Video Lecture On Engineering Fracture Mechanics, Prof. K. Ramesh, IIT Madras 1
Video Lecture On Engineering Fracture Mechanics, Prof. K. Ramesh, IIT Madras 1
In the last class, we had developed the equations for strain energy stored in slender
members subjected to axial load, torsion, then bending. The expressions you had already
learnt in a course in strength of materials. Nevertheless, we looked at what is the nature of
those expressions, and in fact towards the end of the chapter on energy release rate, we
would use them to find out the energy release rate.
initial strain energy, so the change is what you find here. And in this case the strain energy
decreases. (Refer Slide Time: 05:40)
And what is the challenging question, that we have now is, suppose I have a component
with a crack how to find out the strain energy in the presence of a crack? And what did
Griffith do? Griffith took the problem of a central crack following the example of a plate with
the circular hole, a plate with the elliptical hole; he considered the problem of an infinite plate
with the central crack subjected to uniaxial tension.
Once we take up crack-tip stress and displacement fields, we would develop displacement
as well as stress field then come back and do these calculations again, and satisfy ourself
that, we could find out the strain energy in the presence of a crack, from a mathematical
stand point accurately. For the current discussion, we will look at for a central crack, and
keep this information important. Crack with two tips this is also very subtle point, but equally
Video Lecture on Engineering Fracture Mechanics, Prof. K. Ramesh, IIT Madras 4
important. The strain energy is given as U suffix a equal to pi sigma squared a squared
divided by E.
This is for a central crack. Once you say central crack it is crack with two tips. And whatever
the result that you have seen here is for a infinite panel of unit thickness. So it is multiplied
by unity, and if you look at from dimensional point of view, this satisfies the expression for
energy. And the suffix a denotes strain energy in the presence of a crack.
Now, what we will do is, we will look at a relaxation analogy to satisfy ourself how this
expression could be viewed at. But ultimately, you will be able to convince yourself well, only
when you do the actual calculation, this we will postpone it for the moment. Nevertheless we
will have a look at it.
Then, what I do is introduce a crack at the center. Remind you whatever we discuss it is
meant for a brittle material, but for visualization purpose, it is convenient for you to take a
rubber sheet and then look at what will happen. Suppose, I introduce the crack what do you
anticipate? Up to some extent, the crack will remain stable. Beyond a length, the crack will
propagate and the strip will separate itself. This is you can visualize physically.
For the formation of crack, you need to have surfaces to be formed and energy is required to
form the surfaces; this energy would have come from the strain energy of the system. And
what is the kind of situation we are analyzing? We are looking at fixed grips, so there is no
external work done.
We will just look at what way the crack is put and for the formation of these two surfaces,
some energy would have come out of the strain energy of the stretched plate, and we will
have to calculate this quantum of energy. We will make a simple calculation and to make our
life simple, we consider a triangular area in the neighborhood of the crack is what is
releasing the energy to form the two new surfaces. It could be any shape. This is just to
illustrate a simple calculation methodology to find out strain energy in the presence of a
crack.
Video Lecture on Engineering Fracture Mechanics, Prof. K. Ramesh, IIT Madras 5
We can be satisfied only when you do the mathematical calculation completely on the
displacements of the crack faces, as well as the stress field associated with it, that would
require some more time for us to do that. Right now, we only look at an analogy and we also
bring in one more aspect; we have taken a thin plate. We are really considering a plane
stress situation and for plane stress situation, lambda is taken to be pi by 2.
This also will get clarified, when we actually do the calculations based on the displacements
and the stress field. So, when I do this, I get this as Ua equal to pi a squared B sigma
squared divided by E. The earlier expression, we saw was for a plate of unit thickness, now
you have a thickness of B and what do you see here? The strain energy in the presence of a
crack is related to a squared. So, when the crack length changes, the energy would change
as a parabola, second degree curve.
And now we will look at for the formation of two new surfaces what way we require the
energy? If we look at these two natures we can understand, there could be a crack growth
phase and a catastrophic crack
propagation. We would try to plot
them.
is 2 into 2a into B into gamma s. And this we denote as U suffix s. And if you look at this final
expression, it turns out to be 4a into B gamma s. It is a linear function of the crack length. As
the crack length changes this would vary linearly. And we would try to plot this variation in a
graph, where the X axis is taken as the length of the crack, and Y axis is taken as the
energy. This would be obviously a straight line. We will plot that.
So, I have this as a function of crack length you find that this as a straight line, this denotes
the surface energy. We have already got the result of what is the strain energy in the
presence of a crack, which is a second degree curve. You have already got that as U a
equal to pi a squared B sigma squared divided by E. This is plotted in this fashion.
Please note that, this is only a schematic its not drawn with actual values. It is only a
schematic to illustrate the nature of the surface energy, the nature of the strain energy.
Suppose, I find out what is the total, that would change in this fashion. That would change
like this. The slope would become 0 at a particular point, that demarcates the region from
stable crack growth to fracture. You get the critical crack length from that graph. This will
denote the critical crack length, so you have a region where it is incremental crack growth by
crack growth mechanism, and in this zone it is catastrophic crack growth by fracture
mechanism.
This is given pictorially, you know we would like to have this quantitatively. So we will look at
it from a mathematical stand point, but I would like you to make a neat sketch of this
diagram, this is very illustrative of what is really happening. So you are able to find out a
critical crack length beyond which the crack would propagate faster. And what we are really
looking at is the incremental change in crack length to occur. So we are really looking at that
the incremental energy requirement should be satisfied. Mathematically, how will you write
it? The change of strain energy in the presence of a crack should be equal to change of
surface energy with respect to the crack.
If, you differentiate it with respect to a, you will get this as pi sigma squared 2a divided by E
equal to 4 gamma s. If you differentiate with respect to 2a, you will get this as pi sigma
squared a divided by E equal to 2 gamma s. In fact, if you simplify this, you will come to this
Video Lecture on Engineering Fracture Mechanics, Prof. K. Ramesh, IIT Madras 7
step. But dont differentiate with respect to a because I have already emphasized he
considered the problem of a central crack. The moment you look at a central crack you
should recognize, that it has two crack-tips. If you have to interpret these equations properly,
you have to differentiate with respect to 2a. So the caution is dont do it like this, do it like
this. Even in the definition of energy release rate, we will bring in the energy per crack-tip.
That is very important.
So, we get the final expression as pi sigma squared a divided by E equal to 2 gamma s. This
could be rewritten as sigma root a equal to 2E gamma s divided by pi. A very very important
expression. And symbolic step by Griffith in advancing fracture mechanics, because he has
been able to find out the fracture strength in the presence of a crack. You could also look at
this expression differently as sigma equal to root of 2E gamma s by pi a.
I have taken out a in this left hand side, so it appears as sigma root a. And you know once
you develop an expression like this it has to be experimentally justified. Otherwise you know
you will not be satisfied with the mathematical development. And what is the kind of
experiment with Griffith performed to establish this? Griffith carried out a series of
experiments on glass tubes and spherical vessels subjected to internal pressure, p.
some value in the range of 0.25 to 0.28. And this he has got the left hand side expression.
Sigma root a has been able to get it. For you to get the right hand side expression you
need to get the value of surface energy; surface energy in glass has to be estimated. You
know you will have to congratulate Griffith for having taken the pains to perform such difficult
experimentation.
Is it same as what he had got in his experiments? It is not same. When students perform
experiments these days and compare it with theory for fear of experiment being rejected,
they report only a difference in the second decimal place. I have been mentioning, when you
are performing an experiment you have to learn to report the results that you have got
faithfully. Experimentalists have to be extremely honest, because the values may carry more
information than what you had
originally thought of.
This gives you a confidence the way you have developed the expressions. Sigma root a
equal to 2 gamma s E divided by pi is a meaningful step. But we will also look at in what
context Griffith had really approached this problem. He had not in his mind that there is
fracture mechanics that he has to develop. He had not approached the problem from that
perspective. His perspective was different. You know scientists were really looking at how to
estimate strength of the material from a theoretical stand point.
Video Lecture on Engineering Fracture Mechanics, Prof. K. Ramesh, IIT Madras 9
You know, I am going to present this without much of mathematics we will only look at the
final results. In fact, people who work on cohesive zone modeling really play with different
types of law, and then also incorporate it in the finite element modeling and analyze, how the
plastic zone develops near the vicinity of the crack. They do such kind of activity but we will
only look at salient features of it.
Our interest is to see, what kind of result people have got based on lattice property
calculation. So, if I want to have the force-separation law I need to consider the law such
that, the initial slope corresponds to the elastic modulus E of the material, whatever the
material that you are looking forward to. And the second aspect is a total work of separation,
that is area under the curve corresponds to the surface energy gamma s.
So, what you will have to look at is, people who are making theoretical calculation have
already brought in, the role of surface energy in some other form. So, this kind of thinking
was there at that time. So you need to have the force-separation law to have an initial slope
that corresponds to the elastic modulus E.
Total work of separation that is the area under the curve corresponds to the surface energy
gamma s. And finally, you should have the amplitude a maximum value that represents the
inter-atomic cohesive force. So,
you can obtain a force-separation
law if it exhibits the following three
properties.
We will have to see, whether this is so or not. People have also found out what is the value
of gamma s. For many materials it is found to be approximately Eb divided by 40. So when
you substitute this into this expression, you get a value for theoretical strength. And what is
the theoretical strength? Should be like E by 6. In fact it is very high. E by 6 is not a small
value. We will look at for mild steel, because mild steel everybody knows. Although we may
not apply fracture mechanics to mild steel, to dispute the theoretical strength calculations
are way off when you compare with experimental results, mild steel is a good example.
In the case of mild steel, you have Young's modulus as 200 GPa, so 200 by 6 would be
something like 33 GPa. Mild steel is no where near this. If you look at, its yield strength is
something like 220 MPa or so, its ultimate tensile strength would be something like 320 or
400 depends on the material composition. that is all in terms of MegaPascal whereas,
theoretical strength is in terms of GigaPascal, that means the theoretical strength
calculations are way off. Why this is so?
See, there are two issues that we will have to look at. In the case of conventional spring
design if you look at the design books, you will find the yield strength for the same material
when it is used as a wire is much higher than a bulk material, this is one aspect. There is
another aspect, which was noted by da Vinci. What he said was the length, if the length is
shortened, the strength increases. It is an experimental observation.
See we need an explanation for this. Remember experiment is truth. If the experiment is
carefully performed and the results are also honestly reported then, experiment is truth. You
have to find out reasons for such a behavior. And later you find there another person Le
Blanc in 1839, he also established long iron wires to be weaker than short wires of the same
diameter.
So, people are finding there is something to do with the size. And they generalized it as a
size effect, and they reported an increase in strength as the dimensions of the test piece
decreases. A too drastic conclusion. You know, this kind of a conclusion we cannot make it
now. In those days when they were struggling to find out, why the theoretical strength was
not achieved in actual practice. But they observed when the size of the specimen decreases
the strength increases. They have noted this. They needed an explanation for this. And this
was the aspect which Griffith was concerned about and he felt that, the size effect is actually
a crack size length effect.
So, what Griffith looked at was, in the case of glass fiber the presence of cracks have
diminished strength. So he was able to provide a rational explanation why, when the size of
the specimen diminishes, the strength increases. I would like you to make a neat sketch of
this graph. It is very illustrative very important observation by Griffith. And you should also
look at in those days when they were developing the theories if they had a doubt they had
always verified it by suitable experimentation. However, complex it may be.
Video Lecture on Engineering Fracture Mechanics, Prof. K. Ramesh, IIT Madras 12
One is he has related the fracture because of crack length. Other way of looking at it is, in
the case of Inglis solution even a very very small crack with very very small external load,
the stresses would be so high the specimen will break into pieces. But that is not so in actual
practice. So, thats a paradox.
The Inglis solution was useful, but utility of that solution is questioned, because people are
unable to explain how actual structures remain solid with cracks. So that paradox was
resolved by Griffith, that strength of the plate is independent of the size of the crack. That
was the immediate observation from Inglis solution. Whereas Griffith said the crack length
does matter.
I can find out using this expression fracture strength, as well as critical crack length for a
given stress value. So this expression is very important. You know, we have looked at from
very simplistic analysis that, you get energy from strain energy for the formation of two new
surfaces, and we were able to get the expression for fracture strength.
So, you can go and find out for a bulk glass what would be the size of the crack; the size of
the crack is of the order of 0.025 millimeter. Because a bulk glass has this crack length. Its
fracture strength is far below the theoretical strength. See, if you look at the graph we have
noted that theoretical strength is around 11000 MPa, whereas bulk glass strength is around
175 MPa. Because of the expression developed by Griffith, now, we are able to say that, the
crack length for the bulk glass would be around 0.025 millimeter. Now, you take a fiber of
glass thinner than this. What happens when it is much below the size of that crack? You
find the value approaches theoretical strength. Now, the story is complete. People were
struggling why, when the size of the specimen decreases, strength increases.
And people have also looked at another problem. You know when you have a three-
dimensional solid, where you have a penny shaped crack. The moment you have a penny
shaped crack, people found the expression only the numerical factor has changed. the form
remained same; you have Young's modulus, you have surface energy divided by crack that
remains same. only the numerical factor was changing, and this is a very important
observation utilized by Irwin to generalize crack problems.
So from the energy approach, Griffith identified formation of two new surfaces requires
energy to be given from the strain energy of the system helped him to get the expression for
fracture strength, which also explained the paradox of Inglis solution, also explained why
actual strength of the material is far below the theoretical strength, and later on when people
carried out for variety of practical situations, they found a similarity. So, it helped to develop
fracture mechanics in the years to come. that is a very important contribution by Griffith mind
you all his analysis were focused to ideally brittle solids. You should never forget that, thank
you.