各位朋友： z77>W}d
大家好！我通过ansys用有限元强度折减法分析了一个边坡的稳定性，结果不同的工况下得出的塑性区是一样的，我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激！附命令流和塑性区如下：(1)创建物理环境 cS9 jGD92
/COM, Structural ! 指定结构分析 -"dt3$ju
/TITLE, Slope stability Analysis ! 定义工作标题 R8Vf6]s_
/FILNAM,support,1 ! 定义工作文件名 mA+:)?e5~
(2)建立模型 uoI7'
:Nv
!进入前处理器 E6 -*2U)k+
/PREP7 m)Kg6/MV.
!定义单元类型 -;L'Jb>s76
ET,1,PLANE82 !定义边坡围岩单元 9DNp
KEYOPT,1,3,2 !求解类型为平面应变 \>-
M&C
!定义材料属性 r.\L@Y<
!边坡围岩1材料属性 X$ ZV Y2
MP,EX,1,3E10 X&?s:A
MP,PRXY,1,0.25 ;x/eb g
MP,DENS,1,2500 !定义塑性模型 P a{
TB,DP,1 !定义数据表1 F9u?+y-xb
TBDATA,1,9E5,42.8 !给数据表输入数据 `{G&i\"n
!边坡围岩2材料属性 !采用弹性模型 HT .*r6Y>g
MP,EX,2,3.2E10 ！定义弹性模量 KZw~Ch}b9
MP,PRXY,2,0.24 !定义泊松比 UE.kR+1
MP,DENS,2,2700 !定义密度 PVBz~rG
!取强度折减系数F=1.2时边坡围岩材料属性 nOL"6%q
MP,EX,3,3E10 ~ $g:
MP,PRXY,3,0.25 nJ2x;';lA
MP,DENS,3,2500 i9DD)Y<
TB,DP,3 !定义数据表3 nm,LKS7
TBDATA,1,7.5E5,37.7 !给数据表输入数据 I'BhN#GhX
!取强度折减系数F=1.4时边坡围岩材料属性 SscB&{f
MP,EX,4,3E10 .PUp3 X-
MP,PRXY,4,0.25 !定义泊松比 VIP7j(#t_g
MP,DENS,4,2500 Ey|_e3Lf[
TB,DP,4 !定义数据表4 0\ = du
TBDATA,1,6.4E5,33.5 !给数据表输入数据 T"_'sSI>tF
!取强度折减系数F=1.6时边坡围岩材料属性 *(F`NJ 3
MP,EX,5,3E10 73D< wMgZF
MP,PRXY,5,0.25 @a?7D;+<
MP,DENS,5,2500 ！定义密度 *7 jz(iX
TB,DP,5 !定义数据表5 ? 1S sF>|
TBDATA,1,5.6E5,30 !给数据表输入数据 *{DTxEy
!取强度折减系数F=1.8时边坡围岩材料属性 r[a7">n
MP,EX,6,3E10 ！定义弹性模量 pi?MAE*f
MP,PRXY,6,0.25 Gs,:$Im
MP,DENS,6,2500 Y"GNJtsL "
TB,DP,6 !定义数据表6 ;mG*Rad
TBDATA,1,5.0E5,27.2 c>B1cR
!取强度折减系数F=2.0时边坡围岩材料属性 'j6O2=1
MP,EX,7,3E10 o[Jzx2 A<
MP,PRXY,7,0.25 !定义泊松比 ,|({[ 9jA
MP,DENS,7,2500 lb-1z]YwQ
TB,DP,7 !定义数据表7 q*l4h u%3
TBDATA,1,4.5E5,24.8 AAevN3a#nI
!取强度折减系数F=2.2时边坡围岩材料属性 tNi>TkC}`
MP,EX,8,3E10 h^yqrDyJ
MP,PRXY,8,0.25 E[] 5Od5#
MP,DENS,8,2500 ~P1~:AT
TB,DP,8 !定义数据表8 6:7[>|okQ
TBDATA,1,4.09E5,22.8 FvxM
!取强度折减系数F=2.4时边坡围岩材料属性 bP&Q Fc
MP,EX,9,3E10 ?dgyi4J?=`
MP,PRXY,9,0.25 s{Z)<n03
MP,DENS,9,2500 ) CP
TB,DP,9 !定义数据表9 (j&:
TBDATA,1,3.6E5,21.1 !给数据表输入数据 ,yYcjs!=o
!取强度折减系数F=2.6时边坡围岩材料属性 hw^&{x
MP,EX,10,3E10 ！定义弹性模量 52e>f5m.
MP,PRXY,10,0.25 !定义泊松比 MEiP&=gX!
MP,DENS,10,2500 ！定义密度 50^ux:Uv+N
TB,DP,10 !定义数据表10 x;n3 Zr;(
TBDATA,1,3.46E5,19.6 !给数据表10输入数据 m H'~pR>t
!取强度折减系数F=2.8时边坡围岩材料属性 "?I#!t%'
MP,EX,11,3E10 ！定义弹性模量 KUA zJ[>
MP,PRXY,11,0.25 !定义泊松比 U#` e~d t<
MP,DENS,11,2500 bO=|utpk
TB,DP,11 !定义数据表11 ! *N9PUM
TBDATA,1,3.2E5,18.2 !给数据表11输入数据 f2i:I1 p("
!取强度折减系数F=3.0时边坡围岩材料属性 to'7o8Z
MP,EX,12,3E10 ！定义弹性模量 %c }V/v_h
MP,PRXY,12,0.25 <x OX+D
MP,DENS,12,2500 ！定义密度 k {*QU(
TB,DP,12 !定义数据表12 \ H fAKBT
TBDATA,1,3.0E5,17.1 !给数据表12输入数据 dZU #lg
!建立几何模型 J7`mEL>?
!创建边坡线模型线 K3yQ0k
|
K,1,,,, ！ 创建关键点1 .A2$C|a*
K,2,-800,,,
;A*SuFbV
K,3,-800,-800,, H@Z_P p?
K,4,-800,-1200,, zR" c j
K,5,1200,-1200,, 9zu;OK%
K,6,1200,-800,, C$OVN$lL`8
K,7,1200,,, ?Y(
K,8,1200,400,, HL38iXQ(
3
K,9,400,400,, 9P?0D
l,1,2 !连接关键点1、2生成直线L1 m\L`$=eO8
l,2,3 r
)F;8(
l,3,4 EhcJE;S)
l,4,5 52w@.]
l,5,6 `^bP9X_a
l,6,7 bn(N8MFCV
l,7,8 ~]?:v,UIm(
l,8,9 95l)s],
l,9,1 53HA6:Q[
l,1,7 \Z[1m[{
l,3,6 pHV^K v#
!创建边坡面模型 u "[f\l
al,3,4,5,11 !由线L3、L4、L5、L11生成一个面A1 j9p6 rD
al,1,2,11,6,10 !由线L1、L2、L1、L6、L10生成一个面A2 /Z:j:l
al,7,8,9,10 !由线L7、L8、L9、L10生成一个面A3 UVo`jb|>
o
/PNUM, area, on !打开面号开关 Nd~B$venh
!划分网格生成有限元模型 v@ _1V
! 划分边坡围岩2网格 u_(~zs.N]
!设置网格份数 +V3mF_s|z
lsel,s,,,3,5,2 !选择线L3和L5 ^4dE8Ve"@
lesize,all,,,5 !把所选择线分为5段 AUnfhk@$
lsel,s,,,4,11,7 !选择线L4和L11 ".?4 `@7F\
lesize,all,,,26 !把所选择线分为26段 e97G]XLR
mat,2 !给边坡围岩2赋予2号材料特性 3U[O :
type,1 !采用单元类型1 _f1o!4ocx
mshkey,0 !设定自由网格划分 -@pjEI
mshape,0 !设定四边形网格划分 9`tK 9
amesh,1 ！划分面积A1 Pe`(9&iT.
! 划分边坡围岩1网格 &."$kfA+
!设置网格份数 `Ij @;=(
lsel,s,,,2,6,4 !选择线L2和L6 I Dohv[#
lsel,a,,,1 !附加选择线L1 17LhgZs&
lesize,all,,,10 !把所选择线分为10段 "tJ+v*E
lsel,s,,,8,10,2 !选择线L8和L10 ?Nos;_/
lesize,all,,,16 !把所选择线分为16段 *_(X$qfoW
lsel,s,,,7,9,2 &@0~]\,D7
lesize,all,,,12 .ZTvOm'mB^
mat,1 !给边坡围岩1赋予1号材料特性 w'S,{GW
type,1 #[+# bw_6
mshkey,0 !设定自由网格划分 3:)z+#Uk6
mshape,0 !设定四边形网格划分 <DF3! r
amesh,2 !划分面积A2 &d_^k.%y
amesh,3 ！划分面积A3 Ds9)e&yYrb
allsel K@JZ$
!保存网格模型 L5&M@YTH
save,Slope-grid.db kwI``7g8*e
(3)施加约束和荷载 0U '"@A
\
!给边坡模型施加约束 3SVGx< ,2
!边坡两侧施加X方向约束 l/.{F ;3F
nsel,s,loc,x,-800 !选择X=-800线上所有节点 OUX7
*_
nsel,a,loc,x,1200 !选择X=1200线上所有节点 5^B79A"}
d,all,ux !对所选择节点约束X方向位移 ]3KeAJ
allsel C/kf?:j
!边坡底部施加约束 A;oHji#*
nsel,s,loc,y,-1200 !选择Y=--1200线上所有节点 UcRP/LR%C
d,all,uy !对所选择节点约束X、Y方向位移 kK6O ZhLH
d,all,ux %w`d
!施加重力加速度 Hfcpq a
acel,,9.8 yW_yHSx;
(4)求解 +7_qg
i7:
/solu 0=Z[6Q@ :
!求解设置 c+Q'4E0 |
antype,static !设定为静力求解 Qb536RpcTY
nsubst,100 !设定最大子步数为100 r&gvP|W%
pred,on !打开时间步长预测器 ~Amq1KU*Z
nropt,full !设定牛顿-拉普森选项 Q: }]-lJg
nlgeom,on !打开大位移效果 jK|n^5\
lnsrch,on !打开线性搜索 JJ: ku&Mb
outres,all,all !输出所有项 >)>f~ >
cnvtol,f,,0.005,2,0.5 !力收敛准则设定 YpZB-9Krf
cnvtol,u,,0.05,2,1 !位移收敛准则设定 b^o4Q[
!边坡在强度折减系数F=1时求解 djd/QAfSC
allsel {}O~tf_
solve !进行求解 3I6ocj [,
save,F1,db !把F=1时求解结果保存 nm^HL|
!边坡在强度折减系数F=1.2时求解 ;%t Fi
finish JcAsrtrG]
/solu G m~([Ln{
allsel T{{:p\<]_
asel,s,area,,2,3,1 !选择面积A2、A3 77> oQ~q
MPCHG,3,all ! 把所选择单元材料号改为3 7J'%; sH
allsel zl0{lV
solve !进行求解 (3Db}Hnn
save,F1.2,db !把F=1.2时求解结果保存 $8tk|uh
!边坡在强度折减系数F=1.4时求解 <c+.%ka
finish RSX27fb4
/solu {6i|"5_j
allsel c =m'I>A
asel,s,area,,2,3,1 !选择面积A2、A3 1I Xtu
MPCHG,4,all !把所选择单元材料号改为4 zakhJ
allsel X9 x`i
solve !进行求解 O]N /(pe:d
save,F1.4,db !把F=1.4时求解结果保存 57\ 0MQO
!边坡在强度折减系数F=1.6时求解 at uqo3
finish w~lH2U'k}
/solu U4D7@KY +m
allsel Kz HYh
asel,s,area,,2,3,1 !选择面积A2、A3 *z0 !=>(
MPCHG,5,all !把所选择单元材料号改为5 CPz<iU
allsel 9"~ FKM N
solve \!u<)kkyT
save,F1.6,db !把F=1.6时求解结果保存 GppCrQ%Ra|
!边坡在强度折减系数F=1.8时求解 ; !n>
finish mLCD N1UO{
/solu ;ALWL~Xm
allsel I8xdE(o8+
asel,s,area,,2,3,1 !选择面积A2、A3 I)rnF
MPCHG,6,all !把所选择单元材料号改为6 F*72g)hVh
allsel HuhQ|~C+~
solve pMfb(D"
save,F1.8,db !把F=1.8时求解结果保存 (W 1 $+X
!边坡在强度折减系数F=2.0时求解 hD6BP
finish SNj-h>&Mha
/solu *`>BOl+ro
allsel [F-R*}&x
asel,s,area,,2,3,1 !选择面积A2、A3 qO6M5g:
MPCHG,7,all !把所选择单元材料号改为7 (7 I|lf
e
allsel d$!Q6ux;
solve =uP?
?E
save,F2.0,db !把F=2.0时求解结果保存 f$I =o N
!边坡在强度折减系数F=2.2时求解 b:==:d:0s
finish BP/nK.
/solu lM- 9 J?j
allsel !%>RHh[
asel,s,area,,2,3,1 !选择面积A2、A3 Yy *=@qu>g
MPCHG,8,all !把所选择单元材料号改为8 gInh+XZs
allsel =2\2S p
solve h}O tz "
save,F2.2,db !把F=2.2时求解结果保存 80l ei
!边坡在强度折减系数F=2.4时求解 Z?)g' n
finish :i.@d?
/solu ,i:?c
allsel O}M-6!%<,
asel,s,area,,2,3,1 !选择面积A2、A3 ee {ToK
MPCHG,9,all !把所选择单元材料号改为9 j=r1JV
@
allsel {npm9w<;
solve l=DF)#>w
save,F2.4,db !把F=2.4时求解结果保存 2<D | {
!边坡在强度折减系数F=2.6时求解 ! M^O\C)
finish #&Hi0..y
/solu *HoRYCL
allsel 4]o+)d.`(
asel,s,area,,2,3,1 !选择面积A2、A3 ^exU]5nvz
MPCHG,10,all !把所选择单元材料号改为10 (&}[2pb!
allsel x|O7}oj
solve "i0>>@NR'
save,F2.6,db !把F=2.6时求解结果保存
"*d6E}wG
!边坡在强度折减系数F=2.8时求解 ale'-V)5
finish wQ 33Gc
/solu wGyVmC
allsel 25>R^2,LiE
asel,s,area,,2,3,1 !选择面积A2、A3 1!z{{H;W
MPCHG,11,all !把所选择单元材料号改为11 ]az(w&vqg2
allsel e= i X]%^
solve 1 Pk+zBJ$
save,F2.8,db !把F=2.8时求解结果保存 O= S[ n
!边坡在强度折减系数F=3.0时求解 ,L ig6Z`
finish \.L j A_
/solu u)oAQ<w
allsel $d<NN2
asel,s,area,,2,3,1 !选择面积A2、A3 K43 %9=sM
MPCHG,12,all !把所选择单元材料号改为12 M7=,J;@
allsel yqN`R\d
solve (*S<2HN5
save,F3.0,db !把F=3.0时求解结果保存 ng}C$d . I
(5）后处理 ~?\U];l
/post1 ！进入后处理 PUR,r %K`
!边坡在强度折减系数F=1时结果分析 `~XksyT
Resume,'F1','db' !读入边坡在强度折减系数F=1时 ,,Qg"C
set,1,last !读入后一个子步 M9iX_4
pldisp,1 !绘制边坡模型变形图 Xm#W}Y'
plnsol,u,x !绘制边坡模型水平方向位移云图 l 7-lXl"%q
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 q` xc h [H
!边坡在强度折减系数F=1.2时结果分析 &5W;E+Pub
Resume,'F1.2','db' !读入边坡在强度折减系数F=1.2时 A 8-[EBkK
set,1,last !读入后一个子步 wg FX')l:
pldisp,1 !绘制边坡模型变形图 D ,k x B~
plnsol,u,x !绘制边坡模型水平方向位移云图 r8 9o
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 <K CI@
!边坡在强度折减系数F=1.4时结果分析 T5:Q_o]
Resume,'F1.4','db' !读入边坡在强度折减系数F=1.4时 eoiz]L
set,1,last !读入后一个子步 *w0! C:mL&
pldisp,1 !绘制边坡模型变形图 >7W)iwF
plnsol,u,x !绘制边坡模型水平方向位移云图
p9"dm{
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 Yq ]sPE92
!边坡在强度折减系数F=1.6时结果分析 7_\G|Zd
Resume,'F1.6','db' !读入边坡在强度折减系数F=1.6时 TiD#t+g
set,1,last !读入后一个子步 ?I +{S
pldisp,1 !绘制边坡模型变形图 / a\i
plnsol,u,x !绘制边坡模型水平方向位移云图 9]eG |LFD
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 ZiPeP
!边坡在强度折减系数F=1.8时结果分析 2 7 )If E
Resume,'F1.8','db' !读入边坡在强度折减系数F=1.8时 VeQGdyhY
set,1,last !读入后一个子步 :E9 pdx+
pldisp,1 !绘制边坡模型变形图 'I, a 29
plnsol,u,x !绘制边坡模型水平方向位移云图 C^ZoYf8+"m
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 ZID- ~
6
!边坡在强度折减系数F=2.0时结果分析 71{Q#%5U~
Resume,'F2.0','db' !读入边坡在强度折减系数F=2.0时 hO"!q;<eS
set,1,last !读入后一个子步 i^DMnvV.
pldisp,1 !绘制边坡模型变形图 2M x\D
plnsol,u,x !绘制边坡模型水平方向位移云图 uWJJ\
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 2;gv o*k
!边坡在强度折减系数F=2.2时结果分析 v]EMJm6d|
Resume,'F2.2','db' !读入边坡在强度折减系数F=2.2时 _Q+c'q Zkl
set,1,last !读入后一个子步 OlT8pG5Oa
pldisp,1 !绘制边坡模型变形图 CcQc!`YC
plnsol,u,x !绘制边坡模型水平方向位移云图 t4oD> =,92
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 }u;K<<h:
!边坡在强度折减系数F=2.4时结果分析 w"{DLN[Qw
Resume,'F2.4','db' !读入边坡在强度折减系数F=2.4时 0/v]YK.
set,1,last !读入后一个子步 g+B7~Z5,
pldisp,1 !绘制边坡模型变形图 7qqzL_d>
plnsol,u,x !绘制边坡模型水平方向位移云图 "qgwuWb M
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 p8'$@:M\
!边坡在强度折减系数F=2.6时结果分析 p_I^7 $
Resume,'F2.6','db' !读入边坡在强度折减系数F=2.6时 [q|W*[B:@
set,1,last !读入后一个子步 -))S
pldisp,1 !绘制边坡模型变形图 s4fO4.bn m
plnsol,u,x !绘制边坡模型水平方向位移云图 3)Wf BvG
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 Y"lE MY
!边坡在强度折减系数F=2.8时结果分析 ao)';[%9s
Resume,'F2.8','db' !读入边坡在强度折减系数F=2.8时 xX-r<:'tmi
set,1,last !读入后一个子步 tK;xW
pldisp,1 !绘制边坡模型变形图 r R6}
plnsol,u,x !绘制边坡模型水平方向位移云图 wN.S]
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 ) ,yar9C
!边坡在强度折减系数F=3.0时结果分析 e|oMbTZ5m
Resume,'F3.0','db' !读入边坡在强度折减系数F=3.0时 0`"oR3JY
set,1,last !读入后一个子步 fg LY{
pldisp,1 !绘制边坡模型变形图 W5Jw^,iPd
plnsol,u,x !绘制边坡模型水平方向位移云图 @ujwN([I
plnsol,eppl,eqv !绘制边坡模型塑性应变云图
