目前在练习冷喷涂颗粒碰撞的例子,看人家文献里用的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 |