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

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

同求 呵呵 谁知道啊

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

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

我已经自己解决了，哎

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

(1)创建物理环境 .m]}Ba}J$  
/COM, Structural                             BFW b0;+  
/TITLE, Slope stability Analysis                 !vf:mMo  
/FILNAM,support,1                             %N-aLw\  
(2)建立模型 |qMG@  
!进入前处理器 )7*'r@  
/PREP7 E\!X$  
!定义单元类型 cRf F!EV  
ET,1,PLANE82                               !定义边坡围岩单元 S&}7XjY  
KEYOPT,1,3,2                                !求解类型为平面应变 5\pS8<RJ;  
!定义材料属性 eY :"\c3  
!边坡围岩1材料属性   1+Gq<]@G  
MP,EX,1,3E10   P_5aHeiJ  
MP,PRXY,1,0.25     (O+d6oT=Z2  
MP,DENS,1,2500                             !定义塑性模型 hP+4{F*}-  
TB,DP,1                                     !定义数据表1 P#hRqETw  
TBDATA,1,9E5,42.8                           !给数据表输入数据 @~$d4K y<  
!边坡围岩2材料属性                         !采用弹性模型 OY{fxBb  
MP,EX,2,3.2E10                              ！定义弹性模量 {KHI(*r;  
MP,PRXY,2,0.24                              !定义泊松比 A6=Z2i0w>X  
MP,DENS,2,2700                             !定义密度 u*@R`,Y   
!取强度折减系数F=1.2时边坡围岩材料属性 -E|"?  
MP,EX,3,3E10   iKaX8c,zI  
MP,PRXY,3,0.25     k3$'K}=d  
MP,DENS,3,2500 X6 BIZ  
TB,DP,3                                     !定义数据表3 rQ*w3F?:  
TBDATA,1,7.5E5,37.7                         !给数据表输入数据 r[!~~yu/o  
!取强度折减系数F=1.4时边坡围岩材料属性 Nx.9)MjI  
MP,EX,4,3E10   nKGQU,C  
MP,PRXY,4,0.25                              !定义泊松比 f DXK<v)  
MP,DENS,4,2500 Hy3J2p9.  
TB,DP,4                                    !定义数据表4 vE)N6Ss  
TBDATA,1,6.4E5,33.5                        !给数据表输入数据 KhWy  
!取强度折减系数F=1.6时边坡围岩材料属性 6HCg<_j]  
MP,EX,5,3E10   t}R!i-D|HB  
MP,PRXY,5,0.25     YxM\qy{Vr  
MP,DENS,5,2500                             ！定义密度 =Oy&f:s  
TB,DP,5                                    !定义数据表5 3yRvs;nWS  
TBDATA,1,5.6E5,30                          !给数据表输入数据 [UXVL}tk  
!取强度折减系数F=1.8时边坡围岩材料属性 \YzKEYx+  
MP,EX,6,3E10                               ！定义弹性模量 @2"3RmYLo  
MP,PRXY,6,0.25     "J,|),Yd  
MP,DENS,6,2500 9R$$(zB 1;  
TB,DP,6                                    !定义数据表6 W\Pd:t  
TBDATA,1,5.0E5,27.2 q,ie)`  
!取强度折减系数F=2.0时边坡围岩材料属性 /4Wf\ Zu  
MP,EX,7,3E10   :4zPYG o  
MP,PRXY,7,0.25                             !定义泊松比 }T.>p#z  
MP,DENS,7,2500  ="]r{  
TB,DP,7                                    !定义数据表7 `>\ ~y1  
TBDATA,1,4.5E5,24.8  _:\rB  
!取强度折减系数F=2.2时边坡围岩材料属性 C2R"96M7
q  
MP,EX,8,3E10   a.<XJ\  
MP,PRXY,8,0.25     O;4S<N  
MP,DENS,8,2500 6xx.Z3v  
TB,DP,8                                     !定义数据表8 TFb7P/g  
TBDATA,1,4.09E5,22.8 !PoyM[Z"f  
!取强度折减系数F=2.4时边坡围岩材料属性 ]uFJ~:R  
MP,EX,9,3E10   >A D!)&c  
MP,PRXY,9,0.25     ^?fsJ  
MP,DENS,9,2500 e/Z{{FP%6  
TB,DP,9                                      !定义数据表9 *1v3x:pQ'  
TBDATA,1,3.6E5,21.1                          !给数据表输入数据 kq=Htbv7  
!取强度折减系数F=2.6时边坡围岩材料属性 \ &|xMw[  
MP,EX,10,3E10                                ！定义弹性模量   7Q{&L#;  
MP,PRXY,10,0.25                              !定义泊松比 HAE$Np|>a  
MP,DENS,10,2500                           ！定义密度 ^kK")+K  
TB,DP,10                                   !定义数据表10 ])y)]H#{  
TBDATA,1,3.46E5,19.6                        !给数据表10输入数据 R5b,/>^'A  
!取强度折减系数F=2.8时边坡围岩材料属性 8Q/cJ+&  
MP,EX,11,3E10                              ！定义弹性模量 dKPXs-5  
MP,PRXY,11,0.25                             !定义泊松比 9u=]D> kb  
MP,DENS,11,2500 O(evlci  
TB,DP,11                                   !定义数据表11 ZWyf.VJ  
TBDATA,1,3.2E5,18.2                         !给数据表11输入数据 nZj&Ma7R  
!取强度折减系数F=3.0时边坡围岩材料属性 Kc] GE#~g  
MP,EX,12,3E10                              ！定义弹性模量 0KExB{K  
MP,PRXY,12,0.25     TQt[he$O  
MP,DENS,12,2500                            ！定义密度 -d!84_d9  
TB,DP,12                                   !定义数据表12 lz0dt<8eP  
TBDATA,1,3.0E5,17.1                        !给数据表12输入数据 w@P86'< v  
!建立几何模型 umPN=0u6  
!创建边坡线模型线 Ld:-S,2  
K,1,,,,                                      ！ 创建关键点1 pPsTgGai  
K,2,-800,,, MIr[_  
K,3,-800,-800,,   }:?_/$};  
K,4,-800,-1200,, CiU^U|~'L  
K,5,1200,-1200,,   Q.4+"JoG  
K,6,1200,-800,,   r8y,$Mv<)0  
K,7,1200,,, A5fzyG
 
K,8,1200,378,, agoMsxI9  
K,9,430,378,, p/ZgzHyF  
l,1,2                                        !连接关键点1、2生成直线L1 siuDg,uqK5  
l,2,3 "OP$n-*@%  
l,3,4 (8S+-k?  
l,4,5 T;k
h+i  
l,5,6 qoQ,3&<  
l,6,7 !+hw8@A  
l,7,8
qCk`398W  
l,8,9 |.~2C14[  
l,9,1 0a?[@ -Sz  
l,1,7 (F R  
l,3,6 ODCN~7-@  
!创建边坡面模型 4[r:DM|8  
al,3,4,5,11                              !由线L3、L4、L5、L11生成一个面A1 ywjD.od"v  
al,1,2,11,6,10                            !由线L1、L2、L1、L6、L10生成一个面A2 Ov:U3P?%  
al,7,8,9,10                              !由线L7、L8、L9、L10生成一个面A3 (ia+N/$u  
/PNUM, area, on                         !打开面号开关 U&+lw=  
!划分网格生成有限元模型 X=Y>9  
! 划分边坡围岩2网格 YTe8C9eO  
!设置网格份数 $xcU*?=K  
lsel,s,,,3,5,2                               !选择线L3和L5 &fxyY(  
lesize,all,,,5                               !把所选择线分为5段 KmF+3g~#s  
lsel,s,,,4,11,7                              !选择线L4和L11 BYI13jMH+Y  
lesize,all,,,26                              !把所选择线分为26段 :$^sI"hO  
mat,2                                    !给边坡围岩2赋予2号材料特性 SsL>K*t5  
type,1                                   !采用单元类型1 UD'e%IVw  
mshkey,0                                 !设定自由网格划分 q>a/',m  
mshape,0                                 !设定四边形网格划分 P^/e!%UgC  
amesh,1                                  ！划分面积A1 t2-bw6U  
! 划分边坡围岩1网格 P

-N+  
!设置网格份数 Y,%G5X@S<  
lsel,s,,,2,6,4                               !选择线L2和L6 ?P2d 9b  
lsel,a,,,1                                  !附加选择线L1 }PzYt~Z`@  
lesize,all,,,10                              !把所选择线分为10段 N;av  
lsel,s,,,8,10,2                              !选择线L8和L10 BjUz"69  
lesize,all,,,16                               !把所选择线分为16段 }'TZ)=t{J  
lsel,s,,,7,9,2 75H5{#)  
lesize,all,,,12 ZnB|vfL?  
mat,1                                     !给边坡围岩1赋予1号材料特性 WB|SXto%4D  
type,1 #w]:<R^  
mshkey,0                                  !设定自由网格划分 gl6*bB=  
mshape,0                                  !设定四边形网格划分 ;Q8rAsf9  
amesh,2                                   !划分面积A2 TCO^9RP<  
amesh,3                                   ！划分面积A3 JMYM}G  
allsel }qdGS<{  
!保存网格模型 }"9jCxXL  
save,Slope-grid.db =kWm9W<^  
(3)施加约束和荷载 !bD`2m[Q  
!给边坡模型施加约束 A?OaP  
!边坡两侧施加X方向约束 YR/rN,  
nsel,s,loc,x,-800                             !选择X=-800线上所有节点 a@WSIcX*W  
nsel,a,loc,x,1200                             !选择X=1200线上所有节点 KbF,jm5  
d,all,ux                                     !对所选择节点约束X方向位移 E2{SKIUm  
allsel faaFmEC  
!边坡底部施加约束 JS1''^G&.  
nsel,s,loc,y,-1200                           !选择Y=--1200线上所有节点 VFwp .1oa!  
d,all,uy                                   !对所选择节点约束X、Y方向位移 gnSb)!i>z  
d,all,ux                                   /7<l`RSr  
!施加重力加速度 s5ddGiZnBT  
acel,,9.8 cw{[% 7  
(4)求解 GKKf#r74  
/solu q6v%HF-q4  
!求解设置 j_*#"}Lcp  
antype,static                                  !设定为静力求解 {__"Z<  
nsubst,100                                   !设定最大子步数为100 VDq?,4Kb  
pred,on                                      !打开时间步长预测器 'PrrP3lO_~  
nropt,full                                     !设定牛顿-拉普森选项 H{CG/+x  
nlgeom,on                                    !打开大位移效果 d^W
EfH  
lnsrch,on                                     !打开线性搜索 !'^gqaF+  
outres,all,all                                  !输出所有项 9n"
D/NZB  
cnvtol,f,,0.005,2,0.5                            !力收敛准则设定 W[Ro)  
cnvtol,u,,0.05,2,1                              !位移收敛准则设定 EBN'u&zX  
!边坡在强度折减系数F=1时求解 AA$-Lx(UJk  
allsel rOu7r4  
solve                                        !进行求解 jn9KQe\3  
save,F1,db                                    !把F=1时求解结果保存 Fh4kd>1D  
!边坡在强度折减系数F=1.2时求解 QvH=<$  
finish Ojf.D6nY  
/solu TPrqb  
allsel kN*,3)T;}  
asel,s,area,,2,3,1                               !选择面积A2、A3 Rj9YAW$  
MPCHG,3,all                                 ! 把所选择单元材料号改为3 ;X ]+r$_  
allsel {WJ+6!v  
solve                                         !进行求解 @exeHcW61  
save,F1.2,db                                   !把F=1.2时求解结果保存                                                                                   V!F# ek:  
!边坡在强度折减系数F=1.4时求解 ['~B&  
finish |Psi?'4  
/solu #?Z>o16,u  
allsel Oozt&* F  
asel,s,area,,2,3,1                               !选择面积A2、A3 ShdE!q7  
MPCHG,4,all                                 !把所选择单元材料号改为4 =r=YV-D.  
allsel e:E:"elr]  
solve                                        !进行求解 Gvk)H$ni  
save,F1.4,db                                  !把F=1.4时求解结果保存                             c_e2'K:  
!边坡在强度折减系数F=1.6时求解 YzqUOMAt"V  
finish w]hs1vch  
/solu ~5[#c27E9  
allsel SM2
QF
 
asel,s,area,,2,3,1                               !选择面积A2、A3 @Q;%hb  
MPCHG,5,all                                 !把所选择单元材料号改为5 F(J6 XnQ  
allsel )DS|mM)  
solve z %Ty;  
save,F1.6,db                                   !把F=1.6时求解结果保存                                                                                   v"VpE`
z1#  
!边坡在强度折减系数F=1.8时求解 K]{Y >w  
finish Lj"@JF;c  
/solu *UVo>;  
allsel )_eEM1  
asel,s,area,,2,3,1                                !选择面积A2、A3 owDp?Sy}E  
MPCHG,6,all                                  !把所选择单元材料号改为6 n 7Mab  
allsel nm)H\i  
solve K!c "g,S  
save,F1.8,db                                   !把F=1.8时求解结果保存   1PD{m{  
!边坡在强度折减系数F=2.0时求解 }&+b\RE  
finish K}e%E&|>  
/solu uHj"nd13  
allsel .v?x>iV  
asel,s,area,,2,3,1                                !选择面积A2、A3 N8!V%i?  
MPCHG,7,all                                   !把所选择单元材料号改为7 A}Q6DHh26  
allsel s9R#rwI
c  
solve tI42]:z  
save,F2.0,db                                   !把F=2.0时求解结果保存                             f;AQw_{  
!边坡在强度折减系数F=2.2时求解 _Mi`]VSq9  
finish {t<E*5N]a  
/solu TI^W=5W@@  
allsel dux.Z9X?  
asel,s,area,,2,3,1                                !选择面积A2、A3 |$*1!pL-QP  
MPCHG,8,all                                  !把所选择单元材料号改为8 w$Zi'+&*  
allsel P%VEJ5,]b  
solve a_jw4"Sb  
save,F2.2,db                                   !把F=2.2时求解结果保存                                                                                   T?vM\o%i3  
!边坡在强度折减系数F=2.4时求解 [D
r'  
finish )ry7a .39b  
/solu T_\Nvzb}  
allsel
6_Ps*Ed  
asel,s,area,,2,3,1                                 !选择面积A2、A3 =<yMB d\  
MPCHG,9,all                                   !把所选择单元材料号改为9 eh\_;2P  
allsel LqNt.d @  
solve O;zq(/,-l  
save,F2.4,db                                   !把F=2.4时求解结果保存   bSa%?laS  
!边坡在强度折减系数F=2.6时求解 k }amSsE  
finish \9 ,a"g  
/solu _}bs0 kIz  
allsel WC&V9Yk  
asel,s,area,,2,3,1                                 !选择面积A2、A3 i}u,_ }  
MPCHG,10,all                                   !把所选择单元材料号改为10 R[#Np`z  
allsel @477|LO
  
solve O`?qnNmc;  
save,F2.6,db                                   !把F=2.6时求解结果保存                             tQR qQ  
!边坡在强度折减系数F=2.8时求解 F,dPmR  
finish zR@4Z>6   
/solu .rPg  
allsel u[:-^H  
asel,s,area,,2,3,1                                !选择面积A2、A3 6!39t  
MPCHG,11,all                                 !把所选择单元材料号改为11 /penB[1i  
allsel 6Qc *:(GE  
solve qjr:(x/  
save,F2.8,db                                   !把F=2.8时求解结果保存                                                                                  
;AH8/M B9  
!边坡在强度折减系数F=3.0时求解 > J
.q3  
finish $^T
xLv  
/solu Z=+Tw!wR>  
allsel i'cGB5-j  
asel,s,area,,2,3,1                             !选择面积A2、A3 h
5)4Z^n  
MPCHG,12,all                              !把所选择单元材料号改为12 #]rfKHW9  
allsel &0b\E73  
solve fw&cv9X(IU  
save,F3.0,db                                 !把F=3.0时求解结果保存 R$Rub/b6  
/post1                                       ！进入后处理 +X%pUe  
!边坡在强度折减系数F=1时结果分析 po! [Nd&"  
Resume,'F1','db'                              !读入边坡在强度折减系数F=1时 W{$J)iQ  
set,1,last                                    !读入后一个子步 xcQD]"  
pldisp,1                                      !绘制边坡模型变形图 R@df~  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 <)$e*HrI  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 +B '<0  
!边坡在强度折减系数F=1.2时结果分析 Vg^yjP{sv  
Resume,'F1.2','db'                            !读入边坡在强度折减系数F=1.2时 Leu6kPk  
set,1,last                                    !读入后一个子步 7VIfRN{5n  
pldisp,1                                      !绘制边坡模型变形图 \b
;z$P\+*  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 eK
[9wEdn  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 x%yzhIRR  
!边坡在强度折减系数F=1.4时结果分析 6vfut$)[{  
Resume,'F1.4','db'                            !读入边坡在强度折减系数F=1.4时 "8$Muwm  
set,1,last                                    !读入后一个子步 6fm oIK{  
pldisp,1                                      !绘制边坡模型变形图 L9]y~[R:  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 j6d{r\!$4  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 BhUGMK  
!边坡在强度折减系数F=1.6时结果分析 &
s8vmUt  
Resume,'F1.6','db'                            !读入边坡在强度折减系数F=1.6时 A I v  
set,1,last                                    !读入后一个子步 r,<p#4(>_  
pldisp,1                                      !绘制边坡模型变形图 ;uho.)%N`F  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 oe*fgk/o9  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 $ghlrV;:ct  
!边坡在强度折减系数F=1.8时结果分析 m;J'y2h =$  
Resume,'F1.8','db'                            !读入边坡在强度折减系数F=1.8时 Bdcs}Ga  
set,1,last                                    !读入后一个子步 ,`ba?O?*G  
pldisp,1                                      !绘制边坡模型变形图 2"QcjFW%  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 '5lwlF  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 v]'\]U^  
!边坡在强度折减系数F=2.0时结果分析 M_O)w^ '  
Resume,'F2.0','db'                            !读入边坡在强度折减系数F=2.0时 U UhlKV|5  
set,1,last                                    !读入后一个子步 sk5B} -  
pldisp,1                                      !绘制边坡模型变形图 pwvmb\  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 V$_.&S?(Y  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 L@S1C=-/  
!边坡在强度折减系数F=2.2时结果分析 }ww`Y&#  
Resume,'F2.2','db'                            !读入边坡在强度折减系数F=2.2时 Elt"tJ  
set,1,last                                    !读入后一个子步 ,*sKr)9)  
pldisp,1                                      !绘制边坡模型变形图 b%h.>ij?  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 _4 YT2k  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 EeQ2\'t  
!边坡在强度折减系数F=2.4时结果分析 Eb8~i_B-  
Resume,'F2.4','db'                            !读入边坡在强度折减系数F=2.4时 OGcW]i  
set,1,last                                    !读入后一个子步 !iL6/  
pldisp,1                                      !绘制边坡模型变形图 ( m7qc  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 Ri3*au/
Q  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 Q yhu=_&  
!边坡在强度折减系数F=2.6时结果分析 g'9~T8i& ^  
Resume,'F2.6','db'                            !读入边坡在强度折减系数F=2.6时 `NQ;|!  
set,1,last                                    !读入后一个子步 FP")$ ,=s  
pldisp,1                                      !绘制边坡模型变形图 <r)5jf  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 n0)y|B#  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 [[O4_)?el  
!边坡在强度折减系数F=2.8时结果分析 k_nQmU>  
Resume,'F2.8','db'                            !读入边坡在强度折减系数F=2.8时 braI MIQ`  
set,1,last                                    !读入后一个子步 z%nplG'~|  
pldisp,1                                      !绘制边坡模型变形图 L=5Y^f'aU  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 od^o9(.W^  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 4j(*%da  
!边坡在强度折减系数F=3.0时结果分析 vcZ"4%w  
Resume,'F3.0','db'                            !读入边坡在强度折减系数F=3.0时 )1g\v8XT  
set,1,last                                    !读入后一个子步 Y}h&dAr  
pldisp,1                                      !绘制边坡模型变形图 @cQ |`  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 6%V#_]  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图

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

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

好东西,顶一下!

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

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