各位朋友: 3%|&I:tI
大家好!我通过ansys用有限元强度折减法分析了一个边坡的稳定性,结果不同的工况下得出的塑性区是一样的,我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激!附命令流和塑性区如下:(1)创建物理环境 uZYF(Yu
/COM, Structural ! 指定结构分析 :kV#y
/TITLE, Slope stability Analysis ! 定义工作标题 <=&`ZH
/FILNAM,support,1 ! 定义工作文件名 I{&[[7H
(2)建立模型 QL/(72K
!进入前处理器 rXq.DvQ
/PREP7 FxY}m
!定义单元类型 Hio0HL-
ET,1,PLANE82 !定义边坡围岩单元 S+6.ZZ9c
KEYOPT,1,3,2 !求解类型为平面应变 (E3b\lST
!定义材料属性 nW:C/{n2tG
!边坡围岩1材料属性 =%O6:YM
MP,EX,1,3E10 fbvL7*
(
MP,PRXY,1,0.25 ~=LE0. 3[
MP,DENS,1,2500 !定义塑性模型 3M`M
TB,DP,1 !定义数据表1 1.hyCTnI
TBDATA,1,9E5,42.8 !给数据表输入数据 >|=ts
!边坡围岩2材料属性 !采用弹性模型 W*:.Gxv]
MP,EX,2,3.2E10 !定义弹性模量 ;;t yoh~t
MP,PRXY,2,0.24 !定义泊松比 4"ZP 'I;
MP,DENS,2,2700 !定义密度 LOYk9m
!取强度折减系数F=1.2时边坡围岩材料属性 G!##X: 6'
MP,EX,3,3E10 G.B2('
MP,PRXY,3,0.25 Rv>-4@fMJ
MP,DENS,3,2500 Q{>k1$fkV
TB,DP,3 !定义数据表3 wQf-sk#
TBDATA,1,7.5E5,37.7 !给数据表输入数据 ?z+eWL
!取强度折减系数F=1.4时边坡围岩材料属性 yPp9\[+^j
MP,EX,4,3E10 d-ko
^Y0
MP,PRXY,4,0.25 !定义泊松比 j;r-NCBnz
MP,DENS,4,2500 {Xy5pfW
Q
TB,DP,4 !定义数据表4 M3y NAN
TBDATA,1,6.4E5,33.5 !给数据表输入数据 _LnpnL:
!取强度折减系数F=1.6时边坡围岩材料属性 u#~RkY7s
MP,EX,5,3E10 |e&\<LwsP
MP,PRXY,5,0.25 3}1u\(Mf
MP,DENS,5,2500 !定义密度 pki%vRY
TB,DP,5 !定义数据表5 c6]U E@A
TBDATA,1,5.6E5,30 !给数据表输入数据 ^76]0`gS
!取强度折减系数F=1.8时边坡围岩材料属性 re<{
>
MP,EX,6,3E10 !定义弹性模量 N"R]Yp;j
MP,PRXY,6,0.25 ;(%QD
3 >
MP,DENS,6,2500 H?Wya.7
TB,DP,6 !定义数据表6 IOH}x4
TBDATA,1,5.0E5,27.2 ]EAO+x9
!取强度折减系数F=2.0时边坡围岩材料属性 l'qg8
MP,EX,7,3E10 UkC!1Jy
MP,PRXY,7,0.25 !定义泊松比 $PPi5f}HD
MP,DENS,7,2500 Zi
i
TB,DP,7 !定义数据表7 & .j&0WE
TBDATA,1,4.5E5,24.8 :cECRm*
!取强度折减系数F=2.2时边坡围岩材料属性 +sA2WK]
MP,EX,8,3E10 |df Pki{
MP,PRXY,8,0.25 33q}CzK
MP,DENS,8,2500 <lPG=Xt
TB,DP,8 !定义数据表8 _VN?#J)o
TBDATA,1,4.09E5,22.8 6 "sSo j
!取强度折减系数F=2.4时边坡围岩材料属性 N+xP26D8
MP,EX,9,3E10 @. l@\4m
MP,PRXY,9,0.25 "S]TP$O D
MP,DENS,9,2500 jr."I+
TB,DP,9 !定义数据表9 /6*42[r
TBDATA,1,3.6E5,21.1 !给数据表输入数据 e !Y~Qy
!取强度折减系数F=2.6时边坡围岩材料属性 P@B]
MP,EX,10,3E10 !定义弹性模量 x9g#<2w8
MP,PRXY,10,0.25 !定义泊松比 X_h}J=33Q
MP,DENS,10,2500 !定义密度 cI*;k.KU
TB,DP,10 !定义数据表10 2zb"MEOS5
TBDATA,1,3.46E5,19.6 !给数据表10输入数据 ki!0^t:9
!取强度折减系数F=2.8时边坡围岩材料属性 t*u:hex
MP,EX,11,3E10 !定义弹性模量
)10+@d
MP,PRXY,11,0.25 !定义泊松比 * u>\57W
MP,DENS,11,2500 7$=InK
TB,DP,11 !定义数据表11 KpGhQdR#
TBDATA,1,3.2E5,18.2 !给数据表11输入数据 niyV8v
!取强度折减系数F=3.0时边坡围岩材料属性 Hn"RH1Zy
MP,EX,12,3E10 !定义弹性模量 q} >%8;nm
MP,PRXY,12,0.25 F41=b4/
MP,DENS,12,2500 !定义密度 n>YKa)|W`
TB,DP,12 !定义数据表12 H
<l7ZS:
TBDATA,1,3.0E5,17.1 !给数据表12输入数据 TJRCH>E[a
!建立几何模型 0h_|t-9j
!创建边坡线模型线 T8g$uFo
K,1,,,, ! 创建关键点1 `;C V=,M
K,2,-800,,, H{wl% G
K,3,-800,-800,, _D tV
K,4,-800,-1200,, wHy!CP%
K,5,1200,-1200,, :I#V.
K,6,1200,-800,, ez$(c
K,7,1200,,, 1+s;FJ2}
K,8,1200,400,, [W&T(%(W-
K,9,400,400,, k(HUUH_z
l,1,2 !连接关键点1、2生成直线L1 WsB ?C&>x
l,2,3 7[)E>XRE
l,3,4 qE"OB
l,4,5 hVY$;s
l,5,6 2+XAX:YD
l,6,7 })%{AfDRF
l,7,8 MQ2}EY*A
l,8,9 9.#<b|g
l,9,1 AK4t\D)K1
l,1,7 <hyKu
l,3,6 GbI/4<)l}
!创建边坡面模型 59u}W 0
al,3,4,5,11 !由线L3、L4、L5、L11生成一个面A1 g_bLl)g<
al,1,2,11,6,10 !由线L1、L2、L1、L6、L10生成一个面A2 6gDN`e,@
al,7,8,9,10 !由线L7、L8、L9、L10生成一个面A3 L4W5EO$
/PNUM, area, on !打开面号开关 Pm7}"D'/
!划分网格生成有限元模型 /n&&Um\
! 划分边坡围岩2网格 FS O).=#
!设置网格份数 SXh-A1t
lsel,s,,,3,5,2 !选择线L3和L5 wCBplaojJ
lesize,all,,,5 !把所选择线分为5段 %|oym.-I6
lsel,s,,,4,11,7 !选择线L4和L11 ?@x/E&
lesize,all,,,26 !把所选择线分为26段 ;uGv:$([g
mat,2 !给边坡围岩2赋予2号材料特性 :3 mh@[V
type,1 !采用单元类型1 OyIw>Wfv
mshkey,0 !设定自由网格划分 Kg]J/|0\
mshape,0 !设定四边形网格划分 DEgXQ[
amesh,1 !划分面积A1 -9?]IIVb
! 划分边坡围岩1网格 QT}tvm@PMq
!设置网格份数 HzsdHH(J
lsel,s,,,2,6,4 !选择线L2和L6 x8|J-8A(
lsel,a,,,1 !附加选择线L1 y~V(aih}D
lesize,all,,,10 !把所选择线分为10段 *-X[u:
lsel,s,,,8,10,2 !选择线L8和L10 8ao _i=&x
lesize,all,,,16 !把所选择线分为16段 ca9X19NG
lsel,s,,,7,9,2 :`sUt1Fw.
lesize,all,,,12 \;Weizq5
mat,1 !给边坡围岩1赋予1号材料特性 kJR`:J3DJ
type,1 %C]>9."
mshkey,0 !设定自由网格划分 Fr-SvsNFB
mshape,0 !设定四边形网格划分 (8OsGn
amesh,2 !划分面积A2 *Pr )%
amesh,3 !划分面积A3 BA.uw_^4
allsel /-s6<e!
!保存网格模型 |s_GlJV.
save,Slope-grid.db ALHIGJW:6$
(3)施加约束和荷载 6]wIG$j
!给边坡模型施加约束 *
y,v}-
!边坡两侧施加X方向约束 !,PWb3S
nsel,s,loc,x,-800 !选择X=-800线上所有节点 j>kqz>3
nsel,a,loc,x,1200 !选择X=1200线上所有节点 +ZV5o&V>
d,all,ux !对所选择节点约束X方向位移 (@YG~0
allsel :(%5:1W
!边坡底部施加约束 ]M3yLYK/P
nsel,s,loc,y,-1200 !选择Y=--1200线上所有节点 zuCSj~
d,all,uy !对所选择节点约束X、Y方向位移 nk:)j:fr
d,all,ux O6Y0XL
!施加重力加速度 j<$2hiI/?&
acel,,9.8 I-)4YQI
(4)求解 cwL_tq
/solu Q2>gU#
!求解设置 :Dp0?&_
antype,static !设定为静力求解 pR_9NfV{
nsubst,100 !设定最大子步数为100 d;>QhoiL
pred,on !打开时间步长预测器 Bw.i}3UT6
nropt,full !设定牛顿-拉普森选项 Ys7]B9/1O
nlgeom,on !打开大位移效果 Y);=TM6s
lnsrch,on !打开线性搜索 $[|mGae
outres,all,all !输出所有项 Hr C+Yjp
cnvtol,f,,0.005,2,0.5 !力收敛准则设定 Vvo7C!$z
cnvtol,u,,0.05,2,1 !位移收敛准则设定 6u%&<")4HP
!边坡在强度折减系数F=1时求解 pCG}ZKa
allsel i#O SC5ZI
solve !进行求解 <3
uNl
save,F1,db !把F=1时求解结果保存 '%;m?t%q
!边坡在强度折减系数F=1.2时求解 naNghGQ
finish &,/S`ke=
/solu P<-@h1p,
allsel TA\vZGJ('
asel,s,area,,2,3,1 !选择面积A2、A3 MK*r+xfSae
MPCHG,3,all ! 把所选择单元材料号改为3 l-3~K-k<@
allsel xD 7]C|8o
solve !进行求解 ?]5qr?W%
save,F1.2,db !把F=1.2时求解结果保存 |y*c9
!边坡在强度折减系数F=1.4时求解 \7_y%HR
finish EX"yxZ~
/solu ^rz_f{c]-
allsel )%]J>&/0J
asel,s,area,,2,3,1 !选择面积A2、A3 D+c>F5
MPCHG,4,all !把所选择单元材料号改为4 <ZW-QN4
allsel YkADk9fE
solve !进行求解 +|f@^-
save,F1.4,db !把F=1.4时求解结果保存 iDD$pd,e\
!边坡在强度折减系数F=1.6时求解 fV~~J2IK
finish u@444Vzg
/solu QWU-m{@~&
allsel |[lKY+26:{
asel,s,area,,2,3,1 !选择面积A2、A3 UZ";a453r
MPCHG,5,all !把所选择单元材料号改为5 m[2gdJK
allsel ig"L\ C"T
solve tX[WH\(xI
save,F1.6,db !把F=1.6时求解结果保存
#Q5o)x
!边坡在强度折减系数F=1.8时求解 MOC/KNb
finish YZ7.1`8
/solu d=^z`nt !R
allsel v^PO|Z
asel,s,area,,2,3,1 !选择面积A2、A3 {G-kNU
MPCHG,6,all !把所选择单元材料号改为6 cb bFw
allsel d5 -qZ{W
solve WzWXE(
save,F1.8,db !把F=1.8时求解结果保存 yQrD9*t&g
!边坡在强度折减系数F=2.0时求解 (%9$! v{3
finish 0 {mex4
/solu Zd&S@Z
allsel kT=8e;K
asel,s,area,,2,3,1 !选择面积A2、A3 `_h&glMJ,q
MPCHG,7,all !把所选择单元材料号改为7 DwF hK*
allsel #E]59_
solve 4K74=r),i
save,F2.0,db !把F=2.0时求解结果保存 fy$1YI>!Q
!边坡在强度折减系数F=2.2时求解 d5d@k
finish R-Sym8c
/solu -qoH,4w
allsel JI5Dy>u:
asel,s,area,,2,3,1 !选择面积A2、A3 ^@]3R QB
MPCHG,8,all !把所选择单元材料号改为8 LQ% `c
allsel =H~j,K
solve }WV:erg`
save,F2.2,db !把F=2.2时求解结果保存 ^Dx&|UwiZa
!边坡在强度折减系数F=2.4时求解 w
= KPT''!
finish GthYzd:'hJ
/solu mc3"`+o
allsel .(vwIb8\_
asel,s,area,,2,3,1 !选择面积A2、A3 %)wjR/o
MPCHG,9,all !把所选择单元材料号改为9 ?Ob3tUz2
allsel EV?z`jE9
solve <GsuZ
save,F2.4,db !把F=2.4时求解结果保存 PW4q~rc=:
!边坡在强度折减系数F=2.6时求解 0$njMnB2l
finish #;<Y[hR{P
/solu ~K=b\xc^
allsel j!ch5A
asel,s,area,,2,3,1 !选择面积A2、A3 8i pez/
MPCHG,10,all !把所选择单元材料号改为10 Debv4Gr;^
allsel f 1d?.)
solve
bZ6+,J
save,F2.6,db !把F=2.6时求解结果保存 @XVTU
!边坡在强度折减系数F=2.8时求解 cnLro
finish
3CJwj
/solu cNH7C"@GVu
allsel {b{s<@?
asel,s,area,,2,3,1 !选择面积A2、A3 c%&>p||
MPCHG,11,all !把所选择单元材料号改为11 =Sv/IXX\di
allsel <uJ@:oWG7
solve qWw=8Bq
save,F2.8,db !把F=2.8时求解结果保存 YzWz|
!边坡在强度折减系数F=3.0时求解 Y ay?=Y{
finish P* o9a
/solu <}LC~B!
allsel q*KAk{kR(v
asel,s,area,,2,3,1 !选择面积A2、A3 qb4z
T
MPCHG,12,all !把所选择单元材料号改为12 Je{ykL?N
allsel v2?ZQeHr_(
solve h$*!8=M
save,F3.0,db !把F=3.0时求解结果保存 U~8g_*
(5)后处理 F\!
`/4
/post1 !进入后处理 u4j5w
!边坡在强度折减系数F=1时结果分析 Q20%"&Xp]
Resume,'F1','db' !读入边坡在强度折减系数F=1时 h\e.e3/
set,1,last !读入后一个子步 f5r0\7y0
pldisp,1 !绘制边坡模型变形图 D]}G.v1
plnsol,u,x !绘制边坡模型水平方向位移云图 {8OCXus3m
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 :-'qC8C
!边坡在强度折减系数F=1.2时结果分析 26nx`w?j(
Resume,'F1.2','db' !读入边坡在强度折减系数F=1.2时 f<H2-(m
set,1,last !读入后一个子步 ZW}_DT0
pldisp,1 !绘制边坡模型变形图 8_8l.!~
plnsol,u,x !绘制边坡模型水平方向位移云图 c|%6e(g"L
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 ;aBG,dr}i
!边坡在强度折减系数F=1.4时结果分析 ]tD]Wx%
Resume,'F1.4','db' !读入边坡在强度折减系数F=1.4时 v1[29t<I!
set,1,last !读入后一个子步 >/|*DI-HJ
pldisp,1 !绘制边坡模型变形图 o@Oqm> ]SS
plnsol,u,x !绘制边坡模型水平方向位移云图 HDLk>_N_s,
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 "fI6Cpc
!边坡在强度折减系数F=1.6时结果分析 0mnw{fE8_
Resume,'F1.6','db' !读入边坡在强度折减系数F=1.6时 2|L&DF:G
set,1,last !读入后一个子步 Y7aqO5
pldisp,1 !绘制边坡模型变形图 TN.rrop`#g
plnsol,u,x !绘制边坡模型水平方向位移云图 /\Ef%@
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 P'2Qen*
!边坡在强度折减系数F=1.8时结果分析 7"mc+QOp
Resume,'F1.8','db' !读入边坡在强度折减系数F=1.8时 w &(ag$p'
set,1,last !读入后一个子步 P%6~&woF
pldisp,1 !绘制边坡模型变形图 :
'c&,oLY
plnsol,u,x !绘制边坡模型水平方向位移云图 z}@7'_iJ
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 -%~4W?
!边坡在强度折减系数F=2.0时结果分析 NgwbQ7)
Resume,'F2.0','db' !读入边坡在强度折减系数F=2.0时 s>en
set,1,last !读入后一个子步 H. c7Nle
pldisp,1 !绘制边坡模型变形图 u"8yK5!
plnsol,u,x !绘制边坡模型水平方向位移云图 K8~d^G
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 y^k$Us
!边坡在强度折减系数F=2.2时结果分析 KP"+e:a%
Resume,'F2.2','db' !读入边坡在强度折减系数F=2.2时 j6YOKJX
set,1,last !读入后一个子步 G3 m Z($y
pldisp,1 !绘制边坡模型变形图 "zc l|@
plnsol,u,x !绘制边坡模型水平方向位移云图 nEfK53i_
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 <[v[ci
!边坡在强度折减系数F=2.4时结果分析
<Uur^uB
Resume,'F2.4','db' !读入边坡在强度折减系数F=2.4时 9z0p5)]n>
set,1,last !读入后一个子步 phK/
pldisp,1 !绘制边坡模型变形图 d1*<Ll9K
plnsol,u,x !绘制边坡模型水平方向位移云图 pJ=#zsE0
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 }#J/fa9
!
!边坡在强度折减系数F=2.6时结果分析 :Al!1BJQ
Resume,'F2.6','db' !读入边坡在强度折减系数F=2.6时 5bIw?%dk(
set,1,last !读入后一个子步 SKtr tm
pldisp,1 !绘制边坡模型变形图 !<oe=)Iz|
plnsol,u,x !绘制边坡模型水平方向位移云图 2/f}S?@
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 : +u]S2u{
!边坡在强度折减系数F=2.8时结果分析 &L:!VL{I
Resume,'F2.8','db' !读入边坡在强度折减系数F=2.8时 GVz6-T~\>
set,1,last !读入后一个子步 h 0|s
pldisp,1 !绘制边坡模型变形图 D@KlOU{<
plnsol,u,x !绘制边坡模型水平方向位移云图 pw#-_
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 [RL9>n8f
!边坡在强度折减系数F=3.0时结果分析 ,_P-$lB
Resume,'F3.0','db' !读入边坡在强度折减系数F=3.0时 BWNi [^]
set,1,last !读入后一个子步 lFkR=!?=
pldisp,1 !绘制边坡模型变形图 H::bwn`Vc
plnsol,u,x !绘制边坡模型水平方向位移云图 pdMc}=K
plnsol,eppl,eqv !绘制边坡模型塑性应变云图