(1)创建物理环境 qca,a3k
/COM, Structural nM#/uuRl|
/TITLE, Slope stability Analysis lrQNl^K }=
/FILNAM,support,1 `X)y5*##wq
(2)建立模型 2 z E gAc
!进入前处理器 @:s (L]
/PREP7 %JoHc?
!定义单元类型 tx`gXtO$
ET,1,PLANE82 !定义边坡围岩单元 5! -+5TJI
KEYOPT,1,3,2 !求解类型为平面应变 BRSI g]
!定义材料属性 Z P-^10
!边坡围岩1材料属性 X/ Sp!W-H
MP,EX,1,3E10 >L4q>S^v
MP,PRXY,1,0.25 [L(qrAQ2|z
MP,DENS,1,2500 !定义塑性模型 ,!
b9
TB,DP,1 !定义数据表1 bbNN$-S|
TBDATA,1,9E5,42.8 !给数据表输入数据 #w]UP#^io
!边坡围岩2材料属性 !采用弹性模型 1z IX
$A
MP,EX,2,3.2E10 ！定义弹性模量 ~xPU#m<
MP,PRXY,2,0.24 !定义泊松比 U</Vcz
MP,DENS,2,2700 !定义密度 c6-~PKJL
!取强度折减系数F=1.2时边坡围岩材料属性 `-Y8T\
MP,EX,3,3E10 9 n0 ?0mk
MP,PRXY,3,0.25 fj"1TtPq#
MP,DENS,3,2500 ? $$Xg3w_#
TB,DP,3 !定义数据表3 V) xwl vX
TBDATA,1,7.5E5,37.7 !给数据表输入数据 9 4.|l
!取强度折减系数F=1.4时边坡围岩材料属性 U-+o6XX
MP,EX,4,3E10 Y(mnGaV n
MP,PRXY,4,0.25 !定义泊松比 yr vV<}
MP,DENS,4,2500 x_L5NsO:
TB,DP,4 !定义数据表4 AcHr X=O
TBDATA,1,6.4E5,33.5 !给数据表输入数据 t9l]ie{"o.
!取强度折减系数F=1.6时边坡围岩材料属性 aoqG*qh}b
MP,EX,5,3E10 $Iz *W]B!
MP,PRXY,5,0.25 +-2 W{lX
MP,DENS,5,2500 ！定义密度 T"IDCT'z
TB,DP,5 !定义数据表5 '< =77yDg
TBDATA,1,5.6E5,30 !给数据表输入数据 !1m7^3l7j
!取强度折减系数F=1.8时边坡围岩材料属性 z\7-v<ZS
MP,EX,6,3E10 ！定义弹性模量 {gT4Oq__
MP,PRXY,6,0.25 D*0[7:NSO
MP,DENS,6,2500 BcXPgM!Xqz
TB,DP,6 !定义数据表6 Dm{9;Abs%
TBDATA,1,5.0E5,27.2 pgUp1goAU
!取强度折减系数F=2.0时边坡围岩材料属性 p ; ]Qxh
MP,EX,7,3E10 8f`r!/j
MP,PRXY,7,0.25 !定义泊松比 >uLWfk+y1
MP,DENS,7,2500 ?COLjk
TB,DP,7 !定义数据表7 H^ds<I<)
TBDATA,1,4.5E5,24.8 nz2`YyR
!取强度折减系数F=2.2时边坡围岩材料属性 @V}!elV
MP,EX,8,3E10 W79Sz}):
MP,PRXY,8,0.25 CWdpF>En
MP,DENS,8,2500 F H byL\Q
TB,DP,8 !定义数据表8 #M ;j*IBl*
TBDATA,1,4.09E5,22.8 t4d^DZDh!
!取强度折减系数F=2.4时边坡围岩材料属性 >bR oQ 8
MP,EX,9,3E10 Z X0 #I W
MP,PRXY,9,0.25 `_"loPu
MP,DENS,9,2500 0q6xX NAX
TB,DP,9 !定义数据表9 .RW&=1D6
TBDATA,1,3.6E5,21.1 !给数据表输入数据 CXiDe)|<E
!取强度折减系数F=2.6时边坡围岩材料属性 z"%{SI^
MP,EX,10,3E10 ！定义弹性模量 '><I|c}
MP,PRXY,10,0.25 !定义泊松比 \JLea$TM:
MP,DENS,10,2500 ！定义密度 N++ ;}j
TB,DP,10 !定义数据表10 _~Id~b
TBDATA,1,3.46E5,19.6 !给数据表10输入数据 E %%iVFPX
!取强度折减系数F=2.8时边坡围岩材料属性 GHWt3K:*w
MP,EX,11,3E10 ！定义弹性模量 utzf7?nIS
MP,PRXY,11,0.25 !定义泊松比 @ b&_xT
MP,DENS,11,2500 .~ O-
<P#
TB,DP,11 !定义数据表11 um,G^R
TBDATA,1,3.2E5,18.2 !给数据表11输入数据 A'6-E{
!取强度折减系数F=3.0时边坡围岩材料属性 u(W>HVEG
MP,EX,12,3E10 ！定义弹性模量 "UY lC0 S\
MP,PRXY,12,0.25 vC^Ul
MP,DENS,12,2500 ！定义密度 KOhK#t>H@0
TB,DP,12 !定义数据表12 QtHK`f>4#n
TBDATA,1,3.0E5,17.1 !给数据表12输入数据 awB+B8^ s
!建立几何模型 K0A[xkX6
!创建边坡线模型线 U%rEW[ j
K,1,,,, ！ 创建关键点1 u~8=ik n+T
K,2,-800,,, j82x$I*
K,3,-800,-800,, %p;;aZG
K,4,-800,-1200,, `a6AE S'w$
K,5,1200,-1200,, "W^+NeLc
K,6,1200,-800,, :P8X?C63W]
K,7,1200,,, gT_tR_g
K,8,1200,378,, 4|/}~9/
K,9,430,378,, h~pQ
l,1,2 !连接关键点1、2生成直线L1 8hV>Q
l,2,3 6c6 w w"
l,3,4 xp*Wf#BF
l,4,5 )k o[_OJj
l,5,6 A1Es>NK[qW
l,6,7 Bv xLbl}
l,7,8 XOL_vS24
l,8,9 =Jax T90x
l,9,1 S uo%uD
l,1,7 FJD;LpW
l,3,6 PiIP%$72O
!创建边坡面模型 'ws@I?! r
al,3,4,5,11 !由线L3、L4、L5、L11生成一个面A1 NG5k9pJ
al,1,2,11,6,10 !由线L1、L2、L1、L6、L10生成一个面A2 H# H[8#
al,7,8,9,10 !由线L7、L8、L9、L10生成一个面A3 s|vx2-Cu]
/PNUM, area, on !打开面号开关 2WUl8?f2Y
!划分网格生成有限元模型 Egt !N
! 划分边坡围岩2网格 1<G, 0Lt
!设置网格份数 sa?s[
lsel,s,,,3,5,2 !选择线L3和L5 ) vD:
lesize,all,,,5 !把所选择线分为5段 .^xQtnq
lsel,s,,,4,11,7 !选择线L4和L11 XwZR
Kh\>=
lesize,all,,,26 !把所选择线分为26段 {ui{Y c
mat,2 !给边坡围岩2赋予2号材料特性 ,K15KN.'
type,1 !采用单元类型1 bn:74,GeyK
mshkey,0 !设定自由网格划分 " M3S
mshape,0 !设定四边形网格划分 V}h
<,E9
amesh,1 ！划分面积A1 A'aY H`j
! 划分边坡围岩1网格 5fq 4[a
!设置网格份数 O03N$ Jq
A
lsel,s,,,2,6,4 !选择线L2和L6 (M# m BS
lsel,a,,,1 !附加选择线L1 X=@bzL;eq
lesize,all,,,10 !把所选择线分为10段 P"{yV?CNg
lsel,s,,,8,10,2 !选择线L8和L10 NOSL b];
lesize,all,,,16 !把所选择线分为16段 v%muno,
lsel,s,,,7,9,2 Hb3..o:
lesize,all,,,12 .4J7 ^l
mat,1 !给边坡围岩1赋予1号材料特性 }ijFvIHV
type,1 9fy[%M
mshkey,0 !设定自由网格划分 rL,kDSLs
mshape,0 !设定四边形网格划分 7Y.mp9,
amesh,2 !划分面积A2 )mH(Hx
amesh,3 ！划分面积A3 9\Md .>
allsel S&F[\4w5]
!保存网格模型 1\aV4T
save,Slope-grid.db Df@b;-E
(3)施加约束和荷载 K BlJJH`z{
!给边坡模型施加约束 G){A&F
!边坡两侧施加X方向约束 *+TO% {4
nsel,s,loc,x,-800 !选择X=-800线上所有节点 7q&T2?GEN
nsel,a,loc,x,1200 !选择X=1200线上所有节点 h$]nfHi_Q
d,all,ux !对所选择节点约束X方向位移 )i"52!
allsel zeshM8=
!边坡底部施加约束 G:!3X) b
nsel,s,loc,y,-1200 !选择Y=--1200线上所有节点 5cj&D74o
d,all,uy !对所选择节点约束X、Y方向位移 # AH gY.
d,all,ux KFBBqP
!施加重力加速度 1* ?XI
acel,,9.8 *X!+wK-+
(4)求解 ~^/BAc
/solu r ?<?0j
!求解设置 KBDNK_7A
antype,static !设定为静力求解 fQxlYD'peb
nsubst,100 !设定最大子步数为100 ]+5Y\~I
pred,on !打开时间步长预测器 Z|B`n
SzH
nropt,full !设定牛顿-拉普森选项 l0PXU)>C
nlgeom,on !打开大位移效果 cwE ?+vB
lnsrch,on !打开线性搜索 *P\OP'o_
outres,all,all !输出所有项 u
p zBd]
cnvtol,f,,0.005,2,0.5 !力收敛准则设定 =4uO"o
cnvtol,u,,0.05,2,1 !位移收敛准则设定 V]Kk =
!边坡在强度折减系数F=1时求解 *+% $OH,
allsel 0DaKd<Scv
solve !进行求解 ^|%N _ s
save,F1,db !把F=1时求解结果保存 0
s@>e
!边坡在强度折减系数F=1.2时求解 XMF#l]P
finish DfjDw/{U3L
/solu pS "A{k)i
allsel s5 4AM]a{j
asel,s,area,,2,3,1 !选择面积A2、A3 *SYuq)
MPCHG,3,all ! 把所选择单元材料号改为3 +h? Gps
allsel -amNz.`[PR
solve !进行求解 ]u.)6{
save,F1.2,db !把F=1.2时求解结果保存 *JOp)e0 b
!边坡在强度折减系数F=1.4时求解 JM fv|>=
finish FI @kE19
/solu 0HI0/Tvu$<
allsel qzWn l[3
asel,s,area,,2,3,1 !选择面积A2、A3 W[LQ $uj
MPCHG,4,all !把所选择单元材料号改为4 +^q- v-
allsel FwV 5{-(
solve !进行求解 'soll[J
save,F1.4,db !把F=1.4时求解结果保存 I@kMM12>c
!边坡在强度折减系数F=1.6时求解
c`}YL4
finish 8iPA^b|sz{
/solu J ql$
g
allsel S;" $02]
asel,s,area,,2,3,1 !选择面积A2、A3 4}t$Lf_
MPCHG,5,all !把所选择单元材料号改为5 J;k8 a2$_
allsel qFp]jbU
solve E J&w6),d
save,F1.6,db !把F=1.6时求解结果保存 GPrq(
!边坡在强度折减系数F=1.8时求解 Tdxc%'l
finish a+B3`6
/solu )`#SMLMy~
allsel :j [=
asel,s,area,,2,3,1 !选择面积A2、A3 PKZMuEEy,
MPCHG,6,all !把所选择单元材料号改为6 ~
Ofn&[G
allsel f3*SIKi
solve nTE\EZ+=2
save,F1.8,db !把F=1.8时求解结果保存 8CUl |I ~
!边坡在强度折减系数F=2.0时求解 WM0-F@_
finish MSb0J `
/solu D1V^Db Um_
allsel WeIi{<u8R
asel,s,area,,2,3,1 !选择面积A2、A3 lv G+9e3+
MPCHG,7,all !把所选择单元材料号改为7 H on,-<
allsel To;r#h
solve 1QbD ]"=n
save,F2.0,db !把F=2.0时求解结果保存 7g4IAs oD
!边坡在强度折减系数F=2.2时求解 r.' cjUs
finish ?NxaJ^
/solu o, qUf
allsel Xc9NM1bp=
asel,s,area,,2,3,1 !选择面积A2、A3 K8uqLSP '
MPCHG,8,all !把所选择单元材料号改为8 &<F w
allsel 6RfS _
solve Ny$N5/b!!
save,F2.2,db !把F=2.2时求解结果保存 ),@m
3wQ
!边坡在强度折减系数F=2.4时求解 geNvp0
finish 6 u,w
/solu c3 )6{
allsel 8aY}b($*ZI
asel,s,area,,2,3,1 !选择面积A2、A3 f9HoQDFsM
MPCHG,9,all !把所选择单元材料号改为9 @ =? #nB&
allsel n{!=g R.v.
solve 7WHq'R{@
save,F2.4,db !把F=2.4时求解结果保存 >&JS-j Fg
!边坡在强度折减系数F=2.6时求解 !]MGIh#u
finish ^V"0 8
/solu &S[>*+}{+
allsel 2E.D0E Cu
asel,s,area,,2,3,1 !选择面积A2、A3 t;
@T~%
MPCHG,10,all !把所选择单元材料号改为10 z>HM$n`YD
allsel Dc3bG@K*G
solve BO>[\!=y
save,F2.6,db !把F=2.6时求解结果保存 @Ll^z e&HI
!边坡在强度折减系数F=2.8时求解 v 807)JwS
finish .X"\ Mg
/solu dF^`6-K1
allsel ^@$T> SB1
asel,s,area,,2,3,1 !选择面积A2、A3 g{Hb3id9
MPCHG,11,all !把所选择单元材料号改为11 |H%,>r`9S
allsel ahR-^^'$
solve o;v_vCLO
save,F2.8,db !把F=2.8时求解结果保存 p[%B#(]9,
!边坡在强度折减系数F=3.0时求解 -+Z&O?pSH
finish op,mP0b
/solu loD:4e1
allsel #;\tgUQ
asel,s,area,,2,3,1 !选择面积A2、A3 X
CHN'l'
MPCHG,12,all !把所选择单元材料号改为12 in> ?kbaG+
allsel t?FPmbj v
solve N p? / r}
save,F3.0,db !把F=3.0时求解结果保存 &g!yRvM!;Q
/post1 ！进入后处理 #Wt1Ph_;
!边坡在强度折减系数F=1时结果分析 *X 2dS
{
Resume,'F1','db' !读入边坡在强度折减系数F=1时 ~"cqFdnO
set,1,last !读入后一个子步 R aA7 U
pldisp,1 !绘制边坡模型变形图 v {) *P.E
plnsol,u,x !绘制边坡模型水平方向位移云图 H284
]i
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 <%"CQT6g %
!边坡在强度折减系数F=1.2时结果分析 -l < [CI
Resume,'F1.2','db' !读入边坡在强度折减系数F=1.2时 8* sP
set,1,last !读入后一个子步 FXbalQ?^
pldisp,1 !绘制边坡模型变形图 Sr-!-eC
plnsol,u,x !绘制边坡模型水平方向位移云图 }q)dXFL=I#
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 T9AFL;1
!边坡在强度折减系数F=1.4时结果分析 r#c+{yY
Resume,'F1.4','db' !读入边坡在强度折减系数F=1.4时 Sj;:*jk!h
set,1,last !读入后一个子步 V0*9Tnc
pldisp,1 !绘制边坡模型变形图 qSQsY:]j0
plnsol,u,x !绘制边坡模型水平方向位移云图 /< \do 1
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 {'o\#4 Wk
!边坡在强度折减系数F=1.6时结果分析 oYX#VX
Resume,'F1.6','db' !读入边坡在强度折减系数F=1.6时 H6Kt^s<6xu
set,1,last !读入后一个子步 mW#p& {
pldisp,1 !绘制边坡模型变形图 Cp]q>lM"
plnsol,u,x !绘制边坡模型水平方向位移云图 \O\veB8
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 G C@U['
!边坡在强度折减系数F=1.8时结果分析 R}$A>)%dx
Resume,'F1.8','db' !读入边坡在强度折减系数F=1.8时 K>Tv M &
set,1,last !读入后一个子步 ~g &Gi)je
pldisp,1 !绘制边坡模型变形图
?dvcmXR
plnsol,u,x !绘制边坡模型水平方向位移云图 A[Vhy;xz
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 S^)xioKsJ
!边坡在强度折减系数F=2.0时结果分析 ,Z(J; ~
Resume,'F2.0','db' !读入边坡在强度折减系数F=2.0时 \; zix(N[5
set,1,last !读入后一个子步 4x$Ts %]
pldisp,1 !绘制边坡模型变形图 t7/a 5x
plnsol,u,x !绘制边坡模型水平方向位移云图 \7q>4[
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 ~t^'4"K*
!边坡在强度折减系数F=2.2时结果分析 ;h4w<OqcM
Resume,'F2.2','db' !读入边坡在强度折减系数F=2.2时 y<)q;fI7
set,1,last !读入后一个子步 | E FbT>
pldisp,1 !绘制边坡模型变形图 \P9HAz'6
plnsol,u,x !绘制边坡模型水平方向位移云图 9U!#Y%*T
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 $kh6-y@
!边坡在强度折减系数F=2.4时结果分析 +?Y(6$o
Resume,'F2.4','db' !读入边坡在强度折减系数F=2.4时 << aAYkx <
set,1,last !读入后一个子步 eSSv8 [u
pldisp,1 !绘制边坡模型变形图 { pu .l4nk
plnsol,u,x !绘制边坡模型水平方向位移云图 0*:4@go0}i
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 {&h=
!边坡在强度折减系数F=2.6时结果分析 XtIY8wsP
Resume,'F2.6','db' !读入边坡在强度折减系数F=2.6时 @qB1:==@7
set,1,last !读入后一个子步 Tl*FK?)MC^
pldisp,1 !绘制边坡模型变形图 gal.<SVW
plnsol,u,x !绘制边坡模型水平方向位移云图 ;CA7\&L>
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 ETaLE[T%1
!边坡在强度折减系数F=2.8时结果分析 jxZd
=%7Q
Resume,'F2.8','db' !读入边坡在强度折减系数F=2.8时 ~ym-Szo
set,1,last !读入后一个子步 }#E~XlX ^
pldisp,1 !绘制边坡模型变形图 ys9MV%*
plnsol,u,x !绘制边坡模型水平方向位移云图 %loe8yt
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 Es+BV+x[.c
!边坡在强度折减系数F=3.0时结果分析 [4HOWM>\
Resume,'F3.0','db' !读入边坡在强度折减系数F=3.0时 M!iYj+nrP
set,1,last !读入后一个子步 ANd#m9(x
pldisp,1 !绘制边坡模型变形图 T] ?QCf
plnsol,u,x !绘制边坡模型水平方向位移云图 _UUp+Hz
plnsol,eppl,eqv !绘制边坡模型塑性应变云图
