*heading
verification test normal displacement controlled contct
two identical spheres (ADHESION: JKR and JKR shifted)
**
*preprint, model=yes, history=yes
**
**NODES
**
*node,nset=qa_test
1, -5.000, 0.0, -5.5
2, 5.00, 0.0, -5.5
3, -5.00, 0.0, 5.0
4, 5.000, 0.0, 5.0
**
**NODE SETS
**
*nset, nset=JKR
1, 2
*nset, nset=JKR_SHIFTED
3, 4
*nset, nset=dem
1, 2
3, 4
**
**ELEMENTS and ELEMENT SETS
**
*element, type=pd3d, elset=JKR
1, 1
2, 2
*element, type=pd3d, elset=JKR_SHIFTED
3, 3
4, 4
*elset, elset=dem
1, 2
3, 4
**
**SECTION DEFINITIONS 
**
*elset,elset=JKR-1
 1
*elset,elset=JKR-2
 2
*elset,elset=JKR_SHIFTED-1
 3
*elset,elset=JKR_SHIFTED-2
 4
*discrete section, elset=JKR-1,shape=sphere,density=1.4e-09
 5.0
*discrete elasticity
 1.0e+03,0.3
*discrete section, elset=JKR-2,shape=sphere,density=1.4e-09
 5.0
*discrete elasticity
 1.0e+03,0.3 
*discrete section, elset=JKR_SHIFTED-1,shape=sphere,density=1.4e-09
 5.0
*discrete elasticity
 1.0e+03,0.3
*discrete section, elset=JKR_SHIFTED-2,shape=sphere,density=1.4e-09
 5.0
*discrete elasticity
 1.0e+03,0.3 
*section control, name = autodem


**-10.0,-10.0,-10.0,10.0,10.0,10.0


**<fStability>,<fTracking>,<fCollision>,<fRayleigh>,<fSepDist>,<fRotation>
 0.1, 0.1, 0.1, 0.1, 0.1, 0.1 
**
**SURFACES
**
*surface, type=element, name=JKR_surf
 JKR,
*surface, type=element, name=JKR_SHIFTED_surf
 JKR_SHIFTED,
**
**SURFACE INTERACTION 
**
*surface interaction, name=JKR_int
*friction
 0.0
*surface behavior, pressure-overclosure=JKR
 5.0e-02, 170.0
*contact damping, definition=damping coefficient
 0.0
*surface interaction, name=JKR_SHIFTED_int
*friction
 0.0
*surface behavior, pressure-overclosure=JKR
 5.0e-02, 170.0,SHIFTED
*contact damping, definition=damping coefficient
 0.0
**
** displacement controlled test
**
*boundary
1,2,6
2,1,6
3,2,6
4,1,6
*amplitude, name=disp,definition=smooth
0.0,0.0,2.5e-03,0.01,10.0e-03,-0.01
**
**STEP
**
*step
*dynamic, explicit, , element by element
**2.5e-8, 10.0e-03
 , 10.0e-03
**
*boundary, amp=disp
1,1,1,1.0
3,1,1,1.0
**
**CONTACT DEFINITIONS
**
*contact
*contact inclusions
JKR_surf, JKR_surf
JKR_SHIFTED_surf, JKR_SHIFTED_surf
*contact property assignment
JKR_surf, JKR_surf, JKR_int
JKR_SHIFTED_surf, JKR_SHIFTED_surf, JKR_SHIFTED_int
**
**HISTORY OUTPUT
**
*output, history, freq=1
*node output, nset=JKR
U1,RF1
*node output, nset=JKR_SHIFTED
U1,RF1
*incrementation output,var=all
**
**FIELD OUTPUT
**
*output, field, timeinterval=5e-3
*node output, nset=dem
U,V,A
*node output, nset=qa_test
 rf
*contact output
CFORCE
*end step