*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 z=0
** 
** p(z,t) = (K a0 / c^2) (ct - z)  for z < ct
**        = 0  			   for z > ct
** 
*PREPRINT,ECHO=NO,HISTORY=NO,MODEL=YES
*NODE, NSET=ODB_NSET
   1, 0,
  11, 0, .1,
  51, .5E-2,0.
  61, .5E-2, .1
 101, 1.E-2,0.
 111, 1.E-2, .1
*NGEN,NSET=CHANNEL
   1,11
 101,111
*NSET,NSET=SAMPLE,GENERATE
 1,11
*ELEMENT,TYPE=CAX4R
 1,1,101,102,2
*ELGEN,ELSET=CHANNEL
 1,10,1,1
**************************
*SURFACE,NAME=IW_LOAD
1,S1
**************************
** CROSS-SECTION IS Pi*1.E-4
*SOLID SECTION,MATERIAL=SOLID,ELSET=CHANNEL
 1.E-2, 
**
*MATERIAL,NAME=SOLID
*DENSITY
 1.E3,
*ELASTIC
 1.E6,0.0
**
**************************
*INCIDENT WAVE PROPERTY,NAME=IWPROP,TYPE=SPHERE
0.0,0.0,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,0.0015,-1.0
**************************
**
*BOUNDARY
SAMPLE,XSYMM
11,2,2
111,2,2
**************************
*RESTART,WRITE
**
*STEP
*DYNAMIC,EXPLICIT
 ,.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