*HEADING
SUBMERGED CYLINDER UNDEX (KWON & FOX, 1993)
**
**
**  **********************
**  ***  DEFINE PARTS  ***
**  **********************
**
*PART, NAME=WATER
*INCLUDE, INPUT = submerged_cyl_water.inp
*END PART
**
*PART, NAME=CYLINDER
*INCLUDE, INPUT = submerged_cyl_cylinder.inp
*END PART
**
*PART, NAME=PULSE
*NODE, NSET=PT-MASS
*ELEMENT, TYPE=MASS, ELSET=PT-MASS
*MASS, ELSET=PT-MASS
1.0E-7,
*END PART
**
**
**  ******************************
**  ***  DEFINE ASSEMBLY DATA  ***
**  ******************************
**
*ASSEMBLY, NAME=UNDEX
**
*INSTANCE, NAME=WATER-1, PART=WATER
*END INSTANCE
**
*INSTANCE, NAME=CYLINDER-1, PART=CYLINDER
*END INSTANCE
**
*INSTANCE, NAME=SHOCK-PULSE, PART=PULSE
*END INSTANCE
**
**
**  ****************************
**  **  SURFACE DEFINITIIONS  **
**  **    (ASSEMBLY LEVEL)    **
**  ****************************
**
*SURFACE, TYPE=ELEMENT, NAME=CYL-FSI, NO THICK
CYLINDER-1.ALL-ELEM, SPOS
**
*SURFACE, TYPE=ELEMENT, NAME=RAD-MID
WATER-1._RAD-MID_S3, S3
WATER-1._RAD-MID_S1, S1
WATER-1._RAD-MID_S2, S2
WATER-1._RAD-MID_S4, S4
*SURFACE, TYPE=ELEMENT, NAME=RAD-ENDL
WATER-1._RAD-ENDL_S3, S3
WATER-1._RAD-ENDL_S4, S4
WATER-1._RAD-ENDL_S2, S2
WATER-1._RAD-ENDL_S1, S1
*SURFACE, TYPE=ELEMENT, NAME=RAD-ENDR
WATER-1._RAD-ENDR_S3, S3
WATER-1._RAD-ENDR_S4, S4
WATER-1._RAD-ENDR_S2, S2
WATER-1._RAD-ENDR_S1, S1
*SURFACE, TYPE=ELEMENT, NAME=WATER-FSI
WATER-1._WATER-FSI_S3, S3
WATER-1._WATER-FSI_S4, S4
WATER-1._WATER-FSI_S2, S2
WATER-1._WATER-FSI_S1, S1
**
**
**   ***************************************
**   **    ASSEMBLY LEVEL DEFINITION OF   **
**   ** FLUID-STRUCTURE COUPLING VIA *TIE **
**   ***************************************
**
*TIE, NAME=FSI, ADJUST=NO, NO ROTATION
CYL-FSI, WATER-FSI
**
*END ASSEMBLY
**
**
**  *********************
**  ***  MODEL DATA   ***
**  *********************
**
*MATERIAL, NAME=ALUMINIUM
*DENSITY
2784.5,
*ELASTIC
7.56E+10, 0.33
*PLASTIC
  3E+08,0.
3.3E+08,1.
*DAMPING, ALPHA=750.0
**
**  WATER (C = 1463 METERS/SECOND)
**
*MATERIAL, NAME=WATER
*ACOUSTIC MEDIUM
*DENSITY
1000.0,
**
**
**     **********************************
**     ***  INCIDENT WAVE DEFINITION  ***
**     **********************************
**
*INCIDENT WAVE PROPERTY,NAME=SHOCK,TYPE=SPHERE
*INCIDENT WAVE FLUID PROPERTY
**
**
**  ****************************
**  **   SURFACE IMPEDANCES   **
**  ** (RADIATING BOUNDARIES) **
**  ****************************
**
*IMPEDANCE PROPERTY,NAME=IMP-CYL,TYPE=CIRCULAR
*IMPEDANCE PROPERTY,NAME=IMP-SPH,TYPE=SPHERE
**
**
**  ************************************
**  **  RESTRAIN POINT MASS 1 & 2 DOF **
**  ************************************
**
*BOUNDARY, TYPE=VELOCITY
UNDEX.SHOCK-PULSE.PT-MASS, 1, 2, 0.0
**
**
**  **********************************
**  **  EXPLICIT DYNAMIC STEP DATA  **
**  **********************************
**
*NSET, NSET=QA_TEST
UNDEX.CYLINDER-1.ALL-NODE
*********
*STEP, NAME=UNDEX-SHOCK-WAVE
INCIDENT WAVE, CYLINDER ALPHA DAMPING = 750
*DYNAMIC, EXPLICIT
, 0.008
**
*BULK VISCOSITY
**
*RESTART,WRITE,NUMBERINTERVAL=1,TIME MARKS=NO
**
*INCLUDE, INPUT = submerged_cyl_pulse.inp
**
**
**  ***************************************
**  **  POINT MASS LOADING TO REPRODUCE  **
**  **   THE SHOCK PULSE HISTORY CURVE   **
**  ***************************************
**
*CLOAD, AMPLITUDE=PULSE
UNDEX.SHOCK-PULSE.PT-MASS, 3, 1.0e-7
**
**
**   *************************************
**   **    INCIDENT WAVE LOADING AT     **
**   **  THE FLUID-STRUCTURE INTERFACE  **
**   *************************************
**
*INCIDENT WAVE, PROPERTY=SHOCK,
PRESSURE AMPLITUDE=PULSE
UNDEX.WATER-FSI, 1.
UNDEX.CYL-FSI, 1.
**
**
**  *************************************
**  **  ACOUSTIC RADIAITON BOUNDARIES  **
**  *************************************
**
*SIMPEDANCE
UNDEX.RAD-MID,  IMP-CYL
UNDEX.RAD-ENDL, IMP-SPH
UNDEX.RAD-ENDR, IMP-SPH
**
**
**      *************************
**      ***  OUTPUT REQUESTS  ***
**      *************************
**
*OUTPUT, FIELD, NUMBERINTERVAL=40
*ELEMENT OUTPUT,
ELSET=UNDEX.CYLINDER-1.ALL-ELEM
S, LE, PEEQ
*NODE OUTPUT, NSET=UNDEX.CYLINDER-1.ALL-NODE
U,
*NODE OUTPUT, NSET=UNDEX.WATER-1.ALL-NODE
POR,
***************
*OUTPUT, FIELD, NUM=1
*NODE OUTPUT, NSET=QA_TEST
U,
*****************
*OUTPUT, HISTORY, TIMEINTERVAL=1.0E-05
*NODE OUTPUT, NSET=UNDEX.SHOCK-PULSE.PT-MASS
A3,
*NODE OUTPUT, NSET=UNDEX.CYLINDER-1.N-R-CTR
U,V
*NODE OUTPUT, NSET=UNDEX.CYLINDER-1.N-L-CTR
U,V
*NODE OUTPUT, NSET=UNDEX.CYLINDER-1.N-B1
U,V
*NODE OUTPUT, NSET=UNDEX.CYLINDER-1.N-B2
U,V
*NODE OUTPUT, NSET=UNDEX.CYLINDER-1.N-B3
U,V
*NODE OUTPUT, NSET=UNDEX.CYLINDER-1.N-B4
U,V
*ELEMENT OUTPUT, ELSET=UNDEX.CYLINDER-1.A1
LE11, LE22, PEEQ
*ELEMENT OUTPUT, ELSET=UNDEX.CYLINDER-1.A2
LE11, LE22, PEEQ
*ELEMENT OUTPUT, ELSET=UNDEX.CYLINDER-1.B1
LE11, LE22, PEEQ
*ELEMENT OUTPUT, ELSET=UNDEX.CYLINDER-1.B2
LE11, LE22, PEEQ
*ELEMENT OUTPUT, ELSET=UNDEX.CYLINDER-1.B3
LE11, LE22, PEEQ
*ELEMENT OUTPUT, ELSET=UNDEX.CYLINDER-1.C1
LE11, LE22, PEEQ
*ELEMENT OUTPUT, ELSET=UNDEX.CYLINDER-1.C2
LE11, LE22, PEEQ
**
*ENERGY OUTPUT,ELSET=UNDEX.CYLINDER-1.ALL-ELEM
ALLAE, ALLIE, ALLKE, ALLPD, ALLVD, ALLSE
**
*END STEP