Research Article  Open Access
Alexey V. Kornaev, Nickolay V. Kornaev, Elena P. Kornaeva, Leonid A. Savin, "Application of Artificial Neural Networks to Calculation of Oil Film Reaction Forces and Dynamics of Rotors on Journal Bearings", International Journal of Rotating Machinery, vol. 2017, Article ID 9196701, 11 pages, 2017. https://doi.org/10.1155/2017/9196701
Application of Artificial Neural Networks to Calculation of Oil Film Reaction Forces and Dynamics of Rotors on Journal Bearings
Abstract
Increase of energy efficiency and level of information system development of rotor machines in general requires improvement of theoretical approaches to research. In the present paper the problem of highprecision and highperformance computing programs development has been considered to simulate rotor vibrations. Based on twolayer feedforward neural networks, numerical models have been developed to calculate oil film reaction forces to solve the rotor dynamics problems. Comparison has been done of linear and nonlinear approaches to solution of rotor dynamics problems, and a qualitative evaluation has been presented of accuracy and performance of a neural network approach compared to conventional approaches to rotor dynamics.
1. Introduction
Dynamic behavior of a rotating shaft on journal bearings is defined by the resonance vibration at critical speeds, selfexcited vibrations due to the internal damping of the shaft’s material, selfexcited vibrations due to the oil film of journal bearings, oil whip, and flowinduced vibrations [1–4]. Since the authors are closely connected with hydrodynamic lubrication, this paper is dedicated to oil whirl and oil whip problems. Oil film forces are obtained by means of solving the hydrodynamics equations and depend on many design and operational parameters of a bearing as well as on various specific physical effects that may occur, such as vibration, cavitation, and turbulence [2, 4, 5]. The number of input parameters can reach several dozen [6]. At the same time, development of technology in the direction of informatization requires simultaneous increase of accuracy and performance of the numerical models. One of the possible ways to provide both accuracy and high performance is using deterministic and stochastic approaches to modeling. Nowadays the methods and resources of implementation of neural network algorithms are developing quite intensively, which, based on the known data about the subject of study, allows building multifactor models of input and output parameters relations. Three basic levels of developed models could be distinguished.
Models of the first level have no direct interaction with the rotating machine; their training is based on the given numerical or experimental data. In [7] the results of study of an elastic rotor’s dynamics on rigid bearings have been presented. Neural networks that have been applied are feedforward with two hidden layers. High accuracy has been demonstrated for models of rotor systems with rollerelement bearings based on neural networks: accuracy from 82% up to 100% has been presented. Dynamics has been studied in a certain range of speed and for every case of speed separately, so the neural network had to be retrained. Some works are known, where neural network models are used to calculate the distributed characteristics, for example, pressure distribution in the oil film of the journal bearings [8, 9]. It should be noted that application of stochastic methods of rotor dynamics modeling allows skipping the stage of distributed values calculation. Instead, one could instantly determine the components of a resulting oil film force. In [10] dynamics of an imbalanced rigid rotor on fluidfilm bearings has been studied. Authors determined the relation between three kinematic characteristics and two components of the resulting oil film force using a model of a neural network. The numerical tests were carried out for various angular speed and allowed determination of rotor trajectories. The synchronous, periodic, quasiperiodic, and chaotic motions for rotation speed within the range from 1000 to 10000 rpm were presented as the results of calculation. The authors suggest that there is a limitation to using linearized models in solving the problems of dynamics; however, they do not present a quantitative comparison of the numerical results obtained using linear and nonlinear models. Moreover, the reduction of the number of input variables has been done according to Poincare transformation of Reynolds equation, but possible outcomes of such substitution were not discussed in the paper. It is also not quite clear how one introduced kinematic variable allowed accounting for three given variables (two projections of rotor’s lateral vibration speed and its rotating speed). It could be that the neural network had to be retrained for every particular case of the rotating speed. So, the results of this indeed interesting and important work remain in need of additional verification. A number of quite perspective papers should also be highlighted where neural network algorithms are used to account for complex specific phenomena [5, 11] and to optimize the operational conditions of fluidfilm bearings [12]. Models of higher level can learn in real time using the processed measurement data. One of the most popular applications of such models is using neural networks to flaw detection of bearings in rotor machines [13–16]. In [17] authors investigate the inverse model of the squeezefilm damper bearing (force input/displacement output) from empirical data and its application to a nonlinear inverse rotorbearing problem. Comparison of experimental and numerical data showed high accuracy of the model based on recurrent neural networks. Models of highest level of complexity are able to generate some control impact [18, 19]. Today such approach is used to solve the problems of control in magnetic bearings [20–22]. It should be, however, noted that the values of clearance in such bearings are an order of magnitude more than the same value in fluidfilm bearings.
The present paper considers development of neural network algorithms of first level, that is, models with no direct contact to the rotor machine. Matters of quantitative evaluation of accuracy of application of neural network in comparison to known methods of oil film forces calculation and applicability of neural networks during simulation of nonlinear selfexcited oscillations—oil whip—and during simulation of transient regimes of rotor’s operation have been considered.
2. Mathematical Model
A rigid symmetrical imbalanced rotor is considered. Both ends of the rotor are supported by journal bearings of the same specification (see Figure 1). Given the chosen direction of coordinate axes, the pressure distribution and the resulting oil film force can be determined by means of numerically solving the Reynolds equation for a case of an isothermal flow of a viscous incompressible fluid [6]:where is the oil film thickness (radial gap) function, is the dynamic viscosity coefficient (assumed to be constant), is the pressure function, and is accordingly the tangential and normal components of the journal’s surface velocity, .
Radial gap function could be determined geometrically (see Figure 1) (coordinates and velocities of the center of the rotor are indicated with the capital letters, while the coordinates and the velocities of the fluid medium are indicated with the lower case letters):where is the mean oil film thickness, are coordinates of the center of the journal, is the angular coordinate, and is the bearing’s diameter.
The values of over the surface of the journal could be determined using the following kinematic relations:where is angular speed of the rotor and are projections of the velocity vector of the center of the rotor.
Calculation of pressure distribution at some specific point in time is a boundary problem of solution of Reynolds equation (1) considering (2) and (3) with the following boundary conditions: the pressure on each side of the bearing is known to be and ; along the direction the following condition has to be met: ; . Additionally, Gumbel’s boundary condition for film rupture is used [6, 23].
Given the pressure distribution, the components of the resulting oil film reaction force can be obtained as projections on the (see Figure 1):where is the length of the bearing.
The equation for lateral vibrations of the rotor’s center can be presented in simple form [1–4, 6]:where is the mass of the rotor per one bearing, is time, is the bearing load (centrifugal and inertia forces), is the rotor’s imbalance, and is free fall acceleration.
Oscillations’ velocity of the rotor is significantly less than tangential velocity of journal’s surface , so the solution of (5) is linked to solution of the steadystate fluid’s flow problem solution. So, time as variable could be neglected during calculation of oil film force (4). This significantly simplifies the procedure of numerical solution of rotor dynamics problems [1, 6]. Three basic approaches have been considered and are characterized by different methods of oil film forces determination.
The first approach, called trajectory method [1–4, 6], is based on solution of Reynolds equation (1) followed by integration of pressure distribution and determination of (4). Such approach is considered referential due to high precision and however has low performance in computational terms, since the Reynolds equation has to be solved at each time step.
The second approach, called the method of linearization, of the oil film force determination [1–4, 6] is based on the decomposition of reaction force components into the Taylor series in the proximity of the equilibrium position . The matrices of stiffness and damping , calculated in the equilibrium position, shall be introduced as follows:Then the motion equation of the rotor (5) with given operational conditions of the rotor system and given the linearization of the film’s reaction forces will be as follows [1–4, 6, 24]:The equilibrium point can be determined using the trajectory approach. In many practical applications the equilibrium locus curve is calculated to solve problems of rotor dynamics [1–4, 6]. The points on this curve correspond to the equilibrium positions at various rotational speeds. Necessity to determine the equilibrium position or the equilibrium locus curve complicates the method of linearization.
The third approach is based on a neural network program module to approximate the oil film force’s components and to solve (5). A twolayer feedforward neural network with inputs () and outputs () is used (see Figure 2). The functionality of the network is determined by specifying the strengths of the connections, called weights , and the neuron transfer function [25]. The mathematical model of the used neural network is where is a linear transfer function and is a sigmoid symmetric transfer function, and are weights, and are thresholds of output and hidden layers, respectively, and are the intermediate outputs of hidden neurons ().
In the present paper the components of the oil film force are considered as outputs, and four or five kinematic values, depending on the problem, are considered inputs: or , where is the rotor’s speed.
Neural network learning is essentially the solution of optimization problem of a minimum of the mean squared error function search:where is a number of samples of inputs and outputs and is the target value of an th output in the th sample.
The mean squared error function (9) is the function of a number of variables and . The search for an extremum of such function is implemented numerically. One of the criteria of neural network’s quality is the correlation coefficient between outputs and targets [26]. Target data in the present paper is based on the numerical results, obtained using the trajectory approach.
Stability analysis of the obtained results could be implemented based on the rotor’s locus analysis using Poincare maps [10, 27–29]. Moreover, the Routh–Hurwitz criterion could be applied to the linearization approach [1, 10].
3. Simulation Results
Three basic approaches have been considered and characterized by different methods of oil film forces determination: the trajectory method, the linearization method, and the neural network method. Numerical solution of Reynolds equation (1) was implemented using the finite difference method. The obtained set of linear equations was solved using the Gauss method [30]. The LevenbergMarquardt backpropagation algorithm was used in training neural networks [25, 31, 32], which, in its turn, was based on the numerical results, obtained using the trajectory method. The rotor dynamics calculation program was based on the numerical solution of (5) or (7) using the explicit RungeKutta method of 4th and 5th order [30]. During the development of the programs the GNUOctave software and Octave’s Neural Network Package [31] have been used. The GNUOctave products are free software and use the Matlab compatible language. During simulations a PC with following characteristics has been used: Intel® Core™ i54690K CPU 3.5 GHz, RAM 12 GB, x64 OS.
The results presented below correspond to various types of problems solved.
3.1. Rotor Dynamics with Fixed Speed and Low Vibrations
A rotor system with the following parameters has been subjected to the study: rotor with 8.21 kg mass rotates at a constant speed of 4000 rpm under the load and centrifugal force and is supported by two plain journal fluidfilm bearings with 40 mm diameter and 20 mm length with a mean radial gap of 100 μm; lubricant is mineral oil with dynamic viscosity of 13 mPa·s; the pressure on each side of the bearing is 0.1 MPa.
First, trajectory approach has been applied. A preliminary numerical test neglecting the centrifugal force successfully converged to the equilibrium point . Next, linearization approach was used.
Then, a neural network has been developed for training with a number of neurons in the hidden layer . As the input data a matrix of elements was used, which consisted of numerical values of the coordinates and the velocities , with s calculated time of rotor’s motion. The target data matrix of elements consisted of the numerical values of oil film forces . The matrices of input and output data were divided into three groups for training, validation and testing in fractions 0.6, 0.3, and 0.1, respectively. The procedure of adjustment has been implemented using application software [31]. Training was done in 100 iterations, socalled epochs [25], in approximately 10 minutes. The convergence graph of a mean squared error of the network during training is presented in Figure 3. Training quality of a neural network is characterized by the following values of the mean squared error’s value and correlation coefficients: N, .
The numerical tests for comparative calculation of the rotor’s trajectories have been carried out for various values of the relative centrifugal force . Other input parameters remained as during the training procedure.
As a result of numerical tests a calculation time of 1 sec of rotor motion was determined given the action of dimensionless centrifugal force (see Figure 4). It could be seen that the calculation time of rotor dynamics calculation using the linearization approach is the least in the whole range. Application of the neural network module increased the machine time up to 1.5–2 times under given conditions. Application of trajectory approach increased the calculation time up to 60–90 times compared to the linearization approach. It is indicative that calculation time of rotor dynamics using linearization approach and neural network approach is faster than the flow of real time in the whole range. This fact means that the application of linearization or neural network approaches in a real rotational machine control system would make it possible to obtain results of prognostic modeling [20, 33, 34].
As a result of another series of numerical tests, the trajectories have been calculated for a dimensionless centrifugal force (see Figure 5). Calculated rotor motion’s time was set to 1 sec. The simulation started at the equilibrium point , the trajectory approach taken as referential. The relative error of calculation has been chosen as quantitative criteria of accuracy evaluation:where are coordinates, obtained using the trajectory approach, and are coordinates, obtained using the linearization or neural network approach, .
The quality analysis of stability of calculated trajectories has been implemented based on the Poincare maps [27–29]. The point were displayed at moments of time that were multiples of the period of rotation of the rotor . It could be seen from the Poincare maps (see Figure 5) that in the given range of loads the rotor’s trajectories are quasiperiodic. This results from the fact that after a short transient process the trajectory’s mapping becomes a dense pointed line, which closes if the calculation time increases by 100 times or more. This shows the presence of incommensurable frequencies in the oscillations, while the motion is not chaotic [27]. It could also be seen from the Figure 5 that all three approaches to calculation of reaction forces gave the same result. The relative error of the center of mass location calculation with during steadystate regime was 0.4% and 0.3% for the linearization approach and the neural network approach, respectively. With the relative error remains small but accuracy differs by 10 times: 0.3% and 3% for the linearization approach and the neural network approach, respectively. Preservation of accuracy of the neural network approach could be explained by the fact that the model of training was carried out with value of the dimensionless centrifugal force close to 1.6. Finally, with the error was no more than 2% for the neural network approach and no more than 6% for the linearization approach. So, the neural network approach compared to the linearization approach given that the amplitude of oscillations is more than has been proven to be significantly more accurate.
Further increase of the imbalance was not studied as it usually results in unstable trajectories. Increase of imbalance leads to increase of vibration’s amplitude. That is why the linearization method could give inadequate results due to actual nonlinearity of the reaction forces at large eccentricities. The neural network approach gives good interpolation results, that is, within the learning sample. During extrapolation, deviation from the learning conditions could lead to unpredictable results and rapid increase of the error.
3.2. Rotor Dynamics for Various Speed Regimes and Transition from Low to High Vibrations
The same rotor system as in Section 3.1 of the present paper has been studied. Imbalance was constant kg·m, and tests were run for various rotation speeds: 1000 rpm, 3000 rpm, 5000 rpm, and 7000 rpm.
The first series of simulations of rotor’s behavior at various rotation speeds was implemented using the trajectory approach. Calculated motion time was set to 1 sec. Apart from trajectories, the location of equilibrium points has been determined, marked “+” in Figure 6(a). The second series of simulations has been implemented using the linearization approach. For that, matrices of stiffness and damping have been determined at four points on the equilibrium locus curve.
(a)
(b)
Then, four neural networks have been developed for training with neurons in the hidden layer. For every network, matrices and of and elements were used as input and output data, respectively. The matrices of input and output data were divided into three groups for training, validation and testing in fractions 0.6, 0.3, and 0.1, respectively. Training of each network was done in 1000 epochs in approximately 1–3 min. Training quality is characterized by the following values of mean squared errors and correlation coefficients: N, .
The calculation results showed that in the range of speeds from 1000 to 5000 rpm the vibrations of the rotor have the frequency equal to the rotation frequency (see Figure 6(b)). In this range of speed the oil whirl does not appear under given conditions, and the main reason of the occurrence of vibrations is the imbalance of the rotor. Additional calculations have shown that the critical speed of the rotor is 3300 rpm. In the range of speed from 6000 to 7000 rpm some increased vibrations occur that correspond to approximately twice critical speed (see Figures 6(a) and 6(b)). The vibrations of the rotor have the frequency close to the critical frequency (see Figure 6(b)). This is typical of oil whip [6]. Here, linearization approach becomes unstable, and the neural network approach, on the other hand, remains stable with sufficient accuracy (see Table 1).

So, the neural network approach allows modeling of significantly nonlinear vibrations, such as oil whip. The necessity of development of a separate neural network for every rotation speed; however, it is a drawback of the applied approach. Consequently, in the next paragraph a neural network is considered with an additional input value, rotation speed of the rotor.
3.3. Rotor Dynamics during Transient Regimes
Modeling of rotor dynamics during transient start/stop regime is more complex and costly in terms of calculation resources. The main advantages of the linearization approach that are essentially simple implementation and good computational performance manifest themselves weakly in this case. This is due to necessity to implement a number of intermediate calculations of equilibrium position, of stiffness and damping matrices, and due to interpolation procedures. In case of the neural network approach, such additional interpolation procedures are not needed. There is also a possibility of training a neural network based on one experiment and not on a series of them. This is why application of the linearization approach has not been considered in this section.
The transient regime has been studied using the test rig, shown in the Figure 7. The main parameters of the test rig are the same as in the Section 3.2. During the experiment the rotor was steadily accelerated up to 4000 rpm and then the electromotor was shut down and disconnected. During the deceleration the data from tachometer and acceleration sensors was acquired. Deceleration of the rotor took about 1 min. Based on the results of rotation speed measurements, rotor’s speed was approximated as a linear function of time, which describes the deceleration of the rotor during hydrodynamic lubrication regime (down to 500 rpm) accurately enough. Additional experiments at constant rotation speed in the 1000–4000 rpm range showed that in this range the rotor vibrates at frequencies equal to the rotation frequency. This result matches the results in the Section 3.2.
In this case, one neural network has been developed for training with neurons in the hidden layer. Additionally, the rotation speed was taken as an input, so the matrix of inputs was taken as . The combined sample of trajectories of the rotor at eight different rotation speeds rpm was used for training. In other words, the training function had a steplike form. The matrices of input and output data were divided into three groups for training, validation, and testing in fractions 0.7, 0.15, and 0.15, respectively. Training took no more than 10 minutes. Training quality is characterized by the following values of mean squared errors and correlation coefficients: N; .
It has to be noted that one of the alternative ways to train a neural network during monotonous smooth change of speed was not successful. The results obtained based on such neural networks were unstable. This could be explained by fact that every value of speed occurs in a monotonous function only once. In a steplike function, many intermediate values are absent, but the neural network manages to successfully interpolate them.
In Figures 8(a) and 8(b) the results of deceleration of the rotor simulation are shown for trajectory and neural networks approaches. In Figure 8(c) the comparison of numerical and experimental results is shown. It could be seen from the graphs that numerical results match well quantitatively (see Figures 8(a) and 8(b)) and qualitatively the experimental results (see Figure 8(c)). The calculated value of the critical speed matches the experimentally obtained value quite well (see Figure 8(c)).
(a)
(b)
(c)
So, the solution of the rotor dynamics problems during transition regimes could be successfully implemented using neural networks.
4. Conclusion
An approach to neural network program module formation has been presented aimed at calculation of the reaction force’s components of the fluid film depending on the position and oscillation velocities of the center of a rotor. Three approaches used in solving the rotor dynamics problem have been compared: the trajectory approach, the linearization approach, and the neural network approach to approximation of the oil film reaction forces. As a result of quantitative analysis, it has been determined that the use of linear and neural approximation allows decreasing the calculation time up to two orders of magnitude. It has also been determined that if the amplitude of the rotor’s oscillation exceeds 10% of the mean radial gap, the neural network approach is up to 3 times more accurate compared to the linearization approach. Also, neural networks were used to simulate significantly nonlinear regimes of rotor’s vibration with high vibrodisplacement and demonstrated good results in simulation of transient regimes. The considered models of a rotor and the fluid film were relatively simple, which significantly improved performance using conventional methods and aided sampling for neural network’s training. Further complication of the models, as well as training under condition of interaction with a real rotor machine, shall not significantly influence the main characteristics of artificial neural networks.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Authors’ Contributions
Alexey Kornaev developed the fluid flow and rotor’s dynamic mathematical model, scripted the simulation program, and wrote the manuscript. Nickolay Kornaev developed the artificial neural network program module and performed the numerical experiment. Elena Kornaeva and Leonid Savin analyzed the experiment’s results. Leonid Savin was in charge of supervising this work.
Acknowledgments
This work was supported by the Russian Science Foundation under Project no. 161900186 (formulation and solution of the problem in all sections of the present paper) and by the Ministry of Education and Science of the Russian Federation under Project no. 9.2952.2017/4.6 (development and modernization of the test rig in Section 3.3). The authors gratefully acknowledge this support. The authors also express gratitude to Alexander Fetisov for carrying out the experiment and to Alexander Babin for translation.
References
 D. Childs, Turbomachinery Rotordynamics. Phenomena, Modeling, and Analysis, John Wiley and Sons, New York, NY, USA, 1993.
 M. G. Vance, Rotordynamics of Turbomachinery, John Wiley and Sons, New York, NY, USA, 1988.
 T. Yamamoto, Linear and Nonlinear Rotordynamics. A Modern Treatment with Applications, John Wiley and Sons, New York, NY, USA, 2001.
 L. Savin, O. Solomin, and D. Ustinov, “Rotor dynamics on friction bearings with cryogenic lubrication,” in Proceedings of the Tenth World Congress on the Machines and Mechanisms, vol. 4, pp. 1716–1721, June 1999. View at: Google Scholar
 S. P. Asok, K. Sankaranarayanasamy, T. Sundararajan, K. Rajesh, and G. Sankar Ganeshan, “Neural network and CFDbased optimisation of square cavity and curved cavity static labyrinth seals,” Tribology International, vol. 40, no. 7, pp. 1204–1216, 2007. View at: Publisher Site  Google Scholar
 Y. Hori, Hydrodynamic lubrication, Yokendo Ltd., Tokyo, 2006.
 M. Kalkat, S. Yildirim, and I. Uzmay, “Rotor dynamics analysis of rotating machine system using artificial neural networks,” International Journal of Rotating Machinery, p. 9, 2003. View at: Publisher Site  Google Scholar
 M. Karkoub and A. Elkamel, “Modelling pressure distribution in a rectangular gas bearing using neural networks,” Tribology International, vol. 30, no. 2, pp. 139–150, 1997. View at: Publisher Site  Google Scholar
 A. D. Dongare and A. D. Kachare, “Predictive tool: an artificial neural network,” International Journal of Engineering and Innovative Technology (IJEIT), vol. 2, pp. 209–2014, 2012. View at: Google Scholar
 P. Qin, Y. Shen, J. Zhu, and H. Xu, “Dynamic analysis of hydrodynamic bearingrotor system based on neural network,” International Journal of Engineering Science, vol. 43, no. 56, pp. 520–531, 2005. View at: Publisher Site  Google Scholar
 A. O. Kurban and Ş. Yildirim, “Analysis of a hydrodynamic thrust bearing with elastic deformation using a recurrent neural network,” Tribology International, vol. 36, no. 12, pp. 943–948, 2003. View at: Publisher Site  Google Scholar
 J. Ghorbanian, M. Ahmadi, and R. Soltani, “Design predictive tool and optimization of journal bearing using neural network model and multiobjective genetic algorithm,” Scientia Iranica, vol. 18, no. 5, pp. 1095–1105, 2011. View at: Publisher Site  Google Scholar
 M. Simsir, R. Bayir, and Y. Uyaroglu, “Realtime monitoring and fault diagnosis of a low power hub motor using feedforward neural network,” Computational Intelligence and Neuroscience, vol. 2016, Article ID 7129376, 2016. View at: Publisher Site  Google Scholar
 B. Samanta, K. R. AlBalushi, and S. A. AlAraimi, “Bearing fault detection using artificial neural networks and genetic algorithm,” Eurasip Journal on Applied Signal Processing, vol. 2004, no. 3, pp. 366–377, 2004. View at: Publisher Site  Google Scholar
 M. Subrahmanyam and C. Sujatha, “Using neural networks for the diagnosis of localized defects in ball bearings,” Tribology International, vol. 30, no. 10, pp. 739–752, 1997. View at: Publisher Site  Google Scholar
 K. M. Saridakis, P. G. Nikolakopoulos, C. A. Papadopoulos, and A. J. Dentsoras, “Identification of wear and misalignment on journal bearings using artificial neural networks,” Proceedings of the Institution of Mechanical Engineers, Part J: Journal of Engineering Tribology, vol. 226, no. 1, pp. 46–56, 2012. View at: Publisher Site  Google Scholar
 S. G. Tores, C. Bonello, and P. Bonello, “Empirical identification of the inverse model of a squeezefilm damper bearing using neural networks and its application to a nonlinear inverse problem,” Journal of Vibration and Control, p. 22, 2016. View at: Google Scholar
 S. Morishita and T. Ura, “ER fluid applications to vibration control devices and an adaptive neuralnet controller,” Journal of Intelligent Material Systems and Structures, vol. 4, no. 3, pp. 366–372, 1993. View at: Publisher Site  Google Scholar
 R. Jha and C. h. He, “Adaptive neurocontrollers for vibration suppression of nonlinear and time varying structures,” Journal of Intelligent Material Systems and Structures, vol. 15, p. 11, 2004. View at: Google Scholar
 S. Y. Chen, Y. C. Hung, Y. H. Hung, and C. H. Wu, “Application of a recurrent wavelet fuzzyneural network in the control of a magneticbearing mechanism,” Computers Electrical Engineering, vol. 54, pp. 147–158, 2016. View at: Publisher Site  Google Scholar
 P.Y. Couzon and J. der Hagopian, “Neurofuzzy active control of rotor suspended on active magnetic bearing,” Journal of Vibration and Control, vol. 13, no. 4, pp. 365–384, 2007. View at: Publisher Site  Google Scholar
 C. R. Burrows, P. S. Keogh, and M. N. Sahinkaya, “Progress towards smart rotating machinery through the use of active bearings,” Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science, vol. 223, p. 11, 2009. View at: Publisher Site  Google Scholar
 Q. J. Wang and Y. Chung, Encyclopedia of Tribology, Springer, Boston, Mass, USA, 2013. View at: Publisher Site
 A. Kornaev, L. Savin, E. Kornaeva, and A. Fetisov, “Influence of the ultrafine oil additives on friction and vibration in journal bearings,” Tribology International, vol. 101, pp. 131–140, 2016. View at: Publisher Site  Google Scholar
 S. Haykin, Neural Networks. A Comprehensive Foundation, Prentice Hall, NJ, USA, 2nd edition, 1999.
 G. A. Korn and T. M. Korn, Mathematical Handbook for Scientists and Engineers, Dover Publications, New York, NY, USA, 2000. View at: MathSciNet
 F. C. Moon, Chaotic Vibrations. An Introduction for Applied Scientists and Engineers, John Wiley and Sons, New York, NY, USA, 1987.
 L. Cui and J. Zheng, “Nonlinear vibration and stability analysis of a flexible rotor supported on angular contact ball bearings,” JVC/Journal of Vibration and Control, vol. 20, no. 12, pp. 1767–1782, 2014. View at: Publisher Site  Google Scholar
 Y. Wang and X. Wang, “Nonlinear vibration analysis for a Jeffcott rotor with seal and airfilm bearing excitations,” Mathematical Problems in Engineering, vol. 2010, Article ID 657361, 15 pages, 2010. View at: Publisher Site  Google Scholar
 M. G. Salvadori and M. L. Baron, Numerical Methods in Engineering, Englewood Cliffs, NJ, USA, 2nd edition, 1961.
 http://www.gnu.org/software/octave.
 http://www.mathworks.com/help/nnet/ref/feedforwardnet.html.
 A. Harnoy, Bearing Design in Machinery: Engineering Tribology and Lubrication, Marcel Dekker Inc, New York, NY, USA, 2003.
 A. Chasalevris, “An investigation on the dynamics of highspeed systems using nonlinear analytical floating ring bearing models,” International Journal of Rotating Machinery, vol. 2016, Article ID 7817134, 22 pages, 2016. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2017 Alexey V. Kornaev et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.