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

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

同求 呵呵 谁知道啊

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

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

我已经自己解决了，哎

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

(1)创建物理环境 <i_> y~v`  
/COM, Structural                             On{~St'V
 
/TITLE, Slope stability Analysis                 10C 2=  
/FILNAM,support,1                            
-A/ds1=;  
(2)建立模型 xab]q$n]k  
!进入前处理器 `y\*m]:  
/PREP7 9z5"y|$  
!定义单元类型 \8;Qv  
ET,1,PLANE82                               !定义边坡围岩单元 R*s* +I  
KEYOPT,1,3,2                                !求解类型为平面应变 I86e&"40  
!定义材料属性 xn(+G$m  
!边坡围岩1材料属性   X}oj_zsy;^  
MP,EX,1,3E10   Qs|OG  
MP,PRXY,1,0.25     jW-;4e*H=V  
MP,DENS,1,2500                             !定义塑性模型 SU?wFCGT%  
TB,DP,1                                     !定义数据表1 Y
]C;T  
TBDATA,1,9E5,42.8                           !给数据表输入数据 _Ry_K3K  
!边坡围岩2材料属性                         !采用弹性模型 #*!
+b  
MP,EX,2,3.2E10                              ！定义弹性模量 Vh"MKJ'R^  
MP,PRXY,2,0.24                              !定义泊松比 <,jAk
4  
MP,DENS,2,2700                             !定义密度 )46 0Ed  
!取强度折减系数F=1.2时边坡围岩材料属性 "'['(e+7  
MP,EX,3,3E10   6(>3P  
MP,PRXY,3,0.25     FE/$(7rM  
MP,DENS,3,2500 z>4D~HX  
TB,DP,3                                     !定义数据表3 +"]oc{W!  
TBDATA,1,7.5E5,37.7                         !给数据表输入数据 2HcsQ*H]G  
!取强度折减系数F=1.4时边坡围岩材料属性 u-m%=2  
MP,EX,4,3E10   [NYj.#,oR  
MP,PRXY,4,0.25                              !定义泊松比 blJIto'  
MP,DENS,4,2500 DdBxqkh  
TB,DP,4                                    !定义数据表4 hVMYB_<~  
TBDATA,1,6.4E5,33.5                        !给数据表输入数据 CN$I:o04C  
!取强度折减系数F=1.6时边坡围岩材料属性 B{s]juPG  
MP,EX,5,3E10   4\p$4Hs}  
MP,PRXY,5,0.25     h^}_YaT\  
MP,DENS,5,2500                             ！定义密度 yZq?
B  
TB,DP,5                                    !定义数据表5 :/+>e IE  
TBDATA,1,5.6E5,30                          !给数据表输入数据 *"1~b
Pl  
!取强度折减系数F=1.8时边坡围岩材料属性 %$Aqbd  
MP,EX,6,3E10                               ！定义弹性模量 p1+7<Y:  
MP,PRXY,6,0.25     9160L qY  
MP,DENS,6,2500 jvu N  
TB,DP,6                                    !定义数据表6 x+4vss  
TBDATA,1,5.0E5,27.2 E]NY (1  
!取强度折减系数F=2.0时边坡围岩材料属性 {5>3;.  
MP,EX,7,3E10   A:/}`  
MP,PRXY,7,0.25                             !定义泊松比 F&xv z2G  
MP,DENS,7,2500 Q~phGD3!~  
TB,DP,7                                    !定义数据表7 |Eb&}m:E$  
TBDATA,1,4.5E5,24.8 IWQ&6SDW$z  
!取强度折减系数F=2.2时边坡围岩材料属性 y =R aJm  
MP,EX,8,3E10   &GuF\wJ{7  
MP,PRXY,8,0.25     V|TA:&:7  
MP,DENS,8,2500 PNf&@  
TB,DP,8                                     !定义数据表8 \I;cZ>{u"}  
TBDATA,1,4.09E5,22.8 *2`:VFEV  
!取强度折减系数F=2.4时边坡围岩材料属性 8$ic~eJ  
MP,EX,9,3E10   [<yz)<<  
MP,PRXY,9,0.25     WHP;Neb6  
MP,DENS,9,2500 mbyih+amCr  
TB,DP,9                                      !定义数据表9 B%fU'  
TBDATA,1,3.6E5,21.1                          !给数据表输入数据 [m\,+lG?)j  
!取强度折减系数F=2.6时边坡围岩材料属性 (k
8Z=/N~  
MP,EX,10,3E10                                ！定义弹性模量   YZ< NP  
MP,PRXY,10,0.25                              !定义泊松比 ;u,
rtEMy;  
MP,DENS,10,2500                           ！定义密度 ehE-SrkU'  
TB,DP,10                                   !定义数据表10 _lP4}9p  
TBDATA,1,3.46E5,19.6                        !给数据表10输入数据 45)D+  
!取强度折减系数F=2.8时边坡围岩材料属性 JA<~xo[Q9  
MP,EX,11,3E10                              ！定义弹性模量 3w6&&R9  
MP,PRXY,11,0.25                             !定义泊松比 Atd1qJ  
MP,DENS,11,2500 iT.|vr1HG  
TB,DP,11                                   !定义数据表11 G,
]z(%  
TBDATA,1,3.2E5,18.2                         !给数据表11输入数据 1aq2aLx  
!取强度折减系数F=3.0时边坡围岩材料属性 8b7;\C~$p  
MP,EX,12,3E10                              ！定义弹性模量 eQ<xp
A  
MP,PRXY,12,0.25     M6_-f ;.  
MP,DENS,12,2500                            ！定义密度 =:gjz4}_8  
TB,DP,12                                   !定义数据表12 -D#5o,]3  
TBDATA,1,3.0E5,17.1                        !给数据表12输入数据 `E W
!-v)  
!建立几何模型 KzG_ <<  
!创建边坡线模型线 E67XPvo1+@  
K,1,,,,                                      ！ 创建关键点1 i:d`{kJ|[  
K,2,-800,,, $T),DUYO  
K,3,-800,-800,,   xQo~%wW,?  
K,4,-800,-1200,, b/nOdFO@  
K,5,1200,-1200,,   +*C^:^jA  
K,6,1200,-800,,   n'?4.tb  
K,7,1200,,, ^X&)'H  
K,8,1200,378,, UfKkgq#  
K,9,430,378,, 9#Y2`pT  
l,1,2                                        !连接关键点1、2生成直线L1 I8k  
l,2,3 @h#Xix7  
l,3,4 ~M*gsW$  
l,4,5 ?L\z}0#  
l,5,6 q`^3ov^</  
l,6,7 ufPCx|x~  
l,7,8 P>wZ~Hjk  
l,8,9 kwlC[G$j7  
l,9,1 | 3giZ{  
l,1,7 +i=p5d5  
l,3,6 4%7s259%  
!创建边坡面模型 +_~,86  
al,3,4,5,11                              !由线L3、L4、L5、L11生成一个面A1 nLJBq)i  
al,1,2,11,6,10                            !由线L1、L2、L1、L6、L10生成一个面A2 q5p e~  
al,7,8,9,10                              !由线L7、L8、L9、L10生成一个面A3 s@~/x5jwCs  
/PNUM, area, on                         !打开面号开关 X2^_~<I{,  
!划分网格生成有限元模型 TEWAZVE*  
! 划分边坡围岩2网格 mgVML&^  
!设置网格份数 SA+d&H}Fc  
lsel,s,,,3,5,2                               !选择线L3和L5 6V2j*J  
lesize,all,,,5                               !把所选择线分为5段 lR@& Z6lw  
lsel,s,,,4,11,7                              !选择线L4和L11 -!TcQzHUs  
lesize,all,,,26                              !把所选择线分为26段 [4>r6Hqxr  
mat,2                                    !给边坡围岩2赋予2号材料特性 wAh#   
type,1                                   !采用单元类型1 AS|Rd+.  
mshkey,0                                 !设定自由网格划分 h<% U["   
mshape,0                                 !设定四边形网格划分 KO&:06V{  
amesh,1                                  ！划分面积A1
y(v_-6b  
! 划分边坡围岩1网格 V#J"c8n  
!设置网格份数 7|Dn+=  
lsel,s,,,2,6,4                               !选择线L2和L6 }9CrFTbx;  
lsel,a,,,1                                  !附加选择线L1 tS<h8g_  
lesize,all,,,10                              !把所选择线分为10段 -:SIS`0s  
lsel,s,,,8,10,2                              !选择线L8和L10 7Z0/(V.-  
lesize,all,,,16                               !把所选择线分为16段 t',BI
  
lsel,s,,,7,9,2 QNArZ6UQ  
lesize,all,,,12 DA=1KaJ.  
mat,1                                     !给边坡围岩1赋予1号材料特性 O]{*(J/t  
type,1 lFfXWNb  
mshkey,0                                  !设定自由网格划分 gz61FW  
mshape,0                                  !设定四边形网格划分 GNU;jSh5  
amesh,2                                   !划分面积A2 ,I`_F,  
amesh,3                                   ！划分面积A3 cPuHLwwYf  
allsel "8HE^Po/pn  
!保存网格模型 T3z(k la  
save,Slope-grid.db 4AHL3@x  
(3)施加约束和荷载 6Z8l8:r-6  
!给边坡模型施加约束 ZEGd4_ux  
!边坡两侧施加X方向约束 Y<Q\d[3^F  
nsel,s,loc,x,-800                             !选择X=-800线上所有节点 G{o+R]Us  
nsel,a,loc,x,1200                             !选择X=1200线上所有节点 I4ilR$jg  
d,all,ux                                     !对所选择节点约束X方向位移 Q2c*.Y  
allsel S}7>
RHe  
!边坡底部施加约束 &{W^W8,%  
nsel,s,loc,y,-1200                           !选择Y=--1200线上所有节点 x=,8[W#XT  
d,all,uy                                   !对所选择节点约束X、Y方向位移 H]Wp%"L  
d,all,ux                                   U$'y_}V  
!施加重力加速度 >nry0 ;z0,  
acel,,9.8 VAKy^nR5j  
(4)求解 hYb9`0G"2  
/solu t5)J;0/  
!求解设置 V@Wcb$mgk  
antype,static                                  !设定为静力求解 A5XR3$
5P  
nsubst,100                                   !设定最大子步数为100 uT
GcQs}  
pred,on                                      !打开时间步长预测器 [H,u)8)  
nropt,full                                     !设定牛顿-拉普森选项 RrWNJ&o  
nlgeom,on                                    !打开大位移效果 ^~l  $&~  
lnsrch,on                                     !打开线性搜索 }-p,iTm  
outres,all,all                                  !输出所有项 2-v\3voN  
cnvtol,f,,0.005,2,0.5                            !力收敛准则设定 lIh[|]  
cnvtol,u,,0.05,2,1                              !位移收敛准则设定 T"!EK&  
!边坡在强度折减系数F=1时求解 D1Yh,P<CF\  
allsel -t4 [oB  
solve                                        !进行求解 gX,9Gh  
save,F1,db                                    !把F=1时求解结果保存 ow.j+<M  
!边坡在强度折减系数F=1.2时求解 5AOfp2O  
finish bx>i6 R2  
/solu :3?|VE F  
allsel o:UXPAj  
asel,s,area,,2,3,1                               !选择面积A2、A3 g$3>~D  
MPCHG,3,all                                 ! 把所选择单元材料号改为3 r7I B{}>-  
allsel X/+OF'po  
solve                                         !进行求解 gj+3y9  
save,F1.2,db                                   !把F=1.2时求解结果保存                                                                                   U)[ty@zyF  
!边坡在强度折减系数F=1.4时求解 c3k|G<C2  
finish j}RzXJ~t  
/solu XnXb&@Y  
allsel <xn;bp[  
asel,s,area,,2,3,1                               !选择面积A2、A3 \L[i9m|e  
MPCHG,4,all                                 !把所选择单元材料号改为4 s;{K!L@  
allsel U CY2]E  
solve                                        !进行求解 1S%}xsR0  
save,F1.4,db                                  !把F=1.4时求解结果保存                             `|<+  ?  
!边坡在强度折减系数F=1.6时求解 .mPg0  
finish 5 =Z!hQ}  
/solu &qRJceT(  
allsel Dg2uE8k  
asel,s,area,,2,3,1                               !选择面积A2、A3 Fe$/t(  
MPCHG,5,all                                 !把所选择单元材料号改为5 >n!ni(  
allsel )^ <3\e  
solve [xdj6W  
save,F1.6,db                                   !把F=1.6时求解结果保存                                                                                   Coa-8j*R7  
!边坡在强度折减系数F=1.8时求解 @2e2^8X7f  
finish 2c:f<>r0y  
/solu @ yxt($G  
allsel xj
q0D[  
asel,s,area,,2,3,1                                !选择面积A2、A3 )|]Z>>%t  
MPCHG,6,all                                  !把所选择单元材料号改为6 eb!_ie"D  
allsel ^ vbWRG~  
solve \0AiCMX[  
save,F1.8,db                                   !把F=1.8时求解结果保存   ~7SH4Cr  
!边坡在强度折减系数F=2.0时求解 2p:r`THvS5  
finish zk=\lp2  
/solu oJY[{-qW  
allsel ("Dv>&w9  
asel,s,area,,2,3,1                                !选择面积A2、A3  @Fx@5e  
MPCHG,7,all                                   !把所选择单元材料号改为7 #WufZ18#  
allsel !R:y'Y%j  
solve r4eUZ .8R  
save,F2.0,db                                   !把F=2.0时求解结果保存                             C6@t  
!边坡在强度折减系数F=2.2时求解 O9]\Q@M.  
finish `6Hf&u<  
/solu ._<ii2K'  
allsel 40K2uT{cq  
asel,s,area,,2,3,1                                !选择面积A2、A3 77Q}=80GU;  
MPCHG,8,all                                  !把所选择单元材料号改为8 'b[0ci:  
allsel Vxh.<b6&'  
solve h~miP7,c<u  
save,F2.2,db                                   !把F=2.2时求解结果保存                                                                                   K.DXJ UR  
!边坡在强度折减系数F=2.4时求解 \&q=@rJp(z  
finish /f3m)pT  
/solu LprM;Q_  
allsel 9 5 H?{  
asel,s,area,,2,3,1                                 !选择面积A2、A3 ^M Ey,  
MPCHG,9,all                                   !把所选择单元材料号改为9 'G#SLqZy  
allsel Z> 74.r  
solve 7BFN|S_l  
save,F2.4,db                                   !把F=2.4时求解结果保存   V5MO}  
!边坡在强度折减系数F=2.6时求解  5$Kf]ZP  
finish I@./${o  
/solu 6 =gp:I  
allsel . U/k<v<)6  
asel,s,area,,2,3,1                                 !选择面积A2、A3 =A
 n`D  
MPCHG,10,all                                   !把所选择单元材料号改为10 Tbh'_F6  
allsel (9GbG"   
solve &!1}`4$[T  
save,F2.6,db                                   !把F=2.6时求解结果保存                             y1C/v:;  
!边坡在强度折减系数F=2.8时求解 8(A{;9^g  
finish
y/R+$h(%  
/solu \P?A7vuhLs  
allsel 3N<&u   
asel,s,area,,2,3,1                                !选择面积A2、A3 Qpu3(`d<  
MPCHG,11,all                                 !把所选择单元材料号改为11 oI@9}*  
allsel uva\0q  
solve l{%a&/  
save,F2.8,db                                   !把F=2.8时求解结果保存                                                                                   dAxp ,):&J  
!边坡在强度折减系数F=3.0时求解 ,")F[%v  
finish <t}?$1  
/solu K`4GU[ul  
allsel T[2}p=<%  
asel,s,area,,2,3,1                             !选择面积A2、A3 )%mAZk-*;^  
MPCHG,12,all                              !把所选择单元材料号改为12 1xTTJyoq  
allsel fIyPFqf7w)  
solve 6tdI6  
save,F3.0,db                                 !把F=3.0时求解结果保存 2k+16/T  
/post1                                       ！进入后处理 ,- HIFbXx@  
!边坡在强度折减系数F=1时结果分析 Yx1 D)  
Resume,'F1','db'                              !读入边坡在强度折减系数F=1时 Q]\j>>  
set,1,last                                    !读入后一个子步 LK+felL  
pldisp,1                                      !绘制边坡模型变形图 {L9yhY
w  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 ;.s:X  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 ( u f5\}x  
!边坡在强度折减系数F=1.2时结果分析 h5-d;RKE  
Resume,'F1.2','db'                            !读入边坡在强度折减系数F=1.2时 ;|30QUYh  
set,1,last                                    !读入后一个子步 D#v?gPo4  
pldisp,1                                      !绘制边坡模型变形图 B>|5xpZM12  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 }1m_o@{3P  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 KXe ka  
!边坡在强度折减系数F=1.4时结果分析 Gah lS*W  
Resume,'F1.4','db'                            !读入边坡在强度折减系数F=1.4时 SDc" 4g`
 
set,1,last                                    !读入后一个子步 yWHne~!  
pldisp,1                                      !绘制边坡模型变形图 F]5\YYXO  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 re uYTH  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 ;r!\-]5$  
!边坡在强度折减系数F=1.6时结果分析 #F4X}  
Resume,'F1.6','db'                            !读入边坡在强度折减系数F=1.6时 u^=@DO
'  
set,1,last                                    !读入后一个子步 K-4tdC3  
pldisp,1                                      !绘制边坡模型变形图 v@_in(dk  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 kF-TG3  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 QymD-A"P  
!边坡在强度折减系数F=1.8时结果分析 Hsl{rN  
Resume,'F1.8','db'                            !读入边坡在强度折减系数F=1.8时 ragSy8M  
set,1,last                                    !读入后一个子步 =pp:j`B9(  
pldisp,1                                      !绘制边坡模型变形图 7bT /KLU  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 h{/ve`F>@  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 *V|zx#RN  
!边坡在强度折减系数F=2.0时结果分析 S=O$JP79  
Resume,'F2.0','db'                            !读入边坡在强度折减系数F=2.0时 _?b;0{93u  
set,1,last                                    !读入后一个子步 bp%S62Dj  
pldisp,1                                      !绘制边坡模型变形图 AbZ:(+@cP  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 :t>Q:mX(N  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 MPvWCPB  
!边坡在强度折减系数F=2.2时结果分析 h5
!d
  
Resume,'F2.2','db'                            !读入边坡在强度折减系数F=2.2时 Qa?aL  
set,1,last                                    !读入后一个子步 .)`-Hkxa  
pldisp,1                                      !绘制边坡模型变形图 (C\r&N  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 R3G@G  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 a: IwA9!L  
!边坡在强度折减系数F=2.4时结果分析 %#9P?COs&W  
Resume,'F2.4','db'                            !读入边坡在强度折减系数F=2.4时 ~4^p}{  
set,1,last                                    !读入后一个子步 \7 Mq $d  
pldisp,1                                      !绘制边坡模型变形图 S'@=3)  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 la702)N{  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 ^&m?qKN8  
!边坡在强度折减系数F=2.6时结果分析 ;[~^(. f  
Resume,'F2.6','db'                            !读入边坡在强度折减系数F=2.6时 mJ$Htyr  
set,1,last                                    !读入后一个子步 s3< F  
pldisp,1                                      !绘制边坡模型变形图 T6=-hA^A  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 ;
+/NjC1  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 oCA(FQ6  
!边坡在强度折减系数F=2.8时结果分析 x@Vt[}e  
Resume,'F2.8','db'                            !读入边坡在强度折减系数F=2.8时 .K7C-Xn=  
set,1,last                                    !读入后一个子步 V{$Sfmey  
pldisp,1                                      !绘制边坡模型变形图 M.[wKGX(  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 :gaETr  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 G>0)I  
!边坡在强度折减系数F=3.0时结果分析 KpGUq0d@  
Resume,'F3.0','db'                            !读入边坡在强度折减系数F=3.0时 %M6 c0d[9-  
set,1,last                                    !读入后一个子步 yoW> BX  
pldisp,1                                      !绘制边坡模型变形图 i[_|%'p  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 \n(
ROf^'  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图

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

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

好东西,顶一下!

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

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