# Structural energy flow optimization through adaptive shunted piezoelectric metacomposites

Manuel Collet, Morvan Ouisse, Mohamed Ichchou

► To cite this version:

Manuel Collet, Morvan Ouisse, Mohamed Ichchou. Structural energy flow optimization through adaptive shunted piezoelectric metacomposites. Journal of Intelligent Material Systems and Structures, 2012, 23 (15), pp.1661-1677. 10.1177/1045389X12449915. hal-02300118

HAL Id: hal-02300118

https://hal.science/hal-02300118

Submitted on 3 Feb 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Structural energy flow optimization through adaptive shunted
piezoelectric metacomposites

M Collet, M. Ouisse

FEMTO-ST, Applied Mechanics Department,
UMR-CNRS 6174,

24 chemin de l'épitaphe,
25000 Besançon, France

M. N. Ichchou

LTDS UMR5513

Ecole Centrale de Lyon

36 avenue Guy de Collongue
69134 Ecully, France
In this paper, a numerical approach for modeling and optimizing 2D smart metacomposites is presented. The proposed methodology is based on the Floquet-Bloch theorem in the context of elastodynamics including distributed shunted piezoelectric patches. The dedicated numerical technique is able to cope with the multi-modal wave dispersions behaviour over the whole first Brillouin zone for periodically distributed 2D shunted piezo-mechanical systems. Some indicators allowing the optimization of the shunt impedance for specific performance objectives are presented, and applied for illustration purposes on the design of an adaptive metacomposite with specific functionalities. In order to validate the strategy, the designed metacomposite is integrated in a support structure, and a full 3D model is derived to illustrate the efficiency of the approach. This work is a collaborative effort supported by the French Research Agency under grant number NT09 – 617542.

Keywords: *Adaptive Metamaterial*, *Distributed Shunted piezoelectric patches*, *Wave propagation optimization*
# 1 Introduction

The design of innovative structures with high degree of performances is a permanent challenge for both industrial and academic communities. The issue is often posed in many engineering sectors such as civil engineering and transportations (ground and aeronautics). In the transportation sector the question becomes drastic with regards to energy savings concerns and with respect to sustainable developments (see reports of Clean Sky and DREAM EU Projects, AIAA's Emerging Technologies Committee - ETC). This is the case, for instance, in the aeronautic sector where a fantastic effort is made in order to integrate lightweight structures in the design of new aircrafts. Advanced materials and structures such as composite materials are thus more and more used. Such composite materials make the structures lighter and stiffer. They are now intensively included in the design and manufacturing of new aircrafts and space shuttles. However, such composite structures bring new vibration and acoustic problems. Indeed, compared to metallic structures it is noticeable that composites structures evolves specific noise emissions and particular safety issues. The majority of the projects dedicated to smart materials have focused on adding smart functionalities to existing structures - piezoelectric patches to plates and shells for acoustic or structural health monitoring as an example. Nonetheless, these types of solutions tend to address efficiently only one particular problem at the time, and require increased complexity in terms of wiring, controllers and modifications of the baseline structure to tackle successfully structural, acoustic control and structural health monitoring. Hence, among the problems to be considered, the broadband nature (from the low frequency LF up to the high frequency HF and including the mid frequency MF range) of vibration and noise associated to composite structures and their assembly is still an open and challenging concern. This paper is a contribution to the challenges of designing, optimizing and manufacturing a new class of truly integrated and intrinsic smart active metamaterial, structures and systems, also called metacomposite with ground breaking engineering performance in terms of mechanical and vibroacoustics behaviour.

The concept of metacomposite has a specific meaning in this paper. This concept couples two different aspects in vibration control. The first concept is connected to periodic structures theories usually connected to metamaterial developments. In this case, it is well known that the dynamic behaviour is fully connected to periodicity ratios and existing pass bands and blocked bands can be of real use in vibration control. The second concept is associated to vibration control through piezoelectric and smart materials. Specifically, shunted piezoelectric smart materials are employed for the metacomposite achievement by integrating into the metamaterial electronics and numerical components allowing implementation of adaptive and controlled behaviour. We also tend to extend the notion of programmable matter within the meaning of work presented in to vibroacoustic programming. The paper main novelty is then the design through full numerical analysis of a smart structure with broad band control abilities. Wave based methods and numerical simulation tools are adapted to the proposed concept. The realized composite controlled dynamic is in the LF low frequency and mid frequency MF band as well. Development of a wave trap metacomposite can be the ultimate goal of the reported research and results.

Looking at the literature, there are intensive involvements in both the modelling and control communities. Indeed, researchers studied numerous methods to deal with undesirable vibrations and noise. These approaches can be generally classified as three catalogues which are passive control, active control and semi-active or semi-passive control. Several passive control ways were developed. Tuned mass damper (TMD) or harmonic absorber is one of the most popular passive control approach used since many years in many different engineering applications (Ormondroyd & Den Hartog 1928, Den Hartog 1956, Spencer & Nagarajaiah 2003). The corresponding device is normally mounted onto a host structure to absorb the structural mechanical energy when the structure vibrates at its
resonant frequency. Passive damping technologies can also be considered as passive control strategies especially when optimization design is carried out (Milman & Cheng-Chih 1994, Baz 1997). The main disadvantages of these passive strategies are also their narrow band of efficiency and the necessity to introduce massive elements (damping materials or massive absorber) inside the considered structure. Optimal efficiency in low and broad frequency bands can be reached by using active technologies for vibration or noise control (Preumont 2002, Nelson & Elliott 1992, Banks, Smith & Wang 1996) with the drawback of introducing complex and massive numerical and electronics elements and external power supply systems.

Semi-active control or semi-passive control possesses partial advantages of both passive and active control. Piezoelectric shunt technologies are promising passive techniques for vibration control of flexible structures. Contrary to active control, the only external element to be used is a Passive Electrical Network (PEN) that is directly connected to the electrodes of the piezoelectric device. With this approach, the sensing element is not needed and the use of a passive network guarantees the stability of the coupled system. Furthermore, shunted piezoelectric systems offer a very attractive path for electro-mechanical integration and distribution. Numerous works have been published (Park & Palumbo 2004, Becker, Fein, Maess & Gaul 2006, Hurlebaus & Gaul 2006, Hagood & von Flotow 1991, Agnes 1994, dell'Isola, Maurini & Porfiri 2004, Moheimani & Fleming 2006, Cross & Fleeter 2002, Beck, Cunefare & Ruzzene 2008a) that present analyses of the capability and efficiency of a single shunted piezoelectrical patch for structural stabilization and wave cancellation. Hagood and VonFlotow (Hagood & von Flotow 1991) provided the first analytical formulation for passive shunt networks. They demonstrated how a piezoelectric patch shunted through a single resistive-inductive (RL) circuit acts as a vibration absorber tuned at the resonance frequency of the circuit. Since then, more complex shunting circuits have been investigated to extend the effectiveness over broader frequency bands: multi-mode techniques (Wu 1998, Wu & Bicos 1997), state switching (Corr & Clark 2002), synthetic impedance (Fleming, Belirens & Moheimani 2000) or negative impedance circuits (Behrens, Fleming & Moheimani 2001, Park & Palumbo 2004, Park & Baz 2005, Fukada, Date & Kimura 2002, Neubauer, Oleskiewicz, Popp & Krzyzynski 2006). Park and Palumbo decomposed the wave field on a beam to find a minimization of the reflected wave component (Park & Palumbo 2004). By noticing the link between the wave based approach of Park and Palumbo and the power input to the system, Cunefare developed a parameter selection technique based on the reactive input power (Cunefare 2006). Periodic structures have been shown to improve the broadband performance of control systems on structures (Casadei, Beck, Ruzzene & Cunefare 2009, Thorp, Ruzzene & Baz 2001, Thorp, Ruzzene & Baz 2005). On all these works, piezoelectric shunt damping also appears as promising passive or semi-passive technique for vibration control of flexible structures (Beck et al. 2008a, Cunefare 2006). Contrary to active control, the only external element to be used is a passive electrical network (PEN) that is directly connected to the electrodes of the piezoelectric device and can be integrated into the composite material (Forward 1978, Fleming, Behrens & Moheimani 2003, Fleming & Moheimani 2004). In this paper, thus shunted piezoelectric concept will be developed further for wave control purposes mainly. Unlike most of the existing papers, this contribution will concentrate on a deep multi-physical numerical modelling of composite structures coupled to shunted piezoelectric patches based on mathematical approach developed in (Collet, Ouisse, Ruzzene & Ichchou 2011). The numerical investigation is concerned with 2D wave treatment and optimization in such metacomposite. This is also an original aspect of the paper.

Considering numerical simulation of wave behaviours in elastodynamic media, two main strategies can be reported: the semi-analytical finite element (SAFE) and the wave finite element method (WFE). Both are of spectral nature but of different philosophies and practical use. Semi-analytical finite element method was investigated by many researchers in the past decades. This method was first proposed by (Nelson & Dong 1973) for plate waves to analyse natural vibrations of lam-
inated orthotropic plates in which lamina is discretized by finite element. The (cross-sectionally) two-dimensional problem was first developed by (Aalami 1973), which can be extended to the wave guides of arbitrary cross-sectional shape. Recently, Gavric calculated the dispersion relationship in a free rail by using triangular and quadrilateral elements, with those elements obtained from HamiltonŠs principle (Gavric 1995). Hayashi et al. (Hayashi, Song & Rose 2003) derived the SAFE formulation through virtual work principles and proposed a way to calculate the group velocity using the eigensolution at a given frequency. Damljanovic and Weaver (Damljanovic & Weaver 2004) developed the linear triangular elements for SAFE method using Lagrange equations to investigate the elastic waves in wave guides of arbitrary cross-section. Bartoli et al. (Bartoli, Marzani & Viola 2006) examined the wave mode characteristics in damped waveguides by using SAFE method, where the elements were also derived from Hamilton's principle. This method has also been adopted to investigate the wave propagation characteristics for thin-walled structures (Finnveden 2004), where the polynomial interpolation is used in the propagation axis and thus leads to the polynomial eigenvalue problems. The method is extended to curved structures in (Finnveden & Fraggstedt 2008), where an isoparametric element is used.

Another spectral method has attracted many researchers, which is based on the technique dealing with the wave propagation problems of periodic structures. Mead (Mead 1973) proposed a general theory in order to determine harmonic wave propagation characters, where both one-dimensional and two-dimensional periodic systems are considered. Recently, more studies focus on extending this idea to homogeneous structures rather than the generally periodic systems comprised of an arbitrary substructure (see (Houillon, Ichchou & Jezequel 2005, Bocquillet, Ichchou & Jezequel 2003)). This Wave Finite Element (WFE) method regards the homogeneous wave guide structure as a periodic system assembled by identical substructures. Thus the method for periodic systems can be seamlessly transplanted to the WFE method for eigenmode extraction. The identical substructures can be discretized with the aid of some commercial FE procedures, rather than the development of a relatively new FE code for specific elements. This allows the existing element libraries and powerful grid generation procedures to be used for many engineering structures. The WFE method is easy to apply due to its perfect connection with the standard FE method. This method has been used for structural vibration analyses (Mace, Duhamel, Brennan & Hinke 2005, Zhou, Ichchou & Mencik 2009, Hardy, Ichchou, Jezequel & Trentin 2009, Duhamel, Mace & Brennan 2006, Houillon et al. 2005), wave propagation in elastic wave guides (Mencik & Ichchou 2005b, Chen & Wilcox 2007), and calculations of the Green's functions for homogeneous acoustical and elastic media (Duhamel 2007). The method was also implemented for wave propagation and dynamic problems in homogeneous structures with internal fluid (Bocquillet et al. 2003, Mencik & Ichchou 2007, Maess, Wagner & Gaul 2006), where the studies were based on the WFE method for the one-dimensional wave propagation problem and concentrated on the lower frequency problems (far below the frequency of interest for NDT). One of the main problems all these approaches are facing is the difficulty to compute the damped wave numbers in complex media. Precisely, wave properties in the whole Brillouin domain for multiphysics problem incorporating shunted piezoelectric patches is not an easy task. This involvement is necessary for vibroacoustic behaviour optimization purposes of such periodic smart structures (Collet, Ouisse, Ruzzene & Ichchou 2011, Collet, Ouisse, Ichchou & Ruzzene 2011).

The objectives of the work reported in this document were threefold. First, the paper offers a suited numerical approach for computing wave dispersion in smart multiphysics metacomposite with shunted piezoelectric patches. Second, the achieved model is extensively used for electromechanical optimization studies of the metacomposite, and finally, the optimized configuration is used in a finite structure to check the applicability of the concept on a practical application. The paper is structured in four different sections. After recalling the Floquet-Bloch theorems in section 2, we introduce a numerical formulation to compute the multi-modal damped wave numbers dispersion
in the whole first Brillouin domain. Based on this wave modelling, optimization of the electrical impedance of the shunted circuit is made in order to control group velocity of flexural waves or to increase damping induced by the electric circuit. The corresponding optimal impedances are then tested for controlling the mid frequency (MF) response of a semi-distributed finite system. The numerical achievements prove the effectiveness of the proposed metacomposite concept and draw future steps and perspectives.

## 2 Piezo-Elasto Dynamical Application of the Floquet-Bloch Theorem

In this section the application of the celebrated Floquet-Bloch theorem is presented for piezoelastodynamic problems. Based on the well known results obtained by Floquet (Floquet 1883) for monodimensional problems and later rediscovered by Bloch (Bloch 1928) in multidimensional problems, an original application to bi-dimensional piezo-elastodynamical problem has been proposed recently (Collet, Quisse, Ichchou & Ruzzene 2010). This formulation leads to very general numerical implementation for computing waves dispersion for periodically smart distributed mechanical systems incorporating electronic components, damping effects or any frequency-dependent characteristics. The main ideas of the approach are recalled here, together with the specific points related to the inclusion of piezoelectric effects in the model.

### 2.1 The Bloch theorem

The Bloch theorem, in its original version, gives the form of homogeneous states of Schrödinger equation with periodic potential. This theorem can be considered as a multidimensional application of the Floquet theorem (Joannopoulos, Meade & Winn 1995). The periodic medium (or potential) properties satisfy $M(x+R.m) = M(x)$, $m \in \mathbb{Z}^3$ where $R = [r_1, r_2, r_3] \in \mathbb{R}^{3 \times 3}$ is a matrix grouping the three lattice's basis vectors (in 3D). The primitive cell is defined as a convex polyhedron of $\mathbb{R}^3$ called $\Omega_x$. The reciprocal unit cell is denoted by $\Omega_k$ and limited by the reciprocal lattice vector defined by the three vectors $g_j$ so that: $r_i \cdot g_j = 2\pi\delta_{ij}$ ($\delta_{ij}$ being the Kronecker index). If $\Omega_x$ is the irreducible primitive cell, $\Omega_k$ corresponds to the first Brillouin zone of the lattice (Kittel 1986).

The Bloch Theorem stipulates that any function $u(x) \in L^2(\mathbb{R}^3, \mathbb{C}^n)$ can be expressed as

$$u(x) = \int_{\Omega_k} e^{ikx} \tilde{u}(x, k) dk \quad (1)$$

where the Bloch amplitude $\tilde{u}(x, k)$ is $\Omega_x$-periodic and has the representations

$$\begin{aligned} \tilde{u}(x, k) &= \sum_{n \in \mathbb{Z}^3} \hat{u}(k + Gn)e^{iGn \cdot x}, \\ &= \frac{|\Omega_x|}{(2\pi)^3} \sum_{n \in \mathbb{Z}^3} u(x + Rn)e^{ik(x+Rn)} \end{aligned} \quad (2)$$

where $\hat{u}(k)$ stands for the Fourier transform of $u(x)$ and $G = [g_1, g_2, g_3]$ is the reciprocal lattice matrix in the later. It can also be demonstrated that the mean value of the Bloch amplitude is the Fourier amplitude of $u(x)$ for the corresponding wave vector: $\langle \tilde{u}(., k) \rangle_{\Omega_x} = \hat{u}(k)$. Using the
Bloch theorem to represent the solutions of periodical partial derivative equations implies that all derivatives are shifted by $k$ in the sense given by the spatial operator.

Based on that theorem, the expansion functions $v_m(x, k)$ can be defined. They are called the Bloch eigen modes, and can be used to represent the Bloch amplitudes of any solution of the corresponding partial derivative equation as

$$ \tilde{u}(x, k) = \sum_{m} u_m(k) v_m(x, k) \qquad (3) $$

and at the same time diagonalize the partial derivative equations. The expansion coefficients $u_m(k)$ depend on the disturbance and on the induced wave vector (see (Bensoussan, Lions & Pananicolaou 1978) for details).

## 2.2 Application to piezo-elastodynamic

<figure><img src="image_2.png" /><figcaption>Figure 1: Generic 3D piezocomposite periodic cell</figcaption></figure>

Let us consider a piezo-elastodynamic problem made of infinite periodic distribution of unitary cell described in figure 1. The harmonic homogeneous dynamical equilibrium of system is driven by the following partial derivative equations ($t \in \mathbb{R}$):

$$ \begin{cases} \rho \ddot{w}(x, t) - \nabla \sigma(x, t) = 0 & \forall x \in \Omega_x \\ -\nabla D(x, t) = 0 & \forall x \in \Omega_x \end{cases} \qquad (4) $$

where $w(x, t) \in \mathbb{R}^3(\Omega_x, \mathbb{R})$ is the displacement vector, $\sigma$ represents the Cauchy stress tensor, $\epsilon = \nabla_{sym}w = \frac{1}{2}(\nabla w^T(x, t) + w(x, t)\nabla^T)$ the Green strain tensor, $D(x, t)$ the electric displacement. The linear constitutive material behavior relationships can be written as

$$ \sigma = C_E(x)\epsilon - e^T(x)E, \qquad (5) $$

$$ D = e(x)\epsilon + \epsilon_S(x)E, \qquad (6) $$
where $E = -\nabla V$ is the electric field vector ($V$ being the voltage), $C_E$ the elasticity tensor at constant electrical field, $e^T$ the piezoelectric coupling tensor and $\varepsilon_S$ the dielectric permittivity at constant strain. We add to this set of equilibrium equations an output expression

$$q^o = - \int_{S_t} \mathbf{D} \cdot \mathbf{n} dS \qquad (7)$$

allowing the introduction of the charge measurement on the piezoelectric's top electrode and hence the dual counterpart of the imposed electrical Dirichlet boundary condition for applying the shunt impedance operator.

The equations above are consistent for each kind of material to the extent that null piezoelectric and permittivity tensors can be used when passive materials are considered. All of these tensors also depend on the spatial location vector $x$ and are $\Omega_x$ periodic. By applying a Fourier transformation, the piezo-elastodynamic equilibrium can also be written as ($\omega \in \mathbb{R}$):

$$\rho\omega^2\mathbf{w}(\mathbf{x}, \omega) + \nabla C \nabla_{sym}(\mathbf{w}(\mathbf{x}, \omega)) + \nabla e^T(\mathbf{x})\nabla V(\mathbf{x}, \omega) = 0 \quad \forall \mathbf{x} \in \Omega_x \qquad (8)$$

$$-\nabla e(\mathbf{x})\nabla_{sym}(\mathbf{w}(\mathbf{x}, \omega)) + \nabla \varepsilon_S(\mathbf{x}, \omega)\nabla V(\mathbf{x}, \omega) = 0 \quad \forall \mathbf{x} \in \Omega_x \qquad (9)$$

As the problem is 2D infinitely periodic, mechanical boundary conditions are included in the formulation, while electrostatic boundary conditions have to be considered on each cell ($\omega \in \mathbb{R}$):

$$\left\{ \begin{array}{ll} V(\mathbf{x}, \omega) = 0 & \forall \mathbf{x} \in S_b \\ V(\mathbf{x}, \omega) = V^o(\omega) & \forall \mathbf{x} \in S_t \\ \mathbf{D} \cdot \mathbf{n}(\mathbf{x}, \omega) = 0 & \forall \mathbf{x} \in S_l \end{array} \right. \qquad (10)$$

where $S_b$ is the grounded bottom electrode of the piezoelectric layer, $S_t$ is the top electrode connected to the external shunt and $s_l$ the lateral electrode. The top electrode applied feedback voltage $V^o$ depends on the shunt characteristic and on the collected charges $q^o$ given in eq. (7) and can be expressed in the Fourier space by:

$$V^o(\omega) = -Z(\omega)q^o(\omega) \qquad (11)$$

In the following, the frequency dependency on $\omega$ of each functions in the Fourier space is removed for simplifying the notations. Considering a primitive cell of the periodic problem $\Omega_x$, the Bloch eigenmodes (3) and the dispersion functions can be computed by searching the eigen solutions of the homogeneous problem (8) and (9) with mechanical periodic boundary conditions and electric ones given by (10) as:

$$\mathbf{u}(\mathbf{x}) = \begin{bmatrix} \mathbf{w}(\mathbf{x}) \\ V(\mathbf{x}) \end{bmatrix} = \mathbf{u}_{n,k}(\mathbf{x})e^{i\mathbf{k}\cdot\mathbf{x}} \qquad (12)$$

where $\mathbf{u}_{n,k}(\mathbf{x}) = \begin{bmatrix} \mathbf{w}_{n,k}(\mathbf{x}) \\ V_{n,k}(\mathbf{x}) \end{bmatrix}$ is composed of $\Omega_x$ periodic functions. By introducing expression (12) in the piezo-elastodynamic equations (8) and (9), $\mathbf{w}_{n,k}(\mathbf{x})$, $V_{n,k}(\mathbf{x})$ and $\omega_n(\mathbf{k})$ can be found by solving the generalized eigenvalues problem:

$$\begin{aligned} 0 ={}& \rho \omega_n^2(\mathbf{k}) \mathbf{w}_{n,k}(\mathbf{x}) + \nabla C \nabla_{sym} \mathbf{w}_{n,k}(\mathbf{x}) + \nabla e^T \nabla V_{n,k}(\mathbf{x}) \\ & +ik \left[ (C \nabla_{sym} \mathbf{w}_{n,k}(\mathbf{x})).\mathbf{\Phi} + \nabla(C \Xi_{n,k}(\mathbf{x})) \right] \\ & +ik \left[ (\nabla e^T V_{n,k}(\mathbf{x})).\mathbf{\Phi} (e^T \nabla V_{n,k}(\mathbf{x})).\mathbf{\Phi} \right] \\ & -k^2 \left[ (C \Xi_{n,k}(\mathbf{x})).\mathbf{\Phi} + V_{n,k}(\mathbf{x}) (e^T \mathbf{\Phi}).\mathbf{\Phi} \right] \quad \forall \mathbf{x} \in \Omega_x, \end{aligned} \qquad (13)$$

$$\begin{aligned} 0 ={}& -\nabla e \nabla_{sym} \mathbf{w}_{n,k}(\mathbf{x}) + \nabla \varepsilon_S \nabla V_{n,k}(\mathbf{x}) \\ & -ik \left[ \nabla (e \Xi_{n,k}(\mathbf{x})) + (e \nabla_{sym}(\mathbf{w}_{n,k}(\mathbf{x}))).\mathbf{\Phi} \right] \\ & +ik \left[ (\nabla \varepsilon_S V_{n,k}(\mathbf{x})).\mathbf{\Phi} + (\varepsilon_S \nabla V_{n,k}(\mathbf{x})).\mathbf{\Phi} \right] \\ & +k^2 \left[ (e \Xi_{n,k}(\mathbf{x})).\mathbf{\Phi} - (s \mathbf{\Phi} V_{n,k}(\mathbf{x})).\mathbf{\Phi} \right] \quad \forall \mathbf{x} \in \Omega_x, \end{aligned} \qquad (14)$$
with the associated boundary conditions:

$$
\left\{
\begin{array}{ll}
\boldsymbol{w}_{n,k}(\boldsymbol{x}-R.m) = \boldsymbol{w}_{n,k}(\boldsymbol{x}) & \forall \boldsymbol{x} \in S_r, \boldsymbol{m} \in \mathbb{Z}^2, \\
V_{n,k}(\boldsymbol{x}) = 0 & \forall \boldsymbol{x} \in S_b, \\
V_{n,k}(\boldsymbol{x}) = -Z(i\omega)q_{n,k}^o & \forall \boldsymbol{x} \in S_t, \\
\boldsymbol{D}.\boldsymbol{n} = 0 & \forall \boldsymbol{x} \in S_l.
\end{array}
\right.
$$

In these equations $k = k \begin{bmatrix} \cos(\phi) \\ \sin(\phi) \\ 0 \end{bmatrix}$ is the $k \Phi$ where $\Phi$ represents the direction angles into the reciprocal

2D lattice domain and $\Xi_{n,k}(\mathbf{x}) = \frac{1}{2}(\mathbf{w}_{n,k}(\mathbf{x}) \cdot \mathbf{\Phi}^T + \mathbf{\Phi} \cdot \mathbf{w}_{n,k}^T(\mathbf{x}))$ is the symmetric dyadic tensor or the dyadic product of the displacement $\mathbf{w}_{n,k}(\mathbf{x})$ and direction vector $\mathbf{\Phi}$. $S_r$ are the interfaces of the cells continuum, and $R$ is the matrix grouping the two lattice's basis vectors (in 2D in the considered problem). In the electrical boundary conditions, $q_{n,k}^o$ is given by:

$$ q_{n,k}^{o} = \int_{S_t} \left[ \begin{aligned} & -e \nabla_{\text{sym}} \mathbf{w}_{n,k}(\mathbf{x}) + ik \nabla e \Xi_{n,k}(\mathbf{x}) \\ & + \varepsilon_S (\nabla V_{n,k}(\mathbf{x}) + ik V_{n,k}(\mathbf{x}) \mathbf{\Phi}) \end{aligned} \right] .n dS \quad (15) $$

where $\mathbf{n}$ is the surface outpointing unitary normal vector.

The proposed formulation is then based on the computation of the Floquet vectors given in equations (13) and (14)), instead of computing the Floquet propagators commonly used for elastodynamic applications. The full 2D waves dispersions functions can then be obtained, while damping and electrical impedance can clearly be introduced into the piezo-elastodynamic operator. The adopted methodology allows the computation of the complete complex map of the dispersion curves incorporating computation of evanescent waves and allowing the introduction of damping and shunt operator if any.

## 2.3 Weak formulation and computation of waves dispersion functions in periodic piezo-composite lattice

Let us consider the partial derivative equations (13), (14) on a unit cell $\Omega$. It stands for a generalized eigenvalue problem leading to compute the dispersion functions $\omega_n(\mathbf{k})$ and the corresponding Floquet eigenvectors $\mathbf{u}_{n,\mathbf{k}}(\mathbf{x})$. For computing the 2D dispersions curves, we need to introduce a suitable weak formulation.

If $\mathbf{u}_{n,\mathbf{k}}(\mathbf{x})$ is a solution of equations (13), (14), then:

$$ \forall \tilde{\mathbf{w}}_{n,k}(\mathbf{x}) \in \{H_1(\Omega, \mathbb{C}^3)/\tilde{\mathbf{w}}_{n,k}(\mathbf{x} - Rm) = \tilde{\mathbf{w}}_{n,k}(\mathbf{x}) \quad \forall \mathbf{x} \in S_r\} \text{ and} $$

$$ \forall \tilde{V}_{n,k}(\mathbf{x}) \in \{H_1(\Omega, \mathbb{C})/\tilde{V}_{n,k}(\mathbf{x}) = 0 \quad \forall \mathbf{x} \in S_b \text{ and } \tilde{V}_{n,k}(\mathbf{x}) = V^o \quad \forall \mathbf{x} \in S_t\} \text{ one has:} $$

$$ \begin{aligned} 0 = & \int_{\Omega} \left[ \begin{aligned} & \rho \omega_n^2(\mathbf{k}) \tilde{w}_{n,k}(\mathbf{x}) \mathbf{w}_{n,k}(\mathbf{x}) \\ & - (\tilde{\varepsilon}_{n,k}(\mathbf{x}) - ik\tilde{\Xi}_{n,k}(\mathbf{x})) C(\varepsilon_{n,k}(\mathbf{x}) + ik\Xi_{n,k}(\mathbf{x})) \\ & + (\tilde{\varepsilon}_{n,k}(\mathbf{x}) - ik\tilde{\Xi}_{n,k}(\mathbf{x})) e^T(\nabla V_{n,k}(\mathbf{x}) + ikV_{n,k}(\mathbf{x})\mathbf{\Phi}) \\ & - (\nabla \tilde{V}_{n,k}(\mathbf{x}) - ik\tilde{V}_{n,k}(\mathbf{x})\mathbf{\Phi}) e(\varepsilon_{n,k}(\mathbf{x}) + ik\Xi_{n,k}(\mathbf{x})) \\ & + (\nabla \tilde{V}_{n,k}(\mathbf{x}) - ik\tilde{V}_{n,k}(\mathbf{x})\mathbf{\Phi}) \varepsilon_S(\nabla V_{n,k}(\mathbf{x}) + ikV_{n,k}(\mathbf{x})\mathbf{\Phi}) \end{aligned} \right] d\Omega \\ & - \frac{\tilde{V}_{n,k}^t V_{n,k}^t}{Z(i\omega)} \end{aligned} \quad (16) $$
This weak formulation is simply obtained by integrating equation (13), (14) projected onto any test function $\tilde{u}_{n,k}(x)$. The boundary integral vanishes as the test functions are chosen so that $\tilde{w}_{n,k}(x - Rm) = \tilde{w}_{n,k}(x)$ on $S_r$. For a polyhedron cell, each boundary is generally a polyhedral plane sub-domain that can be associated with a parallel opposite one. The symmetry conditions called $\tilde{w}_{n,k}(x - Rm) = \tilde{w}_{n,k}(x)$ explicitly link these associated surfaces.

## 2.4 Numerical computation of the Bloch's waves

The numerical implementation is obtained by using a standard finite elements method to discretize the weak formulation (16). The assembled matrix equation is given by:

$$0 = [K(Z(\omega_n)) + \lambda L(\phi, Z(\omega_n)) - \lambda^2 H(\phi, Z(\omega_n)) - \omega_n^2 M] \mathbf{u}_n(k, \phi), \quad (17)$$

where $\lambda = ik$, $M$ and $K(Z(\omega_n))$ are respectively the standard symmetric semi-definite mass and stiffness matrices (the mass matrix is semi definite because elastostatic equation are condensed into the equation), $L(\phi, Z(\omega_n))$ is a skew-symmetric matrix and $H(\phi, Z(\omega_n))$ is a symmetric semi-definite positive matrix. $\mathbf{u}_n$ is the generalized eigen vector defined on all degrees of freedom of the used finite element model.

When $k$ and $\phi$ are fixed and $Z$ does not depend on $\omega$ the system (17) is a linear eigen value problem allowing us to compute the dispersion functions $\omega_n^2(k, \phi)$ and the associated Bloch eigen-vector $\mathbf{u}_n(k, \phi)$.

This approach has been widely used for developing homogenization techniques and spectral asymptotic analysis like in the work of (Allaire & Congas 1998). It can also be applied for computing wave's dispersion even if Floquet propagators is preferred for 1D or quasi 1D computation, as indicated in (Ichchou, Akrout & Mencik 2007, Mencik & Ichchou 2005a). Nevertheless these approaches have been only developed for undamped or lightly damped mechanical systems. In these cases, most of the previously published works present techniques based on the mesh of a real $k$-space (i.e $k$ or $k$ and $\phi$) following the boundary of the first Brillouin zone for obtaining the corresponding dispersion curves and the associated Floquet vectors. For undamped system, only propagative or evanescent waves exist corresponding to a family of eigen solutions purely real or imaginary. Discrimination between each class of waves is easy. If a highly damped system ($K, L, H$ are complex frequency dependent) and a frequency dependent electrical shunt impedance are considered, the obtained eigenvalue problem is not quadratic and a complex specific numerical methodology has to be implemented. Furthermore, evanescent parts of propagating waves appear as the imaginary part of $\omega_n^2(k, \Phi)$. It then becomes much more difficult to distinguish the propagative and evanescent waves as all solution appear complex.

Another much more suitable possibility for computing damped system, dedicated for time/space deconvolution and for computation of diffusion properties as defined by (Collet, Ouisse, Ruzzene & Ichchou 2011, Collet, Cunefare & Ichchou 2009, Mencik & Ichchou 2005a), is to consider the following generalized eigen value problem:

$$0 = (K(Z(\omega) - \omega^2 M + \lambda_n(\omega, \phi)L(\phi, Z(\omega)) - \lambda_n^2(\omega, \phi)H(\phi, Z(\omega)))\mathbf{u}_n(\omega, \phi). \quad (18)$$

In this problem, the pulsation $\omega$ is a real parameter corresponding to the harmonic frequency. Wave's numbers and Floquet vectors are then computed. An inverse Fourier transformation in the $k$-space domain can lead us to evaluate the physical wave's displacements and energy diffusion operator
when the periodic distribution is connected to another system as in (Collet et al. 2009). Another temporal inverse Fourier transformation can furnish a way to access spatio-temporal response for non-homogeneous initial conditions. As L is skew-symmetric, the obtained eigen values are quadruple $(\lambda, \bar{\lambda}, -\lambda, -\bar{\lambda})$ collapsing into real or imaginary pairs (or a single zero) when all matrices are real (i.e. for an undamped system). In this case a real pair of eigen values correspond to evanescent modes oriented in two opposite directions on the k-space and imaginary values to two traveling waves propagating in opposite direction as described in (Garvic 1995).

As previously mentioned, the real part of $k = k\Phi$ vector is restricted to stand inside the first Brillouin zone. In the quadratic eigen value problem (18) nothing restricts computation to only find eigen values satisfying this condition. For direction vector $\Phi$ orthogonal to the lattice facelets (i.e. for $\Phi_{p1} = [1,0]^T$ and $\Phi_{p2} = [0,1]^T$ in bi-dimensional rectangular cell), the periodical conditions expressed for one dimensional wave guide are still valid: if $\lambda_j(\omega, \Phi_p)$ is an eigen value associated to $w_j(\omega, \Phi_p)$ then $\forall m \in \mathbb{Z}^3$, $\lambda + i \cdot \Phi_p^T(G \cdot m)$ is also an eigen value associated to $w_j(\omega, \Phi_p)e^{-i \cdot \Phi_p^T(G \cdot m)x}$. Thus, for undamped systems, all obtained eigenvalues are periodically distributed in the k-space along its principal directions.

## 2.5 Computation of the evanescence and damped power flow criteria

One aim of this paper is to provide a numerical methodology for optimizing the piezoelectric shunt impedance $Z(\omega)$ for controlling energy flow into the periodically distributed piezo-composite structure. For doing this, we need to define suitable criteria evaluating the efficiency of the approach on the features of interest.

The first criterion which is considered for describing the capability of the metacomposite for transmitting structural energy is based on the computation of the wave group velocities. Indeed, they indicate how energy is transported into the considered system and allow to distinguish the 'propagative' and 'evanescent' waves. If a Bloch eigen solution (i.e $u_n(\omega, \phi)$, $k_n(\omega)$) is considered, the associated group velocity vector (Maysenhölder 1994) is given by:

$$C_{g_n}(\omega, \phi) = \nabla_k \omega = \frac{\langle\langle S \rangle\rangle}{\langle\langle e_{tot} \rangle\rangle} = \frac{\langle I \rangle}{\langle E_{tot} \rangle} \quad (19)$$

where $\langle\langle:\rangle\rangle$ is the spatial and time average respectively on one cell and one period of time, $S$ is the density of energy flow, $I$ the mean intensity and $e_{tot}$, $E_{tot}$ the total piezomechanical energy and its time average on a period (see (Maysenhölder 1994) for details). In this problem, we only consider mechanical energy transportation as the electrostatic coupling is decentralized on each cell and can not induce spatial energy transportation. The electromechanical coupling effect is also condensed into the mechanical behaviour as explicitly introduced in (Hagood & von Flotow 1991) and extended to finite element model in (Collet & Cunefare 2008).

The intensity vector $I$ is expressed as:

$$\langle I_n \rangle = -\frac{\omega}{2} Re \left( \int_{\Omega_x} C(\varepsilon_n(\mathbf{x}) + ik\Xi_n(\mathbf{x})).(\mathbf{w}_n^*(\mathbf{x})) \frac{d\Omega}{V_{\text{ol}}} \right) \quad (20)$$

where $.*$ is the complex conjugate, $Re$ stands for real part and $V_{ol}$ the domain volume.

As the spatio-temporal average of the system Lagragian is null (see (Maysenhölder 1994)), the total energy average is approximated by only computing the kinetic energy average:

$$\langle E_{tot} \rangle = \frac{1}{2V_{ol}} \left( \int_{\Omega_x} \rho \omega^2 \mathbf{w}_n(\mathbf{x}, \omega, \phi) \cdot \mathbf{w}_n^*(\mathbf{x}, \omega, \phi) d\Omega \right). \quad (21)$$
The group velocity vectors $C_{g_n}(\omega, \phi)$ are computed for all wave numbers at each frequency. In order to focus our analysis on only flexural modes (S and SH ones) we introduce an indicator allowing to select them by computing the ratio of kinetic energy average on out of plane displacement as:

$$Ind(n, \omega, \phi) = \frac{\frac{1}{2V_{ol}} \int_{\Omega_x} \rho \omega^2 w z_n(\mathbf{x}) w z_n^*(\mathbf{x}) d\Omega}{\langle E_{tot} \rangle}, \quad (22)$$

$wz_n(\mathbf{x})$ being the (Oz) component of vector $w_n(\mathbf{x})$. Optimization of the shunt impedance $Z(\omega)$ is based on the minimization of the maximal group velocity collinear to the wave number vector (19) for waves having a ratio of transported flexural kinetic energy (22) greater than 0.8. The criterion can then be written as:

$$\mathrm{Crit}_1(Z(\omega), \phi) = \max_{n/\mathrm{Ind}(n,\omega,\phi)>0.8} (C_{g_n}(\omega, \phi) \cdot \mathbf{\Phi}). \quad (23)$$

The second criterion for optimizing the structural energy absorption is based on the maximization of the damped electric power expressed as the active electrical power

$$P_{elec}(n, \omega, \phi) = \frac{1}{2} \mathrm{Re}(i\omega Z(\omega) q_{n,k}^o q_{n,k}^{o*}). \quad (24)$$

In order to increase the damping effects inside the smart metacomposite, this term needs to be sufficiently large. In the second case, the criteria is then

$$\mathrm{Crit}_2(Z(\omega), \phi) = \max_{n/\mathrm{Ind}(n,\omega,\phi)>0.8} \frac{1}{P_{elec}(n, \omega, \phi)}. \quad (25)$$

## 3 Optimization of the flexural energy flow inside the shunted periodic piezo-composite

The considered piezo-composite cell is presented in figure 1. The supporting plate material is standard aluminum with 0.1 % of hysteretic damping ratio and the piezoelectric material characteristics are given in table 1.

### 3.1 Optimization of the waves group velocities by using $\mathrm{Crit}_1$

In this part, optimization of the transmission capability of the designed adaptive metacomposite is considered by using $\mathrm{Crit}_1$ given in equation (23). The objective is to avoid any energy transportation when flexural waves ($A_0, A_1, TH...$) are excited into the periodically distributed shunted piezocomposite cells. The numerical optimization procedure is based on a multidimensional unconstrained nonlinear minimization algorithm (Nelder-Mead). The optimization is done by considering an active/reactive electronic circuit through an impedance $Z(\omega) \in \mathbb{C}$.

The analysis is initialized with an arbitrary complex value of the shunt impedance. Optimization steps are then proceeded using criterion (23) by considering any frequency dependent complex impedance for describing the circuit behaviour. Figure 2 shows the dispersion curves using the propagative parts (here the real one) of the wave numbers $kx_n(\omega)$ along (Ox) axis, for both initial and optimal impedance shunts.
<table><caption>Table 1: *Piezoelectric patch characteristics*</caption><thead><tr><th colspan="3">Piezoelectric Material</th></tr><tr><th>Symbol</th><th>Value</th><th>Property</th></tr></thead><tbody><tr><td>s<sub>11</sub><sup>E</sup> = s<sub>22</sub><sup>E</sup> = s<sub>33</sub><sup>E</sup></td><td>11.6e<sup>-12</sup> Pa<sup>-1</sup></td><td></td></tr><tr><td>s<sub>12</sub><sup>E</sup> = s<sub>13</sub><sup>E</sup> = s<sub>23</sub><sup>E</sup></td><td>-3.33e<sup>-12</sup> Pa<sup>-1</sup></td><td>Compliance matrix</td></tr><tr><td>s<sub>44</sub><sup>E</sup> = s<sub>55</sub><sup>E</sup> = s<sub>66</sub><sup>E</sup></td><td>45.0e<sup>-11</sup> Pa<sup>-1</sup></td><td></td></tr><tr><td>η</td><td>0.1 %</td><td>Hysteretic Damping ratio</td></tr><tr><td>d<sub>31</sub> = d<sub>32</sub></td><td>-6e<sup>-11</sup> C/N</td><td></td></tr><tr><td>d<sub>33</sub></td><td>15.2e<sup>-11</sup> C/N</td><td>Piezoelectric matrix</td></tr><tr><td>d<sub>24</sub> = d<sub>15</sub></td><td>730e<sup>-12</sup> C/N</td><td></td></tr><tr><td>ρ</td><td>7600kg/m<sup>3</sup></td><td>Density</td></tr><tr><td>ε<sub>11</sub><sup>T</sup> = ε<sub>22</sub><sup>T</sup></td><td>504.1 ε<sub>o</sub> C/V/m</td><td>Dielectric Permittivity</td></tr><tr><td>ε<sub>33</sub><sup>T</sup></td><td>270 ε<sub>o</sub> C/V/m</td><td></td></tr></tbody></table>
<figure><img src="image_3.png" /><figcaption>Figure 2: Propagative parts of the wave number $kx_n(i\omega)$ along ($Ox$) (all waves included). o: initial value of shunt; +: optimal value of shunt.</figcaption></figure>

This figure underlines the ability of the procedure to catch all Bloch's solutions including evanescent ones that generally do not appear in literature because of their imaginary character. In the context of structural dynamics, including damping effects is mandatory for real-life applications and all wave numbers become complex, which constitutes a key point in the analysis as underlined in (Collet, Ouisse, Ruzzene & Ichchou 2011, Garvic 1995).

By using a suitable wave filter based on the group velocity, it is possible to retain for presentation only propagative solutions responsible for energy transport. For a sake of clarity, the following results will be presented by exhibiting only propagative waves. Above results will be filtered and presented as in figure 3, which shows the propagative parts of the wave number $kx_n(\omega)$ along ($Ox$). It then becomes easy to observe branches similar to standard $S_0$, $A_0$, $A_1$, TH waves.

The corresponding group velocities along ($Ox$) are presented in figure 4.

It can immediately be observed that the optimization of the shunt impedance leads to a large decrease of the group velocity of the $A_0$ mode while the $A_1$ wave, which becomes propagative at 8.8 kHz, is not controlled by the optimal configuration (see figure 4). The bending waves also propagate energy with a very slow velocity and can be considered as evanescent. Flexural energy is, also, only transported by the $A_1$ mode after the cutting frequency.

The propagative part of the controlled wave is strongly modified, as illustrated in figure 3: from two complex conjugated solutions, the $A_0$ wave changes after control to 4 solutions (2 complex conjugates and their opposites). This situation is described in (Garvic 1995), it corresponds to high order evanescent waves solutions.

It can also be observed on the group velocity diagram 4 that a non dispersive symmetric wave is also affected by the control whereas it was not targeted by the criterion.

The optimization has been performed for several values of $\phi$. In order to show the stability of the optimized configuration with this angle, the figure 5 shows the propagative parts of the wave number $kx_n(\omega)$ along 6 directions forming angles with the ($Ox$) axis of $\frac{n\pi}{20}$, $n$ being all integers from 0 to 5. It is clear on this figure that above comments are valid for all angles. In particular, the
<figure><img src="image_4.png" /><figcaption>Figure 3: Propagative parts of the wave number $kx_n(\omega)$ along ($Ox$) (only propagating waves included). o: initial value of shunt; +: optimal value of shunt. Large marks correspond to flexural waves, small marks correspond to other waves</figcaption></figure>

high order of evanescence of controlled waves in the optimized configurations can be observed in all situations, while some non dispersive symmetric waves are systematically affected by the control whereas they were not targeted by the criterion.

The real and imaginary parts of the optimal impedance are plotted in figure 6 for all angles $\phi$. The optimal impedance values almost correspond to a constant negative capacitance in all directions. The corresponding average value is $-233.66 \, pC.V^{-1}$. Equivalent resistances corresponding to the active part of the shunt impedance (figure 6) are negative which indicates that the optimization leads to provide energy to the system for controlling mechanical damping effects introduced by hysteretic damping in the model. The optimized configuration tends to converge towards a fully conservative system. The obtained mean value of the resistance is $-1.5319 \, \Omega$. This is confirmed in figure 7 where the electrical dissipated energy appears negative when the optimal shunt is connected to the patch.

In order to check this point, an optimization run has been performed on a highly damped system, using an hysteretic coefficient of 5 %. Figure 8 show the corresponding dispersion curves for $\phi = 0$. The optimized configurations leads to a behaviour which is very similar to the one illustrated in figure 3 for lower structural damping value. The same observation can be done for other angle values. This similar behaviour is only related to propagative parts of the wave numbers, while the dissipative effects are modified as illustrated in figure 9 which shows the new optimal impedances. It can be underlined that the electric circuit is fully active with larger resistive parts of average value $-74.88 \, \Omega$. The observed proportionality between the optimal average resistances ($-1.5319$ and $-74.88 \, \Omega$) and the chosen mechanical damping terms (0.1 % and 5 %) indicates that the optimization leads to annihilate mechanical damping in order to reach a total cancellation of the group velocity of $A_0$ mode. The final controlled system is then almost fully reactive.

This last point is of first importance in both physical and technological points of view. Indeed, this results tends to show, on this example, that a pure band gap effect is only reachable if the
<figure><img src="image_5.png" /><figcaption>Figure 4: Group velocities along ($Ox$) direction. o: initial value of shunt; +: optimal value of shunt.</figcaption></figure>

Large marks correspond to flexural waves, small marks correspond to other waves

system is completely reactive, which means that no damping occurs in the system. The technological implementation of such interface allowing cancellation of structural transmissibility is then unrealistic since it leads to the implementation of only 'Lyapunov stable' system: the stability margins have to be drastically minimized by using negative resistance for obtaining transmission control.

## 3.2 Optimization of damped power flow inside the electric shunts by using $Crit_2$

Another strategy for optimizing the adaptive metacomposite consists in focusing on the damped power flow inside the electric shunts by using the criteria $Crit_2$ given in equation (25). In this case the objective is to improve the absorption capability of the smart structure.

The propagative parts of the wave number $kx_n(i\omega)$ along ($Ox$) axis appear as in figure 10. The corresponding group velocities along ($Ox$) are presented in figure 11, while the real and imaginary parts of the optimal impedance are plotted in figure 13 and the electrically damped power in figure 12.

The optimal impedances of the electric shunts are plotted in figure 13 for various values of $\phi$.

The first observation is that the optimization of the shunt impedance for improving the absorption characteristics of the system induces modifications of the group velocities of the controlled waves (see figure 11), while the propagative part of the wave numbers remain unchanged (figure 10). This can be explained by a large improvement of the ratio between the real and imaginary parts of the waves numbers, which physically corresponds to the forcing of propagating effects to increase damping effects: energy can propagate inside the periodically distributed set of active cells for allowing electrical energy conversion. The figure 12 shows that a the dissipated power is largely increased when optimal shunt is connected to the patch.

The optimal impedance values correspond to an almost constant negative capacitance in all directions (see figure 13). The corresponding average value is $-237.43 \ pC.V^{-1}$. Equivalent resis-
<figure><img src="image_6.png" /><figcaption>Figure 5: Propagative parts of the wave number $kx_n(i\omega)$ along 6 directions forming an angle with $(Ox)$ axis of $\frac{n\pi}{20}$ with $n$ being all integers from 0 to 5. o: initial value of shunt; +: optimal value of shunt. Large marks correspond to flexural waves, small marks correspond to other waves</figcaption></figure>

tances corresponding to the active part of the shunt impedance (figure 13) are positive, which is in accordance with the fact that a damping effect is awaited.

## 3.3 Validation on a periodically semi-distributed set of adaptive cells

Up to that point, the optimizations of the smart metacomposite which have been proposed were performed by considering infinite periodic structures. The objective of this section is to illustrate the applicability of the designs for integration into real-life finite structures. The optimal impedance is then applied on a finite set of shunted piezo-composite cells distributed on a part of a totally free plate system submitted to a point force in one corner, as described in figure 14. Material damping is introduced throughout a constant imaginary part of the elasticity tensor of 0.1%. The harmonic response of this system is then computed on the [500, 10000] Hz frequency band when optimal impedances, obtained by criteria $Crit_1$ or $Crit_2$, are connected or not to each patches.

The figure 15 presents the mean kinetic energy on the considered system (i.e on the whole domain $\Omega = \Omega_s \cup \Omega_i \cup \Omega_t$) with open circuit and $Crit_1$, $Crit_2$ optimal shunts connected to each piezoelectric patch.
<figure><img src="image_7.png" /><figcaption>Figure 6: Optimal electric impedance, represented as equivalent resistance and capacitance, obtained along 6 direction forming an angle with (Ox) axis of $\frac{n\pi}{20}$ with n being all integers from 0 to 5.</figcaption></figure>

From the observation of figure 15, the main particularities of the synthesized optimal impedances can be emphasized. The first configuration, obtained by using $Crit_1$, lead the adaptive interface to avoid any flexural energy transfer. The mechanical response of the system is strongly modified, and the resonance peaks are moved far form the open circuit response, which indicates a modification of the internal dynamics. The second configuration, based on $Crit_2$ optimal shunt, forces energy dissipation into the resistive part of the electric shunt. In that case, the conservation of the modal architecture into the energy spectrum can be observed, even if some damping effects limit the peak's amplitudes.

To go further, one can take a closer look at the energy repartition into the considered system in figure 16, which shows the ratio of the kinetic energy finally transmitted to the sub-domain $\Omega_t$. The minimization of the flexural wave group velocity ($Crit_1$) leads to decrease the part of kinetic energy transmitted to the non-excited part of the plate located behind the adaptive interface while the maximization of the damped energy ($Crit_2$) only slightly modify this distribution. The modification of the energy distribution into the system is very important after 6kHz and increase with the frequency when the energy is globally transported by $A_0$ mode. This is not the case in the lower frequency band when the interface works in the near field domain of the applied point force.

The figure 17 shows the corresponding results in terms of velocity root mean square.

These numerical results clearly show a strong influence of the modifications of the wave dispersion characteristics on the standing wave responses plotted in figure 17 especially when $Crit_1$ is used. In particular, at 2550Hz, a large attenuation of the dynamic response of the system can be observed when optimal shunts are connected. In this case, when the shunts are optimized for cancelling flexural energy flow ($Crit_1$), the location of the nodal line are also modified. With an optimization aiming at damping energy in the electric resistance ($Crit_2$), only to a standard attenuation can be observed. At 8000 Hz when $Crit_1$ optimal shunts are connected, a large part of the system energy remains in its left part where the applied force is located. The smart metacomposite interface also changes the system admittance and filter wave diffusion by increasing is reflexivity properties. The energy diffusion is also clearly condensed into the left part of the system with a largely decreased
<figure><img src="image_8.png" /><figcaption>Figure 7: Damped power inside the electric circuit. o: initial value of shunt; +: optimal value of shunt. Large marks correspond to flexural waves, small marks correspond to other waves</figcaption></figure>

amplitude compared to the one obtained with open circuit and $Crit_2$ optimal shunt. With this type of shunt something similar to a wave trap effect can be observed. With $Crit_2$ optimal shunt, an improvement of the damping effect leading to vibration attenuation can be observed, without wave trap behaviour.

## 4 Conclusions

This article presents a numerical procedure able to compute the damped wave's dispersion functions in the whole first Brillouin domain of multi dimensional piezo-elastodynamical wave guides. The method was applied for determining the optimal shunt impedance allowing the minimization of the group velocities of the flexural waves and the damping capability of the distributed system. The smart metacomposite can then be included in a finite structure to confer on it new functionalities. Some numerical tests on a finite dimension system incorporating a semi-distributed set of shunted piezo-composite cells has been performed. A strong influence of the designed shunt circuits has been underlined in the dynamical response of the system and energy distribution into each structural sub domains. Even if the link between the obtained wave properties are not clearly established especially in the low frequency range where the adaptive interface is located in the near field domain of the input force, it has been demonstrated that the numerical procedures can be used for optimizing the energy diffusion operator of such adaptive mechanical interface. To do so, additional work has to be done for optimizing the complete interface scattering and for controlling the evanescent waves playing an important role in the finite system dynamical response in the low frequency band.

The proposed methodology can also be used for studying particular dissipation phenomenon such as those induced by complex shunted piezoelectric patches as proposed by (Beck, Cunefare & Ruzzene 2008b) and (Casadei et al. 2009), or even foams or complex polymers behaviours. The
<figure><img src="image_9.png" /></figure>

Figure 8: Propagative parts of the wave number $kx_n(i\omega)$ along ($Ox$), for high hysteretic damping of 5% (only propagating waves included). o: initial value of shunt; +: optimal value of shunt. Large marks correspond to flexural waves, small marks correspond to other waves

proposed method furnishes an efficient tool for future optimization of distributed smart cells as proposed in the case of 1D wave guide by (Collet et al. 2009).

## References

### References

* Aalami, B. (1973). Waves in prismatic guides of arbitrary cross section, *Journal of Applied Mechanics* **40**(4): 1067–1072.

* Agnes, G. (1994). Active/passive piezoelectric vibration suppression, in C. D. Johnson (ed.), *Proceedings of SPIE Smart Structures and Materials*, pp. 24–34.

* Allaire, G. & Congas, C. (1998). Bloch waves homogenization and spectral asymptotic analysis, *Journal de Mathématiques Pures et Appliquées* **77**: 153–208.

* Banks, H., Smith, R. & Wang, Y. (1996). *Smart material structures: Modeling Estimation and Control*, Masson and Wiley.

* Bartoli, I., Marzani, A. and Lanza Di Scalea, F. & Viola, E. (2006). Modeling wave propagation in damped waveguides of arbitrary cross-section, *Journal of Sound and Vibration* **295**: 685–707.
<figure><img src="image_10.png" /></figure>

Figure 9: Optimal electric impedance for high hysteretic damping of 5 %, represented as equivalent resistances and capacitances, obtained along 6 direction forming an angle with (Ox) axis of $\frac{n\pi}{20}$ with n being all integers from 0 to 5.

Baz, A. (1997). Optimization of energy dissipation characteristics of active constrained layer damping, *Smart Materials and Structures* **6**(3): 360.

Beck, B., Cunefare, K. & Ruzzene, M. (2008a). Broadband vibration suppression assessment of negative impedance shunts, *Proceedings of ASME - SMASIS*.

Beck, B., Cunefare, K. & Ruzzene, M. (2008b). Broadband vibration suppression assessment of negative impedance shunts, *Proceedings of ASME - SMASIS*.

Becker, J., Fein, O., Maess, M. & Gaul, L. (2006). Finite element-based analysis of shunted piezoelectric structures for vibration damping, *Computers and Structures* **84**(31-32): 2340-2350.

Behrens, S., Fleming, A. & Moheimani, S. (2001). New method for multiple-mode shunt damping of structural vibration using a single piezoelectric transducer, in D. Inman (ed.), *Proceedings of SPIE Smart Structures and Materials*, pp. 239-250.

Bensoussan, A., Lions, J. & Pananicolaou, G. (1978). *Asymptotic Analysis for Periodic Structures*, North Holland.

Bloch, F. (1928). Über die Quantenmechanik der Electron in Kristallgittern, *Zeitschrift für Physik* **52**: 550-600.

Bocquillet, A., Ichchou, M. N. & Jezequel, L. (2003). Energetics of axisymmetric fluid-filled pipes up to high frequencies, *Journal of Fluids and Structures* **17**: 491-510.

Casadei, F., Beck, B., Ruzzene, M. & Cunefare, K. (2009). Vibration control of plates featuring periodic arrays of hybrid shunted piezoelectric patches, *Proceedings of SPIE Smart Structures and Materials*.
<figure><img src="image_11.png" /></figure>

Figure 10: Propagative parts of the wave number $kx_n(i\omega)$ along ($Ox$), for high hysteretic damping of 5% (only propagating waves included). o: initial value of shunt; +: optimal value of shunt. Large marks correspond to flexural waves, small marks correspond to other waves

Chen, F. & Wilcox, P. (2007). The effect of load on guided wave propagation, *Ultrasonics* **47**(1-4): 111–122.

Collet, M. & Cunefare, K. (2008). Modal Synthesis and Dynamical Condensation Methods for Ac curate Piezoelectric Systems Impedance Computation, *Journal of Intelligent Material Systems and Structures* **19**(11): 1251–1271.

Collet, M., Cunefare, K. & Ichchou, M. (2009). Wave Motion Optimization in Periodically Dis tributed Shunted Piezocomposite Beam Structures, *Journal of Intelligent Material Systems and Structures* **20**(7): 787–808.

Collet, M., Ouisse, M., Ichchou, M. & Ruzzene, M. (2010). Numerical tools for semi-active opti mization of 2d wave’s dispersion into mechanical system, *Proceedings of ASME-SMASIS*.

Collet, M., Ouisse, M., Ichchou, M. & Ruzzene, M. (2011). Semi-active optimization of 2d wave’s dispersion into mechanical systems by the mean of periodically distributed shunted piezoelectric patches: a new class of adaptive metamaterials, *Proceedings of SPIE Smart Structures/NDE*.

Collet, M., Ouisse, M., Ruzzene, M. & Ichchou, M. (2011). Floquet-bloch decomposition for the com putation of dispersion of two-dimensional periodic, damped mechanical systems, *International Journal of Solids and Structures* **48**(20): 2837–2848.

Corr, L. & Clark, W. (2002). Comparison of lowfrequency piezoelectric switching shunt techniques for structural damping, *Journal of Smart Material and Structures* **11**(3): 370–376.

Cross, C. & Fleeter, S. (2002). Shunted piezoelectrics for passive control of turbomachine blading flow-induced vibrations, *Smart Material Structures* **11**: 239–248.
<figure><img src="image_12.png" /></figure>

Figure 11: Group velocities along (Ox) direction. o: initial value of shunt; +: optimal value of shunt. Large marks correspond to flexural waves, small marks correspond to other waves

Cunefare, K. (2006). Negative capacitance shunts for vibration suppression: wave based tuning and reactive input power, *Proceedings of Active 2006*.

Damljanovic, V. & Weaver, R. L. (2004). Propagating and evanescent elastic waves in cylindrical waveguides of arbitrary cross section, *Journal of the Acoustical Society of America* **115**: 1572–1581.

dell'Isola, F., Maurini, C. & Porfiri, M. (2004). Passive damping of beam vibrations through distributed electric networks and piezoelectric transducers: prototype design and experimental validation, *Smart Material Structures* **13**: 299–308.

Den Hartog, J. (1956). *Mechanical vibrations*, 4th edn, McGraw-Hill, New-York.

Duhamel, D., Mace, B. & Brennan, M. (2006). Finite element analysis of the vibrations of waveguides and periodic structures, *Journal of Sound and Vibration* **294**(1-2): 205–220.

Duhamel, E. (2007). Finite element computation of green's functions, *Engineering Analysis with Boundary Elements* **31**: 919–930.

Finnveden, S. (2004). Evaluation of modal density and group velocity by a finite element method, *Journal of Sound and Vibration* **273**: 51–75.

Finnveden, S. & Fraggstedt, M. (2008). Waveguide finite elements for curved structures, *Journal of Sound and Vibration* **312**(4-5): 644–671.

Fleming, A., Behrens, S. & Moheimani, S. (2003). Reducing the inductance requirements of piezoelectric shunt damping systems, *Smart Materials and Structures* **12**(1): 57–64.

Fleming, A., Belirens, S. & Moheimani, S. (2000). Synthetic impedance for implementation of piezoelectric shunt-damping circuits, *Electronics Letters* **36**(18): 1525–1526.
<figure><img src="image_13.png" /><figcaption>Figure 12: Damped power inside the electric circuit. o: initial value of shunt; +: optimal value of shunt. Large marks correspond to flexural waves, small marks correspond to other waves</figcaption></figure>

Fleming, A. J. & Moheimani, S. O. R. (2004). Improved Current and Charge Amplifiers for Driving Piezoelectric Loads, and Issues in Signal Processing Design for Synthesis of Shunt Damping Circuits, *Journal of Intelligent Material Systems and Structures* **15**(2): 77–92.

Floquet, G. (1883). Sur les équations différentielles linéaires à coefficients périodiques, *Annales de l'Ecole Normale Supérieure* **12**: 47-88.

Forward, R. (1978). Electromechanical transducer-coupled mechanical structure with negative capacitance compensation circuit, Patent US4158787.

Fukada, E., Date, M. & Kimura, K. (2002). Sound isolation by piezoelectric polymer films connected with negative capacitance circuits, *Proceedings of 11th International Symposium on Electrets*.

Garvic, L. (1995). Computation of propagative waves in free rail using a finite element technique, *Journal of Sound and Vibration* **185**: 531–543.

Garvic, L. (1995). Computation of propagatives waves in free rail using a finite element technique, *Journal of Sound and Vibration* **185**(3): 531–543.

Hagood, N. & von Flotow, A. (1991). Damping of structural vibrations with piezoelectric materials and passive electrical networks, *Journal of Sound and Vibration* **146**(2): 243–268.

Hardy, P., Ichchou, M., Jezequel, L. & Trentin, D. (2009). A hybrid local energy formulation for plates mid-frequency flexural vibrations, *European Journal of Mechanics - A/Solids* **28**(1): 121-130.

Hayashi, T., Song, W. & Rose, J. (2003). Guided wave dispersion curves for a bar with an arbitrary cross-section, a rod and rail example, *Ultrasonics* **41**(3): 175-183.
<figure><img src="image_14.png" /><figcaption>Figure 13: Optimal electric impedance, represented as equivalent resistances and capacitances, obtained along 6 direction forming an angle with (Ox) axis of $\frac{n\pi}{20}$ with n being all integers from 0 to 5.</figcaption></figure>

Houillon, L., Ichchou, M. & Jezequel, L. (2005). Wave motion in thin-walled structures, *Journal of Sound and Vibration* **281**(3-5): 483–507.

Hurlebaus, S. & Gaul, L. (2006). Smart structure dynamics, *Mechanical System and Signal Processing* **20**: 255–281.

Ichchou, M., Akrout, S. & Mencik, J. (2007). Guided waves group and energy velocities via finite elements, *Journal of Sound and Vibration* **305**(4-5): 931–944.

Joannopoulos, J., Meade, R. & Winn, J. (1995). *Photonic Crystals: Molding the Flow of Light*, Princeton University Press.

Kittel, C. (1986). *Introduction to Solid State Physics*, John Wiley and Sons, New York.

Mace, B., Duhamel, D., Brennan, M. & Hinke, L. (2005). Finite element prediction of wave motion in structural waveguides, *Journal of Acoustical Society of America* **117**(5): 2835–2843.

Maess, M., Wagner, N. & Gaul, L. (2006). Dispersion curves of fluid filled elastic pipes by standard fe models and eigenpath analysis, *Journal of Sound and Vibration* **296**(1-2): 264 – 276.

Maysenhölder, W. (1994). *Körperschall-energie Grundlagen zur Berechnung von Energiedichten und Intensitäten*, Wissenschaftliche Verlagsgesellschaft, Stuttgart.

Mead, D. (1973). A general theory of harmonic wave propagation in linear periodic systems with multiple coupling, *Journal of Sound and Vibration* **27**: 235–260.

Mencik, J. & Ichchou, J. (2005a). Multi-mode propagation and diffusion in structures through finite elements, *European Journal of Mechanics A-Solids* **24**(5): 877–898.
<figure><img src="image_15.png" /><figcaption>Figure 14: 3D piezocomposite periodic semi-distributed cells as plate interface</figcaption></figure>

Mencik, J. & Ichchou, M. (2005b). Multi-mode propagation and diffusion in structures through finite elements, *European Journal of Mechanics - A/Solids* **24**(5): 877–898.

Mencik, J. & Ichchou, M. (2007). Wave finite elements in guided elastodynamics with internal fluid, *International Journal of Solids and Structures* **44**: 2148–2167.

Milman, M. H. & Cheng-Chih, C. (1994). Optimization methods for passive damper placement and tuning, *J. of guid. cont. and dyn.* **17**(4): 848–856.

Moheimani, S. & Fleming, A. (2006). *Piezoelectric transducers for vibration control & damping, (Advances in industrial control)*, Springer.

Nelson, P. & Elliott, S. (1992). *Active Control of Sound*, Pub. Academic Press, London, San Diego.

Nelson, R. & Dong, S. (1973). High frequency vibrations and waves in laminated orthotropic plates, *Journal of Sound and Vibration* **30**(1): 33–44.

Neubauer, M., Oleskiewicz, R., Popp, K. & Krzyzynski, T. (2006). Optimization of damping and absorbing performance of shunted piezo elements utilizing negative capacitance, *Journal of Sound and Vibration* **298**(1-2): 84–107.

Ormondroyd, J. & Den Hartog, J. (1928). The theory of dynamic vibration absorber, *Trans., ASME APM-50-7*: 9–22.

Park, C. & Baz, A. (2005). Vibration control of beams with negative capacitive shunting of inter-digital electrode piezoceramics, *Journal of Vibration and Control* **11**: 331–346.

Park, J. & Palumbo, D. (2004). A new approach to identify optimal of shunting elements for maximum damping of structural vibration using piezoelectrical patches, *Active 2004*, Institute of Noise Control Engineers.

Preumont, A. (2002). *Vibration control of active structures: an introduction*, Springer.

Spencer, B. F. & Nagarajaiah, S. (2003). State of the art of structural control, *J. Struct. Eng.* **129**(7): 1–10.
<figure><img src="image_16.png" /></figure>

Figure 15: Mean value of the kinetic energy on the whole domain Ω with open circuit (red dashed line), *Crit*₁ (blue line) and *Crit*₂ (black marked line) optimal shunt connected to piezoelectric patches

Thorp, O., Ruzzene, M. & Baz, A. (2001). Attenuation and localization of wave propagation in rods with periodic shunted piezoelectric patches, *Smart Materials and Structures* **10**: 979–989.

Thorp, O., Ruzzene, M. & Baz, A. (2005). Attenuation of wave propagation in fluid-loaded shells with periodic shunted piezoelectric rings, *Smart Materials and Structures* **14**: 594–604.

Wu, S. (1998). Method for multiple-mode shunt damping of structural vibration using a single pzt transducer, *Proceedings of SPIE Smart Structures and Materials*.

Wu, S. & Bicos, A. (1997). Structural vibration damping experiments using improved piezoelectric shunts, in L. P. Davis (ed.), *Proceedings of SPIE Smart Structures and Materials*, pp. 40–50.

Zhou, W., Ichchou, M. & Mencik, J. (2009). Analysis of wave propagation in cylindrical pipes with local inhomogeneities, *Journal of Sound and Vibration* **319**(1-2): 335–354.
<figure><img src="image_17.png" /><figcaption>Figure 16: Ωt ratio of the total kinetic energy mean value in percent with open circuit (red dashed line), Crit1 (blue line) and Crit2 (black marked line) optimal shunt connected to piezoelectric patches</figcaption></figure>
<figure><img src="image_18.png" /><figcaption>Figure 17: Velocity root mean square level of the system submitted to an external point force at 2550 Hz (left) and 8000 Hz (right) for all piezoelectric cells in open circuit (a), connected to optimal shunt impedance obtained with $Crit_1$ (b) and $Crit_2$ (c)</figcaption></figure>