(1)创建物理环境 dw Aju:-H
/COM, Structural ,v)@&1Wh:
/TITLE, Slope stability Analysis P iAA,
/FILNAM,support,1 Y. Gr(]tk
(2)建立模型 Ij{{Z;o3
!进入前处理器 =I7#Vtd^K<
/PREP7 WERK JA
!定义单元类型 M;3uG/E\
ET,1,PLANE82 !定义边坡围岩单元 -?'u"*#1,
KEYOPT,1,3,2 !求解类型为平面应变 y4M<L. RO
!定义材料属性 m= j7 vb
!边坡围岩1材料属性 H> _%ZXL
MP,EX,1,3E10 CS6,mX
MP,PRXY,1,0.25 YSv\T '3
MP,DENS,1,2500 !定义塑性模型 =b !f
TB,DP,1 !定义数据表1 r 97 V X>
TBDATA,1,9E5,42.8 !给数据表输入数据 \Sg&Qv`
!边坡围岩2材料属性 !采用弹性模型 O]lWaiR`
MP,EX,2,3.2E10 ！定义弹性模量 ' +'
MP,PRXY,2,0.24 !定义泊松比 }$iH 3#E8
MP,DENS,2,2700 !定义密度 u49/LtB\
!取强度折减系数F=1.2时边坡围岩材料属性 *qKwu?]?>
MP,EX,3,3E10 8 +uOYNXsA
MP,PRXY,3,0.25 SV8rZWJ
MP,DENS,3,2500 *^" 4 )
TB,DP,3 !定义数据表3 Mx$VAV^\
TBDATA,1,7.5E5,37.7 !给数据表输入数据 f n;7Nf7{
!取强度折减系数F=1.4时边坡围岩材料属性 9\Y j`,i5
MP,EX,4,3E10 *"1]NAz+
MP,PRXY,4,0.25 !定义泊松比 xPsuDi8u
MP,DENS,4,2500 c%i/ '<Afr
TB,DP,4 !定义数据表4 n%Nf \z
TBDATA,1,6.4E5,33.5 !给数据表输入数据 Eiz\Nb
!取强度折减系数F=1.6时边坡围岩材料属性 a.c2ScXG
MP,EX,5,3E10 LFg<j1Gk`
MP,PRXY,5,0.25 52-^HV
MP,DENS,5,2500 ！定义密度 Pme`UcE3H
TB,DP,5 !定义数据表5 W%~ S~wx
TBDATA,1,5.6E5,30 !给数据表输入数据 yrkd#m
!取强度折减系数F=1.8时边坡围岩材料属性 VA2%2g2n{
MP,EX,6,3E10 ！定义弹性模量 +2C:]
MP,PRXY,6,0.25 O!/J2SfuDH
MP,DENS,6,2500 e2/&X;2
TB,DP,6 !定义数据表6 bO^%#<7
TBDATA,1,5.0E5,27.2 >JC.qjA
!取强度折减系数F=2.0时边坡围岩材料属性 =_L"x~0I-
MP,EX,7,3E10 3- LO
MP,PRXY,7,0.25 !定义泊松比 1Qf5H!5vx
MP,DENS,7,2500 ~u}[VP
TB,DP,7 !定义数据表7 w 7=Y_
TBDATA,1,4.5E5,24.8 wm@1jLjrQ
!取强度折减系数F=2.2时边坡围岩材料属性 3 7 M7bB0
MP,EX,8,3E10 (lEWnf=2h
MP,PRXY,8,0.25 QGLfZvTT
MP,DENS,8,2500 7{<t]wQq
TB,DP,8 !定义数据表8 w*Kw#m 'U
TBDATA,1,4.09E5,22.8 "&L<u0KHG
!取强度折减系数F=2.4时边坡围岩材料属性 cWh Aj>?_Q
MP,EX,9,3E10 ; :"~utL7
MP,PRXY,9,0.25 $K;4=zN>t:
MP,DENS,9,2500 ,:;nq> ;
TB,DP,9 !定义数据表9 P1z6 sG G
TBDATA,1,3.6E5,21.1 !给数据表输入数据 u4+)lvt
!取强度折减系数F=2.6时边坡围岩材料属性 !|Vjv}UO
MP,EX,10,3E10 ！定义弹性模量 m}F1sRkdQ
MP,PRXY,10,0.25 !定义泊松比 u%h]k ,(E
MP,DENS,10,2500 ！定义密度 @c7 On)sy
TB,DP,10 !定义数据表10 Ep?a1&b
TBDATA,1,3.46E5,19.6 !给数据表10输入数据 E})PNf;
!取强度折减系数F=2.8时边坡围岩材料属性 ,' 82;oP4
MP,EX,11,3E10 ！定义弹性模量 C{A eud #5
MP,PRXY,11,0.25 !定义泊松比 rmWG9&coW
MP,DENS,11,2500 y>Nlj%XH
TB,DP,11 !定义数据表11 B8[H><)o\y
TBDATA,1,3.2E5,18.2 !给数据表11输入数据 %2/Wy D$U
!取强度折减系数F=3.0时边坡围岩材料属性 jC; XY !d6
MP,EX,12,3E10 ！定义弹性模量 mL3'/3-7:V
MP,PRXY,12,0.25 (Rs <'1+>
MP,DENS,12,2500 ！定义密度 vxbO>c
TB,DP,12 !定义数据表12 \<;/)!Nmw
TBDATA,1,3.0E5,17.1 !给数据表12输入数据 V-J\!CHX
!建立几何模型 ab3" ?.3m
!创建边坡线模型线 B.{0,b W?
K,1,,,, ！ 创建关键点1 ScM2_k`D
K,2,-800,,, C "hc.A&4
K,3,-800,-800,, F"a,[i,[W
K,4,-800,-1200,, gKS^-X{x
K,5,1200,-1200,, hQHV ]xW
K,6,1200,-800,, W&Fa8
K,7,1200,,, h2uO+qEsu
K,8,1200,378,, <8j n_6
K,9,430,378,, B`T9dL[E4
l,1,2 !连接关键点1、2生成直线L1 3H4p$\; C
l,2,3 Q"QrbU
l,3,4 nQ /(*d
l,4,5 5#WZXhlc}
l,5,6 8!:4m"Y
l,6,7 c;(}Ih(#
l,7,8 nLo:\I(
l,8,9 ;k!Ej-(
l,9,1 mN ~;MR;
l,1,7 rQ~%SUM7
l,3,6 C5;"mo-
!创建边坡面模型 63F0 Za}h
al,3,4,5,11 !由线L3、L4、L5、L11生成一个面A1 7xWX:2l*?
al,1,2,11,6,10 !由线L1、L2、L1、L6、L10生成一个面A2 SM0=
al,7,8,9,10 !由线L7、L8、L9、L10生成一个面A3 #4~Ivj
/PNUM, area, on !打开面号开关 PT>,:zY
!划分网格生成有限元模型 bumS>:
! 划分边坡围岩2网格 #pOW2 Uj8\
!设置网格份数 i-tX5Md|
lsel,s,,,3,5,2 !选择线L3和L5 Sy8o/-
lesize,all,,,5 !把所选择线分为5段 xa!@$w=U&
lsel,s,,,4,11,7 !选择线L4和L11 MZ% P(5
lesize,all,,,26 !把所选择线分为26段 J,P7k$t2vv
mat,2 !给边坡围岩2赋予2号材料特性 qK(? \ t$
type,1 !采用单元类型1 (K0FWTmm
mshkey,0 !设定自由网格划分 EY> %#0
mshape,0 !设定四边形网格划分 .F%RW8=Q
amesh,1 ！划分面积A1 E/5w
H/
! 划分边坡围岩1网格 : @YZ6?hf
!设置网格份数 T[ mTA>d
lsel,s,,,2,6,4 !选择线L2和L6 i,b>&V/Y$
lsel,a,,,1 !附加选择线L1 .@ 1\26<
lesize,all,,,10 !把所选择线分为10段 ]lm9D@HMC
lsel,s,,,8,10,2 !选择线L8和L10 ) c+ ZQq
lesize,all,,,16 !把所选择线分为16段 z2 nDD6N
lsel,s,,,7,9,2 ]+fL6"OD/2
lesize,all,,,12 |@V<} 2zCZ
mat,1 !给边坡围岩1赋予1号材料特性 ){8^l0b
type,1 c$ 1ez
mshkey,0 !设定自由网格划分 |%b' L.$4
mshape,0 !设定四边形网格划分 ~EX/IIa {
amesh,2 !划分面积A2 &z% 7Nu
amesh,3 ！划分面积A3 B4U+q|OD#
allsel nA% - <
!保存网格模型 &w LI:x5
save,Slope-grid.db MPM_/dn-
(3)施加约束和荷载 s_E iA _
!给边坡模型施加约束 X0P <ifIv
!边坡两侧施加X方向约束 p?6w/ n
nsel,s,loc,x,-800 !选择X=-800线上所有节点 ~#(bX]+A
nsel,a,loc,x,1200 !选择X=1200线上所有节点 O P``g/x)
d,all,ux !对所选择节点约束X方向位移 5n?fZ?6(
allsel +F+jC9j(<
!边坡底部施加约束 6;5}%
B:#h
nsel,s,loc,y,-1200 !选择Y=--1200线上所有节点 ,4EE9
?J
d,all,uy !对所选择节点约束X、Y方向位移 [&)9|EV
d,all,ux 7j._3'M=Kc
!施加重力加速度 kdg Q -UN$
acel,,9.8 K$f~Fft
(4)求解 3#5sj >
/solu 'nW:2(J
!求解设置 ]Fs Plxk6
antype,static !设定为静力求解 CUY2eQJ{U
nsubst,100 !设定最大子步数为100 5k< HO _]
pred,on !打开时间步长预测器 %Ix^Xb0
nropt,full !设定牛顿-拉普森选项 l|5ss{llR
nlgeom,on !打开大位移效果 ,X9Y/S
l
lnsrch,on !打开线性搜索 INcJ Xlv
outres,all,all !输出所有项 CX\#
|Q8q
cnvtol,f,,0.005,2,0.5 !力收敛准则设定 U_oMR$/Z
cnvtol,u,,0.05,2,1 !位移收敛准则设定 !U=;e ?o
!边坡在强度折减系数F=1时求解 ' 71D:%p
allsel ]6 }|X#_
solve !进行求解 qItj`F)d
save,F1,db !把F=1时求解结果保存 F<G.!Y8!&
!边坡在强度折减系数F=1.2时求解 qR B&R$
finish z[CCgs&vqe
/solu Wp T.25
allsel C 5)G^
asel,s,area,,2,3,1 !选择面积A2、A3 syBYH5
MPCHG,3,all ! 把所选择单元材料号改为3 o5AyJuS-u$
allsel o2F6K*u}
solve !进行求解 MP NBA1s
save,F1.2,db !把F=1.2时求解结果保存 //| 9J(B]
!边坡在强度折减系数F=1.4时求解 bha_b j
finish >&Bg F*mm
/solu ) ~ C)4
allsel hM+nA::w
asel,s,area,,2,3,1 !选择面积A2、A3 dHd{9ftyF
MPCHG,4,all !把所选择单元材料号改为4 +F4xCz7f
allsel 5RZAs63t
solve !进行求解 d]w*fn
save,F1.4,db !把F=1.4时求解结果保存 <R_3; 5J%
!边坡在强度折减系数F=1.6时求解 ]Ns b V
finish rS 6iZp,
/solu s)&"g a
allsel !K6: W1
asel,s,area,,2,3,1 !选择面积A2、A3 +| Cvv]Tx1
MPCHG,5,all !把所选择单元材料号改为5 W99Fb+$I
allsel [+ 1([#
solve E~{-RZNK
save,F1.6,db !把F=1.6时求解结果保存 )mp0k%
!边坡在强度折减系数F=1.8时求解 rK)%n!Z
finish VYlg+MlT0
/solu S(/@.gI:f
allsel WS2TOAya)
asel,s,area,,2,3,1 !选择面积A2、A3
"&v?>
MPCHG,6,all !把所选择单元材料号改为6 5E}0 <&
allsel `M p]iD {
solve KHdj#3<AR
save,F1.8,db !把F=1.8时求解结果保存 8 rnr>Ee@
!边坡在强度折减系数F=2.0时求解 8Ck:c45v
finish .>h|e_E
/solu $6ITa }o
allsel ^VoQGP/cl
asel,s,area,,2,3,1 !选择面积A2、A3 D|_}~T>;&
MPCHG,7,all !把所选择单元材料号改为7 [=.. #y!U
allsel %qN8u Qx
solve 2u#{K9g
save,F2.0,db !把F=2.0时求解结果保存 EMJio\
!边坡在强度折减系数F=2.2时求解 +O9l@X$l=
finish # YhKAG@ |
/solu X @r5^A[9
allsel saYn\o"m
asel,s,area,,2,3,1 !选择面积A2、A3 QWfwoe&;R:
MPCHG,8,all !把所选择单元材料号改为8 ]3 Mm"7`
allsel ^c.b@BE
solve F~<$E *&h@
save,F2.2,db !把F=2.2时求解结果保存 Q_M2!qj
!边坡在强度折减系数F=2.4时求解 I" Y?vj9]
finish *>Om3[D
/solu A}[Lk#|n
allsel ?Tb'J`MO
asel,s,area,,2,3,1 !选择面积A2、A3 /kNr5s
MPCHG,9,all !把所选择单元材料号改为9 eN,m8A`/S
allsel aD0w82s]J
solve )lH?XpfTjm
save,F2.4,db !把F=2.4时求解结果保存 ka" jv"z
!边坡在强度折减系数F=2.6时求解 5.5dB2w
finish 2t#9ih"9
/solu 6n;ew l}
allsel kA\;h|Y3
asel,s,area,,2,3,1 !选择面积A2、A3 @(Q4
MPCHG,10,all !把所选择单元材料号改为10 a08B8
allsel &X +@,!
solve 7r*>? ]y+
save,F2.6,db !把F=2.6时求解结果保存 S.M< (
!边坡在强度折减系数F=2.8时求解 pOw4H67
finish kD"dZ Qx
/solu */?L_\7
allsel A5
8i}G9
asel,s,area,,2,3,1 !选择面积A2、A3 x{RTI#a.
MPCHG,11,all !把所选择单元材料号改为11 z?FZu,h}
allsel 1AA(qE
solve n |.- :Zy
save,F2.8,db !把F=2.8时求解结果保存 >8jDW "Ua
!边坡在强度折减系数F=3.0时求解 Auv/w}zrr
finish 5M*q{kX)
/solu ?Cmb3pX^\
allsel r\_aux^z
asel,s,area,,2,3,1 !选择面积A2、A3 .DCp)&m
l;
MPCHG,12,all !把所选择单元材料号改为12 .L6t3/^
allsel 9l OUE
solve (7-K4j`
save,F3.0,db !把F=3.0时求解结果保存 e`8z1r
/post1 ！进入后处理 %Sxy!gGz%%
!边坡在强度折减系数F=1时结果分析 .h Z =8y9
Resume,'F1','db' !读入边坡在强度折减系数F=1时 \h _hd%'G
set,1,last !读入后一个子步 =a7m^e7
pldisp,1 !绘制边坡模型变形图 0D 0 #*J
plnsol,u,x !绘制边坡模型水平方向位移云图 aLhTaB-va
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 <6- (a;T!7
!边坡在强度折减系数F=1.2时结果分析 vWzNsWPK"{
Resume,'F1.2','db' !读入边坡在强度折减系数F=1.2时 Z&J417buk
set,1,last !读入后一个子步 PMkwY {.u
pldisp,1 !绘制边坡模型变形图 yTbBYx9Bi
plnsol,u,x !绘制边坡模型水平方向位移云图 LOe!qt\&
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 RwT.B+Onuy
!边坡在强度折减系数F=1.4时结果分析 4Mg0 9
Resume,'F1.4','db' !读入边坡在强度折减系数F=1.4时 aLq=%fsV)
set,1,last !读入后一个子步 ;2%3~L 8?V
pldisp,1 !绘制边坡模型变形图 H+_oK
]/
plnsol,u,x !绘制边坡模型水平方向位移云图 xI_WkoI
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 M,y='* \M
!边坡在强度折减系数F=1.6时结果分析 WV?i YX!
Resume,'F1.6','db' !读入边坡在强度折减系数F=1.6时 ]FQ4v.7
set,1,last !读入后一个子步 }{E//o:Ta
pldisp,1 !绘制边坡模型变形图 bCA3w%,kM
plnsol,u,x !绘制边坡模型水平方向位移云图 /sJk[5!z
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 Fv e,& ~
!边坡在强度折减系数F=1.8时结果分析 Cg )#B+
Resume,'F1.8','db' !读入边坡在强度折减系数F=1.8时 QDxL y aL
set,1,last !读入后一个子步 lL%7 lO
pldisp,1 !绘制边坡模型变形图 t_P1a0Zu
plnsol,u,x !绘制边坡模型水平方向位移云图 G{ F>=z"(l
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 kZF\V7k
!边坡在强度折减系数F=2.0时结果分析 !A!zG)Ue<
Resume,'F2.0','db' !读入边坡在强度折减系数F=2.0时 2c9@n9Vx3a
set,1,last !读入后一个子步 OM2|c}]ZQ
pldisp,1 !绘制边坡模型变形图 {zmo7~=
plnsol,u,x !绘制边坡模型水平方向位移云图 uyAhN
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 0#<_:E
!边坡在强度折减系数F=2.2时结果分析 c S{ l2}E
Resume,'F2.2','db' !读入边坡在强度折减系数F=2.2时 EL~s90C
set,1,last !读入后一个子步 h
:NHReMT
pldisp,1 !绘制边坡模型变形图 ;
Sh |6
plnsol,u,x !绘制边坡模型水平方向位移云图 A+ Z3b:}~
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 nE::9Yh8z
!边坡在强度折减系数F=2.4时结果分析 XTHrf'BU
Resume,'F2.4','db' !读入边坡在强度折减系数F=2.4时 mX,#|qLf
set,1,last !读入后一个子步 zM[WbB+"m
pldisp,1 !绘制边坡模型变形图 eYR/kZ %<
plnsol,u,x !绘制边坡模型水平方向位移云图 $+*Zs Io
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 C:gE
!边坡在强度折减系数F=2.6时结果分析 $ #"}g#u
Resume,'F2.6','db' !读入边坡在强度折减系数F=2.6时 nLT]'B]$ +
set,1,last !读入后一个子步 0fNWI
pldisp,1 !绘制边坡模型变形图 LhV4 ^\+
plnsol,u,x !绘制边坡模型水平方向位移云图 KGK8;Q,O
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 ZFrK'BvbR
!边坡在强度折减系数F=2.8时结果分析 #v~zf@<KLB
Resume,'F2.8','db' !读入边坡在强度折减系数F=2.8时 2Uu ,Vv
set,1,last !读入后一个子步 |!IJ/ivEgw
pldisp,1 !绘制边坡模型变形图 `qs'={YtU
plnsol,u,x !绘制边坡模型水平方向位移云图 L{
.r8wSrI
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 QZQ@C# PR;
!边坡在强度折减系数F=3.0时结果分析 9YB~1 M
Resume,'F3.0','db' !读入边坡在强度折减系数F=3.0时 ;|9VPv/
set,1,last !读入后一个子步 ^JY R^X>_
pldisp,1 !绘制边坡模型变形图 AGrGZ7p]
plnsol,u,x !绘制边坡模型水平方向位移云图 t}NxD`8
plnsol,eppl,eqv !绘制边坡模型塑性应变云图
