论坛首页    职业区    学术与写作    工程技术区    软件区    资料区    商务合作区    社区办公室
 您好! 欢迎 登录注册 最新帖子 邀请注册 活动聚焦 统计排行 社区服务 帮助
 
  • 帖子
  • 日志
  • 用户
  • 版块
  • 群组
帖子
  • 2631阅读
  • 1回复

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

上一主题 下一主题
离线fhzhw
 

发帖
426
土币
1699
威望
1245
原创币
0
只看楼主 倒序阅读 使用道具 楼主  发表于: 2008-05-12
各位朋友: 2B*9]AHny  
    大家好!我通过ansys用有限元强度折减法分析了一个边坡的稳定性,结果不同的工况下得出的塑性区是一样的,我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激!附命令流和塑性区如下:(1)创建物理环境 e0aeiG$/0  
/COM, Structural                              ! 指定结构分析 _p90Zm-3X  
/TITLE, Slope stability Analysis                  ! 定义工作标题 {Ynr(J.  
/FILNAM,support,1                            ! 定义工作文件名 Qn6&M  
(2)建立模型 ;[*7UE+#7  
!进入前处理器 bhFzu[B  
/PREP7 0chBw~@*s  
!定义单元类型 VpD9!;S  
ET,1,PLANE82                              !定义边坡围岩单元  9((v.  
KEYOPT,1,3,2                                !求解类型为平面应变 ./[t'dgC  
!定义材料属性 F+}MW/ra@  
!边坡围岩1材料属性  \.}* s]6  
MP,EX,1,3E10  VdYu| w ;v  
MP,PRXY,1,0.25    9GCxF`OB  
MP,DENS,1,2500                            !定义塑性模型 TMq\}k-I5  
TB,DP,1                                    !定义数据表1 i)eub`uMy  
TBDATA,1,9E5,42.8                          !给数据表输入数据 =0O`VSb  
!边坡围岩2材料属性                        !采用弹性模型 ?z?IEj}  
MP,EX,2,3.2E10                              !定义弹性模量 Pxl7zz&pl=  
MP,PRXY,2,0.24                              !定义泊松比 `L0}^ |`9  
MP,DENS,2,2700                            !定义密度 ]"F0"UH,  
!取强度折减系数F=1.2时边坡围岩材料属性 {E>kFeg  
MP,EX,3,3E10  I=. 98v%  
MP,PRXY,3,0.25    W4X=.vr  
MP,DENS,3,2500 w+[r$+z!k  
TB,DP,3                                    !定义数据表3 RW%e%  
TBDATA,1,7.5E5,37.7                        !给数据表输入数据 /}(w{6C  
!取强度折减系数F=1.4时边坡围岩材料属性 #knpZ'  
MP,EX,4,3E10  r"k\G\,%  
MP,PRXY,4,0.25                              !定义泊松比 St'3e<  
MP,DENS,4,2500 z0doL b^!  
TB,DP,4                                    !定义数据表4 yn`H}@`k  
TBDATA,1,6.4E5,33.5                        !给数据表输入数据 XlR.Y~  
!取强度折减系数F=1.6时边坡围岩材料属性 }Vk#w%EJ  
MP,EX,5,3E10  ~-|K5  
MP,PRXY,5,0.25    bT0CQ_g21  
MP,DENS,5,2500                            !定义密度 _0ep[r  
TB,DP,5                                    !定义数据表5 >P/.X^G0  
TBDATA,1,5.6E5,30                          !给数据表输入数据 qm] k (/w  
!取强度折减系数F=1.8时边坡围岩材料属性 ^ WidA-  
MP,EX,6,3E10                              !定义弹性模量 @G^ l`%  
MP,PRXY,6,0.25    yX<Sk q  
MP,DENS,6,2500 n:'Mpux  
TB,DP,6                                    !定义数据表6 c]#F^(-A`  
TBDATA,1,5.0E5,27.2 uq{w1O5  
!取强度折减系数F=2.0时边坡围岩材料属性 e")s1`  
MP,EX,7,3E10  %Q]m6ciAM  
MP,PRXY,7,0.25                            !定义泊松比 \za 0?b  
MP,DENS,7,2500 A6pPx1-&  
TB,DP,7                                    !定义数据表7 t4~Bn<=  
TBDATA,1,4.5E5,24.8 evz{@;.R  
!取强度折减系数F=2.2时边坡围岩材料属性 6|~N5E~SX  
MP,EX,8,3E10  FUiEayM  
MP,PRXY,8,0.25    4f'1g1@$  
MP,DENS,8,2500 Xw_AZ-|1D  
TB,DP,8                                    !定义数据表8 e@w-4G(;  
TBDATA,1,4.09E5,22.8 5Ta<$t  
!取强度折减系数F=2.4时边坡围岩材料属性 "NMX>a,(  
MP,EX,9,3E10  Q6]SsV?x  
MP,PRXY,9,0.25    {OH "d  
MP,DENS,9,2500 /X9Kg  
TB,DP,9                                      !定义数据表9 =0 mf  
TBDATA,1,3.6E5,21.1                          !给数据表输入数据 !^:b?M  
!取强度折减系数F=2.6时边坡围岩材料属性 #2^eGhwnI  
MP,EX,10,3E10                                !定义弹性模量  sWa`-gc  
MP,PRXY,10,0.25                              !定义泊松比 criOJ-  
MP,DENS,10,2500                          !定义密度 zU}Ru&T9  
TB,DP,10                                  !定义数据表10 h3[^uY e  
TBDATA,1,3.46E5,19.6                        !给数据表10输入数据 *@^9 ]$*$  
!取强度折减系数F=2.8时边坡围岩材料属性 ER;?[!  
MP,EX,11,3E10                              !定义弹性模量 NI,i)OSEN  
MP,PRXY,11,0.25                            !定义泊松比 Gx,<|v  
MP,DENS,11,2500 MZZEqsD5[  
TB,DP,11                                  !定义数据表11 W+k SL{0  
TBDATA,1,3.2E5,18.2                        !给数据表11输入数据 3"!h+dXw  
!取强度折减系数F=3.0时边坡围岩材料属性 h/?$~OD  
MP,EX,12,3E10                              !定义弹性模量 4,4S5u[|  
MP,PRXY,12,0.25    Na=.LW-ma=  
MP,DENS,12,2500                            !定义密度 .hPk}B/KV  
TB,DP,12                                  !定义数据表12 14Y_ oH9  
TBDATA,1,3.0E5,17.1                        !给数据表12输入数据 leR-oeSO  
!建立几何模型 1Tm,#o  
!创建边坡线模型线 $F2 A  
K,1,,,,                                      ! 创建关键点1 waT'|9{  
K,2,-800,,, `n)e] dn  
K,3,-800,-800,,  | KY6IGcqV  
K,4,-800,-1200,, -Ri/I4Xj  
K,5,1200,-1200,,  ObnQ,x(  
K,6,1200,-800,,  zZR_&z<  
K,7,1200,,, O*jTrZ(k  
K,8,1200,400,, X\sOeb:]  
K,9,400,400,, !Pd@0n4  
l,1,2                                        !连接关键点1、2生成直线L1 5e1;m6  
l,2,3 vLCyT=OB`  
l,3,4 ||yx?q6\h  
l,4,5 @$mh0K>  
l,5,6 &|MdBJ  
l,6,7 GO4IAUA  
l,7,8 &UrPb%=2H  
l,8,9 ?gYQE&M !  
l,9,1  r`-=<@[  
l,1,7 ^hZZ5(</8P  
l,3,6 ^1`Mz<  
!创建边坡面模型 [L(qrAQ2|z  
al,3,4,5,11                              !由线L3、L4、L5、L11生成一个面A1 <UsFBF  
al,1,2,11,6,10                            !由线L1、L2、L1、L6、L10生成一个面A2 =^ZDP1h/}  
al,7,8,9,10                              !由线L7、L8、L9、L10生成一个面A3 H. o=4[  
/PNUM, area, on                        !打开面号开关 9 n0 ?0mk  
!划分网格生成有限元模型 Y.yiUf/Q  
! 划分边坡围岩2网格 D`hl}  
!设置网格份数 'wyS9^F  
lsel,s,,,3,5,2                              !选择线L3和L5 AcHr X=O  
lesize,all,,,5                              !把所选择线分为5段 @qUgp*+{  
lsel,s,,,4,11,7                              !选择线L4和L11 vd{ban9  
lesize,all,,,26                              !把所选择线分为26段 uSQlE=  
mat,2                                    !给边坡围岩2赋予2号材料特性 " l;=jk]  
type,1                                  !采用单元类型1 ExKyjWAJ  
mshkey,0                                !设定自由网格划分 | -l9Z  
mshape,0                                !设定四边形网格划分 ^ruz-N^Y!  
amesh,1                                  !划分面积A1 ?*AhGza/  
! 划分边坡围岩1网格 \RT3#X+  
!设置网格份数 nQ+$  
lsel,s,,,2,6,4                              !选择线L2和L6 7H])2:)  
lsel,a,,,1                                  !附加选择线L1 "50 c<sZSB  
lesize,all,,,10                              !把所选择线分为10段 :;Xh`br  
lsel,s,,,8,10,2                              !选择线L8和L10 zQ~N(Jj?h  
lesize,all,,,16                              !把所选择线分为16段 tn 38T%  
lsel,s,,,7,9,2 utzf7?nIS  
lesize,all,,,12 3-Bz5sj9  
mat,1                                    !给边坡围岩1赋予1号材料特性 *q-VY[2  
type,1 M!R=&a=Z  
mshkey,0                                  !设定自由网格划分 X*w7q7\8-:  
mshape,0                                  !设定四边形网格划分 l~ Hu#+O  
amesh,2                                  !划分面积A2 .+) AeGh  
amesh,3                                  !划分面积A3 `a6AES'w$  
allsel _s18^7  
!保存网格模型 ]WYddiF  
save,Slope-grid.db u u$Jwn!S  
(3)施加约束和荷载 )ko[_OJj  
!给边坡模型施加约束 )k[{re  
!边坡两侧施加X方向约束  MFyi#nq  
nsel,s,loc,x,-800                            !选择X=-800线上所有节点 {zIcEN$ ~  
nsel,a,loc,x,1200                            !选择X=1200线上所有节点 Og-v][  
d,all,ux                                    !对所选择节点约束X方向位移 ]bP1gV(b-  
allsel }vxRjO,  
!边坡底部施加约束 &|fPskpy  
nsel,s,loc,y,-1200                          !选择Y=--1200线上所有节点 }D]y -BbA.  
d,all,uy                                  !对所选择节点约束X、Y方向位移 laRn![[  
d,all,ux                                  U<|*V5   
!施加重力加速度 qm/#kPlM  
acel,,9.8 r XT6u  
(4)求解 IOddu2.(  
/solu Hb3..o:  
!求解设置 MK.TBv  
antype,static                                  !设定为静力求解 "_0sW3rG  
nsubst,100                                  !设定最大子步数为100 9\Md.>  
pred,on                                      !打开时间步长预测器 MX"M2>"pT  
nropt,full                                    !设定牛顿-拉普森选项 OUhlQq\  
nlgeom,on                                    !打开大位移效果 6 \?GY  
lnsrch,on                                    !打开线性搜索 $sFqMy  
outres,all,all                                  !输出所有项 O/.8;.d;4Y  
cnvtol,f,,0.005,2,0.5                            !力收敛准则设定 1* ?XI  
cnvtol,u,,0.05,2,1                              !位移收敛准则设定 U!jRF  
!边坡在强度折减系数F=1时求解 o'_eLp  
allsel ]tNB^  
solve                                        !进行求解 cwE?+vB  
save,F1,db                                    !把F=1时求解结果保存 $}EI3a  
!边坡在强度折减系数F=1.2时求解 e$`;z%6y  
finish |RH^|2:x9Q  
/solu /Yj; '\3  
allsel ?aQVaw&L!7  
asel,s,area,,2,3,1                              !选择面积A2、A3 W 'PW;.,  
MPCHG,3,all                                ! 把所选择单元材料号改为3 vt#&YXu{A  
allsel =d$m@rc0r  
solve                                        !进行求解 Et'&}NjI  
save,F1.2,db                                  !把F=1.2时求解结果保存                                                                                  RF [81/w]  
!边坡在强度折减系数F=1.4时求解 I@kMM12>c  
finish ~zoZ{YqP  
/solu z(#CO<C.t  
allsel B A i ^t  
asel,s,area,,2,3,1                              !选择面积A2、A3 [5PQrf~Mo  
MPCHG,4,all                                !把所选择单元材料号改为4 Wb5n> *  
allsel 2;7n0LOs}  
solve                                        !进行求解 Bxf&gDwjgr  
save,F1.4,db                                  !把F=1.4时求解结果保存                            nTE\EZ+=2  
!边坡在强度折减系数F=1.6时求解 WM0-F@_  
finish ^YB3$:@$U  
/solu 8w ]'U  
allsel ~X-v@a  
asel,s,area,,2,3,1                              !选择面积A2、A3 F!RP *  
MPCHG,5,all                                !把所选择单元材料号改为5 -+0!Fkt@,  
allsel #iT3 aou  
solve qgxGq(6K  
save,F1.6,db                                  !把F=1.6时求解结果保存                                                                                  *iUR1V Y  
!边坡在强度折减系数F=1.8时求解 C_ W%]8u  
finish m'eM&1Ba  
/solu 2UMX%+ "J  
allsel :Vrj[i-{  
asel,s,area,,2,3,1                                !选择面积A2、A3 L4<=,}KS  
MPCHG,6,all                                  !把所选择单元材料号改为6 z J V>;  
allsel Dc3bG@K*G  
solve 7Q|v5@;pU  
save,F1.8,db                                  !把F=1.8时求解结果保存  s,j=Kym%  
!边坡在强度折减系数F=2.0时求解 +hIMfhF  
finish acI%fYw5p`  
/solu ,Qt2?  
allsel ?:7.3{|Aq  
asel,s,area,,2,3,1                                !选择面积A2、A3 #;\tgUQ  
MPCHG,7,all                                  !把所选择单元材料号改为7 A<-3u  
allsel 36d6KS 7  
solve eMjW^-RgE5  
save,F2.0,db                                  !把F=2.0时求解结果保存                            k^%F4d3z@C  
!边坡在强度折减系数F=2.2时求解 .G>6_n3  
finish AQs_(LR  
/solu FXbalQ?^  
allsel }q)dXFL=I#  
asel,s,area,,2,3,1                                !选择面积A2、A3 !Xzy:  
MPCHG,8,all                                  !把所选择单元材料号改为8 RCMO?CBe  
allsel `8D'r|=`Eh  
solve gFxaUrZA  
save,F2.2,db                                  !把F=2.2时求解结果保存                                                                                  H,)2Ou-Wn  
!边坡在强度折减系数F=2.4时求解 uXdR-@80*  
finish 1ifPc5j}  
/solu qbo W<W<H1  
allsel O*1la/~m  
asel,s,area,,2,3,1                                !选择面积A2、A3 xKXD`-|W  
MPCHG,9,all                                  !把所选择单元材料号改为9 6~Y`<#X5J  
allsel AE4>pzBe  
solve |E FbT>  
save,F2.4,db                                  !把F=2.4时求解结果保存  9U!#Y%*T  
!边坡在强度折减系数F=2.6时求解 `IK3e9QpcA  
finish b@[\+P] "  
/solu rlkg.e6  
allsel ZRYHsl{F+  
asel,s,area,,2,3,1                                !选择面积A2、A3 AA K}t6  
MPCHG,10,all                                  !把所选择单元材料号改为10 ETaLE[T%1  
allsel %M-B"#OB7  
solve AeEF/*  
save,F2.6,db                                  !把F=2.6时求解结果保存                            T0BM:ofx  
!边坡在强度折减系数F=2.8时求解 G=>LW1E|  
finish #po}Y  
/solu L)_L#]Yy  
allsel !{4bC  
asel,s,area,,2,3,1                                !选择面积A2、A3 Q".AmHn  
MPCHG,11,all                                !把所选择单元材料号改为11 UzUt=s!^H  
allsel %|3NCyJ*7  
solve R1\$}ep^  
save,F2.8,db                                  !把F=2.8时求解结果保存                                                                                  #}zL?s^G  
!边坡在强度折减系数F=3.0时求解 lR3JyYY{X  
finish ipJnNy;  
/solu d+g+ {p>?  
allsel H_7X%TvXb  
asel,s,area,,2,3,1                            !选择面积A2、A3 G `JXi/#`  
MPCHG,12,all                              !把所选择单元材料号改为12 >=ng?  
allsel sme!!+Rd  
solve m>-(c=3  
save,F3.0,db                                !把F=3.0时求解结果保存 g}'(V>(  
(5)后处理 }@jT-t]P  
/post1                                      !进入后处理 ,7izrf8  
!边坡在强度折减系数F=1时结果分析 M2y"M,k4  
Resume,'F1','db'                              !读入边坡在强度折减系数F=1时 WAp#[mW.fx  
set,1,last                                    !读入后一个子步 ch]Q%M  
pldisp,1                                      !绘制边坡模型变形图 X" ;ly0Mb  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 \B D'"  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 hG.~[#[&6  
!边坡在强度折减系数F=1.2时结果分析 %/=#8v4*  
Resume,'F1.2','db'                            !读入边坡在强度折减系数F=1.2时 oF#]<Z\  
set,1,last                                    !读入后一个子步 f m'Qif q^  
pldisp,1                                      !绘制边坡模型变形图 /(skIvE|  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 ){XG%nC  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 xnG,1doa  
!边坡在强度折减系数F=1.4时结果分析 ,a?$F1Z-  
Resume,'F1.4','db'                            !读入边坡在强度折减系数F=1.4时 @|J+ f5O  
set,1,last                                    !读入后一个子步 @d=4C{g%o  
pldisp,1                                      !绘制边坡模型变形图 r!J?Lc])8  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 |DMa2}%  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 i:,37INMt  
!边坡在强度折减系数F=1.6时结果分析 [lX3":)  
Resume,'F1.6','db'                            !读入边坡在强度折减系数F=1.6时 '}T6e1#JV  
set,1,last                                    !读入后一个子步 WjvD C"  
pldisp,1                                      !绘制边坡模型变形图 qlfYX8edZ  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 zq&,KZ  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 ~F9WR5}]  
!边坡在强度折减系数F=1.8时结果分析  e{33%5  
Resume,'F1.8','db'                            !读入边坡在强度折减系数F=1.8时 fF} NPl  
set,1,last                                    !读入后一个子步 '74-rL:i  
pldisp,1                                      !绘制边坡模型变形图 C/x<_VJzN/  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 j{u! /FD  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 ckPI^0A!  
!边坡在强度折减系数F=2.0时结果分析 #35@YMF  
Resume,'F2.0','db'                            !读入边坡在强度折减系数F=2.0时 o]|a5. O  
set,1,last                                    !读入后一个子步 #v`J]I)$  
pldisp,1                                      !绘制边坡模型变形图 f6Io|CZWJ  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 b_T?jCyW  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 naoH685R4  
!边坡在强度折减系数F=2.2时结果分析 ce\ F~8y  
Resume,'F2.2','db'                            !读入边坡在强度折减系数F=2.2时 E>3fk  
set,1,last                                    !读入后一个子步 UP}5Eh  
pldisp,1                                      !绘制边坡模型变形图 C) "|sG  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 TGe{NUO  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 GP\Pk/E  
!边坡在强度折减系数F=2.4时结果分析 ^UAL5}CQt  
Resume,'F2.4','db'                            !读入边坡在强度折减系数F=2.4时 -O-qEQd  
set,1,last                                    !读入后一个子步 T5+iX`#M  
pldisp,1                                      !绘制边坡模型变形图 x0!5z1KQh  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 @dj 2#  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 (C EXPf  
!边坡在强度折减系数F=2.6时结果分析 uk~4R@=&H  
Resume,'F2.6','db'                            !读入边坡在强度折减系数F=2.6时 -5 YvtL  
set,1,last                                    !读入后一个子步 XABI2Ex  
pldisp,1                                      !绘制边坡模型变形图 =?/RaK/ w  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 :LwNOuavN  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 r{gJ[%  
!边坡在强度折减系数F=2.8时结果分析 F! ;0eS"xp  
Resume,'F2.8','db'                            !读入边坡在强度折减系数F=2.8时 Z'~yUo=  
set,1,last                                    !读入后一个子步 a}+7MEUmZ/  
pldisp,1                                      !绘制边坡模型变形图 N{<=s]I%x  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 /Ma"a ^  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 &k+'TcWm  
!边坡在强度折减系数F=3.0时结果分析 ~"+Fp&[9f  
Resume,'F3.0','db'                            !读入边坡在强度折减系数F=3.0时 {D jz']  
set,1,last                                    !读入后一个子步 zp2IpYQ,3  
pldisp,1                                      !绘制边坡模型变形图 007SA6xq  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 .V?i3  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图
离线diaohuwang

发帖
273
土币
0
威望
0
原创币
0
只看该作者 1楼 发表于: 2011-05-23
楼主兄弟,Ansys用于岩土工程计算本来就不是很合适,更何况采用强度折减法啊,目前流行的这种思路本人不是很认同。
快速回复
限100 字节
温馨提示:欢迎交流讨论,请勿纯表情、纯引用!
 
上一个 下一个

      https://beian.mps.gov.cn/ 粤公网安备 44010602012919号 广州半山岩土网络科技有限公司 粤ICP备2024274469号

      工业和信息化部备案管理系统网站