各位朋友： TocqoYX{{
大家好！我通过ansys用有限元强度折减法分析了一个边坡的稳定性，结果不同的工况下得出的塑性区是一样的，我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激！附命令流和塑性区如下：(1)创建物理环境 Oqyh{q%]
/COM, Structural ! 指定结构分析 <[Vr(.A
/TITLE, Slope stability Analysis ! 定义工作标题 9Bn
dbS i
/FILNAM,support,1 ! 定义工作文件名 HhO $`YZ%>
(2)建立模型 `#hy'S:e
!进入前处理器 `?:'_K i
/PREP7 0)Z7U$
!定义单元类型 I`|>'$E[ r
ET,1,PLANE82 !定义边坡围岩单元 Y+E@afsKs
KEYOPT,1,3,2 !求解类型为平面应变 -{?Rq'H
!定义材料属性 &nq[Vy0kO4
!边坡围岩1材料属性 ^,3 >}PU
MP,EX,1,3E10
J0`?g6aY
MP,PRXY,1,0.25 8^T' a^Wt
MP,DENS,1,2500 !定义塑性模型 fN9hBC@
TB,DP,1 !定义数据表1
=~)n,5
TBDATA,1,9E5,42.8 !给数据表输入数据 ;Zw28!#Rt
!边坡围岩2材料属性 !采用弹性模型 EpiagCS
MP,EX,2,3.2E10 ！定义弹性模量 &8N\
6K=
MP,PRXY,2,0.24 !定义泊松比 cWi2Sls
MP,DENS,2,2700 !定义密度 {VF p fo
!取强度折减系数F=1.2时边坡围岩材料属性 tU(6%zvR
MP,EX,3,3E10 @U}UC G7+
MP,PRXY,3,0.25 x-"7{@lz
MP,DENS,3,2500 -`( :L[
TB,DP,3 !定义数据表3 I^o^@ C
TBDATA,1,7.5E5,37.7 !给数据表输入数据 *oF { R^
!取强度折减系数F=1.4时边坡围岩材料属性 V1+IqOXAIp
MP,EX,4,3E10 D^pAf/ek@i
MP,PRXY,4,0.25 !定义泊松比 >[#4Pb7_ Y
MP,DENS,4,2500 w v1R
]3}
TB,DP,4 !定义数据表4 Sdn]
f4
TBDATA,1,6.4E5,33.5 !给数据表输入数据 `bi5 #xR
!取强度折减系数F=1.6时边坡围岩材料属性 /w| YNDA]j
MP,EX,5,3E10 ji1A>jepF
MP,PRXY,5,0.25 ,Vogo5~X
MP,DENS,5,2500 ！定义密度 hX^XtIC=
TB,DP,5 !定义数据表5 [wEx jLW
TBDATA,1,5.6E5,30 !给数据表输入数据 _*+M'3&=
!取强度折减系数F=1.8时边坡围岩材料属性 Lczcz"t
MP,EX,6,3E10 ！定义弹性模量 - na]P3 s
MP,PRXY,6,0.25 @8 pRI S"V
MP,DENS,6,2500 V{8mx70
TB,DP,6 !定义数据表6 DP?gozm
TBDATA,1,5.0E5,27.2 ^t| %!r
G
!取强度折减系数F=2.0时边坡围岩材料属性 d00r&Mc
MP,EX,7,3E10 %c&h:7);
MP,PRXY,7,0.25 !定义泊松比 ^Tl|v'
MP,DENS,7,2500 =[P ||
TB,DP,7 !定义数据表7 }5 A?WH_
TBDATA,1,4.5E5,24.8 {5,C W
!取强度折减系数F=2.2时边坡围岩材料属性 2Bx\nLf/
K
MP,EX,8,3E10 {BF$N#7
MP,PRXY,8,0.25 V
Bg\)r[
MP,DENS,8,2500 _H-Lt{k
TB,DP,8 !定义数据表8 MQI6e".
TBDATA,1,4.09E5,22.8 #PiW\Tq
!取强度折减系数F=2.4时边坡围岩材料属性 u+ -}|
MP,EX,9,3E10 !$O +M#
MP,PRXY,9,0.25 v"^G9u
MP,DENS,9,2500 #e=[W) )
TB,DP,9 !定义数据表9 {3l] / X3
TBDATA,1,3.6E5,21.1 !给数据表输入数据 ?&[`=ZVn
!取强度折减系数F=2.6时边坡围岩材料属性 rT x] %{
MP,EX,10,3E10 ！定义弹性模量 m49)c K?
MP,PRXY,10,0.25 !定义泊松比 f>Ge
Em~
MP,DENS,10,2500 ！定义密度 ^mH^cP?/
TB,DP,10 !定义数据表10 |d{4 _o90
TBDATA,1,3.46E5,19.6 !给数据表10输入数据 V%"aU}
!取强度折减系数F=2.8时边坡围岩材料属性 oR5 'g7?
MP,EX,11,3E10 ！定义弹性模量 Y M5;m PR
MP,PRXY,11,0.25 !定义泊松比 J"|o g|Tz
MP,DENS,11,2500 Hi
)n]OE
TB,DP,11 !定义数据表11 EayZ*e ]
TBDATA,1,3.2E5,18.2 !给数据表11输入数据 q$\KE4v "
!取强度折减系数F=3.0时边坡围岩材料属性 1"i/*}M
MP,EX,12,3E10 ！定义弹性模量 vzfWPjpKW
MP,PRXY,12,0.25 fi1tF/ `
MP,DENS,12,2500 ！定义密度 %kF6y_h`
TB,DP,12 !定义数据表12 tY V mB:l
TBDATA,1,3.0E5,17.1 !给数据表12输入数据 `}1 8A.K
!建立几何模型 C}7Sh6
!创建边坡线模型线 9Y0w
SOSW
K,1,,,, ！ 创建关键点1 }e1f kjWk
K,2,-800,,, 59GS:
K,3,-800,-800,, aK 'BC>uFI
K,4,-800,-1200,, ^W;\faG
K,5,1200,-1200,, Wq]^1g_
K,6,1200,-800,, "I
u3&mc
K,7,1200,,, F6yMk%
K,8,1200,400,, oB+drDp8U
K,9,400,400,, }s? 9Hnqa
l,1,2 !连接关键点1、2生成直线L1 Y\s@'UoVN
l,2,3 6?r}bs6Msx
l,3,4 '};pu;GA7
l,4,5 @2V#bK
l,5,6 L\o-z NY
l,6,7 yj-BLR5
l,7,8 X"*pt5B6`
l,8,9 y%f'7YZ4
l,9,1 L~>pSP^a
l,1,7 |
X! d*4
l,3,6 86oa>#opU
!创建边坡面模型 JT!-Q!O}O
al,3,4,5,11 !由线L3、L4、L5、L11生成一个面A1 rQCj^=cf;~
al,1,2,11,6,10 !由线L1、L2、L1、L6、L10生成一个面A2 r"c<15g2'
al,7,8,9,10 !由线L7、L8、L9、L10生成一个面A3 i ,ga2{GnM
/PNUM, area, on !打开面号开关 +hGr2%*0f
!划分网格生成有限元模型 Db@$'
! 划分边坡围岩2网格 ,-D3tleu`
!设置网格份数 z`XX[9 $qm
lsel,s,,,3,5,2 !选择线L3和L5 Kx7s
d i
lesize,all,,,5 !把所选择线分为5段 5*C#~gd& F
lsel,s,,,4,11,7 !选择线L4和L11 h<JV6h :8
lesize,all,,,26 !把所选择线分为26段 4'[/gMUkw
mat,2 !给边坡围岩2赋予2号材料特性 PCtf&U
type,1 !采用单元类型1 uS;N&6;:
mshkey,0 !设定自由网格划分 (}
?")$.
mshape,0 !设定四边形网格划分 @!`x^Tzz
amesh,1 ！划分面积A1 M]
7#
! 划分边坡围岩1网格 ,Tjc\;~%
!设置网格份数 3,t3 \`=
lsel,s,,,2,6,4 !选择线L2和L6 <<4U :
lsel,a,,,1 !附加选择线L1 >We4F2?
lesize,all,,,10 !把所选择线分为10段 E0G"B' x
lsel,s,,,8,10,2 !选择线L8和L10 ;L r]w8 d
lesize,all,,,16 !把所选择线分为16段 {W,&jC
lsel,s,,,7,9,2 @p"m{
lesize,all,,,12 G*^4+^Vz?
mat,1 !给边坡围岩1赋予1号材料特性 ~qT5F)$B-
type,1 -Q9} gaH_
mshkey,0 !设定自由网格划分 EZypqe):/C
mshape,0 !设定四边形网格划分 J#ujI e
amesh,2 !划分面积A2 l$
^LY)i
amesh,3 ！划分面积A3 5o2w)<d!
allsel p?PK8GL
!保存网格模型 @~sJ
((G[5
save,Slope-grid.db bi~1d "j
(3)施加约束和荷载 SJ8Ax_9{q
!给边坡模型施加约束 Y`3V&8X
!边坡两侧施加X方向约束 @L0xU??"|
nsel,s,loc,x,-800 !选择X=-800线上所有节点 }LeizbU
nsel,a,loc,x,1200 !选择X=1200线上所有节点 m9M#)< @*
d,all,ux !对所选择节点约束X方向位移 ^ZP
$(a4
allsel ip8%9fG\>
!边坡底部施加约束 4j_\_:$w<
nsel,s,loc,y,-1200 !选择Y=--1200线上所有节点 B63pu X{u#
d,all,uy !对所选择节点约束X、Y方向位移 :J6 xY y$
d,all,ux 9?]69O
!施加重力加速度 6
);8z!+
acel,,9.8 2+s#5K&i
(4)求解 zl $mt '\y
/solu )xm[m vt
!求解设置 TMJq-u51
antype,static !设定为静力求解 ^pwT8Bp
nsubst,100 !设定最大子步数为100 (P'{A>aHl0
pred,on !打开时间步长预测器 #|8Ia:=s
nropt,full !设定牛顿-拉普森选项 QNWGUg4*&
nlgeom,on !打开大位移效果 > zA*W<g
lnsrch,on !打开线性搜索 [3 D*DyQt
outres,all,all !输出所有项 Y<qWG 8X
cnvtol,f,,0.005,2,0.5 !力收敛准则设定 '- X[T}
cnvtol,u,,0.05,2,1 !位移收敛准则设定 Gn&=<q :H
!边坡在强度折减系数F=1时求解 [8jIu&tJf
allsel q]Vxf!0*>
solve !进行求解 duQ , 6
save,F1,db !把F=1时求解结果保存 +cWo^ d.
!边坡在强度折减系数F=1.2时求解 18nT
Iz_
finish G9_7jX*
/solu 4g+Dp&U
allsel k~Z ;S QyN
asel,s,area,,2,3,1 !选择面积A2、A3 X rl# DN
MPCHG,3,all ! 把所选择单元材料号改为3 Y#[xX2z9
allsel Z2a~1BL
solve !进行求解 vB<9M-sa0
save,F1.2,db !把F=1.2时求解结果保存 ;x.5_Xw{.
!边坡在强度折减系数F=1.4时求解 ?xW,2S
finish QZ&
4W
/solu q\Q'9Rl0(
allsel ij5|P4Eka
asel,s,area,,2,3,1 !选择面积A2、A3 `:.a5
MPCHG,4,all !把所选择单元材料号改为4 +N}yqg E
allsel OMf w#
solve !进行求解 L.2!Q3&
save,F1.4,db !把F=1.4时求解结果保存 7?ICXhu9
!边坡在强度折减系数F=1.6时求解 I
>a K a
finish Ac P d(Pc
/solu >npTUOGL=n
allsel l=UXikx
asel,s,area,,2,3,1 !选择面积A2、A3 l#IN)">1
MPCHG,5,all !把所选择单元材料号改为5 zb3,2D+ P
allsel g o9tvK
solve |_ njN
save,F1.6,db !把F=1.6时求解结果保存 +\-cf,WkI
!边坡在强度折减系数F=1.8时求解 U0=: `G2l
finish %/oeV;D
/solu \AeM=K6q+D
allsel nT
:n>ja
asel,s,area,,2,3,1 !选择面积A2、A3 a(IU Ah*mO
MPCHG,6,all !把所选择单元材料号改为6 ',Y`\X
allsel 2Q%*`
vCuV
solve U4=m>Ty
save,F1.8,db !把F=1.8时求解结果保存 ~jJF&*)
!边坡在强度折减系数F=2.0时求解 .;F%k,!v
finish `ztp u
~?
/solu .WLwAL
allsel j-#h^3l1?
asel,s,area,,2,3,1 !选择面积A2、A3 #(}'G*
MPCHG,7,all !把所选择单元材料号改为7 'i <%kL@
allsel &$bcB]C\3
solve a$K6b5`>Rs
save,F2.0,db !把F=2.0时求解结果保存 ?{L5=X@$$
!边坡在强度折减系数F=2.2时求解 & LhQr-g
finish T ]x]hQ
/solu $px1D$F !
allsel v]\T&w%9
asel,s,area,,2,3,1 !选择面积A2、A3 l6-
n{zG
MPCHG,8,all !把所选择单元材料号改为8 _MW W
allsel (3WK2IM^
solve L 8J] X7
save,F2.2,db !把F=2.2时求解结果保存 3"Zc|Ck <?
!边坡在强度折减系数F=2.4时求解 s14D(:t(
finish :} -VLp4b
/solu rn]F97v@]
allsel y[ZVi5) ,
asel,s,area,,2,3,1 !选择面积A2、A3 vRb7=fXf
MPCHG,9,all !把所选择单元材料号改为9 tIk$4)ZAl
allsel 4C[kj
solve Q!5W x
save,F2.4,db !把F=2.4时求解结果保存 xb2j
|KY7
!边坡在强度折减系数F=2.6时求解 `(r0+Qx
finish 8=)9ZjfD
/solu !?" pnKb}
allsel YJ+l
\Wb}
asel,s,area,,2,3,1 !选择面积A2、A3 9a#Y
D;-p
MPCHG,10,all !把所选择单元材料号改为10 XVF!l> nE
allsel BCO (,k
solve h.\p+Qw.
save,F2.6,db !把F=2.6时求解结果保存 Y() ZM
!边坡在强度折减系数F=2.8时求解 Pv|sPIIB7
finish @O/"s~d-
/solu _hgG F9
allsel ~x[(1
asel,s,area,,2,3,1 !选择面积A2、A3 \!"3 yd
MPCHG,11,all !把所选择单元材料号改为11 F$ckW'V
allsel p+228K ;H
solve ;{Yr|
save,F2.8,db !把F=2.8时求解结果保存 5&.I9}[)j
!边坡在强度折减系数F=3.0时求解 cSj X/%*!m
finish CDY3 +!
/solu fZ(k"*\MZ
allsel X/bu z
asel,s,area,,2,3,1 !选择面积A2、A3 )Y`ybAD d3
MPCHG,12,all !把所选择单元材料号改为12 0+K<;5"63d
allsel vR
(nd
solve k(xB%>ns
save,F3.0,db !把F=3.0时求解结果保存 Kpkpr`:)]
(5）后处理 "EVf1iQ
/post1 ！进入后处理 WD5jO9Oai
!边坡在强度折减系数F=1时结果分析 R[kF( C&
Resume,'F1','db' !读入边坡在强度折减系数F=1时 #7"";"{ z|
set,1,last !读入后一个子步 nt@uV wfQ
pldisp,1 !绘制边坡模型变形图 3vkzN
plnsol,u,x !绘制边坡模型水平方向位移云图 \~r_S
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 ps^["3e
!边坡在强度折减系数F=1.2时结果分析 Ce~Pms]
Resume,'F1.2','db' !读入边坡在强度折减系数F=1.2时 `QZK W
set,1,last !读入后一个子步 s )Xz}QPK.
pldisp,1 !绘制边坡模型变形图 g][n1$%
plnsol,u,x !绘制边坡模型水平方向位移云图 b8 1cq,
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 %_G '#Bn<
!边坡在强度折减系数F=1.4时结果分析 qeC^e}h
Resume,'F1.4','db' !读入边坡在强度折减系数F=1.4时 Md0`/F:+2
set,1,last !读入后一个子步 ,4k3C#!. i
pldisp,1 !绘制边坡模型变形图 #%pY,AK:=
plnsol,u,x !绘制边坡模型水平方向位移云图 vs>Pd |p;
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 _ bRgr
!边坡在强度折减系数F=1.6时结果分析 nkz<t
Resume,'F1.6','db' !读入边坡在强度折减系数F=1.6时 aU/y>Y <k
set,1,last !读入后一个子步 Esm=sPW
pldisp,1 !绘制边坡模型变形图 b~{nS,_Rn
plnsol,u,x !绘制边坡模型水平方向位移云图 C`uL
4r
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 I+`>e*:@W
!边坡在强度折减系数F=1.8时结果分析 r*ziO#[
Resume,'F1.8','db' !读入边坡在强度折减系数F=1.8时 $h}w: AV:
set,1,last !读入后一个子步 ) NTpb
pldisp,1 !绘制边坡模型变形图 qW /&.
plnsol,u,x !绘制边坡模型水平方向位移云图 X5(oL
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 ti 3S'K0t
!边坡在强度折减系数F=2.0时结果分析 #bCUI*N"P
Resume,'F2.0','db' !读入边坡在强度折减系数F=2.0时 ^F- 2tc
set,1,last !读入后一个子步 .J=QWfqt
pldisp,1 !绘制边坡模型变形图
'*EK i
plnsol,u,x !绘制边坡模型水平方向位移云图 b ! `6s
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 UR(i_T&w
!边坡在强度折减系数F=2.2时结果分析 7+Nl)d:C J
Resume,'F2.2','db' !读入边坡在强度折减系数F=2.2时 #G fM^sK
set,1,last !读入后一个子步 x TqP`ljX
pldisp,1 !绘制边坡模型变形图 7jr+jNsowj
plnsol,u,x !绘制边坡模型水平方向位移云图 $GO F'
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 GnC s_[*&r
!边坡在强度折减系数F=2.4时结果分析 C >F5=&
Resume,'F2.4','db' !读入边坡在强度折减系数F=2.4时 I}|E_U1Qj
set,1,last !读入后一个子步 . yHK
pldisp,1 !绘制边坡模型变形图 @LY[kt6o
plnsol,u,x !绘制边坡模型水平方向位移云图 fF2] 7:
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 s) {VU2.ra
!边坡在强度折减系数F=2.6时结果分析 TH55@1W ,[
Resume,'F2.6','db' !读入边坡在强度折减系数F=2.6时 /7ShE-.5#
set,1,last !读入后一个子步 -=n!k^?lK
pldisp,1 !绘制边坡模型变形图 s2F[v:|Wq
plnsol,u,x !绘制边坡模型水平方向位移云图 : GFK
|
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 m E'HRv
!边坡在强度折减系数F=2.8时结果分析 s#aan e
Resume,'F2.8','db' !读入边坡在强度折减系数F=2.8时 x_ <,GE@
set,1,last !读入后一个子步 LoSrXK~0~J
pldisp,1 !绘制边坡模型变形图 O&#>i]*V
plnsol,u,x !绘制边坡模型水平方向位移云图 3Vk\iJ
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 sxdDI?W4
!边坡在强度折减系数F=3.0时结果分析 =L;g:hc<
Resume,'F3.0','db' !读入边坡在强度折减系数F=3.0时 eT?vZH[N
set,1,last !读入后一个子步 K,+z^{Hvh
pldisp,1 !绘制边坡模型变形图 4F<wa s/
plnsol,u,x !绘制边坡模型水平方向位移云图
Y=H_U$
plnsol,eppl,eqv !绘制边坡模型塑性应变云图
