*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=PLATE
 1, 0. , 0. ,
 2, 1. , 0. ,
 3, 1. , 1. ,
 4, 0. , 1. ,
 5, 0.5, 0. ,
 6, 1.0, 0.5,
 7, 0.5, 1.0,
 8, 0. , 0.5,
 9, 0.5, 0.5,
*ELEMENT,TYPE=S8R,ELSET=PLATE
 1,1,2,3,4,5,6,7,8 
**************************
*SURFACE,NAME=IW_LOAD
PLATE,SPOS
**************************
*SHELL SECTION,MATERIAL=SOLID,ELSET=PLATE
 1.E-4, 
**
*MATERIAL,NAME=SOLID
*DENSITY
 1.E3,
*ELASTIC
 1.E6,0.0
**
**************************
*INCIDENT WAVE PROPERTY,NAME=IWPROP,TYPE=SPHERE
0.5,0.5,0.0
0.5,0.5,10.0E04
*INCIDENTWAVE FLUID PROPERTY
0.1424E6,1.21,0.0
*AMPLITUDE,DEFINITION=TABULAR,NAME=PRESSUREVTIME
0.0,0.0,0.0005,1.0
**************************
*BOUNDARY
PLATE,ENCASTRE
**
**************************
*RESTART,WRITE
**
*STEP,INC=1000,NLGEOM
*DYNAMIC,ALPHA=0.0
.0001,.0005
**
*INCIDENTWAVE,PRESSUREAMPLITUDE=PRESSUREVTIME,PROPERTY=IWPROP
IW_LOAD,1.E3
**
*NODE PRINT,NSET=PLATE,FREQ=22
 RF,
*OUTPUT,VAR=PRESELECT,FIELD,FREQ=1
*OUTPUT,VAR=PRESELECT,HISTORY,FREQ=1
*OUTPUT,FIELD, FREQUENCY=1
*NODE OUTPUT
 RF,
*OUTPUT,HISTORY,FREQUENCY=1
*NODE OUTPUT,NSET=PLATE
 RF,
*ENDSTEP
**
**