*HEADING
 PIPE WHIP SIMULATION
**
**  This example simulates a pipe-on-pipe impact resulting from 
**  the rupture of a high-pressure line in a power plant. It is 
**  assumed that a sudden release of fluid could cause one segment 
**  of the pipe to rotate about its support and strike a 
**  neighboring pipe.
**
**  Generate a 92 x 30 mesh for the impacting pipe.
**
*SYSTEM
 0., 0., 3.3125, 0., 1., 3.3125
 0., 0., 4.3125
*NODE, SYSTEM=C
    1, 3.0965, -90., -25.
   31, 3.0965,  90., -25.
*NODE, SYSTEM=R
 100000, 0., 0., -25.
*NGEN, NSET=SECTA, LINE=C
    1, 31, 1, 100000, 0., 0., -25., 0., 0., 1.
*NCOPY, SHIFT, CHANGE NUMBER=558, OLD SET=SECTA, NEW SET=SECTB
  15., 0., 0.
   0., 0., 0., 0., 0., 1., 0.
*NFILL, NSET=END1P1
 SECTA, SECTB, 18, 31
*NCOPY, SHIFT, CHANGE NUMBER=2294, OLD SET=SECTA, NEW SET=SECTC
  35., 0., 0.
   0., 0., 0., 0., 0., 1., 0.
*NFILL, NSET=MIDP1
 SECTB, SECTC, 56, 31
*NCOPY, SHIFT, CHANGE NUMBER=2852, OLD SET=SECTA, NEW SET=SECTD
  50., 0., 0.
   0., 0., 0., 0., 0., 1., 0.
*NFILL, NSET=END2P1
 SECTC, SECTD, 18, 31
*NSET, NSET=PIPE1
 END1P1, MIDP1, END2P1
*ELEMENT, TYPE=S4R
   1,   1,   2,   33,  32
*ELGEN, ELSET=PIPE1
   1, 30, 1, 1, 92, 31, 30
*SHELL SECTION, ELSET=PIPE1, MATERIAL=STEEL, SECTION INTEGRATION=GAUSS
 0.432, 
*NSET, NSET=SECTE, GENERATE
  1,        2853,          31
 31,        2883,          31
*BOUNDARY
 SECTE, YSYMM
 REFPT, 1, 1, 0.
 REFPT, 3, 3, 0.
*NSET,NSET=PLANE,GENERATE
        2853,        2883,           1
*NODE,NSET=REFPT, SYSTEM=C
      200001, 0., 0., 25.
*ELSET, ELSET=MIDP1, GENERATE
         541,        2220, 1
**
**  Generate a 32 x 48 mesh for the impacted pipe.
**
*SYSTEM
 0., 0., -3.3125, -1., 0., -3.3125
 0., 0.,  0.
*NODE, SYSTEM=C
        5001, 3.0965,   0., 0.
        5048, 3.0965,   352.5000    , 0.
*NODE, SYSTEM=R
 100001, 0., 0., 0.
*NGEN, NSET=SECTF, LINE=C
         5001, 5048, 1, 100001, 0., 0., 0., 0., 0., 1.
*NCOPY, SHIFT, CHANGE NUMBER=672, OLD SET=SECTF, NEW SET=SECTG
   0.,  7., 0.
   0.,  0., 0., 0., 0., 1., 0.
*NFILL, NSET=MIDP2
 SECTF, SECTG, 14, 48
*NCOPY, SHIFT, CHANGE NUMBER=1536, OLD SET=SECTF, NEW SET=SECTH
   0., 25., 0.
   0.,  0., 0., 0., 0., 1., 0.
*NFILL, NSET=ENDP2
 SECTG, SECTH, 18,  48
*NSET, NSET=PIPE2
 MIDP2, ENDP2
*ELEMENT, TYPE=S4R
        5001,        5001,        5002,        5050,        5049
*ELGEN, ELSET=OPENLOOP
        5001, 47, 1, 1, 32, 48, 48
*ELEMENT, TYPE=S4R
        5048,        5048,        5001,        5049,        5096
*ELGEN, ELSET=CLOSURE
        5048,          32,          48,          48
*ELSET, ELSET=PIPE2
 OPENLOOP, CLOSURE
*SHELL SECTION, ELSET=PIPE2, MATERIAL=STEEL, SECTION INTEGRATION=GAUSS
 0.432, 
*BOUNDARY
 SECTF, YSYMM
 SECTH, ENCASTRE
*ELSET, ELSET=MIDP2, GENERATE
        5001, 5672, 1
*ELSET, ELSET=ETOP
185, 
*ELSET, ELSET=EBOT
602, 
*ELSET, ELSET=ELOUT
ETOP,EBOT
*NSET, NSET=NOUT
200,602
**
**  Material description.
**
*MATERIAL, NAME=STEEL
*ELASTIC
 30.E6, 0.3
*PLASTIC
 45.E3, 
*DENSITY
 7.324E-4, 
**
**  The impacting pipe is allowed to rotate about a fixed pivot 
**  with an initial angular velocity of 75 radian/sec.
**
*INITIAL CONDITIONS, TYPE=ROTATING VELOCITY
 PIPE1, 75., 0., 0., 0.
 25., 3.0965, 3.3125, 25., -3.0965, 3.3125   
**
**  The simulation will run for 0.015 second
**
*RESTART, WRITE, NUMBER INTERVAL=3, TIMEMARKS=NO
**
*ELGEN, ELSET=PIPE11
   1, 30, 1, 1, 91, 31, 30
*SURFACE,TYPE=ELEMENT, NAME=SURFP1
 PIPE11,
*SURFACE,TYPE=ELEMENT, NAME=SURFP2
 PIPE2,
*RIGID BODY,REF NODE=200001,TIE NSET=PLANE
*ELEMENT, TYPE=MASS, ELSET=MASS
200000, 200001
*ELEMENT, TYPE=ROTARYI, ELSET=ROTARY
200001, 200001
*MASS, ELSET=MASS
1.E-6,
*ROTARY INERTIA, ELSET=ROTARY
1.E-7,1.E-7,1.E-7
*STEP
*DYNAMIC, EXPLICIT
 , 0.015
*CONTACT PAIR
 SURFP1, SURFP2
*OUTPUT,FIELD,VAR=PRESELECT,NUMBER INTERVAL=5
**
** Output requests for qa testing
**
*NSET,NSET=QA_TEST
NOUT,
*ELSET,ELSET=QA_TEST
ELOUT,
*OUTPUT,FIELD,NUMBER INTERVAL=1
*NODE OUTPUT, NSET=QA_TEST
 U,
*ELEMENT OUTPUT, ELSET=QA_TEST
 PEEQ,
*OUTPUT,HISTORY,VAR=PRESELECT,TIME INTERVAL=0.00375
*FILE OUTPUT, NUMBER INTERVAL=1
*EL FILE, ELSET=ELOUT
 PEEQ,
*NODE FILE, NSET=NOUT
 U,
*ENERGY FILE
***OUTPUT,HISTORY,TIME INTERVAL=0.003E-3
***ENERGY OUTPUT
**ALLAE,ALLIE,ALLKE,ALLPD,ALLSE,ALLVD,ALLWK,ETOTAL
***NODE OUTPUT,NSET=NOUT
**U,V
***ELEMENT OUTPUT,ELSET=ELOUT
**SF,STH,PRESS
***ELEMENT OUTPUT,ELSET=ELOUT,SECTIONPOINT=1
**MISES,PEEQ,PRESS
***ELEMENT OUTPUT,ELSET=ELOUT,SECTIONPOINT=2
**MISES,PEEQ,PRESS
***ELEMENT OUTPUT,ELSET=ELOUT,SECTIONPOINT=3
**MISES,PEEQ,PRESS
*END STEP