各位朋友： )AoF-&,w
大家好！我通过ansys用有限元强度折减法分析了一个边坡的稳定性，结果不同的工况下得出的塑性区是一样的，我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激！附命令流和塑性区如下：(1)创建物理环境 \CDzVO 0^
/COM, Structural ! 指定结构分析 ? > 7SZiC`
/TITLE, Slope stability Analysis ! 定义工作标题 JMIS*njq^
/FILNAM,support,1 ! 定义工作文件名 `-QY<STTP9
(2)建立模型 3D*vNVI
!进入前处理器 Tyk\l>S
/PREP7 qp)a`'Pq
!定义单元类型 Q\ k|pg?
ET,1,PLANE82 !定义边坡围岩单元 T>\ r}p
KEYOPT,1,3,2 !求解类型为平面应变 6Lh fb\2?
!定义材料属性 LYvjqNC&4
!边坡围岩1材料属性 H%
"R _[+
MP,EX,1,3E10 NP?hoqeKs
MP,PRXY,1,0.25 "F fIq;
MP,DENS,1,2500 !定义塑性模型 Y,;$RV@g
TB,DP,1 !定义数据表1 8M9}os
TBDATA,1,9E5,42.8 !给数据表输入数据 TU(w>v
!边坡围岩2材料属性 !采用弹性模型 q{*4BL'
MP,EX,2,3.2E10 ！定义弹性模量 g (:%E
MP,PRXY,2,0.24 !定义泊松比 }:f
\!b
MP,DENS,2,2700 !定义密度 A(dWA e,
!取强度折减系数F=1.2时边坡围岩材料属性 =n7 3bm
MP,EX,3,3E10 uBXl ltU
MP,PRXY,3,0.25 7B|
#*IZe
MP,DENS,3,2500 &Y/Myh[P
TB,DP,3 !定义数据表3 y %4G[Dz
TBDATA,1,7.5E5,37.7 !给数据表输入数据 ~:M"JNcs
!取强度折减系数F=1.4时边坡围岩材料属性 JZ M: R
MP,EX,4,3E10 m\O|BMHn
MP,PRXY,4,0.25 !定义泊松比 1deNrmp %
MP,DENS,4,2500 qJ8@A}}8
TB,DP,4 !定义数据表4 Pk 6l*+"r<
TBDATA,1,6.4E5,33.5 !给数据表输入数据 +o7Np| Ou
!取强度折减系数F=1.6时边坡围岩材料属性 ~Vf+@_G8`
MP,EX,5,3E10 map#4\
MP,PRXY,5,0.25 G*x"drP
MP,DENS,5,2500 ！定义密度 JDA :)[;
TB,DP,5 !定义数据表5 "`'' eV3
TBDATA,1,5.6E5,30 !给数据表输入数据 n fU\l<
!取强度折减系数F=1.8时边坡围岩材料属性 \vR&-+8dk
MP,EX,6,3E10 ！定义弹性模量 ]J
aV +b'O
MP,PRXY,6,0.25 vn0}l6n3s
MP,DENS,6,2500 ]l.qp5eQ
TB,DP,6 !定义数据表6 8H})Dq%d 7
TBDATA,1,5.0E5,27.2 H*M )<"X
!取强度折减系数F=2.0时边坡围岩材料属性 Q\ AM]
U
MP,EX,7,3E10 {vL4 :K
MP,PRXY,7,0.25 !定义泊松比 ?-g/hXx;
MP,DENS,7,2500 00d<V:Aoy
TB,DP,7 !定义数据表7 Rl2*oOVz
TBDATA,1,4.5E5,24.8 F7r!zKXZ
!取强度折减系数F=2.2时边坡围岩材料属性 ;E{k+vkqy
MP,EX,8,3E10 !L)~*!+Gf
MP,PRXY,8,0.25 B]xZ
4 Y
MP,DENS,8,2500 zYW+Goz/C
TB,DP,8 !定义数据表8 rnV\O L
TBDATA,1,4.09E5,22.8 GV aIZh<
!取强度折减系数F=2.4时边坡围岩材料属性 C ~&~Ano,
MP,EX,9,3E10 9]DMHA@
MP,PRXY,9,0.25 t|ih{0
MP,DENS,9,2500 jhN]1t /\X
TB,DP,9 !定义数据表9 j~ )GZV
TBDATA,1,3.6E5,21.1 !给数据表输入数据 e t$VR:
!取强度折减系数F=2.6时边坡围岩材料属性 8LMO2Wyq
MP,EX,10,3E10 ！定义弹性模量 37S bF,G
MP,PRXY,10,0.25 !定义泊松比 !&k}YF
MP,DENS,10,2500 ！定义密度 ^lp#j;Df
TB,DP,10 !定义数据表10 S1B^FLe7X
TBDATA,1,3.46E5,19.6 !给数据表10输入数据 {M.OOEcIp
!取强度折减系数F=2.8时边坡围岩材料属性 <+tSTc4>r
MP,EX,11,3E10 ！定义弹性模量 *C*'J7
MP,PRXY,11,0.25 !定义泊松比 gq^j-!Q)Q<
MP,DENS,11,2500 `~hB-Z5dI
TB,DP,11 !定义数据表11 *h `P+_Q7
TBDATA,1,3.2E5,18.2 !给数据表11输入数据 {H5a.+-(bE
!取强度折减系数F=3.0时边坡围岩材料属性 ZB5:FtW4
MP,EX,12,3E10 ！定义弹性模量 XIAHUT5~J
MP,PRXY,12,0.25 aB N^J_
MP,DENS,12,2500 ！定义密度 *x;&fyR
TB,DP,12 !定义数据表12 vl "l
TBDATA,1,3.0E5,17.1 !给数据表12输入数据 4w\@D>@}H
!建立几何模型 K9*vWoP '
!创建边坡线模型线 _=wu>h&7
K,1,,,, ！ 创建关键点1 fXV+aZ
K,2,-800,,, ][jW2;A
K,3,-800,-800,, S {q c1qj
K,4,-800,-1200,, }KB z8M5
K,5,1200,-1200,, J ^y1=PM
K,6,1200,-800,, X &z|im' d
K,7,1200,,, 9uNkd2 #
K,8,1200,400,, Ju"* ;/
K,9,400,400,, h
"MiD
l,1,2 !连接关键点1、2生成直线L1 R8*z}xy{
l,2,3 j0XS12eM
l,3,4 w6RB |^
l,4,5 tTMYqg zUk
l,5,6 * 7<{Xbsj^
l,6,7 r\+0J`
l,7,8 Cm~Pn "K_]
l,8,9 NM ` 5hd{
l,9,1 Zt;dPYq>
l,1,7 r}-si^fo;
l,3,6
(SEE(G3 5
!创建边坡面模型 ? nLlZpZ2v
al,3,4,5,11 !由线L3、L4、L5、L11生成一个面A1 R7q\^Yzo
al,1,2,11,6,10 !由线L1、L2、L1、L6、L10生成一个面A2 *WHQ1geI8
al,7,8,9,10 !由线L7、L8、L9、L10生成一个面A3 j;GH|22
/PNUM, area, on !打开面号开关 Mx3MNX /
!划分网格生成有限元模型 )\q A[rTG
! 划分边坡围岩2网格 t==C dCl
!设置网格份数 X3 bPBv
lsel,s,,,3,5,2 !选择线L3和L5 Z4^O`yS9+
lesize,all,,,5 !把所选择线分为5段 #GJ{@C3H8Q
lsel,s,,,4,11,7 !选择线L4和L11 *t)Y@=k3>
lesize,all,,,26 !把所选择线分为26段 ko2T9NI:S
mat,2 !给边坡围岩2赋予2号材料特性 d3m!34 ml
type,1 !采用单元类型1 mM[!g'*
mshkey,0 !设定自由网格划分 _6tir'z
mshape,0 !设定四边形网格划分 1&,d,<
amesh,1 ！划分面积A1 EDl*UG83G
! 划分边坡围岩1网格 n0ZrgTVJ
!设置网格份数 R;G"LT
lsel,s,,,2,6,4 !选择线L2和L6 ,EE,W0/zzM
lsel,a,,,1 !附加选择线L1 /#lqv)s'
lesize,all,,,10 !把所选择线分为10段 8(6(,WwP}
lsel,s,,,8,10,2 !选择线L8和L10 uuD|%-Ng
lesize,all,,,16 !把所选择线分为16段 Q.]$t
2J
lsel,s,,,7,9,2 kE8\\ }B7
lesize,all,,,12 r+0<A.''a
mat,1 !给边坡围岩1赋予1号材料特性 #cnh
~O
type,1 THVF(M4v
mshkey,0 !设定自由网格划分 sbju3nvk
mshape,0 !设定四边形网格划分 Ug VLHwkvk
amesh,2 !划分面积A2 @ PaOQ@
amesh,3 ！划分面积A3 u/I|<NAC,
allsel dHzQAqb8J
!保存网格模型 SU,#:s(
save,Slope-grid.db k%|Sl>{Ir
(3)施加约束和荷载 ?ufX3yia
!给边坡模型施加约束 SJ4+s4!l
<
!边坡两侧施加X方向约束 `GBa 3
nsel,s,loc,x,-800 !选择X=-800线上所有节点 Q{:5gh
nsel,a,loc,x,1200 !选择X=1200线上所有节点 N<$dbqoT|
d,all,ux !对所选择节点约束X方向位移
W0 &x0
allsel A"V($:>U
!边坡底部施加约束 }Ai_peO0a
nsel,s,loc,y,-1200 !选择Y=--1200线上所有节点 x$:P; #
d,all,uy !对所选择节点约束X、Y方向位移 I*SrK Zb
d,all,ux jDV;tEY#^
!施加重力加速度 (2SmB`g
acel,,9.8 zc`gm~@
(4)求解 -](NMRqfN
/solu F=B[%4q`%
!求解设置 xa=Lu?t%<
antype,static !设定为静力求解 Fad.!%[
nsubst,100 !设定最大子步数为100 ~fht [S?@M
pred,on !打开时间步长预测器 ' 1 }ybSG
nropt,full !设定牛顿-拉普森选项 j Q"z\}Wf
nlgeom,on !打开大位移效果 n&n WY+GEo
lnsrch,on !打开线性搜索 UZV\ ]Y
outres,all,all !输出所有项 |*T`3@R;3
cnvtol,f,,0.005,2,0.5 !力收敛准则设定 *3Vic
cnvtol,u,,0.05,2,1 !位移收敛准则设定 P4VMGP
!边坡在强度折减系数F=1时求解 }MiEbLduN
allsel GBvgVX<
solve !进行求解 C3]"y7
save,F1,db !把F=1时求解结果保存 Qms ,kX
!边坡在强度折减系数F=1.2时求解 EL$"/ptE
finish w<P$)~6
/solu m-v0=+~&
allsel >x3 $Ld
asel,s,area,,2,3,1 !选择面积A2、A3 moOc
G3=9
MPCHG,3,all ! 把所选择单元材料号改为3 I&Z4?K
allsel ;usR=i36b
solve !进行求解 .,9e~6}
save,F1.2,db !把F=1.2时求解结果保存 JQ~[$OGH
!边坡在强度折减系数F=1.4时求解 }t FRl
finish 7 :S4 Ur
/solu }<9IH%sgF
allsel l(<=JU O;
asel,s,area,,2,3,1 !选择面积A2、A3 XzD+#+By
MPCHG,4,all !把所选择单元材料号改为4 wyi%!H
allsel %@~;PS3kd
solve !进行求解 d$gT,+|vu
save,F1.4,db !把F=1.4时求解结果保存 MV?#g -5
!边坡在强度折减系数F=1.6时求解 <Crbc$!OeX
finish Ifu[L&U
/solu k0(_0o
allsel (rE.ft5$9
asel,s,area,,2,3,1 !选择面积A2、A3 \+T U{vr
MPCHG,5,all !把所选择单元材料号改为5 K-2oSS56
allsel c]+uj q
solve $[xS>iuD
save,F1.6,db !把F=1.6时求解结果保存 JGRL&MG4
!边坡在强度折减系数F=1.8时求解 QUO '{;,
finish Q&@e,7]V+
/solu S11ME
allsel _)%Sz"g^Ix
asel,s,area,,2,3,1 !选择面积A2、A3 N87)rhXSo,
MPCHG,6,all !把所选择单元材料号改为6 ec sQshR
allsel +.2O Z3(
solve :AYhBhitC
save,F1.8,db !把F=1.8时求解结果保存 5kx-s6 `!
!边坡在强度折减系数F=2.0时求解 ^B!c L~S*I
finish crbph.0
/solu /7fD;H^*
allsel v 1VH&~e
asel,s,area,,2,3,1 !选择面积A2、A3 'o]k Op@q
MPCHG,7,all !把所选择单元材料号改为7 ) -^(Su(!
allsel o FS2*u
solve L &,&SDr
save,F2.0,db !把F=2.0时求解结果保存 m FgrT
!边坡在强度折减系数F=2.2时求解 0Lx,qZ'
finish R^?9 V=Y<T
/solu gski :C
allsel EpU}~vC9C
asel,s,area,,2,3,1 !选择面积A2、A3 %5 ovW<E:
MPCHG,8,all !把所选择单元材料号改为8 ] J^/`gc
allsel + usB$=kJ
solve {X EX0|TZ
save,F2.2,db !把F=2.2时求解结果保存 Mc9J Fzp
!边坡在强度折减系数F=2.4时求解 &Sdf0"
finish NA!?.zn
/solu 7tP
qez#
allsel wA2^ I70-
asel,s,area,,2,3,1 !选择面积A2、A3 ' +f(9/
MPCHG,9,all !把所选择单元材料号改为9 r7jh)Q;BbR
allsel k!"6mo@rd
solve FI" KJk'
save,F2.4,db !把F=2.4时求解结果保存 *_sSM+S
!边坡在强度折减系数F=2.6时求解 I|<`Er-;58
finish < #zd]t
/solu ;Sl%I+?
allsel VVw5)O1'
asel,s,area,,2,3,1 !选择面积A2、A3 >+9:31p
MPCHG,10,all !把所选择单元材料号改为10 0WSOA[R%[b
allsel gx.\&W b
solve #+ Y%Bxf
save,F2.6,db !把F=2.6时求解结果保存 Ei4Iv#Oi`
!边坡在强度折减系数F=2.8时求解 t" nxny9&
finish "BZL*hHq
/solu <<PXh&wu0
allsel AD~_n ^
asel,s,area,,2,3,1 !选择面积A2、A3 j;J`P H
MPCHG,11,all !把所选择单元材料号改为11 tTb fyI
allsel 3vHkhhYQ
solve O$7cN\Z
save,F2.8,db !把F=2.8时求解结果保存 "$.B@[iY@
!边坡在强度折减系数F=3.0时求解 %)w7t[A2D
finish ! NJGW
/solu "0Z5cQjg
allsel uhU'm@JZ
asel,s,area,,2,3,1 !选择面积A2、A3 5fy{!
MPCHG,12,all !把所选择单元材料号改为12 >mj WC) U
allsel `Y3\R#
solve k'NP +N<M
save,F3.0,db !把F=3.0时求解结果保存 G2&,R{L6w
(5）后处理 48_( 'z*>
/post1 ！进入后处理 {s0!hp
!边坡在强度折减系数F=1时结果分析 ?-'GbOr!
Resume,'F1','db' !读入边坡在强度折减系数F=1时 \olY)b[
set,1,last !读入后一个子步 bPIo9clq
pldisp,1 !绘制边坡模型变形图 2 ^qJ'<2]M
plnsol,u,x !绘制边坡模型水平方向位移云图 (6i4N2
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 m^ /s}WEqp
!边坡在强度折减系数F=1.2时结果分析 \A6 }=
Resume,'F1.2','db' !读入边坡在强度折减系数F=1.2时 ?CldcxM#
set,1,last !读入后一个子步 a- lF}P\
pldisp,1 !绘制边坡模型变形图
Wf~PP;
plnsol,u,x !绘制边坡模型水平方向位移云图 'vV+Wu#[
plnsol,eppl,eqv !绘制边坡模型塑性应变云图
aTk Mg
!边坡在强度折减系数F=1.4时结果分析 K j~!E
H"
Resume,'F1.4','db' !读入边坡在强度折减系数F=1.4时 Bkg./iP5x
set,1,last !读入后一个子步 !{CaW 4
pldisp,1 !绘制边坡模型变形图 )" ]( ?V
plnsol,u,x !绘制边坡模型水平方向位移云图 Rm}G4Pq
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 : o:Z
!边坡在强度折减系数F=1.6时结果分析 !+DhH2;)F
Resume,'F1.6','db' !读入边坡在强度折减系数F=1.6时 b#*"eZj
set,1,last !读入后一个子步 ,ZVhL* "
pldisp,1 !绘制边坡模型变形图
`)>}b 3
plnsol,u,x !绘制边坡模型水平方向位移云图 n(.L=VuXn
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 w,j;XPp
!边坡在强度折减系数F=1.8时结果分析 ts; ^,|h
Resume,'F1.8','db' !读入边坡在强度折减系数F=1.8时 NZ?dJ"eq7
set,1,last !读入后一个子步 [ #fz [U
pldisp,1 !绘制边坡模型变形图 "-\8Y>E
plnsol,u,x !绘制边坡模型水平方向位移云图 $RYOj{1
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 @ +a}O
!边坡在强度折减系数F=2.0时结果分析 _YN
C}PUU
Resume,'F2.0','db' !读入边坡在强度折减系数F=2.0时 [
~kS)
set,1,last !读入后一个子步 `/9I` < y
pldisp,1 !绘制边坡模型变形图 TPLv]$n
plnsol,u,x !绘制边坡模型水平方向位移云图 j=c< Lo`
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 >*\yEH9"
!边坡在强度折减系数F=2.2时结果分析 5=b6B=\*~
Resume,'F2.2','db' !读入边坡在强度折减系数F=2.2时 x`I"%pG
set,1,last !读入后一个子步 ND/oK M+?
pldisp,1 !绘制边坡模型变形图 YxGqQO36
plnsol,u,x !绘制边坡模型水平方向位移云图 )@<HG$#
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 $~\Tl:!#?
!边坡在强度折减系数F=2.4时结果分析 0H x'C^m72
Resume,'F2.4','db' !读入边坡在强度折减系数F=2.4时 K L*+gq0k
set,1,last !读入后一个子步 u>:j$@56
pldisp,1 !绘制边坡模型变形图 Kw2]J)TO
plnsol,u,x !绘制边坡模型水平方向位移云图 { R/e1-;
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 $j \jT
!边坡在强度折减系数F=2.6时结果分析 8 "A0@fNz
Resume,'F2.6','db' !读入边坡在强度折减系数F=2.6时 <sX_hIA^Fx
set,1,last !读入后一个子步 . .w$p-1
pldisp,1 !绘制边坡模型变形图 &*oljGt8
plnsol,u,x !绘制边坡模型水平方向位移云图 m}]QP\
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 h,6> ^A
!边坡在强度折减系数F=2.8时结果分析 _#uRKy<`N
Resume,'F2.8','db' !读入边坡在强度折减系数F=2.8时 -:~ z,F
set,1,last !读入后一个子步 c~}FYO$
pldisp,1 !绘制边坡模型变形图 ;9q$eK%d
plnsol,u,x !绘制边坡模型水平方向位移云图 f0YB y<a
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 ,@I\'os
!边坡在强度折减系数F=3.0时结果分析 MZqHL4<|
Resume,'F3.0','db' !读入边坡在强度折减系数F=3.0时 x
