各位朋友： ! 345
大家好！我通过ansys用有限元强度折减法分析了一个边坡的稳定性，结果不同的工况下得出的塑性区是一样的，我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激！附命令流和塑性区如下：(1)创建物理环境 sXR}#*8p
/COM, Structural ! 指定结构分析 `E5vO1Pl
/TITLE, Slope stability Analysis ! 定义工作标题 G~19Vv*;
/FILNAM,support,1 ! 定义工作文件名 KZI-/H+
(2)建立模型 )B5(V5-!|
!进入前处理器 k^Uk= )9
/PREP7 e%v0EJ},
!定义单元类型 c\N-B,m&
ET,1,PLANE82 !定义边坡围岩单元 FS6I?q#tQ
KEYOPT,1,3,2 !求解类型为平面应变 fR,7l9<%Zp
!定义材料属性 #W[C;f|,
!边坡围岩1材料属性 V6tUijz
MP,EX,1,3E10 2D"\Ox
MP,PRXY,1,0.25 i&zJwUr(<
MP,DENS,1,2500 !定义塑性模型 -"w&g0Z
TB,DP,1 !定义数据表1 ufXU
TBDATA,1,9E5,42.8 !给数据表输入数据 7 w5 L?,a
!边坡围岩2材料属性 !采用弹性模型 ^Z G 3{>
MP,EX,2,3.2E10 ！定义弹性模量 \:_!!
MP,PRXY,2,0.24 !定义泊松比 ziG]BZ
MP,DENS,2,2700 !定义密度 5dEek7wnf
!取强度折减系数F=1.2时边坡围岩材料属性 ~MZ.988:<
MP,EX,3,3E10 fXB64MNo
MP,PRXY,3,0.25 rtk1 8U-
MP,DENS,3,2500 =d1i<iw?-
TB,DP,3 !定义数据表3 m^Rf6O^
TBDATA,1,7.5E5,37.7 !给数据表输入数据 @^K_>s9B
!取强度折减系数F=1.4时边坡围岩材料属性 k4BiH5\hA
MP,EX,4,3E10 [p 8 fg!|
MP,PRXY,4,0.25 !定义泊松比 $6yr:2Xvt
MP,DENS,4,2500 W*Ce1
TB,DP,4 !定义数据表4 ~w}Zv0
TBDATA,1,6.4E5,33.5 !给数据表输入数据 <UBB&}R0
!取强度折减系数F=1.6时边坡围岩材料属性 ZO !)G
MP,EX,5,3E10 AGgL`sP
MP,PRXY,5,0.25 zXT[}J VV
MP,DENS,5,2500 ！定义密度 zK ir
TB,DP,5 !定义数据表5 _|KeB(W
TBDATA,1,5.6E5,30 !给数据表输入数据 Jk=d5B
!取强度折减系数F=1.8时边坡围岩材料属性 )! C |DSw
MP,EX,6,3E10 ！定义弹性模量 nISfRX U;
MP,PRXY,6,0.25 Jkub|w#QH
MP,DENS,6,2500 H^0`YQJ3
TB,DP,6 !定义数据表6 ?KXgG' !!
TBDATA,1,5.0E5,27.2 FW!1 0K?
!取强度折减系数F=2.0时边坡围岩材料属性 & <Jvaf_=
MP,EX,7,3E10 ]?j[P=\
MP,PRXY,7,0.25 !定义泊松比 "jAEZ
MP,DENS,7,2500 =y1/V'2E
TB,DP,7 !定义数据表7 Xd@x(T~'X
TBDATA,1,4.5E5,24.8 GoRSLbCUR
!取强度折减系数F=2.2时边坡围岩材料属性 a3(q;^v
MP,EX,8,3E10 F(+dX4$
MP,PRXY,8,0.25 @P
xX]e
MP,DENS,8,2500 mc}r15:<
TB,DP,8 !定义数据表8 Czt>?8x`
TBDATA,1,4.09E5,22.8 YLe$Vv735
!取强度折减系数F=2.4时边坡围岩材料属性 MDh^ic5
MP,EX,9,3E10 Mf.:y
MP,PRXY,9,0.25 #wL8=QTcNC
MP,DENS,9,2500 =]hPX
TB,DP,9 !定义数据表9 *Q:EICDE7
TBDATA,1,3.6E5,21.1 !给数据表输入数据 =U<6TP]{
!取强度折减系数F=2.6时边坡围岩材料属性 jthGNVZ
MP,EX,10,3E10 ！定义弹性模量 t?cO>4*|
MP,PRXY,10,0.25 !定义泊松比 zoO9N oUHW
MP,DENS,10,2500 ！定义密度 A]mXV4RmI
TB,DP,10 !定义数据表10
O^I%Xk
TBDATA,1,3.46E5,19.6 !给数据表10输入数据 7s'r3}B`
!取强度折减系数F=2.8时边坡围岩材料属性 2ZZF hj
MP,EX,11,3E10 ！定义弹性模量 uY*|bD`6&
MP,PRXY,11,0.25 !定义泊松比 t 4tXLI;'
MP,DENS,11,2500 cT,5xp"a
TB,DP,11 !定义数据表11 2 NrM se
TBDATA,1,3.2E5,18.2 !给数据表11输入数据 Odj4)
!取强度折减系数F=3.0时边坡围岩材料属性 o0 Pc^
MP,EX,12,3E10 ！定义弹性模量 bhc
.UmH
MP,PRXY,12,0.25 +}@6V4BRn
MP,DENS,12,2500 ！定义密度 ]2'{W]m
TB,DP,12 !定义数据表12 Jz' 8|o;^
TBDATA,1,3.0E5,17.1 !给数据表12输入数据 r d4\N2- 6
!建立几何模型 J3#
!创建边坡线模型线 @Z %I g
K,1,,,, ！ 创建关键点1 , K[}Bz
K,2,-800,,, U{~R39
K,3,-800,-800,, 6$"0!fl>
K,4,-800,-1200,, _+x&[^gjP
K,5,1200,-1200,, \`,, r_tO
K,6,1200,-800,, o9D] \PdL>
K,7,1200,,, 'UL"yM
K,8,1200,400,, 'CC;=@J
K,9,400,400,, O(Vi/r2:e
l,1,2 !连接关键点1、2生成直线L1 E/N*n!sV
l,2,3 Q]2v]PJ6"
l,3,4 z\Y-8a.]
l,4,5 bx8 |_K*^
l,5,6 F!qt#Sw!\
l,6,7 !m tX*;b(e
l,7,8 > aV
Q
l,8,9 *Wmn!{\ g
l,9,1 ^q
?xi5 w
l,1,7 YF(TG]?6
l,3,6 hu''"/raM
!创建边坡面模型 UXN!iU)
al,3,4,5,11 !由线L3、L4、L5、L11生成一个面A1 7K}Sk
al,1,2,11,6,10 !由线L1、L2、L1、L6、L10生成一个面A2 7s-ZRb[)1
al,7,8,9,10 !由线L7、L8、L9、L10生成一个面A3 )a'c_ 2 [
/PNUM, area, on !打开面号开关 Bi|XdS$G
!划分网格生成有限元模型 z4[S 02s
! 划分边坡围岩2网格 $l!+SLK
!设置网格份数 % $.]g
lsel,s,,,3,5,2 !选择线L3和L5 D_4UM#Tw
lesize,all,,,5 !把所选择线分为5段 {Tym#
lsel,s,,,4,11,7 !选择线L4和L11 dr8`;$;G*
lesize,all,,,26 !把所选择线分为26段 7SpF&
mat,2 !给边坡围岩2赋予2号材料特性 I Lq"/S.
type,1 !采用单元类型1 pCm |t!,
mshkey,0 !设定自由网格划分 +x" cWOg
mshape,0 !设定四边形网格划分 ]>\!} \R<
amesh,1 ！划分面积A1 JYU0&nZl4
! 划分边坡围岩1网格 tr $~INe
!设置网格份数 =/]d\JSp
lsel,s,,,2,6,4 !选择线L2和L6 #.z`clK#
lsel,a,,,1 !附加选择线L1 ,6FmU$
Kn
lesize,all,,,10 !把所选择线分为10段 YQk<1./}I
lsel,s,,,8,10,2 !选择线L8和L10 ;~5w`F)
lesize,all,,,16 !把所选择线分为16段 +:fr(s!OE
lsel,s,,,7,9,2 }^Kye23
lesize,all,,,12 rezH5d6z62
mat,1 !给边坡围岩1赋予1号材料特性 STH?X]
/
type,1 = ;"$t_t
mshkey,0 !设定自由网格划分 Q g;?C
mshape,0 !设定四边形网格划分 #{u>
amesh,2 !划分面积A2 sVJwe\!
amesh,3 ！划分面积A3 tZn=[X~Vw@
allsel e.:S BXZ
!保存网格模型 y vz2eAXa
save,Slope-grid.db <xWBS/K
(3)施加约束和荷载 F D*w4U5
!给边坡模型施加约束 W2\ Q-4D
!边坡两侧施加X方向约束 ,
,=7deR
nsel,s,loc,x,-800 !选择X=-800线上所有节点 TWFi.w4pY
nsel,a,loc,x,1200 !选择X=1200线上所有节点 -H;p +XAY
d,all,ux !对所选择节点约束X方向位移 ^@0-E@ {c
allsel ]$ gBX=
!边坡底部施加约束 ;wkoQ8FD9
nsel,s,loc,y,-1200 !选择Y=--1200线上所有节点 L?gak@E
d,all,uy !对所选择节点约束X、Y方向位移 r]+N(&q
d,all,ux *K1GX
!施加重力加速度 4,pS C
acel,,9.8 zIjUfgO/M
(4)求解 (0/)vZc
/solu Q\4nduQ
!求解设置 A(>kp=~
antype,static !设定为静力求解 y
2v69nu~q
nsubst,100 !设定最大子步数为100 ]jL`*tI \S
pred,on !打开时间步长预测器 ~Q)137u]P
nropt,full !设定牛顿-拉普森选项 rM?ox
nlgeom,on !打开大位移效果 8!uqR!M<C
lnsrch,on !打开线性搜索 V=g<3R&
outres,all,all !输出所有项 q[w.[]
cnvtol,f,,0.005,2,0.5 !力收敛准则设定 j,c8_;X!
cnvtol,u,,0.05,2,1 !位移收敛准则设定 ntT~_Ba8;u
!边坡在强度折减系数F=1时求解 {_i.IPp~
allsel gAWrn^2L5
solve !进行求解 |p7k2wzN
save,F1,db !把F=1时求解结果保存 h#iFp9N
!边坡在强度折减系数F=1.2时求解 y8.(fil NB
finish ZT;:Hxv0N
/solu ,awp)@VG7
allsel < BNCo5*
asel,s,area,,2,3,1 !选择面积A2、A3 CH/*MA
MPCHG,3,all ! 把所选择单元材料号改为3 P6cc8x9g(
allsel <M4Qc12jP
solve !进行求解 7>Oa, \
save,F1.2,db !把F=1.2时求解结果保存 /K./k!'z
!边坡在强度折减系数F=1.4时求解 |:?JSi0
finish ,wvzY7%
/solu "|:I]ZB
allsel 8+7=yN(
asel,s,area,,2,3,1 !选择面积A2、A3 z|gG%fM
MPCHG,4,all !把所选择单元材料号改为4 fm%1 vM$[J
allsel NhDM h8=$^
solve !进行求解 Cyw
cJ
save,F1.4,db !把F=1.4时求解结果保存 :jp4 !0w
!边坡在强度折减系数F=1.6时求解 ^,.G<2Kx&
finish M;i4ss,}!
/solu d=B
DR^/wA
allsel z
a^s%^:yK
asel,s,area,,2,3,1 !选择面积A2、A3 iqj
ZC80
MPCHG,5,all !把所选择单元材料号改为5 8 O5@FU
3
allsel I3ZbHb-)_,
solve _4VS.~}/R
save,F1.6,db !把F=1.6时求解结果保存 uBe1{Z
!边坡在强度折减系数F=1.8时求解 ) =)=]|3
finish xe3t_y
/solu #n_uELE
allsel "T_OLegdK
asel,s,area,,2,3,1 !选择面积A2、A3 Rr'^l ]
MPCHG,6,all !把所选择单元材料号改为6 "/-T{p;.
allsel /:j9 #kj
solve bl)iji`]
save,F1.8,db !把F=1.8时求解结果保存 a ges-Z_X
!边坡在强度折减系数F=2.0时求解 FGP~^Dr/
finish ped3}i+|]
/solu 4l~0LdYXKm
allsel K&WNtk3hT
asel,s,area,,2,3,1 !选择面积A2、A3 xgeKz^,
MPCHG,7,all !把所选择单元材料号改为7 LFx*_3a
allsel 75pz' Cb
solve gZs UX^%
save,F2.0,db !把F=2.0时求解结果保存 m8|&z{
!边坡在强度折减系数F=2.2时求解 (y xrK
finish H' [#x2
/solu .RNr^ *AQ
allsel +|w-1&-
asel,s,area,,2,3,1 !选择面积A2、A3 *&vySyt
MPCHG,8,all !把所选择单元材料号改为8 6jIW)C
allsel ul',!js?
solve = yH#Iil
save,F2.2,db !把F=2.2时求解结果保存 Gv};mkX[N
!边坡在强度折减系数F=2.4时求解 #!%\97ZR
finish a Dik1Q
/solu }m~2[5q%/
allsel !y>MchNv
asel,s,area,,2,3,1 !选择面积A2、A3 p<@ 0b
MPCHG,9,all !把所选择单元材料号改为9 \5wC&|WEB
allsel |# zznT"
solve :%?\Wj5HW
save,F2.4,db !把F=2.4时求解结果保存 +I?T |Iin
!边坡在强度折减系数F=2.6时求解 ktr l |
finish u$Za hN!
/solu H lw0i a
allsel Z ~-A*{u?
asel,s,area,,2,3,1 !选择面积A2、A3 v<`1z?dch
MPCHG,10,all !把所选择单元材料号改为10 &@dW d
allsel ]eJjffx
solve &x(^=sTHI
save,F2.6,db !把F=2.6时求解结果保存 !:[kS1s>M
!边坡在强度折减系数F=2.8时求解 esM< .
finish til L7
/solu (!nhU
allsel 79>8tOuo
asel,s,area,,2,3,1 !选择面积A2、A3 XVfp* `
MPCHG,11,all !把所选择单元材料号改为11 / !*gH1 s
allsel ?V}AwLX}
solve p?X`f#
save,F2.8,db !把F=2.8时求解结果保存 I o z
rZ
!边坡在强度折减系数F=3.0时求解 G([!(8&2Y
finish Mp V6Vbp
/solu 2_ x~y|<9
allsel myj^c>1Iz
asel,s,area,,2,3,1 !选择面积A2、A3 xCd9b:jG
MPCHG,12,all !把所选择单元材料号改为12 U 6y
;V
allsel 0-^wY8n-=
solve U-$ B"w &
save,F3.0,db !把F=3.0时求解结果保存 90pk
(5）后处理 l|[8'*]r!
/post1 ！进入后处理 hupYiI~
!边坡在强度折减系数F=1时结果分析 2 HNH@K
Resume,'F1','db' !读入边坡在强度折减系数F=1时 GMZj@q
set,1,last !读入后一个子步 # Z*nc0C
pldisp,1 !绘制边坡模型变形图 cN> z`x l
plnsol,u,x !绘制边坡模型水平方向位移云图 a?IL6$z
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 "/)}Cc,L
!边坡在强度折减系数F=1.2时结果分析 Z81{v<c;
Resume,'F1.2','db' !读入边坡在强度折减系数F=1.2时 'S
f
set,1,last !读入后一个子步 ]byj[Gd
pldisp,1 !绘制边坡模型变形图 ZR3x;$I~4
plnsol,u,x !绘制边坡模型水平方向位移云图 q >9F2 1 W
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 C <.t'|
!边坡在强度折减系数F=1.4时结果分析 [p +h b
Resume,'F1.4','db' !读入边坡在强度折减系数F=1.4时 7b_Ihv
set,1,last !读入后一个子步 (\si/&
pldisp,1 !绘制边坡模型变形图 qR~s& SC#
plnsol,u,x !绘制边坡模型水平方向位移云图 fU+A~oL%I
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 6c3+q+#J2
!边坡在强度折减系数F=1.6时结果分析 .g7ebh6D
Resume,'F1.6','db' !读入边坡在强度折减系数F=1.6时 ZcX qH7`r
set,1,last !读入后一个子步 7]q $ sQ
pldisp,1 !绘制边坡模型变形图 U~SOHfZ%(
plnsol,u,x !绘制边坡模型水平方向位移云图 hwmpiyu
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 wNuS'P_(:T
!边坡在强度折减系数F=1.8时结果分析 4 g#pQ
Resume,'F1.8','db' !读入边坡在强度折减系数F=1.8时 p1=sDsLL
set,1,last !读入后一个子步 o y-Qy
pldisp,1 !绘制边坡模型变形图 Ah2%LXdHA
plnsol,u,x !绘制边坡模型水平方向位移云图 a/.O, &3
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 U+!H/R)(
!边坡在强度折减系数F=2.0时结果分析 eTc0u;{V
Resume,'F2.0','db' !读入边坡在强度折减系数F=2.0时 R,hX *yVq
set,1,last !读入后一个子步 k MS[
pldisp,1 !绘制边坡模型变形图 NC 0H5
plnsol,u,x !绘制边坡模型水平方向位移云图 "-N)TIzLX
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 x0])&':!
!边坡在强度折减系数F=2.2时结果分析 9's/~ T
Resume,'F2.2','db' !读入边坡在强度折减系数F=2.2时 8u::f`vi
set,1,last !读入后一个子步 %Nb hR(
pldisp,1 !绘制边坡模型变形图 MR90 }wXE
plnsol,u,x !绘制边坡模型水平方向位移云图 0;-S){
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 RN?z)9 !
!边坡在强度折减系数F=2.4时结果分析 z/7H/~d
Resume,'F2.4','db' !读入边坡在强度折减系数F=2.4时 |4i,Vkfhe
set,1,last !读入后一个子步
")U`W gx
pldisp,1 !绘制边坡模型变形图 $ V"~\h8
plnsol,u,x !绘制边坡模型水平方向位移云图 &8Cuu$T9)
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 sa~.qmqu
!边坡在强度折减系数F=2.6时结果分析 =s P6
Resume,'F2.6','db' !读入边坡在强度折减系数F=2.6时 &/7 AW(?
set,1,last !读入后一个子步 A#mf*]'
pldisp,1 !绘制边坡模型变形图 wR;_x x
plnsol,u,x !绘制边坡模型水平方向位移云图 Kl[WscR
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 x@I@7Pvo3
!边坡在强度折减系数F=2.8时结果分析 /hue]ZaQq
Resume,'F2.8','db' !读入边坡在强度折减系数F=2.8时 Fc;)p88[
set,1,last !读入后一个子步 )> )_>[
pldisp,1 !绘制边坡模型变形图 `A\
!Gn?
plnsol,u,x !绘制边坡模型水平方向位移云图 K%<Z"2!+
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 lA39$oJ
!边坡在强度折减系数F=3.0时结果分析 )T6+}
Resume,'F3.0','db' !读入边坡在强度折减系数F=3.0时 saH +C@_,
set,1,last !读入后一个子步 c FjC
pldisp,1 !绘制边坡模型变形图 rs@,<DV)u
plnsol,u,x !绘制边坡模型水平方向位移云图 /tno`su;
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 wovWEtVBU
4QnJ;&~
- v9V/LJ
n8zh;vuJ
uBk$zs
OC'cP[$ _
jZ < *XX
