各位朋友： Sn[/'V^$a
大家好！我通过ansys用有限元强度折减法分析了一个边坡的稳定性，结果不同的工况下得出的塑性区是一样的，我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激！附命令流和塑性区如下：(1)创建物理环境 a.fdCI]%
/COM, Structural ! 指定结构分析 YhL^kM@c
/TITLE, Slope stability Analysis ! 定义工作标题 k;jXVa
/FILNAM,support,1 ! 定义工作文件名 P
4+}<5
(2)建立模型 Koc5~qUY]
!进入前处理器 $Ao'mT
/PREP7 Oe lf^&m
!定义单元类型 Zhc99 L &K
ET,1,PLANE82 !定义边坡围岩单元 t@q==VHF
KEYOPT,1,3,2 !求解类型为平面应变 aq - |
!定义材料属性 ?vQ:z{BO
!边坡围岩1材料属性 #q6#nfi"
MP,EX,1,3E10 eZ{Ce.lNR
MP,PRXY,1,0.25 k \ \e`=
MP,DENS,1,2500 !定义塑性模型 -!IeP]n#P
TB,DP,1 !定义数据表1 Dw<bLSaW&
TBDATA,1,9E5,42.8 !给数据表输入数据 )oy+-1dE
!边坡围岩2材料属性 !采用弹性模型 ff
6x4t
MP,EX,2,3.2E10 ！定义弹性模量 SZ~lCdWad
MP,PRXY,2,0.24 !定义泊松比 eQ<Vky^SJ
MP,DENS,2,2700 !定义密度 ^;wz+u4^l
!取强度折减系数F=1.2时边坡围岩材料属性 ': ]Hj8t_
MP,EX,3,3E10 t\f[ ->f
MP,PRXY,3,0.25 Av! xI
MP,DENS,3,2500 wxy@XN"/i+
TB,DP,3 !定义数据表3 EF'8-*
TBDATA,1,7.5E5,37.7 !给数据表输入数据 vK$wc ~
!取强度折减系数F=1.4时边坡围岩材料属性 L9d|7.b
MP,EX,4,3E10 A+(+Pf U
MP,PRXY,4,0.25 !定义泊松比 =ba1::18
MP,DENS,4,2500 mQ2=t %
TB,DP,4 !定义数据表4 ubM N
TBDATA,1,6.4E5,33.5 !给数据表输入数据 ?*'0;K13
!取强度折减系数F=1.6时边坡围岩材料属性 -#,4rN#
MP,EX,5,3E10 v1=N?8Hz1
MP,PRXY,5,0.25 sW76RKX8
MP,DENS,5,2500 ！定义密度 oj@=Cq':-
TB,DP,5 !定义数据表5 ,Jf P$HJ
TBDATA,1,5.6E5,30 !给数据表输入数据 yYdh+ x
!取强度折减系数F=1.8时边坡围岩材料属性 +3Z+#nGtk
MP,EX,6,3E10 ！定义弹性模量 8\p"V.o>
MP,PRXY,6,0.25 v
,zD52
MP,DENS,6,2500 JA4}B wn
TB,DP,6 !定义数据表6 X2'XbG 3
TBDATA,1,5.0E5,27.2 9U Hh#
!取强度折减系数F=2.0时边坡围岩材料属性 <,Mf[R2N>
MP,EX,7,3E10 l[[^]__
MP,PRXY,7,0.25 !定义泊松比 i!iG7X)qT
MP,DENS,7,2500 3gtKD9RL:
TB,DP,7 !定义数据表7 ll*Ez"
TBDATA,1,4.5E5,24.8 V>c !V9w
!取强度折减系数F=2.2时边坡围岩材料属性 HhwAzk/G~
MP,EX,8,3E10 1'|gxYT
MP,PRXY,8,0.25 B qLL]%F
MP,DENS,8,2500 *:ErZ UyQM
TB,DP,8 !定义数据表8 2Sk"S/4}Z
TBDATA,1,4.09E5,22.8 j1=su~
!取强度折减系数F=2.4时边坡围岩材料属性 35 L\
MP,EX,9,3E10 yT[=!M
MP,PRXY,9,0.25
I]BhkJ
MP,DENS,9,2500 t:b}Mo0
TB,DP,9 !定义数据表9 ^fmuBe}d{
TBDATA,1,3.6E5,21.1 !给数据表输入数据 N?O^"
!取强度折减系数F=2.6时边坡围岩材料属性 &"7+k5O
MP,EX,10,3E10 ！定义弹性模量 N1"p ;czK
MP,PRXY,10,0.25 !定义泊松比 1fC)&4W
MP,DENS,10,2500 ！定义密度 POf xN.
TB,DP,10 !定义数据表10 a)I>Ns)
TBDATA,1,3.46E5,19.6 !给数据表10输入数据 )TEod!]
!取强度折减系数F=2.8时边坡围岩材料属性 bz.sWBugR
MP,EX,11,3E10 ！定义弹性模量 +.zX ?}
MP,PRXY,11,0.25 !定义泊松比 p{+F{e
MP,DENS,11,2500 >=;hnLu
TB,DP,11 !定义数据表11 9FPqd8(]*V
TBDATA,1,3.2E5,18.2 !给数据表11输入数据 6C\WX(@4
!取强度折减系数F=3.0时边坡围岩材料属性 [z!pm-Ir
MP,EX,12,3E10 ！定义弹性模量 9Ais)Wy%p
MP,PRXY,12,0.25 BCtKxtbS
MP,DENS,12,2500 ！定义密度 4 p%^?L?
TB,DP,12 !定义数据表12 rV} 5&N*c
TBDATA,1,3.0E5,17.1 !给数据表12输入数据 oTri t_@3
!建立几何模型 D
!{e
!创建边坡线模型线 CeM %?fr5
K,1,,,, ！ 创建关键点1 }pGj c_:']
K,2,-800,,, 5rmQ:8_5
K,3,-800,-800,, r! [Qpb-:
K,4,-800,-1200,, /:"^,i\t
K,5,1200,-1200,, 5s NN:m
K,6,1200,-800,, i :jns>E
K,7,1200,,, [f}`reRlZ
K,8,1200,400,, pl$wy}W-
K,9,400,400,, mq (-L
l,1,2 !连接关键点1、2生成直线L1 <[ u(il
l,2,3 & eqqgLz
l,3,4 piqh7u3~
l,4,5 1.TIUH1
l,5,6 T# i~/
l,6,7 '-cayG
l,7,8 cI/ }r Z+
l,8,9
HLQ>
|,9
l,9,1 I!SIy&=W
l,1,7 reM~q-M~o@
l,3,6 !;P[Y"h@r
!创建边坡面模型 | g1Cs
al,3,4,5,11 !由线L3、L4、L5、L11生成一个面A1 %@|)&][hO
al,1,2,11,6,10 !由线L1、L2、L1、L6、L10生成一个面A2 2
/*z5
al,7,8,9,10 !由线L7、L8、L9、L10生成一个面A3 %LD(S * >7
/PNUM, area, on !打开面号开关 ^bfU>02Q6p
!划分网格生成有限元模型 H328I}7
! 划分边坡围岩2网格 y&bZai8WlE
!设置网格份数
V<?0(esgR
lsel,s,,,3,5,2 !选择线L3和L5 2%9 L'-
lesize,all,,,5 !把所选择线分为5段 /Ak\Q5O'3
lsel,s,,,4,11,7 !选择线L4和L11 ,H8M.hbsQ
lesize,all,,,26 !把所选择线分为26段 H9(UzyN>i
mat,2 !给边坡围岩2赋予2号材料特性 4bL? V^@7
type,1 !采用单元类型1 dB S_N/
mshkey,0 !设定自由网格划分 GG-b)64h`
mshape,0 !设定四边形网格划分 y u &mu CA
amesh,1 ！划分面积A1 W\ mgM2p
! 划分边坡围岩1网格 d|)ARRW
!设置网格份数 X*w;6 V
lsel,s,,,2,6,4 !选择线L2和L6 ] w0Y5H "
lsel,a,,,1 !附加选择线L1 =e0MEV#s.
lesize,all,,,10 !把所选择线分为10段 J<4 _<.o(a
lsel,s,,,8,10,2 !选择线L8和L10 b5I 8jPj4c
lesize,all,,,16 !把所选择线分为16段 s@G E(Pu7
lsel,s,,,7,9,2 ~%eE%5!k
lesize,all,,,12 w4_Xby)
mat,1 !给边坡围岩1赋予1号材料特性 7=Pj}x)
type,1 BUV4L5(
mshkey,0 !设定自由网格划分 3<N2 ehi?
mshape,0 !设定四边形网格划分 DY{v@
<3
amesh,2 !划分面积A2 dLeos9M:
amesh,3 ！划分面积A3 m,J
IId%O
allsel lQ4^I^?m
!保存网格模型 _n{_\/A6 f
save,Slope-grid.db ^)a:D KL
(3)施加约束和荷载 R y (<6u0
!给边坡模型施加约束 Fy!s$!\C0
!边坡两侧施加X方向约束 +nim47
nsel,s,loc,x,-800 !选择X=-800线上所有节点 _[i.)8$7
nsel,a,loc,x,1200 !选择X=1200线上所有节点 U2tgBF?)A
d,all,ux !对所选择节点约束X方向位移 uEG4^
allsel QMmZvz\^
!边坡底部施加约束 uFhPNR2l
nsel,s,loc,y,-1200 !选择Y=--1200线上所有节点 1 ?3+>
d,all,uy !对所选择节点约束X、Y方向位移 wZ69W$, p
d,all,ux o|c&$)m
!施加重力加速度 .IVKgQ
B
acel,,9.8 !q$>6P
(4)求解 %++S;#)~
/solu Uovna:"
!求解设置 b'` XFB#V
antype,static !设定为静力求解 y 4aT-^C'
nsubst,100 !设定最大子步数为100 E>gLUMG$
pred,on !打开时间步长预测器 /~7M @`1
nropt,full !设定牛顿-拉普森选项 UV av^<_
nlgeom,on !打开大位移效果 Ag*?>I
lnsrch,on !打开线性搜索 |9Yx`_DF
outres,all,all !输出所有项 \
bWy5/+
cnvtol,f,,0.005,2,0.5 !力收敛准则设定 m=v.< +>
cnvtol,u,,0.05,2,1 !位移收敛准则设定 dM19;R@4
!边坡在强度折减系数F=1时求解 +Z GH
allsel mA_EvzXk\
solve !进行求解 /v| b]Ji
save,F1,db !把F=1时求解结果保存 ;=E}Pb Zt2
!边坡在强度折减系数F=1.2时求解 5|t-CY{?b
finish * C Az_s<
/solu R/Sm
allsel wV{j CQ
asel,s,area,,2,3,1 !选择面积A2、A3 p]?eIovi
MPCHG,3,all ! 把所选择单元材料号改为3 gp~-n7'~O
allsel B|+% ExT7
solve !进行求解
j0O1??
save,F1.2,db !把F=1.2时求解结果保存 -*Qg^ 1]i+
!边坡在强度折减系数F=1.4时求解 'O9Yu{M
finish /R?*i@rvf
/solu 45iO2W uur
allsel h.Sb ds
asel,s,area,,2,3,1 !选择面积A2、A3 \ $t{K
MPCHG,4,all !把所选择单元材料号改为4 9$ VdYw7D
allsel - em3 #V
solve !进行求解 b
j<T`M!
save,F1.4,db !把F=1.4时求解结果保存 p)"EenUK
!边坡在强度折减系数F=1.6时求解 .L^;aL
finish iEy2z+/"^
/solu #)#'^MZX
allsel v<g#/X8
asel,s,area,,2,3,1 !选择面积A2、A3 R`@7f$;wG
MPCHG,5,all !把所选择单元材料号改为5 jv1p'qs4
allsel &9.3-E47*
solve aO9\8\^
save,F1.6,db !把F=1.6时求解结果保存 C:bA:O
!边坡在强度折减系数F=1.8时求解 [xDIK8d:I
finish S"87 <o
/solu ;i+(Q%LO
allsel 0t MzVx S
asel,s,area,,2,3,1 !选择面积A2、A3 VeQ [A?pER
MPCHG,6,all !把所选择单元材料号改为6 # -R]HLW*
allsel 36.mf_AM
solve >SYOtzg%
save,F1.8,db !把F=1.8时求解结果保存 I<xcVY9L
!边坡在强度折减系数F=2.0时求解 !VrBoU 4<d
finish Qt^6w}&
/solu w> Ft5"z
allsel Aqg$q* Y
asel,s,area,,2,3,1 !选择面积A2、A3 |BGzdBm^x:
MPCHG,7,all !把所选择单元材料号改为7 `$3P@SO"
allsel ,pkzNe` F
solve @ e7_&EGR?
save,F2.0,db !把F=2.0时求解结果保存 R\$6_
!边坡在强度折减系数F=2.2时求解 5VfpeA `
finish _nw\ac#*
/solu |Z Cv>8 ?n
allsel x0TE+rf5
asel,s,area,,2,3,1 !选择面积A2、A3 QEx&AT
MPCHG,8,all !把所选择单元材料号改为8 E !I4I'
allsel I4c %>R
solve }2Im?Q
save,F2.2,db !把F=2.2时求解结果保存 +yHzp
!边坡在强度折减系数F=2.4时求解 CyB 1`&G>
finish Rob: W|
/solu kaDn=
={YM
allsel
Ox'K C
asel,s,area,,2,3,1 !选择面积A2、A3 bs?4|#[K
MPCHG,9,all !把所选择单元材料号改为9 u9+)jN<Yh
allsel )S_ %Ip
solve R/^;,.
save,F2.4,db !把F=2.4时求解结果保存 \D[~ 54
!边坡在强度折减系数F=2.6时求解 z:{R4#(Q
finish -**fT?n
/solu ? C6`
allsel 8R;E+B{
asel,s,area,,2,3,1 !选择面积A2、A3 vha@YPC=
MPCHG,10,all !把所选择单元材料号改为10 |txzIc.#
allsel \ nQEvcH
solve )9!ZkZbv_m
save,F2.6,db !把F=2.6时求解结果保存 *I9G" R8
!边坡在强度折减系数F=2.8时求解 a1weTn*
finish !}xRwkN
/solu CR|>?9V
allsel D ,uT#P
asel,s,area,,2,3,1 !选择面积A2、A3 gti=GmL(L
MPCHG,11,all !把所选择单元材料号改为11 s: MJ{r(s
allsel _'s5FlZq
solve $,Y?q n/
save,F2.8,db !把F=2.8时求解结果保存 81wmKqDEs
!边坡在强度折减系数F=3.0时求解 IRI<no
finish + OUM 4y
/solu Gs}lw'pK
allsel a & R,jq
asel,s,area,,2,3,1 !选择面积A2、A3 D/WzYc2h]
MPCHG,12,all !把所选择单元材料号改为12 9Mv4=k^7|4
allsel $I / RN
solve jH4Wu`r;m
save,F3.0,db !把F=3.0时求解结果保存 Szlww
(5）后处理 )v.\4Q4
/post1 ！进入后处理 @{880 5D p
!边坡在强度折减系数F=1时结果分析 0+S'i82=M
Resume,'F1','db' !读入边坡在强度折减系数F=1时 y? 65*lUl
set,1,last !读入后一个子步 V0:d b
pldisp,1 !绘制边坡模型变形图 c<qJs-C4;
plnsol,u,x !绘制边坡模型水平方向位移云图 c-a;n AR
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 9,
792b
!边坡在强度折减系数F=1.2时结果分析 vYG$>*
Resume,'F1.2','db' !读入边坡在强度折减系数F=1.2时 \G1(r=fU
set,1,last !读入后一个子步 *K w/ilI
pldisp,1 !绘制边坡模型变形图 d _koF-7
plnsol,u,x !绘制边坡模型水平方向位移云图 faI4`.i
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 N;e*eMFE
!边坡在强度折减系数F=1.4时结果分析 sk`RaDq@;
Resume,'F1.4','db' !读入边坡在强度折减系数F=1.4时 VnMiZAHR
set,1,last !读入后一个子步 N" oJ3-~
pldisp,1 !绘制边坡模型变形图 k+cHx799
plnsol,u,x !绘制边坡模型水平方向位移云图 <4Cy U
j
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 V{kgD pB
!边坡在强度折减系数F=1.6时结果分析 knZ<V%/e
Resume,'F1.6','db' !读入边坡在强度折减系数F=1.6时 /ec~ ^S8X
set,1,last !读入后一个子步 y]f| U-f:~
pldisp,1 !绘制边坡模型变形图 Ad`jV_z
plnsol,u,x !绘制边坡模型水平方向位移云图 3Kum
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 ^kp u9H
!边坡在强度折减系数F=1.8时结果分析 w+ tO@
Resume,'F1.8','db' !读入边坡在强度折减系数F=1.8时 HnfTj 5J@
set,1,last !读入后一个子步 MMU w+jM4
pldisp,1 !绘制边坡模型变形图 m dHC{sp
plnsol,u,x !绘制边坡模型水平方向位移云图 n? A;'\cK
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 Mxe
!边坡在强度折减系数F=2.0时结果分析 uM~j
Resume,'F2.0','db' !读入边坡在强度折减系数F=2.0时
/=7[Q
set,1,last !读入后一个子步 gG=E2+=uy
pldisp,1 !绘制边坡模型变形图 meV
RdQ
plnsol,u,x !绘制边坡模型水平方向位移云图 aG\B?pn-
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 Z@`HFZJ
!边坡在强度折减系数F=2.2时结果分析 x;R9Gc[5
Resume,'F2.2','db' !读入边坡在强度折减系数F=2.2时 zHCz[jlrMq
set,1,last !读入后一个子步 -vC?bumR%
pldisp,1 !绘制边坡模型变形图 1e^-_Bo6'o
plnsol,u,x !绘制边坡模型水平方向位移云图 [t`QV2um
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 2]*2b{gF,
!边坡在强度折减系数F=2.4时结果分析 c
=i 6
Resume,'F2.4','db' !读入边坡在强度折减系数F=2.4时 =
8%+$vX
set,1,last !读入后一个子步 x_5H_! \#
pldisp,1 !绘制边坡模型变形图 vd <"
G}
plnsol,u,x !绘制边坡模型水平方向位移云图 #VEHyz 6P
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 }mC-SC)oSi
!边坡在强度折减系数F=2.6时结果分析 - gV'z5
Resume,'F2.6','db' !读入边坡在强度折减系数F=2.6时 P1ab2D
set,1,last !读入后一个子步 4m6E~_:F
pldisp,1 !绘制边坡模型变形图 < tg>1,C
plnsol,u,x !绘制边坡模型水平方向位移云图 3J}bI {3
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 j7 D\O
!边坡在强度折减系数F=2.8时结果分析 +d736lLe%
Resume,'F2.8','db' !读入边坡在强度折减系数F=2.8时 f"5vpU^5*
set,1,last !读入后一个子步 H;$O CDRC
pldisp,1 !绘制边坡模型变形图 &u_f:Pog
plnsol,u,x !绘制边坡模型水平方向位移云图 lU.@! rGbw
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 iB 5 Se
!边坡在强度折减系数F=3.0时结果分析 o3l _&?^
Resume,'F3.0','db' !读入边坡在强度折减系数F=3.0时 U.G** v
set,1,last !读入后一个子步 6 l>$N?a
pldisp,1 !绘制边坡模型变形图 j{9sn,<:
plnsol,u,x !绘制边坡模型水平方向位移云图 @vaK-&|#$
plnsol,eppl,eqv !绘制边坡模型塑性应变云图
