| 查看: 911 | 回复: 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 |
回复此楼» 猜你喜欢
博士自荐
已经有6人回复
-
博士推荐
已经有4人回复
-
求环氧树脂研发1名
已经有10人回复
-
280求调剂
已经有5人回复
-
什么是人一生最重要的？
已经有10人回复
-
面上可以超过30页吧？
已经有13人回复
-
为什么中国大学工科教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？
已经有13人回复
-
版面费该交吗
已经有17人回复
-
【博士招生】太原理工大学2026化工博士
已经有8人回复
» 本主题相关商家推荐: (我也要在这里推广)
taiyanghef
铁虫 (著名写手)
- 应助: 5 (幼儿园)
- 金币: 4257.7
- 散金: 1028
- 红花: 7
- 帖子: 1989
- 在线: 356小时
- 虫号: 1176744
- 注册: 2010-12-24
- 性别: GG
- 专业: 机械结构强度学
2楼2017-03-14 18:46:02