【请教】运用修正剑桥模型时的问题

lin182838

小弟在用运用剑桥模型时，提示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

新手路过，完全看不懂，不好意思

lin182838

tinyo1982 发表于 2014-12-9 19:19
新手路过，完全看不懂，不好意思

没关系，一起学习嘛

lin182838

关于这个问题，已经有了解决办法，就是在广义剪应力后面再加上一些项。上面的命令流中用的是sxx，syy，szz，而打的表达式应该是对应s1，s2，s3的，所以有问题。

但现在又有新的问题，就是计算时候会出现“Yield envelope cannot be reached”，目前在各大论坛都没有找到好的解决办法。希望高手不吝赐教，谢谢！

泰勒菲

哥们，问题解决了吗？你的初始应力场计算的是不是有点问题哈!

泰勒菲

“Yield envelope cannot be reached”主要是什么原因造成的！

脱缰的土狗

还是“Yield envelope cannot be reached”，有没有大牛指导一下啊

Icedragonboom

学习了。。。。。。、。。。。。。。。