 Open Access
 Total Downloads : 0
 Authors : Ranaivoarison Zo Nandrianina , Andrianaharison Yvon
 Paper ID : IJERTV7IS090077
 Volume & Issue : Volume 07, Issue 09 (September – 2018)
 Published (First Online): 05012019
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
Shape Control of Piezoelectric Bimorph Cantilever Beam Using Comsol Multiphysics and Genetic Algorithm
Shape Control of Piezoelectric Bimorph Cantilever Beam Using Comsol Multiphysics and Genetic Algorithm
Ranaivoarison Z. N.
Student
University of Antananarivo in STII Antananarivo, Madagascar
Andrianaharison Y.
Full Professor
University of Antananarivo in ESPA Antananarivo, Madagascar
AbstractThis paper presents a numerical study on optimal voltage, optimal load and optimal thickness for static shape control of piezoelectric bimorph cantilever beam. A finite element model using Comsol Multiphysics is developed to characterize the behavior of structure. Genetic Algorithm is used to calculate optimum values.
KeywordsFinite Element; Shape control; Piezoelectric; Comsol Multiphysics; Genetic Algorithm.

INTRODUCTION
Since a few years, investigations on smart structures having integrated selfmonitoring and control capabilities were carried out by many researchers. Among the various smart materials, particular attention is given to piezoelectric. Due to the direct and converse effects of piezoelectric materials, they have found applications in sensor and actuator technologies. One of the most important piezoelectric materials is Lead Zirconate Titanate (PZT). H.S. Kim, J.H. Kim, J. Kim [1] used PZT for vibration energy harvesting. Another application of PZT is shape control. The objective is to optimize some parameters in order to achieve the desired shapes. E.P. Hadjigeorgiou, G.E. Stavroulakis, C.V. Massalas [2] worked on shape control of beam using piezoelectric actuator. Genetic optimization is used to control the voltage applied to the actuator. In the paper of
H.R. Patel, J.R. Mevada [3], size, location and control voltage of actuators are taken as variables. Optimum values of all variables are calculated with Genetic Algorithm. Initial condition of beam is taken as horizontal position and three higher order polynomial curves are taken as desired shapes for beam to achieve. Piezoelectric bimorph cantilever beam was studied by L.N. Sulbhewar, P. Raveendranath [4]. A Timoshenko piezoelectric beam finite element based on First order Shear Deformation Theory (FSDT) was developed. It is an improved model of the FSDT proposed by S. Narayanan, V. Balamurugan [5].
In the present work, piezoelectric bimorph cantilever beam is investigated using Comsol Multiphysics and the results are compared to analytical results based on FSDT. Genetic Algorithm carried out using Matlab works with Comsol Multiphysics to optimize voltage, load and thickness for shape control. Error between desired shape and achieved shape is taken as an objective function.

METHODOLOGY
We firstly introduce the equations which were used to develop the piezoelectric model. This model is validated by
comparing the results got using Comsol Multiphysics with analytical results. In a second part, we take advantage of the power of Comsol Multiphysics combined with Matlab [6]. The developed piezoelectric model is used with Genetic Algorithm for shape control.

Piezoelectric equations and model in Comsol
Using the first law of thermodynamic, the internal energy of an electromechanical system can be expressed as:
dU=Tij dSij+Ei dDi
Where:
U: internal energy;
Tij: stress;
Sij: strain;
Ei: electric field;
Di: electric displacement. Electric enthalpy is defining as:
H=U – Ei Di
H: Electric enthalpy
The total differential of H gives us:
dH= Tij dSij – Di dEi
Equation (3) can be written as:
In linear piezoelectricity theory, the form taken by H is:
Where:
cEijkl: constant elastic;
ekij: constant piezoelectric;
S
ij: constant dielectric.
The piezoelectric constitutive equations are:
Tij =cEijkl Skl – ekijEk
ij
Di =eikl Skl + S Ej
The piezoelectric model is implemented using Comsol Multiphysics. We consider here a bimorph cantilever beam with oppositely poled piezoelectric layers. The structure has clamped end at the lefthand side. We used three types of material while the geometry is fixed: PZT 4, PZT 5A and PZT 5H. Those materials are available in the library of Comsol Multiphysics.
In order to compare Comsolbased simulation results with analytical results based on FSDT, we investigated static analysisactuator configuration and static analysissensor configuration.

Static analysisactuator configuration:
Interface of the bimorph is grounded and the voltages of Â±10 V are applied on the free surfaces.

Static analysissensor configuration:
The piezoelectric bimorph beam is subjected to a tip load of 1000 N.
Displacements of the piezoelectric bimorph cantilever beam are given in the results of Comsol Multiphysics. Tip deflection is compared to analytical results for the three materials.
A schematic view of the piezoelectric bimorph cantilever beam is presented in Fig. 1.
Fig. 1. Piezoelectric bimorph cantilever beam
Dimensions of the top and bottom piezoelectric are the same and given in TABLE I.
TABLE I. DIMENSIONS OF TOP AND BOTTOM PIEZOELECTRIC
Top piezoelectric
Bottom Piezoelectric
Length (Ã—103m)
LT = 100
LB = 100
Thickness(Ã—103m)
hT = 5
hB= 5


Shape control
The model developed in Comsol Multiphysics is used with Genetic Algorithm for shape control. Initial condition is taken as horizontal. Desired shape is taken as polynomial curve. Error between desired shape and achieved shape is to be considered as an objective function to minimize. Error can be expressed as:
(9)
Where:
i: predefined displacement;
qi: achieved displacement;
Genetic Algorithm is a class of probabilistic algorithm based on biological evolution. It starts working on a randomly generated set of solutions, known as population. The solution is called chromosome. Genetic Algorithm generates new population of chromosomes by selecting the better fit solutions from existing population and applying genetic operators (crossover, mutation) to produce new offspring of the solutions.
Here, Genetic Algorithm is carried out using Matlab.
Methodology associating Comsol Multiphysics and Matlab is summarized in Fig. 2.
Fig. 2. Methodology using Comsol Multiphysics and Matlab
For Genetic algorithm, population is set as 100 for all cases.
Number of generation is allowed as 50.
In the shape control, PZT 4 is used as piezoelectric material. We consider three kinds of shape control problems with various boundary conditions: the voltage problem, the tip load problem, the voltage and thickness problem.

Voltage problem:
The first problem consists in finding voltage which minimizes the objective function under the constraint:
Actuation voltage
min max
min: Lower saturation voltage
max: Upper saturation voltage
In order to test the efficiency of methodology associating both Comsol Multiphysics and Matlab, the desired shape that we considered in the voltage problem is displacement corresponding to Â±50V as shown in Fig. 3.
Fig. 3. Initial condition and desired shape
The desired shape used above is the same considered in the tip load problem. It consists to find tip load which minimizes the objective function.
This work shows how tip load are depended on the desired shape.

Voltage and thickness:

In this problem, two variables are considered: voltage and thickness. The thickness of top piezoelectric hT is varied while the thickness of bottom piezoelectric hB is fixed. We investigated two cases:
Thickness is under constraints: 1Ã—103m) hT 4Ã—103m Thickness is under constraint 6Ã—103m) hT 15Ã—103m.
The first case corresponds to thickness of top piezoelectric inferior to the thickness of bottom piezoelectric. The second case corresponds to thickness of top piezoelectric superior to the thickness of bottom piezoelectric.
This work shows how thickness and control voltage are depended on the desired shape.


RESULTS AND DISCUSSIONS

Comparison of tip deflection got in model Comsol and analytical results
The displacement caused by an applied voltage Â± 10V is shown in Fig. 4 and the displacement caused by a tip load of 1000 N is shown in Fig. 5 for PZT 4, PZT 5A, PZT 5H.
Fig. 4. Displacement by Â± 10V
Fig. 5. Displacement by a tip load of 1000N
Absolute tip deflection of different piezoelectric materials got in Comsol is in agreement with analytical results based on FSDT [4] as shown in the TABLE II and TABLE III.
TABLE II. COMPARISON OF TIP DEFLECTION CAUSED BY AN APPLIED VOLTAGE Â±10V
Materials
FSDT (Ã—106m)
Comsol (Ã—106m)
Error percentage (%)
PZT 4
0.3620
0.3559
1.68
PZT 5A
0.5023
0.4927
1.91
PZT 5H
0.7899
0,7801
1.24
TABLE III. COMPARISON OF TIP DEFLECTION CAUSED BY A TIP LOAD OF 1000N
Materials
FSDT (Ã—106m)
Comsol (Ã—106m)
Error percentage (%)
PZT 4
46.061
44.484
3.42
PZT 5A
60.942
58.486
4.03
PZT 5H
57.156
56.625
0.92

Genetic Algorithm and Comsol used for shape control

Voltage problem
This first problem is used to test the efficiency of the methodology associating Comsol Multiphysics and Matlab. The desired shape corresponds to displacement caused by an applied voltage Â±50V.
Constraints: 0
Convergence with error equal to 1.34*104 is achieved at 15th iteration.
Error is plotted with number of generations in Fig.6.
Fig. 6. Error for control voltage

Tip Load problem
Here, the same desired shape used previously is considered. Tip load is taken as variable. The tip load which can produce the desired shape is 37N.

Voltage and thickness problem

Thickness of bottom piezoelectric is fixed. The Thickness of top piezoelectric hT and voltage are variables.
– Case 1: thickness of top piezoelectric is taken as inferior of the thickness of bottom piezoelectric.
Constraints: 1Ã—103m hT 4Ã—103m and 0
Optimum values: hT =4Ã—103m and 40V.
If we compare applied voltage got in this case with results in the control voltage only, we see that decrease in thickness of top piezoelectric needs a smaller applied voltage to achieve the desired shape.

Case 2: thickness of top piezoelectric is taken as superior of the thickness of bottom piezoelectric.
Constraints: 6Ã—103m hT 15Ã—103m and 0
Optimum values: hT =6Ã—103m and 60V
If we compare applied voltage got in this case with results in the control voltage only, we see that increase in thickness of top piezoelectric needs a greater applied voltage to achieve the desired shape.


CONCLUSION
Here, methodology using Comsol Multiphysics and Genetic Algorithm is developed for shape control. To minimize control cost of piezoelectric bimorph cantilever beam, it is required to optimize voltage, tip load and thickness on the structure. Optimization is carried out for the error minimization between desired and actual shape with the use of genetic algorithm.

Tip deflection got with model developed in Comsol Multphysics is in agreement with analytical results based on FSDT;

Efficiency of methodology associating Comsol Multiphysics and Genetic Algorithm was tested by control voltage and the results are satisfactory.

An equivalent tip load can cause the same desired shape got by applied voltage

Decrease in thickness of top piezoelectric needs a smaller applied voltage to achieve the desired shape. Increase in thickness of top piezoelectric needs a greater applied voltage to achieve the desired shape.
REFERENCES
[1] H.S. Kim, J.H. Kim, and J. Kim, A review of piezoelectric energy harvesting based on vibration, International Journal of Precision Engineering and Manufacturing, vol. 12, pp. 11291141, December 2011. [2] E.P. Hadjigeorgiou, G.E. Stavroulakis, and C.V. Massalas, Shape control and damage identification of beams using piezoelectric actuation and genetic optimization, International Journal of Engineering Science, vol. 44, pp. 409421, April 2006. [3] H.R. Patel and J.R. Mevada, Shape control and optimization using cantilever beam, International Journal of Engineering Research and Applications, vol. 3, pp. 155161, June 2013. [4] L.N. Sulbhewar and P. Raveendranath, A Timoshenko piezoelectric beam finite element with consistent performance irrespective of geometric and material configurations, Latin American Journal of Solids and Structures, vol. 13, pp. 9921015, 2016. [5] S. Narayanan and V. Balamurugan, Finite element modelling of piezolaminated smart structures for active vibration control with distributed sensors and actuators, Journal Sound and Vibration, vol. 262, pp.529562, June 2003. [6] S.Ozana, M. Pies and R. Hajovsky, Using matlab and comsol multiphysics for optimization of the model of underground thermal processes at old mining dumps, Applied Mechanics and Materials, vol. 548549, pp. 571578, April 2014.