本人模拟了一个实例工程，但是赋值方面出现了问题，老是提示某节点刚度为零。
我找了好久，始终没有找出究竟哪里的值出现错误了，现在希望FLAC3D学的好的大神们帮我解决一下这个问题吧！！！！！！
跪求！！！！！拜托了。希望大神能将的仔细些！！！！！
命令流如下：
new
;=============土层单元===============
gen zone brick size 100 2 60
group soil1 range z 57.5 60
group soil2 range z 54.5 57.5
group soil3 range z 51.5 54.5
group soil4 range z 43.5 51.5
group soil5 range z 37.5 43.5
group soil6 range z 31.5 37.5
group soil7 range z 21.5 31.5
group soil8 range z 0 21.5
group exc1 range x 0 30 y 0 2 z 60 54.9
group exc2 range x 0 30 y 0 2 z 54.9 52.2
group exc3 range x 0 30 y 0 2 z 52.2 49.5
group exc4 range x 0 30 y 0 2 z 49.5 46.8
group exc5 range x 0 30 y 0 2 z 46.8 44.1
group exc6 range x 0 30 y 0 2 z 44.1 41.4
group exc7 range x 0 30 y 0 2 z 41.4 40
;===============创建桩体=================
gen zone brick p0 30.8 0 87 p1 30.8 0.8 87 p2 30.8 0 61 p3 30 0 87 size 8 8 26
gen zone brick p0 30.8 0.8 87 p1 30.8 1.2 87 p2 30.8 0.8 61 p3 30 0.8 87 size 8 4 26
gen zone brick p0 30.8 1.2 87 p1 30.8 2 87 p2 30.8 1.2 61 p3 30 1.2 87 size 8 8 26
group pile ran z 61 87
ini z add -31.5 ran group pile
;=================连接面的建立===============
interface 1 face ran x 30 30.8 y 0.7 0.9 z 29.5 55.5
interface 2 face ran x 30 30.8 y 1.1 1.3 z 29.5 55.5
interface 3 face ran x 29.9 30.1 y 0 2 z 29.5 55.5
interface 4 face ran x 30.7 30.9 y 0 2 z 29.5 55.5
interface 5 face ran x 30 30.8 y 0 2 z 29.4 29.6
interface 1 prop kn 20e6 ks 20e6 fric 15 cohesion 10e3 tension 1e10
interface 2 prop kn 20e6 ks 20e6 fric 15 cohesion 10e3 tension 1e10
interface 3 prop kn 10e6 ks 10e6 fric 10.8 cohesion 8.5e3 tension 1e10
interface 4 prop kn 10e6 ks 10e6 fric 10.8 cohesion 8.5e3 tension 1e10
interface 5 prop kn 10e6 ks 10e6 fric 10.8 cohesion 8.5e3 tension 1e10
;==================土体和桩体参数=================
model mohr
prop bulk 6.90e6 shear 2.4e6 fric 13.5 cohesion 8.50e3 dilation 0 tension 1e10 density 1880 range group soil8
ini dens 2300 range group soil8
prop bulk 5.30e6 shear 1.8e6 fric 10.3 cohesion 10.7e3 dilation 0 tension 1e10 density 1910 range group soil7
ini dens 2300 range group soil7
prop bulk 6.10e6 shear 2.0e6 fric 13.3 cohesion 16.3e3 dilation 0 tension 1e10 density 1980 range group soil6
ini dens 2300 range group soil6
prop bulk 18.57e6 shear 24.70e6 fric 14.40 cohesion 31.80e3 dilation 0 tension 1e10 density 2080 range group soil5
ini dens 2300 range group soil5
prop bulk 9.70e6 shear 12.91e6 fric 13.50 cohesion 16.40e3 dilation 0 tension 1e10 density 1980 range group soil4
ini dens 2300 range group soil4
prop bulk 8.56e6 shear 11.40e6 fric 10.40 cohesion 10.80e3 dilation 0 tension 1e10 density 1920 range group soil3
ini dens 2300 range group soil3
prop bulk 11.56e6 shear 15.40e6 fric 13.50 cohesion 8.50e3 dilation 0 tension 1e10 density 1880 range group soil2
ini dens 2300 range group soil2
prop bulk 5.92e6 shear 7.88e6 fric 16.50 cohesion 19.00e3 dilation 0 tension 1e10 density 1900 range group soil1
ini dens 2300 range group soil1
def ini_szz
_grad=-1880*(-9.81)
_szz1=-_grad*(60-57.5)
_szz0=-_grad*60
command
ini szz add _szz0 grad 0 0 _grad range z 57.5 60
ini szz add _szz1 range z 0 60
end_command
_grad=-1910*(-9.81)
_szz1=-_grad*(57.5-54.5)
_szz0=-_grad*36.5
command
ini szz add _szz0 grad 0 0 _grad range z 54.5 57.5
ini szz add _szz1 range z 0 54.5
end_command
_grad=-1980*(-9.81)
_szz1=-_grad*(54.5-51.5)
_szz0=-_grad*51.5
command
ini szz add _szz0 grad 0 0 _grad range z 51.5 54.5
ini szz add _szz1 range z 0 51.5
end_command
_grad=-2080*(-9.81)
_szz1=-_grad*(51.5-43.5)
_szz0=-_grad*43.5
command
ini szz add _szz0 grad 0 0 _grad range z 51.5 43.5
ini szz add _szz1 range z 0 43.5
end_command
_grad=-1980*(-9.81)
_szz1=-_grad*(43.5-37.5)
_szz0=-_grad*37.5
command
ini szz add _szz0 grad 0 0 _grad range z 37.5 43.5
ini szz add _szz1 range z 0 37.5
end_command
_grad=-1920*(-9.81)
_szz1=-_grad*(37.5-31.5)
_szz0=-_grad*31.5
command
ini szz add _szz0 grad 0 0 _grad range z 31.5 37.5
ini szz add _szz1 range z 0 31.5
end_command
_grad=-1880*(-9.81)
_szz1=-_grad*(31.5-21.5)
_szz0=-_grad*21.5
command
ini szz add _szz0 grad 0 0 _grad range z 21.5 31.5
ini szz add _szz1 range z 0 21.5
end_command
_grad=-1900*(-9.81)
_szz0=-_grad*21.5
command
ini szz add _szz0 grad 0 0 _grad range z 0 21.5
end_command
end
def ini_conf
pnt=zone_head
loop while pnt # null
val=0.43*z_szz(pnt)
z_sxx(pnt)=val
z_syy(pnt)=val
pnt=z_next(pnt)
end_loop
end
ini_conf
;
model mohr range group pile
prop bulk 2.0e10 shear 1.0e10 fric 30 cohesion 3.0e6 ten 2e6 range group pile id=1
prop bulk 2.5e10 shear 1.15e10 fric 50 cohesion 3.8e6 ten 4e6 range group pile id=2
prop bulk 2.0e10 shear 1.0e10 fric 30 cohesion 3.0e6 ten 2e6 range group pile id=3
ini dens 2500 range group pile
;=============创建梁单元===============
sel beam id=1 begin 30 0 55.8 end 30 2 55.8 nseg=2
sel beam prop emod=30e9 nu=0.2 density=2.5e3 ran sel beam
;=============边界条件=================
fix z range z -0.1 0.1
fix x range x -0.1 0.1
fix x range x 99.9 100.1
fix y range y -0.1 0.1
fix y range y 1.9 2.1
;===============应力条件================
apply szz -20e3 ran z 59.9 60.1
ini szz 0e3
ini sxx 0e3
ini syy 0e3
set grav 0 0 -9.81
hist unbal
set large
solve
save 自重平衡.sav
