[求助]管状建模初始应力生成问题 [复制链接]

本人模拟的是管桩静载试验，现在有一点迷惑的地方就是：建模的时候桩是弹性模型，桩周土是弹塑性模型（m-c），我想在施加第一级荷载前先生成初始应力平衡，但是对陈老师书上介绍的三种应力生成方法都拿不准，不知道该选择哪种来生成初始应力，求高人指点！！！！谢谢 a{xJ#_/6  
i+14!LlI  
下面附带我的命令流： U5C]zswL  
;桩周土 ,\i*vJ#f  
gen zone radcyl p0 0 0 0 p1 8 0 0 p2 0 0 -5 p3 0 8 0 p4 8 0 -5  p5 0 8 -5  p6 8 8 0 p7 8 8 -5  p8 0.25 0 0 p9 0 0.25 0 p10 0.25 0 -5  p11 0 0.25 -5  size 3 10 6 15 ratio 1 1 1 !L.R"8!  
1.15 )B]s.
w  
gen zone radcyl p0 0 0 -5  p1 8 0 -5  p2 0 0 -8  p3 0 8 -5  p4 8 0 -8   p5 0 8 -8   p6 8 8 -5  p7 8 8 -8  p8 0.25 0 -5  p9 0 0.25 -5  p10 0.25 0 -8   p11 0 0.25 -8  size 3 6 6 [;wJM|Z J0  
15   ratio 1 1 1 1.15  fill kTH""h{  
group clay =@d#@  
；桩与桩周土接触面 CcUF)$kz  
interface 1 face range cylinder end1 0 0 0 end2 0 0 -5.1 radius 0.26  cylinder end1 0 0 0 end2 0 0 -5.1 radius 0.24  not w
1I07 (  
；桩端土接触面 FO/cEu  
interface 2 face range cylinder end1 0 0 -4.9 end2 0 0 -5.1 radius 0.26 z%E(o%l8  
; 桩的承台厚10cm Tw';;euw  
gen zone cylinder p0=0 0 6.1  p1=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 ZbC$Fk,,I&  
  size 3 1 6 lG-B)
F  
；桩-柱形壳体网格建立桩的模型 <}lah%4F  
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 ?1L<VL=b  
_GkLspSaU  
group pile ran z 1 6.1 f+9eB
  
；桩下移合并 wn@~80)$  
ini z add -6.0 range group pile 8=$XhC  
QKjn/%l"@  
；桩周土材料模型跟参数 68j1svz9  
model mohr   range group clay ,< g%
}P/  
prop bulk 7.78e6   shear 4.33e6   coh 4400  fric  25  range group clay HN7tIz@Frc  
；桩的模型跟参数 /k/X[/WO  
model elas     range  group   pile m}z6Bbis0  
prop bulk 2.11e10  shear 1.58e10   range group pile -
F?97&G$  
；接触面参数 q;[HUyY,  
interface 1 prop kn  1e8  ks  1e8  fric 12  coh  2200 a]*{!V{$i  
interface 2 prop kn  1e8  ks  1e8  fric 12  coh  2200 x_~_/&X5  
；桩周土跟桩的密度 WOn<JCh]  
ini  dens 14200 range group clay curYD~7  
ini  dens 28000 range group pile oaQW~R`_  
(eF[nfM  
；边界条件-模型底部全约束 Q
crhgR  
fix  z range z -8.1 -7.9 'ge$}L}4  
fix  x range x -0.1 0.1 9C)VW  
fix  x range x 7.9 8.1 f_)#  
fix  y range y -0.1 0.1  el2Wk@*  
fix  y range y 7.9 8.1 &?y@`',a0{  
；监测桩顶的竖向位移 Ub\^3f  
history id=1 gp zdis 0 0 0.1 w<H2#d>5!@  
;设置初始应力 sg=G<50i  
set  grav 0 0 -10 B9|s`o)!  
so dDA8IW![S  
save pile0.sav 4L,wBce;,t  
;输出桩顶的位移 - BWf.  
plot {d<XDx4`  
history 1 0UJ6>Rj  
;施加第一级荷载 yf&_l^!  
ini xdis 0.0 ydis 0.0 zdis 0.0 >>$L vQ  
apply szz -8.15e6   ran z 0 0.1 group pile   cO]w*Hti  
solve 8KJ`+"<=@  
save app1.sav