使用自编强度折减法计算结果为Fs=3.33e8

accbd

使用自编强度折减法，对澳大利亚ACADS ex1(a)算例进行校核，最后结果竟为Fs=3.333e8。问题出现在哪？

ACADS ex1(a):

c(kPa)	phi	γ(kN/m^3)	E(kPa)	泊松比
3.0	19.6	20.0	1.0e4	0.25

以下为命令流：
n
gen zone brick p0 0,0,0   p1 50,0,0    p2 0,0.5,0   p3 0,0,5 size 50 1 5
gen zone brick p0 10,0,5  p1 50,0,5    p2 10,0.5,5  p3 30,0,15   p4 50,0.5,5  &
     p5 30,0.5,15 p6 50,0,15   p7 50,0.5,15 size 40 1 10
plo sur

def SSR
ait1=0.002
k1=0.8
k2=1.2
kr=(k1+k2)/2
loop while (k2-k1)>ait1
co0=3000/kr
fri0=(atan((tan(19.60*pi/180))/kr))*180/pi
dens0=(20.0/10.0)*1000                ;重度γ=20.0,算出密度ρ
E0=1.0e7                                       ;弹模E=1.0e7
Poss0=0.25                                    ;泊松比μ=0.25
G0=E0/(2.0*(1.0+Poss0))  
K0=E0/(3.0*(1.0-2.0*Poss0))
G1=G0*100
K1=K0*100

command
model null

model elas
prop bu K0 sh G0 dens dens0 co co0 fri fri0

fix x range x -0.1 0.1
fix x range x 49.9 50.1
fix x y z range z -0.1 0.1
fix y
set grav 0 0 -10.0
plo con ssi ou on
solve
ini xdis 0 ydis 0 zdis 0 xvel 0 yvel 0 zvel 0

model mohr
prop bu K0 sh G0 dens dens0 co co0 fri fri0
set mech ratio 1.0e-5
solve step 20000
endcommand

if mech_ratio<1.0e-5
k1=kr
k2=k2
else
k2=kr
k1=k1
endif
kr=(k1+k2)/2
endloop

end

SSR
print kr

accbd

该校核题所给的安全系数为1.00，或者在其左右微小浮动。
我利用这个程序段算其他的算例可以成立，但是怎么用到这个上面就不行了？
出现的Kr=3.3333333338833e+008，安全系数不至于这么大啊？！

问题出现在哪里？亟待解决！多谢！

accbd

顶起！希望高人给指点一下！

副热带

期待高人指点一下咯

zhaoti66

设置力不平衡比率时有误

lsj821687340

根据陈育民老师书上273页（第二版）初始应力场的生成采用的体积模型和剪切模量的参数取得都比较大，而不是材料的原始的参数，是不是这个有问题。同时，孙书伟老师书上226页直接采用摩尔库伦模型生成应力场，不知道这两个哪个更合适。

lincolnlee

lsj821687340 发表于 2018-9-1 17:28
根据陈育民老师书上273页（第二版）初始应力场的生成采用的体积模型和剪切模量的参数取得都比较大，而不是 ...

我也在看这两本书的这块内容。陈老师书的初始应力那一章讲到了生成初始应力场的三种方法。兄弟可以看看，希望能一起交流啊。