module no.# 04, lecture no. # 15: airy’s stress function...

14
Video Lecture on Engineering Fracture Mechanics, Prof. K. Ramesh, IIT Madras 1 Module No.# 04, Lecture No. # 15: Airy’s Stress Function for Mode - 1 Let us proceed towards developing crack-tips stress field equations in this class. I had already mentioned for solving crack problems, analytic functions are used as stress functions. You have read what is analytic functions in your earlier courses in mathematics. Nevertheless, we will look at some of the important aspects of it. (Refer Slide Time: 00:43) When you say a function is analytic, where you have a function w given as u of x, y defining the real part plus i into v of x, y. And z is your x plus i y. And when you say that function w is differentiable; we are defining this, dw by dz equal to limit delta z tends to 0. You have to evaluate the quotient, delta w by delta z. What we will have to look at is, for the derivatives to exist, no matter how delta z approaches 0, it is necessary that the limit of the quotient be the same. And we have looked at two possibilities. One is if delta z is real, delta y equal to 0. If delta z is imaginary, delta x equal to 0. And you have a graph that illustrates this. When delta z is imaginary, you will proceed like this. When delta z is real, you will proceed like this, or it could be any generic path, either as a straight line, or as a curve. In whichever way delta z tends to 0, I must have a unique value for delta w by delta z in the limit delta z tends to 0. (Refer Slide Time: 02:28) And what we had looked at in the last class, we have taken up two cases. One is when delta z is real. And we grouped the terms appropriately. And we found dw by dz equal to dou u by dou x plus i dou v by dou x.

Upload: doduong

Post on 07-Feb-2018

214 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: Module No.# 04, Lecture No. # 15: Airy’s Stress Function ...textofvideo.nptel.ac.in/112106065/lec15.pdf · Video Lecture on Engineering Fracture Mechanics, Prof. K. Ramesh, IIT

Video Lecture on Engineering Fracture Mechanics, Prof. K. Ramesh, IIT Madras 1

Module No.# 04, Lecture No. # 15: Airy’s Stress Function for Mode - 1

Let us proceed towards developing crack-tips stress field equations in this class. I had already mentioned for solving crack problems, analytic functions are used as stress functions. You have read what is analytic functions in your earlier courses in mathematics. Nevertheless, we will look at some of the important aspects of it.

(Refer Slide Time: 00:43)

When you say a function is analytic, where you have a function w given as u of x, y defining the real part plus i into v of x, y. And z is your x plus i y. And when you say that function w is differentiable; we are defining this, dw by dz equal to limit delta z tends to 0. You have to evaluate the quotient, delta w by delta z. What we will have to look at is, for the derivatives to exist, no matter how delta z approaches 0, it is

necessary that the limit of the quotient be the same.

And we have looked at two possibilities. One is if delta z is real, delta y equal to 0. If delta z is imaginary, delta x equal to 0. And you have a graph that illustrates this. When delta z is imaginary, you will proceed like this. When delta z is real, you will proceed like this, or it could be any generic path, either as a straight line, or as a curve. In whichever way delta z tends to 0, I must have a unique value for delta w by delta z in the limit delta z tends to 0.

(Refer Slide Time: 02:28)

And what we had looked at in the last class, we have taken up two cases. One is when delta z is real. And we grouped the terms appropriately. And we found dw by dz equal to dou u by dou x plus i dou v by dou x.

Page 2: Module No.# 04, Lecture No. # 15: Airy’s Stress Function ...textofvideo.nptel.ac.in/112106065/lec15.pdf · Video Lecture on Engineering Fracture Mechanics, Prof. K. Ramesh, IIT

Video Lecture on Engineering Fracture Mechanics, Prof. K. Ramesh, IIT Madras 2

(Refer Slide Time: 02:50)

Then we considered what will happen, when delta z is imaginary. Then again you can group the terms. And the final result turns out to be dw by dz equal to dou v by dou y minus i dou u by dou y. In fact in both the cases, we were only evaluating dw by dz.

Cauchy –Riemann Conditions

(Refer Slide Time: 03:28)

So, there has to be some kind of an interrelationship between the partial derivatives. And this is what is summarized in the next slide. These are famously known as Cauchy-Riemann conditions. If the derivative dw by dz, is to exist, the following conditions are to be satisfied. dou u by dou x, where u is the real part of the function w, v is the imaginary part of the function w. dou u by dou x should be equal to dou v by dou y; dou u by dou y should be equal to minus dou v by dou x. Only if these conditions are satisfied, then the function w is

differentiable. That’s very important.

And we have got these conditions by considering only two of the infinitely many ways in which delta z can approach 0. And what is the final conclusion is though I have looked at only two of the many possible ways, whatever the conditions that we have arrived at, are not only necessary but also sufficient conditions for the existence of the derivative of w. You know, I have always been mentioning, whenever you develop a condition mathematically, you must always qualify and see whether those conditions are necessary as well as sufficient. If they are not sufficient, what are the sufficient conditions? This is very important.

Page 3: Module No.# 04, Lecture No. # 15: Airy’s Stress Function ...textofvideo.nptel.ac.in/112106065/lec15.pdf · Video Lecture on Engineering Fracture Mechanics, Prof. K. Ramesh, IIT

Video Lecture on Engineering Fracture Mechanics, Prof. K. Ramesh, IIT Madras 3

(Refer Slide Time: 05:42)

And what you find is for a complex function, w to be differentiable, the Cauchy-Riemann conditions are given like this. If these two conditions are satisfied, the function w is differentiable and it’s an analytic function. And we also have certain generic definitions, we will look these aspects.

If w equal to f of z, possesses a derivative at z equal to z naught. And at every point in some

neighborhood of z naught, then function f is said to be analytic at z naught and z naught is called a regular point. You know, if you take any book in mathematics they would define some of these quantities. And if you take up a research paper, they might mention some of these aspects. So you should not feel uncomfortable, while reading research paper. So it is better that you get to know these terminologies. So you will have to know, when you have an analytic function what is a regular point?

If the function f is not analytic at z naught but if every neighbourhood of z naught contains points at which function f is analytic, then z naught is called a singular point of f of z.

So you define, what is a regular point? And you define, what is a singular point? And you also have one more definition. Suppose you find a function is analytic at every point of a region R. The function is analytic in R, or the function is regular and holomorphic. And this terminology will come across in research papers. They will say it’s a holomorphic function. So that means it is an analytic function where it has a differential available and you can differentiate it comfortably. And if you look at, we have to satisfy the bi-harmonic equation. Unless the function is differentiable, we would not be able to satisfy the bi-harmonic equation.

So we proceed towards solving problems using analytic functions. And we have to have a complete understanding on what aspects of analytic functions are important in this context. So you have to remember the Cauchy-Riemann condition, and also some of these definitions.

(Refer Slide Time: 08:30)

This is just to recapitulate, what kind of problems that we are going to solve in this chapter. We have already looked at, what is a mode I loading, what is a mode II loading, what is a mode III loading. And in each of these cases, our interest is

to find out the stress intensity factors KI KII and KIII. And from a material testing point of

Page 4: Module No.# 04, Lecture No. # 15: Airy’s Stress Function ...textofvideo.nptel.ac.in/112106065/lec15.pdf · Video Lecture on Engineering Fracture Mechanics, Prof. K. Ramesh, IIT

Video Lecture on Engineering Fracture Mechanics, Prof. K. Ramesh, IIT Madras 4 view, you would get KIc, KIIc, and KIIIc. As long as these stress intensity factors are below this, crack will not propagate.

(Refer Slide Time: 09:11)

Now our interest is to find out how to get the stress field for these modes of loading.

So we will look at what are Westergaard stress functions. See the moment you take up theory of elasticity, Airy's stress functions are very popular. We had seen quite a variety of problems that have been solved by that approach. And in 1939, Westergaard published a papers bearing pressures and

cracks in the journal of applied mechanics, where he proposed a simple stress function capital Z. And mind you this symbolism is slightly confusing.

Since these slides are prepared I am able to show the stress function as capital Z. This is a function of small z, which is equal to x plus iy. So when you write it using a software tool, you are able to show the capital Z and small z with sufficient difference. So while you take down the notes, you must find out a way of representing the capital Z and small z appropriately.

May be for the capital Z, you can put a horizontal line. So that will indicate not from the size of your letter but by that line, you can identify that you are referring to the function capital Z. And he has also written the function as real part of Z plus i imaginary part of Z. Though conventionally, we write the real part as u and imaginary part as v, once you come to solid mechanics, you know u and v are reserved as u displacement and v displacement. So, in order to avoid that, he has retained this as real part of Z, and imaginary part of Z. And you will have to live with this kind of notational scheme. There is no escape from it, so even the Cauchy-Riemann conditions; we would look it in terms of real parts of Z and imaginary part of Z. It will be very clumsy, while writing and also reading. But that’s the way, the paper has been written in those early days. It’s a very good paper. He had solved a variety of problems in just five pages. It’s lot of information available in that paper.

So what you have is, we now take up an analytic function capital Z, which is represented as real part of Z plus i imaginary part of Z. And based on Z, we would construct a relevant Airy's stress function for a given problem. And there are also other relationships that are needed for us to carry on further. We have to look at how to represent the derivatives and the integrals of Z. They are defined as follows.

So d of capital Z divided by d of small z equal to Z prime. Similar notation when you have a prime, you say, that this is a first derivative. When you have two primes, you have this as a second derivative. So you have d of capital Z prime divided by dz equal to Z of double prime. These are for differentials. You also have a notational scheme for representing the integrals, Z bar denotes integral capital Z dz.

See normally in complex numbers usual convention is when you put a bar, it’s a conjugate function. It is not referred in that context here. So you should look at the difference in symbolism and based on that interpret whatever the equation that you come across.

Page 5: Module No.# 04, Lecture No. # 15: Airy’s Stress Function ...textofvideo.nptel.ac.in/112106065/lec15.pdf · Video Lecture on Engineering Fracture Mechanics, Prof. K. Ramesh, IIT

Video Lecture on Engineering Fracture Mechanics, Prof. K. Ramesh, IIT Madras 5 So I have Z bar is defined as integral Z dz. Similar to your prime, you have a single prime and you have a double prime. You have a Z double bar that’s equal to integral Z bar dz. And you know in all our mathematical development, we would do it only as capital Z. We would not look at the details of the stress function right away. The advantage in the approach is, once we get the final set of expressions by changing capital Z, you could get solution for a variety of problems. So we will make a very generic mathematical calculations. We will not reduce it to a particular problem. We will develop it in a very generic sense.

(Refer Slide Time: 14:44)

And what are the Cauchy-Riemann conditions in terms of the function Z. We have already seen dou u by dou x equal to dou v by dou y. So which is written as dou of real part of Z divided by dou x equal to dou of imaginary part of Z divided by dou y. So, that is written straight away from your basic definition of Cauchy-Riemann conditions.

And you have to remember these conditions for simplification later.

And what is written on the left hand side? We have also written, what you get by differentiating the real part of Z with respect to x, gives you the real part of Z prime. In fact we would use this identity for simplification. On the other hand, if you look at the imaginary part when I differentiate it with respect to y, the left hand side, it is real part of Z prime. But you have imaginary part. When I do it, with differentiation with respect to x, the real part remains real. Whereas when you differentiate with respect to y, there is a flip over imaginary part turns to real. And the next equation, if you look at the real part will turn to imaginary. This you have to keep in mind. This is nothing but the same Cauchy-Riemann conditions you had written in terms of u and v. This is written in terms of the capital Z, the function itself.

So the next equation gives minus of dou real part of Z divided by dou y equal to dou imaginary part of Z divided by dou x. And this is equal to imaginary part of Z prime. You know you will realize its importance only when you want to satisfy the bi-harmonic equation for the given Airy's stress function. Then you will have to find out the derivatives and substitute. So all those, when you find the derivatives, you have to use these conditions. In fact, you will do that development as part of this class.

Page 6: Module No.# 04, Lecture No. # 15: Airy’s Stress Function ...textofvideo.nptel.ac.in/112106065/lec15.pdf · Video Lecture on Engineering Fracture Mechanics, Prof. K. Ramesh, IIT

Video Lecture on Engineering Fracture Mechanics, Prof. K. Ramesh, IIT Madras 6

Airy’s Stress Function for Mode-I

(Refer Slide Time: 17:18)

Make a neat sketch of this. And what I have here is, I have a central crack in an infinite plate subjected to what. There is a deviation. I have very clearly shown the infinite plate with a

central crack is subjected to biaxial loading. It’s very, very important. See if you look at what people had developed the solution for the plate of a small hole in a tension strip or an elliptical hole in a tension strip, they had only considered uniaxial loading. And even when the energy release rate was developed by Griffith, he considered a problem of a central crack subjected to uniaxial loading. But what you find here is Westergaard, has actually considered only the problem of a central crack in a biaxial loading.

Just because people were accustomed to solving these problems for uniaxial loading, rightly or wrongly whatever the solution which Westergaard provided they simply used it for a uniaxial situation. Then later researchers pointed out, no no this is not for uniaxial loading, it is for a biaxial loading. And if you look at the original paper, the diagram does not show a biaxial loading. The sentences show that he has done it for a biaxial loading. You know this is a misnomer and you have to be very careful about it.

So, you have to recognize that this is subjected to biaxial loading. And another aspect that you have to keep in mind, you know I have shown the central crack and the origin is fixed at the center of the crack, I have this as x axis, coinciding with the crack axis and y axis is perpendicular to that. And when I refer the point in the domain, you know we are not referring it in terms of Cartesian coordinates x comma y.

We are referring it in terms of r comma theta. why do I do? When I have complex numbers, it is easier to represent in terms of r e power i theta. So I will have only Cartesian stress components expressed in terms of r n theta. This is again something new. You know you have a function. The function is expressed in terms of r comma theta, but it is actually giving Cartesian stress components. And what is the stress function which was used, which is constructed based on the Westergaard stress function, was phi equal to real part of Z double bar plus y imaginary part of Z bar.

See in all of our development what we have looked at. If I have a problem once the stress function is specified everything about the problem is known. In fact in our earlier discussion on review of theory of elasticity, I just focused on certain salient aspects. I have not gone into the details of it.

Now, when we take the problem of a crack, we would investigate whether the function phi satisfy the bi-harmonic equation or not in all its completeness. We would see whether it satisfy the bi-harmonic equation. Then I have to do all those derivations. Then later on I

Page 7: Module No.# 04, Lecture No. # 15: Airy’s Stress Function ...textofvideo.nptel.ac.in/112106065/lec15.pdf · Video Lecture on Engineering Fracture Mechanics, Prof. K. Ramesh, IIT

Video Lecture on Engineering Fracture Mechanics, Prof. K. Ramesh, IIT Madras 7 would take up the stress function and ensure whether the boundary conditions of the problem are fully satisfied.

So, all that we will do systematically. So at this stage, what you will have to keep in mind is, I am having a central crack, where the center of the crack is taken as the origin. I had already pointed out, Irwin said don’t focus on the crack, focus at the crack-tip. He shifted the attention to the crack-tip.

So later, we would also shift the origin to the crack-tip and re-modify the equations. And that would really bring in, are we looking at a close formed solution or are we looking at a near field solution. So all those approximations you have to keep track off. And in the case of fracture mechanics, when some equation is specified, you should know whether they are referred with respect to the crack center as the origin or crack tip as the origin. So keep this at the back of your mind.

Satisfaction of B-harmonic equation

(Refer Slide Time: 22:44)

So the first test what we have to do is, when we have defined phi, we have to ensure whether it satisfies the bi-harmonic equation. And what is the bi-harmonic equation? I have this, as dou power 4 phi by dou x power 4 plus 2 times dou power 4 phi by dou x square dou y square plus dou power 4 phi divided by dou y power 4 equal to 0. First thing I have to do is, I have to get dou phi by dou x. Then go on with it until I find out dou power 4 phi by

dou x power 4.

And when I have to do this we have already defined phi. So what we are really looking at is, we have to get the differential of real part of Z double bar. That is dou by dou x of real part of Z double bar plus y into dou by dou x of imaginary part of Z bar. You know I would like you to take 2 minutes of your time and write the expression. If you write for one of this, rest all the derivation is simple and straight forward. Let me see how you take time to develop this. I will go round the class and see.

(No audio from 24:28 to 25:03)

I am happy to see some of you have got it correctly. This is nothing but real part of Z bar plus y imaginary part of Z. If you know how to write this, rest of the derivation is simple and straight forward. In fact you would do this for all the derivatives now. I would give you sufficient time now. We will have to get dou squared phi by dou x square.

So when I have dou squared phi by dou x square, I will have to differentiate this. And that is what is given here. So I have this as dou by dou x of real part of Z bar plus y into dou by dou x of imaginary part of Z. From whatever you have done for this case, what this amounts to, this is nothing but real part of Z plus y imaginary part of Z prime.

Page 8: Module No.# 04, Lecture No. # 15: Airy’s Stress Function ...textofvideo.nptel.ac.in/112106065/lec15.pdf · Video Lecture on Engineering Fracture Mechanics, Prof. K. Ramesh, IIT

Video Lecture on Engineering Fracture Mechanics, Prof. K. Ramesh, IIT Madras 8

(Refer Slide Time: 26:42)

So there is no difficulty in arriving at this. And we would proceed further. We will have to get dou cube phi by dou x cube. So I will have to differentiate this. And that’s what I am going to do. You know when you are differentiating with respect to x, it is simple and straight forward. It’s quite easy for you to handle.

So when I do this, I get this as real part of Z prime plus y imaginary

part of Z double prime. And mind you, these are all the function capital Z. I am not specifically saying that as capital Z. It is shown as capital Z in my slide. So when you are writing it put a horizontal bar to differentiate between capital Z and small z. In some way there should be a differentiation in your notes. In handwritten notes, you have to consciously do that.

If I have to get dou power 4 phi by dou x power 4, I have to do dou by dou x of real part of Z prime plus y into dou by dou x of imaginary part of Z double prime. This is again straight forward. You get this as real part of Z double prime plus y imaginary part of Z triple prime.

(Refer Slide Time: 28:55)

So now we move on to find the differential with respect to y. We want to get dou phi by dou y. You know you have to bring in all your understanding of your chain rule and all that. You know many times if you don’t use the mathematical approaches you tend to forget because I have the function as y imaginary part of something is coming. So when I have to do that, I have recognize that y is also a function. I have to use the chain

rule.

So I will have to have this as imaginary part of Z bar dou by dou y of y plus y into dou by dou y of imaginary part of Z bar. See I have nice slide. I have written this. Is it correct? Just look at, because you know while you are writing it, you tend to copy what is shown on the slides. You know you should reflect, whether what is written in the slide is correct or not. You know when something is written very neatly; even if it is wrong you tend to think that it is correct. What is the mistake here?

Page 9: Module No.# 04, Lecture No. # 15: Airy’s Stress Function ...textofvideo.nptel.ac.in/112106065/lec15.pdf · Video Lecture on Engineering Fracture Mechanics, Prof. K. Ramesh, IIT

Video Lecture on Engineering Fracture Mechanics, Prof. K. Ramesh, IIT Madras 9

(Refer Slide Time: 29:55)

So you have to really look at the Cauchy-Riemann conditions. I have written the first term. I have this as dou by dou y of real part of Z double bar. And if I really look at that there is a minus sign sitting here.

(Refer Slide Time: 30:10)

So I cannot forget that. So that is the problem that was there in that equation. So this is wrong. So

when I refer the Cauchy-Riemann conditions that becomes minus. And what is the advantage? This gets cancelled and your dou phi by dou y is simply y real part of Z, as simple as that. So you have to be careful. And also look at. When I am differentiating the real part, it switches to imaginary part, when I do it with y. When I differentiate with to respect to y, this also you should recognize. This directly comes from your Cauchy-Riemann conditions. See the Cauchy-Riemann conditions are very important. You

are just using them for finding out all this. There is nothing more or nothing less.

So I have y real part of Z. So you should use the chain rule while doing this. So that is what is depicted carefully in this slide. So I have this as real part of Z dou by dou y of y plus y into dou by dou y of real part of Z. Using Cauchy-Riemann conditions, you can write. This turns out to be real part of Z minus y imaginary part of Z prime. See right in the class if you derive them with your involvement rather than writing down from what is written in the slide, while you review the course, you also feel: Yes, I have derived all of them a sense of confidence will develop. These are not Greek and Latin; these were derived step by step. So whatever the final result, I get, I am confident of that. And you feel certain level of closeness to the solution. So try to involve yourself, solve it mentally and then agree with what is written in the slide. So do not take the slide as 100 percent correct. It is deliberately done so that you are alert in the class.

Page 10: Module No.# 04, Lecture No. # 15: Airy’s Stress Function ...textofvideo.nptel.ac.in/112106065/lec15.pdf · Video Lecture on Engineering Fracture Mechanics, Prof. K. Ramesh, IIT

Video Lecture on Engineering Fracture Mechanics, Prof. K. Ramesh, IIT Madras 10

(Refer Slide Time: 32:37)

So now we have to find out dou cube phi by dou y cube. So that is dou by dou y of real part of Z minus y into dou by dou y of imaginary part of Z prime plus imaginary part of Z prime dou by dou y of y. So when you simplify, you get an expression like this. I have this as minus imaginary part of Z prime minus of y real part of Z double prime plus imaginary part of Z prime. So these two add up. So I have this as minus 2 times

imaginary part of Z prime minus y real part of Z double prime.

See our ultimate objective is to see whether the bi-harmonic equation is satisfied. If the bi-harmonic equation is satisfied, what we would conclude? Phi is a valid stress function. Later on we will have to go and investigate what problem it represents. So we are testing out each step in the procedure. Now we evaluate dou power 4 phi by dou y power 4. So that is nothing but minus 2 times dou by dou y of imaginary part of Z prime minus of y dou by dou y of real part of Z double prime plus real part of Z double prime dou by dou y into y.

So upon simplification, you get this as minus 2 times real part of Z double prime or I could also put this. I could also have this as minus 3 times real part of Z double prime plus y imaginary part of Z three prime.

See we have evaluated the dou power 4 phi by dou x power 4 as well as dou power 4 phi divided by dou x power 4. Now we have to get dou power 4 phi dou x square dou y square. I can start from anywhere. I can start from x or I can start from y.

(Refer Slide Time: 35:24)

You know it’s easier because we have seen when you are differentiating with respect to y, you have to be careful. So I would directly start with dou squared phi by dou y square. Then I differentiate it with respect to x once, again differentiate it with respect to x again. That makes my life simpler. So that’s what is depicted here. So what I will do is, I would take the solution which we have already obtained for dou

squared phi by dou y square. Then writing dou square by dou x square is lot more simpler. It’s simple and straight forward. There are less chances of making any calculation error. And that is what is shown here.

Page 11: Module No.# 04, Lecture No. # 15: Airy’s Stress Function ...textofvideo.nptel.ac.in/112106065/lec15.pdf · Video Lecture on Engineering Fracture Mechanics, Prof. K. Ramesh, IIT

Video Lecture on Engineering Fracture Mechanics, Prof. K. Ramesh, IIT Madras 11 So I have to get dou squared by dou x squared of real part of Z minus y imaginary part of Z prime. And when you do this, I have to differentiate dou by dou x real part of Z prime minus y dou by dou x of imaginary part of Z double prime which finally turns out to be real part of Z double prime minus y imaginary part of Z triple prime.

Since you have all these expressions, now what we will have to do? We will have to go and substitute it in the bi-harmonic equation. And see what is it that we get? We have to get that should go to 0. Let us see whether it goes to 0.

(Refer Slide Time: 36:49)

So when I substitute all of these quantities, I get an expression like this. Please take some time to write this down, expression is sufficiently long. And it goes to 0.

So, what is the conclusion? The stress function selected satisfies the bi-harmonic equation. So that means our first step is correct. We are really looking at solving the problem, so we have to select a valid form of stress function. So the

first condition is the bi-harmonic equation should be satisfied. That is satisfied.

Boundary Conditions

(Refer Slide Time: 37:51)

Now we have to look at what are the boundary conditions. See boundary conditions are very, very important. Only to write the boundary condition, I have trained you on the concept of what is a free surface. And if you know how to write the boundary condition carefully and correctly the problem is solved. And in fact in this case, we would write the boundary condition; we would solve. It would appear as if we have satisfied all the boundary conditions

satisfactorily. Nevertheless, the solution will not reflect what is seen in an experiment. It is something very peculiar. In your strength of materials you have not done something like this. Because the theory is very well developed, so you don’t have to do this. But in the case of crack problems, you will have to investigate whether the boundary conditions have been fully satisfied or you have over satisfied some condition. All this have to be looked at. But first let us look at the simple problem. And you should know how to write the boundary condition.

And what I have here? I have the crack that is opening up. And how do you classify the crack surfaces? The crack surfaces are free. You should recognize that. You should

Page 12: Module No.# 04, Lecture No. # 15: Airy’s Stress Function ...textofvideo.nptel.ac.in/112106065/lec15.pdf · Video Lecture on Engineering Fracture Mechanics, Prof. K. Ramesh, IIT

Video Lecture on Engineering Fracture Mechanics, Prof. K. Ramesh, IIT Madras 12 recognize that it is a free surface. And you know we are also looking at small deformation. You have to look at the deformed configuration as well as undeformed configuration for writing out the boundary condition. But you will write the boundary condition only in the undeformed configuration because you say the deformations are small.

So we have to write what happens on the crack surface. Crack is defined, because we have taken the center of the crack as the origin. It is defined between minus a to plus a and y remains 0. So when I have this as a free surface, you have the axis of reference given, which are the stress components should remain 0. See we have sigma x sigma y and tau x y. Of this whether only one stress component should remain 0, or all the stress components should remain 0 or only some stress components should remain 0. This is where you apply your mind because somebody has given you the stress function. Once stress function is given, the mathematical procedure is straight forward that there is nothing special about it. Only thing is you have to remember the mathematical steps, and do it, where you have the creative potential is used, is only in looking at the boundary conditions. That is very, very important. I have this as a free surface.

Shear cannot cross a free surface, so shear stress can be 0. That one condition, I can write. Between sigma x and sigma y, which stress component should go to 0? Sigma y. So if you recognize that then the first condition is completely specified. You know this is, what I have mentioned, that I am showing the deformed configuration for clarity. And another thing we observe whenever I show the deformed configuration, the crack is opened up like an ellipse. There is a reason behind it. You know we have not really looked at the displacement field. The moment you look at the displacement field, the crack faces would open up like an ellipse. Ellipse is very closely linked to crack problem, in many ways. Only the problem of a elliptical hole, in a tension strip alerted the scientists, that cracks are very dangerous. It doesn’t stop there, the ellipse goes with the crack. It opens up like an ellipse only. We would see the equations and reconvince ourself. So whatever the animation you have is, capturing a realistic feature.

(Refer Slide Time: 37:51)

See what you should not confuse is, I have shown this as a curved surface. But I have written the equation for the boundary conditions as minus a; x between minus a plus a on y equal to 0. That means I am considering this still as a straight line for writing the boundary condition. This we are justified, because we are looking at small deformation.

And what are the other boundary conditions? We know what happens at the far field. At the far field, we have taken up a problem of biaxial loading; sigma x equal to sigma y equal to sigma. And no shear stress exists; tau xy equal to 0. This is where I alerted. See people have been used to looking at a small hole in a tension strip. Without looking into the way the solution was developed, if you only look at the final solution there is every reason you wrongly interpret for which problem the solution represents. See in the methodology what we are adapting, we would develop the first term in the series solution. That’s called as a singular solution.

That solution is valid, but that solution is not sufficient to explain experimentally observed fringe patterns. And we will have to investigate the boundary conditions for that. So I have

Page 13: Module No.# 04, Lecture No. # 15: Airy’s Stress Function ...textofvideo.nptel.ac.in/112106065/lec15.pdf · Video Lecture on Engineering Fracture Mechanics, Prof. K. Ramesh, IIT

Video Lecture on Engineering Fracture Mechanics, Prof. K. Ramesh, IIT Madras 13 written two boundary conditions. Are there any more boundary conditions? You have. Along the crack axis, you can write something about it. Along y equal to 0 for any x due to symmetry your shear stress has to be 0.

And you have to distinguish. This is a in-plane shear stress. You are talking about tau xy. See in this course, I said I am going to investigate whether the solution obtained from fracture mechanics is correct or not. Whether some modifications need to be looked at, I said I am going to use photoelasticity. In photoelasticity, you get contours of sigma 1 minus sigma 2. In fact they are contours of maximum shear stress. See the maximum shear stress is different from in-plane shear stress, remember that. And if you recall, I had mentioned along the crack axis, you see several fringe orders. This we had seen when the crack is very long or when the crack is situated in a stress concentration zone, along the crack axis you see a frontal loop and fringe order varies.

(Refer Slide Time: 37:51)

So that means along the crack axis, there is variation of maximum shear stress. And here you are making a statement, along the crack axis in-plane shear stress is 0. This is correct. But what we would show later is by imposing tau xy equal to 0, without our knowledge we have also imposed maximum shear stress is 0 in the process of the solution development. This was not looked at by people in initial stages. This was discussed by Sanford way

back in 1980s after very long time from Westergaard gave the solution. He gave the solution in 1939. But only around 1980s, Sanford identified that there is a deviation. In fact, if you look at the history of fracture mechanics development, in 1957, Williams came up with a different approach to finding the solution for the crack problem. He solved it from polar coordinates. And you have, what are known as William corner functions. That had infinite series solution.

Since the entire fracture mechanics community focused on the singular term, they have not looked at the advantages of the contribution by Williams. It was remained dormant. So the analytical people woke up only around 1977. They said higher order terms also play a role in fracture mechanics. And if you look at the history, when Westergaard proposed the solution in 1939, Wells performed an experiment around 1952. Around that time, for interpreting the results of photoelastic fringes, Irwin recognized the role of higher order terms. So experimentalists knew this 25 years earlier than analytical people. So we will have a look at all of that in the subsequent classes.

So what I want to say is in a normal problem situation, if you write the boundary condition 1, boundary condition 2 boundary condition 3; and if you satisfy the solution and obtain based on all this you would be happy that you have solved the problem. That is what people thought even in the initial stage of fracture mechanics. They thought that they have solved the problem. Only when they were unable to explain, what you find in experiment, they realized that such nuances have to be relooked at. You know it’s very important, you know.

And another aspect all along we have been saying, you know if you look at the boundary condition also. We are looking what happens close to the crack-tip, we are looking what happens at infinity. At this stage of solution development, we are really looking at the

Page 14: Module No.# 04, Lecture No. # 15: Airy’s Stress Function ...textofvideo.nptel.ac.in/112106065/lec15.pdf · Video Lecture on Engineering Fracture Mechanics, Prof. K. Ramesh, IIT

Video Lecture on Engineering Fracture Mechanics, Prof. K. Ramesh, IIT Madras 14 complete domain. See I said closed formed solutions are a luxury in solid mechanics. Only very few problems have closed form solutions. And I also alerted for the problems involving crack, we would not be able to get a close form solution. But we would develop near field solution. Near the crack tip, how the stress field is, how the displacement field so on and so forth. Right at this development, it gives you an impression that you are able to satisfy what happens near the crack-tip? What happens at distances away from the crack-tip?

(Refer Slide Time: 37:51)

So in this class, what we had looked at was for crack problems, analytic functions are used as stress functions. And we need to know certain features of analytic functions for us to proceed further. We looked at what are Cauchy-Riemann conditions. And these Cauchy-Riemann conditions are extensively used, when you want to investigate whether a given stress function is a valid stress function. That is what we had

derived. And I had also alerted whatever the problem that we have taken up is a central crack subjected to biaxial loading. You should never forget that. If you look at the solution and simply jump that for hole in a tension strip or elliptical hole in a tension strip, we have looked at only for uniaxial loading. On similar lines, whatever the solution I get, for the crack problem, I will extrapolate it for uniaxial loading. Don’t look at that. Because the boundary conditions clearly shows, you are really looking at what happens at infinity. And here we specify sigma x equal to sigma y equal to sigma. Keep that in mind. That is essential for you to look at the stress field development in the next class.

Thank you.