接触面应力分布求助

aaronhoog

刚性桩水平位移引起的接触面正应力分布，但是这个程序执行完好像桩的位移没有变化，而且接触面正应力分布貌似也没有变啊，大侠帮帮小弟，感激不尽。

new
;----generate right top
gen zone radcyl p0 (0,0,0) p1 (15,0,0) p2 (0,0,-5.08) p3 (0,15,0) &
p4 (15,0,-5.08) p5 (0,15,-5.08) p6 (15,15,0) p7 (15,15,-5.08) &
p8 (1.27,0,0) p9 (0,1.27,0) p10 (1.27,0,-5.08) p11 (0,1.27,-5.08) &
size 3 10 6 15 ratio 1 1 1 1.15
;------generate right bottom
gen zone radcyl p0 (0,0,-5.08) p1 (15,0,-5.08) p2 (0,0,-7.62) p3 (0,15,-5.08) &
p4 (15,0,-7.62) p5 (0,15,-7.62) p6 (15,15,-5.08) p7 (15,15,-7.62) &
p8 (1.27,0,-5.08) p9 (0,1.27,-5.08) p10 (1.27,0,-7.62) p11 (0,1.27,-7.62) &
size 3 6 6 15 ratio 1 1 1 1.15 fill
gen zone reflect dd 270 dip 90 ; gennerate left top and left bottom
group soil
;------define interface
interface 1 face range cylinder end1 (0,0,0) end2 (0,0,-5.18) radius 1.37 &
cylinder end1 (0,0,0) end2 (0,0,-5.18) radius 1.17 not
interface 2 face range cylinder end1 (0,0,-4.98) end2 (0,0,-5.18) radius 1.27
;-----generate pile
gen zone cyl p0 (0,0,7) p1 (1.27,0,7) p2 (0,0,1) p3 (0,1.27,7) &
p4 (1.27,0,1) p5 (0,1.27,1) &
size 3 10 7
;gen zone cyl p0 (0,0,6.1) p1 (1.27,0,6.1) p2 (0,0,6) p3 (0,1.27,6.1) &
;p4 (1.27,0,6) p5 (0,1.27,6) &
;size 3 1 6
gen zone reflect dd 270 dip 90 range z 1 7
group pile range z 1 7
ini z add -6.08 range group pile
;-----define properties
model mohr range group soil
prop bulk 8.333e7 shear 3.846e7 coh 30000 fric 0 range group soil
model elas range group pile
prop bulk 8.333e77 shear 3.846e77 range group pile
interface 1 prop kn 1e8 ks 1e8 fric 20 coh 30000
interface 2 prop kn 1e8 ks 1e8 fric 20 coh 30000
ini dens 1230 range group soil
ini dens 0 range group pile
;model null range z -0.1 0.15
;-----------apply boundary
fix z range z -7.72 -7.52
fix x range x -15.1 -14.9
fix x range x 14.9 15.1
fix y range y -.1 .1
fix y range y 14.9 15.1
set grav 0 0 -10
solve rat 1.e-4
;---------apply disp
ini xdis 1 range group intruder
cyc 10000