Experimental study on the effect of the milling condition of an aluminum alloy on subsurface residual stress

Aluminum alloys used in monolithic parts for aerospace applications are subjected to distortion and residual stress (RS) generated by milling, affecting the product fatigue life. Particularly, the change in RS with depth (z) has a characteristic distribution with a maximum compressive RS at a z several tens of micrometers from the surface; however, the RS value depends on the measurement method used. In this study, the RS distribution with z from the surface after milling was measured for the AA7050-T7451 aluminum alloy by two-dimensional X-ray diffraction (2D method). The results were compared with those of four prior measurement methods, and the validity of 2D method was verified. The changes in subsurface RS with z showed similar distributions under all measurement conditions except when cos(α)-XRD was employed. The 2D method provides high repeatability. The in-plane RS distribution was also measured using 2D method to investigate the effect of milling conditions on this distribution. The RS values varied markedly depending on the measurement position, particularly at a small collimator diameter of 0.146 mm, allowing detection of localized extreme RS values. The maximum RS at z = 0 mm was − 85.6 MPa at a cutting speed of vc = 200 m/s and feed per tooth of fz = 0.05 mm, while it was − 16 MPa for vc = 450 m/s and 6.8 MPa for fz = 0.2 mm, revealing that the compressive RS changes to tensile RS as vc and fz increase.


Introduction
Monolithic structural components are widely used in aircraft manufacturing and require light products, short assembly times, high production efficiencies, and low production costs [1]. Aluminum alloys are used as monolithic components for aerospace applications owing to their low cost, high specific strength, low density, good thermal conductivity, corrosion resistance, and other material properties that facilitate the formation of complex shapes [2,3]. Such monolithic parts are fabricated using machining techniques, such as milling. Milling is a cutting process that removes workpiece material by pressing a rotating blade against a fixed workpiece, thereby achieving high machining efficiency because a wide area can be removed simultaneously. Moreover, milling can be used for almost all materials [4], and high surface quality can be obtained, even for high-strength, difficult-to-machine materials [5]. However, the thinness and large size of aircraft aluminum alloy parts reduce their rigidity, and high material removal of as much as 90% causes residual stress (RS) on the machined surface, which easily results in distortion [1,6]. The distortion and RS caused by milling can also reduce 1 3 fatigue life [7]. Aircraft manufacturers are estimated to incur hundreds of millions of dollars or more in economic losses owing to scrap materials and rework related to such component distortions [8]. Therefore, to improve aircraft safety, optimize manufacturing, and reduce costs, it is necessary to clarify the causes of RS.
Pioneering studies have attempted to derive the RS distribution with depth (z) in the milled surface region using the electrolytic etching-deflection technique [9]. In recent years, researchers have experimentally confirmed that the maximum compressive RS arises at a z near the surface and that RS gradually decreases with increasing z [10][11][12]. Such a z-directed RS distribution is called a "root" (√) shape [1]. Feng et al. used numerical analysis to obtain the RS distribution with z of milling and showed that the distribution is √-shaped [13,14], similar to prior experimental results [1,9,10,12]. The z range of the compressive RS and the z and magnitude of the maximum compressive RS in the near-surface RS distribution vary depending on milling parameters, such as cutting speed (v c ), feed per tooth (f z ), and z of the cut [15][16][17]. Tensile RS in metallic materials accelerates fatigue crack propagation, whereas compressive RS inhibits crack propagation [18]. For example, in trochoidal milling [19], the cutting tool is moved in a spiral motion to maintain the cutting rate and reduce the cutting resistance and RS. However, because compressive RS plays a role in suppressing fatigue failure, it is not necessarily preferable to reduce the RS. Therefore, clarifying the milling parameters that affect the degree of tensile/compressive RS in the milled material can improve the fatigue life and reduce the product cost. However, when the RS is measured experimentally, it can vary significantly depending on the measurement method, even under the same milling conditions. Several experimental methods exist for measuring RS [20]. Hole-drilling and slotting methods are commonly used as mechanical measurement methods, whereas X-ray diffraction (XRD) methods are X-ray-based and can be further divided into sin 2 (ψ), cos(α), and two-dimensional (2D) methods, depending on the measurement technique. For example, XRD has recently been used to measure the RS in austenitic stainless steels treated at cryogenic temperatures [21]. Vrabeľ et al. measured the RS distribution with z of Inconel 718 after using the sin 2 (ψ) and hole-drilling methods [22]. However, they used the sin 2 (ψ) method for the top surface and the hole-drilling method to measure the RS distribution with z but did not compare measurement methods [22]. Bordinassi et al. compared the surface RS after milling measured using the blind hold method (a mechanical measurement method) with that measured using the cos(α) method; however, they did not compare the z RS distributions [23]. Chighizola et al. mentioned the lack of a systematic comparison of these measurement methods and conducted z RS distribution measurements of milled aluminum using the hole-drilling, slotting, sin 2 (ψ) and cos(α) methods [24]. Consequently, although hole-drilling requires measurement skills, the results with the highest repeatability were obtained at a z of 0.03 mm or greater, whereas the results obtained using the cos(α) method disagreed with the other measurement results. In contrast, the 2D method, which, like the cos(α) method, has a 2D detector, has been shown to have higher measurement accuracy than the sin 2 (ψ) method, which is mostly used for RS measurements in laboratories [25]. RS can be measured in submillimeter areas using a collimator [25]. Therefore, the validity of the 2D method can be investigated by comparing the measurement results with those of the hole-drilling, slotting, sin 2 (ψ) and cos(α) methods [24]. In addition, if large differences resulting from milling exist in the RS distribution in the plane direction, they will distort the machined parts. However, few studies have investigated the distribution of RS on machined surfaces in detail.
In this study, RS measurements using the 2D method were performed on the same samples as those in [24], an aluminum alloy (AA7050-T7451), after milling, to investigate the validity of the 2D method. Furthermore, the in-plane RS distribution, which has not been investigated experimentally or analytically, was evaluated using the 2D method. First, the RS distribution with z was investigated and compared with the major RS measurement methods (hole-drilling, slotting, sin 2 (ψ), and cos(α) methods), and the measurement accuracy was evaluated. Subsequently, the 2D method was used to precisely measure the distribution of RS with z in the inplane direction, and the effect of the milling conditions on the subsurface RS was investigated.

Workpiece material and milling conditions
A common aerospace aluminum sheet (AA7050-T7451) with a thickness of 25 mm was used as the experimental material, and the samples were prepared as previously described [24]. The AA7050-T7451 standard is defined in AMS4050, and AA indicates that this material was provided by Aluminum Association, Inc. Aluminum 7050 has high strength coupled with high resistance to exfoliation corrosion and stress corrosion cracking, high fracture toughness, and fatigue resistance. T7451 indicates that this material was solution-annealed and the RS was removed and further overaged. A milling machine (DMU 70 CNC, DMG Mori Co., Ltd.) and a cemented carbide end mill (F3AA1200AWL, Kennametal, Inc.) were used. Figure 1a shows a schematic of the milling parameters, Fig. 1b shows a photograph of the experimental setup [26], and Fig. 1c shows the layout of the measurement workpiece.
As shown in Fig. 1a and b, the x-, y-, and z-axes are defined as the length, width, and thickness directions of the sample, and the workpiece sizes are 206, 106, and 28 mm. The tool was a square with a diameter D = 12 mm and three cutting edges with a torsional angle of 45° (Fig. 1a). The z of cut a p and width of cut a e were set to 3 and 4 mm, respectively. The workpiece was fixed using a vice, and a piezoelectric dynamometer was used to record the milling forces (Fig. 1b). As shown in Fig. 1c, the milled workpieces were subdivided into area grids of x = 25.4 mm and y = 34.0 mm. RS was measured near the center of the grids. Three measurement areas were randomly selected, and the repeatability of each RS measurement method was investigated. Four conditions were set for v c and f z (Fig. 2).
The effect of f z on RS was examined by changing f z to 0.04, 0.10, and 0.20 mm with v c = 200 m/min for Modes 1-3, and the effect of v c on RS was also examined by changing v c to 200 and 400 m/min with f z = 0.04 mm for Modes 1 and 4. The specific milling procedure was as follows. Beginning at the milling origin, the tool was cut to a p in the positive direction of the z-axis and a p in the positive direction of the y-axis and then scanned 200 mm toward the negative direction of the x-axis. This operation was defined as a single path. After completing one path, the tool was returned to the machining origin along the x-axis and offset by a e along the y-axis; the same path was repeated. Milling was completed when a total of 100 mm was cut in the y-direction. Material from the same lot was used for each milling process. As the milling area is sufficiently large for v c and f z , the repeatability of each measurement method can be evaluated by markedly changing the measurement area, as shown in Fig. 1c. The machined surface was observed via scanning electron microscopy (SEM; JCM-7000, JEOL, Ltd.), and the surface profile was measured using a contact-type surface roughness meter.  1 3

RS measurement methods
RS was measured using 2D method and compared with results of the hole-drilling, slotting, sin 2 (ψ), and cos(α) methods [24]. These measurement methods are briefly described below: (i) The hole-drilling method, as specified in ASTM E837-20, calculates the RS from the three-directional strain generated by the stress release. In this approach, holes are drilled incrementally at the centers of the rosette strain gauges placed on the specimen [27].
(ii) The slotting method is similar to the one-dimensional (1D) hole-drilling method. A single component of RS can be measured from the strain generated in a uniaxial strain gauge by cutting a long, narrow groove perpendicular to the direction of the RS to be measured [28,29].
(iii) The sin 2 (ψ) method is a basic RS measurement technique that uses XRD. This is performed using a zerodimensional or 1D detector [30,31]. First, the lattice plane spacing d is derived from the diffraction angle 2θ when X-rays are incident at various angles ψ. The residual strain calculated from d is plotted against sin 2 ψ to obtain the slope of the regression line. RS can be calculated by multiplying the slope by the elastic constant (1/2 S 2 ) of the lattice plane.
(iv) The cos(α) method employs a 2D detector to detect a single incident X-ray beam and determine the RS from the obtained Debye-Scherrer rings [32]. Based on an unstressed sample, the strain with respect to the center angle α of the Debye-Scherrer ring can be determined, and RS can be calculated from the slope of the strain with respect to cos(α).
The hole-drilling and slotting methods have the advantages of being fast, easy to use, and portable. However, they have the disadvantages of data interpretation and limited strain sensitivity and resolution [33]. The sin 2 (ψ) and cos(α) methods have the advantages of the ability to measure a wide range of materials and the possibility of measuring macro-and micro-RSs; their disadvantages include the use of laboratory-based systems and small components and only providing basic measurements [33]. In each method, the material is incrementally removed from the surface, and measurements are taken at each step to obtain a RS profile from the surface to z.
The 2D method involves using a 2D detector. However, the stress components are obtained by applying X-rays from various directions using the least squares method [34]. The left image in Fig. 3 shows the XRD system configuration, and the right image shows the coordinate system of the measured sample.
The sample was irradiated with X-rays from various incident directions by adjusting φ, ψ, and ω, as shown in Fig. 3, and some of the Debye rings, which were diffraction patterns generated from the sample, were detected using a 2D positionsensitive proportional counter (2D-PSPC). RS was obtained by measuring the diffraction angle 2θ at each position in the χ-direction (χ = 90°ψ). If the stress in the j-direction relative to the i-plane is σ ij , then all stress components in the sample coordinate system are expressed as in Eq. (1), and each stress component can be measured using the least squares method [34]. In Eq. (1), f ij is a function of ψ, φ, ω, χ, and θ and can be obtained from Eqs. (2)-(8) [34].
(1)  As shown in Fig. 3 (right), the y-axis origin was set as the edge of each milling path, and the x-and y-axes were set in the length and width directions of the sample, respectively.

Conditions for measuring the RS distribution with z
An X-ray diffractometer (D8 Discover, Bruker Corp.) and 2D-PSPC (V500, Bruker Corp.) were used for the RS measurements, and the RSs were calculated using LEP-TOS software (version 7.9, Bruker Japan K. K.). The X-ray source was Cr-Kα radiation (wavelength 2.29093 Å, 35 kV, 40 mA). The samples were irradiated with X-rays after narrowing the beam diameter using a collimator diameter (d c ) to 0.8 mm. The exposure time per frame was set to 20 s. Table 1 lists the measurement conditions for the 2D method. Eight conditions were applied in the φ-direction (0-315°) and three in the ψ-direction (0°, 30°, and 60°), and 24 frames were measured while swinging in the ω-direction from 107° to 115°. During the measurement, the sample stage was scanned within ± 2 mm in the x-and y-directions, and diffraction information was obtained for each 4 mm × 4 mm area. Table 2 lists the analyzed areas of the Debye ring.
The RS was calculated using the diffraction angles from the Al [311] plane by analyzing the Debye rings in the range of 135-145° in the 2θ-direction and 65-115° in the χ-direction. The change in the RS distribution with z was obtained by gradually removing the material surface by electropolishing and measuring the stress in the exposed plane. Phosphoric acid with water and aluminum was used as the electropolishing fluid, and the removal z was measured after each electropolishing using a micrometer. The etch pit was rectangular with an in-plane area of x = 3 mm × y = 8 mm. The RS versus z profile was measured at three different locations along x for y = 0, as was done in a previous study [24], and the average of the three z profiles and the standard deviations at each z were compared to the prior results [24]. No z correction was performed, and raw XRD data were used.

Conditions for measuring the RS distribution in plane
Two collimators with d c = 0.8 and 0.146 mm were used to measure the RS distribution in the plane direction of the sample. For the d c = 0.8-and 0.146-mm collimators, the exposure times per frame were set to 40 s and 20 min, respectively. The measurement conditions were the same as those listed in Table 1, and 24 frames were used. During the measurements, the sample stage was not scanned; however, a fixed-point area was measured. These measurements were repeated at different locations in the range of − 2 mm ≤ y ≤ 2 mm, and the RS distribution in the y-axis direction in the xy-plane was measured at 9 and 17 points at equal intervals for the collimators with d c = 0.8 and 0.146 mm, respectively. The analysis range of the Debye ring was set as listed in Table 2, and the RSs were calculated using the procedure described in Sect. 2.3. This measurement was performed for each z within an area of x = 3 mm × y = 8 mm, achieved by electropolishing as described in Sect. 2.3. The maximum height roughness (R z ) of Modes 1-4 were 3.7, 9.9, 25.1, and 7.3 μm, and the arithmetical mean heights (R a ) were 0.7, 1.8, 3.1, 1.1 μm, respectively. In this study, a e was set to 4 mm, and each profile indicated that the surface geometry changed every 4 mm. Figure 4a-c, corresponding to Modes 1-3 with v c = 200 m/min, show that the R z and R a increased as f z increased, and the R z and R a values of Mode 3 (Fig. 4c) were approximately 6.8 and 4.4 times higher, respectively, than those in Mode 1 (Fig. 4a). The f z value of Mode 3 is five times greater than that of Mode 1. R z and R a were used to evaluate the projection and average surface characteristics, respectively. Both roughness values increased significantly with increasing f z . Next, by comparing Mode 1 in Fig. 4a and Mode 4 in Fig. 4d with different values of v c , the v c in Mode 4 was found to be 2.25 times greater than that in Mode 1, and R z and R a were 2.0 and 1.6 times greater, respectively. Considering the increased rate of surface roughness with respect to each parameter value, f z had a greater influence on the surface roughness than v c . Although milling conditions with small f z and v c are desirable for high surface accuracy, the machining efficiency can be improved by increasing f z and v c . Furthermore, the milling

Evaluation of measurement accuracy of the 2D method
The RS of the milled samples with z under the same conditions was compared with the results of the hole-drilling, slotting, sin 2 (ψ), and cos(α) methods [24] to assess the measurement accuracy of the RS from the 2D method. However, the measurement positions differed slightly with z for each sample. Therefore, the constant z RS values were interpolated using the respective measurement results, and the interpolated values were compared. Figure 5 shows the RS distribution with z, as measured using each measurement method. The results show σ yy stress components milled in Mode 3. The three different RS measurements are indicated by the open symbols (triangles, squares, and diamonds). As shown in Fig. 1c, three different points on one sample were randomly selected for RS measurements with z. The average RS obtained by interpolation from each of the three measurement results is indicated by closed-color symbols, and the error bars indicate the standard deviations.
In Fig. 5a-e, the results, other than those of the cos(α) method, show a √-shaped RS distribution. During milling, the tool and workpiece rub against each other, introducing a large RS near the milled surface. Because the material was heat-treated to remove RS before milling, the RS approached zero away from the milled surface. Therefore, the RS distribution with z became the √ shape. The maximum RSs values at the interpolated z obtained using each measurement method are shown in Fig. 5f. The maximum compressive RS was at z = 0.055 mm for the slotting and 2D methods, z = 0.041 mm for the hole-drilling method, z = 0.028 mm for the sin 2 (ψ) method, and z = 0.015 mm for the cos(α) method. Focusing on repeatability, the standard deviations of the mechanical measurement methods-the hole-drilling and slotting methods ( Fig. 5a and b) were larger near the surface but became smaller as the measurement position became deeper. In particular, the hole-drilling method resulted in almost no measurement errors at z > 0.03 mm. Specifically, the standard deviations of the sin 2 (ψ) and 2D methods at z = 0.005 mm were 1% and 8% of those of the hole-drilling method, respectively. In contrast, for the XRD methods (the sin 2 (ψ), cos(α), and 2D methods; see Fig. 5c-e, respectively) the cos(α) method resulted in a large standard deviation at all measured z values, whereas the sin 2 (ψ) and 2D methods resulted in small standard deviations at all z values.
In Fig. 6, the interpolated values obtained from Fig. 5 are compared. The standard deviations are shown as error bars. Figure 6a shows the RS distribution obtained from each measurement method with z; all the interpolated values  In a previous study [24], the highest repeatability was obtained using the hole-drilling method compared with those using the slotting, sin 2 (ψ), and cos(α) methods. Therefore, repeatability was evaluated by comparing the respective measurement results with those of the hole-drilling method. Specifically, Fig. 6b-e show the measurement results of the slotting, sin 2 (ψ), cos(α), and 2D methods, respectively, compared with those of the hole-drilling method. However, the standard deviation from using the hole-drilling method near the surface was large; therefore, the values of z ≥ 0.05334 mm were compared. The correlation coefficients for the results shown in Fig. 6b-e are 0.982, 0.895, 0.657, and 0.973, respectively, and the uncorrelated probabilities are 0.0003, 0.11, 10.9, and 0.0011%, respectively. Thus, the 2D method results are highly positively correlated with the results of the holedrilling method. By contrast, the slopes of the regression lines in Fig. 6b-e are 1.2179, 0.9709, 0.2867, and 0.7965, respectively. Hence, the sin 2 (ψ) method was the most consistent with the absolute RS value of the hole-drilling method. However, the sin 2 (ψ) method provides only partial diffraction information from Debye rings and determines the RS values by linear regression between two parameters obtained from a certain  [35]. If the target material and determination method of the 2θ value are not optimized, the measurement error is extremely large [35], and the 2D method has a smaller repetitive measurement error [25]. Therefore, the sin 2 (ψ) method is effective if the material and 2θ value conditions are satisfied to increase the degree of regression linearity for determining RS. The 2D method obtains diffraction information from a wide range of Debye rings by using a 2D detector to determine the RS and obtain results with high repeatability. Although there is a slight difference (approximately 50 MPa at the maximum compressive RS) from the absolute value measured using the hole-drilling method, it is possible to perform highly accurate measurements with a small error, even in the near-surface area where the mechanical measurement is unsatisfactory. Repeatability is important when evaluating RS distribution trends. The 2D method, in which the measurement area can be selected only by changing the collimator diameter, is extremely effective for measuring the RS in a local area.  Modes 1-4. The milling tool was scanned in the x-axis positive direction, and σ xx , σ yy , and σ xy were the RS components in the x-, y-, and shear directions, respectively. The standard deviations are shown as error bars.

Effects of milling parameters on the RS distribution with z using the 2D method
To investigate the effect of f z , the RS distributions of Modes 1-3 were compared, as shown in Fig. 7a-c, in which v c was constant at 200 m/min and only f z was changed. The z ranges of RS in all RS components were z ≤ 0.06 mm for Mode 1 with f z = 0.04 mm, as shown in Fig. 7a, z ≤ 0.08 mm for Mode 2 with f z = 0.1 mm, as shown in Fig. 7b, and z ≤ 0.13 mm for Mode 3 with f z = 0.2 mm, as shown in Fig. 7c, indicating that the larger the f z , the deeper the RS. The compressive RS peaks in the √ shape for σ xx and σ yy occurred at greater z values as f z increased, with peaks for σ xx in Modes 1-3 at − 112.1 ± 6.7 MPa at z = 0.01 mm, − 124.9 ± 7.2 MPa at z = 0.03 mm, and − 139.9 ± 7.7 MPa at z = 0.05 mm, respectively, and peaks for σ yy at − 152.4 ± 6.7 MPa at z = 0.01 mm, − 164.0 ± 7.4 MPa at z = 0.01 mm, and − 122.4 ± 7.7 MPa at z = 0.05 mm, respectively. In other words, in the σ xx -direction, the z and magnitude of the maximum compressive RS monotonically increased, while, in the σ yy -direction, the z of the maximum compressive RS increased for f z = 0.2 mm, and the magnitude increased to f z = 0.1 mm and then decreased. Tang et al. performed milling on 7050-T7451 and showed that tangential and axial forces increase with increasing f z and that the maximum RS occurs at a greater z for both σ xx and σ yy RS components from XRD methods [36]. Denkena et al. measured the RS of the σ xx component after milling using the sin 2 (ψ) method and found that as f z increased, the tangential and axial forces increased, and the maximum compressive RSs became increasingly deeper and larger [15]. These results are similar to those of the present study. The RSs values in the σ yy -direction were compared using various measurement methods [24]. Although the RSs occurred at greater z values with increasing f z , the z and magnitudes of the RS peaks differed for each measurement method owing to large measurement errors. However, as shown in Fig. 7, the 2D method has a smaller measurement error, enabling the relative evaluation of the f z -dependent compressive RS peaks for the σ xx and σ yy components. This indicates that the 2D method can be used to investigate the RS distribution in greater detail than conventional measurement methods. Comparing the maximum values of σ xx and σ yy measured by the 2D method reveals that the maximum compressive RS in Fig. 7a and b with f z ≤ 0.1 mm was larger for σ yy than for σ xx . As shown in Fig. 7c, with f z = 0.2 mm, σ xx was slightly larger than σ yy , but σ xx and σ yy showed similar RS distributions in terms of peak position and magnitude.
In milling, machining is performed while the tool rotates; therefore, if the tool is scanned in the x-direction, machining is preferentially performed in the y-direction, which is the direction of the rotational circumference of the tool. Furthermore, if v c and the number of flutes Z of the tool remain the same, f z is determined by the number of tool revolutions n, as expressed in Eqs. (9), and n decreases as f z increases.
When f z is small and n is large, the blade is fed in the x-direction after cutting in the y-direction is completed with each rotation of the blade. However, when f z is large and n is small, cutting in the y-direction is not completed when the blade makes one rotation, and cutting in the x-direction is performed simultaneously. Hence, σ xx increases as f z increases. Therefore, in Mode 3, with f z = 0.2 mm, n was insufficiently predicted, resulting in almost the same amount of cutting in the x-and y-directions each time the blade rotated, and σ xx and σ yy had close distributions. This implies that the values of σ xx and σ yy are predicted to be dependent on the amount of cutting in each stress component direction. Thus, σ xx was larger in the deeper range than in Modes 1 and 2 because the amount of cutting in the shear direction was inevitably larger in Mode 3. In other words, the greater the amount of cutting of each stress component, the greater z and magnitude of the √ shape peak, and the deeper the RS distribution. As shown in Fig. 4, the milled surface roughness worsened as f z increased. A compressive RS is preferred because it suppresses crack propagation. However, considering surface roughness, f z should not increase.
To investigate the effect of v c , Modes 1 and 4 were compared, as shown in Fig. 7a and d, where f z was constant at 0.04 mm and only v c was changed. The RS z ranges for all stress components were z ≤ 0.06 mm for Mode 1 with v c = 200 m/min and z ≤ 0.04 mm for Mode 4 with v c = 450 m/min, indicating that the larger v c , the shallower the range of RS. Although the maximum depth of compressive RS in the √ shape for both σ xx and σ yy remained the same at z = 0.01 mm, the values were significantly different: − 112.1 ± 6.7 MPa for σ xx and − 152.4 ± 6.7 MPa for σ yy in Mode 1, and − 77.9 ± 5.4 MPa for σ xx and − 100.0 ± 5.4 MPa for σ yy in Mode 4. Tang et al. showed that the tangential force decreases with increasing v c and that the axial force peaks at a certain v c [36]. As a result, σ xx and σ yy RS components change to tensile RS at z = 0 mm with increasing v c . Denkena also showed that only the tangential force decreases as v c increases, whereas the axial force remains constant [15], suggesting that the compressive RS value may be smaller and its range shallower because of the decrease in the cutting force as v c increases. However, Mode 4 is considered a better milling condition because machining is completed at a speed 2.25 times faster than in Mode 1, resulting in higher machining efficiency and small surface roughness, as shown in Fig. 4.

Effect of the milling parameters on the RS
in-plane distribution using the 2D method  Fig. 7 are indicated by the thin horizontal lines in Fig. 8. Focusing on the difference in collimator d c , a RS distribution with a similar trend was obtained for all stress components for both collimators, clearly indicating the existence of an RS distribution in the plane direction as well as a change with z. Nevertheless, a relatively smooth sinusoidal distribution was observed for the collimator with d c = 0.8 mm, while several extreme values appeared for the collimator with d c = 0.146 mm. All RS measurement methods using XRD provided averaged diffraction information from the grains in the X-ray-irradiated area. A larger d c reflects diffraction information from a larger number of crystals, providing results with less variation, whereas a smaller d c provides a higher spatial resolution. The AA7050-T7451 alloy changed from 1-mm pancake-shaped grains to 1-5-μm fine grains because of the frictional heat from strong machining [37]. Therefore, even with a collimator with d c = 0.146 mm, sufficient diffraction information from the crystals is expected to be obtained because grain refinement caused by strong friction occurs on the machined surface during milling. As evidence, the standard deviation values of the measurement results for d c = 0.146 mm in Fig. 8  at three locations on the milled surface using the sin 2 ψ method [11]. They found that when the axial z of the cut was 0.1 mm, the compressive RS decreased with increasing f z for v c ≤ 350 m/min and increased with increasing f z for v c ≥ 350 m/min [11]. Thus, it was reported that the RS distribution on the milled surface depends on the milling parameters. However, for the Mode 4 condition of v c = 450 m/min in this study, the opposite result was observed. Nevertheless, as shown in Fig. 8, the RS values vary significantly depending on the measurement location, even in the same plane. Therefore, to obtain an accurate trend, it is necessary to show the results for each measurement location rather than the average value. The compressive RS of σ yy tended to be smaller for y > 0 mm and greater for y < 0 mm regardless of the milling conditions. The compressive RS value was particularly small under Modes 3 and 4, indicating that the σ yy stress component decreased as the f z threshold was exceeded or as v c increased. Finally, σ xy showed a minimum value at − 0.5 mm ≤ y ≤ 0 mm for all milling conditions, and the RS decreased as f z increased, but no significant difference in RS distribution occurred with changes in v c . Figure 9a-d show the y-axis direction RS distributions in the xy-plane near z, where the maximum compressive RS occurred, for samples milled in Modes 1-4. The error bars indicate the standard deviations. The evaluated xyplanes were exposed by electropolishing, and the RS in the xy-planes of z = 20, 20, 50, and 14 μm were evaluated for Modes 1-4. The average RS values for each stress component at each z shown in Fig. 7 are indicated by thin horizontal lines in Fig. 9. Compared with Fig. 8, which shows the measurement results at z = 0 mm, many evaluation points have significantly different values owing to the difference in the d c . In the xy-plane at z = 0 mm, the tool was in direct contact with the sample; therefore, the in-plane RS distribution trend remained almost the same, even when d c was different. By contrast, as shown in Fig. 9, the RS distribution changed in a complex manner as the z increased from the milled surface, and the collimator with d c = 0.8 mm could not detect the change, whereas the collimator with d c = 0.146 mm, which had a higher spatial resolution, showed that the RS distribution significantly. The results in Fig. 8 show that the measurement at d c = 0.146 mm is accurate. Comparing the standard deviations in Fig. 9, the standard deviations of d c = 0.146 mm is larger than d c = 0.8 mm, showing that the local RS values in the xy-plane with z > 0 vary. Moreover, some regularities were observed in the RS distribution on the machined surface, as shown in Fig. 8, and it is difficult to detect this regularity, as shown in Fig. 9. Therefore, the degrees of tensile RS that could cause fatigue fractures at each z were compared. Figure 10a and b show the maximum tensile RSs in the xy-plane at z = 0 mm shown in Fig. 8 for collimators with d c = 0.146 and 0.8 mm, respectively, and the error bars indicate the standard deviations. Figure 10c and d show the maximum tensile RSs in the xy-plane near the z-axis, where the maximum compressive RSs are shown in Fig. 9 for collimators with d c = 0.146 and 0.8 mm.
Initially, the effect of v c at z = 0, as in Fig. 10a and b, showed that for both collimators, σ xx remained almost unchanged with increasing v c , while σ yy changed significantly to the tensile side, and σ xy became slightly smaller. The effect of f z was almost the same for σ xx as f z increased, whereas σ yy showed a large tensile change and σ xy showed a slightly decreasing trend. When a collimator with d c = 0.146 mm was used, σ xx showed tensile RS for all conditions, while σ yy showed tensile RS only in Mode 4. For the collimator with d c = 0.8 mm, σ xx and σ yy are compressive RSs. In the near-z region, as shown in Fig. 10c and d, where the maximum compressive RS occurred, the changing trend in values differed for each collimator. For the collimator with d c = 0.146 mm, σ xx remained almost unchanged with increasing v c , while σ yy changed slightly to the tensile side, and σ xy became much smaller. The effects of f z , σ xx , and σ yy showed a maximum compressive RS at f z = 0.1 mm. It then changed significantly to the tensile side at f z = 0.2 mm, and σ xy decreased monotonically. The values of σ xx and σ yy showed tensile RS only in Mode 4. By contrast, for the collimator with d c = 0.8 mm, σ xx changed slightly to the tensile side, and σ yy changed significantly to the tensile side with increasing v c , whereas σ xy became smaller. For the effect of f z, σ xx showed a minimum compressive RS at f z = 0.1 mm and then changed significantly to the compressive side at f z = 0.2 mm, while σ yy changed monotonically to the tensile side with increasing f z , and σ xy remained almost unchanged. As in the case of z = 0 mm, σ xx and σ yy are the compressive RS components. These results show that in the xy-plane at the z where the maximum compressive RS occurred, the effects of the different collimators were significant, especially for σ xx at f z = 0.2 mm. These differences result from the difference in spatial resolution for the collimators, and the collimator with d c = 0.8 mm averaged RSs over a larger In summary, tensile RSs are more likely to occur at σ xx on a surface with z = 0 mm, and tensile RSs that cause fatigue failure are more likely to occur as v c and f z increase, regardless of the z. The variation in the RS distribution in the plane direction is directly related to the distortion of thin-plate materials. The machining efficiency is better with larger v c and f z ; however, it is necessary to determine the milling conditions that consider the surface precision and RS distribution to improve the fatigue performance. If compressive RS can be introduced by adjusting the milling conditions to suppress the initiation and propagation of fatigue cracks without generating tensile RSs that cause fatigue, aircraft safety can be improved.