*HEADING
 FRICTIONAL BRAKING OF A ROTATING RIGID BODY IN 3D
**
**  Define node for rigid circular drum.
**
*NODE, NSET=AXIS
 100,  0.000,  0.000,  0.000
*ELEMENT, TYPE=MASS, ELSET=PMASS
 100, 100
*MASS, ELSET=PMASS
 5., 
*ELEMENT, TYPE=ROTARYI, ELSET=RMASS
 200, 100
*ROTARY INERTIA, ELSET=RMASS
 0.001, 0.001, .175
**
**  Start with it spinning at 60 rad/sec.
**
*BOUNDARY
 100, 1, 5
*INITIAL CONDITIONS, TYPE=VELOCITY
 100, 6, 60.
**
**  Define a deformable brake pad.
**
*NODE
 501, -0.200, -0.050, -0.050
 502, -0.250, -0.050, -0.050
 521, -0.200,  0.050, -0.050
 522, -0.250,  0.050, -0.050
 600, -0.250,  0.000,  0.000
*NGEN, NSET=FSIDE
 501, 521, 2
 502, 522, 2
*NCOPY, SHIFT, CHANGE NUMBER=100, OLD SET=FSIDE, NEW SET=NSIDE
 0.000, 0.000, 0.100
 0.000, 0.000, 0.000,  1.000, 0.000, 0.000,  0.000
*NSET, NSET=PAD
 FSIDE, NSIDE
*ELEMENT, TYPE=C3D8R, ELSET=PAD
 501,  501, 503, 504, 502, 601, 603, 604, 602
*ELGEN, ELSET=PAD
 501,  10, 2, 1
*SOLID SECTION, ELSET=PAD, MATERIAL=STUFF
**
**  Make it a soft rubber pad.
**
*MATERIAL, NAME=STUFF
*HYPERELASTIC, N=1
 552.E3, 138.E3, 1.45E-7
*DENSITY
 5.E3, 
**
**  Use a second rigid body to hold the brake pad.
**
*ELEMENT, TYPE=R3D4, ELSET=PLATE
 601,  502, 504, 604, 602
*ELGEN, ELSET=PLATE
 601,  10, 2, 1
*BOUNDARY
 600, 2, 6, 0.
**
**  Run the simulation for 0.2 seconds.
**
*ELSET,ELSET=EALL1
PAD,
*ELSET,ELSET=EALL2
PMASS,
*ELSET,ELSET=EALL3
RMASS,
*SURFACE, TYPE=CYLIN,NAME=DRUM
0.,0.,0.,1.,0.,0.
0.,0.,-1.
START,0.,0.2
CIRCL,0.2,0.,0.,0.
CIRCL,0.,-0.2,0.,0.
CIRCL,-0.2,0.,0.,0.
CIRCL,0.,0.2,0.,0.
*SURFACE,TYPE=ELEMENT, NAME=PAD_LOAD
 PAD, S5
*SURFACE,TYPE=ELEMENT, NAME=PAD
 PAD, S3
*RIGID BODY, REF NODE=100, ANALYTICAL SURFACE =DRUM
*RIGID BODY, ELSET=PLATE, REF NODE=600
*NSET,NSET=QA_TEST_AXIS
AXIS,
**
*SURFACE INTERACTION,NAME=DRUM_PAD
*FRICTION
 0.15,
*STEP
*DYNAMIC, EXPLICIT
 , .2
*CONTACT 
*CONTACT INCLUSIONS
,
*CONTACT PROPERTY ASSIGNMENT
,,DRUM_PAD
**
**  Apply a constant pressure to the pad.
**
*DSLOAD
 PAD_LOAD, P, 350.E3
**
**  Some typical output requests.
**
*RESTART,TIMEMARKS=YES, WRITE, NUMBER=1
**
*FILE OUTPUT, NUMBER=2, TIMEMARKS=YES
*NODE FILE,NSET=AXIS
 U,
**
***OUTPUT,HISTORY,FREQ=1
***NODE OUTPUT,NSET=AXIS
***RF1,RF2,RM3,UR3,VR3,AR3
***ENERGY OUTPUT
**ALLIE,ALLKE,ALLPD,ALLSE,ALLVD,ALLWK,ETOTAL,ALLFD
***ENERGY OUTPUT,ELSET=EALL1
**ALLIE,ALLKE,ALLSE,
***ENERGY OUTPUT,ELSET=EALL2
**ALLIE,ALLKE,ALLSE,
***ENERGY OUTPUT,ELSET=EALL3
**ALLIE,ALLKE,ALLSE,
**
*OUTPUT, FIELD, NUMBER INTERVAL=2, TIMEMARKS=YES
*NODE OUTPUT, NSET=QA_TEST_AXIS
 U,
*OUTPUT,FIELD,VARIABLE=PRESELECT
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIMEINT=0.1
*END STEP