*HEADING
 1D transient wave propagation, solid medium, incident wave load, spherical
 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=C3D5
 1,   1, 251, 151, 101, 102
11,   1,   2, 252, 251, 102
21, 252, 152, 151, 251, 102
31,   3, 103, 153, 253, 102
41,   2,   3, 253, 252, 102
51, 253, 153, 152, 252, 102
*ELGEN,ELSET=CHANNEL
 1,   5, 2, 1
11,   5, 2, 1
21,   5, 2, 1
31,   5, 2, 1
41,   5, 2, 1
51,   5, 2, 1
**************************
*SURFACE,NAME=IW_LOAD
1,S1
**************************
*SOLID SECTION,MATERIAL=SOLID,ELSET=CHANNEL
**
*MATERIAL,NAME=SOLID
*DENSITY
 1.E3,
*ELASTIC
 1.E6,0.0
**
**************************
*INCIDENT WAVE PROPERTY,NAME=IWPROP,TYPE=SPHERE
0.0,0.0,0.0
-10.E04,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
*DYNAMIC,ALPHA=0.0
.0001,.0015
**
*INCIDENT WAVE,PRESSUREAMPLITUDE=PRESSUREVTIME,PROPERTY=IWPROP
IW_LOAD,415.09517
**
*EL PRINT,FREQ=5
 S,
 E,
*NODE PRINT,NSET=SAMPLE,FREQ=5
 U,
*OUTPUT,VAR=PRESELECT,FIELD,FREQ=1
*OUTPUT,VAR=PRESELECT,HISTORY,FREQ=1
*OUTPUT,FIELD, OP=ADD,FREQUENCY=1
*NODE OUTPUT
 U
*ELEMENT OUTPUT
 S,E
*OUTPUT,HISTORY,OP=ADD,FREQUENCY=1
*NODE OUTPUT,NSET=SAMPLE
 U,
*ENDSTEP