LS-DYNA钢球撞击刚性平面的命令流

/FILNAME,PRACTICE,ON
/TITLE,DROP OF A HOLLOW BALL
/PREP7
ET,1,SHELL163
KEYOPT,1,1,2
R,1,1, ,2.5E-3, , , ,

!建立一个完全固定的刚性材料模型
EDMP,RIGI,1,7,7
MP,DENS,1,7800
MP,EX,1,2E11
MP,NUXY,1,0.3
!建立2号材料为双线性弹塑性模型
MP,DENS,2,2600
MP,EX,2,7E10
MP,NUXY,2,0.3
TB,BKIN,2, , , ,
TBDAT,1,3.5E9
TBDAT,2,1.35E10

!直接建立刚性平面的有限元模型
N, ,-0.6,0.0,-0.6, , , ,
N, ,-0.6,0.0,0.6, , , ,
N, ,0.6,0.0,0.6, , , ,
N, ,0.6,0.0,-0.6, , , ,
E,1,2,3,4

!建立空心球壳的几何模型并进行网格划分
SPH4,0.0,0.3,0.06
VDELE,1
AATT,2,1,1
AESIZE,ALL,0.03
AMESH,ALL

!在球壳表面和刚性平台之间建立面对面自动接触
NSEL,S,LOC,Y,0.1,0.5
CM,SPHERE,NODE
ALLSEL,ALL
NSEL,S,LOC,Y,-0.1,0.1
CM,TABLE,NODE
ALLSEL,ALL
EDCGEN,ASTS,SPHERE,TABLE,0,0,0,0,0, , , , ,0,10000000

!施加初始速度
EDVE,VELO,SPHERE,0,-5.05,0,0,0,0, , , , , ,
!建立时间以及加速度数组并赋值
*DIM,TIME,ARRAY,2,1,1, , ,
*DIM,GY,ARRAY,2,1,1, , ,
*SET,TIME(1,1,1),0.0
*SET,TIME(2,1,1),3.0
*SET,GY(1,1,1),9.8
*SET,GY(2,1,1),9.8
!用EDLOAD命令来施加ACLY载荷
EDLOAD,ADD,ACLY,0,SPHERE,TIME,GY,0, , , , ,

FINISH
/SOL
TIME,2.0,                                !设置求解的结束时间
EDRST,30,                                !设置结果文件Jobname.RST的输出步骤
EDHTIME,30,                              !设置结果文件Jobname.HIS的输出步骤
EDDUMP,1,                                !设置重启文件d3dump的输出步骤
EDENERGY,1,1,1,1                         !设置能量控制选项
SOLVE

Author: Taochun Yang

佐雍得尝，与众分享，收获快乐！