*HEADING
:SUBMODEL [AC2D3]--Global Analysis
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
f_out_hist=T/100
omega=2*3.14159*f
******************
*NODE, Nset=ODB_NSET
1, 
41,,10.
1001,1.
1041,1.,10.
*NSET,NSET=QA_TEST
41,1041
*NGEN,NSET=NALL
1,41
1001,1041
*ELEMENT,TYPE=AC2D3,ELSET=EALL
1,1,1001,1003
2,1,1003,3
*ELGEN,ELSET=EALL
1,20,2,2
2,20,2,2
*ELSET,ELSET=LINER3
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
35,37,2
1001,1041,4
1035,1037,2
*NSET,NSET=DRIV
1,1001
*AMPLITUDE,NAME=SINE,DEF=PERIODIC
1,<omega>,0.0,0.0
0,1,0,0
*STEP
*DYNAMIC,EXPLICIT
,<T>
*IMPEDANCE,PROP=IMPED
LINER3,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=<f_out_hist>
*NODE OUTPUT,NSET=ODB_NSET
POR,
*******************
*END STEP