level 1
大喜剧师
楼主
根据书上的方法，在使用FISH更改修正剑桥模型参数时
def cam_prop
pnt=zone_head
loop while pnt # null
if z_group(pnt)=string(grpnam) then
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+3.0*((z_sxy(pnt))^2+(z_sxz(pnt))^2+(z_syz(pnt))^2))
z_prop(pnt,'mm')=6.0*sin(fai*degrad)/(3.0-sin(fai*degrad)) ;?é¨′?á?騨y??o?￥¨oy?DD¨¨?????¨¨?áa??￥?a???¨¨
temp1=q0/(z_prop(pnt,'mm')*p0_effective)
pc0=p0_effective*(1.0+temp1^2)*OCR ;?¨¨?¨2¨?DD?ì1¨??¨￠?1¨￠|?ê?¨??¨?¨2¨¨?¤??§?¨1?¤t??
v0=1.0+_e0
z_prop(pnt,'cam_cp')=p0_effective ;????¨°a2?¨oy?ê??¤??¨°2??¨1?y¨¨?¤????¨?DD?ì¨?|¨￠|?ê?¨?¨￠¨o?3???¨a"Mean effective pressure is negative"
z_prop(pnt,'mpc')=pc0 ;?é¨2
z_prop(pnt,'poisson')=p_ratio ;?é?
z_prop(pnt,'lambda')=_lambda ;?é¨1
z_prop(pnt,'kappa')=_kappa ;?éY
z_prop(pnt,'mp1')=_mp1 ;?éT
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*40e3 ;?騤
;z_prop(pnt,'bulk_bound')=100*(s_mod+4.0/3.0*s_mod) ;|ì?￥D?¨????ê¨|???Kmax
endif
pnt=z_next(pnt)
endloop
end
提示temp1=q0/(z_prop(pnt,'mm')*p0_effective)行trying to devided by zero
是值q0是零还是p0_effective是0
2017年03月21日 13点03分
1
def cam_prop
pnt=zone_head
loop while pnt # null
if z_group(pnt)=string(grpnam) then
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+3.0*((z_sxy(pnt))^2+(z_sxz(pnt))^2+(z_syz(pnt))^2))
z_prop(pnt,'mm')=6.0*sin(fai*degrad)/(3.0-sin(fai*degrad)) ;?é¨′?á?騨y??o?￥¨oy?DD¨¨?????¨¨?áa??￥?a???¨¨
temp1=q0/(z_prop(pnt,'mm')*p0_effective)
pc0=p0_effective*(1.0+temp1^2)*OCR ;?¨¨?¨2¨?DD?ì1¨??¨￠?1¨￠|?ê?¨??¨?¨2¨¨?¤??§?¨1?¤t??
v0=1.0+_e0
z_prop(pnt,'cam_cp')=p0_effective ;????¨°a2?¨oy?ê??¤??¨°2??¨1?y¨¨?¤????¨?DD?ì¨?|¨￠|?ê?¨?¨￠¨o?3???¨a"Mean effective pressure is negative"
z_prop(pnt,'mpc')=pc0 ;?é¨2
z_prop(pnt,'poisson')=p_ratio ;?é?
z_prop(pnt,'lambda')=_lambda ;?é¨1
z_prop(pnt,'kappa')=_kappa ;?éY
z_prop(pnt,'mp1')=_mp1 ;?éT
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*40e3 ;?騤
;z_prop(pnt,'bulk_bound')=100*(s_mod+4.0/3.0*s_mod) ;|ì?￥D?¨????ê¨|???Kmax
endif
pnt=z_next(pnt)
endloop
end
提示temp1=q0/(z_prop(pnt,'mm')*p0_effective)行trying to devided by zero
是值q0是零还是p0_effective是0
