*HEADING
 1D transient wave propagation, solid  medium, incident wave load, planar
 and defined through pressure amplitude
**
** C = root (K/rho)= 343.05 m/s   K=bulk mod, rho=density.
**
** 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=BEAM
1,0. ,
2,0.5,
3,1.0,
*ELEMENT,TYPE=B23,ELSET=BEAM
1,1,3
**************************
*SURFACE, NAME=IW_LOAD
BEAM,SPOS
**************************
*BEAM SECTION,SECTION=RECT,TEMPERATURE=GRADIENTS,MATERIAL=SOLID,ELSET=BEAM
1.E-2, 1.E-2
0., 0.,  -1.
3,
*MATERIAL,NAME=SOLID
*DENSITY
 1.E3
*ELASTIC
 1.E6,0.0
**
**************************
*INCIDENT WAVE PROPERTY,NAME=IWPROP,TYPE=PLANE
0.5,0.0
0.5,10.0
*INCIDENTWAVE FLUID PROPERTY
0.1424E6,1.21,0.0
*AMPLITUDE,DEFINITION=TABULAR,NAME=PRESSUREVTIME
0.0,0.0,0.0005,1.0
**************************
*BOUNDARY
BEAM,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=BEAM,FREQ=1
 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=BEAM
 RF,
*ENDSTEP
**
**