*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
 401, 0, 0.5E-2, -0.5E-2
 411, .1, 0.5E-2, -0.5E-2
*NGEN,NSET=CHANNEL
   1,11
  51,61
 101,111
 151,161
 201,211
 251,261
 301,311
 351,361
 401,411
*NSET,NSET=SAMPLE,GENERATE
 1,11
*ELEMENT,TYPE=C3D10M
 1,    1, 251, 101,    3, 351, 401,  51,
            2, 352,  52,
 2, 251, 101,    3, 253, 401,  52, 352,
         252, 402, 353,
 3, 101,    3, 253, 103,  52, 353, 402,
         102,  53, 403,
 4, 101, 251, 151, 103, 401, 301, 201,
         102, 402, 202,
 5, 251, 151, 103, 253, 301, 202, 402,
         252, 302, 403,
 6, 151, 103, 253, 153, 202, 403, 302,
         152, 203, 303,
*ELGEN,ELSET=CHANNEL
 1,5,2,6
 2,5,2,6
 3,5,2,6
 4,5,2,6
 5,5,2,6
 6,5,2,6
**************************
*SURFACE,NAME=IW_LOAD
4,S1
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=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
411,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