ansys在公路工程中的应用实例

/施工过程仿真分析材料、实常数和单元类型定义/TITLE，Mechanicalanalysisonrailwaytunnel1ndlining !确定分析标题/NOPR !菜单过滤设置/PMETH，OFF，0KEYW，PR_SET，1KEYW，PR_STRUC，1 !保留结构分析部分菜单/COM，/COM，PreferencesforGUIfilteringhavebeensettodisplay:/COM，Structural/PREP7 !进入前处理器ET,1,PLANE42 !设置实体单元类型KEYOPT,1,1,0KEYOPT,1,2,0KEYOPT,1,3,2 !设置为平面应变模式KEYOPT,1,5,0KEYOPT,1,6,0ET,2,PLANE42 !设置实体单元类型KEYOPT,1,1,0KEYOPT,1,2,0KEYOPT,1,3,2 !设置为平面应变模式KEYOPT,1,5,0KEYOPT,1,6,0ET，3，BEAM3 !设置梁单元类型R,1,0.25,1.5e-3,0.25,,,, !设置梁单元几何常数MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,EX,1,,0.8e9 !输入弹性模量（围岩）MPDATA,PRXY,1,,0.38 !输入泊松比（围岩）MPDATA,DENS,1,,2000 !输入密度（围岩）TB,DP,1,,, !采纳DP准则进行弹塑性分析TBMODIF,1,1,0.2e6 !输入凝合力（围岩）TBMODIF,1,2,30 !输入摩擦角（围岩）TBMODIF,1,3,MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,EX,2,,3.2e9 !输入弹性模量（加固圈）MPDATA,PRXY,2,,0.32 !输入泊松比（加固圈）MPDATA,DENS,2,,2200 !输入密度（加固圈）TB,DP,2,,, !采纳DP准则TBMODIF,1,1,0.5e6 !输入凝合力（加固圈）TBMODIF,1,2,35 !输入摩擦角（加固圈）TBMODIF,1,3,MPDATA,EX,3,,27.5e9 !输入弹性模量（喷射混凝土）MPDATA,PRXY,3,,0.2 !输入泊松比（喷射混凝土）MPDATA,DENS,3,,2500 !输入密度（喷射混凝土）SAVE !保存数据库建立几何模型创建关键点。K,1,0,0,,!创建计算区域轮廓线关键点K,2,65,-22.3,,K,3,-65,-22.3,,K,4,65,37.1,,K,5,19.6,37.1,,K,6,-65,-11.8,,K,7,-44.6,-22.3,,K,8,-44.6,0,,K,9,65,0,,K,10,15,-22.3,,K,11,-15,-22.3,,K,12,-15,17.1,,K,13,15,34.4,,K,14,-15,1.169,, !创建隧道轮廓线关键点，上半圆圆心K,15,-18.408,0.22,, !左下拱脚圆心K,16,-15.022,7.887,, !仰拱圆心K,17,0,25.742,,K,18,0,-22.3,,K,27,-11.636,0.22,, !右下拱脚圆心SAVE !保存数据创建线。Circle,14,5.4 !画上半圆Circle,15,1.9 !画左下拱脚圆Circle,27,1.9 !画右下拱脚圆LDELE,5,6,1 !删除左右下拱脚圆的部分线LDELE,8,11,1LOVLAP,7,3 !作叠加布尔运算LOVLAP,12,4LDELE,6,10,2LDELE,7LDELE,13,16,3NUMMRG,ALL,,,,LOW !合并全部元素NUMCMP,ALL !压缩全部编号FLST,2,3,4,ORDE,3 !作加运算FITEM,2,4FITEM,2,5FITEM,2,6LCOMB,P51X,,0FLST,2,3,4,ORDE,3 !作加运算FITEM,2,3FITEM,2,7FITEM,2,9LCOMB,P51X,,0NUMMRG,ALL,,,,LOW !合并全部元素NUMCMP,ALL !压缩全部编号circle,16,10.178 !画仰拱圆LOVLAP,3,7 !作叠加布尔运算LOVLAP,4,8LDELE,19 !删除部分线LDELE,5,7,1LDELE,10,11,1LDELE,14,16,2FLST,2,4,4,ORDE,4 !做线加布尔运算FITEM,2,9FITEM,2,12FITEM,2,-13FITEM,2,15LCOMB,P51X,,0FLST,2,4,4,ORDE,4 !做线加布尔运算FITEM,2,3FITEM,2,17FITEM,2,-18FITEM,2,20LCOMB,P51X,,0Circle,14,8.9 !画上半圆(加固圈)Circle,15,5.4 !画左下拱脚圆(加固圈)Circle,27,5.4 !画右下拱脚圆(加固圈)ldele,8,12,2ldele,13,15,1LOVLAP,6,11 !作叠加布尔运算LOVLAP,7,16ldele,10,14,2 !删除部分线ldele,17,19,2ldele,11circle,16,13.678 !画仰拱圆(加固圈)LOVLAP,11,13 !作叠加布尔运算LOVLAP,12,18ldele,10,11,1 !删除部分线ldele,23,26,3ldele,7,17,10ldele,14,22,8FLST,2,5,4,ORDE,5 !做线加布尔运算FITEM,2,8FITEM,2,15FITEM,2,-16FITEM,2,19FITEM,2,21LCOMB,P51X,,0FLST,2,5,4,ORDE,5 !做线加布尔运算FITEM,2,6FITEM,2,13FITEM,2,20FITEM,2,24FITEM,2,-25LCOMB,P51X,,0FLST,3,8,4,ORDE,4 !对称操作，由左隧道镜像成右隧道。FITEM,3,1FITEM,3,-6FITEM,3,8FITEM,3,-9LSYMM,X,P51X,,,,0,0NUMMRG,ALL,,,,LOW !合并全部元素NUMCMP,ALL !压缩全部编号SAVE !所画出的隧道轮廓图见图6-7所示。k,56,0,1.169,0, !创建关键点，为画其它直线作打算k,57,30.044,1.169,0,k,58,30.044,-22.3,0,k,59,30.044,37.1,0,k,60,15,1.169,0,LSTR,3,7 !画计算区域线和分块线LSTR,7,8LSTR,8,6LSTR,6,3LSTR,7,11LSTR,11,18LSTR,18,10LSTR,10,58LSTR,58,2LSTR,2,9LSTR,9,4LSTR,4,59LSTR,59,5LSTR,5,13LSTR,13,17LSTR,17,12LSTR,12,8LSTR,21,14LSTR,14,19LSTR,14,20LSTR,14,31LSTR,44,60LSTR,60,46LSTR,60,45LSTR,60,47LSTR,45,53LSTR,46,54LSTR,47,55LSTR,44,52LSTR,52,56LSTR,54,57LSTR,57,59LSTR,57,58LSTR,55,10LSTR,56,17LSTR,56,18LSTR,53,13LSTR,19,32LSTR,32,56LSTR,51,11LSTR,33,12LSTR,20,33LSTR,21,34LSTR,34,8LSTR,31,51LSTR,57,9/PNUM,KP,0 !不显示关键点编号/PNUM,LINE,1 !显示线编号/PNUM,AREA,0 !不显示面编号/PNUM,VOLU,0 !不显示体编号/PNUM,NODE,0 !不显示节点编号!*/PNUM,ELEM,0 !不显示单元编号/REPLOT !重新绘制!*LPLOT !绘制线图，带编号，如图6-8所示。SAVE !保存数据库创建面。本次的计算区域为横向130m，竖向为60m，即左右两侧计算边界为4倍左右双线隧道总跨度，下部边界为2倍隧道总高度。Al,35,1,36, !采纳线创建面，依次创建27个面Al,36,2,34,Al,34,9,37,Al,37,3,35, !左隧道内核心土4个面绘制完成Al,54,4,58,1,Al,58,5,59,2,Al,59,8,61,9,Al,61,6,54,3, !左隧道加固圈土4个面绘制完成Al,39,10,40,Al,40,7,38,Al,38,11,41,Al,41,16,39, !右隧道内核心土4个面绘制完成Al,43,13,42,10,Al,42,12,45,7,Al,45,14,44,11,Al,44,15,43,16,!右隧道加固圈土4个面绘制完成Al,55,51,32,57,4,Al,57,33,60,5,Al,60,18,21,56,8,Al,56,22,52,55,6,!左隧道外围土体4个面绘制完成Al,47,48,29,30,53,13,Al,53,31,51,46,12,Al,46,52,23,50,14,Al,50,24,49,47,15,!右隧道外围土体4个面绘制完成Al,18,19,20,17,!左侧边疆1个面绘制完成Al,49,25,26,62,Al,62,27,28,48,!右侧边疆2个面绘制完成/PNUM,KP,0/PNUM,LINE,0/PNUM,AREA,1APLOT!绘制带编号的面,生成的面如图6-9所示。SAVE!保存数据,创建网格模型设置线被划分的单元数。LESIZE,_Y1,,,8,,,,,1lesize,1,,,8,,,,,1 !设置单元大小,即L1线划分成8个单元TYPE,1 !设置将要创建单元的类型MAT,1 !设置将要创建单元的材料REAL,1 !设置将要创建单元的几何常数/PNUM,LINE,1 !显示线编号LPLOT !绘制线图，带编号和单元数。Save !根据“lesize,1,,,8,,,,,1”设置线单元数NL(NDIV)NL(NDIV)NL(NDIV)NL(NDIV)NL(NDIV)358398544554514368408584574324348388594604224378418614564524187843447418428108424534214981184544642043816844450419448128314301174581382642922446814825104854948815862102752810338234将全部面划分单元，其单元图如图6-11所示。AMAP,1,14,19,20, !划分面1，依次单击关键点14、19、20进行映射划分AMAP,2,14,20,21,AMAP,3,14,21,31,AMAP,4,14,31,19, !面1到4，为左隧道内开挖土体AMAP,9,60,46,45,AMAP,10,60,45,44,AMAP,11,60,44,47,AMAP,12,60,47,46, !面9到12，为右隧道内开挖土体AMAP,7,21,34,51,31AMAP,8,31,51,32,19 !面7和8，为左隧道仰拱下土体AMAP,15,44,52,55,47AMAP,16,47,55,54,46 !面15和16，为右隧道仰拱下土体TYPE,1 !转换为加固圈土体单元类型和材料MAT,2AMAP,5,19,32,33,20 !划分左隧道加固圈土体AMAP,6,20,33,34,21AMAP,13,46,54,53,45 !划分右隧道加固圈土体AMAP,14,45,53,52,44TYPE,1 !转换为一般土体单元类型和材料MAT,1AMAP,17,32,56,12,33 !划分左隧道周边区域面AMAP,18,33,12,8,34AMAP,19,34,8,11,51AMAP,20,51,11,56,32AMAP,21,54,57,13,53 !划分右隧道周边区域面AMAP,22,53,13,56,52AMAP,23,52,56,10,55AMAP,24,55,10,57,54AMAP,25,7,8,6,3 !划分左隧道边界面AMAP,26,58,2,9,57 !划分左隧道边界面AMAP,27,57,9,4,59SAVE !保存数据，本次共划分成1002个四边形单元创建模拟喷射混凝土的梁单元。TYPE,2 !设置将要创建单元的类型MAT,3 !设置将要创建单元的材料REAL,1 !设置将要创建单元的几何常数E,2,11 !通过两个节点创建梁单元E,11,12E,12,13E,13,14E,14,15E,15,16E,16,17E,17,10E,10,63E,63,64E,64,65E,65,66E,66,67E,67,68E,68,69E,69,62E,62,122E,122,123E,123,124E,124,125E,125,126E,126,127E,127,128E,128,114E,114,172E,172,171E,171,170E,170,169E,169,168E,168,167E,167,166E,166,2 !以上梁单元为左隧道喷射混凝土部分Save !保存数据E,211,226E,226,225E,225,224E,224,223E,223,222E,222,221E,221,220E,220,219E,219,278E,278,277E,277,276E,276,275E,275,274E,274,273E,273,272E,272,271E,271,331E,331,332E,332,333E,333,334E,334,335E,335,336E,336,337E,337,323E,323,318E,323,381E,381,380E,380,379E,379,378E,378,377E,377,376E,376,375E,375,211 !以上梁单元为右隧道喷射混凝土部分SAVE !保存数据创建CRD开挖中模拟临时支撑的梁单元，如图6-12所示。TYPE,2 !设置将要创建单元的类型MAT,3 !设置将要创建单元的材料REAL,1 !设置将要创建单元的几何常数E,1,18 !通过两个节点创建梁单元E,18,19E,19,20E,20,21E,21,22E,22,23E,23,24E,24,10E,1,76E,76,75E,75,74E,74,73E,73,72E,72,71E,71,70E,70,62E,1,115E,115,116E,116,117E,117,118E,118,119E,119,120E,120,121E,121,114E,1,3E,3,4E,4,5E,5,6E,6,7E,7,8E,8,9E,9,2 !以上梁单元为左隧道临时支撑喷射混凝土部分Save !保存数据E,210,227E,227,228E,228,229E,229,230E,230,231E,231,232E,232,233E,233,219E,210,285E,285,284E,284,283E,283,282E,282,281E,281,280E,280,279E,279,271E,210,324E,324,325E,325,326E,326,327E,327,328E,328,329E,329,330E,330,323E,210,212E,212,213E,213,214E,214,215E,215,216E,216,217E,217,218E,218,211Finish!返回上级书目SAVE!保存数据以上梁单元为右隧道临时支撑喷射混凝土部分，共128个梁单元。加载和自重应力场求解/SOL !进入求解器NSEL,S,LOC,X,-65.1,-64.9NSEL,A,LOC,X,64.9,65.1d,all,ux,0 !在选择的节点上施加“Ux”约束Allsel !选择全部内容NSEL,S,LOC,Y,-22.5,-22.1d,all,uy,0 !在选择的节点上施加“Uy”约束ACEL,0,10,0, !在Y方向施加重力加速度SAVE !施加了位移约束和自重应力场后见图6-13所示。NROPT,FULL,, !采纳全牛顿－拉普森法进行求解Allsel !选择全部内容ESEL,S,TYPE,,2 !选择2类单元Ekill,all !对选择的单元赐予“死属性”Allsel !选择全部内容Solve !求解计算Finish !求解结束返回MainMenu主菜单SAVE竖向位移（PLNSOL,U,Y,0,1）和应力（PLNSOL,S,Y,0,1）如图6-14和图6-15所示。先左洞后右洞开挖模拟分析左洞左上部分土体开挖施工模拟分析开挖土体。/post1Finish !求解结束返回MainMenu主菜单/SOL !进入求解器Allsel !选择全部内容ANTYPE,,REST,1,13,0 !重新启动求解器ASEL,S,,,2 !选择左洞左上部分土风光 ESLA,R !选择左洞左上部分土体单元Ekill,all !对选择的单元赐予“死属性”Save进行喷锚支护，即激活喷射混凝土衬砌梁单元和临时支撑单元。Allsel !选择全部内容ESEL,S,TYPE,,2 !选择2类单元Ekill,all !对选择的单元赐予“死属性”Allsel !选择全部内容lsel,s,,,2,,, !选择梁单元所在的线lsel,a,,,34,36,2,NSLL,R,1 !选择线上的节点ESLN,R,1 !选择节点上的单元ESEL,R,TYPE,,2 !选择2类单元Ealive,all !对选择的单元赐予“生属性”Save节点力计算。FINISH/POST1 !进入后处理器Allsel !选择全部内容ESEL,S,,, !人工手动选择开挖土体四周的一圈单元图6-16所示NFORCE,ALL !对选择的节点求解节点力SaveFINISH/sol !进入求解理器节点力施加，加上节点力和位移边界条件后的模型如图6-17所示。f, 1 ,fx, 4.37E-09 !施加X方向集中力f, 1 ,fx, -0.4743E+05 !施加Y方向集中力f, 1 ,fy, 0.9367E+05f, 10 ,fx, -0.4789E+05f, 10 ,fy, -0.1099E+06f, 18 ,fx, -0.1116E+06f, 18 ,fy, -0.2069E+05f, 19 ,fx, -0.1041E+06f, 19 ,fy, -0.2041E+05f, 20 ,fx, -0.9727E+05f, 20 ,fy, -0.2031E+05f, 21 ,fx, -0.9131E+05f, 21 ,fy, -0.2034E+05f, 22 ,fx, -0.8623E+05f, 22 ,fy,-0.2040E+05f, 23 ,fx, -0.8191E+05f, 23 ,fy, -0.2042E+05f, 24 ,fx, -0.7817E+05f, 24 ,fy, -0.2032E+05f, 62 ,fx, 0.7342E+05f, 62 ,fy, 0.6471E+05f, 63 ,fx, -5701f, 63 ,fy, -0.1873E+06f, 64 ,fx, 0.2051E+05f, 64 ,fy, -0.1652E+06f, 65 ,fx, 0.4760E+05f, 65 ,fy, -0.1410E+06f, 66 ,fx, 0.7627E+05f, 66 ,fy, -0.1149E+06f, 67 ,fx, 0.1060E+06f, 67 ,fy, -0.8668E+05f, 68 ,fx, 0.1340E+06f, 68 ,fy, -0.5737E+05f, 69 ,fx, 0.1507E+06f, 69 ,fy, -0.2829E+05f, 70 ,fx, 0.1678E+05f, 70 ,fy, 0.1740E+06f, 71 ,fx, 0.1940E+05f, 71 ,fy, 0.1803E+06f, 72 ,fx, 0.2109E+05f, 72 ,fy, 0.1853E+06f, 73 ,fx, 0.2143E+05f, 73 ,fy, 0.1903E+06f, 74 ,fx, 0.2161E+05f, 74 ,fy, 0.1953E+06f, 75 ,fx, 0.2150E+05f, 75 ,fy, 0.2002E+06f, 76 ,fx, 0.2133E+05f, 76 ,fy, 0.2051E+06求解计算。Allsel !选择全部内容Solve !求解计算Finish !求解结束返回MainMenu主菜单SAVE后处理，两个方向的应力和位移，如图6-18～图6-21所示。/post1PLNSOL,U,X,0,1 !X方向位移PLNSOL,U,Y,0,1 !Y方向位移PLNSOL,S,X,0,1 !X方向应力PLNSOL,S,Y,0,1 !Y方向应力Finish !求解结束返回MainMenu主菜单SAVE计算下一步开挖计算所施加的节点力，如图6-22所示。ESEL,S,,, !人工手动选择开挖土体四周的一圈单元NFORCE,ALL !对选择的节点求解节点力FINISH !求解结束返回MainMenu主菜单Save左洞左下部分土体开挖施工模拟分析开挖土体。/SOL !进入求解器Allsel !选择全部内容ANTYPE,,REST,2,13,0 !重新启动求解器ASEL,S,,,2,3,1 !选择左洞左上和左下部分土风光 ESLA,R !选择左洞左上和左下部分土体单元Ekill,all !对选择的单元赐予“死属性”Save进行喷锚支护，即激活喷射混凝土衬砌梁单元和临时支撑单元。Allsel !选择全部内容ESEL,S,TYPE,,2 !选择2类单元Ekill,all !对选择的单元赐予“死属性”Allsel !选择全部内容lsel,s,,,2,9,7, !选择梁单元所在的线lsel,a,,,36,37,1,NSLL,R,1 !选择线上的节点ESLN,R,1 !选择节点上的单元ESEL,R,TYPE,,2 !选择2类单元Ealive,all !对选择的单元赐予“生属性”Save节点力施加，加上节点力和位移边界条件后的模型如图6-23所示。f, 1 ,fx, -0.7585E+05 !施加X方向集中力f, 1 ,fy, 0.5203E+05 !施加Y方向集中力f, 62 ,fx, 0.7840E+05f, 62 ,fy, 219.8f, 114 ,fx, 0.1589E+05f, 114 ,fy, 0.1046E+06f, 115 ,fx, -0.8960E+05f, 115 ,fy, -0.5690E+05f, 116 ,fx, -0.5327E+05f, 116 ,fy, -0.6903E+05f, 117 ,fx, -0.5832E+05f, 117 ,fy, -0.7317E+05f, 118 ,fx, -0.6716E+05f, 118 ,fy, -0.7179E+05f, 119 ,fx, -0.7668E+05f, 119 ,fy, -0.6813E+05f, 120 ,fx, -0.8512E+05f, 120 ,fy, -0.6344E+05f, 121 ,fx, -0.9239E+05f, 121 ,fy, -0.5821E+05f, 122 ,fx, 0.8784E+05f, 122 ,fy, -0.7074E+05f, 123 ,fx, 0.3195E+05f, 123 ,fy, 4682f, 124 ,fx, 0.1113E+05f, 124 ,fy, 0.5868E+05f, 125 ,fx, 0.3296E+05f, 125 ,fy, 0.6390E+05f, 126 ,fx, 0.8549E+05f, 126 ,fy, 0.8524E+05f, 127 ,fx, 0.1220E+06f, 127 ,fy, 0.1426E+06f, 128 ,fx, 0.1310E+06f, 128 ,fy, 0.2152E+06求解计算。Allsel !选择全部内容Solve !求解计算Finish !求解结束返回MainMenu主菜单SAVE后处理，两个方向的应力和位移，如图6-24～图6-27所示。/post1PLNSOL,U,X,0,1!X方向位移PLNSOL,U,Y,0,1!Y方向位移PLNSOL,S,X,0,1!X方向应力PLNSOL,S,Y,0,1!Y方向应力Finish!求解结束返回MainMenu主菜单SAVE计算下一步开挖计算所施加的节点力，如图6-28所示。ESEL,S,,, !人工手动选择开挖土体四周的一圈单元。NFORCE,ALL !对选择的节点求解节点力FINISHSave左洞右上部分土体开挖施工模拟分析开挖土体。/SOL !进入求解器Allsel !选择全部内容ANTYPE,,REST,3,19,0 !重新启动求解器ASEL,S,,,1,3,1 !选择左洞左上、左下和右上部分土风光 ESLA,R !选择左洞左上、左下和右上部分土体单元Ekill,all !对选择的单元赐予“死属性”Save进行喷锚支护，即激活喷射混凝土衬砌梁单元和临时支撑单元。Allsel !选择全部内容ESEL,S,TYPE,,2 !选择2类单元Ekill,all !对选择的单元赐予“死属性”Allsel !选择全部内容lsel,s,,,1,2,1, !选择梁单元所在的线lsel,a,,,35,37,2,lsel,a,,,9,,,NSLL,R,1!选择线上的节点ESLN,R,1!选择节点上的单元ESEL,R,TYPE,,2!选择2类单元Ealive,all!对选择的单元赐予“生属性”Save节点力施加，加上节点力和位移边界条件后的模型如图6-29所示。f, 1 ,fx, 8888!施加X方向集中力f, 1 ,fy, 0.1627E+05!施加Y方向集中力f, 2 ,fx, -0.3130E+05f, 2 ,fy, 0.4226E+05f, 3 ,fx, -4424f, 3 ,fy, 0.2071E+05f, 4 ,fx, 0.1028E+05f, 4 ,fy, 0.2229E+05f, 5 ,fx, 0.2226E+05f, 5 ,fy, 0.2401E+05f, 6 ,fx, 0.2950E+05f, 6 ,fy, 0.2507E+05f, 7 ,fx, 0.3382E+05f, 7 ,fy, 0.2565E+05f, 8 ,fx, 0.3839E+05f, 8 ,fy, 0.2614E+05f, 9 ,fx, 0.4493E+05f, 9 ,fy, 0.2697E+05f, 10 ,fx, -0.1884E+05f, 10 ,fy, -0.2575E+06f, 11 ,fx, -5824f, 11 ,fy, -0.1401E+06f, 12 ,fx, -0.1110E+05f, 12 ,fy, -0.1420E+06f, 13 ,fx, -0.2394E+05f, 13 ,fy, -0.1375E+06f, 14 ,fx, -0.3887E+05f, 14 ,fy, -0.1272E+06f, 15 ,fx, -0.5173E+05f, 15 ,fy, -0.1123E+06f, 16 ,fx, -0.6407E+05f, 16 ,fy, -0.9350E+05f, 17 ,fx, -0.6281E+05f, 17 ,fy, -0.1098E+06求解计算。Allsel !选择全部内容Solve !求解计算Finish !求解结束返回MainMenu主菜单SAVE后处理，两个方向的应力和位移，如图6-30～图6-33所示。/post1PLNSOL,U,X,0,1!X方向位移PLNSOL,U,Y,0,1!Y方向位移PLNSOL,S,X,0,1!X方向应力PLNSOL,S,Y,0,1!Y方向应力Finish !求解结束返回MainMenu主菜单SAVE（6）计算下一步开挖计算所施加的节点力，如图6-34所示。ESEL,S,,, !人工手动选择开挖土体四周的一圈单元。NFORCE,ALL !对选择的节点求解节点力FINISHSave左洞右下部分土体开挖施工模拟分析（1）开挖土体。/SOL !进入求解器Allsel !选择全部内容ANTYPE,,REST,4,13,0 !重新启动求解器ASEL,S,,,1,4,1 !选择左洞全部核心土风光 ESLA,R !选择左洞全部核心土体单元Ekill,all !对选择的单元赐予“死属性”Save提示：在开挖左隧道右下部分土体时，同样要进行第一部分（左上）、其次部分（左下）和第三部分（右上）上全部单元的杀死、梁单元的激活操作，这是重复的，但重新启动后必需进行。（2）进行喷锚支护，即激活喷射混凝土衬砌梁单元和临时支撑单元。Allsel !选择全部内容ESEL,S,TYPE,,2 !选择2类单元Ekill,all !对选择的单元赐予“死属性”Allsel !选择全部内容lsel,s,,,1,2,1, !选择梁单元所在的线lsel,a,,,3,9,6,NSLL,R,1 !选择线上的节点ESLN,R,1 !选择节点上的单元ESEL,R,TYPE,,2 !选择2类单元Ealive,all !对选择的单元赐予“生属性”Save（3）节点力施加，加上节点力和位移边界条件后的模型如图6-35所示。f, 2 ,fx, -0.3948E+05 !施加X方向集中力f, 2 ,fy, -0.1837E+05 !施加Y方向集中力f, 114 ,fx, -0.5205E+05f, 114 ,fy, 0.2791E+06f, 166 ,fx, -0.7718E+05f, 166 ,fy, 0.6841E+05f, 167 ,fx, -0.3975E+05f, 167 ,fy, 0.6883E+05f, 168 ,fx, -0.1029E+06f, 168 ,fy, 0.1820E+06f, 169 ,fx, -0.1153E+06f, 169 ,fy, 0.1805E+06f, 170 ,fx, -0.9114E+05f, 170 ,fy, 0.1538E+06f, 171 ,fx, -0.7583E+05f, 171 ,fy, 0.1488E+06f, 172 ,fx, -0.7772E+05f, 172 ,fy, 0.1858E+06（4）求解计算。Allsel !选择全部内容Solve !求解计算Finish !求解结束返回MainMenu主菜单SAVE（5）后处理，两个方向的应力和位移，如图6-36～图6-42所示。内力和变形如图6-43～图6-47所示。/post1PLDISP,1 !绘制变形图PLNSOL,U,X,0,1 !X方向位移PLNSOL,U,Y,0,1 !Y方向位移PLNSOL,S,X,0,1 !X方向应力PLNSOL,S,Y,0,1 !Y方向应力PLNSOL,S,1,0,1 !第一主应力PLNSOL,S,3,0,1 !第三主应力PLNSOL,EPPL,EQV,0,1 !塑性应变ETABLE,,SMISC,1 !创建梁单元内力表ETABLE,,SMISC,7ETABLE,,SMISC,2ETABLE,,SMISC,8ETABLE,,SMISC,6ETABLE,,SMISC,12ESEL,S,LIVE !选择“属性为活”的单元ESEL,R,TYPE,,2 !选择单元类型为2的全部单元/REPLOT !从新绘制图形PLLS,SMIS1,SMIS7,0.2,0 !画轴力图，显示比例为0.2PLLS,SMIS2,SMIS8,0.2,0 !画剪力图，显示比例为0.2PLLS,SMIS6,SMIS12,-0.2,0 !画弯矩图，显示比例为-0.2SET,LIST !设置计算步骤LCDEF,1,1,13,LCDEF,5,5,13,LCASE,5,LCOPER,SUB,1,,, !将第四步的结果减去第一部（自重应力场）的结果PLDISP,1 !画变形喷射混凝土的图Finish !求解结束返回MainMenu主菜单SAVE（6）计算下一步开挖计算所施加的节点力，如图6-48所示。ESEL,S,,, !人工手动选择开挖土体四周的一圈单元。NFORCE,ALL !对选择的节点求解节点力FINISHSave右洞左上部分土体开挖施工模拟分析（1）开挖土体。/post1Finish !求解结束返回MainMenu主菜单/SOL !进入求解器Allsel !选择全部内容ANTYPE,,REST,5,13,0 !重新启动求解器ASEL,S,,,1,4,1 !选择左洞全部土风光 ASEL,a,,,10 !选择右洞左上部分土风光 ESLA,R !选择左洞全部和右洞左上部分土体单元Ekill,all !对选择的单元赐予“死属性”Save（2）进行喷锚支护，即激活喷射混凝土衬砌梁单元和临时支撑单元。Allsel !选择全部内容ESEL,S,TYPE,,2 !选择2类单元Ekill,all !对选择的单元赐予“死属性”Allsel !选择全部内容lsel,s,,,1,2,1, !选择左隧道梁单元所在的线lsel,a,,,3,9,6,lsel,a,,,7 !选择右隧道左上部分梁单元所在的线lsel,a,,,38,40,2,NSLL,R,1 !选择线上的节点ESLN,R,1 !选择节点上的单元ESEL,R,TYPE,,2 !选择2类单元Ealive,all !对选择的单元赐予“生属性”Save（3）节点力施加，加上节点力和位移边界条件后的模型如图6-49所示。f, 210 ,fy, 0.1999E+06 !施加X方向集中力f, 210 ,fx, -0.5649E+05 !施加Y方向集中力f, 210 ,fy, 0.1999E+05f, 219 ,fx, -0.4300E+05f, 219 ,fy, -0.2487E+05f, 227 ,fx, -0.1101E+06f, 227 ,fy,-1504f, 228 ,fx, -0.1014E+06f, 228 ,fy, -1912f, 229 ,fx, -0.9480E+05f, 229 ,fy, -2394f, 230 ,fx, -0.9065E+05f, 230 ,fy,-2903f, 231 ,fx, -0.8893E+05f, 231 ,fy, -3380f, 232 ,fx, -0.8921E+05f, 232 ,fy, -3784f, 233 ,fx, -0.9096E+05f, 233 ,fy,-4110f, 271 ,fx, 0.3675E+05f, 271 ,fy, 0.9479E+04f, 272 ,fx, 0.1417E+06f, 272 ,fy, -0.1460E+06f, 273 ,fx, 0.1542E+06f, 273 ,fy,-0.1087E+06f, 274 ,fx, 0.1408E+06f, 274 ,fy, -0.1748E+06f, 275 ,fx, 0.1171E+06f, 275 ,fy, -0.1401E+06f, 276 ,fx, 0.8841E+05f, 276 ,fy,-0.1987E+06f, 277 ,fx, 0.5735E+05f, 277 ,fy, -0.1462E+06f, 278 ,fx, 0.2457E+05f, 278 ,fy, -0.1795E+06f, 279 ,fx, 8959f, 279 ,fy,0.3789E+05f, 280 ,fx, 1788f, 280 ,fy, 0.3819E+05f, 281 ,fx, 2381f, 281 ,fy, 0.3851E+05f, 282 ,fx, 3068f, 282 ,fy,0.3888E+05f, 283 ,fx, 3518f, 283 ,fy, 0.3924E+05f, 284 ,fx, 3373f, 284 ,fy,0.3959E+05f, 285 ,fx, 3116f, 285 ,fy, 0.3991E+05（4）求解计算。Allsel !选择全部内容Solve !求解计算Finish !求解结束返回MainMenu主菜单SAVE（5）后处理，两个方向的应力和位移，如图6-50～图6-53所示。/post1PLNSOL,U,X,0,1 !X方向位移PLNSOL,U,Y,0,1!Y方向位移PLNSOL,S,X,0,1 !X方向应力PLNSOL,S,Y,0,1 !Y方向应力Finish !求解结束返回MainMenu主菜单SAVE（6）计算下一步开挖计算所施加的节点力，如图6-54所示。ESEL,S,,,!人工手动选择开挖土体四周的一圈单元。NFORCE,ALL!对选择的节点求解节点力FINISH!求解结束返回MainMenu主菜单Save右洞左下部分土体开挖施工模拟分析（1）开挖土体。/SOL !进入求解器Allsel !选择全部内容ANTYPE,,REST,6,14,0 !重新启动求解器ASEL,S,,,1,4,1 !选择左洞全部土风光 ASEL,a,,,10,11,1 !选择右洞左上和左下部分土风光 ESLA,R !选择左洞全部和右洞左上和左下部分土体单元Ekill,all !对选择的单元赐予“死属性”Save（2）进行喷锚支护，即激活喷射混凝土衬砌梁单元和临时支撑单元。Allsel !选择全部内容ESEL,S,TYPE,,2 !选择2类单元Ekill,all !对选择的单元赐予“死属性”Allsel !选择全部内容lsel,s,,,1,2,1, !选择左隧道梁单元所在的线lsel,a,,,3,9,6,lsel,a,,,7,11,4, !选择右隧道左上和左下部分梁单元所在的线lsel,a,,,40,41,1,NSLL,R,1 !选择线上的节点ESLN,R,1 !选择节点上的单元ESEL,R,TYPE,,2 !选择2类单元Ealive,all !对选择的单元赐予“生属性”Save（3）节点力施加，加上节点力和位移边界条件后的模型如图6-55所示。f, 210 ,fx, -0.1215E+05 !施加X方向集中力f, 210 ,fy, 0.1644E+05 !施加Y方向集中力f, 271 ,fx, 0.3771E+05f, 271 ,fy, 0.6092E+05f, 323 ,fx, -0.1915E+05f, 323 ,fy,-0.2613E+05f, 324 ,fx, -0.5155E+05f, 324 ,fy, -0.6288E+05f, 325 ,fx, -0.2299E+05f, 325 ,fy, -0.7694E+05f, 326 ,fx, -0.3873E+05f, 326 ,fy,-0.7586E+05f, 327 ,fx, -0.5527E+05f, 327 ,fy, -0.7253E+05f, 328 ,fx, -0.7057E+05f, 328 ,fy, -0.6596E+05f, 329 ,fx, -0.8310E+05f, 329 ,fy,-0.5898E+05f, 330 ,fx, -0.9315E+05f, 330 ,fy, -0.5150E+05f, 331 ,fx, 0.1309E+05f, 331 ,fy, 0.1305E+06f, 332 ,fx, 0.2309E+05f, 332 ,fy,0.1460E+06f, 333 ,fx, 0.4951E+05f, 333 ,fy, 0.1202E+06f, 334 ,fx, 0.2678E+05f, 334 ,fy, 0.1517E+06f, 335 ,fx, 0.4138E+05f, 335 ,fy,0.1684E+06f, 336 ,fx, 0.9038E+05f, 336 ,fy, 0.2027E+06f, 337 ,fx, 0.1065E+05f, 337 ,fy, 0.2152E+06（4）求解计算。Allsel !选择全部内容Solve !求解计算Finish !求解结束返回MainMenu主菜单SAVE（5）后处理，两个方向的应力和位移，如图6-56～图6-59所示。/post1PLNSOL,U,X,0,1 !X方向位移PLNSOL,U,Y,0,1 !Y方向位移PLNSOL,S,X,0,1 !X方向应力PLNSOL,S,Y,0,1 !Y方向应力Finish !求解结束返回MainMenu主菜单SAVE（6）计算下一步开挖计算所施加的节点力，如图6-60所示。ESEL,S,,,!人工手动选择开挖土体四周的一圈单元。NFORCE,ALL!对选择的节点求解节点力FINISHSave右洞右上部分土体开挖施工模拟分析（1）开挖土体。/SOL !进入求解器Allsel !选择全部内容ANTYPE,,REST,7,23,0 !重新启动求解器ASEL,S,,,1,4,1 !选择左洞全部土风光 ASEL,a,,,9,11,1 !选择右洞左上、左下和右上部分土风光 ESLA,R !选择所选面上部分土体单元Ekill,all !对选择的单元赐予“死属性”Save（2）进行喷锚支护，即激活喷射混凝土衬砌梁单元和临时支撑单元。Allsel !选择全部内容ESEL,S,TYPE,,2 !选择2类单元Ekill,all !对选择的单元赐予“死属性”Allsel !选择全部内容lsel,s,,,1,2,1, !选择左隧道梁单元所在的线lsel,a,,,3,9,6,lsel,a,,,7,11,4, !选择右隧道左上、左下和右上部分梁单元所在的线lsel,a,,,39,41,2,lsel,a,,,10,NSLL,R,1 !选择线上的节点ESLN,R,1 !选择节点上的单元ESEL,R,TYPE,,2 !选择2类单元Ealive,all !对选择的单元赐予“生属性”Save（3）节点力施加，加上节点力和位移边界条件后的模型如图6-61所示。f, 210 ,fx, -0.1116E+05 !施加X方向集中力f, 210 ,fy, 0.1416E+06 !施加Y方向集中力f, 211 ,fx, -0.5826E+05f, 211 ,fy, 0.1030E+06f, 212 ,fx, -0.3155E+05f, 212 ,fy,0.2814E+06f, 213 ,fx, -0.1581E+05f, 213 ,fy, 0.2795E+06f, 214 ,fx, -4423f, 214 ,fy, 0.2786E+06f, 215 ,fx, -43.02f, 215 ,fy,0.2793E+06f, 216 ,fx, -3788f, 216 ,fy, 0.2765E+06f, 217 ,fx, -9608f, 217 ,fy, 0.2747E+06f, 218 ,fx, -0.1991E+05f, 218 ,fy,0.2742E+06f, 219 ,fx, 0.4335E+05f, 219 ,fy, -0.6617E+06f, 220 ,fx, -0.9339E+05f, 220 ,fy, -0.4809E+06f, 221 ,fx, -0.8952E+05f, 221 ,fy,-0.3837E+06f, 222 ,fx, -0.6087E+05f, 222 ,fy, -0.3509E+06f, 223 ,fx, -0.3349E+05f, 223 ,fy, -0.3138E+06f, 224 ,fx, -0.1001E+05f, 224 ,fy,-0.2732E+06f, 225 ,fx, -0.2350E+05f, 225 ,fy, -0.2306E+06f, 226 ,fx, -0.2150E+05f, 226 ,fy, -0.1873E+06（4）求解计算。Allsel !选择全部内容Solve !求解计算Finish !求解结束返回MainMenu主菜单SAVE（5）后处理，两个方向的应力和位移，如图6-62～图6-65所示。/post1PLNSOL,U,X,0,1 !X方向位移PLNSOL,U,Y,0,1 !Y方向位移PLNSOL,S,X,0,1 !X方向应力PLNSOL,S,Y,0,1 !Y方向应力Finish !求解结束返回MainMenu主菜单SAVE（6）计算下一步开挖计算所施加的节点力，如图6-66所示。ESEL,S,,, !人工手动选择开挖土体四周的一圈单元。NFORCE,ALL !对选择的节点求解节点力FINISHSave右洞右下部分土体开挖施工模拟分析（1）开挖土体。/SOL !进入求解器Allsel !选择全部内容ANTYPE,,REST,8,14,0 !重新启动求解器ASEL,S,,,1,4,1 !选择左洞全部土风光 ASEL,a,,,9,12,1 !选择右洞全部土风光 ESLA,R !选择所选面上部分土体单元Ekill,all !对选择的单元赐予“死属性”Save（2）进行喷锚支护，即激活喷射混凝土衬砌梁单元和临时支撑单元。Allsel !选择全部内容ESEL,S,TYPE,,2 !选择2类单元Ekill,all !对选择的单元赐予“死属性”Allsel !选择全部内容lsel,s,,,1,2,1, !选择左隧道梁单元所在的线lsel,a,,,3,9,6,lsel,a,,,7,11,4, !选择右隧道梁单元所在的线lsel,a,,,10,16,6,NSLL,R,1 !选择线上的节点ESLN,R,1 !选择节点上的单元ESEL,R,TYPE,,2 !选择2类单元Ealive,all !对选择的单元赐予“生属性”Save（3）节点力施加，加上节点力和位移边界条件后的模型如图6-67所示。f, 211 ,fx, -0.4418E+05 !施加X方向集中力f, 211 ,fy, 0.6747E+05 !施加Y方向集中力f, 211 ,fx, -0.4418E+05f, 211 ,fy, 0.6747E+05f, 323 ,fx, -0.1273E+05f, 323 ,fy, 0.1440E+06f, 375 ,fx, -0.1353E+05f, 375 ,fy,0.6025E+05f, 376 ,fx, -0.4306E+05f, 376 ,fy, 0.1044E+06f, 377 ,fx, -0.1824E+05f, 377 ,fy, 0.1097E+06f, 378 ,fx, -0.1861E+05f, 378 ,fy,0.1033E+06f, 379 ,fx, -0.1369E+05f, 379 ,fy, 0.1133E+06f, 380 ,fx, -0.5347E+05f, 380 ,fy, 0.1299E+06f, 381 ,fx, -0.5852E+05f, 381 ,fy,0.1551E+06（4）求解计算。Allsel !选择全部内容Solve !求解计算Finish !求解结束返回MainMenu主菜单SAVE（5）后处理，两个方向的应力和位移，如图6-68～图6-79所示。/post1PLDISP,1 !绘制变形图PLNSOL,U,X,0,1 !X方向位移PLNSOL,U,Y,0,1 !Y方向位移PLNSOL,S,X,0,1 !X方向应力PLNSOL,S,Y,0,1 !Y方向应力PLNSOL,S,1,0,1 !第一主应力PLNSOL,S,3,0,1 !第三主应力PLNSOL,S,EQV,0,1 !等效应力PLNSOL,EPPL,EQV,0,1 !塑性应变ETABLE,,SMISC,1 !创建梁单元内力表ETABLE,,SMISC,7ETABLE,,SMISC,2ETABLE,,SMISC,8ETABL