第6章 高速公路分離式偏壓隧道施工過程仿真分析

1、第1章 大型有限元軟件ANSYS簡介596.1 施工過程仿真分析1 材料、實(shí)常數(shù)和單元類型定義/TITLE，Mechanical analysis on railway tunnel 1nd lining! 確定分析標(biāo)題/NOPR !菜單過濾設(shè)置/PMETH，OFF，0KEYW，PR_SET，1 KEYW，PR_STRUC，1 !保留結(jié)構(gòu)分析部分菜單/COM， /COM，Preferences for GUI filtering have been set to display:/COM， Structural /PREP7 !進(jìn)入前處理器ET,1,PLANE42 !設(shè)置實(shí)體單元類型KEYOP

2、T,1,1,0KEYOPT,1,2,0KEYOPT,1,3,2 !設(shè)置為平面應(yīng)變模式KEYOPT,1,5,0KEYOPT,1,6,0ET,2,PLANE42 !設(shè)置實(shí)體單元類型KEYOPT,1,1,0KEYOPT,1,2,0KEYOPT,1,3,2 !設(shè)置為平面應(yīng)變模式KEYOPT,1,5,0KEYOPT,1,6,0ET，3，BEAM3 !設(shè)置梁單元類型R,1,0.25,1.5e-3,0.25, , , , !設(shè)置梁單元幾何常數(shù)MPTEMP, MPTEMP,1,0MPDATA,EX,1,0.8e9 !輸入彈性模量（圍巖）MPDATA,PRXY,1,0.38 !輸入泊松比（圍巖）MPDATA,D

3、ENS,1,2000 !輸入密度（圍巖）TB,DP,1, !采用DP準(zhǔn)則進(jìn)行彈塑性分析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準(zhǔn)則TBMODIF,1,1,0.5e6 !輸入凝聚力（加固圈）TBMODIF,1,2,35 !輸入摩擦角（加固圈）TBMODIF,1,3, MP

4、DATA,EX,3,27.5e9 !輸入彈性模量（噴射混凝土）MPDATA,PRXY,3,0.2 !輸入泊松比（噴射混凝土）MPDATA,DENS,3,2500 !輸入密度（噴射混凝土）SAVE !保存數(shù)據(jù)庫2 建立幾何模型創(chuàng)建關(guān)鍵點(diǎn)。K,1,0,0, !創(chuàng)建計算區(qū)域輪廓線關(guān)鍵點(diǎn)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

5、,K,14,-15,1.169, !創(chuàng)建隧道輪廓線關(guān)鍵點(diǎn)，上半圓圓心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 !保存數(shù)據(jù)創(chuàng)建線。Circle,14,5.4 !畫上半圓Circle,15,1.9 !畫左下拱腳圓Circle,27,1.9 !畫右下拱腳圓LDELE,5,6,1 !刪除左右下拱腳圓的部分線LDELE,8,11,1 LOVLAP,7,3 !作疊加布爾運(yùn)算LOVLAP,12,4LDELE,6,10,2LDELE

6、,7LDELE,13,16,3NUMMRG,ALL, , , ,LOW !合并所有元素NUMCMP,ALL !壓縮所有編號FLST,2,3,4,ORDE,3 !作加運(yùn)算FITEM,2,4 FITEM,2,5 FITEM,2,6 LCOMB,P51X, ,0 FLST,2,3,4,ORDE,3 !作加運(yùn)算FITEM,2,3 FITEM,2,7 FITEM,2,9 LCOMB,P51X, ,0NUMMRG,ALL, , , ,LOW !合并所有元素NUMCMP,ALL !壓縮所有編號circle, 16,10.178 !畫仰拱圓LOVLAP,3,7 !作疊加布爾運(yùn)算LOVLAP,4,8LDELE,

7、19 !刪除部分線LDELE,5,7,1LDELE,10,11,1LDELE,14,16,2FLST,2,4,4,ORDE,4 !做線加布爾運(yùn)算FITEM,2,9 FITEM,2,12 FITEM,2,-13 FITEM,2,15 LCOMB,P51X, ,0 FLST,2,4,4,ORDE,4 !做線加布爾運(yùn)算FITEM,2,3 FITEM,2,17 FITEM,2,-18 FITEM,2,20 LCOMB,P51X, ,0Circle,14,8.9 !畫上半圓(加固圈)Circle,15,5.4 !畫左下拱腳圓(加固圈)Circle,27,5.4 !畫右下拱腳圓(加固圈)ldele,8,1

8、2,2ldele,13,15,1LOVLAP,6,11 !作疊加布爾運(yùn)算LOVLAP,7,16ldele,10,14,2 !刪除部分線ldele,17,19,2ldele,11circle, 16,13.678 !畫仰拱圓(加固圈)LOVLAP,11,13 !作疊加布爾運(yùn)算LOVLAP,12,18ldele,10,11,1 !刪除部分線ldele,23,26,3ldele,7,17,10ldele,14,22,8FLST,2,5,4,ORDE,5 !做線加布爾運(yùn)算FITEM,2,8 FITEM,2,15 FITEM,2,-16 FITEM,2,19 FITEM,2,21 LCOMB,P51X,

9、 ,0 FLST,2,5,4,ORDE,5 !做線加布爾運(yùn)算FITEM,2,6 FITEM,2,13 FITEM,2,20 FITEM,2,24 FITEM,2,-25 LCOMB,P51X, ,0 FLST,3,8,4,ORDE,4 !對稱操作，由左隧道鏡像成右隧道。FITEM,3,1 FITEM,3,-6 FITEM,3,8 FITEM,3,-9 LSYMM,X,P51X, , , ,0,0 NUMMRG,ALL, , , ,LOW !合并所有元素NUMCMP,ALL !壓縮所有編號SAVE !所畫出的隧道輪廓圖見圖6-7所示。k,56,0,1.169,0, !創(chuàng)建關(guān)鍵點(diǎn)，為畫其它直線作準(zhǔn)

10、備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, 8 LSTR, 8, 6 LSTR, 6, 3 LSTR, 7, 11 LSTR, 11, 18 LSTR, 18, 10 LSTR, 10, 58 LSTR, 58, 2 LSTR, 2, 9 LSTR, 9, 4 LSTR, 4, 59 LSTR, 59, 5 LSTR, 5, 13 LSTR, 13, 17 LSTR, 17, 12 LSTR, 12, 8 LSTR, 21,

11、14 LSTR, 14, 19 LSTR, 14, 20 LSTR, 14, 31 LSTR, 44, 60 LSTR, 60, 46 LSTR, 60, 45 LSTR, 60, 47 LSTR, 45, 53 LSTR, 46, 54 LSTR, 47, 55 LSTR, 44, 52 LSTR, 52, 56 LSTR, 54, 57 LSTR, 57, 59 LSTR, 57, 58 LSTR, 55, 10 LSTR, 56, 17 LSTR, 56, 18 LSTR, 53, 13 LSTR, 19, 32 LSTR, 32, 56 LSTR, 51, 11 LSTR, 33, 1

12、2 LSTR, 20, 33 LSTR, 21, 34 LSTR, 34, 8 LSTR, 31, 51 LSTR, 57, 9/PNUM,KP,0 !不顯示關(guān)鍵點(diǎn)編號/PNUM,LINE,1 !顯示線編號/PNUM,AREA,0 !不顯示面編號/PNUM,VOLU,0 !不顯示體編號/PNUM,NODE,0 !不顯示節(jié)點(diǎn)編號!* /PNUM,ELEM,0 !不顯示單元編號/REPLOT !重新繪制!* LPLOT !繪制線圖，帶編號，如圖6-8所示。SAVE !保存數(shù)據(jù)庫創(chuàng)建面。本次的計算區(qū)域為橫向130m，豎向為60m，即左右兩側(cè)計算邊界為4倍左右雙線隧道總跨度，下部邊界為2倍隧道總高度。

13、Al, 35, 1, 36, !采用線創(chuàng)建面，依次創(chuàng)建27個面Al, 36, 2, 34, Al, 34, 9, 37, Al, 37, 3, 35, !左隧道內(nèi)核心土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, !右隧道內(nèi)核心土4個面繪制完成Al, 43, 13, 42, 10,Al, 42, 12, 45, 7, Al, 45,

14、 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, !左側(cè)邊境1個面繪制完成Al,

15、49, 25, 26, 62, Al, 62, 27, 28, 48, !右側(cè)邊境2個面繪制完成/PNUM,KP,0 /PNUM,LINE,0/PNUM,AREA,1APLOT !繪制帶編號的面, 生成的面如圖6-9所示。SAVE !保存數(shù)據(jù),3 創(chuàng)建網(wǎng)格模型設(shè)置線被劃分的單元數(shù)。LESIZE,_Y1, , ,8, , , , ,1lesize,1,8,1 !設(shè)置單元大小, 即L1線劃分成8個單元TYPE, 1 !設(shè)置將要創(chuàng)建單元的類型MAT, 1 !設(shè)置將要創(chuàng)建單元的材料REAL, 1 !設(shè)置將要創(chuàng)建單元的幾何常數(shù)/PNUM,LINE,1 !顯示線編號LPLOT !繪制線圖，帶編號和單元數(shù)。

16、Save !按照“l(fā)esize,1,8,1”格式，將其他61條線也進(jìn)行被劃分單元數(shù)設(shè)置，其線的編號和被劃分的單元數(shù)如表6-3所示，劃分好后如圖6-10所示。表6-1 設(shè)置線單元數(shù)NL(NDIV)NL(NDIV)NL(NDIV)NL(NDIV)NL(NDIV)358398544554514368408584574324348388594604224378418614564524187843447418428108424534214981184544642043816844450419448128314301174581382642922446814825104854948815862102752

17、810338234將所有面劃分單元，其單元圖如圖6-11所示。AMAP,1,14,19,20, !劃分面1，依次單擊關(guān)鍵點(diǎn)14、19、20進(jìn)行映射劃分AMAP,2,14,20,21,AMAP,3,14,21,31,AMAP,4,14,31,19, !面1到4，為左隧道內(nèi)開挖土體AMAP,9,60,46,45,AMAP,10,60,45,44, AMAP,11,60,44,47, AMAP,12,60,47,46, !面9到12，為右隧道內(nèi)開挖土體AMAP,7,21,34,51,31 AMAP,8,31,51,32,19 !面7和8，為左隧道仰拱下土體AMAP,15,44,52,55,47 AM

18、AP,16,47,55,54,46 !面15和16，為右隧道仰拱下土體TYPE, 1 !轉(zhuǎn)換為加固圈土體單元類型和材料MAT, 2AMAP,5,19,32,33,20 !劃分左隧道加固圈土體AMAP,6,20,33,34,21 AMAP,13,46,54,53,45 !劃分右隧道加固圈土體AMAP,14,45,53,52,44TYPE, 1 !轉(zhuǎn)換為普通土體單元類型和材料MAT, 1AMAP,17,32,56,12,33 !劃分左隧道周邊區(qū)域面AMAP,18,33,12,8,34 AMAP,19,34,8,11,51 AMAP,20,51,11,56,32 AMAP,21,54,57,13,5

19、3 !劃分右隧道周邊區(qū)域面AMAP,22,53,13,56,52 AMAP,23,52,56,10,55 AMAP,24,55,10,57,54 AMAP,25,7,8,6,3 !劃分左隧道邊界面AMAP,26,58,2,9,57 !劃分左隧道邊界面AMAP,27,57,9,4,59 SAVE !保存數(shù)據(jù)，本次共劃分成1002個四邊形單元創(chuàng)建模擬噴射混凝土的梁單元。TYPE, 2 !設(shè)置將要創(chuàng)建單元的類型MAT, 3 !設(shè)置將要創(chuàng)建單元的材料REAL, 1 !設(shè)置將要創(chuàng)建單元的幾何常數(shù)E,2,11 !通過兩個節(jié)點(diǎn)創(chuàng)建梁單元E,11,12E,12,13E,13,14E,14,15E,15,16E

20、,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 !保存數(shù)據(jù)E,211,226E,226,225E,225,224E,224,223E,223,222E,222,221E,221,22

21、0E,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 !保存數(shù)據(jù)創(chuàng)建CRD開挖中模擬臨時支撐的梁單元，如圖6-12所示。TYPE, 2 !

22、設(shè)置將要創(chuàng)建單元的類型MAT, 3 !設(shè)置將要創(chuàng)建單元的材料REAL, 1 !設(shè)置將要創(chuàng)建單元的幾何常數(shù)E,1 ,18 !通過兩個節(jié)點(diǎn)創(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 !以上梁單元為左隧道臨時支

23、撐噴射混凝土部分Save !保存數(shù)據(jù)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,

24、218,211Finish !返回上級目錄SAVE !保存數(shù)據(jù)以上梁單元為右隧道臨時支撐噴射混凝土部分，共128個梁單元。6.1.2 加載與自重應(yīng)力場求解/SOL !進(jìn)入求解器NSEL,S,LOC,X,-65.1,-64.9NSEL,A,LOC,X,64.9,65.1 d,all,ux,0 !在選擇的節(jié)點(diǎn)上施加“Ux”約束Allsel !選擇所有內(nèi)容NSEL,S,LOC,Y,-22.5,-22.1d,all,uy,0 !在選擇的節(jié)點(diǎn)上施加“Uy”約束ACEL, 0, 10, 0, !在Y方向施加重力加速度SAVE !施加了位移約束和自重應(yīng)力場后見圖6-13所示。NROPT, FULL, , !

25、采用全牛頓拉普森法進(jìn)行求解Allsel !選擇所有內(nèi)容ESEL,S,TYPE,2 !選擇2類單元Ekill,all !對選擇的單元給予“死屬性”Allsel !選擇所有內(nèi)容Solve !求解計算Finish !求解結(jié)束返回Main Menu 主菜單SAVE豎向位移（PLNSOL,U,Y,0,1）和應(yīng)力（PLNSOL,S,Y,0,1）如圖6-14和圖6-15所示。6.1.3 先左洞后右洞開挖模擬分析1 左洞左上部分土體開挖施工模擬分析開挖土體。/post1Finish !求解結(jié)束返回Main Menu 主菜單/SOL !進(jìn)入求解器Allsel !選擇所有內(nèi)容ANTYPE,REST,1,13,0

26、!重新啟動求解器ASEL,S, , , 2 !選擇左洞左上部分土體面ESLA,R !選擇左洞左上部分土體單元Ekill,all !對選擇的單元給予“死屬性”Save 進(jìn)行噴錨支護(hù)，即激活噴射混凝土襯砌梁單元和臨時支撐單元。Allsel !選擇所有內(nèi)容ESEL,S,TYPE,2 !選擇2類單元Ekill,all !對選擇的單元給予“死屬性”Allsel !選擇所有內(nèi)容lsel,s,2, !選擇梁單元所在的線lsel,a,34,36,2,NSLL,R,1 !選擇線上的節(jié)點(diǎn)ESLN,R,1 !選擇節(jié)點(diǎn)上的單元ESEL,R,TYPE,2 !選擇2類單元Ealive,all !對選擇的單元給予“生屬性”

27、Save 節(jié)點(diǎn)力計算。FINISH /POST1 !進(jìn)入后處理器Allsel !選擇所有內(nèi)容ESEL,S, , , !人工手動選擇開挖土體周圍的一圈單元圖6-16所示NFORCE,ALL !對選擇的節(jié)點(diǎn)求解節(jié)點(diǎn)力Save FINISH /sol !進(jìn)入求解理器節(jié)點(diǎn)力施加，加上節(jié)點(diǎn)力和位移邊界條件后的模型如圖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.2

28、069E+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.1873

29、E+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,

30、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 !選擇所有內(nèi)容Solve !求解計算Finish !求解結(jié)束返回Main Menu 主菜單SAVE后處理，兩個方向的應(yīng)力和位移，如圖6-18圖6-21所示。/post1PLNSOL,U,X,0,1 ! X方向位移PLNS

31、OL,U,Y,0,1 ! Y方向位移PLNSOL,S,X,0,1 ! X方向應(yīng)力PLNSOL,S,Y,0,1 ! Y方向應(yīng)力Finish !求解結(jié)束返回Main Menu 主菜單SAVE計算下一步開挖計算所施加的節(jié)點(diǎn)力，如圖6-22所示。ESEL,S, , , !人工手動選擇開挖土體周圍的一圈單元NFORCE,ALL !對選擇的節(jié)點(diǎn)求解節(jié)點(diǎn)力FINISH !求解結(jié)束返回Main Menu 主菜單Save 2 左洞左下部分土體開挖施工模擬分析開挖土體。/SOL !進(jìn)入求解器Allsel !選擇所有內(nèi)容ANTYPE,REST,2,13,0 !重新啟動求解器ASEL,S, , , 2,3,1 !選擇

32、左洞左上和左下部分土體面ESLA,R !選擇左洞左上和左下部分土體單元Ekill,all !對選擇的單元給予“死屬性”Save 進(jìn)行噴錨支護(hù)，即激活噴射混凝土襯砌梁單元和臨時支撐單元。Allsel !選擇所有內(nèi)容ESEL,S,TYPE,2 !選擇2類單元Ekill,all !對選擇的單元給予“死屬性”Allsel !選擇所有內(nèi)容lsel,s,2,9,7, !選擇梁單元所在的線lsel,a,36,37,1,NSLL,R,1 !選擇線上的節(jié)點(diǎn)ESLN,R,1 !選擇節(jié)點(diǎn)上的單元ESEL,R,TYPE,2 !選擇2類單元Ealive,all !對選擇的單元給予“生屬性”Save 節(jié)點(diǎn)力施加，加上節(jié)點(diǎn)

33、力和位移邊界條件后的模型如圖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.717

34、9E+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,

35、0.8524E+05f,127,fx,0.1220E+06f,127,fy,0.1426E+06f,128,fx,0.1310E+06f,128,fy,0.2152E+06求解計算。Allsel !選擇所有內(nèi)容Solve !求解計算Finish !求解結(jié)束返回Main Menu 主菜單SAVE后處理，兩個方向的應(yīng)力和位移，如圖6-24圖6-27所示。/post1PLNSOL,U,X,0,1 ! X方向位移PLNSOL,U,Y,0,1 ! Y方向位移PLNSOL,S,X,0,1 ! X方向應(yīng)力PLNSOL,S,Y,0,1 ! Y方向應(yīng)力Finish !求解結(jié)束返回Main Menu 主菜單SAV

36、E計算下一步開挖計算所施加的節(jié)點(diǎn)力，如圖6-28所示。ESEL,S, , , !人工手動選擇開挖土體周圍的一圈單元。NFORCE,ALL !對選擇的節(jié)點(diǎn)求解節(jié)點(diǎn)力FINISH Save3 左洞右上部分土體開挖施工模擬分析開挖土體。/SOL !進(jìn)入求解器Allsel !選擇所有內(nèi)容ANTYPE,REST,3,19,0 !重新啟動求解器ASEL,S, , , 1,3,1 !選擇左洞左上、左下和右上部分土體面ESLA,R !選擇左洞左上、左下和右上部分土體單元Ekill,all !對選擇的單元給予“死屬性”Save 進(jìn)行噴錨支護(hù)，即激活噴射混凝土襯砌梁單元和臨時支撐單元。Allsel !選擇所有內(nèi)容

37、ESEL,S,TYPE,2 !選擇2類單元Ekill,all !對選擇的單元給予“死屬性”Allsel !選擇所有內(nèi)容lsel,s,1,2,1, !選擇梁單元所在的線lsel,a,35,37,2,lsel,a,9,NSLL,R,1 !選擇線上的節(jié)點(diǎn)ESLN,R,1 !選擇節(jié)點(diǎn)上的單元ESEL,R,TYPE,2 !選擇2類單元Ealive,all !對選擇的單元給予“生屬性”Save 節(jié)點(diǎn)力施加，加上節(jié)點(diǎn)力和位移邊界條件后的模型如圖6-29所示。f,1,fx,8888 !施加X方向集中力f,1,fy,0.1627E+05 !施加Y方向集中力f,2,fx,-0.3130E+05f,2,fy,0.4

38、226E+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+05 f,10,fx,-0.1884E+05f,10,fy,-0.2575E+06f,11,fx,-5824f,11,fy,

39、-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 !選擇所有內(nèi)容Solve !求解計算Finish !求解結(jié)束返回Main Menu 主菜單SAVE后處理，兩個方向的應(yīng)力和位移，如圖6-30圖6-33所示。/post1PLNSOL,U,X,0,1 ! X方向位移PLNSOL,U,Y,0,1 ! Y方向位移PLNSOL,S,X,0,1 ! X方向應(yīng)力PLNSOL,S,Y,0,1 ! Y方向應(yīng)力Finish !求解結(jié)束返回Main Menu 主菜單SAVE（6）計算下一步開挖計算所施加的節(jié)