Shock and Vibration
Volume 2020 (2020), Article ID 8816773, 12 pages
https://doi.org/10.1155/2020/8816773

Research Article

Dynamic Modeling and Cutting Stability of Rotating Tapered Composite Cutter Bar considering Material Damping

Yuhuan Zhang, Ren Yongsheng, Bole Ma, and Jinfeng Zhang

College of Mechanical and Electronic Engineering, Shandong University of Science and Technology, Qindao 266510, China

Correspondence should be addressed to Ren Yongsheng; renys@sdust.edu.cn

Received 16 June 2020; Revised 27 July 2020; Accepted 7 August 2020; Published 27 August 2020

Academic Editor: Sundararajan Natarajan

Copyright © 2020 Yuhuan Zhang et al. 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.

Traditional milling cutter bars are generally made up of metals and exhibit poor capacity of chatter suppression. This study proposes an anisotropic composites tapered cutter bar for increasing natural frequency and damping and finally achieves the goal of enhancing chatter stability. Based on Hamilton principle and Euler–Bernoulli beam theory, the partial differential motion equations of the cutting system with a 3D rotating tapered composite cutter bar are established. Next, using the Galerkin method, the equations of motion are discretized so as to derive ordinary differential equations. In the model, damping modeling of the composite cutter bar is achieved theoretically by using damping dissipation constitutive relations for viscoelastic composites. Moreover, by introducing the rotating effect of the 3D cutter bar in the 2-DOF analytical model of stability analysis first proposed for a fixed-type cutter bar, an improved prediction model is developed and used to solve the stability lobes of the cutting system in the frequency domain analytically. Furthermore, the influences of the gyroscopic effect, material, ply angle, stacking sequence, and taper ratio on chatter stability are also discussed.

1. Introduction

Due to low dynamic stiffness and natural frequency, traditional metal cutter tools can hardly perform high-speed boring operations or end-milling operations under chatter. Chatter can be roughly classified into two types, regenerative chatter and mode coupled chatter. Regenerative chatter refers to self-excited vibration under the interaction between cutting force and undulation on the workpiece surface after the cutting of the previous cutter tooth, while mode coupled chatter is induced by the coupling between two natural mode shapes since the vibration system exhibits slight discrepancy in stiffness along two directions. Chatter can inevitably reduce the workpiece’s machining precision and surface processing quality and simultaneously shorten service lives of both the cutting tool and machine tool. This study focuses on regenerative chatter on account of its importance in the cutting process.

Previous research reveals that chatter suppression capability of a cutting system is tightly connected with its damping and static stiffness. The boring bar with higher damping and static stiffness also has a greater dynamic stiffness, and the boring bar with a greater specific stiffness also exhibits a larger natural frequency. Traditional materials cannot ensure great dynamic stiffness and natural frequency at the same time. Owing to a series of advantages including high modulus, great damping and low density, advanced fiber reinforced composites seem to be the best choice for replacing metal materials and developing the boring bars with excellent chatter suppression stability. Fiber reinforced composites now have been widely applied in the structural dynamic design of rotating structures [1–7], such as composite automotive drive shafts, high-speed spindles, and rollers. As a result, the modeling theory and numerical method of rotating composite shaft were developed [8–15]. Meantime, a composite is used in the design of leaf spring for the lightweight of automobiles [16]. Scholars also paid increasing attentions to structural design and dynamics of composites boring bar [17–20]. Lee and Suh [17] performed vibration tests and cutting tests and found that the fundamental frequency and damping ratio of graphite/epoxy composite boring bar are 80% and 50% higher than those of steel boring bar, and the maximum cutting depth using the graphite/epoxy composite boring bar is 5 times greater than that using the steel boring bar. Nagano et al. [19] employed finite element analysis and cutting tests for investigating dynamic properties and chatter stability of the carbon/epoxy composite boring bar with different steel core shapes. Ghorbani et al. [20] utilized an experimental method to study the vibration suppression performance of boring bars filled with epoxy granite, with different cross-sections. The vibro-acoustic signals were measured during a machining process. Results showed that application of this new design of boring bars in the machining process can improve the surface quality up to 30% compared to the conventional boring bar.

Nevertheless, abovementioned studies only examined the simple cases when the boring bar is fixed and the workpiece is rotating. In practical applications, the boring bar is mostly rotates at a high speed while the workpiece is fixed. Lee et al. [18] performed the vibration test and cutting test and optimized the dynamic properties of the carbon/epoxy boring bar with different sizes of cores made up of different materials. Without chatter, the cutting depth can be enhanced by 30% compared with that of the tungsten carbide boring bar. It should be noted that abovementioned studies were performed using a finite element analysis or test. Based on Timoshenko beam theory, Kim et al. [21] established the model of composite cantilever end-milling bar, as well as proposed the chatter analysis model of composite end-milling bar by assuming that regenerative delay cutting force acted on the free end. The research also demonstrates that composite boring bar can greatly improve cutting stability.

It is well known that composite material has higher damping than conventional metallic material. However, it has been shown in previous research [22] that internal damping in a rotating system may lead to whirl instabilities in high-speed rotors. Therefore, accurate prediction of effects of internal damping in composite material rotors, particularly effects of composite damping on the stability of a rotating cutter bar, is essential.

Even though the internal viscous damping constant was included in Kim et al. [21] to investigate effects of composite damping on the cutting stability and the damping ratio of the composite damping cutter bar is given empirically, no composite damping modeling has been performed yet. Ren et al. [22] presented a damping analysis model for the prediction of the variation of the natural frequency and damping with rotating speed. The critical rotating speeds and instability thresholds of composite shaft are evaluated. However, no cutting stability of the composite cutting bar was involved in this study.

Altintas and Budak [23] presented the 2-DOF frequency-domain solutions for the milling process with stationary cutter bars and rotating workpieces in which the machining force directions on the cutter bars in the inertial coordinates are fixed. The cutting stability of the process with the high-speed rotating cutter bar can be conducted in the cutter bar’s rotating coordinate frame, in which case the directional force coefficients are time invariant. However, the two orthogonal coordinates of the tool become dynamically coupled as a function of spindle speed.

A new cutting dynamic model with a rotating tapered composite cutter bar is developed in this paper. The partial differential motion equations of the cutting system with the 3D continuously distributed parameters cutter bar which incorporated internal damping and rotating effects are derived based on Hamilton principle and Euler–Bernoulli beam theory. After discretization using the Galerkin method, ordinary differential equations were obtained. Here, inclusion of viscous damping has been achieved theoretically by using the damping analysis model for composite shaft considering damping dissipation [22]. At the same time, an improving analytical prediction of stability lobes for the cutting system with a 3D cutter bar has been proposed by introducing the rotating effects of the cutter bar in the 2-DOF model of Comak et al. [24] and Eynian and Altintas [25] based on the rotating frame approach and used to investigate the process stability of the cutting system.

By means of Laplace transformation, a time-domain chatter equation was converted into frequency-domain characteristic equation and the stability lobes of the cutting system were plotted. Moreover, the influences of the gyroscopic effect, material, ply angle, and taper ratio on the chatter stability are investigated.

2. Theoretical Formula

2.1. Differential Equation of Motion

As shown in Figure 1, the radius of cross-section of the rotating tapered composite cutter bar with rotating speed of and length of L changes linearly along the axial length, i.e., , where σ denotes the taper of the cutter bar and and denote the cross-section radius on the fixed and free ends of the cutter bar, respectively.

(a)

(b)

Figure 1: Schematic diagram of the rotating tapered composite cutter bar. (a) Schematic diagram of the cutter bar. (b) Cylindrical coordinate system.

The motion differential equation of the cutter bar is established based on the following Hamilton principle:

where T and U denote the variations of kinetic energy and potential energy, respectively, and denotes the virtual work done by the internal damping force.

Based on the Euler–Bernoulli beam theory, the strain-displacement relations can be described as (Ren et al. [22])

where uy and uz denote the displacement of the neutral axis along y-axis and z-axis direction of the points on the reference axis of the cutter bar, while and denote the rotations about y-axis and z-axis, respectively. Here, the x-y-z is a rotating coordinate system attached to the neutral axis of the cutter bar.

In view of the viscoelastic mechanical characteristics of the composite, the stress components can be split into elastic stress and dissipative [22]:

where σ and τ denote normal stress and shear stress, respectively, and the superscripts e and d correspond to elastic and dissipative damping parts. Because the axial displacement ux and torsion angle φ have no effect on the bending deformation, they can be removed from equation (2). In this case, the elastic stress-strain relations can be expressed as

where denotes reduced stiffness. The detailed definitions for the reduced stiffness are given in Appendix.

The damping dissipation stress-strain relations can be written as

where denotes the coefficient of the internal damping of the cutter bar linked to material damping (or call it off-axis stiffness linked to material damping).

The expression for kinetic energy is

The spatial varying flexural stiffness of the composite cutter bar can be expressed as [22]

in which M denotes an M-layered laminate and and are the external and internal radius of layer k in section x (see Figure 1).

According to Zorzi and Nelson [26], the flexural bending moments of the cutter bar may be written as

where I (x) denotes the area moment of inertia.

The first term of equation (8) is the conservative moment and has already been reflected in the potential energy equation (6). The rest of the terms represent the nonconservative moments which are resultant of the internal damping. Thus, the virtual work due to the internal damping can be expressed as

The kinetic energy of a rotating tapered composite cutter bar can be expressed as

where m (x) and Im (x) denote the mass per unit length and the mass moment of inertia, respectively.

The mass per unit length m (x) and the moment of mass inertia Im (x) can be written as

in which is the density of layer k.

By substituting equations (6), (9), and (10) into equation (1) and assuming that the cutting forces Fy and Fz are produced on the free end of the cutter bar along y-axis and z-axis directions, the bending-bending coupling forced vibration equation of the rotating tapered composite cutter bar can be written as

where denotes the Dirac function:

2.2. Galerkin Method

It is assumed that the differential equation (12) has the solutions:

where denotes the mode function and satisfies the boundary condition .

By substituting equation (14) into equation (12), the following expression can be derived with the use of Galerkin method:

where

Equation (15) is the dynamic equation of the cutting system. For analyzing free vibration characteristics in the bending of the cutter bar, it can be first assumed that Fy = Fz = 0.

In order to determine the internal damping coefficient, the following expression can be derived according to the definition of the damping ratio of 1-DOF free vibration systems

The internal damping ratio is related to the modal loss factor :

where denotes the modal loss factor, which can be defined as by using an energy approach [27]:

where is the mode deflection and and are the stiffness and damping matrices of a tapered nonrotating composite cutter bar, respectively, which have the forms

In the previous equations, the elements of damping matrix are given by [22]

in which denotes the off-axis damping dissipation factor of the composite layer. The detailed expressions can be found in Appendix.

2.3. Cutting Force and Stability Analysis

It is assumed that the milling cutter has N teeth that produce vibration displacements along y-axis and z-axis directions. The cutter bar rotates at a rotating speed of (rad/s) as shown in Figure 2.

Figure 2: Cross-section view of the cutting process.

The regenerative dynamic milling force in the rotating coordinate system can be written as Comak et al. [24]

where

in which t denotes the tooth period , b denotes axial cutting depth, and denotes the tangential cutting constant. [A (t)] and [B (t)] denote directional matrix whose coefficients depend on the instantaneous angular immersion of the cutter as , in which is the pitch angle of the cutter with N equally spaced teeth and is the cutter pitch angle and the angle between the local rotating coordinate frame (yj, zj) and principal rotating frame (y, z).

Since [A (t)] and [B (t)] are periodic at tooth passing frequency, they can be expanded into Fourier series. Using the average values and of time varying directional matrices [A (t)] and [B (t)], the frequency-domain solution of stability in rotating coordinates can be obtained analytically.

For two-fluted end mills, when pitch angle is uniform and180 degrees, the dynamic milling force in equation (22) can be expressed as follows:

Average matrix has the form as follows [25]:

where and are start and exit immersion angles of the cutter to and from the cut and is the radial cutting constant.

2.4. Cutting Stability Analysis

It is assumed for simplicity that there is only one dominant vibration mode in each of the two principal axes of the cutter bar. Under the action of regenerative linear cutting force, the chatter equation of the rotating tapered composite cutter bar can be written as

where

where the symbol prime (‘) denotes the first-order derivative with respect to x.

By conducting Laplace transformation on equation (26) and assuming , it can be found that the value of s is connected with system stability.

If , the system is in unstable state.

If , the system is in stable state.

If , the system is in the critical state between the stable and unstable state.

Let , by substituting into the equation after Laplace transformation of equation (26), and the following expression can be written as

in which denotes the chatter frequency.

The transfer function can be written as

Equation (28) can be rewritten as

The system of linear equation (30) has a nontrivial solution if its determinant is zero:

where denotes an identity matrix.

The stability of the cutter bar can be transformed into the characteristic value of the following second-order characteristic equation:

where

Substituting and into equation (33), one has

Since the limit cutting depth is a real number in the cutting process, is a real number, and then the imaginary part of equation (35) must vanish, which will lead to

Finally, the critical axial depth and spindle speed (rpm) can be obtained as follows:

The stability lobes can be calculated as follows:

(1)Calculate the directional matrix in accordance with the milling cutter’s material and geometrical parameters(2)Start a loop incrementing the spindle speed (3)Scan the chatter frequency ωc by using the natural frequency as the reference, and calculate the transfer function(4)Calculate real and imaginary parts according to equation (32), and solve the critical axial cutting depth blim and the corresponding rotating speed (5)If the difference between the rotating speed obtained in the step 4 and input rotating speed in the step 2 satisfies the required precision, plot vs. blim; otherwise, repeat the iterative procedures until the results converge(6)Select the new j and calculate the adjacent lobes

3. Case Study and Discussion

3.1. Natural Frequency and Logarithmic Decrement

By neglecting the cutting forces in equation (15), the eigenvalue can be solved given a rotating speed , whose imaginary part represents the natural frequency and the ratio of imaginary part to the real part defines the logarithmic decrement, respectively. Therefore, the variation curves of both natural frequency and damping with the rotating speed can be plotted.

In order to investigate the variations of the natural frequency and logarithmic decrement with the rotating speed, a case study is then performed on the tapered composite cutter bar with a cross-section thickness of h = 0.01016 m, RT = 0.176 m. The length L can be determined in accordance with the length-to-diameter ratio (L/RT).

Figures 3 and 4 show the variation curves of the first natural frequency and logarithmic decrement of the cutter bar with the rotating speed when the ply angles are set as different values 0°, 30°, and 60°, respectively, with the taper ratio 0.5 and the length-to-diameter ratio 10. The composites have stacking sequence [±θ]₈. Table 1 lists the material property parameters.

Figure 3: Variation of the first natural frequency with rotating speed at different ply angles.

Figure 4: Variation of the first mode logarithmic decrements with rotating speed at different ply angles.

Table 1: Mechanical properties of carbon/epoxy composite material [13].

ρ (kg/m³)
E₁₁ (GPa)
E₂₂ (GP)
G₁₂ (GP)
G₂₃ (GP)
ν₁₂
(%)
(%)
(%)
(%)

1446.2
172.7
7.2
3.76
3.76
0.3
0.45
4.22
7.05
7.05

The numerical results of first natural frequency and logarithmic decrement are shown in Table 2 for an increasing number of mode shape functions. From these tables, it can be seen that to obtain the accurate results of the first natural frequency, no more than three mode shape functions are required.

Table 2: First frequency and modal logarithmic decrement of the tapered composite cutter bar .

Ply angle (degree)
The first natural frequency (Hz)
The first mode logarithmic decrement

N = 1
N = 3
N = 5
N = 1
N = 3
N = 5

0°
441.52
438.27
438.27
0.0095
0.0087
0.0087

30°
337.95
335.51
335.51
0.0143
0.0128
0.0128

60°
146.05
143.24
143.24
0.0585
0.0455
0.0455

As shown in Figures 3 and 4, as the rotating speed increases from 0, both the variation curves of natural frequency and logarithmic decrement split into two curves. In the upper branch which is corresponding to forward whirling, the natural frequency increases with the rotating speed; in the latter branch which is corresponding to backward whirling, natural frequency drops with the rotating speed. The logarithmic decrement in the upper branch is always positive within the whole rotating speed range, suggesting that the forward whirling is stable. Within certain rotating speed range, the logarithmic decrement in the lower branch is positive; as the rotating speed increases to certain value, the logarithmic decrement in the lower branch changes from positive to negative, suggesting that the rotating cutter bar gradually lost motion stability. The value of corresponding to zero logarithmic decrement is called as the instability threshold. It should be noted that, when the external damping equals to 0, the threshold of instability equals to the critical speed.

It can also be observed from Figure 3 that, decreasing of the ply angle increases the natural frequency of the composite cutter bar. This is because the closer the fiber is oriented to 0°, the greater the cutter bar rigidity and the higher the whirl frequencies are (see Table 1).

Figure 4 suggests that the instability thresholds decrease with ply angle. The reason is that the larger the ply angle, the greater the internal damping due to the composite materials is and the more likely the instability appears (see Table 1).

Table 3 shows the first frequency and logarithmic decrement of the tapered cutter bars which are made of carbon fiber/epoxy composite material and steel. The taper ratios 0.75, 0.5, and the length-to-diameter ratio 10 are used. The composite cutter bar has stacking sequence . Table 1 gives the carbon fiber/epoxy composite material property parameters. Table 4 gives the mechanical properties of steel.

Table 3: First frequency and modal logarithmic decrement for the tapered cutter bar made of different materials.

Material
The first natural frequency (Hz)
The first mode logarithmic decrement

Carbon fiber/epoxy
355.82
438.27
0.0089
0.0087

Steel
132.56
180.38
0.0025
0.0024

Table 4: Mechanical properties of steel material.

ρ (kg/m³)
E (GPa)
ν
ς (%)

7850
210
0.3
0.01

It can be observed from Table 3 that the improvement of frequency and logarithmic decrement in using the anisotropic composite tapered cutter bar when compared with an isotropic tapered cutter bar is obvious. For and , the differences in the first natural frequency are 168% and 143%, and the differences in the first mode logarithmic decrement are 256% and 263%.

Figure 5 shows the variation of the modal loss factor of the composite cutter bar with ply angle. It can be seen that the loss factor in the first mode increases with the value of ply angle increasing from 0° to 81°. It is minimum at 0° ply angle, about 0.45%, and has a maximum value, about 4.28%, at 81° ply angle. An insignificant amount of reduction in the modal loss factor is observed in the range of 81°～90° (about 4.22% for 90° ply angle). Such behavior is explained by the fact that the closer the fiber is oriented to 90°, the greater the modal loss factor since the transversal specific damping capacity is greater than that along the longitudinal direction (see Table 2).

Figure 5: Variation of the modal loss factor with ply angle (l/d = 10,  = 0.5).

3.2. Stability Lobes

The stability analytical model is applied to a two-fluted composite end mill cutter for a half immersion down milling (=0° and  = 90°), Kt = 732 MPa and Kr = 0.076 as reported by Comak et al [24]. The composite cutter bar has the stacking sequence , with a mean radius 0.0352 m and a thickness 0.002 m. The length L was determined according to the length-to-diameter ratio.

First, the convergence study of the stability lobes is given in Figure 6. It is seen that the convergence of the stability lobes is good with respect to N, and N = 3 is found to be adequate for convergence. Therefore, in this section, all results have been obtained by three mode shape functions.

Figure 6: The effect of the numbers of mode shape function on the stability lobes (l/d = 10, σ = 0.5, θ = 00).

Figure 7 shows the effect of rotation of the cutter bar on the cutting stability with the length-to-diameter ratio 10, the ply angle 0°, and the taper ratio 0.5, in which the stability lobes without and with the consideration of gyroscopic effect are marked in black solid lines and red dash lines, respectively. Apparently, as the rotating speed increases, the backward frequency drops and forward frequency increases under gyroscopic/rotating effect. The backward frequency is less than the forward one, and thus the position of the stability diagram associated with the backward mode is lower, thereby the backward whirling determines the stability of the cutting system.

Figure 7: Stability lobes with and without the consideration of gyroscopic/rotating effect.

It can also be observed from Figure 7 that the limit critical cutting depth of the cutter bar without gyroscopic effect is a constant, and the envelopes are horizontal lines. By taking the gyroscopic effect into account, the critical cutting depth is no longer a constant but relies on the rotating speed, whose envelopes decrease obliquely. Accordingly, instable region of the cutter bar considering the gyroscopic effect is larger than those of the cutter bar without gyroscopic effect, suggesting that gyroscopic effect plays the role of negative damping in the cutting process. Therefore, the gyroscopic effect imposes nonnegligible effects on the system’s dynamic characteristics and stability limit during high-speed milling process.

Figure 8 shows the stability lobes of the cutter bars made up of carbon fiber/epoxy composite materials and steel. The taper ratio 0.5, the length-to-diameter ratio 10, and the ply angle 0° are used, respectively. Table 4 gives the mechanical properties of steel.

Figure 8: The limit cutting depth versus rotating speed for carbon/epoxy and steel cutter bars.

As shown in Figure 8, the maximum value of the limit cutting depth for the carbon fiber/epoxy cutter bar is 0.412 mm, and the maximum value of the limit cutting depth for the steel cutter bar is 0.188 mm. This is due to the fact that, among three different types of cutter bars, the carbon/epoxy cutter bar has the greater bending stiffness and damping than the steel cutter bar. Also, the cutting speeds of the carbon fiber/epoxy composite cutter bar are larger than those of the steel cutter bar due to higher specific stiffness and in consequence higher fundamental natural frequency.

Figure 9 shows the effect of the ply angle on the stability lobes of the carbon/epoxy resin cutter bar with a taper ratio of 0.5 and a long-to-diameter ratio of 10.

It can be observed from Figure 9 that the stable region of the composite cutter bar increases with the decreasing ply angle. As the ply angle drops, the bending stiffness of the composite cutter bar increases, which can lead to the increase of natural frequency (see Figure 3), thereby raising the rotating speed of machine tool spindle.

Figure 10 shows the effect of stacking sequence of the composite cutter bar on the stability lobes (with a length-to-diameter ratio of 10), from which it can be easily observed that the change of stacking sequence can significantly affect the stability lobes and thereby affect the chatter stability of the cutting system. The greater the number of fibers oriented close to the longitudinal direction of the composite cutter bar (e.g., the stacking sequences unidirectional), the more they contribute to the cutter bar bending stiffness and, consequently, the larger the stable regions will be.

Figure 9: Effect of ply angle on stability lobes of the carbon/epoxy cutter bar.

Figure 10: Effect of stacking sequence of the composite boring bar on the stability lobes.

Figure 11 shows stability diagrams for the tapered composite cutter bar (σ = 0.3) and uniform composite cutter bar (σ = 1). Two cutter bars have same stacking sequence and same volume. The geometrical characteristics of the tapered carbon/epoxy cutter bar are RT = 0.04028 m, RR = 0.03012 m, h = 0.002032 m, and L = 0.704 m. The geometrical characteristics of the uniform carbon/epoxy cutter bar are RT = RR = 0.0352 m, h = 0.002032 m, and L = 0.704 m. As it can be seen in Figure 11 that stable regions of the tapered composite cutter bar are larger than those of the uniform composite cutter bar, this may be due to the higher bending stiffness produced by cutter bar tapering.

Figure 11: The stability lobe diagrams of the tapered composite cutter bar and tapered composite cutter bar and uniform composite cutter bar.

4. Conclusions

An analytical model for predicting free vibration and chatter stability of a rotating tapered composite cutter bar, subject to dynamic regenerative cutting forces, is presented. Based on damping dissipation constitutive relations for viscoelastic composite, damping is introduced to the model. Moreover, the 2-DOF stability analysis model of Comak et al. [24] and Eynian and Altintas [25] for steel cutter bar based on the rotating frame approach is extended here to the case of the cutting process with rotating tapered composite cutter bar, which is used to predict stability lobes analytically. The continuum-based model of the cutter bar is developed by employing the Hamilton principle and Euler–Bernoulli beam theory, and the Galerkin method is used to discretize the partial differential motion equations with bending-bending coupling. The influences of some factors including the rotating gyroscopic effect, materials, ply angle, stacking sequence, and taper on the chatter stability are also examined. Results show that the gyroscopic effect of the composites acts as negative damping in the cutting process. By neglecting the rotating effect, the chatter stability is overestimated. Compared with the maximum limit cutting depth of the steel boring bar, the values of the carbon/epoxy cutter bar is enhanced by 263.8%. As the ply angle drops, the chatter stability increases. The chatter stability of the cutting process can also be improved by changing the taper ratio or stacking sequence of the composites cutter bar.

Appendix

The detailed expressions for and

The formulas for calculating [22]:

in which Cij are known functions of the layer material parameters and represents ply angle.

in which and represent the elastic modulus, represents the shear modulus, and represents the Poisson ratio.

The formulas for calculating [22]:

where can be obtained as follows:

where, represent damping parameters, in which

Nomenclature

: | Rotating speed of the tapered composite cutter bar
: | The cross-section radii on the fixed end of the cutter bar
: | The cross-section radii on the free end of the cutter bar
: | The taper of the cutter bar
: | Kinetic energy
: | Potential energy
: | Virtual work
, : | The displacement of the neutral axis along y-axis and z-axis directions of the points on the reference axis of the cutter bar
, : | The rotations about y-axis and z-axis
: | Reduced stiffness
: | Off-axis stiffness linked to material damping
: | The spatial varying flexural stiffness of the composite cutter bar
, : | The external and internal radius of layer k in section x
, : | The cutting forces produced on the free end of the cutter bar along y-axis and z-axis directions
: | The mode function
: | Internal damping ratio
: | Mode deflection
: | The stiffness matrices of a tapered nonrotating composite cutter bar
: | The damping matrices of a tapered nonrotating composite cutter bar
: | Axial cutting depth
: | The tangential cutting constant
, : | Directional matrix whose coefficients depend on the instantaneous angular immersion of the cutter
: | The pitch angle of the cutter
: | The cutter pitch angle
: | Average matrix
, : | Start and exit immersion angles of the cutter to and from the cut
: | The radial cutting constant
: | The chatter frequency
: | The critical axial depth.

Data Availability

(1) The [1446.2 172.7 7.2 3.76 3.76 0.3 0.45 4.22 7.05 7.05] data used to support the findings of this study have been deposited in the website (https://doi.org/10.1016/j.compscitech.2007.06.019). (2) The [0 90 732 0.076] data used to support the findings of this study have been deposited in the [DOI: 10.1115/1.4032585]. (3) The [7850 210 0.3 0.01] data used to support the findings of this study were measured in the authors experiments.

Conflicts of Interest

The authors have declared that no conflicts of interest exist.

Acknowledgments

The authors want to express their heartfelt thanks to National Natural Science Foundation of China (Grand no. 11672166) for financial supports of this work.

References

A. R. Abu Talib, A. Ali, M. A. Badie, N. Azida Che Lah, and A. F. Golestaneh, “Developing a hybrid, carbon/glass fiber-reinforced, epoxy composite automotive drive shaft,” Materials & Design, vol. 31, no. 1, pp. 514–521, 2010.

K. G. Bang and D. G. Lee, “Design of carbon fiber composite shafts for high speed air spindles,” Composite Structures, vol. 55, no. 2, pp. 247–259, 2002.

H. S. Kim, J. W. Sohn, and S.-B. Choi, “Vibration control of a cylindrical shell structure using macro fiber composite actuators,” Mechanics Based Design of Structures and Machines, vol. 39, no. 4, pp. 491–506, 2011.

D. G. Lee, K. G. Bang, and J. K. Choi, “Development of composite pressing rollers for thin film processing,” Journal of Composite Materials, vol. 34, no. 9, pp. 770–790, 2000.

S. K. S. Nadeem, G. Giridhara, and H. K. Rangavittal, “A Review on the design and analysis of composite drive shaft,” Materials Today: Proceedings, vol. 5, no. 1, pp. 2738–2741, 2018.

M. Yang, Y. Hu, J. Zhang, G. Ding, and C. Song, “Study on forced torsional vibration of CFRP drive-line system with internal damping,” Applied Composite Materials, vol. 25, no. 6, pp. 1307–1322, 2018.

P. S. Filho, J. H. S. A. Jr, and S. C. Amico, “Carbon/epoxy filament wound composite drive shafts under torsion and compression,” Journal of Composite Materials, vol. 52, no. 8, pp. 1103–1111, 2018.

J. R. Banerjee and H. Su, “Dynamic stiffness formulation and free vibration analysis of a spinning composite beam,” Computers & Structures, vol. 84, no. 19-20, pp. 1208–1214, 2006.

M.-Y. Chang, J.-K. Chen, and C.-Y. Chang, “A simple spinning laminated composite shaft model,” International Journal of Solids and Structures, vol. 41, no. 3-4, pp. 637–662, 2004.

H. B. H. Gubran and K. Gupta, “The effect of stacking sequence and coupling mechanisms on the natural frequencies of composite shafts,” Journal of Sound and Vibration, vol. 282, no. 1-2, pp. 231–248, 2005.

J. M. Ma and Y. S. Ren, “Vibration and stability of thin-walled composite shafts of variable cross-section with different end constaints,” Engineering Mechanics, vol. 33, no. 5, pp. 18–33, 2016, (in Chinese).

Y. S. Ren, Q. Y. Dai, and X. Q. Zhang, “Modeling and dynamic analysis of rotating composite shaft,” Journal of Vibroengineering, vol. 15, no. 4, pp. 1790–1806, 2013.

R. Sino, T. N. Baranger, E. Chatelet, and G. Jacquet, “Dynamic analysis of a rotating composite shaft,” Composites Science and Technology, vol. 68, no. 2, pp. 337–345, 2008.

O. Song, N.-H. Jeong, and L. Librescu, “Implication of conservative and gyroscopic forces on vibration and stability of an elastically tailored rotating shaft modeled as a composite thin-walled beam,” The Journal of the Acoustical Society of America, vol. 109, no. 3, pp. 972–981, 2001.

M. Almuslmani and R. Ganesan, “Vibration of tapered composite driveshaft based on the hierarchical finite element analysis,” Composite Structures, vol. 209, pp. 905–927, 2019.

S. X. Yang, W. K. Shi, Z. Y. Chen, Q. Chen, C. H. Yang, and L. L. Hang, “Composite mechanics and energy method based stiffness prediction model for composite leaf springs,” Mechanics Based Design of Structures and Machines, vol. 47, no. 11, pp. 1–12, 2019.

D. G. Lee and N. P. Suh, “Manufacturing and testing of chatter free boring bars,” CIRP Annals, vol. 37, no. 1, pp. 365–368, 1988.

D. G. Lee, H. Yun Hwang, and J. Kook Kim, “Design and manufacture of a carbon fiber epoxy rotating boring bar,” Composite Structures, vol. 60, no. 1, pp. 115–124, 2003.

S. Nagano, T. Koizumi, T. Fujii, N. Tsujiuchi, H. Ueda, and K. steel, “Development of a composite boring bar,” Composite Structures, vol. 38, no. 1-4, pp. 531–539, 1997.

S. Ghorbani, V. A. Rogov, A. Carluccio, and P. S. Belov, “The effect of composite boring bars on vibration in machining process,” The International Journal of Advanced Manufacturing Technology, vol. 105, no. 1-4, pp. 1157–1174, 2019.

W. Kim, A. Argento, and R. A. Scott, “Forced vibration and dynamic stability of a rotating tapered composite Timoshenko shaft: bending motions in end-milling operations,” Journal of Sound and Vibration, vol. 246, no. 4, pp. 583–600, 2001.

Y. S. Ren, Y. Y. Shi, and Y. H. Zhang, “Effects of internal damping on dynamic stability of a rotating composite shaft,” Journal of Vibration and Shock, vol. 36, no. 24, pp. 181–200, 2017, (in Chinese).

Y. Altintas and E. Budak, “Analytical prediction of stability lobes in milling,” Annals of the CIRP, vol. 44, no. 1, pp. 357–362, 1995.

A. Comak, O. Ozsahin, and Y. Altintas, “Stability of milling operations with asymmetric cutter dynamics in rotating coordinates,” Journal of Manufacturing Science and Engineering, vol. 138, no. 8, pp. 081004-1–081004-7, 2016.

M. Eynian and Y. Altintas, “Analytical chatter stability of milling with rotating cutter dynamics at process damping speeds,” Journal of Manufacturing Science and Engineering, vol. 132, pp. 021012-1–021012-14, 2010.

E. S. Zorzi and H. D. Nelson, “Finite element simulation of rotor-bearing systems with internal damping,” Journal of Engineering for Power, vol. 99, no. 1, pp. 71–76, 1977.

D. A. Saravanos, D. Varelis, T. S. Plagianakos, and N. Chrysochoidis, “A shear beam finite element for the damping analysis of tubular laminated composite beams,” Journal of Sound and Vibration, vol. 291, no. 3-5, pp. 802–823, 2006.