24小时热门版块排行榜

返回列表

wanwan1993

新虫 (初入文坛)

应助: 0 (幼儿园)
帖子: 3
在线: 1.3小时
虫号: 4949996
注册: 2016-08-28

[交流] 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

回复此楼