Active Vibration Control of Rotating Machines with Active Machine Foot Mounts on Soft Foundations Based on a 3D-Model ()
1. Introduction
In praxis, large rotating machines are often fixed directly on soft foundations. e.g. elastic steel frame foundations, which sometimes lead to problems regarding vibrations, caused by resonances and instability [1] - [7]. In many technical applications, active vibration control is used to solve vibration problems [8] - [20]. Ushijima and Kumakawa analyzed in [12] Piezo-Actuators, used as active engine mounts for vibration control for high vibrations. Ulbrich investigated in [13] different actuator concepts especially for rotating machinery and compared them regarding their effectiveness. Also Sun (et al.) investigated in [14] vibration control of car active suspensions, but with electrohydraulic actuators. Sohn (et al.) presented in [15] experimental investigations to evaluate an electromagnetic actuator. Zhang (et al.) presented in [16] a review of the current status and developmental progress of active engine mounts, where a theoretical model, a finite-element model, and experimental modelling of electromagnetic active engine mounts are shown. In this paper now, the active vibration system includes active machine foot mounts—actuators, which are put between the machine feet and the foundation, acting only in vertical direction—vibration sensors, which are mounted at each machine foot, detecting the vertical vibrations, and a separate controller for each actuator (Figure 1). This concept was basically investigated in [21], but only for induction motors and it was only based on a plane 2D-model. In [22] a general 3D model for rotating machines mounted on actuators was presented, but only for rigid foundations. In this paper now, the 3D model is enhanced for soft foundations, e.g. steel frame foundations. An additional large enhancement of this paper is, compared to [22], that the excitations are here not only forces and moments with constant amplitudes and arbitrary frequencies, acting separately at the x, y, z-axes, but a rotating unbalance force and a rotating unbalance moment, described by complex rotating vectors. Furthermore, the damping coefficients are here dependent on the whirling angular frequency, whereas in [22] the damping coefficients have been considered constant as a simplification.
2. Vibration Model
The model, which is used here, is a simplified 3D model (Figure 2). The mass and the moments of inertia of the machine are concentrated at the center of gravityS. The rotor and stator have hereby the same center of gravity. The angular frequency Wrepresents the rotational angular frequency of the rotor. The stiffness of the rotor, of the bearings, of the bearing support and of the stator is supposed to be infinitely high. Following definitions are made:
Figure 1. Rotating machine (induction motor) mounted on a soft foundation with actuators.
Figure 2. Simplified 3D-model of rotating machine on a soft foundation with actuators.
· The mass
of the machine consists of the mass
of the stator and of the mass
of the rotor and is included in the mass matrix, as a single mass.
(1)
· The moments of inertia
and
, consisting of
and
(the moments of inertia of the stator at the y-axis and z-axis) and of
(the moments of inertia of the rotor at the y-axis and z-axis), are also included in the mass matrix.
;
(2)
· However, the moment of inertia of the stator at the x-axis
and the moment of inertia of the rotor at the x-axis
have to be used separately, because
has to be included in the mass matrix, and
in the gyroscopic matrix.
· Additionally
and
and
have to be handled separately, because they have to be considered in the excitation vectors, which are described by
and
.
At the centre of gravity S the two excitation vectors
and
are positioned. The vector
represents the rotating force cause by a static unbalance of the rotor, due to a mass eccentricity
with phase shift of
. The vector
represents the rotating moment caused by a moment unbalance of the rotor, due to a tilt—angle
with phase shift of
—of the rotor mass on the shaft. The center of the bearing housings are described by point
(on the drive side) and by point
(on the non-drive side). The boreholes of the machine feet are defined by the points
,
,
and
. Index D stands for drive side, index N for non-drive side and index L for left side and index R for right side:
;
(3)
The mass
,
,
and
present the mass of each amateur of the actuators and
,
,
and
the mass of each stator of the actuators. In the stiffness and damping matrix of each actuator only translational stiffness and damping is considered, referring to [22]:
(4)
The coefficients
,
and
are the stiffness coefficients of the actuators and the coefficients
,
and
are the damping coefficient of the actuators. For each foundation point, also only translational stiffness and damping of the foundation is considered in the stiffness and damping matrix:
(5)
The coefficients
,
and
are the stiffness coefficients of the foundation and the coefficients
,
and
are the damping coefficient of the foundation. Nine global coordinate systems are used here, one at the centre of gravity S and one for each machine foot point (
,
,
,
) and one for each foundation point (
,
,
,
). By using a vibration sensor for each machine foot, the vertical machine foot displacements or velocities or accelerations can be detected and then lead back to separate controllers (
,
,
and
). These controllers create then suitable signals for producing vertical actuator forces (
,
,
and
). The controller structure of each controller may be different. The structure of each controller is described in the Laplace domain, by a transfer function, with the Laplace variable s:
Controller
:
(6)
The constants
and
are hereby the constants of the polynomial functions. To describe different feedback strategies the index
is used, referring to [22]:
(7)
The mechanical damping coefficients of the actuators
and of the foundation
can be derived by the corresponding mechanical loss factor
for the actuators and
for the foundation, by the corresponding stiffness
and
and by the whirling frequency
, referring to [4] [7] and [21]:
(8)
For analysis of the forced vibration due to unbalance, the whirling angular frequency is equal to the rotor angular frequency
. For analysis of natural vibrations with marginal decay, the whirling angular frequency is defined as the correspondent natural angular frequency, as a simplification, referring to [4].
3. Mathematical Description
3.1. Description in the Time Domain
A linearization regarding the machine feet displacements is possible, because only small displacements occur due to the excitations. Therefore, kinematic constraints can be used to express the movement of the motor feet, as well as the displacements of the bearing housing points
and
, by the movement of the center point S, referring to [22]. The vibration system can be described by following differential equation:
(9)
The vectors
and
are the vectors of excitation and the vector
is the vector of the actuator forces. The vector
contains the coordinates for displacements and rotations of the machine centre point S and for displacements of the foundation points.
(10)
The matrix
represents the mass matrix, the matrix
the damping matrix and the matrix
the stiffness matrix:
(11)
The coefficients of the matrices are presented in the Appendix. The gyroscopic matrix
can be written by:
(12)
The vector
, represents the zero-vector with
, the vector
the zero-vector with
and the matrix
the zero-matrix with
. The unbalance force vector
, caused by static unbalance
, can be described by:
(13)
and the unbalance moment vector
, caused by the moment unbalance
, by:
(14)
The actuator force vector
is split into the actuator force vector for each machine foot:
(15)
with the actuator force vectors:
(16)
(17)
(18)
(19)
The vector
represents the zero-vector with
. The vectors
,
,
and
are the actuator force transmission vectors. Now, a state space formulation is used (Figure 3), based on [17] [18] [19] [20].
For avoiding mix-up with the stiffness matrix
and damping matrix
, index “st” is used for the matrices of the state space. According to Figure 3, the state space vector
and the output vector
can be written as follows, referring to [22]:
(20)
Figure 3. State space model for active vibration control with negative feedback of the output vector, considering excitation due to static unbalance and moment unbalance.
The system matrix
, the input matrix
, the output matrix
, and the straight-way matrix
can be written as:
(21)
with the zero-matrix
and the unit-matrix
. Therefore, the state space equations can be described by:
(22)
(23)
By using a controller matrix
(Figure 3)—which will be defined later— the actuator force vector
can be described by:
(24)
3.2. Description in the Fourier-Domain for Excitation by Unbalance
For deriving the frequency response matrix of the system, Equations (22), (23) and (24) are transferred in the Fourier-domain:
(25)
(26)
(27)
With these equations follows the output vector:
(28)
with the unit-matrices
and
. With the definition
(29)
it is highlighted, that the output vector
dependents on the different feedback strategies(7), represented by the controller matrix
. The controller matrix
can now be derived, based on [22].
(30)
with the matrix
, also described in the Fourier-domain with
, which are the frequency response functions of the controllers, based on the controller transfer functions (6) with
:
(31)
With:
and
and with the zero-matrix
. Now, the output vector
can be calculated by:
(32)
and the frequency response matrix
can be written by:
(33)
The frequency response vector for each single kind of excitation can now be derived, with following Fouier-transformations, with the Dirac-delta function
:
(34)
(35)
Further on, the index
is used for both kinds of unbalance, representing the excitations:
(36)
The output vector
for each kind of unbalance can be calculated by:
(37)
Using the sifting property of the Dirac delta function
, the Fourier-transformed output vector for each kind of unbalance follows:
(38)
Now, the inverse transformation of this Fourier-transformed output vector
back into the time-domain can be done:
(39)
It is useful to relate the amplitude output vector on the respective unbalance:
(40)
At this point it is important to mention, that beside the gyroscopic matrix
also the damping matrix
depends on the rotary angular frequency W. The reason is that for forced vibration due to the unbalance, the whirling angular frequency
for the mechanical damping coefficients (8), is equal to the rotary angular frequency W (
).Therefore, the matrices
and
are also functions of the rotor angular frequency W. The frequency response matrix
can now be formulated by:
(41)
With this matrix, the frequency response vector
for both kinds of unbalance can be described:
(42)
With this frequency response vector
and the kinematic constraints in [22], the response functions for vibration velocities of the bearing housings, and of the foundation points can be derived, as well as the response functions for the actuator forces, based on [22]. The equations are shown in the Appendix.
3.3. Stability Analysis
Due to the use of active vibration control, a stability analysis is very important, therefore the poles of the system have to be derived. Based on [22], the poles can be calculated by solving the following equation:
(43)
However, this direct procedure is only possible if the matrices
and
are independent of the whirling angular frequency
, which is here not the case, because the damping coefficients depend here on
. The whirling angular frequency
corresponds here—when calculating the pols—to the natural angular frequency and therefore to the imaginary part of the complex poles. Therefore, the use of the damping coefficient lead here to a causality problem. Therefore, a worst case procedure is here derived to investigate the stability of the system. The basic requirements are:
· The mechanical damping of the actuators and the foundation is low, which should be usually the case, so that the mechanical loss factors for the actuators and the foundation fulfill following conditions:
and
(44)
· The modal mass of the foundation is so low, that only the first six natural vibration modes have to be taken into account.
· The first six natural vibration modes have conjugate complex poles.
Therefore, following worst case procedure is derived to calculate the pols
, which is shown in Figure 4. In the first step, the pols are calculated without damping (
) and with open control loops (
).
Then, the natural angular frequency of the 6th conjugate complex pole pair
—which is set to be equal to the whirling angular frequency
—is considered, for calculating the damping coefficients
and
. Afterwards, the pols are calculated again, considering damping (
) and the closed control loop operation (
). The corresponding conjugate complex pol pair
is again taken into account for deriving the new damping coefficients. Finally the pols are calculated with the modified damping matrix
. If the change of the whirling angular frequency is too large (e.g. more than 5%) this approach can be repeated in loops. This procedure presents a worst case scenario regarding instability, because the natural angular frequency of the highest considerable mode (here mode 6) is used for calculating the damping coefficients, which lead to the
lowest damping coefficients.
Figure 4. Flow diagram for calculating the pols of the vibration system.
3.4. Mathematical Description for Special Controllers
In this section, special controllers—standard controllers (P-, I-, PI-, PD-(ideal), PID-(ideal) controllers)—are analyzed. Therefore, the control parameter can now be directly implemented into the mass matrix and/or the damping matrix and/or the stiffness matrix, which can be seen in Figure 5, referring to [22].
Therefore, the differential equation can be described now by:
(45)
By implementing the control parameters into the mass matrix,
becomes
. When they are implemented into the damping matrix,
becomes
, and if they are implemented into the stiffness matrix,
becomes
. The structure of these matrices is shown in (49). Following definitions are used, referring to [22]:
For the mass matrix
:
and
(46)
For the damping matrix
:
and
(47)
For the stiffness matrix
:
and
(48)
The matrices with the integrated controller parameters are formulated by:
Figure 5. Matrices, depending on standard controller structures and different feedback strategies, referring to [22].
(49)
The controller parameters
depend on the controller structure and the chosen feedback strategy (Figure 5). The combination in Figure 5 marked with “¾” cannot be deduced, if the differential equation system (45) shall be used. If the control parameters are identical, the cells with identical color lead to the same differential equation. The inhomogeneous differential equation (45) can now be solved, based on the common approach.
4. Numerical Example
Now, a numerical example is presented, where the vibration stability, as well as the forced vibrations of a rotating machine are analyzed.
4.1. Boundary Conditions
The data of the rotating machine are shown in Table 1. As it can be seen, the rotating machine has a symmetrical machine design (bL = bR; aD = aN).
The foundation is a simplified steel frame foundation, which consists of two I-beams (green), which are stiffened by additional welded steel sheets (brown) in the area of the machine feet (cyan), and fixed (red) to ground (Figure 6). No
Table 1. Data of the rotating machine.
Figure 6. Simplified model of the steel frame foundation.
stiffening effect in the area of the machine feet is considered, so that a low stiffness of the steel frame foundation is considered here, as a worst case.
The substitute foundation stiffness, which is necessary to use the model in Table 2, is derived by a finite element analysis of the foundation, but still presenting a simplification.
In this example only three actuators are used, two on drive-side, left and right (DL, DR) and only one actuator on non-drive side left (NL) (Table 3). The fourth actuator on the non-drive side right (NR) is exchanged by a stiff element.
The data of the stiff element is shown in Table 4.
For the control system, I-controllers with feedbacks of the vertical machine feet accelerations are chosen, so that following controller transfer function is used.
(50)
Table 3. Data of the actutaors on position DL, DR and NL.
Table 4. Data of the stiff element on position NR.
Table 5. Data of the controller parameters.
The data of the control parameters are shown in Table 5. All three controllers—for the three actuators—have identical parameters. The control parameter for the fourth controller on NR has to be set to zero, because here a stiff element is used instead of an actuator.
Afterwards three different cases are analyzed:
· Case 1: The rotating machine is directly mounted on the steel frame foundation.
· Case 2: Three actuators at DL, DR, NL and a stiff element at NR are placed between machine feet and foundation. The actuators are only operating passively (open control loops).
· Case 3: Same setting as case 2, but the actuators are now operating actively (closed control loops).
4.2. Stability Analysis
To analyze the stability of the system, the poles are calculated for the three different cases, based on the procedure, described in Figure 4. It can be clearly shown, that for all three cases vibration stability exists, because all characteristic poles have no positive real parts (Figure 7, Figure 8 and Figure 9).
When comparing the figures, it is obvious that the damping of the poles can be strongly increased by the active vibration control system (case 3), compared to case 1 and case 2.
4.3. Frequency Response Analysis
Now, the amplitudes of the frequency response functions of bearing housing vibration velocities and of foundation vibration velocities are computed, as well as the frequency response functions of the actuator forces, all related to the respective unbalance. The amplitudes of the frequency response functions are calculated in [dB] with the reference gauge in Si-Unit, in a frequency range from 1 Hz to 250 Hz. An additional index d is now introduced for case 1, so that now following definitions are used:
(51)
Related bearing housing vibration velocities:
Figure 7. Poles for case 1, calculated at rotational frequency
and
.
Figure 8. Poles for case 2, calculated at rotational frequency
and
.
Figure 9. Poles for case 3, calculated at rotational frequency
and
.
(52)
with:
Related foundation vibration velocities:
(53)
with:
Related actuator forces:
(54)
with:
4.3.1. Excitation by Static Unbalance
The amplitude response functions for the bearing housing vibration velocities at the drive side and at the non-drive side, related to static unbalance are shown in Figure 10.
The amplitude response functions for the foundation vibration velocities at the drive side, related to static unbalance are shown in Figure 11.
The amplitude response functions for the foundation vibration velocities at the non-drive side, related to static unbalance are shown in Figure 12.
The frequency response functions for the actuator forces, related to static
Figure 10. Amplitude response functions for the bearing housing vibration velocities at the drive side and at the non-drive side, related to static unbalance
for case 1, case 2 and case 3.
Figure 11. Amplitude response functions for the foundation vibration velocities at the drive side (left and right), related to static unbalance
for case 1, case 2 and case 3.
Figure 12. Amplitude response functions for the foundation vibration velocities at the non-drive side (left and right), related to static unbalance
for case 1, case 2 and case 3.
Figure 13. Frequency response functions for the actuator forces, related to static unbalance
(case 3) at drive side left (DL), drive side right (DR), non-drive side left (NL) and non-drive side right (NR).
unbalance are shown in Figure 13.
4.3.2. Excitation by Moment Unbalance
The amplitude response functions for the bearing housing vibration velocities at the drive side and at the non-drive side, related to moment unbalance are shown in Figure 14.
The amplitude response functions for the foundation vibration velocities at the drive side, related to moment unbalance are shown in Figure 15.
The amplitude response functions for the foundation vibration velocities at the non-drive side, related to moment unbalance are shown in Figure 16.
The frequency response functions for the actuator forces, related to moment unbalance are shown in Figure 17.
4.3.3. Discussion of the Frequency Response Analysis
The amplitude response functions for the bearing housing vibration velocities
Figure 14. Amplitude response functions for the bearing housing vibration velocities on the drive side and on the non-drive side, related to moment unbalance
for case 1, case 2 and case 3.
Figure 15. Amplitude response functions for the foundation vibration velocities at the drive side (left and right), related to moment unbalance
for case 1, case 2 and case 3.
Figure 16. Amplitude response functions for the foundation vibration velocities at the non-drive side (left and right), related to moment unbalance
for case 1, case 2 and case 3.
Figure 17. Frequency response functions for the actuator forces, related to moment unbalance
(case 3) at drive side left (DL), drive side right (DR), non-drive side left (NL) and non-drive side right (NR).
for excitation by a static unbalance (Figure 10) show, that in the analyzed rotational frequency range two resonance frequencies occur for case 1—where the machine is directly mounted on the steel frame foundation—at about
and at
. Figure 10 shows additionally, that no axial vibrations occur for case 1, if only a static unbalance is considered. For a moment unbalance, three resonance frequencies occur for case 1, at about
,
and
(Figure 14). The natural mode, where the machine mass is mainly rotating at the x-axis—with a natural angular frequency of about 1088 rad/s (Figure 7), which means a natural frequency of about 173, 2 Hz—leads not to a resonance regarding bearing housing vibrations, because this mode shape is hardly excited by these two kinds of unbalance. Because of the symmetrical machine design and the identical foundation stiffness and damping coefficient under each machine feet, the frequency response functions for the bearing housing vibration velocities are identical for the drive side and the non-drive side for case 1. If now the machine is mounted on a stiff element and on three actuators, operating with open control loops (case 2), the resonance frequencies for the bearing housing vibrations occur at lower frequencies and six different resonance frequencies are obvious, when considering both kind of unbalance (Figure 10 and Figure 14). The amplitudes in the resonance frequencies are clearly lower for case 2 compared to amplitudes in the resonance frequencies for case 1, because the amplitudes of excitation force and of the excitation moment are proportional to
. When now the control loops are closed— which means case 3—most of the resonances regarding the related bearing housing vibrations (Figure 10 and Figure 14) get well damped—as well as most of the related foundation vibrations (Figure 11, Figure 12, Figure 15, Figure 16), when comparing the red curves with the blue curves. Therefore, an operation in the rotational frequency range from 0 Hz to 80 Hz seems now possible for case 3, without off-limits areas regarding rotational frequency. If the machine would be mounted directly on the foundation (case 1), this would not be possible, because a sharp resonance at 38.1 Hz occurs here due to static unbalance, regarding the related bearing housing vibrations (Figure 10) and regarding the related foundation vibrations (Figure 11 and Figure 12). Figure 13 and Figure 17 show the frequency response functions regarding the active actuator forces. It can be seen, that with increasing rotational frequency also the amplitudes of the actuator force frequency response functions increase, because the unbalance force and unbalance moment also increase with the rotational frequency. Of course, the actuator forces are also strongly influenced by the vibration mode shapes. No active actuator force occurs on position NR, because a stiff element is positioned there instead of an actuator. If the stiff element would be exchanged by a fourth actuator, so that four identical actuators would be used—case 2 gets case 2* and case 3 gets case 3* (with four identical control parameters)—a symmetrical system would exist. With such a symmetrical system, a vibration mode exists, where the machine makes a pure rotating at the z-axis. This mode cannot be influenced by the actuators anymore, because they are only acting in vertical direction. This is shown in Figure 18 regarding the amplitude response functions for the bearing housing vibration velocities in horizontal direction. The resonance frequency at about 18.6 Hz cannot be damped by the control system. Only a marginal difference between the blue curve and the red curve for horizontal direction
Figure 18. Amplitude response functions for the bearing housing vibration velocities on the drive side and on the non-drive side, related to moment unbalance
for case 1, case 2* and case 3* (*: Four identical actuators are used).
is obvious in Figure 18, which is caused by the gyroscopic effect, leading here to a marginal coupling of the vibration modes, so that no absolute pure vibration mode at the z-axis occurs. Because of the symmetrical system, the frequency response functions for the bearing housing vibration velocities are now identical for drive side and non-drive side (Figure 18).
Therefore, it could be shown, that the method “vibration mode coupling by asymmetry”, which was developed and mathematically described in [22] for a rigid foundation, is also transferable for a soft foundation.
5. Conclusion
The paper presents a simplified 3D-model for active vibration control of rotating machines with active machine foot mounts on soft foundations, considering static and moment unbalance. After the model was mathematical described in the time domain, it was transferred into the Fourier domain, where the frequencies response functions regarding bearing housing vibrations, foundation vibrations and actuator forces have been derived. Afterwards, the mathematical coherences have been described in the Laplace domain and a worst case procedure was derived to analyze the vibration stability. For special controller structures in combination with certain feedback strategies, a calculation method was shown, where the controller parameters can directly be implemented into the stiffness matrix, damping matrix and mass matrix. Additionally a numerical example was presented, where the vibration stability and the frequency response functions have been analyzed. It could be shown, that with the active vibration control system all vibration modes can be damped well, so that an operation in the rotational frequency range from 0 Hz to 80 Hz is possible without off-limits areas.
Finally it could be demonstrated, that the method “vibration mode coupling by asymmetry”, which was developed and mathematically described in [22] for a rigid foundation, is also transferable for soft foundations. The numerical example was also used, to validate the presented model with a finite element analysis by comparing the results. However also a validation by measurement is planned. Therefore, a test bed for an 11 kW induction motor, mounted on a soft steel frame foundation with actuators, has been built and first experimental investigations have been deduced [23], showing good congruence with the presented 3D-model.
Appendix:
A1. Coefficients of the Matrices
Coefficients of the mass matrix M:
(55)
(56)
(57)
(58)
(59)
(60)
(61)
(62)
(63)
(64)
(65)
(66)
(67)
(68)
(69)
(70)
(71)
The other coefficients are zero.
Coefficients of the stiffness matrix C:
(72)
(73)
(74)
(75)
(76)
(77)
(78)
(79)
(80)
(81)
(82)
(83)
(84)
(85)
(86)
(87)
(88)
(89)
(90)
(91)
(92)
(93)
(94)
(95)
(96)
(97)
(98)
(99)
(100)
(101)
(102)
(103)
(104)
(105)
(106)
(107)
(108)
(109)
(110)
(111)
(112)
(113)
(114)
(115)
(116)
(117)
(118)
(119)
(120)
(121)
(122)
(123)
(124)
(125)
(126)
(127)
(128)
(129)
(130)
(131)
(132)
(133)
(134)
The other coefficients are zero.
Coefficients of the damping matrix D:
The coefficients of the damping matrix can be easily derived by just replacing “c” in the formulas (72)-(134) by “d”.
A2. Response Functions for the Bearing Housing Vibration Velocities
· In vertical direction:
Drive side:
(135)
Non-drive side:
(136)
· In horizontal direction:
Drive side:
(137)
Non-drive side:
(138)
· In axial direction:
Drive side:
(139)
Non-drive side:
(140)
The element
is the nth element of the frequency response vector
.
A3. Response Functions for the Foundation Points Velocities
· In vertical direction:
Drive side, left:
(141)
Drive side right:
(142)
Non-drive side, left:
(143)
Non-drive side, right:
(144)
· In horizontal direction:
Drive side, left:
(145)
Drive side right:
(146)
Non-drive side, left:
(147)
Non-drive side, right:
(148)
· In axial direction:
Drive side, left:
(149)
Drive side right:
(150)
Non-drive side, left:
(151)
Non-drive side, right:
(152)
A4. Response Functions for the Actuator Forces
Feedback of the motor feet displacements
:
· For drive side, left:
(153)
· For drive side, right:
(154)
· For non-drive side, left:
(155)
· For non-drive side, right:
(156)
For feedback of the motor feet velocities, index
gets v (
) and the index number n of the elements
in (153)-(156) has to be changed to
. For feedback of motor feet accelerations index
gets a (
) and the index number n of the elements
in (153)-(156) has to be changed to
. For open loop operation index
gets 0 (
), and the response functions of the actuator forces are zero.