Vibration Control of Double-Beam System with Multiple Smart Damping Members

A controlmethod to stabilize vibration of a double cantilever systemwith a set of smart damping blocks is designed and numerically evaluated.The externally controlled magnetorheological sheared elastomer damping block is considered, but other smart materials can be used as well.The robust bang-bang control law for stabilization the bilinear system is elaborated.The key feature of the closed loop controller is the efficiency for different types of initial excitement. By employing the finite element model, the performance of the controller is validated for strongwind blow load and concentrated impact excitement of the particular point of one of the beams. For each of the excitations, the closed loop control outperforms the optimal passive damping case by over 27% for the considered energy metric.


Introduction
Passive damping solutions are appreciated for simplicity and robust performance, but they may be insufficient when resistance to incidental extreme loads is considered. When introducing additional damping to improve the system, it is pragmatic to incorporate some vibration control possibility. Semiactive techniques allow adjusting damping parameters of the system while being easy to operate without high power consumption. Application of smart materials for adaptivepassive or semiactive operation regimes is recently in favour of scientists. In these cases, damping forces are developed due to the variation of the system motion and not due to direct application of the control force as in fully active solutions. The system requires less power to operate and is less susceptible to instability and malfunctions [1].
In this paper, the authors focus on the application of control regime, which is based on mistuning the system motion according to a state-feedback control law. The strategy aims at absorbing energy from one characteristic mode and releasing it at a particular instant of the vibration, increasing higher modes of the double-beam system. A stepwise parametric modification of the system allows achieving the desired level of vibration damping and is often referred to as the switching times control [2], adaptive elasticity [3], or prestress-accumulation-release [4]. Usually, the switching operates between two extreme control values, as in the bangbang type of control, or between two states of the actuator, as in an on-off logic [5,6].
When beam systems are considered, different smart materials may be employed to perform semiactive damping strategy. The switchable stiffness is frequently obtained with piezoelectric patches. Ramaratnam and Jalili [7] applied switched piezoelectrics to stabilize flexible beam with a moving base, intended for robotic application. The phase transformation of shape memory alloys and polymers was used in [8] for morphing of aircraft wings. Ostachowicz et al. [9] used shape memory to stabilize smart beams. In [3], the authors varied the flexural bending stiffness of a multilayered beam with sheared polymer layers. The shear 2 Shock and Vibration modulus switching of the polymer layer was obtained by heating it with ultrathin electric blankets. Recently, the present authors in [10,11] proposed a pneumatic damping member based on switching the jamming transition of the granular material. Interesting practical examples of robust switched control based on a variable stiffness was applied to civil structures by means of on-off hydraulic devices [12][13][14]. Winthrop et al. in [15] provided a good review of the most notable solutions based on the periodic switching of the system parameters.
The paper [16] devoted to the theoretical and experimental investigation of free vibrations of a layered beam with a single magnetorheological elastomer member proved the concept of switch stiffness to be effective, giving 20% improvement of the damping performance compared to the passive system. The research advances motivated us to study the performance of a more complex system, where several smart blocks are collaborating to counteract the vibrations introduced by distributed or impact loads, introducing higher structural modes.
The contribution of the work lies in the derivation, verification, and characterization of semiactive control strategy for independently operated smart blocks coupling two beams. The control logic is proposed, and theoretical effectiveness is demonstrated in numerical simulations for multidegree of freedom structure with multiple variable-stiffness members. As far as different types of smart materials can be used as actuators, we do not suggest any particular type to be controlled with the derived algorithm.

Mathematical and Numerical Model
The double-beam structure to be investigated consists of two elastic face sheets coupled by several viscoelastic smart damping elements (see Figure 1). Both beams have length , width , thickness , and mass per unit length and are mounted in a clamped support. We are going to investigate vibrations of the structure in the plane of symmetry -. The displacements of beams from their equilibria indirection are denoted by 1,2 ( , ). The beams are slender and their deflections are small, so the linear Bernoulli-Euler beam theory is valid. The external damping of the beams is neglected. The effect of gravity, which acts perpendicularly to the plane of vibrations, is not taken into account either, since the equilibrium state is taken as a neutral one.
We consider identical damping blocks of rectangular shape 2 × × 2ℎ. Their centers are located at distances 1 , . . . , from the support. The rheological parameters of these elements, namely, the storage and loss moduli, are controllable by means of external devices embedded in the structure or placed nearby. The masses of the damping blocks as well as control devices are assumed to be negligible compared to the mass of the beams. Alternatively, they might be incorporated into the mass of the beams, as the damping blocks are stuck uniformly across the length of the structure.
To be more specific, we are going to consider a magnetorheological elastomer as the material of the damping blocks. Exposed to a magnetic field, the MRE increases its storage and loss moduli. Still, other smart materials, for instance, vacuum packed granules, can be applied as well. The overall elastic and dissipative characteristics of the structure are to a large extent determined by the rheological parameters of the damping elements. Thus, modification of these parameters allows controlling the system. Note that the control changes both parameters, that is, the storage and the loss modulus.
Every element is assumed to be made of a Kelvin-Voigt material. The shear stress-strain relation is given by where , are the shear stress and angular deformation, ℎ 1 = ℎ 1 ( , ) = ( , ) is the storage modulus, ℎ 2 = ℎ 2 ( , ) = ( , )/ is the loss modulus divided by the angular frequency of vibrations , and is the magnitude of the applied magnetic field. Although the values of and depend on the frequency of vibrations, a typical change of storage or loss modulus due to the applied magnetic field yields only a slight change of this frequency. The values of and are nearly constant within such a narrow range of , so we may treat the values of both rheological parameters as independent of it. The relations between these parameters and the magnetic field are generally nonlinear. However, the external devices are able to control these parameters in a limited range only, where linear approximations are justified. Additionally, for the purpose of control design, we assume these relations to be monotonically growing.
The damping blocks undergo shear deformations when the beams are deflected. This yields bending moments as well as axial forces acting on both beams. Moreover, relative motion of the beams generates transverse forces. To develop the formulas for all these forces and moments, it turned out to be convenient to use the truss model proposed in [17]. Having assumed that the blocks are not too long (2 ≪ ) and treating the parts of the beams over the damping blocks as rigid bodies, we may reduce the forces to the pole, which is the point over the center of the block (see Figure 2) where the two-link diagonal truss replaced the damping block attached at tip ( = ). Let us present the formulas for the bending moments and forces acting on Beam 1. Assume that the angle of deflection of both beams at point is the same and equal to , and the relative displacement of the beams at amounts to Δ = 1 ( , ) − 2 ( , ). Then, we come up with the following expressions for bending moments : transverse forces : and axial forces : The influence of the axial forces appears to be of minor importance and it will be disregarded in the study. This is mainly because the axial forces act in the opposite directions; that is, when Beam 1 is compressed, then Beam 2 is stretched, and partly because these forces are too small to disturb the dynamics of the rigid face sheets.
The dynamics of the system are governed by the following coupled partial differential equations: where and are given by formulas (2)  We build FE model of the system to design the control strategy. The original computer code for numerical simulation was elaborated for this purpose. The characteristic matrices describing the finite element of length of discretized continuous system (5) were developed. In order to obtain the virtual work of the system, we multiply the motion equations of Beam 1 and Beam 2 by virtual displacements * 1 ( ) and * 2 ( ), respectively. After rearrangements, the total virtual energy in the element containing both beams is equal to the integrals as follows: where the prime and the dot denote differentiation with respect to spatial variable and time variable . As a result of the integration of the terms with distributions, (6) can be written in the following form: In the case of the Euler beam, the nodal displacements are associated with the rotation angles. Hence, the displacements approximation for both beams in the finite element are given by Subscripts L and R denote left and right node of the finite element, respectively. The nodal displacements and angles in the element are interpolated by the 3rd-order polynomials as follows: More details can be found in [18]; for example, taking into account the virtual hat-shaped function in the form * we can obtain the characteristic matrices of the considered system. The classical energy minimization leads to the following local inertia matrix M and the local stiffness matrix K of both beams:  ] , Shock and Vibration  ] .
Coupling matrix Y is presented below: where, according to the shape functions (9), Let be the number of nodes of each beam and T = [ 1 , . . . , 4 ] T be the vector of displacements and angles of the nodes of both beams. Finally, according to (7), the dynamics are governed by the following second-order differential equation: where M is the global mass matrix, K is global stiffness matrix, and Y is matrix that establishes the effect of the damping block. Global matrices M and K are assembled from elemental matrices M and K . Matrices Y are assembled in matrix places related to positions of the coupling blocks.

Control Design
For the sake of the further study, we will now introduce the state-space first-order ordinary differential representation of the dynamics (14). For state vector x ∈ R 8 defined as we introduce system matrices A (8 ×8 ) , given by Then, the system dynamics given by (14) can be represented by the following system of the ordinary differential equations: 6 Shock and Vibration Here, x 0 stands for a nonzero initial condition. For the control inputs, we assume that each is bounded by two nonnegative values corresponding to the minimum and maximum admissible magnitude of the magnetic field denoted by min and max , respectively; namely, Note that control variable appears linearly in control system (17). Moreover, due to the middle term in the right-hand side of (17), the system is classified as a bilinear [19]. By Picard-Lindelöf theorem for given function , the initial boundary value problem (17) possesses unique solution x( ). For the numerical results presented in the sequel of this work, the dynamics (17) will be solved by employing the 4th-order Runge-Kutta scheme. Before we design a controller, let us first recall the result from [1,20]. It has been demonstrated that the optimal policy to stabilize the bilinear systems as (17) is realized by the socalled bang-bang control laws, where the control parameter is switched between two extreme values of the admissible set. The result is also valid for optimal control of any bilinear system with a finite time integral objective functional that does not explicitly depend on the control parameter, as in the cases of energy-related functionals. Therefore, unless there are strong practical reasons to reduce the amount of control input, the usage of controls other than bang-bang controls has no rational argument. To determine an exact structure of the optimal bang-bang control, one needs to estimate the number of switches and the corresponding time instants (referred to as the switching times). Practically, this requires solving a two-point boundary value problem that, for the multidimensional system, is concerned with employment of significant amount of computational burden. Another inconvenience of the optimal bang-bang approach is due to the fact that the structures of the optimal controls strictly depend on the system's initial conditions. Since we are aiming at a controller that is robust for any initial excitement, we will opt for an alternative method.
In the present work, a switched control strategy will be designed based on Lyapunov's classical second method for stability. It will enable us to derive an explicit, simple, and robust state-feedback control law capable of real-time implementation with the use of accessible state information. For this purpose, let us first introduce the energy-related Lyapunov function corresponding to system (17): The time derivative of the energy function iṡ Our goal is to design controls * providing the best instantaneous energy dissipation; namely, ∀ : * ( ) = arg min ∈U̇( x ( )) , = 1, 2, . . . , .
To find a solution to problem (21), we first insert dynamical equation (17) into (20). That yieldṡ For simplicity, let us define the following state-dependent functions: For such a representation, the solution to problem (21) can be determined by solving the following mathematical programming problems: Assuming that, for all time , state x( ) is accessible, we will now consider three cases depending on the values of functions 1 (x( )) and 2 (x( )).

Numerical Study
In this section, we will perform a set of numerical experiments to validate the designed controller. Our primary goal will be to examine its robustness when the structure is subjected to different types of the impact excitations. For that purpose, the simulations will be performed for different initial condition scenarios.
For the considered structure, we assume two identical aluminum beams coupled by five identical controlled damping blocks (see Figure 3). For a such system, we select a set of parameters listed as follows. In the simulations, we will use the finite element model introduced in Section 2, where, for each of the beams, we will assume ten equal elements placed among eleven nodes. The controlled blocks are located at positions (1/10) , (3/10) , (5/10) , (7/10) , and (9/10) . The simulations will be performed for the time horizon of = 3 s. To establish the comparative results, we will utilize the total energy metric defined as follows:

Beam Parameters
Here, K and M are defined according to (14). For the extreme values of the admissible control (see (18)), we will assume min = 0 and max = 1. Control variable denotes the intensity of the field in the dimensionless form; that is, this is the actual value divided by the maximum one.
For the assumed finite element representation, we introduce I 1 = {1, . . . , 11} and I 2 = {12, . . . , 22} standing for the set of nodes associated with the lower and upper beam, respectively. Let stand for the horizontal position of th node. We will consider the following three variants of the initial condition for system (5).

Variant 1.
In this case, we will be interested in how the designed control can stabilize the structure subjected to a single earth tremor or a strong wind blow. For that purpose, we will impose the initial velocity distributed over the set of (31) In this case, we will examine the performance of the controller when the structure is subjected to concentrated impact excitement. We will impose an initial velocity to the node placed at the end of the lower beam. As before, for the upper beam, we will assume the zero initial condition.  Figure 4. We can clearly observe that the increase of the value of the control parameter results in higher rates of decrease of the vibration amplitude. The same result applies to the upper beam. For the assumed metric (30), over 46% benefit is achieved for  = 1 compared to the case of = 0. A similar conclusion is observed for the initial conditions of variants 2 and 3 ( Figures 5 and 6) where the total energy is reduced by 48% and 44%, respectively (see Table 1 for the full set of the energy metric values). The presented results confirm the expected dynamical behavior since the value of the control parameter is related to the structural damping parameters (see (1)). To examine the performance of the designed controller, for each   initial condition variant, we will use the passive strategies = 1, later referred to as the optimal passive cases.

Switched Controller versus the Optimal Passive Cases.
For the testing of the designed controller, we begin with the simulations under the initial condition of variant 1. The comparison of the deflection of the ending point of lower beam 1 ( , ) for the controlled case (referred to as closed loop case) and the optimal passive case is depicted in Figure 7.
At the first stage of the simulations, we observe that both trajectories are almost identical which is due to the fact that, for each of the blocks, the switching controller mainly operates on maximum admissible value max = 1. After reaching the first maximal deflection (at time ≈ 0.18 s), we notice a gradual divergence of the curves. After this time, we can observe a decrease of the velocity (inspected by the lower slope) in the closed loop case. This velocity decrease results in 16% reduction of the second extreme deflection. For the rest of the simulation, we clearly observe the tendency of faster stabilization for the closed loop case. Switching of the control parameter influences both the damping and the stiffness of the blocks according to (1). As the control is applied, due to the temporary stiffness loss, the peak frequencies slightly shift towards the lower values ( Figure 10). The amplitude is presented as a relative strength of the harmonic component of the source signal. For each frequency, the damping ratio was estimated using half-power bandwidth method and plotted on the frequency response graph. The damping ratio for the optimal passive case is denoted , while is used for the controlled, closed loop case. One can notice that the controlled case outperforms the optimal passive case for each variant. Still, the damping effectiveness differs.
In Figure 11, we illustrate the switching control modes of the damping blocks at time instants = 0.002 s, = 0.02 s, = 0.2 s, and = 2s. Here, for each of the blocks, red and blue colors correspond to switching modes max and min , respectively. Note that the initial condition of variant 1 enforces the approximate shapes of the beams' first mode, so single frequency peak is observed. As a result, for all of the damping blocks, the switches triggered by the sign of the functions (23) take place nearly at the same time (due to the coherence in the changes of signs of the state vector terms), and thus their switching modes remain consistent. In terms of the assumed energy metric, the closed loop case outperforms the optimal passive case by over 31% (see Table 2). As a result of the control, the damping ratio increases from = 0.164 to = 0.186, that is, by 19%.
The comparison of the deflection of the end point of lower beam 1 ( , ) in the case of the initial condition of variant 2 is demonstrated in Figure 8. At the first stage of the simulation, for both trajectories, we observe a minor presence of the higher vibration modes which is due to the point impact. At this stage, the switching modes of the damping blocks are not consistent (see Figure 12). Thanks to the control policy, the first extreme deflection of the closed loop trajectory is slightly reduced. As for the previous initial condition scenario, we observe a reduction of the velocity and the second maximal deflection for the closed loop case. Afterwards, the dynamics become dominated by the first mode, and the switching pattern returns to that observed in the case of variant 1. The overall reduction of the energy is 33% compared to the optimal passive case. When the closed loop switching control is applied, the damping ratio for the first-peak frequency increases by 22%, from 1 = 0.158 up to 1 = 0.193. When the frequency of the second mode is considered, the damping increases by the same percentage, from 2 = 0.133 up to 2 = 0.162.
In Figure 9, we present the comparison of deflection 1 ( , ) in the case of the initial condition of variant 3. As in variant 2, at the first stage of the simulation, we observe the presence of higher vibration modes resulting in irregular switching control patterns (see Figure 13). We can notice a minor drop of the peak amplitude for the closed loop case. After vanishing of the higher vibration modes, the dynamics and the control pattern return to those of variant 1. Compared to the optimal passive case, the closed loop system resulted in 28% reduction of the total energy metric. The passive damping ratio at the first-peak frequency 1 = 0.161 is increased by the switching control policy to 1 = 0.185, which is 15% growth. When the second mode is considered, the damping ratio for the optimal passive case is 2 = 0.215. As a result of the designed closed loop controller, this value was raised by 16% to 2 = 0.250.
Finally, we will analyze the character of damping in the system. We define the damping decrement as the ratio of consecutive extreme deflections of the lower beam (Note that multiple extremes may occur at the first cycles of vibrations of the dominating mode, which results from the fact that several natural modes are induced, and this is clearly visible in Figure 9. In such a case only the extremes of a dominating mode are taken into account.). Specifically, if denotes a time instant of th extreme deflection ( = 1, . . . , , where is the number of extremes), then the damping decrement at is equal to The values of damping decrements for optimal passive and closed loop cases are depicted in Figure 14. For every variant of initial conditions, the closed loop control almost doubles the decrease rate of the amplitude, from approximately 25% to 50%. For most of the simulation, the values of decrements are constant which is a consequence of the presence of linear viscous damping. The exception is observed at the beginning of the simulation. Here, the values of the initial 2-3 decrements differ because more than one of the natural  Figure 11: Activity of the controlled blocks for the selected time instants in the case of variant 1. Blue and red colors correspond to = max and = min , respectively. modes are induced. It is particularly observable in variants 2 and 3, when the initial conditions induce several modes at once, whereas in variant 1 the initial condition given in a third-order polynomial form is closely related to the first eigenfunction of a cantilever beam.

Conclusions
In this paper, a control method to stabilize adjacent structures by means of incorporated controllable damping blocks has been designed and numerically evaluated. The controller operates on two switching modes triggered by a practical and simple implementation state-feedback control law. It provides high rates of the energy dissipation and is robust to a change in the initial condition. For each of the considered scenarios, the closed loop system outperformed the optimal passive cases by over 27%. In this work, the designed controller was suggested for implementation to elastomer based damping blocks that change their mechanical properties due to the magnetic field. The analogous operation, however, can be performed by a number of smart material actuators, and it is applicable to structures as an additional system supporting the existing damping solutions if specific constructional factors, such as the uncertainty of the material properties, nonlinearity, and fail-safe design, are taken into account.
Currently, the authors are preparing an experimental setup to validate the proposed method on a real structure. For the future study, the authors address the problem of designing a distributed structure of the proposed controller. Relying on a limited amount of state information, it can benefit functionality and robustness.