*heading
Oblique impact of a sphere with a rigid plane at a constant normal velocity but varying tangential velocities
**
*preprint, model=yes, history=yes
**
*parameter
tangential_speed=20
**
**NODES
**
*node,nset=qa_test
1, -0.000025, 0.0, 0.0
2, 0.000025, 0.0, 0.0
3, 0.00005, 0.00005, -0.00001001
4, -0.00005, 0.00005, -0.00001001
5, -0.000050, -0.00005, -0.00001001
6, 0.00005, -0.00005,-0.00001001
7, 0.0, 0.0, 0.0
**
**NODE SETS
**
*nset, nset=Steel
1,
*nset, nset=Polyethylene
2,
*nset ,nset=dem
1, 2
*nset,nset=refnode
7,
**
**ELEMENTS and ELEMENT SETS
**
*element, type=pd3d, elset=Steel
1, 1
*element, type=pd3d, elset=Polyethylene
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=Steel,shape=sphere,density=7850
1e-5
*discrete elasticity
208e9, 0.3
*discrete section, elset=Polyethylene,shape=sphere,density=1400
1e-5
*discrete elasticity
1.0e9, 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=Steel_surf
Steel,
*surface, type=element, name=Polyethylene_surf
Polyethylene,
*surface, type=element, name=rigid_plane_surf
rigid_plane, SPOS
**
**SURFACE INTERACTION 
**
*surface interaction, name=Steel_int
*friction
0.3
*surface behavior, pressure-overclosure=hertz
99999999999.9
*contact damping, definition=damping coefficient
0.0
*surface interaction, name=Polyethylene_int
*friction
0.3
*surface behavior, pressure-overclosure=hertz
99999999999.9
*contact damping, definition=damping coefficient
0.0
**
**INITIAL CONDITIONS
**
*initial conditions, type=velocity
1, 3, -5
2, 3, -5
1, 2, <tangential_speed>
2, 2, <tangential_speed>
**
** MATERIAL
**
*material,name=mat-1
*elastic
2e11,0.3
*density
8000
**
**BOUNDARY CONDITIONS
**
*boundary
7, 1, 6
**
**STEP
**
*step
*dynamic, explicit, element by element
, 1.e-7
**
**CONTACT DEFINITIONS
**
*contact
*contact inclusions
Steel_surf, rigid_plane_surf
Polyethylene_surf, rigid_plane_surf
*contact property assignment
Steel_surf, rigid_plane_surf, Steel_int
Polyethylene_surf, rigid_plane_surf, Polyethylene_int
*contact controls assignment, rotational terms=structural
**
**HISTORY OUTPUT
**
*output, history, freq=1
*node output, nset=qa_test
V, VR
*node output, nset=Steel
V, VR
*node output, nset=Polyethylene
V, VR
*energy output, variable=all
*INCREMENTATION OUTPUT,var=all
**
**FIELD OUTPUT
**
*output, field, timeinterval=5e-8
*node output, nset=dem
U,V,A
*output, field, timeinterval=5e-8
*contact output
CTHICK, CFORCE, CSTRESS
*end step