(1)创建物理环境 |Z "h q
/COM, Structural )<Yy.Z_:DC
/TITLE, Slope stability Analysis '7=*n_l
/FILNAM,support,1 _ wqFKj
(2)建立模型 RhDa`kV%t
!进入前处理器 ~MQN&
/PREP7 1yK f=LZ^
!定义单元类型 ?Ts
Z_
ET,1,PLANE82 !定义边坡围岩单元 x'
KEYOPT,1,3,2 !求解类型为平面应变 WK<pZ *x
!定义材料属性 c76^x
!边坡围岩1材料属性 @yek6E&9
MP,EX,1,3E10 uZ'5&k96T
MP,PRXY,1,0.25 "/\:Fdc^
MP,DENS,1,2500 !定义塑性模型 Nd6N:1 -
TB,DP,1 !定义数据表1 g6*}& .&
TBDATA,1,9E5,42.8 !给数据表输入数据 ;N,7#l|wi
!边坡围岩2材料属性 !采用弹性模型 VLOyUt~O#
MP,EX,2,3.2E10 ！定义弹性模量 mV'd9(s?
MP,PRXY,2,0.24 !定义泊松比 f| apk ,o_
MP,DENS,2,2700 !定义密度 SE/@ li
!取强度折减系数F=1.2时边坡围岩材料属性 Q2#)Jx\6!
MP,EX,3,3E10 +~[19'GH
MP,PRXY,3,0.25 $hN!DHz
MP,DENS,3,2500 <4>6k7 W
TB,DP,3 !定义数据表3 T &0tW"r?
TBDATA,1,7.5E5,37.7 !给数据表输入数据 JUXK}0d%eN
!取强度折减系数F=1.4时边坡围岩材料属性 eq/s8]uM
MP,EX,4,3E10 o= 8yp2vG
MP,PRXY,4,0.25 !定义泊松比 5N[ Y2
MP,DENS,4,2500 ',CcL N
TB,DP,4 !定义数据表4 M.l;!U!}
TBDATA,1,6.4E5,33.5 !给数据表输入数据 h
knobk
!取强度折减系数F=1.6时边坡围岩材料属性 ?dTz?C.w
MP,EX,5,3E10 FEP\5d >
MP,PRXY,5,0.25 .}0Cg2W
MP,DENS,5,2500 ！定义密度 K+ GjJ8
TB,DP,5 !定义数据表5 x t-;7
TBDATA,1,5.6E5,30 !给数据表输入数据 Dljq
!取强度折减系数F=1.8时边坡围岩材料属性 B$lbp03z
MP,EX,6,3E10 ！定义弹性模量 #p [= iP
MP,PRXY,6,0.25 m.4y=69 &
MP,DENS,6,2500 >MhkNy
TB,DP,6 !定义数据表6 Q.8Jgel1
TBDATA,1,5.0E5,27.2 f.8L<<5 c
!取强度折减系数F=2.0时边坡围岩材料属性 * dNMnZ@Y
MP,EX,7,3E10 @r
.K>+1
MP,PRXY,7,0.25 !定义泊松比 ,Y&kW'2
MP,DENS,7,2500 kPxrI=
TB,DP,7 !定义数据表7 p<J/J.E
TBDATA,1,4.5E5,24.8 {fS/ZG"5<t
!取强度折减系数F=2.2时边坡围岩材料属性 "fmJ;W;#1
MP,EX,8,3E10 Rn TP U`
MP,PRXY,8,0.25 ?c43cYb
MP,DENS,8,2500 O=+C Kx@
TB,DP,8 !定义数据表8 |Rab'9U^
TBDATA,1,4.09E5,22.8 *]H ./a:1
!取强度折减系数F=2.4时边坡围岩材料属性 t
Y^:C[
MP,EX,9,3E10 R.RCa$
MP,PRXY,9,0.25 ksK
lw_%o
MP,DENS,9,2500 &0o&!P8CB
TB,DP,9 !定义数据表9 \K)q$E <!
TBDATA,1,3.6E5,21.1 !给数据表输入数据 -BjB>Vt
!取强度折减系数F=2.6时边坡围岩材料属性 v/m6 (z
MP,EX,10,3E10 ！定义弹性模量 ! AMPA*
MP,PRXY,10,0.25 !定义泊松比 i\ PN
MP,DENS,10,2500 ！定义密度 $MR{3-
TB,DP,10 !定义数据表10 j5 RM S V
TBDATA,1,3.46E5,19.6 !给数据表10输入数据 $4eogI7N>w
!取强度折减系数F=2.8时边坡围岩材料属性 g| T' oK
MP,EX,11,3E10 ！定义弹性模量 f< '~K
MP,PRXY,11,0.25 !定义泊松比 u{_T,k<!
MP,DENS,11,2500 ,qr)}s-
TB,DP,11 !定义数据表11 ; aMMI p
TBDATA,1,3.2E5,18.2 !给数据表11输入数据 iE&`F hf?
!取强度折减系数F=3.0时边坡围岩材料属性 WFh!re%Z
MP,EX,12,3E10 ！定义弹性模量 ] #J ]f
MP,PRXY,12,0.25 | e pe;/
MP,DENS,12,2500 ！定义密度 ao,LP,_
TB,DP,12 !定义数据表12 2Hq!YsJ4]
TBDATA,1,3.0E5,17.1 !给数据表12输入数据 *X)OdU
!建立几何模型 c(eu[vj:
!创建边坡线模型线 B)c.`cfr*\
K,1,,,, ！ 创建关键点1 K1 C#
K,2,-800,,, #6YNgJNk
K,3,-800,-800,, CBF>157B
K,4,-800,-1200,, 2[:`w),.
K,5,1200,-1200,, %b_zUFHPp
K,6,1200,-800,, h<QXr '4+
K,7,1200,,, z2 4-h C
K,8,1200,378,, }ShZ4 xMz
K,9,430,378,, V&f3>#n\
l,1,2 !连接关键点1、2生成直线L1 g:*yjj
l,2,3 sB"]R%`_
l,3,4 EoY570PN
l,4,5 Y${ $7+@
l,5,6 T&{EqsI=B
l,6,7 1o$<pZZ
l,7,8 6Es-{u(,
l,8,9 fNlU c
l,9,1 lc'Jn$O@
l,1,7 jbg@ CA*=C
l,3,6 }L E/{]A
!创建边坡面模型 6DExsB~@
al,3,4,5,11 !由线L3、L4、L5、L11生成一个面A1 k r/[|.bq
al,1,2,11,6,10 !由线L1、L2、L1、L6、L10生成一个面A2 eH6#'M4+\
al,7,8,9,10 !由线L7、L8、L9、L10生成一个面A3 CE+\|5u
W
/PNUM, area, on !打开面号开关 |E}N 8 \Gr
!划分网格生成有限元模型 vu*08<M~i|
! 划分边坡围岩2网格 N,;Bl&EU
!设置网格份数 SWm dU]
lsel,s,,,3,5,2 !选择线L3和L5 K3@ UoR
lesize,all,,,5 !把所选择线分为5段 `@:^(sMo
lsel,s,,,4,11,7 !选择线L4和L11 t[ DXG2&
lesize,all,,,26 !把所选择线分为26段 0_j! t
mat,2 !给边坡围岩2赋予2号材料特性 uS&bfx2
type,1 !采用单元类型1 '7xY ,IY
mshkey,0 !设定自由网格划分 v8Wo V*
mshape,0 !设定四边形网格划分 S]9xqiJ W
amesh,1 ！划分面积A1 T0Gu(c`1d
! 划分边坡围岩1网格 &-{4 JSII
!设置网格份数 *=ALns?y
lsel,s,,,2,6,4 !选择线L2和L6 <ZnAPh
lsel,a,,,1 !附加选择线L1 j0 OxR.S
lesize,all,,,10 !把所选择线分为10段 ?vk&k (FT
lsel,s,,,8,10,2 !选择线L8和L10 {X< t Uco
lesize,all,,,16 !把所选择线分为16段 OgzPX^q/=
lsel,s,,,7,9,2 ;=E3f ^'s
lesize,all,,,12 DG&
kY+
mat,1 !给边坡围岩1赋予1号材料特性 KQ 2]VN"?_
type,1 MqNp*n 2
mshkey,0 !设定自由网格划分 %f>V\z_C
mshape,0 !设定四边形网格划分 O!|:ZMjF
amesh,2 !划分面积A2 Ml)WY#7
amesh,3 ！划分面积A3 :oZ 30}
allsel >=B 8PK+<
!保存网格模型 Lu <'A4Q1
save,Slope-grid.db k!! o!r BS
(3)施加约束和荷载 gfp#G,/B
!给边坡模型施加约束 3_ D$6/i
!边坡两侧施加X方向约束 p2cKtk+
nsel,s,loc,x,-800 !选择X=-800线上所有节点 3 =S.-
nsel,a,loc,x,1200 !选择X=1200线上所有节点 i,V~5dE[I<
d,all,ux !对所选择节点约束X方向位移 f:= ?"MX7
allsel /
AW]12_
!边坡底部施加约束 %f^TZ,q$
nsel,s,loc,y,-1200 !选择Y=--1200线上所有节点 %i96@ 6O
d,all,uy !对所选择节点约束X、Y方向位移 .]jKuTC\<
d,all,ux |M+ !O93
!施加重力加速度 {0[qERj"z
acel,,9.8 K~X t`
(4)求解 *W0`+#Dcv
/solu XCk \#(VSE
!求解设置 DsP +#PX
antype,static !设定为静力求解 xo]|m\#k5E
nsubst,100 !设定最大子步数为100 >uI|S
pred,on !打开时间步长预测器 kdv>QZ
nropt,full !设定牛顿-拉普森选项 Kj}}O2
nlgeom,on !打开大位移效果 Uyv FR @
lnsrch,on !打开线性搜索 /njN*rhx&Z
outres,all,all !输出所有项 <7)@Jds\
cnvtol,f,,0.005,2,0.5 !力收敛准则设定 \75%[;.
cnvtol,u,,0.05,2,1 !位移收敛准则设定 kgA')]
!边坡在强度折减系数F=1时求解 Q#vur o
allsel ++FMkeHZ
solve !进行求解 |D<~a(0
save,F1,db !把F=1时求解结果保存 gjwp' GN
!边坡在强度折减系数F=1.2时求解 xvW+;3;
finish .m4K ]^m
/solu )! OEa]
allsel \BS^="AcpP
asel,s,area,,2,3,1 !选择面积A2、A3 6 .*=1P*?
MPCHG,3,all ! 把所选择单元材料号改为3 4w\')@`[jk
allsel ZOU$do>O
solve !进行求解 $ A ( #^&
save,F1.2,db !把F=1.2时求解结果保存 d_OHQ pfK
!边坡在强度折减系数F=1.4时求解 V%3K")
finish Y pp>7J/
/solu nGg>lRL
allsel v/(< fI^
asel,s,area,,2,3,1 !选择面积A2、A3 ;[*7UE+#7
MPCHG,4,all !把所选择单元材料号改为4 |}#Rn`*2y
allsel F 02Nn F
solve !进行求解 3ldOOQW%
save,F1.4,db !把F=1.4时求解结果保存 iHR?]]RF
!边坡在强度折减系数F=1.6时求解 PDtaL
finish WSh+5](:
/solu <Z}2A8mjY
allsel gB&'MA!
asel,s,area,,2,3,1 !选择面积A2、A3 @90)
MPCHG,5,all !把所选择单元材料号改为5 ?6a:!^eL
allsel Hm*n ,8_
solve 6 @ nE cr
save,F1.6,db !把F=1.6时求解结果保存 +nZx{d,wt
!边坡在强度折减系数F=1.8时求解 =b 6Q2s,i
finish !,I}2,1%k
/solu \.}* s] 6
allsel ,(]hykbXp
asel,s,area,,2,3,1 !选择面积A2、A3 5Rc
5/ m
MPCHG,6,all !把所选择单元材料号改为6 F*(<`V
allsel I|08 [
mO
solve m 'a3}vRV(
save,F1.8,db !把F=1.8时求解结果保存 yA6"8fr
!边坡在强度折减系数F=2.0时求解 UW40Y3W0
finish K 0b(D8!
/solu "&>$/b$
allsel fRq2sK;+
asel,s,area,,2,3,1 !选择面积A2、A3 }7UE
MPCHG,7,all !把所选择单元材料号改为7 ;Qw>&24h[
allsel tcmG>^YM
solve F_@PSA+
save,F2.0,db !把F=2.0时求解结果保存 {@({po
!边坡在强度折减系数F=2.2时求解 E 5Z,4B
finish wG -X833\(
/solu I V!&jL
allsel zg "<N
asel,s,area,,2,3,1 !选择面积A2、A3 5z=;q !3
MPCHG,8,all !把所选择单元材料号改为8 Xd{"+'29
allsel obY5taOw
solve gx #TRp}-
save,F2.2,db !把F=2.2时求解结果保存 }yX a1#3
!边坡在强度折减系数F=2.4时求解 QQ pe.oF
finish k(V#{
YP
/solu ;K`qSX;;c(
allsel S 3.Pqp_<
asel,s,area,,2,3,1 !选择面积A2、A3 TqzkF7;k4
MPCHG,9,all !把所选择单元材料号改为9 $n>.;CV
allsel MQLa+I,S4
solve 8+lM6O ~!
save,F2.4,db !把F=2.4时求解结果保存 3'IF? ](]U
!边坡在强度折减系数F=2.6时求解
rq+E"Uj ?
finish (L
q^C=
/solu T?Hs_u {
allsel @m# 7E4 +
asel,s,area,,2,3,1 !选择面积A2、A3 /} (w{6C
MPCHG,10,all !把所选择单元材料号改为10 #lF8"@)a-$
allsel <w 8*Ly:L
solve s,lrw~17
save,F2.6,db !把F=2.6时求解结果保存 ^e)KEkh
!边坡在强度折减系数F=2.8时求解 R5| c4v{B
finish R ]HHbD&;
/solu &i6WVNGy
allsel & [4 Gv61
asel,s,area,,2,3,1 !选择面积A2、A3 z0doL b^!
MPCHG,11,all !把所选择单元材料号改为11 D4[t@* m>7
allsel poQY X5
solve iN=-N =
save,F2.8,db !把F=2.8时求解结果保存 @ VVBl I
!边坡在强度折减系数F=3.0时求解 N^:)U"9*e
finish v=@Z,-
/solu \z 0"
allsel \V}?K0#bt
asel,s,area,,2,3,1 !选择面积A2、A3 ~-|K5
MPCHG,12,all !把所选择单元材料号改为12 #
bP1rQ0
allsel Bg Uf:PT
solve P T|t6V"wd
save,F3.0,db !把F=3.0时求解结果保存 h _f A
/post1 ！进入后处理 +Z /Pj_.o
!边坡在强度折减系数F=1时结果分析 c:4 i&|n
Resume,'F1','db' !读入边坡在强度折减系数F=1时 td%EbxJK]`
set,1,last !读入后一个子步 `WX @1 ]m
pldisp,1 !绘制边坡模型变形图 V"k * PLt
plnsol,u,x !绘制边坡模型水平方向位移云图 :+Y+5:U]
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 ^ WidA-
!边坡在强度折减系数F=1.2时结果分析 s [@I I ]
Resume,'F1.2','db' !读入边坡在强度折减系数F=1.2时 0~)cAK us
set,1,last !读入后一个子步 ]B\H
pldisp,1 !绘制边坡模型变形图 D1#fy=u69|
plnsol,u,x !绘制边坡模型水平方向位移云图 B`9'COw
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 w
{6k U
!边坡在强度折减系数F=1.4时结果分析 ? !cUAa>iH
Resume,'F1.4','db' !读入边坡在强度折减系数F=1.4时 vz/.* u
set,1,last !读入后一个子步 f)/Yru. ;
pldisp,1 !绘制边坡模型变形图 3 /oVl
6
plnsol,u,x !绘制边坡模型水平方向位移云图 ub7|'+5
plnsol,eppl,eqv !绘制边坡模型塑性应变云图
_6xC4@~h*
!边坡在强度折减系数F=1.6时结果分析 O~trv,?)
Resume,'F1.6','db' !读入边坡在强度折减系数F=1.6时 e=Kf<ZQt
set,1,last !读入后一个子步 -NH c~ =m
pldisp,1 !绘制边坡模型变形图 sBB>O@4
plnsol,u,x !绘制边坡模型水平方向位移云图 #r >
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 5Z1b9.;.,
!边坡在强度折减系数F=1.8时结果分析 xyBWV]Y
Resume,'F1.8','db' !读入边坡在强度折减系数F=1.8时 A6pPx1 -&
set,1,last !读入后一个子步 ?/5WM%
pldisp,1 !绘制边坡模型变形图 I 6X_DPY
plnsol,u,x !绘制边坡模型水平方向位移云图 3~%9;.I3!
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 m.Yj{u8zX
!边坡在强度折减系数F=2.0时结果分析 6f{K j)
Resume,'F2.0','db' !读入边坡在强度折减系数F=2.0时 SW#
5px`
set,1,last !读入后一个子步 A^*0{F?,)
pldisp,1 !绘制边坡模型变形图 SfEgmp-m
plnsol,u,x !绘制边坡模型水平方向位移云图 &Z#g/Hc
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 %h( J+_"L6
!边坡在强度折减系数F=2.2时结果分析 gT?:zd=;
Resume,'F2.2','db' !读入边坡在强度折减系数F=2.2时 #]cO]
I
set,1,last !读入后一个子步 X\V1c$13CK
pldisp,1 !绘制边坡模型变形图 | *K AqTO0
plnsol,u,x !绘制边坡模型水平方向位移云图 e@w-4G(;
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 IP9mv`[
!边坡在强度折减系数F=2.4时结果分析 %?@N-$j
Resume,'F2.4','db' !读入边坡在强度折减系数F=2.4时 bw&8"k>D?
set,1,last !读入后一个子步 g >u { H:
pldisp,1 !绘制边坡模型变形图 Jvgx+{Xu
plnsol,u,x !绘制边坡模型水平方向位移云图 =c[9:&5Q
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 M9M~[[
!边坡在强度折减系数F=2.6时结果分析 Hq ]f$Q6:
Resume,'F2.6','db' !读入边坡在强度折减系数F=2.6时 R:fERj<s
set,1,last !读入后一个子步 .\" .}4qQ
pldisp,1 !绘制边坡模型变形图 SI^!e1@M[
plnsol,u,x !绘制边坡模型水平方向位移云图 M Zl6 J
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 #px74EeI\
!边坡在强度折减系数F=2.8时结果分析 A )tGB&
Resume,'F2.8','db' !读入边坡在强度折减系数F=2.8时 y)C nH4{
set,1,last !读入后一个子步 1 cvoI
pldisp,1 !绘制边坡模型变形图 -Uh3A\#(
plnsol,u,x !绘制边坡模型水平方向位移云图 #2^eGhwnI
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 LJ\uRfs
!边坡在强度折减系数F=3.0时结果分析 [x{'NwP?
Resume,'F3.0','db' !读入边坡在强度折减系数F=3.0时 JPR o<jt=
set,1,last !读入后一个子步 ko2 ?q
pldisp,1 !绘制边坡模型变形图 Z vM~]8m
plnsol,u,x !绘制边坡模型水平方向位移云图 luY#l!mx3
plnsol,eppl,eqv !绘制边坡模型塑性应变云图
