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

各位朋友： uF<\|y rFt  
    大家好！我通过ansys用有限元强度折减法分析了一个边坡的稳定性，结果不同的工况下得出的塑性区是一样的，我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激！附命令流和塑性区如下：(1)创建物理环境 BD&AtOj[,  
/COM, Structural                              ! 指定结构分析 Fz^5cxmw  
/TITLE, Slope stability Analysis                  ! 定义工作标题 V5S6?V\  
/FILNAM,support,1                            ! 定义工作文件名 8QN/D\uq  
(2)建立模型 i?|b:lcV  
!进入前处理器 =[IKwmCX  
/PREP7 -'RD%_  
!定义单元类型 \bv JZ_  
ET,1,PLANE82                              !定义边坡围岩单元 8o[+>W  
KEYOPT,1,3,2                                !求解类型为平面应变 9[Xe|5?c  
!定义材料属性 :[bpMP<bz;  
!边坡围岩1材料属性  drh,=M\F  
MP,EX,1,3E10  0<42\ya  
MP,PRXY,1,0.25    gutf[Ksu  
MP,DENS,1,2500                            !定义塑性模型 ~ ve  
TB,DP,1                                    !定义数据表1 r,cK#!<%  
TBDATA,1,9E5,42.8                          !给数据表输入数据 _QiGrC  
!边坡围岩2材料属性                        !采用弹性模型 4\(|V fy  
MP,EX,2,3.2E10                              ！定义弹性模量 BaOPtBYA:  
MP,PRXY,2,0.24                              !定义泊松比 AqjEz+TVt  
MP,DENS,2,2700                            !定义密度
y.?Q  
!取强度折减系数F=1.2时边坡围岩材料属性 \\$wg   
MP,EX,3,3E10  0t.v  
MP,PRXY,3,0.25    p@%H. 5&&  
MP,DENS,3,2500 uAv'%/  
TB,DP,3                                    !定义数据表3 <M M(Z  
TBDATA,1,7.5E5,37.7                        !给数据表输入数据 I0(nRu<  
!取强度折减系数F=1.4时边坡围岩材料属性 VpWpC&  
MP,EX,4,3E10  `&g1`vg  
MP,PRXY,4,0.25                              !定义泊松比 HQ`A.E2  
MP,DENS,4,2500 `lN Z|U  
TB,DP,4                                    !定义数据表4 f^ 6da6Z  
TBDATA,1,6.4E5,33.5                        !给数据表输入数据 );L+)UV  
!取强度折减系数F=1.6时边坡围岩材料属性 ^LAdN8Cbb  
MP,EX,5,3E10  ;#S4$wISw`  
MP,PRXY,5,0.25    <k7q9"\4  
MP,DENS,5,2500                            ！定义密度 LGPg\g`  
TB,DP,5                                    !定义数据表5 HOlMj!.  
TBDATA,1,5.6E5,30                          !给数据表输入数据 `g:bvIV5x>  
!取强度折减系数F=1.8时边坡围岩材料属性 |5me }!C  
MP,EX,6,3E10                              ！定义弹性模量 A&=`?4>  
MP,PRXY,6,0.25    onF?;>[  
MP,DENS,6,2500 Pc=:j(  
TB,DP,6                                    !定义数据表6  ){xMMQ5  
TBDATA,1,5.0E5,27.2 & 6~AY:0r  
!取强度折减系数F=2.0时边坡围岩材料属性 ~ ]^<*R  
MP,EX,7,3E10  +V/mV7FK  
MP,PRXY,7,0.25                            !定义泊松比 }
BLT2]y0  
MP,DENS,7,2500 6Y/TqI[   
TB,DP,7                                    !定义数据表7 }7Y@u
@R  
TBDATA,1,4.5E5,24.8 lBfG#\rdW~  
!取强度折减系数F=2.2时边坡围岩材料属性 6x"|,,&MD0  
MP,EX,8,3E10  &z%DX   
MP,PRXY,8,0.25    u
U#e54^  
MP,DENS,8,2500 0
xpE+GY  
TB,DP,8                                    !定义数据表8 e(Ub7
L#  
TBDATA,1,4.09E5,22.8 rI4N3d;C  
!取强度折减系数F=2.4时边坡围岩材料属性 L+TM3*a*  
MP,EX,9,3E10  zq4)Uab*  
MP,PRXY,9,0.25    <PXnR\  
MP,DENS,9,2500 JURJN+)z  
TB,DP,9                                      !定义数据表9 xZp`Ke
!  
TBDATA,1,3.6E5,21.1                          !给数据表输入数据 #(d/A<  
!取强度折减系数F=2.6时边坡围岩材料属性 #{|
F2AM  
MP,EX,10,3E10                                ！定义弹性模量  @4%x7%+[c  
MP,PRXY,10,0.25                              !定义泊松比
HD9+4~8  
MP,DENS,10,2500                          ！定义密度 Q)l~?Fx  
TB,DP,10                                  !定义数据表10 #GA6vJ4^s  
TBDATA,1,3.46E5,19.6                        !给数据表10输入数据 H"%SzU  
!取强度折减系数F=2.8时边坡围岩材料属性 :qO)^~x  
MP,EX,11,3E10                              ！定义弹性模量 6%2\bI.#  
MP,PRXY,11,0.25                            !定义泊松比 `&qeSEs\  
MP,DENS,11,2500 J7s\  
TB,DP,11                                  !定义数据表11 C_ (s  
TBDATA,1,3.2E5,18.2                        !给数据表11输入数据 4dUr8]BkG  
!取强度折减系数F=3.0时边坡围岩材料属性 vm`\0VGSW  
MP,EX,12,3E10                              ！定义弹性模量 E>w|i  
MP,PRXY,12,0.25    v#Y
9O6g]T  
MP,DENS,12,2500                            ！定义密度 r`!S*zK  
TB,DP,12                                  !定义数据表12 ,P$Crs[  
TBDATA,1,3.0E5,17.1                        !给数据表12输入数据 lr&O@ 5"oy  
!建立几何模型 7;H P_oAu  
!创建边坡线模型线 $Y_v X 2  
K,1,,,,                                      ！ 创建关键点1 ulxy 4] h  
K,2,-800,,, s14;\  
K,3,-800,-800,,  \_PD@A9  
K,4,-800,-1200,, &g\?znF]H  
K,5,1200,-1200,,  WU<C7  
K,6,1200,-800,,  b5d;_-~d  
K,7,1200,,, r[y3@SE5  
K,8,1200,400,, 50^T\u  
K,9,400,400,, iJ 8I# j+N  
l,1,2                                        !连接关键点1、2生成直线L1 \[;Q
qn0  
l,2,3 3M<T}>  
l,3,4 7P7d[KP<  
l,4,5 i 79;;9M  
l,5,6 8WL*Pr1I  
l,6,7 ,?Nc\Q<:  
l,7,8 5sK1rDN  
l,8,9 8i'EO6  
l,9,1 DJ<F8-sb2r  
l,1,7 %!QY:[  
l,3,6  *"K7<S[  
!创建边坡面模型 'Z ,T,zW  
al,3,4,5,11                              !由线L3、L4、L5、L11生成一个面A1 JB
vP {5  
al,1,2,11,6,10                            !由线L1、L2、L1、L6、L10生成一个面A2 Z*Jp?[##  
al,7,8,9,10                              !由线L7、L8、L9、L10生成一个面A3 +q@g  
/PNUM, area, on                        !打开面号开关 Yeb-u+23  
!划分网格生成有限元模型 ctWH?b/ua  
! 划分边坡围岩2网格 x\2N @*I:  
!设置网格份数 u!K5jqP  
lsel,s,,,3,5,2                              !选择线L3和L5 xlU:&=|  
lesize,all,,,5                              !把所选择线分为5段 =}Xw}X+[WY  
lsel,s,,,4,11,7                              !选择线L4和L11 #Yr/GNN  
lesize,all,,,26                              !把所选择线分为26段 3rQ;}<*M  
mat,2                                    !给边坡围岩2赋予2号材料特性 CXfPC[o  
type,1                                  !采用单元类型1 3QO*1P@q  
mshkey,0                                !设定自由网格划分
-p2 =?a  
mshape,0                                !设定四边形网格划分 f+j-M|A  
amesh,1                                  ！划分面积A1 hp)k[|u;  
! 划分边坡围岩1网格 g$$j:U*-  
!设置网格份数 {[Vkht}  
lsel,s,,,2,6,4                              !选择线L2和L6 b<?A  
lsel,a,,,1                                  !附加选择线L1 ? {vY3~  
lesize,all,,,10                              !把所选择线分为10段 Ve\=By-a|  
lsel,s,,,8,10,2                              !选择线L8和L10 1!`B8y)  
lesize,all,,,16                              !把所选择线分为16段 ums*EKjs97  
lsel,s,,,7,9,2 j,i> 1|J  
lesize,all,,,12  {]=oOy1  
mat,1                                    !给边坡围岩1赋予1号材料特性 b^I(>l-  
type,1 8ECBi(  
mshkey,0                                  !设定自由网格划分 8WvQ[cd  
mshape,0                                  !设定四边形网格划分 %44Z7  
amesh,2                                  !划分面积A2 WjsE#9D!of  
amesh,3                                  ！划分面积A3 g_F-PT>($  
allsel *^b<CZd9  
!保存网格模型 ;fnE"}  
save,Slope-grid.db lH8e?zJ  
(3)施加约束和荷载 8{iFxTz  
!给边坡模型施加约束 u*i[A\Y  
!边坡两侧施加X方向约束 N J_#;t#j  
nsel,s,loc,x,-800                            !选择X=-800线上所有节点 wSP'pM{#2  
nsel,a,loc,x,1200                            !选择X=1200线上所有节点 [j^c&}0  
d,all,ux                                    !对所选择节点约束X方向位移 _ BUD~'Q5  
allsel G0VbW-`O  
!边坡底部施加约束 _ZU.;0  
nsel,s,loc,y,-1200                          !选择Y=--1200线上所有节点 = 7TK&  
d,all,uy                                  !对所选择节点约束X、Y方向位移 Fi!XaO  
d,all,ux                                  lf%Ju$H   
!施加重力加速度 /6Vn WrN_  
acel,,9.8 ]v{TSP^/  
(4)求解 js<}>wD7<  
/solu  y+.E}  
!求解设置 q"sD>Yh&  
antype,static                                  !设定为静力求解 BavGirCp  
nsubst,100                                  !设定最大子步数为100 R" '=^  
pred,on                                      !打开时间步长预测器 :k*3?*'K  
nropt,full                                    !设定牛顿-拉普森选项 7y2-8eL  
nlgeom,on                                    !打开大位移效果 L-v-KO6  
lnsrch,on                                    !打开线性搜索 fPLi8`r  
outres,all,all                                  !输出所有项 _onEXrM  
cnvtol,f,,0.005,2,0.5                            !力收敛准则设定 >s+TD4OfY  
cnvtol,u,,0.05,2,1                              !位移收敛准则设定 mrvPzoF,]  
!边坡在强度折减系数F=1时求解 iYHCa }  
allsel +rA:/!b)Y  
solve                                        !进行求解 3A5:D#  
save,F1,db                                    !把F=1时求解结果保存 a="\?L5  
!边坡在强度折减系数F=1.2时求解 4RYvI!
 
finish :i>/aRNh1  
/solu t<QSp6n""  
allsel 6EeK5XLf,  
asel,s,area,,2,3,1                              !选择面积A2、A3 tQ > IJ  
MPCHG,3,all                                ! 把所选择单元材料号改为3 ",&c"r4c  
allsel g=)djXW  
solve                                        !进行求解 AJ`R2 $  
save,F1.2,db                                  !把F=1.2时求解结果保存                                                                                  UAi]hUq  
!边坡在强度折减系数F=1.4时求解 540,A,>:tb  
finish Sd0y=!Pj=  
/solu 7,![oY[  
allsel 5o dtYI%L  
asel,s,area,,2,3,1                              !选择面积A2、A3 wmf#3"n  
MPCHG,4,all                                !把所选择单元材料号改为4 jLLZZPBK  
allsel +S3r]D3v/  
solve                                        !进行求解 +,BJ4``*k  
save,F1.4,db                                  !把F=1.4时求解结果保存                            n-Qpg  
!边坡在强度折减系数F=1.6时求解 _x ;fTW0  
finish OY>0qj  
/solu 'K0=FPB/@  
allsel Rw9 *!<Izt  
asel,s,area,,2,3,1                              !选择面积A2、A3 `LID*uD;_  
MPCHG,5,all                                !把所选择单元材料号改为5 R?K[O   
allsel [)&(zJHX  
solve > l@o\  
save,F1.6,db                                  !把F=1.6时求解结果保存                                                                                  wK[Xm'QTPJ  
!边坡在强度折减系数F=1.8时求解 U;Ne"Jh  
finish Q:4euhz*  
/solu Q|`sYm'.  
allsel ;0!rq^JG  
asel,s,area,,2,3,1                                !选择面积A2、A3 {_{&t>s2  
MPCHG,6,all                                  !把所选择单元材料号改为6 cqyrao3;  
allsel Ao/KB_4f*Q  
solve aAX(M=3  
save,F1.8,db                                  !把F=1.8时求解结果保存  s[HQq;S  
!边坡在强度折减系数F=2.0时求解 u(`,7 o "  
finish Ea7LPHE#  
/solu 4xE [S  
allsel 7P2?SW^
 
asel,s,area,,2,3,1                                !选择面积A2、A3 z2GT9  
MPCHG,7,all                                  !把所选择单元材料号改为7 MCcWRbE5#  
allsel ,n&e,I  
solve B- VhUS  
save,F2.0,db                                  !把F=2.0时求解结果保存                            qAF.i^  
!边坡在强度折减系数F=2.2时求解 b&$sY!iU  
finish [qEd`8V(  
/solu ~!Q\\_  
allsel lN-[2vT<  
asel,s,area,,2,3,1                                !选择面积A2、A3 *ydU3LG7  
MPCHG,8,all                                  !把所选择单元材料号改为8 Vu`O%[Q/  
allsel XtV=Gr8"  
solve ]1XJQW@gF  
save,F2.2,db                                  !把F=2.2时求解结果保存                                                                                  H)
${"  
!边坡在强度折减系数F=2.4时求解 #ir~v>J||  
finish 0R0j7\{  
/solu v'QmuMWF  
allsel jPjFp35;zb  
asel,s,area,,2,3,1                                !选择面积A2、A3 I>jDM  
MPCHG,9,all                                  !把所选择单元材料号改为9 ?\l@k(w4[x  
allsel ]5=
C3Y  
solve l]GUQcN=  
save,F2.4,db                                  !把F=2.4时求解结果保存  ?z2k74&M^  
!边坡在强度折减系数F=2.6时求解 qL03iV#h*V  
finish 8@f=GJf  
/solu e
{dYLQd  
allsel )|`#BC  
asel,s,area,,2,3,1                                !选择面积A2、A3 ny.YkN2  
MPCHG,10,all                                  !把所选择单元材料号改为10 4X5Tyv(Dp  
allsel EZ.|6oug\  
solve y_=},a  
save,F2.6,db                                  !把F=2.6时求解结果保存                            u\JYxNj1  
!边坡在强度折减系数F=2.8时求解 MJ
)aY2  
finish qrj:H4#VB  
/solu %z_PEqRj  
allsel fs=W(~"  
asel,s,area,,2,3,1                                !选择面积A2、A3 -0{"QhdE%  
MPCHG,11,all                                !把所选择单元材料号改为11 \R0&*cnmo  
allsel Y_}DF.>I P  
solve 9Xu O\+z  
save,F2.8,db                                  !把F=2.8时求解结果保存                                                                                  ;U^7]JO;  
!边坡在强度折减系数F=3.0时求解 abVz/R/o  
finish `XJU$c  
/solu V_SZp8  
allsel 
jd&kak
 
asel,s,area,,2,3,1                            !选择面积A2、A3 A{!D7kwTz~  
MPCHG,12,all                              !把所选择单元材料号改为12 !P6\-.  
allsel v/Z!Wp1LV  
solve Zu$30&U  
save,F3.0,db                                !把F=3.0时求解结果保存 j;|rI`67~  
(5）后处理 f0^DsP  
/post1                                      ！进入后处理 iYyJq;
S   
!边坡在强度折减系数F=1时结果分析 G%V*+Ond  
Resume,'F1','db'                              !读入边坡在强度折减系数F=1时 ^@&RJa-kb  
set,1,last                                    !读入后一个子步 U#jbii6e  
pldisp,1                                      !绘制边坡模型变形图 d`_X$P4y  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 wjr1?c  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 "t{|e6   
!边坡在强度折减系数F=1.2时结果分析 fgg;WXcT ~  
Resume,'F1.2','db'                            !读入边坡在强度折减系数F=1.2时  -<'&"-  
set,1,last                                    !读入后一个子步 >4zH\T!  
pldisp,1                                      !绘制边坡模型变形图 Dm"@59x  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 *W#_W]Tu  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 AvN\^ &G  
!边坡在强度折减系数F=1.4时结果分析 V ?10O  
Resume,'F1.4','db'                            !读入边坡在强度折减系数F=1.4时 fFHT`"bD:  
set,1,last                                    !读入后一个子步 ~;f,Ad`Q  
pldisp,1                                      !绘制边坡模型变形图 } h.]sF  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 fh1rmet&Ts  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 t/=xY'7  
!边坡在强度折减系数F=1.6时结果分析 UY
_'F5X  
Resume,'F1.6','db'                            !读入边坡在强度折减系数F=1.6时 !1:364  
set,1,last                                    !读入后一个子步 {hr+ENgV  
pldisp,1                                      !绘制边坡模型变形图 U(.3[x  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 0;b%@_E  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 Felu`@b  
!边坡在强度折减系数F=1.8时结果分析 9Okb)K95  
Resume,'F1.8','db'                            !读入边坡在强度折减系数F=1.8时 oWZbfR9R  
set,1,last                                    !读入后一个子步 483BrFV  
pldisp,1                                      !绘制边坡模型变形图 U-P\F-  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 gU
oL8~  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 pMB~Lt9  
!边坡在强度折减系数F=2.0时结果分析 5df~] -=0Y  
Resume,'F2.0','db'                            !读入边坡在强度折减系数F=2.0时 {~"&$DY2  
set,1,last                                    !读入后一个子步 7h4"5GlO0  
pldisp,1                                      !绘制边坡模型变形图 kT!Y~c  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 * 7
ki$f!  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 RaS7IL:e  
!边坡在强度折减系数F=2.2时结果分析 | 'SqG}h  
Resume,'F2.2','db'                            !读入边坡在强度折减系数F=2.2时 uKI2KWU?2  
set,1,last                                    !读入后一个子步 6QCU:2IiL  
pldisp,1                                      !绘制边坡模型变形图 `XwFH#_  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 KT
)A{i  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 (Ut)APM  
!边坡在强度折减系数F=2.4时结果分析 FQbF)K~e  
Resume,'F2.4','db'                            !读入边坡在强度折减系数F=2.4时 +$eEZ;4  
set,1,last                                    !读入后一个子步 f$lf(brQ:  
pldisp,1                                      !绘制边坡模型变形图 Ol,Tw=?  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 .,C8ASfh  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 }}";)}C`  
!边坡在强度折减系数F=2.6时结果分析 y] Io`w(>  
Resume,'F2.6','db'                            !读入边坡在强度折减系数F=2.6时 -NG`mfu  
set,1,last                                    !读入后一个子步 BwN
65_5p  
pldisp,1                                      !绘制边坡模型变形图 PQmgv&!DP  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 K%) K$/A  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 _?M71>3$.  
!边坡在强度折减系数F=2.8时结果分析
'NM$<<0  
Resume,'F2.8','db'                            !读入边坡在强度折减系数F=2.8时 +v 9@du  
set,1,last                                    !读入后一个子步 n]/7UH}(<&  
pldisp,1                                      !绘制边坡模型变形图 (z}q6Lfa  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 DQ{Yr>
J  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 >f [Lb|t  
!边坡在强度折减系数F=3.0时结果分析 6#/Riu%  
Resume,'F3.0','db'                            !读入边坡在强度折减系数F=3.0时 L}bS"=B[&W  
set,1,last                                    !读入后一个子步 ,qj  
pldisp,1                                      !绘制边坡模型变形图 3H0~?z_  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 9B
lc  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图

如果你这个能算出正确的结果的话,才真的叫出错,......... ^^[MDjNy@  
>&K1+FSmyJ  
x)M=_u2 _  
2k,!P6fgl  
Mf0XQ3n`H  
)q?z"F|  
#!wL0p  
~ {sRK  
     

好强啊，楼主能帮忙告诉我一下怎样才能输出网格上的节点的坐标呢？谢谢楼主，我写论文急用

不知道二楼的说什么

高人现身啊

做的不错嘛

顶！！！！！！！！！！！

晕  错误很明显呀 比如：Resume,'F3.0','db'                            !读入边坡在强度折减系数F=3.0时 nfck3h  
set,1,last    你每次读入的都是第一荷载步，结果肯定是一样的  (9x8,f0z  
有没有人具体用ANSYS算过郑院士论文上的实例？结果相同吗？欢迎加入边破加固群：64321843

你这个命令流好像是错的，在李围那本书里有个例子是ANSYS强度折减法的，界面操作，推荐你看看

    

谢谢前辈，刚入行，多指教

是对整体都折减吗，还是只对滑体折减啊???????