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

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

同求 呵呵 谁知道啊

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

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

我已经自己解决了，哎

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

(1)创建物理环境 `7+j0kV)  
/COM, Structural                             ogeL[7  
/TITLE, Slope stability Analysis                 m-V02's  
/FILNAM,support,1                             @!Hr|k|  
(2)建立模型 y69J%/c ra  
!进入前处理器 Y\+KoR';  
/PREP7 :XV} c(+d  
!定义单元类型 MS<SAD>w  
ET,1,PLANE82                               !定义边坡围岩单元 LS.r%:$mb  
KEYOPT,1,3,2                                !求解类型为平面应变 E-ZRG!)[v  
!定义材料属性 ,Vd7V}t  
!边坡围岩1材料属性   =y%rG :!  
MP,EX,1,3E10   9l+{OA  
MP,PRXY,1,0.25     zzZg$9PT[  
MP,DENS,1,2500                             !定义塑性模型 BRbV7&  
TB,DP,1                                     !定义数据表1 XTq
m]  
TBDATA,1,9E5,42.8                           !给数据表输入数据 #T~&]|{,  
!边坡围岩2材料属性                         !采用弹性模型 +&ZX$  
MP,EX,2,3.2E10                              ！定义弹性模量 b X)|MiWI  
MP,PRXY,2,0.24                              !定义泊松比 8v yG*UK  
MP,DENS,2,2700                             !定义密度 $_j1kx$  
!取强度折减系数F=1.2时边坡围岩材料属性 Az,- C
q  
MP,EX,3,3E10   Y:x/!-  
MP,PRXY,3,0.25     GlR~%q-jiQ  
MP,DENS,3,2500 , fb( WY  
TB,DP,3                                     !定义数据表3 h ?#@~  
TBDATA,1,7.5E5,37.7                         !给数据表输入数据 :#="%  
!取强度折减系数F=1.4时边坡围岩材料属性 Jm(ixekp  
MP,EX,4,3E10   I__4I{nI  
MP,PRXY,4,0.25                              !定义泊松比 We`'>'W0  
MP,DENS,4,2500 m-1?\bs  
TB,DP,4                                    !定义数据表4 L2-^!'  
TBDATA,1,6.4E5,33.5                        !给数据表输入数据 Jp;k+"<q  
!取强度折减系数F=1.6时边坡围岩材料属性 j%& IL0  
MP,EX,5,3E10   qeaA&(|5  
MP,PRXY,5,0.25     d*;$AYI#R  
MP,DENS,5,2500                             ！定义密度 eZF'Ck
y  
TB,DP,5                                    !定义数据表5 oEzDMImJ5  
TBDATA,1,5.6E5,30                          !给数据表输入数据 D\]&8w6&  
!取强度折减系数F=1.8时边坡围岩材料属性 }++5_Z_  
MP,EX,6,3E10                               ！定义弹性模量 r(uP!n1+  
MP,PRXY,6,0.25     f\ "`7  
MP,DENS,6,2500 4/wwn6I}G  
TB,DP,6                                    !定义数据表6 M%=P)cC  
TBDATA,1,5.0E5,27.2 \~"#ld(x7  
!取强度折减系数F=2.0时边坡围岩材料属性 U6WG?$x  
MP,EX,7,3E10   dc^Vc{26Z  
MP,PRXY,7,0.25                             !定义泊松比 RR>G]#k  
MP,DENS,7,2500 EI`vVI  
TB,DP,7                                    !定义数据表7 %"7WXOv&z  
TBDATA,1,4.5E5,24.8 gM_:l  
!取强度折减系数F=2.2时边坡围岩材料属性 mf*Nr0L;J  
MP,EX,8,3E10   9[>Lp9l'  
MP,PRXY,8,0.25     6S*zzJ.0K  
MP,DENS,8,2500 t.3\/  
TB,DP,8                                     !定义数据表8 ha'qIT3&  
TBDATA,1,4.09E5,22.8 B1w0cS%%:  
!取强度折减系数F=2.4时边坡围岩材料属性 =g0*MZ;
"  
MP,EX,9,3E10   ;?im(9h"v!  
MP,PRXY,9,0.25     rycJyiw<-  
MP,DENS,9,2500 5upShtC  
TB,DP,9                                      !定义数据表9 Xu<FDjr  
TBDATA,1,3.6E5,21.1                          !给数据表输入数据 xw%)rm<t  
!取强度折减系数F=2.6时边坡围岩材料属性 _O52ai><b  
MP,EX,10,3E10                                ！定义弹性模量   r%!FmS<  
MP,PRXY,10,0.25                              !定义泊松比 z<QIuq  
MP,DENS,10,2500                           ！定义密度 3}M\c)
 
TB,DP,10                                   !定义数据表10 /Nqrvy=  
TBDATA,1,3.46E5,19.6                        !给数据表10输入数据 iwz` x  
!取强度折减系数F=2.8时边坡围岩材料属性 `Yogq)G}  
MP,EX,11,3E10                              ！定义弹性模量 `!kL1oUYE  
MP,PRXY,11,0.25                             !定义泊松比 eJf>"IF-  
MP,DENS,11,2500 wF;B@  
TB,DP,11                                   !定义数据表11 T#e4":A&x  
TBDATA,1,3.2E5,18.2                         !给数据表11输入数据 Dh8(HiXf:  
!取强度折减系数F=3.0时边坡围岩材料属性 ~j}J<4&OvC  
MP,EX,12,3E10                              ！定义弹性模量 JFJIls  
MP,PRXY,12,0.25     s*{mT6s+T  
MP,DENS,12,2500                            ！定义密度 x(yX0 ,P/7  
TB,DP,12                                   !定义数据表12 E8;TLk4\  
TBDATA,1,3.0E5,17.1                        !给数据表12输入数据 -y<x!61  
!建立几何模型 I3ho(Kdi
  
!创建边坡线模型线 x#o?>5Qg?  
K,1,,,,                                      ！ 创建关键点1 U$@83?O{iM  
K,2,-800,,, o~}1oN  
K,3,-800,-800,,   ~GYpat  
K,4,-800,-1200,, (l;C%O7*  
K,5,1200,-1200,,   _3KZME  
K,6,1200,-800,,   !SN WB  
K,7,1200,,, ox JGJ  
K,8,1200,378,, t_Z _!Qy  
K,9,430,378,, MyM+C}  
l,1,2                                        !连接关键点1、2生成直线L1 $*aE$O6l  
l,2,3 Ahd\TH  
l,3,4 7n [12:  
l,4,5 jXH0BPa,  
l,5,6 ~\-r  
l,6,7 %,Q;<axzi  
l,7,8 dd19z%  
l,8,9 afV P-m4L  
l,9,1 W
kcH5[  
l,1,7 ?bn;{c;E  
l,3,6 p( HyRCH  
!创建边坡面模型 xf[zEEt  
al,3,4,5,11                              !由线L3、L4、L5、L11生成一个面A1 kYmo7  
al,1,2,11,6,10                            !由线L1、L2、L1、L6、L10生成一个面A2 +0dQORo  
al,7,8,9,10                              !由线L7、L8、L9、L10生成一个面A3 1wU=WE(kKZ  
/PNUM, area, on                         !打开面号开关 @>gD1Q7v b  
!划分网格生成有限元模型 o9~h%&  
! 划分边坡围岩2网格 &s2#1  
!设置网格份数 HUD7{6
}4  
lsel,s,,,3,5,2                               !选择线L3和L5 3[,wM
y"  
lesize,all,,,5                               !把所选择线分为5段 S]5VEn;pV  
lsel,s,,,4,11,7                              !选择线L4和L11 c-ttds  
lesize,all,,,26                              !把所选择线分为26段 L;4[ k;5  
mat,2                                    !给边坡围岩2赋予2号材料特性 J]$er0`LY  
type,1                                   !采用单元类型1 li1v 4  
mshkey,0                                 !设定自由网格划分 $8EV,9^U  
mshape,0                                 !设定四边形网格划分 d"LoK,p#  
amesh,1                                  ！划分面积A1 n=;';(wR[  
! 划分边坡围岩1网格
LAizx^F  
!设置网格份数 ;K>{_kf  
lsel,s,,,2,6,4                               !选择线L2和L6 0P3j+? N%  
lsel,a,,,1                                  !附加选择线L1 J&n ^y  
lesize,all,,,10                              !把所选择线分为10段 5#yJK>a7  
lsel,s,,,8,10,2                              !选择线L8和L10 uT'-B7N  
lesize,all,,,16                               !把所选择线分为16段 ?,D>+::  
lsel,s,,,7,9,2 s+(l7xH$  
lesize,all,,,12 0>I]=M]@  
mat,1                                     !给边坡围岩1赋予1号材料特性 y*y`t6D  
type,1 uDpf2(>s
 
mshkey,0                                  !设定自由网格划分 AB/,S  
mshape,0                                  !设定四边形网格划分 }g%KvYB_  
amesh,2                                   !划分面积A2 `4$4bXrP'  
amesh,3                                   ！划分面积A3 ( [K2:n\  
allsel %"1` NT  
!保存网格模型 NvEm,E\|  
save,Slope-grid.db %kFTnX
HK  
(3)施加约束和荷载 .vCY%0oE  
!给边坡模型施加约束 L"Qh_+  
!边坡两侧施加X方向约束 j|3g(_v4W  
nsel,s,loc,x,-800                             !选择X=-800线上所有节点 7Qm;g-)f  
nsel,a,loc,x,1200                             !选择X=1200线上所有节点 Y\+(rC27  
d,all,ux                                     !对所选择节点约束X方向位移 ?ZSXoy-kr  
allsel Eqz4{\   
!边坡底部施加约束 J-@o@!o  
nsel,s,loc,y,-1200                           !选择Y=--1200线上所有节点 z Y$X|=f  
d,all,uy                                   !对所选择节点约束X、Y方向位移 U:s}/to  
d,all,ux                                   va`l*N5  
!施加重力加速度 3bT6W,J4T  
acel,,9.8 2P;%P]~H  
(4)求解 i6meY$l  
/solu BFzcoBu-  
!求解设置 :t2 9`x  
antype,static                                  !设定为静力求解 f_Wkg)g  
nsubst,100                                   !设定最大子步数为100 zbF:R[)  
pred,on                                      !打开时间步长预测器 N/`TrWVF  
nropt,full                                     !设定牛顿-拉普森选项 FYx `o\  
nlgeom,on                                    !打开大位移效果 hLDch5J5~  
lnsrch,on                                     !打开线性搜索 ~ LJ>WA  
outres,all,all                                  !输出所有项 h|(ZXCH  
cnvtol,f,,0.005,2,0.5                            !力收敛准则设定 &cv@Kihq(  
cnvtol,u,,0.05,2,1                              !位移收敛准则设定 PK:o}IWn~x  
!边坡在强度折减系数F=1时求解 U}A|]vi@  
allsel F20wf1^  
solve                                        !进行求解 b7W=HR  
save,F1,db                                    !把F=1时求解结果保存 y
(aAp.S>  
!边坡在强度折减系数F=1.2时求解 )[@YHE5g  
finish :Y}Y&mA4  
/solu )vE
HLp.  
allsel cJ7{4YK_#/  
asel,s,area,,2,3,1                               !选择面积A2、A3 4~MJ4:  
MPCHG,3,all                                 ! 把所选择单元材料号改为3 7/p J6>  
allsel m>Yo9/XpZ  
solve                                         !进行求解 z*NC?\  
save,F1.2,db                                   !把F=1.2时求解结果保存                                                                                   #Lhj0M;a  
!边坡在强度折减系数F=1.4时求解 '%n<MTL  
finish rm5T=fNJ  
/solu YJ$Vn>6Z  
allsel ex@,F,u>o  
asel,s,area,,2,3,1                               !选择面积A2、A3 /el["l
 
MPCHG,4,all                                 !把所选择单元材料号改为4 ] dm1Qm  
allsel &6,Yjs:T m  
solve                                        !进行求解 z^a6%N  
save,F1.4,db                                  !把F=1.4时求解结果保存                             ]RJb;  
!边坡在强度折减系数F=1.6时求解 g|v1qfK  
finish #Vn>ue+?  
/solu %ojR?=ON  
allsel r{/
G\  
asel,s,area,,2,3,1                               !选择面积A2、A3 9QX
~aX  
MPCHG,5,all                                 !把所选择单元材料号改为5 zDvP7hl  
allsel [0tfY0  
solve .%J?T5D  
save,F1.6,db                                   !把F=1.6时求解结果保存                                                                                   U'Ja\Ek/f  
!边坡在强度折减系数F=1.8时求解 (A]m=  
finish c5wkzY h  
/solu <.AC=4@V  
allsel iXG>j.w{79  
asel,s,area,,2,3,1                                !选择面积A2、A3 #pDGaqeX  
MPCHG,6,all                                  !把所选择单元材料号改为6 8XH|T^5  
allsel XRz%KVysp  
solve :cmI"Bo  
save,F1.8,db                                   !把F=1.8时求解结果保存   R+kZLOE  
!边坡在强度折减系数F=2.0时求解 }`<>$2b  
finish 53,,%Ue  
/solu },Y; (n'  
allsel ~f]r>jQM  
asel,s,area,,2,3,1                                !选择面积A2、A3 Qa7S'(  
MPCHG,7,all                                   !把所选择单元材料号改为7 v]`A_)[  
allsel ].<sAmL^  
solve &JXHDpd$a^  
save,F2.0,db                                   !把F=2.0时求解结果保存                             $cJN9|$6  
!边坡在强度折减系数F=2.2时求解 |qy"%W@  
finish yf#%)-7(  
/solu ]\D6;E8P-~  
allsel YPff)0Nh  
asel,s,area,,2,3,1                                !选择面积A2、A3 F~Z0  
MPCHG,8,all                                  !把所选择单元材料号改为8 Br ^rK}|l  
allsel T956L'.+G  
solve -t~B@%  
save,F2.2,db                                   !把F=2.2时求解结果保存                                                                                   4_m /_Z0x  
!边坡在强度折减系数F=2.4时求解 Zs5I?R1e8  
finish @R OY}CZ{/  
/solu SX?$H~A  
allsel &TkbnDuYd~  
asel,s,area,,2,3,1                                 !选择面积A2、A3 Lx:N!RDw  
MPCHG,9,all                                   !把所选择单元材料号改为9 ap'kxOf"1  
allsel hRy}G'0  
solve = C'e1=]  
save,F2.4,db                                   !把F=2.4时求解结果保存   }D
UDA%U  
!边坡在强度折减系数F=2.6时求解 1=q?#PQ  
finish 5KH'|z  
/solu 0h@%q;g  
allsel :{xN33@6\X  
asel,s,area,,2,3,1                                 !选择面积A2、A3 Qqhb]<z  
MPCHG,10,all                                   !把所选择单元材料号改为10 mN'9|`>V>  
allsel u<y\iZ[   
solve k>0cTBY&  
save,F2.6,db                                   !把F=2.6时求解结果保存                             kN9sug^  
!边坡在强度折减系数F=2.8时求解 P7x?!71?L  
finish iUG/   
/solu ?n9$,-^v  
allsel + $a:X  
asel,s,area,,2,3,1                                !选择面积A2、A3 W];6u  
MPCHG,11,all                                 !把所选择单元材料号改为11 4Ig{#}<  
allsel ~"rwP=<}  
solve K!: ,l  
save,F2.8,db                                   !把F=2.8时求解结果保存                                                                                   !Ub?eJp  
!边坡在强度折减系数F=3.0时求解 8rG&CxI  
finish [rQ(ae  
/solu wCTR-pL^  
allsel {kRC!}  
asel,s,area,,2,3,1                             !选择面积A2、A3 <K[Zl/7I  
MPCHG,12,all                              !把所选择单元材料号改为12 /qq&'}TZP  
allsel ,/b/O4`;y  
solve 3F'dT[;  
save,F3.0,db                                 !把F=3.0时求解结果保存 @[b:([  
/post1                                       ！进入后处理 8
}#Lo9:,d  
!边坡在强度折减系数F=1时结果分析 D_ZBx+/_?  
Resume,'F1','db'                              !读入边坡在强度折减系数F=1时 8`~3MsE"  
set,1,last                                    !读入后一个子步 5NF&LM;i(  
pldisp,1                                      !绘制边坡模型变形图 bSmaE7  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 H6JMN1#t$  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 Okg8Ve2  
!边坡在强度折减系数F=1.2时结果分析 l`%} {3r9  
Resume,'F1.2','db'                            !读入边坡在强度折减系数F=1.2时 >(6\ C  
set,1,last                                    !读入后一个子步 UuqnL{  
pldisp,1                                      !绘制边坡模型变形图 e/Oj T  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 z57papo  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 dDn4nwH  
!边坡在强度折减系数F=1.4时结果分析 BG8)bhk;/  
Resume,'F1.4','db'                            !读入边坡在强度折减系数F=1.4时 ;)]zv\fC  
set,1,last                                    !读入后一个子步
;yER V  
pldisp,1                                      !绘制边坡模型变形图 WAa1H60VkS  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 :=7'1H  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 rLVc<595  
!边坡在强度折减系数F=1.6时结果分析 gwz _b  
Resume,'F1.6','db'                            !读入边坡在强度折减系数F=1.6时 u R%R]X  
set,1,last                                    !读入后一个子步 r~2@#gTbl  
pldisp,1                                      !绘制边坡模型变形图 KC-aLq/  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 \fR:+rbQ&|  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 4$2
T zJE  
!边坡在强度折减系数F=1.8时结果分析 eC`f8=V  
Resume,'F1.8','db'                            !读入边坡在强度折减系数F=1.8时 F}.TT=((8  
set,1,last                                    !读入后一个子步 *Z2Q]?:{ i  
pldisp,1                                      !绘制边坡模型变形图 m.a
1  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 +sluu!~  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 d lLk4a+  
!边坡在强度折减系数F=2.0时结果分析 _r|ytQ)  
Resume,'F2.0','db'                            !读入边坡在强度折减系数F=2.0时 >`|uc  
set,1,last                                    !读入后一个子步 =l'_*B8  
pldisp,1                                      !绘制边坡模型变形图 HPdwx V  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 &8i{'k,l  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 \ g(#)f  
!边坡在强度折减系数F=2.2时结果分析 wP1VQUL  
Resume,'F2.2','db'                            !读入边坡在强度折减系数F=2.2时 nJ})6/gK  
set,1,last                                    !读入后一个子步 BF[?* b  
pldisp,1                                      !绘制边坡模型变形图 (a
!,)  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 wXGFq3`  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 }UdqX1jz  
!边坡在强度折减系数F=2.4时结果分析 VEEeQy
  
Resume,'F2.4','db'                            !读入边坡在强度折减系数F=2.4时 *dTf(J  
set,1,last                                    !读入后一个子步 3}1+"? s  
pldisp,1                                      !绘制边坡模型变形图 RX\O'Zwlj  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 p% %Y^=z  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 76r s)J[*w  
!边坡在强度折减系数F=2.6时结果分析 :Qp/3(g e  
Resume,'F2.6','db'                            !读入边坡在强度折减系数F=2.6时 k$_]b0D{4  
set,1,last                                    !读入后一个子步 FN`kSTm*0!  
pldisp,1                                      !绘制边坡模型变形图 h#ot)m|I  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 &.4_4"l(  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 m8^2k2  
!边坡在强度折减系数F=2.8时结果分析 %9Fg1LH42r  
Resume,'F2.8','db'                            !读入边坡在强度折减系数F=2.8时 QaBXzf   
set,1,last                                    !读入后一个子步 L/LNX{|  
pldisp,1                                      !绘制边坡模型变形图 '="){  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 ,T
oEKId  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 aKtTx~$@  
!边坡在强度折减系数F=3.0时结果分析 *)D*iU&  
Resume,'F3.0','db'                            !读入边坡在强度折减系数F=3.0时 deY
<+!  
set,1,last                                    !读入后一个子步 W3rvKqdw5  
pldisp,1                                      !绘制边坡模型变形图 "Oko|3  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 S$On$]~\"  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图

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

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

好东西,顶一下!

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

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