*HEADING
 STRAIGHT LINE ROLLING FRICTION. 
 C3D20R elements; Elastic material; R_i = 1.7;
 Delta = 0.05; Friction = 0.02;   
 Uses 2D dependency file  PSTCA8SREFA.RES 
*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
 0.0, 1.0, 0.0
 157.5, 4
  45.0, 10
 157.5, 4
*rigid body,analytical surface=SROAD,REFNODE=1
*surface,name=SROAD,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
*surface,NAME=STREAD
 ETREAD, S4
*CONTACT PAIR, INTERACTION=SRIGID
 STREAD, SROAD
*SURFACE INTERACTION,NAME=SRIGID
*FRICTION
 0.0, 
*****************************
*STEP,INC=100,NLGEOM
 FOOTPRINT
*STATIC
*BOUNDARY
 NSIDE1, 3
 NSIDE2, 3
 NRIM,  1, 3
 NROAD, 1
 NROAD, 2, , 0.05
 NROAD, 3, 6
*CONTACT PRINT,FREQ=100
*CONTACT FILE,FREQ=100
*OUTPUT,FIELD,FREQ=100
*CONTACT OUTPUT,VARIABLE=PRESELECT
*NODE PRINT,NSET=NROAD,FREQ=1
 U2, RF1, RF2, RM3
*EL PRINT,FREQ=0
*EL FILE,FREQ=0
*OUTPUT,FIELD,FREQ=0
*ELEMENT OUTPUT
*END STEP
*****************************
*STEP,INC=100,NLGEOM,UNSYMM=YES
 STEADY STATE ROLLING
*STEADY STATE TRANSPORT
*CHANGE FRICTION,INTERACTIO=SRIGID
*FRICTION,ELAS=0.01
 0.02, 
*MOTION,TYPE=VELOCITY,TRANSLATION
 NTIRE, 1, 1, 2.0
*TRANSPORT VELOCITY
 NTIRE, -1.002834
*EL FILE,FREQ=100,ELSET=ETREAD
 S, E
*OUTPUT,FIELD,FREQ=100
*ELEMENT OUTPUT,ELSET=ETREAD
S,E
*NODE PRINT,NSET=NTREAD,FREQ=100
 V1, V2, V3
*NODE PRINT,NSET=NROAD,FREQ=100
 U2, RF1, RF2, RM3
*NODE FILE,NSET=NTREAD,FREQ=1
 V,
*OUTPUT,FIELD,FREQ=1
*NODE OUTPUT,NSET=NTREAD
V,
*NODE FILE,NSET=NROAD,FREQ=1
 U, RF
*OUTPUT,FIELD,FREQ=1
*NODE OUTPUT,NSET=NROAD
U,RF
*END STEP