*HEADING
DERA UNDEX Benchmark Problem 4
Cylindrical Elastic Shell with Planar Exponentially Decaying 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
**
** Water Nodes
**
*NODE
*NGEN,NSET=BWN,LINE=C
*NGEN,NSET=FWN,LINE=C
*NCOPY,CHANGE NUMBER=15000,OLD SET=BWN,POLE,NEW SET=BWF
*NCOPY,CHANGE NUMBER=15000,OLD SET=FWN,POLE,NEW SET=FWF
*NFILL,NSET=BW
 BWN,BWF,150,100
 FWN,FWF,150,100
**
** Shell Elements
**
*ELEMENT,TYPE=S4R
*ELGEN,ELSET=SHELL
**
** Water Elements
**
** This mesh density is 150 elements per meter, resulting in an element length of 0.00667.
** Note that tau*c=0.208331, so we have approximately 30 elements per decay length.
**
*ELEMENT,TYPE=AC3D8R
*ELGEN,ELSET=WATER
*ELSET,ELSET=INNER,GENERATE
*ELSET,ELSET=OUTER,GENERATE
**
**
*SHELL SECTION,ELSET=SHELL,MATERIAL=STEEL
 0.02899,
*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,DEFINITION=TABULAR,TIME=TOTAL TIME
** Rise time of Tau/Twenty
.00000D+00,.00000D+00,.85437D-07,.16014D+03,.17087D-06,.64035D+03,.25631D-06,.14403D+04
.34175D-06,.25598D+04,.42719D-06,.39985D+04,.51263D-06,.57560D+04,.59806D-06,.78321D+04
.68350D-06,.10226D+05,.76894D-06,.12939D+05,.85437D-06,.15969D+05,.93981D-06,.19316D+05
.10253D-05,.22981D+05,.11107D-05,.26962D+05,.11961D-05,.31260D+05,.12816D-05,.35874D+05
.13670D-05,.40804D+05,.14524D-05,.46049D+05,.15379D-05,.51610D+05,.16233D-05,.57486D+05
.17087D-05,.63677D+05,.17942D-05,.70182D+05,.18796D-05,.77001D+05,.19651D-05,.84134D+05
.20505D-05,.91580D+05,.21359D-05,.99340D+05,.22214D-05,.10741D+06,.23068D-05,.11580D+06
.23922D-05,.12450D+06,.24777D-05,.13351D+06,.25631D-05,.14283D+06,.26486D-05,.15246D+06
**
**
*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
*NODE OUTPUT
 U,
*OUTPUT,HISTORY,TIME INTERVAL=3.28E-5
*NODE OUTPUT,NSET=SYM
 U,V
*NODE OUTPUT,NSET=TOP
 U,V
*END STEP