*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,
  51, 0,.5E-2
  61, .1,.5E-2
 101, 0, 1.E-2
 111, .1, 1.E-2
 151, 0, 1.E-2, -1.E-2
 161, .1, 1.E-2, -1.E-2
 201, 0, 1.E-2, -0.5E-2
 211, .1, 1.E-2, -0.5E-2
 251, 0, 0, -1.E-2
 261, .1, 0, -1.E-2
 301, 0,.5E-2, -1.E-2
 311, .1,.5E-2, -1.E-2
 351, 0, 0, -0.5E-2
 361, .1, 0, -0.5E-2
*NGEN,NSET=CHANNEL
   1,11
  51,61
 101,111
 151,161
 201,211
 251,261
 301,311
 351,361
*NSET,NSET=SAMPLE,GENERATE
 1,11
*ELEMENT,TYPE=C3D20
 1,1,3,253,251,101,103,153,151,
 2,353,252,351,102,203,152,201,
 51,53,303,301
*ELGEN,ELSET=CHANNEL
 1,5,2
**************************
*SURFACE,NAME=IW_LOAD
1,S6
**************************
*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
 61,1,1
111,1,1
161,1,1
211,1,1
261,1,1
311,1,1
361,1,1
**************************
*RESTART,WRITE
**
*STEP,INC=1000,NLGEOM
*DYNAMIC,ALPHA=0.0
.0001,.0015
**
*INCIDENT WAVE,PRESSUREAMPLITUDE=PRESSUREVTIME,PROPERTY=IWPROP
IW_LOAD,415.09517
**
*EL PRINT,FREQ=5
 S,
*NODE PRINT,NSET=SAMPLE,FREQ=5
 U,
*OUTPUT,VAR=PRESELECT,FIELD,FREQ=1
*OUTPUT,VAR=PRESELECT,HISTORY,FREQ=1
*OUTPUT,FIELD,FREQUENCY=1
*NODE OUTPUT
 U
*ELEMENT OUTPUT
 S,
*OUTPUT,HISTORY,FREQUENCY=1
*NODE OUTPUT,NSET=SAMPLE
 U,
*ENDSTEP