| 查看: 983 | 回复: 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 |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有17人回复
-
基金不中，共勉
已经有10人回复
-
青基如何破局
已经有6人回复
-
国自然面上复盘~欢迎讨论
已经有11人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有7人回复
-
基金未中，这种答复是模板吗？
已经有6人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有6人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有6人回复
-
怎么看青基中了没有啊
已经有5人回复
-
梦想
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
taiyanghef
铁虫 (著名写手)
- 应助: 5 (幼儿园)
- 金币: 4183.2
- 散金: 1142
- 红花: 7
- 帖子: 2065
- 在线: 359.3小时
- 虫号: 1176744
- 注册: 2010-12-24
- 性别: GG
- 专业: 机械结构强度学
2楼2017-03-14 18:46:02
