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

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

同求 呵呵 谁知道啊

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

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

我已经自己解决了，哎

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

(1)创建物理环境 *r$.1nke  
/COM, Structural                             6'!{0 5=m  
/TITLE, Slope stability Analysis                 R9G)X]  
/FILNAM,support,1                             K/m)f#  
(2)建立模型 u@u.N2H
.%  
!进入前处理器 UVUO}B@[S  
/PREP7 })vr*[  
!定义单元类型 i9U_r._qj;  
ET,1,PLANE82                               !定义边坡围岩单元 {b
qKb=nyZ  
KEYOPT,1,3,2                                !求解类型为平面应变 %
ab)Gs  
!定义材料属性 D#}Yx]Q1  
!边坡围岩1材料属性   B/kn&^z$|~  
MP,EX,1,3E10   [iy;}5XK  
MP,PRXY,1,0.25     ,GP4I3D  
MP,DENS,1,2500                             !定义塑性模型 _k.bGYldk  
TB,DP,1                                     !定义数据表1 \&)W#8V  
TBDATA,1,9E5,42.8                           !给数据表输入数据 lILtxVBO2o  
!边坡围岩2材料属性                         !采用弹性模型 h!
CX`pBM  
MP,EX,2,3.2E10                              ！定义弹性模量 ~U6YN_W  
MP,PRXY,2,0.24                              !定义泊松比 $: m87cR~  
MP,DENS,2,2700                             !定义密度 :";D.{||  
!取强度折减系数F=1.2时边坡围岩材料属性 o 0ivja  
MP,EX,3,3E10   E wsq0D  
MP,PRXY,3,0.25     rGGepd  
MP,DENS,3,2500 h^$
}1[  
TB,DP,3                                     !定义数据表3 %kFELtx  
TBDATA,1,7.5E5,37.7                         !给数据表输入数据 ey\{C`(__y  
!取强度折减系数F=1.4时边坡围岩材料属性 3N0X?* (x|  
MP,EX,4,3E10   G5y  
MP,PRXY,4,0.25                              !定义泊松比 f<altz_\q  
MP,DENS,4,2500 rtmt3  
TB,DP,4                                    !定义数据表4 k&iScMgCTH  
TBDATA,1,6.4E5,33.5                        !给数据表输入数据 ^|i\d\  
!取强度折减系数F=1.6时边坡围岩材料属性 sxwW9_C  
MP,EX,5,3E10   &dDI*v+  
MP,PRXY,5,0.25     TM}'XZ&  
MP,DENS,5,2500                             ！定义密度 P`IG9  
TB,DP,5                                    !定义数据表5 3b\8907  
TBDATA,1,5.6E5,30                          !给数据表输入数据 fB,1s}3Hn  
!取强度折减系数F=1.8时边坡围岩材料属性 :_,]?n  
MP,EX,6,3E10                               ！定义弹性模量 q}v04Yy,o  
MP,PRXY,6,0.25     =C>`}%XT}  
MP,DENS,6,2500 |abst&yp  
TB,DP,6                                    !定义数据表6 +xBK^5/x  
TBDATA,1,5.0E5,27.2 1(7.V-(G  
!取强度折减系数F=2.0时边坡围岩材料属性 s_Oh >y?Aq  
MP,EX,7,3E10   vcJb\LW  
MP,PRXY,7,0.25                             !定义泊松比 >VUQTg  
MP,DENS,7,2500 &
W<>^C2v  
TB,DP,7                                    !定义数据表7 'L|GClc6)
 
TBDATA,1,4.5E5,24.8 V0 +k3H  
!取强度折减系数F=2.2时边坡围岩材料属性 Q>a7Ps@~  
MP,EX,8,3E10   J$9xC{L4  
MP,PRXY,8,0.25     [^}>AC*im  
MP,DENS,8,2500 zo7XmUI3P  
TB,DP,8                                     !定义数据表8 &Yf#O*  
TBDATA,1,4.09E5,22.8 \M"^Oe{Dy?  
!取强度折减系数F=2.4时边坡围岩材料属性 skDk/-*R  
MP,EX,9,3E10   :`u&TXsu  
MP,PRXY,9,0.25     M:UB>-`bW  
MP,DENS,9,2500 cm 9oG  
TB,DP,9                                      !定义数据表9 C6V&R1"s  
TBDATA,1,3.6E5,21.1                          !给数据表输入数据 V*7Z,nA  
!取强度折减系数F=2.6时边坡围岩材料属性 iVi3 :7*  
MP,EX,10,3E10                                ！定义弹性模量   KD"&_PX  
MP,PRXY,10,0.25                              !定义泊松比 Xm=^\K3  
MP,DENS,10,2500                           ！定义密度 f,HzrHax  
TB,DP,10                                   !定义数据表10 ~{,vg4L  
TBDATA,1,3.46E5,19.6                        !给数据表10输入数据 j YIV^o 0  
!取强度折减系数F=2.8时边坡围岩材料属性 *(&ClUQQ  
MP,EX,11,3E10                              ！定义弹性模量 Lr}b,  
MP,PRXY,11,0.25                             !定义泊松比 syW9Hlm  
MP,DENS,11,2500 q?-3^z%u  
TB,DP,11                                   !定义数据表11 :{i$2\DH6  
TBDATA,1,3.2E5,18.2                         !给数据表11输入数据 eMl]td rI  
!取强度折减系数F=3.0时边坡围岩材料属性 +fP/|A8P  
MP,EX,12,3E10                              ！定义弹性模量 v;bP8)mI  
MP,PRXY,12,0.25     P`3s\8[Q  
MP,DENS,12,2500                            ！定义密度 <r+!hJ[s'  
TB,DP,12                                   !定义数据表12 P-OPv%jyi  
TBDATA,1,3.0E5,17.1                        !给数据表12输入数据 &QOWW}  
!建立几何模型 IgiF,{KE,  
!创建边坡线模型线 Op/79]$  
K,1,,,,                                      ！ 创建关键点1 Ce5w0&VlS  
K,2,-800,,, <V:<x  
K,3,-800,-800,,   _95}ifSVm  
K,4,-800,-1200,, <K#'3&*$s  
K,5,1200,-1200,,   "H="Ip!s  
K,6,1200,-800,,   Jm(&G  
K,7,1200,,, hLytKPgt  
K,8,1200,378,, 0gOrW=  
K,9,430,378,, *)`kx  
l,1,2                                        !连接关键点1、2生成直线L1 s\Pt,I@Y_  
l,2,3 lKV\1(`  
l,3,4 kBiBXRt  
l,4,5 X=p3KzzX  
l,5,6 42J';\)oP  
l,6,7 Y7kb1UG  
l,7,8 ( ztim  
l,8,9 k$-~_^4m  
l,9,1 Rg?{?qK\K  
l,1,7 aq0iNbv@  
l,3,6 U%^eIXV|  
!创建边坡面模型 LE<u&
9I\  
al,3,4,5,11                              !由线L3、L4、L5、L11生成一个面A1 q1TW?\pjb:  
al,1,2,11,6,10                            !由线L1、L2、L1、L6、L10生成一个面A2 ^#BGA|j  
al,7,8,9,10                              !由线L7、L8、L9、L10生成一个面A3 y  KYP  
/PNUM, area, on                         !打开面号开关 ;N$
0)2w  
!划分网格生成有限元模型 A`x -L  
! 划分边坡围岩2网格 9o`7Kc/g  
!设置网格份数 @k+%y'Y?  
lsel,s,,,3,5,2                               !选择线L3和L5 n-hvh-ZO  
lesize,all,,,5                               !把所选择线分为5段 ]/o12pI  
lsel,s,,,4,11,7                              !选择线L4和L11 Jny)uo8  
lesize,all,,,26                              !把所选择线分为26段 O@r.>  
mat,2                                    !给边坡围岩2赋予2号材料特性 zY1s7/$i  
type,1                                   !采用单元类型1 ovDJ{3L6O  
mshkey,0                                 !设定自由网格划分  z _O,Y  
mshape,0                                 !设定四边形网格划分 Y%fVt|  
amesh,1                                  ！划分面积A1 4z9#M;qT  
! 划分边坡围岩1网格 CP]S-o}yd  
!设置网格份数 KM/U?`6>:  
lsel,s,,,2,6,4                               !选择线L2和L6 (zs4#ja2,  
lsel,a,,,1                                  !附加选择线L1 0eqi1;$b]  
lesize,all,,,10                              !把所选择线分为10段 (aJP: ^  
lsel,s,,,8,10,2                              !选择线L8和L10 q[)q|R|  
lesize,all,,,16                               !把所选择线分为16段 &SjHrOG?  
lsel,s,,,7,9,2 H}
sS4[z  
lesize,all,,,12 S$jV|xKB  
mat,1                                     !给边坡围岩1赋予1号材料特性 BSfm?ku"!  
type,1 bB:r]*_ s]  
mshkey,0                                  !设定自由网格划分 *^@#X-NG  
mshape,0                                  !设定四边形网格划分 5?5-;H  
amesh,2                                   !划分面积A2 {K45~ha9!m  
amesh,3                                   ！划分面积A3 .0 s[{x  
allsel zNV!@Yr  
!保存网格模型 ?E+:]j_  
save,Slope-grid.db
$(@o$%d  
(3)施加约束和荷载 cebs.sF:  
!给边坡模型施加约束 MegE--h  
!边坡两侧施加X方向约束 lh\ICN\O  
nsel,s,loc,x,-800                             !选择X=-800线上所有节点 X-)RU?  
nsel,a,loc,x,1200                             !选择X=1200线上所有节点 fO^e+Mz  
d,all,ux                                     !对所选择节点约束X方向位移 r=~WMDCz@  
allsel pHen>BA[  
!边坡底部施加约束 RIy5ww}3|  
nsel,s,loc,y,-1200                           !选择Y=--1200线上所有节点 %@lV-(5q  
d,all,uy                                   !对所选择节点约束X、Y方向位移 h"%|\o+3  
d,all,ux                                   29Gwv  
!施加重力加速度 ~!]&>n;=G  
acel,,9.8 0%xR<<gir  
(4)求解 sK`~Csb iB  
/solu G$;] ?g  
!求解设置 \~@[QGKN  
antype,static                                  !设定为静力求解 'yPCZ`5H(  
nsubst,100                                   !设定最大子步数为100 ~A`&/U  
pred,on                                      !打开时间步长预测器  <+AIt  
nropt,full                                     !设定牛顿-拉普森选项 N5 SLF4R1  
nlgeom,on                                    !打开大位移效果 {W5ydHXy  
lnsrch,on                                     !打开线性搜索 eg"=H50  
outres,all,all                                  !输出所有项 >=-GD2WK  
cnvtol,f,,0.005,2,0.5                            !力收敛准则设定 1B)Y;hg6&  
cnvtol,u,,0.05,2,1                              !位移收敛准则设定 TL},Unq  
!边坡在强度折减系数F=1时求解 hrs#ZZ:E  
allsel YgE]d?_h  
solve                                        !进行求解 pk-yj~F}  
save,F1,db                                    !把F=1时求解结果保存 < /;Q8;
0  
!边坡在强度折减系数F=1.2时求解 <w0NPrS]  
finish A1T<  
/solu vKNt$]p
m=  
allsel 1Pk mg%+  
asel,s,area,,2,3,1                               !选择面积A2、A3 E rop9T1  
MPCHG,3,all                                 ! 把所选择单元材料号改为3 @br@[RpB  
allsel FI
]P<)*r  
solve                                         !进行求解 }P<Qz^sr_  
save,F1.2,db                                   !把F=1.2时求解结果保存                                                                                   }>MP{67Dm  
!边坡在强度折减系数F=1.4时求解 q>_vE{UB  
finish _:35d1[  
/solu B{7Kzwh;  
allsel u0m5JD0/  
asel,s,area,,2,3,1                               !选择面积A2、A3 -VS9`7k  
MPCHG,4,all                                 !把所选择单元材料号改为4 kOs(?=  
allsel '[Mlmgc5  
solve                                        !进行求解 #yW.o'S+  
save,F1.4,db                                  !把F=1.4时求解结果保存                             J\het2?\  
!边坡在强度折减系数F=1.6时求解 [YC=d1F5  
finish EY+/.=$x  
/solu I jZ]_*^!  
allsel k?3S  
asel,s,area,,2,3,1                               !选择面积A2、A3 t)-*.qZh  
MPCHG,5,all                                 !把所选择单元材料号改为5 }.0Bl&\UK  
allsel ^)&Ly_xrU  
solve eeuZUf+~]  
save,F1.6,db                                   !把F=1.6时求解结果保存                                                                                   [Q4_WKI0T  
!边坡在强度折减系数F=1.8时求解 -Y*VgoK%  
finish ^"3\iA:  
/solu PX;Vo~6  
allsel ;% 2wGT  
asel,s,area,,2,3,1                                !选择面积A2、A3 L
nP3z5d(  
MPCHG,6,all                                  !把所选择单元材料号改为6 SArfczoB  
allsel P!kw;x  
solve e8$l0gzaD  
save,F1.8,db                                   !把F=1.8时求解结果保存   6!P`XTTE  
!边坡在强度折减系数F=2.0时求解 P
DRnW  
finish N>
+P WE$  
/solu :w@F?:C  
allsel <,\ `Psa)N  
asel,s,area,,2,3,1                                !选择面积A2、A3 &^V~cJ  
MPCHG,7,all                                   !把所选择单元材料号改为7 gRBSt M&hU  
allsel y Wpi|  
solve YiD-F7hf.*  
save,F2.0,db                                   !把F=2.0时求解结果保存                            
}$o*  
!边坡在强度折减系数F=2.2时求解 !CVuw  
finish z0#-)AeS  
/solu >; aCf#q  
allsel |#{-.r6Y]  
asel,s,area,,2,3,1                                !选择面积A2、A3 3pvYi<<D'  
MPCHG,8,all                                  !把所选择单元材料号改为8 U
1OFDXHG  
allsel  s[3e=N  
solve k!5m@'f  
save,F2.2,db                                   !把F=2.2时求解结果保存                                                                                   $tI]rU  
!边坡在强度折减系数F=2.4时求解 2ck4C/ h  
finish ujU=JlJ7dl  
/solu g %f*ofb  
allsel <"W?<VjO  
asel,s,area,,2,3,1                                 !选择面积A2、A3 dXmV@ Noo  
MPCHG,9,all                                   !把所选择单元材料号改为9 Ng\/)^  
allsel ~A6"sb=  
solve _@Y"$V]=Vt  
save,F2.4,db                                   !把F=2.4时求解结果保存   ?a S%  
!边坡在强度折减系数F=2.6时求解 [`d$X^<y;  
finish :z]}ZZ  
/solu O7%2v@j|8  
allsel
!<&m]K  
asel,s,area,,2,3,1                                 !选择面积A2、A3 Mtm/}I  
MPCHG,10,all                                   !把所选择单元材料号改为10 ^$!987"  
allsel (ab{F5  
solve ;x-(kIiE  
save,F2.6,db                                   !把F=2.6时求解结果保存                             OA[&Za#w  
!边坡在强度折减系数F=2.8时求解 P}0*{%jB  
finish 0
o;O`/x  
/solu $f#agq_  
allsel g?B3!,!9  
asel,s,area,,2,3,1                                !选择面积A2、A3 X 7=fX~s  
MPCHG,11,all                                 !把所选择单元材料号改为11 3{KR {B#L  
allsel ['z!{Ez  
solve (:5G#?6,  
save,F2.8,db                                   !把F=2.8时求解结果保存                                                                                   -T7%dLHY  
!边坡在强度折减系数F=3.0时求解 Y$g}XN*)E  
finish Z%m-HE:k  
/solu 9M_(He -  
allsel 9:5:`'b  
asel,s,area,,2,3,1                             !选择面积A2、A3 +iy7e6P  
MPCHG,12,all                              !把所选择单元材料号改为12 Y2C9(Zk U  
allsel XAPYpBgm  
solve 3t)v%S|k  
save,F3.0,db                                 !把F=3.0时求解结果保存 1"}cdq.  
/post1                                       ！进入后处理 77V .["=7  
!边坡在强度折减系数F=1时结果分析 .e@> 
 
Resume,'F1','db'                              !读入边坡在强度折减系数F=1时 9Y/L?km_(  
set,1,last                                    !读入后一个子步 Cp#)wxi6[y  
pldisp,1                                      !绘制边坡模型变形图 A3HF,EG  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 $J.T$0pFa  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 Z

U^IH9  
!边坡在强度折减系数F=1.2时结果分析 '_dzcN,z  
Resume,'F1.2','db'                            !读入边坡在强度折减系数F=1.2时 ~]
BMrgn  
set,1,last                                    !读入后一个子步 3J{vt"dS  
pldisp,1                                      !绘制边坡模型变形图 Rs*vm  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 JWROYED  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 p~2UUmV  
!边坡在强度折减系数F=1.4时结果分析 m*Lo|F  
Resume,'F1.4','db'                            !读入边坡在强度折减系数F=1.4时 q_ %cbAcD  
set,1,last                                    !读入后一个子步 m#f{]+6U  
pldisp,1                                      !绘制边坡模型变形图 z%1{  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 -I":Z2.fR  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 (DDyK[t+VX  
!边坡在强度折减系数F=1.6时结果分析 P}V=*g
  
Resume,'F1.6','db'                            !读入边坡在强度折减系数F=1.6时 k;I &.H  
set,1,last                                    !读入后一个子步 +E/y ~s  
pldisp,1                                      !绘制边坡模型变形图 1DE@N1l  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 dA_YL?or  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 S_v(S^x6  
!边坡在强度折减系数F=1.8时结果分析 F[saP0 *  
Resume,'F1.8','db'                            !读入边坡在强度折减系数F=1.8时 fTqC:r|st  
set,1,last                                    !读入后一个子步 *f5l=lDOB  
pldisp,1                                      !绘制边坡模型变形图 HSN8O@dy  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 Q$ri=uB;+  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 >`'O7.R  
!边坡在强度折减系数F=2.0时结果分析 Ri#H.T<'  
Resume,'F2.0','db'                            !读入边坡在强度折减系数F=2.0时 &3Z?UhH  
set,1,last                                    !读入后一个子步 xY\0zQ  
pldisp,1                                      !绘制边坡模型变形图 N8v'70  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 -kpswP  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 "CQw/qZw  
!边坡在强度折减系数F=2.2时结果分析 dRI^@n  
Resume,'F2.2','db'                            !读入边坡在强度折减系数F=2.2时 ow6*Xr8eQ  
set,1,last                                    !读入后一个子步 w8iR|TV  
pldisp,1                                      !绘制边坡模型变形图
RKZ6}q1n  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 `TtXZ[gP}  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 6_gnEve h  
!边坡在强度折减系数F=2.4时结果分析 :B|Dr v  
Resume,'F2.4','db'                            !读入边坡在强度折减系数F=2.4时 PWB(5 f?  
set,1,last                                    !读入后一个子步 w~Ff%p@9  
pldisp,1                                      !绘制边坡模型变形图 ZDx@^P y  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 &<5+!cV=  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 AW,OHSXh6  
!边坡在强度折减系数F=2.6时结果分析 9!U@"~yB  
Resume,'F2.6','db'                            !读入边坡在强度折减系数F=2.6时 -?6MU~"GK  
set,1,last                                    !读入后一个子步 PXz
T6)  
pldisp,1                                      !绘制边坡模型变形图  U4
7}QDh  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 j>Ag\@2ME  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 T*~H m  
!边坡在强度折减系数F=2.8时结果分析  npp[@*~  
Resume,'F2.8','db'                            !读入边坡在强度折减系数F=2.8时 W`z 0"  
set,1,last                                    !读入后一个子步 M
+`Hg_#Q  
pldisp,1                                      !绘制边坡模型变形图 .>pgU{C`!  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 zf[`~g  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图 i)ASsYG!  
!边坡在强度折减系数F=3.0时结果分析 w$|l{VI  
Resume,'F3.0','db'                            !读入边坡在强度折减系数F=3.0时 dQb.BOI)h  
set,1,last                                    !读入后一个子步 1tMQqI`N  
pldisp,1                                      !绘制边坡模型变形图 R?66b{O  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 k(%QIJH  
plnsol,eppl,eqv                               !绘制边坡模型塑性应变云图

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

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

好东西,顶一下!

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

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