各位朋友： #;2Ju'e#z
大家好！我通过ansys用有限元强度折减法分析了一个边坡的稳定性，结果不同的工况下得出的塑性区是一样的，我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激！附命令流和塑性区如下：(1)创建物理环境 oIX]9~
/COM, Structural ! 指定结构分析 }1 Q]C"hY
/TITLE, Slope stability Analysis ! 定义工作标题 zFwp$K>{QY
/FILNAM,support,1 ! 定义工作文件名 WG(t t.
(2)建立模型 n T\ W|
!进入前处理器 :`oY D
/PREP7 glPOW
!定义单元类型 ~@TNVkw
ET,1,PLANE82 !定义边坡围岩单元 kS3wa3bT
KEYOPT,1,3,2 !求解类型为平面应变 O`~T:N| D
!定义材料属性 wnha c}
!边坡围岩1材料属性 \q%li)
MP,EX,1,3E10 }U%2)M
MP,PRXY,1,0.25 ^;Ap-2Ww
MP,DENS,1,2500 !定义塑性模型 /]iv9e{uh(
TB,DP,1 !定义数据表1 '$4o,GA8
TBDATA,1,9E5,42.8 !给数据表输入数据 P\N$TYeH
!边坡围岩2材料属性 !采用弹性模型 B9Y " J
MP,EX,2,3.2E10 ！定义弹性模量 "4}{Z)&R2
MP,PRXY,2,0.24 !定义泊松比 *\:u}'[
MP,DENS,2,2700 !定义密度 :+m|K C(Z
!取强度折减系数F=1.2时边坡围岩材料属性 gb8nST$r
MP,EX,3,3E10 VO[s:e9L
MP,PRXY,3,0.25 rhwY5FD ?
MP,DENS,3,2500 }-YD_Pm
K-
TB,DP,3 !定义数据表3 vs HY; [
TBDATA,1,7.5E5,37.7 !给数据表输入数据 FZb\VUmnV
!取强度折减系数F=1.4时边坡围岩材料属性 l9? ] t;
MP,EX,4,3E10 >\3=h8zw
MP,PRXY,4,0.25 !定义泊松比 WX]kez{<uP
MP,DENS,4,2500 >*{\N^:z
TB,DP,4 !定义数据表4 pH' #v]"
TBDATA,1,6.4E5,33.5 !给数据表输入数据
8X[G)J;
!取强度折减系数F=1.6时边坡围岩材料属性 |jF)~k6
MP,EX,5,3E10 F;5.nKo
MP,PRXY,5,0.25 B k#68p
MP,DENS,5,2500 ！定义密度 X^r HugQ
TB,DP,5 !定义数据表5 l 0jjLqm :
TBDATA,1,5.6E5,30 !给数据表输入数据 IHMZ E42
!取强度折减系数F=1.8时边坡围岩材料属性 Xn9TQ"[4
MP,EX,6,3E10 ！定义弹性模量 ~z%K9YcyU
MP,PRXY,6,0.25 B T gL:
MP,DENS,6,2500 69TQH J[
TB,DP,6 !定义数据表6 uZ\ >
TBDATA,1,5.0E5,27.2 8cF-kfbfZ
!取强度折减系数F=2.0时边坡围岩材料属性 Bp>%'L
MP,EX,7,3E10 T6%*t#8r
MP,PRXY,7,0.25 !定义泊松比 T&X*[kP
MP,DENS,7,2500 4 Z1- RS
TB,DP,7 !定义数据表7 N-4LdC
TBDATA,1,4.5E5,24.8 fG2\p&z
!取强度折减系数F=2.2时边坡围岩材料属性 K zKHC
MP,EX,8,3E10 (]vHW+'
MP,PRXY,8,0.25 F+v? 2|03
MP,DENS,8,2500 # ]&=]K1V
TB,DP,8 !定义数据表8 Xc"S"a^\%
TBDATA,1,4.09E5,22.8 f[!N]*
!取强度折减系数F=2.4时边坡围岩材料属性 \y"!`.E7\d
MP,EX,9,3E10 r,!7TuBl
MP,PRXY,9,0.25 Iy6 "2$%a
MP,DENS,9,2500 G*3O5m
TB,DP,9 !定义数据表9 z?Hvh
TBDATA,1,3.6E5,21.1 !给数据表输入数据 I[W W1P5
!取强度折减系数F=2.6时边坡围岩材料属性 rwv_
RN
MP,EX,10,3E10 ！定义弹性模量 Dwe_ytjpc
MP,PRXY,10,0.25 !定义泊松比 srw5&s(3X
MP,DENS,10,2500 ！定义密度 w}K<,5I>
TB,DP,10 !定义数据表10 `BdZqXKG
TBDATA,1,3.46E5,19.6 !给数据表10输入数据 .~4%TsBa Y
!取强度折减系数F=2.8时边坡围岩材料属性 g\OPidY
MP,EX,11,3E10 ！定义弹性模量 n$F&gx'^
MP,PRXY,11,0.25 !定义泊松比 v|<Dc8i+
MP,DENS,11,2500 =YE"6 iU
TB,DP,11 !定义数据表11 -r_z,h|
TBDATA,1,3.2E5,18.2 !给数据表11输入数据 ,A h QA
!取强度折减系数F=3.0时边坡围岩材料属性 ? D2:'gg
MP,EX,12,3E10 ！定义弹性模量 aM, >LKNbQ
MP,PRXY,12,0.25 b 6'%nR*f
MP,DENS,12,2500 ！定义密度 =Bu d!
TB,DP,12 !定义数据表12 My&h{Qk
TBDATA,1,3.0E5,17.1 !给数据表12输入数据 Z; r}G m
!建立几何模型 ?9i
7 w1`
!创建边坡线模型线 XTXRC$B
K,1,,,, ！ 创建关键点1 xbxU`2/
K,2,-800,,, gmJiKuAL5
K,3,-800,-800,, 6,ylk f3
K,4,-800,-1200,, GX4# I Rq
K,5,1200,-1200,, O|O #T.Tg
K,6,1200,-800,, ZUVk~X3
K,7,1200,,, %($sj| _l
K,8,1200,400,, A=Q"IdK
K,9,400,400,, vbXuT$
l,1,2 !连接关键点1、2生成直线L1 #O7|&DqF{
l,2,3 `[ 0.G0i
l,3,4 l[fU0;A
l,4,5 V)^ nVD)e
l,5,6 w[qWr@
l,6,7 f5IO<(:E^
l,7,8 Sx5r u?$.
l,8,9 XwOj`N{!H
l,9,1 !rlN|HB
l,1,7 d/k&f5
l,3,6 JVD#wwic
!创建边坡面模型 3gQQ,V..
al,3,4,5,11 !由线L3、L4、L5、L11生成一个面A1 =m xmJFA
al,1,2,11,6,10 !由线L1、L2、L1、L6、L10生成一个面A2 ,ic.b
@u1
al,7,8,9,10 !由线L7、L8、L9、L10生成一个面A3 22a$//}E
/PNUM, area, on !打开面号开关 qh Rs5QXL
!划分网格生成有限元模型 -4m UGh1dy
! 划分边坡围岩2网格 {*bXO8vi((
!设置网格份数 QHUoAa`6v
lsel,s,,,3,5,2 !选择线L3和L5 M~t S
*
lesize,all,,,5 !把所选择线分为5段 Ratg!l|'-
lsel,s,,,4,11,7 !选择线L4和L11 N J9H=
lesize,all,,,26 !把所选择线分为26段 #*:y2W%H
mat,2 !给边坡围岩2赋予2号材料特性 g<-x" $(C&
type,1 !采用单元类型1 CL5u{ i5
mshkey,0 !设定自由网格划分 3 (Kj| u
mshape,0 !设定四边形网格划分 DB-4S-2
amesh,1 ！划分面积A1 oC#@9>+@+"
! 划分边坡围岩1网格 2_6x2I a4
!设置网格份数 ] O3[Te
lsel,s,,,2,6,4 !选择线L2和L6 Af pB=3
lsel,a,,,1 !附加选择线L1 +e. bO5Y
lesize,all,,,10 !把所选择线分为10段 -8X* (7
lsel,s,,,8,10,2 !选择线L8和L10 lwIU|T<4
lesize,all,,,16 !把所选择线分为16段 US 9cuah1/
lsel,s,,,7,9,2
S =!3t`
lesize,all,,,12 7Q7z6p/\v
mat,1 !给边坡围岩1赋予1号材料特性 m/r4f279
type,1 i9[=x(-@
mshkey,0 !设定自由网格划分 C S"2Sd 1`
mshape,0 !设定四边形网格划分 y,v*jE
amesh,2 !划分面积A2 'oBT*aL
amesh,3 ！划分面积A3 DJr{;t$7~
allsel a*o#,T5A
!保存网格模型 qP-*
save,Slope-grid.db t5G@M&d4Eo
(3)施加约束和荷载 5K| 1Y#X
!给边坡模型施加约束 W2J"W=:z
!边坡两侧施加X方向约束 j
y R 9a!
nsel,s,loc,x,-800 !选择X=-800线上所有节点 ]e@'9`G-'
nsel,a,loc,x,1200 !选择X=1200线上所有节点 J'`,];su
d,all,ux !对所选择节点约束X方向位移 8q{
%n
allsel w+(bkqz]
!边坡底部施加约束 % zHs h
nsel,s,loc,y,-1200 !选择Y=--1200线上所有节点 ' \Z54$
d,all,uy !对所选择节点约束X、Y方向位移 (?Fz{
d,all,ux 6":=p:PT.
!施加重力加速度 Rb\M63q
acel,,9.8 N'w ;1,c+
(4)求解 m!L&_ Z|j
/solu c&z@HEzV7
!求解设置 "l6 v[yv
antype,static !设定为静力求解 U@ x5cw:
nsubst,100 !设定最大子步数为100 \^or l9
pred,on !打开时间步长预测器 ]f=108|8
nropt,full !设定牛顿-拉普森选项 N|Mzj|i.
nlgeom,on !打开大位移效果 H]K(`)y}4
lnsrch,on !打开线性搜索 @2c Gx/1#
outres,all,all !输出所有项 D6@c&
cnvtol,f,,0.005,2,0.5 !力收敛准则设定 *6=2UJcJ
cnvtol,u,,0.05,2,1 !位移收敛准则设定 :KGUO{_u
!边坡在强度折减系数F=1时求解 #Y`GWT1==
allsel [ GJ_]w^}j
solve !进行求解 :
&>PN,q>
save,F1,db !把F=1时求解结果保存 eSHyA+ F
!边坡在强度折减系数F=1.2时求解 p "EQ 6_f
finish '+LC.l M
/solu RI9 &KS
allsel ,r!_4|\
asel,s,area,,2,3,1 !选择面积A2、A3 l?L wQmq6
MPCHG,3,all ! 把所选择单元材料号改为3 w K0vKd i
allsel MX2 Zm
solve !进行求解 ^cB83%<Z
save,F1.2,db !把F=1.2时求解结果保存 c L}} ^
!边坡在强度折减系数F=1.4时求解 tP8>0\$)
finish o5)lTVQ~~
/solu J@Z IW%5
allsel :1NF#-2\f
asel,s,area,,2,3,1 !选择面积A2、A3 % H Pwu &
MPCHG,4,all !把所选择单元材料号改为4 aE BQx
allsel Xg E\q
solve !进行求解 kBg,U 8|S
save,F1.4,db !把F=1.4时求解结果保存 ]JF>a_2wG
!边坡在强度折减系数F=1.6时求解 / @- !JF#g
finish VgLrufJ
/solu &6V[@gmD
allsel C)66 ^l!x
asel,s,area,,2,3,1 !选择面积A2、A3 -,+zA.{+W
MPCHG,5,all !把所选择单元材料号改为5 },zP,y:cH
allsel l0N~mes
solve 1{{z[ w#
save,F1.6,db !把F=1.6时求解结果保存 P oTJ4z
!边坡在强度折减系数F=1.8时求解 8i"v7}
finish ] /&qv6D*d
/solu 1'ts>6b
allsel FY1
>{Bn
asel,s,area,,2,3,1 !选择面积A2、A3 0<Q['l4Ar
MPCHG,6,all !把所选择单元材料号改为6 Q |,(C0<G
allsel ?Co)7} N
solve uL| Wuq
save,F1.8,db !把F=1.8时求解结果保存 nKpXR uFn\
!边坡在强度折减系数F=2.0时求解 KG7 ~)g
finish SbS*z:
/solu %<c2jvn+k
allsel ~5S[S l
asel,s,area,,2,3,1 !选择面积A2、A3 :98<dQIG
MPCHG,7,all !把所选择单元材料号改为7 H8@1Kt
allsel *l\wl @{
solve ot Tv,T182
save,F2.0,db !把F=2.0时求解结果保存 #pDWwnP[rt
!边坡在强度折减系数F=2.2时求解 g{0a]'ph
finish .=@xTJh
/solu TKJs'%Q7F6
allsel 5qZebD2a
asel,s,area,,2,3,1 !选择面积A2、A3 UqHk2h-
MPCHG,8,all !把所选择单元材料号改为8 eQK}J]S<
allsel !V/Vy/'` *
solve C{Y0}ZrmlF
save,F2.2,db !把F=2.2时求解结果保存 0m8mHJ<&
!边坡在强度折减系数F=2.4时求解 ]=Im 0s
finish Xm#rkF[,
/solu W&MZ5t ,k=
allsel V
,#
|\
asel,s,area,,2,3,1 !选择面积A2、A3 ]k0Pe;<
MPCHG,9,all !把所选择单元材料号改为9 rBY)rUDd4
allsel 2;&K*>g&.
solve Fl8 *dXG&
save,F2.4,db !把F=2.4时求解结果保存 Jkf%k3H3I*
!边坡在强度折减系数F=2.6时求解 F_ C7S
finish I$yFCd Xr
/solu ha Tmfh_|
allsel v/ $~ifY"
asel,s,area,,2,3,1 !选择面积A2、A3 ?xkw~3Yfi
MPCHG,10,all !把所选择单元材料号改为10 ~O|g~H5;
allsel 8kn> ?
solve w~+C.4=7
save,F2.6,db !把F=2.6时求解结果保存 P_7QZ0k/
!边坡在强度折减系数F=2.8时求解 R9 +0ZoS
finish V c 2 (R^
/solu YCMXF#1
allsel G#uD CF,O
asel,s,area,,2,3,1 !选择面积A2、A3 [Q0n-b,Q
MPCHG,11,all !把所选择单元材料号改为11 (%N=7?
allsel VPC7Dh%.
solve w^B F.Nu
save,F2.8,db !把F=2.8时求解结果保存 goa@ e
!边坡在强度折减系数F=3.0时求解 5f#N$ mh
finish c\P,ct
}>
/solu D5Z@6RVt
allsel E }qW'
asel,s,area,,2,3,1 !选择面积A2、A3
ORCG(N
MPCHG,12,all !把所选择单元材料号改为12
EU5^"\
allsel V=
wWY*C
solve [ :Upn)9
save,F3.0,db !把F=3.0时求解结果保存 \}Hi\k+h':
(5）后处理 d+m}Z>iQ1O
/post1 ！进入后处理 e@j&c:p(Y
!边坡在强度折减系数F=1时结果分析 | XGj97#M
Resume,'F1','db' !读入边坡在强度折减系数F=1时 Fy6(N{hql
set,1,last !读入后一个子步 p+.xye U(
pldisp,1 !绘制边坡模型变形图
]S2F9
plnsol,u,x !绘制边坡模型水平方向位移云图 }F
B]LLi
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 4LY$;J ;2
!边坡在强度折减系数F=1.2时结果分析 Ld H23\
Resume,'F1.2','db' !读入边坡在强度折减系数F=1.2时 UR{OrNg*
set,1,last !读入后一个子步 ]cmq
pldisp,1 !绘制边坡模型变形图 V7S [rI<<r
plnsol,u,x !绘制边坡模型水平方向位移云图 GZq~P l
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 z<I@SI^>
!边坡在强度折减系数F=1.4时结果分析 +hZ {/
Resume,'F1.4','db' !读入边坡在强度折减系数F=1.4时 +sY8<y@%
set,1,last !读入后一个子步 fK
4,k:YC
pldisp,1 !绘制边坡模型变形图 uUIjntSF(
plnsol,u,x !绘制边坡模型水平方向位移云图 ._X|Ye9/
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 ow<z @^ 3'
!边坡在强度折减系数F=1.6时结果分析 \Bvy~UeE)>
Resume,'F1.6','db' !读入边坡在强度折减系数F=1.6时 vB?(|
set,1,last !读入后一个子步 >LAhc 7I
pldisp,1 !绘制边坡模型变形图 [Dq@(Q s'
plnsol,u,x !绘制边坡模型水平方向位移云图 E*"-U!?)l2
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 d Eu\ }y|
!边坡在强度折减系数F=1.8时结果分析 px=k&|l
Resume,'F1.8','db' !读入边坡在强度折减系数F=1.8时 -J&
b~t@
set,1,last !读入后一个子步 AqZ( )p*z
pldisp,1 !绘制边坡模型变形图 j6&q6C X
plnsol,u,x !绘制边坡模型水平方向位移云图 ?PSJQ3BC|
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 B] nu \!
!边坡在强度折减系数F=2.0时结果分析 d xa[9>V
Resume,'F2.0','db' !读入边坡在强度折减系数F=2.0时 p=T6Ix'_2e
set,1,last !读入后一个子步 N5F+h94z]
pldisp,1 !绘制边坡模型变形图 $/)0iL{0
plnsol,u,x !绘制边坡模型水平方向位移云图 77/j}Pxh
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 v(GT+i)|
!边坡在强度折减系数F=2.2时结果分析 D:56>%y@
Resume,'F2.2','db' !读入边坡在强度折减系数F=2.2时 Ly~s84k_po
set,1,last !读入后一个子步 Xw'Y
&!z
pldisp,1 !绘制边坡模型变形图 $5r1Si)
plnsol,u,x !绘制边坡模型水平方向位移云图 D,,$
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 C"cBlru8B
!边坡在强度折减系数F=2.4时结果分析 &0Bs ?oq_
Resume,'F2.4','db' !读入边坡在强度折减系数F=2.4时 Aq*,cOF+
set,1,last !读入后一个子步 1 i=p5,|
pldisp,1 !绘制边坡模型变形图 +\ GuZ 5`
plnsol,u,x !绘制边坡模型水平方向位移云图 [+3~wpU(p
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 3d;w\#? L;
!边坡在强度折减系数F=2.6时结果分析 [ApAd
Resume,'F2.6','db' !读入边坡在强度折减系数F=2.6时 yyv<MSU8
set,1,last !读入后一个子步 Yw6d-5=:
pldisp,1 !绘制边坡模型变形图 }n t,DG!r
plnsol,u,x !绘制边坡模型水平方向位移云图 i}_d&.DbF
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 O|e/(s?$
!边坡在强度折减系数F=2.8时结果分析 # n\|Q\W
Resume,'F2.8','db' !读入边坡在强度折减系数F=2.8时 >U#j\2!Sg
set,1,last !读入后一个子步 oFDJwOJ'Bj
pldisp,1 !绘制边坡模型变形图 eFz!`a^dX
plnsol,u,x !绘制边坡模型水平方向位移云图 0N`'a?x
plnsol,eppl,eqv !绘制边坡模型塑性应变云图 =| O ><O|
!边坡在强度折减系数F=3.0时结果分析 ApotRr$)
Resume,'F3.0','db' !读入边坡在强度折减系数F=3.0时 5X}O Un8
set,1,last !读入后一个子步 '(fCi
pldisp,1 !绘制边坡模型变形图 ZK?:w^Z
plnsol,u,x !绘制边坡模型水平方向位移云图 E>L_$J -A-
plnsol,eppl,eqv !绘制边坡模型塑性应变云图
