Shock and Vibration

Volume 2014 (2014), Article ID 787518, 10 pages

http://dx.doi.org/10.1155/2014/787518

Research Article

On the Free Vibration Modeling of Spindle Systems: A Calibrated Dynamic Stiffness Matrix

Omar Gaber and Seyed M. Hashemi

Department of Aerospace Engineering, Ryerson University, 350 Victoria Street, Toronto, ON, Canada, M5B 2K3

Received 17 October 2013; Accepted 17 March 2014; Published 12 May 2014

Academic Editor: Valder Steffen Jr.

Copyright © 2014 Omar Gaber and Seyed M. Hashemi. 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.

Abstract.

The effect of bearings on the vibrational behavior of machine tool spindles is investigated. This is done through the development of a calibrated dynamic stiffness matrix (CDSM) method, where the bearings flexibility is represented by massless linear spring elements with tuneable stiffness. A dedicated MATLAB code is written to develop and to assemble the element stiffness matrices for the system’s multiple components and to apply the boundary conditions. The developed method is applied to an illustrative example of spindle system. When the spindle bearings are modeled as simply supported boundary conditions, the DSM model results in a fundamental frequency much higher than the system’s nominal value. The simply supported boundary conditions are then replaced by linear spring elements, and the spring constants are adjusted such that the resulting calibrated CDSM model leads to the nominal fundamental frequency of the spindle system. The spindle frequency results are also validated against the experimental data. The proposed method can be effectively applied to predict the vibration characteristics of spindle systems supported by bearings.

1. Introduction

The booming aerospace industry and high levels of competition have forced companies to constantly look for ways to optimize their machining processes. Cycle time has been a major concern at various industries dealing with manufacturing of airframe parts and subassemblies. When trying to machine a part as quickly as possible, spindle speed or metal removal rates are no longer the limiting factors but rather the chatter that occurs during the machining process. It is well established that chatter is directly linked to the natural frequency of the cutting system, which includes the spindle, shaft, tool, and hold combination.

The first mention of chatter can be credited to Taylor in 1907 [1]; however, the first comprehensive investigation of the problem was published by Arnold, in 1946 [2]. Arnold’s experiments were conducted on the turning process, where he theorized that the machine could be modeled as a simple oscillator and that the force on the tool decreased as the speed of the tool in relation to the work piece increased. Gurney and Tobias later theorized the now widely accepted belief that chatter is caused by wave patterns traced onto the surface of the work piece by preceding tool passes [3]. The phase shift of the preceding wave to the wave currently being traced determines whether there is any amplification in the tool head movement. If there exists a phase shift between the two tool paths, then the uncut chip cross-sectional area is varied over the pass. Thomas and Shabana [4] showed that the cutting force is dependent on the chips cross-sectional area, and, as a result, a varying cutting force is produced. This was demonstrated through an illustrative example of a grinding machine, modeled as a one-degree-of-freedom mass-spring system, as opposed to an oscillating system. The spring-mass system is also the widely used modeling theory for how a vibrating tool should be characterized today.

Prior to 1961, the machining procedures were regarded as steady state and discrete processes [5]. This erroneous idea led to the creation of overly heavy and thick-walled machines, believed to provide high damping to the tooltip forces (assumed to be static). However, one must note that machining is a continuous, dynamic process, with constantly fluctuating tooltip forces. Moreover, while being identical in construction, the dynamics and response of similar machines may be slightly different due to structural imperfections, imbalances, and so forth. Therefore, the specific characteristics and dynamics of each machine must also be taken into consideration when performing chatter-related calculations [6]. The structural modes of the machine determine the frequency and the direction that the tool is going to vibrate at [7]. Eisele [5] and Peng et al. [8] stressed that designers must investigate the mode-shapes, weak points, bearing clearances, and self-inducing vibratory components of the machine to reduce chatter. Butlin and Woodhouse also investigated the number of structural modes required to produce accurate results [9]. They showed that low order models, incorporating only two modes, are sufficiently accurate to model the machines.

In 1981, Tlusty and Ismail published one of the first papers that documented the nonlinearity of the vibratory system occurring during chatter [10]. Kondo et al. [11] added to Tlusty and Ismail’s investigation the behaviour of the tool after the onset of chatter. An important study, carried out by Lee et al. in 1989, helped bring a more complete theory of chatter to light [12], where they looked into the response of the work piece to chatter. They modeled the tool as a typical two-degree-of-freedom spring-mass system and the work piece as an Euler beam. When the responses of both the tool and the work piece were incorporated into the typical spring-mass equations, the results obtained far surpassed the accuracy of previously used equations. They observed a large increase in surface finish quality while machining, even with the use of long end mills.

While a comprehensive set of mathematical models have already been derived to explain chatter using spring-mass and Euler beam systems, which have led to dramatic increases in surface finish quality, chatter is still an issue in the machining industry as the surface tolerances keep increasing. From the 1990s on, this drove researchers to use more refined mathematical models to characterize the machines [13]. With the rise of computer aided calculations and simulations, it has become increasingly easier to perform more complex analyses of the machining process [14]. Tool wear is an often-overlooked factor that contributes to chatter. With the aid of more powerful computers, this variable can now be included in simulations. Clancy and Shin [15] investigated the wear of turning tools and how that affects the stability of the system. It was found that as the tool becomes worn, its limits of stability increase. As a result, stability lobes have now become a function of tool wear. Li et al. [16] determined that the coherence function of two crossed accelerations in the bending vibration of the tool shank approaches unity at the onset of chatter. A threshold needs to be set [17] and then detected using simple mechanism to alert the operator to change the machining conditions and avoid increased tool wear.

In most of the existing stability prediction methods, a frequency response function (FRF)—acquired using an impact test—is required to extract the system fundamental frequency, used to perform the calculations. In order to perform accurate impact tests, the machine would need to be shut off and allowed to cool down and the coolant drip off. This process could lead to several hours of down time. An offline method of obtaining the FRF could greatly benefit manufacturing companies by eliminating the downtime needed for the impact tests. While stability lobe graphs have become more accurate, the industry is trending towards computer simulations and numerical analyses. With computing power becoming relatively cheap and readily available, this method has become feasible. Providing the most accurate results due to its constantly updated stability calculations, this is the area of research that engineers will most likely progress towards. Lee et al. [18] suggested that the whole system comprising machine, tool and work piece could be modeled using finite element method (FEM) analysis. A computer simulation would be able to predict the FRF during all phases of the machining process. As the part is machined and becomes thinner, its response to vibration changes dramatically. Most existing research papers assume that the FRF remains constant throughout the whole process. However, a constantly updated FRF would allow for accurate, real time stability calculations.

In an attempt to numerically model the entire spindle systems for the vibration modeling and analysis and to incorporate the effects of bearings, Cao and Altintas [19] modeled spindle bearing and machine tool systems for virtual simulation of milling operations, where they presented a general, integrated model of the spindle bearing and machine tool system, consisting of a rotating shaft, tool holder, angular contact ball bearings, housing, and the machine tool mounting. The proposed model, also verified experimentally, allows virtual cutting of a work material with the numerical model of the spindle during the design stage and predicts bearing stiffness, mode shapes, frequency response function (FRF), static and dynamic deflections along the cutter and spindle shaft, and contact forces on the bearings with simulated cutting forces before physically building and testing the spindles. The study showed that the accuracy of predicting the performance of the spindles requires integrated modeling of all spindle elements and mounting on the machine tool. Jones’ bearing model, which considers the bearing balls and rings as elastic parts, and the Hertzian contact theory were used to calculate the contact force and displacement. Patwari et al. [20] presented a detailed systematic procedure for experimental and analytical modal analysis techniques, to evaluate structural dynamics of a vertical machining centre. Abuthakeer et al. [21] carried out the dynamic characteristics analysis of high speed motorized spindles. More recently, Delgado et al. [22] used FEM to study thenonlinear dynamic behavior of a high speed machine tool spindle.

It is known that FEM predictions and the numerical results are often different from test and experimental results. However, based on the model updating, corrections can be made on the FEM models by processing records of test results. Motthershead and Friswell [23] presented a survey of model updating in structural dynamics. Later, Mares et al. [24] investigated the model updating using robust estimation, and more recently Mottershead et al. [25] presented a basic introduction to the most important procedures of computational model updating, including tutorial examples and a large scale model updating example of a helicopter airframe and discussed the sensitivity method in FEM model updating.

To the authors’ best knowledge, except authors’ previous works [26, 27], no studies on the calibrated/updated dynamic stiffness matrix (DSM) models of spindle systems have appeared in the open literature. The DSM method [28] provides an analytical solution to free vibration problems, achieved by combining the coupled governing differential equations of motion of the system into a single higher-order ordinary differential equation. Seeking the most general solution to the resulting higher-order ODE and enforcing the end displacements and loads through extensive matrix manipulations lead to the DSM of the structural element. Unlike limited nodal DOFs in conventional FEM [29], the DSM-based element is modeled as a continuous system with infinite number of DOFs within the element. The DSM formulation has been proven to provide exact frequency (within the limits of the theory) for various beam configurations and beam-structures, with the use of a single continuous element per segment (see, e.g., [28, 30–33]). The objective of the present study is to determine the natural frequencies/vibration characteristics of machine tool spindle systems by developing the corresponding dynamic stiffness matrix (DSM) [31, 32] and the proper boundary conditions. These results would then be compared to the existing results for a common cutting system to validate/tune the model developed. The coupled bending-bending (B-B) vibration of a spinning beam is investigated. A MATLAB code is developed to assemble the DSM element matrices for multiple components and to apply the boundary conditions (BC). The bearings are first modeled as simply supported (S-S) frictionless pins, which are then replaced by linear spring elements to incorporate the flexibility of bearings into the model. In comparison with the existing data on spindle’s fundamental frequency, the bearing stiffness coefficient, , is then varied to achieve a calibrated (or updated) dynamics stiffness matrix (CDSM) vibrational model. The proposed formulation can also be extended to include torsional degree-of-freedom (DOF) for further modeling purposes.

2. Problem Description and Governing Equations

Computer numeric control (CNC) machines are becoming the norm in manufacturing plants worldwide. These machines are generally 3-, 4-, or 5-axis, depending on the number of degree-of-freedoms the device has. Having the tool translating in the , , and directions accounts for the first three degrees of freedom. Rotation about the spindle axes account for any further DOF. A typical spindle contains the motors that rotate the shaft connected to the tool holder, the bearings, the tools, and all the mechanisms that hold the tool in place (see Figure 1 for a sample spindle configuration). As mentioned earlier in this paper, the fundamental natural frequency of this system dictates the limitation of the cutting parameters. If these limits are surpassed, the system tends to vibrate which causes chatter. Also, as the bearings wear, the system characteristics and spindle’s natural frequency change.

Figure 1: Typical spindle configuration.

In what follows, the derivation of the governing differential equations for the coupled bending-bending vibration of a spinning beam for a rotating shaft is first briefly discussed. Bearings restrain the spinning motion of the spindle system in five degrees of freedom, that is, displacement in , , and directions and the rotation in and directions. The beam torsional rigidity, , is large enough so that the torsional vibrations can be ignored [31, 32]. Figure 2 shows a cylindrical beam in a right-handed rectangular Cartesian coordinates system. The beam length is denoted by , mass per unit length is , polar mass moment of inertia per unit length is , and the principal axes bending rigidities are for both planes. At an arbitrary cross section along -axis, and are displacement of a point in the and directions. The cross section is allowed to rotate or twist about the axis. The position vector of the point after deformation is given by (refer to Figures 2 and 3)

where and are unit vectors in the and directions. The velocity of point is given by

The kinetic and potential energies of the beam ( and ) are given by

Using the Hamilton principle in the usual notation state,

where and are the time intervals in the dynamic trajectory and is the variational operator. Substituting the kinetic and potential energies of the beam in the Hamilton principle, collecting like terms and integrating by parts, the following differential equations are obtained:

which govern the bending-bending (B-B) vibration of a beam, coupled by the spinning speed. The resulting loads, also force boundary conditions at free ends, are then found to be in the following forms, written for shear forces as

and for bending moments as

Figure 2: Spinning beam.

Figure 3: Degrees of freedom of the system.

Assuming simple harmonic motion,

and substituting (8) into (5), they can be rewritten as

where is frequency of oscillation and , , and are the amplitudes of , and . Introducing and , which are nondimensional length and the differential operator, and back substitution into (9) lead to

3. Solution of the Governing Equations

3.1. Dynamic Stiffness Matrix (DSM) Solution

As suggested by Banerjee and Su [31, 32], the above equation (10) can be combined into one 8th-order differential equation written as

where “” represents both lateral displacements and and

The general solution of the differential equation is sought in the form

Substituting (13) into (11) leads to

where , , , , and

Based on Euler-Bernoulli beam bending theory, the corresponding bending rotation, that is, slope, about -axis and -axis, and , respectively, is given by (see the appendix for further details)

By doing similar substitutions in load (3) and (4), one finds

Substituting the boundary conditions into the governing equations, enforced at () and (), leads to

where

Substituting similarly for the force equation one finds

where

The frequency-dependent, dynamic stiffness matrix (DSM) of the spinning beam, , can then be derived by eliminating . We also know that the force amplitude is related to the displacement vector by

Explicit matrix forms of , , and are given in the appendix.

3.2. Classical Finite Element Method (FEM)

For validation and comparison purposes, a conventional finite element formulation of the problem and a MATLAB-based FEM code were developed to carry out the free vibration analysis of the spindle system [33]. The Galerkin method of weighted residuals [29, 33] was employed to develop the integral form of the governing equation (10). Performing integration by parts on the resulting integral equations leads to the weak integral form of equations, which also satisfy the principle of virtual work. The system is then discretized using beam elements with four DOF per node (i.e., two lateral displacements, and , and two slopes, and , eight degrees of freedom in total). Introducing cubic Hermite-type polynomial approximations [29, 33] for bending displacements in the element weak integral form, the classical finite element formulation is developed. This formulation results in constant (or static) element stiffness, mass, and coupling matrices, which are then assembled to find the global matrices. Finally, enforcing the boundary conditions leads to a linear eigenvalue problem, which is solved to extract the system natural frequencies and modes. There has also been an attempt to develop a dynamic finite element (DFE) [34] formulation for the problem in hand, using frequency-dependent interpolation functions, and the preliminary results were found to be promising [35].

4. Application of the Theory and Numerical Results

The DSM formulation was first validated against various analytical and numerical frequency data available in the literature [31, 32] for various uniform simple beam configurations. The DSM frequency data were also compared with those obtained from commercial FEM-based software ABAQUS and the effect of various element types on the convergence was also examined. In addition, the conventional FEM formulation [29, 33] was used for further comparisons; however, as the in-house FEM-based code and commercial software both led to similar results, the FEM results have been omitted here for the brevity.

A typical spindle system model (see Figure 4) was then constructed using multiple segments. As the nominal fundamental frequency was available for nonspinning spindle, therefore the effect of spinning speed was taken out by simply setting to zero in FEM models. However, one should note that setting the spinning speed exactly to zero would cause the DSM formulation to collapse. This is due to the fact that the DSM formulation pivots around combining the two expressions (10) into a single on (11), through the coupling terms involving [31, 32]. However, an extremely small number can be used to model the stationary spindle.

Figure 4: Modified boundary conditions.

The entire system was assumed to be made from the same material, that is, tooling steel. Initially it was also assumed that the system is simply supported at the bearing locations, that is, the bearings modeled as frictionless pins, rigidly attached to the spindle casing with zero flexibility, allowing only for rotations in all directions. This system was modeled using the DSM method described above, as well as finite element analysis (FEA) in ABAQUS commercial software to evaluate the first six natural frequencies of the spindle system. Several different element types and mesh sizes were examined to ensure the convergence (Table 1); beam element B33 uses a 2-node cubic beam (similar to the in-house FEM code), C3D20R is a general-purpose quadratic brick element, and C3D10M is a general-purpose tetrahedral element (see Table 1 for results). Unlike the conventional FEM formulation, the DSM method results in exact solution for all natural frequencies of a spinning uniform beam [31, 32]. Therefore, as the spindle system in hand is made of 12 uniform sections, only a 12-element piecewise-uniform (stepped) DSM model was required. As expected, the results using beam element type B33 yielded the best agreement with the DSM results, as the 3D and shear effects are neglected in both. However, the fundamental frequency values were found to be different from the nominal value, that is, 1393 Hz compared to  Hz, as provided by the manufacturer.

Table 1: Spindle natural frequencies; DSM versus FEA for a nonspinning ( RPM), with simply supported boundary conditions.

Mode
DSM
FEA B33 element (Hz)
FEA C3D20R (Hz)
FEA C3D10M (Hz)

12 sections
600 elements
1200 elements
16965 elements
20122 elements
34465 elements
45645 elements

1
1392.92
1393.00
1393.00
1308.00
1308.00
1301.00
1301.00

2
1393.24
1393.00
1393.00
1308.00
1308.00
1301.00
1301.00

3
4071.02
4071.00
4071.00
3208.00
3212.00
3335.00
3322.00

4
4071.34
4071.00
4071.00
3210.00
3212.00
3335.00
3324.00

5
5668.14
5669.00
5669.00
4339.00
4311.00
4322.00
4291.00

6
5668.46
5669.00
5669.00
4341.00
4316.00
4323.00
4292.00

In order to address this issue, the DSM and FEM models were calibrated to take into account the flexibility of the bearing systems. To this end, the simply supported boundary conditions were modified and replaced by simple linear spring elements, attached to the spindle shaft in and directions, representing general flexibility of the bearing systems (see Figure 4). The bearings were all assumed to be identical, with the same spring coefficient, . The spring stiffness value, , was then varied and a calibrated model equipped with spring elements of constant  N/m was found to result in a fundamental bending frequency of 1004 Hz, closely equivalent to the system’s nominal value (Figure 5). The effect of spinning speed, , on the system’s fundamental frequency was also investigated. In the absence of spinning speed, the system exhibits uncoupled behavior (see (7) and (8)), undergoing two separate bending displacements in and directions. As a result, for a circular shaft with , the modal analysis results in repeated frequency values, appearing in pairs (Table 1). By introducing the spinning speed, , the bending natural frequencies in and directions become coupled, changing the repeated natural frequencies to two distinct ones (see Table 2), with one decreasing and approaching zero (Figure 6) and the other increasing from their uncoupled values. Further increase in the spinning speed would make one of the bending frequencies go to zero; that is, whirling instability, also reported and discussed in [31, 32]. The first critical spindle speed, in the case of calibrated model, was found to be  RPM, which is well above the operating speed of the spindle, that is, 3.5 × 10⁴ RPM (see Figure 7).

Table 2: Spindle natural frequencies versus various spinning speeds, for flexible boundary conditions;  N/m.

Mode
9.55

(RPM)
4784.2

(RPM)
9558.9

(RPM)
14334

(RPM)
19108

(RPM)
23883

(RPM)
28657

(RPM)
33432

(RPM)

1
1004.27
924.69
845.11
765.54
685.96
606.38
526.80
447.23

2
1004.59
1084.16
1143.21
1063.63
984.06
904.48
824.90
745.32

3
1302.37
1222.79
1163.741
1134.14
1054.56
974.98
895.41
815.83

4
1302.68
1293.29
1213.716
1243.32
1322.90
1402.47
1468.05
1388.47

5
1372.87
1382.26
1461.838
1541.42
1541.42
1547.62
1482.05
1561.63

6
1373.19
1452.77
1532.344
1611.92
1627.20
1700.57
1780.15
1859.73

Figure 5: Varying spring constant versus system natural frequency using DSM and FEA.

Figure 6: Spindle speed versus natural frequency.

Figure 7: Spindle natural frequency versus spindle RPM; simply supported and flexible bearing models.

In order to further validate the proposed CDSM model, the tap testing method was used to experimentally collect the frequency response function (FRF) of the spindle system in hand, while being installed on its corresponding machine tool. A 1-inch diameter blank tool, with a 2-inch protrusion, in a typical shrink fit tool holder was used. The tool and holder were inserted into the spindle. The spindle was set to be horizontal (i.e., neutral) position and an accelerometer (352A21, 0.8 g) was then attached to the edge of the tool holder. The spindle was struck 10 times in the direction using a hammer (086C04, 5000 N) and the average FRF graph was generated (SIM3 Module Photon + Data Acquisition). Similarly, the spindle was also tested along the direction. The nonspinning spindle’s natural frequencies extracted from the FRF graphs are presented in Table 3, alongside the CDSM results. As can be seen, the maximum difference between the first three CDSM flexural natural frequencies and the experimentally evaluated data are, respectively, 8%, 0.07%, and 14.3%, in direction, and 8.4%, 0.06%, and 13.3%, in direction. The average difference is found to be around 7.4%, which is within the acceptable range.

Table 3: Spindle natural frequencies; experimental versus CDSM for a nonspinning spindle system, with flexible boundary conditions;  N/m.

Nonspinning bending modes (Hz)
CDSM ( RPM)
Experimental (x-axis/y-axis)
Difference (%) (x-axis/y-axis)

1004.27/1004.59
929.81/926.51
8.0%/8.4%

1302.37/1302.68
1301.51/1310.67
0.07%/0.06%

1372.87/1373.19
1601.44/1583.86
14.3%/13.3%

5. Discussion and Conclusion

It has been established that the service life changes the vibrational behavior of spindle systems, that is, reduced natural frequency over the time, associated with the bearings wear. Therefore, in order to generate an updated spindle model, it is essential to include the effect of the bearings in the system. An analytical model of a multisegment spinning spindle, based on the dynamic stiffness matrix (DSM) formulation and exact within the limits of the Euler-Bernoulli beam bending theory, was developed. The beam exhibits coupled bending-bending (B-B) vibration and, as expected, its fundamental coupled frequency was found to decrease with increasing spinning speed. The bearings were included in the model using two different models: rigid, simply supported, frictionless pins, and flexible linear spring elements. The simply supported bearing model was found to yield an excessively rigid system with a much higher natural frequency than the nominal value of the system. Replacing the simply supported boundary condition with linear spring elements, the spring stiffness, , was then calibrated so that the fundamental frequency of the system matched the nominal value, that is, a more realistic and accurate representation of the spindle system. The accuracy of the proposed calibrated dynamic stiffness matrix (CDSM) method was confirmed through comparison with both numerical and experimental results. When compared to the conventional finite element method (FEM), the CDSM frequencies were found to best match with those obtained from a cubic beam finite element (B33 in ABAQUS). Moreover, in comparison with experimental data, the CDSM results were found to be within acceptable accuracy; the average error for the first three bending natural frequencies is less than 8%.

It is worth noting that the proposed CDSM can also be fine-tuned to match spindle’s experimental fundamental frequency, instead of the nominal value. Further research is underway to improve the CDSM by using more accurate material properties for different segments of the spindle. The developed calibrated model could then be exploited in an attempt to take into account the effects of bearing wear in terms of spindle’s service time, numerically evaluate the equivalent cutting force, and calculate updated stability lobes, used to achieve chatter-free cuts.

Appendix

From the solutions of and (see expressions (13)) the corresponding bending rotation about -axis and -axis, and , are, respectively, given by

By doing similar substitutions we find (see expressions (3), (4), and (14))

where

To obtain the dynamic stiffness matrix we introduce the boundary conditions into the governing equations. For displacements,

 at , 

 at , 

And for forces we have the following:

 at ,

 at ,

Substituting the boundary conditions into the governing equations we find

where

Also

Substituting similarly for the force equation,

where

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors wish to acknowledge the supports provided by Natural Sciences and Engineering Research Council of Canada (NSERC) and Ryerson University.

References

F. W. Taylor, “On the art of cutting metal,” Transactions of the American Society of Mechanical Engineers, vol. 28, pp. 31–350, 1907.

R. N. Arnold, “Discussion on the mechanism of tool vibration in the cutting of steel,” Proceedings of the Institution of Mechanical Engineers, vol. 154, pp. 429–432, 1946.

J. P. Gurney and S. A. Tobias, “A graphical method for the determination of the dynamic stability of machine tools,” International Journal of Machine Tool Design and Research, vol. 1, no. 1-2, pp. 148–156, 1961.

B. Thomas and A. Shabana, “Chatter vibration of flexible multibody machine tool mechanisms,” Mechanism and Machine Theory, vol. 22, no. 4, pp. 359–369, 1987.

F. Eisele, “Machine tool research at the Technological University of Munich,” International Journal of Machine Tool Design and Research, vol. 1, no. 4, pp. 249–274, 1961.

C. Andrew and S. A. Tobias, “Vibration in horizontal milling,” International Journal of Machine Tool Design and Research, vol. 2, no. 4, pp. 369–378, 1962.

T. Insperger, G. Stépán, P. V. Bayly, and B. P. Mann, “Multiple chatter frequencies in milling processes,” Journal of Sound and Vibration, vol. 262, no. 2, pp. 333–345, 2003.

Z. K. Peng, M. R. Jackson, L. Z. Guo, R. M. Parkin, and G. Meng, “Effects of bearing clearance on the chatter stability of milling process,” Nonlinear Analysis: Real World Applications, vol. 11, no. 5, pp. 3577–3589, 2010.

T. Butlin and J. Woodhouse, “Friction-induced vibration: should low-order models be believed?” Journal of Sound and Vibration, vol. 328, no. 1-2, pp. 92–108, 2009.

J. Tlusty and F. Ismail, “Basic non-linearity in machining chatter,” CIRP Annals—Manufacturing Technology, vol. 30, no. 1, pp. 299–304, 1981.

Y. Kondo, O. Kawano, and H. Sato, “Behaviour of self excited chatter due to multiple regenerative effect,” Journal of Engineering for Industry, vol. 103, no. 3, pp. 324–329, 1981.

A.-C. Lee, C.-S. Liu, and S.-T. Chiang, “Analysis of chatter vibration in a cutter-workpiece system,” International Journal of Machine Tools and Manufacture, vol. 31, no. 2, pp. 221–234, 1991.

I. E. Minis, E. B. Magrab, and I. O. Pandelidis, “Improved methods for the prediction of chatter in turning, part 3. A generalized linear theory,” Journal of Engineering for Industry, vol. 112, no. 1, pp. 28–35, 1990.

H. Li and X. Li, “Modelling and simulation of chatter in milling using a predictive force model,” International Journal of Machine Tools and Manufacture, vol. 40, no. 14, pp. 2047–2071, 2000.

B. E. Clancy and Y. C. Shin, “A comprehensive chatter prediction model for face turning operation including tool wear effect,” International Journal of Machine Tools and Manufacture, vol. 42, no. 9, pp. 1035–1044, 2002.

X. Q. Li, Y. S. Wong, and A. Y. C. Nee, “Tool wear and chatter detection using the coherence function of two crossed accelerations,” International Journal of Machine Tools and Manufacture, vol. 37, no. 4, pp. 425–435, 1996.

M. Rahman and Y. Ito, “Detection of the onset of chatter vibration,” Journal of Sound and Vibration, vol. 109, no. 2, pp. 193–205, 1986.

A.-C. Lee, C.-S. Liu, and S.-T. Chiang, “Analysis of chatter vibration in a cutter-workpiece system,” International Journal of Machine Tools and Manufacture, vol. 31, no. 2, pp. 221–234, 1991.

Y. Cao and Y. Altintas, “Modeling of spindle-bearing and machine tool systems for virtual simulation of milling operations,” International Journal of Machine Tools and Manufacture, vol. 47, no. 9, pp. 1342–1350, 2007.

A. U. Patwari, W. F. Faris, A. K. M. Nurul Amin, and S. K. Loh, “Dynamic modal analysis of vertical machining centre components,” Advances in Acoustics and Vibration, vol. 2009, Article ID 508076, 10 pages, 2009.

S. S. Abuthakeer, P. V. Mohanram, and G. Mohan Kumar, “Dynamic characteristics analysis of high speed motorized spindle,” Annals of Faculty Engineering Hunedorar—International Journal of Engineering, 2011, Tome IX, Fascicule 2

A. S. Delgado, E. Ozturk, and N. Sims, “Analysis of non-linear machine tool synamic behaviour,” in Proceedings of the 5th Manufacturing Engineering Society International Conference, Zaragoza, Spain, June 2013.

J. E. Mottershead and M. I. Friswell, “Model updating in structural dynamics: a survey,” Journal of Sound and Vibration, vol. 167, no. 2, pp. 347–375, 1993.

C. Mares, M. I. Friswell, and J. E. Mottershead, “Model updating using robust estimation,” Mechanical Systems and Signal Processing, vol. 16, no. 1, pp. 169–183, 2002.

J. E. Mottershead, M. Link, and M. I. Friswell, “The sensitivity method in finite element model updating: a tutorial,” Mechanical Systems and Signal Processing, vol. 25, no. 7, pp. 2275–2296, 2011.

S. M. Hashemi and O. Gaber, “Free vibration analysis of spinning spindles: a calibrated dynamic stiffness matrix method,” in Advances in Vibration Engineering and Structural Dynamics, F. Beltran-Carbajal, Ed., chapter 4, p. 369, InTech, 2012.

O. Gaber and S. M. Hashemi, “Vibration modeling of machine tool spindles: a calibrated dynamic stiffness matrix method,” Advanced Materials Research, vol. 651, pp. 710–716, 2013.

J. R. Banerjee, “Dynamic stiffness formulation for structural elements: a general approach,” Computers & Structures, vol. 63, no. 1, pp. 101–103, 1997.

K.-J. Bathe, Element Procedures, Prentice Hall, Englewood Cliffs, NJ, USA, 1996.

J. R. Banerjee and F. W. Williams, “Exact Bernoulli-Euler dynamic stiffness matrix for a range of tapered beams,” International Journal for Numerical Methods in Engineering, vol. 21, no. 12, pp. 2289–2302, 1985.

J. R. Banerjee and H. Su, “Dynamic stiffness formulation and free vibration analysis of spinning beams,” in Proceedings of the 6th International Conference on Computational Structures Technology (ICCST '02), pp. 77–78, Civil-Comp Press, Prague, Czech Republic, September 2002.

J. R. Banerjee and H. Su, “Development of a dynamic stiffness matrix for free vibration analysis of spinning beams,” Computers and Structures, vol. 82, no. 23–26, pp. 2189–2197, 2004.

M. Marciniak, Modal analysis of a tool holder system for a five axis milling machine [BEng thesis], Department of Aerospace Eng., Ryerson University, Toronto, Canada, 2012.

S. M. Hashemi, Free vibrational analysis of rotating beam-like structures: a dynamic finite element approach [Ph.D. thesis], Department of Mechanical Engineering, Laval University, Québec, Canada, 1998.

S. Samiezadeh, “Free vibratial analysis of spinning composite beams, using Dynamic Finite Element (DFE) method,” Internal Report, Department of Aerospace Eng., Ryerson University, Toronto, Canada, 2012.