Flac3d 计算不收敛，求大神指导，实在不知道问题在哪
new
generate zone brick p0 0 0 0 p1 20 0 0 p2 0 0.1 0 p3 0 0 20 size 100 1 100
model elas
pro bulk 1.71e10 shear 8.80e9
apply sxx -5.23e7 range x 0
apply sxx -5.23e7 range x 20
apply syy -5.47e7 range y 0
apply syy -5.47e7 range y 0.1
apply szz -6.59e7 range z 0
apply szz -6.59e7 range z 20
apply sxy -2.11e6 range x 0
apply sxz 3.43e6 range x 0
apply sxy -2.11e6 range x 20
apply sxz 3.43e6 range x 20
apply syx -2.11e6 range y 0
apply syz -4.23e6 range y 0
apply syx -2.11e6 range y 0.1
apply syz -4.23e6 range y 0.1
apply szx 3.43e6 range z 0
apply szy -4.23e6 range z 0
apply szx 3.43e6 range z 20
apply szy -4.23e6 range z 20
ini sxx -5.23e7
ini syy -5.47e7
ini szz -6.59e7
ini sxy -2.11e6
ini syz -4.23e6
ini sxz 3.43e6
solve
model mohr
pro bulk 1.71e10 shear 8.80e9 friction 32.27 coh 7.8e5 ten 0.76e5
solve
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
hist id 1 unbal
model null range x 9 11 y 0 0.1 z 9 11
step 4000
