Next Article in Journal
Prediction and Compensation of Color Deviation by Response Surface Methodology for PolyJet 3D Printing
Next Article in Special Issue
A Methodology for Tribo-Mechanical Characterization of Metallic Alloys under Extreme Loading and Temperature Conditions Typical of Metal Cutting Processes
Previous Article in Journal
Development of a Cost Model for Vertical Milling Machines to Assess Impact of Lightweighting
Previous Article in Special Issue
Minimisation of Pose-Dependent Regenerative Vibrations for 5-Axis Milling Operations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Investigation of Gyroscopic Effect on the Stability of High Speed Micromilling via Bifurcation Analysis

1
Department of Mechanical Engineering, Indian Institute of Science, Bengaluru 560012, India
2
Department of Mechanical Engineering, Indian Institute of Technology Bombay, Mumbai 400076, India
*
Author to whom correspondence should be addressed.
J. Manuf. Mater. Process. 2021, 5(4), 130; https://doi.org/10.3390/jmmp5040130
Submission received: 8 September 2021 / Revised: 18 November 2021 / Accepted: 29 November 2021 / Published: 2 December 2021
(This article belongs to the Special Issue Advances in Modelling of Machining Operations)

Abstract

:
Catastrophic tool failure due to the low flexural stiffness of the micro-tool is a major concern for micromanufacturing industries. This issue can be addressed using high rotational speed, but the gyroscopic couple becomes prominent at high rotational speeds for micro-tools affecting the dynamic stability of the process. This study uses the multiple degrees of freedom (MDOF) model of the cutting tool to investigate the gyroscopic effect in machining. Hopf bifurcation theory is used to understand the long-term dynamic behavior of the system. A numerical scheme based on the linear multistep method is used to solve the time-periodic delay differential equations. The stability limits have been predicted as a function of the spindle speed. Higher tool deflections occur at higher spindle speeds. Stability lobe diagram shows the conservative limits at high rotational speeds for the MDOF model. The predicted stability limits show good agreement with the experimental limits, especially at high rotational speeds.

1. Introduction

Micromilling can create complex miniature 3D features with sub-micron accuracy in various engineering materials, such as metals, ceramics, polymers, and composites. The micromilling process is extensively used in many high-technology sectors like biomedical, defense, and aerospace. However, in micromilling, self-excited vibration or chatter may cause catastrophic tool failure and damage the spindle and workpiece. Chatter vibration is detrimental to the dimensional accuracy and surface finish of the workpiece.
The dynamic interaction between the tool and cutting process causes chatter in machining. Regenerative chatter occurs during the cutting process when a flexible tool leaves behind a wavy surface on the workpiece in the preceding cut. The phase difference between the waviness of the current and the preceding passes can result in a condition where the amplitude of the vibration grows rapidly. In the micromilling process, the key reason behind the occurrence of the regenerative chatter is the limited flexural stiffness of the small end mills. Very high spindle speeds can counter the limited flexural stiffness resulting in lower chip loads and cutting forces. However, the spindle dynamics can influence the chatter at high rotational speeds due to gyroscopic effects, especially for low stiffness tools (≤100 µm in diameter and high length to diameter ratio). Therefore, it is essential to investigate the impact of the gyroscopic couple on stability in high-speed machining.
The adverse effects caused by chatter vibration have inspired researchers to develop different models for predicting and avoiding the chatter in conventional milling and micromilling. Tobias [1] and Tlusty [2] presented the mechanism of chip regeneration and stability charts for single-point metal cutting that provides stable depths of cut as a function of spindle speed. Minis and Yanushevsky [3] proposed a theoretical approach based on Fourier analysis and the transfer function. They have developed a comprehensive mathematical model for a milling process with two orthogonal degrees-of-freedom (DOF). Budak and Altintas [4,5] analytically predicted the stability limits in the frequency domain. They also proposed a multi-frequency solution to improve the prediction at small radial immersions where the process is highly intermittent. Insperger et al. [6,7] performed stability analysis in the discrete time domain, including the periodically varying parameters. It is important to note that these studies have considered single or two DOF systems to predict the stability limits.
Several researchers have studied the occurrence of bifurcation during the machining process. Bifurcation is the qualitative behavior of the solutions as the control parameters in a system are varied [8,9]. Davies et al. [10] used the Poincare sectioning technique to characterize and identify chatter in high-speed milling. Zhao and Balachandran [11] developed a time-domain model for both regenerative and loss of contact effects in a milling system. They studied the bifurcations diagrams for parameters, such as axial depth of cut. Schmitz and Honeycutt [12] presented the extended milling bifurcation diagram using time-domain simulation and once-per-revolution sampling. They found both stability behavior such as Hopf and period-n bifurcations in bifurcation diagrams. Consequently, it is crucial to understand the system behavior with a control parameter of interest using bifurcation analysis.
It may be noted that the gyroscopic effect on dynamic stability becomes more pronounced in high-speed micromachining. Movahhedy and Mosaddegh [13] used finite element analysis to determine the modal properties of the spindle system and its frequency response function (FRF). They assumed rigid connections between the tool, tool holder, and spindle. They predicted speed-dependent stability boundaries and observed that the gyroscopic effects become significant only at very high speeds. Bediz et al. [14] used a custom-made impact excitation system to obtain speed-dependent dynamics of high-speed spindles through experimental modal analysis. They noticed that the spindle speed affected the dynamic response significantly, especially at higher spindle speeds. Cao et al. [15] investigated the effect of the gyroscopic moment and centrifugal forces on spindle dynamics. They used the Nyquist stability criterion to generate the stability lobe diagram of high-speed milling. They observed a speed-dependent shift in stability lobes. Xiong et al. [16] investigated the gyroscopic effect of the rotating spindle on the stability characteristics of the conventional milling system and found that it makes the stability prediction less conservative. Wan et al. [17] investigated the effect of gyroscopic couple on the stability. They used speed-dependent FRFs of the tool. They considered the effect of runout on the stability. Mokhtari et al. [18] investigated the effects of rotary inertial dynamics on prediction of chatter in the milling operation. They considered a 3D rotating cantilever Timoshenko beam for modeling the cutting tool. They found that without the rotary inertial dynamics, significant errors are in the prediction of chatter limits at high rotational speeds. Lu et al. [19] obtained the FRF of the spindle system based on the rotating Timoshenko beam theory and the receptance coupling substructure analysis (RCSA). They investigated the centrifugal force and gyroscopic effect caused by the high-speed rotation of the micro-milling spindle. Jalili et al. [20] modeled a cutting tool as a rotating clamped-free beam which is excited by cutting forces. They investigated various kinds of resonances in milling operation using a 3-D nonlinear dynamic model of the milling process including both structural and cutting force nonlinearities, gyroscopic moment, and rotary inertia. Totis [21] introduced additional rotational degrees of freedom in the stability analysis of the conventional milling process. However, the effect of the gyroscopic couple was not investigated in this study. Very few studies have investigated the gyroscopic effect on stability limits in high-speed machining. These models have not explicitly considered the rotational degrees of freedom along with the translational degrees of freedom.
Hence, this study is focused on understanding and modeling the gyroscopic effect in a compliant micro-tool (at high rotational speeds) on the stability boundaries via a rotor dynamics-based analysis with the cutting forces at one end. A MDOF model includes rotational degrees of freedom with the translational degrees of freedom to investigate the effect of gyroscopic moment. Bifurcation analysis was performed to understand the system behavior for a control parameter in the MDOF system. The time-periodic delay differential equations governing the dynamics of the MDOF milling system were numerically solved using a linear multistep (LMS) method. The stability charts and bifurcation diagrams were constructed along with a comparison of displacement profiles at different depths of cut and spindle speeds. Poincare maps were also developed to observe phase space trajectories. The predicted stability limits were validated at different rotational speeds and depths of cut. This study will be beneficial to understand the long-term dynamic behavior of the micromilling system, which is affected by gyroscopic moment especially, at high rotational speeds.

2. Equations of Motion for MDOF Stability Model

For a low flexural stiffness end mill, the effect of the gyroscopic couple will be prominent at high rotational speeds. Both rotational and translational degrees of freedom have been considered in the stability model in the authors’ previous work [22] to capture the effect of gyroscopic couple on stability. This model considers a micro end mill with mass ‘m’ a rotor, as shown in Figure 1. Flexible supports at the ends (clamping end and tooltip end) of the micro end mill are assumed. The center of mass is located at distance ‘a’ and ‘b’ from the clamping end and the tool tip-workpiece interface, respectively. The micro end mill has a fixed frame of reference ‘XYZ’ and a rotating frame of reference ‘xyz’. The micro end mill is rotating around the z-axis with spindle speed ‘ω’. The center of mass has orthogonal translational degrees of freedom in x G and   y G , and rotational degrees of freedom in   θ x and   θ y rotations about the x and y axes, respectively. x 1 , y 1 and x 2 , y 2 are the displacements at the clamping end and tooltip end, respectively, which can be defined in terms of x G ,   y G ,   θ x and   θ y by,
x 1 = x G a θ y ;     x 2 = x G + b θ y y 1 = y G + a θ x ;     y 2 = y G b θ x
The dynamic equations of motion of the center of mass in the translational and rotational directions can be determined using a Lagrangian formulation described in Mittal et al. [22] as,
m x ¨ G + α x x ˙ G + α x x G γ x θ ˙ y γ x θ y = F x
m y ¨ G + α y y ˙ G + α y y G + γ y θ ˙ x + γ y θ x = F y
I x θ ¨ x + I z ω θ ˙ y + γ y y ˙ G + γ y y G + δ y θ ˙ x + δ y θ x = b F y
I x θ ¨ y I z ω θ ˙ x γ x x ˙ G γ x x G + δ x θ ˙ y + δ x θ y = b F x
where
α q = k q , 1 + k q , 2 ,   α q = c q , 1 + c q , 2 γ q = k q , 1 a k q , 2 b ,   γ q = c q , 1 a c q , 2 b δ q = k q , 1 a 2 + k q , 2 b 2 ,   δ q = c q , 1 a 2 + c q , 2 b 2 ;   q = x ,   y
where I x and I z are the transverse and the axial moments of inertia, respectively, k x 1 ,   k x 2 ,   k y 1 ,   k y 2 are the stiffness values of the flexible bearing supports in x and y directions, c x 1 ,   c x 2 ,   c y 1 ,   c y 2 are the damping coefficients of the flexible bearing supports in x and y directions, and Fx and Fy are the cutting forces in the x and y directions. It may be noted that in Equations (4) and (5), I z ω θ ˙ y and I z ω θ ˙ x are the gyroscopic moment acting about x and y axes, respectively.
The derived set of equations of motion for a higher-order multiple degrees of freedom stability model consists of translational and rotational degrees of freedom can be expressed in the following matrix form,
M Φ ¨ + C Φ ˙ + K Φ = F
where M, C, and K are the mass matrix, damping matrix with the gyroscopic component, and stiffness matrix, respectively. Φ is the displacement vector and F is the force matrix, determined using dynamic force modeling. M, C, K, Φ , and F are defined as,
M = m 0 0 0 0 m 0 0 0 0 I x 0 0 0 0 I x ;   C = α x 0 0 γ x 0 α y γ y 0 0 γ y δ y I z ω γ x 0 I z ω δ x ; K = α x 0 0 γ x 0 α y γ y 0 0 γ y δ y 0 γ x 0 0 δ x   ;   Φ = x G y G θ x θ y :   F = F x F y b F y b F x
A mechanistic force was used to determine the dynamic cutting forces Fx and Fy in terms of dynamic chip thickness (as shown in Figure 2) as described by Budak and Altintas [4].
The dynamic chip thickness is given by,
h ϕ j = Δ x 2 t   s i n ϕ j + Δ y 2 t   c o s ϕ j
where Δ x 2 t and Δ y 2 t are the difference between the dynamic displacement during the present (j) and previous (j − 1) cut.
The cumulative cutting forces in matrix form are expressed as,
F t = F x F y = 1 2 a p K t c a x x a x y a y x a y y Δ x 2 Δ y 2
where a x x , a x y ,   a y x , and a y y are the directional factors.
By substituting Fx and Fy and Δ x 2 and Δ y 2 in terms of   x G ,   y G ,   θ x , and   θ y , the force matrix F can be expressed as,
F = F x F y b F y b F x = 1 2 a p K t c a x x a x y b a x y b a x x a y x a y y b a y y b a y x b a y x b a y y b 2 a y y b 2 a y x b a x x b a x y b 2 a x y b 2 a x x Δ x G Δ y G Δ θ x Δ θ y
Equation (10) can be rewritten as
F t = 1 2 a p K t c A t Δ t
where Δ t   is given as Δ t = Φ Φ 0 , K t c is cutting force coefficient in tangential direction and A(t) is the directional factor matrix.
By substituting force matrix in Equation (6), the dynamic equations of motion for a micro end mill system can be expressed in matrix form as,
M Φ ¨ + C ω Φ ˙ + K Φ = 1 2 a p K t c A t Φ t Φ t T

3. Stability Analysis Using Linear Multistep Method (LMM)-Based Numerical Scheme

The multiple degrees of freedom model has been expressed in terms of the periodic delay-differential equation where A(t) is a time-periodic term ( A t = A t + T . The frequency-domain approach will not predict the stability limits as accurately as the time-domain approach. Therefore, in this study, the periodic delay differential equations have been solved numerically using the fourth-order Explicit Adams–Bashforth method, a linear multistep method (LMM) [23,24]. Further, the stability limits have been identified using numerical solutions at different rotational speeds. The dynamic equations of motion for a micro end mill system can be expressed in matrix form as second-order delay differential equation (Equation (12)),
M Φ ¨ + C ω Φ ˙ + K Φ = 1 2 a p K t c A t Φ t Φ t T
where M , C , K , A t , and Φ are the matrices of mass, damping coefficients, stiffness, directional factors, and displacement variables, respectively.
By defining the state variables as,
x G 1 = x G ,     x G 2 = d x G d t = x G 1 ˙ y G 1 = y G ,     y G 2 = d y G d t = y G 1 ˙ θ x 1 = θ x ,     θ x 2 = d θ x d t = θ x 1 ˙ θ y 1 = θ y ,     θ y 2 = d θ y d t = θ y 1 ˙
The second-order delay differential equation for the periodic milling dynamics shown in Equation (12) is organized in the first-order delay differential equation using state variables as follows,
M q ˙ t = L t q t + R t q t T
where
M = I 4 × 4 0 4 × 4 0 4 × 4 M ,   L t = 0 4 × 4 I 4 × 4 δ A t K C ,   R t = 0 4 × 4 0 4 × 4 δ A t 0 4 × 4 , δ = 1 2 a p K t c ,   q t = x G 1     y G 1     θ x 1     θ y 1       x G 2     y G 2     θ x 2     θ y 2 T and q ˙ t = x G 1 ˙     y G 1 ˙     θ x 1 ˙       θ y 1 ˙     x G 2 ˙     y G 2 ˙     θ x 2 ˙     θ y 2 ˙   T
By multiplying the Equation (14) on both sides with the inverse of matrix   M , Equation (14) becomes
q ˙ t = L t q t + R t q t T
where
L t = M 1 L t   and   R t = M 1 R t
The values of state variables in the previous cut are required to solve Equation (15). The first-order delay differential equations as shown in Equation (15) can be written as
q ˙ t = F t , q t , q t T     t 0 q t = φ t     t 0
where φ t is the initial function, T is the time period, q t is called the delay argument, and q t T   is the solution of the delay term.
The linear multistep methods are extensively used for the numerical solution of initial value problems without the delay term. The delay period T is divided into r discrete time intervals of Δ t . Assume the value to q t and F t , q t , q t T   at the current time t i is q i and F i , respectively. According to the fourth-order explicit Adams–Bashforth method, the next value q i + 1 from previous values of q i can be computed as
q i + 1 = q i + Δ t 4 55 F i 59 F i 1 + 37 F i 2 9 F i 3
While solving the delay equation (Equation (15)) using the Adams–Bashforth method, the solution for first interval t   ϵ 0 ,   T will depend on the value of function φ t . The values of q at previous steps will be determined using φ t . For subsequent time intervals, the values of q in previous steps will be used from determined values of q in previous intervals. Using Equation (17), the motion of the micro end mill system (see Equation (12)) can be numerically simulated for different spindle speeds and axial depths of cut.
A schematic to generate the stability plots is shown in Figure 3. The stability limits are predicted by progressively increasing the axial depth of cut for a spindle speed given by ω = 2π/T until the numerical simulations for motion show instability. Root mean square (RMS) values of first half and second half of displacements data were calculated as RMS1 and RMS2, respectively. If the RMS1 is greater than the RMS2, then the system is assumed as stable, otherwise it will be unstable system. This procedure is repeated for the next spindle speeds to construct the stability lobe diagrams.

4. Results and Discussion

The numerical modeling is carried out for slot milling operation. The simulations were performed for different spindle speeds and depths of cut to construct the stability lobe diagrams. In this study, cutting force coefficients were assumed to be constant. To determine the tangential and radial cutting force coefficients, slot milling experiments were conducted on the high-speed micromachining center developed in the Machine Tools Lab, IIT Bombay, as shown in Figure 4. The two-fluted micro end has a diameter of 100 µm. The length of flute is 200 µm. The overhang length of tool (length of tool outside the clamping) is 20 mm. The workpiece material is Ti6Al4V which finds application in the aerospace and biomedical industries. The cutting forces in x and y directions were measured using a 3-axis tabletop dynamometer (Minidyn 9256C1, Kistler, Pune, India). The experimental conditions as shown in Table 1 were used to carry out experiments to determine the cutting force coefficients. Experiments were carried out at feed per rev of 2 µm/flute and depth of cut of 10 µm. The mechanistic force modeling was used to determine the radial and tangential cutting force coefficients from experimentally measured cutting forces. Average value of measured cutting force coefficients for different speeds was used. The value of radial and tangential cutting force coefficients were 9182 MPa and 7205 MPa, respectively.
The modal parameters (natural frequencies, stiffnesses, and damping ratios) at the clamping end and tooltip end are required to predict the stability limits. Impact testing is performed in which excitation is done using a miniature impact hammer (Dytran Chatsworth, CA, USA, Model 5800SL) and displacement measured using a laser displacement sensor (Keyence, Pune, India, Model LK-H020). The stiffness and damping ratios of the machine tool structure at the ends of the rotor, which are assumed as flexible support bearings, are shown in Table 2. The dynamic parameters in x and y directions are assumed to be isotropic because of symmetrical collet and tool.

4.1. Bifurcation Analysis

A bifurcation diagram is a tool to study the long-term qualitative behavior of the dynamic system with the variation of the control parameter. In the stability analysis of milling, the control parameters are axial depth of cut and spindle speed during the machining operation. To construct the bifurcation diagram, simulations were performed for the tooltip displacement at a given depth of cut. Initial values for state variables x G 1   ,   y G 1 ,   θ x 1 ,   θ y 1 ,   x G 2 ,   y G 2 ,   θ x 2   a n d   θ y 2 are 10−6, 10−6, 10−4, 10−4, 0, 0, 0, and 0, respectively. The displacement values are collected for a particular time period and the process is repeated for the next depths of cut. The bifurcation diagrams are numerically constructed for a slot milling operation, as shown in Figure 5, for a higher spindle speed of 120,000 rpm. The control parameter for these bifurcation diagrams is the axial depth of cut. It can be seen that bifurcation occurs at depth of cut ap = 4.1 µm. The system is stable at depths of cut below 4.1 µm and the system shows unstable behavior at depths of cut greater than 4.1 µm.
To find the dynamic behavior of Equation (12), the numerical simulations are performed for different values of bifurcation parameter depth of cut (ap). Figure 6 and Figure 7 show the waveform and Poincare maps for pre-bifurcation (depth of cut ap = 4.1 µm) and post-bifurcation (depth of cut ap = 4.6 µm) motions at a spindle speed of 120,000 rpm. Figure 6 shows that if the axial depth of cut is less than the critical value, then the zero-equilibrium point is asymptotically stable. The amplitude of displacements (see Figure 6a,b) and tool deflections (see Figure 6c,d) approaches zero with time. It can also be seen from the orbital map (see Figure 6e) and the Poincare maps (see Figure 6f) that trajectory moves towards the equilibrium point, and the system is stable. When the axial depth of cut increases beyond 4.1 µm, a Hopf bifurcation occurs at ap = 4.4 µm, as shown in Figure 5. For post-bifurcation axial depth of cut ap = 4.6 µm, the displacement profiles (see Figure 7a,b) and tool deflection profiles (see Figure 7c,d) diverge with time, indicating unstable behavior. The orbital map and Poincare maps as shown in Figure 7e,f respectively, show that trajectory moves away from the equilibrium point, and the system is unstable.
Phase portrait for pre-bifurcation, bifurcation point, and post-bifurcation motions are shown in Figure 8. It can be seen that at the bifurcation point (ap = 4.4 µm), the vibration is not growing or decaying. A depth of cut below 4.4 µm, trajectory moves towards zero (shown in the center of the plot), and the system is stable, but depth of cut greater than 4.4 µm indicates unstable behavior. Figure 9 shows a comparison of the orbital maps ( θ x vs. θ y ) for spindle speeds of 60,000 and 120,000 rpm at one fixed depth of cut ap = 4.1 µm. Variation in θ x and θ y are significantly higher for spindle speed of 120,000 rpm as compared to spindle speed of 60,000 rpm. This can also be seen from the waveform plots as shown in Figure 10a,b. The amplitude of deflection in x and y direction is significantly higher for a higher spindle speed of 120,000 rpm. Root mean square (RMS) of θ x and θ y increases by 37% and 45%, respectively, for spindle speed of 120,000 rpm as compared to 60,000 rpm. The higher deflection at higher speeds shows the presence of the gyroscopic effect in the machining process.

4.2. Stability Lobe Diagram and Validation

The stability lobe diagrams were constructed for slot milling operation (100% immersion) on Ti6Al4V using LMS based numerical scheme. Table 2 shows the modal parameters at tool tip and clamping ends used in the numerical simulations. The feed per flute is fixed at 2 µm for the prediction of stability limits. The delay period T is divided in r = 100, which is a parameter for the discretization of the tooth period. Figure 11 shows the stability limits for the 2DOF and MDOF stability models. In the stability chart, the loss of stability is associated with the Neimark bifurcation of periodic motions, also known as Hopf bifurcation. The MDOF chatter model shows conservative stability limits, especially at high rotational speeds (>100,000 rpm). This conservative behavior of stability limits at high rotational speeds can be attributed to the inclusion of additional rotational degrees of freedom with gyroscopic components. The higher deflection at higher spindles (as shown in Figure 10) shows the effect of gyroscopic couple. Due to effect of gyroscopic couple at higher rotational speeds, there is a reduction in stability limits at higher spindle speeds compared to the lower spindle speeds.
Experiments were carried out to verify the stability predictions at different depths of cut and spindle speeds selected from the stability chart shown in Figure 11. The FFT spectrum of velocity and machined surface data were used to identify the stable and unstable depths of cut experimentally. Figure 12 shows the FFT spectrum of the velocity data at spindle speeds of 50,000 and 120,000 and corresponding depths of cut of 20 µm and 5 µm. It can be seen that there is no high amplitude peak at a rotational speed of 50,000 rpm and a depth of cut of 20 µm which shows stable machining. At a spindle speed of 120,000 rpm and a depth of cut of 5 µm, the amplitude peaks of the velocity spectrum are higher at around 6000 Hz which is in the vicinity of one of the natural frequencies of the tool tip. The critical limits for stability obtained from numerical simulations and experiments are compared in Figure 13. It can be seen that the predicted stability limits from the MDOF model are in better agreement with experimental limits at high rotational speeds (>100,000 rpm). Difference between predicted and experimental limits at high rotational speeds is less than 16%. Hence, using an MDOF model is imperative to predict stability limits at very rotational speeds as the effect of the gyroscopic couple becomes pronounced, especially for miniature tools (~100 µm in diameter).

5. Conclusions

This paper uses a multiple degrees of freedom model based on rotor dynamics to investigate the gyroscopic effect. A numerical scheme based on the linear multistep method has been presented to predict the stability limits. Bifurcation analysis has been performed to understand the dynamic behavior of the system under high rotational speeds. The loss of stability is associated with Hopf bifurcation. Numerical simulations were presented to compare the pre-bifurcation and post-bifurcation motions. It was found that vibration does not grow and decay at the bifurcation point at depth of cut of 4.1 µm at spindle speed of 120,000 rpm. RMS of θ x and θ y increases by 37% and 45%, respectively, for spindle speed of 120,000 rpm as compared to 60,000 rpm. This increase in deflection in x and y direction ( θ x vs. θ y ) is significantly higher for higher spindle speeds which shows the presence of gyroscopic effect at high rotational speeds. It has been found that at high rotational speeds, stability limits are conservative due to the gyroscopic effect at high rotational speeds.

Author Contributions

Conceptualization, R.K.M. and R.K.S.; methodology, R.K.M.; software, R.K.M.; validation, R.K.M.; formal analysis, R.K.M.; investigation, R.K.M.; resources, R.K.M.; data curation, R.K.M.; writing—original draft preparation, R.K.M.; writing—review and editing, R.K.M. and R.K.S.; visualization, R.K.M.; supervision, R.K.S.; project administration, R.K.M.; funding acquisition, R.K.M. All authors have read and agreed to the published version of the manuscript.

Funding

The authors wish to gratefully acknowledge the support of Department of Science and Technology, Govt. of India—Inspire Faculty Fellowship award (DST/INSPIRE/04/2019/000795).

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tobias, S.A. Vibration of Machine Tools. Prod. Eng. 1964, 43, 599. [Google Scholar] [CrossRef]
  2. Tlusty, J.; Polacek, M. The Stability of Machine Tools against Self-Excited Vibrations in Machining. Int. Res. Prod. Eng. ASME 1963, 1, 465–474. [Google Scholar]
  3. Minis, I.; Yanushevsky, R. A New Theoretical Approach for the Prediction of Machine Tool Chatter in Milling. J. Manuf. Sci. Eng. Trans. ASME 1993, 115, 1–8. [Google Scholar] [CrossRef]
  4. Budak, E.; Altintaş, Y. Analytical Prediction of Chatter Stability in Milling—Part I: General Formulation. J. Dyn. Syst. Meas. Control Trans. ASME 1998, 120, 22–30. [Google Scholar] [CrossRef]
  5. Budak, E.; Altintaş, Y. Analytical Prediction of Chatter Stability in Milling—Part II: Application of the General Formulation to Common Milling Systems. J. Dyn. Syst. Meas. Control Trans. ASME 1998, 120, 31–36. [Google Scholar] [CrossRef]
  6. Insperger, T.; Stépán, G. Semi-Discretization Method for Delayed Systems. Int. J. Numer. Methods Eng. 2002, 55, 503–518. [Google Scholar] [CrossRef]
  7. Insperger, T.; Stépán, G.; Bayly, P.V.; Mann, B.P. Multiple Chatter Frequencies in Milling Processes. J. Sound Vib. 2003, 262, 333–345. [Google Scholar] [CrossRef]
  8. Guo, S.; Li, J. Bifurcation Theory of Functional Differential Equations: A Survey. J. Appl. Anal. Comput. 2015, 5, 751–766. [Google Scholar] [CrossRef]
  9. Guo, S.; Wu, J. Bifurcation Theory of Functional Differential Equations; Springer: New York, NY, USA, 2013; Volume 10, p. 978. [Google Scholar]
  10. Davies, M.A.; Dutterer, B.; Pratt, J.R.; Schaut, A.J.; Bryan, J.B. On the Dynamics of High-Speed Milling with Long, Slender Endmills. CIRP Ann.-Manuf. Technol. 1998, 47, 55–60. [Google Scholar] [CrossRef]
  11. Zhao, M.X.; Balachandran, B. Dynamics and Stability of Milling Process. Int. J. Solids Struct. 2001, 38, 2233–2248. [Google Scholar] [CrossRef]
  12. Schmitz, T.L.; Honeycutt, A. The Extended Milling Bifurcation Diagram. Procedia Manuf. 2015, 1, 466–476. [Google Scholar] [CrossRef] [Green Version]
  13. Movahhedy, M.R.; Mosaddegh, P. Prediction of Chatter in High Speed Milling Including Gyroscopic Effects. Int. J. Mach. Tools Manuf. 2006, 46, 996–1001. [Google Scholar] [CrossRef]
  14. Bediz, B.; Ozdoganlar, O.B. Rotational Dynamics of Micro-Scale Cutting Tools. Precis. Eng. 2019, 60, 1–11. [Google Scholar] [CrossRef]
  15. Cao, H.; Li, B.; He, Z. Chatter Stability of Milling with Speed-Varying Dynamics of Spindles. Int. J. Mach. Tools Manuf. 2012, 52, 50–58. [Google Scholar] [CrossRef]
  16. Xiong, G.L.; Yi, J.M.; Zeng, C.; Guo, H.K.; Li, L.X. Study of the Gyroscopic Effect of the Spindle on the Stability Characteristics of the Milling System. J. Mater. Process. Technol. 2003, 138, 379–384. [Google Scholar] [CrossRef]
  17. Wan, M.; Dong, Z.Y.; Yang, Y.; Zhang, W.H. Stability Analysis of Milling Process by Combining the Gyroscopic Effect with the Symmetry and Runout of the Cutter. Mech. Syst. Signal Process. 2021, 161, 107977. [Google Scholar] [CrossRef]
  18. Mokhtari, A.; Mazidi, A.; Jalili, M.M. Investigation of Rotary Inertial Dynamic Effects on Chatter Boundary in Milling Process Using Three-Dimensional Timoshenko Tool Model. Proc. Inst. Mech. Eng. Part K J. Multi-Body Dyn. 2018, 233, 93–110. [Google Scholar] [CrossRef]
  19. Lu, X.; Jia, Z.; Liu, S.; Yang, K.; Feng, Y.; Liang, S.Y. Chatter Stability of Micro-Milling by Considering the Centrifugal Force and Gyroscopic Effect of the Spindle. J. Manuf. Sci. Eng. Trans. ASME 2019, 141, 1–32. [Google Scholar] [CrossRef]
  20. Jalili, M.M.; Hesabi, J.; Abootorabi, M.M. Simulation of Forced Vibration in Milling Process Considering Gyroscopic Moment and Rotary Inertia. Int. J. Adv. Manuf. Technol. 2017, 89, 2821–2836. [Google Scholar] [CrossRef]
  21. Totis, G. Breakthrough of Regenerative Chatter Modeling in Milling by Including Unexpected Effects Arising from Tooling System Deflection. Int. J. Adv. Manuf. Technol. 2017, 89, 2515–2534. [Google Scholar] [CrossRef]
  22. Mittal, R.K.; Kulkarni, S.S.; Singh, R.K. Multiple Degree of Freedom Rotordynamics Based Stability Modeling in High-Speed Micromilling of Ti-6Al-4V. Procedia Manuf. 2018, 26, 607–616. [Google Scholar] [CrossRef]
  23. Jankowska, M.; Marciniak, A. On Explicit Interval Methods of Adams-Bashforth Type. Comput. Methods Sci. Technol. 2002, 8, 46–57. [Google Scholar] [CrossRef]
  24. Marciniak, A.; Jankowska, M.A. Interval Methods of Adams-Bashforth Type with Variable Step Sizes. Numer. Algorithms 2020, 84, 651–678. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Micro end mill as a rotor with flexible bearing supports.
Figure 1. Micro end mill as a rotor with flexible bearing supports.
Jmmp 05 00130 g001
Figure 2. Schematic of milling process.
Figure 2. Schematic of milling process.
Jmmp 05 00130 g002
Figure 3. Schematic diagram to generate stability plots.
Figure 3. Schematic diagram to generate stability plots.
Jmmp 05 00130 g003
Figure 4. High-speed micromachining center.
Figure 4. High-speed micromachining center.
Jmmp 05 00130 g004
Figure 5. Bifurcation diagram for a spindle speed of 120,000 rpm.
Figure 5. Bifurcation diagram for a spindle speed of 120,000 rpm.
Jmmp 05 00130 g005
Figure 6. Waveform plots: (a) XG vs. time; (b) YG vs. time; (c) θ x vs. time; (d) θ y vs. time, orbital map: (e) YG vs. XG and Poincare map: (f) X G 1 ˙ vs. XG for pre bifurcation motion (depth of cut ap = 4.1 µm and spindle speed of 120,000 rpm).
Figure 6. Waveform plots: (a) XG vs. time; (b) YG vs. time; (c) θ x vs. time; (d) θ y vs. time, orbital map: (e) YG vs. XG and Poincare map: (f) X G 1 ˙ vs. XG for pre bifurcation motion (depth of cut ap = 4.1 µm and spindle speed of 120,000 rpm).
Jmmp 05 00130 g006
Figure 7. Waveform plots: (a) XG vs. time; (b) YG vs. time; (c) θ x vs. time; (d) θ y vs. time, orbital map: (e) YG vs. XG and Poincare map: (f) X G ˙ vs. XG for post bifurcation motion (depth of cut ap = 4.6 µm and spindle speed of 120,000 rpm).
Figure 7. Waveform plots: (a) XG vs. time; (b) YG vs. time; (c) θ x vs. time; (d) θ y vs. time, orbital map: (e) YG vs. XG and Poincare map: (f) X G ˙ vs. XG for post bifurcation motion (depth of cut ap = 4.6 µm and spindle speed of 120,000 rpm).
Jmmp 05 00130 g007
Figure 8. Poincare maps for pre-bifurcation (depth of cut ap = 4.1 µm), bifurcation point (depth of cut ap = 4.4 µm), and post-bifurcation motion (depth of cut ap = 4.6 µm) (a) X G ˙ vs. XG and (b) Y G ˙ vs. YG.
Figure 8. Poincare maps for pre-bifurcation (depth of cut ap = 4.1 µm), bifurcation point (depth of cut ap = 4.4 µm), and post-bifurcation motion (depth of cut ap = 4.6 µm) (a) X G ˙ vs. XG and (b) Y G ˙ vs. YG.
Jmmp 05 00130 g008
Figure 9. Orbital maps at a depth of cut ap = 4.1 µm for spindle speeds of 60,000 and 120,000 rpm.
Figure 9. Orbital maps at a depth of cut ap = 4.1 µm for spindle speeds of 60,000 and 120,000 rpm.
Jmmp 05 00130 g009
Figure 10. Comparison of waveform plots (a) θ x vs. time; (b) θ y vs. time at a depth of cut ap = 4.1 µm for spindle speeds of 60,000 and 120,000 rpm.
Figure 10. Comparison of waveform plots (a) θ x vs. time; (b) θ y vs. time at a depth of cut ap = 4.1 µm for spindle speeds of 60,000 and 120,000 rpm.
Jmmp 05 00130 g010aJmmp 05 00130 g010b
Figure 11. Stability lobe diagram for 2DOF and MDOF stability models.
Figure 11. Stability lobe diagram for 2DOF and MDOF stability models.
Jmmp 05 00130 g011
Figure 12. FFT spectrum of the velocity data.
Figure 12. FFT spectrum of the velocity data.
Jmmp 05 00130 g012
Figure 13. Validation of stability lobe diagram.
Figure 13. Validation of stability lobe diagram.
Jmmp 05 00130 g013
Table 1. Experimental conditions for determination of cutting force coefficients.
Table 1. Experimental conditions for determination of cutting force coefficients.
Experimental ConditionsValue
Tool Diameter (µm)100
Cutting edge radius (µm)2–3
Speed (RPM)40,000, 60,000, 80,000
Feed (µm/flute)2
Axial depth of cut (µm)10
Table 2. Stiffness, natural frequency, and damping ratio at the clamping end and the tooltip end.
Table 2. Stiffness, natural frequency, and damping ratio at the clamping end and the tooltip end.
Clamping End
Stiffness (kx1, ky1)9 × 106 N/m
Natural Frequency (Hz)4030 Hz
Damping ratio (ζx1, ζy1)0.03
Tooltip end
Stiffness (kx2, ky2)0.65 × 106 N/m
Natural Frequency (Hz)4000, 4800, 6440 Hz
Damping ratio (ζx2, ζy2)0.03
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mittal, R.K.; Singh, R.K. Investigation of Gyroscopic Effect on the Stability of High Speed Micromilling via Bifurcation Analysis. J. Manuf. Mater. Process. 2021, 5, 130. https://doi.org/10.3390/jmmp5040130

AMA Style

Mittal RK, Singh RK. Investigation of Gyroscopic Effect on the Stability of High Speed Micromilling via Bifurcation Analysis. Journal of Manufacturing and Materials Processing. 2021; 5(4):130. https://doi.org/10.3390/jmmp5040130

Chicago/Turabian Style

Mittal, Rinku K., and Ramesh K. Singh. 2021. "Investigation of Gyroscopic Effect on the Stability of High Speed Micromilling via Bifurcation Analysis" Journal of Manufacturing and Materials Processing 5, no. 4: 130. https://doi.org/10.3390/jmmp5040130

APA Style

Mittal, R. K., & Singh, R. K. (2021). Investigation of Gyroscopic Effect on the Stability of High Speed Micromilling via Bifurcation Analysis. Journal of Manufacturing and Materials Processing, 5(4), 130. https://doi.org/10.3390/jmmp5040130

Article Metrics

Back to TopTop