关于二分法计算边坡稳定系数
在二分法计算边坡稳定系数的时候，会有这样的语句：“if mech_ratio<1.0e-5”，这个选择语句，当模型运行收敛时，mech_ratio会小于1e-5；但是当不收敛，运行的时候，mech_ratio会一直在某一个大于1e-5的数值附近徘徊，而不减小，这样程序则会一直运行下去？？？请教高手们，这样的问题该怎么解决？？？:):):)下面是相关的程序片段。
def SSR
ait1=0.02
floww=1.0
fup=2.0
n=0
ff=(floww+fup)/2
loop while (fup-floww)>ait1
flag=0
coh1=1e5/ff
fri1=(atan((tan(45*pi/180))/ff))*180/pi
command
model null
model mohr
pro bulk 2e8 she 1e8 fric fri1 coh coh1 ten 1e5
so ratio 9.8e-6
endcommand
if mech_ratio<1.0e-5
floww=ff
fup=fup
else
fup=ff
floww=floww
endif
ff=(floww+fup)/2
n=n+1
endloop
end
SSR
页:
[1]