Research Article  Open Access
LowFrequency Broadband Sound Transmission Loss of Infinite Orthogonally RibStiffened Sandwich Structure with Periodic Subwavelength Arrays of Shunted Piezoelectric Patches
Abstract
This paper studies lowfrequency sound transmission loss (STL) of an infinite orthogonally ribstiffened sandwich structure flexibly connected with periodic subwavelength arrays of finite shunted piezoelectric patches. A complete theoretical model is proposed by three steps. First, the panels and piezoelectric patches on both sides are equivalent to two homogeneous facesheets by effective medium method. Second, we take into account all inertia terms of the ribstiffeners to establish the governing equations by space harmonic method, separating the amplitude coefficients of the equivalent facesheets through virtual work principle. Third, the expression of STL is reduced. Based on the two prerequisites of subwavelength assumption and convergence criterion, the accuracy and validity of the model are verified by finite element simulations, cited experiments, and theoretical values. In the end, parameters affecting the STL performance of the structure are studied. All of these results show that the sandwich structure can improve the lowfrequency STL effectively and broaden the sound insulation bandwidth.
1. Introduction
Classical sandwich structure consists of two parallel panels (as the facesheets) reinforced by sets of periodic ribstiffeners (as the sound bridges) and filled with air, which can be found in a wide range of engineering applications, such as lightweight aircraft sidewall panels, submarine/ship constructions, and express trains [1–11]. In the lowfrequency range, the classical sandwich structure can be modeled as an orthotropic plate when the ribstiffener spacing is much smaller than a wavelength of the panel flexural wave. However, in the relatively high frequency range, the exact model of the ribstiffeners should be taken into account in the case where the wavelength is comparable with ribstiffener spacing [2, 6, 10, 11]. In view of the vibroacoustic characteristics of sandwich structure, a large number of valuable literatures cover the corresponding theoretical models, simulations, or experiments.
For the study of sandwich structure, more than 20 theoretical models have been presented. More intuitively, only by the detailed comparison of Hongisto [12] can it be found that the prevalent models for the prediction of STL are as high as 16. According to the classical sandwich structure, the sound bridges will significantly reduce the sound insulation in the middle and high frequency bands, and three methods are commonly used in predicting the vibroacoustic response. The first method is the decoupled approach, which divides the sound transmission paths into a fluidborne through the air cavity and a structureborne through the bridges. The sound transmission coefficient of the fluidborne path can usually be analyzed by the transfer matrix method (TMM) [8, 13]; the pressures and velocities at the top and bottom facesheets of the sandwich structure form a fourterminal network. To solve the sound transmission coefficient of the structureborne path, Sharp [14] assumed that the velocity of the first facesheet was unaffected by the connections, and the bridges were rigid and massless; then a theoretical model was proposed. But the shortcoming is that the STL curve was required to add a positive empirical correction of 5 dB. Gu and Wang [15] extended Sharp’s theory by modeling the bridge connections as equivalent translational springs, and then the sandwich structure was simplified as a classical massspringmass system. The model developed by Fahy and Gardonio [16] has similarities and differences with Sharp’s model, it is based on the assumption that the dynamic properties of each stud are independent of each other, and their motions are restricted to translation, caused by incident bending waves normal to the line connections between the facesheets and studs. Davy [17] relied on Fahy’s theory; he supposed that the bridges have flexible properties and facesheets have different line impedances. Legault and Atalla [8] summed up the previous four theoretical models and obtained the conclusion that Davy’s model is good over the whole frequency range 100 Hz–10 kHz by comparing with the experimental values.
The second method is the Fourier transform technique, which has a wide range of applications. For example, Rumerman [18] proposed a general solution for the forced vibration problem of an infinite thin plate with periodically uniform ribstiffeners, but this method only considers the tension forces and bending moments between the facesheet and ribstiffeners, and the impact of the torsional moments is ignored. Lin and Garrelick [19] calculated the sound transmission of a sandwich structure with rigid bridges. Takahashi [20] discussed the sound radiation of the stiffened sandwich structure with multipoint connections excited by mechanical point forces, but the governing equations did not consider the fluid loads. Recently, Maxit [21] studied the vibration and sound radiation of a periodically ribbed plate to a mechanical point force drive by an efficient method based on the Fourier transform technique. Xin and Lu [11] proposed a theoretical model to calculate the sound radiation from an infinite sandwich structure reinforced by orthogonal ribstiffeners; it is worth mentioning that the motions of equally spaced ribstiffeners are determined by considering their tensional forces, bending moments, and torsional moments.
The third method is a hybrid of space harmonic method and virtual work principle. In advance, Mead [22, 23] analyzed the responses of the structure by the space harmonic method and then studied the dynamics of a periodic unit of the system based on the virtual work principle. Further, Lee and Kim [6] proposed a theoretical analysis model based on the hybrid method, to calculate STL performances of a thin plate stiffened by equally spaced line stiffeners; each stiffener was simplified as a combination of a lumped mass and rotational and translational springs. Extending this model to sandwich structure stiffened with parallelly placed studs, Wang et al. [7] developed a periodic model and took account of the discrete placing of the studs. However, this model only included the effect of tensional forces and bending moments between the studs and facesheets and ignored the torsional moments. On the basis of the above two models, Xin and Lu [10] perfected a deterministic theoretical model by coupling the acoustic and structure vibration and studied the STL of an infinite sandwich structure reinforced by orthogonal ribstiffeners. All possible motions of the ribstiffeners, such as tensional forces, bending moments, and torsional moments as well as the corresponding inertial terms, were taken into account in the governing equations of the two facesheets.
When the experiment of STL through sandwich structure was mentioned, the experimental results measured by the National Research Council of Canada (NRCC) were more comprehensive and authoritative, its internal report IRCIR761 [24] contained 350 sets of measurements, and the experimental specimens were mostly single and doubleply gypsum board attached parallel to steel/wood/nonloadbearing steel studs. In the existing literatures on the STL of stiffened plate structure, the experimental results in some literatures [25, 26] were directly referenced to the measurement data in the internal report IRCIR761 of NRCC [24] and in some documents [8] are measured by the relevant research institutions. In addition, there are parts of the literatures based on experimental study; for example, a series of articles by Hongisto [27–30] measured the structural transmission through the door leaf, the leak transmission through the slits, and flanking transmission via double structures by sound intensity method or sound pressure method. Among them, the test laboratory, measurement equipment, test specimens, and test methods were described in detail.
Although there are many studies on the STL of sandwich structure, only the middle and high frequency bands have a good sound insulation effect, while the STL performance in lowfrequency range is still not improved. Therefore, this work first introduces periodic subwavelength arrays of finite shunted piezoelectric patches into an infinite orthogonally ribstiffened sandwich structure and aims to enhance the STL in the lowfrequency range. However, the existing research on periodic subwavelength arrays of shunted piezoelectric patches is basically focused on the vibration control, as well as noise radiation from metamaterialbased thin plates [31–35]; study into STL remains to be accomplished, not to mention the lowfrequency broadband STL. For example, Xiao et al. [36] were concerned with STL of metamaterialbased thin plates consisting of multiple subwavelength arrays of springmass resonators and completed the prediction by two analytical wave approaches: the plane wave expansion (PWE) and effective medium (EM) method. Zhang et al. [37] introduced the piezoelectric shunt damping into the periodic subwavelength structure and developed EM method for calculation of STL, but it was insufficient that the connections between the piezoelectric patches and substrate are rigid.
In view of the infinite orthogonally ribstiffened sandwich structure flexibly connected with periodic subwavelength arrays of finite shunted piezoelectric patches (hereafter termed research structure), the studies on STL prediction of the research structure are blank. Therefore, this work focuses on the lowfrequency broadband STL prediction of the research structure and carries out studies from the following four aspects. First, a complete theoretical model of STL for solving the studied sandwich structure can be built based on the EM method, space harmonic method, and virtual work principle. Second, two prerequisites for the application of the theoretical model are given. Next, the accuracy of the theoretical model is verified by finite element simulation and cited experimental data, and the validity is verified by the comparison of the STLs before and after adding the shunt piezoelectric patches. Afterwards, the parameters affecting the STL performance of the research structure are studied.
2. Theoretical Modeling of Vibroacoustic Response
2.1. Preliminary Considerations
Considering the research model is an infinitely lager twodimensional (2D) sandwich structure as shown in Figure 1. It consists of two infinite thin plates connected with periodic subwavelength arrays of shunted piezoelectric patches on the inner side and the two parallel facesheets reinforced by an orthogonally ribstiffened core which has two periodic uniform spacing of and in the  and directions, respectively. Its geometrical parameters are thickness of top and bottom plates and , depth of orthogonal ribstiffeners , and thickness of  and wise stiffeners and , respectively. A righthanded Cartesian coordinate system is established on the upper surface of top plate, the positive direction of axis is vertical downward, and the positive directions of the  and axis, are, respectively, along the two orientations of the orthogonal ribstiffeners. The research structure divides the spatial sound field into three parts: the top sound field (semiinfinite sound field) in the area of < 0, the middle sound field (divided into periodic cavities by the orthogonal ribstiffeners) in between the area of , and the bottom sound field (semiinfinite sound field) in the area of .
An acoustic plane wave varying harmonically in time impinging on the upper surface of the top plate is included in the sandwich structure. The incident plane wave induces motion of the top plate, and the motion creates a pressure inside the middle sound field and drives the deformations of the orthogonal ribstiffeners. Therefore, the acoustic energy can travel from the top plate to the bottom plate via two separate paths: the fluidborne path through the middle sound field and the structureborne through the orthogonal ribstiffeners. The incident wave makes an elevation angle with the axis and its projection in the plane makes an azimuth angle with the axis. Its amplitude is and the sound wavenumber vector can be decomposed in the , , and directions [8]:where in (2) is the sound speed in air.
2.2. Equivalent Characteristics of Periodic Facesheets
In view of the spatial periodicity of the research structure as shown in Figure 1, it is necessary to consider only the vibroacoustic response contribution from one period of element which is shown in the red cube box in Figure 2(a). Among them, assuming the ribstiffener spacing , that is, the top and bottom facesheets of one period are square structure, their inner surfaces each contain pieces of piezoelectric patches, with the formation of a periodic subwavelength structure of rows and columns (see Figure 2(b)). The yellow box is the basic unit of periodic subwavelength structure in Figure 2(a); its top view is shown in Figure 3(a), which is divided into a multilayer region (m) and a single layer region (s); the lattice constant and size of piezoelectric patch are, respectively, and .
(a)
(b)
(a)
(b)
The flexible connections between the piezoelectric patches and facesheets can be initially equivalent to a spring and a damper in a unit cell. The equivalent spring, damper, and piezoelectric patch form a resonator. Among them, the stiffness coefficient is (), where subscript represents the th resonator, and the damping coefficient can be introduced by means of a complex stiffness coefficient constant using the symbol to represent it, and it satisfies the relationship: ), where is the loss factor. Therefore, the flexible connection unit can be finally equivalent to the structure shown in Figure 3(b). In addition, an external shunt circuit is connected to the piezoelectric patch of a unit cell as shown in Figure 3(b); its impedance is .
To achieve the equivalent surface mass density and equivalent dynamic bending stiffness of the top and bottom facesheets of one period, it is based on EM method and can be completed by three steps. The first step is to determine the effective material parameters of a single piezoelectric patch with external shunt circuit. The second step is to derive the effective material parameters of regions m and s of the unit cell. And the last step is to determine the effective material parameters of the top and bottom periodic facesheet.
Assuming that the displacement of each section of the spring has a linear distribution, the equivalent mass of the individual resonator can be obtained as follows:where , , and , respectively, denote the volume mass density, thickness, and area of the shunted piezoelectric patch of the th resonator, is the angular frequency, is the mass of the th equivalent spring, and represents the angular frequency of th resonator corresponding to the resonance frequencyNote that ; it can be regarded as a massless spring.
In addition, the piezoelectric patch of th resonator can be equivalent to an isotropic plate, whose equivalent Young’s modulus and Poisson ratio are, respectively, as follows:where is the Laplace operator () and and are, respectively, the th piezoelectric patch compliance coefficient and piezoelectric constant. Subscripts 1, 2, and 3 are the three components of the corresponding parameter; its vector points along the direction of the three coordinate axes. is the inherent capacitance of th piezoelectric patch at constant stress, which can be expressed as where is the dielectric constant of the th piezoelectric patch at constant stress.
Therefore, the equivalent surface density and equivalent dynamic bending stiffness of regions and of the periodic facesheet can be obtained aswhere subscripts and 2 represent the top periodic facesheet and bottom periodic facesheet. , , , and are, respectively, the volume mass density, thickness, Young’s modulus, and Poisson ratio of the basic material of facesheet.
Furthermore, the expressions of the equivalent surface mass density and equivalent dynamic bending stiffness of the periodic facesheets are obtained by EM method as follows [37]:where denotes the ratio between the total area of region and the total area of region in one periodic facesheet, given by
2.3. Dynamic Motions of Periodic Facesheets and RibStiffeners
To a twodimensional sandwich structure, reinforced by an orthogonally ribstiffened core repeats in the  and direction and excited by a harmonic plane sound wave, the top and bottom facesheet responses can be expressed using space harmonic expansion as follows [10, 23, 38]:where subscripts = 1 and 2 also represent the top periodic facesheet and bottom periodic facesheet; the th harmonic wave components in the top and bottom facesheets have the same wavenumbers (, ) but different amplitudes: that is,Note that the wavenumbers , denote positivegoing harmonic waves in the  and direction and vice versa.
The orthogonal ribstiffeners impose a strong constraint on the motions of the facesheets and then affect the sound transmission and wave propagation, where the influences of ribstiffeners exist in the form of tensional forces, bending moments, and torsional moments, which satisfy the following expressions: tensional forces = general force + inertial force, bending moments = general bending moment + inertial bending moment, and torsional moments = general torsional moment + inertial torsional moment. The impacts at the interfaces between the top facesheet and /wise ribstiffeners are marked by (, , ) and (, , ) as shown in Figures 4(a) and 4(b), respectively. Similarly, the impacts at the interfaces between the bottom facesheet and x/ywise ribstiffeners are, respectively, marked by (, , ) and (, , ) as shown in Figure 5.
(a)
(b)
(a)
(b)
Taking all the impacts of the ribstiffeners into consideration, and applying Hooke’s law and the Newton’s second law, the expressions of the tensional forces, bending moments, and torsional moments can be obtained as follows [10]:(1)Tensional forces:(2)Bending moments: (3)Torsional moments:whereIn (12) to (14), , , , and are tensional force coefficients, , , , and are bending moment coefficients, and , , , and are torsional moment coefficients. Their calculation expressions are as follows:(1)Tensional force coefficients:(2)Bending moment coefficients:(3)Torsional moment coefficients:where (, ), (, ), and (, ) are, respectively, the tensional, bending, and torsional stiffness of half ribstiffeners per unit length. and are the line and volume mass density for the ribstiffeners, respectively. and are polar and torsional moment of inertia of the ribstiffeners, where subscripts and denote the corresponding orientations of the orthogonal ribstiffeners. The geometrical properties of orthogonal ribstiffeners are given by [10, 11]
2.4. Continuity Conditions
The sound pressure in the top, middle, and bottom sound field is marked by , , and , which can be expressed by space harmonic series as [7, 10, 23] where is the th space harmonic wavenumber in the direction, which can by calculated from [6, 39]when , the corresponding pressure waves become evanescent, and then the expression is replaced by
Moreover, the sound pressures in the three sound fields meet the wave equation [40]:
The continuity condition of fluidstructure coupling is to ensure that the facesheet velocity is equal to the fluid velocity in the facesheet surface, which is based on the momentum equation and satisfies the relationships as follows [19, 41]:where is the density of fluid field.
By substituting (10) and (20) into (24), the relevant coefficients can be obtained as
2.5. Separation of Facesheet Response Amplitudes
As can be seen from (25) to (28), the facesheet response amplitudes and determine the motions of the facesheets and ribstiffeners. Based on one periodic element, the amplitudes and can be derived by applying the virtual work principle [6, 7, 10, 23]. Therefore, the virtual placements of the top and bottom facesheets can be obtained aswhere is the Dirac delta function.
The equations governing the vibration responses of the periodic facesheets are given by
Moreover, the virtual works of periodic facesheets and wise ribstiffeners are shown as follows:(1)Periodic facesheets:(2)xwise ribstiffeners:(3)ywise ribstiffeners:where and are the modal amplitudes and
Finally, the virtual work principle satisfies the fact that
Substituting (10) and (20) into (31)–(33), then putting the simplified results into (35), it can be derived that where (36) describe an infinite set of coupled equations. The solution becomes convergent when the single mode number is greater than a constant, that is, the sumindex in the finite range of and . Therefore, the simultaneous equations can be written as follows in matrix form of finite order (i.e., 2KL, where , )
The detailed derivation of (37) is similar to that in the literature [6, 7, 10, 11].
2.6. Sound Transmission Loss
From the above it can be known that the coefficients ( = 1, 2) have been solved; then the remaining unknown coefficients , , , and can be found by using (25)~(28), respectively. The sound transmission coefficient can then be calculated aswhich is dependent on the elevation angle and azimuth angle . The diffusefield transmission coefficient is finally calculated by integrating over and : where = 90° for random incidence and ≈ 78° for field incidence transmission.
The sound transmission loss (STL) may be defined as [6–8, 10, 11, 39]where is calculated from (38) or (39).
3. Fundamental Conditions of Research
3.1. Basic Parameters
To the STL problem of the periodic sandwich structure as stated above, the geometrical and characteristic parameters have a direct influence on the sound insulation performance, and the STL performance of periodic structure can be characterized by one periodic element.
According to the periodic element as shown in Figure 2(a), the piezoelectric patches of each facesheet are selected as 25 pieces (i.e., = 5). The corresponding external shunt circuits are the equivalent negativecapacitance (NC) circuit as shown in Figure 6, which consists of a fixed resistor (), an adjustable resistor (), a capacitor (C), and an operational amplifier (OA). Among them, the fixed resistor is composed of the internal resistor of piezoelectric patch and the external shunt resistor. The equivalent negative capacitance is given by
Therefore, the impedance of each external shunted circuit can be obtained as
In addition, the properties of the facesheet material (aluminum) and piezoelectric patch (PZT_5H) are represented in Table 1; the meaning of each symbol is as stated above.

3.2. Verification of Subwavelength Hypothesis
Since the EM theory is valid only in the subwavelength regime, it is necessary to verify whether the unit cell in Figure 3(a) satisfies the subwavelength hypothesis, which requires that the lattice constant () is smaller than the flexural wavelength () of the substrate, and is defined by
The ratio () of the lattice constant to the flexural wavelength is plotted in Figure 7 over the frequency range of 20–2000 Hz. Figure 7 describes that the ratio increases with frequency; that is, the calculation accuracy of the EM model decreases with the increasing frequency. It is can be seen that the ratio is smaller than 0.5 and confirms the validity of the hypothesis.
3.3. Criteria for Convergence
For the research sandwich structure, there must be enough truncations about infinite and to ensure the convergence and accuracy of the solution, but when the number of terms is too large, it needs high hardware configuration and low computational efficiency as the cost. The description of a feasible criterion [6] is when the solution converges at a given frequency; it is convergent in the region below that frequency. Therefore, it only needs to determine the necessary number of terms corresponding to the highest frequency (i.e., 2 kHz). Moreover, the convergence check of the sandwich structure is determined by the following three steps. The first step is to calculate the STL performance at the highest frequency. The second step is to add one term to further calculate the corresponding STL. In the last step, in order to find whether the difference between the two results is within the preset error band (0.008 dB in this work), the solution can be considered convergent and the corresponding number of terms is suitable for calculating the STL performances at all other frequencies below 2 kHz.
Since the infinite orthogonally ribstiffened sandwich structure is symmetrical about  and axis, all the equations have been evolved into a finite set of equations with and , and an assumption have been made where . Figure 8 illustrates that the STL decreases gradually with the single mode number ( or ) increasing at the highest frequency of 2 kHz, and finally it tends to be stable. In addition, the results of Figure 8 show that the solution achieved convergence conditions when single mode number (or ) ≥ 11; that is, it requires at least 529 terms to ensure solution convergence.
4. Calculation of Sound Transmission Loss
4.1. Accuracy of the Theoretical Model
To represent the accuracy of the theoretical model, we compare the theoretical values with the finite element simulation results based on VA One software and the reference experimental data is this section, respectively.
4.1.1. Simulation Comparison
For a theoretical model, its feasibility and correctness are the key to development. In this section, the finite element analysis of the research structure is carried out by using VA One, vibroacoustics simulation software developed by the ESI Group. The finite element results are then compared with the theoretical values to verify the accuracy of the theoretical model.
Since VA One can perform finite element simulation of periodic subsystems, it is only necessary to analyze the periodic element in Figure 2(a). For the orthogonally ribstiffened sandwich structure with no shunt piezoelectric patches, the corresponding finite element model is shown in Figure 9(a), which includes a periodic element, two semiinfinite fluids (SIFs), and a diffuse acoustic field (DAF). Consider the finite periodic element with fluid on both sides, the acoustic transmission coefficient through the structure can be derived by the diffusefield reciprocity result [42]. For an eigensolution associated with a pair of phase constants (), the radiation stiffness in the SIF on side is written as follows:where is the spatial Fourier transform of the freefield Green function [43] and is the spatial Fourier transform of displacement field over the entire periodic structure (obtained from the spatial Fourier transform of the eigenvector over a single research element). When denotes the acoustic wavenumber in the source cavity and A is the surface area per periodic element coupled to the fluid on the source side, the transmission loss is obtained from the diffusefield reciprocity result given by the following equation:where the sum runs over the eigensolutions associated with the phase constants pair (, ) and denotes the imaginary part. For a given frequency band , the contribution to the transmission from the modes with natural frequency within the band (i.e., ) defines the resonant transmission and is analytically averaged over the band. Alternatively, the nonresonant transmission is obtained from the contribution of the modes with natural frequency below and above the frequency band.
(a)
(b)
For the finite element model, it requires that the minimum wavelength corresponding to the highest frequency contains at least four grids. The geometrical and material parameters of the periodic element and the fluid properties are shown in Table 1. To improve the solving efficiency, the finite element model is used, 1/6 octave in the frequency range of 10 Hz–2000 Hz. Furthermore, the simulation results obtained by the finite element method (FEM) are compared with the theoretical values obtained by integration when is limited to 78° as shown in Figure 9(b).
As can be seen from Figure 9(b), the STL curve simulated by the FEM has a similar trend with the theoretical calculation curve. Although the simulation curve is not as smooth as the theoretical curve, the corresponding valley frequencies (such as around 425 Hz, 865 Hz, and 1275 Hz) and peak frequencies (such as around 555 Hz, 1115 Hz, and 1725 Hz) can effectively match. In addition, the difference between the two STL curves is within acceptable limits. Therefore, the theoretical model established for the orthogonally ribstiffened sandwich structure can predict the sound insulation.
4.1.2. Cited Experimental Comparison
The feasibility and accuracy of the theoretical model can also be verified by a comparison, between the theoretical prediction values and the experimental results in the internal report IRCIR761 [24], completed by National Research Council of Canada (NRCC) in March 1998. The experimental specimen of number G16_SS90(610)_G16 is made up from two parallel gypsum boards, which are separated by an air cavity and linked with periodically distributed “Csection” steel studs. For this reason, the research sandwich structure is simplified as a typically infinite sandwich structure with unidirectionally parallel ribstiffeners; that is, the theoretical model is not restricted by the subwavelength assumption. In addition, the material and geometrical properties of the structure adopted by the experiment G16_SS90(610)_G16 are completely followed in the comparison.
Figure 10 shows a comparison of the predictions from theoretical model and test data from internal report IRCIR761 [24]. Among them, the STL measurements are taken in 1/3 octave bands from 50 Hz to 6300 Hz, and the STL is calculated in 1/27 octave bands with a diffusefield integration limit of 78° in the same frequency range.
It can be seen from Figure 10 that overall the comparison of theoretical prediction values and experimental results is satisfactory, except for a certain discrepancy in narrowband 200 Hz–300 Hz and in the vicinity of 1000 Hz. Since the influence factors considered in the theoretical model are not comprehensive in the experimental measurement, the theory and experiment must be more or less different. As Hongisto et al. [28] did the STL experiments on nine steel passage doors and nine timber passage doors using the sound intensity method, the average difference between the results and the corresponding theory is in the range of −3 to 12 dB. Therefore, the comparison shows that the theoretical model can be used to analyze the structure characteristics and acoustical properties of the sandwich structure with ribstiffeners as a core accurately.
4.2. Validity of Theoretical Model
To explore the validity of the research sandwich structure composed of 25 piezoelectric patches in each unit cell, the orthogonally ribstiffened sandwich structure with no shunt piezoelectric patches is introduced as a contrast model. The material parameters and geometrical properties of the two calculation models are equal as shown in Table 1. Since the equivalent facesheet parameters of each model are identical, the natural frequencies of sandwich structure with ribstiffened core can be calculated by the following formula:where and are, respectively, the equivalent dynamic bending stiffness and Poisson’s ratio of the facesheet. is the total equivalent surface density of the studied structure. is shear stiffness of the ribstiffener.
The STLs of the two calculation models are compared in the frequency range 10 Hz–2000 Hz as shown in Figure 11. It illustrates that the STL of the research structure is better than that of the contrast structure in the whole frequency band, except for the vicinity of the second valley frequency 1146.9 Hz, which is generated by the (3, 0)order modal of the research structure. More concretely, as the number of piezoelectric patches on the top and bottom facesheet in a typically periodic element increases from 0 to 25, the increment of STL exceeds 20 dB over the frequency range 10 Hz–120 Hz. Moreover, the first valley frequency is moved from 437.3 Hz to 133.1 Hz; it corresponds to the (1, 5)order modal of the contrast structure and the (1, 1)order modal of the research structure, respectively. In addition, it can be concluded from Figure 11 that the multimode resonant frequencies are generated by the resonators, which are formed by the flexible connections and the corresponding piezoelectric patches; they can improve the STL effectively and broaden the sound insulation bandwidth of the research structure.
4.3. Parameter Studies
4.3.1. Effect of Elevation Angle
According to the research structure, the azimuth angle has no influence on the STL in the whole frequency band when the elevation angle is certain. Therefore, only the effect of the elevation angle on the STL is studied under constant azimuth angle in the frequency range of 10 Hz–2000 Hz.
Under the condition of the azimuth angle being identical to 30° (i.e., = 30°), the STLs calculated for three different elevation angles (i.e., = 0°, 30°, and 60°) are plotted in Figure 12, which illustrates that the STL decreases with increasing elevation angle due to the waveform conversion. There is no waveform conversion when = 0° (i.e., normal incident), but with the increases of elevation angle (i.e., oblique incident) the waveform conversion occurs, resulting in the fact that the space harmonic wavenumber in direction changes, and finally the sound transmission coefficient increases as the elevation angle increases. Figure 12 also indicates that the first valley frequency is fixed at 133.1 Hz; this phenomenon is independent of the incident angle and is determined by the inherent characteristics of the research structure.
4.3.2. Effect of RibStiffener Spacing
Ribstiffener spacing and spacing have a great impact on the properties of structureborne and fluidborne transmission. Therefore, on the basis that the ratio between the total area of region m and the total area of region s in one periodic facesheet is fixed at 9/16, the ribstiffener spacing of two different studied sandwich structures is, respectively, selected as = = 200 mm and = = 250 mm. The STL curves are obtained as shown in Figure 13, wherein the group with spacing of 250 mm is a contrast structure, the lattice constant and width of piezoelectric patch are, respectively, 50 mm and 37.5 mm, and the remaining parameters are unchanged.
As can be seen from Figure 13, the two STL curves have similar trends, but the first valley frequency of the contrast structure is lower than that of the research structure; the exact values are 87.0 Hz and 133.1 Hz, respectively. Therefore, increasing the ribstiffener spacing improves the STL and broadens sound transmission bandwidth to a certain extent. This is because the natural frequency decreases as the ribstiffener spacing increases, and the modal amplitude coefficient at the resonant frequency increases rapidly.
4.3.3. Effect of Facesheet Thickness
For the research structure, the increase of total mass to some extent can improve the STL performance, but the location of the valley frequencies limit the sound insulation and bandwidth of the structure, thereby reducing its engineering application value. Therefore, in order to minimize the total mass of the structure, the rational design of structural geometry to improve its sound insulation performance is critical. To quantify the effects of facesheet thickness, the three STL versus frequency curves of different studied sandwich structures under the condition of the thicknesses (, ), respectively, selected as (1.0, 1.0), (1.5, 1.5), and (2.0, 2.0) mm, are presented in Figure 14.
Figure 14 indicates that, with the increase of the facesheet thicknesses (, ), the first valley frequency moves to a higher frequency, and the second valley frequency moves to a lower frequency. That is, the sound insulation bandwidth gradually narrows, because the equivalent surface mass density and equivalent dynamic bending stiffness of the facesheet change as the thicknesses increase. For the three studied models with (1.0, 1.0), (1.5, 1.5), and (2.0, 2.0) mm, the natural frequency is 57.4 Hz, 91.4 Hz, and 133.1 Hz, which are generated by the (1, 1)order modal of each model. And the second valley values were 1312.6 Hz, 1266.6 H, and 1146.9 Hz, respectively, by the role of the (5, 7), (4, 2), and (3, 0)order modal. In addition, the increments of sound insulation caused by the equivalent resonators tend to deteriorate when the thicknesses (, ) decrease, since the multimode local resonance frequencies are close to other modal frequencies between the first valley frequency and the second valley frequency of the studied sandwich structure.
4.3.4. Effect of External Shunt Circuit
In view of the equivalent negativecapacitance circuit shown in Figure 6, the contrast curves depicted in Figure 15 are the STLs corresponding to the negative capacitance of −0.1C_{p}, −0.3C_{p}, and −0.5C_{p}, respectively. This figure shows that, as the negative capacitance increases, the second valley frequency changes slightly, with the exact values being 1143.2 Hz, 1146.9 Hz, and 1155.0 Hz, respectively. In addition, the change in negative capacitance has no effect on the STL of the research sandwich structure in the range of less than about 1000 Hz, but beyond this frequency, the amplitudes of the STLs and all extreme frequencies are affected. This is because the negative capacitance directly affects the equivalent Young’s modulus and equivalent Poisson’s ratio of the facesheets in the research structure.
5. Conclusions
This paper is related to the lowfrequency STL of an infinite orthogonally ribstiffened sandwich structure, with the inner sides of the upper and bottom facesheet flexibly connected with periodic subwavelength arrays of shunted piezoelectric patches. For exploring its lowfrequency vibroacoustic characteristics, a complete theoretical model of STL for solving the research sandwich structure can be built in three steps. The first step is that the upper and bottom facesheet with periodic subwavelength arrays of shunted piezoelectric patches are equivalent to two homogeneous facesheets in a periodic element by effective medium (EM) method. The second step is to establish the dynamic motion governing equations by space harmonic method, and the tensional forces, bending moments, and torsional moments of the ribstiffeners as well as the inertial terms are taken into account; then the amplitude coefficients of the equivalent facesheets are separated by the virtual work principle under the velocity continuity conditions. The third step is to deduce the calculation formula of STL. But most importantly, the application of the theoretical model has two prerequisites. One is that the unit cells in facesheets satisfy the subwavelength assumption; that is, the lattice constant is smaller than the flexural wavelength of the substrate. The other is that the infinite sets of coupled equations must be truncated to be finite by the convergence criterion; that is, at the highest frequency, the difference between the STLs corresponding to two adjacent terms is within the preset error band (0.008 dB in this work).
To verify the validity and feasibility of the theoretical model, the finite element simulation of the orthogonally ribstiffened sandwich structure is carried out. The corresponding valley frequencies and peak frequencies can effectively match. Then the complete model is simplified to a typically infinite sandwich structure with unidirectionally parallel ribstiffeners and compares with the experimental results in the frequency range 50 Hz–6300 Hz of the experimental specimen numbered G16_SS90(610)_G16 in the internal report of the NRCC. The theoretical values are in good agreement with the experimental values, except for a certain discrepancy in narrowband 200 Hz–300 Hz and in the vicinity of 1000 Hz. Furthermore, the STL of the complete model is compared with that of the typical orthogonally ribstiffened sandwich structure, which proves that the research sandwich structure has good sound insulation performance.
At the end of the work, the influences of the parameters such as the elevation angle, ribstiffener spacing, facesheet thicknesses, and negative capacitance on the research structure are studied, and four main conclusions are drawn. First, the STL decreases with increasing elevation angle, and it has no effect on the bandwidth. Second, appropriately increasing the ribstiffener spacing can improve the sound insulation and broaden the bandwidth to a certain extent. Third, as the thickness of the facesheets increases from 1 mm to 2 mm, the sound insulation bandwidth gradually narrows, and the increments of sound insulation caused by the equivalent resonators become better. Finally, the change of negative capacitance has a slight influence on the bandwidth of STL, but it has a certain effect on the amplitude of STL above the second valley frequency.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (nos. 51575201 and 51405093), Guangxi Key Laboratory of Manufacturing Systems and Advanced Manufacturing Technology, China (no. 1404515006Z), and Natural Science Foundation of Guangxi, China (no. 2015GXNSFBA139232).
References
 F. J. Fahy and E. Lindqvist, “Wave propagation in damped, stiffened structures characteristic of ship construction,” Journal of Sound and Vibration, vol. 45, no. 1, pp. 115–138, 1976. View at: Publisher Site  Google Scholar
 B. R. Mace, “Periodically stiffened fluidloaded plates, I: response to convected harmonic pressure and free wave propagation,” Journal of Sound and Vibration, vol. 73, no. 4, pp. 473–486, 1980. View at: Publisher Site  Google Scholar
 B. R. Mace, “Periodically stiffened fluidloaded plates, II: response to line and point forces,” Journal of Sound and Vibration, vol. 73, no. 4, pp. 487–504, 1980. View at: Publisher Site  Google Scholar
 B. R. Mace, “The vibration of plates on twodimensionally periodic point supports,” Journal of Sound and Vibration, vol. 192, no. 3, pp. 629–643, 1996. View at: Publisher Site  Google Scholar
 D. J. Mead, “Wave propagation in continuous periodic structures: research contributions from Southampton, 1964–1995,” Journal of Sound and Vibration, vol. 190, no. 3, pp. 495–524, 1996. View at: Publisher Site  Google Scholar
 J.H. Lee and J. Kim, “Analysis of sound transmission through periodically stiffened panels by spaceharmonic expansion method,” Journal of Sound and Vibration, vol. 251, no. 2, pp. 349–366, 2002. View at: Publisher Site  Google Scholar
 J. Wang, T. J. Lu, J. Woodhouse, R. S. Langley, and J. Evans, “Sound transmission through lightweight doubleleaf partitions: theoretical modelling,” Journal of Sound and Vibration, vol. 286, no. 45, pp. 817–847, 2005. View at: Publisher Site  Google Scholar
 J. Legault and N. Atalla, “Numerical and experimental investigation of the effect of structural links on the sound transmission of a lightweight double panel structure,” Journal of Sound and Vibration, vol. 324, no. 35, pp. 712–732, 2009. View at: Publisher Site  Google Scholar
 L. Dozio and M. Ricciardi, “Free vibration analysis of ribbed plates by a combined analyticalnumerical method,” Journal of Sound and Vibration, vol. 319, no. 12, pp. 681–697, 2009. View at: Publisher Site  Google Scholar
 F. X. Xin and T. J. Lu, “Analytical modeling of fluid loaded orthogonally ribstiffened sandwich structures: sound transmission,” Journal of the Mechanics and Physics of Solids, vol. 58, no. 9, pp. 1374–1396, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 F. X. Xin and T. J. Lu, “Analytical modeling of wave propagation in orthogonally ribstiffened sandwich structures: sound radiation,” Computers and Structures, vol. 89, no. 56, pp. 507–516, 2011. View at: Publisher Site  Google Scholar
 V. Hongisto, “Airborne sound insulation of wall structures—measurement and prediction methods,” Tech. Rep. 56, Laboratory of Acoustics and Audio Signal Processing, Helsinki University of Technology, Espoo, Finland, 2000. View at: Google Scholar
 Y. Hamada, Y. Hidaka, and H. Tachibana, “Estimation of random incidence sound transmission loss of multiple structure walls by fourterminal network theory,” Journal of Structural and Construction Engineering (Transactions of Architectural Institute of Japan), vol. 362, pp. 1–7, 1986 (Japanese). View at: Google Scholar
 B. H. Sharp, “Prediction methods for the sound transmission of building elements,” Noise Control Engineering, vol. 11, no. 2, pp. 53–63, 1978. View at: Publisher Site  Google Scholar
 Q. Gu and J. Wang, “Effect of resilient connection on sound transmission loss of metal stud double panel partitions,” Chinese Journal of Acoustics, vol. 2, pp. 113–126, 1983. View at: Google Scholar
 F. Fahy and P. Gardonio, Sound and Structural Vibration, Elsevier, Amsterdam, The Netherlands, 2nd edition, 2007.
 J. L. Davy, “The sound transmission of cavity walls due to studs,” in Proceedings of the INTERNOISE and NOISECON Congress and Conference (InterNoise '93), vol. 4, pp. 975–978, Leuven, Belgium, 1993. View at: Google Scholar
 M. L. Rumerman, “Vibration and wave propagation in ribbed plates,” Journal of the Acoustical Society of America, vol. 57, no. 2, pp. 370–373, 1975. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 G. F. Lin and J. A. M. Garrelick, “Sound transmission through periodically framed parallel plates,” Journal of the Acoustical Society of America, vol. 61, no. 4, pp. 1014–1018, 1977. View at: Publisher Site  Google Scholar
 D. Takahashi, “Sound radiation from periodically connected doubleplate structures,” Journal of Sound and Vibration, vol. 90, no. 4, pp. 541–557, 1983. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 L. Maxit, “Wavenumber space and physical space responses of a periodically ribbed plate to a point drive: a discrete approach,” Applied Acoustics, vol. 70, no. 4, pp. 563–578, 2009. View at: Publisher Site  Google Scholar
 D. J. Mead, “Free wave propagation in periodically supported, infinite beams,” Journal of Sound and Vibration, vol. 11, no. 2, pp. 181–197, 1970. View at: Publisher Site  Google Scholar
 D. J. Mead and K. K. Pujara, “Spaceharmonic analysis of periodically supported beams: response to convected random loading,” Journal of Sound and Vibration, vol. 14, no. 4, pp. 525–541, 1971. View at: Publisher Site  Google Scholar
 R. E. Halliwell, T. R. T. Nightingale, A. C. C. Warnock, and J. A. Birta, “Gypsum board walls: transmission loss data,” Internal Report IRCIR761, Institute for Research in Construction, National Research Council of Canada, Ottawa, Canada, 1998. View at: Google Scholar
 J. L. Davy, “Predicting the sound insulation of walls,” Building Acoustics, vol. 16, no. 1, pp. 1–20, 2009. View at: Publisher Site  Google Scholar
 J. L. Davy, “Sound transmission of cavity walls due to structure borne transmission via point and line connections,” Journal of the Acoustical Society of America, vol. 132, no. 2, pp. 814–821, 2012. View at: Publisher Site  Google Scholar
 V. Hongisto, “Sound insulation of doors—Part 1: prediction models for structural and leak transmission,” Journal of Sound and Vibration, vol. 230, no. 1, pp. 133–148, 2000. View at: Publisher Site  Google Scholar
 V. Hongisto, J. Keränen, and M. Lindgren, “Sound insulation of doors—Part 2: comparison between measurement results and predictions,” Journal of Sound and Vibration, vol. 230, no. 1, pp. 149–170, 2000. View at: Publisher Site  Google Scholar
 V. Hongisto, M. Lindgren, and J. Keränen, “Enhancing maximum measurable sound reduction index using sound intensity method and strong receiving room absorption,” Journal of the Acoustical Society of America, vol. 109, no. 1, pp. 254–265, 2001. View at: Publisher Site  Google Scholar
 V. Hongisto, “Case study of flanking transmission through double structures,” Applied Acoustics, vol. 62, no. 5, pp. 589–599, 2001. View at: Publisher Site  Google Scholar
 N. W. Hagood and A. von Flotow, “Damping of structural vibrations with piezoelectric materials and passive electrical networks,” Journal of Sound and Vibration, vol. 146, no. 2, pp. 243–268, 1991. View at: Publisher Site  Google Scholar
 J. J. Hollkamp, “Multimodal passive vibration suppression with piezoelectric materials and resonant shunts,” Journal of Intelligent Material Systems and Structures, vol. 5, no. 1, pp. 49–57, 1994. View at: Publisher Site  Google Scholar
 O. Thorp, M. Ruzzene, and A. Baz, “Attenuation and localization of wave propagation in rods with periodic shunted piezoelectric patches,” Smart Materials and Structures, vol. 10, no. 5, pp. 979–989, 2001. View at: Publisher Site  Google Scholar
 O. Thorp, M. Ruzzene, and A. Baz, “Attenuation of wave propagation in fluidloaded shells with periodic shunted piezoelectric rings,” Smart Materials and Structures, vol. 14, no. 4, pp. 594–604, 2005. View at: Publisher Site  Google Scholar
 S. Chen, G. Wang, J. Wen, and X. Wen, “Wave propagation and attenuation in plates with periodic arrays of shunted piezopatches,” Journal of Sound and Vibration, vol. 332, no. 6, pp. 1520–1532, 2013. View at: Publisher Site  Google Scholar
 Y. Xiao, J. Wen, and X. Wen, “Sound transmission loss of metamaterialbased thin plates with multiple subwavelength arrays of attached resonators,” Journal of Sound and Vibration, vol. 331, no. 25, pp. 5408–5423, 2012. View at: Publisher Site  Google Scholar
 H. Zhang, J. Wen, Y. Xiao, G. Wang, and X. Wen, “Sound transmission loss of metamaterial thin plates with periodic subwavelength arrays of shunted piezoelectric patches,” Journal of Sound and Vibration, vol. 343, pp. 104–120, 2015. View at: Publisher Site  Google Scholar
 B. R. Mace, “Sound radiation from fluid loaded orthogonally stiffened plates,” Journal of Sound and Vibration, vol. 79, no. 3, pp. 439–452, 1981. View at: Publisher Site  Google Scholar
 G. P. Mathur, B. N. Tran, J. S. Bolton, and N.M. Shiau, “Sound transmission through stiffened doublepanel structures lined with elastic porous materials,” in Proceedings of the 14th DGLR/AIAA Aeroacoustics Conference, pp. 102–105, Aachen, Germany, May 1992. View at: Google Scholar
 A. Spadoni and M. Ruzzene, “Structural and acoustic behavior of chiral trusscore beams,” Journal of Vibration and Acoustics, Transactions of the ASME, vol. 128, no. 5, pp. 616–626, 2006. View at: Publisher Site  Google Scholar
 F. X. Xin, T. J. Lu, and C. Q. Chen, “External mean flow influence on noise transmission through doubleleaf aeroelastic plates,” AIAA Journal, vol. 47, no. 8, pp. 1939–1951, 2009. View at: Publisher Site  Google Scholar
 E. E. Ungar and J. A. Zapfe, Structural Damping, Noise and Vibration Control Engineering: Principles and Applications, John Wiley & Sons, 2nd edition, 2007.
 E. G. Williams, Fourier Acoustics: Sound Radiation and Nearfield Acoustical Holography, Academic Press, 1999.
Copyright
Copyright © 2017 Zhifu Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.