麻烦问下哪位做过抗滑桩啊，有命令流的发上来参考下。。
现在做一个抗滑桩不知道哪里有问题。。感觉很怪异
我把我自己做的命令流传上来大家帮我看下
new
ge zo b p0 0 0 0 p1 8 0 0 p2 0 32 0 p3 0 0 11 p4 8 32 0 p5 0 32 11 p6 8 0 11 p7 8 32 11 size 4 16 11 gro 01
ge zo b p0 -24 0 0 p1 0 0 0 p2 -24 32 0 p3 -24 0 32 p4 0 32 0 p5 -24 32 32 p6 0 0 11 p7 0 32 11 size 8 16 11 gro 02
gro 1 ran gro 01 any gro 02 any
ge zo b p0 -24 0 32 p1 0 0 11 p2 -24 32 32 p3 -24 0 48 p4 0 32 11 p5 -24 32 48 p6 0 0 27 p7 0 32 27 size 8 16 16 gro 2
fix z range z -0.1 0.1
fix y range y -0.1 0.1
fix y range y 31.9 32
fix x range x -24.1 -23.9
fix x range x 7.9 8.1
; 属性
model ela range z 0 48
prop density 2600 bulk 1.515e9 shear 7.813e8 coh 400e3 fric 38 range group 1
prop density 1900 bulk 1.818e7 shear 9.375e6 coh 14e3 fric 18 range group 2
set large
set grav 0,0,-9.81
solve
save 初始平衡.sav
rest 初始平衡.sav
ini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0
plot block group
rest 初始平衡.sav
ini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0
plot block group
;del ran x -3 0 y 10 12
;del ran x -3 0 y 22 24
gro 1_1 ran gro 1 x -3 0 y 10 12
gro 1_2 ran gro 1 x -3 0 y 22 24
gro 2_1 ran gro 2 x -3 0 y 10 12
gro 2_2 ran gro 2 x -3 0 y 22 24
ge zo b p0 -3 10 29.43 p1 0 10 27 p2 -3 12 29.43 p3 -3 10 32 p4 0 12 27 p5 -3 12 32 p6 0 10 32 p7 0 12 32 size 1 1 6 gro 00
ge zo b p0 -3 22 29.43 p1 0 22 27 p2 -3 24 29.43 p3 -3 22 32 p4 0 24 27 p5 -3 24 32 p6 0 22 32 p7 0 24 32 size 1 1 6 gro 00
gro 3 ran gro 00 any gro 2_1 any gro 2_2 any gro 1_1 any gro 1_2 any
ge merge 0.2
ge separate 1
int 1 wrap 1 2
Int 1 maxedge 1
ge separate 3
int 2 wrap 1 3
int 2 wrap 2 3
Int 2 maxedge 1
plot int yellow
int 1 prop ks 2.57e8 kn 2.57e8 coh 7e3 fric 9
int 2 prop ks 5.46e8 kn 5.46e8 coh 7e3 fric 9
model mo range group 1
model mo range group 2
prop density 2600 bulk 4e9 shear 9e8 coh 400e3 fric 38 range group 1
prop density 1900 bulk 5e7 shear 9.375e7 coh 14e3 fric 18 range group 2
model elas range group 3
prop density 2500 bulk 4.012e9 shear 1.082e9 range group 3
pl bl gro
pl add int yel
set large
solve
save 加桩平衡.sav
模型的思想：
（1）先求解初始地应力，然后将地应力加载在模型上，求解抗滑桩作用下滑体的最大水平位移：
问题：
（1）如何将初始地应力加载到模型上？
（2）为什么计算出来的滑体会顺着基岩一直下滑，导致滑体网格变形过大，无法计算下去。我个人认为是网格接触面设置有问题。
