ANSYS隧洞蠕变分析求助

eureka

问一下各位老师，我想在隧道开挖过程中考虑围岩蠕变的影响，同时隧道开挖所需时间按“真实的时间”来进行，为什么蠕变计算的结果和不考虑蠕变效应的计算结果一模一样呢？是我的设置出了问题？这段命令流有什么问题吗？命令流有点长，谢谢各位老师了。

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

eureka

我把问题再详细的描述一下，一般情况下的隧道开挖过程，模拟是考虑围岩的非线性，混凝土为线弹性，在一般的开挖分析中，时间步是一个人为的数量，仅仅是作为计算从迭代到稳态的一种手段，而我现在需要考虑围岩的蠕变本构关系和时间效应，求解时间和时间步都应该是“真实的时间”，其中的蠕变参数由于暂时没有做相关的流变试验，采用的是所查阅的文献中的数据，应用的是隐式蠕变的3号方程。但是不知道为什么计算出来的结果还是和一般性（不考虑蠕变本构关系和时间效应）计算出来的结果一模一样，不知道哪里出了问题，请各位老师指导一下。谢谢了，这个问题困扰了我很久了。