(1)创建物理环境 xqG[~)~
/COM, Structural ;RC{<wBTx
/TITLE, Slope stability Analysis )f[
B6Y
/FILNAM,support,1 :*Z4yx
(2)建立模型 = C8 ?M
!进入前处理器 4gz
H8sF
/PREP7 L\xk:j1[
!定义单元类型 "!ug_'VW
ET,1,PLANE82 !定义边坡围岩单元 Ez
fN&8E
KEYOPT,1,3,2 !求解类型为平面应变 [6 %VR qY
!定义材料属性 v4`"1Ss,K
!边坡围岩1材料属性 ^cP!\E-^
MP,EX,1,3E10 AQ ,'
6F9
MP,PRXY,1,0.25 Mb|a+,:>3
MP,DENS,1,2500 !定义塑性模型 [S9K6%w_!
TB,DP,1 !定义数据表1 :toh0oB[
TBDATA,1,9E5,42.8 !给数据表输入数据 ;5S9y7[i|
!边坡围岩2材料属性 !采用弹性模型 p93r'&Q
MP,EX,2,3.2E10 ！定义弹性模量 =J K@z
MP,PRXY,2,0.24 !定义泊松比 ^6 MU
0Q2
MP,DENS,2,2700 !定义密度 g9}DnCT*.
!取强度折减系数F=1.2时边坡围岩材料属性 t4zkt!`B
MP,EX,3,3E10 7XTkX"zKj
MP,PRXY,3,0.25 9=8iy
w
MP,DENS,3,2500 8hOk{xs8
TB,DP,3 !定义数据表3 i1NY9br
TBDATA,1,7.5E5,37.7 !给数据表输入数据 t(NI-UXBp
!取强度折减系数F=1.4时边坡围岩材料属性 D%OQ e#!
MP,EX,4,3E10 4U+xb>
MP,PRXY,4,0.25 !定义泊松比 r%yvOF\>
MP,DENS,4,2500 7vrl'^ 1
TB,DP,4 !定义数据表4 { GK;63`1
TBDATA,1,6.4E5,33.5 !给数据表输入数据 P3x= 8_#
!取强度折减系数F=1.6时边坡围岩材料属性 j<V Fn~*_
MP,EX,5,3E10 '
V^6XI
MP,PRXY,5,0.25 _VRpI)mu
MP,DENS,5,2500 ！定义密度 Q
Nh|Wz
TB,DP,5 !定义数据表5 Vt %bI0#
TBDATA,1,5.6E5,30 !给数据表输入数据 1{ ~#H<K
!取强度折减系数F=1.8时边坡围岩材料属性 l1OE!W W
MP,EX,6,3E10 ！定义弹性模量 p.v0D:@&
MP,PRXY,6,0.25 P2BWuh F
MP,DENS,6,2500 ao1(]64X"
TB,DP,6 !定义数据表6 |,gc_G
TBDATA,1,5.0E5,27.2 8 *#R]9
!取强度折减系数F=2.0时边坡围岩材料属性 2Mc3|T4)U
MP,EX,7,3E10 j,lT>/
MP,PRXY,7,0.25 !定义泊松比 M "p
MP,DENS,7,2500 A@1W}8qY:
TB,DP,7 !定义数据表7 ;=eDO (Ij
TBDATA,1,4.5E5,24.8 bLi j7K 2H
!取强度折减系数F=2.2时边坡围岩材料属性 dJeN bVd
MP,EX,8,3E10 7Bzq,2s
MP,PRXY,8,0.25 VM
GS[qrG
MP,DENS,8,2500 pfA|I*`XV
TB,DP,8 !定义数据表8
- D
TBDATA,1,4.09E5,22.8 v &Yi
!取强度折减系数F=2.4时边坡围岩材料属性 !;Yg/'vD-
MP,EX,9,3E10 fk6%XO
MP,PRXY,9,0.25 cl=EA6P\X
MP,DENS,9,2500 A+ZK4]xb
TB,DP,9 !定义数据表9 [!HEQ8 2g
TBDATA,1,3.6E5,21.1 !给数据表输入数据 la0BiLzb]
!取强度折减系数F=2.6时边坡围岩材料属性 "GMBjT8
MP,EX,10,3E10 ！定义弹性模量 AN8`7F1
MP,PRXY,10,0.25 !定义泊松比 B%)%
MP,DENS,10,2500 ！定义密度 |:nOp(A\*
TB,DP,10 !定义数据表10 O`x;, 6Vr
TBDATA,1,3.46E5,19.6 !给数据表10输入数据 5cL83F Qh
!取强度折减系数F=2.8时边坡围岩材料属性 1PVtxL?1P
MP,EX,11,3E10 ！定义弹性模量 1 d}Z(My
MP,PRXY,11,0.25 !定义泊松比 zZPuha8
MP,DENS,11,2500 p*4':TFuD;
TB,DP,11 !定义数据表11 e6R}0w~G
TBDATA,1,3.2E5,18.2 !给数据表11输入数据 ]@j*/IP
!取强度折减系数F=3.0时边坡围岩材料属性 _~IR6dKE
MP,EX,12,3E10 ！定义弹性模量 %Gz0^[+
MP,PRXY,12,0.25 B( LWdap~
MP,DENS,12,2500 ！定义密度 )t0$qd ]
TB,DP,12 !定义数据表12 ~:kZgUP_f
TBDATA,1,3.0E5,17.1 !给数据表12输入数据 {5U{ 8b]k
!建立几何模型 42{Ew8
!创建边坡线模型线 o{* e'4
K,1,,,, ！ 创建关键点1 m ZtCL
K,2,-800,,, QdH\LL^8R4
K,3,-800,-800,, ZRh~`yy
K,4,-800,-1200,, V:In>u$QJ!
K,5,1200,-1200,, 5[k/s}g
K,6,1200,-800,, );
!eow
K,7,1200,,, Xx. "$l
K,8,1200,378,, 9Y IM'q>`v
K,9,430,378,, M -cTR d-i
l,1,2 !连接关键点1、2生成直线L1 :~e>Ob[,"
l,2,3 ww\C Q6/h
l,3,4 Neq+16*u
l,4,5 l&OKBUG
l,5,6 D/Z6C&/I
l,6,7 y~ AVei&
l,7,8 X$
0 ?j 1
l,8,9 VRWAm>u
l,9,1 u]<,,
l,1,7 f HE <(
l,3,6 OE_XCZ!5P
!创建边坡面模型 *}F3M\
al,3,4,5,11 !由线L3、L4、L5、L11生成一个面A1 S!jTyY7e
al,1,2,11,6,10 !由线L1、L2、L1、L6、L10生成一个面A2 z1PBMSG
al,7,8,9,10 !由线L7、L8、L9、L10生成一个面A3 p4. wh|n
/PNUM, area, on !打开面号开关 -LK
B$
!划分网格生成有限元模型 Se :.4<
! 划分边坡围岩2网格 m[? E
!设置网格份数 2,$ 8icM
lsel,s,,,3,5,2 !选择线L3和L5 |o H,
lesize,all,,,5 !把所选择线分为5段 >kj`7GA
lsel,s,,,4,11,7 !选择线L4和L11 #%a;"w
lesize,all,,,26 !把所选择线分为26段 qON|4+~u%
mat,2 !给边坡围岩2赋予2号材料特性 D.B.7-_8
type,1 !采用单元类型1 R&8Iz
yM
mshkey,0 !设定自由网格划分 s @&`f{
mshape,0 !设定四边形网格划分 R}
eN@#"D
amesh,1 ！划分面积A1 3 }TaF~
! 划分边坡围岩1网格 kO.%9wFbz
!设置网格份数 >Ea8G,
lsel,s,,,2,6,4 !选择线L2和L6 =x%dNf$e{W
lsel,a,,,1 !附加选择线L1 ~
-4{B
lesize,all,,,10 !把所选择线分为10段 nhB1D-
lsel,s,,,8,10,2 !选择线L8和L10 :~b3^xhc^
lesize,all,,,16 !把所选择线分为16段 gp };D
lsel,s,,,7,9,2 Y;dz,}re
lesize,all,,,12 8;b( 0^
mat,1 !给边坡围岩1赋予1号材料特性 2iY3L sna
type,1 GY6`JWk
mshkey,0 !设定自由网格划分 [YRz*5
mshape,0 !设定四边形网格划分 .b3Qfxc>
amesh,2 !划分面积A2 f=(?JT
amesh,3 ！划分面积A3 nrL9
E'F'
allsel q@QksAq
!保存网格模型 /\ y?Y
save,Slope-grid.db Y_;#UU689
(3)施加约束和荷载 ~P*6ozSYpY
!给边坡模型施加约束 s:>Va GC
!边坡两侧施加X方向约束 3m ]4=
nsel,s,loc,x,-800 !选择X=-800线上所有节点 ~("5y G
nsel,a,loc,x,1200 !选择X=1200线上所有节点 bR*-Ht+wd
d,all,ux !对所选择节点约束X方向位移 XX7{-Y y
allsel KyVQh8
!边坡底部施加约束 {@H6HqD
nsel,s,loc,y,-1200 !选择Y=--1200线上所有节点 8 ##-EN;ag
d,all,uy !对所选择节点约束X、Y方向位移 {C 6=[
d,all,ux *g:4e3Iy
!施加重力加速度 Qw Nly4
acel,,9.8 <C1H36p
(4)求解 !O+) sbd<
/solu C]O(T2l{l
!求解设置 mq aHwID
antype,static !设定为静力求解 R kH W
nsubst,100 !设定最大子步数为100 rHC>z7+z.
pred,on !打开时间步长预测器 3c#BKHNC
nropt,full !设定牛顿-拉普森选项 -[ *,^Ti`
nlgeom,on !打开大位移效果 %+ @O#P
lnsrch,on !打开线性搜索 SN9kFFIPb=
outres,all,all !输出所有项 `YFtL
cnvtol,f,,0.005,2,0.5 !力收敛准则设定 UldXYtGe
cnvtol,u,,0.05,2,1 !位移收敛准则设定 ^W`<gR
!边坡在强度折减系数F=1时求解 2 Wt> Mi
allsel 5A)2} D]
solve !进行求解 zvYq@Mhr
save,F1,db !把F=1时求解结果保存 (Mo*^pVr
!边坡在强度折减系数F=1.2时求解 yh Yb'GK
finish K SbK EA
/solu -V-RP;">
allsel w j*,U~syB
asel,s,area,,2,3,1 !选择面积A2、A3 v~f HYa>
MPCHG,3,all ! 把所选择单元材料号改为3 04LI]'
allsel $:=A'd2
solve !进行求解 Pu7_
v
save,F1.2,db !把F=1.2时求解结果保存 B -?6M6#
!边坡在强度折减系数F=1.4时求解 F3N?Nk/
finish yCd-9zb=
/solu aGrIQq/k)%
allsel *rM^;4Zt
asel,s,area,,2,3,1 !选择面积A2、A3 9=v MgW
MPCHG,4,all !把所选择单元材料号改为4 oI0M%/aM
allsel ; kFDMuuO
solve !进行求解 [>+ 4^&
save,F1.4,db !把F=1.4时求解结果保存 *;l] 8.
!边坡在强度折减系数F=1.6时求解 (yu/l 6[
finish r54&XE]O
/solu N|8P)
allsel %824Cqdc
asel,s,area,,2,3,1 !选择面积A2、A3 S @jQX
MPCHG,5,all !把所选择单元材料号改为5 6*PYFf`
allsel K,Ef9c/+K
solve ;X8y Fq
save,F1.6,db !把F=1.6时求解结果保存 hEA<o67
!边坡在强度折减系数F=1.8时求解 EY^1Y3D w0
finish l7XUXbYp&=
/solu j#X.KM
allsel PXK7b2fE.
asel,s,area,,2,3,1 !选择面积A2、A3 s [M?as
MPCHG,6,all !把所选择单元材料号改为6 6_J$UBT
allsel i2@VB6]?
solve 6 CV*
Z\b
save,F1.8,db !把F=1.8时求解结果保存 y~B h
!边坡在强度折减系数F=2.0时求解 #ZJ _T`l
finish n &{Dq}q
/solu h%o%fH&F!
allsel W 3]_m8,Z
asel,s,area,,2,3,1 !选择面积A2、A3 xHUsFm s
MPCHG,7,all !把所选择单元材料号改为7 MuYk};f
allsel l Q'I
solve ;+e }aER&9
save,F2.0,db !把F=2.0时求解结果保存 Nh8Q b/::
!边坡在强度折减系数F=2.2时求解 j | v%)A
finish NTdi xfR
/solu v0
nj M
allsel (_niMQtF}
asel,s,area,,2,3,1 !选择面积A2、A3 Upc+Ukw
MPCHG,8,all !把所选择单元材料号改为8 Ee)T1 ~;W
allsel K$&s=Hm
solve >QjAoDVX ?
save,F2.2,db !把F=2.2时求解结果保存 $ yn];0$J
!边坡在强度折减系数F=2.4时求解 )bS~ 1n_0
finish ^`*9QjY
/solu X^T:8npxt
allsel Y'c>:;JEe
asel,s,area,,2,3,1 !选择面积A2、A3 (X $=Q6
MPCHG,9,all !把所选择单元材料号改为9 j|4<i9^}
allsel WzPTFw[
solve <L mIK
save,F2.4,db !把F=2.4时求解结果保存 nim*/LC[:
!边坡在强度折减系数F=2.6时求解 O}+.U<V
finish 3p3 9`"~
/solu C\S3Gs
allsel @KWb+?_H{<
asel,s,area,,2,3,1 !选择面积A2、A3 _K`wG}YIE
MPCHG,10,all !把所选择单元材料号改为10 v/Py "hQ
allsel 4iqoR$3Fc
solve 1{r3#MVL
save,F2.6,db !把F=2.6时求解结果保存 [=E
!边坡在强度折减系数F=2.8时求解 4H,`]B8(D
finish &R[ M c-2
/solu n(b(yXYm]
allsel ?+_Gs;DGVE
asel,s,area,,2,3,1 !选择面积A2、A3 4~k\ j
MPCHG,11,all !把所选择单元材料号改为11
txJr;
allsel GQ t8p[!
solve 8e*,jH 3
save,F2.8,db !把F=2.8时求解结果保存 gD,1 06%
!边坡在强度折减系数F=3.0时求解 xAqb\|$^
finish %b%-Ogz;4
/solu YNLV9.P6
allsel vL| SY_:4
asel,s,area,,2,3,1 !选择面积A2、A3 OglEt[ "
MPCHG,12,all !把所选择单元材料号改为12 OP~HdocB
allsel V^7V[(~`
solve ) T /0S$@
save,F3.0,db !把F=3.0时求解结果保存 bt"W(m&f
/post1 ！进入后处理 cQ/5qg
!边坡在强度折减系数F=1时结果分析 ~T ]m>A!
Resume,'F1','db' !读入边坡在强度折减系数F=1时 R{W E\T '
set,1,last !读入后一个子步 88VZR&v
pldisp,1 !绘制边坡模型变形图 tR}MrM
plnsol,u,x !绘制边坡模型水平方向位移云图 $ }<PL}+
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 I~q#eO)
!边坡在强度折减系数F=1.2时结果分析 E .r>7`E
Resume,'F1.2','db' !读入边坡在强度折减系数F=1.2时 1 9&<|qTz
set,1,last !读入后一个子步 /,89p&h
pldisp,1 !绘制边坡模型变形图 j.C`U(n}`
plnsol,u,x !绘制边坡模型水平方向位移云图 "w ^! /
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 %@wJ`F2a_
!边坡在强度折减系数F=1.4时结果分析 J,V9 k[88
Resume,'F1.4','db' !读入边坡在强度折减系数F=1.4时 pf%=h
|
set,1,last !读入后一个子步 h XfQ)$J
pldisp,1 !绘制边坡模型变形图 NgADKrDU
plnsol,u,x !绘制边坡模型水平方向位移云图 H(R1o~
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 $LKIT0
!边坡在强度折减系数F=1.6时结果分析 M=%l}FSTw(
Resume,'F1.6','db' !读入边坡在强度折减系数F=1.6时 o}$XH,-9&
set,1,last !读入后一个子步 a; rdQ>
pldisp,1 !绘制边坡模型变形图 5A%w 8Qv
plnsol,u,x !绘制边坡模型水平方向位移云图 @ >d*H75
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 b1^vd@(l x
!边坡在强度折减系数F=1.8时结果分析 >/:" D$
Resume,'F1.8','db' !读入边坡在强度折减系数F=1.8时 Ozw;(fDaU
set,1,last !读入后一个子步 JI? rL
pldisp,1 !绘制边坡模型变形图 y-T| #
plnsol,u,x !绘制边坡模型水平方向位移云图 *}]Nf
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 EqyeJq .
!边坡在强度折减系数F=2.0时结果分析 jq-p;-i
Resume,'F2.0','db' !读入边坡在强度折减系数F=2.0时 K-e9>fmB#
set,1,last !读入后一个子步 ZO!I .
pldisp,1 !绘制边坡模型变形图 sc|_Q/`\.
plnsol,u,x !绘制边坡模型水平方向位移云图 Qt iDTr
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 z8jk[5z
!边坡在强度折减系数F=2.2时结果分析 <A [E:*`*
Resume,'F2.2','db' !读入边坡在强度折减系数F=2.2时 `{eyvW[Ks
set,1,last !读入后一个子步 |
h`0u'#
pldisp,1 !绘制边坡模型变形图 &)8:h+&Z
plnsol,u,x !绘制边坡模型水平方向位移云图 2>|dF~"
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 8I-u2Y$Sr
!边坡在强度折减系数F=2.4时结果分析 e>7]w,*|
Resume,'F2.4','db' !读入边坡在强度折减系数F=2.4时 ]=.\-K
set,1,last !读入后一个子步 e&sim X;W
pldisp,1 !绘制边坡模型变形图 ?i)f^O
plnsol,u,x !绘制边坡模型水平方向位移云图 *v;!-F&8>
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 " W!M[qBW
!边坡在强度折减系数F=2.6时结果分析 c]$i\i#
Resume,'F2.6','db' !读入边坡在强度折减系数F=2.6时 Fw/6?:C}O6
set,1,last !读入后一个子步 Fgk a jig
pldisp,1 !绘制边坡模型变形图 k' gh
plnsol,u,x !绘制边坡模型水平方向位移云图 6FUw"|\u{
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 B\<Q ;RI2;
!边坡在强度折减系数F=2.8时结果分析 N96jJk
Resume,'F2.8','db' !读入边坡在强度折减系数F=2.8时 Ao&\E cIOT
set,1,last !读入后一个子步
*-+&[P]m
pldisp,1 !绘制边坡模型变形图 Y&[1` :-~-
plnsol,u,x !绘制边坡模型水平方向位移云图 R ? ,an2
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 ~res V
!边坡在强度折减系数F=3.0时结果分析 VGfMN|h
Resume,'F3.0','db' !读入边坡在强度折减系数F=3.0时 ?Eg (Gu.J
set,1,last !读入后一个子步 ?Y)vGlWDW<
pldisp,1 !绘制边坡模型变形图 Q~814P 8]
plnsol,u,x !绘制边坡模型水平方向位移云图 tkVbo.[8K
plnsol,eppl,eqv !绘制边坡模型塑性应变云图
