[讨论]有限元强度折减法 [复制链接]

各位朋友： Ul_M3"Z  
    大家好！我通过ansys用有限元强度折减法分析了一个边坡的稳定性，结果不同的工况下得出的塑性区是一样的，我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激！附命令流和塑性区如下：(1)创建物理环境 Qd_Y\PzS  
/COM, Structural                              ! 指定结构分析 R g?1-|Tj  
/TITLE, Slope stability Analysis                  ! 定义工作标题 Ja]?&j  
/FILNAM,support,1                            ! 定义工作文件名 $B#6tk~u  
(2)建立模型 +B'8|5tPX  
!进入前处理器 .fi/I  
/PREP7 L>!8YUz7p$  
!定义单元类型 rrnNn'  
ET,1,PLANE82                              !定义边坡围岩单元 Zes+/.sA}]  
KEYOPT,1,3,2                                !求解类型为平面应变 yJ
sH=5A  
!定义材料属性 T/c<23i  
!边坡围岩1材料属性  p1D[YeF4  
MP,EX,1,3E10  @v}B6j b;  
MP,PRXY,1,0.25    y8s!M  
MP,DENS,1,2500                            !定义塑性模型 IcrL   
TB,DP,1                                    !定义数据表1 4D$;KokZ  
TBDATA,1,9E5,42.8                          !给数据表输入数据 Uz608u  
!边坡围岩2材料属性                        !采用弹性模型 X53mzs  
MP,EX,2,3.2E10                              ！定义弹性模量 H{?9CxY
a  
MP,PRXY,2,0.24                              !定义泊松比 ?f*Q>3
S)  
MP,DENS,2,2700                            !定义密度 T# lP!c  
!取强度折减系数F=1.2时边坡围岩材料属性 ~1*A  
MP,EX,3,3E10 
yoRU_%xA  
MP,PRXY,3,0.25    hH
->%*  
MP,DENS,3,2500 Um)>2|rp}  
TB,DP,3                                    !定义数据表3 (f_g7B2&y  
TBDATA,1,7.5E5,37.7                        !给数据表输入数据 \dlph  
!取强度折减系数F=1.4时边坡围岩材料属性 4Hb $0l  
MP,EX,4,3E10  '~ 4pl0TWc  
MP,PRXY,4,0.25                              !定义泊松比 vBvNu<v7te  
MP,DENS,4,2500 [EY`am8[  
TB,DP,4                                    !定义数据表4 "o&HE@t  
TBDATA,1,6.4E5,33.5                        !给数据表输入数据 wNNB;n`l  
!取强度折减系数F=1.6时边坡围岩材料属性 yMc:n"-[  
MP,EX,5,3E10  6*n
Ao8gl  
MP,PRXY,5,0.25    3Wcy)y>2Ap  
MP,DENS,5,2500                            ！定义密度 PF`rWw  
TB,DP,5                                    !定义数据表5 'r} zY-FM`  
TBDATA,1,5.6E5,30                          !给数据表输入数据 D6~+Y~R  
!取强度折减系数F=1.8时边坡围岩材料属性 1/ZR*fa  
MP,EX,6,3E10                              ！定义弹性模量 ]DvO:tM  
MP,PRXY,6,0.25    b@t5`Y-+K  
MP,DENS,6,2500 4KI
[D{  
TB,DP,6                                    !定义数据表6 ' )-M\'S$E  
TBDATA,1,5.0E5,27.2 C}?0`!Cc%  
!取强度折减系数F=2.0时边坡围岩材料属性 \OC6M` /  
MP,EX,7,3E10  tfkr+ /  
MP,PRXY,7,0.25                            !定义泊松比 >p])it[q&$  
MP,DENS,7,2500 B|%tE{F  
TB,DP,7                                    !定义数据表7 !r+IXuqV,!  
TBDATA,1,4.5E5,24.8 t` 8!AhOgc  
!取强度折减系数F=2.2时边坡围岩材料属性 Aaw(Ed  
MP,EX,8,3E10  _M,lQ~  
MP,PRXY,8,0.25    4'Xgk
8)  
MP,DENS,8,2500 ;pJ2V2 g8  
TB,DP,8                                    !定义数据表8 !umEyd@ "  
TBDATA,1,4.09E5,22.8 64j|}wJ$  
!取强度折减系数F=2.4时边坡围岩材料属性 ~%#mK:+  
MP,EX,9,3E10  i3mAfDF  
MP,PRXY,9,0.25    K?u(1  
MP,DENS,9,2500 P20|RvE  
TB,DP,9                                      !定义数据表9 *zVvQ=  
TBDATA,1,3.6E5,21.1                          !给数据表输入数据 la$%H<,7  
!取强度折减系数F=2.6时边坡围岩材料属性 [U\
(G  
MP,EX,10,3E10                                ！定义弹性模量  O'NW Ebl/  
MP,PRXY,10,0.25                              !定义泊松比  >0Ev#cX4  
MP,DENS,10,2500                          ！定义密度 'GJVWpvUU  
TB,DP,10                                  !定义数据表10 Vv*NFJ|  
TBDATA,1,3.46E5,19.6                        !给数据表10输入数据 XD-^w_  
!取强度折减系数F=2.8时边坡围岩材料属性 /.V0ag'G  
MP,EX,11,3E10                              ！定义弹性模量 Uh|>Skic4  
MP,PRXY,11,0.25                            !定义泊松比 2GWDEgI1o  
MP,DENS,11,2500 5V-jMB  
TB,DP,11                                  !定义数据表11 kII7z;<^`  
TBDATA,1,3.2E5,18.2                        !给数据表11输入数据 kGN||h  
!取强度折减系数F=3.0时边坡围岩材料属性 +`bC%\T8?  
MP,EX,12,3E10                              ！定义弹性模量 X1A<$Am1  
MP,PRXY,12,0.25    NvtM3  
MP,DENS,12,2500                            ！定义密度 jN/C'\QL  
TB,DP,12                                  !定义数据表12 {UH9i'y:t  
TBDATA,1,3.0E5,17.1                        !给数据表12输入数据 Di=9m
HC  
!建立几何模型 uCzii o`S  
!创建边坡线模型线 9 1P4:6  
K,1,,,,                                      ！ 创建关键点1 t7F.[uWD  
K,2,-800,,, ^m7~:=K7WG  
K,3,-800,-800,,  ivrXwZ7jT  
K,4,-800,-1200,, 6ri#Lw  
K,5,1200,-1200,,  :#="%  
K,6,1200,-800,,  FA>1x*;c  
K,7,1200,,, tSv0" L  
K,8,1200,400,, A8?[6^%O|  
K,9,400,400,, ~U}Mv{y  
l,1,2                                        !连接关键点1、2生成直线L1 q42FPq  
l,2,3 _MYx%Z  
l,3,4 aW{L7N%  
l,4,5 3@5p"X  
l,5,6 ~!+h?[miV  
l,6,7 ARL  
l,7,8 :kw0y  
l,8,9 EychR/s  
l,9,1 [@{0o+.]'H  
l,1,7 <9@7,2  
l,3,6 IftxSaP  
!创建边坡面模型 9 Aivf+  
al,3,4,5,11                              !由线L3、L4、L5、L11生成一个面A1 ;G
m>O7"|@  
al,1,2,11,6,10                            !由线L1、L2、L1、L6、L10生成一个面A2 A['uD<4b  
al,7,8,9,10                              !由线L7、L8、L9、L10生成一个面A3 R
Td^
ImV  
/PNUM, area, on                        !打开面号开关 l+ T,2sd  
!划分网格生成有限元模型 9Z!lmfnJ  
! 划分边坡围岩2网格  pbB2wt  
!设置网格份数 )teFS%  
lsel,s,,,3,5,2                              !选择线L3和L5 17py).\  
lesize,all,,,5                              !把所选择线分为5段 rS~qi}4X  
lsel,s,,,4,11,7                              !选择线L4和L11 BpT"~4oV5  
lesize,all,,,26                              !把所选择线分为26段 'M%5v'$y  
mat,2                                    !给边坡围岩2赋予2号材料特性 n@B{vyy  
type,1                                  !采用单元类型1 w@N{@tG  
mshkey,0                                !设定自由网格划分 zS% m_,t  
mshape,0                                !设定四边形网格划分 /)1v9<vM"  
amesh,1                                  ！划分面积A1 6S*zzJ.0K  
! 划分边坡围岩1网格 P"4Mm
, C  
!设置网格份数 kEK[\f VE  
lsel,s,,,2,6,4                              !选择线L2和L6 lL2-.!]R  
lsel,a,,,1                                  !附加选择线L1 k\(4sY M  
lesize,all,,,10                              !把所选择线分为10段 %yw=[]Vjze  
lsel,s,,,8,10,2                              !选择线L8和L10 N+>'J23d!  
lesize,all,,,16                              !把所选择线分为16段 #)i&DJ^Y  
lsel,s,,,7,9,2 'yT`ef  
lesize,all,,,12 p~>_T7ze  
mat,1                                    !给边坡围岩1赋予1号材料特性 Xu<FDjr  
type,1 d)*(KhYie@  
mshkey,0                                  !设定自由网格划分 sC6r.@[u8t  
mshape,0                                  !设定四边形网格划分 P06.1  
amesh,2                                  !划分面积A2 Pe,;MP\2  
amesh,3                                  ！划分面积A3 T8q[7Zn  
allsel ?Y`zg`  
!保存网格模型 ggCr-  
save,Slope-grid.db u[`v&e  
(3)施加约束和荷载 YeIe\3x!N  
!给边坡模型施加约束 'jbMTI  
!边坡两侧施加X方向约束 G^"H*a  
nsel,s,loc,x,-800                            !选择X=-800线上所有节点 n2na9dX)w  
nsel,a,loc,x,1200                            !选择X=1200线上所有节点 3Wj,}  
d,all,ux                                    !对所选择节点约束X方向位移 r } 7:#XQ  
allsel UT-ewXh  
!边坡底部施加约束 ;qVG \wQq  
nsel,s,loc,y,-1200                          !选择Y=--1200线上所有节点 ZKsQ2"8{M  
d,all,uy                                  !对所选择节点约束X、Y方向位移 8on[%Vk  
d,all,ux                                  /Hk07:"c  
!施加重力加速度 -RCv7U`  
acel,,9.8 m3%
ef  
(4)求解 9L
=;KtE1  
/solu (1y='L2rj  
!求解设置 zg0%>iqO  
antype,static                                  !设定为静力求解 qj;l,Kua  
nsubst,100                                  !设定最大子步数为100 8DL hk  
pred,on                                      !打开时间步长预测器 +QupM  
nropt,full                                    !设定牛顿-拉普森选项 0Q*-g}wXfS  
nlgeom,on                                    !打开大位移效果 p[;8  
lnsrch,on                                    !打开线性搜索 [#zE. TW  
outres,all,all                                  !输出所有项 b60[({A\s&  
cnvtol,f,,0.005,2,0.5                            !力收敛准则设定 yr{5Rp05=  
cnvtol,u,,0.05,2,1                              !位移收敛准则设定 D3]@i&^B  
!边坡在强度折减系数F=1时求解 {*nEKPq(_*  
allsel iiehrK&T!  
solve                                        !进行求解 =u[rOU{X"W  
save,F1,db                                    !把F=1时求解结果保存 u mqKFM$  
!边坡在强度折减系数F=1.2时求解 <!hpfTz*  
finish g?$e^ls  
/solu ]g7HEB.Y  
allsel $*aE$O6l  
asel,s,area,,2,3,1                              !选择面积A2、A3 >EeAPO4  
MPCHG,3,all                                ! 把所选择单元材料号改为3 /9TL&_A-T  
allsel B^Bbso'{1  
solve                                        !进行求解 G#g{3}dcK  
save,F1.2,db                                  !把F=1.2时求解结果保存                                                                                  ~\-r  
!边坡在强度折减系数F=1.4时求解 yj]ML:n  
finish ylT6h_z1[Y  
/solu $KH@,;Xz  
allsel sk\_[p  
asel,s,area,,2,3,1                              !选择面积A2、A3 W
kcH5[  
MPCHG,4,all                                !把所选择单元材料号改为4 l92!2$]b  
allsel +Gp!cGaAm  
solve                                        !进行求解 &:C{
/QnA  
save,F1.4,db                                  !把F=1.4时求解结果保存                            <Q\H
  
!边坡在强度折减系数F=1.6时求解 kYmo7  
finish JgxA^>|9;  
/solu s|BX>1  
allsel 8_+vb#M  
asel,s,area,,2,3,1                              !选择面积A2、A3 oFwG+W/  
MPCHG,5,all                                !把所选择单元材料号改为5 ,?t}NZY&  
allsel `6n!$Cxo  
solve Y+OYoI  
save,F1.6,db                                  !把F=1.6时求解结果保存                                                                                  Kyyih|{  
!边坡在强度折减系数F=1.8时求解 '[n)N@h  
finish u% r!?-z  
/solu N!.kq4$.  
allsel giQ{Xrj  
asel,s,area,,2,3,1                                !选择面积A2、A3 @OBHAoz%/  
MPCHG,6,all                                  !把所选择单元材料号改为6 #%:`p9p.S  
allsel //6^+-he  
solve e2qpJ4i  
save,F1.8,db                                  !把F=1.8时求解结果保存  ./';P<)  
!边坡在强度折减系数F=2.0时求解 IsWcz+1n  
finish BQ5_s,VM  
/solu u~ VswXc4  
allsel [}jj<!9A_;  
asel,s,area,,2,3,1                                !选择面积A2、A3 ?'TA!MR  
MPCHG,7,all                                  !把所选择单元材料号改为7 kT>r<`rt  
allsel ;`;G/1]#9  
solve 5#yJK>a7  
save,F2.0,db                                  !把F=2.0时求解结果保存                            A8#.1uEgNb  
!边坡在强度折减系数F=2.2时求解 RCoeJ|  
finish Kp$_0
 
/solu B#| Z`mZ  
allsel T-hU+(+hg  
asel,s,area,,2,3,1                                !选择面积A2、A3 NW
}>pb9  
MPCHG,8,all                                  !把所选择单元材料号改为8 xYT#!K1*  
allsel AlAh S<  
solve d3\OHkM0^  
save,F2.2,db                                  !把F=2.2时求解结果保存                                                                                  DN4fP-m-  
!边坡在强度折减系数F=2.4时求解 _ .-o%6  
finish TEH*@~P"  
/solu MT;SRAmUr  
allsel J(3gT}z-  
asel,s,area,,2,3,1                                !选择面积A2、A3 3D]2$a_d  
MPCHG,9,all                                  !把所选择单元材料号改为9 }C_G0'"F  
allsel 8=lHUn9l  
solve E`A<]dAoK  
save,F2.4,db                                  !把F=2.4时求解结果保存  ./7&_9|<  
!边坡在强度折减系数F=2.6时求解 |LXrGyk^  
finish o+]Y=r2  
/solu 7KvXTrN!9  
allsel fF8a 1XV  
asel,s,area,,2,3,1                                !选择面积A2、A3 ?ZSXoy-kr  
MPCHG,10,all                                  !把所选择单元材料号改为10 SJ1w1^#Pz  
allsel e6tH/`Uln  
solve TBoM{s=.  
save,F2.6,db                                  !把F=2.6时求解结果保存                            n25irCD`  
!边坡在强度折减系数F=2.8时求解 U:s}/to  
finish Z/XM`Cy  
/solu xn%l  
allsel

q[PD  
asel,s,area,,2,3,1                                !选择面积A2、A3 Sb@{f<3E  
MPCHG,11,all                                !把所选择单元材料号改为11 NqQM!B]  
allsel N
#<zEAB  
solve ZXF
AuF  
save,F2.8,db                                  !把F=2.8时求解结果保存                                                                                  5`1(}  
!边坡在强度折减系数F=3.0时求解 t(Q&H!~e   
finish

B[) [fE  
/solu m;;0 Cl  
allsel $|`t9-EA/  
asel,s,area,,2,3,1                            !选择面积A2、A3 .g6PrhzFbk  
MPCHG,12,all                              !把所选择单元材料号改为12 hLDch5J5~  
allsel H<r
nJ  
solve $V`KrA~]  
save,F3.0,db                                !把F=3.0时求解结果保存 .}L-c>o"o  
(5）后处理 >!HfH(is\  
/post1                                      ！进入后处理 iBGS
BSeL&  
!边坡在强度折减系数F=1时结果分析 1q}u?7nnSG  
Resume,'F1','db'                              !读入边坡在强度折减系数F=1时 4y
J*85e]  
set,1,last                                    !读入后一个子步 [H6X2yjj|  
pldisp,1                                      !绘制边坡模型变形图 2uJNc!&  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 (>!]A6^L~  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 l|R
<F;|  
!边坡在强度折减系数F=1.2时结果分析 yV 9]_k  
Resume,'F1.2','db'                            !读入边坡在强度折减系数F=1.2时 yr/]xc$  
set,1,last                                    !读入后一个子步 )vE
HLp.  
pldisp,1                                      !绘制边坡模型变形图 rL!_&|  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 ain#_H  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 VuX>  
!边坡在强度折减系数F=1.4时结果分析 !I~C\$^U  
Resume,'F1.4','db'                            !读入边坡在强度折减系数F=1.4时 1sIy*z  
set,1,last                                    !读入后一个子步 :{TmR3.  
pldisp,1                                      !绘制边坡模型变形图 /ZIJ<#o[
 
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 p]1yd;Jt  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 ;Sx'O  
!边坡在强度折减系数F=1.6时结果分析 avHD'zU}N  
Resume,'F1.6','db'                            !读入边坡在强度折减系数F=1.6时 T!^?d5uW#  
set,1,last                                    !读入后一个子步 'XZI{q2i  
pldisp,1                                      !绘制边坡模型变形图 Tdk2436=  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 KG4#BY&^  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 |1<]o;:  
!边坡在强度折减系数F=1.8时结果分析 *MnG-\{j  
Resume,'F1.8','db'                            !读入边坡在强度折减系数F=1.8时 ]2aYi9)  
set,1,last                                    !读入后一个子步 Oet#wp/I  
pldisp,1                                      !绘制边坡模型变形图 yKe*<\  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 cKi^C  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 *x*,I,03  
!边坡在强度折减系数F=2.0时结果分析 r{/
G\  
Resume,'F2.0','db'                            !读入边坡在强度折减系数F=2.0时 9rpg10/T  
set,1,last                                    !读入后一个子步 EL 8N[]RF  
pldisp,1                                      !绘制边坡模型变形图 z'\}/k+  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 6`s%%v  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 .%J?T5D  
!边坡在强度折减系数F=2.2时结果分析 oqd N5+xt  
Resume,'F2.2','db'                            !读入边坡在强度折减系数F=2.2时 AihL>a%  
set,1,last                                    !读入后一个子步 YvxMA#  
pldisp,1                                      !绘制边坡模型变形图 1{hoO<CJ  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 <.AC=4@V  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 PcN
fTB{  
!边坡在强度折减系数F=2.4时结果分析 f'zU^/$rf  
Resume,'F2.4','db'                            !读入边坡在强度折减系数F=2.4时 fzkCI  
set,1,last                                    !读入后一个子步 @sgT[P*ut  
pldisp,1                                      !绘制边坡模型变形图 XRz%KVysp  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 D_ Bx>G9  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 zT4ulXN  
!边坡在强度折减系数F=2.6时结果分析 /BjGAa(  
Resume,'F2.6','db'                            !读入边坡在强度折减系数F=2.6时 JK:mQ_  
set,1,last                                    !读入后一个子步 qUoMg%Z%l  
pldisp,1                                      !绘制边坡模型变形图 k8x&aH  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 xT=kxyu  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 @F3-Ugm  
!边坡在强度折减系数F=2.8时结果分析 %!HnGwv-  
Resume,'F2.8','db'                            !读入边坡在强度折减系数F=2.8时 }n2-*{)x  
set,1,last                                    !读入后一个子步 aG8D%i0  
pldisp,1                                      !绘制边坡模型变形图 Y#SmZ*zok  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 fa+W9  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 $!l2=^\3  
!边坡在强度折减系数F=3.0时结果分析 <igx[2X  
Resume,'F3.0','db'                            !读入边坡在强度折减系数F=3.0时 _;J9q}X  
set,1,last                                    !读入后一个子步 5|~r{w)9  
pldisp,1                                      !绘制边坡模型变形图 Hh/ -^G  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 XoxR5arj  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图

如果你这个能算出正确的结果的话,才真的叫出错,......... etK,zEd  
/L|}Y242  
>3 .ep},  
hL{B9?  
g/X=#!  
]qza*ba  
{
5tEsv  
vS t=Ax3]  
     

好强啊，楼主能帮忙告诉我一下怎样才能输出网格上的节点的坐标呢？谢谢楼主，我写论文急用

不知道二楼的说什么

高人现身啊

做的不错嘛

顶！！！！！！！！！！！

晕  错误很明显呀 比如：Resume,'F3.0','db'                            !读入边坡在强度折减系数F=3.0时 L/LNX{|  
set,1,last    你每次读入的都是第一荷载步，结果肯定是一样的  EoCwS  
有没有人具体用ANSYS算过郑院士论文上的实例？结果相同吗？欢迎加入边破加固群：64321843

你这个命令流好像是错的，在李围那本书里有个例子是ANSYS强度折减法的，界面操作，推荐你看看

    

谢谢前辈，刚入行，多指教

是对整体都折减吗，还是只对滑体折减啊???????