1、读位移双侧壁导坑法静态模拟(三维)/prep7et,1,solid45et,2,shell63mp,ex,1,1.5e9 !所有岩石的 mpmp,dens,1,1850mp,prxy,1,0.0.4tb,dp,1tbdata,1,0.125e6,45mp,ex,2,3.65e9 !加固圈岩石的 mpmp,dens,2,2150mp,prxy,2,0.32tb,dp,2tbdata,1,0.45e6,55mp,ex,3,3.1e10 !隧道钢筋混凝土的 mpmp,dens,3,2500mp,prxy,3,0.2real,1,0.25,0.25,0.25,0.25mp,ex,4,1.5e9 !所有
2、岩石的 mpmp,dens,4,1850mp,prxy,4,0.4tb,dp,4tbdata,4,0.125e6,45mp,ex,5,1.5e9 !所有岩石的 mpmp,dens,5,1850mp,prxy,5,0.4tb,dp,5tbdata,5,0.125e6,45mp,ex,6,1.5e9 !所有岩石的 mpmp,dens,6,1850mp,prxy,6,0.4tb,dp,6tbdata,6,0.125e6,45mp,ex,7,3.1e10 !地表混凝土 mpmp,dens,7,2500mp,prxy,7,0.2real,2,0.85,0.85,0.85,0.85/pnum,kp,1/p
3、num,line,1/pnum,area,1/pnum,volu,1k,1,-47.16,12.17641,0k,2,-31,12.17641,0k,3,-17.16,12.17641,0k,4,-9.44112,12.17641,0k,5,-6.654410,12.17641,0k,6,-1.75516,12.17641,0k,7,0,12.17641,0k,8,1.75516,12.17641,0k,9,6.654410,12.17641,0k,10,9.44112,12.17641,0k,11,17.16,12.17641,0k,12,31,12.17641,0k,13,47.16,12
4、.17641,0KGEN,2,1,13,1,0,-1.16949,0KGEN,2,14,26,1,0,-10.8,0KGEN,2,27,39,1,0,-3.64187,0KGEN,2,40,52,1,0,-7.17813,0KGEN,2,53,65,1,0,-24,0!生成线*do,i,0,5,1*do,j,1,12,1l,i*13+j,i*13+j+1*enddo*enddo*do,i,1,65,1l,i,i+13*enddo!线划分网格*do,i,1,137,1lsel,s,line,ilesize,all,3*enddoALLSEL,ALLldele,76ldele,77ldele,78
5、ldele,80ldele,81ldele,82lcomb,3,4lcomb,3,5lcomb,3,6lcomb,7,8lcomb,7,9lcomb,7,10ldele,89ldele,90ldele,91ldele,92ldele,93ldele,94ldele,95ldele,102ldele,103ldele,104ldele,105ldele,106ldele,107ldele,108ldele,115ldele,116ldele,117ldele,118ldele,119ldele,120ldele,121ldele,27ldele,28ldele,29ldele,30ldele,3
6、1ldele,32ldele,33ldele,34ldele,39ldele,40ldele,41ldele,42ldele,43ldele,44ldele,45ldele,46lcomb,15,16lcomb,15,17lcomb,15,18lcomb,19,20lcomb,19,21lcomb,19,22lsel,s,line,3lsel,a,line,7lsel,a,line,15lsel,a,line,19lesize,all,3!生成面ALLSEL,ALL FLST,2,4,4 FITEM,2,1 FITEM,2,73 FITEM,2,13 FITEM,2,74 AL,P51X FL
7、ST,2,4,4 FITEM,2,2 FITEM,2,74 FITEM,2,14 FITEM,2,75 AL,P51X FLST,2,4,4 FITEM,2,75 FITEM,2,3 FITEM,2,79 FITEM,2,15 AL,P51X FLST,2,4,4 FITEM,2,79 FITEM,2,7 FITEM,2,19 FITEM,2,83 AL,P51X FLST,2,4,4 FITEM,2,11 FITEM,2,83 FITEM,2,23 FITEM,2,84 AL,P51X FLST,2,4,4 FITEM,2,84 FITEM,2,12 FITEM,2,85 FITEM,2,2
8、4 AL,P51X FLST,2,4,4 FITEM,2,98 FITEM,2,24 FITEM,2,97 FITEM,2,36 AL,P51X FLST,2,4,4 FITEM,2,23 FITEM,2,96 FITEM,2,97 FITEM,2,35 AL,P51X FLST,2,4,4 FITEM,2,35 FITEM,2,109 FITEM,2,47 FITEM,2,110 AL,P51X FLST,2,4,4 FITEM,2,110 FITEM,2,36 FITEM,2,111 FITEM,2,48 AL,P51X FLST,2,4,4 FITEM,2,48 FITEM,2,124
9、FITEM,2,60 FITEM,2,123 AL,P51X FLST,2,4,4 FITEM,2,47 FITEM,2,122 FITEM,2,59 FITEM,2,123 AL,P51X FLST,2,4,4 FITEM,2,60 FITEM,2,137 FITEM,2,136 FITEM,2,72 AL,P51X FLST,2,4,4 FITEM,2,59 FITEM,2,135 FITEM,2,136 FITEM,2,71 AL,P51X FLST,2,4,4 FITEM,2,58 FITEM,2,134 FITEM,2,135 FITEM,2,70 AL,P51X FLST,2,4,
10、4 FITEM,2,57 FITEM,2,134 FITEM,2,133 FITEM,2,69 AL,P51X FLST,2,4,4 FITEM,2,56 FITEM,2,132 FITEM,2,133 FITEM,2,68 AL,P51X FLST,2,4,4 FITEM,2,55 FITEM,2,131 FITEM,2,132 FITEM,2,67 AL,P51X FLST,2,4,4 FITEM,2,54 FITEM,2,130 FITEM,2,131 FITEM,2,66 AL,P51X FLST,2,4,4 FITEM,2,53 FITEM,2,130 FITEM,2,129 FIT
11、EM,2,65 AL,P51X FLST,2,4,4 FITEM,2,52 FITEM,2,128 FITEM,2,129 FITEM,2,64 AL,P51X FLST,2,4,4 FITEM,2,51 FITEM,2,127 FITEM,2,128 FITEM,2,63 AL,P51X FLST,2,4,4 FITEM,2,50 FITEM,2,127 FITEM,2,126 FITEM,2,62 AL,P51X FLST,2,4,4 FITEM,2,49 FITEM,2,126 FITEM,2,125 FITEM,2,61 AL,P51X FLST,2,4,4 FITEM,2,112 F
12、ITEM,2,37 FITEM,2,113 FITEM,2,49 AL,P51X FLST,2,4,4 FITEM,2,38 FITEM,2,114 FITEM,2,113 FITEM,2,50 AL,P51X FLST,2,4,4 FITEM,2,26 FITEM,2,101 FITEM,2,100 FITEM,2,38 AL,P51X FLST,2,4,4 FITEM,2,25 FITEM,2,99 FITEM,2,37 FITEM,2,100 AL,P51X FLST,2,4,4 FITEM,2,13 FITEM,2,86 FITEM,2,25 FITEM,2,87 AL,P51X FL
13、ST,2,4,4 FITEM,2,14 FITEM,2,88 FITEM,2,87 FITEM,2,26 AL,P51X k,100,0,7.17641k,101,6.21496,3.58821k,102,6.4108,3.22396k,103,7.11501,-0.27906k,104,7.05067,-0.86479k,105,6.65441,-3.44187k,106,0,10.17641k,107,8.81303,5.08821k,108,10.11959,-0.27906k,109,9.44112,-3.44187k,110,4.23861,3.17641k,111,1.75516,
14、-0.27906k,112,1.75516,-4.68728k,113,0,-4.77679k,114,0,-0.27906k,115,-0.34113,3.17641k,120,0,0k,121,0.3,0.17321k,122,5.40086,-0.27906k,123,0,12.47633larc,100,101,120,7.17641larc,101,102,120,6.83larc,102,103,121,6.83larc,103,104,121,6.83larc,106,107,120,10.17641larc,107,108,121,9.83larc,108,109,121,9.
15、83larc,105,112,123,17.25312larc,112,113,123,17.25312larc,102,110,122,3.6457larc,110,111,122,3.6457l,111,112l,103,111l,111,114l,105,109l,104,105FLST,3,16,4,ORDE,10 FITEM,3,4 FITEM,3,-6 FITEM,3,8 FITEM,3,-10 FITEM,3,16 FITEM,3,-18 FITEM,3,20 FITEM,3,-22 FITEM,3,27 FITEM,3,-30 LSYMM,X,P51X, , , ,0,0 lc
16、omb,4,5lcomb,31,32lcomb,34,80lcomb,8,30l,110,115l,115,79l,100,115l,115,114l,114,113LSTR, 16, 17 LSTR, 20, 106 LSTR, 24, 107 LSTR, 37, 108 LSTR, 50, 109 LSTR, 109, 62 LSTR, 105, 61 LSTR, 112, 60 LSTR, 113, 59 LSTR, 22, 58 LSTR, 21, 57 LSTR, 19, 56 LSTR, 19, 42 LSTR, 29, 18 LSTR, 18, 8 LSTR, 17, 6 LST
17、R, 106, 100 LSTR, 107, 102 LSTR, 108, 103 NUMMRG,NODE, , , ,LOW NUMCMP,NODE NUMMRG,NODE, , , ,LOW NUMCMP,ALL NUMMRG,ALL, , , ,LOWNUMCMP,ALL ALLSEL,ALL !隧道线的网分ALLSEL,ALL FLST,5,52,4,ORDE,22 FITEM,5,4 FITEM,5,-6 FITEM,5,8 FITEM,5,-10 FITEM,5,16 FITEM,5,-18 FITEM,5,20 FITEM,5,-22 FITEM,5,27 FITEM,5,-34
18、 FITEM,5,39 FITEM,5,-46 FITEM,5,76 FITEM,5,-78 FITEM,5,80 FITEM,5,-82 FITEM,5,89 FITEM,5,-95 FITEM,5,102 FITEM,5,-108FITEM,5,115 FITEM,5,-118CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,_Y LESIZE,_Y1, , ,3, , , , ,1 !隧道线组成面FLST,2,4,4 FITEM,2,82 FITEM,2,15 FITEM,2,89 FITEM,2,39 AL,P51X FLST,2,4,4 FIT
19、EM,2,89 FITEM,2,9 FITEM,2,90 FITEM,2,19 AL,P51X FLST,2,4,4 FITEM,2,90 FITEM,2,96 FITEM,2,91 FITEM,2,10 AL,P51X FLST,2,4,4 FITEM,2,91 FITEM,2,109 FITEM,2,92 FITEM,2,16 AL,P51X FLST,2,4,4 FITEM,2,92 FITEM,2,119 FITEM,2,58 FITEM,2,93 AL,P51X FLST,2,4,4 FITEM,2,93 FITEM,2,94 FITEM,2,29 FITEM,2,57 AL,P51
20、X FLST,2,4,4 FITEM,2,17 FITEM,2,94 FITEM,2,56 FITEM,2,95 AL,P51X FLST,2,4,4 FITEM,2,18 FITEM,2,95 FITEM,2,55 FITEM,2,102 AL,P51X FLST,2,4,4 FITEM,2,43 FITEM,2,102 FITEM,2,103 FITEM,2,54 AL,P51X FLST,2,4,4 FITEM,2,42 FITEM,2,103 FITEM,2,53 FITEM,2,104 AL,P51X FLST,2,4,4 FITEM,2,78 FITEM,2,104 FITEM,2
21、,52 FITEM,2,105 AL,P51X FLST,2,4,4 FITEM,2,106 FITEM,2,105 FITEM,2,51 FITEM,2,114 AL,P51X FLST,2,4,4 FITEM,2,101 FITEM,2,106 FITEM,2,107 FITEM,2,41 AL,P51X FLST,2,4,4 FITEM,2,88 FITEM,2,107 FITEM,2,40 FITEM,2,82 AL,P51X FLST,2,4,4 FITEM,2,78 FITEM,2,41 FITEM,2,34 FITEM,2,108 AL,P51X FLST,2,4,4 FITEM
22、,2,40 FITEM,2,108 FITEM,2,33 FITEM,2,115 AL,P51X FLST,2,4,4 FITEM,2,39 FITEM,2,31 FITEM,2,115 FITEM,2,116 AL,P51X FLST,2,4,4 FITEM,2,116 FITEM,2,9 FITEM,2,117 FITEM,2,4 AL,P51X FLST,2,4,4 FITEM,2,117 FITEM,2,10 FITEM,2,118 FITEM,2,6 AL,P51X FLST,2,4,4 FITEM,2,118 FITEM,2,16 FITEM,2,29 FITEM,2,8 AL,P
23、51X FLST,2,4,4 FITEM,2,44 FITEM,2,45 FITEM,2,76 FITEM,2,33 AL,P51X FLST,2,4,4 FITEM,2,76 FITEM,2,46 FITEM,2,34 FITEM,2,42 AL,P51X FLST,2,4,4 FITEM,2,21 FITEM,2,6 FITEM,2,20 FITEM,2,27 AL,P51X FLST,2,4,4 FITEM,2,27 FITEM,2,8 FITEM,2,17 FITEM,2,22 AL,P51X FLST,2,4,4 FITEM,2,31 FITEM,2,32 FITEM,2,30 FI
24、TEM,2,44 AL,P51X FLST,2,4,4 FITEM,2,4 FITEM,2,32 FITEM,2,5 FITEM,2,20 AL,P51X FLST,2,4,4 FITEM,2,30 FITEM,2,45 FITEM,2,80 FITEM,2,77 AL,P51X FLST,2,4,4 FITEM,2,80 FITEM,2,5 FITEM,2,21 FITEM,2,28 AL,P51X FLST,2,4,4 FITEM,2,77 FITEM,2,46 FITEM,2,81 FITEM,2,43 AL,P51X FLST,2,4,4 FITEM,2,81 FITEM,2,28 F
25、ITEM,2,22 FITEM,2,18 AL,P51X k,160,k,161,-5l,160,161vdra,all,135!隧道围岩单元allstype,1 mat,1 vsel,s,volu,allvmesh,all !隧道衬砌的壳单元allstype,2MAT,3real,1asel,s,216 !1号圆asel,a,214 !2号圆asel,a,200asel,a,225 !3号圆asel,a,228 !4号圆asel,a,191asel,a,219 !5号圆加底板asel,a,222 !6号圆加底板asel,a,197 !8号圆加底板asel,a,194amesh,allR,1,
26、0.25,0.25,0.25,0.25, , ,RMORE, , , ,RMORE RMORE, ,!施加边界条件及重力加速度nsel,s,loc,y,-34.61308 d,all,uy,0 allsnsel,s,loc,x,-47.16 nsel,a,loc,x,47.16 d,all,ux,0 allsnsel,s,loc,z,0 nsel,a,loc,z,-5d,all,uz,0 allsacel,9.8!施加水平地应力及房屋荷载FINISH /SOLFLST,2,2,5,ORDE,2 FITEM,2,66 FITEM,2,78 /GO SFA,P51X,1,PRES,70000 FL
27、ST,2,1,5,ORDE,1 FITEM,2,64 /GO SFA,P51X,1,PRES,100000 FLST,2,1,5,ORDE,1 FITEM,2,169 /GO SFA,P51X,1,PRES,200000 FLST,2,1,5,ORDE,1 FITEM,2,166 /GO SFA,P51X,1,PRES,300000 FLST,2,1,5,ORDE,1 FITEM,2,154 /GO SFA,P51X,1,PRES,400000 FLST,2,1,5,ORDE,1 FITEM,2,152 /GO SFA,P51X,1,PRES,500000 FLST,2,1,5,ORDE,1
28、FITEM,2,83 /GO SFA,P51X,1,PRES,100000 FLST,2,1,5,ORDE,1 FITEM,2,86 /GO SFA,P51X,1,PRES,200000 FLST,2,1,5,ORDE,1 FITEM,2,98 /GO SFA,P51X,1,PRES,300000 FLST,2,1,5,ORDE,1 FITEM,2,100 /GO SFA,P51X,1,PRES,400000 FLST,2,1,5,ORDE,1 FITEM,2,107 /GO SFA,P51X,1,PRES,500000 !杀死隧道壳单元esel,s,mat,3ekill,allalls!求初
29、始应力场/solutime,1E-6antype,static nsub,5,100,1nropt,fullsolve!求第一、二次开挖所加节点力vsel,s,45vsel,a,46vsel,a,55vsel,a,57vsel,a,59vsel,a,40eslv,esel,u,mat,3nsle,!nforce,all!工况 1-1/solutime,2fdele,allvsel,s,51vsel,a,52eslv,ekill,allalls!释放 50%荷载(加节点反力)solve!工况 1-2/solutime,3fdele,allvsel,s,46!置换加固圈参数vsel,a,45esl
30、v,esel,u,mat,3mpchg,2,allallsasel,s,216!激活壳单元asel,a,214asel,a,200esla,ealive,allallssolve!求第三、四次开挖所加节点力vsel,s,56vsel,a,58vsel,a,60vsel,a,37vsel,a,49vsel,a,50eslv,esel,u,mat,3nsle,!nforce,all!工况 2-1/solutime,4fdele,allvsel,s,53vsel,a,54eslv,ekill,allalls!释放 50%荷载(加节点反力)solve!工况 2-2/solutime,5fdele,al
31、lvsel,s,49!置换加固圈参数vsel,a,50eslv,esel,u,mat,3mpchg,2,allallsasel,s,225!激活壳单元asel,a,228asel,a,191esla,ealive,allallssolve!求第五、六次开挖所加节点力vsel,s,47vsel,a,48vsel,a,57vsel,a,58eslv,esel,u,mat,3nsle,!nforce,all!工况 3-1/solutime,6fdele,allvsel,s,55vsel,a,56eslv,ekill,allalls!释放 50%荷载(加节点反力)solve!工况 3-2/soluti
32、me,7fdele,allvsel,s,47!置换加固圈参数vsel,a,48eslv,esel,u,mat,3mpchg,2,allallsasel,s,219!激活壳单元asel,a,222esla,ealive,allallssolve!求第七、八次开挖所加节点力vsel,s,39vsel,a,38eslv,esel,u,mat,3nsle,!nforce,all!工况 3-1/solutime,8fdele,allvsel,s,57vsel,a,58vsel,a,59vsel,a,60eslv,ekill,allalls!释放 50%荷载(加节点反力)solve!工况 3-2/solutime,9fdele,allasel,s,197!激活壳单元asel,a,194esla,ealive,allalls
