*HEADING
EC26AFAW: Nonreflective end, complex material via volumetric
drag and direct 
Compare IMPEDANCE WITH SIMPEDANCE, and STEADY STATE DYNAMICS, DIRECT
with SUBSPACE PROJECTION.
********************
*part,name=duct1
*NODE, Nset=ODB_NSET
1,    0
41,  .5
101,  0,1.25e-2
141, .5,1.25e-2 
51,   0,6.25e-3 
91,  .5,6.25e-3 
*NGEN,NSET=NALL
1,41
51,91
101,141
*NSET,NSET=LOAD
1,51,101
*NSET,NSET=QA_TEST,GENERATE
1,41,2
*ELEMENT,TYPE=AC2D6,ELSET=EALL
1,1,3,103,2,53,52
2,1,103,101,52,102,51
*ELGEN,ELSET=EALL
1,20,2,2
2,20,2,2
*ELSET,ELSET=NONREF
20,
*SOLID SECTION,MATERIAL=ABSORB,ELSET=EALL
1.,
*end part
**
*part,name=duct2
*NODE, Nset=ODB_NSET
1,    0
41,  .5
101,  0,1.25e-2
141, .5,1.25e-2 
51,   0,6.25e-3 
91,  .5,6.25e-3 
*NGEN,NSET=NALL
1,41
51,91
101,141
*NSET,NSET=LOAD
1,51,101
*NSET,NSET=QA_TEST,GENERATE
1,41,2
*ELEMENT,TYPE=AC2D6,ELSET=EALL
1,1,3,103,2,53,52
2,1,103,101,52,102,51
*ELGEN,ELSET=EALL
1,20,2,2
2,20,2,2
*ELSET,ELSET=NONREF
20,
*SOLID SECTION,MATERIAL=complex,ELSET=EALL
1.,
*end part
**
*assembly,name=AcousticMaterialTest
*instance,part=duct1,name=duct1
*end instance
*instance,part=duct2,name=duct2
0,0.1,0
*end instance
*Elset,Elset=ODB_ESET
DUCT1.EALL,DUCT2.EALL
*Nset, Nset=QA_TEST
DUCT1.qa_test,DUCT2.QA_TEST
*NSET,NSET=LOAD
DUCT1.LOAD,DUCT2.LOAD
*ELSET,ELSET=NONREF
DUCT1.NONREF,DUCT2.NONREF
*surface,NAME=NRS
DUCT1.NONREF,S2
DUCT2.NONREF,S2
*Nset,Nset=corners
duct1.101,duct1.1,duct2.101,duct2.1
*nset,nset=mid
duct1.51,duct2.51
*end assembly
**
*MATERIAL,NAME=ABSORB
*DENSITY
 1.21,
*ACOUSTIC MEDIUM,BULK MODULUS
 0.1424E6,
*ACOUSTIC MEDIUM,VOLUMETRIC DRAG
1400.,
**
*MATERIAL,NAME=complex
*ACOUSTIC MEDIUM,COMPLEX DENSITY
 1.21, -2.2282,100.0
*ACOUSTIC MEDIUM,BULK MODULUS
 0.1424E6,
**
*AMPLITUDE,NAME=OMEGA
0.,0., 1E6,6.2831853E6
************************************************
**       ** ANALYTIC SOLUTION  **
**  
**  at freq=100,
**  Impedance: Z = root( K_f (rho  - i r/omega))
**               = 516.411 -307.208 i
**  Propagation constant 
**             b = i (omega/K_f) Z
**               = 1.35551 + 2.27859 i
**
**  Solution: with nonreflective end,
**  Complex pressure
**             p(x) = Z v(0) exp( - b x)
**  i.e. |p(x)| = |Z v(0)| exp ( - Re(b) x)
**              = 600.8805  exp ( -1.35551 x)
**
**       Arg(p) = Arg(Z) - Im(b) x
**       ppor   = -30.74799 - 130.5534 x
**
************************************************  
*physi,spl=1.23456
**
*STEP
*STEADY STATE DYNAMIC,DIRECT
100,
*IMPEDANCE
NONREF,I2
**
** Conditions at x=0 : Define the source inward accel
**  to be  a =  omega i, so v (inward) = 1.
**
**  Area of face=A. Total inward VOLUME acceleration  
**  is A*a.  This should be distributed in the ratio  
**  (1/6, 2/3, 1/6) over (corner, midedge, corner)
**  nodes respectively:
*CLOAD,LOAD CASE=2, AMP=OMEGA
corners,  8, 2.08333333e-3
mid, 8, 8.33333333e-3
*ENERGYPRINT
*NODE PRINT, NSET=qa_test
COORD,POR,PPOR
*OUTPUT,FIELD
*NODE OUTPUT
POR,SPL
*ELEMENTOUTPUT
ENER,ELEN,ELEDEN
*OUTPUT,HISTORY
*NODE OUTPUT,NSET=QA_TEST
POR,
*ENERGYOUTPUT
ALLQB,ALLWK,ALLVD,ALLSE,ALLKE
*END STEP
**
**
*STEP
*STEADY STATE DYNAMIC,DIRECT
100,
*SIMPEDANCE
NRS,
*CLOAD,LOAD CASE=2, AMP=OMEGA
corners,  8, 2.08333333e-3
mid, 8, 8.33333333e-3
*ENERGYPRINT
*NODE PRINT, NSET=QA_TEST
COORD,POR,PPOR
*OUTPUT,FIELD
*NODE OUTPUT
POR,SPL
*ELEMENTOUTPUT
ENER,ELEN,ELEDEN
*OUTPUT,HISTORY
*NODE OUTPUT,NSET=QA_TEST
POR,
*ENERGYOUTPUT
ALLQB,ALLWK,ALLVD,ALLSE,ALLKE
*END STEP
**
**
*STEP
*FREQUENCY,NORM=MASS,EIGENSOLVER=LANCZOS, SIM=NO
22, , ,
*NODE PRINT,F=0
*EL PRINT,F=0
*OUTPUT,FIELD,Frequency=0
*END STEP
**
**
*STEP
*STEADY STATE DYNAMIC,SUBSPACE PROJECTION
100,1000,20
*IMPEDANCE
NONREF,I2
*CLOAD,LOAD CASE=2, AMP=OMEGA
corners,  8, 2.08333333e-3
mid, 8, 8.33333333e-3
*ENERGYPRINT
*OUTPUT,FIELD
*NODE OUTPUT
POR,SPL
*ELEMENTOUTPUT
ENER,ELEN,ELEDEN
*OUTPUT,HISTORY
*ENERGYOUTPUT
ALLQB,ALLWK,ALLVD,ALLSE,ALLKE
*END STEP
**
*STEP
*STEADY STATE DYNAMIC,SUBSPACE PROJECTION
100,
*SIMPEDANCE
NRS, 
*CLOAD,LOAD CASE=2, AMP=OMEGA
corners,  8, 2.08333333e-3
mid, 8, 8.33333333e-3
*ENERGYPRINT
*OUTPUT,FIELD
*NODE OUTPUT
POR,SPL
*ELEMENTOUTPUT
ENER,ELEN,ELEDEN
*OUTPUT,HISTORY
*NODE OUTPUT,NSET=QA_TEST
POR,
*ENERGYOUTPUT
ALLQB,ALLWK,ALLVD,ALLSE,ALLKE
*END STEP