【正文】 30 plot cont zdisp outl on417不規(guī)則自由面應(yīng)力初始化gen zone brick size 15 15 10 p0 0,0,0 edge=model elasticprop shear 3e8 bulk 5e8def mountaingp = gp_headloop while gp nullzz = sqrt(gp_xpos(gp)?2 + gp_ypos(gp)?2)dz = * sin( * zz + ) 。網(wǎng)格三維顯示418非均網(wǎng)格內(nèi)部壓實(shí)new gen zone brick size 8 8 10 ratio 1 1model mohrini dens 2000prop bulk 2e8 shear 1e8prop fric 30fix x range x .1fix x range x fix y range y .1fix y range y fix z range z .1。pause。step 網(wǎng)格三維顯示plot set rotation 30 0 30 plot cont zdisp outl on420應(yīng)力與孔隙壓力的初始化new config fluidgen zone brick size 8 5 10model elasmodel fl_isoini dens 2000prop bulk 1e9 shear 5e8prop poros perm 1e10ini fmod 2e9ini fdensity 1e3ini sat 0ini sat 1 range z set grav 0 0 10fix x range x .1fix x range x fix y range y .1fix y range y fix z range z .1ini pp grad 0,0, range z 5.ini szz 20e4 grad 0,0,20e3ini szz add grad 0,0,.25e4 range z 5,6ini szz add grad 0,0,.5e4 range z 0,5ini szz add range z 0,5solve。size 4 10 6 4 dim 4 4 4 4 rat 1 1 1 1 fillgen zon reflect dip 90 dd 90 orig 15 0 0gen zon reflect dip 0 ori 0 0 0。p5 7 50 0 p6 7 0 8 p7 7 50 8 p9 point 1 p11 point 2 amp。p8 23 0 0 p9 7 0 8 p10 23 50 0 p11 7 50 8 amp。end2 , rad save rest 。 symmetry cond.sel node fix x yr zr range x 。 size 1 2 10 4 dim rat 1 1 1 1 group concretlinersy fill group insiderocksygen zon radcyl p0 0 0 0 p1 0 0 9 p2 0 1 0 p3 0 0 amp。size 5 2 12 group outsiderock2gen zon brick p0 0 0 p1 50 0 0 p2 1 0 p3 0 amp。size 20 2 12 group outsiderock5gen zon reflect dip 90 dd 90 orig 0 0 0plot set rotation 30 0 30 。設(shè)置左邊界fix x range x , 。莫爾～庫侖模型ini density 2000 。plot cont zdisp 。繪制豎向應(yīng)力場，如圖617所示 。5右上半斷面開挖rest plot block groupplot add axes redmodel null range group insiderocksy any group concretlinersy anystep 2000save plot cont zdisp rest sel shell id=1 range x z cyl end1 0 0 0 end2 0 1 0 rad sel shell id=1 prop iso=(,) thick=step 2000plot cont zdisp save 。 dim rat 1 1 1 1 group concretlinersy fill group insiderocksy。1/2繪制上下地層網(wǎng)格gen zon brick p0 0 0 p1 0 p2 0 1 p3 0 0 size 5 2 8 group outsiderock1
gen zon brick p0 0 0 p1 0 p2 0 1 p3 0 0 9 size 5 2 12 group outsiderock2。設(shè)置重力加速度set gravity 0 0 10。體積、剪切、摩擦角、凝聚力、抗拉強(qiáng)度step 6000save 。 group insiderocksy any group concretlinersy anystep 2000save plot cont zdisp 。4二次襯砌rest model elas range group concretlinerxymodel elas range group concretlinerxzmodel elas range group concretlinersymodel elas range group concretlinerszini density 2600 range group concretlinerxyini density 2600 range group concretlinerxzini density 2600 range group concretlinersyini density 2600 range group concretlinerszprop bulk=, shear= range group concretlinerxyprop bulk=, shear= range group concretlinerxzprop bulk=, shear= range group concretlinersyprop bulk=, shear= range group concretlinerszstep 2000save III級(jí)圍巖施工過程模擬new。1/4 下半部分gen zon radcyl p0 0 0 0 p1 0 0 9 p2 0 1 0 p3 0 0 amp。1/2右側(cè)中間部分土體網(wǎng)格gen zon brick p0 0 0 p1 50 0 0 p2 1 0 p3 0 size 20 2 5 group outsiderock3gen zon brick p0 0 9 p1 50 0 9 p2 1 9 p3 0 0 size 20 2 5 group outsiderock3。給定邊界條件fix z range z ,fix x range x ,fix x range x ,fix y range y fix y range y ,。1全斷面開挖rest ini xdis ydis zdis model null range group insiderocksz any group concretlinersz any amp。3二次襯砌rest model elas range group concretlinerxymodel elas range group concretlinerxzmodel elas range group concretlinersymodel elas range group concretlinerszini density 2600 range group concretlinerxyini density 2600 range group concretlinerxzini density 2600 range group concretlinersyini density 2600 range group concretlinerszprop bulk=, shear= range group concretlinerxyprop bulk=, shear= range group concretlinerxzprop bulk=, shear= range group concretlinersyprop bulk=, shear= range group concretlinerszstep 2000save plot block group第七章 分離式公路隧道施工過程分析72 IV級(jí)圍巖施工過程模擬new。1/4 下半部分gen zon radcyl p0 0 0 p1 0 p2 1 0 p3 0 0 amp。 size 5 2 10 4 dim rat 1 1 1 group outrockysgen zone cshell p0 0 0 p1 0 0 p2 1 0 p3 0 amp。 group insiderockxy any group concretlinerxy anystep 2000save 。圍巖的密度prop bulk=, shear= fric=45, coh= tension=2e6。對稱得出完整模型gen zon reflect dip 90 dd 90 orig 0 0 0。 dim rat 1 1 1 1 group concretlinerxy fill group insiderockxy。size 5 2 10 4 dim rat 1 1 1 group outsiderocksygen zone cshell p0 0 0 0 p1 0 0 p2 0 1 0 p3 0 0 size 1 2 10 4 amp。3下半斷面開挖rest model null range group insiderockxz any group concretlinerxz any amp。new。0求解自重應(yīng)力場model mohr ini density 2200 。1/2右上側(cè)和下側(cè)土體網(wǎng)格gen zon brick p0 0 p1 50 0 p2 1 p3 0 size 20 2 8 group outsiderock4
gen zon brick p0 0 p1 50 0 p2 1 p3 0 9 size 20 2 12 group outsiderock5
。 size 5 2 10 4 dim 4 4 rat 1 1 1 group outsiderockxygen zone cshell p0 0 0 0 p1 0 0 4 p2 0 1 0 p3 0 0 size 1 2 10 4amp。1/4 上半部分gen zon radcyl p0 0 0 0 p1 0 0 p2 0 1 0 p3 0 0 size 5 2 10 4amp。 tension=5e5, range group outsiderocksy any group outsiderocksz anymodel null range group insiderocksz any group concretlinersz anystep 2000save plot cont zdisprest sel shell id=1 range x z cyl end1 0 0 0 end2 0 1 0 rad sel shell id=1 range x z sel shell id=1 prop iso=(,) thick=step 2000plot cont zdisp save 。繪制豎向位移場，如圖615所示plot cont szz 。求解6000次Solve 。設(shè)置前邊界fix y range y , 。設(shè)置重力加速度為z方向10fix z range z , 。size 20 2 5 group outsiderock3gen zon brick p0 0 p1 50 0 p2 1 p3 0 amp。 size 1 2 10 4 dim rat 1 1 1 1 group concretlinerxy fill group insiderockxy