Research Article  Open Access
Weitao Jia, Feng Gao, Yan Li, Wenwu Wu, Zhongwei Li, "Nonlinear Dynamic Analysis and Chaos Prediction of Grinding Motorized Spindle System", Shock and Vibration, vol. 2019, Article ID 9643124, 10 pages, 2019. https://doi.org/10.1155/2019/9643124
Nonlinear Dynamic Analysis and Chaos Prediction of Grinding Motorized Spindle System
Abstract
The paper determines the impact factors of dynamics of a motorized spindle rotor system due to high speed: centrifugal force and bearing stiffness softening. A nonlinear dynamic model of the grinding motorized spindle system considering the above impact factors is constructed. Through system simulation including phase portraits and Poincaré map, the periodic behavior and chaotic behavior of the nonlinear grinding motorized spindle system are revealed. The threshold curve of chaos motion is obtained through the Melnikov method. The conclusion can provide a theoretical basis for researching deeply the dynamic behaviors of the grinding motorized spindle system.
1. Introduction
Although grinding motorized spindle machining technology has been widely used in ultraprecision and highspeed machine tools that can realize highspeed ultraprecision machining at a low cost, improvement of stability is still an important subject [1]. The unbalance of rotor mass is the main factor affecting the dynamic behavior of the motorized spindle system. The instability of the system results in abnormal operation, serious damage, and even “axle holding” [2, 3]. Therefore, analyzing the stability and predicting the threshold curve of chaos motion of motorized spindle system are necessary.
In engineering practice, there is still the unbalance of rotor mass. Many scholars have conducted indepth research on this. The studies showed that dynamic characteristics of the system are related to the unbalance of rotor mass [4–7]. Huihui and Shuyun [8, 9], taking water lubricated motorized spindle as the research object, analyzed the dynamic behavior of the motorized spindle system under the coupling action of the unbalance of rotor mass and bearing tilt effect are studied. Bo et al. [10] investigated the dynamic performance of motorized spindle system under centrifugal force and bearing stiffness softening. The results display that the dynamic behavior of the system is affected seriously by the two factors and must be taken into account when modeling the motorized spindle system.
Under the highspeed running state, the rigidity of the motorized spindle bearing is softened to some extent. The centrifugal force of the spindle is the main influencing factor of the stiffness softening of the motorized spindle in the existing literature [11–15].
The dynamic model established from practical problem is generally nonlinear, and the system has complex dynamic behavior. Chaos is one of the unique phenomena of nonlinear systems. Along with the progress of nonlinear dynamics, nonlinear behaviors of the system such as stability solution, periodic solution, bifurcation, and chaos, have become hot internet topic. Many scholars have focused on rotor bearing system models, and the main methods for analyzing the dynamic behavior of nonlinear systems include axis orbit, power spectra, phase portraits, Poincare map, and bifurcation diagrams. The modeling results can predict the stability parameters of the rotor bearing system and avoid the range of unstable motion of the rotor bearing system [16–20].
Li et al. [21] & Sheu et al. [22] investigated the nonlinear dynamics of Duffing system. The nonlinear dynamic characteristics of the Duffing system are researched by timedomain diagram, phase portraits, Poincaré map, and bifurcation diagrams.
At present, numerical analysis is mainly used to study nonlinear dynamics, but it cannot provide the analytic solution of the nonlinear system. Melnikov method proposed a new and more applicable method of chaos criterion, and proved its effectiveness by solving nonlinear system described through the Duffing equation. Because of the effectiveness of Melnikov method, many scholars are still applying it nowadays. Yagasaki [23] & Zhang et al. [24] utilized the Melnikov method, and the research show that abundant bifurcation behavior and chaos motions occur in vibrations of the microcantilever. Yu et al. [25, 26] investigated nonlinear dynamics of a flexible rotor system, and obtained the parameter range of the occurrence of chaos motion in the sense of Smale horseshoes. Siewe et al. [27, 28] et al. obtained the condition of heteroclinic and homoclinic chaos of Rayleigh and Van der Pol oscillator by applying Melnikov method. Zhang and Zhou [29] investigated the chaotic behaviors of a nanoplate by the Melnikov method. The threshold curve of chaotic motion is given. Jibin and Fengjuan [30], Liu et al. [31], and Li et al. [32] gave the critical parameter curve of the Duffing oscillator for the existence of chaos in the Smale horse sense by Melnikov method.
At high speed, bearing softening will occur, which strongly affects the dynamic behavior of motorized spindle system. However, in previous literatures, the study of bearing softening effect is relatively rare, mainly considering the system as a hard spring system. The novelty of this paper is to consider the effect of bearing softening on the motorized spindle system.
This article established the nonlinear dynamic model of the nonlinearly supported motorized spindle rotor bearing system under centrifugal force and bearing stiffness softening. The nonlinear characteristics of motorized spindle system are shown using phase portraits and Poincaré map by means of numerical analysis. The threshold curve of motorized spindle system is obtained under periodic perturbation through Melnikov method. The results could be used to design and optimize motorized spindle system, and it is of practical significance for the purpose of guaranteeing system stability.
2. Mathematical Formulation
The stiffness and damping of motorized spindle system may be influenced such as bearing, pedestal, and so on. The model of motorized spindle system is as follows:
The rotor bearing model with nonlinear support is established by equivalent spring and damping. Then, the nonlinear dynamics equation of motorized spindle system under centrifugal force is as follows:where m, e, and ω are mass of the rotor, mass eccentricity, and rotating speed, respectively. c, k_{1}, and k_{3} denote damping, linear, and nonlinear stiffness of bearing.
Since there is no coupling term between the first and the second formulas in equation (1), the vibration of each degree of freedom is relatively independent, and the vibration in x_{2} direction is the same as that in x_{1}, but there is only a phase difference, so in this paper, only one of the equations needs to be investigated.
The highspeed grinding motorized spindle system is a kind of electromechanical coupling system. The rotor bearing system exhibits soft characteristics; therefore, a < 0.
The first of equation (1) is rewritten aswhere ξ = c/m and .
One can rewrite the equation of motion (2) in the dimensionless form as [30]where the dimensionless parameters are , , , , , and .
For the rotor bearing system of highspeed grinding motorized spindle, r << 1, f << 1, so r, f is regarded as a small parameter. Equation (3) is the Duffing equation with weak periodic forced oscillation equation.
3. Simulation Analysis
Phase portraits and Poincaré map are obtained by using displacement and velocity of system response as horizontal and longitudinal coordinates. That is to say horizontal coordinates is displacement, and longitudinal coordinates is velocity. Poincare diagram is obtained by sampling every rotation cycle. If the response of the system is periodic n, there are n isolated points on Poincaré map. If the system is periodic motion, it presents a closed curve. If the system is chaotic motion, it is a disordered point.
The Runge–Kutta method is used to solve the differential equation of the system. Taking 150MD24Z7.5 highspeed grinding motorized spindle for example, the motorized spindle system is found to display intricate dynamics characteristics.
3.1. Response of the System at Different Rotational Speeds for e = 0.001 mm, k_{1} = 2.2488 × 10^{9} N/m
Phase portraits and Poincaré map of the system are shown in Figures 2–6, respectively, ω as the control parameter.
(a)
(b)
(a)
(b)
(a)
(b)
(a)
(b)
(a)
(b)
When the speed is 3000 r/min, it can be seen from Figure 2(a) that the phase portraits form a closed curve, and there are three isolated points in the Poincaré map; then the system is period3(P3) motion.
In the same way, when the speed reaches 5000 r/min, there are eleven isolated points in Poincaré map the system is period11(P11) motion.
When the speed is 10000 r/min, phase portraits form a complex curve, and there are a lot of disordered points in the Poincaré map. Thus, the system is in chaotic motion.
In the same way, when the speed is 15000 r/min, the system is still in chaotic motion.
As the rotational speed continues to increase, when the speed reaches 20000 r/min, and there are two points in Poincaré map, the system returns to period2(P2) motion.
It can be seen from Figures 2–6 that when the rotation speed is 3000 r/min, 5000 r/min, or 20000 r/min, phase portraits form a closed figure. When the rotating speed is 10000 r/min or 15000 r/min, phase portraits are irregular curves. From the Poincaré map, it takes the general form (P3 ⟶ P11 ⟶ chaos ⟶ P2) with the change of rotational speed between the values (3000 r/min∼20000 r/min).
3.2. Response of the System at Different Rotational Speeds for e = 0.001 mm, k_{1} = 2.2488 × 10^{10} N/m
Phase portraits and Poincaré map of the system are shown in Figures 7–11, respectively, with ω as the control parameter.
(a)
(b)
(a)
(b)
(a)
(b)
(a)
(b)
(a)
(b)
When the speed is 3000 r/min, it can be seen from Figure 7(a) that the phase portraits is a closed curve, and there is one point in the Poincaré map; then the dynamic response undergoes synchronous vibration with period1(P1) motion.
In the same way, when the speed is 5000 r/min, the system is still in period1(P1) motion.
When the speed increases to 10000 r/min, as can be seen from Figure 7, and there are seven points in Poincaré map, the system is in period7(P7) motion.
Similarly, when the speed is 15000 r/min, the system is still in period3(P3) motion.
It is obvious that when the speed is 20000 r/min, there are countable points, so the system is in periodn(Pn) motion.
It can be seen from Figures 7–11, no matter what the speed is, phase portraits form a closed figure. Thus, the system has gone through a process (P1 ⟶ P7 ⟶ P3 ⟶ Pn) from the Poincaré map when the rotating speed varied from 3000 r/min to 20000 r/min.
In short, the nonlinear characteristics of the motorized spindle system are shown in Table 1 synthetically.

In Table 1, when k_{1} = 2.2488 × 10^{9}, the nonlinear characteristics of motorized spindle system have undergone different stages as the rotating speed varied, from periodic to chaotic, and finally back to periodic; when k_{1} = 2.2488 × 10^{10}, the system is periodic motion. As shown in Figures 2–11, with the increase of stiffness, it can be seen that the stiffness suppresses the occurrence of chaos, improving the vibration behavior for rotor bearing system of grinding motorized spindle.
4. Chaotic Prediction
The research and application of chaos theory has been extended to many fields such as natural science and social science and has become an important part o f modern nonlinear system theory. With the deepening of scientific research, the status of prediction is particularly prominent. In this paper, the Melnikov method is used to predict the threshold of the occurrence of chaos motion [27, 28].
Equation (3) is the Duffing equation with weak periodic forced oscillation equation. Then, the chaotic threshold of rotor bearing system of grinding motorized spindle is studied based on the Melnikov method.
The systems shown in equation (3) are equivalent to
Equation (4) is written in matrix form:where , , and .
When ε = 0, the unperturbed system is
The system of (6) is a Hamilton system with a potential function
Obviously, this is a singlewell system. The potential energy surface of the system in the threedimensional phase space and the potential energy function on the plane are shown in Figure 12.
The Hamilton function of the unperturbed is given by
The Hamilton function of the unperturbed is shown in Figure 13.
The characteristic equation of the linear approximation system = 0, its characteristic root is = ±i; therefore, (0, 0) is the center, while (1, 0) and (−1, 0) are singularity. The characteristic equation is:
So = ±, (1, 0) and (−1, 0) are saddle points. When h = 1/4, there are two heteroclinic trajectories connected to (1, 0) and (−1, 0), and these two trajectories make up a heteroclinic circle. When 0 < h < 1/4, there is a closed loop surrounded by a family of (0, 0), as shown in Figure 14.
By integrating, the parameter equations of two heteroclinic trajectories are obtained:
The corresponding Melnikov functions are as follows:where
Therefore,where I_{1} = 2/3 and I_{2} = .
Because = 0,
If the equation is solved, then
Because = , then has its single zero point; thus, the system is chaotic. Therefore, when the ratio of the parameter f to the r exceeds the threshold , the system has chaos in the sense of Smale horseshoe.
The chaotic threshold of the rotor bearing system of grinding motorized spindle is shown in Figure 15.
5. Conclusion
In this article, a nonlinear dynamic model of grinding motorized spindle system considering centrifugal force and bearing stiffness softening has been proposed to research dynamical behavior and predict the threshold of chaos motion of system. Several results are as follows:(1)Through the system simulation including phase portraits and Poincaré map, the periodic and chaotic behavior of grinding motorized spindle system are revealed as rotating speed increases.(2)Comparing the dynamic behavior of the motorized spindle system under different stiffnesses, it is found that the system exhibits good nonlinearity as stiffness increases. At the same time, it is shown that the increase of system stiffness can suppress the occurrence of chaos.(3)Using the Melnikov method, the threshold curve of chaos motion has been predicted. The results could be used to provide the theory basis for practical production.
The results could help engineers and researchers design and research the grinding motorized spindle system and provide a theoretical basis for further research about complex nonlinear characteristics, selecting the range of system stable operating parameters.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
The authors acknowledge support by the Shaanxi Provincial Key Research and Development Program (grant no. 2018ZDXMGY074) and Natural Science Basic Research Plan in Shaanxi Province of China (program no. 2019JM148).
References
 X. Wanli, Y. Xuebing, L. Lang, and Y. Julong, “Review on key technology of hydrodynamic and hydrostatic highfrequency motor spindles,” Journal of Mechanical Engineering, vol. 45, no. 9, pp. 1–18, 2009. View at: Publisher Site  Google Scholar
 X. Wanli, H. Zhiquan, L. Lang, and Y. Xuebing, “Review on the dynamic characteristics of aerostatic motorized spindles,” Journal of Mechanical Engineering, vol. 47, no. 5, pp. 40–58, 2011. View at: Publisher Site  Google Scholar
 J. S. Larsen and I. F. Santos, “On the nonlinear steadystate response of rigid rotors supported by air foil bearingstheory and experiments,” Journal of Sound and Vibration, vol. 346, pp. 284–297, 2015. View at: Publisher Site  Google Scholar
 L. Shusen, “Analysis of the dynamics of precision centrifuge spindle system with the externally pressurized gas bearing,” Chinese Journal of Mechanical Engineering, vol. 41, no. 2, pp. 28–32, 2005. View at: Publisher Site  Google Scholar
 C. Dongju, F. Jinwei, and Z. Feihu, “Dynamic and static characteristics of a hydrostatic spindle for machine tools,” Journal of Manufacturing Systems, vol. 31, no. 1, pp. 26–33, 2012. View at: Publisher Site  Google Scholar
 H. M. Navazi and M. Hojjati, “Nonlinear vibrations and stability analysis of a rotor on highstaticlowdynamicstiffness supports using method of multiple scales,” Aerospace Science and Technology, vol. 63, pp. 259–265, 2017. View at: Publisher Site  Google Scholar
 M. R. Rashidi, M. A. Karami, and B.N. Firooz, “Numerical analysis of a rigid rotor supported by aerodynamic fourlobe journal bearing system with mass unbalance,” Communications in Nonlinear Science & Numerical Simulation, vol. 17, no. 1, pp. 454–471, 2012. View at: Publisher Site  Google Scholar
 F. Huihui and J. Shuyun, “Dynamic analysis of waterlubricated motorized spindle considering tilting effect of thrust bearing,” Proceedings of the Institution of Mechanical Engineers, Part C Journal of Mechanical Engineering Science, vol. 231, no. 20, pp. 3780–3790, 2016. View at: Publisher Site  Google Scholar
 F. Huihui and J. Shuyun, “Dynamics of a motorized spindle supported on waterlubricated bearings,” Proceedings of the Institution of Mechanical Engineers, Part C Journal of Mechanical Engineering Science, vol. 231, no. 3, pp. 459–472, 2015. View at: Publisher Site  Google Scholar
 W. Bo, S. Wei, and W. BangChun, “The effect of high speeds on dynamic characteristics of motorized spindle system,” Engineering Mechanics, vol. 32, no. 6, pp. 231–237, 2015. View at: Google Scholar
 J.S. Chen and Y.W. Hwang, “Centrifugal force induced dynamics of a motorized highspeed spindle,” The International Journal of Advanced Manufacturing Technology, vol. 30, no. 12, pp. 10–19, 2006. View at: Publisher Site  Google Scholar
 C.W. Lin, J. F. Tu, and J. Kamman, “An integrated thermomechanicaldynamic model to characterize motorized machine tool spindles during very high speed rotation,” International Journal of Machine Tools & Manufacture, vol. 43, no. 10, pp. 1035–1050, 2003. View at: Publisher Site  Google Scholar
 J.J. Sinou, “Nonlinear dynamics and contacts of an unbalanced flexible rotor supported on ball bearings,” Mechanism and Machine Theory, vol. 44, no. 9, pp. 1713–1732, 2009. View at: Publisher Site  Google Scholar
 C. Li and Z. Jianrong, “Nonlinear vibration and stability analysis of a flexible rotor supported on angular contact ball bearings,” Journal of Vibration and Control, vol. 20, no. 12, pp. 1767–1782, 2014. View at: Publisher Site  Google Scholar
 C. Xing, Y. Shihua, and P. Zengxiong, “Nonlinear vibration for PMSM used in HEV considering mechanical and magnetic coupling effects,” Nonlinear Dynamics, vol. 80, no. 12, pp. 541–552, 2015. View at: Publisher Site  Google Scholar
 C.W. ChangJian, “Nonlinear dynamic analysis of dual flexible rotors supported by long journal bearings,” Mechanism and Machine Theory, vol. 45, no. 6, pp. 844–866, 2010. View at: Publisher Site  Google Scholar
 C.W. ChangJian and C.K. Chen, “Bifurcation and chaos analysis of a flexible rotor supported by turbulent long journal bearings,” Chaos, Solitons & Fractals, vol. 34, no. 4, pp. 1160–1179, 2007. View at: Publisher Site  Google Scholar
 C.W. ChangJian and C.K. Chen, “Chaos and bifurcation of a flexible rotor supported by porous squeeze couple stress fluid film journal bearings with nonlinear suspension,” Chaos, Solitons & Fractals, vol. 35, no. 2, pp. 358–375, 2008. View at: Publisher Site  Google Scholar
 C.W. ChangJian and C. O.K. Chen, “Chaotic response and bifurcation analysis of a flexible rotor supported by porous and nonporous bearings with nonlinear suspension,” Nonlinear Analysis: Real World Applications, vol. 10, no. 2, pp. 1114–1138, 2009. View at: Publisher Site  Google Scholar
 L. Xiang, A. Hu, L. Hou, Y. Xiong, and J. Xing, “Nonlinear coupled dynamics of an asymmetric doubledisc rotorbearing system under rubimpact and oilfilm forces,” Applied Mathematical Modelling, vol. 40, no. 78, pp. 4505–4523, 2016. View at: Publisher Site  Google Scholar
 Z. Li, D. Chen, J. Zhu, and Y. Liu, “Nonlinear dynamics of fractional order Duffing system,” Chaos, Solitons & Fractals, vol. 81, pp. 111–116, 2015. View at: Publisher Site  Google Scholar
 L.J. Sheu, H.K. Chen, J.H. Chen, and L.M. Tam, “Chaotic dynamics of the fractionally damped Duffing equation,” Chaos, Solitons & Fractals, vol. 32, no. 4, pp. 1459–1468, 2007. View at: Publisher Site  Google Scholar
 K. Yagasaki, “Bifurcations and chaos in vibrating microcantilevers of tapping mode atomic force microscopy,” International Journal of Nonlinear Mechanics, vol. 42, no. 4, pp. 658–672, 2007. View at: Publisher Site  Google Scholar
 W. Zhang, M. H. Yao, and J. H. Zhang, “Using the extended Melnikov method to study the multipulse global bifurcations and chaos of a cantilever beam,” Journal of Sound and Vibration, vol. 319, no. 12, pp. 541–569, 2009. View at: Publisher Site  Google Scholar
 T.J. Yu, S. Zhou, X.D. Yang, and W. Zhang, “Global dynamics of a flexible asymmetrical rotor,” Nonlinear Dynamics, vol. 91, no. 2, pp. 1041–1060, 2018. View at: Publisher Site  Google Scholar
 T.J. Yu, W. Zhang, and X.D. Yang, “Global bifurcations and chaotic motions of a flexible multibeam structure,” International Journal of Nonlinear Mechanics, vol. 95, pp. 264–271, 2017. View at: Publisher Site  Google Scholar
 M. S. Siewe, H. Cao, and M. A. F. Sanjuán, “On the occurrence of chaos in a parametrically driven extended Rayleigh oscillator with threewell potential,” Chaos, Solitons & Fractals, vol. 41, no. 2, pp. 772–782, 2009. View at: Publisher Site  Google Scholar
 M. S. Siewe, F. M. M. Kakmeni, C. Tchawoua, and P. Woafo, “Bifurcations and chaos in the triplewell Van der Pol oscillator driven by external and parametric excitations,” Physica A: Statistical Mechanics and its Applications, vol. 357, no. 34, pp. 383–396, 2005. View at: Publisher Site  Google Scholar
 X. Zhang and L. Zhou, “Melnikov’s method for chaos of the nanoplate postulating nonlinear foundation,” Applied Mathematical Modelling, vol. 61, pp. 744–749, 2018. View at: Publisher Site  Google Scholar
 L. Jibin and C. Fengjuan, Chaos, Melnikov Method and New Development, Science Press, Beijing, China, 2012.
 Z. R. Liu, W. G. Yao, and Z. X. Zhu, “Road to chaos for a soft spring system under weak periodic disturbance,” Applied Mathematics & Mechanics, vol. 7, no. 2, pp. 111–116, 1986. View at: Publisher Site  Google Scholar
 H. Li, X. Liao, S. Ullah, and L. Xiao, “Analytical proof on the existence of chaos in a generalized Duffingtype oscillator with fractionalorder deflection,” Nonlinear Analysis: Real World Applications, vol. 13, no. 6, pp. 2724–2733, 2012. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2019 Weitao Jia 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.