Embrittlement of an elasto-plastic medium by an inclusion

A mathematical model for the embrittlement of a long elastic-plastic crack by a relatively small, misfitting inclusion is presented. The model makes direct contact with the Dugdale–Bilby–Cottrell–Swinden model as a limiting case. The particular case of an oxide inclusion with a triangular cross-section at the tip of an intergranular crack in the Ni-based superalloy RR1000 at 650∘C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$650\,^{\circ }\hbox {C}$$\end{document} is considered. The positive misfit of the intrusion provides an additional tensile load on the crack tip and on the plastic zone, raising the local stress intensity factor kI\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$k_I$$\end{document} and the crack tip opening displacement Δu\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varDelta } u$$\end{document} above those when the inclusion is replaced by a dislocation-free zone of the same length. It is shown that for a given misfit strain and inclusion shape, the enhancement of kI\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$k_I$$\end{document} and Δu\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varDelta } u$$\end{document} is controlled by a dimensionless parameter ω=(σ/σ1)c/(2l)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega = (\sigma /\sigma _1)\sqrt{c/(2l)}$$\end{document} where σ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma $$\end{document} is the applied stress, σ1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma _1$$\end{document} is the yield stress, c is the crack length and l is the length of the inclusion. The anti-shielding effect of the intrusion is significant only when ω≲6\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega \lesssim 6$$\end{document}. As a result of the anti-shielding effect of the intrusion, the stress singularity at the crack tip always exceeds the compressive normal stress that exists within the thickest part of the intrusion when it is isolated. It is also shown that the gradient of the hydrostatic stress within the intrusion subjected to different applied stresses drives the oxygen diffusion and, hence, assists the oxidation at the grain boundary. The fracture toughness is considerably greater than that of a bulk sample of the oxide particle, which we attribute to the plastic zone.

by a dimensionless parameter ω = (σ/σ 1 ) √ c/(2l) where σ is the applied stress, σ 1 is the yield stress, c is the crack length and l is the length of the inclusion. The anti-shielding effect of the intrusion is significant only when ω 6. As a result of the anti-shielding effect of the intrusion, the stress singularity at the crack tip always exceeds the compressive normal stress that exists within the thickest part of the intrusion when it is isolated. It is also shown that the gradient of the hydrostatic stress within the intrusion subjected to different applied stresses drives the oxygen diffusion and, hence, assists the oxidation at the grain boundary. The fracture toughness is considerably greater than that of a bulk sample of the oxide particle, which we attribute to the plastic zone.

Introduction
When a brittle inclusion forms at a loaded crack tip it may enhance the growth of the crack by providing a path of relatively low fracture energy. An important example is the embrittlement of Ni-based superalloys due to environmental oxygen at temperatures around 700 • C in rotor discs of gas turbine engines. At low temperatures, the fatigue crack growth is transgranular and is caused by localised plastic deformation at the crack tip (Liu and Oshida 1986;Andrieu et al. 1992;Bache et al. 1999;Knowles and Hunt 2002;Li et al. 2015;Cruchley et al. 2015a, b;Németh et al. 2016;Hörnqvist et al. 2014). But at temperatures around 700 • C in an oxidising environment, the mode of fracture changes to intergranular. In this case oxygen enters the crack tip, diffuses along a γ /γ or γ /γ interface and reacts to produce an oxide particle called an intrusion (McMahon and Coffin 1970;Andrieu et al. 1992;Bache et al. 1999;Knowles and Hunt 2002;Li et al. 2015;Cruchley et al. 2015aCruchley et al. , b, 2016Németh et al. 2016Németh et al. , 2017Hörnqvist et al. 2014;Viskari et al. 2013;Pineau and Antolovich 2009). Crack growth is then enhanced through fracture of the oxide intrusion (McMahon and Coffin 1970;Andrieu et al. 1992;Cruchley et al. 2015a, b;Németh et al. 2016Németh et al. , 2017. The oxide materials of which the intrusion is composed are brittle in comparison to the γ -matrix, with fracture toughnesses typically between 0.4 and 2 MPa m 1/2 (Evans 1995).
The stress field of a loaded intergranular crack has been widely assumed to enhance diffusion of oxygen along the grain boundary. This mechanism of failure of Ni-based superalloys at high temperatures has been called stress assisted grain boundary oxidation (SAGBO). It is quite different from the mechanism of Krupp (2005), which assumes the ingress of atomic oxygen weakens bonds at the grain boundary, resulting in decohesion. Németh et al. (2017) studied the temperature dependence of the ductility of a Ni-based superalloy in air between room temperature and 1000 • C. They observed a pronounced minimum in the ductility at 750 • C resulting from oxidation along the γ grain boundaries and incoherent interfaces of γ particles. This experimental study, together with the experimental work of Kitaguchi et al. (2013), provide compelling evidence that the embrittlement is a result of the formation of oxide particles at grain boundaries rather than grain boundary decohesion due to the ingress of atomic oxygen. In this paper we assume the SAGBO mechanism of failure applies to Ni-based superalloys.
As might be expected in a multi-component superalloy the oxide intrusion is a complex agglomerate of nano-scale crystals in a layered structure (Kitaguchi et al. 2013), reflecting the decreasing chemical potential of oxygen with increasing distance normal to the grain boundary. The middle of the intrusion comprises fast growing oxides such as NiO and CoO, while the more thermodynamically stable oxides such as Cr 2 O 3 and Al 2 O 3 , which can form at lower oxygen chemical potentials, form in the outer layers of the intrusion. The intrusion is wedge-shaped. It is some tens of nanometres thick adjacent to the crack tip, and it tapers along the grain boundary over a distance of a few microns (McMahon and Coffin 1970;Andrieu et al. 1992;Li et al. 2015;Németh et al. 2016Németh et al. , 2017Viskari et al. 2013;Kitaguchi et al. 2013Kitaguchi et al. , 2015. The length of the intrusion is much less than the length of the crack which is hundreds of microns (Andrieu et al. 1992;Viskari et al. 2013;Németh et al. 2017). In the matrix surrounding the intrusion and crack, slip bands have been observed (Németh et al. 2016;Hörnqvist et al. 2014), indicating that the failure is not purely brittle. Karabela et al. (2013) simulated using finite elements, fatigue crack growth and oxygen diffusion in the Ni-based superalloy RR1000. They treated the grain microstructure explicitly and adopted a power law relationship for the inelastic strain rate. The diffusion of oxygen was assumed to be primarily along grain boundaries and assisted by stresses calculated using a crystal plasticity finite element model. They optimised the parameters of the simulation using experimental data for RR1000. They then obtained a fatigue failure curve for crack growth depending on both the oxygen concentration and accumulated inelastic strain at the crack tip, which compared well with experimental results. In contrast, the approach we have taken in this paper is theoretical, providing physical insight into the influence of the oxide intrusion on the fracture of the elastic-plastic medium through explicit functional relationships. Evans et al. (2013) argued that diffusion of oxygen along the grain boundary is enhanced principally by creating a difference in the chemical potential of oxygen vacancies between the ends of the intrusion, and that the gradient of stress within the intrusion has a negligible effect. But to estimate the enhanced diffusion they assumed the variation of the normal stress within the intrusion ahead of the crack. In this paper we find the assumed variation in Evans et al. (2013) differs significantly from the profile we derive, most notably we find the intrusion adjacent to the crack tip is always loaded in tension for experimentally observed intrusion geometries. In addition, we find the gradient of the hydrostatic stress within the intrusion is consistent with experimental observations (Kitaguchi et al. 2013) of the growth of oxide intrusions under different loading conditions. Chan (2014) argued that an oxide intrusion, with a larger volume per Ni atom than the Ni-matrix, induces a compressive stress in the matrix and thus screens the crack tip. We find the opposite. In addition, Chan's use of the superposition principle of linear elasticity does not satisfy the boundary condition that the crack faces are free of tractions. In this paper we use the superposition principle in a self-consistent approach that satisfies the boundary conditions on the stress field.
The outline of the paper is as follows. In the next section we describe the geometry of the mathematical model, which includes a long crack, a misfitting intrusion at the crack tip, and a plastic zone ahead of the intrusion tip. In Sect. 3 we calculate the stress field of the isolated misfitting intrusion with a triangular crosssection in an infinite elastic medium with no applied stress. The superposition principle is used in Sect. 4 to add the stress field of the intrusion to the applied normal stress on the long crack. This is done by replacing the long crack by a semi-infinite crack and setting up an integral equation for the Burgers vector density of dislocations representing the crack and the plastic zone. The condition for the normal stress to be bounded at either end of the plastic zone enables the size of the plastic zone to be determined. This condition leads to the solution of the integral equation for the Burgers vector density, from which the crack tip opening displacement, the stress field, the local stress intensity factor at the crack tip, and the associated J-integrals are determined. A dimensionless parameter, which we call ω, emerges from the analysis. It shows how the applied normal stress, the yield stress, the crack length, and the length of the intrusion conspire to control the influence of the intrusion on the embrittlement of the medium.

The geometry of the model
We consider an infinite, homogeneous, isotropic linear elastic medium in which there is a stationary crack of length c in the plane y = 0, between x = −c and x = 0, and which is infinitely extended along the z-axis. There is a uniform applied stress σ yy = σ at y = ± ∞. The stress intensity factor of the crack is K I = σ √ π c. Ahead of the crack there is a grain boundary occupying the half-plane y = 0, x ≥ 0. We assume there is an oxide intrusion at the crack tip with a triangular crosssection, and infinitely long along the z-axis, as shown in Fig. 1. The base of the intrusion is at the crack-tip and has width t, and its tip is at x = l ≥ 0. We assume a plastic zone between x = l and x = a > l, and infinitely extended along the z-axis, in which the stress Fig. 1 Geometry of the model comprising a long stationary crack along the negative x-axis, an intrusion with a triangular cross-section and a plastic zone. The model is subjected to an applied stress σ yy = σ at y = ± ∞, so the crack is loaded in mode I. The geometry is invariant normal to the page and therefore the elastic field is one of plane strain is constant and equal to the yield stress, σ 1 . In line with experimental observations we assume the crack is much longer than the intrusion, c l. This is a plane strain model because the strain tensor has non-zero components e x x , e yy , e xy only. As the length l of the intrusion tends to zero the plastic zone approaches the crack tip and the model becomes equivalent to the DBCS-model (Dugdale 1960;Bilby et al. 1963) with a long crack. As the angle at the tip of the intrusion tends to zero the intrusion vanishes and it is replaced by a dislocationfree zone (DFZ) occupying the region 0 ≤ x ≤ l. We call this limit the DFZ-model Ohr 1981, 1983;Ohr and Chang 1982;Majumdar and Burns 1983;Weertman et al. 1983). These models are instructive limiting cases for comparison to our model.

The intrusion
In this section we calculate the stress field in an infinite medium arising from the intrusion only with no applied stress.

The stress tensor
For simplicity we assume the elastic constants of the intrusion are those of the surrounding isotropic elastic medium. The intrusion is assumed to be homogeneous and to be characterised only by its volumetric misfit. Assuming the intrusion is composed of NiO the volumetric misfit is Φ = 1.65 (Evans 1995). That is, when the intrusion forms its stress-free volume is 1.65 times greater than the volume occupied by the same number of nickel atoms in the superalloy. Assuming the dilatation is isotropic the linear misfit strain is then 18. Let the interface between the intrusion and the matrix be labelled Γ . We associate with Γ a contin- The intrusion semi-angle θ is enlarged for clarity. The 18% misfit along y is represented by a continuous distribution of interfacial edge dislocations. It is assumed the same misfit along x and z cannot be sustained elastically because l t and the intrusion has infinite length along z, so that the misfit is assumed to be fully relaxed along x and z. Note the extra half planes of the dislocations are inside the intrusion because it has a larger volume when it is stress-free. The dislocations form interstitial loops, closed at z = ± ∞ uous distribution of dislocations, which generates the elastic field of the intrusion. In Fig. 2 the dimensions of the intrusion are those of a triangular prism of the Ni-matrix before it has been oxidised. The interface Γ comprises a vertical interface, Γ 0 , and two inclined interfaces Γ + and Γ − (see Fig. 2). The intrusion semiangle is θ = arctan [t/2l]. TEM observations show that θ 1 (Kitaguchi et al. 2013). Since t l we assume the 18% misfit along y is accommodated elastically. That is, the interface Γ is in a state of forced elastic coherence along the y-axis. This is modelled by a continuous distribution of edge dislocations with infinitesimal Burgers vectors along the negative y-axis, and with their lines parallel to z, consistent with the FS/RH convention. But since the interface along the z-axis is infinitely long we assume it is fully incoherent along z, that is there are no dislocations with Burgers vectors along z. Similarly, since l t we assume the interfaces cannot sustain coherency along the x-direction and therefore these interfaces do not have dislocations with Burgers vectors parallel to x. Figure 2 shows schematically the continuous distribution of interfacial edge dislocations with Burgers vectors along −y and infinitely long lines parallel to z forming loops closed at z = ± ∞. The elastic field of the intrusion is generated by this continuous distribution of dislocations.
To model the volumetric misfit of the intrusion we assign a positive sign to the line sense along Γ 0 and a negative sign to the line sense along Γ + and Γ − . Let γ be the position vector of a point on the interface Γ , and dγ be an infinitesimal element of Γ . Let b(γ )dγ be the magnitude of the Burgers vector of the continuously distributed dislocations between γ and γ + dγ . Then the Burgers vector density b(γ ) is as follows: The total Burgers vector of the continuously distributed dislocations on Γ vanishes: The stress tensor for the intrusion, σ I αβ (x), is then: where x = (x, y) and μ and ν are the shear modulus and Poisson's ratio, and the functions ζ αβ (x) are determined by the stress field of a positive edge dislocation: It is convenient to introduce a stress σ 0 = μ m /2π(1− ν). In RR1000 the shear modulus and Poisson's ratio are μ = 74 GPa and ν = 0.285 (Karabela et al. 2013) at 650 • C, for which σ 0 = 2.96 GPa. These elastic constants are used throughout the paper whenever stresses in RR1000 are calculated.

The hydrostatic stress along the x-axis
The hydrostatic stress is of interest because oxygen vacancies experience a force to drift from regions of hydrostatic tension to regions of hydrostatic compression. We obtain the following expression for the hydrostatic stress on the x-axis, σ I (x, 0) = 1 Fig. 3 The hydrostatic stress σ I (x, 0) (5) in GPa generated by the intrusion along the x-axis for RR1000, plotted as a function of x for four values of the intrusion semi-angle θ where H (x) denotes the Heaviside step function. We note that as θ → 0, σ I (x, 0) → 0, as it should because in this limit the intrusion vanishes. But when θ > 0, the hydrostatic stress is discontinuous at x = 0. The discontinuity does not depend on the intrusion angle throughout the range 0 < θ ≤ π/2: The jump in the hydrostatic stress at x = 0 is 15.9 GPa for RR1000. The hydrostatic stress along the x-axis is plotted in Fig. 3 for RR1000 as an illustration of (5). The hydrostatic stress is tensile on the x-axis at all intrusion angles just outside the intrusion and compressive just inside at x = 0. As expected the thickest part of the intrusion is always in compression. At the apex of the intrusion the hydrostatic stress is always tensile both inside and outside the intrusion, diverging logarithmically as the tip is approached from both sides along the x-axis. Fig. 4 The normal stress σ I yy (x, 0) (7) in GPa generated by the intrusion along the x-axis for RR1000, plotted as a function of x for four values of the intrusion semi-angle θ

The normal stress along the x-axis
The normal stress σ I yy (x, 0) is of interest because it will contribute to tractions on the crack faces in the next section. We obtain the following expression for it: This function is plotted in Fig. 4 for RR1000. The normal stress is predominantly compressive inside the intrusion, becoming more so as the intrusion angle increases. The normal stress becomes tensile inside the intrusion near the tip at x = l, and σ I yy (x, 0) diverges as − sin(2θ)(1 + 2 sin 2 θ) log |x − l| as x → l. This is a geometrical effect arising from the tapered shape of the intrusion. Outside the intrusion σ I yy (x, 0) is everywhere tensile and it becomes negligible at distances of more than two intrusion lengths from the intrusion.

The intrusion at the crack tip
In this section we use the stress field of the intrusion to calculate the stress field and related properties of the whole model portrayed in Fig. 1.

The approximation of a semi-infinite crack
The crack and the plastic zone are modelled as continuous distributions of climb edge dislocations in the usual way, e.g. Lardner (1974). These dislocations have infinitesimal Burgers vectors along the negative y-axis and lines parallel to the z-axis, consistent with the FS/RH convention. The total traction at any point on the faces of the crack must be zero. The sources of the traction are the stress fields of the dislocations comprising the crack, the applied load σ , and the stresses created by the intrusion and dislocations in the plastic zone. The x y-component of the total stress tensor on the crack faces is zero by virtue of the mirror symmetry in the plane y = 0. The yz-component of the total stress tensor on the crack faces is zero because the strain field is confined to the x − y plane. Therefore, only the yy component of the total stress tensor on the crack faces has to be considered. For the plastic zone we assume the yy-component of the total stress tensor is a constant, σ 1 , which we interpret as the yield stress. Let L c denote the set of points along the x-axis occupied by the crack, −c ≤ x ≤ 0, and the plastic zone In L c the Burgers vector density of the dislocations is b c (x). Then, at any point x in the crack or plastic zone we require the following singular integral equation to be satisfied: where A = μ/2π(1 − ν). The integral is a Cauchy principal value integral because a straight dislocation does not exert a stress upon itself. The term involving two Heaviside step functions is equal to σ 1 if x is in the plastic zone, otherwise it is zero. In this equation the addition of the intrusion stress tensor, σ I yy , to the applied load σ is an application of the superposition principle.
In (8) we see that the integral extends over the entire crack length, c, and the plastic zone. But provided c is much larger than the intrusion and the plastic zone it is otherwise arbitrary. This dependence on c complicates the analysis unnecessarily and it may be removed, as we now show.
We have already noted in Fig. 4 that the stress σ I yy (x, 0) generated by the intrusion is negligible at distances along the x-axis from the intrusion of more than twice its length. In small-scale yielding the field of the plastic zone is also short ranged compared to the crack length. Therefore, the intrusion and plastic zone modify the distribution of Burgers vectors in the crack only in the vicinity of its tip at x = 0. Conversely, the intrusion and plastic zone are in the near-tip stress field of the crack. Summarising, the Burgers vector density inside the crack is modified by the intrusion and plastic zone near the tip at x = 0 only, and the intrusion and plastic zone are located in the stress field σ yy (x, 0) of the crack only slightly ahead of this tip.
Since most of the crack surfaces are unaffected by the presence of the intrusion and plastic zone, we may view their presence as a perturbation localised near the crack tip. However, we will see below that they affect the stress intensity factor of the crack significantly. Let b e (x) be the Burgers vector density of an elastic crack of length 2c, occupying −2c ≤ x ≤ 0 and y = 0, in an infinite elastic medium. The crack is subjected to an applied normal load σ at y = ± ∞, and the intrusion and plastic zone are absent. Our crack of length c will be modelled as half of this crack of length 2c, occupying −c ≤ x ≤ 0. We view b e (x) as an unperturbed Burgers vector density of an elastic crack because it is not disturbed by the plastic zone or intrusion. It is well known, e.g. Lardner (1974), that b e (x) is as follows: This distribution ensures the crack faces are free of tractions along their entire length. Ahead of the crack, at x > 0, we obtain the stress σ e yy (x, 0): The same stress field σ e yy (x, 0) just ahead of the crack may be obtained by extending the crack to x = −∞, so that it occupies −∞ ≤ x ≤ 0, and by introducing the following Burgers vector density into the semi-infinite crack: where K I is the same stress intensity σ √ π c. By comparing with (9) it is seen that b e (x) = b (0) e (x)+O(x/c). Therefore it is a good approximation to replace the unperturbed Burgers vector density b e (x) with b (0) e (x) to model (a) the stress field just ahead of the crack tip at x = 0 and (b) the Burgers vector density just inside the crack near x = 0. However, the traction on the crack faces within the semi-infinite crack with the Burgers vector density b (0) e (x) is then zero, not −σ as it should be. It follows that if we replace b e (x) in the finite crack of length 2c in (8) e (x) in a semi-infinite crack we must remove the applied stress σ from the right hand side of the equation. Equation (8) then becomes: where L extends from x = −∞ to x = 0 and from x = l to x = a, and b(x) is the Burgers vector density of the dislocations in L. Thus, the removal of the dependence on the crack length c from the equilibrium condition in (8) has been achieved without compromising the accuracy of the model in the crucial region near the crack tip. Attention can now be focused on the Burgers vector density b(x) in this region, where the interaction with the intrusion and plastic zone dominate.

The solution of the integral equation
To solve the integral equation (12) we introduce the complex function g(z) =i √ z/[(z − l)(z − a)], where z = x + iy. Then, using the techniques described by Lardner (1974) the formal solution is as follows: This solution exists when the following condition is satisfied: This condition ensures the Burgers vector density is finite at the ends x = l and x = a of the plastic zone, as required by the finite value of σ 1 throughout the plastic zone. The condition in (14) provides a relationship for the size of the plastic zone. Let η = l/a. The plastic zone  (14) provides the following implicit relationship for η in terms of the dimensionless variable In this equation χ = 3θ/2 + π/4 + arctan [η tan θ ] /2, and E (m) and Π (n; φ|m) are respectively the complete elliptic integral of the second kind and the incomplete elliptic integral of the third kind, where m is the parameter, n is the characteristic, and the amplitude φ = arctan √ (i tan θ − 1) / (1 −i η tan θ ) . There is ambiguity in the definitions of the elliptic integrals in the literature. In the Appendix we state the definitions used in this paper for the incomplete and complete elliptic integral of the first, second and third kinds. Here and elsewhere in this paper we assume σ 0 /σ 1 = 3 as in RR1000 at 650 • C. The variable ω ≡ (σ/σ 1 ) √ c/(2l) is a ratio of stress intensity factors, combining the applied stress σ , the yield stress σ 1 , the crack length c and the length of the intrusion l into a single dimensionless variable. We will show that the conditions under which the intrusion embrittles the material are described by a limited range of values of ω.
The solution of (15) for η is plotted in Fig. 5 for four values of the intrusion semi-angle. When θ = 0 there is no intrusion. The presence of the intrusion increases the size of the plastic zone. As θ increases the size of the plastic zone increases for a given value of ω. We note that even when the applied stress is compressive (so that ω < 0) there is still a plastic zone ahead of the intrusion. All these features stem from the tensile normal stress ahead of the intrusion seen in Fig. 4, which increases with θ .
Using (15) the solution of (13) for the Burgers vector density b(x) becomes: In this expression we follow the convention (Abramowitz and Stegun 1972) that whenever the phase of the elliptic integral of the third kind is not specified-as in the first occurrence of Π on the right hand side of this equation-it is understood the phase is π/2 and the integral is then a complete elliptic integral of the third kind. Also,

The crack tip opening displacement
We may use Eq. (16) for the Burgers vector density to calculate the crack tip opening displacement as follows: In this expression Z (φ|m) is the Jacobi Zeta function, K (m) is the complete elliptic integral of the first kind, and F(φ|m) and E(φ|m) are the incomplete elliptic integrals of the first and second kind, respectively. This expression for Δu is plotted in Fig. 6, where it is normalised to the crack tip opening displacement Δu 0 = 2aσ 1 /(Aπ 2 ) in the DBCS-model. The length of the plastic zone is a c. When θ = 0, taking the limit as l → 0, we recover the DBCS-model and the crack tip opening displacement approaches Δu 0 . But in the presence of an intrusion (i.e. θ > 0), the crack tip opening displacement increases with θ at a given value of ω. This is consistent with Fig. 5 where the size of the plastic zone increases with θ for a given value of ω. The asymptotic behaviour of all the curves in Fig. 6 at large values of ω reflects the diminishing significance of the intrusion as the intrusion length l decreases relative to the crack length c, and as the applied stress σ increases relative to the yield stress σ 1 .

The normal stress along the x-axis
The stress field in the medium is generated by edge dislocations with density b(x) (16) in the crack and plastic zone, and by the intrusion. The total stress tensor, σ αβ (x), is given by the superposition principle of linear elasticity: The functions ζ αβ (x) are defined in (4). The integrals in (18) may be evaluated by complex variable techniques, and the resulting components of the total stress tensor are stated in the Supplementary Material. The normal stress σ yy (x, 0) is of interest because it governs the fracture of the grain boundary and, hence, of the intrusion. In the crack this normal stress is zero, in the plastic zone it is σ 1 , and at other positions along the x-axis it is as follows: In Fig. 7 we plot σ yy (x, 0) for an intrusion with a semi-angle θ = 5 • , for four values of ω from − 1 to 40. It is seen that the normal stress at the tip of the intrusion approaches σ 1 for all values of ω and this small value of θ . This is also obtained in (19): This equation shows that for experimentally observed intrusion geometries the normal stress at the intrusion tip is finite and equal to σ 1 . In contrast, at the crack tip there is always a tensile square root singularity. It is interesting to compare the solid curve in Fig.  7, showing σ yy (x, 0) when there is no applied stress (ω = 0), with the curve in Fig. 4 for θ = 5 • , showing the stress generated by the intrusion only. The presence of the crack tip singularity in Fig. 7, and its absence in Fig. 4, puts both ends of the intrusion into normal tension in Fig. 7, although at the thicker end of the intru- Fig. 7 The yy-component of the total stress tensor on the x-axis for RR1000 as a function of x for four values of ω from − 1 to 40. The intrusion semi-angle, θ, is constant and equal to 5 • . The normal stress σ yy (x, 0) = 0 in the crack, for which x ≤ 0. The intrusion lies between x = 0 and x = l. In the plastic zone σ yy (x, 0) = σ 1 and it is between x = l and x = a. Note that a varies with ω Fig. 8 The yy-component of the stress tensor plotted against x for RR1000 for four values of the intrusion angle θ, and ω = 4 sion (x = 0 + ) the tensile region is very narrow at the smallest values of ω. This feature persists even when the applied load is slightly compressive (ω = − 1).
It is seen that as ω increases the spatial extent of the compressive region inside the intrusion diminishes and eventually all the intrusion is in normal tension. This behaviour changes the gradient of σ yy (x, 0) inside the intrusion significantly-a point we return to in Sect. 5.
In Fig. 8 we show the normal stress σ yy (x, 0) when ω = 4 for four values of intrusion semi-angle θ . As the intrusion semi-angle decreases the normal stress σ yy (x, 0) inside the intrusion becomes increasingly tensile; when θ = 5 • there is only a very small part of the x-axis where the normal stress is compressive. Fig. 9 The hydrostatic stress on the x-axis for RR1000 as a function of x for four values of ω from − 1 to 40. The intrusion semi-angle, θ, is constant and equal to 5 • . The hydrostatic stress σ (x, 0) < 0 in the crack, for which x ≤ 0. The intrusion lies between x = 0 and x = l. In the plastic zone σ (x, 0) ≈ σ 1 and it is between x = l and x = a

The hydrostatic stress along the x-axis
We obtain the following expression for the hydrostatic stress on the x-axis, σ (x, 0) = 1 3 α σ αα (x, 0): In Fig. 9 we plot σ (x, 0) for an intrusion with a semiangle θ = 5 • , for four values of ω from − 1 to 40. It is seen that the hydrostatic stress is similar to the normal stress illustrated in Fig. 7. However, the crack faces are under compression and the compressive stress inside the intrusion is smaller than that in Fig. 7.

The local stress intensity factor
Let k I be the local stress intensity factor of the crack tip. It is defined by k I = lim x→0 + √ 2π x σ yy (x, 0). It differs from the applied stress intensity factor K I owing to the interaction with the intrusion and the plastic zone. Using (19) we obtain the following expression for k I : In Fig. 10 we plot k I /K I as a function of ω for four values of the intrusion semi-angle θ . Taking the limit l → 0 in (22) we obtain k I = 0 because the model Fig. 10 Plot of the local stress intensity factor, k I , normalised to K I , as a function of ω for four values of the intrusion semi-angle then becomes equivalent to the DBCS-model where the crack tip is completely screened by the plastic zone. This is also evident in Fig. 10 where k I /K I → 0 as ω → ∞.
It is seen in Fig. 10 that when θ > 0 the presence of the intrusion raises k I above the local stress intensity factor of the DFZ-model. This shows that the existence of the intrusion, as compared with its replacement by a DFZ of the same length, increases the intensity of the stress singularity at the crack tip. k I increases with the length l of the intrusion. The size of the plastic zone increases with increasing ω, as seen in Fig. 5, which increasingly screens and blunts the crack tip. A more ductile material, characterised by a smaller value of the yield stress σ 1 , or a larger applied stress σ , reduces the stress intensification at the crack tip caused by the intrusion. When θ = 30 • the larger plastic zone screens the crack more effectively and k I is lower than that for θ = 5 • and 15 • , as seen in Fig. 10. Eshelby (1951Eshelby ( , 1956 derived the force on an elastic singularity as an integral taken around a contour enclosing the singularity, a result that was discovered independently by Rice (1968) who called it the J-integral. We consider three contours, S 1 , S 2 and S 3 , illustrated in Fig.  11. The radius of S 1 is very large compared to the crack tip, intrusion and plastic zone which it encloses. S 2 is a contour around the crack tip only with an infinitesimal radius. S 3 encloses the plastic zone only. The non-zero component, J x , of the J-integral for each of these contours is then: Fig. 11 The contours S 1 , S 2 and S 3 comprise the singularities of the medium (i.e. the crack tip and the plastic zone) Fig. 12 The ratio of J x (S 2 ) + J x (S 3 ) to J x (S 1 ), given by (23) as a function of ω for four values of the intrusion semi-angle θ, and σ 0 = 2.96 GPa for RR1000. For ω 6 the presence of the intrusion raises the ratio significantly above one

J-integrals
When . This is illustrated in Fig. 12. In this figure it is seen that for experimentally observed intrusion geometries the force J x (S 2 )+ J x (S 3 ) does not change significantly. Also, when ω increases the force J x (S 2 )+ J x (S 3 ) tends to J x (S 1 ) and, hence, the presence of the intrusion becomes increasingly negligible. In fact, only when ω 6 the intrusion raises significantly the force on the crack tip and on the plastic zone.

Discussion
By comparing Figs. 4 and 7 it is seen that the normal stresses within the intrusion differ significantly, even when there is no applied stress, from a linear superposition of the stress fields of the isolated intrusion (7), and the unloaded crack. It highlights the point made in the Introduction, with reference to the work of Chan (2014), of the necessity of satisfying the boundary conditions of traction free surfaces in the crack and a constant normal stress in the plastic zone, when the normal stress within the intrusion is calculated. Figures 7 and 8 show the variation of the normal stress within the intrusion subjected to different values of ω and to different intrusion semi-angles θ , respectively. It is seen that for experimentally observed intrusion geometries there is always a positive stress concentration at the crack tip, which generates a tensile region inside the thickest part of the intrusion. This tensile stress decreases with distance from the crack tip and, when ω is sufficiently small, it becomes compressive, gradually becoming tensile again towards the intrusion tip. At the intrusion tip the normal stress is tensile and is approximately equal to σ 1 . This behaviour of the normal stress indicates that at moderate values of ω (ω < 4) the crack will penetrate the intrusion during crack growth. Figure 9 shows the variation of the hydrostatic stress within the intrusion subjected to different applied loads. It is seen that when there is no applied load (ω = 0) the gradient of the hydrostatic stress within the intrusion encourages the drift of oxygen vacancies from the intrusion tip towards the crack tip, and hence the diffusion of oxygen atoms in the opposite direction. The gradient is reduced to almost zero when ω = 4 and when ω = 40 its sign is reversed. This suggests that the intrusion will grow to a larger size when there is no applied load, which is supported by the experimental observations of Kitaguchi et al. (2013). In their work on RR1000, the intrusions at static crack tips were examined following two loading conditions each for 5 h at 650 • C in an oxidising environment. One of the loads was close to zero, and the other was slightly less than the threshold for crack-growth. The intrusion formed under almost zero load was significantly longer. In contrast, Evans et al. (2013) suggested that oxygen diffusion is dominated by the stress at the ends of the intrusion (i.e. x = 0 and x = l). But in Fig. 9 we see that if that were true the growth of the intrusion would be independent of the applied load. We conclude that the origin of the enhancement of oxygen diffusion in SAGBO is primarily the gradient of the hydrostatic stress within the intrusion, and not the hydrostatic stresses at the ends of the intrusion.
The local stress intensity factor k I , normalised to the applied stress intensity factor K I , is shown in Fig.  10 as a function of ω for four values of the intrusion semi-angle θ . Since an intrusion raises the local stress intensity factor above that of the DFZ-model, which corresponds to θ = 0, it follows that the intrusion always provides an additional tensile load on the crack tip. This additional tensile load is clearly visible in the region −2l x ≤ 0 of Fig. 4. The force acting to advance the crack tip, as determined by J x (S 2 ) in (23), is proportional to k 2 1 , and therefore it is also raised by the presence of the intrusion.
If we adopt the fracture criterion that k I has to exceed a critical value then as we increase the applied load from zero, fracture of the intrusion will occur so long as k I exceeds the critical value. As the crack grows from c to c + δc, it penetrates the intrusion by δc and the effective length l of the intrusion decreases to l − δc. At a given applied stress σ and yield stress σ 1 , the variable ω then increases because δω ≈ (ω/(2l))δc, where the approximation assumes c l. Thus if δc > 0 then δω > 0 for constant values of σ and σ 1 .
That k I decreases with increasing ω may be seen by writing k I in (22) as k I = σ 1 √ 2πlk I (ω). Then we obtain: where the partial derivatives are evaluated at constant semi-angle θ of the intrusion, constant applied load σ and constant yield stress σ 1 . In Fig. 13 we see that the term in the square brackets of (24) is negative for ω ≥ 0 for each value of θ . Since the pre-factor σ 1 √ 2πl/ω is always positive it follows that k I always decreases with increasing ω, that is with increasing crack length at constant applied load, yield stress and intrusion semiangle. The crack may be arrested inside the intrusion if k I falls below the critical value.
It follows that if the crack is arrested inside the intrusion further crack growth will resume only when Fig. 13 To show that the term in square brackets in (24) is always negative for ω ≥ 0 the intrusion has grown by further oxidation, thereby increasing l and raising k I above the critical value. Thus, we may envisage steady-state crack growth with the crack tip always inside the intrusion, where the rate of crack growth is equal to the rate of growth of the oxide intrusion ahead of the crack. Evans et al. (2013) stated that fracture of the intrusion in RR1000 occurs when k I reaches approximately 10 MPa m 1/2 , about an order of magnitude larger than the typical critical stress intensity factors for oxides of the metals present in the superalloy. They attributed this observation to compressive stresses inside the intrusion just ahead of the crack. However, we have seen in Fig.  4 that the intrusion loads the crack in tension. And, as seen in Fig. 7 this tensile loading creates a tensile stress singularity just ahead of the crack, even when there is no externally applied load. We suggest the increased fracture toughness of the oxide intrusion, as compared with bulk samples of oxide material, originates from the plastic zone, which exists even when there is no externally applied load.
The extent to which the intrusion embrittles the material, by raising the force on the crack tip and on the plastic zone, is governed by ω. As ω increases the plastic zone increases in size, screening and blunting the crack tip, which brings about the decrease of k I /K I with increasing ω in Fig. 10. The intrusion then has a smaller influence on the stresses at the crack tip. Conversely, as the intrusion grows through oxidation, so that ω decreases, it begins to increase the force on the crack tip and on the plastic zone. As seen in Fig. 12 the intrusion begins to have a marked influence on both the crack tip and the plastic zone when ω 6. This requirement leads to the following approximate relation for the length of the intrusion when it begins to have a signif-icant influence on the crack and plastic zone: (25)

Conclusions
A mathematical model has been developed to provide physical insight into the embrittlement of a long elastic-plastic crack by a relatively small, misfitting inclusion through functional relationships. In particular, we addressed the important example of the embrittlement of the Ni-based superalloys RR1000 at 650 • C caused by an oxide intrusion. The dimensionless parameter ω = (σ/σ 1 )/ √ c/(2l) emerges from the analysis, where σ is applied normal stress, σ 1 is the yield stress, c is the crack length, and l is the length of the intrusion. This variable is key to understanding the embrittlement caused by the intrusion. We find the intrusion has a marked influence only when ω 6. At larger values of ω the intrusion is too small to have a detrimental effect on the fracture toughness of the alloy.
The intrusion is always 'anti-shielding', adding to a tensile applied load on the crack tip. Even in the absence of an externally applied load the tensile load on the crack tip arising from the intrusion creates a square root singularity in the stress field at the crack. In turn this singularity exceeds the compressive normal stress that exists within the thickest part of the intrusion when it is isolated, making it tensile. As a result, the intrusion fails near the crack tip because it is the region where the tensile stress is highest. When the crack advances into the intrusion the local stress intensity factor k I on the crack tip decreases with increasing crack length, because the effective length l of the intrusion decreases raising ω and lowering k I . Thus, it is possible that this reduction arrests the crack advance inside the intrusion. In that case we may envisage a steady state where the rate of crack growth is equal to the rate of growth of the intrusion ahead of the crack, and the plastic zone is continuously regenerated ahead of the intrusion. As a result of the anti-shielding effect of the intrusion both the local stress intensity factor k I and the length of the plastic zone, are increased compared to the case where the intrusion is replaced by a dislocation free zone of the same length. Similarly, the crack tip opening displacement Δu is increased. Although the intrusion embrittles the medium, the fracture toughness is larger than that of the bulk oxide because the plastic zone ahead of the intrusion shields the crack tip.
The experimental observations of Kitaguchi et al. (2013) and Fig. 9 support the view that it is the gradient of the hydrostatic stress within the intrusion that is responsible for the enhancement of oxygen diffusion through the intrusion, and hence the increased rate of growth of the intrusion.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http:// creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Appendix: Definitions of the elliptic integrals
Let m be the parameter, n the characteristic, and φ the amplitude. The definitions used in this paper of the incomplete elliptic integrals of the first kind, F (φ|m), second kind, E (φ|m), and third kind, Π (n; φ|m), are as follows: These definitions correspond to those in Abramowitz and Stegun (1972), where properties of the elliptic integrals may be found. m and n are usually real numbers between 0 and 1. The amplitude φ can be a complex number. When φ = π/2, the elliptic integrals are called complete. The complete elliptic integrals of the first kind, K (m), second kind, E (m), and third kind, Π (n|m), are then respectively: