Chapter 9d Fracture
Chapter 9d Fracture
T.L. Anderson
CRC Press, Boca Raton, USA (1995).
Fracture Mechanics
C.T. Sun & Z.-H. Jin
Academic Press, Oxford (2012).
Theoretical fracture strength and cracks
Let us consider a perfect crystalline material loaded in tension. Failure by fracture can occur
if bonds are broken and fresh surfaces are created.
If two atomic planes are to be separated the force required initially increases to a maximum
(Fmax) and then decreases. The maximum stress corresponding to F max is the theoretical
strength t . This stress is given by:
By Energy By atomistic E → Young’s modulus of the crystal
E E → Surface energy
consideration
TFS t approach TFS a0 → Equilibrium distance between atomic
a0 For many metals ~ 2 centres
0.01Ea0
This implies the theoretical fracture strength is in the range of E/10 to E/6*.
The strength of real materials is of the order of E/100 to E/1000 (i.e. much lower in
magnitude). Tiny cracks are responsible for this (other weak regions in the crystal could also be responsible for this).
*For Al:
Fma
Applied Force (F) →
As we have seen crack is an amplifier of ‘far-field’ mean stress. The sharper the crack-tip,
the higher will be the stresses at the crack-tip. It is a region where atoms are ‘debonded’ and
an internal surface exists (this internal surface may be connected to the external surface).
Cracks can be sharp in brittle materials, while in ductile materials plastic deformation at the
crack-tip blunts the crack (leading to a lowered stress at the crack tip and further alteration
of nature of the stress distribution).
Even void or a through hole in the material can be considered a crack. Though often a crack is considered to
be a discontinuity in the material with a ‘sharp’ feature (i.e. the stress amplification factor is large).
A second phase (usually hard brittle phase) in a lens/needle like geometry can lead to stress
amplification and hence be considered a crack. Further, (in some cases) debonding at the
interface between the second phase and matrix can lead to the formation of an interface
cracks.
As the crack propagates fresh (internal) surface area is created. The fracture surface energy
required for this comes from the strain energy stored in the material (which could further
come from the work done by externally applied loads). In ductile materials energy is also
expended for plastic deformation at the crack tip.
A crack reduces the stiffness of the structure (though this may often be ignored).
Hard second phase in
Though often in figures the crack is shown to have a large
A crack in a material the material
lateral extent, it is usually assumed that the crack does not lead
to an appreciable decrease in the load bearing area [i.e. crack is
2a
a local stress amplifier, rather than a ‘global’ weakener by
decreasing the load bearing area].
Characterization of Cracks
Cracks can be characterized looking into the following aspects.
Its connection with the external free surface: (i) completely internal, (ii) internal cracks with
connections to the outer surfaces, (iii) Surface cracks.
Cracks with some contact with external surfaces are exposed to outer media and hence
may be prone to oxidation and corrosion (cracking). We will learn about stress corrosion cracking later.
Crack length (the deleterious effect of a crack further depends on the type of crack (i, ii or
iii as above).
Crack tip radius (the sharper the crack, the more deleterious it is). Crack tip radius is
dependent of the type of loading and the ductility of the material.
Crack orientation with respect to geometry and loading. We will see modes of loading in
this context soon.
~
2a a
Modes of deformation of a cracked body (modes of fracture)
How many ways are there to load a cracked body?
Three ideal cases of loading of a cracked body can be considered, which are called the
modes of deformation:
Mode I: Opening mode
Mode II: Sliding mode
Mode III: Tearing mode
In the general case (for a crack in an arbitrarily shaped body, under an arbitrary loading),
the mode is not pure (i.e. is mixed mode). The essential aspects of fracture can be understood by considering mode I.
Mode I
Modes of Deformation /
fracture of a cracked body
Mode II
Mode III
Important note: the loading specified and the geometry of the specimen illustrated for Mode II & III above do not give rise
to pure Mode II and II deformation (other constraints or body shapes are required).
Fracture: Important Points
One of the goals of fracture mechanics is to derive a material property (the fracture
toughness), which can characterize the mechanical behaviour of a material with flaws
(cracks) in it.
Fracture can broadly be classified into Brittle and Ductile fracture. This is usually done
using the macroscopic ductility observed and usually not taking into account the microscale
plasticity, which could be significant. A ductile material is one, which yields before fracture.
Further, one would like to avoid brittle fracture, wherein crack propagation leading to
failure occurs with very little absorption of energy (in brittle fracture the crack may grow
unstably, without much predictability).
Three factors have a profound influence on the nature of fracture:
(i) temperature, (ii) strain rate, (iii) the state of stress.
Materials which behave in a brittle fashion at low temperature may become ductile at high
temperatures. When strain rate is increased (by a few orders of magnitude) a ductile
material may start to behave in a brittle fashion.
Temperature
Ductile material : y < f Factors affecting
Strain rate
(the nature of) fracture
State of stress
Ductile Low Temperature
Fracture Promoted by High Strain rate
Brittle Triaxial state of State of stress
Funda Check Why do high strain rate, low temperature and triaxial state of stress promote
brittle fracture?
High strain rate (by not giving sufficient time) and low temperature essentially have a
similar effect of not allowing thermally activated motion of dislocations (i.e. ‘not helping’
plastic deformation by slip).
In specific cases some of the slip systems being active at high temperatures may become
inactive at low temperatures.
By triaxial state of stress (SoS) we mean tensile stresses of same sign along ‘y’ and ‘z’ also.
Triaxial SoS does not promote crack propagation, but suppresses plastic deformation (click
on link below to know more). Since plastic deformation is suppressed the crack tip remains
sharp, thus promoting brittle fracture. Click here to know more about which state of stress is good for plastic deformation
So for plastic deformation the following order is better: tri-axial < bi-axial < uni-axial.
s2 s2
s1 s1 s1
s3
* The Scanning Electron Microscope (SEM) with a large depth of field is an ideal tool to do fractography.
Classification of Fracture (based on various features)
0.66 nm
IGF
grain boundaries along
(ii) Appearance of Fracture surface, which crack can propagate
Grain-2
(iii) Strain to fracture,
(iv) Crack Path, etc. (As in the table below).
Presence of chemical species at the crack tip can lead to reduced fracture stress and
enhanced crack propagation.
Presence of brittle phase along the grain boundaries (Fe3C along GB in steel, glassy phase
at GB in Si3N4 ceramics) can lead to inter-granular crack propagation. This preferred ‘weak’
path along grain boundaries implies low energy expenditure during fracture (i.e. low
fracture toughness).
e
an
Cleavage plane Pl
ip
Sl
[1] C.E. Inglis, Stresses in a plate due to the presence of cracks and sharp corners, Trans. Inst. Naval Architechts 55 (1913) 219-230.
[2] A.A. Griffith, The phenomena of rupture and flow in solids, Philos. Trans. R. Soc. Lond. A221 (1920) 163-198. → Fat paper!
Crack growth and failure Brittle Materials
Initially we try to understand crack propagation$ in brittle materials (wherein the cracks are
sharp and there is very little crack-tip plasticity). The is the domain of Linear Elastic
Fracture Mechanics (LEFM).
For crack to propagate the necessary global criterion (due to Griffith) and the sufficient
local criterion (due to Inglis) have to be satisfied (as in figure below). Global vs. Local
The kind of loading/stresses also matters. Tensile stresses* tend to open up cracks, while
compressive stresses tend to close cracks.
Griffith
Energy based Global
Crack growth criteria
Stress based Local
Inglis
$ Note: the crack propagation we will study in this chapter will be quasi-static (i.e. elastic wave propagation due to crack growth is ignored)
* More on this later.
Stress based criterion for crack propagation (Inglis criterion)
In 1913 Inglis observed that the stress concentration around a hole (or a ‘notch’) depended
on the radius of curvature of the notch. I.e. the far field stress (0) is amplified near the
hole. [(max / 0) is the stress concentration factor ()].
hole crack
A ‘flattened’ (elliptical) hole can be thought of as a crack.
One way of understanding this formula is that if max For a circular hole
exceeds t (the theoretical fracture stress), then the c
σ max σ 0 1 2
material fails. c
Notes:
The units of Us is [J/m] (Joules per meter depth of the crack→ as
this is a through crack).
Though Us has a symbol of energy, it is actually a difference
between two energies
(i.e. two states of a body→ one with a crack and one without).
Half crack length ‘c’ appears in the formula.
E is assumed constant in the process (the apparent modulus will decrease
slightly).
a is the ‘far field’ stress (this may result from displacements
rather than from applied forces see note later).
The formula for Us can be appreciated by considering the energy released from a circular
region of diameter 2c as in the figure below. (The region is cylindrical in 3D).
The energy released is:
2
1
Elastic energy released from a circular region U scircular region
2 E
a
c 2 (1)
The computation of the actual energy released is more involved and is given by the formula
as noted before:
c 2 a2
Us
E Energy released from this circular region
is given by the formula (1) as above
Plane stress (not a true value, but to get a feel of the
condition predominant region involved).
U s 2 c a2
Henc
e c E
For a body in plane strain condition (i.e. ~ thick in the z-direction, into the plane of the
page), E is replaced with E/(12):
This is the difference in the energy between a body with a crack and one without a crack.
As pointed out before, the surface energy is the fracture surface energy and not just the
surface free energy. The origin of this energy is contributions from dissipative mechanisms
like plastic deformation, micro-cracking & phase transformation, in addition to the energy
of the ‘broken bonds’.
The units are Joules per meter depth of the body: [J/m].
Important note
The “Griffith experiment” is easily understood in displacement control mode (i.e. apply a
constant displacement and ‘see’ what happens to the crack) and is more difficult to
comprehend it in the force control mode (by applying constant ‘far-field’ forces).
In force control mode, the forces do work on the system and hence the ‘energy accounting’
process is more involved.
Hence, it is better to visualize as arising from ‘far field’ applied displacements.
Now we have the formulae for Us & U (which are required to write down the Griffith’s condition):
dU s dU c a2 U s 2 c a2 U
2 f c
E c
2 f c
dc dc E
LHS increases linearly with c, while RHS is constant.
The ‘equal to’ (=) represents the bare minimum requirement (i.e. the critical condition) →
the minimum crack size, which will propagate with a ‘balance’ in energy (i.e. between
elastic energy released due to crack extension and the penalty in terms of the fracture
surface energy). * 2 E f
c
The critical crack size (c*): a2(Note that ‘c’ is half the crack length internal)
A crack below this critical size will not propagate under a constant stress a.
Weather a crack of size greater than or equal to c* will propagate will depend on the Inglis
condition being satisfied at the crack-tip.
This stress a now becomes the fracture stress (f)→ cracks of length c* will grow
(unstably) if the stress exceeds f (= a)
2 E 2 E
At constant c (= c*) f
Griffith f when exceeds f then specimen fails c* (1 2 )
c*
Plane strain
An alternate way of understanding the Griffith’s criterion (energy based), though personally I
prefer the previous method
0 c 2 a2
Change in energy on the introduction of a crack U 4 f c
E
This change in energy (U) should be negative with an increase in crack
length (or at worst equal to zero). I.e. (dU/dc) ≤ 0.
At c* the slope of U vs c curve is zero [(dU/dc)c* = 0]. This is a point
of unstable equilibrium.
0
With increasing stress the value of c* decreases (as expected→ more
elastic strain energy stored in the material).
dU
0
Positive slope dc c* re ss
st
si ng
r ea
c0 Inc
U →
c1* c 2*
U →
c*
Negative slope c →
Stable
cracks Unstable cracks
c 2 a2
U 4 γ f c Us
c → Equations for ready reference
E
Griffith versus Inglis criteria
For very sharp cracks, the available elastic energy near the crack-tip, will determine if the
crack will grow.
On the other hand if available energy is sufficient, then the ‘sharpness’ of the crack-tip will
determine if the crack will grow.
A sharp crack is limited by availability of energy, while a blunt crack is limited by stress
concentration.
2 E E
f f
c* 4a 0 c
Griffith Inglis
8a0
If Griffith' s and Inglis criterion give the same result
If 3a 0 Griffith's and Inglis criterion give the same result
the 'Dieter' cross-over criterion
2 E E
c* c
*
2f 4a 2
0 f
‘Modern’ Fracture Mechanics
G.R. Irwin[1]
Stress Intensity Factor (K)
Material Parameter Fracture Toughness (KC)
Energy Release Rate (G)
Material Parameter Critical Energy Release Rate (GC)
J-integral
Material Parameter: JC
[1] G.R. Irwin, “Fracture Dynamics”, in: “Fracture of Metals”, ASM, Cleaveland, OH, 1948, pp.147-166.
[2] G.R. Irwin, “Analysis of stresses and strains near the end of a crack traversing a plate, J. Appl. Mech 24 (1957) 361-364.
Fracture Mechanics
Historically (in the ‘old times’ ~1910-20) fracture was studied using the Inglis and Griffith
criteria.
The birth of fracture mechanics (~1950+) led to the concepts of stress intensity factor (K)
and energy release rate (G). Due to Irwin and others.
Concept of Energy Release Rate (G)
G is defined as the total potential energy () decrease during unit crack extension (dc):
G GC For perfectly brittle solids: GC = 2f (i.e. this is equivalent to Griffith’s criterion).
* It is better to understand the basics of fracture with fixed boundary conditions (without any surface tractions).
Stress fields at crack tips
For a body subjected far field biaxial stress 0, with a double ended crack of length 2c, the
stress state is given by (this is mode-I loading):
KI 3
xx Cos 1 Sin Sin (1)
2r 2 2 2
KI 3
yy Cos 1 Sin Sin (2)
2r 2 2 2
KI 3
xy Cos Sin Sin (3)
2r 2 2 2
Fig.1
Note the inverse square root (of r) singularity at the crack tip. The intensity of the
singularity is captured by KI (the Stress Intensity Factor). I.e. K is the scaling factor for the singularity.
I
As no material can withstand infinite stresses (in ductile materials plasticity will intervene),
clearly the solutions are not valid exactly at (& ‘very near’) the crack tip.
At = 0 and r → the stresses (xx & yy) should tend to 0. This is not the case, as seen
from the equations ((1) & (2)). This implies that the equations should be used only close to
crack tip (with little errors) or additional terms must be used.
Understanding the stress field equation
KI 3 f ( )
xx Cos 1 Sin Sin → xx K I (1)
2r 2 2 2 2r
xx K I f ( r, )
K I Y 0 c Half the crack length
“KI (the Stress Intensity Factor) quantifies the magnitude of the effect of stress singularity
at the crack tip”[1].
Quadrupling the crack length is equivalent to doubling the stress ‘applied’. Hence, K
captures the combined effect of crack length and loading. The remaining part in equation(1)
is purely the location of a point in (r, ) coordinates (where the stress has to be computed).
Note that there is no crack tip radius () in the equation! The assumptions used in the
derivation of equations (1-3) are: = 0, infinite body, biaxial loading.
The factor ‘Y’ is considered in the next page.
KI > KIC
Irwin [K] Concept of stress intensity factor.
(in mode I)
Energy release rate based. Same as K based
- [G]
criterion for elastic bodies.
J-integral
Region of K-Dominance
The crack tip fields consists of two parts: (i) singular part (which blow up near the crack
tip) and (ii) the non-singular part.
The region near the crack tip, where the singular part can describe the stress fields is the K-
Dominance region. This is the region where the stress intensity factor can be used to
characterize the crack tip stress fields.
Fracture Toughness (Irwins’s K- Based)
One of the important goals of fracture mechanics is to derive a material parameter, which
characterizes cracks in a material. This will be akin to yield stress (y) in a uniaxial tension
test (i.e. y is the critical value of stress, which if exceeded ( y) then yielding occurs).
The criterion for fracture in mode-I can be written as:
K I K IC Where, KIC is the critical value of stress intensity factor (K) and is known
as Fracture Toughness
KIC is a material property (like yield stress) and can be determined for different materials
using standard testing methods. KIC is a microstructure sensitive property.
The focus here is the ‘local’ crack tip region and not ‘global’, as in the case of Griffith’s
approach.
All the restrictions/assumptions on K will apply to K IC: (i) material has a liner elastic
behaviour (i.e. no plastic deformation or other non-linear behaviour), (ii) inverse square
root singularity exists at crack tip (eq. (1)), (iii) the K-dominance region characterizes the
crack tip.
f ( )
xx K I (1)
2r
Fracture Toughness* (KIC) for some typical materials [1]
* We have already noted that fracture toughness is a microstructure sensitive property and hence to get ‘true’ value the
microstructure has to be specified.
** Note the strange units for fracture toughness!
[1] Fracture Mechanics, C.T. Sun & Z.-H. Jin, Academic Press, Oxford (2012).
Funda Check Is KIC really a material property like
y?
Ideally, we would like KIC (in mode-I loading) (KIIC & KIIIC will be the corresponding
material properties under other modes of loading $) to be a material property, independent of
the geometry of the specimen*. In reality, KIC depends on the specimen geometry and
loading conditions.
The value KIC is especially sensitive to the thickness of the specimen. A thick specimen
represents a state that is closer to plane strain condition, which tends to suppress plastic
deformation and hence promotes crack growth (i.e. the experimentally determined value of
KIC will be lower for a body in plane strain condition). On the other hand, if the specimen is
thin (small value ‘t’ in the figure), plastic deformation can take place and hence the
measured KIC will be higher (in this case if the extent of plastic deformation is large then K I
will no longer be a parameter which characterizes the crack tip accurately).
To use KIC as a design parameter, we have to use its ‘conservative
value’. Hence, a minimum thickness is prescribed in the standard
sample for the determination of fracture toughness.
This implies that KIC is the value determined from ‘plane strain tests’.
Essentially there are two approaches: global (energy based) and local (stress based).
For linear elastic materials the energy and stress field approaches can be considered
equivalent.
Ductile fracture
In ductile materials:
Crack-tip stresses lead to plastic deformation at the crack-tip, which further leads to
crack tip blunting.
Energy is consumed due to plastic deformation at the crack-tip (which comes from
elastic strain energy). This implies less energy is available for crack growth (& creation of new
surfaces).
Crack-tip blunting leads to a reduced stress amplification at the crack-tip. Blunting will
avoid ‘stress singularity’ at the crack tip and may lead to a maximum stress at a certain
distance from the crack-tip (as in the figure below).
Crack-tip blunting will lead to an increased resistance to crack propagation (i.e. increased
fracture toughness).
High magnitude of crack tip stresses can cause yielding at the crack tip (plastic
deformation).
This leads to crack tip blunting, which reduces the stress amplification.
There develops a zone ahead of the crack tip known as the process zone.
Funda Check
What else can happen at the crack tip due to high stresses?
i
Orowan’s modification to the Griffith’s equation to include “plastic energy”
2 2
Change in energy U 4 ( s p ) c c
E
2 ( s p ) E s ~ (1 2) J/m 2
f
c* p ~ (10 2 10 3 ) J/m 2
2 p E
f
c*
Ductile – brittle transition
Certain materials which are ductile at a given temperature (say room temperature), become
brittle at lower temperatures. The temperature at which this happens is terms as the Ductile
Brittle Transition Temperature (DBTT).
As obvious, DBT can cause problems in components, which operate in ambient and low
temperature conditions.
Typically the phenomena is reported in polycrystalline materials. Deformation should be
continuous across grain boundary in polycrystals for them to be ductile. This implies that
five independent slip systems should be operative (this is absent in HCP and ionic
materials).
This phenomenon (ductile to brittle transition) is not observed in FCC metals (they remain
ductile to low temperatures).
Common BCC metals become brittle at low temperatures (as noted before a decrease in
temperature can be visualized as an increase in strain rate, in terms of the effect on the
mechanical behaviour).
As we have noted before a ductile material is one which yields before fracture (i.e. its yield
strength is lower in magnitude than its fracture strength).
What causes the ductile to brittle transition phenomenon?
Both the fracture stress (f) and the yield stress (y) are temperature dependent. However,
as slip is a thermally activated process, the yield stress is a stronger function of temperature
as compared to the fracture stress.
If one looks at the Griffith’s criterion of fracture, f has a slight dependence on temperature
as E increases with decreasing the temperature ( also has a slight temperature dependence,
which is ignored here). y on the other hand has a steeper increase with decreasing
temperature.
f
Brittle Ductile 2 E
Griffith’s criterion f
T → c*
Inglis
E
f
4a 0 c
f , y →
f
y (BCC)
y (FCC)
T →
DBTT
No DBTT
Griffith versus Hall-Petch
Griffith Hall-Petch
2 E k
f y i
c* d
2 E 1 k'
f *
c c*
Grain size dependence of DBTT
T2 > T1
T1 f T2
T1
y
T2
f , y →
d-½ →
DBT
T2 > T1
T1 f
T1
y
T2
f , y →
Finer size
d-½ →
DBT
Torsion
Fatigue
Conditions of fracture
Creep
Temper embrittlement
Hydrogen embrittlement
Funda Check Why do we need a large ductility (say more than 10% tensile elongation)
material, while ‘never’ actually in service component is going to see/need such
large plastic deformation (without the component being classified as ‘failed’).
Let us take a gear wheel for an example. The matching tolerances between gears are so
small that this kind of plastic deformation is clearly not acceptable.
In the case of the case carburized gear wheel, the surface is made hard and the interior is
kept ductile (and tough).
The reason we need such high values of ductility is so that the crack tip gets blunted and the
crack tip stress values are reduced (thus avoiding crack propagation).
* 2 E
c
2
c →
Fracture
* stable
c
0
0 →
Rajesh Prasad’s Diagrams Validity domains for brittle fracture criteria
Sharp
cracks
>c
a0 3a0 →
Sharpest possible crack Approximate border for changeover of criterion
Safety regions applying Griffith’s criterion alone
c →
Unsafe
* 2 E
c
c* 2f
Safe
a0 →
c → Safety regions applying Inglis’s criterion alone
Safe E
Unsafe c
*
4a 2
0 f
a0 →
Griffith unsafe
Inglis unsafe Griffith
c → unsafe unsafe
Inglis safe
safe
c*
Griffith safe
Inglis
unsafe
Griffith safe unsafe
Inglis Griffith safe
unsafe Inglis safe
safe safe
a0 3a0 →
Role of Environment in Fracture
Similar to a critical value of the stress intensity factor (K IC) in normal fracture mechanics,
we can define a critical stress intensity factor in the presence of a corrosive environment (at
the crack tip) (KISCC). This value as seen from the table below can be much lower than K .
IC
* Metals are considered here, although other materials are also prone to such effects.
KIC SCC KISCC
Alloy
(MN/m3/2) environment (MN/m3/2)
13Cr steel 60 3% NaCl 12
http://en.wikipedia.org/wiki/Stress_corrosion_cracking
Another related phenomenon, which can be classified under the broad ambit of SCC is
hydrogen embrittlement.
Hydrogen may be introduced into the material during processing (welding, pickling,
electroplating, etc.) or in service (from nuclear reactors, corrosive environments, etc.).
Q&A What are the characteristics of brittle
fracture
Extreme case scenario is considered here:
Cracks are sharp & no crack tip blunting.
No energy spent in plastic deformation at the crack tip.
Fracture surfaces are flat.
Q&A What is the difference between plane stress and plane strain as far as fracture
goes?
C
END
Ductile fracture
Ductile fracture →
► Crack tip blunting by plastic deformation at tip
► Energy spent in plastic deformation at the crack tip
y y
Schematic
→
→
r → r →
Sharp crack Blunted crack