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

ScienceSoft
ScienceSoft
ScienceSoft
240
文章
67
评论
6月 9, 201816:41:29
评论
1
;修正剑桥模型
;cam-clay模型则不需定义弹性模量(E、G、K)等参数,自动计算
;cam-clay模型中需确定8个模型参数(①-⑧)和一个分组名称,手册property中的初始比体积cv(v0)和shear无须给定
def cam_prop
pnt=zone_head
loop while pnt # null
if z_group(pnt)=string(grpnam) then    ;note:新增控制分组参数,grpnam=分组名称。
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
flac3d多层土设置修正剑桥模型
from: http://blog.sina.com.cn/s/blog_6c64ec090101ipvh.html
weinxin
我的微信公众号
分享科研软件、科研方法,为你的科研助力。
ScienceSoft
  • 本文由 发表于 6月 9, 201816:41:29
  • 转载请务必保留本文链接:https://www.sciencesoft.cn/use-cam-clay-model-for-multi-layer-in-flac3d/
FLAC3D有哪些值得一看的参考书? 科研软件

FLAC3D有哪些值得一看的参考书?

FLAC3D 是一种专业的沿途分析软件,也是最具流行度和影响力的数值软件之一。我们可以通过阅读书籍来帮助掌握FLAC3D的相关知识。这里列出了几本值得一看的参考书。 1、FLAC3D 5.0说明手册(...
在Notepad++中一键调用FLAC3D 科研软件

在Notepad++中一键调用FLAC3D

由于在Notepad++中通过自定义语法高亮可以很方便的编写代码,因此很多平时的代码编写都是在Notepad++中完成的。但代码编写完后进行保存再在Flac3D中点击菜单File->call,这...
给FLAC/FLAC3D软件初学者的一些意见 科研方法

给FLAC/FLAC3D软件初学者的一些意见

一、认识和了解FLAC/FLAC3D 有关这两个软件的介绍,在这里不再赘述,网上有大量而详细的材料。需要提示大家的是,这两者都是基于连续介质及有限差分原理的2D(或3D)程序,一般不用于断裂、崩塌、侵...
Office终于原生支持SVG格式图片 办公软件

Office终于原生支持SVG格式图片

最近发现,Office终于原生支持SVG格式的图片了,真是太好了。(我用的Office 365,可以随时更新到最新版本) 以前,我们在Excel中作图之后,如果要插入到Word中,想要得到清晰的图片,...
匿名

发表评论

匿名网友 填写信息

:?: :razz: :sad: :evil: :!: :smile: :oops: :grin: :eek: :shock: :???: :cool: :lol: :mad: :twisted: :roll: :wink: :idea: :arrow: :neutral: :cry: :mrgreen: