Spinning Bar

A bar is given an initial angular velocity and its free rotation is computed. To eliminate hourglassing, a fully-integrated solid element formulation is used.

*CONTROL_ACCURACY
*CONTROL_ENERGY
*CONTROL_OUTPUT
*CONTROL_TERMINATION
*CONTROL_TIMESTEP
*DATABASE_BINARY_D3PLOT
*DATABASE_GLSTAT
*DATABASE_HISTORY_NODE
*DATABASE_MATSUM
*DATABASE_NODOUT
*ELEMENT_SOLID
*END
*INITIAL_VELOCITY_GENERATION
*KEYWORD
*MAT_ELASTIC
*NODE
*PART
*SECTION_SOLID
*TITLE

  *KEYWORD
  *TITLE
  spinning bar
  $
  $  J.D. Reid    2/11/00 - original model obtained from Jim Day (LSTC)
  $
  $  Units: unknown
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $$$$  Control Ouput
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
  $
  *CONTROL_TERMINATION 
  $   endtim    endcyc     dtmin    endeng    endmas
       0.010
  $
  *CONTROL_TIMESTEP
  $   dtinit      scft      isdo    tslimt     dt2ms      lctm     erode     ms1st
                  0.70
  $
  *CONTROL_ACCURACY
  $      osu       inn
           1         2
  $
  *CONTROL_ENERGY
  $     hgen      rwen    slnten     rylen
           2         2         2
  $
  *CONTROL_OUTPUT
  $    npopt    neecho    nrefup    iaccop     opifs    ipnint    ikedit
           1         3
  $
  $
  *DATABASE_BINARY_D3PLOT
  $       dt      lcdt
      0.0002
  $
  $
  *DATABASE_GLSTAT
  $       dt
      0.0001
  $
  *DATABASE_MATSUM
  $       dt
      0.0001
  $
  *DATABASE_NODOUT
  $       dt
      0.0001
  $
  *DATABASE_HISTORY_NODE
  $      id1       id2       id3       id4       id5       id6       id7       id8
          88       101
  $
  $   88 - node at end of beam (ctr of end)
  $  101 - node at center of model
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $$$$  Initial Conditions
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
  $
  *INITIAL_VELOCITY_GENERATION
  $      sid      styp     omega        vx        vy        vz
           1         2     1000.
  $       xc        yc        zc        nx        ny        nz     phase
         0.0       0.0       0.0       0.0       0.0       1.0
  $
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $$$$  Define Parts and Materials
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
  $
  *PART
  $      pid       sid       mid     eosid      hgid      grav    adpopt      tmid
  bar
           1         1         1
  $
  *MAT_ELASTIC
  $      mid        ro         e        pr
           1  3.00E-04  1.00E+07      0.30
  $
  *SECTION_SOLID
  $      sid    elform
           1         2
  $
  $   - fully-integrated solid element formulation to eliminate hourglassing
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $$$$  Define Nodes and Elements
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
  $
  *END