*HEADING
 DISK ROLLING ON DRUM.
 C3D8H elements; Hyperelastic material; R_i = 1.0;
 Delta = 0.2; Friction = 0.2;
 Uses 2D dependency file  pstca4shhfa
*PREPRINT,MODEl=YES
*NODE,NSET=NDRUM
*SYMMETRIC MODEL GENERATION,REVOLVE,NODEOFF=500,ELEMENT=500
*RIGID BODY,ANALYTICAL SURFACE=SDRUM,REFNODE=NDRUM
*SURFACE,NAME=SDRUM,TYPE=REV
 START, 1000.0,  2.0
 LINE,  1000.0, -2.0
*SURFACE,NAME=STREAD
 ETREAD, S4
*CONTACT PAIR, INTERACTION=TREADDRUM
 STREAD, SDRUM
*SURFACE INTERACTION,NAME=TREADDRUM
*FRICTION
 0.0,
*****************************
*STEP,INC=100,NLGEOM
 FOOTPRINT
*STATIC
*BOUNDARY
 NSIDE1, 3
 NSIDE2, 3
 NRIM, 1, 3
 NDRUM, 1
 NDRUM, 2, , 0.1
 NDRUM, 3, 6
*NODE PRINT,NSET=NDRUM,FREQ=1
 U2, RF1, RF2, RM3
*CONTACT PRINT,FREQ=100
*EL PRINT,FREQ=100,ELSET=ETREAD
 S11, S22, S33, CE11, CE22, CE33
*OUTPUT, FIELD,FREQ=100
*NODE OUTPUT
 U
*ELEMENT OUTPUT, ELSET=ETREAD
 S
*CONTACT OUTPUT,MASTER=RIM,VARIABLE=ALL
*OUTPUT,HIST,FREQ=1
*NODE OUTPUT,NSET=NDRUM
 U, RF
*END STEP
*****************************
*STEP,INC=100,NLGEOM,UNSYMM=YES
 STEADY STATE ROLLING
*STEADY STATE TRANSPORT
*CHANGE FRICTION,INTERACTIO=TREADDRUM
*FRICTION,SLIPTOL=0.0025
 0.02,
*TRANSPORT VELOCITY
 NTIRE, 0.946345
*TRANSPORT VELOCITY
 NDRUM, -0.002
*EL FILE,FREQ=100,ELSET=ETREAD
 S, E
*NODE FILE,NSET=NTREAD,FREQ=1
 V,
*NODE FILE,NSET=NDRUM,FREQ=1
 U, RF
*END STEP