*heading
Impact of two different spheres with a constant normal velocity but varying angular velocities
*preprint, model=yes, history=yes
**
*parameter
angular_speed=20
**
**NODES
**
*node,nset=qa_test
1, -0.1000005, 0.0, -0.8
2, 0.5000005, 0.0, -0.8
3, -0.1000005, 0.0, 0.8
4, 0.5000005, 0.0, 0.8
*nset, nset=Al_alloy
1, 2
*nset, nset=Nylon
3, 4
*nset, nset=dem
1, 2, 3, 4
**
**ELEMENTS and ELEMENT SETS
**
*element, type=pd3d, elset=Al_alloy1
1, 1
*element, type=pd3d, elset=Al_alloy2
2, 2
*element, type=pd3d, elset=Nylon1
3, 3
*element, type=pd3d, elset=Nylon2
4, 4
*elset, elset=dem
1, 2, 3, 4
*elset, elset=Al_Alloy
1, 2
*elset, elset=Nylon
3, 4
**
**SECTIONS 
**
*discrete section, elset=Al_alloy1,shape=sphere,density=2700
0.1
*discrete elasticity
70.0e9, 0.33
*discrete section, elset=Al_alloy2,shape=sphere,density=2.7E6
0.5
*discrete elasticity
70.0e9, 0.33
*discrete section, elset=Nylon1,shape=sphere,density=1000
0.1
*discrete elasticity
2.5e9, 0.4
*discrete section, elset=Nylon2,shape=sphere,density=1e6
0.5
*discrete elasticity
2.5e9, 0.4
*section control, name=autodem




**<fStability>,<fTracking>,<fCollision>,<fRayleigh>,<fSepDist>,<fRotation>
 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=Nylon_surf
Nylon,
**
**SURFACE INTERACTION 
**
*surface interaction, name=Al_alloy_int
*friction
0.4
*surface behavior, pressure-overclosure=Hertz
1e10
*contact damping, definition=damping coefficient
0.0
*surface interaction, name=Nylon_int
*friction
0.4
*surface behavior, pressure-overclosure=Hertz
1e10
*contact damping, definition=damping coefficient
0.0
**
**INITIAL CONDITIONS
**
*initial conditions, type=velocity
1, 1, 0.2
1, 5, <angular_speed>
3, 1, 0.2
3, 5, <angular_speed>
**
**STEP
**
*step
*dynamic, explicit, element by element
, 3e-3
**
**CONTACT DEFINITIONS
**
*contact
*contact inclusions
Al_alloy_surf, Al_alloy_surf
Nylon_surf, Nylon_surf
*contact property assignment
Al_alloy_surf, Al_alloy_surf, Al_alloy_int
Nylon_surf, Nylon_surf, Nylon_int
**
**HISTORY OUTPUT
**
*output, history, freq=1
*node output, nset=Al_alloy
V, VR
*node output, nset=Nylon
V, VR
*energy output, variable=all
*INCREMENTATION OUTPUT,var=all
**
**FIELD OUTPUT
**
*output, field, timeinterval=1e-3
*node output, nset=qa_test
V
*node output, nset=dem
U,V,A
*contact output
CFORCE
*end step