matlab 颗粒碰撞,Lsdyna颗粒碰撞图 - 仿真模拟 - 小木虫 - 学术 科研 互动社区

目前在练习冷喷涂颗粒碰撞的例子，看人家文献里用的shell163单元我也尝试着做了一下，但是感觉做的结果不对不知道错在哪现贴出语句供大家修改，目前我模拟出来如图所示。这个明显是不对的

!*

/PREP7

!*

ET,1,SHELL163

!*

KEYOPT,1,1,2

KEYOPT,1,2,0

KEYOPT,1,3,0

KEYOPT,1,4,0,

!*

*SET,_RC_SET,1,

*SET,_RC_SET,1,

R,1

RMODIF,1,1,0.83333333333333,1,100e-6,100e-6,100e-6,100e-6,

!*

!*

MP,DENS,1,8900

MP,EX,1,117e9

MP,NUXY,1,0.35

TB,PLAW,1,,,1,

TBDAT,1,100e6

TBDAT,2,400e6

TBDAT,3,0.5

TBDAT,4,100

TBDAT,5,100

TBDAT,6,

TBDAT,7,

MP,DENS,2,8900

MP,EX,2,117e9

MP,NUXY,2,0.35

TB,PLAW,2,,,1,

TBDAT,1,100e6

TBDAT,2,400e6

TBDAT,3,0.5

TBDAT,4,100

TBDAT,5,100

TBDAT,6,

TBDAT,7,

PCIRC,10e-6, ,90,-90,

FLST,3,1,5,ORDE,1

FITEM,3,1

AGEN, ,P51X, , , ,10e-6, , , ,1

RECTNG,0,100e-6,0,-100e-6,

CM,_Y,AREA

ASEL, , , ,       1

CM,_Y1,AREA

CMSEL,S,_Y

!*

CMSEL,S,_Y1

AATT,       1,       1,   1,       0,

CMSEL,S,_Y

CMDELE,_Y

CMDELE,_Y1

!*

CM,_Y,AREA

ASEL, , , ,       2

CM,_Y1,AREA

CMSEL,S,_Y

!*

CMSEL,S,_Y1

AATT,       2,       1,   1,       0,

CMSEL,S,_Y

CMDELE,_Y

CMDELE,_Y1

!*

FLST,2,1,5,ORDE,1

FITEM,2,1

AESIZE,P51X,5e-7,

FLST,2,1,5,ORDE,1

FITEM,2,2

AESIZE,P51X,10e-007,

MSHAPE,0,2D

MSHKEY,0

!*

FLST,5,2,5,ORDE,2

FITEM,5,1

FITEM,5,-2

CM,_Y,AREA

ASEL, , , ,P51X

CM,_Y1,AREA

CHKMSH,'AREA'

CMSEL,S,_Y

!*

AMESH,_Y1

!*

CMDELE,_Y

CMDELE,_Y1

CMDELE,_Y2

!*

/UI,MESH,OFF

!*

EDPART,CREATE

!*

EDCGEN,ASS2D, , ,0,0,0,0,0, , , , ,0,10000000,0,0

EDCONTACT,0.9,0,2,0,1,1,1,4,0

EDSP,OFF,' ',' ',1,

!*

FLST,2,3,4,ORDE,3

FITEM,2,2

FITEM,2,-3

FITEM,2,7

!*

/GO

DL,P51X, ,UX,0

FLST,2,1,4,ORDE,1

FITEM,2,4

!*

/GO

DL,P51X, ,ALL,0

EDPV,VELO,       1,0,-600,0,0,0,0,

FINISH

/SOL

EDENERGY,1,1,1,1

EDMP,HGLS,1,4,0.1,1.5,0.06, , ,

EDMP,HGLS,2,4,0.1,1.5,0.06, , ,

TIME,50e-9,

EDOPT,ADD,blank,LSDYNA

EDRST,10,

EDHTIME,10,

EDDUMP,1,

!*

/STATUS,SOLU

SOLVE

FINISH

QQ截图20130804204606.jpg