ansys在公路工程中的應(yīng)用實例

1、6.1 施工過程仿真分析1 材料、實常數(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 !進入前處理器ET,1,PLANE42 !設(shè)置實體單元類型KEYOPT,1,1,0KEYOPT,1,2,0K

2、EYOPT,1,3,2 !設(shè)置為平面應(yīng)變模式KEYOPT,1,5,0KEYOPT,1,6,0ET,2,PLANE42 !設(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 !輸入密度（圍巖）

3、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 !輸入

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

5、隧道輪廓線關(guā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 !作疊加布爾運算LOVLAP,12,4LDELE,6,10,2LDELE,7LDELE,13,16,3NUMMR

6、G,ALL, , , ,LOW !合并所有元素NUMCMP,ALL !壓縮所有編號FLST,2,3,4,ORDE,3 !作加運算FITEM,2,4 FITEM,2,5 FITEM,2,6 LCOMB,P51X, ,0 FLST,2,3,4,ORDE,3 !作加運算FITEM,2,3 FITEM,2,7 FITEM,2,9 LCOMB,P51X, ,0NUMMRG,ALL, , , ,LOW !合并所有元素NUMCMP,ALL !壓縮所有編號circle, 16,10.178 !畫仰拱圓LOVLAP,3,7 !作疊加布爾運算LOVLAP,4,8LDELE,19 !刪除部分線LDELE,5,7,1

7、LDELE,10,11,1LDELE,14,16,2FLST,2,4,4,ORDE,4 !做線加布爾運算FITEM,2,9 FITEM,2,12 FITEM,2,-13 FITEM,2,15 LCOMB,P51X, ,0 FLST,2,4,4,ORDE,4 !做線加布爾運算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,1LOVL

8、AP,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,8 FITEM,2,15 FITEM,2,-16 FITEM,2,19 FITEM,2,21 LCOMB,P51X, ,0 FLST,2,5,4,ORDE,

9、5 !做線加布爾運算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)鍵點，為畫其它直線作準備k,57,30.044,1.169,0

10、,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, 14 LSTR, 14, 19 LSTR

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

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

13、建面，依次創(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, 44, 11, Al, 44,

14、 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, 49, 25, 26, 62, Al,

15、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ù)。Save !按照“l(fā)esize,1,8,

16、1”格式，將其他61條線也進行被劃分單元數(shù)設(shè)置，其線的編號和被劃分的單元數(shù)如表6-3所示，劃分好后如圖6-10所示。表6-1 設(shè)置線單元數(shù)NL(NDIV)NL(NDIV)NL(NDIV)NL(NDIV)NL(NDIV)358398544554514368408584574324348388594604224378418614564524187843447418428108424534214981184544642043816844450419448128314301174581382642922446814825104854948815862102752810338234將所有面劃分單元，其單

17、元圖如圖6-11所示。AMAP,1,14,19,20, !劃分面1，依次單擊關(guān)鍵點14、19、20進行映射劃分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 !面

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

19、,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é)點創(chuàng)建梁單元E,11,12E,12,13E,13,14E,14,15E,15,16E,16,17E,17,10E,10,63

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

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

22、置將要創(chuàng)建單元的材料REAL, 1 !設(shè)置將要創(chuàng)建單元的幾何常數(shù)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 !保存數(shù)據(jù)E,

23、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 !返回上級目

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

25、!選擇所有內(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 !進入求解器Allsel !選擇所有內(nèi)容ANTYPE,REST,1,13,0 !重新啟動求解器ASEL,S, , ,

26、2 !選擇左洞左上部分土體面ESLA,R !選擇左洞左上部分土體單元Ekill,all !對選擇的單元給予“死屬性”Save 進行噴錨支護，即激活噴射混凝土襯砌梁單元和臨時支撐單元。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é)點ESLN,R,1 !選擇節(jié)點上的單元ESEL,R,TYPE,2 !選擇2類單元Ealive,all !對選擇的單元給予“生屬性”Save 節(jié)點力計算。FINISH /P

27、OST1 !進入后處理器Allsel !選擇所有內(nèi)容ESEL,S, , , !人工手動選擇開挖土體周圍的一圈單元圖6-16所示NFORCE,ALL !對選擇的節(jié)點求解節(jié)點力Save FINISH /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.10

28、41E+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+

29、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

30、,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方向位移PLNSOL,U,Y,0,1 ! Y方向位移PL

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

32、選擇左洞左上和左下部分土體單元Ekill,all !對選擇的單元給予“死屬性”Save 進行噴錨支護，即激活噴射混凝土襯砌梁單元和臨時支撐單元。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é)點ESLN,R,1 !選擇節(jié)點上的單元ESEL,R,TYPE,2 !選擇2類單元Ealive,all !對選擇的單元給予“生屬性”Save 節(jié)點力施加，加上節(jié)點力和位移邊界條件后的模型如圖6-23所示

33、。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.766

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

35、.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 主菜單SAVE計算下一步開挖計算所施加的節(jié)點力，如圖

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

37、元Ekill,all !對選擇的單元給予“死屬性”Allsel !選擇所有內(nèi)容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

38、,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,-0.1401E+06f,12,fx,-

39、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é)點力，如圖6-34所示。ESEL,S, , , !人工手動選擇開挖土體周圍