Part 2 - Section Analysis

Section analysis

As I mentioned in the introduction, I was taught many different methods of section design, one for each material type, and then other versions at serviceability and at ultimate, and other variations with axial forces thrown in. But there is only actually one method.
Surprised?

The only difference between the various situations is the behaviour of the particular combination of materials.

Stress-strain diagrams


The behaviour of any particular material is best understood by looking at a plot of stress against strain. Different codes often use different stress strain curves, but in each case they are just the code author’s particular way of simplifying actual test results.




Figure 4 - BS5400 stress strain diagram for reinforcement (modified as described)






Figure 5 - BS5400 stress strain diagram for concrete (modified as described)

The stress strain curve for structural steel can be quite hard to find, and I could not find one in BS 5400 or in the Eurocodes. The following is therefore an unofficial modification of the graph for reinforcing steel.


Figure 6 - idealised stress strain curve for structural steel (tension side only shown)

For simplicity’s sake, the main graph has been changed to only have two lines, a linear elastic part and a plastic part. I have, however, also drawn in some subtleties in thin lines to aid your understanding of a few other matters.
Naturally, you will have heard the terms linear and elastic, but let us assume that you have just started and define them for you. Linear behaviour means that the graph of any particular effect, in this case stress vs. strain, is a straight line. For this particular example, the line also goes through zero, so we can say that doubling the input (strain) will double the output (stress). Elastic means that the stress state (i.e. the relationship between stress and strain) will follow the same path back down the curve (or straight line in our case) as it did on the way up: think of it as the behaviour of a spring which always goes back to its original length when it is unloaded.

On the linear elastic section:
stress = strain x Young’s modulus, which is a rearrangement of the definition of Young’s modulus.

On the plastic part of the graph, the strain in the material keeps on increasing without a corresponding increase in stress, i.e. it stretches and stretches, without any extra force required.

This is obviously a simplification of the more realistic graph, shown with the thin grey line, but it is the normal assumption for working with structural steel.

Subtleties: a more realistic stress stain curve would not just be a pair of lines with a sharp break. When the steel starts to yield, the corner of the graph is not sharp, but a gradual curve, as shown in grey. So how do we determine the yield stress ? Think back to the description, above, of elastic behaviour. If we consider the stress strain curve to be straight, followed by a short curve, and followed then by another straight section, then, in the elastic range, loading and unloading a sample will go up the straight line and then back down the same way. But one we get onto the curved section, unloading a sample will immediately start back down a different line, parallel to the linear section. This is where the sample has started to go plastic and when the sample is unloaded, there will be a small permanent strain. Yield is normally taken at a particular “Proof Stress”. As an example, if the critical value is the 0.2% proof stress, then the value will be the stress (partially around the curved bit) where fully unloading the sample will return to a permanent strain of 0.002.

Likewise, the flat plastic section of the graph will probably have some strain hardening, which is a posh phrase for the fact that the stress will actually increase a bit with increasing strain, although the rate will be orders of magnitude less than for the linear elastic part of the curve.

Using stress diagrams in section analysis


I am afraid that the next bit will require a small amount of visual gymnastics. It is something that some people see straight away but it is so simple lot of people cannot get their heads around it: but it is just a way that I use to visualise how things work. The method is shown in Figure 7via three examples with progressively higher peak strains.

In case you can see through the principle straight away, or just as an introduction, the following demonstrates that the relevant part of the stress-strain diagram, from zero to the maximum strain, can be transposed on to the stress diagram.


Figure 7 - steel in compression



During any section analysis, and here we are just looking at the compression side to prove a point, the strain will vary linearly from zero at the neutral axis (by definition) to a particular value of strain (to be determined in a later part of the calculation) at the tip of the section. The compression strain diagrams for the three examples are shown as (1). Variations on the stress strain diagram for the material are drawn to the same scale in (2), directly above (1). Since the strain scales are the same, the stresses at the tip of section will be the stresses corresponding to the points directly above the extreme values of the strain diagram. If we “Wipe out” the rest of the stress strain diagram (shown light and dotted to give a visual reference) then we can see how much of the stress strain diagram is relevant for the particular maximum values of strain: these sections have been coloured in a tasteful shade of pink. It should not be too hard to see that the pattern of stress that is highlighted should be mirrored, rotated (look at the order of the letters in the corners of each shape: the x and y axes have merely been swapped to suit the way that we normally look at sections) and “Stretched” to fit the plot of stress versus height. In the examples above, the pink parts of the stress strain curve have just been repeated on the stress diagram at exactly the same scale: the first example is purely elastic, the second is elastic plastic and the third is mainly plastic: I will explain later why the plastic section is actually almost completely plastic but I need to keep some teasers.

This principle of superimposing the relevant part of the stress strain curve onto the stress diagram is the basis of all of the various methods of section analysis and means that the same principle can be used for all forms of section analysis.

Balancing the forces


I have introduced the complicated bit (not) and there is only one very simple bit left. Phrase it however you want – resolved forces must total zero – what goes in comes out – there ain’t no such thing as a free lunch – but the idea is that the total of all the internal forces must equal the externally applied force. In real money that means that in a reinforced concrete section in pure bending, i.e. with no axial force, the sum of the compressive force in the concrete (positive) and the tension force in the steel (negative) equals the applied force (zero). If there is an applied force as well, e.g. in a column, then the sum of the compressive and tension forces equals the applied axial force.

Reinforced concrete at serviceability under pure bending


This is not the easiest of the strength analyses, but it does demonstrate almost everything discussed above in a clear cut form.

When working at serviceability, the stress in the concrete is relatively low and it will remain in the straight line (linear) part of the stress/strain curve. Because the behaviour is linear, the stress is simply equal to strain x Young’s modulus.

The diagram of strain against depth, as always, is a straight line, hence it can be defined by any two points. The easiest point to assume is a particular value for the compressive strain at the “Top” of the section. I use the word “Top” because that is the case for the most common case of a simply supported beam, but the orientation is irrelevant and we could just as easily be talking about bottom for a beam that is hogging, or side for a beam that is being bent sideways.

Fundamental: this is a classic case of the difficulty with some explanations, demonstrated by the continuous bridge aside discussed earlier. It is very hard to describe something without putting in an orientation, but it is vital to be able to keep everything general in your mind, as far as possible. Remember that the value of shares can go up or down but gravity always sucks.

The second point we need to assume is the depth in the section where the stress is zero, i.e. the depth from the top to the neutral axis (which is not the same as the centroid if you remember).

The fundamental assumption with the behaviour of concrete is that it cannot take any stress in tension. This is, in fact, a slight exaggeration since concrete can take a small tensile stress before it cracks, but since concrete can crack without any externally applied strain, e.g. through shrinkage cracking, then no tension stress can be relied on so it is ignored: and purely by coincidence, it makes the arithmetic easier too so it is a win win assumption.

Aside: one of my favourite instructions for my grads is to be lazy in everything we do – this is nothing to do with being slapdash, but is more about being efficient. The above is a nice example of simple being efficient at the same time as being rigorous – yeah us !

The stress diagram for the concrete forms a triangle “Above”the neutral axis. The tension reinforcement is (hopefully) on the tension side of the neutral axis. Once we have assumed a peak concrete strain and a depth to the neutral axis, the strain in the steel can be calculated using similar triangles and hence the steel stress can be calculated. The force in the concrete, and the force in the steel can now be calculated, and the two assumed values varied until the total force is zero. It works out that, for a rectangular section, the depth of the neutral axis only depends on a series of knowns, i.e. the depth to the reinforcement, the proportion of steel and the ratio between the Young’s moduli of steel and concrete. The peak strain in the concrete, and hence the stress, as well as the stress in the steel, follow on from the depth of the neutral axis and then depends on the applied moment.

The simplest hand calc for this case is called the modular ratio method. The way I was taught used a chart from Reynold's Reinforced Concrete Designer Handbook. This involved working out any one of a number of values, drawing a line across a chart, and reading off the remaining values. The way it was used depended on whether we knew the amount of steel and wanted to calculate the stresses, or whether we knew the moment and wanted to design the steel.

But it all boiled down to one equation


where x is the depth to the neutral axis

d is the depth to the centroid of the reinforcement (please note that d is measured to the steel centroid and that the total depth, when used, is usually referred to as h. One common usage of h is in BS 8110 charts for columns under combined bending and compression)

m is the modular ratio Esteel / Econc
ρ is the proportion of steel i.e. Asteel / (b.d)

b is the width of the concrete.

This can easily be fitted into a spreadsheet as the basis for any calculation. Since the concrete stress block is triangular, its centroid is 1/3rd of the depth from the top, so the lever arm


which is then divided into the moment to calculate the force in the steel or concrete.

Having told you to question everything (and if I haven’t done so then I am doing so now) you should now be asking me to prove this. What is that I hear ? Ok then …

Let strain at top of concrete = εconc and strain in reinforcement = εsteel .

By similar triangles:

The force in concrete has a triangular pattern so the average is half the peak:

Force in steel:

The two forces must be equal, and substituting for εsteel from above:


Multiplying the x by d/d (aka 1), to make d . (x/d), cancelling out the εc’s and moving the Ec and 2 to the right hand side:

The area, down to the tension steel (which is all we are interested in for a pure bending case) equals b.d, and dividing As by b.d we get the proportion of steel ρ. Also, Esteel divided by Econcis the modular ratio m so:


Almost there.

If we now multiply each term by (x/d) then we get a quadratic equation in (x/d):


Those of you who are racing ahead will recall the solution of:

So:

Cancelling out the 2’s, including the 4’s in the square root we get:

Which then rearranges to:
Round of applause please.

The nice thing is that you only have to check my proof and then remember that the equation exists … and that the method boils down to a straight line strain graph, application of the relevant stress strain curves and a balancing of the forces.

From the calculation of x/d, you can carry this on to work out the forces in the concrete or steel, using the lever arm z defined earlier, and check your material stresses against the relevant code limits. If you need to do a crack width calculation, you can extend the strain calcs down to the level that you are interested in and go on from there. Etcetera etcetera etcetera.

The formulae above only work for a rectangular section. But if you have a complex shape, like a voided deck, exactly the same rules are used, you just have to use a tool like a spreadsheet to work out the stresses, and hence forces, over a series of (say) horizontal strips based on assumed strains, and then vary those strains until the forces balance and the moment from the internal stresses equal the externally applied moment.

Simples.

Hint: you can use any tools that you want, but I usually find that the spin button very useful. This is a control that can increase or decrease the integer value in a control cell between specified values. The value limits can be set to anything between 0 and 30,000, but those are the outer limits. If I wanted to set a value for the depth to the neutral axis accurately then I would choose a maximum of 2/3rds of the depth (even halfway down is probably too deep since it means that there is too much reinforcement: ultimate failure of an over reinforced section is sudden and catastrophic – bad – as opposed to a lightly reinforced beam which forms large cracks and fails with warning). I would then do the calculation using a depth that was x/30,000 of that 3/5ths of the section depth and vary the controlling value using the spin box. The example does this with the depth set to x/50,000. You can set as many controls as you want, e.g. for the depth to NA and for the concrete strain, and if I have a lot of variations to do I often create three spin boxes to control the same cell and use them for coarse, medium and fine adjustments since each button can be set to change the control cell by a different integer number.
If you have never come across the spin control, add the developer toolbar in file, option, customise ribbon and then click on the developer tools to add it to the toolbar.

Reinforced concrete at serviceability with bending moment and axial force

The principal is exactly the same if an axial force is applied as well. The only difference is that there is no easy formula and you either need to numerically find the answer by varying the two controlling strain values (at the two extremes of the section) until the force and moment balance out.

Fundamental: this is where the difference between the neutral axis and the centroid is very important. The neutral axis can move up and down depending on the amount of rebar and the combination of moment and axial force. But any calculations that you have done to determine the axial force will have assumed that it runs down the centroid, and should therefore be applied there during any section capacity calculation. It is therefore easier to calculate moments about the centroid with the axial force not coming directly into the answer.

Reinforced concrete at ultimate under pure bending


This one is either very easy, or a bit detailed, depending on how rigorous you want/need to be. The normal behaviour for reinforced concrete has the concrete failing first, although it is good practice for the rebar to reach yield first (and some codes reduce the calculated capacity by a percentage if the rebar does not yield).

A typical code approach specifies a particular crushing yield strain at the extreme compression fibre. The stress diagram follows the full range from zero at the neutral axis to compressive yield at the extreme fibre. As can be seen in the concrete stress strain diagram in Figure 5, there is usually a short parabolic part ending in a larger rectangular part at full yield.

In BS5400, the maximum stress equals 0.45 fcu and can be applied rigorously, by equation or by slices in a spreadsheet, or a reduced value of 0.40 fcu can be applied over the full depth. For a rectangular section there is a simple formula using the uniform stress pattern: it gives a very similar result for a rectangular section to the rigorous method. There is also a second check which does not allow the neutral axis to be taken below halfway. However, on a bridge with a voided deck, like a box girder, it is usually better to take the higher 0.45 fcu over the top flange, since the loss in stress over the parabolic section is limited to the webs, which are much narrower. BS5400 also has rules which require an additional 15% capacity, i.e. you divide by 1.15, if the steel has not reached the plateau of the stress strain diagram.

In the American bridge code AASHTO, the full peak value is applied over a rectangle, but the rectangle is assumed to stop at a particular proportion of the depth to the neutral axis which depends on the strength of the concrete. For lower concrete strengths, this is 85%of the depth so effectively they are taking a uniform stress equal to the peak, but with nothing over the first 15% from the neutral axis. You can still apply the stress routine rigorously and, again, this is more important for voided sections.

In the Eurocode, you seem to be able to do whatever you want, with either a parabola/rectangular profile, triangular/rectangular or a rectangle over most of the depth (similar to AASHTO).

Reinforced concrete at ultimate with bending moment and axial force


This is another one of the nasty ones where you have to do it the hard way, by finding the strains at the relevant locations. Again, the force is applied along the centroid, so that is the best place to take moments about. Otherwise it just follows the general methods for concrete.

Structural steel at serviceability under pure bending


At serviceability, most codes will limit the steel to some relatively high proportion of the yield stress, but not actually as far as yield. We are therefore working on the straight section of the stress strain curve from the compression side to the tension side.

Elastic analysis of steel sections is probably the easiest of all the section analysis techniques that you have been taught. You work out the Second Moment of Area, aka I, divide this into the applied moment, and then calculate the stress at any level by multiplying by the distance from the centroid.
Easy isn’t it?

But have you ever wondered why the centroid and why the Second Moment of Area?

Aside: this is where the camera in your computer has just been automatically turned on so that I can take a picture of the puzzled look on your face as you try to remember if you ever knew.

It is actually quite easy and follows on from three facts:

· the strain diagram is straight,

· since the stress strain curve for steel at serviceability is straight, then the stress versus height plot is also straight (or we will fail the serviceability check), and

· there is no external axial force, so the total force in the section is zero.

Are we detecting a common thread here by any chance?

Let us start from an hypothesis that the stress at the centroid actually is zero. If this does not conflict with any of the rules then we can say that the hypothesis is correct.

With a straight strain line, and zero strain at the centroid, the strain at a level y above the centroid is proportional to y, and equal to k.y. The stress σ is therefore E.k.y. Please remember that signs are important and that y is negative below the centroid and the stress is therefore negative. The small amount of force δF on a small area δA is the product of stress x area so

If we add up all the little bits of force then the total force

Dividing both sides by E . k (which are non-zero constants) we get:

Since the total applied force F = 0 then F / (E . k) = 0 and

So what? Well if you think about it, the bit on the left is called the First Moment of Area and the only place about which the First Moment of Area is equal to zero is the centroid. If you have problems getting your head round that one, then use an analogy and change to the centre of gravity theory and you only get all of the products of weight times distance equal to zero at the balance point, or centre of gravity.

Hence, to have a zero total force the stress must be zero at the centroid, is proved.

That is not a circular argument, honestly, it is putting forward an hypothesis and checking to see if all the known criteria are met. Since they are, the hypothesis is proved to be true.

The next, and actually simpler, stage is to calculate the moments that arise from the calculation above.

If we now sum all of the little bits of (force x lever arm from the centroid) then

Those of you who are still (miraculously) awake, might have noticed that the last integral is the Second Moment of Area aka I, so the above can be replaced by

But from above:

Substituting into the moment equation above:

Which is the fundamental equation that we all use to get stress in steel at serviceability.

Ta Daaa.



TO BE CONTINUED …

Plastic behaviour of steel and steel with axial forces need to be added.

This would be followed by prestressed concrete
AKP 28/10/2012

No comments:

Post a Comment