# Optimal Shape Control of Piezoelectric Intelligent Structure Based on Genetic Algorithm

Authors: Zhanxi Wang; Xiansheng Qin; Shunqi Zhang; Jing Bai; Jing Li; Genjie Yu
Journal: Advances in Materials Science and Engineering
 (2017)

---

## Abstract

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.


---

## Body

## 1. 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 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.


## 2. Piezoelectric Intelligent Structure Shape Control Theory

### 2.1. Shape Function and Finite Element Equation

The top and bottom layers of the piezoelectric intelligent plate are symmetrically bonded with piezoelectric actuators, which are applied to produce control bending moments and detect shape change. Because the thickness of these piezoelectric actuators is assumed to be very small as compared with the thickness of the plates, the structure changes due to these piezoelectric layers can be considered negligible. The piezoelectric intelligent plate is modeled by using the four-corner-node rectangular plate elements, which include four bending degrees of freedomw,wx,wy,wxy at each corner node. According to the Hermitian interpolation polynomial relationship between dimensionless natural coordinates ξ and η, the relative displacement in any point inside the rectangular element can be expressed as follows: (1)wξ,η=∑i=14GiξGiηwi+HiξGiη∂w∂xi+GiξHiη∂w∂yi+GiξGiη∂w∂zi,where (2)Giξ=14-ξiξ3+3ξiξ+2,Hiξ=14ξ3+ξiξ2-ξ-ξi.The relationship between natural coordinates and geometric coordinates is (3)ξ=2xa-1,η=2yb-1,where a and b are the length and width of rectangular element, respectively.The shape function is expressed according to(4)f1x=G1ξ=G4ξ,f1y=G1η=G2η,f2x=G2ξ=G3ξ,f2y=G3η=G4η,f3x=H1ξ=H4ξ,f3y=H1η=H2η,f4x=H2ξ=H3ξ,f4y=H3η=H4η.After substituting (2), (3), and (4) into (1), the displacement of any point inside the finite element is expressed as follows: (5)wx,y=∑i=116fix,yqi,where qi=w1,w1x,w1y,w1xy,…,w4,w4x,w4y,w4xy, and fi(x,y) represents the shape functions of a mode in every degree of freedom, which are listed below:(6)f1x,y=f1xf1y=1a3b3a3+2x3-3ax2b3+2y3-3by2,f2x,y=f2xf1y=1a3b3ax3-2a2x2+a3xb3+2y3-3by2,f3x,y=f1xf2y=1a3b3a3+2x3-3ax2by3-2b2y2+b3y,f4x,y=f2xf2y=1a3b3ax3-2a2x2+a3xby3-2b2y2+b3y,f5x,y=f3xf1y=1a3b33ax2-2x3b3+2y3-3by2,f6x,y=f4xf1y=1a3b3ax3-a2x2b3+2y3-3by2,f7x,y=f3xf2y=1a3b33ax2-2x3by3-2b2y2+b3y,f8x,y=f4xf2y=1a3b3ax3-a2x2by3-2b2y2+b3y,f9x,y=f3xf3y=1a3b33ax2-2x33by2-2y3,f10x,y=f4xf3y=1a3b3ax3-a2x23by2-2y3,f11x,y=f3xf4y=1a3b33ax2-2x3by3-b2y2,f12x,y=f4xf4y=1a3b3ax3-a2x2by3-b2y2,f13x,y=f1xf3y=1a3b3a3+2x3-3ax23by2-2y3,f14x,y=f2xf3y=1a3b3ax3-2a2x2+a3x3by2-2y3,f15x,y=f1xf4y=1a3b3a3+2x3-3ax2by3-b2y2,f16x,y=f2xf4y=1a3b3ax3-2a2x2+a3xby3-b2y2.The linear constitutive equations of the piezoelectric material can be written as [15](7)σ=cEε-eTE,D=eε+єSE,where cE is the elastic stiffness matrix, the superscript T denotes transpose, E is the electric field vector, D is the electric displacement vector, єS is the dielectric constant matrix, and e is the piezo stress/charge constant.The finite element dynamic equation of the piezoelectric intelligent structure is expressed as follows:(8)M000q¨ϕ¨+C000q˙ϕ˙+KqqKqϕKϕqKϕϕqϕ=FQ,where q and ϕ are the displacement vector and potential vector of each mode in the overall structure, respectively; M, C, and Kqq are the overall weight, resistance, and stiffness matrix of the integrated piezoelectric intelligent structure, respectively; Kqϕ=KϕqT 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):(9)KqqKqϕKϕqKϕϕqϕ=FQ,(10)Kqqq+Kqϕϕ=F,Kϕqq+Kϕϕϕ=Q.


### 2.2. Shape Control Mechanism of the Intelligent Structure

The static condensation method was applied to simplify (10) to obtain (11)Kq=F+Fϕ, where (12)K=Kqq-KqϕKϕϕ-1Kϕq,Fϕ=-KqϕKϕϕ-1Q.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): (13)ϕ=-Kϕϕ-1Kϕqq,where -Kϕϕ-1Kϕq is defined as the displacement sensitivity matrix.Expressions for stiffness matrixKqqe, force-electric coupling stiffness matrix Kqϕe, and unit dielectric matrix Kϕϕe in the finite unit of the piezoelectric intelligent structure are as follows: (14)Kqqe=∫VBqTcEBqdV,Kqϕe=∫VBqTeTBqdV,Kϕϕe=∫VBϕTєSBϕdV,where matrices Bq and Bϕ consist of partial derivatives of the displacement shape function fq and potential shape function fϕ of the unit along the displacement and potential directions, respectively,(15)Bq=∂2fq∂x2∂2fq∂y2∂2fq∂z2T,Bϕ=∂fϕ∂x∂fϕ∂y∂fϕ∂zT.According to the boundary conditions, (14) are integrated to obtain the overall stiffness matrix Kqq, force-electric coupling stiffness matrix Kqϕ, and dielectric matrix Kϕϕ. From (14), the force-electric coupling stiffness matrix Kqϕ 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 displacementq 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.Figure 1
Active static shape control system of a piezoelectric intelligent plate.


## 2.1. Shape Function and Finite Element Equation

The top and bottom layers of the piezoelectric intelligent plate are symmetrically bonded with piezoelectric actuators, which are applied to produce control bending moments and detect shape change. Because the thickness of these piezoelectric actuators is assumed to be very small as compared with the thickness of the plates, the structure changes due to these piezoelectric layers can be considered negligible. The piezoelectric intelligent plate is modeled by using the four-corner-node rectangular plate elements, which include four bending degrees of freedomw,wx,wy,wxy at each corner node. According to the Hermitian interpolation polynomial relationship between dimensionless natural coordinates ξ and η, the relative displacement in any point inside the rectangular element can be expressed as follows: (1)wξ,η=∑i=14GiξGiηwi+HiξGiη∂w∂xi+GiξHiη∂w∂yi+GiξGiη∂w∂zi,where (2)Giξ=14-ξiξ3+3ξiξ+2,Hiξ=14ξ3+ξiξ2-ξ-ξi.The relationship between natural coordinates and geometric coordinates is (3)ξ=2xa-1,η=2yb-1,where a and b are the length and width of rectangular element, respectively.The shape function is expressed according to(4)f1x=G1ξ=G4ξ,f1y=G1η=G2η,f2x=G2ξ=G3ξ,f2y=G3η=G4η,f3x=H1ξ=H4ξ,f3y=H1η=H2η,f4x=H2ξ=H3ξ,f4y=H3η=H4η.After substituting (2), (3), and (4) into (1), the displacement of any point inside the finite element is expressed as follows: (5)wx,y=∑i=116fix,yqi,where qi=w1,w1x,w1y,w1xy,…,w4,w4x,w4y,w4xy, and fi(x,y) represents the shape functions of a mode in every degree of freedom, which are listed below:(6)f1x,y=f1xf1y=1a3b3a3+2x3-3ax2b3+2y3-3by2,f2x,y=f2xf1y=1a3b3ax3-2a2x2+a3xb3+2y3-3by2,f3x,y=f1xf2y=1a3b3a3+2x3-3ax2by3-2b2y2+b3y,f4x,y=f2xf2y=1a3b3ax3-2a2x2+a3xby3-2b2y2+b3y,f5x,y=f3xf1y=1a3b33ax2-2x3b3+2y3-3by2,f6x,y=f4xf1y=1a3b3ax3-a2x2b3+2y3-3by2,f7x,y=f3xf2y=1a3b33ax2-2x3by3-2b2y2+b3y,f8x,y=f4xf2y=1a3b3ax3-a2x2by3-2b2y2+b3y,f9x,y=f3xf3y=1a3b33ax2-2x33by2-2y3,f10x,y=f4xf3y=1a3b3ax3-a2x23by2-2y3,f11x,y=f3xf4y=1a3b33ax2-2x3by3-b2y2,f12x,y=f4xf4y=1a3b3ax3-a2x2by3-b2y2,f13x,y=f1xf3y=1a3b3a3+2x3-3ax23by2-2y3,f14x,y=f2xf3y=1a3b3ax3-2a2x2+a3x3by2-2y3,f15x,y=f1xf4y=1a3b3a3+2x3-3ax2by3-b2y2,f16x,y=f2xf4y=1a3b3ax3-2a2x2+a3xby3-b2y2.The linear constitutive equations of the piezoelectric material can be written as [15](7)σ=cEε-eTE,D=eε+єSE,where cE is the elastic stiffness matrix, the superscript T denotes transpose, E is the electric field vector, D is the electric displacement vector, єS is the dielectric constant matrix, and e is the piezo stress/charge constant.The finite element dynamic equation of the piezoelectric intelligent structure is expressed as follows:(8)M000q¨ϕ¨+C000q˙ϕ˙+KqqKqϕKϕqKϕϕqϕ=FQ,where q and ϕ are the displacement vector and potential vector of each mode in the overall structure, respectively; M, C, and Kqq are the overall weight, resistance, and stiffness matrix of the integrated piezoelectric intelligent structure, respectively; Kqϕ=KϕqT 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):(9)KqqKqϕKϕqKϕϕqϕ=FQ,(10)Kqqq+Kqϕϕ=F,Kϕqq+Kϕϕϕ=Q.


## 2.2. Shape Control Mechanism of the Intelligent Structure

The static condensation method was applied to simplify (10) to obtain (11)Kq=F+Fϕ, where (12)K=Kqq-KqϕKϕϕ-1Kϕq,Fϕ=-KqϕKϕϕ-1Q.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): (13)ϕ=-Kϕϕ-1Kϕqq,where -Kϕϕ-1Kϕq is defined as the displacement sensitivity matrix.Expressions for stiffness matrixKqqe, force-electric coupling stiffness matrix Kqϕe, and unit dielectric matrix Kϕϕe in the finite unit of the piezoelectric intelligent structure are as follows: (14)Kqqe=∫VBqTcEBqdV,Kqϕe=∫VBqTeTBqdV,Kϕϕe=∫VBϕTєSBϕdV,where matrices Bq and Bϕ consist of partial derivatives of the displacement shape function fq and potential shape function fϕ of the unit along the displacement and potential directions, respectively,(15)Bq=∂2fq∂x2∂2fq∂y2∂2fq∂z2T,Bϕ=∂fϕ∂x∂fϕ∂y∂fϕ∂zT.According to the boundary conditions, (14) are integrated to obtain the overall stiffness matrix Kqq, force-electric coupling stiffness matrix Kqϕ, and dielectric matrix Kϕϕ. From (14), the force-electric coupling stiffness matrix Kqϕ 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 displacementq 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.Figure 1
Active static shape control system of a piezoelectric intelligent plate.


## 3. 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 Figure1 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 variablesX1 and X2 represent the location and driving voltage of the piezoelectric plate, respectively. The binary-encoded method and the float-encoded method were applied to encode X1 and X2, respectively, and then(16)X1=γ11,γ12,γ13,…,γ1n,X2=γ21,γ22,γ23,…,γ2n,where X1 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, n represents the number of finite elements in the piezoelectric intelligent structure plate. Assuming the intelligent structure is equipped with l pairs of actuators/sensors, then (17)∑i=1nγ1i=l.Here,X2 is a continuous variable, float-encoded variable that directly reflects the actual value of the design variable and has the same length as variable X1. In addition, γ2i is any real number within the given range of driving voltage: -a,a. Figure 2 illustrates the individual string structure that applies the parameter cascade code.Figure 2
Individual bit string structure of parameter cascade coding.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:(18)minGX1,X2S.t.RX1,X2≤QFit.FG,where GX1,X2 is the target function of the genetic algorithm optimization, RX1,X2 is the optimized restriction function, and FG 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:(19)minG1X1,X2=∑i=1nX1iminG2X1,X2=∑i=1nX2iS.t.RX1,X2=∑j=1mvOj-vGj2≤RpFit.F1G1=n-G1X1,X2Fit.F2G2=na-G2X1,X2,where n is the number of structure units; m is the number of structure nodes; vOj is the optimal displacement of each node; vGj is the targeted displacement for each node; and [-a,a] is the space defined by the driving voltage. Furthermore, the target function of the genetic optimization algorithm G1X1,X2 represents the configuration number of piezoelectric plates, G2X1,X2 is the sum of absolute driving voltages of each plate, the restriction function RX1,X2 provides variance between the optimal and targeted shapes that satisfies control precision Rp, and the fitness functions F1G1 and F2G2 are the differences between target functions G1, and G2 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:(20)minG1X1,X2=∑j=1mvOj-vGj2≤RpminG2X1,X2=∑i=1nX2iS.t.RX1,X2=∑i=1nX1i=lFit.F1G1=11+G1X1,X2Fit.F2G2=na-G2X1,X2,where l 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 G1X1,X2 becomes the variance of the optimal shape and target shape, the restriction function RX1,X2 becomes the number of piezoelectric plates equal to the given value, and the fitness function F1G1 becomes the reciprocal of 1+G1X1,X2.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.


## 4. 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 Figure3. 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 P, and the shape control target was defined as zero displacement of the free end. Structure optimization was then performed.Table 1
Property parameters of cantilever intelligent structural plate.
Main structure: Al
Elastic module:E = 71.0 GPa; the Poison parameter: v=0.3; density: ρ = 2690 kg/m3
Intelligent material: PZT-5H
Elastic parameter:E = 126.0 GPa; Poison parameter: v = 0.35; density: ρ = 7500 kg/m3
Elastic rigid parameter:c11 = 12.6 × 1010 N/m2, c12 = 7.95 × 1010 N/m2, c13 = 8.41 × 1010 N/m2, c33 = 11.7 × 1010 N/m2, c44 = 2.33 × 1010 N/m2;
Piezoelectric stress constant:e31 = −6.5 C/m2, e32 = 23.3 C/m2, e33 = 17 C/m2;
Dielectric constant:ϵ11 = 1.503 × 10−8 C·V/m, ϵ33 = 1.503 × 10−8 C⋅V/mFigure 3
Finite element model of cantilever intelligent structural plate.


### 4.1. Configuration Analysis of the Genetic Algorithm: Optimizing the Number of Piezoelectric Actuators

According to (19), a minimum number of actuators should be used to satisfy structure shape control. Based on the finite element model and structure symmetry, 32 binary-encoded discrete design parameters G1 were set in the parameter cascade coding to represent configuration information regarding the location and number of piezoelectric actuators. The 32 float-encoded continuous design parameters G2 represented the working voltage of piezoelectric actuators. Working voltages were chosen within the range of -100,100, and the shape control precision (variance of actual displacement and expected displacement at the free end) was set as the restriction function. The shape control precision requirement was set as less than 1 × 10−3 mm2, the maximum number of generation as 1000, the crossover possibility of the genetic algorithm as Pc=0.625, and the mutation probability as Pm=1/64. Then, the optimization calculation of the genetic algorithm was performed.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 RX1,X2 provides variance of actual displacement and expected displacement at the free end. If the shape control precision requirement is more than 1 × 10−3 mm2, 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. Figure4 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 mm2. 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.Table 2
Optimal working voltage of piezoelectric actuators under optimal configuration.
Unit number
Working voltage (V)
6
59.4
10
72.1
11
−24.8
15
42.5
17
61.3
20
−36.7Figure 4
Variation of restriction function value.Figure 5
Variation of target function value.Figure 6
Variation of the fitness function value.Figure 7
Configuration diagrams of piezoelectric actuators at different generations: (a) 40 generations; (b) 125 generations; (c) 200 generations; (d) 1000 generations.
(a)
(b)
(c)
(d)Figure 8
Initial shape without control.Figure 9
Controlled shape when evolution is at 125 generations.Figure 10
Controlled shape under optimal configuration.


### 4.2. 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 mm2. 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 Pc=0.625, and the mutation probability as Pm=1/64. Then, the optimization calculation of the genetic algorithm was performed.Figure11 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 mm2. 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.Table 3
Optimal working voltage of piezoelectric actuator.
Unit number
Working voltage (V)
6
67.2
7
56.8
9
−19.2
11
35.2
14
41.2
15
46.5
16
21.6
17
−27.7
18
41.5
19
56.4
22
61.4
23
47.2Figure 11
Variation of restriction function value.Figure 12
Variation of fitness function value.Figure 13
Configuration diagrams of piezoelectric actuators at different evolution periods: (a) 20 generations; (b) 130 generations; (c) 320 generations; (d) 2000 generations.
(a)
(b)
(c)
(d)Figure 14
Controlled shape when evolution is at 320 generations.Figure 15
Controlled shape under optimal configuration.


## 4.1. Configuration Analysis of the Genetic Algorithm: Optimizing the Number of Piezoelectric Actuators

According to (19), a minimum number of actuators should be used to satisfy structure shape control. Based on the finite element model and structure symmetry, 32 binary-encoded discrete design parameters G1 were set in the parameter cascade coding to represent configuration information regarding the location and number of piezoelectric actuators. The 32 float-encoded continuous design parameters G2 represented the working voltage of piezoelectric actuators. Working voltages were chosen within the range of -100,100, and the shape control precision (variance of actual displacement and expected displacement at the free end) was set as the restriction function. The shape control precision requirement was set as less than 1 × 10−3 mm2, the maximum number of generation as 1000, the crossover possibility of the genetic algorithm as Pc=0.625, and the mutation probability as Pm=1/64. Then, the optimization calculation of the genetic algorithm was performed.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 RX1,X2 provides variance of actual displacement and expected displacement at the free end. If the shape control precision requirement is more than 1 × 10−3 mm2, 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. Figure4 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 mm2. 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.Table 2
Optimal working voltage of piezoelectric actuators under optimal configuration.
Unit number
Working voltage (V)
6
59.4
10
72.1
11
−24.8
15
42.5
17
61.3
20
−36.7Figure 4
Variation of restriction function value.Figure 5
Variation of target function value.Figure 6
Variation of the fitness function value.Figure 7
Configuration diagrams of piezoelectric actuators at different generations: (a) 40 generations; (b) 125 generations; (c) 200 generations; (d) 1000 generations.
(a)
(b)
(c)
(d)Figure 8
Initial shape without control.Figure 9
Controlled shape when evolution is at 125 generations.Figure 10
Controlled shape under optimal configuration.


## 4.2. 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 mm2. 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 Pc=0.625, and the mutation probability as Pm=1/64. Then, the optimization calculation of the genetic algorithm was performed.Figure11 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 mm2. 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.Table 3
Optimal working voltage of piezoelectric actuator.
Unit number
Working voltage (V)
6
67.2
7
56.8
9
−19.2
11
35.2
14
41.2
15
46.5
16
21.6
17
−27.7
18
41.5
19
56.4
22
61.4
23
47.2Figure 11
Variation of restriction function value.Figure 12
Variation of fitness function value.Figure 13
Configuration diagrams of piezoelectric actuators at different evolution periods: (a) 20 generations; (b) 130 generations; (c) 320 generations; (d) 2000 generations.
(a)
(b)
(c)
(d)Figure 14
Controlled shape when evolution is at 320 generations.Figure 15
Controlled shape under optimal configuration.


## 5. 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 float-encoded 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 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 high precision can be realized through optimization when using a certain number of piezoelectric actuators.



---

*Source: 6702183-2017-10-11.xml*