各位朋友： /`)>W :
大家好！我通过ansys用有限元强度折减法分析了一个边坡的稳定性，结果不同的工况下得出的塑性区是一样的，我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激！附命令流和塑性区如下：(1)创建物理环境 !~Am1\02
/COM, Structural ! 指定结构分析 4Z"JC9As
/TITLE, Slope stability Analysis ! 定义工作标题 3$E\B=7/U
/FILNAM,support,1 ! 定义工作文件名 XX@@tzN
(2)建立模型 p~h)@
!进入前处理器 afJ`1l
/PREP7 iCK p"(kf
!定义单元类型 GNIZHyT(O
ET,1,PLANE82 !定义边坡围岩单元 TGe)%jZ
KEYOPT,1,3,2 !求解类型为平面应变 l[)ZEEP
!定义材料属性 '=^$ ;3Z
!边坡围岩1材料属性 <&:OSd:%
MP,EX,1,3E10 lB9 9J"A
MP,PRXY,1,0.25 zi l^^wT0J
MP,DENS,1,2500 !定义塑性模型 R{"Kh2q_
TB,DP,1 !定义数据表1 c3,YA,skb!
TBDATA,1,9E5,42.8 !给数据表输入数据 F)_zR
!边坡围岩2材料属性 !采用弹性模型 +6 tj
w 6
MP,EX,2,3.2E10 ！定义弹性模量 k^ B<t'
MP,PRXY,2,0.24 !定义泊松比 x"q! =&>f
MP,DENS,2,2700 !定义密度 z+IBy+
!取强度折减系数F=1.2时边坡围岩材料属性 M:6Yy@#T.
MP,EX,3,3E10 K*^3FO}JG
MP,PRXY,3,0.25 NuZi LtC
MP,DENS,3,2500 (Tt\6-
TB,DP,3 !定义数据表3 qyzmjV6J2
TBDATA,1,7.5E5,37.7 !给数据表输入数据 z\wY3pIr2
!取强度折减系数F=1.4时边坡围岩材料属性 |P!7T.
MP,EX,4,3E10 qKu/~0a/
MP,PRXY,4,0.25 !定义泊松比 r:fMd3;gq
MP,DENS,4,2500 SPfz/ q{
TB,DP,4 !定义数据表4 *Ii_d pJ
TBDATA,1,6.4E5,33.5 !给数据表输入数据 PhTMXv<cE
!取强度折减系数F=1.6时边坡围岩材料属性 [?z`XY_-
MP,EX,5,3E10 ?`ETlFtD4
MP,PRXY,5,0.25 \
vf&Ldk
MP,DENS,5,2500 ！定义密度 Z4 y9d?g%b
TB,DP,5 !定义数据表5 KQG dV{VFs
TBDATA,1,5.6E5,30 !给数据表输入数据 vz'/] E
!取强度折减系数F=1.8时边坡围岩材料属性 }ZaZPB/_}P
MP,EX,6,3E10 ！定义弹性模量 (JV [7u -
MP,PRXY,6,0.25 3%gn:.9N
MP,DENS,6,2500 <y4hK3wP
TB,DP,6 !定义数据表6 [t #xX59
TBDATA,1,5.0E5,27.2 ?wM{NVt#-
!取强度折减系数F=2.0时边坡围岩材料属性 GHeucG} ?
MP,EX,7,3E10 4yy9m8 /
MP,PRXY,7,0.25 !定义泊松比 gE r @L
MP,DENS,7,2500
X0a)6HZ{
TB,DP,7 !定义数据表7 h=:Ls]ZU
TBDATA,1,4.5E5,24.8 j{ QzD^t
!取强度折减系数F=2.2时边坡围岩材料属性 /j%(Z/RM
MP,EX,8,3E10 5dwC~vn}c
MP,PRXY,8,0.25 :(7icHa
MP,DENS,8,2500 KV&4Ep#
TB,DP,8 !定义数据表8 9B W"^$
TBDATA,1,4.09E5,22.8 2z*EamF
!取强度折减系数F=2.4时边坡围岩材料属性 _,r2g8qm
MP,EX,9,3E10 cX~J6vNy5
MP,PRXY,9,0.25 VINb9W}G[
MP,DENS,9,2500 4Wsp PHj
TB,DP,9 !定义数据表9 VY9o}J>,w
TBDATA,1,3.6E5,21.1 !给数据表输入数据 E}-Y!,v^
!取强度折减系数F=2.6时边坡围岩材料属性 -d'swx2aZ!
MP,EX,10,3E10 ！定义弹性模量 M/<>'%sj
MP,PRXY,10,0.25 !定义泊松比 lqaOLZH
MP,DENS,10,2500 ！定义密度 $)or{Z$&
TB,DP,10 !定义数据表10 o'Rr2,lVi
TBDATA,1,3.46E5,19.6 !给数据表10输入数据 hB[VU
";
!取强度折减系数F=2.8时边坡围岩材料属性 dV/ ^@[
MP,EX,11,3E10 ！定义弹性模量 Twq/ Y07M
MP,PRXY,11,0.25 !定义泊松比 t[ =-4;
MP,DENS,11,2500 4$W}6 v
TB,DP,11 !定义数据表11 P]h-**O
TBDATA,1,3.2E5,18.2 !给数据表11输入数据 zDK"Y{
!取强度折减系数F=3.0时边坡围岩材料属性 <PMQ$s>KK
MP,EX,12,3E10 ！定义弹性模量 J_-K"T|f
MP,PRXY,12,0.25 qnO>F^itF
MP,DENS,12,2500 ！定义密度 #!yW)RG
TB,DP,12 !定义数据表12 UO~Xzx!e
TBDATA,1,3.0E5,17.1 !给数据表12输入数据 )\O;Rt(
!建立几何模型 "+dByaY
!创建边坡线模型线 *3R3C +
L
K,1,,,, ！ 创建关键点1 A!GQ4.~%
K,2,-800,,, X}f u $2
K,3,-800,-800,, BQ o$c~
K,4,-800,-1200,, H;D C kVL
K,5,1200,-1200,, bz'#YM
K,6,1200,-800,, TmE J!)*
K,7,1200,,, g.,_E4L
K,8,1200,400,, 2=]Xe#5J=
K,9,400,400,, Y%zYO
l,1,2 !连接关键点1、2生成直线L1 Y'&rSHI"
l,2,3 vRI0fDu
l,3,4 z79c30 y]"
l,4,5 BmhIKXE{*
l,5,6 k`(Cwp{Oc
l,6,7 c9/&A
l,7,8 #_p QS}$
l,8,9 9YQYg@+R
l,9,1 +eUWf{(_
l,1,7 8~!9bg6C
l,3,6 fUKi@*^ZUa
!创建边坡面模型 ^.g BHZ
al,3,4,5,11 !由线L3、L4、L5、L11生成一个面A1 Oaj$Z-
f
al,1,2,11,6,10 !由线L1、L2、L1、L6、L10生成一个面A2 ;${_eab ]
al,7,8,9,10 !由线L7、L8、L9、L10生成一个面A3 dTTC6?yPXf
/PNUM, area, on !打开面号开关 (8In f_ 59
!划分网格生成有限元模型 gt \O
! 划分边坡围岩2网格 Ge_Gx*R
!设置网格份数 xf|vz|J?y
lsel,s,,,3,5,2 !选择线L3和L5 JIOeDuw+
lesize,all,,,5 !把所选择线分为5段 c,L{Qv"n{
lsel,s,,,4,11,7 !选择线L4和L11 ]1)#Y
lesize,all,,,26 !把所选择线分为26段 ~ TfN*0
mat,2 !给边坡围岩2赋予2号材料特性 N~!
G AaD
type,1 !采用单元类型1 a<C J#B2K
mshkey,0 !设定自由网格划分 9%iv?/ o*L
mshape,0 !设定四边形网格划分 Ha41Wn'tZ
amesh,1 ！划分面积A1 -qs
R,H
! 划分边坡围岩1网格 ?#0m [k&`
!设置网格份数 *|*6 q /
lsel,s,,,2,6,4 !选择线L2和L6 ,t|qhJF
lsel,a,,,1 !附加选择线L1 v/G)E_
lesize,all,,,10 !把所选择线分为10段 @Wl2E.)K;
lsel,s,,,8,10,2 !选择线L8和L10 8{B]_:
-:
lesize,all,,,16 !把所选择线分为16段 :pw6#yi8`
lsel,s,,,7,9,2 ~{00moN"m
lesize,all,,,12 :8LK}TY7
mat,1 !给边坡围岩1赋予1号材料特性 Oh W o
type,1 [kKg?I$D@B
mshkey,0 !设定自由网格划分 Dd' 4W
mshape,0 !设定四边形网格划分 aB ^`3J
amesh,2 !划分面积A2 LUKt!I0l
amesh,3 ！划分面积A3 =u<jxV9
allsel M
$\!SXL
!保存网格模型 ,sZ)@?e
save,Slope-grid.db m#|h2 2^H
(3)施加约束和荷载 @!KG;d:l
!给边坡模型施加约束 3R ?6{.
!边坡两侧施加X方向约束 #q9cjEd_7
nsel,s,loc,x,-800 !选择X=-800线上所有节点 QdK
PzjA
nsel,a,loc,x,1200 !选择X=1200线上所有节点 b{(= C
3
d,all,ux !对所选择节点约束X方向位移 bFB.hkTP
allsel Xf%wW [~
!边坡底部施加约束 h
{ M=V
nsel,s,loc,y,-1200 !选择Y=--1200线上所有节点 K06 x7W
d,all,uy !对所选择节点约束X、Y方向位移 E HH+)mlo
d,all,ux KYM%U" j D
!施加重力加速度 _)a! g-Do7
acel,,9.8 F x3 X
(4)求解 b~Un=-@5a
/solu /D,<2>o
!求解设置 *T-v^ndJh
antype,static !设定为静力求解 uZXG"
nsubst,100 !设定最大子步数为100 ?F^O7\rw
pred,on !打开时间步长预测器 V2o1~R~
nropt,full !设定牛顿-拉普森选项 hwi$: [
nlgeom,on !打开大位移效果 !n`Y^
lnsrch,on !打开线性搜索 =y]F cxF
outres,all,all !输出所有项 q vGP$g
cnvtol,f,,0.005,2,0.5 !力收敛准则设定 owe6ge7m
cnvtol,u,,0.05,2,1 !位移收敛准则设定 kn= fW1
!边坡在强度折减系数F=1时求解 'tJ@+(tqw
allsel WKB
K)=
solve !进行求解 9/dI 6 P7
save,F1,db !把F=1时求解结果保存 3Bbd2[<W
!边坡在强度折减系数F=1.2时求解 n0vhc; d
finish ve*m\DU
/solu WVdF/H
allsel 0I7 r{T
asel,s,area,,2,3,1 !选择面积A2、A3 V]6CHE:BS
MPCHG,3,all ! 把所选择单元材料号改为3 D[Kq`
allsel .2x`Fj;o1
solve !进行求解 5 NN`tv
save,F1.2,db !把F=1.2时求解结果保存 4%>2 >5
!边坡在强度折减系数F=1.4时求解 Rd ,5 &X$
finish ij&T \):d
/solu G~{# %i
allsel wvPS0]
asel,s,area,,2,3,1 !选择面积A2、A3 f19'IH$n{
MPCHG,4,all !把所选择单元材料号改为4 d_ )VeuE2
allsel {(_> A\zi
solve !进行求解 IP+. L]S
save,F1.4,db !把F=1.4时求解结果保存 `s ~[q
!边坡在强度折减系数F=1.6时求解 ~!#2s'
finish OYj~"-3y)
/solu !"bU|a
allsel ~pPj
asel,s,area,,2,3,1 !选择面积A2、A3 66l+ cb
MPCHG,5,all !把所选择单元材料号改为5 =S< E[D{V`
allsel QRz5eGpW
solve c+ .?+g
save,F1.6,db !把F=1.6时求解结果保存 ; {Tf:j'g
!边坡在强度折减系数F=1.8时求解 arb'.:[z^
finish pg/SYEvs V
/solu SxW.dT8{
allsel 9k3RC}dEr
asel,s,area,,2,3,1 !选择面积A2、A3 KCi0v
MPCHG,6,all !把所选择单元材料号改为6 {LqahO*
allsel MxL i'R=
solve qV:TuR-|w
save,F1.8,db !把F=1.8时求解结果保存 I0v4T jHH
!边坡在强度折减系数F=2.0时求解 ^wMZG'/
finish g$^I/OK?
/solu B<!WAw+
allsel A*TO0L
asel,s,area,,2,3,1 !选择面积A2、A3 2A(IsUtqO:
MPCHG,7,all !把所选择单元材料号改为7 3<X*wVi)NN
allsel W7?f_E\>W
solve ^9
Pae)
save,F2.0,db !把F=2.0时求解结果保存 !xz{X ?
!边坡在强度折减系数F=2.2时求解 `+ 5,=S
finish ?.6fVSa
/solu p$}1V2h;
allsel C>Cb
asel,s,area,,2,3,1 !选择面积A2、A3 D UW SY?^c
MPCHG,8,all !把所选择单元材料号改为8 R[f@g;h
allsel \ 3F OI
solve (laVm U?I7
save,F2.2,db !把F=2.2时求解结果保存 d/l>~%bR
!边坡在强度折减系数F=2.4时求解 cp L '
finish BB3wG*q
/solu CNrK]+>
allsel o`~ %}3
asel,s,area,,2,3,1 !选择面积A2、A3 2M*84 oh8P
MPCHG,9,all !把所选择单元材料号改为9 4j}uVGi{e
allsel IpKpj"eoLy
solve /}R*'y
save,F2.4,db !把F=2.4时求解结果保存 _OvIi~KW+
!边坡在强度折减系数F=2.6时求解 f<Xi/ (
finish DW0UcLO
/solu qExmf%q:q
allsel zz1]6B*eX
asel,s,area,,2,3,1 !选择面积A2、A3
:qrCqFl
MPCHG,10,all !把所选择单元材料号改为10 u8Oo@xf0Fr
allsel TUM7(-,9
solve {O (@}
save,F2.6,db !把F=2.6时求解结果保存 V J]S"
!边坡在强度折减系数F=2.8时求解 =6< Am
finish ,]Q
i/m
/solu S|/Za".Gr
allsel M }qrF~
asel,s,area,,2,3,1 !选择面积A2、A3 #6F/:j;
MPCHG,11,all !把所选择单元材料号改为11 t{o&$s93
allsel N^xk.O_TO
solve
Vc?=cQ'c
save,F2.8,db !把F=2.8时求解结果保存 hV d PO
!边坡在强度折减系数F=3.0时求解 9P\R ?~3
finish jZP~!q
/solu tP0!TkTo9
allsel ^eY% T5K
asel,s,area,,2,3,1 !选择面积A2、A3 q+?>shqsZ
MPCHG,12,all !把所选择单元材料号改为12 ))eQZ3ap9
allsel Bxs0m]
solve $p~X"f? 0
save,F3.0,db !把F=3.0时求解结果保存 0mY KzJi
(5）后处理 \. _TOE9L
/post1 ！进入后处理 I'a&n}j x
!边坡在强度折减系数F=1时结果分析 )&wJ _ (z
Resume,'F1','db' !读入边坡在强度折减系数F=1时 >%Nqgn$V
set,1,last !读入后一个子步 2/.I6IbL
pldisp,1 !绘制边坡模型变形图 8|H^u6+yz
plnsol,u,x !绘制边坡模型水平方向位移云图 z<z \)
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 V;%DS)-
!边坡在强度折减系数F=1.2时结果分析 `cn}}1Lg]
Resume,'F1.2','db' !读入边坡在强度折减系数F=1.2时 .|x" '3 #
set,1,last !读入后一个子步 8D1+["&
pldisp,1 !绘制边坡模型变形图 oC >l|?h,
plnsol,u,x !绘制边坡模型水平方向位移云图 /d
prs(*K
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 O&ZVu>`g
!边坡在强度折减系数F=1.4时结果分析 m5\ /7 VC
Resume,'F1.4','db' !读入边坡在强度折减系数F=1.4时 U7le> d;L
set,1,last !读入后一个子步 >TB Rp,;r
pldisp,1 !绘制边坡模型变形图 cH{[\F"Eb
plnsol,u,x !绘制边坡模型水平方向位移云图 Mxk0XFA
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 " UxKG+
!边坡在强度折减系数F=1.6时结果分析 o}VW%G"
Resume,'F1.6','db' !读入边坡在强度折减系数F=1.6时 $hE,BeQ
set,1,last !读入后一个子步 04P! l
pldisp,1 !绘制边坡模型变形图 D vz}sQZ
plnsol,u,x !绘制边坡模型水平方向位移云图 HYWKx><
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 8V(-S,
!边坡在强度折减系数F=1.8时结果分析 Az4a| .
Resume,'F1.8','db' !读入边坡在强度折减系数F=1.8时 <WGx
6{
set,1,last !读入后一个子步 o%V%@q H
pldisp,1 !绘制边坡模型变形图 3*;S%1C^
plnsol,u,x !绘制边坡模型水平方向位移云图 L"ob ))GF
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 &H IG776
!边坡在强度折减系数F=2.0时结果分析 gRQV)8uh
Resume,'F2.0','db' !读入边坡在强度折减系数F=2.0时 wTK>U`o
set,1,last !读入后一个子步 &I=F4 z
pldisp,1 !绘制边坡模型变形图 V/CZcMY_
plnsol,u,x !绘制边坡模型水平方向位移云图 XX;MoE~MM
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 Sz3Tp5b
!边坡在强度折减系数F=2.2时结果分析 z3 Ro*yJU
Resume,'F2.2','db' !读入边坡在强度折减系数F=2.2时 0DgEOW9H
set,1,last !读入后一个子步 JtGBNz!"
pldisp,1 !绘制边坡模型变形图 6O#
xV:Uc<
plnsol,u,x !绘制边坡模型水平方向位移云图 aH+n]J]
=)
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 c</u]TD
!边坡在强度折减系数F=2.4时结果分析 (J,^)!g7
Resume,'F2.4','db' !读入边坡在强度折减系数F=2.4时 Gp5[H}8K
set,1,last !读入后一个子步 KQv97#n1
pldisp,1 !绘制边坡模型变形图 jo ^ +
plnsol,u,x !绘制边坡模型水平方向位移云图 %.`<ud
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 y`Km96 Ui
!边坡在强度折减系数F=2.6时结果分析 ~n ) <L7
Resume,'F2.6','db' !读入边坡在强度折减系数F=2.6时 3IHA+ Zz
set,1,last !读入后一个子步 I@9 [
pldisp,1 !绘制边坡模型变形图 3,4m|Z2)
plnsol,u,x !绘制边坡模型水平方向位移云图 R<FW?z*
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 b am*&E%0K
!边坡在强度折减系数F=2.8时结果分析 Y*q_>kps"
Resume,'F2.8','db' !读入边坡在强度折减系数F=2.8时 t9(sSl
set,1,last !读入后一个子步 g W(7jFl
pldisp,1 !绘制边坡模型变形图 T3#K uiwU9
plnsol,u,x !绘制边坡模型水平方向位移云图 Y[WL}:"93
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 zy!mP
!边坡在强度折减系数F=3.0时结果分析 .Bb86Y=3
Resume,'F3.0','db' !读入边坡在强度折减系数F=3.0时 8
DE%ot
set,1,last !读入后一个子步 "2a&G3}t"
pldisp,1 !绘制边坡模型变形图 Wp0L!X=0
plnsol,u,x !绘制边坡模型水平方向位移云图 T>\ r}p
plnsol,eppl,eqv !绘制边坡模型塑性应变云图
