*HEADING
DERA UNDEX Benchmark Problem 3
Cylindrical Elastic Shell with Planar Step Shockwave
**
** Shell Nodes
**
*NODE
  1,  1.0,0.0,0.0
 65, -1.0,0.0,0.0
998,  0.0,0.0,0.0
999,  0.0,0.0,0.0490873849509
*NGEN,NSET=B,LINE=C
 1,65,1,998,,,,0,0,1
*NCOPY,CHANGE NUMBER=65,OLD SET=B,SHIFT,NEW SET=F
 0.0,0.0,0.0490873849509

*NSET,NSET=BACK,GENERATE
 2,64,1
*NSET,NSET=FRONT,GENERATE
 67,129,1
*NSET,NSET=SYM
 1,65,66,130
*NSET,NSET=TOP
 33,98
*NORMAL
 1,1,1.0,0.0,0.0
 1,66,1.0,0.0,0.0
**
** Water Nodes
**
*NODE
 201,    1.0,0.0,0.0
 265,   -1.0,0.0,0.0
 20066,  1.0,0.0,0.0490873849509
 20130, -1.0,0.0,0.0490873849509
*NGEN,NSET=BWN,LINE=C
 201,265,1,998,,,,0,0,1
*NGEN,NSET=FWN,LINE=C
 20066,20130,1,999,,,,0,0,1
*NODE
 15201,  3.0,0.0,0.0
 15265, -3.0,0.0,0.0
 35066,  3.0,0.0,0.0490873849509
 35130, -3.0,0.0,0.0490873849509
*NGEN,NSET=BWF,LINE=C
 15201,15265,1,998,,,,0,0,1
*NGEN,NSET=FWF,LINE=C
 35066,35130,1,999,,,,0,0,1
*NSET, NSET=QA_TEST
201,202,203,20066,20067,20068,
*NFILL,NSET=BW
 BWN,BWF,25,600
 FWN,FWF,25,600
**
** Shell Elements
**
*ELEMENT,TYPE=S4R
 1,66,1,2,67
*ELGEN,ELSET=SHELL
 1,64,1,1,
**
** Water Elements
**
*ELEMENT,TYPE=AC3D8R
 101,20066,201,202,20067,20666,801,802,20667
*ELGEN,ELSET=WATER
 101,64,1,1,25,600,100
*ELSET,ELSET=INNER,GENERATE
 101,164,1
*ELSET,ELSET=OUTER,GENERATE
 2501,2664,1
**
**
*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
 2.316E9,
*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
 1.0145,0.0,0.0,
 3.0,0.0,0.0
*INCIDENT WAVE FLUID PROPERTY
 2.316E9,997.0
**
**
*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,non=c
 QUIET,2.0
*BULK VISCOSITY
 0.75,10.0
*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