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

各位朋友： ;j/ur\37  
    大家好！我通过ansys用有限元强度折减法分析了一个边坡的稳定性，结果不同的工况下得出的塑性区是一样的，我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激！附命令流和塑性区如下：(1)创建物理环境 8-SVg
o(  
/COM, Structural                              ! 指定结构分析 wkd591d*  
/TITLE, Slope stability Analysis                  ! 定义工作标题 `I'=d4  
/FILNAM,support,1                            ! 定义工作文件名 ^kXDEKm  
(2)建立模型 B
B|{VwN  
!进入前处理器 @
cjhri|vH  
/PREP7 crhck'?0  
!定义单元类型 r_2  
ET,1,PLANE82                              !定义边坡围岩单元 DF&C7+hO  
KEYOPT,1,3,2                                !求解类型为平面应变 9r8{9h:  
!定义材料属性
j,G/[V  
!边坡围岩1材料属性   |u$AzI  
MP,EX,1,3E10  }"cb^3  
MP,PRXY,1,0.25    WO}JIExy  
MP,DENS,1,2500                            !定义塑性模型 \8t g7Sdq  
TB,DP,1                                    !定义数据表1 <MfB;M  
TBDATA,1,9E5,42.8                          !给数据表输入数据 [e&$4l IS  
!边坡围岩2材料属性                        !采用弹性模型 +Oyt  
MP,EX,2,3.2E10                              ！定义弹性模量 q_ 5xsTlTR  
MP,PRXY,2,0.24                              !定义泊松比 Rw. Uz&  
MP,DENS,2,2700                            !定义密度 ipb
VQ7  
!取强度折减系数F=1.2时边坡围岩材料属性 r/{VL3}F_e  
MP,EX,3,3E10  F_I! +  
MP,PRXY,3,0.25    'Y&yt"cs  
MP,DENS,3,2500 {) .=G  
TB,DP,3                                    !定义数据表3 |r /}r,t}  
TBDATA,1,7.5E5,37.7                        !给数据表输入数据 Gwvs~jN  
!取强度折减系数F=1.4时边坡围岩材料属性 !'[sV^ds  
MP,EX,4,3E10  M_XZOlW5  
MP,PRXY,4,0.25                              !定义泊松比 3cfkJ|fuwe  
MP,DENS,4,2500 Hh qx)u  
TB,DP,4                                    !定义数据表4 uvmNQg   
TBDATA,1,6.4E5,33.5                        !给数据表输入数据 |"@E"Za^  
!取强度折减系数F=1.6时边坡围岩材料属性 ]X-ZRmB`  
MP,EX,5,3E10  x?j&Jn_@w  
MP,PRXY,5,0.25    HV?awc  
MP,DENS,5,2500                            ！定义密度 _5X}&>>lhF  
TB,DP,5                                    !定义数据表5 zD?$O7 |ZK  
TBDATA,1,5.6E5,30                          !给数据表输入数据 -l,ib=ne  
!取强度折减系数F=1.8时边坡围岩材料属性 <T7y85  
MP,EX,6,3E10                              ！定义弹性模量 pCC0:  
MP,PRXY,6,0.25    >?tcL*  
MP,DENS,6,2500 uZQ)A,#n;  
TB,DP,6                                    !定义数据表6 'Lm.`U  
TBDATA,1,5.0E5,27.2 5Y,e}+I>  
!取强度折减系数F=2.0时边坡围岩材料属性 ;N/c5+  
MP,EX,7,3E10  YobIbpo  
MP,PRXY,7,0.25                            !定义泊松比 Ri-I+7(n!  
MP,DENS,7,2500 c?opVbJB\  
TB,DP,7                                    !定义数据表7 @FbzKHd
V/  
TBDATA,1,4.5E5,24.8 Nf;vUYP  
!取强度折减系数F=2.2时边坡围岩材料属性 2cmqtlW"  
MP,EX,8,3E10  l"cO@.T3  
MP,PRXY,8,0.25    `i.f4]r  
MP,DENS,8,2500 -iW[cj R`$  
TB,DP,8                                    !定义数据表8 cJhf{{_oR  
TBDATA,1,4.09E5,22.8 8XY4  
!取强度折减系数F=2.4时边坡围岩材料属性 Aa1 |{^$:L  
MP,EX,9,3E10  Ik}*7D  
MP,PRXY,9,0.25    v =u|D$  
MP,DENS,9,2500
a.Mp1W  
TB,DP,9                                      !定义数据表9 h{h=',o1  
TBDATA,1,3.6E5,21.1                          !给数据表输入数据 Ia(A&Za  
!取强度折减系数F=2.6时边坡围岩材料属性 I&cb5j]C  
MP,EX,10,3E10                                ！定义弹性模量  V*)6!N[5  
MP,PRXY,10,0.25                              !定义泊松比 D&oC1  
MP,DENS,10,2500                          ！定义密度 .<j8>1  
TB,DP,10                                  !定义数据表10 3Ys|M%N  
TBDATA,1,3.46E5,19.6                        !给数据表10输入数据 fO:*85%}7  
!取强度折减系数F=2.8时边坡围岩材料属性 jV7q)\uu^  
MP,EX,11,3E10                              ！定义弹性模量 Sqb#U{E  
MP,PRXY,11,0.25                            !定义泊松比 [@\f 0R  
MP,DENS,11,2500 "@f`O  
TB,DP,11                                  !定义数据表11 HSUr  
TBDATA,1,3.2E5,18.2                        !给数据表11输入数据 n@%'Nbc>b  
!取强度折减系数F=3.0时边坡围岩材料属性 ;=n7 Z  
MP,EX,12,3E10                              ！定义弹性模量 Et- .[  
MP,PRXY,12,0.25    [t6)M~&e:_  
MP,DENS,12,2500                            ！定义密度 v:vA=R2  
TB,DP,12                                  !定义数据表12 ^X_%e|  
TBDATA,1,3.0E5,17.1                        !给数据表12输入数据 )Z^(+  
!建立几何模型 Q \{\uJ x  
!创建边坡线模型线 j"NqNv  
K,1,,,,                                      ！ 创建关键点1 dIf Jr}ih  
K,2,-800,,, SS`\,%aog  
K,3,-800,-800,,  )pvZM?  
K,4,-800,-1200,, leD?yyjw7  
K,5,1200,-1200,,  @lJGdp  
K,6,1200,-800,,  .*zQ\P  
K,7,1200,,, m?m,w$K  
K,8,1200,400,, %r>vZ/>a  
K,9,400,400,, D~%cf  
l,1,2                                        !连接关键点1、2生成直线L1 >kB?C!\  
l,2,3 UGN. ]#"#  
l,3,4 +V'Z%;/  
l,4,5 +g.WO5A  
l,5,6 fe Q%L  
l,6,7 U}5uy9A  
l,7,8 TQEZ<B$
  
l,8,9 Y0b.utR&  
l,9,1 V~+Unn  
l,1,7 -\I0*L'$|\  
l,3,6 +fwq9I>L  
!创建边坡面模型 |4'E&(BU-  
al,3,4,5,11                              !由线L3、L4、L5、L11生成一个面A1 X;0EgIqh3  
al,1,2,11,6,10                            !由线L1、L2、L1、L6、L10生成一个面A2 UtWoSFZ'o!  
al,7,8,9,10                              !由线L7、L8、L9、L10生成一个面A3 nBD7  
/PNUM, area, on                        !打开面号开关 %v|,-B7Yx  
!划分网格生成有限元模型 \(z)]D  
! 划分边坡围岩2网格 *".7O*jjV  
!设置网格份数 J]~3{Mi  
lsel,s,,,3,5,2                              !选择线L3和L5 g=%W"v  
lesize,all,,,5                              !把所选择线分为5段 77 `/YE#M  
lsel,s,,,4,11,7                              !选择线L4和L11 c^=:]^  
lesize,all,,,26                              !把所选择线分为26段 eIEcj<f  
mat,2                                    !给边坡围岩2赋予2号材料特性 lS,Hr3Lz  
type,1                                  !采用单元类型1 %|u"0/  
mshkey,0                                !设定自由网格划分 :N[2*
.c[  
mshape,0                                !设定四边形网格划分 \UE9Ff+{  
amesh,1                                  ！划分面积A1 MdyH/.Te  
! 划分边坡围岩1网格 Zkz:h7GUG-  
!设置网格份数 Y#Z&$&n  
lsel,s,,,2,6,4                              !选择线L2和L6 hyBSS,I  
lsel,a,,,1                                  !附加选择线L1 '}
OrFN  
lesize,all,,,10                              !把所选择线分为10段 @-7K~in?^  
lsel,s,,,8,10,2                              !选择线L8和L10 S\f^y8*<  
lesize,all,,,16                              !把所选择线分为16段 6[CX[=P30  
lsel,s,,,7,9,2 vw!i)JO8M  
lesize,all,,,12 SZ0Zi\W  
mat,1                                    !给边坡围岩1赋予1号材料特性 {<_}[} XY  
type,1 )Q\;N C=4  
mshkey,0                                  !设定自由网格划分 oZ=e/\[K  
mshape,0                                  !设定四边形网格划分 $b8>SSz  
amesh,2                                  !划分面积A2 VJNPs6  
amesh,3                                  ！划分面积A3 G%>[I6G  
allsel 5
6v<!L5%  
!保存网格模型 &?9.Y,  
save,Slope-grid.db ybB}|4d&  
(3)施加约束和荷载 IGp-`%9
 
!给边坡模型施加约束 &zm5s*yNt  
!边坡两侧施加X方向约束 G%6wk=IH  
nsel,s,loc,x,-800                            !选择X=-800线上所有节点 7.=u:PK7kM  
nsel,a,loc,x,1200                            !选择X=1200线上所有节点 R<wb8iir  
d,all,ux                                    !对所选择节点约束X方向位移 l/1uP  
allsel wbpxJtJB  
!边坡底部施加约束 +6>Pp[%  
nsel,s,loc,y,-1200                          !选择Y=--1200线上所有节点 o3`Z@-.G  
d,all,uy                                  !对所选择节点约束X、Y方向位移 k'\RS6M`L  
d,all,ux                                  ?e ~*,6  
!施加重力加速度 BujWql  
acel,,9.8 pu^1s#g8w  
(4)求解 Oq.) 8E.  
/solu jrW7AT)\  
!求解设置 B 51LZP  
antype,static                                  !设定为静力求解 kiECJ@5p  
nsubst,100                                  !设定最大子步数为100 bhgh ]{  
pred,on                                      !打开时间步长预测器 aRV!0?fS  
nropt,full                                    !设定牛顿-拉普森选项 q
aQ
 
nlgeom,on                                    !打开大位移效果 ZuE0'9  
lnsrch,on                                    !打开线性搜索 G9gvOEI/  
outres,all,all                                  !输出所有项 !7w-?1?D  
cnvtol,f,,0.005,2,0.5                            !力收敛准则设定 Sj*H4ZHD<&  
cnvtol,u,,0.05,2,1                              !位移收敛准则设定 Kzmgy14o  
!边坡在强度折减系数F=1时求解 1"e=Zqn$)  
allsel ~7=,)Q  
solve                                        !进行求解 XrvrN^'  
save,F1,db                                    !把F=1时求解结果保存 ?@u &3/&  
!边坡在强度折减系数F=1.2时求解 R@EFG%|`_  
finish ar{e<&Bny  
/solu {!="PnB  
allsel E6G^?k~q  
asel,s,area,,2,3,1                              !选择面积A2、A3 y3xP~]n  
MPCHG,3,all                                ! 把所选择单元材料号改为3 1"k +K~:  
allsel &fTCY-W[  
solve                                        !进行求解 l%)XPb2$J  
save,F1.2,db                                  !把F=1.2时求解结果保存                                                                                  $~2Ao[  
!边坡在强度折减系数F=1.4时求解 DnTM#i:  
finish oNdO@i%.q4  
/solu a#{a{>  
allsel =8tK]lb  
asel,s,area,,2,3,1                              !选择面积A2、A3 #prYZcHv:_  
MPCHG,4,all                                !把所选择单元材料号改为4 xYv;l\20.  
allsel oxeu%wj_  
solve                                        !进行求解 ,:J[|9  
save,F1.4,db                                  !把F=1.4时求解结果保存                            ic!% }S?  
!边坡在强度折减系数F=1.6时求解 yl@Nyu  
finish
KGc.YUoE  
/solu J %A=  
allsel l
(_|CkcZ  
asel,s,area,,2,3,1                              !选择面积A2、A3 Vq-Kl[-|  
MPCHG,5,all                                !把所选择单元材料号改为5 {Y2J:x  
allsel J%r:"Jm[y1  
solve '-N5F  
save,F1.6,db                                  !把F=1.6时求解结果保存                                                                                  H? z~V-8  
!边坡在强度折减系数F=1.8时求解 VO {z)_  
finish yevJA?C4 v  
/solu d(^HO~p  
allsel DMDtry?1:  
asel,s,area,,2,3,1                                !选择面积A2、A3 >\:GFD{z  
MPCHG,6,all                                  !把所选择单元材料号改为6 &v;fK$=2C  
allsel BGj!/E  
solve L?|
}!  
save,F1.8,db                                  !把F=1.8时求解结果保存  (j@c946z""  
!边坡在强度折减系数F=2.0时求解 &0bq3JGW  
finish d6[' [
dG  
/solu Q=DMfJ"  
allsel oidK_mU9q  
asel,s,area,,2,3,1                                !选择面积A2、A3 -u"|{5? '  
MPCHG,7,all                                  !把所选择单元材料号改为7 @uV]7d"z(  
allsel }{@y]DcdM4  
solve ZD|F"v.  
save,F2.0,db                                  !把F=2.0时求解结果保存                            2NjgLXP  
!边坡在强度折减系数F=2.2时求解 2-
Q5l*  
finish zd$?2y8  
/solu _fY9u2Y  
allsel Hq<4G:#  
asel,s,area,,2,3,1                                !选择面积A2、A3 VKYljY0#  
MPCHG,8,all                                  !把所选择单元材料号改为8 M"p%CbcI]  
allsel $")Gd@aR  
solve svxjad@l/  
save,F2.2,db                                  !把F=2.2时求解结果保存                                                                                  C|IHRw`[  
!边坡在强度折减系数F=2.4时求解 [
$d]U.  
finish GKF!GbGR@  
/solu u$5.GmKm  
allsel `~VL&o1>  
asel,s,area,,2,3,1                                !选择面积A2、A3 IOC$jab@  
MPCHG,9,all                                  !把所选择单元材料号改为9 Rm)hgmZ  
allsel VrfEa d  
solve v_U/0 0
  
save,F2.4,db                                  !把F=2.4时求解结果保存  *$R9'Yo}F  
!边坡在强度折减系数F=2.6时求解 Wvu1?  
finish yN`&oya  
/solu 45[,LJaMd  
allsel ~J~R.r/  
asel,s,area,,2,3,1                                !选择面积A2、A3 9yTDuhJ6  
MPCHG,10,all                                  !把所选择单元材料号改为10 r 20!  
allsel ) ] Ro  
solve os;94yd)  
save,F2.6,db                                  !把F=2.6时求解结果保存                            r|uR!=*|?  
!边坡在强度折减系数F=2.8时求解 wX6-WQR  
finish [Rq|;p  
/solu II _CT=  
allsel ?-3G5yy  
asel,s,area,,2,3,1                                !选择面积A2、A3 ,/V~T<FI  
MPCHG,11,all                                !把所选择单元材料号改为11 a[j]f
v*6  
allsel `#!>}/m  
solve q\jq9)  
save,F2.8,db                                  !把F=2.8时求解结果保存                                                                                  1~E4]Ef:W  
!边坡在强度折减系数F=3.0时求解 bk=;=K  
finish dD39?K/  
/solu 8q6Le{G  
allsel ^jS1g*nrN  
asel,s,area,,2,3,1                            !选择面积A2、A3 t(PA+~sIp  
MPCHG,12,all                              !把所选择单元材料号改为12 CXz9bhn<4  
allsel Tyaqa0  
solve (gvnIoDl0  
save,F3.0,db                                !把F=3.0时求解结果保存 <Pe'&u  
(5）后处理 k^;/@:  
/post1                                      ！进入后处理 Y^d#8^cP  
!边坡在强度折减系数F=1时结果分析 0i8\Lu6  
Resume,'F1','db'                              !读入边坡在强度折减系数F=1时 AEWrrE
 
set,1,last                                    !读入后一个子步 Z]2z*XD  
pldisp,1                                      !绘制边坡模型变形图 (Lh!7g/0N  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 q2`mu4B  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 W;!V_-:  
!边坡在强度折减系数F=1.2时结果分析 :iE`=( o  
Resume,'F1.2','db'                            !读入边坡在强度折减系数F=1.2时 @~7au9.V=X  
set,1,last                                    !读入后一个子步 :wRfk*Ly  
pldisp,1                                      !绘制边坡模型变形图 I(&N2L$-  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 /GJL&RMx  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 r+#g  
!边坡在强度折减系数F=1.4时结果分析 7tJ#0to  
Resume,'F1.4','db'                            !读入边坡在强度折减系数F=1.4时 ~^'t70 :D  
set,1,last                                    !读入后一个子步 g%^/^<ei  
pldisp,1                                      !绘制边坡模型变形图 8 ks\-38n1  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 zZ51jA9x  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 7iv g3*  
!边坡在强度折减系数F=1.6时结果分析 "9hD4R  
Resume,'F1.6','db'                            !读入边坡在强度折减系数F=1.6时 Mkt_pr  
set,1,last                                    !读入后一个子步 s{X+0_@Q  
pldisp,1                                      !绘制边坡模型变形图 !P=L0A`  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 'ju_l)(R  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 ei1
;@k/  
!边坡在强度折减系数F=1.8时结果分析 fo>_*6i74  
Resume,'F1.8','db'                            !读入边坡在强度折减系数F=1.8时 =1!.g"0  
set,1,last                                    !读入后一个子步 wM;=^br  
pldisp,1                                      !绘制边坡模型变形图 g.;2N9  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 C~.\2D`zy  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 9X#]Lg?b  
!边坡在强度折减系数F=2.0时结果分析 <pG 4g  
Resume,'F2.0','db'                            !读入边坡在强度折减系数F=2.0时 bYhG`1,$-a  
set,1,last                                    !读入后一个子步 j3+ hsA/(k  
pldisp,1                                      !绘制边坡模型变形图 "5N$u(: b  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 G'wW-|  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 2t?>0)*m  
!边坡在强度折减系数F=2.2时结果分析 9]G~i`QQ  
Resume,'F2.2','db'                            !读入边坡在强度折减系数F=2.2时 -<T>paE9  
set,1,last                                    !读入后一个子步 +m~3InWq  
pldisp,1                                      !绘制边坡模型变形图 4v5qK  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 c;yp}k]\  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 /9D mK%d  
!边坡在强度折减系数F=2.4时结果分析 `8\"3S  
Resume,'F2.4','db'                            !读入边坡在强度折减系数F=2.4时 ?TW?2+  
set,1,last                                    !读入后一个子步 "_BWUY  
pldisp,1                                      !绘制边坡模型变形图 |VyN>
&r~6  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 v-!Spf  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 &i`(y>\  
!边坡在强度折减系数F=2.6时结果分析 # M3d=  
Resume,'F2.6','db'                            !读入边坡在强度折减系数F=2.6时 ^ rO}'~(  
set,1,last                                    !读入后一个子步 [;D1O;c'W.  
pldisp,1                                      !绘制边坡模型变形图 H#nJWe_9A  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 T|(w-)mv  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 .]<iRf[\[  
!边坡在强度折减系数F=2.8时结果分析 c5KciTD^  
Resume,'F2.8','db'                            !读入边坡在强度折减系数F=2.8时 _]3#C[1L  
set,1,last                                    !读入后一个子步 4/S3hH  
pldisp,1                                      !绘制边坡模型变形图 PuBE=9,  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 I
>-1kFma;  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 ,CqGO %DY  
!边坡在强度折减系数F=3.0时结果分析 CA$|3m9)NM  
Resume,'F3.0','db'                            !读入边坡在强度折减系数F=3.0时 g~=- ,j|  
set,1,last                                    !读入后一个子步 !{g>g%2!  
pldisp,1                                      !绘制边坡模型变形图 ~@}n}aV'!  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 dd6l+z  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图

如果你这个能算出正确的结果的话,才真的叫出错,......... CYhSCT!-?  
 KzZRFEA_  
_6\"U5*Y  
FMqes5\ 3  
mfgUf  
zm=|#f  
#8B4*gAM  
Rw`64L_
  
     

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

不知道二楼的说什么

高人现身啊

做的不错嘛

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

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

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

    

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

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