我的命令流是
new
gen zone brick p0 0 0 0 p1 18 0 0 p2 0 5 0 p3 0 0 8 size 17 10 8 group 1
gen zone brick p0 18 0 0 p1 24 0 0 p2 18 5 0 p3 18 0 8 size 9 10 8 group 2
gen zone brick p0 0 0 8 p1 18 0 8 p2 0 5 8 p3 0 0 18 p4 18 5 8 p5 0 5 18 p6 8 0 18 p7 8 5 18 size 17 10 17 group 3
group 0 range group 1 any group 2 any group 3 any
plo blo gro
;====================================
;设置边界条件
fix x y z range z -0.1 0.1
fix x range x 23.9 24.1
fix x range x -0.1 0.1
fix y range y -0.1 0.1
fix y range y 4.9 5.1
;====================================
;初始地应力的生成
model elas
prop density 2000 bulk 3e9 shear 1e9 range group 0
set gravity 0 0 -10.0
solve
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zdisp 0
model mohr
prop density 2000 bulk 1.0e8 shear 3.0e7 coh 12100 tens 1.0e6 fric 20 dila 20 range group 0
plot cont ssi outline on
solve
save 初始应力.sav
; =======================================================
sel pile id=1 begin=(13 2.5 13) end=(13 2.5 4) nseg=12
sel pile prop Emod=8.0e10 Nu=0.30 XCArea=12.56 XCJ=0 XCIy=0 XCIz=0 Per=3.14 CS_sK=1.3e11 CS_sCoh=10e3 CS_sFric=10.0 &
CS_nK=1.3e11 CS_nCoh=10e3 CS_nFric= 20 CS_nGap=off
plo con ssi
solve
; =======================================================
sel cable id 1 beg 13 2.5 13 end 4 2.5 9 nseg 8
sel cable id 1 prop emod 200e9 xcarea 5.54e-3 ytens 15.34e5 gr_k 0 gr_coh 0 gr_fric 0 range cid 1 3
sel cable id 1 prop emod 200e9 xcarea 5.54e-3 ytens 30e5 gr_k 0.56e9 gr_coh 1.5e10 gr_fric 25 gr_per 0.264 range cid 4 8
sel delete link range id 1
sel link id 1 1 target node tgt_num 14
sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id 1
sel cable id 1 pretension 10e3 range cid 1 3 ;在自由段施加预应力
solve
问题是：（1）加抗滑桩之前和之后的位移以及应力等基本未发生变化。这是什么原因？
（2）模拟抗滑桩治理边坡时是否需要给抗滑桩施加一定的荷载？
（3）如何判断抗滑桩是否发生破坏？
