#### Abstract

The modeling as well as application of process damping in milling of thin-walled workpiece made of titanium alloy is investigated. Titanium alloy used commonly in aviation industry is one typical difficult-to-machine material. Chatter usually occurs in cutting of titanium alloy, which results in poor surface quality and damaged tool. Thus, chatter is one important restriction for the quality and efficiency of titanium alloy manufacture, especially for the thin-walled workpiece made of titanium alloy due to poor structural stiffness. Process damping results from interference between flank face and machined surface, which is critical but usually ignored in chatter analysis for difficult-to-machine material. The paper presents one nonlinear dynamic model considering process damping for milling of thin-walled workpiece made of titanium alloy and designs antivibration clearance angle to suppress chatter based on the model. The experimental and computational results indicate that the presented methods for chatter stability analysis are reasonable, and the antivibration clearance angle designed is effective in suppressing chatter and improving machining quality.

#### 1. Introduction

Titanium alloy is used widely in aerospace industry, which possesses excellent overall performances owing to high specific strength, low density, strong heat, and corrosion resistance; corresponding thin-walled structural parts made of titanium alloy also play one essential role in aerospace industry. However, titanium alloy is also one typical type of difficult-to-cutting material, which possesses low machinability owing to poor thermal conductivity, high chemical activity, and the severe hardening. Particularly, due to large unit cutting force and low elastic modulus, chatter usually occurs in machining process. Vibration marks left on the machined surface are usually required to be removed manually, which result in low machining efficiency. Worse, it may result in scraped workpiece and destroyed tool. Chatter problem is one of the main restrictions to the machining quality and efficiency of titanium alloy, which is much more troublesome in milling of thin-walled workpiece made of titanium alloy due to the poor structural stiffness.

The methods to suppress chatter can all be attributed to increase the damping of system. Damping of cutting system includes structural damping of machine tool and process damping.

Process damping (PD) results from interference between flank face of tool and undulations left on the workpiece. Tlusty and Ismail [1] showed that process damping has a significant effect on chatter stability decreasing with cutting speed. Structural damping can be identified using modal analysis easily. However, the identification as well as modeling of process damping has proved to be very difficult, and currently there is no practical method for its measurement or estimation. Altintas and Weck have addressed the process damping as the most challenging unsolved problems [2]. Therefore, in most of the studies, the effect of process damping is ignored, leading to significant errors in estimation of the stability limit [3, 4].

In one of the early works, the regenerative chip thickness effect on process stability is modeled by Das and Tobias, where the contact between the tool flank face and the undulations left on the surface was not considered. In this work [5], a velocity term in addition to the structural damping is introduced into the system equations increasing process stability. Later on, the focus of the research on dynamic cutting force coefficients shifted toward the identification and modeling of dynamic cutting force coefficients. The results of the CIRP efforts on dynamic cutting are summarized by Tlusty [6], where the difficulty of the measurements and the inconsistency of the test data from different labs are discussed. Sisson and Kegg [7], Wu [8], Lee et al. [9], and Shawky and Elbestawi [10] showed that the contact forces due to flank-wave interaction contribute to the dynamic of the cutting process by increasing the overall damping acting on the system. Altintas et al. [11] developed a dynamic force model which includes chip thickness, velocity, and acceleration terms. They identified dynamic cutting force coefficients from a series of dynamic cutting tests, where the cutting tool is oscillated by a fast tool servo at the desired frequency and amplitude.

Compared to turning, the reports on the process damping in milling, especially for the thin-walled workpiece, are much fewer, because milling is one discontinuity and time-invariant cutting process and the descriptions of process damping are much more complicated. For thin-walled workpiece, the modal parameters are also time-invariant, which increase the difficulties in stability analysis furthermore. Huang and Junz Wang [12, 13] extended analytical modeling of the milling process to include process damping effects. In their model, two cutting mechanisms (shearing and plowing) and two process damping effects (direction and magnitude) are included. The analytical nature of this model makes it possible to determine unknown process damping coefficients from measured vibration signal during milling. Budak and Tunc [14–17] considered the effect of process damping as an additional damper in the zero order MFS. They identified the additional damping coefficients by measuring stability limits experimentally and fitting the results to the stability model and studied the effects of tool geometry and cutting conditions on process damping. Ahmadi and Ismail [18, 19] simplified process damping to a piecewise linear viscous damper by assuming small amplitude of vibration and integrated the equivalent viscous model of process damping into the multifrequency solution and the semidiscrete method to establish the stability lobes in milling.

It can be known from above that higher stable cutting depth can be achieved under the effect of process damping, and high efficiency cutting can be achieved at low speeds. This can be important to increase machining productivity of titanium alloy, especially for the thin-walled workpiece with weak structural stiffness, because low cutting speeds have to be used for longer tool life due to low machinability. Besides, process damping has a lot to do with tool geometry; design of tool geometry is very important in suppressing chatter. However, very few cases have been reported for the application of process damping in designing of tool for milling of titanium alloy.

In this paper, one nonlinear milling dynamic model considering process damping for thin-walled workpiece made of titanium alloy is presented. In the model, indentation area arised from face-wave interference is calculated and stability limits are predicted by time domain simulations. Based on the model, antivibration clearance angle is designed to increase process damping and expand stable region further.

#### 2. Modeling of Process Damping

##### 2.1. Formation Mechanism of Process Damping

In finish machining, the stiffness of a thin-walled workpiece is far lower than that of the cutting tool. It is usually assumed that spindle tool is a rigid and the thin-walled workpiece can be considered to have a degree of freedom as shown in direction in Figure 1(a). The 1 DoF model shown in Figure 1(a) is widely used in the dynamic analysis for milling of thin-walled workpiece [20, 21]. is the angular position of the th cutting tool. is the spindle speed. is the vibration displacements of workpiece determined in the normal () direction, the wave on machined surface left behind in radial direction (i.e., chip thickness or radial direction), and direction represents the cutting speed direction. The interference between flank face and wave will occur when the chatter amplitude grows larger. The flank face of the tool indents into the wave, as shown in Figure 1. As a result, indentation forces (i.e., and ) arise in normal, , and frictional, , directions on the tool flank face as shown in Figure 1(b), which create a damping effects (i.e., process damping effects). The damping forces in and directions acting on each tooth are oriented and summed up as a resultant damping force in directions.

**(a)**

**(b)**

Indentation forces (i.e., and ) are proportional to indentation area; the calculation method of is presented in Section 2.2. In Figure 2(a), is the clearance angle of tool, and angle represents the slope of wave left in radial direction. shown in Figure 2(a) is the radial coordinate values of wave, which represent the radial vibration displacement. The instantaneous slope of wave represents the radial vibration speed (i.e., ). It is assumed that the tool does not indent the previous wave as shown in Figure 2(a). Therefore, the interference between flank face and workpiece occurs only as the tool slides down the wave, where . The occurrence of interference is determined by the relationship between and . It can be seen from Figure 2(a) that the interference will occur when . If satisfies the conditions of (1), the interference will occur:

**(a)**

**(b)**

In (1) is the cutting speed, represents tangential arc length of tool edge, and is time here (). It can be known from (1) that interference is easy to occur; as and decrease, indentation area and process damping will also increase. The law is also shown in Figure 2(b).

##### 2.2. Calculation of Indentation Area

The calculation of indentation area shown in Figure 1(a) is rather complicated, which is also the basis of calculation of indentation force (as seen in (8)). As shown in Figure 3, the cutting tool is divided into a number of axial elements; elemental axial width is . In axial level , point represents the current position of the tool edge, which can be also seen in Figure 4(a); correspondingly, is the current radial coordinate values of wave. As shown in Figure 4(a), represent the point on the surface wave generated by the tool edge in the previous time step; corresponding radial coordinate value of wave calculated in previous time steps is . For every , is the corresponding point on the flank face; is its radial coordinate values. The distance between and is . is zero. The indentation area in every axial level is bounded by the tool flank face and the workpiece surface wave. This area is calculated at each simulation time step numerically.

**(a)**

**(b)**

The calculation of is calculated according to (2); here, is the tangential arc length between adjacent two time points. is the simulation time step:

The intersection point of flank face and wave, , needs to be searched in the calculation process. As shown in Figure 4(b), when , the calculation of indentation area can stop. The index of front and back the intersection can be recorded as and , respectively; here , ; thus, indentation area can be calculated by summing up the discrete areas according to

##### 2.3. Effects of Antivibration Clearance Angle on Process Damping

It can be seen from (1) that interference is much easier to occur as cutting speed and clearance angle decreases, therefore the lower speed and smaller clearance angle can increase the process damping. However, in the tool design, severe extrusion is easy to occur and the surface quality will become very poor if the clearance angle is too small. Therefore, transition edge needs to be adopted in design. As shown in Figure 5, the tools modified have two clearance angles, and ; is designed to be smaller to suppress vibration, which can be named antivibration angle; the length of its corresponding transition edge is .

The calculation of indentation area follows the methods provided in Section 2.2; the main parameters in Figure 6 are the same as ones in Figure 4(a). Minor modification needs to be done because of the adoption of antivibration clearance angle. The calculation of follows (4); then the calculation of indentation area also follows (3):

#### 3. Time Domain Simulations

##### 3.1. Equations of Milling Dynamics

Dynamic cutting force model is of great importance in analysis of milling dynamics. At each axial level , the location of the th cutting tooth is determined by its angular position (as shown in Figure 1):

In (5), is the number of teeth, is the helix angle (as shown in Figure 3), and is the tool radius. is the angular position of tool bottom (axial level ), calculated from angular velocity and the time .

The dynamic cutting force is modeled using the linear-edge force model [22], where elemental tangential and radial forces, and , at flute, , are expressed in terms of cutting force coefficients, (), edge force coefficients, (), instantaneous chip thickness, , and elemental axial width, (as shown in Figure 3):

The differential cutting forces are then oriented in directions using the forces in and directions as follows:

For the th cutting flute, differential indentation forces acting on the axial level , in normal and in frictional directions are given as

In (8), is the indentation area at axial level , which is calculated by methods provided in Sections 2.2 and 2.3. is the coefficient of contact friction between the flank face and the workpiece surface; is indentation coefficient. The process damping forces in and directions are expressed in terms of the forces in direction and the friction coefficient (as shown in Figure 1(b)):

The differential process damping forces are then oriented in directions using the process damping forces in and directions as follows:

Thus, the equations of milling dynamic can be generated; it should be mentioned that the direction of cutting force exerted on the workpiece is opposite to the forces exerted on the tool: where is the numbers of axial elements (; is axial cutting depth), milling at small radial immersion is adopted in the paper, and, generally, only one tooth is cutting workpiece in one tooth passing period. , are the modal mass, modal stiffness of workpiece, respectively (as shown in Figure 1(a)), and is the coefficient of vibration mode, which is set to 1 at the tool bottom (). , , and can be all obtained from FEM analysis. is the corresponding modal damping ratio, which can be obtained from modal test; then the modal damping can be also generated. is a unit step function which determines whether the tooth is in or out of the cutting, which will be discussed in Appendix. Besides, the calculation of dynamic chip thickness in (6) and in (2) are rather complicated, which will be also discussed in detail in Appendix.

##### 3.2. Determining Chatter Stability Lobes from Time Domain Simulations

Equation (11) is solved using the classical 4th order Runge-Kutta method, where time increment is adopted as 10^{−5} seconds. Using the following steps, the stability limits are evaluated by finding the limiting axial cutting depth through a range of spindle speeds.

Modal parameters of workpiece, cutting tool geometry, feed rate, tool orientation, entry angle and exit angle, cutting force coefficients, and a starting axial cutting depth, , are specified.

A range of spindle speeds and a spindle speed step size are specified.

For a given spindle speed, a static time domain simulation is run (suppressing any tool and workpiece vibrations) and the maximum static uncut chip thickness, , is stored.

A second time domain simulation is run. For given spindle speed and cutting depth, (11) is solved, the vibration displacements and dynamic chip thickness can be obtained during the process. The largest dynamic chip thickness of all cutting points on the tool, , is stored for the last few revolutions of the simulation.

The nondimensional chatter parameter, , is evaluated as

If is greater than a predetermined limit (1.25 is used in this paper), the process is unstable; otherwise the process is stable. The criterion to judge is provided in [23].

If the process is stable, is set to the current value of , is doubled, and steps 3–6 are repeated until chatter occurs. Then is set to the value of when chatter occurred.

If the process is unstable, is set to the current value of , is halved, and steps 3–6 are repeated until the process is stable. Then is set to the stable axial cutting depth, .

Once the range to is found between which the limiting axial cutting depth lies, a bisection search is performed with , repeating steps 3–6 until the limiting axial cutting depth, , is found within a given tolerance.

#### 4. Comparison of Computational and Experimental Results

To verify the nonlinear model provided in Sections 1–3, the computations and experiments are all achieved. Computations are achieved in MATLAB according to time domain simulation provided in Section 3. Test scenes are shown in Figure 7, experimental equipments are as follows.

*Machine*. Mikron UCP DURO710.

*Tool 1*. Carbide end mill, 80 mm overall length, 40 mm overhang length, 20 mm edge length, 16 mm diameter, 4 flutes, 35° helix angle, and 9° clearance angle.

*Tool 2 (with Antivibration Clearance Angle)*. Clearance angle: is 4°; is 9°, length of transition edge: *μ*m, and the other parameters are the same as tool 1.

*Workpiece Rectangular Thin-Walled Workpiece Made of Typical Titanium Alloy (Ti6AL4V), with Physical Properties*. Young’s modulus of elasticity Gpa, density kg/m^{3}, and Poisson’s ratio . The dimension of workpiece is 70 × 40 × 5 mm, 5 mm wall thickness. The clamping height is 30 mm; thus the workpiece model shown in Figure 8 can be seen as one cantilever plate, the dimension of which is 40 × 40 × 5 (mm).

The modal stiffness and mass are different with the change of tool position; particularly the axial position of tool bottom has great effects on modal stiffness; therefore, the base line shown in Figure 8 is adopted. Axial cutting depths in experimental conditions are all chosen based on the line. The radial position of tool also has effects on the modal parameters. Therefore, only one thin coat of materials (thickness is equal to radial cutting depth) is removed for every workpiece. At last, the modal parameters of thin-walled workpiece are time-varying along feed direction in cutting, but the effects on stability prediction are much lighter than the first two. The modal stiffness and mass are calculated based on the point in Figure 8.

FEM is used to obtain modal stiffness and mass. The modal parameters (the first bending mode of the cantilever plate is adopted) in (11) are: natural frequency: Hz; modal damping ratio: ; modal stiffness: N/m.

The modal damping ratio is obtained from modal test. Figure 9 shows the measured frequency response function; modal parameters are identified by means of rational fraction method. The corresponding frequency response function fitted through these parameters is also shown in Figure 9. The identified from modal test is 2680 Hz, which is very close to the ones computed by FE model.

*Cutting Force Coefficients*. Tangential force coefficients Mpa, N/mm, and radial force coefficients MPa, N/mm, which are reported in [22].

*Process Damping Coefficients*. Indentation coefficient and friction coefficient , for titanium alloy Ti6AL4V, were reported at 30,000 N/mm^{3} and 0.3 in [15]. To show the effects of process damping, the stability limits predicted without process damping are also generated.

*Sensors*. Piezoelectric acceleration sensors used to measure accerleration signals of workpiece; dynamometer used to measure cutting force. Acceleration sensor is attached to the corner point of workpiece.

*Data Acquisition Card*. It is NIUSB9233.

*Main Parameters*. Radial cutting depth mm, feed per tooth mm/z, down milling, and mist cooling are considered.

The comparisons of stability limits generated by computation and experiments, respectively, are shown in Figure 10, from which it can be seen that the experimental results match with computational results well. The vibration signals measured in partial experimental conditions are shown in Figure 11. From Figure 11, it can be seen that vibration signals become divergent in unstable conditions, which is very steady in stable conditions.

**(a)**

**(b)**

**(c)**

The spectral analysis is also very important in chatter analysis. The corresponding spectral analysis results of vibration signals measured in partial experimental conditions are shown in Figure 12 (the sampling time is 0.5 s and the sampling region is near to the point in Figure 8). To see the spectrum clearly (spectral lines are very dense due to low rotating speed), the region from 1800 Hz to 2400 Hz is adopted to observe the milling stability in spectral analysis, the higher magnitude in frequency spectral concentrate on the region. From Figure 12, it can be seen that the main frequencies are all spindle rotating frequency (recorded as “★”) and tooth passing frequency (recorded as “×”) in stable conditions, while the chatter frequency (recorded as “∘”) becomes prominent in unstable conditions. In unstable conditions (chatter), the phase and chatter frequency do not coincide with spindle rotating frequency; the unstable types belong to Hopf bifurcation.

**(a)**

**(b)**

**(c)**

From Figures 10–12, it can be seen that the process damping increases the stability at low speed significantly and the method to predict stability limit is valid. For tool 1, the effect of process damping lasts until the spindle speed is up to 750 r/min (cutting speed m/min). For tool 2, stable regions are expanded remarkably relative to tool 1. The effects of process damping increase due to antivibration clearance angle, which lasts until the spindle speed exceeds 1200 r/min (cutting speed m/min).

Figure 13 shows the comparison of machined surface using two types of tool in the same condition (1000 r/min, 7 mm). In Figure 13(a), a lot of inclined vibration marks are left on the machined surface when using tool 1, which is the symbol of the occurrence of regenerative chatter. Whereas in Figure 13(b) the surface is very smooth and only some straight marks are left (it indicates the slight imposed vibrations) when using tool 2, antivibration clearance angles have significant effects on suppressing chatter.

**(a) Tool 1**

**(b) Tool 2**

Figure 14 shows the dynamic chip thickness, , calculated in different working conditions adopted in experiments, which can indicate the criterion provided in Section 3 is valid.

**(a)**

**(b)**

**(c)**

The comparative analysis of computational and experimental results shows that the nonlinear model considering process damping can predict stability limits at low speed in milling of thin-walled workpiece made of typical titanium alloy accurately, including the working conditions using tool with antivibration clearance angle. For tool 1, process damping can expand the stable regions significantly; limits can be improved from 2.5 mm to over 7 mm at 500 r/min (cutting speed m/min). The tool 2 can increase the process damping further and expand the speed range of process damping effect; the laws presented by (1) are proved fully.

#### 5. Discussions and Conclusions

Chatter problem in milling of titanium alloy, especially for thin-walled structural workpiece, is one of the main restrictions to the quality and efficiency of aeronautical manufacture. In the machining of titanium alloy, quite low cutting speed is always adopted for ensuring the tool life. Stability limits predicted by regular model ignoring process damping are much lower than actual limits, which may lead to terrible machining efficiency. In the paper, one milling dynamic model for thin-walled workpiece considering process damping is presented to predict stability limit. The computations and experiments indicate that the nonlinear model could provide accurate results. Besides the antivibration clearance angle designed could increase process damping to suppress chatter and improve machining surface quality effectively. Furthermore, the antivibration angle can expand speed range of process damping effects significantly.

However, it must be pointed out that the calculation model of process damping proposed in paper is a simplification of a much complicated interaction between the cutting edge, the flank face, and the workpiece. The interaction is highly dependent on the edge condition and in particular on tool wear. In the current work, only sharp cutters were utilized.

In summary, although lots of simplifications have been adopted in computational model, the importance and effect of model can be also shown in paper. In future work, more complex tool geometry will be tried out, and the effect of honed radius and flank wear will be investigated.

#### Appendix

#### A. The Calculation of in (6) and in (2)

We adsorbed some ideas in [23] and provide one useful method to calculate dynamic chip thickness and indentation area in milling of thin-walled workpiece. It considers the nonlinear effects of engagement and disengagement from the workpiece and could describe the surface wave caused by vibration of workpiece accurately.

##### A.1. Determined Function

As shown in Figure 15, the tool is assumed to be rotating anticlockwise, the workpiece surface is updated as material is removed, and the chip thickness is evaluated at each discretized point on the cutting edge. For downmilling, the workpiece surface is divided into two regions: the cutting arc surface and finished surfaces, represented by arrays of coordinates, , , as shown in Figure 15. Here, is the entry angle, is the exit angle, and is radial cutting depth, for downmilling:

When angular position , the tool is out of contact with workpiece, and the cutting force and process damping force are zero. According to common analytic methods, when , the cutting force is also regarded as zero [4]. The unit step function which determines whether the tooth is in or out of the cutting, (mentioned in Section 3.1), is always defined as 1 when ; otherwise is 0. However, in practice, the finished surface is affected probably by a tooth at the instant it passes through. Therefore, the scope of should be expanded. In the paper, the definition of is where is the angular position of tooth bottom (axial level ), which is shown in (5), and is axial cutting depth. When , the calculation of cutting force and process damping force at every axial level should continue.

Additionally, we should mention that the calculation of indentation area should also stop when corresponding angular position at th time point, , is smaller than entry angle.

##### A.2. Calculation of Instantaneous Chip Thickness

As mentioned in Section 3.1, the calculation of instantaneous chip thickness is of great importance in the whole computation process, which is required to describe the wave in workpiece surface accurately. As shown in Figure 16, two surface arrays are stored, one of the surface left by the previous tooth and one being created by the current tooth. Figure 16 shows how the surface is updated. At each time interval, a new point is added to the current surface at the angular position of the tooth,. If the tooth is cutting, the instantaneous position of tooth edge and vibration displacement of workpiece are used; otherwise a point is found on the previous surface at the tooth angle.

**(a) Tooth submersed in workpiece**

**(b) Tooth separated from workpiece**

We provide one method to represent the wave in workpiece surface. As shown in Figure 16, this method generates a point on the surface at each time interval corresponding to the instantaneous position of the cutting edge and vibration displacement of workpiece.

In order to determine if the tooth is submersed in the workpiece, the intersection point of the tooth and the previous surface is calculated. The point on the previous surface, , is found which has an angular position, , immediately preceding the angular position of the tooth, ; the intersection point is then found by linear interpolation.

At angular position, , the coordinate values of cutting point on workpiece surface are

Firstly, the radial distance, , is obtained by linear interpolation of the radial distance between points and on the previous surface: The coordinate values of intersection point are

As shown in Figure 16(a), it indicates that the tool is cutting workpiece if ; otherwise, it indicates the tool has separated from workpiece as shown in Figure 16(b). Thus, the dynamic chip thickness can be obtained:

Then the coordinate values of surface wave caused by current tooth path can also be generated. If the tooth is cutting workpiece, the coordinate values of surface wave, , are taken as the current coordinate values of cutting point, . Otherwise, the coordinate value of wave is taken as the value of :

##### A.3. Calculation of Coordinates of Surface Wave in Radial Direction

For calculating indentation area and process damping force in accordance with (2)–(4), the radial coordinate values of surface wave, , should be also obtained. If tool is cutting workpiece, is equal to the radial vibration displacement of workpiece, which can be obtained using and angular position , so does the corresponding radial vibration speed. Otherwise, if tool has separated from workpiece, is taken as the ones caused by the previous tooth pass, which is also generated by linear interpolation. In (A.8), is radial coordinate values of wave caused by previous pass:

##### A.4. Feeding the Workpiece

The workpiece is fed into the tool at the rate [mm/tooth]. With a spindle speed, , a time interval of , and number of teeth, the motion of the workpiece along the -axis during a single time step is

The components of all points in all surface arrays are incremented by at each time step before the surface is updated due to the teeth cutting. Correspondingly, and in (A.8) are also incremented by .

##### A.5. Special Case

When tooth is just entering the arc surface as shown in Figure 17, tooth may first strike the upper uncut workpiece surface before cutting the surface cut by previous teeth. The chip thickness is found using (A.6) with the point found as intersection between the tooth edge and the line () as shown in (A.10). Then chip thickness is obtained using the above methods:

#### Conflict of Interests

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

#### Acknowledgment

The research is sponsored by the National Natural Science Funds (51475234).