应用修正剑桥模型进行flac计算时，出现Degenerate yield envelope detected for Cam-clay material negative square root encountered.的错误 ，我感觉是参数赋的不对，但不知道怎么修改，望大神们指导呀。我之前的参数如下：den 2000 bulk 3.3e7 shear 1.2e6 mm 0.96 lambda 0.09377 kappa 0.02344 mp1 1 mv_l 3.27
mpc赋值函数：
def set_prop
loop i (1,izones)
loop j (1,jzones) s1=-sxx(i,j)
s2=-syy(i,j) s3=-szz(i,j)
p0=(s1+s2+s3)/3-pp(i,j) q0=sqrt(((s1-s2)^2+(s2-s3)^2+(s3-s1)^2)*0.5)
e0=q0/(mm(i,j)*p0) pc=p0*(1+e0*e0)*1.2
command prop mpc pc ;bulk bm
endcommand end_loop
end_loopend
set_prop
