(1)创建物理环境 Di=9mHC
/COM, Structural qJ8-9^E,L
/TITLE, Slope stability Analysis oP,9#FC|(
/FILNAM,support,1 t7F.[uWD
(2)建立模型 `_ (~ Ud
!进入前处理器 > %*B`oqo
/PREP7 Vm8D"I5i
!定义单元类型 Mth6-^g5
ET,1,PLANE82 !定义边坡围岩单元 dL;HV8z^
KEYOPT,1,3,2 !求解类型为平面应变 FN
)d1q(~
!定义材料属性 kJ=L2g>W<.
!边坡围岩1材料属性 3gfimD$_E
MP,EX,1,3E10 ~U}Mv{y
MP,PRXY,1,0.25 noA-)
MP,DENS,1,2500 !定义塑性模型 Ie'P#e'
TB,DP,1 !定义数据表1 X;fy\HaU
TBDATA,1,9E5,42.8 !给数据表输入数据 QLbMPS
!边坡围岩2材料属性 !采用弹性模型 @qK<T
MP,EX,2,3.2E10 !定义弹性模量 6~5$s1Yc
MP,PRXY,2,0.24 !定义泊松比 ARL
MP,DENS,2,2700 !定义密度 `1p 8C%
!取强度折减系数F=1.2时边坡围岩材料属性 tfiqr|z
MP,EX,3,3E10 Rt=
X%[YL
MP,PRXY,3,0.25 hSqMaX%G
MP,DENS,3,2500 2HOe__Ns
TB,DP,3 !定义数据表3 's@MQ!
*
TBDATA,1,7.5E5,37.7 !给数据表输入数据 9 Aivf+
!取强度折减系数F=1.4时边坡围岩材料属性 o;J;*~g
MP,EX,4,3E10 [{F%LRCo-
MP,PRXY,4,0.25 !定义泊松比 %!.M~5mCd
MP,DENS,4,2500 t6u-G+}
TB,DP,4 !定义数据表4 ~v:#zU
TBDATA,1,6.4E5,33.5 !给数据表输入数据 {^&@gkYY
!取强度折减系数F=1.6时边坡围岩材料属性 pbB2wt
MP,EX,5,3E10 &v#`t~
MP,PRXY,5,0.25 :d'65KMi
MP,DENS,5,2500 !定义密度 K&pM o.
TB,DP,5 !定义数据表5 dc^Vc{26Z
TBDATA,1,5.6E5,30 !给数据表输入数据 izt^Wi|
!取强度折减系数F=1.8时边坡围岩材料属性 9NIy#
MP,EX,6,3E10 !定义弹性模量 rFXSO=P?Z
MP,PRXY,6,0.25 {-*\w-~G
MP,DENS,6,2500 c%<2z
TB,DP,6 !定义数据表6 IUhp;iH
TBDATA,1,5.0E5,27.2 Ao`_",E
!取强度折减系数F=2.0时边坡围岩材料属性 b>q6:=((
MP,EX,7,3E10 ]XrE
MP,PRXY,7,0.25 !定义泊松比 6$B'Q30}r
MP,DENS,7,2500 LZ&uj{ <
TB,DP,7 !定义数据表7 ha'qIT3&
TBDATA,1,4.5E5,24.8 2uu[52H8d%
!取强度折减系数F=2.2时边坡围岩材料属性 kfpm=dKL
MP,EX,8,3E10 %yw=[]Vjze
MP,PRXY,8,0.25 8[\79|
MP,DENS,8,2500 ]Ti$ztJ
TB,DP,8 !定义数据表8 sX'U|)/pD
TBDATA,1,4.09E5,22.8 1*R_"#
!取强度折减系数F=2.4时边坡围岩材料属性 \j:gr>4
MP,EX,9,3E10 E\e]K
!
MP,PRXY,9,0.25 ATO
5
MP,DENS,9,2500 sC6r.@[u8t
TB,DP,9 !定义数据表9 Z>{*ISvpq
TBDATA,1,3.6E5,21.1 !给数据表输入数据 x*mc -&N
!取强度折减系数F=2.6时边坡围岩材料属性 }|He?[TR
MP,EX,10,3E10 !定义弹性模量 ib50LCm
MP,PRXY,10,0.25 !定义泊松比 <kc]L x
MP,DENS,10,2500 !定义密度 0_V*B[V
TB,DP,10 !定义数据表10 u[`v&e
TBDATA,1,3.46E5,19.6 !给数据表10输入数据 iwz`
x
!取强度折减系数F=2.8时边坡围岩材料属性 @aB9%An1
MP,EX,11,3E10 !定义弹性模量 j:?N!*r=
MP,PRXY,11,0.25 !定义泊松比 `!kL1oUYE
MP,DENS,11,2500 Fg@ ACv'@
TB,DP,11 !定义数据表11 3Wj,}
TBDATA,1,3.2E5,18.2 !给数据表11输入数据 wF;B@
!取强度折减系数F=3.0时边坡围岩材料属性 U(A4v0T
MP,EX,12,3E10 !定义弹性模量 9 x [X<
MP,PRXY,12,0.25 LD@7(?mlU
MP,DENS,12,2500 !定义密度 7ti<
TB,DP,12 !定义数据表12 CveWl$T12
TBDATA,1,3.0E5,17.1 !给数据表12输入数据 /Hk07:"c
!建立几何模型 1nXqi)&?;
!创建边坡线模型线 {_ 6t4h}
K,1,,,, ! 创建关键点1 QJM(UfHUD
K,2,-800,,, (wlfMiO
K,3,-800,-800,, z\h,SX<U
K,4,-800,-1200,, W8uVd zQ
K,5,1200,-1200,, %QE5<2k
K,6,1200,-800,, {3SdX
K,7,1200,,, {fElto
K,8,1200,378,, )v-Cj_W5]"
K,9,430,378,, ;Bnr='[
l,1,2 !连接关键点1、2生成直线L1 x?>!UqgkY
l,2,3 P7Z<0Dt\}
l,3,4 o~}1oN
l,4,5 yr{5Rp05=
l,5,6 ?k
w/S4
l,6,7 (l;C%O7*
l,7,8 YZ{jP?x
l,8,9 \vs%U}IrO
l,9,1 !SN WB
l,1,7 u
mqKFM$
l,3,6 wV
%8v\
!创建边坡面模型 t_Z _!Qy
al,3,4,5,11 !由线L3、L4、L5、L11生成一个面A1 >~>{;Wq(p+
al,1,2,11,6,10 !由线L1、L2、L1、L6、L10生成一个面A2 -}AE\qXs/
al,7,8,9,10 !由线L7、L8、L9、L10生成一个面A3 Ku&*`dME
/PNUM, area, on !打开面号开关 {SHqW5VX
!划分网格生成有限元模型
$Gd5wmb!
! 划分边坡围岩2网格 iZu:uMoc
!设置网格份数 o>$|SU!a
lsel,s,,,3,5,2 !选择线L3和L5 8q{1E];:q
lesize,all,,,5 !把所选择线分为5段 ${CYDD"mdy
lsel,s,,,4,11,7 !选择线L4和L11 %,Q;<axzi
lesize,all,,,26 !把所选择线分为26段 Yg|l?d"
mat,2 !给边坡围岩2赋予2号材料特性 mj,qQ=n;p
type,1 !采用单元类型1 kYTOldfY2
mshkey,0 !设定自由网格划分 E.U0qK],
mshape,0 !设定四边形网格划分 XzlIW&"uC
amesh,1 !划分面积A1 ^h"n03VFA
! 划分边坡围岩1网格 ->Q`'@'|P
!设置网格份数 )MMhlcNC
lsel,s,,,2,6,4 !选择线L2和L6 (,d/JnP
lsel,a,,,1 !附加选择线L1 JgxA^>|9;
lesize,all,,,10 !把所选择线分为10段 VEr 6uvB
lsel,s,,,8,10,2 !选择线L8和L10 j&
<tdORT
lesize,all,,,16 !把所选择线分为16段 d{iL?>'?^
lsel,s,,,7,9,2 +H?<}N*T
lesize,all,,,12
}Olr
mat,1 !给边坡围岩1赋予1号材料特性 Qlf
9]ug)
type,1 g8rp|MOH
mshkey,0 !设定自由网格划分 Kyyih|{
mshape,0 !设定四边形网格划分 Sn+FV+D
amesh,2 !划分面积A2 }^IwQm*i
amesh,3 !划分面积A3 f>?^uSpWH
allsel IMw
"eV
!保存网格模型 dp33z"<3
save,Slope-grid.db 5_PWGaQa
(3)施加约束和荷载 nP5d?
!给边坡模型施加约束 //6^+-he
!边坡两侧施加X方向约束 zL6
\p)y
nsel,s,loc,x,-800 !选择X=-800线上所有节点 y`\mQ48V
nsel,a,loc,x,1200 !选择X=1200线上所有节点 }ty"fI3&iY
d,all,ux !对所选择节点约束X方向位移 Vx}Yl&*D
allsel A>J1B(up
!边坡底部施加约束 LAizx^F
nsel,s,loc,y,-1200 !选择Y=--1200线上所有节点 [}jj<!9A_;
d,all,uy !对所选择节点约束X、Y方向位移 !~ BZHi6\
d,all,ux 3^j~~"2,w
!施加重力加速度 3"f)*w7d
acel,,9.8 V^9$t/c&
(4)求解 'MSEki67
/solu A8#.1uEgNb
!求解设置 /0Rt+`
antype,static !设定为静力求解 (QA-"9v#i,
nsubst,100 !设定最大子步数为100 Y1m}@k,+M
pred,on !打开时间步长预测器 >a?OXqYP
nropt,full !设定牛顿-拉普森选项 J2)-cY5G
nlgeom,on !打开大位移效果 Wk0>1 rlu
lnsrch,on !打开线性搜索 x:=0.l#
outres,all,all !输出所有项 AlAh
S<
cnvtol,f,,0.005,2,0.5 !力收敛准则设定 xI-=tib
cnvtol,u,,0.05,2,1 !位移收敛准则设定 FGV}5L
!边坡在强度折减系数F=1时求解 ',L{CQA?c
allsel s$js5
ou
solve !进行求解 k,
$I59
save,F1,db !把F=1时求解结果保存 97['VOh0
!边坡在强度折减系数F=1.2时求解 J(3gT}z-
finish k'6<jEbk
/solu Fl8w7LcF7
allsel 2]?w~qjWm
asel,s,area,,2,3,1 !选择面积A2、A3 / c4;3>IS
MPCHG,3,all ! 把所选择单元材料号改为3 HVtr,jg
allsel iMP*]K-O
solve !进行求解 |LXrGyk^
save,F1.2,db !把F=1.2时求解结果保存 Ufm(2`FQ
!边坡在强度折减系数F=1.4时求解 bbE bf !E
finish KyuA5jQ7
/solu 4.,KEt'H
allsel g,A.Y,})
asel,s,area,,2,3,1 !选择面积A2、A3 [K"U_b}w
MPCHG,4,all !把所选择单元材料号改为4 DBqg_v
allsel ~E^yM=:h
solve !进行求解 ckH$E%j
save,F1.4,db !把F=1.4时求解结果保存 KK&<Vw|O\
!边坡在强度折减系数F=1.6时求解 [Ihp\!xqI
finish I}6DoLbV
/solu xn%l
allsel Qx6,>'Qk'
asel,s,area,,2,3,1 !选择面积A2、A3 /}h71V!
MPCHG,5,all !把所选择单元材料号改为5 "R@$Wu53|
allsel m_{%tU;N
solve BFzcoBu-
save,F1.6,db !把F=1.6时求解结果保存 $[HcHnf
!边坡在强度折减系数F=1.8时求解 5`1(}
finish f_Wkg)g
/solu +YGw4{\EL
allsel [HC8-N^.}
asel,s,area,,2,3,1 !选择面积A2、A3 6Tm
Rc
MPCHG,6,all !把所选择单元材料号改为6 \;3B?8wbIl
allsel z5|e\Z
solve Pg!;o=
{M
save,F1.8,db !把F=1.8时求解结果保存 n"^/UQ|#j
!边坡在强度折减系数F=2.0时求解 h,!G7V
finish h|(ZXCH
/solu e>])m3xvn
allsel WHpUjyBP
asel,s,area,,2,3,1 !选择面积A2、A3 PK:o}IWn~x
MPCHG,7,all !把所选择单元材料号改为7 1q}u?7nnSG
allsel 3{2^G@j
solve r`&2-]
save,F2.0,db !把F=2.0时求解结果保存 vF*^xhh
!边坡在强度折减系数F=2.2时求解 0?J|C6XM#4
finish ? 6yF{!F*
/solu 0)6i~MglY
allsel yV 9]_k
asel,s,area,,2,3,1 !选择面积A2、A3 Z@>=&
MPCHG,8,all !把所选择单元材料号改为8 7G<KrKal
allsel I]uOMWZs
solve +
d+hvwEM
save,F2.2,db !把F=2.2时求解结果保存 5 WN`8?
!边坡在强度折减系数F=2.4时求解 #{)r*"%
finish !I~C\$^U
/solu Z1)jRE2dl
allsel cuV8#:
i
asel,s,area,,2,3,1 !选择面积A2、A3 F#!@}K8
MPCHG,9,all !把所选择单元材料号改为9 =|qt!gY)Y
allsel LK
solve ?$)x$nS`
save,F2.4,db !把F=2.4时求解结果保存 Tc'{i#%9j
!边坡在强度折减系数F=2.6时求解 T!^?d5uW#
finish Vid{6?7kh
/solu tdw\Di#m
allsel E1U4v&P
asel,s,area,,2,3,1 !选择面积A2、A3 yL.PGF1(
MPCHG,10,all !把所选择单元材料号改为10 -H ac^4uF
allsel EMVoTW)z
solve =ELDJt
save,F2.6,db !把F=2.6时求解结果保存 *MnG-\{j
!边坡在强度折减系数F=2.8时求解 D^N#E>,
finish K#JabT
/solu Cu
['&_@
allsel dIBKE0`
asel,s,area,,2,3,1 !选择面积A2、A3 cKi^C
MPCHG,11,all !把所选择单元材料号改为11 p,[XT`q^
allsel E&s'uE=w+
solve 4BduUH
save,F2.8,db !把F=2.8时求解结果保存 =x='<{jtgW
!边坡在强度折减系数F=3.0时求解 y'0dl "Dy\
finish @~!-a
s7
/solu iSZctsqE
allsel -A-hxK*^
asel,s,area,,2,3,1 !选择面积A2、A3 OUIUgej
MPCHG,12,all !把所选择单元材料号改为12 m! '1$G
allsel %X0NHta~@
solve
Z}'F"}QI
save,F3.0,db !把F=3.0时求解结果保存 d0H
/post1 !进入后处理 Z3abem<Q
!边坡在强度折减系数F=1时结果分析 YP$*;l
Resume,'F1','db' !读入边坡在强度折减系数F=1时 @LWxz
set,1,last !读入后一个子步 ]JqkC4|
pldisp,1 !绘制边坡模型变形图 %0~wtZH_!
plnsol,u,x !绘制边坡模型水平方向位移云图 Q~b M
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 XRz%KVysp
!边坡在强度折减系数F=1.2时结果分析 fbzKO^Ub
Resume,'F1.2','db' !读入边坡在强度折减系数F=1.2时 UpszCY4
set,1,last !读入后一个子步 Hl3XqR
pldisp,1 !绘制边坡模型变形图 j
J`Zz
plnsol,u,x !绘制边坡模型水平方向位移云图 C\a:eSgaC
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 53,,%Ue
!边坡在强度折减系数F=1.4时结果分析 k8x&aH
Resume,'F1.4','db' !读入边坡在强度折减系数F=1.4时 d=4f`q0k
set,1,last !读入后一个子步 8~[C'+r
pldisp,1 !绘制边坡模型变形图 syC"eH3{
plnsol,u,x !绘制边坡模型水平方向位移云图 2l[A=Z
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 iw~V_y4
!边坡在强度折减系数F=1.6时结果分析 /_VRO9R\V
Resume,'F1.6','db' !读入边坡在强度折减系数F=1.6时 Y#SmZ*zok
set,1,last !读入后一个子步 'wB Huq
pldisp,1 !绘制边坡模型变形图 K9I,Q$&xX
plnsol,u,x !绘制边坡模型水平方向位移云图 ot>EnHfV
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 eMm~7\
R
!边坡在强度折减系数F=1.8时结果分析 U$/Hp#~X
Resume,'F1.8','db' !读入边坡在强度折减系数F=1.8时 Kt4\&l-De
set,1,last !读入后一个子步 z:i X]df
pldisp,1 !绘制边坡模型变形图 w
/W
Cj4`
plnsol,u,x !绘制边坡模型水平方向位移云图 fN"oa>X
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 A9qO2kq7_
!边坡在强度折减系数F=2.0时结果分析 Y)4Nydq
Resume,'F2.0','db' !读入边坡在强度折减系数F=2.0时 {Hp}F!X$
set,1,last !读入后一个子步 NBg>i7KQ
pldisp,1 !绘制边坡模型变形图 !6tC[W`
plnsol,u,x !绘制边坡模型水平方向位移云图 8SCW.;0
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 PkCeV]`w
!边坡在强度折减系数F=2.2时结果分析 Zs5I?R1e8
Resume,'F2.2','db' !读入边坡在强度折减系数F=2.2时 CI~;B
set,1,last !读入后一个子步 SJ~I
r#
pldisp,1 !绘制边坡模型变形图 =@Nv:1:r
plnsol,u,x !绘制边坡模型水平方向位移云图 ^;k _
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 l5y#i7q
!边坡在强度折减系数F=2.4时结果分析 DKVt8/vq
Resume,'F2.4','db' !读入边坡在强度折减系数F=2.4时 {DXZ}7w:v
set,1,last !读入后一个子步 H-cBXp5z
pldisp,1 !绘制边坡模型变形图 R
!%m5Q?5
plnsol,u,x !绘制边坡模型水平方向位移云图 ?k:])^G5
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 Er/5 ,
!边坡在强度折减系数F=2.6时结果分析 'd.@4 9
Resume,'F2.6','db' !读入边坡在强度折减系数F=2.6时
oRbYna?J
set,1,last !读入后一个子步 }DUDA%U
pldisp,1 !绘制边坡模型变形图 j]?0}Z*
plnsol,u,x !绘制边坡模型水平方向位移云图 PRk%C0`
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 ^; V>}08
!边坡在强度折减系数F=2.8时结果分析 4h_4jqf=pU
Resume,'F2.8','db' !读入边坡在强度折减系数F=2.8时 CF}Nom)
set,1,last !读入后一个子步 7f\^VG
pldisp,1 !绘制边坡模型变形图 zloaU
plnsol,u,x !绘制边坡模型水平方向位移云图 J2rLsNC]0
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 =<'iLQb1
!边坡在强度折减系数F=3.0时结果分析 f`9rTc
Resume,'F3.0','db' !读入边坡在强度折减系数F=3.0时 -SY:qG3?
set,1,last !读入后一个子步 w[A3;]la
pldisp,1 !绘制边坡模型变形图 rIFC#Jd/
plnsol,u,x !绘制边坡模型水平方向位移云图 /6+%(f}7l
plnsol,eppl,eqv !绘制边坡模型塑性应变云图