请问修正剑桥模型中提示参数错误？麻烦大神帮看看代码！

科费卡

RT请问修正剑桥模型中提示参数错误？
总是这样*** Bad conversion of parameter number 2 (p_ratio).

麻烦大神帮看看代码！

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

科费卡

另外想请教一下，我在看相关资料时，没有发现cam_cp是修正剑桥模型的必要参数，但是所有的实例中都出现这个参数，且注明是重要参数，否则不能正确计算有效应力。而fish中也没有找到他是预留的变量。请问cam_cp到底代表什么？为什么必须要有，是否可以用自己随意定义的变量代替？谢谢！

zaiyunduan

楼主 这些你解决了吗