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

各位朋友： Q)6va}2ai  
    大家好！我通过ansys用有限元强度折减法分析了一个边坡的稳定性，结果不同的工况下得出的塑性区是一样的，我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激！附命令流和塑性区如下：(1)创建物理环境 |+EKF.K  
/COM, Structural                              ! 指定结构分析 w6,*9(;$Pk  
/TITLE, Slope stability Analysis                  ! 定义工作标题 LdTIR]  
/FILNAM,support,1                            ! 定义工作文件名 {/UhUG  
(2)建立模型 ,?b78_,2  
!进入前处理器 I"Q<n[g0'  
/PREP7 P1d,8~;  
!定义单元类型 ua& @GXvZ  
ET,1,PLANE82                              !定义边坡围岩单元 03E3cp"  
KEYOPT,1,3,2                                !求解类型为平面应变 ;~tsF.=
 
!定义材料属性 C!UEXj`l9  
!边坡围岩1材料属性  xUj2]Q>R+  
MP,EX,1,3E10  O2#S: ~h  
MP,PRXY,1,0.25    #vnT&FN0[  
MP,DENS,1,2500                            !定义塑性模型 :I
/  
TB,DP,1                                    !定义数据表1 {OxWcK\2@h  
TBDATA,1,9E5,42.8                          !给数据表输入数据 el@XK}<dr  
!边坡围岩2材料属性                        !采用弹性模型 C6k4g75U2  
MP,EX,2,3.2E10                              ！定义弹性模量 (wL3 +  
MP,PRXY,2,0.24                              !定义泊松比 ?n*fy  
MP,DENS,2,2700                            !定义密度 X5E '*W  
!取强度折减系数F=1.2时边坡围岩材料属性 O3_B<Em  
MP,EX,3,3E10  i-13~Dk  
MP,PRXY,3,0.25    co]Gmg6p  
MP,DENS,3,2500 Zq?_dIX %  
TB,DP,3                                    !定义数据表3 Va9q`X
byO  
TBDATA,1,7.5E5,37.7                        !给数据表输入数据 KRk~w]  
!取强度折减系数F=1.4时边坡围岩材料属性 h(:<(o@<  
MP,EX,4,3E10  OlcP(  
MP,PRXY,4,0.25                              !定义泊松比 IRB& j%LA  
MP,DENS,4,2500 4]BJ0+|mT  
TB,DP,4                                    !定义数据表4 %-^}45](q  
TBDATA,1,6.4E5,33.5                        !给数据表输入数据 i,OKfXp  
!取强度折减系数F=1.6时边坡围岩材料属性 9/;{>RL=  
MP,EX,5,3E10  q+/l"&j.  
MP,PRXY,5,0.25    @\!!t{y  
MP,DENS,5,2500                            ！定义密度 gX/NtO%  
TB,DP,5                                    !定义数据表5 @InJ_
9E
 
TBDATA,1,5.6E5,30                          !给数据表输入数据 [@.B4p  
!取强度折减系数F=1.8时边坡围岩材料属性 KS! iL=i  
MP,EX,6,3E10                              ！定义弹性模量 k:0P+d  
MP,PRXY,6,0.25    Mvof%I  
MP,DENS,6,2500 mU]s7` %<>  
TB,DP,6                                    !定义数据表6 swZpWC  
TBDATA,1,5.0E5,27.2 r{"uv=,`  
!取强度折减系数F=2.0时边坡围岩材料属性 5#u.pu  
MP,EX,7,3E10  'Oy5e@G+?  
MP,PRXY,7,0.25                            !定义泊松比 3X'W
R]  
MP,DENS,7,2500 rt.[,m  
TB,DP,7                                    !定义数据表7 0eA5zFU7  
TBDATA,1,4.5E5,24.8 {E~l>Z88  
!取强度折减系数F=2.2时边坡围岩材料属性 b>=7B6 Aw  
MP,EX,8,3E10  9.8,q  
MP,PRXY,8,0.25    m3?e]nL4W  
MP,DENS,8,2500 DT? m/*  
TB,DP,8                                    !定义数据表8 8t``NZ[  
TBDATA,1,4.09E5,22.8 hDtKnF  
!取强度折减系数F=2.4时边坡围岩材料属性 %|?1B$s0  
MP,EX,9,3E10  _7 `E[&v  
MP,PRXY,9,0.25    !GNXt4D  
MP,DENS,9,2500 =HPu{K$  
TB,DP,9                                      !定义数据表9 1o#vhk/"+  
TBDATA,1,3.6E5,21.1                          !给数据表输入数据 a/e\vwHLv  
!取强度折减系数F=2.6时边坡围岩材料属性 R~c1)[[E  
MP,EX,10,3E10                                ！定义弹性模量  p":@>v?  
MP,PRXY,10,0.25                              !定义泊松比 Jk*QcEE=  
MP,DENS,10,2500                          ！定义密度 )k%M.{&bji  
TB,DP,10                                  !定义数据表10 #:W%,$9\P  
TBDATA,1,3.46E5,19.6                        !给数据表10输入数据 u9}!Gq  
!取强度折减系数F=2.8时边坡围岩材料属性 |Y{
PO&-?r  
MP,EX,11,3E10                              ！定义弹性模量 SKxe3  
MP,PRXY,11,0.25                            !定义泊松比 B!`\L!  
MP,DENS,11,2500 /+P5)q TKL  
TB,DP,11                                  !定义数据表11 QN #)F  
TBDATA,1,3.2E5,18.2                        !给数据表11输入数据 D=9}|b/  
!取强度折减系数F=3.0时边坡围岩材料属性 33M}>$ZH  
MP,EX,12,3E10                              ！定义弹性模量 V_M@g;<o  
MP,PRXY,12,0.25    q%.bnF/Yd  
MP,DENS,12,2500                            ！定义密度 MgnE-6_c  
TB,DP,12                                  !定义数据表12 4<yK7x  
TBDATA,1,3.0E5,17.1                        !给数据表12输入数据 w a.f![  
!建立几何模型 0WS|~?OR@  
!创建边坡线模型线 %MNk4UsV  
K,1,,,,                                      ！ 创建关键点1 BGpk&.J  
K,2,-800,,, ~^7  
K,3,-800,-800,,  uHrb:X!q  
K,4,-800,-1200,, ((9YG  
K,5,1200,-1200,,  TolrEcI  
K,6,1200,-800,,  [tN` :}?  
K,7,1200,,, 9Z9l:}bO  
K,8,1200,400,, bA+[{  
K,9,400,400,, .\4l'THn,0  
l,1,2                                        !连接关键点1、2生成直线L1 V85.DK!  
l,2,3 _S<?t9mS  
l,3,4 yM17H\=  
l,4,5 '?k' 6R$'\  
l,5,6 C38XQLC  
l,6,7 >Fh#DmQ  
l,7,8 NEg>lIu<~  
l,8,9 8_awMVAy  
l,9,1 IDmsz  
l,1,7 .nPL2zO  
l,3,6 ^je528%H  
!创建边坡面模型 ylim/`u}6  
al,3,4,5,11                              !由线L3、L4、L5、L11生成一个面A1 l'2H4W_+  
al,1,2,11,6,10                            !由线L1、L2、L1、L6、L10生成一个面A2 k!c7a\">{  
al,7,8,9,10                              !由线L7、L8、L9、L10生成一个面A3 y*|L:!   
/PNUM, area, on                        !打开面号开关 &?}1AQAYg  
!划分网格生成有限元模型 x~(y "^ph  
! 划分边坡围岩2网格 thQ J(w  
!设置网格份数 &G=0  
lsel,s,,,3,5,2                              !选择线L3和L5 )8]3kQffJ=  
lesize,all,,,5                              !把所选择线分为5段 0%;MVMH  
lsel,s,,,4,11,7                              !选择线L4和L11 kpT>G$s~gy  
lesize,all,,,26                              !把所选择线分为26段 W^|J/Y48  
mat,2                                    !给边坡围岩2赋予2号材料特性 g 2#F_  
type,1                                  !采用单元类型1 KDaN-r^{%
 
mshkey,0                                !设定自由网格划分 M\jB)@)  
mshape,0                                !设定四边形网格划分 4g'}h`kh  
amesh,1                                  ！划分面积A1 G|V\^.f<  
! 划分边坡围岩1网格 TMtI^mkB:  
!设置网格份数 (olLB  
lsel,s,,,2,6,4                              !选择线L2和L6 LH.%\TMN$  
lsel,a,,,1                                  !附加选择线L1 TPqvp|~2  
lesize,all,,,10                              !把所选择线分为10段 i0i`k^bA  
lsel,s,,,8,10,2                              !选择线L8和L10 %%x0w^  
lesize,all,,,16                              !把所选择线分为16段 .' IeHh  
lsel,s,,,7,9,2 r4S=I  
lesize,all,,,12 !3?~#e{_  
mat,1                                    !给边坡围岩1赋予1号材料特性 k) 3s
?  
type,1 6'vi68  
mshkey,0                                  !设定自由网格划分 uf?b%:A  
mshape,0                                  !设定四边形网格划分 R}.3|0  
amesh,2                                  !划分面积A2 Wa}"SqYr h  
amesh,3                                  ！划分面积A3 J;HkTT  
allsel TRCI\  
!保存网格模型 S]b xQa+  
save,Slope-grid.db HYFN?~G  
(3)施加约束和荷载 cS 4T\{B;  
!给边坡模型施加约束 g`.{K"N>!  
!边坡两侧施加X方向约束 u!u5g.Q  
nsel,s,loc,x,-800                            !选择X=-800线上所有节点 $$~a=q,P[  
nsel,a,loc,x,1200                            !选择X=1200线上所有节点 EFv4=OWB  
d,all,ux                                    !对所选择节点约束X方向位移 1!s!wQgS  
allsel :'ihE\j  
!边坡底部施加约束 x {vIT- f  
nsel,s,loc,y,-1200                          !选择Y=--1200线上所有节点 N=
<=dp(  
d,all,uy                                  !对所选择节点约束X、Y方向位移 +<B|qcT!  
d,all,ux                                  w?/f Zx  
!施加重力加速度 nO}$ 76*'0  
acel,,9.8 omT(3)TP  
(4)求解 *sAOpf@M  
/solu F%y{% C7l  
!求解设置
ytob/tc  
antype,static                                  !设定为静力求解 QP<FCmt8  
nsubst,100                                  !设定最大子步数为100 PuU
*vs3  
pred,on                                      !打开时间步长预测器 6.UKB<sV  
nropt,full                                    !设定牛顿-拉普森选项 Ir>2sTrm  
nlgeom,on                                    !打开大位移效果 1::LN(`<  
lnsrch,on                                    !打开线性搜索 U>1b9G"_  
outres,all,all                                  !输出所有项 ] EV`dIk  
cnvtol,f,,0.005,2,0.5                            !力收敛准则设定 (6o:4|xl0  
cnvtol,u,,0.05,2,1                              !位移收敛准则设定 cyHhy_~R  
!边坡在强度折减系数F=1时求解 @oA0{&G{  
allsel ( w5f(4  
solve                                        !进行求解 ,aYU$~o#  
save,F1,db                                    !把F=1时求解结果保存 t@r#b67WJe  
!边坡在强度折减系数F=1.2时求解 E7gL
~4I  
finish V:QdQ;c  
/solu ,-!2 5G  
allsel
`M6YblnJZ  
asel,s,area,,2,3,1                              !选择面积A2、A3 ef,F[-2^o  
MPCHG,3,all                                ! 把所选择单元材料号改为3 1zR/HT  
allsel Ki63Ox^O  
solve                                        !进行求解 P_mi)@  
save,F1.2,db                                  !把F=1.2时求解结果保存                                                                                  n8Q* _?Z/  
!边坡在强度折减系数F=1.4时求解 fYs?D+U;PF  
finish Bf+^O)Ns^  
/solu p&m ^
IWD  
allsel YjL t&D:IZ  
asel,s,area,,2,3,1                              !选择面积A2、A3 :Ip
~)n9t  
MPCHG,4,all                                !把所选择单元材料号改为4 W`5a:"Vg  
allsel b+_hI)T  
solve                                        !进行求解 GGp{b>E+ #  
save,F1.4,db                                  !把F=1.4时求解结果保存                            J[MVE4&  
!边坡在强度折减系数F=1.6时求解 0hb/`[Q  
finish 6w@,I;  
/solu 5C*?1& !  
allsel N@}gLBf  
asel,s,area,,2,3,1                              !选择面积A2、A3 CJ:uYXJJ:z  
MPCHG,5,all                                !把所选择单元材料号改为5 ]p}#NPe5  
allsel /xF 9:r  
solve " C&x,Ic  
save,F1.6,db                                  !把F=1.6时求解结果保存                                                                                  p^Ag
h  
!边坡在强度折减系数F=1.8时求解 >
07shN
X  
finish fvO;lA>`  
/solu >waN;&>/  
allsel M!-q}5';  
asel,s,area,,2,3,1                                !选择面积A2、A3 k5g@myb-  
MPCHG,6,all                                  !把所选择单元材料号改为6 "s> >V,  
allsel =_uol8v  
solve m.Twgin  
save,F1.8,db                                  !把F=1.8时求解结果保存  ?|)rv  
!边坡在强度折减系数F=2.0时求解 %L28$c3p  
finish VZIKjrKs  
/solu %db3f z  
allsel r /^'X
j'(  
asel,s,area,,2,3,1                                !选择面积A2、A3 ?&XzW+(X  
MPCHG,7,all                                  !把所选择单元材料号改为7 J/= +r0c  
allsel mUiOD$rO  
solve 3 ^pYCK%  
save,F2.0,db                                  !把F=2.0时求解结果保存                            8Y7 @D$=w  
!边坡在强度折减系数F=2.2时求解 :K:f^o]s  
finish {DSyV:  
/solu jB`7T^bU  
allsel 6G$/NW=L  
asel,s,area,,2,3,1                                !选择面积A2、A3 l-Fmn/V  
MPCHG,8,all                                  !把所选择单元材料号改为8 DYkC'+TEX  
allsel m
_(E(_  
solve ^b:Xo"q#H  
save,F2.2,db                                  !把F=2.2时求解结果保存                                                                                  wJ/~q)  
!边坡在强度折减系数F=2.4时求解 i5Eeg`NMl  
finish GIK u  
/solu F],TG&>5  
allsel >b3@>W  
asel,s,area,,2,3,1                                !选择面积A2、A3 2S'AIuIew
 
MPCHG,9,all                                  !把所选择单元材料号改为9 cu:-MpE  
allsel [:Xn6)qz  
solve {GAsFnZk  
save,F2.4,db                                  !把F=2.4时求解结果保存  ` v>
/  
!边坡在强度折减系数F=2.6时求解 $>EqH?EQ  
finish  w}"!l G  
/solu \A ;^ UxG  
allsel |E? ,xWN  
asel,s,area,,2,3,1                                !选择面积A2、A3 l,^xX=,  
MPCHG,10,all                                  !把所选择单元材料号改为10 |c=d;+  
allsel pAMo XJ`  
solve j_(?=7Y3g  
save,F2.6,db                                  !把F=2.6时求解结果保存                            >2nF"?"=  
!边坡在强度折减系数F=2.8时求解 \?r$&K]4  
finish  &Q<EfB  
/solu a4
:`2  
allsel Rnz8 f}  
asel,s,area,,2,3,1                                !选择面积A2、A3 f<^ScFVR  
MPCHG,11,all                                !把所选择单元材料号改为11 yg`E22  
allsel #0jSZg^,"  
solve /%-o.hT  
save,F2.8,db                                  !把F=2.8时求解结果保存                                                                                  M&eQ=vew.  
!边坡在强度折减系数F=3.0时求解 9A(n_Rs7?  
finish *1i?6$[ "  
/solu G]at{(^Vz  
allsel n}+wd9J*!2  
asel,s,area,,2,3,1                            !选择面积A2、A3 {a[&#Uv  
MPCHG,12,all                              !把所选择单元材料号改为12 ?-4OfGN  
allsel ?{?
Vy9'B  
solve B%)zGTp6  
save,F3.0,db                                !把F=3.0时求解结果保存 d8D yv#gT  
(5）后处理 QXsfp  
/post1                                      ！进入后处理 k:`a+LiZ  
!边坡在强度折减系数F=1时结果分析 +BU0 6lLD  
Resume,'F1','db'                              !读入边坡在强度折减系数F=1时 8u/3?Kc  
set,1,last                                    !读入后一个子步 |e~u!V\m  
pldisp,1                                      !绘制边坡模型变形图 LPb]mC6
#  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 >}70]dN7b  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 j1W bD7*8  
!边坡在强度折减系数F=1.2时结果分析 6|%^pjX5  
Resume,'F1.2','db'                            !读入边坡在强度折减系数F=1.2时 33O)k*g  
set,1,last                                    !读入后一个子步 ;B*L1'FF%t  
pldisp,1                                      !绘制边坡模型变形图 @Ap@m6K?q  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 =z+-l5Gu"  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 ^NY+wR5Sn  
!边坡在强度折减系数F=1.4时结果分析 JN-D/s  
Resume,'F1.4','db'                            !读入边坡在强度折减系数F=1.4时 <\+Po<)3j  
set,1,last                                    !读入后一个子步 bJANZn|H  
pldisp,1                                      !绘制边坡模型变形图 4$..r4@  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 H&w(]PDh  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 w4NZt|>5j;  
!边坡在强度折减系数F=1.6时结果分析 zI1(F67d`  
Resume,'F1.6','db'                            !读入边坡在强度折减系数F=1.6时 Q3(hK<Qh;  
set,1,last                                    !读入后一个子步 G,+xT}@wu  
pldisp,1                                      !绘制边坡模型变形图 d$4WK)U  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 FYxUOO  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 :fl*w"
"V@  
!边坡在强度折减系数F=1.8时结果分析 7 FEzak'  
Resume,'F1.8','db'                            !读入边坡在强度折减系数F=1.8时 md.*  
set,1,last                                    !读入后一个子步 PpW A f\  
pldisp,1                                      !绘制边坡模型变形图 "sdcP8])d  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 RA!x  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 <.;@ksCPW{  
!边坡在强度折减系数F=2.0时结果分析 sV]i/B  
Resume,'F2.0','db'                            !读入边坡在强度折减系数F=2.0时 #W_i{
bdO  
set,1,last                                    !读入后一个子步 i?lX,9%  
pldisp,1                                      !绘制边坡模型变形图 }<E sS  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 )J{.Cx<E  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 b?sAEU;  
!边坡在强度折减系数F=2.2时结果分析 \a\^(`3a[  
Resume,'F2.2','db'                            !读入边坡在强度折减系数F=2.2时 ZCj>MA  
set,1,last                                    !读入后一个子步 aeLBaS  
pldisp,1                                      !绘制边坡模型变形图 $_ST:h&C  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 lkI8{  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 \MfR #k0  
!边坡在强度折减系数F=2.4时结果分析 [^h/(a`  
Resume,'F2.4','db'                            !读入边坡在强度折减系数F=2.4时 |:~("rA+v  
set,1,last                                    !读入后一个子步 5dbX%e_OP  
pldisp,1                                      !绘制边坡模型变形图 e6Wl7&@6  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 t)YFT
O"Jj  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 b(g_.1[  
!边坡在强度折减系数F=2.6时结果分析 MoZU(j  
Resume,'F2.6','db'                            !读入边坡在强度折减系数F=2.6时 Ar\IZ_Q  
set,1,last                                    !读入后一个子步 GH[ U!J  
pldisp,1                                      !绘制边坡模型变形图 Bi'I18<  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 jn0t-":  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 J11dqj  
!边坡在强度折减系数F=2.8时结果分析 |G[
{{qZM5  
Resume,'F2.8','db'                            !读入边坡在强度折减系数F=2.8时 Pw0{.W~r  
set,1,last                                    !读入后一个子步 pHq{S;R2G  
pldisp,1                                      !绘制边坡模型变形图 pZyb  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 YhEiN. ~  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 Gj
G{qR  
!边坡在强度折减系数F=3.0时结果分析 @ARAX\F  
Resume,'F3.0','db'                            !读入边坡在强度折减系数F=3.0时 f<Va<TL6-  
set,1,last                                    !读入后一个子步 Sr4dY`V*:z  
pldisp,1                                      !绘制边坡模型变形图 z?8zFP  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 ]!n*V
/g  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 ['Hp?Q|k  
hz&^_G6`  
?IL! X-xx  
e{c._zr,  
;JMOsn}
8  

,)0/Ec  
/%2:+w