*HEADING
 1D transient wave propagation, solid medium, incident wave load, planar
 and defined through pressure amplitude
**
** 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,
  11, .1,
 101, 0, 1E-2
 111, .1, 1E-2
 151, 0, 1E-2, -1E-2
 161, .1, 1E-2, -1E-2
 251, 0, 0, -1E-2
 261, .1, 0, -1E-2
*NGEN,NSET=CHANNEL
   1,11
 101,111
 151,161
 251,261
*NSET,NSET=SAMPLE,GENERATE
 1,11
*ELEMENT,TYPE=C3D4
    1,   1, 102,   2, 251
    2, 102,   2, 251, 152
    3,   2, 251, 152, 252
    4,   1, 101, 102, 251
    5, 101, 102, 251, 151
    6, 102, 251, 151, 152 
*ELGEN,ELSET=CHANNEL
 1,10,1,6
 2,10,1,6
 3,10,1,6
 4,10,1,6
 5,10,1,6
 6,10,1,6
**************************
*SURFACE,NAME=IW_LOAD
4,S2
5,S4
**************************
*SOLID SECTION,MATERIAL=SOLID,ELSET=CHANNEL
**
*MATERIAL,NAME=SOLID
*DENSITY
 1.E3,
*ELASTIC
 1.E6,0.0
**
**************************
*INCIDENT WAVE PROPERTY,NAME=IWPROP,TYPE=PLANE
0.0,0.0,0.0
-10,0.0,0.0
*INCIDENT WAVE FLUID PROPERTY
0.1424E6,1.21,0.0
*AMPLITUDE,DEFINITION=TABULAR,NAME=PRESSUREVTIME
0.0,0.0,0.0015,-1.0
**************************
**
*BOUNDARY
CHANNEL,2,3
 11,1,1
111,1,1
161,1,1
261,1,1
**************************
*RESTART,WRITE
**
*STEP,INC=1000,NLGEOM
*DYNAMIC,ALPHA=0.0
.0001,.0015
**
*INCIDENT WAVE,PRESSUREAMPLITUDE=PRESSUREVTIME,PROPERTY=IWPROP
IW_LOAD,415.09517
**
*EL PRINT,FREQ=5
 S,
*NODE PRINT,NSET=SAMPLE,FREQ=5
 U,
*OUTPUT,VAR=PRESELECT,FIELD,FREQ=1
*OUTPUT,VAR=PRESELECT,HISTORY,FREQ=1
*OUTPUT,FIELD, FREQUENCY=1
*NODE OUTPUT
 U
*ELEMENT OUTPUT
 S,
*OUTPUT,HISTORY,FREQUENCY=1
*NODE OUTPUT,NSET=SAMPLE
 U,
*ENDSTEP