flac3D隧道模型例题命令流TXT

; main tunnel
gen zon radcyl p0 15 0 0 p1 23 0 0 p2 15 50 0 p3 15 0 8 &
size 4 10 6 4 dim 4 4 4 4 rat 1 1 1 1 fill
gen zon reflect dip 90 dd 90 orig 15 0 0
gen zon reflect dip 0 ori 0 0 0
; service tunnel
gp id 3000 (2.969848,0.0,-0.575736)
gp id 3001 (2.969848,50.0,-0.575736)
gen zon radcyl p0 0 0 -1 p1 7 0 0 p2 0 50 -1 p3 0 0 8 p4 7 50 0 &
p5 0 50 8 p6 7 0 8 p7 7 50 8 p8 gp 3000 p10 gp 3001 &
size 3 10 6 4 dim 3 3 3 3 rat 1 1 1 1 fill
gen zon radcyl p0 0 0 -1 p1 0 0 -8 p2 0 50 -1 p3 7 0 0 p4 0 50 -8 &
p5 7 50 0 p6 7 0 -8 p7 7 50 -8 p9 gp 3000 p11 gp 3001 &
size 3 10 6 4 dim 3 3 3 3 rat 1 1 1 1 fill
; outer boundary
gen zone radtun p0 7 0 0 p1 50 0 0 p2 7 50 0 p3 15 0 50 p4 50 50 0 &
p5 15 50 50 p6 50 0 50 p7 50 50 50 &
p8 23 0 0 p9 7 0 8 p10 23 50 0 p11 7 50 8 &
size 6 10 3 10 rat 1 1 1 1.1

gen zone brick p0 0 0 8 p1 7 0 8 p2 0 50 8 p3 0 0 50 &
p4 7 50 8 p5 0 50 50 p6 15 0 50 p7 15 50 50 &
size 3 10 10 rat 1 1 1.1
gen zon reflect dip 0 ori 0 0 0 range x 0 23 y 0 50 z 8 50
gen zon reflect dip 0 ori 0 0 0 range x 23 50 y 0 50 z 0 50
group zone main1 range cyl end1 15 0 0 end2 15 25 0 rad 4
group zone main2 range cyl end1 15 25 0 end2 15 50 0 rad 4
group zone service range cyl end1 0.0,0.0,-0.575736 &
end2 0.0,50.0,-0.575736 rad 3.0
save tun0

rest tun0
; mohr-coulomb model
model mech mohr
prop shear 0.36e9 bulk 0.6e9 coh 1e5 fric 20 tens 1e5
; boundary and initial conditions
apply szz -1.4e6 range z 49.9 50.1
fix z range z -50.1 -49.1
fix x range x -.1 .1
fix x range x 49.9 50.1
fix y range y -.1 .1
fix y range y 49.9 50.1
ini sxx -1.4e6 syy -1.4e6 szz -1.4e6

hist add unbal
hist add gp xdis 3,0,-1
hist add gp zdis 0,0,2
hist add gp xdis 3,25,-1
hist add gp zdis 0,25,2

solve
save tun1


; exacavate 25 m section of service tunnel
rest tun1
model mech null range group service y 0,25
solve
save tun2

restore tun2
ini xdis 0.0 ydis 0.0 zdis 0.0
hist purge
;
model mech null range group service y 25 50
; liner
sel shell id=1 range cyl end1 0 0 -1 end2 0 25 -1 rad 3
sel shell prop iso=(25.3e9, 0.266) thick = 0.5
sel node fix y xr zr range y -0.1 0.1 ; symmetry cond.
sel node fix x yr zr range x -0.1 0.1 ; symmetry cond.
;
model mech mohr range group service y 25 50
prop shear 0.36e9 bulk 0.6e9 coh 1e5 fric 20 tens 1e5
;
; excavate main tunnel - main1
model mech null range group main1
solve
sel recover surface surfx -1 0 0
sel recover stress depth_fac 1.0
save tun3