隧道及地下工程ANSYS實(shí)例分析雙線鐵路隧道襯砌結(jié)構(gòu)設(shè)計(jì)力學(xué)分析

1、CM CO 寸 9N- 009if)51二次襯砌結(jié)構(gòu)力學(xué)分析/TITLE , Mechanical analysis on railway tunnel 2nd lining ! 確定分析標(biāo)題/NOPR!菜單過濾設(shè)置/PMETH , OFF, 0KEYW, PR_SET, 1KEYW , PR_STRUC , 1!保留結(jié)構(gòu)分析部分菜單/COM ,/COM , Preferences for GUI filtering have been set to display:/COM , Structural/PREP7!進(jìn)入前處理器ET, 1, BEAM3!設(shè)置梁?jiǎn)呜ｎ愋虴T, 2, COMBIN1

2、4!設(shè)置彈簧單兀類型R, 1, 0.4, 0.0053333, 0.4,!設(shè)置梁?jiǎn)卧獛缀纬?shù)R, 2, 400e6,!設(shè)置彈簧單元幾何常數(shù)MPTEMP,!設(shè)置材料模型MPTEMP ,1,0MPDATA , EX, 1, , 30.0e9!輸入彈性模量MPDATA , PRXY , 1 , , 0.2!輸入泊松比MPTEMP,!設(shè)置材料模型MPTEMP ,1,0MPDATA , DENS, 1 , , 2500!輸入密度SAVE!保存數(shù)據(jù)庫(kù)1.建立幾何模型K, 1,!創(chuàng)建關(guān)鍵點(diǎn)（隧道二次襯砌）K, 2, 4.71,-1.82,K, 3,-4.71,-1.82,K, 4, 0, -3.75,K,

3、5, 0, 5.05,K, 20, 5.6, -2.162,!創(chuàng)建關(guān)鍵點(diǎn)（地層彈簧）K, 30, -5.6, -2.16,K, 40, 0, -4.75,K, 50, 0, 6,!創(chuàng)建隧道襯砌線LARC, 2, 5, 1, 5.05,!創(chuàng)建圓弧線（拱頂部）LARC, 5, 3, 1, 5.05,LARC,2,3,4!創(chuàng)建圓弧線（仰拱部）!創(chuàng)建地層彈簧線LARC, 20, 50, 1,6,!創(chuàng)建圓弧線（拱頂部）LESIZE,ALL, , ,10, ,1, , ,1,! 設(shè)置單元大小, 每條弧線劃分成TYPE,1! 設(shè)置將要?jiǎng)?chuàng)建單元的類型MAT,1! 設(shè)置將要?jiǎng)?chuàng)建單元的材料REAL,1! 設(shè)置將要

4、創(chuàng)建單元的幾何常數(shù)LMESH,1,6,1! 將所有直線劃分單元/PNUM, KP, 0/PNUM, ELEM, 1! 以下為顯示單元編號(hào)和顏色/REPLOT SAVE! 重新顯示TYPE,2! 設(shè)置將要?jiǎng)?chuàng)建單元的類型MAT,1! 設(shè)置將要?jiǎng)?chuàng)建單元的材料REAL,2! 設(shè)置將要?jiǎng)?chuàng)建單元的幾何常數(shù)E, 1,31! 通過兩個(gè)節(jié)點(diǎn)創(chuàng)建彈簧單元2 單元網(wǎng)格劃分設(shè)置單元大小并將所有直線劃分單元，其單元圖如圖5-8 所示。10個(gè)單元LARC,50, 30, 1, 6,LARC,20,SAVE30,40 !創(chuàng)建圓弧線（仰拱部）! 保存數(shù)據(jù)E, 2,32 E,3,33 E,4,34 E,5,35 E,6,36

5、E,7,37 E,8,38 E,9,39 E,10,40 E,11,41 E,12,42 E,13,43 E,14,44 E,15,45 E,16,46 E,17,47 E,18,48 E,19,49 E,20,50 E,21,51 E,22,52 E,23,53 E,24,54 E,25,55 E,26,56 E,27,57E,28,58E,29,59E,30,60Lclear,4,6,1LDELE,4, 6,1 ,1 Finish! 清除用于創(chuàng)建地層彈簧單元的外層梁?jiǎn)卧? 返回 Main Menu 主菜單5.18.2加載 與 求 解/SOLd,31,ux,0,0,60,1,uy,ACEL,

6、 0, 10, 0,F,2,Fx,0F,13,Fx,43554.43356F,14,Fx,86574.64508F,15,Fx,128464.6048F,16,Fx,168285.8813F,17,Fx,204799.5937F,18,Fx,236448.4444F,19,Fx,261396.9365F,20,Fx,277795.0781F,21,Fx,283894.0611F,12,Fx,247397.3915F,30,Fx,190308.0663F,29,Fx,149683.6899F,28,Fx,103226.8339F,27,Fx,52665.0399F,26,Fx,0F,25,Fx,

7、-52665.0399F,24,Fx,-103226.8339F,23,Fx,-149683.6899F,22,Fx,-190308.0663F,1,Fx,-247397.3915F,3,Fx,-283894.0611F,4,Fx,-277795.0781F,5,Fx,-261396.9365F,6,Fx,-236448.4444F,7,Fx,-204799.5937F,8,Fx,-168285.8813F,9,Fx,-128464.6048F,10,Fx,-86574.64508F,11,Fx,-43554.43356F,2,Fy,-449636.88F,13,Fy,-441210F,14,

8、Fy,-416229.66!在“Ux”和“ Uy”兩個(gè)方向的施加約束!在 Y 方向施加重力加速度! 在節(jié)點(diǎn)上施加X 方向集中力Y 方向集中力F,15,Fy,-375629.1F,16,Fy,-320951.4F,17,Fy,-254238.6F,18,Fy,-178008.6F,19,Fy,-95102.7F,20,Fy,0F,21,Fy,82772.64F,12,Fy,259324.92F,30,Fy,413398.46F,29,Fy,454320.04F,28,Fy,484308F,27,Fy,502604.08F,26,Fy,508719.08F,25,Fy,502604.08F,24,

9、Fy,484308F,23,Fy,454320.04F,22,Fy,413398.46F,1,Fy,259324.92F,3,Fy,82772.64F,4,Fy,0F,5,Fy,-95102.7F,6,Fy,-178008.6F,7,Fy,-254238.6F,8,Fy,-320951.4F,9,Fy,-375629.1F,10,Fy,-416229.66F,11,Fy,-441210SaveNROPT, FULL,AllselOutres, all, allSolveFinishSAVE! 采用全牛頓拉普森法進(jìn)行求解! 選擇所有內(nèi)容! 輸出所有內(nèi)容! 求解計(jì)算! 求解結(jié)束返回Main Men

10、u 主菜單5.18.3后 處 理1 初次查看內(nèi)力和變形結(jié)果/POST1! 進(jìn)入后處理器PLDISP, 1! 繪制變形和未變形圖ETABLE, SMISC,6! 6 、 12 表示彎矩ETABLE, SMISC,12ETABLE, SMISC,1! 1 、 7 表示軸力ETABLE, SMISC,7ETABLE, SMISC,2! 2、 8 表示剪力ETABLE, , SMISC, 8查看內(nèi)力, 包括彎矩、軸力和剪力，如圖 5-12、 圖 5-13 和圖 5-14 所示。 路徑： GeneralPostproc Plot Results Contour Plot Line Elem Res 。!

11、 繪制彎矩圖! 僅顯示單元類型1! 繪制軸力圖! 繪制剪力圖PLLS, SMIS6, SMIS12, -1, 0ESEL,R,TYPE,1PLLS, SMIS1, SMIS7, 1, 0PLLS, SMIS2, SMIS8, 1, 02 去除受拉彈簧再計(jì)算! 結(jié)束后處理器操作! 進(jìn)入求解器! 殺死 62 號(hào)地層彈簧單元Finish /solEkill,62 Ekill,71 Ekill,70 Ekill,73 Ekill,74 Ekill,88 Ekill,86Ekill,87 Ekill,85 Ekill,86 SAVEAllsel!選擇所有內(nèi)容Solve!求解計(jì)算Finish!求解結(jié)束返回

12、Main Menu主菜單SAVE3 查看最后計(jì)算結(jié)果/POST1ETABLE, REFL!更新單元表數(shù)據(jù)PRETAB, SMIS6, SMIS12, SMIS1, SMIS7, SMIS2,SMIS8!打印單元表數(shù)據(jù)PRNSOL, DOF,!打印節(jié)點(diǎn)位移5.19臺(tái)階法施工過程仿真分析1. /TITLE , /NOPR /PMETH KEYW , KEYW , /COM ,材料、實(shí)常數(shù)和單元類型定義Mechanical analysis on railway tunnel 1nd liningOFF 0PR_SET, 1PR_STRUC, 1Preferences for GUI filteri

13、ng have been set to display:!確定分析標(biāo)題!菜單過濾設(shè)置!保留結(jié)構(gòu)分析部分菜單/COM , Structural /PREP7 ET,1,PLANE42 KEYOPT,1,1,0 KEYOPT,1,2,0 KEYOPT,1,3,2 KEYOPT,1,5,0 KEYOPT,1,6,0 ET, 2, BEAM3 ET,3,LINK1 R,1,0.2,0.0008,0.2, R20.0005, MPTEMP, MPTEMP,1,0 MPDATA,EX,1,3.2e9 MPDATA,PRXY,1,0.32 MPDATA,DENS,1,2200 MPDATA,EX,2,27.

14、5e9 MPDATA,PRXY,2,0.2 MPDATA,DENS,2,2500 MPDATA,EX,3,200e9 MPDATA,PRXY,3,0.3 MPDATA,DENS,3,7800!進(jìn)入前處理器!設(shè)置實(shí)體單元類型!設(shè)置為平面應(yīng)變模式!設(shè)置梁?jiǎn)卧愋?設(shè)置桿單元類型!設(shè)置梁?jiǎn)卧獛缀纬?shù)!設(shè)置桿單元幾何常數(shù)SAVE2.建立幾何模型!輸入彈性模量（圍巖）!輸入泊松比（圍巖）!輸入密度（圍巖）!輸入彈性模量（噴射混凝土） !輸入泊松比（噴射混凝土） !輸入密度（噴射混凝土）!輸入彈性模量（錨桿）!輸入泊松比（錨桿）!輸入密度（錨桿）!保存數(shù)據(jù)庫(kù)K,1,0,0, K,2,5.45,0, K,3

15、,0,5.45, K,4,-5.45,0, K,5,-5.0838,-1.964, K,6,0,-4.15,!創(chuàng)建隧道輪廓線關(guān)鍵點(diǎn)K,7,5.0838,-1.964,K,8,8.45,0,!創(chuàng)建錨桿加固區(qū)域關(guān)鍵點(diǎn)K,9,0,8.45,K,10,-8.45,0,K,11,13,0,K,12,0,13, K,13,-13,0, K,14,-13,13, K,15,13,13, K,16,8.45,-8.45, K,17,-8.45,-8.45, K,18,-13,-13, K,19,13,-13,K,20,0,26,!創(chuàng)建計(jì)算區(qū)域關(guān)鍵點(diǎn)K,21,13,26, K,22,-13,26, K,23,-

16、13,-34, K,24,13,-34, K,25,0,-34, K,26,50,-34, K,27,-50,-34, K,28,-50,26, K,29,50,26, K,30,50,13, K,31,50,0, K,32,50,-13, K,33,-50,-13, K,34,-50,0, K,35,-50,13, K,36,0,-8.45K,37,0,-13K,38,2.8653,-3.5371 K,39,-2.8653,-3.5371 SAVELARC,2,3,1,5.45,!* 畫隧道輪廓線LARC,3,4,1,5.45,!*LARC,4,5,1,5.45,LARC,5,6,39!*L

17、ARC,7,2,1,5.45,LARC,8,9,1,8.45,!* 畫錨桿加固區(qū)域線LARC,9,10,1,8.45,LSTR,1,2!* 畫其它分割線LSTR,1,3LSTR,1,4LSTR,1,6LSTR,2,8LSTR,4,10LSTR,8,11LSTR,11,15LSTR,15,12LSTR,12,9LSTR,12,14LSTR,14,13LSTR,13,10LSTR,10,17LSTR,13,18LSTR,8,16LSTR,11,19LSTR,17,36LSTR,36,16LSTR,6,36LSTR,36,37LSTR,18,37LSTR,37,19LSTR,37,25LSTR,19

18、,24LSTR,18,23LSTR,14,22LSTR,12,20LSTR,15,21LSTR,22,20LSTR,20,21LSTR,23,25LSTR,25,24LSTR,24,26LSTR,19,32LSTR,11,31LSTR,15,30LSTR,21,29LSTR,22,28LSTR,14,35LSTR,13,34LSTR,18,33LSTR,23,27LSTR,26,32LSTR,32,31LSTR,31,30LSTR,30,29LSTR,27,33LSTR,33,34LSTR,34,35LSTR,35,28LSTR,3,9LARC,6,7,38!*SAVE!保存數(shù)據(jù)，生成的直線如

19、圖5-20所示。Al,8,1,9,!采用線創(chuàng)建面，依次創(chuàng)建 24個(gè)面Al,9,2,10,Al,10,3,4,11,Al,11,60,5,8,Al,12,6,59,1,Al,59,7,13,2,Al,27,26,23,12,5,60,Al,27,25,21,13,3,4,Al,14,15,16,17,6,Al,17,7,20,19,18,Al,20,21,25,28,29,22,Al,30,24,14,23,26,28,Al,39,31,29,33,Al,40,32,30,31,Al,41,51,42,32,Al,42,52,43,24,Al,43,53,44,15,Al,44,54,45,36

20、,Al,16,35,38,36,Al,18,34,37,35,Al,47,58,46,34,Al,48,57,47,19,Al,49,56,48,22,Al,50,55,49,33,SAVE!保存數(shù)據(jù)3.創(chuàng)建網(wǎng)格模型lesize,1,8,1!設(shè)置單兀大小，即L1線劃分成8個(gè)單兀TYPE,1!設(shè)置將要?jiǎng)?chuàng)建單兀的類型MAT,1!設(shè)置將要?jiǎng)?chuàng)建單元的材料REAL,1!設(shè)置將要?jiǎng)?chuàng)建單元的幾何常數(shù)NL(NDIV)N L(NDIV)N L(NDIV)N L(NDIV)N L(NDIV)N L(NDIV)18118214315418515281232243254285243213323433543853446

21、1432443444485445215425435445855568164264364468564781732733744785748818428338448858498194294394P 498593108203304404508606!按照“ lesize,1,8,1”格式，將其他59條線也進(jìn)行被劃分單元數(shù)設(shè)置AMAP,1,1,2,3, AMAP,2,1,3,4, AMAP,3,1,4,6, AMAP,4,6,2,1, AMAP,5,2,8,9,3 AMAP,6,3,9,10,4 AMAP,7,6,36,8,2AMAP,8,4,10,36,6AMAP,9,8,11,12,9 AMAP,1

22、0,9,12,13,10 AMAP,11,10,13,37,36 AMAP,12,36,37,11,8劃分面 1 ，依次單擊關(guān)鍵點(diǎn)1 ， 2， 3 進(jìn)行映射面 1 到4，為隧道內(nèi)開挖土體劃分錨固加固去圍巖單元，包括面5 和 6采用映射劃分中心區(qū)域其它面amesh,13,24,1! 劃分周邊區(qū)域面SAVE! 存儲(chǔ)數(shù)據(jù)TYPE,2! 設(shè)置將要?jiǎng)?chuàng)建單元的類型MAT,2! 設(shè)置將要?jiǎng)?chuàng)建單元的材料REAL,1! 設(shè)置將要?jiǎng)?chuàng)建單元的幾何常數(shù)! 通過兩個(gè)節(jié)點(diǎn)創(chuàng)建梁?jiǎn)卧狤, 2,11E,10,63E,62,115E,116,167E,11,12E,63,64E,115,114E,167,168E,12,13E

23、,64,65E,114,117E,168,169E,13,14E,65,66E,117,118E,169,170E,14,15E,66,67E,118,119E,170,171E,15,16E,67,68E,119,120E,171,166E,16,17E,68,69E,120,121E,166,172E,17,10E,69,62E,121,116E,172,2TYPE,3! 設(shè)置將要?jiǎng)?chuàng)建單元的類型MAT,3! 設(shè)置將要?jiǎng)?chuàng)建單元的材料REAL,2! 設(shè)置將要?jiǎng)?chuàng)建單元的幾何常數(shù)! 通過兩個(gè)節(jié)點(diǎn)創(chuàng)建梁?jiǎn)卧狤,2,211E,233,217E,247,254E,68,252E,211,212E,15,

24、227E,254,238E,252,259E,212,210E,227,234E,64,248E,259,243E,11,223E,234,218E,248,255E,69,253E,223,230E,16,228E,255,239E,253,260E,230,214E,228,235E,65,249E,260,244E,12,224E,235,219E,249,256E,62,245E,224,231E,17,229E,256,240E,245,246E,231,215E,229,236E,66,250E,246,237E,13,225E,236,220E,250,257FinishE,22

25、5,232E,10,221E,257,241E,232,216E,221,222E,67,251E,14,226E,222,213E,251,258E,226,233E,63,247E,258,2425.19.2加 載 與 自 重 應(yīng) 力 場(chǎng) 求 解/SOLNSEL,S,LOC,X,-50.1,-49.9NSEL,A,LOC,X,49.9,50.1 d,all,ux,0AllselNSEL,S,LOC,Y,-34.1,-33.9 d,all,uy,0 （ 2）施加重力加速度。ACEL, 0, 10, 0, SAVENROPT, FULL, ,AllselESEL,S,TYPE,2 ESEL,A

26、,TYPE,3 Ekill,all Allsel Solve Finish SAVE! 進(jìn)入求解器! 在選擇的節(jié)點(diǎn)上施加“Ux ”約束! 選擇所有內(nèi)容! 在選擇的節(jié)點(diǎn)上施加“Uy ”約束!在 Y 方向施加重力加速度! 施加了位移約束和自重應(yīng)力場(chǎng)后見圖5-25 所示。! 采用全牛頓拉普森法進(jìn)行求解! 選擇所有內(nèi)容! 選擇2 類單元! 選擇3 類單元! 對(duì)選擇的單元給予“死屬性”! 選擇所有內(nèi)容! 求解計(jì)算! 求解結(jié)束返回Main Menu 主菜單1 自重應(yīng)力場(chǎng)后處理豎向位移( PLNSOL,U,Y,0,1 ) 和應(yīng)力 ( PLNSOL,S,Y,0,1 ) 如圖 5-26 和圖 5-27 所示。5

27、.19.3 上 臺(tái) 階 開 挖 模 擬 分 析1 上臺(tái)階開挖施工/post1Finish! 求解結(jié)束返回Main Menu 主菜單/SOL! 進(jìn)入求解器Allsel! 選擇所有內(nèi)容ANTYPE,REST,1,1,0! 重新啟動(dòng)求解器ASEL,S, , ,1,2,1! 選擇上臺(tái)階面ESLA,R! 選擇上臺(tái)階土體單元Ekill,all! 對(duì)選擇的單元給予“死屬性”Allsel! 選擇所有內(nèi)容lsel,s,1,2,1,NSLL,R,1ESLN,RESEL,R,TYPE,2! 選擇2 類單元ESEL,A,TYPE,3! 選擇3 類單元Ealive,all! 對(duì)選擇的單元給予“生屬性”口下點(diǎn)福裁2 施加

28、節(jié)點(diǎn)荷載FINISH/POST1! 進(jìn)入后處理器Allsel! 選擇所有內(nèi)容ESEL,S, , ,! 人工手動(dòng)選擇開挖土體周圍的一圈單元，如圖NFORCE,ALL! 對(duì)選擇的節(jié)點(diǎn)求解節(jié)點(diǎn)力FINISH /sol! 進(jìn)入求解理器節(jié)點(diǎn)力施加，加上節(jié)點(diǎn)力和位移邊界條件后的模型如圖5-29 所示。f, 1 ,fx,4.37E-09! 施加X 方向集中力f, 2 ,fx,-1.41E+05f, 3 ,fx,-67.92f, 4 ,fx,-201.3f, 5 ,fx,-357.9f, 6 ,fx,-499.4f, 7 ,fx,-616.4f, 8 ,fx,-664.6f, 9 ,fx,-617.6f, 1

29、0 ,fx,6.73E-05f, 11 ,fx,-2.70E+05f, 12 ,fx,-2.44E+05f, 13 ,fx,-2.11E+05f, 14 ,fx,-1.73E+05f, 15 ,fx,-1.32E+05f, 16 ,fx,-8.90E+04f, 17 ,fx,-4.47E+04f, 62 ,fx,1.41E+05f, 63 ,fx,4.47E+04f, 64 ,fx,8.90E+04f, 65 ,fx,1.32E+05f, 66 ,fx,1.73E+05f, 67 ,fx,2.11E+05f, 68 ,fx,2.44E+05f, 69 ,fx,2.70E+05f, 70 ,fx

30、,67.92f, 71 ,fx,201.3f, 72 ,fx,357.9f, 73 ,fx,499.4f, 74 ,fx,616.4f, 75 ,fx,664.6f, 76 ,fx,617.6f, 1 ,fy,3.90E+05! 施加Y 方向集中力f, 2 ,fy,1.66E+05f, 3 ,fy,3.90E+055-28 所示。f,4,fy,3.90E+05f,5,fy,3.90E+05f,6,fy,3.90E+05f,7,fy,3.90E+05f,8,fy,3.90E+05f,9,fy,3.90E+05f,10,fy,-4.82E+05f,11,fy,-1.12E+05f,12,fy,-2

31、.13E+05f,13,fy,-2.98E+05f,14,fy,-3.67E+05f,15,fy,-4.18E+05f,16,fy,-4.54E+05f,17,fy,-4.75E+05f,62,fy,1.66E+05f,63,fy,-4.75E+05f,64,fy,-4.54E+05f,65,fy,-4.18E+05f,66,fy,-3.67E+05f,67,fy,-2.98E+05f,68,fy,-2.13E+05f,69,fy,-1.12E+05f,70,fy,3.90E+05f,71,fy,3.90E+05f,72,fy,3.90E+05f,73,fy,3.90E+05f,74,fy,3

32、.90E+05f,75,fy,3.90E+05f,76,fy,3.90E+053上臺(tái)階開挖求解Allsel! 選擇所有內(nèi)容Solve! 求解計(jì)算Finish! 求解結(jié)束返回Main Menu 主菜單SAVE4上臺(tái)階開挖模擬分析后處理/post1PLNSOL,U,Y,0,1! Y 方向位移PLNSOL,S,Y,0,1! Y 方向應(yīng)力Finish! 求解結(jié)束返回Main Menu 主菜單SAVE5.19.4 下 臺(tái) 階 開 挖 模 擬 分 析/SOL! 進(jìn)入求解器Allsel! 選擇所有內(nèi)容ANTYPE,REST,2,1,0! 重新啟動(dòng)求解器ASEL,S, , ,3,4,1! 選擇下臺(tái)階面ESLA

33、,R! 選擇下臺(tái)階土體單元Ekill,all! 對(duì)選擇的單元給予“死屬性”進(jìn)行噴錨支護(hù)，即激活噴射混凝土襯砌梁?jiǎn)卧湾^桿桿單元。Allsel! 選擇所有內(nèi)容ESEL,S,TYPE,2 ESEL,A,TYPE,3 Ealive,all上臺(tái)階土體開挖。Allsel! 選擇所有內(nèi)容ASEL,S, , ,1,2,1! 選擇下臺(tái)階面ESLA,R! 選擇下臺(tái)階土體單元Ekill,all! 對(duì)選擇的單元給予“死屬性”2 施加節(jié)點(diǎn)荷載 節(jié)點(diǎn)力計(jì)算。FINISH/POST1! 進(jìn)入后處理器Allsel! 選擇所有內(nèi)容ESEL,S, , ,! 人工手動(dòng)選擇開挖土體周圍的一圈單元。NFORCE,ALL! 對(duì)選擇的節(jié)點(diǎn)求解節(jié)點(diǎn)力FINISH/sol! 進(jìn)入求解理器節(jié)點(diǎn)力施加，加上節(jié)點(diǎn)力和位移邊界條件后的模型如圖1 下臺(tái)階開挖施工 下臺(tái)階土體開挖。圈單兀。5-32 所示。f,114,fx,6.46E+04f,115,fx,1.59E+05f,116,fx,0.00E+00f,117,fx,1.03E+05f,118,fx,8.16E+04f,119,fx,6.19E+04f,120,fx,4.06E+04f,121,fx,2.02E+04f,166,fx,-6.46E+04f,167,fx,-2.02E+04f,168,f