我把命令放在下面了，请高手帮忙看一下，为什么一个简单的不考虑重力的 均匀应力场，居然平衡不了！愁死了，其中有两种材料，上面为泥岩，下面是岩盐，中间加了个界面！考虑受25MP的均匀应力场！万分感谢!
gen zone radcylinder p0 0 0 0 p1 60 0 0 p2 0 30 0 p3 0 0 40 dimension 30 30 30 30 size 10 30 60 15 fill group null
gen zone brick p0 60 0 0 p1 120 0 0 p2 60 30 0 p3 60 0 40 size 30 15 30
gen zone reflect normal 0 0 -1 origin 0 0 0
gen zone brick p0 0 0 -100 p1 120 0 -100 p2 0 30 -100 p3 0 0 -40 size 60 15 30
gen zone brick p0 0 0 50 p1 120 0 50 p2 0 30 50 p3 0 0 120 size 60 15 35
gen zone reflect normal -1 0 0 origin 0 0 0
pl bl group
group rocksalt range z -200 45
group mudstone range z 45 200
group null range cylinder end1 0 0 0 end2 0 30 0 radius 30
int 1 face range plane origin (0,0,40) norm (0,0,1) dist 0.001
int 1 prop coh 0 fric 30 kn 19.258e11 ks 19.258e5 tension 0 bslip=on
plo int yel ske
pl bl group
ini z add -10 range group mudstone
model mohr range group rocksalt
model mohr range group mudstone
model mohr range group null
def derive ; caculate bulk and shear by E and u
s_mod=y_mod/(2.0*(1.0+p_ratio))
b_mod=y_mod/(3.0*(1.0-2.0*p_ratio))
end
set y_mod=21.8e9 p_ratio=0.23
derive
prop bulk b_mod shear s_mod coh 0 tens 1e6 fri 35 range group mudstone
def derive2
a_mod=l_mod/(2.0*(1.0+q_ratio))
c_mod=l_mod/(3.0*(1.0-2.0*q_ratio))
end
set l_mod=5.15e9 q_ratio=0.3
derive2
prop bulk c_mod shear a_mod coh 0 tens 1e6 fri 25 range group rocksalt
prop bulk c_mod shear a_mod coh 0 tens 1e6 fri 25 range group null
set grav 0 0 0
fix y range y -0.1 0.1
fix y range y 29.9 30.1
apply nstress -25e6 range z -99.9 -100.1
apply nstress -25e6 range z 109.9 110.1
apply nstress -25e6 range x 119.9 120.1
apply nstress -25e6 range x -119.9 -120.1
set mech force 0.001
so
