Homogenization of the vibro-acoustic transmission on perforated plates

by   Eduard Rohan, et al.
University of West Bohemia

The paper deals with modelling of acoustic waves which propagate in inviscid fluids interacting with perforated elastic plates. The plate can be replaced by an interface on which transmission conditions are derived by homogenization of a problem describing vibroacoustic fluid-structure interactions in a transmission layer in which the plate is embedded. The Reissner-Mindlin theory of plates is adopted for periodic perforations designed by arbitrary cylindrical holes with axes orthogonal to the plate midplane. The homogenized model of the vibroacoustic transmission is obtained using the two-scale asymptotic analysis with respect to the layer thickness which is proportional to the plate thickness and to the perforation period. The nonlocal, implicit transmission conditions involve a jump in the acoustic potential and its normal one-side derivatives across the interface which represents the plate with a given thickness. The homogenized model was implemented using the finite element method and validated using direct numerical simulations of the non-homogenized problem. Numerical illustrations of the vibroacoustic transmission are presented.



There are no comments yet.



Homogenization of the vibro-acoustic transmission on periodically perforated elastic plates with arrays of resonators

Based on our previous work, we propose a homogenized model of acoustic w...

Modelling wave dispersion in fluid saturating periodic scaffolds

Acoustic waves in a slightly compressible fluid saturating porous period...

An acoustic model of a Helmholtz resonator under a grazing turbulent boundary layer

Acoustic models of resonant duct systems with turbulent flow depend on f...

Diffraction of acoustic waves by a wedge of point scatterers

This article considers the problem of diffraction by a wedge consisting ...

A numerical algorithm based on probing to find optimized transmission conditions

Optimized Schwarz Methods (OSMs) are based on optimized transmission con...

Homogenization of the fluid-saturated piezoelectric porous media

The paper is devoted to the homogenization of porous piezoelectric mater...

Acoustic passive cloaking using thin outer resonators

We consider the propagation of acoustic waves in a 2D waveguide unbounde...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.

1 Introduction

The noise and vibration reduction belongs to important issues in design of structures used in the automotive industry, or civil engineering. The engine silencer used to reduce the noise emitted by the exhaust gas presents an important and well known example. However, there are many similar solid structures which can influence the acoustic wave propagation in fluid. Usually they involve porous, or perforated plates, or panels, such that they are permeable for the gas flow. The straightforward approach to modelling the acoustic wave propagation through vibrating perforated plates consists in solving directly the vibroacoustic problem with a 3D elastic structure describing the plate. However, its numerical treatment using the finite element method can lead to an intractable problem because of the prohibitive number of DOFs corresponding to the geometric complexity of the perforated structure. Therefore, it is reasonable to replace the elastic plate by an interface on which coupling transmission conditions are prescribed.

In this paper, we consider the acoustic wave propagation in an inviscid fluid interacting with elastic structures designed as periodically perforated plates. The aim is to derive non-local vibro-acoustic transmission conditions using the periodic homogenization method. Although similar problems have been treated in the literature, cf. DELOURME201228, in this context, the plate elasticity has not been considered yet. As for the rigid structures, semi-empirical formulae for the acoustic impedance exist which were tuned by experiments, or developed using the electro-acoustic equivalent circuit theory Jung-etal-2007-JKPS; Sakagami2010; Stremtan2012, or the Helmholtz-Kirchhoff integral theory Zhou2013. During the last decade, a number of works appeared which are based on a homogenization strategy. For a thin rigid perforated plate represented by interface and characterized by the thickness it has been shown in bb2005; DELOURME201228 that this interface is totally transparent for the acoustic field at the zero order terms of the model which describes the limit behaviour for , cf. Dorlemann2017. For a higher order approximation, an approach based on the so-called inner and outer asymptotic expansions has been developed, such that two associated acoustic fields are coupled, one being relevant in the proximity of the perforations, the other at a distance from the limit interface, see e.g. Clayes-Delourm-AA2013; Marigo-Maurel-JASA2016; Marigo-Maurel-PRSA2016. In contrast with bb2005 dealing with thin perforated interfaces only, in rohan-lukes-waves07 we were concerned with homogenization of a fictitious layer in which rigid periodically distributed obstacles were placed. In particular, a rigid plate perforated by arbitrary shaped pores could be considered. Therein nonlocal transmission conditions were obtained as the two-scale homogenization limit of a standard acoustic problem imposed in the layer.

Here we follow the approach reported in rohan-lukes-waves07 to develop vibroacoustic transmission conditions which substitute the vibroacoustic interaction on an elastic perforated plate immersed in the acoustic fluid. Up to our knowledge, despite some numerical studies, see e.g. . Takahashi-Tanaka2002, a rigorous treatment of such a problem has not been treated using the homogenization method so far. As the result we obtain vibroacoustic transmission conditions in a form of an implicit Dirichlet-to-Neumann operator. Due to this operator, the elastic perforated plate can be replaced by an interface on which a jump of the global acoustic pressure is linked to the acoustic momenta associated with two faces of the homogenized plate. It allows us to obtain an efficient numerical model which takes into account geometrical details of the periodic perforation without need of discretizing the vibroacoustic problem at the global level. In other words, the homogenized interface provides a reduced model in which a complex 3D elastic structure is replaced by a 2D perforated plate model whose coefficients retain information about the perforation geometry. To do so, we rely on the homogenized Reissner-Mindlin plate tailor-made for the “simple” perforation represented by general cylindrical holes with axes orthogonal to the mid-plane of the plate. Elastic strongly heterogeneous plates were treated in rohan-miara-CRAS2011; rohan-miara-ZAMM2015 where the framework of the Reissner-Mindlin theory was used to derive a model of phononic plates, cf. Rohan2015bg-plates, but without the interaction with an exterior acoustic field.

The proposed modelling conception based on the problem decomposition and using the homogenization provides an alternative framework for modelling of microporous panels which are known for their capabilities of acoustic attenuation Toyoda-JSV2005; Zhou2013; LIU2016149. In Maxit-JASA2012 the so-called patch transfer functions were developed for numerical modelling of compliant micro-perforated panels.

The plan of the paper is as follows. In Section 2 the vibroacoustic problem of the wave propagation in a waveguide containing the perforated plate is decomposed into the problem in a fictitious transmission layer (the “in-layer” problem) and the “outer” problem governing the acoustic field out of the layer. The “in-layer” vibroacoustic problem is treated using the homogenization method in Section 3, where the local problems imposed in the representative periodic cell are introduced and formulae for the homogenized coefficients are given. In Section LABEL:sec-global, as the main result of this paper, the global acoustic problem is established using the limit “in-layer” and the “outer” problems which are coupled using additional conditions derived by an additional integration and averaging procedure. The limit two-scale model of the homogenized layer is validated in Section LABEL:sec-valid using direct numerical simulations of the original problem. Finally, in Section LABEL:sec-simul, the proposed model is employed to simulate wave propagation in a waveguide equipped with the perforated plate. Some technical auxiliary derivations are presented in the Appendix.


In the paper, the mathematical models are formulated in a Cartesian coordinate system where is the origin of the space and is a orthonormal basis for this space. The spatial position in the medium is specified through the coordinates with respect to a Cartesian reference frame

. The boldface notation for vectors,

, and for tensors,

, is used. The gradient and divergence operators applied to a vector a are denoted by and , respectively. By we denote the symmetrized gradient , i.e. the strain tensor. When these operators have a subscript which is space variable, it is for indicating that the operator acts relatively at this space variable, for instance . The symbol dot ‘’ denotes the scalar product between two vectors and the symbol colon ‘’ stands for scalar (inner) product of two second-order tensors. Throughout the paper, denotes the global (“macroscopic”) coordinates, while the “local” coordinates describe positions within the representative unit cell where is the set of real numbers. By latin subscripts we refere to vectorial/tensorial components in , whereas subscripts are reserved for the tangential components with respect to the plate midsurface, i.e. coordinates of vector represented by are associated with directions . Moreover, is the “in-plane” gradient. The gradient in the so-called dilated configuration with coordinates is denoted by . We also use the jump w.r.t. the transversal coordinate, .

Figure 1: Left: The global domain decomposition into and separated by the homogenized perforated plate represented by the interface . Center: detail of the layer structure; the layer thickness is proportional to the plate thickness and to the perforation period. Right: perforated interface and the representative periodic cell .

2 Formulation and decomposition of the vibroacoustic transmission problem

The aim of the paper is to find a representation of the vibro-acoustic interaction on a perforated plate. For this, homogenized vibroacoustic transmission conditions are derived using the asymptotisc analysis w.r.t. a scale parameter which has a double role: on one hand it deals with the thickness of an elastic plate when considered as a 3D object, on the other hand it describes the size and spacing of holes periodically drilled in the plate structure.

The flowchart of deriving the transmission conditions for a limit global problem consists of the following steps:

  • The vibro-acoustic problem (later called the “global problem”) is formulated in a domain in which the perforated elastic plate is embedded, being represented by a planar surface – the plate midsurface.

  • A transmission layer of the thickness is introduced in terms of which constitutes its midsurface. This will allow to decompose the global problem into two subproblems: the vibroacoustic interaction in the layer and the outer acoustic problems in . The two subproblems are coupled by natural transmission conditions on the “fictitious” interfaces .

  • We consider the layer thickness being proportional to the scale parameter, thus, , where is fixed. The asymptotic analysis considered for the problem in with the Neumann type boundary conditions on leads to the homogenized vibroacoustic transmission problem defined on . In this analysis, has the double role announced above and the plate is described using its 2D representation in the framework of the Reissner-Mindlin plate theory. In Remark 2.4 we explain the dual interpretation of the plate thickness used in the asymptotic analysis of the vibro-acoustic problem.

  • The final step is to derive the limit global problem for the acoustic waves in the fluid interacting with the homogenized perforated plate represented by . For this, with a few modifications we follow the approach used in rohan-lukes-waves07, where the rigid plate was considered; a given plate thickness corresponds to given finite thickness of the transmission layer. Then the continuity of the acoustic fields on interfaces yields the homogenized vibroacoustic transmission conditions which hold on .

Figure 2: Left: Transmission layer of thickness embedded in the global domain . Right: Interface representing the homogenized transmission layer.

2.1 Global problem with transmission layer

In this section, we introduce the problem of acoustic waves in a domain with embedded perforated elastic plate , see Fig. 2 and Fig. 3. The acoustic fluid occupies the domain . We consider a fictitious transimssion layer with a thickness , such that . The plate thickness is , while the layer thickness for a given fixed .

For a fixed parameter , correspondingly to and , we use a simplified notation and . The acoustic harmonic wave with the frequency is described by the acoustic potential in the fluid, the corresponding wave in the elastic body is described by the displacement field . Assuming the body is fixed to a rigid frame on the boundary and interacting with the fluid on , these fields satisfy the following equalities:


Above, is the sound speed in the acoustic fluid, is the stress in the lienar elastic solid, is reference fluid density, and by we denote the normal vector. The constants and are defined to describe incident, reflected, or absorbed acoustic waves in the fluid, according to a selected part of the boundary.

2.2 Geometry of the perforated layer

Given a bounded 2D manifold representing the plate mid-plane, we introduce , an open domain representing the transmission layer. This enables to decompose into three nonoverlapping parts, as follows: . Thus, the transmission layer is bounded by which splits into three parts:


where is the layer thickness and , see Fig. 1. In the context of the transmission layer definition, we consider the plate as a 3D domain defined in terms of the perforated midsurface ; the following definitions are employed:


where is the surface where the plate is clamped.

The midsurface representing the perforated plate is generated using a representative cell , as a periodic lattice. Let , where are given (usually ) and consider the hole , whereas its complement defines the solid plate segment. Then


Further we introduce the representative periodic cell and define its solid part ,


so that is the fluid part. Obviously, in the transmission layer , the fluid occupies the part


where and .

For completeness, by virtue of (3) we can introduce the decomposition of boundary . For this we need the boundary , where , so that the closed curve generates the cylindrical boundary :


For the sake of simplicity, by we shall refere to .

Figure 3: Scheme of the transmission layer in which the perforated plate (dark gray) is embedded. The complementary domain is occupied by the acoustic fluid (light gray).

2.3 Problem decomposition

The domain split allows us to decompose problem (1) into three parts. By we denote the acoustic potential in and , whereas is the acoustic potential in the transmission layer , see Fig. 3. Further, by we denote the acoustic fluid velocity projected into the normal of the interfaces . The following subproblems are considered:

  • 1,2: Given on , find defined in , such that


    where is the “external” boundary. As in problem (1), and are constants attaining values 0, or 1, whereas is the amplitude of an incident wave. In the context of a waveguide, we consider to be decomposed into three parts, , , and , denoting the input, the output and walls, respectively. By the constants in (8) different conditions on are respected: on the walls , whereas on and , on , which accounts for the non-reflection condition.

  • 3: Given on , find in and in , such that


    where is the surface of the elastic structure in contact with the fluid, thus, .

  • 4. For a fixed and , solutions to problems (8) and (9) are equivalent to the solution of (1), if the coupling conditions hold:


    where referes to normals outer to domains .

2.4 Plate model

The 3D model of an elastic plate involved in problem (1) can be replaced by a plate model which describes a thin structure. We assume a small for which the limit model of acoustic transmission can be interpreted. In this paper we shall approximate behaviour of the thin elastic structure by the Reissner-Mindlin (R-M) plate model, which allows us to consider the effects related to shear stresses induced by rotations of the plate crossections w.r.t. the mid-plane.

The R-M plate model can be obtained by the asymptotic analysis of the corresponding 3D elastic structure while its thickness . However, the obtained limit model is then interpreted in terms of a given thickness . We shall discus this point in Remark 2.4.

The plate is represented by its perforated mean surface , therefore all involved variables depend on . However, for a while we drop the superscript related to these variables. The plate deflections are described by amplitude of the membrane elastic wave , of the transverse wave and of the rotation wave . Two linear constitutive laws are involved, which depend upon the second order tensor , where is the shear coefficient, and the fourth order elasticity tensor which is given by the Hooke law adapted for the plane stress constraint; we define (all indices )


The Reissner-Mindlin plate model is derived using the following kinematic ansatz confining the displacement in a plate with the actual thickness ,


where is “membrane-mode” displacements, i.e. vector involves also the “transversal mode” (the deflection). The vector fields satisfy the following equations in ,


where describes the perforations. Above the applied forces ,

and moments

, depend on the acoustic potential . The crucial step in deriving the model of vibroacoustic transmission consists in describing these forces in terms of imposed on surface in the 3D plate representation.

Remark 1. In our asymptotic analysis of the acoustic transmission layer, we shall use the plate thickness in two contexts:

  • The periodically perforated plate model defined in terms of the 2D domain representing the mid-plane and the thickness with being fixed. In fact, for a given thickness and the perforation design (a given size of the holes yielding ) we can obtain .

  • To describe the interaction between the 3D elastic structure and the acoustic fluid, the thickness must be proportional to which is also related to the transmission layer thickness , thus, we consider and the elastic body occupying domain , see (3).

Thus, the homogenization of the periodically perforated plate is done by pursuing the asymptotic analysis applied to the 2D plate model (13) divided by . Whereas is fixed in the plate equation operator, beeing independent of , at the r.h.s. terms we get which is coherent with the dilation operation applied when dealing with fluid equation, see Section 2.7.

2.5 Variational formulation of the vibroacoustic problem in the layer

In order to derive the homogenized model of the transmission layer, we shall need the variational formulation of problem (9) with the plate model (13).

Find and such that


for all , where n is outward normal to domain , and


for all test functions . In (14), the displacements defined on the surface are expressed using the mid-plane kinematic fields. Due to (12), it holds that


where , . In analogy, the test displacements can bentroduced in terms of the test functions involved in (15); in this equation, the r.h.s. integrals express the virtual power


where the traction stress is induced by the acoustic pressure in the fluid.

2.6 Fluid structure interaction on the plate surface

The forces and moments involved in the r.h.s. of (15) can be identified using the 3D representation of the plate surface decomposed according to (3). The actual surface traction is given by the acoustic potential and by the surface normal ; note that , on , whereas on . Hence, it can be shown that the following expressions hold:


Then we consider the fluid equation. In (14), in the integral on , the displacement field w must be expressed in terms of the mid-plane displacements and rotations and , as introduced in (16). This yields


2.7 Dilated formulation

We can now state the vibro-acoustic problem in the dilated layer , where the fluid occupies domain , see (6). Using , while , with new coordinates , the gradients are , thus ; to simplify the notation, we shall use the same notation for functions depending on , but expressed in terms of .

By virtue of the dilation and the periodic unfolding, the vibroacoustic problem can be transformed in the domain which does not change with . Consequently the standard means of convergence can be used to obtain the limit model.

Equation (14) with the substitution (19) can now be transformed by the dilatation (the same notation for all variables is adhered, but should be interpreted in this new context of this dilated formulation):


Further we employ (18) to rewrite (15) which is divided by ; by virtue of Remark 2.4, the plate thickness is given, i.e. , however, when dealing with the r.h.s. interaction terms, in accordance with the dilation transformation. Thus we get


It is worth noting that, in (20) and (21), the r.h.s. integrals provide a symmetry of the following formulation.

The problem formulation

The vibroacoustic interaction in the dilated layer is described by which satisfy equations (20)-(21) for any test fields .

Let and , whereby being -periodic in the second variable; we define


For any and defined according to (22), the vibroacoustic interaction problem constituted by equations (20)-(21) possesses a unique solution . As an essential step of the proof, the a prioriestimates are derived in the Appendix LABEL:sec-appendixA.

3 Homogenization of the transmission layer

In this section, we introduce the convergence result which yields the limit acoustic pressure and the plate displacements and rotations. These are involved in the limit two-scale equations of the vibroacoustic problem imposed in the transmission layer. The asymptotic analysis is based on the unfolding method which was inaugurated in the seminal paper Cioranescu-etal-2008 and elaborated further for thin structures in Cioranescu2008-Neumann-sieve. In our setting, the unfolding operator transforms a function defined in into a function of two variables, and . For any , the cell average involved in all unfolding intergartion formulae will be abreviated by


whatever the domain of the the integral is (i.e. volume, or surface).

3.1 The convergence results

Based on the a priori estimates derived in the Apendix A, the following theorem holds.

Theorem 1. Let us assume


then the folloving estimates can be obtained: