*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
*DYNAMIC,EXPLICIT, DIRECT USER
0.6E-4 ,.0015,
**
*INCIDENT WAVE,PRESSUREAMPLITUDE=PRESSUREVTIME,PROPERTY=IWPROP
IW_LOAD,415.09517
**
*NSET,NSET=QA_TEST
SAMPLE,
*OUTPUT,FIELD,NUMBER INTERVAL=2
*NODE OUTPUT
U,
*NODE OUTPUT,NSET=QA_TEST
U,
*ELEMENT OUTPUT
S,
*OUTPUT,HISTORY,TIME INTERVAL=0.1E-3
*NODE OUTPUT,NSET=SAMPLE
 U,
*ENDSTEP