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

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

上一主题 下一主题
离线fhzhw
 

发帖
426
土币
1699
威望
1245
原创币
0
只看楼主 倒序阅读 使用道具 楼主  发表于: 2008-05-12
各位朋友: 3%|&I:tI  
    大家好!我通过ansys用有限元强度折减法分析了一个边坡的稳定性,结果不同的工况下得出的塑性区是一样的,我不知道问题出在哪里,请各位看后给予解决问题的方案。不甚感激!附命令流和塑性区如下:(1)创建物理环境 uZYF(Yu  
/COM, Structural                              ! 指定结构分析 :kV#y  
/TITLE, Slope stability Analysis                  ! 定义工作标题 <=&`ZH   
/FILNAM,support,1                            ! 定义工作文件名 I{&[[7H  
(2)建立模型 QL/(72K  
!进入前处理器 rXq.DvQ  
/PREP7 FxY}m  
!定义单元类型 Hio0HL-  
ET,1,PLANE82                              !定义边坡围岩单元 S+6.ZZ9c  
KEYOPT,1,3,2                                !求解类型为平面应变 (E 3b\lST  
!定义材料属性 nW:C/{n2tG  
!边坡围岩1材料属性  =%O6:YM   
MP,EX,1,3E10  fbvL7* (  
MP,PRXY,1,0.25    ~=LE0.3[  
MP,DENS,1,2500                            !定义塑性模型 3M`M  
TB,DP,1                                    !定义数据表1 1.hyCTnI  
TBDATA,1,9E5,42.8                          !给数据表输入数据  > |=ts  
!边坡围岩2材料属性                        !采用弹性模型 W*:.Gxv]  
MP,EX,2,3.2E10                              !定义弹性模量 ;;t yoh~t  
MP,PRXY,2,0.24                              !定义泊松比 4"ZP 'I;  
MP,DENS,2,2700                            !定义密度 LOYk9m  
!取强度折减系数F=1.2时边坡围岩材料属性 G!##X: 6'  
MP,EX,3,3E10  G.B2('  
MP,PRXY,3,0.25    Rv>-4@fMJ  
MP,DENS,3,2500 Q{>k1$fkV  
TB,DP,3                                    !定义数据表3 wQf-sk#  
TBDATA,1,7.5E5,37.7                        !给数据表输入数据 ?z+eWL  
!取强度折减系数F=1.4时边坡围岩材料属性 y Pp9\[+^j  
MP,EX,4,3E10  d-ko ^Y0  
MP,PRXY,4,0.25                              !定义泊松比 j;r-NCBnz  
MP,DENS,4,2500 {Xy5pfW Q  
TB,DP,4                                    !定义数据表4 M3y NAN  
TBDATA,1,6.4E5,33.5                        !给数据表输入数据 _LnpnL:  
!取强度折减系数F=1.6时边坡围岩材料属性 u#~RkY7s  
MP,EX,5,3E10  |e&\<LwsP  
MP,PRXY,5,0.25    3}1u\(Mf  
MP,DENS,5,2500                            !定义密度 pki%vRY  
TB,DP,5                                    !定义数据表5 c6]U E@A  
TBDATA,1,5.6E5,30                          !给数据表输入数据 ^76]0`gS  
!取强度折减系数F=1.8时边坡围岩材料属性 re<{ >  
MP,EX,6,3E10                              !定义弹性模量 N"R]Yp;j  
MP,PRXY,6,0.25    ;(%QD 3>  
MP,DENS,6,2500 H?Wya.7  
TB,DP,6                                    !定义数据表6 IOH}x4  
TBDATA,1,5.0E5,27.2 ]EAO+x9  
!取强度折减系数F=2.0时边坡围岩材料属性 l'qg8  
MP,EX,7,3E10  UkC!1Jy  
MP,PRXY,7,0.25                            !定义泊松比 $PPi5f}HD  
MP,DENS,7,2500 Zi i   
TB,DP,7                                    !定义数据表7 & .j&0WE  
TBDATA,1,4.5E5,24.8 :cECRm*  
!取强度折减系数F=2.2时边坡围岩材料属性 +sA2WK]  
MP,EX,8,3E10  |df Pki{  
MP,PRXY,8,0.25    33q}CzK  
MP,DENS,8,2500 <lPG=Xt  
TB,DP,8                                    !定义数据表8 _VN?#J)o  
TBDATA,1,4.09E5,22.8 6 "sSoj  
!取强度折减系数F=2.4时边坡围岩材料属性 N+xP26D8  
MP,EX,9,3E10  @.l@\4m  
MP,PRXY,9,0.25    "S]TP$O D  
MP,DENS,9,2500 jr. "I+  
TB,DP,9                                      !定义数据表9 /6* 42[r  
TBDATA,1,3.6E5,21.1                          !给数据表输入数据 e !Y~Qy  
!取强度折减系数F=2.6时边坡围岩材料属性 P@B]  
MP,EX,10,3E10                                !定义弹性模量  x9g#<2w8  
MP,PRXY,10,0.25                              !定义泊松比 X_h}J=33Q  
MP,DENS,10,2500                          !定义密度 cI*;k.KU  
TB,DP,10                                  !定义数据表10 2zb"MEOS5  
TBDATA,1,3.46E5,19.6                        !给数据表10输入数据 ki!0^t:9  
!取强度折减系数F=2.8时边坡围岩材料属性 t*u:hex  
MP,EX,11,3E10                              !定义弹性模量 )10+@d  
MP,PRXY,11,0.25                            !定义泊松比 * u>\57W  
MP,DENS,11,2500 7$=In K  
TB,DP,11                                  !定义数据表11 KpGhQdR#  
TBDATA,1,3.2E5,18.2                        !给数据表11输入数据 niyV8v  
!取强度折减系数F=3.0时边坡围岩材料属性 Hn"RH1Zy  
MP,EX,12,3E10                              !定义弹性模量 q} >%8;nm  
MP,PRXY,12,0.25    F41=b4/  
MP,DENS,12,2500                            !定义密度 n>YKa)|W`  
TB,DP,12                                  !定义数据表12 H <l7ZS:  
TBDATA,1,3.0E5,17.1                        !给数据表12输入数据 TJRCH>E[a  
!建立几何模型 0h_|t-9j  
!创建边坡线模型线 T8g$uFo  
K,1,,,,                                      ! 创建关键点1 `;C  V=,M  
K,2,-800,,, H{wl% G  
K,3,-800,-800,,  _DtV  
K,4,-800,-1200,, wHy!CP%  
K,5,1200,-1200,,  :I#V.  
K,6,1200,-800,,  ez$(c  
K,7,1200,,, 1+s;FJ2}  
K,8,1200,400,, [W&T(%(W-  
K,9,400,400,, k(HUUH_z  
l,1,2                                        !连接关键点1、2生成直线L1 WsB?C&>x  
l,2,3 7[)E>XRE  
l,3,4 qE"OB  
l,4,5 hVY$;s  
l,5,6 2+XA X:YD  
l,6,7 })%{AfDRF  
l,7,8 MQ2}EY*A  
l,8,9 9.#<b |g  
l,9,1 AK4t\D)K1  
l,1,7 <hyKu  
l,3,6 GbI/4<)l}  
!创建边坡面模型 59u }W 0  
al,3,4,5,11                              !由线L3、L4、L5、L11生成一个面A1 g_bLl)g<  
al,1,2,11,6,10                            !由线L1、L2、L1、L6、L10生成一个面A2 6gDN`e,@  
al,7,8,9,10                              !由线L7、L8、L9、L10生成一个面A3 L4W5EO$  
/PNUM, area, on                        !打开面号开关 Pm7}"D'/  
!划分网格生成有限元模型 /n&&Um\  
! 划分边坡围岩2网格 FSO).=#  
!设置网格份数 SXh-A1t  
lsel,s,,,3,5,2                              !选择线L3和L5 wCBplaojJ  
lesize,all,,,5                              !把所选择线分为5段 %|oym.-I6  
lsel,s,,,4,11,7                              !选择线L4和L11 ?@x/E&  
lesize,all,,,26                              !把所选择线分为26段 ;uGv:$([g  
mat,2                                    !给边坡围岩2赋予2号材料特性 :3 mh@[V  
type,1                                  !采用单元类型1 OyIw>Wfv  
mshkey,0                                !设定自由网格划分 Kg]J/|0\  
mshape,0                                !设定四边形网格划分 DEgXQ[  
amesh,1                                  !划分面积A1 -9?]IIVb  
! 划分边坡围岩1网格 QT}tvm@PMq  
!设置网格份数 HzsdHH(J  
lsel,s,,,2,6,4                              !选择线L2和L6 x8|J-8A(  
lsel,a,,,1                                  !附加选择线L1 y~V(aih}D  
lesize,all,,,10                              !把所选择线分为10段 *-X[u:  
lsel,s,,,8,10,2                              !选择线L8和L10 8ao_i=&x  
lesize,all,,,16                              !把所选择线分为16段 ca9X19NG  
lsel,s,,,7,9,2 :`sUt1Fw.  
lesize,all,,,12 \;Weizq5  
mat,1                                    !给边坡围岩1赋予1号材料特性 kJR`:J3DJ  
type,1 %C]>9."  
mshkey,0                                  !设定自由网格划分 Fr-SvsNFB  
mshape,0                                  !设定四边形网格划分 (8OsGn  
amesh,2                                  !划分面积A2 *Pr )%  
amesh,3                                  !划分面积A3 BA.uw_^4  
allsel /-s6<e!  
!保存网格模型 |s_GlJV.  
save,Slope-grid.db ALHIGJW:6$  
(3)施加约束和荷载 6]wIG$j  
!给边坡模型施加约束 * y,v}-  
!边坡两侧施加X方向约束 !,PWb3S  
nsel,s,loc,x,-800                            !选择X=-800线上所有节点 j>kqz>3  
nsel,a,loc,x,1200                            !选择X=1200线上所有节点 +ZV5o&V>  
d,all,ux                                    !对所选择节点约束X方向位移 (@YG~ 0  
allsel :(%5:1W  
!边坡底部施加约束 ]M3yLYK/P  
nsel,s,loc,y,-1200                          !选择Y=--1200线上所有节点 zuCSj~  
d,all,uy                                  !对所选择节点约束X、Y方向位移 nk:)j:fr  
d,all,ux                                  O6Y0XL  
!施加重力加速度 j<$2hiI/?&  
acel,,9.8 I-)4YQI  
(4)求解 cwL_tq  
/solu Q2> gU#  
!求解设置 : Dp0?&_  
antype,static                                  !设定为静力求解 pR_9NfV{  
nsubst,100                                  !设定最大子步数为100 d;>QhoiL  
pred,on                                      !打开时间步长预测器 Bw.i}3UT6  
nropt,full                                    !设定牛顿-拉普森选项 Ys7]B9/1O  
nlgeom,on                                    !打开大位移效果 Y);=TM6s  
lnsrch,on                                    !打开线性搜索 $[|mGae  
outres,all,all                                  !输出所有项 Hr C+Yjp  
cnvtol,f,,0.005,2,0.5                            !力收敛准则设定 Vvo 7C!$z  
cnvtol,u,,0.05,2,1                              !位移收敛准则设定 6u%&<")4HP  
!边坡在强度折减系数F=1时求解 pCG}Z Ka  
allsel i#O SC5ZI  
solve                                        !进行求解 <3 uNl  
save,F1,db                                    !把F=1时求解结果保存 '%;m?t% q  
!边坡在强度折减系数F=1.2时求解 naNghGQ  
finish &,/ S`ke=  
/solu P<-@h1p,  
allsel TA\vZGJ('  
asel,s,area,,2,3,1                              !选择面积A2、A3 MK*r+xfSae  
MPCHG,3,all                                ! 把所选择单元材料号改为3 l-3~K-k<@  
allsel xD7]C|8o  
solve                                        !进行求解 ?]5qr?W%  
save,F1.2,db                                  !把F=1.2时求解结果保存                                                                                  |y*c9  
!边坡在强度折减系数F=1.4时求解 \7_y%HR  
finish EX"yxZ~  
/solu ^rz_f{c]-  
allsel )%]J>&/0J  
asel,s,area,,2,3,1                              !选择面积A2、A3 D+c>F5  
MPCHG,4,all                                !把所选择单元材料号改为4 <ZW-QN4  
allsel YkADk9fE  
solve                                        !进行求解 +|f@^-  
save,F1.4,db                                  !把F=1.4时求解结果保存                            iDD$pd,e\  
!边坡在强度折减系数F=1.6时求解 fV~~J2IK  
finish u@444Vzg  
/solu QWU-m{@~&  
allsel |[lKY+26:{  
asel,s,area,,2,3,1                              !选择面积A2、A3 UZ";a453r  
MPCHG,5,all                                !把所选择单元材料号改为5 m[2gdJK  
allsel ig"L\ C"T  
solve tX[WH\(xI  
save,F1.6,db                                  !把F=1.6时求解结果保存                                                                                  #Q5o)x  
!边坡在强度折减系数F=1.8时求解 MOC/KNb  
finish YZ7.1`8  
/solu d=^z`nt !R  
allsel v^P O|Z  
asel,s,area,,2,3,1                                !选择面积A2、A3 {G-kNU  
MPCHG,6,all                                  !把所选择单元材料号改为6 cb bFw  
allsel d5-qZ{W  
solve WzWX E(  
save,F1.8,db                                  !把F=1.8时求解结果保存  yQrD9*t&g  
!边坡在强度折减系数F=2.0时求解 (% 9$!v{3  
finish 0{mex4  
/solu Zd&S@Z  
allsel kT=8e;K  
asel,s,area,,2,3,1                                !选择面积A2、A3 `_h&glMJ,q  
MPCHG,7,all                                  !把所选择单元材料号改为7 DwF hK*  
allsel #E]59_  
solve 4K74=r),i  
save,F2.0,db                                  !把F=2.0时求解结果保存                            fy$1YI>!Q  
!边坡在强度折减系数F=2.2时求解 d5d@k  
finish R-Sym8c  
/solu -qoH,4w  
allsel JI5Dy>u:  
asel,s,area,,2,3,1                                !选择面积A2、A3 ^@]3R QB  
MPCHG,8,all                                  !把所选择单元材料号改为8 LQ% `c  
allsel =H~j,K  
solve }WV:erg`  
save,F2.2,db                                  !把F=2.2时求解结果保存                                                                                  ^Dx&|UwiZa  
!边坡在强度折减系数F=2.4时求解 w =KPT''!  
finish GthYzd:'hJ  
/solu mc3"`+o  
allsel .(vwIb8\_  
asel,s,area,,2,3,1                                !选择面积A2、A3 %)wjR/o  
MPCHG,9,all                                  !把所选择单元材料号改为9 ?Ob3tUz2  
allsel EV?z`jE9  
solve <Gsu Z  
save,F2.4,db                                  !把F=2.4时求解结果保存  PW4q~rc=:  
!边坡在强度折减系数F=2.6时求解 0$njMnB2l  
finish #;<Y[hR{P  
/solu ~ K=b\xc^  
allsel j!ch5A  
asel,s,area,,2,3,1                                !选择面积A2、A3 8ipez/  
MPCHG,10,all                                  !把所选择单元材料号改为10 Debv4Gr;^  
allsel f 1d?.)  
solve bZ6+,J  
save,F2.6,db                                  !把F=2.6时求解结果保存                            @XVTU  
!边坡在强度折减系数F=2.8时求解 cnLro  
finish  3CJwj  
/solu cNH7C"@GVu  
allsel {b{s<@?  
asel,s,area,,2,3,1                                !选择面积A2、A3 c%&>p||  
MPCHG,11,all                                !把所选择单元材料号改为11 =Sv/IXX\di  
allsel <uJ@:oWG7  
solve qWw=8Bq  
save,F2.8,db                                  !把F=2.8时求解结果保存                                                                                  YzWz|  
!边坡在强度折减系数F=3.0时求解 Y ay?=Y{  
finish P*o9a  
/solu <}LC~B!  
allsel q*KAk{kR(v  
asel,s,area,,2,3,1                            !选择面积A2、A3 qb4z T  
MPCHG,12,all                              !把所选择单元材料号改为12 Je{ykL?N  
allsel v2?ZQeHr_(  
solve h$*!8=M  
save,F3.0,db                                !把F=3.0时求解结果保存 U~8g_*  
(5)后处理 F\! `/4  
/post1                                      !进入后处理 u4j5w  
!边坡在强度折减系数F=1时结果分析 Q20 %"&Xp]  
Resume,'F1','db'                              !读入边坡在强度折减系数F=1时 h\e.e3/  
set,1,last                                    !读入后一个子步 f5r0\7y0  
pldisp,1                                      !绘制边坡模型变形图 D]}G.v1  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 {8OCXus3m  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 :-'qC8C  
!边坡在强度折减系数F=1.2时结果分析 26nx`w?j(  
Resume,'F1.2','db'                            !读入边坡在强度折减系数F=1.2时 f<H2-(m  
set,1,last                                    !读入后一个子步 ZW}_DT0  
pldisp,1                                      !绘制边坡模型变形图 8_8l.!~  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 c|%6e(g"L  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 ;aBG,dr}i  
!边坡在强度折减系数F=1.4时结果分析 ]tD]Wx%  
Resume,'F1.4','db'                            !读入边坡在强度折减系数F=1.4时 v1[29t<I!  
set,1,last                                    !读入后一个子步 >/|*DI-HJ  
pldisp,1                                      !绘制边坡模型变形图 o@Oqm>]SS  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 HDLk>_N_s,  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 "fI6Cpc  
!边坡在强度折减系数F=1.6时结果分析 0mnw{fE8_  
Resume,'F1.6','db'                            !读入边坡在强度折减系数F=1.6时 2|L&DF:G  
set,1,last                                    !读入后一个子步 Y 7aqO5  
pldisp,1                                      !绘制边坡模型变形图 TN.rrop`#g  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 /\Ef%@  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 P'2Qen*  
!边坡在强度折减系数F=1.8时结果分析 7"mc+QOp  
Resume,'F1.8','db'                            !读入边坡在强度折减系数F=1.8时 w &(ag$p'  
set,1,last                                    !读入后一个子步 P%6~&woF  
pldisp,1                                      !绘制边坡模型变形图 : 'c&,oLY  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 z}@7'_iJ  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 -%~4W?  
!边坡在强度折减系数F=2.0时结果分析 Ngwb Q7)  
Resume,'F2.0','db'                            !读入边坡在强度折减系数F=2.0时 s>en  
set,1,last                                    !读入后一个子步 H.c7Nle  
pldisp,1                                      !绘制边坡模型变形图 u"8yK5!  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 K8~d^G  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 y^k$Us  
!边坡在强度折减系数F=2.2时结果分析 KP"+e:a%  
Resume,'F2.2','db'                            !读入边坡在强度折减系数F=2.2时 j6YOKJX  
set,1,last                                    !读入后一个子步 G3 m Z($y  
pldisp,1                                      !绘制边坡模型变形图 "z c l|@  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 nEfK53i_  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 <[v[ci  
!边坡在强度折减系数F=2.4时结果分析  <Uur^uB  
Resume,'F2.4','db'                            !读入边坡在强度折减系数F=2.4时 9z0p5)]n>  
set,1,last                                    !读入后一个子步 phK/   
pldisp,1                                      !绘制边坡模型变形图 d1*<Ll9K  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 pJ=#zsE0  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 } #J/fa9 !  
!边坡在强度折减系数F=2.6时结果分析 :Al!1BJQ  
Resume,'F2.6','db'                            !读入边坡在强度折减系数F=2.6时 5bIw?%dk(  
set,1,last                                    !读入后一个子步 SKtrtm  
pldisp,1                                      !绘制边坡模型变形图 !<oe=)Iz|  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 2/f}S?@   
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 : +u]S2u{  
!边坡在强度折减系数F=2.8时结果分析 &L:!VL{I  
Resume,'F2.8','db'                            !读入边坡在强度折减系数F=2.8时 GVz6-T~\>  
set,1,last                                    !读入后一个子步 h 0|s  
pldisp,1                                      !绘制边坡模型变形图 D@KlOU{<  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 pw#-_  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图 [RL9>n8f  
!边坡在强度折减系数F=3.0时结果分析 ,_P-$lB  
Resume,'F3.0','db'                            !读入边坡在强度折减系数F=3.0时 BWNi [^]  
set,1,last                                    !读入后一个子步 lFk R=!?=  
pldisp,1                                      !绘制边坡模型变形图 H::bwn`Vc  
plnsol,u,x                                    !绘制边坡模型水平方向位移云图 pdMc}=K  
plnsol,eppl,eqv                              !绘制边坡模型塑性应变云图
离线diaohuwang

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

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

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