各位朋友： i8EMjLBUR
大家好！我通过ansys用有限元强度折减法分析了一个边坡的稳定性，结果不同的工况下得出的塑性区是一样的，我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激！附命令流和塑性区如下：(1)创建物理环境 |>d5 6
/COM, Structural ! 指定结构分析 (?$}Vp
/TITLE, Slope stability Analysis ! 定义工作标题 (L
q^C=
/FILNAM,support,1 ! 定义工作文件名 3d
\bB !
(2)建立模型 A5 /Q:8b
!进入前处理器 l'mgjv~
/PREP7 q d(`~a
!定义单元类型 & [4 Gv61
ET,1,PLANE82 !定义边坡围岩单元 )6Hc Pso6
KEYOPT,1,3,2 !求解类型为平面应变 F~&bgl[YZ
!定义材料属性 U
Tw\_s
!边坡围岩1材料属性 ECQ>VeP
MP,EX,1,3E10 `@-H
;
MP,PRXY,1,0.25 -2bu`oD
`
MP,DENS,1,2500 !定义塑性模型 ;CF I*Wfp
TB,DP,1 !定义数据表1 Pij*?qmeQ
TBDATA,1,9E5,42.8 !给数据表输入数据 V"k * PLt
!边坡围岩2材料属性 !采用弹性模型 ^ WidA-
MP,EX,2,3.2E10 ！定义弹性模量 @G^
l`%
MP,PRXY,2,0.24 !定义泊松比 L%I@HB9-Q0
MP,DENS,2,2700 !定义密度 p
0R)Yc+;
!取强度折减系数F=1.2时边坡围岩材料属性 /tDwgxJ
MP,EX,3,3E10 P**h\+M>{
MP,PRXY,3,0.25 /+iU1m'(
MP,DENS,3,2500 -NH c~ =m
TB,DP,3 !定义数据表3 #r >
TBDATA,1,7.5E5,37.7 !给数据表输入数据 xyBWV]Y
!取强度折减系数F=1.4时边坡围岩材料属性 ?/5WM%
MP,EX,4,3E10 vm
1vX;
MP,PRXY,4,0.25 !定义泊松比 K+2<{qwh
MP,DENS,4,2500 ):kDWc
TB,DP,4 !定义数据表4 &Z#g/Hc
TBDATA,1,6.4E5,33.5 !给数据表输入数据 gT?:zd=;
!取强度折减系数F=1.6时边坡围岩材料属性 auK*\Wjm?
MP,EX,5,3E10 &O7]e3E j
MP,PRXY,5,0.25 !S$LRm\ '
MP,DENS,5,2500 ！定义密度 "NMX>a,(
TB,DP,5 !定义数据表5 tg.[.v Ks
TBDATA,1,5.6E5,30 !给数据表输入数据 o@XhL9
!取强度折减系数F=1.8时边坡围岩材料属性 MB%yC]w8
MP,EX,6,3E10 ！定义弹性模量 #px74EeI\
MP,PRXY,6,0.25 OXT 5
y)
MP,DENS,6,2500 H0LEK(K
TB,DP,6 !定义数据表6 |W:xbtPNy
TBDATA,1,5.0E5,27.2 ot8UuBq
!取强度折减系数F=2.0时边坡围岩材料属性 <My4 )3
MP,EX,7,3E10 M4~^tML>Ey
MP,PRXY,7,0.25 !定义泊松比 j.= VZ
MP,DENS,7,2500 R m>AU=
TB,DP,7 !定义数据表7 F ^f L
TBDATA,1,4.5E5,24.8 g\ilK:r}
!取强度折减系数F=2.2时边坡围岩材料属性 w|lA%H7`J
MP,EX,8,3E10 e5W 8YNA
MP,PRXY,8,0.25 WzDL(~m+Z
MP,DENS,8,2500 3"!h+dXw
TB,DP,8 !定义数据表8 h/?$~OD
TBDATA,1,4.09E5,22.8 RoS&oGYqR
!取强度折减系数F=2.4时边坡围岩材料属性 }%x 2Z{VF
MP,EX,9,3E10 i Glg@
MP,PRXY,9,0.25 mu@He&w"
MP,DENS,9,2500 Bi:%}8STH
TB,DP,9 !定义数据表9 k[#<=G_=/E
TBDATA,1,3.6E5,21.1 !给数据表输入数据 aQzx^%B1
!取强度折减系数F=2.6时边坡围岩材料属性 KxhMPvN'
MP,EX,10,3E10 ！定义弹性模量 mmRxs1 0$
MP,PRXY,10,0.25 !定义泊松比 %{I b
MP,DENS,10,2500 ！定义密度 (_1(<Jw
TB,DP,10 !定义数据表10 "tdF #>x
TBDATA,1,3.46E5,19.6 !给数据表10输入数据 A6eIf
!取强度折减系数F=2.8时边坡围岩材料属性 0 w Q'~8
MP,EX,11,3E10 ！定义弹性模量 }$
C;ccWL
MP,PRXY,11,0.25 !定义泊松比 =#L\fe)q)
MP,DENS,11,2500 5e1 ;m6
TB,DP,11 !定义数据表11 vLCyT=OB`
TBDATA,1,3.2E5,18.2 !给数据表11输入数据 ||yx?q6\h
!取强度折减系数F=3.0时边坡围岩材料属性 w5Lev}Rb
MP,EX,12,3E10 ！定义弹性模量 \USl 9*E
MP,PRXY,12,0.25 UZvF5Hoe+O
MP,DENS,12,2500 ！定义密度 uC$!|I
TB,DP,12 !定义数据表12 #$n >+ lc
TBDATA,1,3.0E5,17.1 !给数据表12输入数据 )seeBm-`
!建立几何模型 BRSI g]
!创建边坡线模型线 X/ Sp!W-H
K,1,,,, ！ 创建关键点1 _Fe%Ek1Yy
K,2,-800,,, <UsFB F
K,3,-800,-800,, =^ZDP1h/}
K,4,-800,-1200,, kZ 5; Fe\*
K,5,1200,-1200,, KJ (|skO
K,6,1200,-800,, ? $$Xg3w_#
K,7,1200,,, 9 4.|l
K,8,1200,400,, _"l2UDx
K,9,400,400,, l/xpAx
l,1,2 !连接关键点1、2生成直线L1 qL2!\zt>g
l,2,3 =Vi e0TV&h
l,3,4 -< 0xS.^
l,4,5 )>"|<h.2]
l,5,6 D*0[7:NSO
l,6,7 ]Xm+-{5?!R
l,7,8 U\y:\+e l
l,8,9 }vK8P r%
l,9,1 \gC h'3
l,1,7 BFnp[93N
l,3,6 1L; 3e@G
!创建边坡面模型 w 3kX!%a:
al,3,4,5,11 !由线L3、L4、L5、L11生成一个面A1 >bR oQ 8
al,1,2,11,6,10 !由线L1、L2、L1、L6、L10生成一个面A2 3kz O
VZ
al,7,8,9,10 !由线L7、L8、L9、L10生成一个面A3 V
|cPAT%
/PNUM, area, on !打开面号开关 2p %j@O
!划分网格生成有限元模型 zQ~N(Jj?h
! 划分边坡围岩2网格 z&wJ"[nOC
!设置网格份数 kY?w] lS)t
lsel,s,,,3,5,2 !选择线L3和L5 WBN3:Y7
lesize,all,,,5 !把所选择线分为5段 A'6-E{
lsel,s,,,4,11,7 !选择线L4和L11 > q&X#E<w
lesize,all,,,26 !把所选择线分为26段 d D351 !-
mat,2 !给边坡围岩2赋予2号材料特性 P(xgIMc H
type,1 !采用单元类型1 tqD=)0Uzs
mshkey,0 !设定自由网格划分 %p;;aZG
mshape,0 !设定四边形网格划分 "W^+NeLc
amesh,1 ！划分面积A1 w<LV5w+
! 划分边坡围岩1网格 imc1rY!~'
!设置网格份数 `c)[aP{vN
lsel,s,,,2,6,4 !选择线L2和L6 O>y*u 8
lsel,a,,,1 !附加选择线L1 XOL_vS24
lesize,all,,,10 !把所选择线分为10段 [Smqe>U 1
lsel,s,,,8,10,2 !选择线L8和L10 7u`:e,'
lesize,all,,,16 !把所选择线分为16段 +aQM %~
lsel,s,,,7,9,2 ]bP1gV(b-
lesize,all,,,12 }v xRjO,
mat,1 !给边坡围岩1赋予1号材料特性 ) vD:
type,1 XwZR
Kh\>=
mshkey,0 !设定自由网格划分 j4wsDtmAU
mshape,0 !设定四边形网格划分 PR3i }y>
amesh,2 !划分面积A2 a9_KoOa.H
amesh,3 ！划分面积A3 cv5+[;(b
allsel H0\' , X
!保存网格模型 =d BK,/
save,Slope-grid.db .4J7 ^l
(3)施加约束和荷载 ujW C!*W(Q
!给边坡模型施加约束 HDi_|{2^
!边坡两侧施加X方向约束 cO8yu`4!e
nsel,s,loc,x,-800 !选择X=-800线上所有节点 >H5 _,A}f
nsel,a,loc,x,1200 !选择X=1200线上所有节点 m1D,#=C,_
d,all,ux !对所选择节点约束X方向位移 OUhlQq\
allsel T@xaa\bzg
!边坡底部施加约束 4(? Z1S
nsel,s,loc,y,-1200 !选择Y=--1200线上所有节点 R$x(3ey x
d,all,uy !对所选择节点约束X、Y方向位移 LFAefl\
d,all,ux Ca["tks
!施加重力加速度 KBDNK_7A
acel,,9.8 ]+5Y\~I
(4)求解 Mtaky=l8~I
/solu 1RauI0d*
!求解设置 m KJO?7tj
antype,static !设定为静力求解 p
~pl|
nsubst,100 !设定最大子步数为100 I6i qC"BK
pred,on !打开时间步长预测器 D}rnp wp{
nropt,full !设定牛顿-拉普森选项 s5 4AM]a{j
nlgeom,on !打开大位移效果 XnKf<|j6k
lnsrch,on !打开线性搜索 P N(<=v&E
outres,all,all !输出所有项
qN'%q+ n
cnvtol,f,,0.005,2,0.5 !力收敛准则设定 06 an(& a9
cnvtol,u,,0.05,2,1 !位移收敛准则设定 ; EsfH Ci)
!边坡在强度折减系数F=1时求解 pmiC|F83!8
allsel V#+M lN
solve !进行求解 MKh}2B#S
save,F1,db !把F=1时求解结果保存 4}t$Lf_
!边坡在强度折减系数F=1.2时求解 S;SI#Vg@
finish JSoInR1E
/solu ~H4Tr[8a
allsel YRu/KUT$ 7
asel,s,area,,2,3,1 !选择面积A2、A3 * $|9e
MPCHG,3,all ! 把所选择单元材料号改为3 g*W Y kv
allsel ~X-v@a
solve !进行求解 %~\I*v04
save,F1.2,db !把F=1.2时求解结果保存 :cA8[!
!边坡在强度折减系数F=1.4时求解 bwK1XlfD.s
finish &r!jjT
/solu $'w>doUlA
allsel ;y;UgwAM
asel,s,area,,2,3,1 !选择面积A2、A3 fM3ZoH/
MPCHG,4,all !把所选择单元材料号改为4 gMPvzBpP
allsel ^V"0 8
solve !进行求解 POQ4&ChA
save,F1.4,db !把F=1.4时求解结果保存 +vYV x<uTQ
!边坡在强度折减系数F=1.6时求解 c^~R %Bx
finish b~ ;M&Y
/solu {A\y 4D@
allsel |H%,>r`9S
asel,s,area,,2,3,1 !选择面积A2、A3 Nkx W*w%}l
MPCHG,5,all !把所选择单元材料号改为5 2 U3WH.o
allsel vv D515i
solve in> ?kbaG+
save,F1.6,db !把F=1.6时求解结果保存
A/OGF>
!边坡在强度折减系数F=1.8时求解 eMjW^-RgE5
finish k^%F4d3z@C
/solu .G>6_n3
allsel [ z { }?
asel,s,area,,2,3,1 !选择面积A2、A3 FXbalQ?^
MPCHG,6,all !把所选择单元材料号改为6 }q)dXFL=I#
allsel !Xzy:
solve RCMO?CBe
save,F1.8,db !把F=1.8时求解结果保存 `8D'r|=`Eh
!边坡在强度折减系数F=2.0时求解 +2m\Sv V
finish 4EJ6Zy![0*
/solu J6J;
!~>_
allsel K>Tv M &
asel,s,area,,2,3,1 !选择面积A2、A3
?dvcmXR
MPCHG,7,all !把所选择单元材料号改为7 960rbxKy3
allsel u:>*~$f
solve Gu%}B@ 4^
save,F2.0,db !把F=2.0时求解结果保存 dLw,dg
!边坡在强度折减系数F=2.2时求解 vl5r ~F
finish @|}=W Q
/solu +?Y(6$o
allsel eSSv8 [u
asel,s,area,,2,3,1 !选择面积A2、A3 /&RS+By(i
MPCHG,8,all !把所选择单元材料号改为8 =
$6pL
allsel 2w:cdAv$
solve #+;0=6+SM
save,F2.2,db !把F=2.2时求解结果保存 ^S^7 u
!边坡在强度折减系数F=2.4时求解 ys9MV%*
finish bAL!l\&2
/solu G=>LW1E|
allsel #po}Y
asel,s,area,,2,3,1 !选择面积A2、A3 lLZ?&z$
MPCHG,9,all !把所选择单元材料号改为9 k^\>=JTq=
allsel Q".AmHn
solve Nk shJ2
save,F2.4,db !把F=2.4时求解结果保存 xJ)vfo
!边坡在强度折减系数F=2.6时求解 +WSM< S2 U
finish yQ+C}8r5
/solu 6E.64+P Jw
allsel oBz jEv
asel,s,area,,2,3,1 !选择面积A2、A3 lHSu T2)x;
MPCHG,10,all !把所选择单元材料号改为10 H_7X%TvXb
allsel G `JXi/#`
solve >=n g ?
save,F2.6,db !把F=2.6时求解结果保存 sme!!+Rd
!边坡在强度折减系数F=2.8时求解 m>-(c=3
finish g}'(V>(
/solu } @jT-t]P
allsel ,7izrf8
asel,s,area,,2,3,1 !选择面积A2、A3 M2y"M ,k4
MPCHG,11,all !把所选择单元材料号改为11 WAp#[mW.fx
allsel ch]Q% M
solve X"
;ly0Mb
save,F2.8,db !把F=2.8时求解结果保存 \ B
D' "
!边坡在强度折减系数F=3.0时求解 hG.~[#[&6
finish uM9RlI 5
/solu oF#]<Z\
allsel 9h
0^_|"
asel,s,area,,2,3,1 !选择面积A2、A3 x 0x/2re
MPCHG,12,all !把所选择单元材料号改为12 ){XG%nC
allsel xnG,1doa
solve ,a?$F1Z-
save,F3.0,db !把F=3.0时求解结果保存 yX;v
(5）后处理 @d=4C{g%o
/post1 ！进入后处理 6"U )d7^
!边坡在强度折减系数F=1时结果分析 4) ~ GHb
Resume,'F1','db' !读入边坡在强度折减系数F=1时 b7? 2Pu
set,1,last !读入后一个子步 U^{'"x+
pldisp,1 !绘制边坡模型变形图 :1\QM'O
plnsol,u,x !绘制边坡模型水平方向位移云图 6GtXM3qtS
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 pZe:U; bb
!边坡在强度折减系数F=1.2时结果分析 4).>b3OhX
Resume,'F1.2','db' !读入边坡在强度折减系数F=1.2时 ~85Pgb<
set,1,last !读入后一个子步 W]DZ '
pldisp,1 !绘制边坡模型变形图 PEX(*GS
plnsol,u,x !绘制边坡模型水平方向位移云图 bol #[_~
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 Q d$!?h
!边坡在强度折减系数F=1.4时结果分析 hh>mX6A
Resume,'F1.4','db' !读入边坡在强度折减系数F=1.4时
mH*6Q>
set,1,last !读入后一个子步 _<1uO=km6
pldisp,1 !绘制边坡模型变形图 @AU<'?k
plnsol,u,x !绘制边坡模型水平方向位移云图 O8% Y .SK
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 A1(=7ZKz
!边坡在强度折减系数F=1.6时结果分析 T'nQj<dBt:
Resume,'F1.6','db' !读入边坡在强度折减系数F=1.6时 SmYY){AQ/
set,1,last !读入后一个子步 BKQI |i
pldisp,1 !绘制边坡模型变形图 zg83->[
plnsol,u,x !绘制边坡模型水平方向位移云图 _oJq32
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 z`[q$H7?
!边坡在强度折减系数F=1.8时结果分析 _P*<T6\J>
Resume,'F1.8','db' !读入边坡在强度折减系数F=1.8时 7I _lTu (
set,1,last !读入后一个子步 79 ZBVe(}
pldisp,1 !绘制边坡模型变形图 =D2x@ank[
plnsol,u,x !绘制边坡模型水平方向位移云图 X#*|_(^
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 s`:- 6{E
!边坡在强度折减系数F=2.0时结果分析 KW.QVBuVO#
Resume,'F2.0','db' !读入边坡在强度折减系数F=2.0时 +aWI"d--h
set,1,last !读入后一个子步 Ge]2g0
pldisp,1 !绘制边坡模型变形图 }$ZcC_
plnsol,u,x !绘制边坡模型水平方向位移云图 T7{Z0-
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 - 6KGQc}U
!边坡在强度折减系数F=2.2时结果分析 aGY F\7
Resume,'F2.2','db' !读入边坡在强度折减系数F=2.2时 L l
!J!{
set,1,last !读入后一个子步 BI,j/SRK
pldisp,1 !绘制边坡模型变形图 n(#159pZ
plnsol,u,x !绘制边坡模型水平方向位移云图 dM1)wkbET
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 @i[z4)"S
!边坡在强度折减系数F=2.4时结果分析 &[hq ! v
Resume,'F2.4','db' !读入边坡在强度折减系数F=2.4时 -S7rOq2Li
set,1,last !读入后一个子步 3 ./4] _p
pldisp,1 !绘制边坡模型变形图 `2@t) :
plnsol,u,x !绘制边坡模型水平方向位移云图 L27i_4 E,
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 .iB ?:
!边坡在强度折减系数F=2.6时结果分析 0=r.I}x
Resume,'F2.6','db' !读入边坡在强度折减系数F=2.6时 \e/'d~F
set,1,last !读入后一个子步 yjbqby7
pldisp,1 !绘制边坡模型变形图 `_6!nk q8
plnsol,u,x !绘制边坡模型水平方向位移云图 1cyX9X
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 |? !Ew# w
!边坡在强度折减系数F=2.8时结果分析 Ta38/v;S
Resume,'F2.8','db' !读入边坡在强度折减系数F=2.8时 U >z8gdzu
set,1,last !读入后一个子步 q.OkZI0n
pldisp,1 !绘制边坡模型变形图 I
=t{ u;
plnsol,u,x !绘制边坡模型水平方向位移云图 ': f q
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 y#3mc#)k
!边坡在强度折减系数F=3.0时结果分析 5CxD ys&<
Resume,'F3.0','db' !读入边坡在强度折减系数F=3.0时 &r6VF/
set,1,last !读入后一个子步 b0 Ct Qe
pldisp,1 !绘制边坡模型变形图 uOqWMRsoi
plnsol,u,x !绘制边坡模型水平方向位移云图 rZDlPp>BPZ
plnsol,eppl,eqv !绘制边坡模型塑性应变云图
