Quantitative Fracture Mechanics Approach to Fatigue Life Improvement by Shot Peening Application

Stress intensity factor has been a favorite tool in fracture mechanics due to its simplicity. Here, it is used to study the influence of the shot peening application in fatigue life improvement. Using a CT-specimen, the study of shot peening effect on fatigue life was performed. Here, the compressive residual stress is converted into residual stress intensity factor, Kres.. Together with the applied stress intensity factor, Kapl., the superposition principle can be used. The result of the investigation clarifies quantitatively what has been widely believed that one of the easiest ways to improve fatigue life is by perfecting the surface via shot peening. Compressive stresses are beneficial not only in postponing fatigue crack initiation but also in preventing the crack propagation, which lead to increasing fatigue strength.


Introduction
Shot peening has been believed to be beneficial's for quite some time. It has been used to improve fatigue resistance by introducing compressive stresses onto surface [1]. However, it has been treated too practical in industry without academic and scientific supports. The influence of shot peening to fatigue strength is generally underrated by engineering design codes. The difficulty of defining a general fatigue design approach is the reason, to correctly and quantitatively consider the effect of shot peening, in terms of residual stresses, surface

Cracked Body
To explore the characteristics of the stress field surrounding a crack in a loaded body, one can start from the Westergaard function [7]. The Westergaard function is a complex solution to the Airy stress functions. To do this, consider a coordinate system X, Y , Z in a stressed solid. At each point (x, y, z) one can define the stresses σ x , σ y , σ z , τ xy , τ yz , τ zx . Neglecting the body forces, for two dimensional problems, equilibrium equations are: If the displacements in x and y are u and v respectively, the expressions for the strains are: and the stress-strain relations: Where E is Young's modulus and v is Poisson's ratio. The above equations are satisfied if: The function ψ is called the Airy stress function. Substitution of the above equations and differentiating twice leads to the compatibility equation: One can define a complex function by: For Z to be an analytic function, the derivative dZ/dz must be defined unambiguously. This leads to the Cauchy-Riemann conditions: In the case of Mode I cracks it is convenient to use a function proposed by Westergaard. The Westergaard function is:

dZ Z dz
With Cauchy-Riemann equations, it follows that Which means that Equation (17) automatically satisfies the compatibility equation, which is Equation (12). By using equilibrium equation the stresses can be determined as: Any analytic function ( ) Z z will result in stresses defined by the equations above. For a crack problem of Figure 1, representing an infinite plate under biaxial stress, the stress function is: Where . z x iy = + The function is analytic except for -a ≤ x ≤ a, y = 0. The boundary stresses follow from Equations (19, 20 and 21). At infinity, where , z → ∞ the result is σ x = σ y = σ and τ xy = 0, and on the crack surface σ y = τ xy = 0, which means that the boundary conditions are satisfied.
It is more convenient to convert to a coordinate system with the origin at the crack tip, hence z should be replaced by (z + a). Turning to the general problem, which is shown in Figure 2, that has the form: Where r and θ are the polar coordinates of a point with respect to the crack tip. The parameter K I in these equations is known as the stress intensity factor. The stress intensity factor is then a measure for the stress singularity at the crack tip. Since the stresses are elastic they must be proportional to the external load. For the case of uniaxial tension with σ → ∞, it means that K I must be proportional to σ. In order to give the proper dimension to the stresses in Equation (25), K I must also be proportional to the square root of a length. For an infinite plate the only characteristic length is the crack size, hence K I must take form of: Comparison of Equation (24) and the above equation shows that The stress system parallel to the crack is not disturbed by the crack. Therefore, the solution for the uniaxial case must be the same as for the biaxial case.
However, for finite sized plates there are no closed form solutions available. For practical engineering use, the general equation for the stress intensity factor at the crack tip is Where β is a dimensionless constant, which depends upon geometry as The factor β has been calculated for many geometries and can often be obtained from various fracture mechanics literatures and standards.

Shot Peening
Shot peening has been used for decades as a measure to overcome surface engineering problems in various industries. Its improvements are produced mainly by combinations of compressive residual stress and cold work. Compressive residual stresses are known to be beneficial in increasing resistances to fatigue failures and corrosion fatigue, while the cold work effects of shot peening treatments can increase the surface hardness According to our boundary conditions, both σ y and τ xy are zero at the crack surface. The required real and constant value of f (z) at the crack tip is given by the notation K I , hence Taking polar coordinates from the origin with the stress at the crack tip can be calculated from compatibility equations to be: The general expression of these is usually also written as:  been known well that many parameters influence the efficiency of shot peening process. These are the peening coverage, saturation, shot material, shot size, speed, and peening time [8,10].

Analytical Solution and FEA Modeling for Shot Peening Application
While the whole complete report on this topics is already published previously by our research group [3,11], here the main result is represented. We developed a model by solving the problem for a single impact of a ball on a half-space. The FEA modeling was conducted using the nonlinear finite element (FE) code LS-DYNA 971. A model of the workpiece surfaces was developed using a constant stress sol-id element formulation which suitably models the impacted surface response of the workpiece. In a way it is similar to that of Purohit's [3,12] with mul-tiple shot balls.
The finite element discretization of the shot peening simulation is depicted in Figure 3. Shot peening loading was simulated by defining an im-pact velocity of 90 m/s to the impacting balls over a three repeated impact loading. Figure 4 shows the result of the FE models. One can see that as the number of the pass increases, the maximum com-pressive residual stress also increases.

Experiment
The shot peening application was conducted by shot peening on two different CT-specimens, namely sample SP-1, which has larger shot and sample SP-2, which has smaller shot size. Figure 5 shows the sample. The shot peening was only given at crack prospective line.
For the residual stress measurement, an X-ray diffraction (XRD) was used. The macroscopic resid-ual stress and information related to the degree of cold working can be obtained simultaneously by XRD methods. XRD is applicable to most poly-crystalline materials, metallic or ceramic, and is non-destructive at the sample surface. XRD meth-ods are well established, having been developed and standardized, E.g. [13,14].
Because the depth of penetration of the X-ray beam is extremely shallow, the diffracting volume can be considered to represent a free surface under plane stress. As shown in Figure 6, the biaxial sur-face stress field is defined by the principal stress-es, σ 1 and σ 2 , with stress normal to the surface = 0. of many materials [8]. It is believed to be the most economical and effective method of producing and making surface residual compressive stresses to in-crease the product life of treated metal parts. The increased strength of treated parts allows for light-er-weight parts that exhibit high wear and fatigue resistance. In some cases, shot peening induced grain refinement, martensitic phase transformation and compressive residual stresses and therefore, enhances mechanical properties in terms of hard-ness and work hardening. The induced compres-sive residual stresses are expected to play major role in reducing crack propagation rate, while the surface work hardening related to grain refinement and phase transformation can promote a higher re-sistance towards crack initiation enhancing crack growth threshold [9].
The process can be defined as work hardening to the surface of components by propelling streams of spherical shots to the surface. The surface layer of material yields plastically to generate residual compressive stress. Among the practitioners, it has Consider the strain vector, ε φψ , lying in the plane defined by the surface normal and the stress, σ φ , to be determined. ε φψ is at an angle ψ, to the surface normal, and can be expressed in terms of the stress of interest and the sum of the principal strains as, The only crystals which diffract X-rays are those The stress to be determined is the stress, σ φ , lying in the plane of the surface at an angle, φ, to the maximum principal stress, σ 1 . The direction of measurement is determined by the plane of diffraction. The stress in any direction (for any angle, φ) can be determined by rotating the specimen in the X-ray beam. If the stress is measured in at least three different directions, the principal stresses and their orientation can be calculated.    cept of the shot peening to academic world is not unique [10,12,[16][17][18]. Our approach is unique due to the work that we did was from very fundamental solid mechanics approach [3], continued to fracture mechanics. Figure 7 and Figure 8 show the results of actual residual stresses after shot peening with different shot sizes. It is easy to understand that the larger shot size produces deeper and larger compressive residual stress as we discussed in our previous publication [3]. Analytically, we proved that residual stress on the material is

On the residual stress creation
Where H ij σ -load stresses and sph ij σ -discharged stresses, obtained from the elastic solution. The first term of load stresses is indeed caused by the shot impact, which is needless to say influenced by the shot size. Figure 9 shows the result of the fatigue crack initiation and fatigue crack propagation test. In the figure, the time when the crack start to deviate from the base line is an indication of the crack initiation. The figure shows that the sample with no shot peening already cracked in less than 10,000 which are properly oriented relative to the incident and diffracted X-ray beam to satisfy Bragg's Law,

On the crack initiation
Where λ is the X-ray wavelength, n is an integer (typically 1), θ is the diffraction angle, and d is the lattice spacing. Subsequently, the measured strain is Where d 0 is the stress-free lattice spacing.
Using the concept highlighted above, the samples were then profiled with an X-ray diffraction device. Figure 7 and Figure 8 show the results. This trend is similar with both the analytical and finite element results as well as the finding of other researchers [15].
The CT-specimen were tested with the frequency of 40 Hz and the lengths were recorded electronically with occasional scope verification. Figure  9 shows the result. It is clear that both SP-1 and SP-2 have much longer time to break. The same data was than processed to construct Figure 10. In this figure, it is more obvious that SP-1 and SP-2 also have lower crack propagation rate.

Discussions
The attempts to bring the mostly industrial con-  fact that it lower the slope. In here, since mathematically, since the constant C in Paris law is simpler to be maintained, based on the following reason, the data is presented in term of m reduction.
Where m is directly related to the slope.

On quantitative fracture mechanics
For this, the concept of the weight function is used. The method is simple and can be applied to any shape of sample with any type of stress. In here, the application of the weight function in fracture mechanics is briefly described. Further discussion is available elsewhere [19]. Bueckner and Rice [19][20][21][22][23] introduced the concept of the weight function (WF), which enable one to calculate the stress intensity factor (SIF) for certain loading system by using a reference SIF for different loading system. A WF exists for any crack problem specified by the geometry of the component and a crack type. If this function is known, the SIF can be obtained by simply multiplying this function by the stress distribution and integrating it along the crack length. Basic Relations: A crack of length a in a body may be loaded by tractions T acting normal to a curve Γ, see Figure 12.
The tractions are responsible for a stress field cycles, while the sample shot peened with small shot size developed crack in approximately 15,000 cycle. This means that the shot peening delays the crack initiation significantly. The sample that was shot peened with larger shot size in this experiment was found to even postpone the crack initiation until 20,000 cycle. This is also reflected as increasing the threshold stress intensity factor as summarized Figure 11.
The shot peening has a very localized effect near the surface, which results in the introduction of compressive residual stresses, the propagation of cracks will be affected only in these areas. The shot peening increases the retardation of the surface crack propagation observed for small crack length that can be caused by residual stresses profile. This effect seems to increase near the threshold condition has also been verified by other researcher [6].

On the crack propagation
Fatigue crack propagation is typically the closest to be evaluated with stage II in the Paris curve. Figure 10 show the result of the fatigue crack propagation test. It is also summarized in Figure 11. Fundamental Paris law consists of C and m constants. Clearly the slope of the Paris law for shot peened samples are lower than that of non-shot peened sample, assuming the same value of m. Likewise, the value of the m reduces as the constant C is assumed constant. Since Paris law is purely derived based on empirical model, it does not have any significant metallurgical meaning in here besides the work directly like compressive applied stress. Consequently, the superposition principle also applies.
Where the subscripts use here are self-explanatory.

Conclusions
Fracture mechanics was used to quantitatively analyze the effect of shot peening on fatigue life.
The finding is that shot peening extend the fatigue life by both increase the stress intensity factor threshold and also lower the slope on stage II. In-creasing the threshold means delaying the crack ini-tiation while lowering the slope on stage II means lowering fatigue crack propagation until the crack becomes unstable. Within the limit of our evalua-tion, we also concluded that larger shot size is more effective than that of smaller shot size.
at the crack tip, which can be characterized by a SIF K T , where the superscript T refers to the loading system. As Bueckner and Rice [19][20][21][22][23] has suggest-ed one can write Where m is the vector of the WF, and WF is re-lated to the displacement field u, under an arbi-trary reference load by With H = E for plane stress and E/(1-ν 2 ) for plane strain. K ref is the SIF for the chosen reference load-ing case. In most practical cases of the mode I, the stresses along the crack prospective line are of in-terest, see Figure 13.
Therefore, referring to Figure 13, the distribu-tion of the stress perpendicular to the crack area in the uncracked component along the location of the crack, σ(x), the SIF for this stress distribution is given by Therefore, the WF depends only on the m, which is essentially independent of the stress state and depends only on the geometry. In our situation here, the value of the residual stress is constant.
This makes the relation to stress intensity factor simple. Therefore, it is not difficult to imagine that compressive residual stress in this case indeed