*HEADING
 1D transient wave propagation, acoustic medium, incident wave load, spherical
 and defined through pressure amplitude
**
** C = root (K/rho)= 343.05 m/s   K=bulk mod, rho=density.  
**
** Analytical solution to prescribed constant acceleration of a0 at x=0
** 
** p(x,t) = (K a0 / c^2) (ct - x)  for x < ct
**        = 0  			   for x > ct
** 
*PREPRINT,ECHO=NO,HISTORY=NO,MODEL=YES
*NODE, NSET=ODB_NSET
    1, 0,
  101, 1,
  501, 0,.5E-2
  601, 1,.5E-2
 1001, 0, 1.E-2
 1101, 1, 1.E-2
 1501, 0, 1.E-2, -1.E-2
 1601, 1, 1.E-2, -1.E-2
 2001, 0, 1.E-2, -0.5E-2
 2101, 1, 1.E-2, -0.5E-2
 2501, 0, 0, -1.E-2
 2601, 1, 0, -1.E-2
 3001, 0,.5E-2, -1.E-2
 3101, 1,.5E-2, -1.E-2
 3501, 0, 0, -0.5E-2
 3601, 1, 0, -0.5E-2
 4001, 0, 0.5E-2, -0.5E-2
 4101, 1, 0.5E-2, -0.5E-2
*NGEN,NSET=CHANNEL
    1,101
  501,601
 1001,1101
 1501,1601
 2001,2101
 2501,2601
 3001,3101
 3501,3601
 4001,4101
*NSET,NSET=SAMPLE,GENERATE
 1,101
*ELEMENT,TYPE=AC3D10
 1,    1, 2501, 1001,    3, 3501, 4001,  501,
            2, 3502,  502,
 2, 2501, 1001,    3, 2503, 4001,  502, 3502,
         2502, 4002, 3503,
 3, 1001,    3, 2503, 1003,  502, 3503, 4002,
         1002,  503, 4003,
 4, 1001, 2501, 1501, 1003, 4001, 3001, 2001,
         1002, 4002, 2002,
 5, 2501, 1501, 1003, 2503, 3001, 2002, 4002,
         2502, 3002, 4003,
 6, 1501, 1003, 2503, 1503, 2002, 4003, 3002,
         1502, 2003, 3003,
*ELGEN,ELSET=CHANNEL
 1,50,2,6
 2,50,2,6
 3,50,2,6
 4,50,2,6
 5,50,2,6
 6,50,2,6
**************************
*SURFACE,NAME=IW_LOAD
 4,S1
 1,S1
**************************
*SOLID SECTION,MATERIAL=AIR,ELSET=CHANNEL
**
*MATERIAL,NAME=AIR
*DENSITY
 1.21,
*ACOUSTIC MEDIUM,BULK MODULUS
 0.1424E6,
**
**************************
*INCIDENT WAVE PROPERTY,NAME=IWPROP,TYPE=SPHERE
0.0,0.0,
-10.0E04,0.0,
*INCIDENT WAVE FLUID PROPERTY
0.1424E6,1.21,0.0
*AMPLITUDE,DEFINITION=TABULAR,NAME=PRESSUREVTIME
0.0,0.0,1.0,-1.0
**************************
*RESTART,WRITE
**
*STEP,INC=1000
*DYNAMIC
.0001,.0044
** DEFINE NONREFLECTIVE BOUNDARY AT X=1
*IMPEDANCE
 297,I3
 300,I3
**
*INCIDENT WAVE,PRESSUREAMPLITUDE=PRESSUREVTIME,PROPERTY=IWPROP
IW_LOAD,415.09517
**
*EL PRINT,FREQ=0
*NODE PRINT,NSET=SAMPLE,FREQ=22
 POR,
*OUTPUT,FIELD,FREQUENCY=1
*NODE OUTPUT
 POR,
*OUTPUT,HISTORY,FREQUENCY=1
*NODE OUTPUT,NSET=SAMPLE
 POR,
*ENDSTEP