Nonlinear dynamic characteristic of the spindle-cutter system

Xiangsheng Gao1, Yidu Zhang2, Min Wang3 1, 3Beijing Key Laboratory of Advanced Manufacturing Technology, College of Mechanical Engineering and Applied Electronics Technology, Beijing University of Technology, Beijing, China 2State Key Laboratory of Virtual Reality Technology and Systems, Beihang University, Beijing, China 2Corresponding author E-mail: 1gaoxsh@bjut.edu.cn, 2ydzhang@buaa.edu.cn, 3wangm@bjut.edu.cn


Introduction
With the development of manufacturing industry, the machining speed has increased considerably over the past years and the spindle tool system is one of the key components in ensuring machining precision and accuracy.Sometimes the spindle system is in unusual working condition.Higher demands for nonlinear dynamic characteristic of spindle system are required.Therefore, research on nonlinear dynamic characteristic of spindle-cutter system is necessary.
Yu et al. [1] conducted nonlinear dynamic analysis of ultra-high speed grinding spindle system by using digital analysis methods, and dynamic characteristics of liquid hybrid bearing were analyzed by using Fluent software.Based on Hertz contact force formulation, Zhang et al. [2] established a six-degree-of-freedom model for investigating a machine-tool spindle system supported with ball bearings.The effects of negative bearing clearance on the dynamic response and stability of the machine tool spindle system with unbalanced force were investigated.Gao et al. [3] conducted unstable islands and bifurcation analysis in a spindle milling system, and dynamic milling forces were considered [3,4].Floquet theory was applied to investigate the stability of the periodic solution in [5].A non-uniform continuous Timoshenko beam model for spindle system with nonlinear and nonsmooth boundaries was established in [6].Galerkin method was used to discretize the governing nonlinear partial differential equations.Zhang et al. [7] considered effects caused by bearings and the interface of spindle-holder on the vibration response of spindle system.Si et al. [8] considered the nonlinear bearing force based on Hertz theory.Nonlinear bearing force caused by surface waviness on raceway was considered in [9].Nonlinear dynamic characteristic of the spindle system was studied via numerical simulation in [8,9], respectively.Lioulios et al. [10] studied the effects of rotational speed changes on dynamic characteristic of rotor-bearing system.In practical engineering applications, the vibration of cutter tip is always concerned.However, in the aforementioned researches, the vibration of cutter in a spindle system has received little consideration.Gao et al. [11] established a dynamic model of spindle-cutter coupling system including the gyroscopic effects and shear deformation.Based on Riccati transformation, the dynamic characteristic of the system was analyzed by using the whole transfer matrix method; however, the research focused on the vibration mode of the system and did not involve transient response of cutter tip.
In this paper, a nonlinear dynamic model of spindle-tool coupling system based on the second Newton's law is established considering milling cutter stiffness, cutting force, and clearance of bearing.According to the classic method of elastic theory, based on cross-sectional analysis, the cutter stiffness is determined and it is added into the nonlinear dynamic model.The influence of bearing clearance to the system is discussed.The displacement of cutter tip is accumulated into cutting width in the model of cutting force.Transient response of cutter tip is predicted and analyzed under cutting force as well.

Nonlinear dynamic model of the spindle-cutter coupling system
Fig. 1 shows the typical spindle-cutter coupling system for machine tools.The rotor of spindle is supported by two pairs of angular contact ball bearings.An integral induction motor is located between the front and rear bearings to drive the spindle system.Since the main vibration affecting machining precision is that normal to the machining plane, the proposed model only considers the vibration in this direction.Therefore, spindle-cutter system is simplified to a five-degree-offreedom system as shown in Fig. 2. Where, , and are discrete masses of the spindle; and m are discrete masses of the milling cutter; and are stiffness of spindle. is coupling stiffness at spindle-cutter interface determined by the semi analytical method combining classic elasticity theory with Yoshimura integral method [12]; is stiffness of milling cutter; and are damping coefficient of spindle; is damping coefficient of milling cutter; and are bearing load.
In this paper, the rotation speed = 20044 r/min.The equivalent mass ( ( =1, 2,…, 5)) and equivalent stiffness ( ( = 1, 2, 3)) are identified by the equivalent method described in [13].The damping coefficient is assumed to be proportional to the corresponding stiffness, and = 0.001 ⁄ ( = 1, 2, 3).According to the second Newton's law, dynamic differential equation of spindle-cutter system is as follows: where , , , and are displacements of nodes corresponding to the mass shown in Fig. 2; the downward direction is the positive direction.= 3.21×10 -4 mm is the offset of unbalance mass of spindle, and is the gravity acceleration.

Cutter stiffness
The cutter is usually simplified as cantilever beam [14][15][16][17][18][19], and each cross section along the axial direction is assumed as the same.In fact the error between each cross section exists because of an angle of Δ , which is depicted in Fig. 3-4.A certain error will occur when the cross section is assumed to determine the cutter stiffness.To avoid the error, it is thus necessary to develop a method to evaluate the cutter stiffness considering the difference between the cross sections.All the cross sections of cutter have the same shape and each section rotates an angle of along the cutting edge relative to the neutral plane shown in Fig. 3. Defining the inertia moment relative to the neutral plane at each cross section is thus required to calculate the stiffness of the whole milling cutter.
In this research, cutter stiffness is determined by the analysis of the cross-section.The cutter has helicoidal shape, the cross section is assumed according to [20] and it is shown in Fig. 4.
A two-flute end mill is chosen as the research object, whose cross section rotates an angle of θ.The distance between the section and the initial section is tan .The cross section of cutter can be simplified as in Fig. 4. The parameter is = + , where and are the radius depicted in Fig. 4. The inertia moment of cross section that rotates an angle of relative to the neutral plane is given by: From Eq. ( 2), it can be concluded that the inertia moment relative to the neutral plane is constant, although the cross section rotates an angle of along the cutting edge.The inertia moment of cross section is independent of the angle rotated from the initial section.
The end mill consists of cutting edge and shank, and it can be simplified as the non-uniform beam shown in Fig. 5.According to classic theories of material mechanics [21], the displacement at cutter tip under a unit force applied on the cutter tip can be obtained.The stiffness at cutter tip (i.e.cutter stiffness) is the reciprocal of the displacement, and it can be determined by: where is length of cutting edge, is length of shank, is inertia moment of cutting edge, = ( − ) 4 ⁄ , I is inertia moment of shank, = ( + ) 4 ⁄ .The parameters of the end mill adopted in the model are listed in Table 1.

Nonlinearity of the bearing clearance
Angular contact bearings are generally used in machine tool spindle.A micro-gap may exist between the balls and the outer or inner ring because of the installation error or wear during spindle operation.
The support force of the bearing is expressed as [22]: where is support force of the bearing, is displacement of bearing, is clearance of bearing, is number of balls, is diameter of ball, is preload, and α is contact angle of bearing.The front bearing of the spindle is assumed as loose fit.

Modeling of transient cutting force
According to [23], the transient cutting force can be predicted as follows: where and is transient cutting force in and direction, respectively.is feed of each tooth, is diameter of end mill, is helix angle of end mill, and are coefficient of cutting force in tangential and radial direction, respectively.The coefficients represent the cutting force per unit thickness and unit width, respectively and they are experimentally evaluated.is number of tooth of the cutter.
is lag angle in the maximum cutting deep, is lag angle of micro element relative to cutter tooth on the bottom, is transient position angle of the th tooth, which can be expressed by: where ( ) is a function considering effects of continuous cutting on the cutting force: where and is angle of cutting in and cutting out, respectively.During the down milling = −arccos(1 − 2 ⁄ ) and = 0.The displacement of cutter tip is accumulated into the cutting width.
= + , where is initial cutting width.In this way, the nonlinear factor of cutting force can be considered in this research.
The coefficients of cutting force and have been measured in our previous research [24] and their values are the following: =1880.0N/mm 2 , = 893.3N/mm 2 .

Results and discussion
Dimensionless variables: The dimensionless governing equation of motion of spindle-cutter system is a set of ordinary differential equations with nonlinear and nonsmooth excitation.For this system, the analytical method to study nonlinear dynamic characteristic is not available.The fourth-order Runge-Kutta method is adopted to solve Eq. (8).The relatively shorter time step is chosen to reduce the calculation error; meanwhile a certain periods are abandoned to obtain stable transient solution of the system.In this research, the time step /100 is chosen and the first 1000 periods are abandoned.
To study the effect of bearing clearance on nonlinear dynamic characteristic of the system, responses of cutter tip under different bearing clearances (i.e.= 0, = 0.5 , = , = 5 ) with the same unbalanced force and no cutting force are compared and the results are shown in Fig. 6.
The comparison highlights that each Poincare section under different bearing clearances is a single point; therefore, the responses of cutter tip under different bearing clearances are all periodic motions.The dimensionless displacement of cutter tip ranges [0.787, 0.824] when = 0.It ranges [0.918, 0.955] when = 0.5 , and it ranges [1.003, 1.081] when = 5 .The equilibrium position varies with the increase of bearing clearance.It can be concluded that the bearing clearance mainly has effects on equilibrium position of vibration, and a greater bearing clearance leads to more serious deviation of the equilibrium position from its original position.
When = , bearing clearance changes not only the equilibrium position but also the vibration mode of cutter tip.Although the vibration is still a periodic motion, three peaks occur in a period.
To investigate the effects of clearance on the motion of cutter tip, the dimensionless bearing clearance = / and frequency ratio = / are introduced to obtain the bifurcation diagram, where is the critical rotation speed of the system.In the process of numerical simulation, the cutting force is not considered.The numerical results with = 0.8 is shown in Similar investigation conducted with = 1.0.The numerical results are depicted in Fig. 8.The parameter is chosen as the control parameter to generate the bifurcation diagram, which is shown in Fig. 8(a).From the diagram, the conclusion can be drawn that the motion of cutter tip is a periodic orbit in the range ∈ [0.01, 0.31], which is confirmed by Poincare section depicted in From Fig. 7 and Fig. 8, it is concluded that the responses of cutter tip exhibit periodic, quasi-periodic and chaotic characteristics with different values of bearing clearance and rotation speed.The dynamic characteristic of spindle-cutter system strongly depends on the bearing clearance and rotation speed.Therefore, the proper bearing clearance and rotation speed should be chosen to ensure a stable cutting and high cutting rate according to the bifurcation diagram.According to Fig. 9, the response of the cutter tip is similar to a periodic motion when the cutting force is considered.The amplitude of time-domain response is in a certain range; the phase trajectory is a curve with a certain width.The time-domain response of cutter tip predicted by nonlinear dynamic analysis can provide the basis for the machining error prediction.The dynamic XIANGSHENG GAO, YIDU ZHANG, MIN WANG system strongly depends on the bearing clearance and rotation speed.The proper bearing clearance and rotation speed should be chosen to ensure a stable cutting and high cutting rate according to the bifurcation diagram.The response of the cutter tip is a quasi-periodic motion when the cutting force is considered.The time-domain response of cutter tip predicted by nonlinear dynamic analysis can provide the basis for machining error prediction.

Fig. 5 .
Fig. 5. Model of cutter for stiffness calculation DYNAMIC CHARACTERISTIC OF THE SPINDLE-CUTTER SYSTEM.XIANGSHENG GAO, YIDU ZHANG, MIN WANG

6 .
DYNAMIC CHARACTERISTIC OF THE SPINDLE-CUTTER SYSTEM.XIANGSHENG GAO, YIDU ZHANG, MIN WANG g) Time-domain response h) Poincare section i) Phase trajectory j) Time-domain response k) Poincare section l) Phase trajectory Fig.Comparison of responses of cutter tip under different bearing clearances: a), b)

Fig. 7 . 3 Fig. 7 .
Fig.7.The dimensionless bearing clearance γ is chosen as the control parameter to generate the bifurcation diagram, which is depicted in Fig.7(a).From the bifurcation diagram, it is concluded that the motion of cutter tip is a periodic orbit in the range ∈ [0.6, 1.3], which is confirmed by the result in Fig.7(b).With the increasing of the bearing clearance, the response of cutter tip is chaotic in the range ∈ [1.35, 1.5].Poincare section with = 1.45 is shown in Fig.7(c), which can confirm this response is chaotic.As the bearing clearance increases further, the response regains a periodic orbit in the range ∈ [1.55,1.6].This conclusion can be confirmed by Poincare section depicted in Fig.7(d).As increases beyond 1.65, the response of cutter tip becomes chaotic again, which is confirmed by the result in Fig.7(e).
Fig.8(b).As the increasing of , the response of cutter tip develops into a quasi-periodic orbit.Poincare section depicted in Fig.8(c)-(e) can confirm the conclusion.As increases beyond 3.61, the response of cutter tip becomes a periodic orbit again, which is confirmed by Poincare section shown in Fig.8(f).From Fig.7and Fig.8, it is concluded that the responses of cutter tip exhibit periodic, quasi-periodic and chaotic characteristics with different values of bearing clearance and rotation speed.The dynamic characteristic of spindle-cutter system strongly depends on the bearing clearance and rotation speed.Therefore, the proper bearing clearance and rotation speed should be chosen to ensure a stable cutting and high cutting rate according to the bifurcation diagram.

4 Fig. 8 .
a) Bifurcation diagram b) Poincare section with = 0.16 c) Poincare section with = 0.6 d) Poincare section with = 1.6 e) Poincare section with = 3.2 f) Poincare section with = Nonlinear dynamics of cutter tip with = 1.0

Table 1 .
Parameters of the end mill