各位朋友： 3LEN~ N}
大家好！我通过ansys用有限元强度折减法分析了一个边坡的稳定性，结果不同的工况下得出的塑性区是一样的，我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激！附命令流和塑性区如下：(1)创建物理环境 @gZ%>qe
/COM, Structural ! 指定结构分析 KmMt:^9
/TITLE, Slope stability Analysis ! 定义工作标题 G",+jR ]
/FILNAM,support,1 ! 定义工作文件名 n44j]+ P
(2)建立模型 &oT]ycz%
!进入前处理器 2+| [e_
/PREP7 ;N)qNiJY
!定义单元类型 0hPm,H*Y]
ET,1,PLANE82 !定义边坡围岩单元 5 bgx;z9
KEYOPT,1,3,2 !求解类型为平面应变 n |(Y?`(
!定义材料属性 2^|*M@3r
!边坡围岩1材料属性 A0'Yfuie
MP,EX,1,3E10 _mcD* V
MP,PRXY,1,0.25 `/0FXb
8h
MP,DENS,1,2500 !定义塑性模型 -1 fT2e
TB,DP,1 !定义数据表1 AR&u9Y)I
TBDATA,1,9E5,42.8 !给数据表输入数据 ,#s}nJ4
!边坡围岩2材料属性 !采用弹性模型 y M>c**9
MP,EX,2,3.2E10 ！定义弹性模量 FQ );el'_V
MP,PRXY,2,0.24 !定义泊松比 :
G< 1
MP,DENS,2,2700 !定义密度 IuRKj8J)o
!取强度折减系数F=1.2时边坡围岩材料属性 e\\ I,
MP,EX,3,3E10 dD#A.C,Rz
MP,PRXY,3,0.25 w@hm>6j
MP,DENS,3,2500 `2GHB@S"k
TB,DP,3 !定义数据表3 *`|F?wF
TBDATA,1,7.5E5,37.7 !给数据表输入数据 :c
c#e&BO
!取强度折减系数F=1.4时边坡围岩材料属性 b8>r UGA{
MP,EX,4,3E10 s$mcIMqs
MP,PRXY,4,0.25 !定义泊松比 '#^ONn STn
MP,DENS,4,2500 C] { 43
TB,DP,4 !定义数据表4 ,*Sj7qb#
TBDATA,1,6.4E5,33.5 !给数据表输入数据 T'FRnC^~
!取强度折减系数F=1.6时边坡围岩材料属性 yY8q{\G
MP,EX,5,3E10 E{h
MP,PRXY,5,0.25 z~Gi/Ln
MP,DENS,5,2500 ！定义密度 Fz-Bd*uS
TB,DP,5 !定义数据表5 R\)pW9)
TBDATA,1,5.6E5,30 !给数据表输入数据 ]>&au8
!取强度折减系数F=1.8时边坡围岩材料属性 `ez_
{
MP,EX,6,3E10 ！定义弹性模量 I9H+ $Wjd
MP,PRXY,6,0.25 c}lUP(Ss
MP,DENS,6,2500 vWwp'q
TB,DP,6 !定义数据表6 $:SHZe
TBDATA,1,5.0E5,27.2 1)H+iN|im/
!取强度折减系数F=2.0时边坡围岩材料属性 =ejk E;
%L
MP,EX,7,3E10 `}"*i_0-5'
MP,PRXY,7,0.25 !定义泊松比 zuXJf+]
MP,DENS,7,2500 _r&`[@m
TB,DP,7 !定义数据表7 e5C560
TBDATA,1,4.5E5,24.8 N EJxd%-
!取强度折减系数F=2.2时边坡围岩材料属性 | M4_@P
MP,EX,8,3E10 d XrLeoK
MP,PRXY,8,0.25 Wfz&:J#
MP,DENS,8,2500 y?;&(Tcbt8
TB,DP,8 !定义数据表8 *JUP~/N r
TBDATA,1,4.09E5,22.8 < OCy
!取强度折减系数F=2.4时边坡围岩材料属性 b v~"_)C
MP,EX,9,3E10 cd#@"&r
MP,PRXY,9,0.25 vEk
jd#
MP,DENS,9,2500 V~ %!-7?
TB,DP,9 !定义数据表9 {|bf`
TBDATA,1,3.6E5,21.1 !给数据表输入数据 LDx1@a|83
!取强度折减系数F=2.6时边坡围岩材料属性 D!+d]A[r
MP,EX,10,3E10 ！定义弹性模量 QVsOB$
MP,PRXY,10,0.25 !定义泊松比 )u`q41!
MP,DENS,10,2500 ！定义密度 ]:8:|*w
TB,DP,10 !定义数据表10 w/#k .YE
TBDATA,1,3.46E5,19.6 !给数据表10输入数据 d@JjqE[
!取强度折减系数F=2.8时边坡围岩材料属性 QGs\af
MP,EX,11,3E10 ！定义弹性模量 >S,yqKp37~
MP,PRXY,11,0.25 !定义泊松比 a"4j9cO
MP,DENS,11,2500 &82Za %
TB,DP,11 !定义数据表11 QJ\
o"c
TBDATA,1,3.2E5,18.2 !给数据表11输入数据 tU.Y$%4
!取强度折减系数F=3.0时边坡围岩材料属性 $OldHe[p
MP,EX,12,3E10 ！定义弹性模量 &;DK^ta*P
MP,PRXY,12,0.25 }vgeQh-G
MP,DENS,12,2500 ！定义密度 |>Qj ]
TB,DP,12 !定义数据表12 bk44 qL;8
TBDATA,1,3.0E5,17.1 !给数据表12输入数据 [< Bk% B 5
!建立几何模型 Y/?V%X
!创建边坡线模型线 UOC>H%r~M?
K,1,,,, ！ 创建关键点1
^"STM'Zh
K,2,-800,,, uS`XWn<CSD
K,3,-800,-800,, w3WBgH
K,4,-800,-1200,, 7j| ^ZuI+
K,5,1200,-1200,, ZxF RE#y~2
K,6,1200,-800,, <F{ EZ Ii
K,7,1200,,, xp7 `[.
K,8,1200,400,, roz p
K,9,400,400,, | vL 0}e
l,1,2 !连接关键点1、2生成直线L1 )O*\ }6:S
l,2,3 4+"2K-]
l,3,4 GH[ATL
l,4,5 [|. IXdJ!
l,5,6 H0r@dn
l,6,7 z3x /Y/X$S
l,7,8 %38HGjS
l,8,9 )_B Q@5NK
l,9,1 0h=NbLr|S-
l,1,7 yq]= +X>(
l,3,6 jMr [ UZ
!创建边坡面模型 kCRf O}wt3
al,3,4,5,11 !由线L3、L4、L5、L11生成一个面A1 !Lkk1z o
al,1,2,11,6,10 !由线L1、L2、L1、L6、L10生成一个面A2 |Lf>Z2E
al,7,8,9,10 !由线L7、L8、L9、L10生成一个面A3 Pfi|RTX$'*
/PNUM, area, on !打开面号开关 >%LZ|*U
!划分网格生成有限元模型 q(xr5iuP_
! 划分边坡围岩2网格 DXA<m2&64N
!设置网格份数 F#Y9 @E
lsel,s,,,3,5,2 !选择线L3和L5 ]t0]fb[J
lesize,all,,,5 !把所选择线分为5段 4~<78r5m
lsel,s,,,4,11,7 !选择线L4和L11 U1 nObA
lesize,all,,,26 !把所选择线分为26段 c[VVCN8dA
mat,2 !给边坡围岩2赋予2号材料特性 0+S ;0
type,1 !采用单元类型1 6)=`&>9
mshkey,0 !设定自由网格划分 w]1hoYuV
mshape,0 !设定四边形网格划分 69< <pm,m
amesh,1 ！划分面积A1 FJ-X~^
! 划分边坡围岩1网格 <OgwA$abl%
!设置网格份数 Ql>bsr}
lsel,s,,,2,6,4 !选择线L2和L6 %[S-"k
lsel,a,,,1 !附加选择线L1 &FrUj>i
lesize,all,,,10 !把所选择线分为10段 |Yb]@9 >vn
lsel,s,,,8,10,2 !选择线L8和L10 ,i RUR 8
lesize,all,,,16 !把所选择线分为16段 *{y({J
lsel,s,,,7,9,2 ;Qn)~b~
lesize,all,,,12 )O+V ft&#
mat,1 !给边坡围岩1赋予1号材料特性 xB4}9zN s
type,1 nIZ;N!r=i
mshkey,0 !设定自由网格划分 0nr 5(4h
mshape,0 !设定四边形网格划分 POXd ,ON9
amesh,2 !划分面积A2 pTeN[Yu?
amesh,3 ！划分面积A3 @=<B8VPJd
allsel d!57`bVOd
!保存网格模型 Q&5s,)w-
save,Slope-grid.db x qSoE[<v
(3)施加约束和荷载 ~cv322N
!给边坡模型施加约束 mY= Q#nG
!边坡两侧施加X方向约束 M:M"7>:
nsel,s,loc,x,-800 !选择X=-800线上所有节点 _w}l,
nsel,a,loc,x,1200 !选择X=1200线上所有节点 MD9 8N{+[|
d,all,ux !对所选择节点约束X方向位移 mP*Ct6628n
allsel PQj 'D <G
!边坡底部施加约束 3QI. |;X
nsel,s,loc,y,-1200 !选择Y=--1200线上所有节点 i2P:I A|@
d,all,uy !对所选择节点约束X、Y方向位移 Talmc|h
d,all,ux H|)F-aL[
!施加重力加速度 I3qT SX-
acel,,9.8 ;~,)6UX 7
(4)求解 P9
w);jp;
/solu j
D kBe-`
!求解设置 Nn#;Kjul.
antype,static !设定为静力求解 yQcIfl]f
nsubst,100 !设定最大子步数为100 2WK c;?
pred,on !打开时间步长预测器 DuQW?9^232
nropt,full !设定牛顿-拉普森选项 3ncN) E/@
nlgeom,on !打开大位移效果 XjXz#0nR
lnsrch,on !打开线性搜索 nV:RL|p2jw
outres,all,all !输出所有项 t
i&!_
cnvtol,f,,0.005,2,0.5 !力收敛准则设定 b($9gre>mI
cnvtol,u,,0.05,2,1 !位移收敛准则设定 H-,p.$3}
!边坡在强度折减系数F=1时求解 dL]wu! wE
allsel 7v&>d,
solve !进行求解 !aB~G} '
save,F1,db !把F=1时求解结果保存 yGY:EvH^?
!边坡在强度折减系数F=1.2时求解 M3G ecjR
finish vw6 >eT
/solu ~KQiNkA\|l
allsel B 3 |G&Kg
asel,s,area,,2,3,1 !选择面积A2、A3 Q7#t #XM
MPCHG,3,all ! 把所选择单元材料号改为3 [*J?TNk
allsel dY{qdQQ}
solve !进行求解 `m thzc3W
save,F1.2,db !把F=1.2时求解结果保存 !bIE% cq
!边坡在强度折减系数F=1.4时求解 Mt4*`CxtH;
finish EK& ";(x2(
/solu
RnSll-
allsel *C);IdhK%y
asel,s,area,,2,3,1 !选择面积A2、A3 /_X`i[
MPCHG,4,all !把所选择单元材料号改为4 bc gXpP
allsel Zi?:< H}
solve !进行求解 ,8.$!Zia
save,F1.4,db !把F=1.4时求解结果保存 v|MT^.
!边坡在强度折减系数F=1.6时求解 j&u/T
finish LGK&& srJs
/solu | |L^yI~_d
allsel fK _uuw4
asel,s,area,,2,3,1 !选择面积A2、A3 *xjP^y":
MPCHG,5,all !把所选择单元材料号改为5 3 a|pk4M
allsel v\@pZw=x
solve :='I>Gn
save,F1.6,db !把F=1.6时求解结果保存 dT0>\9ZNr
!边坡在强度折减系数F=1.8时求解 nYR#Q|
finish l"1*0jgBw
/solu O!t=, F1j
allsel _p mo
6O
asel,s,area,,2,3,1 !选择面积A2、A3 R$p(5>#\5
MPCHG,6,all !把所选择单元材料号改为6 /'/i? 9:
allsel &e78xtA {
solve on;>iKta9
save,F1.8,db !把F=1.8时求解结果保存 $<9u:.9xf
!边坡在强度折减系数F=2.0时求解 -
~4na{6x
finish JZK93R
/solu S['cX ~
allsel / ykc`E?f
asel,s,area,,2,3,1 !选择面积A2、A3 1?yj<^"
MPCHG,7,all !把所选择单元材料号改为7 m[z $y
allsel +v~x_E5FP
solve qyAnq%B}
save,F2.0,db !把F=2.0时求解结果保存 a`8]TD
!边坡在强度折减系数F=2.2时求解 sjOv!|]A
finish G3 |x%/Fbp
/solu UM`{V5NG#
allsel O c.fvP^ZD
asel,s,area,,2,3,1 !选择面积A2、A3 puLgc$?
MPCHG,8,all !把所选择单元材料号改为8 1]9w9! j
allsel -k@1# c+z
solve EDuH+/:n
save,F2.2,db !把F=2.2时求解结果保存 w5^k84vye
!边坡在强度折减系数F=2.4时求解 +h Ks
finish ,
@!X! L
/solu I:HrBhI)wP
allsel Dw.I<fns^B
asel,s,area,,2,3,1 !选择面积A2、A3 "h #/b}/
MPCHG,9,all !把所选择单元材料号改为9 )&O6d .
allsel [?hv x}
solve 9 ^9-\DG
save,F2.4,db !把F=2.4时求解结果保存 ~ g-(
!边坡在强度折减系数F=2.6时求解 0b/@QgJ
finish LF
@_|o I
/solu Alo L+eN@
allsel alB' l
asel,s,area,,2,3,1 !选择面积A2、A3 0w:
3/WO
MPCHG,10,all !把所选择单元材料号改为10 `qYiic%
allsel wdAKU +tM
solve (w{T[ ~6
save,F2.6,db !把F=2.6时求解结果保存 W9pY =9]p+
!边坡在强度折减系数F=2.8时求解 ,Tu.cg
finish ;c>"gW8
/solu -EE'xh-zD
allsel w@&z0ODJ
asel,s,area,,2,3,1 !选择面积A2、A3 Y 9|!= T%
MPCHG,11,all !把所选择单元材料号改为11 jf- XVk5q
allsel o& &`_"18
solve Yku6\/^
save,F2.8,db !把F=2.8时求解结果保存 \|9B:y'y
!边坡在强度折减系数F=3.0时求解 NJ+$3n om
finish _"Z?O)d*
/solu +7o1&D*v
allsel (9oo8& GG
asel,s,area,,2,3,1 !选择面积A2、A3 p"c6d'qe
MPCHG,12,all !把所选择单元材料号改为12 s9a` 2Wm
allsel H la?\
solve Zs{7km
save,F3.0,db !把F=3.0时求解结果保存 [k>{q+MWK
(5）后处理 1Ml<>
/post1 ！进入后处理 FZn1$_Svr
!边坡在强度折减系数F=1时结果分析 &6C]| 13;
Resume,'F1','db' !读入边坡在强度折减系数F=1时 9 " t; 6
set,1,last !读入后一个子步 -seLa(8F
pldisp,1 !绘制边坡模型变形图 6)ibXbH
plnsol,u,x !绘制边坡模型水平方向位移云图 VBQAkl?(}4
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 Xz^k.4 Y{4
!边坡在强度折减系数F=1.2时结果分析 jrFPd
Resume,'F1.2','db' !读入边坡在强度折减系数F=1.2时 k(pJVez
set,1,last !读入后一个子步 n@ SUu7o
pldisp,1 !绘制边坡模型变形图 05k'TqT{c
plnsol,u,x !绘制边坡模型水平方向位移云图 >\^oCbqF}~
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 Z{BK@Q4z
!边坡在强度折减系数F=1.4时结果分析 @--"u_[
Resume,'F1.4','db' !读入边坡在强度折减系数F=1.4时 @'}2xw[eU
set,1,last !读入后一个子步 =.;ib6M
pldisp,1 !绘制边坡模型变形图 C4$P#D ZT^
plnsol,u,x !绘制边坡模型水平方向位移云图 S*=^I2;
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 l^ay* H
!边坡在强度折减系数F=1.6时结果分析 O|+ZEBP
Resume,'F1.6','db' !读入边坡在强度折减系数F=1.6时 =1,!EkG
set,1,last !读入后一个子步 i+mU(/l2{
pldisp,1 !绘制边坡模型变形图 JZ`SV}\`
plnsol,u,x !绘制边坡模型水平方向位移云图 Nwt [)\W `
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 &gruYZGK
!边坡在强度折减系数F=1.8时结果分析 i8k} B
o
Resume,'F1.8','db' !读入边坡在强度折减系数F=1.8时 ]T'8O`
set,1,last !读入后一个子步 :0Jn`Ds4o
pldisp,1 !绘制边坡模型变形图 S+~;PmN9qL
plnsol,u,x !绘制边坡模型水平方向位移云图 MymsDdQ]
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 eu" m0Q
!边坡在强度折减系数F=2.0时结果分析 i_Q1\_m !
Resume,'F2.0','db' !读入边坡在强度折减系数F=2.0时 p@% Pdx
set,1,last !读入后一个子步 lAM)X&}0
pldisp,1 !绘制边坡模型变形图 7Z`4K dh .
plnsol,u,x !绘制边坡模型水平方向位移云图 |v%$Q/zp&
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 -rI7ihr*
!边坡在强度折减系数F=2.2时结果分析 8tR6.09'
Resume,'F2.2','db' !读入边坡在强度折减系数F=2.2时 C_/eNu\I
set,1,last !读入后一个子步 SQ> Yf\
pldisp,1 !绘制边坡模型变形图 &4|]VOf
plnsol,u,x !绘制边坡模型水平方向位移云图 .<> t2,Af
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 /klo),|&
!边坡在强度折减系数F=2.4时结果分析 zA6C{L G3
Resume,'F2.4','db' !读入边坡在强度折减系数F=2.4时 =23B9WT
set,1,last !读入后一个子步 7F>]zrbK
pldisp,1 !绘制边坡模型变形图 >^`# %$+
plnsol,u,x !绘制边坡模型水平方向位移云图 [&*irk
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 Z0zEX?2mb
!边坡在强度折减系数F=2.6时结果分析 JUA%l
Resume,'F2.6','db' !读入边坡在强度折减系数F=2.6时 9ThsR&h3
set,1,last !读入后一个子步 Cj 3Xp~
pldisp,1 !绘制边坡模型变形图 (*\&xR Y|C
plnsol,u,x !绘制边坡模型水平方向位移云图 fJ.=,9:<
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 >S` =~4
!边坡在强度折减系数F=2.8时结果分析 #aC&!Rei{
Resume,'F2.8','db' !读入边坡在强度折减系数F=2.8时 6wB
! dl
set,1,last !读入后一个子步 cJq {;~
pldisp,1 !绘制边坡模型变形图 *!"T^ 4DEg
plnsol,u,x !绘制边坡模型水平方向位移云图 M5+R8ttc
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 }0E@eL
!边坡在强度折减系数F=3.0时结果分析 85io %>&0
Resume,'F3.0','db' !读入边坡在强度折减系数F=3.0时 4`#%<G
set,1,last !读入后一个子步 2i!R>`
pldisp,1 !绘制边坡模型变形图 2& Hl
wpx
plnsol,u,x !绘制边坡模型水平方向位移云图 x C _3&.
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 2 N &B
;WN% tI)
bt=D<YZk
w/o8R3 F
tQ; Fgv8Y!
u"v$[8
|AvsT{2
