[讨论]有限元强度折减法 [复制链接]

各位朋友： 4>(
rskl_  
    大家好！我通过ansys用有限元强度折减法分析了一个边坡的稳定性，结果不同的工况下得出的塑性区是一样的，我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激！附命令流和塑性区如下：(1)创建物理环境 -$cO0RSY  
/COM, Structural                              ! 指定结构分析 5O"$'iL  
/TITLE, Slope stability Analysis                  ! 定义工作标题 w7QYWf'  
/FILNAM,support,1                            ! 定义工作文件名 #7p!xf^  
(2)建立模型 oR'u&\mB  
!进入前处理器 D7v_<  
/PREP7 ^D A<=C-[!  
!定义单元类型 s-JS[  
ET,1,PLANE82                              !定义边坡围岩单元 lHc9D  
KEYOPT,1,3,2                                !求解类型为平面应变 /G= ?E]^  
!定义材料属性 !p{
CsR8c  
!边坡围岩1材料属性  28LYGrB  
MP,EX,1,3E10  B PG&R  
MP,PRXY,1,0.25    Pd;Gc@'~  
MP,DENS,1,2500                            !定义塑性模型 0@kL<\u  
TB,DP,1                                    !定义数据表1 o4nDjFhh  
TBDATA,1,9E5,42.8                          !给数据表输入数据 :*WiswMFm  
!边坡围岩2材料属性                        !采用弹性模型 tYCVVs`?  
MP,EX,2,3.2E10                              ！定义弹性模量 Z^_gS&nDa~  
MP,PRXY,2,0.24                              !定义泊松比 [Lq9lw&   
MP,DENS,2,2700                            !定义密度 q0VR&b`?>D  
!取强度折减系数F=1.2时边坡围岩材料属性 F $1f8U8  
MP,EX,3,3E10  EjP9/VG@=  
MP,PRXY,3,0.25    ZhY03>X  
MP,DENS,3,2500 >- U+o.o  
TB,DP,3                                    !定义数据表3 ~
;ObT=  
TBDATA,1,7.5E5,37.7                        !给数据表输入数据 |X;|=
.  
!取强度折减系数F=1.4时边坡围岩材料属性 )F,z pGG  
MP,EX,4,3E10  %`}nP3  
MP,PRXY,4,0.25                              !定义泊松比 U[W &D%'  
MP,DENS,4,2500 dK>sHUu  
TB,DP,4                                    !定义数据表4 v:]z-zU  
TBDATA,1,6.4E5,33.5                        !给数据表输入数据 l;}3J3/qq]  
!取强度折减系数F=1.6时边坡围岩材料属性 W}@IUCRs  
MP,EX,5,3E10  7R$O~R3p  
MP,PRXY,5,0.25    jR
>`Xz  
MP,DENS,5,2500                            ！定义密度 Y]bS=*q  
TB,DP,5                                    !定义数据表5 #M@~8dAH}M  
TBDATA,1,5.6E5,30                          !给数据表输入数据 5Kw?#  
!取强度折减系数F=1.8时边坡围岩材料属性 ~{-9qOGw;  
MP,EX,6,3E10                              ！定义弹性模量 U;t1 K  
MP,PRXY,6,0.25    w$"^)EG,7  
MP,DENS,6,2500 kbZp
i`w  
TB,DP,6                                    !定义数据表6 .Ky)Co  
TBDATA,1,5.0E5,27.2 I %|;M%B  
!取强度折减系数F=2.0时边坡围岩材料属性 in`|.#  
MP,EX,7,3E10  ^o4](l  
MP,PRXY,7,0.25                            !定义泊松比 &1ZUMc  
MP,DENS,7,2500 'PWA  
TB,DP,7                                    !定义数据表7 @S1Z"%S  
TBDATA,1,4.5E5,24.8 NiD_v  
!取强度折减系数F=2.2时边坡围岩材料属性 'zOB!QqA`v  
MP,EX,8,3E10   Lr0:yo  
MP,PRXY,8,0.25    k5)a|  
MP,DENS,8,2500 G%viWWTY  
TB,DP,8                                    !定义数据表8 (@V_47o  
TBDATA,1,4.09E5,22.8 b*1yvkX5  
!取强度折减系数F=2.4时边坡围岩材料属性 q1Mt5O}  
MP,EX,9,3E10  m~-O
}i~)  
MP,PRXY,9,0.25    GI6]Ecc  
MP,DENS,9,2500 \&[(PNl  
TB,DP,9                                      !定义数据表9 LZ RP}|  
TBDATA,1,3.6E5,21.1                          !给数据表输入数据 ic}mru  
!取强度折减系数F=2.6时边坡围岩材料属性 k%VYAON  
MP,EX,10,3E10                                ！定义弹性模量  p4D.nB8  
MP,PRXY,10,0.25                              !定义泊松比 {@hJPK8  
MP,DENS,10,2500                          ！定义密度 8J:=@X^}  
TB,DP,10                                  !定义数据表10 % _nmv  
TBDATA,1,3.46E5,19.6                        !给数据表10输入数据 kLc@U~M  
!取强度折减系数F=2.8时边坡围岩材料属性 R]3j6
\  
MP,EX,11,3E10                              ！定义弹性模量 aNP\Q23D  
MP,PRXY,11,0.25                            !定义泊松比 d|>/eb.R  
MP,DENS,11,2500 `R!Q(rePx  
TB,DP,11                                  !定义数据表11 '3?-o|v@D  
TBDATA,1,3.2E5,18.2                        !给数据表11输入数据 nf1O8FwRb  
!取强度折减系数F=3.0时边坡围岩材料属性 wV-9T*QrM  
MP,EX,12,3E10                              ！定义弹性模量 <!F".9c@A  
MP,PRXY,12,0.25    #n]K$k>  
MP,DENS,12,2500                            ！定义密度 [:+f Y[4==  
TB,DP,12                                  !定义数据表12 TjHt:%7.  
TBDATA,1,3.0E5,17.1                        !给数据表12输入数据 MV]`[^xQ5  
!建立几何模型 2D/bMq  
!创建边坡线模型线 Xyjd7"  
K,1,,,,                                      ！ 创建关键点1 ),Hr  
K,2,-800,,, 3^5h:OaT  
K,3,-800,-800,,  pog   
K,4,-800,-1200,, NS-0-o|4#  
K,5,1200,-1200,,  ZsSW{ffZ77  
K,6,1200,-800,,  FmSE]et  
K,7,1200,,, _qk yU)z  
K,8,1200,400,, #x`K

4f)  
K,9,400,400,, |AS~sjWSJ  
l,1,2                                        !连接关键点1、2生成直线L1 b[<L l%K  
l,2,3 /B)2L]6p  
l,3,4 "z*.Bk  
l,4,5 _}F_Q5)  
l,5,6 %xr'96d  
l,6,7 _0UE*l$t  
l,7,8 =J|jCK[r  
l,8,9 ) ]DqK<-  
l,9,1 -[}Aka,f!  
l,1,7 d0R;|p''Z  
l,3,6 (,KzyR=*'  
!创建边坡面模型 6\k~q.U@XI  
al,3,4,5,11                              !由线L3、L4、L5、L11生成一个面A1 &hrMpD6z6i  
al,1,2,11,6,10                            !由线L1、L2、L1、L6、L10生成一个面A2 +\$c_9|C+  
al,7,8,9,10                              !由线L7、L8、L9、L10生成一个面A3 !{fu(E  
/PNUM, area, on                        !打开面号开关 ;YSe:m*  
!划分网格生成有限元模型 T}
/|nOu 5  
! 划分边坡围岩2网格 c-_1tSh}  
!设置网格份数 P+BGCc%);B  
lsel,s,,,3,5,2                              !选择线L3和L5 Kp^"<%RT  
lesize,all,,,5                              !把所选择线分为5段 ZMLN ;.{Na  
lsel,s,,,4,11,7                              !选择线L4和L11 s\<UDW  
lesize,all,,,26                              !把所选择线分为26段 {{?MO{Mh*  
mat,2                                    !给边坡围岩2赋予2号材料特性 |=07n K2  
type,1                                  !采用单元类型1 LcB+L](  
mshkey,0                                !设定自由网格划分 -xbs'[  
mshape,0                                !设定四边形网格划分 rT\~VJ>+i  
amesh,1                                  ！划分面积A1 ]>1`Fa6_  
! 划分边坡围岩1网格 ug.|ag'R  
!设置网格份数 g/}d> 6  
lsel,s,,,2,6,4                              !选择线L2和L6 "?<(-,T  
lsel,a,,,1                                  !附加选择线L1 vKWi?}1  
lesize,all,,,10                              !把所选择线分为10段 K1o>>388G  
lsel,s,,,8,10,2                              !选择线L8和L10 l(Dr@LB~  
lesize,all,,,16                              !把所选择线分为16段 ay|jq"a  
lsel,s,,,7,9,2 iJj!-a:z.  
lesize,all,,,12 R!yh0y}Z  
mat,1                                    !给边坡围岩1赋予1号材料特性 "a9j2+9  
type,1 @,7r<6E  
mshkey,0                                  !设定自由网格划分 EV-sEl8ki  
mshape,0                                  !设定四边形网格划分 _>BYUPY  
amesh,2                                  !划分面积A2 { ^^5FE)%  
amesh,3                                  ！划分面积A3 #!E`%' s]  
allsel &n6L;y-  
!保存网格模型 dm]g:KWg  
save,Slope-grid.db JdLPIfI^  
(3)施加约束和荷载 pL!,1D!  
!给边坡模型施加约束 v 2p  
!边坡两侧施加X方向约束 (P;TM1k  
nsel,s,loc,x,-800                            !选择X=-800线上所有节点 bRrSd:e  
nsel,a,loc,x,1200                            !选择X=1200线上所有节点 `JY+3d,Ui  
d,all,ux                                    !对所选择节点约束X方向位移 v_Df+  
allsel }V*
?~.R  
!边坡底部施加约束 #Hz9@H  
nsel,s,loc,y,-1200                          !选择Y=--1200线上所有节点 zA}JVB  
d,all,uy                                  !对所选择节点约束X、Y方向位移 Z6rZAwy  
d,all,ux                                  1zC
u1'Wv  
!施加重力加速度 Wp+lI1t  
acel,,9.8 @$!6u0x  
(4)求解 P3-O)m]jv  
/solu o.w/?  
!求解设置 SP/b4  
antype,static                                  !设定为静力求解 ?iV}U  
nsubst,100                                  !设定最大子步数为100 m mZP;  
pred,on                                      !打开时间步长预测器 h Ypj  
nropt,full                                    !设定牛顿-拉普森选项 {&XTa`C  
nlgeom,on                                    !打开大位移效果 tzfyS#E  
lnsrch,on                                    !打开线性搜索 B9[vv;lzu  
outres,all,all                                  !输出所有项 M$.bC0}T  
cnvtol,f,,0.005,2,0.5                            !力收敛准则设定 60]VOQku  
cnvtol,u,,0.05,2,1                              !位移收敛准则设定 YtKT3u:x  
!边坡在强度折减系数F=1时求解 pUS:HJk|  
allsel 7)[2Ud8  
solve                                        !进行求解 uF1 4;  
save,F1,db                                    !把F=1时求解结果保存 UJQTArf  
!边坡在强度折减系数F=1.2时求解 I'^XEl?  
finish }st~$JsV1  
/solu I\1"E y  
allsel mtkZF{3Jx  
asel,s,area,,2,3,1                              !选择面积A2、A3 M$Ui=GGq  
MPCHG,3,all                                ! 把所选择单元材料号改为3 ]kx<aQ^  
allsel ']fyD3N  
solve                                        !进行求解 G6qZ>-GiL  
save,F1.2,db                                  !把F=1.2时求解结果保存                                                                                  j,;f#+O`g  
!边坡在强度折减系数F=1.4时求解 J%|;  
finish )/JVp>  
/solu 8t=O=l\  
allsel /4O
Qx0Xmm  
asel,s,area,,2,3,1                              !选择面积A2、A3  B9y5NX  
MPCHG,4,all                                !把所选择单元材料号改为4 9H;Os:"\|  
allsel }yn%_KQ0  
solve                                        !进行求解 gK;dfrU.8Y  
save,F1.4,db                                  !把F=1.4时求解结果保存                            qoH:_o8ClO  
!边坡在强度折减系数F=1.6时求解 kTfRm^  
finish X@}7 #Vt  
/solu -?:8sv*X  
allsel 1Az&BZU[  
asel,s,area,,2,3,1                              !选择面积A2、A3 5+!yXkE^e  
MPCHG,5,all                                !把所选择单元材料号改为5 Pv,PS.,-  
allsel E2d
S@!]V  
solve jD"nEp-  
save,F1.6,db                                  !把F=1.6时求解结果保存                                                                                  p7Zeudmj  
!边坡在强度折减系数F=1.8时求解 1%vE7a>{  
finish Sz<:WY/(x  
/solu Gey-8  
allsel p/Q< VV  
asel,s,area,,2,3,1                                !选择面积A2、A3 V"(5U(v{~  
MPCHG,6,all                                  !把所选择单元材料号改为6 -T1R}ew*t  
allsel l3BN,HNv+  
solve u/wX7s  
save,F1.8,db                                  !把F=1.8时求解结果保存  s.rQiD  
!边坡在强度折减系数F=2.0时求解 xzA!,75@U  
finish #o[n.  
/solu h$$JXf  
allsel .sQV0jF{  
asel,s,area,,2,3,1                                !选择面积A2、A3 !`7evV:  
MPCHG,7,all                                  !把所选择单元材料号改为7 FN/l/OSb  
allsel k$m'ebrS.~  
solve l l*g *zt3  
save,F2.0,db                                  !把F=2.0时求解结果保存                            +PWm=;tcC  
!边坡在强度折减系数F=2.2时求解 :|S[i('  
finish yK"\~t[@X:  
/solu \'u+iB g  
allsel 3I)oqS@q'  
asel,s,area,,2,3,1                                !选择面积A2、A3 I4w``""c  
MPCHG,8,all                                  !把所选择单元材料号改为8  0%,W5w  
allsel FZ<6kk4  
solve ^`dMjeF  
save,F2.2,db                                  !把F=2.2时求解结果保存                                                                                  *oIIcE4g7  
!边坡在强度折减系数F=2.4时求解 0S;Ipg  
finish Fw(  
/solu ]AM*9!  
allsel ws,?I
mA  
asel,s,area,,2,3,1                                !选择面积A2、A3 tj00xYY  
MPCHG,9,all                                  !把所选择单元材料号改为9 H|aC(c  
allsel (zy|>u  
solve G7,v:dlK   
save,F2.4,db                                  !把F=2.4时求解结果保存  7b-[# g  
!边坡在强度折减系数F=2.6时求解 YqXN|&  
finish }j1;0kb?  
/solu 4IB`7QJq  
allsel 9;vES^  
asel,s,area,,2,3,1                                !选择面积A2、A3 i$3#/*Y7_L  
MPCHG,10,all                                  !把所选择单元材料号改为10 jqj}j2 9  
allsel 8KigGhY'ms  
solve +/%4E %  
save,F2.6,db                                  !把F=2.6时求解结果保存                            G.iQ\'1_h  
!边坡在强度折减系数F=2.8时求解 MFO%F) 5  
finish )>b1%x} =  
/solu Sh-B!  
allsel Z ]ZUK  
asel,s,area,,2,3,1                                !选择面积A2、A3 K*'AjT9wX+  
MPCHG,11,all                                !把所选择单元材料号改为11 NcwUK\  
allsel XPq`;<G  
solve [:e>FXV  
save,F2.8,db                                  !把F=2.8时求解结果保存                                                                                  y6sY?uu  
!边坡在强度折减系数F=3.0时求解 w^HI lA  
finish bOrE86v:  
/solu bT9:9LP  
allsel S\sy^Kt~4:  
asel,s,area,,2,3,1                            !选择面积A2、A3 y|*4XF<b  
MPCHG,12,all                              !把所选择单元材料号改为12 y,Bj,zw  
allsel L
{&1w  
solve K)`R?CZ:s  
save,F3.0,db                                !把F=3.0时求解结果保存 =? q&/ cru  
(5）后处理 <?8cVLW}O  
/post1                                      ！进入后处理 d/3&3>/  
!边坡在强度折减系数F=1时结果分析 wod{C!  
Resume,'F1','db'                              !读入边坡在强度折减系数F=1时 >.C$2bW<L  
set,1,last                                    !读入后一个子步 r z@%rOWV  
pldisp,1                                      !绘制边坡模型变形图 RiZ}cd  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 Qd% (]L[N.  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 jL4"FTcE]3  
!边坡在强度折减系数F=1.2时结果分析 P&5vVA6K7  
Resume,'F1.2','db'                            !读入边坡在强度折减系数F=1.2时
#q0xlF@  
set,1,last                                    !读入后一个子步 GO][`zZJ]  
pldisp,1                                      !绘制边坡模型变形图 XM?c*,=fu  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 i ^N}avO  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 Ly, ];  
!边坡在强度折减系数F=1.4时结果分析 {O!;cI~  
Resume,'F1.4','db'                            !读入边坡在强度折减系数F=1.4时 ^dxy%*Z/  
set,1,last                                    !读入后一个子步 Kb5}M/8  
pldisp,1                                      !绘制边坡模型变形图 "4smW>f:%  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 j`3IizN2  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 o0b\<}  
!边坡在强度折减系数F=1.6时结果分析 @N>rOA  
Resume,'F1.6','db'                            !读入边坡在强度折减系数F=1.6时 UQ^ )t ]  
set,1,last                                    !读入后一个子步 jl]p e7-  
pldisp,1                                      !绘制边坡模型变形图 >/@Q7V99{  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 B1i'Mzm-4  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 A"+t[0$.  
!边坡在强度折减系数F=1.8时结果分析 (lit^v,9  
Resume,'F1.8','db'                            !读入边坡在强度折减系数F=1.8时 )F'hn+(B|G  
set,1,last                                    !读入后一个子步 a
hM?;p  
pldisp,1                                      !绘制边坡模型变形图 L>Y3t1=  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 YNB7`:  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 j"s7P%  
!边坡在强度折减系数F=2.0时结果分析 j8G$,~v  
Resume,'F2.0','db'                            !读入边坡在强度折减系数F=2.0时 l
u?:1V-  
set,1,last                                    !读入后一个子步 Y3\EX  
pldisp,1                                      !绘制边坡模型变形图 s&4&\Aq}x#  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 *Fg)`M3g  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 7w<e^H?  
!边坡在强度折减系数F=2.2时结果分析 i5,yrPF  
Resume,'F2.2','db'                            !读入边坡在强度折减系数F=2.2时 iYf)FPET  
set,1,last                                    !读入后一个子步 
#Dea$  
pldisp,1                                      !绘制边坡模型变形图 3L CT-rp  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 L)n_  Q  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 TVM19
)9  
!边坡在强度折减系数F=2.4时结果分析 <N:)Xf9`  
Resume,'F2.4','db'                            !读入边坡在强度折减系数F=2.4时 ?Rk[P cX<  
set,1,last                                    !读入后一个子步 SsEpuEn  
pldisp,1                                      !绘制边坡模型变形图 ICEyz| C  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 }BUm}.-{u,  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 P6I<M}p  
!边坡在强度折减系数F=2.6时结果分析 Yr*!
T= z  
Resume,'F2.6','db'                            !读入边坡在强度折减系数F=2.6时 S"t\LB*'Ls  
set,1,last                                    !读入后一个子步 1=h5Z3/fj  
pldisp,1                                      !绘制边坡模型变形图 KO\-|#3y>  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 ' GUCXx  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 v5 @9  
!边坡在强度折减系数F=2.8时结果分析
wmA TV/  
Resume,'F2.8','db'                            !读入边坡在强度折减系数F=2.8时 :}R,a=N  
set,1,last                                    !读入后一个子步 m1e Sn |)7  
pldisp,1                                      !绘制边坡模型变形图 )<f4F!?,A  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 
0R,.  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 ["#H/L]3  
!边坡在强度折减系数F=3.0时结果分析 *10qP?0H  
Resume,'F3.0','db'                            !读入边坡在强度折减系数F=3.0时 Om*(dK]zHQ  
set,1,last                                    !读入后一个子步 c
*y*UG  
pldisp,1                                      !绘制边坡模型变形图 D4N(FZ
0~  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 73_=CP"t  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图

楼主兄弟，Ansys用于岩土工程计算本来就不是很合适，更何况采用强度折减法啊，目前流行的这种思路本人不是很认同。