*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,,0.0,0.0 0,1,0,0 *STEP *DYNAMIC,EXPLICIT , *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= *NODE OUTPUT,NSET=ODB_NSET POR, ******************* *END STEP