flac3d多层土设置修正剑桥模型

;修正剑桥模型

;cam-clay模型则不需定义弹性模量(E、G、K)等参数，自动计算

;cam-clay模型中需确定8个模型参数(①-⑧)和一个分组名称，手册property中的初始比体积cv(v0)和shear无须给定

def cam_prop

pnt=zone_head

loop while pnt # null

abs_sxx=abs(z_sxx(pnt)) ;|sxx|

abs_syy=abs(z_syy(pnt)) ;|syy|

abs_szz=abs(z_szz(pnt)) ;|szz|

p0=(abs_sxx+abs_syy+abs_szz)/3.0

;cam-clay模型中p、q均须为正值，p0由初应力场确定，故cam-clam定义模型参数前须先已知初应力

p0_effective=p0-z_pp(pnt) ;p0'

;q0=sqrt(((abs_sxx-abs_syy)^2+(abs_syy-abs_szz)^2+(abs_szz-abs_sxx)^2)*0.5)

        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))  ;①注三角函数中需将角度转化为弧度

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 ;★重要参数，否则不能正确计算有效应力，提示出错"Mean effective pressure is negative"

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 ;⑧

;z_prop(pnt,'bulk_bound')=100*(s_mod+4.0/3.0*s_mod) ;弹性体模上界Kmax

;自动确定Kmax时会出现“property bad”错误提示

;因为弹性上界对计算结果无影响，在不提示Kmax太小的性况下，取值越小计算收敛越快

endif

pnt=z_next(pnt)

endloop

end

model cam-clay

set grpnam='soil_3_x_out' p_ratio=0.24 fai=30 _lambda=0.14 _kappa=0.012 _mp1=1e3 _e0=0.8 _cu=10e3 OCR=1.0  ;模型所需参数

cam_prop

set grpnam='soil_1_x_out' p_ratio=0.24 fai=30 _lambda=0.14 _kappa=0.90  _mp1=1e3 _e0=0.8 _cu=10e3 OCR=1.0  ;模型所需参数

cam_prop

from: http://blog.sina.com.cn/s/blog_6c64ec090101ipvh.html