(1)创建物理环境 ^$Dpdz I
/COM, Structural ;=Jj{FoG%
/TITLE, Slope stability Analysis 2:Q9 g ru
/FILNAM,support,1 M;={] w@ n
(2)建立模型 Zi[ @xG8dm
!进入前处理器 `9Zoq=/
/PREP7 a0C f.[L
!定义单元类型 0ejx; Mum
ET,1,PLANE82 !定义边坡围岩单元 /Ws@YP
KEYOPT,1,3,2 !求解类型为平面应变 y iE[^2Pv
!定义材料属性 bN@V=C 3
!边坡围岩1材料属性 4r `u@
MP,EX,1,3E10 sXkWs2!
MP,PRXY,1,0.25 f*7/O |Gp
MP,DENS,1,2500 !定义塑性模型 xxgS! J
TB,DP,1 !定义数据表1 f@+[-y F
TBDATA,1,9E5,42.8 !给数据表输入数据 Xz&Hfs"/J
!边坡围岩2材料属性 !采用弹性模型 dU\,>3tG
MP,EX,2,3.2E10 ！定义弹性模量 .{LFc|Z[
MP,PRXY,2,0.24 !定义泊松比 /cUcfe#X
MP,DENS,2,2700 !定义密度 :f?\ mVS+
!取强度折减系数F=1.2时边坡围岩材料属性 C] q Y
MP,EX,3,3E10 z _~f/
MP,PRXY,3,0.25 UMuqdLaT9
MP,DENS,3,2500 ~cVFCM
TB,DP,3 !定义数据表3 2j` x^
TBDATA,1,7.5E5,37.7 !给数据表输入数据 k5]s~* ,0
!取强度折减系数F=1.4时边坡围岩材料属性 Mb=vIk{B f
MP,EX,4,3E10 [WnX'R R
MP,PRXY,4,0.25 !定义泊松比 <Vm+Lt9
MP,DENS,4,2500 V03U"eI="
TB,DP,4 !定义数据表4 |\%[e@u
TBDATA,1,6.4E5,33.5 !给数据表输入数据 A9ia[2[
!取强度折减系数F=1.6时边坡围岩材料属性 e3UGYwQ
MP,EX,5,3E10 t&^9o $
MP,PRXY,5,0.25 >xws
MP,DENS,5,2500 ！定义密度 -z./6dQ
TB,DP,5 !定义数据表5 Lnbbv
*
TBDATA,1,5.6E5,30 !给数据表输入数据 1w/1 k6`0
!取强度折减系数F=1.8时边坡围岩材料属性 ifl`QZp_
MP,EX,6,3E10 ！定义弹性模量 RrRrB"!8nR
MP,PRXY,6,0.25 & WE| 9
MP,DENS,6,2500 -gefdx6ES
TB,DP,6 !定义数据表6 $]K gs6=r
TBDATA,1,5.0E5,27.2 k{(R.gLZG
!取强度折减系数F=2.0时边坡围岩材料属性 N}b/; Y
MP,EX,7,3E10 YwyP+S r\
MP,PRXY,7,0.25 !定义泊松比 $Mm=5 K %
MP,DENS,7,2500 @ |v4 B[/
TB,DP,7 !定义数据表7 n]?Yv E
TBDATA,1,4.5E5,24.8 Ikql
!取强度折减系数F=2.2时边坡围岩材料属性 bO>q`%&
MP,EX,8,3E10 #\[h.4i
MP,PRXY,8,0.25 ,3 &XV %1
MP,DENS,8,2500 j}3Avu%
TB,DP,8 !定义数据表8 i,1=5@rw5
TBDATA,1,4.09E5,22.8 ] /hF!eO
!取强度折减系数F=2.4时边坡围岩材料属性 k'E3{8<!
MP,EX,9,3E10 R7}=k)U?d@
MP,PRXY,9,0.25 A94ZG:
MP,DENS,9,2500 nfET;:{
TB,DP,9 !定义数据表9 htX;"R&
TBDATA,1,3.6E5,21.1 !给数据表输入数据 |1wfLJ4--l
!取强度折减系数F=2.6时边坡围岩材料属性 c[J(H,m t/
MP,EX,10,3E10 ！定义弹性模量 16ahU$@-
MP,PRXY,10,0.25 !定义泊松比 =/Pmi_
MP,DENS,10,2500 ！定义密度 X ptb4]
TB,DP,10 !定义数据表10 M3 ihtY
TBDATA,1,3.46E5,19.6 !给数据表10输入数据 A5cx !h
!取强度折减系数F=2.8时边坡围岩材料属性 5[}3j1
MP,EX,11,3E10 ！定义弹性模量 EjW3_ %
MP,PRXY,11,0.25 !定义泊松比 S7]\tw_L)
MP,DENS,11,2500 /~[+'
TB,DP,11 !定义数据表11 ;[TC`DuNj0
TBDATA,1,3.2E5,18.2 !给数据表11输入数据 ivDmPHj{
!取强度折减系数F=3.0时边坡围岩材料属性 "S)2<tV
MP,EX,12,3E10 ！定义弹性模量 nf=*KS\v
MP,PRXY,12,0.25 =&G<^7
MP,DENS,12,2500 ！定义密度 oU`8\ n](
TB,DP,12 !定义数据表12 m}&cX Y
TBDATA,1,3.0E5,17.1 !给数据表12输入数据 a,g3 /
!建立几何模型 f.8Jp<S2K
!创建边坡线模型线 e^2e[rp0
K,1,,,, ！ 创建关键点1 2u9O+]EP
K,2,-800,,, CZ* #FY
K,3,-800,-800,, @v{lH&K:;
K,4,-800,-1200,, )@+lfIE(l
K,5,1200,-1200,, IlQNo 1
K,6,1200,-800,, DT *'r;
K,7,1200,,, %;ZWYj`]n
K,8,1200,378,, ^f*}]`S
K,9,430,378,, h!.#r*vV
l,1,2 !连接关键点1、2生成直线L1 ^*JpdmVhu
l,2,3 vM )2F
l,3,4 !lQGoXQ'4
l,4,5 /Iht,@%E
l,5,6 8axz`2 `
l,6,7 V@krw" vW
l,7,8 ETU.v*HT]
l,8,9 uG<VQ2LM
l,9,1 Tk~RT<\Ab+
l,1,7 5@l[!Jl0k
l,3,6 y.,S}7l:
!创建边坡面模型 lO=+V 6
al,3,4,5,11 !由线L3、L4、L5、L11生成一个面A1 HQPb
al,1,2,11,6,10 !由线L1、L2、L1、L6、L10生成一个面A2 rHe*/nN%*
al,7,8,9,10 !由线L7、L8、L9、L10生成一个面A3 M]9oSi
/PNUM, area, on !打开面号开关 #m>Rt~(,S
!划分网格生成有限元模型 T}')QC&wQ
! 划分边坡围岩2网格 6G2s^P1Dl@
!设置网格份数 4HX qRFUD
lsel,s,,,3,5,2 !选择线L3和L5 j 7 r! N^
lesize,all,,,5 !把所选择线分为5段 a (Q4*XH4
lsel,s,,,4,11,7 !选择线L4和L11 T.ZPpxY
lesize,all,,,26 !把所选择线分为26段 >! u@>
mat,2 !给边坡围岩2赋予2号材料特性 vwy10PlqL
type,1 !采用单元类型1 r*$$82s
mshkey,0 !设定自由网格划分 },1**_#<Br
mshape,0 !设定四边形网格划分 i>=d7'oR
amesh,1 ！划分面积A1 : (IP rQ
! 划分边坡围岩1网格 _<Hx1l~
!设置网格份数 Twqkd8[
lsel,s,,,2,6,4 !选择线L2和L6 zx)^!dEMM
lsel,a,,,1 !附加选择线L1 Ve/"9 ?Y_
lesize,all,,,10 !把所选择线分为10段 j5hM |\]
lsel,s,,,8,10,2 !选择线L8和L10 QG{).|pm
lesize,all,,,16 !把所选择线分为16段 J6m`XC
lsel,s,,,7,9,2 nT7{`aaQl
lesize,all,,,12 <eWGvIEP[
mat,1 !给边坡围岩1赋予1号材料特性 Xxl>,QUA
type,1 VrG 4wLpLs
mshkey,0 !设定自由网格划分 1X-K u GaD
mshape,0 !设定四边形网格划分 P
"S=RX#+
amesh,2 !划分面积A2 /48W]a}JS
amesh,3 ！划分面积A3 \Xy]z
allsel {8L)F w
!保存网格模型 #$ \ cRLPg
save,Slope-grid.db >Hh8K<@NL
(3)施加约束和荷载 >&uG1q0p.
!给边坡模型施加约束 R X]x3 -
!边坡两侧施加X方向约束 oECM1'=Bf
nsel,s,loc,x,-800 !选择X=-800线上所有节点 Ub1? dk
nsel,a,loc,x,1200 !选择X=1200线上所有节点 X`,4pSQ;
d,all,ux !对所选择节点约束X方向位移 !A R$JUn X
allsel wLa^pI4p ^
!边坡底部施加约束 sU7>q }!
nsel,s,loc,y,-1200 !选择Y=--1200线上所有节点 [,GXA)j
d,all,uy !对所选择节点约束X、Y方向位移 F[aow$",+}
d,all,ux ?(8z O"
!施加重力加速度 {HgW9N(
acel,,9.8 o;#{N~4[$
(4)求解 x .>E 7
+
/solu IKJ~sw~AQ
!求解设置 bYzBe\^3q3
antype,static !设定为静力求解 GPqF>
nsubst,100 !设定最大子步数为100 doR'E= Z4h
pred,on !打开时间步长预测器 $Cu/!GA4.>
nropt,full !设定牛顿-拉普森选项 %}z/_QZ
nlgeom,on !打开大位移效果 E\M{/.4 4
lnsrch,on !打开线性搜索 >.R6\ >N%
outres,all,all !输出所有项 3YY<2<
cnvtol,f,,0.005,2,0.5 !力收敛准则设定 " mKMym2
cnvtol,u,,0.05,2,1 !位移收敛准则设定 >M` swEj
!边坡在强度折减系数F=1时求解 cQ4TYr;?
allsel
Q@3.0Hf|{
solve !进行求解 nMU#g])y)
save,F1,db !把F=1时求解结果保存 IB wqu w+
!边坡在强度折减系数F=1.2时求解 vFfvvRda4x
finish Z &mV1dxR
/solu niO(>
allsel -KRHcr \
asel,s,area,,2,3,1 !选择面积A2、A3 V/5hEo Dt
MPCHG,3,all ! 把所选择单元材料号改为3 nG#lr YZw
allsel ;qI5GQ {
solve !进行求解 Df_W>QC
save,F1.2,db !把F=1.2时求解结果保存 0F'75
!边坡在强度折减系数F=1.4时求解 d}f| HOFq
finish =goZI6 7
/solu 0/.#V*KM
allsel c$f|a$$b
asel,s,area,,2,3,1 !选择面积A2、A3 c%.f|/.k
MPCHG,4,all !把所选择单元材料号改为4 .?SClTqg
allsel F" a^ `E&
solve !进行求解 =q|fe%#
save,F1.4,db !把F=1.4时求解结果保存 lwH&4K
!边坡在强度折减系数F=1.6时求解 UT [9ER S
finish A!v-[AI[
/solu (4V1%0
allsel m[^;HwJ
asel,s,area,,2,3,1 !选择面积A2、A3 da@y*TO#i
MPCHG,5,all !把所选择单元材料号改为5 A>L(#lz#ek
allsel *D1fSu!
solve #SY 8Zv
save,F1.6,db !把F=1.6时求解结果保存 )Jx!VJ^Y
!边坡在强度折减系数F=1.8时求解 IidZ -Il
finish V,qZF=} S
/solu :C42yQAP
allsel )!dELS \ix
asel,s,area,,2,3,1 !选择面积A2、A3 C# r_qn
MPCHG,6,all !把所选择单元材料号改为6 hoC}@8_
allsel + N[dYm
solve gb:Cc,F, %
save,F1.8,db !把F=1.8时求解结果保存 81a&99k #
!边坡在强度折减系数F=2.0时求解 U]sU
b3
finish }wR)p
/solu BO\l>\)Ir
allsel |PN-,f{ -
asel,s,area,,2,3,1 !选择面积A2、A3 hBsjO3n
MPCHG,7,all !把所选择单元材料号改为7 c : *wev
allsel rOY^w9!
solve BN4dr9T
save,F2.0,db !把F=2.0时求解结果保存 Kw'Dzz%kN
!边坡在强度折减系数F=2.2时求解 ?HIc=
finish *DkA$Eu3u
/solu }jU{RR%6B
allsel )I+1 b
!U
asel,s,area,,2,3,1 !选择面积A2、A3 P7 \(D`
MPCHG,8,all !把所选择单元材料号改为8 V`69%35*@
allsel ?l,i(I
solve *Oz5I
save,F2.2,db !把F=2.2时求解结果保存 wri[#D {
!边坡在强度折减系数F=2.4时求解 }sy3M rb
finish !W$3p'8Tu
/solu 62EJ# q[
allsel 1 ,PFz
asel,s,area,,2,3,1 !选择面积A2、A3 -lL*WA `
MPCHG,9,all !把所选择单元材料号改为9 9+QLc b
allsel 3 :Co K#
solve q1dYiG.-Z
save,F2.4,db !把F=2.4时求解结果保存 !xo@i XL
!边坡在强度折减系数F=2.6时求解 q\#3 G
finish qw87B!D
/solu B_c n[?M
allsel 4BEVG&Ks
asel,s,area,,2,3,1 !选择面积A2、A3 4v2( YJ%u
MPCHG,10,all !把所选择单元材料号改为10 k5-mK{RZ
allsel .M>u:,v
solve |[34<tIN
save,F2.6,db !把F=2.6时求解结果保存 6lT'%ho}B
!边坡在强度折减系数F=2.8时求解 GWE`'V
finish ]VJcV.7`
/solu n] Dq
allsel `*oLEXYN
asel,s,area,,2,3,1 !选择面积A2、A3 \f66ipZK*
MPCHG,11,all !把所选择单元材料号改为11 \yG_wZs
allsel !.t'3~dUf$
solve 1Y9Ye?~jd
save,F2.8,db !把F=2.8时求解结果保存 @oRYQ |.R
!边坡在强度折减系数F=3.0时求解 HTqik w5X
finish [ JAd1%$3
/solu 9: ~,T H
allsel Zl, K#
asel,s,area,,2,3,1 !选择面积A2、A3 _ux 6SIyp`
MPCHG,12,all !把所选择单元材料号改为12 #gN{8Yk>
allsel X_!mZ\H7
solve s?Q `#qD
save,F3.0,db !把F=3.0时求解结果保存 &\p :VF.
/post1 ！进入后处理 2:DpnLU5
!边坡在强度折减系数F=1时结果分析 B 8s|VI
Resume,'F1','db' !读入边坡在强度折减系数F=1时 3al5Vu2:
set,1,last !读入后一个子步 5 #kvb$97
pldisp,1 !绘制边坡模型变形图 M.OWw#?p:_
plnsol,u,x !绘制边坡模型水平方向位移云图 RTl7vzG
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 &tRnI $D
!边坡在强度折减系数F=1.2时结果分析 |9\Lv $VJ
Resume,'F1.2','db' !读入边坡在强度折减系数F=1.2时 >a4Bfnf"eI
set,1,last !读入后一个子步 },Z -w_H
pldisp,1 !绘制边坡模型变形图 V GvOwd)E
plnsol,u,x !绘制边坡模型水平方向位移云图 :m$%D]WY
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 A`N;vq,
!边坡在强度折减系数F=1.4时结果分析 S'@"a%EV
Resume,'F1.4','db' !读入边坡在强度折减系数F=1.4时 q6G( [h7
set,1,last !读入后一个子步 *<0g/AL
pldisp,1 !绘制边坡模型变形图 +$
-#V
plnsol,u,x !绘制边坡模型水平方向位移云图 ' xi..
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 O3BU.X1'%
!边坡在强度折减系数F=1.6时结果分析 fGe{7p6XV*
Resume,'F1.6','db' !读入边坡在强度折减系数F=1.6时 F 1zc4l6
set,1,last !读入后一个子步 v <c Hx/
pldisp,1 !绘制边坡模型变形图 *(k=!`4(
plnsol,u,x !绘制边坡模型水平方向位移云图 U73`HDJ
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 73<iK]*c
!边坡在强度折减系数F=1.8时结果分析 ZR QPOy
Resume,'F1.8','db' !读入边坡在强度折减系数F=1.8时 {(MG:
B
set,1,last !读入后一个子步 SQ1&n;M}f
pldisp,1 !绘制边坡模型变形图 WI~%n
plnsol,u,x !绘制边坡模型水平方向位移云图 /gd o~
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 & {/ u>,
!边坡在强度折减系数F=2.0时结果分析 h =U 4
Resume,'F2.0','db' !读入边坡在强度折减系数F=2.0时 RCX4;,DHx
set,1,last !读入后一个子步 ~Igo
8ykl
pldisp,1 !绘制边坡模型变形图 -"nYCF
plnsol,u,x !绘制边坡模型水平方向位移云图 mn4;$1~e>H
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 %4-pw|':
!边坡在强度折减系数F=2.2时结果分析 `fV$'u
Resume,'F2.2','db' !读入边坡在强度折减系数F=2.2时 y ( z U:.
set,1,last !读入后一个子步 dk]ro~ [
pldisp,1 !绘制边坡模型变形图 :|A db\b
plnsol,u,x !绘制边坡模型水平方向位移云图 0ND7F
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 S8cFD):q
!边坡在强度折减系数F=2.4时结果分析 H5%I?ZXw4
Resume,'F2.4','db' !读入边坡在强度折减系数F=2.4时 u J y@
set,1,last !读入后一个子步 vSL{W T]m
pldisp,1 !绘制边坡模型变形图 e^an` </{
plnsol,u,x !绘制边坡模型水平方向位移云图 r 1a{Y8?
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 Arv8P
P^'
!边坡在强度折减系数F=2.6时结果分析 YO oP] 0'L
Resume,'F2.6','db' !读入边坡在强度折减系数F=2.6时 .KLuGb 3JJ
set,1,last !读入后一个子步 ?d0Dfqh_
pldisp,1 !绘制边坡模型变形图 1lQ1 0J
plnsol,u,x !绘制边坡模型水平方向位移云图 #P@r[VZ{6
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 ;7A, 'y4f
!边坡在强度折减系数F=2.8时结果分析 t8+93,*B
Resume,'F2.8','db' !读入边坡在强度折减系数F=2.8时 ujh4cp
set,1,last !读入后一个子步 CKur$$B
pldisp,1 !绘制边坡模型变形图 bDNd
m-
plnsol,u,x !绘制边坡模型水平方向位移云图 #e=^[E-yE
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 J|q_&MX/
!边坡在强度折减系数F=3.0时结果分析 O:x%!-w
Resume,'F3.0','db' !读入边坡在强度折减系数F=3.0时 e_;6U Z+
set,1,last !读入后一个子步 H}~^,B2;
pldisp,1 !绘制边坡模型变形图 CUG<v3\
plnsol,u,x !绘制边坡模型水平方向位移云图 8UcT? Zp
plnsol,eppl,eqv !绘制边坡模型塑性应变云图
