This models a drum rotating with initial angular velocity. The drum is brought to rest by a rubber pad pressed against it. From this example you can learn about friction, rubber, rigid material and angular velocity.
*BOUNDARY_SPC_SET *CONSTRAINED_EXTRA_NODES_NODE *CONSTRAINED_NODE_SET *CONTACT_NODES_TO_SURFACE *CONTROL_CONTACT *CONTROL_ENERGY *CONTROL_TERMINATION *CONTROL_TIMESTEP *DATABASE_BINARY_D3PLOT *DATABASE_GLSTAT *DATABASE_HISTORY_NODE_SET *DATABASE_MATSUM *DATABASE_NODOUT *DATABASE_RCFORC *DEFINE_CURVE *ELEMENT_SHELL *ELEMENT_SOLID *END *INITIAL_VELOCITY_GENERATION *KEYWORD *LOAD_NODE_SET *MAT_MOONEY-RIVLIN_RUBBER *MAT_RIGID *NODE *PART *SECTION_SHELL *SECTION_SOLID *SET_NODE_LIST *TITLE
*KEYWORD
*TITLE
BRAKING PROBLEM (N2S, RUBBER PAD)
$ learn friction, rubber and rigid material, and angular velocity
$ drum with radius of 200 mm width of 150 mm and thickness of 5 mm
$ is rotating with initial angulay velocity. The drum is to be brought to
$ rest by rubber brake pad 50x100x100 mm. The normal force pressing the pad
$ against the drum is 3498 N. The friction coeffecient is 0.15.
$ units: m, s, Kg, N
$ A. Tabiei 3/3/98
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*CONTROL_TIMESTEP
0.0000 0.9000
*CONTROL_TERMINATION
0.2500
*CONTROL_CONTACT
$ SLSFAC RWPNAL ISLCHK SHLTHK PENOPT THKCHG ORIEN
2 2
$ USRSTR USRFAC NSBCS INTERM XPENE
*CONTROL_ENERGY
$ hgen rwen slnten rylen
2 2
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*DATABASE_BINARY_D3PLOT
0.01
*DATABASE_MATSUM
0.100E-04
*DATABASE_NODOUT
0.100E-04
*DATABASE_RCFORC
0.100E-02
*DATABASE_GLSTAT
0.100E-04
*DATABASE_HISTORY_NODE_SET
10
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*MAT_RIGID
$ MID RO E PR
1 7.890E+03 2.100E+11 3.000E-01
$ CMO CON CON2
1.000E+00 7.000E+00 6.000E+00
$ alternative method for specifying local system below (optional)
*MAT_MOONEY-RIVLIN_RUBBER
$ MID RO PR A B
2 5.000E+03 4.900E-01 5.520E+05 1.380E+05
$ if A=B=0, then a least square fit is computed from test based on data below
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*SECTION_SHELL
1 2 0.000E+00 3.000E+00 0.000E+00 0.000E+00
5.000E-03 5.000E-03 5.000E-03 5.000E-03 0.000E+00
*SECTION_SOLID
2 1
*PART
DRUM:
1 1 1
*PART
BRAKE-PAD:
2 2 2
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*CONSTRAINED_NODE_SET
$ NSID DOF
11 1
*CONSTRAINED_EXTRA_NODES_NODE
$ PID NID/NSID
1 259
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*INITIAL_VELOCITY_GENERATION
$ ID STYP OMEGA
1 2 3.000E+02
$ XC YC ZC NX NY NZ PHASE
0.000E+00 0.000E+00 0.000E+00 0.000E+00 1.000E+00 0.000E+00 0
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*SET_NODE_LIST
10
259
*SET_NODE_LIST
11
215 216 219 220 223 224 227 228
231 232 235 236 239 240 243 244
247 248 251 252 255 256
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*CONTACT_nodes_TO_surface
2 1 3 3 0 0
0.1500000 0.1500000 0.0000000 0.0000000 0.0000000 0 0.0000000 1.000E+20
1.0000000 1.0000000 0.0000000 0.0000000 1.0000000 1.0000000 1.0000000 1.0000000
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*DEFINE_CURVE
1
0.00000000000000E+00 1.0000000000000E+00
1.00000000000000E+04 1.0000000000000E+00
*LOAD_NODE_SET
$ NID/NSID DOF LCID SF CID
11 1 1 -1.59E+02 0
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*BOUNDARY_SPC_SET
$ NID/NSID CID DOFX DOFY DOFZ DOFRX DOFRY DOFRZ
11 0 0 1 1 0 0 0
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*NODE
$
*ELEMENT_SHELL
$
*ELEMENT_SOLID
$
*END
