管状建模初始应力生成问题
本人模拟的是管桩静载试验，现在有一点迷惑的地方就是：建模的时候桩是弹性模型，桩周土是弹塑性模型（m-c），我想在施加第一级荷载前先生成初始应力平衡，但是对陈老师书上介绍的三种应力生成方法都拿不准，不知道该选择哪种来生成初始应力，求高人指点！！！！谢谢下面附带我的命令流：
;桩周土
gen zone radcyl p0 0 0 0 p1 8 0 0 p2 0 0 -5 p3 0 8 0 p4 8 0 -5p5 0 8 -5p6 8 8 0 p7 8 8 -5p8 0.25 0 0 p9 0 0.25 0 p10 0.25 0 -5p11 0 0.25 -5size 3 10 6 15 ratio 1 1 1
1.15
gen zone radcyl p0 0 0 -5p1 8 0 -5p2 0 0 -8p3 0 8 -5p4 8 0 -8 p5 0 8 -8 p6 8 8 -5p7 8 8 -8p8 0.25 0 -5p9 0 0.25 -5p10 0.25 0 -8 p11 0 0.25 -8size 3 6 6
15 ratio 1 1 1 1.15fill
group clay
；桩与桩周土接触面
interface 1 face range cylinder end1 0 0 0 end2 0 0 -5.1 radius 0.26cylinder end1 0 0 0 end2 0 0 -5.1 radius 0.24not
；桩端土接触面
interface 2 face range cylinder end1 0 0 -4.9 end2 0 0 -5.1 radius 0.26
; 桩的承台厚10cm
gen zone cylinder p0=0 0 6.1p1=0.25 0 6.1 p2=0 0 6 p3=0 0.25 6.1 p4=0.25 0 6 p5=0 0.25 6
size 3 1 6
；桩-柱形壳体网格建立桩的模型
gen zone cshell p0=0 0 6 p1=0.25 0 6 p2=0 0 1 p3=0 0.25 6 p4=0.25 0 1 p5=0 0.25 1 p8=0.15 0 6 p9=0 0.15 6 p10=0.15 0 1 p11=0 0.15 1 size 1 10 6 3 ratio 1 1 1.15 1
group pile ran z 1 6.1
；桩下移合并
ini z add -6.0 range group pile
；桩周土材料模型跟参数
model mohr range group clay
prop bulk 7.78e6 shear 4.33e6 coh 4400fric25range group clay
；桩的模型跟参数
model elas rangegroup pile
prop bulk 2.11e10shear 1.58e10 range group pile
；接触面参数
interface 1 prop kn1e8ks1e8fric 12coh2200
interface 2 prop kn1e8ks1e8fric 12coh2200
；桩周土跟桩的密度
inidens 14200 range group clay
inidens 28000 range group pile
；边界条件-模型底部全约束
fixz range z -8.1 -7.9
fixx range x -0.1 0.1
fixx range x 7.9 8.1
fixy range y -0.1 0.1
fixy range y 7.9 8.1
；监测桩顶的竖向位移
history id=1 gp zdis 0 0 0.1
;设置初始应力
setgrav 0 0 -10
so
save pile0.sav
;输出桩顶的位移
plot
history 1
;施加第一级荷载
ini xdis 0.0 ydis 0.0 zdis 0.0
apply szz -8.15e6 ran z 0 0.1 group pile
solve
save app1.sav 我觉得初始应力前都设成弹性的,初始完再改回实际的模型,然后ini 位移和速度,进行计算
页:
[1]