LS-DANA后处理中位移-时间图获取方法

请教： 我是做应力波在混凝土中传播检测其缺陷，不知道什么原因，我在LS-DANA后处理中提取不出来节点位移-时间图，提取出来的位移-时间图显示位移是0，但是能够提取出速度-时间图和加速度-时间图。 模型的命令流附下表
/prep7
et,1,plane162
mp,dens,1,2500
mp,ex,1,40e9
mp,nuxy,1,0.2
mp,ex,2,2e11
mp,dens,2,7800
mp,nuxy,2,0.3
RECTNG,0.028,0.23,0,0.028,
RECTNG,0.23,0.25,0,0.25,
RECTNG,0,0.23,0.028,0.15,
RECTNG,0,0.23,0.15,0.151,
RECTNG,0,0.23,0.151,0.152,
RECTNG,0,0.23,0.152,0.25,
FLST,3,1,8
FITEM,3,0,0,0
K, ,P51X
FLST,2,3,3
FITEM,2,25
FITEM,2,9
FITEM,2,4
A,P51X
CYL4,0,0,0.024,45, ,90
ASBA, 7, 8
NUMMRG,ALL, , , ,LOW
NUMCMP,ALL
FLST,5,15,4,ORDE,9
FITEM,5,1
FITEM,5,3
FITEM,5,-9
FITEM,5,11
FITEM,5,-12
FITEM,5,14
FITEM,5,17
FITEM,5,19
FITEM,5,-21
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1,0.002, , , , , , ,1
!*
MSHAPE,0,2D
MSHKEY,0
!*
FLST,5,6,5,ORDE,2
FITEM,5,1
FITEM,5,-6
CM,_Y,AREA
ASEL, , , ,P51X
CM,_Y1,AREA
CHKMSH,'AREA'
CMSEL,S,_Y
!*
AMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
FLST,5,6,4,ORDE,4
FITEM,5,1
FITEM,5,4
FITEM,5,8
FITEM,5,-11
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1,0.002, , , , , , ,1
!*
FLST,5,2,5,ORDE,2
FITEM,5,1
FITEM,5,3
CM,_Y,AREA
ASEL, , , ,P51X
CM,_Y1,AREA
CHKMSH,'AREA'
CMSEL,S,_Y
!*
!*
ACLEAR,_Y1
AMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
FLST,5,4,4,ORDE,2
FITEM,5,1
FITEM,5,-4
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1,0.002, , , , , , ,1
!*
CM,_Y,AREA
ASEL, , , , 1
CM,_Y1,AREA
CHKMSH,'AREA'
CMSEL,S,_Y
ACLEAR,_Y1
AMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
FLST,5,2,4,ORDE,2
FITEM,5,15
FITEM,5,18
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1,0.002, , , , , , ,1
FLST,5,2,4,ORDE,2
FITEM,5,13
FITEM,5,16
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1,0.001, , , , , , ,1
FLST,5,2,5,ORDE,2
FITEM,5,4
FITEM,5,-5
CM,_Y,AREA
ASEL, , , ,P51X
CM,_Y1,AREA
CHKMSH,'AREA'
CMSEL,S,_Y
ACLEAR,_Y1
AMESH,_Y1
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
FLST,5,2,4,ORDE,2
FITEM,5,22
FITEM,5,-23
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,14, , , , ,1
FLST,5,2,4,ORDE,2
FITEM,5,24
FITEM,5,-25
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,8,1.1, , , ,0
CM,_Y,AREA
ASEL, , , , 7
CM,_Y1,AREA
CHKMSH,'AREA'
CMSEL,S,_Y
!*
AMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
TYPE, 1
MAT, 2
REAL,
ESYS, 0
FLST,5,2,5,ORDE,2
FITEM,5,4
FITEM,5,-5
CM,_Y,AREA
ASEL, , , ,P51X
CM,_Y1,AREA
CHKMSH,'AREA'
CMSEL,S,_Y
ACLEAR,_Y1
AMESH,_Y1
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
FLST,3,6,5,ORDE,2
FITEM,3,1
FITEM,3,-6
ARSYM,X,P51X, , , ,0,0
FLST,3,12,5,ORDE,4
FITEM,3,1
FITEM,3,-6
FITEM,3,8
FITEM,3,-13
ARSYM,Y,P51X, , , ,0,0
FLST,3,1,5,ORDE,1
FITEM,3,7
ARSYM,X,P51X, , , ,0,0
CSYS,1
FLST,3,2,5,ORDE,2
FITEM,3,7
FITEM,3,26
AGEN,4,P51X, , , ,90, , ,0
NUMMRG,ALL, , , ,LOW
NUMCMP,ALL
!*
EDPART,CREATE
NSEL,S, , , 41764
NPLOT
CM,jiedian,NODE
ALLSEL,ALL
*SET,calu_t,3.0e-4 !计算时间
*SET,t_c,1.0e-4 !脉冲力宽度时间
*SET,t_inc,0.5e-6 !计算时间步长
*SET,j,0
*SET,pi,3.1415926
*SET,m_hang,t_c/t_inc+2
*DIM,time,ARRAY,m_hang,1,1, , ,
*DIM,force,ARRAY,m_hang,1,1, , ,
*do,i,0,t_c,t_inc
*SET,j,j+1
*SET,time(j,1,1),i !时间
!!!!!原函数force(j,1,1)=-sin(pi*i/t_c)
!!!!!修改后函数,其中freq为周期，Num为周期，需要自行定义
*SET,freq,100000
*SET,Num,10
*SET,force(j,1,1),0.5*cos(2*pi*freq*i)*(1-sin(2*pi*freq*i/Num))
*enddo
*SET,time(m_hang,1,1),calu_t ! 计算时间
EDLOAD,ADD,Fy,0,jiedian,TIME,FORCE, 0, , , , ,
EPLOT
EDFPLOT, 1
ALLSEL,ALL
FINISH
/SOL
EDENERGY,1,1,1,1
EDBVIS,1.5,0.06
EDHGLS,0.1
TIME,calu_t !时间
EDCTS,0,0.9
EDOPT,ADD,blank,LSDYNA !输出lsdyna
!EDOPT,ADD,blank,ANSYS !输出
!EDOPT,ADD,blank,BOTH !输出全部
EDRST,800,
EDHTIME,1000,
EDDUMP,1,
EDWRITE,LSDYNA,'mg3k','k',' ' !输出lsdyna
!EDWRITE,ANSYS,'mg3k','k',' ' !输出
!EDWRITE,BOTH,'mg3k','k',' ' !输出全部
SAVE