*HEADING
 EC36AFAV: wave propagation with CLOAD on dof 8. AC3D6
 Verification problem: Transient Dynamics Tests
**
** C = root (K/rho)= 343.05 m/s   K=bulk mod, rho=density.
**
** Analytic 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=no
*NODE, Nset=ODB_NSET
*NGEN,NSET=CHANNEL_NODES
*NSET,NSET=SAMPLE_NODES,GENERATE
*NSET,NSET=LOADED_NODES
*ELEMENT,TYPE=AC3D6
*ELGEN,ELSET=CHANNEL
*SOLID SECTION,MATERIAL=AIR,ELSET=CHANNEL
** Cross-section is 3.0
**
*MATERIAL,NAME=AIR
*DENSITY
 1.21,
*ACOUSTIC MEDIUM,BULK MODULUS
**
** Define nonreflective boundary at x=1
** Using user calculated constants:
** 1/k1=0, 1/c1=1/(rho*c)
*IMPEDANCE PROPERTY,NAME=NONREF
**
*RESTART,WRITE
**
**
** Sets added by addodb
**
*Elset,Elset=ODB_ESET
CHANNEL,
**
**
*Nset, Nset=ODB_NSET
ODB_NSET, CHANNEL_NODES, SAMPLE_NODES, LOADED_NODES
**
**
*STEP,INC=1000
*DYNAMIC
.0001,.0044
*IMPEDANCE
 199,I4,NONREF
**
** Prescribe  acceleration inwards of  a0=1.
** volume acc = a0 * cross - sectional area = 1e-4
**  Distribute 1/4 to each corner node:
*CLOAD
 LOADED_NODES,8,.25e-4
**
*EL PRINT,FREQ=0
*NODE PRINT,NSET=SAMPLE_NODES,FREQ=22
 POR,
*NODE FILE,NSET=SAMPLE_NODES,FREQ=22
 POR,
*OUTPUT,FIELD,VARIABLE=PRESELECT,OP=ADD,FREQUENCY=9999
*OUTPUT,FIELD, OP=ADD,Frequency=22
*NODE OUTPUT,NSET=SAMPLE_NODES
 POR,
*OUTPUT,HISTORY,OP=ADD,Frequency=22
*NODE OUTPUT,NSET=SAMPLE_NODES
 POR,
*ENDSTEP