3DEC算例1

3DEC算例1

poly brick 0,20 0,20 0,20

pl bl

pl reset

pl set dip 70 dd 210

gen edge 4

prop mat=1 dens 2500 bulk 5e9 g 3e9 phi 35

change cons 2

gravity 0 0 -10

bound stress -2.75e6,-2.75e6,-5.5e6 0,0,0 zgrad 1.25e4,1.25e4,2.5e4 0,0,0 range x 0

bound stress -2.75e6,-2.75e6,-5.5e6 0,0,0 zgrad 1.25e4,1.25e4,2.5e4 0,0,0 range x 20

bound stress -2.75e6,-2.75e6,-5.5e6 0,0,0 zgrad 1.25e4,1.25e4,2.5e4 0,0,0 range y 0

bound stress -2.75e6,-2.75e6,-5.5e6 0,0,0 zgrad 1.25e4,1.25e4,2.5e4 0,0,0 range y 20

bound stress 0,0,-5e6 0,0,0 range z 20

bound zvel=0 range z 0

insitu stress -2.75e6,-2.75e6,-5.5e6 0,0,0 zgrad 1.25e4,1.25e4,2.5e4 0,0,0

step 500

plot block fill off stress colorby maximum reverse

ret

！！！！！！！！！！！！！

new

poly brick -1.5 1.5 -1.5 1.5 -1.5 1.5

pl bl color by material

3DEC算例1

pl set dip 70 dd 200

tunnel radial region 1 &

a (-0.30 -1.5 0)(-0.3 -1.5 0.4)(-0.25 -1.5 0.47) &

(-0.15 -1.5 0.52)( 0 -1.5 0.55)( 0.15 -1.5 0.52) &

( 0.25 -1.5 0.47)( 0.3 -1.5 0.4)( 0.3 -1.5 0) &

b (-0.30 1.5 0)(-0.3 1.5 0.4)(-0.25 1.5 0.47) &

(-0.15 1.5 0.52)( 0 1.5 0.55)( 0.15 1.5 0.52) &

( 0.25 1.5 0.47)( 0.3 1.5 0.4)( 0.3 1.5 0)

remove region 1

ret

！！！！！！！！！！！！

new

poly brick -1.5 1.5 -1.5 1.5 -1.5 1.5

pl bl colorby material

pl reset

pl set dip 70 dd 200

tunnel radial region 1 &

a (-0.30 -1.5 0)(-0.3 -1.5 0.4)(-0.25 -1.5 0.47) &

(-0.15 -1.5 0.52)( 0 -1.5 0.55)( 0.15 -1.5 0.52) &

( 0.25 -1.5 0.47)( 0.3 -1.5 0.4)( 0.3 -1.5 0) &

b (-0.30 1.5 0)(-0.3 1.5 0.4)(-0.25 1.5 0.47) &

(-0.15 1.5 0.52)( 0 1.5 0.55)( 0.15 1.5 0.52) &

( 0.25 1.5 0.47)( 0.3 1.5 0.4)( 0.3 1.5 0)

remove region 1

ret

3DEC算例1

！！！！！！！

new

poly brick 0 10 0 10 0 10

plot block axes

pl reset

plot set dip 70 dd 150

jset dip 0 dd 180 origin 0 0 5

gen quad ndiv 4 4 4

change cons 1

prop mat=1 dens 1000 bulk 8e9 g 5e9

prop jmat=1 jkn 1e10 jks 1e10

gravity 0 0 -10

bound stress 0 0 -1e5 0 0 0 range x 0 10 y 0 10 z -.1 .1

bound xvel=0 range x -.1 .1

bound xvel=0 range x 9.9 10.0

bound yvel=0 range y -.1 .1

bound yvel=0 range y 9.9 10.1

insitu stress 0 0 -1e5 0 0 0 zgrad 0 0 1e4 0 0 0

hist zdisp 5 5 2.5

step 300

excavate range x 0 10 y 0 10 z 5 10

step 100

3DEC算例1

plot block colorby mat fill off axes

plot add vel line color red

ret

！！！！！！！！！！！！

step 300

excavate range x 0 10 y 0 10 z 5 10

step 100

plot block colorby mat fill off axes

plot add vel line color red

ret

！！！！！！！！！！

new

poly brick 0 10 0 10 0 10

plot block

plot reset

gen edge 2.0

prop mat=1 dens 1000 bulk 8e9 g 5e9

bound stress -5e6 0 0 0 0 0 range x 0

bound stress -5e6 0 0 0 0 0 range x 10

bound stress 0 -5e6 0 0 0 0 range y 0

bound stress 0 -5e6 0 0 0 0 range y 10

bound stress 0 0 -1e7 0 0 0 range z 0

3DEC算例1

bound stress 0 0 -1e7 0 0 0 range z 10

insitu stress -5e6 -5e6 -1e7 0 0 0

step 1

plot block fill off stress colorby maximum

ret

new

poly brick 0 20 0 20 0 20

pl bl

pl reset

gen edge 4.0

prop mat=1 dens 2500 bulk 5e9 g 3e9 phi 35

change cons 2

gravity 0 0 -10

bound stress -2.75e6 -2.75e6 -5.5e6 0 0 0 &

zgrad 1.25e4 1.25e4 2.5e4 0 0 0 range x 0

bound stress -2.75e6 -2.75e6 -5.5e6 0 0 0 &

zgrad 1.25e4 1.25e4 2.5e4 0 0 0 range x 20

bound stress -2.75e6 -2.75e6 -5.5e6 0 0 0 &

zgrad 1.25e4 1.25e4 2.5e4 0 0 0 range y 0

bound stress -2.75e6 -2.75e6 -5.5e6 0 0 0 &

zgrad 1.25e4 1.25e4 2.5e4 0 0 0 range y 20

bound stress 0 0 -5e6 0 0 0 range z 20

bound zvel=0 range z 0

insitu stress -2.75e5 -2.75e6 -5.5e6 0 0 0 &

zgrad 1.25e4 1.25e4 2.5e4 0 0 0

3DEC算例1

step 500

plot block fill off stress colorby maximu

ret

！！！！！！！！

new

poly brick 0 20 0 20 0 25

pl bl

pl reset

jset dip 0 origin 0 0 10

jset dip 0 origin 0 0 15

gen edge 2

change mat 1 range z 0 10

change mat 2 range z 10 15

change mat 3 range z 15 25

prop mat=1 dens 1600 bulk 5e9 g 3e9

prop mat=2 dens 2000 bulk 5e9 g 3e9

prop mat=3 dens 2200 bulk 5e9 g 3e9

change jmat 1

prop jmat 1 jkn 1e10 jks 1e10 coh 1e10

gravity 0 0 -10

insitu stress 0 0 -4.8e5 0 0 0 &

zgrad 0 0 2.2e4 0 0 0 range z 0 10

3DEC算例1

insitu stress 0 0 -4.6e5 0 0 0 &

zgrad 0 0 2e4 0 0 0 range z 10 15

insitu stress 0 0 -4e5 0 0 0 &

zgrad 0 0 1.6e4 0 0 0 range z 15 25

bound xvel=0 range x 0

bound xvel=0 range x 20

bound yvel=0 range y 0

bound yvel=0 range y 20

bound zvel=0 range z 0

hist unbal

step 500

plot set dip 90 dd 180 cent 10 11 12.5

plot x bcont szz minimum -6e5 maximum 0 interval 1.25e5 wireframe off

本文来源：https://www.bwwdw.com/article/ugti.html