On the Simulation of the Influence of Defects on Immersed Plane Periodic Multilayer Viscoelastic Media ()
1. Introduction
Periodic multi-layered media are widely studied because of their very interesting properties such as exhibiting absolute or local band gaps in their dispersion curves. The width and position of these bands depend on the physical and geometric characteristics of the materials and their frequency. Indeed, by the modification of the physical or geometric properties of a material constituting one or more periods, the width and the position of the forbidden or passbands can be modified, and very fine bandwidths can appear in the prohibited bands. It is thus possible to perform more or less selective wave filtering depending on the intended application [1] [2]. Numerous theoretical methods can be used to obtain the forbidden and passbands of the periodic media. Our theoretical approach is based on the stiffness matrix [3] method which ensures unconditional numerical stability, whatever the frequency range and the incidence angle. This method allows to study the forbidden and passbands from the calculation of the reflection and transmission coefficients. In the literature, several studies can be found concerning the study of periodic multilayer media [4] [5] [6]. For example, Lenoir et al. [4] [5], studied two cases of insonation and showed that, if we consider attenuation in one of the constitutive plates (polyethylene), the reflection coefficients are different depending on the insonation side while the transmission coefficients are identical. Maréchal et al. [7] showed that the number of periods necessary to obtain stable forbidden and passbands depends on the geometrical and physical characteristics of the viscoelastic plate. For their configurations studied, they demonstrated that three periods (N = 3) ensure the convergence of the reflection spectrum (or transmission spectrum). They have also shown that [8] the importance that should be given to the thickness measurement and the characterization of the viscoelastic properties of the materials used as input data for the multilayer structure. Siryabe et al. [6] compared results of multilayer without or with hole configurations and showed that throughout holes in the polyethylene plate allow the rapid observation of forbidden bands. Lenoir et al. [9] studied the influence of defect layers on the reflection coefficients on a plane multilayered fluid-loaded structure composed of N periods. Two types of inserted flaws are considered, either inside a period by varying the polyethylene layer thickness (without taking into account the attenuation), or by inserting a defect layer. The former case is based on experimental ascertainments which imply to take into account thickness variations. In the latter case, a layer whose acoustic impedance is different from those of the Al and PE plates plays the role of the defect. The defect layer is either in copper or in glass. Particular attention has been paid to the reflection minimums at low frequency, these ones being linked to the so-called vertical modes. They showed the linearity of the sensitivity of the vertical modes to thickness variations. In the case of the insertion of a defect layer, a very narrow passband emerge in a wide stopband, at a frequency depending on the thickness and acoustical properties of the defect layer. More recently, Khaled et al. [10] investigated the influence of the insertion of two geometric defect layers in a 1D periodic structure composed of N elementary cells, through the analysis of the band gaps and passbands of the transmission coefficient. The concept behind the introduction of a defect into the multilayer stacks is a way to create the so-called defect modes or resonant modes. The numerical study is carried out in the case of N-periods made of a Glass isotropic layer and another one of fluid water. The variations of the resonant modes in terms of the defect locations and thicknesses are shown as graphs.
The present study is the contribution to the previous work on multilayer periodic structure in which defects were introduced [9] [10]. We consider plane periodic multilayers composed of N periods (1 ≤ N ≤ 3). The elementary stack is made of two plates, one in aluminum and the other one in polyethylene, both having the same thickness, and presenting a high acoustic impedance contrast. Simulations are made at normal incidence. Contrary to the previous studies [6] [9] [10] concerning periodic structures, we will consider firstly, the presence of inclusion matter in the polyethylene plate, secondly the effect of boundary conditions different from perfect ones between the two plates composing a period and finally the presence of more than two layers in a stack. Theoretical basis used to obtain effective acoustic properties and reflection, transmission and absorption coefficients are described. Then, the results for the case of polyethylene containing grease inclusions configurations are compared with those of air porosities [6] and without porosities [4]. The effect of disbond simulated by the insertion of a thin Teflon layer and the number of layers composing a period are discussed. Finally, conclusion and some prospects are presented.
2. Theoretical Basis
2.1. Context
The geometry of the problem (Figure 1(a)) consists with a stacking of N periods (1 ≤ N ≤ 3), each one made up of two isotropic plates of different acoustic properties. As a case of study (Figure 1(b)), Siryabe et al. investigated the AL/PEp stack, where the polyethylene plate (PE) was perforated (air porosities) [6]. In this study, as an extension, we consider the case of grease inclusions. The aim is to observe the effect of the presence of grease on the reflection and transmission coefficients of the considered multilayer. Here, one will note (PEg) for polyethylene with grease. In both cases, the longitudinal and transverse waves attenuations are taken into account in the polyethylene plate but neglected in the aluminum one [4].
2.2. Homogenization Methods
Since the polyethylene layer was considered as containing grease inclusions after manufacturing, thus, it is suitable to determine its new effective properties using homogenization models as was done in [6]. In order to retain the appropriate model to our study, a comparison of the approaches suggested by Voigt [11] [12], Gaunaurd (G-sph) [13] and Hashin-Shtrikman (HS-iso, HS-ort) [14] [15]
Figure 1. (a) Geometry of the problem, (b) PE with spherical or cylindrical grease inclusions and (c) position of the insertion of thin Teflon plate.
is made. For a plate with inclusions, the estimation of effective properties can be done using the following relations [6]:
(1)
were {Keff, Geff} are given by homogenization models in Table 1.
The longitudinal and transverse effective velocities {cL,eff, cT,eff} which constitute the effective acoustic properties of a plate with inclusion, the effective density ρeff and the effective acoustic impedance ZL,eff are respectively expressed by [6]:
(2)
and
(3)
In a general way, depending on the homogenization model, the effective mechanical properties Xeff of a heterogeneous material can be written as follows:
(4)
Table 1. Effective mechanical properties according to the homogenization models.
ρ: Polyethylene’s density; K: Polyethylene bulk modulus; G: Polyethylene shear modulus; ρg: Grease’s density; Kg: Grease bulk modulus; Gg: Grease shear modulus; Vf: Volume fraction.
where {x, a, b, m} are either a constant or function depending on the homogenization model, associated to the elastic properties of the matrix X = {K, G, ρ} and Vf is the volume fraction of inclusions having the mechanical properties Xi = {Ki, Gi, ρi}.
Since the effective density ρeff is considered as that given in Equation (3), the effective properties {Keff, Geff} are summarized in Table 1 according to the different homogenization models for a polyethylene plate containing grease inclusions.
2.3. Comparison of Homogenization Models
In this part, the equivalent modulus curves according to the volume fraction Vf are plotted for the four models. The normalized {Keff, Geff} coefficients, i.e. {Keff/K, Geff/G} are shown to decrease when Vf increases in Figure 2(a) and Figure 2(b), respectively.
This reduction depends on the homogenization model. Moreover, we observed that the Hashin-Shtrikman model for orthotropic inclusions (HS-ort) and isotropic inclusions (HS-iso) gives values of effective mechanical properties lower than those obtained with other models. Consequently, the effective acoustic properties {cL,eff, cT,eff} required also will be lowest possible: what will make it possible as well as possible to describe the behavior of the wave in the structure with porosities (see Figure 3(a) and Figure 3(b), respectively).
For the polyethylene plate containing grease inclusions, the new values of the effective properties for longitudinal and transverse velocities (Equation (2)), density and acoustic impedance (Equation (3) for the four models are calculated. The results are presented in Table 2.
As we can see, the effective acoustic properties for the HS-ort model are very
Figure 2. Normalized elastic properties deduced from (Voigt, G-sph, HS-Iso, HS-ort) homogenization models according to the volume fraction Vf: (a) normalized effective bulk modulus Keff/K and (b) normalized effective shear modulus Geff/G.
Figure 3. Effective elastic and acoustic properties deduced from HS-ort homogenization model as a function of the volume fraction Vf: (a) normalized Young modulus, Poisson’s ratio and density {Eeff, neff, ρeff } and (b) normalized longitudinal and transverse velocities {cL,eff, cT,eff}.
Table 2. Effective acoustic properties according to the homogenization models for the polyethylene plate with 20% grease inclusion (PEg).
ρeff (kg/m3): Density; {cL,eff, cT,eff} (m/s): Longitudinal and transverse wave velocities; ZL,eff (MRa): Longitudinal acoustic impedance.
lower than effective acoustic properties for the HS-iso model. Thus, in the following, since orthotropic inclusions are considered, effective acoustic properties calculated using Hashin-Shtrikman (HS-ort) model will be retained and will be useful as input data for simulations.
2.5. Stiffness Matrix Method (SMM)
We consider two configurations of insonation: direct insonation (made from the aluminum layer side) and reverse insonation (from polyethylene side). For the direct case, the obtained reflection coefficient is noted Rd, and for the reverse one, the reflection coefficient is noted Rr. When attenuation is considered in the multilayer structure, those reflection coefficients are different according to the insonation side, either in the direct case Rd or in the reverse case Rr, whereas the transmission coefficient T is identical whatever the insonation side. These coefficients are obtained from a procedure suggested by Rokhlin et al. [3] [16]. For each period N, one defines a (4 × 4) flexibility matrix:
(5)
where Sij are 2 × 2 submatrices.
The Sij submatrices for N periods are deduced from the ones for (N − 1) periods and those of an additional n layer by the following recursive relationships:
(6)
It is shown that the expressions of Rd, Rrand T coefficients are expressed according to the elements (2, 2) of the
submatrices extracted from global compliance matrix SN, denoted as
into the forms:
(7)
where
.
3. Results and Discussions
In this section, the theoretical reflection, transmission and attenuation spectra are plotted according to the insonation side. The multilayer structure is immerged in a coupling medium (water in our case). In Table 3 are summarized the characteristics of each layer constituting the period and those of the coupling medium. Among acoustic properties, the effective properties of the air porosities polyethylene plate and those of the grease inclusions polyethylene plate are
Table 3. Acoustic properties of each layer constituting a period and the propagation medium.
ρ (kg/m3): Density; {cL, cT} (m/s): Longitudinal and transverse wave velocities; ZL (MRa): Longitudinal acoustic impedance.
calculated using the HS-ort homogenization model presented in Table 2. The considered volume fraction was evaluated on the basis of grease inclusions, distributed in a random way into the thickness of the polyethylene plate. The diameter of each inclusion was considered negligible compared to the plate thickness. The Poisson’s ratio of PE material is ν = 1/3 and the chosen volume fraction is Vf = 20%. The chosen bandwidth frequency is ranging in [0.75; 3.25] MHz.
The constitutive layers of periods have each one 4 mm thickness and are denoted by AL, PE, PEp and PEg for aluminum, polyethylene, perforated polyethylene and polyethylene with grease inclusions, respectively. A computation carried out with plates without losses will give an identical reflection spectrum whatever the insonation side [4] [5]. For that, we take into account longitudinal absorption in the polyethylene, by considering that the attenuation is linearly dependent on the frequency:
. This implies that the longitudinal wave velocity as a complex can be expressed as:
, where the losses are estimated around
in the polyethylene plate.
3.1. Grease Inclusions Effect
In previous works, Lenoir et al. [4] for (AL/PE) periodic structures and Siryabe et al. [6] for (AL/PEp) periodic structures showed that reflection coefficients reach their final shape as from the second period (N = 2). It is also the case of our configuration (AL/PEg) (Figure 4(a) and Figure 4(b)). In addition, the reflection coefficients minima tend towards zero and the undulations in both AL/PEp and AL/PEg configurations disappear for N = 1 period. A small shift and interchange is observed on local maxima and minima of reflection coefficients.
In the case of PEg/AL configuration (Figure 4(b)), undulations related to the vibrations of the polyethylene layer are obtained at the period N = 1. We also observed a very great decreasing in the amplitude of the reflection coefficients as the frequency increases. This result was also obtained by Siryabe et al. [6] in the case of PEp/AL. Concerning the minima of reflection; they tend towards zero for the three configurations whatever the frequency and the number of periods.
The difference of spectrum according to the insonation side is also observed
Figure 4. Reflection coefficients obtained with a (AL/PEg) stack, for N = 3 periods: (a) direct reflection coefficient |Rd|2, or (b) reverse reflection coefficient |Rr|2, and associated (c) transmission coefficient |T|2 and (d) reverse absorption coefficient |Ar|2.
in Figure 4(a) and Figure 4(b) for the configurations (AL/PEg) and (PEg/AL), respectively. From N = 1 period, the reflection minima linked to the propagating Lamb modes in the aluminum plate keep their positions for all configurations. These minima of reflection are observed at the frequencies 1.6, 2.4 and 3.2 MHz, in the bandwidth of the transducer.
The transmission coefficient illustrated by Figure 4(c), shows the forbidden band gaps as from N = 2 periods, similarly as for the previously cited works [4] [6]. More precisely, the transmission coefficients (Figure 4(c)) exhibit a clear reduction of its amplitudes when the frequency and the number of periods increase. Thus, most of the wave energy is attenuated as the frequency and the number of periods increases. For example, if we take the case of the first peak at 1.6 MHz, the amplitude is equal to 0.63, 0.41 and 0.25 for N = 1, 2 and 3 periods, respectively. This leads us to deduce the absorption by the following relationship:
(8)
As a result, the absorption coefficients |A{d,r}|2 also depend on the insonation side. The estimated absorption on the reverse side is equal to 0.33, 0.57 and 0.73 at 1.6 MHz (Figure 4(d)).
3.2. Volume Fraction Effect for N = 3
Siryabe et al., [6] observed an increase in the oscillations amplitude of the reflection coefficients when the volume fraction ratio and frequency increase (Figure 5(a)). In the case of our configuration (PEg/AL), for a fixed period N = 3 (Figure 5(b)), a decrease of the oscillations amplitude of reflection coefficients with the volume fraction and the frequency is observed since they are directly linked to the acoustic impedance contrast. This divergence can be explained by the fact that in the case of our configuration, the acoustic impedance contrast is lower than the case of study carried out by Siryabe et al., where air properties were neglected. In addition, we can observe that the minima are periodically spaced out by a shift of frequency Δf = cL/(2d) such as ΔfPEc = [267, 242, 221] kHz for PEg/AL configuration, corresponding to the longitudinal velocities cL,PEc = [2133, 1933, 1769] m/s, respectively associated to the volume fraction Vf = 20%, 40% and 60%.
3.3. Effect of Disbond Simulated by a Thin Teflon Plate
The insertion of a thin Teflon layer of which properties are given in Table 3, within two adjacent plates (AL and PE), increases the layer number of the period from two to three. The presence of 100 to 500 µm Teflon thickness is analyzed. We observed a shift to low frequencies on the local maxima and minima of reflection. The maximum shift is around 250 - 300 kHz. As illustrated by Figure 6, there are not significant changes depicted on the minima of reflection. By increasing the Teflon thickness, undulations on the transmission spectra increase too. Furthermore, reflection coefficients still have approximately same amplitude, whatever the number of period increases (Figure 6). About the location of the forbidden bands in the periodic structure, as we can see in Figure 6(c), three forbidden bands in the frequency range [0.75; 3.25] MHz are observed: from 1 to 1.5 MHz, 1.7 to 2.3 MHz and 2.5 to MHz. As an illustration, the influence of the variations of Teflon thickness layer on these passbands is depicted in Figure 6(c) where narrow passbands are observed in the frequency range from 1 to 1.5 MHz, reducing the forbidden bands.
Figure 5. Effects of volume fraction on the reflection coefficients (a) (PEp/AL) configuration [6] and (b) (PEg/AL) configuration for N = 3; Vf = 20%, 40%, 60%.
Figure 6. Effects of thin Teflon plate between AL and PE plate for N = 3 (a) direct (AL/PE) and (b) reverse (PE/AL) configuration [6].
3.4. Number of Layer of Period Effect: Reference Case of (AL/PEg/PEg/AL) × N and Symmetric Case of (PEg/AL/AL/PEg) × N
In the structure of type AL/PEg, the reflection coefficients reach their maximum from the second period (N = 2). It is not the case of the structure of type AL/PEg/PEg/AL, where the reflection coefficients do not evolve anymore and this from the period N = 1 (Figure 7(a) and Figure 7(b)). This is due to the fact that, the number of layers forming one period has been increased, but also to their arrangements. On the one hand, this result means that the reflection coefficient is related only to the first layers of the multilayer structure, the following ones being dissimulated by the attenuation effect. On the other hand, in the reverse configuration PEg/AL/AL/PEg (Figure 7(b)), the more the frequency increases, the more the amplitude of the reflection coefficients decrease. We also observe a duplication of the spectrum width with the apparition of the news minima of reflection, due to the double thickness of the aluminum layer.
Concerning the transmission coefficients, depending on the reference stack configuration (AL/PEg/PEg/AL) × N (Figure 7(a)) or the symmetric one (PEg/AL/AL/PEg) × N (Figure 7(b)), the insonation side influences significantly the evolution of the transmission coefficients. The forbidden band gap appears at a period N = 1 for the reference configuration (Figure 7(c)), and from N ≥ 2 periods in the symmetric case (Figure 7(d)). This leads us to conclude that, in the reference configuration case (AL/PEg/PEg/AL), the more the layer number
Figure 7. Coefficients obtained in (a), (b) reflection and (c), (d) transmission for the (a), (c) AL/PEg/PEg/AL and (b), (d) PEg/AL/AL/PEg configurations for N = 1, 2, 3; Vf = 20%.
of a period increases, the more the forbidden bands are observed in the multilayer. Moreover, amplitudes of the maxima of transmission coefficients decrease considerably when the frequency and the number of period increase. Regarding the number of peaks, it shifted from npeaks to 2npeaks for the reference configuration versus the symmetric one, respectively. In the reference configuration, the position of the maximum of transmission is the same as for the case of AL/Peg, what is to be understood as the resonances of the first AL layer. In the symmetric configuration, additional maxima are observed at 1.2, 2.0 and 2.8 MHz, in between to the three maxima of transmission of the reference stack, respectively. Moreover, in the reference configuration, the positions and widths of the forbidden bands are identical to those of the AL/PEg configuration.
4. Conclusion
In this study, the effect of grease inclusion, Teflon, volume fraction of grease and the number of layers composing one period were highlighted by the analysis of the reflection and transmission spectra of an immersed multilayered media. As in the porous case of AL/PEp, reflection coefficients reach their maximum as from the second period (N = 2), in the case of AL/PEg, the forbidden bands gap appears in the multilayer stack as from N = 2 periods. In addition, in the case of the studied (PEg/AL) configuration, for a fixed number of periods N = 3, a decrease of the oscillations amplitude of the reflection coefficients with the volume fraction and the frequency is observed, since they are directly linked to the acoustic impedance contrast. The presence of Teflon defects of 100 to 500 µm thick was analyzed. As a result, a shift is observed to low frequencies on the local maxima and minima of reflection. The maximum shift was around 250 - 300 kHz. Finally, concerning the transmission coefficients, unlike to the structures of type (AL/PE or PE/AL) [4], (AL/PEp or PEp/AL) [6] and (AL/PEg or PEg/AL) where the transmission spectra are identical whatever the insonation side (direct or reverse). In the case where the elementary stack consists of four layers, the forbidden band gaps appear at a period N = 1 for the reference configuration and from N ≥ 2 periods in the symmetric case. The number of peaks, shifted from npeaks to 2npeaks for the reference and symmetric configurations, respectively. Further studies are to be led experimentally in order to validate those numerical results. It would also be interesting to include piezoelectric layers as dynamic modulators of the reflection and transmission coefficients. As a perspective, such piezo-active and controllable multilayers have many potential applications as actuators and sensors in the industry. We could plan to evaluate the electromechanical response of the obtained global structure.