Vibration Control Using Heat Actuators

This paper focuses on theoretical investigation of active vibration control of a cantilever beam using heat actuation. The actuator is a thin metal bar rigidly bonded to the beam on one face and subject to heat input on the opposite face. The actuator then works like a piezoelectric actuator, and expands and contracts in response to applied heat. We assume that the actuator is insulated so that no heat is transferred to the beam, ensuring that the heat does not alter the beam’s thermal state. To avoid necessity of cooling, we consider two actuators working together at the same spanwise location, one on the upper and one on the lower face of the beam. Then, the beam can be bent up and down by applying heat to the lower and upper actuators, respectively. The governing equations are partial differential equations for one-dimensional heat conduction of the actuators and the bending vibration of the beam with attached actuators. For an approximate solution, Rayleigh-Ritz method replaces the partial differential equations with a system of ordinary differential equations. A control model is obtained from a low-dimensional representation of the system, and used to design feedback control and observer by means of LQR and Kalman-Bucy filtering techniques. The control signal obtained is introduced into the plant model, a high-dimensional representation of the system, to mimic the true system as closely as possible. In a numerical application, the response of the beam to an initial excitation is simulated, which demonstrates that the heat actuators are in fact effective in active vibration control of the beam.


Introduction
Suppression of vibration in flexible structures remains to be a challenging problem in the related technical community. The current trend is to use active means to suppress the vibration. The active control of vibration in structures is implemented by means of sensors, actuators, an accurate mathematical model of the whole system (flexible structure + sensors + actuators), and a computer to carry out all required computations in real time.
Sensors do all measurements, but usually the states of the system cannot be measured directly. Hence, the mathematical model of the system is used to estimate all of the states based on the available measurements, provided that the system is observable. The mathematical model is also used to compute the control signal to be inputed to the actuator based on these estimates of the states. The actuator produces what is needed to affect the true system to suppress the vibration. This can be force, moment, strain, etc. The actuators widely used in vibration control are hydraulic, pneumatic, and piezo-electric, with each having its own benefits and disadvantages depending on the specific applications in which they are used [1]. A special attention has been given to piezoelectric actuators in the recent years. These actuators are bonded on surfaces of structures at some discrete locations to alter the local strain therein by the application of proper voltage signals to actuators [2].
Another type of actuation that has started receiving some attention in vibration control is simple heat. Considering a structure initially at room temperature, an applied load causes temperature to vary over the body of the structure nonuniformly. Areas of the structure under tension are now cooler and areas under compression are hotter than the room temperature, creating a temperature gradient in the structure. This in turn causes a heat flow by conduction from hotter areas to cooler areas. Hence, a portion of the mechanical energy introduced to the structure is converted to heat through this process. Since the heat flow created is irreversible, the portion of the mechanical energy converted to heat is dissipated. This phenomenon is known as thermoelastic damping [3]. With this fact, question remains as to whether vibration in a structure can be suppressed by artificially introducing a thermal gradient in the body by means of some external heat actuation. Some earlier studies had shown that, when applied properly, external heat actuation could in fact introduce significant amount of thermoelastic damping in structures [4]- [9].
Pioneering research by Edberg [4] considers control of vibration in a beam using a thermoelectric actuator, which has an ability to heat or cool the beam to produce a thermal gradient across the thickness of the beam. Experimental results reported in the paper show that with a simple lead compensation, damping ratio for the first vibrational mode of a cantilever beam is increased by more than nine times. Another work by Murozono and Sumi [5] also investigates active vibration control of a cantilever beam, but this time by heat input into the beam through the upper and lower surfaces, solely for creating a temperature gradient across the beam section. Foil strain gauges bonded on the surfaces are used as heat actuators. The heat input is through the upper or lower surface, depending on the direction of the gradient needed. Experiments reveal that the damping ratio for the first vibration mode is increased by about ten times. A computational study by Tuzcu et al. [6] considers a uniform rod in longitudinal vibration and shows that control of both vibration and temperature can be achieved by heat input through thermoelectric actuators bonded on the lower and upper surfaces of the rod. Tuzcu and Gonzalez-Rocha [7] extended the work of Tuzcu et al. [6] to a beam in bending. Both vibration and temperature are controlled simultaneously using thermoelectric actuators. Zhang et al. [8] proposes a vibration control of thermally induced vibration in space structures. Actuators are heaters bounded on the spacecraft structure at some optimal positions. Tuzcu et al. [9] explores the idea of using heat patches as actuators working in pairs on the upper and lower surfaces at some discrete spanwise locations of a uniform beam to simultaneously control vibration and temperature in the beam. A numerical example shows the presence of thermoelastic damping and demonstrates simultaneous control of elastic displacement and temperature gradient using a set of displacement measurements.
References [4]- [9] have a major disadvantage: all use heat actuation that requires heating the flexible structure to be controlled, which has a potential to alter the thermal state of the structure significantly, which in turn can prevent the structure from functioning properly. To overcome this problem, this paper proposes a simple design of an actuator which includes a thin metal bar bonded to the structure to be controlled on one face and subject to heat input on the opposite face, as seen in Figure 1. Then, the actuators will work somewhat like piezoelectric actuators, expanding and contracting in response to applied heat, instead of electric voltage in the case of piezoelectric actuators. For more effective control of the structure, multiple of these actuators can be rigidly bonded to the upper and lower surfaces of the beam in pairs at a desired number of spanwise locations. It is assumed that the actuators are insulated so that no heat is transfered to the flexible structure. Since each actuator and the structure are rigidly attached to each other, extension in an actuator due to the application of heat is transferred directly to the host structure. The beam will then bend down if heat is applied to one or more upper actuators, and it will bend up if heat is applied to one or more lower actuators. Hence, the beam's deflection can be controlled by choosing proper heat inputs to the actuators.

Modeling
We will demonstrate our modeling and control approach on a cantilever beam, which is a fundamental element commonly encountered in structures. To make the presentation of our approach more clear and focus, we limit ourselves to a simple model by assuming the beam is uniform and there is only one pair of actuators, as shown in Figure 1. The beam is slender so that Euler-Bernoulli beam theory is valid. The upper and lower actuators are assumed to be rigidly bonded to the beam, which means the displacements of the beam and the actuators are equal at every point on the bonding surfaces. We define the displacement field in the beam as where , , u v w are the displacement components in , , x y z -directions, respectively, and subscript x denotes a partial derivative with respect to x. Displacement components of a typical point P on the beam are illustrated in Figure 1. The beam's left-end is fixed and its right-end is free so that w is subject to the boundary conditions: where each subscript x indicates a partial derivative with respect to x. Then, the strain components of the beam can be expressed as 2 2 , 0 Since the ends of the actuators are free, the normal strains of the actuators will be constant at beam's strains at its interface with the actuators [10], which means Assuming a linear stress-strain relation, we can express the only nonzero stress component in beam as where E is the Young's modulus of the beam. On the other hand, the normal stresses in the actuators are where a E is the Young's modulus and α the coefficient of thermal expansion for the actuator material, and L T and U T are the temperature deviations from the ambient temperature 0 T in the lower and upper actuators, T T z t = . We approximate w, U T , as well as L T using the expansions are the matrices consisting of m and n number of shape functions and 3 7.8548 L β = and 4 10.9955 L β = .
Let us now express the total potential energy, which is the summation of potential energies of the beam and the lower and upper actuators: where a A bh = is the cross-sectional area of the each actuator, V is the volume of the beam, L V and U V are the volumes of the actuators, and is the area moment of inertia of the beam. Substituting Equation (7) into Equation (11), we rewrite the potential energy expression as We note that K is the stiffness matrix. Kinetic energy of the beam is where the mass matrix is where ρ is the density of the beam, m mass per unit length of the beam, and a ρ density of the actuator material. Generalized force is obtained from the virtual work: is the generalized force vector. The equations of motion of the beam can be derived by means of Lagrange's equations of motion, which can be written in the vector as d dt (12) and (14) into the Lagrangian, and the result into Equation (18), we get Equation (19) governs the elastic deflection of the beam. Since this equation is coupled with the temperature difference θ , it must be solved simultaneously with the heat equations governing the temperatures in the upper and lower actuators. Hence, we next derive these heat equations. The heat conduction in the upper actuator is described by the boundary-value problem defined by the following partial differential equation and boundary conditions: is the surface area of the upper face. We choose the shape functions ( ) ( ) Substituting the second of Equation (7), premultiplying both sides by T κΨ , and integrating the result over the domain of the upper actuator, we get Similarly, substituting the second of Equations (7) into the boundary conditions, we get Using the integration by parts, the second integral in Equation (22) can be written as Substituting Equation (24), Equation (22) can be rewritten in the compact form Equation (25) completely describes the temperature in the upper actuator for prescribed heat input ( ) U Q t . Note that by neglecting the thermoelasticity of the actuator, we assume the temperature is not affected by the elastic deformation of the beam. Equation (25) can be written in the state-space form 1 1 which then leads to the following scalar equations for the individual generalized coordinates: Hence, the eigenvalues of the heat conduction are ( ) We note that the eigenvalue corresponding to 1 i = is zero and all the remaining ones are negative real so that U T approaches to a constant, instead of zero, as 0 U Q → . The heat conduction in the lower actuator is similarly described by the boundary-value problem Following the same steps, Equation (22)

Plant Model
Equations (19), (25) which can be rewritten in the state-space form are the state and input vectors, and 1 1 1 1 are the coefficient matrices. Note that the dimensions of x , u , A, and B are ( ) 2 m n m n + × + , and ( ) ( ) 2 2 m n m + × + , respectively. We will call the state-space model of Equation (32) the plant model, which, in absence of an experimental model, will be used to simulate the behavior of the actual system and to test the performance of our control design. In order for this model to approximate the behavior of the actual system better, we will choose m and n (the number of shape functions used in the expansions of Equations (7)) as high as possible. Control design, on the other hand, will be based on a control model that uses a less number of shape functions, as explained next.
Our objective in this paper is to control the deflection of the beam using

Control Model
Control design must be as simple as possible for easy implementation and fast real-time computations. This in turn requires a low-order control model. To that end, we first note that the second integral of the expression for a C defined in Equation (13) (38) are the control state vector and the control input, and The dimensions of c x , c u , c A and c B in this case are ( ) To keep the dimension even smaller, we will use an m for the control model several times smaller than for the plant model. Note that this model allows us to control not only the deflection in the beam, but also 1 θ , which is the most dominant temperature component among all 1 , , n θ θ .

Output Feedback Control Design
The system to be controlled is inherently underactuated because cantilever beam is a distributed-parameter system with infinite degrees of freedom while the number of control inputs that can be used must always be finite. This inherent property of the system remains a fact even though we use a finite-order approximation of the actual system, namely the control model, because the order of the control model is 1 2m + while the x is the estimate of the actual output ( ) t y [12]. Using the Kalman-Bucy Filter, ∆ is obtained from where W is the sensor noise intensity matrix and Σ is the solution of the algebraic matrix Riccati equation in which V is the state excitation noise intensity matrix. We assume deflection measurements at the points 1 P and 2 P shown in Figure 1 are available as output signals. Since we use the plant model to approximate the actual system, the output of the plant can be written as is the output matrix relating the actual state to the output. The estimate of the output, on the other hand, is written as is the output matrix relating the estimated state to the estimated output. Once

Numerical Application
The modeling approach and the control design presented above will be numerically demonstrated for a uniform cantilever beam. We will use the same cantilever beam and the related parameters used in [9] so that we can compare the results from this paper directly to the ones from [9]. Physical parameters of the beam and the actuators are listed in Table 1. For most effective actuators, we choose 0 l x = and r x l = so that the actuators are closed to the root where the strain is high. As revealed in Table 1, the material of the beam is chosen to be steel. We choose steel also for the actuators because it has high a E α , which makes the actuators more effective.
Using the values of the parameters listed in Table 1, we construct the matrices M, K and a C using the expressions given in Equations (13) The matrices  ,  and  are as given in Equation (26). A and B matrices given in (34) can be constructed by using the numerical forms of the submatrices: To demonstrate the simulation of the closed-loop system, we choose 7 m = and 10 n = for the plant model and 3 m = and 1 n = for the control model so that A is 34 34 The corresponding eigenvalues of open-loop plant and control models, closed-loop control model, and observer dynamics are listed in Table 2. The open-loop plant model has 34 eigenvalues while the control model has only 7. The zero and the negative real eigenvalues of the plant correspond to the heat conduction of the actuators, as explained earlier, and purely complex ones correspond to the vibration of the beam. The control model has a zero eigenvalue and three pairs of purely complex conjugate eigenvalues. The third and the fourth columns of Table 2 show how these eigenvalues are stabilized by the control and observer designs, respectively. However, the question remains as to whether the control and observer designs stabilize the plant. This question   Figure 2 completely describes how the control of the plant is implemented. We will test the performance of our control design by simulating the response of the closed-loop plant to the same initial conditions used in [9], namely , we conclude that the damping achieved in the present paper is about the same as that of [9]; in both, tip displacement reduces to 5 mm in approximately 16 cycles. Figure 4 shows the heat inputs to the lower and upper actuators needed to achieve this damping. Considering that the amplitude of the vibration in this numerical example is high, these heat inputs seem to be reasonable. Even though [9] uses heat actuators in two spanwise locations, the required heat inputs are quite similar to those of this paper.
In this paper, we assume each actuator is insulated in all faces except the face that is subject to the heat, including the faces normal to x and y axes. This assumption clearly means the temperatures everywhere in the actuators increase monotonically until the heat inputs are ceased after vibration amplitude reduces to a desired value. As an example, temperatures Since the heat is confined to a small volume in the actuators, the temperatures get very high, but still much less than the material of the actuators can tolerate. To prevent the temperatures in the actuators to get very high values, some sort of cooling, such as thermoelectric cooling, can be used to cool the actuators. In that case, cooling must be introduced into the control design as a disturbance. Contrary to the actuators, no temperature changes occur in the cantilever beam. In [9], since the heat is directly applied to the beam itself, which has a much larger volume than the actuators of this paper, the temperature increase in the beam is limited to few degrees.
We can also compare the total heat injected into the volumes in the two papers. we compute the total heat input to the actuators within t seconds as the integral of total heat rate: q t versus t is plotted for 100 seconds in Figure 6, from which we read the total heat input to the actuators in 100 s as ( ) On the other hand, in [9], ( ) 100 2.555 kJ q = , which is much larger than the heat input in the present paper.

Conclusions
This paper investigates active vibration control of a structure using a new type of actuator. The structure of interest is a cantilever beam. The actuator, on the other hand, is a thin metal bar rigidly bonded to the beam on one face and subject to heat input on the opposite face. The actuator then can work like a piezoelectric actuator, and expand and contract in response to applied heat. We assume that the actuator is insulated so that no heat is transferred to the beam. This ensures that the heat applied to the actuator does not alter the thermal state of the beam. One problem that needs to be addressed is that expansion of the actuator requires heating while contraction requires cooling. However, cooling is not practical, and to avoid it, we consider two actuators working together at the same spanwise location, one on the upper and one on the lower face of the beam. Then, the beam can be bent down by applying heat to the upper actuator, and bent up by applying heat to the lower actuator.   In the first part of the paper, we obtain a comprehensive model of the system that can be used to design control and simulate the response of the closed-loop system. Euler-Bernoulli beam theory is assumed to be valid to describe the motion of the beam with the attached actuators. Then, the Lagrange's equations of motion are used to derive the governing equations for the vibration of the beam with the actuators. Since we assume the elastic motions of the actuators only depend on the elastic motion of the beam, the actuators do not bring any additional elastic degrees of freedom into the governing equations. One-dimensional heat conduction equation is used as the equation governing the temperature in each actuator. We also assume the temperatures of the actuators are not affected by the elastic motion. The governing equations are a set of three partial differential equations, each of which has an infinite dimension. For an approximate solution, we use Rayleigh-Ritz method to replace the partial differential equations with a desired number of ordinary differential equations, obtaining a finite, but otherwise high-dimensional system. The higher the dimension, the more accurate the model is.
We obtain the plant model by using a high-dimensional representation of the system so that it can be represented as accurately as possible. On the other hand, we obtain the control model by using a low-dimensional representation of the system so that estimates of the states as well as control inputs can easily be computed in real time. Since we do not have an experimental model, we use the plant model to mimic the true system as closely as possible, which is the reason to choose a high dimensional model. The control model is used to design feedback control and observer by means of LQR and Kalman-Bucy filtering techniques. The observer design is based on only two displacement measurements, one in the middle and one at the free end of the beam.
In a numerical application, control signal is computed using the control and observer gains based on the control model, and inputted into the plant. Then, the response of the beam to an initial excitation is simulated. The numerical simulation demonstrates that the actuators are in fact effective in active vibration control of the beam. The results obtained are comparable to some of the results of [9]. For comparison reason, the feedback control is designed so that approximately the same damping factor as that of [9] is achieved. However, we show that the total heat input needed to achieve this damping is much smaller than that of [9], from which we conclude that the actuators of this paper are more effective than those of [9].
For simplicity of the models, we assume that each actuator is insulated in all faces except the face it is subject to heat. With this assumption, the heat conduction is one-dimensional (only in z-direction), and also we do not have to make assumptions about the convection and/or radiation heat transfer from the faces normal to x and y axes to outside. However, in real-life applications, this assumption is hard to enforce and it must be relaxed. Our future work will aim to make design alterations to relax this assumption to make the work more applicable to real-life situations.