各位朋友： A<*tn?M]
大家好！我通过ansys用有限元强度折减法分析了一个边坡的稳定性，结果不同的工况下得出的塑性区是一样的，我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激！附命令流和塑性区如下：(1)创建物理环境 0
6G[^
/COM, Structural ! 指定结构分析 l,FG:"`Z@
/TITLE, Slope stability Analysis ! 定义工作标题 9b=^"K
/FILNAM,support,1 ! 定义工作文件名 u0 'pR#
m|
(2)建立模型 7p"~:1hU
!进入前处理器 M%Kx{*aw&
/PREP7 G3^n_]Jb
!定义单元类型 ZCOuv6V+
ET,1,PLANE82 !定义边坡围岩单元 MIM PJXT#.
KEYOPT,1,3,2 !求解类型为平面应变 zxynEdO
!定义材料属性 %(wsG Nd
!边坡围岩1材料属性 =:R[gdA#1
MP,EX,1,3E10 yU&g|MV_
MP,PRXY,1,0.25 ?6l,
MP,DENS,1,2500 !定义塑性模型 bju,p"J1-E
TB,DP,1 !定义数据表1 L}sm R,
TBDATA,1,9E5,42.8 !给数据表输入数据 1PT_1[eAR
!边坡围岩2材料属性 !采用弹性模型 [7 Kj$PB3
MP,EX,2,3.2E10 ！定义弹性模量 `- uZv
MP,PRXY,2,0.24 !定义泊松比 :8GxcqvCWq
MP,DENS,2,2700 !定义密度 *7D$;?"
!取强度折减系数F=1.2时边坡围岩材料属性 /v5g;x_T
MP,EX,3,3E10 y0mg}N1
MP,PRXY,3,0.25 9Av- ;!]
MP,DENS,3,2500 \5L 4*
TB,DP,3 !定义数据表3 uhN(`E@
TBDATA,1,7.5E5,37.7 !给数据表输入数据
vG Lb2Q
!取强度折减系数F=1.4时边坡围岩材料属性 yKk,);
MP,EX,4,3E10 Ul~}@^m]4}
MP,PRXY,4,0.25 !定义泊松比 /Z9`u K
MP,DENS,4,2500 fN? Lz%z3
TB,DP,4 !定义数据表4 +
1%^c (3
TBDATA,1,6.4E5,33.5 !给数据表输入数据 n*8RYm)?
!取强度折减系数F=1.6时边坡围岩材料属性 #lVl?F+~
MP,EX,5,3E10 mYs->mg1
MP,PRXY,5,0.25 3E2.v5*
MP,DENS,5,2500 ！定义密度 NB6h/0*v
TB,DP,5 !定义数据表5 ~>-;(YU"t
TBDATA,1,5.6E5,30 !给数据表输入数据 MAhPO!e5.
!取强度折减系数F=1.8时边坡围岩材料属性 6s"bstc{
MP,EX,6,3E10 ！定义弹性模量 }mS0{rxD4
MP,PRXY,6,0.25 Jtk.v49Ad>
MP,DENS,6,2500 <?va)
ou
TB,DP,6 !定义数据表6 p-03V"^&
TBDATA,1,5.0E5,27.2 kGTc~p(
!取强度折减系数F=2.0时边坡围岩材料属性 Q8/0C b/
MP,EX,7,3E10 k~EPVJh"
MP,PRXY,7,0.25 !定义泊松比 R"kE5 :
MP,DENS,7,2500 vYm& AD
TB,DP,7 !定义数据表7 s$y#Ufz
TBDATA,1,4.5E5,24.8 !iz vY
!取强度折减系数F=2.2时边坡围岩材料属性 g1 !L.
On
MP,EX,8,3E10 xTawG?"D
MP,PRXY,8,0.25 M>z7H"jCu
MP,DENS,8,2500 o %J IJ7M
TB,DP,8 !定义数据表8 V$F.`O!hfi
TBDATA,1,4.09E5,22.8 5$G??="K
!取强度折减系数F=2.4时边坡围岩材料属性 e ~,'|~
C5
MP,EX,9,3E10 g4qdm{BL
MP,PRXY,9,0.25 #3[b|cL
MP,DENS,9,2500 x r-`i
TB,DP,9 !定义数据表9 EYD{8Fw-
TBDATA,1,3.6E5,21.1 !给数据表输入数据 ML"P"&~u6
!取强度折减系数F=2.6时边坡围岩材料属性 4 AmF ^H
MP,EX,10,3E10 ！定义弹性模量 O)Dw<j)
MP,PRXY,10,0.25 !定义泊松比 R3! vS+5rR
MP,DENS,10,2500 ！定义密度 G(3 la3\(
TB,DP,10 !定义数据表10 w3cK:
C 0
TBDATA,1,3.46E5,19.6 !给数据表10输入数据 tNf" X !
!取强度折减系数F=2.8时边坡围岩材料属性 KIl. ?_61O
MP,EX,11,3E10 ！定义弹性模量 Z!l!3(<G.f
MP,PRXY,11,0.25 !定义泊松比 Q>c6ouuJ
MP,DENS,11,2500 1E_ Ui1 [
TB,DP,11 !定义数据表11 Qi]Z)v{ ^
TBDATA,1,3.2E5,18.2 !给数据表11输入数据 *xVAm 7_v
!取强度折减系数F=3.0时边坡围岩材料属性 s ^@Cq=
MP,EX,12,3E10 ！定义弹性模量 ]TprPU39
MP,PRXY,12,0.25 *K^O oS
MP,DENS,12,2500 ！定义密度 l5nDt$Ex
TB,DP,12 !定义数据表12 O,Gn2Do
TBDATA,1,3.0E5,17.1 !给数据表12输入数据 ?
Z
fhz
!建立几何模型 R2ZQBwB
!创建边坡线模型线 2 3 w{h d
K,1,,,, ！ 创建关键点1 0BD((oNg
K,2,-800,,, O}!@28|3"
K,3,-800,-800,, 0lh6b3tdP
K,4,-800,-1200,, 1)r _h(
K,5,1200,-1200,, :R x"WY
K,6,1200,-800,, O +}EE^*a
K,7,1200,,, fzG1<Gem
K,8,1200,400,, 8J U~Q
K,9,400,400,, ov>L-
l,1,2 !连接关键点1、2生成直线L1 M.^A`
l,2,3 zm"g,\.d
l,3,4 s""8V_,;
l,4,5 2/N*Uk 0
l,5,6 fsuvg jlE
l,6,7 .g?,:$`0D?
l,7,8 [
[CXMbD`*
l,8,9 ]arskmB]
l,9,1 S ylsp%A
l,1,7 Y &+/[ [
l,3,6 N0hU~| /
!创建边坡面模型 ?tQv|x
al,3,4,5,11 !由线L3、L4、L5、L11生成一个面A1 A6.'1OD
al,1,2,11,6,10 !由线L1、L2、L1、L6、L10生成一个面A2 _oVA0@#n
al,7,8,9,10 !由线L7、L8、L9、L10生成一个面A3 74Wg@! P
/PNUM, area, on !打开面号开关 3mz>Y*^?0
!划分网格生成有限元模型 ` :O.g9
! 划分边坡围岩2网格 M X\-)e#
!设置网格份数 S6X<3L`FfH
lsel,s,,,3,5,2 !选择线L3和L5 )KQ um`pO
lesize,all,,,5 !把所选择线分为5段 a[ l5k
lsel,s,,,4,11,7 !选择线L4和L11 XnOl*#P
lesize,all,,,26 !把所选择线分为26段 qEz' l'%(
mat,2 !给边坡围岩2赋予2号材料特性 TvwIro
type,1 !采用单元类型1 H E '8
mshkey,0 !设定自由网格划分 ^k##a-t<_>
mshape,0 !设定四边形网格划分 1 `\kXaG
amesh,1 ！划分面积A1 z 59J=?|
! 划分边坡围岩1网格 ^SCWT\E
!设置网格份数 ~C
x2Q4E
lsel,s,,,2,6,4 !选择线L2和L6 qN L~m'
lsel,a,,,1 !附加选择线L1 ^u,x~nPXg
lesize,all,,,10 !把所选择线分为10段 O r>[_3
lsel,s,,,8,10,2 !选择线L8和L10 hZJqo + s
lesize,all,,,16 !把所选择线分为16段 \M@8# k|
lsel,s,,,7,9,2 f'U]Ik;Jy
lesize,all,,,12 t Ya8I/HpT
mat,1 !给边坡围岩1赋予1号材料特性 eO G%6C%a
type,1 WQL`;uIX
mshkey,0 !设定自由网格划分 &X(-C9'j
mshape,0 !设定四边形网格划分 &N ;6G`3
amesh,2 !划分面积A2 #9Z-Hd<
amesh,3 ！划分面积A3 %L.+r!.
allsel TD7ON a-,
!保存网格模型 a-y5 \x
save,Slope-grid.db UkV] F]
(3)施加约束和荷载
(5_(s`q.
!给边坡模型施加约束 )_?$B6hf,&
!边坡两侧施加X方向约束 .`].\Zyk f
nsel,s,loc,x,-800 !选择X=-800线上所有节点 :g";p.~=
nsel,a,loc,x,1200 !选择X=1200线上所有节点 &pz8v WCk
d,all,ux !对所选择节点约束X方向位移 ZeZwzH)B D
allsel _jz=BRO$
!边坡底部施加约束 ?2;G_P+
nsel,s,loc,y,-1200 !选择Y=--1200线上所有节点 )' #(1
,1k
d,all,uy !对所选择节点约束X、Y方向位移 _3S{n=9
d,all,ux 1 Y& d%AA
!施加重力加速度 P0=F9`3 wb
acel,,9.8 kpu^:N &
(4)求解 ,L G&sa"
/solu cj
*4 XYu
!求解设置 LAeX e!y
antype,static !设定为静力求解 Q<pM
tW
nsubst,100 !设定最大子步数为100 X;NTz75
pred,on !打开时间步长预测器 ew g WzB9c
nropt,full !设定牛顿-拉普森选项 GZo4uwG@a
nlgeom,on !打开大位移效果 !A o?bs'
lnsrch,on !打开线性搜索 ?da 3Azp
outres,all,all !输出所有项 $Vzfhj-if
cnvtol,f,,0.005,2,0.5 !力收敛准则设定 ]u ';zJ.
cnvtol,u,,0.05,2,1 !位移收敛准则设定 ,+&j/0U
!边坡在强度折减系数F=1时求解 D~:fn|/Brp
allsel rGWTp N
solve !进行求解 _?ZT[t<
save,F1,db !把F=1时求解结果保存 U<b!$"P9
!边坡在强度折减系数F=1.2时求解 9':$!Eoq
finish A-FwNo2"%
/solu n C^'2z
allsel 2 O TpGl
asel,s,area,,2,3,1 !选择面积A2、A3 jo~Pr
MPCHG,3,all ! 把所选择单元材料号改为3 d}Om?kn
allsel MR}\fw$(.
solve !进行求解 rL=$Wxd PU
save,F1.2,db !把F=1.2时求解结果保存 FcA0 \`0M
!边坡在强度折减系数F=1.4时求解 V|)>
finish V_* ^2c)
/solu @w0[5ZA j
allsel jVz1`\Nje
asel,s,area,,2,3,1 !选择面积A2、A3 6D],275`J
MPCHG,4,all !把所选择单元材料号改为4 2wd(0K}b
allsel v6:DA#0
solve !进行求解 he Wb(E&
save,F1.4,db !把F=1.4时求解结果保存 ,n* .Yq
!边坡在强度折减系数F=1.6时求解 MYDSkW
finish Tx5L
/solu Gm. hBNgp
allsel C6n4OU
asel,s,area,,2,3,1 !选择面积A2、A3 EB#z \
MPCHG,5,all !把所选择单元材料号改为5 5IMH G%W7
allsel d,77 L
solve P*SXfb"HC
save,F1.6,db !把F=1.6时求解结果保存 A SSoKrFL
!边坡在强度折减系数F=1.8时求解 s@$AYZm_
finish c nzPq\
/solu }/VHeHd
allsel Ij+zR>P8=\
asel,s,area,,2,3,1 !选择面积A2、A3 pqe**`z@y
MPCHG,6,all !把所选择单元材料号改为6 fHigLL0B
allsel luMNi^FQ
solve md0=6<
}P
save,F1.8,db !把F=1.8时求解结果保存 AS7!FD6b
!边坡在强度折减系数F=2.0时求解 xn"g_2Hi
finish f As:[
/solu 1} 6pq 2
allsel j !s& yHE1
asel,s,area,,2,3,1 !选择面积A2、A3 ? _W*7<
MPCHG,7,all !把所选择单元材料号改为7
)nY/ RO
allsel W.a/k7 p
solve "#7 i-?=
save,F2.0,db !把F=2.0时求解结果保存 CqoL5qt
!边坡在强度折减系数F=2.2时求解 g _=Q=y@,
finish C$ bK!]a
/solu DB0xIP~i,?
allsel `G'V9Xs(
asel,s,area,,2,3,1 !选择面积A2、A3 Bg`b*(Q
MPCHG,8,all !把所选择单元材料号改为8 zE [c$KPP
allsel SCjACQ}-
solve *M"wH_cd
save,F2.2,db !把F=2.2时求解结果保存 rnr7t \a~]
!边坡在强度折减系数F=2.4时求解 h2q]!01XP
finish MiC&av
/solu s+[=nau('w
allsel d>%_<pw
asel,s,area,,2,3,1 !选择面积A2、A3 & zG=
MPCHG,9,all !把所选择单元材料号改为9 E|>I/!{u7`
allsel P
,i)A
solve U0rz 4fxc
save,F2.4,db !把F=2.4时求解结果保存 en
!边坡在强度折减系数F=2.6时求解 te|VKYN%}[
finish 8{ep`$(K@
/solu {\V)bizY;
allsel 8[DD=[&
asel,s,area,,2,3,1 !选择面积A2、A3 lY 1m%
MPCHG,10,all !把所选择单元材料号改为10 / nr DU*
allsel B0Z@ Cf
solve _ehU:3L`s
save,F2.6,db !把F=2.6时求解结果保存 |9y &;3
!边坡在强度折减系数F=2.8时求解
pK3cg|}
finish -X=f+4 j
/solu ~DJ/sY2/
allsel DgdW.Kj|IL
asel,s,area,,2,3,1 !选择面积A2、A3 <W?,n%
MPCHG,11,all !把所选择单元材料号改为11 78X;ZMY
allsel +q4T] ;<
solve 65AOFH
save,F2.8,db !把F=2.8时求解结果保存 4uz\Me(
!边坡在强度折减系数F=3.0时求解 "-hgeQX
finish }$7Hf+G
/solu ]8mBFr5E9
allsel `jV0;sPd;
asel,s,area,,2,3,1 !选择面积A2、A3 [p <L*3<
MPCHG,12,all !把所选择单元材料号改为12 $] Q_x?
allsel 8\yH 7H
solve PC9:n ee
save,F3.0,db !把F=3.0时求解结果保存 w qsPGkJJ7
(5）后处理 I@hC$o
/post1 ！进入后处理 oreS u;`$
!边坡在强度折减系数F=1时结果分析 9Kqr9U--v
Resume,'F1','db' !读入边坡在强度折减系数F=1时 q)P<lKi
set,1,last !读入后一个子步 %PG0PH4?
pldisp,1 !绘制边坡模型变形图 mb&b = &
plnsol,u,x !绘制边坡模型水平方向位移云图 OxD\e5r
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 mU3UQ
j
!边坡在强度折减系数F=1.2时结果分析 pG"
4qw
Resume,'F1.2','db' !读入边坡在强度折减系数F=1.2时 H[Qh* pq2
set,1,last !读入后一个子步 >uQ!B/C!
pldisp,1 !绘制边坡模型变形图 BjB2YO& /
plnsol,u,x !绘制边坡模型水平方向位移云图 eSvu:euv
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 , vky
!边坡在强度折减系数F=1.4时结果分析 H
fRxgA @
Resume,'F1.4','db' !读入边坡在强度折减系数F=1.4时 &o?pZ(\C
set,1,last !读入后一个子步 _-D(N/
pldisp,1 !绘制边坡模型变形图 5[GX
plnsol,u,x !绘制边坡模型水平方向位移云图 A1:Fe9q
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 t'z] <7
!边坡在强度折减系数F=1.6时结果分析 M0 ^r!f>O
Resume,'F1.6','db' !读入边坡在强度折减系数F=1.6时 M~@\x]p >
set,1,last !读入后一个子步 uZqo"
pldisp,1 !绘制边坡模型变形图 >6aCBS? 2
plnsol,u,x !绘制边坡模型水平方向位移云图 xUG|@xIwc
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 X=DJOepH'
!边坡在强度折减系数F=1.8时结果分析 vbSz&+52;
Resume,'F1.8','db' !读入边坡在强度折减系数F=1.8时 M HgS5 b2
set,1,last !读入后一个子步 wi9fYfuv3R
pldisp,1 !绘制边坡模型变形图 k_!z=6?[:
plnsol,u,x !绘制边坡模型水平方向位移云图 YKk%lZ.8
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 mX<D]Z< k
!边坡在强度折减系数F=2.0时结果分析 :?60pu=
Resume,'F2.0','db' !读入边坡在强度折减系数F=2.0时 >s1HQSe66
set,1,last !读入后一个子步 wngxVhu8Ld
pldisp,1 !绘制边坡模型变形图 !x,3k\M
plnsol,u,x !绘制边坡模型水平方向位移云图 #|'8 O
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 % D&FnTa
!边坡在强度折减系数F=2.2时结果分析 :}E*u^v K
Resume,'F2.2','db' !读入边坡在强度折减系数F=2.2时 :\.v\.wm
set,1,last !读入后一个子步 s_RYYaM
pldisp,1 !绘制边坡模型变形图 mhcJ0\@_
plnsol,u,x !绘制边坡模型水平方向位移云图 ag] nVE/
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 wv1?v_4
!边坡在强度折减系数F=2.4时结果分析 a5&wS@)
;
Resume,'F2.4','db' !读入边坡在强度折减系数F=2.4时 SBYRN##n_
set,1,last !读入后一个子步 iA%3cpIc(Z
pldisp,1 !绘制边坡模型变形图 q6osRK*20
plnsol,u,x !绘制边坡模型水平方向位移云图 I{ M2nQi
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 \0b ",|"3
!边坡在强度折减系数F=2.6时结果分析 ?Ww ',e
Resume,'F2.6','db' !读入边坡在强度折减系数F=2.6时 ma-|L3 #
set,1,last !读入后一个子步 S;CT:kG6Y{
pldisp,1 !绘制边坡模型变形图 h>\}- |Ek
plnsol,u,x !绘制边坡模型水平方向位移云图 T sR20P @
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 ,PY<AI^59
!边坡在强度折减系数F=2.8时结果分析 Qi[T!1
Resume,'F2.8','db' !读入边坡在强度折减系数F=2.8时 6_9w1
,W E
set,1,last !读入后一个子步 JhuK W>7
pldisp,1 !绘制边坡模型变形图 S?D2`b
plnsol,u,x !绘制边坡模型水平方向位移云图 %,>,J`
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 $bMeL7CN
!边坡在强度折减系数F=3.0时结果分析 uz*C`T0:rj
Resume,'F3.0','db' !读入边坡在强度折减系数F=3.0时 *aTM3k)Zs
set,1,last !读入后一个子步 d *H-l3N
pldisp,1 !绘制边坡模型变形图 |px 4a"
plnsol,u,x !绘制边坡模型水平方向位移云图 2*Gl|@~N
plnsol,eppl,eqv !绘制边坡模型塑性应变云图
