Volume Fraction B

This is an Eulerian mesh and a rigid plate. The keyword *INITIAL_VOLUME_FRACTION_GEOMETRY is used to define a metal cylinder inside the mesh. The cylinder is assigned an initial velocity and impacts the rigid plate.

*ALE_MULTI-MATERIAL_GROUP
*CONSTRAINED_LAGRANGE_IN_SOLID
*CONTROL_ALE
*CONTROL_TERMINATION
*CONTROL_TIMESTEP
*DATABASE_BINARY_D3PLOT
*DEFINE_CURVE
*ELEMENT_SHELL
*ELEMENT_SOLID
*END
*INITIAL_VELOCITY_GENERATION
*INITIAL_VOLUME_FRACTION_GEOMETRY
*KEYWORD
*MAT_PLASTIC_KINEMATIC
*MAT_VACUUM
*NODE
*PART
*SECTION_SHELL
*SECTION_SOLID

  *KEYWORD
  *CONTROL_TERMINATION
  1.0e-3
  *CONTROL_TIMESTEP
  0.0, 0.8
  *DATABASE_BINARY_D3PLOT
  1.0e-4
  *CONTROL_ALE
  0, 1, 2, -1.0
  
  *ALE_MULTI-MATERIAL_GROUP
  1, 1
  2, 1
  *MAT_VACUUM
  1, 1.0
  *MAT_PLASTIC_KINEMATIC
  2, 8000.0, 200.0e9, 0.3, 200.0e6
  
  *MAT_PLASTIC_KINEMATIC
  3, 8000.0, 200.0e9, 0.3, 200.0e6
  
  *SECTION_SOLID
  1, 11
  *SECTION_SOLID
  2, 11
  *SECTION_SHELL
  3, 0
  0.1, 0.1, 0.1, 0.1
  *PART
  
  1, 1, 1 
  *PART
  
  2, 2, 2 
  *PART
  
  3, 3, 3 
  *INITIAL_VOLUME_FRACTION_GEOMETRY
  1, 1, 1
  4, 0, 2
  0.03, 0.0, 0.0, 0.18, 0.0, 0.0,  0.1, 0.1
  *INITIAL_VELOCITY_GENERATION
  1, 2, 0.0, 200.0, 0.0, 0.0
  
  *CONSTRAINED_LAGRANGE_IN_SOLID
  3, 1, 1, 1, 4, 4, 2, 1
  0, 0, -1
  
  *DEFINE_CURVE
  1 
  0.0, 0.0
  1.0e-3, 1.0e9
  *NODE
  $
  *ELEMENT_SOLID
  $
  *ELEMENT_SHELL
  $
  *END