求助!!计算坡体安全系数!!急!!

simsky

安全度!!求救我在计算一个边坡安全系数时总是出错,老是出现这样的提示" 'cohesion' not found in zone 1"。求高手指导下！
计算代码如下：
网格及代码见附件
model mohr  
pro bulk 3.3e7 shear 1.5e7 fric 25 cohesion 15000 ten 0 dens 1800  
fix x y z range y -10.87 -10.85
fix x range x -10.83 -10.81
fix x range x 29.17 29.19
fix z range z -0.01 0.01
fix z range z 0.99 1.01
def zonesafety
p_z=zone_head
loop while p_z # null
sigma1=z_sig1(p_z)
sigma3=z_sig3(p_z)
c=z_prop(p_z,' cohesion')
phi=z_prop(p_z,' friction')
sigmat=z_prop(p_z,' tension')
sigmac=2*c*sqrt((1+sin(phi*pi/180))/(1-sin(phi*pi/180)))
if sigma3>0
if sigma1>0
f=sigmat/sigma3
else
if abs(sigma1)<=abs(sigma3)
f=sigmat/sigma3
endif
if abs(sigma1)>abs(sigma3)
f=(sigma*sigmac)/(sigmac*sigma3-sigmat*(sigma1+sigma3))
endif
endif
endif
if sigma3<=0
f=-1*sigmac/sigma1
endif
z_extra(p_z,1)=f
p_=z_next(p_z)
endloop
end
zonesafety
hist unbal
set grav 0 -10 0
solve
求斑竹们指点下！小弟万分感谢!!

Stphenche

既然是摩尔库伦模型用solve fos 就行了

成吉思汗5578

直接用 flac slope 出来的就是安全系数

li_2001

直接用 flac slope 出来的就是安全系数