各位朋友： _L@2_#h!
大家好！我通过ansys用有限元强度折减法分析了一个边坡的稳定性，结果不同的工况下得出的塑性区是一样的，我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激！附命令流和塑性区如下：(1)创建物理环境 ET[k pL
/COM, Structural ! 指定结构分析 glDh( [
/TITLE, Slope stability Analysis ! 定义工作标题 bd<m%OM""
/FILNAM,support,1 ! 定义工作文件名 lhtZaU~V
(2)建立模型 CYKr\DA
!进入前处理器 JqMDqPIQ
/PREP7 bzvh%RsW
!定义单元类型 8
O 67
ET,1,PLANE82 !定义边坡围岩单元 vmLxkjUm#
KEYOPT,1,3,2 !求解类型为平面应变 ?g wUwOV"
!定义材料属性 :OW ;?{ ~j
!边坡围岩1材料属性 7{xh 8#m
MP,EX,1,3E10 lRv eHB&V
MP,PRXY,1,0.25 3 q`)*
MP,DENS,1,2500 !定义塑性模型 YSj+\Z$(
TB,DP,1 !定义数据表1 h68]=KyK
TBDATA,1,9E5,42.8 !给数据表输入数据 P(;?kg}0
!边坡围岩2材料属性 !采用弹性模型 I8QjKI (
MP,EX,2,3.2E10 ！定义弹性模量 > ]8a 3x
MP,PRXY,2,0.24 !定义泊松比 {U84 _Pi
MP,DENS,2,2700 !定义密度 HU9Sl*/
!取强度折减系数F=1.2时边坡围岩材料属性 {&"L ~>/o
MP,EX,3,3E10 $DZH QH
MP,PRXY,3,0.25 ?_ v{|
YI=
MP,DENS,3,2500 eYC ^4g%l(
TB,DP,3 !定义数据表3 uio@r^Xz
TBDATA,1,7.5E5,37.7 !给数据表输入数据 = +x yI
!取强度折减系数F=1.4时边坡围岩材料属性 Q7g>4GZC
MP,EX,4,3E10 11,!XD*"
MP,PRXY,4,0.25 !定义泊松比 JkGnKm9G
MP,DENS,4,2500 JKX_q&bUw
TB,DP,4 !定义数据表4 fPR$kc h
TBDATA,1,6.4E5,33.5 !给数据表输入数据 e5OsI Vtjr
!取强度折减系数F=1.6时边坡围岩材料属性 qVdwfT{1J
MP,EX,5,3E10 ]I L;`>Gp
MP,PRXY,5,0.25 4&HXkRs:
MP,DENS,5,2500 ！定义密度 fd #j Y}
TB,DP,5 !定义数据表5 (?H0+zws^
TBDATA,1,5.6E5,30 !给数据表输入数据 b ]7GmRekl
!取强度折减系数F=1.8时边坡围岩材料属性 PC qZNBN
MP,EX,6,3E10 ！定义弹性模量 <&!v1yR
MP,PRXY,6,0.25 OKHX)"j\\
MP,DENS,6,2500 E d:eGm }
TB,DP,6 !定义数据表6 Cn_$l>
TBDATA,1,5.0E5,27.2 5BR2?hO4
!取强度折减系数F=2.0时边坡围岩材料属性 jn: NYJv
MP,EX,7,3E10 o=RM-tR`v
MP,PRXY,7,0.25 !定义泊松比 9m-)Xdoy
MP,DENS,7,2500 cS. -7
TB,DP,7 !定义数据表7 1 &/FG(*/
TBDATA,1,4.5E5,24.8 DaH?@Q
!取强度折减系数F=2.2时边坡围岩材料属性 ~h" /Tce
MP,EX,8,3E10 4Td{;Y="yF
MP,PRXY,8,0.25 (o J9k[(
MP,DENS,8,2500 "7 v-` i
TB,DP,8 !定义数据表8 )~V4+*<
TBDATA,1,4.09E5,22.8 fRfn2jA)d
!取强度折减系数F=2.4时边坡围岩材料属性 &h~Xq^
MP,EX,9,3E10 )I\=BPo|B
MP,PRXY,9,0.25 .@"q$\
MP,DENS,9,2500 @ w>zF/
TB,DP,9 !定义数据表9 = 0-
$W5E
TBDATA,1,3.6E5,21.1 !给数据表输入数据 'SYo_!
!取强度折减系数F=2.6时边坡围岩材料属性 Z--@.IYoJ
MP,EX,10,3E10 ！定义弹性模量 @
a4/ELx
MP,PRXY,10,0.25 !定义泊松比 q\rC5gk >
MP,DENS,10,2500 ！定义密度 9
C {;h
TB,DP,10 !定义数据表10 OgcHS?
TBDATA,1,3.46E5,19.6 !给数据表10输入数据 ?DH"V7bs
!取强度折减系数F=2.8时边坡围岩材料属性 HAYMX:%
MP,EX,11,3E10 ！定义弹性模量 [I:KpAd/
MP,PRXY,11,0.25 !定义泊松比 huudBc
A[
MP,DENS,11,2500 2qi'g:qe
TB,DP,11 !定义数据表11 {}o>{&X
TBDATA,1,3.2E5,18.2 !给数据表11输入数据 P"lBB8\eku
!取强度折减系数F=3.0时边坡围岩材料属性 O'GG Ti]e
MP,EX,12,3E10 ！定义弹性模量 ,+;:3gRk9
MP,PRXY,12,0.25 !m(4F(!"h
MP,DENS,12,2500 ！定义密度 \*\ R1_+
TB,DP,12 !定义数据表12 h C=:q
TBDATA,1,3.0E5,17.1 !给数据表12输入数据 u4"SH(
!建立几何模型 ^>?CMcN4*
!创建边坡线模型线 F/1#l@q N
K,1,,,, ！ 创建关键点1 < p TTo
K,2,-800,,, Wy]^Ub gW
K,3,-800,-800,, Y :l dR
K,4,-800,-1200,, l/y]nw
K,5,1200,-1200,, Kd3?I5 t
K,6,1200,-800,, $ DBGLmw
K,7,1200,,, AOp/d(vx5i
K,8,1200,400,, $hR)i
K,9,400,400,, ^2o dr \
l,1,2 !连接关键点1、2生成直线L1 " 7d-z<^n
l,2,3 OquAql:
l,3,4 Gq#~vr
l,4,5 !T Ap+b
l,5,6 ki;UY~
l,6,7 zNSu
l,7,8 L5IbExjV
l,8,9 rC_*sx
r^
l,9,1 +fKtG]$
l,1,7 uE=$p)
l,3,6 >\<*4J$PZ
!创建边坡面模型 $ XBAZ<"hd
al,3,4,5,11 !由线L3、L4、L5、L11生成一个面A1 L?Yoh<
al,1,2,11,6,10 !由线L1、L2、L1、L6、L10生成一个面A2 fWGOP~0
al,7,8,9,10 !由线L7、L8、L9、L10生成一个面A3 CJaKn z
/PNUM, area, on !打开面号开关 H ftxS
!划分网格生成有限元模型 @^ ik[9^H
! 划分边坡围岩2网格 hIJ) MZU|
!设置网格份数 1x~U*vbh Q
lsel,s,,,3,5,2 !选择线L3和L5 Q 6 C-4ja
lesize,all,,,5 !把所选择线分为5段 5~XN>>hp
lsel,s,,,4,11,7 !选择线L4和L11 'j%F]C K
lesize,all,,,26 !把所选择线分为26段 'rb'7=z5
mat,2 !给边坡围岩2赋予2号材料特性
he+#Q 6
type,1 !采用单元类型1 r C[6lIP
mshkey,0 !设定自由网格划分 f DgD@YC D
mshape,0 !设定四边形网格划分 K 3GSOD>
amesh,1 ！划分面积A1 s,*c@1f?
! 划分边坡围岩1网格 1on'^8]0
!设置网格份数 m+$ @'TbP
lsel,s,,,2,6,4 !选择线L2和L6 I-NN29 Sk
lsel,a,,,1 !附加选择线L1 .wrL3z_
lesize,all,,,10 !把所选择线分为10段 }i!pL(8;
lsel,s,,,8,10,2 !选择线L8和L10 )< 4 _ :
lesize,all,,,16 !把所选择线分为16段 1U\$iy8}
lsel,s,,,7,9,2 1+y" i<3)
lesize,all,,,12 02JL *
mat,1 !给边坡围岩1赋予1号材料特性 /DS?}I.*]
type,1 % qsvtc`
mshkey,0 !设定自由网格划分 !wLg67X$
-
mshape,0 !设定四边形网格划分 Lb=W;9;
amesh,2 !划分面积A2 >n.z)ZJ
amesh,3 ！划分面积A3 oh#N
0
0X
allsel zGcqzYbuA
!保存网格模型 R@T6U: 1
save,Slope-grid.db j%Y\A ~DV
(3)施加约束和荷载 J#3{S]* v_
!给边坡模型施加约束 @nV5.r0W}B
!边坡两侧施加X方向约束 Q}2[hB
nsel,s,loc,x,-800 !选择X=-800线上所有节点 ZC^C
nsel,a,loc,x,1200 !选择X=1200线上所有节点
'3l$al:H^
d,all,ux !对所选择节点约束X方向位移 b 1^n KB
allsel pF=g||gS
!边坡底部施加约束 10sK]XI
nsel,s,loc,y,-1200 !选择Y=--1200线上所有节点 {#Q\z>
d,all,uy !对所选择节点约束X、Y方向位移 `kN #4p
d,all,ux u>T76,8|\
!施加重力加速度 SjcL#S($&Y
acel,,9.8 j;P+_Hfe/E
(4)求解 ? S=W&
/solu `kBnSi o~
!求解设置 |@ VF.)_
antype,static !设定为静力求解 GSVdb/+
nsubst,100 !设定最大子步数为100 eTE2J~\
pred,on !打开时间步长预测器 &b C}3D
nropt,full !设定牛顿-拉普森选项 Il ~01|3+m
nlgeom,on !打开大位移效果 :9d \Uj,
lnsrch,on !打开线性搜索 dX u { p
outres,all,all !输出所有项 \x\.
cnvtol,f,,0.005,2,0.5 !力收敛准则设定 /nB'kg[h\
cnvtol,u,,0.05,2,1 !位移收敛准则设定 -
HOnB =
!边坡在强度折减系数F=1时求解 Fb Sa ~uN
allsel ?KN:r E
solve !进行求解 \MYU<6{u
save,F1,db !把F=1时求解结果保存 '?/&n8J\
!边坡在强度折减系数F=1.2时求解 ~\_T5/I%
finish {/M\Q@j
/solu E)' 8U
allsel o$bUY7_
asel,s,area,,2,3,1 !选择面积A2、A3 B1T5f1;uY
MPCHG,3,all ! 把所选择单元材料号改为3 w^VSj%XH!
allsel Uxl(9 6
solve !进行求解 ,a(O`##Bn
save,F1.2,db !把F=1.2时求解结果保存 s^C;>
!边坡在强度折减系数F=1.4时求解 >2- F2E ,
finish =&NOHT>
/solu .;Y
x*]
allsel z& 'f/w8
asel,s,area,,2,3,1 !选择面积A2、A3 }ymvC
MPCHG,4,all !把所选择单元材料号改为4 K r3];(w{
allsel nmE5]Pcg
solve !进行求解 $i Jnxqn
save,F1.4,db !把F=1.4时求解结果保存 P1d,8~ ;
!边坡在强度折减系数F=1.6时求解 ILG?r9 x
finish 1Kc^m\
/solu ; h + q
allsel ?n*fy
asel,s,area,,2,3,1 !选择面积A2、A3 =:5<{J OG
MPCHG,5,all !把所选择单元材料号改为5 m 6so]xr
allsel dK # h<q1
solve X ]s"5ju|t
save,F1.6,db !把F=1.6时求解结果保存 4]BJ0+|mT
!边坡在强度折减系数F=1.8时求解 G,DOBA
finish ep?:;98|t
/solu $N \+,?
allsel wu^q`!ml
asel,s,area,,2,3,1 !选择面积A2、A3 ~/ c5 hyTx
MPCHG,6,all !把所选择单元材料号改为6 {!K;`I[] v
allsel (|0b7 |'T
solve NW ISS
save,F1.8,db !把F=1.8时求解结果保存 ,C!n}+27
!边坡在强度折减系数F=2.0时求解 xii$e
finish o>u!CL<
/solu .~<]HAwq
allsel )fCMITq.|
asel,s,area,,2,3,1 !选择面积A2、A3 M.k|bh8
MPCHG,7,all !把所选择单元材料号改为7 8eww7k^R
allsel R~c1)[[E
solve TzY!D *%z
save,F2.0,db !把F=2.0时求解结果保存 s`Vf+ l0
!边坡在强度折减系数F=2.2时求解 S Kx e3
finish +u#Sl)F
/solu 3:" &Z6t#
allsel %c0z) R~
asel,s,area,,2,3,1 !选择面积A2、A3 uYY=~o[
Tw
MPCHG,8,all !把所选择单元材料号改为8 >z5Oy
allsel /xF 9:r
solve p^Agh
save,F2.2,db !把F=2.2时求解结果保存 c0SX]4}
G
!边坡在强度折减系数F=2.4时求解 @,GjeF]!
finish }oV3EIH
/solu ;i}i5yv2
allsel gDM Ac/V`l
asel,s,area,,2,3,1 !选择面积A2、A3 +!G4tA$g
MPCHG,9,all !把所选择单元材料号改为9 +{'lZa
allsel h2AGEg'g2[
solve =J`gGDhGY-
save,F2.4,db !把F=2.4时求解结果保存 jB` 7T^bU
!边坡在强度折减系数F=2.6时求解 l- Fmn/ V
finish h k/ +
/solu c'xUJhEL
allsel <TL])@da
asel,s,area,,2,3,1 !选择面积A2、A3 _J` |<}?t;
MPCHG,10,all !把所选择单元材料号改为10 *J.c $1#h
allsel ` v>/
solve w}"!l G
save,F2.6,db !把F=2.6时求解结果保存 4YBf ~Pp
!边坡在强度折减系数F=2.8时求解 x }_r nf_
finish rAL1TU(vm
/solu R2SBhs,+R
allsel 3V"dG1?
asel,s,area,,2,3,1 !选择面积A2、A3 yg`E2 2
MPCHG,11,all !把所选择单元材料号改为11 sN;(/O
allsel IC\E , m
solve FF8WTuzB+
save,F2.8,db !把F=2.8时求解结果保存 EQ6l:[
!边坡在强度折减系数F=3.0时求解 2$iw/ r
finish Q X sfp
/solu &^AzIfX}Gw
allsel MkkA {p
asel,s,area,,2,3,1 !选择面积A2、A3 %bEGv:88s
MPCHG,12,all !把所选择单元材料号改为12 m\l51}xz
allsel Io2,% !D
solve +yt 6.L
save,F3.0,db !把F=3.0时求解结果保存 <\+Po<)3 j
(5）后处理 vY4WQbz(
/post1 ！进入后处理 tsB.oDMP
!边坡在强度折减系数F=1时结果分析 mf+K{y,L
Resume,'F1','db' !读入边坡在强度折减系数F=1时 N'I?fWN!;R
set,1,last !读入后一个子步 t; h+Cf4
pldisp,1 !绘制边坡模型变形图 gQu\[e%mVo
plnsol,u,x !绘制边坡模型水平方向位移云图 RA! x
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 sV]i /B
!边坡在强度折减系数F=1.2时结果分析 i?lX,9%
Resume,'F1.2','db' !读入边坡在强度折减系数F=1.2时 3O #~dFnp
set,1,last !读入后一个子步 ? Ozk^#H[
pldisp,1 !绘制边坡模型变形图 >3<&V{<K
plnsol,u,x !绘制边坡模型水平方向位移云图 "r:H5) !
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 B8`R(vu;
!边坡在强度折减系数F=1.4时结果分析 @ysJt
Resume,'F1.4','db' !读入边坡在强度折减系数F=1.4时 f S(^["*G
set,1,last !读入后一个子步 Ar\IZ_Q
pldisp,1 !绘制边坡模型变形图 E=$7ieW
plnsol,u,x !绘制边坡模型水平方向位移云图 ,oC= {^l{
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 Pw0{.W~r
!边坡在强度折减系数F=1.6时结果分析 H't `Q&]a
Resume,'F1.6','db' !读入边坡在强度折减系数F=1.6时 7Ntjx(b$"h
set,1,last !读入后一个子步 0 =3FO}[u
pldisp,1 !绘制边坡模型变形图 "(9=h@@Y"
plnsol,u,x !绘制边坡模型水平方向位移云图 W ]a7&S
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 8 h55$j
!边坡在强度折减系数F=1.8时结果分析 &z7N \n
Resume,'F1.8','db' !读入边坡在强度折减系数F=1.8时 cpP.7ZR
set,1,last !读入后一个子步 kxcgOjrmI
pldisp,1 !绘制边坡模型变形图 se 4w~\/
plnsol,u,x !绘制边坡模型水平方向位移云图 >^M! @=/?J
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 U@ 1#!ZZ6
!边坡在强度折减系数F=2.0时结果分析 B7NmET4
Resume,'F2.0','db' !读入边坡在强度折减系数F=2.0时 #p11D=
@[
set,1,last !读入后一个子步 B iVd
ka
pldisp,1 !绘制边坡模型变形图 v'mJ~tz
plnsol,u,x !绘制边坡模型水平方向位移云图 %#[r_QQ^
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 m0dFA <5-
!边坡在强度折减系数F=2.2时结果分析 i KP\/LR<n
Resume,'F2.2','db' !读入边坡在强度折减系数F=2.2时 ]L5Z=.z&
set,1,last !读入后一个子步 < P? 3GT/
pldisp,1 !绘制边坡模型变形图 'XC&BWJ
plnsol,u,x !绘制边坡模型水平方向位移云图 S-H-tFy\\
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 YzforM^F
!边坡在强度折减系数F=2.4时结果分析
4WBo ZJ
Resume,'F2.4','db' !读入边坡在强度折减系数F=2.4时 V\e13cL]
set,1,last !读入后一个子步 ' lQ
pldisp,1 !绘制边坡模型变形图 oyi7Y Rvwd
plnsol,u,x !绘制边坡模型水平方向位移云图 EZgxSQaPH
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 T%N~oa
!边坡在强度折减系数F=2.6时结果分析
75Q XkJu
Resume,'F2.6','db' !读入边坡在强度折减系数F=2.6时 e<: 4czh8
set,1,last !读入后一个子步 -oaG|
pldisp,1 !绘制边坡模型变形图 vV$hGS(f~
plnsol,u,x !绘制边坡模型水平方向位移云图 =R"Eb1
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 mR!&.R?
!边坡在强度折减系数F=2.8时结果分析 '#pMEVP
Resume,'F2.8','db' !读入边坡在强度折减系数F=2.8时 E@\d<c.
set,1,last !读入后一个子步 \4]zNV ~x
pldisp,1 !绘制边坡模型变形图 ^Yg|P&e(;
plnsol,u,x !绘制边坡模型水平方向位移云图 EU?&
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 %_%f# S
!边坡在强度折减系数F=3.0时结果分析 |wyJh"4!
Resume,'F3.0','db' !读入边坡在强度折减系数F=3.0时 yhJA;&}>
set,1,last !读入后一个子步 ;r0|_mnf
pldisp,1 !绘制边坡模型变形图 Y ;~~?[6
plnsol,u,x !绘制边坡模型水平方向位移云图 ILu0J`;}
plnsol,eppl,eqv !绘制边坡模型塑性应变云图
