- 积分
- 0
- 注册时间
- 2016-1-5
- 仿真币
-
- 最后登录
- 1970-1-1
|
建了一个基坑的模型,1到7是坑内土,8-16是坑外土,其中2和9是软土。我先用摩尔库伦本构模型计算了地应力,并将速度位移都清零了。
现在想将2和9改为修正剑桥模型计算基坑内外变形,其他底层继续采用摩尔库伦模型。
在论坛上借鉴了丁勇春博士的成果。但是总提示我property bad
在地应力清零并加上结果单元后的命令流:
def cam_prop
pnt=zone_head
loop while pnt # null
if z_group(pnt)=string(grpnam) then
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)
z_prop(pnt,'mm')=6.0*sin(fai*degrad)/(3.0-sin(fai*degrad)) ;①注三角函数中需将角度转化为弧度
temp1=q0/(z_prop(pnt,'mm')*p0_effective+1.0)
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*10e6 ;⑧
;z_prop(pnt,'bulk_bound')=100*(b_mod+4.0/3.0*s_mod) ;弹性体模上界Kmax
;自动确定Kmax时会出现“property bad”错误提示
;因为弹性上界对计算结果无影响,在不提示Kmax太小的性况下,取值越小计算收敛越快
endif
pnt=z_next(pnt)
endloop
end
model cam-clay
set grpnam='2' p_ratio=0.37 fai=8.972 _lambda=0.0620 _kappa=0.0271 _mp1=1e3 _e0=1.25 _cu=10e3 OCR=1.0
cam_prop
set grpnam='9' p_ratio=0.37 fai=8.972 _lambda=0.1305 _kappa=0.0281 _mp1=1e3 _e0=1.25 _cu=10e3 OCR=1.0
cam_prop
|
本帖子中包含更多资源
您需要 登录 才可以下载或查看,没有账号?注册
×
|