各位朋友： (8s]2\/Ar
大家好！我通过ansys用有限元强度折减法分析了一个边坡的稳定性，结果不同的工况下得出的塑性区是一样的，我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激！附命令流和塑性区如下：(1)创建物理环境 5 TD "
/COM, Structural ! 指定结构分析 dXh@E 7
/TITLE, Slope stability Analysis ! 定义工作标题 N(?yOB4gt
/FILNAM,support,1 ! 定义工作文件名 qYFol# =%
(2)建立模型 %iI0JF*E z
!进入前处理器 GLb}_-|
/PREP7 LnJ/t(KV
!定义单元类型 ;G.m;5A
ET,1,PLANE82 !定义边坡围岩单元 DA
oOs}D
KEYOPT,1,3,2 !求解类型为平面应变 %PG::b
!定义材料属性 :):=KowI
!边坡围岩1材料属性 y( : hN)
MP,EX,1,3E10 FhB^E$r%
MP,PRXY,1,0.25 sBIqee'T
MP,DENS,1,2500 !定义塑性模型 Vgs( feGs
TB,DP,1 !定义数据表1 oHm U|
TBDATA,1,9E5,42.8 !给数据表输入数据 JF*JF Ob
!边坡围岩2材料属性 !采用弹性模型 x8T5aS
MP,EX,2,3.2E10 ！定义弹性模量 ^zGgvFf>
MP,PRXY,2,0.24 !定义泊松比 ]{OEU]I@
MP,DENS,2,2700 !定义密度 " 7!K'i
!取强度折减系数F=1.2时边坡围岩材料属性 k=[!{I
MP,EX,3,3E10 |}*k |
MP,PRXY,3,0.25 -[#Mx}%
MP,DENS,3,2500 {6<7 M
TB,DP,3 !定义数据表3 L%a ni}V
TBDATA,1,7.5E5,37.7 !给数据表输入数据 )o[ O%b
!取强度折减系数F=1.4时边坡围岩材料属性 tg~&kaz
MP,EX,4,3E10 G'(
%8\
MP,PRXY,4,0.25 !定义泊松比 66=6;77
MP,DENS,4,2500 6| #^4D)
TB,DP,4 !定义数据表4 (
$3j
TBDATA,1,6.4E5,33.5 !给数据表输入数据 vi6EI
wZG
!取强度折减系数F=1.6时边坡围岩材料属性 i:WHql"Kw_
MP,EX,5,3E10 wLD/#Hfi7
MP,PRXY,5,0.25 V/+r"le
MP,DENS,5,2500 ！定义密度 [;VNuF
TB,DP,5 !定义数据表5 1h\: Lj
TBDATA,1,5.6E5,30 !给数据表输入数据 (j u-r*0
!取强度折减系数F=1.8时边坡围岩材料属性 oKTIoTb
MP,EX,6,3E10 ！定义弹性模量 RR:m <9l
MP,PRXY,6,0.25 0D>~uNcT}
MP,DENS,6,2500 [pb X_
TB,DP,6 !定义数据表6 }H{{ @RU
TBDATA,1,5.0E5,27.2 DR9M 8E
!取强度折减系数F=2.0时边坡围岩材料属性 1vu4}%nD
MP,EX,7,3E10 M[_~7~4
MP,PRXY,7,0.25 !定义泊松比 QEM") (
MP,DENS,7,2500 xIF
z@9+k
TB,DP,7 !定义数据表7 9AJ!7J#v"
TBDATA,1,4.5E5,24.8 RlX;c!K
!取强度折减系数F=2.2时边坡围岩材料属性 gFJ& t^yL
MP,EX,8,3E10 K0]'v>AWr
MP,PRXY,8,0.25 -e%=Mpq.
MP,DENS,8,2500 w\; =3C`
TB,DP,8 !定义数据表8 IxCEE 5+`%
TBDATA,1,4.09E5,22.8 ?ZSG4La\
!取强度折减系数F=2.4时边坡围岩材料属性 .i/]1X*;r^
MP,EX,9,3E10 &a8#q v"l
MP,PRXY,9,0.25 (0W%Y Z!&
MP,DENS,9,2500 d@"eWvnlZ
TB,DP,9 !定义数据表9 `QR2!W70o3
TBDATA,1,3.6E5,21.1 !给数据表输入数据 -!MDYj +U
!取强度折减系数F=2.6时边坡围岩材料属性 N_L&!%s
MP,EX,10,3E10 ！定义弹性模量 8Es]WR5
^
MP,PRXY,10,0.25 !定义泊松比 !23W=N}82
MP,DENS,10,2500 ！定义密度 _]tR1T5e
TB,DP,10 !定义数据表10 }i/ &m&VU
TBDATA,1,3.46E5,19.6 !给数据表10输入数据 .jr1<LE
!取强度折减系数F=2.8时边坡围岩材料属性 %H- [u}s
MP,EX,11,3E10 ！定义弹性模量 Ta!.oC[
MP,PRXY,11,0.25 !定义泊松比 *|Re,cY
MP,DENS,11,2500 #hd<5+$U}l
TB,DP,11 !定义数据表11 ~0fT*lp
TBDATA,1,3.2E5,18.2 !给数据表11输入数据 JBE'B Q@
!取强度折减系数F=3.0时边坡围岩材料属性 j:|60hDz^
MP,EX,12,3E10 ！定义弹性模量 /,5`#Gte_
MP,PRXY,12,0.25 mf@YmKbp
MP,DENS,12,2500 ！定义密度 >w 9)c|
TB,DP,12 !定义数据表12 -3Vx jycY
TBDATA,1,3.0E5,17.1 !给数据表12输入数据 NV;5T3
!建立几何模型 | qHWM
!创建边坡线模型线 y wk;
K,1,,,, ！ 创建关键点1 U:MkA(S%c
K,2,-800,,, Qd!;CoOmZs
K,3,-800,-800,, <_ */
K,4,-800,-1200,, rK~362|mo
K,5,1200,-1200,, _\"P<+!
K,6,1200,-800,, K 3&MR=#^
K,7,1200,,, hzPx8sO
K,8,1200,400,, 1l\O9D +$
K,9,400,400,, 5vY h~|
l,1,2 !连接关键点1、2生成直线L1 nl5K1!1
l,2,3 scqG$~O)
l,3,4 yQhrPw> m
l,4,5 1q~U3'l:$
l,5,6 a-Cp"pKlVY
l,6,7 !j4C: L3F
l,7,8 5"y
p|Yl
l,8,9 P1
(8foZA
l,9,1 svyC(m)'
l,1,7 +M@G 8l
l,3,6 5S$HDO&
!创建边坡面模型 m[oe$yH
al,3,4,5,11 !由线L3、L4、L5、L11生成一个面A1 t2OXm
al,1,2,11,6,10 !由线L1、L2、L1、L6、L10生成一个面A2 _89
_*t(
al,7,8,9,10 !由线L7、L8、L9、L10生成一个面A3 N>Y`>5
/PNUM, area, on !打开面号开关
$7)O&T*q'
!划分网格生成有限元模型 Dt1{]~30
! 划分边坡围岩2网格 /Fh"Gl^
!设置网格份数 #X"\:yN
lsel,s,,,3,5,2 !选择线L3和L5 qPE(Lt1
lesize,all,,,5 !把所选择线分为5段 dWD,iO_"@
lsel,s,,,4,11,7 !选择线L4和L11 =Gd[Qn83.%
lesize,all,,,26 !把所选择线分为26段 D5T\X-+]O
mat,2 !给边坡围岩2赋予2号材料特性 2O9dU 5b
type,1 !采用单元类型1 ~<v.WP<:
mshkey,0 !设定自由网格划分 FTCp3 g
mshape,0 !设定四边形网格划分 wXZ . D}d
amesh,1 ！划分面积A1 -ihF)^"a
! 划分边坡围岩1网格 ' A !Dg
!设置网格份数 p>K'6lCa
lsel,s,,,2,6,4 !选择线L2和L6 uA!T@>vl
lsel,a,,,1 !附加选择线L1 :M|c,SQK
lesize,all,,,10 !把所选择线分为10段 o^&;
`XOd
lsel,s,,,8,10,2 !选择线L8和L10 30>TxL=&
lesize,all,,,16 !把所选择线分为16段 P>pkLP}
Vo
lsel,s,,,7,9,2 K;6K!6J:[
lesize,all,,,12 R_ vZh|
mat,1 !给边坡围岩1赋予1号材料特性 l$,l3
type,1 ) 0AE*S
mshkey,0 !设定自由网格划分 2t [c^J
mshape,0 !设定四边形网格划分 j#&sZ$HQ4
amesh,2 !划分面积A2 Z6>:k,-Ot
amesh,3 ！划分面积A3 4>Uo0NfL
allsel )\^o<x2S
!保存网格模型 7L;yN..0
save,Slope-grid.db :`3b|u=KZ
(3)施加约束和荷载 ~uC4>+dk
!给边坡模型施加约束 }jiqUBn%
!边坡两侧施加X方向约束 /l+x&xYD
nsel,s,loc,x,-800 !选择X=-800线上所有节点 ADv
a@P
nsel,a,loc,x,1200 !选择X=1200线上所有节点 l,*Q?q
d,all,ux !对所选择节点约束X方向位移 Q(k $HP
allsel
UUb!2sO
!边坡底部施加约束 Yc?t aL)
nsel,s,loc,y,-1200 !选择Y=--1200线上所有节点 bq[Q
d,all,uy !对所选择节点约束X、Y方向位移 2y_rsu\
d,all,ux OM!ES%c,
!施加重力加速度 J~gf Mp.
acel,,9.8 Kz3u
(4)求解 f ` A
/solu &O0+\A9tP
!求解设置 ~8 pf.^,fi
antype,static !设定为静力求解 5B2x#
m|8
nsubst,100 !设定最大子步数为100 a4 `@z:l
pred,on !打开时间步长预测器 bHS2;K~
nropt,full !设定牛顿-拉普森选项 @aCg1Rm
nlgeom,on !打开大位移效果 8ZG'?A+{
lnsrch,on !打开线性搜索 m1F<L
outres,all,all !输出所有项 #4na>G|
cnvtol,f,,0.005,2,0.5 !力收敛准则设定 y13=y}dyDH
cnvtol,u,,0.05,2,1 !位移收敛准则设定 gyCb\y+\a
!边坡在强度折减系数F=1时求解 O|y-nAZgU
allsel $o]zNW;X
solve !进行求解 tx5_e [
save,F1,db !把F=1时求解结果保存 ;S`N q%,
!边坡在强度折减系数F=1.2时求解 308w 0eP
finish .j}u'!LKul
/solu ?]9uHrdsN}
allsel Rdt8jY6F/
asel,s,area,,2,3,1 !选择面积A2、A3 B"KsYB79t
MPCHG,3,all ! 把所选择单元材料号改为3 ;%dkwKO
allsel *$# r %
solve !进行求解 i'e^[oZ
save,F1.2,db !把F=1.2时求解结果保存 9d[0i#` :q
!边坡在强度折减系数F=1.4时求解 xA!o"VZPq7
finish 6^{ hY^Z
/solu $Q{1^
allsel lBG* P>;
asel,s,area,,2,3,1 !选择面积A2、A3 D<++6HN&#
MPCHG,4,all !把所选择单元材料号改为4 ?783LBe
allsel f!LZT! y
solve !进行求解 ni y@'
save,F1.4,db !把F=1.4时求解结果保存 crgYr$@s?
!边坡在强度折减系数F=1.6时求解 4#2iL+
finish #(aROTV5a
/solu ~BS*x+M
allsel p6Z]oL q
asel,s,area,,2,3,1 !选择面积A2、A3 ~iwEhF
MPCHG,5,all !把所选择单元材料号改为5 i $I|JJ J
allsel ~d5"<`<^o
solve ^Y'J0v2
save,F1.6,db !把F=1.6时求解结果保存 _\]D<\St
!边坡在强度折减系数F=1.8时求解 R X2=
iO"
finish o4~ft!>
/solu 1v2wP2]|;
allsel 3sp*.dk
asel,s,area,,2,3,1 !选择面积A2、A3 sgX}`JH?z
MPCHG,6,all !把所选择单元材料号改为6 {f^30F w
allsel w,}}mC)\*
solve m\Tq0cT$
save,F1.8,db !把F=1.8时求解结果保存 n"FOCcTIs
!边坡在强度折减系数F=2.0时求解 $d8A_CUU
finish [7I| 8
/solu &2xYG{Z
allsel ljt1:@SN(
asel,s,area,,2,3,1 !选择面积A2、A3 RTY$oUqlZ
MPCHG,7,all !把所选择单元材料号改为7 3:Z(tM&-O
allsel o=`9JKB~
solve 4` 8.\
save,F2.0,db !把F=2.0时求解结果保存 (
?/0$DB
!边坡在强度折减系数F=2.2时求解 _a<PUdP
finish ['R=@.
/solu
Gzp)OHgJ
allsel hLm9"N 'Pf
asel,s,area,,2,3,1 !选择面积A2、A3 :!} zdeRJ
MPCHG,8,all !把所选择单元材料号改为8 B. P64"w
allsel lC _zSmT
solve hq,;H40%/
save,F2.2,db !把F=2.2时求解结果保存 Cg{$$&_(Hj
!边坡在强度折减系数F=2.4时求解 [tD*\\IA
finish W^=89I4]
/solu iBo-ANnK9
allsel $\^]MxI
asel,s,area,,2,3,1 !选择面积A2、A3 Uw&+zJ
MPCHG,9,all !把所选择单元材料号改为9 V'mpl
allsel 4uftx1o
solve 2{ V|
save,F2.4,db !把F=2.4时求解结果保存 t&P5Zw*B
!边坡在强度折减系数F=2.6时求解 GU[ Cq=k
finish Tl%4L %
bE
/solu rM=Hd/ki5
allsel LWQ BGiJj
asel,s,area,,2,3,1 !选择面积A2、A3 {eZ j[* P
MPCHG,10,all !把所选择单元材料号改为10 f "&q~V4?
allsel #[KwR\b{:+
solve rXm!3E6JL
save,F2.6,db !把F=2.6时求解结果保存 :X4\4B*~
!边坡在强度折减系数F=2.8时求解 A\# ? rK
finish \~?s= LT
/solu <BU|?T 6~
allsel E?9_i
:IX
asel,s,area,,2,3,1 !选择面积A2、A3 (B$FX<K3
MPCHG,11,all !把所选择单元材料号改为11 1MahFeQ[
allsel /pk; E$qv
solve 4 l $8 lYi
save,F2.8,db !把F=2.8时求解结果保存 zf o.S[R@
!边坡在强度折减系数F=3.0时求解 ycE<7W
finish _-!6@^+
/solu "5 y^s !/
allsel CIaabn
asel,s,area,,2,3,1 !选择面积A2、A3 FBY~Z$o0.
MPCHG,12,all !把所选择单元材料号改为12 6wu/6DO
allsel l&|{uk
solve w5*18L=O\
save,F3.0,db !把F=3.0时求解结果保存 !k s<VJh
(5）后处理 = oh6;Ojt
/post1 ！进入后处理 c=b\9!hr_E
!边坡在强度折减系数F=1时结果分析 XdS< 51 C
Resume,'F1','db' !读入边坡在强度折减系数F=1时 ^_=0.:QaW
set,1,last !读入后一个子步 /GeS(xzQ
pldisp,1 !绘制边坡模型变形图 GUp51*#XR
plnsol,u,x !绘制边坡模型水平方向位移云图 ZDDwh&h
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 H~G=0_S
!边坡在强度折减系数F=1.2时结果分析 ,@!d%rL:4]
Resume,'F1.2','db' !读入边坡在强度折减系数F=1.2时 CqX%V":2
set,1,last !读入后一个子步 F`- [h )e.
pldisp,1 !绘制边坡模型变形图 aZ0H)
plnsol,u,x !绘制边坡模型水平方向位移云图 kcOp O<oE
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 ix#
!边坡在强度折减系数F=1.4时结果分析 % qE#^ U
Resume,'F1.4','db' !读入边坡在强度折减系数F=1.4时 D$mrnm4d
set,1,last !读入后一个子步 ?x[ >g!r
pldisp,1 !绘制边坡模型变形图 " LxJPt\
plnsol,u,x !绘制边坡模型水平方向位移云图 kW:!$MX!
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 @2$8o]et
!边坡在强度折减系数F=1.6时结果分析 a<o0B{7{BM
Resume,'F1.6','db' !读入边坡在强度折减系数F=1.6时 }`M6+.z3F
set,1,last !读入后一个子步 y]CJOC)/K
pldisp,1 !绘制边坡模型变形图 FN,uD:a
plnsol,u,x !绘制边坡模型水平方向位移云图 {(q U n
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 zp9 ?Ia
!边坡在强度折减系数F=1.8时结果分析 ` q@ ~78`
Resume,'F1.8','db' !读入边坡在强度折减系数F=1.8时 =f|>7m.p
set,1,last !读入后一个子步 EV(/@kN2
pldisp,1 !绘制边坡模型变形图 hy]AH)?pR
plnsol,u,x !绘制边坡模型水平方向位移云图 |1<Z3\+_/
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 HkV/+ {;S~
!边坡在强度折减系数F=2.0时结果分析 *
;M?R?+
Resume,'F2.0','db' !读入边坡在强度折减系数F=2.0时 CTP%
set,1,last !读入后一个子步 8/F2V?iT
pldisp,1 !绘制边坡模型变形图 #-f^; =7
plnsol,u,x !绘制边坡模型水平方向位移云图 5Y&@
:Y
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 lL DHx3+
!边坡在强度折减系数F=2.2时结果分析 (qG$u&
Resume,'F2.2','db' !读入边坡在强度折减系数F=2.2时 mPD'"
set,1,last !读入后一个子步 Qdn:4yk
pldisp,1 !绘制边坡模型变形图 (i? ^g &
plnsol,u,x !绘制边坡模型水平方向位移云图 j?mJ1J5
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 6h,'#|:d
!边坡在强度折减系数F=2.4时结果分析 b
g D Dys
Resume,'F2.4','db' !读入边坡在强度折减系数F=2.4时 [ho 'Pc3A<
set,1,last !读入后一个子步 K5p h x
pldisp,1 !绘制边坡模型变形图 *7v PU:Q[
plnsol,u,x !绘制边坡模型水平方向位移云图 *AA1e}R{B
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 "Jwz.,Y\
!边坡在强度折减系数F=2.6时结果分析 aLq; a
Resume,'F2.6','db' !读入边坡在强度折减系数F=2.6时 Grub1=6l
set,1,last !读入后一个子步 0=5i\*5 p
pldisp,1 !绘制边坡模型变形图 +]e4c;`ko}
plnsol,u,x !绘制边坡模型水平方向位移云图 qB:`tHy
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 e=R}
4`
!边坡在强度折减系数F=2.8时结果分析 .F6#s
Resume,'F2.8','db' !读入边坡在强度折减系数F=2.8时 mzw*6e2T
set,1,last !读入后一个子步 <5#e.w
pldisp,1 !绘制边坡模型变形图 kL3=7t^ 1
plnsol,u,x !绘制边坡模型水平方向位移云图 v6n(<0:
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 xD#/@E1'Y
!边坡在强度折减系数F=3.0时结果分析 DVS7N_cx2o
Resume,'F3.0','db' !读入边坡在强度折减系数F=3.0时 lz* 2wGI9
set,1,last !读入后一个子步 ri^yal<'
pldisp,1 !绘制边坡模型变形图 jFc{$#g-
plnsol,u,x !绘制边坡模型水平方向位移云图 8xv\Zj +
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 kj#yG"3+
TSk6Q'L\v
~k%\ LZ3s
l
)4 OV>
P OG5 x
Uy_= #&jg
+O H."4Z
