高速公路隧道施工力学数值模拟分析实例

/NOPR/NOPRKEYW,PR_STRUC,1/PREP7ET,1,BEAM3ET,2,PLANE42KEYOPT,2,3,3R,1,0.15,0.000281,0.15MP,EX,1,23E9MP,PRXY,1,0.2MP,DENS,1,2500MP,EX,2,8E9MP,PRXY,2,0.28MP,DENS,2,2400TB,DP,2TBMODIF,1,1,1.2E6TBMODIF,1,2,42MP,EX,3,16E9MP,PRXY,3,0.25MP,DENS,3,2500TB,DP,3TBMODIF,1,1,2.4E6TBMODIF,1,2,50SAVEK,1,21,0,0K,2,15,-1,0K,3,27,-1,0K,4,0,0,0K,5,0,21,0K,6,21,21,0K,7,42,21,0高速公路隧道力学数值分析!菜单过滤设置!保留结构分析部分菜单!进入前处理器!设置了梁单元类型，模拟喷射混凝土!设置平面实体单元类型，模拟围岩!设置平面应变模式!设置梁单元几何常数!输入弹性模量!输入泊松比!输入密度!输入弹性模量!输入泊松比!输入密度!输入DP模型的相应参数!输入粘聚力!输入内摩擦角!以下输入加固圈（锚杆加固范围）参数K,8,42,0,0K,9,55,0,0K,10,55,21,0K,11,55,56,0K,12,0,56,0K,13,21,56,0K,14,42,56,0K,15,55,-14,0K,16,42,-14,0K,17,27,-14,0K,18,0,-14,0K,19,14.9896,6.0104K,20,27.0104,6.0104!对称复制所有关键点!对称复制所有关键点!合并所有关键点!压缩所有关键点编号NUMMRG,ALL,,,,LOWNUMCMP,ALLSAVEK,,15,-14K,,-15,-14CIRCLE,1,6CIRCLE,21,6CIRCLE,21,8.5CIRCLE,1,8.5LDELE,3,4,1,1LDELE,7,8,1,1LDELE,11,12,1,1LDELE,15,16,1,1NUMMRG,ALL,,,,LOWNUMCMP,ALLSAVE!创建将要划分六段的线L,31,30L,30,12L,12,13L,13,14L,29,28L,31,25L,30,24L,12,5L,13,6L,14,7L,11,10L,25,24L,24,5L,5,6L,6,7L,28,27L,25,26L,25,36L,24,46L,43,21L,35,5L,5,4L,5,19L,6,49L,40,1L,20,7L,7,8L,9,10L,8,48L,4,50L,4,45L,26,47L,27,32L,26,33L,4,18L,8,16L,9,15!创建将要划分三段的线L,29,31L,14,11L,28,25L,7,10L,27,26L,8,9L,33,32L,15,16!绘制将要划分五段的线L,23,34L,22,38L,2,37L,3,17!绘制将要划分二段的线L,44,47L,43,46L,42,45L,50,41L,40,49L,39,48!绘制将要划分一段的线L,23,44L,22,42L,41,2L,39,3!绘制将要划分十二段的线L,44,21L,21,42L,41,1L,1,39!绘制将要划分二十四段的线L,22,23L,34,38L,2,3L,37,17!绘制将要划分八段的线L,33,34L,38,18L,18,37L,17,16SAVE!通过线创建面AL,22,23,35,37,7,8,38,30LSEL,S,LINE,,31,38,7LSEL,A,LINE,,8ASBL,1,ALLLSEL,S,LINE,,31,34,3LSEL,A,LINE,,7,81,74ASBL,3,ALLALLSEL,ALLAL,20,21,30,39,5,6,40,25LSEL,S,LINE,,26,40,14LSEL,A,LINE,,6ASBL,3,ALLLSEL,S,LINE,,26,29,3LSEL,A,LINE,,5,83,78ASBL,6,ALLNUMMRG,ALL,,,,LOWNUMCMP,ALLSAVEALLSEL,ALLAL,2,32,69AL,32,70,1AL,79,80,61,2,60AL,61,81,7,62,1AL,37,60,65,55,77,42AL,69,70,66,73,65AL,73,56,74,55AL,62,36,43,78,56,66AL,50,44,52,43AL,48,35,50,34AL,46,18,48,17AL,11,17,22,16AL,10,16,21,15AL,9,15,20,14AL,8,14,19,13AL,45,13,47,12AL,47,24,49,23AL,4,27,67AL,27,68,3AL,6,82,58,4,57AL,83,5,59,3,58AL,49,41,51,40AL,39,57,63,53,75,41AL,67,68,64,71,63AL,71,54,72,53AL,64,59,38,42,76,54SAVE!单元网格划分ALLSEL,ALLLSEL,S,LINE,,5,44,1LSEL,A,LINE,,79,83,1LSEL,U,LINE,,27,32,5LESIZE,ALL,,,6LSEL,S,LINE,,45,52,1LESIZE,ALL,,,3LSEL,S,LINE,,53,56,1LESIZE,ALL,,,5LSEL,S,LINE,,57,62,1LESIZE,ALL,,,2LSEL,S,LINE,,63,66,1LESIZE,ALL,,,1LSEL,S,LINE,,67,70,1LSEL,A,LINE,,1,4,1LSEL,A,LINE,,27,32,5LESIZE,ALL,,,12LSEL,S,LINE,,71,74,1LESIZE,ALL,,,24LSEL,S,LINE,,75,78,1LESIZE,ALL,,,8SAVE!加固圈单元划分ALLSEL,ALLTYPE,2MAT,3REAL,1LCCAT,79,80LCCAT,81,7LCCAT,6,82LCCAT,5,83AMAP,26,47,46,43,44AMAP,27,46,45,42,43AMAP,9,50,49,40,41AMAP,10,49,48,39,40LSEL,S,LCCATLDELE,ALL!以下划分开挖土体网格ALLSEL,ALLTYPE,2MAT,2REAL,1AMAP,24,44,43,21AMAP,25,43,42,21AMAP,7,41,40,1AMAP,8,40,39,1LCCAT,67,68LCCAT,69,70AMAP,30,44,42,22,23AMAP,12,41,39,3,2LSEL,S,LCCATLDELE,ALL!划分未开挖土体网格ALLSEL,ALLLCCAT,19,20LCCAT,82,83LCCAT,21,22LCCAT,80,81AMAP,6,25,5,35,36AMAP,4,5,7,20,19LSEL,S,LCCATLSELE,ALLAMAP,2,5,19,50,4AMAP,1,20,7,8,48AMAP,13,2,3,17,37AMAP,15,8,9,15,16AMAP,16,7,10,9,8AMAP,17,14,11,10,7AMAP,18,13,14,7,6AMAP,19,12,13,6,5AMAP,20,30,12,5,24AMAP,21,31,30,24,25AMAP,22,29,31,25,28AMAP,23,28,25,26,27AMAP,5,25,36,47,26AMAP,3,35,5,4,45AMAP,28,27,26,33,32AMAP,31,23,22,38,34LCCAT,39,57LCCAT,63,53LCCAT,54,64LCCAT,59,38LCCAT,37,60LCCAT,55,65LCCAT,56,66LCCAT,62,36AMAP,29,26,44,34,33AMAP,32,42,4,18,38AMAP,11,4,41,37,18AMAP,14,39,8,16,17LSEL,S,LCCATLDELE,ALLSAVE!生成初期支护喷射混凝土梁单元TYPE,1MAT,1REAL,1!右隧道E,608,127E,127,129E,129,130E,130,131E,131,132E,132,133E,133,134E,134,135E,135,136E,136,137E,137,138E,138,139E,139,91E,91,92E,92,93E,93,94E,94,95E,95,96E,96,97E,97,98E,98,99E,99,100E,100,101E,101,102E,102,89E,89,609!左隧道E,583,52E,52,53E,53,54E,54,55E,55,56E,56,57E,57,58E,58,59E,59,60E,60,61E,61,62E,62,63E,63,16E,16,17E,17,18E,18,19E,19,20E,20,21E,21,22E,22,23E,23,24E,24,25E,25,26E,26,27E,27,14E,14,584SAVE

!加载与初始地应力的求解!对左右两侧边界施加UX方向的约束，在下边界施加UY方向的约束/SOLALLSEL,ALLNSEL,S,LOC,X,-55.1,-54.9NSEL,A,LOC,X,54.9,55.1D,ALL,UX,0NSEL,S,LOC,Y,-14.1,-13.9D,ALL,UY,0SAVE!施加重力加速度ACEL,0,10,0SAVE!杀死梁单元ESEL,S,TYPE,,1EKILL,ALLSAVE!求解!求解前设置NROPT,FULLDELTIM,100,20,1000ALLSEL,ALLOUTRES,ALL,ALL!初始应力求解SOLVEFINISHSAVE!自重应力场处理/POST1!采用全牛顿—拉普森法进行求解!设置求解次数!采用全牛顿—拉普森法进行求解!设置求解次数!竖向位移PLNSOL,S,2,0,1 !第一主应力PLNSOL,S,3,0,1 !第三主应力!左隧道土体开挖施工!开挖左侧隧道土体FINISH/SOLANTYPE,,REST,1,1,0ASEL,S,AREA,,24,25,1ASEL,A,AREA,,30ESLA,SEKILL,ALL!进行初期支护，即激活喷射混凝土衬砌梁单元ALLSEL,ALLESEL,S,ELEM,,1549,1574,1EKILL,ALL!杀死右隧道喷射混凝土梁单元ALLSEL,ALLESEL,S,ELEM,,1575,1600,1EALIVE,ALLSAVEFINISH!进入后处理器!选择所有内容!进入后处理器!选择所有内容!人工手动选择开挖土体周围的一圈单元!对选择的节点求解节点力/POST1ALLSEL,ALLESEL,ALLESEL,S,LIVEESEL,S,PNSEL,S,PNFORCE,ALLFINISH/SOL

!施加节点力F,14,FX,2.83E+04F,16,FX,7.48E+01F,17,FX,3.27E+03F,18,FX,6.42E+03F,19,FX,9.48E+03F,20,FX,1.24E+04F,21,FX,1.51E+04F,22,FX,1.75E+04F,23,FX,1.99E+04F,24,FX,2.19E+04F,25,FX,2.38E+04F,26,FX,2.51E+04F,27,FX,2.47E+04F,52,FX,-2.82E+04F,53,FX,-2.47E+04F,54,FX,-2.50E+04F,55,FX,-2.36E+04F,56,FX,-2.18E+04F,57,FX,-1.97E+04F,58,FX,-1.74E+04F,59,FX,-1.49E+04F,60,FX,-1.22E+04F,61,FX,-9.30E+03F,62,FX,-6.30E+03F,63,FX,-3.12E+03F,583,FX,-1.39E+04F,584,FX,1.39E+04F,585,FX,4.18E+03F,586,FX,2.19E+03F,587,FX,1.54E+03F,588,FX,1.07E+03F,589,FX,7.81E+02F,590,FX,5.70E+02F,591,FX,4.13E+02F,592,FX,2.88E+02F,593,FX,1.86E+02F,594,FX,9.81E+01F,595,FX,1.85E+01F,596,FX,-5.62E+01F,597,FX,-1.31E+02F,598,FX,-2.11E+02F,599,FX,-2.99E+02F,600,FX,-4.02E+02F,601,FX,-5.28E+02F,602,FX,-6.86E+02F,603,FX,-8.98E+02F,604,FX,-1.19E+03F,605,FX,-1.66E+03F,606,FX,-2.31E+03F,607,FX,-4.31E+03F,14,FY,-2.06E+04F,16,FY,-8.10E+04F,17,FY,-8.04E+04F,18,FY,-7.86E+04F,19,FY,-7.56E+04F,20,FY,-7.08E+04F,21,FY,-6.54E+04F,22,FY,-5.88E+04F,23,FY,-5.12E+04F,24,FY,-4.28E+04F,25,FY,-3.40E+04F,26,FY,-2.55E+04F,27,FY,-1.85E+04F,52,FY,-2.04E+04F,53,FY,-1.83E+04F,54,FY,-2.53E+04F,55,FY,-3.38E+04F,56,FY,-4.27E+04F,57,FY,-5.10E+04F,58,FY,-5.87E+04F,59,FY,-6.54E+04F,60,FY,-7.08E+04F,61,FY,-7.56E+04F,62,FY,-7.86E+04F,63,FY,-8.04E+04F,583,FY,2.86E+04F,584,FY,2.85E+04F,585,FY,6.60E+04F,586,FY,6.24E+04F,587,FY,6.18E+04F,588,FY,6.12E+04F,589,FY,6.12E+04F,590,FY,6.06E+04F,591,FY,6.06E+04F,592,FY,6.06E+04F,593,FY,6.06E+04F,594,FY,6.06E+04F,595,FY,6.06E+04F,596,FY,6.06E+04F,597,FY,6.06E+04F,598,FY,6.06E+04F,599,FY,6.06E+04

F,600,FY,6.06E+04F,601,FY,6.06E+04F,602,FY,6.12E+04F,603,FY,6.12E+04F,604,FY,6.12E+04F,605,FY,6.18E+04F,606,FY,6.24E+04F,607,FY,6.60E+04!左隧道开挖求解ALLSEL,ALLSOLVEFINISHSAVE!左隧道开挖模拟分析后处理/POST1PLNSOL,U,Y,0,1!Y方向的位移PLNSOL,S,Y,0,1!Y方向的应力PLNSOL,S,1,0,1!第一主应力PLNSOL,S,3,0,1!第三主应力PLNSOL,EPPL,EQV,0,1!塑性区ESEL,S,LIVEESEL,R,TYPE,,1ETABLE,,SMISC,6ETABLE,,SMISC,12!6和12表示弯矩ETABLE,,SMISC,1ETABLE,,SMISC,7!1和7表示轴力ETABLE,,SMISC,2ETABLE,,SMISC,8!2和8表示剪力PLLS,SMIS6,SMIS12,1,0!绘制弯矩图PLLS,SMIS1,SMIS7,1,0!绘制轴力图PLLS,SMIS2,SMIS8,1,0!绘制剪力图!右隧道土体开挖施工!开挖右隧道土体SAVE/SOLANTYPE,,REST,2,1,0ASEL,S,AREA,,7,8,1ASEL,A,AREA,,12!重新启动求解器!选择右隧道面ESLA,SEKILL,ALLASEL,S,AREA,,24,25,1ASEL,A,AREA,,30ESLA,SEKILL,ALL!选择右隧道土体单元!对选择的单元给予“死单元”!选择左隧道面，每次计算都要重新杀死单元!选择左隧道土体单元对选择的单元给予“死单元”!进行初期支护，即激活喷射混凝土衬砌梁单元ESEL,S,ELEM,,1549,1574,1EALIVE,ALLSAVEFINISH!施加节点荷载!节点力计算/POST1ALLSEL,ALLESEL,ALL!人工手动选择开挖土体周围的一圈单元!!人工手动选择开挖土体周围的一圈单元!进入求解器ESEL,S,PNFORCE,ALLFINISH/SOL!节点力施加F,89,FX, 2.82E+04F,91,FX, -7.48E+01

F,92,FX,3.12E+03F,93,FX,6.27E+03F,94,FX,9.31E+03F,95,FX,1.22E+04F,96,FX,1.49E+04F,97,FX,1.74E+04F,98,FX,1.97E+04F,99,FX,2.18E+04F,100,FX,2.37E+04F,101,FX,2.50E+04F,102,FX,2.47E+04F,127,FX,-2.83E+04F,129,FX,-2.47E+04F,130,FX,-2.51E+04F,131,FX,-2.37E+04F,132,FX,-2.19E+04F,133,FX,-1.98E+04F,134,FX,-1.75E+04F,135,FX,-1.51E+04F,136,FX,-1.24E+04F,137,FX,-9.46E+03F,138,FX,-6.42E+03F,139,FX,-3.27E+03F,608,FX,-1.39E+04F,609,FX,1.39E+04F,610,FX,-4.18E+03F,611,FX,-2.19E+03F,612,FX,-1.54E+03F,613,FX,-1.07E+03F,614,FX,-7.81E+02F,615,FX,-5.70E+02F,616,FX,-4.13E+02F,617,FX,-2.88E+02F,618,FX,-1.86E+02F,619,FX,-9.81E+01F,620,FX,-1.85E+01F,621,FX,5.62E+01F,622,FX,1.31E+02F,623,FX,2.11E+02F,624,FX,2.99E+02F,625,FX,4.02E+02F,626,FX,5.28E+02F,627,FX,6.86E+02F,628,FX,8.98E+02F,629,FX,1.19E+03F,630,FX,1.66E+03F,631,FX,2.31E+03F,632,FX,4.31E+03F,89,FY,-2.04E+04F,91,FY,-8.12E+04F,92,FY,-8.05E+04F,93,FY,-7.86E+04F,94,FY,-7.54E+04F,95,FY,-7.09E+04F,96,FY,-6.53E+04F,97,FY,-5.87E+04F,98,FY,-5.10E+04F,99,FY,-4.26E+04F,100,FY,-3.38E+04F,101,FY,-2.53E+04F,102,FY,-1.83E+04F,127,FY,-2.06E+04F,129,FY,-1.85E+04F,130,FY,-2.55E+04F,131,FY,-3.40E+04F,132,FY,-4.28E+04F,133,FY,-5.12E+04F,134,FY,-5.88E+04F,135,FY,-6.55E+04F,136,FY,-7.10E+04F,137,FY,-7.54E+04F,138,FY,-7.87E+04F,139,FY,-8.06E+04F,608,FY,2.85E+04F,609,FY,2.86E+04F,610,FY,6.60E+04F,611,FY,6.25E+04F,612,FY,6.16E+04F,613,FY,6.11E+04F,614,FY,6.10E+04F,615,FY,6.08E+04F,616,FY,6.08E+