/COM, Structural ! 指定結構分析
/TITLE, Slope stability Analysis ! 定義工作標題
/FILNAM,support,1 ! 定義工作文件名
!進入前處理器
/PREP7
!定義單元類型
ET,1,PLANE82 !定義邊坡圍巖單元
KEYOPT,1,3,2 !求解類型為平面應變
!定義材料屬性
!邊坡圍巖1材料屬性
MP,EX,1,3E10
MP,PRXY,1,0.25
MP,DENS,1,2500 !定義塑性模型
TB,DP,1 !定義數據表1
TBDATA,1,9E5,42.8 !給數據表輸入數據
!邊坡圍巖2材料屬性 !采用彈性模型
MP,EX,2,3.2E10 !定義彈性模量
MP,PRXY,2,0.24 !定義泊松比
MP,DENS,2,2700 !定義密度
!取強度折減系數F=1.2時邊坡圍巖材料屬性
MP,EX,3,3E10
MP,PRXY,3,0.25
MP,DENS,3,2500
TB,DP,3 !定義數據表3
TBDATA,1,7.5E5,37.7 !給數據表輸入數據
!取強度折減系數F=1.4時邊坡圍巖材料屬性
MP,EX,4,3E10
MP,PRXY,4,0.25 !定義泊松比
MP,DENS,4,2500
TB,DP,4 !定義數據表4
TBDATA,1,6.4E5,33.5 !給數據表輸入數據
!取強度折減系數F=1.6時邊坡圍巖材料屬性
MP,EX,5,3E10
MP,PRXY,5,0.25
MP,DENS,5,2500 !定義密度
TB,DP,5 !定義數據表5
TBDATA,1,5.6E5,30 !給數據表輸入數據
!取強度折減系數F=1.8時邊坡圍巖材料屬性
MP,EX,6,3E10 !定義彈性模量
MP,PRXY,6,0.25
MP,DENS,6,2500
TB,DP,6 !定義數據表6
TBDATA,1,5.0E5,27.2
!取強度折減系數F=2.0時邊坡圍巖材料屬性
MP,EX,7,3E10
MP,PRXY,7,0.25 !定義泊松比
MP,DENS,7,2500
TB,DP,7 !定義數據表7
TBDATA,1,4.5E5,24.8
!取強度折減系數F=2.2時邊坡圍巖材料屬性
MP,EX,8,3E10
MP,PRXY,8,0.25
MP,DENS,8,2500
TB,DP,8 !定義數據表8
TBDATA,1,4.09E5,22.8
!取強度折減系數F=2.4時邊坡圍巖材料屬性
MP,EX,9,3E10
MP,PRXY,9,0.25
MP,DENS,9,2500
TB,DP,9 !定義數據表9
TBDATA,1,3.6E5,21.1 !給數據表輸入數據
!取強度折減系數F=2.6時邊坡圍巖材料屬性
MP,EX,10,3E10 !定義彈性模量
MP,PRXY,10,0.25 !定義泊松比
MP,DENS,10,2500 !定義密度
TB,DP,10 !定義數據表10
TBDATA,1,3.46E5,19.6 !給數據表10輸入數據
!取強度折減系數F=2.8時邊坡圍巖材料屬性
MP,EX,11,3E10 !定義彈性模量
MP,PRXY,11,0.25 !定義泊松比
MP,DENS,11,2500
TB,DP,11 !定義數據表11
TBDATA,1,3.2E5,18.2 !給數據表11輸入數據
!取強度折減系數F=3.0時邊坡圍巖材料屬性
MP,EX,12,3E10 !定義彈性模量
MP,PRXY,12,0.25
MP,DENS,12,2500 !定義密度
TB,DP,12 !定義數據表12
TBDATA,1,3.0E5,17.1 !給數據表12輸入數據
!建立幾何模型
!創建邊坡線模型線
K,1,,,, ! 創建關鍵點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 !連接關鍵點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
!創建邊坡面模型
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 !打開面號開關
!劃分網格生成有限元模型
! 劃分邊坡圍巖2網格
!設置網格份數
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 !采用單元類型1
mshkey,0 !設定自由網格劃分
mshape,0 !設定四邊形網格劃分
amesh,1 !劃分面積A1
! 劃分邊坡圍巖1網格
!設置網格份數
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 !設定自由網格劃分
mshape,0 !設定四邊形網格劃分
amesh,2 !劃分面積A2
amesh,3 !劃分面積A3
allsel
!保存網格模型
save,Slope-grid.db
!給邊坡模型施加約束
!邊坡兩側施加X方向約束
nsel,s,loc,x,-800 !選擇X=-800線上所有節點
nsel,a,loc,x,1200 !選擇X=1200線上所有節點
d,all,ux !對所選擇節點約束X方向位移
allsel
!邊坡底部施加約束
nsel,s,loc,y,-1200 !選擇Y=--1200線上所有節點
d,all,uy !對所選擇節點約束X、Y方向位移
d,all,ux
!施加重力加速度
acel,,9.8
/solu
!求解設置
antype,static !設定為靜力求解
nsubst,100 !設定最大子步數為100
pred,on !打開時間步長預測器
nropt,full !設定牛頓-拉普森選項
nlgeom,on !打開大位移效果
lnsrch,on !打開線性搜索
outres,all,all !輸出所有項
cnvtol,f,,0.005,2,0.5 !力收斂準則設定
cnvtol,u,,0.05,2,1 !位移收斂準則設定
!邊坡在強度折減系數F=1時求解
allsel
solve !進行求解
save,F1,db !把F=1時求解結果保存
!邊坡在強度折減系數F=1.2時求解
finish
/solu
allsel
asel,s,area,,2,3,1 !選擇面積A2、A3
MPCHG,3,all ! 把所選擇單元材料號改為3
allsel
solve !進行求解
save,F1.2,db !把F=1.2時求解結果保存
!邊坡在強度折減系數F=1.4時求解
finish
/solu
allsel
asel,s,area,,2,3,1 !選擇面積A2、A3
MPCHG,4,all !把所選擇單元材料號改為4
allsel
solve !進行求解
save,F1.4,db !把F=1.4時求解結果保存
!邊坡在強度折減系數F=1.6時求解
finish
/solu
allsel
asel,s,area,,2,3,1 !選擇面積A2、A3
MPCHG,5,all !把所選擇單元材料號改為5
allsel
solve
save,F1.6,db !把F=1.6時求解結果保存
!邊坡在強度折減系數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時求解結果保存
!邊坡在強度折減系數F=2.0時求解
finish
/solu
allsel
asel,s,area,,2,3,1 !選擇面積A2、A3
MPCHG,7,all !把所選擇單元材料號改為7
allsel
solve
save,F2.0,db !把F=2.0時求解結果保存
!邊坡在強度折減系數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時求解結果保存
!邊坡在強度折減系數F=2.4時求解
finish
/solu
allsel
asel,s,area,,2,3,1 !選擇面積A2、A3
MPCHG,9,all !把所選擇單元材料號改為9
allsel
solve
save,F2.4,db !把F=2.4時求解結果保存
!邊坡在強度折減系數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時求解結果保存
!邊坡在強度折減系數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.8時求解結果保存
!邊坡在強度折減系數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時求解結果保存
/post1 !進入后處理
!邊坡在強度折減系數F=1時結果分析
Resume,'F1','db' !讀入邊坡在強度折減系數F=1時
set,1,last !讀入后一個子步
pldisp,1 !繪制邊坡模型變形圖
plnsol,u,x !繪制邊坡模型水平方向位移云圖
plnsol,eppl,eqv !繪制邊坡模型塑性應變云圖
!邊坡在強度折減系數F=1.2時結果分析
Resume,'F1.2','db' !讀入邊坡在強度折減系數F=1.2時
set,1,last !讀入后一個子步
pldisp,1 !繪制邊坡模型變形圖
plnsol,u,x !繪制邊坡模型水平方向位移云圖
plnsol,eppl,eqv !繪制邊坡模型塑性應變云圖
!邊坡在強度折減系數F=1.4時結果分析
Resume,'F1.4','db' !讀入邊坡在強度折減系數F=1.4時
set,1,last !讀入后一個子步
pldisp,1 !繪制邊坡模型變形圖
plnsol,u,x !繪制邊坡模型水平方向位移云圖
plnsol,eppl,eqv !繪制邊坡模型塑性應變云圖
!邊坡在強度折減系數F=1.6時結果分析
Resume,'F1.6','db' !讀入邊坡在強度折減系數F=1.6時
set,1,last !讀入后一個子步
pldisp,1 !繪制邊坡模型變形圖
plnsol,u,x !繪制邊坡模型水平方向位移云圖
plnsol,eppl,eqv !繪制邊坡模型塑性應變云圖
!邊坡在強度折減系數F=1.8時結果分析
Resume,'F1.8','db' !讀入邊坡在強度折減系數F=1.8時
set,1,last !讀入后一個子步
pldisp,1 !繪制邊坡模型變形圖
plnsol,u,x !繪制邊坡模型水平方向位移云圖
plnsol,eppl,eqv !繪制邊坡模型塑性應變云圖
!邊坡在強度折減系數F=2.0時結果分析
Resume,'F2.0','db' !讀入邊坡在強度折減系數F=2.0時
set,1,last !讀入后一個子步
pldisp,1 !繪制邊坡模型變形圖
plnsol,u,x !繪制邊坡模型水平方向位移云圖
plnsol,eppl,eqv !繪制邊坡模型塑性應變云圖
!邊坡在強度折減系數F=2.2時結果分析
Resume,'F2.2','db' !讀入邊坡在強度折減系數F=2.2時
set,1,last !讀入后一個子步
pldisp,1 !繪制邊坡模型變形圖
plnsol,u,x !繪制邊坡模型水平方向位移云圖
plnsol,eppl,eqv !繪制邊坡模型塑性應變云圖
!邊坡在強度折減系數F=2.4時結果分析
Resume,'F2.4','db' !讀入邊坡在強度折減系數F=2.4時
set,1,last !讀入后一個子步
pldisp,1 !繪制邊坡模型變形圖
plnsol,u,x !繪制邊坡模型水平方向位移云圖
plnsol,eppl,eqv !繪制邊坡模型塑性應變云圖
!邊坡在強度折減系數F=2.6時結果分析
Resume,'F2.6','db' !讀入邊坡在強度折減系數F=2.6時
set,1,last !讀入后一個子步
pldisp,1 !繪制邊坡模型變形圖
plnsol,u,x !繪制邊坡模型水平方向位移云圖
plnsol,eppl,eqv !繪制邊坡模型塑性應變云圖
!邊坡在強度折減系數F=2.8時結果分析
Resume,'F2.8','db' !讀入邊坡在強度折減系數F=2.8時
set,1,last !讀入后一個子步
pldisp,1 !繪制邊坡模型變形圖
plnsol,u,x !繪制邊坡模型水平方向位移云圖
plnsol,eppl,eqv !繪制邊坡模型塑性應變云圖
!邊坡在強度折減系數F=3.0時結果分析
Resume,'F3.0','db' !讀入邊坡在強度折減系數F=3.0時
set,1,last !讀入后一個子步
pldisp,1 !繪制邊坡模型變形圖
plnsol,u,x !繪制邊坡模型水平方向位移云圖
plnsol,eppl,eqv !繪制邊坡模型塑性應變云圖