/CLEAR FINISH !定義工作名稱和分析標題 /FILNAME,Slope,1 /TITLE,Slope stability Analysis !定義單元類型 /prep7 ET,1,PLANE183 KEYOPT,1,3,2 !平面應變選項 !定義模型中的材料參數 MP,EX,1,31E9 !圍巖1彈性模量 MP,PRXY,1,0.24 !圍巖1泊松比 MP,DENS,1,2700 !圍巖1密度 *DO,I,2,13 XS=1.0+(I-2)*0.2 MP,EX,I,30E9 !圍巖2彈性模量 MP,PRXY,I,0.25 !圍巖2泊松比 MP,DENS,I,2500 !圍巖2密度 TB,DP,I TBMODIF,1,1,0.9E6/XS !圍巖2折減XS后的內聚力 TBMODIF,1,2,42/XS !圍巖2折減XS后的摩擦角 *ENDDO SAVE !存儲數據庫 !建立幾何模型并劃分單元 K,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 SAVE *DO,I,1,8 !由關鍵點創建線 L,I,I+1 *ENDDO L,9,1 L,1,7 L,3,6 SAVE AL,3,4,5,11 AL,1,2,11,6,10 AL,7,8,9,10 SAVE !劃分網格 LSEL,S,,,3,5,2 LESIZE,ALL,,,5 LSEL,S,,,4,11,7 LESIZE,ALL,,,26 MAT,1 TYPE,1 MSHAPE,0 !網格形狀為四邊形，1為三角形網格 MSHKEY,0 !網格劃分采用掃略方式；0為自由劃分；2為以掃略為主，自由為輔 AMESH,1 LSEL,S,,,2,6,4 LSEL,A,,,1 LESIZE,ALL,,,10 LSEL,S,,,8,10,2 LESIZE,ALL,,,16 LSEL,S,,,7,9,2 LESIZE,ALL,,,12 MAT,2 TYPE,1 MSHAPE,0 MSHKEY,0 AMESH,2 AMESH,3 ALLSEL SAVE,Slope_Grid.db !指定邊界條件 NSEL,S,LOC,X,-800 NSEL,A,LOC,X,1200 D,ALL,UX ALLSEL NSEL,S,LOC,Y,-1200 D,ALL,UX D,ALL,UY ACEL,9.8 SAVE !求解設置 /SOLU ANTYPE,STATIC !靜力求解 NSUBST,5,100,1 !荷載步5，最大荷載子步數100，最小荷載子步數1 PRED,ON !打開時間步長預測器 NROPT,FULL !設定牛頓-拉普森選項 NLGEOM,ON !打開大變形 LNSRCH,ON !打開線性搜索 OUTRES,ALL,ALL !輸出所有項 CNVTOL,F,,0.005,2,0.5 !力收斂準則 CNVTOL,U,,0.005,2,1 !位移收斂準則 !邊坡在強度折減系數XS=1.0時求解 ALLSEL SOLVE SAVE,F1.0,DB