*HEADING
 CORNERING ROLLING FRICTION. 
 C3D8I elements; Elastic material; R_i = 0.2;
 Delta = 0.05; Friction = 0.02;   
 Uses 2D dependency file  PSTCA4SIEFA.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, 7
  45.0, 20
 157.5, 7
*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
*NODE PRINT,NSET=NROAD,FREQ=1
 U2, RF1, RF2, RM3
*EL PRINT,FREQ=0
*EL FILE,FREQ=0
*END STEP
*****************************
*STEP,INC=100,NLGEOM,UNSYMM=YES
 STEADY STATE ROLLING
*STEADY STATE TRANSPORT
*CHANGE FRICTION,INTERACTIO=SRIGID
*FRICTION,SLIPTOL=0.0025
 0.02, 
*MOTION,TYPE=VELOCITY,ROTATION
 NTIRE, -0.02,   0.0, 0.0, 100.0,  0.0, 1.0, 100.0
*TRANSPORT VELOCITY
 NTIRE, -1.00673
*EL FILE,FREQ=100,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,
*NODE FILE,NSET=NROAD,FREQ=1
 U, RF
*END STEP