*HEADING
DERA UNDEX Benchmark Problem 4
Cylindrical Elastic Shell with Planar Exponentially Decaying 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
**
** Water Nodes
**
*NODE
 201,    1.014495,0.0,0.0
 265,   -1.014495,0.0,0.0
 20066,  1.014495,0.0,0.0490873849509
 20130, -1.014495,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
*NCOPY,CHANGE NUMBER=15000,OLD SET=BWN,POLE,NEW SET=BWF
 998,,,,
*NCOPY,CHANGE NUMBER=15000,OLD SET=FWN,POLE,NEW SET=FWF
 999,,,,
*NFILL,NSET=BW
 BWN,BWF,150,100
 FWN,FWF,150,100
**
** Shell Elements
**
*ELEMENT,TYPE=S4R
 1,66,1,2,67
*ELGEN,ELSET=SHELL
 1,64,1,1,
**
** 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
 101,20066,201,202,20067,20166,301,302,20167
*ELGEN,ELSET=WATER
 101,64,1,1,150,100,100
*ELSET,ELSET=INNER,GENERATE
 101,164,1
*ELSET,ELSET=OUTER,GENERATE
 15001,15064,1
**
**
*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
 2.316E9,
*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
.27340D-05,.16240D+06,.28194D-05,.17266D+06,.29049D-05,.18322D+06,.29903D-05,.19410D+06
.30757D-05,.20529D+06,.31612D-05,.21678D+06,.32466D-05,.22859D+06,.33321D-05,.24070D+06
.34175D-05,.25312D+06,.35029D-05,.26586D+06,.35884D-05,.27890D+06,.36738D-05,.29224D+06
.37592D-05,.30590D+06,.38447D-05,.31986D+06,.39301D-05,.33413D+06,.40156D-05,.34871D+06
.41010D-05,.36359D+06,.41864D-05,.37878D+06,.42719D-05,.39428D+06,.43573D-05,.41008D+06
.44427D-05,.42619D+06,.45282D-05,.44260D+06,.46136D-05,.45932D+06,.46991D-05,.47634D+06
.47845D-05,.49366D+06,.48699D-05,.51129D+06,.49554D-05,.52922D+06,.50408D-05,.54746D+06
.51263D-05,.56600D+06,.52117D-05,.58484D+06,.52971D-05,.60399D+06,.53826D-05,.62343D+06
.54680D-05,.64318D+06,.55534D-05,.66323D+06,.56389D-05,.68359D+06,.57243D-05,.70424D+06
.58097D-05,.72519D+06,.58952D-05,.74645D+06,.59806D-05,.76800D+06,.60661D-05,.78986D+06
.61515D-05,.81201D+06,.62369D-05,.83447D+06,.63224D-05,.85722D+06,.64078D-05,.88027D+06
.64933D-05,.90362D+06,.65787D-05,.92727D+06,.66641D-05,.95121D+06,.67496D-05,.97546D+06
.68350D-05,.10000D+07,.41010D-04,.77880D+06,.75185D-04,.60653D+06,.10936D-03,.47237D+06
.14354D-03,.36788D+06,.17771D-03,.28650D+06,.21188D-03,.22313D+06,.24606D-03,.17377D+06
.28023D-03,.13534D+06,.31441D-03,.10540D+06,.34858D-03,.82085D+05,.38276D-03,.63928D+05
.41693D-03,.49787D+05,.45111D-03,.38774D+05,.48528D-03,.30197D+05,.51946D-03,.23518D+05
.55363D-03,.18316D+05,.58781D-03,.14264D+05,.62198D-03,.11109D+05,.65616D-03,.86517D+04
.69033D-03,.67379D+04,.72451D-03,.52475D+04,.75868D-03,.40868D+04,.79286D-03,.31828D+04
.82703D-03,.24788D+04,.86121D-03,.19305D+04,.89538D-03,.15034D+04,.92956D-03,.11709D+04
.96373D-03,.91188D+03,.99791D-03,.71017D+03,.10321D-02,.55308D+03,.10663D-02,.43074D+03
.11004D-02,.33546D+03,.11346D-02,.26126D+03,.11688D-02,.20347D+03,.12030D-02,.15846D+03
.12371D-02,.12341D+03,.12713D-02,.96112D+02,.13055D-02,.74852D+02,.13397D-02,.58295D+02
.13738D-02,.45400D+02,.14080D-02,.35358D+02,.14422D-02,.27536D+02,.14764D-02,.21445D+02
.15105D-02,.16702D+02,.15447D-02,.13007D+02,.15789D-02,.10130D+02,.16131D-02,.78893D+01
.16472D-02,.61442D+01,.16814D-02,.47851D+01,.17156D-02,.37267D+01,.17498D-02,.29023D+01
.17839D-02,.22603D+01,.18181D-02,.17603D+01,.18523D-02,.13710D+01,.18865D-02,.10677D+01
.19206D-02,.83153D+00,.19548D-02,.64760D+00,.19890D-02,.50435D+00,.20232D-02,.39279D+00
.20573D-02,.30590D+00,.20915D-02,.23824D+00,.21257D-02,.18554D+00,.21599D-02,.14450D+00
.21940D-02,.11254D+00,.22282D-02,.87642D-01,.22624D-02,.68256D-01,.22966D-02,.53158D-01
.23307D-02,.41399D-01,.23649D-02,.32242D-01,.23991D-02,.25110D-01,.24333D-02,.19556D-01
.24674D-02,.15230D-01,.25016D-02,.11861D-01,.25358D-02,.92374D-02,.25700D-02,.71941D-02
.26041D-02,.56028D-02,.26383D-02,.43635D-02,.26725D-02,.33983D-02,.27067D-02,.26466D-02
.27408D-02,.20612D-02,.27750D-02,.16052D-02,.28092D-02,.12502D-02,.28434D-02,.97362D-03
.28775D-02,.75826D-03,.29117D-02,.59053D-03,.29459D-02,.45991D-03,.29801D-02,.35817D-03
.30142D-02,.27895D-03,.30484D-02,.21724D-03,.30826D-02,.16919D-03,.31168D-02,.13177D-03
.31509D-02,.10262D-03,.31851D-02,.79920D-04,.32193D-02,.62241D-04,.32535D-02,.48474D-04
.32876D-02,.37751D-04,.33218D-02,.29401D-04,.33560D-02,.22897D-04,.33902D-02,.17832D-04
.34243D-02,.13888D-04,.34585D-02,.10816D-04,.34927D-02,.84235D-05,.35269D-02,.65602D-05
**
**
*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.014495,0.0,0.0,
 2.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=circular
 QUIET,2.0
*BULK VISCOSITY
 0.25,10.0
*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