ansys邊坡穩(wěn)定命令流

1、(1創(chuàng)建物理環(huán)境 /COM, Structural ! 指定結(jié)構(gòu)分析 /TITLE, Slope stability Analysis ! 定義工作標(biāo)題 /FILNAM,support,1 ! 定義工作文件名 (2建立模型 !進(jìn)入前處理器 /PREP7 !定義單元類型 ET,1,PLANE82 !定義邊坡圍巖單元 KEYOPT,1,3,2 !求解類型為平面應(yīng)變 !定義材料屬性 !邊坡圍巖1材料屬性 MP,EX,1,3E10 MP,PRXY,1,0.25 MP,DENS,1,2500 !定義塑性模型 TB,DP,1 !定義數(shù)據(jù)表1 TBDATA,1,9E5,42.8 !給數(shù)據(jù)表輸入數(shù)據(jù) !邊坡圍

2、巖2材料屬性 !采用彈性模型 MP,EX,2,3.2E10 ！定義彈性模量 MP,PRXY,2,0.24 !定義泊松比 MP,DENS,2,2700 !定義密度 !取強(qiáng)度折減系數(shù)F=1.2時邊坡圍巖材料屬性 MP,EX,3,3E10 MP,PRXY,3,0.25 MP,DENS,3,2500 TB,DP,3 !定義數(shù)據(jù)表3 TBDATA,1,7.5E5,37.7 !給數(shù)據(jù)表輸入數(shù)據(jù) !取強(qiáng)度折減系數(shù)F=1.4時邊坡圍巖材料屬性 MP,EX,4,3E10 MP,PRXY,4,0.25 !定義泊松比 MP,DENS,4,2500 TB,DP,4 !定義數(shù)據(jù)表4 TBDATA,1,6.4E5,33.

3、5 !給數(shù)據(jù)表輸入數(shù)據(jù) !取強(qiáng)度折減系數(shù)F=1.6時邊坡圍巖材料屬性 MP,EX,5,3E10 MP,PRXY,5,0.25 MP,DENS,5,2500 ！定義密度 TB,DP,5 !定義數(shù)據(jù)表5 TBDATA,1,5.6E5,30 !給數(shù)據(jù)表輸入數(shù)據(jù) !取強(qiáng)度折減系數(shù)F=1.8時邊坡圍巖材料屬性 MP,EX,6,3E10 ！定義彈性模量 MP,PRXY,6,0.25 MP,DENS,6,2500 TB,DP,6 !定義數(shù)據(jù)表6 TBDATA,1,5.0E5,27.2 !取強(qiáng)度折減系數(shù)F=2.0時邊坡圍巖材料屬性 MP,EX,7,3E10 MP,PRXY,7,0.25 !定義泊松比 MP,D

4、ENS,7,2500 TB,DP,7 !定義數(shù)據(jù)表7 TBDATA,1,4.5E5,24.8 !取強(qiáng)度折減系數(shù)F=2.2時邊坡圍巖材料屬性 MP,EX,8,3E10 MP,PRXY,8,0.25 MP,DENS,8,2500 TB,DP,8 !定義數(shù)據(jù)表8 TBDATA,1,4.09E5,22.8 !取強(qiáng)度折減系數(shù)F=2.4時邊坡圍巖材料屬性 MP,EX,9,3E10 MP,PRXY,9,0.25 MP,DENS,9,2500 TB,DP,9 !定義數(shù)據(jù)表9 TBDATA,1,3.6E5,21.1 !給數(shù)據(jù)表輸入數(shù)據(jù) !取強(qiáng)度折減系數(shù)F=2.6時邊坡圍巖材料屬性 MP,EX,10,3E10 ！

5、定義彈性模量 MP,PRXY,10,0.25 !定義泊松比 MP,DENS,10,2500 ！定義密度 TB,DP,10 !定義數(shù)據(jù)表10 TBDATA,1,3.46E5,19.6 !給數(shù)據(jù)表10輸入數(shù)據(jù) !取強(qiáng)度折減系數(shù)F=2.8時邊坡圍巖材料屬性 MP,EX,11,3E10 ！定義彈性模量 MP,PRXY,11,0.25 !定義泊松比 MP,DENS,11,2500 TB,DP,11 !定義數(shù)據(jù)表11 TBDATA,1,3.2E5,18.2 !給數(shù)據(jù)表11輸入數(shù)據(jù) !取強(qiáng)度折減系數(shù)F=3.0時邊坡圍巖材料屬性 MP,EX,12,3E10 ！定義彈性模量 MP,PRXY,12,0.25 MP

6、,DENS,12,2500 ！定義密度 TB,DP,12 !定義數(shù)據(jù)表12 TBDATA,1,3.0E5,17.1 !給數(shù)據(jù)表12輸入數(shù)據(jù) !建立幾何模型 !創(chuàng)建邊坡線模型線 K,1, ！ 創(chuàng)建關(guān)鍵點(diǎn)1 K,2,-800, K,3,-800,-800, K,4,-800,-1200, K,5,1200,-1200, K,6,1200,-800, K,7,1200, K,8,1200,378, K,9,430,378, l,1,2 !連接關(guān)鍵點(diǎn)1、2生成直線L1 l,2,3 l,3,4 l,4,5 l,5,6 l,6,7 l,7,8 l,8,9 l,9,1 l,1,7 l,3,6 !創(chuàng)建邊坡面模

7、型 al,3,4,5,11 !由線L3、L4、L5、L11生成一個面A1 al,1,2,11,6,10 !由線L1、L2、L1、L6、L10生成一個面A2 al,7,8,9,10 !由線L7、L8、L9、L10生成一個面A3 /PNUM, area, on !打開面號開關(guān) !劃分網(wǎng)格生成有限元模型 ! 劃分邊坡圍巖2網(wǎng)格 !設(shè)置網(wǎng)格份數(shù) lsel,s,3,5,2 !選擇線L3和L5 lesize,all,5 !把所選擇線分為5段 lsel,s,4,11,7 !選擇線L4和L11 lesize,all,26 !把所選擇線分為26段 mat,2 !給邊坡圍巖2賦予2號材料特性 type,1 !采用

8、單元類型1 mshkey,0 !設(shè)定自由網(wǎng)格劃分 mshape,0 !設(shè)定四邊形網(wǎng)格劃分 amesh,1 ！劃分面積A1 ! 劃分邊坡圍巖1網(wǎng)格 !設(shè)置網(wǎng)格份數(shù) lsel,s,2,6,4 !選擇線L2和L6 lsel,a,1 !附加選擇線L1 lesize,all,10 !把所選擇線分為10段 lsel,s,8,10,2 !選擇線L8和L10 lesize,all,16 !把所選擇線分為16段 lsel,s,7,9,2 lesize,all,12 mat,1 !給邊坡圍巖1賦予1號材料特性 type,1 mshkey,0 !設(shè)定自由網(wǎng)格劃分 mshape,0 !設(shè)定四邊形網(wǎng)格劃分 amesh,

9、2 !劃分面積A2 amesh,3 ！劃分 面積A3 allsel !保存網(wǎng)格模型 save,Slope-grid.db (3施加約束和荷載 !給邊坡模型施加約束 !邊坡兩側(cè)施加X方向約束 nsel,s,loc,x,-800 !選擇X=-800線上所有節(jié)點(diǎn) nsel,a,loc,x,1200 !選擇X=1200線上所有節(jié)點(diǎn) d,all,ux !對所選擇節(jié)點(diǎn)約束X方向位移 allsel !邊坡底部施加約束 nsel,s,loc,y,-1200 !選擇Y=-1200線上所有節(jié)點(diǎn) d,all,uy !對所選擇節(jié)點(diǎn)約束X、Y方向位移 d,all,ux !施加重力加速度 acel,9.8 (4求解 /s

10、olu !求解設(shè)置 antype,static !設(shè)定為靜力求解 nsubst,100 !設(shè)定最大子步數(shù)為100 pred,on !打開時間步長預(yù)測器 nropt,full !設(shè)定牛頓-拉普森選項(xiàng) nlgeom,on !打開大位移效果 lnsrch,on !打開線性搜索 outres,all,all !輸出所有項(xiàng) cnvtol,f,0.005,2,0.5 !力收斂準(zhǔn)則設(shè)定 cnvtol,u,0.05,2,1 !位移收斂準(zhǔn)則設(shè)定 !邊坡在強(qiáng)度折減系數(shù)F=1時求解 allsel solve !進(jìn)行求解 save,F1,db !把F=1時求解結(jié)果保存 !邊坡在強(qiáng)度折減系數(shù)F=1.2時求解 finis

11、h /solu allsel asel,s,area,2,3,1 !選擇面積A2、A3 MPCHG,3,all ! 把所選擇單元材料號改為3 allsel solve !進(jìn)行求解 save,F1.2,db !把F=1.2時求解結(jié)果保存 !邊坡在強(qiáng)度折減系數(shù)F=1.4時求解 finish /solu allsel asel,s,area,2,3,1 !選擇面積A2、A3 MPCHG,4,all !把所選擇單元材料號改為4 allsel solve !進(jìn)行求解 save,F1.4,db !把F=1.4時求解結(jié)果保存 !邊坡在強(qiáng)度折減系數(shù)F=1.6時求解 finish /solu allsel as

12、el,s,area,2,3,1 !選擇面積A2、A3 MPCHG,5,all !把所選擇單元材料號改為5 allsel solve save,F1.6,db !把F=1.6時求解結(jié)果保存 !邊坡在強(qiáng)度折減系數(shù)F=1.8時求解 finish /solu allsel asel,s,area,2,3,1 !選擇面積A2、A3 MPCHG,6,all !把所選擇單元材料號改為 6 allsel solve save,F1.8,db !把F=1.8時求解結(jié)果保存 !邊坡在強(qiáng)度折減系數(shù)F=2.0時求解 finish /solu allsel asel,s,area,2,3,1 !選擇面積A2、A3 MP

13、CHG,7,all !把所選擇單元材料號改為7 allsel solve save,F2.0,db !把F=2.0時求解結(jié)果保存 !邊坡在強(qiáng)度折減系數(shù)F=2.2時求解 finish /solu allsel asel,s,area,2,3,1 !選擇面積A2、A3 MPCHG,8,all !把所選擇單元材料號改為8 allsel solve save,F2.2,db !把F=2.2時求解結(jié)果保存 !邊坡在強(qiáng)度折減系數(shù)F=2.4時求解 finish /solu allsel asel,s,area,2,3,1 !選擇面積A2、A3 MPCHG,9,all !把所選擇單元材料號改為9 allsel

14、 solve save,F2.4,db !把F=2.4時求解結(jié)果保存 !邊坡在強(qiáng)度折減系數(shù)F=2.6時求解 finish /solu allsel asel,s,area,2,3,1 !選擇面積A2、A3 MPCHG,10,all !把所選擇單元材料號改為10 allsel solve save,F2.6,db !把F=2.6時求解結(jié)果保存 !邊坡在強(qiáng)度折減系數(shù)F=2.8時求解 finish /solu allsel asel,s,area,2,3,1 !選擇面積A2、A3 MPCHG,11,all !把所選擇單元材料號改為11 allsel solve save,F2.8,db !把F=2.

15、8時求解結(jié)果保存 !邊坡在強(qiáng)度折減系數(shù)F=3.0時求解 finish /solu allsel asel,s,area,2,3,1 !選擇面積A2、A3 MPCHG,12,all !把所選擇單元材料號改為12 allsel solve save,F3.0,db !把F=3.0時求解結(jié)果保存 (5）后處理 /post1 ！進(jìn)入后處理 !邊坡在強(qiáng)度折減系數(shù)F=1時結(jié)果分析 Resume,'F1','db' !讀入邊坡在強(qiáng)度折減系數(shù)F=1時 set,1,last !讀入后一個子步 pldisp,1 !繪制邊坡模型變形圖 plnsol,u,x !繪制邊坡模型水平方向位移

16、云圖 plnsol,eppl,eqv !繪制邊坡模型塑性應(yīng)變云圖 !邊坡在強(qiáng)度折減系數(shù)F=1.2時結(jié)果分析 Resume,'F1.2','db' !讀入邊坡在強(qiáng)度折減系數(shù)F=1.2時 set,1,last !讀入后一個子步 pldisp,1 !繪制邊坡模型變形圖 plnsol,u,x !繪制邊坡模型水平方向位移云圖 plnsol,eppl,eqv !繪制邊坡模型塑性應(yīng)變云圖 !邊坡在強(qiáng)度折減系數(shù)F=1.4時結(jié)果分析 Resume,'F1.4','db' !讀入邊坡在強(qiáng)度折減系數(shù)F=1.4時 set,1,last !讀入后一個子步

17、pldisp,1 !繪制邊坡模型變形圖 plnsol,u,x !繪制邊坡模型水平方向位移云圖 plnsol,eppl,eqv !繪制邊坡模型塑性應(yīng)變云圖 !邊坡在強(qiáng)度折減系數(shù)F=1.6時結(jié)果分析 Resume,'F1.6','db' !讀入邊坡在強(qiáng)度折減系數(shù)F=1.6時 set,1,last !讀入后一個子步 pldisp,1 !繪制邊坡模型變形圖 plnsol,u,x !繪制邊坡模型水平方向位移云圖 plnsol,eppl,eqv !繪制邊坡模型塑性應(yīng)變云圖 !邊坡在強(qiáng)度折減系數(shù)F=1.8時結(jié)果分析 Resume,'F1.8','db&

18、#39; !讀入邊坡在強(qiáng)度折減系數(shù)F=1.8時 set,1,last !讀入后一個子步 pldisp,1 !繪制邊坡模型變形圖 plnsol,u,x !繪制邊坡模型水平方向位移云圖 plnsol,eppl,eqv !繪制邊坡模型塑性應(yīng)變云圖 !邊坡在強(qiáng)度折減系數(shù)F=2.0時結(jié)果分析 Resume,'F2.0','db' !讀入邊坡在強(qiáng)度折減系數(shù)F=2.0時 set,1,last !讀入后一個子步 pldisp,1 !繪制邊坡模型變形圖 plnsol,u,x !繪制邊坡模型水平方向位移云圖 plnsol,eppl,eqv !繪制邊坡模型塑性應(yīng)變云圖 !邊坡在強(qiáng)度折減系數(shù)F=2.2時結(jié)果分析 Resume,'F2.2','db' !讀入邊坡在強(qiáng)度折減系數(shù)F=2.2時 set,1,last !讀入后一個子步 pldisp,1 !繪制邊坡模型變形圖 plnsol,u,x !繪制邊坡模型水平方向位移云圖 plnsol,eppl,eqv !繪制邊坡模型塑性應(yīng)變云圖 !邊坡在強(qiáng)度折減系數(shù)F=2.4時結(jié)果分析 Resume,'F2.4','db' !讀入邊坡在強(qiáng)