3DEC经典算例.docx
《3DEC经典算例.docx》由会员分享,可在线阅读,更多相关《3DEC经典算例.docx(25页珍藏版)》请在冰豆网上搜索。
3DEC经典算例
3DEC-经典算例,书上内容
仔细学习,融会贯通
P47 ----------------------------------------------
poly brick 0,80 -30,80 0,50
plot block
jset dip 90 dd 180 origin 0,0,0
jset dip 90 dd 180 origin 0,50,0
hide range y -30,0
hide range y 50,80
mark region 1
jset dip 2.5 dd 235 or 30,0,12.5
jset dip 2.5 dd 315 or 35,0,30
jset dip 76 dd 270 spacing 4 number 5 or 38,0,12.5
hide range x 30,80 y 0,50 z 0,50
jset dip 0 dd 0 or 0,0,10
hide range z 0,10
mark region 2 seek
hide range region 0
hide range z 0,10
hide range x 55,80
hide range x 0,30
jset dip 70 dd 200 or 0,35,0
jset dip 60 dd 330 or 50,15,50 seek
hide region 0 2
plot set dip 70 dd 210 mag 2
plot axes block colorby material seek
fix range z 0 10
fix range x 55 80
fix region 0
hide region 0
delete region 2
gravity 0 0 -10 seek
prop mat=1 dens=2000
prop jmat=1 kn=1e9 ks=1e9 f=89
prop jmat=2 kn=1e9 ks=1e9 f=0.0
change dip 90 dd 180 jmat=2
hist zvel (30,30,30) type 1
step 500
plot hist 1 title
new title>rock slope stability
plot hardcopy
save slope.sav
prop jmat=1 f=6.0
cycle 2000
hide reg 0 title
new title> ROCK SLOPE STABILITY -- WEDGE FAILURE
plot dip 70 dd 210 mag 2
plot clear
plot set orien (90,180,0) center (38,25,35) mag 1.5
plot cut add plane origin (0,25,0) name plane normal (0,1,0)
plot add displacement colorby on plane onplane on front off behind off
plot add block fill off plane onp on front off behind off
rest slope.sav
P57 ---------------------------------------------- new
poly brick 0,10 0,10 0,10 pl bl
pl reset
jset origin 5 5 5
prop jmat 1 jkn 1.33e7 jks 1.33e7 fric 20.0
prop mat 1 dens 2000
plot block colorby mat
plot set dip 70 dd 210
cycle 1
list block
list contact
list contact location
gen edge 20
list zone location
P64 ---------------------------------------------- new
poly brick -1 1 -1 1 -1 1
plot block colorby mat
pl reset
plot set dip 70 dd 200
jset dip 65 dd 270 origin .3 0 0
jset dip 40 dd 230 origin 0 -.3 0
jset dip 50 dd 320 origin 0 .3 0
ret
tunnel a (-.3 -1.5 -.3) (-.3 -1.5 .3) (.3 -1.5 .3) (.3 -1.5 -.3) & b (-.3 1.5 -.3) (-.3 1.5 .3) (.3 1.5 .3) (.3 1.5 -.3)
remove range x -.3 .3 y -1.5 1.5 z -.3 .3
plot excavate joint ret
gen edge 1.0
prop mat=1 dens 2000 bulk 1.5e9 g .6e9
prop jmat=1 jkn 1e9 jks 1e9 coh 1e9 ten 1e9 ret
bound stress 0.0 0.0 -1.0e6 0.0 0.0 0.0 range z 0.9 1.1
bound xvel 0.0 range x -1.1 -0.9
bound xvel 0.0 range x 0.9 1.1
bound yvel 0.0 range y -1.1 -0.9
bound yvel 0.0 range y 0.9 1.1
bound zvel 0.0 range z -1.1 -0.9
gravity 0,0,-10.0
insitu stress -0.5e6 -0.5e6 -1.0e6 0.0 0.0 0.0 ret
hist unbal
hist zdis .3 0 .3
hist ty 1
pl hist 1
step 500
save tun0.sav ret
rest tun0.sav
prop jmat 1 fric 6.0 coh 0.0 ten 0.0
reset time hist disp
hist unbal
hist zdisp .3 -0.1 .3
hist ty 2
hide range x -0.4 0.5 y -1 -.5 z .3 .8
plot block disp line color cyan
cycle 5000
save tun1.sav ret
rest tun0.sav
prop jmat 1 fric 6.0 coh 0.0 ten 0.0
struct cable .3 0 .3 .7 0 .7 prop 1 seg 4
struct prop 1 area 5e-4 e 1e9 yield 1e6 kbond 15e8 sbond 1e9
reset time hist disp
hist unbal
hist zdisp .3 -.1 .3
hist ty 2
hide range x -.4 .5 y -1 -.5 z .3 .8
plot block disp line color cyan
plot add hist 2 xaxis label 'Step' yaxis label 'Z-Displacement'
cycle 5000
save tun2.sav ret
P95 ---------------------------------------------- new
poly &
face 0,0,0 1,0,0 1,1,0 0,1,0 &
face 0,0,0 0,0,1 1,0,1 1,0,0 &
face 0,0,0 0,1,0 0,1,1 0,0,1 &
face 1,1,1 1,1,0 1,0,0 1,0,1 &
face 1,1,1 1,0,1 0,0,1 0,1,1 &
face 1,1,1 0,1,1 0,1,0 1,1,0 pl bl
pl reset
pl set dip 60 dd 210
ret
P95 ---------------------------------------------- new
poly brick 0,1 0,1 0,1 pl bl
pl reset
pl set dip 60 dd 210 ret
P97 ---------------------------------------------- new
poly prism &
a(0,0,0)(-.5,.87,0)(-.5,1.87,0)(0,2.74,0) & (1,2.74,0)(1.5,1.87,0)(1.5,.87,0)(1,0,0) &
b(0,0,4)(-.5,.87,4)(-.5,1.87,4)(0,2.74,4) & (1,2.74,4)(1.5,1.87,4)(1.5,.87,4)(1,0,4) pl bl
pl reset
pl set dip 60 dd 200 ret
P98 ----------------------------------------------
New
poly tunnel rad=2 leng=-10,10 ratr=3.0 dip=0 dd=0 nr=2 nt=1 nx=3 pl bl
pl reset
delete -2 2,-10 10 -2 2 ret
join on
P101 ---------------------------------------------- new
poly brick 0,1 0,1 0,1
pl block
pl reset
jset dip 0 dd 0 or 0,0,.5
hide dip 0 dd 0 or 0,0,.5 below
jset dip 90 dd 90 or .5,0,0 seek ret
hide (0.5,1.0) (0,1) (0.5,1.0)
join on
P104 ---------------------------------------------- new
poly brick 0 80 -30 80 0 50
pl block
pl reset
plot set dip 70 dd 210
; shallow-dipping fracture planes (continuous)
jset dip 2.45 dd 235 org 30 0 12.5
jset dip 2.45 dd 315 org 35 0 30
; high angle foliation planes (continuous)
jset dip 76 dd 270 spac 16 num 3 org 30,0,12.5
; intersecting discontinuities (non-continuous)
hide range x 0,80 y 0,50 z 0,10
hide range x 55,80 y 0,50 z 0,50
jset dip 70 dd 200 org 0 35 0
jset dip 60 dd 330 org 50 15 50 seek
hide range x 0,30 y -30,80 z 13,50 ret
P104 ---------------------------------------------- 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.00) (-0.3,-1.5,0.40) (-0.25,-1.5,0.47) & (-0.15,-1.5,0.52) ( 0.0,-1.5,0.55) ( 0.15,-1.5,0.52) & ( 0.25,-1.5,0.47) ( 0.3,-1.5,0.40) ( 0.30,-1.5,0.00) &
b (-0.30, 1.5,0.00) (-0.3, 1.5,0.40) (-0.25, 1.5,0.47) & (-0.15, 1.5,0.52) ( 0.0, 1.5,0.55) ( 0.15, 1.5,0.52) & ( 0.25, 1.5,0.47) ( 0.3, 1.5,0.40) ( 0.30, 1.5,0.00)
remove region 1 ret
P104 ------------------跳过下列例子--------------P124
算例算例算例算例算例算例算例算例算例算例 算例算例算例算例算例算例算例算例算例算例 new
poly brick 0 20 0 25 0 20
jset dip 0.0 or 0 10 0
jset dip 0.0 or 0 15 0
gen edge 5.0
change 0 20 0 10 0 20 mat 1
change 0 20 10 15 0 20 mat 2
change 0 20 15 25 0 20 mat 3
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 kn 1e10 ks 1e10 coh 1e10
gravity 0 -10 0
insitu 0 20 0 10 0 20 stress 0 -4.8e5 0 0 0 0 & ygrad 0 2.2e4 0 0 0 0
insitu 0 20 10 15 0 20 stress 0 -4.6e5 0 0 0 0 & ygrad 0 2.0e4 0 0 0 0
insitu 0 20 15 25 0 20 stress 0 -4.0e5 0 0 0 0 & ygrad 0 1.6e4 0 0 0 0
bound 0,20 -0.1,0.1 0,20 yvel=0.0
bound -0.1,0.1 0,20 0,20 xvel=0.0
bound 19.9,20.1 0,20 0,20 xvel=0.0
bound 0,20 0,20 -0.1,0.1 zvel=0.0
bound 0,20 0,20 19.9,20.1 zvel=0.0
hist unbal
step 500 ret
算例算例算例算例算例算例算例算例算例算例 算例算例算例算例算例算例算例算例算例算例
new
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
3dec汉语翻译3-豆丁网--3dec汉语翻译3-豆丁网--3dec汉语翻译3-豆丁网--3dec汉语翻译3-豆丁网
3dec汉语翻译3-豆丁网--3dec汉语翻译3-豆丁网--3dec汉语翻译3-豆丁网--3dec汉语翻译3-豆丁网
3dec汉语翻译3-豆丁网--3dec汉语翻译3-豆丁网--3dec汉语翻译3-豆丁网--3dec汉语翻译3-豆丁网
3dec汉语翻译3-豆丁网:
P18 3dec汉语翻译3-豆丁网:
P18 3dec汉语翻译3-豆丁网:
P18
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汉语翻译3-豆丁网:
P21-36 3dec汉语翻译3-豆丁网:
P21-36 3dec汉语翻译3-豆丁网:
P21-36 new ;
;example of use of poly cube command ;
config LH
poly br 0 100 0 40 0 100
plot block colorby reg
pl reset
poly cube dip 90 dd 180 num 10 5 10 spac 4 8 4 top 50 0 50 seek
hide reg 1 ret
P124 ---------------------------------------------- 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
plot block colorby mat fill off axes
plot add vel line color red ret
P126 ---------------------------------------------- new
poly brick 0 10 0 10 0 10 pl bl
plot reset
gen quad ndiv 4 4 4
prop mat=1 dens 1000 bulk 8e9 g 5e9
bound zvel=0 range z 0
bound stress -1e5 0 0 0 0 0 range x 0
bound stress -1e5 0 0 0 0 0 range x 10
bound stress 0 -1e5 0 0 0 0 range y 0
bound stress 0 -1e5 0 0 0 0 range y 10
bound stress 0 0 -2e5 0 0 0 range z 10
hist zdisp 0 0 0
step 100
plot axes block fill off vel line color blue