Next Article in Journal
Research on Wave Energy Harvesting Technology of Annular Triboelectric Nanogenerator Based on Multi-Electrode Structure
Next Article in Special Issue
Force Prediction and Material Removal Mechanism Analysis of Milling SiCp/2009Al
Previous Article in Journal
Demonstration of 144-Gbps Photonics-Assisted THz Wireless Transmission at 500 GHz Enabled by Joint DBN Equalizer
Previous Article in Special Issue
Study on the Formation Mechanism of Cutting Dead Metal Zone for Turning AISI4340 with Different Chamfering Tools
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Ti-6Al-4V Milling Force Prediction Model Based on the Taylor Factor Model and Microstructure Evolution of the Milling Surface

1
School of Mechanical and Automotive Engineering, Shanghai University of Engineering Science, Shanghai 201620, China
2
School of Mechanical Engineering, Georgia Institute of Technology, North Ave NW, Atlanta, GA 30332, USA
*
Author to whom correspondence should be addressed.
Micromachines 2022, 13(10), 1618; https://doi.org/10.3390/mi13101618
Submission received: 10 September 2022 / Revised: 23 September 2022 / Accepted: 25 September 2022 / Published: 27 September 2022
(This article belongs to the Special Issue Advances in Micro-Milling)

Abstract

:
In this paper, a milling force prediction model considering the Taylor factor is established, and the Ti-6Al-4V milling force predicted by the model under different milling parameters is presented. In the study, the milling experiment of Ti-6Al-4V was carried out, the milling force was collected by the dynamometer, and the microstructure evolution of the milling surface before and after milling was observed by EBSD. Through the comparative analysis of the experimental results and the model prediction results, the reliability of the prediction model proposed in this study was verified, and the influences of the milling parameters on the milling force were further analyzed. Finally, based on the EBSD observation results, the effects of the milling parameters on the microstructure evolution of the milling surface were studied. The results show that both the tangential milling force and normal milling force increase with the increase in the milling depth and feed rate. Among the milling parameters selected in this study, the milling depth has the greatest influence on the milling force. The average errors of the tangential milling force and normal milling force predicted by the milling force model are less than 10%, indicating that the milling force prediction model established in this paper considering Taylor factor is suitable for the prediction of the Ti-6Al-4V milling force. With the change in the milling parameters, the grain structure, grain size, grain boundary distribution, phase distribution, and micro-texture of the material surface change to varying degrees, and the plastic deformation of the milling surface is largely coordinated by the slip.

1. Introduction

Titanium is an allotrope, and α -titanium has an HCP (hexagonal close-packed) structure, while β -titanium has a BCC (body-centered cubic) structure. Titanium alloy is an alloy based on titanium and is composed of other elements. Common titanium alloys are divided into the following three categories: α -titanium alloy, β -titanium alloy, and ( α + β ) titanium alloy [1]. Ti-6Al-4V is a typical ( α + β ) titanium alloy, and Ti-6Al-4V is also one of the materials widely used in the field of aerospace engineering [2,3,4]. Because Ti-6Al-4V has the physical properties of a high hardness, high strength, and low thermal conductivity, it produces a very high temperature in the cutting process, and the material is easily attached to the tool surface. At present, there are certain common problems regarding the machining process of titanium alloy, such as a low machining efficiency and serious tool wear. In the field of aerospace engineering, the material removal per unit time in the actual processing of titanium alloy workpieces is much less than that of other materials. Therefore, the question of how to select a reasonable processing method and improve the machining efficiency and quality has become an important research topic [5].
The cutting force is a very important parameter in the machining process, which has direct impacts on the stability and efficiency of the machining process and further affects the quality of the machined surface and tool wear. Over the last century, many scholars have carried out a great deal of research on the prediction of the cutting force, and their research methods can be divided into the empirical model, finite element model, and analytical model [6]. Based on the machining principle, the analytical method uses geometric and mechanical methods to analyze the cutting process, which can reflect the real machining state and predict the results accurately. The classical analytical model includes the right angle cutting model proposed by Merchant [7] and Oxley [8], which assumes that the shear plane is located at the place where the shear energy is minimal and calculates the shear angle according to this assumption. Based on the classical model, Endres et al. [9] and Moufki et al. [10] introduced the ploughing force model and force–thermal coupling model, respectively, which further improved the prediction accuracy of the analytical model. By using the combination of the analytical model and finite element simulation, Yang [11] established a constitutive model considering the large deformation, large strain rate, and high temperature characteristics and predicted the cutting force and cutting deformation in the milling process of titanium alloy. Das et al. [12] proposed a new method of grey relational analysis of the Ti6Al4V machining parameters based on Taguchi. The effects of the cutting speed, tool feed, and cutting depth on the matching response were studied. Kim [13] carried out the milling force and thermal analysis of Ti-6Al-4V by the finite element method, predicted the temperature distribution, and evaluated the surface roughness, tool wear, and cutting force. Compared with the results of traditional face milling, laser-assisted face milling improves the surface roughness. In addition, in terms of tool damage, laser-assisted face milling has been demonstrated to be more effective than traditional face milling. Perez et al. [14] proposed a Taylor model based on oblique cutting so as to quantify the effect of the crystal structure on the shear strength. The microstructure, crystal texture, and particle morphology of different samples were analyzed by EBSD, and the tool geometry, tool position, and laser scanning path were taken into account. The analysis showed that there is a significant interaction between the plane direction of the shear band and the grain orientation on the spindle, and there is a strong correlation between the cutting force and the prediction model. Based on the J–C (Johnson–Cook) flow stress model and mechanical–thermal coupling model, Ji [15] established a cutting force prediction model considering the cutting parameters, tool parameters, material properties, and surface lubrication parameters. Through the coupled iterative calculation method, this model iteratively determines the cutting force and cutting heat in the cutting process and is better than other models in terms of the scope of its application and prediction accuracy.
The above studies on the prediction models of the cutting force do not account for the microstructure of the material. In the process of material deformation, the evolution of the microstructure usually has an important influence on the physical properties of the material [16]. Many scholars use EBSD (electron backscattering diffraction) to analyze the microstructure changes in materials (including the microstructure, grain size, slip, twin, etc.). Based on the new method of external friction theory, Tong et al. [17] analyzed the friction-reducing and anti-wear mechanisms of micro-texture cutting tools. The optimum area ratio of the micro-texture in the contact zone between the chip and tool was studied when the cutting force was considered. Su et al. [18] used PCD tools to study the effect of ultrasonic vibration-assisted machining (UVAM) on the machinability of the Ti6Al4V titanium alloy and its mechanism. The cutting force, machined surface, surface bonding force, and wear morphology were analyzed. The results showed that UVAM can effectively improve the cutting performance. Wu et al. [19] carried out a comparative experiment on the dry milling of the Ti-6Al-4V alloy with micro-texture cutters and common cutters. The results showed that the micro-texture tool can reduce the cutting force, cutting temperature, and power consumption by about 15%, 10%, and 5%, respectively. The results showed that the cutting radial width, cutting speed, and axial cutting depth have obvious influences on the dry grinding of the Ti-6Al-4V alloy, and the mechanism for improving the performance of the microstructure tool is mainly its self-lubrication. Dabwan et al. [20] studied the effects of the feed speed, radial cutting depth, and cutting speed on the surface roughness, cutting force, microhardness, microstructure, chip morphology, and surface morphology of Ti6Al4V. It was found that different orientations have different effects on the machined surface. The results showed that the EBM parts machined in the TLP direction have a good surface quality and surface integrity. To sum up, the microstructure of the material has an important influence on the mechanical properties, and it also affects the cutting force and cutting temperature in the material processing.
In the above studies, the effects of the grain orientation on the properties of the materials were not considered. In this paper, the Taylor factor model is introduced, and the Taylor factor is taken as a variable. Based on the milling force prediction model proposed by Ji [15], an improved milling force prediction model considering the effect of the grain orientation is established. Combined with the milling parameters of a small titanium alloy workpiece in actual production, the orthogonal grinding experiment of Ti-6Al-4V was established. The milling force in the milling process was measured and recorded using a Kistler three-dimensional dynamometer (model: 9272), and the microstructure of the milling surface was observed using an EBSD detector (model: Nordly-Max-2). Through the analysis of the experimental results, the applicability of the prediction model established in this paper was verified, and the effects of the milling parameters on the milling force and milling surface microstructure were further studied.

2. Materials and Methods

2.1. Materials and Experimental Equipments

The titanium alloy bars used in this study are obtained by rolling and machined into cuboid samples of 15 × 15 × 30 mm size by wire cutting. In order to ensure that all the milling samples have the same surface microstructure before machining, all wire cutting samples are obtained from the same Ti-6Al-4V bar. The samples and the experimental equipment used in the milling experiment are shown in Figure 1. The milling force in the milling process is transmitted from the dynamometer to the portable USB signal collector by means of electrical signals. The electrical signals are processed by the charge signal amplifier and then transmitted to the PC.
The components of Ti-6Al-4V (data provided by the supplier) are shown in Table 1. The milling cutter used in the experiment was a YG6X cemented carbide end milling cutter with 4 edges, 10 mm diameter, 9 ° front angle, 8 ° back angle, and 35 ° helix angle. The milling method is climb milling, and the milling process is dry milling (without adding coolant). In order to verify the accuracy of the prediction model, dry milling is used in the experimental design, which is consistent with the mechanical–thermal coupling mode in the prediction model. Before the milling experiment, all milling surfaces are finished with the same face milling cutters.
The material removal rate per unit time is low in the actual milling process of small titanium alloy workpieces. In this study, with reference to the range of processing parameters of Ti-6Al-4V small workpieces, an orthogonal experimental table was designed, as shown in Table 2. In order to investigate the evolution of the micro-texture after milling, EBSD samples of Ti-6Al-4V were prepared after mechanical polishing. In mechanical polishing, the sample surface is polished with 120 μ m sandpaper, 3 μ m and 0.3 μ m polishing cloth until several shallow scratches are formed. The schematic diagram of the EBSD observation is shown in Figure 2.

2.2. Prediction Model of the Milling Force

2.2.1. Taylor Factor Model

Schmidt’s law is widely used to study the slip and twin mechanisms of materials. Schmidt’s law holds that there is a critical value of the shear stress in the slip system and twin system, and the slip system and twin system that reach the critical value start first. The Schmidt factor is calculated by the following formula:
S c h m i d t   f a c t o r = c o s ( Φ ) · c o s ( Ω )
Here, Φ represents the angle between the slip plane and twin plane and macro-stress, and Ω represents the angle between the slip direction, twin direction, and macro-stress. The greater the Schmidt factor is, the greater the shear stress is, and the material is more likely to slip or twin. According to von Mises law [21], polycrystals usually require at least five independent or systematic systems during plastic deformation. In polycrystal deformation, the complex distribution of the microscopic stress field leads to the anisotropy of the flow stress, and there is a difference between the local stress and macroscopic stress of each grain. The expression of the flow stress is as follows:
σ s = τ c c o s ( Φ ) · c o s ( Ω ) = m c · τ c
Here, τ c represents the critical shear stress of the activated slip system and twin system, and m c is the reciprocal of the maximum Schmidt factor. For single crystals, the normal stress σ = M · τ c , and M is called the Taylor factor [22]. In polycrystals, the Taylor factor changes with the change in the crystal structure and force direction of each grain. Mao et al. [23] and Zhao [24] assumed that the Taylor factor is still applicable in a high-strain-rate machining environment and multi-slip system, and the Taylor factor of a polycrystal is the linear superposition of the single crystal Taylor factor; thus, they proposed the Taylor factor model of the FCC crystal structure and BCC crystal structure. The validity of the Taylor factor model is verified by micro-grinding experiments on 3J33b alloy steel. The single crystal Taylor factor model of the FCC crystal structure and BCC crystal structure is as follows:
M F C C , B C C = i = 1 N 1 c o s ( Φ i + γ 1 ) · c o s ( Ω i + γ 1 )
In the above formula, N represents the number of slip systems in the FCC or BCC crystal structure, and γ 1 and γ 1 are correction angles of 0 ° and 10 ° , respectively. By regarding the Taylor factor model of the polycrystal as the linear superposition of the single crystal, the Taylor factor model of the polycrystal is obtained as follows:
M P o l y = j = 1 N 1 f j · M j F C C , B C C
In the formula, f j represents the distribution function value of the grain orientation j, and N 1 represents the number of the grain orientation.
In this study, on the basis of the above research, the Taylor factor model of the HCP crystal structure for Ti-6Al-4V is established. Ti-6Al-4V is a typical ( α + β ) titanium alloy, in which α -Ti has an HCP structure, while β -Ti has a BCC crystal structure. Figure 3 shows the crystal structure diagrams of HCP structure (Figure 3a) and BCC structure (Figure 3b).
In the HCP crystal structure, the slip types can be divided into three types: a slip, c slip, and a + c slip. Among them, a slips and a + c slips are easily opened during plastic deformation. Because the a slip can only provide four independent slip systems, which cannot meet the von Mises criterion, the HCP crystal structure also needs to turn on the a + c slip and twin in order to meet the five independent slip system conditions for plastic deformation. The common slip systems of HCP metals are listed in Table 3 [25].
Although the a + c slip can coordinate the deformation in the direction of the c-axis to some extent, the a + c slip is difficult to open and unstable in HCP metal. It is inevitable that the twinning mode must be activated to coordinate the deformation on the c-axis [26,27]. The common twin systems of HCP materials are shown in Table 4 [28]. Many studies have shown that the c/a axial ratio ( γ c / a ) has an effect on the selection of deformation twins in HCP materials [29,30,31], and the twins that can be activated in HCP materials with a different γ c / a are different. The γ c / a of titanium alloy is 1.587, which is less than the ideal axial ratio (1.633). In the five common twinning modes of HCP materials, the { 10 1 ¯ 2 } 10 1 ¯ 1 ¯ extension twins and { 11 2 ¯ 2 } 11 2 ¯ 3 ¯ contraction twins mainly occur in titanium alloys.
In the BCC crystal structure, the number of independent slip systems according to von Mises criterion can satisfy the uniform plastic deformation, so that only the role of the slip systems is generally considered. The slip plane families of BCC materials are { 101 } , and the slip direction families are 111 . All slip systems are shown in Table 5 [32].
Based on the plastic deformation mechanism, the Taylor factor model considering the slip and twin in α -Ti and slip in β -Ti is established in order to describe the effect of the grain orientation of Ti-6Al-4V during plastic deformation. The Taylor factor model of the HCP structure of α -Ti constructed in this study is as follows:
M H C P = i = 1 N 1 c o s ( Φ i + γ 1 ) · c o s ( Ω i + γ 1 ) + j = 1 N 1 1 cos ( Φ j + α 1 ) · c o s ( Ω j + α 1 )
In the formula, N represents the number of slip systems, N 1 represents the number of twin systems, γ 1 and γ 1 are the correction angles, which are 0° and 10°, respectively, [24], and α 1 and α 1 are deviation angles of 0° and 3°, respectively [33]. The Taylor factor model of a single crystal with a BCC structure refers to Formula (3). In this study, the Taylor factor model of polycrystals with an HCP structure and a BCC structure is regarded as the linear superposition of the Taylor factors of single crystals, as follows:
M P o l y = i = 1 N 3 f i · M i H C P , B C C  
In the formula, f i represents the distribution function value of the grain orientation, and N 3 refers to the number of the grain orientation. The HCP structure and BCC structure belong to α -Ti and β -Ti, respectively, in Ti-6Al-4V. Thus, the Taylor factor model of the two-phase materials was constructed in this study as follows:
M = f α T i · M H C P + f β T i · M B C C  
Here, f α T i and f β T i are the proportions of α -Ti and β -Ti in Ti-6Al-4V, respectively.

2.2.2. Milling Force Prediction Model Considering the Taylor Factor

When establishing the analytical model of the milling force, the milling process needs to be simplified using the right angle cutting model or bevel cutting model. On the basis of Oxley’s prediction model, Ji [15] established a cutting force prediction model with a wider range of parameters, considering the thermo-mechanical coupling effects of the material and tool parameters. In this study, the cutting force prediction model considering the Taylor factor was established by combining the cutting force prediction model proposed by Ji [15] with the Johnson–Cook flow stress model, considering the influence of the microstructure.
In the prediction model of the cutting force, the main sources of the cutting force are the chip-forming force and ploughing force. The chip-forming force analytical model and ploughing force model [34] used in this study are shown in Figure 4. According to the analytical model of the chip-forming force, the chip-forming force ( F C h i p ) can be divided into the cutting force along the cutting direction ( F C u t t i n g ) and the cutting force perpendicular to the machined surface ( F T h r u s t ). The components F C u t t i n g and F T h r u s t can be calculated by the following formula:
{ F C u t t i n g = R · c o s ( λ α ) F T h r u s t = R · s i n ( λ α )  
R = k A B t 1 w s i n ϕ · c o s θ
In the formula, R is the resultant force of the normal force and friction force, θ is the angle between R and the shear plane A-B, α is the tool rake angle, ϕ is the shear angle, λ is the angle of friction, k A B is the flow stress on the shear plane A-B, t 1 is the undeformed chip thickness, and w is the cutting width. Among these, k A B must be calculated by the Johnson–Cook flow stress constitutive equation. The milling process of titanium alloy is a process comprising a large strain, high strain rate, and high temperature, which is consistent with the scope of the application of the Johnson–Cook model. According to the research results of Perez et al. [14], the Zerilli–Armstrong model can also effectively describe the milling process of titanium alloy.
When the material is deformed during processing, the microstructure of the material also evolves, including elements such as the phase distribution, grain size, and grain orientation. The effect of the grain size on the flow stress during machining can be quantified by Hall–Petch formula [35]. According to the Hall–Petch formula, the flow stress is inversely proportional to the square root of the grain size. On the other hand, the Taylor factor model established in this study can quantitatively describe the effects of the phase distribution and grain orientation on the plastic deformation process of materials. Combined with the Hall–Petch formula and the influence model of the microstructure on the flow stress [36,37,38], the Johnson–Cook flow stress model, considering the initial microstructure of the material, is used to calculate the flow stress in the cutting process as follows:
σ = ( A + B ε p n ) ( 1 + C l n ε p ˙ ε 0 ˙ ) ( 1 ( T T w T m T w ) m   ) + M α G b ρ 1 + K H P D d  
K H P = M · τ b · 4 G · b ( 1 ν ) π
In the formula, A, B, C, n, and m are all material constants in the Johnson–Cook model, ε p is the effective plastic strain, ε p ˙ is the effective strain rate, and ε 0 ˙ is the reference strain rate (usually 1 s 1 ). T is the current temperature, T w is the ambient temperature, and T m is the material melting point. M is the Taylor factor model, G is the material shear modulus, b is the Burgers vector, ρ 1 is the dislocation density, ν is Poisson’s ratio, τ b = 0.057 G , α is the dislocation density interaction coefficient ( α = 2 ( 1 υ ) / ( 2 υ ) ) [37], and D d is the diameter of the grain. According to the Johnson–Cook flow stress model, the flow stress of the shear plane AB ( k A B ) can be expressed as:
k A B = σ A B 3  
In this study, the analytical model of the ploughing force proposed by Waldorf et al. [34] was used to predict the ploughing force caused by the fillet of the knife tip. The ploughing force can also be divided into two components, including the ploughing force P C u t t i n g along the cutting direction and the ploughing force P T h r u s t perpendicular to the machined surface, as follows:
P C u t t i n g = k A B · w · [ c o s ( 2 η p l o w ) c o s ( ϕ γ p l o w + η p l o w ) + ( 1 + 2 θ f a n + 2 γ p l o w + s i n ( 2 η p l o w ) ) · s i n ( ϕ γ p l o w + η p l o w ) ] · C A  
P T h r u s t = k A B · w · [ ( 1 + 2 θ f a n + 2 γ p l o w + s i n ( 2 η p l o w ) ) · c o s ( ϕ γ p l o w + η p l o w ) c o s ( 2 η p l o w ) s i n ( ϕ γ p l o w + η p l o w ) ] · C A
C A = R f a n s i n η p l o w
The fan-shaped angles θ f a n , γ p l o w , ρ p r o w , and η p l o w can be solved according to the geometric relation and friction relation [34]. Here, r is the radius of the tool edge and ϕ is the shear angle. R f a n is solved by the following equation:
R f a n = ( r · t a n ( π + 2 α 4 ) + 2 · R f a n · s i n ( ρ p r o w ) t a n ( π + 2 α 2 ) ) 2 + 2 [ R f a n · s i n ( ρ p r o w ) ] 2 · s i n ( η p l o w )  
After the calculation of the chip-forming force and ploughing force is completed, the total cutting force in the cutting process can be expressed as:
{ F t = F C u t t i n g + P C u t t i n g F n = F T h r u s t + P T h r u s t  

3. Results and Discussion

3.1. Analysis of the Milling Force and Verification of the Prediction Model

The physical parameters of the workpiece material Ti-6Al-4V are shown in Table 6. The Johnson–Cook flow stress model, considering the microstructure of materials, is required in the cutting force prediction model. The Johnson–Cook material parameters of Ti-6Al-4V are shown in Table 7 [39]. The material parameters of the tool are shown in Table 8, below.
Microstructure maps of the material surface before milling are shown in Figure 5. From the figure, we can observe that the grains on the material surface are mainly equiaxed grains, the average grain size is about 4.5 μ m , the proportion of the α phase is about 98.6%, and the proportion of the β phase is about 1.4%. Figure 5d shows the KAM map of the surface. According to the KAM diagram, the degree of deformation can be analyzed intuitively. The higher the KAM value is, the greater the degree of plastic deformation is.
In this study, the geometric dislocation density ( ρ G N D ) was calculated according to the KAM value, which is expressed by the following formula:
{ ρ G N D = 2 K A M a v g μ b K A M a v g = exp ( 1 N 1 i l n K A M i )  
In the formula, μ is the scanning step, b is the Burgers vector, K A M a v g is the average KAM value of the selected area, K A M i represents the KAM value at point i , and N represents the number of points in the test area.
In this study, the milling force was collected by the Kistler-9272 three-dimensional dynamometer, which can collect and record FX, FY, FZ, and the corresponding torque. In order to eliminate the influence of the noise signal caused by the spindle vibration, the collected data must be transformed by FFT (fast Fourier transform). Figure 6 shows the signal curve of the normal milling force in the first group of experiments before and after low frequency filtering, with a fixed frequency of 42 Hz.
The experimental parameters and results under the condition of dry milling are shown in Table 9.
In order to intuitively analyze the effect of each milling parameter on the actual milling force F t and F n , the range analysis of the milling force was carried out in this study. In the analysis, the K value of each level represents the sum of the experimental results of the selected factors at the current level. The K-avg value indicates the average value of K corresponding to the selected factor at the current level, and the R value represents the range of each factor, which is obtained by subtracting the minimum value from the maximum value in the K-avg. The higher the R value is, the more obvious the influence of this factor on the result is. The range analysis tables of the milling force measured in the experiment are shown in Table 10 and Table 11, respectively. By analyzing the results of the orthogonal experiment, we can see that, for the tangential component of the cutting force F t , the ranges of the linear speed, feed rate, and cutting depth are 21.90, 9.73, and 22.20, respectively. The ranges of the linear velocity, feed rate, and cutting depth of the normal component of the cutting force F n are 3.10, 5.40, and 11.13, respectively. Therefore, the influence of the cutting depth on F t is the greatest, the influence of the linear velocity is the second greatest, and the influence of the feed rate is the least. For F n , the influence of the cutting depth is the greatest, the influence of the feed rate is the second greatest, and the influence of the linear velocity is the least. The above analysis shows that the milling depth has the greatest influence on the overall milling force in the process of milling.
In the actual milling process, because the surface of the machined material is not completely horizontal, there is a difference between the actual machining depth and the theoretical value of each sample in the milling process. Moreover, the tool wears during the machining process, which changes the geometric parameters of the tool (rake angle, edge radius, etc.) and then affects the milling force. The milling force prediction model used in this study is based on the assumption that the milling parameters (cutting speed, feed rate, cutting depth) and tool geometric parameters do not change in the machining process; thus, there is an error between the milling force prediction model and the actual experimental results. Generally, when the actual milling depth is greater than the parameters set in the model, or when the tool is worn, the actual milling force will be greater than the predicted value of the model. When the actual milling depth is less than the parameters in the model, or a new tool is replaced, the actual milling force will be less than the predicted value of the model.
Figure 7 shows that, within the given range of the milling parameters, both the tangential milling force and the normal milling force show an overall increasing trend with the increase in the milling depth and feed rate. The tangential milling force decreases with the increase in the milling speed, and the normal milling force first decreases and then increases slightly with the increase in the milling speed. In the cutting direction, both the chip-forming force and the ploughing force increase with the increase in the feed rate. This is due to the increase in the feed, which is the thickness of the cut, and, therefore, the increase in the cutting force. However, when the feed rate increases, the deformation coefficient of the chip decreases, and the friction coefficient also decreases, thereby reducing the friction force. Due to these effects, the increasing speed of the cutting force is not proportional to the increasing speed of the feed. When the cutting speed increases, the cutting force in the cutting direction decreases. As the cutting speed increases, the cutting temperature also increases, and the strength and hardness of the processed metal decrease, resulting in a decrease in the cutting force. The cutting speed and the feed amount have almost no effect on the cutting force in the direction perpendicular to the surface of the workpiece.
Many researchers have found that the milling force decreases with the increase in the milling speed in their experiments. The research results of Klocke et al. [40] indicate that, on the one hand, intermittent cutting increases with the increase in the milling speed, and the hardening of the material will lead to the increase in the deformation resistance in the area affected by the shear force. On the other hand, the increase in the milling speed leads to the increase in the milling temperature and the softening of the materials, which leads to the decrease in the milling force. The influence of the milling speed on the milling force is the result of the comprehensive action of the above two factors. The increase in the feed rate and milling depth affect the geometry of the milling area, resulting in the increase in the cross-sectional area of the milling area and the amount of material removal per unit time, which leads to the change in the milling force. Many scholars have studied the milling force of Ti-6Al-4V, and most of the analytical results concerning the influence of the milling parameters on the milling force correspond to the results of this study [41,42].
Figure 8 shows a comparison of all the experimental results with those predicted by the model. According to the analysis of Figure 8, it is found that there is a small error between the predicted and experimental values of the cutting force in both the tangential and normal directions. The reason may be that the change in the milling force caused by tool wear and the change in the milling depth caused by the horizontal degree of the machined surface are not considered in theoretical research.
Table 12 shows the relative errors of the milling forces and predictions for each group. From the table, we can see that the relative error range of the tangential milling force obtained from the prediction model in this study is 3.4–20.7%, and the average error is 8.3%. And the relative error range of the normal milling force is 3.7–14.3%, and the average error is 8.5%.
From Table 12, we can see that when the milling speed and feed rate are at a high level but the milling depth is small, the relative error between the predicted results and the experimental results increases (14.2%, 20.7%, etc.). This may be caused by the following factors: (1) Firstly, when the milling depth decreases, the influence of the chip-forming force on the milling force decreases, which leads to the greater influence of the ploughing force on the milling force. The ploughing force is also affected by the material hardness, but the material hardness is not considered in the prediction model established in this study, which may lead to the increase in the relative error. (2) Secondly, when the milling depth is very low, the change in the actual milling depth caused by the uneven surface of the material has a great influence on the results. The prediction model cannot simulate the change in the milling depth in the machining process, which also leads to the increase in the relative error.
The average prediction errors obtained by Ji [15] and Li et al. [43] through the cutting force prediction model, without considering the Taylor factor, are 13.4% and 10.2%, respectively. This indicates that the prediction accuracy of the milling force prediction model considering the Taylor factor proposed in this study is better than the original model.

3.2. Effects of the Milling Parameters on the Microstructure Evolution of the Milling Surface

3.2.1. Microstructure Shape and Grain Size

In the milling process, the change in the milling parameters lead to changes in the temperature, strain, and strain rate of the machined surface, thus affecting the grain structure of the machined surface. Under the action of milling, some of the grains in the cutting are distorted along the milling direction of the tool. The grain orientation distribution maps and grain size statistics of Ti-6Al-4V under different milling parameters are shown in Figure 9 and Figure 10 (subfigures 1–9 correspond to 1–9 groups of samples in the orthogonal experiment table, respectively). The research object is the milling surface of the sample, and the EBSD calibration area is 60 μ m × 50 μ m , which characterizes the local microstructure of the milling surface of the sample. From Figure 9, we can see that the grains in the surface of the milled materials are refined, and most of them are equiaxed grains.
By analyzing the statistical diagram of the grain size, it is found that the grain size decreases with the increase in the milling speed and milling depth, and the effect of the milling depth on the grain size is the most significant. Further analysis shows that the grain size of the subsurface layer is smaller than that of the surface layer, and the grain size of the machined surface shows a gradient change in the depth direction. The change in the grain size is closely related to the temperature field and strain rate field distribution of the material surface. On the machined surface, which is in direct contact with the tool, the plastic deformation caused by the friction and shear between the tool and the material and the high temperature caused by milling has obvious effects on the grain size. Due to the high temperature produced by the machined surface during milling, the refined grains become larger under the action of the high temperature. Meanwhile, most of the heat generated in the milling process is carried away by the chip, and in the subsurface layer of the material surface, the high strain rate plays a dominant role in shaping the grain size. When the milling speed increases from 60 m/min to 120 m/min, the milling temperature of the material does not change obviously in this range of the milling speed because the contact time between the tool and the material becomes shorter and more chips reduce the heat generated in the milling. Therefore, in this study, with the increase in the milling speed, the effect of plastic deformation on the grain size plays a major role, and the grain size decreases with the increase in the milling speed. Similarly, with the increase in the milling depth, the plastic deformation increases significantly, and the strain rate field of the machined surface layer is transferred to the subsurface layer, causing grain refinement.
The above results show that the change in the grain size follows the law of grain refinement on the micron scale, as follows: with the enhancement of the plastic deformation, the grain is stretched, distorted, and broken, and the grain refinement is enhanced, while with the increase in the temperature, the grain easily recrystallizes and grows, and the grain refinement is weakened. Moreover, there are no twins in most grains and flat twins in a few grains, indicating that the deformation twinning behavior occurs in the grains to coordinate the plastic deformation. In the milling process of Ti-6Al-4V, dislocation slip plays a dominant role in the plastic deformation of the surface, while the deformation twins coordinate only a small part of the plastic deformation process, and the grains of almost all the samples are refined to varying degrees, which can be considered as the result of the complex interaction between the twins and dislocations.

3.2.2. Grain Boundary Distribution and Phase Distribution

The grain boundary distribution maps and phase distribution maps of Ti-6Al-4V under different milling parameters are shown in Figure 11 and Figure 12 (subfigures 1–9 correspond to 1–9 groups of samples in the orthogonal experiment table, respectively). According to the misorientation angle difference of the adjacent grains, the grain boundary is divided into the high-angle grain boundary (misorientation angle 15 ° , HAGB) and low-angle grain boundary (misorientation angle < 15 ° , LAGB). In Figure 11, the HAGB is represented by blue lines, the LAGB is represented by green lines, and the grain boundaries of the sub-grains are represented by red lines. Based on the analysis of the proportion of each grain boundary, it was found that the proportion of the LAGB increases with the increase in the milling depth. With the increase in the milling speed and feed rate, the proportion of LAGB first increases and then decreases. Further analysis showed that the influence of the milling depth on the LAGB distribution is the greatest.
Grain refinement caused by plastic deformation is characterized by the appearance of a large LAGB within the grains, while the grains undergoing dynamic recrystallization usually have an HAGB. The percentage of the LAGB can reflect the intensity of the dislocation activity. When the milling speed, feed rate, and milling depth increase, the amount of strain and the strain rate increase, and the percentage of the LAGB increases. This indicates that the proportion of refined grains caused by plastic deformation gradually increases, the slip system inside the grains is activated, and the dislocation slip is the main deformation mode. With further increases in the milling speed and feed, the percentage of the LAGB tends to decrease, which may be due to the shortened time given for dislocations to occur inside the grains due to the increase in the strain rate. The time given for dislocation slip to occur inside the grain is reduced, and the grain does not have sufficient time to open the slip system so as to coordinate the deformation. Moreover, with the increase in the strain amount and strain rate, the heat generated by the milling surface also increases, and the grains are more likely to undergo dynamic recrystallization under the action of the high temperature, resulting in an increase in the percentage of the HAGB. Some studies have also pointed out that with the increase in the degree of grain deformation, a large number of dislocations in the grain will offset each other and rearrange, which also reduces the percentage of the LAGB to a certain extent.
As mentioned above, Ti-6Al-4V is a two-phase ( α + β ) titanium alloy. The phase distributions of each group of the sample milling surface are shown in Figure 12, with the α phase in red and the β phase in green. The proportion of the α phase in the milling surface of each group is higher than 90%. Compared with the phase distribution of the surface of the material before milling, this indicates that the proportion of the α phase and β phase changed to varying degrees under the action of milling. With the increase in the milling speed and milling depth, the proportion of β phase increases, and the milling depth has the most significant effect on the beta phase distribution.
The cutting temperature and local stress are the main factors causing Ti-6Al-4V phase transformation. In the milling process, the material passes through two stages; the first stage is the heat generation stage, caused by instantaneous pressure, and the second stage is the cooling stage. In the first stage, there is a transition from the α phase to the β phase, while in the second stage, there is a transition from the β phase to the α phase. The milling method used in this study was dry milling, without the effect of a coolant, and the effect of the cooling stage is not obvious. The above analysis only considers the effect of the temperature on the phase transformation, but in the actual milling process, the severe plastic deformation and the mechanical load imposed on the material by the cutting tool also transform the α phase into the   β phase. Due to the fact that the   α -Ti has an HCP structure and the β -Ti has a BCC structure, the number of independent slip systems of α -Ti is much less than that of   β -Ti. In other words, β -Ti is more prone to plastic deformation than α -Ti under the same milling force, which can also explain the increased proportion of the β phase in the milling process.

3.2.3. Texture Evolution and Twinning Behavior

Figure 13 (subfigures 1–9 correspond to 1–9 groups of samples in the orthogonal experiment table, respectively) shows the pole figures and twin grain boundary contrast maps of the Ti-6Al-4V surface under different milling parameters. The orientation distribution of the α Ti grains perpendicular to the ND direction under different milling conditions is represented by the { 0001 } pole figure, { 10 1 ¯ 0 } pole figure, and { 1120 } pole figure. In the IPF (Figure 9), most of the areas are blue, indicating that the grain orientation distribution of the { 10 1 ¯ 0 } crystallographic planes in the vertical ND direction accounts for a large proportion. From the corresponding pole figure ( { 10 1 ¯ 0 } pole figure), the maximum intensity of most textures appears in the ND direction, indicating that there is a large number of prismatic textures on this cross section. Further analysis indicated that with the increase in the milling depth, the intensity of the prismatic surface texture increases (in Figure 9(3),(4),(7)). Moreover, with the increase in the milling depth, the texture of the { 10 1 ¯ 0 } crystallographic plane aggregates to the [ 0001 ] orientation, and the main texture distribution is in the range of ± 45 ° inclination from the [ 0001 ] orientation towards the TD and RD direction.
According to the analysis of the { 0001 } pole figure, there is no obvious basal texture of the milling surface, and the main component is the [ 0001 ] orientation, with an inclination of ± 45 ° ~ 60 ° towards TD and RD. When the feed rate reaches 240 mm / min , a small degree of low-intensity basal texture can be observed in the { 0001 } pole figure (in Figure 9(6), Figure 9(9)). The results show that the texture intensity of the { 0001 } pole figure decreases with the increase in the milling speed, and with the increase in the milling depth, the texture orientation of the { 0001 } crystallographic plane is inclined towards the RD and TD direction (in Figure 9(3),(4),(7)).
In addition, some textures with a high intensity are also observed in the { 1120 } pole figure—that is, the   [ 0001 ] orientation with a tilt of ± 20 ° ~ 30 ° towards TD and RD. With the change in the milling parameters, the texture distribution of the {1120} plane shows no obvious regularity, but the texture intensity of the {1120} plane decreases with the increase in the milling speed. The texture intensity also increases slightly with the increase in the milling depth (in Figure 9(3),(4)). Compared with the texture intensity before milling (Figure 5), it was found that the texture intensity changes in varying degrees with the change in the milling speed and milling depth. Among all the milling parameters, the effects of the milling speed on texture strength are the most significant, and the overall intensity of the texture decreases with the increase in the milling speed.
In Figure 14, the red line, green line, yellow line, and blue line represent the twin boundaries of { 10 1 ¯ 2 } , { 11 2 ¯ 1 } , { 11 2 ¯ 2 } , and { 10 1 ¯ 1 } , respectively. By observing the selected region, it is found that the content of all the twin grain boundaries is very low, and the volume fraction amounts to less than 1%, so that only two twin grain boundary distribution maps were selected to be illustrated. In all twin modes, the main activated twins are { 10 1 ¯ 1 } and { 10 1 ¯ 2 } , and the number of { 10 1 ¯ 1 } twins is much higher than that of { 10 1 ¯ 2 } , while a small number of { 11 2 ¯ 1 } twins can also be observed in Figure 14b. On the whole, when milling under the milling parameters selected in this study, the strain in the grain is coordinated mainly by the dislocation slip, and the contribution of the twins is very minimal.

4. Conclusions

1. The influence of the cutting depth on F t is the greatest, the influence of the linear velocity is the second greatest, and the influence of the feed rate is the least. For F n , the influence of the cutting depth is the greatest, the influence of the feed rate is the second greatest, and the influence of the linear velocity is the least. In the process of milling, the milling depth has the greatest influence on the overall grinding force.
2. Within the given range of the milling parameters, both the tangential milling force and the normal milling force show an overall increasing trend with the increase in the milling depth and feed rate. The tangential milling force decreases with the increase in the milling speed, and the normal milling force decreases first and then increases slightly with the increase in the milling speed.
3. The average error of the tangential milling force obtained from the prediction model in this study is 8.3%, and the average error of the normal milling force is 8.5%. This indicates that the prediction accuracy of the milling force prediction model considering the Taylor factor proposed in this study is better than the original model. The error may be caused by the tool wear and the change in the milling depth (caused by the horizontal degree of the machined surface), which are not considered in theoretical research.
4. The grain size decreases with the increase in the milling speed and milling depth, and the effect of the milling depth on grain size is the most significant. The proportion of the LAGB increases with the increase in the milling depth. Further analysis indicated that the influence of the milling depth on the LAGB distribution is the greatest. With the increase in the milling speed and milling depth, the proportion of the β phase increases, but with the increase in the feed rate, the proportion of the β phase decreases. Moreover, the milling depth has the most significant effect on the beta phase distribution. The texture intensity decreases with the increase in the milling speed and increases with the increase in the feed rate. Among all the milling parameters, the effect of the milling speed on the texture strength is the most significant. When milling under the milling parameters selected in this study, the strain in the grain is coordinated mainly by the dislocation slip, and the contribution of the twins is very minimal.

Author Contributions

Conceptualization, S.Y.L.; investigation, M.Z.; writing—original draft preparation, S.Z.; funding acquisition, J.M. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported and sponsored by the Shanghai Pujiang Program (20PJ1404700) and Young Scientific Research Team Cultivation Program of SUES (QNTD202112).

Data Availability Statement

The author declares that all data are available.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Li, S.; Han, J.; Yu, H.; Wang, J.; Lu, M.; Tian, Y.; Lin, J. Finite Element Investigation on Cutting Force and Residual Stress in 3D Elliptical Vibration Cutting Ti6Al4V. Micromachines 2022, 13, 1278. [Google Scholar] [CrossRef] [PubMed]
  2. Shen, L.; Chen, Y.; Zhu, H.; Lei, Y.; Qiu, C. Ti6Al4V Alloy Remelting by Modulation Laser: Deep Penetration, High Compactness and Metallurgical Bonding with Matrix. Micromachines 2022, 13, 1107. [Google Scholar] [CrossRef] [PubMed]
  3. Cui, Y.; Cai, J.; Li, Z.; Jiao, Z.; Hu, L.; Hu, J. Effect of Porosity on Dynamic Response of Additive Manufacturing Ti-6Al-4V Alloys. Micromachines 2022, 13, 408. [Google Scholar] [CrossRef] [PubMed]
  4. Wang, D.; Wang, H.; Chen, X.; Liu, Y.; Lu, D.; Liu, X.; Han, C. Densification, Tailored Microstructure, and Mechanical Properties of Selective Laser Melted Ti–6Al–4V Alloy via Annealing Heat Treatment. Micromachines 2022, 13, 331. [Google Scholar] [CrossRef]
  5. Ninpetch, P.; Kowitwarangkul, P.; Chalermkarnnon, P.; Promoppatum, P.; Chuchuay, P.; Rattanadecho, P. Numerical Modeling of Distortion of Ti-6Al-4V Components Manufactured Using Laser Powder Bed Fusion. Metals 2022, 12, 1484. [Google Scholar] [CrossRef]
  6. Liu, L.; Lu, M.; Wu, W.; Zhu, X. Experimental study on chip morphology of titanium alloy Ti-6Al-4V in High-speed milling. J. Mech. Eng. 2015, 196, 200–205. [Google Scholar]
  7. Merchant, M. Mechanics of the metal cutting process II. Plasticity conditions in orthogonal cutting. J. Appl. Phys. 1945, 16, 318–324. [Google Scholar] [CrossRef]
  8. Oxley, P. The Mechanics of Machining: An Analytical Approach to Assessing Machinability, 1st ed.; Ellis Horwood Ltd.: New York, NY, USA, 1989. [Google Scholar]
  9. Endres, W.; DeVor, R.; Kapoor, S. A Dual-Mechanism Approach to the Prediction of Machining Forces, Part 2: Calibration and Validation. J. Eng. Ind. 1995, 117, 534–541. [Google Scholar] [CrossRef]
  10. Moufki, A.; Devillez, A.; Dudzinski, D.; Molinari, A. Thermomechanical modelling of oblique cutting and experimental validation. Int. J. Mach. Tools Manuf. 2004, 44, 971–989. [Google Scholar] [CrossRef]
  11. Yang, Y. Study on Prediction Theory and Method of Milling Deformation of Titanium Alloy Aeronautical Integral Structure; Zhejiang University: Hangzhou, China, 2007. [Google Scholar]
  12. Das, L.; Nayak, R.; Saxena, K.K.; Nanda, J.; Jena, S.P.; Behera, A.; Sehgal, S.; Prakash, C.; Dixit, S.; Abdul-Zahra, D.S. Determination of Optimum Machining Parameters for Face Milling Process of Ti6A14V Metal Matrix Composite. Materials 2022, 15, 4765. [Google Scholar] [CrossRef]
  13. Kim, D.-H.; Lee, C.-M. Experimental Investigation on Machinability of Titanium Alloy by Laser-Assisted End Milling. Metals 2021, 11, 1552. [Google Scholar] [CrossRef]
  14. Perez, R.; Luis, N.; Gorka, U.; Octavio, P.; Silvia, M.; Jorge, B. On the relationship between cutting forces and anisotropy features in the milling of LPBF Inconel 718 for near net shape parts. Int. J. Mach. Tools Manuf. 2021, 170, 103801. [Google Scholar] [CrossRef]
  15. Ji, X. Prediction and Modeling of Residual Stress on Micro-Lubricated Cutting Surface; Shanghai Jiaotong University: Shanghai, China, 2014. [Google Scholar]
  16. Ren, Y. Study on Microstructure Evolution and Fine Structure of High Purity Titanium during Dynamic Plastic Deformation; Chongqing University: Chongqing, China, 2017. [Google Scholar]
  17. Tong, X.; Liu, X.; Yu, S. Anti-Friction and Anti-Wear Mechanisms of Micro Textures and Optimal Area Proportion in the End Milling of Ti6Al4V Alloy. J. Manuf. Mater. Process. 2019, 3, 91. [Google Scholar] [CrossRef]
  18. Su, Y.; Li, L. An Investigation of Cutting Performance and Action Mechanism in Ultrasonic Vibration-Assisted Milling of Ti6Al4V Using a PCD Tool. Micromachines 2021, 12, 1319. [Google Scholar] [CrossRef]
  19. Wu, Z.; Xing, Y.; Chen, J. Improving the Performance of Micro-Textured Cutting Tools in Dry Milling of Ti-6Al-4V Alloys. Micromachines 2021, 12, 945. [Google Scholar] [CrossRef] [PubMed]
  20. Dabwan, A.; Anwar, S.; Al-Samhan, A.M.; Nasr, M.M. On the Effect of Electron Beam Melted Ti6Al4V Part Orientations during Milling. Metals 2020, 10, 1172. [Google Scholar] [CrossRef]
  21. Von, M. Mechanics of the ductile form changes of crystals. Z. Fur Angew. Math. Und Mech. 1928, 8, 161–185. [Google Scholar]
  22. Taylor, G.I. Plastic strain in metals. J. Inst. Met. 1938, 62, 307–324. [Google Scholar]
  23. Mao, J.; Zhao, M.; Zhang, L. Mechanism and experimental study of grain orientation on grinding force in micromachining. J. Mech. Eng. 2021, 757, 262–272. [Google Scholar]
  24. Zhao, M. Study on the Mechanism and Prediction Modeling of Grain Orientation on Micro-Grinding Surface Residual Stress; Donghua University: Shanghai, China, 2019. [Google Scholar]
  25. Groves, G. Independent slip systems in crystals. Philos. Mag. 1963, 8, 877–887. [Google Scholar] [CrossRef]
  26. Rosi, F.; Dube, C.; Alexander, B. Mechanism of plastic flow in titanium-determination of slip and twinning elements. Trans. Aime. 1953, 197, 257. [Google Scholar]
  27. Rosi, F.; Perkins, F.; Seigle, L. Mechanism of plastic flow in titanium at low and high temperatures. JOM 1956, 8, 115–122. [Google Scholar] [CrossRef]
  28. Zhong, Y.; Yin, F.; Nagai, K. Role of deformation twin on texture evolution in cold-rolled commercial-purity Ti. J. Mater. Res. 2011, 23, 2954–2966. [Google Scholar] [CrossRef]
  29. Xiao, L. Twinning behavior in the Ti–5at% Al single crystals during cyclic loading along. Mater. Sci. Eng. A 2005, 394, 168–175. [Google Scholar] [CrossRef]
  30. Yi, N.; Hama, T.; Kobuki, A. Anisotropic deformation behavior under various strain paths in commercially pure titanium Grade 1 and Grade 2 sheets. Mater. Sci. Eng. A 2016, 655, 70–85. [Google Scholar] [CrossRef] [Green Version]
  31. Yoo, M. Slip, twinning, and fracture in hexagonal close-packed metals. Metall. Trans. A 1981, 12, 409–418. [Google Scholar] [CrossRef]
  32. Jackson, A. Slip Systems; Springer: New York, NY, USA, 1991; pp. 83–88. [Google Scholar]
  33. Meng, L.; Yang, P.; Cui, F.; Bian, J.; Xie, Q. Orientation imaging analysis of deformation twins in magnesium alloys. In Proceedings of the 11th China Symposium on Stereology and Image Analysis, Ningbo, China, 11 May 2006; pp. 293–297. [Google Scholar]
  34. Waldorf, D.; DeVor, R.; Kapoor, S. A slip-line field for ploughing during orthogonal cutting. J. Manuf. Sci. Eng. 1998, 120, 693–698. [Google Scholar] [CrossRef]
  35. Hall, E. The Deformation and Ageing of Mild Steel: III Discussion n of Results. Proc. Phys. Society. Sect. B 1951, 64, 747–753. [Google Scholar] [CrossRef]
  36. Hansen, N.; Huang, X.; Winther, G. Grain orientation, deformation microstructure and flow stress. Mater. Sci. Eng. A 2008, 494, 61–67. [Google Scholar] [CrossRef]
  37. Hughes, G.; Smith, S.; Pande, C. Hall-Petch strengtheninig for the microhardness of twelve nanometer grain diameter electrodeposited nickel. Scr. Metall. 1986, 20, 93–97. [Google Scholar] [CrossRef]
  38. Hughes, D.; Hansen, N. Microstructure and strength of nickel at large strains. Acta Mater. 2000, 48, 2985–3004. [Google Scholar] [CrossRef]
  39. Wang, X. Simulation Study on Microstructure Evolution of Hot Deformation Dynamic Recrystallization of Ti-6Al-4V-DT Titanium Alloy; Nanjing University of Aeronautics and Astronautics: Nanjing, China, 2012. [Google Scholar]
  40. Klocke, F.; Hoppe, S.; Schulz, H. Mechanisms of Chip Formation in High Speed Cutting; Scientific Fundamentals of HSC, Hanser Press: Butler, OH, USA, 2001; p. 1. [Google Scholar]
  41. Luo, H.; Zhao, J.; Li, A.; Cui, X. Experimental study on milling force of high speed milling titanium alloy ti6a4v. Modul. Mach. Tool Autom. Mach. Technol. 2011, 5, 18–20. [Google Scholar]
  42. Yi, J. High Speed Milling Force Test and Finite Element Numerical Analysis of Titanium Alloy; Nanjing University of Aeronautics and Astronautics: Nanjing, China, 2009. [Google Scholar]
  43. Li, K.; Liang, S. Modeling of cutting forces in near dry machining under tool wear effect. Int. J. Mach. Tools Manuf. 2007, 47, 1292–1301. [Google Scholar] [CrossRef]
Figure 1. The samples and the experimental equipment used in the milling experiment.
Figure 1. The samples and the experimental equipment used in the milling experiment.
Micromachines 13 01618 g001aMicromachines 13 01618 g001b
Figure 2. Schematic representation of the EBSD of single crystal planes.
Figure 2. Schematic representation of the EBSD of single crystal planes.
Micromachines 13 01618 g002
Figure 3. Schematic diagram of the crystal structures of α -Ti and β -Ti. (a) α -Ti. (b) β -Ti.
Figure 3. Schematic diagram of the crystal structures of α -Ti and β -Ti. (a) α -Ti. (b) β -Ti.
Micromachines 13 01618 g003
Figure 4. Geometry of the milling force model. (a) Chip-forming force model [15]. (b) Ploughing force model [34].
Figure 4. Geometry of the milling force model. (a) Chip-forming force model [15]. (b) Ploughing force model [34].
Micromachines 13 01618 g004
Figure 5. Microstructure maps of the Ti-6Al-4V surface before milling. (a) The grain orientation distribution map. (b) The grain size statistic map. (c) The phase distribution map. (d) The KAM map.
Figure 5. Microstructure maps of the Ti-6Al-4V surface before milling. (a) The grain orientation distribution map. (b) The grain size statistic map. (c) The phase distribution map. (d) The KAM map.
Micromachines 13 01618 g005
Figure 6. Low frequency filtering of the milling force signal. (a) Before FFT. (b) After FFT.
Figure 6. Low frequency filtering of the milling force signal. (a) Before FFT. (b) After FFT.
Micromachines 13 01618 g006
Figure 7. Effect of the level of the milling parameters on the milling force. (a) Tangential milling force F t . (b) Normal milling force F n .
Figure 7. Effect of the level of the milling parameters on the milling force. (a) Tangential milling force F t . (b) Normal milling force F n .
Micromachines 13 01618 g007
Figure 8. Comparison between the experimental results and predicted values. (a) Tangential milling force F t . (b) Normal milling force F n .
Figure 8. Comparison between the experimental results and predicted values. (a) Tangential milling force F t . (b) Normal milling force F n .
Micromachines 13 01618 g008
Figure 9. The grain orientation distribution maps of Ti-6Al-4V under different milling parameters.
Figure 9. The grain orientation distribution maps of Ti-6Al-4V under different milling parameters.
Micromachines 13 01618 g009
Figure 10. The grain size statistics of Ti-6Al-4V under different milling parameters.
Figure 10. The grain size statistics of Ti-6Al-4V under different milling parameters.
Micromachines 13 01618 g010
Figure 11. The grain boundary distribution maps of Ti-6Al-4V under different milling parameters.
Figure 11. The grain boundary distribution maps of Ti-6Al-4V under different milling parameters.
Micromachines 13 01618 g011
Figure 12. The grain boundary distribution maps of Ti-6Al-4V under different milling parameters.
Figure 12. The grain boundary distribution maps of Ti-6Al-4V under different milling parameters.
Micromachines 13 01618 g012
Figure 13. Pole figures of the Ti-6Al-4V surface under different milling parameters.
Figure 13. Pole figures of the Ti-6Al-4V surface under different milling parameters.
Micromachines 13 01618 g013aMicromachines 13 01618 g013b
Figure 14. Pole figures of the Ti-6Al-4V surface under different milling parameters. (a) Sample 1. (b) Sample 9.
Figure 14. Pole figures of the Ti-6Al-4V surface under different milling parameters. (a) Sample 1. (b) Sample 9.
Micromachines 13 01618 g014
Table 1. The main chemical composition of Ti-6Al-4V.
Table 1. The main chemical composition of Ti-6Al-4V.
AlVFeSiCNHOTi
wt%6.133.930.1030.110.0070.0350.0130.12Bal.
Table 2. Orthogonal experimental parameters.
Table 2. Orthogonal experimental parameters.
Factors V c   ( m / m i n ) f   ( m m / m i n ) a p   ( m m )
Levels
1601600.1
2902000.2
31202400.3
Table 3. Main slip systems and corresponding parameters of α-Ti [25].
Table 3. Main slip systems and corresponding parameters of α-Ti [25].
DirectionTypesCrystallographic ElementsIndependent Slip System
a Basal slip ( 0002 ) 11 2 ¯ 0 2
a Prismatic slip { 10 1 ¯ 0 } 11 2 ¯ 0 2
a Pyramidal slip { 10 1 ¯ 1 } 11 2 ¯ 0 4
a + c Pyramidal slip { hkil } 11 2 ¯ 3 ¯ 5
Table 4. Main twinning systems and corresponding parameters of α-Ti [28].
Table 4. Main twinning systems and corresponding parameters of α-Ti [28].
TypesTwin ModesRotation Angle ( ° ) Rotation Axis
Extension twin { 10 1 ¯ 2 } 10 1 ¯ 1 ¯ 84.78 11 2 ¯ 0
{ 11 2 ¯ 1 ¯ } 11 2 ¯ 6 ¯ 35.10 10 1 ¯ 0
Contraction twin { 10 1 ¯ 1 } 10 2 ¯ 0 57.22 11 2 ¯ 0
{ 11 2 ¯ 2 } 11 2 ¯ 3 ¯ 64.62 10 1 ¯ 0
{ 11 2 ¯ 4 } 22 4 ¯ 3 ¯ 76.66 10 1 ¯ 0
Table 5. Main slip systems and corresponding parameters of β -Ti [32].
Table 5. Main slip systems and corresponding parameters of β -Ti [32].
Slip PlaneSlip Direction
( 101 ) [ 1 ¯ 1 ¯ 1 ] [ 1 ¯ 11 ]
( 10 1 ¯ ) [ 1 1 ¯ 1 ] [ 111 ]
( 011 ) [ 1 ¯ 1 ¯ 1 ] [ 1 1 ¯ 1 ]
( 0 1 ¯ 1 ) [ 1 ¯ 11 ] [ 111 ]
( 1 1 ¯ 0 ) [ 1 ¯ 1 ¯ 1 ] [ 111 ]
( 110 ) [ 1 1 ¯ 1 ] [ 1 ¯ 11 ]
Table 6. The physical parameters of Ti-6Al-4V.
Table 6. The physical parameters of Ti-6Al-4V.
Materials ρ   ( k g / m 3 ) E   ( G p a ) υ λ   ( W / m · K ) C   ( J / k g · K )
Ti-6Al-4V4430 109   ( 50   ) 0.34 6.8   ( 20   ) 611   ( 20   )
Table 7. The Johnson–Cook parameters of Ti-6Al-4V [39].
Table 7. The Johnson–Cook parameters of Ti-6Al-4V [39].
Materials A   ( M p a ) B   ( M p a ) C n m T m   ( ) T w   ( ) α b(m) G   ( N / m 2 )
Ti-6Al-4V4183940.0350.4711650250.8 2.86 × 10 10 2.05 × 10 10
Table 8. The material parameters of the tool.
Table 8. The material parameters of the tool.
Materials ρ   ( k g / m 3 ) E   ( M p a ) υ λ   ( W / m · K ) C   ( J / k g · K )
YG6X14,600640,0000.2279.6176
Table 9. Experimental results of Ti-6Al-4V under dry milling.
Table 9. Experimental results of Ti-6Al-4V under dry milling.
No. V c
( m / m i n )
f
( m m / m i n )
a p
( m m )
F t
( N )
F n
( N )
P r e d   F t
( N )
Pred   F n
( N )
1601600.125.313.426.412.9
2602000.244.123.042.624.3
3602400.362.832.959.734.4
4901600.229.917.428.718.9
5902000.343.625.839.827.5
6902400.123.216.819.914.4
71201600.325.622.828.024.6
81202000.116.917.913.415.7
91202400.224.020.122.717.5
Table 10. Range analysis of the tangential milling force F t .
Table 10. Range analysis of the tangential milling force F t .
ItemLevel V c   ( m / m i n ) f   ( m m / m i n ) a p   ( mm )
K1132.280.865.4
296.7104.698
366.5110132
K-avg144.0726.9321.80
232.2334.8732.67
322.1736.6744.00
R21.909.7322.20
Table 11. Range analysis of the normal milling force F n .
Table 11. Range analysis of the normal milling force F n .
ItemLevel V c   ( m / m i n ) f   ( m m / m i n ) a p ( mm )
K169.353.648.1
260.066.760.5
360.869.881.5
K-avg123.1017.8716.03
220.0022.2320.17
320.2723.2727.17
R3.105.4011.13
Table 12. The relative errors of the milling forces and predictions for each group.
Table 12. The relative errors of the milling forces and predictions for each group.
No. F t
( N )
P r e d   F t
( N )
Relative Error (%) F n
( N )
Pred   F n
( N )
Relative Error (%)
125.326.44.313.412.93.7
244.142.63.423.024.35.7
362.859.74.932.934.44.6
429.928.74.017.418.98.6
543.639.88.725.827.56.6
623.219.914.216.814.414.3
725.628.09.422.824.67.9
816.913.420.717.915.712.3
924.022.75.420.117.512.9
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhu, S.; Zhao, M.; Mao, J.; Liang, S.Y. A Ti-6Al-4V Milling Force Prediction Model Based on the Taylor Factor Model and Microstructure Evolution of the Milling Surface. Micromachines 2022, 13, 1618. https://doi.org/10.3390/mi13101618

AMA Style

Zhu S, Zhao M, Mao J, Liang SY. A Ti-6Al-4V Milling Force Prediction Model Based on the Taylor Factor Model and Microstructure Evolution of the Milling Surface. Micromachines. 2022; 13(10):1618. https://doi.org/10.3390/mi13101618

Chicago/Turabian Style

Zhu, Siyuan, Man Zhao, Jian Mao, and Steven Y. Liang. 2022. "A Ti-6Al-4V Milling Force Prediction Model Based on the Taylor Factor Model and Microstructure Evolution of the Milling Surface" Micromachines 13, no. 10: 1618. https://doi.org/10.3390/mi13101618

APA Style

Zhu, S., Zhao, M., Mao, J., & Liang, S. Y. (2022). A Ti-6Al-4V Milling Force Prediction Model Based on the Taylor Factor Model and Microstructure Evolution of the Milling Surface. Micromachines, 13(10), 1618. https://doi.org/10.3390/mi13101618

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop