In this paper the problem of modeling, analysis and unbalance response control of a rotor system with two disks in an asymmetrical configuration is treated. The Finite Element Method (FEM) is used to get the system model including the gyroscopic effects and then, the obtained model is experimentally validated. Rotordynamic analysis is carried out using the finite element model obtaining the Campbell diagram, the natural frequencies and the critical speeds of the rotor system. An asymptotic observer is designed to estimate the full state vector which is used to synthesize a Linear Quadratic Regulator (LQR) to reduce the vibration amplitudes when the system passes through the first critical speed. Some numerical simulations are carried out to verify the closed-loop system behavior. The active vibration control scheme is experimentally validated using an active suspension with electromechanical linear actuators, obtaining significant reductions in the resonant peak.
Rotating machinery plays an important role in modern industry due to the wide range of applications of these kind of machines (i.e. turbomachinery, compressors, generators, etc.). The presence of mechanical vibrations is an inherent phenomenon in rotating machinery. Mass unbalance and dynamic interaction between the stator and rotating parts are the main causes for vibrations. The trend in industry has been to move towards high speed, high power, lighter and more compact machinery, which has resulted in machines operating above their first critical speed and increasing the vibrations problem [1-4]. In spite of the fact that the total elimination of the mass unbalance is impossible, it is very important to control the vibration amplitudes within acceptable limits for a safe operation of the machines. To do this, machinery designers have to understand and predict the rotordynamic behavior of machinery taking advantage of the improvements in computing technologies and modeling techniques, nowadays finite element techniques are widely used to model and to analyze rotor systems (see e.g. [2,3,5,6,7]). However, we must keep in mind that any model is only an approximation of a real system and its experimental validation has to carry out. The prototypes in laboratories are very important in the study of rotordynamics because by these, the analytical models and the control schemes proposed can be experimentally validated.
Nowadays, many machines must operate to high rotating speeds which yields a supercritical operation, that is, during the running up the machine has to cross one or more critical speeds causing large vibration amplitudes and more complex dynamic effects [2,3]. Some authors have presented important results about rotordynamics analysis such as [8,9], or fault diagnosis in rotor-bearing systems [10-12], but they do not include control schemes to attenuate the vibration amplitudes. On the other hand, rotating machines increasingly incorporate transducers, actuators and control systems, in order to attenuate vibration amplitudes, to improve the machines performance and to extend their useful life. One of the most important issues of this tendency is the use of active bearings. Active vibration control has been an area of theoretical and experimental research in rotordynamics, providing many advantages for the attenuation of vibration amplitude, during run-up and coast-down through critical speeds, and minimization of sudden transient behavior due to rotor unbalance or parametric uncertainty. In literature many control strategies to reduce the unbalance response in the classical Jeffcott rotor have been published, but using very simple models [13,14]. Other authors have proposed the use of active bearings to control the unbalance response in rotor systems with more than a disk. Nicoletti and Santos [15] show a methodology to design an active control of vibrations in a rotor system using actively lubricated bearings, they employ two control strategies, in the first one the control system reduces the vibration amplitudes at the desired position of the shaft, at the expense of increasing displacement at the ends of the flexible rotor and, in the second one, significant reduction of vibration amplitude along the whole rotor is achieved, and finally, they found that the actuator must operate in a linear range and the excitation frequency must be less than the natural frequency of the servo valve, otherwise, the control system does not properly operate, the results are presented by numerical simulations. Active Magnetic Bearings (AMB) and piezoelectric actuators are other devices used to control actively vibrations in rotor systems. Bi et al. [16] propose an automatic learning control for unbalance compensation using two AMB in a small rotor of around 1.5kg obtaining an important reduction in vibration amplitudes crossing the first critical speed. Couzon and Der Hagopian [17] combine neural networks to identify and fuzzy logic to control the flexible modes of a rotor suspended on active magnetic bearings and working under the first critical speed, they obtain the system finite element model by software to design the control scheme and then, they validate experimentally this control law proving that it is possible to use artificial intelligence in the active control of a flexible structure. Lei and Palazzolo [18] present the design and analysis of a rotating system with magnetic bearings including dynamics, control and simulation results, by rotordynamic analysis they find that the first bending mode is above the maximum operating speed which permits the controller design to ignore the flexible mode; they do not present experimental results, but they show the importance of a multidisciplinary design. Simoes et al. [19] show an active vibration control for unbalance compensation using piezoelectric actuators and show that important reductions in the vibration amplitudes are possible using active bearings with linear actuators. However, AMB and piezoelectric actuators are rather expensive, exhibit physical limitations and sometimes they are not easily available (see e.g. [20]). All of the works mentioned above use a reduced order model obtained by finite element techniques, but only some of them use experimental results to validate the mathematical models and control schemes.
In this work a reduced order model for a rotor system with two disks is obtained by the finite element method. The system model includes the gyroscopic effects and it is experimentally validated and used to get the Campbell diagram, the critical speeds and the mode shapes. An asymptotic observer is designed to estimate the full vector of system states and then, with the estimated states, a Linear Quadratic Regulator (LQR) with state feedback controller is synthesized to attenuate the resonant peak in vibration amplitudes when the rotor passes through the first critical speed. Some numerical simulations to validate the control scheme proposed are presented. The experimental test rig described in [21] was used to carry out some experiments to validate the closed loop system behavior. This experimental results show important reductions in the unbalance response of the overall system.
2System modeling and analysis2.1Description and finite element model of the rotor systemFirst, we consider a rotor-bearing system shown in Figure 1. The rotor system is driven by an AC motor, which is connected to a steel shaft by a flexible couple. The shaft is supported by a traditional journal bearing on its ends. Two unbalanced disks are mounted in an asymmetrical way along the shaft, and the disks are fixed to the shaft by bridles. The motor and both bearings are fixed on a thick metallic base in order to give high stiffness to the system.
In order to obtain a finite element model, the rotor system was divided in three elements and four nodes, one in each journal bearing and another one in each disk location. It is important to mention that such elements have different length. A finite element type Euler beam is considered with the following assumptions: the system is isotropic such as we analyze only a plane of motion, the shaft presents a linear elastic behavior and the disk masses are lumped in the corresponding nodes. This element has two degrees of freedom in each node, a radial displacement denoted by Rx, and an angle of deflection denoted by βy. The mass, stiffness and gyroscopic matrices for this type of element are given in [3]. To assembly the mass and stiffness matrices, the consistent approach presented by Genta [22] was employed.
Considering the vector of generalized displacements
we can obtain dynamics of undamped system behavior in free vibrations by
Where MT and MR are mass global matrices linked with translational and rotational inertia respectively, GGYR is the global gyroscopic matrix, K is the global stiffness matrix and Ω is the rotational speed.
When the boundary conditions are applied (rigid supports in shaft ends), the vector of generalized displacements becomes
2.2Natural frequencies and critical speedsAs the rotational speed appears explicitly in the equation of motion (Eq. 2), the natural frequencies ω of the system depend on that rotational speed. When this occurs, the behavior of the system is summarized by a plot of the natural frequencies as functions of the rotational speed Ω. Because in many cases the frequencies of the exciting forces depend on the speed, they can be reported on the same plot, obtaining the Campbell diagram.
In the particular case of the disturbance due to mass unbalance, the forcing frequency can be represented on the Campbell diagram by a straight line ω = Ω. The critical speeds occur when the frequency of the forcing function coincides with one of the natural frequencies of the system. It can be seen on the Campell diagram by the intersection of the curves of the natural frequencies with those related with the forcing frequencies. The natural frequencies can be found as eigenvalues of the matrix (see [3]):
where M = MT + MR.
The parameters used to obtain the Campbell diagram and for numerical simulation are presented in Table 1, which correspond with the physical prototype described in [21]. It is important to mention that the eccentricity parameters were experimentally estimated by peak peaking techniques.
System parameters.
Parameter | Value |
Shaft density, ρ | 7850kgm3 |
Shaft diameter, d | 0.020m |
Mass per length unit, m | 2.466kgm |
Disks mass, md1, md2 | 3.9kg |
Length, L1 | 0.25m |
Length, L2 | 0.35m |
Length, L3 | 0.15m |
Young’s modulus, E | 211 GPa |
Moment of inertia of shaft, I | 7.854 × 10-9 m4 |
Diametral moment of inertia per lenght unit of shaft, j | 6.1654 × 10-5kg.m |
Polar moment of inertia per lenght unit of shaft, j¯ | 1.2331 × 10-4kg.m |
Diametral moment of inertia disk,Jd | 0.0057kg.m2 |
Polar moment of inertia of disk, J¯Jd | 0.0113kg.m2 |
Eccentricity of disk 1, a1 | 7.3 × 10-5m |
Eccentricity of disk 2, a2 | 22.2 × 10-5m |
In Figure 2 the Campbell diagram for the two first natural frequencies is depicted. We can observe two branches stemming from a same point, this point is a natural frequency at standstill and the branches can be interpreted as the frequencies of two circular whirling motions, one occurring in the same direction of the rotational speed (forward whirling) and another one in the opposite direction (backward whirling). If the natural frequencies are independent of the spin speed, the Campbell diagram is represented by a straight line parallel to the Ω axis, as in the case in which the gyroscopic effects are neglected and the critical speeds coincide with the natural frequencies at a standstill.
The inclusion of gyroscopic effects in the model tends to increase the value of the forward critical speeds, as we can observe in the Campbell diagram.
The first two critical speed obtained by Figure 2 are ω1 = 33Hz = 1980rpm and ω2 = 109.8Hz = 6592rpm.
When the eigen-problem for matrix (4) is solved, complex numbers are obtained, the complex frequencies. The natural frequencies are given by the imaginary part and they are used to plot the Campbell diagram, while the real part represents the decay rate which can be used to prove system stability [3]. In Figure 3 the real part of the first two complex frequencies is plotted. We can observe negative values for decay rate which implies the system is stable for the range of rotational speed shown (30000rpm).
2.3Mode shapesA modal system analysis was carried out in ANSYS finite element software to obtain the first two natural frequencies and their corresponding mode shapes. The mode shapes represent the deformation pattern induced in each natural frequency. In Figures 4 and 5, mode shapes for the first two natural frequencies are shown.
3Experimental model validationValidation of model is carried out by comparing simulation and experimental results in time and in frequency domains. Experimental data were obtained by proximity sensors in the test rig (shown in Figure 6) operating from 0 to 3000rpm. System response in time is presented in Figure 7, where we can observe that the model reasonable tracks the experimental results.
Then, the Fast Fourier Transforms (FFT) of the data of Figure 7 were obtained in order to get the system frequency response, and these results are presented in Figure 8, where we can see that the resonant peak in simulation corresponds with the critical speed obtained by the Campbell diagram 33Hz, while the resonant peak in experimental data is at 32.9Hz, that is only 0.3 percent of difference.
The amplitude peak observed in Figure 7 between 5 and 6 seconds may be caused by no modeled dynamic effects such as misalignment or bearing characteristics (they are not completely rigid). However, in the frequency response shown in Figure 8 does not present a resonant peak which means it is not a critical speed so that it will be neglected in the design of the active vibration control.
4System with active suspensionThe system with the active suspension is presented in Figure 9.
When the active suspension is incorporated to the system, one more degree of freedom must be added because the radial displacement at the right end is not constrained anymore, so the new generalized displacements are
The rotor system dynamics is given by
where Mcl, Gcl, Dcl and Kcl are the mass, gyroscopic, damping and stiffness matrices, respectively, for the generalized displacements given by Eq. 5, f(t) is the force control input given by the linear actuators, ϖ1 and ϖ2 are the disturbance forces due to mass unbalance of disks, b, e1 and e2 are vectors indicating the degrees of freedom where the forces are acting.
The disturbance forces in transient state are
with dφdt=ΩIt is important to remark that critical speeds change when the active suspension is incorporated to rotor system. Solving the eigen-problem with the new boundary conditions, the first critical speed is at 39.99Hz = 2400rpm.
5Active control scheme5.1Controllability and observabilityBy defining the state vector z = [u u]T, dynamics system given by Eq. 6 can be described in state space form
with
From the state space representation, controllability and observability matrices can be computed as
which are of full rank, proving that the system of Eq. 9 is completely controllable from the control force f(t) and completely observable from the output Rx2 (radial displacement of the first disk).
5.2Asymptotic observerIn order to design a LQR controller based on a full state feedback, first, we design a full order asymptotic (Luenberger) observer (see e.g. [23]) to estimate the state vector considering that only the system output Rx2 and the control input f(t) are available. The asymptotic observer is expressed by the dynamics
where zˆ is the estimated state vector and L is a feedback gain matrix.
By defining the estimation error as
and properly selecting the matrix L, the asymptotic convergence of the estimation error to zero can be guaranteed limt→∞z˜;=0, the estimated state zˆ converges to the actual state z In order to satisfy this condition, the matrix L is computed via pole-placement technique, such that the observer dynamics described by Eq. 13 be sufficiently faster than the system dynamics given by Eq. 9.
5.3LQR controlUsing the estimated state vector given by the observer, a Linear Quadratic Regulator (LQR) based on full state feedback to attenuate the vibration amplitudes is proposed as follows
where G is a gain matrix that minimizes a conventional LQR performance index with with Q > 0 and r > 0 (see, e.g., [23]).
The linear electromechanical actuators are DC motors with geared reducer and a gear-screw type mechanism to transform the rotational in translational motion. Therefore, in order to take into account the actuators dynamics, the control force obtained by Eq. 15 is transformed and considered as a reference current for a PID current control for the linear actuator, whose output is the applied force to the system by the active bearing. A schematic diagram of the closed loop system is depicted in Figure 10.
The actuator dynamics in state space is given by
where w1 is the armature current, w2 the angular position, w3 the angular speed, Ra the armature resistance, La the armature inductance, kem the back emf constant, u the input voltage, kT the torque constant, Jo the rotor inertia, bo the viscous damping coefficient, T the output torque, T2 the torque after the geared reducer, n the reduction constant, F the linear force provided by the actuator and r the radius in the gear-screw mechanism. The motor parameters are presented in Table 2. Some parameters were obtained from the manufacturer and other ones were estimated by experimental tests.
The PID current control action is obtained by
with yˆ2 the reference current, kp, kd and ki the proportional, derivative and integral gains, respectively.
6Numerical and experimental resultsAs we mentioned above, many rotating machines must operate to high rotating speeds which yields a supercritical operation, that is, during the running up the machine has to cross one or more critical speeds causing large vibration amplitudes [2,3]. The proposed controller was tuned to the critical case, when rotor system passes through the first resonant peak.
In Figures11 and 12, the numerical results of the rotor system with the active suspension in open and in closed loop are presented. System response is presented in time and in frequency domain from the standstill to operation speed (3000rpm). A substantial reduction in the vibration amplitudes at the resonant peaks can be observed. In [24, 25], the response of the closed loop rotor system at constant speed is presented showing that the proposed control scheme can work in this condition too.
In Figures11 and 12 we can appreciate a magnification in closed loop system response in relation to open loop system in some low operation speeds, particularly at disk 1 location (Rx2). This may be due to the fact that the controller was tuned for attenuation in the critical case (resonant peak corresponding to the first critical speed). Moreover, in a specific application, if those amplitudes were excessive, the controller would be triggered just before the system reaches the resonant condition to operate only in a time interval.
In order to validate the control scheme proposed, some experiments were carried out. In Figure 13 the test rig with the active suspension used to carry out the experiments is presented.
The results are shown in Figures 14 and 15, in time and in frequency domain, respectively. The system operates from standstill to 3000rpm. We can observe important reductions in the vibration amplitudes in both disks during the resonant peak.
7ConclusionsIn this work an active control scheme to reduce the unbalance response in a rotor system was presented. To do this a reduced order model for the rotor system was obtained by finite element techniques. This model includes the gyroscopic effects which allowed us to obtain the Campbell diagram to represent the dynamic behavior of the rotor system and it was experimentally validated. The inclusion of gyroscopic effects permits us to differentiate between the concepts of natural frequency and critical speed as well as to predict more precisely the resonance condition. Also, the decay rate plot was obtained from the finite element model, showing the system is stable in a wide range of rotational speed. The model was experimentally validated and used to design an asymptotic observer to estimate the full state vector in order to propose a linear quadratic regulator with full state feedback to attenuate the vibration amplitudes when the rotor pass through its first critical speed. This scheme was proved by numerical simulations. A novel active suspension using electromechanical actuators was designed and constructed in order to validate experimentally the closed loop system performance. Experimental results demonstrate reductions in the unbalance response about 66 percent in the first disk and 44 percent in the second one. Furthermore, the applied control efforts acquire reasonable values even when the system goes through its critical speed. It is important to mention that the power consumption of the actuator is around 2 W while the rotor system is driven by a motor of 1HP.