【请教】运用修正剑桥模型时的问题
小弟在用运用剑桥模型时，提示trying to divide by zero，但是不知道那里出现了0做除数……n
ge zo cshell p0 0 0 15 p1 4 0 15 p2 0 0 9 p3 0 4 15 p4 4 0 9 p5 0 4 9 p8 0.1155 0 15 &
p9 0 0.1155 15 p10 0.1155 0 9 p11 0 0.1155 9 size 5 6 5 2 ratio 1.2 1 1 1
ge zo cshell p0 0 0 9 p1 4 0 9 p2 0 0 0 p3 0 4 9 p4 4 0 0 p5 0 4 0 p8 0.1155 0 9 &
p9 0 0.1155 9 p10 0.1155 0 0 p11 0 0.1155 0 size 5 6 5 2 ratio 1.2 1 1 1 fill
plo sur lbl
;pause
;initial state
mo el
prop y=40e6 p=0.35 den=1.8e3 ; c=10e3 f=10 t=1e3
set grav 0 0 -10
ini szz 270e3 grad 0 0 -18e3
ini sxx 135e3 grad 0 0 -9e3
ini syy 135e3 grad 0 0 -9e3
mo cam-clay
def install_prop
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
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*pi/180)/(3.0-sin(fai*pi/180))
temp1=q0/(z_prop(pnt,'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')*_mpl))+(_lambda-_kappa)*ln(2.0)
z_prop(pnt,'bulk_bound')=100*40e6
;z_prop(pnt,'bulk_bound')=100*(s_mod+4.0/3.0*s_mod)
pnt=z_next(pnt)
endloop
end
set p_ratio=0.35 fai=34.5 _lambda=0.14 _kappa=0.01 _mp1=1e3 _e0=1.2 _cu=13.9e3 OCR=1.0
install_prop
以上是剑桥模型的设置，希望大家多多指教，谢谢！ 新手路过，完全看不懂，不好意思 tinyo1982 发表于 2014-12-9 19:19
新手路过，完全看不懂，不好意思
没关系，一起学习嘛 关于这个问题，已经有了解决办法，就是在广义剪应力后面再加上一些项。上面的命令流中用的是sxx，syy，szz，而打的表达式应该是对应s1，s2，s3的，所以有问题。
但现在又有新的问题，就是计算时候会出现“Yield envelope cannot be reached”，目前在各大论坛都没有找到好的解决办法。希望高手不吝赐教，谢谢！ 哥们，问题解决了吗？你的初始应力场计算的是不是有点问题哈! “Yield envelope cannot be reached”主要是什么原因造成的！ 还是“Yield envelope cannot be reached”，有没有大牛指导一下啊 学习了。。。。。。、。。。。。。。。
页:
[1]