Acessibilidade / Reportar erro

A Note on the Inverse Reconstruction of Residual Fields in Surface Peened Plates

Abstract

A modified stress function approach is developed here to reconstruct induced stress, residual stress and eigenstrain fields from limited experimental measurements. The present approach is successfully applied to three experimental measurements set in surface peened plates with shallow shot peening affected zone. The well-rehearsed advantage of the proposed approach is that it not only minimizes the deviation of measurements from its approximations but also will result in an inverse solution satisfying a full range of continuum mechanics requirements. Also, the effect of component thickness as a geometric parameter influencing the residual stress state is comprehensively studied. A key finding of present study is that the plate thickness has no influence on the maximum magnitude of eigenstrain profile and compressive residual stresses within the shot peening affected zone while having a great influence on the magnitude of tensile residual stress and the gradient of linear residual stresses present in deeper regions.

Keywords:
Surface peening; Residual stress; Induced Stress; Inverse problems; Eigenstrain

1 INTRODUCTION

Shot peening is a mechanical surface treatment generally utilized to improve the stress corrosion resistance and enhance the fatigue life of engineering structures (Meguid, 1991Meguid, S.A., (1991). Effect of partial coverage upon the fatigue fracture behavior of peened components, Fatigue and Fracture of Engineering Materials and Structures 14(5):515-30. DOI: 10.1111/j.1460-2695.1991.tb00680.x.
https://doi.org/10.1111/j.1460-2695.1991...
). It is well known that the shot peening process induces compressive residual stresses and work hardening in near-surface regions of treated components in order to suppress and decelerate crack initiation and propagation in engineering components subjected to cyclic loading. Therefore predicting the distribution and magnitude of residual stresses generated by shot peening is of utmost importance in understanding the influence of physical parameters of the process (Fathallah et al., 1998Fathallah, R., Inglebert, G., Castex, L., (1998). Prediction of plastic deformation and residual stresses induced in metallic parts by shot peening, Journal of Materials Science & Technology 14(7):631-639. DOI: 10.1179/mst.1998.14.7.631.
https://doi.org/10.1179/mst.1998.14.7.63...
).

A remarkable set of publications could be found in the literature dedicated to the modeling of residual stress state of shot peening process utilizing different types of finite element modeling (Guagliano, 2001Guagliano, M., (2001). Relating Almen intensity to residual stresses induced by shot peening: a numerical approach, Journal of Materials Processing Technology 110(3):277-286. DOI: 10.1016/S0924-0136(00)00893-1.
https://doi.org/10.1016/S0924-0136(00)00...
; Zimmermann et al., 2009Zimmermann, M., Klemenz, M., Schulze, V., Lohe, D., (2009). Numerical studies on the influence of thickness on the residual stress development during shot peening, High Performance Computing in Science and Engineering '08. Springer, Berlin, pp. 481-492. DOI: 10.1007/978-3-540-88303-6_34.
https://doi.org/10.1007/978-3-540-88303-...
; Gonzales et al., 2010Gonzales, M.A.C., Barrios, D.B., de Lima, N.B., Goncalves, E., (2010). Importance of considering a material micro-failure criterion in the numerical modeling of the shot peening process applied to parabolic leaf springs, Latin American Journal of Solids and Structures 7(1):21-40.), analytical methods (Shen and Atluri, 2006Shen, S., Atluri, S.N., (2006). An analytical model for shot-peening induced residual stresses, CMC: Computers, Materials, & Continua 4(2):75-86. DOI: 10.3970/cmc.2006.004.075.
https://doi.org/10.3970/cmc.2006.004.075...
; Franchim et al., 2009Franchim, A.S., de Campos, V.S., Travessa, D.N., de Moura Neto, C., (2009). Analytical modeling for residual stresses produced by shot peening, Materials & Design 30(5):1556-60. DOI: 10.1016/j.matdes.2008.07.040.
https://doi.org/10.1016/j.matdes.2008.07...
; Miao et al., 2010Miao, H.Y., Larose, S., Perron, C., Lévesque, M., (2010). An analytical approach to relate shot peening parameters to Almen intensity, Surface and Coatings Technology 205(7):2055-2066. DOI: 10.1016/j.surfcoat.2010.08.105.
https://doi.org/10.1016/j.surfcoat.2010....
) or experimental methods (Niku-Lari et al., 1985Niku-Lari, A., Lu, J., Flavenot, J.F., (1985). Measurement of Residual-Stress Distribution by the Incremental Hole-Drilling Method, Journal of Mechanical Working Technology 11(2):167-188. DOI: 10.1016/0378-3804(85)90023-3.
https://doi.org/10.1016/0378-3804(85)900...
; Fathallah et al., 1998Fathallah, R., Inglebert, G., Castex, L., (1998). Prediction of plastic deformation and residual stresses induced in metallic parts by shot peening, Journal of Materials Science & Technology 14(7):631-639. DOI: 10.1179/mst.1998.14.7.631.
https://doi.org/10.1179/mst.1998.14.7.63...
; Menig et al., 2001Menig, R., Pintschovius, L., Schulzea, V., Vohringer, O., (2001). Depth profiles of macro residual stresses in thin shot peened steel plates determined by X-ray and neutron diffraction, Scripta Materialia 45(8):977-983. DOI: 10.1016/S1359-6462(01)01063-6.
https://doi.org/10.1016/S1359-6462(01)01...
). However it is currently difficult and prohibitively time-consuming to measure the complete patterns of residual stresses throughout the most structural materials (Krawitz, 2011Krawitz, A.D., (2011). Neutron strain measurement, Materials Science and Technology 27(3):589-603. DOI: 10.1179/1743284710Y.0000000029.
https://doi.org/10.1179/1743284710Y.0000...
). Consequently most measurements are restricted and hence resulted in relatively incomplete experimental information. So, there is great practical value in inferring the complete residual stress field from a limited number of measurements (Coules et al., 2014Coules, H.E., Smith, D.J., Venkata, K.A., Truman, C.E., (2014). A method for reconstruction of residual stress fields from measurements made in an incompatible region, International Journal of Solids and Structures 51(10):1980-1990. DOI: 10.1016/j.ijsolstr.2014.02.008.
https://doi.org/10.1016/j.ijsolstr.2014....
).

Theory of inverse method for determination of residual stress field from limited experimental measurements has been developed over the past decade and gained a great attention. One such method known as variational eigenstrain method is based on the theory of eigenstrains. The variational eigenstrain method utilized a combination of finite element analysis and distributed basis eigenstrains that combines experimental characterization in terms of residual elastic strains (Korsunsky et al., 2007Korsunsky, A.M., Regino, G.M., Nowell, D., (2007). Variational eigenstrain analysis of residual stresses in a welded plate, International Journal of Solids and Structures 44(13):4574-4591. DOI: 10.1016/j.ijsolstr.2006.11.037.
https://doi.org/10.1016/j.ijsolstr.2006....
; Korsunsky, 2009Korsunsky, A.M., (2009). Eigenstrain analysis of residual strains and stresses, Journal of Strain Analysis for Engineering Design 44(1):29-43. DOI: 10.1243/03093247JSA423.
https://doi.org/10.1243/03093247JSA423...
; Jun and Korsunsky, 2010Jun, T.-S., Korsunsky, A.M., (2010). Evaluation of residual stresses and strains using the eigenstrain reconstruction method, International Journal of Solids and Structures 47(13):1678-1686. DOI: 10.1016/j.ijsolstr.2010.03.002.
https://doi.org/10.1016/j.ijsolstr.2010....
). A smoothed inverse eigenstrain method is also developed for reconstruction of residual field from limited strain measurements that allows suppressing fluctuations that are contrary to the physics of the problem (Faghidian, 2014Faghidian, S.A., (2014). A smoothed inverse eigenstrain method for reconstruction of the regularized residual fields, International Journal of Solids and Structures 51(25-26):4427-4434. DOI: 10.1016/j.ijsolstr.2014.09.012.
https://doi.org/10.1016/j.ijsolstr.2014....
). An alternative method is the stress function approach which does not utilize assumed eigenstrain distribution. Stress function approach approximate the residual stress field using a series of stress functions which directly solve the stress equilibrium equations together with the traction free boundary conditions (Farrahi et al., 2009aFarrahi, G.H., Faghidian, S.A., Smith, D.J., (2009a) Reconstruction of residual stresses in autofrettaged thick-walled tubes from limited measurements, International Journal of Pressure Vessels and Piping 86(11):777-784. DOI: 10.1016/j.ijpvp.2009.03.010.
https://doi.org/10.1016/j.ijpvp.2009.03....
, 2010Farrahi, G.H., Faghidian, S.A., Smith, D.J., (2010). An inverse method for reconstruction of the residual stress field in welded plates, Journal of Pressure Vessel Technology, Trans. ASME 132(6):612051-612059. DOI: 10.1115/1.4001268.
https://doi.org/10.1115/1.4001268...
; Faghidian et al., 2012Faghidian, S.A., Goudar, D., Farrahi, G.H., Smith, D.J., (2012). Measurement, analysis and reconstruction of residual stresses, Journal of Strain Analysis for Engineering Design 47(4):254-264. DOI: 10.1177/0309324712441146.
https://doi.org/10.1177/0309324712441146...
). Reconstruction of residual stress distribution due to surface peening was first attempted by (Korsunsky, 2005Korsunsky, A.M., (2005). On the modeling of residual stresses due to surface peening using eigenstrain distributions, Journal of Strain Analysis for Engineering Design 40(8):817-824. DOI: 10.1243/030932405X30984.
https://doi.org/10.1243/030932405X30984...
) utilizing the variational eigenstrain approach. Afterward based on the stress function approach, axi-symmetric residual stress field was reconstructed in shot-blasted steel round bars (Farrahi et al., 2009bFarrahi, G.H., Faghidian, S.A., Smith, D.J., (2009b). An inverse approach to determination of residual stresses induced by shot peening in round bars, International Journal of Mechanical Sciences 51(9-10):726-731. DOI: 10.1016/j.ijmecsci.2009.08.004.
https://doi.org/10.1016/j.ijmecsci.2009....
). Also (Jun et al., 2011Jun, T.-S., Venter, A.M., Korsunsky, A.M., (2011). Inverse eigenstrain analysis of the effect of non-uniform sample shape on the residual stress due to shot peening, Experimental Mechanics 51(2):165-174. DOI: 10.1007/s11340-010-9346-2.
https://doi.org/10.1007/s11340-010-9346-...
) analyzed residual elastic strains in non-uniformly shaped shot-peened samples based on the variational eigenstrain approach. Recently based on the stress function approach, (Faghidian, 2015Faghidian, S.A., (2015). Inverse determination of the regularized residual stress and eigenstrain fields due to surface peening, Journal of Strain Analysis for Engineering Design 50(2):84-91. DOI: 10.1177/0309324714558326.
https://doi.org/10.1177/0309324714558326...
) reconstruct the residual fields form experimental measurements in surface peened plates with deep shot peening affected zone utilizing Cauchy type modulation functions. Also (Faghidian, 2015Faghidian, S.A., (2015). Inverse determination of the regularized residual stress and eigenstrain fields due to surface peening, Journal of Strain Analysis for Engineering Design 50(2):84-91. DOI: 10.1177/0309324714558326.
https://doi.org/10.1177/0309324714558326...
) performed a full discussion on various mathematical aspects of numerical reconstruction including the regularity issue of the approximated solution by Tikhonov-Morozov regularization method.

In present study, the stress function approach is modified to reconstruct the residual stress and eigenstrain field form experimental residual stress measurements in surface peened plates with shallow shot peening affected zone. The proposed approach results in a continuous reconstructed residual stresses and eigenstrains while satisfying all of the continuum mechanics requirements. In contrast to stress function approach where the residual stress field are directly determined from Airy stress function, in present approach the induced stress field would be first reconstructed utilizing Gaussian type modulation function. Subsequently the residual stress distribution and eigenstrain field are inversely determined satisfying the static balance across the domain and strain compatibility equations, respectively. Furthermore, the effect of component thickness as a geometric parameter influencing the residual stress state and eigenstrain distribution would be examined.

2 DETERMINATION OF RESIDUAL FIELDS

2.1 Residual stresses

Residual stresses are the stress field supported by a continuum in a fixed reference configuration where there is no external forces and thermal gradients (Hoger, 1986Hoger, A., (1986). On the determination of residual stress in an elastic body, Journal of Elasticity 16(3):303-324. DOI: 10.1007/BF00040818.
https://doi.org/10.1007/BF00040818...
). The domain of interest containing the distribution of residual stresses due to surface peening is a plate with depth of h in the z direction, while coordinates origin coincides on the un-peened surface. The inherent symmetry of the shot peening process under the conditions of normal impact, results in equal residual stress distribution in x and y directions (Cao et al., 1995Cao, W., Fathallah, R., Castex, L., (1995). Correlation of Almen arc height with residual stresses in shot peening process, Materials Science and Technology 11(9):967-973. DOI: 10.1179/mst.1995.11.9.967.
https://doi.org/10.1179/mst.1995.11.9.96...
). Thus residual stress field is supposed to be independent of the x and y coordinates and its unknown variation in the z direction is assumed to be,

The residual stress field must satisfy equilibrium equations in absence of body forces and traction free boundary conditions at all boundaries (Timoshenko and Goodier, 1970Timoshenko, S.P. and Goodier, J.N., (1970). Theory of elasticity, third ed. McGraw-Hill (New York).). Furthermore Cartesian components of mean residual stress must be always zero (Mura, 1987Mura, T., (1987). Micromechanics of defects in solids, Kulwer Academic Publishers (Dordrecht).). So zero mean residual stress for the introduced domain is reduced to,

In present approach, rather than direct determination of residual stress field from Airy stress function, the residual stress distribution is supposed to be,

where represents the depth of the shot peening affected zone and would be determined later based on the continuity condition of the reconstructed residual stress profile. Also σInd (z) and σinr (z) are the induced stress and linear residual stress distribution defined within and outside of shot peening affected zone, respectively. In shot peening process, induced stresses are those encountered in a fully constrained component that create stretching and bending of the part to reach a balanced residual stress state (Van Luchene et al., 1995Van Luchene, R.D., Johnson, J., Carpenter, R.G., (1995). Induced stress relationships for wing skin forming by shot peening, Journal of Materials Engineering and Performance 4(3):283-290. DOI: 10.1007/BF02649065.
https://doi.org/10.1007/BF02649065...
). In the proposed approach, both of the stretching and bending effects are regarded as the linear residual stresses.

The residual stresses measured in a thick component can be used to estimate induced stresses since the large thickness does not allow significant bending and stretching (Gariépy et al., 2011Gariépy, A., Larose, S., Perron, C., Lévesque, M., (2011). Shot peening and peen forming finite element modeling-Towards a quantitative method, International Journal of Solids and Structures 48(20):2859-2877. DOI: 10.1016/j.ijsolstr.2011.06.003.
https://doi.org/10.1016/j.ijsolstr.2011....
). Therefore induced stresses should have force balance across the domain of thick parts and its distribution is supposed to be,

Although the force equilibrium across the entire domain are obviously satisfied by the proposed form of the induced stresses as Eq. (4) but the moment static balance are not satisfied yet. To control the numerical behavior of the approximate solution and to ensure its existence and uniqueness, an arbitrary function of F(z)is assumed to have the following asymptotic expansion,

Where λm are unknown coefficients to be determined later based on the minimization of the deviation of measurements from its approximations and the smooth shape functions of fm(z) should exhibit the expected distribution of induced stresses in surface peening process. However it is well-known that the compressive residual stresses in the near-surface region have to be balanced by tensile residual stresses in the bulk of the specimen. Also the induced stress distribution is supposed to be concentrated in a relatively thin layer adjacent to the layer directly influenced by shot peening (Gariépy et al., 2011Gariépy, A., Larose, S., Perron, C., Lévesque, M., (2011). Shot peening and peen forming finite element modeling-Towards a quantitative method, International Journal of Solids and Structures 48(20):2859-2877. DOI: 10.1016/j.ijsolstr.2011.06.003.
https://doi.org/10.1016/j.ijsolstr.2011....
). Thus to impose the physical conditions of the shot peening induced stresses to shape functions, a family of shape functions is derived from a modulation function of Γ(z) by translation and dilation as fm(z) = Γ(ζm(z - h)). The Gaussian type modulation function of Γ(z) = exp(yz) could be an appropriate choice that results in the expected distribution of induced stresses in the surface peening process. Therefore the explicit form of the shape functions are given by,

where the positive real constants of 1 < ζ < y govern the rate of convergence of the asymptotic solution and would be determined numerically based on the well-posedness of the solution considering relatively small condition number of the system. Accordingly, the final form of the induced stress would be expressed as a truncated series with M terms by,

The proposed form of shape functions will result in induced stress distribution exhibiting a compressive stress with the maximum value occurring at a distinct depth below the shot peened surface while concentrated in a relatively thin layer adjacent to the layer directly influenced by shot peening. Also the linear residual stress σInr (z) = c0 + c1z is defined outside the shot peening affected zone and unknown coefficients of c0, c1would be determined by satisfying the requirements of force and moment balance across the thickness of the surface peened plate as,

Moreover in contrast to conventional reconstruction methods that assume the depth of the shot peening zone, it would be determined here based on the continuity condition of the reconstructed residual stress profile at the boundary of shot peening affected zone, i.e. σlnr() = σind().

Thus, the proposed approach results in a continuous reconstructed residual stress field while clearly satisfies stress equilibrium equations and traction free boundary conditions. Also the requirement of vanishing mean residual stress is guaranteed by satisfying the force and moment balance across the thickness of the surface peened plate.

2.2 Eigenstrains

The residual stress field in a continuum is not uniquely defined by the geometry, constitutive relations and boundary conditions as is the case for stress fields which arise solely from external loading (Timoshenko and Goodier, 1970Timoshenko, S.P. and Goodier, J.N., (1970). Theory of elasticity, third ed. McGraw-Hill (New York).). Instead, the residual stresses are dependent on the internal incompatibility within the continuum where the distribution of incompatible strain field in such a continuum is known as the eigenstrain distribution. Thus total strain tensor εij can be expressed as the sum of elastic eij and eigenstrain terms as εij = eij+ (Mura, 1987Mura, T., (1987). Micromechanics of defects in solids, Kulwer Academic Publishers (Dordrecht).). Moreover kinematic analysis of continuous deformation requires that the total strain field to be compatible. Hence the total strain compatibility equation for equi-biaxial strains in the x and y directions, εxx(z) = εyy (z) = ε(z), is reduced to (Timoshenko and Goodier, 1970Timoshenko, S.P. and Goodier, J.N., (1970). Theory of elasticity, third ed. McGraw-Hill (New York).),

As a result, total strain should be a linear function of z to satisfy the strain compatibility equation. Motivated from fundamental assumptions of variational eigenstrain approach (Korsunsky, 2005Korsunsky, A.M., (2005). On the modeling of residual stresses due to surface peening using eigenstrain distributions, Journal of Strain Analysis for Engineering Design 40(8):817-824. DOI: 10.1243/030932405X30984.
https://doi.org/10.1243/030932405X30984...
); it is assumed that residual stress components are proportional to elastic strains and also the eigenstrain field outside the shot peening affected zone vanishes. Therefore the difference between the residual stress profile and linear residual stress is proportional to the eigenstrain field,

where the elastic coefficient of = E/(1 - υ) may be readily found for biaxial stress state of shot peening process (Fathallah et al., 1998Fathallah, R., Inglebert, G., Castex, L., (1998). Prediction of plastic deformation and residual stresses induced in metallic parts by shot peening, Journal of Materials Science & Technology 14(7):631-639. DOI: 10.1179/mst.1998.14.7.631.
https://doi.org/10.1179/mst.1998.14.7.63...
) and E and v are Young's modulus and Poisson's ratio of the material, respectively. Furthermore in the case of extremely thick plates where the shot peening affected zone is extremely small compared to the thickness of the plate, linear residual stress approaches to zero and eigenstrain distribution is given by,

The direct proportionality of residual stress and permanent plastic strain in extremely thick plates subjected to uniform shot peening is in total agreement with variational eigenstrain approach (Korsunsky, 2005Korsunsky, A.M., (2005). On the modeling of residual stresses due to surface peening using eigenstrain distributions, Journal of Strain Analysis for Engineering Design 40(8):817-824. DOI: 10.1243/030932405X30984.
https://doi.org/10.1243/030932405X30984...
). In Section 3, the effects of component thickness as a geometric parameter influencing the residual stress state would be examined comprehensively. As a final point it should be emphasized that linearity of total strains is due to satisfy strain compatibility equations and is valid independent of plate bending theory assumptions.

2.3 Least squares analysis

Provided a limited set of experimental measurements consisting of the values of residual stress as Tn are measured at positions zn, then evaluating the shape function of φm(z) at each of n measurement positions, results in predicted values of φmn = φm (zn). However the unknown coefficients λm appeared in the asymptotic expansion of induced stress field cannot be directly evaluated by solving the over-determined set of linear algebraic equations of Tn = λmφmn. Since the number of truncated series M used to approximate the induced stress field should be less than the number of measurement points N, to attain a fast convergent series.

A standard alternative approach, as described earlier in (Faghidian, 2014Faghidian, S.A., (2014). A smoothed inverse eigenstrain method for reconstruction of the regularized residual fields, International Journal of Solids and Structures 51(25-26):4427-4434. DOI: 10.1016/j.ijsolstr.2014.09.012.
https://doi.org/10.1016/j.ijsolstr.2014....
, 2015Faghidian, S.A., (2015). Inverse determination of the regularized residual stress and eigenstrain fields due to surface peening, Journal of Strain Analysis for Engineering Design 50(2):84-91. DOI: 10.1177/0309324714558326.
https://doi.org/10.1177/0309324714558326...
), is to find the best values of coefficients λm that minimizes the deviation of the experimental measurements from its predicted values. Therefore the following norm of error should be minimized,

It may be shown that unique values of coefficients Λ = {λm}can be determined by setting the gradient of error norm with respect to the coefficients equal to zero, ∂e/∂λm = 0 (Faghidian, 2015Faghidian, S.A., (2015). Inverse determination of the regularized residual stress and eigenstrain fields due to surface peening, Journal of Strain Analysis for Engineering Design 50(2):84-91. DOI: 10.1177/0309324714558326.
https://doi.org/10.1177/0309324714558326...
). Therefore final form of the coefficients may be determined as,

where the matrix notation of Φ = [φmn] and T = {Tn} are introduced over the n=1...N measurement points and m=1...M the number of truncated series used to approximate the induced stress field.

Furthermore, a comprehensive discussion on the various mathematical aspects of the numerical reconstruction consisting of invertibility, uniqueness, well-posedness and convergence of asymptotic residual fields could be found in former works (Faghidian, 2014Faghidian, S.A., (2014). A smoothed inverse eigenstrain method for reconstruction of the regularized residual fields, International Journal of Solids and Structures 51(25-26):4427-4434. DOI: 10.1016/j.ijsolstr.2014.09.012.
https://doi.org/10.1016/j.ijsolstr.2014....
). It should be noted that the continuous dependence of the approximate solution on measured data is shown for shape functions derived from Gaussian type modulation function and the non-linear smooth shape functions are known to be numerically stabilized (Faghidian, 2014Faghidian, S.A., (2014). A smoothed inverse eigenstrain method for reconstruction of the regularized residual fields, International Journal of Solids and Structures 51(25-26):4427-4434. DOI: 10.1016/j.ijsolstr.2014.09.012.
https://doi.org/10.1016/j.ijsolstr.2014....
).

3 ANALYSIS OF EXPERIMENTAL RESULTS AND DISCUSSIONS

The presented framework of inverse determination of induced stress, residual stress and eigenstrain fields would be applied to three experimental data sets here. The experimental data examined here are measured in specimens with shallow shot peening affected zone where the depth of shot peening affected zone is h / 10.

To reconstruct the residual fields, it is assumed that the magnitude of residual stresses is solely known at N measurement points within the shot peening affected zone, reported in Table 1 for every set of measurements. Also in each case studies, values of the parameters ζ and γ appeared in explicit form of shape functions, Eq. (6), are numerically determined based on the well-posedness of the system and are given in Table 1. Once the coefficients λm are evaluated, the induced and linear residual stresses and consequently the eigenstrain distribution over the entire domain would be determined. Moreover for each reconstructed residual field, number of the terms M utilized in truncated series form of residual field is tabulated in Table 1 to examine the rate of convergence issue.

Table 1
Values of shape function parameters, number of measurement points and terms used for residual stress reconstruction.

As the first case study and to illustrate the strength of the present approach, the residual stress and eigenstrain profile determined by variational inverse eigenstrain analysis is examined. The aluminum alloy thick coupon with effective elastic modulus of =103 GPa has dimensions of 150×50×4 mm and in tests is peened to full coverage using typical 0.5 mm shots (Levers and Prior, 1998Levers, A., Prior, A., (1998). Finite element analysis of shot peening, Journal of Materials Processing Technology 80-81:304-308. DOI: 10.1016/S0924-0136(98)00188-5.
https://doi.org/10.1016/S0924-0136(98)00...
). The residual stress and eigenstrain profile determined by (Korsunsky, 2005Korsunsky, A.M., (2005). On the modeling of residual stresses due to surface peening using eigenstrain distributions, Journal of Strain Analysis for Engineering Design 40(8):817-824. DOI: 10.1243/030932405X30984.
https://doi.org/10.1243/030932405X30984...
) utilizing variational inverse eigenstrain analysis are also shown in Fig. 1 and 2 and compared with the results of present approach, respectively. It is obviously inferred from Fig. 1 and Fig. 2 that the reconstructed residual fields by present approach are in excellent agreement with residual stress and eigenstrain profile of variational inverse eigenstrain analysis (Korsunsky, 2005Korsunsky, A.M., (2005). On the modeling of residual stresses due to surface peening using eigenstrain distributions, Journal of Strain Analysis for Engineering Design 40(8):817-824. DOI: 10.1243/030932405X30984.
https://doi.org/10.1243/030932405X30984...
) in the entire shot peening domain.

Figure 1
Residual stress data (Levers and Prior, 1998Levers, A., Prior, A., (1998). Finite element analysis of shot peening, Journal of Materials Processing Technology 80-81:304-308. DOI: 10.1016/S0924-0136(98)00188-5.
https://doi.org/10.1016/S0924-0136(98)00...
), reconstructed induced and residual stress distribution compared with variational inverse eigenstrain results (Korsunsky, 2005Korsunsky, A.M., (2005). On the modeling of residual stresses due to surface peening using eigenstrain distributions, Journal of Strain Analysis for Engineering Design 40(8):817-824. DOI: 10.1243/030932405X30984.
https://doi.org/10.1243/030932405X30984...
).

Figure 2
Reconstructed eigenstrain field compared with eigenstrain profile by variational inverse eigenstrain approach (Korsunsky, 2005Korsunsky, A.M., (2005). On the modeling of residual stresses due to surface peening using eigenstrain distributions, Journal of Strain Analysis for Engineering Design 40(8):817-824. DOI: 10.1243/030932405X30984.
https://doi.org/10.1243/030932405X30984...
).

Also to study the effects of plate thickness on the distribution of eigenstrain and residual stress fields, the X-ray diffraction measurements and FE simulated results of (Zimmermann et al., 2009Zimmermann, M., Klemenz, M., Schulze, V., Lohe, D., (2009). Numerical studies on the influence of thickness on the residual stress development during shot peening, High Performance Computing in Science and Engineering '08. Springer, Berlin, pp. 481-492. DOI: 10.1007/978-3-540-88303-6_34.
https://doi.org/10.1007/978-3-540-88303-...
) is considered. The age hardened IN718 superalloy specimen with Young's modulus of E=200 GPa and Poisson's ratio of ν=0.32 is treated with an air blast machine. The test specimen featuring two sections of different thickness of 5 mm and 20 mm that is shot peened to 98 % coverage with a shot medium of type CCW31 with a measured mean diameter of 0.89 mm, a hardness of 550 HV, an impact angle of 80˚and mean shot velocity of 23 m/s (Zimmermann et al., 2009Zimmermann, M., Klemenz, M., Schulze, V., Lohe, D., (2009). Numerical studies on the influence of thickness on the residual stress development during shot peening, High Performance Computing in Science and Engineering '08. Springer, Berlin, pp. 481-492. DOI: 10.1007/978-3-540-88303-6_34.
https://doi.org/10.1007/978-3-540-88303-...
). The residual stress profiles simulated by FE analysis using an elasto-viscoplastic material model with combined isotropic and kinematic hardening (Zimmermann et al., 2009Zimmermann, M., Klemenz, M., Schulze, V., Lohe, D., (2009). Numerical studies on the influence of thickness on the residual stress development during shot peening, High Performance Computing in Science and Engineering '08. Springer, Berlin, pp. 481-492. DOI: 10.1007/978-3-540-88303-6_34.
https://doi.org/10.1007/978-3-540-88303-...
) are also shown in Fig. 3(a) and 3(b) for thickness of 5 mm and 20 mm, respectively. The reconstructed induced stress, residual stress and eigenstrain distribution by present approach are also exhibited together with X-ray measurements and simulated residual stress depth profiles for different thicknesses in Fig. 3 and Fig. 4, respectively.

Figure 3
Reconstructed induced and residual stress distribution compared with X-ray measurements (Zimmermann et al., 2009Zimmermann, M., Klemenz, M., Schulze, V., Lohe, D., (2009). Numerical studies on the influence of thickness on the residual stress development during shot peening, High Performance Computing in Science and Engineering '08. Springer, Berlin, pp. 481-492. DOI: 10.1007/978-3-540-88303-6_34.
https://doi.org/10.1007/978-3-540-88303-...
) and simulated residual stress profile (Zimmermann et al., 2009Zimmermann, M., Klemenz, M., Schulze, V., Lohe, D., (2009). Numerical studies on the influence of thickness on the residual stress development during shot peening, High Performance Computing in Science and Engineering '08. Springer, Berlin, pp. 481-492. DOI: 10.1007/978-3-540-88303-6_34.
https://doi.org/10.1007/978-3-540-88303-...
) for thickness of (a) 5 mm, (b) 20 mm.

Figure 4
Reconstructed eigenstrain depth profile for X-ray measurements of (Zimmermann et al., 2009Zimmermann, M., Klemenz, M., Schulze, V., Lohe, D., (2009). Numerical studies on the influence of thickness on the residual stress development during shot peening, High Performance Computing in Science and Engineering '08. Springer, Berlin, pp. 481-492. DOI: 10.1007/978-3-540-88303-6_34.
https://doi.org/10.1007/978-3-540-88303-...
) with thickness of 5 mm and 20 mm.

It is clearly inferred from Fig. 3 that decreasing the component thickness does not affect the maximum magnitude of compressive residual stresses in shot peening affected region but leads to higher tensile residual stresses in greater depths. Also the difference in magnitude of tensile stresses and gradient of linear residual stress can be attributed to the different thicknesses available to balance the compressive residual stresses. It should be noted that increasing the thickness of surface peened plates would result in lower magnitudes of linear residual stress which is sometimes less than the resolution of measurement methods; as a result the linear residual stress is not constant through the thickness and still varies linearly even if with small gradient.

Furthermore it is obviously deduced from Fig. 4 that the maximum magnitude of eigenstrain field is approximately the same and identically independent of the plate thickness but the eigenstrain field is distributed over the larger zone for smaller thickness. The results are in complete agreement with FE simulated plastic deformation which the maximum magnitudes are always identically independent of FE model thickness (Zimmermann et al., 2009Zimmermann, M., Klemenz, M., Schulze, V., Lohe, D., (2009). Numerical studies on the influence of thickness on the residual stress development during shot peening, High Performance Computing in Science and Engineering '08. Springer, Berlin, pp. 481-492. DOI: 10.1007/978-3-540-88303-6_34.
https://doi.org/10.1007/978-3-540-88303-...
).

However it worth to note that, although the induced stresses have the same profile as the residual stresses in shot peening affected zone but induced stresses cannot represent the residual field in deeper regions since not only they do not satisfy the moment balance across the plate thickness but also would result in negative eigenstrain distribution near to shot peening affected zone which is not acceptable (Korsunsky, 2009Korsunsky, A.M., (2009). Eigenstrain analysis of residual strains and stresses, Journal of Strain Analysis for Engineering Design 44(1):29-43. DOI: 10.1243/03093247JSA423.
https://doi.org/10.1243/03093247JSA423...
). Thus the compressive residual stresses generated by surface peening procedure in the shot peening affected region could be only balanced by rather low residual stresses varying linearly over the whole cross-section as opposed to relatively high tensile stresses in a near-surface layer.

In reconstruction of residual fields in surface peened plates, a full range of continuum mechanics requirements are satisfied. Furthermore an optimal agreement between experimental measurements and model predictions is attained in the least squares sense. The eigenstrain field is also inversely determined by satisfying strain compatibility equations while in conventional reconstruction methods generally a direct approach would be utilized to determine the residual stresses from an assumed eigenstrain distribution.

The open mathematical problem concerns the development of present approach for appropriate choice of unknown parameters ζ, γ of shape functions fm(z). Although a variety of examples exists for numerical selection of such variable parameters in stress function approach (Faghidian, 2014Faghidian, S.A., (2014). A smoothed inverse eigenstrain method for reconstruction of the regularized residual fields, International Journal of Solids and Structures 51(25-26):4427-4434. DOI: 10.1016/j.ijsolstr.2014.09.012.
https://doi.org/10.1016/j.ijsolstr.2014....
, 2015; Farrahi et al. 2009aFarrahi, G.H., Faghidian, S.A., Smith, D.J., (2009a) Reconstruction of residual stresses in autofrettaged thick-walled tubes from limited measurements, International Journal of Pressure Vessels and Piping 86(11):777-784. DOI: 10.1016/j.ijpvp.2009.03.010.
https://doi.org/10.1016/j.ijpvp.2009.03....
, 2009bFarrahi, G.H., Faghidian, S.A., Smith, D.J., (2009b). An inverse approach to determination of residual stresses induced by shot peening in round bars, International Journal of Mechanical Sciences 51(9-10):726-731. DOI: 10.1016/j.ijmecsci.2009.08.004.
https://doi.org/10.1016/j.ijmecsci.2009....
; Faghidian et al., 2012Faghidian, S.A., Goudar, D., Farrahi, G.H., Smith, D.J., (2012). Measurement, analysis and reconstruction of residual stresses, Journal of Strain Analysis for Engineering Design 47(4):254-264. DOI: 10.1177/0309324712441146.
https://doi.org/10.1177/0309324712441146...
) based on the well-posedness of the reconstructed solution considering relatively small condition number of the system of linear equations (Allaire and Kaber, 2008Allaire, G. and Kaber, S.M., (2008). Numerical linear algebra. Springer (New York).) appearing in the equation of ∂e/∂λm = 2Φ(ΦTΛ - T) = 0 and the fast convergence rate of the solution; but no theoretical results are available at present time.

Nonetheless numerical determination of appropriate values for shape function parameters results in a truncated series with six terms that has a fast convergence. It should be noted that, the convergence criterion for the asymptotic series would be defined as the norm of the discrepancy of successive approximations is of the order of the noise in data. Also a good estimate for the measurement error could be the unbiased estimator for variance of the measured data (Faghidian, 2015Faghidian, S.A., (2015). Inverse determination of the regularized residual stress and eigenstrain fields due to surface peening, Journal of Strain Analysis for Engineering Design 50(2):84-91. DOI: 10.1177/0309324714558326.
https://doi.org/10.1177/0309324714558326...
).

4 CONCLUSIONS

An examination of the eigenstrain field, induced and residual stresses in surface peened specimens with shallow shot peening affected zone, has revealed the points outlined below.

  1. 1

    . The well-rehearsed advantage of the proposed approach is that it not only minimizes the deviation of the residual measurements from its predicted values but also will result in an inverse solution satisfying a full range of continuum mechanics requirements.

  2. 2

    . In contrast to conventional reconstruction methods, the eigenstrain field is analytically determined by satisfying the strain compatibility equation

  3. 3

    . surface region but also depends on the component thickness as an influencing geometric parameter.

  4. 4

    . The plate thickness has no influence on the maximum magnitude of eigenstrain profile within the shot peening affected region but would effectively modify the magnitude and gradient of linear residual stresses present in deeper regions.

  5. 5

    . Although lower linear residual stress would be expected by increasing the thickness of surface peened plate but the residual stress is not constant in the core regions and still varies linearly even if with small gradient.

References

  • Allaire, G. and Kaber, S.M., (2008). Numerical linear algebra. Springer (New York).
  • Cao, W., Fathallah, R., Castex, L., (1995). Correlation of Almen arc height with residual stresses in shot peening process, Materials Science and Technology 11(9):967-973. DOI: 10.1179/mst.1995.11.9.967.
    » https://doi.org/10.1179/mst.1995.11.9.967
  • Coules, H.E., Smith, D.J., Venkata, K.A., Truman, C.E., (2014). A method for reconstruction of residual stress fields from measurements made in an incompatible region, International Journal of Solids and Structures 51(10):1980-1990. DOI: 10.1016/j.ijsolstr.2014.02.008.
    » https://doi.org/10.1016/j.ijsolstr.2014.02.008
  • Faghidian, S.A., Goudar, D., Farrahi, G.H., Smith, D.J., (2012). Measurement, analysis and reconstruction of residual stresses, Journal of Strain Analysis for Engineering Design 47(4):254-264. DOI: 10.1177/0309324712441146.
    » https://doi.org/10.1177/0309324712441146
  • Faghidian, S.A., (2014). A smoothed inverse eigenstrain method for reconstruction of the regularized residual fields, International Journal of Solids and Structures 51(25-26):4427-4434. DOI: 10.1016/j.ijsolstr.2014.09.012.
    » https://doi.org/10.1016/j.ijsolstr.2014.09.012
  • Faghidian, S.A., (2015). Inverse determination of the regularized residual stress and eigenstrain fields due to surface peening, Journal of Strain Analysis for Engineering Design 50(2):84-91. DOI: 10.1177/0309324714558326.
    » https://doi.org/10.1177/0309324714558326
  • Farrahi, G.H., Faghidian, S.A., Smith, D.J., (2009a) Reconstruction of residual stresses in autofrettaged thick-walled tubes from limited measurements, International Journal of Pressure Vessels and Piping 86(11):777-784. DOI: 10.1016/j.ijpvp.2009.03.010.
    » https://doi.org/10.1016/j.ijpvp.2009.03.010
  • Farrahi, G.H., Faghidian, S.A., Smith, D.J., (2009b). An inverse approach to determination of residual stresses induced by shot peening in round bars, International Journal of Mechanical Sciences 51(9-10):726-731. DOI: 10.1016/j.ijmecsci.2009.08.004.
    » https://doi.org/10.1016/j.ijmecsci.2009.08.004
  • Farrahi, G.H., Faghidian, S.A., Smith, D.J., (2010). An inverse method for reconstruction of the residual stress field in welded plates, Journal of Pressure Vessel Technology, Trans. ASME 132(6):612051-612059. DOI: 10.1115/1.4001268.
    » https://doi.org/10.1115/1.4001268
  • Fathallah, R., Inglebert, G., Castex, L., (1998). Prediction of plastic deformation and residual stresses induced in metallic parts by shot peening, Journal of Materials Science & Technology 14(7):631-639. DOI: 10.1179/mst.1998.14.7.631.
    » https://doi.org/10.1179/mst.1998.14.7.631
  • Franchim, A.S., de Campos, V.S., Travessa, D.N., de Moura Neto, C., (2009). Analytical modeling for residual stresses produced by shot peening, Materials & Design 30(5):1556-60. DOI: 10.1016/j.matdes.2008.07.040.
    » https://doi.org/10.1016/j.matdes.2008.07.040
  • Gariépy, A., Larose, S., Perron, C., Lévesque, M., (2011). Shot peening and peen forming finite element modeling-Towards a quantitative method, International Journal of Solids and Structures 48(20):2859-2877. DOI: 10.1016/j.ijsolstr.2011.06.003.
    » https://doi.org/10.1016/j.ijsolstr.2011.06.003
  • Gonzales, M.A.C., Barrios, D.B., de Lima, N.B., Goncalves, E., (2010). Importance of considering a material micro-failure criterion in the numerical modeling of the shot peening process applied to parabolic leaf springs, Latin American Journal of Solids and Structures 7(1):21-40.
  • Guagliano, M., (2001). Relating Almen intensity to residual stresses induced by shot peening: a numerical approach, Journal of Materials Processing Technology 110(3):277-286. DOI: 10.1016/S0924-0136(00)00893-1.
    » https://doi.org/10.1016/S0924-0136(00)00893-1
  • Hoger, A., (1986). On the determination of residual stress in an elastic body, Journal of Elasticity 16(3):303-324. DOI: 10.1007/BF00040818.
    » https://doi.org/10.1007/BF00040818
  • Jun, T.-S., Korsunsky, A.M., (2010). Evaluation of residual stresses and strains using the eigenstrain reconstruction method, International Journal of Solids and Structures 47(13):1678-1686. DOI: 10.1016/j.ijsolstr.2010.03.002.
    » https://doi.org/10.1016/j.ijsolstr.2010.03.002
  • Jun, T.-S., Venter, A.M., Korsunsky, A.M., (2011). Inverse eigenstrain analysis of the effect of non-uniform sample shape on the residual stress due to shot peening, Experimental Mechanics 51(2):165-174. DOI: 10.1007/s11340-010-9346-2.
    » https://doi.org/10.1007/s11340-010-9346-2
  • Korsunsky, A.M., (2005). On the modeling of residual stresses due to surface peening using eigenstrain distributions, Journal of Strain Analysis for Engineering Design 40(8):817-824. DOI: 10.1243/030932405X30984.
    » https://doi.org/10.1243/030932405X30984
  • Korsunsky, A.M., Regino, G.M., Nowell, D., (2007). Variational eigenstrain analysis of residual stresses in a welded plate, International Journal of Solids and Structures 44(13):4574-4591. DOI: 10.1016/j.ijsolstr.2006.11.037.
    » https://doi.org/10.1016/j.ijsolstr.2006.11.037
  • Korsunsky, A.M., (2009). Eigenstrain analysis of residual strains and stresses, Journal of Strain Analysis for Engineering Design 44(1):29-43. DOI: 10.1243/03093247JSA423.
    » https://doi.org/10.1243/03093247JSA423
  • Krawitz, A.D., (2011). Neutron strain measurement, Materials Science and Technology 27(3):589-603. DOI: 10.1179/1743284710Y.0000000029.
    » https://doi.org/10.1179/1743284710Y.0000000029
  • Levers, A., Prior, A., (1998). Finite element analysis of shot peening, Journal of Materials Processing Technology 80-81:304-308. DOI: 10.1016/S0924-0136(98)00188-5.
    » https://doi.org/10.1016/S0924-0136(98)00188-5
  • Meguid, S.A., (1991). Effect of partial coverage upon the fatigue fracture behavior of peened components, Fatigue and Fracture of Engineering Materials and Structures 14(5):515-30. DOI: 10.1111/j.1460-2695.1991.tb00680.x.
    » https://doi.org/10.1111/j.1460-2695.1991.tb00680.x
  • Menig, R., Pintschovius, L., Schulzea, V., Vohringer, O., (2001). Depth profiles of macro residual stresses in thin shot peened steel plates determined by X-ray and neutron diffraction, Scripta Materialia 45(8):977-983. DOI: 10.1016/S1359-6462(01)01063-6.
    » https://doi.org/10.1016/S1359-6462(01)01063-6
  • Miao, H.Y., Larose, S., Perron, C., Lévesque, M., (2010). An analytical approach to relate shot peening parameters to Almen intensity, Surface and Coatings Technology 205(7):2055-2066. DOI: 10.1016/j.surfcoat.2010.08.105.
    » https://doi.org/10.1016/j.surfcoat.2010.08.105
  • Mura, T., (1987). Micromechanics of defects in solids, Kulwer Academic Publishers (Dordrecht).
  • Niku-Lari, A., Lu, J., Flavenot, J.F., (1985). Measurement of Residual-Stress Distribution by the Incremental Hole-Drilling Method, Journal of Mechanical Working Technology 11(2):167-188. DOI: 10.1016/0378-3804(85)90023-3.
    » https://doi.org/10.1016/0378-3804(85)90023-3
  • Shen, S., Atluri, S.N., (2006). An analytical model for shot-peening induced residual stresses, CMC: Computers, Materials, & Continua 4(2):75-86. DOI: 10.3970/cmc.2006.004.075.
    » https://doi.org/10.3970/cmc.2006.004.075
  • Timoshenko, S.P. and Goodier, J.N., (1970). Theory of elasticity, third ed. McGraw-Hill (New York).
  • Van Luchene, R.D., Johnson, J., Carpenter, R.G., (1995). Induced stress relationships for wing skin forming by shot peening, Journal of Materials Engineering and Performance 4(3):283-290. DOI: 10.1007/BF02649065.
    » https://doi.org/10.1007/BF02649065
  • Zimmermann, M., Klemenz, M., Schulze, V., Lohe, D., (2009). Numerical studies on the influence of thickness on the residual stress development during shot peening, High Performance Computing in Science and Engineering '08. Springer, Berlin, pp. 481-492. DOI: 10.1007/978-3-540-88303-6_34.
    » https://doi.org/10.1007/978-3-540-88303-6_34

Publication Dates

  • Publication in this collection
    Dec 2015

History

  • Received
    28 Dec 2014
  • Reviewed
    10 May 2015
  • Accepted
    27 May 2015
Individual owner www.lajss.org - São Paulo - SP - Brazil
E-mail: lajsssecretary@gmsie.usp.br