- 积分
- 0
- 注册时间
- 2016-10-7
- 仿真币
-
- 最后登录
- 1970-1-1
|
很多例子中先进行的开挖后计算平衡,这样是不是存在问题?陈育民老师书中关于修正剑桥模型中,为什么照他的代码运行,却出现了*** 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
|
|