各位朋友： vS$_H<;P
大家好！我通过ansys用有限元强度折减法分析了一个边坡的稳定性，结果不同的工况下得出的塑性区是一样的，我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激！附命令流和塑性区如下：(1)创建物理环境 \ &S-lsLY
/COM, Structural ! 指定结构分析 _Ka6! 9
/TITLE, Slope stability Analysis ! 定义工作标题 '"/Yk=EmlU
/FILNAM,support,1 ! 定义工作文件名 keYvscRBI
(2)建立模型 bim
82<F
!进入前处理器 P}k p_l27
/PREP7 dmkd.aP4
!定义单元类型 g;Lk 'Ky6
ET,1,PLANE82 !定义边坡围岩单元 10 D6fkjf
KEYOPT,1,3,2 !求解类型为平面应变 O0Vtvbj
!定义材料属性 ]U8VU
!边坡围岩1材料属性 :172I1|7
MP,EX,1,3E10 |Qe #[Q7
MP,PRXY,1,0.25 U(jZf{ `Mz
MP,DENS,1,2500 !定义塑性模型 #cRw0bn:
TB,DP,1 !定义数据表1 =(^-s Jk
TBDATA,1,9E5,42.8 !给数据表输入数据 3i c6!T#t"
!边坡围岩2材料属性 !采用弹性模型 wCwJ#-z.=
MP,EX,2,3.2E10 ！定义弹性模量 w>q:&Q
MP,PRXY,2,0.24 !定义泊松比 GkT:7`|C
MP,DENS,2,2700 !定义密度 dD Zds
k+!
!取强度折减系数F=1.2时边坡围岩材料属性 'K9 {xI@N
MP,EX,3,3E10 P})Iwk|Z
MP,PRXY,3,0.25 O&Ws*k
MP,DENS,3,2500 Hik :Sqpox
TB,DP,3 !定义数据表3 \
M_}V[1+
TBDATA,1,7.5E5,37.7 !给数据表输入数据 *61+Fzr
!取强度折减系数F=1.4时边坡围岩材料属性 EMV< PshW=
MP,EX,4,3E10 [=
G VK
MP,PRXY,4,0.25 !定义泊松比 Y<vsMf_U
MP,DENS,4,2500 opKk#40
TB,DP,4 !定义数据表4 #00k7y>OyD
TBDATA,1,6.4E5,33.5 !给数据表输入数据 P9\!JH!
!取强度折减系数F=1.6时边坡围岩材料属性 6U`<+[K7
MP,EX,5,3E10 EP|OKXRltA
MP,PRXY,5,0.25 Y)-)owx7
MP,DENS,5,2500 ！定义密度 dHtbl\6
TB,DP,5 !定义数据表5 ZDI%?.U
TBDATA,1,5.6E5,30 !给数据表输入数据 Wx i|(}
!取强度折减系数F=1.8时边坡围岩材料属性 oB m^RHTZ
MP,EX,6,3E10 ！定义弹性模量 *,UD&N_)*6
MP,PRXY,6,0.25 1%;o-F@
MP,DENS,6,2500 )< kI d4E
TB,DP,6 !定义数据表6 ^zMME*G
TBDATA,1,5.0E5,27.2 L6xLD X7y
!取强度折减系数F=2.0时边坡围岩材料属性 *7g gw[~
MP,EX,7,3E10 qJQ!e
MP,PRXY,7,0.25 !定义泊松比 wQ4 IQ!
MP,DENS,7,2500 } +@H&}u
TB,DP,7 !定义数据表7 XM+o e0:[
TBDATA,1,4.5E5,24.8 "K"]/3`k-
!取强度折减系数F=2.2时边坡围岩材料属性 JVoW*uA
MP,EX,8,3E10 _msV3JBr
MP,PRXY,8,0.25 Z ]1=nSv
MP,DENS,8,2500 @ -~
)M_
TB,DP,8 !定义数据表8 7.r}9 8V
TBDATA,1,4.09E5,22.8 Ls<.&3X2
!取强度折减系数F=2.4时边坡围岩材料属性 ~1NK@=7T
MP,EX,9,3E10 C {i9~80n
MP,PRXY,9,0.25 me#?1 r
MP,DENS,9,2500 g `)5g5
TB,DP,9 !定义数据表9 x}`]9XQ
TBDATA,1,3.6E5,21.1 !给数据表输入数据 [Ipg",Su;f
!取强度折减系数F=2.6时边坡围岩材料属性 ?9_RI(a.}
MP,EX,10,3E10 ！定义弹性模量 5i+0GN3nd
MP,PRXY,10,0.25 !定义泊松比 j=%^CRum
MP,DENS,10,2500 ！定义密度 )]=1W
TB,DP,10 !定义数据表10 ;Jn"^zT
TBDATA,1,3.46E5,19.6 !给数据表10输入数据 -7^A_!.
!取强度折减系数F=2.8时边坡围岩材料属性 \8~P3M":c
MP,EX,11,3E10 ！定义弹性模量 g=*`6@_=
MP,PRXY,11,0.25 !定义泊松比 8 f n 7!
MP,DENS,11,2500 fI/ ?2ZH
TB,DP,11 !定义数据表11 T[z]~MJL
TBDATA,1,3.2E5,18.2 !给数据表11输入数据 gbf-3KSp^
!取强度折减系数F=3.0时边坡围岩材料属性 (|6Y1``
MP,EX,12,3E10 ！定义弹性模量 >d`XR" _e
MP,PRXY,12,0.25 bN ,>,hj
MP,DENS,12,2500 ！定义密度 _?vh#6F
TB,DP,12 !定义数据表12 M*r/TT
TBDATA,1,3.0E5,17.1 !给数据表12输入数据 [|=#~(yYQ
!建立几何模型 d ,Fj|}S
!创建边坡线模型线 a w0;
K,1,,,, ！ 创建关键点1 92@/8,[
K,2,-800,,, sG/mmZHYzr
K,3,-800,-800,, pzo9?/-
K,4,-800,-1200,,
g 1je':
K,5,1200,-1200,, ~\4B 1n7
K,6,1200,-800,, ^1a/)Be{_
K,7,1200,,, G2bZl%
,D
K,8,1200,400,, X>[x7t:
K,9,400,400,, X#qm wcF
l,1,2 !连接关键点1、2生成直线L1 Nh I&wl
l,2,3 6I$laHx?
l,3,4
wkn r^A
l,4,5 '`-W!g[
>
l,5,6 `]U u` b
l,6,7 ^J?ExMu
l,7,8 2k
}:)] m
l,8,9 5 /v,|
l,9,1 (1
"unP-
l,1,7 K tT.WHr(m
l,3,6 m3apeIEi[
!创建边坡面模型 fp jy[$8
al,3,4,5,11 !由线L3、L4、L5、L11生成一个面A1 ~ wMdk9RQ
al,1,2,11,6,10 !由线L1、L2、L1、L6、L10生成一个面A2 uhyw?#f
al,7,8,9,10 !由线L7、L8、L9、L10生成一个面A3 [j6EzMN
/PNUM, area, on !打开面号开关 ho1Mo
!划分网格生成有限元模型 X{n7)kgL
! 划分边坡围岩2网格 pjC2jlwm*
!设置网格份数 m])!'Pa( =
lsel,s,,,3,5,2 !选择线L3和L5 N1 lhlw6
lesize,all,,,5 !把所选择线分为5段 Dq#/Uw#
lsel,s,,,4,11,7 !选择线L4和L11 WIytgM
lesize,all,,,26 !把所选择线分为26段 SIJ:[=5!7
mat,2 !给边坡围岩2赋予2号材料特性 p.r \ |
type,1 !采用单元类型1 thoAEG 80
mshkey,0 !设定自由网格划分 z7[TgL7
mshape,0 !设定四边形网格划分 Yy 4EM
amesh,1 ！划分面积A1 -BY'E$]4
! 划分边坡围岩1网格 z (N3oBW
!设置网格份数 *%G$[=
lsel,s,,,2,6,4 !选择线L2和L6 vybQ}dscn
lsel,a,,,1 !附加选择线L1 KGMX >t'
lesize,all,,,10 !把所选择线分为10段 6GxQ<
lsel,s,,,8,10,2 !选择线L8和L10 R^ }}-Dv r
lesize,all,,,16 !把所选择线分为16段 O*hd@2h d
lsel,s,,,7,9,2 6^W6As0
lesize,all,,,12 Cm~h\+"
mat,1 !给边坡围岩1赋予1号材料特性 zW`Zmt\T2
type,1 9;Z 2.P"w
mshkey,0 !设定自由网格划分 +3HukoR(
mshape,0 !设定四边形网格划分 R&w2y$
amesh,2 !划分面积A2 _(CuuP$`I
amesh,3 ！划分面积A3 $jt UQ1
allsel oYYns%r}{
!保存网格模型 |4UU`J9M
save,Slope-grid.db i\P?Y(-{
(3)施加约束和荷载 g$&uD
!给边坡模型施加约束 ) V!9/d
!边坡两侧施加X方向约束 8 1 \$X
nsel,s,loc,x,-800 !选择X=-800线上所有节点 u/ Os
nsel,a,loc,x,1200 !选择X=1200线上所有节点
5^Gv! XW
d,all,ux !对所选择节点约束X方向位移 R&z)
allsel h)%}O.ueB
!边坡底部施加约束 87/!u]q
nsel,s,loc,y,-1200 !选择Y=--1200线上所有节点 u9WQ0.
d,all,uy !对所选择节点约束X、Y方向位移 :<-,[(@bR
d,all,ux (nhv#&Fd+
!施加重力加速度 =)56]ki}
acel,,9.8 r ]XXN2[jO
(4)求解 ?e? mg
/solu .YvE
!求解设置
b>N) H
antype,static !设定为静力求解 -IB f;"8f
nsubst,100 !设定最大子步数为100 aR }|^ex
pred,on !打开时间步长预测器 Rp~#zt9:
nropt,full !设定牛顿-拉普森选项 vYh_<Rp5
nlgeom,on !打开大位移效果 "W_C%elg
lnsrch,on !打开线性搜索 s]=s2.=
outres,all,all !输出所有项 \z}/=Qgc
cnvtol,f,,0.005,2,0.5 !力收敛准则设定 sGNVZx
cnvtol,u,,0.05,2,1 !位移收敛准则设定 ZE#f{qF(
!边坡在强度折减系数F=1时求解 ,/. U'{
allsel 9m2_zfO[ w
solve !进行求解 aBPaC=g{HO
save,F1,db !把F=1时求解结果保存 ]J^ 9iDTTA
!边坡在强度折减系数F=1.2时求解 b \/:-][
finish | v?
pS
/solu t 6.hg3Y
allsel RB5fn+FiZ
asel,s,area,,2,3,1 !选择面积A2、A3 uV]4C^k;`[
MPCHG,3,all ! 把所选择单元材料号改为3 JHY0 J
&4s
allsel zj!&12w%3
solve !进行求解 %G%##wv:
save,F1.2,db !把F=1.2时求解结果保存 f%LzWXA
!边坡在强度折减系数F=1.4时求解 ;4jRsirx9
finish :Zy7h7P,lT
/solu cD-.thHO
allsel <1(:W[M
asel,s,area,,2,3,1 !选择面积A2、A3 U
D9 &k^
MPCHG,4,all !把所选择单元材料号改为4 T-Od|T@[
allsel } z4=3 '
solve !进行求解 #;bpxz1lR9
save,F1.4,db !把F=1.4时求解结果保存 dzY B0vut@
!边坡在强度折减系数F=1.6时求解 r7tN(2;5
finish ?^z!yD \
/solu :kycIM]s
allsel Mis t,H7
asel,s,area,,2,3,1 !选择面积A2、A3 ICNS+KsI
MPCHG,5,all !把所选择单元材料号改为5 Rd5ni2-nve
allsel CroI,=a&,
solve ETP}mo
save,F1.6,db !把F=1.6时求解结果保存 (>qX>
!边坡在强度折减系数F=1.8时求解 Az. k6)~
finish RUTlwTdv
/solu m1 7 8S3
allsel K
7)1wiEj
asel,s,area,,2,3,1 !选择面积A2、A3 KlgPDV9mg
MPCHG,6,all !把所选择单元材料号改为6 *|n::9
allsel ; 6W lu3I
solve %K^l]tWa@
save,F1.8,db !把F=1.8时求解结果保存 2e}${NZN
!边坡在强度折减系数F=2.0时求解 8?S32Gdu
finish QM I&?Q:=
/solu } d / 5_X
allsel +w Oa
asel,s,area,,2,3,1 !选择面积A2、A3 Ti/t\'6
MPCHG,7,all !把所选择单元材料号改为7 fI%+
allsel DPvM|n`TW
solve QDlEby m
save,F2.0,db !把F=2.0时求解结果保存 !FweXFl
!边坡在强度折减系数F=2.2时求解 :sRV]!Iw
finish `BKo`@
/solu 7%Zl^c>q
allsel LOyL:~$
asel,s,area,,2,3,1 !选择面积A2、A3 2 Mc/ah
MPCHG,8,all !把所选择单元材料号改为8 Sf>R7.lpP
allsel \1 4"B gj1
solve \F%5T RoC
save,F2.2,db !把F=2.2时求解结果保存 V\zf yH\~
!边坡在强度折减系数F=2.4时求解 i^ 9Pi P|U
finish .BGM1ph}~
/solu 37 )Dx
allsel *F+t`<2
asel,s,area,,2,3,1 !选择面积A2、A3 -bHQy:
MPCHG,9,all !把所选择单元材料号改为9 >_QC_UX>4i
allsel q u[ ~#
solve -UaUFJa8K&
save,F2.4,db !把F=2.4时求解结果保存 7<[p1C*B
!边坡在强度折减系数F=2.6时求解 @Md%gEh;&
finish &H!3]
/solu }7CMXw
[
allsel !RLg[_'
asel,s,area,,2,3,1 !选择面积A2、A3 6#XB'PR2p
MPCHG,10,all !把所选择单元材料号改为10 ]!J 6S.@#+
allsel 27*u^N*z@
solve 0 uzm@'^
save,F2.6,db !把F=2.6时求解结果保存 (#* 7LdZ
!边坡在强度折减系数F=2.8时求解 Efi@hdEV
finish o.t$hv|
/solu &LE ,.Q34
allsel <,d .`0:y
asel,s,area,,2,3,1 !选择面积A2、A3 ;bVC7D~~4w
MPCHG,11,all !把所选择单元材料号改为11 @cPb*
allsel I0 a,mO;m
solve >N>WOLbb7(
save,F2.8,db !把F=2.8时求解结果保存 SU; PmG4
!边坡在强度折减系数F=3.0时求解 v_f8zk
finish gB@Xi*
/solu '3Yci(t+
allsel ;cl\$TDL
asel,s,area,,2,3,1 !选择面积A2、A3 Uw^`_\si
MPCHG,12,all !把所选择单元材料号改为12 ~dfc
allsel `_sc_Y|C!
solve >*H>'O4
save,F3.0,db !把F=3.0时求解结果保存 }}XYV eI
(5）后处理 ?Y2ZqI
/post1 ！进入后处理 |du@iA]dP
!边坡在强度折减系数F=1时结果分析 v4!z B9d
Resume,'F1','db' !读入边坡在强度折减系数F=1时 zVe@`gc
set,1,last !读入后一个子步 N2uxiXpQZ=
pldisp,1 !绘制边坡模型变形图 {dJC3/ Rf
plnsol,u,x !绘制边坡模型水平方向位移云图 T;jp2 #
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 /W'GX n
!边坡在强度折减系数F=1.2时结果分析 U'zW; Lt
Resume,'F1.2','db' !读入边坡在强度折减系数F=1.2时 mD- qJ6AM
set,1,last !读入后一个子步 ?P[uf
pldisp,1 !绘制边坡模型变形图 S (] (C
plnsol,u,x !绘制边坡模型水平方向位移云图 ^,' )1r,
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 'EJ8)2
!边坡在强度折减系数F=1.4时结果分析 ~U`o ew
Resume,'F1.4','db' !读入边坡在强度折减系数F=1.4时 p+U}oC
set,1,last !读入后一个子步 Tx?,]c,(u
pldisp,1 !绘制边坡模型变形图 SCE5|3j
plnsol,u,x !绘制边坡模型水平方向位移云图 n'=-bj`
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 R *IO%9O
!边坡在强度折减系数F=1.6时结果分析 a'f"Zdh%w
Resume,'F1.6','db' !读入边坡在强度折减系数F=1.6时 Rk.YnA_J6
set,1,last !读入后一个子步 6 lO]V=+
pldisp,1 !绘制边坡模型变形图 g1VdP[ Y#
plnsol,u,x !绘制边坡模型水平方向位移云图 VO3& ! uOd
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 + 9G
GC
!边坡在强度折减系数F=1.8时结果分析 `E>HpRcxD
Resume,'F1.8','db' !读入边坡在强度折减系数F=1.8时 ![3#([>4>
set,1,last !读入后一个子步 BL<.u
pldisp,1 !绘制边坡模型变形图 xaSvjc\
plnsol,u,x !绘制边坡模型水平方向位移云图 zQ9"i
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 9K~2!<
!边坡在强度折减系数F=2.0时结果分析 xy$73K6
Resume,'F2.0','db' !读入边坡在强度折减系数F=2.0时 QlJ
cj+_h
set,1,last !读入后一个子步 \bqIe}3V7
pldisp,1 !绘制边坡模型变形图 W-q2|NK
plnsol,u,x !绘制边坡模型水平方向位移云图 b {<qt} )
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 D'Tb=
!边坡在强度折减系数F=2.2时结果分析 o9ZHa
Resume,'F2.2','db' !读入边坡在强度折减系数F=2.2时 B1 T:c4:N
set,1,last !读入后一个子步 24PEt%2
pldisp,1 !绘制边坡模型变形图 k@U8K(:x
plnsol,u,x !绘制边坡模型水平方向位移云图 F,Q\_H##x4
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 5G=< 2;
!边坡在强度折减系数F=2.4时结果分析 mZ
39 s
Resume,'F2.4','db' !读入边坡在强度折减系数F=2.4时 q 65KxOf`
set,1,last !读入后一个子步 K:
g_M
pldisp,1 !绘制边坡模型变形图 BDSZ '
plnsol,u,x !绘制边坡模型水平方向位移云图 \$Y Kw0K
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 )orVI5ti
!边坡在强度折减系数F=2.6时结果分析 !s pp*Q)#\
Resume,'F2.6','db' !读入边坡在强度折减系数F=2.6时 w2!G"oD
set,1,last !读入后一个子步 gR v5l3k
pldisp,1 !绘制边坡模型变形图 sm5\> L3V
plnsol,u,x !绘制边坡模型水平方向位移云图 $X8(OS5d'
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 C( 8i0(1
!边坡在强度折减系数F=2.8时结果分析 n 1MZHa,
Resume,'F2.8','db' !读入边坡在强度折减系数F=2.8时 =r"8 J5 [f
set,1,last !读入后一个子步 &C<K|F!j!
pldisp,1 !绘制边坡模型变形图 D7|[:``
plnsol,u,x !绘制边坡模型水平方向位移云图 <+ UEM ~)
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 bh|M]*Pq
!边坡在强度折减系数F=3.0时结果分析 0X@! i3eu
Resume,'F3.0','db' !读入边坡在强度折减系数F=3.0时 j78WPG
set,1,last !读入后一个子步 ~$n4Yuu2[
pldisp,1 !绘制边坡模型变形图 > Jh*S`e
plnsol,u,x !绘制边坡模型水平方向位移云图 ifj %!*
plnsol,eppl,eqv !绘制边坡模型塑性应变云图
