Optimal Shape Control of Piezoelectric Intelligent Structure Based on Genetic Algorithm

Shape variation induced by mismachining tolerance, humidity and temperature of the working environment, material wear and aging, and unknown external load disturbances have a relatively large influence on the dynamic shape of a mechanical structure. When integrating piezoelectric elements into the main mechanical structure, active control of the structural shape is realized by utilizing the inverse piezoelectric effect. This paper presents a mathematical model regarding piezoelectric intelligent structure shape control. We also applied a genetic algorithm, and given a piezoelectric intelligent cantilever plate with both ends affected by a certain load, optimal shape control results of piezoelectric materials were analyzed from different perspectives (precision reference or cost reference). The mathematical model and results indicate that, by optimizing a certain number of piezoelectric actuators, high-precision active shape control can be realized.


Introduction
A variety of high-speed and high-precision mechanical structures require extremely precise shape and position when in operation, and slight variations greatly influence the dynamic performances of these structures [1,2].For instance, mismachining tolerance, humidity and temperature of the working environment, material wear and aging, and unknown external load disturbances may affect the shape of a structure.Consequently, active shape control and compensation are required in real operation.Piezoelectric elements are integrated into the main structure, the inverse piezoelectric effect is utilized, and piezoelectric elements are deformed via an external voltage, which is passed on to the main structure, thus realizing active control of structural shape.
Due to structure complexity and the inherent electromechanical coupling effect of a piezoelectric intelligent structure, studies on structural shape control are typically based on the mathematical mode of the intelligent structure using a finite element method.Following conclusions toward the theory of intelligent structure, Wada et al. [3] proposed a shape-variable and self-adaptive intelligent structure framework to meet future requirements of structures like space stations, aircrafts, and satellites.Based on a hierarchical theory, Donthireddy and Chandrashekhara [4] established a finite element numerical model of laminated piezoelectric beams, studied the shape control problem of these beam, and analyzed the influence of stacking sequence and boundary conditions on shape control.Varadarajan et al. [5] performed shape control of laminated composite plates with piezoelectric actuators, and using the minimum error function of ideal shape and control shape, voltage distribution was optimized and an active feedback control algorithm was established.Conversely, Chee et al. [6,7] established the static shape control equation of intelligent composite piezoelectric beams and plates based on a high-order displacement function, and classic beam and plate theory, respectively.Furthermore, they applied the generalized function to analyze the specific control effect.Lastly, Lin and Nien [8] established the finite element model of laminated plate shape control using piezoelectric actuators, considered the influence of actuator and sensor locations on shape control, calculated structural stress, and concluded that internal stress of the piezoelectric actuator has huge influence on shape control.
It is important to optimize the design of embedded depth, thickness, and location of piezoelectric elements, as well as 2 Advances in Materials Science and Engineering control voltage, for realizing the shape control effect.Based on a nonlinear constitutive piezoelectric equation and according to the error function and control energy, Sun and Tong [9] utilized the finite element method and Lagrange multiplier method to define an optimization algorithm for piezoelectric actuators on the control voltage of static nonlinear deformation control, and further demonstrated the effectiveness of this algorithm through simulation experiments.On the other hand, Barboni et al. [10] used a method combining a dynamic influence function and closed-loop feedback to study the optimal location of a pair of piezoelectric actuators, thus maximizing piezoelectric intelligent beam displacement.Another study considered the influence of the location, dimension, and voltage of piezoelectric actuators on shape control, and performed multi-objective optimization of cantilever shape control, minimizing beam deflection under the external load effect [11].Conversely, Zhang et al. [12] studied the integrated optimization configuration problem of piezoelectric actuators and sensors in a flexible structure system by applying the genetic algorithm, while Da Mota Silva et al. [13] used the variance between the preset and real displacement of a specific node to study the static shape control of self-adaptable structure, established an infinite model of the structure system, and determined the optimal driving voltage of piezoelectric elements by applying the genetic algorithm.Experiments verified the applicability and effectiveness of this algorithm.In contrast, based on the shear deformation beam theory and linear piezoelectric theory, Hadjigeorgiou et al. [14] established a numerical model of piezoelectric beams, and utilized the genetic optimization algorithm to optimize the voltage of piezoelectric actuator elements.Numerical results indicated that fewer piezoelectric driving elements achieved the same expected deformation effect following optimization.
In this paper, the piezoelectric intelligent structure was applied to perform active shape control, and the effect location, geometric parameters, and working voltage were solved based on the expected shape or displacement value.Thus, this is a typical multiobjective optimization problem.First, according to the shape function of the piezoelectric intelligent structure and the finite dynamic function, the relationship between the shape or displacement of the structure and its piezoelectric properties, function location, geometric parameters, or input voltage under the influence of constant mechanical load was deduced.Then, from different perspectives (precision preference or cost preference), proper object function was selected to perform multiobjective optimization of the number, function location, and working voltage of the piezoelectric driving elements, thus providing optimal shape control results that satisfy precision requirements.

Piezoelectric Intelligent Structure
Shape Control Theory where The relationship between natural coordinates and geometric coordinates is where  and  are the length and width of rectangular element, respectively.The shape function is expressed according to After substituting (2), (3), and (4) into (1), the displacement of any point inside the finite element is expressed as follows: where   = { 1 ,  1 ,  1 ,  1 , . . .,  4 ,  4 ,  4 ,  4 }, and   (, ) represents the shape functions of a mode in every degree of freedom, which are listed below: The linear constitutive equations of the piezoelectric material can be written as [15] where c E is the elastic stiffness matrix, the superscript  denotes transpose, E is the electric field vector, D is the electric displacement vector, M S is the dielectric constant matrix, and  is the piezo stress/charge constant.
The finite element dynamic equation of the piezoelectric intelligent structure is expressed as follows: where q and  are the displacement vector and potential vector of each mode in the overall structure, respectively; M, C, and K  are the overall weight, resistance, and stiffness matrix of the integrated piezoelectric intelligent structure, respectively; K  = K   is the overall force-electricity coupled stiffness matrix of the integrated piezoelectric material; K  is the overall dielectric stiffness matrix of the piezoelectric material; and F and Q are the structure's external force load vector and electric load vector, respectively.The detailed derivation of these element matrices can be found in [16].
Taking into consideration the active static shape control of the piezoelectric intelligent structure, the dynamic response term in the equation is neglected, and (8) then becomes (9).This can then be rewritten as (10):

Shape Control Mechanism of the Intelligent
Structure.The static condensation method was applied to simplify (10) to obtain where From ( 11), the external load of the piezoelectric intelligent structure system consists of 2 parts: the mechanical load F and electric load F  .The electric load is determined by the dielectric properties, geometric parameters, and input voltage of the piezoelectric driving element.Under the condition that the mechanical load remains constant, because of the inverse piezoelectric effect, the displacement of the structure under control can be changed by changing the driving element input voltage changes.
The external electric field of the piezoelectric intelligent structure's sensing unit is zero, and therefore, the potential output is obtained by the second term in (10): where −K −1  K  is defined as the displacement sensitivity matrix.
Expressions for stiffness matrix K   , force-electric coupling stiffness matrix K   , and unit dielectric matrix K   in the finite unit of the piezoelectric intelligent structure are as follows: where matrices [  ] and [  ] consist of partial derivatives of the displacement shape function f  and potential shape function f  of the unit along the displacement and potential directions, respectively, According to the boundary conditions, ( 14) are integrated to obtain the overall stiffness matrix K  , force-electric coupling stiffness matrix K  , and dielectric matrix K  .From ( 14), the force-electric coupling stiffness matrix K  and dielectric matrix K  of the system are not only influenced by parameters of the material itself, but also closely related to the location and number of piezoelectric driving units.Therefore, by changing the location, number, and input voltage of piezoelectric driving units, the overall displacement q in (11) may be changed to achieve shape control.A basic piezoelectric intelligent structure consists of a main structural layer, an upper layer, and a lower layer.If 1 piezoelectric layer is used as the sensing unit, the other layer is used as the driving unit, the output potential of the sensing unit is enlarged using the feedback control system, and the output voltage required by the driving unit is calculated via the feedback gain and control law.Then, active static shape control of the piezoelectric intelligent structure is fulfilled.The static shape system of a simple piezoelectric intelligent structure is shown in Figure 1.

Specific Configuration of the Active Shape Controlled Genetic Algorithm
The location and working voltage of piezoelectric actuators and sensors are the optimization goal of static shape control in a piezoelectric intelligent structure.Based on the finite element model of the overall structure, the finite elements shown in Figure 1 are discretized.Assume the piezoelectric actuators/sensors in the upper and lower surfaces of the finite elements in a unit are independently controllable, with the minimum variance of the actual and expected displacement of the structure being the target function, and with the location, number, and driving voltage of the piezoelectric plate being the variables of optimal design.Then, a computational program for the genetic optimization algorithm was compiled to seek the optimal solution that satisfies the expected shape of the piezoelectric intelligent structure.
Among the design variables, the location and number of piezoelectric plates are binary variables, while the driving voltage is a continuous variable.In this paper, a parameter cascade method was applied to code the location and driving voltage of the piezoelectric plate, and further compile codes to form the individual coding that represents the entire parameters.
Assume variables X 1 and X 2 represent the location and driving voltage of the piezoelectric plate, respectively.The

Discrete variable
Binary coding

Continuous variable
Float coding binary-encoded method and the float-encoded method were applied to encode X 1 and X 2 , respectively, and then where X 1 is a discrete variable.Using the binary-encoded method,  11 = 0 represents a situation where the unit finite element is not equipped with actuators/sensors, and  11 = 1 represents a situation where the unit finite element is equipped with actuators/sensors.In addition,  represents the number of finite elements in the piezoelectric intelligent structure plate.Assuming the intelligent structure is equipped with  pairs of actuators/sensors, then Here, X 2 is a continuous variable, float-encoded variable that directly reflects the actual value of the design variable and has the same length as variable X 1 .In addition,  2 is any real number within the given range of driving voltage: [−, ]. Figure 2 illustrates the individual string structure that applies the parameter cascade code.
The genetic strategy was determined according to the form of the individual parameter cascade code, and the optimization model of piezoelectric intelligent structure shape control was established.The optimization algorithm system expression is as follows: where (X 1 , X 2 ) is the target function of the genetic algorithm optimization, (X 1 , X 2 ) is the optimized restriction function, and () is the individual fitness function.
In real optimization of static shape control with different optimization objects, the corresponding optimization configuration differs.Taking into consideration the economy of controlling the cost, number, and voltage of actuators/sensors requires optimization.A minimum number of piezoelectric actuators/sensors and minimum driving voltage were used to obtain an ideal control shape.The genetic algorithm configuration is then described using a bi-objective optimization model: where  is the number of structure units;  is the number of structure nodes; V  is the optimal displacement of each node; V  is the targeted displacement for each node; and [−, ] is the space defined by the driving voltage.Furthermore, the target function of the genetic optimization algorithm  1 (X 1 , X 2 ) represents the configuration number of piezoelectric plates,  2 (X 1 , X 2 ) is the sum of absolute driving voltages of each plate, the restriction function (X 1 , X 2 ) provides variance between the optimal and targeted shapes that satisfies control precision   , and the fitness functions  1 ( 1 ) and  2 ( 2 ) are the differences between target functions  1 , and  2 and their theoretical maximum values, respectively.Another optimization configuration method uses a given number of piezoelectric actuators/sensors, and optimizing the configuration location and driving voltage of the piezoelectric plate to achieve the optimal shape control effect within the precision requirement range.The specific configuration method of the genetic algorithm is as follows: where  is the given number of piezoelectric actuators/sensors, while the other variables are the same as those in (19).Under this configuration, the target function of the genetic algorithm  1 (X 1 , X 2 ) becomes the variance of the optimal shape and target shape, the restriction function (X 1 , X 2 ) becomes the number of piezoelectric plates equal to the given value, and the fitness function  1 ( 1 ) becomes the reciprocal of 1 +  1 (X 1 , X 2 ).
Referring to the specific operational procedure of the genetic algorithm and based on the finite element model of piezoelectric intelligent structure, a computational procedure was compiled using Matlab, and a series of genetic algorithm optimization procedures were designed to realize active control optimization of the piezoelectric intelligent structure.

Numerical Analysis
Given a static shape control of a cantilever plate with the upper and lower surfaces painted with piezoelectric  intelligent layers, the control theory of the finite shape and optimization control of the genetic algorithm were verified.The finite element model of the structure, which was divided into 32 finite rectangular plate elements, is shown in Figure 3.The upper and lower surfaces of the structure were symmetrically attached with piezoelectric thin film PZT-5H (1 mm thickness).The main structural material in the middle layer was Al, which had a thickness of 2.5 mm.The property parameters of each material in the piezoelectric intelligent structure are listed in Table 1.Lastly, the middle point of the structural free end was affected by a constant load , and the shape control target was defined as zero displacement of the free end.Structure optimization was then performed.function.The shape control precision requirement was set as less than 1 × 10 −3 mm 2 , the maximum number of generation as 1000, the crossover possibility of the genetic algorithm as   = 0.625, and the mutation probability as   = 1/64.Then, the optimization calculation of the genetic algorithm was performed.

Configuration Analysis of the Genetic
Firstly, we give random 32 binary-encoded discrete design parameters as the initial value.The displacement of the free end can be calculated through (10), and the restriction function (X 1 , X 2 ) provides variance of actual displacement and expected displacement at the free end.If the shape control precision requirement is more than 1 × 10 −3 mm 2 , mutation will be applied to change the design parameters G1 and G2 to the new generation according to the mutation probability.The new generation design parameters will be used to start the new iteration.
After 1000 iterations, the optimization configuration that satisfied the shape control precision requirement was obtained.Figure 4 shows the performance tracking process of the restriction function value.The restriction function value rapidly converged and approached the required precision value after 350 iterations.After 1000 iterations, precision reached 9 × 10 −4 mm 2 .Figures 5 and 6 show the iterative evolution of the target function value and the fitness function value, respectively.Similarly, in the early stage of genetic optimization, the target function value immediately converged, and after 645 iterations, the number of configuration design parameters evolved from 32 to 6, less than one-fifth of the initial design, which greatly saved the control cost.A configuration information diagram of each piezoelectric actuator  in the iterative evolution process is shown in Figure 7, and Table 2 shows the working voltage of each piezoelectric actuator under the best configuration.It should be noted that variance of actual displacement and desired displacement at the free end was set as the restriction function, and then the symmetry of numerical results along the central axis of the beam could not be expected in genetic optimization.For example, the values of voltages between elements 10 and 11, or 17 and 20 are derived to be positive or negative in Figure 7 and Table 2. Figures 8, 9, and 10 show shape diagrams of the piezoelectric intelligent cantilever at initial status, during evolution, and after optimization, respectively.

Configuration Analysis of Genetic Algorithm: Optimizing Shape Control Precision.
With the number of actuators fixed, the optimum values of actuator location and working voltage were obtained via the optimization of the genetic algorithm to achieve the best control effect.The optimized genetic algorithm was performed according to (19), the parameter cascade encoding method was adopted, and 32 binary-encoded discrete design variables were chosen to represent configuration information regarding the location and number of piezoelectric actuators.The 32 float-encoded continuous design variables represented the working voltage within the range of [−100, 100].The shape control precision was selected as the target function (sum of squares of displacement and target displacement of 5 nodes in the free end after optimization), and the shape control precision requirement was set as less than 1 × 10 −5 mm 2 .Given 10 piezoelectric actuators, the restriction function value was 16.Lastly, we assumed the maximum number of generations as 2000, the crossover possibility of the genetic algorithm as   = 0.625, and the mutation probability as   = 1/64.Then, the optimization calculation of the genetic algorithm was performed.
Figure 11 shows the target function over 2000 iterations.The function rapidly converged during the initial iterative stages and approached the required precision after 400 iterations.Figure 12 shows the iterative variance of the corresponding fitness function.For a certain number of actuators, a relatively high control precision can be reached with proper optimization configuration.After 2000 iterations, the shape control precision of the piezoelectric actuators reached 5.21 × 10 −6 mm 2 .Figure 13 shows the configuration diagram of the piezoelectric actuators at different evolution periods during iteration, while Figures 14 and 15 provide shape deformation diagrams of the piezoelectric intelligent plates during the optimization control process.Lastly, Table 3 shows the working voltage of each piezoelectric actuator under optimal configuration following iteration.

Conclusions
(1) Based on the dynamic finite element equation of piezoelectric intelligent structure, the static shape control principle of treating the piezoelectric material as an actuator was analyzed, demonstrating that the shape control effect is not only influenced by the parameters of the material itself, but also closely related to the location and number of piezoelectric driving elements.(2) By utilizing operation characteristics of the genetic algorithm (such as the embarrassingly parallel, random, and self-adaptive), each finite element of the piezoelectric intelligent structure was separated, and the computational and basic operational procedures of the optimal shape control genetic algorithm were designed to find the optimal solution that satisfied the expected shape.(3) Based on the modified parameter cascade encoding method, we used the binary-encoded method and floatencoded method to code the location and driving voltage of the piezoelectric plates, respectively.We proposed a genetic algorithm optimization mode for shape control, and from the perspective of cost control economy and shape control precision, we then designed a computational program for the genetic algorithm.(4) Given a piezoelectric intelligent

Figure 1 :
Figure 1: Active static shape control system of a piezoelectric intelligent plate.

Figure 2 :
Figure 2: Individual bit string structure of parameter cascade coding.

Figure 3 :
Figure 3: Finite element model of cantilever intelligent structural plate.

Figure 4 :
Figure 4: Variation of restriction function value.

Figure 5 :
Figure 5: Variation of target function value.

Figure 6 :
Figure 6: Variation of the fitness function value.

Table 1 :
Property parameters of cantilever intelligent structural plate.

Table 2 :
Optimal working voltage of piezoelectric actuators under optimal configuration.

Table 3 :
Optimal working voltage of piezoelectric actuator.
cantilever with a load on 1 end, static shape optimization control results of piezoelectric materials based on the genetic algorithm were analyzed.Conclusively, our mathematical model and results indicated that active shape control with