# Nonlinear Dynamic Analysis and Chaos Prediction of Grinding Motorized Spindle System

Authors: Weitao Jia; Feng Gao; Yan Li; Wenwu Wu; Zhongwei Li
Journal: Shock and Vibration
 (2019)

---

## 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.


---

## Body

## 1. Introduction

Although grinding motorized spindle machining technology has been widely used in ultraprecision and high-speed machine tools that can realize high-speed 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 in-depth 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 high-speed 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 time-domain 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:Figure 1
Diagram of the motorized spindle system.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:(1)mx¨1+cx˙1+k1x1+k3x13=meω2cosωt,mx¨2+cx˙2+k1x2+k3x23=meω2sinωt,where m, e, and ω are mass of the rotor, mass eccentricity, and rotating speed, respectively. c, k1, and k3 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 x2 direction is the same as that in x1, but there is only a phase difference, so in this paper, only one of the equations needs to be investigated.The high-speed 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 as(2)x¨1+ξx˙1+ωn2x1−a^x13=eω2cosωt,where ξ = c/m and a^=k3/m.One can rewrite the equation of motion (2) in the dimensionless form as [30](3)y¨+ry˙+y−y3=fcosλτ,where the dimensionless parameters are y=kx1, k=a^/ωn2, τ=ωnt, r=ξ/ωn, f=eω2a^/ωn3, and λ=ω/ωn.For the rotor bearing system of high-speed 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 periodicn, 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 high-speed 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 fore = 0.001 mm, k1 = 2.2488 × 109 N/m

Phase portraits and Poincaré map of the system are shown in Figures2–6, respectively, ω as the control parameter.Figure 2
System response forω = 3000 r/min. (a) Phase portraits. (b) Poincaré map.
(a)
(b)Figure 3
System response forω = 5000 r/min. (a) Phase portraits. (b) Poincaré map.
(a)
(b)Figure 4
System response forω = 10000 r/min. (a) Phase portraits. (b) Poincaré map.
(a)
(b)Figure 5
System response forω = 15000 r/min. (a) Phase portraits. (b) Poincaré map.
(a)
(b)Figure 6
System response forω = 20000 r/min. (a) Phase portraits. (b) Poincaré map.
(a)
(b)When the speed is 3000 r/min, it can be seen from Figure2(a) that the phase portraits form a closed curve, and there are three isolated points in the Poincaré map; then the system is period-3(P3) motion.In the same way, when the speed reaches 5000 r/min, there are eleven isolated points in Poincaré map the system is period-11(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 period-2(P2) motion.It can be seen from Figures2–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 fore = 0.001 mm, k1 = 2.2488 × 1010 N/m

Phase portraits and Poincaré map of the system are shown in Figures7–11, respectively, with ω as the control parameter.Figure 7
System response forω = 3000 r/min. (a) Phase portraits. (b) Poincaré map.
(a)
(b)Figure 8
System response forω = 5000 r/min. (a) Phase portraits. (b) Poincaré map.
(a)
(b)Figure 9
System response forω = 10000 r/min. (a) Phase portraits. (b) Poincaré map.
(a)
(b)Figure 10
System response forω = 15000 r/min. (a) Phase portraits. (b) Poincaré map.
(a)
(b)Figure 11
System response forω = 20000 r/min. (a) Phase portraits. (b) Poincaré map.
(a)
(b)When the speed is 3000 r/min, it can be seen from Figure7(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 period-1(P1) motion.In the same way, when the speed is 5000 r/min, the system is still in period-1(P1) motion.When the speed increases to 10000 r/min, as can be seen from Figure7, and there are seven points in Poincaré map, the system is in period-7(P7) motion.Similarly, when the speed is 15000 r/min, the system is still in period-3(P3) motion.It is obvious that when the speed is 20000 r/min, there are countable points, so the system is in period-n(P-n) motion.It can be seen from Figures7–11, no matter what the speed is, phase portraits form a closed figure. Thus, the system has gone through a process (P1 ⟶ P7 ⟶ P3 ⟶ P-n) 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 Table1 synthetically.Table 1
Motion characteristics of the motorized spindle.
r/min
k
1 (N/m)
2.2488 × 109
2.2488 × 1010
3000
P3
P1
5000
P11
P1
10000
Chaos
P7
15000
Chaos
P3
20000
P2
P-nIn Table1, when k1 = 2.2488 × 109, 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 k1 = 2.2488 × 1010, 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.


## 3.1. Response of the System at Different Rotational Speeds fore = 0.001 mm, k1 = 2.2488 × 109 N/m

Phase portraits and Poincaré map of the system are shown in Figures2–6, respectively, ω as the control parameter.Figure 2
System response forω = 3000 r/min. (a) Phase portraits. (b) Poincaré map.
(a)
(b)Figure 3
System response forω = 5000 r/min. (a) Phase portraits. (b) Poincaré map.
(a)
(b)Figure 4
System response forω = 10000 r/min. (a) Phase portraits. (b) Poincaré map.
(a)
(b)Figure 5
System response forω = 15000 r/min. (a) Phase portraits. (b) Poincaré map.
(a)
(b)Figure 6
System response forω = 20000 r/min. (a) Phase portraits. (b) Poincaré map.
(a)
(b)When the speed is 3000 r/min, it can be seen from Figure2(a) that the phase portraits form a closed curve, and there are three isolated points in the Poincaré map; then the system is period-3(P3) motion.In the same way, when the speed reaches 5000 r/min, there are eleven isolated points in Poincaré map the system is period-11(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 period-2(P2) motion.It can be seen from Figures2–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 fore = 0.001 mm, k1 = 2.2488 × 1010 N/m

Phase portraits and Poincaré map of the system are shown in Figures7–11, respectively, with ω as the control parameter.Figure 7
System response forω = 3000 r/min. (a) Phase portraits. (b) Poincaré map.
(a)
(b)Figure 8
System response forω = 5000 r/min. (a) Phase portraits. (b) Poincaré map.
(a)
(b)Figure 9
System response forω = 10000 r/min. (a) Phase portraits. (b) Poincaré map.
(a)
(b)Figure 10
System response forω = 15000 r/min. (a) Phase portraits. (b) Poincaré map.
(a)
(b)Figure 11
System response forω = 20000 r/min. (a) Phase portraits. (b) Poincaré map.
(a)
(b)When the speed is 3000 r/min, it can be seen from Figure7(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 period-1(P1) motion.In the same way, when the speed is 5000 r/min, the system is still in period-1(P1) motion.When the speed increases to 10000 r/min, as can be seen from Figure7, and there are seven points in Poincaré map, the system is in period-7(P7) motion.Similarly, when the speed is 15000 r/min, the system is still in period-3(P3) motion.It is obvious that when the speed is 20000 r/min, there are countable points, so the system is in period-n(P-n) motion.It can be seen from Figures7–11, no matter what the speed is, phase portraits form a closed figure. Thus, the system has gone through a process (P1 ⟶ P7 ⟶ P3 ⟶ P-n) 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 Table1 synthetically.Table 1
Motion characteristics of the motorized spindle.
r/min
k
1 (N/m)
2.2488 × 109
2.2488 × 1010
3000
P3
P1
5000
P11
P1
10000
Chaos
P7
15000
Chaos
P3
20000
P2
P-nIn Table1, when k1 = 2.2488 × 109, 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 k1 = 2.2488 × 1010, 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(4)y˙=z,z˙=−y+y3+ε−ry˙+fcosλτ.Equation (4) is written in matrix form:(5)Y˙=fY+εgY,t,where Y=y,zT, fY=z,−y+y3T, and gY,t=0,−ry′+fcosλτT.Whenε = 0, the unperturbed system is(6)y˙=z,z˙=−y+y3.The system of (6) is a Hamilton system with a potential function(7)Ux=y22−y44.Obviously, this is a single-well system. The potential energy surface of the system in the three-dimensional phase space and the potential energy function on the plane are shown in Figure12.Figure 12
Potential function of the unperturbed system.The Hamilton function of the unperturbed is given by(8)Hy,z=12z2+12y2−14y2=h.The Hamilton function of the unperturbed is shown in Figure13.Figure 13
Hamilton function of the unperturbed system.The characteristic equation of the linear approximation systemλ′2+1 = 0, its characteristic root is λ1,2′ = ±i; therefore, (0, 0) is the center, while (1, 0) and (−1, 0) are singularity. The characteristic equation is:(9)0−λ′1−1+3y20−λ′=0.Soλ1,2′ = ±2, (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.Figure 14
Phase portraits of rotor bearing system of grinding motorized spindle.By integrating, the parameter equations of two heteroclinic trajectories are obtained:(10)y±0t=±th22t,z±0t=±22sech222t.The corresponding Melnikov functions are as follows:(11)M±t0=∫−∞+∞fqi0t∧gqi0t,t+t0dt,where(12)qi0t=y±0tz±0t,fqi0t=z±0t−y±0t+y±0t3,gqi0t,t+t0=0−rz±0t+fcosλt+t0,fqi0t∧gqi0t,t+t0=y±0t−y±0t3z±0t⋅0−rz±0t+fcosλt+t0=−rz±0t+fcosλt+t0z±0t.Therefore,(13)M±t0=∫−∞+∞−rz±0t+fcosλt+t0z±0tdt=rI1±fI2cosλt0,where I1 = 2/3 and I2 = πλcoshπλ/2.BecauseM±t0 = 0,(14)cosωt0=±rI1fI2.If the equation is solved, then(15)rI1fI2<1,that isfr>I1I2,BecauseM±′t0 = −λfI2sinλt0±⋅≠0, then M±t0 has its single zero point; thus, the system is chaotic. Therefore, when the ratio of the parameter f to the r exceeds the threshold 2sinhπλ/2/3πλ, 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 Figure15.Figure 15
Melnikov function threshold curve of rotor bearing system of motorized spindle.


## 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.



---

*Source: 9643124-2019-09-22.xml*