new
gen zone brick &
p0 0 0 0 p1 30 0 0 p2 0 5 0 p3 0 0 20 size 3 1 3 rat 0.9 1 0.9
gen zone brick &
p0 30 0 0 p1 150 0 0 p2 30 5 0 p3 30 0 20 &
size 17 1 3 rat 1.05 1 0.9
gen zone brick &
p0 30 0 20 p1 150 0 20 p2 30 5 20 p3 90 0 80 &
p4 150 5 20 p5 90 5 80 p6 150 0 80 &
p7 150 5 80 size 17 1 17 rat 1.05 1 1
;===============================================
;*************强度折减法*******************
;===============================================
def SSR
;------------------
ait1=0.02 ;精度
k11=1.0 ;上限值
k12=2.0 ;下限值
ks=(k11+k12)/2
loop while (k12-k11)>ait1
coh1=12380/ks
fril=(atan((tan(20*pi/180))/ks))*180/pi
dila1=20.0
ten1=1e6
grav0=-10
dens1=2000
K1=1e8
G1=3e7
;==================
command
model null ;挖空
model elas
pro bu 1e10 sh 3e9 dens dens1
fix x y z range z -0.1 0.1
fix x range x -0.1 0.1
fix x range x 149.9 15.1
fix y
set grav 0 0 grav0
solve
ini xd 0 yd 0 zd 0
ini xv 0 yv 0 zv 0
mo mo ;采用莫尔库仑模型
pro bulk K1 she G1 dens dens1 coh coh1 &
friction fri1 dil dila1 tens ten1
set mech ratio 9.8e-6
solve step 10000
endcommand
;==================
;收敛(二分法的实现过程）
if mech_ratio<1.0e-5
k11=ks
k12=k12
else
k12=ks
k11=k11
endif
ks=(k11+k12)/2
endloop
fosfile0=string(size1)+'_'+string(size2)+'_fos'+'.sav'
command
save fosfile0
endcommand
end
SSR
pr ks

运行之后会出现这种情况，这怎么回事

是参数转换错误吗

大神在哪里

fril改成fri1

你用的是强度折减法吧