(1)创建物理环境 Z^6#4Q]YC
/COM, Structural b& t[S[P.V
/TITLE, Slope stability Analysis f~gSJ< t4
/FILNAM,support,1 t3F?>G#y
(2)建立模型 LdTIR ]
!进入前处理器 c;V D}UD'
/PREP7 $ AwZ2HY
!定义单元类型 LF=c^9t
ET,1,PLANE82 !定义边坡围岩单元 Kmry=`=A
KEYOPT,1,3,2 !求解类型为平面应变 N~#D\X^t.
!定义材料属性 {OxWcK\2@h
!边坡围岩1材料属性 E [S' :Q
MP,EX,1,3E10 L{PH0Jf
MP,PRXY,1,0.25 Gp1EJ2d8
MP,DENS,1,2500 !定义塑性模型 !UNNjBBP7
TB,DP,1 !定义数据表1 KR k~w]
TBDATA,1,9E5,42.8 !给数据表输入数据 $pu3Ig$^
!边坡围岩2材料属性 !采用弹性模型 p{H0dj ^|
MP,EX,2,3.2E10 ！定义弹性模量 qC
j*> D
MP,PRXY,2,0.24 !定义泊松比 !k h{9I>M
MP,DENS,2,2700 !定义密度 Qb@eK$wo}
!取强度折减系数F=1.2时边坡围岩材料属性 wu^q`!ml
MP,EX,3,3E10 ~/ c5 hyTx
MP,PRXY,3,0.25 fPE ? hG<x
MP,DENS,3,2500 (|0b7 |'T
TB,DP,3 !定义数据表3 NW ISS
TBDATA,1,7.5E5,37.7 !给数据表输入数据 5#u.pu
!取强度折减系数F=1.4时边坡围岩材料属性 kMS5h~D[
MP,EX,4,3E10 BvJ=iB<E
MP,PRXY,4,0.25 !定义泊松比 FGVb@=TO>
MP,DENS,4,2500 y&rY0bm
TB,DP,4 !定义数据表4 h DtK nF
TBDATA,1,6.4E5,33.5 !给数据表输入数据 wznn #j
!取强度折减系数F=1.6时边坡围岩材料属性 @&:VKpu\
MP,EX,5,3E10 5'9.np F)
MP,PRXY,5,0.25 ?'+8[OHiF^
MP,DENS,5,2500 ！定义密度 DcU C,
TB,DP,5 !定义数据表5 A}4t9|/K6
TBDATA,1,5.6E5,30 !给数据表输入数据 7|~:P $M
!取强度折减系数F=1.8时边坡围岩材料属性 Y(z }[`2
MP,EX,6,3E10 ！定义弹性模量 twv
lQ|
MP,PRXY,6,0.25 w b+<a
MP,DENS,6,2500 8nu> gA
TB,DP,6 !定义数据表6 Ki 3_N*z
TBDATA,1,5.0E5,27.2 ^Jtl;Q
!取强度折减系数F=2.0时边坡围岩材料属性 ( (9 YG
MP,EX,7,3E10 Kw*~W
i
MP,PRXY,7,0.25 !定义泊松比 QZ0R :TY
MP,DENS,7,2500 ohTd' +Lm
TB,DP,7 !定义数据表7 |8;?
*s`H
TBDATA,1,4.5E5,24.8 Qdm(q:w
!取强度折减系数F=2.2时边坡围岩材料属性 `r.N
MP,EX,8,3E10 ~h|m&XK+Q
MP,PRXY,8,0.25 ylim/`u}6
MP,DENS,8,2500 vTcZ8|3 e
TB,DP,8 !定义数据表8 Qg{WMlyOP
TBDATA,1,4.09E5,22.8 '_E c_ F
!取强度折减系数F=2.4时边坡围岩材料属性 J(hA^;8:
MP,EX,9,3E10 GWh|FEqUbf
MP,PRXY,9,0.25 #XL`S
MP,DENS,9,2500 /fC\K_<N
TB,DP,9 !定义数据表9 <|Iyt[s
TBDATA,1,3.6E5,21.1 !给数据表输入数据 mrRea st
!取强度折减系数F=2.6时边坡围岩材料属性 pg5&=
MP,EX,10,3E10 ！定义弹性模量 r$?Vx_f`Q
MP,PRXY,10,0.25 !定义泊松比 k) 3s?
MP,DENS,10,2500 ！定义密度 uf?b%:A
TB,DP,10 !定义数据表10 UhA _1A'B
TBDATA,1,3.46E5,19.6 !给数据表10输入数据 500>
CBL0O
!取强度折减系数F=2.8时边坡围岩材料属性 j #es2;
MP,EX,11,3E10 ！定义弹性模量 M)tv;!eQ
MP,PRXY,11,0.25 !定义泊松比 ,N;v~D$Y
MP,DENS,11,2500 2b~
HHVruX
TB,DP,11 !定义数据表11 N=<=dp(
TBDATA,1,3.2E5,18.2 !给数据表11输入数据 2moIgJ
!取强度折减系数F=3.0时边坡围岩材料属性 G)4SWu0<t
MP,EX,12,3E10 ！定义弹性模量 Wa{%0inZ
MP,PRXY,12,0.25 'M
lXnHxt
MP,DENS,12,2500 ！定义密度
"$Y(NFb
TB,DP,12 !定义数据表12 z^9E;
TBDATA,1,3.0E5,17.1 !给数据表12输入数据 (6o:4|xl0
!建立几何模型 @ oA0{&G{
!创建边坡线模型线 lk N'uZ
K,1,,,, ！ 创建关键点1 Q.>/*8R;
K,2,-800,,, [{/$9k-aF?
K,3,-800,-800,, <YSg~T
K,4,-800,-1200,, ?}v}U^
K,5,1200,-1200,, c-|kv[\a
K,6,1200,-800,, uYY=~o[
Tw
K,7,1200,,, uh1S
7!^
K,8,1200,378,, ]p}#N Pe5
K,9,430,378,, CY5w$E
l,1,2 !连接关键点1、2生成直线L1 >0 7shNX
l,2,3 k|l5 "&K~.
l,3,4 Qv\bLR
l,4,5 !&\meS{
l,5,6 Slo9#26
l,6,7 4"z;CGE7
l,7,8 G6<HO7\
l,8,9 H[Cj7{ V
l,9,1 &6Ns7w6*z
l,1,7 nc- Qz
l,3,6 -#daBx
?
!创建边坡面模型 ;i/"$K
al,3,4,5,11 !由线L3、L4、L5、L11生成一个面A1 hO%Y{Gg
al,1,2,11,6,10 !由线L1、L2、L1、L6、L10生成一个面A2 i5Eeg`NMl
al,7,8,9,10 !由线L7、L8、L9、L10生成一个面A3 MHGj vSx
/PNUM, area, on !打开面号开关 kO jEY
!划分网格生成有限元模型 8KZ$ F>T]>
! 划分边坡围岩2网格 qih 6me8 C
!设置网格份数 ]|t9B/()i
lsel,s,,,3,5,2 !选择线L3和L5 C1n? ?Y[
lesize,all,,,5 !把所选择线分为5段 p AMo
XJ`
lsel,s,,,4,11,7 !选择线L4和L11 aQxe)
lesize,all,,,26 !把所选择线分为26段 g&q^.7c}
mat,2 !给边坡围岩2赋予2号材料特性 6(,ItMbI
type,1 !采用单元类型1 iY}QgB< M
mshkey,0 !设定自由网格划分 >Sh0dFqeT
mshape,0 !设定四边形网格划分 *1i?6$[
"
amesh,1 ！划分面积A1 x
R idc^
! 划分边坡围岩1网格 o}d2N/T
!设置网格份数 }WA<=9e
lsel,s,,,2,6,4 !选择线L2和L6 _dJp
3D
lsel,a,,,1 !附加选择线L1 _d/GdeLs
lesize,all,,,10 !把所选择线分为10段 LPb]mC6#
lsel,s,,,8,10,2 !选择线L8和L10 j1W
bD7*8
lesize,all,,,16 !把所选择线分为16段 <` #,AVH
lsel,s,,,7,9,2 M PqY?KF
lesize,all,,,12 *h >OW
mat,1 !给边坡围岩1赋予1号材料特性 N&x@_t""
type,1 H&w(]PDh
mshkey,0 !设定自由网格划分 zI1(F67d`
mshape,0 !设定四边形网格划分 %f_FGh
amesh,2 !划分面积A2 is64)2F](
amesh,3 ！划分面积A3 5sG ] 3z+1
allsel U`:l AG
!保存网格模型 z w K g
save,Slope-grid.db "$# $f
(3)施加约束和荷载 (H'_KPK
!给边坡模型施加约束 5%EaX ?0h+
!边坡两侧施加X方向约束 ="/R5fp
nsel,s,loc,x,-800 !选择X=-800线上所有节点 $_ST:h&C
nsel,a,loc,x,1200 !选择X=1200线上所有节点 \MfR #k0
d,all,ux !对所选择节点约束X方向位移 Dm&lSWW`/
allsel M0 Lon/%
!边坡底部施加约束 3S;>ki4(0
nsel,s,loc,y,-1200 !选择Y=--1200线上所有节点 yjeqv-7
d,all,uy !对所选择节点约束X、Y方向位移 \ x3^
d,all,ux =ltbS f7
!施加重力加速度 n iXHK$@5
acel,,9.8 .WxFm@]/\
(4)求解 =c
:lS&B
/solu "K9vm^xP
!求解设置 Uyz;U34 oI
antype,static !设定为静力求解 $0S.@wUG
nsubst,100 !设定最大子步数为100 B.h0" vJ
pred,on !打开时间步长预测器 -HO6K) ur
nropt,full !设定牛顿-拉普森选项 }w;Q^EU
nlgeom,on !打开大位移效果 \Mobq
lnsrch,on !打开线性搜索 Uh<H*o6e 9
outres,all,all !输出所有项 cP/F| uG5
cnvtol,f,,0.005,2,0.5 !力收敛准则设定 AaJ,=eQ
cnvtol,u,,0.05,2,1 !位移收敛准则设定 |:!E HFr
!边坡在强度折减系数F=1时求解 Tb>IHoil
allsel 9{au leu
R
solve !进行求解 *?EjYI
save,F1,db !把F=1时求解结果保存 v'mJ~tz
!边坡在强度折减系数F=1.2时求解 %#[r_QQ^
finish CUAg{]
/solu gt] . rwo"
allsel pZni,< Q
asel,s,area,,2,3,1 !选择面积A2、A3 AJJ%gxqGq
MPCHG,3,all ! 把所选择单元材料号改为3 | 7t= \
allsel nPQZI6>
solve !进行求解 (ouRf;\6$8
save,F1.2,db !把F=1.2时求解结果保存 V\e13cL]
!边坡在强度折减系数F=1.4时求解 ' lQ
finish oyi7Y Rvwd
/solu RZcx4fL }x
allsel T%N~oa
asel,s,area,,2,3,1 !选择面积A2、A3 CDwFVR'_Af
MPCHG,4,all !把所选择单元材料号改为4 f(@"[-[
allsel 7]<F>97
solve !进行求解 d()zW7}W
save,F1.4,db !把F=1.4时求解结果保存 !W+p<F1i
!边坡在强度折减系数F=1.6时求解 .3S\Rrv
finish GN~[xX JU
/solu C[Y%=\6'0
allsel Z7m GC`>
asel,s,area,,2,3,1 !选择面积A2、A3 >*<6 zQf
MPCHG,5,all !把所选择单元材料号改为5 a:(: :m
allsel s$2l"|h>B
solve Q]2sj:
save,F1.6,db !把F=1.6时求解结果保存 UH1S_:6
!边坡在强度折减系数F=1.8时求解 "4Wp>B
finish URmAI8fq*M
/solu rU2YMg hE
allsel }uV?
asel,s,area,,2,3,1 !选择面积A2、A3 R-2F Nl
MPCHG,6,all !把所选择单元材料号改为6 KV) Hywl`
allsel 5kRwSOG%'
solve YokZar2a0
save,F1.8,db !把F=1.8时求解结果保存 x):k#cu[L
!边坡在强度折减系数F=2.0时求解 E'Fv *UA
finish 1OfSq1G>v$
/solu 3<~2"@J
allsel E9 Qd>o
asel,s,area,,2,3,1 !选择面积A2、A3 Jm*M7g j
MPCHG,7,all !把所选择单元材料号改为7 <?=mLOo =
allsel +*r** (-Dm
solve ~_# Y,)S!z
save,F2.0,db !把F=2.0时求解结果保存 %Mb(
c+7
!边坡在强度折减系数F=2.2时求解 y . ivz
finish FJwZo}<6E
/solu }
K hq
allsel K@$L~G
asel,s,area,,2,3,1 !选择面积A2、A3 XeIUdg4>R
MPCHG,8,all !把所选择单元材料号改为8 Zh fD`@>&
allsel J 4*:.8Ki
solve >TY5ZRB
save,F2.2,db !把F=2.2时求解结果保存 &ttv4BC^r
!边坡在强度折减系数F=2.4时求解 _L `N^I.
finish 9wYtOQ{g
/solu F`ZIc7(.{
allsel "Wn?8vR
asel,s,area,,2,3,1 !选择面积A2、A3 w $U/;C
MPCHG,9,all !把所选择单元材料号改为9 #)h
~ .D{
allsel U_
?elz\
solve yBD2
save,F2.4,db !把F=2.4时求解结果保存 =p&'_a^$
!边坡在强度折减系数F=2.6时求解 DESV iQM
finish #pm-nU%|_j
/solu GJ^]ER-K
allsel mC[U)` ey
asel,s,area,,2,3,1 !选择面积A2、A3 9,EaN{GM
MPCHG,10,all !把所选择单元材料号改为10 Uf2v$Jl+Yh
allsel ,GXfy9x7U
solve SJ W P8+
save,F2.6,db !把F=2.6时求解结果保存 |au q j2
!边坡在强度折减系数F=2.8时求解 M_e$l`"G
finish [K4+G]6
/solu rO?x/{;ai
allsel ~&j`9jdO j
asel,s,area,,2,3,1 !选择面积A2、A3 2dd:5L,
MPCHG,11,all !把所选择单元材料号改为11 BHJ'[{U*w
allsel 0@FM^ejA#
solve :C:N]6_{SZ
save,F2.8,db !把F=2.8时求解结果保存 BtNW5'^
!边坡在强度折减系数F=3.0时求解 ,b.4uJg'
finish gt t$O
/solu C +?@iMh
allsel 4gkV]"
H!
asel,s,area,,2,3,1 !选择面积A2、A3 s
eZ<52f2
MPCHG,12,all !把所选择单元材料号改为12 Vw;ldEdx
allsel |C)UZ4A/p
solve <K=B(-~
save,F3.0,db !把F=3.0时求解结果保存 &fd4IO/O
/post1 ！进入后处理 64 \5 v?C
!边坡在强度折减系数F=1时结果分析 :rs\ydDUF
Resume,'F1','db' !读入边坡在强度折减系数F=1时 Pdm6u73
set,1,last !读入后一个子步 yL3<X w|
pldisp,1 !绘制边坡模型变形图 ? "8A^
^
plnsol,u,x !绘制边坡模型水平方向位移云图 h3d\MYO)B
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 }jY[| >z
!边坡在强度折减系数F=1.2时结果分析 R;TEtu7
Resume,'F1.2','db' !读入边坡在强度折减系数F=1.2时 M}oFn}-T9a
set,1,last !读入后一个子步 n-b<vEZw#
pldisp,1 !绘制边坡模型变形图 ">vYEkZ3
plnsol,u,x !绘制边坡模型水平方向位移云图 Z/05 wB
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 hte9l)
!边坡在强度折减系数F=1.4时结果分析 0-~6}
r$
Resume,'F1.4','db' !读入边坡在强度折减系数F=1.4时 G>ptwB81KM
set,1,last !读入后一个子步 mv%:[+!
pldisp,1 !绘制边坡模型变形图 k]/6 /s\
plnsol,u,x !绘制边坡模型水平方向位移云图 })0 7u
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 ,nChwEn
!边坡在强度折减系数F=1.6时结果分析 7" STS7_
Resume,'F1.6','db' !读入边坡在强度折减系数F=1.6时 E#2k|TpH4
set,1,last !读入后一个子步 GWqY$YT
pldisp,1 !绘制边坡模型变形图 J_[[BJ&}x
plnsol,u,x !绘制边坡模型水平方向位移云图 ?5YmE(v7
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 yDyeP{
!边坡在强度折减系数F=1.8时结果分析 1DlcO>#@
Resume,'F1.8','db' !读入边坡在强度折减系数F=1.8时 hHl-;%#
set,1,last !读入后一个子步 vB YT)S
pldisp,1 !绘制边坡模型变形图 UrcN?
plnsol,u,x !绘制边坡模型水平方向位移云图 >'TD?@sr
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 c'#J{3d
!边坡在强度折减系数F=2.0时结果分析 "`3 ^M vC
Resume,'F2.0','db' !读入边坡在强度折减系数F=2.0时 s^?sJUj
set,1,last !读入后一个子步 ^YJ^+:D(
pldisp,1 !绘制边坡模型变形图 @oc%4~zl
plnsol,u,x !绘制边坡模型水平方向位移云图 2SPFjpG8n
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 )4_6\VaM
!边坡在强度折减系数F=2.2时结果分析 rWfurB5f
Resume,'F2.2','db' !读入边坡在强度折减系数F=2.2时 YG1`%,OW`
set,1,last !读入后一个子步 '-] BSU
pldisp,1 !绘制边坡模型变形图 (DMnwqr
plnsol,u,x !绘制边坡模型水平方向位移云图 2l(j
4~ g
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 vgW1hWmHJ
!边坡在强度折减系数F=2.4时结果分析 >.wd)
Resume,'F2.4','db' !读入边坡在强度折减系数F=2.4时 y3[)zv
set,1,last !读入后一个子步 [O+^e E6h
pldisp,1 !绘制边坡模型变形图 7C?mD75j
plnsol,u,x !绘制边坡模型水平方向位移云图 a={qA4N
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 $~@096`QL<
!边坡在强度折减系数F=2.6时结果分析 - QI`npsnV
Resume,'F2.6','db' !读入边坡在强度折减系数F=2.6时 IS
9q 5/]
set,1,last !读入后一个子步 v?TJ!o
pldisp,1 !绘制边坡模型变形图 wYMX1=
plnsol,u,x !绘制边坡模型水平方向位移云图 L8tLW09
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 y'O{8 Q8T
!边坡在强度折减系数F=2.8时结果分析 W>wi;Gf#
Resume,'F2.8','db' !读入边坡在强度折减系数F=2.8时 tMBy
^@p
set,1,last !读入后一个子步 2T%f~yQ^
pldisp,1 !绘制边坡模型变形图 ,Ve@=<
plnsol,u,x !绘制边坡模型水平方向位移云图 RrpF i'R
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 HWfX>Vf>}k
!边坡在强度折减系数F=3.0时结果分析 , Y8X"~{A
Resume,'F3.0','db' !读入边坡在强度折减系数F=3.0时 jsaCnm>&
set,1,last !读入后一个子步 wVJFA1
pldisp,1 !绘制边坡模型变形图 J'c9577$
plnsol,u,x !绘制边坡模型水平方向位移云图 yIf}b
plnsol,eppl,eqv !绘制边坡模型塑性应变云图
