各位朋友： =}PdH`S
大家好！我通过ansys用有限元强度折减法分析了一个边坡的稳定性，结果不同的工况下得出的塑性区是一样的，我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激！附命令流和塑性区如下：(1)创建物理环境 L. ]$6Q0
/COM, Structural ! 指定结构分析 z"*/mP2
/TITLE, Slope stability Analysis ! 定义工作标题 XT ;u<aJs
/FILNAM,support,1 ! 定义工作文件名 ^( Rvk
(2)建立模型 =V"(AuCVE
!进入前处理器 ]0L&v7[
/PREP7 t'm;:J1
!定义单元类型 b=3H
ET,1,PLANE82 !定义边坡围岩单元 Gn;@{x6
KEYOPT,1,3,2 !求解类型为平面应变 _,</1~.
!定义材料属性 Dde]I_f}
!边坡围岩1材料属性 nNXgW
MP,EX,1,3E10 M4xi1M#%
MP,PRXY,1,0.25 D /h/Y) Y
MP,DENS,1,2500 !定义塑性模型 M`6y@<
TB,DP,1 !定义数据表1 Jjl`_X$CB
TBDATA,1,9E5,42.8 !给数据表输入数据 h5 yzwj:C?
!边坡围岩2材料属性 !采用弹性模型 "(qw-kil
MP,EX,2,3.2E10 ！定义弹性模量 vky . ^
MP,PRXY,2,0.24 !定义泊松比 m[k_>e\ u
MP,DENS,2,2700 !定义密度 uIU5.\"s
!取强度折减系数F=1.2时边坡围岩材料属性 b* ef);
MP,EX,3,3E10 ki>~H!zB
MP,PRXY,3,0.25 ':R,53tjl
MP,DENS,3,2500 x nE |Umz
TB,DP,3 !定义数据表3 %p
X6QRt?
TBDATA,1,7.5E5,37.7 !给数据表输入数据 HNL42\Kz!
!取强度折减系数F=1.4时边坡围岩材料属性 gNG r!3*)w
MP,EX,4,3E10 cWajrLw
MP,PRXY,4,0.25 !定义泊松比 g R
nOd
MP,DENS,4,2500 1,5E `J
TB,DP,4 !定义数据表4 kp\\"+,VC
TBDATA,1,6.4E5,33.5 !给数据表输入数据 h=_mNG>R)
!取强度折减系数F=1.6时边坡围岩材料属性 t\$U`V)
MP,EX,5,3E10 Et_V,s<|
MP,PRXY,5,0.25 R-^96fFBy
MP,DENS,5,2500 ！定义密度 0| ;
.6\
TB,DP,5 !定义数据表5 fu$R7
TBDATA,1,5.6E5,30 !给数据表输入数据 K!,<7[MBg
!取强度折减系数F=1.8时边坡围岩材料属性 M@W[Bz
MP,EX,6,3E10 ！定义弹性模量 S .!UPkW H
MP,PRXY,6,0.25 _w*}\~`=^
MP,DENS,6,2500 :$+-3_oLMQ
TB,DP,6 !定义数据表6 ;\MW$/[JCy
TBDATA,1,5.0E5,27.2 @ |'5 n
!取强度折减系数F=2.0时边坡围岩材料属性 Hi]cxD*`
MP,EX,7,3E10 6 =N`wi
MP,PRXY,7,0.25 !定义泊松比 7)%+=@
MP,DENS,7,2500 :rP#I#,7w
TB,DP,7 !定义数据表7 67y Tvr@a
TBDATA,1,4.5E5,24.8 $H$j-)\D
!取强度折减系数F=2.2时边坡围岩材料属性 US
MP,EX,8,3E10 -|rLs$V1r
MP,PRXY,8,0.25 -[L!3jU
MP,DENS,8,2500 d,$d~alY
TB,DP,8 !定义数据表8 ;l$ \6T
TBDATA,1,4.09E5,22.8 ,.gQ^^+=
!取强度折减系数F=2.4时边坡围岩材料属性 TY(bPq
MP,EX,9,3E10 'EFyIVezg9
MP,PRXY,9,0.25 r]ShZBAbYp
MP,DENS,9,2500 JMdP wI
TB,DP,9 !定义数据表9 U.{l;EL:T
TBDATA,1,3.6E5,21.1 !给数据表输入数据 r <
cVp^
!取强度折减系数F=2.6时边坡围岩材料属性 \
u_ui
MP,EX,10,3E10 ！定义弹性模量 5{ $LsL
MP,PRXY,10,0.25 !定义泊松比 z#F.xVg'
MP,DENS,10,2500 ！定义密度 OxGE%R,
TB,DP,10 !定义数据表10 ,Z S6jZ
TBDATA,1,3.46E5,19.6 !给数据表10输入数据 GA`PY-Vs)
!取强度折减系数F=2.8时边坡围岩材料属性 !a$ D4(`v
MP,EX,11,3E10 ！定义弹性模量 e *j .
MP,PRXY,11,0.25 !定义泊松比 {-BRt)L[
MP,DENS,11,2500 ZtHm\VTS
TB,DP,11 !定义数据表11 f3|@|'
;
TBDATA,1,3.2E5,18.2 !给数据表11输入数据 Pl-5ncb\
!取强度折减系数F=3.0时边坡围岩材料属性 FYS/##r
MP,EX,12,3E10 ！定义弹性模量 )J?{+3
MP,PRXY,12,0.25 upvS|KUil
MP,DENS,12,2500 ！定义密度 @x c ', I
TB,DP,12 !定义数据表12 >&!RWH9 *q
TBDATA,1,3.0E5,17.1 !给数据表12输入数据 :R.&`4=X
!建立几何模型 vy,&N^P
!创建边坡线模型线 >3u ]OSb
K,1,,,, ！ 创建关键点1 w{k)XY40sW
K,2,-800,,, Dz./w
K,3,-800,-800,, dJ?XPo"Cm=
K,4,-800,-1200,, }h 3K@R
K,5,1200,-1200,,
y< C<_2
K,6,1200,-800,, }K hjlPhx
K,7,1200,,, k=GG>]<i
K,8,1200,400,, _ bD/D !|
K,9,400,400,, RmrL^asg
l,1,2 !连接关键点1、2生成直线L1 =+w*gDr
l,2,3 -)vEWn$3<
l,3,4 BnRN;bu
l,4,5 2YuN~-
l,5,6 NzKUtwnIz
l,6,7 n4lu tnF
l,7,8 Ej7 /X ~
l,8,9 |j3'eW& =
l,9,1 +y 87~]]
l,1,7 0j(M*
sl
l,3,6 WL+]4Wiz
!创建边坡面模型 hXGwP4
al,3,4,5,11 !由线L3、L4、L5、L11生成一个面A1 L#)(H^[
al,1,2,11,6,10 !由线L1、L2、L1、L6、L10生成一个面A2 /*Qq[C
al,7,8,9,10 !由线L7、L8、L9、L10生成一个面A3 e|4&b@
/PNUM, area, on !打开面号开关 XlI!{qj|
!划分网格生成有限元模型 *._|- L
! 划分边坡围岩2网格 7 h y&-<
!设置网格份数 Dup;e&9g
lsel,s,,,3,5,2 !选择线L3和L5
rxO2QQ%V
lesize,all,,,5 !把所选择线分为5段
b3 YO!cJ
lsel,s,,,4,11,7 !选择线L4和L11 _Jv
9F8v
lesize,all,,,26 !把所选择线分为26段 |y<),j6
mat,2 !给边坡围岩2赋予2号材料特性 &Z?ut *%S
type,1 !采用单元类型1 #.[AK_S5&
mshkey,0 !设定自由网格划分 )etmE
mshape,0 !设定四边形网格划分 8.bKb<y
amesh,1 ！划分面积A1 s( <u o{
! 划分边坡围岩1网格 f&D]anf33
!设置网格份数 |y20Hi':
lsel,s,,,2,6,4 !选择线L2和L6 8}w6z7e|{
lsel,a,,,1 !附加选择线L1 m5 G \}8|
lesize,all,,,10 !把所选择线分为10段 >yJ9U,Y
lsel,s,,,8,10,2 !选择线L8和L10 %zU`XVNN+
lesize,all,,,16 !把所选择线分为16段 dz>;<&2Z
lsel,s,,,7,9,2 =uDgzdDyE
lesize,all,,,12 * Ei|fe$sa
mat,1 !给边坡围岩1赋予1号材料特性 'B0{U4?
type,1 0q\7C[R_
mshkey,0 !设定自由网格划分 |w}xl '>q
mshape,0 !设定四边形网格划分 `"@ X.}\
amesh,2 !划分面积A2 _tr<}PnZ
amesh,3 ！划分面积A3 ;WL1B
allsel U}SXJH&&E
!保存网格模型 6WoAs)ZF
save,Slope-grid.db 'Pvm8t
(3)施加约束和荷载 7*DMVok:
!给边坡模型施加约束 - y9>;6
!边坡两侧施加X方向约束 5X.e*;
nsel,s,loc,x,-800 !选择X=-800线上所有节点 n}xhW'3hU=
nsel,a,loc,x,1200 !选择X=1200线上所有节点 fJZp?e"
d,all,ux !对所选择节点约束X方向位移 {G*A.$-d
allsel E2L(wt}^
!边坡底部施加约束 |# l=
nsel,s,loc,y,-1200 !选择Y=--1200线上所有节点 q2:K 4
d,all,uy !对所选择节点约束X、Y方向位移 Z>)][pL
d,all,ux *)]"27^
!施加重力加速度 aCV4Ay G
acel,,9.8 F_uY{bg
(4)求解 {A|T owBN
/solu ;IK[Y{W /
!求解设置 K\XyZ
antype,static !设定为静力求解 Jx#k, Z4
nsubst,100 !设定最大子步数为100 1{_A:<VBl
pred,on !打开时间步长预测器 v +"rZ
nropt,full !设定牛顿-拉普森选项 \Ep0J $ #o
nlgeom,on !打开大位移效果 /J)l /oI
lnsrch,on !打开线性搜索 #}^-C&~
outres,all,all !输出所有项 Jw~( G9G
cnvtol,f,,0.005,2,0.5 !力收敛准则设定 6mH/ m&
cnvtol,u,,0.05,2,1 !位移收敛准则设定 ``ekR6[ 8c
!边坡在强度折减系数F=1时求解 W0jZOP5_.$
allsel *Ywpz^2?:
solve !进行求解 7kKy\W
save,F1,db !把F=1时求解结果保存 oPc\<$
!边坡在强度折减系数F=1.2时求解 L}#0I+M l7
finish 4 (l?uU$
/solu 4lKVY<
allsel
htY=w}>
asel,s,area,,2,3,1 !选择面积A2、A3 vILy>QS)
MPCHG,3,all ! 把所选择单元材料号改为3 *c[2C
allsel x_|F|9
solve !进行求解 S]sk7
save,F1.2,db !把F=1.2时求解结果保存 ~7!7\i,Y8\
!边坡在强度折减系数F=1.4时求解 |lH;Fq {\
finish 5BR 5X\f0
/solu j'i0*"x
allsel juBw5U <
asel,s,area,,2,3,1 !选择面积A2、A3 D }3XFuZs_
MPCHG,4,all !把所选择单元材料号改为4 ;d$qc<2uA
allsel 6a}"6d/sTL
solve !进行求解 B5X sGLV
save,F1.4,db !把F=1.4时求解结果保存 v6\F
Q9|t
!边坡在强度折减系数F=1.6时求解 J/);"bg_O
finish p1c3Q$>i
/solu
wiX ~D
allsel X}P$e mr7
asel,s,area,,2,3,1 !选择面积A2、A3 9{j66
MPCHG,5,all !把所选择单元材料号改为5
>ds%].$-\
allsel ENh!N4vbO
solve 0tk#Gs[
save,F1.6,db !把F=1.6时求解结果保存 @xsCXCRWVV
!边坡在强度折减系数F=1.8时求解 V Cy5JH
finish Z['\61
/solu l:]Nn%U(>
allsel M\b")Tu{0
asel,s,area,,2,3,1 !选择面积A2、A3 ~8|t*@ D
MPCHG,6,all !把所选择单元材料号改为6
QH]M
allsel gqib:q ;r
solve ~tB;@ e
save,F1.8,db !把F=1.8时求解结果保存 W\f9jfD
!边坡在强度折减系数F=2.0时求解 #<( = }?
finish (yo;NKq,@
/solu c^8o~K>w84
allsel <ktz T&A
asel,s,area,,2,3,1 !选择面积A2、A3 +*oS((0s
MPCHG,7,all !把所选择单元材料号改为7 -oz`"&%
allsel d +iR/Ssc
solve ^BZkHAp
save,F2.0,db !把F=2.0时求解结果保存 SvvUkQ#1w
!边坡在强度折减系数F=2.2时求解 bU 63X={
finish TgU**JN)
/solu C} IbxKl
allsel 6B$q,"%S@
asel,s,area,,2,3,1 !选择面积A2、A3 n 3MWs );5
MPCHG,8,all !把所选择单元材料号改为8 >*twTlb{
allsel vFrt|JC_{
solve 6{d6 s#|%
save,F2.2,db !把F=2.2时求解结果保存 ac d:r%y
!边坡在强度折减系数F=2.4时求解 U-wLt(Y<
finish
t4Z
/solu t)oa pIeIe
allsel
O?EB8RB
asel,s,area,,2,3,1 !选择面积A2、A3 6pE :A@
MPCHG,9,all !把所选择单元材料号改为9 4\. V
allsel ^ 0W(hA
solve EPW7+Ve
save,F2.4,db !把F=2.4时求解结果保存 /RLq>#:h**
!边坡在强度折减系数F=2.6时求解 c':ezEaC
finish `nR %Cav,U
/solu C9S@v D+
allsel t<:D@J]a
asel,s,area,,2,3,1 !选择面积A2、A3 W&:[r/8wA
MPCHG,10,all !把所选择单元材料号改为10 #0b&^QL
allsel cD1o"bq
solve 3>z[PPw
save,F2.6,db !把F=2.6时求解结果保存 &$`hQgi
!边坡在强度折减系数F=2.8时求解 ;evCW$G=
finish q ^rl)
/solu 0e ["]Tlnm
allsel k&hc m
asel,s,area,,2,3,1 !选择面积A2、A3 l6[lJ0Y
MPCHG,11,all !把所选择单元材料号改为11 2Ha5yaTL
allsel *p ? e.%nd
solve 1gO2C $
save,F2.8,db !把F=2.8时求解结果保存 $3=:E 36K
!边坡在强度折减系数F=3.0时求解 KRsAv^']
finish a=G M[{og
/solu =R*Gk4<Y
allsel "%8A :^1
asel,s,area,,2,3,1 !选择面积A2、A3 v;y0jD#b
MPCHG,12,all !把所选择单元材料号改为12 A{o ' z_zC
allsel xa( m5P
solve 3-40'$lE
save,F3.0,db !把F=3.0时求解结果保存 tNs~M4TVVH
(5）后处理 {mE! Vf
/post1 ！进入后处理 &K^MN d
!边坡在强度折减系数F=1时结果分析 p<WFqLe(":
Resume,'F1','db' !读入边坡在强度折减系数F=1时 <I;*[ ;AK
set,1,last !读入后一个子步 7=4 A;Ybq
pldisp,1 !绘制边坡模型变形图 U3vEdw<lV
plnsol,u,x !绘制边坡模型水平方向位移云图 =j0x.f Se
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 puJ#w1!x`
!边坡在强度折减系数F=1.2时结果分析 ~#:e *:ro
Resume,'F1.2','db' !读入边坡在强度折减系数F=1.2时 M
ixwK,
set,1,last !读入后一个子步 .V6-(d
pldisp,1 !绘制边坡模型变形图 0@1:M
plnsol,u,x !绘制边坡模型水平方向位移云图 E&
36 H
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 ZA#y)z8!E
!边坡在强度折减系数F=1.4时结果分析 C-P06Q]
Resume,'F1.4','db' !读入边坡在强度折减系数F=1.4时 cd ;NpN
set,1,last !读入后一个子步 c.H?4j7ga
pldisp,1 !绘制边坡模型变形图 bAxTLIf
plnsol,u,x !绘制边坡模型水平方向位移云图 ,gnQa
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 +?RGta'%k
!边坡在强度折减系数F=1.6时结果分析 LE?u`i,e=+
Resume,'F1.6','db' !读入边坡在强度折减系数F=1.6时 @E`?<| B}
set,1,last !读入后一个子步 !a1i Un9
pldisp,1 !绘制边坡模型变形图 -U 2mfW
plnsol,u,x !绘制边坡模型水平方向位移云图 T#:F]=
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 !Q-wdzsp?
!边坡在强度折减系数F=1.8时结果分析 p@&R0>6j
Resume,'F1.8','db' !读入边坡在强度折减系数F=1.8时 ]Gc3Ea;4
set,1,last !读入后一个子步 j_?cpm{~ml
pldisp,1 !绘制边坡模型变形图 g ( 0;[#@
plnsol,u,x !绘制边坡模型水平方向位移云图 FgA//)1
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 xSpC'"
!边坡在强度折减系数F=2.0时结果分析 $7I] `Jt
Resume,'F2.0','db' !读入边坡在强度折减系数F=2.0时 k7_I$ <YDj
set,1,last !读入后一个子步 _8K%`6!"Z
pldisp,1 !绘制边坡模型变形图 |c -LSs'\
plnsol,u,x !绘制边坡模型水平方向位移云图 %R0v5=2'
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 Oi:JiD=
!边坡在强度折减系数F=2.2时结果分析 qUhRu>
Resume,'F2.2','db' !读入边坡在强度折减系数F=2.2时 D1,O:+[;.
set,1,last !读入后一个子步 .
,NB( s`
pldisp,1 !绘制边坡模型变形图 Kn+=lCk
plnsol,u,x !绘制边坡模型水平方向位移云图 IKvBf'%-
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 wCZO9sU:6=
!边坡在强度折减系数F=2.4时结果分析 ST1Ts 5I
Resume,'F2.4','db' !读入边坡在强度折减系数F=2.4时 QL"gWr`R
set,1,last !读入后一个子步 *2u
E
pldisp,1 !绘制边坡模型变形图 U6 82 Th
plnsol,u,x !绘制边坡模型水平方向位移云图 _J?SIm
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 ?SY<~i<K-
!边坡在强度折减系数F=2.6时结果分析 zW{ 6Eg
Resume,'F2.6','db' !读入边坡在强度折减系数F=2.6时 FYPz 4K
set,1,last !读入后一个子步 #`GbHxd
pldisp,1 !绘制边坡模型变形图 E(+T*
plnsol,u,x !绘制边坡模型水平方向位移云图 }wt%1v-10U
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 7x. %hRk
!边坡在强度折减系数F=2.8时结果分析 a j|5 #
Resume,'F2.8','db' !读入边坡在强度折减系数F=2.8时 pt:;9hA
set,1,last !读入后一个子步 Q 882B1H
pldisp,1 !绘制边坡模型变形图 Pn TZ/|
plnsol,u,x !绘制边坡模型水平方向位移云图 r
-f
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 jeN1eM8 WI
!边坡在强度折减系数F=3.0时结果分析 SA6hbcYk
Resume,'F3.0','db' !读入边坡在强度折减系数F=3.0时 !>;w!^U
set,1,last !读入后一个子步 FyD.>ot7M
pldisp,1 !绘制边坡模型变形图 %|3e.1oX
plnsol,u,x !绘制边坡模型水平方向位移云图 2xmk,&s
plnsol,eppl,eqv !绘制边坡模型塑性应变云图
