找回密码
 注册
Simdroid-非首页
查看: 296|回复: 3

[计算分析] 请教修正剑桥模型的问题

[复制链接]
发表于 2016-1-5 20:10:23 | 显示全部楼层 |阅读模式 来自 重庆沙坪坝区
建了一个基坑的模型,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


本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

×
发表于 2016-2-25 14:48:16 | 显示全部楼层 来自 福建福州
Simdroid开发平台
问问题最好把模型和命令流给出来,或者给个类似的简单的例子,光靠这点描述,估计很难有人可以帮到你。
回复 不支持

使用道具 举报

发表于 2018-4-23 19:43:17 | 显示全部楼层 来自 甘肃兰州
楼主,我也遇到这个问题了,实在求教啊
回复 不支持

使用道具 举报

发表于 2018-7-15 18:35:28 | 显示全部楼层 来自 美国
我喜欢,请大家鼓掌支持,谢谢












四虎影库
辣文
四虎影库
伦理片
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-3-28 17:55 , Processed in 0.031147 second(s), 11 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表