求助刚性球撞板

greenhh

想知道刚性球撞板，板的情况。
不知道哪里出问题了。麻烦大家帮我看看。

命令流如下，
!定义单元类型
/PREP7  

ET,1,SOLID168         !小球
  
ET,2,SOLID168         !链板
  
EDMP,RIGI,1,7,7      !小球刚性体
MP,DENS,1,2300      
MP,EX,1,4.8e10  
MP,NUXY,1,0.2
  
MP,DENS,2,7850       !链板弹性体
MP,EX,2,2.1e11  
MP,NUXY,2,0.3   


BLOCK,-1.450,1.450,0,-0.02,-0.1,0.1,     
K,9,-0.1,-0.02,-0.1,
K,10,-0.1,-0.02,0.1,
K,11,0.1,-0.02,0.1,  
K,12,0.1,-0.02,-0.1,
K,13,0.5,-0.02,-0.1,
K,14,0.5,-0.02,0.1,  
K,15,-0.5,-0.02,0.1,
K,16,-0.5,-0.02,-0.1,
K,17,-0.7,-0.02,-0.1,
K,18,-0.7,-0.02,0.1,
K,19,0.7,-0.02,0.1,  
K,20,0.7,-0.02,-0.1,
K,21,1.1,-0.02,-0.1,
K,22,1.1,-0.02,0.1,  
K,23,-1.1,-0.02,0.1,
K,24,-1.1,-0.02,-0.1,
K,25,-1.3,-0.02,-0.1,
K,26,-1.3,-0.02,0.1,
K,27,1.3,-0.02,0.1,  
K,28,1.3,-0.02,-0.1,
LSTR,25,26
LSTR,24,23
LSTR,17,18
LSTR,16,15
LSTR,9,10
LSTR,12,11
LSTR,13,14
LSTR,20,19
LSTR,21,22
LSTR,28,27
LPLOT   
FLST,3,10,4,ORDE,2  
FITEM,3,13  
FITEM,3,-22
ASBL,       3,P51X  
VPLOT
SPH4,0,1,0.175

VSEL,S,VOLU,,2       !划分网格
VATT,1,,1,0,
VSEL,S,VOLU,,1
VATT,2,,2,0,
ALLSEL,ALL
MSHKEY,0                 !采用自由网格化分
MSHAPE,1,3D
LESIZE,ALL,0.05,
VMESH,ALL

EDPART,CREATE          !定义PART号,球2，板1
EDCGEN,ASTS,1,2,0,0,0,0,0,,,,,0,10000000


DA,9,UY,0                     !加载约束
DA,11,UY,0
DA,13,UY,0
DA,15,UY,0
DA,17,UY,0

EDVE,VELO,2,0, -11.7 ,0,0,0,0, , , , , ,    !加初速度


*DIM,TIME,ARRAY,2,1,1, , ,                     !定义加速度数组
!*  
*SET,TIME(2,1,1) , 1
*DIM,GY,ARRAY,2,1,1, , ,  
!*  
*SET,GY(1,1,1) , 9.8  
*SET,GY(2,1,1) , 9.8  
EDLOAD,ADD,ACLY,0,2,TIME,GY, 0, , , , ,
FINISH  

/SOL
TIME,0.1,   
EDOPT,ADD,blank,BOTH
EDRST,30,   
EDHTIME,30,
EDDUMP,1,
SAVE
SOLVE