Contact Pipes II

Two pipes consisting of plastic material are modeled. The upper one is given an initial angular velocity and crushes into the stiffer lower one. Surface to surface contact is employed. Adaptive process is controlled by a total angle change criteria. If total angle change in the surrounding of an element exceeds ten degrees, the element is refined to the second level and if it exceeds twenty degrees, the element is refined to the third and last level.

*BOUNDARY_SPC_NODE
*CONTACT_SURFACE_TO_SURFACE
*CONTROL_ADAPT
*CONTROL_CONTACT
*CONTROL_ENERGY
*CONTROL_TERMINATION
*CONTROL_TIMESTEP
*DATABASE_BINARY_D3DUMP
*DATABASE_BINARY_D3PLOT
*DATABASE_BINARY_D3THDT
*DATABASE_ELOUT
*DATABASE_GLSTAT
*DATABASE_MATSUM
*DATABASE_NODOUT
*DATABASE_RCFORC
*ELEMENT_SHELL
*END
*INITIAL_VELOCITY_GENERATION
*KEYWORD
*MAT_PLASTIC_KINEMATIC
*NODE
*PART
*SECTION_SHELL
*SET_PART
*TITLE

  *KEYWORD
  *TITLE
  pipe whip simulation
  $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
  $ mm, ton, s, N
  $ stationary pipe (material 1), Steel
  $ rotation pipe (materia 2), Alum.
  $ A. Tabiei, 2-1-98
  $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
  *CONTROL_TERMINATION
  $   ENDTIM    ENDCYC     DTMIN    ENDNEG    ENDMAS
   0.002E+00         0     0.000     0.000     0.000
  *CONTROL_TIMESTEP
  $   DTINIT      SCFT 
       0.000     0.900
  *CONTROL_CONTACT
  $   SLSFAC    RWPNAL    ISLCHK    SHLTHK
       0.100                   2         2
  $   USRSTR    USRFAC     NSBCS    INTERM
   
  *CONTROL_ENERGY
  $     HGEN      RWEN    SLNTEN     RYLEN
           2         1         1         1
  $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
  *DATABASE_MATSUM
   0.100E-05
  *DATABASE_NODOUT
   0.100E-05
  *DATABASE_ELOUT
   0.100E-05
  *DATABASE_GLSTAT
   0.100E-05
  *DATABASE_RCFORC
   0.100E-05
  $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
  *DATABASE_BINARY_D3PLOT
  $  DT/CYCL
  0.0001E+00
  *DATABASE_BINARY_D3THDT
  $  DT/CYCL
   0.100E+03
  *DATABASE_BINARY_D3DUMP
  $  DT/CYCL
   1.000E+05 
  $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
  *PART 
  LOWER, stationary
  $      PID       SID       MID
           1         1         1
  *PART 
  UPPER, rotating 
  $      PID       SID       MID                                  ADPOPT
           2         2         2                                       1
  $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
  *MAT_Plastic_kinematic
  $      MID        RO         E        PR      SIGY      ETAN      BETA
           1 7.890E-09 2.100E+05 3.000E-01 4.400E+02 1.700E+03 1.000E+00
  $      SRC       SRP        FS
  
  *MAT_Plastic_kinematic
  $      MID        RO         E        PR      SIGY      ETAN      BETA
           2 2.730E-09 7.2408E+4 3.000E-01 9.500E+01 7.2408E+3 0.000E+00
  $      SRC       SRP        FS
  
  $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
  *SECTION_SHELL                                                         
  $      SID    ELFORM      SHRF       NIP
           1         2 0.000E+00       3.0 
  $       T1        T2        T3        T4
   1.000E+01 1.000E+01 1.000E+01 1.000E+01
  *SECTION_SHELL 
  $      SID    ELFORM      SHRF       NIP
           2         2 0.000E+00       3.0
  $       T1        T2        T3        T4
   1.000E+01 1.000E+01 1.000E+01 1.000E+01
  $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
  *INITIAL_VELOCITY_GENERATION
           2         2     200.0 
       635.0     -635.   168.275        -1         0         0
  $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
  *CONTACT_surface_to_surface
           5         4         2         2
   0.0000000 0.0000000 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
  *set_part
  5
  2
  *set_part
  4
  1
  $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
  *CONTROL_ADAPT
  $     freq       tol       opt       lvl     birth     death
       0.001      10.0         2         3       0.0       0.0
  $             ADPASS
                     0
  $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
  *NODE
  $
  *ELEMENT_SHELL
  $
  $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
  *BOUNDARY_SPC_NODE
  $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
  *END