Thermoelastic and Pyroelectric Couplings Effects on Dynamics and Active Control of Smart Piezolaminated Beam Modeled by Finite Element Method

Smart structures with integrated sensors, actuators, and control electronics are of importance to the next generation highperformance structural systems. In this study, thermopiezoelastic characteristics of piezoelectric beam continua are studied and applications of the theory to active structures in sensing and optimal control are discussed. Using linear thermopiezoelastic theory and Timoshenko assumptions, a generic thermopiezoelastic theory for piezolaminated composite beam is derived. Finite element equations for the thermopiezoelastic media are obtained by using the linear constitutive equations in Hamilton’s principle together with the finite element approximations. The structure consists of a modeling of cantilevered piezolaminated Timoshenko beam with integrated thermopiezoelectric elements between two aluminium layers. 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 accompanied by the Kalman filter is applied. The effects of thermoelastic and pyroelectric couplings on the dynamics of the structure and on the control procedure are studied and 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 beam.


Introduction
In the development of distributed sensors, actuators, and thin-film devices, thin layer piezoelectrics (either laminated, deposited, or embedded) are of importance in many applications, for example, dynamic measurement, control, actuation, and so forth.In case studies of electromechanical coupling and implementation of the finite element method, several studies have been the subject of research on the topic covered in this paper.Aldraihem and Khdeir [1] have studied smart beams with extension and thickness-shear piezoelectric actuators.Trindade et al. [2] have investigated the piezoelectric active vibration control of damped sandwich beams.Gabbert et al. [3] have implemented the modeling, control, and simulation of piezoelectric smart structures using finite element method and optimal LQ control.Raja et al. [4] have analysed the active vibration control of composite sandwich beams with piezoelectric extension-bending and shear actuators.Moita et al. [5] have studied the active control of adaptative laminated structures with bonded piezoelectric sensors and actuators.Manjunath and Bandyopadhyay have used the technique of fast output sampling feedback in the control of vibrations in SISO based Timoshenko structures [6,7].Trindade and Benjeddou [8] have evaluated and optimized the effective electromechanical coupling coefficients of piezoelectric adaptive structures.
Besides mechanical and electric couplings and interactions, temperature can also influence the performance of piezoelectric devices and its variation can introduce voltage/charge generation in piezoelectric sensors.In addition, control voltage can cause temperature rise in piezoelectric actuators.Temperature can introduce the pyroelectric effect and the thermal strain effect to the distributed sensors and also thermal deflection in dynamic oscillations.Aouadi [9] has discussed the generalized thermopiezoelectric problems with temperature-dependent properties.Ganesan and Sethuraman [10] have studied the thermally induced vibrations of piezothermoviscoelastic composite beam with relaxation times and system response.Sadek and Abukhaled [11] have implemented the optimal control of thermoelastic beam vibrations by piezoelectric actuation.In this study, distributed sensing and control of a piezolaminated composite beam under sudden and intense thermal gradient have been studied and sensing/control demonstrated.The present work investigates the influence of the thermal and pyroelectric coupling on the dynamic behavior of the flexible composite piezolaminated beam and on the control procedure by the application of a thermal gradient on the faces of the structure.For this purpose, we consider a piezolaminated Timoshenko's beam with sandwiched thermopiezoelectric sensors and actuators placed at different positions.The structure is modelled by finite element method where the linear constitutive equations in Hamilton's principle together with the finite element approximations are used.We look first for the effects of changes in temperature and position sensor on quality control by varying the sensor position along the beam.We are secondly looking if there is a decrease in the control quality when it is applied after the start of vibration of the beam.We demonstrate that the active control of the beam is influenced by the variation of temperature.

Basic Equations and FE Method Implementation
In this study, we consider a structure consisting of a cantilever composite piezolaminated beam ( ×  ×  = 0.2 × 0.03 × 0.001) with shear thermopiezoelectric layers (actuator or sensor) and a rigid foam both placed at the core of the structure and sandwiched between two relatively thick aluminum  layers (Figure 1).The foam is introduced to fill the space between aluminium layers in order to obtain a compact beam.To obtain a better coupling between the main structure and piezoelectric layers, the center layers (thermopiezoelectric + foam) are considered perfectly bonded.Thickness, mass, and stiffness of the adhesive are considered relatively negligible.The beam is devised in 5 two-node finite beam elements (FE) where the actuator is fixed at the FE1 while the sensor occupies successively the four other finite elements (Figure 1).The physical, piezoelectric, pyroelectric, and thermal characteristics of the used beam and thermopiezoelectric elements (sensor and actuator) are given in the Table 1.
We start by giving equations of displacement, strain, and stress; then we give the constitutive equations and Hamilton's principle to the structure.In this work the following linear constitutive relations for thermopiezoelectric materials are employed [12]: where the superscript  means that the values are measured at constant strain and the superscript  means that the values are measured at constant electric field, {} is the stress tensor, {} is the electric displacement vector, {Θ} is temperature, {} is the entropy, and {} is the electric field.{} is the strain tensor, [  ] is the elastic constants at constant electric field, [] denotes the piezoelectric stress coefficients, [] is the dielectric tensor at constant mechanical strain, [] is the thermoelastic tensor, [] is the pyroelectric tensor, and α is the expansion coefficient with α =    0 /Θ 0 where  0 and Θ 0 are the specific heat and initial temperature, respectively.
Hamilton's principle is employed here to derive the finite element equations: where  1 and  2 are two arbitrary instants,  is the kinetic energy,  is the potential energy,   denotes the work done by electrical forces, and  th is the work done by thermal forces.The total kinetic energy  and the potential energy  of the composite beam are described by the following relations: The work done by electrical forces and thermal forces and the element virtual works done by appliqued surface forces {  } are given by The mass moments characteristics of the cross-section of the beam are defined as with   and ℎ 1,2,3 being the mass density and the height of the beam + piezoelectric patches (the thickness of the total structure), respectively.The dynamic equation can be found in [13], and the displacements () and () are written as where  0 ,  0 (), and () are the axial, transverse midplane displacements and -rotation, respectively [14].Assuming that there is no compressibility in the  direction, the normal and transverse components of strain are: The beam constitutive equations can be written as where where   =  11 is the extension stress,   =  55 ,   is the shear stress,  is the beam width,  is the distance measured between the plane of the structure and that of the th layer laminate, ℎ is the total thickness of the structure (beam + actuator/sensor + beam),   is the axial force,   is the bending moment,   is the shear force,  11 ,  11 ,  11 ,  55 are the extension, extension-bending, bending, and transverse shear stiffness coefficients given by [13,15] where   is the distance of the th layer relative to the -axis,  is the number of layers, K is the correction shear deformation factor generally taken to be 5/6, and  11 ,  55 are calculated based on the physical properties of piezoelectric material [13,15,16]: where is the vector of nodal displacements and [  ], [  ], [  ] are the mode shape functions due to the axial displacement, transverse displacement, and the slop, which are defined as [17] [ The elements of the shape functions are given in [17].The inertial forces vector N can be written as The mass matrix of the regular beam element is given by where is the inertia matrix.Similarly, the stiffness matrix of the regular beam element can be written as where   is the area of cross-section of the beam element and Temperature flux Application of the thermal gradient

Control box
Actuator Sensor The mass and stiffness element matrices of the piezoelectric element are obtained in the same manner with respect to the physical characteristics of the material.The mass and stiffness of the finite elements which together contain the element of the beam and the piezoelectric element are given by  =   + 2  and  =   + 2  .The sensor output voltage, due to thermal and mechanical deformations and temperature (pyroelectric effect), can be written as where the control force developed by the thermopiezoelectric actuator is written where G is the transverse module, ℎ = (ℎ  + ℎ  )/2 is the distance between neutral axes of the beam and the thermopiezoelectric layer,   is the shape function of rotations, and   () is the actuator input voltage.
Similarly, the forces due to thermoelastic and pyroelectric couplings f eth and f pth are given by where  Θ is thermal shape function, [] is the thermoelastic tensor, and [p] is pyroelectric tensor.If an external force f ext is applied, the total force acting on the beam is Therefore, the equation of motion can be written as where M * , K * , C * , and g are the generalized mass, stiffness, damping matrices, and generalized displacement.The above with where T is the modal matrix, () is the command vector, and () is the external force vector.

Results and Discussion
The graphs presented below correspond to the responses of the free end of the beam.We compare noncontrolled and controlled responses for a disturbance in impulse (2N) before, during, and after applying a thermal gradient so that the faces of the structure are under two different temperatures as illustrated in Figure 2. We assume that the material with which the material bonded the thermopiezoelectric elements is resistant to the temperature change (the patches remain perfectly bonded to the surfaces of the beam).A thermal gradient will be applied to visualize the thermal and pyroelectric effect on the vibration behavior of the structure under thermal perturbations.Figures 3 and 4 show the three first modes.The different responses, step, impulse, and sinusoidal are illustrated in the Figures 5, 6, and 7. Figure 8 shows 3D control visualisation for the three first modes and for (2N) pulse excitation at free end of the beam.
Figure 10 shows in 3D the deflection of the beam under a 5 ∘ C thermal gradient.In addition to the sensor voltage created by the expansion of the beam under the thermal field, another voltage is created by the temperature increase (pyroelectric effect).In fact, respectively, Figures 9 and 10 show the variation of amplitude and voltage at the input of the actuator for the beam before and during the control applied 0.5 sec after the application of the pulse.The control method can be applied during the vibration without diminishing its quality or effectiveness.The same above figures illustrate the thermal compensation of responses in amplitude and voltage of the actuator.The gradient is applied at the beginning of the pulse after 0.5 sec.We noticed that the application of such a gradient did not alter the quality and effectiveness of control.The effectiveness of control will depend on the intensity of the thermal gradient and its duration.Indeed, the intense gradients can infect the pyroelectric effect in the thermopiezoelectric elements and deteriorate their polarization.We finally conclude that the control method used is better adapted to sudden changes in disturbance types of vibration or for small variations in the thermal gradient.To show the influence of thermal effect on the used control method, we maintain a gradient for 3 sec.We note that the control procedure is effective in applying or removing the thermal gradient.However, the application of a constant gradient for a long duration can not be monitored or controlled by the method LQG.Instead, the sudden application of a gradient is controllable.This is one of the advantages of this method, since all the thermal disturbances are unpredictable.We conclude, therefore, that the control method used is effective for sudden changes in temperature.

Conclusion
Modeling by the finite element method according to the theory of Timoshenko, of a cantilever piezolaminated composite beam with embedded thermopiezoelectric element, is presented.The optimal control based on the method of LQG-Kalman is applied and discussed.The influence of the location of the sensor and of the application of thermal gradient on the effectiveness of control is analyzed.In fact, our analysis shows that the more the sensor is near the free end of the beam, the more the control is effective.This is due to the the increase of the deformation amplitudes in the sensor, which affects the control voltage of the actuator.We have also 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 beam.Moreover, the deformations produced by the sudden application of a thermal gradient can be controlled.We reported that the application of an intense thermal gradient, or of long duration, can infect the pyroelectric effect in the sensor or may deteriorate the polarization of the actuator.

Figure 1 :
Figure 1: Cantilever beam with embedded thermopiezoelectric elements and the different models.

Figure 2 :
Figure 2: Application of a thermal gradient.

Figure 3 :
Figure 3: The first three modes of the beam subjected to a pulse of 2N at its free end.

Figure 4 :
Figure 4: Spectra of the first three modes for the sensor positions EF5 and EF2 in SISO model, respectively.

Figure 5 :Figure 6 :
Figure 5: Influence of the sensor location on the quality of control in the case of a pulse excitation ( = 10 7 ,  = 1, and controller gain   = 200) for the two first modes.

Figure 9 :Figure 10 :
Figure 9: Application of control during the vibration of the beam.

Table 1 :
Properties of the aluminium and piezoelectric layers.
11 cos 4  +  22 sin 4  + 2 ( 12 +  66 ) sin 2 cos 2 ,  55 =  13 cos 2  +  23 sin 2 .(11) The angle  is the angle between the direction of the fibers and the longitudinal axis of the beam.The physical constants  11 ,  22 ,  12 ,  66 ,  13 ,  23 relative to the foam, aluminum, and the piezoelectric material are  11 , and  55 are the induced piezoelectric axial force, bending moment due to deformation of the actuator, and the shear strength given by [13, 15]  11 =  11 = 0 since the piezoelectric layers are polarized longitudinally.  (, ) = the voltage applied to the th actuator with thickness (  + −   − ) and piezoelectric constants   31 and   15 .  = the number of actuators.Consider where ] is the poisson coefficient and  is the rigidity transverse modulus.Respectively,  11 ,