各位朋友： 'w/ S6j
大家好！我通过ansys用有限元强度折减法分析了一个边坡的稳定性，结果不同的工况下得出的塑性区是一样的，我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激！附命令流和塑性区如下：(1)创建物理环境 Oq}7q!H
/COM, Structural ! 指定结构分析 t="nmjQs
/TITLE, Slope stability Analysis ! 定义工作标题 [ 'B u
/FILNAM,support,1 ! 定义工作文件名 ~xJr|_,gp
(2)建立模型 ]h`d>#Hw!
!进入前处理器 c|iTRco
/PREP7 Vhn Ir#L+
!定义单元类型 v[DxWs8q
ET,1,PLANE82 !定义边坡围岩单元 egWfKL&iy
KEYOPT,1,3,2 !求解类型为平面应变 ( enOj0
!定义材料属性 cp`ZeLz2^
!边坡围岩1材料属性 k
$# ,^)T
MP,EX,1,3E10 C(xsMO'k,,
MP,PRXY,1,0.25 uE %2kB*]
MP,DENS,1,2500 !定义塑性模型 #>z !ns
TB,DP,1 !定义数据表1 J'&K
TBDATA,1,9E5,42.8 !给数据表输入数据 4 4WyfpTJ*
!边坡围岩2材料属性 !采用弹性模型 4^ 0CHy
MP,EX,2,3.2E10 ！定义弹性模量 NUtKT~V
MP,PRXY,2,0.24 !定义泊松比 \^'-=8<* >
MP,DENS,2,2700 !定义密度 bXM&VW?OP
!取强度折减系数F=1.2时边坡围岩材料属性 `lbRy($L
MP,EX,3,3E10 \4fuC6d2
MP,PRXY,3,0.25 %w!x \U V
MP,DENS,3,2500 LS-_GslE7\
TB,DP,3 !定义数据表3 G8Ow;:Ro
TBDATA,1,7.5E5,37.7 !给数据表输入数据 F+D
e"^As
!取强度折减系数F=1.4时边坡围岩材料属性 KfC{/J\
MP,EX,4,3E10 e!k4 Ij-]
MP,PRXY,4,0.25 !定义泊松比 mZnsr@KF
MP,DENS,4,2500 }4,[oD
TB,DP,4 !定义数据表4 >V%.=})K
TBDATA,1,6.4E5,33.5 !给数据表输入数据 zSOZr2-
^a
!取强度折减系数F=1.6时边坡围岩材料属性 ^%U`|GBZp
MP,EX,5,3E10 \z@ :OR,
MP,PRXY,5,0.25 +t]Ge
>S
MP,DENS,5,2500 ！定义密度 Wrm3U/>e
TB,DP,5 !定义数据表5 R_:lp \S&
TBDATA,1,5.6E5,30 !给数据表输入数据 =|c7#GaiF
!取强度折减系数F=1.8时边坡围岩材料属性 au+:-Khm
MP,EX,6,3E10 ！定义弹性模量
z' 5
MP,PRXY,6,0.25 ]% G# x
MP,DENS,6,2500 ?cK67|% W
TB,DP,6 !定义数据表6 f$?`50D"1
TBDATA,1,5.0E5,27.2 x.I?)x!C'
!取强度折减系数F=2.0时边坡围岩材料属性 9zLeyw\
MP,EX,7,3E10 @ RdNAP_6
MP,PRXY,7,0.25 !定义泊松比 pG v*{.
MP,DENS,7,2500 gEgd/Le
TB,DP,7 !定义数据表7 &pv* TL 8
TBDATA,1,4.5E5,24.8 5RF*c,cNq
!取强度折减系数F=2.2时边坡围岩材料属性 \SJX;7 ST
MP,EX,8,3E10 .\
vrBf
MP,PRXY,8,0.25 3?+t%_[
MP,DENS,8,2500 K 'K/}q<
TB,DP,8 !定义数据表8 S[l z>I
TBDATA,1,4.09E5,22.8 O^3XhTW^\~
!取强度折减系数F=2.4时边坡围岩材料属性 ],zp~yVU&
MP,EX,9,3E10 aOUTKyR ~
MP,PRXY,9,0.25 95/; II
MP,DENS,9,2500 R3#| *)q
TB,DP,9 !定义数据表9 j0=H6 Y
TBDATA,1,3.6E5,21.1 !给数据表输入数据 ZxCXru1
!取强度折减系数F=2.6时边坡围岩材料属性 9`&sZ|"3
MP,EX,10,3E10 ！定义弹性模量 O /&Qzt
MP,PRXY,10,0.25 !定义泊松比 4jVd
MP,DENS,10,2500 ！定义密度 #! (2@N8
TB,DP,10 !定义数据表10 3]&le[.
TBDATA,1,3.46E5,19.6 !给数据表10输入数据 16~5 ;u
!取强度折减系数F=2.8时边坡围岩材料属性 J'wJe,
MP,EX,11,3E10 ！定义弹性模量 kFfNDM#D
MP,PRXY,11,0.25 !定义泊松比 @cG+ D
MP,DENS,11,2500 zvv/|z2(r
TB,DP,11 !定义数据表11 x|Ms2.!
TBDATA,1,3.2E5,18.2 !给数据表11输入数据 dL1{ i ,M
!取强度折减系数F=3.0时边坡围岩材料属性 (r1"!~d@
MP,EX,12,3E10 ！定义弹性模量 A(+V{1 L'
MP,PRXY,12,0.25 SEM- t
MP,DENS,12,2500 ！定义密度 b>}
)G7b}
TB,DP,12 !定义数据表12 Pn ?gB}l
TBDATA,1,3.0E5,17.1 !给数据表12输入数据 i\K88B&24
!建立几何模型 }JUc!cH8z
!创建边坡线模型线 ,n UovWN07
K,1,,,, ！ 创建关键点1 ')WS :\J
K,2,-800,,, Q[T)jo,j%
K,3,-800,-800,, 2UBAk')O}
K,4,-800,-1200,, D~2n8h"2ye
K,5,1200,-1200,, T- js*
K,6,1200,-800,, g6][N{xW0
K,7,1200,,, A#F6~QX(.9
K,8,1200,400,, c11;(
K,9,400,400,, u3jLe=Y'\
l,1,2 !连接关键点1、2生成直线L1 raMtTL+
l,2,3 BY$L[U;@T
l,3,4 4Le{| B
l,4,5 I5Rd ~-="G
l,5,6 qzu(4*Gk6
l,6,7 6>b#nFVJ
l,7,8 R0GD9
l,8,9 sei%QE]!/
l,9,1 '^' PdB
l,1,7 [E9_ZdB T
l,3,6 ?uF3Q)rCk
!创建边坡面模型 cNy*< Tv
al,3,4,5,11 !由线L3、L4、L5、L11生成一个面A1 R@IwmJxX
al,1,2,11,6,10 !由线L1、L2、L1、L6、L10生成一个面A2 l5enlYH
al,7,8,9,10 !由线L7、L8、L9、L10生成一个面A3 nbDjoZZ4
/PNUM, area, on !打开面号开关 @k-GyV-v
!划分网格生成有限元模型 [a#*%H{OC
! 划分边坡围岩2网格 ny<D1>{90
!设置网格份数 C5X! H_p
lsel,s,,,3,5,2 !选择线L3和L5 H<*n5r(c
lesize,all,,,5 !把所选择线分为5段 K j-zEl
lsel,s,,,4,11,7 !选择线L4和L11 5VGZ5,+<<
lesize,all,,,26 !把所选择线分为26段 +N|t:8 qaf
mat,2 !给边坡围岩2赋予2号材料特性 7e)j|a-!<
type,1 !采用单元类型1 ndvt
$*
mshkey,0 !设定自由网格划分 &~=d;llkT
mshape,0 !设定四边形网格划分 AFsYP/ g]
amesh,1 ！划分面积A1 LO%OH
u}]
! 划分边坡围岩1网格 W7\UZPs5t
!设置网格份数
_akpW
lsel,s,,,2,6,4 !选择线L2和L6 *4Z! 5iOs
lsel,a,,,1 !附加选择线L1 }C
J K9*Z
lesize,all,,,10 !把所选择线分为10段 r aR=k!3i
lsel,s,,,8,10,2 !选择线L8和L10 {J;(K~>?m
lesize,all,,,16 !把所选择线分为16段 aMxM3"
lsel,s,,,7,9,2 F]RZP/D`
lesize,all,,,12 ABq#I'H#@2
mat,1 !给边坡围岩1赋予1号材料特性 Yg;7TK y
type,1 :{-/b
mshkey,0 !设定自由网格划分 ;;432^jD
mshape,0 !设定四边形网格划分 HoZsDs.XZ
amesh,2 !划分面积A2 LS<*5 HWX
amesh,3 ！划分面积A3 x*:"G'z T
allsel ,jy9\n*<t9
!保存网格模型 u*T#? W?
save,Slope-grid.db $A98h -*x
(3)施加约束和荷载 8;3I:z&muQ
!给边坡模型施加约束 k+eeVy
!边坡两侧施加X方向约束 h,MaF<~
nsel,s,loc,x,-800 !选择X=-800线上所有节点
1<0Z @D~F
nsel,a,loc,x,1200 !选择X=1200线上所有节点 & i *e&{L7
d,all,ux !对所选择节点约束X方向位移 B2)5Z]
allsel B\~(:(OPM]
!边坡底部施加约束 6 ziBGU#.-
nsel,s,loc,y,-1200 !选择Y=--1200线上所有节点 QC1\Sn /
d,all,uy !对所选择节点约束X、Y方向位移 [E qZj /
d,all,ux WvcPOt8Bp>
!施加重力加速度 H00iy$R
acel,,9.8 :;& 3"-
(4)求解 QghL=
/solu 7lzmAih
!求解设置 {@ tO9pc`8
antype,static !设定为静力求解 , Mn`kL<F
nsubst,100 !设定最大子步数为100 t+Qx-sW
pred,on !打开时间步长预测器 %)o;2&aD
nropt,full !设定牛顿-拉普森选项 qt. =
nlgeom,on !打开大位移效果 LP? *RrM
lnsrch,on !打开线性搜索 \.=,}s V2Z
outres,all,all !输出所有项 z
E\~Oa;
cnvtol,f,,0.005,2,0.5 !力收敛准则设定 L~Xzo
cnvtol,u,,0.05,2,1 !位移收敛准则设定 1E$^ul-v
!边坡在强度折减系数F=1时求解 :M@#.
allsel V'l9fj*E
solve !进行求解 X09i +/ICK
save,F1,db !把F=1时求解结果保存 "Q[?W( SA
!边坡在强度折减系数F=1.2时求解 <4"Bb_U
finish ;F /w&u.n
/solu pO N@
allsel }l5Q0'
asel,s,area,,2,3,1 !选择面积A2、A3 Z ..s /K {
MPCHG,3,all ! 把所选择单元材料号改为3 T^2 o' _:
allsel 7K24sHw;%
solve !进行求解 q9nQ/]rkHF
save,F1.2,db !把F=1.2时求解结果保存 w!Ii
!边坡在强度折减系数F=1.4时求解 MX|@x~9W
finish `pd+ as
/solu _u#r;h[
allsel J
c:j7}OOV
asel,s,area,,2,3,1 !选择面积A2、A3 5^N` ~
MPCHG,4,all !把所选择单元材料号改为4 jZ<f-Ff 0
allsel WG&WPV/p
solve !进行求解 bZgFea_>i
save,F1.4,db !把F=1.4时求解结果保存 'pT8S
!边坡在强度折减系数F=1.6时求解 .ITTY QHv)
finish c:-n0m'i
/solu f Q f5%
allsel V~ QOl=`K:
asel,s,area,,2,3,1 !选择面积A2、A3 2:1
kSR^Ky
MPCHG,5,all !把所选择单元材料号改为5 L,sXJ23.
allsel A-u}&}l<
solve sQO>1bh
save,F1.6,db !把F=1.6时求解结果保存 07vzVsQ}p
!边坡在强度折减系数F=1.8时求解 yk2XfY
finish ?|GwuG8g
/solu W: 3fLXk+
allsel 0)9n${P7d
asel,s,area,,2,3,1 !选择面积A2、A3
&/)To
MPCHG,6,all !把所选择单元材料号改为6 $$T a
allsel @CA{uP ;
solve ZP<OyX?
save,F1.8,db !把F=1.8时求解结果保存 /Em6+DN>
!边坡在强度折减系数F=2.0时求解 sGGi7 %
finish 6D4 j];~X
/solu c u4 |!s`#
allsel 6PMu*-Nv!j
asel,s,area,,2,3,1 !选择面积A2、A3
3nx*M=
MPCHG,7,all !把所选择单元材料号改为7 [tJp^?6 *
allsel 58PL@H~@0
solve 6^z):d #u
save,F2.0,db !把F=2.0时求解结果保存 yDi'@Z9R?
!边坡在强度折减系数F=2.2时求解 !*,m=*[3
finish +"VXw2R_ e
/solu N1dM,H
allsel rpL]5e!
asel,s,area,,2,3,1 !选择面积A2、A3 E$4Ik.k
MPCHG,8,all !把所选择单元材料号改为8 d. y-R#F_]
allsel b Kr73S9
solve KN.WT aO
save,F2.2,db !把F=2.2时求解结果保存 0E^S!A 7
!边坡在强度折减系数F=2.4时求解 pH396GFIW
finish |_16IEJ
/solu 4B J w+EV8
allsel X D \;|
asel,s,area,,2,3,1 !选择面积A2、A3 V"A* B
MPCHG,9,all !把所选择单元材料号改为9 q)RTy|N J^
allsel i MF-TR
solve %)y-BdSp .
save,F2.4,db !把F=2.4时求解结果保存 w#>CYP`0k6
!边坡在强度折减系数F=2.6时求解 fLuOxYQb f
finish OB+QVYk"
/solu 6KX/Y j~B
allsel J/c5)IB|
asel,s,area,,2,3,1 !选择面积A2、A3 2))p B/
MPCHG,10,all !把所选择单元材料号改为10 [h0)V(1KR
allsel 1HeE$
solve Sh u=o weJ
save,F2.6,db !把F=2.6时求解结果保存 JiX-t\V ~
!边坡在强度折减系数F=2.8时求解 bG]?AiW r
finish Ol X
otp8
/solu 3Io7!:+
allsel wkD"EuW(
asel,s,area,,2,3,1 !选择面积A2、A3 xp]_>WGq
MPCHG,11,all !把所选择单元材料号改为11 I:] Pd
allsel :MF+`RpL
solve -g4 {:!*D
save,F2.8,db !把F=2.8时求解结果保存 9i!|w kx
!边坡在强度折减系数F=3.0时求解 Ka8Bed3
finish W'5c%SI
/solu 9gETWz(3I
allsel KW n .
asel,s,area,,2,3,1 !选择面积A2、A3 A3Vj3em
MPCHG,12,all !把所选择单元材料号改为12 .:Zb ~
allsel
^{64b
solve (l) r.Vj
save,F3.0,db !把F=3.0时求解结果保存 JzkI!5c< j
(5）后处理 Jwbb>mB!
/post1 ！进入后处理 nO8e'&|
!边坡在强度折减系数F=1时结果分析 1sXVuto
Resume,'F1','db' !读入边坡在强度折减系数F=1时 Ot s] y
set,1,last !读入后一个子步 > NtJ)N*
pldisp,1 !绘制边坡模型变形图 S\6.vw!'
plnsol,u,x !绘制边坡模型水平方向位移云图 ohPDknHp
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 8q|T`ac+N
!边坡在强度折减系数F=1.2时结果分析 KK /siG~O
Resume,'F1.2','db' !读入边坡在强度折减系数F=1.2时 D|/
4) ,v
set,1,last !读入后一个子步 2J t *s$
pldisp,1 !绘制边坡模型变形图 ^7Z.~ A y
plnsol,u,x !绘制边坡模型水平方向位移云图 F2',3
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 Y-]Ne"+v f
!边坡在强度折减系数F=1.4时结果分析 7"Q;Yi2(
Resume,'F1.4','db' !读入边坡在强度折减系数F=1.4时 vgKdhN2kI
set,1,last !读入后一个子步 b5l;bXp]
pldisp,1 !绘制边坡模型变形图 Gp1?drF6
plnsol,u,x !绘制边坡模型水平方向位移云图 <1kK@m -E
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 eMU t%zvb
!边坡在强度折减系数F=1.6时结果分析 I= 7 YAm[W
Resume,'F1.6','db' !读入边坡在强度折减系数F=1.6时 x#'v}(v
set,1,last !读入后一个子步 E&\ 0+-Dw
pldisp,1 !绘制边坡模型变形图 G@,XUP
plnsol,u,x !绘制边坡模型水平方向位移云图
R7Z!
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 b25C[C5C
!边坡在强度折减系数F=1.8时结果分析 piAFxS<6
Resume,'F1.8','db' !读入边坡在强度折减系数F=1.8时 ynZfO2kf
set,1,last !读入后一个子步 v. > 95|8
pldisp,1 !绘制边坡模型变形图 dK7BjZTJo
plnsol,u,x !绘制边坡模型水平方向位移云图 P?<G:]W
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 !eD
f}~
!边坡在强度折减系数F=2.0时结果分析 E7@m& R
Resume,'F2.0','db' !读入边坡在强度折减系数F=2.0时 :\|<7n
set,1,last !读入后一个子步 B \qu XE)
pldisp,1 !绘制边坡模型变形图 y~py+:_
plnsol,u,x !绘制边坡模型水平方向位移云图 1j!{?t ?
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 <p#+('N`
!边坡在强度折减系数F=2.2时结果分析 dz )(~@tgz
Resume,'F2.2','db' !读入边坡在强度折减系数F=2.2时 :a)` iJnb
set,1,last !读入后一个子步 #$ ,b )Uy
pldisp,1 !绘制边坡模型变形图 W9jxw 4)
plnsol,u,x !绘制边坡模型水平方向位移云图 W>E|I v[o
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 Vd A!tL
!边坡在强度折减系数F=2.4时结果分析 ,UW!?}@
Resume,'F2.4','db' !读入边坡在强度折减系数F=2.4时 ?Nl@K/
set,1,last !读入后一个子步 |G P1[Q{
pldisp,1 !绘制边坡模型变形图 4l_~-Peh
plnsol,u,x !绘制边坡模型水平方向位移云图 +LaR_n[
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 D3C3_
@ *
!边坡在强度折减系数F=2.6时结果分析 (CY#B%*
Resume,'F2.6','db' !读入边坡在强度折减系数F=2.6时 ;x-] 1 xx_
set,1,last !读入后一个子步 P:m6:F@hO
pldisp,1 !绘制边坡模型变形图 $kY ]HI
plnsol,u,x !绘制边坡模型水平方向位移云图 N[sJ5o F
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 +w (B9rH
!边坡在强度折减系数F=2.8时结果分析 R rp-SR?O
Resume,'F2.8','db' !读入边坡在强度折减系数F=2.8时 6f;20dn 6
set,1,last !读入后一个子步 A 7zL\U4
pldisp,1 !绘制边坡模型变形图 m@g9+ 7
plnsol,u,x !绘制边坡模型水平方向位移云图 G4J)o?:m@
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 o= C 'u
!边坡在强度折减系数F=3.0时结果分析 uVzvUz{b
Resume,'F3.0','db' !读入边坡在强度折减系数F=3.0时 4u7^v1/
set,1,last !读入后一个子步 yzyK$WN\[3
pldisp,1 !绘制边坡模型变形图 2h`Tn{&1/
plnsol,u,x !绘制边坡模型水平方向位移云图 U;FJSy
plnsol,eppl,eqv !绘制边坡模型塑性应变云图
