*HEADING
 1D transient wave propagation in an acoustic medium with CLOAD on dof 8. AC3D10
 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=AC3D10
    1,    1, 2501, 1001,    3, 3501, 4001,  501,
            2, 3502,  502,
    2, 2501, 1001,    3, 2503, 4001,  502, 3502,
         2502, 4002, 3503,
    3, 1001,    3, 2503, 1003,  502, 3503, 4002,
         1002,  503, 4003,
    4, 1001, 2501, 1501, 1003, 4001, 3001, 2001,
         1002, 4002, 2002,
    5, 2501, 1501, 1003, 2503, 3001, 2002, 4002,
         2502, 3002, 4003,
    6, 1501, 1003, 2503, 1503, 2002, 4003, 3002,
         1502, 2003, 3003,
*ELGEN,ELSET=CHANNEL
 1,50,2,6
 2,50,2,6
 3,50,2,6
 4,50,2,6
 5,50,2,6
 6,50,2,6
*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
 300,I3
 297,I3
**
**  Prescribe unit _particle_ acceleration a0 inwards;
**  _volume_ acceleration ("force" for acoustics) is a0*cross-sec.
**  Distribute (1/6,-1/6, 1) to each (corner,midside,center) 
**  node, respectively (appropriate for triangles):
*CLOAD
 CRNRLOAD,8,0.16666666666e-4
 EDGELOAD,8,-.16666666666e-4
 4001,8,1.0e-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