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

本人模拟的是管桩静载试验，现在有一点迷惑的地方就是：建模的时候桩是弹性模型，桩周土是弹塑性模型（m-c），我想在施加第一级荷载前先生成初始应力平衡，但是对陈老师书上介绍的三种应力生成方法都拿不准，不知道该选择哪种来生成初始应力，求高人指点！！！！谢谢 ~mi4V  
wArfnB&  
下面附带我的命令流： y%kZ##  
;桩周土 qE]e+S?57a  
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 ReA-.j_2@  
1.15 ijzwct#.  
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 HSAr6h  
15   ratio 1 1 1 1.15  fill LGX+_"  
group clay %:;g|PC  
；桩与桩周土接触面 !7MRHI/0C  
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 P*VZ$bUe5@  
；桩端土接触面 ~(GNY5  
interface 2 face range cylinder end1 0 0 -4.9 end2 0 0 -5.1 radius 0.26 ~<_WYSzS  
; 桩的承台厚10cm *vvm8ik  
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 ae]6F_Qtc*  
  size 3 1 6 YgS,5::SU  
；桩-柱形壳体网格建立桩的模型 Z|ZB6gP>h1  
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 >'.: Acn  
e+{lf*"3  
group pile ran z 1 6.1 rzLW@k  
；桩下移合并 ;ny9q  
ini z add -6.0 range group pile ?.H]Y&XF  
kY{$[+-jR  
；桩周土材料模型跟参数 d2-oy5cEB  
model mohr   range group clay L
NHi}P~  
prop bulk 7.78e6   shear 4.33e6   coh 4400  fric  25  range group clay lmL$0{Yr  
；桩的模型跟参数 ^M;#x$Y?  
model elas     range  group   pile K' <[kh:cl  
prop bulk 2.11e10  shear 1.58e10   range group pile qc\D=3#Yp  
；接触面参数 v=cQ`nou  
interface 1 prop kn  1e8  ks  1e8  fric 12  coh  2200 O7uCTB+  
interface 2 prop kn  1e8  ks  1e8  fric 12  coh  2200 3T4HX|rC  
；桩周土跟桩的密度 jiLJiYMg  
ini  dens 14200 range group clay WY!4^<|w"  
ini  dens 28000 range group pile ?#z<<FR  
dh&>

E  
；边界条件-模型底部全约束 <1<xSr  
fix  z range z -8.1 -7.9 K4YD}[  
fix  x range x -0.1 0.1 A=p'`]Yld  
fix  x range x 7.9 8.1 7\R"RH-  
fix  y range y -0.1 0.1 \4C[<Gbx$(  
fix  y range y 7.9 8.1 .q[}e);)  
；监测桩顶的竖向位移 W-/V5=?   
history id=1 gp zdis 0 0 0.1 }N$f=:i
I  
;设置初始应力 {>~9?Xwh   
set  grav 0 0 -10 EUQtl_h/H  
so ]yZ%wU9!  
save pile0.sav [K'gvLt1  
;输出桩顶的位移 *)6\V}`  
plot k6RVP:V  
history 1 lm
D[Cn  
;施加第一级荷载 N#t`ZC&m'  
ini xdis 0.0 ydis 0.0 zdis 0.0 n9`]}bnX  
apply szz -8.15e6   ran z 0 0.1 group pile   MtN!Xx  
solve woBx609Aak  
save app1.sav