*IMPORT,UPDATE=NO,STATE=YES,STEP=1
 FOAM, PUNCH
*SURFACE,TYPE=ELEMENT,NAME=PUNCH
 PUNCH, SPOS
*SURFACE,TYPE=ELEMENT,NAME=FOAM
 TOP, S3
*CONTACT PAIR,INTERACTION=INTPROP
 FOAM, PUNCH
*SURFACE INTERACTION,NAME=INTPROP
*FRICTION,TAUMAX=1.15E5
 0.2,
** 
** BOUNDARY CONDITIONS
*BOUNDARY
 LFT, XSYMM
 BOT, ENCASTRE
 PUNCHREF, 1
 PUNCHREF, 6
*AMPLITUDE,NAME=DISP,TIME=TOTAL TIME,DEFINITION=SMOOTH STEP
 0., 0., 0.05, 0.027
** 
** STEP: Step-1
*STEP,INC=500,NLGEOM,UNSYMM=YES
 Step-1: Initialization - ramp down initial out-of-balance forces
*STATIC
 1.0e-6, 1.0e-6
**
** APPLIED LOAD
*BOUNDARY,AMPLITUDE=DISP
 PUNCHREF, 2, 2, -1.0
**BOUNDARY
** PUNCHREF, 2, 2, -0.02
**
** Relax tolerance from a default of 0.005 to 0.01 because the time and
** space averaging produce values that are very small, thus requiring
** convergence with unnecessarily small residuals. Relax displacement
** tolerance from 0.01 to 0.1
*CONTROLS,PARAMETERS=FIELD
 0.01, 0.1
, 0.0, 
**
*MONITOR,NODE=PUNCHREF,DOF=2
*ENERGY PRINT,FREQUENCY=999
*NODE FILE,NSET=PUNCHREF
 U, RF
*OUTPUT,FIELD,FREQ=100
*ELEMENT OUTPUT
 S, PE, PEEQ, LE
*NODE OUTPUT
 U,
*OUTPUT,HISTORY,FRE=1
*ELEMENT OUTPUT,ELSET=OUT
 MISES, PRESS, PEEQ
*NODE OUTPUT,NSET=PUNCHREF
 U, RF1, RF2
*END STEP
** STEP: Step-2
*STEP,INC=500,NLGEOM,UNSYMM=YES
 Step-2
*STATIC
 0.005, 0.025
**
** APPLIED LOAD
*BOUNDARY,AMPLITUDE=DISP
 PUNCHREF, 2, 2, -1.0
**-----
*OUTPUT,FIELD,FREQ=100
*ELEMENT OUTPUT
 S, PE, PEEQ, LE
*NODE OUTPUT
 U,
*OUTPUT,HISTORY,FRE=1
*ELEMENT OUTPUT,ELSET=OUT
 MISES, PRESS, PEEQ
*NODE OUTPUT,NSET=PUNCHREF
 U, RF
*OUTPUT,FIELD, time interval=0.025
*ELEMENT OUTPUT, elset=qa_test
 MISES,
*OUTPUT,HISTORY, time interval=1.25e-3
*ELEMENT OUTPUT,ELSET=qa_test
 MISES, PRESS,
*NODE OUTPUT,NSET=qa_test
 U, RF1, RF2
**----
*END STEP