【正文】 !定義數(shù)據(jù)表7TBDATA,1,!取強(qiáng)度折減系數(shù)F=MP,EX,8,3E10 MP,PRXY,8, MP,DENS,8,2500TB,DP,8 !定義數(shù)據(jù)表8TBDATA,1,!取強(qiáng)度折減系數(shù)F=MP,EX,9,3E10 MP,PRXY,9, MP,DENS,9,2500TB,DP,9 !定義數(shù)據(jù)表9TBDATA,1, !給數(shù)據(jù)表輸入數(shù)據(jù)!取強(qiáng)度折減系數(shù)F=MP,EX,10,3E10 ！定義彈性模量 MP,PRXY,10, !定義泊松比MP,DENS,10,2500 ！定義密度TB,DP,10 !定義數(shù)據(jù)表10TBDATA,1, !給數(shù)據(jù)表10輸入數(shù)據(jù)!取強(qiáng)度折減系數(shù)F=MP,EX,11,3E10 ！定義彈性模量MP,PRXY,11, !定義泊松比MP,DENS,11,2500TB,DP,11 !定義數(shù)據(jù)表11TBDATA,1, !給數(shù)據(jù)表11輸入數(shù)據(jù)!取強(qiáng)度折減系數(shù)F=MP,EX,12,3E10 ！定義彈性模量MP,PRXY,12, MP,DENS,12,2500 ！定義密度TB,DP,12 !定義數(shù)據(jù)表12TBDATA,1, !給數(shù)據(jù)表12輸入數(shù)據(jù)!建立幾何模型!創(chuàng)建邊坡線模型K,1, ！ 創(chuàng)建關(guān)鍵點(diǎn)1K,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)2生成直線L1l,2,3l,3,4l,4,5l,5,6l,6,7l,7,8l,8,9l,9,1l,1,7l,3,6!創(chuàng)建邊坡面模型al,3,4,5,11 !由線LLLL11生成一個(gè)面A1al,1,2,11,6,10 !由線LLLLL10生成一個(gè)面A2al,7,8,9,10 !由線LLLL10生成一個(gè)面A3/PNUM, area, on !打開面號(hào)開關(guān)!劃分網(wǎng)格生成有限元模型! 劃分邊坡圍巖2網(wǎng)格!設(shè)置網(wǎng)格份數(shù)lsel,s,3,5,2 !選擇線L3和L5lesize,all,5 !把所選擇線分為5段lsel,s,4,11,7 !選擇線L4和L11lesize,all,26 !把所選擇線分為26段mat,2 !給邊坡圍巖2賦予2號(hào)材料特性type,1 !采用單元類型1mshkey,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和L6lsel,a,1 !附加選擇線L1lesize,all,10 !把所選擇線分為10段lsel,s,8,10,2 !選擇線L8和L10lesize,all,16 !把所選擇線分為16段lsel,s,7,9,2lesize,all,12mat,1 !給邊坡圍巖1賦予1號(hào)材料特性type,1mshkey,0 !設(shè)定自由網(wǎng)格劃分mshape,0 !設(shè)定四邊形網(wǎng)格劃分amesh,2 !劃分面積A2 amesh,3 ！劃分面積A3allsel!保存網(wǎng)格模型save, (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 !對(duì)所選擇節(jié)點(diǎn)約束X方向位移allsel!邊坡底部施加約束nsel,s,loc,y,1200 !選擇Y=1200線上所有節(jié)點(diǎn)d,all,uy !對(duì)所選擇節(jié)點(diǎn)約束X、Y方向位移d,all,ux !施加重力加速度acel,(4)求解/solu!求解設(shè)置antype,static !設(shè)定為靜力求解nsubst,100 !設(shè)定最大子步數(shù)為100pred,on !打開時(shí)間步長預(yù)測(cè)器nropt,full !設(shè)定牛頓拉普森選項(xiàng)nlgeom,on !打開大位移效果lnsrch,on !打開線性搜索outres,all,all !輸出所有項(xiàng)vtol,f,2, !力收斂準(zhǔn)則設(shè)定vtol,u,2,1 !位移收斂準(zhǔn)則設(shè)定!邊坡在強(qiáng)度折減系數(shù)F=1時(shí)求解allselsolve !進(jìn)行求解save,F1,db !把F=1時(shí)求解結(jié)果保存!邊坡在強(qiáng)度折減系數(shù)F=finish/soluallselasel,s,area,2,3,1 !選擇面積AA3MPCHG,3,all ! 把所選擇單元材料號(hào)改為3allselsolve !進(jìn)行求解save,db !把F= !邊坡在強(qiáng)度折減系數(shù)F=finish/soluallselasel,s,area,2,3,1 !選擇面積AA3MPCHG,4,all !把所選擇單元材料號(hào)改為4allselsolve !進(jìn)行求解save,db !把F= !邊坡在強(qiáng)度折減系數(shù)F=finish/soluallselasel,s,area,2,3,1 !選擇面積AA3MPCHG,5,all !把所選擇單元材料號(hào)改為5allselsolvesave,db !把F= !邊坡在強(qiáng)度折減系數(shù)F=finish/soluallselasel,s,area,2,3,1 !選擇面積AA3MPCHG,6,all !把所選擇單元材料號(hào)改為6allselsolvesave,db !把F= !邊坡在強(qiáng)度折減系數(shù)F=finish/soluallselasel,s,area,2,3,1 !選擇面積AA3MPCHG,7,all !把所選擇單元材料號(hào)改為7allselsolvesave,db !把F= !邊坡在強(qiáng)度折減系數(shù)F=finish/soluallselasel,s,area,2,3,1 !選擇面積AA3MPCHG,8,all !把所選擇單元材料號(hào)改為8allselsolvesave,db !把F= !邊坡在強(qiáng)度折減系數(shù)F=finish/soluallselasel,s,area,2,3,1 !選擇面積AA3MPCHG,9,all !把所選擇單元材料號(hào)改為9allselsolvesave,db !把F= !邊坡在強(qiáng)度折減系數(shù)F=finish/soluallselasel,s,area,2,3,1 !選擇面積AA3MPCHG,