3DEC经典算例.doc

1、 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

2、.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

3、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 titlenew titlerock

4、 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 dis

5、placement 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.0prop mat 1 dens 2000 plot block colorby mat plot set dip

6、 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 rettunnel a (

7、-.3 -1.5 -.3) (-.3 -1.5 .3) (.3 -1.5 .3) (.3 -1.5 -.3) 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-contin

8、uous) 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 resetpl set dip 70 dd 200 tunnel radial region 1

9、 ;example of use of poly cube command ; config LH poly br 0 100 0 40 0 100 plot block colorby reg pl resetpoly 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

10、 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.1insitu stre

11、ss 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

12、=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 10hist zdisp 0 0 0 step 100 plot axes block fill off vel line color blue ret P132 - new poly bric

13、k 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 bound stress 0 0 -1e

14、7 0 0 0 range z 10 insitu stress -5e6 -5e6 -1e7 0 0 0 step 1 plot block fill off stress colorby maximum ret P133 -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 gen edge .3 g

15、en quad ndiv 1 5 10 seek gen quad ndiv 1 2 10 ;gen edge 1 change cons 2prop mat=1 dens 2000 bulk 2e8 g 1e8 phi 30 prop jmat 1 jkn 1e8 jks 1e8 coh 1e10 ten 1e10 gravity 0 0 -10 bound xvel=0 range x -0.1 0.1 bound xvel=0 range x 9.9 10.1 bound yvel=0 range y -0.1 0.1 bound yvel=0 range y 1.9 2.1 bound

16、 zvel=0 range z -0.1 0.1 hist unbal sav ex3_14.sav step 500 plot bcont szz max 0 min -200000 int 20000 pause key ;optional method 1 restore ex3_14.sav insitu stress -1.5e5 -1.5e5 -2.0e5 0 0 0 optional method 2 restore ex3_14.sav prop mat=1 bcoh 1e10 bten 1e10 step 750 prop mat=1 bcoh 0 bten 0step 25

17、0 ret P139 - new poly brick 0 5 0 5 0 5 plot block colorby model plot reset gen edge 1.0 prop mat=1 dens 2000 bulk 3e8 g 2e8 gravity 0 0 -10 bound xvel=0 zvel=0 range z -0.1 0.1 hist unbal step 250 pause key change cons 2 prop mat=1 dens 2000 bulk 3e8 g 2e8 phi 34 ret P140 - new poly brick -10 10 -1

18、0 10 -20 0plot block plot reset jset dip 60 dd 90 origin 0 0 -10 jset dip 20 dd 90 origin 0 0 -8 jset dip 20 dd 90 origin 0 0 -12 join on range z -20 -14 join on range z -7 0 gen edge 2 prop mat=1 dens 2000 bulk 8e9 g 5e9 prop jmat=1 jkn 5e11 jks 2.5e11 fric 30 insitu stress -2.5e6 -2.5e6 -1e7 0 0 0

19、 bound stress -2.5e6 0 0 0 0 0 range x -10 bound stress -2.5e6 0 0 0 0 0 range x 10 bound stress 0 -2.5e6 0 0 0 0 range y -10 bound stress 0 -2.5e6 0 0 0 0 range y 10 bound stress 0 0 -1e7 0 0 0 range z 0 bound zvel 0 range z -20 hist unbal hist zdis 0 0 -10 step 1000 plot clear plot cut add plane o

20、rigin 0 0 -10 name “Plane 1“ normal 0 1 0plot add jointvector shear shearoffset 0.1 plane on point size 0 - tunz: FISH function to define tunnel geometry parameters -; . tunnel along Y axis, from YYA to YYB ; . semi-circular roof, centered at (TXC,TZC) ; def tunz ; yya = -10.0 yyb = 10.0 ; ; - outer

21、 surface - txb1 = -4.0 tzb1 = -4.0 txb2 = 4.0 tzb2 = -4.0 ; txc = 0.0 tzc = 0.0 tr = 4.0 tx1 = txc + tr * cos(180*degrad) tz1 = tzc + tr * sin(180*degrad) tx2 = txc + tr * cos(135*degrad) tz2 = tzc + tr * sin(135*degrad) tx3 = txc + tr * cos(90*degrad) tz3 = tzc + tr * sin(90*degrad) tx4 = txc + tr

22、* cos(45*degrad)tz4 = tzc + tr * sin(45*degrad) tx5 = txc + tr * cos(0*degrad) tz5 = tzc + tr * sin(0*degrad) ; ; - inner surface - ; thickness th th = 0.5 txb1i = -4.0 + th tzb1i = -4.0 + th txb2i = 4.0 - th tzb2i = -4.0 + th ; txc = 0.0 tzc = 0.0 tri = tr - th tx1i = txc + tri * cos(180*degrad) tz

23、1i = tzc + tri * sin(180*degrad) tx2i =txc + tri * cos(135*degrad) tz2i = tzc + tri * sin(135*degrad) tx3i = txc +tri * cos(90*degrad) tz3i = tzc + tri * sin(90*degrad) tx4i = txc + tri * cos(45*degrad) tz4i = tzc + tri * sin(45*degrad) tx5i = txc + tri * cos(0*degrad) tz5i = tzc + tri * sin(0*degra

24、d) ; end ; ; -; (execute function) tunz ; ; create outer surface tunnel radial ; create inner surfacetunnel radial ; -NOTE: region inside inner surface is REG 7 ; region between surface (to be liner) is REG 5 ; save tun_a.sav ; ; - joints - 3 joints to form a wedge in the roof ; jset dd 270 dip 70 o

25、rigin 0 0,5.7 id 10 jset dd 40 dip 65 origin 0 0,5.7 id 10 jset dd 130 dip 60 origin 0 0,5.7 id 10 ; save tun_b.sav ; ; - mesh generation - ; rock blocks hide reg 5 7 gen edge 5 ; ; liner find reg 5 gen edge 2 ; find reg 7 gen edge 2 ; save tun_z.sav plot block colorby matplot set dip 70 dd 210 ret

26、P150 - rest tun_z.sav ; ; - properties - ; ; - MAT 1 : rock - ; density = 2700 kg/m3 = 0.0027e6 kg/m3 ; E=50 GPa, Poissons ratio=0.2 prop mat 1 dens 0.0027 k 27778 g 20833 ; ; - MAT=5 : concrete liner - ; density = 2400 kg/m3 = 0.0024e6 kg/m3 ; E=30 GPa, Poissons ratio=0.2 prop mat 5 dens 0.0025 k 1

27、6667 g 12500 ; ; - JMAT=1 : rock joints - prop mat 1 jkn 10000 jks 2000 fric 25 ; ; - JMAT=5 : concrete-concrete joints (elastic) -prop mat 5 jkn 30000 jks 12000 coh 1e6 tens 1e6 ; ; - JMAT=6 : concrete-rock interface -prop mat 6 jkn 10000 jks 2000 fric 0.001 ; ; - assign material numbers -; initial

28、ly all materials are rock change mat 1 change jmat 1 ; ; - insitu stress state - ; assume tunnel at 200 m depth ; vertical stress: syy=(0.0027*g)*(y-200) ; at z=0: szz=-5.4 ; z-gradient of szz: 0.027 ; (positive: less compression going up) ; horizontal sxx=szz=0.5*syy ; insitu stress -2.7 -2.7 -5.4

29、0 0 0 gravity grav 0 0 -10 ; ; - boundary conditions for insitu stress state -; top of model (z=10): szz=-0.027*190=-5.13 bound stress 0 0 -5.13 0 0 0 range z 10.0 ; bottom bound zvel 0 range z -10.0 ; sides bound xvel 0 range x -10.0 bound xvel 0 range x 10.0 bound yvel 0 range y -10.0 bound yvel 0

30、 range y 10.0 ; ; - histories to monitor convergence - hist nc=1 unbal ; top of model hist xdis 0 0 10 ydis 0 0 10 zdis 0 0 10 ; save tun_c0.sav plot hist 2 3 4 yaxis label Displacement cycle 1000 save tun_c.sav ret P152 - rest tun_c.sav pl block pl reset ; delete interior blocks remove reg 5 7 paus

31、e key ; ; history point at tunnel roof reset disp time hist hide seek range xr .44 .46 hist zdis 1.8 0 3.72 hist label 1 Vertical Displacement seek ; pl hist 1 yaxis label Displacement cycle 5000 ; save tun_x.sav pause key hide dip 60 dd 130 or 0 0 5.7 above pl block velocity line color cyan plot se

32、t dip 80 dd 180 mag 3 ret P155 - rest tun_c.sav ; ; delete interior blocks remove region 7 ; ; delete liner blocks remove reg 5 ; ; - install axial elements - struct axial -8 -5 -2 -3.9 -5 -2 prop 7 struct axial -8 0 -2 -3.9 0 -2 prop 7 struct axial -8 5 -2 -3.9 5 -2 prop 7 struct axial -6.8 -5 6.8

33、-2.8 -5 2.8 prop 7 struct axial -6.8 0 6.8 -2.8 0 2.8 prop 7 struct axial -6.8 5 6.8 -2.8 5 2.8 prop 7 ; struct axial 8 -5 -2 3.9 -5 -2 prop 7 struct axial 8 0 -2 3.9 0 -2 prop 7 struct axial 8 5 -2 3.9 5 -2 prop 7 struct axial 6.8 -5 6.8 2.8 -5 2.8 prop 7 struct axial 6.8 0 6.8 2.8 0 2.8 prop 7 str

34、uct axial 6.8 5 6.8 2.8 5 2.8 prop 7 ; struct axial 0 -5 4 0 -5 8 prop 7 struct axial 0 0 4 0 0 8 prop 7struct axial 0 5 4 0 5 8 prop 7 ; struct prop 7 rkax 250 rlen 0.10 rult 0.55 ; reset disp time hist ; history point at tunnel roof hide seek range xr .44 .46 hist zdis 1.8 0 3.72 hist label 1 Vert

35、ical Displacement seek ; pl hist 1 yaxis label Displacement cy 2000 ; save tun_lr.sav ret P156 - rest tun_c.sav ; ; delete interior blocks plot block remove region 7 ; ; delete liner blocks remove reg 5 ; ; - install cable elements - struct cable -8 -5 -2 -4.05 -5 -2 prop 8 seg 4 struct cable -8 0 -

36、2 -4.05 0 -2 prop 8 seg 4 struct cable -8 5 -2 -4.05 5 -2 prop 8 seg 4 struct cable -6.8 -5 6.8 -2.85 -5 2.85 prop 8 seg 4 struct cable -6.8 0 6.8 -2.85 0 2.85 prop 8 seg 4 struct cable -6.8 5 6.8 -2.85 5 2.85 prop 8 seg 4 ; struct cable 8 -5 -2 4.05 -5 -2 prop 8 seg 4 struct cable 8 0 -2 4.05 0 -2

37、prop 8 seg 4 struct cable 8 5 -2 4.05 5 -2 prop 8 seg 4 struct cable 6.8 -5 6.8 2.85 -5 2.85 prop 8 seg 4 struct cable 6.8 0 6.8 2.85 0 2.85 prop 8 seg 4 struct cable 6.8 5 6.8 2.85 5 2.85 prop 8 seg 4 ; struct cable 0 -5 4.1 0 -5 8 prop 8 seg 4 struct cable 0 0 4.1 0 0 8 prop 8 seg 4 struct cable 0

38、 5 4.10 5 8 prop 8 seg 4 ; ; start with high SBOND struct prop 8 area 5e-4 e 100000 yield 0.55 kbond 15e4 sbond 1e6 ; reset disp time hist ; history point at tunnel roof hide seek range xr .44 .46 hist zdis 1.8 0 3.72 hist label 1 Vertical Displacement seek ; cycle 500 ; ; set real SBOND struct prop

39、 8 sbond 0.8 ; pl hist 1 yaxis label Displacement cy 1500 ; save tun_cab.sav pause key pl exc colorby material cable pl set dip 80 dd 190 ret P160 - rest tun_c.sav plot block plot set dip 80 dd 200 ; delete interior blocks delete region 7 ; ; excavate liner blocks (not deleted)excavate reg 5 ; ; his

40、tory point at tunnel roof hide seek range xr .44 .46 hist xdis 1.8 0 3.72 ydis 1.8 0 3.72 zdis 1.8 0 3.72 hist label 7 Vertical Displacement seek ; ; simulate the removal of approximately 50% of insitu stress ; applying at liner-rock interface a stress state ; szz=-2.7 sxx=syy=-1.35 ; bound str -1.3

41、5 -1.35 -2.7 0 0 0 range x -4.1,-3.9 y -11 11 z -4.1 0.1 bound str -1.35 -1.35 -2.7 0 0 0 range x 3.9,4.1 y -11 11 z -4.1 0.1 bound str -1.35 -1.35 -2.7 0 0 0 range x -4.1,4.1 y -11 11 z -4.1 -3.9 ; note: need to include all faces on tunnel surface ; (inner radius must be a bit smaller than 4.0) bou

42、nd str -1.35 -1.35 -2.7 0 0 0 ; must fix again end-surfaces that were freed by BOU STRESS bou yvel 0 range yr -10.1,-9.9 bou yvel 0 range yr 9.9,10.1 ; ; check that sum of applied forces on tunnel surface is zero list boun forcerange x -5,5 y -11,11 z -5,5 ; check that sum of applied x and z forces

43、on tunnel surface is near zero pause ; cycle 2000 ; save tun_l1.sav ; ; - insert liner - ; remove loads from tunnel surface ; bound xfree yfree zfree range x -4.1,-3.9 y -11,11 z -4.1,0.1 bound xfree yfree zfree range x 3.9,4.1 y -11,11 z -4.1,0.1 bound xfree yfree zfree range x -4.1,4.1 y -11,11 z

44、-4.1,-3.9 bound xfree yfree zfree range zr -0.1,4.1 ; must fix again end-surfaces bou yvel 0 range yr -10.1 -9.9 bou yvel 0 range yr 9.9 10.1 ; ; insert liner fill mat 5 jmat 5 range reg 5 ; ; join liner blocks join on range reg 5 ; ; assign rock-liner interface material change jmat 6 range rint 0 5 ; pl hist 7 yaxis label Displacement cy 2000 ;