Mod-01 Lec-04 Introduction to Finite Element Method

Video Statistics and Information

Video
Captions Word Cloud
Reddit Comments
Captions
The last class we were looking at a very simple problem. We were trying to do this problem by using the finite element analysis. We had done the preliminary work which we will summarize later in this class, but let us pick the threads from where we left. We were here. We had defined the complete matrix which we called as stiffness matrix for this problem. We defined the matrices for the elements and then assembled them and we got to this stage and applied the boundary conditions. It is quite clear from this stage as to how this problem can be solved. It is simple that these two matrix gives a set of simultaneous equations, two equations, simple simultaneous equations which can be written as say K1 plus K2 into u2 minus K2 into u3 is equal to zero. What is the next equation? Minus K2 into u2 plus K2 into u3 is equal to P. There are two unknowns, u2 and u3. It is quite straight forward to solve them. The first equation gives rise to a simple equation in an explicit form by defining u2 is equal to K2 divided by K1 plus K2 into u3. Now what we have to do is simple. We have to substitute this value for u2 in the second equation and then write it down as minus K2 into K2 by K1 plus K2 into u3 plus K2 into u3 is equal to P. Once this is done, I have to just sort out or determine what is u3? Solving it is very straight forward. I can write u3 to be equal to P divided by K1 plus P divided K2; P divided by K1 plus P divided K2 and u2 can be written as P divided by K1. We have determined or we have found out what are the displacements for the system? Once I know the displacements, what is my next step? I have to determine the strain because ultimately I am interested in stress. So, what we do is we solve these equations, determine displacement. Then we go to strain through what are called as strain displacement relationship and then ultimately we determine stress through stress strain relation. So, we first determine strain through strain displacement relationship. What would be the strain? There are two strains, epsilon1 and we can say epsilon1 and epsilon2 as the strains that exist in the two bars. In this particular case, we can use a simple formula of change in length by original length. So, I can determine epsilon1. What is epsilon1 in this case? Look at this. Remember that we had called this as node 1 and node 3. So, what is strain in 1? It is nothing but u2 divided by L1. From this relationship, I can define the strain in 1 as u2 divided by L1. So, I can also determine epsilon2 and epsilon2 is written as, now there are two nodes here, so u3 minus u2 divided by L2; very good. Once I know this, I determine the stress by a simple relationship. Again for this problem We have to expand this relationship later for this problem as sigma1 is equal to E1 epsilon1 and sigma2 is equal to E2 epsilon2, which on simplifying, by substituting the values for u2 or u3 from whatever we have derived here in this case and substituting for K1 and K2, a very straight forward substitutions, this is the relationship for u2 and this is the relationship for u1. Let me just do only one of the two. The other I can just write down, it is quite straight forward. So, if we do not waste much time on it, remember what K1 is and what is K2? K1 is equal to A1 E1 by L1 and K2 is equal to A2 E2 by L2. By doing that, it is very easy to show that sigma1 is equal to P by A1 and you can just substitute it and see it is equal to P by A2. You can straight away substitute it and you can see that for example if you have a doubt I can say that we can work it out for this; how we got it. Though it is very straight forward, it is useful to look at that quite closely. Remember that epsilon1 is equal to I am substituting from here; so epsilon1 is equal to u2 by L1 and u2 is equal to P by K1 and that is equal to E1 by L1 into P by K1, K1 being A1 E1 by L1. Substituting that A1 E1 into L1, striking them off, I get sigma1 is equal to P by A1. Same way, you can substitute it, it will take two minutes but you can see that sigma2 is equal to P by A2. This particular procedure which I will summarize in a minute is in a nutshell what we do in finite element analysis. There are very interesting things we have not talked about. Look at that result. There is always a misconception, mostly among students, that when I have a bar like this with two different E’s, E1 and E2, there is always a misconception that when I apply a force to this particular bar, to this bar for example what is here, that since E2 say for example let us assume that E2 is greater than E1, just for an argument. Since E2 is greater than E1, the students feel that the strains or displacements are lower; strains are lower and hence stress is lower, when I have here higher E2. It is absolutely not correct. Why that is, because ultimately we are looking at equilibrium equations. What does this actually indicate, in one sense if these two come out from the equilibrium of the bar? Sigma1 A1, if I bring the A1 to the other side, in a sense if these two come out from the equilibrium of the bar; now, sigma1 A1, if I bring the A1 to the other side, sigma1 A1 signifies the internal forces due to stresses; stresses developed due to internal forces and in one hand they are nothing but an equilibrium equation. So, the force, whatever be the material when I apply a force, this force has to be equilibrated by the internal stresses and when they develop, when this internal stresses develop, they develop in such a fashion that they are independent of E1 and E2. That is a very important thing to realize, because people usually make a mistake and that is why I have pointed out this particular case. Now let me summarize. If we look at what all the steps we have gone through, you have gone through quite a few steps, though every step was simple but you have gone through quite a few steps, what is the first step? We looked at the problem and said that this behavior is a simple bar behavior. We said that this is simple bar. Please note there is a difference between geometry and its behavior. The same bar if I apply another type of load, for example let us assume that this bar has a circular cross section. The same bar if I apply a load in this direction perpendicular to it, would behave or what would it undergo? It would undergo bending. It would undergo bending and hence the behavior of the body would be like a beam. So, the first thing that we did was to immediately recognize from this problem and with the engineering knowledge we have at the back of our mind, the first thing we said that look this behavior is a bar behavior. In fact I told you that if this is a bar behavior, there is an element called bar element. So, I discretise this using the knowledge of engineering which I had gained and correspondingly looking at a thing which is close to the engineering behavior in this finite element library. What does it look like or what is the analogy? You remember that the first day we had a problem of determining the area. When I wanted you to determine the area what did you do? You said that look there is a triangle there. I know what is the area of this triangle? Immediately you discretised that body into a triangle. Similarly you know what the area of a circle is and you apply that knowledge; that is what you did that day. You applied that knowledge and you knew that there is a circle whose area is known and so on and this is exactly what you are doing here. Knowing how the material or the geometry is going to behave under a combination of the geometry as well as loading, you said let us choose bar element first thing. So, after choosing the element, what did you do? You discretised it; you divided it into two elements. Why did you do it? Before you did this, you immediately had an idea of what is stiffness of a bar? That is what we did. What is the stiffness of a bar? Then you realized that there is going to be a term say A1 E1 by L1 which we called as K1 and A2 E2 by L2. After having realized this you also felt that if I divide this into two elements it is nice, because this guy has one area, this guy has another area. So, I can have two elements so that my calculation of stiffness matrix does not become very messy. So, you choose this element again looking at geometry and having knowledge of the stiffness matrix or stiffness determination. That is the second step you did; you discretised it. As a third step what did you do? You wrote down what is the stiffness of each of these elements? In that process, you also defined K1 and K2. What is the next step you did? You assembled them, you assembled them, you put them together so that I get a global stiffness matrix. So, you got what is called as global stiffness matrix. Next is the assembly procedure. After having assembled it, you applied boundary condition. This is a very crucial step in finite element modeling, finite element analysis. Let us take a simple example now apart from this to show the importance of boundary conditions which most not only students, but practicing engineers also get confused. The question they ask is, let me keep this figure because we will explain that again. The question they ask is suppose I have a bar; let us stick to the bar for time being and if I apply a load say P on either side of the bar which is typical, very simple thing; I take this bar go to some testing machine and usually this is the type of free body diagram you would usually put, simple thing. The question that people ask is where is the boundary condition here? Can I not solve it straight away? Can I not solve it straight away is the question. The other confusion is that look why do you want to apply boundary condition because anyway I have u1 u2 and u3 in this previous problem. Similarly I have u1 and u2 in this problem and I have three equations. So, why not I go and solve it? Can I solve this without applying boundary conditions? Because, that is what it looks like from the free body diagram. But unfortunately that is not true. It is Just because you have written down three equations and just because you have three unknowns it does not mean you can solve this equation. because These equations have to be linearly independent, then only it would be possible to solve this problem because you have got these particular equations for the stiffnesses from this the equilibrium, internal equilibrium. So, the solution for this problem as far as say for example displacement is concerned will not be unique. Physically, what does it mean? It means that I have to get the correct positions of these, for example, ends as I apply this by as i apply this bar as I apply the load in this bar. That is not unique; it is something like adding a constant to integration; something like adding a constant to integration. For example, I take this bar and then determine the stresses for this bar in this position. I will get a particular u1 and u2 with respect to say coordinate system. On the other hand, if I take the same bar here, same bar and apply the same loads and determine the displacement and stress it has to be the same. But I would get a different value. This can be say at a distance of around 5 units from the origin. This can be at a distance of 3 units from the ……. So, you would not get a unique u or unique solution. In other words a rigid body mode can be applied. There is a rigid body translation, a rigid body translation that is applied here and these kinds of rigid body translation will not affect the result physically as you see it, but will make it impossible to get a unique solution. We have to arrest this kind of rigid body motions so that we get a unique solution. Look at that again. Then the question that you may ask is how do I solve this problem? This is free body diagram. We are all people with mechanics background, so, we would like to put the free body diagram and then start this problem. Then how do I solve this problem? Where do I put a boundary condition? Yes; one solution, as you say that immediately, why not I fix it here and then solve this problem or these two problems from which I started. Are these two problems the same? You feel that both the solutions are the same; you are just fixing it and when I fix it the reaction is anyway going to be P. That is what you are saying and hence these two problems are the same. Yes, to certain extent what you say is right but the purist will not agree to it for a simple reason that if I take a small element there, in this particular problem corresponding element being here, in this particular problem, when I fix it that element will not move. Can you say the same thing of this element? No; this element would start moving. But what you say is right because the stress anyway will not be affected, stress should be the same. On one hand in this problem if I plot the displacement right across, it would start from zero. On the other hand the displacement here would start from negative; negative value, if I consider this as the negative x1 direction, it will start from the negative value. How do you solve this problem? One of the easiest ways to solve this problem is to exploit what is called as symmetry of the problem. What is meant by symmetry? Now let us look at a material particle. Let us sit here in this material particle and see what is going to happen to us as this load is applied; go and sit at this place. When you sit at this place, this guy here to the left is going to pull me and this guy to the right is going to pull me. I am going to have a horrible stress state. Both the guys are going to pull me. They are going to pull me with equal force. What happens? Equal force; you know both the guys are going to pull me. I am going to stand there. Please note again when I am going to stand there, it does not mean that I have no stress. I am stressed; both the guys are pulling me. All these material particles which are sitting on this line are not going to have any motion which means that instead of solving this problem I can say take half of the problem. I can take that problem fix this side and then apply a load here and then solve this particular problem. We will come. What is the boundary condition, what is the load to be applied, all these things we will come later; just note that the whole lesson here is not symmetry. Symmetry will come bit later in the course but what is important to realize is that we cannot solve a problem without giving boundary condition. That is the lesson I want to look at from this example. We will talk about symmetry which is slightly more complex later in the course, but note this carefully. Hence boundary conditions are important. Once you apply the boundary conditions you can solve the problem, determine the displacements. What is the result you are going to get from finite element? The result that you get from this finite element analysis is displacements. You do not get stress; please note you get displacements. Then you go back and through strain displacement relationship, you find out the strain and then through stress strain relationship you determine stress. Is that clear? This in a nutshell completes all the steps that are required in order to solve this problem. Yes, any questions? Yes, so that is a good question. The first thing is what happens if the bar is inclined? Suppose I want to solve the same problem, this same problem which is sitting here, if I have to solve this same problem with the bar being in an inclined position or for that matter an element being inclined, how do you now solve this problem? Is there going to be any difference in the procedure? This problem becomes important because of two or three reasons or rather most important reason is this. There are number of softwares that are used to do finite element analysis. We will again look at where we use softwares? What is that we use as with respect to different steps that we have seen but immediately these softwares have a tendency to give names. Looking at this they will say three dimensional bar element or generalized bar, three dimensional bar element. This is the name they would give; three dimensional beam, three dimensional bar. These are the names that are usually given. People get confused because as we had just now talked the bar can take only degree of freedom along its axis. So, what do you mean by 3D or displacement of a node? There are three degrees of freedom for the displacement of the node; that is what they would say. What does it mean or how am I suddenly going to get three displacements? Does it mean that this bar element can also take into account bending? No. What essentially it means is that for example if this is the original bar and if this is the displacement and this displacement happens to be say ui and that displacement happens to be uj, this displacement, what essentially you do as far as a stiffness matrix is concerned is to resolve these displacements along the direction of the beam or sorry of the bar. If this happens to be, if this angle happens to be beta and let me define C is equal to say cos beta, then if this is going to be ui then this shortening, keeping say for example let us say uj is equal to zero so that things become easier to understand, then this shortening of length happens to be what? Cui and the force that is developed due to this shortening of length now is to be calculated. Let me call this force as F that will be same as what? AE by L into Cui; if I call that force as F, let me write it here, F is equal to AE by L into Cui and this force is now resolved again in the x and y direction which means that if I call this force that is acting in this direction as say pi and the force that is acting say in this direction as a qi, then pi is what? It is equal to AE by L into C square; C square because Cui into C. So, this becomes AE by L into Cui into C. Please remember C is cos beta. Same way you can get q to be defined as qi is equal to AE by L Cui into, correct, let us call that as sin beta and s or yes? It is possible to develop the stiffness matrix of an inclined bar by following the same steps but by resolving the forces as well as displacement along the axis of the bar. So, my physical behavior of the bar does not change but I change a little, the mathematics that is involved so that I am able to define a more general inclined bar. The first column of the stiffness matrix would now look like, the corresponding stiffness matrix would like C square, CS, minus C square, minus CS corresponding to four degrees of freedom. Please note now that we call this as four degrees of freedom, two in the x and the y direction for this node and two in the x and y direction for this node, though physically, in actuality the bar only shortens or lengthens. So this displacement along the axis of the bar has been resolved and that is what you get. So, the stiffness matrix now becomes four by four and this would result in K which is a four by four matrix. So, that is clear. Fine; now, so any other question? Yeah! Right. See one that is a good question. Again we will discuss that point. Let me repeat that question what he says in a minute. The problem which we have done so far involves only a uniaxial problem, in other words uniaxial stress. Now, there are lots of things which we have to change; the way we are going to think further. What is it that we are going to change or why we should change our thinking? In the uniaxial problem, we did not have much difficulty in understanding stress because stress we said is just E into epsilon, finished it. There is only one stress called sigma and it happens to be or it turns out to be the same as load per unit area P by A1, P by A2 and so on. There are two very simple definitions we used. One is to define strain; just change in length by original length. The other one which we used was just E into epsilon which turns out to be P by A. What happens to a complex situation? Many of the examples which we saw the other day are definitely not in the class of example which we have seen. None of them are one dimensional bar type of elements. All of them were 3D; remember the railway wheel which we saw, it was a 3D problem. Now a three dimensional case where the deformation is not also homogenous in sense that the stress as you had seen varies from one point to another point. How do you define stress? Load per unit area or P by A becomes quite difficult to define. What is P by A? What is that area that we have to take? Pardon; can you take projected area? See, the problem here is suppose I have a sheet with the hole and so on and if I start applying loads, what we call as biaxial loads like this and say like this and so on then, we know that the stresses near the hole may be different from the other region. There can be one hole; it can be circular hole or elliptic hole or whatever it is, they are going to be different from other positions. Now in other words, you know probably you have all studied that there is stress concentration region near this hole; probably all of you would say; a very simple answer you would get. What do you mean by stress concentration region? You would immediately say that look when there is a hole in this part then the stresses near this hole is going to be high; is going to be high. What is this stress? Is it just P by, suppose I say this is say P2 and this is say P1. The question that comes to our mind is how do I now define stress in this part? What is A now, P2 by what A? If you say it is projected area that is what many people say and so the stresses is P2 by A then you mean to say that at all these regions the stresses would be higher? Then it becomes a problem because immediately you would again come and tell me that no, no it is not like that. Very near the ends the stresses would be higher. In fact when the stresses or when the hole is very sharp, stresses may even reach infinity, the singular point and so on. Let us not worry about those things now, but it is important to realize that now the stresses would change as I move away from this corner. If I use one area as you said if I use just the projected area, the normal area, then I would get only one answer. For example in a more complex case of the railway wheel or for that matter the connecting rod or the tyre and so on stresses vary from point to point, from one position to another position. Then what is the area that I am going to define there. Is it just load per unit area? So, we have to put stress in a much more proper perspective and what is it that we are going to do for that? We are going to define what is called as state of stress at a point ultimately. But, before we go to state of stress at a point, we have to define first of all what stress is, in a much more rigorous fashion than just load per unit area. See, in your earlier classes you would have used a coordinate systems; usually you will call this as x y z or zee. American’s call the z as zee; you know why? Because in English z is never defined or never pronounced as z. Take for example z is used in zebra; but, you have z to be, the phonetics is z to be zee. So, the Americans define x y zee. We are not going to get into that complication in language. We are going to define for example these coordinate systems to be x1 x2 and x3. Let us now take a complicated body and then let us say that I apply loads on the body and I have some boundary conditions. Before we go further let us look at the type of loads that happen to be present in engineering components. Loads are due to certain external influence. Right? Loads are due to external influence. What are the types of external influence? Can you classify them? Yes, we can classify them into two categories, one which produces what are called as body forces and the other which produces what are called as surface forces. What are body forces and what are surface forces? There are external influences like gravity or centrifugal force which happens to be present in the external or sorry in the total volume of the body. Though they are external influences, they are present in the volume elements of the body. Their influence is ever pervading throughout the volume of the body. Such forces which act on the volume of the body are called as body forces. On the other hand we have what are called as surface forces. The surface forces, for example contact forces, you have a connecting rod; there is a …… pin which is in contact in one end called small end and the crankshaft, the other end and so on. So, there are contact forces. These contact forces, as an example, gives rise to surface forces or when body is in water, it is the force that acts on the surface and so on. So, the forces that act on the body can be classified into body forces and surface forces. Now, there can be a question. Usually in all our earlier classes I have applied force, something like this. I have put an arrow and said that look this is the force and probably many of you would have called this force as a concentrated force. Where do I fit in this category? Where do I fit in this concentrated force? Is it body force? Surface forces; correct. So, these forces are surface forces. In actuality there is no such thing as a concentrated force and forces do not act just at one point. They act in a very, very, very small area, infinitesimal area; whatever it is, it acts on an area and that has been modified to be a concentrated force for engineering purposes and hence we have this term concentrated loads or concentrated force and that happens to be a part of the surface forces. Stresses are produced due to this kind of forces. Then, what is a stress, how do you define a stress? First thing that you have to realize is that stress is defined say at a point say P. Stress is defined at a point P. How is it defined at a point P? We will come back to our original definition, but before that what we have to do is to imagine that the body can be split. I am not going to Please note that when I cut the body Before I go further we should realize that when I cut the body it does not mean that I am taking an axe or something and cutting the body and what about all the damage you are going to do to the body? No, we just imagine. So, the first thing is that the stress can be defined at any point P of the body. In order to determine the stress P, I have to cut the body by an imaginary plane. I will say that you want the stress at P, so, let me cut it by an imaginary plane and then take say a part of the body out. Let me take part of the body out. You want me to concentrate on P; yes, let me concentrate on P. This is the point P. So, the next thing I will do is now what is the influence of the other body or the other part of the body which you have cut onto the body which you are now having in your hand and looking at it? Please note the body still has all the forces that are acting. It is not that I have removed everything and cut the body but I have just cut the body, taken it imaginarily and looking at this point, P. We will go through that in the next class; to define what stress is from this point of view in the next class. Before that if there are any questions I will answer. Yeah, yeah very good. So, that is a next question; that is a very important question. It is good that we stopped here and you asked the question. The technique that we used in order to arrive at the stiffness matrix is a direct method. There the question is that is this the only way or in other words if I have a tapered bar, if I have some other shape and if it is not going to be bar, but any other shape how do I get a stiffness matrix. Can I use this kind of same approach? The approach that we used is what is called as a direct method. But this direct method may not work for all cases, but we have to resort to other engineering principles. It may be based on energy theorems, it may be based on virtual work principle and so on. We are going to answer that question later in the course. But in order to do that we have to first understand stress and hence we are now going to look at stress very closely. Once we do that and we understand what is stress, what is strain in a much more general setting then, it would be possible for us to go into more involved theorems, define them clearly and then define stiffness matrices.
Info
Channel: nptelhrd
Views: 107,563
Rating: undefined out of 5
Keywords: Introduction to Finite Element Method
Id: _-5W9RvrjvA
Channel Id: undefined
Length: 49min 45sec (2985 seconds)
Published: Tue Jan 07 2014
Related Videos
Note
Please note that this website is currently a work in progress! Lots of interesting data and statistics to come.