各位朋友: 2B*9]AHny
大家好!我通过ansys用有限元强度折减法分析了一个边坡的稳定性,结果不同的工况下得出的塑性区是一样的,我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激!附命令流和塑性区如下:(1)创建物理环境 e0aeiG$/0
/COM, Structural ! 指定结构分析 _p90Zm-3X
/TITLE, Slope stability Analysis ! 定义工作标题 {Ynr(J.
/FILNAM,support,1 ! 定义工作文件名 Qn6&M
(2)建立模型 ;[*7UE+#7
!进入前处理器 bhFzu[B
/PREP7 0chBw~@*s
!定义单元类型 VpD9!;S
ET,1,PLANE82 !定义边坡围岩单元
9((v.
KEYOPT,1,3,2 !求解类型为平面应变 ./[t'dgC
!定义材料属性 F+}MW/ra@
!边坡围岩1材料属性 \.}* s]6
MP,EX,1,3E10 VdYu| w;v
MP,PRXY,1,0.25 9GCxF`OB
MP,DENS,1,2500 !定义塑性模型 TMq\}k-I5
TB,DP,1 !定义数据表1 i)eub`uMy
TBDATA,1,9E5,42.8 !给数据表输入数据 =0O`VSb
!边坡围岩2材料属性 !采用弹性模型 ?z?IEj}
MP,EX,2,3.2E10 !定义弹性模量 Pxl7zz&pl=
MP,PRXY,2,0.24 !定义泊松比 `L0}^|`9
MP,DENS,2,2700 !定义密度 ]"F0"UH,
!取强度折减系数F=1.2时边坡围岩材料属性 {E>kFeg
MP,EX,3,3E10 I=.98v%
MP,PRXY,3,0.25 W4X=.vr
MP,DENS,3,2500 w+[r$+z!k
TB,DP,3 !定义数据表3 RW%e%
TBDATA,1,7.5E5,37.7 !给数据表输入数据 /}(w{6C
!取强度折减系数F=1.4时边坡围岩材料属性 #knpZ'
MP,EX,4,3E10 r"k\G\,%
MP,PRXY,4,0.25 !定义泊松比 St'3e<
MP,DENS,4,2500 z0doLb^!
TB,DP,4 !定义数据表4 yn`H }@`k
TBDATA,1,6.4E5,33.5 !给数据表输入数据 XlR.Y~
!取强度折减系数F=1.6时边坡围岩材料属性 }Vk#w%EJ
MP,EX,5,3E10 ~-|K5
MP,PRXY,5,0.25 bT0CQ_g21
MP,DENS,5,2500 !定义密度 _0ep[r
TB,DP,5 !定义数据表5 >P/.X^G0
TBDATA,1,5.6E5,30 !给数据表输入数据 qm]k
(/w
!取强度折减系数F=1.8时边坡围岩材料属性 ^ WidA-
MP,EX,6,3E10 !定义弹性模量 @G^
l`%
MP,PRXY,6,0.25 yX<Sk q
MP,DENS,6,2500 n:' Mpux
TB,DP,6 !定义数据表6 c]#F^(-A`
TBDATA,1,5.0E5,27.2 uq{w1O5
!取强度折减系数F=2.0时边坡围岩材料属性 e")s1`
MP,EX,7,3E10 %Q]m6ciAM
MP,PRXY,7,0.25 !定义泊松比 \za 0?b
MP,DENS,7,2500 A6pPx1-&
TB,DP,7 !定义数据表7 t4~Bn<=
TBDATA,1,4.5E5,24.8 evz{@;.R
!取强度折减系数F=2.2时边坡围岩材料属性 6|~N5E~SX
MP,EX,8,3E10 FUiEayM
MP,PRXY,8,0.25 4f'1g1@$
MP,DENS,8,2500 Xw_AZ-|1D
TB,DP,8 !定义数据表8 e@w-4G(;
TBDATA,1,4.09E5,22.8
5Ta<$t
!取强度折减系数F=2.4时边坡围岩材料属性 "NMX>a,(
MP,EX,9,3E10 Q6]SsV?x
MP,PRXY,9,0.25 {OH"d
MP,DENS,9,2500 /X9K g
TB,DP,9 !定义数据表9 =0
mf
TBDATA,1,3.6E5,21.1 !给数据表输入数据 !^:b?M
!取强度折减系数F=2.6时边坡围岩材料属性 #2^eGhwnI
MP,EX,10,3E10 !定义弹性模量 sWa`-gc
MP,PRXY,10,0.25 !定义泊松比 criOJ-
MP,DENS,10,2500 !定义密度 zU}Ru&T9
TB,DP,10 !定义数据表10 h3[^uYe
TBDATA,1,3.46E5,19.6 !给数据表10输入数据 *@^9]$*$
!取强度折减系数F=2.8时边坡围岩材料属性 ER;?[!
MP,EX,11,3E10 !定义弹性模量 NI,i)OSEN
MP,PRXY,11,0.25 !定义泊松比 Gx,<|v
MP,DENS,11,2500 MZZEqsD5[
TB,DP,11 !定义数据表11 W+k SL{0
TBDATA,1,3.2E5,18.2 !给数据表11输入数据 3"!h+dXw
!取强度折减系数F=3.0时边坡围岩材料属性 h/?$~OD
MP,EX,12,3E10 !定义弹性模量 4,4S5u[|
MP,PRXY,12,0.25 Na=.LW-ma=
MP,DENS,12,2500 !定义密度 .hPk}B/KV
TB,DP,12 !定义数据表12 14Y_ oH9
TBDATA,1,3.0E5,17.1 !给数据表12输入数据 leR-oeSO
!建立几何模型 1Tm,#o
!创建边坡线模型线 $F2A
K,1,,,, ! 创建关键点1 waT'|9{
K,2,-800,,, `n)e]
dn
K,3,-800,-800,, |KY6IGcqV
K,4,-800,-1200,, -Ri/I4Xj
K,5,1200,-1200,, ObnQ,x(
K,6,1200,-800,, zZR_&z<
K,7,1200,,, O*jTrZ(k
K,8,1200,400,, X\sO eb:]
K,9,400,400,, !Pd@0n4
l,1,2 !连接关键点1、2生成直线L1 5e1;m6
l,2,3 vLCyT=OB`
l,3,4 ||yx?q6\h
l,4,5 @$mh0K>
l,5,6 &|MdBJ
l,6,7 GO4IAUA
l,7,8 &UrPb%=2H
l,8,9 ?gYQE&M !
l,9,1 r`-=<@[
l,1,7 ^hZZ5(</8P
l,3,6 ^1`Mz<
!创建边坡面模型 [L(qrAQ2|z
al,3,4,5,11 !由线L3、L4、L5、L11生成一个面A1 <UsFB F
al,1,2,11,6,10 !由线L1、L2、L1、L6、L10生成一个面A2 =^ZDP1h/}
al,7,8,9,10 !由线L7、L8、L9、L10生成一个面A3 H.o=4[
/PNUM, area, on !打开面号开关 9 n0?0mk
!划分网格生成有限元模型 Y.yiUf/Q
! 划分边坡围岩2网格 D`hl}
!设置网格份数 'wyS9^F
lsel,s,,,3,5,2 !选择线L3和L5 AcHr X=O
lesize,all,,,5 !把所选择线分为5段 @qUgp*+{
lsel,s,,,4,11,7 !选择线L4和L11 vd{ban9
lesize,all,,,26 !把所选择线分为26段 uSQlE=
mat,2 !给边坡围岩2赋予2号材料特性 "
l;=jk]
type,1 !采用单元类型1 ExKyjWAJ
mshkey,0 !设定自由网格划分 |-l9 Z
mshape,0 !设定四边形网格划分 ^ruz-N^Y!
amesh,1 !划分面积A1 ?*AhGza/
! 划分边坡围岩1网格 \RT3#X+
!设置网格份数
nQ +$
lsel,s,,,2,6,4 !选择线L2和L6 7H])2:)
lsel,a,,,1 !附加选择线L1 "50c<sZSB
lesize,all,,,10 !把所选择线分为10段 :;Xh`br
lsel,s,,,8,10,2 !选择线L8和L10 zQ~N(Jj?h
lesize,all,,,16 !把所选择线分为16段 tn38T%
lsel,s,,,7,9,2 utzf7?nIS
lesize,all,,,12 3-Bz5sj9
mat,1 !给边坡围岩1赋予1号材料特性 *q-VY[2
type,1 M!R=&a=Z
mshkey,0 !设定自由网格划分 X*w7q7\8-:
mshape,0 !设定四边形网格划分 l~Hu#+O
amesh,2 !划分面积A2 .+)
AeGh
amesh,3 !划分面积A3 `a6AES'w$
allsel _s18^7
!保存网格模型 ]WYddiF
save,Slope-grid.db u u$Jwn!S
(3)施加约束和荷载 )ko[_OJj
!给边坡模型施加约束 )k[{re
!边坡两侧施加X方向约束
MFyi#nq
nsel,s,loc,x,-800 !选择X=-800线上所有节点 {zIcEN$ ~
nsel,a,loc,x,1200 !选择X=1200线上所有节点 Og-v][
d,all,ux !对所选择节点约束X方向位移 ]bP1gV(b-
allsel }v xRjO,
!边坡底部施加约束 &|fPskpy
nsel,s,loc,y,-1200 !选择Y=--1200线上所有节点 }D]y-BbA.
d,all,uy !对所选择节点约束X、Y方向位移 laRn![[
d,all,ux U<|*V5
!施加重力加速度 qm/#kPlM
acel,,9.8 r XT6u
(4)求解 IOddu2.(
/solu Hb3..o:
!求解设置 MK.TBv
antype,static !设定为静力求解 "_0sW3rG
nsubst,100 !设定最大子步数为100 9\Md.>
pred,on !打开时间步长预测器 MX"M2>" pT
nropt,full !设定牛顿-拉普森选项 OUhlQq\
nlgeom,on !打开大位移效果 6 \?GY
lnsrch,on !打开线性搜索 $sFqMy
outres,all,all !输出所有项 O/.8;.d;4Y
cnvtol,f,,0.005,2,0.5 !力收敛准则设定 1*?XI
cnvtol,u,,0.05,2,1 !位移收敛准则设定 U!jRF
!边坡在强度折减系数F=1时求解 o'_eLp
allsel ]tNB^
solve !进行求解 cwE?+vB
save,F1,db !把F=1时求解结果保存 $}EI3a
!边坡在强度折减系数F=1.2时求解 e$`;z%6y
finish |RH^|2:x9Q
/solu /Yj; '\3
allsel ?aQVaw&L!7
asel,s,area,,2,3,1 !选择面积A2、A3 W'PW;.,
MPCHG,3,all ! 把所选择单元材料号改为3 vt#&YXu{A
allsel =d$m@rc0r
solve !进行求解 Et'&}NjI
save,F1.2,db !把F=1.2时求解结果保存 RF
[81/w]
!边坡在强度折减系数F=1.4时求解 I@kMM12>c
finish ~zoZ{YqP
/solu z(#CO<C.t
allsel B A
i ^t
asel,s,area,,2,3,1 !选择面积A2、A3 [5PQrf~Mo
MPCHG,4,all !把所选择单元材料号改为4 Wb5n> *
allsel 2;7n0LOs}
solve !进行求解 Bxf&gDwjgr
save,F1.4,db !把F=1.4时求解结果保存 nTE\EZ+=2
!边坡在强度折减系数F=1.6时求解 WM0-F@_
finish ^YB3$:@$U
/solu 8w ]'U
allsel ~X-v@a
asel,s,area,,2,3,1 !选择面积A2、A3 F!RP *
MPCHG,5,all !把所选择单元材料号改为5 -+0!Fkt@,
allsel #iT3aou
solve qgxGq(6K
save,F1.6,db !把F=1.6时求解结果保存 *iUR1V Y
!边坡在强度折减系数F=1.8时求解 C_ W%]8u
finish m' eM&1Ba
/solu 2UMX%+ "J
allsel :Vrj[i-{
asel,s,area,,2,3,1 !选择面积A2、A3 L4<=,}KS
MPCHG,6,all !把所选择单元材料号改为6 z
J V>;
allsel Dc3bG@K*G
solve 7Q|v5@;pU
save,F1.8,db !把F=1.8时求解结果保存 s,j=Kym%
!边坡在强度折减系数F=2.0时求解 +hIMfhF
finish acI%fYw5p`
/solu ,Qt2 ?
allsel ?:7.3{|Aq
asel,s,area,,2,3,1 !选择面积A2、A3 #;\tgUQ
MPCHG,7,all !把所选择单元材料号改为7 A<-3u
allsel 36d6KS 7
solve eMjW^-RgE5
save,F2.0,db !把F=2.0时求解结果保存 k^%F4d3z@C
!边坡在强度折减系数F=2.2时求解 .G>6_n3
finish AQs_(LR
/solu FXbalQ?^
allsel }q)dXFL=I#
asel,s,area,,2,3,1 !选择面积A2、A3 !Xzy:
MPCHG,8,all !把所选择单元材料号改为8 RCMO?CBe
allsel `8D'r|=`Eh
solve gFxa UrZA
save,F2.2,db !把F=2.2时求解结果保存 H,)2Ou-Wn
!边坡在强度折减系数F=2.4时求解 uXdR-@80*
finish 1ifPc5j}
/solu qbo
W<W<H1
allsel O*1la/~m
asel,s,area,,2,3,1 !选择面积A2、A3 xKXD`-|W
MPCHG,9,all !把所选择单元材料号改为9 6~Y`<#X5J
allsel AE4>pzBe
solve | EFbT>
save,F2.4,db !把F=2.4时求解结果保存 9U!#Y%*T
!边坡在强度折减系数F=2.6时求解 `IK3e9QpcA
finish b@[\+P] "
/solu rlkg.e6
allsel ZRYHsl{F+
asel,s,area,,2,3,1 !选择面积A2、A3 AAK}t6
MPCHG,10,all !把所选择单元材料号改为10 ETaLE[T%1
allsel %M-B"#OB7
solve AeEF/*
save,F2.6,db !把F=2.6时求解结果保存 T0BM:ofx
!边坡在强度折减系数F=2.8时求解 G=>LW1E|
finish #po}Y
/solu L)_L#]Yy
allsel !{4bC
asel,s,area,,2,3,1 !选择面积A2、A3 Q".AmHn
MPCHG,11,all !把所选择单元材料号改为11 UzUt=s!^H
allsel %|3NCyJ*7
solve R1\$}ep^
save,F2.8,db !把F=2.8时求解结果保存 #}zL?s^G
!边坡在强度折减系数F=3.0时求解 lR3JyYY{X
finish ipJnNy;
/solu d+g+{p>?
allsel H_7X%TvXb
asel,s,area,,2,3,1 !选择面积A2、A3 G `JXi/#`
MPCHG,12,all !把所选择单元材料号改为12 >=ng?
allsel sme!!+Rd
solve m>-(c=3
save,F3.0,db !把F=3.0时求解结果保存 g}'(V>(
(5)后处理 } @jT-t]P
/post1 !进入后处理 ,7izrf8
!边坡在强度折减系数F=1时结果分析 M2y"M ,k4
Resume,'F1','db' !读入边坡在强度折减系数F=1时 WAp#[mW.fx
set,1,last !读入后一个子步 ch]Q% M
pldisp,1 !绘制边坡模型变形图 X"
;ly0Mb
plnsol,u,x !绘制边坡模型水平方向位移云图 \B
D'"
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 hG.~[#[&6
!边坡在强度折减系数F=1.2时结果分析 %/=#8v4*
Resume,'F1.2','db' !读入边坡在强度折减系数F=1.2时 oF#]<Z\
set,1,last !读入后一个子步 fm'Qifq^
pldisp,1 !绘制边坡模型变形图 /(skIvE|
plnsol,u,x !绘制边坡模型水平方向位移云图 ){XG%nC
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 xnG,1doa
!边坡在强度折减系数F=1.4时结果分析 ,a?$F1Z-
Resume,'F1.4','db' !读入边坡在强度折减系数F=1.4时 @|J+f5O
set,1,last !读入后一个子步 @d=4C{g%o
pldisp,1 !绘制边坡模型变形图 r!J?Lc])8
plnsol,u,x !绘制边坡模型水平方向位移云图 |DMa2}%
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 i:,37INMt
!边坡在强度折减系数F=1.6时结果分析 [l X3":)
Resume,'F1.6','db' !读入边坡在强度折减系数F=1.6时 '}T6e1#JV
set,1,last !读入后一个子步 WjvD C"
pldisp,1 !绘制边坡模型变形图 qlfYX8edZ
plnsol,u,x !绘制边坡模型水平方向位移云图 zq&,KZ
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 ~F9WR5}]
!边坡在强度折减系数F=1.8时结果分析 e{33%5
Resume,'F1.8','db' !读入边坡在强度折减系数F=1.8时 fF} NPl
set,1,last !读入后一个子步 '74-rL:i
pldisp,1 !绘制边坡模型变形图 C/x<_VJzN/
plnsol,u,x !绘制边坡模型水平方向位移云图 j{u!/FD
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 ckPI^0A!
!边坡在强度折减系数F=2.0时结果分析 #35@YMF
Resume,'F2.0','db' !读入边坡在强度折减系数F=2.0时 o]|a5.O
set,1,last !读入后一个子步 #v`J]I)$
pldisp,1 !绘制边坡模型变形图 f6Io|CZWJ
plnsol,u,x !绘制边坡模型水平方向位移云图 b_T?jCyW
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 naoH685R4
!边坡在强度折减系数F=2.2时结果分析 ce\ F~8y
Resume,'F2.2','db' !读入边坡在强度折减系数F=2.2时 E>3fk
set,1,last !读入后一个子步 UP}5E h
pldisp,1 !绘制边坡模型变形图 C) "|sG
plnsol,u,x !绘制边坡模型水平方向位移云图 TGe{NUO
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 GP\Pk/E
!边坡在强度折减系数F=2.4时结果分析 ^UAL5}CQt
Resume,'F2.4','db' !读入边坡在强度折减系数F=2.4时 -O-qEQd
set,1,last !读入后一个子步 T5+iX`#M
pldisp,1 !绘制边坡模型变形图 x0!5z1KQh
plnsol,u,x !绘制边坡模型水平方向位移云图 @dj2#
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 (C
EXPf
!边坡在强度折减系数F=2.6时结果分析 uk~4R@=&H
Resume,'F2.6','db' !读入边坡在强度折减系数F=2.6时 -5 YvtL
set,1,last !读入后一个子步 XABI2Ex
pldisp,1 !绘制边坡模型变形图 =?/RaK/
w
plnsol,u,x !绘制边坡模型水平方向位移云图 :LwNOuavN
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 r{gJ[%
!边坡在强度折减系数F=2.8时结果分析 F!;0eS"xp
Resume,'F2.8','db' !读入边坡在强度折减系数F=2.8时 Z'~yUo=
set,1,last !读入后一个子步 a}+7MEUmZ/
pldisp,1 !绘制边坡模型变形图 N{<=s]I%x
plnsol,u,x !绘制边坡模型水平方向位移云图 /Ma"a
^
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 &k+'TcWm
!边坡在强度折减系数F=3.0时结果分析 ~"+Fp&[9f
Resume,'F3.0','db' !读入边坡在强度折减系数F=3.0时 {D
jz']
set,1,last !读入后一个子步 zp2IpYQ,3
pldisp,1 !绘制边坡模型变形图 007SA6xq
plnsol,u,x !绘制边坡模型水平方向位移云图 .V?i 3
plnsol,eppl,eqv !绘制边坡模型塑性应变云图