各位朋友： :jCaDhK
大家好！我通过ansys用有限元强度折减法分析了一个边坡的稳定性，结果不同的工况下得出的塑性区是一样的，我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激！附命令流和塑性区如下：(1)创建物理环境 KPr xw }P
/COM, Structural ! 指定结构分析 l$@lk?dc
/TITLE, Slope stability Analysis ! 定义工作标题 zL}DL fy>R
/FILNAM,support,1 ! 定义工作文件名 :O=Vr]Y8K
(2)建立模型 7hE=+V8
!进入前处理器 Jk{2!uP
/PREP7 } z:=b8}
!定义单元类型 NV ~i4R*#
ET,1,PLANE82 !定义边坡围岩单元
<T9m.:l
KEYOPT,1,3,2 !求解类型为平面应变 i I RigW
!定义材料属性 &?pAt30K:
!边坡围岩1材料属性 b8Sl3F?-~
MP,EX,1,3E10 ~Dh}E9E:
MP,PRXY,1,0.25 w N.Jy b
MP,DENS,1,2500 !定义塑性模型 x=+H@YO\
TB,DP,1 !定义数据表1 <<6#Uz.1
TBDATA,1,9E5,42.8 !给数据表输入数据 WJ,ON-v
!边坡围岩2材料属性 !采用弹性模型 lsf?R'1
MP,EX,2,3.2E10 ！定义弹性模量 )^3655mb
MP,PRXY,2,0.24 !定义泊松比 _ 3JTHf<+
MP,DENS,2,2700 !定义密度 G9%4d;uFT
!取强度折减系数F=1.2时边坡围岩材料属性 ZmF32 Ir
MP,EX,3,3E10 / d;l:
MP,PRXY,3,0.25 Tl-B[CT
MP,DENS,3,2500 cV,URUD
TB,DP,3 !定义数据表3 j(N9%/4u
TBDATA,1,7.5E5,37.7 !给数据表输入数据 }7^*%$
!取强度折减系数F=1.4时边坡围岩材料属性 JE!Xf}nEi
MP,EX,4,3E10 ~<-h# B
MP,PRXY,4,0.25 !定义泊松比 cFo-NI2
MP,DENS,4,2500 KO7cZME
TB,DP,4 !定义数据表4 o^J&c_U\3'
TBDATA,1,6.4E5,33.5 !给数据表输入数据 eGg #=l=
!取强度折减系数F=1.6时边坡围岩材料属性 1Tkz!
MP,EX,5,3E10 3iwZUqyq
MP,PRXY,5,0.25 ObnB6ShKi
MP,DENS,5,2500 ！定义密度 j'cS_R
TB,DP,5 !定义数据表5 M1\/ueOe
TBDATA,1,5.6E5,30 !给数据表输入数据 0}!lN{m?
!取强度折减系数F=1.8时边坡围岩材料属性 ^hJ ,1{o
MP,EX,6,3E10 ！定义弹性模量 Gd 9B
MP,PRXY,6,0.25 /qr8
MP,DENS,6,2500 UNHHzTsr?
TB,DP,6 !定义数据表6 ponvi42u
TBDATA,1,5.0E5,27.2 n_Dhq (.
!取强度折减系数F=2.0时边坡围岩材料属性 r 6<}S(
MP,EX,7,3E10 \U8Vsx1tl
MP,PRXY,7,0.25 !定义泊松比 sIe(;%[`
MP,DENS,7,2500 \ Ld7fP
TB,DP,7 !定义数据表7 C7:Ry)8'I
TBDATA,1,4.5E5,24.8 pj`-T"Q
!取强度折减系数F=2.2时边坡围岩材料属性 '-_PO|}
MP,EX,8,3E10 nz+KA\iW
MP,PRXY,8,0.25 wq#3f# 3V
MP,DENS,8,2500 9 R1]2U $|
TB,DP,8 !定义数据表8 otx7J \4
TBDATA,1,4.09E5,22.8 e)8iPu ..
!取强度折减系数F=2.4时边坡围岩材料属性 I_Qnq4Sk(
MP,EX,9,3E10 ,GlK_-6>
MP,PRXY,9,0.25 -W ,b*U
MP,DENS,9,2500 1-fz564
TB,DP,9 !定义数据表9 oB@ C-(M
TBDATA,1,3.6E5,21.1 !给数据表输入数据 C_/o ORvK
!取强度折减系数F=2.6时边坡围岩材料属性 d29HEu
MP,EX,10,3E10 ！定义弹性模量 {DR+sE
MP,PRXY,10,0.25 !定义泊松比 * G 4;
MP,DENS,10,2500 ！定义密度 <R $|J|
TB,DP,10 !定义数据表10 H'.d'OE:I
TBDATA,1,3.46E5,19.6 !给数据表10输入数据 fy eS )
!取强度折减系数F=2.8时边坡围岩材料属性 ]Ea6Z
MP,EX,11,3E10 ！定义弹性模量 6=k^gH[g
MP,PRXY,11,0.25 !定义泊松比 _!ed.h.r:
MP,DENS,11,2500 xUU p ?]9y
TB,DP,11 !定义数据表11 y{<7OTA)
TBDATA,1,3.2E5,18.2 !给数据表11输入数据 W *2SlS7
!取强度折减系数F=3.0时边坡围岩材料属性 {lA@I*_lj
MP,EX,12,3E10 ！定义弹性模量 fi)ypv*
MP,PRXY,12,0.25 $Z4p$o
dk
MP,DENS,12,2500 ！定义密度 kW9STN
TB,DP,12 !定义数据表12 g`1i[ Iu2
TBDATA,1,3.0E5,17.1 !给数据表12输入数据 p%_TbH3j`
!建立几何模型 iGIaZ!j aW
!创建边坡线模型线 z& 8#1'
K,1,,,, ！ 创建关键点1 }Z|a?J@CZm
K,2,-800,,, P'+*d#*S
K,3,-800,-800,, /(w5S',EL
K,4,-800,-1200,, Fei$94 a
K,5,1200,-1200,, ,>Q,0bVhH0
K,6,1200,-800,, 4{_5z7ody
K,7,1200,,, ]z,?{S
K,8,1200,400,, ScD9Ct*):C
K,9,400,400,, G[mYx[BTz
l,1,2 !连接关键点1、2生成直线L1 4^r}&9C ~
l,2,3 M E.LS2'n
l,3,4 /[p4. F L
l,4,5 E4Sp^,
l,5,6 e) M1$
l,6,7 7%tn+
l,7,8 6ZE] 7~X
l,8,9 N78Ev7PN
l,9,1 T r8+E;;
l,1,7 >_o_&;=`v
l,3,6 Ni`qU(I'|
!创建边坡面模型 k`d
al,3,4,5,11 !由线L3、L4、L5、L11生成一个面A1 \n5,!,A
al,1,2,11,6,10 !由线L1、L2、L1、L6、L10生成一个面A2 \Z$*8z=
al,7,8,9,10 !由线L7、L8、L9、L10生成一个面A3 F {[Q
/PNUM, area, on !打开面号开关 8[k-8h|
!划分网格生成有限元模型 9M6&+1XE
! 划分边坡围岩2网格 vc<8ApK3V
!设置网格份数 -(jcsqDk
lsel,s,,,3,5,2 !选择线L3和L5 ; O CYx[|
lesize,all,,,5 !把所选择线分为5段 lyy i?/W%
lsel,s,,,4,11,7 !选择线L4和L11 GM5::M]fS
lesize,all,,,26 !把所选择线分为26段 ??F{Gli"C`
mat,2 !给边坡围岩2赋予2号材料特性 TS;MGi 0`}
type,1 !采用单元类型1 <MYD`,$yu
mshkey,0 !设定自由网格划分 q&vr;f B2
mshape,0 !设定四边形网格划分 b(yY.L=K
amesh,1 ！划分面积A1 Bn.8wMB
! 划分边坡围岩1网格 *(rq AB0~
!设置网格份数 + $%o#~
lsel,s,,,2,6,4 !选择线L2和L6 8yd OS
lsel,a,,,1 !附加选择线L1 /N $T[
lesize,all,,,10 !把所选择线分为10段 m4EkL
lsel,s,,,8,10,2 !选择线L8和L10 lr1i DwZV
lesize,all,,,16 !把所选择线分为16段 7}Gy%SJ`
lsel,s,,,7,9,2 bV"0}|A~K
lesize,all,,,12 S&a 44i
mat,1 !给边坡围岩1赋予1号材料特性 ;RR)C@n1
type,1 `
p) #!
mshkey,0 !设定自由网格划分 iZ{D_uxq
mshape,0 !设定四边形网格划分 nPKj%g3h
amesh,2 !划分面积A2 /+rHy7(\
amesh,3 ！划分面积A3 UZyo:*yB
allsel c 9Cp!.#*E
!保存网格模型 (qaY,>je]D
save,Slope-grid.db fE(rDQI
(3)施加约束和荷载 'b- }K DP
!给边坡模型施加约束 ]8RcZn
!边坡两侧施加X方向约束 3V-pLs|
nsel,s,loc,x,-800 !选择X=-800线上所有节点 $I_aHhKt
nsel,a,loc,x,1200 !选择X=1200线上所有节点 [P^ .=F
d,all,ux !对所选择节点约束X方向位移 P%1s6fjU
allsel Eq|_> f@@8
!边坡底部施加约束 'p<(6*,"
nsel,s,loc,y,-1200 !选择Y=--1200线上所有节点 pvX\k X3}
d,all,uy !对所选择节点约束X、Y方向位移 z=!xN5
d,all,ux [j5L}e!T
!施加重力加速度 C&Rv) j
acel,,9.8 hE>Mo$Q(
(4)求解 ~($h9* \
/solu km *$;Nli
!求解设置 $ED<:[3N
antype,static !设定为静力求解 5[0n'uH
nsubst,100 !设定最大子步数为100
hp)3@&T
pred,on !打开时间步长预测器 ;lP)
nropt,full !设定牛顿-拉普森选项 HH>]"mv
nlgeom,on !打开大位移效果 *b> ~L
lnsrch,on !打开线性搜索 =!Ok 079{[
outres,all,all !输出所有项 ( `T:b1
cnvtol,f,,0.005,2,0.5 !力收敛准则设定 n84*[d}t
cnvtol,u,,0.05,2,1 !位移收敛准则设定 _KK G^
u<
!边坡在强度折减系数F=1时求解 :}Z+K*%o-
allsel f@Db._ E
solve !进行求解 'E6)6N
save,F1,db !把F=1时求解结果保存 =l
{>-`:
!边坡在强度折减系数F=1.2时求解 LdA&F&
pI
finish vrH/Z.WD
/solu ,<%],-Lt[
allsel CYz ]tv}g:
asel,s,area,,2,3,1 !选择面积A2、A3 ZMbv1*Vt
MPCHG,3,all ! 把所选择单元材料号改为3 p 5P<3(
allsel Pj^6.f+
solve !进行求解 B&"fPi
save,F1.2,db !把F=1.2时求解结果保存 cBU>/
zI p
!边坡在强度折减系数F=1.4时求解 [N{Rd[{QTL
finish |Kb-oM&^#
/solu I"r*p?
allsel uA,K}sNRZ
asel,s,area,,2,3,1 !选择面积A2、A3 lE /"
MPCHG,4,all !把所选择单元材料号改为4 |06G )r&
allsel CB |Z~_Bm
solve !进行求解 H$Q_K<V
save,F1.4,db !把F=1.4时求解结果保存 >*EcX 3
!边坡在强度折减系数F=1.6时求解 E: #VS~
finish f.Jz]WXw,
/solu w
J; y 4
allsel TTS }, `
asel,s,area,,2,3,1 !选择面积A2、A3 i,,U D
MPCHG,5,all !把所选择单元材料号改为5 hV_0f_Og
allsel xCGvLvFn
solve 6:@tHUm
save,F1.6,db !把F=1.6时求解结果保存 uYF _sf
!边坡在强度折减系数F=1.8时求解 {~ VgXkjsC
finish V*LpO 8=
/solu aabnlOVw
allsel &vo]l~.
asel,s,area,,2,3,1 !选择面积A2、A3 F8.Fp[_tM
MPCHG,6,all !把所选择单元材料号改为6 #TRPq>Xz D
allsel s <td n[d
solve jf@#&%AC9
save,F1.8,db !把F=1.8时求解结果保存 *SL v$A
!边坡在强度折减系数F=2.0时求解 ` O-$qT, _
finish [/6IEt3}B
/solu njX!Ez
allsel ^~?VD
asel,s,area,,2,3,1 !选择面积A2、A3 y~jTI[kS
MPCHG,7,all !把所选择单元材料号改为7 B usxg?=
allsel Vf,t=$.[Q
solve T-uI CMEf
save,F2.0,db !把F=2.0时求解结果保存 NJwcb=*
!边坡在强度折减系数F=2.2时求解 MX]<tR `
finish u<:R Sg
/solu [{L4~(uU8
allsel ?=-18@:.ss
asel,s,area,,2,3,1 !选择面积A2、A3 Ub%5# <k|-
MPCHG,8,all !把所选择单元材料号改为8 v~9PS2
allsel > }Za)
solve 4 cl}ouG
save,F2.2,db !把F=2.2时求解结果保存 vQ$ FMKz7
!边坡在强度折减系数F=2.4时求解 `#E1FB2M
finish AKej Wh
/solu vaW, O/F
allsel 8W9kd"=U
asel,s,area,,2,3,1 !选择面积A2、A3 %TP0i#J
MPCHG,9,all !把所选择单元材料号改为9 _}8O15B|
allsel D`V03}\-
solve \GEFhM4)
save,F2.4,db !把F=2.4时求解结果保存 ULNU '6
!边坡在强度折减系数F=2.6时求解 >iDV8y
finish
T32C=7
/solu =LC:1zn4
allsel q",n:=PL
asel,s,area,,2,3,1 !选择面积A2、A3 De6WC*trq
MPCHG,10,all !把所选择单元材料号改为10 B~G ?&"]
allsel /q0[T{Wz$
solve n>" 0 y^v
save,F2.6,db !把F=2.6时求解结果保存 c= 0S]_
!边坡在强度折减系数F=2.8时求解 !"Jne' f
finish [2"a~o\
/solu An{`'U(l
allsel qk<( iVUO
asel,s,area,,2,3,1 !选择面积A2、A3 @2na r<
MPCHG,11,all !把所选择单元材料号改为11 N`h, 2!(j
allsel tt $DWmm
solve 9@9(zUS|
save,F2.8,db !把F=2.8时求解结果保存 [Ab q("9p\
!边坡在强度折减系数F=3.0时求解 w^6rgCl
finish &!SdO<agZ
/solu b3N1SC:Wn
allsel kj{z;5-dl
asel,s,area,,2,3,1 !选择面积A2、A3 J^m#984
MPCHG,12,all !把所选择单元材料号改为12 i 3?=up!
allsel ?); 6]"k:3
solve i7 _Nv
save,F3.0,db !把F=3.0时求解结果保存 }6 *+>?
(5）后处理 |3<tDq@+
/post1 ！进入后处理 gdPv,p19L
!边坡在强度折减系数F=1时结果分析 hd^?mZ
Resume,'F1','db' !读入边坡在强度折减系数F=1时 .7 6T<j_
set,1,last !读入后一个子步 roSdcQTeT
pldisp,1 !绘制边坡模型变形图 z~\ a]MB
plnsol,u,x !绘制边坡模型水平方向位移云图 ?%-VSL>$w=
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 XHKiz2Pc1
!边坡在强度折减系数F=1.2时结果分析 SVB> 1s9F
Resume,'F1.2','db' !读入边坡在强度折减系数F=1.2时 SaceIV %(
set,1,last !读入后一个子步 1zqIB")s>
pldisp,1 !绘制边坡模型变形图 ce7$r*@!
plnsol,u,x !绘制边坡模型水平方向位移云图 +i HZ*
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 Vb yGr~t
!边坡在强度折减系数F=1.4时结果分析 a~ sU
Resume,'F1.4','db' !读入边坡在强度折减系数F=1.4时 W5z<+8R
set,1,last !读入后一个子步 GrPKJ~{6
pldisp,1 !绘制边坡模型变形图 k<(G)7'gm
plnsol,u,x !绘制边坡模型水平方向位移云图 #; ~`+[y?\
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 "*UN\VV+s
!边坡在强度折减系数F=1.6时结果分析 }-4@EC>
Resume,'F1.6','db' !读入边坡在强度折减系数F=1.6时 j Gy%O3/
set,1,last !读入后一个子步 Hz~?"ts@;
pldisp,1 !绘制边坡模型变形图 .
7*k}@k
plnsol,u,x !绘制边坡模型水平方向位移云图 %Z_/MNI
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 ?&A)%6` ~
!边坡在强度折减系数F=1.8时结果分析 Lu?MRF
f
Resume,'F1.8','db' !读入边坡在强度折减系数F=1.8时 {,xI|u2R
set,1,last !读入后一个子步 utOATjB.z
pldisp,1 !绘制边坡模型变形图 *9T a0e*
plnsol,u,x !绘制边坡模型水平方向位移云图 w{TZN {Y
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 M~/7thP{
!边坡在强度折减系数F=2.0时结果分析 b* (~8JxZ
Resume,'F2.0','db' !读入边坡在强度折减系数F=2.0时 m03D+@F
set,1,last !读入后一个子步 w(Jf;[o
pldisp,1 !绘制边坡模型变形图 XK#~w:/fB
plnsol,u,x !绘制边坡模型水平方向位移云图 jEU`ko_
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 @ByD=
!边坡在强度折减系数F=2.2时结果分析 W;L<zFFbU)
Resume,'F2.2','db' !读入边坡在强度折减系数F=2.2时 \"k[y+O],4
set,1,last !读入后一个子步 tV.qdy/]}
pldisp,1 !绘制边坡模型变形图 d ? Uj3G
plnsol,u,x !绘制边坡模型水平方向位移云图 1'Sr0
oEd3
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 ZW+M<G
!边坡在强度折减系数F=2.4时结果分析 GO&~)Vh&7
Resume,'F2.4','db' !读入边坡在强度折减系数F=2.4时 :DWvH,{+&
set,1,last !读入后一个子步 }OZut! _
pldisp,1 !绘制边坡模型变形图 3T1t !q4/5
plnsol,u,x !绘制边坡模型水平方向位移云图 R%Xhd cn7
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 f
wE
b
!边坡在强度折减系数F=2.6时结果分析 r7|_Fm Qf
Resume,'F2.6','db' !读入边坡在强度折减系数F=2.6时 }Br=eaY
set,1,last !读入后一个子步 sk aPC#u
pldisp,1 !绘制边坡模型变形图 ({& \~"
plnsol,u,x !绘制边坡模型水平方向位移云图 r{Lr Q
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 py|ORVN(Z
!边坡在强度折减系数F=2.8时结果分析 =bKz$
_W
Resume,'F2.8','db' !读入边坡在强度折减系数F=2.8时 ;@ <E
set,1,last !读入后一个子步 7 +]+S`p
pldisp,1 !绘制边坡模型变形图 2bv/ -^
plnsol,u,x !绘制边坡模型水平方向位移云图 =: fN
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 mH`K~8pRg
!边坡在强度折减系数F=3.0时结果分析 &,*G}6wa;&
Resume,'F3.0','db' !读入边坡在强度折减系数F=3.0时 dMd2a4
set,1,last !读入后一个子步 &JUHm_wd&S
pldisp,1 !绘制边坡模型变形图 !m {d6C[
plnsol,u,x !绘制边坡模型水平方向位移云图 LOO<)XFJ
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 ~?6M4!u
K%jh 6c8
O Ko)p`BX
y V.E+~y
CU`yi.)T{
WNnB
s
u B~C8}
