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

gaoyuan780

主要现在要做一个有叠层角度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                          !冲击高度ballsize=0.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  

TB,PLAW,,,,1

TBDATA,1,270E6                   !*Yield Strength

TBDATA,2,1.35E8                  !*E2

TBDATA,4,6500.0                  !*! C (s-1)

TBDATA,5,4.0                     !*P  

TBDATA,6,0.17                    !*Failure Strain

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

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的叠层问题能解决，在程式中有什么问

题麻烦帮我解答，感激

速度与激情6

发帖辛苦了,楼主 支持