(1)创建物理环境 j'7FTVmJ
/COM, Structural PY^Yx$t9
/TITLE, Slope stability Analysis PC9:n ee
/FILNAM,support,1 ';1
c
(2)建立模型 lu >>~vy6
!进入前处理器 FSM~Rl
/PREP7 cZwQ{9>
!定义单元类型 HsO=%bb
ET,1,PLANE82 !定义边坡围岩单元 v7ae^iU
KEYOPT,1,3,2 !求解类型为平面应变 tHGK<r b
!定义材料属性 0nv3JX^l]
!边坡围岩1材料属性 "evLI?
MP,EX,1,3E10 I*#~@:4*
MP,PRXY,1,0.25 |BXq8Erh
MP,DENS,1,2500 !定义塑性模型 A.y$.(
TB,DP,1 !定义数据表1 `Q{ki y
TBDATA,1,9E5,42.8 !给数据表输入数据 ^ =t yf&"
!边坡围岩2材料属性 !采用弹性模型 ;w1h)
MP,EX,2,3.2E10 ！定义弹性模量 @}FRiPo6
MP,PRXY,2,0.24 !定义泊松比 UL oTPx@N
MP,DENS,2,2700 !定义密度 cJqPcCq(wn
!取强度折减系数F=1.2时边坡围岩材料属性
W6a2I
MP,EX,3,3E10 D`LwW` 9
MP,PRXY,3,0.25 c`iSe$ eS
MP,DENS,3,2500 J)R2O4OEd
TB,DP,3 !定义数据表3 *kDXx&7B$
TBDATA,1,7.5E5,37.7 !给数据表输入数据 D.HAp+lx
!取强度折减系数F=1.4时边坡围岩材料属性 Jus )cO#I
MP,EX,4,3E10 _ p?q/-[4
MP,PRXY,4,0.25 !定义泊松比 sG\=_-"v(
MP,DENS,4,2500 !Zs;m`j&9
TB,DP,4 !定义数据表4 s# d>yx_b
TBDATA,1,6.4E5,33.5 !给数据表输入数据 Pg`^EJ+
!取强度折减系数F=1.6时边坡围岩材料属性 H(0d(c1s
MP,EX,5,3E10 =B; )h
MP,PRXY,5,0.25 AD8~
MP,DENS,5,2500 ！定义密度 A\=:h AQ
TB,DP,5 !定义数据表5 &AoWT:Ea
TBDATA,1,5.6E5,30 !给数据表输入数据 Y(&phv&
!取强度折减系数F=1.8时边坡围岩材料属性 \TKv3N
MP,EX,6,3E10 ！定义弹性模量 UboOIx5:
MP,PRXY,6,0.25 fBQZ=zh
MP,DENS,6,2500 >s1HQSe66
TB,DP,6 !定义数据表6 au GN~"n^
TBDATA,1,5.0E5,27.2 /
{A]('t
!取强度折减系数F=2.0时边坡围岩材料属性 Uqkh@-6-
MP,EX,7,3E10 -5E<B mM
MP,PRXY,7,0.25 !定义泊松比 YN\
Q wV
MP,DENS,7,2500 p,xM7 V"O)
TB,DP,7 !定义数据表7 0a(*/u
TBDATA,1,4.5E5,24.8 | o+ vpy
!取强度折减系数F=2.2时边坡围岩材料属性 6z/8n f +u
MP,EX,8,3E10 ag] nVE/
MP,PRXY,8,0.25 EmLPq!C
MP,DENS,8,2500 /1O6;'8He
TB,DP,8 !定义数据表8 M T0}MMr
TBDATA,1,4.09E5,22.8 Vv zd>yII
!取强度折减系数F=2.4时边坡围岩材料属性 bi,%QZZ
MP,EX,9,3E10 6jKM,%l
MP,PRXY,9,0.25 YpuA ,r;"
MP,DENS,9,2500 PZ"xW0"-
TB,DP,9 !定义数据表9 Ue8_Q8q 5
TBDATA,1,3.6E5,21.1 !给数据表输入数据 u]zb<)'_
!取强度折减系数F=2.6时边坡围岩材料属性 wl! 'Bck=
MP,EX,10,3E10 ！定义弹性模量 F`Pu$>8C
MP,PRXY,10,0.25 !定义泊松比 xE + Go
MP,DENS,10,2500 ！定义密度 of(Nq@
TB,DP,10 !定义数据表10 hz Ppw.
TBDATA,1,3.46E5,19.6 !给数据表10输入数据 zO! `sPP
!取强度折减系数F=2.8时边坡围岩材料属性 p4t(xm2T
MP,EX,11,3E10 ！定义弹性模量 PY` V]|J
MP,PRXY,11,0.25 !定义泊松比 RzyEA3L'
MP,DENS,11,2500 Cu!4ha.e`
TB,DP,11 !定义数据表11 V,2O `D%
TBDATA,1,3.2E5,18.2 !给数据表11输入数据 Gk!v-h9cq
!取强度折减系数F=3.0时边坡围岩材料属性 L*FnFRhU
MP,EX,12,3E10 ！定义弹性模量 YXBS!89m
MP,PRXY,12,0.25 NeNKOW#X
MP,DENS,12,2500 ！定义密度 2*Gl|@~N
TB,DP,12 !定义数据表12 8fdOV&&D~i
TBDATA,1,3.0E5,17.1 !给数据表12输入数据 5|AZ/!rb
!建立几何模型 |>[w $
!创建边坡线模型线 LxMOs Nv
K,1,,,, ！ 创建关键点1 IMV oNKW-
K,2,-800,,, :s8,i$Ex
K,3,-800,-800,, -"(e*&TJ#
K,4,-800,-1200,, rPQ$e!m1Ee
K,5,1200,-1200,, Th"7p:SE?
K,6,1200,-800,, CNP?i(Rk
K,7,1200,,, ph69u #Og
K,8,1200,378,, !lHsJ)t
K,9,430,378,, OD*DHC2rN]
l,1,2 !连接关键点1、2生成直线L1 y 5=J6a2.
l,2,3 Z3[,X w
l,3,4 m
81\cg
l,4,5 Q,jlKgB 5:
l,5,6 CIV6 Qe"<
l,6,7 }P^ n /
l,7,8 ,kLe K{
l,8,9 G-rN?R.
l,9,1 &9^c-;Vs
l,1,7 1f~_# EIC
l,3,6 +!w?g/dV
!创建边坡面模型 O 89BN6p
al,3,4,5,11 !由线L3、L4、L5、L11生成一个面A1 G|H\(3hHLZ
al,1,2,11,6,10 !由线L1、L2、L1、L6、L10生成一个面A2 p|W:;(
al,7,8,9,10 !由线L7、L8、L9、L10生成一个面A3 1 t g
/PNUM, area, on !打开面号开关 4CNK ]2
!划分网格生成有限元模型 }LE.kd&
! 划分边坡围岩2网格 smk0 * m4
!设置网格份数 f!1K GP
lsel,s,,,3,5,2 !选择线L3和L5 =#1/<q)L
lesize,all,,,5 !把所选择线分为5段 KX*e2 /0
lsel,s,,,4,11,7 !选择线L4和L11 D4`7,JC}<
lesize,all,,,26 !把所选择线分为26段 bv|v9_i
mat,2 !给边坡围岩2赋予2号材料特性 xXLKL6F(\
type,1 !采用单元类型1 Ih"f98l V
mshkey,0 !设定自由网格划分 )Q9 J,
mshape,0 !设定四边形网格划分 CuDU~)`
amesh,1 ！划分面积A1 'zRd?Z>%
! 划分边坡围岩1网格 qf
]ax!bK
!设置网格份数 NH|v`r O
lsel,s,,,2,6,4 !选择线L2和L6 .o
fYFK
lsel,a,,,1 !附加选择线L1 ;/ p)vR
lesize,all,,,10 !把所选择线分为10段 > a^H7kp
lsel,s,,,8,10,2 !选择线L8和L10 s;$
eq);
lesize,all,,,16 !把所选择线分为16段 M~3(4,
lsel,s,,,7,9,2 LG51e7_gFi
lesize,all,,,12 o q;}q
mat,1 !给边坡围岩1赋予1号材料特性 GfV Mj 7{
type,1 UIzk-.<
mshkey,0 !设定自由网格划分 Qa"R?d fr
mshape,0 !设定四边形网格划分 "%0R R?
amesh,2 !划分面积A2 `(1 6 _a
amesh,3 ！划分面积A3 r\L:JTZ$
allsel r?H {Y3 ,
!保存网格模型 HaLEQ73
save,Slope-grid.db S0kH/ A
(3)施加约束和荷载 gSC8qip
!给边坡模型施加约束 =(v/pLLK?
!边坡两侧施加X方向约束 g[!sGa &
nsel,s,loc,x,-800 !选择X=-800线上所有节点 X/'B*y'=U
nsel,a,loc,x,1200 !选择X=1200线上所有节点 $^`hu%s,~
d,all,ux !对所选择节点约束X方向位移 yUBic~S
allsel drF"kTD"7
!边坡底部施加约束 ~>)cY{wE_
nsel,s,loc,y,-1200 !选择Y=--1200线上所有节点 ,{YC|uB
d,all,uy !对所选择节点约束X、Y方向位移 NW^}u ~-f
d,all,ux [<OMv9(l'o
!施加重力加速度 P]INYH
acel,,9.8 E-n!3RQ(w
(4)求解 n j1 cqh
/solu j^:\a\-1
!求解设置 {dCk iF
antype,static !设定为静力求解 92eS*x2@
nsubst,100 !设定最大子步数为100
%lEPFp
pred,on !打开时间步长预测器 jj.)$|&#`
nropt,full !设定牛顿-拉普森选项 ./.E=,j
nlgeom,on !打开大位移效果 ]+,Z()
lnsrch,on !打开线性搜索 CYG'W FvZZ
outres,all,all !输出所有项 -cJ(iz9!
cnvtol,f,,0.005,2,0.5 !力收敛准则设定 V@T G"YF
cnvtol,u,,0.05,2,1 !位移收敛准则设定 "PnYa)?1
!边坡在强度折减系数F=1时求解 -3ha LdRk6
allsel }| Cw]GW
solve !进行求解 QZ
h|6&yI
save,F1,db !把F=1时求解结果保存 )W uuU [(
!边坡在强度折减系数F=1.2时求解 m;8_A|$A
finish |0m h*+i
/solu \:^$ZBQr<n
allsel DR^mT$
asel,s,area,,2,3,1 !选择面积A2、A3 Vs m06Rj{
MPCHG,3,all ! 把所选择单元材料号改为3 jhK&Z7;
allsel XWk/S $-d
solve !进行求解 7dN]OUdi
save,F1.2,db !把F=1.2时求解结果保存 P0UMMn\-#
!边坡在强度折减系数F=1.4时求解 F;`es%8
finish vPpbm
/solu B*K%&w10~
allsel &V iIxJZ1$
asel,s,area,,2,3,1 !选择面积A2、A3 nhQ.U>&-M
MPCHG,4,all !把所选择单元材料号改为4 k'
pu%nWN
allsel 0jE,=<W0>
solve !进行求解 !0E$9 Xon
save,F1.4,db !把F=1.4时求解结果保存 ^ !E;+o' t
!边坡在强度折减系数F=1.6时求解 ;I]TM#qGF
finish mqwN<:
/solu R 3%T}^;f
allsel ?=f\oH$
asel,s,area,,2,3,1 !选择面积A2、A3 u=k\]W-
MPCHG,5,all !把所选择单元材料号改为5 S ^2'O7uj
allsel cXt L3T+
solve 69\0$ O
save,F1.6,db !把F=1.6时求解结果保存 c[J#Hc8;
!边坡在强度折减系数F=1.8时求解 SO8Ej) m
finish C[ <&% =
/solu TkmN.@w _C
allsel `F u|50_@V
asel,s,area,,2,3,1 !选择面积A2、A3 WY0u9M4
MPCHG,6,all !把所选择单元材料号改为6 aD 24)?db-
allsel Ow+7o@$"/
solve >8 VfijK
save,F1.8,db !把F=1.8时求解结果保存 < _c84,[V
!边坡在强度折减系数F=2.0时求解 6WI_JbT~
finish @q[-, EA9
/solu T\. 8og
allsel j&[63XSe
asel,s,area,,2,3,1 !选择面积A2、A3 I3s'44
MPCHG,7,all !把所选择单元材料号改为7 SAly~(r?/
allsel ` hI1
solve v;jrAND
save,F2.0,db !把F=2.0时求解结果保存 Ff%m.A8d,4
!边坡在强度折减系数F=2.2时求解 !>g c!8Y'o
finish < P$b$fh/
/solu -&@[]/
allsel e Hd{'J<
asel,s,area,,2,3,1 !选择面积A2、A3 v>7t J[s
MPCHG,8,all !把所选择单元材料号改为8 q>!T*BQ
allsel tgc&DT; E
solve D Eqk9Exk`
save,F2.2,db !把F=2.2时求解结果保存 $aG]V-M>
!边坡在强度折减系数F=2.4时求解 nolTvqMT
finish .# rI9op
/solu uh`~K6&*\w
allsel iTqv=
asel,s,area,,2,3,1 !选择面积A2、A3 *E+) mB"~
MPCHG,9,all !把所选择单元材料号改为9 8t0i
j
allsel ]:m*7p\uk
solve ",Ek| z
save,F2.4,db !把F=2.4时求解结果保存 ur9 -F^$
!边坡在强度折减系数F=2.6时求解 7A3 e-51 >
finish K xh)'aal
/solu f5 mk\^
allsel ,Sghi&Ky
asel,s,area,,2,3,1 !选择面积A2、A3 g) u2
MPCHG,10,all !把所选择单元材料号改为10 y)Ip\.KV\
allsel \`|,wLgH
solve ^
chlAQz(
save,F2.6,db !把F=2.6时求解结果保存 5o rA# B
!边坡在强度折减系数F=2.8时求解 Scxf5x-
finish C
Hyb{:<
/solu q d!$ nr
allsel .E[k}{k,
asel,s,area,,2,3,1 !选择面积A2、A3 '00DUUa
MPCHG,11,all !把所选择单元材料号改为11 d=N5cCqq
allsel aH, 0 + |
solve Eu1 s
save,F2.8,db !把F=2.8时求解结果保存 k&yy_r
!边坡在强度折减系数F=3.0时求解 B YNOgB1
finish Zff-Hl
/solu |VaXOdD`&
allsel 93[`1_q7\
asel,s,area,,2,3,1 !选择面积A2、A3 _aa3Qw x
MPCHG,12,all !把所选择单元材料号改为12 h9g 5W'.#
allsel \A keC 6[D
solve @2kt6
W
save,F3.0,db !把F=3.0时求解结果保存 w4RP*Da?:
/post1 ！进入后处理 D#^v=U
!边坡在强度折减系数F=1时结果分析 Lr40rLx;u
Resume,'F1','db' !读入边坡在强度折减系数F=1时 _EjS(.e/=
set,1,last !读入后一个子步 *kZJ
pldisp,1 !绘制边坡模型变形图 |p .mA-81
plnsol,u,x !绘制边坡模型水平方向位移云图 H)t8 d_^|j
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 'X@j
!边坡在强度折减系数F=1.2时结果分析 A -c3B+
Resume,'F1.2','db' !读入边坡在强度折减系数F=1.2时 z }u
set,1,last !读入后一个子步 EC?!%iO`
pldisp,1 !绘制边坡模型变形图 -%%2Pz0I
plnsol,u,x !绘制边坡模型水平方向位移云图 U%w ?muJW
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 /of,4aaK7
!边坡在强度折减系数F=1.4时结果分析 &eK8v]|"W
Resume,'F1.4','db' !读入边坡在强度折减系数F=1.4时 R"=G?d)
set,1,last !读入后一个子步 8%vk"h:u:
pldisp,1 !绘制边坡模型变形图 ,B!u*
plnsol,u,x !绘制边坡模型水平方向位移云图 GS< ,adD
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 CNf eHMT
!边坡在强度折减系数F=1.6时结果分析 |]\qI
Resume,'F1.6','db' !读入边坡在强度折减系数F=1.6时 E83{4A4
set,1,last !读入后一个子步 {IqbO>|"O_
pldisp,1 !绘制边坡模型变形图 N ?0V0B
plnsol,u,x !绘制边坡模型水平方向位移云图 (fY (-
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 sE-"TNONZ
!边坡在强度折减系数F=1.8时结果分析 jc)D*Cf
Resume,'F1.8','db' !读入边坡在强度折减系数F=1.8时 }g>&l.2X
set,1,last !读入后一个子步 |/YT.c%
pldisp,1 !绘制边坡模型变形图 } 3M\&}=8
plnsol,u,x !绘制边坡模型水平方向位移云图 ny TfTn
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 IjJ3CJ<
!边坡在强度折减系数F=2.0时结果分析 'rq
[P ",
Resume,'F2.0','db' !读入边坡在强度折减系数F=2.0时 g0_8:Gs}^
set,1,last !读入后一个子步 YI&^j2
pldisp,1 !绘制边坡模型变形图 eX;"kO
plnsol,u,x !绘制边坡模型水平方向位移云图 "d%":F(
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 "m/0>UU0
!边坡在强度折减系数F=2.2时结果分析 R &T(S
Resume,'F2.2','db' !读入边坡在强度折减系数F=2.2时 Rz*%(2Vz
set,1,last !读入后一个子步 T9N /;3
pldisp,1 !绘制边坡模型变形图 eUx|_*`
plnsol,u,x !绘制边坡模型水平方向位移云图 J"I{0> @
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 @ ;rU #
!边坡在强度折减系数F=2.4时结果分析 3mr9}P9;
Resume,'F2.4','db' !读入边坡在强度折减系数F=2.4时 L[##w?Xf.
set,1,last !读入后一个子步 Gi Zy C
pldisp,1 !绘制边坡模型变形图 k Af2g
plnsol,u,x !绘制边坡模型水平方向位移云图 /d8PDc "
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 v=9:N/ sW
!边坡在强度折减系数F=2.6时结果分析 :;)K>g,b
Resume,'F2.6','db' !读入边坡在强度折减系数F=2.6时 * J~N
set,1,last !读入后一个子步 ^EM##Ss_
pldisp,1 !绘制边坡模型变形图 :);GeZ
plnsol,u,x !绘制边坡模型水平方向位移云图 pPeS4$Y
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 sq&$
!边坡在强度折减系数F=2.8时结果分析 CbI[K|
Resume,'F2.8','db' !读入边坡在强度折减系数F=2.8时 VY Va8[}
set,1,last !读入后一个子步 o!\Vk~Vi&
pldisp,1 !绘制边坡模型变形图 V
mKMj'
plnsol,u,x !绘制边坡模型水平方向位移云图 %'k^aq FL
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 ks:Z=%o
!边坡在强度折减系数F=3.0时结果分析 NfDg=[FN[
Resume,'F3.0','db' !读入边坡在强度折减系数F=3.0时 v3 4XcA
set,1,last !读入后一个子步 vSPkm)O0)
pldisp,1 !绘制边坡模型变形图 ]NG`MZ
plnsol,u,x !绘制边坡模型水平方向位移云图 NC@OmSR\0
plnsol,eppl,eqv !绘制边坡模型塑性应变云图
