Summary
Summary
The SIF was developed by Irwin to describe the stress state at a crack tip, is
related to the rate of crack growth, and is used as a failure criterion. Irwin defined
SIF as a near-crack-tip approximation to Westergaard's solution for the stress field
surrounding a crack [40, 41].
If the applied stress at the crack tip exceeds a critical value known as the fracture
toughness (KIC), crack propagation occurs. The fracture toughness represents the
material's resistance to crack propagation and is a material property that
characterizes its toughness. For a simple case, the failure in an infinite plate with
central crack can be expressed in terms of critical energy release rate in relation
with critical stress (σc) [39].
σc2 ∗𝜋𝑎
Gc = (2.1)
𝐸
Griffith's energy release rate (G) and Irwin's SIF (K) is related at the plane stress
K2
condition by the relationship G =
𝐸
K2 𝐸
And for the plane strain condition by G = where E′=
𝐸′ (1−ʋ2 )
From the fundamental definition, the potential energy P is related to the crack
𝑑𝑃
growth 𝑑𝑎 proposed by Irwin G = - [43]. From the previous equations, for
𝑑𝑎
different crack modes the following equations is valid for plane stress and plane
strain conditions [44].
1 (1+ʋ)
G = (KI2 + KII2) + KIII2 (2.3) for the plane stress.
𝐸 𝐸
(1−ʋ2 ) (1+ʋ)
G= (KI2 + KII2) + KIII2 (2.4) for the plane strain.
𝐸 𝐸
Paris and Erdogan studied the crack growth behavior experimentally to find a
relation between the crack growth rate 𝑑𝑎/𝑑𝑁 and the stress intensity factor range
𝛥𝐾 [45]. The fatigue crack growth rate in metals can be described by the following
empirical relationship:
da
= C(∆K)m (2.5)
dN
Where: da/dN is the crack growth rate, ∆K is the stress intensity range,
When: 𝑑𝑎/𝑑𝑁 is plotted against the 𝛥𝐾 in log-log scale it was found the two
extreme asymptotes indicate the start and the end of the crack life. The left value
occurs at 𝛥𝐾 = 𝛥𝐾𝑡ℎ signifies that the values of K below this value cause no crack
growth. The second value occurs at 𝐾𝑚𝑎𝑥 = 𝐾𝑐 signifies for the complete fracture as
shown in figure 2.2.
𝐝𝐚
𝐝𝐍
G
Gthresh Gpl Gc
1) They take the value of 1 at the nodes corresponding to the element and 0 at all
other nodes within the element.
While FEM is very effective for many problems, it faces challenges when
dealing with discontinuities in the solution, such as cracks as mentioned above.
Multiple strategies have been put forward for dealing crack-related problems.
These include methods like the quarter-point finite element method, the enriched
finite element method, the boundary collocation method, the integral equation
method, the body force method, the boundary elements method, the dislocation
method, as well as mesh-free techniques like the element-free Galerkin method. In
order to eliminate the need for re-meshing during crack modeling, various methods
have been suggested, including the incorporation of a discontinuous mode at the
element level, the utilization of a moving mesh technique, and an enrichment
approach based on the partition-of-unity X-FEM.
Where: the first term is relevant to all nodes within the model, the second term
applies to nodes enriched by the discontinuous jump function (crack interior), and
the third term associates to crack-tip.
H(x) is the discontinuous jump function. The Heaviside jump function exhibits a
discontinuous behavior along the crack surface, maintaining a constant value on
each side of the crack, with a value of +1 on one side and -1 on the other.
𝐹𝛼 (𝑥) is the asymptotic crack-tip functions that assigned to nodes of the element
where the crack tip locates.
𝜃 𝜃 𝜃 𝜃
𝐹𝛼 (𝑥) = (√𝑟 sin , √𝑟cos , √𝑟 sin θ sin , √𝑟 sin θ cos )
2 2 2 2
Where 𝑟and 𝜃are the polar coordinates found at the crack tip at 𝜃= 0 tangent to the
crack face.
STATUSXFEM:
The XFEM technology was introduced through the development of the level set
method (LSM), which led to the establishment of nodal signed distance functions
Ψ and Φ. These serve to identify the crack surface and the orthogonal plane, as
shown in Figure 3.2. Φ defines the geometric distance from a node to the crack
surface, while Ψ describes the plane that is perpendicular to the crack surface.
Nodes with Φ (Ψ) values approaching zero correspond to the crack surface (the
plane orthogonal to the crack surface).
XFEM in ABAQUS is used with the direct cyclic approach in combination with
Paris law to simulate crack propagation. The direct cyclic defines the cyclic
loading condition by defining of a periodic function, which will be used by the
software to apply the amount of load at its corresponding time. The crack growth is
defined using Paris law, which is based on relative fracture energy release rates.
The initiation of the crack propagation condition;
𝑁
f= >1
C1 ∆𝐺 C2
Where: ΔG is the relative fracture energy rate between its maximum and minimum
values.
Three common mixed mode model for evaluating the equivalent fracture energy
release rate GequivC : the BK law, the Power law and the Reeder law.
To define this model, provide GIC , GIIC , GIIIC , am, an and ao.
To define this model, some values must be provided GIC , GIIC and η. This model
provides a power law relationship combining energy release rates in Mode I, Mode
II, and Mode III into a single scalar fracture criterion.
GIII G +GIII
GIIC ) (
GII + GIII
) (G +IIG )η
I II + GIII
1.5 Contour integral method
The J-integral is a very important fracture parameter for linear and nonlinear
material response. It is related to the energy release associated with crack
extension. If the material behavior is linear, the J-integral can be linked to the
stress intensity factor. ABAQUS provides a procedure to get the J-integral based
on studies [147, 148], with the virtual crack extension/domain. The method is very
attractive due to its simplicity, its low cost for the analysis, and high accuracy,
even if a coarse mesh is used. The J integral is defined as the rate of energy release
per unit step for crack propagation. The equation representing the energy release
rate is given as follows:
𝐽 = lim ∮ 𝑛. 𝐻. 𝑞 𝑑г (3.2)
г⇾0
Where:
Г is a contour beginning from the bottom to the top of the crack surface, as shown
in figure 3.3a.
H is given by H= WI – σ * ∂u/∂x
Where W is the elastic strain energy if the material behavior is linear; and is
defined as the elastic strain energy density plus the plastic dissipation for elastic-
plastic or elastic-viscoplastic materials.
Following Shih et al. (1986), equation 3.2 can be written in the form
. .
𝜕𝑢
𝐽= ∫ 𝑚. 𝐻. 𝑞′ 𝑑г − ∫ t. . 𝑞′ 𝑑г (3.3)
𝐶+ 𝐶++ г+ 𝐶− 𝐶+ + 𝐶− 𝜕𝑥
Where:
m is the outward normal to the domain enclosed by the closed contour, as shown in
Figure 3.3b.
m = -n on г
Using the divergence theorem, the closed contour integral can be converted into
the domain integral
. .
𝜕 ′
𝜕𝑢
𝐽 = ∫ ( ) . (𝐻. 𝑞 ) 𝑑г − ∫ t. . 𝑞′ 𝑑г
𝐴 𝜕𝑥 𝐶+ + 𝐶− 𝜕𝑥
(a) (b)
Fig3.3. (a) Contour for evaluation of the J-integral and (b) Closed contour C+ C++
г+ C- encloses a domain A that includes the crack-tip region as г=0
Once the J integral is calculated which is the energy release rate so the SIF can be
evaluated by:
𝐸
Where: E′=
(1−ʋ2 )