运行这段命令的时候会出现如下的提示,请问各位大神是什么情况,该怎么解决啊?
new
config fluid
gen zone radcyl p0 0 0 0 p1 5 0 0 p2 0 20 0 p3 0 0 5 &
dim 3 3 3 3 size 2 10 4 2 ratio 1.0 1.0 1.0 1.2 fill &
group tunnel
gen z b p0 0 0 5 p1 5 0 5 p2 0 20 5 p3 0 0 15 &
rat 1 1 1 size 2 10 4
gen z b p0 5 0 0 p1 25 0 0 p2 5 20 0 p3 5 0 15 &
rat 1 1 1 size 8 10 6
group soil range group tunnel not
gen zone reflect normal 0 0 -1
fix x range x -.1 .1
fix x range x 24.9 25.1
fix y range y -.1 .1
fix y range y 19.9 20.1
fix z range z -14.9 -15.1
;cam clay
model cam-clay
prop shear 150000 bulk_bound 20e6
prop mm 0.73 lambda 0.0938 kappa 0.0234
prop mpc 0.395e6 mp1 1.0 mv_l 3.32
ini dens 1270
;fluid para
model fl_iso
prop perm 2.66e-13
ini fdens 1000
ini sat 1.0
set grav 0 0 -10
;ini stress and pp
ini szz -190500 grad 0 0 1.27e4
ini szz add -0.795e5 grad 0 0 0.53e4
ini sxx -128509.5 grad 0 0 0.857e4
ini sxx add -0.795e5 grad 0 0 0.53e4
ini syy -128509.5 grad 0 0 0.857e4
ini syy add -0.795e5 grad 0 0 0.53e4
ini pp 1.5e5 grad 0 0 -10000
fix pp 0 range z 14.9 15.1
;xiuzeng jianqiaomoxing
def camclay_ini_p
pnt = zone_head
loop while pnt # null
OCR=1.2
s1=-z_sxx(pnt)
s2=-z_syy(pnt)
s3=-z_szz(pnt)
p0=(s1+s2+s3)/3.0-z_pp(pnt)
z_prop(pnt,'cam_cp')=p0
q0=sqrt(((s1-s2)*(s1-s2)+(s2-s3)*(s2-s3)+(s3-s1)*(s3-s1))*0.5)
temp1=q0/(z_prop(pnt,'mm')*p0)
pc=p0*(1.0+temp1*temp1)*OCR
z_prop(pnt,'mpc')=pc
pnt=z_next(pnt)
end_loop
end
@camclay_ini_p
;chushi pingheng
set fl off
set mech on
solve
ini xdisp 0 ydisp 0 zdisp 0
save ini_equivalence
hist fltime
;kaiwamian parameter
def sup_stress
lumda=0.9
o_press=-208009.55
o_grad=1.39e4
s_press=o_press*lumda
s_grad=o_grad*lumda
end
@sup_stress
;first huan kaiwa parameter
def excate_step1
n=1
a1=2*n-2
a2=2*n
b1=2*n-0.01
b2=2*n+0.01
n1=16*n-1
n2=16*n
t=4*3600*n
end
@excate_step1
model null range cyl end1 0 a1 0 end2 0 a2 0 rad 3.0