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

各位朋友： 0]HIc  
    大家好！我通过ansys用有限元强度折减法分析了一个边坡的稳定性，结果不同的工况下得出的塑性区是一样的，我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激！附命令流和塑性区如下：(1)创建物理环境 o%iTYR:x  
/COM, Structural                              ! 指定结构分析 b"V-!.02  
/TITLE, Slope stability Analysis                  ! 定义工作标题 TjdYCk]'  
/FILNAM,support,1                            ! 定义工作文件名 X35hLp8 M  
(2)建立模型 S7*:eo  
!进入前处理器 j.*}W4`Q_  
/PREP7 w4mL/j  
!定义单元类型 bx>D  
ET,1,PLANE82                              !定义边坡围岩单元 .nJGxz+X"  
KEYOPT,1,3,2                                !求解类型为平面应变 fQh!1R  
!定义材料属性 J)_42Z  
!边坡围岩1材料属性  x$*OglaS  
MP,EX,1,3E10  #e[5O|V~  
MP,PRXY,1,0.25    ho. a93  
MP,DENS,1,2500                            !定义塑性模型 7y&`H  
TB,DP,1                                    !定义数据表1 *np|PyLP:  
TBDATA,1,9E5,42.8                          !给数据表输入数据 B'yrXa|P  
!边坡围岩2材料属性                        !采用弹性模型 (A_9;uL^_  
MP,EX,2,3.2E10                              ！定义弹性模量 5q[@N J  
MP,PRXY,2,0.24                              !定义泊松比 P? n`n!qZ  
MP,DENS,2,2700                            !定义密度 c"r( l~fc  
!取强度折减系数F=1.2时边坡围岩材料属性 2vLn#  
MP,EX,3,3E10  +$C5V,H~  
MP,PRXY,3,0.25    ~yV0SpL  
MP,DENS,3,2500 SV >EB;<  
TB,DP,3                                    !定义数据表3 z#b6 aP  
TBDATA,1,7.5E5,37.7                        !给数据表输入数据 .'M]cN~  
!取强度折减系数F=1.4时边坡围岩材料属性 TaH9Nu  
MP,EX,4,3E10  bVym  
MP,PRXY,4,0.25                              !定义泊松比 FrS>.!OFn  
MP,DENS,4,2500 4T"P#)z  
TB,DP,4                                    !定义数据表4 9b>a<Z  
TBDATA,1,6.4E5,33.5                        !给数据表输入数据 L '=3y$"],  
!取强度折减系数F=1.6时边坡围岩材料属性 |O
NOF  
MP,EX,5,3E10  W=drp>Uj  
MP,PRXY,5,0.25    VlH9ap  
MP,DENS,5,2500                            ！定义密度 }j x{Cw  
TB,DP,5                                    !定义数据表5 y)2]:nD`B  
TBDATA,1,5.6E5,30                          !给数据表输入数据 mb/Y  
!取强度折减系数F=1.8时边坡围岩材料属性 ul E\>5O4h  
MP,EX,6,3E10                              ！定义弹性模量  *hba>LZ  
MP,PRXY,6,0.25    9S}PCAA;  
MP,DENS,6,2500 >*&[bW'}?  
TB,DP,6                                    !定义数据表6 !#QD;,SE+  
TBDATA,1,5.0E5,27.2 HDYoM  
!取强度折减系数F=2.0时边坡围岩材料属性 6k-]2,\#  
MP,EX,7,3E10  rz(DZV  
MP,PRXY,7,0.25                            !定义泊松比 YAD9'h]d\  
MP,DENS,7,2500 G8@%)$A  
TB,DP,7                                    !定义数据表7 RAvV[QkT  
TBDATA,1,4.5E5,24.8 ]M.ufbguq  
!取强度折减系数F=2.2时边坡围岩材料属性 hM-qC|!  
MP,EX,8,3E10  Xv9CD  
MP,PRXY,8,0.25    nZ]d[  
MP,DENS,8,2500 sy~mcH:%+  
TB,DP,8                                    !定义数据表8 xJ4T7 )*  
TBDATA,1,4.09E5,22.8 ke0Vy(3t{h  
!取强度折减系数F=2.4时边坡围岩材料属性 /dIiFr"e}G  
MP,EX,9,3E10  ~be&T:7.  
MP,PRXY,9,0.25    &A>J>b  
MP,DENS,9,2500 {"s8X(#_sC  
TB,DP,9                                      !定义数据表9 7J\I%r  
TBDATA,1,3.6E5,21.1                          !给数据表输入数据 I>:'
5V  
!取强度折减系数F=2.6时边坡围岩材料属性 x9a0J1Nb-h  
MP,EX,10,3E10                                ！定义弹性模量  =6>mlI>i  
MP,PRXY,10,0.25                              !定义泊松比 j&F&wRD%r  
MP,DENS,10,2500                          ！定义密度 9,4a?.*4~  
TB,DP,10                                  !定义数据表10 <=p"ck@  
TBDATA,1,3.46E5,19.6                        !给数据表10输入数据 FBe1f1 sm  
!取强度折减系数F=2.8时边坡围岩材料属性 v+Hu=RZE  
MP,EX,11,3E10                              ！定义弹性模量 ]C_6I\Z#=W  
MP,PRXY,11,0.25                            !定义泊松比 18~j>fN  
MP,DENS,11,2500 'O CVUF,  
TB,DP,11                                  !定义数据表11 zk?lNs  
TBDATA,1,3.2E5,18.2                        !给数据表11输入数据 QJxcH$  
!取强度折减系数F=3.0时边坡围岩材料属性 ~*&_zPTN  
MP,EX,12,3E10                              ！定义弹性模量 Mb-C DPT  
MP,PRXY,12,0.25    Gz`Zp "i%0  
MP,DENS,12,2500                            ！定义密度 ~{xY{qL  
TB,DP,12                                  !定义数据表12 mIOx)`
$  
TBDATA,1,3.0E5,17.1                        !给数据表12输入数据 ^#;2 Pd>  
!建立几何模型 >t.I,Zn  
!创建边坡线模型线 @qeI4io-n  
K,1,,,,                                      ！ 创建关键点1 3?C$Tl2G8  
K,2,-800,,, P
)`-cfg  
K,3,-800,-800,,  UJO+7h'  
K,4,-800,-1200,, rLp0VKPe  
K,5,1200,-1200,,  gZ 9<H q  
K,6,1200,-800,,  (UpSi6?\  
K,7,1200,,,
{R(q7ALR  
K,8,1200,400,, kIQMIL0+  
K,9,400,400,, o|*,<5t  
l,1,2                                        !连接关键点1、2生成直线L1 =r`>tWs
  
l,2,3 WFTTBUoH  
l,3,4 'M&`l%dIPf  
l,4,5 5'@}8W3b  
l,5,6 `CW=*uBH  
l,6,7 VF;%Z  
l,7,8 
TJ#<wIiX  
l,8,9 vAX %i(4  
l,9,1 S-My6'ar  
l,1,7 J-=&B5"O>  
l,3,6 pj'Yv  
!创建边坡面模型 S);bcowf_  
al,3,4,5,11                              !由线L3、L4、L5、L11生成一个面A1 (Ys0|I3  
al,1,2,11,6,10                            !由线L1、L2、L1、L6、L10生成一个面A2 V;9.7v  
al,7,8,9,10                              !由线L7、L8、L9、L10生成一个面A3 *PD7H9m  
/PNUM, area, on                        !打开面号开关 (T,ST3{*k  
!划分网格生成有限元模型 :Z
,zWk1|  
! 划分边坡围岩2网格 1_%3cN.  
!设置网格份数 I/t2c=f  
lsel,s,,,3,5,2                              !选择线L3和L5 "H6DiPh.E  
lesize,all,,,5                              !把所选择线分为5段 a&9+<  
lsel,s,,,4,11,7                              !选择线L4和L11 L_Om<LO2  
lesize,all,,,26                              !把所选择线分为26段 5s8k^n"A  
mat,2                                    !给边坡围岩2赋予2号材料特性 [i\K#O +f  
type,1                                  !采用单元类型1 Ymx/N+Jl  
mshkey,0                                !设定自由网格划分 *&!&Y*Jzg  
mshape,0                                !设定四边形网格划分 ?z p$Wz;k  
amesh,1                                  ！划分面积A1 zGA#7W2?0  
! 划分边坡围岩1网格 dQTJC %]O  
!设置网格份数 Iq(B
H^K  
lsel,s,,,2,6,4                              !选择线L2和L6 uPt({H  
lsel,a,,,1                                  !附加选择线L1 p%8y!^g  
lesize,all,,,10                              !把所选择线分为10段 K%)u zP  
lsel,s,,,8,10,2                              !选择线L8和L10 G8H=xr#  
lesize,all,,,16                              !把所选择线分为16段 y3kXfSe  
lsel,s,,,7,9,2 RWK|?FD\<  
lesize,all,,,12 eWE7>kwh  
mat,1                                    !给边坡围岩1赋予1号材料特性 =56O-l7T*w  
type,1 rMUn ~  
mshkey,0                                  !设定自由网格划分 wm_xH_{F  
mshape,0                                  !设定四边形网格划分 ??)IPRv?yF  
amesh,2                                  !划分面积A2 sBt,y_LW  
amesh,3                                  ！划分面积A3 g=Rl4F]  
allsel ?`lIsd  
!保存网格模型 s&TPG0W  
save,Slope-grid.db 274F+X  
(3)施加约束和荷载 [<rV "g  
!给边坡模型施加约束 ch!/k  
!边坡两侧施加X方向约束  YH@p\#Y  
nsel,s,loc,x,-800                            !选择X=-800线上所有节点 &( aw  
nsel,a,loc,x,1200                            !选择X=1200线上所有节点 Gg$4O
8  
d,all,ux                                    !对所选择节点约束X方向位移 y38x^fuYJ~  
allsel }I`"$2   
!边坡底部施加约束 S8 .1%sw  
nsel,s,loc,y,-1200                          !选择Y=--1200线上所有节点 rUmaKh?v|X  
d,all,uy                                  !对所选择节点约束X、Y方向位移 ,Y ./9F  
d,all,ux                                  bW-9YXj%  
!施加重力加速度 tFn_{fCc>  
acel,,9.8 ~nDbWv"  
(4)求解 k>SPtiAs  
/solu N"L@  
!求解设置 ~U}0=lRVS  
antype,static                                  !设定为静力求解 b8@?fC+tm  
nsubst,100                                  !设定最大子步数为100 5bBY[qp  
pred,on                                      !打开时间步长预测器 {Dupk0'(  
nropt,full                                    !设定牛顿-拉普森选项 :"{("!x  
nlgeom,on                                    !打开大位移效果 zsOOx% +  
lnsrch,on                                    !打开线性搜索 h6`v%7H?  
outres,all,all                                  !输出所有项 wkc)2z   
cnvtol,f,,0.005,2,0.5                            !力收敛准则设定 Nz1u:D]  
cnvtol,u,,0.05,2,1                              !位移收敛准则设定 8Mg wXH  
!边坡在强度折减系数F=1时求解 \I"n~h^_  
allsel |$  
solve                                        !进行求解 2D;2QdO  
save,F1,db                                    !把F=1时求解结果保存 Pg4go10|  
!边坡在强度折减系数F=1.2时求解 yI$KBx/]n  
finish @` 5P^H7  
/solu Q!.JV.(  
allsel =e*S h0dK  
asel,s,area,,2,3,1                              !选择面积A2、A3 %1h%#/#[  
MPCHG,3,all                                ! 把所选择单元材料号改为3 q|h#J}\  
allsel J@$KF GUs  
solve                                        !进行求解 ;XagLy  
save,F1.2,db                                  !把F=1.2时求解结果保存                                                                                 
M5c$  
!边坡在强度折减系数F=1.4时求解 '8Wu9 phT  
finish Zmx[:-  
/solu `"Lk@  
allsel Qv1<)&Ft<  
asel,s,area,,2,3,1                              !选择面积A2、A3 ((Uw[8#2`  
MPCHG,4,all                                !把所选择单元材料号改为4 JIL(\d  
allsel  |pgrR7G'  
solve                                        !进行求解 cfSQqH  
save,F1.4,db                                  !把F=1.4时求解结果保存                           
@;;3B  
!边坡在强度折减系数F=1.6时求解 Ndmki 7A  
finish nOTe 3?i>  
/solu ;x,yGb`  
allsel ,)?!p_*@:  
asel,s,area,,2,3,1                              !选择面积A2、A3 5LaF'>1yY  
MPCHG,5,all                                !把所选择单元材料号改为5 ;
r?s7b/>  
allsel SR>(GQ,m0;  
solve LiyR,e  
save,F1.6,db                                  !把F=1.6时求解结果保存                                                                                  tY@+d*u  
!边坡在强度折减系数F=1.8时求解 $&Lw 2 c0  
finish '"C& dia  
/solu !(A<  
allsel kltorlH  
asel,s,area,,2,3,1                                !选择面积A2、A3 [[s^rC<d  
MPCHG,6,all                                  !把所选择单元材料号改为6 FQ_a=v  
allsel T|k_$LH  
solve o#uhPUZ
 
save,F1.8,db                                  !把F=1.8时求解结果保存  <H,E1kGw9  
!边坡在强度折减系数F=2.0时求解 9s4>hw@u  
finish C[R|@9NI  
/solu s7?d_+O  
allsel TdeHs{|  
asel,s,area,,2,3,1                                !选择面积A2、A3 SDu%rr7sQ  
MPCHG,7,all                                  !把所选择单元材料号改为7 z"\<GmvB  
allsel _J&IL!S2  
solve &0x;60b  
save,F2.0,db                                  !把F=2.0时求解结果保存                            d~qZ;uw  
!边坡在强度折减系数F=2.2时求解 \ORE;pG  
finish @Ta0v:Y  
/solu AB+HyZ*//  
allsel B,b^_4XX$  
asel,s,area,,2,3,1                                !选择面积A2、A3 Ve\.7s  
MPCHG,8,all                                  !把所选择单元材料号改为8 %6Y\4Fe  
allsel W4#E&8g%  
solve VXPsYR&  
save,F2.2,db                                  !把F=2.2时求解结果保存                                                                                  sjy/[.4-  
!边坡在强度折减系数F=2.4时求解 =xjtPmZ5X  
finish ,5& Rra/  
/solu BD2Gv)?g  
allsel ^XT;n  
asel,s,area,,2,3,1                                !选择面积A2、A3 &8HJ4Vj2  
MPCHG,9,all                                  !把所选择单元材料号改为9 ZFC&&[%-sG  
allsel B-p ].  
solve '$h0l-mQ  
save,F2.4,db                                  !把F=2.4时求解结果保存  II]-mb  
!边坡在强度折减系数F=2.6时求解 %Y=  
finish /u~L3Cp(  
/solu .A F94OlE/  
allsel Ri/D>[  
asel,s,area,,2,3,1                                !选择面积A2、A3 0`4Fa^o]h  
MPCHG,10,all                                  !把所选择单元材料号改为10 eI
cIl2  
allsel \SooIEl@  
solve 1}q(Pn2  
save,F2.6,db                                  !把F=2.6时求解结果保存                            Iuve~ugO  
!边坡在强度折减系数F=2.8时求解 J9&#);(  
finish 8Xa{.y"  
/solu H2{&da@D5  
allsel :gM_v?sy  
asel,s,area,,2,3,1                                !选择面积A2、A3 +{\b&q_  
MPCHG,11,all                                !把所选择单元材料号改为11 Og1Hg B3v  
allsel Wb#ON|.2  
solve ~p{ fl?  
save,F2.8,db                                  !把F=2.8时求解结果保存                                                                                  u}-)ywX  
!边坡在强度折减系数F=3.0时求解 $r1{Nh  
finish _N"c,P0  
/solu %'MR;hQsd8  
allsel obK6GG?ZE  
asel,s,area,,2,3,1                            !选择面积A2、A3 IAd[_<9D  
MPCHG,12,all                              !把所选择单元材料号改为12 PEMuIYm$  
allsel mr+8[0  
solve +7Qj%x\  
save,F3.0,db                                !把F=3.0时求解结果保存 u>.>hQ  
(5）后处理 Ai~d  
/post1                                      ！进入后处理 DfZ)gqp/Av  
!边坡在强度折减系数F=1时结果分析 Zr|\T7w 3  
Resume,'F1','db'                              !读入边坡在强度折减系数F=1时 G(hzW%P  
set,1,last                                    !读入后一个子步 ]#\/1!W
  
pldisp,1                                      !绘制边坡模型变形图 7#sb},J{  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 eY\!}) 5  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 /mQ9}E4X  
!边坡在强度折减系数F=1.2时结果分析 K 5AArI  
Resume,'F1.2','db'                            !读入边坡在强度折减系数F=1.2时 vPsf{[Kr  
set,1,last                                    !读入后一个子步 m88[(l  
pldisp,1                                      !绘制边坡模型变形图 2 W Wr./q  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 ]g]~!":  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 -^,wQW:o)  
!边坡在强度折减系数F=1.4时结果分析 2+C8w%F8  
Resume,'F1.4','db'                            !读入边坡在强度折减系数F=1.4时 *40Z}1ng  
set,1,last                                    !读入后一个子步 KV|ywcGhT  
pldisp,1                                      !绘制边坡模型变形图 4EFP*7X  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 p><DA fB  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 =UV=F/Af^  
!边坡在强度折减系数F=1.6时结果分析 8O.5ML{  
Resume,'F1.6','db'                            !读入边坡在强度折减系数F=1.6时 TnLblkX  
set,1,last                                    !读入后一个子步 -O6\!Wo=-  
pldisp,1                                      !绘制边坡模型变形图 ?j$*a7[w  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 H8$";T(I  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 v]m#+E  
!边坡在强度折减系数F=1.8时结果分析 1-}$sO c  
Resume,'F1.8','db'                            !读入边坡在强度折减系数F=1.8时 +||[H)qym  
set,1,last                                    !读入后一个子步 Dl_SEf6b  
pldisp,1                                      !绘制边坡模型变形图 oI9-jW  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 Y#f
iJ  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 8_uh2`+Bvb  
!边坡在强度折减系数F=2.0时结果分析 7WN$ rl5/  
Resume,'F2.0','db'                            !读入边坡在强度折减系数F=2.0时 m@y_Wt  
set,1,last                                    !读入后一个子步 p/?TU  
pldisp,1                                      !绘制边坡模型变形图 fVxRK\a\\  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 g&L $5  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 {ud^+I&  
!边坡在强度折减系数F=2.2时结果分析 y:dwx*Q9I  
Resume,'F2.2','db'                            !读入边坡在强度折减系数F=2.2时 l]mn4cn3  
set,1,last                                    !读入后一个子步 @v=A)L  
pldisp,1                                      !绘制边坡模型变形图 etL)T":XV  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 *Bx'g| u  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 N0
&#fXO  
!边坡在强度折减系数F=2.4时结果分析 23gJD8i8  
Resume,'F2.4','db'                            !读入边坡在强度折减系数F=2.4时 }9>X M  
set,1,last                                    !读入后一个子步 .iH#8Z  
pldisp,1                                      !绘制边坡模型变形图 @-BgPDi.Z  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 hionR)R4  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 WyJXT.  
!边坡在强度折减系数F=2.6时结果分析 k?'PCV  
Resume,'F2.6','db'                            !读入边坡在强度折减系数F=2.6时 \:-N<[   
set,1,last                                    !读入后一个子步 d5
tpw$A  
pldisp,1                                      !绘制边坡模型变形图 .E0*lem'hE  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 0:k ~lz  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 ;AE%f.Y  
!边坡在强度折减系数F=2.8时结果分析 :]-? l4(%  
Resume,'F2.8','db'                            !读入边坡在强度折减系数F=2.8时 2s8(r8AI  
set,1,last                                    !读入后一个子步 ]@7]mu:oL  
pldisp,1                                      !绘制边坡模型变形图 &Qq4xn+J  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 nwAx47>{  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 p$&6E\#7  
!边坡在强度折减系数F=3.0时结果分析 P%xz"l i  
Resume,'F3.0','db'                            !读入边坡在强度折减系数F=3.0时 aoS]Qp  
set,1,last                                    !读入后一个子步 o!M*cyq  
pldisp,1                                      !绘制边坡模型变形图 0Z>oiBr4  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 Abf=b<bu  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图

如果你这个能算出正确的结果的话,才真的叫出错,......... md/Z[du:'  
{G%!M+n<  
fE/8;v!
=  
mckrR$>  
\HOOWaapN  
Lay+)S.ta[  
8_tMiIE-pS  
-zVa[&  
     

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

不知道二楼的说什么

高人现身啊

做的不错嘛

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

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

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

    

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

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