请问为什么修正剑桥模型中提示参数错误？麻烦大神帮看看代码！ [复制链接]

RT请问修正剑桥模型中为什么提示参数错误？ 总是这样*** Bad conversion of parameter number 2 (p_ratio). 图片:QQ图片20160801173410.png 命令流如下请大神帮看看！ new gen z b p0 0 0 188.35 p1 2000 0 188.35 p2 0 2000 188.35 p3 0 0 200.5 size 100 100 2 ratio 1 1 1 group diceng1 gen z b p0 0 0 175 p1 2000 0 175 p2 0 2000 175 p3 0 0 188.35 size 100 100 2 ratio 1 1 1 group diceng2 gen z b p0 0 0 165.1 p1 2000 0 165.1 p2 0 2000 165.1 p3 0 0 175 size 50 50 1 ratio 1 1 1 group diceng3 gen z b p0 0 0 154.8 p1 2000 0 154.8 p2 0 2000 154.8 p3 0 0 165.1 size 50 50 1 ratio 1 1 1 group diceng4 gen z b p0 0 0 146 p1 2000 0 146 p2 0 2000 146 p3 0 0 154.8 size 50 50 1 ratio 1 1 1 group diceng5 gen z b p0 0 0 91 p1 2000 0 91 p2 0 2000 91 p3 0 0 146 size 10 10 3 ratio 1 1 1 group diceng6 gen z b p0 0 0 84.7 p1 2000 0 84.7 p2 0 2000 84.7 p3 0 0 91 size 10 10 1 ratio 1 1 1 group diceng7 gen z b p0 0 0 60.5 p1 2000 0 60.5 p2 0 2000 60.5 p3 0 0 84.7 size 5 5 1 ratio 1 1 1 group diceng8 gen z b p0 0 0 49.3 p1 2000 0 49.3 p2 0 2000 49.3 p3 0 0 60.5 size 5 5 1 ratio 1 1 1 group diceng9 gen z b p0 0 0 44.5 p1 2000 0 44.5 p2 0 2000 44.5 p3 0 0 49.3 size 5 5 1 ratio 1 1 1 group diceng10 gen z b p0 0 0 16.1 p1 2000 0 16.1 p2 0 2000 16.1 p3 0 0 44.5 size 1 1 1 ratio 1 1 1 group diceng11 gen z b p0 0 0 10.5 p1 2000 0 10.5 p2 0 2000 10.5 p3 0 0 16.1 size 1 1 1 ratio 1 1 1 group diceng12 gen z b p0 0 0 0 p1 2000 0 0 p2 0 2000 0 p3 0 0 10.5 size 1 1 1 ratio 1 1 1 group diceng13 attach face range z 0 201 Initial density 2000 Fix x range x=-0.1 0.1 Fix x range x=2000 fix y range y=-0.1 0.1 fix y range y=2000 fix z range z=0 set gravity 0 0 -9.8 model cam-clay def install_prop pnt=zone_head loop while pnt # null   OCR=1.2   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   p0_effective=p0-z_pp(pnt)   q0=sqrt(((abs_sxx-abs_syy)^2+(abs_syy-abs_szz)^2+(abs_szz-abs_sxx)^2)*0.5)   z_prop(pnt,'mm')=6.0*sin(fai*degrad)/(3.0-sin(fai*degrad))   temp1=q0/(z_prop(pnt,'mm')*p0_effective)   ;temp1=q0/(mm*p0_effective)   pc0=p0_effective*(1.0+temp1^2)*OCR   v0=1.0+_e0   z_prop(pnt,'cam_cp')=p0_effective   z_prop(pnt,'mpc')=pc0   z_prop(pnt,'poisson')=p_ratio   z_prop(pnt,'lambda')=_lambda   z_prop(pnt,'kappa')=_kappa   z_prop(pnt,'mp1')=_mp1   z_prop(pnt,'mv_l')=v0+_lambda*ln(2.0*_cu/(z_prop(pnt,'mm')*_mp1))+(_lambda-_kappa)*ln(2.0)   z_prop(pnt,'bulk_bound')=100*40e6  pnt=z_next(pnt) endloop end set p_ratio=0.24 fai=30 _lambda=0.14 _kappa=0.012 _mp1=1e3 _e0=0.8 _cu=10e3 OCR=1.0 @install_prop solve save diceng4.sav

还有个问题是，flac3d 是不是不允许Z为负？我把地面作为0时，向下地层为负时，就不能赋本构模型了，请问这是怎么回事？

楼主，请问你这个问题怎么解决的，麻烦分享一下方法，谢谢