各位朋友： ~4iI G}Y<
大家好！我通过ansys用有限元强度折减法分析了一个边坡的稳定性，结果不同的工况下得出的塑性区是一样的，我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激！附命令流和塑性区如下：(1)创建物理环境 U {>eE8l
/COM, Structural ! 指定结构分析 (&$|R\W.
/TITLE, Slope stability Analysis ! 定义工作标题
0fP qO2
/FILNAM,support,1 ! 定义工作文件名 ^eEj
5Rh
(2)建立模型 g ~@0p7]Y
!进入前处理器 xW,(d5RtZ
/PREP7 , &HZvU&
!定义单元类型 3Ec m Nwr
ET,1,PLANE82 !定义边坡围岩单元 Zh.fv-Ecp
KEYOPT,1,3,2 !求解类型为平面应变 Al`[Iu&
!定义材料属性 sy` : wp
!边坡围岩1材料属性 hqY9\,.C
MP,EX,1,3E10 [hS?d.D
MP,PRXY,1,0.25
Xze
MP,DENS,1,2500 !定义塑性模型 05MtQB
TB,DP,1 !定义数据表1 3Fn26Ri j
TBDATA,1,9E5,42.8 !给数据表输入数据 v|Tg %
!边坡围岩2材料属性 !采用弹性模型 wDZFOx0#8
MP,EX,2,3.2E10 ！定义弹性模量 G1~|$X@@
MP,PRXY,2,0.24 !定义泊松比 WU1 o4& OF
MP,DENS,2,2700 !定义密度 kI3-G~2
!取强度折减系数F=1.2时边坡围岩材料属性 b~b(Ed{r
MP,EX,3,3E10 l1N{ ujM
MP,PRXY,3,0.25 lq_W;L
MP,DENS,3,2500 z ux+ooU
TB,DP,3 !定义数据表3 W &4`eB/4}
TBDATA,1,7.5E5,37.7 !给数据表输入数据 v-) eT
!取强度折减系数F=1.4时边坡围岩材料属性 XI\aZ\v
MP,EX,4,3E10 7-ba-[t#A
MP,PRXY,4,0.25 !定义泊松比 9,'5~+7
MP,DENS,4,2500 *<U&DOYV:
TB,DP,4 !定义数据表4
jl2nRo
TBDATA,1,6.4E5,33.5 !给数据表输入数据 FV$= l
%
!取强度折减系数F=1.6时边坡围岩材料属性 lrh6lt)
MP,EX,5,3E10 *4qsM, t
MP,PRXY,5,0.25 &bj :,$@
MP,DENS,5,2500 ！定义密度 Z=!*7@QY
TB,DP,5 !定义数据表5 A0rdQmrOL
TBDATA,1,5.6E5,30 !给数据表输入数据 ?/)5U}*M0T
!取强度折减系数F=1.8时边坡围岩材料属性 "`"j2{9|e!
MP,EX,6,3E10 ！定义弹性模量 {zUc*9
MP,PRXY,6,0.25 gqD`1/
MP,DENS,6,2500 @\!wW-:A
TB,DP,6 !定义数据表6 '{ cN~A2b4
TBDATA,1,5.0E5,27.2 0v0Y(
Mo@
!取强度折减系数F=2.0时边坡围岩材料属性 vo`wYJ3W
MP,EX,7,3E10 `z` `d*_
MP,PRXY,7,0.25 !定义泊松比 y=vH8D]%X
MP,DENS,7,2500 R C (v# G
TB,DP,7 !定义数据表7 R /^JyL
TBDATA,1,4.5E5,24.8 |7`Vw Z
!取强度折减系数F=2.2时边坡围岩材料属性 -t6d`p;dR
MP,EX,8,3E10 [l #WS
MP,PRXY,8,0.25 @LL&ggV?
MP,DENS,8,2500 pN0c'COy^
TB,DP,8 !定义数据表8 45MK|4\Y_
TBDATA,1,4.09E5,22.8 : 6>H\
!取强度折减系数F=2.4时边坡围岩材料属性 P>@ `hZ9
o
MP,EX,9,3E10 ;wQWt_OtuJ
MP,PRXY,9,0.25 <YeF?$S}
MP,DENS,9,2500 }7=a,1T
TB,DP,9 !定义数据表9 ,uKvE`H
TBDATA,1,3.6E5,21.1 !给数据表输入数据 7Xad2wXn
!取强度折减系数F=2.6时边坡围岩材料属性 HqXo;`Yy}
MP,EX,10,3E10 ！定义弹性模量 \;7 DS:d@
MP,PRXY,10,0.25 !定义泊松比
eN>
(IW
MP,DENS,10,2500 ！定义密度 HnP ;1Gi
TB,DP,10 !定义数据表10 eF8`an5S
TBDATA,1,3.46E5,19.6 !给数据表10输入数据 |}M']V z
!取强度折减系数F=2.8时边坡围岩材料属性 7-o=E =
MP,EX,11,3E10 ！定义弹性模量 (u G4W|?p
MP,PRXY,11,0.25 !定义泊松比 R\.huOJh
MP,DENS,11,2500 Q"hI !PO+
TB,DP,11 !定义数据表11 i 9!Urq-
TBDATA,1,3.2E5,18.2 !给数据表11输入数据 "j a0,%3
!取强度折减系数F=3.0时边坡围岩材料属性 Ve8!
MP,EX,12,3E10 ！定义弹性模量 j^%i?BWw
MP,PRXY,12,0.25 |O4A+S
MP,DENS,12,2500 ！定义密度 @eT sS%f2
TB,DP,12 !定义数据表12 ax^${s|{-
TBDATA,1,3.0E5,17.1 !给数据表12输入数据 K%pmE?%,8
!建立几何模型 Tc88U8Gc
!创建边坡线模型线 Y$ChMf
K,1,,,, ！ 创建关键点1 z}Z`kq+C
K,2,-800,,, 5R=lTx/Hj
K,3,-800,-800,, #h3+T*5} 6
K,4,-800,-1200,, xeI ,Kz."
K,5,1200,-1200,, +1)C&:
K,6,1200,-800,, 9wq%Fnt
K,7,1200,,, @&Nvb.5nT
K,8,1200,400,, @^` <iTK&p
K,9,400,400,, o` 77gkLO
l,1,2 !连接关键点1、2生成直线L1 G wW#Ww;Oc
l,2,3 y2 +a2
l,3,4 H~ZSw7!M8
l,4,5 !%u#J:z2
l,5,6 li%-9Jd
l,6,7 42wcpSp
l,7,8 i9xv`Ev=R
l,8,9 nUd\4;J#
l,9,1 6 gKOpa
l,1,7 R82Y&s;
l,3,6 } ) e`0)
!创建边坡面模型 %UQ{'JW?K
al,3,4,5,11 !由线L3、L4、L5、L11生成一个面A1 OpL 6Y+<
al,1,2,11,6,10 !由线L1、L2、L1、L6、L10生成一个面A2 uWWv`bI>x
al,7,8,9,10 !由线L7、L8、L9、L10生成一个面A3 P+b^;+\1s
/PNUM, area, on !打开面号开关 s{k\1 P(G}
!划分网格生成有限元模型 q9
;\ B&
! 划分边坡围岩2网格 t<|s &
!设置网格份数 lXPn]iLJ
lsel,s,,,3,5,2 !选择线L3和L5 zl:D|h77
lesize,all,,,5 !把所选择线分为5段 vv0zUvmT
lsel,s,,,4,11,7 !选择线L4和L11 5{g9Wh[
lesize,all,,,26 !把所选择线分为26段 Pu^~]^W)
mat,2 !给边坡围岩2赋予2号材料特性 DrEtnt
type,1 !采用单元类型1 $0*47+f
mshkey,0 !设定自由网格划分 $dnHUBB
mshape,0 !设定四边形网格划分 {6n \532@
amesh,1 ！划分面积A1 V 1:3
! 划分边坡围岩1网格 ;:|KfXiC8
!设置网格份数 jkP70Is
lsel,s,,,2,6,4 !选择线L2和L6 W<LaR,7
lsel,a,,,1 !附加选择线L1 6J
5)4^bk
lesize,all,,,10 !把所选择线分为10段 j?,*fp8
lsel,s,,,8,10,2 !选择线L8和L10 V[I<9xaE
lesize,all,,,16 !把所选择线分为16段 0l6iv[qu5w
lsel,s,,,7,9,2 K%} I}8M
lesize,all,,,12 rl^LS z
mat,1 !给边坡围岩1赋予1号材料特性 v 2 GhR*
type,1 d*>k
]X@G
mshkey,0 !设定自由网格划分 ziycyf.d
mshape,0 !设定四边形网格划分 "? {yVu~9
amesh,2 !划分面积A2 /a(zLHyz)
amesh,3 ！划分面积A3 BLH=:zb5
allsel Dq?HUb^X
!保存网格模型 u}1 vn} F{
save,Slope-grid.db ~yi&wbTjM
(3)施加约束和荷载 0w ['jh|,
!给边坡模型施加约束 D%idlL2%J
!边坡两侧施加X方向约束 BZ <z@DJp
nsel,s,loc,x,-800 !选择X=-800线上所有节点 k@aP&Z~
nsel,a,loc,x,1200 !选择X=1200线上所有节点 %~ 5Q^3$O
d,all,ux !对所选择节点约束X方向位移 4) 3pa *
allsel ? 03Zy3 /
!边坡底部施加约束 D&r8V;G[[
nsel,s,loc,y,-1200 !选择Y=--1200线上所有节点 b(q& }60
d,all,uy !对所选择节点约束X、Y方向位移 7h]R{ _
d,all,ux Qg>GW
!施加重力加速度 "|'`'W
acel,,9.8 &i Otw0E
(4)求解 x#0 @ $
/solu 6 h_ k`z
!求解设置 <e&QTyb
antype,static !设定为静力求解 'T rrOq4
nsubst,100 !设定最大子步数为100 _ <a)\UR
pred,on !打开时间步长预测器 i 9tJHeSm
nropt,full !设定牛顿-拉普森选项 0o| ,& K
nlgeom,on !打开大位移效果 D:ugP ,
lnsrch,on !打开线性搜索 `7%eA9*.m
outres,all,all !输出所有项 *w B-lg7%
cnvtol,f,,0.005,2,0.5 !力收敛准则设定 DiQkT R
cnvtol,u,,0.05,2,1 !位移收敛准则设定 (1,#=e+
!边坡在强度折减系数F=1时求解 G^ ZkY
allsel /mwr1GU
solve !进行求解 HCIF9{o1j>
save,F1,db !把F=1时求解结果保存 +Jka :]MW!
!边坡在强度折减系数F=1.2时求解 O|&TL9:
finish fF > qU-
/solu m9 0R8 V
allsel r)5 \3j[P
asel,s,area,,2,3,1 !选择面积A2、A3 $ Xs Q e
MPCHG,3,all ! 把所选择单元材料号改为3 FFeRE{,
allsel 0bx vM
solve !进行求解 A4Q)YY9~
save,F1.2,db !把F=1.2时求解结果保存 ZU.) K>'
!边坡在强度折减系数F=1.4时求解 Kc[u}
. U
finish xnQGCw?S&}
/solu a61eH )a
allsel DVoV:pk
asel,s,area,,2,3,1 !选择面积A2、A3 f76|
MPCHG,4,all !把所选择单元材料号改为4 ,L{o, qzC
allsel a^T 4\
solve !进行求解 <SM&VOiaOz
save,F1.4,db !把F=1.4时求解结果保存 Z c3:9
!边坡在强度折减系数F=1.6时求解 q AVfbcb
finish inv{dg/2
/solu C $RAJ
allsel SxW}Z_8x
asel,s,area,,2,3,1 !选择面积A2、A3 !bV(VRbu
MPCHG,5,all !把所选择单元材料号改为5 mYjiiql~
allsel K@ sP~('
solve Rf(x^J {
save,F1.6,db !把F=1.6时求解结果保存 YeExjC
!边坡在强度折减系数F=1.8时求解 e N< pU%7
finish mz VuQ
/solu 0HJqsSZ$mW
allsel 0"<;You
asel,s,area,,2,3,1 !选择面积A2、A3 3Ra\2(bR
MPCHG,6,all !把所选择单元材料号改为6 Cb%.C;q
allsel D@( M+u9/%
solve ( {_:^$E\
save,F1.8,db !把F=1.8时求解结果保存 Ul713Bjz
!边坡在强度折减系数F=2.0时求解 FX|lhwmc(
finish V0'p1J tD
/solu LJNie*
allsel d+m6-4[_k
asel,s,area,,2,3,1 !选择面积A2、A3 ;p#Z :6
MPCHG,7,all !把所选择单元材料号改为7 tD~PvUJ
allsel Xq^y<[
solve M`C~6Mf+
save,F2.0,db !把F=2.0时求解结果保存 e7bT%h9i
!边坡在强度折减系数F=2.2时求解 p]V-<
finish .f !]@"\
/solu (LT\
IJSM
allsel ^8-~@01.`_
asel,s,area,,2,3,1 !选择面积A2、A3 X |1_0
MPCHG,8,all !把所选择单元材料号改为8 4'{hI;&a&
allsel :&1=8^B Y
solve >.O*gv/ _
save,F2.2,db !把F=2.2时求解结果保存 ngkeJ)M0$
!边坡在强度折减系数F=2.4时求解 O^R:_vb3I
finish Ss<_K>wk
/solu 3dxnh,]&@
allsel emkMR{MY
asel,s,area,,2,3,1 !选择面积A2、A3 ny
MPCHG,9,all !把所选择单元材料号改为9 l\sS?
allsel >0<KkBH
solve L=$?q/=-
save,F2.4,db !把F=2.4时求解结果保存 @RoRNat
!边坡在强度折减系数F=2.6时求解 nT@6g|!
finish 43u PH1
)
/solu (w#s lTFT
allsel / 8WpX
asel,s,area,,2,3,1 !选择面积A2、A3 @ O<@f8-
MPCHG,10,all !把所选择单元材料号改为10 Y( V3P nH
allsel "doU.U&u
solve d#vS E.&
save,F2.6,db !把F=2.6时求解结果保存 iFwyh`Bcg
!边坡在强度折减系数F=2.8时求解 Oa.f~|
finish *_ {l
/solu |>IUtUg\
allsel $||ns@F+
asel,s,area,,2,3,1 !选择面积A2、A3 i!.I;@
MPCHG,11,all !把所选择单元材料号改为11 oc3dd "8}@
allsel BS-:dyBw
solve Rg8m4x w
save,F2.8,db !把F=2.8时求解结果保存 @`IXu$ Wm(
!边坡在强度折减系数F=3.0时求解 5!SoN}$
finish \`P2Yq
/solu (pT(&/\8
allsel `b Z 2x@
asel,s,area,,2,3,1 !选择面积A2、A3 o% /-5-
MPCHG,12,all !把所选择单元材料号改为12 /jZaU`
allsel cXOje"5i
solve 3X A p Y '
save,F3.0,db !把F=3.0时求解结果保存 '&.)T 2Kw
(5）后处理 z<yNG/M1>U
/post1 ！进入后处理 SLQ\Y%F
!边坡在强度折减系数F=1时结果分析 C-a*EG
Resume,'F1','db' !读入边坡在强度折减系数F=1时 3WO#^}t
set,1,last !读入后一个子步 L Xh@o1
pldisp,1 !绘制边坡模型变形图 F=-uDtQ <N
plnsol,u,x !绘制边坡模型水平方向位移云图 T bKP8zw{
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 &J lpA<^s;
!边坡在强度折减系数F=1.2时结果分析 j.=UI-&m
Resume,'F1.2','db' !读入边坡在强度折减系数F=1.2时 l50|`
6t
set,1,last !读入后一个子步 z)58\rtz
pldisp,1 !绘制边坡模型变形图 e(O wS?K
plnsol,u,x !绘制边坡模型水平方向位移云图 ( Lok
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 x9x #'H3
!边坡在强度折减系数F=1.4时结果分析 2Y E;m&
Resume,'F1.4','db' !读入边坡在强度折减系数F=1.4时 DD12pL{QA
set,1,last !读入后一个子步 C=oM,[ESQ0
pldisp,1 !绘制边坡模型变形图 :rz9M@7
plnsol,u,x !绘制边坡模型水平方向位移云图 9*}iBs
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 h*;g0QBkl
!边坡在强度折减系数F=1.6时结果分析 L;1$xI8tx
Resume,'F1.6','db' !读入边坡在强度折减系数F=1.6时 #D|!
.I)
set,1,last !读入后一个子步 " 8wf.nZ
pldisp,1 !绘制边坡模型变形图 ] ?DDCew
plnsol,u,x !绘制边坡模型水平方向位移云图 E3 ~,+68U
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 w+Z- -@\
!边坡在强度折减系数F=1.8时结果分析 `y8
? =
Resume,'F1.8','db' !读入边坡在强度折减系数F=1.8时 ZdE>C
set,1,last !读入后一个子步 M >*0r<qn
pldisp,1 !绘制边坡模型变形图 /<3;0~#){
plnsol,u,x !绘制边坡模型水平方向位移云图 54
M!Fq -
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 EX`"z(L
!边坡在强度折减系数F=2.0时结果分析 0t# NMW
Resume,'F2.0','db' !读入边坡在强度折减系数F=2.0时 u @{E{
set,1,last !读入后一个子步 P &;y]
,)E
pldisp,1 !绘制边坡模型变形图 iGG6Myp-
plnsol,u,x !绘制边坡模型水平方向位移云图 ;;ED N45
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 9">zdFC'
!边坡在强度折减系数F=2.2时结果分析 ioB|*D<U2
Resume,'F2.2','db' !读入边坡在强度折减系数F=2.2时 0K=Qf69Y
set,1,last !读入后一个子步 Ys"|</;dbj
pldisp,1 !绘制边坡模型变形图 B#HV20\?v
plnsol,u,x !绘制边坡模型水平方向位移云图 !Gln Q`T
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 l,*yEkU
!边坡在强度折减系数F=2.4时结果分析 =)(3Dp
Resume,'F2.4','db' !读入边坡在强度折减系数F=2.4时 89LpklD
set,1,last !读入后一个子步 !MZ+- dpK
pldisp,1 !绘制边坡模型变形图 1a7!4)\
plnsol,u,x !绘制边坡模型水平方向位移云图 ZWf-X
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 qSFc=Wwc
!边坡在强度折减系数F=2.6时结果分析 lhLnyg Uk
Resume,'F2.6','db' !读入边坡在强度折减系数F=2.6时 <U@P=G< t
set,1,last !读入后一个子步 >Y 7r \
pldisp,1 !绘制边坡模型变形图 "KMLk
plnsol,u,x !绘制边坡模型水平方向位移云图 ;EgzC^2e
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 mND XzT&
!边坡在强度折减系数F=2.8时结果分析 z=ML(1c=
Resume,'F2.8','db' !读入边坡在强度折减系数F=2.8时 6BDt.bG
set,1,last !读入后一个子步 |0tg:\.
pldisp,1 !绘制边坡模型变形图 UppBnw
plnsol,u,x !绘制边坡模型水平方向位移云图 `d75@ 0:
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 cqeR<len
!边坡在强度折减系数F=3.0时结果分析 w-
UKMW9"
Resume,'F3.0','db' !读入边坡在强度折减系数F=3.0时 :=rA Yc3]
set,1,last !读入后一个子步 Q Oz9\,C
pldisp,1 !绘制边坡模型变形图 m_c O<LB
plnsol,u,x !绘制边坡模型水平方向位移云图 CD^CUbGk
plnsol,eppl,eqv !绘制边坡模型塑性应变云图
