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

各位朋友： ^f,4=-  
    大家好！我通过ansys用有限元强度折减法分析了一个边坡的稳定性，结果不同的工况下得出的塑性区是一样的，我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激！附命令流和塑性区如下：(1)创建物理环境 tQ9%rb  
/COM, Structural                              ! 指定结构分析 (C. $w  
/TITLE, Slope stability Analysis                  ! 定义工作标题 bX$z)]KKu  
/FILNAM,support,1                            ! 定义工作文件名 tLOGj?/r  
(2)建立模型 t,9+G<)>H  
!进入前处理器 9F807G\4Qt  
/PREP7 \< .BN;t{  
!定义单元类型 c6[m'cy  
ET,1,PLANE82                              !定义边坡围岩单元 mEJ7e#  
KEYOPT,1,3,2                                !求解类型为平面应变 Z^V;B _  
!定义材料属性 {}$rN@OM$  
!边坡围岩1材料属性  "0!h-bQN  
MP,EX,1,3E10  ATkd#k%S  
MP,PRXY,1,0.25    iW6O9~  
MP,DENS,1,2500                            !定义塑性模型 t+KW=eW  
TB,DP,1                                    !定义数据表1 tKt}]KHV  
TBDATA,1,9E5,42.8                          !给数据表输入数据 *W8n8qG%T  
!边坡围岩2材料属性                        !采用弹性模型 #1%@R<`  
MP,EX,2,3.2E10                              ！定义弹性模量 ^# $IoW  
MP,PRXY,2,0.24                              !定义泊松比 zdwQpB,+^  
MP,DENS,2,2700                            !定义密度 U:*rlA@_.  
!取强度折减系数F=1.2时边坡围岩材料属性 <~dfp  
MP,EX,3,3E10  4+hNP'e  
MP,PRXY,3,0.25    h-PJC/>  
MP,DENS,3,2500 eF%M2:&c;  
TB,DP,3                                    !定义数据表3 7"Xy8]i{z  
TBDATA,1,7.5E5,37.7                        !给数据表输入数据 &MpLm&  
!取强度折减系数F=1.4时边坡围岩材料属性 3g;Y  
MP,EX,4,3E10  c$ /.Xp  
MP,PRXY,4,0.25                              !定义泊松比 c \??kQH  
MP,DENS,4,2500 =wX;OK|U(^  
TB,DP,4                                    !定义数据表4 tCrEcjT-  
TBDATA,1,6.4E5,33.5                        !给数据表输入数据  b*Qd9  
!取强度折减系数F=1.6时边坡围岩材料属性 X}Q4;='C-  
MP,EX,5,3E10  US)i"l7:H*  
MP,PRXY,5,0.25    iOZ9A~Ywy  
MP,DENS,5,2500                            ！定义密度 Hpix:To  
TB,DP,5                                    !定义数据表5 Qp<*or@  
TBDATA,1,5.6E5,30                          !给数据表输入数据 KRxJ2  
!取强度折减系数F=1.8时边坡围岩材料属性 2-'_Nwkl*  
MP,EX,6,3E10                              ！定义弹性模量 ug]2wftlQ  
MP,PRXY,6,0.25    7;p/S#P:  
MP,DENS,6,2500 DPf].i#  
TB,DP,6                                    !定义数据表6 OFr"RGW"  
TBDATA,1,5.0E5,27.2 F qJ`d2E  
!取强度折减系数F=2.0时边坡围岩材料属性 vT#R>0@mi  
MP,EX,7,3E10  2@MpWj4  
MP,PRXY,7,0.25                            !定义泊松比 qVh?%c1.Y  
MP,DENS,7,2500 O.FTToh<  
TB,DP,7                                    !定义数据表7 8d*S9p,/  
TBDATA,1,4.5E5,24.8 +}m`$B}mJ  
!取强度折减系数F=2.2时边坡围岩材料属性 fL|9/sojz  
MP,EX,8,3E10  6H0kY/quL|  
MP,PRXY,8,0.25    C1{Q 4(K%  
MP,DENS,8,2500 oqvu8"  
TB,DP,8                                    !定义数据表8 ky|kg@n{  
TBDATA,1,4.09E5,22.8 sFZdj0tQ4  
!取强度折减系数F=2.4时边坡围岩材料属性 HBp??.r  
MP,EX,9,3E10  x9F*$G  
MP,PRXY,9,0.25    >q;| dn9  
MP,DENS,9,2500 0 @#Jz#?  
TB,DP,9                                      !定义数据表9 ^JxVs 7  
TBDATA,1,3.6E5,21.1                          !给数据表输入数据 <,DMD  
!取强度折减系数F=2.6时边坡围岩材料属性 }N9PV/a  
MP,EX,10,3E10                                ！定义弹性模量  J <z ^C  
MP,PRXY,10,0.25                              !定义泊松比 $|7"9W}m*  
MP,DENS,10,2500                          ！定义密度 $zOV*O2  
TB,DP,10                                  !定义数据表10 06HU6d,  
TBDATA,1,3.46E5,19.6                        !给数据表10输入数据 IsT}T}p,t  
!取强度折减系数F=2.8时边坡围岩材料属性 NCg("n,jx  
MP,EX,11,3E10                              ！定义弹性模量 ;L:UYhDbUx  
MP,PRXY,11,0.25                            !定义泊松比 rWWpP<  
MP,DENS,11,2500 NdJ]\>5oN,  
TB,DP,11                                  !定义数据表11 r2f%E:-0G  
TBDATA,1,3.2E5,18.2                        !给数据表11输入数据 t~qSiHw  
!取强度折减系数F=3.0时边坡围岩材料属性 kDJYEI9j>  
MP,EX,12,3E10                              ！定义弹性模量 d0T 8Cwcb  
MP,PRXY,12,0.25    15_"U+O(/  
MP,DENS,12,2500                            ！定义密度 B[mZQ&Gz`a  
TB,DP,12                                  !定义数据表12 1?)h-aN  
TBDATA,1,3.0E5,17.1                        !给数据表12输入数据 yUcU-pQ  
!建立几何模型 j(I(0Yyh  
!创建边坡线模型线 KOR*y(*8  
K,1,,,,                                      ！ 创建关键点1 X7*ossv  
K,2,-800,,, MX*4d{l  
K,3,-800,-800,,  W/RB|TMT  
K,4,-800,-1200,, DPY+{5
q2  
K,5,1200,-1200,,  >Ch2Ep  
K,6,1200,-800,,  1WaQWZ:=  
K,7,1200,,, \gT({XU?  
K,8,1200,400,, x%)oL:ue  
K,9,400,400,, &iD&C>;pf  
l,1,2                                        !连接关键点1、2生成直线L1 -^_^ByJe  
l,2,3 xX%{i0E  
l,3,4 
~rrl"a>  
l,4,5 .HTRvE`X  
l,5,6 k3LHLJZ#  
l,6,7 %a8e_  
l,7,8 p}~Sgi  
l,8,9 YA+R!t:F{  
l,9,1 , MqoX-+  
l,1,7 "w 4^i!\  
l,3,6 ~Q^.7.-T  
!创建边坡面模型 \*5z0A9)5)  
al,3,4,5,11                              !由线L3、L4、L5、L11生成一个面A1 w[vccARQ  
al,1,2,11,6,10                            !由线L1、L2、L1、L6、L10生成一个面A2 yEB#*}K?  
al,7,8,9,10                              !由线L7、L8、L9、L10生成一个面A3 n2o)K;wW+  
/PNUM, area, on                        !打开面号开关 u=PLjrB~}  
!划分网格生成有限元模型 +bso4 }
rS  
! 划分边坡围岩2网格 "8)%XSb  
!设置网格份数 iqoMQ7%  
lsel,s,,,3,5,2                              !选择线L3和L5 ?w'03lr%  
lesize,all,,,5                              !把所选择线分为5段 18~jUYMV  
lsel,s,,,4,11,7                              !选择线L4和L11 g9Dynm5  
lesize,all,,,26                              !把所选择线分为26段 ?W dY{;&  
mat,2                                    !给边坡围岩2赋予2号材料特性 J10/pS  
type,1                                  !采用单元类型1 7<:o4\q?m  
mshkey,0                                !设定自由网格划分 eF0FQlMe[  
mshape,0                                !设定四边形网格划分 <2O#!bX1  
amesh,1                                  ！划分面积A1 d8Cd4qIXX  
! 划分边坡围岩1网格 Bvj  
!设置网格份数 A]%*ye"NT  
lsel,s,,,2,6,4                              !选择线L2和L6 Ei
2M~/  
lsel,a,,,1                                  !附加选择线L1 |BtFT  
lesize,all,,,10                              !把所选择线分为10段 CE (zt  
lsel,s,,,8,10,2                              !选择线L8和L10 .LM|@OeaD!  
lesize,all,,,16                              !把所选择线分为16段 x*j eCD,  
lsel,s,,,7,9,2 N.|zz)y  
lesize,all,,,12 @G>eCj  
mat,1                                    !给边坡围岩1赋予1号材料特性 -oZw+ge}  
type,1 Fv(FRZ)  
mshkey,0                                  !设定自由网格划分 3:w_49~:~  
mshape,0                                  !设定四边形网格划分 r^E(GmW  
amesh,2                                  !划分面积A2 1#|lt\T  
amesh,3                                  ！划分面积A3 a&kt!%p:  
allsel @}r s6 G  
!保存网格模型 o=zl{tZV  
save,Slope-grid.db r_8;aPL  
(3)施加约束和荷载 _/ 5  
!给边坡模型施加约束 $WRRCB/A6  
!边坡两侧施加X方向约束 i$.!8AV6  
nsel,s,loc,x,-800                            !选择X=-800线上所有节点 6c>:h)?  
nsel,a,loc,x,1200                            !选择X=1200线上所有节点 # dUi['  
d,all,ux                                    !对所选择节点约束X方向位移 );S8`V  
allsel 00-2u~D&  
!边坡底部施加约束 >,Zn~8&Z  
nsel,s,loc,y,-1200                          !选择Y=--1200线上所有节点 Wj)v,v2&  
d,all,uy                                  !对所选择节点约束X、Y方向位移 1JOoICjB  
d,all,ux                                  M)JozD%  
!施加重力加速度 all2?neK  
acel,,9.8 Tqs|2at<t  
(4)求解 benqm ~{\  
/solu U.Z5;E0:  
!求解设置 0Bkc93  
antype,static                                  !设定为静力求解 zII^Ny8D  
nsubst,100                                  !设定最大子步数为100 /,<s9 :  
pred,on                                      !打开时间步长预测器 jW^]N$>  
nropt,full                                    !设定牛顿-拉普森选项 vXM{)
  
nlgeom,on                                    !打开大位移效果 /nh3/[u  
lnsrch,on                                    !打开线性搜索 Rb^G~82d?  
outres,all,all                                  !输出所有项 NTGWI$  
cnvtol,f,,0.005,2,0.5                            !力收敛准则设定 M
4 }))  
cnvtol,u,,0.05,2,1                              !位移收敛准则设定 [T3%Xt'4  
!边坡在强度折减系数F=1时求解 dtG>iJ  
allsel Gy6x.GX  
solve                                        !进行求解 k5}Qx'/l  
save,F1,db                                    !把F=1时求解结果保存 l%u8Lq  
!边坡在强度折减系数F=1.2时求解 lQRtsmZ0  
finish ^mut-@ N9  
/solu DHvZ:)aT}  
allsel ^%\MOjSN  
asel,s,area,,2,3,1                              !选择面积A2、A3 Fl(j,B6Z  
MPCHG,3,all                                ! 把所选择单元材料号改为3 XQOM6$~,  
allsel +T,0,^*  
solve                                        !进行求解 s|[qq7  
save,F1.2,db                                  !把F=1.2时求解结果保存                                                                                  ]Ee$ulJ02  
!边坡在强度折减系数F=1.4时求解 >C5u>@%
9O  
finish &:}WfY!hX  
/solu #g2&x sU  
allsel ssH[\i  
asel,s,area,,2,3,1                              !选择面积A2、A3 z63
y8  
MPCHG,4,all                                !把所选择单元材料号改为4 # xX  
allsel @L 6)RF  
solve                                        !进行求解 S*76V"")  
save,F1.4,db                                  !把F=1.4时求解结果保存                            r0pwKRE~t  
!边坡在强度折减系数F=1.6时求解 L@?3E`4/v  
finish W]>%*n  
/solu (*T$:/zIS  
allsel Nn?$}
g  
asel,s,area,,2,3,1                              !选择面积A2、A3 ZP0D)@8  
MPCHG,5,all                                !把所选择单元材料号改为5 YPKB4p#  
allsel Zv8GrkK  
solve 0q]0+o*%  
save,F1.6,db                                  !把F=1.6时求解结果保存                                                                                  =L?2[a$2;  
!边坡在强度折减系数F=1.8时求解 `Hu2a]e9  
finish >$a;+v  
/solu ^nFP#J)_5  
allsel M[u6+`  
asel,s,area,,2,3,1                                !选择面积A2、A3 00r7trZW^  
MPCHG,6,all                                  !把所选择单元材料号改为6 d*U<Ww^q  
allsel iG=Di)O  
solve {DUtdu[  
save,F1.8,db                                  !把F=1.8时求解结果保存  oDt{;S8|]  
!边坡在强度折减系数F=2.0时求解 mwZ)PySm)  
finish Uaj_,qb(  
/solu Wn?),=WQ{  
allsel lO0}  
asel,s,area,,2,3,1                                !选择面积A2、A3 >V\^oh)t]t  
MPCHG,7,all                                  !把所选择单元材料号改为7 ]9W7]$  
allsel b9v<Jk  
solve ##alzC  
save,F2.0,db                                  !把F=2.0时求解结果保存                            rrphOG  
!边坡在强度折减系数F=2.2时求解 /cvMp#<]  
finish "];@N!dA  
/solu 8NJT:6Q7l  
allsel
O{U j  
asel,s,area,,2,3,1                                !选择面积A2、A3 wAL}c(EHO  
MPCHG,8,all                                  !把所选择单元材料号改为8 -]n%+,3L  
allsel {r5OtYmpR  
solve H
9!q)qlK  
save,F2.2,db                                  !把F=2.2时求解结果保存                                                                                  $
1m}lXk  
!边坡在强度折减系数F=2.4时求解 :s-9@Yl|  
finish V"}Jsr  
/solu RU>qj *e  
allsel U9 mK^  
asel,s,area,,2,3,1                                !选择面积A2、A3 m4'jTC$  
MPCHG,9,all                                  !把所选择单元材料号改为9 ip1gCH/?_+  
allsel ',rK\&lL6  
solve MGre_=Dm_  
save,F2.4,db                                  !把F=2.4时求解结果保存  {:Kr't<XzF  
!边坡在强度折减系数F=2.6时求解 MY}K.^4^  
finish ;&O *KhLH  
/solu }-u%6KZ   
allsel Oi-%6&}J  
asel,s,area,,2,3,1                                !选择面积A2、A3 f=S2O_Ee  
MPCHG,10,all                                  !把所选择单元材料号改为10 _8\B~;0  
allsel 1<*U:W $g  
solve Zp9kxm'  
save,F2.6,db                                  !把F=2.6时求解结果保存                            K5jeazasp  
!边坡在强度折减系数F=2.8时求解 R-wz+j#  
finish ]M'~uTf  
/solu x]3[0K5;  
allsel {2i8]Sp1d/  
asel,s,area,,2,3,1                                !选择面积A2、A3 !aL=R)G&e  
MPCHG,11,all                                !把所选择单元材料号改为11 i` ay9J8N  
allsel n.C5w8f  
solve /?Mr2!3N  
save,F2.8,db                                  !把F=2.8时求解结果保存                                                                                  :)&vf<JL  
!边坡在强度折减系数F=3.0时求解 QBN\wL8g  
finish $ \o)-3  
/solu TNK1E  
allsel TZ!@IBu  
asel,s,area,,2,3,1                            !选择面积A2、A3 Z/|=@gpw  
MPCHG,12,all                              !把所选择单元材料号改为12 >[~7fxjK-  
allsel dep"$pys>  
solve
-*Th=B-  
save,F3.0,db                                !把F=3.0时求解结果保存 <$otBC/%  
(5）后处理 GS}JyU  
/post1                                      ！进入后处理 >Q?
8tGfB  
!边坡在强度折减系数F=1时结果分析 }Gy M<!:  
Resume,'F1','db'                              !读入边坡在强度折减系数F=1时 n1:q:qMR1  
set,1,last                                    !读入后一个子步 Z_Y'#5o#  
pldisp,1                                      !绘制边坡模型变形图 ,H^!G\  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 |{_>H'  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 UP%6s:>:  
!边坡在强度折减系数F=1.2时结果分析 ["4Tn0g ;  
Resume,'F1.2','db'                            !读入边坡在强度折减系数F=1.2时 ^Xu4N"@  
set,1,last                                    !读入后一个子步 ~-B+7  
pldisp,1                                      !绘制边坡模型变形图 (mKH,r  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 Ndgx@LTQQ  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 S kB*w'k  
!边坡在强度折减系数F=1.4时结果分析 ,LxkdV  
Resume,'F1.4','db'                            !读入边坡在强度折减系数F=1.4时 %/5Wj_|p  
set,1,last                                    !读入后一个子步 49c-`[d L  
pldisp,1                                      !绘制边坡模型变形图 ~!cxRd5;F  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 e0otr_)3F  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 EbVva{;#$;  
!边坡在强度折减系数F=1.6时结果分析 .aVtd [  
Resume,'F1.6','db'                            !读入边坡在强度折减系数F=1.6时 W8;!rFW  
set,1,last                                    !读入后一个子步 *c&|2EsZ  
pldisp,1                                      !绘制边坡模型变形图 M]/DKo  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 bSz7?NAp  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 $$"G1<EZ  
!边坡在强度折减系数F=1.8时结果分析 Dg{d^>T!_x  
Resume,'F1.8','db'                            !读入边坡在强度折减系数F=1.8时 c"pOi&  
set,1,last                                    !读入后一个子步 F[`dX  
pldisp,1                                      !绘制边坡模型变形图 Y:, rN  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 R^P>yk8  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 vLyazVj..  
!边坡在强度折减系数F=2.0时结果分析 a7453
s  
Resume,'F2.0','db'                            !读入边坡在强度折减系数F=2.0时 bovAFdHW  
set,1,last                                    !读入后一个子步 BMuEfa^  
pldisp,1                                      !绘制边坡模型变形图 cCNRv$IO\  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 !\9^|Ef?  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 Kfr1k  
!边坡在强度折减系数F=2.2时结果分析 o0|Ex\  
Resume,'F2.2','db'                            !读入边坡在强度折减系数F=2.2时 2CneR
KQy  
set,1,last                                    !读入后一个子步 eF9GhwE=  
pldisp,1                                      !绘制边坡模型变形图 { I\og  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 evYn}  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 7]ysvSM  
!边坡在强度折减系数F=2.4时结果分析 SgehOu  
Resume,'F2.4','db'                            !读入边坡在强度折减系数F=2.4时 Q~VM.G  
set,1,last                                    !读入后一个子步 K7K/P{@9[9  
pldisp,1                                      !绘制边坡模型变形图 Hbd>sS  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 ( N};.DB1Y  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 fb>$p_s]  
!边坡在强度折减系数F=2.6时结果分析 Jww#zEK  
Resume,'F2.6','db'                            !读入边坡在强度折减系数F=2.6时 L/`1K_\l  
set,1,last                                    !读入后一个子步 ahy6a,)K~  
pldisp,1                                      !绘制边坡模型变形图 a#;;0R $  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 -)VjjKz]8  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 n-yUt72  
!边坡在强度折减系数F=2.8时结果分析 lw.[qP  
Resume,'F2.8','db'                            !读入边坡在强度折减系数F=2.8时 K,,) FM  
set,1,last                                    !读入后一个子步 *kg->J  
pldisp,1                                      !绘制边坡模型变形图 PPIO<K 3`  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 Jou*e%  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 MGt>:&s(]  
!边坡在强度折减系数F=3.0时结果分析 6=;(~k&x9:  
Resume,'F3.0','db'                            !读入边坡在强度折减系数F=3.0时 +>K&zS  
set,1,last                                    !读入后一个子步 6w,"i#E!  
pldisp,1                                      !绘制边坡模型变形图 *4}NLUVX
 
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 ZDgT"53  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图

如果你这个能算出正确的结果的话,才真的叫出错,.........  YYYF a  
uH%b rbrU  
<KU0K  
eTvjo(Lvx  
vjx'yh|  
`y.4FA4"8  
}%/mPbd#  
z}u`45W+  
     

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

不知道二楼的说什么

高人现身啊

做的不错嘛

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

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

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

    

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

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