ansys路基沉陷分析

ANSYS路基沉陷的彈塑性分析（命令流）問題描述：如圖所示路基結(jié)構(gòu)簡圖,各層材料均為Drucker-Prager彈塑性材料,但材料參數(shù)不相同。第一層截面為梯形,上下底邊長分別為14m和30m、高為6m,材料參數(shù)：材料彈性模量2.1GPa、泊松比0.25、密度1910kg/m3、粘聚力35kPa、內(nèi)摩擦角21°、膨脹角0°。第二層截面為矩形,長80m。寬8m,材料參數(shù)：彈性模量85MPa、泊松比0.30、密度1890kg/m3、粘聚力25kPa、內(nèi)摩擦角22°、膨脹角0°。第三層截面為矩形,長80m、寬12m,材料參數(shù)：彈性模量45MPa、泊松比0.30、密度1880kg/m3、粘聚力30kPa、內(nèi)摩擦角14°、膨脹角0°。第四層截面為矩形,長80m,寬6m,材料參數(shù)：彈性模量65MPa、泊松比0.30、密度2010kg/m3、粘聚力21kPa、內(nèi)摩擦角15°、膨脹角0°。！進入前處理模塊,定義平面應(yīng)變單元,定義各層的材料參數(shù)。/FILNAME,luji,1/PREP7ET,1,PLANE42 !定義平面單元PLANE42KEYOPT,1,3,2 !設(shè)置為平面應(yīng)變單元！第1層材料參數(shù)MP,EX,1,2.1E9 !第1層材料彈性模量2.1GPaMP,PRXY,1,0.25 !第1層材料泊松比0.25MP,DENS,1,1910 !第1層材料密度1910kg/m3TB,DP,1TBMODIF,1,1,35E3 !第1層材料粘聚力35kpaTBMODIF,1,2,21 !第1層材料摩擦角21度TBMODIF,1,3,0 !第1層材料膨脹角0度!第2層材料參數(shù)MP,EX,2,85E6 !第2層材料彈性模量85MPaMP,PRXY,2,0.30 !第2層材料泊松比0.30MP,DENS,2,1890 !第2層材料密度1890kg/m3TB,DP,2TBMODIF,1,1,25E3 !第2層材料粘聚力25kpaTBMODIF,1,2,22 !第2層材料摩擦角22度TBMODIF,1,3,0 !第2層材料膨脹角0度!第3層材料參數(shù)MP,EX,3,45E6 !第3層材料彈性模量45MPaMP,PRXY,3,0.30 !第3層材料泊松比0.30MP,DENS,3,1880 !第3層材料密度1880kg/m3TB,DP,3TBMODIF,1,1,30E3 !第3層材料粘聚力30kpaTBMODIF,1,2,14 !第3層材料摩擦角14度TBMODIF,1,3,0 !第3層材料膨脹角0度!第4層材料參數(shù)MP,EX,4,65E6 !第4層材料彈性模量65MPaMP,PRXY,4,0.30 !第4層材料泊松比0.30MP,DENS,4,2010 !第4層材料密度2010kg/m3TB,DP,4TBMODIF,1,1,21E3 !第4層材料粘聚力21kpaTBMODIF,1,2,15 !第4層材料摩擦角15度TBMODIF,1,3,0 !第4層材料膨脹角0度！先定義描述第1層梯形截面4個頂點的4個關(guān)鍵點,然后通過這4個關(guān)鍵點建立第1層的梯形截面,直接建立第2到第4層的矩形截面,將建立的4個截面粘結(jié)在1起,得到路基有限元模型。K,1,-15,0,0 !建立梯形頂點K,2,15,0,0 !建立梯形頂點K,3,7,5,0 !建立梯形頂點K,4,-7,5,0 !建立梯形頂點A,1,2,3,4 !建立第1層梯形截面RECTNG,-40,40,-8,0 !建立第2層矩形截面RECTNG,-40,40,-20,-8 !建立第3層矩形截面RECTNG,-40,40,-26,-20 !建立第4層矩形截面AGLUE,ALL !粘結(jié)所有面！設(shè)置總體單元尺寸,逐層設(shè)置單元材料類型、單元劃分方式,并逐層劃分有限元網(wǎng)格,得到有限元模型。ESIZE,1 !設(shè)置總體單元尺寸MSHKEY,1 !設(shè)置映射劃分網(wǎng)格方式MAT,1 !設(shè)置第1層材料類型AMESH,1 !劃分第1層有限元網(wǎng)格MAT,2 !設(shè)置第2層材料類型MSHKEY,0 !設(shè)置自由劃分網(wǎng)格方式AMESH,7 !劃分第2層有限元網(wǎng)格MAT,3 !設(shè)置第3層材料類型MSHKEY,1 !設(shè)置映射劃分網(wǎng)格方式AMESH,5 !劃分第3層有限元網(wǎng)格MAT,4 !設(shè)置第4層材料類型MSHKEY,1 !設(shè)置映射劃分網(wǎng)格方式AMESH,6 !劃分第4層有限元網(wǎng)格！進入計算求解模塊,定義靜力學分析類型,設(shè)置幾何非線性選項,設(shè)置邊界條件、并施加重力加速度,得到施加載荷后的路基幾何實體模型。設(shè)置計算結(jié)果輸出選項,定義計算時間和載荷子步數(shù),進行計算求解。/SOL !進入計算求解模塊ANTYPE,0 !定義靜力學分析類型NLGEOM,1 !考慮幾何非線性LSEL,S,LOC,X,-40 !選擇左邊界LSEL,A,LOC,X,40 !選擇右邊界DL,ALL,,UX !約束左右兩側(cè)水平位移LSEL,ALL !選擇所有線LSEL,S,LOC,Y,-26 !選擇底邊界DL,ALL,ALL !約束底邊界所有位移LSEL,ALL !選擇所有線ACEL,0,9.8 !施加重力加速度LPLOT !繪制線OUTRES,ALL,ALL !輸出所有計算結(jié)果TIME,1 !設(shè)置計算時間AUTOTS,1 !設(shè)置自動時間和載荷步選項NSUBST,100 !定義求解子步數(shù)KBC,0 !定義斜坡加載方式SOLVE