2维轴对称运行出错。

haolpk37

有做过2维轴对称准静态分析没，麻烦帮我看一下，我的程序错在哪，提升错误是part与边界不匹配，接触那里也有问题。谢谢解答！

问题是一个刚体穿透三层材料。
下面是我的命令流

fini
/clear
/prep7
!element type
et,1,plane162
KEYOPT,1,3,1
KEYOPT,1,5,1
KEYOPT,1,2,0
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!material for head
EDMP,RIGI,1,6,7
mp,dens,1,7.8e-6
mp,ex,1,2e5
mp,nuxy,1,0.2

!material for bone
MPTEMP,,,,,,,   
MP,DENS,2,3.17e-6   
MP,EX,2,1200
MP,NUXY,2,.2
TB,PLAW,2,,,1,  
TBDAT,1,200
TBDAT,2,0   
TBDAT,3,0   
TBDAT,4,0   
TBDAT,5,0   
TBDAT,6,.2  
TBDAT,7,

!materialfor collage
MP,DENS,3,1.3e-6
MP,EX,3,448
MP,NUXY,3,.49   
TB,PLAW,3,,,1,  
TBDAT,1,65  
TBDAT,2,0   
TBDAT,3,00  
TBDAT,4,0   
TBDAT,5,0   
TBDAT,6,0.3
TBDAT,7,

!materials for sub
MP,DENS,4,1.33e-9   
MP,NUXY,4,.49
TB,MOONEY,4,,,0,
TBDAT,1,0.000001   
TBDAT,2,.3  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! modeling for head
k,1,0,0
k,2,0.0215,0.0125
k,3,0.075,0.2152
k,4,0,0.2152
k,1000,0,0.025
larc,1,2,1000,0.025
l,2,3
l,3,4
l,4,1
al,1,2,3,4

!modeling for bone and collage
k,5,0,0
k,6,4.5,0
k,7,0,-0.1
k,8,4.5,-0.1
k,9,0,-0.2
k,10,4.5,-0.2
l,5,6
l,6,8
l,8,7
l,7,5
al,5,6,7,8
l,7,9
l,9,10
l,10,8
al,10,11,7,9
aglue,2,3
!modeling for sub
k,11,0,-0.2
k,12,6,-0.2
k,13,0,-4.7
k,14,6,-4.7
l,11,12
l,12,14
l,14,13
l,13,11

al,12,13,14,15

allsel,all

!meshing
lesize,5,,,160,1
lesize,7,,,160,1/1
lesize,10,,,160,1

lesize,6,,,10
lesize,8,,,10
lesize,9,,,10
lesize,11,,,10

TYPE,1        ! meshing for bone
MAT,2
MSHAPE,0,2D
MSHKEY,1
AMESH,2
   
MAT,3        ! meshing for collage
MSHAPE,0,2D
MSHKEY,1
AMESH,3


MAT,4       ! meshing for sub
lesize,15,,,90,1/2
lesize,13,,,90,2
lesize,12,,,120
lesize,14,,,120
MSHAPE,0,2D
MSHKEY,1
AMESH,4

MAT,1
AESIZE,1,0.01   
CHKMSH,'AREA'   
AMESH,1
!create part
EDPART,CREATE   

!contacting
EDCGEN,ASTS,1,2,0,0,0,0,0, , , , ,0,10000000,0,0  

!EDCGEN,TDSS,2,3,0,0,0,0,0, , , , ,0,10000000,0,0  

EDCGEN,ASTS,3,4,0,0,0,0,0, , , , ,0,10000000,0,0  


!loading
*DIM,TIME,ARRAY,2,1,1, , ,  
*SET,TIME(2,1,1) , 1e-005   

*DIM,DISP,ARRAY,2,1,1, , ,  
*SET,TIME(2,1,1) , -1

EDLOAD,ADD,RBUY,0,1,TIME,DISP, 0, , , , ,
!
!compents for constraint
NSEL,S,LOC,Y,-4.7,-4.699
!CM,under_y,NODE
D,all,,,,,,,UY
allsel,all

NSEL,S,LOC,X,-0.1,0.0001
NSEL,R,LOC,Y,-4.7,0
!CM,left_x,NODE
D,all,,,,,,UX
allsel,all

/SOL
TIME,1e-5,  
EDCTS,0,0.9,
EDOPT,ADD,blank,ANSYS   
EDRST,100,  
EDHTIME,1000,   
EDDUMP,1,   
EDOUT,RCFORC
EDOUT,DEFGEO
EDOUT,RBDOUT
EDOUT,GCEOUT
!*  
EDWRITE,ANSYS,'fish_scales_ls-dyna','k',' '
FINISH