水平荷载下桩的位移问题

哼歌而梦

小弟写了一段命令流，模拟承台受水平推力下桩的位移问题，承台和土用的实体单元，桩用的pile单元，计算求解后发现桩的位移不对，请教一下各位大牛是什么原因？代码如下：
new
def para
    mo=200e5
    bo=0.3
    bulk_1=mo/(3*(1-2*bo))
    shear_1=mo/(2*(1+bo))
end
para
gen zone brick p0=(0,0,0),p1=(8,0,0),p2=(0,8,0),p3=(0,0,8) size 16,16,16 group 1
gen zone brick p0=(3.5,3.5,8),p1=(4.5,3.5,8),p2=(3.5,4.5,8),p3=(3.5,3.5,10) &
               size 2,2,4 group 2
sav jianmo.sav
;--------------赋材料参数------------------------------------
mod elastic
pro bulk bulk_1 shear shear_1 density=2000 range group 1
pro bulk=13.9e9 she=10.4e9 density=2500 range group 2
;---------------边界条件--------------------------------
fix z ran z -0.1 0.1
fix x ran x -0.1 0.1
fix x ran x 7.9 8.1
fix y ran y -0.1 0.1
fix y ran y 7.9 8.1
ini sxx -2.20e5 grad=(0,0,2.0e4)
ini syy -1.32e5 grad=(0,0,1.2e4)
ini szz -0.88e5 grad=(0,0,8.0e3)
set gravity=(0.0,0.0,-10.0)
solve
;----------------加桩--------------------------------
sel pile id=1 begin=(4, 4, 8.5) end=(4, 4, 4.0) nseg=8
sel pile id=1 prop Emod=8.0e10 Nu=0.30 XCArea=0.7854 &
        XCJ=9.82e-2 XCIy=4.91e-2 XCIz=4.91e-2 &
        per=3.14 &
        cS_sK=1.3e11 CS_sCoh=1.0e10 CS_sFric=0.0 &
        cS_nK=1.3e09 CS_nCoh=1.0e04 CS_nFric=0.0 &
        cS_nGap=off
sel delete link range id 9
sel link id 100 2 target zone
sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid &
yrdir=rigid zrdir=rigid range id 100
sel delete link range id 1
sel link id 101 1 target zone
sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid &
yrdir=rigid zrdir=rigid range id 101
app nstress 10e5 ran x 3.49 3.51 y 3.5 4.5 z 8 10

solve ratio=1e-5
save 22.sav
计算后的位移矢量图

桩的位移

显示的桩单元的位移数据

yinwen0809

期待高手解答

awk520

谢谢分享！学习后再说！