地鐵明挖隧道襯砌結構設計力學分析

1、襯砌結構力學有限元分析過程/TITLE,Mechanical analysis on cut and cover tunnel lining ! 確定分析標題/NOPR !菜單過濾設置/PMETH,OFF,0KEYW,PR_SET,1 KEYW,PR_STRUC,1 !保留結構分析部分菜單/COM, /COM,Preferences for GUI filtering have been set to display:/COM, Structural /PREP7 !進入前處理器ET,1,BEAM3 !設置梁單元類型ET,2, COMBIN14 !設置彈簧單元類型R,1,0.5,0.01041

2、67,0.5, , , , !設置梁單元幾何常數(shù)R,2,50e6, , , !設置彈簧單元幾何常數(shù)MPTEMP, !設置材料模型MPTEMP,1,0 MPDATA,EX,1,30e9 !輸入彈性模量MPDATA,PRXY,1,0.2 !輸入泊松比MPTEMP, !設置材料模型MPTEMP,1,0 MPDATA,DENS,1,2500 !輸入密度SAVE !保存數(shù)據(jù)庫K,1, !創(chuàng)建關鍵點K,2 ,4.7, K,3 ,4.7,5, K, 4,0,5, LSTR, 1, 2 !創(chuàng)建直線LSTR, 2, 3 LSTR, 3, 4 LSTR, 4, 1 SAVE !保存數(shù)據(jù)，生成的直線如圖4-7所示。

3、LESIZE,ALL,0.5, , , ,1, , ,1, !設置單元大小，本次劃分成0.5m長TYPE, 1 !設置將要創(chuàng)建單元的類型MAT, 1 !設置將要創(chuàng)建單元的材料REAL, 1 !設置將要創(chuàng)建單元的幾何常數(shù)FLST,2,4,4,ORDE,2 FITEM,2,1 FITEM,2,-4LMESH,P51X !將所有直線劃分單元/PNUM,KP,0 !以下為顯示單元編號和顏色/PNUM,ELEM,1/REPLOT !重新顯示SAVEFLST,4,9,1,ORDE,2 !選擇要復制的節(jié)點數(shù)量和起始點FITEM,4,23 FITEM,4,-31 NGEN,2,100,P51X, , , ,1

4、, ,1, !復制生成新的節(jié)點，Y坐標增加1，節(jié)點號增加100FLST,4,9,1,ORDE,2 FITEM,4,3 FITEM,4,-11 NGEN,2,100,P51X, , , ,-1, ,1, FLST,4,9,1,ORDE,2 FITEM,4,13 FITEM,4,-21 NGEN,2,100,P51X, , ,1, , ,1, FLST,4,9,1,ORDE,2 FITEM,4,32 FITEM,4,-40 NGEN,2,100,P51X, , ,-1, , ,1, FLST,4,1,1,ORDE,1 FITEM,4,12 NGEN,2,140,P51X, , ,0.71,0.71

5、, ,1, FLST,4,1,1,ORDE,1 FITEM,4,2 NGEN,2,152,P51X, , ,0.71,-0.71, ,1, FLST,4,1,1,ORDE,1 FITEM,4,22 NGEN,2,154,P51X, , ,-0.71,0.71, ,1, FLST,4,1,1,ORDE,1 FITEM,4,1 NGEN,2,176,P51X, , ,-0.71,-0.71, ,1, SAVETYPE, 2 !設置將要創(chuàng)建單元的類型MAT, 1 !設置將要創(chuàng)建單元的材料REAL, 2 !設置將要創(chuàng)建單元的幾何常數(shù)FLST,2,2,1 !選擇要創(chuàng)建單元的兩個節(jié)點FITEM,2,12

6、FITEM,2,152 E,P51X !創(chuàng)建彈簧單元FLST,2,2,1 FITEM,2,23 FITEM,2,123 E,P51X !依次繼續(xù)直到所有的彈簧單元創(chuàng)建完成FLST,2,2,1 FITEM,2,2 FITEM,2,154 E,P51X Finish !返回Main Menu 主菜單/SOL !進入求解器FLST,2,40,1,ORDE,10 !選擇要施加約束的節(jié)點FITEM,2,103 FITEM,2,-111FITEM,2,113 FITEM,2,-121FITEM,2,123 FITEM,2,-140FITEM,2,152 FITEM,2,154 FITEM,2,176 FI

7、TEM,2,-177D,P51X, , , , , ,UX,UY, , , , !在“Ux”和“Uy”兩個方向的施加約束ACEL,0,10,0, !在Y方向施加重力加速度F,1,FY,41125 !在節(jié)點上施加Y方向集中力F,2,FY,41125 !在底邊上施加力F,3,FY,82250F,4,FY,82250F,5,FY,82250F,6,FY,82250F,7,FY,82250F,8,FY,82250F,9,FY,82250F,10,FY,82250F,11,FY,82250F,12,FY,-32700 !在頂邊上施加力F,22,FY,-32700F,23,FY,-65400F,24,FY

8、,-65400F,25,FY,-65400F,26,FY,-65400F,27,FY,-65400F,28,FY,-65400F,29,FY,-65400F,30,FY,-65400F,31,FY,-65400F,2,FX,-32700 !在右邊上施加力F,13,FX,-62400 !在X方向施加力F,14,FX,-59400F,15,FX,-56400F,16,FX,-53400F,17,FX,-50400F,18,FX,-47400F,19,FX,-44400F,20,FX,-41400F,21,FX,-38400F,12,FX,-17700F,22,FX,17700 !在左邊上施加力F,

9、32,FX,38400F,33,FX,41400F,34,FX,44400F,35,FX,47400F,36,FX,50400F,37,FX,53400F,38,FX,56400F,39,FX,59400F,40,FX,62400F,41,FX,32700NROPT,FULL, , !采用全牛頓拉普森法進行求解Allsel !選擇所有內容Outres,all,all !輸出所有內容Solve !求解計算Finish !求解結束返回Main Menu 主菜單SAVE/POST1 !進入后處理器PLDISP,1 !繪制變形和未變形圖ETABLE, ,SMISC, 6 ! 6、12表示彎矩ETABL

10、E, ,SMISC, 12 ETABLE, ,SMISC, 1 ! 1、7表示軸力ETABLE, ,SMISC, 7 ETABLE, ,SMISC, 2 ! 2、8表示剪力ETABLE, ,SMISC, 8 查看內力，PLLS,SMIS6,SMIS12,-1,0 !繪制彎矩圖ESEL,U,TYPE,2 !僅顯示單元類型2PLLS,SMIS1,SMIS7,1,0 !繪制軸力圖PLLS,SMIS2,SMIS8,1,0 !繪制剪力圖Finish !結束后處理器操作/sol !進入求解器FLST,2,36,2 !選擇單元FITEM,2,50 FITEM,2,49 FITEM,2,48 FITEM,2,

11、47 FITEM,2,46 FITEM,2,45 FITEM,2,44 FITEM,2,43 FITEM,2,42 FITEM,2,80 FITEM,2,79 FITEM,2,78 FITEM,2,77 FITEM,2,76 FITEM,2,75 FITEM,2,74 FITEM,2,73 FITEM,2,72 FITEM,2,70 FITEM,2,69 FITEM,2,68 FITEM,2,67 FITEM,2,66 FITEM,2,65 FITEM,2,64 FITEM,2,63 FITEM,2,62 FITEM,2,60 FITEM,2,59 FITEM,2,58 FITEM,2,57 FITEM,2,56 FITEM,2,55 FITEM,2,54 FITEM,2,53 FITEM,2,52 EKILL,P5