各位朋友： Z5Ihc%J^
大家好！我通过ansys用有限元强度折减法分析了一个边坡的稳定性，结果不同的工况下得出的塑性区是一样的，我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激！附命令流和塑性区如下：(1)创建物理环境
_)E8XyzF
/COM, Structural ! 指定结构分析 DH9?2)aR
/TITLE, Slope stability Analysis ! 定义工作标题 vB5iG|b}
/FILNAM,support,1 ! 定义工作文件名 ~Ls I<z
(2)建立模型 +&,\ J9'B
!进入前处理器 9Nu#&_2R
/PREP7 PAw g&._K
!定义单元类型 |V\.[F2Fe
ET,1,PLANE82 !定义边坡围岩单元 z8|9WZ:
KEYOPT,1,3,2 !求解类型为平面应变 rO$>zdmYHs
!定义材料属性 5 "am>$rh
!边坡围岩1材料属性 va(9{AXI
MP,EX,1,3E10 ZPsY0IzLo
MP,PRXY,1,0.25 [\9(@Bx
MP,DENS,1,2500 !定义塑性模型 ?0NSjK5ma
TB,DP,1 !定义数据表1 p qN[G=0
TBDATA,1,9E5,42.8 !给数据表输入数据 Ro]IE|Fv
!边坡围岩2材料属性 !采用弹性模型 uS#Cb+*F
MP,EX,2,3.2E10 ！定义弹性模量 %`QsX {?,
MP,PRXY,2,0.24 !定义泊松比 KN:dm!A
MP,DENS,2,2700 !定义密度 0juIkN#
!取强度折减系数F=1.2时边坡围岩材料属性 )_e"N d4
MP,EX,3,3E10 )m8>w6"
MP,PRXY,3,0.25 `^-Be
MP,DENS,3,2500 /Q8A"'Nk
TB,DP,3 !定义数据表3 T DIOK
TBDATA,1,7.5E5,37.7 !给数据表输入数据 1K9?a;.
!取强度折减系数F=1.4时边坡围岩材料属性 W<91m*
MP,EX,4,3E10 [ |n-x3h
MP,PRXY,4,0.25 !定义泊松比 &PuJV + y
MP,DENS,4,2500 -Y =c g;
TB,DP,4 !定义数据表4 M5kHD]b
TBDATA,1,6.4E5,33.5 !给数据表输入数据 d:pm|C|F
!取强度折减系数F=1.6时边坡围岩材料属性 ^3|$wB=
MP,EX,5,3E10 1vs>2` DLa
MP,PRXY,5,0.25 bM^A9BxD
MP,DENS,5,2500 ！定义密度 W lQ=CRY
TB,DP,5 !定义数据表5 s66XdM
TBDATA,1,5.6E5,30 !给数据表输入数据 !* Z)[[
!取强度折减系数F=1.8时边坡围岩材料属性 6r=)V$K <
MP,EX,6,3E10 ！定义弹性模量 jL+}F /~r
MP,PRXY,6,0.25 pE/3-0;}N
MP,DENS,6,2500 'uAC oME@
TB,DP,6 !定义数据表6 hS{
*l9v7
TBDATA,1,5.0E5,27.2 hav?mnVJ
!取强度折减系数F=2.0时边坡围岩材料属性 eBTedSM?t
MP,EX,7,3E10 N#['fg'
MP,PRXY,7,0.25 !定义泊松比 7(8
MP,DENS,7,2500 sLa)~To
TB,DP,7 !定义数据表7 zXPJ;^Xxa
TBDATA,1,4.5E5,24.8 *rz(}(r
!取强度折减系数F=2.2时边坡围岩材料属性 !VX_'GyK
MP,EX,8,3E10 q>(u>z!
MP,PRXY,8,0.25 G=!bM(]R~
MP,DENS,8,2500 oHXW])[
TB,DP,8 !定义数据表8 \G=R hx f
TBDATA,1,4.09E5,22.8 yr9A0F0
!取强度折减系数F=2.4时边坡围岩材料属性 o>;0NF| }
MP,EX,9,3E10 |C6(0fgWd
MP,PRXY,9,0.25 sQAc"S
MP,DENS,9,2500 ICbdKgLz
TB,DP,9 !定义数据表9 ,XG|oo -
TBDATA,1,3.6E5,21.1 !给数据表输入数据 Zmbz-##HQ
!取强度折减系数F=2.6时边坡围岩材料属性 M(zY[O
MP,EX,10,3E10 ！定义弹性模量 $ t# ,'M
MP,PRXY,10,0.25 !定义泊松比 qb> r\bc
MP,DENS,10,2500 ！定义密度 XjZao<?u
TB,DP,10 !定义数据表10 }0*ra37z>
TBDATA,1,3.46E5,19.6 !给数据表10输入数据 kqigFcz!Y
!取强度折减系数F=2.8时边坡围岩材料属性 sq(Ar(L<
MP,EX,11,3E10 ！定义弹性模量 &@utAuI
MP,PRXY,11,0.25 !定义泊松比 }t
D!xI;
MP,DENS,11,2500 X,EYa>RSy_
TB,DP,11 !定义数据表11 liw 9:@+V
TBDATA,1,3.2E5,18.2 !给数据表11输入数据 csX*XiDWm
!取强度折减系数F=3.0时边坡围岩材料属性 +'j*WVE%5
MP,EX,12,3E10 ！定义弹性模量 gQd=0"MV
MP,PRXY,12,0.25 *<1x:PR
MP,DENS,12,2500 ！定义密度 8+|V!q
TB,DP,12 !定义数据表12 `V):V4!j),
TBDATA,1,3.0E5,17.1 !给数据表12输入数据 p5;,/
|Ft
!建立几何模型 `J#xyDL6?
!创建边坡线模型线 ws2 j:B
K,1,,,, ！ 创建关键点1 l[ ": tG
K,2,-800,,, ENXW#{N.v
K,3,-800,-800,, R^M (fC
K,4,-800,-1200,, 6a]f&={E
K,5,1200,-1200,, \1`DaQp7
K,6,1200,-800,, oB06{/6
K,7,1200,,, W/r?0 E
K,8,1200,400,, *!]Epb
K,9,400,400,, |z|)r"*\4
l,1,2 !连接关键点1、2生成直线L1 199h Qxib:
l,2,3 @VC9gd O/
l,3,4 _2X6bIE
l,4,5 Qv0>Pf
l,5,6 P tQ#
l,6,7 @52=3
l,7,8 renmz,dJ,
l,8,9 iC|6roO!jk
l,9,1 Be>c)90bO_
l,1,7 QjjJtKz
l,3,6 5f&{ !N
!创建边坡面模型 y~c4:*L3
al,3,4,5,11 !由线L3、L4、L5、L11生成一个面A1 , HI%Xn
al,1,2,11,6,10 !由线L1、L2、L1、L6、L10生成一个面A2 0SCW2/o8
al,7,8,9,10 !由线L7、L8、L9、L10生成一个面A3 ym*#ZE`B!
/PNUM, area, on !打开面号开关 (z J$oRq
!划分网格生成有限元模型 0?Tk* X
! 划分边坡围岩2网格 o*wC{VP_
!设置网格份数 o% ^k T&
lsel,s,,,3,5,2 !选择线L3和L5 JeTrMa 2
lesize,all,,,5 !把所选择线分为5段 }Q r0T
lsel,s,,,4,11,7 !选择线L4和L11 Hrg=sR
lesize,all,,,26 !把所选择线分为26段 h[je _^5
mat,2 !给边坡围岩2赋予2号材料特性 -~ O;tJF2
type,1 !采用单元类型1 B,vHn2W
mshkey,0 !设定自由网格划分 D0y,TF
mshape,0 !设定四边形网格划分 JNM@Q
amesh,1 ！划分面积A1 `-K)K<
! 划分边坡围岩1网格 <aSLm=
!设置网格份数 /zG-\e U
lsel,s,,,2,6,4 !选择线L2和L6 _h=< _Z
lsel,a,,,1 !附加选择线L1 wdcryejCkr
lesize,all,,,10 !把所选择线分为10段 AV[P QI
lsel,s,,,8,10,2 !选择线L8和L10 h/0- Mrk;e
lesize,all,,,16 !把所选择线分为16段 E}b> 7L&w
lsel,s,,,7,9,2 lmtQr5U
lesize,all,,,12 W3{<e"
mat,1 !给边坡围岩1赋予1号材料特性 .`Old{<
type,1 iWN.3 |r
mshkey,0 !设定自由网格划分 qe6C|W~n
mshape,0 !设定四边形网格划分 ?^48Zq6wM
amesh,2 !划分面积A2 _
U8OIXN
amesh,3 ！划分面积A3 N7$DRG/<b
allsel X :SzkkVl7
!保存网格模型 Z_V&IQo-7
save,Slope-grid.db 18p3
(3)施加约束和荷载 o(X90X
!给边坡模型施加约束 U??f<
!边坡两侧施加X方向约束 v`+n`DT
nsel,s,loc,x,-800 !选择X=-800线上所有节点 4`!
nsel,a,loc,x,1200 !选择X=1200线上所有节点 F {*9[jY
d,all,ux !对所选择节点约束X方向位移 ~?`9i>3W~
allsel %<J(lC9,C
!边坡底部施加约束 ^W [B[Y<k
nsel,s,loc,y,-1200 !选择Y=--1200线上所有节点 K jn&
d,all,uy !对所选择节点约束X、Y方向位移 ghobu}wuF
d,all,ux Q0q)n=i }]
!施加重力加速度 vi! r8k
acel,,9.8 snTJe[^d
(4)求解 w] 5U
/solu IJ_ 'w[k
!求解设置 (m3I#L
antype,static !设定为静力求解 mam5 G!$
nsubst,100 !设定最大子步数为100 :S99}pgY
pred,on !打开时间步长预测器 *Nf4bH%MN
nropt,full !设定牛顿-拉普森选项 7Ysy\gZ&wp
nlgeom,on !打开大位移效果 4 &]To @>
lnsrch,on !打开线性搜索 "Yfr"1RmO
outres,all,all !输出所有项 X \p`pw$
cnvtol,f,,0.005,2,0.5 !力收敛准则设定 AYPf)K;%
cnvtol,u,,0.05,2,1 !位移收敛准则设定 3
!> L?
!边坡在强度折减系数F=1时求解 @[>+Dzn[6
allsel 'tV"^K QHI
solve !进行求解 uU[[[LQq
save,F1,db !把F=1时求解结果保存 d JQ }{,+6
!边坡在强度折减系数F=1.2时求解 z@%/r~?|
finish kf^-m/
/solu 0 bPJEEd
allsel ^= izqh5S
asel,s,area,,2,3,1 !选择面积A2、A3 Wr HgF*[
MPCHG,3,all ! 把所选择单元材料号改为3 3<)@ll
allsel [Z5}2gB&
solve !进行求解 g"Q}h
save,F1.2,db !把F=1.2时求解结果保存 ;!0.Kk
4
!边坡在强度折减系数F=1.4时求解 3h[:0W!C]
finish g=oeS%>E
/solu 'x45E.wYw
allsel 76IALJ00V
asel,s,area,,2,3,1 !选择面积A2、A3 {GX
&)c4
MPCHG,4,all !把所选择单元材料号改为4 yNqm]H3<MP
allsel ndKvJH 4
solve !进行求解 # McK46B z
save,F1.4,db !把F=1.4时求解结果保存 @u"kX2>Eq
!边坡在强度折减系数F=1.6时求解 (ju
aDn)
finish C?/r}ly<\
/solu q]iKz%|Z/
allsel SD|4ybK>d
asel,s,area,,2,3,1 !选择面积A2、A3 %KJhtd"q
MPCHG,5,all !把所选择单元材料号改为5 c5iormb"#
allsel w*6!?=jP
solve m.HX2(&\3
save,F1.6,db !把F=1.6时求解结果保存 ,p*ntj{
!边坡在强度折减系数F=1.8时求解 _!C)r*0(
finish N}}PlGp$
/solu ^Z-. [Y
allsel l NA' M&
asel,s,area,,2,3,1 !选择面积A2、A3 TM/ |K|_
MPCHG,6,all !把所选择单元材料号改为6 fV A=<:
allsel i B}LnC:
solve cFI7}#,5
save,F1.8,db !把F=1.8时求解结果保存 S4 k^&$;
!边坡在强度折减系数F=2.0时求解 ^`TKvcgIc
finish > G4HZE
/solu 3D$\y~HU
allsel 5}X<(q(
asel,s,area,,2,3,1 !选择面积A2、A3 !4 4mT'Y
MPCHG,7,all !把所选择单元材料号改为7 anz9lGG#
allsel vSYun I
solve "S B%02
save,F2.0,db !把F=2.0时求解结果保存 RP`GG+K
!边坡在强度折减系数F=2.2时求解 DX& lBV
finish hk"9D<&i>b
/solu zO).<xIq+
allsel a_ 9 |xI
asel,s,area,,2,3,1 !选择面积A2、A3 l?@MUsg+
MPCHG,8,all !把所选择单元材料号改为8 6_9:Eb=^v!
allsel "
g0-u(Y
solve 6cQeL$,SQ
save,F2.2,db !把F=2.2时求解结果保存 O{") i;v @
!边坡在强度折减系数F=2.4时求解 oA!5dpNhU
finish y?Hj %,
/solu -
5o<Q'(
allsel ~.4W,QLuD
asel,s,area,,2,3,1 !选择面积A2、A3 k}I5x1>&
MPCHG,9,all !把所选择单元材料号改为9 u"#6_-0 y
allsel j h1 bn
solve o&hKg#nO83
save,F2.4,db !把F=2.4时求解结果保存 Y @XkqvX
!边坡在强度折减系数F=2.6时求解 *3.yumcv{L
finish B{OW}D$P#
/solu B~V<n&<
allsel _>jrlIfc
asel,s,area,,2,3,1 !选择面积A2、A3 75\RG+kQ
MPCHG,10,all !把所选择单元材料号改为10 ;9p#xW 6
allsel ~ F?G5cN5
solve =q" w2b&
save,F2.6,db !把F=2.6时求解结果保存 t-eKr uj+
!边坡在强度折减系数F=2.8时求解 KD%xo/Z.
finish _#J_$CE#
/solu EU^}NZW&v:
allsel cYq' ]$]
asel,s,area,,2,3,1 !选择面积A2、A3 cwM#X;FGq
MPCHG,11,all !把所选择单元材料号改为11 >Bh)7>`3c
allsel !!-}ttFA
solve +
4V1>e+
save,F2.8,db !把F=2.8时求解结果保存 X ]pR,\B
!边坡在强度折减系数F=3.0时求解 Z0( ) pT
finish ) 8x:x7 ?
/solu ;"d ,~nLn
allsel .y %pGi
asel,s,area,,2,3,1 !选择面积A2、A3 @pqY9_:P1
MPCHG,12,all !把所选择单元材料号改为12 XQ0#0<
allsel J+3\2D?
solve u5cV z_S
save,F3.0,db !把F=3.0时求解结果保存 vB{; N
(5）后处理 To# E@Nw
/post1 ！进入后处理 .-('C> @
!边坡在强度折减系数F=1时结果分析 qEK4I}Q-=
Resume,'F1','db' !读入边坡在强度折减系数F=1时 k7yv>iN
set,1,last !读入后一个子步 /`4v"f0V
pldisp,1 !绘制边坡模型变形图 y"|K
|QT
plnsol,u,x !绘制边坡模型水平方向位移云图 r&%gjqt
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 t`<}UWAH+
!边坡在强度折减系数F=1.2时结果分析 >} dTO/
Resume,'F1.2','db' !读入边坡在强度折减系数F=1.2时 C}(<PN T
set,1,last !读入后一个子步 ]HJ{dcF
pldisp,1 !绘制边坡模型变形图 #51 4a(6
plnsol,u,x !绘制边坡模型水平方向位移云图 vDK:v$g
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 pIZLGsu[
!边坡在强度折减系数F=1.4时结果分析 v2M"b?Q
Resume,'F1.4','db' !读入边坡在强度折减系数F=1.4时 r6F{
set,1,last !读入后一个子步 u_}`y1Xu#
pldisp,1 !绘制边坡模型变形图 e?^ \r)1
plnsol,u,x !绘制边坡模型水平方向位移云图 S.Wh4kMUe
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 3r~>~ueZ
!边坡在强度折减系数F=1.6时结果分析 P S _3Oq)
Resume,'F1.6','db' !读入边坡在强度折减系数F=1.6时 PmPyb>HK=P
set,1,last !读入后一个子步 gtaV6sD
pldisp,1 !绘制边坡模型变形图 b!4N)t>gl
plnsol,u,x !绘制边坡模型水平方向位移云图 Qm35{^p+
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 ;PfeP ;z
!边坡在强度折减系数F=1.8时结果分析 zqHpT^ B?
Resume,'F1.8','db' !读入边坡在强度折减系数F=1.8时 R
"/xne
set,1,last !读入后一个子步 pIID= 8RJ.
pldisp,1 !绘制边坡模型变形图 n=<NFkeX
plnsol,u,x !绘制边坡模型水平方向位移云图 Wz6]*P`qv
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 |dl0B26x
!边坡在强度折减系数F=2.0时结果分析 xe cieC
Resume,'F2.0','db' !读入边坡在强度折减系数F=2.0时 "t (1tWO1o
set,1,last !读入后一个子步 ||{T5E-.F
pldisp,1 !绘制边坡模型变形图 ! F0r d9
plnsol,u,x !绘制边坡模型水平方向位移云图 5YTb7M
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 RsqRR`|X?
!边坡在强度折减系数F=2.2时结果分析 *}
*!+C3
Resume,'F2.2','db' !读入边坡在强度折减系数F=2.2时 )XmCy"xx
set,1,last !读入后一个子步 c W81
pldisp,1 !绘制边坡模型变形图 AkYupP2]v
plnsol,u,x !绘制边坡模型水平方向位移云图 _" ?c9
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 iXUWIgr
!边坡在强度折减系数F=2.4时结果分析 >zPO>.?h7T
Resume,'F2.4','db' !读入边坡在强度折减系数F=2.4时 ^f^-.X
set,1,last !读入后一个子步
K;<NBnH
pldisp,1 !绘制边坡模型变形图 Xk {!' 0
plnsol,u,x !绘制边坡模型水平方向位移云图 P[Y{LKAbb
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 D rF
!边坡在强度折减系数F=2.6时结果分析 0%;N9\
Resume,'F2.6','db' !读入边坡在强度折减系数F=2.6时 X[GIOPDx
set,1,last !读入后一个子步 rQu
pldisp,1 !绘制边坡模型变形图 VZT6;1TD$8
plnsol,u,x !绘制边坡模型水平方向位移云图 +Fc ET
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 #Acon7R p
!边坡在强度折减系数F=2.8时结果分析 (*Gi~?-
Resume,'F2.8','db' !读入边坡在强度折减系数F=2.8时 )CC? vV
set,1,last !读入后一个子步 ysSjc
pldisp,1 !绘制边坡模型变形图 HZ QDe&
plnsol,u,x !绘制边坡模型水平方向位移云图 38V $ <w
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 Hk <X
!边坡在强度折减系数F=3.0时结果分析 4c5^7";P
Resume,'F3.0','db' !读入边坡在强度折减系数F=3.0时 _L&n&y1+%
set,1,last !读入后一个子步 [+hy_Nc$
pldisp,1 !绘制边坡模型变形图 8N=%X-R%
plnsol,u,x !绘制边坡模型水平方向位移云图 Y=P9:unG
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 b{~fVil$y
Mv/IMO0rR
%+AS0 JhB
9 0
7>nhIp))
1KeJd&e
+8LM~voB
