[求助]求教“在ANSYS里如何用强度折减法算安全系数？” [复制链接]

各位高手： 0eP ]  
    小弟初学ANSYS，看到不少文章在ANSYS里面用强度折减法求边坡的安全系数，但不知具体过程如何实现，不知哪位可传一个例子给小弟看看，GUI或命令流的都可以啊

同求 呵呵 谁知道啊

可以参考一下徐干成的有关论文，上面有相关内容

不断降低c和摩擦角，直到不收敛

我已经自己解决了，哎

拜托，请楼主换张照片！好嘿人哦！

(1)创建物理环境 6KRC_-  
/COM, Structural                             Tpd|+60g  
/TITLE, Slope stability Analysis                 F+SqJSa  
/FILNAM,support,1                             4~K%,K+Du  
(2)建立模型 LG+2?+tE"  
!进入前处理器 0sA+5*mdM  
/PREP7 0g`$Dap  
!定义单元类型 fpa~~E-  
ET,1,PLANE82                               !定义边坡围岩单元 :OFs"bC  
KEYOPT,1,3,2                                !求解类型为平面应变 FTQNS8  
!定义材料属性 mz|p=[lR|  
!边坡围岩1材料属性   !kS/Ei  
MP,EX,1,3E10   k.ttrKy<q/  
MP,PRXY,1,0.25     Q@ Ze+IhK`  
MP,DENS,1,2500                             !定义塑性模型 `oU|U!|  
TB,DP,1                                     !定义数据表1 dLfB){>S  
TBDATA,1,9E5,42.8                           !给数据表输入数据 0NF=7 j  
!边坡围岩2材料属性                         !采用弹性模型 VTwDa*]AhB  
MP,EX,2,3.2E10                              ！定义弹性模量 |JLXgwML  
MP,PRXY,2,0.24                              !定义泊松比 bgYUsc*uR  
MP,DENS,2,2700                             !定义密度 NXCvS0/h  
!取强度折减系数F=1.2时边坡围岩材料属性 %6W%-`  
MP,EX,3,3E10   bs&>QsI?j  
MP,PRXY,3,0.25     8Drz i!}  
MP,DENS,3,2500 Ii&\LJ  
TB,DP,3                                     !定义数据表3 RG.wu6Av  
TBDATA,1,7.5E5,37.7                         !给数据表输入数据 ]Zz.n5c  
!取强度折减系数F=1.4时边坡围岩材料属性 ;
Rljx3!N  
MP,EX,4,3E10   {SkE`u4Sz  
MP,PRXY,4,0.25                              !定义泊松比 x']Fe7nv  
MP,DENS,4,2500 0T7
t.  
TB,DP,4                                    !定义数据表4 Rc vp@  
TBDATA,1,6.4E5,33.5                        !给数据表输入数据 )*XD"-9  
!取强度折减系数F=1.6时边坡围岩材料属性 ni85Ne$  
MP,EX,5,3E10   "WzKJwFr  
MP,PRXY,5,0.25     c;Gf$9?iC  
MP,DENS,5,2500                             ！定义密度 c`@";+|r  
TB,DP,5                                    !定义数据表5 al"=ld(  
TBDATA,1,5.6E5,30                          !给数据表输入数据 f~10 iD  
!取强度折减系数F=1.8时边坡围岩材料属性 [jv+Of IZ  
MP,EX,6,3E10                               ！定义弹性模量 )|=4H
>?%  
MP,PRXY,6,0.25     I.[Lv7U-  
MP,DENS,6,2500 Fs3 :NH  
TB,DP,6                                    !定义数据表6 w>o/)TTJL  
TBDATA,1,5.0E5,27.2 G*f\ /  
!取强度折减系数F=2.0时边坡围岩材料属性 h343$,))u  
MP,EX,7,3E10   2FcNzAaV  
MP,PRXY,7,0.25                             !定义泊松比 w{*PZb4  
MP,DENS,7,2500 `&9iC 4P  
TB,DP,7                                    !定义数据表7 E&N~h|CL  
TBDATA,1,4.5E5,24.8 3$_JNF`  
!取强度折减系数F=2.2时边坡围岩材料属性 p,.6sk  
MP,EX,8,3E10   N%F4ug@i   
MP,PRXY,8,0.25     suS[P?4  
MP,DENS,8,2500 2){O&8
A  
TB,DP,8                                     !定义数据表8 ob;O,&e0>  
TBDATA,1,4.09E5,22.8 oOHY+'V  
!取强度折减系数F=2.4时边坡围岩材料属性 7`f%?xVn0  
MP,EX,9,3E10   Q5b9q$L$  
MP,PRXY,9,0.25     e%lxRN"b  
MP,DENS,9,2500 ;0U*N& f  
TB,DP,9                                      !定义数据表9 /i~x.i3  
TBDATA,1,3.6E5,21.1                          !给数据表输入数据 zI0d  
!取强度折减系数F=2.6时边坡围岩材料属性 c'>_JlG~  
MP,EX,10,3E10                                ！定义弹性模量   f`)*bx  
MP,PRXY,10,0.25                              !定义泊松比 BwkY;Ur/AL  
MP,DENS,10,2500                           ！定义密度 Wu"1M^a  
TB,DP,10                                   !定义数据表10 g4u6#.m(  
TBDATA,1,3.46E5,19.6                        !给数据表10输入数据 >=4('  
!取强度折减系数F=2.8时边坡围岩材料属性
W7. +  
MP,EX,11,3E10                              ！定义弹性模量 la}cGZ; p.  
MP,PRXY,11,0.25                             !定义泊松比 fy$?~Ji&  
MP,DENS,11,2500 ?N(<w?Gat  
TB,DP,11                                   !定义数据表11 G c\^Kg^#  
TBDATA,1,3.2E5,18.2                         !给数据表11输入数据 UwxszEHC  
!取强度折减系数F=3.0时边坡围岩材料属性 }<YU4EW  
MP,EX,12,3E10                              ！定义弹性模量 P{x6e/  
MP,PRXY,12,0.25     d N$,AOT  
MP,DENS,12,2500                            ！定义密度 !S%0#d2  
TB,DP,12                                   !定义数据表12 W4,'?o  
TBDATA,1,3.0E5,17.1                        !给数据表12输入数据
('{aOiSH  
!建立几何模型 CBv0fQtL  
!创建边坡线模型线 PXyv);#Q`  
K,1,,,,                                      ！ 创建关键点1 ):[}NDmC  
K,2,-800,,, p|(SR~;6  
K,3,-800,-800,,   _RIlGs\.  
K,4,-800,-1200,, i),bAU!+m  
K,5,1200,-1200,,   ap8q`a{j^  
K,6,1200,-800,,   4l7 Ny\J  
K,7,1200,,, K iEmvC  
K,8,1200,378,, d@p#{ -  
K,9,430,378,, Wb>;L@jB7  
l,1,2                                        !连接关键点1、2生成直线L1 1_b*j-j  
l,2,3 14"+ctq  
l,3,4 +4  h!;i  
l,4,5 i)'tt9f$  
l,5,6 3vKTCHbk9  
l,6,7 IJ~j(.W  
l,7,8 8ok=&Gq4  
l,8,9 g60k R7;\  
l,9,1 l2kGFgc  
l,1,7 P@keg*5@  
l,3,6 |;7mDhj=  
!创建边坡面模型 b8_F2  
al,3,4,5,11                              !由线L3、L4、L5、L11生成一个面A1 ;*$e8y2  
al,1,2,11,6,10                            !由线L1、L2、L1、L6、L10生成一个面A2 Jt[,V*:#  
al,7,8,9,10                              !由线L7、L8、L9、L10生成一个面A3 Y!8FW|  
/PNUM, area, on                         !打开面号开关 yIcTc  
!划分网格生成有限元模型 c6
lCF &  
! 划分边坡围岩2网格 <=V{tl  
!设置网格份数 `KN>0R2k  
lsel,s,,,3,5,2                               !选择线L3和L5 d]Y;rqjue  
lesize,all,,,5                               !把所选择线分为5段 0-[naGz  
lsel,s,,,4,11,7                              !选择线L4和L11 Lg~C:BNF  
lesize,all,,,26                              !把所选择线分为26段 0QT
:@v2R  
mat,2                                    !给边坡围岩2赋予2号材料特性 Fu
zb4Df  
type,1                                   !采用单元类型1 ALy7D*Z]w  
mshkey,0                                 !设定自由网格划分 /`l;u7RD  
mshape,0                                 !设定四边形网格划分 }W'4(V;:  
amesh,1                                  ！划分面积A1 2lO(f+  
! 划分边坡围岩1网格 ^86M94k  
!设置网格份数 zPc"r$'0U  
lsel,s,,,2,6,4                               !选择线L2和L6 h=
0a9vIXF  
lsel,a,,,1                                  !附加选择线L1 P%)r4+at  
lesize,all,,,10                              !把所选择线分为10段 Ix6\5}.c9  
lsel,s,,,8,10,2                              !选择线L8和L10 cFt&Efj  
lesize,all,,,16                               !把所选择线分为16段 XPU>} 4{  
lsel,s,,,7,9,2 P1Z"}Qw  
lesize,all,,,12 /OWwC%tM/  
mat,1                                     !给边坡围岩1赋予1号材料特性 BvsSrse  
type,1 1f<R,>  
mshkey,0                                  !设定自由网格划分 #G.eiqh$a  
mshape,0                                  !设定四边形网格划分 &92/qRh7  
amesh,2                                   !划分面积A2 +]nIr'V  
amesh,3                                   ！划分面积A3 oX8EY l  
allsel mEbI\!}H0  
!保存网格模型 ^?Mp(o  
save,Slope-grid.db @lF?+/=$  
(3)施加约束和荷载 D*ZjoU  
!给边坡模型施加约束 jLVG=rO
n  
!边坡两侧施加X方向约束 yKoZj   
nsel,s,loc,x,-800                             !选择X=-800线上所有节点 a_V\[V{R=  
nsel,a,loc,x,1200                             !选择X=1200线上所有节点 _FYA? d}  
d,all,ux                                     !对所选择节点约束X方向位移 s*[ I"iE  
allsel q~b# ml2QS  
!边坡底部施加约束 ":8\2Qp  
nsel,s,loc,y,-1200                           !选择Y=--1200线上所有节点 2 4+  
d,all,uy                                   !对所选择节点约束X、Y方向位移 ^8;MY5Wbs  
d,all,ux                                    =aZ d>{Y  
!施加重力加速度 X!qK[b@Z  
acel,,9.8 CNefk$/cR  
(4)求解 nj'5iiV`]  
/solu fj;ZGbg-O  
!求解设置 OemY'M?ZQ  
antype,static                                  !设定为静力求解 5, ,~k=  
nsubst,100                                   !设定最大子步数为100 C@6:uiT$  
pred,on                                      !打开时间步长预测器 mLqqo2u  
nropt,full                                     !设定牛顿-拉普森选项 j(Lz& *4  
nlgeom,on                                    !打开大位移效果 t\hnnu`Pq  
lnsrch,on                                     !打开线性搜索 Yu\$Y0 {]  
outres,all,all                                  !输出所有项 fJ[ ^_,O  
cnvtol,f,,0.005,2,0.5                            !力收敛准则设定 -dixiJ
=  
cnvtol,u,,0.05,2,1                              !位移收敛准则设定 U8Zb&6  
!边坡在强度折减系数F=1时求解 @k&6\1/U  
allsel Vf&U`K  
solve                                        !进行求解 Yr/$92(  
save,F1,db                                    !把F=1时求解结果保存 Jgv Mx  
!边坡在强度折减系数F=1.2时求解 88~Nrl=co  
finish n82tZpn  
/solu zPa2fS8  
allsel ~c35Y9-5  
asel,s,area,,2,3,1                               !选择面积A2、A3 12z!{k7N  
MPCHG,3,all                                 ! 把所选择单元材料号改为3 e{@RBYX@+c  
allsel !i)!|9e  
solve                                         !进行求解 9L:wfg}8s  
save,F1.2,db                                   !把F=1.2时求解结果保存                                                                                   xeSch?}  
!边坡在强度折减系数F=1.4时求解 iRnjN  
finish 7h}gIm7e"  
/solu &liON1GLM  
allsel ud r\\5  
asel,s,area,,2,3,1                               !选择面积A2、A3 ]JjS$VMauX  
MPCHG,4,all                                 !把所选择单元材料号改为4 x[)-h/&Fh  
allsel K2
zln_W  
solve                                        !进行求解 Qdq;C,}Ai.  
save,F1.4,db                                  !把F=1.4时求解结果保存                             SK\@w9#&$  
!边坡在强度折减系数F=1.6时求解 .Qx5,)@
9  
finish .DhI3'Jrl  
/solu (vO3vCYeQ  
allsel |y U!d %  
asel,s,area,,2,3,1                               !选择面积A2、A3 Q38+`EhLA  
MPCHG,5,all                                 !把所选择单元材料号改为5 A.vAk''(}+  
allsel bRx
I7 '  
solve EL$DvJ~  
save,F1.6,db                                   !把F=1.6时求解结果保存                                                                                   P!W%KobZ7|  
!边坡在强度折减系数F=1.8时求解 UHZ&7jfl  
finish ; R}>SS'  
/solu ZPmqoR[  
allsel TQT3]h6  
asel,s,area,,2,3,1                                !选择面积A2、A3 Xx{| [2`  
MPCHG,6,all                                  !把所选择单元材料号改为6 VGc*aQYa  
allsel qU!dg  
solve ;[}OZt  
save,F1.8,db                                   !把F=1.8时求解结果保存   rJpr;QKf%  
!边坡在强度折减系数F=2.0时求解 ~i"=:D  
finish j~k,d.17M  
/solu E+lR&~mK=  
allsel ?.SGn[  
asel,s,area,,2,3,1                                !选择面积A2、A3 (Lgea  
MPCHG,7,all                                   !把所选择单元材料号改为7 IQPu%n{0v  
allsel ?Fl}@EA#M  
solve N l@G\_  
save,F2.0,db                                   !把F=2.0时求解结果保存                             -U/)y:k!%  
!边坡在强度折减系数F=2.2时求解 N.
JR($N$  
finish /P*XB%y  
/solu ,2fi`9=\  
allsel [t?tLUg|6  
asel,s,area,,2,3,1                                !选择面积A2、A3 x#EE_i/W  
MPCHG,8,all                                  !把所选择单元材料号改为8 KSPa2>lz?  
allsel .D 4G;=Q  
solve |reA`&<q  
save,F2.2,db                                   !把F=2.2时求解结果保存                                                                                   09-8Xzz  
!边坡在强度折减系数F=2.4时求解 'Gn>~
m  
finish Y1-dpML   
/solu <{kPa_`'  
allsel r >'tE7W9  
asel,s,area,,2,3,1                                 !选择面积A2、A3 L;RE5YrH%6  
MPCHG,9,all                                   !把所选择单元材料号改为9 FMVAXOO  
allsel OJcS%-~  
solve m[}k]PB>  
save,F2.4,db                                   !把F=2.4时求解结果保存   LeL
Ut<4~  
!边坡在强度折减系数F=2.6时求解 LL#REK|lm8  
finish JD$g%hcVZa  
/solu .t_t)'L  
allsel PY5&Fwjc  
asel,s,area,,2,3,1                                 !选择面积A2、A3 &s+l/;3  
MPCHG,10,all                                   !把所选择单元材料号改为10 4=^_VDlpd  
allsel ]o<]A[<  
solve Kz"3ba}KH  
save,F2.6,db                                   !把F=2.6时求解结果保存                             =y][j+WH  
!边坡在强度折减系数F=2.8时求解 /*qRbN  
finish TmG)
;B}  
/solu hik.qK  
allsel KV_/fa~Ry  
asel,s,area,,2,3,1                                !选择面积A2、A3 ddfGR/1X  
MPCHG,11,all                                 !把所选择单元材料号改为11 =wh[D$n$~  
allsel Ab@G^SLX  
solve NfvPE]S  
save,F2.8,db                                   !把F=2.8时求解结果保存                                                                                   :*}Q/]N  
!边坡在强度折减系数F=3.0时求解 i//H5D3  
finish >9{?&#]x  
/solu SY+0~5E  
allsel L+8ar9es  
asel,s,area,,2,3,1                             !选择面积A2、A3 5skN'*oG  
MPCHG,12,all                              !把所选择单元材料号改为12 EB[T 5{  
allsel N(7XILC  
solve }3{eVct#|  
save,F3.0,db                                 !把F=3.0时求解结果保存 k_7m[o  
/post1                                       ！进入后处理 y!=,u  
!边坡在强度折减系数F=1时结果分析 VmqJMU>.  
Resume,'F1','db'                              !读入边坡在强度折减系数F=1时 +l7)7qKx  
set,1,last                                    !读入后一个子步 l(Rn=
?  
pldisp,1                                      !绘制边坡模型变形图 u"HGT=Nl  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 |{N{VK  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 PR@6=[|d  
!边坡在强度折减系数F=1.2时结果分析 "N}t =3i$  
Resume,'F1.2','db'                            !读入边坡在强度折减系数F=1.2时 JY"jj}H]|  
set,1,last                                    !读入后一个子步 ,.<mj !YE  
pldisp,1                                      !绘制边坡模型变形图 2<@2_wSJ  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 f;{Q ~  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 1C
B&z@  
!边坡在强度折减系数F=1.4时结果分析 ?@rd,:'dE  
Resume,'F1.4','db'                            !读入边坡在强度折减系数F=1.4时 i(j/C  
set,1,last                                    !读入后一个子步 J~2SGXH)^?  
pldisp,1                                      !绘制边坡模型变形图 9hA`I tS  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 t1?aw<  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 Z mJ<h&  
!边坡在强度折减系数F=1.6时结果分析 sLr47 NC  
Resume,'F1.6','db'                            !读入边坡在强度折减系数F=1.6时 TF8#I28AD  
set,1,last                                    !读入后一个子步 p3M!H2W  
pldisp,1                                      !绘制边坡模型变形图 j9+4},>>CU  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 WQ9e~D"  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 fQfn7FaW_\  
!边坡在强度折减系数F=1.8时结果分析 VE+H! ob A  
Resume,'F1.8','db'                            !读入边坡在强度折减系数F=1.8时 e$~[\ w  
set,1,last                                    !读入后一个子步 zgwez$  
pldisp,1                                      !绘制边坡模型变形图 T?7u [D[[  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 'S_i6K  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 RCWmdR#}V  
!边坡在强度折减系数F=2.0时结果分析 )pHtsd.eP  
Resume,'F2.0','db'                            !读入边坡在强度折减系数F=2.0时 x"b'Pmw  
set,1,last                                    !读入后一个子步 DG;7+2U  
pldisp,1                                      !绘制边坡模型变形图 P 2WAnm  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 l!tR
<$|  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 296}L
W  
!边坡在强度折减系数F=2.2时结果分析 ["3dr@T9Z  
Resume,'F2.2','db'                            !读入边坡在强度折减系数F=2.2时 &&
&-P\3  
set,1,last                                    !读入后一个子步 A8m06  
pldisp,1                                      !绘制边坡模型变形图 f!'i5I]  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 UY(T>4H+h  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 ,Pa*; o\  
!边坡在强度折减系数F=2.4时结果分析 J'%i?cuV  
Resume,'F2.4','db'                            !读入边坡在强度折减系数F=2.4时 <A,V
/']  
set,1,last                                    !读入后一个子步 m Q9
dF,  
pldisp,1                                      !绘制边坡模型变形图 -Uo11'{  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 {:1j>4m2  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 OH t)z.  
!边坡在强度折减系数F=2.6时结果分析 v #+ECx
 
Resume,'F2.6','db'                            !读入边坡在强度折减系数F=2.6时 K7RAmX  
set,1,last                                    !读入后一个子步 P6v ANL-B  
pldisp,1                                      !绘制边坡模型变形图 ~+Gh{,f  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 ~[%CUc"  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 joe)b  
!边坡在强度折减系数F=2.8时结果分析 ,CqWm9  
Resume,'F2.8','db'                            !读入边坡在强度折减系数F=2.8时 %66="1z0@  
set,1,last                                    !读入后一个子步
Xd(^7~i  
pldisp,1                                      !绘制边坡模型变形图 `dvg5qQ  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 hN3FH#YO  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 yx]9rD1cz  
!边坡在强度折减系数F=3.0时结果分析 P{o)Ir8Tt  
Resume,'F3.0','db'                            !读入边坡在强度折减系数F=3.0时 uBlPwb,V  
set,1,last                                    !读入后一个子步 /*#o1W?wQZ  
pldisp,1                                      !绘制边坡模型变形图 tpJA~!mG3  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 w/6X9d  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图

楼上的是一个试算过程，最后求出的安全系数是多少呢？人工判读？？   

  随着折减系数的不断增大并达到某一数值时, 边坡内广义剪应变自坡底向坡顶贯通, 边坡达到极限状态, 此时的折减系数即为安全系数。

好东西,顶一下!

正在学习，判断不收敛有几种办法，比如计算不收敛，位移逐渐增大后又突然急剧降低，塑性区贯通，此时均可以看作失稳，则前一级的即为安全系数！

就是不断降低土体强度，也就是用把土体强度参数如内摩擦角除以F，得到的新内摩擦角用于ansys分析，计算不收敛就说明安全系数为假定的F值