FLAC3D流固耦合中隧道开挖的问题？

脱缰的土狗

很多例子中先进行的开挖后计算平衡，这样是不是存在问题？陈育民老师书中关于修正剑桥模型中，为什么照他的代码运行，却出现了*** Cam-clay: Yield envelope cannot be reached？求大神帮部啊！另附代码：
new
set fish safe_conversion off
config fluid
gen zon 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 rat 1 1 1 1.2  group soil fill group tunnel
      
gen zone brick p0 0 0 5 p1 5 0 5 p2 0 20 5 p3 0 0 15 size 2 10 4 ratio 1 1 1 group soil
gen zone brick p0 5 0 0  p1 25 0 0 p2 5 20 0 p3 5 0 15 size 8 10 6 ratio 1 1 1 group soil

gen zon reflect normal 0 0 -1  

;plot create displ
;plot set back white
;plot set rot 0 0 40
;plot set mag 1
;plot add axes red
;plot add contour disp out on
;plor show

fix z range z -14.9,-15.1
fix x range x -0.1 0.1
fix x range x 24.9,25.1
fix y range y -0.1 0.1
fix y range y 19.9,20.1

model cam-clay

prop bulk_bound=20e6 shear=1.5e5  
prop mm 0.73 lambda 0.0938 kappa 0.0234
prop mpc 0.395e6 mp1 1.0 mv_l 3.32
ini dens 1270


model fl_iso
prop perm 2.66e-13
ini fdensity 1e3
ini fmod 2e9
ini sat 1.0
set gravity 0 0 -10

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 -1e4
fix pp 0 range z 14.9 15.1

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)
endloop
end
camclay_ini_p

set fl off mech on
solve
ini xdis 0 ydis 0 zdis 0
save ini.sav

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

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
sel shell id =1 range cyl end1 0 a1 0 end2 0 a2 0 rad 3.0 group soil
sel node local xdir 1 0 0 ydir 0 -1 0 range x -0.1 0.1
sel node fix lsys                     range x -0.1 0.1
sel node fix x yr zr                  range x -0.1 0.1
sel node local xdir 1 0 -1 ydir 0 -1 0 range y -0.1 0.1
sel node fix lsys                     range y -0.1 0.1
sel node fix y xr zr                  range y -0.1 0.1
sel shell id=1 prop iso=(3.45e10,0.3)thick=0.35 dens 2450
sel shell apply press -160e3 range cid n1 n2

apply nstress 187208.595 grad 0 0 1.251e4 range cyl end1 0 b1 0 end2 0 b2 0 rad 3.0

set fl off
solve
set fl on
solve age t

zhuqiankun2016

:D:D:D:D:D:D:D:D:D

chengdu815cn

谢谢分享。