plot set rotation 30 0 30 plot block group plot grid vel
4-13不考虑重力影响的均匀应力
new ;set grav 0 0 0
gen zone brick size 6 6 6 model elas
ini sxx=-5e6 syy=-1e7 szz=-2e7 apply sxx=-5e6 range x -0.1 0.1 apply sxx=-5e6 range x 5.9 6.1 apply syy=-1e7 range y -0.1 0.1 apply syy=-1e7 range y 5.9 6.1 apply szz=-2e7 range z -0.1 0.1 apply szz=-2e7 range z 5.9 6.1 prop bulk 8e9 shear 5e9 ini dens 2000 fix x range z -.1 .1 step 10
plot set rotation 30 0 30 ;网格三维显示 plot cont zdisp outl on
4-14考虑应力梯度的均匀材料
new
gen zone brick size 10 10 10 p1 20,0,0 p2 0,20,0 p3 0,0,20 model mohr
prop bulk 5e9 shear 3e9 fric 35 ini density 2500 set gravity 0,0,-10
fix x y z range z -0.1 0.1 ini szz -5.0e6 grad 0,0,2.5e4 ini syy -2.5e6 grad 0,0,1.25e4 ini sxx -2.5e6 grad 0,0,1.25e4 apply szz -4.5e6 range z 19.9 20.1 apply szz -5.0e6 range z -0.1 0.1
apply sxx -2.5e6 grad 0,0,1.25e4 range x -0.1 0.1 apply sxx -2.5e6 grad 0,0,1.25e4 range x 19.9 20.1 apply syy -2.5e6 grad 0,0,1.25e4 range y -0.1 0.1 apply syy -2.5e6 grad 0,0,1.25e4 range y 19.9 20.1 step 10
;网格三维显示
plot set rotation 30 0 30
plot cont zdisp outl on
4-15考虑应力梯度的非均匀材料
new
gen zone brick size 10 10 10 &
p0 0,-25,0 p1 20,-25,0 p2 0,0,0 p3 0,-25,20 model elas
prop bulk 5e9 shear 3e9
ini density 1600 range y -10,0 ini density 2000 range y -15,-10 ini density 2200 range y -25,-15 set gravity 0,-10,0 fix x range x -.1 .1 fix x range x 19.9 20.1 fix z range z -.1 .1 fix z range z 19.9 20.1 fix y range y -25.1 -24.9
ini syy 0.0 grad 0,1.6e4,0 range y -10,0 ini syy 4e4 grad 0,2.0e4,0 range y -15,-10 ini syy 7e4 grad 0,2.2e4,0 range y -15,-25 step 10
;网格三维显示
plot set rotation 30 0 30 plot cont zdisp outl on
4-16非均匀网格应力初始化
new
gen zone radcyl size 3 8 4 5 fill p1 10,0,0 p2 0,10,0 p3 0,0,10 mode elastic
prop shear 3e8 bulk 5e8 fix x range x -.1 .1 fix x range x 9.9 10.1 fix y range y -.1 .1 fix y range y 9.9 10.1 fix z range z -.1 .1
ini szz = -2.5e5 grad 0,0,2.5e4 ini density 2500 set grav 0,0,-10 step 10
;网格三维显示
plot set rotation 30 0 30 plot cont zdisp outl on
4-17不规则自由面应力初始化
gen zone brick size 15 15 10 p0 0,0,0 edge=100.0 model elastic
prop shear 3e8 bulk 5e8 def mountain gp = gp_head
loop while gp # null
zz = sqrt(gp_xpos(gp)?2 + gp_ypos(gp)?2)
dz = 0.06 * sin(0.2 * zz + 100.0) ; Sum Fourier terms for dz = dz + 0.06 * sin(0.22 * zz - 20.3) ; quasi-random surface dz = dz - 0.04 * sin(0.33 * zz + 33.3) ; topology. gp_zpos(gp) = 0.5 * gp_zpos(gp) * (1.0 + dz) gp = gp_next(gp) end_loop end
mountain
fix x range x -.1 .1
fix x range x 99.9 100.1 fix y range y -.1 .1
fix y range y 99.9 100.1 fix z range z -.1 .1 set grav 0,0,-10 ini density=2000
ini szz=-2.0e6 (grad 0,0,2.0e4) sxx=-4.0e6 (grad 0,0,4.0e4) syy=-4.0e6 (grad 0,0,4.0e4) step 100
plot set rotation 30 0 30
plot cont zdisp outl on;网格三维显示
4-18非均网格内部压实
new
gen zone brick size 8 8 10 ratio 1.2 1 1 model mohr ini dens 2000
prop bulk 2e8 shear 1e8 prop fric 30
fix x range x -.1 .1 fix x range x 7.9 8.1 fix y range y -.1 .1 fix y range y 7.9 8.1 fix z range z -.1 .1
;ini szz -2.0e5 grad 0,0,2e4 ;ini sxx -1.5e5 grad 0,0,1.5e4 ;ini syy -1.5e5 grad 0,0,1.5e4
set grav 10 step 1000 ;pause
;prop tens 1e10 coh 1e10 ;step 750
;prop tens 0 coh 0 ;step 250
;网格三维显示
plot set rotation 30 0 30 plot cont zdisp outl on
4-19模型改变后初始应力变化
new
gen zone brick size 5 5 5 model elas
prop sh 2e8 bu 3e8 fix x y z range z -.1 .1 set grav 0 0 -10 ini dens 2000 Solve
;网格三维显示
plot set rotation 30 0 30 plot cont zdisp outl on ;模型改变
model mohr range x 0 2 y 0 5 z 0 2
prop sh 2e8 bu 3e8 fric 35 range x 0 2 y 0 5 z 0 2 Step 2000
;网格三维显示
plot set rotation 30 0 30 plot cont zdisp outl on
4-20应力与孔隙压力的初始化
new
config fluid
gen zone brick size 8 5 10 model elas model fl_iso ini dens 2000
prop bulk 1e9 shear 5e8 prop poros 0.5 perm 1e-10 ini fmod 2e9 ini fdensity 1e3 ini sat 0
ini sat 1 range z -.1 5.1
set grav 0 0 -10 fix x range x -.1 .1 fix x range x 7.9 8.1 fix y range y -.1 .1 fix y range y 4.9 5.1 fix z range z -.1 .1
ini pp 5.e4 grad 0,0,-1.e4 range z 0.0 5. ini szz -20e4 grad 0,0,20e3
ini szz add -1.5e4 grad 0,0,.25e4 range z 5,6 ini szz add -2.5e4 grad 0,0,.5e4 range z 0,5 ini szz add -.25e4 range z 0,5 solve
;网格三维显示
plot set rotation 30 0 30 plot cont zdisp outl on
new
gen zone brick size 1 1 10 model elas
ini dens 2500 range z 0 5 ini dens 2250 range z 5 6 ini dens 2000 range z 6 10 prop bulk 1e9 shear 5e8 set grav 10
water dens 1000
water table ori 0 0 5 normal 0 0 1 fix x range x -.1 .1 fix x range x 7.9 8.1 fix y range y -.1 .1 fix y range y 4.9 5.1 fix z range z -.1 .1
ini szz -20e4 grad 0,0,20e3
ini szz add -1.5e4 grad 0,0,.25e4 range z 5,6 ini szz add -2.5e4 grad 0,0,.5e4 range z 0,5 ini szz add -.25e4 range z 0,5 solve
;网格三维显示
plot set rotation 30 0 30 plot cont zdisp outl on
4-21加载顺序
new
; 主隧道建模