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

各位朋友： #kC~qux^  
    大家好！我通过ansys用有限元强度折减法分析了一个边坡的稳定性，结果不同的工况下得出的塑性区是一样的，我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激！附命令流和塑性区如下：(1)创建物理环境 Q VTL}AT2:  
/COM, Structural                              ! 指定结构分析 ;_cTrjMv\  
/TITLE, Slope stability Analysis                  ! 定义工作标题 [inlxJD  
/FILNAM,support,1                            ! 定义工作文件名 >-MnB  
(2)建立模型 N!K%aH~O  
!进入前处理器
Jp=qPG|  
/PREP7 TWRP|i!i  
!定义单元类型 z0 9Gp}^;  
ET,1,PLANE82                              !定义边坡围岩单元 P~\rP6 ;  
KEYOPT,1,3,2                                !求解类型为平面应变 60+zoL'  
!定义材料属性 I0}.!  
!边坡围岩1材料属性  ztO)~uL  
MP,EX,1,3E10  +KTfGwKt  
MP,PRXY,1,0.25    7%^G]AFi  
MP,DENS,1,2500                            !定义塑性模型 /I7V\  
TB,DP,1                                    !定义数据表1 ='m$O  
TBDATA,1,9E5,42.8                          !给数据表输入数据 /z-rBfdy^  
!边坡围岩2材料属性                        !采用弹性模型 k)b{UFRW  
MP,EX,2,3.2E10                              ！定义弹性模量 ]\M{Abqd{  
MP,PRXY,2,0.24                              !定义泊松比 x4q}xwH  
MP,DENS,2,2700                            !定义密度 # ?2*I2_  
!取强度折减系数F=1.2时边坡围岩材料属性 s>>&3jfM  
MP,EX,3,3E10  roS" q~GS,  
MP,PRXY,3,0.25    c]9gf\WW  
MP,DENS,3,2500 mo| D  
TB,DP,3                                    !定义数据表3 2)=whnFS  
TBDATA,1,7.5E5,37.7                        !给数据表输入数据 ,t`V^(PEq  
!取强度折减系数F=1.4时边坡围岩材料属性 6F(z6_<  
MP,EX,4,3E10  ^N={4'G)  
MP,PRXY,4,0.25                              !定义泊松比 =!\Nh,\eQ  
MP,DENS,4,2500 xTAfVN  
TB,DP,4                                    !定义数据表4 %%NoXW  
TBDATA,1,6.4E5,33.5                        !给数据表输入数据 );0  
!取强度折减系数F=1.6时边坡围岩材料属性 Dd?G4xUG  
MP,EX,5,3E10  u nv:sV#b  
MP,PRXY,5,0.25    JQM_9
6\  
MP,DENS,5,2500                            ！定义密度 Usf@kVQ  
TB,DP,5                                    !定义数据表5
{"wF;*U.V  
TBDATA,1,5.6E5,30                          !给数据表输入数据 ZG=]b%  
!取强度折减系数F=1.8时边坡围岩材料属性 UdO8KD#r3  
MP,EX,6,3E10                              ！定义弹性模量 hV7EjQp  
MP,PRXY,6,0.25    ,j%\3g`  
MP,DENS,6,2500 lM\dK)p21O  
TB,DP,6                                    !定义数据表6 IO\1nB$0nb  
TBDATA,1,5.0E5,27.2 KTm^}')C8  
!取强度折减系数F=2.0时边坡围岩材料属性 ^OV; P[  
MP,EX,7,3E10  /]U$OP*0  
MP,PRXY,7,0.25                            !定义泊松比  |#yu  
MP,DENS,7,2500 if'=W6W  
TB,DP,7                                    !定义数据表7 m^ xTV-#l@  
TBDATA,1,4.5E5,24.8 hY4#4A`I  
!取强度折减系数F=2.2时边坡围岩材料属性 #&|"t<}  
MP,EX,8,3E10  v<1@"9EH  
MP,PRXY,8,0.25    iV{_?f1jo  
MP,DENS,8,2500 oywiX@]~7  
TB,DP,8                                    !定义数据表8 [piK"N  
TBDATA,1,4.09E5,22.8 fV"Y/9}(  
!取强度折减系数F=2.4时边坡围岩材料属性 I1]YT  
MP,EX,9,3E10  t1Ts!Q2  
MP,PRXY,9,0.25    d'_q9uf'  
MP,DENS,9,2500 f+:iz'b#U  
TB,DP,9                                      !定义数据表9 0C<\m\|~k  
TBDATA,1,3.6E5,21.1                          !给数据表输入数据 [(n5-#1S  
!取强度折减系数F=2.6时边坡围岩材料属性 JO|j?%6YY  
MP,EX,10,3E10                                ！定义弹性模量  6(E4l5%  
MP,PRXY,10,0.25                              !定义泊松比 K
&[0`sH!  
MP,DENS,10,2500                          ！定义密度 )la3GT*1mS  
TB,DP,10                                  !定义数据表10 +-!3ruwSn  
TBDATA,1,3.46E5,19.6                        !给数据表10输入数据 q-z1ElrN7u  
!取强度折减系数F=2.8时边坡围岩材料属性 &y_t,8>5  
MP,EX,11,3E10                              ！定义弹性模量 Wet0qt]  
MP,PRXY,11,0.25                            !定义泊松比 ;#Po}8Y=  
MP,DENS,11,2500 ?T/4 =  
TB,DP,11                                  !定义数据表11 WM+8<|)n  
TBDATA,1,3.2E5,18.2                        !给数据表11输入数据 {7e(0QK  
!取强度折减系数F=3.0时边坡围岩材料属性 Q`bXsH  
MP,EX,12,3E10                              ！定义弹性模量 /O[6PG  
MP,PRXY,12,0.25    :I#.d7`uk  
MP,DENS,12,2500                            ！定义密度 08ZvRy(Je<  
TB,DP,12                                  !定义数据表12 V[.{cY?6  
TBDATA,1,3.0E5,17.1                        !给数据表12输入数据 NO*,}aeG  
!建立几何模型 u$JAjA
 
!创建边坡线模型线 "(mJupI  
K,1,,,,                                      ！ 创建关键点1 ;2kQ)Bq"  
K,2,-800,,, kQ=bd{a6  
K,3,-800,-800,,  12^uu)6Xm,  
K,4,-800,-1200,, 0/QDfA?  
K,5,1200,-1200,,  >v,X:B?+FL  
K,6,1200,-800,,  
g]f<k2  
K,7,1200,,, Szzj9K  
K,8,1200,400,, [4yHXZxza  
K,9,400,400,, B
e{@ L  
l,1,2                                        !连接关键点1、2生成直线L1 ' #K@%P  
l,2,3 XoyxS:=>|[  
l,3,4 |lLe^FM  
l,4,5 g=td*S  
l,5,6 xC< )]  
l,6,7 O8Mypv/C  
l,7,8  m}yu4  
l,8,9 Qy:yz  
l,9,1 _h1bVd-  
l,1,7 9'(_*KSH  
l,3,6 X0j\nXk  
!创建边坡面模型 8ve-g\C8 H  
al,3,4,5,11                              !由线L3、L4、L5、L11生成一个面A1 /o4_rzR?  
al,1,2,11,6,10                            !由线L1、L2、L1、L6、L10生成一个面A2 j"jssbu}  
al,7,8,9,10                              !由线L7、L8、L9、L10生成一个面A3 8J,^O04<  
/PNUM, area, on                        !打开面号开关 ^$oa`B^2JM  
!划分网格生成有限元模型 k)knyEUi  
! 划分边坡围岩2网格
r8(oTx  
!设置网格份数 Jz''UJY/O  
lsel,s,,,3,5,2                              !选择线L3和L5 O+-+=W  
lesize,all,,,5                              !把所选择线分为5段 w^L`"  
lsel,s,,,4,11,7                              !选择线L4和L11 ,i*rHMe  
lesize,all,,,26                              !把所选择线分为26段 `)O9 '568  
mat,2                                    !给边坡围岩2赋予2号材料特性 `6rLd>=R  
type,1                                  !采用单元类型1 wQ(DX!   
mshkey,0                                !设定自由网格划分 z"6o|]9I  
mshape,0                                !设定四边形网格划分 \0|x<~#j'  
amesh,1                                  ！划分面积A1 }e7/F[c.U  
! 划分边坡围岩1网格 Q^{TcL8  
!设置网格份数 .EhC\QpP  
lsel,s,,,2,6,4                              !选择线L2和L6 f?Ex$gnI  
lsel,a,,,1                                  !附加选择线L1 2@(+l*.Q  
lesize,all,,,10                              !把所选择线分为10段 9?Bh8%$  
lsel,s,,,8,10,2                              !选择线L8和L10 hEjvtfM9\-  
lesize,all,,,16                              !把所选择线分为16段 \WE/#To  
lsel,s,,,7,9,2 UusAsezm:  
lesize,all,,,12 Z(:\Vj"  
mat,1                                    !给边坡围岩1赋予1号材料特性 jpi,BVTI-X  
type,1 5JOfJ$(n  
mshkey,0                                  !设定自由网格划分 :/6:&7s  
mshape,0                                  !设定四边形网格划分 bN?*p($/  
amesh,2                                  !划分面积A2 E`N`  
amesh,3                                  ！划分面积A3 *effDNE!  
allsel igOX0  
!保存网格模型 0^{Tq0Ri[  
save,Slope-grid.db !
o| ex+z;  
(3)施加约束和荷载 f.ua,,P.  
!给边坡模型施加约束 G$zY&  
!边坡两侧施加X方向约束 x6>WvFZ  
nsel,s,loc,x,-800                            !选择X=-800线上所有节点 <2*+Y|Lk2  
nsel,a,loc,x,1200                            !选择X=1200线上所有节点 G,A?yM'Vw  
d,all,ux                                    !对所选择节点约束X方向位移 ,pcyU\68v  
allsel M]V j  
!边坡底部施加约束 pYCMJK-H  
nsel,s,loc,y,-1200                          !选择Y=--1200线上所有节点 O#tmB?n*  
d,all,uy                                  !对所选择节点约束X、Y方向位移 y2%[/L:u~  
d,all,ux                                  i.9}bw 9u@  
!施加重力加速度 bIT[\Q  
acel,,9.8 fe4/[S{a   
(4)求解 OY"BaSEOw}  
/solu 1mgw0Q
O  
!求解设置 {{A=^rr%C  
antype,static                                  !设定为静力求解 `
mkOjsj &  
nsubst,100                                  !设定最大子步数为100 '!X`X=  
pred,on                                      !打开时间步长预测器 qw4wg9w5p  
nropt,full                                    !设定牛顿-拉普森选项 wB8548C}-  
nlgeom,on                                    !打开大位移效果 {(-TWh7V  
lnsrch,on                                    !打开线性搜索 K{y`Sb~k  
outres,all,all                                  !输出所有项 CEk[&39"  
cnvtol,f,,0.005,2,0.5                            !力收敛准则设定 Iv7BIK^0  
cnvtol,u,,0.05,2,1                              !位移收敛准则设定 |x.[*'X@  
!边坡在强度折减系数F=1时求解 J{Ij  
allsel XPYf1
H  
solve                                        !进行求解 H|s Iw:  
save,F1,db                                    !把F=1时求解结果保存 W*H%\Y:N  
!边坡在强度折减系数F=1.2时求解 j.Y!E<e4]  
finish d;%~\+)x4  
/solu (|W6p%(  
allsel 2aivc,m{r  
asel,s,area,,2,3,1                              !选择面积A2、A3 pC4uar  
MPCHG,3,all                                ! 把所选择单元材料号改为3 ]mBlXE:Z  
allsel 2P57C;N8|  
solve                                        !进行求解 L%7WHtU*#  
save,F1.2,db                                  !把F=1.2时求解结果保存                                                                                  +L?;g pVE&  
!边坡在强度折减系数F=1.4时求解 k;umLyz  
finish g3n>}\xG>  
/solu 6mZpyt  
allsel O,@QGUoA  
asel,s,area,,2,3,1                              !选择面积A2、A3 3rh t5n2-  
MPCHG,4,all                                !把所选择单元材料号改为4 ,vi6<C\  
allsel L#vk77  
solve                                        !进行求解 W[!bF'-10  
save,F1.4,db                                  !把F=1.4时求解结果保存                            -}qay@cDt  
!边坡在强度折减系数F=1.6时求解 ;).QhHeg>  
finish `5t~ Vlp  
/solu 1%.CtTi  
allsel .Xta;Py|J  
asel,s,area,,2,3,1                              !选择面积A2、A3 ld~*w  
MPCHG,5,all                                !把所选择单元材料号改为5 N}bZdE9F  
allsel g[2[ zIB=  
solve w2"]%WS%  
save,F1.6,db                                  !把F=1.6时求解结果保存                                                                                  A}!D&s&UH  
!边坡在强度折减系数F=1.8时求解 1g!%ej jd  
finish 1\f8-:C  
/solu AxJf\B8  
allsel 0} \;R5a<  
asel,s,area,,2,3,1                                !选择面积A2、A3 <Dnv=)Rq  
MPCHG,6,all                                  !把所选择单元材料号改为6 blV'-Al  
allsel bxz6 >>  
solve 7Il /+l(  
save,F1.8,db                                  !把F=1.8时求解结果保存  ng2yZ @$  
!边坡在强度折减系数F=2.0时求解 %'F[(VB  
finish [:Odb?+`F  
/solu >48)@sS  
allsel x@@k_'~t%  
asel,s,area,,2,3,1                                !选择面积A2、A3 wW/wvC-  
MPCHG,7,all                                  !把所选择单元材料号改为7 NLWj5K)1P  
allsel 'vIVsv<p  
solve mnMY)-6C  
save,F2.0,db                                  !把F=2.0时求解结果保存                            i#lO{ ]  
!边坡在强度折减系数F=2.2时求解 =>C3IR/  
finish ~xZ)b
tf  
/solu ?IG+U TI  
allsel xI=[=;L  
asel,s,area,,2,3,1                                !选择面积A2、A3 !]f:dWSLB  
MPCHG,8,all                                  !把所选择单元材料号改为8 kZ_5R#xK  
allsel ~o;*{ Q  
solve j|% C?N  
save,F2.2,db                                  !把F=2.2时求解结果保存                                                                                  G_p13{"IM  
!边坡在强度折减系数F=2.4时求解 e3&.RrA  
finish ZONe
}tv:  
/solu n]JfdI  
allsel D/zp_9B  
asel,s,area,,2,3,1                                !选择面积A2、A3 QEL3b4Vm  
MPCHG,9,all                                  !把所选择单元材料号改为9 wmv/?g  
allsel M\9p-%"L  
solve 0Ioa;XgOn  
save,F2.4,db                                  !把F=2.4时求解结果保存  $uNYus^vS  
!边坡在强度折减系数F=2.6时求解 W$Z""
 
finish g|3FJA/  
/solu V,9UOC,Gn  
allsel DOo34l6#  
asel,s,area,,2,3,1                                !选择面积A2、A3 F[|aDj@q e  
MPCHG,10,all                                  !把所选择单元材料号改为10 |w^nCsv  
allsel l< |)LDq~  
solve r+l3J>:K  
save,F2.6,db                                  !把F=2.6时求解结果保存                            0Z[8d0  
!边坡在强度折减系数F=2.8时求解 == i?lbj  
finish dJg72?"ka  
/solu Z"<tEOs/En  
allsel oY%NDTVN  
asel,s,area,,2,3,1                                !选择面积A2、A3 s2+s1%^Ll  
MPCHG,11,all                                !把所选择单元材料号改为11 qxwD4L`S  
allsel :] +D+[c)  
solve k!,&L$sG  
save,F2.8,db                                  !把F=2.8时求解结果保存                                                                                  blB00   
!边坡在强度折减系数F=3.0时求解 n47v5.Wn  
finish b{d@:"  
/solu t?kb
N\,  
allsel ;,]Wtmu)7  
asel,s,area,,2,3,1                            !选择面积A2、A3 ~); 7D'[  
MPCHG,12,all                              !把所选择单元材料号改为12 ;i&'va$  
allsel 3
>sA_  
solve Qjh @oWT  
save,F3.0,db                                !把F=3.0时求解结果保存 |4FvPR[  
(5）后处理 hbdM}"&]  
/post1                                      ！进入后处理 ZgI1Byf  
!边坡在强度折减系数F=1时结果分析 j1,ir  
Resume,'F1','db'                              !读入边坡在强度折减系数F=1时  bR5+({yH  
set,1,last                                    !读入后一个子步 D7x"P-ie  
pldisp,1                                      !绘制边坡模型变形图 M>g\Y  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 t7DT5SrR  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 "(d7:!%  
!边坡在强度折减系数F=1.2时结果分析 Go_~8w0<  
Resume,'F1.2','db'                            !读入边坡在强度折减系数F=1.2时 djcCm5m  
set,1,last                                    !读入后一个子步 oW/ #/;|`  
pldisp,1                                      !绘制边坡模型变形图 J4ltHk.|  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 jkIgEF2d*  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 +lqX;*a=N  
!边坡在强度折减系数F=1.4时结果分析 ;/Dp  
Resume,'F1.4','db'                            !读入边坡在强度折减系数F=1.4时 @(A[H^E  
set,1,last                                    !读入后一个子步 Dos`lh  
pldisp,1                                      !绘制边坡模型变形图 "o[j'  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 HI30-$9  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 Nu'T0LPNq(  
!边坡在强度折减系数F=1.6时结果分析 ;HeUD5Nt6F  
Resume,'F1.6','db'                            !读入边坡在强度折减系数F=1.6时 RI_3X5.KQ  
set,1,last                                    !读入后一个子步 WY%'ps_]<  
pldisp,1                                      !绘制边坡模型变形图 'e>0*hF[  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 9rmOf Jo:  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 oUBn:Ir@  
!边坡在强度折减系数F=1.8时结果分析 vN%SN>=L<  
Resume,'F1.8','db'                            !读入边坡在强度折减系数F=1.8时 ceR zHq=  
set,1,last                                    !读入后一个子步 +H~})PeQ  
pldisp,1                                      !绘制边坡模型变形图 3Ga!)  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 y\&`A:^[ A  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 @nS+!t{  
!边坡在强度折减系数F=2.0时结果分析 V}kZowWD  
Resume,'F2.0','db'                            !读入边坡在强度折减系数F=2.0时 G? "6[w/p  
set,1,last                                    !读入后一个子步 5l"v:Px  
pldisp,1                                      !绘制边坡模型变形图 /_P5UE(  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 dw!cDfT+  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 rBZ0(XSZQ  
!边坡在强度折减系数F=2.2时结果分析 i7w>Nvj]  
Resume,'F2.2','db'                            !读入边坡在强度折减系数F=2.2时 E(oI0*S.5  
set,1,last                                    !读入后一个子步 qJN2\e2~f  
pldisp,1                                      !绘制边坡模型变形图 /r Hd9^Y  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 3R[5prE<  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 O?9&6x   
!边坡在强度折减系数F=2.4时结果分析 1^zpO~@S  
Resume,'F2.4','db'                            !读入边坡在强度折减系数F=2.4时 AVA hS}*t  
set,1,last                                    !读入后一个子步 \]W*0t>s  
pldisp,1                                      !绘制边坡模型变形图
f6ad@2  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图
y-93 >Y  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 >I3#ALF  
!边坡在强度折减系数F=2.6时结果分析 <OcD[5  
Resume,'F2.6','db'                            !读入边坡在强度折减系数F=2.6时 S]3t{s#JW7  
set,1,last                                    !读入后一个子步 y#Ao6Od6  
pldisp,1                                      !绘制边坡模型变形图 iii2nmiK  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 q(J3fjY)  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 nDSmr  
!边坡在强度折减系数F=2.8时结果分析 COa"zg  
Resume,'F2.8','db'                            !读入边坡在强度折减系数F=2.8时 _kb $S  
set,1,last                                    !读入后一个子步 .n
s1;8  
pldisp,1                                      !绘制边坡模型变形图 >c>f6  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 Nj_h+=UE!  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 T^ ^o  
!边坡在强度折减系数F=3.0时结果分析 S&% GB  
Resume,'F3.0','db'                            !读入边坡在强度折减系数F=3.0时 %klC& _g~_  
set,1,last                                    !读入后一个子步 nTweQ  
pldisp,1                                      !绘制边坡模型变形图 &JM|u ww?1  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 *;wPAQE  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图

如果你这个能算出正确的结果的话,才真的叫出错,......... 0 V
:z(r  
PK@hf[YHe  
s88lN=;  
x8xSA*@k  
ML!Zm[I9  
X|)Ox ,(  
8S[`(] )  
|T~C($9  
     

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

不知道二楼的说什么

高人现身啊

做的不错嘛

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

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

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

    

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

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