# Study on Acoustic Characteristics of a Flexible Plate Strongly Coupled with Rectangular Cavity

Authors: Jinlong Liao; Haichao Zhu; Jiuxiao Hou; Suwei Yuan
Journal: Shock and Vibration
 (2021)

---

## Abstract

This paper presents a method to predict the acoustic characteristics and steady-state responses of a flexible plate strongly coupled with rectangular cavity based on energy principle theory and Legendre polynomial series. First, the displacement of the plate and the sound pressure in the cavity are constructed in the form of two-dimensional and three-dimensional Legendre polynomial series, respectively. The unknown expansion coefficients are obtained using the Rayleigh–Ritz technique based on the energy expressions for the strongly coupled plate-cavity system. The accuracy, convergence, and efficiency of the present method are verified by comparing with the results available in the FEM and literature. Finally, the effects of the structural boundary conditions, cavity depth, and structural length-width ratio on the coupling natural frequency and the steady-state responses under three excitation conditions are analyzed.


---

## Body

## 1. Introduction

The structure-acoustic interaction system consisting of a flexible vibrating plate and a closed acoustic cavity is widely used in daily life and practical engineering, for instance, bedrooms, aircraft, ship cabins, and sonar platform cavity. In recent years, with the growing requirements for acoustical performance of products as well as comfortable working and living conditions, a clear understanding of acoustic characteristics of structure-acoustic interaction system is the key to acoustical product design and noise control. However, most research is carried out with the assumption of weak coupling between structure and acoustics. With the development of engineering field, the strong coupling situation where, for example, the contained acoustic medium has large mass density, the plate is thin, and the cavity is shallow is becoming more and more common. Therefore, it is of great value to deeply study the acoustic characteristics of the strongly coupled structure-acoustic system.In practice, a flexible vibrating plate backed by a closed rectangular cavity is a classical representation of numerous engineering backgrounds and has been extensively studied for a long time. In 1963, Lyon [1] and Dowell and Voss [2] firstly studied the acoustic characteristic of a flexible plate coupled with a rectangular cavity and theoretically analyzed the attenuation of internal noise in the cavity. Dowell gave the theoretical solution of structural displacement and sound pressure in form of Fourier series. On this basis, Pretlove [3, 4] proposed to calculate the vibration response of the structure by using the structural modes in vacuum and analyzed the effect of the cavity depth on the vibration and acoustical response. In order to express the effect of acoustic cavity on structural vibration, Fahy [5] proposed the concept of modal coupling, which uses modal coupling coefficients to characterize the interaction between structure and cavity. By combining numerical and analytical methods, Qaisi [6] obtained the mass and stiffness matrices of the coupled system and analyzed the effect of cavity depth and boundary conditions on the natural frequencies of the coupled system. Pan and Bies [7] gave a comprehensive literature review of structure-acoustic interaction system and studied the coupling mechanism. Kim and Brennan [8] extended the classical modal coupling theory of structure-acoustic interaction system. A compact matrix form of a plate-cavity coupled system was established based on the impedance and mobility method, which can be easily solved in numerical computer software. Using the variation of modal coupling theory, Geng et al. [9] analyzed the response of irregular structure-acoustic interaction system. Luo et al. [10] derived the compact matrix formulation by using Green’s function theorem and proposed a numerical simulation method to calculate the response of sound pressure and velocity. Zhang et al. [11] studied vibration characteristics of the laminated plate-cavity coupling system using an improved Fourier method (IFM). Furthermore, they [12] analyzed vibro-acoustic property of a laminated rotary plate coupled with an impedance cavity under complex boundary conditions by combining the simplified plate theory (SPT) and the improved Fourier–Ritz method. Lee [13] proposed to examine the effect of cavity length on the natural frequency of a nonlinear panel coupled with extended cavity by using the multilevel residue harmonic balance method.In addition to the research on weakly coupled panel-cavity system, the study of strongly coupled panel-cavity system has attracted the interest of many scholars. Tournour and Atalla [14] analyzed the effect of the truncation value on the response of the strongly coupled plate-cavity system and pointed out that the vibration and sound response of the system could not be accurately predicted by increasing the truncation order in the strongly coupled case and proposed a pseudostatic correction to improve its convergence. David and Menelle [15] proposed a decoupled modal projection method for analyzing the characteristic frequency of the strongly coupled plate-cavity system. In order to investigate the acoustic characteristic of an elastically restrained plate coupled with an enclosure, Du et al. [16–19] proposed an improved Fourier series by introducing auxiliary function to ensure the continuity of vibration velocity on the plate-cavity interface. Shi et al. [20] used coordinate transformation to map irregular cavities to the unit square domain and studied its acoustic properties based on the three-dimensional improved Fourier series. On the basis, Chen et al. [21, 22] modeled the coupling system by replacing the improved Fourier series with Chebyshev polynomial series (CPSM), which was more concise and required less computation. Kim et al. [23] presented a new formulation of the coupled reduced-order modeling technique for fluid-structure interaction problems.In the view of the existing literatures, different from the research on the weakly coupled panel-cavity system, the research on the strongly coupled plate-cavity system focuses on the system characteristics and less on the steady-state response. Meanwhile, the computational cost of IFM and CPSM is high [20, 21]. Inspired by these limitations, based on the principle of energy, the acoustical characteristics of the strongly coupled plate-cavity system with arbitrary structure boundary constraints is analyzed by using Legendre polynomial series in this paper. A general expression of structural boundary conditions is established. The structural displacement functions of the plate and sound pressure functions of the cavity are expanded by Legendre polynomial series, respectively. Based on the energy principle, the Lagrange functions of structure and acoustics in the strongly coupled plate-cavity system are developed, respectively. The acoustic characteristics of the system are obtained by the Rayleigh–Ritz technique. Compared with FEM and literatures, the accuracy, convergence, and efficiency of the present method are verified. Finally, the effects of the structural boundary conditions, cavity depth, and structural length-width ratio on the natural frequency are analyzed. On this basis, the steady-state response of the system under force source excitation, internal acoustic source excitation, and a combination of both is studied.


## 2. Theoretical Formulations

### 2.1. Model Description

As shown in Figure1, a plate-cavity coupled system consisting of a plate with arbitrary boundary conditions and an acoustic cavity is established to analyze the acoustic characteristics. The rectangular cavity has a length, width, and height of Lx, Ly, and Lz, respectively, and is filled with a heavy medium (e.g., water) with a mass density of ρ0 and a speed of sound propagation of c0. The plate has the same length and width as the cavity and is coupled to it at the top surface of the cavity. In addition, the thickness of the plate is denoted as h, and the mass density is ρh. A point force source excitation Fp on the plate xp,yp and a point sound source excitation Qc at the cavity xc,yc,zc are introduced to analyze the steady-state response of the coupling system.Figure 1
Schematic of the considered flexible plate strongly coupled with a rectangular cavity.The strongly coupled plate-cavity system can be divided into the acoustic cavity system and flexible plate system. For the acoustic cavity system, the internal sound pressure satisfies the Helmholtz equation and boundary conditions:(1)∇2px,y,z,t+ω2c02px,y,z,t=0,∂px,y,z,t∂n=ρ0ω2wx,y,t,on the coupling surface,0,on the rigid wall,where ∇ is the Laplace operator, ω=2πf is the excitation frequency, wx,y,t is the structural displacement of the flexible plate, and n is the external normal direction of the flexible plate.For the flexible plate system, the effect of sound pressure at the interface should be considered as well as the point force source excitation. Based on the theory of structural dynamics, the differential equation for the elastic structural plate takes the following form:(2)D∇4wx,y,t+m∂2wx,y,t∂t2=px,y,Lz,t−Fp,where D and m are the structural bending stiffness and density of plate, respectively. Each boundary constraint of the plate can be described by the linear spring stiffness coefficient k and rotational spring stiffness coefficient K; take x=0 as an example [16]:(3)kx=0wx,y,t=−D∂3wx,y,t∂x3+2−η∂3wx,y,t∂x∂y2,Kx=0∂wx,y,t∂x=D∂2wx,y,t∂x2+η∂2wx,y,t∂y2,where η is Poisson’s ratio and kx=0 and Kx=0 are the linear spring stiffness coefficient and rotational spring stiffness coefficient, respectively. The variation of the elastic boundary constraints can be realized by setting the stiffness coefficient to zero or infinity [16]. For example, if the linear spring stiffness coefficient is set to infinity (which is represented by 1011 in following numerical calculations) and the rotational spring stiffness coefficient is set to zero, the simply supported boundary condition can be obtained. If the linear spring stiffness coefficient and the rotational spring stiffness coefficient are both set to infinity, the clamped boundary conditions can be obtained.


### 2.2. Series Representation of Sound Pressure and Structural Displacement Functions

The velocity continuity will not be satisfied at the coupling interface with the strong coupling between plate and cavity if the rigidly walled cavity mode is used to represent the sound pressure as well as the vacuum structure mode is used to represent the displacement. In order to overcome the problem, the structural displacement and sound pressure are expanded by Legendre polynomial series; in addition, coordinate transformation is required since the Legendre polynomials define the interval as [−1, 1]:(4)α=2xLx−1,β=2yLy−1,γ=2zLz−1.The sound pressure and structural displacement functions can then be expressed as(5)pα,β,γ,t=e−jwt∑mα=0∞∑mβ=0∞∑mγ=0∞AmαmβmγPmααPmββPmγγ,(6)wα,β,t=e−jwt∑nα=0∞∑nβ=0∞BnαnβPnααPnββ,where Amαmβmγ is the mα,mβ,mγth term of the expansion coefficient vectors of sound pressure A and Bnαnβ is the nα,nβth term of the expansion coefficient vectors of structural displacement B. The time item e−jwt will be omitted below for simplicity; Pmαα is Legendre polynomial [24] with the form as(7)Pmαα=12mαmα!dmαdαmαα2−1mα.It should be noted that although the infinite polynomial expansion terms are included in equations (5) and (6), only the first finite polynomials will be used in numerical calculations. In addition, Legendre polynomials are not only complete and recursive as other polynomials, but also their L2 inner product satisfies orthogonality:(8)∫−11PmααPmα′α=22mα+1δmαmα′,where δmαmα′ is the Kronecker function, only when mα=mα′, δmαmα′ is 1, otherwise, δmαmα′ is 0.


### 2.3. Energy Expression of Acoustic Cavity System and Flexible Plate System

Since the strongly coupled plate-cavity system can be divided into the acoustic cavity system and the flexible plate system, based on the energy principle, the Lagrange functions [18, 19] of the acoustic cavity system and the flexible plate system can be obtained, respectively. For the acoustic cavity system:(9)LC=UC−TC−WQ−WP2C,where UC is the total acoustic potential energy of the acoustic cavity, TC is the total kinetic energy of the acoustic cavity, WQ is the work done by the point sound source, and WP2C represents the work done by structural vibration to the medium in the cavity [16]. Their specific expressions are(10)UC=LxLyLz16ρc2∫−11∫−11∫−11p2α,β,γdαdβdγ,(11)TC=LxLyLz4ρω2∫−11∫−11∫−11∂pα,β,γLx∂α2+∂pα,β,γLy∂β2+∂pα,β,γLz∂γ2dαdβdγ,(12)WQ=−1jωQcpαc,βc,γc,(13)WP2C=LxLy4∫−11∫−11wα,βpα,β,1dαdβ.For the flexible plate system, the Lagrange function expression [18, 19] is shown as below:(14)LP=UP−TP−WP−WC2P,where UP is the total acoustic potential energy of the flexible plate, TP is the total kinetic energy of the flexible plate, WP is the work done by the point force source, and WC2P represents the work done by sound pressure to the flexible plate [20]. Their specific expressions are(15)UP=D2∫−11∫−114LyLx3∂2wα,β∂α22+4LxLy3∂2wα,β∂β22+8ηLxLy∂2wα,β∂α2∂2wα,β∂β2+81−ηLxLy∂2wα,β∂α2∂β22dαdβ+Lx4∫−11ky0w2α,β+4Ly2Ky0∂wα,β∂β2β=−1dα+Lx4∫−11kyLyw2α,β+4Ly2KyLy∂wα,β∂β2β=1dα+Ly4∫−11kx0w2α,β+4Lx2Kx0∂wα,β∂α2α=−1dβ+Ly4∫−11kxLxw2α,β+4Lx2KxLx∂wα,β∂α2α=1dβ,(16)TP=LxLy8ρhω2∫−11∫−11wα,βdαdβ,(17)WP=wαp,βpFP,(18)WC2P=−WP2C.


### 2.4. Solution Method

Substitute equation (5) and equations (10)∼(13) into equation (9). Substitute equation (6) and equations (15)∼(18) into equation (14). According to the Rayleigh–Ritz method [25, 26], set the partial derivative of the Lagrange functions LC with respect to expansion coefficient vectors of sound pressure A and the Lagrange functions LP with respect to expansion coefficient vectors of structural displacement B as zero, respectively,(19)∂LC∂A=0,∂LP∂B=0.The linear algebraic equations of the expansion coefficient vectorsA and B can be obtained:(20)KP−ω2MPB+KC2PA=F,(21)KC−ω2MCA+ω2KP2CB=Q.The overall form can be developed by combining equations (20) and (21):(22)KPKC2P0KC−ω2MP0−KP2CMCBA=FQ,where KP and KC are the stiffness matrices of the acoustic cavity system and the flexible plate system, MP and MC are the mass matrices of the acoustic cavity system and the flexible plate system, KC2P represents the coupling effect matrix of the acoustic cavity system to the flexible plate system, and KP2C represents the coupling effect matrix of the flexible plate system to the acoustic cavity system.According to the principle of action and reaction, we can get the equality equation ofKP2C=KC2PT. F is is the force vector which acted on the flexible plate and Q is the sound source vector of the internal point sound source. The acoustic characteristics of the strongly coupled plate-cavity system can be obtained by setting F and Q to zero. On the contrary, when considering the external force vectors F and Q, the steady-state responses of the coupling system can be developed. Since the L2 inner product of Legendre polynomial satisfies the orthogonality, the high-multiple integral term in each matrix element can be simplified to low-multiple integral or direct product term. According to equation (8), the matrix element can be simplified into the following formulae:(23)kpnαnβ,nα′nβ′=4DLxLy22nβ+1δnβnβ′Ξnαnα′2,2+22nα+1δnαnα′Ψnβnβ′2,2+ηΞnαnα′2,0Ψnβnβ′0,2+ηΞnαnα′0,2Ψnβnβ′2,0+21−ηΞnαnα′1,1Ψnβnβ′1,1+Lx2nα+1δnαnα′ky0Pnβ−1Pnβ′−1+4Ly2Ky0Pnβ′−1Pnβ′′−1+kyLyPnβ1Pnβ′1+4Ly2KyLyPnβ′′1Pnβ′′1+Ly2q+1δnβnβ′kx0Pnα−1Pnα′−1+4Lx2Kx0Pnα′−1Pnα′′−1+kxLxPnα1Pnα′1+4Lx2KxLxPnα′1Pnα′′1,(24)kcmαmβmγ,mα′mβ′mγ′=LxLyLzρ02Lx22mβ+12mγ+1δmβmβ′δmγmγ′Ξmαmα′1,1+2Ly22mα+12mγ+1δmαmα′δmγmγ′Ψmβmβ′1,1+2Lz22mα+12mβ+1δmαmα′δmβmβ′ϒmγmγ′1,1,(25)mpnαnβ,nα′nβ′=LxLyρhh2nα+12nβ+1δnαnα′δnβnβ′,(26)mcmαmβmγ,mα′mβ′mγ′=LxLyLzρ0c022mα+12mβ+12mγ+1δmαmα′δmβmβ′δmγmγ′,(27)kp2cmαmβmγ,nαnβ=LxLy2nα+12nβ+1δnαmαδnβmβPmγ1,(28)fnαnβ=FPPnααpPnββp,(29)qmαmβmγ=−jωQcPmααcPmββcPmγγc,where(30)Ξabv,v¯=∫−11dvPavαLxvdαvdv¯Pav¯αLxv¯dαv¯a=nα,mα;b=nα′,mα′,(31)Ψabv,v¯=∫−11dvPavβLyvdβvdv¯Pav¯βLyv¯dβv¯a=nβ,mβ;b=nβ′,mβ′,(32)ϒa,bv,v¯=∫−11dvPavγLzvdγvdv¯Pav¯γLzv¯dγv¯a=mγ;b=mγ′.


## 2.1. Model Description

As shown in Figure1, a plate-cavity coupled system consisting of a plate with arbitrary boundary conditions and an acoustic cavity is established to analyze the acoustic characteristics. The rectangular cavity has a length, width, and height of Lx, Ly, and Lz, respectively, and is filled with a heavy medium (e.g., water) with a mass density of ρ0 and a speed of sound propagation of c0. The plate has the same length and width as the cavity and is coupled to it at the top surface of the cavity. In addition, the thickness of the plate is denoted as h, and the mass density is ρh. A point force source excitation Fp on the plate xp,yp and a point sound source excitation Qc at the cavity xc,yc,zc are introduced to analyze the steady-state response of the coupling system.Figure 1
Schematic of the considered flexible plate strongly coupled with a rectangular cavity.The strongly coupled plate-cavity system can be divided into the acoustic cavity system and flexible plate system. For the acoustic cavity system, the internal sound pressure satisfies the Helmholtz equation and boundary conditions:(1)∇2px,y,z,t+ω2c02px,y,z,t=0,∂px,y,z,t∂n=ρ0ω2wx,y,t,on the coupling surface,0,on the rigid wall,where ∇ is the Laplace operator, ω=2πf is the excitation frequency, wx,y,t is the structural displacement of the flexible plate, and n is the external normal direction of the flexible plate.For the flexible plate system, the effect of sound pressure at the interface should be considered as well as the point force source excitation. Based on the theory of structural dynamics, the differential equation for the elastic structural plate takes the following form:(2)D∇4wx,y,t+m∂2wx,y,t∂t2=px,y,Lz,t−Fp,where D and m are the structural bending stiffness and density of plate, respectively. Each boundary constraint of the plate can be described by the linear spring stiffness coefficient k and rotational spring stiffness coefficient K; take x=0 as an example [16]:(3)kx=0wx,y,t=−D∂3wx,y,t∂x3+2−η∂3wx,y,t∂x∂y2,Kx=0∂wx,y,t∂x=D∂2wx,y,t∂x2+η∂2wx,y,t∂y2,where η is Poisson’s ratio and kx=0 and Kx=0 are the linear spring stiffness coefficient and rotational spring stiffness coefficient, respectively. The variation of the elastic boundary constraints can be realized by setting the stiffness coefficient to zero or infinity [16]. For example, if the linear spring stiffness coefficient is set to infinity (which is represented by 1011 in following numerical calculations) and the rotational spring stiffness coefficient is set to zero, the simply supported boundary condition can be obtained. If the linear spring stiffness coefficient and the rotational spring stiffness coefficient are both set to infinity, the clamped boundary conditions can be obtained.


## 2.2. Series Representation of Sound Pressure and Structural Displacement Functions

The velocity continuity will not be satisfied at the coupling interface with the strong coupling between plate and cavity if the rigidly walled cavity mode is used to represent the sound pressure as well as the vacuum structure mode is used to represent the displacement. In order to overcome the problem, the structural displacement and sound pressure are expanded by Legendre polynomial series; in addition, coordinate transformation is required since the Legendre polynomials define the interval as [−1, 1]:(4)α=2xLx−1,β=2yLy−1,γ=2zLz−1.The sound pressure and structural displacement functions can then be expressed as(5)pα,β,γ,t=e−jwt∑mα=0∞∑mβ=0∞∑mγ=0∞AmαmβmγPmααPmββPmγγ,(6)wα,β,t=e−jwt∑nα=0∞∑nβ=0∞BnαnβPnααPnββ,where Amαmβmγ is the mα,mβ,mγth term of the expansion coefficient vectors of sound pressure A and Bnαnβ is the nα,nβth term of the expansion coefficient vectors of structural displacement B. The time item e−jwt will be omitted below for simplicity; Pmαα is Legendre polynomial [24] with the form as(7)Pmαα=12mαmα!dmαdαmαα2−1mα.It should be noted that although the infinite polynomial expansion terms are included in equations (5) and (6), only the first finite polynomials will be used in numerical calculations. In addition, Legendre polynomials are not only complete and recursive as other polynomials, but also their L2 inner product satisfies orthogonality:(8)∫−11PmααPmα′α=22mα+1δmαmα′,where δmαmα′ is the Kronecker function, only when mα=mα′, δmαmα′ is 1, otherwise, δmαmα′ is 0.


## 2.3. Energy Expression of Acoustic Cavity System and Flexible Plate System

Since the strongly coupled plate-cavity system can be divided into the acoustic cavity system and the flexible plate system, based on the energy principle, the Lagrange functions [18, 19] of the acoustic cavity system and the flexible plate system can be obtained, respectively. For the acoustic cavity system:(9)LC=UC−TC−WQ−WP2C,where UC is the total acoustic potential energy of the acoustic cavity, TC is the total kinetic energy of the acoustic cavity, WQ is the work done by the point sound source, and WP2C represents the work done by structural vibration to the medium in the cavity [16]. Their specific expressions are(10)UC=LxLyLz16ρc2∫−11∫−11∫−11p2α,β,γdαdβdγ,(11)TC=LxLyLz4ρω2∫−11∫−11∫−11∂pα,β,γLx∂α2+∂pα,β,γLy∂β2+∂pα,β,γLz∂γ2dαdβdγ,(12)WQ=−1jωQcpαc,βc,γc,(13)WP2C=LxLy4∫−11∫−11wα,βpα,β,1dαdβ.For the flexible plate system, the Lagrange function expression [18, 19] is shown as below:(14)LP=UP−TP−WP−WC2P,where UP is the total acoustic potential energy of the flexible plate, TP is the total kinetic energy of the flexible plate, WP is the work done by the point force source, and WC2P represents the work done by sound pressure to the flexible plate [20]. Their specific expressions are(15)UP=D2∫−11∫−114LyLx3∂2wα,β∂α22+4LxLy3∂2wα,β∂β22+8ηLxLy∂2wα,β∂α2∂2wα,β∂β2+81−ηLxLy∂2wα,β∂α2∂β22dαdβ+Lx4∫−11ky0w2α,β+4Ly2Ky0∂wα,β∂β2β=−1dα+Lx4∫−11kyLyw2α,β+4Ly2KyLy∂wα,β∂β2β=1dα+Ly4∫−11kx0w2α,β+4Lx2Kx0∂wα,β∂α2α=−1dβ+Ly4∫−11kxLxw2α,β+4Lx2KxLx∂wα,β∂α2α=1dβ,(16)TP=LxLy8ρhω2∫−11∫−11wα,βdαdβ,(17)WP=wαp,βpFP,(18)WC2P=−WP2C.


## 2.4. Solution Method

Substitute equation (5) and equations (10)∼(13) into equation (9). Substitute equation (6) and equations (15)∼(18) into equation (14). According to the Rayleigh–Ritz method [25, 26], set the partial derivative of the Lagrange functions LC with respect to expansion coefficient vectors of sound pressure A and the Lagrange functions LP with respect to expansion coefficient vectors of structural displacement B as zero, respectively,(19)∂LC∂A=0,∂LP∂B=0.The linear algebraic equations of the expansion coefficient vectorsA and B can be obtained:(20)KP−ω2MPB+KC2PA=F,(21)KC−ω2MCA+ω2KP2CB=Q.The overall form can be developed by combining equations (20) and (21):(22)KPKC2P0KC−ω2MP0−KP2CMCBA=FQ,where KP and KC are the stiffness matrices of the acoustic cavity system and the flexible plate system, MP and MC are the mass matrices of the acoustic cavity system and the flexible plate system, KC2P represents the coupling effect matrix of the acoustic cavity system to the flexible plate system, and KP2C represents the coupling effect matrix of the flexible plate system to the acoustic cavity system.According to the principle of action and reaction, we can get the equality equation ofKP2C=KC2PT. F is is the force vector which acted on the flexible plate and Q is the sound source vector of the internal point sound source. The acoustic characteristics of the strongly coupled plate-cavity system can be obtained by setting F and Q to zero. On the contrary, when considering the external force vectors F and Q, the steady-state responses of the coupling system can be developed. Since the L2 inner product of Legendre polynomial satisfies the orthogonality, the high-multiple integral term in each matrix element can be simplified to low-multiple integral or direct product term. According to equation (8), the matrix element can be simplified into the following formulae:(23)kpnαnβ,nα′nβ′=4DLxLy22nβ+1δnβnβ′Ξnαnα′2,2+22nα+1δnαnα′Ψnβnβ′2,2+ηΞnαnα′2,0Ψnβnβ′0,2+ηΞnαnα′0,2Ψnβnβ′2,0+21−ηΞnαnα′1,1Ψnβnβ′1,1+Lx2nα+1δnαnα′ky0Pnβ−1Pnβ′−1+4Ly2Ky0Pnβ′−1Pnβ′′−1+kyLyPnβ1Pnβ′1+4Ly2KyLyPnβ′′1Pnβ′′1+Ly2q+1δnβnβ′kx0Pnα−1Pnα′−1+4Lx2Kx0Pnα′−1Pnα′′−1+kxLxPnα1Pnα′1+4Lx2KxLxPnα′1Pnα′′1,(24)kcmαmβmγ,mα′mβ′mγ′=LxLyLzρ02Lx22mβ+12mγ+1δmβmβ′δmγmγ′Ξmαmα′1,1+2Ly22mα+12mγ+1δmαmα′δmγmγ′Ψmβmβ′1,1+2Lz22mα+12mβ+1δmαmα′δmβmβ′ϒmγmγ′1,1,(25)mpnαnβ,nα′nβ′=LxLyρhh2nα+12nβ+1δnαnα′δnβnβ′,(26)mcmαmβmγ,mα′mβ′mγ′=LxLyLzρ0c022mα+12mβ+12mγ+1δmαmα′δmβmβ′δmγmγ′,(27)kp2cmαmβmγ,nαnβ=LxLy2nα+12nβ+1δnαmαδnβmβPmγ1,(28)fnαnβ=FPPnααpPnββp,(29)qmαmβmγ=−jωQcPmααcPmββcPmγγc,where(30)Ξabv,v¯=∫−11dvPavαLxvdαvdv¯Pav¯αLxv¯dαv¯a=nα,mα;b=nα′,mα′,(31)Ψabv,v¯=∫−11dvPavβLyvdβvdv¯Pav¯βLyv¯dβv¯a=nβ,mβ;b=nβ′,mβ′,(32)ϒa,bv,v¯=∫−11dvPavγLzvdγvdv¯Pav¯γLzv¯dγv¯a=mγ;b=mγ′.


## 3. Method Validation

The accuracy and convergence of the method are verified by comparing with the results available in the FEM and literature. The FEM model is established in COMSOL, and the sound-shell coupling module is used to model the strongly coupled plate-cavity system. To obtain reliable results, the maximum mesh size is set to 1/6 of the minimum wavelength of the analysis band. A water-filled cavity bounded by a simply supported plate is considered as an analysis model which has been studied in Ref [21]. The specific parameters of the plate-cavity system are shown in Table 1.Table 1
The model parameters.
CavityLength (m)0.29Width (m)0.35Depth (m)0.14Water density (kg⋅m−3)1000Sound velocity (m⋅s−1)1500PlateLength (m)0.29Width (m)0.35Thickness (m)0.0015Poisson’s ratio0.3Young’s modulus (Pa)7.2×1010Mass density (kg⋅m−3)2700


### 3.1. Accuracy and Convergence Analysis

Table2 shows the first eight natural frequencies of the strongly coupled plate-cavity system. Different truncation values are selected for sound pressure expression (represented by mα, mβ, mγ) and structural displacement expression (represented by nα, nβ) to analyze the accuracy and convergence of the present method. It can be seen from Table 2 that with the increase of the truncation values of the sound pressure expression and structural displacement expression, the results of present method converge rapidly and are in good agreement with the results available in the FEM and literature.Table 2
The first eight natural frequencies of the strongly coupled plate-cavity system.
mα×mβ×mγnα×nβMode number123456788×8×86×633.56446.52378.589112.209163.658171.812207.975269.7278×833.54846.58078.57585.232125.997137.665162.189175.86310×1033.54746.57978.57284.802125.100137.282161.111172.60612×1233.54746.57978.57284.797125.090137.275161.095172.55214×1433.54746.57978.57284.797125.090137.275161.094172.55210×10×106×633.56446.62278.587112.179163.590171.759207.902269.5888×833.54846.58078.57485.211125.963137.635162.141175.74210×1033.54746.57878.57184.764125.033137.231160.992172.26812×1233.54746.57878.57184.757125.020137.220160.965172.17514×1433.54746.57878.57184.757125.019137.220160.965172.17412×12×126×633.56446.62278.587112.176163.582171.754207.894269.5718×833.54746.58078.57385.211125.962137.634162.140175.74110×1033.54746.57878.57184.764125.032137.229160.991172.26512×1233.54746.57878.57184.757125.018137.218160.962172.16714×1433.54746.57878.57084.757125.018137.218160.962172.166FEM33.50946.51378.41384.577124.670136.810160.420171.570Ref [21]33.57446.59678.60884.102124.680137.112161.082172.334Figure2 shows the first five modes of displacement distribution of the plate calculated by FEM and the present method with the truncation values values mα=mβ=mγ=10, nα=nβ=10. In addition to accurately expressing the displacement distribution of the structure, we can also find that the strong coupling between plate and cavity caused by the incompressibility of the heavy fluid medium changes the asymmetric mode displacement distribution of the plate. Different from the structural displacement distribution in vacuum [27], where the first piston-like mode in vacuum is replaced by the symmetric mode, the fourth and fifth modes (the fifth and sixth modes in vacuum) have obvious changes, and the second and third modes (the third and fourth modes in vacuum) basically have no effect. The incompressibility of the heavy fluid medium leads to significant changes in the asymmetrically distributed modes, while the symmetrically distributed modes counteract the effects of coupling.Figure 2
The first five structural modes of the strongly coupled panel-cavity system. (a) FEM. (b) Present method.
(a)(b)


### 3.2. Computational Efficiency Analysis

In contrast to Ref [21], which proposed to expand the structural displacement function and the sound pressure function with Chebyshev polynomial series, this paper uses Legendre polynomial series to describe the sound pressure and displacement. Since the L2 inner product of Legendre polynomial satisfies the orthogonality, it will simplify the calculation. For code implementations of the two methods on computers with the same computing power, the main influence on the computational efficiency is the dimensionality of one-dimension integration terms. The appendix gives expressions for each element of the stiffness matrix, mass matrix, and coupling interaction matrix from Ref [21]. For the same analytical model with the same truncation values, the total dimensionality of the one-dimension integral terms for present method and Ref [21] can be obtained according to equations (23)–(32) and the Appendix:(33)Dimpre=nα2+nβ2+3nα2nβ2+mα2+mβ2+mγ2,(34)Dimref=6nα2nβ2+4mα2mβ2mγ2+nα2+nβ2+nαmαnβmβ.It can be seen from equations (33) and (34) that the total dimensionality of the one-dimension integral term is greatly reduced, making the code easier to handle and improving computational efficiency when using Legendre polynomial series to describe sound pressure and structural displacement.


## 3.1. Accuracy and Convergence Analysis

Table2 shows the first eight natural frequencies of the strongly coupled plate-cavity system. Different truncation values are selected for sound pressure expression (represented by mα, mβ, mγ) and structural displacement expression (represented by nα, nβ) to analyze the accuracy and convergence of the present method. It can be seen from Table 2 that with the increase of the truncation values of the sound pressure expression and structural displacement expression, the results of present method converge rapidly and are in good agreement with the results available in the FEM and literature.Table 2
The first eight natural frequencies of the strongly coupled plate-cavity system.
mα×mβ×mγnα×nβMode number123456788×8×86×633.56446.52378.589112.209163.658171.812207.975269.7278×833.54846.58078.57585.232125.997137.665162.189175.86310×1033.54746.57978.57284.802125.100137.282161.111172.60612×1233.54746.57978.57284.797125.090137.275161.095172.55214×1433.54746.57978.57284.797125.090137.275161.094172.55210×10×106×633.56446.62278.587112.179163.590171.759207.902269.5888×833.54846.58078.57485.211125.963137.635162.141175.74210×1033.54746.57878.57184.764125.033137.231160.992172.26812×1233.54746.57878.57184.757125.020137.220160.965172.17514×1433.54746.57878.57184.757125.019137.220160.965172.17412×12×126×633.56446.62278.587112.176163.582171.754207.894269.5718×833.54746.58078.57385.211125.962137.634162.140175.74110×1033.54746.57878.57184.764125.032137.229160.991172.26512×1233.54746.57878.57184.757125.018137.218160.962172.16714×1433.54746.57878.57084.757125.018137.218160.962172.166FEM33.50946.51378.41384.577124.670136.810160.420171.570Ref [21]33.57446.59678.60884.102124.680137.112161.082172.334Figure2 shows the first five modes of displacement distribution of the plate calculated by FEM and the present method with the truncation values values mα=mβ=mγ=10, nα=nβ=10. In addition to accurately expressing the displacement distribution of the structure, we can also find that the strong coupling between plate and cavity caused by the incompressibility of the heavy fluid medium changes the asymmetric mode displacement distribution of the plate. Different from the structural displacement distribution in vacuum [27], where the first piston-like mode in vacuum is replaced by the symmetric mode, the fourth and fifth modes (the fifth and sixth modes in vacuum) have obvious changes, and the second and third modes (the third and fourth modes in vacuum) basically have no effect. The incompressibility of the heavy fluid medium leads to significant changes in the asymmetrically distributed modes, while the symmetrically distributed modes counteract the effects of coupling.Figure 2
The first five structural modes of the strongly coupled panel-cavity system. (a) FEM. (b) Present method.
(a)(b)


## 3.2. Computational Efficiency Analysis

In contrast to Ref [21], which proposed to expand the structural displacement function and the sound pressure function with Chebyshev polynomial series, this paper uses Legendre polynomial series to describe the sound pressure and displacement. Since the L2 inner product of Legendre polynomial satisfies the orthogonality, it will simplify the calculation. For code implementations of the two methods on computers with the same computing power, the main influence on the computational efficiency is the dimensionality of one-dimension integration terms. The appendix gives expressions for each element of the stiffness matrix, mass matrix, and coupling interaction matrix from Ref [21]. For the same analytical model with the same truncation values, the total dimensionality of the one-dimension integral terms for present method and Ref [21] can be obtained according to equations (23)–(32) and the Appendix:(33)Dimpre=nα2+nβ2+3nα2nβ2+mα2+mβ2+mγ2,(34)Dimref=6nα2nβ2+4mα2mβ2mγ2+nα2+nβ2+nαmαnβmβ.It can be seen from equations (33) and (34) that the total dimensionality of the one-dimension integral term is greatly reduced, making the code easier to handle and improving computational efficiency when using Legendre polynomial series to describe sound pressure and structural displacement.


## 4. Numerical Results and Discussion

### 4.1. Natural Characteristics Analysis

In this section, the effects of the structural boundary conditions, cavity depth, and structural length-width ratio on the natural frequency of simply supported strongly coupled plate-cavity system are studied. The model parameters are shown in Table1, which are used in subsequent calculation. Considering the correctness of the results and the efficiency of the calculation, the truncation values can be set to mα=mβ=mγ=10 and nα=nβ=10 in the following numerical calculations.


#### 4.1.1. Effects of Structural Boundary Conditions

In this section, the natural frequencies under different boundary conditions are analyzed by changing the spring stiffness coefficient. “C” represents clamped boundary condition and “S” represents simply supported boundary condition. For example, “CSSC” means that x=0 is clamped supported, x=Lx is simply supported, y=0 is simply supported, and y=Ly is clamped supported. Table 3 shows the natural frequency under different boundary conditions, and the finite element results are selected as the reference results. For different boundary conditions, the results of the present method and FEM are in good agreement, which on the one hand further demonstrates the accuracy of present method in predicting the acoustic characteristics of strongly coupled systems. On the other hand, the present method is more flexible than traditional methods because it only needs to change the spring stiffness coefficients to achieve different combinations of boundary conditions. In addition, we can also find that the natural frequency of the coupled system under the clamped supported boundary is higher than that of the simply supported boundary. With the increase of the number of the clamped supported boundary, the natural frequency gradually increases.Table 3
Natural frequencies under different boundary conditions.
Boundary conditionsResultsMode number12345678SSSSPresent method33.54746.57878.57184.764125.033137.231160.992172.268FEM33.50946.51378.41384.577124.670136.810160.420171.570SSSCPresent method40.38948.41085.34096.514127.324150.292168.131191.673FEM40.30948.33085.10996.178126.930149.650167.440190.370SSCCPresent method48.62550.75593.262108.600130.969164.836175.427213.522FEM48.48050.65192.932108.070130.500163.880174.590211.580SCCCPresent method51.43962.058104.229111.989150.289175.091194.538217.899FEM51.27061.869103.770111.400149.510173.960193.280215.930CCCCPresent method55.19875.789115.662117.261172.088187.087210.256228.636FEM54.99375.459115.010116.620170.820185.720208.430226.560


#### 4.1.2. Effects of the Cavity Depth

The natural frequencies in various cavity depths in the range of 0.01–1.0 m are studied. It should be noted that since the cavity depth is variable, the truncation value in the direction of cavity depth is selected as follows: when the cavity depth is less than 0.5 m,mγ=10; otherwise, mγ=15. Figure 3 shows the first eight natural frequency curves in different cavity depths. It can be seen that the frequency variation trend in different order is consistent. When the cavity depth is shallow, the frequency gradually increases with the increase of the cavity depth. As the cavity depth deepens further, the frequency does not change. In addition, it can be found that the frequency variation of the higher order is larger than that of the lower order under the same cavity depth variation.Figure 3
Variation of the first eight-order coupling frequency with various cavity depths.


#### 4.1.3. Effects of Structural Length-Width Ratio

Keeping cavity depth and the length of plate constant, Figure4 shows the variation of the first eight natural frequencies of the plate-cavity system under different length-width ratios Ly/Lx. As the width of the plate changes, the truncation value in this direction should also change as follows: when the length is less than 0.5 m, mβ=nβ=10, and when the length is less than 1 m but is greater than 0.5 m, mβ=nβ=15; otherwise, mβ=nβ=20. It can be seen from Figure 4 that the different order natural frequency variation rules of the coupled system under different length-width ratios are relatively consistent. With the increase of length-width ratio, the coupling frequency decreases gradually until it remains constant. In addition, when the length-width ratio is small, the coupling frequency changes more with the length-width ratio. We can slightly adjust the natural frequency of the coupling system by changing the length-width ratio.Figure 4
Variation of the first eight-order coupling frequency with various length-width ratios.


### 4.2. Steady-State Response Analysis

On the basis of the above results, the steady-state responses of the plate-cavity system are studied in this section under three excitation conditions, which are point force source excitation, internal point acoustic source excitation, and a combination of both. The model parameters are shown in Table1. The point force is applied at the thin rectangular plate at (0.04, 0.08, 0.14) with the magnitude F = 1 N, and the internal point acoustic source is applied at (0.17, 0.20, 0.10) with the volume velocity Qc=2×10−5m3⋅s−1. Two observation points are set at (0.10, 0.12, 0.14) and (0.15, 0.30, 0.14) of the plate to obtain the vibrating velocity response. Similarly, two observation points are set in the cavity at (0.08, 0.11, 0.07) and (0.20, 0.33, 0.06) to obtain the sound pressure response. In response calculation, structural damping ξ and acoustic medium damping ζ are both 0.01, which can be achieved by setting complex Young’s modulus as 7.2×1010×1−ξ and complex sound velocity as 1500/1−ζ, respectively. Figures 5–7 show the vibrating velocity response and the sound pressure response under three excitation conditions in the frequency range of 0–350 Hz. It can be seen that no matter the point force excitation, internal acoustic excitation, or a combination of both, the response calculation results of the proposed method are in good agreement with the FEM, which also demonstrates the accuracy of the proposed method in predicting the response. Meanwhile, it is feasible to introduce structural damping and acoustic damping by means of complex sound velocity and complex Young’s modulus.Figure 5
Response under point force source excitation. (a) Velocity response at (0.10, 0.12, 0.14). (b) Velocity response at (0.15, 0.30, 0.14). (c) Sound pressure response at (0.08, 0.11, 0.07). (d) Sound pressure response at (0.20, 0.33, 0.06).
(a)(b)(c)(d)Figure 6
Response under internal point acoustic source excitation. (a) Velocity response at (0.10, 0.12, 0.14). (b) Velocity response at (0.15, 0.30, 0.14). (c) Sound pressure response at (0.08, 0.11, 0.07). (d) Sound pressure response at (0.20, 0.33, 0.06).
(a)(b)(c)(d)Figure 7
Response under the combined excitation of point force source and internal point acoustic source. (a) Velocity response at (0.10, 0.12, 0.14). (b) Velocity response at (0.15, 0.30, 0.14). (c) Sound pressure response at (0.08, 0.11, 0.07). (d) Sound pressure response at (0.20, 0.33, 0.06).
(a)(b)(c)(d)


## 4.1. Natural Characteristics Analysis

In this section, the effects of the structural boundary conditions, cavity depth, and structural length-width ratio on the natural frequency of simply supported strongly coupled plate-cavity system are studied. The model parameters are shown in Table1, which are used in subsequent calculation. Considering the correctness of the results and the efficiency of the calculation, the truncation values can be set to mα=mβ=mγ=10 and nα=nβ=10 in the following numerical calculations.


### 4.1.1. Effects of Structural Boundary Conditions

In this section, the natural frequencies under different boundary conditions are analyzed by changing the spring stiffness coefficient. “C” represents clamped boundary condition and “S” represents simply supported boundary condition. For example, “CSSC” means that x=0 is clamped supported, x=Lx is simply supported, y=0 is simply supported, and y=Ly is clamped supported. Table 3 shows the natural frequency under different boundary conditions, and the finite element results are selected as the reference results. For different boundary conditions, the results of the present method and FEM are in good agreement, which on the one hand further demonstrates the accuracy of present method in predicting the acoustic characteristics of strongly coupled systems. On the other hand, the present method is more flexible than traditional methods because it only needs to change the spring stiffness coefficients to achieve different combinations of boundary conditions. In addition, we can also find that the natural frequency of the coupled system under the clamped supported boundary is higher than that of the simply supported boundary. With the increase of the number of the clamped supported boundary, the natural frequency gradually increases.Table 3
Natural frequencies under different boundary conditions.
Boundary conditionsResultsMode number12345678SSSSPresent method33.54746.57878.57184.764125.033137.231160.992172.268FEM33.50946.51378.41384.577124.670136.810160.420171.570SSSCPresent method40.38948.41085.34096.514127.324150.292168.131191.673FEM40.30948.33085.10996.178126.930149.650167.440190.370SSCCPresent method48.62550.75593.262108.600130.969164.836175.427213.522FEM48.48050.65192.932108.070130.500163.880174.590211.580SCCCPresent method51.43962.058104.229111.989150.289175.091194.538217.899FEM51.27061.869103.770111.400149.510173.960193.280215.930CCCCPresent method55.19875.789115.662117.261172.088187.087210.256228.636FEM54.99375.459115.010116.620170.820185.720208.430226.560


### 4.1.2. Effects of the Cavity Depth

The natural frequencies in various cavity depths in the range of 0.01–1.0 m are studied. It should be noted that since the cavity depth is variable, the truncation value in the direction of cavity depth is selected as follows: when the cavity depth is less than 0.5 m,mγ=10; otherwise, mγ=15. Figure 3 shows the first eight natural frequency curves in different cavity depths. It can be seen that the frequency variation trend in different order is consistent. When the cavity depth is shallow, the frequency gradually increases with the increase of the cavity depth. As the cavity depth deepens further, the frequency does not change. In addition, it can be found that the frequency variation of the higher order is larger than that of the lower order under the same cavity depth variation.Figure 3
Variation of the first eight-order coupling frequency with various cavity depths.


### 4.1.3. Effects of Structural Length-Width Ratio

Keeping cavity depth and the length of plate constant, Figure4 shows the variation of the first eight natural frequencies of the plate-cavity system under different length-width ratios Ly/Lx. As the width of the plate changes, the truncation value in this direction should also change as follows: when the length is less than 0.5 m, mβ=nβ=10, and when the length is less than 1 m but is greater than 0.5 m, mβ=nβ=15; otherwise, mβ=nβ=20. It can be seen from Figure 4 that the different order natural frequency variation rules of the coupled system under different length-width ratios are relatively consistent. With the increase of length-width ratio, the coupling frequency decreases gradually until it remains constant. In addition, when the length-width ratio is small, the coupling frequency changes more with the length-width ratio. We can slightly adjust the natural frequency of the coupling system by changing the length-width ratio.Figure 4
Variation of the first eight-order coupling frequency with various length-width ratios.


## 4.1.1. Effects of Structural Boundary Conditions

In this section, the natural frequencies under different boundary conditions are analyzed by changing the spring stiffness coefficient. “C” represents clamped boundary condition and “S” represents simply supported boundary condition. For example, “CSSC” means that x=0 is clamped supported, x=Lx is simply supported, y=0 is simply supported, and y=Ly is clamped supported. Table 3 shows the natural frequency under different boundary conditions, and the finite element results are selected as the reference results. For different boundary conditions, the results of the present method and FEM are in good agreement, which on the one hand further demonstrates the accuracy of present method in predicting the acoustic characteristics of strongly coupled systems. On the other hand, the present method is more flexible than traditional methods because it only needs to change the spring stiffness coefficients to achieve different combinations of boundary conditions. In addition, we can also find that the natural frequency of the coupled system under the clamped supported boundary is higher than that of the simply supported boundary. With the increase of the number of the clamped supported boundary, the natural frequency gradually increases.Table 3
Natural frequencies under different boundary conditions.
Boundary conditionsResultsMode number12345678SSSSPresent method33.54746.57878.57184.764125.033137.231160.992172.268FEM33.50946.51378.41384.577124.670136.810160.420171.570SSSCPresent method40.38948.41085.34096.514127.324150.292168.131191.673FEM40.30948.33085.10996.178126.930149.650167.440190.370SSCCPresent method48.62550.75593.262108.600130.969164.836175.427213.522FEM48.48050.65192.932108.070130.500163.880174.590211.580SCCCPresent method51.43962.058104.229111.989150.289175.091194.538217.899FEM51.27061.869103.770111.400149.510173.960193.280215.930CCCCPresent method55.19875.789115.662117.261172.088187.087210.256228.636FEM54.99375.459115.010116.620170.820185.720208.430226.560


## 4.1.2. Effects of the Cavity Depth

The natural frequencies in various cavity depths in the range of 0.01–1.0 m are studied. It should be noted that since the cavity depth is variable, the truncation value in the direction of cavity depth is selected as follows: when the cavity depth is less than 0.5 m,mγ=10; otherwise, mγ=15. Figure 3 shows the first eight natural frequency curves in different cavity depths. It can be seen that the frequency variation trend in different order is consistent. When the cavity depth is shallow, the frequency gradually increases with the increase of the cavity depth. As the cavity depth deepens further, the frequency does not change. In addition, it can be found that the frequency variation of the higher order is larger than that of the lower order under the same cavity depth variation.Figure 3
Variation of the first eight-order coupling frequency with various cavity depths.


## 4.1.3. Effects of Structural Length-Width Ratio

Keeping cavity depth and the length of plate constant, Figure4 shows the variation of the first eight natural frequencies of the plate-cavity system under different length-width ratios Ly/Lx. As the width of the plate changes, the truncation value in this direction should also change as follows: when the length is less than 0.5 m, mβ=nβ=10, and when the length is less than 1 m but is greater than 0.5 m, mβ=nβ=15; otherwise, mβ=nβ=20. It can be seen from Figure 4 that the different order natural frequency variation rules of the coupled system under different length-width ratios are relatively consistent. With the increase of length-width ratio, the coupling frequency decreases gradually until it remains constant. In addition, when the length-width ratio is small, the coupling frequency changes more with the length-width ratio. We can slightly adjust the natural frequency of the coupling system by changing the length-width ratio.Figure 4
Variation of the first eight-order coupling frequency with various length-width ratios.


## 4.2. Steady-State Response Analysis

On the basis of the above results, the steady-state responses of the plate-cavity system are studied in this section under three excitation conditions, which are point force source excitation, internal point acoustic source excitation, and a combination of both. The model parameters are shown in Table1. The point force is applied at the thin rectangular plate at (0.04, 0.08, 0.14) with the magnitude F = 1 N, and the internal point acoustic source is applied at (0.17, 0.20, 0.10) with the volume velocity Qc=2×10−5m3⋅s−1. Two observation points are set at (0.10, 0.12, 0.14) and (0.15, 0.30, 0.14) of the plate to obtain the vibrating velocity response. Similarly, two observation points are set in the cavity at (0.08, 0.11, 0.07) and (0.20, 0.33, 0.06) to obtain the sound pressure response. In response calculation, structural damping ξ and acoustic medium damping ζ are both 0.01, which can be achieved by setting complex Young’s modulus as 7.2×1010×1−ξ and complex sound velocity as 1500/1−ζ, respectively. Figures 5–7 show the vibrating velocity response and the sound pressure response under three excitation conditions in the frequency range of 0–350 Hz. It can be seen that no matter the point force excitation, internal acoustic excitation, or a combination of both, the response calculation results of the proposed method are in good agreement with the FEM, which also demonstrates the accuracy of the proposed method in predicting the response. Meanwhile, it is feasible to introduce structural damping and acoustic damping by means of complex sound velocity and complex Young’s modulus.Figure 5
Response under point force source excitation. (a) Velocity response at (0.10, 0.12, 0.14). (b) Velocity response at (0.15, 0.30, 0.14). (c) Sound pressure response at (0.08, 0.11, 0.07). (d) Sound pressure response at (0.20, 0.33, 0.06).
(a)(b)(c)(d)Figure 6
Response under internal point acoustic source excitation. (a) Velocity response at (0.10, 0.12, 0.14). (b) Velocity response at (0.15, 0.30, 0.14). (c) Sound pressure response at (0.08, 0.11, 0.07). (d) Sound pressure response at (0.20, 0.33, 0.06).
(a)(b)(c)(d)Figure 7
Response under the combined excitation of point force source and internal point acoustic source. (a) Velocity response at (0.10, 0.12, 0.14). (b) Velocity response at (0.15, 0.30, 0.14). (c) Sound pressure response at (0.08, 0.11, 0.07). (d) Sound pressure response at (0.20, 0.33, 0.06).
(a)(b)(c)(d)


## 5. Conclusions

In the paper, a strongly coupled structure-acoustic model for a rectangular cavity and its flexible wall is proposed, and the displacement of the plate and the sound pressure in the cavity are developed based on energy principle theory and Legendre polynomial series. The acoustic characteristics and the forced response of the coupled system are obtained using the Rayleigh–Ritz technique. In addition, the influence of geometric parameters on the acoustic characteristics of the coupled system is analyzed. The results show the following:(1)
The acoustic characteristics of a strongly coupled plate-cavity system can be obtained accurately and quickly, and since the L2 inner product of the Legendre polynomials satisfies orthogonality, the efficiency of the calculation is greatly enhanced.(2)
The change of the structural boundary from clamped supported to simply supported, the reduction of the cavity depth, and the increase in the length-width ratio of the structure all result in a shift of the coupled natural frequencies towards lower frequencies. The property can guide the design of acoustic products. Meanwhile, when the cavity depth is shallow, the coupled natural frequency is sensitive to the change of cavity depth, Therefore, special attention should be paid to cavity depth when designing and applying cavity structures.(3)
The steady-state response under the point force excitation, internal acoustic excitation, and a combination of both can be predicted accurately by the present method, which can provide a new idea for active noise control of the strongly coupled plate-cavity system.



---

*Source: 5565136-2021-05-31.xml*