我加抗滑桩后用自编强度折减法求安全系数，求得安全系数是
麻烦大家给我看看这个问题在哪？命令流如下：
impgrid 7.5m.Flac3D
gen merge 1.0e-7
group 抗滑桩 range group 2 ;抗滑桩
group 桩侧岩土体 range group 3 ;桩侧岩土体
group 桩下岩土体 range group 4 ;桩下岩土体
group 桩前后岩土体 range group 5 ;桩前后岩土体
plo blo group
; =======================================================建立桩土接触面
gen sep 抗滑桩
;桩侧接触面
int 1 wrap 抗滑桩 桩侧岩土体
;桩底接触面
int 2 wrap 抗滑桩 桩下岩土体
; =======================================================强度折减法
;自定义强度折减法
def SSR
step1=10000
ait1=0.001 ;精度
k11=0.0 ;lower bound
k12=2.0 ;upper bound
ks=(k11+k12)/2
loop while (k12-k11)>ait1
coh1=1.0e4/ks ; 粘结力(Pa)
fri1=(atan((tan(20.0*pi/180.0))/ks))*180.0/pi ; 内摩擦角(度)
ten1=0 ; 抗拉强度
grav0=-9.80 ; 重力加速度(N/kg)
dens1=2000 ; 密度(kg/m3)
dens2=2400 ; 密度(kg/m3)
K1=1.3e7 ; 岩土体体积模量(Pa)
K2=16.67e9 ; 抗滑桩体积模量(Pa)
G1=8e6 ; 岩土体剪切模量(Pa)
G2=12.50e9 ; 抗滑桩剪切模量(Pa)
;折减的实现过程
command
hist reset
model null
;初始地应力场的形成
model elas range group 桩前后岩土体 any group 桩侧岩土体 any group 桩下岩土体 any group 抗滑桩 any
pro bulk K1 she G1 dens dens1 range group 桩前后岩土体 any group 桩侧岩土体 any group 桩下岩土体 any
pro bulk K2 she G2 dens dens2 range group 抗滑桩
;接触面参数赋值
interface 1 prop Kn 4.7e8 Ks 4.7e8 coh 0.8e4 fric 16
interface 2 prop Kn 4.7e8 Ks 4.7e8 coh 0.8e4 fric 16
fix x y z range z 0
fix x range x 0
fix x range x 40
fix y range y 0
fix y range y 3
set grav 0 0 grav0
solve
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
;塑性阶段的求解
model mohr range group 桩前后岩土体 any group 桩侧岩土体 any group 桩下岩土体 any ;采用莫尔库仑模型
pro dens dens1 bulk K1 she G1 coh coh1 friction fri1 tens ten1 range group 桩前后岩土体 any group 桩侧岩土体 any group 桩下岩土体 any
hist unbal
set mech ratio 1e-5
solve step step1
endcommand
; =======================================================求解
;收敛
if mech_ratio<1.0e-5
k11=ks
k12=k12
else
k12=ks
k11=k11
endif
ks=(k11+k12)/2
endloop
;=====================================
;计算结果的保存
fosfile0='_fos'+'.sav'
command
save fosfile0
endcommand
end
;**********************************************
;程序执行及结果显示
SSR
pr ks
自编强度折减求安全系数.txt (3 K) 下载次数:91 7.5mslope.txt (5 K) 下载次数:75 
自编强度折减法求解.zip (4640 K) 下载次数:231
