Analysis of Cutting Stability of a Composite Variable-Section Boring Bar with a Large Length-to-Diameter Ratio Considering Internal Damping

Chattering in composite deep-hole boring can directly affect surface processing quality and efficiency and has always been a research hotspot in machining mechanics. In this study, based on Euler–Bernoulli beam theory, the fine control equations for the cutting stability of composite variable-section boring bars were established using the Hamilton principle, in which the sectional change and internal damping of the material were considered. Next, using the Galerkin method and semi-discrete method, the effects of the taper ratio, damping ratio, length-to-diameter ratio, and ply angle on the free vibration characteristics and cutting stability were analyzed in detail. The results show that at a low damping ratio, both the first-order inherent frequency and boring stability can be enhanced with the increase in the taper ratio; at a large damping ratio, increasing the taper ratio can reduce the first-order inherent frequency and boring stability. Finally, the effects of the sectional change on the inherent frequency, displacement response, and convergence were analyzed. A numerical simulation was performed for the model reliability validation. The present research results can provide a theoretical basis and technical guidance for analyzing the cutting stability and fine control of composite variable-section boring bars with large length-to-diameter ratios.


Introduction
Some deep holes, such as cylinder holes, axial oil holes of the shaft, rocket projectiles, and barrels of various artilleries, require high machining precision and surface quality, whereas some materials to be processed show poor cutting machinability. These have become a great challenge in production. Meanwhile, the boring bar of deep-hole processing is restricted by the hole diameter, showing a great length-to-diameter ratio, poor first-order bending rigidity, and low strength. Under the vortex motion of the boring bar, some problems including vibration, ripples, and tapers are easily generated, which can affect the linearity and surface roughness of deep holes, thereby inducing a series of processing problems such as noise [1]. A good-quality boring bar with a large length-to-diameter ratio should have high static rigidity, high dynamic rigidity, and high first-order natural bending frequency so as to avoid the chattering and vortex motion generated in high-speed rotation [2,3]. The suitable selection of the cutting parameters and size of the boring bar to enhance stability in cutting are key to enhancing the processing efficiency and quality in deep-hole boring as well as deep-hole processing capabilities. Chatter prediction and control in deep-hole boring has become a hot topic in machining dynamics.

Constitutive Equation of the Boring Bar after the Consideration of Material Dampi
The composite boring bar in this study is anisotropic and the constitutive can be expressed as denotes the damping stiffness matrix of the mate

Constitutive Equation of the Boring Bar after the Consideration of Material Damping
The composite boring bar in this study is anisotropic and the constitutive equation can be expressed as where Given the viscoelastic properties of the composite materials, the stress can be divided into elastic stress {σ e } and dissipative stress {σ d }: The dissipation characteristics of the various layers can be described by the following damping ratio matrix Therefore, the damping stiffness matrix [Q] ψ can be described by the damping ratio matrix [ψ] as Assume that α denotes the intersection angle between any fiber ply in the laminated plate and the axial x-axis direction, as shown in Figure 1. The relationship between stress and strain in the global coordinate system can be expressed as ε (8) where -Q and -Q ψ denote the stiffness matrix and the damping stiffness matrix in the global coordinate system that have been converted by the matrices in the local coordinate system, respectively, which can be expressed as The matrices for the coordinate transformation can be written as Considering the plane stress state, the elastic stress and dissipative stress can be expressed as where K ij and K ij ψ are shown in Appendix A and k s denotes the shear correction factor.

Elastic Strain Energy of the Boring Bar
Under the plane stress state, the elastic strain energy of the boring bar can be written as (σ xx e ε xx + τ xz e γ xz + τ xy e γ xy )dSdx (12) Ma et al. confirmed that when the length-to-diameter ratio (L/2R) exceeded 5, the Euler-Bernoulli beam theory showed high precision for anisotropic steel bars [19]; for the anisotropic composite bar, the value of L/2R should be greater than 8 to achieve the error of below 5%. Since the length-to-diameter ratio of the boring bar in this study exceeded 8, the Euler-Bernoulli beam theory can satisfy the actual requirements in engineering. Based on the Euler-Bernoulli beam theory, the shear stress induced by the bending and torsion on the cross-section can be ignored, i.e., τ e xz = τ e xy = 0. Simultaneously, considering small deformations, the axial strain induced by bending can be expressed as where v and w denote the lateral displacements of the boring bar along the y-axis and z-axis and θ y and θ z are the rotation angles of the boring bar on the xand x-z planes. Substituting Equations (13) and (14) into the strain energy and stress equations, the strain energy of the boring bar can be expressed as where In which R i−1 and R i are the inner diameter and outer diameter of the i-th layer. R i−1 and R i are also the functions of the section position x.

Kinetic Energy of the Boring Bar
Taking both the deformation and rotation of the boring bar into account, the kinetic energy of the boring bar in the global coordinate system can be written as (17) where I m denotes the mass of the boring bar per unit length and I d and I P denote the inertia moment and the polar inertial moment of the cross-section, which can be expressed as where n denotes the total number of plies and ρ i denotes the mass density of the i-th layer, as shown in Figure 2.
layer, as shown in Figure 2. Substituting the relationship between the rotation angle and the deflec Equation (17), it can be seen that

Dissipated Virtual Work of the Boring Bar
The dissipated virtual work of the boring bar can be expressed as: Substituting the expressions of σ d xx and ε xx into Equation (20), the d virtual work of the boring bar can be expressed as

Vibration Control Equation of the Boring Bar
Based on the Hamilton principle, the relationship between strain energy kinetic energy density, and the dissipated virtual work can be obtained as Substituting the relationship between the rotation angle and the deflection into Equation (17), it can be seen that

Dissipated Virtual Work of the Boring Bar
The dissipated virtual work of the boring bar can be expressed as: Substituting the expressions of σ d xx and ε xx into Equation (20), the dissipated virtual work of the boring bar can be expressed as where

Vibration Control Equation of the Boring Bar
Based on the Hamilton principle, the relationship between strain energy density, kinetic energy density, and the dissipated virtual work can be obtained as Substituting the above relation into the Hamilton principle, the vibration differential equation of the system can be expressed as In above equations, the single-underlined terms indicate the effect of the sectional change on the vibration, and the double-underlined terms indicate the effects of both damping and the sectional change.

Discrete Solution to Vibration Equation
Using the Galerkin approximate solution method, which is commonly used for continuous system vibration and stability analysis, the partial differential equation of free vibration can be transformed into the ordinary differential equation (ODE) with respect to time.
For the cantilever bar, assuming the bending deformation can be expressed as The vibration mode functions of the bending deformation of standard non-rotary, non-coupling, and uniform cantilever beams are as follows: Substituting Equation (25) into the vibration differential equation, the weighted integral was performed on the vibration mode functions using the Galerkin method and then the 2N ordinary differential equations of the boring bar can be obtained as From Equation (28), it can be seen that the non-diagonal terms in the damping matrix G reflect the influence of the gyroscopic effect, and the diagonal items reflect the influence of the material resistance. The vortex motion appeared when the boring bar was influenced by the gyroscopic effect.

Boring Stability Analysis
This study adopted the model of boring force established by Kapoor et al., as shown in Figure 3 [20]. The boring force is in direct proportion to the cutting cross-sectional area on the r-x plane. Assuming that the endpoint of the bar C shows no displacement on the y-z plane, the boring force F can be expressed as where K s denotes the cutting force on the unit cutting area; f r denotes the feed per rotation during the boring process; and d r denotes the cutting radial depth. During the boring process, due to the deformation of the boring bar, the endpoint C will deviate along the radial depth direction. Assuming that u r denotes the offset, the boring force after the consideration of the deviation can be expressed as  During single-tool boring machining, the boring force refers to the stress imposed on the endpoint at x = L of the boring bar. Combined with the research results in Refs. [22] and [23], and assuming that y f and z f denote the cutting forces along the y-axis and zaxis direction per unit length, respectively, Equation (27) can be rewritten as in which δ is the Dirac-delta function and can be expressed as Substituting Equations (36)-(38) into Equation (35) and combined with Equation (25), the vibration equation of the boring bar imposed by cutting forces can be obtained as where M, G, and K are identical to those in Equation (28). The other parameters are as follows: ; ; According to Equation (40), the off-diagonal elements in the f K matrix are opposite in sign, which can induce chattering instability during the cutting process. The elements In combination with Equation (29), the cutting force in Figure 3a can be obtained as where β 0 denotes the intersection angle between the normal direction of the cutting face and the y-z plane and β 0 = 60 • [20]. With reference to the research results by Subramani et al. [21], the axial cutting force can be assumed as From Figure 3, through the transformation of coordinates, the boring forces along the y-axis and z-axis can be obtained as In which θ t = Ωt denotes the angular displacement during the rotation of the boring bar; K s is the coefficient of cutting force during high-speed boring on steel; according to the test data in [20], K s = 1330 × 10 6 N/m 2 ; and f r denotes the boring feed, with the unit of mm/rotation.
During single-tool boring machining, the boring force refers to the stress imposed on the endpoint at x = L of the boring bar. Combined with the research results in refs. [22] and [23], and assuming that f y and f z denote the cutting forces along the y-axis and z-axis direction per unit length, respectively, Equation (27) can be rewritten as where In which δ is the Dirac-delta function and can be expressed as where M, G, and K are identical to those in Equation (28). The other parameters are as follows: According to Equation (40), the off-diagonal elements in the K f matrix are opposite in sign, which can induce chattering instability during the cutting process. The elements in the K f matrix change in a periodic pattern, which can also induce unstable boring parameters.
During the boring process, a constant force will not induce the chatter of the boring bar. Accordingly, the constant force on the right side of Equation (39) can be ignored. Only considering the effect of the time-delay force on the boring bar, the boring vibration equation can be expressed as Equation (41) is a standard 2-dof delay differential equation during the cutting vibration. The time-delay terms related to the bending and deformation of the boring bar can induce chattering.

Natural Vibration Analysis
The first-order inherent frequency of the boring bar imposes great significance on the cutting stability. In this study, in order to reveal the effects of the various factors on the boring bar's inherent vibration characteristics, the inherent vibration characteristics of the hollow circular-section boring bar made up of epoxy-resin-based graphite fibers were analyzed. The boring bar is shown in Figure 1, and Table 1 lists the performance parameters of the material. From the table, N = 10, t c = 0.01321 m, d = 2R = r 1 + r 2 = 0.127 m, ψ 1 = 0.0045, ψ 2 = 0.0422, and ψ 12 = 0.0705, where N denotes the number of plies, whereas ψ 1 , ψ 2 , and ψ 12 are the internal damping parameters. For convenience for the subsequent description, we used ψ = 1 to represent this group of values and used ψ = 0.5 to represent the parameters of ψ 1 = 0.00225, ψ 2 = 0.0211, and ψ 12 = 0.03525. The rest were treated in the same manner. If there was no specified explanation, the boring bars all adopted the [0 • ] 10 -ply mode. The taper ratio (TR) is to reveal the linear sectional change of the boring bar, which can be defined as

Graphite Epoxy Steel
1610 7700 Figure 4 shows the variation in the first-order vortex frequency of the boring bars with different length-to-diameter ratios when TR = 1 and ψ = 0, in which FW and BW denote the forward and backward procession, respectively. In the figure, it can be seen that as the length-to-diameter ratio increased, the first-order vortex frequency decreased gradually, and simultaneously, the critical rotation speed decreased. As the length-to-diameter ratio increased, the boring bar became thinner and longer, thereby reducing the first-order natural bending frequency and the critical rotation speed. Accordingly, the chattering was more likely to occur.  Figure 5 displays the variation in the first-order inherent frequency with the internal damping coefficient under a fixed length-to-diameter ratio and a changing taper ratio. As the internal damping coefficient increased, the first-order inherent frequency of the boring bar decreased gradually, and the larger the taper ratio, the more significant the effect influenced by the inner damping of the boring bar. This can be attributed to the existence of the double-underlined terms as shown in Equation (24) in the vibration equation. A larger taper ratio is indicative of a greater effect on the change in damping. Figure 6 shows the changes in the inherent frequency of uniform-section boring bars with different length-to-diameter ratios with the damping coefficient. For the boring bars with different length-to-diameter ratios, the inherent frequency decreased gradually with the increase in the damping coefficient, and the smaller the length-to-diameter ratio, the more significant the effect of the internal damping coefficient on the first-order inherent frequency of the boring bar. This suggests that the effect of A ψ 11 in Equation (24) was weakened for thinner and longer bars. Accordingly, the effect of the inherent frequency on damping weakened and the curve flattened. For the boring bars with large length-todiameter ratios, changing the inherent frequency with damping showed a limited effect.
damping coefficient under a fixed length-to-diameter ratio and a changing taper ratio the internal damping coefficient increased, the first-order inherent frequency of the bo bar decreased gradually, and the larger the taper ratio, the more significant the e influenced by the inner damping of the boring bar. This can be attributed to the exist of the double-underlined terms as shown in Equation (24) in the vibration equatio larger taper ratio is indicative of a greater effect on the change in damping.  Figure 6 shows the changes in the inherent frequency of uniform-section boring with different length-to-diameter ratios with the damping coefficient. For the boring with different length-to-diameter ratios, the inherent frequency decreased gradually the increase in the damping coefficient, and the smaller the length-to-diameter ratio more significant the effect of the internal damping coefficient on the first-order inhe frequency of the boring bar. This suggests that the effect of ψ 11 A in Equation (24) weakened for thinner and longer bars. Accordingly, the effect of the inherent frequ on damping weakened and the curve flattened. For the boring bars with large lengt diameter ratios, changing the inherent frequency with damping showed a limited ef  Figure 7 shows the variations in the first-order inherent frequencies of the boring with the taper ratio under a fixed length-to-diameter ratio and different dam coefficients. In the figure, it can be seen that when the damping ratio equaled 0, 0.5 1, the first-order inherent frequency of the boring bar increased gradually wit increasing taper ratio. When the damping ratio equaled 1.5 and 2, the inherent frequ first increased and then decreased with the increase in the taper ratio. This can al explained by the existence of the double-underlined terms in Equation (24). The inh frequency of the boring bar with low damping can be enhanced by increasing the ratio. For the materials with large damping, the taper ratio imposed a complex effec effects of both the taper ratio and damping should be considered.   Figure 7 shows the variations in the first-order inherent frequencies of the boring bars with the taper ratio under a fixed length-to-diameter ratio and different damping coefficients. In the figure, it can be seen that when the damping ratio equaled 0, 0.5, and 1, the first-order inherent frequency of the boring bar increased gradually with the increasing taper ratio. When the damping ratio equaled 1.5 and 2, the inherent frequency first increased and then decreased with the increase in the taper ratio. This can also be explained by the existence of the double-underlined terms in Equation (24). The inherent frequency of the boring bar with low damping can be enhanced by increasing the taper ratio. For the materials with large damping, the taper ratio imposed a complex effect. The effects of both the taper ratio and damping should be considered.
first increased and then decreased with the increase in the taper ra explained by the existence of the double-underlined terms in Equatio frequency of the boring bar with low damping can be enhanced by ratio. For the materials with large damping, the taper ratio imposed a effects of both the taper ratio and damping should be considered.  In seen that the greater the number of 0° ply layers, the greater th frequency and the larger the critical rotation speed of the boring critical rotation speed of the boring bar decreased with the decreasi layers. The first-order vortex frequency of the material with a ply minimum. Through comparison, the first-order vortex frequency o mode of [90°/0°]5 was slightly higher than that of the bar with a ply m suggests that the closer the 0° ply to the outer surface of the boring critical rotation speed. This is because the 0° ply showed higher rig with other angles. Meanwhile, the cross-section where the 0° ply was   5 . This suggests that the closer the 0 • ply to the outer surface of the boring bar, the higher the critical rotation speed. This is because the 0 • ply showed higher rigidity than the plies with other angles. Meanwhile, the cross-section where the 0 • ply was located was greater in radius, perimeter, and volume, and thus can bear greater the bending moment induced by elastic stress and dissipative stress. Therefore, compared with the arrangement of 0 • ply on the inner layer, arranging the 0 • ply on the outer layer can enhance the stability of the boring bar.
Materials 2021, 14, x FOR PEER REVIEW in radius, perimeter, and volume, and thus can bear greater the bending moment by elastic stress and dissipative stress. Therefore, compared with the arrangem ply on the inner layer, arranging the 0° ply on the outer layer can enhance the sta the boring bar. Compared with existing models, the established stability analysis mode variable-stability boring bar considered the effect of the sectional change. In analyze the sectional change in the inherent vibration characteristics, the cal Compared with existing models, the established stability analysis model of the variable-stability boring bar considered the effect of the sectional change. In order to analyze the sectional change in the inherent vibration characteristics, the calculation results before and after the consideration of the sectional change were compared for further analysis. Since both A ψ 11 and A ψ 11 included the effects of damping on the boring bar's inherent vibration characteristics, the internal damping coefficient was set as 0 only to analyze the effect of the sectional change on the boring bar's vibration characteristics. Then, the free vibration characteristics were analyzed by setting the terms of I d , I P , A 11 , and A 11 in Equations (24a) and (24b) as 0. Figure 9 shows the effect of the sectional change on the first-order inherent frequency of the boring bar when TR = 1, 2, and 3, respectively. The larger the taper ratio of the boring bar, the greater the effect. Table 2 lists the calculation results of the first-order inherent frequencies of the boring bar with different TR values using different models. The ANSYS model of the boring bar is shown in Figure 10. Solid 185 element is used to simulate the composite laminated structure. From the table, the calculation results using the established model were closer to the data using the ANSYS model, showing more favorable reliability.

Cutting Stability Analysis Validation of Model Reliability
The Galerkin method is a numerical discrete algorithm, whose numerical calculation results are related to precision and the number of intercepted modes. In the case of insufficient intercepted modes, the calculation results cannot be convergent with great errors. However, a larger number of intercepted modes is also not good. In the case of a larger number of modes, the calculation time is long and the calculation results are divergent. Therefore, it is necessary to validate the convergence of the selected number of modes used in this study. Figure 11 shows the calculated stability lobes of different boring bars with different mode numbers. In the figure, it can be seen that the stability lobe diagrams were quite similar when the model numbers were set as 1, 2, and 3. Table 3 shows the limit feeds of the unconditional stable region at a rotation speed of 0~1500 rad/s using different mode numbers. The numerical results were quite close, suggesting the favorable convergence of the solution to the stability analysis of the boring system using the time-domain semi-discrete method.

Validation of Model Reliability
The Galerkin method is a numerical discrete algorithm, whose results are related to precision and the number of intercepted m insufficient intercepted modes, the calculation results cannot be c errors. However, a larger number of intercepted modes is also not larger number of modes, the calculation time is long and the c divergent. Therefore, it is necessary to validate the convergence of t modes used in this study. Figure 11 shows the calculated stability lo bars with different mode numbers. In the figure, it can be seen diagrams were quite similar when the model numbers were set a shows the limit feeds of the unconditional stable region at a rotation using different mode numbers. The numerical results were quite favorable convergence of the solution to the stability analysis of the the time-domain semi-discrete method.     Table 3. Effect of the number of vibration mode functions on the limit boring feed of the rotated composite material boring system (ψ = 1 and L/d = 10). Cutting stability analysis curves can be obtained via feature extraction. It is necessary to validate the reliability of cutting stability. This study used the Runge-Kutta method to perform the displacement response analysis on the control equation and validated the results by combining the displacement response results and stability analysis curves. Figure 12 shows the cutting stability lobes of the uniform section of the boring bars with different damping coefficients. A total of 3 points, denoted as Point A, Point B, and Point C, were selected for the displacement response analysis, and the results are shown in Figure 13.  Table 4 lists the rotation speeds and feeds at different coordinates. Th located in the stable and unstable regions of the lobe diagrams of the b different damping coefficients. In Figure 12 it can be seen that Point A w stable regions of the boring bars with damping coefficients of 0.5 and displacement response should be convergent. These results were con results shown in Figure 13. Point B was located in the unstable region o with a damping coefficient of 0.5 and the stable region of the boring bar coefficient of 1, indicating its displacement response should be divergent These results were also consistent with the divergent and converge responses at Point B in Figure 13. Point C was located in the unstab damping coefficients of 0.5 and 1, indicating its displacement resp divergent. These results were consistent with the divergent displacem Point C in Figure 13. The displacement response results fit well with drawn from the lobe diagrams, suggesting the reliability of the lobe diag  Table 4 lists the rotation speeds and feeds at different coordinates. These points were located in the stable and unstable regions of the lobe diagrams of the boring bars with different damping coefficients. In Figure 12 it can be seen that Point A was located in the stable regions of the boring bars with damping coefficients of 0.5 and 1, indicating its displacement response should be convergent. These results were consistent with the results shown in Figure 13. Point B was located in the unstable region of the boring bar with a damping coefficient of 0.5 and the stable region of the boring bar with a damping coefficient of 1, indicating its displacement response should be divergent and convergent. These results were also consistent with the divergent and convergent displacement responses at Point B in Figure 13. Point C was located in the unstable regions with damping coefficients of 0.5 and 1, indicating its displacement response should be divergent. These results were consistent with the divergent displacement response at Point C in Figure 13. The displacement response results fit well with the conclusions drawn from the lobe diagrams, suggesting the reliability of the lobe diagrams.  Figure 14 shows the boring stability lobe diagrams of the boring bars with different length-to-diameter ratios. In the figure it can be seen that the limit feed of the boring bar at the same rotation speed decreased with the increasing length-to-diameter ratio, suggesting that the increase in the length-to-diameter ratio increased the instability. When the length-to-diameter ratio exceeded 10, the limit cutting depth dropped rapidly, indicating great difficulties in deep-hole processing. A greater length-to-diameter ratio is indicative of thinner pores. At that moment, the boring bar with a larger length-todiameter ratio is featured by quite low rigidity with a very high risk of chattering. The limit cutting depth is quite low with a fairly low processing efficiency. In addition, it can also be observed from the figure that the length-to-diameter ratio of the boring bar imposed a slight effect on the initiation speed of new chatter instability in the high-speed rotation region.   Table 4. Rotation speeds and cutting depths of the different points in Figure 11.

Point Rotation Speed (rad/s) Feed (mm/rev)
A 600 0.25 B 600 0.75 C 600 1.5 Figure 14 shows the boring stability lobe diagrams of the boring bars with different length-to-diameter ratios. In the figure it can be seen that the limit feed of the boring bar at the same rotation speed decreased with the increasing length-to-diameter ratio, suggesting that the increase in the length-to-diameter ratio increased the instability. When the lengthto-diameter ratio exceeded 10, the limit cutting depth dropped rapidly, indicating great difficulties in deep-hole processing. A greater length-to-diameter ratio is indicative of thinner pores. At that moment, the boring bar with a larger length-to-diameter ratio is featured by quite low rigidity with a very high risk of chattering. The limit cutting depth is quite low with a fairly low processing efficiency. In addition, it can also be observed from the figure that the length-to-diameter ratio of the boring bar imposed a slight effect on the initiation speed of new chatter instability in the high-speed rotation region.

Analysis of the Influencing Factors of the Boring Stability of the Boring Bar
suggesting that the increase in the length-to-diameter ratio increased the instabilit the length-to-diameter ratio exceeded 10, the limit cutting depth dropped indicating great difficulties in deep-hole processing. A greater length-to-diamete indicative of thinner pores. At that moment, the boring bar with a larger le diameter ratio is featured by quite low rigidity with a very high risk of chatter limit cutting depth is quite low with a fairly low processing efficiency. In additio also be observed from the figure that the length-to-diameter ratio of the bo imposed a slight effect on the initiation speed of new chatter instability in the hig rotation region.  Figure 15 shows the effect of the change in internal damping on the cutting stability of the boring bar when TR = 1 and L = 2 m (L/d = 15.7). Through comparison, when the rotation speed was low, increasing the internal damping of the boring bar can enhance the cutting stability of the boring bar, thereby leading to greater stable feed. However, at a high rotation speed, a new boring chatter unstable region appeared due to the existence of damping. The greater the damping, the smaller the initiation speed of the unstable region and the larger the range of the unstable region enclosed by the high-speed lobe curve and the horizontal axis. Generally, a higher inherent frequency suggests a larger limit cutting depth, whereas the inner damping shows an opposite effect. Under higher inner damping, the inherent vibration frequency can be reduced but still enhance the cutting stability at a low rotation speed, thereby reducing the cutting stability at a high rotation speed.  Figure 15 shows the effect of the change in internal damping on the cutting stability of the boring bar when TR = 1 and L = 2 m (L/d = 15.7). Through comparison, when the rotation speed was low, increasing the internal damping of the boring bar can enhance the cutting stability of the boring bar, thereby leading to greater stable feed. However, at a high rotation speed, a new boring chatter unstable region appeared due to the existence of damping. The greater the damping, the smaller the initiation speed of the unstable region and the larger the range of the unstable region enclosed by the high-speed lobe curve and the horizontal axis. Generally, a higher inherent frequency suggests a larger limit cutting depth, whereas the inner damping shows an opposite effect. Under higher inner damping, the inherent vibration frequency can be reduced but still enhance the cutting stability at a low rotation speed, thereby reducing the cutting stability at a high rotation speed.  Figure 16 shows the effect of the taper ratio of the boring bar on the cutting stability. From the figure it can be seen that ψ = 1 at the same rotation speed and increasing the taper ratio of the boring bar can obtain a greater stable feed, thereby enhancing the cutting efficiency. However, as the taper ratio increased, the initialization speed of the new highspeed chattering unstable region changed slightly. By combining Figures 14-16, it can be concluded that the initialization speed of the cutting unstable region at a high rotation speed was most significantly affected by the damping coefficient but slightly affected by the taper ratio. The length-to-diameter ratio imposed almost no effect.   Figure 16 shows the effect of the taper ratio of the boring bar on the cutting stability. From the figure it can be seen that ψ = 1 at the same rotation speed and increasing the taper ratio of the boring bar can obtain a greater stable feed, thereby enhancing the cutting efficiency. However, as the taper ratio increased, the initialization speed of the new highspeed chattering unstable region changed slightly. By combining Figures 14-16, it can be concluded that the initialization speed of the cutting unstable region at a high rotation speed was most significantly affected by the damping coefficient but slightly affected by the taper ratio. The length-to-diameter ratio imposed almost no effect.
taper ratio of the boring bar can obtain a greater stable feed, thereby enhancing the cutting efficiency. However, as the taper ratio increased, the initialization speed of the new highspeed chattering unstable region changed slightly. By combining Figures 14-16, it can be concluded that the initialization speed of the cutting unstable region at a high rotation speed was most significantly affected by the damping coefficient but slightly affected by the taper ratio. The length-to-diameter ratio imposed almost no effect.     Figure 17d shows the stability lobe diagrams of the boring bars with ply modes of [0°/90°]5 and [90°/0°]5, respectively. The limit stable cutting depth increased slightly when 0° plies were arranged from inside to outside. It should be noted that more 0° plies were more favorable for bending deformation; however, the composite boring bar with a ply angle of 0° was the weakest in torsion resistance.

Conclusions
Aiming at overcoming the limitations in the current variable-section boring bar   5 , respectively. The limit stable cutting depth increased slightly when 0 • plies were arranged from inside to outside. It should be noted that more 0 • plies were more favorable for bending deformation; however, the composite boring bar with a ply angle of 0 • was the weakest in torsion resistance.

Conclusions
Aiming at overcoming the limitations in the current variable-section boring bar models, this study introduced both sectional change and inner damping to establish the kinetic control equations of composite boring bars with large length-to-diameter ratios. Moreover, the effects of the various factors on the inherent vibration frequency and cutting stability were qualitatively and quantitatively analyzed in theory. This study innovatively explored the variation rules of the inherent frequency and cutting stability by taking both the sectional change and damping into account. The main conclusions are described below.
(1) The increase in the length-to-diameter ratio can reduce the first-order inherent frequency of free vibration and limit the boring feed while imposing no effect on the initiation speed of the new chattering unstable region. (2) The increase in the internal damping coefficient can reduce the first-order free vibration inherent frequency of the boring bar and simultaneously enhance the limit boring feed. However, at a high rotation speed, the existence of damping will generate a new boring chattering unstable region. Moreover, the unstable region will appear earlier under a greater damping condition. (3) At a damping ratio of 1 (i.e., ψ = 1), as the taper ratio of the boring bar increased, the first-order inherent frequency and the limit feed during stable boring were enhanced. The taper ratio imposed a slight effect on the initiation speed of the new chattering unstable region during high-speed boring. Meanwhile, both the taper ratio of the boring bar and the damping coefficient had interaction effects on the free vibration characteristics and cutting stability. For the composite variable-section boring bar, the effects of the internal damping and taper ratio should be considered. This is an innovative point in this study. (4) The ply angle can significantly affect the stable region during the boring process.