*heading
normal impact of a sphere with a rigid plane
**
*preprint, model=yes, history=yes
**
**NODES
**
*node
1, -0.25, 0.0, 0.0
2, 0.25, 0.0, 0.0
3, 0.5, 0.5, -0.1000001
4, -0.5, 0.5, -0.1000001
5, -0.5, -0.5, -0.1000001
6, 0.5, -0.5, -0.1000001
7, 0.0, 0.0, 0.0
**
**NODE SETS
**
*Nset, nset=Al_Alloy
1,
*Nset, nset=Mg_Alloy
2,
*Nset,nset=qa_test
1,2
*Nset,nset=dem
1,2
*nset,nset=refnode
7,
**
**ELEMENTS and ELEMENT SETS
**
*element, type=pd3d, elset=Al_alloy
1, 1
*element, type=pd3d, elset=Mg_alloy
2, 2
*element, type=s4r, elset=rigid_plane
3, 3, 4, 5, 6
*Elset, elset=dem
1, 2
**
**SECTIONS
**
*shell section, elset=rigid_plane, material=mat-1, offset=SPOS
0.1, 5
*rigidbody, refnode=refnode, elset=rigid_plane
*discrete section, elset=Al_Alloy,shape=sphere,density=2699
0.1
*discrete elasticity
70e9, .3
*discrete section, elset=Mg_Alloy,shape=sphere,density=1800
0.1
*discrete elasticity
40e9, .35
*section control, name=autodem
 
 
 
 
 0.1, 0.1, 0.1, 0.1, 0.1, 0.1 
**
**SURFACES
**
*surface, type=element, name=Al_alloy_surf
Al_alloy,
*surface, type=element, name=Mg_alloy_surf
Mg_Alloy,
*surface, type=element, name=rigid_plane_surf
rigid_plane, SPOS
**
**SURFACE INTERACTION 
**
*surface interaction, name=Al_alloy_int
*friction
0.0
*surface behavior, pressure-overclosure=hertz
99999999999.9
*contact damping, definition=damping coefficient
0.0
*surface interaction, name=Mg_alloy_int
*friction
0.0
*surface behavior, pressure-overclosure=hertz
99999999999.9
*contact damping, definition=damping coefficient
0.0
**
**INITIAL CONDITIONS
**
*initial conditions, type=velocity
1, 3, -0.2
2, 3, -0.2
**
** MATERIAL
**
*material,name=mat-1
*elastic
2e11,0.3
*density
8000
**
**BOUNDARY CONDITIONS
**
*boundary
7, 1, 6
**
**STEP
**
*step
*dynamic, explicit, element by element
 , 8e-4
**
**CONTACT DEFINITIONS
**
*contact
*contact inclusions
Al_alloy_surf, rigid_plane_surf
Mg_alloy_surf, rigid_plane_surf
*contact property assignment
Al_alloy_surf, rigid_plane_surf, Al_alloy_int
Mg_alloy_surf, rigid_plane_surf, Mg_alloy_int
*contact controls assignment, rotational terms=structural
**
**HISTORY OUTPUT
**
*output, history, freq=1
*node output, nset=qa_test
V
*node output, nset=Al_alloy
V, VR
*node output, nset=Mg_alloy
V, VR
*energy output, variable=all
*INCREMENTATION OUTPUT,var=all
**
**FIELD OUTPUT
**
*output, field, timeinterval=4e-4
*node output, nset=dem
U,V,A
*contact output
CFORCE
*end step