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 !保存數據庫建立幾何模型創(chuàng)建關鍵點。K,1,0,0,,!創(chuàng)建計算區(qū)域輪廓線關鍵點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,, !創(chuàng)建隧道輪廓線關鍵點，上半圓圓心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 !保存數據創(chuàng)建線。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, !創(chuàng)建關鍵點，為畫其它直線作打算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 !畫計算區(qū)域線和分塊線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 !不顯示節(jié)點編號!*/PNUM,ELEM,0 !不顯示單元編號/REPLOT !重新繪制!*LPLOT !繪制線圖，帶編號，如圖6-8所示。SAVE !保存數據庫創(chuàng)建面。本次的計算區(qū)域為橫向130m，豎向為60m，即左右兩側計算邊界為4倍左右雙線隧道總跨度，下部邊界為2倍隧道總高度。Al,35,1,36, !采納線創(chuàng)建面，依次創(chuàng)建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!保存數據,創(chuàng)建網格模型設置線被劃分的單元數。LESIZE,_Y1,,,8,,,,,1lesize,1,,,8,,,,,1 !設置單元大小,即L1線劃分成8個單元TYPE,1 !設置將要創(chuàng)建單元的類型MAT,1 !設置將要創(chuàng)建單元的材料REAL,1 !設置將要創(chuàng)建單元的幾何常數/PNUM,LINE,1 !顯示線編號LPLOT !繪制線圖，帶編號和單元數。Save !根據“l(fā)esize,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 !劃分左隧道周邊區(qū)域面AMAP,18,33,12,8,34AMAP,19,34,8,11,51AMAP,20,51,11,56,32AMAP,21,54,57,13,53 !劃分右隧道周邊區(qū)域面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個四邊形單元創(chuàng)建模擬噴射混凝土的梁單元。TYPE,2 !設置將要創(chuàng)建單元的類型MAT,3 !設置將要創(chuàng)建單元的材料REAL,1 !設置將要創(chuàng)建單元的幾何常數E,2,11 !通過兩個節(jié)點創(chuàng)建梁單元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 !保存數據創(chuàng)建CRD開挖中模擬臨時支撐的梁單元，如圖6-12所示。TYPE,2 !設置將要創(chuàng)建單元的類型MAT,3 !設置將要創(chuàng)建單元的材料REAL,1 !設置將要創(chuàng)建單元的幾何常數E,1,18 !通過兩個節(jié)點創(chuàng)建梁單元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 !在選擇的節(jié)點上施加“Ux”約束Allsel !選擇全部內容NSEL,S,LOC,Y,-22.5,-22.1d,all,uy,0 !在選擇的節(jié)點上施加“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 !選擇線上的節(jié)點ESLN,R,1 !選擇節(jié)點上的單元ESEL,R,TYPE,,2 !選擇2類單元Ealive,all !對選擇的單元賜予“生屬性”Save節(jié)點力計算。FINISH/POST1 !進入后處理器Allsel !選擇全部內容ESEL,S,,, !人工手動選擇開挖土體四周的一圈單元圖6-16所示NFORCE,ALL !對選擇的節(jié)點求解節(jié)點力SaveFINISH/sol !進入求解理器節(jié)點力施加，加上節(jié)點力和位移邊界條件后的模型如圖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計算下一步開挖計算所施加的節(jié)點力，如圖6-22所示。ESEL,S,,, !人工手動選擇開挖土體四周的一圈單元NFORCE,ALL !對選擇的節(jié)點求解節(jié)點力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 !選擇線上的節(jié)點ESLN,R,1 !選擇節(jié)點上的單元ESEL,R,TYPE,,2 !選擇2類單元Ealive,all !對選擇的單元賜予“生屬性”Save節(jié)點力施加，加上節(jié)點力和位移邊界條件后的模型如圖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計算下一步開挖計算所施加的節(jié)點力，如圖6-28所示。ESEL,S,,, !人工手動選擇開挖土體四周的一圈單元。NFORCE,ALL !對選擇的節(jié)點求解節(jié)點力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!選擇線上的節(jié)點ESLN,R,1!選擇節(jié)點上的單元ESEL,R,TYPE,,2!選擇2類單元Ealive,all!對選擇的單元賜予“生屬性”Save節(jié)點力施加，加上節(jié)點力和位移邊界條件后的模型如圖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）計算下一步開挖計算所施加的節(jié)點力，如圖6-34所示。ESEL,S,,, !人工手動選擇開挖土體四周的一圈單元。NFORCE,ALL !對選擇的節(jié)點求解節(jié)點力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 !選擇線上的節(jié)點ESLN,R,1 !選擇節(jié)點上的單元ESEL,R,TYPE,,2 !選擇2類單元Ealive,all !對選擇的單元賜予“生屬性”Save（3）節(jié)點力施加，加上節(jié)點力和位移邊界條件后的模型如圖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 !創(chuàng)建梁單元內力表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）計算下一步開挖計算所施加的節(jié)點力，如圖6-48所示。ESEL,S,,, !人工手動選擇開挖土體四周的一圈單元。NFORCE,ALL !對選擇的節(jié)點求解節(jié)點力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 !選擇線上的節(jié)點ESLN,R,1 !選擇節(jié)點上的單元ESEL,R,TYPE,,2 !選擇2類單元Ealive,all !對選擇的單元賜予“生屬性”Save（3）節(jié)點力施加，加上節(jié)點力和位移邊界條件后的模型如圖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）計算下一步開挖計算所施加的節(jié)點力，如圖6-54所示。ESEL,S,,,!人工手動選擇開挖土體四周的一圈單元。NFORCE,ALL!對選擇的節(jié)點求解節(jié)點力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 !選擇線上的節(jié)點ESLN,R,1 !選擇節(jié)點上的單元ESEL,R,TYPE,,2 !選擇2類單元Ealive,all !對選擇的單元賜予“生屬性”Save（3）節(jié)點力施加，加上節(jié)點力和位移邊界條件后的模型如圖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）計算下一步開挖計算所施加的節(jié)點力，如圖6-60所示。ESEL,S,,,!人工手動選擇開挖土體四周的一圈單元。NFORCE,ALL!對選擇的節(jié)點求解節(jié)點力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 !選擇線上的節(jié)點ESLN,R,1 !選擇節(jié)點上的單元ESEL,R,TYPE,,2 !選擇2類單元Ealive,all !對選擇的單元賜予“生屬性”Save（3）節(jié)點力施加，加上節(jié)點力和位移邊界條件后的模型如圖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）計算下一步開挖計算所施加的節(jié)點力，如圖6-66所示。ESEL,S,,, !人工手動選擇開挖土體四周的一圈單元。NFORCE,ALL !對選擇的節(jié)點求解節(jié)點力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 !選擇線上的節(jié)點ESLN,R,1 !選擇節(jié)點上的單元ESEL,R,TYPE,,2 !選擇2類單元Ealive,all !對選擇的單元賜予“生屬性”Save（3）節(jié)點力施加，加上節(jié)點力和位移邊界條件后的模型如圖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 !創(chuàng)建梁單元內力表ETABLE,,SMISC,7ETABLE,,SMISC,2ETABLE,,SMISC,8ETABL