This study examines sound transmission of thin membranes with square frame-shaped masses. Numerical results indicate that multiple transmission loss peaks can be generated by adding more frame mass inclusions. The number and the location of the peaks are controlled by the number of frames, the frame distribution, and the frame width. Near the transmission loss peak frequencies, the dynamic effective mass density turns from positive to negative. The validity of the present model has been verified by comparing the analytical results with FE results. Two types of cell arrangements are also considered in this study, namely, cells in series and cells in array. It is seen that either the stacked or array configurations can produce better sound attenuation than single-celled structures. Moreover, the frequency band where sound wave is blocked can be broadened by stacking more layers with different mass magnitudes. Furthermore, additional frequency bands due to the periodicity of the structure are found in the stacked configurations.
1. Introduction
Broadband noise reduction has been investigated over decades. Many researchers have demonstrated that the addition of foams or fiber batting can effectively attenuate sound at certain frequencies [1, 2]. Impedance mismatch of gap layer and the use of microperforated panels/mass inclusions also have been feasible tools for enhancing sound insulation [3–6]. Recently, negative dynamic mass density structures become popular candidates for attenuating elastic waves [7–12] or acoustic waves [13–22] over a certain band of frequencies. For the attenuation of elastic waves, the first proposed promising candidate was a one-dimensional mass-in-mass system usually comprised of a host mass connected to an auxiliary mass by massless and elastic springs. This class of metamaterials is able to create a resonant-type bandgap where waves cannot propagate freely. To acquire multiple bandgaps at desired frequency range, two-dimensional designs of metamaterials, such as metacomposites [11] and multiscale mass-in-mass systems [12], were subsequently introduced. For sound isolation, the basic unit of the pioneering composite is usually comprised of an elastic membrane and a centrally loaded mass. The desired frequency band where sound waves cannot propagate can be achieved by tailoring the mass weight. To eliminate sounds with multiple frequencies, multicelled structures were also proposed [16–18]. Different magnitudes of mass attached to each of the cells in the array can create a multipeak transmission loss profile. In addition to using four-celled array, geometric variation of inclusions is an alternative solution to effectively increase frequency bandwidth and produce multiple transmission loss peaks [19]. Their findings show that the characteristics of transmission loss profile can be manipulated by the number of rings and the distribution of mass. Phononic crystals constituted by periodical repetition of inclusions in a background material also can create frequency gaps which are known as Bragg bandgaps and usually occurred at high frequencies [23, 24]. Because Bragg-type bandgaps are not practice for filtering low-frequency waves, a new design of phononic crystals with locally resonant structures, which can generate flexible resonant-type bandgaps, was developed [25–28]. By altering the properties of the resonant units, the opening and closing of the low-frequency gaps and the gap frequency ranges can be controlled or shifted. Other types of acoustic metamaterials, such as an array of side holes on a tube or an array of thin elastic membranes, have demonstrated their potential capability of filtering low-frequency sound [29, 30]. Those structures exhibit either a negative effective modulus or a negative effective density in a certain frequency range where acoustic waves produced by the excitation are completely blocked.
Sandwich panels containing lightweight honeycomb core are widely used in aircraft and satellite launch vehicles [6, 31]. Although sandwich panels have good benefit of static strength, they are relatively weak in sound insulation. To improve their performance in acoustic absorption, the use of the aforementioned efficient tools becomes necessary. Since the honeycomb core usually has a shape of square form or hexagon form, membrane-frame acoustic metamaterials might be more suitable selections. In this study, the acoustic response of square membranes with multiple frame masses is presented. The governing equation of the proposed vibroacoustic problem is obtained by using the energy method. Parameter studies including frame mass density, frame location, and frame width on transmission loss are emphasized. Two types of structure arrangements, a series-type structure and an array-type structure, are also introduced and compared with single-celled structures. The dynamic effective mass density of the system is also studied.
2. Model
Figure 1 shows a novel design of a basic unit that can be used to construct a membrane-type locally resonant acoustic metamaterial. This structure unit consists of a thin tensioned membrane with one/multiple frame mass (mass m) embedded in the membrane. Consider a plane sound wave of frequency ω and amplitude A normally incident upon the membrane from the region z<0. Here we assume that the bending stiffness of the membrane is overlooked and the frame mass does not hinder the bending of the membrane segment.
A sketch of membrane-frame structure.
The governing equations for the present structure can be derived by using Hamilton’s principle [32] which states(1)∫t1t2δT-U+Vdt=0,where U, T, V are deformation, kinetic, and potential energy, respectively. We have(2)T=12∫0b∫0aρmembrane+∑i=1Qρmassiħi-ħi′∂w∂t2dxdy,U=12∫0b∫0aF∂w∂x2+∂w∂y2dxdy,V=-∫0b∫0aP--P+wdxdy,in which a, b, ρmembrane denote the length, width, and density per unit area of the membrane, ρmassi denotes the density per unit area of the ith frame mass, Q is the number of frame masses, F is the intensity of tension, P- and P+ are acoustic pressure field in the negative and positive z direction, w=w(x,y,t) is the transverse displacement in the z direction of the membrane-frame structure, and ħi-ħi′ represents the area of the ith frame mass as shown as follows:(3)ħi-ħi′=Hx-x0i-Hx-x0i-ci·Hy-y0i-Hy-y0i-di-Hx-x0i′-Hx-x0i′-ci′·Hy-y0i′-Hy-y0i′-di′,where H is Heaviside function. Substituting (2) into (1) gives the equation of motion; namely,(4)ρmembrane∂2w∂t2+∑i=1Qρmassiħi-ħi′∂2w∂t2-F∇2w=P--P+,where P--P+ is the pressure difference between the left-hand side and right-hand side of the membrane at z=0. It can be simply expressed as [33](5)P--P+=2Aejωt-2ρaircair∂w∂t,in which A is plane sound wave pressure amplitude, ω is the wave frequency, ρair is the density of the air, and cair is the speed of sound. To apply Galerkin discretization procedure [34] to (4), the transverse displacement function w is expressed as(6)wx,y,t=∑n=1NWnx,yunt,where Wnx,y is the mode function and un(t) is the harmonic time function excited by the acoustic load; that is,(7)unt=u^nejωt.The mode shape of the membrane simply supported on all four sides can be expressed as [35](8)Wnx,y=sinkπxasinlπybk,l=1,2,….Substituting (6) and (7) into (4), multiplying each term by Wm(x,y), and integrating all terms in the equation over the domain (0≤x≤a,0≤y≤b) give(9)-ω2ρmembrane∫0b∫0aWm∑n=1NWndxdyu^m-ω2∑i=1Qρmassi∑n=1N∫y0iy0i+di∫x0ix0i+ciWmWndxdy-∫y0i′y0i′+di′∫x0i′x0i′+ci′WmWndxdyu^n+jω2ρaircair∫0b∫0aWm∑n=1NWndxdyu^m+-F∫0b∫0aWm∇2∑n=1NWndxdyu^m=2A∫0b∫0aWmdxdy.Let (10)Mm=ρmembrane∫0b∫0aWm∑n=1NWndxdy,Smni=∫y0iy0i+di∫x0ix0i+ciWmWndxdy-∫y0i′y0i′+di′∫x0i′x0i′+ci′WmWndxdy,Cm=2ρaircair∫0b∫0aWm∑n=1NWndxdy,Km=-F∫0b∫0aWm∇2∑n=1NWndxdy,Dm=∫0b∫0aWmdxdy.Equation (9) can be simplified as(11)-ω2Mmu^m+∑i=1Qρmassi∑n=1NSmniu^n+jωCmu^m+Kmu^m=2ADmm=1,2,…,N.Rewriting (11) in the matrix form gives(12)-ω2M+R+jωC+Ku^=2AD,where M, K, and C are N×N diagonal matrices, D, W, and u^ are N×1 matrices, and R is also N×N matrix and has a form as follows:(13)R=∑i=1QρmassiS11iS12i⋯S1NiS12iS22i⋯S2Ni⋮⋮⋱⋮SN1iSN2i⋯SNNi.The solution of displacement amplitude w^(x,y) also can be written in the matrix form as(14)w^x,y=∑n=1NWnx,yu^n=WTu^,where u^ is acquired from (12). Hence we can get the sound power amplitude transmission coefficient τ defined as the ratio of transmitted to incident intensities [33]:(15)τ=ρaircairv^A2,where v^ is the average velocity amplitude of the entire membrane. Then the sound transmission loss TL can be obtained by(16)TL=10log101τ.To treat the whole structure as a homogeneous membrane, the dynamic effective mass density is introduced and acquired by(17)ρeff=p^t-p^ia^,where a^ is the averaged out-of-plane acceleration of the structure, which is normal to the membrane surface at rest; p^t-p^i is the averaged pressure difference between the transmitted wave and the incident wave.
3. Numerical Results
Consider a single-celled structure comprised of a square membrane with an embedded frame mass as shown in Figure 1. The characteristics of sound transmission of the present structure were evaluated using the analytical approach and finite element method. In the finite element analysis, the commercial code COMSOL Multiphysics is used to capture the dynamic behavior of the membrane-frame structure. All the physical properties of the membrane material are displayed in Table 1. Table 2 lists the properties (mass, width, and location) of frame-shaped inclusions embedded in the membrane.
Physical properties of membrane material [16].
Property
Polyetherimide (PEI)
Mass density (kg/m3)
1200
Young’s modulus (Pa)
3.6×109
Poisson ratio
0.36
Thickness (mm)
0.076
Length (mm)
27.4
Tension (Pa)
6.4×106
Properties of frame-shaped mass inclusions.
Configuration
Central mass (g)
Frame I (g)
Frame II (g)
Distance from the center: Frame I (mm)
Distance from the center: Frame II (mm)
Frames I and II width (mm)
Model I
0.32
—
—
—
—
—
Model II
—
0.32
—
2.14
—
0.856
Model III
0.32
0.32
—
2.14
—
0.856
Model IV
—
0.32
0.32
2.14
3.85
0.856
Figure 2 shows the transmission loss of four illustrative examples. Analytical results are obtained by using (16); FE results are calculated from the ratio of the incident power to the transmitted power. It is found that the agreement between analytical and finite element methods is quite good. In Figures 2(a) and 2(b), only one TL peak and two TL dips occur while in Figures 2(c) and 2(d), two TL peaks and three TL dips are found. It is evident that the addition of frame masses can produce more TL peaks. The TL dips correspond to resonance behavior where sound transmission reaches a maximum, while the TL peaks correspond to antiresonance where sound transmission is minimized.
Transmission loss of (a) Model I, (b) Model II, (c) Model III, and (d) Model IV.
Figures 3(a)–3(d) exhibit the dynamic effective mass density of four membrane structures (Model I, Model II, Model III, and Model IV), respectively, as a function of the frequency. Analytical results are acquired by using (17); FE results are obtained by dividing the average pressure on the membrane surface by the average membrane acceleration. It is seen that near the TL peak frequencies, negative effective mass occurs. Observation of the results shown in Figures 3(c) and 3(d) presents that the addition of a frame mass leads to an additional jump where the effective mass turns from positive to negative.
Dynamic effective mass of (a) Model I, (b) Model II, (c) Model III, and (d) Model IV.
Figure 4 illustrates the effect of the material and geometrical properties of the embedded frame on the TL plot of the membrane-frame structure. The results of Figure 4(a) indicate that the first TL peak and valley have a strong dependency on the magnitude of the frame mass. With an increase of the frame mass, the first peak and valley move toward low-frequency region. As the frame width increases, the first TL peak and valley slightly shift to higher frequencies (Figure 4(b)). The first TL valley is from the frame mass and membrane vibrating coherently while the second TL valley is from the membrane vibration with the motionless frame mass. The first TL peak is viewed as the result of superposition of two resonant eigenmodes and has a strong dependency on the frame mass. Hence, as shown in Figure 4(a), the change in frame mass magnitude greatly influences the first TL peak and valley frequencies but has no effect on the second TL valley. On the contrary, with the same mass magnitude and frame location, the change in frame width significantly influences the second TL valley but has a relatively small impact on the first TL peak and valley (as shown in Figure 4(b)). Figure 4(c) exhibits the transmission loss varying with the locations of the frame masses. The frame widths for three cases remain the same. The distances from the center for Frame I, Frame II, and Frame III are 2.14 mm, 3.85 mm, and 5.56 mm, respectively. As the frame moves away from the center, a notable increase in the first TL peak and valley frequencies is observed. The reason is possibly that different frame mass locations result in different eigenmodes of the membrane-frame structure. The characteristic frequencies of the first resonance and TL peak varying with the location of the frame mass are presented in Figure 5. It seems that the first TL peak and resonance frequencies are almost proportional to the mass location which is measured from the membrane center. Also, it is seen that the former rising rate is higher than the latter.
Transmission loss of membrane-frame structure with varying (a) frame mass, (b) frame width, and (c) frame location.
The dependence of the first TL peak and TL valley upon the mass location.
The schematics of the stacked membrane-central-mass structure and stacked membrane-frame structure are displayed in Figures 6(a) and 6(b), respectively. Four different mass magnitudes are considered: 0.08 g, 0.16 g, 0.24 g, and 0.32 g. Figure 7 shows the comparison of sound transmission of the stacked structure, obtained by FEA analysis, with that of four single-celled structures, obtained by analytical models. It is seen that the TL peaks and valley frequencies of the stacked structure coincide with the ones of single-celled structures. In other words, the characteristic frequencies of the series-type structure can be precisely predicted by single-layer structures. It is also seen that the width of frequency band in which sound wave cannot propagate is effectively enlarged. Compared with the results shown in Figures 7(a) and 7(b), it seems that the acoustic performances of two stacked structures present a similar behavior. Moreover, more frequency bands due to the periodicity of the structure are created in both cases. Figure 8 exhibits the effects of air spacing on sound transmission of two structures. It is evident that air spacing does affect the frequency band caused by the periodic assembly but has no effect on the TL peak and valley frequencies, caused by resonance and antiresonance of the local resonator. With larger air spacing, the band width becomes wider and the beginning frequency of the band moves toward lower frequency regime. Figure 9 depicts two sketches of the array-type structures. Two unit cells are considered and shown in Figures 9(a) and 9(b), respectively, a membrane with a central mass and the one with a frame mass. The results show that a multipeak TL profile can be created by using different masses at adjacent cells. Sound transmission of the array-type structure and four single-celled structures are compared and displayed in Figure 10. From the results of Figures 10(a) and 10(b), it is seen that the TL valleys for four-celled array correspondingly match the ones for single-celled structures. In other words, the resonances of the local resonators do not depend on the cell arrangement whereas the antiresonance frequencies are affected. As expected, series-type structures can produce stronger attenuation in sound wave than array-type structures and single-celled structures.
Schematics of the structures stacked in series with (a) membrane-central-mass units and (b) membrane-frame units.
(a) Comparison of membrane-central-mass structures stacked in series with four single-celled structures and (b) comparison of membrane-frame structures stacked in series with four single-celled structures.
Effect of air spacing on sound transmission for (a) stacked membrane-central-mass structure and (b) stacked membrane-frame structure.
Schematics of the structures arranged in array with (a) membrane-central-mass units and (b) membrane-frame units.
(a) Comparison of membrane-central-mass structures arranged in array with four single-celled structures and (b) comparison of membrane-frame structures arranged in array with four single-celled structures.
4. Conclusions
This study addresses the acoustic behavior of thin membranes with square frame-shaped masses. Analytical results are consistent with FE results in transmission loss and effective mass. The results indicate that a multipeak TL profile can be created by adding more frame masses in the membrane. Frequencies where negative effective mass occurs coincide with the TL peak frequencies. The TL peak and valley frequencies can be tailored by altering the geometrical and material properties of the frame. Two types of structure arrangements, a series-type structure and an array-type structure, are also evaluated and compared with single-celled structures. It is found that structures stacked in series can produce a frequency band with wide bandwidth and enhance attenuation in sound wave.
Competing Interests
The authors declare that they have no competing interests.
Acknowledgments
The authors would like to thank Dr. Y. T. Chiu and Dr. T. C. Lu from Industrial Technology Research Institute for their constructive comments. This work was supported by the Ministry of Science and Technology of Taiwan under Grant MOST 103-2221-E-006-032 and in part by the Headquarters of University Advancement at the National Cheng Kung University, which is sponsored by the Ministry of Education, Taiwan.
BoltonJ. S.ShiauN.-M.KangY. J.Sound transmission through multi-panel structures lined with elastic porous materials1996191331734710.1006/jsvi.1996.01252-s2.0-0030568501LaneS. A.GriffinS.RichardR. E.Fairing noise mitigation using passive vibroacoustic attenuation devices2006431314410.2514/1.139262-s2.0-33644760427MaaD. Y.Theory and design of microperforated panel sound-absorbing constructions1975185571KidnerM. R. F.FullerC. R.GardnerB.Increase in transmission loss of single panels by addition of mass inclusions to a poro-elastic layer: experimental investigation2006294346647210.1016/j.jsv.2005.11.0222-s2.0-33646518989IdrisiK.JohnsonM. E.TosoA.CarnealJ. P.Increase in transmission loss of a double panel system by addition of mass inclusions to a poro-elastic layer: a comparison between theory and experiment20093231-2516610.1016/j.jsv.2008.12.0262-s2.0-64449087855NaifyC. J.HuangC.SneddonM.NuttS.Transmission loss of honeycomb sandwich structures with attached gas layers2011722-3717710.1016/j.apacoust.2010.09.0052-s2.0-78649463986YaoS.ZhouX.HuG.Experimental study on negative effective mass in a 1D mass-spring system20081004302010.1088/1367-2630/10/4/0430202-s2.0-46349109609HuangH. H.SunC. T.Wave attenuation mechanism in an acoustic metamaterial with negative effective mass density20091101300310.1088/1367-2630/11/1/0130032-s2.0-59849106849HuangH. H.SunC. T.HuangG. L.On the negative effective mass density in acoustic metamaterials200947461061710.1016/j.ijengsci.2008.12.0072-s2.0-60449109154YaoS.ZhouX.HuG.Investigation of the negative-mass behaviors occurring below a cut-off frequency20101210302510.1088/1367-2630/12/10/1030252-s2.0-77958580696HuangH. H.SunC. T.Continuum modeling of a composite material with internal resonators20124611010.1016/j.mechmat.2011.11.0062-s2.0-84055191147AhmedR. U.BanerjeeS.Wave propagation in metamaterial using multiscale resonators by creating local anisotropy20131325159YangZ.MeiJ.YangM.ChanN. H.ShengP.Membrane-type acoustic metamaterial with negative dynamic mass20081012020430110.1103/PhysRevLett.101.2043012-s2.0-56849128459NaifyC. J.ChangC.-M.McKnightG.NuttS.Transmission loss and dynamic response of membrane-type locally resonant acoustic metamaterials20101081111490510.1063/1.35140822-s2.0-78751508080YangZ.DaiH. M.ChanN. H.MaG. C.ShengP.Acoustic metamaterial panels for sound attenuation in the 50–1000 Hz regime201096404190610.1063/1.32990072-s2.0-75749151973NaifyC. J.ChangC.-M.McKnightG.ScheulenF.NuttS.Membrane-type metamaterials: transmission loss of multi-celled arrays20111091010490210.1063/1.35836562-s2.0-79958844508ZhangY.WenJ.ZhaoH.YuD.CaiL.WenX.Sound insulation property of membrane-type acoustic metamaterials carrying different masses at adjacent cells2013114606351510.1063/1.48184352-s2.0-84883430150NaifyC. J.ChangC.-M.McKnightG.NuttS. R.Scaling of membrane-type locally resonant acoustic metamaterial arrays201213242784279210.1121/1.47449412-s2.0-84867391742NaifyC. J.ChangC.-M.McKnightG.NuttS.Transmission loss of membrane-type acoustic metamaterials with coaxial ring masses20111101212490310.1063/1.36652132-s2.0-84855301489ZhangY.WenJ.XiaoY.WenX.WangJ.Theoretical investigation of the sound attenuation of membrane-type acoustic metamaterials2012376171489149410.1016/j.physleta.2012.03.0102-s2.0-84862784968VaranasiS.BoltonJ. S.SiegmundT. H.CipraR. J.The low frequency performance of metamaterial barriers based on cellular structures201374448549510.1016/j.apacoust.2012.09.0082-s2.0-84878063774ChenY. Y.HuangG. L.ZhouX. M.HuG. K.SunC.-T.Analytical coupled vibroacoustic modeling of membrane-type acoustic metamaterials: membrane model2014136396997910.1121/1.48928702-s2.0-84911424374SigalasM. M.EconomouE. N.Elastic and acoustic wave band structure1992158237738210.1016/0022-460X(92)90059-72-s2.0-0026941042KushwahaM. S.HaleviP.DobrzynskiL.Djafari-RouhaniB.Acoustic band structure of periodic elastic composites199371132022202510.1103/PhysRevLett.71.20222-s2.0-3242887087LiuZ.ChanC. T.ShengP.Analytic model of phononic crystals with local resonances200571101410310.1103/physrevb.71.0141032-s2.0-16844373388OudichM.LiY.AssouarB. M.HouZ.A sonic band gap based on the locally resonant phononic plates with stubs20101208304910.1088/1367-2630/12/8/0830492-s2.0-77956588005HsuJ.-C.Local resonances-induced low-frequency band gaps in two-dimensional phononic crystal slabs with periodic stepped resonators201144505540110.1088/0022-3727/44/5/0554012-s2.0-79251573858BanerjeeS.AhmedR. U.Phonon confinement using spirally designed elastic resonators in discrete continuum20143161310.11648/j.ijmsa.20140301.12LeeS. H.ParkC. M.SeoY. M.WangZ. G.KimC. K.Acoustic metamaterial with negative modulus2009211717570410.1088/0953-8984/21/17/1757042-s2.0-65149090707LeeS. H.ParkC. M.SeoY. M.WangZ. G.KimC. K.Acoustic metamaterial with negative density2009373484464446910.1016/j.physleta.2009.10.0132-s2.0-74249112401QianZ.FuY.YuP.DengZ.AnZ.ChenJ.ZhouZ.ChenW.A new method for calculating noise reduction of sandwich structures based on modal testingProceedings of the 21st International Congress on Sound and VibrationJuly 2014Beijing, China18WellsteadP. E.1979London, UKAcademic PressMR570825FahyF.GardonioP.1985Orlando, Fla, USAAcademic PressKopmazO.TelliS.Free vibrations of a rectangular plate carrying a distributed mass20022511395710.1006/jsvi.2001.39772-s2.0-0037149204MukhopadhyayM.2000Brookfield, Vt, USAA. A. Balkema Publishers