文档库 最新最全的文档下载
当前位置:文档库 › kaiwa

kaiwa


block 0 0 140 0 140 -99 0 -99
round 0.5
crack 0 -50 140 -50
crack 0 -60 140 -60
crack 0 -70 140 -70
crack 0 -74 140 -74
crack 0 -79 140 -79
crack 0 -99 140 -99
gen edge 5.0
bound stress -1.02e5 0 -1.02e5 ygrad 1.34e4 0 1.34e4
insit stres -1.02e5 0 -1.02e5 szz -1.02e5 ygra 1.34e4 0 1.34e4 zgra 0 1.34e4
gravity 0 -10
change cons=3
prop mat=1 d=1000 s=5.0e7 co=6.00e6 dil=4 bul=1.5e7 fric=15
change mat=1 range 0 140 -50 0
prop mat=2 d=2300 s=8.0e8 co=3.5e7 dil=9 bul=4.5e7 fric=45
change mat=2 range 0 140 -60 -50
prop mat=3 d=2500 s=1.52e9 co=5e7 dil=10 bul=5.5e7 fric=70
change mat=3 range 0 140 -70 -60
prop mat=4 d=2200 s=4.00e8 co=2.0e7 dil=7 bul=3.0e7 fric=40
change mat=4 range 0 140 -74 -70
prop mat=5 d=1600 s=7.00e7 co=1.00e7 dil=8 bul=2.7e7 fric=27
change mat=5 range 0 140 -79 -74
prop mat=6 d=2000 s=8.00e8 co=3.9e7 dil=9 bul=4.5e7 fric=40
change mat=6 range 0 140 -99 -79
change jcons=2
prop jmat=1 jkn=1e9 jks=1e9 jcoh=1e10 jtens=1e10
solve
del 50 90 -79 -74
solve

材料性质
*material 1 粗砂岩-砾岩
unmark i 1 96 j 1 33
mark i 1 96j 1 4
prop d 2500 s 1.52e9 co 5e7 dil 10 bul 5.5e7 fric 70 &
ctab 1 dtab 2 ftab 3 i 1 95 j 1 4
table 1 0 5.000e7 0.001 3e7 0.1 3e7
table 2 0 10 0.001 5 0.1 5
table 3 0 70.000 0.001 40 0.1 40

*material 2 砂质泥岩
unmark i 1 96 j 1 33
mark i 1 96 j 4 5
prop d 2200 s 4.00e8 co 2.0e7 dil 7 bul 3.0e7 fric 40 &
ctab 4 dtab 5 ftab 6 i 1 95 j 4 5
table 4 0 2.00e7 0.001 3e7 0.1 3e7
table 5 0 7 0.001 4.0 0.1 4.0
table 6 0 40.000 0.001 16 0.1 16
*1,2,3,4,5,6 material 3 2-3煤
unmark i 1 96 j 1 33
mark i 1 96 j 5 10
prop d 1600 s 7.00e7 co 1.000e7 dil 8 bul 2.7e7 fric 27 &
ctab 7 dtab 8 ftab 9 i 1 95 j 5 11
table 7 0 1.000e7 0.001 5e6 0.1 5e6
table 8 0 8 0.001 4.0 0.1 4.0
table 9 0 27.000 0.001 14 0.1 14

*4material 3 2-3煤
*unmark i 1 95 j 1 32
* mark i 1 95 j 5 6
*prop d 1600 s 7.00e7 co 1.000e7 dil 8 bul 2.7e7 fric 27 &
* ctab 7 dtab 8 ftab 9 i 1 95 j 5 6
*table 7 0 1.000e7 0.001 5e6 0.1 5e6
*table 8 0 8 0.001 4.0 0.1 4.0
*table 9 0 27.000 0.001 14 0.1 14

*3material 3 2-3煤
*unmark i 1 95 j 1 32
* mark i 1 95 j 5 6
*prop d 1600 s 7.00e7 co 1.000e7 dil 8 bul 2.7e7 fric 27 &
* ctab 7 dtab 8 ftab 9 i 1 95 j 6 7
*table 7 0 1.000e7 0.001 5e6 0.1 5e6
*table 8 0 8 0.001 4.0 0.1 4.0
*table 9 0 27.000 0.001 14 0.1 14


*2material 3 2-3煤
*unmark i 1 95 j 1 32
* mark i 1 95 j 6 7
*prop d 1600 s 7.00e7 co 1.000e7 dil 8 bul 2.7e7 fric 27 &
* ctab 10 dtab 11 ftab 12 i 1 95 j 7 8
*table 10 0 1.000e7 0.001 5e6 0.1 5e6
*table 11 0 8 0.001 4.0 0.1 4.0
*table 12 0 27.000 0.001 14 0.1 14
*1material 3 2-3煤
*unmark i 1 95 j 1 32
* mark i 1 95 j 7 8
*prop d 1600 s 7.00e7 co 1.000e7 dil 8 bul 2.7e7 fric 27 &
* ctab 13 dtab 14 ftab 15 i 1 95 j 8 9
*table 13

0 1.000e7 0.001 5e6 0.1 5e6
*table 14 0 8 0.001 4.0 0.1 4.0
*table 15 0 27.000 0.001 14 0.1 14
*material 4 细砂岩
unmark i 1 96 j 1 33
mark i 1 96 j 10 14
prop d 2000 s 8.00e8 co 3.9e7 dil 9 bul 4.5e7 fric 40 &
ctab 16 dtab 17 ftab 18 i 1 95 j 11 15
table 16 0 3.900e7 0.001 1e7 0.1 1e7
table 17 0 9.0 0.001 4.0 0.1 4.0
table 18 0 40.000 0.001 22 0.1 22
*material 5 2-1煤层
unmark i 1 96 j 1 33
mark i 1 96 j 14 16
prop d 1600 s 7.00e7 co 1.000e7 dil 8 bul 2.7e7 fric 27 &
ctab 19 dtab 20 ftab 21 i 1 95 j 15 17
table 19 0 1.00e7 0.001 5e6 0.1 5e6
table 20 0 8 0.001 4.0 0.1 4.0
table 21 0 27.00 0.001 14 0.1 14
*material 6 泥岩
unmark i 1 96 j 1 33
mark i 1 96 j 16 19
prop d 1600 s 7e7 co 1.20e7 dil 8 bul 2.7e7 fric 25 &
ctab 22 dta 23 ftab 24 i 1 95 j 17 20
table 22 0 1.20e7 0.001 5e6 0.1 5e6
table 23 0 8 0.001 4.0 0.1 4.0
table 24 0 25.00 0.001 18 0.1 18
*material 7 1-1 1-2煤
unmark i 1 96 j 1 33
mark i 1 96 j 19 21
prop d 1600 s 7.00e7 co 1.000e7 dil 8 bul 2.7e7 fric 27 &
ctab 25 dta 26 ftab 27 i 1 95 j 20 22
table 25 0 1.00e7 0.001 5e6 0.1 5e6
table 26 0 8 0.001 4.0 0.1 4.0
table 27 0 27.00 0.001 13 0.1 13
*material 8 砂质泥岩
unmark i 1 96 j 1 33
mark i 1 96 j 21 27
prop d 2200 s 4.00e8 co 2.0e7 dil 7 bul 3.0e7 fric 35 &
ctab 28 dtab 29 ftab 30 i 1 95 j 22 28
table 28 0 2.00e6 0.001 1e6 0.1 1e6
table 29 0 7 0.001 4.0 0.1 4.0
table 30 0 35.000 0.001 22 0.1 22
*material 9 砂岩
unmark i 1 96 j 1 33
mark i 1 96 j 27 28
prop d 2300 s 8.00e8 co 3.5e7 dil 9 bul 4.5e7 fric 45 &
ctab 31 dtab 32 ftab 33 i 1 95 j 28 29
table 31 0 3.5000e7 0.001 10e6 0.1 10e6
table 32 0 9 0.001 4.0 0.1 4.0
table 33 0 45.000 0.001 22 0.1 22
*material 10 黄土层
unmark i 1 96 j 1 33
mark i 1 96 j 28 33
prop d 1000 s 5.0e7 co 6.00e6 dil 4 bul 1.5e7 fric 15 &
ctab 34 dtab 35 ftab 36 i 1 95 j 29 34
table 34 0 6.00e6 0.001 3e6 0.1 3e6
table 35 0 4 0.001 2.0 0.1 2.0
table 36 0 15.000 0.001 7 0.1 7
Tunnel Support Loading 3 - 7
; tunnel support loading
;block 0 -30 60 -30 60 -90 0 -90
round 0.1
crack 0 -70 60 -70
crack 30 0 30 -90
crack 42 0 42 -90
tun 30 -70 4.11 12
tun 42 -70 2.62 8
tun 30 -70 5.5 12
tun 42 -70 5.5 12
gen edge 2.0
; initial stress state
bound stress 1.02e5 0 1.02e5 ygrad 1.34e4 0 1.34e4
insit stres 1.02e5 0 1.02e5 szz 1.02e5 ygra 1.34e4 0 1.34e4 zgra 0 1.34e4
gravity 0 -10
; rock properties
prop mat=1 d=1340
zone model mohr
zone shear=.33e9 bulk=.99e9 coh=1e6 fric=30.0
;prop mat=1 d=1340 g=.33e9 k=.99e9 coh=1e6 fric=30.0
;change con=3
; elastic joint properties
prop jmat=1 jkn=1e9 jks=1e9 jcoh=1e10 jtens=1e10
; cycle to initial equilibrium
hist solve type 1
solve rat 1e-5
save tun1.sav
; excavate service tunnel
del 40 44 -72 -68
bound -1 1 -91 0 xvel=0

.0
bound -1 90 -91 -89 yvel=0.0
bound 59 61 -91 0 xvel=0.0
; histories around tunnel 1
hist ydis 42 -67 sxx 42 -67
hist ydis 42 -73 sxx 42 -73
hist xdis 39 -70 syy 39 -70
hist xdis 45 -70 syy 39 -70
; histories around tunnel 2
hist ydis 30 -65.0 sxx 30 -65.0
hist ydis 30 -75.0 sxx 30 -75.0
hist xdis 25.0 -70 syy 25.0 -70
hist xdis 35.0 -70 syy 35.0 -70
reset disp jdisp
solve rat 1e-5
save tun2.sav
; line service tunnel
struct gen xc=42 yc=-70 npoint=16 mat=5 thick=0.37 fang=-11.25 theta 360
prop mat=5 st_d=2400 st_ymod=24.0e9 st_prat=0.20 st_yield=1e10
prop mat=5 if_kn=1e8 if_ks=1e7 if_coh=1.0e10
; excavate main tunnel
del 28 32 -72 -68
reset disp
solve rat 1e-5
save tun3.sav
; line main tunnel
struc gen xc=30 yc=-70 npoint=8 mat=5 thick=0.46 fang=22.5 theta -360
; add additional load representing raised water level
bound stress 0.0 0.0 -1.0e6 range -1 91 -31 -29
; add hydrostatic loads to tunnel liners
struct apply press 0.0 2.06e6
reset disp
solve rat 1e-5
save tun4.sav
ret
开挖程序
round 0.1
block 0 0 0 100 150 100 150 0
crack 0 30 150 30
crack 0 40 150 40
crack 0 60 150 60
crack 70 30 70 33
crack 70 33 74 33
crack 74 33 74 30
crack 74 30 70 30
change cons=3
jregion id=1 0 0 0 30 150 30 150 0
jregion id=2 0 30 0 40 150 40 150 30
jregion id=3 0 40 0 60 150 60 150 40
jregion id=4 0 60 0 100 150 100 150 60
jset 0,0 150,0 0,0 6,0 range jreg 1
jset 90,0 6,0 6,0 15,0 (15 0) range jreg 1
jset 90,0 6,0 6,0 15,0 (7.5 6) range jreg 1
block being split
jset 0,0 150,0 0,0 2,0 range jreg 2
jset 90,0 2,0 2,0 4,0 (4 30) range jreg 2
jset 90,0 2,0 2,0 4,0 (2 32) range jreg 2
jset 0,0 150,0 0,0 5,0 range jreg 3
jset 90,0 5,0 5,0 10,0 (10 40) range jreg 3
jset 90,0 5,0 5,0 10,0 ( 5 45) range jreg 3
jset 0,0 150,0 0,0 10,0 range jreg 4
jset 90,0 10,0 10,0 15,0 ( 15 60) range jreg 4
jset 90,0 10,0 10,0 15,0 (7.5 70) range jreg 4
change mat=1 range 0 150 0 30
change mat=2 range 0 150 30 40
change mat=3 range 0 150 40 60
change mat=4 range 0 150 60 100
generate edge 5 range 0 150 0 30
generate edge 2 range 0 150 30 40
generate edge 4 range 0 150 40 60
generate edge 10 range 0 150 60 100
delete range 70 74 30 33
bound stress -10 0 -10
insitu stress -10e6 0 -10e6 szz-10e6
; rock blocks ; joints (no cohesion)
prop m=1 dens=2500 k=5e9 g=3e9
prop jm=1 jkn=10e10 jks=10e10 jfric=30
prop m=2 dens=1000 k=5e8 g=3e8
prop jm=2 jkn=10e9 jks=10e9 jfric=20
prop m=3 dens=1500 k=1e9 g=1e9
prop jm=3 jkn=10e10 jks=10e10 jfric=25
prop m=4 dens=2500 k=5e9 g=3e9
prop jm=4 jkn=10e10 jks=10e10 jfric=32
change jcons=5
set jcondf=5
prop jmat jkn=40000 jks=40000 jfric=35 jrfric=10 &
jcoh=2.0 jrescoh=0 &
jdil=6 zdil=4e-4
bound xvel 0 range -0.01

0.01 -0.01 100.01
bound xvel 0 range 149.99 150.01 -0.01 100.01
bound yvel 0 range -0.01 150.01 -0.01 0.01
;set gravity ; set in-situ stresses (total)
grav 0 -9.8
hist type=2
hist unbal solve_ratio
hist xdis=-40,100 ydis=-40,100
damp auto
solve

相关文档