# Vibration Attenuation in Rotating Machines Using Smart Spring Mechanism

Authors: Aldemir Ap. Cavalini; Thiago Vianna Galavotti; Tobias Souza Morais; Edson Hideki Koroishi; Valder Steffen
Journal: Mathematical Problems in Engineering
 (2011)

---

## Abstract

This paper proposes a semiactive vibration control technique dedicated to a rotating machine passing by its critical speed during the transient rotation, by using a Smart Spring Mechanism (SSM). SSM is a patented concept that, using an indirect piezoelectric (PZT) stack actuation, changes the stiffness characteristics of one or more rotating machine bearings to suppress high vibration amplitudes. A Genetic Algorithm (GA) optimization technique is used to determine the best design of the SSM parameters with respect to performance indexes associated with the control efficiency. Additionally, the concept of ecologically correct systems is incorporated to this work including the PZT stack energy consumption in the indexes considered for the optimization process. Simulation carried out on Finite Element Method (FEM) model suggested the feasibility of the SSM for vibration attenuation of rotors for different operating conditions and demonstrated the possibility of incorporating SSM devices to develop high-performance ecologic control systems.


---

## Body

## 1. Introduction

Commonly, rotating machines cross critical speeds during their transient rotation leading the system to undesirable vibration amplitudes. Under this condition, the occurrence of catastrophic failures caused by crack propagation due to the fatigue process is intensified. Therefore, significant research effort has been devoted to the development and improvement of mechanisms capable to attenuate undesirable vibrations in rotating machines [1–3].The control approaches for rotating machines are clustered into three main categories, namely passive, active, and semi-active techniques. Passive techniques are normally performed by devices known as absorbers or isolators. These techniques perform typically over a limited frequency bandwidth and, consequently, are unable to adapt their characteristics to changes in the system. Differently, active approaches promise vibration suppression over a broadband of frequencies in which the suppression is performed by incorporating active actuators, such as PZT stacks, magnetic bearings, and electromagnetic actuators to the machine to act directly against the vibratory loads. Unfortunately, successful implementation of these approaches has been limited by displacement capabilities of the piezoelectric actuators and the expensive costs of magnetic bearings [4]. The semi-active techniques represent an alternative solution to these problems. In semi-active approaches, the vibration is attenuated through an indirect manner by changing the structural parameters of the machine, such as damping and/or stiffness. In our days, the implementation of this technique in rotor dynamics is made possible by techniques such as magnetorheological and electrorheological dampers and SSM.SSM uses an indirect PZT stack actuation to change the stiffness characteristics of the rotating machine to attenuate high vibration amplitudes. This mechanism can be implemented to attenuate a wide range of vibratory loads such as axial, bending, torsion, or a combination of these ones [4]. For bending control purpose, which is the case studied in the present paper, the SSM must be orthogonally arranged in a plane located at the one or more bearings of the rotating machine.However, this paper presents a numerical simulation to evaluate the effectiveness of the SSM to attenuate high vibration amplitudes of a rotating machine passing by a critical speed during its transient rotation. The system is composed by horizontal flexible shaft with two rigid discs and three bearings. The SSM parameters are optimized using a GA multiobjective optimization technique, so that the objective function includes the simultaneous minimization of the norm and the maximum absolute value of the outputs measured at all the bearings together with the PZT stack energy consumption, in order to improve the control performance. The inclusion of the energy consumption in the minimization process aims at considering a new worldwide trend in machine design, that is, to obtain control systems that are ecologically correct. This relevant concept has attracted the attention of many researches, being Skladanek et al. [5], Guyomar et al. [6], Matichard and Gaudiller [7], and Maslen et al. [8] the examples of recently published papers involving ecologically correct control systems design, in some of them, for rotating machines, namely ecorotors.


## 2. Rotor Model

The FEM model of a flexible rotor in transient motion is represented by a matrix differential equation that describes the dynamic behavior of the system(2.1)[M]{ẍ}+[D+∅̇G]{ẋ}+[K+∅̈Kst]{x}={Fext},
where M is the inertia matrix, D is the damping matrix, G is the gyroscopic matrix, K is the stiffness matrix, Kst is a stiffness matrix resulting from the transient motion, x is the generalized displacement vector, Fext is the external forces vector, and ∅̇ is the angular velocity of the rotor [9].In this paper, the shaft was modeled by using the Timoshenko’s beam element with two nodes and four degrees of freedom per node, two displacements, and two rotations. Due to the size of the matrices involved in the equation of motion, the pseudomodal method was used to reduce the dimension of the FEM model. For this aim, the reduction is achieved by changing from the physical coordinatesxto modal coordinates q as follows:   (2.2)x=Φq,
where Φ is the modal matrix containing the m first vibration modes of the nongyroscopic, symmetric and undamped associated rotor.Substituting (2.2) into (2.1) and multiplying the resulting expression by ΦT, the reduced equation of motion of the rotor is given by (2.3)[M̃]{q̈}+[D̃+∅̇G̃]{q̇}+[K̃+∅̈K̃st]{q}={Q},
where(2.4)[M̃]=ΦT[M]Φ,[D̃]=ΦT[D]Φ,[G̃]=ΦT[G]Φ;[K̃]=ΦT[K]Φ,[K̃st]=ΦT[Kst]Φ,[Q]=ΦT[Fext].The solution of (2.3) results in a response vector described in modal coordinates. By applying (2.2), it is possible to convert the dynamic response to physical coordinates.


## 3. Smart Spring Mechanism

Several applications of the SSM have been explored for vibration suppression in the last years. Daley et al. [10] developed a system based upon an electromagnet combined in parallel with passive elements for vibration suppression in marine structures, namely, the smart spring mounting system. The results demonstrated that only the rigid body modes of the machinery were controlled. Yong et al. [11] obtained successful results using a different SSM-patented concept based on two springs arranged in parallel. In their paper, it was demonstrated the concept ability for multiple harmonics vibration control of helicopter blades. Nitzsche et al. [12] presented the control efficiency of the same concept through numerical and experimental investigations. Furthermore, a feedback control system was used to improve the control efficiency. The architecture of the SSM used in the last two presented papers is conceptually shown in Figure 1.Figure 1
Smart spring concept (adapted from [4]).If the friction between the PZT stack and the host structure is disregarded (normal forceN(t) = infinite and the friction f(t)=0), the SSM behavior can be represented by two distinctly dynamic systems, according to the PZT stack status. The first one is active when the PZT is turned off, and the equivalent SSM stiffness is given by kps+kas. The second status is obtained when the PZT is turned on. In this situation, the active spring is unattached leading the SSM to operate only through the primary spring kps. Note that if the friction is considered, the SSM becomes able to vary combinations of stiffness and damping at a particular location of the structure. Thus, the SSM mechanisms can actuate over a broadband of frequencies with low voltage requirement and displacement capability of the piezoelectric actuator [4]. It is worth mentioning that most semi-active vibration control approaches are capable to change only a single structural property.


## 4. Applications

The proposed methodology was numerically applied to a rotor system composed of a horizontal flexible steel shaft, represented by 20 Timoshenko’s beam elements, two rigid steel discs (D1 and D2), and three asymmetric bearings (Figure 2). The physic and geometric properties of the shaft, discs, and bearings are given in Table 1. The matrix equation of motion of the studied rotor was solved by using a MATLAB/SIMULINK code. The obtained responses were compared with the commercial software ROTORINSA for validation purposes. In all analyses performed in this study, only the displacement responses generated by the first ten vibration modes of the rotor measured along the x and z directions at the bearings positions were considered.Table 1
Physics and geometric properties of the rotor elements.
Rotor elementsPropertiesValuesShaftLength (m)0.588Diameter (m)0.010Young’s modulus (Pa)2.0× 1011Density (Kg/m3)7800DiscD1Thickness (m)0.005Diameter (m)0.100Young’s modulus (Pa)2.0× 1011Density (Kg/m3)7800DiscD2Thickness (m)0.010Diameter (m)0.150Young’s modulus (Pa)2.0× 1011Density (Kg/m3)7800BearingsB1, B2 and B3kxx (N/m)49.0×103kzz (N/m)60.0×103Dxx (Ns/m)5.0Dzz (Ns/m)7.0Proportional dampingD = αM+ βKα1.0×10-1β1.0×10-5Figure 2
Rotor model.The proposed SSM on-off control strategy is based on instantly reductions in the stiffness of the bearings when the rotor gets close to a critical speed during its transient rotation. This procedure is performed in order to change the position of the critical speed in the Campbell diagram and, consequently, to attenuate the vibration amplitudes generated at that particular rotation speed. For this aim, in this paper, two smart springs were installed along thex and z directions at the location of the bearing B3(Figure 3)—the highest outputs amplitudes were measured at B3—to control the vibration generated by the rotor in transient rotation from 0 to 3000 RPM, during 3.5 seconds. In this application, the target section of each SSM (see Figure 1) was attached to the ground, and only the stiffnesses kxx and kzz of the bearing B3 were used as varying parameters, being the friction between the PZT stack and the host structure disregarded. Additionally, for design purposes the equivalent stiffness of each SSM, namely, kxps+kxas related with the SSM installed along the x direction, and kzps+kzas related with the SSM installed along the z direction, are equal to the stiffnesses kxx and kzzof the bearing B3, respectively.Figure 3
SSM locations atx and zdirections of the bearing B3.Figure4 shows the regions where the PZT stack of the SSM installed along the z direction is either turned off, thus producing an equivalent stiffness associated with the zdirection of the bearing B3(kzps+kzas), or turned on, so that the equivalent stiffness is simply kzps. The same strategy is applied to the SSM along the x direction. Note that for the output amplitudes between 1.0×10-3 m and   -1.0×10-3 m (outputs measured along the x and z directions of the bearing B3) the PZT stack of each direction is in the off status. Outside of this range, the region delimited by the two lines, the PZT stack remains on. An important remark is that the choice of the SSM operation region is defined by the user, considering some key points such as, for instance, the desired vibration attenuation and the available electric power for the PZT stack.Figure 4
Equivalent SSM stiffness along thez direction related with the output amplitude.The best choice of the SSM parameterskxps and kxas for the x direction and kzps and kzas for the z direction, which are associated with the control performance, was determined by using a multiobjective optimization procedure, as performed by the so-called compromise programming (CP). The unconstrained minimization of the scalar objective function was achieved by using a GA technique with an initial population of 200 individuals and 40 generations.According to Vanderplaats [13], CP is able to combine various objective functions in order to obtain a reasonable compromise solution to many objectives. The compromise objective function is given by(4.1)F(x)={∑k=1K[Wk{Fk(x)-Fk*(x)}Fkw(x)-Fk*(x)]2}1/2,
where F(x) is the CP objective function, Wk is the weighting factor assigned to the kth objective function, Fk(x) is the kth objective function, Fk*(x) is the kth objective function target, and Fkw(x) is the worst known value of the kth objective function. In this paper, the Fk*(x) was determined by the optimization of each objective function considered independently, and Fkw(x) was associated to the worst configuration. CP was performed combining the norm and the maximum absolute value of the responses measured at the bearings positions, measured along the xand zdirections, and the total number of activation times, Nb, which each PZT stack is turned on. The weighting factors Wkrelated with the two first indexes were chosen as equal to 1.0 and the third one was chosen as equal to 1.5.These indexes were minimized considering the design variable Px and Pz (0<Px<1 and 0<Pz<1), being kxas=Pxkxx and kzas=Pzkzz. Table 2 presents the SSM parameters found in the minimization process.Table 2
Results obtained by the minimization processes.
DirectionObjectivePkas (N/m)kps (N/m)xNorm0.639531335.517664.5Maximum0.626030674.018326.0Nb0.958946986.12013.9CP0.650731884.317115.7zNorm0.809548570.011430.0Maximum0.483128986.031014.0Nb0.910154606.05394.0CP0.705242312.017688.0Note that the results for CP show that the stiffness reduction provided by the SSM installed along thex direction was as big as 65.07% and 70.52% along the z direction, abrupt reductions that could generate instability in the rotor motion. However, Figure 5 shows through the outputs measured along the x and z directions at the position of the bearings that the rotor keeps stable despite the large stiffness reductions imposed to the rotor. This behavior can be explained by the amount of damping in the system. Additionally, Figure 5 compares the dynamic responses for the rotor without and with the SSMs. The control performance appears to be efficient along the both directions at all the positions where a bearing is placed to support the rotor. This means that the SSMs are efficient not only at the point to which the SSMs were attached (bearing B3).Dynamic responses calculated for the rotor without and with the SSM.
(a)
Output measured along thex direction of the bearing B1(b)
Output measured along thez direction of the bearing B1(c)
Output measured along thex direction of the bearing B2(d)
Output measured along thez direction of the bearing B2(e)
Output measured along thex direction of the bearing B3(f)
Output measured along thez direction of the bearing B3In order to demonstrate quantitatively the SSM control efficiency, Table3 shows the RMS value (Root Mean Square) of the outputs obtained for the rotor without and with the SSMs. One can note that the vibration attenuation is efficient along both directions of the bearing B3 (greater than 30%), as shown in Figure 5. Expressive vibration attenuation is observed also at the bearings B1 and B2.Table 3
SSM control efficiency (response RMS values).
Directions/bearingsWithout SSMs  (10-3)With SSMs (10-3)Difference (%)x/B10.12850.0984-23.42x/B20.38180.2255-40.93x/B30.41090.2581-37.17z/B10.13690.1008-26.38z/B20.45760.2919-36.21z/B30.47490.3210-32.41Another aspect of the SSM technology presented in this paper is related to the design of ecologically correct machines. RemovingNb from the optimization process, the PZT stack of the z direction is turned on 50 times. When the energy consumption is considered in the process (case presented in Figure 5 and Table 3), this number decreases to 36 (28% of reduction). For the x direction, this number remained approximately constant (reduction from 34 to 33 times).Additionally, supporting the adoption of the ecorotor concept in the design of machines and systems, Table4 shows the RMS values of the rotor response. This means that including the energy consumption in the optimization process does not alter significantly the efficiency of the control. However, the amount of energy consumed by the system is reduced for this configuration.Table 4
SSM control efficiency of the ecorotor.
Directions/bearingsEcorotor (10-3)Norm and maximum (10-3)Difference (%)x/B10.09840.09661.8394x/B20.22550.2266-0.4750x/B30.25810.2595-0.5051z/B10.10080.09525.9292z/B20.29190.27326.8444z/B30.32100.30714.5213Previous results showed the efficiency of the SSM with respect to vibration reduction of a rotor that crosses critical speeds during its transient rotation. However, rotating machines are also susceptible to external disturbances, such as impact. Therefore, it is necessary to evaluate the SSM influence for this kind of excitation. Figure6 shows the behavior of the rotor subjected to impacts applied along the x and zdirections at the position of the disc D2 when t =1.75 sec. The dynamic responses of the rotor with the SSMs are then presented in Figure 6. The rotor is accelerating from 0 to 3000 RPM, similarly to the previous case. It is possible to observe significant differences between the outputs obtained with and without impact (refer to Figure 5).Output responses measured for the rotor with and without impact during the transient rotation.
(a)
Output measured along thexdirection of the bearing B1(b)
Output measured along thezdirection of the bearing B1(c)
Output measured along thexdirection of the bearing B2(d)
Output measured along thezdirection of the bearing B2(e)
Output measured along thexdirection of the bearing B3(f)
Output measured along thezdirection of the bearing B3However, one can see that the rotor motion is kept stable when the SSMs are operating, which means that the behavior of the system with the SSMs is robust to external disturbances during the transient motion of the rotor. The same evaluation was performed for the rotor under steady state condition at 3000 RPM as depicted in Figure7. Here again, impacts were applied along the xand z directions at the position of the disc D2 when        t =1.75 sec. It is possible to observe significant differences between the outputs obtained with and without the impacts disturbances. However, as in the previous cases, the rotor motion remains stable when the SSMs are operating due to the amount of damping in the system. This indicates that the system with SSMs operating under steady-state condition is robust to external disturbances.Output responses measured for the rotor with and without impact for steady state rotation.
(a)
Output measured along thexdirection of the bearing B1.(b)
Output measured along thezdirection of the bearing B1.(c)
Output measured along thexdirection of the bearing B2.(d)
Output measured along thezdirection of the bearing B2.(e)
Output measured along thexdirection of the bearing B3.(f)
Output measured along thezdirection of the bearing B3.Finally, the influence of the SSMs was evaluated for the rotor at rest. Figure8 shows the output measured when the rotor was excited by impacts applied along the x and z directions at the position of the disc D2 when           t =1.75 sec. Note that when the SSMs are operating, the maximum output amplitude is 2.0 × 10-3 m (Figure 8(e)).Output responses measured for the rotor at rest excited by an impact force.
(a)
Output measured along thexdirection of the bearingB1.(b)
Output measured along thezdirection of the bearing B1.(c)
Output measured along thexdirection of the bearing B2.(d)
Output measured along thezdirection of the bearing B2.(e)
Output measured along thexdirection of the bearing B3.(f)
Output measured along thezdirection of the bearing B3.However, one can conclude that the difference founded between the cases with the rotor in rotation and the one with the rotor at rest is associated with the damping.It is worth mentioning that other ranges can be used to activate the smart spring (see Figure4). For example, if the value of the amplitude response is reduced to 0.5×10-3m the optimal stiffness values result very much, that is, the system becomes very flexible (the first critical speed is reduced significantly). Also, the strategy to cross the critical speed safely could be reached by increasing the stiffness of the system so that the first critical speed would move to a higher rotation. As this is an exploratory study, the results shown aim at demonstrating the potential use of smart springs in rotor dynamics.


## 5. Conclusions

In this paper, a vibration attenuation technique based on SSM for rotating machines that cross critical speeds during its transient rotation was evaluated. The proposed approach presented promising results, leading to reductions as large as 30% for the case in which the parameters of the SSM were optimized. The optimization process was conducted using a GA technique that was applied to a multi-objective problem. Despite the stiffness reduction (around 65%) imposed to the rotor by the SSM, mathematical simulation results indicate that the system did not become unstable during the test. The robustness of the approach with respect to external disturbances was analysed. First, an impact force was applied to the rotor disc during the transient rotation. In the second case, the rotor was evaluated for steady state conditions (constant rotation) and, finally, the system vibration was analysed for the rotor at rest. For the three cases considered, the results showed that the on-off control technique is robust to external disturbances. Additionally, it was possible to conclude that the inclusion of the concept of ecorotors in the design of the controller did not reduce significantly the SSM efficiency while reducing the amount of energy in the controller. The proposed approach seems to be very effective to reduce vibrations for different operating conditions of the rotor. Further work includes the modelling of the friction of the PZT stack inside the SSM so that the stiffness and damping of the system can be altered simultaneously. Also, the technique will be implemented experimentally.



---

*Source: 340235-2010-11-08.xml*