*HEADING
 FRICTIONAL BRAKING OF A ROTATING RIGID BODY IN 2D
 TEST 9-DIGIT LABELS    
*preprint, echo=yes,hist=yes,model=yes  
**
**  Define a rigid circular drum.
**
*NODE, NSET=AXIS
900000100,  0.000,  0.000
*NODE
900000101,  0.200,  0.000
900000119,  0.000,  0.200
900000137, -0.200,  0.000
900000155,  0.000, -0.200
900000173,  0.200,  0.000
*NGEN, LINE=C
900000101,900000119, 1,900000100
900000119,900000137, 1,900000100
900000137,900000155, 1,900000100
900000155,900000173, 1,900000100
*ELEMENT, TYPE=R2D2, ELSET=DRUM
900000101,900000101,900000102
900000172,900000172,900000101
*ELGEN, ELSET=DRUM
900000101, 71, 1, 1
*ELEMENT, TYPE=MASS, ELSET=PMASS
900000100,900000100
*MASS, ELSET=PMASS
 5., 
*ELEMENT, TYPE=ROTARYI, ELSET=RMASS
900000200,900000100
*ROTARY INERTIA, ELSET=RMASS
 0., 0., .175
**
**  Start with it spinning at 60 rad/sec.
**
*BOUNDARY
900000100, PINNED
*INITIAL CONDITIONS, TYPE=VELOCITY
 900000100, 6, 60.
**
**  Define a deformable brake pad.
**
*NODE
900000501, -0.200, -0.050
900000502, -0.250, -0.050
900000521, -0.200,  0.050
900000522, -0.250,  0.050
900000600, -0.250,  0.000
*NGEN
900000501,900000521, 2
900000502,900000522, 2
*ELEMENT, TYPE=CPE4R, ELSET=PAD
900000501,900000501,900000503,900000504,900000502
*ELGEN, ELSET=PAD
900000501,  10, 2, 1
*SOLID SECTION, ELSET=PAD, MATERIAL=STUFF
 0.100, 
**
**  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=R2D2, ELSET=PLATE
900000601,900000502,900000504
*ELGEN, ELSET=PLATE
900000601,  10, 2, 1
*BOUNDARY
900000600, 2, 2, 0.
900000600, 6, 6, 0.
**
**  Run the simulation for 0.2 seconds.
**
*ELSET,ELSET=EALL1
PAD
*ELSET,ELSET=EALL2
PMASS
*ELSET,ELSET=EALL3
RMASS
**
*SURFACE,TYPE=ELEMENT, NAME=DRUM
 DRUM, SNEG
*SURFACE,TYPE=ELEMENT, NAME=PAD_LOAD
 PAD, S3
*SURFACE,TYPE=ELEMENT, NAME=PAD
 PAD, S1
*RIGID BODY,ELSET=DRUM,REFNODE=900000100
*RIGID BODY,ELSET=PLATE,REFNODE=900000600
*NSET,NSET=QA_TEST_AXIS
AXIS
**
*STEP
*DYNAMIC, EXPLICIT
 , .2
**
**  Define contact between the drum and the pad.
**
*SURFACE INTERACTION,NAME=DRUM_PAD
*FRICTION
 0.15,
*CONTACT PAIR,INTERACTION=DRUM_PAD
 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
**
**
***OUTPUT,HISTORY,TIME INTERVAL=0.
***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