各位朋友： m`q&[:
大家好！我通过ansys用有限元强度折减法分析了一个边坡的稳定性，结果不同的工况下得出的塑性区是一样的，我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激！附命令流和塑性区如下：(1)创建物理环境 ew dTsgt'
/COM, Structural ! 指定结构分析 gd qBT]j
/TITLE, Slope stability Analysis ! 定义工作标题 .b<W*4{j0H
/FILNAM,support,1 ! 定义工作文件名 ]yqE6Lf9
(2)建立模型 :wg =H
!进入前处理器 ^=5y;
/PREP7 <,/ 7:n
!定义单元类型 : WQlpL n
ET,1,PLANE82 !定义边坡围岩单元 z6d0Y $A G
KEYOPT,1,3,2 !求解类型为平面应变 1t^9.!$@y
!定义材料属性 %3t;[$n#
!边坡围岩1材料属性 4 J(-~
MP,EX,1,3E10 ln8NcAEx
MP,PRXY,1,0.25 A! HJ
MP,DENS,1,2500 !定义塑性模型 P*|=Z>%[0
TB,DP,1 !定义数据表1 Kj3Gm>B<y
TBDATA,1,9E5,42.8 !给数据表输入数据
w >2G@
!边坡围岩2材料属性 !采用弹性模型 Ac|dmu
MP,EX,2,3.2E10 ！定义弹性模量 I"3C/ pU2
MP,PRXY,2,0.24 !定义泊松比 W+"^! p|
MP,DENS,2,2700 !定义密度 6H U*,
!取强度折减系数F=1.2时边坡围岩材料属性 0MxK+8\y
MP,EX,3,3E10 Lr8|S
MP,PRXY,3,0.25 SVd@-
'-K
MP,DENS,3,2500 (>x05nh
TB,DP,3 !定义数据表3 >35w"a7S
TBDATA,1,7.5E5,37.7 !给数据表输入数据 :KXI@)M
!取强度折减系数F=1.4时边坡围岩材料属性 OUGkam0UK
MP,EX,4,3E10 ~$[fG}C.K
MP,PRXY,4,0.25 !定义泊松比 ;]> )6
MP,DENS,4,2500 q^zG+FN
TB,DP,4 !定义数据表4 8c9<kGm$E
TBDATA,1,6.4E5,33.5 !给数据表输入数据 -D=Sj@G
!取强度折减系数F=1.6时边坡围岩材料属性 aL9 0:,V
MP,EX,5,3E10 kRX?o'U~C
MP,PRXY,5,0.25 M,li\) J!&
MP,DENS,5,2500 ！定义密度 {~Jk (c~I
TB,DP,5 !定义数据表5 )YAU|sCAi$
TBDATA,1,5.6E5,30 !给数据表输入数据 8{i}^.p
!取强度折减系数F=1.8时边坡围岩材料属性 h2Th)&Fb>
MP,EX,6,3E10 ！定义弹性模量 @0t[7Nv-1
MP,PRXY,6,0.25 &^HVuYa.0
MP,DENS,6,2500 $)9 |"q6
TB,DP,6 !定义数据表6 C^B $_?
TBDATA,1,5.0E5,27.2 %
C2Vga#
!取强度折减系数F=2.0时边坡围岩材料属性 +0Q +0:
MP,EX,7,3E10 NR
k~
MP,PRXY,7,0.25 !定义泊松比 > 3(,s^
MP,DENS,7,2500 Nz`v+sp
TB,DP,7 !定义数据表7 $^?VyHXvY
TBDATA,1,4.5E5,24.8 oZ t Cx
!取强度折减系数F=2.2时边坡围岩材料属性 |Z2"pV
MP,EX,8,3E10 whHuV*K}
MP,PRXY,8,0.25 #Cu$y8~as
MP,DENS,8,2500 z;<~j=lP
TB,DP,8 !定义数据表8 q%$p56\?3
TBDATA,1,4.09E5,22.8 &Q}%b7
!取强度折减系数F=2.4时边坡围岩材料属性 E7@Gpu,o
MP,EX,9,3E10 PO6yE r
MP,PRXY,9,0.25 ~UO} PI`C
MP,DENS,9,2500 u1#(~[.
TB,DP,9 !定义数据表9 :@-yK8q's
TBDATA,1,3.6E5,21.1 !给数据表输入数据 ?(K=du
!取强度折减系数F=2.6时边坡围岩材料属性 m$v >r\*X
MP,EX,10,3E10 ！定义弹性模量 y6[ le*T
MP,PRXY,10,0.25 !定义泊松比 \>lA2^E f
MP,DENS,10,2500 ！定义密度 i +~BVb
TB,DP,10 !定义数据表10 =l*xM/S
TBDATA,1,3.46E5,19.6 !给数据表10输入数据 2?Jw0Wq5D
!取强度折减系数F=2.8时边坡围岩材料属性 +s8R]3NJ_H
MP,EX,11,3E10 ！定义弹性模量 .S/zxf~h
MP,PRXY,11,0.25 !定义泊松比 Xfqin4/jC
MP,DENS,11,2500 rrj.]^E_~
TB,DP,11 !定义数据表11 3^ y<Db
TBDATA,1,3.2E5,18.2 !给数据表11输入数据 m }RZ )c
!取强度折减系数F=3.0时边坡围岩材料属性 "gm5 DE
MP,EX,12,3E10 ！定义弹性模量 Z~-N'L t{
MP,PRXY,12,0.25 m9:ah<
MP,DENS,12,2500 ！定义密度 Y(kf<Wo
TB,DP,12 !定义数据表12 SvvNk
TBDATA,1,3.0E5,17.1 !给数据表12输入数据 1%N*GJlwJ
!建立几何模型 w <"mS*Q
!创建边坡线模型线 'OP0#`6`
K,1,,,, ！ 创建关键点1 ?f q!BV
K,2,-800,,, 4Nt4(3Kf
K,3,-800,-800,, u|AMqS
K,4,-800,-1200,, es #6/
K,5,1200,-1200,, Zxql hq/)
K,6,1200,-800,, /w5*R5B{
K,7,1200,,, Dr%wa b"yy
K,8,1200,400,, Qb/:E}h]$
K,9,400,400,, %3#C0%{x
l,1,2 !连接关键点1、2生成直线L1 8uH8)
l,2,3 "Z,T%]
l,3,4 T=M##`jP%
l,4,5 vo*oCfm
l,5,6 CZeZk
l,6,7 zSfUM.fM
l,7,8 =4SXntU!e
l,8,9 ` W~
l,9,1 9609
l,1,7 R0tT4V +
l,3,6 x`@`y7(
!创建边坡面模型 ~ |A0*
al,3,4,5,11 !由线L3、L4、L5、L11生成一个面A1 $)o0{HsL+
al,1,2,11,6,10 !由线L1、L2、L1、L6、L10生成一个面A2 Xz)F-C27h
al,7,8,9,10 !由线L7、L8、L9、L10生成一个面A3 Mz2TwU_
/PNUM, area, on !打开面号开关 Kn@#5MC
rU
!划分网格生成有限元模型 JJbd h \
! 划分边坡围岩2网格 2=8PA/
!设置网格份数 5<v1v&
lsel,s,,,3,5,2 !选择线L3和L5 Q25VG5 G
lesize,all,,,5 !把所选择线分为5段 ^5TVm>F@3
lsel,s,,,4,11,7 !选择线L4和L11 F4z#u2 ~TC
lesize,all,,,26 !把所选择线分为26段 6<fG; :
mat,2 !给边坡围岩2赋予2号材料特性 Vym 0|cW
type,1 !采用单元类型1 MO7R3PP
mshkey,0 !设定自由网格划分 HZ Wt >f
mshape,0 !设定四边形网格划分 $m*Gu:#xm&
amesh,1 ！划分面积A1 D^ .
c:
! 划分边坡围岩1网格 GCO: !,1
!设置网格份数 a*.#Zgy:lK
lsel,s,,,2,6,4 !选择线L2和L6 URbB2
Bi
lsel,a,,,1 !附加选择线L1 7[qL~BT+
lesize,all,,,10 !把所选择线分为10段 Jx}-Y*
o
lsel,s,,,8,10,2 !选择线L8和L10 j_<!y(W
lesize,all,,,16 !把所选择线分为16段 $^{#hYq)o
lsel,s,,,7,9,2 ~:f..|JM
lesize,all,,,12 ]|,}hsN
mat,1 !给边坡围岩1赋予1号材料特性 R"P-+T=7M
type,1 Lkx~>U
mshkey,0 !设定自由网格划分 R*lq7n9
mshape,0 !设定四边形网格划分 )&>W/56/
amesh,2 !划分面积A2 nC% qdzT
amesh,3 ！划分面积A3 YMK ![ q-
allsel C<(oaeQY
!保存网格模型 mS![J69(
save,Slope-grid.db Fih
pp<
(3)施加约束和荷载 {xov8 M
!给边坡模型施加约束 0 \ (:y^X
!边坡两侧施加X方向约束 3Xd:LDZ{
nsel,s,loc,x,-800 !选择X=-800线上所有节点 E JuTv%Y8
nsel,a,loc,x,1200 !选择X=1200线上所有节点 3Z*o5@RI
d,all,ux !对所选择节点约束X方向位移 <y^_&9
allsel =IQ}Y_xr
!边坡底部施加约束 @/^mFqr2
nsel,s,loc,y,-1200 !选择Y=--1200线上所有节点 BYM6cp+S
d,all,uy !对所选择节点约束X、Y方向位移 sHk>ek]2I
d,all,ux {9V.l.Q
!施加重力加速度 P3|s}&
acel,,9.8 O]@#53)Tz
(4)求解 h
k a_Fo
/solu -6NoEmb)\'
!求解设置 a <?~1pWtc
antype,static !设定为静力求解 ZM v\j|{8
nsubst,100 !设定最大子步数为100 &b5(Su
pred,on !打开时间步长预测器 vVa|E#
[
nropt,full !设定牛顿-拉普森选项 0 ^o/c SF
nlgeom,on !打开大位移效果 ZnI15bsDx
lnsrch,on !打开线性搜索 jED.0,+K !
outres,all,all !输出所有项 id5 ` YA$
cnvtol,f,,0.005,2,0.5 !力收敛准则设定 y||RK` H
cnvtol,u,,0.05,2,1 !位移收敛准则设定 gz[3 xH~
!边坡在强度折减系数F=1时求解 _Q
I!UQdW
allsel 81s
} 4
solve !进行求解 *. |%uf.
save,F1,db !把F=1时求解结果保存 YT(Eh3ID
!边坡在强度折减系数F=1.2时求解 ;xzaW4(3
finish C]5 kQ1Og
/solu [
fzYC'A=
allsel kV?fie<\)
asel,s,area,,2,3,1 !选择面积A2、A3 bl^Ihza
MPCHG,3,all ! 把所选择单元材料号改为3 Bz-jy.
allsel . yXqa"p
solve !进行求解 $>O~7Nfst7
save,F1.2,db !把F=1.2时求解结果保存 F/> \ uzu
!边坡在强度折减系数F=1.4时求解 !R\FCAW[ x
finish Y01! D"{\
/solu lbIPtu
allsel e]88 4FP
asel,s,area,,2,3,1 !选择面积A2、A3 XJ3sqc S
MPCHG,4,all !把所选择单元材料号改为4 o#f"wQH;p
allsel \#dacQ2E@
solve !进行求解 pUqC 88*j
save,F1.4,db !把F=1.4时求解结果保存 jLVD37 P^
!边坡在强度折减系数F=1.6时求解 _r\M }lDh*
finish =%I yR
/solu QNU~G3
allsel 6Nn+7z<*&z
asel,s,area,,2,3,1 !选择面积A2、A3 fpo{`;&F
MPCHG,5,all !把所选择单元材料号改为5 =VuSi(d;e{
allsel 7(.Z8AO
solve p5or"tK
save,F1.6,db !把F=1.6时求解结果保存 3lpx h_
!边坡在强度折减系数F=1.8时求解 M;AD L|
finish x@ rQ7K >
/solu 2m yxwA5
allsel -2J37
asel,s,area,,2,3,1 !选择面积A2、A3 `DgK$ QM
MPCHG,6,all !把所选择单元材料号改为6 0g|5s
allsel ~ B JE~
solve 9A B U^ig
save,F1.8,db !把F=1.8时求解结果保存
-NN=(p!<
!边坡在强度折减系数F=2.0时求解 HV/:OCK
finish (iir,Ks2C
/solu ^OWG9`p+
allsel k"& o)*d
asel,s,area,,2,3,1 !选择面积A2、A3 wxh\CBxG
MPCHG,7,all !把所选择单元材料号改为7 V27RK-.N!
allsel QtKcv7:4
solve S}%z0 g<
save,F2.0,db !把F=2.0时求解结果保存 Ta?}n^V?;
!边坡在强度折减系数F=2.2时求解 Wmcd{MOS
finish N2A6 C$s
/solu EC,`t*<
allsel j`RG Moq
asel,s,area,,2,3,1 !选择面积A2、A3 *1`X}
MPCHG,8,all !把所选择单元材料号改为8 Z8xB
a0
allsel b1 w@toc
solve }p2 iF2g9`
save,F2.2,db !把F=2.2时求解结果保存 1s =Q~*f~d
!边坡在强度折减系数F=2.4时求解 Gg9MAK\ C9
finish ~ d]v{<3
/solu wT" :
allsel SU~.baP?
asel,s,area,,2,3,1 !选择面积A2、A3 a!: N
C
MPCHG,9,all !把所选择单元材料号改为9 |JYb4J4Ni
allsel V)/J2 -w
solve LiT%d
save,F2.4,db !把F=2.4时求解结果保存 Jk\-e` eE
!边坡在强度折减系数F=2.6时求解 A2M(
ad
finish #d\&6'O
/solu =#W:z.w
allsel S5 q1M n
asel,s,area,,2,3,1 !选择面积A2、A3 ..u{v}4&
MPCHG,10,all !把所选择单元材料号改为10 {E Ay~lo
allsel 9_:"`)] 3B
solve H2R3I<j
save,F2.6,db !把F=2.6时求解结果保存 7mMGH(
!边坡在强度折减系数F=2.8时求解 ;Hk3y+&]a
finish "*t6KXVaM
/solu (wZ!OLY%}
allsel uWdF7|PN7
asel,s,area,,2,3,1 !选择面积A2、A3 P%e7 c,
MPCHG,11,all !把所选择单元材料号改为11 04|ZwX$>+
allsel = N* Jis
solve <.4(#Ebd
save,F2.8,db !把F=2.8时求解结果保存 *
CR#D}F
!边坡在强度折减系数F=3.0时求解 s~ 8 g
finish PZ8,E{V
/solu 2Wluc37
allsel LPt9+sauf1
asel,s,area,,2,3,1 !选择面积A2、A3 ,k4pW&A
MPCHG,12,all !把所选择单元材料号改为12 oH x :["F
allsel oxc; DfJ_
solve % mT/y%&:
save,F3.0,db !把F=3.0时求解结果保存 PJN9[Y{^3
(5）后处理 <L qJg
/post1 ！进入后处理 n
A b~
!边坡在强度折减系数F=1时结果分析 BK%B[f*[OA
Resume,'F1','db' !读入边坡在强度折减系数F=1时 ?}s;, _GH
set,1,last !读入后一个子步 $]E+E .P
pldisp,1 !绘制边坡模型变形图 \/3(>g?4
plnsol,u,x !绘制边坡模型水平方向位移云图 g[pU5%|" [
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 nI6om pTX
!边坡在强度折减系数F=1.2时结果分析 \vT~2Y(K
Resume,'F1.2','db' !读入边坡在强度折减系数F=1.2时 9Tt%~m^
set,1,last !读入后一个子步 z&d.YO_W
pldisp,1 !绘制边坡模型变形图 pK3A/ry<
plnsol,u,x !绘制边坡模型水平方向位移云图 ktTP~7UVi
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 66eJp-5e8
!边坡在强度折减系数F=1.4时结果分析 aHW34e@ebL
Resume,'F1.4','db' !读入边坡在强度折减系数F=1.4时 K}@rte
set,1,last !读入后一个子步 $Xlr@)%
pldisp,1 !绘制边坡模型变形图 r]p3DQ
plnsol,u,x !绘制边坡模型水平方向位移云图 !X-\;3kC0
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 VM\R- [
!边坡在强度折减系数F=1.6时结果分析 C'$}{%Cc@$
Resume,'F1.6','db' !读入边坡在强度折减系数F=1.6时 " E2 0Y"[h
set,1,last !读入后一个子步 aH _c 84DS
pldisp,1 !绘制边坡模型变形图 P5_Ajb(@'
plnsol,u,x !绘制边坡模型水平方向位移云图 lY
tt|J
plnsol,eppl,eqv !绘制边坡模型塑性应变云图
{ % X2K
!边坡在强度折减系数F=1.8时结果分析 uM[|>t
Resume,'F1.8','db' !读入边坡在强度折减系数F=1.8时 lF! PiL
set,1,last !读入后一个子步 tp cB}HUv
pldisp,1 !绘制边坡模型变形图 lF}@ @e)N
plnsol,u,x !绘制边坡模型水平方向位移云图 J Ah!#S(
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 @L!^2 v
!边坡在强度折减系数F=2.0时结果分析 ;Y*K!iFWH
Resume,'F2.0','db' !读入边坡在强度折减系数F=2.0时 `~u=[} w
set,1,last !读入后一个子步 iXn XZ|M
pldisp,1 !绘制边坡模型变形图 cHF W"g78
plnsol,u,x !绘制边坡模型水平方向位移云图 ftPps -
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 LsERcjwwK
!边坡在强度折减系数F=2.2时结果分析 I&La0g _E
Resume,'F2.2','db' !读入边坡在强度折减系数F=2.2时 ^ l]!'"
set,1,last !读入后一个子步 ekyCZ8iai
pldisp,1 !绘制边坡模型变形图 G:$kGzhJ
plnsol,u,x !绘制边坡模型水平方向位移云图 {\1bWr8!U
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 15j5F5P
!边坡在强度折减系数F=2.4时结果分析 C
6
\
Resume,'F2.4','db' !读入边坡在强度折减系数F=2.4时 &]n }fq
set,1,last !读入后一个子步 C][hH?.
pldisp,1 !绘制边坡模型变形图 ,6g{-r-2
plnsol,u,x !绘制边坡模型水平方向位移云图 L4/ns@e
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 %[*-aA
!边坡在强度折减系数F=2.6时结果分析 'M"z3j]m-,
Resume,'F2.6','db' !读入边坡在强度折减系数F=2.6时 0@ zJa;z'
set,1,last !读入后一个子步 St%x \[D
pldisp,1 !绘制边坡模型变形图 a4^hC[a
plnsol,u,x !绘制边坡模型水平方向位移云图 +-|""`I 1I
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 [6mK<A,/
!边坡在强度折减系数F=2.8时结果分析 ,#ZPg_x?1
Resume,'F2.8','db' !读入边坡在强度折减系数F=2.8时 ru eaP
set,1,last !读入后一个子步 q\o#<'F1J
pldisp,1 !绘制边坡模型变形图 <7J3tn B
plnsol,u,x !绘制边坡模型水平方向位移云图 /OztkThx=
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 2w7$"N
!边坡在强度折减系数F=3.0时结果分析 x7zc3%T's
Resume,'F3.0','db' !读入边坡在强度折减系数F=3.0时 3O$l;|SX
set,1,last !读入后一个子步 ]z^jz#>um&
pldisp,1 !绘制边坡模型变形图 ?; W"=I*3
plnsol,u,x !绘制边坡模型水平方向位移云图 cl^UFl f[
plnsol,eppl,eqv !绘制边坡模型塑性应变云图
