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

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

同求 呵呵 谁知道啊

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

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

我已经自己解决了，哎

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

(1)创建物理环境 X#>:9  
/COM, Structural                             az1#:Go  
/TITLE, Slope stability Analysis                 ]++,7Z\AU  
/FILNAM,support,1                             8ZF!}kb0F  
(2)建立模型 JT!
Cb$!  
!进入前处理器 7J;\&q'  
/PREP7 R > [2*o"  
!定义单元类型 V!Pe%.>  
ET,1,PLANE82                               !定义边坡围岩单元 tR\cS)  
KEYOPT,1,3,2                                !求解类型为平面应变 j#jwK(:]  
!定义材料属性 pkf$%{"e  
!边坡围岩1材料属性   '2Lx>nByk  
MP,EX,1,3E10   (?xR<]~g*  
MP,PRXY,1,0.25     +bGO"*  
MP,DENS,1,2500                             !定义塑性模型 Z<iK(?@O  
TB,DP,1                                     !定义数据表1 Y?6}r;<  
TBDATA,1,9E5,42.8                           !给数据表输入数据 .lAPlJOO  
!边坡围岩2材料属性                         !采用弹性模型 _;,"!'R`f  
MP,EX,2,3.2E10                              ！定义弹性模量 VGf&'nL@,  
MP,PRXY,2,0.24                              !定义泊松比 VXnWY8\  
MP,DENS,2,2700                             !定义密度 j+PLtE  
!取强度折减系数F=1.2时边坡围岩材料属性 ZPG~@lU  
MP,EX,3,3E10   }X6w"  
MP,PRXY,3,0.25     zZW5M^z8  
MP,DENS,3,2500 Nwo*tb:  
TB,DP,3                                     !定义数据表3 P(.XB`  
TBDATA,1,7.5E5,37.7                         !给数据表输入数据 yAFt|<  
!取强度折减系数F=1.4时边坡围岩材料属性 +G6 Ge;  
MP,EX,4,3E10   eH V#Mey[  
MP,PRXY,4,0.25                              !定义泊松比 1QPz|3f@\  
MP,DENS,4,2500 mM r$~^P:  
TB,DP,4                                    !定义数据表4 mT/^F{c  
TBDATA,1,6.4E5,33.5                        !给数据表输入数据 Ob/i_  
!取强度折减系数F=1.6时边坡围岩材料属性 x];i? 4  
MP,EX,5,3E10   |\Q2L;4C  
MP,PRXY,5,0.25     vm7ag 7@O  
MP,DENS,5,2500                             ！定义密度 PXG@]$~3  
TB,DP,5                                    !定义数据表5 o!lKP>  
TBDATA,1,5.6E5,30                          !给数据表输入数据 jwg*\HO,s  
!取强度折减系数F=1.8时边坡围岩材料属性 h$#PboLd  
MP,EX,6,3E10                               ！定义弹性模量 ;H3~r^>c  
MP,PRXY,6,0.25     ;jJ4H+8  
MP,DENS,6,2500 #uV J  
TB,DP,6                                    !定义数据表6 !"ir}Y%  
TBDATA,1,5.0E5,27.2 H.;2o(vD  
!取强度折减系数F=2.0时边坡围岩材料属性 kCALJRf~d  
MP,EX,7,3E10   g~2=he\C  
MP,PRXY,7,0.25                             !定义泊松比 CkRilS<  
MP,DENS,7,2500 icQQLSU5  
TB,DP,7                                    !定义数据表7 =w3cF)&  
TBDATA,1,4.5E5,24.8 (f2r4Io|}  
!取强度折减系数F=2.2时边坡围岩材料属性
@ULd~  
MP,EX,8,3E10   Xf{9rZ+  
MP,PRXY,8,0.25     _f~$iY  
MP,DENS,8,2500 bk1.H@8  
TB,DP,8                                     !定义数据表8 [Op^l%BC  
TBDATA,1,4.09E5,22.8 )5`~WzA  
!取强度折减系数F=2.4时边坡围岩材料属性 o@hj.)u  
MP,EX,9,3E10   H& $M/`  
MP,PRXY,9,0.25     \<hHZS  
MP,DENS,9,2500 XSCcumde!  
TB,DP,9                                      !定义数据表9 oWx^_wQ-=  
TBDATA,1,3.6E5,21.1                          !给数据表输入数据 +^jm_+  
!取强度折减系数F=2.6时边坡围岩材料属性 ?-~<Vc*  
MP,EX,10,3E10                                ！定义弹性模量   wA"d?x  
MP,PRXY,10,0.25                              !定义泊松比 Qp7F3,/#  
MP,DENS,10,2500                           ！定义密度 )W6l/  
TB,DP,10                                   !定义数据表10 |S~$IFN4  
TBDATA,1,3.46E5,19.6                        !给数据表10输入数据 (5A8#7a  
!取强度折减系数F=2.8时边坡围岩材料属性 8;"9A  
MP,EX,11,3E10                              ！定义弹性模量 v0TbQ  
MP,PRXY,11,0.25                             !定义泊松比 s)?GscPG!  
MP,DENS,11,2500 N{}8Zh4op  
TB,DP,11                                   !定义数据表11 r/@Wn  
TBDATA,1,3.2E5,18.2                         !给数据表11输入数据 P#MUS_x  
!取强度折减系数F=3.0时边坡围岩材料属性 NHQoP&OG  
MP,EX,12,3E10                              ！定义弹性模量 m?gGFxo  
MP,PRXY,12,0.25     >5E1y!  
MP,DENS,12,2500                            ！定义密度 E,dUO;  
TB,DP,12                                   !定义数据表12 lWH#/5`h  
TBDATA,1,3.0E5,17.1                        !给数据表12输入数据 =7#"}%4Q  
!建立几何模型 VT-%o7%N  
!创建边坡线模型线 LqO=wK~  
K,1,,,,                                      ！ 创建关键点1 r,3\32[?  
K,2,-800,,, XwfR/4  
K,3,-800,-800,,   .'.bokl/  
K,4,-800,-1200,, D>^g2!b:  
K,5,1200,-1200,,   ?~IdPSY  
K,6,1200,-800,,   h4=mGJpm
 
K,7,1200,,, U<
r!G;^`  
K,8,1200,378,, obN8+ j  
K,9,430,378,, _;%l~q/  
l,1,2                                        !连接关键点1、2生成直线L1 y\4/M6  
l,2,3
cs_  
l,3,4 [c -|`d^  
l,4,5 ?bu=QV@  
l,5,6 .S1MxZhbP  
l,6,7 u\MxQIo'u  
l,7,8 '@ p464  
l,8,9 db=S*LUbl  
l,9,1 i`Qa7  
l,1,7 GY %$7
  
l,3,6 a#G3dY>  
!创建边坡面模型 6nJQPa  
al,3,4,5,11                              !由线L3、L4、L5、L11生成一个面A1 6tjV^sjs  
al,1,2,11,6,10                            !由线L1、L2、L1、L6、L10生成一个面A2 a\~118 !  
al,7,8,9,10                              !由线L7、L8、L9、L10生成一个面A3 miTff[hsMa  
/PNUM, area, on                         !打开面号开关 2#00<t\  
!划分网格生成有限元模型 $RB p!7  
! 划分边坡围岩2网格 %!AzFL J|Z  
!设置网格份数 \rATmjsKzS  
lsel,s,,,3,5,2                               !选择线L3和L5 DSWmQQ  
lesize,all,,,5                               !把所选择线分为5段 ._t1eb`m{  
lsel,s,,,4,11,7                              !选择线L4和L11 @v:Eh  
lesize,all,,,26                              !把所选择线分为26段 \YFM5l;IU  
mat,2                                    !给边坡围岩2赋予2号材料特性 lFduX D  
type,1                                   !采用单元类型1 xX9snSGz  
mshkey,0                                 !设定自由网格划分 yTn<5T[H  
mshape,0                                 !设定四边形网格划分 ZR-64G=L,  
amesh,1                                  ！划分面积A1 ycwkF$7  
! 划分边坡围岩1网格 s&'FaqE  
!设置网格份数 |lZJt  
lsel,s,,,2,6,4                               !选择线L2和L6 00s)=A_  
lsel,a,,,1                                  !附加选择线L1 WbBd<^Q  
lesize,all,,,10                              !把所选择线分为10段 #_`qbIOAj  
lsel,s,,,8,10,2                              !选择线L8和L10 2`%a[t@M.  
lesize,all,,,16                               !把所选择线分为16段 : 2$*'{mM  
lsel,s,,,7,9,2 a1Q%Gn@R  
lesize,all,,,12 >qOj^WO~  
mat,1                                     !给边坡围岩1赋予1号材料特性 l!KPgRw  
type,1 #,#:{&H  
mshkey,0                                  !设定自由网格划分 pq\N2d  
mshape,0                                  !设定四边形网格划分 ywkRH  
amesh,2                                   !划分面积A2 Sv0?_3C  
amesh,3                                   ！划分面积A3 $.:x3TsA  
allsel #U7pT!Fx  
!保存网格模型 * hS6F  
save,Slope-grid.db {IvA 5^  
(3)施加约束和荷载 ~t/i0pKq.  
!给边坡模型施加约束 aG]^8`~>'  
!边坡两侧施加X方向约束 KZ%us6  
nsel,s,loc,x,-800                             !选择X=-800线上所有节点 g]c6_DMfb1  
nsel,a,loc,x,1200                             !选择X=1200线上所有节点 1 ?Zw  
d,all,ux                                     !对所选择节点约束X方向位移 >dJ~  
allsel (fr=N5   
!边坡底部施加约束 {o5^nd  
nsel,s,loc,y,-1200                           !选择Y=--1200线上所有节点 ZBFn  
d,all,uy                                   !对所选择节点约束X、Y方向位移 9$U4x|n  
d,all,ux                                   W}2!~ep!  
!施加重力加速度 f[.'V1  
acel,,9.8 Ctn?O~u  
(4)求解 /3tErc'  
/solu c#`IF6qj  
!求解设置 S-ZN}N{,6  
antype,static                                  !设定为静力求解 R".$x{{  
nsubst,100                                   !设定最大子步数为100 ;UUgqX#  
pred,on                                      !打开时间步长预测器 ='=4tj=z  
nropt,full                                     !设定牛顿-拉普森选项 ~tV7yY|zr  
nlgeom,on                                    !打开大位移效果 I{WP:]"Yf  
lnsrch,on                                     !打开线性搜索 @V:Y%#%  
outres,all,all                                  !输出所有项 sKsMF:|OT  
cnvtol,f,,0.005,2,0.5                            !力收敛准则设定 osdoL  
cnvtol,u,,0.05,2,1                              !位移收敛准则设定 m
k^,{D  
!边坡在强度折减系数F=1时求解 INp:;  
allsel @D=%J!!*  
solve                                        !进行求解 -*~ @?  
save,F1,db                                    !把F=1时求解结果保存 ;hV|W{=w  
!边坡在强度折减系数F=1.2时求解
W,[QK~  
finish )8eb(!}7  
/solu (=fLWK{8  
allsel &WL::gy_S  
asel,s,area,,2,3,1                               !选择面积A2、A3 nCp_RJu  
MPCHG,3,all                                 ! 把所选择单元材料号改为3 ,EVPnH[F~  
allsel L6i|5 P  
solve                                         !进行求解 aaw[ia_EL  
save,F1.2,db                                   !把F=1.2时求解结果保存                                                                                   9wGsHf8]  
!边坡在强度折减系数F=1.4时求解 [hpkE lE  
finish #gT"G1
8/!  
/solu )ZBNw{nh  
allsel aPK:k$.  
asel,s,area,,2,3,1                               !选择面积A2、A3 ]];pWlo!  
MPCHG,4,all                                 !把所选择单元材料号改为4 &7>]# *  
allsel N-&ZaK  
solve                                        !进行求解
+:=FcsY  
save,F1.4,db                                  !把F=1.4时求解结果保存                             S@S4<R1{\  
!边坡在强度折减系数F=1.6时求解 #Ha"rr46p  
finish jAue+tB  
/solu -sZ'<(3  
allsel XqUQ{^;aI  
asel,s,area,,2,3,1                               !选择面积A2、A3 c2,g%(  
MPCHG,5,all                                 !把所选择单元材料号改为5 jF j'6LT9/  
allsel
O"EL3$9V  
solve j[FB*L1!D  
save,F1.6,db                                   !把F=1.6时求解结果保存                                                                                   WQze|b%  
!边坡在强度折减系数F=1.8时求解 Uf]$I`T#  
finish GYiL}itD=3  
/solu DTM(SN8R+n  
allsel "aI)LlyCY  
asel,s,area,,2,3,1                                !选择面积A2、A3 `GY3H3B  
MPCHG,6,all                                  !把所选择单元材料号改为6 t']/2m.&p  
allsel |2n*Ds'  
solve o#1Ta7Ro  
save,F1.8,db                                   !把F=1.8时求解结果保存   "v`q%(TA  
!边坡在强度折减系数F=2.0时求解 wbcip8<t  
finish sJ^Ff  
/solu U
BM#~~sM  
allsel rw3tU0j  
asel,s,area,,2,3,1                                !选择面积A2、A3 do.>Y}d  
MPCHG,7,all                                   !把所选择单元材料号改为7 4&%H;Q  
allsel =q)+_@24>d  
solve ]v),[]Xs  
save,F2.0,db                                   !把F=2.0时求解结果保存                             rP2h9Cb  
!边坡在强度折减系数F=2.2时求解 5S!j$_(  
finish L;?F^RK{U  
/solu dTCLE t.  
allsel `Npo|.?=  
asel,s,area,,2,3,1                                !选择面积A2、A3 $joGda  
MPCHG,8,all                                  !把所选择单元材料号改为8 1@ &J"*  
allsel D
O-M0L  
solve a73VDQr I  
save,F2.2,db                                   !把F=2.2时求解结果保存                                                                                   ZsXw]Wa  
!边坡在强度折减系数F=2.4时求解 WN%,  
finish *{k
{  
/solu }fkdv6mz  
allsel 4g2`[<S  
asel,s,area,,2,3,1                                 !选择面积A2、A3 %saP>]o  
MPCHG,9,all                                   !把所选择单元材料号改为9 }`H{;A h  
allsel "eOl(TSu/  
solve <pa-C2Ky  
save,F2.4,db                                   !把F=2.4时求解结果保存   IAMa  
!边坡在强度折减系数F=2.6时求解 50Gu~No6  
finish sb3z8:r  
/solu y
(22m+B  
allsel ,X:3w3nr^  
asel,s,area,,2,3,1                                 !选择面积A2、A3 wme#8/eUk  
MPCHG,10,all                                   !把所选择单元材料号改为10 LEtGrA/%@b  
allsel N}NKQ]=  
solve /ar0K9`c  
save,F2.6,db                                   !把F=2.6时求解结果保存                             ')yYpWO  
!边坡在强度折减系数F=2.8时求解 {dxl8~/I  
finish 7G;1n0m-T  
/solu u$zRm(!RB  
allsel $M0l (htR  
asel,s,area,,2,3,1                                !选择面积A2、A3 .57Fh)Y  
MPCHG,11,all                                 !把所选择单元材料号改为11 1pG|jT+Bi  
allsel oMLs22Do?  
solve 4:Ju|g]O  
save,F2.8,db                                   !把F=2.8时求解结果保存                                                                                   }Rh%bf7,  
!边坡在强度折减系数F=3.0时求解 [\<#iRcP  
finish |.yS~XFJS  
/solu A&B|n!;b  
allsel $8&Y(`  
asel,s,area,,2,3,1                             !选择面积A2、A3 NvTK7? v  
MPCHG,12,all                              !把所选择单元材料号改为12 >b9J!'G,(  
allsel *,:2O&P  
solve 2vG X\W%3  
save,F3.0,db                                 !把F=3.0时求解结果保存 &nVekE:!  
/post1                                       ！进入后处理 Qnt}:M+  
!边坡在强度折减系数F=1时结果分析 ntPj9#lf  
Resume,'F1','db'                              !读入边坡在强度折减系数F=1时 C{nk,j L  
set,1,last                                    !读入后一个子步 J?dz>3Rhx9  
pldisp,1                                      !绘制边坡模型变形图 V6_":L"!  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 :eD-'#@$u  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 Hf1b&8&:K  
!边坡在强度折减系数F=1.2时结果分析 YqWNp  
Resume,'F1.2','db'                            !读入边坡在强度折减系数F=1.2时 !t~tIJ>6  
set,1,last                                    !读入后一个子步 <*^|Aj|#  
pldisp,1                                      !绘制边坡模型变形图 +_*NY~  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 &J|I&p   
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 F;_o`h  
!边坡在强度折减系数F=1.4时结果分析 Wtk|}>Pf  
Resume,'F1.4','db'                            !读入边坡在强度折减系数F=1.4时 *H
HL a  
set,1,last                                    !读入后一个子步 ;vPFRiFK  
pldisp,1                                      !绘制边坡模型变形图 sUmpf4/  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 'WmjQsf  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 9ei'o
Z  
!边坡在强度折减系数F=1.6时结果分析 T|h!06   
Resume,'F1.6','db'                            !读入边坡在强度折减系数F=1.6时 '-n Iy$>  
set,1,last                                    !读入后一个子步 B:UPSX)A  
pldisp,1                                      !绘制边坡模型变形图 ecH7")  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 :8}Qt^p  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 rwW"B  
!边坡在强度折减系数F=1.8时结果分析 '~f*O0_  
Resume,'F1.8','db'                            !读入边坡在强度折减系数F=1.8时 &FOq c  
set,1,last                                    !读入后一个子步 ;=^J_2ls  
pldisp,1                                      !绘制边坡模型变形图 /g9^g(  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 3Rsrb  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 (6mw@gzr  
!边坡在强度折减系数F=2.0时结果分析 9e
P*N(m<  
Resume,'F2.0','db'                            !读入边坡在强度折减系数F=2.0时 L>WxAeyu1K  
set,1,last                                    !读入后一个子步 |E$q S)y  
pldisp,1                                      !绘制边坡模型变形图 >c Tt2v  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 lFcHE c  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 a_[+id  
!边坡在强度折减系数F=2.2时结果分析 =MM+(mD  
Resume,'F2.2','db'                            !读入边坡在强度折减系数F=2.2时 0LzS #J+  
set,1,last                                    !读入后一个子步 tBZ?UAe;  
pldisp,1                                      !绘制边坡模型变形图 ,|?#+O{  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 |?^N@  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 .=G3wox3  
!边坡在强度折减系数F=2.4时结果分析 :@E^oNKa0  
Resume,'F2.4','db'                            !读入边坡在强度折减系数F=2.4时 Pj g#  
set,1,last                                    !读入后一个子步 nTs\zikP  
pldisp,1                                      !绘制边坡模型变形图 <_*5BO  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 ^2%)Nq;O  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 `}uM91;  
!边坡在强度折减系数F=2.6时结果分析 fU%Ys9:wU  
Resume,'F2.6','db'                            !读入边坡在强度折减系数F=2.6时 =s'H o  
set,1,last                                    !读入后一个子步 d Efk~V\  
pldisp,1                                      !绘制边坡模型变形图 #n.v#FyNx  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 b">"NvlB  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 0;=]MEk?  
!边坡在强度折减系数F=2.8时结果分析 )>Z@')Uk:  
Resume,'F2.8','db'                            !读入边坡在强度折减系数F=2.8时 o.|36#Fa  
set,1,last                                    !读入后一个子步 Er$&}9G+-  
pldisp,1                                      !绘制边坡模型变形图 yLXIjR  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 uX{g4#eG  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 ?aBAmyxm  
!边坡在强度折减系数F=3.0时结果分析 u6Ux nqNc  
Resume,'F3.0','db'                            !读入边坡在强度折减系数F=3.0时 i 
Pl/I  
set,1,last                                    !读入后一个子步 n4XkhY|  
pldisp,1                                      !绘制边坡模型变形图 }7{(o-  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 glM42s  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图

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

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

好东西,顶一下!

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

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