*HEADING
 PIPE WHIP SIMULATION for S4RSW
**
**  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=S4RSW
   1,   1,   2,   18,   17
*ELGEN, ELSET=PIPE1
   1, 15, 1, 1, 46, 16, 15
*SECTION CONTROLS, NAME=HG, HOURGLASS=STIFFNESS
*SHELL SECTION, ELSET=PIPE1, MATERIAL=STEEL, SECTION INTEGRATION=GAUSS,CONTROLS=HG
 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=S4RSW
        5001,        5001,        5002,        5026,        5025
*ELGEN, ELSET=OPENLOOP
        5001, 23, 1, 1, 16, 24, 24
*ELEMENT, TYPE=S4RSW
        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,CONTROLS=HG
 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
*SURFACE,TYPE=ELEMENT, NAME=SURFP1
 PIPE11,
*SURFACE,TYPE=ELEMENT, NAME=SURFP2
 PIPE2,
*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
*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