各位朋友： $_-f } E
大家好！我通过ansys用有限元强度折减法分析了一个边坡的稳定性，结果不同的工况下得出的塑性区是一样的，我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激！附命令流和塑性区如下：(1)创建物理环境 x+zz:^yHYf
/COM, Structural ! 指定结构分析 6! } @vp![
/TITLE, Slope stability Analysis ! 定义工作标题 N|c;Qzl
/FILNAM,support,1 ! 定义工作文件名 [&n2 yt
(2)建立模型 5Sd+ Cc
!进入前处理器 Egv (n@1
/PREP7 {g?$u
!定义单元类型 uBLI!N-G
ET,1,PLANE82 !定义边坡围岩单元 X>w( ^L*>
KEYOPT,1,3,2 !求解类型为平面应变 /TTmMx*
!定义材料属性 '>AOJ aA
!边坡围岩1材料属性 zbdmz
MP,EX,1,3E10 :E/]Bjq$;
MP,PRXY,1,0.25 8kQ
>M
MP,DENS,1,2500 !定义塑性模型 $m`?x5rL8
TB,DP,1 !定义数据表1 xz g81sV7
TBDATA,1,9E5,42.8 !给数据表输入数据 _tZT
!边坡围岩2材料属性 !采用弹性模型 =[!&&,c=
MP,EX,2,3.2E10 ！定义弹性模量 A4. Q \0
MP,PRXY,2,0.24 !定义泊松比 *TY?*H
MP,DENS,2,2700 !定义密度 d\}r.pD
!取强度折减系数F=1.2时边坡围岩材料属性 4I7}
MP,EX,3,3E10 B&&:A4
MP,PRXY,3,0.25 -QL _a8NL
MP,DENS,3,2500 ahNpHTPa
TB,DP,3 !定义数据表3 q.0a0 /R
TBDATA,1,7.5E5,37.7 !给数据表输入数据 M] YK]VyG
!取强度折减系数F=1.4时边坡围岩材料属性 * 8n0
MP,EX,4,3E10 K_
P08
MP,PRXY,4,0.25 !定义泊松比 V3WHp'1
MP,DENS,4,2500 Zr, :i
MPZ
TB,DP,4 !定义数据表4 w-l:* EV8
TBDATA,1,6.4E5,33.5 !给数据表输入数据 : T(3!}4
!取强度折减系数F=1.6时边坡围岩材料属性 PSz|I8
c
MP,EX,5,3E10 VfOm#Ue0 q
MP,PRXY,5,0.25 smQ<lwA
MP,DENS,5,2500 ！定义密度 4S>A}rWz
TB,DP,5 !定义数据表5 0R&$P 6
TBDATA,1,5.6E5,30 !给数据表输入数据 [(U:1&x &
!取强度折减系数F=1.8时边坡围岩材料属性 gT$`a
MP,EX,6,3E10 ！定义弹性模量 (j;s6g0
MP,PRXY,6,0.25 hER]%)#r
MP,DENS,6,2500 s^F6sXhyPi
TB,DP,6 !定义数据表6 ~G.'pyW
TBDATA,1,5.0E5,27.2 v'(p."g
!取强度折减系数F=2.0时边坡围岩材料属性 -@{5
u d
MP,EX,7,3E10 %EA|2O.D
MP,PRXY,7,0.25 !定义泊松比 jovI8Dw
>
MP,DENS,7,2500 &U%AV D[
TB,DP,7 !定义数据表7 OnE#8*8
TBDATA,1,4.5E5,24.8 LIYj__4=|
!取强度折减系数F=2.2时边坡围岩材料属性 n} cjVH5
MP,EX,8,3E10 [F4] pR(
MP,PRXY,8,0.25 H'E(gc)>)
MP,DENS,8,2500 (L z VWz m
TB,DP,8 !定义数据表8 Zq7Y('=`t@
TBDATA,1,4.09E5,22.8 _]btsv\)f
!取强度折减系数F=2.4时边坡围岩材料属性 $_N Yu
MP,EX,9,3E10 "w. gP8`
MP,PRXY,9,0.25 $d 2mcwh\
MP,DENS,9,2500 Mz,G; x}
TB,DP,9 !定义数据表9 aI]EwVz-q
TBDATA,1,3.6E5,21.1 !给数据表输入数据 'E FP/(2J
!取强度折减系数F=2.6时边坡围岩材料属性 55 5j@
MP,EX,10,3E10 ！定义弹性模量 Y=+pz^/"
MP,PRXY,10,0.25 !定义泊松比 1sgI,5liUs
MP,DENS,10,2500 ！定义密度 Hd
H ,
TB,DP,10 !定义数据表10 t]LOBy-Kv
TBDATA,1,3.46E5,19.6 !给数据表10输入数据 I;xrw?=\L
!取强度折减系数F=2.8时边坡围岩材料属性 JGl0
(i*|
MP,EX,11,3E10 ！定义弹性模量 4X tIMa28
MP,PRXY,11,0.25 !定义泊松比 cKt=_4Lf
MP,DENS,11,2500 53efF bo
TB,DP,11 !定义数据表11 Q-<N)K$F(4
TBDATA,1,3.2E5,18.2 !给数据表11输入数据 )Y.H*ca
!取强度折减系数F=3.0时边坡围岩材料属性 hKo& ZWPq
MP,EX,12,3E10 ！定义弹性模量 lezX-5Z
MP,PRXY,12,0.25 %mKM9>lf#
MP,DENS,12,2500 ！定义密度 s`Z |
A
TB,DP,12 !定义数据表12 F
uYjrzmx
TBDATA,1,3.0E5,17.1 !给数据表12输入数据 g{{D C )>
!建立几何模型 SVKjhZK
!创建边坡线模型线 ,*_=w^;Rr
K,1,,,, ！ 创建关键点1 jpijnz{M
K,2,-800,,, QP HibPP:
K,3,-800,-800,, A=8%2U wI
K,4,-800,-1200,, MvV\?Lzj
K,5,1200,-1200,, G`1!SEae
K,6,1200,-800,, ! R@v\Eu
K,7,1200,,, Sep/N"7~t
K,8,1200,400,, (R^X3
K,9,400,400,, &c[.&L,w4
l,1,2 !连接关键点1、2生成直线L1 8SH&b8k<<
l,2,3 FfEP@$
l,3,4 yR?S]
l,4,5 5dwC~vn}c
l,5,6 :(7icHa
l,6,7 KV&4Ep#
l,7,8 O^/z 7,
l,8,9 2z*EamF
l,9,1 _,r2g8qm
l,1,7 ,j3Yvn W
l,3,6 VINb9W}G[
!创建边坡面模型 4Wsp PHj
al,3,4,5,11 !由线L3、L4、L5、L11生成一个面A1 hbfN 1 "z
al,1,2,11,6,10 !由线L1、L2、L1、L6、L10生成一个面A2
mE1m
al,7,8,9,10 !由线L7、L8、L9、L10生成一个面A3 D1G9^7:^E
/PNUM, area, on !打开面号开关 _=Y?' gHH
!划分网格生成有限元模型 ":ig Yh
! 划分边坡围岩2网格 p;nRxi7'
!设置网格份数 J ZH~ {
lsel,s,,,3,5,2 !选择线L3和L5 e'nhP
lesize,all,,,5 !把所选择线分为5段 pAdx 6
lsel,s,,,4,11,7 !选择线L4和L11 faTp|T`nY
lesize,all,,,26 !把所选择线分为26段 Xg <R+o
mat,2 !给边坡围岩2赋予2号材料特性 sKu/VAh
x
type,1 !采用单元类型1 wlC7;u
mshkey,0 !设定自由网格划分 Cpg>5N~;L
mshape,0 !设定四边形网格划分
(zIW JJw
amesh,1 ！划分面积A1 #E`wqI\'
! 划分边坡围岩1网格 =[_=y=G
!设置网格份数 YB`1S
lsel,s,,,2,6,4 !选择线L2和L6 WR :I2-1
lsel,a,,,1 !附加选择线L1 /+
yIcE(&3
lesize,all,,,10 !把所选择线分为10段 Kg8n3pLAX
lsel,s,,,8,10,2 !选择线L8和L10 LZ'Y3 *
lesize,all,,,16 !把所选择线分为16段 p)&Yr
lsel,s,,,7,9,2 6j 2mr6o
lesize,all,,,12 -d+o\qp"#
mat,1 !给边坡围岩1赋予1号材料特性 .:wo
ARW!
type,1 Pl/}`H: R&
mshkey,0 !设定自由网格划分 !_i;6UVG
mshape,0 !设定四边形网格划分 lEXER^6
amesh,2 !划分面积A2 eN,6p '&
amesh,3 ！划分面积A3 9_Ws8nE
allsel ' lWgHmE
!保存网格模型 .X.,.vHx
save,Slope-grid.db @sPuc.
(3)施加约束和荷载 ]=Tle&yM+T
!给边坡模型施加约束 YP 4lizs.
!边坡两侧施加X方向约束 *!5X!\e_
nsel,s,loc,x,-800 !选择X=-800线上所有节点 IIh \ d.o
nsel,a,loc,x,1200 !选择X=1200线上所有节点 jA' 7@/F/
d,all,ux !对所选择节点约束X方向位移 br3r!Vuz/-
allsel (qyT,K8
!边坡底部施加约束 H$M{thW
nsel,s,loc,y,-1200 !选择Y=--1200线上所有节点 UlD]!5NO
d,all,uy !对所选择节点约束X、Y方向位移 ^l8&y;-T
d,all,ux YJ m64H,[
!施加重力加速度 ]tsp}M@
acel,,9.8 & @ U)
(4)求解 wg}rMJoG|
/solu [sc4ULS &
!求解设置 U32&"&";c
antype,static !设定为静力求解 LW#$%}
nsubst,100 !设定最大子步数为100 ]1)#Y
pred,on !打开时间步长预测器 ~ TfN*0
nropt,full !设定牛顿-拉普森选项 N~!
G AaD
nlgeom,on !打开大位移效果 a<C J#B2K
lnsrch,on !打开线性搜索 Fi8#r)G.
outres,all,all !输出所有项 P9f,zM-
cnvtol,f,,0.005,2,0.5 !力收敛准则设定 -qs
R,H
cnvtol,u,,0.05,2,1 !位移收敛准则设定 E``\Jre@
!边坡在强度折减系数F=1时求解 YZ(tjIgQ
allsel Nc_Qd4<[@G
solve !进行求解 EC8 Fapy
save,F1,db !把F=1时求解结果保存 ^ 3O `8o
!边坡在强度折减系数F=1.2时求解 N@x5h8
finish V2oXg
/solu T`W37fz0
allsel I&Z+FL&@f
asel,s,area,,2,3,1 !选择面积A2、A3 ?' /#Gt`
MPCHG,3,all ! 把所选择单元材料号改为3 oR@1/lV
allsel
% G>
solve !进行求解 LXq0hI
save,F1.2,db !把F=1.2时求解结果保存 #>_5PdO
!边坡在强度折减系数F=1.4时求解 Q&`if
O
finish @g%^H )T
/solu f7v|N)
allsel AYHefAF<w
asel,s,area,,2,3,1 !选择面积A2、A3 j`l'Mg
MPCHG,4,all !把所选择单元材料号改为4 "?Yf3G: \0
allsel t2_ pwd*B
solve !进行求解 hOM#j
save,F1.4,db !把F=1.4时求解结果保存 >k9W+mk
!边坡在强度折减系数F=1.6时求解 $*qQ/h i
finish !V27ln KP+
/solu zW:r7
P.
allsel 7*C>4Gs
asel,s,area,,2,3,1 !选择面积A2、A3 (X/dP ~
MPCHG,5,all !把所选择单元材料号改为5 X*<
!_3
allsel r}M2t$nv
solve Sej( jJX1
save,F1.6,db !把F=1.6时求解结果保存 /;;$9O9
!边坡在强度折减系数F=1.8时求解 ['_W <
finish Z@gEJ^"yA"
/solu Cm%xI& Y
allsel I*EHZctH
asel,s,area,,2,3,1 !选择面积A2、A3 ~v V+)KI
MPCHG,6,all !把所选择单元材料号改为6 zOn % \
allsel xY@< <
solve a "!r]=r
save,F1.8,db !把F=1.8时求解结果保存 [W Ud9fUL
!边坡在强度折减系数F=2.0时求解 V 7)<MY
finish 2'-o'z<
/solu vC%Hc/&.}
allsel 2@>#?c7
asel,s,area,,2,3,1 !选择面积A2、A3 |*y'H*
MPCHG,7,all !把所选择单元材料号改为7 4;)aGN{e
allsel RCY}JH>}
solve & d@N3y
save,F2.0,db !把F=2.0时求解结果保存 v {O(}@
!边坡在强度折减系数F=2.2时求解 c^8csQ fG
finish 7CSd}@71\
/solu l5<&pb#b
allsel qrkJ:
asel,s,area,,2,3,1 !选择面积A2、A3 1JY3c
M
MPCHG,8,all !把所选择单元材料号改为8 Z ItS(o J.
allsel LDY k\[81
solve {33B%5n"
save,F2.2,db !把F=2.2时求解结果保存 ;GZ/ V;S
!边坡在强度折减系数F=2.4时求解 @%gt h@8
finish (3Q G
/solu ;hCU y=m.
allsel ?vg|; Q
asel,s,area,,2,3,1 !选择面积A2、A3 T3In0LQ
MPCHG,9,all !把所选择单元材料号改为9 W3y9>]{x^
allsel }]+k
solve Z> _F:1x
save,F2.4,db !把F=2.4时求解结果保存 w3 K>IDWI7
!边坡在强度折减系数F=2.6时求解 2T3 b6
finish }HxC ~J"
/solu L%31>)8
allsel cb`ik)=K%
asel,s,area,,2,3,1 !选择面积A2、A3 VL/KC -6
MPCHG,10,all !把所选择单元材料号改为10 \PM5B"MDZ
allsel N 6w!V]b
solve I0v4T jHH
save,F2.6,db !把F=2.6时求解结果保存 ^wMZG'/
!边坡在强度折减系数F=2.8时求解 F+2 85JK
finish fea4Ul{ib
/solu ]Nb~-)t%B
allsel ;C~:C^Q\H
asel,s,area,,2,3,1 !选择面积A2、A3 >36>{b<'$*
MPCHG,11,all !把所选择单元材料号改为11 IT f4PxF
allsel vhL/L?NB$
solve @9\L|O'~?
save,F2.8,db !把F=2.8时求解结果保存 jI4 5X22j
!边坡在强度折减系数F=3.0时求解 Y%#r&de
finish xM D]b
/solu 4nU+Wj?T
allsel f# |
wb~
asel,s,area,,2,3,1 !选择面积A2、A3 (T1d!v"~"
MPCHG,12,all !把所选择单元材料号改为12 3$h yV{
allsel ;\rKkH"K8n
solve LHHDD\X
save,F3.0,db !把F=3.0时求解结果保存 u8b ^DB#+W
(5）后处理 ' qN"!\
/post1 ！进入后处理 #GIjU1-
!边坡在强度折减系数F=1时结果分析 'I5~<"E
Resume,'F1','db' !读入边坡在强度折减系数F=1时 QO <.l`F
set,1,last !读入后一个子步 VD+y4t'^
pldisp,1 !绘制边坡模型变形图 LNI]IITx /
plnsol,u,x !绘制边坡模型水平方向位移云图 G&d z<f
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 Oi, :q&
!边坡在强度折减系数F=1.2时结果分析 ;Ff5ooL{
Resume,'F1.2','db' !读入边坡在强度折减系数F=1.2时 H\<^p",`
set,1,last !读入后一个子步 1Gh3o}z
pldisp,1 !绘制边坡模型变形图 #Y<(7
plnsol,u,x !绘制边坡模型水平方向位移云图 }D*5PV%d
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 S*S @a4lV7
!边坡在强度折减系数F=1.4时结果分析 ,&0iFUwN_
Resume,'F1.4','db' !读入边坡在强度折减系数F=1.4时 on)$y&lu
set,1,last !读入后一个子步 Ou wEO
pldisp,1 !绘制边坡模型变形图 s#%P 9A
plnsol,u,x !绘制边坡模型水平方向位移云图 y({ EF~w
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 t[H A86X
!边坡在强度折减系数F=1.6时结果分析 2PG= T/
Resume,'F1.6','db' !读入边坡在强度折减系数F=1.6时 (VB-5&b
set,1,last !读入后一个子步 #6F/:j;
pldisp,1 !绘制边坡模型变形图 :WAFBK/x
plnsol,u,x !绘制边坡模型水平方向位移云图 m$ 80D,3
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 ? ?h4qJ
!边坡在强度折减系数F=1.8时结果分析 F8>J (7On
Resume,'F1.8','db' !读入边坡在强度折减系数F=1.8时 /r Z`e'}
set,1,last !读入后一个子步 ]n^TN
r7
pldisp,1 !绘制边坡模型变形图 ,n/^;. _1
plnsol,u,x !绘制边坡模型水平方向位移云图 LRK l3"M
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 BM_Rl cx~
!边坡在强度折减系数F=2.0时结果分析 } TsND6Ws3
Resume,'F2.0','db' !读入边坡在强度折减系数F=2.0时 G ,? l
o=m
set,1,last !读入后一个子步 q
!N b-O{
pldisp,1 !绘制边坡模型变形图
&b!|Y
plnsol,u,x !绘制边坡模型水平方向位移云图 3FE=?Q
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 K4j2xSGeo
!边坡在强度折减系数F=2.2时结果分析 [ @`K i
Resume,'F2.2','db' !读入边坡在强度折减系数F=2.2时
hp!. P1b
set,1,last !读入后一个子步 UQ.DKU g
pldisp,1 !绘制边坡模型变形图 mxkv{;ad
plnsol,u,x !绘制边坡模型水平方向位移云图
:JfT&YYi"
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 $p~X"f? 0
!边坡在强度折减系数F=2.4时结果分析 0mY KzJi
Resume,'F2.4','db' !读入边坡在强度折减系数F=2.4时 \. _TOE9L
set,1,last !读入后一个子步 y5$AAas
pldisp,1 !绘制边坡模型变形图 sq1v._^s
plnsol,u,x !绘制边坡模型水平方向位移云图 *|^}=ioj*
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 *;X,yEK[
!边坡在强度折减系数F=2.6时结果分析 /4bHN:I]M
Resume,'F2.6','db' !读入边坡在强度折减系数F=2.6时 KdU&q+C^
set,1,last !读入后一个子步 dBM> ;S;v
pldisp,1 !绘制边坡模型变形图
U!Eo*?LU$
plnsol,u,x !绘制边坡模型水平方向位移云图 Fa^I 1fk
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 {aVL3QU
!边坡在强度折减系数F=2.8时结果分析 '1[Bbs
Resume,'F2.8','db' !读入边坡在强度折减系数F=2.8时 mYw9lM
set,1,last !读入后一个子步 U[WR?J4~LX
pldisp,1 !绘制边坡模型变形图 bbA J5EqL
plnsol,u,x !绘制边坡模型水平方向位移云图 jp viX#\S_
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 ?*: mR|=
!边坡在强度折减系数F=3.0时结果分析 /Lt Lu
Resume,'F3.0','db' !读入边坡在强度折减系数F=3.0时 +RN|ZG&
set,1,last !读入后一个子步 Nx^r&pr
pldisp,1 !绘制边坡模型变形图 BY!M(X
jrZ
plnsol,u,x !绘制边坡模型水平方向位移云图 O\ph!?L
plnsol,eppl,eqv !绘制边坡模型塑性应变云图
2%gLq
d|RDx;r l8
v+qHH 8
$<v{$UOh
fd /?x^Z
rIo)'L$uU
3*;S%1C^
