各位朋友： <Pqv;WI |R
大家好！我通过ansys用有限元强度折减法分析了一个边坡的稳定性，结果不同的工况下得出的塑性区是一样的，我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激！附命令流和塑性区如下：(1)创建物理环境 *E>.)B i
/COM, Structural ! 指定结构分析 ofc.zwH
/TITLE, Slope stability Analysis ! 定义工作标题 i;\s.wrzH
/FILNAM,support,1 ! 定义工作文件名 # jYpVc{]
(2)建立模型 h SqY$P
!进入前处理器 {B$2"q/~
/PREP7 #~:P}<h
!定义单元类型 <x%M3BTx
ET,1,PLANE82 !定义边坡围岩单元 )\/
=M*
KEYOPT,1,3,2 !求解类型为平面应变 +\`vq"e
!定义材料属性 Ar9nBJ`
!边坡围岩1材料属性 KHiFJ_3
MP,EX,1,3E10 = Je>`{J
MP,PRXY,1,0.25 jX=lAs~6
MP,DENS,1,2500 !定义塑性模型 zbw7U'jk
TB,DP,1 !定义数据表1 L6:h.1 U$
TBDATA,1,9E5,42.8 !给数据表输入数据 )9"oL!2h
!边坡围岩2材料属性 !采用弹性模型 `ue[q!Qq
MP,EX,2,3.2E10 ！定义弹性模量 khd5 Cf[
MP,PRXY,2,0.24 !定义泊松比 6y+b5- {'
MP,DENS,2,2700 !定义密度
GrJ#.
!取强度折减系数F=1.2时边坡围岩材料属性 BWPP5X9
MP,EX,3,3E10 4 |J[Jdj
MP,PRXY,3,0.25 AG" l1wz
MP,DENS,3,2500
{E9v`u\
TB,DP,3 !定义数据表3 A//?6O Jx?
TBDATA,1,7.5E5,37.7 !给数据表输入数据 =y ]Jl,_.
!取强度折减系数F=1.4时边坡围岩材料属性 O>r-]0DI[
MP,EX,4,3E10 6o3 T;h
MP,PRXY,4,0.25 !定义泊松比 hAR?
t5c
MP,DENS,4,2500 Os),;W0w 4
TB,DP,4 !定义数据表4 V*n==Nb5L
TBDATA,1,6.4E5,33.5 !给数据表输入数据 >G)qns9
!取强度折减系数F=1.6时边坡围岩材料属性 `{<f rB@
MP,EX,5,3E10 *3{J#Q6fk3
MP,PRXY,5,0.25 M.%shrJ/
MP,DENS,5,2500 ！定义密度 I %_MV
TB,DP,5 !定义数据表5 J07O:cjyu
TBDATA,1,5.6E5,30 !给数据表输入数据 54p tP
!取强度折减系数F=1.8时边坡围岩材料属性 <)m%*9{
MP,EX,6,3E10 ！定义弹性模量 =iz,S:[
MP,PRXY,6,0.25 =WZ%H_oxi
MP,DENS,6,2500 =Z0t :{
TB,DP,6 !定义数据表6 Z_z#QX>=D
TBDATA,1,5.0E5,27.2 _j#SpL' P
!取强度折减系数F=2.0时边坡围岩材料属性 WnxEu3U
MP,EX,7,3E10 ,\ldz(D?+
MP,PRXY,7,0.25 !定义泊松比 9w ~cvlv[
MP,DENS,7,2500 q|#MB7e/
TB,DP,7 !定义数据表7 t\y-T$\\
TBDATA,1,4.5E5,24.8 TYS\95<
!取强度折减系数F=2.2时边坡围岩材料属性 Rq)BssdF
MP,EX,8,3E10 IhonnLLW
MP,PRXY,8,0.25 WO </Q6+
MP,DENS,8,2500 GG-[`!>.pw
TB,DP,8 !定义数据表8 Ub)I66
TBDATA,1,4.09E5,22.8 s${_K* g6
!取强度折减系数F=2.4时边坡围岩材料属性 <&#+ E%E4
MP,EX,9,3E10 lglYJ,
MP,PRXY,9,0.25 m@qqVRn#)
MP,DENS,9,2500 e1 a*'T$z
TB,DP,9 !定义数据表9 <! x+e E`
TBDATA,1,3.6E5,21.1 !给数据表输入数据 ~t/JCxa
!取强度折减系数F=2.6时边坡围岩材料属性 s OC&Q&eg
MP,EX,10,3E10 ！定义弹性模量 %]4-{%v
MP,PRXY,10,0.25 !定义泊松比 r2eQ{u{nX
MP,DENS,10,2500 ！定义密度 Dx+ K +(
TB,DP,10 !定义数据表10 1 p\Ak
TBDATA,1,3.46E5,19.6 !给数据表10输入数据 z]#hWfM4B:
!取强度折减系数F=2.8时边坡围岩材料属性 pG yRX_;
MP,EX,11,3E10 ！定义弹性模量 }C6@c1myq-
MP,PRXY,11,0.25 !定义泊松比 *3Nn +T
MP,DENS,11,2500 H~9=&p[ Q
TB,DP,11 !定义数据表11 vZjZb(jlN
TBDATA,1,3.2E5,18.2 !给数据表11输入数据 :EHQ .^
!取强度折减系数F=3.0时边坡围岩材料属性 &TT":FPR
MP,EX,12,3E10 ！定义弹性模量 ?ApRJm:T
MP,PRXY,12,0.25 Qlz Q]:dWC
MP,DENS,12,2500 ！定义密度 ,seFkG@1
TB,DP,12 !定义数据表12 W"sr$K2m|
TBDATA,1,3.0E5,17.1 !给数据表12输入数据 I-:` cON=G
!建立几何模型 it}-^3 A M
!创建边坡线模型线 J/Q|uRpmqr
K,1,,,, ！ 创建关键点1 m6n hC
K,2,-800,,, qi=3L
K,3,-800,-800,, &hN,xpC
K,4,-800,-1200,, 'U )8rR
K,5,1200,-1200,, n(&*kfk
K,6,1200,-800,, 9 81!2*
K,7,1200,,, +WF.wP?y
K,8,1200,400,, pMN<p[MB
K,9,400,400,, *8\(FVyG^
l,1,2 !连接关键点1、2生成直线L1 r z6jx
l,2,3 'j7 9GC0
l,3,4 o/JPYBhdl
l,4,5 g _ M-F
l,5,6 ]h@{6N'oNS
l,6,7 &5q{ viI
l,7,8 h~h a
l,8,9 lO 2k <
l,9,1 {
/<4'B
l,1,7 f5?hnt`m
l,3,6 XZO<dhZX:
!创建边坡面模型 G[6=u| (M
al,3,4,5,11 !由线L3、L4、L5、L11生成一个面A1 60PYCqWc
al,1,2,11,6,10 !由线L1、L2、L1、L6、L10生成一个面A2 yiT{+;g^
al,7,8,9,10 !由线L7、L8、L9、L10生成一个面A3 `pYE[y+
/PNUM, area, on !打开面号开关 3 7>MJ
!划分网格生成有限元模型 .BL:h&h|y
! 划分边坡围岩2网格 Zr"dOj$Jf
!设置网格份数 G%OpO.Wf
lsel,s,,,3,5,2 !选择线L3和L5 jOl 1_
lesize,all,,,5 !把所选择线分为5段 A?Sm-#n{
lsel,s,,,4,11,7 !选择线L4和L11 M[,^KJ!
lesize,all,,,26 !把所选择线分为26段 ZjD2u 8e
mat,2 !给边坡围岩2赋予2号材料特性 oNSz& )LP
type,1 !采用单元类型1 "jGe^+9uT
mshkey,0 !设定自由网格划分 F*<Ws;j
mshape,0 !设定四边形网格划分 M(S{1|,V
amesh,1 ！划分面积A1 Oxn'bh6R0
! 划分边坡围岩1网格 Gg +Y fY_
!设置网格份数 -A@U0=o
lsel,s,,,2,6,4 !选择线L2和L6 2{};6{yz
lsel,a,,,1 !附加选择线L1 Wf1-" Q
lesize,all,,,10 !把所选择线分为10段 UCF[oO>v
lsel,s,,,8,10,2 !选择线L8和L10 jr1Se9u D
lesize,all,,,16 !把所选择线分为16段 WguV{#=H
lsel,s,,,7,9,2 w:R]!e_6\9
lesize,all,,,12 YHh u^}|jQ
mat,1 !给边坡围岩1赋予1号材料特性 EG1x
type,1 YPQCOG
mshkey,0 !设定自由网格划分 ^ K>pT}u
mshape,0 !设定四边形网格划分 ddMM74
amesh,2 !划分面积A2 = +Tsknq
amesh,3 ！划分面积A3 2V< # Y
allsel 3mA/Nu_
!保存网格模型 },3R%?8 9%
save,Slope-grid.db :+A; TV
(3)施加约束和荷载 Bcm=G""
!给边坡模型施加约束 fW[.r== Kf
!边坡两侧施加X方向约束 +}]wLM}\UF
nsel,s,loc,x,-800 !选择X=-800线上所有节点 d.1Q~&`
nsel,a,loc,x,1200 !选择X=1200线上所有节点 F\P!NSFZV
d,all,ux !对所选择节点约束X方向位移 KqY>4tb
allsel s.e y!ew
!边坡底部施加约束 <M>#qd@c
nsel,s,loc,y,-1200 !选择Y=--1200线上所有节点 B2"+Hwbk
d,all,uy !对所选择节点约束X、Y方向位移 % NwoU%q
d,all,ux mZ &]
!施加重力加速度 Q,TaJ ]
acel,,9.8 -CZ-l;5
(4)求解 >,w P! ;dh
/solu qZc)Sa. S
!求解设置 Cj=J;^vf
antype,static !设定为静力求解 .=aMjrME
nsubst,100 !设定最大子步数为100 Xa6qvg7/
pred,on !打开时间步长预测器 h@/>?Va
nropt,full !设定牛顿-拉普森选项 "p2u+ 8?
nlgeom,on !打开大位移效果 l?N`V2SuR
lnsrch,on !打开线性搜索 p!V) 55J*
outres,all,all !输出所有项 6P6 Jx;
cnvtol,f,,0.005,2,0.5 !力收敛准则设定 $q=h c u
cnvtol,u,,0.05,2,1 !位移收敛准则设定 Ut=0~x.=<
!边坡在强度折减系数F=1时求解 ~M(pCSJ[
allsel {.%0@{Y
solve !进行求解 \ !6t
save,F1,db !把F=1时求解结果保存 q .)^B@}_
!边坡在强度折减系数F=1.2时求解 :4b- sg#
finish [/n' @cjNZ
/solu wf^cyCR0
allsel r6/<&1[
asel,s,area,,2,3,1 !选择面积A2、A3 JC BnFrP
MPCHG,3,all ! 把所选择单元材料号改为3 _2x uzmz0
allsel 0c{Gr 0[>
solve !进行求解 13]y) (
save,F1.2,db !把F=1.2时求解结果保存 [kQ"6wh8
!边坡在强度折减系数F=1.4时求解 HJ(=?TU
finish LE Jlo%M
/solu sV;qpD XX
allsel E ,Dlaq
asel,s,area,,2,3,1 !选择面积A2、A3 M |aQ)ivh3
MPCHG,4,all !把所选择单元材料号改为4 LIc*tsl
allsel F\l!A'Q+t
solve !进行求解 ROcY'-
save,F1.4,db !把F=1.4时求解结果保存 IrUpExJ
!边坡在强度折减系数F=1.6时求解 jR }*bIzv
finish $::51# ^Wg
/solu ],WwqD=
allsel `
C/fF_YA
asel,s,area,,2,3,1 !选择面积A2、A3 z bdOCfA;
MPCHG,5,all !把所选择单元材料号改为5 _p ? I{1O
allsel 6 YB-}>?
solve ',<{X (#(
save,F1.6,db !把F=1.6时求解结果保存 !G@V<'F
!边坡在强度折减系数F=1.8时求解 an #jZ[
finish nT(!HDH
/solu UX+?0 K
allsel iOE9FW|e
asel,s,area,,2,3,1 !选择面积A2、A3 GWVE IZ
MPCHG,6,all !把所选择单元材料号改为6 BD<rQ mfA^
allsel 6B4hSqjh
solve i] V
F'tG
save,F1.8,db !把F=1.8时求解结果保存 4d,qXSKty
!边坡在强度折减系数F=2.0时求解 HE4`9$kVLr
finish *s
1D\/H
/solu i@sCMCu6
allsel |v,}%UN 2
asel,s,area,,2,3,1 !选择面积A2、A3 7 B4w.P,B
MPCHG,7,all !把所选择单元材料号改为7 _
+u sn.
allsel 1KNkl,E
solve n C\(+K1%
save,F2.0,db !把F=2.0时求解结果保存 8JFnB(3xU
!边坡在强度折减系数F=2.2时求解 * <Nk%`
finish 2*b# + b
/solu KI^ q 5D ?
allsel lU >)n
asel,s,area,,2,3,1 !选择面积A2、A3 + ;{rU&
MPCHG,8,all !把所选择单元材料号改为8 K]lb8q}Z~
allsel #h}IUR
solve U<#$w{d:
save,F2.2,db !把F=2.2时求解结果保存 i|:: v l
!边坡在强度折减系数F=2.4时求解 _S[Rvb1e
finish F=P+;%.
/solu Y,p2eAss
allsel 50O7=
asel,s,area,,2,3,1 !选择面积A2、A3 je$R\7B<
MPCHG,9,all !把所选择单元材料号改为9 il 8A&`%
allsel P W0q71
solve 96^1Ivd
save,F2.4,db !把F=2.4时求解结果保存 KR+ a Y.
!边坡在强度折减系数F=2.6时求解 VhJyWH%(
finish (dw3'W
/solu yAZ.L/j yr
allsel 0Te)s 3X
asel,s,area,,2,3,1 !选择面积A2、A3 !IQfe o T
MPCHG,10,all !把所选择单元材料号改为10 |>m# m*{S
allsel }-k_?2"A
solve ^H y)< P
save,F2.6,db !把F=2.6时求解结果保存 {g23[$X]N
!边坡在强度折减系数F=2.8时求解 #Q!c42}M
finish !rN#PF>
/solu FN8=YUYK%
allsel D:S6Mu
asel,s,area,,2,3,1 !选择面积A2、A3 N23 +1 h
MPCHG,11,all !把所选择单元材料号改为11 X8NO;w@z#
allsel . T
N`p*
solve 8z3I~yL_`+
save,F2.8,db !把F=2.8时求解结果保存 )U3 H1 5
!边坡在强度折减系数F=3.0时求解 E : LQ!
finish -hJ>wGI
/solu 9 Le/'o vq
allsel >V ]*mS %K
asel,s,area,,2,3,1 !选择面积A2、A3 %^)Ja EUC
MPCHG,12,all !把所选择单元材料号改为12 @pEO@bbg >
allsel qJAv=D
solve Sj=69>m]5
save,F3.0,db !把F=3.0时求解结果保存 W6!4Q yn
(5）后处理 1' @lg*^9
/post1 ！进入后处理 : p*ojl|
!边坡在强度折减系数F=1时结果分析 E?;T:7.%
Resume,'F1','db' !读入边坡在强度折减系数F=1时 c=U$$|qHV
set,1,last !读入后一个子步 bl+@}+A
pldisp,1 !绘制边坡模型变形图 i32_ZB Z?y
plnsol,u,x !绘制边坡模型水平方向位移云图 iM Y0xf8l
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 nP3;<*T P0
!边坡在强度折减系数F=1.2时结果分析 Z 1MJ!{@6
Resume,'F1.2','db' !读入边坡在强度折减系数F=1.2时 3
MI ) E
set,1,last !读入后一个子步 6=`m
pldisp,1 !绘制边坡模型变形图 xs83S.fHg
plnsol,u,x !绘制边坡模型水平方向位移云图 W@uH!n>k
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 9^[5!SMzCj
!边坡在强度折减系数F=1.4时结果分析 }G:uzud10
Resume,'F1.4','db' !读入边坡在强度折减系数F=1.4时 4q 2=:"z4
set,1,last !读入后一个子步
GwIfGixqH
pldisp,1 !绘制边坡模型变形图 W;)F NP|MT
plnsol,u,x !绘制边坡模型水平方向位移云图 qyG636i
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 |
9\7xT
!边坡在强度折减系数F=1.6时结果分析 q4#f
*]
Resume,'F1.6','db' !读入边坡在强度折减系数F=1.6时 Q i'WV9ke
set,1,last !读入后一个子步 M
b /X@ 51
pldisp,1 !绘制边坡模型变形图 VTgbJ {?
plnsol,u,x !绘制边坡模型水平方向位移云图 .]sIoB-5 4
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 ?H86Wbz
!边坡在强度折减系数F=1.8时结果分析 "3K Smb
Resume,'F1.8','db' !读入边坡在强度折减系数F=1.8时 ?}lp o; $
set,1,last !读入后一个子步 *Yk8Mj^_h
pldisp,1 !绘制边坡模型变形图 >cr _^(UW&
plnsol,u,x !绘制边坡模型水平方向位移云图 >[P7Zlwv4
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 ^J}$y7
!边坡在强度折减系数F=2.0时结果分析 ff5
Lwf{{
Resume,'F2.0','db' !读入边坡在强度折减系数F=2.0时 _XO)`D~
set,1,last !读入后一个子步 q MT.7n:
pldisp,1 !绘制边坡模型变形图 -aeo7C
plnsol,u,x !绘制边坡模型水平方向位移云图 h9OL%n 7m'
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 =QKgsgLh
!边坡在强度折减系数F=2.2时结果分析 ;
K
6Fe)
Resume,'F2.2','db' !读入边坡在强度折减系数F=2.2时 $r QFM[
set,1,last !读入后一个子步 A>1 p]#
pldisp,1 !绘制边坡模型变形图 J7xT6Q=
plnsol,u,x !绘制边坡模型水平方向位移云图 :`"T Eif
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 }+RF~~H/
!边坡在强度折减系数F=2.4时结果分析 'rq#q)1MT
Resume,'F2.4','db' !读入边坡在强度折减系数F=2.4时 kI[O {<kQ
set,1,last !读入后一个子步 _2Xu1q.6~5
pldisp,1 !绘制边坡模型变形图 sz2 SWk^&
plnsol,u,x !绘制边坡模型水平方向位移云图 8iH;GFNJ7'
plnsol,eppl,eqv !绘制边坡模型塑性应变云图
I~5fz4Q
!边坡在强度折减系数F=2.6时结果分析 V1fvQ= 9
Resume,'F2.6','db' !读入边坡在强度折减系数F=2.6时 !7a^8
set,1,last !读入后一个子步 _Ag/gu2-?
pldisp,1 !绘制边坡模型变形图 a]=vq(N'r
plnsol,u,x !绘制边坡模型水平方向位移云图 "i<3}6/*
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 -O> mY)
!边坡在强度折减系数F=2.8时结果分析 "d rh+oo.
Resume,'F2.8','db' !读入边坡在强度折减系数F=2.8时 08n%%
F
set,1,last !读入后一个子步 eYSGxc x
pldisp,1 !绘制边坡模型变形图 @p` CAB
plnsol,u,x !绘制边坡模型水平方向位移云图 !DUg"o3G>
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 !}Ou|r4_
!边坡在强度折减系数F=3.0时结果分析 *Wo$ $T
Resume,'F3.0','db' !读入边坡在强度折减系数F=3.0时 VOK$;s'9}
set,1,last !读入后一个子步 YU0HySP:
pldisp,1 !绘制边坡模型变形图 Y l3[~S
plnsol,u,x !绘制边坡模型水平方向位移云图 RDZh>K
PG
plnsol,eppl,eqv !绘制边坡模型塑性应变云图
