各位朋友： :X`Bc"
大家好！我通过ansys用有限元强度折减法分析了一个边坡的稳定性，结果不同的工况下得出的塑性区是一样的，我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激！附命令流和塑性区如下：(1)创建物理环境 +P~E54
/COM, Structural ! 指定结构分析 @ a1+
/TITLE, Slope stability Analysis ! 定义工作标题 ?'_Q^O>
/FILNAM,support,1 ! 定义工作文件名 z5CWgN
(2)建立模型 cXO_g!&2A
!进入前处理器 c !ybz{L
/PREP7 "/)}Cc,L
!定义单元类型 z.9
#AN=&[
ET,1,PLANE82 !定义边坡围岩单元 EuAJ .n
KEYOPT,1,3,2 !求解类型为平面应变 "KY9MBzPD
!定义材料属性 'ErtiD
!边坡围岩1材料属性 (\si/&
MP,EX,1,3E10 fU+A~oL%I
MP,PRXY,1,0.25 {GS7J
MP,DENS,1,2500 !定义塑性模型 ` NC{+A
TB,DP,1 !定义数据表1 }xl
@:Qo
TBDATA,1,9E5,42.8 !给数据表输入数据 nJTV@m XVq
!边坡围岩2材料属性 !采用弹性模型 ?^F#}>C
MP,EX,2,3.2E10 ！定义弹性模量 Ah2%LXdHA
MP,PRXY,2,0.24 !定义泊松比 * n)3y.s
MP,DENS,2,2700 !定义密度 T#1>p ED
!取强度折减系数F=1.2时边坡围岩材料属性 r"a4 ;&mf
MP,EX,3,3E10 } 3 1z
35
MP,PRXY,3,0.25 7^bO`
MP,DENS,3,2500 w@P c7$EP
TB,DP,3 !定义数据表3 5@+8*Fdk
TBDATA,1,7.5E5,37.7 !给数据表输入数据 Uv6#d":f;
!取强度折减系数F=1.4时边坡围岩材料属性 W`C&$v#
MP,EX,4,3E10 h-1eDxK6
MP,PRXY,4,0.25 !定义泊松比 _ "ysJ&
MP,DENS,4,2500 \jdp L1
TB,DP,4 !定义数据表4 : B=p%C
TBDATA,1,6.4E5,33.5 !给数据表输入数据 '\:?FQ
C
!取强度折减系数F=1.6时边坡围岩材料属性 XV2 f|8d>
MP,EX,5,3E10 f N8| 4
MP,PRXY,5,0.25 6 m 5\f
MP,DENS,5,2500 ！定义密度 ms=I lz
TB,DP,5 !定义数据表5 3ySP *J5
TBDATA,1,5.6E5,30 !给数据表输入数据 ;6o p|
!取强度折减系数F=1.8时边坡围岩材料属性 877>=Tp |
MP,EX,6,3E10 ！定义弹性模量 <R:KR(bT
MP,PRXY,6,0.25 n8zh;vuJ
MP,DENS,6,2500 w+ cI0lj
TB,DP,6 !定义数据表6 1rV?^5
TBDATA,1,5.0E5,27.2 ^P-!pK*
!取强度折减系数F=2.0时边坡围岩材料属性 3<x_[0v`K1
MP,EX,7,3E10 {Ex0mw)T
MP,PRXY,7,0.25 !定义泊松比 'M\ou}P
MP,DENS,7,2500 xA nAW
TB,DP,7 !定义数据表7 %S22[;v{N
TBDATA,1,4.5E5,24.8 cI3uH1;#
!取强度折减系数F=2.2时边坡围岩材料属性 z(^p@&r)F
MP,EX,8,3E10 V#W(c_g
MP,PRXY,8,0.25 |WeLmy%9
MP,DENS,8,2500 ,\5]n&T;r
TB,DP,8 !定义数据表8 ?-O(EY1E
TBDATA,1,4.09E5,22.8 ^/HE_keY
!取强度折减系数F=2.4时边坡围岩材料属性 uU`zbh}]L.
MP,EX,9,3E10 M i\f?
MP,PRXY,9,0.25 a pUV6h-v
MP,DENS,9,2500 mp~\ioI*d
TB,DP,9 !定义数据表9 17G7r\iNYq
TBDATA,1,3.6E5,21.1 !给数据表输入数据 C ,Je>G
!取强度折减系数F=2.6时边坡围岩材料属性 d]h[]Su/?
MP,EX,10,3E10 ！定义弹性模量 d}b# "A
MP,PRXY,10,0.25 !定义泊松比 n<7#?X7
MP,DENS,10,2500 ！定义密度 M`umfw T
TB,DP,10 !定义数据表10 `S Wf)1K
TBDATA,1,3.46E5,19.6 !给数据表10输入数据 +MOUO$;fGt
!取强度折减系数F=2.8时边坡围岩材料属性 kX {c+qH M
MP,EX,11,3E10 ！定义弹性模量 ^!|BKH8>f%
MP,PRXY,11,0.25 !定义泊松比 WKp Hb:H
MP,DENS,11,2500 6^['g-\2
TB,DP,11 !定义数据表11 K/A xojo
TBDATA,1,3.2E5,18.2 !给数据表11输入数据 ;W#G<M&n'
!取强度折减系数F=3.0时边坡围岩材料属性 8 k+Ctk
MP,EX,12,3E10 ！定义弹性模量 *#EyfMz-B
MP,PRXY,12,0.25 !.iA^D//]
MP,DENS,12,2500 ！定义密度 SZ c6=^$
TB,DP,12 !定义数据表12 _y `'T;~OY
TBDATA,1,3.0E5,17.1 !给数据表12输入数据 C,-V>bx g
!建立几何模型 1K,bmb xRt
!创建边坡线模型线 c*!bT$]~\
K,1,,,, ！ 创建关键点1 bdqo2ZO
K,2,-800,,, lN1T\
K,3,-800,-800,, z@&_3 Gl
K,4,-800,-1200,, bn^^|i
K,5,1200,-1200,, ;4XX8W1
K,6,1200,-800,, 33couAP#
K,7,1200,,, xJ%b<y{@
K,8,1200,400,, 8Md*9E#J("
K,9,400,400,, wmY6&^ ?uS
l,1,2 !连接关键点1、2生成直线L1 9VkuYm,3
l,2,3 x!!: jL'L
l,3,4 H5/%"1Q
l,4,5 l4u`R(!n5
l,5,6 &cDnZ3 Q;
l,6,7 RXgi> Hz
l,7,8 Q= ~e|
l,8,9 a&ZH
l,9,1 Czb@:l%sc
l,1,7 P 2;j>=W
l,3,6 w0moC9#$?
!创建边坡面模型 _}`iLA!$I
al,3,4,5,11 !由线L3、L4、L5、L11生成一个面A1 "xS",6Sy
al,1,2,11,6,10 !由线L1、L2、L1、L6、L10生成一个面A2 \I#2Mq?
al,7,8,9,10 !由线L7、L8、L9、L10生成一个面A3 Lt H;#Q
/PNUM, area, on !打开面号开关 XXmtpM8
!划分网格生成有限元模型 Aye!@RjM8
! 划分边坡围岩2网格 ^`=Z=C$fj
!设置网格份数 /5suyM=U
lsel,s,,,3,5,2 !选择线L3和L5 "S#0QH%5
lesize,all,,,5 !把所选择线分为5段 |!I#T
lsel,s,,,4,11,7 !选择线L4和L11 y"!+Fus9
lesize,all,,,26 !把所选择线分为26段 ykl./uY'
mat,2 !给边坡围岩2赋予2号材料特性 3,tKqR7g
type,1 !采用单元类型1 ]=q?= %H
mshkey,0 !设定自由网格划分 |...T
4:^Y
mshape,0 !设定四边形网格划分 e|AJxn]
amesh,1 ！划分面积A1 jbC7U9t7
! 划分边坡围岩1网格 HnioB=fc
!设置网格份数 v"_hWJ)
lsel,s,,,2,6,4 !选择线L2和L6 &hd +x5
lsel,a,,,1 !附加选择线L1 YG?W8)T
lesize,all,,,10 !把所选择线分为10段 <+sv7"a
lsel,s,,,8,10,2 !选择线L8和L10 #)=P/N1
lesize,all,,,16 !把所选择线分为16段 lGjmw"/C
lsel,s,,,7,9,2 Hc^b}A y7
lesize,all,,,12 Uh?SDay
mat,1 !给边坡围岩1赋予1号材料特性 T
-C2V$1
type,1 GIo&zPx
mshkey,0 !设定自由网格划分 5x4JDaG2
mshape,0 !设定四边形网格划分 H
<F6o-*
amesh,2 !划分面积A2 J9I!d.U
amesh,3 ！划分面积A3 v2JC{X qrI
allsel im%'S6_X4
!保存网格模型 B4[o nYU
save,Slope-grid.db -UPdgZ_Vxz
(3)施加约束和荷载 6/s#'#jh
!给边坡模型施加约束 R S;r
!边坡两侧施加X方向约束 x1$fkNu
nsel,s,loc,x,-800 !选择X=-800线上所有节点 aQ]C`9k
nsel,a,loc,x,1200 !选择X=1200线上所有节点 #= 7~.Y
d,all,ux !对所选择节点约束X方向位移 sqJ?dIBH
allsel #\@*C=
!边坡底部施加约束 ./[%%"
nsel,s,loc,y,-1200 !选择Y=--1200线上所有节点 cRT@Cu
d,all,uy !对所选择节点约束X、Y方向位移 2@:Go`mg
d,all,ux 5"^$3&)
!施加重力加速度 l5D8DvJCj
acel,,9.8 #Cvjv;
QwY
(4)求解 vy 1:>N?#5
/solu Po(9BRd7
!求解设置 \ tF><
antype,static !设定为静力求解 rMfp%DMA
nsubst,100 !设定最大子步数为100 8>6+]]O
pred,on !打开时间步长预测器 o}7`SYn
nropt,full !设定牛顿-拉普森选项 :s$ rD
nlgeom,on !打开大位移效果 %@kmuz??
lnsrch,on !打开线性搜索 V8`t7[r
outres,all,all !输出所有项 kVy%y"/
cnvtol,f,,0.005,2,0.5 !力收敛准则设定 >F!2 ib8
cnvtol,u,,0.05,2,1 !位移收敛准则设定 g G~UsA
!边坡在强度折减系数F=1时求解 4[Hf[ .
allsel qL ,!
solve !进行求解 \@GA;~x.b
save,F1,db !把F=1时求解结果保存 vM 1f-I-
!边坡在强度折减系数F=1.2时求解
. sg V
finish ,>` wz^ z
/solu @FdtM<X
allsel Ngi$y> {Sq
asel,s,area,,2,3,1 !选择面积A2、A3 k[gO>UGB;
MPCHG,3,all ! 把所选择单元材料号改为3 l`~*" 4|/
allsel mcpM<vY/H
solve !进行求解 c3Y\XzV3v
save,F1.2,db !把F=1.2时求解结果保存 68+ 9^
!边坡在强度折减系数F=1.4时求解 ^4_.5~(
finish P<
O [S
/solu Qg<_te)\
allsel ujmO'blO
asel,s,area,,2,3,1 !选择面积A2、A3 |.)oV;9
MPCHG,4,all !把所选择单元材料号改为4 a 5:YP
allsel kDS4 t?Ig
solve !进行求解 j0Kj>
save,F1.4,db !把F=1.4时求解结果保存 m/Q@-
!边坡在强度折减系数F=1.6时求解 AWi~qzTZ
finish |t!kD(~r
/solu k3}|^/bHJ
allsel ]j~V0 1p/e
asel,s,area,,2,3,1 !选择面积A2、A3 5|9,S
MPCHG,5,all !把所选择单元材料号改为5 u qXvN'Jr
allsel 4! XB?-.
solve M0S}-eXc5
save,F1.6,db !把F=1.6时求解结果保存 SQ1.jcWW[
!边坡在强度折减系数F=1.8时求解 JC_Y#kN@z
finish S c_*L<$
/solu @vCPX=c
allsel gieT kZ
asel,s,area,,2,3,1 !选择面积A2、A3 &BFW`5N
MPCHG,6,all !把所选择单元材料号改为6 !\z :S?V
allsel
B ;9^
solve ^j0Mu.+_
save,F1.8,db !把F=1.8时求解结果保存 ~kD/dXt
!边坡在强度折减系数F=2.0时求解 UMma|9l(i
finish /![S 3Ol
/solu *rXESw]BR
allsel ?7 6 Wg::
asel,s,area,,2,3,1 !选择面积A2、A3 0 gL]^_+7
MPCHG,7,all !把所选择单元材料号改为7 cWA9n}Z
allsel M-e!F+d{od
solve ^}8(o
save,F2.0,db !把F=2.0时求解结果保存 gah3d*d7
!边坡在强度折减系数F=2.2时求解 ) ~rf x
finish |ITp$ _S
/solu 4askQV &hj
allsel "
2Dz5L1v
asel,s,area,,2,3,1 !选择面积A2、A3 dpDVEEs84
MPCHG,8,all !把所选择单元材料号改为8 N&]v\MjI62
allsel M$B9?N6
solve lQ<2Vw#Yl
save,F2.2,db !把F=2.2时求解结果保存 +\fr3@Yc
!边坡在强度折减系数F=2.4时求解 =!*e; L
finish j #f+0
/solu /XeDN-{
allsel 0 k@4;BYu
asel,s,area,,2,3,1 !选择面积A2、A3 &B Y%<h0c
MPCHG,9,all !把所选择单元材料号改为9 V}. uF,>V
allsel d(3F:dbk
solve &na#ES $X,
save,F2.4,db !把F=2.4时求解结果保存 =;W"Pi;*
!边坡在强度折减系数F=2.6时求解 .0:BgM
finish pEuZsQ
/solu D^baXp8
allsel .{1G"(z
asel,s,area,,2,3,1 !选择面积A2、A3 zH0%;
o}
MPCHG,10,all !把所选择单元材料号改为10 yM}}mypS
allsel 9z$]hl
solve WS/^WxRY
save,F2.6,db !把F=2.6时求解结果保存 z\Rs?v"
!边坡在强度折减系数F=2.8时求解 GpMKOjVm|
finish `MA ee8u'
/solu HgvgO\`]
allsel 0&mo1 k_U
asel,s,area,,2,3,1 !选择面积A2、A3 @zL)R b%P$
MPCHG,11,all !把所选择单元材料号改为11 %0fF _OU
allsel `KqMcAW
solve Dd-;; Y1C
save,F2.8,db !把F=2.8时求解结果保存 +FfT)8@W
!边坡在强度折减系数F=3.0时求解 d rnqX-E;
finish 5+vCuVZ
/solu |NJe4lw+?
allsel iS&~oj_-%
asel,s,area,,2,3,1 !选择面积A2、A3 jV]'/X<
MPCHG,12,all !把所选择单元材料号改为12 ZM K"3c9
allsel }zY)H9J~
solve #s$b\"4
save,F3.0,db !把F=3.0时求解结果保存 q#:, 6HDd
(5）后处理 H% t/-'U?
/post1 ！进入后处理 O$k;p<?M
!边坡在强度折减系数F=1时结果分析 LZch7Xe3
Resume,'F1','db' !读入边坡在强度折减系数F=1时 1t_$pDF}
set,1,last !读入后一个子步 hb9e6Cc
pldisp,1 !绘制边坡模型变形图 Gtd!Y
x
plnsol,u,x !绘制边坡模型水平方向位移云图 )xX(Et6+`
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 9I0/KuZd
O
!边坡在强度折减系数F=1.2时结果分析 :y== O4
Resume,'F1.2','db' !读入边坡在强度折减系数F=1.2时 3 $ cDC8
set,1,last !读入后一个子步 =2] .G Gg
pldisp,1 !绘制边坡模型变形图 a*REx_gLG
plnsol,u,x !绘制边坡模型水平方向位移云图 ]W7(}~m
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 J~eY,n.6]
!边坡在强度折减系数F=1.4时结果分析 jb~a z
Resume,'F1.4','db' !读入边坡在强度折减系数F=1.4时 BF@(`D&>
set,1,last !读入后一个子步 )z&0 g2Am
pldisp,1 !绘制边坡模型变形图 \HLI
y
plnsol,u,x !绘制边坡模型水平方向位移云图 5LbU'5
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 A%> Ir`I
!边坡在强度折减系数F=1.6时结果分析 j / 4N
Resume,'F1.6','db' !读入边坡在强度折减系数F=1.6时 )8kcOBG^L
set,1,last !读入后一个子步 },KY9w
pldisp,1 !绘制边坡模型变形图 b Gq 0k&
plnsol,u,x !绘制边坡模型水平方向位移云图 @=,2 {JF*6
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 pJr c\`D
!边坡在强度折减系数F=1.8时结果分析 X&o!xV -+
Resume,'F1.8','db' !读入边坡在强度折减系数F=1.8时 7Fw`s@/%
set,1,last !读入后一个子步 sDT(3{)L7
pldisp,1 !绘制边坡模型变形图 RIOR%~ U
plnsol,u,x !绘制边坡模型水平方向位移云图 F, Y@
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 et(/`
!边坡在强度折减系数F=2.0时结果分析 zWoPa,
Resume,'F2.0','db' !读入边坡在强度折减系数F=2.0时 3v)v92;
set,1,last !读入后一个子步 vCyvy^s-I
pldisp,1 !绘制边坡模型变形图 R$' 4 d
plnsol,u,x !绘制边坡模型水平方向位移云图 m^rgzx19?
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 _ I8L#4\(=
!边坡在强度折减系数F=2.2时结果分析 E/|]xKG
Resume,'F2.2','db' !读入边坡在强度折减系数F=2.2时 5tT-[mQ*
set,1,last !读入后一个子步 s\i=-`
pldisp,1 !绘制边坡模型变形图 &sR=N60n
plnsol,u,x !绘制边坡模型水平方向位移云图 ;j])h !8X
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 k@JDG]R<{
!边坡在强度折减系数F=2.4时结果分析 <MZ$b aK
Resume,'F2.4','db' !读入边坡在强度折减系数F=2.4时 &dF$:$'s
set,1,last !读入后一个子步 Rn~FCj,-
pldisp,1 !绘制边坡模型变形图 vZj^&/F$=g
plnsol,u,x !绘制边坡模型水平方向位移云图 nv1'iSEeOl
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 oJe9H <
!边坡在强度折减系数F=2.6时结果分析 P1;T- .X~&
Resume,'F2.6','db' !读入边坡在强度折减系数F=2.6时 0* < gGC
set,1,last !读入后一个子步 L@2%a'
pldisp,1 !绘制边坡模型变形图 #c@Dn.W
plnsol,u,x !绘制边坡模型水平方向位移云图 ^prseO?A
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 6k uN)
!边坡在强度折减系数F=2.8时结果分析 &o{I9MD
Resume,'F2.8','db' !读入边坡在强度折减系数F=2.8时 La48M'u
set,1,last !读入后一个子步 pW7vY)hj
pldisp,1 !绘制边坡模型变形图 K&0op 4&
plnsol,u,x !绘制边坡模型水平方向位移云图 [R CUP.
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 Gc>bli<-
!边坡在强度折减系数F=3.0时结果分析 LWHP31{ R
Resume,'F3.0','db' !读入边坡在强度折减系数F=3.0时 5%"${ywI
set,1,last !读入后一个子步 ?z %@ ;&
pldisp,1 !绘制边坡模型变形图 9 P_`IsVK
plnsol,u,x !绘制边坡模型水平方向位移云图 hO(8v&ns3
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 vA@K b3 ,
s:lar4>kM
]2(vO0~
JIvVbI
QLH& WF
3dfG_a61y
qb(#{Sw0
