*HEADING
 PIPE WHIP SIMULATION
*preprint, model=yes, contact=yes, history=yes
**
**  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 46 x 15 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.
   16, 3.0965,  90., -25.
*NODE, SYSTEM=R
 100000, 0., 0., -25.
*NGEN, NSET=SECTA, LINE=C
    1, 16, 1, 100000, 0., 0., -25., 0., 0., 1.
*NCOPY, SHIFT, CHANGE NUMBER=144, OLD SET=SECTA, NEW SET=SECTB
  15., 0., 0.
   0., 0., 0., 0., 0., 1., 0.
*NFILL, NSET=END1P1
 SECTA, SECTB,           9,          16
*NCOPY, SHIFT, CHANGE NUMBER=592, OLD SET=SECTA, NEW SET=SECTC
  35., 0., 0.
   0., 0., 0., 0., 0., 1., 0.
*NFILL, NSET=MIDP1
 SECTB, SECTC, 28, 16
*NCOPY, SHIFT, CHANGE NUMBER=736, OLD SET=SECTA, NEW SET=SECTD
  50., 0., 0.
   0., 0., 0., 0., 0., 1., 0.
*NFILL, NSET=END2P1
 SECTC, SECTD, 9, 16
*NSET, NSET=PIPE1
 END1P1, MIDP1, END2P1
*ELEMENT, TYPE=S4R
   1,   1,   2,   18,   17
*ELGEN, ELSET=PIPE1
   1, 15, 1, 1, 46, 16, 15
*SHELL SECTION, ELSET=PIPE1, MATERIAL=STEEL, SECTION INTEGRATION=GAUSS
 0.432, 
*NSET, NSET=SECTE, GENERATE
  1,         737,          16
 16,         752,          16
*BOUNDARY
 SECTE, YSYMM
 REFPT, 1, 1, 0.
 REFPT, 3, 3, 0.
*NSET,NSET=PLANE,GENERATE
737, 752, 1
*NODE,NSET=REFPT, SYSTEM=C
 1001, 0., 0., 25.
*ELSET, ELSET=MIDP1, GENERATE
            136, 555, 1
**
** Generate a 16 x 24 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.
        5024, 3.0965,   345.0000    , 0.
*NODE, SYSTEM=R
 100001, 0., 0., 0.
*NGEN, NSET=SECTF, LINE=C
         5001, 5024, 1, 100001, 0., 0., 0., 0., 0., 1.
*NCOPY, SHIFT, CHANGE NUMBER=168, OLD SET=SECTF, NEW SET=SECTG
   0.,  7., 0.
   0.,  0., 0., 0., 0., 1., 0.
*NFILL, NSET=MIDP2
 SECTF, SECTG, 7, 24
*NCOPY, SHIFT, CHANGE NUMBER=384, OLD SET=SECTF, NEW SET=SECTH
   0., 25., 0.
   0.,  0., 0., 0., 0., 1., 0.
*NFILL, NSET=ENDP2
 SECTG, SECTH, 9, 24
*NSET, NSET=PIPE2
 MIDP2, ENDP2
*ELEMENT, TYPE=S4R
        5001,        5001,        5002,        5026,        5025
*ELGEN, ELSET=OPENLOOP
        5001, 23, 1, 1, 16, 24, 24
*ELEMENT, TYPE=S4R
        5024,        5024,        5001,        5025,        5048
*ELGEN, ELSET=CLOSURE
        5024,          16,          24,          24
*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, 5168, 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, 15, 1, 1, 45, 16, 15
*RIGID BODY,REF NODE=1001,TIE NSET=PLANE
*ELEMENT, TYPE=MASS, ELSET=MASS
10000, 1001
*ELEMENT, TYPE=ROTARYI, ELSET=ROTARY
10001, 1001
*MASS, ELSET=MASS
1.E-6,
*ROTARY INERTIA, ELSET=ROTARY
1.E-7,1.E-7,1.E-7
*ELSET,ELSET=QA_TEST_ELOUT
ELOUT,
*NSET,NSET=QA_TEST_NOUT
NOUT,
**
*STEP
*DYNAMIC, EXPLICIT
 , 0.015
*CONTACT
*CONTACT INCLUSIONS, ALL ELEMENT BASED
*FILE OUTPUT, NUMBER INTERVAL=2, TIMEMARKS=YES
*EL FILE, ELSET=ELOUT
 PEEQ,
*NODE FILE, NSET=NOUT
 U,
*ENERGY FILE
*OUTPUT,FIELD,OP=NEW,NUMBER INTERVAL=6,TIMEMARKS=NO
*ELEMENT OUTPUT,ELSET=ELOUT
 S,PEEQ
*NODE OUTPUT
 U, 
*OUTPUT,HISTORY,OP=NEW,TIME INTERVAL=3.E-5
*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
 MISES,PEEQ,PRESS
*OUTPUT, FIELD, NUMBER INTERVAL=2, TIMEMARKS=YES
*ELEMENT OUTPUT, ELSET=QA_TEST_ELOUT
 PEEQ,
*NODE OUTPUT, NSET=QA_TEST_NOUT
 U,
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIMEINT=0.0075
*END STEP