做管线的动力时程分析，建模完成后进行瞬态分析显示没问题，但当加地震波之后不收敛，请大神帮忙看看什么问题，命令流如下，地震波文件见附件

/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,7800

k,1

k,2,0,0,30

l,1,2

circle,1,0.5

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

lsel,r,loc,z,0

lesize,all,,,4

lsel,s,length,,30

lesize,all,,,60

mshape,0,2d

mshkey,1

amesh,all

NGEN,2,1000,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,30 

TYPE, 5  

MAT,1

REAL, 5  

ESYS,  0   

SECNUM, 

TSHAP,LINE  

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

alls

NSEL,S, , ,1001,1976,1

D,all,all

nsel,s,loc,z,0

NSEL,A,LOC,Z,30

nsel,r,,,1,1000

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,10,2

kbc,1

acel,0,-9.8

solve

timint,on

/input,dizheng,txt

*do,i,1,NT

acel,0,AC(i),0

time,i*DT

nsubst,4,10,2

solve

*enddo

dizheng.txt