各位朋友： |DoD.?v
大家好！我通过ansys用有限元强度折减法分析了一个边坡的稳定性，结果不同的工况下得出的塑性区是一样的，我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激！附命令流和塑性区如下：(1)创建物理环境 dxWG+S
/COM, Structural ! 指定结构分析 brt` oR
/TITLE, Slope stability Analysis ! 定义工作标题 4 = hz4(5a
/FILNAM,support,1 ! 定义工作文件名 x*F_XE1#M
(2)建立模型 zLV k7u{e
!进入前处理器 to 'O;f">n
/PREP7 G/}nwj\
!定义单元类型 jR,3 -JQ
ET,1,PLANE82 !定义边坡围岩单元 Ns `:=
KEYOPT,1,3,2 !求解类型为平面应变 \F),SL
!定义材料属性 e&XJK*Wf
!边坡围岩1材料属性 .^?Z3iA",
MP,EX,1,3E10 JuXuS
MP,PRXY,1,0.25 zzyD'n7D
MP,DENS,1,2500 !定义塑性模型 vtf`+q
TB,DP,1 !定义数据表1 uB3Yl =P
TBDATA,1,9E5,42.8 !给数据表输入数据 *
n>YS
!边坡围岩2材料属性 !采用弹性模型 .?-]+ -J?`
MP,EX,2,3.2E10 ！定义弹性模量 )'t&LWS~
MP,PRXY,2,0.24 !定义泊松比 \EeK <)4:
MP,DENS,2,2700 !定义密度 VxGR[kq$]
!取强度折减系数F=1.2时边坡围岩材料属性 o7s<G8;?
MP,EX,3,3E10 u=!n9W~"
MP,PRXY,3,0.25 EoHrXv
MP,DENS,3,2500 Y> 't)PK
TB,DP,3 !定义数据表3 C Z/:(sOJ
TBDATA,1,7.5E5,37.7 !给数据表输入数据 w<=-n ;2
!取强度折减系数F=1.4时边坡围岩材料属性 aboA9pwH
MP,EX,4,3E10 hk+"c^g:j<
MP,PRXY,4,0.25 !定义泊松比 NpxND0
MP,DENS,4,2500 @RVj~J.A
TB,DP,4 !定义数据表4 KH76Vts
TBDATA,1,6.4E5,33.5 !给数据表输入数据 ^o6&|q
!取强度折减系数F=1.6时边坡围岩材料属性 L@LT *M
MP,EX,5,3E10 9TOq A4
MP,PRXY,5,0.25 uz e5u \
MP,DENS,5,2500 ！定义密度 [5jXYqD=vj
TB,DP,5 !定义数据表5 #jR1ti)p
TBDATA,1,5.6E5,30 !给数据表输入数据 &<S]= \
!取强度折减系数F=1.8时边坡围岩材料属性 U10: @Wzh
MP,EX,6,3E10 ！定义弹性模量 H}&4#CQ'!
MP,PRXY,6,0.25 Qx}hiv/
MP,DENS,6,2500 -Mufo.Jz1o
TB,DP,6 !定义数据表6 2o9IP>#u
TBDATA,1,5.0E5,27.2 \"hP*DJ"
!取强度折减系数F=2.0时边坡围岩材料属性 '9q:gFO
MP,EX,7,3E10 ,zH\&D $>u
MP,PRXY,7,0.25 !定义泊松比 Z "g6z#L&
MP,DENS,7,2500 $iAd)2LT
TB,DP,7 !定义数据表7 %(n^re uP
TBDATA,1,4.5E5,24.8 ^Dw18gqr=@
!取强度折减系数F=2.2时边坡围岩材料属性 8AV G pL
MP,EX,8,3E10 q#!c6 lG
MP,PRXY,8,0.25 m^ [VM &%
MP,DENS,8,2500 dRnO5
7+{
TB,DP,8 !定义数据表8 0@PI=JZ%
TBDATA,1,4.09E5,22.8 Hqm1[G)
!取强度折减系数F=2.4时边坡围岩材料属性 c9 N5 c
MP,EX,9,3E10 g|V0[Hnq6
MP,PRXY,9,0.25 RiM!LX
MP,DENS,9,2500 p3z%Y$!Tm
TB,DP,9 !定义数据表9 [^E{Yz=8,
TBDATA,1,3.6E5,21.1 !给数据表输入数据 XYsU)(;j
!取强度折减系数F=2.6时边坡围岩材料属性 v?(9ZY]
MP,EX,10,3E10 ！定义弹性模量 ^B7C8YP
MP,PRXY,10,0.25 !定义泊松比 0w
]
pDj
MP,DENS,10,2500 ！定义密度 D<XRu4^;
TB,DP,10 !定义数据表10 ,O ]AB
TBDATA,1,3.46E5,19.6 !给数据表10输入数据 t:.ZvA3
!取强度折减系数F=2.8时边坡围岩材料属性 n!t][d/g+
MP,EX,11,3E10 ！定义弹性模量 ? ;)F_aHp
MP,PRXY,11,0.25 !定义泊松比 s~MC t|a
MP,DENS,11,2500 lC(g&(\{
TB,DP,11 !定义数据表11 05Go*Qv V
TBDATA,1,3.2E5,18.2 !给数据表11输入数据 .LHzaeJCX
!取强度折减系数F=3.0时边坡围岩材料属性 d<B=p&~
MP,EX,12,3E10 ！定义弹性模量 14;lB.$p
MP,PRXY,12,0.25 p*S;4+>#
MP,DENS,12,2500 ！定义密度 [I7([l1Wvd
TB,DP,12 !定义数据表12 DANndXQLH
TBDATA,1,3.0E5,17.1 !给数据表12输入数据 H?<c eK'e
!建立几何模型 n aB`@
!创建边坡线模型线 W8{zV_TBm
K,1,,,, ！ 创建关键点1 @je vY81)
K,2,-800,,, ds UY[X-<6
K,3,-800,-800,, ,}_uk]AQ
K,4,-800,-1200,, (dVrGa54
K,5,1200,-1200,, 1brKs- z
K,6,1200,-800,, ~`Xu 6+1o
K,7,1200,,, I &I
q
K,8,1200,400,, maTZNzy
K,9,400,400,, meNz0v e
l,1,2 !连接关键点1、2生成直线L1 &C7HG^;W9
l,2,3 4 Z<
l,3,4 `h|>;u
l,4,5 bv8GJ #
l,5,6 i7(\i2_P
l,6,7 G`r*)pdm
l,7,8 PFnq:G^L
l,8,9 v&*}O
l,9,1 hz _F^gF
l,1,7 Q.Ljz
Z
l,3,6 /Jc i 1o
!创建边坡面模型 gR:21*&cz
al,3,4,5,11 !由线L3、L4、L5、L11生成一个面A1 t.E3Fh!o
al,1,2,11,6,10 !由线L1、L2、L1、L6、L10生成一个面A2 esVZ2_eL
al,7,8,9,10 !由线L7、L8、L9、L10生成一个面A3 (N25.}8 Y
/PNUM, area, on !打开面号开关 9F"Q2^l'
!划分网格生成有限元模型 qo" _w%{
! 划分边坡围岩2网格
@s@ 67\
!设置网格份数 =c(3EI'w
lsel,s,,,3,5,2 !选择线L3和L5 <E.$4/T
lesize,all,,,5 !把所选择线分为5段 [v$NxmRu
lsel,s,,,4,11,7 !选择线L4和L11 :Y ;\1J<b1
lesize,all,,,26 !把所选择线分为26段 gC qQ~lWZ
mat,2 !给边坡围岩2赋予2号材料特性 vs~lyM/
type,1 !采用单元类型1 lKF<] 25
mshkey,0 !设定自由网格划分 %!e;sL~&
mshape,0 !设定四边形网格划分 gELb(Y\ak
amesh,1 ！划分面积A1 'uOzC"_yF
! 划分边坡围岩1网格 rCa2$#Z
!设置网格份数 &k2nt
lsel,s,,,2,6,4 !选择线L2和L6 uuMHD{}?}
lsel,a,,,1 !附加选择线L1 wk"zpI7L
lesize,all,,,10 !把所选择线分为10段 3V`.<
lsel,s,,,8,10,2 !选择线L8和L10 hPHrq{YZ
lesize,all,,,16 !把所选择线分为16段 @|GKNW#
lsel,s,,,7,9,2 ^ M4-O~
lesize,all,,,12 15 ^5y RXC
mat,1 !给边坡围岩1赋予1号材料特性 N>}2&'I
type,1 Ho:X.Z9A^
mshkey,0 !设定自由网格划分 Y2,\ WKa
mshape,0 !设定四边形网格划分 8^Ov.$rP
amesh,2 !划分面积A2 L/R ES
amesh,3 ！划分面积A3 /4Lmu+G4
allsel HGjGV] N5
!保存网格模型 f:-dw6a=s
save,Slope-grid.db 9&`e jeD
(3)施加约束和荷载 M :Aik&
!给边坡模型施加约束 2i3& 3oz]O
!边坡两侧施加X方向约束 c_r&)8
nsel,s,loc,x,-800 !选择X=-800线上所有节点 ;c_p a0L
nsel,a,loc,x,1200 !选择X=1200线上所有节点 z.7'yJIP#
d,all,ux !对所选择节点约束X方向位移 H^no&$2`1
allsel rX<gcntv
!边坡底部施加约束 MjHjL~ Tg
nsel,s,loc,y,-1200 !选择Y=--1200线上所有节点 TGJ\f
d,all,uy !对所选择节点约束X、Y方向位移 axW4 cS ?
d,all,ux )d |hIW]7(
!施加重力加速度 ?w[M{
acel,,9.8 *t+ E8)qL
(4)求解 MYb^ILz H3
/solu >O{/%(9
!求解设置 KVrK:W--p
antype,static !设定为静力求解 =':,oz^|
nsubst,100 !设定最大子步数为100 _~ipO 1*
pred,on !打开时间步长预测器 G$D6#/rR
nropt,full !设定牛顿-拉普森选项 {PCf' n
nlgeom,on !打开大位移效果 %aszZP
lnsrch,on !打开线性搜索 f1RfNiW.
outres,all,all !输出所有项 .{|AHW&0<
cnvtol,f,,0.005,2,0.5 !力收敛准则设定 > WW5A py[
cnvtol,u,,0.05,2,1 !位移收敛准则设定 tyLR_@i%%
!边坡在强度折减系数F=1时求解 5E+k}S]M$
allsel 7o4 vf~
solve !进行求解 -^JGa{9*
save,F1,db !把F=1时求解结果保存 3[jk}2R';p
!边坡在强度折减系数F=1.2时求解 |e@Bi#M[
finish Q b|.;_
/solu Um$a9S8b&
allsel Q4;br ?2H
asel,s,area,,2,3,1 !选择面积A2、A3 *ppb 4R;CW
MPCHG,3,all ! 把所选择单元材料号改为3 d"> Ya !W
allsel 6GD Uo} .
solve !进行求解 ] cv|A^
save,F1.2,db !把F=1.2时求解结果保存 D? w-uR%Y
!边坡在强度折减系数F=1.4时求解 Ht 'jm (
finish '*KP{"3\
/solu U&|$B |[
allsel Nv
iPrp>c
asel,s,area,,2,3,1 !选择面积A2、A3 CKTD27})
MPCHG,4,all !把所选择单元材料号改为4 _ z!0 ab
allsel ~"*;lT5KX
solve !进行求解 LEJ7. 82
save,F1.4,db !把F=1.4时求解结果保存 nAp7X-t
!边坡在强度折减系数F=1.6时求解 afZPju"-
finish Z0~ }'K
/solu kPVP+}cA
allsel ib OXh U
asel,s,area,,2,3,1 !选择面积A2、A3 RhJL`>W`
MPCHG,5,all !把所选择单元材料号改为5 Wsr #YNhx|
allsel GZN@MK*co
solve R<!WW9IM
save,F1.6,db !把F=1.6时求解结果保存 |7CH
!边坡在强度折减系数F=1.8时求解 lE54RX}e4
finish TLL. Ch|#Y
/solu h ,@x5q>g
allsel Ccocv>=Q&J
asel,s,area,,2,3,1 !选择面积A2、A3 &\C{,:[
MPCHG,6,all !把所选择单元材料号改为6 w^*jhvV%kW
allsel ^j2:fJOU#
solve L?^C\g6u]
save,F1.8,db !把F=1.8时求解结果保存 QsmG(1=
!边坡在强度折减系数F=2.0时求解 f-\l<o(
finish )W@H
/solu "*@iXJxv5
allsel _e_]$G/TM
asel,s,area,,2,3,1 !选择面积A2、A3 Oed &B
MPCHG,7,all !把所选择单元材料号改为7 "z }bgy
allsel pg~`NN
solve uREc9z `Q'
save,F2.0,db !把F=2.0时求解结果保存 *TuoC5
!边坡在强度折减系数F=2.2时求解 S)"5X)mq
finish ^yRCR] oT
/solu n#N<zC/
allsel b)u9#%Q
asel,s,area,,2,3,1 !选择面积A2、A3 u%O^hc fb
MPCHG,8,all !把所选择单元材料号改为8 O D" eB?
allsel b>%I=H %g
solve K@{jY\AZNx
save,F2.2,db !把F=2.2时求解结果保存 q\[31$i$
!边坡在强度折减系数F=2.4时求解 (D8'qx-M
finish !qH=l-7A
/solu na
FZ<'t>&
allsel U6=m4]~Z
asel,s,area,,2,3,1 !选择面积A2、A3 uF T\a=
MPCHG,9,all !把所选择单元材料号改为9 WII_s|YSt%
allsel f tq~AF
solve i[ BR"(
save,F2.4,db !把F=2.4时求解结果保存 }}qR~.[
!边坡在强度折减系数F=2.6时求解 [,5clR=F
finish pAmTwe
/solu _<.R \rX&
allsel e7L;{+XI
asel,s,area,,2,3,1 !选择面积A2、A3 e_ 6
i896
MPCHG,10,all !把所选择单元材料号改为10 7$GP#V1r/
allsel ej<z]{`05
solve #% e`OA (b
save,F2.6,db !把F=2.6时求解结果保存 d1j v>tu
!边坡在强度折减系数F=2.8时求解 O)5-6lm
finish PyH E >C%
/solu 1&Y P}sg)
allsel ]dDyz[NuvD
asel,s,area,,2,3,1 !选择面积A2、A3 Pr<?E [
MPCHG,11,all !把所选择单元材料号改为11 Us Qv!Cwu^
allsel gr=ke #
solve q0y ?$XS
save,F2.8,db !把F=2.8时求解结果保存 &F#X0h/m=
!边坡在强度折减系数F=3.0时求解 %h%r6EB1F
finish ]?)zH :2)
/solu "_ )
allsel wCvD4C.WH
asel,s,area,,2,3,1 !选择面积A2、A3 *OF7 {^~&
MPCHG,12,all !把所选择单元材料号改为12 %
wRJ"T`Tt
allsel `T9<}&=!
solve qrX6FI
save,F3.0,db !把F=3.0时求解结果保存 xs &vgel>
(5）后处理 Gz *U?R-T
/post1 ！进入后处理 F{T|lTl
!边坡在强度折减系数F=1时结果分析 bM^'q
Resume,'F1','db' !读入边坡在强度折减系数F=1时 {/7'uD\
H
set,1,last !读入后一个子步 v2@M,xbxF:
pldisp,1 !绘制边坡模型变形图 6}ax~wYct
plnsol,u,x !绘制边坡模型水平方向位移云图 l:@.D|(o3
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 cDTDim1F
!边坡在强度折减系数F=1.2时结果分析 & 8&WY1cU
Resume,'F1.2','db' !读入边坡在强度折减系数F=1.2时 0/K NXz
set,1,last !读入后一个子步 g1J]z<&
pldisp,1 !绘制边坡模型变形图 EoY#D '[
plnsol,u,x !绘制边坡模型水平方向位移云图 Vxgc|E^J
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 "Pys3=h
!边坡在强度折减系数F=1.4时结果分析 OClY ,@
Resume,'F1.4','db' !读入边坡在强度折减系数F=1.4时 >Ft jrEB
set,1,last !读入后一个子步 w-nkf
M~
pldisp,1 !绘制边坡模型变形图 SwP h-6
plnsol,u,x !绘制边坡模型水平方向位移云图 <1jiU %!w
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 # 3C A
!边坡在强度折减系数F=1.6时结果分析 E$G "R =
Resume,'F1.6','db' !读入边坡在强度折减系数F=1.6时 t%dPj8~
set,1,last !读入后一个子步 R1q04Zj{2
pldisp,1 !绘制边坡模型变形图 SyYa_=En
plnsol,u,x !绘制边坡模型水平方向位移云图 d&'6l"${
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 {{DW P-v4
!边坡在强度折减系数F=1.8时结果分析 *gT
TI;:
Resume,'F1.8','db' !读入边坡在强度折减系数F=1.8时 'ZDa *9nkF
set,1,last !读入后一个子步 \O/=g6w|t}
pldisp,1 !绘制边坡模型变形图 <b Ta88,)
plnsol,u,x !绘制边坡模型水平方向位移云图 mU& J,C
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 xU rfH$$!`
!边坡在强度折减系数F=2.0时结果分析 Y66 vJ<lM
Resume,'F2.0','db' !读入边坡在强度折减系数F=2.0时 AARhGx |L<
set,1,last !读入后一个子步 L7xTAFe
pldisp,1 !绘制边坡模型变形图 < g %xo"
plnsol,u,x !绘制边坡模型水平方向位移云图 E_ y h 9lk
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 2d I:],7
!边坡在强度折减系数F=2.2时结果分析 YT:5J %"
Resume,'F2.2','db' !读入边坡在强度折减系数F=2.2时 rz.`$b
set,1,last !读入后一个子步 3 S 5QqAm
pldisp,1 !绘制边坡模型变形图 Q'=!1^&
plnsol,u,x !绘制边坡模型水平方向位移云图 i([A8C_A
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 zR<jZwo]#
!边坡在强度折减系数F=2.4时结果分析 OQsH,'
Resume,'F2.4','db' !读入边坡在强度折减系数F=2.4时 03#_ (
set,1,last !读入后一个子步 1%=,J'AH
pldisp,1 !绘制边坡模型变形图 gVl#pVO`N
plnsol,u,x !绘制边坡模型水平方向位移云图 WD )[Ac[
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 nEn2!)$
!边坡在强度折减系数F=2.6时结果分析 6 % ,Q
Resume,'F2.6','db' !读入边坡在强度折减系数F=2.6时 ?4~lA
L1
set,1,last !读入后一个子步 d~_OWCg`
pldisp,1 !绘制边坡模型变形图 6a ,YxR\
plnsol,u,x !绘制边坡模型水平方向位移云图 tIvtiN6[|l
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 W@0(Y9jdg
!边坡在强度折减系数F=2.8时结果分析 ?NwFpSB2
Resume,'F2.8','db' !读入边坡在强度折减系数F=2.8时 <n`|zQ
set,1,last !读入后一个子步 ?QMs<
pldisp,1 !绘制边坡模型变形图 !{3pp
plnsol,u,x !绘制边坡模型水平方向位移云图 7xd}J(l
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 0 s 4j>
!边坡在强度折减系数F=3.0时结果分析 eeX)JC0A
Resume,'F3.0','db' !读入边坡在强度折减系数F=3.0时 $ ` ho+
set,1,last !读入后一个子步 n_+Iw,a'm
pldisp,1 !绘制边坡模型变形图 {"y 6l
plnsol,u,x !绘制边坡模型水平方向位移云图 0+H"$2/
plnsol,eppl,eqv !绘制边坡模型塑性应变云图
