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

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

同求 呵呵 谁知道啊

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

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

我已经自己解决了，哎

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

(1)创建物理环境 FEw51a+V  
/COM, Structural                             2I(0EBW  
/TITLE, Slope stability Analysis                 RD_&m?d  
/FILNAM,support,1                             Ku*@4#<L6h  
(2)建立模型 gfw,S;  
!进入前处理器 4\x'$G  
/PREP7 dY68wW>d|  
!定义单元类型 :Sk0?WU  
ET,1,PLANE82                               !定义边坡围岩单元 THWT\3~,  
KEYOPT,1,3,2                                !求解类型为平面应变 OqtQLqN  
!定义材料属性 =|bM|8,  
!边坡围岩1材料属性   t=NPo+fm  
MP,EX,1,3E10   9 bYoWw  
MP,PRXY,1,0.25     9kuL1tcY  
MP,DENS,1,2500                             !定义塑性模型 *TVr| to  
TB,DP,1                                     !定义数据表1 XL>Vwd  
TBDATA,1,9E5,42.8                           !给数据表输入数据 s\1h=V)!H  
!边坡围岩2材料属性                         !采用弹性模型 mm'n#%\G  
MP,EX,2,3.2E10                              ！定义弹性模量 <33[qt~  
MP,PRXY,2,0.24                              !定义泊松比 K_bF)6"  
MP,DENS,2,2700                             !定义密度 ^E8&!s  
!取强度折减系数F=1.2时边坡围岩材料属性 ~;QO`I=0P  
MP,EX,3,3E10   L
peQx\  
MP,PRXY,3,0.25     PQ<""_S||  
MP,DENS,3,2500 l|^p;z:d  
TB,DP,3                                     !定义数据表3 I1BVqIt1i  
TBDATA,1,7.5E5,37.7                         !给数据表输入数据 $)vljM<<  
!取强度折减系数F=1.4时边坡围岩材料属性 *L%HH@] %_  
MP,EX,4,3E10   FF6[qSV  
MP,PRXY,4,0.25                              !定义泊松比 YTpSR~!Rj  
MP,DENS,4,2500 =Bi>$Ly  
TB,DP,4                                    !定义数据表4 G$}\~dD  
TBDATA,1,6.4E5,33.5                        !给数据表输入数据 
]8*g%  
!取强度折减系数F=1.6时边坡围岩材料属性 )TFaG[tj  
MP,EX,5,3E10   +'2Mj|d@p  
MP,PRXY,5,0.25     VZ'[\3J  
MP,DENS,5,2500                             ！定义密度 DbP!wU lqR  
TB,DP,5                                    !定义数据表5 YwQxN"  
TBDATA,1,5.6E5,30                          !给数据表输入数据 hf^,  
!取强度折减系数F=1.8时边坡围岩材料属性 *4Y1((1k
 
MP,EX,6,3E10                               ！定义弹性模量 _-x|g~pV*  
MP,PRXY,6,0.25     vXLiYWo  
MP,DENS,6,2500 }RYr)  
TB,DP,6                                    !定义数据表6 63QMv[`,  
TBDATA,1,5.0E5,27.2 $GNN*WmHw  
!取强度折减系数F=2.0时边坡围岩材料属性 v#@"Evh7  
MP,EX,7,3E10   ~dC)EG  
MP,PRXY,7,0.25                             !定义泊松比 [!,&A{.!  
MP,DENS,7,2500 (Ybc~M)z  
TB,DP,7                                    !定义数据表7 z
*>"I  
TBDATA,1,4.5E5,24.8 X,:^})]  
!取强度折减系数F=2.2时边坡围岩材料属性 SN(:\|f 2  
MP,EX,8,3E10    @D^y<7(  
MP,PRXY,8,0.25     AcyiP   
MP,DENS,8,2500 @bOhnd#W  
TB,DP,8                                     !定义数据表8 6A;V[3  
TBDATA,1,4.09E5,22.8 EA|*|o4)  
!取强度折减系数F=2.4时边坡围岩材料属性 HsGXb\  
MP,EX,9,3E10   }L>}_NV\  
MP,PRXY,9,0.25     #Z)e]4{!l  
MP,DENS,9,2500 @X?DHLM  
TB,DP,9                                      !定义数据表9 m{x[q  
TBDATA,1,3.6E5,21.1                          !给数据表输入数据 OGh9^,v  
!取强度折减系数F=2.6时边坡围岩材料属性 S=~[6;G  
MP,EX,10,3E10                                ！定义弹性模量   v*<hE>J0  
MP,PRXY,10,0.25                              !定义泊松比 h^D?G2O  
MP,DENS,10,2500                           ！定义密度 jxL}tS{j  
TB,DP,10                                   !定义数据表10 RmY5/IYR|:  
TBDATA,1,3.46E5,19.6                        !给数据表10输入数据 a&6 3[p.<}  
!取强度折减系数F=2.8时边坡围岩材料属性 g2JNa?z  
MP,EX,11,3E10                              ！定义弹性模量 w<<>XIL  
MP,PRXY,11,0.25                             !定义泊松比 G`
9Ud  
MP,DENS,11,2500 /T]2ZX>  
TB,DP,11                                   !定义数据表11 Dz
:A.x@$*  
TBDATA,1,3.2E5,18.2                         !给数据表11输入数据 H ifKa/}P8  
!取强度折减系数F=3.0时边坡围岩材料属性 21bvSK  
MP,EX,12,3E10                              ！定义弹性模量 G)]'>m<y  
MP,PRXY,12,0.25     Ii%^z?'  
MP,DENS,12,2500                            ！定义密度 K>l$Y#x}k  
TB,DP,12                                   !定义数据表12 B BbGq8p  
TBDATA,1,3.0E5,17.1                        !给数据表12输入数据 wUW^ O  
!建立几何模型 A&jkc'  
!创建边坡线模型线 rS\j9@=Y4  
K,1,,,,                                      ！ 创建关键点1 E'j>[C:U  
K,2,-800,,, fPZt*A__  
K,3,-800,-800,,   ZZ?
0
%9  
K,4,-800,-1200,, 0z #'=XWk  
K,5,1200,-1200,,   E?z3 D*U  
K,6,1200,-800,,   _?M34&.X  
K,7,1200,,, F$ShhZgi  
K,8,1200,378,, tisSj?+  
K,9,430,378,, V$VqYy9 *  
l,1,2                                        !连接关键点1、2生成直线L1 "}!|V)K
 
l,2,3 9cp-Rw<tI  
l,3,4 ci0)kxUBF  
l,4,5 Urj8v2k  
l,5,6 >N62t9Ll[  
l,6,7 Xt^ldW  
l,7,8 a$yAF4HR
<  
l,8,9 c [sydl  
l,9,1 aTuD|s  
l,1,7 UBzX%:A  
l,3,6 9u^PM  
!创建边坡面模型 Z,)4(#b =  
al,3,4,5,11                              !由线L3、L4、L5、L11生成一个面A1 ~m8".Z"  
al,1,2,11,6,10                            !由线L1、L2、L1、L6、L10生成一个面A2 &YGd!Q  
al,7,8,9,10                              !由线L7、L8、L9、L10生成一个面A3 0f&B;?)!  
/PNUM, area, on                         !打开面号开关 ;e415T  
!划分网格生成有限元模型 .LhIB?  
! 划分边坡围岩2网格 9+nB;vA  
!设置网格份数 z85%2Apd  
lsel,s,,,3,5,2                               !选择线L3和L5 Ci4`,  
lesize,all,,,5                               !把所选择线分为5段 juG?kL.  
lsel,s,,,4,11,7                              !选择线L4和L11 d&4ve Lu  
lesize,all,,,26                              !把所选择线分为26段 }pdn-#  
mat,2                                    !给边坡围岩2赋予2号材料特性 M(KsLu1   
type,1                                   !采用单元类型1 P}29wrIZ  
mshkey,0                                 !设定自由网格划分 ,
7KP  
mshape,0                                 !设定四边形网格划分 8om6wALXB  
amesh,1                                  ！划分面积A1 F&%@p&  
! 划分边坡围岩1网格 ,Md8A`7x~  
!设置网格份数 ztTj2M"  
lsel,s,,,2,6,4                               !选择线L2和L6 $wg5q\Rv  
lsel,a,,,1                                  !附加选择线L1 &l"/G%W  
lesize,all,,,10                              !把所选择线分为10段 N4I`6uDgD  
lsel,s,,,8,10,2                              !选择线L8和L10 jzI70+E  
lesize,all,,,16                               !把所选择线分为16段 V zuW]"  
lsel,s,,,7,9,2 >!848J  
lesize,all,,,12 :m]~o3KRy  
mat,1                                     !给边坡围岩1赋予1号材料特性 rn $a)^!  
type,1 f6vhW66:?x  
mshkey,0                                  !设定自由网格划分 h:-ZXIv?  
mshape,0                                  !设定四边形网格划分 njtz,qt_;G  
amesh,2                                   !划分面积A2 &a5UQ>  
amesh,3                                   ！划分面积A3 B(7
oHj.i2  
allsel O;z:?  
!保存网格模型 "XfCLc1 T  
save,Slope-grid.db Xyy;BO:  
(3)施加约束和荷载 y$|%K3  
!给边坡模型施加约束 i'OFun+-,  
!边坡两侧施加X方向约束 kd)Q$RA(  
nsel,s,loc,x,-800                             !选择X=-800线上所有节点 px8988X  
nsel,a,loc,x,1200                             !选择X=1200线上所有节点 >lQ@" U  
d,all,ux                                     !对所选择节点约束X方向位移 XLb lVi@
 
allsel c[J?`8  
!边坡底部施加约束 g>-pC a  
nsel,s,loc,y,-1200                           !选择Y=--1200线上所有节点 *`Swv`  
d,all,uy                                   !对所选择节点约束X、Y方向位移 3O7]~5 j1  
d,all,ux                                   `ltc)$  
!施加重力加速度 b3F)$UQ  
acel,,9.8 FM;NA{  
(4)求解 -0r0M)  
/solu _8A  
!求解设置 v/*}M&vo  
antype,static                                  !设定为静力求解 g"ha1<y<  
nsubst,100                                   !设定最大子步数为100 h/5|3  
pred,on                                      !打开时间步长预测器 r*HbglB  
nropt,full                                     !设定牛顿-拉普森选项 I{tY;b'w  
nlgeom,on                                    !打开大位移效果 #%N v\g;  
lnsrch,on                                     !打开线性搜索 `-fWN
Hs  
outres,all,all                                  !输出所有项 ^MIF+/bQ  
cnvtol,f,,0.005,2,0.5                            !力收敛准则设定 G d~ v _  
cnvtol,u,,0.05,2,1                              !位移收敛准则设定 N;4bEcWjp  
!边坡在强度折减系数F=1时求解 %c"PMTq(  
allsel S"'0lS   
solve                                        !进行求解 7rQwn2XD{  
save,F1,db                                    !把F=1时求解结果保存 @&?E3?5ll  
!边坡在强度折减系数F=1.2时求解 mivb}cKM  
finish `|coA2$rw  
/solu rV84?75(Y  
allsel u^|c_5J(  
asel,s,area,,2,3,1                               !选择面积A2、A3 <}
t~^E,  
MPCHG,3,all                                 ! 把所选择单元材料号改为3 fb-Lp#!T39  
allsel J9eOBom8e<  
solve                                         !进行求解 q
;Tdqv!Ju  
save,F1.2,db                                   !把F=1.2时求解结果保存                                                                                   YbtsJ <w  
!边坡在强度折减系数F=1.4时求解 .DsdQ4Y  
finish g xY6M4  
/solu 1/+d@s#t  
allsel ]xC#XYE:dy  
asel,s,area,,2,3,1                               !选择面积A2、A3  9uR+  
MPCHG,4,all                                 !把所选择单元材料号改为4 w\,N}'
G  
allsel 3.+TM]RYN  
solve                                        !进行求解 \ne1Xu:hM  
save,F1.4,db                                  !把F=1.4时求解结果保存                             .7&V@A7  
!边坡在强度折减系数F=1.6时求解 g%Bh-O9\  
finish dp
#JvZb  
/solu ve($l"T  
allsel 7f|8SB  
asel,s,area,,2,3,1                               !选择面积A2、A3 ${m;x:'  
MPCHG,5,all                                 !把所选择单元材料号改为5 
?lq  
allsel M2H +1ic  
solve lC/1,Z/M  
save,F1.6,db                                   !把F=1.6时求解结果保存                                                                                   uonCD
8  
!边坡在强度折减系数F=1.8时求解 "@^Pb$BLY  
finish 2t 6m#  
/solu %]
7'2  
allsel DmU,}]#:  
asel,s,area,,2,3,1                                !选择面积A2、A3 'L?e)u.  
MPCHG,6,all                                  !把所选择单元材料号改为6 K+MSjQS"
 
allsel fh_+M"Y0`  
solve r5 tn'  
save,F1.8,db                                   !把F=1.8时求解结果保存   -!;2?6R9{  
!边坡在强度折减系数F=2.0时求解 >*s_)IH2  
finish O+x"c3@Z)D  
/solu EP,j+^RVf  
allsel $`j%z@[g  
asel,s,area,,2,3,1                                !选择面积A2、A3 X3e&c  
MPCHG,7,all                                   !把所选择单元材料号改为7 ,1/O2aQ%\0  
allsel kvcDa+#  
solve ]H|O  
save,F2.0,db                                   !把F=2.0时求解结果保存                             scX'>\w&c  
!边坡在强度折减系数F=2.2时求解 <G|i!Pm  
finish g>G+?PY  
/solu j5m KJC  
allsel m}A|W[p<  
asel,s,area,,2,3,1                                !选择面积A2、A3 !yqez  
MPCHG,8,all                                  !把所选择单元材料号改为8 TOapq9B]  
allsel "Vh3hnS~  
solve ^POHQQ  
save,F2.2,db                                   !把F=2.2时求解结果保存                                                                                   \]C_ul'  
!边坡在强度折减系数F=2.4时求解 V%h,JA  
finish "uCO?hv0  
/solu XSk*w'xO  
allsel #1*
#3p9UL  
asel,s,area,,2,3,1                                 !选择面积A2、A3 $mxm?7ZVR  
MPCHG,9,all                                   !把所选择单元材料号改为9 QO2@K1Y  
allsel `#HtVI  
solve q9Zp8&<EqH  
save,F2.4,db                                   !把F=2.4时求解结果保存   3LGX ^J<f  
!边坡在强度折减系数F=2.6时求解 L$L/5/  
finish  _
U.|$pU  
/solu yPY}b_W  
allsel A0hKzj  
asel,s,area,,2,3,1                                 !选择面积A2、A3 /2^"c+/'p  
MPCHG,10,all                                   !把所选择单元材料号改为10 C7*n<+e  
allsel ]%M&pc3U  
solve :I_p4S.)  
save,F2.6,db                                   !把F=2.6时求解结果保存                             ]1hyvm3  
!边坡在强度折减系数F=2.8时求解 r$[`A_  
finish /pY-how%!  
/solu Gyx4}pV  
allsel GDF/0-/Z  
asel,s,area,,2,3,1                                !选择面积A2、A3 /tm2b<G  
MPCHG,11,all                                 !把所选择单元材料号改为11 Kzwbr?&z  
allsel m)AF9#aT2  
solve 1WjNFi  
save,F2.8,db                                   !把F=2.8时求解结果保存                                                                                   4Vtug>  
!边坡在强度折减系数F=3.0时求解 (#kKL??W  
finish 1lo.X_  
/solu Hjhgu=  
allsel Q$+6f,m#W  
asel,s,area,,2,3,1                             !选择面积A2、A3 '8J!(+  
MPCHG,12,all                              !把所选择单元材料号改为12 +YZ*>ki  
allsel YRg"{[+#]k  
solve #[f]-c(!  
save,F3.0,db                                 !把F=3.0时求解结果保存 I3gl+)Q  
/post1                                       ！进入后处理 .~$!BWP  
!边坡在强度折减系数F=1时结果分析 Q g~cYwX  
Resume,'F1','db'                              !读入边坡在强度折减系数F=1时 {p\ll  
set,1,last                                    !读入后一个子步 |RjAp.pm  
pldisp,1                                      !绘制边坡模型变形图  I~T   
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 )<Fq}Q86  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 Ri>4:V3K  
!边坡在强度折减系数F=1.2时结果分析 4)"S/u  
Resume,'F1.2','db'                            !读入边坡在强度折减系数F=1.2时 nTsKJX%\  
set,1,last                                    !读入后一个子步 %xt;&HE  
pldisp,1                                      !绘制边坡模型变形图 PRJ  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 >?<S(  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 cWtuI(.  
!边坡在强度折减系数F=1.4时结果分析 Tp46K\}Uf  
Resume,'F1.4','db'                            !读入边坡在强度折减系数F=1.4时 /!Ay12lKE}  
set,1,last                                    !读入后一个子步 +*RaX (&  
pldisp,1                                      !绘制边坡模型变形图 i<0_sxfUD  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 mR|L'[l  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 ^H"o=K8=  
!边坡在强度折减系数F=1.6时结果分析 Ml_Hq>
\U  
Resume,'F1.6','db'                            !读入边坡在强度折减系数F=1.6时 &F- \t5X=i  
set,1,last                                    !读入后一个子步 D:] QBA)C  
pldisp,1                                      !绘制边坡模型变形图 QPX&P{!g  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 wE[gp+X~  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 a\m_Q{:  
!边坡在强度折减系数F=1.8时结果分析 O22Q g  
Resume,'F1.8','db'                            !读入边坡在强度折减系数F=1.8时 P~ : N  
set,1,last                                    !读入后一个子步 6VUs:iO1j5  
pldisp,1                                      !绘制边坡模型变形图 9xi nX-x;n  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 KH$|
wv  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 5P Zzaz<  
!边坡在强度折减系数F=2.0时结果分析 s&hJ[$i  
Resume,'F2.0','db'                            !读入边坡在强度折减系数F=2.0时 E5aRTDLq  
set,1,last                                    !读入后一个子步 QyghNImp  
pldisp,1                                      !绘制边坡模型变形图 K;z$~;F  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图
(}g4}A@x  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 R + ~b@  
!边坡在强度折减系数F=2.2时结果分析 H;nEU@>"Z  
Resume,'F2.2','db'                            !读入边坡在强度折减系数F=2.2时 = N&5]Z  
set,1,last                                    !读入后一个子步 'C4cS[1  
pldisp,1                                      !绘制边坡模型变形图 kEgpF{"%n  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 s$%t2UaV  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 M*!WXQlud  
!边坡在强度折减系数F=2.4时结果分析 =yF]#>Ah  
Resume,'F2.4','db'                            !读入边坡在强度折减系数F=2.4时 0$A^ .M;  
set,1,last                                    !读入后一个子步 :V3z`}Rl  
pldisp,1                                      !绘制边坡模型变形图 Hf/ZaBn  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 }y%c.  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 JDJ"D\85  
!边坡在强度折减系数F=2.6时结果分析 J>l?HK  
Resume,'F2.6','db'                            !读入边坡在强度折减系数F=2.6时 l1bkhA b  
set,1,last                                    !读入后一个子步 fmq''1u  
pldisp,1                                      !绘制边坡模型变形图 Y~xo=v(  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 K| dI'TnW  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 (> W\Nf  
!边坡在强度折减系数F=2.8时结果分析 44NMof8N  
Resume,'F2.8','db'                            !读入边坡在强度折减系数F=2.8时 l~]D|92  
set,1,last                                    !读入后一个子步 Um}f7^fp^l  
pldisp,1                                      !绘制边坡模型变形图 +0}z3T1L  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 LZ34x: ,C  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 SR$ 'JGfp  
!边坡在强度折减系数F=3.0时结果分析 ;NOmI+t0w&  
Resume,'F3.0','db'                            !读入边坡在强度折减系数F=3.0时 7s0\`eXo/  
set,1,last                                    !读入后一个子步 ;,8 )%
[  
pldisp,1                                      !绘制边坡模型变形图 =cpUc]~  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 r<$o [,W  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图

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

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

好东西,顶一下!

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

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