[求助]一个有限元折减法分析边坡的实例 [复制链接]

各位朋友： CK|AXz+EN  
    大家好！我通过ansys用有限元强度折减法分析了一个边坡的稳定性，结果不同的工况下得出的塑性区是一样的，我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激！附命令流和塑性区如下：(1)创建物理环境 m J$[X  
/COM, Structural                              ! 指定结构分析 #cW:04  
/TITLE, Slope stability Analysis                  ! 定义工作标题 {%P2.:  
/FILNAM,support,1                            ! 定义工作文件名 xX{Zh;M&[  
(2)建立模型 9AQ,@xP|  
!进入前处理器 ]mNsG0r6  
/PREP7 `m#G'E I  
!定义单元类型 e0ni  
ET,1,PLANE82                              !定义边坡围岩单元 L})*ck
 
KEYOPT,1,3,2                                !求解类型为平面应变 zLg$|@E&  
!定义材料属性 $
T66%wX  
!边坡围岩1材料属性  5.oY$tb(  
MP,EX,1,3E10  o /1+ }f  
MP,PRXY,1,0.25    :J x%K  
MP,DENS,1,2500                            !定义塑性模型 =WZ9|
e  
TB,DP,1                                    !定义数据表1 8cMX
=P  
TBDATA,1,9E5,42.8                          !给数据表输入数据 aMkuyqPf{  
!边坡围岩2材料属性                        !采用弹性模型 `)KGajB  
MP,EX,2,3.2E10                              ！定义弹性模量 ySDo(EI4  
MP,PRXY,2,0.24                              !定义泊松比 MF*4E9Ue.  
MP,DENS,2,2700                            !定义密度 N'l
2$8  
!取强度折减系数F=1.2时边坡围岩材料属性 ,z`D}<3  
MP,EX,3,3E10  (]&B'1b  
MP,PRXY,3,0.25    <}c7E3Uc  
MP,DENS,3,2500 "cjD-42  
TB,DP,3                                    !定义数据表3 %BMlcm7Ec  
TBDATA,1,7.5E5,37.7                        !给数据表输入数据 "
;T a8  
!取强度折减系数F=1.4时边坡围岩材料属性 :f_oN3F p  
MP,EX,4,3E10  HFFrS%  
MP,PRXY,4,0.25                              !定义泊松比 #uC}IX2n  
MP,DENS,4,2500 luac  
TB,DP,4                                    !定义数据表4 FzCXA=m  
TBDATA,1,6.4E5,33.5                        !给数据表输入数据 |f1^&97=+  
!取强度折减系数F=1.6时边坡围岩材料属性 }y J,&N'p  
MP,EX,5,3E10  2>9..c  
MP,PRXY,5,0.25    p0l.f
`B  
MP,DENS,5,2500                            ！定义密度 FjiIB1 T  
TB,DP,5                                    !定义数据表5 VQ2'a/s  
TBDATA,1,5.6E5,30                          !给数据表输入数据 -ULgVGYKK  
!取强度折减系数F=1.8时边坡围岩材料属性 GiK,+M"d  
MP,EX,6,3E10                              ！定义弹性模量 ![vy{U.:`  
MP,PRXY,6,0.25    0
8k  
MP,DENS,6,2500 g3Hi5[-H  
TB,DP,6                                    !定义数据表6 Qgf|obrEi6  
TBDATA,1,5.0E5,27.2 W >}T$a}\  
!取强度折减系数F=2.0时边坡围岩材料属性 &m9= q|;m  
MP,EX,7,3E10  %1{O  
MP,PRXY,7,0.25                            !定义泊松比 BXxJra/V  
MP,DENS,7,2500 ''!j:49  
TB,DP,7                                    !定义数据表7 +7 j/.R  
TBDATA,1,4.5E5,24.8 q@VIFmqY!  
!取强度折减系数F=2.2时边坡围岩材料属性 7(C)vtEO:  
MP,EX,8,3E10  nox-)e  
MP,PRXY,8,0.25    KjF8T7%  
MP,DENS,8,2500 saQo]6#  
TB,DP,8                                    !定义数据表8 %gSmOW2.c^  
TBDATA,1,4.09E5,22.8 &t_TLV 8T  
!取强度折减系数F=2.4时边坡围岩材料属性
rf'A+q  
MP,EX,9,3E10  e}7!A  
MP,PRXY,9,0.25    Vu4LC&q  
MP,DENS,9,2500 mF4OLG3L0  
TB,DP,9                                      !定义数据表9 ePaC8sd0
  
TBDATA,1,3.6E5,21.1                          !给数据表输入数据 )$a6l8  
!取强度折减系数F=2.6时边坡围岩材料属性 eOXu^M>:F  
MP,EX,10,3E10                                ！定义弹性模量  EKN<KnU%  
MP,PRXY,10,0.25                              !定义泊松比 :=!6w  
MP,DENS,10,2500                          ！定义密度 55] MRv  
TB,DP,10                                  !定义数据表10 q;f L@L@-  
TBDATA,1,3.46E5,19.6                        !给数据表10输入数据 u WdKG({][  
!取强度折减系数F=2.8时边坡围岩材料属性 'gD./|Z0  
MP,EX,11,3E10                              ！定义弹性模量 cG@Wo8+  
MP,PRXY,11,0.25                            !定义泊松比 H.]<fvP  
MP,DENS,11,2500 Nr)DU.f  
TB,DP,11                                  !定义数据表11 ux6)K= ]  
TBDATA,1,3.2E5,18.2                        !给数据表11输入数据 3i4m!g5Z?  
!取强度折减系数F=3.0时边坡围岩材料属性 MU `!sb*  
MP,EX,12,3E10                              ！定义弹性模量 >f-RzQ k  
MP,PRXY,12,0.25    u$ci{<  
MP,DENS,12,2500                            ！定义密度 ER[$TH&  
TB,DP,12                                  !定义数据表12 'IVC!uL,%  
TBDATA,1,3.0E5,17.1                        !给数据表12输入数据 z^4+Un  
!建立几何模型 0@EI@X;q  
!创建边坡线模型线 60e{]}Z  
K,1,,,,                                      ！ 创建关键点1 1~2+w]-kU  
K,2,-800,,, 'Ce?!UO  
K,3,-800,-800,,  P%vouC0W  
K,4,-800,-1200,, K$:btWSm  
K,5,1200,-1200,,  T?jN/}qg  
K,6,1200,-800,,  >){}nlQf  
K,7,1200,,, `qCL&(`%  
K,8,1200,400,, $/(``8li_  
K,9,400,400,, RX
^8`}N  
l,1,2                                        !连接关键点1、2生成直线L1 9asA-'fZ  
l,2,3 I1dOMu9  
l,3,4 nG?Z* n  
l,4,5 @<4U &  
l,5,6 s,/C^E  
l,6,7 l>BM}hS  
l,7,8 ;<+Z}d/g9  
l,8,9  =>Qd  
l,9,1 4R8Qn^  
l,1,7 i=rA;2>  
l,3,6 rTJqw@]#WH  
!创建边坡面模型 @(N} {om  
al,3,4,5,11                              !由线L3、L4、L5、L11生成一个面A1 'iwT
vkf{  
al,1,2,11,6,10                            !由线L1、L2、L1、L6、L10生成一个面A2 At[Sk
G}b  
al,7,8,9,10                              !由线L7、L8、L9、L10生成一个面A3 Z?9G2<i  
/PNUM, area, on                        !打开面号开关 9oP  
!划分网格生成有限元模型 ma
QxU(  
! 划分边坡围岩2网格 a%6=sqxE  
!设置网格份数 e8xNZG;  
lsel,s,,,3,5,2                              !选择线L3和L5 *J':U>p  
lesize,all,,,5                              !把所选择线分为5段 jJ2{g> P0P  
lsel,s,,,4,11,7                              !选择线L4和L11 gA1j'!\6l9  
lesize,all,,,26                              !把所选择线分为26段 cf"!U+x
  
mat,2                                    !给边坡围岩2赋予2号材料特性 \S?-[v*{  
type,1                                  !采用单元类型1 ,Tx38  
mshkey,0                                !设定自由网格划分 j!;E>`g  
mshape,0                                !设定四边形网格划分 )u]=^  
amesh,1                                  ！划分面积A1 0A8G8^T  
! 划分边坡围岩1网格 ]+w 27!  
!设置网格份数 $DnJ/hg;qD  
lsel,s,,,2,6,4                              !选择线L2和L6 $`vXI%|.  
lsel,a,,,1                                  !附加选择线L1 e85E+S%  
lesize,all,,,10                              !把所选择线分为10段 m@L>6;*  
lsel,s,,,8,10,2                              !选择线L8和L10 gOE?  
lesize,all,,,16                              !把所选择线分为16段 }"%mP 4]&  
lsel,s,,,7,9,2 i/!KUbt  
lesize,all,,,12 < %<nh`D  
mat,1                                    !给边坡围岩1赋予1号材料特性 WHLTJ]OB  
type,1 <1x u&Z7  
mshkey,0                                  !设定自由网格划分 d#ab"&$bv  
mshape,0                                  !设定四边形网格划分 :8N by$#V  
amesh,2                                  !划分面积A2 #KHj.Vg  
amesh,3                                  ！划分面积A3 w6lx
&K-  
allsel B !
rb*"[  
!保存网格模型 ^Mhh2v  
save,Slope-grid.db VtU2&  
(3)施加约束和荷载 vJ 28A  
!给边坡模型施加约束 M-+!z5q~d  
!边坡两侧施加X方向约束 8)Vl2z  
nsel,s,loc,x,-800                            !选择X=-800线上所有节点 *qm>py`O  
nsel,a,loc,x,1200                            !选择X=1200线上所有节点 qAlX#]  
d,all,ux                                    !对所选择节点约束X方向位移 =0;njL(7;  
allsel 3Y +;8ld  
!边坡底部施加约束 zc,X5R1  
nsel,s,loc,y,-1200                          !选择Y=--1200线上所有节点 JLu$UR4  
d,all,uy                                  !对所选择节点约束X、Y方向位移 <RH%FhT  
d,all,ux                                  !Bg
^-F:N  
!施加重力加速度 ~rXLb:  
acel,,9.8 Y[rCF=ZVH  
(4)求解 0Am\02R.C,  
/solu od,,2pwK+  
!求解设置 B_8JwMJu3  
antype,static                                  !设定为静力求解 ! z5c+JqN  
nsubst,100                                  !设定最大子步数为100 rF{,]U9`  
pred,on                                      !打开时间步长预测器 +~ Y.m8  
nropt,full                                    !设定牛顿-拉普森选项 .]<gm9l  
nlgeom,on                                    !打开大位移效果 #BH]`A J  
lnsrch,on                                    !打开线性搜索 x1Gc|K/-  
outres,all,all                                  !输出所有项 X_rv}  
cnvtol,f,,0.005,2,0.5                            !力收敛准则设定 30sA\TZ  
cnvtol,u,,0.05,2,1                              !位移收敛准则设定 eE\T,u5:  
!边坡在强度折减系数F=1时求解 AxO.adQE%  
allsel {S@, ,  
solve                                        !进行求解 qzZ;{>_f  
save,F1,db                                    !把F=1时求解结果保存 h+YPyeAs  
!边坡在强度折减系数F=1.2时求解 oGbh*  
finish !g|[A7<|  
/solu 8W7ET@`  
allsel :qShP3^  
asel,s,area,,2,3,1                              !选择面积A2、A3 dg+"G|nr  
MPCHG,3,all                                ! 把所选择单元材料号改为3
7[R`52pP  
allsel X%;4G^%ZI  
solve                                        !进行求解 ALInJ{X  
save,F1.2,db                                  !把F=1.2时求解结果保存                                                                                  j~9Y0jz_  
!边坡在强度折减系数F=1.4时求解 vYo~36  
finish }y(cv}8
Y  
/solu ~Xr[d07bC  
allsel w=ufJRj  
asel,s,area,,2,3,1                              !选择面积A2、A3 /%{CJ0Y  
MPCHG,4,all                                !把所选择单元材料号改为4 Ia{t/IX\[  
allsel 0dD.xuor  
solve                                        !进行求解 ?a?4;Y!  
save,F1.4,db                                  !把F=1.4时求解结果保存                            20Jlf?  
!边坡在强度折减系数F=1.6时求解 S~|\bnE  
finish L$,Kdpj  
/solu #W_-S0>&  
allsel cmd7-2  
asel,s,area,,2,3,1                              !选择面积A2、A3 j!0-3YKv  
MPCHG,5,all                                !把所选择单元材料号改为5 "s`#`'  
allsel x%W~@_  
solve *kj+6`:CPs  
save,F1.6,db                                  !把F=1.6时求解结果保存                                                                                  ds{)p<LpT  
!边坡在强度折减系数F=1.8时求解 I.SMn,N  
finish ?01ru5ys/o  
/solu GFnwj<V+{  
allsel 6vU%Y_n=y]  
asel,s,area,,2,3,1                                !选择面积A2、A3 m5P@F@  
MPCHG,6,all                                  !把所选择单元材料号改为6 ;{e'q?Y  
allsel :Z83*SPc  
solve tm_\(  
save,F1.8,db                                  !把F=1.8时求解结果保存  u2I@ fH/  
!边坡在强度折减系数F=2.0时求解 }<'ki ;  
finish Aaz:C5dtU  
/solu R.GDCGAL  
allsel "/taatcH  
asel,s,area,,2,3,1                                !选择面积A2、A3 6aWnj*dF  
MPCHG,7,all                                  !把所选择单元材料号改为7 B~O<?@]d  
allsel `Uvc^  
solve *N6sxFs  
save,F2.0,db                                  !把F=2.0时求解结果保存                            ,Vz-w;oDn  
!边坡在强度折减系数F=2.2时求解 \"5p)(  
finish "N}MhcdS  
/solu =dWqB&  
allsel mWLiXKnb  
asel,s,area,,2,3,1                                !选择面积A2、A3 Vy=+G~  
MPCHG,8,all                                  !把所选择单元材料号改为8 M3JV^{O/DV  
allsel 7MKZ*f@x;  
solve `bLJwJ7  
save,F2.2,db                                  !把F=2.2时求解结果保存                                                                                  k@>(sXs  
!边坡在强度折减系数F=2.4时求解 9"M-nH*<  
finish )hVn/*mH  
/solu -&%! 4(Je  
allsel o?#-Tkb  
asel,s,area,,2,3,1                                !选择面积A2、A3 XN3'k[  
MPCHG,9,all                                  !把所选择单元材料号改为9 n%QWs1 b  
allsel 9%MgAik(  
solve X@A8~kj1  
save,F2.4,db                                  !把F=2.4时求解结果保存  $}0\sj%  
!边坡在强度折减系数F=2.6时求解 0juP"v$C>  
finish ybIqn0&[  
/solu i"4;{C{s  
allsel iUqD>OV  
asel,s,area,,2,3,1                                !选择面积A2、A3 ]\ZmK0q<:  
MPCHG,10,all                                  !把所选择单元材料号改为10 9El{>&Fs4  
allsel =#{q#COK$  
solve ,WoV)L'?  
save,F2.6,db                                  !把F=2.6时求解结果保存                            :#N]s  
!边坡在强度折减系数F=2.8时求解 "b)EH/s  
finish i0?/\@gd  
/solu K
z]\o"K
 
allsel E429<LQI/  
asel,s,area,,2,3,1                                !选择面积A2、A3 d@Wze[M?0  
MPCHG,11,all                                !把所选择单元材料号改为11 3_{rXtT)'  
allsel }p8iq  
solve YCdtf7P=q  
save,F2.8,db                                  !把F=2.8时求解结果保存                                                                                  mK^E@uxN  
!边坡在强度折减系数F=3.0时求解 Y|KT3  
finish 6J-=6t|  
/solu Cw5B p9  
allsel \t=#MzjR  
asel,s,area,,2,3,1                            !选择面积A2、A3 nLrCy5R:  
MPCHG,12,all                              !把所选择单元材料号改为12 .^ba*qb`{  
allsel &$~irI  
solve 85A7YraL  
save,F3.0,db                                !把F=3.0时求解结果保存 yi-0CHo  
(5）后处理 dtV7YPz4+  
/post1                                      ！进入后处理 ClPE_Cfw~  
!边坡在强度折减系数F=1时结果分析 /(%!txSNEt  
Resume,'F1','db'                              !读入边坡在强度折减系数F=1时 T w1&<S  
set,1,last                                    !读入后一个子步 .`&k`  
pldisp,1                                      !绘制边坡模型变形图 wRX#^;O9?>  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 7WNUHLEt  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 f]Rh<N$  
!边坡在强度折减系数F=1.2时结果分析 G) 37?A)  
Resume,'F1.2','db'                            !读入边坡在强度折减系数F=1.2时 _0iV6Bj  
set,1,last                                    !读入后一个子步 rfh`;G5s  
pldisp,1                                      !绘制边坡模型变形图 <e@4;Z(h04  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图  -f<}lhmQ  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 l
pbcpB  
!边坡在强度折减系数F=1.4时结果分析 =C7<I   
Resume,'F1.4','db'                            !读入边坡在强度折减系数F=1.4时 * COC&  
set,1,last                                    !读入后一个子步 "837b/>/  
pldisp,1                                      !绘制边坡模型变形图 .GCJA`0h  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 }+)q/]%  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 nH+wU;M  
!边坡在强度折减系数F=1.6时结果分析 e%=SgXl2t  
Resume,'F1.6','db'                            !读入边坡在强度折减系数F=1.6时 R|*Eg,1g -  
set,1,last                                    !读入后一个子步 |`
AJ
P  
pldisp,1                                      !绘制边坡模型变形图 IfP?+y
Pa  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 w,<n5dMv  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 C+Wa(K  
!边坡在强度折减系数F=1.8时结果分析 7eF
FKl  
Resume,'F1.8','db'                            !读入边坡在强度折减系数F=1.8时 6r h#ATep  
set,1,last                                    !读入后一个子步 6_;n b
qY&  
pldisp,1                                      !绘制边坡模型变形图 idjk uB(6  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图  |vBy=:  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 v++&%  
!边坡在强度折减系数F=2.0时结果分析 ~*tn|?%  
Resume,'F2.0','db'                            !读入边坡在强度折减系数F=2.0时 F$YT4414  
set,1,last                                    !读入后一个子步 |2jA4C2L}  
pldisp,1                                      !绘制边坡模型变形图 #3FsK  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 !_iv~Q zv  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 O6\c1ha  
!边坡在强度折减系数F=2.2时结果分析 sWVapup?  
Resume,'F2.2','db'                            !读入边坡在强度折减系数F=2.2时 @$~;vS  
set,1,last                                    !读入后一个子步 &hM7y
7  
pldisp,1                                      !绘制边坡模型变形图 ~svea>Fmr  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 hij 9r z  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 +Je%8jH  
!边坡在强度折减系数F=2.4时结果分析 Zb<D%9  
Resume,'F2.4','db'                            !读入边坡在强度折减系数F=2.4时 }YU\}T-P  
set,1,last                                    !读入后一个子步 Vw&#
Lo  
pldisp,1                                      !绘制边坡模型变形图 W*u$e8i7  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 mWMtz]M}  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 "_+8z_  
!边坡在强度折减系数F=2.6时结果分析 $?;)uoAg  
Resume,'F2.6','db'                            !读入边坡在强度折减系数F=2.6时 Dh68=F0  
set,1,last                                    !读入后一个子步 @F|pKf:M+  
pldisp,1                                      !绘制边坡模型变形图 I 6L3M\+-  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 -AB0uMot  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 iBY16_q  
!边坡在强度折减系数F=2.8时结果分析 dvAz}3p0]  
Resume,'F2.8','db'                            !读入边坡在强度折减系数F=2.8时
7oCY@>(f  
set,1,last                                    !读入后一个子步 .kfx\,lgm  
pldisp,1                                      !绘制边坡模型变形图 py%:,hi  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 OgzKX>N`A  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 X'/'r
.b6  
!边坡在强度折减系数F=3.0时结果分析 gA]3h8%w  
Resume,'F3.0','db'                            !读入边坡在强度折减系数F=3.0时 `"#hhKG  
set,1,last                                    !读入后一个子步 "z*?#&?,  
pldisp,1                                      !绘制边坡模型变形图 60xL.Z   
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 B~YOU3  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 B@8lD\  
/3;]e3x  
(De>k8  
qtz~Y~h|>  
i G%R'/*  
q0nIJ(  
:=:m4UJ
b