*HEADING
DERA UNDEX Benchmark Problem 10
Fluid Coupled Concentric Cylinders, Inner Cylinder Air Backed
**
** Shell Nodes
**
*NODE,NSET=PLANE
*NODE
*NGEN,NSET=IB,LINE=C
*NGEN,NSET=IF,LINE=C
*NGEN,NSET=OB,LINE=C
*NGEN,NSET=OF,LINE=C
*NSET,NSET=SYMM,GENERATE
*NSET, NSET=QA_TEST
**
** Water Nodes
**
*NODE,NSET=BP
**
*NGEN,NSET=IFN,LINE=C
*NGEN,NSET=IFF,LINE=C
*NGEN,NSET=OFN,LINE=C
*NGEN,NSET=OFF,LINE=C
**
*NFILL,NSET=BP
 IFN,IFF,15,100
 OFN,OFF,75,100
**
** Front Plane
**
*NCOPY,CHANGE NUMBER=50000,OLD SET=BP,SHIFT,NEW SET=FP
**
** Shell Elements
**
*ELEMENT,TYPE=S4R
*ELGEN,ELSET=CYLINDER_I
*ELGEN,ELSET=CYLINDER_E
**
** Water Elements
**
** This mesh density is 15 elements per 0.2 meter, resulting in an element length of 0.013.
**
*ELEMENT,TYPE=AC3D8R
*ELGEN,ELSET=IWATER
*ELGEN,ELSET=OWATER
*ELSET,ELSET=INEAR,GENERATE
*ELSET,ELSET=IFAR,GENERATE
*ELSET,ELSET=ONEAR,GENERATE
*ELSET,ELSET=OFAR,GENERATE
**
**
*SHELL SECTION, ELSET=CYLINDER_E, MATERIAL=STEEL_E
 0.00581,
*SHELL SECTION, ELSET=CYLINDER_I, MATERIAL=STEEL_I
 0.02324,
*SOLID SECTION, ELSET=IWATER, MATERIAL=FLUID_E
*SOLID SECTION, ELSET=OWATER, MATERIAL=FLUID_E
*MATERIAL, NAME=STEEL_E
*DENSITY
 7766.0,0.0
*ELASTIC, TYPE=ISOTROPIC
 2.0640E+11,0.3,0.0
*MATERIAL,NAME=STEEL_I
 *DENSITY
 7766.0,0.0
*ELASTIC,TYPE=ISOTROPIC
 2.0640E+11,0.3,0.0
**
** NOTE THAT WE HAVE TAKEN THE FLUID PROPERTIES DIRECTLY
** FROM THE DERA DECK.  THESE VALUES DO NOT SATISFY K=C*C*RHO
**
*MATERIAL,NAME=FLUID_I
*ACOUSTIC MEDIUM
*DENSITY
 997.0,
*MATERIAL,NAME=FLUID_E
*ACOUSTIC MEDIUM
*DENSITY
 997.0,
**
**
*AMPLITUDE,NAME=UNDEX,DEFINITION=TABULAR,TIME=TOTAL TIME
0.0,0.0,1.0E-6,1.0,10.0,1.0
**
**
*SURFACE,NAME=FLUID1
 INEAR,S1
*SURFACE,NAME=FLUID2
 IFAR,S2
*SURFACE,NAME=FLUID3
 ONEAR,S1
*SURFACE,NAME=QUIET
 OFAR,S2
*SURFACE,NAME=SOLID1
 CYLINDER_I,SPOS
*SURFACE,NAME=SOLID2
 CYLINDER_E,SNEG
*SURFACE,NAME=SOLID3
 CYLINDER_E,SPOS
**
**
*TIE,NAME=COUPLE
 FLUID1,SOLID1
 FLUID2,SOLID2
 FLUID3,SOLID3
**
**
*BOUNDARY
 SYMM,ZSYMM
 PLANE,2,6
**
**
*INCIDENT WAVE PROPERTY,NAME=WAVE,TYPE=PLANE
*INCIDENT WAVE FLUID PROPERTY
*IMPEDANCE PROPERTY,TYPE=CIRCULAR,NAME=OUTER
**
**
*STEP
*DYNAMIC,EXPLICIT,DIRECT USER CONTROL
 1.0E-7,0.00538
*INCIDENT WAVE,PROPERTY=WAVE,PRESSURE AMPLITUDE=UNDEX
 SOLID3,1.0E+06
 FLUID3,1.0E+06
*SIMPEDANCE
 QUIET,OUTER
*BULK VISCOSITY
*OUTPUT,FIELD
*NODE OUTPUT
 U,
*OUTPUT,FIELD,NUMBER=2
*NODE OUTPUT, NSET=QA_TEST
 U,
*OUTPUT,HISTORY,TIME INTERVAL=1.076E-6
*NODE OUTPUT,NSET=PLANE
 U,V
*END STEP