# Modeling of Self-Vibratory Drilling Head-Spindle System for Predictions of Bearings Lifespan

Authors: F. Forestier; V. Gagnol; P. Ray; H. Paris
Journal: Advances in Acoustics and Vibration
 (2011)

---

## Abstract

The machining of deep holes is limited due to inadequate chip evacuation, which induces tool breakage. To limit this drawback, retreat cycles and lubrication are used. An alternative response to the evacuation problem is based on high-speed vibratory drilling. A specific tool holder induces axial self-maintained vibration of the drill, which enables the chips to be split. The chips are thus of a small size and can be evacuated. To anticipate the potential risk of decreased spindle lifespan associated with these vibrations, a model of the behavior of the system (spindle—self-vibrating drilling head—tool) is elaborated. In order to assess the dynamic behavior of the system, this study develops a rotor-based finite element model, integrated with the modelling of component interfaces. The current results indicate that the simulations are consistent with the experimental measurements. The influence of spindle speed and feed rate on bearing lifespan is highlighted.


---

## Body

## 1. Introduction

High-speed vibratory drilling allows chips to be split thanks to self-maintained vibration during cutting [1]. When these vibrations have a magnitude greater than the advance per tooth, the drill continuously enters and exits the material, which allows the fragmentation of the chips, as shown by Tichkiewitch et al. [2]. The chips are small and can thus be removed easily without retreat cycles or lubrication. High-speed vibratory drilling enables productivity to increase by a factor of three compared to traditional techniques. Currently, high-speed vibratory drilling is the only process capable of producing deep holes using high-speed machining centers, without lubricants, and with high productivity [3].However, the excitation generated by the cutting process can also be a source of damage to the self-vibrating drilling head and the spindle. So the industrialization of high-speed vibratory drilling requires the effects of vibration on the machine, and more particularly on the spindle, to be predicted. For this purpose, a model of the dynamic behavior of the system (self-vibrating drilling head spindle) is elaborated from a realistic assembly of the components, in order to predict spindle bearing lifespan.The SVDH is composed of an axial vibrating system, consisting of the SVDH vibrating subsystem, mounted on a specific HSK 63 taper, and called the SVDH body. The SVDH-vibrating subsystem is composed of the vibrating parts of the drill holder. The SVDH body guides the axial vibrating subsystem through a ball retainer and a classic HSK63 taper connexion with the spindle head. The self-excited vibrations must be tuned and controlled in order to have a magnitude greater than the advance per tooth. Mathematical models of the SVDH dynamics appearing in most of previous works lead to a one-dimensional linear or nonlinear model in the axial vibration direction governed by the mass, damping, and stiffness of the considered system [3]. In these studies, damping, which plays an important role in SVDH dynamic behaviour, is estimated but not experimentally identified. Moreover, the spindle behaviour and the interfaces between the spindle and the SVDH are assumed to be rigid and have not been taken into account. However, many works have showed that tool point dynamics can be significantly influenced by spindle dynamics, as well as the spindle-holder—tool interfaces. This creates a demand for predictive knowledge models that are capable of investigating the influence of cutting conditions on high-speed spindle-SVDH system.Many authors have investigated the dynamic behaviour of machine tool spindle-bearing systems, both analytically and experimentally. They show that spindle dynamics are influenced by a large number of factors, including holder characteristics [4], spindle shaft geometry and drawbar force [5, 6], and the stiffness and damping provided by the bearings [7]. Most of these factors are independent of spindle speed, contrary to bearing properties [8, 9] and the spindle rotor dynamics, which change according to spindle speed. Such rotating systems have been successfully modelled through rotor-dynamics studies [10, 11]. In previous works [12, 13], a dynamic high-speed spindle-bearing system model based on rotor dynamics prediction was presented. Element kinematics were formulated in a corotational coordinate frame and enabled a special rotor beam element to be developed. Model results showed that spindle dynamics are influenced by the gyroscopic coupling and the spin softening of the rotating shaft due to high rotation speeds.The literature on the modelling and analysis of spindles shows that the tool tip FRF is also greatly influenced by the contact dynamics of the spindle-holder-tool interfaces. The flexibility of the afore-mentioned interfaces can dominate the dynamics of the spindle. The tool tip FRF is usually obtained using experimental measurements, which require significant testing time to take into account the large number of spindle, holder, and tool combinations. Thus, semi-analytical approaches have been proposed to minimise experimental approaches. Erturk et al. [14] use a receptance coupling (RC) and structural modification method to connect the spindle shaft and the tool holder. Schmitz et al. [15] model the spindle holder experimentally and couple it with an analytical model of the overhang portion of the tool. Recent approaches consider distributed springs and dampers between the tool and the holder along the interface contact. [15, 16]. Contact stiffness and damping values alter the frequencies and peak values respectively of dominant tool tip vibration modes. The fast and accurate identification of contact dynamics in spindle-tool assemblies has become an important issue. Ren and Beards [17], Schmitz et al. [15], and Movahhedy and Gerami [18] treat identification as a nonlinear optimisation problem. Movahhedy and Gerami [18] uses a genetic algorithm to find the global minimum of the optimisation function(1)|g|=|{Re(gp-gm)Im⁡(gp-gm)}|,
where gp and gmare, respectively, predicted and measured receptances of the assembly. Ozsahin et al. [19] present an original identification procedure based on experimental measurements. In their work, the elastic RC equations allow the stiffness and damping parameters of the spindle-tool assemblies to be obtained in closed-form expressions. In order to predict the drill dynamics and the adequate cutting conditions that lead to controlled self-excited vibrations, an accurate comprehensive dynamic model of the cutting process and spindle-SVDH dynamics is required.In this paper, a hybrid model based on numerical and experimental approaches of the dynamic behavior of the system is proposed. This finite element model takes into account rotor dynamics effects, bearing stiffness and the real behavior of interfaces between different system components. The simulation allows achieving the best possible cutting parameters without damaging the bearings of the spindle.In the second section, the spindle-SVDH rotor dynamics model is presented. A special rotor-beam element is implemented. The rolling bearing stiffness matrices are calculated around a static function point on the basis of Lim and Singh’s [20] formulation and then integrated into the global finite element model. The identification of contact dynamics in tool-SVDH-spindle assemblies is carried out using the RC method on the basis of experimental substructure characterisation. The identified models are then integrated into the global spindle-SVDH-tool model. Finally, numerical and experimental tool tip FRF, in radial and axial directions, is compared in order to validate the global assembled model.Section3 is dedicated to the definition of optimal cutting conditions with respect to industrial objectives. The studied industrial context requires the maximum material cutting rate and a rational use of the tool-SVDH-spindle set in order to guarantee adequate rolling bearing lifespan. As a result, a recommendation for the use of a spindle-SVDH-tool set which respects the defined zone of interest, combining reliability and productivity constraints, is proposed.Finally, conclusions are presented.


## 2. Model Building

The vibratory drilling system is composed of a SVDH body clamped to the spindle by a standard HSK63A tool-holder interface. A SVDH-vibrating subsystem is jointed to the SVDH body with a specific spring and axially guided by a ball retainer. Finally, a long drill is held in the SVDH-vibrating subsystem with a standard ER25 collet chuck. The spindle has four angular bearings in overall back-to-back configuration (Figure1).Figure 1
The spindle-SVDH-tool system and its finite element model.The spindle-SVDH-tool finite element model is restricted to the rotating system composed of the spindle shaft, the SVDH and the drill. An experimental modal identification procedure was carried out on the different spindle substructures and showed that spindle behaviour can be restricted to rotating structure behaviour [12]. The interfaces represented by the HSK63 taper, spring and ball retainer, and collet chuck were taken into account in the model. The CNC milling machine structure was assumed to be infinitely rigid compared to the other parts of the system.The numerical model of the spindle SVDH tool is based on the integration of the rotating system’s finite element model (FEM), the rolling bearing model, and the interface model. Figure2 summarises the various stages of model development. The system substructures were modelled through rotor-dynamics formulations. A readjustment procedure was carried out on undefined FEM material properties in order to fit model results to experimental ones. The receptance coupling method was used to identify the dynamic parameters of the system’s interfaces. Once each structural subsystem model was validated, the identified interface behaviour parameters were integrated to obtain the assembled global model.Figure 2
Development of the numerical model.


### 2.1. Modelling of Structural Subsystems

The spindle SVDH tool system is composed of four structural subsystems: the drill, the SVDH-vibrating subsystem, the SVDH body, and the spindle. The motion of the rotating structure is considered as the superposition of rigid and elastic body displacements. Dynamic equations were obtained using Lagrange formulation associated with a finite element method. Due to the size of the rotor sections, shear deformations had to be taken into account. Then, the rotating substructure was derived using Timoshenko beam theory. The relevant shape functions were cubic in order to avoid shear locking. A special three-dimensional rotor-beam element with two nodes and six degrees of freedom per node was developed in the corotational reference frame. The damping model used draws onRayleigh viscous equivalent damping, which makes it possible to regard the damping matrix D as a linear combination of the mass matrix M and the spindle rigidity matrix K(2)D=aM+bK,
where a and b are damping coefficients.The set of differential equations can be written as(3)M(qN)q̈N+(C(qN,q̇N)+D)q̇N+KqN=F(t),
where M is the mass matrix, and Cmatrix contains the rotational dynamics effects. qN and F(t) are the nodal displacement and force vectors. An accelerating rotor gives rise to previous time-variant equations (3), but treatment of the rotor using a pseudoconstant speed approach can still be described by means of linear time-invariant models and is valid in many cases (4):
(4)M0q̈N+(2ΩG+D)q̇N+(K-Ω2N)qN=F(t),
where M0 is the constant part of matrix M, and G and N proceed from the decomposition of matrix C.


### 2.2. Modelling Angular Contact Ball Bearings

The rotating system is supported by four (two front and two rear) hybrid angular contact bearings. The rolling bearing stiffness matrices were calculated using in-house software developed on the basis of Lim and Singh’s [20] formulation. The bearing stiffness model represents the load-displacement relation combined with the Hertzian contact stress principle and was calculated around a static function point characterised by the bearing preload: δ. Based on Rantatalo’s prediction [9], the initially calculated bearing stiffness is spindle speed dependent because of the gyroscopic and centrifugal force Fc, which acts on each ball (Figure 3(a)). As the speed increased, the load conditions between the balls and the rings in the bearing changed because of the centrifugal force (Figure 3(a)). Then, speed-dependent bearing stiffness was integrated into the global spindle FEM and influenced the natural frequencies of the spindle-SVDH-tool unit under consideration.Bearing stiffness variation depending on spindle speed [12].
(a)
Angular contact ball bearing Hertzian and centrifugal force(b)
Bearing stiffness variation


### 2.3. Structural Systems Dynamic Readjustment

While detailed knowledge of the spindle-SVDH-tool system is in general not available in a manufacturing environment, models need to be readjusted in order to fit experimental results. The readjustment parameters are the Young modulus:E, the damping coefficients: (a,b), and the rolling bearing preload: δ. The readjustment procedure is proposed to tune the previous variable in order to fit the model results to the experimental frequency response function (FRF). These parameters are readjusted by minimizing the gap between the measured and the modelled tool tip node FRF for nonrotating components, using an optimisation routine and a least squares type objective function defined as(5)r=∑ω=ωmin⁡ωmax⁡(Hnum(ω,(E,a,b,δ))-Hexp⁡(ω))2,
where Hnum(ω,(E,a,b,δ)) and Hexp(ω) are, respectively, the numerical and experimental FRF.


### 2.4. Modelling and Identification of Spindle-SVDH-Tool Interfaces

The dynamic behaviour of the interfaces represented by the HSK63 taper, spring and ball retainer, and collet chuck was taken into account. The identification procedure of the interface models was based on the receptance coupling method.


#### 2.4.1. Receptance Coupling Background

In this section, the RC equations are established. At the top of Figure4, the substructure A and the substructure B are represented, connected by interface H. At the bottom of Figure 4, the assembled structure is represented. This figure enables the excitation point and the measurement points used in the receptance coupling approach to be located.Figure 4
Receptance coupling notation.The notationAij(ω)=(xai(ω)/Faj,Ext→A(ω))∣Fak,Ext→A=0,∀k≠j refers to the spatial receptance vector, whose output is the translation of point ai and whose input is the force at point aj, when all other forces applied to substructure A are zero. The receptance of the assembled structure ABij can be expressed according to the receptance of substructures Aij and Bij and the receptance of interface H. The movement equations of each end point of substructures A and B and of structure AB are written as shown in (6) to (8). The behaviour of the interface, assuming zero mass, is given in (9), and the compatibility conditions are written in  (10).SubstructureA:(6)xa1=A11Fa1,Ext⁡→A+A12Fa2,H→A,xa2=A21Fa1,Ext⁡→A+A22Fa2,H→A.SubstructureB:
(7)xb1=B11Fb1,H→B+B12Fb2,Ext⁡→B,xb2=B21Fb1,H→B+B22Fb2,Ext⁡→B.SubstructureAB:(8)xab1=AB11Fab1,Ext⁡→AB+AB12Fab2,Ext⁡→AB,xab2=AB21Fab1,Ext⁡→AB+AB22Fab2,Ext⁡→AB.Interface:(9)xh2-xh1=HFh2,B→H,Fh2,B→H=-Fh1,A→H.
Compatibility conditions:(10)xa1=xab1,xb2=xab2,xh1=xa2,xh2=xb1,Fa2,H→A=Fh1,A→H,Fb1,H→B=Fh2,B→H,Fa1,Ext⁡→A=Fab1,Ext⁡→AB,Fb2,Ext⁡→B=Fab2,Ext⁡→AB.
By integrating the compatibility and interface equations ((9) and (10)) into (6) to (8), system (11) is obtained(11)Fab1,Ext⁡→AB(A11-AB11)+Fa2,H→A(A12)+Fab2,Ext⁡→AB(-AB12)=0,Fab1,Ext⁡→AB(A21)+Fa2,H→A(A22+B11)+Fab2,Ext⁡→AB(-B12)=0,Fab1,Ext⁡→AB(-AB21)+Fa2,H→A(-B21)+Fab2,Ext⁡→AB(B22-AB22)=0,AB11 and AB21 are elaborated from the definition of ABij(ω)=(xabi(ω)/Fabj,Ext→AB(ω))∣Fabk,Ext→AB=0,∀k≠j, when Fab2,Ext→AB=0 in (11). The receptance equations are obtained ((12)-(13))
(12)AB11=A11-A12[B11+A22+H]-1A21,(13)AB21=A21[B11+A22+H]-1B21.
Similarly, to obtain AB22 and, AB12 it is sufficient to impose Fab1,Ext→AB=0 in (11). Then the standard receptance equations are obtained(14)(14)AB22=B22-B12[B11+A22+H]-1B21,AB12=A12[B11+A22+H]-1B12.From the receptance of the system substructuresAij, Bij, associated with the interface model H, the reconstruction of the assembled structure FRF: ABij is possible.


#### 2.4.2. Application to the SVDH-Vibrating Subsystem/SVDH Body Interface

The interface between the SVDH-vibrating subsystem and the SVDH body is a prismatic joint whose axial stiffness is controlled by a spring. The SVDH-vibrating subsystem is guided in axial translation through a ball retainer, which controls the radial stiffness of the interface. Figure5(a) shows these components. The axial dynamic behaviour of the interface was modelled using a spring damper: Hmodel=1/(k+icω). This interface was assumed to be rigid in the radial direction.(a) Components (SVDH-body, SVDH vibrating subsystem) and interfaces (spring and ball retainer), (b) identification of the stiffness:k and the damping factor: c by minimizing the gap between AB11-reconstructed and AB11-direct.
(a)(b)Thek and c values were determined by minimizing the following r criteria:
(15)r=∑ω=ω0-εω0+ε(AB11-reconstructed(ω,k,c)-AB11-direct(ω))2,
where AB11-reconstructed and AB11-direct represent, respectively, the RC-constructed FRF obtained by (12) and the measured FRF on the assembled system. ω0 is the interface mode pulsation. The optimisation procedure was carried out on a 3 dB bandwidth around ω0. In Figure 5(b), it can be noticed that the AB11-reconstructedand AB11-direct curves are in good agreement, which enables the spring-damper interface model and the identified values (k=1.13×106 N/m and c=24 Ns/m) to be validated.


#### 2.4.3. Application to the Collet Chuck Interface

The interface between the drill and the SVDH-vibrating subsystem is a collet chuck joint. ReceptanceAB21-direct is preferred to AB11-direct in order to facilitate the experimental identification procedure. Moreover, the mass of the accelerometer is not negligible compared to the mass of the drill. In addition, the SVDH-vibrating subsystem is not in the same material configuration when the collet chuck is tightened on a drill as when it is not. For these reasons, receptance AB21-reconstructed (13) was obtained using numerical receptance for components A and B. The identification procedure provided collet chuck stiffness and damping factors of, respectively, 14.8 × 106 N/m and 4 Ns/m.


#### 2.4.4. Application to the HSK 63 Taper Interface

The experimental modal analyses carried out on the spindle/SVDH body system in the axial direction allow the HSK 63 interface to be considered as a rigid connection (Figure6). Indeed, no specific mode for the HSK 63 taper interface appears between 0 and 3000 Hz.HSK 63 experimental modal analyses and associated axial FRF.
(a)(b)(c)


### 2.5. Model Assembly and Experimental Validation

As in a classic finite element procedure, dynamic equations of the overall system, composed of the drill, the SVDH-vibrating subsystems, the SVDH body, and the spindle, were obtained by assembling element matrices. The spring-damper connection parameters between the drill and the SVDH-vibrating subsystems and between the SVDH-vibrating subsystems and SVDH body, identified by the receptance coupling method, enabled the rotor-beam models of the components to be assembled.The spindle-SVDH-tool assembled model was validated by comparison between numerical and experimental FRF, as shown in Figure7.Numerical versus experimental system FRF.
(a)(b)Figure7(a) represents experimental and numerical axial FRF of the assembled system. The 60 Hz and 4700 Hz modes are, respectively, due to the spring-ball retainer and collet chuck interfaces as mentioned in the previous subsystem identification procedure. Some parasitical experimental bending modes at 193 Hz, 1237 Hz, and 3433 Hz are present in the experimental FRF.Figure7(b) represents experimental and numerical radial FRF of the assembled system. The 193, 1237, and 3433 Hz modes are related to the drill’s bending modes. The 376 Hz mode is controlled by the bearings. Additional numerical modes are present at 1830 and 3740 Hz. These frequency peaks are related to the dynamics behaviour of the rear side of the spindle. They do not appear in the experimental FRF since the displacements of the rear side of the real spindle are blocked by the motor. For the axial and radial FRF, a good correspondence between the numerical and experimental curves enabled the numerical model to be used for further investigations.


## 2.1. Modelling of Structural Subsystems

The spindle SVDH tool system is composed of four structural subsystems: the drill, the SVDH-vibrating subsystem, the SVDH body, and the spindle. The motion of the rotating structure is considered as the superposition of rigid and elastic body displacements. Dynamic equations were obtained using Lagrange formulation associated with a finite element method. Due to the size of the rotor sections, shear deformations had to be taken into account. Then, the rotating substructure was derived using Timoshenko beam theory. The relevant shape functions were cubic in order to avoid shear locking. A special three-dimensional rotor-beam element with two nodes and six degrees of freedom per node was developed in the corotational reference frame. The damping model used draws onRayleigh viscous equivalent damping, which makes it possible to regard the damping matrix D as a linear combination of the mass matrix M and the spindle rigidity matrix K(2)D=aM+bK,
where a and b are damping coefficients.The set of differential equations can be written as(3)M(qN)q̈N+(C(qN,q̇N)+D)q̇N+KqN=F(t),
where M is the mass matrix, and Cmatrix contains the rotational dynamics effects. qN and F(t) are the nodal displacement and force vectors. An accelerating rotor gives rise to previous time-variant equations (3), but treatment of the rotor using a pseudoconstant speed approach can still be described by means of linear time-invariant models and is valid in many cases (4):
(4)M0q̈N+(2ΩG+D)q̇N+(K-Ω2N)qN=F(t),
where M0 is the constant part of matrix M, and G and N proceed from the decomposition of matrix C.


## 2.2. Modelling Angular Contact Ball Bearings

The rotating system is supported by four (two front and two rear) hybrid angular contact bearings. The rolling bearing stiffness matrices were calculated using in-house software developed on the basis of Lim and Singh’s [20] formulation. The bearing stiffness model represents the load-displacement relation combined with the Hertzian contact stress principle and was calculated around a static function point characterised by the bearing preload: δ. Based on Rantatalo’s prediction [9], the initially calculated bearing stiffness is spindle speed dependent because of the gyroscopic and centrifugal force Fc, which acts on each ball (Figure 3(a)). As the speed increased, the load conditions between the balls and the rings in the bearing changed because of the centrifugal force (Figure 3(a)). Then, speed-dependent bearing stiffness was integrated into the global spindle FEM and influenced the natural frequencies of the spindle-SVDH-tool unit under consideration.Bearing stiffness variation depending on spindle speed [12].
(a)
Angular contact ball bearing Hertzian and centrifugal force(b)
Bearing stiffness variation


## 2.3. Structural Systems Dynamic Readjustment

While detailed knowledge of the spindle-SVDH-tool system is in general not available in a manufacturing environment, models need to be readjusted in order to fit experimental results. The readjustment parameters are the Young modulus:E, the damping coefficients: (a,b), and the rolling bearing preload: δ. The readjustment procedure is proposed to tune the previous variable in order to fit the model results to the experimental frequency response function (FRF). These parameters are readjusted by minimizing the gap between the measured and the modelled tool tip node FRF for nonrotating components, using an optimisation routine and a least squares type objective function defined as(5)r=∑ω=ωmin⁡ωmax⁡(Hnum(ω,(E,a,b,δ))-Hexp⁡(ω))2,
where Hnum(ω,(E,a,b,δ)) and Hexp(ω) are, respectively, the numerical and experimental FRF.


## 2.4. Modelling and Identification of Spindle-SVDH-Tool Interfaces

The dynamic behaviour of the interfaces represented by the HSK63 taper, spring and ball retainer, and collet chuck was taken into account. The identification procedure of the interface models was based on the receptance coupling method.


### 2.4.1. Receptance Coupling Background

In this section, the RC equations are established. At the top of Figure4, the substructure A and the substructure B are represented, connected by interface H. At the bottom of Figure 4, the assembled structure is represented. This figure enables the excitation point and the measurement points used in the receptance coupling approach to be located.Figure 4
Receptance coupling notation.The notationAij(ω)=(xai(ω)/Faj,Ext→A(ω))∣Fak,Ext→A=0,∀k≠j refers to the spatial receptance vector, whose output is the translation of point ai and whose input is the force at point aj, when all other forces applied to substructure A are zero. The receptance of the assembled structure ABij can be expressed according to the receptance of substructures Aij and Bij and the receptance of interface H. The movement equations of each end point of substructures A and B and of structure AB are written as shown in (6) to (8). The behaviour of the interface, assuming zero mass, is given in (9), and the compatibility conditions are written in  (10).SubstructureA:(6)xa1=A11Fa1,Ext⁡→A+A12Fa2,H→A,xa2=A21Fa1,Ext⁡→A+A22Fa2,H→A.SubstructureB:
(7)xb1=B11Fb1,H→B+B12Fb2,Ext⁡→B,xb2=B21Fb1,H→B+B22Fb2,Ext⁡→B.SubstructureAB:(8)xab1=AB11Fab1,Ext⁡→AB+AB12Fab2,Ext⁡→AB,xab2=AB21Fab1,Ext⁡→AB+AB22Fab2,Ext⁡→AB.Interface:(9)xh2-xh1=HFh2,B→H,Fh2,B→H=-Fh1,A→H.
Compatibility conditions:(10)xa1=xab1,xb2=xab2,xh1=xa2,xh2=xb1,Fa2,H→A=Fh1,A→H,Fb1,H→B=Fh2,B→H,Fa1,Ext⁡→A=Fab1,Ext⁡→AB,Fb2,Ext⁡→B=Fab2,Ext⁡→AB.
By integrating the compatibility and interface equations ((9) and (10)) into (6) to (8), system (11) is obtained(11)Fab1,Ext⁡→AB(A11-AB11)+Fa2,H→A(A12)+Fab2,Ext⁡→AB(-AB12)=0,Fab1,Ext⁡→AB(A21)+Fa2,H→A(A22+B11)+Fab2,Ext⁡→AB(-B12)=0,Fab1,Ext⁡→AB(-AB21)+Fa2,H→A(-B21)+Fab2,Ext⁡→AB(B22-AB22)=0,AB11 and AB21 are elaborated from the definition of ABij(ω)=(xabi(ω)/Fabj,Ext→AB(ω))∣Fabk,Ext→AB=0,∀k≠j, when Fab2,Ext→AB=0 in (11). The receptance equations are obtained ((12)-(13))
(12)AB11=A11-A12[B11+A22+H]-1A21,(13)AB21=A21[B11+A22+H]-1B21.
Similarly, to obtain AB22 and, AB12 it is sufficient to impose Fab1,Ext→AB=0 in (11). Then the standard receptance equations are obtained(14)(14)AB22=B22-B12[B11+A22+H]-1B21,AB12=A12[B11+A22+H]-1B12.From the receptance of the system substructuresAij, Bij, associated with the interface model H, the reconstruction of the assembled structure FRF: ABij is possible.


### 2.4.2. Application to the SVDH-Vibrating Subsystem/SVDH Body Interface

The interface between the SVDH-vibrating subsystem and the SVDH body is a prismatic joint whose axial stiffness is controlled by a spring. The SVDH-vibrating subsystem is guided in axial translation through a ball retainer, which controls the radial stiffness of the interface. Figure5(a) shows these components. The axial dynamic behaviour of the interface was modelled using a spring damper: Hmodel=1/(k+icω). This interface was assumed to be rigid in the radial direction.(a) Components (SVDH-body, SVDH vibrating subsystem) and interfaces (spring and ball retainer), (b) identification of the stiffness:k and the damping factor: c by minimizing the gap between AB11-reconstructed and AB11-direct.
(a)(b)Thek and c values were determined by minimizing the following r criteria:
(15)r=∑ω=ω0-εω0+ε(AB11-reconstructed(ω,k,c)-AB11-direct(ω))2,
where AB11-reconstructed and AB11-direct represent, respectively, the RC-constructed FRF obtained by (12) and the measured FRF on the assembled system. ω0 is the interface mode pulsation. The optimisation procedure was carried out on a 3 dB bandwidth around ω0. In Figure 5(b), it can be noticed that the AB11-reconstructedand AB11-direct curves are in good agreement, which enables the spring-damper interface model and the identified values (k=1.13×106 N/m and c=24 Ns/m) to be validated.


### 2.4.3. Application to the Collet Chuck Interface

The interface between the drill and the SVDH-vibrating subsystem is a collet chuck joint. ReceptanceAB21-direct is preferred to AB11-direct in order to facilitate the experimental identification procedure. Moreover, the mass of the accelerometer is not negligible compared to the mass of the drill. In addition, the SVDH-vibrating subsystem is not in the same material configuration when the collet chuck is tightened on a drill as when it is not. For these reasons, receptance AB21-reconstructed (13) was obtained using numerical receptance for components A and B. The identification procedure provided collet chuck stiffness and damping factors of, respectively, 14.8 × 106 N/m and 4 Ns/m.


### 2.4.4. Application to the HSK 63 Taper Interface

The experimental modal analyses carried out on the spindle/SVDH body system in the axial direction allow the HSK 63 interface to be considered as a rigid connection (Figure6). Indeed, no specific mode for the HSK 63 taper interface appears between 0 and 3000 Hz.HSK 63 experimental modal analyses and associated axial FRF.
(a)(b)(c)


## 2.4.1. Receptance Coupling Background

In this section, the RC equations are established. At the top of Figure4, the substructure A and the substructure B are represented, connected by interface H. At the bottom of Figure 4, the assembled structure is represented. This figure enables the excitation point and the measurement points used in the receptance coupling approach to be located.Figure 4
Receptance coupling notation.The notationAij(ω)=(xai(ω)/Faj,Ext→A(ω))∣Fak,Ext→A=0,∀k≠j refers to the spatial receptance vector, whose output is the translation of point ai and whose input is the force at point aj, when all other forces applied to substructure A are zero. The receptance of the assembled structure ABij can be expressed according to the receptance of substructures Aij and Bij and the receptance of interface H. The movement equations of each end point of substructures A and B and of structure AB are written as shown in (6) to (8). The behaviour of the interface, assuming zero mass, is given in (9), and the compatibility conditions are written in  (10).SubstructureA:(6)xa1=A11Fa1,Ext⁡→A+A12Fa2,H→A,xa2=A21Fa1,Ext⁡→A+A22Fa2,H→A.SubstructureB:
(7)xb1=B11Fb1,H→B+B12Fb2,Ext⁡→B,xb2=B21Fb1,H→B+B22Fb2,Ext⁡→B.SubstructureAB:(8)xab1=AB11Fab1,Ext⁡→AB+AB12Fab2,Ext⁡→AB,xab2=AB21Fab1,Ext⁡→AB+AB22Fab2,Ext⁡→AB.Interface:(9)xh2-xh1=HFh2,B→H,Fh2,B→H=-Fh1,A→H.
Compatibility conditions:(10)xa1=xab1,xb2=xab2,xh1=xa2,xh2=xb1,Fa2,H→A=Fh1,A→H,Fb1,H→B=Fh2,B→H,Fa1,Ext⁡→A=Fab1,Ext⁡→AB,Fb2,Ext⁡→B=Fab2,Ext⁡→AB.
By integrating the compatibility and interface equations ((9) and (10)) into (6) to (8), system (11) is obtained(11)Fab1,Ext⁡→AB(A11-AB11)+Fa2,H→A(A12)+Fab2,Ext⁡→AB(-AB12)=0,Fab1,Ext⁡→AB(A21)+Fa2,H→A(A22+B11)+Fab2,Ext⁡→AB(-B12)=0,Fab1,Ext⁡→AB(-AB21)+Fa2,H→A(-B21)+Fab2,Ext⁡→AB(B22-AB22)=0,AB11 and AB21 are elaborated from the definition of ABij(ω)=(xabi(ω)/Fabj,Ext→AB(ω))∣Fabk,Ext→AB=0,∀k≠j, when Fab2,Ext→AB=0 in (11). The receptance equations are obtained ((12)-(13))
(12)AB11=A11-A12[B11+A22+H]-1A21,(13)AB21=A21[B11+A22+H]-1B21.
Similarly, to obtain AB22 and, AB12 it is sufficient to impose Fab1,Ext→AB=0 in (11). Then the standard receptance equations are obtained(14)(14)AB22=B22-B12[B11+A22+H]-1B21,AB12=A12[B11+A22+H]-1B12.From the receptance of the system substructuresAij, Bij, associated with the interface model H, the reconstruction of the assembled structure FRF: ABij is possible.


## 2.4.2. Application to the SVDH-Vibrating Subsystem/SVDH Body Interface

The interface between the SVDH-vibrating subsystem and the SVDH body is a prismatic joint whose axial stiffness is controlled by a spring. The SVDH-vibrating subsystem is guided in axial translation through a ball retainer, which controls the radial stiffness of the interface. Figure5(a) shows these components. The axial dynamic behaviour of the interface was modelled using a spring damper: Hmodel=1/(k+icω). This interface was assumed to be rigid in the radial direction.(a) Components (SVDH-body, SVDH vibrating subsystem) and interfaces (spring and ball retainer), (b) identification of the stiffness:k and the damping factor: c by minimizing the gap between AB11-reconstructed and AB11-direct.
(a)(b)Thek and c values were determined by minimizing the following r criteria:
(15)r=∑ω=ω0-εω0+ε(AB11-reconstructed(ω,k,c)-AB11-direct(ω))2,
where AB11-reconstructed and AB11-direct represent, respectively, the RC-constructed FRF obtained by (12) and the measured FRF on the assembled system. ω0 is the interface mode pulsation. The optimisation procedure was carried out on a 3 dB bandwidth around ω0. In Figure 5(b), it can be noticed that the AB11-reconstructedand AB11-direct curves are in good agreement, which enables the spring-damper interface model and the identified values (k=1.13×106 N/m and c=24 Ns/m) to be validated.


## 2.4.3. Application to the Collet Chuck Interface

The interface between the drill and the SVDH-vibrating subsystem is a collet chuck joint. ReceptanceAB21-direct is preferred to AB11-direct in order to facilitate the experimental identification procedure. Moreover, the mass of the accelerometer is not negligible compared to the mass of the drill. In addition, the SVDH-vibrating subsystem is not in the same material configuration when the collet chuck is tightened on a drill as when it is not. For these reasons, receptance AB21-reconstructed (13) was obtained using numerical receptance for components A and B. The identification procedure provided collet chuck stiffness and damping factors of, respectively, 14.8 × 106 N/m and 4 Ns/m.


## 2.4.4. Application to the HSK 63 Taper Interface

The experimental modal analyses carried out on the spindle/SVDH body system in the axial direction allow the HSK 63 interface to be considered as a rigid connection (Figure6). Indeed, no specific mode for the HSK 63 taper interface appears between 0 and 3000 Hz.HSK 63 experimental modal analyses and associated axial FRF.
(a)(b)(c)


## 2.5. Model Assembly and Experimental Validation

As in a classic finite element procedure, dynamic equations of the overall system, composed of the drill, the SVDH-vibrating subsystems, the SVDH body, and the spindle, were obtained by assembling element matrices. The spring-damper connection parameters between the drill and the SVDH-vibrating subsystems and between the SVDH-vibrating subsystems and SVDH body, identified by the receptance coupling method, enabled the rotor-beam models of the components to be assembled.The spindle-SVDH-tool assembled model was validated by comparison between numerical and experimental FRF, as shown in Figure7.Numerical versus experimental system FRF.
(a)(b)Figure7(a) represents experimental and numerical axial FRF of the assembled system. The 60 Hz and 4700 Hz modes are, respectively, due to the spring-ball retainer and collet chuck interfaces as mentioned in the previous subsystem identification procedure. Some parasitical experimental bending modes at 193 Hz, 1237 Hz, and 3433 Hz are present in the experimental FRF.Figure7(b) represents experimental and numerical radial FRF of the assembled system. The 193, 1237, and 3433 Hz modes are related to the drill’s bending modes. The 376 Hz mode is controlled by the bearings. Additional numerical modes are present at 1830 and 3740 Hz. These frequency peaks are related to the dynamics behaviour of the rear side of the spindle. They do not appear in the experimental FRF since the displacements of the rear side of the real spindle are blocked by the motor. For the axial and radial FRF, a good correspondence between the numerical and experimental curves enabled the numerical model to be used for further investigations.


## 3. Bearing Lifespan Predictions

The industrial context of the proposed paper is to realize high-speed vibratory drilling operation, with a drill diameter of 5 mm, with a drill depth of 100 mm, without retreat cycles or lubrication, in a mass production system.In this section, the numerical model will be used to predict bearings lifespan for various cutting conditions, in this context. The calculated bearing lifespan, thanks to experimental data, was compared to industrial recommendations in order to give rules of uses to obtain the maximum removal rate with respect to system lifespan.For bearing lifespan calculations, experiments were carried out to measure the cutting force for different cutting conditions, using a three-component dynamometer (Kistler dynamometer type 9257B). High-speed vibratory drilling operations, representative of the industrial context, were performed on 35MnV7 steel, with a drill of 116 mm length and of 5 mm diameter. The experimental setup is represented in Figure8(a). Rotation speed and feed rate were tested, respectively, between 5500 and 15500 rpm and between 0.05 and 0.15 mm/rev. Only the drilling operations that led to regular chips fragmentation were used to predict bearing lifespan, as shown in Figure 8(b). The time-variant value of the radial forces, which depends on the angular orientation of the cutter as it rotates through the cut, was expanded into a Fourier series and then truncated to include only the fundamental frequency. These forces, combined with the numerical FRF of the model, between the excitation node and the bearings nodes, allow the resulting rolling bearings solicitation to be determined. The rolling bearings lifespan is
(16)Lh=10660N(CP)3,
where C is the basic dynamic load and P the equivalent dynamic load in Newtons. N is the spindle speed in rpm. The values obtained are expressed in hours and can be compared to an industrial objective of 20000 hours.Experimental setup.
(a)(b)


### 3.1. Effect of Cutting Conditions on Bearing Lifespan

First, the bearings with the shortest lifespan is sought, because it is the element upon which the lifespan of the entire spindle depends. For all cutting conditions, the bearings with the shortest lifespan is always the front bearings. From Figure9, which represents the extreme lifespan of the different bearings, the strong influence of the cutting conditions on the bearings lifespan can be noticed. Industrial recommendations for spindle lifespan are 20000 hours. The shortest calculated lifespan is 55000 and corresponds to the cutting conditions given rise to the best material removal rate (N=12000 rpm, f=0.125 mm/rev). Thus, a high-speed vibratory drilling operation is always compatible with industrial standards, even under the best material removal rate cutting conditions.Figure 9
Influence of location and cutting conditions on bearings lifespan.Figure10(a) shows the influence of the feed rate on bearings lifespan, for a spindle speed of 12000 rpm. Only feed rates higher than 0.05 mm/rev were tested, because for feed rates below this limit, the productivity is too low and high-speed vibratory drilling loses its relevance. For feed rates higher than 0.125 mm/rev, no stable vibrations of the drill were obtained. The bigger the feed rate, the lower the bearings Lifespan. However, even the shortest bearings Lifespan is compatible with industrial standard. Hence, the optimal feed rate of 0.125 mm/rev is retained. Figure 10(b) shows that cutting forces are modified by the feed rate. Similarities between the curves in Figures 10(a) and 10(b) indicate that the bearing lifespan is mainly influenced by the feed rate.(a) Influence of feed rate on bearings lifespan, for a spindle speed of 12000 rpm, and (b) influence of the feed rate on the inverse of experimental axial cutting force magnitude.
(a)(b)Figure11(b) shows the variations in bearing lifespan depending on spindle speed, for a feed rate of 0.05 mm per revolution. The curve enables the determination of a spindle speed at 8500 rpm which maximizes bearing Lifespan and a spindle speed at 11500 rpm which minimizes it. In order to optimize productivity, as shown on Figure 8(b), the recommended cutting spindle speed is 12000 rpm.Figure 11
(a) Campbell diagram. (b) Influence of spindle speed on bearings lifespan, for a feed rate of 0.05 mm/rev.The curves of Figure11 illustrate the dynamic effects due to high rotational speed, such as gyroscopic coupling and spin softening, on system behaviour, and hence on bearing lifespan. Figure 11(a) represents the effect of spindle speed on the first two radial modes of the system. The dotted line indicates the drill oscillation frequency and is plotted from experimental data measured during high-speed vibratory drilling operations. The critical speed, at 11500 rpm, corresponds to the intersection of the radial mode frequencies with the excitation line due to the drill oscillations.


## 3.1. Effect of Cutting Conditions on Bearing Lifespan

First, the bearings with the shortest lifespan is sought, because it is the element upon which the lifespan of the entire spindle depends. For all cutting conditions, the bearings with the shortest lifespan is always the front bearings. From Figure9, which represents the extreme lifespan of the different bearings, the strong influence of the cutting conditions on the bearings lifespan can be noticed. Industrial recommendations for spindle lifespan are 20000 hours. The shortest calculated lifespan is 55000 and corresponds to the cutting conditions given rise to the best material removal rate (N=12000 rpm, f=0.125 mm/rev). Thus, a high-speed vibratory drilling operation is always compatible with industrial standards, even under the best material removal rate cutting conditions.Figure 9
Influence of location and cutting conditions on bearings lifespan.Figure10(a) shows the influence of the feed rate on bearings lifespan, for a spindle speed of 12000 rpm. Only feed rates higher than 0.05 mm/rev were tested, because for feed rates below this limit, the productivity is too low and high-speed vibratory drilling loses its relevance. For feed rates higher than 0.125 mm/rev, no stable vibrations of the drill were obtained. The bigger the feed rate, the lower the bearings Lifespan. However, even the shortest bearings Lifespan is compatible with industrial standard. Hence, the optimal feed rate of 0.125 mm/rev is retained. Figure 10(b) shows that cutting forces are modified by the feed rate. Similarities between the curves in Figures 10(a) and 10(b) indicate that the bearing lifespan is mainly influenced by the feed rate.(a) Influence of feed rate on bearings lifespan, for a spindle speed of 12000 rpm, and (b) influence of the feed rate on the inverse of experimental axial cutting force magnitude.
(a)(b)Figure11(b) shows the variations in bearing lifespan depending on spindle speed, for a feed rate of 0.05 mm per revolution. The curve enables the determination of a spindle speed at 8500 rpm which maximizes bearing Lifespan and a spindle speed at 11500 rpm which minimizes it. In order to optimize productivity, as shown on Figure 8(b), the recommended cutting spindle speed is 12000 rpm.Figure 11
(a) Campbell diagram. (b) Influence of spindle speed on bearings lifespan, for a feed rate of 0.05 mm/rev.The curves of Figure11 illustrate the dynamic effects due to high rotational speed, such as gyroscopic coupling and spin softening, on system behaviour, and hence on bearing lifespan. Figure 11(a) represents the effect of spindle speed on the first two radial modes of the system. The dotted line indicates the drill oscillation frequency and is plotted from experimental data measured during high-speed vibratory drilling operations. The critical speed, at 11500 rpm, corresponds to the intersection of the radial mode frequencies with the excitation line due to the drill oscillations.


## 4. Conclusions

In this paper, a comprehensive approach to developing a hybrid model of the dynamic behavior of the spindle self-vibratory drilling head—tool system has been proposed. This approach has resulted in a numerical model enriched with physical data. The various components of the system are modelled using a specific beam element, taking into account the gyroscopic effects, centrifugal forces, and shear deflection. The receptance coupling method is used to identify the dynamic behavior of the interface. The complete system is then obtained by assembling the beam model of each component using spring-damper elements. Finally, the assembled model is validated by comparisons between numerical and experimental FRFs. The model is used to predict the influence of a high-speed vibratory drilling operation on the bearing lifespan. The predictions of bearing lifespan are used to give rules of uses of the high-speed vibratory drilling head.



---

*Source: 606087-2011-08-21.xml*