Contact Friction II

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