Final Thesis Report Rajeesh SK
Final Thesis Report Rajeesh SK
Chapter 1
Introduction
The common sources of mechanical stresses on aerospace structures are dynamic
loads and in fact dynamical load cycles can damage or cause a reduction in the service
life of aerospace structures. Therefore the investigation of vibration characteristics is an
important design phase of aerospace structures which are frequently experiencing
dynamic loading conditions. In essence, there are tremendous amounts of numerical and
experimental studies focused on investigation of the vibration characteristics and
attenuation of vibration levels of aerospace structures. When the frequency of the
dynamic loading matches with the natural frequency of the structure, the resonance
occurs, and it may cause severe structural vibrations. In this situation, severe vibrations
may damage components of aerospace vehicles, as aerospace structures are mostly light
weight and have low-stiffness characteristics.
The undesirable effects of induced-vibration in aerospace vehicles can be
exemplified by research studies for a fighter-jet, a helicopter and a satellite. Over the past
decade, research studies showed that severe vibrations in the form of buffet can damage
the components of a fighter-jet. Since flight envelope of fighter-jets includes many highly
acrobatic maneuvers and certain speeds higher than the speed of sound, severe vibrations
occur and may damage their components. The cracks in the components of fighter-jets
may cost millions to be replaced and maintained. On the other hand, helicopters are the
aerial vehicles whose structures are under dynamic loading in all flight envelopes because
of their rotary elements such as main rotor, tail rotor and transmission units. Their cabin
crew is exposed to the high levels of vibration in all flight zones and therefore there are
researches focusing on the investigation of vibration characteristics of a helicopter seat
Department of Electrical Engineering, College of Engineering Trivandrum
and effects on cabin crew health. Vibration at the helicopter seat causes excitation at the
natural frequency of the spine and abdominal of the cabin crew and exposition to these
types of vibrations for long time causes variety of health problems on cabin crew.
Satellites are other type of aerospace vehicles which are under dynamic loading during
launch and in-orbit operations. The induced-vibrations may cause both reduction of the
precision of pointing accuracy and cracks on the components of small satellites.
It is obvious that vibration suppression of structures is very crucial for better, safer
and easier life. Therefore, engineers are attempting to suppress such vibrations of
structures by using passive and active methodologies. However, passive vibration
suppression techniques are generally not suitable for low frequency applications.
Recently, active and adaptive vibration Control is receiving considerable attention as
alternative solutions to those passive methods. Passive vibration suppression
methodologies have some drawbacks such as not suitable for low frequency application,
increases weight of system and once designed their design parameters cannot be varied
easily The technological advances in piezoelectric materials also motivate scientists and
engineers to use these materials for the active vibration control as well. Piezoelectric
elements such as piezoceramics have excellent electric-mechanical conversion
characteristics. Therefore, they are widely used as sensors, which utilize the voltage
generated by the strain to which they are subjected, i.e., the piezoelectric effect. They are
also used as actuators, which utilize the strain due to the applied voltage, i.e. the inverse
piezoelectric effect.
The research and production of piezoelectric materials and piezoelectric devices
are rapidly developing. In recent years, active vibration control based on smart material
and structures especially piezoelectric smart structures has been dramatically developed
Department of Electrical Engineering, College of Engineering Trivandrum
for both research and engineering applications because of its good electro mechanical
coupling characteristics, preferable dynamic performance and higher sensitivity.
Now a new concept for the application of a piezoelectric element has emerged
from the field of material engineering. A Self-sensing actuator (SSA) is one in which a
single piezoelectric element functions simultaneously as a sensor and an actuator.
Therefore some researchers have tried to apply this idea to vibration control of flexible
structures. SSA is applied to bimorph cantilever beam which is an example of flexible
structure for vibration suppression.
vibration of a cantilever beam. Also studied the effects of different types of controller on
vibration. [3]
The modern technology demands the system to be light and reliable. However,
with the conventional vibration control techniques it is difficult to keep the system light
with all the damping mechanisms. This work deals with the Active Vibration Control
(AVC) of cantilever beam using piezoelectric (PZT) transducers. Active control involves
the use of sensors to sense the vibratory motion of the structure, a controller to generate a
control signal and an amplifier to amplify the control signal and an actuator which exerts
control force on the structure to reduce the vibrations. An experimental setup is made,
consisting of the aluminum cantilever beam with the PZT patches mounted on both the
sides of the beam. A proportional-integral derivative (PID) controller is designed to
generate the required control signal. A high frequency switch mode power converter is
designed to generate high voltage required for the actuator to produce the control force.
Experiments are performed for the active control of vibration [4]
utilization of LQR (Linear Quadratic Regulator) design. The LQR method is based on the
minimization of a quadratic performance index that is associated with energy of the state
variable and control signals. The goal of LQR controller design is to establish a
compromise between the energy state and control by minimizing a cost function [9].
Inverted pendulum has been the subject of numerous studies in automatic control
system. Since it the system is inherently non linear,it is useful to illustrate some of the
ideas in non linear control system. Wheeled mobile robots have in the recent years
become increasingly important in industry, since they provide a large degree of flexibility
and efficiency with respect to transportation and operation. The objective of this paper is
to design linear quadratic controllers for a system with an inverted pendulum on a mobile
robot.to this goal, it has to be determined which control strategy delivers better
performance with respect to pendulums angle and robots position, since it continually
moves towards an uncontrolled state.[10]
Linear Quadratic Gaussian (LQG) design problem is rooted in optimal stochastic
control theory and has many applications in the modern world which ranges from flight
and missile navigation control systems, nuclear power plants etc. It combines both the
concepts linear quadratic regulator for full state feedback and Kalman filters for state
estimation. Thus the configuration and design of a LQG controller for linear systems are
principally involved in both determining optimal process estimation by a linear quadratic
estimator and making an optimal control strategy by a linear quadratic regulator [11].
1.3 Objectives
The thesis report is organized in 9 chapters. After a brief introduction of the topic in chapter 1,
chapter 2 describes about piezoelectric materials. Piezoceramics and its properties are discussed
in chapter 3. Chapter 4 deals with the basics of cantilever beam. The concept of Self Sensing
Actuator is explained in chapter 5. Chapter 6 deals with the designing of controllers . Matlab
simulation results are included in chapter 7.Thesis conclusion is given in chapter 8. Finally, the
future scope of thesis is mentioned in chapter 9.
Chapter 2
generated and to apply the electric field to the piezoelectric element. When an electric
field is applied to the poled piezoelectric ceramic through electrodes on its surfaces, the
piezoelectric material gets strained (converse effect). The converse effect property is used
for actuator purposes. Figure 2.1 shows the converse piezoelectric effect.
Based on the converse and direct effects, a piezoelectric material can act as a
transducer to convert mechanical to electrical or electrical to mechanical energy. When
piezoelectric transducer converts the electrical energy to mechanical energy it is called as
piezo-motor/ actuator, and when it converts the mechanical energy to electrical energy it
is called as piezo-generator/ sensor. The sensing and the actuation capabilities of the
piezoelectric materials depend mostly on the coupling coefficient, the direction of the
polarization, and on the charge coefficients (d31 and d33). Figure 2.2 in the form of block
diagrams shows the transducer characteristics of the piezoelectric materials.
Department of Electrical Engineering, College of Engineering Trivandrum
Some of the typical piezoelectric materials include quartz, barium titanante, lead
titanate, cadmium sulphide, lead zirconate titanate (PZT), lead lanthanum zirconate
titanate, lead magnesium niobate, piezoelectric polymer polyvinylidene fluoride (PVDF),
polyvinyl fluoride (PVF). The piezoelectric ceramics are highly brittle and they have
better electromechanical properties when compared to the piezoelectric polymers. This
section gives in brief introduction about the various classes of piezoelectric materials:
single crystal materials, piezo-ceramics, piezo-polymers, piezo-composites, and piezofilms.
10
of a simple cubic cell that has a large cation A at the corner, a smaller cation B in the
body center, and oxygen O in the centers of the faces. The structure is a network of
corner-linked oxygen octahedral surroundings B cations.
For the case of Barium Titanate ceramic, the large cation A is Barium ion, smaller
cation B is Titanium ion. The unit cell of perovskite cubic structure of Barium Titanate is
shown in figure 2.3. The piezoelectric properties of the perovskite-structured materials
can be easily tailored for applications by incorporating various cations in the perovskite
structure. Barium Titanate
2.3.3 Polymers
The polymers like polypropylene, polystyrene, poly (methyl methacrylate), vinyl
acetate, and odd number nylons are known to possess piezoelectric properties. However,
strong piezoelectric effects have been observed only in polyvinylidene fluoride (PVDF or
PVF2) and PVDF copolymers. The molecular structure of PVDF consists of a repeated
monomer unit (-CF2-CH2-)n . The permanent dipole polarization of PVDF is obtained
Department of Electrical Engineering, College of Engineering Trivandrum
11
through a technological process that involves stretching and poling of extruded thin
sheets of polymer. These piezoelectric polymers are mostly used for directional
microphones and ultrasonic hydrophones application.
2.3.4 Composites
Piezo-composites comprised piezoelectric ceramics and polymers are promising
materials because of excellent tailored properties. These materials have many advantages
including high coupling factors, low acoustic impedance, mechanical flexibility, a broad
bandwidth in combination with low mechanical quality factor. They are especially useful
for underwater sonar and medical diagnostic ultrasonic transducers.
12
Chapter 3
Piezo Ceramics
A piezoelectric ceramic is a mass of perovskite crystals. Each crystal is composed
of a small, tetravalent metal ion placed inside a lattice of larger divalent metal ions and
oxygen. To prepare a piezoelectric ceramic, fine powders of the component metal oxides
are mixed in specific proportions. This mixture is then heated to form a uniform powder.
The powder is then mixed with an organic binder and is formed into specific shapes, e.g.
discs, rods, plates, etc. These elements are then heated for a specific time, and under a
predetermined temperature. As a result of this process the powder particles sinter and the
material forms a dense crystalline structure. The elements are then cooled and, if needed,
trimmed into specific shapes. Finally, electrodes are applied to the appropriate surfaces of
the structure.
Figure 3.1: Crystalline structure of a piezoelectric ceramic, before and after polarization
13
Figure 3.2. Poling process: (a) Prior to polarization polar domains are oriented randomly; (b) A very large DC
electric field is used for polarization; (c) After the DC field is removed, the remnant polarization remains
.
14
() =
(, )
(1)
15
where
Here
= ( + )/2
(2)
constant, t is the thickness, b is the width, (, ) the difference in the slopes of both
ends of the piezoelectric element. In this study (, ) is the slope of end of the
piezoelectric element.
Since the slope is proportional to the strain in the beam, is also proportional to
the strain.
The applied voltage causes the distributed bending moment M of the beam in a
proportional manner a s follow,
(, ) = () (0 )
= ( + )/2
(3)
(4)
16
17
Chapter 4
Beam Model
Fig 4.1 shows structure of bimorph cantilever beam fixed at x=0. The beam is
assumed to be a smart structure .In order to make beam bimorph type, piezoceramics are
bonded to the beam on both sides. This beam is composed of a metal shim (phosphor
bronze) between two piezoceramics sheets. The term y(x,t) denotes deflection of beam
and =
, =
(, ) =
() ()
(5)
18
+ + = ( () (0)
(6)
(7)
where M and denote the distributed bending moment and the natural frequency of ith mode
respectively
Where =
State space expression of beam the first 2 modes of vibration is given by [7].
= +
(8)
y=
(9)
State variable vector and output vector with states and respectively deflection and
slope are given by
[ ]
(10)
[ () ()]
(11)
Where
0
0
=[
0
0
0
()
= [ ()
0
1
0
0
1
]
0
0
0
= [ ()]
()
0
0
()]
()
19
This system is called BEAM with only one input, actuator voltage ( ) and four
outputs , , , which represent deflection and slopes of first and second modes of
vibration respectively. All the constants are given in Table 1.
20
Chapter 5
complete collocation of the sensor and the actuator, which is very advantageous
i n terms of control. Therefore some researchers have tried to apply this idea to
vibration control of flexible structures. Furthermore, it can reduce the size and
weight of the system because individual sensors and actuators are not needed.
An electric bridge circuit is essential for self-sensing actuation. Two
voltages, the sensor voltage due to the piezoelectric effect and the actuation
voltage due to the inverse piezoelectric effect, are mixed in a piezoceramics when
it is used as a SSA. It is therefore impossible to measure the sensor voltage
directly. The bridge circuit that includes the piezoceramics as one of the four
elements makes it possible to derive the sensor voltage that indicates the strain in
the structure, if the bridge is well balanced.
21
output the effects due to control voltage, only the voltage due to strain of the material
remains in the output signal. Therefore a RC bridge circuit (figure 4) which includes
piezoceramics is introduced to discriminate from . In RC bridge circuit, strain rate
(velocity) is sensed and strain rate feedback is effective in suppressing vibration.
( ) =
( () + ())
( )
(12)
(13)
Here the difference between and is defined as sensor voltage which is equal to -
And if the two RC time constants are equivalent ( =
22
( ) =
( )
(14)
1
, + 11. So can be expressed as
( ) = ( )
(15)
But
- =
( )
(16)
( )
(17)
The last term in the above equation is the same as the sensor signal - .Simplifying above
equation gives the actuator dynamics
=
For frequencies
( )
(18)
1
, will be equal to
23
[ ]=[
] * + + [
1
* + = *
10
][ ]
0
+* +
10
(19)
(20)
This system is called SSA with two inputs and two outputs. The inputs are time
derivative of slope of beam tip, and control voltage from the controller. The outputs
are sensor voltage to controller and actuating voltage applied to piezoceramics .The
two states in SSA are and + .
24
Figure 5 shows the plant including the systems SSA, BEAM and the blocks for
control specification. The dashed box represents the control object, and K is the
controller to be designed. Input disturbance to the beam results in the vibration of the
beam. The strain rate from the beam is sensed by SSA. Sensing voltage developed on
SSA is given to controller to obtain control voltage. The control voltage generated is fed
to SSA which will develop actuator voltage. This actuator voltage is used to suppress
vibration of the beam.
25
Chapter 6
Design of Controller
In this section, three control strategies are proposed and described in detail. The
control strategies are PID, LQR and LQG. The main objective of control strategy is to
suppress the vibration of beam due to disturbances.
( )= +
and
is given by,
(21)
= Proportional gain
= Integral gain
= Derivative gain
26
PID operations. The optimum values are obtained by tuning these parameters and it
depends on system and its response. Fig 6.1 shows the block diagram of the system with
a PID controller. The control signal from PID controller is again fed to SSA to generate
the actuating voltage to suppress the vibration. The controller parameters values are
0.9,
= 0.05 and
= 0.001.
27
In LQR, a control method is Linear Quadratic Control methods are best choice for
MIMO systems for effective vibration suppression and for damping out disturbances as
quickly as possible. The control voltage for SSA is determined by optimal control
solution of linear quadratic regulator (LQR) provided the full state vector is observable
LQR design is a part of what in the control area is called optimal control. For a LTI
system, this technique involves choosing a control law () = () which stabilizes
the origin (i.e., regulates x(t) to zero) while minimizing the quadratic cost function
= ( ()
where
() + () ())dt
(22)
system dynamics and the quadratic cost function. Q and R are weighting matrices for
states and control voltages respectively which are positive semi definite matrices. The
matrices Q and R are also called the state and control penalty matrices, respectively. If
the components of Q are chosen large relative to those of R, then deviations of x from
zero will be penalized heavily relative to deviations of u from zero. On the other hand, if
28
the components of R are large relative to those of Q , then control effort will be more
costly and the state will not converge to zero as quickly.
Hence, in an optimal control problem the control system seeks to maximize the
return from the system with minimum cost. It is crucial that Q must be chosen in
accordance to the emphasize that wanted to be given the response of certain states, or in
other words; how we will penalize the states. Likewise, the chosen value of R will
penalize the control effort u. As an example, if Q is increased while keeping R at the
same value, the settling time will be reduced as the states approach zero at a faster rate.
This means that more importance is being placed on keeping the states small at the
expense of increased control effort. On the other side, if R is very large relative to Q, the
control energy is penalized very heavily. Hence, in an optimal control problem the
control system seeks to maximize the return from the system with minimum cost. Since
the objective is to suppress the beam vibration, The weighting values corresponding to
these states are kept high that is the first two states in Q matrix corresponds to beam
vibration which is kept at high value.
The solution to above problem, the control voltage is
() = ()
where =
(23)
=0
(22)
The state matrix of plant (Beam and SSA) and weighting matrixes Q and R obtained by
tuning are given by
29
30
processes controllers and even nuclear power plants. This is simply a combination of
Kalman filter and a LQR controller. The separation principle guarantees that these can be
designed and computed independently. LQG controller can be used both in linear time
invariant as well as linear time variant systems. Figure 6.3 illustrates schematic of LQG
control approach
In most cases, all states of system will not be available for feedback. This problem
can be solved by Kalman filter (linear quadratic estimator).Thus in LQG control
technique, the control voltage for the actuators are determined by optimal control solution
of LQR problem of system with state estimated by a Kalman filter. It also takes in to
account the disturbances affecting the system.. The system with noise is given by
= + +
= + +
(23)
(24)
where A is the plant state matrix , B is the plant input matrix , C is the plant output matrix
, D is the plant feed forward, and both w and v are modeled as white noises associated
with process and measurement respectively.
() + () ()) dt
(25)
Where Q and R are weighting matrices which are positive semi definite matrices. Q
and R are weighting matrices for states and control voltages respectively. Hence, in an
optimal control problem the control system seeks to maximize the return from the system
31
with minimum cost. So in order to suppress the beam vibration, the weighting values
corresponding to these states are kept high.
The control problem is to determine the optimal control input;
u(t) = Kx(t)
(26)
=0
(27)
In real world control design problems, it is rarely possible to have access to all states
of the system which are needed for full state feedback. Instead, access is only possible to
specific measured outputs of the system. If these measurements carry enough information
about the states of the system, then a state observer using Kalman Filter could be
implemented to estimate all states of the system. This observer is capable of rejecting
disturbances of the system by acting as a low pass filter. The main inputs to the observer
are the control input (u(t)) and the system output (y(t)).
The state space equations of the Kalman Filter are shown in equation (28). It should
be noticed that it uses the same state space matrices (A, B and C) as the main system and
the estimated states are used as the system states.
= ++L ( )
(28)
The estimator gain must minimize the estimation error and is given by
where X is the positive semi definite solution to the algebraic Riccati equation:
Department of Electrical Engineering, College of Engineering Trivandrum
32
=0
(29)
From figure 6.3, it should be noticed that LQG is formed by connecting the system
and the Kalman Filter through the optimal state estimation gain (L) and then creating full
state feedback by using the estimated states ( ()) which passed through the optimal
feedback gain (-K). Thus Kalman filter can be applied to estimate state vector and
output vector y by using inputs u and measurements y. Because of the stochastic
separation principle, the previously mentioned gain could be designed individually.
The design process starts with checking controllability and observability of the pairs
(A, B) and (A, C), respectively. These criteria are necessary for the existence of the
solutions for the equations used to find the optimal gains. Then, the optimal state
estimation gain L (
and the solution of the Filter Algebraic Riccati Equation (FARE) shown in equation (29).
This solution ensures a minimum value of the cost function shown by equation (25).
After that, the optimal state feedback gain (K ) is calculated, where X is a positive semidefinite matrix and the solution of the Control Algebraic Riccati Equation (CARE) shown
in equation(27). This solution ensures a minimum value of the cost function in equation
(25).
The estimator and controller gain matrices can be designed separately by separation
principle. The state matrix of plant (Beam and SSA) and weighting matrixes Q and R
obtained by tuning are given by
33
The optimal state estimation gain L is obtained by considering 1% noise in the state variables
is
7.
10
06.45 0.07]
controller
is
designed. Also discussed in detail on designing of LQR and LQG optimal control
techniques.
Department of Electrical Engineering, College of Engineering Trivandrum
34
Chapter 7
Simulation Results
In this chapter, response of the system is analyzed in detail. The performance of the
system in open loop, closed loop, system with PID controller, system with LQR and
system analysis with LQG controller is done using commercial software package,
MATLAB/Simulink R2013b. Disturbance signal is given as input to the system in the
form of chirp signal. The frequency range of signal varies from 10Hz 100Hz. The beam
under analysis is having first mode resonant frequency of vibration at 91 Hz.
35
36
controller
37
LQG is an optimal controller that take in to account the noises affecting the
system. Figure 7.5 shows the tip vibration of the beam with LQG controller. The system
is showing best performance in terms of initial vibration(0.07mm) and steady state
vibration (table:2). With LQG controller beam tip vibration settles to
0.05
of final
Initial Tip
Displacement(mm)
Tip
displacement at
10 sec (mm)
0.20
0.10
0.15
0.04
0.10
0.03
0.07
0.025
0.02
0.006
38
Chapter 8
Conclusion
This paper presents a theoretical analysis of vibration suppression of bimorph
cantilever beam using SSA. Simple beam theory is used to model the beam. RC bridge
circuit which includes piezoceramics is used as SSA. SSA with conventional PID, LQR
and LQG controls are applied to suppress the vibration of beam. Simulation results
indicate that vibration of beam has been actively suppressed by applying control voltage
to SSA. From the results it is clear that performance of the system with LQG controller is
better than that of the system with LQR controller and conventional controllers.
39
Chapter 9
Future scope
The use of Piezoceramic for structural actuation and sensing is a developing area.
But piezoceramics imposes certain restrictions for its practical use in real world
applications. For instance brittle nature of Piezoceramic requires extra attention during
handling and bonding procedures. In addition the conformability to curved surface is
extremely poor requiring extra treatment of the surfaces. Active fiber composites like
macro fiber composites (MFC) provide not only the required durability and flexibility,
but also higher electro-mechanical coupling characteristics. So the effectiveness of MFC
on vibration suppression can be studied further.
40
References
[1] Dong Jingshi, Guo Kang, Zheng Wei, Study on Active Vibration Control Technique
of a
Cantilever
Vibration
[6]
control of Flexible
using Self-Sensing Actuator in JSME Intl. Journal, 1997, vol.40, No. 4, pp. 681-687,
41
42