各位朋友： 0>`69&;g|
大家好！我通过ansys用有限元强度折减法分析了一个边坡的稳定性，结果不同的工况下得出的塑性区是一样的，我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激！附命令流和塑性区如下：(1)创建物理环境 smU+:~
/COM, Structural ! 指定结构分析 rij[ZrJ
/TITLE, Slope stability Analysis ! 定义工作标题 |H 4'*NP"
/FILNAM,support,1 ! 定义工作文件名 T,fz/ 5w
(2)建立模型 iI1t
P
!进入前处理器 $3ILVT
/PREP7 9s)oC$\
!定义单元类型 1:t>}[Y
ET,1,PLANE82 !定义边坡围岩单元 `jHGN i
KEYOPT,1,3,2 !求解类型为平面应变 Qi61(lK
!定义材料属性 'FhnSNT(4=
!边坡围岩1材料属性 3C2 >
MP,EX,1,3E10 bsm,lx]bH^
MP,PRXY,1,0.25 &M!:,B
MP,DENS,1,2500 !定义塑性模型 qrkT7f
TB,DP,1 !定义数据表1 Q'VS]n
TBDATA,1,9E5,42.8 !给数据表输入数据 Cw}\t!* !
!边坡围岩2材料属性 !采用弹性模型 8\9EDgT
MP,EX,2,3.2E10 ！定义弹性模量 \) ; rOqh
MP,PRXY,2,0.24 !定义泊松比 8 f.L a
MP,DENS,2,2700 !定义密度 X@)lPr$a
!取强度折减系数F=1.2时边坡围岩材料属性 ?1uAY.~ZZB
MP,EX,3,3E10 2$91+N*w9
MP,PRXY,3,0.25 O2e "TH3
MP,DENS,3,2500 f/x "yUq
TB,DP,3 !定义数据表3 &g!/@*[Nhh
TBDATA,1,7.5E5,37.7 !给数据表输入数据 1 W u
!取强度折减系数F=1.4时边坡围岩材料属性 C 0%%@
2+
MP,EX,4,3E10 SMyg=B\x?7
MP,PRXY,4,0.25 !定义泊松比 ?2TH("hV$
MP,DENS,4,2500 ;k8}D*?8
TB,DP,4 !定义数据表4 i@*
^]'
TBDATA,1,6.4E5,33.5 !给数据表输入数据 }0(
Na
!取强度折减系数F=1.6时边坡围岩材料属性 9& j]
MP,EX,5,3E10 Z:_y,( 1Q
MP,PRXY,5,0.25 \abl|;fj
MP,DENS,5,2500 ！定义密度 ?zEF?LJoK
TB,DP,5 !定义数据表5 ^~ Sn{e sA
TBDATA,1,5.6E5,30 !给数据表输入数据 Q [{v U
!取强度折减系数F=1.8时边坡围岩材料属性 f+V':qz
MP,EX,6,3E10 ！定义弹性模量 F*4+7$E0B
MP,PRXY,6,0.25 7E95" B&w
MP,DENS,6,2500 E'G>'c W;x
TB,DP,6 !定义数据表6 R;o_ *
TBDATA,1,5.0E5,27.2 63Hk N4D4
!取强度折减系数F=2.0时边坡围岩材料属性 dc)Gk
MP,EX,7,3E10 {E/ TC%
MP,PRXY,7,0.25 !定义泊松比 X3[!x Mij
MP,DENS,7,2500 kXr%73s
TB,DP,7 !定义数据表7 :dzU]pk%0
TBDATA,1,4.5E5,24.8 bv/b<N@4?$
!取强度折减系数F=2.2时边坡围岩材料属性 1N+#(<x@,
MP,EX,8,3E10 TRZ^$<AG
MP,PRXY,8,0.25 ^n/uY94E)p
MP,DENS,8,2500 vF&b|V+,
TB,DP,8 !定义数据表8 l_ c?q"X
TBDATA,1,4.09E5,22.8 Y"qY@`
!取强度折减系数F=2.4时边坡围岩材料属性 I,eyL$x
MP,EX,9,3E10 |@BN+o;`Om
MP,PRXY,9,0.25 DtZm|~)a
MP,DENS,9,2500 .IXwa,
TB,DP,9 !定义数据表9 g&>Hy!v,
TBDATA,1,3.6E5,21.1 !给数据表输入数据 &; \v_5N6
!取强度折减系数F=2.6时边坡围岩材料属性 g8Z14'Ke
MP,EX,10,3E10 ！定义弹性模量 v,&2 !Zv
MP,PRXY,10,0.25 !定义泊松比 Eg*3**gTO
MP,DENS,10,2500 ！定义密度 qpl5n'qHUc
TB,DP,10 !定义数据表10 hgK=fHJ k
TBDATA,1,3.46E5,19.6 !给数据表10输入数据 w%;'uN_
!取强度折减系数F=2.8时边坡围岩材料属性 pQNFH)=nw
MP,EX,11,3E10 ！定义弹性模量 5[_8N{QC;
MP,PRXY,11,0.25 !定义泊松比 o__q)"^~-
MP,DENS,11,2500 lNb\^b
TB,DP,11 !定义数据表11 L
~w=O!
TBDATA,1,3.2E5,18.2 !给数据表11输入数据
={^#E?
!取强度折减系数F=3.0时边坡围岩材料属性 8; 8}Oq
MP,EX,12,3E10 ！定义弹性模量 oK6lCGM5
MP,PRXY,12,0.25 d3GK.8y_z
MP,DENS,12,2500 ！定义密度 V1.F`3h~
TB,DP,12 !定义数据表12 meR2"JN'
TBDATA,1,3.0E5,17.1 !给数据表12输入数据 )a\h5nQI)
!建立几何模型 ~|rkt` 8p
!创建边坡线模型线 +b+sQ<w?.
K,1,,,, ！ 创建关键点1 5WT\0]RUa
K,2,-800,,, N6;Z\\&0^q
K,3,-800,-800,, ' T]oV~H
K,4,-800,-1200,, j,XKu5w)Oi
K,5,1200,-1200,, `?x$J
6p
K,6,1200,-800,, {rZ "cUm
K,7,1200,,, }H=OVbQor
K,8,1200,400,, )^TQedF
K,9,400,400,, (Y( [^N q
l,1,2 !连接关键点1、2生成直线L1 PS6`o
l,2,3 9f\8oJQ
l,3,4 l i-YkaP
l,4,5 ^v-'=1ub?
l,5,6 O 0#Jl8
l,6,7 919g5f`
l,7,8 9f,:j
l,8,9 AX+d? M
l,9,1 YW<2:1A|
l,1,7 zwtsw [.
l,3,6 [Jo TWouNU
!创建边坡面模型 ]B4mm__
al,3,4,5,11 !由线L3、L4、L5、L11生成一个面A1 .TC
`\m V
al,1,2,11,6,10 !由线L1、L2、L1、L6、L10生成一个面A2 w>u Z$/
al,7,8,9,10 !由线L7、L8、L9、L10生成一个面A3 `Tj}4f
/PNUM, area, on !打开面号开关 >{a,]q*
!划分网格生成有限元模型 3;NRW+
! 划分边坡围岩2网格 L])w-
!设置网格份数 7VcVI? ?
lsel,s,,,3,5,2 !选择线L3和L5 jhv1 D' >6
lesize,all,,,5 !把所选择线分为5段 _GqE'VX
lsel,s,,,4,11,7 !选择线L4和L11
y6}): |
lesize,all,,,26 !把所选择线分为26段 \]xYV}(FO
mat,2 !给边坡围岩2赋予2号材料特性 ak `) >
type,1 !采用单元类型1 U[l%oLra
mshkey,0 !设定自由网格划分 $4
Uy3C+6
mshape,0 !设定四边形网格划分 ItADO'M
amesh,1 ！划分面积A1 !\1 W*6U8;
! 划分边坡围岩1网格 5H+k_U
!设置网格份数 E,7~k d~y`
lsel,s,,,2,6,4 !选择线L2和L6 lIg2iun[n
lsel,a,,,1 !附加选择线L1 l{9h8]^
lesize,all,,,10 !把所选择线分为10段 ;L2bC3
lsel,s,,,8,10,2 !选择线L8和L10 9#ft;c
lesize,all,,,16 !把所选择线分为16段 @'@6vC
lsel,s,,,7,9,2 I?>T"nV +'
lesize,all,,,12 pm\X*t}L
mat,1 !给边坡围岩1赋予1号材料特性 )\vHIXnfJ1
type,1 moR2iyO_
mshkey,0 !设定自由网格划分 yU,xcq~l
mshape,0 !设定四边形网格划分 s.' \&B[
amesh,2 !划分面积A2 p'~5[JR:
amesh,3 ！划分面积A3 p;$9W+H0
allsel x[UO1% _o-
!保存网格模型 : !3 y>bP)
save,Slope-grid.db <q2nZI^
(3)施加约束和荷载 I })t
!给边坡模型施加约束 <R>z;2c
!边坡两侧施加X方向约束 #~;8#!X
nsel,s,loc,x,-800 !选择X=-800线上所有节点 K/`RZ!
nsel,a,loc,x,1200 !选择X=1200线上所有节点 *GP_ut%
d,all,ux !对所选择节点约束X方向位移 .5;LL,S-
allsel GDp p`' \
!边坡底部施加约束 Jr)`shJ"
nsel,s,loc,y,-1200 !选择Y=--1200线上所有节点 a}fClI-u
d,all,uy !对所选择节点约束X、Y方向位移 YS~t d+*
d,all,ux +o]BjgG
!施加重力加速度 I!gj; a?R
acel,,9.8 5Q` n6 x|
(4)求解 9
w1ONw8v
/solu (JW?azU
!求解设置 9^ p{/Io
antype,static !设定为静力求解 -P>=WZu
nsubst,100 !设定最大子步数为100 |+- i'N 9
pred,on !打开时间步长预测器 /T )n5X
nropt,full !设定牛顿-拉普森选项 e"#D ){k#
nlgeom,on !打开大位移效果 acQN pT
lnsrch,on !打开线性搜索 4Z9 wzQ>
outres,all,all !输出所有项 4iqmi<[("
cnvtol,f,,0.005,2,0.5 !力收敛准则设定 nM-SDVF M
cnvtol,u,,0.05,2,1 !位移收敛准则设定 F |TMpH/
!边坡在强度折减系数F=1时求解 DWQQ615i
allsel "R@N|Qx'
solve !进行求解 e oSM@Isu
save,F1,db !把F=1时求解结果保存 mthl?,I|
!边坡在强度折减系数F=1.2时求解 |SKG 4_wGe
finish !m( L0YH
/solu z \>X[yNpA
allsel I^(#\vRW
asel,s,area,,2,3,1 !选择面积A2、A3 x9l0UD*+g
MPCHG,3,all ! 把所选择单元材料号改为3 Aq%^>YAp
allsel mo[<4U ks
solve !进行求解 NxO^VUD
save,F1.2,db !把F=1.2时求解结果保存 }Xj_Y]T
!边坡在强度折减系数F=1.4时求解 <0)ud)~ u
finish d~-p;i
/solu D|j \ nQ
allsel *)1Vs'!-
asel,s,area,,2,3,1 !选择面积A2、A3 u3m T
l
MPCHG,4,all !把所选择单元材料号改为4 g ypq`F
allsel -WvgK"k
solve !进行求解 ^alZ\!B8
save,F1.4,db !把F=1.4时求解结果保存 BWdc^
!边坡在强度折减系数F=1.6时求解 2F g t)`{!
finish GA.bRN2CI2
/solu + <9
e N
allsel orH0M!OtS!
asel,s,area,,2,3,1 !选择面积A2、A3 ,$ zl w\
MPCHG,5,all !把所选择单元材料号改为5 ApYud?0b
allsel H~ n~5 sF"
solve x ;,xd
save,F1.6,db !把F=1.6时求解结果保存 D1 ~ x
!边坡在强度折减系数F=1.8时求解 e=sJMzm~
finish aGb.
Lh9
/solu F*t_lN5{
allsel ~h0BT(p/
asel,s,area,,2,3,1 !选择面积A2、A3 Xj ~EVD
MPCHG,6,all !把所选择单元材料号改为6 ([b!$o<v
allsel L~y t AZ,
solve y*h1W4:^-
save,F1.8,db !把F=1.8时求解结果保存 'h> 5 &=r
!边坡在强度折减系数F=2.0时求解 V9u\;5oL
finish lc7a@qnw
/solu 9zYiG3 d
allsel bDBO+qA
asel,s,area,,2,3,1 !选择面积A2、A3 NjN ?RB/5
MPCHG,7,all !把所选择单元材料号改为7 /I'u/{KB
allsel L 8wcH
solve 9+
l3 $
save,F2.0,db !把F=2.0时求解结果保存 V}fKV6 v9
!边坡在强度折减系数F=2.2时求解 e~.?:7t
finish > '
0 ][~
/solu k_>Fw>Y
allsel 6h6?BQSE
asel,s,area,,2,3,1 !选择面积A2、A3 r/hyW6e_
MPCHG,8,all !把所选择单元材料号改为8 1] Q 2qs
allsel P3tG #cJ
solve #0hNk%X=
save,F2.2,db !把F=2.2时求解结果保存 U!?gdX
!边坡在强度折减系数F=2.4时求解 9< h]OXv
finish 5}bZs ` C
/solu ds;cfj[
allsel ;'QY<,p[e
asel,s,area,,2,3,1 !选择面积A2、A3 s%[GQQ-N
MPCHG,9,all !把所选择单元材料号改为9 e ]o'i;I
allsel UXPegK!
solve ~X< cG=p~u
save,F2.4,db !把F=2.4时求解结果保存 Wk#h ,p3
!边坡在强度折减系数F=2.6时求解 7[v@*/W@
finish o 7W Kh=
/solu t-*|Hfp*^
allsel Y%:0|utQC
asel,s,area,,2,3,1 !选择面积A2、A3 s^YTI\ L
\
MPCHG,10,all !把所选择单元材料号改为10 5b1uD>,;y
allsel aq+IC@O
solve rjHIQC C
save,F2.6,db !把F=2.6时求解结果保存 E\~ KVn
!边坡在强度折减系数F=2.8时求解 `k]!6osZo
finish ITIj=!F*
/solu E? eWv)//
allsel %M# ?cmt
asel,s,area,,2,3,1 !选择面积A2、A3 }?]yxa ~
MPCHG,11,all !把所选择单元材料号改为11 br o
allsel [~c'|E 8Q
solve 3'*%R48P`
save,F2.8,db !把F=2.8时求解结果保存 H9VXsFTW
!边坡在强度折减系数F=3.0时求解 IAF;mv}'
finish ()M@3={R
/solu Secq^#]8
allsel 7k =F6k0)
asel,s,area,,2,3,1 !选择面积A2、A3 5 k%9>U%$
MPCHG,12,all !把所选择单元材料号改为12 MiH}VfI
allsel S=H_9io
solve 6w"( y~c1
save,F3.0,db !把F=3.0时求解结果保存 FaE # \Q
(5）后处理 !Nl.Vb
/post1 ！进入后处理 DwmU fZp
!边坡在强度折减系数F=1时结果分析 M*|VLOo=v
Resume,'F1','db' !读入边坡在强度折减系数F=1时 wK'! xH^
set,1,last !读入后一个子步 }"?nU4q;S
pldisp,1 !绘制边坡模型变形图 OssR[$69
plnsol,u,x !绘制边坡模型水平方向位移云图 ) HX:U0
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 ,ah*!Zm.kk
!边坡在强度折减系数F=1.2时结果分析 (e>Rot0
Resume,'F1.2','db' !读入边坡在强度折减系数F=1.2时 fA_%8CjI
set,1,last !读入后一个子步 I+"?,Ej$K
pldisp,1 !绘制边坡模型变形图 ct}%Mdg
plnsol,u,x !绘制边坡模型水平方向位移云图
$.Q>M]xH
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 qJ+52U|z
!边坡在强度折减系数F=1.4时结果分析 q@!'R{fu
Resume,'F1.4','db' !读入边坡在强度折减系数F=1.4时 (;pi"/x [
set,1,last !读入后一个子步 "WbVCT'i
pldisp,1 !绘制边坡模型变形图 guOSO@
plnsol,u,x !绘制边坡模型水平方向位移云图 r}0C8(o q
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 MziZN^(
!边坡在强度折减系数F=1.6时结果分析 AR~$MCR]"k
Resume,'F1.6','db' !读入边坡在强度折减系数F=1.6时 Np<&#s[dQ
set,1,last !读入后一个子步 #CA%]*l*F
pldisp,1 !绘制边坡模型变形图 ur<eew@8@i
plnsol,u,x !绘制边坡模型水平方向位移云图 y (nsyA
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 P B (
!边坡在强度折减系数F=1.8时结果分析 VP %i1|XZJ
Resume,'F1.8','db' !读入边坡在强度折减系数F=1.8时 mPfUJ#rS
set,1,last !读入后一个子步 AwXt @!(
pldisp,1 !绘制边坡模型变形图 - *%!q$:
plnsol,u,x !绘制边坡模型水平方向位移云图 a,'Ncg
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 /MqXwUbO
!边坡在强度折减系数F=2.0时结果分析 {(z(NgXG/
Resume,'F2.0','db' !读入边坡在强度折减系数F=2.0时 /pN'K 5@
set,1,last !读入后一个子步 hkwa ""-
pldisp,1 !绘制边坡模型变形图 a We Bav}_
plnsol,u,x !绘制边坡模型水平方向位移云图 {!}F
:~*r
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 hzQ+9-qA
!边坡在强度折减系数F=2.2时结果分析 w^])(
Resume,'F2.2','db' !读入边坡在强度折减系数F=2.2时 /}$T38
set,1,last !读入后一个子步 3p3WDL7
pldisp,1 !绘制边坡模型变形图 tpa^k
plnsol,u,x !绘制边坡模型水平方向位移云图 2Za ,4'
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 J, 0pe\5
!边坡在强度折减系数F=2.4时结果分析 u-QO>3oY6
Resume,'F2.4','db' !读入边坡在强度折减系数F=2.4时 @>G&7r:U
set,1,last !读入后一个子步 Y62u%':X
pldisp,1 !绘制边坡模型变形图 >bA$SN
plnsol,u,x !绘制边坡模型水平方向位移云图 =,4
' "
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 ZEj!jWP2m
!边坡在强度折减系数F=2.6时结果分析 ;EJPrDHTk
Resume,'F2.6','db' !读入边坡在强度折减系数F=2.6时 /MKNv'5&!%
set,1,last !读入后一个子步 inPE/Ux
pldisp,1 !绘制边坡模型变形图 #PslrA.
E
plnsol,u,x !绘制边坡模型水平方向位移云图 wD6!#t k
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 ]A]Ft!`6 z
!边坡在强度折减系数F=2.8时结果分析 Lo9G4Cu
Resume,'F2.8','db' !读入边坡在强度折减系数F=2.8时 n^AP"1l8?0
set,1,last !读入后一个子步 z^rhgs?4
pldisp,1 !绘制边坡模型变形图 uKbHFF
plnsol,u,x !绘制边坡模型水平方向位移云图 vX WsF\g
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 ~mK-8U4>K,
!边坡在强度折减系数F=3.0时结果分析 slge+xq\J
Resume,'F3.0','db' !读入边坡在强度折减系数F=3.0时 +~
3w5.8
set,1,last !读入后一个子步 k%]DT.cE
pldisp,1 !绘制边坡模型变形图 UXwB$@8
plnsol,u,x !绘制边坡模型水平方向位移云图 _}zo
/kDA
plnsol,eppl,eqv !绘制边坡模型塑性应变云图
