*HEADING
 STRAIGHT LINE ROLLING WITH MATERIAL HISTORY EFFECTS INCLUDED. 
 C3D8RH elements
 USES DEPENDENCY FILE: spinningdisk_axi
*RESTART,WRITE,FREQ=100
*NODE,NSET=NROAD
 1, 0.0, 0.0, 0.5
*SYMMETRIC MODEL GENERATION,REVOLVE,NODEOFF=500,ELEMENT=500
 0.0, 0.0, 0.0,    0.0, 0.0, 1.0
 1.0, 0.0, 0.0
 360.0,90
*rigid body,analytical surface=SURFACE_ROAD,REFNODE=1
*surface,name=SURFACE_ROAD,TYPE=CYLINDER
 0.0, 0.0, 0.0,   1.0, 0.0, 0.0
 0.0, 0.0, -1.0
 START, -2.0, -2.0
 LINE,   2.0, -2.0
*COUPLING,CONSTRAINT=XX,REFNODE=AXLE,SURFACE=RIM
*KINEMATIC
 1,6
*surface,NAME=SURFACE_TREAD
 ETREAD, S4
*CONTACT PAIR, INTERACTION=SURFACE_RIGID
 SURFACE_TREAD, SURFACE_ROAD
*SURFACE INTERACTION,NAME=SURFACE_RIGID
*FRICTION
 0.0, 
*ELSET,ELSET=SECT,GEN
 1,499
*ELSET,ELSET=EOUT
 SECT,ETREAD
*NSET,NSET=SECT,GEN
 1, 499
*NSET,NSET=NOUT
 AXLE,NROAD
*****************************
*STEP,INC=100,NLGEOM
 FOOTPRINT
*STATIC, LONG TERM
 0.25, 1.0
*BOUNDARY
 NSIDE1, 3
 NSIDE2, 3
 AXLE,  1, 6
 NROAD, 1
 NROAD, 2, , 0.3
 NROAD, 3, 6
*EL PRINT,FREQ=100,ELSET=ETREAD
 S11, S22, S33, CE11, CE22, CE33
*NODE PRINT,NSET=NOUT,FREQ=1
 RF
*OUTPUT,FIELD,FREQ=10
*NODE OUTPUT
 U
*ELEMENT OUTPUT
 S
*NSET,NSET=QA_TEST_NOUT
 NOUT
*OUTPUT,HIST
*NODE OUTPUT,NSET=NOUT
 RF
*NODE OUTPUT,NSET=QA_TEST_NOUT
 RF1,RF2,RF3
*END STEP
*****************************
*STEP,INC=1000,NLGEOM
 LAGRANGIAN ROLLING
*VISCO,CETOL=0.001
**Set total step time to 2.5 to get 2 revolutions.
 0.2E-4, 0.5
*BOUNDARY,OP=NEW
 NSIDE1, 3
 NSIDE2, 3
 AXLE,  1, 5
 NROAD, 1
 NROAD, 2, , 0.3
 NROAD, 3, 6
*BOUNDARY,OP=NEW,TYPE=VELOCITY
 AXLE, 6, ,2.5
*END STEP