Modelling of transient and steady-state modes of a vertical rotor with an automatic balancing device
Guntis Strautmanis^{1} , Gennadiy Filimonikhin^{2} , Mareks Mezitis^{3} , Alexander Gorbenko^{4} , Valentina Strautmane^{5} , Zura Sansyzbajeva^{6}
^{1}Riga Technical University, Riga, Latvia
^{2}Central Ukrainian National Technical University, Kropyvnytskyi, Ukraine
^{3}Transport Academy, Riga, Latvia
^{4}Kerch State Maritime Technological University, Kerch, Russia
^{5}Railway Transport Institute, Daugavpils, Latvia
^{6}L. N. Gumilev Eurasian National University, Astana, Kazakhstan
^{1}Corresponding author
Journal of Vibroengineering, Vol. 23, Issue 3, 2021, p. 759-769.
https://doi.org/10.21595/jve.2021.21804
Received 15 November 2020; received in revised form 19 December 2020; accepted 5 January 2021; published 31 January 2021
The authors of the study work out the differential equations of motion of a vertical rotor model on an elastic-dissipative suspension, balanced by a ball-type automatic balancing device. Often, the cross-section of the cavity of the body of the automatic balancing device is rectangular and during rolling the balls have two points of contact, in one of which the balls slide along the surface of the cavity. To prevent the balls from sliding, the inner surface of the cavity of the automatic balancing device is made in the shape of a torus, which provides one point of contact. The forces of gravity and the forces of resistance to the movement of the correcting weights are taken into account, and the model is drawn up for both viscous and dry friction forces inside the body of the automatic balancing device. The obtained mathematical model of the rotor makes it possible to study the transient and steady-state modes of motion of the rotor system.
- The authors obtained a system of differential equations for a rotor system with a ball-type ABD, in which the inner space of the body is made in the shape of a torus, which provides a minimum resistance force to the motion of the ball and the best conditions for the location of balls in the auto-balancing mode.
- The proposed system of differential equations makes it possible to proceed to the study of simpler rotor systems with a vertical axis of rotation, equipped with an ABD, in which the balls move in both longitudinal and transverse directions.
- The system of differential equations makes it possible to investigate both steady-state modes of rotor motion and transient modes.
Keywords: rotor, automatic balancing device, equations of motion, elastic element, dissipation coefficient, transient mode.
1. Introduction
One of the design stages of a rotor system with a ball-type automatic balancing device (ABD) is the development and substantiation of the design model of the rotor and its mathematical model. The advantage of this type of the ABD is that the balls themselves are set opposite to the rotor unbalance in the superresonance zone of the rotor motion frequency [1-3].
Most of the studies and publications are concerned with the consideration of the steady-state mode of rotor systems with ABD. At the same time, transient modes, although relatively short-lived, are accompanied by increased oscillation amplitudes. The role of the ABD in transient modes of the rotor system has not been studied to the full. This is due, inter alia, to the complexity of mathematical modelling not only of steady-state [5], but especially of transient processes [6].
In [15], it is indicated that in transient modes at insufficient engine power, the Sommerfeld effect may occur, in this case the ABD aggravates increased vibrations. In the scientific literature there are publications [21], where the issue of reducing resonant vibrations is considered by varying the parameters of the rotor system when installing the ABD. However, the influence of the ABD parameters on the possibility of reducing resonance amplitudes is insufficiently disclosed in the scientific literature. The elaboration of a mathematical model of a rotor with an ABD for a transient mode will allow further deeper and comprehensive study of this issue.
In scientific publications axisymmetric rotors on two isotropic [10] or anisotropic elastic supports [11, 13] are mainly considered. At the same time, a large number of rotor systems have a cantilever layout. In this case, in the motion of such rotor systems, both translational and angular oscillations are realized, which complicates both the derivation of the differential equations of motion and the analysis of the calculation results. At the same time, the elaboration of new and improvement of existing rotor systems requires the solution of such complex mathematical models. As a rule, such mathematical models imply their numerical integration, which requires the use of special software, in particular, the SPRING software described in [12] and allowing implementation of various research tasks.
In [16-19], it is mentioned that weights in the form of balls in the acceleration mode of the rotor and in the stationary mode roll along the ABD treadmill, at the same time it is indicated that the movement of weights in the ABD body occurs in the same plane (along the treadmill). When the ball rolls on the treadmill, the contact with the ABD body is provided at one point, which reduces the resistance force. In real conditions, during the acceleration of the rotor, the ball, for various reasons, inside the ABD body begins to move both along the treadmill and in the transverse direction, which can lead to the emergence of side steady modes that do not allow the ball to accelerate to the operating speed of the rotor. This is confirmed by mathematical modelling of the conditions for acceleration of the weight, considered in [6]. Thus, it becomes necessary to take into account the transverse motion of the ball in the mathematical model.
Cyclically alternating weights in rotor systems also require ensuring their fatigue strength. To check the strength of the rotor and ABD parts, one can, in particular, use the approach mentioned in [14].
When creating a mathematical model of a rotor system with an ABD, the laws of theoretical mechanics [4] and the theory of rotors with automatic balancing devices [7-9] were applied.
The aim of this work is to obtain a mathematical model of a cantilever rotor system with a ball-type ABD in transient and steady-state modes taking into account the longitudinal and transverse movements of the balls inside the ABD body.
2. Mechanical model of the rotor system
The rotor system consists of a platform, a rotor, and a ball-type ABD mounted on the rotor. The ABD contains $n$ identical balls of radius $r$. The balls are inside a toroidal cavity (torus). The radius of the generating circumference of the torus is ${R}_{1}$, and of the guide in the circumferential direction – $R$.
The platform is fixed on elastic-viscous supports. The rotor oscillates along with the platform and rotates relative to the platform. The masses of the rotor, platform and one ball are, respectively, equal to ${m}_{r}$, ${m}_{p}$, $m$.
The design model of the rotor system is shown in Fig. 1. The positions of the rotor system are determined by the following generalized coordinates:
– linear coordinates $x$, $y$, $z$ of the origin of the moving coordinate system ${0}_{1}{x}_{1}{y}_{1}{z}_{1}$ relative to the fixed system $0xyz$;
– angular coordinates $\alpha $, $\beta $, $\theta $ (angles of inclination of the rotor axis in plane ${0}_{1}{x}_{1}{z}_{1}$ and ${0}_{1}{y}_{1}{z}_{1}$, and the angle of twist around ${z}_{1}$ axis, respectively);
– angular coordinates ${\psi}_{j}$, ${\gamma}_{j}$, $/j=\stackrel{-}{1,n}/$, determining the position of the center of $j$ ball in the ABD.
When deriving the equations of motion of the system, the assumptions are made that:
– linear and angular coordinates of the platform are sufficiently small $\left|x\right|$, $\left|y\right|$, $\left|z\right|$, $\left|\alpha \right|$, $\left|\beta \right|$, $\left|\theta \right|\ll 1$ and their time derivatives are small, too;
– balls roll on the surface of the torus without sliding.
The movement of compensating weights (balls) inside the torus is considered in two sections – in the circumferential and transverse directions.
The study of the steady-state mode of the rotor system is carried out at a constant angular velocity of the rotor rotation, i.e. $\omega =$ const. The angle of the rotor rotation in this case is determined by the expression $\phi =\omega t$. In the study of transient modes, the angular velocity of the rotor rotation changes and is a given function of time, i.e. $\omega =f\left(t\right)$.
Fig. 1. Design model of the rotor system; a) general view; b) torus cross section
a)
b)
The tensors of the moments of inertia, characterizing the correspondingly distributed rotating $i$masses of the rotor and the ABD body, as well as ${m}_{p}$ of the platform, have the form:
and ${I}_{xy}={I}_{yx}$; ${I}_{xz}={I}_{zx}$; ${I}_{yz}={I}_{zy}.$
The radius vectors determining the position of the centers of the rotor masses (point ${C}_{1}$) and platform (point ${C}_{2}$, not shown in Fig. 1), as well as $j$ ball in the ABD (points ${C}_{j}$) in the moving coordinate system ${0}_{1}{x}_{1}{y}_{1}{z}_{1}$, have the following form:
where $\phi $ – angle of the rotor rotation around its axis; $h$ – distance from the center of the rotor masses to the suspension plane (not shown in Fig. 1); $e$ – value of the static eccentricity of the rotor.
The differential equations of motion of the rotor system are obtained using the Lagrange equations of the 2nd kind, which have the form:
where $T$ – kinetic energy of the rotor system; $\u041f$ – potential energy of the system; ${R}_{d}$ – dissipative function in generalized coordinates; ${Q}_{i}$ – generalized strength.
The kinetic energy of the rotor system is made up of the energy of rotating ${T}_{1}$ and non-rotating parts of the rotor ${T}_{2}$, as well as the kinetic energy of the balls ${T}_{m}$:
where ${T}_{1}$ – kinetic energy of the rotating parts of the rotor; ${T}_{2}$ – non-rotating parts of the rotor; ${T}_{m}$ – kinetic energy of the balls.
The kinetic energy of the rotating parts of the rotor has two components, caused by the translational movement of the rotor with the center of masses (point ${C}_{1}$) and the rotational movement of the rotor around the center of masses [6], and is expressed by the equation:
where ${V}_{{C}_{1}}$ – absolute velocity of the center of masses of the rotor; ${T}_{1rot}$ – kinetic energy of the rotational motion of the rotor relative to its center of masses.
The position of the center of masses of rotor ${C}_{1}$ in the fixed coordinate system $0xyz$ is determined by the expression:
where $\stackrel{-}{r}$ – radius vector that determines the position of the origin of the moving coordinate system ${0}_{1}{x}_{1}{y}_{1}{z}_{1}$; ${\stackrel{-}{r}}_{1}$ – radius vector that determines the position of the center of masses of the rotor in the moving coordinate system ${0}_{1}{x}_{1}{y}_{1}{z}_{1}$ in relative motion.
The absolute velocity of the center of the rotor masses after transformation is determined by:
where:
where $\stackrel{-}{\mathrm{\Omega}}$ – vector of the angular velocity of the rotor in the moving coordinate system ${0}_{1}{x}_{1}{y}_{1}{z}_{1}$_{}when it moves relative to pole ${0}_{1}$.
The transfer matrix for the transition to the fixed coordinate system $0xyz$ is obtained by multiplying the rotation matrices of the moving coordinate system ${0}_{1}{x}_{1}{y}_{1}{z}_{1}$ relative to the fixed one:
Taking into account the smallness of the rotor vibration angles, rotation matrix ${V}_{1}$ takes the form:
Multiplying Eqs. (10) and (12) by the rotation matrix ${V}_{1}$ Eqs. (14), taking into account Eqs. (9), we obtain in the final form the velocity vector of the center of masses of the rotor:
$+\stackrel{-}{k}\left(\dot{z}+e\dot{\alpha}\mathrm{sin}\phi -e\dot{\beta}\mathrm{cos}\phi \right).$
The expression for the kinetic energy of the rotor in relative motion is obtained taking into account the transition matrix from the central axes of the rotor ${C}_{1}{x}_{2}{y}_{2}{z}_{2}$ to the main central axes ${C}_{1}{x}_{2}^{\text{'}}{y}_{2}^{\text{'}}{z}_{2}^{\text{'}}$. Due to the fact that the inertia moments of the rotor relative to the axes ${C}_{1}{x}_{2}$ and ${C}_{1}{y}_{2}$ are equal to each other, the couple unbalance of the rotor is determined by angle $\delta $, and the proper rotation of the rotor is determined by angle $\phi $, the rotation matrices have the form:
Transition matrix ${V}_{2}$ is obtained by multiplying rotation matrices ${V}_{\phi}$ and ${V}_{\delta}$ and, taking into account the smallness of angle $\delta $, it takes the following form:
The angular velocity vector, taking into account its own rotation in the moving coordinate system ${C}_{1}{x}_{2}{y}_{2}{z}_{2}$, has the following form:
Multiplying Eqs. (18) by (17), we obtain the vector of the angular rotor velocity when it moves relative to the center of masses in a moving coordinate system associated with the main central axes of inertia ${C}_{1}{x}_{2}^{\text{'}}{y}_{2}^{\text{'}}{z}_{2}^{\text{'}}$, i.e. ${\stackrel{~}{\mathrm{\Omega}}}_{1}={V}_{2}{\stackrel{-}{\mathrm{\Omega}}}_{1}.$ Then, the kinetic energy of the rotor in relative motion is determined by the expression:
The total kinetic energy ${T}_{1}$ of the rotor, taking into account Eqs. (6) and (19), takes the form:
$\left.+{\left(\dot{z}+e\dot{\alpha}\mathrm{sin}\phi -e\dot{\beta}\mathrm{cos}\phi \right)}^{2}\right]+\mathrm{0,5}{I}_{x1}\left({\dot{\alpha}}^{2}\right.+{\dot{\beta}}^{2}+2\delta \dot{\phi}\dot{\alpha}\mathrm{cos}\phi +2\delta \dot{\theta}\dot{\alpha}\mathrm{cos}\phi $
$\left.-2\delta \dot{\phi}\dot{\beta}\mathrm{sin}\phi -2\delta \dot{\theta}\dot{\beta}\mathrm{sin}\phi +2\beta \dot{\alpha}\dot{\theta}-2\alpha \dot{\beta}\dot{\theta}\right)+\mathrm{0,5}{I}_{z1}\left({\dot{\phi}}^{2}\right.+{\dot{\theta}}^{2}+2\delta \dot{\phi}\dot{\beta}\mathrm{sin}\phi $
$\left.-2\delta \dot{\phi}\dot{\alpha}\mathrm{cos}\phi +2\delta \dot{\theta}\dot{\beta}\mathrm{sin}\phi -2\delta \dot{\theta}\dot{\alpha}\mathrm{cos}\phi +2\alpha \dot{\beta}\dot{\phi}-2\beta \dot{\alpha}\dot{\phi}+2\dot{\phi}\dot{\theta}+2\alpha \dot{\beta}\dot{\theta}-2\beta \dot{\alpha}\dot{\theta}\right),$
where ${I}_{x1}={I}_{y1}$ – moment of inertia of the rotor relative to axes ${x}_{1}$ or ${y}_{1}$; ${I}_{z1}$ – moment of inertia of the rotor relative to ${z}_{1}$_{}axis.
The kinetic energy of ball ${T}_{{m}_{j}}$ is determined in a similar way. The ball participates simultaneously in two motions in the circumferential and cross sections [3]:
- in a translational motion of the ball with a center of masses (point ${C}_{{m}_{j}}$)
– in a rotational motion around the center of masses.
The design model for determining the kinetic energy of the weights is shown in Fig. 2:
where ${\dot{\omega}}_{j}=\dot{\phi}+{\dot{\psi}}_{j}$ – absolute velocity of the ball in the circumferential plane; ${\dot{\omega}}_{j}={\dot{\gamma}}_{j}$– absolute velocity of the ball in the transverse plane.
Fig. 2. Rolling pattern without sliding the ball on the ball raceway; a) in the circumferential plane, b) in the transverse plane
a)
b)
If we neglect the displacements of rotor $z$ along the axis of rotation and rotational vibration $\theta $ around the same axis, due to their smallness and negligible effect on the transverse vibrations of the rotor, the kinetic energy of the entire rotor system takes the following form:
$\left.-2\delta \dot{\phi}\dot{\beta}\mathrm{sin}\phi \right)+\mathrm{0,5}{I}_{z1}\left({\dot{\phi}}^{2}\right.+2\delta \dot{\phi}\dot{\beta}\mathrm{sin}\phi -2\delta \dot{\phi}\dot{\alpha}\mathrm{cos}\phi \left.+2\alpha \dot{\beta}\dot{\phi}-2\beta \dot{\alpha}\dot{\phi}\right)$
$+\mathrm{0,5}{m}_{2}\left({\dot{x}}^{2}+{\dot{y}}^{2}\right)+\mathrm{0,5}{m}_{m}{\sum}_{1}^{n}\left\{{\dot{x}}^{2}+{\dot{y}}^{2}+{R}^{2}\left[{\left(\dot{\phi}+{\dot{\psi}}_{j}\right)}^{2}+{\dot{\beta}}^{2}{\mathrm{cos}}^{2}\left(\phi +{\psi}_{j}\right)\right.\right.$
$+{\dot{\alpha}}^{2}{\mathrm{sin}}^{2}\left(\phi +{\psi}_{j}\right)-2\beta \dot{\beta}\left(\dot{\phi}+{\dot{\psi}}_{j}\right)\mathrm{sin}\left(\phi +{\psi}_{j}\right)\mathrm{cos}\left(\phi +{\psi}_{j}\right)-2\dot{\alpha}\dot{\beta}\mathrm{sin}\left(\phi +{\psi}_{j}\right)$
$\times \mathrm{cos}\left(\phi +{\psi}_{j}\right)+2\alpha \dot{\alpha}\left(\dot{\phi}+{\dot{\psi}}_{j}\right)\mathrm{sin}\left(\phi +{\psi}_{j}\right)\mathrm{cos}\left(\phi +{\psi}_{j}\right)+2\beta \dot{\alpha}\left(\dot{\phi}+{\dot{\psi}}_{j}\right)$
$\times {\mathrm{sin}}^{2}\left(\phi +{\psi}_{j}\right)\left.-2\alpha \dot{\beta}\left(\dot{\phi}+{\dot{\psi}}_{j}\right){\mathrm{cos}}^{2}\left(\phi +{\psi}_{j}\right)\right]+2R\left(\dot{\phi}+{\dot{\psi}}_{j}\right)\left[\dot{y}\mathrm{cos}\left(\phi +{\psi}_{j}\right)\right.$
$\left.-\dot{x}\mathrm{sin}\left(\phi +{\psi}_{j}\right)\right]+2\left({R}_{1}-r\right){\dot{\gamma}}_{j}\left[\dot{x}\mathrm{cos}\left(\phi +{\psi}_{j}\right)+\dot{y}\mathrm{sin}\left(\phi +{\psi}_{j}\right)\right]\mathrm{cos}{\gamma}_{j}$
$\left.+{\left({R}_{1}-r\right)}^{2}{\dot{\gamma}}_{j}^{2}\right\}+\mathrm{0,2}{m}_{m}{\sum}_{1}^{n}\left[{\dot{\psi}}_{j}^{2}{R}^{2}+{\left({R}_{1}-r\right)}^{2}{\dot{\gamma}}_{j}^{2}\right].$
When determining the potential energy of the system, its change due to the vertical displacement of the center of the rotor masses is not taken into account. In addition, when determining the potential energy of the balls, the angle of rotation of rotor axis $\alpha $ and $\beta $ is not taken into account due to its smallness in comparison with angle ${\gamma}_{j}$. The value of the potential energy of the rotor system is determined by the deformation of the elastic elements of the rotor suspension and the change in the position of the center of masses of the compensating weights (balls):
where ${c}_{i}$ – stiffness of $i$ elastic element along its axis; ${\mathrm{\Delta}}_{i}$ – value of deformation of $i$ elastic element; $r$ – radius of the ball; ${R}_{1}$ – radius of the torus cross-section.
The impact of dissipative forces in the elastic-dissipative suspension of the rotor system is taken into account by introducing equivalent viscous friction. According to the hypothesis of dissipative forces, the resistance of the medium is taken proportional to the generalized velocities of translational and angular displacements. In accordance with the hypothesis, the dissipative function in the elastic suspension has the following form:
where ${b}_{i}$– dissipation coefficient of $i$ suspension element; ${\dot{q}}_{i}$ – value of the generalized velocity of $i$ element.
Dissipative forces in the ABD can be modelled through the forces of dry and / or viscous rolling friction. It is assumed that the ball, when moving relative to the body, moves without slipping and without separation from the body. In this case, the moment of friction forces acting on the ball is equal to:
– dry rolling friction:
where ${N}_{j}$ – reaction force from the pressure of $i$ ball on the inner surface of the torus; $k$ – rolling friction coefficient;
– viscous friction:
where ${b}_{r}$_{}– dissipation coefficient for viscous friction in the ABD.
Force ${N}_{j}$ without taking into account the accelerations, which are small in comparison with the centripetal acceleration of the ball, is determined from the expression:
After transformations, in accordance with Eqs. (4), we obtain a system of differential equations for a rotor system with a vertical axis of rotation with a ball-type ABD:
1) For $x,y$:
$\left.+\left(\ddot{\phi}+{\ddot{\psi}}_{j}\right)\mathrm{cos}\left(\phi +{\psi}_{j}\right)\right]+\left({R}_{1}-r\right)\left({\ddot{\gamma}}_{j}\mathrm{cos}\gamma -{\dot{\gamma}}_{j}^{2}\mathrm{sin}{\gamma}_{j}\right)\mathrm{cos}\left(\phi +{\psi}_{j}\right)$
$\left.-2\left({R}_{1}-r\right){\dot{\gamma}}_{j}\left(\dot{\phi}+{\dot{\psi}}_{j}\right)\mathrm{sin}\left(\phi +{\psi}_{j}\right)\mathrm{cos}{\gamma}_{j}\right\}={m}_{1}e\left({\dot{\phi}}^{2}\mathrm{cos}\phi +\ddot{\phi}\mathrm{sin}\phi \right),$
$\left.-\left(\ddot{\phi}+{\ddot{\psi}}_{j}\right)\mathrm{sin}\left(\phi +{\psi}_{j}\right)\right]+\left({R}_{1}-r\right)\left({\ddot{\gamma}}_{j}\mathrm{cos}{\gamma}_{j}-{\dot{\gamma}}_{j}^{2}\mathrm{sin}{\gamma}_{j}\right)\mathrm{sin}\left(\phi +{\psi}_{j}\right)$
$\left.+2\left({R}_{1}-r\right){\dot{\gamma}}_{j}\left(\dot{\phi}+{\dot{\psi}}_{j}\right)\mathrm{cos}\left(\phi +{\psi}_{j}\right)\mathrm{cos}{\gamma}_{j}\right\}={m}_{1}e\left({\dot{\phi}}^{2}\mathrm{sin}\phi -\ddot{\phi}\mathrm{cos}\phi \right)$
2) For $\alpha $, $\beta $:
$\times \mathrm{sin}\left(\phi +{\psi}_{j}\right)\mathrm{cos}\left(\phi +{\psi}_{j}\right)-{\left(\dot{\phi}+{\dot{\psi}}_{j}\right)}^{2}\left[\alpha \mathrm{cos}2\left(\phi +{\psi}_{j}\right)+\beta \mathrm{sin}2\left(\phi +{\psi}_{j}\right)\right]$
$\left.-\left[\left(\ddot{\phi}+{\ddot{\psi}}_{j}\right)\beta +\left(\dot{\phi}+{\dot{\psi}}_{j}\right)\dot{\beta}-\ddot{\alpha}\right]{\mathrm{sin}}^{2}\left(\phi +{\psi}_{j}\right)\right\}=\delta \left({I}_{x1}-{I}_{z1}\right)\left(\ddot{\phi}\mathrm{sin}\phi +{\dot{\phi}}^{2}\mathrm{cos}\phi \right),$
$\times \mathrm{sin}\left(\phi +{\psi}_{j}\right)\mathrm{cos}\left(\phi +{\psi}_{j}\right)-{\left(\dot{\phi}+{\dot{\psi}}_{j}\right)}^{2}\left[\alpha \mathrm{sin}2\left(\phi +{\psi}_{j}\right)-\beta \mathrm{cos}2\left(\phi +{\psi}_{j}\right)\right]$
$\left.+\left[\left(\ddot{\phi}+{\ddot{\psi}}_{j}\right)\alpha +\left(\dot{\phi}+{\dot{\psi}}_{j}\right)\dot{\alpha}+\ddot{\beta}\right]{\mathrm{cos}}^{2}\left(\phi +{\psi}_{j}\right)\right\}=\delta \left({I}_{y1}-{I}_{z1}\right)\left(\ddot{\phi}\mathrm{cos}\phi -{\dot{\phi}}^{2}\mathrm{sin}\phi \right).$
3) For ${\psi}_{j}$, ${\gamma}_{j}$:
$\left.+\ddot{y}\mathrm{cos}\left(\phi +{\psi}_{j}\right)\right]+R\left[\left(\alpha \ddot{\alpha}-\beta \ddot{\beta}\right)\mathrm{sin}\left(\phi +{\psi}_{j}\right)\mathrm{cos}\left(\phi +{\psi}_{j}\right)+\beta \ddot{\alpha}{\mathrm{sin}}^{2}\left(\phi +{\psi}_{j}\right)\right.$
$\left.\left.-\alpha \ddot{\beta}{\mathrm{cos}}^{2}\left(\phi +{\psi}_{j}\right)\right]\right\}+{M}_{{\psi}_{j}}=0,j=1,n,$
$-R{\left(\dot{\phi}+{\dot{\psi}}_{j}\right)}^{2}\mathrm{cos}{\gamma}_{j}+g\mathrm{sin}{\gamma}_{j}-R\left[{\dot{\beta}}^{2}{\mathrm{cos}}^{2}\left(\phi +{\psi}_{j}\right)-{\dot{\alpha}}^{2}{\mathrm{sin}}^{2}\left(\phi +{\psi}_{j}\right)\right]\mathrm{cos}{\gamma}_{j}$
$+2R\left(\dot{\phi}+{\dot{\psi}}_{j}\right)\left[\dot{\beta}\mathrm{cos}\left(\phi +{\psi}_{j}\right)-\dot{\alpha}\mathrm{sin}\left(\phi +{\psi}_{j}\right)\right]\left[\beta \mathrm{sin}\left(\phi +{\psi}_{j}\right)-\alpha \mathrm{cos}\left(\phi +{\psi}_{j}\right)\right]$
$\left.\times \mathrm{cos}{\gamma}_{j}+2R\dot{\alpha}\dot{\beta}\mathrm{sin}\left(\phi +{\psi}_{j}\right)\mathrm{cos}\left(\phi +{\psi}_{j}\right)\mathrm{cos}{\gamma}_{j}\right\}+{M}_{{\gamma}_{j}}=0,j=1,n.$
4) For $\phi $. With a sufficiently large value of the engine torque ${M}_{rot}$ on the rotor shaft, in a first approximation, the following conditions for changing the angular velocity can be accepted:
where ${I}_{z1}$ – moment of inertia of the rotor system relative to the axis of rotation; ${\omega}_{0}$ – operating angular velocity of the rotor.
The system of differential Eqs. (28-34) describes the motion of the rotor system in transient and steady-state modes.
This system of equations can be used to study other rotor systems with ball-type ABDs, in particular, a symmetric rotor on two supports with a symmetric or asymmetric suspension [8, 13]. At the same time, comparing the system of differential equations of a symmetric vertical rotor with the system of differential equations proposed in this work, one can state their complete coincidence provided there is exclusion from the system of angular coordinates for the rotor axis.
As an example, a system of differential equations for the design model of the rotor from [8] is given, which is obtained from the system of differential Eqs. (28-34) under the following conditions: $\alpha =\beta =\dot{\alpha}=\dot{\beta}=\ddot{\alpha}=\ddot{\beta}=0$; $\ddot{\phi}=0$; $j=1.$ Fig. 3 shows a design model of a vertical symmetrical rotor with a torus-shaped automatic balancer:
$+\left({R}_{1}-r\right)\left(\ddot{\gamma}\mathrm{cos}\gamma -{\dot{\gamma}}^{2}\mathrm{sin}\gamma \right)\mathrm{cos}\left(\phi +\psi \right)\left.-2\left({R}_{1}-r\right)\dot{\gamma}\left(\dot{\phi}+\dot{\psi}\right)\mathrm{sin}\left(\phi +\psi \right)\mathrm{cos}\gamma \right\}$
$={m}_{1}e\left({\dot{\phi}}^{2}\mathrm{cos}\phi \right),$
$+\left({R}_{1}-r\right)\left(\ddot{\gamma}\mathrm{cos}\gamma -{\dot{\gamma}}^{2}\mathrm{sin}\gamma \right)\mathrm{sin}\left(\phi +\psi \right)\left.+2\left({R}_{1}-r\right)\dot{\gamma}\left(\dot{\phi}+\dot{\psi}\right)\mathrm{cos}\left(\phi +\psi \right)\mathrm{cos}\gamma \right\}$
$={m}_{1}e\left({\dot{\phi}}^{2}\mathrm{sin}\phi \right),$
$+{M}_{\psi}=0,$
$-{m}_{m}R\left({R}_{1}-r\right){\left(\dot{\phi}+\dot{\psi}\right)}^{2}\mathrm{cos}\gamma +{m}_{m}\left({R}_{1}-r\right)g\mathrm{sin}\gamma +{M}_{\gamma}=0.$
On the basis of the system of differential Eqs. (35-38), in particular, numerical calculations of the effect of the ball size on its acceleration in the autobalancing mode were carried out using the SPRING software [12], which is confirmed by field experiments [6, 22].
Fig. 3. Design model of the rotor system; a) general view; b) torus cross section
a)
b)
3. Conclusions
The authors obtained a system of differential equations for a rotor system with a ball-type ABD, in which the inner space of the body is made in the shape of a torus, which provides a minimum resistance force to the motion of the ball and the best conditions for the location of balls in the auto-balancing mode.
The proposed system of differential equations makes it possible to proceed to the study of simpler rotor systems with a vertical axis of rotation, equipped with an ABD, in which the balls move in both longitudinal and transverse directions. The authors carried out a physical experiment with such a rotor system and a vertical axis of rotation. The results of the experiment completely coincided with mathematical modelling in the SPRING environment, thereby confirming the correctness of the proposed differential equations.
The system of differential equations makes it possible to investigate both steady-state modes of rotor motion and transient modes, as well as to study the stability areas of various modes of movement of ABD weights and select the parameters of the rotor system.
It is possible to model the force of resistance to the motion of the ball both by dry rolling friction and by viscous friction, which depends on the velocity of the ball.
References
- Gorbenko A., Strautmanis G., Filimonikhin G., Mezītis M. Motion Modes of the Nonlinear Mechanical System of the Rotor Autobalancer. Vibroengineering Procedia, Vol. 25, Issues 1, 2019, p. 1-6. [Publisher]
- Nesterenko V. P. Automatic elimination by balls of static imbalance of the rotor of a two-mass system. University News. Engineering, Vol. 3, 1983, p. 46-50, (in Russian). [Search CrossRef]
- Strautmanis G., Mezitis M., Strautmane V., Gorbenko A. Model of a vertical rotor with an automatic balancer with two compensating masses. Vibroengineering Procedia, Vol. 21, Issue 1, 2018, p. 202-207. [Publisher]
- Strauch D. Classical Mechanics. An Introduction. Springer-Verlag, Berlin Heidelberg, 2009. [Search CrossRef]
- Gorbenko A. N., Shmelev S. Kh. Necessary self-balancing robustness conditions for a two-bearing rotor taking unbalance mass into account. Herald of the Bauman Moscow State Technical University. Series Mechanical Engineering, Vol. 5, Issue 122, 2018, p. 36-50. [Publisher]
- Gorbenko A., Mezītis M., Strautmane V., Strautmanis G. The impact of an elastic rotor suspender and the size of the compensating mass on the acceleration of the automatic balancer. Procedia Computer Science, Vol. 149, Issue 2019, 2018, p. 301-306. [Search CrossRef]
- Goncharov V. V., Filimonikhin G. B. Form and structure of differential equations of motion and process of auto-balancing in the rotor machine with auto-balancers. Bulletin of the Tomsk Polytechnic University, Geo Assets Engineering, Vol. 326, Issue 12, 2015, p. 20-30. [Search CrossRef]
- Strautmanis G., Mezitis M., Strautmane V. Model of a vertical rotor with a ball-type automatic balancer. Vibroengineering Procedia, Vol. 8, 2016, p. 57-62. [Search CrossRef]
- Gorbenko A. N., Klimenko N. P., Strautmanis G. Influence of rotor unbalance increasing on its autobalancing stability. Procedia Engineering, Vol. 206, 2017, p. 266-271. [Publisher]
- Sperling L., Ryzhik B., Duckstein H. Single-plain auto-balancing of rigid rotors. Technische Mechanik, Vol. 24, Issue 1, 2004, p. 1-24. [Search CrossRef]
- Sperling L., Ryzhik B., Duckstein H. Single-plain auto-balancing of anisotropically supported rigid rotors. Technische Mechanik, Vol. 24, Issue 1, 2004, p. 37-50. [Search CrossRef]
- Ščukins I., Zakrževskis M., Ivanov Y., et. al. Application of software SPRING and method of complete bifurcation groups for the bifurcation analysis of nonlinear dynamical system. Journal of Vibroengineering, Vol. 10, Issue 4, 2008, p. 510-518. [Search CrossRef]
- Strautmanis G., Mezītis M., Strautmane V., Gorbenko A. On the issue of impact of anisotropy of the rotor elastic suspension on the performance of the automatic balancer. Vibroengineering Procedia, Vol. 17, Issue 1, 2019, p. 1-6. [Publisher]
- Evseev D., Medvedev B., Medvedev P., Strautmanis G., Samoshkin S. Acoustic emission approach to determining survivability in fatigue tests. Procedia Computer Science, Vol. 149, 2019, p. 282-287. [Publisher]
- Yatsun V., Filimonikhina I., Podoprygora N., Hurievska O. Motion equations of the single-mass vibratory machine with a rotary-oscillatory motion of the platform and the form of a passive auto-balancer. Eastern-European Journal of Enterprise Technologies, Vol. 6, Issues 7(96), 2018, p. 58-67. [Publisher]
- Goncharov V. V., Filimonikhin G. B. Form and structure of differential equations of motion and process of autobalancing in the rotor machine with auto-balancers. Bulletin of the Tomsk Polytechnic University. Georesource Engineering, Vol. 326, Issue 12, 2015, p. 19-30, (in Russian). [Search CrossRef]
- Filimonikhin G., Yatsun V., Filimonikhina I., Ienina I., Munshtukov I. Studyng the load jam modes within the frame work of a flat model of the rotor with an autobalancer. Eastern-European Journal of Enterprise Technologies, Vol. 5, Issues 7(101), 2019, p. 51-61. [Publisher]
- Bykov V. G., Kovachev A. S. Dynamic of statically unbalanced rotor with eccentric ball autobalancer. Bulletin of St. Petersburg State University, Vol. 1, Issue 1, 2014, p. 579-588, (in Russian). [Search CrossRef]
- Green K., Friswell M. I., Champneys A. R., Lieven N. A. J. The stability of automatic ball balancers. 7th IFToMM-Conference on Rotor Dynamics, Vienna, Austria, 2006. [Search CrossRef]
- Adolfsson J. Passive Control of Mechanical Systems. Bipedal Walking and Autobalancing. Department of Mechanics. Royal Institute of Technology, S-10044 Stockholm, Sweden, 2001. [Search CrossRef]
- Ding H.,Li M., Fu H., Li M. The impact of rotor elastic suspension setting on the acceleration of the automatic balancer compensating mass. Recent Patents on Mechanical Engineering, Vol. 11, Issue 3, 2018, p. 232-241. [Publisher]
- Yatsun V., Filimonikhin G., Haleeva A., Krivoblotsky L., Machok Y., Mezitis M., Podoprygora N., Sadovyi M., Strautmanis G. Searching for the two frequency motion modes of a three-mass vibratory machine with a vibration exciter in the form of a passive auto-balancer. Eastern-European Journal of Enterprise Technologies, Vol. 4, Issues 7(106), 2020, p. 103-111. [Publisher]