Next Article in Journal
Simple Approach to the Fabrication of Lanthanum Orthoniobates and Nanocomposites with Ni, Cu, and Co Metal Nanoparticles Using Supercritical Isopropanol
Next Article in Special Issue
Localization of Vibration Weak Position of Composites Based on Weighted Modal Strain Energy Summation
Previous Article in Journal
Special Issues on Composite Carbon Fibers
Previous Article in Special Issue
A New Study on the Structure, and Phase Transition Temperature of Bulk Silicate Materials by Simulation Method of Molecular Dynamics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Robust Control and Thermal Analysis of a Reduced Model of Kirchhoff Composite Plate with Random Distribution of Thermopiezoelectric Sensors and Actuators

1
Sciences and Advanced Technologies Team, National School of Applied Sciences ENSATE, Abdelmalek Essaadi University UAE, Tetouan 93000, Morocco
2
Laboratory for the Study of Advanced Materials and Applications, FSM-ESTM, Moulay Ismail University, Meknes 50000, Morocco
*
Author to whom correspondence should be addressed.
J. Compos. Sci. 2022, 6(8), 242; https://doi.org/10.3390/jcs6080242
Submission received: 4 July 2022 / Revised: 7 August 2022 / Accepted: 12 August 2022 / Published: 18 August 2022
(This article belongs to the Special Issue Dynamic Performance and Response of Composite Based Structures)

Abstract

:
This paper presents an implementation of a robust control LQG-Kalman model applied to composite Kirchhoff plate dynamics. A reduced model of a finite element method and control procedure is considered in the modeling of a structure because of the important number of piezoelectric patches used in control. Replacing the full model with a short model reduces the computational and time costs, especially when the number of degrees of freedom is significant. In robust control, the measurement of all states is not necessary and the observability and estimability criteria can be exploited, while conventional LQR control assumes that the data accessibility of all states is available. For this reason, robust control is proposed to control the random external disturbances and is compared to LQR control to illustrate its practicability and efficiency. The sensors and actuators in the thermo-piezoelectric material are randomly distributed on both sides of the plate to establish the control procedure. A Monte Carlo simulation is used in the selection of the degrees of freedom of sensors presenting high electrical outputs. Numerical simulations are performed to demonstrate the effectiveness of the proposed control procedure in a reduced model and under mechanical and thermal disturbances in comparison with the LQR control.

1. Introduction

The progress of modeling of smart composite structures offers great improvement in advanced engineering systems. The employment of piezoelectric materials as actuators and/or sensors in the modeling of these structures is generally very useful. In analysis of smart structures, the coupling issues associated with thermal, piezoelectric and mechanical fields play an important role. In the literature, the piezo-thermo-elastic behavior of a laminated plate using classical lamination theory (CLT) was evaluated by [1], and, using layer wise theory, by [2]. First-order shear deformation theory was applied to the active control of plates by [3,4]. The control of plates using 3D solid elements was carried out by [5,6], considering the computation cost due to an increased number of degrees of freedom. A linear quadratic regulator (LQR) controller which does not require collocated actuator/sensor pairs was employed for vibration control design for beams and plates in [7,8]. Tzou et al. [9] used first-order shear deformation theory (FSDT) to develop a dynamic model of a thermo-piezo-elastic shell laminate. Sanbi et al. [10,11] studied the pyroelectric effect on the quality of control for beams and plates modeled using a finite element method. A hybrid laminate theory was applied to model smart composite laminates by Mitchell et al. [12]. For thin and relatively thick laminates, Gu et al. [13] developed a higher-order theory which included transverse shear effects. Many of the studies conducted are concerned with model reduction. In this respect, nonlinear systems with time-varying and uncertain variables with modal/balanced truncation are described in [14]. For nonlinear models, an analysis and comparison of linear principal component analysis (PCA) and nonlinear principal component analysis (NLPCA) are provided in [15,16]. For linear structural models, the system equivalent reduction expansion process (SEREP) can be used, some descriptions of which are provided in [17].
With respect to control implementation in recent decades, much research has been undertaken in the field of active vibration control. Many studies have investigated the procedures to make control algorithms better and faster [7,18,19]. Various procedures have been proposed in the literature to analyze the optimal placement of sensors and actuators [20,21]. Some papers have considered the effect of environmental variables on smart structures [10,11,22,23].
The application of reduced-order models is more advantageous in the design of controllers. Generally, increasing the order of the system increases the order of the controller; a relatively high order of controller is complicated to design because of the large number of states and can cause system instability. Reducing the order of a system is very useful to ensure feasible and robust control and to reduce the required computational effort. It is also used to reduce the computational cost and the complexity of computations in modeling and the control of structures without change to their accuracy. The disadvantage of the LQR scheme is that it requires the measurement of all the state variables while in the linear quadratic Gaussian (LQG) controller—this measurement of all states is not necessary. The LQG-Kalman design is implemented in three steps: first, construction of a linear quadratic regulator (LQR) to obtain the control input relative to a measured state; second, construction of a Kalman estimator to provide a best estimation of the states; and third, combination of the optimal regulator and the Kalman filter in one optimal controller/estimator to obtain the input estimated state vector and the measured output vector.
This paper presents an analysis of a reduced model process applied to an FE model and to the active vibration control of a thermo-piezo-elastic composite plate. In the process of model reduction, effort is made to conserve, in reduced space, the modal and dynamical representations of the full space of the studied structure. The classical theory used in the modeling of the plate is based on Kirchhoff’s assumptions, neglecting the effects of transverse shear deformation. Hamilton’s principle is then used to derive the governing equations of motion. A finite element method is established using a four-node rectangular element which is based on first-order theory. The resulting FE equations in nodal variables are transformed into modal form and then introduced into the state space to design the controller. Piezoelectric patches of actuators and sensors, covering one or more finite elements, are distributed randomly along the upper and bottom surfaces of the plate. Sensors are placed at the bottom surface of the plate to recover the maximum number of deformations and deflections. The location of the sensors along the structure is performed using a Monte Carlo simulation (MCS). The optimal LQR and LQG-Kalman controllers are designed and their gain is calculated using output state law. Numerical results are presented to illustrate the performances of a robust control approach represented using a reduced model in a thermal environment. Some results are compared with the full model and with LQR control in the same mechanical and thermal circumstances.

2. Mathematical Modeling

2.1. Kirchhoff Plate Model

In this formulation, the piezoelectric layers are assumed to be perfectly bonded together and the formulation is restricted to linear elastic material behavior. A quadrilateral Kirchhoff plate element with four nodes is considered. Each node has three mechanical degrees of freedom (dofs), the displacement in the z-direction and two rotations in the x- and y-directions, one electric dof ϕ , and one temperature dof θ . Using the Kirchhoff assumptions for thin plates, the displacement fields u, v, and w can be expressed as
u = z w x ; v = z w y ; w = w ( x , y ) .
where x and y are the axes at the mid-planes of the plate, and z is along the plate thickness direction, as seen Figure 1. u and v are the displacements along the x- and y-axes, respectively, and w is the transverse displacement along the z-axis. Neglecting the shear deformations, the strains can be written as a function of the displacements as
ε = ε x x ε y y ε x y T = z 2 w x 2 2 w y 2 2 w x y T .

2.2. Constitutive Equations of Thermopiezoelectric Material

In this paper, the following linear constitutive relations for thermo-piezoelectric materials are employed [24]:
σ = C E ε e E λ Θ ,
D = e T ε ϵ E + p Θ ,
s = λ σ ε + p T E + α ˜ Θ .
where σ is the stress tensor, D is the electric displacement vector, Θ is the temperature, s is the entropy, ε is the strain tensor, E is the electric field, C E are the elastic constants, e denotes the piezoelectric stress coefficients, ϵ is the dielectric tensor, λ is the thermo-elastic tensor, p is the pyroelectric tensor, and α ˜ is the expansion coefficient with α ˜ = r h o p c 0 / Θ 0 where c 0 and Θ 0 are the specific heat and initial temperature, respectively.

2.3. Dynamic Equation

To derive the FEM equations, Hamilton’s principle is employed here and can be written as
t 1 t 2 δ ( T U + W e W t h ) + δ W d t = 0 .
where T is the kinetic energy, U is the potential energy, W e represents the work of electric forces and W t h is the work induced by thermal forces. The total kinetic T and the potential U energies of the structure are described as follows
T = 1 2 V ρ { q ˙ } T { q ˙ } d V ,
U = 1 2 V { ε } T { σ } d V V { ε } T { λ } { Θ } d V .
The element works produced by electric forces and by external volume/surface forces { f v } and { f A } are given by
W e = 1 2 V ( { E } T { D } + { p } { Θ } ) d V ,
δ W = V { δ q } T { f v } d V + A { δ q } T { f A } δ Φ σ q d A .
The generalized dynamic equation of the plate and the distributed sensors equation in the global form are
M g ¨ + C q q g ˙ + K q q g + K q Φ Φ K q Θ Θ F = 0 ,
K Φ q g + K Φ Θ Θ + K Φ Φ Φ + Q a = 0 .
where g is the generalized displacement, M is the mass matrix, C q q is the damping matrix, K q q is the stiffness matrix, K q Φ is the electrical-mechanical coupling stiffness matrix, K q Θ is the thermo-elastic matrix, K Φ Θ is the pyroelectric matrix, K Φ Φ the dielectric stiffness matrix, F is the external force and Q a are the electrical charges.

3. Finite Element Formulation

3.1. Full Order Process

In this section, the host plate and the piezoelectric elements are modeled by proceeding to a full-order FE method. Using a finite element with four nodes and Kirchhoff’s classical plate theory, the strains developing in the plate can be written as
{ ε } 3 × 1 = z [ β u ] 3 × 12 { u e } 12 × 1 .
where
[ β u ] 3 × 12 = z 2 x 2 z 2 y 2 2 z 2 x y 3 × 1 T × [ N ] 1 × 12 .
in which [ N ] 1 × 12 is the interpolation function and u e is the element displacement vector. Since the control patches are made using the piezoelectric material PZT-5H for the actuators and PVDF for the sensors, and assuming that the piezoelectric patches are poled in the z-direction, the electric field can be written as
{ E } 3 × 1 = { β ϕ } V = [ 0 0 1 / h p ] T V .
Applying Hamilton’s variational principle, the dynamic equation of one selected finite element of a structure (plate-piezoelectric patches) is given by
[ M e ] { u ¨ e } + [ K e ] { u e } = { F e } .
where [ M e ] and [ K e ] are the elemental mass and stiffness matrices. The total force vector { F e } acting on the finite element is
[ F e ] = F S e + F T , S e + F T , P e + F E , P e .
where the mechanical, the plate material thermal expansion, the piezoelectric material thermal expansion and the electric forces, respectively, are expressed by
{ F S e } = A S [ N ] T { f S e } d A S ,
{ F T , S e } = 1 / 2 S z [ β u ] T [ C S ] { α S } Δ T d τ ,
{ F T , P e } = 1 / 2 P z [ β u ] T [ C P ] { α P } Δ T d τ ,
{ F E , P e } = ( [ K u ϕ e ] + [ K ˜ u ϕ e ] ) Z s 1 ( Q e x t e + Q p y r e ) .
where
[ K u ϕ e ] = [ K ϕ u e ] = P z { β u } [ e T ] [ β ϕ ] d τ ,
[ K ˜ u ϕ e ] = [ K ˜ ϕ u e ] = P z { β u } [ e ˜ ] [ β ϕ ] Δ T d τ .
here, [ K u ϕ e ] is the electromechanical interaction matrix and [ K ˜ u ϕ e ] is the change in the electromechanical interaction matrix at the applied temperature. The external charge applied on the piezoelectric surface Q e x t e , the produced charge on the piezoelectric patch due to the pyroelectric effect Q p y r e and the sensor capacitance at a reference temperature Z s , are given by
Q e x t e = A P q d A P ,
Q p y r e = 1 2 P { P } T { β ϕ } Δ T d τ ,
Z s = ( [ K ϕ ϕ e ] + [ K ˜ ϕ ϕ e ] ) .
where [ K ϕ ϕ e ] is the capacitance of the piezoelectric patches at the reference temperature and [ K ˜ ϕ ϕ e ] is the temporal relative change in this capacitance at a given temperature:
[ K ϕ ϕ e ] = P { β ϕ } T [ ε T ] { β ϕ } d τ ,
[ K ˜ ϕ ϕ e ] = P { β ϕ } T [ ε ˜ ] { B ϕ } Δ T d τ .
The electrical forces on a piezoelectric patch are due to the thermal strain effect, the pyroelectric effect, and external charges. The voltage produced across these patches can be written as
V = [ K ϕ ϕ e ] + [ K ˜ ϕ ϕ e ] 1 Q e x t e + Q p y r e + [ K ϕ u e ] + [ K ˜ ϕ u e ] { u e } .

3.2. Modal Analysis

Using uncoupled equations, the analysis of the coupled equations System (16) becomes easier and the multi-degree freedom system can be replaced by a single-degree freedom system using modal analysis. The process consists in adopting the following transformation
{ X } n × 1 = [ H ] n × r { η } r × 1 .
where [ H ] is the eigen-vector matrix, { η } is the modal displacement vector, n is the total number of dofs of the structure and r is the considered number of modes. Assuming a linear damping as a combination of the mass and stiffness matrices, the uncoupled dynamic equations can be rewritten as
η ¨ m + c m η ˙ + ω m 2 η m = f m .
where c m = 2 ζ m ω m = μ + λ ω m 2 with the given modal damping ratio ζ m = μ 2 ω m + λ ω m 2 . Here, “m” denote a mode number and the ( μ , λ ) correspond to the mass and stiffness proportional Rayleigh damping coefficients.
The uncoupled Equation (31) represents the equation of motion of the m t h mode of a smart plate. The three first five modal frequencies are calculated at a reference temperature using commercial software for the derived FEM model. Applying the boundary conditions of the plate with one cantilevered edge, the Equation (17) can be transformed as
{ F } 270 × 1 = { F S } 270 × 1 + { F T , S } 270 × 1 + { F T , P } 270 × 1 + { F E , P } 270 × 1 .
Assuming that the external surface force F S = 0 and the thermal expansion force F T S = 0 , the modal force vector { f } = [ H ] T { F } , depending on the pyroelectric effect, can be evaluated. Substituting { F T , P } = F p t h and the value of { F E , P } from Equation (17) into Equation (32), we can obtain the following modal force
{ f } = [ H ] T F p t h [ K u ϕ e ] + [ K ˜ u ϕ e ] Z s 1 ( Q e x t + Q p y r ) .
where Q e x t = Z a c t V e x t is the total external charge applied on the actuator surface. The component of the modal force { f } which depends upon the applied external voltage ( V e x t ) on the actuator is known as the ‘modal control force’ and is given by
{ f c } = H T [ K u ϕ e ] + [ K ˜ u ϕ e ] Z s 1 Z a c t V e x t .
when sensor and actuator patches are subjected to the same temperature, we have Z s = Z a c t , and the modal control force becomes
{ f c } = H T [ K u ϕ e ] + [ K ˜ u ϕ e ] Z s 1 Z a c t V e x t = H T [ K u ϕ e ] + [ K ˜ u ϕ e ] V e x t .

3.3. General Reduction Model

The reduced-order modeling is a global reduction process based on the eigenvalue analysis of the FEM model. There are two reduction process levels; the modal reduction and the dofs reduction. The transformation matrix obtained using this process transforms the full model to a short model using selected eigen-modes with some arbitrary dofs.
Consider an equivalent Kirchhoff plate model element of the above section with four nodes, two mechanical dofs ( ω , θ ), one electric dof ( ϕ ) and one temperature dof ( θ ) at each node (Figure 1). The FEM model of the plate is derived using Hamilton’s principle; the compact equation of motion for an undamped system with electromechanical coupling can be rewritten as
M q q q ¨ + K q q q + K q ϕ ϕ F = 0
K ϕ q q + K ϕ ϕ ϕ Q a = 0
where K q ϕ , K ϕ ϕ , ϕ and Q a are the electromechanical stiffness matrix, the stiffness matrix due to electrical dofs, the piezoelectric actuator’s electrical vector, and the vector of the electrical charge, respectively.
Substituting ϕ from Equation (37) into Equation (36), the dynamic equation in damped and compact form is given by
M X ¨ ( t ) + C X ˙ ( t ) + K X ( t ) = F ( t ) .
where M, C, and K are, respectively, the mass, damping, and stiffness matrices of the full-order model. X ( t ) , X ˙ ( t ) and X ¨ ( t ) represent the displacement, velocity and acceleration vectors. To apply the reduced-order model, Equation (38) will be used in the next formulation. Using a modal matrix, the system can be written in terms of modal coordinates and for “j” eigen-modes as
X ( t ) = ψ ( n × j ) q ^ ( j × 1 ) ( t ) .
If only “a” eigen-modes are considered, the solution can be written as
X ( t ) = ψ ( n × a ) q ( a × 1 ) ( t ) .
The full-order model can be decomposed into master and slave degrees of freedom. For only “j” selected eigen-modes in Equation (40), we can write in terms of active dofs
X ( t ) = X m ( t ) X s ( t ) ψ m ( t ) ψ s ( t ) . q ^ j ( t ) .
where q ^ j is the modal matrix having “j” modes, and “m” and “s” denote the superscripts corresponding to the master and slave dofs. Moreover, here, ψ m ϵ m × j and ψ s ϵ ( n m ) × j . Considering only the master dofs, the responses in modal form can be given as
q ^ j ( t ) = ψ T m X m ( t ) .
Since the matrix ψ T m may not be a square, the number of modes and number of dofs in the analysis may be the same—this depends on the problem and hence the matrix. The different selection possibilities of modes and dofs are cited in [25]. Substituting Equation (42) in Equation (41), the modal vector can be written as
X ( t ) = ψ ψ T m X m ( t ) = T r X m .
where ( T r ) is the transformation matrix that transforms the full-order model into a reduced-order model. Using this transformation matrix, the full-order model matrices in reduced dimensions can be written as
M r = T r T M T r ; K r = T r T K T r ; C r = T r T C T r .
here, M r , K r , C r and F r ( t ) = T r T F ( t ) represent the mass, stiffness, damping matrices, and force vector expressed in the reduced-order model, respectively. The selection of modes and dofs are governed by the type of system and the conditions of the external charges.

3.4. Monte Carlo Method

A Monte Carlo method is used here to randomly select the number of sensors and their discrete locations along the plate. This selection is established by minimizing the variance and depends on a random points choice of the degrees of freedom of the sensor locations. Applying this distribution, the measurements are randomly allocated by the selected points. MCS requires the specification of dimensions, the FEM grid map, the density of nodes for each region in the plate and, finally, the computational number of samples of repetition n m c . For each point in the grid, the generation of the grid map, depending on prefixed dimensions and data, is randomly established from the defined statistical properties of the studied model. The selected random points are chosen from the grid map and the spatial analysis procedure is performed as a function of the different values of location. The optimized variables are determined from the mean-variance values and stored. The above steps are repeated until the minimal mean-variance is found. Therefore, the sensor locations that generated the minimum mean-variance are the most optimal sensor placements based on these repetitions.

4. Control Design

4.1. State Space Model

To control the response of a composite plate using an active vibration control procedure, the uncoupled equations of motion have to be converted into a state-space (SS) model. Therefore, in this section, the uncoupled equations are rewritten in an SS form to represent the structure in a first differential order. The state variables of the structure, hence, are given by
{ s ˙ } = [ A ] { s } + [ B ] u .
where ( s ˙ ) is the first time derivative of the state variables and u is the input of the control box. The sensor outputs are a function of { s } and are given by
{ y } = C { s } .
Converting the first three modal equations of motion into the state space form by η ˙ ( j = 1 , 2 , 3 ) = x ( j = 1 , 2 , 3 ) and using Equation (31), the first three modes can be given by
x ˙ j + c j x j + ω j 2 η j = f j , ( j = 1 , 2 , 3 ) .
or in matrix form
0 0 0 1 0 0 0 0 0 0 1 0 1 0 0 0 0 1 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 η ˙ 1 η ˙ 2 η ˙ 3 x ˙ 1 x ˙ 2 x ˙ 3 + ω 1 2 0 0 c 1 0 0 0 ω 2 2 0 0 c 2 0 0 0 ω 3 2 0 0 c 3 1 0 0 1 0 0 0 1 0 0 1 0 0 0 1 0 0 1 η 1 η 2 η 3 x 1 x 2 x 3 = f 1 f 2 f 3 0 0 0 .
where
1 0 0 0 1 0 0 0 1 = [ H ] T [ M ] [ H ] , c 1 0 0 0 c 2 0 0 0 c 3 = [ H ] T [ C ] [ H ] , ω 1 2 0 0 0 ω 2 2 0 0 0 ω 3 2 = [ H ] T [ K ] [ H ] .
are the modal mass, modal damping, and modal stiffness sub-matrices, respectively. [ H ] is the normalized eigenvector matrix. (45) can be rewritten as
{ s ˙ } = [ A ] { s } + [ B ] V e x t + [ G ] .
with { s } = { η 1 η 2 η 3 x 1 x 2 x 3 } T and [ G ] is the matrix corresponding to the thermal strain and to the pyroelectric effects. Equation (50) is the state space model equation. From Equation (29), the voltage produced by the sensors due to deformations of the plate is given by
V s e n s = Z s e n s 1 ( Q p y r + ( [ K u ϕ e ] + [ K ˜ u ϕ e ] ) u e ) .

4.2. LQR Control

A linear quadratic regulator (LQR) controller was used and implemented here for comparison with the LQG-Kalman scheme. Therefore, the SS is converted to the reduced-order state model, as shown in Figure 2. The actuator inputs are obtained through the LQR algorithm by minimizing the cost function ( J ) defined by
J = lim τ 1 τ 0 ( X m ) T Q ( X m ) + ( u m ) T R ( u m ) d t .
where Q and R are the weight matrices of the state space variable X m and the input u m . By minimization of the quadratic cost function defined in Equation (52) and using the dynamic equation, we obtain the state feedback in the form of the required control input u m as follows
u m = K s X m .
where K s is the ( m × r ) gain matrix resulting from the reduced model. The gain of a closed loop is determined by solving the Riccati reduced-matrix. With this gain matrix, the control input needed to control the full-order model dynamics of the plate can, hence, be estimated. The capacitance of sensor patch Z s can be evaluated at an initial temperature by
Z s = K ˘ ϕ ϕ e = ( [ K ϕ ϕ e ] + [ K ˜ ϕ ϕ e ] ) .
where [ K ϕ ϕ e ] is the capacitance of a piezoelectric patch at reference temperature and [ K ˜ ϕ ϕ e ] are the changes in the capacitance at a given temperature. For an initial temperature, the voltage generated across a piezoelectric patch and a modal force at reference temperature can be written as
V = [ K ϕ ϕ e ] 1 Q e x t e + [ K ϕ ϕ e ] { u e } ,
{ F c } = [ H ] T [ K u ϕ e ] [ K ϕ ϕ e ] 1 Q e x t .
Using identical sensor and actuator patches at the same temperature, the modal control force becomes
{ F c } = H T [ K u ϕ e ] V e x .

4.3. LQG-Kalman Estimator

The filter of Kalman is a computational estimator based on two cycles: time and measurement updates. It is used to identify the mechanical displacements and the parameters of a system. Its advantage is that the dynamic noise can be included in the dynamical state model. It is a real time estimator supplying the estimates for an available measurement at a given instant. The filter algorithm is generally a set of analytical equations that implement a prediction (correction-estimation) that is optimal in the sense of minimizing the estimated error covariance in some control conditions. For a studied plate, the feedback control system (LQG-Kalman filter) requires the measurement of tip deflection and its rate. By adding the filter of Kalman, the control design gives the required potential noise control. For this purpose, in the proposed scheme in Figure 3, the noise has been introduced as an auxiliary gain to be highly sensible to noise. In addition, sometimes, due to the nature of the system input to the non-linearity or modal frequency, the system response becomes unstable. This can be solved by filtering the output response in the controller feedback (using measurement data). According to [26,27], the system can be represented as a linear stochastic difference relation as
{ x ˙ } = [ A ] { x } + [ B ] u + w .
Using the state equation of the structure (50), a first modal control force, a control voltage and a voltage in the sensor are given by
{ F c } = [ H ] T [ K u ϕ e ] { V e x } ,
V e x t = [ H ] T [ K u ϕ e ] 1 k η ˙ e ,
V s e n s = [ K ϕ ϕ e ] 1 [ K u ϕ e ] { u e } .
where η ˙ e is estimated using the Kalman observer. The modal displacements and velocities are estimated using the Kalman observer/estimator and can be constructed as
{ η ˙ e } = [ A d ] { η e } + [ B d ] V e x t + [ L ] V s e n s [ X ] η e , { η ˙ e } = { η e } + [ M ] ( V s e n s [ X ] { η e } ) . .
where η e is the estimated vector in full state and [ A d ] , [ B d ] are discretized forms of [ A ] and [ B ] . [ L ] and [ M ] are the Kalman filter gain matrices. V s e n is the sensor voltage and X is the sensor location vector. The matrices [ B ] , [ X ] , [ L ] and [ M ] and the state vector η ˙ e all depend on the electromechanical and permittivity coefficients of the used piezoelectric material. The external control voltage to be applied to the actuator patch and corresponding to the modal control force is given by [28]
V c t r l = G η ˙ e [ H ] T Z a c t .
where g is the velocity gain of the controller, H is the modal matrix of the first three modes, Z a c t is the capacitance of the actuator and f c is the corresponding control force.

5. Discussion

This paper presents an analysis of the behavior of a smart piezocomposite plate, actively controlled in a reduced model and under different values of temperature. Two types of controllers are used: a robust LQG-Kalman controller and an optimal LQR controller. The Kalman filter is based on an analytical model deduced from FEM formulation and the LQR is a regulator based on velocity feedback of the first modal state.
To perform this comparison, to verify the robustness of the LQG-Kalman control, and to analyze the effects of temperature changes and the distribution of sensors, a thin one-side cantilevered square plate of 18 cm × 18 cm × 0.08 cm size is considered, as shown in Figure 1. A range of piezoelectric actuators (PZT) and sensors (PVDF) are randomly distributed along the plate. They are assumed to be perfectly symmetrically bonded upon the top and the bottom sides of the plate. The material properties of both the piezoelectric elements and the host plate are presented in Table 1. The plate is divided into 81 quadrilateral finite elements of the same size. The plate element is modeled using Kirchhoff’s theory and has four nodes (Figure 1). Each node has three mechanical dofs (one translation and two rotations), one electrical dof, and one thermal dof.
To establish the control procedure, the sensors are distributed over the entire surface of the plate and in three areas. The first area is near the cantilever edge (element numbers 1 to 27); the second area corresponds to the middle of the plate (elements 28 to 54); and the third concerns the end-edge of the plate (elements 55 to 81). For each distribution, the MCS method is used to select the dofs of elements having the maximum electrical responses, i.e., those which correspond to the significant mechanical deformations (including deformation due to temperature increase, i.e., the pyroelectric effect), which allows the actuator to act locally and not over the entire surface of the plate. The comparison of the control response is, then, accomplished between all proposed distributions. The numeric formulation and control implementation of the structure is achieved using commercial software.
The first vibration mode shapes of the plate are visualized. Impulse responses at the end of the plate are traced. The sensor voltage output, the control voltage applied to the actuator input, the total control deformations, and the state vector estimated by the Kalman observer at different temperatures are evaluated and simulated. The control is next compared to the LQR algorithm to evaluate of his efficiency. Figure 4 and Figure 5 present the first 3D mode shapes of vibration with relative frequencies of the one-edge-cantilevered plate. Figure 6 shows deformation responses using LQG-Kalman control in full and reduced order models for an ambient temperature 25 ° C. The MCS is applied here to select 20 active dofs from the distributed sensors and with 100 computational number of repetition samples. The Figure 7 compare, in the reduced order model, the impulse response of the strategies of control (LQR and LQG-Kalman) at a temperature 80 ° C with the same situation of MCS distribution and sensor positions. The comparison between full and reduced order model using LQG-Kalman procedure has been established in Figure 8. the actuator input voltages are given in Figure 9 in both full and reduced models.
Figure 10 and Figure 11 illustrate, in the reduced and full models, the sensor outputs for multiple temperatures: (25 ° C, 45 ° C, 65 ° C, and 85 ° C). Figure 12 shows a comparison of the responses of the control according to the selection of the nodes and which belong to the sensors distributed on the three parts of the plate (the embedded, middle, and free parts). It is noted that it is preferable to multiply the number of sensors to ensure a very good selection through the MCS and to obtain the maximum amount of information on all parts of the plate undergoing the deformations. These are not filtered and contribute to forming the output of all the sensors, which are then introduced into the control box, processed and transmitted to the actuators, which act on their part on the plate.

6. Conclusions

In this study, the active vibration control and thermal effects of a piezo-composite plate are analyzed. Representations of the two procedures of control, i.e., LQR and LQG-Kalman, are established in a reduced model and performed. Based on the problem formulation, the control performance depends on the sensor voltage output, the actuator control voltage (proportional to the control effort), the Kalman filter algorithm and changes in temperature. Further, to maintain the control performance for all temperatures and noise perturbations, the Kalman observer is generally added to the LQG regulator. Therefore, it is concluded that the control is more effective using the reduced model formulation compared to the full model. The simulations demonstrated the robustness of the control law with respect to temperature variations and unpredictable perturbations. The quality of control can be improved by choosing adequate controls to be used as a function of the nature of the perturbations of the structures. The control voltage applied to the inputs of the actuators depends on the number of modes to be controlled, the type of control law, and the piezoelectric material of the actuators. The selection of active dofs produced by the MCS distribution helps in the selection of the nodes with high amplitude of mechanical displacement and neglects the dofs with low amplitudes. All displacements having a low amplitude are ignored under some criteria. The simulations, relative to the reduced model, also show the feasibility, and the effectiveness, of the LQG-Kalman control under thermal conditions.

Author Contributions

Conceptualization, M.S., R.S. and M.R.; methodology, M.S.; software, L.E.K. and M.S.; validation, M.S., M.R. and R.S.; formal analysis, R.S.; investigation, L.E.K.; resources, L.E.K.; data curation, L.E.K. and M.S.; writing—original draft preparation, L.E.K. and M.S.; writing—review and editing, M.S.; visualization, L.E.K.; supervision, M.R.; project administration, M.S.; funding acquisition, L.E.K., M.S., R.S. and M.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tauchert, T. Piezotheelmoelastic behavior of a laminated plate. J. Therm. Stress. 1992, 15, 25–37. [Google Scholar] [CrossRef]
  2. Lee, H.; Saravanos, D. Generalized finite element formulation for smart multilayered thermal piezoelectric composite plates. Int. J. Solids Struct. 1997, 34, 3355–3371. [Google Scholar] [CrossRef]
  3. Caruso, G.; Galeani, S.; Menini, L. Active vibration control of an elastic plate using multiple piezoelectric sensors and actuators. Simul. Model. Pract. Theory 2003, 11, 403–419. [Google Scholar] [CrossRef]
  4. Kumar, K.; Narayanan, S. The optimal location of piezoelectric actuators and sensors for vibration control of plates. Smart Mater. Struct. 2007, 16, 2680–2691. [Google Scholar] [CrossRef]
  5. Dong, X.; Meng, G.; Peng, J. Vibration control of piezoelectric smart structures based on system identification technique: Numerical simulation and experimental study. J. Sound Vib. 2006, 297, 680–693. [Google Scholar] [CrossRef]
  6. Xu, S.; Koko, T. Finite element analysis and design of actively controlled piezoelectric smart structure. Finite Elem. Anal. Des. 2004, 40, 241–262. [Google Scholar] [CrossRef]
  7. Kim, H.; Sohn, J.; Choi, S. Vibration control of a cylindrical shell structure using macro fibre composite actuators. Mech. Based Des. Struct. Mach. 2011, 39, 491–506. [Google Scholar] [CrossRef]
  8. Kusculuoglu, Z.; Royston, T. Finite element formulation for composite plates with piezoceramic layers for optimal vibration control applications. Smart Mater. Struct. 2005, 14, 1139–1153. [Google Scholar] [CrossRef]
  9. Tzou, H.; Bao, Y. Nonlinear piezothermoelasticity and multi-field actuations. J. Vib. Acoust. 1997, 119, 374–381. [Google Scholar] [CrossRef]
  10. Sanbi, M.; Saadani, R.; Sbai, K.; Rahmoune, M. Thermoelastic and Pyroelectric Couplings Effects on Dynamics and Active Control of Smart Piezolaminated Beam Modeled by Finite Element Method. Smart Mater. Res. 2014, 2014, 145087.1–145087.10. [Google Scholar] [CrossRef]
  11. Sanbi, M.; Saadani, R.; Sbai, K.; Rahmoune, M. Thermal Effects on Vibration and Control of Piezocomposite Kirchhoff Plate Modeled by Finite Elements Method. Smart Mater. Res. 2015, 2015, 1–15. [Google Scholar] [CrossRef]
  12. Mitchell, J.; Reddy, H. A refined hybrid platetheoryfor composite laminates with piezoelectric laminae. Int. J. Solids Struct. 1995, 32, 2345–2367. [Google Scholar] [CrossRef]
  13. Gu, H.; Chattopadhyay, A.; Li, J.; Zhou, X. A higher order temperature theory for coupled thermopiezoelectric-mechanical modeling of smart composites. Int. J. Solids Struct. 2000, 37, 6479–6497. [Google Scholar] [CrossRef]
  14. Hughes, P.; Skelton, R. Modal truncation for flexible spacecraft. J. Guid. Control 1981, 4, 291–297. [Google Scholar] [CrossRef]
  15. Willcox, K.; Peraire, J. Balanced model reduction via the proper orthogonal decomposition. AIAA J. 2002, 40, 2323–2330. [Google Scholar] [CrossRef]
  16. Hsieh, W. Nonlinear principal component analysis by neural networks. Tellus A 2002, 53, 2323–2330. [Google Scholar] [CrossRef]
  17. Callahan, J.; Avitabile, P.; Riemer, R. System equivalent reduction expansion process (SEREP). In Proceedings of the 7th International Modal Analysis Conference, Las Vegas, NV, USA, 30 January–2 February 1989; Union College: Schenectady, NY, USA, 1989; pp. 29–37. [Google Scholar]
  18. Sharma, M.; Singh, S.; Sachdeva, B. Modal control of a plate using fuzzy logic controller. Smart Mater. Struct. 2007, 16, 1331–1341. [Google Scholar] [CrossRef]
  19. Tanaka, N.; Sanada, T. Modal control of a rectangular plate using smart sensors and smart actuators. Smart Mater. Struct. 2007, 16, 36. [Google Scholar] [CrossRef]
  20. Rader, A.; Afagh, F.; Yousefi-Koma, A.; Zimcik, D. Optimization of piezoelectric actuator configuration on a flexible fin for vibration control using genetic algorithms. J. Intell. Mater. Syst. Struct. 2007, 18, 1015. [Google Scholar] [CrossRef]
  21. Spier, C.; Bruch, J.; Sloss, J.; Adali, S.; Sadek, I. Placement of multiple piezo patch sensors and actuators for a cantilever beam to maximize frequencies and frequency gaps. J. Vib. Control 2009, 15, 643. [Google Scholar] [CrossRef]
  22. Sharma, S.; Vig, R.; Kumar, N. Temperature compensation in a smart structure by the application of DC bias on piezoelectric patches. J. Intell. Mater. Struct. 2016, 27, 2524–2535. [Google Scholar] [CrossRef]
  23. Gupta, V.; Sharma, M.; Thakur, N. Active structural vibration control: Robust to temperature variations. Mech. Syst. Signal Process. 2012, 33, 80. [Google Scholar] [CrossRef]
  24. Nowacki, W. Some general theorems of thermopiezoelectricity. J. Therm. Stress. 1978, 1, 171–182. [Google Scholar] [CrossRef]
  25. Lal, H.; Sarkar, S.; Gupta, S. Stochastic model order reduction in randomly parametered linear dynamical systems. Appl. Math. Model. 2017, 51, 744–763. [Google Scholar] [CrossRef]
  26. Anderson, B.; Moore, J. Optimal Filtering; Thomas Kailath Editor: Piscataway, NJ, USA, 1979. [Google Scholar]
  27. Simon, D.; Curadelli, O.; Ambrosini, D. Kalman Filtering for fuzzy discrete time dynamic systems. Appl. Soft Comput. 2003, 3, 191–207. [Google Scholar] [CrossRef]
  28. Garrido, H.; Curadelli, O.; Ambrosini, D. A straightforward method for tuning of Lyapunov based controllers in semi active vibration control applications. J. Sound Vib. 2014, 333, 1119–1131. [Google Scholar] [CrossRef]
Figure 1. The smart plate with distributed thermo-piezoelectric patches and the four-noded Kirchhoff element with mechanical degrees of freedom.
Figure 1. The smart plate with distributed thermo-piezoelectric patches and the four-noded Kirchhoff element with mechanical degrees of freedom.
Jcs 06 00242 g001
Figure 2. A schematic representation corresponding to the reduced model control LQR.
Figure 2. A schematic representation corresponding to the reduced model control LQR.
Jcs 06 00242 g002
Figure 3. A schematic representation corresponding to the reduced model of LQG-Kalman control.
Figure 3. A schematic representation corresponding to the reduced model of LQG-Kalman control.
Jcs 06 00242 g003
Figure 4. 3D first mode shapes 1, 2, 3 and 4 of the plate vibration.
Figure 4. 3D first mode shapes 1, 2, 3 and 4 of the plate vibration.
Jcs 06 00242 g004
Figure 5. 3D mode shapes 5, 6, 7 and 8 of the plate vibration.
Figure 5. 3D mode shapes 5, 6, 7 and 8 of the plate vibration.
Jcs 06 00242 g005
Figure 6. Impulse responses of control using reduced LQG-Kalman model for various temperatures with 5 distributed pairs of sensor/actuator and 10 selected sensors nodes.
Figure 6. Impulse responses of control using reduced LQG-Kalman model for various temperatures with 5 distributed pairs of sensor/actuator and 10 selected sensors nodes.
Jcs 06 00242 g006
Figure 7. Comparison of impulse responses of LQG-Kalman and LQR controls in reduced model for T = 80 ° C with 4 distributed pairs of sensor/actuator and 8 selected sensors nodes.
Figure 7. Comparison of impulse responses of LQG-Kalman and LQR controls in reduced model for T = 80 ° C with 4 distributed pairs of sensor/actuator and 8 selected sensors nodes.
Jcs 06 00242 g007
Figure 8. Impulse responses of control corresponding to the full and reduced LQG-Kalman models for T = 25 ° C with 8 distributed pairs of sensor/actuator and 20 selected sensors nodes.
Figure 8. Impulse responses of control corresponding to the full and reduced LQG-Kalman models for T = 25 ° C with 8 distributed pairs of sensor/actuator and 20 selected sensors nodes.
Jcs 06 00242 g008
Figure 9. Actuator input Voltage using full and reduced LQG-Kalman models for T = 25 ° C with 8 distributed pairs of sensor/actuator and 20 selected sensors nodes.
Figure 9. Actuator input Voltage using full and reduced LQG-Kalman models for T = 25 ° C with 8 distributed pairs of sensor/actuator and 20 selected sensors nodes.
Jcs 06 00242 g009
Figure 10. Sensor output voltage using LQG-Kalman control for full model for multiple temperatures with 5 distributed pairs of sensor/actuator and 10 selected sensor nodes.
Figure 10. Sensor output voltage using LQG-Kalman control for full model for multiple temperatures with 5 distributed pairs of sensor/actuator and 10 selected sensor nodes.
Jcs 06 00242 g010
Figure 11. Sensor output voltage using LQG-Kalman control in reduced model for multiple temperatures with 5 distributed pairs of sensor/actuator and 10 selected sensor nodes.
Figure 11. Sensor output voltage using LQG-Kalman control in reduced model for multiple temperatures with 5 distributed pairs of sensor/actuator and 10 selected sensor nodes.
Jcs 06 00242 g011
Figure 12. Impulse responses Corresponding to the LQG-Kalman control for different selections of dofs and 10 distributed sensors along the surface of the plate.
Figure 12. Impulse responses Corresponding to the LQG-Kalman control for different selections of dofs and 10 distributed sensors along the surface of the plate.
Jcs 06 00242 g012
Table 1. Material properties of the piezoelectric elements and the host plate.
Table 1. Material properties of the piezoelectric elements and the host plate.
PropertiesActuators
PZT
Sensors
PVDF
Host
Plate
Length (cm) l a l s 18
Width (cm) b a b s 18
Thickness (cm) 0 , 6 0 , 5 0 , 08
Piezoelectric constant (pm/V)
d 31 , d 32 , d 33 , d 24 , d 15
−220, −220,
374, 670, 670
−22, 22, 0, 0, 0
Permittivity (nF/m)
ϵ 11 = ϵ 22 ϵ 33
15, 30, 1062
Density (Kg/m 3 )760027001800
Young’s modulus E (GPa)63 2.0 70
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

EI Khaldi, L.; Sanbi, M.; Saadani, R.; Rahmoune, M. Robust Control and Thermal Analysis of a Reduced Model of Kirchhoff Composite Plate with Random Distribution of Thermopiezoelectric Sensors and Actuators. J. Compos. Sci. 2022, 6, 242. https://doi.org/10.3390/jcs6080242

AMA Style

EI Khaldi L, Sanbi M, Saadani R, Rahmoune M. Robust Control and Thermal Analysis of a Reduced Model of Kirchhoff Composite Plate with Random Distribution of Thermopiezoelectric Sensors and Actuators. Journal of Composites Science. 2022; 6(8):242. https://doi.org/10.3390/jcs6080242

Chicago/Turabian Style

EI Khaldi, Loukmane, Mustapha Sanbi, Rachid Saadani, and Miloud Rahmoune. 2022. "Robust Control and Thermal Analysis of a Reduced Model of Kirchhoff Composite Plate with Random Distribution of Thermopiezoelectric Sensors and Actuators" Journal of Composites Science 6, no. 8: 242. https://doi.org/10.3390/jcs6080242

APA Style

EI Khaldi, L., Sanbi, M., Saadani, R., & Rahmoune, M. (2022). Robust Control and Thermal Analysis of a Reduced Model of Kirchhoff Composite Plate with Random Distribution of Thermopiezoelectric Sensors and Actuators. Journal of Composites Science, 6(8), 242. https://doi.org/10.3390/jcs6080242

Article Metrics

Back to TopTop