*HEADING
 EC35AFAD_AMS: AC3D5 ELEMENTS
 Acoustic test - 1D Problem, SS-DIRECT.
 A 4m long tube; from left to right,
 a sound source, 3m is air, then 1m of a highly
 absorbent material, followed by a rigid
 termination. 
*PREPRINT,ECHO=NO,HISTORY=NO,MODEL=NO
*RESTART,WRITE
***********************************************************
**
**
**   Mesh Generation
**
*NODE, Nset=ODB_NSET
    1, 0, 0,  0
  401, 4, 0,  0
  501, 0, 0, -2e-2
  901, 4, 0, -2e-2
 1001, 0, 2e-2,  0
 1401, 4, 2e-2,  0
 1501, 0, 2e-2, -2e-2
 1901, 4, 2e-2, -2e-2
*NGEN,NSET=ALL
    1,401
  501,901
 1001,1401
 1501,1901
*NSET,NSET=QA_TEST,GENERATE
 1,401
*NSET,NSET=LOAD
 1,1001,1501,501
*ELEMENT,TYPE=AC3D5
   1,1,501,1501,1001,2
1001,1001,1501,1502,1002,2
2001,1501,501,502,1502,2
   2,3,1003,1503,503,2
1002,1002,1502,1503,1003,2
2002,1503,1502,502,503,2
*ELGEN,ELSET=CHANNEL
   1,150,2,2
1001,150,2,2
2001,150,2,2
   2,150,2,2
1002,150,2,2
2002,150,2,2
*ELEMENT,TYPE=AC3D5
 301,301,801,1801,1301,302
1301,1301,1801,1802,1302,302
2301,1801,801,802,1802,302
 302,303,1303,1803,803,302
1302,1303,1302,1802,1803,302
2302,1803,1802,802,803,302
*ELGEN,ELSET=ABSORB
  301,50,2,2
 1301,50,2,2
 2301,50,2,2
  302,50,2,2
 1302,50,2,2
 2302,50,2,2
**
**
************************************************
**    ** ELEMENT PROPERTIES **
**
*SOLID SECTION,MATERIAL=AIR,ELSET=CHANNEL
*SOLID SECTION,MATERIAL=ABSORB,ELSET=ABSORB
**
**
*MATERIAL,NAME=AIR
*DENSITY
 1.21,
*ACOUSTIC MEDIUM,BULK MODULUS
 0.1424E6,
**
**
*MATERIAL,NAME=ABSORB
*DENSITY
 1.21,
*ACOUSTIC MEDIUM,BULK MODULUS
 0.1424E6,
*ACOUSTIC MEDIUM,VOLUMETRIC DRAG
 1e4,
**
**
**
*************************************************
**
*Elset,Elset=ODB_ESET
CHANNEL, ABSORB
**
**
*Nset, Nset=ODB_NSET
ODB_NSET, ALL, QA_TEST
**
**
*STEP
*STEADY STATE DYNAMICS,DIRECT,FREQUENCY SCALE=LINEAR
 10,1100,13,1
**
** Prescribe outward acceleration of 100 at node LOAD.
** The tributary area of each node is 1, so the 
** outward volume acceleration is 100.
*CLOAD
 LOAD,8, -1E2
**
*NODE PRINT, NSET=QA_TEST
 POR,PPOR,
*EL PRINT,FREQUENCY=0
*ENERGYPRINT
*OUTPUT,FIELD
*NODE OUTPUT,NSET=QA_TEST
 POR,
*ELEMENT OUTPUT
ACV,INTEN,GRADP,ENER,ELEN,ELEDEN
*OUTPUT,HISTORY
*NODE OUTPUT,NSET=QA_TEST
 POR,
*ENERGY OUTPUT
ALLSE,ALLKE,ALLWK,ALLVD
*ENDSTEP
**
**
*STEP
***FREQUENCY,NORM=MASS,EIGENSOLVER=LANCZOS
*FREQUENCY,EIGENSOLVER=AMS
200,,300.
*OUTPUT,FIELD,FREQUENCY=999
*NODE OUTPUT
U,
*ELEMENT OUTPUT
S,
**50,,,50
*NODE PRINT,F=0
*EL PRINT,F=0
*OUTPUT,FIELD,Frequency=0
*END STEP
**
**
*STEP
*STEADY STATE DYNAMICS,FREQUENCY SCALE=LINEAR,INTERVAL=RANGE
 10,1100,13,1
**
** Prescribe outward acceleration of 100 at node LOAD.
** The tributary area of each node is 1, so the 
** outward volume acceleration is 100.
*SELECT EIGENMODES,GENERATE
1,47
*CLOAD
 LOAD,8, -1E2
**
*NODE PRINT, NSET=QA_TEST
 POR,PPOR,
*EL PRINT,FREQUENCY=0
*energyprint
*OUTPUT,FIELD,VARIABLE=PRESELECT,OP=ADD,FREQUENCY=9999
*OUTPUT,FIELD, OP=ADD
*NODE OUTPUT,NSET=QA_TEST
 POR,
*ELEMENT OUTPUT,ELSET=ABSORB
ACV,INTEN,GRADP
*OUTPUT,HISTORY
*NODE OUTPUT,NSET=QA_TEST
 POR,
***ENERGY OUTPUT
**ALLSE,ALLKE,ALLWK,ALLVD
*ENDSTEP