各位朋友： dn& s *
大家好！我通过ansys用有限元强度折减法分析了一个边坡的稳定性，结果不同的工况下得出的塑性区是一样的，我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激！附命令流和塑性区如下：(1)创建物理环境
{y)=eX9
/COM, Structural ! 指定结构分析 W8G,=d}6
/TITLE, Slope stability Analysis ! 定义工作标题 CT&|QH{
/FILNAM,support,1 ! 定义工作文件名 'lH|eU&-
(2)建立模型 b !+hH Hv:
!进入前处理器 U gr!"Q#M
/PREP7 ncaT?~u j
!定义单元类型 S;Fi? M
ET,1,PLANE82 !定义边坡围岩单元 atj (eg
KEYOPT,1,3,2 !求解类型为平面应变 {B~QQMEow
!定义材料属性 ?al'F q
!边坡围岩1材料属性 9=s<Ld
MP,EX,1,3E10 4VHn \
MP,PRXY,1,0.25 N:^n('U&j
MP,DENS,1,2500 !定义塑性模型 ><4<yj1
TB,DP,1 !定义数据表1 kXViWOXU^
TBDATA,1,9E5,42.8 !给数据表输入数据 jVEGj5F;N
!边坡围岩2材料属性 !采用弹性模型 YX!iL6 ?~
MP,EX,2,3.2E10 ！定义弹性模量 y#`tgJ:
MP,PRXY,2,0.24 !定义泊松比 Q-(zwAaE
MP,DENS,2,2700 !定义密度 q v-8)MSr
!取强度折减系数F=1.2时边坡围岩材料属性 ~]sc^[
MP,EX,3,3E10 t$`r4Lb9/
MP,PRXY,3,0.25 irZ])a
MP,DENS,3,2500 `~cqAs}6]Q
TB,DP,3 !定义数据表3 49eD1h3'X[
TBDATA,1,7.5E5,37.7 !给数据表输入数据 ez7A4>/
!取强度折减系数F=1.4时边坡围岩材料属性 |44Ploz2b
MP,EX,4,3E10 R8K& R\
MP,PRXY,4,0.25 !定义泊松比 |NlO7aQ>2H
MP,DENS,4,2500 %:i7s-0w
TB,DP,4 !定义数据表4 ~?l |
[
TBDATA,1,6.4E5,33.5 !给数据表输入数据 ;x y"\S]
!取强度折减系数F=1.6时边坡围岩材料属性 zOJ% }
MP,EX,5,3E10 ${DUCud,kY
MP,PRXY,5,0.25 )7hqJa-V
MP,DENS,5,2500 ！定义密度 QRw"H 8nW
TB,DP,5 !定义数据表5 L7l
FtX+b
TBDATA,1,5.6E5,30 !给数据表输入数据 VM ZMG$C
!取强度折减系数F=1.8时边坡围岩材料属性 kj Jn2c:y
MP,EX,6,3E10 ！定义弹性模量 n3WlZ!$
MP,PRXY,6,0.25 }H53~@WP>
MP,DENS,6,2500 aHD]k8 m z
TB,DP,6 !定义数据表6 9p]QM)M
TBDATA,1,5.0E5,27.2 ,Co|-DYf}
!取强度折减系数F=2.0时边坡围岩材料属性 HVRZ[Y<^
MP,EX,7,3E10 ldf\;Qk
MP,PRXY,7,0.25 !定义泊松比 U svl}{L[
MP,DENS,7,2500 [DuttFX^x
TB,DP,7 !定义数据表7 d z|or9&
TBDATA,1,4.5E5,24.8 P1!qbFDv8
!取强度折减系数F=2.2时边坡围岩材料属性 -uS!\
MP,EX,8,3E10 )705V|v
MP,PRXY,8,0.25 &bS ,hbDt
MP,DENS,8,2500 Zj( AJ*r
TB,DP,8 !定义数据表8 <|HV. O/!
TBDATA,1,4.09E5,22.8 X;$+,&M"
!取强度折减系数F=2.4时边坡围岩材料属性 `Gs9Xmc|
MP,EX,9,3E10 _T60;ZI+^
MP,PRXY,9,0.25 ?4Y GT
MP,DENS,9,2500 ' B |JAi?
TB,DP,9 !定义数据表9 a,,exi
TBDATA,1,3.6E5,21.1 !给数据表输入数据 ?d*z8 w
!取强度折减系数F=2.6时边坡围岩材料属性 yNPVOp*
MP,EX,10,3E10 ！定义弹性模量 @@f"%2ZR[
MP,PRXY,10,0.25 !定义泊松比 -abt:or
MP,DENS,10,2500 ！定义密度 0S"mVZ*P
TB,DP,10 !定义数据表10 "69s) ~
TBDATA,1,3.46E5,19.6 !给数据表10输入数据 8|gIhpO?^
!取强度折减系数F=2.8时边坡围岩材料属性 t5Sy V:fP
MP,EX,11,3E10 ！定义弹性模量 [+Iz@0 q
MP,PRXY,11,0.25 !定义泊松比 /'SNw?&
MP,DENS,11,2500 Zpt\p7WQ
TB,DP,11 !定义数据表11 R*, MfV
TBDATA,1,3.2E5,18.2 !给数据表11输入数据 Cp\6W[2+B
!取强度折减系数F=3.0时边坡围岩材料属性 @NR>{Eg
MP,EX,12,3E10 ！定义弹性模量 poE0{HOU
MP,PRXY,12,0.25 {qk1_yP
MP,DENS,12,2500 ！定义密度 ~g91Pr
TB,DP,12 !定义数据表12 sJKI!
TBDATA,1,3.0E5,17.1 !给数据表12输入数据 #<fRE"v:Q
!建立几何模型 _;"il%l=1
!创建边坡线模型线 p%ki>p )E|
K,1,,,, ！ 创建关键点1 #m xP w
K,2,-800,,, &$+ AXzn
K,3,-800,-800,, q])K,)
K,4,-800,-1200,, ,~U>'&M;
K,5,1200,-1200,, N%@Qf~
K,6,1200,-800,, !|(-=2`
K,7,1200,,, -OV&Md:~
K,8,1200,400,, 1er
TldX
K,9,400,400,, g b1V ~
l,1,2 !连接关键点1、2生成直线L1 6jaEv#
l,2,3 2Ah#<k-gC;
l,3,4 /|}EL%a
l,4,5 {p2!|A&a
l,5,6 iqsCB%;5
l,6,7 l$KA)xbI
l,7,8 cVv=*81 \
l,8,9 t 9lPb_70
l,9,1 v&\Q8!r_
l,1,7 FaAC&F@u
l,3,6 w7L{_aom
!创建边坡面模型 MpT8" /.]A
al,3,4,5,11 !由线L3、L4、L5、L11生成一个面A1 \
#F
al,1,2,11,6,10 !由线L1、L2、L1、L6、L10生成一个面A2 70d1ReQ
al,7,8,9,10 !由线L7、L8、L9、L10生成一个面A3 +Ze} B*0
/PNUM, area, on !打开面号开关 [g |_~h
!划分网格生成有限元模型 hPkp;a #
! 划分边坡围岩2网格 :
$1?i)
!设置网格份数 =IZT(8
lsel,s,,,3,5,2 !选择线L3和L5 qZd QD
lesize,all,,,5 !把所选择线分为5段 ,)cM3 nu
lsel,s,,,4,11,7 !选择线L4和L11 M/f<A$xx_
lesize,all,,,26 !把所选择线分为26段 @?s Rj&w
mat,2 !给边坡围岩2赋予2号材料特性 s S+MqBh&I
type,1 !采用单元类型1 AYBns]!
mshkey,0 !设定自由网格划分 'ms-*c &
mshape,0 !设定四边形网格划分 #^0R&) T
amesh,1 ！划分面积A1 = jN.1}
! 划分边坡围岩1网格 VD*6g%p
!设置网格份数 b=C*W,Q_#
lsel,s,,,2,6,4 !选择线L2和L6 C O/]wS
lsel,a,,,1 !附加选择线L1 zpn9,,~u
lesize,all,,,10 !把所选择线分为10段 `v!urE/gg%
lsel,s,,,8,10,2 !选择线L8和L10 , >a&"V^k
lesize,all,,,16 !把所选择线分为16段 ^d Wa;m]l
lsel,s,,,7,9,2 WCZjXDiwJ
lesize,all,,,12 jVe1b1rt~3
mat,1 !给边坡围岩1赋予1号材料特性 :U|1xgB
type,1 iQ67l\{R
mshkey,0 !设定自由网格划分 B`)BZ,#p
mshape,0 !设定四边形网格划分 )MVz$h{c.]
amesh,2 !划分面积A2 >58YjLXb
amesh,3 ！划分面积A3 Pm6p v;WK
allsel [>I<#_^~
!保存网格模型 j^sg6.Z*
save,Slope-grid.db +fB5w?Rg
(3)施加约束和荷载 (XTG8W sN
!给边坡模型施加约束 J3V=
46Yc
!边坡两侧施加X方向约束 Oi.C(@^(
nsel,s,loc,x,-800 !选择X=-800线上所有节点 uo9B 9"&
nsel,a,loc,x,1200 !选择X=1200线上所有节点 tAd%#: K
d,all,ux !对所选择节点约束X方向位移 /xBb[44z8
allsel I9A~Ye
5O&
!边坡底部施加约束
z _$%-6
nsel,s,loc,y,-1200 !选择Y=--1200线上所有节点 P8:dU(nlW
d,all,uy !对所选择节点约束X、Y方向位移 Y(y kng
d,all,ux 5vQHhwO50k
!施加重力加速度 5Pc;5
o0C
acel,,9.8 s[>,X#7 y
(4)求解 au(D66VO
/solu mthA4sz
!求解设置 r8?gD&c}
antype,static !设定为静力求解 n&4N[Qlv,
nsubst,100 !设定最大子步数为100 8
/]S^'>
pred,on !打开时间步长预测器 C} j "Qi`
nropt,full !设定牛顿-拉普森选项 ?FZ HrA
nlgeom,on !打开大位移效果 N{!i =A
lnsrch,on !打开线性搜索 l 'rja.\
outres,all,all !输出所有项 {lzWrUGO
cnvtol,f,,0.005,2,0.5 !力收敛准则设定 P= BZ+6DS
cnvtol,u,,0.05,2,1 !位移收敛准则设定 QW~E&B%
!边坡在强度折减系数F=1时求解 o'aEY<mZ7
allsel =ZznFVJ`={
solve !进行求解 QE+g
j 8
save,F1,db !把F=1时求解结果保存 dES"@?!^
!边坡在强度折减系数F=1.2时求解 1ba~ SHi
finish &J]K3w1p
/solu :`#d:.@]o@
allsel P bn*_/H
asel,s,area,,2,3,1 !选择面积A2、A3 QO:!p5 ^:
MPCHG,3,all ! 把所选择单元材料号改为3 x;.Jw 6g
allsel )u&|_&g{}J
solve !进行求解 9.M4o[
save,F1.2,db !把F=1.2时求解结果保存 L<cx:Vz
!边坡在强度折减系数F=1.4时求解 )
w5SUb
finish k9R4Y\8P
/solu *8 A
allsel NN{?z!
asel,s,area,,2,3,1 !选择面积A2、A3 C3f' {}
MPCHG,4,all !把所选择单元材料号改为4 }0*@fO
allsel ! I:% 0D
solve !进行求解 L[fiU 0^o
save,F1.4,db !把F=1.4时求解结果保存 Tk[ $5u*,
!边坡在强度折减系数F=1.6时求解 9<?M 8_
finish p$c6<'UqH
/solu W+c<2? d:
allsel e)k9dOR
asel,s,area,,2,3,1 !选择面积A2、A3 x j)F55e?
MPCHG,5,all !把所选择单元材料号改为5 [85spub&}
allsel HyQJXw?A:
solve ( $MlXBI
save,F1.6,db !把F=1.6时求解结果保存 u@) U "FZ
!边坡在强度折减系数F=1.8时求解 D /gw .XYL
finish x8B}ZIbT9
/solu .hb:s,0mP
allsel
Mx?d
asel,s,area,,2,3,1 !选择面积A2、A3 3pROf#M
MPCHG,6,all !把所选择单元材料号改为6 hh%-(HaLX3
allsel n38p!oS
solve B"w?;EeV.
save,F1.8,db !把F=1.8时求解结果保存 ub0.J#j@
!边坡在强度折减系数F=2.0时求解 a5^] 20Fa
finish Z clQ
/solu jRlYU`?
allsel <$$yw=ef
asel,s,area,,2,3,1 !选择面积A2、A3 7aRi 5
MPCHG,7,all !把所选择单元材料号改为7 ` $IK`O
allsel p`dU2gV
solve fplow
save,F2.0,db !把F=2.0时求解结果保存 2a)xTA#
!边坡在强度折减系数F=2.2时求解 FX&~\kmV'j
finish Lg+Ac5 y}`
/solu |^I0dR/w:
allsel +)om^e@.
asel,s,area,,2,3,1 !选择面积A2、A3
_"yh.N&
MPCHG,8,all !把所选择单元材料号改为8 qA7>vi%
allsel pU}(@oy
solve k"%~" 9
save,F2.2,db !把F=2.2时求解结果保存 &t@jl\ND
!边坡在强度折减系数F=2.4时求解 2zA4vZkbcw
finish S3%FHS
/solu s c,Hq\$&
allsel -);Wfs
asel,s,area,,2,3,1 !选择面积A2、A3 4Z=_,#h4.
MPCHG,9,all !把所选择单元材料号改为9 \:'/'^=#|
allsel tS5hv@9cWx
solve Rok7n1gW
save,F2.4,db !把F=2.4时求解结果保存
DPxM'7
!边坡在强度折减系数F=2.6时求解 Ug SB>V<?
finish r,3DT Be
/solu O6 3<AY@
allsel ?3,:-"(@p
asel,s,area,,2,3,1 !选择面积A2、A3 2wg5#i
MPCHG,10,all !把所选择单元材料号改为10 .VJMz4$]O
allsel )EuvRLo{S7
solve CsR$c,8X.
save,F2.6,db !把F=2.6时求解结果保存 uAq~=)F>,
!边坡在强度折减系数F=2.8时求解 Kk0g0C:"EO
finish {]!mrAjD
/solu &{hL&BLr
allsel i# / Jr=
asel,s,area,,2,3,1 !选择面积A2、A3 YlQ=5u^+
MPCHG,11,all !把所选择单元材料号改为11 { lDd.Fn
allsel d"mkL -
solve ah &D%8E
save,F2.8,db !把F=2.8时求解结果保存 [b%D3-}'
!边坡在强度折减系数F=3.0时求解 /Iy]DU8
finish >8^
$ [}w
/solu SM#]H-3
allsel X7 MM2V
asel,s,area,,2,3,1 !选择面积A2、A3 !Pvf;rNI1T
MPCHG,12,all !把所选择单元材料号改为12 n)-$e4u2
allsel gfd"v
solve {6|G@ ""O
save,F3.0,db !把F=3.0时求解结果保存 g)[V(yWu
(5）后处理 %XDc,AR[
/post1 ！进入后处理 LmrfN?5
!边坡在强度折减系数F=1时结果分析 HZB>{O
Resume,'F1','db' !读入边坡在强度折减系数F=1时 myQagqRx
set,1,last !读入后一个子步 P ) "m0Lu<
pldisp,1 !绘制边坡模型变形图 5lmHotj#
plnsol,u,x !绘制边坡模型水平方向位移云图 2;`1h[, -^
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 kCF>nt@
!边坡在强度折减系数F=1.2时结果分析 #Y `~(K47
Resume,'F1.2','db' !读入边坡在强度折减系数F=1.2时 ZF8 yw(z
set,1,last !读入后一个子步 ?
(Oy\
pldisp,1 !绘制边坡模型变形图 7IH@oMvE
plnsol,u,x !绘制边坡模型水平方向位移云图 %N6A+ 5H
plnsol,eppl,eqv !绘制边坡模型塑性应变云图
(N6 i4
g6
!边坡在强度折减系数F=1.4时结果分析 2#]#sZmk
Resume,'F1.4','db' !读入边坡在强度折减系数F=1.4时 z,
) 6"/;
set,1,last !读入后一个子步 xh,qNnGG i
pldisp,1 !绘制边坡模型变形图 l/ GGC nO/
plnsol,u,x !绘制边坡模型水平方向位移云图 ^zmG0EH,
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 6vo;! V6
!边坡在强度折减系数F=1.6时结果分析 k,6f
&#x
Resume,'F1.6','db' !读入边坡在强度折减系数F=1.6时 }OR@~V{Gj
set,1,last !读入后一个子步 /4V#C-
pldisp,1 !绘制边坡模型变形图 G6P?2@
plnsol,u,x !绘制边坡模型水平方向位移云图 N^G
Mp,8
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 xC:L)7#aw
!边坡在强度折减系数F=1.8时结果分析 IqHV)A
Resume,'F1.8','db' !读入边坡在强度折减系数F=1.8时 qJs<#MQ2
set,1,last !读入后一个子步 A RuA<vQ
pldisp,1 !绘制边坡模型变形图 #U4F0BdA
plnsol,u,x !绘制边坡模型水平方向位移云图
Y_IF;V\
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 Gr'
CtO
!边坡在强度折减系数F=2.0时结果分析 sqwGsO$#
Resume,'F2.0','db' !读入边坡在强度折减系数F=2.0时 bHYy}weZ
set,1,last !读入后一个子步 BO ;tCEV?
pldisp,1 !绘制边坡模型变形图 X/!o\yyT
plnsol,u,x !绘制边坡模型水平方向位移云图 D,*3 w'X!K
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 6 7.+
.2
!边坡在强度折减系数F=2.2时结果分析 orpriO|qD
Resume,'F2.2','db' !读入边坡在强度折减系数F=2.2时 [Td4K.c
set,1,last !读入后一个子步 -HbC!w v
pldisp,1 !绘制边坡模型变形图 {X+3;&@
plnsol,u,x !绘制边坡模型水平方向位移云图 bdrg(d6
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 O,
wJR
!边坡在强度折减系数F=2.4时结果分析 S~bOUdV
Z
Resume,'F2.4','db' !读入边坡在强度折减系数F=2.4时 K(rWNO
set,1,last !读入后一个子步 .t-4o<7 3
pldisp,1 !绘制边坡模型变形图 [wOn|)&
&
plnsol,u,x !绘制边坡模型水平方向位移云图 TDKki(o=~
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 z+wA
rPxc
!边坡在强度折减系数F=2.6时结果分析 BLdvyVFx
Resume,'F2.6','db' !读入边坡在强度折减系数F=2.6时 !u[9a;Sa#
set,1,last !读入后一个子步 ItVWO:x&v
pldisp,1 !绘制边坡模型变形图 } 5[qo`M
plnsol,u,x !绘制边坡模型水平方向位移云图 %6,SKg p
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 ".V$~n(
!边坡在强度折减系数F=2.8时结果分析 +F` S>U
Resume,'F2.8','db' !读入边坡在强度折减系数F=2.8时 k68T`Ub\W6
set,1,last !读入后一个子步 # e1>H1eU
pldisp,1 !绘制边坡模型变形图 'Cfl*iNb
plnsol,u,x !绘制边坡模型水平方向位移云图 z&)A,ryW0
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 faX #**r
!边坡在强度折减系数F=3.0时结果分析 . B9iLI
Resume,'F3.0','db' !读入边坡在强度折减系数F=3.0时 29" 'K.r
set,1,last !读入后一个子步 LVfF[
pldisp,1 !绘制边坡模型变形图 qp }Cqi
plnsol,u,x !绘制边坡模型水平方向位移云图 DB|Y
plnsol,eppl,eqv !绘制边坡模型塑性应变云图
