希望各位高手看看下面的程序什么地方不对，加桩后计算边坡安全系数出问题，恳请给予指出和修改，小弟感激不尽！n
;===================================
;建立网格模型
gen zone brick &
p0 0 0 0 p1 15 0 0 p2 0 6 0 p3 0 0 20 size 10 6 16 ratio 1 1 1.03
gen zone brick &
p0 15 0 0 p1 50 0 0 p2 15 6 0 p3 15 0 20 &
size 16 6 16 ratio 1.03 1 1.03
gen zone brick &
p0 15 0 20 p1 50 0 20 p2 15 6 20 p3 25 0 30 &
p4 50 6 20 p5 25 6 30 p6 50 0 30 &
p7 50 6 30 size 16 6 8 ratio 1.03 1 1
;====================================
plot block group
pause
;设置边界条件
fix x y z range z -0.1 0.1
fix x range x 49.9 50.1
fix x range x -0.1 0.1
fix y
;====================================
;初始地应力的生成
model elas
prop density 2000 bulk 3e9 shear 1e9
set gravity 0 0 -10.0
solve
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
;===================================
;安全系数求解
model mohr
prop density 2000 bulk 1.0e8 shear 3.0e7 &
coh 12380.0 tens 1.0e6 fric 20 dila 20
;=================================================================
;建立桩单元
sel pile id=1 begin=(15 3 20) end=(15 3 4) nseg=16
sel pile id=1 prop emod=8E10 nu=0.3 xcarea=6 &
xcj=1.94E2 xciy=1.62E2 xciz=3.2E1 &
per=10 &
cs_sk=1.3E11 cs_scoh=1.0E10 cs_sfric=0.0 &
cs_nk=1.3E9 cs_ncoh=1.0E4 cs_nfric=0.0 &
cs_ngap=off
hist id=1 unbal
hist id=2 sel pilesel moment mz (19 3 16)
hist id=3 sel pilesel force fy end1 cid 1
hist id=4 gp xdisp 19 3 24
solve fos file slope3dfos12-3.sav associated
