《隧道及地下工程FLAC解析方法》第7章 命令流--自己手输.doc

第七章 命令流按照顺序进行new;1/4 上半部分gen zon radcyl p0 15.7 0 0 p1 24.5 0 0 p2 15.7 1 0 p3 15.7 0 8.8 & size 5 2 10 4 dim 5.8 5.8 5.8 5.8 rat 1 1 1 1.2 group outrockysgen zone cshell p0 15.7 0 0 p1 21.5 0 0 p2 15.7 1 0 p3 15.7 0 5.8 & size 1 2 10 4 dim 5.4 5.4 5.4 5.4 rat 1 1 1 1 group liningys fill group inrockysgen zon reflect dip 90 dd 270 orig 15.7 0 0;1/4 下半部分gen zon radcyl p0 15.7 0 0 p1 15.7 0 -6.9 p2 15.7 1 0 p3 24.5 0 0 & size 5 2 10 4 dim 3.9 5.8 3.9 5.8 rat 1 1 1 1.2 group outrockyxgen zone cshell p0 15.7 0 0 p1 15.7 0 -3.9 p2 15.7 1 0 p3 21.5 0 0 & size 1 2 10 4 dim 3.5 5.4 3.5 5.4 rat 1 1 1 1 group liningyx fill group inrockyxgen zon reflect dip 90 dd 270 orig 15.7 0 0 range group liningyx any & group inrockyx any group outrockyx any;右隧道中地层网格gen zon brick p0 0 0 -6.9 p1 6.9 0 -6.9 p2 0 1 -6.9 p3 0 0 0 size 6 2 5 group outrock1gen zon brick p0 24.5 0 -6.9 p1 65 0 -6.9 p2 24.5 1 -6.9 p3 24.5 0 0 size 20 2 5 group outrock1gen zon brick p0 0 0 0 p1 6.9 0 0 p2 0 1 0 p3 0 0 8.8 size 6 2 5 group outrock2gen zon brick p0 24.5 0 0 p1 65 0 0 p2 24.5 1 0 p3 24.5 0 8.8 size 20 2 5 group outrock2;右隧道上侧土体网格gen zon brick p0 0 0 8.8 p1 6.9 0 8.8 p2 0 1 8.8 p3 0 0 36 size 6 2 10 group outrock3gen zon brick p0 24.5 0 8.8 p1 65 0 8.8 p2 24.5 1 8.8 p3 24.5 0 36 size 20 2 10 group outrock3gen zon brick p0 6.9 0 8.8 p1 24.5 0 8.8 p2 6.9 1 8.8 p3 6.9 0 36 size 10 2 10 group outrock3;右隧道下侧土体网格gen zon brick p0 0 0 -24 p1 6.9 0 -24 p2 0 1 -24 p3 0 0 -6.9 size 6 2 8 group outrock4gen zon brick p0 24.5 0 -24 p1 65 0 -24 p2 24.5 1 -24 p3 24.5 0 -6.9 size 20 2 8 groupoutrock4gen zon brick p0 6.9 0 -24 p1 24.5 0 -24 p2 6.9 1 -24 p3 6.9 0 -6.9 size 10 2 8 group outrock4gen zon reflect dip 90 dd 270 orig 0 0 0plot add axes redplot block group;设置重力加速度set gravity 0 0 -10;给定边界条件fix z range z -24.01,-23.99fix x range x -65.01,-64.99fix x range x 64.99,65.01fix y range y -0.01 0.01fix y range y 0.99,1.01;0-求解自重应力场model mohr ini density 2200 ;围岩的密度prop bulk=1.923e9, shear=4.167e9 fric=35, coh=0.5e6 tension=5e5;体积、剪切、摩擦角、凝聚力、抗拉强度step 6000save Gsol.sav plot cont zdispplot cont szz;1-左隧道上半断面开挖rest Gsol.savini xdis 0.0 ydis 0.0 zdis 0.0ini density 2400 range group outrockzs any prop bulk=6.0e9, shear=10.0e9 fric=45, coh=1.2e6 tension=2e6 &range group outrockzs anymodel null range group inrockzs any group liningzs any step 2000save stepp1.sav;2-初期支护rest stepp1.savsel shell id=1 range x -21.51 -9.89 cyl end1 -15.7 0 0 end2 -15.7 1 0 rad 5.8sel shell id=1 prop iso=(21.0e9,0.2) thick=0.20step 2000save stepp2.sav;3-左隧道下半断面开挖rest stepp2.savmodel null range group inrockzx any group liningzx any step 2000save stepp3.savplot add axes redplot block groupplot cont zdispplot cont szz;4-左隧道二次衬砌rest stepp3.savmodel elas range group liningzsmodel elas range group liningzxini density 2600 range group liningzsini density 2600 range group liningzxprop bulk=12.5e9, shear=16.667e9 range group liningzsprop bulk=12.5e9, shear=16.667e9 range group liningzxstep 2000save stepp4.sav plot add axes redplot block groupplot cont zdispplot cont szz;5-右隧道上半断面开挖rest stepp4.savini density 2400 range group outrockys any prop bulk=6.0e9, shear=10.0e9 fric=45, coh=1.2e6 tension=2e6 &range group outrockys anymodel null range group inrockys any group liningys any step 2000save stepp5.savplot cont zdisp;6-初期支护rest stepp5.savsel shell id=1 range x 9.89 21.51 cyl end1 15.7 0 0 end2 15.7 1 0 rad 5.8sel shell id=1 prop iso=(21.0e9,0.2) thick=0.20step 2000plot cont zdisp save stepp6.savplot add axes redplot block groupplot cont zdispplot cont szz;7-左隧道下半断面开挖rest stepp6.savmodel null range group inrockyx any group liningyx any step 2000save stepp7.savplot cont zdisp;8-左隧道二次衬砌rest stepp7.savmodel elas range group liningysmodel elas range group liningyxini density 2600 range group liningysini density 2600 range group li