修正剑桥模型出现错误“yield envelope cannot be reached”

zhouzhou8246384

最近在flac3d中运用修正剑桥模型出现这个错误，不知道什么原因及解决办法，求赐教。
restore inidis0.sav
model cam-clay
prop shear 150000 bulk_bound 20e7
prop mm 0.73 lambda 0.0938 kappa 0.0234
prop mpc 0.395e6 mp1 1.0 mv_l 3.32
def install_prop
OCR=2
pnt=zone_head
loop while pnt # null
abs_sxx=abs(z_sxx(pnt))
abs_syy=abs(z_syy(pnt))
abs_szz=abs(z_szz(pnt))
p0=(abs_sxx+abs_syy+abs_szz)/3.0-z_pp(pnt)
q0=sqrt(((abs_sxx-abs_syy)^2+(abs_sxx-abs_szz)^2+(abs_syy-abs_szz)^2)*0.5)
temp1=q0/(z_prop(pnt,'mm')*p0)
pc=p0*(1.0+temp1^2)*OCR
z_prop(pnt,'cam_cp')=p0
z_prop(pnt,'mpc')=pc
pnt=z_next(pnt)
endloop
end
install_prop
cyc 0
solve

zhouzhou8246384

看论坛上的 调整ocr，我试了几次不管用，另个人认为调节ocr后计算结果会不对。

zhouzhou8246384

求各位同行不吝赐教，谢谢