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

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

同求 呵呵 谁知道啊

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

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

我已经自己解决了，哎

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

(1)创建物理环境 _G|6xlO  
/COM, Structural                             9^tyjX2  
/TITLE, Slope stability Analysis                 p pq#5t^[)  
/FILNAM,support,1                             C#R9Hlb  
(2)建立模型 14R))Dz"  
!进入前处理器 n9^zAcUbAW  
/PREP7 r[~$  
!定义单元类型 o%a$m9I  
ET,1,PLANE82                               !定义边坡围岩单元 &AVi4zV  
KEYOPT,1,3,2                                !求解类型为平面应变 A5b}G
 
!定义材料属性 qz&)|~,\C  
!边坡围岩1材料属性   8TZe=sD~cr  
MP,EX,1,3E10   nj  
MP,PRXY,1,0.25     g d-fJ._1  
MP,DENS,1,2500                             !定义塑性模型 4]GyuY  
TB,DP,1                                     !定义数据表1 5{> cfN\q  
TBDATA,1,9E5,42.8                           !给数据表输入数据 h! <8=V(  
!边坡围岩2材料属性                         !采用弹性模型 m[f\I^\%8  
MP,EX,2,3.2E10                              ！定义弹性模量 q'q{M-U<  
MP,PRXY,2,0.24                              !定义泊松比 vY6
|V$  
MP,DENS,2,2700                             !定义密度 )j]RFt  
!取强度折减系数F=1.2时边坡围岩材料属性 xjpW<-)MLf  
MP,EX,3,3E10   Lnzhs;7L  
MP,PRXY,3,0.25     By9/tB  
MP,DENS,3,2500 r[x7?cXsW  
TB,DP,3                                     !定义数据表3 `*a,8M%  
TBDATA,1,7.5E5,37.7                         !给数据表输入数据 5tL6R3  
!取强度折减系数F=1.4时边坡围岩材料属性 m78PQx H  
MP,EX,4,3E10   .z"[z^/uF  
MP,PRXY,4,0.25                              !定义泊松比 n|.;g!QDA  
MP,DENS,4,2500 T"jl;,gr]J  
TB,DP,4                                    !定义数据表4 C0M{zGT>}  
TBDATA,1,6.4E5,33.5                        !给数据表输入数据 LFC k
6 R  
!取强度折减系数F=1.6时边坡围岩材料属性 YS*t
7  
MP,EX,5,3E10   !IS,[  
MP,PRXY,5,0.25     oS4ag  
MP,DENS,5,2500                             ！定义密度 c LJCLKJ  
TB,DP,5                                    !定义数据表5 >/*\xg&J  
TBDATA,1,5.6E5,30                          !给数据表输入数据 'zaB5d~l  
!取强度折减系数F=1.8时边坡围岩材料属性 <#UvLll  
MP,EX,6,3E10                               ！定义弹性模量 R)=<q]Ms  
MP,PRXY,6,0.25     `t -3(>P  
MP,DENS,6,2500 ?:E;C<Ar  
TB,DP,6                                    !定义数据表6 h| Ih4  
TBDATA,1,5.0E5,27.2 vuf|2!kh/  
!取强度折减系数F=2.0时边坡围岩材料属性 ~U|te_l  
MP,EX,7,3E10   P_gQ-pF.  
MP,PRXY,7,0.25                             !定义泊松比 >Vvc55z  
MP,DENS,7,2500 !ktr|9Bl  
TB,DP,7                                    !定义数据表7 Evc 9k  
TBDATA,1,4.5E5,24.8 ;g9+*$Gw  
!取强度折减系数F=2.2时边坡围岩材料属性 &}r932  
MP,EX,8,3E10   ;#due  
MP,PRXY,8,0.25     qA30G~S  
MP,DENS,8,2500 ,x/j&S9!  
TB,DP,8                                     !定义数据表8 O_ cK4  
TBDATA,1,4.09E5,22.8 "'Q:%_;  
!取强度折减系数F=2.4时边坡围岩材料属性 0U<9=[~q7@  
MP,EX,9,3E10   ]x|sTKv2  
MP,PRXY,9,0.25     2+.m44>Ti  
MP,DENS,9,2500 jcj)9;n=!  
TB,DP,9                                      !定义数据表9 z!%}0  
TBDATA,1,3.6E5,21.1                          !给数据表输入数据 >&R|t_ypw  
!取强度折减系数F=2.6时边坡围岩材料属性 e#wn;w
o?  
MP,EX,10,3E10                                ！定义弹性模量   .JqIAC~  
MP,PRXY,10,0.25                              !定义泊松比 $f+
9svq  
MP,DENS,10,2500                           ！定义密度 .o>QBYpTw/  
TB,DP,10                                   !定义数据表10 ":Q^/;D}U  
TBDATA,1,3.46E5,19.6                        !给数据表10输入数据 RwE]t$T/  
!取强度折减系数F=2.8时边坡围岩材料属性 <bH>\@p7}  
MP,EX,11,3E10                              ！定义弹性模量 [o~w>,a  
MP,PRXY,11,0.25                             !定义泊松比 Z&%61jGK  
MP,DENS,11,2500 ,<BTv;4p  
TB,DP,11                                   !定义数据表11 3-05y!vbcE  
TBDATA,1,3.2E5,18.2                         !给数据表11输入数据 |Bid(`t.  
!取强度折减系数F=3.0时边坡围岩材料属性 +vP1DXtj(  
MP,EX,12,3E10                              ！定义弹性模量 0czy:d,M%  
MP,PRXY,12,0.25     X:A\{^~  
MP,DENS,12,2500                            ！定义密度 LYX+/@OU2  
TB,DP,12                                   !定义数据表12 >nxtQ  
TBDATA,1,3.0E5,17.1                        !给数据表12输入数据 a'Z"Yz^Eo  
!建立几何模型 d={}a,3?  
!创建边坡线模型线 ktCh*R[`  
K,1,,,,                                      ！ 创建关键点1 V;!D:N8<  
K,2,-800,,, ~VOmMw4HV
 
K,3,-800,-800,,   aF:I]]TfK~  
K,4,-800,-1200,, G4i&:0  
K,5,1200,-1200,,   1\McsX4  
K,6,1200,-800,,   4{Iz\:G:{/  
K,7,1200,,, G9 !1Wz
s  
K,8,1200,378,, )JPcSy*  
K,9,430,378,, }7V/(K  
l,1,2                                        !连接关键点1、2生成直线L1 Wg[`H=)Q  
l,2,3 ;8@A7`^  
l,3,4 t`?FSV  
l,4,5 ,oCr6 ]  
l,5,6 Q7C'O @  
l,6,7 i< ih :  
l,7,8 &Wba2fD  
l,8,9 _ |; bh  
l,9,1 mXr)lA  
l,1,7 nT>?}/S  
l,3,6 &zZSWNW  
!创建边坡面模型 Oj:`r*z43  
al,3,4,5,11                              !由线L3、L4、L5、L11生成一个面A1 ^%L$$V nG  
al,1,2,11,6,10                            !由线L1、L2、L1、L6、L10生成一个面A2 AXyuXB  
al,7,8,9,10                              !由线L7、L8、L9、L10生成一个面A3 (n<xoV[e  
/PNUM, area, on                         !打开面号开关 &^EkM  
!划分网格生成有限元模型 46vz=# ,6L  
! 划分边坡围岩2网格 X7G6y|4;w  
!设置网格份数 0ode&dB  
lsel,s,,,3,5,2                               !选择线L3和L5 {XVSHUtw  
lesize,all,,,5                               !把所选择线分为5段 Q_dFZ  
lsel,s,,,4,11,7                              !选择线L4和L11 eg3{sDv,  
lesize,all,,,26                              !把所选择线分为26段 P|\,kw>l  
mat,2                                    !给边坡围岩2赋予2号材料特性 (w.B_9#  
type,1                                   !采用单元类型1 Y4_i=}\*vf  
mshkey,0                                 !设定自由网格划分 pO^ 6p%  
mshape,0                                 !设定四边形网格划分 5XhV+t g.  
amesh,1                                  ！划分面积A1 (<ejJPWT  
! 划分边坡围岩1网格 4&\m!s  
!设置网格份数 vq{:=:5'P  
lsel,s,,,2,6,4                               !选择线L2和L6 @*oi1
_q  
lsel,a,,,1                                  !附加选择线L1 ,FTF@h-Cs  
lesize,all,,,10                              !把所选择线分为10段 loZfzN&6A  
lsel,s,,,8,10,2                              !选择线L8和L10 */1z=  
lesize,all,,,16                               !把所选择线分为16段 Na=q(OKN  
lsel,s,,,7,9,2 JwO+Dd  
lesize,all,,,12 ukw'$Yt2  
mat,1                                     !给边坡围岩1赋予1号材料特性 m*'#`vIbb  
type,1 V{n
pK(  
mshkey,0                                  !设定自由网格划分 %63<Iz"  
mshape,0                                  !设定四边形网格划分 ?$ 3=m)s  
amesh,2                                   !划分面积A2 43eGfp'  
amesh,3                                   ！划分面积A3 b7$?'neH/.  
allsel gnv4.f:  
!保存网格模型 CB~&!MdMr  
save,Slope-grid.db [L8gG.wy  
(3)施加约束和荷载 cX*^PSM  
!给边坡模型施加约束 3laSPih[.  
!边坡两侧施加X方向约束 u^ T2  
nsel,s,loc,x,-800                             !选择X=-800线上所有节点 PtHT>  
nsel,a,loc,x,1200                             !选择X=1200线上所有节点 T:si?7CR  
d,all,ux                                     !对所选择节点约束X方向位移 ]?kf;A@  
allsel 0<Y)yNsV  
!边坡底部施加约束 ':Te#S  
nsel,s,loc,y,-1200                           !选择Y=--1200线上所有节点 0iYe>u  
d,all,uy                                   !对所选择节点约束X、Y方向位移 K= 69z  
d,all,ux                                   xZkLN5I{  
!施加重力加速度 ~"-wS
Am  
acel,,9.8 csC3Wm{v  
(4)求解 sB6UlX;b:  
/solu Z5+0?X0i  
!求解设置 ''Hq-Ng  
antype,static                                  !设定为静力求解 ISl'g'o  
nsubst,100                                   !设定最大子步数为100 6ul34\;  
pred,on                                      !打开时间步长预测器 /XXW4_>  
nropt,full                                     !设定牛顿-拉普森选项 I=1tf;Bsi  
nlgeom,on                                    !打开大位移效果 `);`E_'U k  
lnsrch,on                                     !打开线性搜索 AOTI&v  
outres,all,all                                  !输出所有项 D@2Tx  
cnvtol,f,,0.005,2,0.5                            !力收敛准则设定 +h =lAHn&  
cnvtol,u,,0.05,2,1                              !位移收敛准则设定 o5)U3U1|  
!边坡在强度折减系数F=1时求解 {DpZg",H-  
allsel A`@we  
solve                                        !进行求解 FOZqN K  
save,F1,db                                    !把F=1时求解结果保存 y]/{W}D  
!边坡在强度折减系数F=1.2时求解 ^}WeBU  
finish os>|LPv4  
/solu 6)TFb,  
allsel 9TF[uC)-2  
asel,s,area,,2,3,1                               !选择面积A2、A3 n_aKciF  
MPCHG,3,all                                 ! 把所选择单元材料号改为3 QoqdPk#1  
allsel !u^(<.xJ   
solve                                         !进行求解 xTe?*  
save,F1.2,db                                   !把F=1.2时求解结果保存                                                                                   0q/g:"|j  
!边坡在强度折减系数F=1.4时求解 OvFZ&S[  
finish ,xGlWH wrY  
/solu O6`@'N>6P  
allsel M?
_VYK  
asel,s,area,,2,3,1                               !选择面积A2、A3 *P_TG"^{W  
MPCHG,4,all                                 !把所选择单元材料号改为4 03MB,  
allsel
"8-]6p3u  
solve                                        !进行求解 s~GO-v7  
save,F1.4,db                                  !把F=1.4时求解结果保存                             r6^DD$X  
!边坡在强度折减系数F=1.6时求解 ON=xn|b4  
finish 0c]Lm?&  
/solu TPkm~>zD.  
allsel 6gp3n;D  
asel,s,area,,2,3,1                               !选择面积A2、A3
xT@\FwPr  
MPCHG,5,all                                 !把所选择单元材料号改为5 l_8t[  
allsel c<e$6:|xM  
solve s?=J#WV1y  
save,F1.6,db                                   !把F=1.6时求解结果保存                                                                                   y"7?]#$9/  
!边坡在强度折减系数F=1.8时求解 -i58FJ`B  
finish 6rRPqO j  
/solu _-EHG  
allsel jtZ@`io  
asel,s,area,,2,3,1                                !选择面积A2、A3 t+vn.X+&  
MPCHG,6,all                                  !把所选择单元材料号改为6 5!?><{k=%  
allsel q* m%Fv  
solve 6Up,B=sX0  
save,F1.8,db                                   !把F=1.8时求解结果保存   )q#b^( v  
!边坡在强度折减系数F=2.0时求解 w_9:gprf  
finish %1#5 7-  
/solu 0s4%22  
allsel hX;xbl  
asel,s,area,,2,3,1                                !选择面积A2、A3 tUtl>>6Iu  
MPCHG,7,all                                   !把所选择单元材料号改为7 HMBxj($eR  
allsel g
SP|;Gy  
solve r+) A)a,  
save,F2.0,db                                   !把F=2.0时求解结果保存                             xbIxtZm  
!边坡在强度折减系数F=2.2时求解 <uA|nYpp  
finish c=ZX7U  
/solu 7OT}V}iP  
allsel E;h#3 B9  
asel,s,area,,2,3,1                                !选择面积A2、A3 3i7n"8\$  
MPCHG,8,all                                  !把所选择单元材料号改为8 rtY0?  
allsel Jx'p\*  
solve n&@\[,B  
save,F2.2,db                                   !把F=2.2时求解结果保存                                                                                   Q<"zpwHR  
!边坡在强度折减系数F=2.4时求解 Qd@`jwjS  
finish f$P pFSY4  
/solu vHao y  
allsel g6N{Z e Wg  
asel,s,area,,2,3,1                                 !选择面积A2、A3 50CU|  
MPCHG,9,all                                   !把所选择单元材料号改为9 4Jw0m#UN1  
allsel TIS}'c'C  
solve t.]oLG22r  
save,F2.4,db                                   !把F=2.4时求解结果保存   ><$hFrR!  
!边坡在强度折减系数F=2.6时求解 qD%Jf4.0j  
finish f~E'0f_  
/solu W1Ht8uYG3  
allsel M'*
 Y  
asel,s,area,,2,3,1                                 !选择面积A2、A3 hG3b7!^#g  
MPCHG,10,all                                   !把所选择单元材料号改为10 gllXJM^ -  
allsel u%&zY97/  
solve = uOFaZ4  
save,F2.6,db                                   !把F=2.6时求解结果保存                             w;X-i.%`  
!边坡在强度折减系数F=2.8时求解 qgu.c`GmW  
finish T9u/|OP  
/solu .>&kAf.  
allsel B=9|g1e  
asel,s,area,,2,3,1                                !选择面积A2、A3 "5vFa7y  
MPCHG,11,all                                 !把所选择单元材料号改为11 |vzGFfRI  
allsel #w#B'  
solve Fm*O&6W\@A  
save,F2.8,db                                   !把F=2.8时求解结果保存                                                                                   w
KwireOs  
!边坡在强度折减系数F=3.0时求解 s7=]!7QGS!  
finish '*22j ]  
/solu p#3P`I>ZrT  
allsel rQ/S|gG  
asel,s,area,,2,3,1                             !选择面积A2、A3 lGs fs(  
MPCHG,12,all                              !把所选择单元材料号改为12 *
F&C`]  
allsel %[RLc[pB  
solve O10h(Wg  
save,F3.0,db                                 !把F=3.0时求解结果保存 pTcm2-J  
/post1                                       ！进入后处理 #.) qQ8*(  
!边坡在强度折减系数F=1时结果分析 aG,N>0k8  
Resume,'F1','db'                              !读入边坡在强度折减系数F=1时 wWwY.}j  
set,1,last                                    !读入后一个子步 Y(<>[8S m  
pldisp,1                                      !绘制边坡模型变形图 lqu1H&  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 N2C^'dFj  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 nE,"3X"  
!边坡在强度折减系数F=1.2时结果分析 XO\P4x:c  
Resume,'F1.2','db'                            !读入边坡在强度折减系数F=1.2时 _w(SHWh2  
set,1,last                                    !读入后一个子步 8on2BC2  
pldisp,1                                      !绘制边坡模型变形图 (zUERw\aX  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 p7|~x@q+  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 YT#
"HYO  
!边坡在强度折减系数F=1.4时结果分析 :U?Kwv8s  
Resume,'F1.4','db'                            !读入边坡在强度折减系数F=1.4时 [_${N,1  
set,1,last                                    !读入后一个子步 `57ffQR9  
pldisp,1                                      !绘制边坡模型变形图 r]2}S=[  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 Dtelr=/s  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 V" I+E  
!边坡在强度折减系数F=1.6时结果分析 p=eSHs{>A  
Resume,'F1.6','db'                            !读入边坡在强度折减系数F=1.6时 xAsbP$J:  
set,1,last                                    !读入后一个子步 c>I^SY(r%  
pldisp,1                                      !绘制边坡模型变形图 KO:o GUR  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 mw.9cDf  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 h4ZrD:D0\  
!边坡在强度折减系数F=1.8时结果分析 .4U::j}  
Resume,'F1.8','db'                            !读入边坡在强度折减系数F=1.8时 2c}>}A4  
set,1,last                                    !读入后一个子步 #VD[\#  
pldisp,1                                      !绘制边坡模型变形图 MA"DP7e?v  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 A`IE8@&Z'  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 M7En%sBp  
!边坡在强度折减系数F=2.0时结果分析 !30BZM^  
Resume,'F2.0','db'                            !读入边坡在强度折减系数F=2.0时 yay{lP}b"  
set,1,last                                    !读入后一个子步 BLfoU_Z  
pldisp,1                                      !绘制边坡模型变形图 g,9o'
fs`x  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 J5I
Q  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 J8(v65  
!边坡在强度折减系数F=2.2时结果分析 is`le}$^y  
Resume,'F2.2','db'                            !读入边坡在强度折减系数F=2.2时 U2!9Tl9".  
set,1,last                                    !读入后一个子步 5y@JMQSO  
pldisp,1                                      !绘制边坡模型变形图 #Ey!
?Z  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 3TO$J  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 EpS"NQEe  
!边坡在强度折减系数F=2.4时结果分析 !x|Ok'izDL  
Resume,'F2.4','db'                            !读入边坡在强度折减系数F=2.4时 YwEXTy>0  
set,1,last                                    !读入后一个子步 Ao 1*a%-.  
pldisp,1                                      !绘制边坡模型变形图 )x#^fN~ 7`  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 DaaLRMQ=  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 ?-J\~AXL  
!边坡在强度折减系数F=2.6时结果分析 :tNH Cx  
Resume,'F2.6','db'                            !读入边坡在强度折减系数F=2.6时 w,D(zk$   
set,1,last                                    !读入后一个子步 M" %w9)@  
pldisp,1                                      !绘制边坡模型变形图 GtbIw  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 @TsOc0?-  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 entO"~*EX  
!边坡在强度折减系数F=2.8时结果分析 }F**!%4d  
Resume,'F2.8','db'                            !读入边坡在强度折减系数F=2.8时 Q;SMwCB0M  
set,1,last                                    !读入后一个子步 _aq3G9C_  
pldisp,1                                      !绘制边坡模型变形图 HJM-;C](  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 8L.Y0_x  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 Vhv<w O Ct  
!边坡在强度折减系数F=3.0时结果分析 ]M>mwnt+  
Resume,'F3.0','db'                            !读入边坡在强度折减系数F=3.0时 ]{Iy<  
set,1,last                                    !读入后一个子步 ><5tnBP|+L  
pldisp,1                                      !绘制边坡模型变形图 2{g~6U.  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 2,'m]`;GNr  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图

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

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

好东西,顶一下!

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

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