埋地管道动力时程分析没有问题和警告，计算结果位移非常小，应力也很小，感觉不太对，我用的是m，pa为单位，有大神能帮忙检查下，看是什么问题不，非常感谢，以下是完整的命令流，附件为地震数据文件。

/prep7

et,1,shell63

r,1,0.016,0.016,0.016,0.016,,,,,

et,2,combin39

keyopt,2,3,1

et,3,combin39

keyopt,3,3,2

et,4,combin39

keyopt,4,3,3

et,5,combin39

keyopt,5,3,3

r,2,-5e-3,-5e4,-1.5e-3,-5e4,0,0!定义横向弹簧本构

rmore,1.5e-3,5e4,5e-3,5e4  

r,3,-5e-3,-5e4,-1.5e-3,-5e4,0,0!定义轴向弹簧本构

rmore,1.5e-3,5e4,5e-3,5e4  

r,4,-4e-3,-6.9e4,-2.5e-3,-6.9e4,0,0!定义垂直弹簧本构

rmore,1.5e-3,1.667e4,5e-3,1.667e4  

r,5,0,0,0.05,600000,0.3,1.4e6!定义等效弹簧本构

mp,ex,1,2.1e11

mp,prxy,1,0.3

mp,dens,1,15600

k,1

k,2,0,0,15

l,1,2

circle,1,0.25

adrag,2,3,4,5,,,1

lsel,r,loc,z,0

lesize,all,,,4

lsel,s,length,,15

lesize,all,,,25

mshape,0,2d

mshkey,1

amesh,all

NGEN,2,10000,all, , , , , ,1,   

TYPE,   2   

MAT,       1

REAL,       2   

ESYS,       0   

SECNUM, 

TSHAP,LINE  

EINTF,0.0001, ,LOW,0,0,0,0,0

TYPE,3   

MAT,       1

REAL,3  

ESYS,       0   

SECNUM, 

TSHAP,LINE  

EINTF,0.0001, ,LOW,0,0,0,0,0

TYPE,4  

MAT, 1

REAL,4   

ESYS,0   

SECNUM, 

TSHAP,LINE  

EINTF,0.0001, ,LOW,0,0,0,0,0

alls

!两端弹簧

NSEL,S,LOC,Z,0

NSEL,A,LOC,Z,15

MAT,1

REAL, 5  

ESYS,  0   

SECNUM, 

TSHAP,LINE  

EINTF,0.0001, ,LOW,0,0,0,0,0

alls

NSEL,S, , ,10001,15000,1

D,all,all

nsel,s,loc,z,0

NSEL,A,LOC,Z,15

nsel,r,,,1,5000

d,all,ux,,,,,uy,roty,rotz

alls

*SET,NT,1000

*SET,DT,0.02

*DIM,AC,,NT,1

/solu

antype,trans

alphad,

betad,

timint,off

nlgeom,on

time,1e-6

nsubst,4

kbc,1

acel,0,9.8

solve

timint,on

NROPT,FULL,,

CNVTOL,F,,0.05,,0.01

CNVTOL,U,,0.05,,0.01

/input,dizhengbo,txt

*do,i,1,NT

acel,0,AC(i)*9.8,0

time,i*DT

nsubst,4

solve

*enddo

dizhengbo.txt