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

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

同求 呵呵 谁知道啊

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

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

我已经自己解决了，哎

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

(1)创建物理环境 .|T2\M  
/COM, Structural                             *Y8XP8u/  
/TITLE, Slope stability Analysis                 >HE,'  
/FILNAM,support,1                             'eqiYY|  
(2)建立模型 tY{; U#9  
!进入前处理器 i4h

JE  
/PREP7 ,/~[S  
!定义单元类型 g93I+  
ET,1,PLANE82                               !定义边坡围岩单元 'b8R#R\P  
KEYOPT,1,3,2                                !求解类型为平面应变 66oK3%[  
!定义材料属性 ?sV[MsOsC  
!边坡围岩1材料属性   O x{Q.l  
MP,EX,1,3E10   Kn']n91m  
MP,PRXY,1,0.25     |kId8WtA  
MP,DENS,1,2500                             !定义塑性模型 0r1g$mKb  
TB,DP,1                                     !定义数据表1 joe9.{  
TBDATA,1,9E5,42.8                           !给数据表输入数据 -Bj.hx*  
!边坡围岩2材料属性                         !采用弹性模型 2*+3RrJ  
MP,EX,2,3.2E10                              ！定义弹性模量 W&<g} N+  
MP,PRXY,2,0.24                              !定义泊松比 g<f
P:/  
MP,DENS,2,2700                             !定义密度 fCLcU@3W?  
!取强度折减系数F=1.2时边坡围岩材料属性 p?Z(rCp  
MP,EX,3,3E10   Gu2_dT  
MP,PRXY,3,0.25     3f_i1|>)'  
MP,DENS,3,2500 SV&kWbS  
TB,DP,3                                     !定义数据表3 S,lxM,DL&  
TBDATA,1,7.5E5,37.7                         !给数据表输入数据 !d\t:0;  
!取强度折减系数F=1.4时边坡围岩材料属性 doLkrEm&
 
MP,EX,4,3E10   ^rL,&rk  
MP,PRXY,4,0.25                              !定义泊松比 #-,g&)`]  
MP,DENS,4,2500 v#zPH5xo  
TB,DP,4                                    !定义数据表4 %>i@F=O2<  
TBDATA,1,6.4E5,33.5                        !给数据表输入数据 O'fk&&l  
!取强度折减系数F=1.6时边坡围岩材料属性 1XG$ z@NN  
MP,EX,5,3E10   DYX-5~;!  
MP,PRXY,5,0.25     /v5qyR7an  
MP,DENS,5,2500                             ！定义密度 /E)9v$!  
TB,DP,5                                    !定义数据表5 =KV@&Y^x4  
TBDATA,1,5.6E5,30                          !给数据表输入数据 /4yOs@#  
!取强度折减系数F=1.8时边坡围岩材料属性 ?~!tM}X0:3  
MP,EX,6,3E10                               ！定义弹性模量 0[.3Es:_  
MP,PRXY,6,0.25     H\ 3M  
MP,DENS,6,2500 8GY.){d!l  
TB,DP,6                                    !定义数据表6 _HwpPRVP/  
TBDATA,1,5.0E5,27.2 pP3U,n   
!取强度折减系数F=2.0时边坡围岩材料属性 ]22C)<  
MP,EX,7,3E10   iu+3,]7Fm  
MP,PRXY,7,0.25                             !定义泊松比 (/'h4KS@  
MP,DENS,7,2500 HCJ;&C73&  
TB,DP,7                                    !定义数据表7  KZ]r8  
TBDATA,1,4.5E5,24.8 p:B ]Ft  
!取强度折减系数F=2.2时边坡围岩材料属性 h9)RJSF4  
MP,EX,8,3E10   ~u!gUJ:  
MP,PRXY,8,0.25     F@9Y\. ,  
MP,DENS,8,2500 6{Ks`Af  
TB,DP,8                                     !定义数据表8 pqJ)G;%9  
TBDATA,1,4.09E5,22.8 +Z ><  
!取强度折减系数F=2.4时边坡围岩材料属性 jM*AL X  
MP,EX,9,3E10   9J?W '8s5  
MP,PRXY,9,0.25     |Td_S|:d  
MP,DENS,9,2500 PCtkjd  
TB,DP,9                                      !定义数据表9 ; VBpp<  
TBDATA,1,3.6E5,21.1                          !给数据表输入数据 w5,6$#  
!取强度折减系数F=2.6时边坡围岩材料属性 m`'=)x|  
MP,EX,10,3E10                                ！定义弹性模量   RYt6=R+f  
MP,PRXY,10,0.25                              !定义泊松比 U M@naU  
MP,DENS,10,2500                           ！定义密度 J=):+F=  
TB,DP,10                                   !定义数据表10 K${}r0   
TBDATA,1,3.46E5,19.6                        !给数据表10输入数据 /M:H9Z8!  
!取强度折减系数F=2.8时边坡围岩材料属性 9X3yp:>V  
MP,EX,11,3E10                              ！定义弹性模量 V7P6zAJy  
MP,PRXY,11,0.25                             !定义泊松比 \
4aKLr  
MP,DENS,11,2500 =:4?>2)  
TB,DP,11                                   !定义数据表11 Y:wF5pp;  
TBDATA,1,3.2E5,18.2                         !给数据表11输入数据 N*f^Z#B]  
!取强度折减系数F=3.0时边坡围岩材料属性 M2dm
G<  
MP,EX,12,3E10                              ！定义弹性模量 )i @1XH"D  
MP,PRXY,12,0.25     q?yMa9ZZky  
MP,DENS,12,2500                            ！定义密度 &RWM<6JP  
TB,DP,12                                   !定义数据表12 Lj iI+NJ  
TBDATA,1,3.0E5,17.1                        !给数据表12输入数据 IK3qE!,&U  
!建立几何模型 .?f:Nb.O  
!创建边坡线模型线 @.k5MOn  
K,1,,,,                                      ！ 创建关键点1 )`<6taKx@n  
K,2,-800,,, 0F^]A"kF  
K,3,-800,-800,,   @YCv  
K,4,-800,-1200,, aRX  
K,5,1200,-1200,,   zHV|-R  
K,6,1200,-800,,   3x![8 x  
K,7,1200,,, Hr6wgYPi  
K,8,1200,378,, )6G"*  
K,9,430,378,, H"O$&  
l,1,2                                        !连接关键点1、2生成直线L1 t~ -J %$  
l,2,3 '|&,E#`  
l,3,4 ^PC\E}  
l,4,5 HM):"  
l,5,6 $:e)$Xnn-  
l,6,7 y<|)'(  
l,7,8 ?s%v 3T  
l,8,9 x:qr\Rz  
l,9,1 dsK/6yu  
l,1,7 H-Pq!9[DB  
l,3,6 mY`@'  
!创建边坡面模型 AQe!Sqg'  
al,3,4,5,11                              !由线L3、L4、L5、L11生成一个面A1 3q
"7K  
al,1,2,11,6,10                            !由线L1、L2、L1、L6、L10生成一个面A2 9Fk4|+OJ  
al,7,8,9,10                              !由线L7、L8、L9、L10生成一个面A3 b{BaQ>.(`  
/PNUM, area, on                         !打开面号开关 %lV@:"G  
!划分网格生成有限元模型 ` VwN!B:  
! 划分边坡围岩2网格 [7RheXO<  
!设置网格份数 Ae6("Oid  
lsel,s,,,3,5,2                               !选择线L3和L5 b"t")U==  
lesize,all,,,5                               !把所选择线分为5段 ?ZaD=nh$mK  
lsel,s,,,4,11,7                              !选择线L4和L11 \BUqDd!  
lesize,all,,,26                              !把所选择线分为26段 Wk}D]o0^@  
mat,2                                    !给边坡围岩2赋予2号材料特性 Mm.Ql  
type,1                                   !采用单元类型1 O] H=s  
mshkey,0                                 !设定自由网格划分 %]#VdS|N
  
mshape,0                                 !设定四边形网格划分 b~j~  
amesh,1                                  ！划分面积A1 ) oxIzF  
! 划分边坡围岩1网格 847 R   
!设置网格份数 QNb>rLj52  
lsel,s,,,2,6,4                               !选择线L2和L6 %[XY67A3I  
lsel,a,,,1                                  !附加选择线L1 dhW<p5  
lesize,all,,,10                              !把所选择线分为10段 P% Q@9kO>  
lsel,s,,,8,10,2                              !选择线L8和L10 dnwdFsf  
lesize,all,,,16                               !把所选择线分为16段 .liyC~YW  
lsel,s,,,7,9,2 O4E(R?wd  
lesize,all,,,12 *="m3:c'J  
mat,1                                     !给边坡围岩1赋予1号材料特性 WO
@H*  
type,1 ".E5t@ }?m  
mshkey,0                                  !设定自由网格划分 9k}<Fz"^.  
mshape,0                                  !设定四边形网格划分 ) \
4 |  
amesh,2                                   !划分面积A2 dgslUg9z3g  
amesh,3                                   ！划分面积A3 jXWNHIl)@  
allsel Mb/L~gd"  
!保存网格模型 kxh 5}eB  
save,Slope-grid.db fe';b[q)#  
(3)施加约束和荷载 "7d.i(vw  
!给边坡模型施加约束 o!utZmk$  
!边坡两侧施加X方向约束 O~6Q;qP  
nsel,s,loc,x,-800                             !选择X=-800线上所有节点 7 5|pp  
nsel,a,loc,x,1200                             !选择X=1200线上所有节点 UQ~4c,  
d,all,ux                                     !对所选择节点约束X方向位移 *0~M
 
allsel AFm,CINa  
!边坡底部施加约束 UW/N MjK  
nsel,s,loc,y,-1200                           !选择Y=--1200线上所有节点 E+z18Lf?  
d,all,uy                                   !对所选择节点约束X、Y方向位移 =-qf;5[|  
d,all,ux                                   =53bLzr  
!施加重力加速度 q`[K3p   
acel,,9.8 C3)|<E  
(4)求解 y:,9I`aW  
/solu /VO^5Dnb  
!求解设置 8?1o<8hV  
antype,static                                  !设定为静力求解 ESIzGaM  
nsubst,100                                   !设定最大子步数为100 Mn@$;\:  
pred,on                                      !打开时间步长预测器 5U
~OP  
nropt,full                                     !设定牛顿-拉普森选项 5z~\5x  
nlgeom,on                                    !打开大位移效果 lbS?/f  
lnsrch,on                                     !打开线性搜索 ]26m
B  
outres,all,all                                  !输出所有项 H} 6CKP}  
cnvtol,f,,0.005,2,0.5                            !力收敛准则设定 wyzOcx>M  
cnvtol,u,,0.05,2,1                              !位移收敛准则设定 +2Ql~w@$^l  
!边坡在强度折减系数F=1时求解 |!Fk2Je,  
allsel waCboK'  
solve                                        !进行求解 /[iG5~G  
save,F1,db                                    !把F=1时求解结果保存 q & b5g !  
!边坡在强度折减系数F=1.2时求解 69/?7r  
finish TP{Gt.e  
/solu G'9{a'  
allsel T(V8;!  
asel,s,area,,2,3,1                               !选择面积A2、A3 JOHRmfqR  
MPCHG,3,all                                 ! 把所选择单元材料号改为3 .'/l'>  
allsel (]XbPW  
solve                                         !进行求解 b_=8!Q.:  
save,F1.2,db                                   !把F=1.2时求解结果保存                                                                                   KmL$M  
!边坡在强度折减系数F=1.4时求解 sPy2/7Wqd  
finish 87<9V.s2  
/solu xs%LRF#u  
allsel 5oJ Dux }  
asel,s,area,,2,3,1                               !选择面积A2、A3 7]s%rya  
MPCHG,4,all                                 !把所选择单元材料号改为4 {5-zyE  
allsel !}5*?k g  
solve                                        !进行求解 [O_^MA,z  
save,F1.4,db                                  !把F=1.4时求解结果保存                             :+:6_x  
!边坡在强度折减系数F=1.6时求解 w;+ br  
finish YjLe(+WQ  
/solu AW/wI6[T  
allsel q@kOTkHv)  
asel,s,area,,2,3,1                               !选择面积A2、A3 x'qWM/  
MPCHG,5,all                                 !把所选择单元材料号改为5 sT`^ljp4  
allsel sAYV)w3u"  
solve &K *X)DAs  
save,F1.6,db                                   !把F=1.6时求解结果保存                                                                                   g4wZvra6%)  
!边坡在强度折减系数F=1.8时求解 LN^UC$[tk  
finish ik5|,#}m&  
/solu {zP#woz2Q  
allsel LwOJ|jA(,  
asel,s,area,,2,3,1                                !选择面积A2、A3 0[)VO[  
MPCHG,6,all                                  !把所选择单元材料号改为6 > :Ze4}(  
allsel O:^'x*}  
solve `hzrfum4  
save,F1.8,db                                   !把F=1.8时求解结果保存   j#VIHCzlr  
!边坡在强度折减系数F=2.0时求解
?PH/?QP  
finish 8*s7m   
/solu VFSz-<L  
allsel %iJ|H(P  
asel,s,area,,2,3,1                                !选择面积A2、A3 5m7b\Mak  
MPCHG,7,all                                   !把所选择单元材料号改为7 *,lh:  
allsel $O9Xx  
solve Us-A+)r*!  
save,F2.0,db                                   !把F=2.0时求解结果保存                             W2eAhz&  
!边坡在强度折减系数F=2.2时求解 Q]rqD83((  
finish EZib1g&:R/  
/solu )o jDRJ&  
allsel 7~b!4x|Z  
asel,s,area,,2,3,1                                !选择面积A2、A3 hwVAXsF~  
MPCHG,8,all                                  !把所选择单元材料号改为8 ru{f]|  
allsel 2x
n<E>]  
solve *d>vR1  
save,F2.2,db                                   !把F=2.2时求解结果保存                                                                                   ^i'y6J  
!边坡在强度折减系数F=2.4时求解 J>%t<xYf4  
finish 94{)"w]  
/solu aD ESr?  
allsel XV=S)  
asel,s,area,,2,3,1                                 !选择面积A2、A3 s4$X  
MPCHG,9,all                                   !把所选择单元材料号改为9 M
3(k'q7&:  
allsel /.$L"u  
solve T4r5s  
save,F2.4,db                                   !把F=2.4时求解结果保存   6>,# 6{?jl  
!边坡在强度折减系数F=2.6时求解 NR4Jn?l{  
finish C),7- ?  
/solu v"RiPHLT  
allsel a4&:@`=  
asel,s,area,,2,3,1                                 !选择面积A2、A3 k|FSz#Y  
MPCHG,10,all                                   !把所选择单元材料号改为10 FZXyfZw!|  
allsel Jq .L:>x  
solve OJ/SYZ.r  
save,F2.6,db                                   !把F=2.6时求解结果保存                             0^#DNq*NQ  
!边坡在强度折减系数F=2.8时求解 {155b0  
finish p7C!G1+z  
/solu n+;vjVS%  
allsel CCqT tp  
asel,s,area,,2,3,1                                !选择面积A2、A3 P+Z\3re  
MPCHG,11,all                                 !把所选择单元材料号改为11 O-PdM`mqW  
allsel M&y5AB0  
solve [bjN f2  
save,F2.8,db                                   !把F=2.8时求解结果保存                                                                                   2*u.3,aW  
!边坡在强度折减系数F=3.0时求解 <'&F;5F3V  
finish RtC'v";6  
/solu hS:jBp,  
allsel [M:S`{SbY  
asel,s,area,,2,3,1                             !选择面积A2、A3 +.@c{5J<  
MPCHG,12,all                              !把所选择单元材料号改为12 :c7CiP  
allsel V"#Jk!k9k  
solve ia4k:\  
save,F3.0,db                                 !把F=3.0时求解结果保存 ,w2WS\`%  
/post1                                       ！进入后处理 =nU/ [T.  
!边坡在强度折减系数F=1时结果分析 b/<mRQ{  
Resume,'F1','db'                              !读入边坡在强度折减系数F=1时 h/<=u9J  
set,1,last                                    !读入后一个子步 3SNL5  
pldisp,1                                      !绘制边坡模型变形图 R#qI(V  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 \84v-VK  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 1b3(  
!边坡在强度折减系数F=1.2时结果分析 .X
`C^z]+  
Resume,'F1.2','db'                            !读入边坡在强度折减系数F=1.2时 OVR?*"N_  
set,1,last                                    !读入后一个子步 |s=`w8p  
pldisp,1                                      !绘制边坡模型变形图 mW4%2fD[  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 4x,hj  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 >SDpuG&>  
!边坡在强度折减系数F=1.4时结果分析 T'1gy}  
Resume,'F1.4','db'                            !读入边坡在强度折减系数F=1.4时 f^
9&WT  
set,1,last                                    !读入后一个子步 `FJ|W6%  
pldisp,1                                      !绘制边坡模型变形图 XoItV  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 {Q~7M$  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 VVuR+=.&  
!边坡在强度折减系数F=1.6时结果分析 vZkXt!%)  
Resume,'F1.6','db'                            !读入边坡在强度折减系数F=1.6时 i8~r  
set,1,last                                    !读入后一个子步 |nY~ZVTt/  
pldisp,1                                      !绘制边坡模型变形图 mxJXL":|  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 &U"X$aFc  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 G{b:i8}l  
!边坡在强度折减系数F=1.8时结果分析 B4yh3cf  
Resume,'F1.8','db'                            !读入边坡在强度折减系数F=1.8时 Y"!uU.=xJ  
set,1,last                                    !读入后一个子步 N:x0w+Ca  
pldisp,1                                      !绘制边坡模型变形图 7petHi  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 S&?7K-F>_o  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 D.*>;5:0'  
!边坡在强度折减系数F=2.0时结果分析 hz\WZ^  
Resume,'F2.0','db'                            !读入边坡在强度折减系数F=2.0时 eko]H!Ov(  
set,1,last                                    !读入后一个子步 l67KJ  
pldisp,1                                      !绘制边坡模型变形图 maC>LBa2/  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 %0XvJF)s  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 z[9UQU~x?  
!边坡在强度折减系数F=2.2时结果分析 \c7>:DH  
Resume,'F2.2','db'                            !读入边坡在强度折减系数F=2.2时 6E|S  
set,1,last                                    !读入后一个子步 tln1eN((q  
pldisp,1                                      !绘制边坡模型变形图 *)>do L  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 IU!Ht>  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 _%ZP{5D>  
!边坡在强度折减系数F=2.4时结果分析 kus}WJ  
Resume,'F2.4','db'                            !读入边坡在强度折减系数F=2.4时 V1utUGJV  
set,1,last                                    !读入后一个子步 Yc`<S
 
pldisp,1                                      !绘制边坡模型变形图 q*8lnk  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 BU6Jyuwn  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 2 9#]Vr  
!边坡在强度折减系数F=2.6时结果分析 y3IWfiz>/d  
Resume,'F2.6','db'                            !读入边坡在强度折减系数F=2.6时 =
;{^"#r\  
set,1,last                                    !读入后一个子步 6y Wc1  
pldisp,1                                      !绘制边坡模型变形图 r{[OJc!  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 (oaYF+T  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 mqFq_UX/T  
!边坡在强度折减系数F=2.8时结果分析 6sB$<#  
Resume,'F2.8','db'                            !读入边坡在强度折减系数F=2.8时 ;&f1vi4  
set,1,last                                    !读入后一个子步 \%5MAQS  
pldisp,1                                      !绘制边坡模型变形图 %fhNxR  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 r]LCvsVa  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 !/
hsJ9  
!边坡在强度折减系数F=3.0时结果分析 X:=c5*0e  
Resume,'F3.0','db'                            !读入边坡在强度折减系数F=3.0时 2P9J' L  
set,1,last                                    !读入后一个子步 2o5;Uz1{  
pldisp,1                                      !绘制边坡模型变形图 kX2d7yQZz  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 }1QF+Cf  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图

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

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

好东西,顶一下!

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

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