*HEADING
 ECA64AFAT: ACAX6 AND CAX8 ELEMENTS
 Tube with nonreflective end-condition and *TIE constraint.
 Part of the air in the tube is modelled with continuum elements
 with air properties.
************************************************
*NODE, Nset=ODB_NSET
1,1000.0,
41,1000.5,
101,1000.0,0.0125
141,1000.5,0.0125
201,1000.0,0.025
241,1000.5,0.025
*NGEN,NSET=NALL
1,41
101,141
201,241
*NCOPY,OLDSET=NALL,NEWSET=NSOL,CHANGENUMBER=1000,SHIFT
0.0,0.0,0.0
1.0,0.0,0.0,0.0,1.0,0.0
*NSET,NSET=SAMPLE,GENERATE
1,41
*NSET,NSET=LOAD1
1,201
*NSET,NSET=LOAD2
101, 
** First 10 elements are acoustic
*ELEMENT,TYPE=ACAX6
1,1,3,203,2,103,102
2,203,201,1,202,101,102
*ELGEN,ELSET=ACOU
1,5,2,2
2,5,2,2
** Then 5 continuum elements
*ELEMENT,TYPE=CAX4
1006,1011,1013,1213,1211
*ELGEN,ELSET=STRUCT
1006,5,2,1
** Next 20 elements are acoustic
*ELEMENT,TYPE=ACAX6
31,21,23,223,22,123,122
32,223,221,21,222,121,122
*ELGEN,ELSET=ACOU
31,10,2,2
32,10,2,2
*ELSET,ELSET=NONREF
49,
************************************************
*SURFACE,TYPE=ELEMENT,NAME=ACR
9,S2
*SURFACE,TYPE=ELEMENT,NAME=SOL
1006,S4
*SURFACE,TYPE=ELEMENT,NAME=SOR
1010,S2
*SURFACE,TYPE=ELEMENT,NAME=ACL
32,S2
*TIE,NAME=CRAVAT
ACR,SOL
ACL,SOR
************************************************
**
**  Acoustic Element Properties
**
*SOLID SECTION,MATERIAL=AIR,ELSET=ACOU
*MATERIAL,NAME=AIR
*DENSITY
 1.21,
*ACOUSTIC MEDIUM,BULK MODULUS
 0.1424E6,
************************************************
**
**  Continuum Element Properties
**
*SOLID SECTION,MATERIAL=AIRSOLID,ELSET=STRUCT
*MATERIAL,NAME=AIRSOLID
*DENSITY
 1.21,
*ELASTIC
 0.1424E6,
************************************************
*AMPLITUDE,NAME=OMEGA
0.,0., 1E6,6.2831853E6
************************************************
**       ** ANALYTIC SOLUTION  **
**  
**  at freq=100,
**  Impedance: Z = root( K_f rho )
**               = 415.095
**  Propagation constant 
**             b = i (omega/K_f) Z
**               = 1.8315  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)
**              =  415.095
**
**       Arg(p) = Arg(Z) - Im(b) x
**       ppor   = - 104.94 x
**
************************************************  
*Elset,Elset=ODB_ESET
ACOU, STRUCT
*Nset, Nset=ODB_NSET
ODB_NSET, NALL, SAMPLE, LOAD1
LOAD2, 
************************************************  
*STEP
*STEADY STATE DYNAMIC,DIRECT
100,
**  For reference, these IMPEDANCE PROPERTY
**  constants could be used to specify the nonreflecting 
**  condition at the given frequency:
**  1/k1= 0, 1/c1= 2.40908E-03
**        0,       1/Z
**
*IMPEDANCE
 NONREF,I2
**
** Conditions at x=0 : Define the source inward accel
** to be  a =  omega i, so v (inward) = 1.
**
*CLOAD,LOAD CASE=2, AMP=OMEGA
LOAD1, 8, 26.179939
LOAD2, 8, 104.71976
** simulation of uniaxial condition
*BOUNDARY
NSOL,2,2
*NODE PRINT, NSET=SAMPLE
COORD,POR,PPOR
*OUTPUT,FIELD,OP=ADD,FREQUENCY=9999
*ELEMENT OUTPUT,ELSET=STRUCT
S
*NODE OUTPUT,NSET=SAMPLE
POR,
*OUTPUT,HISTORY,OP=ADD,FREQUENCY=9999
*NODE OUTPUT,NSET=SAMPLE
POR,
*ELEMENT OUTPUT,ELSET=STRUCT
COORD,S11
*END STEP