各位朋友： j +Ro?
大家好！我通过ansys用有限元强度折减法分析了一个边坡的稳定性，结果不同的工况下得出的塑性区是一样的，我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激！附命令流和塑性区如下：(1)创建物理环境 Dos';9Uq
/COM, Structural ! 指定结构分析 Bx0=D:j
/TITLE, Slope stability Analysis ! 定义工作标题 2z9N/ SyN
/FILNAM,support,1 ! 定义工作文件名 mx Tuwx
(2)建立模型 aTE;Gy,W
!进入前处理器 R5fZ }C7
/PREP7 ;et(Yi;9
!定义单元类型 |/!RN[<
ET,1,PLANE82 !定义边坡围岩单元 I6[=tB
KEYOPT,1,3,2 !求解类型为平面应变 6jz~q~ I
!定义材料属性 @0iXqM#jH
!边坡围岩1材料属性 5>rjL ;
MP,EX,1,3E10 T@Th?
MP,PRXY,1,0.25 je=XZ's,i~
MP,DENS,1,2500 !定义塑性模型 iG=XRctgj)
TB,DP,1 !定义数据表1 +}NQ |y V
TBDATA,1,9E5,42.8 !给数据表输入数据 A "~Oi
!边坡围岩2材料属性 !采用弹性模型 i/F ]. Sag
MP,EX,2,3.2E10 ！定义弹性模量 H2E'i\
MP,PRXY,2,0.24 !定义泊松比 .$ 5*v
MP,DENS,2,2700 !定义密度 6 ?= ^8
!取强度折减系数F=1.2时边坡围岩材料属性 }*56 DX
MP,EX,3,3E10 %(g!,!l)
MP,PRXY,3,0.25
_c QhT
MP,DENS,3,2500 tj[-|h
TB,DP,3 !定义数据表3 ;ZR^9%+y9
TBDATA,1,7.5E5,37.7 !给数据表输入数据 W";Po)YC
!取强度折减系数F=1.4时边坡围岩材料属性 i&JI"Dd7
MP,EX,4,3E10 CV *
MP,PRXY,4,0.25 !定义泊松比 (0zYS_m A
MP,DENS,4,2500 (z X&feq
TB,DP,4 !定义数据表4 > `+lEob
TBDATA,1,6.4E5,33.5 !给数据表输入数据 2G }@s. iE
!取强度折减系数F=1.6时边坡围岩材料属性 =}[m_rp&
MP,EX,5,3E10 }62Q{>`
MP,PRXY,5,0.25 +<f !#4T
MP,DENS,5,2500 ！定义密度 u 1/q8'RW
TB,DP,5 !定义数据表5 B {lL}"++0
TBDATA,1,5.6E5,30 !给数据表输入数据 9, 0}}3J
!取强度折减系数F=1.8时边坡围岩材料属性 ]VK9d;0D
MP,EX,6,3E10 ！定义弹性模量 "oT]_WHqo
MP,PRXY,6,0.25 LH>h]OTQF
MP,DENS,6,2500 qnf\ K}
TB,DP,6 !定义数据表6 FO?I}G22
TBDATA,1,5.0E5,27.2 ~D! Y]
SK
!取强度折减系数F=2.0时边坡围岩材料属性 D0 rqte
MP,EX,7,3E10 l6HT}x7OiH
MP,PRXY,7,0.25 !定义泊松比 )iCg,?SSw=
MP,DENS,7,2500 Owm 2/
TB,DP,7 !定义数据表7 >I9|N}I
TBDATA,1,4.5E5,24.8 3E:+DF-Z\
!取强度折减系数F=2.2时边坡围岩材料属性 7kQ,D,c'
MP,EX,8,3E10 H^8 t/h
MP,PRXY,8,0.25 v QDkZ
MP,DENS,8,2500 gYfO a`k
TB,DP,8 !定义数据表8 8:}$L)[V
TBDATA,1,4.09E5,22.8 6>- Gi
!取强度折减系数F=2.4时边坡围岩材料属性 w m/>_
MP,EX,9,3E10 L'}^Av_+
MP,PRXY,9,0.25 j'&a)-Wx_
MP,DENS,9,2500 $ q$\
TB,DP,9 !定义数据表9 E0"DHjR
TBDATA,1,3.6E5,21.1 !给数据表输入数据 a: [m;
!取强度折减系数F=2.6时边坡围岩材料属性 {;?bC'
MP,EX,10,3E10 ！定义弹性模量 W[.UM
MP,PRXY,10,0.25 !定义泊松比 DP4l
%2m0
MP,DENS,10,2500 ！定义密度 j0IuuJ+
TB,DP,10 !定义数据表10
< zF /at
TBDATA,1,3.46E5,19.6 !给数据表10输入数据 ,:"c"
!取强度折减系数F=2.8时边坡围岩材料属性 N;XJMk_ H
MP,EX,11,3E10 ！定义弹性模量 '2c4
4F)i
MP,PRXY,11,0.25 !定义泊松比 u^&,~n@n7
MP,DENS,11,2500 ;#6j9M0
TB,DP,11 !定义数据表11 ~n!7 ?4%U
TBDATA,1,3.2E5,18.2 !给数据表11输入数据 `jH 0FJQ
!取强度折减系数F=3.0时边坡围岩材料属性 ({p @Ay
MP,EX,12,3E10 ！定义弹性模量 9/L J tM
MP,PRXY,12,0.25 %'<m[wf^ o
MP,DENS,12,2500 ！定义密度 dG'5: ,n/
TB,DP,12 !定义数据表12 i@`T_&6l
TBDATA,1,3.0E5,17.1 !给数据表12输入数据 K6z-brvw "
!建立几何模型 *saO~.-;4
!创建边坡线模型线 Nyo ,6 AA
K,1,,,, ！ 创建关键点1 ?kM53zbT#
K,2,-800,,, z}*74lhF
K,3,-800,-800,, vCtag]H2@
K,4,-800,-1200,, f=*xdOB3
K,5,1200,-1200,, XW9
[VUW~
K,6,1200,-800,, O{`r.H1',
K,7,1200,,, OL)M`eVQ'
K,8,1200,400,, 3ZU<u;
K,9,400,400,, 'aW<C>
l,1,2 !连接关键点1、2生成直线L1 /-jk_8@a
l,2,3 <IL$8a
l,3,4 `jB2'
l,4,5 p&ml$N9fd
l,5,6 ~",,&>#[K
l,6,7 [sH[bmLR
l,7,8 Gm1vVHAxv
l,8,9 Ws(#T hA
l,9,1 XDi[ Iyj
l,1,7 &-1;3+#w
l,3,6 u6( >?r-
!创建边坡面模型 ;l 6tZ]-"
al,3,4,5,11 !由线L3、L4、L5、L11生成一个面A1 27Ve $Q8]v
al,1,2,11,6,10 !由线L1、L2、L1、L6、L10生成一个面A2 O <v9i4*
al,7,8,9,10 !由线L7、L8、L9、L10生成一个面A3 ,ZL g=
/PNUM, area, on !打开面号开关 oZ (T` 5
!划分网格生成有限元模型 " | l-NUe
! 划分边坡围岩2网格 ~
-hH#5
!设置网格份数 /N&)r wc
lsel,s,,,3,5,2 !选择线L3和L5 zWh[U'6
lesize,all,,,5 !把所选择线分为5段 #7A_p8
lsel,s,,,4,11,7 !选择线L4和L11 =a .avOZ
lesize,all,,,26 !把所选择线分为26段 va/$dD9
mat,2 !给边坡围岩2赋予2号材料特性 icPg<>TQ
type,1 !采用单元类型1 *r k!`n&
mshkey,0 !设定自由网格划分 :3x |U,wC
mshape,0 !设定四边形网格划分 @gZ<!g/vza
amesh,1 ！划分面积A1 8Dq;QH}
! 划分边坡围岩1网格 g.d%z
!设置网格份数 0-W{( xy@4
lsel,s,,,2,6,4 !选择线L2和L6 75V?K
lsel,a,,,1 !附加选择线L1 OxraaN`
lesize,all,,,10 !把所选择线分为10段 @*gm\sU4
lsel,s,,,8,10,2 !选择线L8和L10 kG7q4jFwP
lesize,all,,,16 !把所选择线分为16段 |[VtYV _{
lsel,s,,,7,9,2 $(6 .K- D
lesize,all,,,12 |E)Es!dr
mat,1 !给边坡围岩1赋予1号材料特性 ui :
type,1 Uaho.(_GP
mshkey,0 !设定自由网格划分 qi\!<clv
mshape,0 !设定四边形网格划分 0} UJP
amesh,2 !划分面积A2 Tej-mr3P
amesh,3 ！划分面积A3 [~ sXjaL8
allsel l ^*GqP5
!保存网格模型 S3=M k~_&
save,Slope-grid.db l*|^mx^Q
(3)施加约束和荷载 o[bG(qHZ
!给边坡模型施加约束 D %`64R
!边坡两侧施加X方向约束 6N&S3<c4JO
nsel,s,loc,x,-800 !选择X=-800线上所有节点 wR?M2*ri
nsel,a,loc,x,1200 !选择X=1200线上所有节点 -T4{PM
d,all,ux !对所选择节点约束X方向位移 bqXCe\#
allsel $Wy7z^ t
!边坡底部施加约束 o?f7_8fG
nsel,s,loc,y,-1200 !选择Y=--1200线上所有节点 a<%WFix
d,all,uy !对所选择节点约束X、Y方向位移 {g2cm'hD
d,all,ux x7X"'1U
!施加重力加速度 5O
Ob(
acel,,9.8 fn zj@_{|
(4)求解 2B8p3 A
/solu HM])m>KeT
!求解设置 pCz;km
antype,static !设定为静力求解 !<! sB)
nsubst,100 !设定最大子步数为100 q wu++9BM
pred,on !打开时间步长预测器 )>atoA
nropt,full !设定牛顿-拉普森选项 }ZEh^zdz8
nlgeom,on !打开大位移效果 *J!oV0# 1
lnsrch,on !打开线性搜索 ,|_ewye
outres,all,all !输出所有项 w| x=^
cnvtol,f,,0.005,2,0.5 !力收敛准则设定 FDaHsiI:
cnvtol,u,,0.05,2,1 !位移收敛准则设定 z= -u89]
!边坡在强度折减系数F=1时求解 XnQd(B`M
allsel (}fbs/8\p
solve !进行求解 p3:x\P<|
save,F1,db !把F=1时求解结果保存 ~ X8U @f
!边坡在强度折减系数F=1.2时求解 PKR0y%Ar
finish cI4%z eR
/solu xDu11 W+g
allsel D>,]EE-
asel,s,area,,2,3,1 !选择面积A2、A3 |2X Et\P
MPCHG,3,all ! 把所选择单元材料号改为3 #lQbMuR
allsel PG9won5_
solve !进行求解 bCk_ZA
save,F1.2,db !把F=1.2时求解结果保存 |H5){ 2V>K
!边坡在强度折减系数F=1.4时求解 9U6y<X
finish -8qLshQ
/solu StWDNAf)
allsel MC/$:PV
asel,s,area,,2,3,1 !选择面积A2、A3 } 1w[ G;$
MPCHG,4,all !把所选择单元材料号改为4 =/J4(# Xb
allsel kf0zL3|
solve !进行求解 RMDs~
save,F1.4,db !把F=1.4时求解结果保存 a=gT GG"9
!边坡在强度折减系数F=1.6时求解 3 S:}fPR
finish JGSeu =)
/solu fyx-VX u
allsel WUAjb ,eo
asel,s,area,,2,3,1 !选择面积A2、A3 " j3Yu4_ks
MPCHG,5,all !把所选择单元材料号改为5 bOux8OHt*
allsel [ X)+(-J
solve 3j0/&ON
save,F1.6,db !把F=1.6时求解结果保存 b]~X
U
!边坡在强度折减系数F=1.8时求解 Gp$[u4-6M6
finish e}f!zA
/solu |]DZc/
allsel b#P ,
asel,s,area,,2,3,1 !选择面积A2、A3 n<\^&_ a
MPCHG,6,all !把所选择单元材料号改为6 AZorz Q]s
allsel (iGk]Rtzt
solve W+a>*#*
save,F1.8,db !把F=1.8时求解结果保存 A2 Rr*e
!边坡在强度折减系数F=2.0时求解 dnNC
=
siY
finish Uk-^n~y
/solu !li Q;R&
allsel H]"Z_n_
asel,s,area,,2,3,1 !选择面积A2、A3 6*k Y7
MPCHG,7,all !把所选择单元材料号改为7 )& %X
AW{
allsel I5F oh|)
solve f `p`c*
save,F2.0,db !把F=2.0时求解结果保存 O\XN/R3
!边坡在强度折减系数F=2.2时求解 +uKlg#wqc
finish &B#H gWud
/solu ; _of'
allsel 9Q^>.^~ ^
asel,s,area,,2,3,1 !选择面积A2、A3 rcq(p (!
MPCHG,8,all !把所选择单元材料号改为8 ;QW6Tgt11
allsel "X>Z!>
solve :Py/d6KK
save,F2.2,db !把F=2.2时求解结果保存 nBaY|
!边坡在强度折减系数F=2.4时求解 B~lrd#qC
finish ai`fP{WlX
/solu ??aOr*%
allsel ]Tf.KUm
asel,s,area,,2,3,1 !选择面积A2、A3 =! 9+f
MPCHG,9,all !把所选择单元材料号改为9 so8isDC'9
allsel 6fY(u7m|p
solve aoc o'BR F
save,F2.4,db !把F=2.4时求解结果保存 y{s?]hLk
!边坡在强度折减系数F=2.6时求解 N3MMxm_u
finish ' >`?T}a,
/solu b@Ik
c<
allsel ^EjZ.#2l;
asel,s,area,,2,3,1 !选择面积A2、A3 E[Tz%x=P
MPCHG,10,all !把所选择单元材料号改为10 8K] fw{-$L
allsel e~ W35Y>A
solve 7sC$hm]
save,F2.6,db !把F=2.6时求解结果保存 :'f#0 ox
!边坡在强度折减系数F=2.8时求解 +yVz )
X
finish ^uM y|d
/solu =KmjCz:
allsel Dsm_T1X
asel,s,area,,2,3,1 !选择面积A2、A3 3 {\b/NL$
MPCHG,11,all !把所选择单元材料号改为11 T#|Qexz6 @
allsel fR~_5 pt7
solve jL I (Z
save,F2.8,db !把F=2.8时求解结果保存 ;B^G<
!边坡在强度折减系数F=3.0时求解 i t ,i^32|
finish ?2
u_E "
/solu gJ+MoAM"
allsel gPT-zul
asel,s,area,,2,3,1 !选择面积A2、A3 G=jdb@V/?
MPCHG,12,all !把所选择单元材料号改为12 8,a&i:C
allsel yj `b-^$?
solve c` #E#
save,F3.0,db !把F=3.0时求解结果保存 (NFq/w%
(5）后处理 b DeHU$
/post1 ！进入后处理 !\}Dxt
!边坡在强度折减系数F=1时结果分析 ?#< Fxme
Resume,'F1','db' !读入边坡在强度折减系数F=1时
PQa {5"
set,1,last !读入后一个子步 ^IjKT
pldisp,1 !绘制边坡模型变形图 ]sd|u[:k
plnsol,u,x !绘制边坡模型水平方向位移云图 yXh=~:1~
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 D}SRr, 4v
!边坡在强度折减系数F=1.2时结果分析 gLsl/G
Resume,'F1.2','db' !读入边坡在强度折减系数F=1.2时 5 4JZEc
set,1,last !读入后一个子步 WQ \' z?P
pldisp,1 !绘制边坡模型变形图 0<>I\UN0b
plnsol,u,x !绘制边坡模型水平方向位移云图 d0Py[37V
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 X{| 1E85fl
!边坡在强度折减系数F=1.4时结果分析 o*_[3{FU
Resume,'F1.4','db' !读入边坡在强度折减系数F=1.4时 e-ljwCD
set,1,last !读入后一个子步 3=o3VGZP
pldisp,1 !绘制边坡模型变形图 glH&v8
plnsol,u,x !绘制边坡模型水平方向位移云图 `4'v)!?
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 MSUkCWt!
!边坡在强度折减系数F=1.6时结果分析 UAFl+d!
Resume,'F1.6','db' !读入边坡在强度折减系数F=1.6时 yc8iT`
set,1,last !读入后一个子步 ]c tUl #j
pldisp,1 !绘制边坡模型变形图 <*V%!pwIG
plnsol,u,x !绘制边坡模型水平方向位移云图 L p 5LRw
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 >)R7*^m{'
!边坡在强度折减系数F=1.8时结果分析 2o] V q
Resume,'F1.8','db' !读入边坡在强度折减系数F=1.8时 .5 p"o- :D
set,1,last !读入后一个子步 M9.jJf
pldisp,1 !绘制边坡模型变形图 Om0Z\GP=
plnsol,u,x !绘制边坡模型水平方向位移云图 T\:Vu{|
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 |#D3~au
!边坡在强度折减系数F=2.0时结果分析 "o*(i7T =n
Resume,'F2.0','db' !读入边坡在强度折减系数F=2.0时 ` {c %d
set,1,last !读入后一个子步 a9j
f7r1
pldisp,1 !绘制边坡模型变形图 \l?\%aqm
plnsol,u,x !绘制边坡模型水平方向位移云图 g>12!2}
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 R+Dx#Wn I
!边坡在强度折减系数F=2.2时结果分析 Jt\?,~,
Resume,'F2.2','db' !读入边坡在强度折减系数F=2.2时 1\z5[
_
set,1,last !读入后一个子步 mDEO$:A
pldisp,1 !绘制边坡模型变形图 )[|TxXz
d
plnsol,u,x !绘制边坡模型水平方向位移云图 8~YhT]R=
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 2.{:PM4Z4
!边坡在强度折减系数F=2.4时结果分析 '`p#%I@
Resume,'F2.4','db' !读入边坡在强度折减系数F=2.4时 kO{A]LnAH
set,1,last !读入后一个子步 tV%:sk^d
pldisp,1 !绘制边坡模型变形图 7m~.V[l1
plnsol,u,x !绘制边坡模型水平方向位移云图 ;{[&&qMwU
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 [ G>8N5@*
!边坡在强度折减系数F=2.6时结果分析 .t*MGUg
Resume,'F2.6','db' !读入边坡在强度折减系数F=2.6时 . :A9*,
set,1,last !读入后一个子步 es 8%JTi
pldisp,1 !绘制边坡模型变形图 Eyv%"+>
plnsol,u,x !绘制边坡模型水平方向位移云图 sgu#`@o
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 |=fa`8m G
!边坡在强度折减系数F=2.8时结果分析 tGDsZ;3Yr
Resume,'F2.8','db' !读入边坡在强度折减系数F=2.8时 L 'B=
=#
set,1,last !读入后一个子步 ]&w8"q
pldisp,1 !绘制边坡模型变形图 89e<,f`h
plnsol,u,x !绘制边坡模型水平方向位移云图 o *)>aw
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 7TD%vhbiwi
!边坡在强度折减系数F=3.0时结果分析 F6J,:
Resume,'F3.0','db' !读入边坡在强度折减系数F=3.0时 ;) ].Dj9
set,1,last !读入后一个子步 LI1OocY.]
pldisp,1 !绘制边坡模型变形图 <Dojl
#
plnsol,u,x !绘制边坡模型水平方向位移云图 luXcr
H+w
plnsol,eppl,eqv !绘制边坡模型塑性应变云图
