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

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

同求 呵呵 谁知道啊

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

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

我已经自己解决了，哎

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

(1)创建物理环境 Ud$Q0m&  
/COM, Structural                             DG3[^B  
/TITLE, Slope stability Analysis                 I"B8_  
/FILNAM,support,1                             [w,(EE   
(2)建立模型 FX7M4t#<  
!进入前处理器 a|U}Ammr  
/PREP7 n;%y  
!定义单元类型 8Bj4_!g  
ET,1,PLANE82                               !定义边坡围岩单元 -<xyC8$^$  
KEYOPT,1,3,2                                !求解类型为平面应变 YbzM6u2  
!定义材料属性 R_7 6W&  
!边坡围岩1材料属性   dl:-k r8  
MP,EX,1,3E10   mU50pM~/i  
MP,PRXY,1,0.25     +a #lofhv  
MP,DENS,1,2500                             !定义塑性模型 q1STRYb  
TB,DP,1                                     !定义数据表1 =PV/`I_h  
TBDATA,1,9E5,42.8                           !给数据表输入数据 4ffU;6~l'  
!边坡围岩2材料属性                         !采用弹性模型 \WcB9  
MP,EX,2,3.2E10                              ！定义弹性模量 ]\7lbLv  
MP,PRXY,2,0.24                              !定义泊松比 5qC:yI  
MP,DENS,2,2700                             !定义密度 2%5?Fn=  
!取强度折减系数F=1.2时边坡围岩材料属性 `N$!s7M  
MP,EX,3,3E10   XjxPIdX_H  
MP,PRXY,3,0.25     %06vgjOa (  
MP,DENS,3,2500 AfN&n= d K  
TB,DP,3                                     !定义数据表3 UkZ\cc}aC/  
TBDATA,1,7.5E5,37.7                         !给数据表输入数据 >'N!dM.+9  
!取强度折减系数F=1.4时边坡围岩材料属性 8[oYZrg  
MP,EX,4,3E10   #Gd7M3  
MP,PRXY,4,0.25                              !定义泊松比 )AJ=an||5  
MP,DENS,4,2500 \9'!"-i  
TB,DP,4                                    !定义数据表4 7^Na9]PY  
TBDATA,1,6.4E5,33.5                        !给数据表输入数据 -/zp&*0gcx  
!取强度折减系数F=1.6时边坡围岩材料属性 R+d< fe  
MP,EX,5,3E10   _AprkI_  
MP,PRXY,5,0.25     BMdSf(l  
MP,DENS,5,2500                             ！定义密度 t=Jm|wJnUA  
TB,DP,5                                    !定义数据表5 l~rb]6E  
TBDATA,1,5.6E5,30                          !给数据表输入数据 Qo.Uqz.C  
!取强度折减系数F=1.8时边坡围岩材料属性 Vg{Zv4+t  
MP,EX,6,3E10                               ！定义弹性模量 -3` "E%9  
MP,PRXY,6,0.25     G)8ChnJa!m  
MP,DENS,6,2500 n08; <  
TB,DP,6                                    !定义数据表6 []"=]f{1};  
TBDATA,1,5.0E5,27.2 , |l@j%  
!取强度折减系数F=2.0时边坡围岩材料属性 U-eI\Lu  
MP,EX,7,3E10   qc
YNtEs*c  
MP,PRXY,7,0.25                             !定义泊松比 [A]Ca$':  
MP,DENS,7,2500 mpAHL(  
TB,DP,7                                    !定义数据表7 2 Kla8  
TBDATA,1,4.5E5,24.8 J,(7.+`~#  
!取强度折减系数F=2.2时边坡围岩材料属性 S.W^7Ap  
MP,EX,8,3E10   >ukQ, CE~  
MP,PRXY,8,0.25     ;gUXvx~~r  
MP,DENS,8,2500 'PpZ/ry$  
TB,DP,8                                     !定义数据表8 t7bqk!6hM\  
TBDATA,1,4.09E5,22.8 IM=+3W;ak  
!取强度折减系数F=2.4时边坡围岩材料属性 =e
{.yggE  
MP,EX,9,3E10   >fH*XP>(  
MP,PRXY,9,0.25     5}XvL'  
MP,DENS,9,2500 q? 9GrwL8F  
TB,DP,9                                      !定义数据表9 tFiR!f)  
TBDATA,1,3.6E5,21.1                          !给数据表输入数据 Ig9d#c  
!取强度折减系数F=2.6时边坡围岩材料属性 84^[/d;!  
MP,EX,10,3E10                                ！定义弹性模量   Tm_8<$ 7  
MP,PRXY,10,0.25                              !定义泊松比 v,jhE9_O0  
MP,DENS,10,2500                           ！定义密度 f?^S bp
 
TB,DP,10                                   !定义数据表10 6{.J:S9n   
TBDATA,1,3.46E5,19.6                        !给数据表10输入数据 hg8Be6G<  
!取强度折减系数F=2.8时边坡围岩材料属性 &\^rQi/tf  
MP,EX,11,3E10                              ！定义弹性模量 {fU?idY)c  
MP,PRXY,11,0.25                             !定义泊松比 m2O&2[g  
MP,DENS,11,2500 W$'0Dc  
TB,DP,11                                   !定义数据表11 B?3juyB`--  
TBDATA,1,3.2E5,18.2                         !给数据表11输入数据 cIb4-TeV  
!取强度折减系数F=3.0时边坡围岩材料属性 Yb6q))Y  
MP,EX,12,3E10                              ！定义弹性模量 ZdH1nX(Yh3  
MP,PRXY,12,0.25     "r..  
MP,DENS,12,2500                            ！定义密度 CW-Ae  
TB,DP,12                                   !定义数据表12 L M /Ga  
TBDATA,1,3.0E5,17.1                        !给数据表12输入数据 \I:.<2i  
!建立几何模型 Za%LAyT_s  
!创建边坡线模型线 MMgx|"  
K,1,,,,                                      ！ 创建关键点1 2+:'0Krc  
K,2,-800,,, \[ M_\&GC  
K,3,-800,-800,,   Q'K$L9q  
K,4,-800,-1200,, _mk5^u/u  
K,5,1200,-1200,,   |dk[cX>  
K,6,1200,-800,,   Ri @`a  
K,7,1200,,, 3D_"y
Z  
K,8,1200,378,, q@ZlJ3%l,  
K,9,430,378,, smup,RNZRX  
l,1,2                                        !连接关键点1、2生成直线L1 O%g\B8;  
l,2,3  iSiDSeW8  
l,3,4 ;LQ9#M?  
l,4,5 /M@PO"  
l,5,6 X#ZgS!Mn  
l,6,7 Q["t eo]DQ  
l,7,8 s
i(cOCj/  
l,8,9 6P*O&1hv  
l,9,1 [/$N!2'5  
l,1,7 :1;Q(9:v  
l,3,6 d!V$Y}n  
!创建边坡面模型 QDE$E.a  
al,3,4,5,11                              !由线L3、L4、L5、L11生成一个面A1 y@;%Uv&  
al,1,2,11,6,10                            !由线L1、L2、L1、L6、L10生成一个面A2 Jhy(x1%  
al,7,8,9,10                              !由线L7、L8、L9、L10生成一个面A3 HnU Et/  
/PNUM, area, on                         !打开面号开关 d~Mg vh'  
!划分网格生成有限元模型 Mu2`ODe]  
! 划分边坡围岩2网格 1'O0`Me>#  
!设置网格份数 ;Z9IZ~  
lsel,s,,,3,5,2                               !选择线L3和L5 {Uq:Xw   
lesize,all,,,5                               !把所选择线分为5段 ;GOz>pg  
lsel,s,,,4,11,7                              !选择线L4和L11 Gx'TkU=  
lesize,all,,,26                              !把所选择线分为26段 nIBFk?)6  
mat,2                                    !给边坡围岩2赋予2号材料特性 N wtg%;  
type,1                                   !采用单元类型1 ;LMWNy4  
mshkey,0                                 !设定自由网格划分 :E'P7A  
mshape,0                                 !设定四边形网格划分 %Q~CB7ILK  
amesh,1                                  ！划分面积A1 rMp9jG@3   
! 划分边坡围岩1网格 u|fXP)>.  
!设置网格份数 _~M^ uW^l  
lsel,s,,,2,6,4                               !选择线L2和L6 wC` R>)  
lsel,a,,,1                                  !附加选择线L1 /!T> b:0  
lesize,all,,,10                              !把所选择线分为10段 o~1 Kp!U  
lsel,s,,,8,10,2                              !选择线L8和L10 > 0.W`j(s  
lesize,all,,,16                               !把所选择线分为16段 Cq\I''~8  
lsel,s,,,7,9,2 j*3}1L4P  
lesize,all,,,12 d8OL!Rk  
mat,1                                     !给边坡围岩1赋予1号材料特性 &leK}je [  
type,1 A-uB\ L  
mshkey,0                                  !设定自由网格划分
v\%G|8+]  
mshape,0                                  !设定四边形网格划分 RiR:69xwR*  
amesh,2                                   !划分面积A2 /3B6Mtb  
amesh,3                                   ！划分面积A3 T3pdx~66  
allsel XsL#;a C  
!保存网格模型 o#D.9K(  
save,Slope-grid.db s 3r=mp{  
(3)施加约束和荷载 o2uj =Gnx  
!给边坡模型施加约束 ,Z~;U  
!边坡两侧施加X方向约束 s9j7Psd  
nsel,s,loc,x,-800                             !选择X=-800线上所有节点 FxRXPt FK  
nsel,a,loc,x,1200                             !选择X=1200线上所有节点 ax|1b`XUr"  
d,all,ux                                     !对所选择节点约束X方向位移 *b)Q5dw@1  
allsel 0LL c 1t>}  
!边坡底部施加约束  }roG(  
nsel,s,loc,y,-1200                           !选择Y=--1200线上所有节点 {{!Y]\2S  
d,all,uy                                   !对所选择节点约束X、Y方向位移 NlF}{   
d,all,ux                                   2R.2D'4)`  
!施加重力加速度 Em^(  
acel,,9.8 \x+3f  
(4)求解 gEJi[E@  
/solu P Sx304  
!求解设置 76H!)={  
antype,static                                  !设定为静力求解 -iN.Iuc{b_  
nsubst,100                                   !设定最大子步数为100 .}`hCt08  
pred,on                                      !打开时间步长预测器 k6**u  
nropt,full                                     !设定牛顿-拉普森选项 :&w{\-0{  
nlgeom,on                                    !打开大位移效果 tIuoD+AW  
lnsrch,on                                     !打开线性搜索 %XTA;lrz  
outres,all,all                                  !输出所有项 lE5v-z? &|  
cnvtol,f,,0.005,2,0.5                            !力收敛准则设定 DQXS$uBT  
cnvtol,u,,0.05,2,1                              !位移收敛准则设定 :}q\tNY<  
!边坡在强度折减系数F=1时求解 ,
)fkr]`<  
allsel ]axh*J3`i  
solve                                        !进行求解 B'weok  
save,F1,db                                    !把F=1时求解结果保存 v:|(8Y  
!边坡在强度折减系数F=1.2时求解 q/gB<p9  
finish p{Sh F.  
/solu 5VE2@Fn}  
allsel w4:\N U  
asel,s,area,,2,3,1                               !选择面积A2、A3 hOk00az  
MPCHG,3,all                                 ! 把所选择单元材料号改为3 G|u3UhyB  
allsel 0l\y.   
solve                                         !进行求解 LE=k  
save,F1.2,db                                   !把F=1.2时求解结果保存                                                                                   q4k@l  
!边坡在强度折减系数F=1.4时求解 "Ty/k8?  
finish 39TT{>?`w  
/solu & >JDPB?5  
allsel b~0N^p[&%  
asel,s,area,,2,3,1                               !选择面积A2、A3 7Z,/g|s}z  
MPCHG,4,all                                 !把所选择单元材料号改为4 kQv*eZ~  
allsel \%A%s*1  
solve                                        !进行求解 A74920X`W  
save,F1.4,db                                  !把F=1.4时求解结果保存                             l!~
mxUb  
!边坡在强度折减系数F=1.6时求解 , c;eN  
finish Z)"61) )  
/solu bGXR7u&K  
allsel =VP=|g  
asel,s,area,,2,3,1                               !选择面积A2、A3 [_6&N.  
MPCHG,5,all                                 !把所选择单元材料号改为5 pW?&J>\6  
allsel #D%ygh=  
solve #-#NqX:  
save,F1.6,db                                   !把F=1.6时求解结果保存                                                                                   *nTU#U  
!边坡在强度折减系数F=1.8时求解 ,^'R_efY  
finish wts:65~  
/solu 2>PH8  
allsel 8cWZ"v  
asel,s,area,,2,3,1                                !选择面积A2、A3 ,?-\ x6  
MPCHG,6,all                                  !把所选择单元材料号改为6 v3RcwySk  
allsel  saZ>?Owz  
solve yy2I2Bv  
save,F1.8,db                                   !把F=1.8时求解结果保存   C!x/ ^gw  
!边坡在强度折减系数F=2.0时求解 } :?.>#  
finish 9w4sSj`  
/solu [`_io>*g  
allsel /Z%>ArAx  
asel,s,area,,2,3,1                                !选择面积A2、A3 !u;>Wyd W  
MPCHG,7,all                                   !把所选择单元材料号改为7 =(, ^du'  
allsel m <k!^jp  
solve `NNf&y)y  
save,F2.0,db                                   !把F=2.0时求解结果保存                             3Rv7Qx  
!边坡在强度折减系数F=2.2时求解 _YHu96H;  
finish ,^s  
/solu *~VxC{  
allsel S"2qJ!.u  
asel,s,area,,2,3,1                                !选择面积A2、A3 n \G Ry'  
MPCHG,8,all                                  !把所选择单元材料号改为8 8Jr?ZDf`  
allsel Ca>&  
solve (?&=T.*^  
save,F2.2,db                                   !把F=2.2时求解结果保存                                                                                   +Ov2`O8?  
!边坡在强度折减系数F=2.4时求解 Qvoqx>2p5  
finish #gbH^a'  
/solu SW?p?<  
allsel `$RA< 3  
asel,s,area,,2,3,1                                 !选择面积A2、A3 e$/B_o7(  
MPCHG,9,all                                   !把所选择单元材料号改为9 lPP,`  
allsel uOi&G:
=  
solve L2z2}U=<  
save,F2.4,db                                   !把F=2.4时求解结果保存   `)>7)={  
!边坡在强度折减系数F=2.6时求解 V64L,u#`
l  
finish NX6nQ  
/solu x5si70BKC/  
allsel  Qo
0H  
asel,s,area,,2,3,1                                 !选择面积A2、A3 "5|Lz)=  
MPCHG,10,all                                   !把所选择单元材料号改为10 LLAa1Wq  
allsel `iG,H[t+j  
solve |;vi*u  
save,F2.6,db                                   !把F=2.6时求解结果保存                             [J#1Ff;  
!边坡在强度折减系数F=2.8时求解 woOy*)@  
finish =_TaA(79  
/solu b\55,La  
allsel @9G- m(?*  
asel,s,area,,2,3,1                                !选择面积A2、A3 x8#ODuH  
MPCHG,11,all                                 !把所选择单元材料号改为11 xK%=  
allsel JiS5um=(.  
solve (jWss V1  
save,F2.8,db                                   !把F=2.8时求解结果保存                                                                                   :
os
z  
!边坡在强度折减系数F=3.0时求解 wzWbB2Mb5  
finish |$lwkC)O  
/solu ZU&"73   
allsel %:C ]7gQ  
asel,s,area,,2,3,1                             !选择面积A2、A3 P>QpvSd_#  
MPCHG,12,all                              !把所选择单元材料号改为12 .m4;^S2cO  
allsel l5t2\Fl  
solve 3tS~:6-/  
save,F3.0,db                                 !把F=3.0时求解结果保存 =bt]JRU  
/post1                                       ！进入后处理 IO$z%r7  
!边坡在强度折减系数F=1时结果分析 {{yt*7k{  
Resume,'F1','db'                              !读入边坡在强度折减系数F=1时 hZ\+FOx;  
set,1,last                                    !读入后一个子步 '_0]vupvY  
pldisp,1                                      !绘制边坡模型变形图 CC >=UF  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 (&\aA 0-}H  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 2ef;NC.&n  
!边坡在强度折减系数F=1.2时结果分析 _!V%fw  
Resume,'F1.2','db'                            !读入边坡在强度折减系数F=1.2时 $a;]_Y  
set,1,last                                    !读入后一个子步 dfj\RIV8  
pldisp,1                                      !绘制边坡模型变形图 f<jb=\}x  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 vP-M,4c  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 v Y|!  
!边坡在强度折减系数F=1.4时结果分析 JcO08
n  
Resume,'F1.4','db'                            !读入边坡在强度折减系数F=1.4时 k!XhFWb  
set,1,last                                    !读入后一个子步 T6#"8qz<  
pldisp,1                                      !绘制边坡模型变形图 &lc8G  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 6gO9 MQY  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图  q" @  
!边坡在强度折减系数F=1.6时结果分析 <I*x0BM=  
Resume,'F1.6','db'                            !读入边坡在强度折减系数F=1.6时 D$eB ,~  
set,1,last                                    !读入后一个子步 h rN%  
pldisp,1                                      !绘制边坡模型变形图 3ha|0[r9  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 XAOak$(j  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 E">T*ao  
!边坡在强度折减系数F=1.8时结果分析 rJfqA@  
Resume,'F1.8','db'                            !读入边坡在强度折减系数F=1.8时 pb;")Q'  
set,1,last                                    !读入后一个子步 #&!G"x
7  
pldisp,1                                      !绘制边坡模型变形图 %i{;r35M;9  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 Yn51U6_S  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 3NK ^AaTK  
!边坡在强度折减系数F=2.0时结果分析 "E!mva*NU  
Resume,'F2.0','db'                            !读入边坡在强度折减系数F=2.0时 P1zK2sL_  
set,1,last                                    !读入后一个子步 QGnxQ{ko  
pldisp,1                                      !绘制边坡模型变形图 eARk QV  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 j0-McLc  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 USML~]G z  
!边坡在强度折减系数F=2.2时结果分析 #TJk-1XM*q  
Resume,'F2.2','db'                            !读入边坡在强度折减系数F=2.2时 LX+5|u  
set,1,last                                    !读入后一个子步 e,1u  
pldisp,1                                      !绘制边坡模型变形图 ik1tidw  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 obClBO)@Y  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 3HV%4nZLf  
!边坡在强度折减系数F=2.4时结果分析 qB6dFl\ (  
Resume,'F2.4','db'                            !读入边坡在强度折减系数F=2.4时 /ONV5IkPy  
set,1,last                                    !读入后一个子步 /)|X.D  
pldisp,1                                      !绘制边坡模型变形图 9|r* pK[  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 MLVB^<qkeH  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 eH>#6R1-  
!边坡在强度折减系数F=2.6时结果分析 m7]hJ,0  
Resume,'F2.6','db'                            !读入边坡在强度折减系数F=2.6时 .q`{Dgc~  
set,1,last                                    !读入后一个子步 SqPtWEq@P  
pldisp,1                                      !绘制边坡模型变形图 / l>.mK()  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 i\}:hU-U  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 <_S@6?  
!边坡在强度折减系数F=2.8时结果分析 QOFvsJ<s  
Resume,'F2.8','db'                            !读入边坡在强度折减系数F=2.8时 'w/qcD-  
set,1,last                                    !读入后一个子步 UV7%4xM5v  
pldisp,1                                      !绘制边坡模型变形图 3SG?W_  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 ?^!,vh  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 nY-*
i!H  
!边坡在强度折减系数F=3.0时结果分析 m
S=r(3#  
Resume,'F3.0','db'                            !读入边坡在强度折减系数F=3.0时 TvhJVVQ+?  
set,1,last                                    !读入后一个子步 0OZMlt%z  
pldisp,1                                      !绘制边坡模型变形图 [2Mbk~  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 ReB7vpd  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图

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

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

好东西,顶一下!

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

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