Smart Materials Research
Volume 2015 (2015), Article ID 748459, 15 pages
http://dx.doi.org/10.1155/2015/748459

Research Article

Thermal Effects on Vibration and Control of Piezocomposite Kirchhoff Plate Modeled by Finite Elements Method

M. Sanbi,¹ R. Saadani,² K. Sbai,² and M. Rahmoune²

¹Equipe Sciences et Technologies Avancées, Ecole Nationale des Sciences Appliquées, Université Abdelmalek Essaadi, 93030 Tétouan, Morocco
²Laboratoire d’Etude des Matériaux Avancés et Applications, Faculté des Sciences et Ecole Supérieure de Technologie, Université Moulay Ismail, 50040 Meknès, Morocco

Received 5 December 2014; Revised 2 April 2015; Accepted 3 April 2015

Academic Editor: Weihua Li

Copyright © 2015 M. Sanbi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Theoretical and numerical results of the modeling of a smart plate are presented for optimal active vibration control. The smart plate consists of a rectangular aluminum piezocomposite plate modeled in cantilever configuration with surface bonded thermopiezoelectric patches. The patches are symmetrically bonded on top and bottom surfaces. A generic thermopiezoelastic theory for piezocomposite plate is derived, using linear thermopiezoelastic theory and Kirchhoff assumptions. Finite element equations for the thermopiezoelastic medium are obtained by using the linear constitutive equations in Hamilton’s principle together with the finite element approximations. The structure is modelled analytically and then numerically and the results of simulations are presented in order to visualize the states of their dynamics and the state of control. The optimal control LQG-Kalman filter is applied. By using this model, the study first gives the influences of the actuator/sensor pair placement and size on the response of the smart plate. Second, the effects of thermoelastic and pyroelectric couplings on the dynamics of the structure and on the control procedure are studied and discussed. It is shown that the effectiveness of the control is not affected by the applied thermal gradient and can be applied with or without this gradient at any time of plate vibrations.

1. Introduction

In the piezoelectric sensors applications, mechanically or thermally induced deformations can be determined from measurement of the induced electrical potential, whereas in piezoelectric actuator applications deformation of strain can be controlled through the introduction of appropriate electric potential. By integrating distributed piezoelectric sensors/actuators and advanced composites, the potential exists for forming high-strength, high stiffness, lightweight structures capable of self-monitoring and self-controlling. Typical applications of such structures are envisioned in the thermal distortion management of propulsion components and space structures. Before they can be utilized in these applications, the performance of piezoelectric structures in thermal environment must be quantified. Thermal effects become important when the piezoelectric structure has to operate in either extremely hot or cold temperature environments. These extreme conditions may severely affect the response of piezoelectric elements by induction of thermal stresses resulting from thermoelastic and pyroelectric coefficients.

Recently, there have been various mathematical models developed to describe the behavior of the piezocomposite or laminated plates actuated and sensed by piezoelectric materials. Görnandt and Gabbert [1] have analysed the finite element method of thermopiezoelectric smart structures. Vel and Batra [2] presented a generalized plane strain thermopiezoelectric analysis of multilayered plates. de Abreu et al. [3] have implemented finite element modeling of a plate with localized piezoelectric sensors and actuators. Deü and Benjeddou [4] have studied free-vibration analysis of laminated plates with embedded shear-mode piezoceramic layers. Tanaka and Sanada [5] have detailed modal control of a rectangular plate using smart sensors and smart actuators. Trindade and Benjeddou [6] have discussed the critical evaluation and optimization of the effective electromechanical coupling coefficients of piezoelectric adaptive structures. Sanbi et al. [7] have analyzed thermoelastic and pyroelectric couplings effects on dynamics and active control of piezolaminated beam. A general solution for piezothermoelasticity of transversely isotropic piezoelectric materials and its applications have been deduced by Ding et al. [8].

In this work, fundamental equations governing the behavior of smart structures are at first given, and, based on the work of de Abreu et al. [3], the generalized linear finite element formulation of the problem is described. Thermopiezoelastic characteristics of aluminium rectangular plate and the theory of active structures in sensing and optimal control are established. A generic thermopiezoelastic theory for piezocomposite plate is derived, using linear thermopiezoelasticity and Kirchhoff assumptions. Generalized finite element equations for the thermopiezoelasticity are obtained by using the linear constitutive equations in Hamilton’s principle together with the finite element approximations. A Kirchhoff four-node rectangular element with one electrical, one temperature, and three mechanical degrees of freedom is considered. The structure consists of a modeling of cantilevered piezocomposite plate with perfectly surface bonded thermopiezoelectric elements. The structure is modelled analytically and then numerically and the results of simulations are presented in order to visualize the states of their dynamics and the state of control. The optimal control LQG-Kalman accompanied is applied. The effects of thermoelastic and pyroelectric couplings coefficients on the dynamics of the structure and on the control effectiveness are discussed. We show that the control procedure cannot be perturbed by applying a thermal gradient and the control can be applied at any time during the period of vibration of the plate.

2. Basic Equations and FEM Implementation

The physical and dimensional characteristics of the plate material are given in Table 1. Properties for the two piezoelectric materials used in the simulation can be found in the literature as in [9, 10]. The piezoelectric and elastic properties of the sensor are chosen smaller compared to those of the actuator. This is justified by the fact that the sensor must have high sensitivity to capture the deformation of the structure, while the actuator is required to act as a secondary source on the plate to deaden its vibrations. In order to model the structure by finite element method, we consider an element with four-node quadrilateral flexible element type “plate.” Each node has three mechanical degrees of freedom (dofs), the displacement in -direction and two rotations in - and -direction, one electric dof , and one temperature dof (Figure 1).

Table 1: Properties of the plate, sensor, and actuator.

Properties/labels/units
Plate
Sensor
Actuator

Al
PVDF
PZT-5H32

Length (m)
0.6
0.25
0.25

Width (m)
0.4
0.25
0.25

Thickness (m)
1
0.2
0.2

Density
8.03
1.78
7.7

Young modulus (GPa)
68
5.04
62

Reduced elastic constants (GPa)

—
5.04
66.2

—
3.25
19.2

—
—
23.5

Piezoelectric stress constant
—
0.15
−9.11

(Vm/N)

Piezoelectric strain constant
—
4.34
−274

  ×10⁻¹² (m/V)

Reduced piezoelectric constant
—
0.018
−23.4

(C/m²)

Pyroelectric constant
—
−1.25
—

  ×10⁻⁵ Cm⁻²K⁻¹

Dielectric constant
—
12
3300

Thermal conductivity
—
1.5
1.5

(W/m·K)

Thermal expansion
—
140
3

  ×10⁻⁶ (K⁻¹)

Specific heat J·kg⁻¹·K⁻¹
—
350
1.5

Damping   ×10⁻³
1
—
—

Constants   ×10⁻³
0.1
—
—

Figure 1: The finite element type “plate” used in the formulation and different degrees of freedom.

In the present formulation, we assume that the piezoelectric layers are perfectly bonded and the formulation is restricted to linear elastic material behavior with small displacement and strains. This formulation uses the Kirchhoff assumption for thin plate in which the transverse normal remains straight and normal to the midsurface after deformation. The thickness of the plate does not change during a deformation. Depending on these assumptions and as shown in Figure 1, displacements , , and can be expressed by

where and are the longitudinal and transverse axes located on the medial surface of the plate and is the thickness direction (Figure 1). The displacements , , and are taken along the axes , , and , respectively. Assuming that the transverse deformation is negligible, the deformation in terms of displacements can be written as

For an isotropic material, the constraints are related to deformation by

where and the matrix is given by

with being the Poisson coefficient and the plate Young modulus.

In the finite element modeling, a four-node bending plate rectangular element based on the classical theory of plates is considered [11]. Each node element has three mechanical degrees of freedom: the displacement in -direction, -rotation , and -rotation (Figure 1). The function of the displacement can be written as

with

 ; ; ; ; ; ; ; ; .

The transverse displacement can be expressed by

where

We define the vector that contains the nodal displacements of the rectangular element by

where

By combining (5) and (8) with (9) for the four nodes, we obtain the following matrix expression:

where is a matrix given by

Therefore, the coefficient vector can be calculated from (10) as follows:

Substituting (12) into (6) leads to

where is the -direction shape function matrix given by

Substituting (3) into (2) yields

By manipulating (15), we can obtain

with

The displacements , , and can be expressed by the vector as

The substitution of (1) into (18) yields

and substituting (13) into (19), we get

where and are given by

Considering thermopiezoelastic material, the constitutive equations of linear thermopiezoelasticity and temperature are written as [12]

where the exponent means that the values are measured at constant strain and the exponent means they are measured at a constant electric field, is the stress tensor, is the vector of electric displacement, is the entropy, is the temperature vector, is the tensor deformation, is the electric field, is the matrix of elastic coefficients at constant electric field, is the matrix of piezoelectric coefficients, is the tensor dielectric constant mechanical strain, is the pyroelectric tensor, and is the thermoelastic tensor. is the expansion coefficient given by , where and are the specific heat and the initial temperature, respectively.

On the other hand, for the relation between and , the coefficients of the piezoelectric strain are written as

If the polarization direction is held in, the electrical potential applied or captured through the elements sensor or actuator is given by the following equation [13]:

where and are the thickness and the maximum electric potential at the outer surface of the corresponding thermopiezoelectric element (actuator or sensor) and ( or ) is defined in the following two intervals:

Assuming that the electric field is constant along the thickness of the actuator and sensor, we can therefore write

Hamilton’s principle is used here to derive finite element equations and is given by

where and are two arbitrary instants, is the kinematic energy, is the energy potential, is the work due to electrical forces, and is the work due to thermal stresses. The total kinematic energy and energy potential are given by the following relations:

with being the derivative with respect to time and a volume element given by

where the indexes , , and refer to the plate, the actuator, and the sensor, respectively, and , , and are given by

The energy of electric forces and elementary virtual work produced by the volume forces , surface forces , and the applied electrical loads are written as

with being the electric displacement vector, the volume force, the surface force, and the electric surface stress.

Substituting (3) and (22) into (31) and (23) into (34), respectively, leads to

The substitution of (30), (35), and (36) into (29) gives

Substituting (20), (16), and (28) into (37) leads to

where

Considering that the variation of and is arbitrary, we obtain two equilibrium equations rewritten in generalized coordinates for the th element:

with being the extended stiffness matrix and the elementary mass matrix. By integrating in the direction of , we obtain

where are given by

and , , and , for , are calculated from (4) for the thermopiezoelectric material and the plate, respectively, and .

The integration in of the elementary mass matrix gives

with , , and and (for ) are

By integrating the direction , the thermoelastic coupling matrix for the three materials is given by

with being the matrix of each thermoelastic layer given in Table 1. The elementary stiffness matrix is due to electromechanical coupling; and pyroelectric are integrated along the direction with respect to the thickness of each thermopiezoelectric layer ( or ). Thus,

Equations (41), (43), (46), and (48) are integrated numerically using the quadratic integration method [11] as follows:

where () are the coordinates of Gaussian point integration and , are the associated weight factors with and .

All elementary matrices constructed may be assembled into global matrices as follows:

where is the number of finite elements and is a distribution matrix defined by

for and with being the number of degrees of freedom of the entire structure. is the indicator vector containing the degrees of freedom 3(dof) of th node (1, 2, 3, or 4; Figure 1) in th finite element given by

Considering that there are actuators and sensors distributed on the surface of the plate, (40) including damping matrix can be written in the following global shape:

where is the distribution matrix indicating the position of the th element on the plate and is composed of numbers 1 or 0. The figure 0 means that no thermopiezoelectric element (actuator or sensor) is present while the figure indicates the presence of the element. is the number of finite elements of th thermopiezoelectric actuator/sensor and is the vector of nodal displacement of the entire structure.

We have no voltage applied to the sensor (), so the electrical potential generated is calculated through (55) for as follows:

The total voltage is composed of tension sensed by the sensor, tension sensed by the actuator, and finally the applied voltage . So, the total voltage is given by

The global dynamic equation can be formed by substitution of (56) into (57) and then into (54). The introduction of electrical forces due to the actuator forces, pyroelectric thermoelastic forces, and external mechanical forces leads to ()

where and are the thermoelastic and pyroelectric forces, respectively, given by

and and (electrical force) are given, respectively, by

for . The electric stiffness matrix can be written as

The static equation can be written using (33) as follows:

For the transformation into a state space model, equations of motion can be written as

where is the matrix structure, is the control matrix, is the matrix of disturbances, and is the thermal disturbance, given, respectively, by

3. Results and Discussions

Thermopiezoelastic behavior and the active control of cantilever aluminium plate are presented. The model considered is a flexible piezocomposite rectangular plate size (0.6 m × 0.4 m × 0.001 m) whose geometric and physical properties are shown in Table 1 and in which two piezoelectric layers (size 0.25 m × 0.25 m × 0.0002 m) are perfectly bonded (Figure 2). In the finite element modeling, the plate is divided into 60 × 40 finite elements of equal size and each element has three mechanical, one electrical, and one temperature dofs. The study consists of three phases; the first is to verify the actuation in the static case by exciting the actuator by constant and harmonic inputs; the second phase corresponds to analyse the responses of the structure dynamics under a center plate pulse and study the effect of the location of piezoelectric elements on the effectiveness of control; and the third phase was dedicated to analyzing the effect of thermoelastic and pyroelectric coupling, by applying thermal gradient, on the dynamics of the structure and the control procedure.

Figure 2: Actuation by application of constants or harmonic tensions.

3.1. Case of Absence of Thermal Effect

To validate the mathematical formulation and the finite element model in the static case, the configuration shown in Figure 2 gives an illustration of the mode of actuation of the actuator(s) bonded on the surface of the plate. This allows us to check the actuation states of the actuators and their acting on the plate. To this end, we apply a voltage (constant or harmonic) at the surfaces of the actuators causing a vibration plate.

We have given in Figure 3 the spectra of the three first modes 11, 12, and 21 and in Figure 4, for illustration, some examples of the first modes shape of vibration of the plate, that is, modes 11, 22, 33, 44, 12, 21, 13, and 14. Figure 5 shows the response to the actuator excitations by tensions with constant or different amplitudes. Figure 6 illustrates the state of the plate after a harmonic input voltage whose amplitudes and frequencies are chosen arbitrarily.

Figure 3: Spectra of first modes 11, 12, and 21 of plate dynamics in open and closed loop. Mode 11 at 3.31 Hz, mode 12 at 34.7 Hz, and mode 21 at 56.4 Hz.

Figure 4: First modes shape of a plate drawn in 3D.

Figure 5: Static test of the used finite element and control methods by giving to the actuator in different sizes pulses with different amplitudes.

Figure 6: Plate dynamics under actuator harmonic input with change in position and size of actuator patches.

In order to verify the characteristics of the actuators, and by varying their size and position, and applying the two types of voltage excitations, it was shown that the plate can generate dynamics which can be controllable. In Figure 7, the symmetrical distributed pairs of patches (actuator/sensor) in several positions and with different sizes on both surfaces of the plate are placed. The center plate displacement in mode 11 with the corresponding actuator voltage and sensor output is shown in Figures 8 and 9. The LQR-Kalman method parameters (weight matrices and with the gain control ) are taken equal to 10⁶, 1, and 200, respectively. Figure 10 gives the first modes of plate vibration with and without control. Moreover, we conclude that the control can be applied during the vibration of the plate without diminishing its quality or effectiveness as shown in Figure 11. The control is applied 0.5 seconds after the start of the oscillation. We note that the plate is quickly damped and consequently the actuator input tension begins to decrease after the application of control. In another side, the position and size of piezoelectric patches may increase or decrease the quality control. Indeed, we deduce from Figure 12 that the wider and nearer to the center of plate the piezoelectric patches are (point undergoing large deformations), the better the structure is controlled. This is justified by the fact that the sensor (mounted in front with the actuator) captures the highest deformations in the plate vibration so that the patches near to the sides of the plate capture the low vibration amplitudes, which affects the output voltage sensor, and therefore the responses of the actuator would be less effective. This is explained by the fact that the patches must be placed as possible near to the the pulse application point. However, if external perturbations are unpredictable, trials must be established to obtain the best positions to the control.

Figure 7: Closed loop of the control (sensor-controller-actuator) used in simulation.

Figure 8: Response of a pulse (3 N) of displacement of center of plate and corresponding actuator input voltage (, , and ).

Figure 9: Response of sensor output voltage in open and closed loop (, , and ).

Figure 10: Shape of first modes 11 and 33 with and without control in 3D (, , and ).

Figure 11: Application of control after 0.5 seconds from the vibration beginning (, , and ).

Figure 12: Influence of the position distribution of actuators and sensors on the quality of control (, , and ).

3.2. Case of Presence of the Thermal Effect

We assume at first that the plate and the layers thermopiezoelectric actuator/sensor are placed in a thermal field and the temperature of the plate reached quickly the studied case. From (23), the change in temperature can generate a voltage in the thermopiezoelectric sensor and the mechanical deformation created under the effect of temperature also produces an electrical signal in the sensor. In order to verify these changes in the deformation, we apply successively, in static case, a thermal gradient in the surface of the plate so that the top surface is hotter than the bottom surface. Because of this gradient, the centerline of the plate is deformed upward. Figure 13 gives the deformation of the plate for different values of the gradient of temperature and illustrates the thermoelastic effect. We recall that the sensor must be more sensitive than the actuator (the coefficient of rigidity of the actuator must be higher compared to the sensor; see Table 1).

Figure 13: Plate neutral plan deformation for different values of temperature.

We apply, now, at the center of the flexible plate a pulse with amplitude 2 N. After a period of 0.5 seconds of applying the impulse, a control based on the method LQG-Kalman is applied. The thermal effect due to a temperature gradient of 4∘C is considered in this case (Figure 14). The weight matrices and and gain control are taken equal to 10¹⁰, 1, and 200, respectively. We assume that the material with which the piezoelectric patches are bonded has a good resistance to temperature change (the patches remain perfectly bonded to the surfaces of the plate). Figure 15 illustrates the compensation due to the application of the gradient without control. The amplitudes of vibration of the neutral plane centre stabilize at the amplitude −2.5 mm. This is logical because the temperature above the plate is high compared to that below. The plate continues to oscillate in the same amplitude. Consequently, the tension in the actuator increases at the application of the gradient due to the increased mechanical deformations by thermal expansion in the sensor and also the tension due to pyroelectric effect. Figure 15 gives also the corresponding actuator voltage that increases when the gradient is applied.

Figure 14: Application of a thermal gradient between surfaces of the plate.

Figure 15: Displacement and the corresponding actuator input voltage of plate center and the actuator tension after the application of a 4∘C thermal gradient.

We take the same position as before, and we apply the control (Figure 16). We note that the plate is checked just to the deviation amplitude (−2.5 mm). The corresponding tension in the actuator is shown in Figure 16. The only precaution is to not let the gradient be applied for a long time; otherwise, we lose the controllability monitoring. It should be noted here that the tension needed to be applied to piezoelectric layers is very important to offset a small gradient of 4∘C, for example. In addition, excessively high voltages applied to the actuator must be avoided in order to not depolarize the thermopiezoelectric material (PZT). This means that thermal gradients applied must be limited by the used material of the actuator. We also note that the application of a constant gradient for a long period cannot be monitored or controlled by the LQG method, but a sudden applied gradient is controllable. This is one of the advantages of this method, since all the thermal disturbances are unpredictable. Figure 17 shows this effect. After a period of 0.5 seconds of applying the pulse (2 N), the plate is subjected suddenly to a temperature gradient of 4°C. The control is applied a second later. The structure is controlled during and after the application of the gradient. This shows that the used control method is effective for sudden changes in temperature. The tension with which the actuator acts is due to the pyroelectric effect (by increasing the temperature of the actuator itself) and the amplified voltage from the controller. It receives a voltage sensor that is the result of the thermal effect (pyroelectric effect) and the effect of increasing mechanical and thermal deformations. The voltage then increases when a thermal gradient is applied and the intensity of control increases accordingly (Figure 17).

Figure 16: The noninfluence of the application of thermal gradient on control effectiveness (, , and ).

Figure 17: Application of sudden thermal gradient during the application of control (, , and ).

4. Conclusion

The modeling by the finite element method according to the theory of Kirchhoff of a rectangular cantilever piezocomposite plate with bonded thermopiezoelectric elements at its surfaces is presented and implemented. The optimal control based on the method of LQG-Kalman is applied and discussed. The influence of the location of the patches sensor/actuator and of the application of thermal gradient on the effectiveness of control is demonstrated. It is shown that the application of control during the vibration of the structure does not diminish the control quality; that is, the control can be applied at any time during the vibration of the plate. Moreover, the deformations produced by the sudden application of a thermal gradient can be controlled, to point out that the simulation programs diverge when an intense, or of long duration, thermal gradient is applied. Otherwise, this can infect the pyroelectric effect in the sensor or may deteriorate the polarization of the actuator.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

A. Görnandt and U. Gabbert, “Finite element analysis of thermopiezoelectric smart structures,” Acta Mechanica, vol. 154, no. 1–4, pp. 129–140, 2002.

S. S. Vel and R. C. Batra, “Generalized plane strain thermopiezoelectric analysis of multilayered plates,” Journal of Thermal Stresses, vol. 26, no. 4, pp. 353–377, 2003.

G. L. C. M. de Abreu, J. F. Ribeiro, and V. Steffen Jr., “Finite element modeling of a plate with localized piezoelectric sensors and actuators,” Journal of the Brazilian Society of Mechanical Sciences and Engineering, vol. 26, no. 2, pp. 117–128, 2004.

J.-F. Deü and A. Benjeddou, “Free-vibration analysis of laminated plates with embedded shear-mode piezoceramic layers,” International Journal of Solids and Structures, vol. 42, no. 7, pp. 2059–2088, 2005.

N. Tanaka and T. Sanada, “Modal control of a rectangular plate using smart sensors and smart actuators,” Smart Materials and Structures, vol. 16, no. 1, pp. 36–46, 2007.

M. A. Trindade and A. Benjeddou, “Effective electromechanical coupling coefficients of piezoelectric adaptive structures: critical evaluation and optimization,” Mechanics of Advanced Materials and Structures, vol. 16, no. 3, pp. 210–223, 2009.

M. Sanbi, R. Saadani, K. Sbai, and M. Rahmoune, “Thermoelastic and pyroelectric couplings effects on dynamics and active control of smart piezolaminated beam modeled by finite element method,” Smart Materials Research, vol. 2014, Article ID 145087, 10 pages, 2014.

H.-J. Ding, F.-L. Guo, and P.-F. Hou, “A general solution for piezothermoelasticity of transversely isotropic piezoelectric materials and its applications,” International Journal of Engineering Science, vol. 38, no. 13, pp. 1415–1440, 2000.

J. Erhart, E. Kittinger, J. Privratska, and J. Tichy, Fundamentals of Piezoelectric Sensorics: Mechanical, Dielectric, and Thermodynamical Properties of Piezoelectric Materials, Springer, 2010.

A. Erturk, Electromechanical modeling of piezoelectric energy harvesters [Ph.D. thesis], Virginia Polytechnic Institute and State University, 2009.

K. J. Bathe, Finite Element Procedures in Engineering Analysis, Prentice-Hall, 1982.

J. S. Yang and R. C. Batra, “Free vibrations of a linear thermopiezoelectric body,” Journal of Thermal Stresses, vol. 18, no. 2, pp. 247–262, 1995.

V. Lopes Jr., J. A. Pereira, and D. J. Inman, “Structural FRF acquisition via electric impedance measurement applied to damage location,” in Proceedings of the 18th International Modal Analysis Conference (IMAC '00), European University Association (EUA), San Antonio, Tex, USA, February 2000.