*HEADING
GLOBAL MODEL [AC3D6]
ABAQUS/Explicit
**********************************************
** The speed of sound c=sqrt(K/rho)=330m/s  **
** The pressure wave will reach the end of  **
** tube (node 41) in t=L_tube/c=0.0301 s    **
** Plot P_1 and P_41 to check the results   **
**********************************************
*PARAMETER
f=25.0
T=1.0/f
out_hist=T/100
omega=2*3.14159*f
******************
*NODE, NSET=ODB_NSET
1, 
41,,10.
1001,1.
1041,1.,10.
2001,,,1.
2041,,10.,1.
3001,1.,,1.
3041,1.,10.,1.
*NSET,NSET=QA_TEST
ODB_NSET,
*NGEN,NSET=NALL
1,41
1001,1041
2001,2041
3001,3041
*ELEMENT,TYPE=AC3D6,ELSET=EALL
1,1,1001,1003,2001,3001,3003
2,1,1003,3,2001,3003,2003
*ELGEN,ELSET=EALL
1,20,2,2
2,20,2,2
*ELSET,ELSET=LINER5
40, 
*MATERIAL,NAME=AIR
*DENSITY
1.293, 
*ACOUSTIC MEDIUM,BULK MODULUS
1.42176E5, 
*SOLID SECTION,MATERIAL=AIR,ELSET=EALL
*IMPEDANCE PROPERTY,NAME=IMPED
0.0,2.3323E-3
*NSET,NSET=RECRD,GEN
1,41,4
1001,1041,4
2001,2041,4
3001,3041,4
*NSET,NSET=DRIV
1,1001,2001,3001
*AMPLITUDE,NAME=SINE,DEF=PERIODIC
1,<omega>,0.0,0.0
0,1,0,0
*STEP
*DYNAMIC,EXPLICIT
,<T>
*IMPEDANCE
LINER5,I4,IMPED
*BOUNDARY,AMP=SINE
DRIV,8,,1.
*********************
*FILE OUTPUT,NUM=10
*NODE FILE
POR,
*********************
*OUTPUT,FIELD,NUM=1
*NODE OUTPUT,NSET=QA_TEST
POR,
*******************
*OUTPUT,FIELD,NUM=4
*NODE OUTPUT
POR,
*******************
*OUTPUT,HISTORY,TIME INTERVAL=<out_hist>
*NODE OUTPUT,NSET=ODB_NSET
POR,
*******************
*END STEP