*HEADING
--
 Same as iw_1d_c3d8r_xpl_s_pp.inp but now the total wave formulation 
 option is chosen.  Note however that this option is meaningful only
 when the wave propagation in fluids is modeled as well. Hence,  
 the presence of the total wave formulation option should not 
 affect the current analysis.  This input deck checks to make
 sure that load application on solid surfaces is not affected. 
--
 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
** 
*ACOUSTIC WAVE FORMULATION, TYPE=TOTAL
*PREPRINT,ECHO=NO,HISTORY=NO,MODEL=YES
*NODE, NSET=ODB_NSET
   1, 0,
  11, .1,
 101, 0, 1E-2
 111, .1, 1E-2
 151, 0, 1E-2, -1E-2
 161, .1, 1E-2, -1E-2
 251, 0, 0, -1E-2
 261, .1, 0, -1E-2
*NGEN,NSET=CHANNEL
   1,11
 101,111
 151,161
 251,261
*NSET,NSET=SAMPLE,GENERATE
 1,11
*ELEMENT,TYPE=C3D8R
 1,1,2,252,251,101,102,152,151
*ELGEN,ELSET=CHANNEL
 1,10
**************************
*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=SPHERE
0.0,0.0,0.0
-10.E04,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
111,1,1
161,1,1
261,1,1
**************************
*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