3DEC经典算例.docx

1、3DEC经典算例3DEC-经典算例，书上内容仔细学习，融会贯通P47-polybrick0,80-30,800,50plotblockjsetdip90dd180origin0,0,0jsetdip90dd180origin0,50,0hiderangey-30,0hiderangey50,80markregion1jsetdip2.5dd235or30,0,12.5jsetdip2.5dd315or35,0,30jsetdip76dd270spacing4number5or38,0,12.5hiderangex30,80y0,50z0,50jsetdip0dd0or0,0,10hideran

2、gez0,10markregion2seekhiderangeregion0hiderangez0,10hiderangex55,80hiderangex0,30jsetdip70dd200or0,35,0jsetdip60dd330or50,15,50seekhideregion02plotsetdip70dd210mag2plotaxesblockcolorbymaterialseekfixrangez010fixrangex5580fixregion0hideregion0deleteregion2gravity00-10seekpropmat=1dens=2000propjmat=1k

3、n=1e9ks=1e9f=89propjmat=2kn=1e9ks=1e9f=0.0changedip90dd180jmat=2histzvel(30,30,30)type1step500plothist1titlenewtitlerockslopestabilityplothardcopysaveslope.savpropjmat=1f=6.0cycle2000hidereg0titlenewtitleROCKSLOPESTABILITY-WEDGEFAILUREplotdip70dd210mag2plotclearplotsetorien(90,180,0)center(38,25,35)

4、mag1.5plotcutaddplaneorigin(0,25,0)nameplanenormal(0,1,0)plotadddisplacementcolorbyonplaneonplaneonfrontoffbehindoffplotaddblockfilloffplaneonponfrontoffbehindoffrestslope.savP57-newpolybrick0,100,100,10plblplresetjsetorigin555propjmat1jkn1.33e7jks1.33e7fric20.0propmat1dens2000plotblockcolorbymatplo

5、tsetdip70dd210cycle1listblocklistcontactlistcontactlocationgenedge20listzonelocationP64-newpolybrick-11-11-11plotblockcolorbymatplresetplotsetdip70dd200jsetdip65dd270origin.300jsetdip40dd230origin0-.30jsetdip50dd320origin0.30rettunnela(-.3-1.5-.3)(-.3-1.5.3)(.3-1.5.3)(.3-1.5-.3)&b(-.31.5-.3)(-.31.5.

6、3)(.31.5.3)(.31.5-.3)removerangex-.3.3y-1.51.5z-.3.3plotexcavatejointretgenedge1.0propmat=1dens2000bulk1.5e9g.6e9propjmat=1jkn1e9jks1e9coh1e9ten1e9retboundstress0.00.0-1.0e60.00.00.0rangez0.91.1boundxvel0.0rangex-1.1-0.9boundxvel0.0rangex0.91.1boundyvel0.0rangey-1.1-0.9boundyvel0.0rangey0.91.1boundz

7、vel0.0rangez-1.1-0.9gravity0,0,-10.0insitustress-0.5e6-0.5e6-1.0e60.00.00.0rethistunbalhistzdis.30.3histty1plhist1step500savetun0.savretresttun0.savpropjmat1fric6.0coh0.0ten0.0resettimehistdisphistunbalhistzdisp.3-0.1.3histty2hiderangex-0.40.5y-1-.5z.3.8plotblockdisplinecolorcyancycle5000savetun1.sa

8、vretresttun0.savpropjmat1fric6.0coh0.0ten0.0structcable.30.3.70.7prop1seg4structprop1area5e-4e1e9yield1e6kbond15e8sbond1e9resettimehistdisphistunbalhistzdisp.3-.1.3histty2hiderangex-.4.5y-1-.5z.3.8plotblockdisplinecolorcyanplotaddhist2xaxislabelStepyaxislabelZ-Displacementcycle5000savetun2.savretP95

9、-newpoly&face0,0,01,0,01,1,00,1,0&face0,0,00,0,11,0,11,0,0&face0,0,00,1,00,1,10,0,1&face1,1,11,1,01,0,01,0,1&face1,1,11,0,10,0,10,1,1&face1,1,10,1,10,1,01,1,0plblplresetplsetdip60dd210retP95-newpolybrick0,10,10,1plblplresetplsetdip60dd210retP97-newpolyprism&a(0,0,0)(-.5,.87,0)(-.5,1.87,0)(0,2.74,0)&

10、(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)plblplresetplsetdip60dd200retP98-Newpolytunnelrad=2leng=-10,10ratr=3.0dip=0dd=0nr=2nt=1nx=3plblplresetdelete-22,-1010-22retjoinonP101-newpolybrick0,10,10,1plblockplresetjsetdip0

11、dd0or0,0,.5hidedip0dd0or0,0,.5belowjsetdip90dd90or.5,0,0seekrethide(0.5,1.0)(0,1)(0.5,1.0)joinonP104-newpolybrick080-3080050plblockplresetplotsetdip70dd210;shallow-dippingfractureplanes(continuous)jsetdip2.45dd235org30012.5jsetdip2.45dd315org35030;highanglefoliationplanes(continuous)jsetdip76dd270sp

12、ac16num3org30,0,12.5;intersectingdiscontinuities(non-continuous)hiderangex0,80y0,50z0,10hiderangex55,80y0,50z0,50jsetdip70dd200org0350jsetdip60dd330org501550seekhiderangex0,30y-30,80z13,50retP104-newpolybrick-1.5,1.5-1.5,1.5-1.5,1.5plblcolorbymaterialplresetplsetdip70dd200tunnelradialregion1&a(-0.30

13、,-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)removeregion1retP10

14、4-跳过下列例子-P124算例算例算例算例算例算例算例算例算例算例算例算例算例算例算例算例算例算例算例算例newpolybrick020025020jsetdip0.0or0100jsetdip0.0or0150genedge5.0change020010020mat1change0201015020mat2change0201525020mat3propmat=1dens1600bulk5e9g3e9propmat=2dens2000bulk5e9g3e9propmat=3dens2200bulk5e9g3e9changejmat1propjmat1kn1e10ks1e10coh1e10gr

15、avity0-100insitu020010020stress0-4.8e50000&ygrad02.2e40000insitu0201015020stress0-4.6e50000&ygrad02.0e40000insitu0201525020stress0-4.0e50000&ygrad01.6e40000bound0,20-0.1,0.10,20yvel=0.0bound-0.1,0.10,200,20xvel=0.0bound19.9,20.10,200,20xvel=0.0bound0,200,20-0.1,0.1zvel=0.0bound0,200,2019.9,20.1zvel=

16、0.0histunbalstep500ret算例算例算例算例算例算例算例算例算例算例算例算例算例算例算例算例算例算例算例算例newpolybrick0,200,200,20plblplresetplsetdip70dd210genedge4propmat=1dens2500bulk5e9g3e9phi35changecons2gravity00-10boundstress-2.75e6,-2.75e6,-5.5e60,0,0zgrad1.25e4,1.25e4,2.5e40,0,0rangex0boundstress-2.75e6,-2.75e6,-5.5e60,0,0zgrad1.25e4,

17、1.25e4,2.5e40,0,0rangex20boundstress-2.75e6,-2.75e6,-5.5e60,0,0zgrad1.25e4,1.25e4,2.5e40,0,0rangey0boundstress-2.75e6,-2.75e6,-5.5e60,0,0zgrad1.25e4,1.25e4,2.5e40,0,0rangey20boundstress0,0,-5e60,0,0rangez20boundzvel=0rangez0insitustress-2.75e6,-2.75e6,-5.5e60,0,0zgrad1.25e4,1.25e4,2.5e40,0,0step500p

18、lotblockfilloffstresscolorbymaximumreverseret3dec汉语翻译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-豆丁网：P183dec汉语翻译3-豆丁网：P183dec汉语翻译3-豆丁网：P18newpolybrick-1.51.5-1.51.5-1.51.5plblco

19、lorbymaterialplresetplsetdip70dd200tunnelradialregion1&a(-0.30-1.50)(-0.3-1.50.4)(-0.25-1.50.47)&(-0.15-1.50.52)(0-1.50.55)(0.15-1.50.52)&(0.25-1.50.47)(0.3-1.50.4)(0.3-1.50)&b(-0.301.50)(-0.31.50.4)(-0.251.50.47)&(-0.151.50.52)(01.50.55)(0.151.50.52)&(0.251.50.47)(0.31.50.4)(0.31.50)removeregion1re

20、t3dec汉语翻译3-豆丁网：P21-363dec汉语翻译3-豆丁网：P21-363dec汉语翻译3-豆丁网：P21-36new;exampleofuseofpolycubecommand;configLHpolybr01000400100plotblockcolorbyregplresetpolycubedip90dd180num10510spac484top50050seekhidereg1retP124-newpolybrick010010010plotblockaxesplresetplotsetdip70dd150jsetdip0dd180origin005genquadndiv44

21、4changecons1propmat=1dens1000bulk8e9g5e9propjmat=1jkn1e10jks1e10gravity00-10boundstress00-1e5000rangex010y010z-.1.1boundxvel=0rangex-.1.1boundxvel=0rangex9.910.0boundyvel=0rangey-.1.1boundyvel=0rangey9.910.1insitustress00-1e5000zgrad001e4000histzdisp552.5step300excavaterangex010y010z510step100plotbl

22、ockcolorbymatfilloffaxesplotaddvellinecolorredretP126-newpolybrick010010010plblplotresetgenquadndiv444propmat=1dens1000bulk8e9g5e9boundzvel=0rangez0boundstress-1e500000rangex0boundstress-1e500000rangex10boundstress0-1e50000rangey0boundstress0-1e50000rangey10boundstress00-2e5000rangez10histzdisp000step100plotaxesblockfilloffvellinecolorblue