This study includes the first work about the absorption of a panel absorber under the effects of microperforation, air pumping, and linear and nonlinear vibrations. In practice, thin perforated panel absorber is backed by a flexible wall to enhance the acoustic performance within the room. The panel is easily excited to vibrate nonlinearly and the wall can vibrate linearly. However, the assumptions of linear panel vibration and no wall vibration are adopted in many research works. The development of the absorption formula is based on the classical approach and the electroacoustic analogy, in which the impedances of microperforation, air pumping, and linear and nonlinear vibrations are in parallel and connected to that of the air cavity in series. Unlike those finite element, numerical integration, and multiscale solution methods and so forth, the analytic formula to calculate the absorption of a panel absorber does not require heavy computation effort and is suitable for engineering calculation purpose. The theoretical result obtained from the proposed method shows reasonable agreement with that from a previous numerical integration method. It can be concluded that the overall absorption bandwidth of a panel absorber with an appropriate configuration can be optimized and widened by making use of the positive effects of microperforation, air pumping, and panel vibration.
1. Introduction
Numerous studies of structural-acoustic interaction have been carried out in recent decades (e.g., [1–6]). Particularly, the topics of microperforated panel and panel absorber have attracted the interest of many researchers. It is because (1) microperforated panel and panel absorber can be made of any durable materials, (2) their maintenance is relatively simple, and (3) microperforated panel requires small space to achieve high sound absorption when compared with typical foam or porous materials. Acoustic theories for these two treatments were developed by Mulholland and Parbrook [7], Ford and McCormick [8], and Maa [9, 10]. According to the assumptions of the previous studies, the wall behind a panel absorber was considered as rigid panel (i.e., no vibration) and the panel absorber was considered to vibrate linearly (in fact, thin panels are easy to vibrate nonlinearly). Although many linear structural-acoustic or nonlinear structural vibration problems had been solved (e.g. [11–17]), very limited studies considered the nonlinear vibration effect on the sound absorption performance of a panel absorber. In fact, it is more appropriate to employ nonlinear approach to analyze the vibration of a thin panel. There are five of those limited nonlinear structural-acoustic research works [18–22], which focused on the natural frequencies of nonlinear structural-acoustic systems, and the sound absorption, radiation, and transmission of nonlinear flat/curved panels. These works are relevant to the nonlinear structural-acoustic problem handled in this paper, but the researchers adopted the solution methods which required heavy computational efforts (e.g., finite element method, numerical integration method, and harmonic balance method). Other solution methods for large amplitude structural vibrations and nonlinear oscillations (e.g., the method of multiple scales, the method of normal forms, the method of Shaw and Pierre, and the method of King and Vakakis [23–27]) also require relatively tedious solution implementation and heavy computational efforts. In [28], the simple and straightforward electroacoustic analogy was adopted for modelling of the sound absorption of a panel absorber. Besides, the air pumping effect (or Helmholtz resonance effect) has been mentioned in many textbooks and articles (e.g., [29]), but it has not been incorporated into the research problem of nonlinear panel absorber. In this paper, it focuses on the development of the analytic absorption formula of a panel absorber which considers the four effects (i.e., microperforation, air pumping, and linear and nonlinear vibrations).
2. Analytic Formulation2.1. Nonlinear Panel Vibration and Liner Wall Vibration
Figure 1 shows the proposed panel absorber subject to the effects of microperforation, air pumping, and linear and nonlinear vibrations. Firstly, the impedances of the nonlinear panel and linear wall are derived in this subsection. According to the formulation of Chu and Herrmann [30], the governing equation for the nonlinear vibration of a rectangular panel subject to uniformly distributed harmonic excitation can be expressed in the following form (see Appendix A for more):
(1)ρd2εdτ2+ρωo2ε+βε3+F(τ)=0,
where
r=a/b is the aspect ratio; E is Young’s modulus; ρ is the surface density; ν is Poisson’s ratio; and t, a, and b are the thickness, length, and width.
Consider the approximation of ε=Acos(ωτ) in (1). Then
(4)34βA3+ρ(ωo2-ω2)A+Fo,mn=0,⇒(5)A3+43βρ(ωo2-ω2)A+43βFo,mn=0,
where A= modal vibration amplitude.
A panel absorber subject to the effects of microperforation, air pumping, and linear and nonlinear vibrations.
For simplicity, let p=(4/3β)ρ(ωo2-ω2),q=(4/3β)Fo,mn, and A=B-(p/3B) (i.e., Vieta’s substitution [31]). Equation (5) can be rewritten as
(6)B3+q-p327B3=0,⇒(7)B6+qB3-p327=0.
If B3 is considered as an independent unknown, (7) is a “modified” quadratic equation. Then, the analytic formulas for the cubic solutions are given below:(8a)A1=-1312(Θ+Ω)3-1312(Θ-Ω)3,(8b)A2=1+j3612(Θ+Ω)3+1-j3612(Θ-Ω)3,(8c)A3=1-j3612(Θ+Ω)3+1+j3612(Θ-Ω)3,where Θ=(27q)2+4(3p)3, Ω=27q.
In the next step, the nonlinear stiffness can be defined as (3/4)βA2 and overall panel stiffness is ρωo2+(3/4)βA2. The panel impedance of the (m,n) mode is given by
(9)ZP,mn=Fo,mnVP,mn=2ρξωωN+j[ρ(ω2-ωo2)-(3/4)βA2]ω,
where j=-1= complex number; ξ is the panel damping ratio. VP,mn is the modal panel velocity. Note that ωN is the nonlinear peak frequency, instead of the linear natural frequency, ωo in [32] (see Appendix B for the derivation of ωN). When it is the case of small amplitude vibration or linear vibration, ωN is equal to ωo.
Consider the overall panel velocity contributed from the modal velocities and find the overall panel impedance:
(10)VP=∑mM∑nNFoΛmnΛmn′ZP,mn,⇒(11)ZP=FoVP=(∑mM∑nNΛmnΛmn′ZP,mn)-1,
where Λmn′=(∫0b∫0asin(mπx/a)sin(nπy/b)dxdy)/ab. According to [33], the sound absorption of two linearly vibrating plates with different boundary conditions was very close. It was concluded that the boundary condition would not affect the sound absorption significantly. That is why only one boundary condition case is considered here.
Consider the linear case (or small amplitude case) in the model representing the panel vibration and the corresponding modal and overall panel impedances (i.e., (9), (11)). It is noted that, for the linear vibration, the fundamental harmonic balance can lead to the exact solution (no higher harmonic terms are required). Hence the modal and overall wall impedances are given below:
(12)ZW,mn=2ρWξWωωW+j[ρ(ω2-ωW2)]ω,(13)ZW=(∑mM∑nNΛmnΛmn′ZW,mn)-1,
where ξW,ρW, and ωW are the wall damping factor, density, and natural frequency
2.2. Air Pumping (Helmholtz Resonance)
In Figure 1, the impedance of the air pumping along the tube mounted on the panel absorber is obtained here. The differential equation for the air mass pumping in/out along the tube is given by [29]
(14)ρaL′d2Λdτ2+(Rr+Rw)SdΛdτ+P(τ)=0,
where Λ= air mass displacement, P(τ)= acoustic pressure excitation, Rr=ρacak2S2/4π= damping due to radiation resistance, Rw=2ρaL′Sαw= damping due to thermoviscous resistance, L′=L+1.4at= effective tube length, αw=(1/atca)(ηω/2ρa)0.5(1+(γ/Pr))= absorption coefficient for wall losses, and ρa is the air density; L′ is the effective length of the tube with unflanged outer end; k is the wave number; L is the tube length; at is the radius of the tube; S is the cross section area of the tube; η is the air viscosity; γ is the ratio of specific heat; Pr is Prandtl number.
Similar to the overall panel impedance in (9), the impedance of air pumping along the tube is given below:
(15)ZH′=((Rr+Rw)/S)ω+jρaL′ω2ω.
Using (15), the air mass velocity is given by
(16)ZH′=PoVH′⟹VH′=PoZH′,
where Po = acoustic pressure magnitude.
Consider the average velocity over the panel surface
(17)VH=VH′σH,
where σH= the ratio of tube area to panel area.
Then the overall impedance of air pumping is obtained by substituting (17) into (16):
(18)ZH=ZH′σH.
2.3. Microperforation
In Figure 1, the impedance of the microholes over the panel surface is obtained here. According to [9, 10], the real and imaginary parts of the acoustic impedance of a microhole are given below:(19a)ZM,R′=ρaca(0.147td2)(9+100d2f32+1.768fd2t),(19b)ZM,I′=ρaca1.847ft(1+19+50d2f+0.85dt),where ZM′=ZM,R′+jZM,I′= the impedance of a microhole, f is frequency, t is panel thickness, and d is the perforation diameter. ρa and ca are air density and sound speed.
Then, the overall impedance of all microholes over the panel is given by
(20)ZM=ZM′σ,
where σ is the perforation ratio.
2.4. Air Cavity
In Figure 1, the impedance of the air cavity is obtained here. The acoustic velocity potential within a rectangular cavity is given by the following homogeneous wave equation [1]:
(21)∇2ϕ-1ca2∂2ϕ∂τ2=0,
where ϕ is the velocity potential function.
The air particle velocities in the x, y, and z directions and pressures within the air cavity can be given by ∂ϕ/∂x,∂ϕ/∂y,∂ϕ/∂z, and -ρa(∂ϕ/∂τ), respectively.
The boundary conditions of the rectangular cavity to be satisfied are(22a)∂ϕ∂x|x=0=∂ϕ∂x|x=a=0;(22b)∂ϕ∂y|y=0=∂ϕ∂y|y=b=0,(22c)∂ϕ∂z|z=0=0;(22d)∂ϕ∂z|z=-D=VD,where a and b are the length and width of the rectangular cavity or the panel and VD is the average velocity at z=D= cavity depth. Note that the cavity impedance of no wall vibration is derived here and used into the sound absorption formulation in the next subsection which can consider both the panel and wall vibration effects and is derived using the electroacoustic analogy.
By applying the boundary conditions in (22a) and (22b), the solution of (21) can be expressed as
(23)ϕ=∑u=0U∑w=0W[Luwcosh(μuwz)+Nuwsinh(μuwz)]×cos(uπxa)cos(wπyb)eiωτ,
where μuw=((uπ/a)2+(wπ/b)2)-(ω/ca)2, u and w are the acoustic mode numbers, U and W are the numbers of acoustic modes used in the x and y directions, and Luw and Nuw are coefficients to be determined in the following step.
Substituting (23) into (22c) and (22d) gives(24a)∂ϕ∂z|z=0=0⟹Nuw=0,(24b)∂ϕ∂z|z=D=VDeiωτ⟹Luw=-∫0a∫0bVDcos(uπx/a)cos(wπy/b)dxdyκuwμuwsinh(μuwD),where κuw=∫0a∫0bcos(uπx/a)2cos(wπy/b)2dxdy.
By inputting (24a) and (24b) into (23), the pressure at z=D is obtained and given below:
(25)FD=∑u=0U∑w=0W[jρaωcoth(μuwD)κuwμuwaaaaaaaaa×coth(μuwD)κuwμuw∫0a∫0bVDcos(uπxa)cos(wπyb)dxdy]×cos(uπxa)cos(wπyb).
The average pressure force at z=D is then given by taking integration over the panel area
(26)F¯D=∫0a∫0bFDcos(uπx/a)cos(wπy/b)dxdyab=VDjρaω∑u=0U∑w=0W[κuw′κuwcoth(μuwD)μuw],
where κuw′=(∫0a∫0bcos(uπx/a)cos(wπy/b)dxdy)2/ab.
Finally, the impedance of the air cavity is given by
(27)Zcav=F¯DVD=jρaω∑u=0U∑w=0W[κuw′κuwcoth(μuwD)μuw].
Note that since even pressure distribution is considered, then u,w=0,2,4,…. If a and b are less than 0.3 m, then the frequencies of the fundamental resonances parallel to the cavity length or width are higher than 1000 Hz (note that 1000 Hz is not within the frequency range concerned in this study). Then, in this case, only the mode numbers, u=0 and w=0, are considered. The impedance can be approximated by
(28)Zcav≈-jρaωcot(ωcaD).
2.5. Electroacoustic Analogy
According to [34], the acoustic impedance of a double panel absorber can be analogous to that of the electrical circuit. The acoustic impedance is given below:
(29)Zdou=Zp,1+Zcav,1+(Zcav,1)2Zp,2+Zcav,1+Zcav,2,
where Zp,1 and Zp,2= acoustic impedance of the 1st and 2nd panels; Zcav,1 and Zcav,2= acoustic impedance of the 1st and 2nd cavities.
If the impedance of the 2nd cavity is set to zero, that is, equivalent to the case of a panel absorber backed by a flexible wall, the overall impedance in (29) is rewritten as
(30)Zo=Zcom+Zcav+(Zcav)2ZW+Zcav,
where Zcom=(ZMZPZH)/(ZMZP+ZPZH+ZHZM)= combined impedance. ZM= impedance of the microperforation; ZP= acoustic impedance of the panel absorber; ZH= impedance of the air pumping; ZW= impedance of the wall.
If the wall is perfectly rigid (i.e., no vibration), the natural frequency of wall ωW in (12)→∞. Hence, ZW,mn, ZW→∞ and (Zcav)2/(ZW+Zcav)≈0. Equation (30) can be rewritten as
(31)Zo=Zcom+Zcav.
Then the absorption coefficient of the absorber αo for normal incidence can be calculated by the well-known formula (see [29]):
(32)αo=4Re(Zo)(1+Re(Zo))2+(Im(Zo))2.
3. Results and Discussions
Using (32), the sound absorption is obtained for the cases considered in this section. Each panel absorber is made of aluminum panel. The material properties are Young’s modulus E=7×1010N/m2, Poisson’s ratio ν=0.3, and mass density ρ=2700kg/m3. The sound speed and air density are 340 m/s and 1.2 kg/m3, respectively. Air viscosity η=1.85×10-5Pa·s; the ratio of specific heat γ=1.402; Prandtl number Pr=0.71. Various excitation levels, damping factors, panel sizes, and so forth are considered and their effects on the sound absorption performances are studied in the following paragraphs.
Figures 2 and 3 present comparisons between the frequency-absorption curves obtained from the proposed method and numerical integration method [32]. The results show that there is an absorption peak due to the microperforation around 178 Hz. The well-known “jump phenomenon” can be seen on both curves. The results obtained from the two methods are generally in reasonably good agreement, although some deviations are observed at the absorption peaks around 422 Hz in Figure 2 and 390 Hz in Figure 3. Besides, small differences between the jump frequencies (i.e., 464 and 470 Hz in Figure 2 and 410 and 420 Hz in Figure 3) are also found. The deviations are caused by that the damping term in the proposed formulation is proportional to the nonlinear peak frequency ωN, instead of the linear natural frequency ωo, which was adopted in [32]. The excitation level in Figure 3 is lower than that in Figure 2. Therefore, the structural vibration in Figure 3 is less nonlinear (i.e., ωN is closer to ωo) and the deviation between the peak values from the two methods is smaller. Note that in [32], the numerical integration method did not deliver the unstable solutions. Therefore, only stable results obtained by the proposed method in Figure 3 are compared with the numerical results in [32]. Besides, the unstable results are not obtainable in practical environments and not very useful for practical engineering purpose.
Comparison between the results from the proposed classical and numerical methods (t=0.8 mm, a=b=0.1 m, D=250 mm, ξ=0.01, σ=0.5%, d=0.4 mm, and κ=5.8).
Comparison between the results from the proposed classical and numerical methods (t=0.8 mm, a=b=0.1 m, D=250 mm, ξ=0.01, σ=0.5%, d=0.4 mm, and κ=2.9).
Figures 4(a)–4(c) present the sound absorption plotted against the excitation frequency for various panel resonant frequencies. In each of these three figures, there are two nonlinear solution types. The “jump up” or “jump down” solution is obtained by sweeping the excitation frequency from “high to low” or “low to high”. In Figures 4(a)–4(c), the linear (1,1) mode panel resonant frequencies of the three cases are 387 Hz, 172 Hz, and 97 Hz, respectively. The absorption peak frequency due to the microperforation effect is 178 Hz. It can be seen that (1) if the panel resonant frequency is higher than the absorption peak frequency due to the microperforation effect, the panel vibration can enhance and widen the absorption bandwidth and (2) if the panel resonant frequency is close to or lower than the absorption peak frequency, the panel vibration degrades the absorption performance and the peak absorption is deteriorated.
(a) The effect of the (1,1) mode vibration (ωo= 387.25 × 2π Hz, t=0.8 mm, a=b=0.1 m, D=250 mm, ξ=0.01, σ=0.5%, d=0.4 mm, and κ=10). (b) The effect of the (1,1) mode vibration (ωo= 172.11 × 2π Hz, t=0.8 mm, a=b=0.15 m, D=250 mm, ξ=0.01, σ=0.5%, d=0.4 mm, and κ=10). (c) The effect of the (1,1) mode vibration (ωo= 96.81 × 2π Hz, t=0.8 mm, a=b=0.2 m, D=250 mm, ξ= 0.01, σ=0.5%, d=0.4 mm, and κ=10).
Figures 5(a) and 5(b) present the sound absorption plotted against the excitation frequency for various excitation levels. The “jump down” solutions are shown in Figure 5(a) (i.e., sweeping the excitation frequency from “high to low”) while the “jump up” solutions are in Figure 5(b) (i.e., sweeping the excitation frequency from “low to high”). The absorption peak around 178 Hz is due to the microperforation effect and is not significantly affected by the excitation level. It is noted that in the cases of κ=1, there is no jump phenomenon as the excitation force is quite small (i.e., linear case). Besides, it can be seen that (1) in Figure 5(a) the higher the excitation force, the wider the absorption peak due to the panel vibration and (2) in Figure 5(b) if the panel vibration is linear (or low excitation level), the absorption peak is very narrow and not very useful for widening the absorption bandwidth. Figures 6(a) and 6(b) present the sound absorption plotted against the excitation frequency for various damping ratios. The “jump down” solutions are shown in Figure 6(a) while the “jump up” solutions are in Figure 6(b). The absorption peak due to the microperforation effect is not significantly affected by the panel damping ratio. Besides, it can be seen that, (1) in Figure 6(a), the higher the damping ratio, the narrower the absorption peak due to the “jump down” panel vibration and the lower the “jump down” frequency and, (2) in Figure 6(b), the absorption peak due to the “jump up” panel vibration is not significantly affected by the damping ratio.
(a) The effect of excitation level (low to high frequency excitation, ωo= 387.25 × 2π Hz, t=0.8 mm, a=b=0.1 m, D=250 mm, ξ=0.01, σ=0.5%, and d=0.4 mm). (b) The effect of excitation level (high to low frequency excitation, ωo=387.25×2πHz, t=0.8 mm, a=b=0.1 m, D=250 mm, ξ=0.01, σ=0.5%, and d=0.4 mm).
(a) The effect of damping ratio (low to high frequency excitation, ωo= 387.25 × 2π Hz, t=0.8 mm, a=b=0.1 m, D=250 mm, σ=0.5%, d=0.4 mm, and κ=10). (b) The effect of damping ratio (high to low frequency excitation, ωo= 387.25 × 2π Hz, t=0.8 mm, a=b=0.1 m, D=250 mm, σ=0.5%, d=0.4 mm, and κ=10).
Figures 7(a) and 7(b) present the sound absorption and vibration amplitude plotted against the excitation frequency. Note that this case does not consider any perforation effect and other effects. There are two peaks in each of the two figures, which are due to the (1,1) and (1,3) mode resonances, respectively. It can be seen that the peaks due to the (1,1) mode resonance are more nonlinear. As the two peaks are far from each other, no strong mode coupling effect is induced. In Figure 7(b) the vibration amplitude is almost linearly increasing from 100 Hz to 286 Hz while the sound absorption is almost linearly increasing from 100 Hz to 150 Hz only and slowly and nonlinearly increasing from 150 Hz to 286 Hz. Besides, the absorption peak due to the (1,3) mode resonance is much more significant than the vibration peak. It can be implied that higher mode responses are important for absorption but not for vibration.
(a) The absorption peaks due to the (1,1) and (1,3) mode vibrations (low to high frequency excitation, t=0.5 mm, a=b=0.15 m, D=250 mm, and k=5). (b) The vibration peaks due to the (1,1) and (1,3) mode vibrations (low to high frequency excitation, t=0.5 mm, a=b=0.15 m, D=250 mm, and κ=5).
Figures 8(a) and 8(b) present the sound absorption plotted against the excitation frequency and show the effects of air pumping (Helmholtz resonance) and microperforation. In Figure 8(a), the peak frequencies of air pumping and microperforation are around 100 Hz. It can be seen that the absorption bandwidth of the case with both effects is much wider and its peak absorption is much higher than the cases with air pumping or microperforation only. The two effects can interact with each other and significantly enhance the absorption performance. In Figure 8(b), the peak frequencies of air pumping and microperforation are quite far from each other. The enhancement of the absorption performance due to the two effects is not as good as that in Figure 8(a).
(a) The effects of the air pumping and microperforation (t=0.3 mm, a=b=0.1 m, D=150 mm, σ=0.5%, d=0.3 mm, L=10 mm, and at=5 mm). (b) The effects of the air pumping and microperforation (t=0.3 mm, a=b=0.1 m, D=100 mm, σ=0.28%, d=0.3 mm, L=10 mm, and at=5 mm).
Figures 9(a) and 9(b) present the sound absorption plotted against the excitation frequency for various damping ratios and show the effect of linear wall vibration. The peak frequencies of wall panel and microperforation are around 80 Hz and 178 Hz, respectively. It can be seen that the wall vibration does not affect the absorption performance significantly and just induces an abrupt change around 80 Hz. The magnitude of the abrupt change depends on the damping ratio (i.e., the higher the damping ratio, the smaller the abrupt change). Figures 10(a) and 10(b) present the sound absorption plotted against the excitation frequency and show the effects of microperforation and linear and nonlinear vibrations. It can be seen that the effects of microperforation and nonlinear panel vibration are positive and significant while the effect of linear wall vibration is negative but not significant.
(a) The effects of linear wall vibration and microperforation (ωW=80×2πHz, ρw=240/m2, t=0.8 mm, a=b=0.1 m, D=250 mm, ξW=0.02, σ=0.5%, and d=0.4 mm). (b) The effects of linear wall vibration and microperforation (ωW=80×2πHz, ρw=240/m2, t=0.8 mm, a=b=0.1 m, D=250 mm, ξW=0.01, σ=0.5%, and d=0.4 mm).
(a) The effects of microperforation and linear and nonlinear vibrations (low to high frequency excitation, ωo=387.25×2πHz, ωW=80×2πHz, ρw=240kg/m2, t=0.8 mm, a=b=0.1 m, D=250 mm, ξ=ξW=0.01, σ=0.5%, d=0.4 mm, and κ=10). (b) The effects of microperforation and linear and nonlinear vibrations (high to low frequency excitation, ωo=387.25×2πHz, ωW=80×2πHz, ρw=240kg/m2, t=0.8 mm, a=b=0.1 m, D=250 mm, ξ=ξW=0.01, σ=0.5%, d=0.4 mm, and κ=10).
4. Conclusions
In this study, the analytic absorption formulation has been developed for the panel absorber under the effects of microperforation, air pumping, and linear and nonlinear vibration. From the predictions, it can be concluded that (1) if the panel resonant frequency is higher than the absorption peak frequency due to the microperforation effect, the panel vibration can enhance and widen the absorption bandwidth; (2) if the panel resonant frequency is close to or lower than the absorption peak frequency, the panel vibration degrades the absorption performance or the peak absorption is deteriorated; (3) the higher the excitation force, the wider the absorption peak due to the panel vibration; (4) if the panel vibration is linear (or low excitation level), the absorption peak is very narrow and not very useful for widening the absorption bandwidth; (5) the absorption bandwidth of the case with the air pumping and microperforation effects is much wider and its peak absorption is much higher than the cases with any one of these two effects. If the peak frequencies of air pumping and microperforation are far from each other, the enhancement of the absorption performance is small; and (6) the wall vibration does not affect the absorption performance significantly and just induces a small abrupt change on the absorption bandwidth around the wall resonant frequency.
AppendicesA. Derivation of the Nonlinear Differential Equation
The explanation of how the formulation of a nonlinear panel is expressed into the form of Duffing equation, which is abstracted from [30], is shown here. The flexible panel vibration is governed by the Von-Karman plate theory. The well-known partial differential equation of the nonlinear panel can be presented in terms of its displacement and the Airy stress function as follows:
(A.1)112(d2wdξ2+∇4w)=(∂2F¯∂η2∂2w∂ξ2+∂2F¯∂ξ2∂2w∂η2-2∂2F¯∂η∂ξ∂2w∂η∂ξ).
Note that (A.1) is abstracted [30] and the notations are the same as those in [30] but the definitions of some notations are different from those in the main text (see the following definitions). w= panel displacement; ξ and η are the dimensionless space variables of the x and y directions; ξ is the dimensionless time variable; ν is Poisson’s ratio. F¯ is the Airy stress function. According to [30], the derivatives of F¯ can be expressed in terms of the in-plane displacements along the x and y directions (i.e., u and v) and transverse displacement (i.e., w).Consider
(A.2a)∂2F¯∂ξ2=∂v∂η+12(∂w∂η)2+ν[∂u∂ξ+12(∂w∂ξ)2],(A.2b)∂2F¯∂η2=∂u∂ξ+12(∂w∂ξ)2+ν[∂v∂η+12(∂w∂η)2],(A.2c)-∂2F¯∂ξ∂η=1-ν2[∂u∂η+∂v∂ξ+∂w∂η∂w∂ξ].As it is assumed that the panel is simply supported, the following displacement functions which can satisfy the boundary conditions are adopted:(A.3a)w=εsin(mπξ)sin(nπrη),(A.3b)u=ε2π16(cos(2πnrη)-1+νn2r2m2)msin(2πmξ),(A.3c)v=ε2π16(nrcos(2πmξ)-nr+νm2nr)sin(2πnrη).Again, (A.3a)–(A.3c) originate from [30] and the definitions of some notations are different from those in the main text (see the following definitions). m and n are the modal number. ε is the modal amplitude. Consider r=a/b.
Then, consider putting (A.2a)–(A.2c) and (A.3a)–(A.3c) into the right terms of (A.1):(A.4a)∂2F¯∂ξ2∂2w∂η2=-∂2F¯∂ξ2π2n2r2Msin(mπξ)sin(πnrη)=-ε3π48n4r4(cos(2πmξ)-1+νm2n2r2)×cos(2πnrη)sin(mπξ)sin(πnrη)-12ε3π4n4r4sin2(πmξ)×cos2(πnrη)sin(mπξ)sin(πnrη)-νε3π48n2r2m2(cos(2πrη)-1+νn2r2m2)×cos(2πmξ)sin(mπξ)sin(πnrη)-ν12ε3π4m2n2r2×cos2(πmξ)sin2(πnrη)sin(mπξ)sin(πnrη),(A.4b)∂2F¯∂η2∂2w∂ξ2=-∂2F¯∂η2m2π2εsin(πmξ)sin(πnrη)=-ε3π48m4(cos(2πrη)-1+νn2r2m2)×cos(2πmξ)sin(πmξ)sin(πnrη)-12ε3π4m4×cos2(πmξ)sin2(πnrη)sin(πmξ)sin(πnrη)-νε3π48n2r2m2(cos(2πmξ)-1+νm2n2r2)×cos(2πnrη)sin(πmξ)sin(πnrη)-ν12ε3π4m2n2r2sin2(πmξ)×cos2(πnrη)sin(πmξ)sin(πnrη).(A.4c)-∂2F¯∂ξ∂η=1-ν2[∂u∂η+∂v∂ξ+∂w∂η∂w∂ξ]=1-ν2[-ε2mnrπ24sin(2πnrη)sin(2πmξ)aaaaaaaa+ε2mnrπ2sin(πnrη)aaaaaaaad7dεd2mnrπ24×cos(πnrη)sin(πmξ)cos(πmξ)]=0.Hence, (A.1) can be rewritten in the following form:
(A.5)[C1d2εdξ2+C2ε+C3ε3]sin(πmξ)sin(πnrη)+highermodeterms=0,
where C1, C2, and C3 are the constants dependent of m, n, r, ν, and the material and dimension parameters and so forth. According to [30], the higher mode terms in (A.5) are neglected. For forced vibration cases, a forcing term P is considered. That is,
(A.6)C1d2εdξ2+C2ε+C3ε3+P=0,
where (A.6) is a Duffing equation which form is the same as that of (1).
B. Derivation of the Peak Frequency
Consider a damping term and rewrite (1):
(B.1)ρd2εdτ2+2ξρωodεdτ+ρωo2ε+βε3+Fo,mnsin(ωτ)=0.
Take the approximation of ε=Acoscos(ωt)+Asinsin(ωt) and consider the harmonic balance of sin(ωt) and cos(ωt): (B.2a)-ρω2Asin-2ξρωoωAcos+ρωo2Asin+34βAsin3+34βAsinAcos2+Fo,mn=0,(B.2b)-ρω2Acos+2ξρωoωAsin+ρωo2Acos+34βcAcos3+34βcAsin2Acos=0. Let Ao2=Acos2+Asin2. Multiplying (B.2a) and (B.2b) by Acos and Asin, respectively, and then subtracting the new (B.2a) by the new (B.2b) yield(B.3a)-2ξρωoωAo2Fo,mn=Acos,(B.3b)Ao1-(2ξρωoωAoFo,mn)2=Asin.Then inputting (B.3a)-(B.3b) into (B.2b) yields
(B.4)ρ(ω2-ωo2)2ξρωoωAo2Fo,mn+2ξρωoωAo1-(2ξρωoωAoFo,mn)2-34β(Ao)42ξρωoωFo,mn=0.
Consider the equation representing the backbone curve by taking the forcing and damping terms out of (B.1):
(B.5)ρ(-ω2+ωo2)Ao+34βAo3=0,⇒(B.6)Ao2=ρ(ω2-ωo2)(3/4)β.
Inputting (B.6) into (B.4) yields
(B.7)1-(2ξρωoωFo,mn)2ρ(ω2-ωo2)(3/4)β=0,
where (B.7) contains only one unknown, ω. It is obtainable and represents the frequency at which the vibration amplitude peaks and is equal to ωN in (9).
Conflict of Interests
The author declares that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The work in this paper was fully supported by the CityU grant “Nonlinear Dropped Ceiling for Noise and Vibration Reduction, Project no. 7004039.” The author would like to express the sincere gratitude and appreciation to Dr. WY Poon for his advice about this research.
LeeY. Y.NgC. F.Sound insertion loss of stiffened enclosure plates using the finite element method and the classical approach1998217223926010.1006/jsvi.1998.17482-s2.0-0032183976LeeY. Y.The response frequency conversion characteristic of a nonlinear curved panel with a centre mass and the sound radiations2012201211298413MR298383310.1155/2012/298413HsiaS. Y.ChouY. T.Traffic noise propagating from vibration of railway wagon20132013765179710.1155/2013/651797FengJ.ZhengX. P.WangH. T.ZouY. J.LiuY. H.YaoZ. H.Low-frequency acoustic-structure analysis using coupled FEM-BEM method201320138583079MR312158210.1155/2013/583079ChenS. M.WangD. F.ZanJ. M.Interior noise prediction of the automobile based on hybrid FE-SEA method201120112032717010.1155/2011/3271702-s2.0-81555219665CrasterR. V.SmithS. G. L.A class of expansion functions for finite elastic plates in structural acoustics199910663128313410.1121/1.4281662-s2.0-0040313664MulhollandK. A.ParbrookH. D.Transmission of sound through apertures of negligible thickness19675349950810.1016/0022-460X(67)90194-02-s2.0-22344446941FordR. D.McCormickM. A.Panel sound absorbers196910341142310.1016/0022-460X(69)90219-32-s2.0-0014599658MaaD.Microperforated-panel wideband absorbers1987293778410.3397/1.28276942-s2.0-0023452729MaaD.Potential of microperforated panel absorber199810452861286610.1121/1.4238702-s2.0-0031765146WuY.HanG.Infinitely many sign-changing solutions for some nonlinear fourth-order beam equations2013201311635265MR305597110.1155/2013/635265LiuX.HuoB.ZhangS.Nonlinear dynamic analysis on the rain-wind-induced vibration of cable considering the equilibrium position of rivulet201320131092763210.1155/2013/927632MR3143559LinX.LiF.Asymptotic energy estimates for nonlinear Petrovsky plate model subject to viscoelastic damping2012201225419717MR299100310.1155/2012/419717Esquivel-AvilaJ. A.Dynamic analysis of a nonlinear Timoshenko equation20112011 3672481510.1155/2011/724815MR28028462-s2.0-80052671790SantosM. L.FerreiraJ.RaposoC. A.Existence and uniform decay for a nonlinear beam equation with nonlinearity of Kirchhoff type in domains with moving boundary2005890191910.1155/AAA.2005.901MR2201923ZBL1092.350682-s2.0-33749832856HeX.Modal decoupling using the method of weighted residuals for the nonlinear elastic dynamics of a clamped laminated composite200920091997293010.1155/2009/972930MR25300742-s2.0-68949122486WangD.ZhangJ.WangY.Strong attractor of beam equation with structural damping and nonlinear damping201320138769514MR309097310.1155/2013/769514LeeY. Y.Structural-acoustic coupling effect on the nonlinear natural frequency of a rectangular box with one flexible plate200263111157117510.1016/S0003-682X(02)00033-62-s2.0-0036833099LeeY. Y.Analysis of the nonlinear structural-acoustic resonant frequencies of a rectangular tube with a flexible end using harmonic balance and homotopy perturbation methods201220121339158410.1155/2012/3915842-s2.0-84871447341HuiC. K.LeeY. Y.ReddyJ. N.Approximate elliptical integral solution for the large amplitude free vibration of a rectangular single mode plate backed by a multi-acoustic mode cavity20114991191119410.1016/j.tws.2011.03.0022-s2.0-79959327417LeeY. Y.LiQ. S.LeungA. Y. T.SuR. K. L.The jump phenomenon effect on the sound absorption of a nonlinear panel absorber and sound transmission loss of a nonlinear panel backed by a cavity2012691-29911610.1007/s11071-011-0249-2MR29298582-s2.0-84861734877LeeY. Y.SuR. K. L.NgC. F.HuiC. K.The effect of modal energy transfer on the sound radiation and vibration of a curved panel: theory and experiment20093243–51003101510.1016/j.jsv.2009.02.0422-s2.0-67349204803MamandiA.KargarnovinM. H.FarsiS.Dynamic analysis of a simply supported beam resting on a nonlinear elastic foundation under compressive axial load using nonlinear normal modes techniques under three-to-one internal resonance condition20127021147117210.1007/s11071-012-0520-1MR29921232-s2.0-84871010630NayfehA. H.20112ndWeinheim, GermanyWiley-VCH10.1002/9783527635801MR2975979GhayeshM. H.BalarS.Non-linear parametric vibration and stability of axially moving visco-elastic Rayleigh beams20084525-266451646710.1016/j.ijsolstr.2008.08.0022-s2.0-54149118730ShawS. W.PierreC.Normal modes of vibration for nonlinear continuous systems1994169331934710.1006/jsvi.1994.1021MR12622642-s2.0-0028203123KingM. E.VakakisA. F.Energy-based formulation for computing nonlinear normal modes in undamped continuous systems1994116333234010.1115/1.29304332-s2.0-0028461550KangJ.FuchsH. V.Predicting the absorption of open weave textiles and micro-perforated membranes backed by an air space1999220590592010.1006/jsvi.1998.19772-s2.0-0000782333KinslerL. E.FreyA. R.CoppensA. B.SandersJ. V.19994thJohn Wiley & SonsChuH. N.HerrmannG.Influence of large amplitudes on free flexural vibrations of rectangular elastic plates195623532540MR0085015Weisstein EricW.Vieta's SubstitutionFrom MathWorld—A Wolfram Web Resource, http://mathworld.wolfram.com/VietasSubstitution.htmlLeeY. Y.GuoX.LeeE. W. M.Effect of the large amplitude vibration of a finite flexible micro-perforated panel absorber on sound absorption20078141442-s2.0-34250023426LeeY. Y.LeeE. W. M.NgC. F.Sound absorption of a finite flexible micro-perforated panel backed by an air cavity20052871-222724310.1016/j.jsv.2004.11.0242-s2.0-22344456634LeeY. Y.LeeE. W. M.Widening the sound absorption bandwidths of flexible micro-perforated curved absorbers using structural and acoustic resonances200749892593410.1016/j.ijmecsci.2007.01.0082-s2.0-34250870174