*HEADING
 EC3FAFAV: AC3D15 ELEMENTS
 1D transient wave propagation in an acoustic medium with CLOAD on dof 8. AC3D15
 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
 1501, 0, 1.e-2, -1.e-2
 1601, 1, 1.e-2, -1.e-2
 2001, 0, 1.e-2, -0.5e-2
 2101, 1, 1.e-2, -0.5e-2
 2501, 0, 0, -1.e-2
 2601, 1, 0, -1.e-2
 3001, 0,.5e-2, -1.e-2
 3101, 1,.5e-2, -1.e-2
 3501, 0, 0, -0.5e-2
 3601, 1, 0, -0.5e-2
 4001, 0, 0.5e-2, -0.5e-2
 4101, 1, 0.5e-2, -0.5e-2
*NGEN,NSET=CHANNEL
    1,101
  501,601
 1001,1101
 1501,1601
 2001,2101
 2501,2601
 3001,3101
 3501,3601
 4001,4101
*NSET,NSET=SAMPLE,GENERATE
 1,101
*NSET,NSET=CRNRLOAD
 1,1001,1501,2501
*NSET,NSET=EDGELOAD
 501,2001,3001,3501
*ELEMENT,TYPE=AC3D15
    1,    1, 1001,    3, 2501, 1501, 2503,
          501,  502,    2, 3001, 3002, 2502,
         3501, 2001, 3503,
    2,    3, 1001, 1003, 2503, 1501, 1503,
          502, 1002,  503, 3002, 1502, 3003,
         3503, 2001, 2003,
*ELGEN,ELSET=CHANNEL
 1,50,2,2
 2,50,2,2
*SOLID SECTION,MATERIAL=AIR,ELSET=CHANNEL
** Cross-section is 3.0  
**
*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, 
**
**
*STEP,INC=1000
*DYNAMIC
.0001,.0044
** Define nonreflective boundary at x=1
*IMPEDANCE
 100,I5
**
** Prescribe volume acceleration inwards
**  of  a0 * cross-section= 3., therefore a0=1.
**  Distribute (-1/12,1/3) to each (corner,
**  midside) node, respectively:
*CLOAD
 CRNRLOAD,8,-.08333333333333e-4
 EDGELOAD,8, .33333333333333e-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