*HEADING
--
 Same as iw_1d_s4r_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=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=S4R,ELSET=PLATE
 1,1,2,3,4
**************************
*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
*DYNAMIC,EXPLICIT, DIRECT USER
0.2E-4,.0005,
**
*INCIDENTWAVE,PRESSUREAMPLITUDE=PRESSUREVTIME,PROPERTY=IWPROP
IW_LOAD,1.E3
**
*OUTPUT,FIELD,NUMBER INTERVAL=2
*NODE OUTPUT,NSET=PLATE
 RF,
*OUTPUT,HISTORY,TIME INTERVAL=0.1E-3
*NODE OUTPUT,NSET=PLATE
 RF,
*ENDSTEP
**
**