[求助]FLAC边坡安全系数的计算 [复制链接]

使用FLAC3d计算某一边坡的安全系数，但是再计算出初始应力状态以后，在计算安全系数时出现计算错误 +2=N#LM  

图片:1111111.jpg

`#hdb=3  
程序 [^B04
x@  
new KwOn<0P  
gen zone brick p0 0 0 0 p1 140 0 0 p2 0 20 0 p3 104 0 120 p4 140 20 0 p5 104 20 120 p6 140 0 120 p7 140 20 120 size 10 5 12 'u$e2^  
gen zone brick p0 -14 0 -20 p1 140 0 -20 p2 -14 20 -20 p3 -14 0 0 size 11 5 2 5An|#^]  
fix x y z range z -20 RQ#gn  
fix x range x 140 "q4c[dna  
fix x range x -14 ++-\^'&1  
fix y GcQO&oq|  
;初始条件设置 8 ~.|^no  
model e ;JMd(\+-  
prop density 2800 bulk 4.72e12 shear 3.54e12 KFBo1^9N  
set gravity 0 0 -10 @a) x^d  
solve %zQME6WELz  
save solpe1-chu-1.sav r'M|mQ$s>  
;位移、速度清零 >(Jy=m?  
ini xdisp 0 ydisp 0 zdisp 0 QA\eXnR  
ini xvel 0 yvel 0 zvel 0 /u{ 9UR[g  
;安全系数求解 ,JyE7h2%i  
model mohr u1O?
`  
prop density 2800 bulk 4.72e10 shear 3.54e10 coh 4.5e7 tens 1e6 fric 35 dila 35 dhpEBJ
  
solve fos file slopepoumian1fos2201.sav associated 3!&lio+<  
K:fK!/  
zr%lBHuW  
请各位大神帮忙提一下修改意见，谢谢 .asHFT7]9  
]Qe;+p9vU  

参数设置有问题，不收敛