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

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

同求 呵呵 谁知道啊

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

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

我已经自己解决了，哎

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

(1)创建物理环境 t*H|*L#YR  
/COM, Structural                             z4$9,p `  
/TITLE, Slope stability Analysis                 #9~,d<H  
/FILNAM,support,1                             Mc!LC .8  
(2)建立模型 c27(en(  
!进入前处理器 b#R$P]dr=  
/PREP7 fNfa.0s  
!定义单元类型 dnQ6Ras  
ET,1,PLANE82                               !定义边坡围岩单元 6a_U[-a9;  
KEYOPT,1,3,2                                !求解类型为平面应变 .GL@`7"  
!定义材料属性 S\Le;,
5Z  
!边坡围岩1材料属性   XXA]ukj;r  
MP,EX,1,3E10   9Oo`4  
MP,PRXY,1,0.25     5FF28C)>/  
MP,DENS,1,2500                             !定义塑性模型 +{C)^!zBK  
TB,DP,1                                     !定义数据表1 NJVkn~<  
TBDATA,1,9E5,42.8                           !给数据表输入数据 hZ!kh3@:`  
!边坡围岩2材料属性                         !采用弹性模型 # ,eC&X45  
MP,EX,2,3.2E10                              ！定义弹性模量 rm$dv%q  
MP,PRXY,2,0.24                              !定义泊松比 <5P*uZ  
MP,DENS,2,2700                             !定义密度 r/"^{0;F{W  
!取强度折减系数F=1.2时边坡围岩材料属性 \NEk B&^n  
MP,EX,3,3E10   pl)?4[`LUc  
MP,PRXY,3,0.25     ,2+d+Zuh  
MP,DENS,3,2500 8#- Nx]VM  
TB,DP,3                                     !定义数据表3 $5&~gHc,  
TBDATA,1,7.5E5,37.7                         !给数据表输入数据 V0W4M%  
!取强度折减系数F=1.4时边坡围岩材料属性 56.JBBZZ  
MP,EX,4,3E10   9!Jt}n?!g  
MP,PRXY,4,0.25                              !定义泊松比 svb7-.!  
MP,DENS,4,2500 %<[{zd1C-  
TB,DP,4                                    !定义数据表4 =&qH%S
6  
TBDATA,1,6.4E5,33.5                        !给数据表输入数据 cAot+N+9|]  
!取强度折减系数F=1.6时边坡围岩材料属性 gA}<Y  
MP,EX,5,3E10   ^;.u}W  
MP,PRXY,5,0.25     z@pa;
_  
MP,DENS,5,2500                             ！定义密度 /!A"[Tyt  
TB,DP,5                                    !定义数据表5 ZGsd cnz  
TBDATA,1,5.6E5,30                          !给数据表输入数据 M_.Jmh<&&  
!取强度折减系数F=1.8时边坡围岩材料属性 1#Hr{&2  
MP,EX,6,3E10                               ！定义弹性模量 }#%YeCA?  
MP,PRXY,6,0.25     IYXN}M.=  
MP,DENS,6,2500 @:#J^CsM+'  
TB,DP,6                                    !定义数据表6 jYFmL_{  
TBDATA,1,5.0E5,27.2 /h.{g0Xc  
!取强度折减系数F=2.0时边坡围岩材料属性 =Y6W Qf  
MP,EX,7,3E10   hsS&|7Pt  
MP,PRXY,7,0.25                             !定义泊松比 D %)L"5C  
MP,DENS,7,2500 gsM^Pu09ud  
TB,DP,7                                    !定义数据表7 .=t:Uy  
TBDATA,1,4.5E5,24.8 Ip=QtNW3\  
!取强度折减系数F=2.2时边坡围岩材料属性 }1~9i'o%Z  
MP,EX,8,3E10   q5-i=lw  
MP,PRXY,8,0.25     6NhGTLI  
MP,DENS,8,2500 W*CRxGyZCl  
TB,DP,8                                     !定义数据表8 B)1(  
TBDATA,1,4.09E5,22.8 gO"G/  
!取强度折减系数F=2.4时边坡围岩材料属性 u a_w5o7  
MP,EX,9,3E10   $x]/|u/9  
MP,PRXY,9,0.25     Ol+D"k~<C  
MP,DENS,9,2500 L?N-uocT  
TB,DP,9                                      !定义数据表9 !K}W.yv,  
TBDATA,1,3.6E5,21.1                          !给数据表输入数据 8wOscL f:  
!取强度折减系数F=2.6时边坡围岩材料属性 ,+6u6  
MP,EX,10,3E10                                ！定义弹性模量   qGEp 6b H  
MP,PRXY,10,0.25                              !定义泊松比 Hi
U)q  
MP,DENS,10,2500                           ！定义密度 -WF((s;<#  
TB,DP,10                                   !定义数据表10 j|K;Yi  
TBDATA,1,3.46E5,19.6                        !给数据表10输入数据 ha=2isq  
!取强度折减系数F=2.8时边坡围岩材料属性 j4=iHnE;  
MP,EX,11,3E10                              ！定义弹性模量 oVn&L*H  
MP,PRXY,11,0.25                             !定义泊松比 Ng*-Bw)p]  
MP,DENS,11,2500 PuGs%{$(h  
TB,DP,11                                   !定义数据表11 \R]2YY`EP  
TBDATA,1,3.2E5,18.2                         !给数据表11输入数据 $L6R,%c  
!取强度折减系数F=3.0时边坡围岩材料属性 %X %zK1  
MP,EX,12,3E10                              ！定义弹性模量 V@
-Q&K#  
MP,PRXY,12,0.25     ~M} K]Li  
MP,DENS,12,2500                            ！定义密度 tp7$t#  
TB,DP,12                                   !定义数据表12 S2*sh2-&6  
TBDATA,1,3.0E5,17.1                        !给数据表12输入数据 &?H$-r1/?V  
!建立几何模型 GWPBP-)0  
!创建边坡线模型线 S\:+5}  
K,1,,,,                                      ！ 创建关键点1 I#O"<0 *r  
K,2,-800,,, E&y)`>Nq{  
K,3,-800,-800,,   .IdbaH _a  
K,4,-800,-1200,, 2sk7E'2(  
K,5,1200,-1200,,   ln=fq:  
K,6,1200,-800,,   TaHcvjhR  
K,7,1200,,, 2z;3NUL$n  
K,8,1200,378,, 'r%(,=L  
K,9,430,378,, M80Q6K  
l,1,2                                        !连接关键点1、2生成直线L1 Z.Z;p/4F  
l,2,3 (&/4wI^M  
l,3,4 wLqj<ot  
l,4,5 ppvlU H5;  
l,5,6 Q*ELMib  
l,6,7 ?z l<"u  
l,7,8 s=jYQ5nv  
l,8,9 h`N2M,  
l,9,1 ny1O- `!1  
l,1,7 0LL0
\ly]  
l,3,6 Er~KX3vF  
!创建边坡面模型 Um4zI>  
al,3,4,5,11                              !由线L3、L4、L5、L11生成一个面A1
C{uT1`  
al,1,2,11,6,10                            !由线L1、L2、L1、L6、L10生成一个面A2 {fHor  
al,7,8,9,10                              !由线L7、L8、L9、L10生成一个面A3 8_HBcZWs  
/PNUM, area, on                         !打开面号开关 qwj7CIc(  
!划分网格生成有限元模型 (P-^ PNz&  
! 划分边坡围岩2网格 ^$lZ  
!设置网格份数 hZ e{Ri  
lsel,s,,,3,5,2                               !选择线L3和L5 _1G/qHf^S  
lesize,all,,,5                               !把所选择线分为5段 j;
_  
lsel,s,,,4,11,7                              !选择线L4和L11 S 4 17.n  
lesize,all,,,26                              !把所选择线分为26段 ly*v|(S&  
mat,2                                    !给边坡围岩2赋予2号材料特性 J.4U;A5  
type,1                                   !采用单元类型1 mKO~`Wq%@  
mshkey,0                                 !设定自由网格划分 HOWpTu(  
mshape,0                                 !设定四边形网格划分 H2: Zda#  
amesh,1                                  ！划分面积A1 \l/}` w  
! 划分边坡围岩1网格 SIrNZ^I  
!设置网格份数 C@`#@1X  
lsel,s,,,2,6,4                               !选择线L2和L6 p$1 'e,G  
lsel,a,,,1                                  !附加选择线L1 ^t gjs$M|  
lesize,all,,,10                              !把所选择线分为10段 oydP}X  
lsel,s,,,8,10,2                              !选择线L8和L10 _p0Yhju?  
lesize,all,,,16                               !把所选择线分为16段 ^"?a)KC  
lsel,s,,,7,9,2 q%kCTw  
lesize,all,,,12 >_yL@^  
mat,1                                     !给边坡围岩1赋予1号材料特性 _^GBfM.  
type,1 :K2 X~Ty  
mshkey,0                                  !设定自由网格划分 o<!H/PN  
mshape,0                                  !设定四边形网格划分 '{ =F/q  
amesh,2                                   !划分面积A2 CS~onf<xz  
amesh,3                                   ！划分面积A3 IL.bwtpQD  
allsel -Z,r\9d  
!保存网格模型 V9"R8*@-  
save,Slope-grid.db 0$Rn|yqf%  
(3)施加约束和荷载 v\?J=|S+  
!给边坡模型施加约束 ms}o[Z@n  
!边坡两侧施加X方向约束 #?S"y:  
nsel,s,loc,x,-800                             !选择X=-800线上所有节点 ,S`n?.&& 7  
nsel,a,loc,x,1200                             !选择X=1200线上所有节点 t`Z3*?UqI  
d,all,ux                                     !对所选择节点约束X方向位移 |Sjy   
allsel p!:oT1U  
!边坡底部施加约束 pIK:$eN!/  
nsel,s,loc,y,-1200                           !选择Y=--1200线上所有节点 ?o+%ckH  
d,all,uy                                   !对所选择节点约束X、Y方向位移 "QXnE^  
d,all,ux                                   Y3[KS;_fr9  
!施加重力加速度 ]m 3cm  
acel,,9.8 `H:`JBe=+[  
(4)求解 AA,/AKikd  
/solu +:jonN9d  
!求解设置 oA-,>:}g{  
antype,static                                  !设定为静力求解 KM[0aXOtv  
nsubst,100                                   !设定最大子步数为100 @uRJl$3  
pred,on                                      !打开时间步长预测器 7=?!B#hm!  
nropt,full                                     !设定牛顿-拉普森选项 nrev!h  
nlgeom,on                                    !打开大位移效果 --l UEo~  
lnsrch,on                                     !打开线性搜索 <CO_JWD  
outres,all,all                                  !输出所有项 G9okl9;od  
cnvtol,f,,0.005,2,0.5                            !力收敛准则设定 sNpA!!\PM  
cnvtol,u,,0.05,2,1                              !位移收敛准则设定 |T"vF`Kr(>  
!边坡在强度折减系数F=1时求解 d]sqj\Q57  
allsel T:{&eWH  
solve                                        !进行求解 V^Z5i]zT  
save,F1,db                                    !把F=1时求解结果保存 =Y|( }92  
!边坡在强度折减系数F=1.2时求解 ^TZmc{i  
finish hq{{XQ  
/solu %b*N.v1+  
allsel $-UVN0=  
asel,s,area,,2,3,1                               !选择面积A2、A3 /YZMP'v  
MPCHG,3,all                                 ! 把所选择单元材料号改为3 ELnUpmv\  
allsel Yn [ F:Z  
solve                                         !进行求解 &D{!zF  
save,F1.2,db                                   !把F=1.2时求解结果保存                                                                                   00i MU  
!边坡在强度折减系数F=1.4时求解 "TaLvworb4  
finish 9b=0 4aWHm  
/solu 9kWI2cLzQt  
allsel zT)cg$8%fY  
asel,s,area,,2,3,1                               !选择面积A2、A3 !?ZR_=Y%  
MPCHG,4,all                                 !把所选择单元材料号改为4 :X}SuM?c  
allsel 8k:^( kByF  
solve                                        !进行求解 :atd_6   
save,F1.4,db                                  !把F=1.4时求解结果保存                             uZL,%pF3A  
!边坡在强度折减系数F=1.6时求解 oZ/"^5  
finish S[$9_Jf  
/solu 7)<Ib j<M  
allsel {"~[F2qR  
asel,s,area,,2,3,1                               !选择面积A2、A3 Xux[  
MPCHG,5,all                                 !把所选择单元材料号改为5 :
35h0;8+  
allsel XZ.D<T"  
solve M/.M~/~  
save,F1.6,db                                   !把F=1.6时求解结果保存                                                                                   q j21#q .  
!边坡在强度折减系数F=1.8时求解 E)m \KSwh  
finish ]8}2  
/solu {9Qc\Ij  
allsel \,/ozfJ7dT  
asel,s,area,,2,3,1                                !选择面积A2、A3 )?ra
dg  
MPCHG,6,all                                  !把所选择单元材料号改为6 Jxe5y3* (  
allsel k]:`<`/I_  
solve .?g=mh79(  
save,F1.8,db                                   !把F=1.8时求解结果保存   3,G|oR{D  
!边坡在强度折减系数F=2.0时求解 Y1r'\@L w  
finish ~b#<
HG\,,  
/solu gh[q*%#  
allsel rz|Sjtq  
asel,s,area,,2,3,1                                !选择面积A2、A3 Z;~
%!  
MPCHG,7,all                                   !把所选择单元材料号改为7 AaB1H7r-  
allsel lt4IoE`tk?  
solve XN t` 4$L  
save,F2.0,db                                   !把F=2.0时求解结果保存                             r0k:RJP  
!边坡在强度折减系数F=2.2时求解 R?lTB3"  
finish S"Vr+x?  
/solu Z)xcxSo  
allsel d!<>Fh^6,  
asel,s,area,,2,3,1                                !选择面积A2、A3 & ;5f/  
MPCHG,8,all                                  !把所选择单元材料号改为8 YY tVp_)  
allsel GNmP_N  
solve e^\(bp+83  
save,F2.2,db                                   !把F=2.2时求解结果保存                                                                                   q'H6oD`  
!边坡在强度折减系数F=2.4时求解 qOpwl*?x+  
finish w-)JCdS6Tb  
/solu W^[QEmyn  
allsel VWdTnu  
asel,s,area,,2,3,1                                 !选择面积A2、A3 r`FTiPD.C  
MPCHG,9,all                                   !把所选择单元材料号改为9 -ys/I,}<  
allsel V+>RF  
solve 3_;=y\F  
save,F2.4,db                                   !把F=2.4时求解结果保存   0[ "CP:u  
!边坡在强度折减系数F=2.6时求解 n@kJ1ee'  
finish .XK3o .ZhW  
/solu U<XfO'XJ  
allsel aW|=|K  
asel,s,area,,2,3,1                                 !选择面积A2、A3 }3Ke  
MPCHG,10,all                                   !把所选择单元材料号改为10 J Z@sk2  
allsel wb#[&2
i  
solve M^l%*QF[,q  
save,F2.6,db                                   !把F=2.6时求解结果保存                             A`vRUl,c=  
!边坡在强度折减系数F=2.8时求解 jZ5ac=D&I  
finish A9Ea}v9:  
/solu ||cI~qg  
allsel Qt'3v"S>)  
asel,s,area,,2,3,1                                !选择面积A2、A3 jsV1~1:83  
MPCHG,11,all                                 !把所选择单元材料号改为11 7Y.yl F:  
allsel F@t\D?  
solve
=Ldf#8J  
save,F2.8,db                                   !把F=2.8时求解结果保存                                                                                   'dQGb-<_<  
!边坡在强度折减系数F=3.0时求解 0.!vp?  
finish mDwuJf8}  
/solu (NnE\2  
allsel Y~( 8<`^  
asel,s,area,,2,3,1                             !选择面积A2、A3 }|=Fnyj  
MPCHG,12,all                              !把所选择单元材料号改为12 P%w!4v~"  
allsel ^N ;TCn  
solve .gRb'  
save,F3.0,db                                 !把F=3.0时求解结果保存 hQzT =0  
/post1                                       ！进入后处理 g:c?%J  
!边坡在强度折减系数F=1时结果分析 J~
@W":v  
Resume,'F1','db'                              !读入边坡在强度折减系数F=1时 4yMi9Ri4H  
set,1,last                                    !读入后一个子步 XI ><;#  
pldisp,1                                      !绘制边坡模型变形图 0 >(hiTy<  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 QRdtr  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 Ttb?x<)+8  
!边坡在强度折减系数F=1.2时结果分析 n]l3 )u  
Resume,'F1.2','db'                            !读入边坡在强度折减系数F=1.2时 a(-t"OL\  
set,1,last                                    !读入后一个子步 SZxnYVY  
pldisp,1                                      !绘制边坡模型变形图 [8DPZU@  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 ij_5=4aZ
-  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 .UK`~17!  
!边坡在强度折减系数F=1.4时结果分析 9j;!4AJ1t  
Resume,'F1.4','db'                            !读入边坡在强度折减系数F=1.4时 7|ACJv6%9  
set,1,last                                    !读入后一个子步 h^H)p`[Gme  
pldisp,1                                      !绘制边坡模型变形图 Yv
{$
XI7  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 |D;I>
O^"R  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 "[_gRe*2  
!边坡在强度折减系数F=1.6时结果分析 +WxD=|p;  
Resume,'F1.6','db'                            !读入边坡在强度折减系数F=1.6时 70R_O&f-k  
set,1,last                                    !读入后一个子步 ^|h.B$_F,  
pldisp,1                                      !绘制边坡模型变形图 o|s|Wmx>u  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 >n&+<06  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 Q`=d5Uvw  
!边坡在强度折减系数F=1.8时结果分析 6%c]{eTd9  
Resume,'F1.8','db'                            !读入边坡在强度折减系数F=1.8时 1NN#-U  
set,1,last                                    !读入后一个子步 7;#o?6!7  
pldisp,1                                      !绘制边坡模型变形图 y?*Y=,"  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 }mk z_P(Z  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 %NhZTmWm  
!边坡在强度折减系数F=2.0时结果分析  </Dv?  
Resume,'F2.0','db'                            !读入边坡在强度折减系数F=2.0时 p@YbIn  
set,1,last                                    !读入后一个子步
!|,djo!N  
pldisp,1                                      !绘制边坡模型变形图 eN TKX  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 )RT:u)N  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 *@ S+J$  
!边坡在强度折减系数F=2.2时结果分析 A22'qgKm@  
Resume,'F2.2','db'                            !读入边坡在强度折减系数F=2.2时 @Rqn&tA8  
set,1,last                                    !读入后一个子步 sdQ"[`~2R  
pldisp,1                                      !绘制边坡模型变形图 ]PH'G>x  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 DL '{ rK  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 Uo:=-NNI  
!边坡在强度折减系数F=2.4时结果分析 .Jc<Gg  
Resume,'F2.4','db'                            !读入边坡在强度折减系数F=2.4时 s<LYSrd  
set,1,last                                    !读入后一个子步 Q"GZh.m  
pldisp,1                                      !绘制边坡模型变形图 O 0P4uq  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 m8*)@e  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 =kW7|c5Z  
!边坡在强度折减系数F=2.6时结果分析 7t% |s!~  
Resume,'F2.6','db'                            !读入边坡在强度折减系数F=2.6时 eyByAT~W,  
set,1,last                                    !读入后一个子步 bqZ5GKUo  
pldisp,1                                      !绘制边坡模型变形图 Z~g7^,-t  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 ~Y7:08  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 @w;&:J9m
 
!边坡在强度折减系数F=2.8时结果分析 u5
P2*  
Resume,'F2.8','db'                            !读入边坡在强度折减系数F=2.8时 HhqNpU
 
set,1,last                                    !读入后一个子步 &
UX:KW`=  
pldisp,1                                      !绘制边坡模型变形图 @_"Z]Y ,D0  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 T^nOv2@,  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 T5TAkEVl  
!边坡在强度折减系数F=3.0时结果分析 5k<0>6;XH  
Resume,'F3.0','db'                            !读入边坡在强度折减系数F=3.0时 wvEdZGO8!  
set,1,last                                    !读入后一个子步 CGZ3-OW@E  
pldisp,1                                      !绘制边坡模型变形图 WP?AQD  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 U;Q?Rh-W  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图

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

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

好东西,顶一下!

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

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