#### Abstract

This paper proposes different design strategies of robust controllers for high-order plants. The design is tailored on the structure of the equations resulting from modeling flexible structures by using modal coordinates. Moreover, the control laws have some characteristics which make them specially suited for active vibration reduction, such as strong stabilization property and bandpass frequency shape. The approach is also targeted the case of more sensors than actuators, which is very frequent in practical applications. Indeed, actuators are often rather heavy and bulky, while small and light sensors may be placed more freely. In such cases, sensors can be usefully placed in the locations where the primary force fields act on the structure, so as to provide the controller with a direct information on the disturbance effects in terms of structural vibrations. Eventually, this approach may lead to uncolocated control strategies. The design problem is here solved by resorting to a Linear Matrix Inequality technique, which allows also to select the performance weights based on different design requirements, for example, a suitable bandpass frequency shape. Experimental results are presented for a vibration reduction problem of a stiffened aeronautical panel controlled by piezoelectric actuators.

#### 1. Introduction

Reduction of vibrations in flexible systems has been considered for years a challenging problem, due to some specific aspects. First, as the model of the structure is described by partial differential equations that can be treated as an infinite set of ordinary differential equations, any finite-order controller must rely on a limited (truncated) description of the plant. Thus, neglected dynamics actions must be taken into account by robustness properties of the controller, otherwise, spillover effects [1] may lead to control problems, even to the instability of the closed-loop system, in the most unfavorable situation. A second problem is that, in order to base the design technique for the controller on a meaningful model of the structure, many structural modes must be included in the model, especially for structures with high modal density, so that control problems for flexible structures are intrinsically large-scale problems. Only in the simpler case of low modal density nonoptimal control strategies (e.g., sliding manifold-based control [2] may be used). Finally, stability of the controller (in order to reduce sensitivity to disturbances) and limited bandwidth, for example, to avoid spillover, due to interaction of unmodelled dynamics with control action, must be enforced. From a practical point of view, vibrations are often measured by using accelerometers, which are prone to offsets which can even saturate actuators if not suitably filtered out in the feedback loop. Thus, the need for a filtering action both at high and at low frequency arises. Hence, a limited-bandwidth strongly stabilizing controller is needed. Further discussions on this topic can be found in [3, 4]. Furthermore, in [5] a broader list of references on the topic is presented and discussed.

In this paper the selection of a stable stabilizing controller with bandpass frequency shape is discussed. First, optimal controllers are discussed, with analytical expression in closed form that is useful when high-order models may result into increased computational burden. The main limitation of this approach is that it can deal only with square plants, that is, as many inputs as outputs must be present. Next, a further technique, able to deal with rectangular plants (i.e., with more outputs than inputs), is proposed. The motivation is provided by the requirements of many applications, where the rejection of a disturbance acting on the structure is desired in locations where the installation of actuator/sensor couples is prohibited. Therefore, the possibility to exploit the vibration measurement in those locations appears the only viable option. The second control strategy proposed in this paper is specifically designed to usefully exploit the knowledge of the disturbance entry point on the structure to place additional sensors able to provide the controller with the necessary information on the vibration status caused by the disturbance. The resulting controller thus has more inputs than outputs, but it is still guaranteed to be strongly stabilizing and with bandpass frequency response. Moreover, selection of free control parameters in order to improve effectiveness of the control, by emphasizing the action of the most controllable and observable modes, is discussed. To achieve such an objective, the control design procedure, still based on an approach, adopts an LMI (linear matrix inequality) technique to compute the control parameters to allow the user to select the design weights based on various design requirements or by taking into account modal controllability/observability indices. The paper reports experimental results of broadband vibration reduction for a stiffened aeronautical panel controlled with piezoelectric actuators.

A final consideration is in order in this section. It is clear that the placement of the actuators and sensors is of paramount importance in the whole response of the flexible structure to external disturbances. Thus, an efficient control system should in general address also the problem of the optimal allocation of actuators and sensors, for example, based on maximizing vibrational energy reduction. However, in practical situations the position of sensors and actuators is often constrained by different considerations (geometry, weight limitations, safety, and cost); thus the control designer has very limited degrees of freedom in deciding sensors and actuators position and even number. For this reason, in this paper the location of sensors and actuators is assumed to be given.

#### 2. The Model

In this paper a flexible structure with control inputs and measured (control) outputs is considered. The structure is assumed to be infinite-dimensional, linear, and time invariant. Since in general , the system is not necessarily square, that is, it may have more outputs than inputs. This assumption is quite reasonable in flexible systems control applications, where actuators are heavy and expensive, while sensors (accelerometers) are light, small, and relatively cheap.

As extensively discussed in [5], the mathematical model of a flexible structure can be written as where is the state vector, the control input, and the measured (control) output. Obviously, since we are approximating an infinite-order system, the number of retained modes is generally very high (with systems of order of hundreds). By resorting to the “modal coordinates’’ the matrices of the state space model assume the form where with and being the natural frequency and damping coefficient of the th mode (). In this case are the modal displacements and the modal velocities. Moreover, assuming that the velocity measurement is available, the output matrix has the form (4). By following the gray-box approach proposed by Cavallo et al. [5], the model above can be accurately identified through a measurement campaign carried out on the structure. The gray-box approach assumes that the structure of the system to identify is exactly the one in (2)–(4), and this assumption is of paramount importance in the design of the controller, as it will be illustrated in Section 3.

#### 3. The Control Strategy

In order to design a stabilizing controller for the flexible system it is convenient to resort to the standard TITO (two-input-two-output) control problem framework depicted in Figure 1. In this figure, is the system to control, with measured output , control input , generalized disturbance , and performance output . Augmenting the flexible system with suitable fictitious “disturbance’’ and “performance’’ matrices, the following generalized plant is obtained: where is a fictitious disturbance, a fictitious performance output, and with and being the weighting matrices to be selected. The matrices and are and real matrices, respectively, with the following structure: where the matrices and are to be selected. Finally, for the sake of simplicity, the matrices and are chosen as to satisfy the following conditions: Although not essential, conditions (13) are standard in or control [6]. It is possible to remove them by slightly increasing the complexity of the control strategy. A possible selection is where

Finally, let denote the lower LFT (linear fractional transformation) [6].

Before computing the controller, a discussion on the required properties of the controller is in order. As stated in the introduction, stable bandpass controllers are preferred options in the control of flexible systems, due to their ability to filter out accelerometer biases. The following lemma is useful for characterizing the state space representation of a general class of bandpass systems.

Lemma 1. *Let be a linear time invariant (LTI) -(McMillan) degree system with minimal state space representation
**
where
**
where all the matrices are partitioned assuming the state to be split into two -dimensional substates and the square matrix is invertible. Then, the system is bandpass, in the sense that
*

*Proof. *Although the lemma is an extension of Theorem 1 in Cavallo et al. [3], it can be easily proved observing that
Thus, direct computation shows that the static gain is zero. The high-frequency behavior is trivial.

The above lemma has a key role in our context. Indeed, it guarantees that a controller with a bandpass frequency shape, as required in the introduction, can be imposed by satisfying the rather mild conditions (18). Thus, the lemma is useful to generate different families of bandpass controllers.

The last issue to address in order to design a useful controller for flexible structures is stability of the controller itself. It is well known that inserting unstable elements into the control loop increases sensitivity to disturbances and calls for high-bandwidth controller, thus also increasing the effect of measurement noise [7]. Moreover, in the case of flexible structures the high-bandwidth requirement conflicts with the need to avoid spillover effects.

A first possible family of controllers is obtained in the case of square plants, that is, when . Indeed, in this case it is possible to compute the explicit solution of a stable stabilizing controller satisfying an index. Preliminarily, the following lemma is needed.

Lemma 2. *Consider the system (6)–(8), assume, , and let and be positive semidefinite matrices satisfying
**
for a given . Then, a necessary condition for the solution of the AREs (21) is to be of the form
**
where is a real scalar, that is, and are block diagonal.*

Although rather technical, Lemma 2 is useful to characterize the closed-form solution to an control problem, as stated by the following theorem.

Theorem 3. *Consider the system (6)–(8), with matrices given by (2), (10), and (11), assume, , let the pair be stabilizable and the pair detectable, and let
**
where
**
and is a scalar parameter.**Then, the following controller is strongly stabilizing and has zeros at the origin
**
where
**
and
**
Moreover, if
**
the LFT satisfies
**
Finally, and the controller
**
where denotes the pseudoinverse, achieves the optimal value
**
where is the set of stable, stabilizing controllers with zeros at the origin .*

The proofs of Lemma 2 and Theorem 3 can be found in Cavallo et al. [5]. They are based on the classic solution to the problem involving two coupled Riccati equations [6], while the lower bound for is obtained by resorting to the theory of singularly perturbed systems. It is worth noticing that the above theorem yields as a byproduct the closed-form solution of the problem, obtained by letting .

Corollary 4. *With the same hypotheses as those in Theorem 3, the controller
**
minimizes the norm .*

*Proof. *It is sufficient to let in the controller in Theorem 3.

Finally, by resorting to the equations for the parameterization of all the stabilizing controllers [8], it is possible to give a general expression for the whole family of optimal controllers and suboptimal controllers as follows.

Let be a stable transfer matrix with . Then, the family of controllers such that is where

Analogously, an (simpler) expression for the suboptimal controller may be obtained as follows: where, it has been, obviously, assumed that .

In the general case of the above formulation does not hold any longer, and the design of the controller has to be redefined. The main reason is that the Riccati equations to be solved for the problem are strongly coupled. The choice of a stabilizing controller with the structure in Lemma 1 can be carried out by exploiting the following theorem.

Theorem 5. *Any system of the form (6)–(8) with matrices given by (2), (10), and (11) is closed-loop stabilized by a controller with state space representation
**
where
**
where and are diagonal matrices with nonnegative entries.*

*Proof. *Let
It is clear that both and are positive semidefinite for any nonnegative matrices and . Moreover, by defining the matrices
it is easy to show that and , where, for a given Hamiltonian matrix , denotes the *stabilizing solution* of the Riccati equation associated to the Hamiltonian matrix [6], and the Hamiltonian matrices and are
Finally, the matrices and are
Thus, for any matrices and the controller (38) solves an problem and is therefore stabilizing.

Controller stability is not guaranteed by the above theorem; thus it must be enforced as follows. Preliminarily, define From the above formulation it is easy to deduce the following theorem.

Theorem 6. *Consider the system with minimal state space representation
**
where the matrices , , and are given in (38). If the diagonal matrices and are chosen so as to satisfy the LMI
**
then the system strongly stabilizes the system (6)–(8) and is bandpass.*

*Proof. *The theorem can be easily proved by considering the Lyapunov function
The rest of the theorem follows from Theorem 5 and Lemma 1.

The above theorem shows that a parametrization of stable, stabilizing, and bandpass controller is obtained, with free parameters, namely, the diagonal entries of the matrices and . Since these matrices define the controller gain, a further constraint in the form of LMI where is the trace of the matrix and is an upper bound on the controller gain depending on the problem, can be imposed. Moreover, also the matrices and have still to be selected. Expressions (39) and (40) can suggest how to choose the degrees of freedom. According to them, and can weigh modes on which to operate, while the matrices and define the authority of the control.

The above considerations can be summarized in a procedure for controller design, namely,

(1)preliminarily, collect the diagonal elements of and in the vector , define a vector of performance weights , whose th element is the weight of the th mode,(2) solve the optimization problem where and is the -dimensional column vector with all elements equal to 1;(3)compute the controller according to (37)–(38).In the above procedure, the scalar is chosen to limit the controller gain since (54) can be written as Obviously, in the above algorithm , the th element of can assume any real value. For instance, choosing , , and would pose different penalties on different modes, thus focussing the control action more on mode one than on mode two and on mode three. Cavallo et al. [9] showed that the optimal choice of the input scaling matrices is performed based on the controllability Gramian of the open-loop and closed-loop systems. However, it is clear that the efficiency of the control action depends on both the controllability and the observability Gramians. As an example, a well-known strategy for model order reduction [10] consists of the computation of a balanced realization and a reduction based on neglecting the modes associated with the Hankel singular values, that is, square roots of the eigenvalues of the product of the controllability Gramian and the observability Gramian, with the lowest magnitude. This kind of consideration leads to the idea of selecting the modal performance weights in (50) based on the Hankel singular values. Another aspect of the presented control strategy to highlight is the selection of the matrices and, even more importantly, in (15) and (16), respectively. The former is responsible for the authority of the control. The latter is responsible for the relative weights of colocated measurements with respect to noncolocated ones. By varying the magnitude of relative to the case of noncolocated feedback can be considered with stability guaranteed.

#### 4. Experimental Results

Control procedures described in the previous section have been experimentally tested on a very complex flexible structure. In particular, the selected test article is a fuselage skin panel of a BOEING 717 (depicted in Figure 2). The panel is stiffened by two bulkheads and three orthogonal stiffeners riveted to the rear of the panel itself. In order to reproduce the conditions usually adopted in the structural testing phase, the panel has been suspended by a couple of soft springs to simulate free boundary conditions. In the experiments presented in this section piezoelectric ceramic actuators have been used both for producing the primary disturbance field and for generating the control forces on the structure aimed at counteracting the effects of the primary field.

##### 4.1. Square Case

In detail, the control inputs are produced using three of the four piezoelectric patches bonded on the panel surface. The output acceleration is measured by means of the three colocated accelerometers placed on the rear part of the panel behind the piezoelectric actuators. The structure has a high modal density in the considered frequency range Hz and a model with 148 modes has been identified as proposed by Cavallo et al. [11]. The experiments were performed by using a dSPACE rapid prototyping real-time control system with 16-bit A/D channels and 14-bit D/A channels at a sampling frequency of kHz.

A fourth piezo is used only as a disturbance input so as to introduce a disturbance contribution not belonging to the range of the input matrix in (6).

The effectiveness of the control strategy has been evaluated comparing the frequency response functions (FRFs) from disturbance input to acceleration outputs and in terms of an energy-based index defined as follows. First, the energy of the measured acceleration in the given frequency range has been estimated for the th of the five used sensors as where is the acceleration spectrum measured at the th point at the frequency , a set of linearly spaced frequencies within . Then, the relative energy reduction at the th point has been quantified as where the subscripts “ol’’ and “cl’’ stand for “open loop’’ and “closed loop,’’ respectively.

In a first experimental campaign, the optimal and controllers are tested. One of the known drawbacks in the application of optimal feedback control to vibration reduction problems is that the control action may focus only on the performance objectives, reducing vibrations only at the controlled points, while increasing vibration levels at other points of the structure. In order to check that the controller actually has a beneficial effect on the *whole* structure a matrix of 35 points evenly spaced on the panel has been considered, as shown in Figure 2.

A band-limited white noise in the range Hz has been chosen as primary field, to excite all the system natural modes in the frequency range . The noise has been applied to the structure on all four input channels so that the disturbance has components both in the range of and off range. Figures 3 and 4 present results for both and cases. The controller has been designed fixing in (25). The figures show the acceleration spectrum measured in a control point and in a performance test point. Obviously, performance is better at the control points than at the test points, where a vibration reduction is still apparent, anyway.

Figures 5 and 6 report the values of the relative energy reduction index as defined in expression (57) at each grid point for the and cases, respectively, showing the effectiveness of both controllers on the whole panel. The average vibration reductions are with the controller and 55.1% with the controller almost uniformly on the whole panel, despite the small number of control actuators, the high modal density, and the broadband nature of the disturbance.

Figures 7 and 8 further clarify the effectiveness of the proposed controllers reporting the values of the energy index in all the grid points with the same color scale in order to compare open loop with closed loop.

**(a)**

**(b)**

**(a)**

**(b)**

Finally, Figure 9 presents the time history of a control input (piezoelectric actuator voltage) in both the and the cases, showing that the control voltage with the controller is generally higher than the case (this explains why the controller behaves better). Nevertheless, in both cases the input voltage is far below the saturation limits of the actuators (V).

**(a)**

**(b)**

##### 4.2. Rectangular Case

The next campaign of experiments addresses the case of more measurements than controls available. In this case, a fifth accelerometer is used to measure the acceleration in a performance point far from both control and disturbance inputs. Moreover, the noise has been applied to the structure *only through the fourth piezo actuator* not used as control input, thus the disturbance is totally out of the range of the input matrix . Although this makes the two sets of experiments not comparable, this choice allows better highlighting of the feature of the controller when used with additional sensors located near the entry point of the primary field. Since the controller computation depends on the choice of modal performance weights and on the number of feedback sensors, the guidelines described at the end of Section 3 have been followed. As a consequence, three controllers have been computed and compared:A: colocated controller with modal performance weights based on the Henkel singular values,B:rectangular controller (with an additional sensor) with modal performance weights based on the Henkel singular values,C:noncolocated controller with modal performance weights based on the Henkel singular values.

All these controllers use the same three piezoelectric actuators as control inputs and the four accelerometers as measured outputs and are designed by selecting the modal performance weights in (50) according to the Henkel singular values. However, controller A practically uses, as output accelerations, only the measurements provided by the three colocated sensors, since weighting matrices are selected as , , and . The large value of makes the information of the fourth sensor meaningless for the control action. Controller B takes into account the problem of vibration reduction when information on the location of the application point of the disturbance force is available. Of course, in such a case, the optimal solution would be to locate an actuator/sensor pair in this entry point. However, assume that mounting the actuator in such a location is not possible for some reason, while placement of the sensor is still allowed. Under this assumption, controller B uses the information provided by the fourth sensor by selecting the weighting matrices , , and , where a lower value for has been set. The controller C has been considered to assess the behaviour of the active vibration control system when the application at hand does not allow to place actuators and sensors in colocated pairs. In particular, it practically uses only the signal of the fourth sensor placed in the disturbance entry point, by selecting the weighting matrices , , and , where a large value for has been set. The rectangular noncolocated controller C has no information on what happens to the structure in the control points.

The performance obtained with the controllers above are evaluated both in terms of vibration reduction quantified according to the index in (57) and by comparing the FRFs from disturbance input point to a control point, the disturbance point, and the performance point. In the first experiment controller A, used as reference, has been compared to controller B. Table 1 shows that the reduction energy index for controller B is significantly higher in the disturbance and performance points. Moreover, the effects of controller B, compared with controller A and open loop, on the entire frequency range are shown in Figure 10, where over dB of vibration reduction peaks are visible and no spillover appears.

**(a)**

**(b)**

**(c)**

The controllers have been also compared in terms of required control energy according to the index defined below. First, the energy of a generic discrete-time signal has been evaluated as where is a fixed number of samples equal for all experiments. Taking into account that the primary field has been randomly generated, the energy of the th control signal in the experiment ( denotes the control label as listed above) has been normalized with respect to the actual energy of the disturbance in the same experiment, that is, The control energy required by controllers B and C has been then compared to the one required by controller A, by considering the control energy reduction index The results are summarized in Table 2, which shows that controller B requires less control energy than controller A despite its higher performance. It is apparent that the information from the additional sensor improves the control action.

In a subsequent experiment controller C has been tested. Although in the design phase the weighting matrices have been selected positive definite according to (15) and (16), in the real experiment, only the fourth sensor placed in the disturbance point has been physically used, in order to demonstrate the possibility to use the proposed controller as a true non-colocated one. The performance in terms of reduction energy index is reported in Table 3. The reduction energy index values are lower than in the previous cases, but the same applies to the required control energy indicated in Table 2. Interestingly enough, even in this actually non-colocated control no spillover phenomenon occurs, as evident from the FRFs reported in Figure 11, still obtaining over dB of vibration reduction peaks.

**(a)**

**(b)**

**(c)**

#### 5. Conclusions

The paper reported a procedure for designing closed-form optimal ( and ) controllers for large-scale square systems, together with an extension to nonsquare plants to tackle the vibration reduction problem of flexible structure with the strong requirement of stable and bandpass controller. The original control law is here extended to the case of more sensors than actuators. The more complex design is tackled by solving an LMI problem, which allows also to select the performance weights based on controllability/observability indices. The experimental results obtained for a vibration reduction problem of a stiffened aeronautical panel controlled by piezoelectric actuators confirm the applicability of the strategy in actual systems and its advantages with respect to a controller assuming only colocated actuator/sensor pairs.