newtitle基坑开挖模拟—未穿通;---建立壳模型gen zone brick p0 -22.5 15 -0.6 p1 -20.5 13 -0.6 p2 0 17.5 -0.6 p3 -22.5 15 0 &p4 0 15.5 -0.6 p5 0 17.5 0 p6 -20.5 13 0 p7 0 15.5 0gen zone brick p0 -22.5 -1 -0.6 p1 -20.5 0 -0.6 p2 -22.5 15 -0.6 p3 -22.5 -1 0 &p4 -20.5 13 -0.6 p5 -22.5 15 0 p6 -20.5 0 0 p7 -20.5 13.0 0gen zone brick p0 -12.5 -15.5 -0.6 p1 -10.5 -13.5 -0.6 p2 -22.5 -1 -0.6 &p3 -12.5 -15.5 0 p4 -20.5 0 -0.6 p5 -22.5 -1 0 p6 -10.5 -13.5 0 p7 -20.5 0 0 gen zone brick p0 -12.5 -15.5 -0.6 p1 0 -15.5 -0.6 p2 -10.5 -13.5 -0.6 &p3 -12.5 -15.5 0 p4 0 -13.5 -0.6 p5 -10.5 -13.5 0 p6 0 -15.5 0 p7 0 -13.5 0 ;---建立对称单元gen zone reflect dd 90 dip 90 origin 0 0 0group shell;---冠梁参数model elassel shell id=1 range group shellsel shell id=1 elemtype=dkt crossdiag range group shell z -0.6 0sel shell prope density 1600 iso 2e11 0.0 thick 0.6 range group shell;---建立基坑模型generate zone radbric size 5 5 5 5 &p0 0 0 0 p1 50 0 0 p2 0 0 -20 p3 0 40 0 p4 50 0 -20 &p5 0 40 -20 p6 50 40 0 p7 50 40 -20 p8 20 0 0 p9 0 0 -8.5 &p10 0 15 0 p11 20 0 -8.5 p12 0 15 -8.5 p13 20 13 0 p14 20 13 -8.5 fill group kengnei generate zone radbrick size 5 5 5 5 &p0 0 0 0 p1 0 -40 0 p2 0 0 -20 p3 50 0 0 &p4 0 -40 -20 p5 50 0 -20 p6 50 -40 0 p7 50 -40 -20 &p8 0 -13 0 p9 0 0 -8.5 p10 20 0 0 p11 0 -13 -8.5 p12 20 0 -8.5 &p13 10 -13 0 p14 10 -13 -8.5 fill group kengnei;---建立对称单元gen zone reflect dd 90 dip 90 origin 0 0 0;---连接attach face rang z -8.6 -8.4;---基坑参数设定莫尔库仑参数attach face ran y -.1 .1attach face ran x -.1 .1group section1 range z -2,0 group kengneigroup section2 range z -4,-2 group kengneigroup section3 range z -6,-4 group kengneigroup section4 range z -8,-6 group kengneigroup section5 range z -8.5,-8 group kengnei;---基坑土体设置材料模型参数model mohrprop bulk 8e9 shear 5e9 fric 35 coh 1e10 tens 1e10;---设置初始应力状态ini xvel=0.0ini density 1600 range z -1,0ini density 2000 range z -6,-1ini density 2200 range z -20,-6ini szz 0. grad 0 0 16000 range z -1 0ini szz 0. grad 0 0 20000 range z -6 -1ini szz 0. grad 0 0 22000 range z -20 -6set gravity 0,0,-10.0;---设置边界条件fix x range x -50.1 -49.9fix x range x 50.1 49.9fix z range z -19.9 -20.1fix y range y 40.1 39.9fix y range y -40.1 -39.9fix z range group shellfix x range group shell x -0.1 0.1fix x range group shell x -22.6 -22.4fix x range group shell x -20.6 -20.4fix x range group shell x 22.6 22.4fix x range group shell x 20.4 20.6fix x y z range group shell y -13.6 -13.4fix x y z range group shell y 17.4 17.6;---监测变量hist nstep=5hist unbalset mechanical force 50solveplothistory 1plotcreate trenchadd contour dispadd axes blackshowfree x range x -0.1 0.1free y range y -0.1 0.1set large;---建立前排桩模型sel pile id 1 begin 20.5 0 -0.6 end 20.5 0 -17.1 nseg 10sel pile id 2 begin 20.5 2 -0.6 end 20.5 2 -17.1 nseg 10sel pile id 3 begin 20.5 4 -0.6 end 20.5 4 -17.1 nseg 10sel pile id 4 begin 20.5 6 -0.6 end 20.5 6 -17.1 nseg 10sel pile id 5 begin 20.5 8 -0.6 end 20.5 8 -17.1 nseg 10sel pile id 6 begin 20.5 10 0.6 end 20.5 10 -17.1 nseg 10sel pile id 7 begin 20.5 12 -0.6 end 20.5 12 -17.1 nseg 10sel pile id 8 begin 20.5 13 -0.6 end 20.5 13 -17.1 nseg 10sel pile id 9 begin 19.0 13.2 -0.6 end 19.0 13.2 -17.1 nseg 10sel pile id 10 begin 17.5 13.4 -0.6 end 17.5 13.4 -17.1 nseg 10sel pile id 11 begin 16.0 13.6 -0.6 end 16.0 13.6 -17.1 nseg 10sel pile id 12 begin 14.5 13.8 -0.6 end 14.5 13.8 -17.1 nseg 10sel pile id 13 begin 13.0 14.0 -0.6 end 13.0 14.0 -17.1 nseg 10sel pile id 14 begin 11.5 14.2 -0.6 end 11.5 14.2 -17.1 nseg 10sel pile id 15 begin 10.0 14.4 -0.6 end 10.0 14.4 -17.1 nseg 10sel pile id 16 begin 8.5 14.6 -0.6 end 8.5 14.6 -17.1 nseg 10sel pile id 17 begin 7.0 14.8 -0.6 end 7.0 14.8 -17.1 nseg 10sel pile id 18 begin 5.5 15.0 -0.6 end 5.5 15.0 -17.1 nseg 10sel pile id 19 begin 4.0 15.2 -0.6 end 4.0 15.2 -17.1 nseg 10sel pile id 20 begin 2.5 15.4 -0.6 end 2.5 15.4 -17.1 nseg 10sel pile id 21 begin 1.0 15.5 -0.6 end 1.0 15.5 -17.1 nseg 10sel pile id 22 begin 0.0 15.5 -0.6 end 0.0 15.5 -17.1 nseg 10sel pile id 23 begin -1.0 15.5 -0.6 end -1.0 15.5 -17.1 nseg 10 sel pile id 24 begin -2.5 15.4 -0.6 end -2.5 15.4 -17.1 nseg 10 sel pile id 25 begin -4.0 15.2 -0.6 end -4.0 15.2 -17.1 nseg 10 sel pile id 26 begin -5.5 15.0 -0.6 end -5.5 15.0 -17.1 nseg 10 sel pile id 27 begin -7.0 14.8 -0.6 end -7.0 14.8 -17.1 nseg 10sel pile id 28 begin -8.5 14.6 -0.6 end -8.5 14.6 -17.1 nseg 10sel pile id 29 begin -10.0 14.4 -0.6 end -10.0 14.4 -17.1 nseg 10 sel pile id 30 begin -11.5 14.2 -0.6 end -11.0 14.2 -17.1 nseg 10 sel pile id 31 begin -13.0 14.0 -0.6 end -13.0 14.0 -17.1 nseg 10 sel pile id 32 begin -14.5 13.8 -0.6 end -14.5 13.8 -17.1 nseg 10sel pile id 33 begin -16.0 13.6 -0.6 end -16.0 13.6 -17.1 nseg 10 sel pile id 34 begin -17.5 13.4 -0.6 end -17.5 13.4 -17.1 nseg 10sel pile id 35 begin -19.0 13.2 -0.6 end -19.0 13.2 -17.1 nseg 10 sel pile id 36 begin -20.5 13.0 -0.6 end -20.5 13.0 -17.1 nseg 10 sel pile id 37 begin -20.5 12.0 -0.6 end -20.5 12.0 -17.1 nseg 10 sel pile id 38 begin -20.5 10.0 -0.6 end -20.5 10.0 -17.1 nseg 10 sel pile id 39 begin -20.5 8.0 -0.6 end -20.5 8.0 -17.1 nseg 10 sel pile id 40 begin -20.5 6.0 -0.6 end -20.5 6.0 -17.1 nseg 10 sel pile id 41 begin -20.5 4.0 -0.6 end -20.5 4.0 -17.1 nseg 10sel pile id 42 begin -20.5 2.0 -0.6 end -20.5 2.0 -17.1 nseg 10 sel pile id 43 begin -20.5 0.0 -0.6 end -20.5 0.0 -17.1 nseg 10 sel pile id 44 begin -19.2 -1.7 -0.6 end -19.2 -1.7 -17.1 nseg 10 sel pile id 45 begin -17.9 -3.4 -0.6 end -17.9 -3.4 -17.1 nseg 10 sel pile id 46 begin -16.6 -5.1 -0.6 end -16.6 -5.1 -17.1 nseg 10 sel pile id 47 begin -15.3 -6.8 -0.6 end -15.3 -6.8 -17.1 nseg 10 sel pile id 48 begin -14.0 -8.5 -0.6 end -14.0 -8.5 -17.1 nseg 10 sel pile id 49 begin -12.7 -10.2 -0.6 end -12.7 -10.2 -17.1 nseg 10 sel pile id 50 begin -11.4 -11.9 -0.6 end -11.4 -11.9 -17.1 nseg 10 sel pile id 51 begin -10.5 -13.5 -0.6 end -10.5 -13.5 -17.1 nseg 10 sel pile id 52 begin -8.5 -13.5 -0.6 end -8.5 -13.5 -17.1 nseg 10 sel pile id 53 begin -6.5 -13.5 -0.6 end -6.5 -13.5 -17.1 nseg 10 sel pile id 54 begin -4.5 -13.5 -0.6 end -4.5 -13.5 -17.1 nseg 10 sel pile id 55 begin -2.5 -13.5 -0.6 end -2.5 -13.5 -17.1 nseg 10 sel pile id 56 begin -0.5 -13.5 -0.6 end -0.5 -13.5 -17.1 nseg 10 sel pile id 57 begin 0.5 -13.5 -0.6 end 0.5 -13.5 -17.1 nseg 10sel pile id 58 begin 2.5 -13.5 -0.6 end 2.5 -13.5 -17.1 nseg 10sel pile id 59 begin 4.5 -13.5 -0.6 end 4.5 -13.5 -17.1 nseg 10sel pile id 60 begin 6.5 -13.5 -0.6 end 6.5 -13.5 -17.1 nseg 10sel pile id 61 begin 8.5 -13.5 -0.6 end 8.5 -13.5 -17.1 nseg 10sel pile id 62 begin 10.5 -13.5 -0.6 end 10.5 -13.5 -17.1 nseg 10 sel pile id 63 begin 11.4 -11.9 -0.6 end 11.4 -11.9 -17.1 nseg 10 sel pile id 64 begin 12.7 -10.2 -0.6 end 12.7 -10.2 -17.1 nseg 10 sel pile id 65 begin 14.0 -8.5 -0.6 end 14.0 -8.5 -17.1 nseg 10sel pile id 66 begin 15.3 -6.8 -0.6 end 15.3 -6.8 -17.1 nseg 10sel pile id 67 begin 16.6 -5.1 -0.6 end 16.6 -5.1 -17.1 nseg 10sel pile id 68 begin 17.9 -3.4 -0.6 end 17.9 -3.4 -17.1 nseg 10sel pile id 69 begin 19.2 -1.7 -0.6 end 19.2 -1.7 -17.1 nseg 10sel pile id 70 begin 20.5 0.0 -0.6 end 20.5 0.0 -17.1 nseg 10sel pile id 71 begin 22.5 1 -0.6 end 22.5 1 -10.6 nseg 10sel pile id 72 begin 22.5 3 -0.6 end 22.5 3 -10.6 nseg 10sel pile id 73 begin 22.5 5 -0.6 end 22.5 5 -10.6 nseg 10sel pile id 74 begin 22.5 7 -0.6 end 22.5 7 -10.6 nseg 10sel pile id 75 begin 22.5 9 -0.6 end 22.5 9 -10.6 nseg 10sel pile id 76 begin 22.5 11 -0.6 end 22.5 11 -10.6 nseg 10sel pile id 77 begin 22.5 13 -0.6 end 22.5 13 -10.6 nseg 10sel pile id 78 begin 22.5 15 -0.6 end 22.5 15 -10.6 nseg 10sel pile id 79 begin 20.7 15.2 -0.6 end 20.7 15.2 -10.6 nseg 10 sel pile id 80 begin 19.9 15.4 -0.6 end 19.9 15.4 -10.6 nseg 10 sel pile id 81 begin 18.1 15.6 -0.6 end 18.1 15.6 -10.6 nseg 10 sel pile id 82 begin 16.3 15.8 -0.6 end 16.3 15.8 -10.6 nseg 10 sel pile id 83 begin 14.5 16.0 -0.6 end 14.5 16.0 -10.6 nseg 10 sel pile id 84 begin 12.7 16.2 -0.6 end 12.7 16.2 -10.6 nseg 10sel pile id 85 begin 10.9 16.4 -0.6 end 10.9 16.4 -10.6 nseg 10sel pile id 86 begin 9.1 16.6 -0.6 end 9.1 16.6 -10.6 nseg 10sel pile id 87 begin 7.3 16.8 -0.6 end 7.3 16.8 -10.6 nseg 10sel pile id 88 begin 5.5 17.0 -0.6 end 5.5 17.0 -10.6 nseg 10sel pile id 89 begin 3.7 17.2 -0.6 end 3.7 17.2 -10.6 nseg 10sel pile id 90 begin 1.9 17.4 -0.6 end 1.9 17.4 -10.6 nseg 10sel pile id 91 begin 0 17.5 -0.6 end 0 17.5 -10.6 nseg 10sel pile id 92 begin -1.9 17.4 -0.6 end -1.9 17.4 -10.6 nseg 10sel pile id 93 begin -3.7 17.2 -0.6 end -3.7 17.2 -10.6 nseg 10sel pile id 94 begin -5.5 17.0 -0.6 end -5.5 17.0 -10.6 nseg 10sel pile id 95 begin -7.3 16.8 -0.6 end -7.3 16.8 -10.6 nseg 10sel pile id 96 begin -9.1 16.6 -0.6 end -9.1 16.6 -10.6 nseg 10sel pile id 97 begin -10.9 16.4 -0.6 end -10.9 16.4 -10.6 nseg 10 sel pile id 98 begin -12.7 16.2 -0.6 end -12.7 16.2 -10.6 nseg 10 sel pile id 99 begin -14.5 16.0 -0.6 end -14.5 16.0 -10.6 nseg 10 sel pile id 100 begin -16.3 15.8 -0.6 end -16.3 15.8 -10.6 nseg 10 sel pile id 101 begin -18.1 15.6 -0.6 end -18.1 15.6 -10.6 nseg 10 sel pile id 102 begin -19.9 15.4 -0.6 end -19.9 15.4 -10.6 nseg 10 sel pile id 103 begin -20.7 15.2 -0.6 end -20.7 15.2 -10.6 nseg 10 sel pile id 104 begin -22.5 15.0 -0.6 end -22.5 15.0 -10.6 nseg 10 sel pile id 105 begin -22.5 13.0 -0.6 end -22.5 13.0 -10.6 nseg 10 sel pile id 106 begin -22.5 11.0 -0.6 end -22.5 11.0 -10.6 nseg 10 sel pile id 107 begin -22.5 9.0 -0.6 end -22.5 9.0 -10.6 nseg 10sel pile id 108 begin -22.5 7.0 -0.6 end -22.5 7.0 -10.6 nseg 10sel pile id 109 begin -22.5 5.0 -0.6 end -22.5 5.0 -10.6 nseg 10sel pile id 110 begin -22.5 3.0 -0.6 end -22.5 3.0 -10.6 nseg 10sel pile id 111 begin -22.5 1.0 -0.6 end -22.5 1.0 -10.6 nseg 10sel pile id 112 begin -22.5 -1.0 -0.6 end -22.5 -1.0 -10.6 nseg 10 sel pile id 113 begin -21.3 -2.8 -0.6 end -21.3 -2.8 -10.6 nseg 10 sel pile id 114 begin -20.1 -4.6 -0.6 end -20.1 -4.6 -10.6 nseg 10 sel pile id 115 begin -18.9 -6.4 -0.6 end -18.9 -6.4 -10.6 nseg 10 sel pile id 116 begin -17.7 -8.2 -0.6 end -17.7 -8.2 -10.6 nseg 10 sel pile id 117 begin -16.5 -10.0 -0.6 end -16.5 -10.0 -10.6 nseg 10 sel pile id 118 begin -15.3 -11.8 -0.6 end -15.3 -11.8 -10.6 nseg 10 sel pile id 119 begin -14.1 -13.6 -0.6 end -14.1 -13.6 -10.6 nseg 10 sel pile id 120 begin -12.5 -15.5 -0.6 end -12.5 -15.5 -10.6 nseg 10 sel pile id 121 begin -11.5 -15.5 -0.6 end -11.5 -15.5 -10.6 nseg 10 sel pile id 122 begin -9.5 -15.5 -0.6 end -9.5 -15.5 -10.6 nseg 10 sel pile id 123 begin -7.5 -15.5 -0.6 end -7.5 -15.5 -10.6 nseg 10 sel pile id 124 begin -5.5 -15.5 -0.6 end -5.5 -15.5 -10.6 nseg 10 sel pile id 125 begin -3.5 -15.5 -0.6 end -3.5 -15.5 -10.6 nseg 10 sel pile id 126 begin -1.5 -15.5 -0.6 end -1.5 -15.5 -10.6 nseg 10 sel pile id 127 begin 1.5 -15.5 -0.6 end 1.5 -15.5 -10.6 nseg 10 sel pile id 128 begin 3.5 -15.5 -0.6 end 3.5 -15.5 -10.6 nseg 10sel pile id 129 begin 5.5 -15.5 -0.6 end 5.5 -15.5 -10.6 nseg 10 sel pile id 130 begin 7.5 -15.5 -0.6 end 7.5 -15.5 -10.6 nseg 10 sel pile id 131 begin 9.5 -15.5 -0.6 end 9.5 -15.5 -10.6 nseg 10 sel pile id 132 begin 11.5 -15.5 -0.6 end 11.5 -15.5 -10.6 nseg 10 sel pile id 133 begin 12.5 -15.5 -0.6 end 12.5 -15.5 -10.6 nseg 10 sel pile id 134 begin 14.1 -13.6 -0.6 end 14.1 -13.6 -10.6 nseg 10 sel pile id 135 begin 15.3 -11.8 -0.6 end 15.3 -11.8 -10.6 nseg 10 sel pile id 136 begin 16.5 -10.0 -0.6 end 16.5 -10.0 -10.6 nseg 10 sel pile id 137 begin 17.7 -8.2 -0.6 end 17.7 -8.2 -10.6 nseg 10 sel pile id 138 begin 18.9 -6.4 -0.6 end 18.9 -6.4 -10.6 nseg 10sel pile id 139 begin 20.1 -4.6 -0.6 end 20.1 -4.6 -10.6 nseg 10sel pile id 140 begin 21.3 -2.8 -0.6 end 21.3 -2.8 -10.6 nseg 10sel pile id 141 begin 22.5 -1.0 -0.6 end 22.5 -1.0 -10.6 nseg 10sel pile prop emod=8.0e10 nu=0.30 xcarea=0.7854 &xcj=0.0 xciy=0.0 xciz=0.0 &per=3.14 &cs_sk=1.3e11 cs_scoh=0.0 cs_sfric=10.0 &cs_nk=1.3e11 cs_ncoh=0.0 cs_nfric=0.0 cs_ngap=off;---建立后排桩模型attach face ran z -0.5 -0.7;-----第一阶段：开挖隧道区域1mod null range group section1set largeinit xdisp 0.0 ydisp 0.0 zdisp 0.0step 2000init xdisp 0.0 ydisp 0.0 zdisp 0.0;hist id=2 gp zvel 0 1 1;hist id=3 gp zvel 0 1 1;;mod null range group section2;solve;save shell-tun2.sav;mod null range group section3;solve;save shell-tun3.sav;mod null range group section4;solve;save shell-tun4.sav;mod null range group section5;solve;save shell-tun5.sav;set large;init xdisp 0.0 ydisp 0.0 zdisp 0.0;step 2000plot create dispcontplot copy gravV dispcont settingsplot add contour disp plane behind shade onplot add axesplot show;plot set rot 40 0 20;plot sh;plot add blo mo me blue ;--显示几何模型;plo add blo group green blue green ;--显示group ;plo add sel geom shell black ;---显示结构体;plo add sel geom pile red ;---显示结构体;plot add axes blackhist unbal。