| 查看: 967 | 回复: 1 | |||
[交流]
ansys隧洞蠕变分析 已有1人参与
|
|
问一下各位老师,我想在隧道开挖过程中考虑围岩蠕变的影响,同时隧道开挖所需时间按“真实的时间”来进行,为什么蠕变计算的结果和不考虑蠕变效应的计算结果一模一样呢?是我的设置出了问题?这段命令流有什么问题吗?命令流有点长,谢谢各位老师了。 FINISH /CLEAR !修改自dxf to ansys !几何建模 /PREP7 K, 1 ,-5.422,3.015,.000, K,2 ,-5.949,-1.820,.000, K, 3 ,-5.931,1.874,.000, LARC,1 ,2 ,3 K, 4 ,-4.510,-3.299,.000, K, 5 ,-5.730,-2.307,.000, LARC,2 ,4 ,5 K, 6 ,5.949,-1.820,.000, K, 7 ,5.422,3.015,.000, K, 8 ,6.200,-.596,.000, LARC,6 ,7 ,8 K, 9 ,4.510,-3.299,.000, K, 10 ,4.990,-3.066,.000, LARC,9 ,6 ,10 K, 11 ,.000,6.145,.000, K, 12 ,4.427,4.312,.000, LARC,7 ,11 ,12 K, 13 ,.000,-4.059,.000, K, 14 ,-3.410,-3.630,.000, LARC,4 ,13 ,14 K, 15 ,-50.000,45.000,.000, K, 16 ,-50.000,13.161,.000, L, 15 ,16 K, 17 ,50.000,45.000,.000, K, 18 ,50.000,13.161,.000, L,17 ,18 K, 19 ,50.000,-13.161,.000, K, 20 ,50.000,-45.000,.000, L, 19 ,20 K, 21 ,-50.000,-13.161,.000, K, 22 ,-50.000,-45.000,.000, L, 21 ,22 K, 23 ,14.624,13.161,.000, K, 24 ,14.624,-13.161,.000, L, 23 ,24 K, 25 ,-14.624,-13.161,.000, K, 26 ,-14.624,13.161,.000, L, 25 ,26 L, 21 ,25 K, 27 ,-50.000,45.000,.000, K, 28 ,-14.624,45.000,.000, L, 27 ,28 K, 29 ,-14.624,-45.000,.000, L, 22 ,29 L, 26 ,28 K, 30 ,14.624,45.000,.000, L, 23 ,30 L, 26 ,16 L, 23 ,18 L, 24 ,19 L, 25 ,29 K, 31 ,14.624,-45.000,.000, L, 24 ,31 L, 30 ,17 L, 31 ,20 L, 16 ,21 L, 18 ,19 L, 7 ,23 L, 1 ,26 L, 2 ,25 L, 6 ,24 K, 32 ,.000,13.161,.000, L, 11 ,32 L, 26 ,32 L, 32 ,23 K, 33 ,.000,-13.161,.000, L, 13 ,33 L, 25 ,33 L, 33 ,24 K, 34 ,.000,45.000,.000, L, 32 ,34 K, 35 ,.000,-45.000,.000, L, 33 ,35 L, 28 ,34 L, 34 ,30 L, 29 ,35 L, 35 ,31 K, 36 ,-1.620,5.932,.000, LARC,11 ,1 ,36 K, 37 ,1.147,-4.011,.000, LARC,13 ,9 ,37 L, 11 ,13 ALLSEL LPLOT /REPLOT !合并线 FLST,2,2,4,ORDE,2 FITEM,2,2 FITEM,2,6 LCOMB,P51X FLST,2,2,4,ORDE,2 FITEM,2,4 FITEM,2,44 LCOMB,P51X NUMMRG,KP LPLOT /REPLOT !生成面 AL,1,2,43,45 AL,3,4,5,45 AL,1,12,28,29 AL,43,28,31,32 AL,5,27,31,33 AL,3,27,30,11 AL,4,30,34,36 AL,2,29,34,35 AL,7,14,18,16 AL,16,39,37,32 AL,33,17,40,37 AL,19,8,23,17 AL,25,18,12,13 AL,11,19,26,20 AL,10,13,21,15 AL,35,21,41,38 AL,36,38,42,22 AL,20,22,24,9 !几何参数设定 *AFUN,DEG !定义单元类型 ET,1,MESH200,2 !辅助画网格(3D线单元2结点) ET,2,MESH200,6 !辅助画网格(3D面单元4结点) ET,3,SHELL181 ET,4,SOLID45 R,1,0.4 MP,EX,1,2.8E10 !衬砌材料(C25混凝土) MP,PRXY,1,0.2 MP,DENS,1,2400 MP,EX,2,2.5E10 !土层 MP,PRXY,2,0.3 MP,DENS,2,2500 TB,DP,2,,, TBMODIF,1,1,1.2E6 TBMODIF,1,2,42 TBMODIF,1,3,21 TB,CREE,2,1,5,11 TBTEMP,293 TBDATA,,1.3E-12,1.04,-0.8,0,1.3E-19,1.04,0 MP,EX,3,2.5E10 !被挖土层 MP,PRXY,3,0.3 MP,DENS,3,2500 TB,DP,3,,, TBMODIF,1,1,1.2E6 TBMODIF,1,2,42 TBMODIF,1,3,21 TB,CREE,3,1,5,3 TBTEMP,293 TBDATA,,1.3E-12,1.04,-0.8,0,1.3E-19,1.04,0 SAVE !网格划分 ALLSEL FLST,2,6,4,ORDE,6 FITEM,2,1 FITEM,2,3 FITEM,2,11 FITEM,2,12 FITEM,2,25 FITEM,2,26 LESIZE,P51X,,,12 LESIZE,45,,,15 ALLSEL FLST,2,7,4,ORDE,7 FITEM,2,1 FITEM,2,3 FITEM,2,11 FITEM,2,12 FITEM,2,25 FITEM,2,26 FITEM,2,45 LSEL,U,,,P51X LESIZE,ALL,,,10 !进行面划分 TYPE,2 AMESH,1 AMESH,2 ALLSEL TYPE,2 MSHKEY,1 MSHAPE,0,2D FLST,2,16,5,ORDE,2 FITEM,2,3 FITEM,2,-18 AMESH,P51X ALLSEL !将模型进行三维拓展 K,1000,-5.949,-1.820,-30 L,2,1000 /VIEW,1,1,1,1 /REP LPLOT !衬砌结构 EXTOPT,ESIZE,10,0, FLST,2,6,4,ORDE,3 FITEM,2,1 FITEM,2,-5 FITEM,2,43 ADRAG,P51X,,,,,,6 APLOT GPLOT !显示全部元素 TYPE,3 REAL,1 MAT,1 ASEL,S,LOC,Z,-15 APLOT LSEL,S,LOC,Z,-15 LESIZE,ALL,,,10 MSHAPE,0,2D MSHKEY,1 AMESH,ALL !创建地层单元 !未挖地层 ASEL,INVE APLOT EXTOPT,ESIZE,10,0 EXTOPT,ACLEAR,1 !拉伸之后将面单元清除 TYPE,4 MAT,2 ASEL,R,AREA,,3,18,1 VDRAG,ALL,,,,,,6 ALLSEL MAT,3 FLST,2,2,5,ORDE,2 FITEM,2,1 FITEM,2,2 VDRAG,P51X,,,,,,6 EPLOT NUMMRG,NODE NUMMRG,KP NUMMRG,ELEM FINISH !设置边界条件及重力加速度 /SOLU ANTYPE,0 ASEL,S,LOC,X,-50 ASEL,A,LOC,X,50 DA,ALL,UX,0 ALLSEL ASEL,S,LOC,Y,-45 DA,ALL,UY,0 !UY=0 ALLS ASEL,S,LOC,Z,-30 ASEL,A,LOC,Z,0 DA,ALL,UZ,0 ALLSEL ACEL,,10 !设置分析选项 RATE,OFF NSUBST,100,1000,5,1 AUTOTS,ON PRED,ON LNSRCH,ON NLGEOM,ON NROPT,FULL CNVTOL,F,,0.02,2,0.5 !进行初始地应力的计算 TIME,1E-8 ESEL,S,TYPE,,3 EKILL,ALL ESEL,ALL ESEL,S,LIVE NSLE,S NSEL,INVERT D,ALL,ALL NSEL,ALL ESEL,ALL /PBC,ALL,,1 GPLOT SOLVE !第一步开挖 /SOLU RATE,ON TIME,23400 NSUBST,100,1000,5,1 !杀死被挖土体单元 ALLSEL ESEL,S,MAT,,3 NSLE,S NSEL,R,LOC,Z,0.001,-3.001 ESLN,R,1 EKILL,ALL !约束杀死单元的节点的全部自由度 ESEL,ALL ESEL,S,LIVE NSLE,S NSEL,INVE D,ALL,ALL NSEL,ALL ESEL,ALL /PBC,ALL,,1 GPLOT SOLVE |
» 猜你喜欢
售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急
已经有11人回复
售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急
已经有14人回复
面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇!
已经有26人回复
面上提前没消息,有中的吗
已经有18人回复
纯娱乐,不喜欢勿喷
已经有9人回复
售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急
已经有6人回复
售SCI文章,我:8O.5.5.1O.54,科目全,可十急
已经有7人回复
售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急
已经有8人回复
» 本主题相关商家推荐: (我也要在这里推广)
taiyanghef
铁虫 (著名写手)
- 应助: 5 (幼儿园)
- 金币: 4183.2
- 散金: 1142
- 红花: 7
- 帖子: 2065
- 在线: 359.3小时
- 虫号: 1176744
- 注册: 2010-12-24
- 性别: GG
- 专业: 机械结构强度学
2楼2017-03-14 18:46:02










回复此楼
50