各位朋友： \ i<7Lk
大家好！我通过ansys用有限元强度折减法分析了一个边坡的稳定性，结果不同的工况下得出的塑性区是一样的，我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激！附命令流和塑性区如下：(1)创建物理环境 }2Y`Lr
/COM, Structural ! 指定结构分析 (''w$qq"D
/TITLE, Slope stability Analysis ! 定义工作标题 *,=8x\Shp
/FILNAM,support,1 ! 定义工作文件名 ? Y luX
(2)建立模型 ~b4kV)[ q
!进入前处理器 u
Qg$h S
/PREP7 8CH9&N5W5t
!定义单元类型 [Ov/&jD"
ET,1,PLANE82 !定义边坡围岩单元 :0bjP Qj
KEYOPT,1,3,2 !求解类型为平面应变 P=s3&NDD
!定义材料属性 4 `Jf_C
!边坡围岩1材料属性 ] 8 <`&~a
MP,EX,1,3E10 F4ylD5Y!
MP,PRXY,1,0.25 -av=5hm
MP,DENS,1,2500 !定义塑性模型 <KE%|6oER
TB,DP,1 !定义数据表1 / neY2D6
TBDATA,1,9E5,42.8 !给数据表输入数据 :X}Ie P
!边坡围岩2材料属性 !采用弹性模型 kX)*:~*
MP,EX,2,3.2E10 ！定义弹性模量 0+.<BOcW5
MP,PRXY,2,0.24 !定义泊松比 =p:~sn#
MP,DENS,2,2700 !定义密度 lc]cs D
!取强度折减系数F=1.2时边坡围岩材料属性 W[>qiYf^b
MP,EX,3,3E10 e-VGJxR
MP,PRXY,3,0.25 7=&+0@R#/d
MP,DENS,3,2500 5s> >]
.%
TB,DP,3 !定义数据表3 B ^{~,'
TBDATA,1,7.5E5,37.7 !给数据表输入数据 b%0@nu4
!取强度折减系数F=1.4时边坡围岩材料属性 b7gN|Hw5 H
MP,EX,4,3E10 ]]y,FQ,r
MP,PRXY,4,0.25 !定义泊松比 Zvra > %
MP,DENS,4,2500 K b'4W-&u!
TB,DP,4 !定义数据表4 LX =cx$ K
TBDATA,1,6.4E5,33.5 !给数据表输入数据 !HeQMz
!取强度折减系数F=1.6时边坡围岩材料属性 }U%T6~_wR
MP,EX,5,3E10 c}H}fyu%n
MP,PRXY,5,0.25 j'lC]}kH
MP,DENS,5,2500 ！定义密度 K1&t>2=%
TB,DP,5 !定义数据表5 _3#_6>=M
TBDATA,1,5.6E5,30 !给数据表输入数据 ",aEN=+|hV
!取强度折减系数F=1.8时边坡围岩材料属性 X}ihYM3y /
MP,EX,6,3E10 ！定义弹性模量 YHx Qb$v)
MP,PRXY,6,0.25 qt4%=E;[
MP,DENS,6,2500 :lK8i{o
TB,DP,6 !定义数据表6 f >&*%[fw
TBDATA,1,5.0E5,27.2 6^2='y~e
!取强度折减系数F=2.0时边坡围岩材料属性 aEun *V^,
MP,EX,7,3E10 ]Z52L`k
MP,PRXY,7,0.25 !定义泊松比 S@TfZ3Go|
MP,DENS,7,2500 <Ynrw4[)t
TB,DP,7 !定义数据表7 ~n(LBA
TBDATA,1,4.5E5,24.8 `\/\C[Gg
!取强度折减系数F=2.2时边坡围岩材料属性 VA %lJ!$
MP,EX,8,3E10 (CAkzgTfc
MP,PRXY,8,0.25 ?D(aky#cyc
MP,DENS,8,2500 ~aXJ5sY"f&
TB,DP,8 !定义数据表8 ,F+,A].wG
TBDATA,1,4.09E5,22.8 *)vy%\
!取强度折减系数F=2.4时边坡围岩材料属性 R0|4KT-i
MP,EX,9,3E10 7$ 8DMBqq
MP,PRXY,9,0.25 -M 4VC^_
MP,DENS,9,2500 =-qYp0sVP
TB,DP,9 !定义数据表9 U r8J G&,
TBDATA,1,3.6E5,21.1 !给数据表输入数据 ,| j\x
!取强度折减系数F=2.6时边坡围岩材料属性 KTeR;6oZn"
MP,EX,10,3E10 ！定义弹性模量 w @\4ft6d
MP,PRXY,10,0.25 !定义泊松比 Yj l:i*u/
MP,DENS,10,2500 ！定义密度 $I<\Yuy-M9
TB,DP,10 !定义数据表10 D u_ ;!E
TBDATA,1,3.46E5,19.6 !给数据表10输入数据 {!! 8 *ix
!取强度折减系数F=2.8时边坡围岩材料属性 (`R
heEg@f
MP,EX,11,3E10 ！定义弹性模量 _ x$\ E
MP,PRXY,11,0.25 !定义泊松比 j+>N&.zs
MP,DENS,11,2500 R 0G!5>1i
TB,DP,11 !定义数据表11 >X5RRSo
TBDATA,1,3.2E5,18.2 !给数据表11输入数据 X?(R!=a
!取强度折减系数F=3.0时边坡围岩材料属性 @4Q /J$
MP,EX,12,3E10 ！定义弹性模量 F;Q'R |HQ
MP,PRXY,12,0.25 _ D"S
MP,DENS,12,2500 ！定义密度 :8N{;aui
TB,DP,12 !定义数据表12 IYr}%:P)
TBDATA,1,3.0E5,17.1 !给数据表12输入数据 s{42_O?,c
!建立几何模型 H;[?8h(
!创建边坡线模型线 $+,kibk*R
K,1,,,, ！ 创建关键点1 g!)*CP#;
K,2,-800,,, ]O0u.=1k
K,3,-800,-800,, 'aS: Azb
K,4,-800,-1200,, |:SIyXGbY
K,5,1200,-1200,, Ix~rBD9
K,6,1200,-800,, Ds{DVdqA$c
K,7,1200,,, ME.!l6lm\
K,8,1200,400,, G+\&8fi0
K,9,400,400,, Ow v}lJ
l,1,2 !连接关键点1、2生成直线L1 8v=t -GJW
l,2,3 3{TE6&HIa
l,3,4 .&h|r>*|J
l,4,5 [ fs.D /
l,5,6 L [U?{
l,6,7 Fc|N6I'o
l,7,8 Nk1p)V SC
l,8,9 O5:[]vIn
l,9,1 \FnR'ne
l,1,7 dqi31e{*2\
l,3,6 \3hFb,/4k
!创建边坡面模型
=NWzsRl,
al,3,4,5,11 !由线L3、L4、L5、L11生成一个面A1 tJm1Q#||
al,1,2,11,6,10 !由线L1、L2、L1、L6、L10生成一个面A2 zQ(`pld
al,7,8,9,10 !由线L7、L8、L9、L10生成一个面A3 lg{M\
+
/PNUM, area, on !打开面号开关 u)%/df qzZ
!划分网格生成有限元模型 R~;8v 1>K
! 划分边坡围岩2网格 PtGFLM9R
!设置网格份数 ~y/qm
[P
lsel,s,,,3,5,2 !选择线L3和L5 ^S(QvoaQ
lesize,all,,,5 !把所选择线分为5段 DU-dIq i
lsel,s,,,4,11,7 !选择线L4和L11 HHOqJb{8S
lesize,all,,,26 !把所选择线分为26段 Z@AN0?,`~o
mat,2 !给边坡围岩2赋予2号材料特性 m;q qjzy
type,1 !采用单元类型1 AE Abny
q
mshkey,0 !设定自由网格划分 <L2z| %`
mshape,0 !设定四边形网格划分 =dp`4N
amesh,1 ！划分面积A1 R{6 M(!x
! 划分边坡围岩1网格 <@y(ikp>
!设置网格份数 `X B$t?xi
lsel,s,,,2,6,4 !选择线L2和L6 3ik~PgGoKQ
lsel,a,,,1 !附加选择线L1 }|nEbM]#
lesize,all,,,10 !把所选择线分为10段 at\$
IK_
lsel,s,,,8,10,2 !选择线L8和L10 urQ<r{$x0
lesize,all,,,16 !把所选择线分为16段 v9Lf |FXo&
lsel,s,,,7,9,2 k4` %.;
lesize,all,,,12 iT+t
mat,1 !给边坡围岩1赋予1号材料特性 lbh7`xCR
type,1 /XdLdA!v
mshkey,0 !设定自由网格划分 (%9J( 4
mshape,0 !设定四边形网格划分 zKh <zj
amesh,2 !划分面积A2 ViUx^e\
amesh,3 ！划分面积A3 L^2wEF
allsel hI*6f3Vn(n
!保存网格模型
lk=[Xo
save,Slope-grid.db Yqv!ZJ6
(3)施加约束和荷载 FS*J8)
!给边坡模型施加约束 "
^!=e72
!边坡两侧施加X方向约束 %H3 iX^}*
nsel,s,loc,x,-800 !选择X=-800线上所有节点 UgOhx- 8
nsel,a,loc,x,1200 !选择X=1200线上所有节点 []?*}o5&>T
d,all,ux !对所选择节点约束X方向位移 3@1$y`SN
allsel | `+ (O
!边坡底部施加约束 sXY{g0%
nsel,s,loc,y,-1200 !选择Y=--1200线上所有节点 ^Of\l:q*
d,all,uy !对所选择节点约束X、Y方向位移 lQ4$d{m`
d,all,ux 3(:?Z-iKe
!施加重力加速度 pezfB{x?
acel,,9.8 PK&X |
h
(4)求解 7lr;S(C
/solu .g.g lQ_~=
!求解设置 th5UzpB4
antype,static !设定为静力求解 Dj 3,SJ*x
nsubst,100 !设定最大子步数为100 Gc;-zq
pred,on !打开时间步长预测器 GKG:iR)
nropt,full !设定牛顿-拉普森选项 zXx A"
nlgeom,on !打开大位移效果 Ym $`EN
lnsrch,on !打开线性搜索 V8Z@y&ny
outres,all,all !输出所有项 l .m #
cnvtol,f,,0.005,2,0.5 !力收敛准则设定 ?iL-2I3*
cnvtol,u,,0.05,2,1 !位移收敛准则设定 0;%\L :,O
!边坡在强度折减系数F=1时求解 ly@%1
allsel }s8xr>
solve !进行求解 Wxi;Tq9C@_
save,F1,db !把F=1时求解结果保存 L\"eE'A
!边坡在强度折减系数F=1.2时求解 {#&D=7LP
finish uI3oPP> $
/solu fr8';Jm
allsel @[Wf!8_
asel,s,area,,2,3,1 !选择面积A2、A3 cVSns\QO
MPCHG,3,all ! 把所选择单元材料号改为3 INR RA
allsel },O7NSG<o
solve !进行求解 Qh]k)]+*|
save,F1.2,db !把F=1.2时求解结果保存 ]|[mwC4
!边坡在强度折减系数F=1.4时求解 7(H?3)%0
finish }$* z:E
/solu 46H@z =5
allsel [lz H%0
V
asel,s,area,,2,3,1 !选择面积A2、A3 }T5 3y6J#
MPCHG,4,all !把所选择单元材料号改为4 <d{>[R)
allsel 8sq0 BH
solve !进行求解 upq3)t_
save,F1.4,db !把F=1.4时求解结果保存 T`c:16I
!边坡在强度折减系数F=1.6时求解 -.Z y(
finish ft!D2M
/solu <<9|*Tz
allsel ) [=C@U
asel,s,area,,2,3,1 !选择面积A2、A3 M -8d*#_P
MPCHG,5,all !把所选择单元材料号改为5 WWLf'89It
allsel ;h#Q!M&e#
solve dx.J v/Mb
save,F1.6,db !把F=1.6时求解结果保存 %mOQIXr1s
!边坡在强度折减系数F=1.8时求解 dd4^4 X`j
finish <W*6=HZ'
/solu C k/DV
allsel WI SK-z
asel,s,area,,2,3,1 !选择面积A2、A3 s1X?]A
MPCHG,6,all !把所选择单元材料号改为6 ^xr & E
allsel X& XD2o"rt
solve Q{~;4+ZD
save,F1.8,db !把F=1.8时求解结果保存 gU?M/i2
!边坡在强度折减系数F=2.0时求解 B.) ;Ju
finish -y/Y% ]%0
/solu qporH]J-E
allsel Z e?H
asel,s,area,,2,3,1 !选择面积A2、A3 [5GzY`/m
MPCHG,7,all !把所选择单元材料号改为7 S5cs(}Bq
allsel
7uzc1}r
solve 0bu!(Tpg7
save,F2.0,db !把F=2.0时求解结果保存 :aD_>,n
!边坡在强度折减系数F=2.2时求解 s2#}@b6'.
finish w ; :{
/solu }G" bD8+
allsel :2~2j-m
asel,s,area,,2,3,1 !选择面积A2、A3 $`L
|
MPCHG,8,all !把所选择单元材料号改为8 ^ JU#_
allsel v}@Uc-(
solve "a<:fEsSE
save,F2.2,db !把F=2.2时求解结果保存 C~M,N|m+^
!边坡在强度折减系数F=2.4时求解 6hHMxS^o
finish ^vI`#}?
/solu O1oh,~W
allsel 41+@!`z7
asel,s,area,,2,3,1 !选择面积A2、A3 Yv[<c!\
MPCHG,9,all !把所选择单元材料号改为9 pQ8f$I#v
allsel =
jTC+0u
solve g c<Y?a-
save,F2.4,db !把F=2.4时求解结果保存 | }d+BD
!边坡在强度折减系数F=2.6时求解 c Hnd
gUW]
finish |"}rC >+
/solu <JW %h :\t
allsel r6nWrO>y
asel,s,area,,2,3,1 !选择面积A2、A3 }'uV{$
MPCHG,10,all !把所选择单元材料号改为10 ];u nR<H
allsel m-S e-aF
solve 6-\M }xq?
save,F2.6,db !把F=2.6时求解结果保存 (7X^ z&2
!边坡在强度折减系数F=2.8时求解 j<h0`v
finish ];QX&";Z
/solu j% '~l#nw
allsel ?VyiR40-Cx
asel,s,area,,2,3,1 !选择面积A2、A3 EEH TlqvR
MPCHG,11,all !把所选择单元材料号改为11 3+! G 9T!
allsel iq5-eJmq
solve p)~EG=p
save,F2.8,db !把F=2.8时求解结果保存 [] R8VC>Ah
!边坡在强度折减系数F=3.0时求解 4v`;D,dIu
finish 6L-3cxqf\
/solu o\nFSG kn
allsel Bey9P)_Of
asel,s,area,,2,3,1 !选择面积A2、A3 Z{_YH7_
MPCHG,12,all !把所选择单元材料号改为12 (?P\;yDG
allsel X$_z"t
solve Z_Y gV:jc
save,F3.0,db !把F=3.0时求解结果保存 _ ujhD
(5）后处理 yz%o?%@
/post1 ！进入后处理 mC'<Ov<eJ
!边坡在强度折减系数F=1时结果分析 8P}
a
Resume,'F1','db' !读入边坡在强度折减系数F=1时 RuOse9
set,1,last !读入后一个子步 <"7Wb"+
pldisp,1 !绘制边坡模型变形图 x,
'KI?TyQ
plnsol,u,x !绘制边坡模型水平方向位移云图 |doG}C
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 35fj-J$8
!边坡在强度折减系数F=1.2时结果分析 2>xEE
Resume,'F1.2','db' !读入边坡在强度折减系数F=1.2时 Qgq VbJP"
set,1,last !读入后一个子步 5zF$Q {3
pldisp,1 !绘制边坡模型变形图 ,F=FM> o
plnsol,u,x !绘制边坡模型水平方向位移云图
~vMJ?P@
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 ZD<e$PxxCd
!边坡在强度折减系数F=1.4时结果分析 O
2+taB
Resume,'F1.4','db' !读入边坡在强度折减系数F=1.4时 f~f)6XU|
set,1,last !读入后一个子步 =@d->d
pldisp,1 !绘制边坡模型变形图 _ F2ofB'
plnsol,u,x !绘制边坡模型水平方向位移云图 ~@"H\):/
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 5W09>C>OC
!边坡在强度折减系数F=1.6时结果分析 D+Z 2y1
Resume,'F1.6','db' !读入边坡在强度折减系数F=1.6时 i8 =+ <d
set,1,last !读入后一个子步 <qBM+m$|)
pldisp,1 !绘制边坡模型变形图 *~zB {
plnsol,u,x !绘制边坡模型水平方向位移云图 w`F'loUEt
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 OK
\9 `
!边坡在强度折减系数F=1.8时结果分析 >Xxi2V y
Resume,'F1.8','db' !读入边坡在强度折减系数F=1.8时 Sj vSnb_3
set,1,last !读入后一个子步 43!E> mq
pldisp,1 !绘制边坡模型变形图 R vd'uIJ
plnsol,u,x !绘制边坡模型水平方向位移云图 BfD C[(n`
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 L!Gpk) }[i
!边坡在强度折减系数F=2.0时结果分析 a@C}0IP)
Resume,'F2.0','db' !读入边坡在强度折减系数F=2.0时 CZ kmd
set,1,last !读入后一个子步 QH k jxj
pldisp,1 !绘制边坡模型变形图 Yd<9Y\W%?
plnsol,u,x !绘制边坡模型水平方向位移云图 perhR!#J
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 9e;:(jl^
!边坡在强度折减系数F=2.2时结果分析 eo&G@zwN
Resume,'F2.2','db' !读入边坡在强度折减系数F=2.2时 $kxu-
set,1,last !读入后一个子步 m =60a@o]
pldisp,1 !绘制边坡模型变形图 g2YE^EKU~
plnsol,u,x !绘制边坡模型水平方向位移云图 c7]0 >nU;
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 9x#T j/5%
!边坡在强度折减系数F=2.4时结果分析 ?:+p#&I
Resume,'F2.4','db' !读入边坡在强度折减系数F=2.4时 Am >b 7Z!
set,1,last !读入后一个子步 r>6FJ:T x
pldisp,1 !绘制边坡模型变形图 ]#W9l\
plnsol,u,x !绘制边坡模型水平方向位移云图 }eSrJgF4M
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 &3\3wcZ,q
!边坡在强度折减系数F=2.6时结果分析 jEL"Q?#
Resume,'F2.6','db' !读入边坡在强度折减系数F=2.6时 3s#/d,+
set,1,last !读入后一个子步 {v2[x W
pldisp,1 !绘制边坡模型变形图 Ys<z%
plnsol,u,x !绘制边坡模型水平方向位移云图 `KieN/d%
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 s@*i
!边坡在强度折减系数F=2.8时结果分析 X#ud_+6x
Resume,'F2.8','db' !读入边坡在强度折减系数F=2.8时 B_"P FWwg
set,1,last !读入后一个子步 @u:q#b
pldisp,1 !绘制边坡模型变形图 ~bgM*4GW
plnsol,u,x !绘制边坡模型水平方向位移云图 6|1*gl1_LD
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 `6v24?z
!边坡在强度折减系数F=3.0时结果分析 Tzfk_h3hE
Resume,'F3.0','db' !读入边坡在强度折减系数F=3.0时 |W5lhx0U
set,1,last !读入后一个子步 i({MID)/_
pldisp,1 !绘制边坡模型变形图 cEK#5
plnsol,u,x !绘制边坡模型水平方向位移云图 P9M%B2DQ6f
plnsol,eppl,eqv !绘制边坡模型塑性应变云图
