隧道及地下工程ANSYS分析 第6章 高速公路分離式偏壓隧道

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í)體單元類型KEYOPT,1,1,0KEYOPT,1,2,0KEYOPT,1,3

2、,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,DENS,1,2000 !輸入密度（圍巖）TB,DP,1,

3、!采用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, MPDATA,EX,3,27.5e9 !輸入彈性模量（噴射

4、混凝土）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)建計(jì)算區(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,K,14,-15,1.169, !創(chuàng)建隧道輪廓線關(guān)鍵點(diǎn)

5、，上半圓圓心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,7LDELE,13,16,3NUMMRG,ALL, ,

6、 , ,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,19 !刪除部分線LDELE,5,7,1LDELE,10

7、,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,12,2ldele,13,15,1LOVLAP,6,11 !作疊

8、加布爾運(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, ,0 FLST,2,5,4,ORDE,5 !做線加布爾運(yùn)算FI

9、TEM,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)備k,57,30.044,1.169,0,k,58,30.044

10、,-22.3,0,k,59,30.044,37.1,0,k,60,15,1.169,0,LSTR, 3, 7 !畫計(jì)算區(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, 14 LSTR, 14, 19 LSTR, 14, 20 LST

11、R, 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, 12 LSTR, 20, 33 LSTR, 21, 34 LSTR

12、, 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)建面。本次的計(jì)算區(qū)域?yàn)闄M向130m，豎向?yàn)?0m，即左右兩側(cè)計(jì)算邊界為4倍左右雙線隧道總跨度，下部邊界為2倍隧道總高度。Al, 35, 1, 36, !采用線創(chuàng)建面，依次創(chuàng)建27個面A

13、l, 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, 44, 11, Al, 44, 15, 43, 16,

14、 !右隧道加固圈土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, 49, 25, 26, 62, Al, 62, 27, 28,

15、48, !右側(cè)邊境2個面繪制完成/PNUM,KP,0 /PNUM,LINE,0/PNUM,AREA,1APLOT !繪制帶編號的面, 生成的面如圖6-9所示。SAVE !保存數(shù)據(jù),1創(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ù)。Save !按照“l(fā)esize,1,8,1”格式，將其他61條線也

16、進(jìn)行被劃分單元數(shù)設(shè)置，其線的編號和被劃分的單元數(shù)如表6-3所示，劃分好后如圖6-10所示。表6-1 設(shè)置線單元數(shù) NL(NDIVNL(NDIVNL(NDIVNL(NDIVNL(NDIV358398544554514368408584574324348388594604224378418614564524187843447418428108424534214981184544642043816844450419448128314301174581382642922446814825104854948815862102752810338234將所有面劃分單元，其單元圖如圖6-11所示。AMAP,1

17、,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 AMAP,16,47,55,54,46 !面15和16，為右隧道仰拱下土體TY

18、PE, 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,53 !劃分右隧道周邊區(qū)域面AMAP,22,53,13,56,52 AMAP

19、,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,16,17E,17,10E,10,63E,63,64E,64,65E,6

20、5,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,220E,220,219E,219,278E,278,277E,277,276

21、E,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 !設(shè)置將要創(chuàng)建單元的類型MAT, 3 !設(shè)置將要創(chuàng)建單元的材料REAL, 1

22、 !設(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 !以上梁單元為左隧道臨時支撐噴射混凝土部分Save !保存數(shù)據(jù)E,210,227E,227,228E

23、,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 !保存數(shù)據(jù)以上梁單元為

24、右隧道臨時支撐噴射混凝土部分，共128個梁單元。6.1.1 加載與自重應(yīng)力場求解/SOL !進(jìn)入求解器NSEL,A,LOC,X,64.9,65.1 d,all,ux,0 !在選擇的節(jié)點(diǎn)上施加“Ux”約束Allsel !選擇所有內(nèi)容d,all,uy,0 !在選擇的節(jié)點(diǎn)上施加“Uy”約束ACEL, 0, 10, 0, !在Y方向施加重力加速度SAVE !施加了位移約束和自重應(yīng)力場后見圖6-13所示。NROPT, FULL, , !采用全牛頓拉普森法進(jìn)行求解Allsel !選擇所有內(nèi)容ESEL,S,TYPE,2 !選擇2類單元Ekill,all !對選擇的單元給予“死屬性”Allsel !選擇所有內(nèi)

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

26、錨支護(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 !對選擇的單元給予“生屬性”Save 節(jié)點(diǎn)力計(jì)算。FINISH /POST1 !進(jìn)入后處理器Allsel !選擇所有內(nèi)容ESEL,S, , , !人工手動選擇開挖土體周圍的一圈單元圖6-16所示NFO

27、RCE,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.2069E+05f, 19 ,fx, -0.1041E+06f, 19 ,fy, -0.2041E+05f, 20 ,fx, -0.

28、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.2

29、051E+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.1

30、940E+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求解計(jì)算。Allsel !選擇所有內(nèi)容Solve !求解計(jì)算Finish !求解結(jié)束返回Main Menu 主菜單SAVE后

31、處理，兩個方向的應(yīng)力和位移，如圖6-18圖6-21所示。/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計(jì)算下一步開挖計(jì)算所施加的節(jié)點(diǎn)力，如圖6-22所示。ESEL,S, , , !人工手動選擇開挖土體周圍的一圈單元NFORCE,ALL !對選擇的節(jié)點(diǎn)求解節(jié)點(diǎn)力FINISH !求解結(jié)束返回Main Menu 主菜單Save 2左洞左下部分土體開挖施工模擬分析開挖土體。/SOL !進(jìn)入求解器All

32、sel !選擇所有內(nèi)容ANTYPE,REST,2,13,0 !重新啟動求解器ASEL,S, , , 2,3,1 !選擇左洞左上和左下部分土體面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)上的單元ES

33、EL,R,TYPE,2 !選擇2類單元Ealive,all !對選擇的單元給予“生屬性”Save 節(jié)點(diǎn)力施加，加上節(jié)點(diǎn)力和位移邊界條件后的模型如圖6-23所示。f, 1 ,fx, -0.7585E+05 !施加X方向集中力f, 1 ,fy, 0.5203E+05 !施加Y方向集中力f, 62 ,fx, 0.7840E+05f, 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,

34、 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, 12

35、3 ,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求解計(jì)算。Allsel !選擇所有內(nèi)容Solve !求解計(jì)算Finish !求解結(jié)束返回Main Menu 主菜單SAVE后處理，兩個方向

36、的應(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 主菜單SAVE計(jì)算下一步開挖計(jì)算所施加的節(jié)點(diǎn)力，如圖6-28所示。ESEL,S, , , !人工手動選擇開挖土體周圍的一圈單元。NFORCE,ALL !對選擇的節(jié)點(diǎn)求解節(jié)點(diǎn)力FINISH Save1左洞右上部分土體開挖施工模擬分析開挖土體。/SOL !進(jìn)入求解器Allsel !選擇所有內(nèi)容ANTYPE,REST,3,19

37、,0 !重新啟動求解器ASEL,S, , , 1,3,1 !選擇左洞左上、左下和右上部分土體面ESLA,R !選擇左洞左上、左下和右上部分土體單元Ekill,all !對選擇的單元給予“死屬性”Save 進(jìn)行噴錨支護(hù)，即激活噴射混凝土襯砌梁單元和臨時支撐單元。Allsel !選擇所有內(nèi)容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

38、!選擇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.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,

39、 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, -0.1401E+06f, 12 ,fx, -0.1110E+05f, 12 ,fy, -0.1420E+06f, 13 ,fx, -0.2394E+05f, 13 ,fy, -0

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

41、方向位移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）計(jì)算下一步開挖計(jì)算所施加的節(jié)點(diǎn)力，如圖6-34所示。ESEL,S, , , !人工手動選擇開挖土體周圍的一圈單元。NFORCE,ALL !對選擇的節(jié)點(diǎn)求解節(jié)點(diǎn)力FINISH Save2左洞右下部分土體開挖施工模擬分析（1）開挖土體。/SOL !進(jìn)入求解器Allsel !選擇所有內(nèi)容ANTYPE,REST,4,13,0 !重新啟動求解器ASEL,S, , , 1,4,1 !選擇左洞全部核心土

42、體面ESLA,R !選擇左洞全部核心土體單元Ekill,all !對選擇的單元給予“死屬性”Save 提示：在開挖左隧道右下部分土體時，同樣要進(jìn)行第一部分（左上）、第二部分（左下）和第三部分（右上）上所有單元的殺死、梁單元的激活操作，這是重復(fù)的，但重新啟動后必須進(jìn)行。（2）進(jìn)行噴錨支護(hù)，即激活噴射混凝土襯砌梁單元和臨時支撐單元。Allsel !選擇所有內(nèi)容ESEL,S,TYPE,2 !選擇2類單元Ekill,all !對選擇的單元給予“死屬性”Allsel !選擇所有內(nèi)容lsel,s,1,2,1, !選擇梁單元所在的線lsel,a,3,9,6,NSLL,R,1 !選擇線上的節(jié)點(diǎn)ESLN,R,1

43、 !選擇節(jié)點(diǎn)上的單元ESEL,R,TYPE,2 !選擇2類單元Ealive,all !對選擇的單元給予“生屬性”Save （3）節(jié)點(diǎn)力施加，加上節(jié)點(diǎn)力和位移邊界條件后的模型如圖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 ,

44、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）求解計(jì)算。Allsel !選擇所有內(nèi)容Solve !求解計(jì)算Finish !求解結(jié)束返回Main Menu 主菜單SAVE（5）后處理，兩個方向的應(yīng)力和位移，

45、如圖6-36圖6-42所示。內(nèi)力和變形如圖6-43圖6-47所示。/post1PLDISP,1 !繪制變形圖PLNSOL,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)力PLNSOL,S,1,0,1 !第一主應(yīng)力PLNSOL,S,3,0,1 !第三主應(yīng)力p PLNSOL,EPPL,EQV,0,1 ! 塑性應(yīng)變p ETABLE, ,SMISC,1 !span創(chuàng)建梁單元內(nèi)力表p ETABLE, ,SMISC, 7 p ETABLE, ,SMISC, 2 ETABLE, ,SMISC,8

46、ETABLE, ,SMISC,6ETABLE, ,SMISC,12 ESEL,S,LIVE !選擇“屬性為活”的單元ESEL,R,TYPE,2 !選擇單元類型為2的所有單元/REPLOT !從新繪制圖形PLLS,SMIS1,SMIS7,0.2,0 !畫軸力圖，顯示比例為PLLS,SMIS2,SMIS8,0.2,0 !畫剪力圖，顯示比例為PLLS,SMIS6,SMIS12,-0.2,0 !畫彎矩圖，顯示比例為SET,LIST !設(shè)置計(jì)算步驟LCDEF,1,1,13,LCDEF,5,5,13,LCASE,5,LCOPER,SUB,1, , , !將第四步的結(jié)果減去第一部（自重應(yīng)力場）的結(jié)果PLDI

47、SP,1 ! 畫變形噴射混凝土的圖Finish !求解結(jié)束返回Main Menu 主菜單SAVE（6）計(jì)算下一步開挖計(jì)算所施加的節(jié)點(diǎn)力，如圖6-48所示。ESEL,S, , , !人工手動選擇開挖土體周圍的一圈單元。NFORCE,ALL !對選擇的節(jié)點(diǎn)求解節(jié)點(diǎn)力FINISH Save2右洞左上部分土體開挖施工模擬分析（1）開挖土體。/post1Finish !求解結(jié)束返回Main Menu 主菜單/SOL !進(jìn)入求解器Allsel !選擇所有內(nèi)容ANTYPE,REST,5,13,0 !重新啟動求解器ASEL,S, , , 1,4,1 !選擇左洞全部土體面ASEL,a, , , 10 !選擇右洞

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

49、單元Ealive,all !對選擇的單元給予“生屬性”Save（3）節(jié)點(diǎn)力施加，加上節(jié)點(diǎn)力和位移邊界條件后的模型如圖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+0

50、5f, 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, -3784項(xiàng)f, 233 ,fx, -0.9096E+05 f, 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 ,f

51、y, -0.1087E+06f, 274 ,fx, 0.1408E+06f, 274 ,fy, -0.1748E+06工程名稱：萬科天譽(yù)花園二期f, 275 ,fx, 10f, 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,

52、 280 ,fy, 0.3819E+05f, 281 ,fx, 2381f, 281 ,fy, 0.3851E+05f, 282 ,fx,f, 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）求解計(jì)算。Allsel !選擇所有內(nèi)容Solve !求解計(jì)算Finish !求解結(jié)束返回Main Menu 主菜單SAVE（5）后處理，兩個方向的應(yīng)力和位移，如圖圖所示。頁/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）計(jì)算下一步開挖計(jì)算所施加的節(jié)點(diǎn)力，如圖6-54所示。ESEL,S, , , !人工手動選擇開挖土體周圍的一圈單元。萬科天譽(yù)花園二期FINISH !求解結(jié)束返回Main Menu 主菜單Save 3右洞左下部分土體開挖施工模擬分析（1）開挖土體。/SOL !進(jìn)入求解器Allsel !選擇所有內(nèi)容ANTYPE,REST,6,14,0