Next Article in Journal
Aluminium Matrix Composite Materials Reinforced by 3D-Printed Ceramic Preforms
Previous Article in Journal
Non-Destructive Evaluation of Thermal Aging in EPDM Rubber Using Electromagnetic Techniques
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Discrete-Finite Element Analysis Model Based on Almen Intensity Test for Evaluation of Real Shot Peening Residual Stress

1
Graduate School, Department of Mechanical and Aeronautical System Engineering, Cheongju University, 298 Daeseong-ro, Cheongju-si 28503, Republic of Korea
2
Department of Aeronautical and Mechanical Engineering, Cheongju University, 298 Daeseong-ro, Cheongju-si 28503, Republic of Korea
*
Author to whom correspondence should be addressed.
Materials 2023, 16(15), 5472; https://doi.org/10.3390/ma16155472
Submission received: 4 July 2023 / Revised: 30 July 2023 / Accepted: 1 August 2023 / Published: 4 August 2023

Abstract

:
In this study, a combined discrete-finite element model based on the Almen intensity measurement test was proposed to evaluate the real shot peening residual stress. The discrete element analysis was utilized to simulate the random behavior of numerous shot balls, while the finite element analysis was employed to quantitatively predict the residual stress induced by shot peening. Moreover, the Almen intensity, an essential factor in the actual shot peening process, was taken into account. Initially, an Almen strip analysis model was established, and the multi-random impact analysis was performed to validate the good agreement between the analytical Almen curve and experimental Almen curve. Subsequently, the unit cell discrete-finite element analysis model was expanded for predicting the peening residual stress, incorporating the Almen intensity. The analysis results showed a significant correlation between the predicted peening residual stress and the XRD (X-ray diffraction) experimental residual stress. Therefore, it was confirmed that the proposed discrete-finite element random impact analysis model in this study could serve as an effective analytical technique capable of substituting for the actual shot peening process.

1. Introduction

Shot peening technology is widely used in various industrial fields such as aircraft, industrial power equipment components, railway vehicles, and automotive durable parts [1,2,3]. Shot peening is a surface treatment technique that induces compressive residual stress by subjecting the surface of a component or metal material to the continuous multiple impacts of numerous small shot particles, causing plastic deformation. Therefore, when shot peening is applied to the surface of a mechanical component subjected to repetitive tensile stress, the compressive residual stress resists the tensile stress and extends the fatigue life of the component [4,5,6]. To quantitatively analyze the effects of such compressive residual stress, it is necessary to predict the magnitude of the surface and maximum compressive residual stress, as well as the depth of the residual stress layer. Typically, peening residual stress generated on the surface of fatigue-prone components is measured using X-ray diffraction [7,8,9]. However, this method poses potential harm to the operator due to the use of radiation. Moreover, the preparation and measurement times for the measurement vary depending on the operator’s skill level, and the measurement results are not consistent, making it difficult to quantitatively predict residual stress. For these reasons, various studies based on finite element analysis have been conducted to quantitatively predict the multi-impact shot peening residual stress.
Recently, numerical analysis studies utilizing three-dimensional random impact have been introduced to simulate the phenomenon of actual shot particles randomly impacting the material surface [10,11,12,13,14,15,16,17,18,19,20,21]. Kim et al. and Meguid et al. proposed a four-fold symmetric cell analysis model that replaces a wide peening area and allows for various impact sequences perpendicular to the four vertices of the symmetric cell [10,11]. Miao et al. proposed a hexahedral analysis model that focuses on localized areas with impact areas ranging from 1/6 to 1/25 of the total area for multi-impact impact [12]. Additionally, Kim et al. and Gangaraj et al. proposed a multi-impact analysis model in the form of a disk-type plate, taking into account the effect of peening coverage. They simulated the phenomenon of random impact by concentrating the impact of shot balls in the central region of the disk [13,14]. Xiao et al. established a finite element model with shot balls randomly impacting an overall area that is approximately 2.7 times the length of a single shot ball’s diameter [15]. However, these previous studies all had in common that the positions of the shot balls were pre-determined to impact the material surface. In other words, they used random number generation functions in software such as Excel, Matlab, or Python to obtain the impact and height coordinates of the shot balls and then arranged them accordingly. Therefore, these papers have limitations in implementing the random impact positions and sequences of shot balls. Furthermore, due to the increase in analysis time caused by the deformation of shot balls upon impact with the material, the peening area and the quantity of shot balls are often limited in most cases. Therefore, recently, a new analysis technique combining discrete element analysis and finite element analysis has been introduced to simulate actual random phenomena [22,23,24,25,26,27,28]. Tu et al. first obtained the random velocities and impact positions of shot balls via discrete element analysis and then transferred them to dynamic finite element analysis to implement the multi-impact phenomenon [22]. They obtained residual stress at different depths from the surface at the impact centers of the shot balls and averaged them arithmetically. However, actual peening residual stress measured by X-ray diffraction exhibits variations in stress not only at the impact center but also in the surrounding area and the entire surface. Therefore, analytically evaluated residual stress based solely on the impact center differs from XRD experimental residual stress. Bhuvaraghan et al. performed single-impact finite element analysis to obtain pressure data from the indentation generated on the material surface [23]. They incorporated this pressure into a multi-shot random impact discrete element analysis to obtain equivalent pressure exerted on the material at the impact locations, which were also randomly obtained. Then, they obtained the peening residual stress from a finite element analysis model that incorporated only the equivalent pressure and randomly obtained impact positions instead of using a multi-shot model. However, these studies did not include experimental verification to establish the validity of the analytical techniques. Additionally, they did not consider the essential factor of peening intensity, which is crucial in predicting the peening residual stress of actual fatigue-prone components. Edward et al. first performed dynamic single-shot impact finite element analysis to calculate the restoration coefficient, which was then incorporated into the discrete element analysis [24]. After the discrete element analysis, they obtained the contact load exerted by the shot ball on the material surface and incorporated it back into the finite element analysis model to simulate the peening phenomenon even in the absence of contact between the shot ball and the material. However, they did not simulate the multi-shot impact phenomenon based on actual peening processes, and process variables such as Almen intensity and peening coverage were not considered. Furthermore, experimental residual stress that could confirm the validity of the analytical residual stress was excluded from their study. Indeed, the recent studies mentioned earlier have incorporated advanced analytical techniques in simulating actual peening phenomena. However, it is true that they have not considered the crucial factor of Almen intensity, which represents the peening intensity and encompasses variables such as the impact velocity and quantity of shot balls. Almen intensity is a key process parameter that must be taken into account during the peening process. Fatigue-prone components have diverse fatigue characteristics depending on their usage and operating environment, so it is important to apply the optimal peening intensity [29,30,31,32,33]. Therefore, the multi-random impact peening analysis model in this study, which considers Almen intensity, stands out from the existing analytical models. In this study, a new analysis technique is proposed to evaluate the real shot peening residual stress using the discrete-finite element analysis model based on the Almen intensity measurement test. Here, the discrete element analysis is used for the random behavior of numerous shot balls, and the impact load is obtained when the shot balls impact the surface of the model, while finite element analysis is used for the deformation behavior of the model, and residual stress is indicated. Ultimately, this technique will be verified by comparing the analytical residual stress with the X-ray diffraction (XRD). This approach can be effectively utilized as a valuable analytical tool for exploring optimal shot peening conditions in various durability components.

2. Almen Intensity Simulation and Experimental Verification

2.1. Almen Intensity by DE-FE Analysis

In this section, an Almen strip analysis model was developed using discrete-finite element analysis to obtain the Almen intensity corresponding to the actual experimental results and assess its suitability. In shot peening processes, the Almen intensity is a crucial factor that must be taken into account. It combines various variables such as peening time, peening coverage, shot diameter, shot velocity, and impact angle into a single parameter that represents the Almen intensity [29,30,31,32,33]. Therefore, even when simulating the actual shot peening phenomenon, a random multi-impact analysis based on the Almen intensity must be performed. To secure the Almen intensity, an Almen strip analysis model was established to obtain the analytical Almen intensity. The Almen strip analysis model targeted an A-type strip with dimensions of 76 mm in length, 19 mm in width, and 1.29 mm in thickness, made of SAE1070 material, as shown in Figure 1. The mechanical properties of SAE1070 material were considered in the analysis model, as shown in Table 1.
In this study, the Johnson–Cook model, as shown in Equation (1), was incorporated to simulate the deformation behavior that occurs rapidly due to plastic deformation when the shot impacts the material surface. This is because the material experiences a high strain rate during the impact process. The Johnson–Cook model offers the advantages of easy simulation of deformation behavior and simple structure, making it computationally efficient [34].
σ = [ A + B ( ε p ) n ] 1 + C l n ε ˙ p ε ˙ 0 1 T T r T m T r m
The Johnson–Cook model includes several parameters: A represents the initial yield stress, B is the strain hardening coefficient, n is the strain hardening exponent, C is the strain rate sensitivity coefficient, ε p is equivalent plastic strain, and m is the thermal softening exponent. ε ˙ p is the equivalent plastic strain rate, ε ˙ 0 is the initial equivalent plastic strain rate, T is working temperature, Tr is the room temperature, and Tm is the melting temperature. The values for these parameters are listed in Table 1.
To reduce the analysis time, a quarter-symmetric Almen strip analysis model, as shown in Figure 2a, was used. The boundary conditions for the model were as follows: the endpoints were restrained in the x, y, and z directions (Ux = Uy = Uz = 0), and displacement was constrained in the direction perpendicular to the symmetry plane (Ux = 0, Uy = 0) [32]. To validate the quarter-symmetric analysis model, a full analysis model was also established and analyzed for comparison. The full analysis model, depicted in Figure 2b, had all four endpoints restrained in the x, y, and z directions (Ux = Uy = Uz = 0).
Figure 3 illustrates the phenomenon of random impact of shot balls in the Almen strip analysis model. To simulate the random impact of shot balls, a probability density function was employed using Abaqus software. The shot balls were assumed to be rigid bodies with a diameter of 0.8 mm and a density of 7800 kg/m3. Furthermore, it was assumed that the shot balls undergo no deformation during the impacts.
Figure 4a shows the deformation of the 1/4 symmetry model and the full model of the Almen strip after multiple impacts. Figure 4b displays the deformation in the z-direction at points P1 to P5 along the length (y-direction) of the Almen strip, as indicated in Figure 4. P1 represents the endpoint in the length direction (y-direction) of the Almen strip, while P5 represents the central region along the length direction (y-direction) of the strip. The shot ball’s impact velocity was set at 60 m/s, and the impact time was 0.3 ms.
After analysis, it was confirmed that the deformation obtained from the 1/4 symmetric model closely approximated the deformation obtained from the full model in Figure 5. Furthermore, the analysis time required for the 1/4 symmetric model was 82 min, while the analysis time for the full model was 475 min, resulting in a reduction of 393 min, or approximately 82%, in analysis time using the 1/4 symmetric model. Based on these findings, the 1/4 symmetric model was used as the multi-random impact analysis model in this study.
The minimum element sizes are 0.09 mm, 0.08 mm, 0.07 mm, 0.06 mm, and so on. These element sizes were used in Figure 6 to examine the convergence of the deformation measured at the center position of the Almen strip analysis model. In Figure 7, the average deformation obtained after analysis was approximately 1.314 mm, and all four element sizes provided values that were close to the average deformation. Considering the reduction in analysis time, a minimum element size (Le) of 0.09 mm was chosen.
Figure 8 shows the deformations of the Almen strip in a 1/4 symmetric analysis model with a minimum element size (Le) of 0.09 mm after multiple shot balls impact randomly with velocities of 40 m/s, 50 m/s, 60 m/s, and 70 m/s for a duration of 0.3 ms. From Figure 8, it can be observed that the maximum deformations are 0.674 mm at an impact velocity of 40 m/s, 0.907 mm at 50 m/s, 1.306 mm at 60 m/s, and 1.407 mm at 70 m/s.
Table 2 presents the arc heights obtained from the Almen strip analysis model as a function of the analysis time. These analytical arc heights are calculated based on the maximum deformation of the Almen strip obtained from the analysis, following the experimental measurement method for determining arc height. Specifically, the height of the deflection within the 31.75 mm gauge length of the Almen gauge is calculated to represent the analytical arc height. These analytical arc heights are compared with the experimental arc heights to validate the effectiveness of the proposed discrete-finite element multiple-random shot peening impact analysis in this study. The analysis times are set as 0.3 ms, 0.6 ms, 0.9 ms, and 1.2 ms, with corresponding impact velocities of 40 m/s, 50 m/s, 60 m/s, and 70 m/s for each analysis time.

2.2. Experimental Effectiveness of the Almen Intensity Simulation

In order to validate the accuracy of the previously performed discrete-finite element (DE-FE) analysis based on the Almen intensity analysis model, an Almen intensity measurement experiment was conducted in this section. An impeller-type shot peening equipment was used in the experiment as shown in Figure 9, with a diameter of 250 mm. The rotation speed of the impeller was adjusted to achieve shot velocities ranging from 40 m/s to 70 m/s. The Almen strips used in the experiment were of type A, identical to the ones used in the analysis model. The shot used was a round-cut wire shot with an average diameter of 0.8 mm. The shot peening duration was set to 60 s, 120 s, 180 s, and 240 s, with a projection rate of 42 kg/min.
Figure 10 illustrates the process of measuring the deflection height of the curved Almen strip using the Almen gauge after shot peening. The Almen intensity is measured by placing the Almen strip on the upper surface of the Almen gauge and using a dial gauge to measure the deflection height within the 31.75 mm span (L) of the gauge [32,35].
Table 3 presents the measured values of Almen intensity for the Almen strip. It can be observed that the Almen intensity increases with the shot peening time for each impact velocity, and it also increases with the increasing impact velocity for the same peening time. This can be attributed to the fact that an increase in impact velocity results in higher kinetic energy, and an increase in peening time leads to the accumulation of kinetic energy, both of which contribute to the deformation of the Almen strip in proportion.
By comparing the analytical Almen intensity values with the experimental Almen intensity values, we can validate the validity of the analysis technique. Figure 11 shows the comparison between the analytical and experimental Almen intensities. The analytical Almen intensity curve closely matches the experimental Almen intensity curve. For all impact velocities, the analysis time of 0.3 ms corresponds well to the experimental peening time of 60 s, 0.6 ms corresponds to 120 s, 0.9 ms corresponds to 180 s, and 1.2 ms corresponds to 240 s. This means that the analytical Almen intensity curves can be used as a substitute for the experimentally obtained Almen intensity curves, which are essential for actual peening processes. These analytical Almen intensity curves are employed for determining the Almen intensity and residual stress prediction in the shot peening process of real metal materials, as well as in the unit cell analysis models.

3. Analytical Peening Residual Stress and Its Experimental Verification

3.1. Peening Residual Stress by the DE-FE Analysis

Figure 12 represents a DE-FE multi-random impact unit cell analysis model used to predict residual stresses induced in the metal material after actual shot peening processes. The model is capable of representing the entire surface area. The material used in the model is AISI4340, with dimensions of 1 mm × 1 mm × 3 mm. The same constant probability density function and Johnson–Cook model, as introduced in Section 2, were applied to the Almen strip analysis models. The shot ball model reflects a diameter of 0.8 mm and a density of 7800 kg/m3, assuming it is an undeformed rigid body. These shot balls impact at a 90° angle, with impact velocities of 40 m/s, 50 m/s, 60 m/s, and 70 m/s. Contact or impacts between the shot balls are not considered. Table 4 presents the mechanical properties and constant values for the AISI 4340 material and the shot balls.
During the shot peening process, the continuous impacts of shot balls create surface indentations, and significant deformation occurs beneath the impact area. To analyze the deformation of the elasto-plastic specimen, the NLGEOM (non-linear Geometry) option in the ABAQUS Explicit code was utilized. Additionally, a three-dimensional 8-node reduced integration element (C3D8R) was used for meshing [35,36,37]. As for the boundary conditions, the bottom surface of the unit cell was fully constrained (Ux = Uy = Uz = 0), and the four side surfaces were restrained in the direction perpendicular to the surface (Ux = 0, Uz = 0). To incorporate the penalty algorithm for contact, contact surface elements were placed on both sides of the material and shot balls. The analysis involved approximately 12,650 elements and 13,968 nodes. The compressive residual stresses generated by shot peening are primarily embedded in the surface and subsurface regions. Therefore, the finite element mesh was refined in the areas with significant plastic deformation on the surface and coarser in the subsurface regions with smaller plastic deformation. The minimum element size (Le) adopted for the analysis was 0.02 mm. To evaluate the residual stresses from the proposed analytical model in this study, an analytical residual stress evaluation method based on XRD experiments, commonly used for measuring residual stresses, is required. Generally, XRD residual stresses are obtained as averaged stresses over a certain area where X-rays are probed. Taking this into consideration, in Figure 12, the evaluation regions were defined into three groups with different numbers of nodes to obtain the analytical average solutions that closely match the experimental field [38]. These evaluation regions consist of 3, 9, and 25 nodes, respectively, and are centered around the analytical model. The evaluation region composed of three nodes is formed by consecutive nodes aligned with the direction of residual stress measurement. The evaluation region composed of nine nodes is a 3 × 3 area, horizontally and vertically, resulting in a total of nine nodes. The evaluation region composed of 25 nodes is a 5 × 5 area, horizontally and vertically, representing the entire set of nodes.

3.2. Experimental Verification of Peening Residual Stress

In order to validate the effectiveness of the proposed unit cell analysis model shown in Figure 12, an analysis of multiple random peening residual stresses was performed, considering the actual measured Almen intensity. The analysis considered the average over the evaluation regions, which refers to the multi-node average for the three evaluation regions shown in Figure 12. After the analysis, the obtained peening residual stresses were compared with the experimentally measured residual stresses by XRD in the depth direction from the surface of specimens peened with four different Almen intensities [39] in Figure 13. In case (a), the Almen intensity is 0.0027 inchA (=0.06858 mmA), and the impact velocity is 40 m/s, with an analytical peening time of 0.13188 ms. After the analysis, it can be observed that both the 9-point and 25-point averaged solutions are closer to the XRD experimental residual stress than the 3-point one. Case (a) represents the scenario where the maximum compressive residual stress occurs at the surface. Considering both the surface and maximum compressive residual stress, the 25-point average solution is slightly closer to the XRD experimental results than the 9-point average solution. In case (b), the Almen intensity is 0.0063 inchA (=0.16002 mmA), and the impact velocity is 50 m/s, with an analytical peening time of 0.25949 ms. After the analysis, it is observed that the 25-point averaged solution is closer to the XRD experimental residual stress than the 3-point and 9-point averaged solutions. In case (c), the Almen intensity is 0.0083 inchA (=0.21082 mmA), and the impact velocity is 60 m/s, with an analytical peening time of 0.25502 ms. After the analysis, it is found that the 3-point, 9-point, and 25-point averaged solutions are close to the XRD experimental residual stress. In case (d), the Almen intensity is 0.0141 inchA (=0.35814 mmA), and the impact velocity is 70 m/s, with an analytical peening time of 0.60531 ms. After the analysis, it is observed that both the 9-point and 25-point averaged solutions are close to the XRD experimental residual stress, but the 25-point average solution is closer to the XRD experimental residual stress at different depths. Overall, based on the comparison with XRD experimental residual stress [39], the 25-point averaged solution is found to be more effective in all cases. Based on the results shown in Figure 13, the discrete-finite element model that incorporates the actual Almen intensity and provides the 25-point averaged solution is chosen as the optimal peening analysis model.
The effectiveness of shot peening is typically evaluated based on the magnitude of surface residual stress and maximum compressive residual stress. Table 5 presents the magnitudes of analytical and experimental SRCS, along with the corresponding errors. The error ranges from 0.75% to 6.55%, indicating a good agreement between them. Additionally, Table 6 shows the magnitudes of analytical and experimental MRCS, along with the corresponding errors. The error ranges from 1.58% to 7.93%, demonstrating an even better agreement compared to the surface residual stress. The reason for the larger error in surface residual stress compared to maximum compressive residual stress is attributed to the sensitivity of the experimental measurement positions and directions to the irregular shape of surface indentation generated by randomly impacting shot-peening balls on the material surface. In addition, it is believed that the analysis model’s surface grid, which constitutes the surface region of the model, undergoes irregular deformation due to the random impacts of shot peening balls, resulting in a greater stress interference effect compared to the interior grid of the analysis model. Figure 14 illustrates the distribution of residual stress within the analysis model after shot peening analysis. It shows a trend where the compressive residual stress increases and then decreases as the depth from the surface increases. This trend is consistent with the distribution of measured residual stress typically observed after conventional shot peening processes.

4. Conclusions

In this study, a new combined discrete-finite element analysis model incorporating both discrete elements and finite elements was used to simulate the residual stress distribution in actual shot peening processes. The conclusion is as follows.
A discrete-finite element analysis model was established for the Almen strips used in measuring Almen intensity, allowing us to obtain the deformation quantities via the analysis of multiple random impacts. The discrete element analysis simulated the random behavior of numerous shot-peening balls, while the finite element analysis verified the deformation behavior of the Almen strips. It was also observed that the deformation increased proportionally with the increasing impact velocity of the shot balls. By converting these deformations into Arc heights, the analytical Almen curve was obtained. Additionally, the experimental Almen curve based on the actual peening time was measured and matched with the analytical curve. The results showed a good agreement between the analytical and experimental Almen curves for different peening times. This confirmed that the analytical Almen curve could effectively replace the experimental Almen curve based on the peening time. Afterward, the analysis model was expanded to predict the actual residual stresses in shot peening. This analysis model was constructed using unit cells, and an analysis technique was applied where the stresses obtained at multiple points within the effective area were averaged to evaluate the area-averaged stresses. In this study, the analysis was performed on AISI 4340 material that was shot peened under four different Almen intensities, obtaining the residual stresses. Then, these analysis results were compared with the experimental residual stresses obtained from specimens shot peened with the same Almen intensities. The results showed that the surface residual stresses with an error range from 0.75% to 6.55% and the maximum compressive residual stresses with an error range from 1.58% to 7.93% obtained from the analysis were in close agreement with the experimental residual stresses.
Ultimately, the proposed discrete-finite element random impact analysis model, which incorporates Almen intensity-dependent arc height, has been validated for its effectiveness and validity as an analysis technique that can replace actual shot peening processes. This analysis technique can be effectively used in various industries to apply durable mechanical parts with optimal peening intensity.

Author Contributions

Conceptualization and methodology, C.W.; investigation and review, C.W. and Y.P.; computational analysis, C.W.; validation, Y.P. and T.K.; Writing—original draft, C.W.; Writing-review, Y.P. and T.K. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the National Research Foundation of Korea (NRF) grant funded by the Korean government (MSIT) (NRF-2022R1FA1074178).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hu, D.; Tian, T.; Wang, X.; Mao, J.; Xiao, Z.; Wang, R. Surface hardening analysis for shot peened GH4720Li superalloy using a DEM-FEM coupling RV simulation method. Int. J. Mech. Sci. 2021, 209, 106689. [Google Scholar] [CrossRef]
  2. Xia, B.; Wang, B.; Zhang, P.; Ren, C.; Duan, Q.; Li, X.; Zhang, Z. Improving the high-cycle fatigue life of a high-strength spring steel for automobiles by suitable shot peening and heat treatment. Int. J. Fatigue 2022, 161, 106891. [Google Scholar] [CrossRef]
  3. Gan, J.; Gao, Z.A.; Wang, Y.; Wang, Z.; Wu, W. Small-Scale Experimental Investigation of Fatigue Performance Improvement of Ship Hatch Corner with Shot Peening Treatments by Considering Residual Stress Relaxation. J. Mar. Sci. Eng. 2021, 9, 419. [Google Scholar] [CrossRef]
  4. Xie, X.; Zhang, L.; Zhu, L.; Li, Y.; Hong, T.; Yang, W.; Shan, X. State of the Art and Perspectives on Surface-Strengthening Process and Associated Mechanisms by Shot Peening. Coatings 2023, 13, 859. [Google Scholar] [CrossRef]
  5. Bagherifard, S. Enhancing the Structural Performance of Lightweight Metals by Shot Peeing. Adv. Eng. Mater. 2019, 21, 1801140. [Google Scholar] [CrossRef]
  6. Bianchetti, C.; Delbergue, D.; Bocher, P.; Lévesque, M.; Brochu, M. Analytical fatigue life prediction of shot peened AA 7050-T7451. Int. J. Fatigue 2019, 118, 271–281. [Google Scholar] [CrossRef]
  7. Lin, J.; Ma, N.; Lei, Y.; Murakawa, H. Measurement of residual stress in arc welded lap joints by cosα X-ray diffraction method. J. Mater. Process. Technol. 2017, 243, 387–394. [Google Scholar] [CrossRef]
  8. Maleki, E.; Farrahi, G.H.; Reza Kashyzadeh, K.; Unal, O.; Gugaliano, M.; Bagherifard, S. Effects of Conventional and Severe Shot Peening on Residual Stress and Fatigue Strength of Steel ALSI 1060 and Residual Stress Relaxation Due to Fatigue Loading: Experimnetal and Numerical Simulation. Met. Mater. Int. 2021, 27, 2575–2591. [Google Scholar] [CrossRef]
  9. Zhan, K.; Wu, Y.; Li, J.; Zhao, B.; Yan, Y.; Xie, L.; Ji, V. Investigation on surface layer characteristics of shot peened graphene reinforced Al composite by X-ray diffraction method. Appl. Surf. Sci. 2018, 435, 1257–1264. [Google Scholar] [CrossRef]
  10. Kim, T.; Lee, H.; Hyun, H.C.; Jung, S. A simple but effective FE model with plastic shot for evaluation of peening residual stress and its experimental validation. Mater. Sci. Eng. A 2011, 528, 5945–5954. [Google Scholar] [CrossRef]
  11. Meguid, S.A.; Shagal, G.; Stranart, J.C. 3D FE analysis of peening of strain-rate sensitive materials using multiple impingement model. Int. J. Impact Eng. 2002, 27, 119–134. [Google Scholar]
  12. Miao, H.Y.; Larose, S.; Perron, C.; Lévesque, M. On the potential applications of a 3D random finite element model for the simulation of shot peening. Adv. Eng. Softw. 2009, 40, 1023–1038. [Google Scholar] [CrossRef]
  13. Kim, T.; Lee, H.; Jung, S.; Lee, J.H. A 3D model with plastic shot for evaluation of equi-biaxial peening residual stress due to multi-impacts. Surf. Coat. Technol. 2012, 206, 3125–3136. [Google Scholar] [CrossRef]
  14. Gangaraj, S.M.H.; Guagliano, M.; Farrahi, G.H. An approach to relate shot peening finite element simulation to the actual coverage. Surf. Coat. Technol. 2014, 243, 39–45. [Google Scholar] [CrossRef]
  15. Xiao, X.; Tong, X.; Liu, Y.; Zhao, R.; Gao, G.; Li, Y. Prediction of shot peen forming effects with single and repeated impacts. Int. J. Mech. Sci. 2018, 137, 182–194. [Google Scholar]
  16. Lin, Q.; Liu, H.; Zhu, C.; Parker, R.G. Investigation on the effect of shot peening coverage on the surface integrity. Appl. Surf. Sci. 2019, 489, 66–72. [Google Scholar]
  17. Majzoobi, G.H.; Azizi, R.; Nia, A.A. A three-dimensional simulation of shot peening process using multiple shot impacts. J. Mater. Process. Technol. 2005, 164, 1226–1234. [Google Scholar] [CrossRef]
  18. Bagherifard, S.; Ghelichi, R.; Guagliano, M. On the shot peening surface coverage and its assessment by mens of finite element simulation: A critical review and some original developments. Appl. Surf. Sci. 2012, 259, 186–194. [Google Scholar] [CrossRef]
  19. Wang, C.; Wang, L.; Wang, X.; Xu, Y. Numerical study of grain refinement induced by severe shot peeing. Int. J. Mech. Sci. 2018, 146, 280–294. [Google Scholar] [CrossRef]
  20. Bagherifard, S.; Ghelichi, R.; Guagliano, M. Numerical and experimental analysis of surface roughness by shot peening. Appl. Surf. Sci. 2012, 258, 6831–6840. [Google Scholar] [CrossRef]
  21. Zhao, J.; Tang, J.; Ding, H.; Shao, W.; Zhao, X.; Liu, H.; Jiang, T. A numerical and experimental investigation on the evolution of three-dimensional surface topography of 12Cr2Ni4A steel in shot peeing. J. Manuf. Process. 2021, 70, 259–270. [Google Scholar] [CrossRef]
  22. Tu, F.; Delbergue, D.; Miao, H.; Klotz, T.; Brochu, M.; Bocher, P.; Levesque, M. A sequential DEM-FEM coupling method for shot peening simulation. Surf. Coat. Technol. 2017, 319, 200–212. [Google Scholar] [CrossRef]
  23. Bhuvaraghan, B.; Srinivasan, S.M.; Maffeo, B.; McCLain, R.D.; Potdar, Y.; Prakash, O. Shot peening simulation using discrete and finite element methods. Adv. Eng. Softw. 2010, 41, 1266–1276. [Google Scholar] [CrossRef]
  24. Edward, A.B.; Heyns, P.S.; Kok, S. A Numerical Investigation of a Single-Shot in DEM-FEM Approach to Shot Peening Simulation. Metals 2019, 9, 1183. [Google Scholar] [CrossRef] [Green Version]
  25. Jebahi, M.; Gakwaya, A.; Lévesque, J.; Mechri, O.; Ba, K. Robust methodology to simulate real shot peeing precess using discrete-continuum coupling method. Int. J. Mech. Sci. 2016, 107, 21–33. [Google Scholar]
  26. Zhao, J.; Tang, J.; Zhou, W.; Jiang, T.; Liu, H.; Xing, B. Numerical modeling and experimental verification of residual stress distribution evolution of 12Cr2Ni4A steel generated by shot peening. Surf. Coat. Technol. 2022, 430, 127993. [Google Scholar]
  27. Murugaratnam, K.; Utili, S.; Petrinic, N. A combined DEM-FEM numerical method for Shot Peening parameter optimization. Adv. Eng. Softw. 2015, 79, 13–26. [Google Scholar] [CrossRef]
  28. Zhang, J.; Lu, S.; Wu, T.; Zhou, Z.; Zhang, W. An evaluation on SP surface property by means of combined FEM-DEM shot dynamics simulation. Adv. Eng. Softw. 2018, 115, 283–296. [Google Scholar] [CrossRef]
  29. Zhang, Y.; Lai, F.; Qu, S.; Ji, V.; Liu, H.; Li, X. Effect of shot peening on residual stress distribution and tribological behaviors of 17Cr2Ni2MoVNb steel. Surf. Coat. Technol. 2020, 386, 125497. [Google Scholar] [CrossRef]
  30. Guagliano, M. Relating Almen intensity to residual stresses induced by shot peening: A numerical approach. J. Mater. Process. Technol. 2001, 110, 277–286. [Google Scholar] [CrossRef]
  31. Gariépy, A.; Miao, H.Y.; Lévesque, M. Simulation of the shot peeing process with variable shot diameters and impacting velocities. Adv. Eng. Softw. 2017, 114, 121–133. [Google Scholar] [CrossRef]
  32. Wang, C.; Li, W.; Jiang, J.; Chao, X.; Zeng, W.; Xu, J.; Yang, J. An Improved Approach to Direct Simulation of an Actual Almen Shot Peening Intensity Test with a Large Number of Shots. Materials 2020, 13, 5088. [Google Scholar] [CrossRef] [PubMed]
  33. Yang, Z.; Lee, Y.; He, S.; Jia, W.; Zhao, J. Analysis of the Influence of High Peening Coverage on Almen Intensity and Residual Compressive Stress. Appl. Sci. 2019, 10, 105. [Google Scholar] [CrossRef] [Green Version]
  34. Huang, H.; Wang, Z.; Gan, J.; Yang, Y.; Wang, X.; He, J.; Gan, X. The study of universality of a method for predicting surface nanocrystallization after high energy shot peening based on finite element analysis. Surf. Coat. Technol. 2019, 358, 617–627. [Google Scholar] [CrossRef]
  35. Gariépy, A.; Larose, S.; Perron, C.; Lévesque, M. Shot peening and peen forming finite element modelling—Towards a quantitatie mothod. Int. J. Solids Struct. 2011, 48, 2859–2877. [Google Scholar] [CrossRef]
  36. AMahmoudi, A.H.; Ghasemi, A.; Farrahi, G.H.; Sherafatnia, K. A comprehensive experimental and numerical study on redistribution of residual stresses by shot peeing. Mater. Des. 2016, 90, 478–487. [Google Scholar] [CrossRef]
  37. Lin, Q.; Liu, H.; Zhu, C.; Chen, D.; Zhou, S. Effects of different shot peening parameters on residual stress, surface roughness and cell size. Surf. Coat. Technol. 2020, 398, 126054. [Google Scholar] [CrossRef]
  38. Kim, T.; Lee, J.H.; Lee, H.; Cheong, S.K. An area-average approach to peening residual stress under multi-impacts using a three-dimensional symmetry-cell finite element model with plastic shots. Mater. Des. 2010, 31, 50–59. [Google Scholar] [CrossRef]
  39. Torres, M.A.S.; Voorwald, H.J.C. An evaluation of shot peening, residual stress and stress relaxation on the fatigue life of AISI 4340 steel. Int. J. Fatigue 2002, 24, 877–886. [Google Scholar] [CrossRef]
Figure 1. Sample. (a) geometric dimensions; (b) an Almen strip of type A.
Figure 1. Sample. (a) geometric dimensions; (b) an Almen strip of type A.
Materials 16 05472 g001
Figure 2. Boundary conditions for the quarter Almen strip model and the whole Almen strip model. (a) quarter model; (b) full model.
Figure 2. Boundary conditions for the quarter Almen strip model and the whole Almen strip model. (a) quarter model; (b) full model.
Materials 16 05472 g002
Figure 3. The 3D discrete element analysis model with multi-shot random impact on the surface of the Almen strip. (a) detail view; (b) overall view.
Figure 3. The 3D discrete element analysis model with multi-shot random impact on the surface of the Almen strip. (a) detail view; (b) overall view.
Materials 16 05472 g003
Figure 4. Deflection of the Almen strip. (a) deflection of the quarter Almen strip; (b) deflection of the full Almen strip.
Figure 4. Deflection of the Almen strip. (a) deflection of the quarter Almen strip; (b) deflection of the full Almen strip.
Materials 16 05472 g004
Figure 5. Deflections for the full model and quarter model of the Almen strip.
Figure 5. Deflections for the full model and quarter model of the Almen strip.
Materials 16 05472 g005
Figure 6. The various element minimum sizes of the Almen strip DE-FE analysis model.
Figure 6. The various element minimum sizes of the Almen strip DE-FE analysis model.
Materials 16 05472 g006
Figure 7. Convergence of arc height with various mesh sizes.
Figure 7. Convergence of arc height with various mesh sizes.
Materials 16 05472 g007
Figure 8. The deformation of the Almen strip in the z-direction, when (a) the peening velocity is 40 m/s; (b) the peening velocity is 50 m/s; (c) the peening velocity is 60 m/s; (d) the peening velocity is 70 m/s.
Figure 8. The deformation of the Almen strip in the z-direction, when (a) the peening velocity is 40 m/s; (b) the peening velocity is 50 m/s; (c) the peening velocity is 60 m/s; (d) the peening velocity is 70 m/s.
Materials 16 05472 g008
Figure 9. Shot peening equipment.
Figure 9. Shot peening equipment.
Materials 16 05472 g009
Figure 10. Measurement devices. (a) Almen gage; (b) Arc height calculation principle.
Figure 10. Measurement devices. (a) Almen gage; (b) Arc height calculation principle.
Materials 16 05472 g010
Figure 11. Arc height of the Almen strip at different impact velocities (a) with an impact velocity of 40 m/s; (b) with an impact velocity of 50 m/s; (c) with an impact velocity of 60 m/s; (d) with an impact velocity of 70 m/s.
Figure 11. Arc height of the Almen strip at different impact velocities (a) with an impact velocity of 40 m/s; (b) with an impact velocity of 50 m/s; (c) with an impact velocity of 60 m/s; (d) with an impact velocity of 70 m/s.
Materials 16 05472 g011
Figure 12. The concept of the area-averaged (=full node averaged) residual stress from unit cell FE model.
Figure 12. The concept of the area-averaged (=full node averaged) residual stress from unit cell FE model.
Materials 16 05472 g012
Figure 13. Comparison of the DE-FE analytical residual stress and XRD test residual stress [39] with four different Almen intensities with (a) the Almen intensity is 0.069 mmA (v = 40 m/s); (b) the Almen intensity is 0.160 mmA (v = 50 m/s); (c) the Almen intensity is 0.211 mmA (v = 60 m/s); (d) the Almen intensity is 0.358 mmA (v = 70 m/s).
Figure 13. Comparison of the DE-FE analytical residual stress and XRD test residual stress [39] with four different Almen intensities with (a) the Almen intensity is 0.069 mmA (v = 40 m/s); (b) the Almen intensity is 0.160 mmA (v = 50 m/s); (c) the Almen intensity is 0.211 mmA (v = 60 m/s); (d) the Almen intensity is 0.358 mmA (v = 70 m/s).
Materials 16 05472 g013aMaterials 16 05472 g013b
Figure 14. Different surface morphologies at different peening velocities (a) with an impact velocity of 40 m/s; (b) with an impact velocity of 50 m/s; (c) with an impact velocity of 60 m/s; (d) with an impact velocity of 70 m/s.
Figure 14. Different surface morphologies at different peening velocities (a) with an impact velocity of 40 m/s; (b) with an impact velocity of 50 m/s; (c) with an impact velocity of 60 m/s; (d) with an impact velocity of 70 m/s.
Materials 16 05472 g014
Table 1. Coefficients of Johnson-Cook model adopted for deformation analysis of Almen strip [10].
Table 1. Coefficients of Johnson-Cook model adopted for deformation analysis of Almen strip [10].
MaterialE
(GPa)
Poisson’s Ratioρ
(kg/m3)
A
(MPa)
B
(MPa)
nmTm
(K)
Tr
(K)
C
SAE 10702050.2978001408600.80.2341.017932980.0134
Shot ball2100.307800
Table 2. Arc height by DE-FE analysis.
Table 2. Arc height by DE-FE analysis.
Analysis Time
t (ms)
Arc Height
(mmA)
40 m/s50 m/s60 m/s70 m/s
0.3 ms0.1240.1610.2430.266
0.6 ms0.1980.2210.3220.371
0.9 ms0.2330.2540.3350.401
1.2 ms0.2430.2930.3610.438
Table 3. Arc heights by Almen test.
Table 3. Arc heights by Almen test.
Exposure Time
t (s)
Arc Height (mmA)
v = 40 m/sv = 50 m/sv = 60 m/sv = 70 m/s
600.1430.1850.2460.275
1200.1970.2380.3120.355
1800.2350.2740.3360.390
2400.2520.3040.3720.416
Table 4. Coefficients of Johnson-Cook model adopted for residual stress analysis of the unit cell DE-FE model. [10].
Table 4. Coefficients of Johnson-Cook model adopted for residual stress analysis of the unit cell DE-FE model. [10].
MaterialE
(GPa)
Poisson’s Ratioρ
(kg/m3)
A
(MPa)
B
(MPa)
nmTm
(K)
Tr
(K)
C
AISI43402100.3078001498943.80.2601.0317932980.014
Shot ball2100.307800
Table 5. Comparison of the numerical and experimental results of the surface residual compressive stress (SRCS).
Table 5. Comparison of the numerical and experimental results of the surface residual compressive stress (SRCS).
Shot Ball Velocity
(m/s)
Arc Height
(mmA)
Analysis
σsrcs/MPa
XRD Test [39]
σsrcs/MPa
Error
(%)
400.069−922.986−929.9610.75
500.160−874.086−917.514.73
600.211−911.864−975.8286.55
700.358−730.179−781.256.53
Table 6. Comparison of the numerical and experimental results of the maximum residual compressive stress (MRCS).
Table 6. Comparison of the numerical and experimental results of the maximum residual compressive stress (MRCS).
Shot Ball Velocity
(m/s)
Arc Height
(mmA)
Analysis
σsrcs/MPa
XRD Test [39]
σsrcs/MPa
Error
(%)
400.069−944.677−929.9611.58
500.160−1041.03−1010.892.98
600.211−1067.78−1159.847.93
700.358−1211.59−1256.253.55
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wang, C.; Park, Y.; Kim, T. A Discrete-Finite Element Analysis Model Based on Almen Intensity Test for Evaluation of Real Shot Peening Residual Stress. Materials 2023, 16, 5472. https://doi.org/10.3390/ma16155472

AMA Style

Wang C, Park Y, Kim T. A Discrete-Finite Element Analysis Model Based on Almen Intensity Test for Evaluation of Real Shot Peening Residual Stress. Materials. 2023; 16(15):5472. https://doi.org/10.3390/ma16155472

Chicago/Turabian Style

Wang, Chengan, Yujin Park, and Taehyung Kim. 2023. "A Discrete-Finite Element Analysis Model Based on Almen Intensity Test for Evaluation of Real Shot Peening Residual Stress" Materials 16, no. 15: 5472. https://doi.org/10.3390/ma16155472

APA Style

Wang, C., Park, Y., & Kim, T. (2023). A Discrete-Finite Element Analysis Model Based on Almen Intensity Test for Evaluation of Real Shot Peening Residual Stress. Materials, 16(15), 5472. https://doi.org/10.3390/ma16155472

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