*HEADING
DERA UNDEX Benchmark Problem 3
Cylindrical Elastic Shell with Planar Step Shockwave
**
** Shell Nodes
**
*NODE
*NGEN,NSET=B,LINE=C
*NCOPY,CHANGE NUMBER=65,OLD SET=B,SHIFT,NEW SET=F
*NSET,NSET=BACK,GENERATE
*NSET,NSET=FRONT,GENERATE
*NSET,NSET=SYM
*NSET,NSET=TOP
*NORMAL
**
** Water Nodes
**
*NODE
*NGEN,NSET=BWN,LINE=C
*NGEN,NSET=FWN,LINE=C
*NODE
*NGEN,NSET=BWF,LINE=C
*NGEN,NSET=FWF,LINE=C
*NSET, NSET=QA_TEST
*NFILL,NSET=BW
 BWN,BWF,25,600
 FWN,FWF,25,600
**
** Shell Elements
**
*ELEMENT,TYPE=S4R
*ELGEN,ELSET=SHELL
**
** Water Elements
**
*ELEMENT,TYPE=AC3D8R
*ELGEN,ELSET=WATER
*ELSET,ELSET=INNER,GENERATE
*ELSET,ELSET=OUTER,GENERATE
**
**
*SHELL SECTION,ELSET=SHELL,MATERIAL=STEEL
 0.02899,3
*SOLID SECTION,ELSET=WATER,MATERIAL=WATER
 1.0,
**
** The material properties for steel and water are taken from
** DERA input deck cylinder_ps_64.inp
** Water properties in that deck given as rho and c
** Calculate the bulk modulus of water as: K=c*c*rho
**
*MATERIAL,NAME=STEEL
*ELASTIC
 2.064E11,0.3
*DENSITY
 7766.0,
*MATERIAL,NAME=WATER
*ACOUSTIC MEDIUM
*DENSITY
 997.0,
**
**
*AMPLITUDE,NAME=UNDEX,TIME=TOTAL TIME,DEFINITION=TABULAR
 0.0,0.0,1.0E-6,1.0,1.0,1.0
**
**
*SURFACE,NAME=FLUID
 INNER,S1
*SURFACE,NAME=SOLID
 SHELL,SPOS
*SURFACE,NAME=QUIET
 OUTER,S2
*TIE,NAME=COUPLE
 FLUID,SOLID
**
**
*BOUNDARY
 FRONT,3,5
 BACK,3,5
 SYM,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.00656
*INCIDENT WAVE,PROPERTY=WAVE,PRESSURE AMPLITUDE=UNDEX
 SOLID,1.0
 FLUID,1.0
*SIMPEDANCE
 QUIET,OUTER
*BULK VISCOSITY
*OUTPUT,FIELD,VARIABLE=PRESELECT
*NODE OUTPUT
 POR,
*OUTPUT, FIELD
*NODE OUTPUT, NSET=QA_TEST
 U,
*OUTPUT,HISTORY,TIME INTERVAL=6.56E-6
*NODE OUTPUT,NSET=SYM
 U,V
*NODE OUTPUT,NSET=TOP
 U,V
*END STEP