*HEADING
GLOBAL MODEL [AC3D4]
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=AC3D4,ELSET=EALL
    1,    1, 1001, 1003, 2001
    2, 1001, 1003, 2001, 3001
    3, 1003, 2001, 3001, 3003
    4,    1, 1003,    3, 2001
    5, 1003,    3, 2001, 3003
    6,    3, 2001, 3003, 2003
*ELGEN,ELSET=EALL
1,20,2,6
2,20,2,6
3,20,2,6
4,20,2,6
5,20,2,6
6,20,2,6
*MATERIAL,NAME=AIR
*DENSITY
1.293, 
*ACOUSTIC MEDIUM,BULK MODULUS
1.42176E5, 
*SOLID SECTION,MATERIAL=AIR,ELSET=EALL
*NSET,NSET=RECRD,GEN
1,41,2
1001,1041,2
2001,2041,2
3001,3041,2
*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
120,I4
119,I2
*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=10
*NODE OUTPUT
POR,
*******************
*OUTPUT,HISTORY,TIME INTERVAL=<out_hist>
*NODE OUTPUT,NSET=ODB_NSET
POR,
*******************
*END STEP