在LS-DYNA的SHELL163叠层问题请教

主要现在要做一个有叠层角度90/90/0/0] S的叠层，

冲击模型试片长宽高为200毫米* 50MM * 1.2MM，做一种钢球自由落下

冲击分析，使用元素是SOLID164跟SHELL163

在SHELL163不知道如何设定叠层角度及设定方法

以下是我所写的程式：

FINISH

/CLEAR,START

/CONTOUR,1,128

/eshape,1

!*****DefineParameters*******

pi =ACOS(-1) !pi=3.14

!***********************************************

Leng=0.2 ! 冲击试片尺寸

Wide=0.05

High=0.0012

Th=0.00015

BBB=0.138 !冲击高度

ball.03808 ! 钢球直径

!**************************************

/PREP7

!**************************************

!***** ELEMENTTYPE *****

ET,1,shell163

ET,2,solid164

!*******TC12K35/EPO-AIMP材料性质********

R,1,5/6,,high

MPTEMP,1,0 !AL

MP,EX,1,129.86e9

MP,EY,1,7.71e9

MP,EZ,1,7.71e9

MP,PRXY,1,0.33

MP,PRYZ,1,0.33

MP,PRXZ,1,0.33

MP,GXY,1,3.45e9

MP,GYZ,1,3.45e9

MP,GXZ,1,3.45e9

MP,DENS,1,1700

!*************** 钢球材料性值******************

MP,EX,2,200e9

MP,PRXY,2,0.3

MP,DENS,2,7847

!******************* 板建模*********************

k,1,0,0,0

k,2,leng,0,0

k,3,leng,0,-Wide

k,4,0,0,-wide

a,1,2,3,4

TYPE,1

MAT,1

REAL,1

LSEL,S,LOC,x,0

LSEL,A,LOC,x,leng

LEsize,all,,,10

ALLSEL

Amesh,all,all

LOCAL,15,0,,,,,-90

wpcsys,15

!*******************球建模*********************

WPOFFS,Leng/2,wide/2,BBB+ballsize/2

SPH4,0,0,ballsize/2

WPROTA,,90,

VSEL,S,,,1

type,2

MAT,2

ESIZE,,5

ASEL,ALL

MSHAPE,1,3D

MSHKEY,0

VMESH,1

wpcsys,15

WPROTA,,90,

ALLSEL

CM,allnode,NODE

ASEL,S,,,1

NSLV,S,1

CM,BEAM,NODE

ALLSEL,ALL

ASEL,S,,,2,3

NSLV,S,1

CM,CUBE,NODE

ALLSEL,ALL

!*** Contact***

EDCGEN,ASTS,CUBE,BEAM,,,,,,,,,,0,1E7

FINISH

!**************************************

/SOLU

!**************************************

ALLSEL

!***BOUND ***

NSEL,all

NSEL,S,LOC,x,Leng

NSEL,A,LOC,x,0

D,ALL,ALL

NSEL,all

ALLSEL

!*** InitialVelocity ***

EDVEL,VELO,CUBE,0,0,0,0,0 !EDVEL or EDPVEL

!*** 时间及加速度之定义 ***

*DIM,TIME,ARRAY,2,1,1

*DIM,GY,ARRAY,2,1,1

*SET,TIME(1,1,1),0.0 !*time_start

*SET,TIME(2,1,1),0.7 !*time_end

*SET,GY(1,1,1),9.81

*SET,GY(2,1,1),9.81

!*** Specifiesloads for an explicit dynamics analysis. ***

EDLOAD,ADD,ACLY,0,allnode,TIME,GY,0

!*** LS-DynaSET ***

TIME,0.7

EDRST,100

EDHTIME,100

!EDHIST,ALLNODE

EDDUMP,1

EDENERGY,1,1,1,1

EDOPT,ADD,BLANK,BOTH

EDOUT,ALL

EDWRITE,BOTH,'FILE','K',''

SAVE

EDRUN,SMP,1,2

SOLVE

FINISH

n_number=NODE(leng/2,Wide/2,0)

!*** POSTSetting ***

/POST1

/DSCALE,1,1.0

SET,LAST

PLNSOL,U,Y,0,1.0

/POST26

NSOL,2,81,U,Y,Uy

PLVAR,2,

*GET,UYNODE,NODE,n_number,u,y

!/POST1

!/DSCALE,1,1.0

!SET,LAST

!PLNSOL,S,EQV,0,1

!

!/POST26

!NSOL,2,NODE(1,0.5,0),U,Y,Uy

!NSOL,3,NODE(1.5,0,0),U,Y,Uy

!DERIV,4,2,1,,Vy,,,1

!DERIV,5,3,1,,Ay,,,1

!

!PLVAR,3

希望使用SHELL163的叠层问题能解决，在程式中有什么问

题麻烦帮我解答，感激