*HEADING
 Same as iwr_1d_1_xpl_p_pp.inp except that the total wave formulation 
 is specified.  However, the solution should not change as there are 
 no acoustic elements in this problem. 
*** 
 1D transient wave propagation, solid medium, incident wave 
 load and incident wave reflection, spherical and defined through 
 pressure amplitude
**
*ACOUSTIC WAVE FORMULATION, TYPE=TOTAL WAVE
*PREPRINT,ECHO=NO,HISTORY=yes,MODEL=YES
*NODE,NSET=PLATE
 1,-0.5,-0.5
 2,0.5,-0.5
 3,0.5,0.5
 4,-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=PLANE
0.0,0.0,0.
0.0,10.0,10.
*INCIDENTWAVE FLUID PROPERTY
0.16e6,1.0
*AMPLITUDE,DEFINITION=periodic,NAME=PRESSUREVTIMESINE
1,350.0,0.0,0.0
0.0,1.0
*AMPLITUDE,DEFINITION=tabular,NAME=PRESSUREVTIMESTEP
0.0,0.0,0.01,1.0
**************************
*BOUNDARY
PLATE,ENCASTRE
**************************
*STEP
*DYNAMIC,EXPLICIT,DIRECTUSERCONTROL
0.0001,0.06
**
*INCIDENTWAVE,PRESSUREAMPLITUDE=PRESSUREVTIMESTEP,PROPERTY=IWPROP
IW_LOAD,1.E3
**
*INCIDENT WAVE REFLECTION
20.0,0,-1,0, -1.0
**
*OUTPUT,FIELD,NUMBERINTERVAL=1
*NODE OUTPUT,NSET=PLATE
RF,
*OUTPUT,HISTORY,FREQUENCY=1
*NODE OUTPUT,NSET=PLATE
RF,
**
*END STEP
*STEP
*DYNAMIC,EXPLICIT,DIRECTUSERCONTROL
0.0001,0.09
**
*INCIDENTWAVE,PRESSUREAMPLITUDE=PRESSUREVTIMESINE,PROPERTY=IWPROP
IW_LOAD,1.E3
**
*INCIDENT WAVE REFLECTION
20.0,0,-1,0, -1.0
**
*OUTPUT,FIELD,NUMBERINTERVAL=1
*NODE OUTPUT,NSET=PLATE
RF,
*OUTPUT,HISTORY,FREQUENCY=1
*NODE OUTPUT,NSET=PLATE
RF,
**
*END STEP