*HEADING
 EC26AFAV: AC2D6
 1D transient wave propagation in an acoustic medium with CLOAD DOF 8. 
 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
    1, 0,
  101, 1,
  501, 0,.5e-2
  601, 1,.5e-2
 1001, 0, 1.e-2
 1101, 1, 1.e-2
*NGEN,NSET=CHANNEL
    1,101
  501,601
 1001,1101
*NSET,NSET=SAMPLE,GENERATE
 1,101
*NSET,NSET=CRNRLOAD
 1,1001
*NSET,NSET=EDGELOAD
 501, 
*NSET,NSET=LOAD
 CRNRLOAD,EDGELOAD
*ELEMENT,TYPE=AC2D6
 1,1,3,1003,2,503 , 502
 2,1003,1001,1,1002,501,502
*ELGEN,ELSET=CHANNEL
 1,50,2,2
 2,50,2,2
*SOLID SECTION,MATERIAL=AIR,ELSET=CHANNEL
 1.e-2, 
** Cross-section is 1e-4
**
*MATERIAL,NAME=AIR
*DENSITY
 1.21,
*ACOUSTIC MEDIUM,BULK MODULUS
 0.1424E6,
**
**
*RESTART,WRITE
**
**
** Sets added by addodb
**
*Elset,Elset=ODB_ESET
CHANNEL,
**
**
*Nset, Nset=ODB_NSET
ODB_NSET, CHANNEL, SAMPLE, CRNRLOAD
EDGELOAD, LOAD
**
**
*STEP,INC=1000
*DYNAMIC
.0001,.0044,,0.0001
** Define nonreflective boundary at x=1
*IMPEDANCE
 99,I2
**
** Prescribe volume acceleration inwards
**  of a0=1.
**  Distribute (1/6,2/3,1/6) to (corner,midside,
**  corner) nodes, respectively:
*CLOAD
 CRNRLOAD,8,.16666666e-4
 EDGELOAD,8,.66666666e-4
**
*EL PRINT,FREQ=0
*NODE PRINT,NSET=SAMPLE,FREQ=22
 POR,
*NODE FILE,NSET=SAMPLE,FREQ=22
 POR,
*OUTPUT,FIELD,VARIABLE=PRESELECT,OP=ADD,FREQUENCY=9999
*OUTPUT,FIELD, OP=ADD,Frequency=22
*NODE OUTPUT,NSET=SAMPLE
 POR,
*OUTPUT,HISTORY,OP=ADD,Frequency=22
*NODE OUTPUT,NSET=SAMPLE
 POR,
*ENDSTEP