#### Abstract

In the paper, we propose distributed feedback control laws for active damping of one-dimensional mechanical structures equipped with dense arrays of force actuators and position and velocity sensors. We consider proportional position and velocity feedback from the neighboring nodes with symmetric gains. Achievable control performance with respect to stability margin and damping ratio is discussed. Compared to full-featured complex controllers obtained by modern design methods like LQG, H-infinity, or mu-synthesis, these simplistic controllers are more suitable for experimental fine tuning and are less case-dependent, and they shall be easier to implement on the target future smart-material platforms.

#### 1. Introduction

The established paradigm in past and current active damping projects is as follows: the mechanical object is defined first (plate, beam, car door, wing panel, etc.). Systems detailed design and modelling phases follow the methods in [1, 2] giving rise to very accurate FEM models with tens of thousands of degrees of freedom. Alternatively for existing prototypes, the experimental identification approach can be applied to get the mathematical models directly via experimental modal analysis [3]. Model order reduction [1, 4] then gives accurate enough yet tractable models for optimal actuators and sensors placement [5–7]. Finally, a very limited number of them are considered (say up to twenty) for the design of the control laws [1, 8]. Finally, validation and verification of the solutions by high-fidelity simulations is performed, followed by laboratory experiments and final deployment of the product. For any new project—or even a relatively mild modification of a previously accomplished project—all these steps must be performed (or re-visited) again. Implications towards requested research and development costs are significant and obvious.

Therefore, there is a need to use other type of control methodologies, and recent advances in MEMS sensors and microactuators, ongoing intensive research on new smart materials, and progress in computational power pave the way to massive development of heavily distributed control in this context.

Distributed control is now a very active field of research, thanks to potential applications which require high scalability and reliability. The main advantage of using distributed control is the locality of the necessary measurement and actuation—the measurements are collected and processed in a distributed manner. This kind of control can be applied for automated highway systems [9], car formations [10], and also flexible structures. The work in [11], for instance, studies a flexible beam model with bending and torsion motions, and a distributed arrangement with two force-actuators and three moment-actuators paired with rate gyros was elaborated. In [12], a dense network of piezoelectric patch actuators was proposed to realize the distributed actuation. In [13], a distributed piezoelectric actuation was involved and applied to the placement problem of patches so that the deformations are suppressed at preselected locations. Multipositive feedback approach for flexible structure control was presented in [14]. Since the flexible systems are passive by nature, one can also employ a lot of results available for distributed control of the passive system [15, 16]. Completely passive solutions can be obtained using piezostructures, as reported in [17].

One of the natural goals when dealing with control of flexible mechanical structures is vibration suppression. One standard approach relies on application of a large number of neutralizers placed in prespecified locations along the structure composed from masses and springs. The goal is not only to design the neutralizers’ parameters but specify their locations as well since vibrations can be eliminated only at the attachment point of the vibrating beam while amplification of vibration may occur in other parts of the beam. Dynamic vibration absorbers using magnetorheological elastomers were used in [18]. In [19, 20], a set of optimum conditions for global control of the kinetic energy based on the fixed-points theory was proposed. Dynamic transfer matrices using mobility or impedance were used in [21]. In [22], an iterative procedure was developed to find the required resonance frequencies of variable stiffness neutralizers to create nodes at selected locations. Wide-band frequency passive vibration attenuation design for the absorbers was introduced in [23]. In [24], explicit model predictive vibration control was tested. A different approach consists in control and attenuation of multiple travelling waves propagated in a one-dimensional structure [25–28]. Sliding mode control on seat vibration reduction problem was applied in [29].

#### 2. Structured Control Laws for Smart Materials

The paper presents an attempt to systematic proportional decentralized position-velocity feedback for active damping of mechanical structures equipped with dense arrays of force actuators and position and velocity sensors. Such a control law is characterized by a very small number of parameters and simple procedures for their tuning compared to centralized approach. Although the results are presented for a one-dimensional structure model only, it is believed that a generalization to two-dimensional mechanical structures will be possible.

The research is motivated by vehicular platoon control where relative position and relative or absolute velocity feedback related to the preceding and succeeding vehicle is often considered [30–32]. Nevertheless, the measure of control performance in both applications is different. For vehicular platooning, the main goal consists in preserving a prescribed spacing between the vehicles and in keeping the leader’s velocity, whereas when dealing with mechanical structures, a fast and adequate damping of the oscillating modes is required. Hence, the presented control design is focused on investigation of feasible damping ratio of the least damped mode and achievable stability margin of all modes.

Throughout the paper, the superscript denotes transpose, stands for identity matrix, and denotes real and imaginary part, respectively, denotes the Kronecker product, and denotes the spectrum of a matrix.

#### 3. One-Dimensional Structure Longitudinal Model

Let us consider a one-dimensional structure composed from the masses , springs , and dampings , each of the same value. Let us assume that the input forces may act on each individual mass independently and we are able to measure positions and velocities of each mass, i.e., actuators and sensors are placed in the same positions. Longitudinal vibrations of such a structure can be described by a state-space model:wherewhere , are the positions of the masses, is the number of nodes, and is the vector of the input forces [32].

The matrices and are given bywithwhere , . The matrix is the Laplacian of the graph corresponding to the structure which is in this case given as

Since for practical reasons at least one of the nodes has to be fixed, the equation (1) becomeswhere comes from and by omitting the entries corresponding to the fixed nodes, and is the number of the nonfixed nodes.

The matrices and are then given bywhere the matrix is called the *grounded* Laplacian that results from Laplacian (5) by omitting the rows and columns corresponding to the fixed nodes. In the sequel, we will assume without loss of generality that the fixed nodes are the first and last one, , and

The eigenvalues of are all positive and are given by

We will see later that from the eigenvalues (9), the minimum and maximum ones are of special interest. Those can be determined as

In this paper, we will use distributed control law where the control action applied to each node depends symmetrically on relative positions and velocities with respect to its neighbors, i.e.,with and being fixed (Figure 1).

The control law (11) can be written asAfter substituting (12) into (6), one obtains the description of the closed-loop system:wherewith

For determination of the eigenvalues of matrix , we will use the following lemma.

Lemma 1 [32]. *Using Lemma 1, we immediately obtain the following result:It turns out that the closed-loop eigenvalues are given as the roots of the characteristic equation:i.e.,*

#### 4. Control Strategies

There are many options where to place the closed-loop eigenvalues. Nevertheless, from the vibration suppression point of view, the following two are the most interesting ones.

##### 4.1. Prescribed Damping Ratio

A quite natural option is to force all modes to be damped with a prescribed minimum damping ratio . From (19), one can see that if is fixed then with increasing value of , the least damped mode is that corresponding to . The ratio of imaginary and real part of the least damped mode that corresponds to is given byHence, the minimum value of satisfying this condition is given bywhere is set arbitrarily.

Since the corresponding damping ratio is given asafter substitution into (21), we obtain

Let us define stability margin as

Stability margin of such a control law is determined by the distance of the least and most damped closed-loop eigenvalues, and , respectively, from the imaginary axis for . The distance of complex conjugate corresponding to whose position is given by the prescribed damping ratio is given bywhereas distance of real corresponding to can be obtained as

Stability margin is then given as minimum of (25) and (26):which after some algebraic manipulations yields

It should be noted that, for very small damping , the eigenvalue is not real and the second term in (28) becomes complex and should not be considered. Nevertheless, considering such a small damping is highly impractical.

The dependence of stability margin on prescribed minimum damping ratio given by (28) for and is depicted in Figure 2.

##### 4.2. Maximum Stability Margin

Another interesting option is to find the control parameters that maximize the stability margin (24). From (19), it follows that the stability margin of the eigenvalues lying on real axis is determined by that corresponding to and stability margin of the eigenvalues lying out of real axis is determined by those corresponding to . Thus, the maximum stability margin is achieved ifwhere and are the eigenvalues corresponding to and , respectively.

Condition (29) can be written asfrom which we obtainfor arbitrarily chosen .

By substituting (31) in (29), we obtain the maximum stability margin as

The ratio of imaginary and real part of the least damped mode that corresponds to is given bythat yields after substitution from (31) and some simplificationsFrom that, the damping ratio of the least damping mode follows as

The achievable stability margins for different values of for both approaches are shown in Figure 3, whereas the minimum damping ratio corresponding to maximum stability margin is depicted in Figure 4. Both figures are plotted for .

#### 5. Example

Let us illustrate the results derived in the previous section on an example. We will consider the following parameters: , and . We set the control parameter .

The minimum and maximum eigenvalues of the grounded Laplacian become

The minimum value of damping all the eigenvalues with minimum damping ratio is given by

Such a control law guarantees the stability margin:

The position of dominant open- and closed-loop eigenvalues is plotted in Figure 5.

The value of control parameter guaranteeing maximum stability margin yieldscorresponding to stability margin:

Damping ratio of the least damped mode is given by

The corresponding position of dominant open- and closed-loop eigenvalues is depicted in Figure 6.

To demonstrate the presented design, we compare time and frequency responses of a point lying in the middle of the beam for different values of damping ratios. The open-loop responses to initial condition are shown in Figure 7.

The initial condition response for different prescribed minimum damping ratios is depicted in Figure 8. The Bode plots are compared in Figure 9.

Let us compare the achieved results with other two standard design methods typically used by the control community. At first, we design an LQ controller with relative positions and velocities considered as measurable state variables, i.e.,

The difference between the LQ and presented control law is that the LQ control law uses relative positions and velocities between all neighboring nodes and not between the closest neighbors only as in (11). To force the LQ control to use relative positions and velocities, we introduce a modified state vector:

The criterion to be minimized is then given byresulting in LQ feedback control lawwith

By tuning the weighting matrices *Q* and *R* to guarantee minimum damping ratio we obtained the corresponding stability margin , see dominant poles in Figure 10. Tuning the weighting matrices to maximize stability margin, we arrived to with corresponding that can be seen from dominant poles in Figure 11. Hence, the dominant poles configuration is very similar to the proposed design (Figures 5 and 6).

The control gains and for the former case are depicted in Figures 12 and 13, respectively. One can see that the control law uses the relative positions and velocities to the closest neighbors only and that the gains are almost the same for all nodes.

To compare our methodology with another control design approach, we formulated the task as an design for the fixed structure controller (11). It can be easily done with *hinfstruct()* function in Matlab Robust Control Toolbox. This tuning minimizes the norm of the closed-loop transfer function modeled by the closed-loop control system with tunable components and weighting filters. In our case, the high-pass filter with cutoff frequency 8 rad/s has been used to penalize all system modes.

The *H*-infinity design methodology offers efficient algorithms how to obtain multivariable control laws by specifying closed-loop frequency response requirements. This approach was used, for e.g., in [33], where authors compare classical single-input single-output controllers with *H*-infinity approaches in terms of robustness and performance. The order of the *H*-infinity control system is however equal to the so-called augmented plant containing the model of the controlled system along with the weighting filters defining performance and robustness requirements. This leads to excessively high-order control laws typically, with strong negative impact on implementation and experimental fine-tuning. For this reason, e.g., in [34], there was a method presented for the controller order reduction which is one possible way how to get control laws with reasonable complexity. Nevertheless, loss or deterioration of closed-loop performance and/or stability is often an unwanted effect associated with this approach. Thanks to recent structured H-infinity control synthesis results (see, e.g., [35, 36]), it is possible to receive the parameters of such reduced-order controllers directly, minimizing the H-infinity norm under the controller complexity constraints.

The Bode plots of original system, LQ controller tuned for prescribed minimum damping ratio , the controller, and the proposed design for are compared in Figure 14. The plots confirm that all designs give very similar results.

#### 6. Conclusion

In the paper, we presented an active approach of one-dimensional structures with dense array of collocated sensors and actuators using proportional position and velocity feedback control laws. The control law was formulated in a distributive manner, i.e., each actuator uses information from its closest neighbors only. The achievable stability margins and damping ratios were analyzed based on the properties of Laplacian matrix describing the corresponding information graph. Comparison with LQ controller and designs shows that the presented approach achieves similar results yet with much lower computational and actuator complexity.

#### Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Acknowledgments

The presented research has been supported by the Czech Science Foundation under the project No. 16-21961S.