*HEADING
 PIPE WHIP SIMULATION
**  pipewhip_std.inp
**
**  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
*NODE, SYSTEM=C
*NODE, SYSTEM=R
*NGEN, NSET=SECTA, LINE=C
*NCOPY, SHIFT, CHANGE NUMBER=144, OLD SET=SECTA, NEW SET=SECTB
*NFILL, NSET=END1P1
 SECTA, SECTB,           9,          16
*NCOPY, SHIFT, CHANGE NUMBER=592, OLD SET=SECTA, NEW SET=SECTC
*NFILL, NSET=MIDP1
 SECTB, SECTC, 28, 16
*NCOPY, SHIFT, CHANGE NUMBER=736, OLD SET=SECTA, NEW SET=SECTD
*NFILL, NSET=END2P1
 SECTC, SECTD, 9, 16
*NSET, NSET=PIPE1
 END1P1, MIDP1, END2P1
*ELEMENT, TYPE=S4R
*ELGEN, ELSET=PIPE1
*SHELL SECTION, ELSET=PIPE1, MATERIAL=STEEL, SECTION INTEGRATION=GAUSS
 0.432,
*NSET, NSET=SECTE, GENERATE
*BOUNDARY
 SECTE, YSYMM
 REFPT, 1, 1, 0.
 REFPT, 3, 3, 0.
*NSET,NSET=PLANE,GENERATE
*NODE,NSET=REFPT, SYSTEM=C
*ELSET, ELSET=MIDP1, GENERATE
**
** Generate a 16 x 24 mesh for the impacted pipe.
**
*SYSTEM
*NODE, SYSTEM=C
*NODE, SYSTEM=R
*NGEN, NSET=SECTF, LINE=C
*NCOPY, SHIFT, CHANGE NUMBER=168, OLD SET=SECTF, NEW SET=SECTG
*NFILL, NSET=MIDP2
 SECTF, SECTG, 7, 24
*NCOPY, SHIFT, CHANGE NUMBER=384, OLD SET=SECTF, NEW SET=SECTH
*NFILL, NSET=ENDP2
 SECTG, SECTH, 9, 24
*NSET, NSET=PIPE2
 MIDP2, ENDP2
*ELEMENT, TYPE=S4R
*ELGEN, ELSET=OPENLOOP
*ELEMENT, TYPE=S4R
*ELGEN, ELSET=CLOSURE
*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
*ELSET, ELSET=ETOP
*ELSET, ELSET=EBOT
*ELSET, ELSET=ELOUT
ETOP,EBOT
*NSET, NSET=NOUT
**
**  Material description.
**
*MATERIAL, NAME=STEEL
*DAMPING,BETA=1.E-5
*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
*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
*ELEMENT, TYPE=ROTARYI, ELSET=ROTARY
*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,
*SURFACE INTERACTION,NAME=INT1
1.,
*SURFACE BEHAVIOR,PENALTY
*CONTACT PAIR,INTERACTION=INT1
 SURFP1, SURFP2
**
*STEP,NLGEOM,INC=1000,EXTRAPOL=PARABOLIC
*DYNAMIC
1.e-4 , 0.015,, 5.e-4
*OUTPUT,FIELD,freq=1
*NODE OUTPUT
U, V
*ELEMENT OUTPUT
S,PEEQ
*CONTACT OUTPUT
CSTRESS,
*OUTPUT,HISTORY,variable=preselect,freq=1
*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
*END STEP