Updated 7 June 2020
Reading time: 10 mins

Yielding, Plastic Deformation and Moment Redistribution in Beams (2/2)

Learn how calculate plastic moment capacities and how moment redistribution occurs in a structure
[object Object]
by Dr Seán Carroll

In the previous post in this series we introduced the idea of elastoplastic bending. We showed how the normal/bending stress distribution transitioned from linearly elastic to fully plastic. We saw that a fully plastic cross-section behaved just like a hinge allowing free rotation if further loading was applied to the structure. We finished that post with the clear understanding that a cross-section has additional plastic moment capacity (ductile materials only!)

In this tutorial, we’re going to work out exactly how to determine the plastic moment capacity of a cross section. We’ll also explore the concept of moment redistribution with an illustrative example. By the end of this post you’ll be able to calculate the plastic moment capacity of any cross-section and understand in detail how moment redistribution occurs in a structure.

Plastic Moment Capacity and the Shape Factor

Once we’ve established the fact that a cross-section can go beyond its elastic limit and sustain further loading after it has begun to yield, the obvious question is how much more load can it withstand? We can start by considering a cross-section which has reached its plastic limit.

Plastic moment capacity | EngineeringSkills.com

We know that to satisfy equilibrium, the total force acting on the cross-section must be zero, therefore,

Ft=Fc(1)F_t = F_c \tag{1}

We know force is simply the product of stress and area (over which the stress acts), so,

σyA1=σyA2(2)\sigma_yA_1 = \sigma_yA_2 \tag{2}

From this we can conclude that when a cross-section has fully yielded, the neutral axis divides the cross-section into two equal areas. So, if the z-axis is not an axis of symmetry (and it’s not in the case of the trapezoidal cross-section shown above), the neutral axis under plastic conditions will be in a different location to the neutral axis under elastic conditions. We can now state an expression for the plastic moment as,

Mp=Fcyˉ1+Ftyˉ2(3)M_p=F_c \bar{y}_1 + F_t \bar{y}_2 \tag{3}

Now since,

Ft=Fc=σyA2(4)F_t = F_c = \frac{\sigma_y A}{2} \tag{4}

We can state that the plastic moment is,

Mp=σyA2(yˉ1+yˉ2)(5)M_p = \frac{\sigma_y A}{2}(\bar{y}_1 + \bar{y}_2) \tag{5}

We now define the plastic section modulus as,

S=A2(yˉ1+yˉ2)(6)S = \frac{A}{2}(\bar{y}_1 + \bar{y}_2) \tag{6}

And we can state the plastic moment as,

Mp=σyS(7)M_p = \sigma_yS \tag{7}

It should be noted that the plastic modulus is purely a function of the cross-section shape.

Shape Factor

We know from basic beam bending theory that the elastic moment is given by,

My=σyIy(8)M_y =\frac{\sigma_y I}{y} \tag{8}

and that the elastic modulus is given by,

Z=Iy(9)Z=\frac{I}{y} \tag{9}

Now we define the ratio of plastic moment, MpM_p to yield moment MyM_y  as the shape factor, ff and note that this ratio is ultimately a function of the cross-section shape,

f=MpMy=SZ(10)f=\frac{M_p}{M_y} = \frac{S}{Z} \tag{10}

💡 The shape factor is a measure of how much reserve strength is in the cross-section after yielding begins.

The shape factor is a helpful indicator of how much more loading can be sustained after yielding but before a plastic hinge is formed. Remember, it is simply a function of the shape of the cross-section and is independent of the material properties of the cross-section.

Calculating the Plastic Moment Capacity

Now we’ve covered all of the theory, let’s actually calculate the plastic moment capacity and shape factor of a cross-section. To make things a little more interesting we’ll consider a compound cross-section consisting of an I-beam with a steel plate welded to the bottom flange. We’ll assume a steel yield stress of σy=350N/mm2\sigma_y = 350 \:N/mm^2.

Plastic moment capacity example | EngineeringSkills.com

You can watch me walk through the full solution to this question in the video below (the thought of typesetting the full thing in this post was too much to take!) But the outline steps are as follows:

  • We need to calculate the second moment of area using the parallel axis theorem. But as this is a compound section and not symmetrical about the neutral axis, we’ll need to first locate the neutral axis, by taking area moments about an arbitrary axis.
  • Armed with the second moment of area and neutral axis location, we can plug these into the engineer’s bending equation (equation 8 above) to determine the moment required to cause yielding, MyM_y .
  • Then we can move on and determine the location of the equal area axis. Remember the neutral axis becomes an equal area axis (and may shift location) in the fully plastic condition.
  • Once we have the equal area axis, we can work out the plastic moment capacity and apply equation (10) to determine the shape factor.

It’s worth you trying to solve this yourself first before watching me solve it in the video below.

Plastic Hinges in Beams and Moment Redistribution

Now we’re going to expand our study of elastoplastic behaviour beyond the cross-section to look at hinge formation in a beam. Consider the simply supported beam below subject to an increasing point load at mid-span. By examining the bending moment diagram, we can map out and visualise the formation of a plastic hinge under the point load.

Hinge formation in a beam | EngineeringSkills.com

We can clearly see how the bending moment diagram can be used to identify the extent of the elastoplastc region in the beam. As we now know, once the entire section becomes plastic (as indicated at the mid-span in the beam above) a plastic hinge forms.

Now we’ll consider how hinge formation allows moment to be redistributed in a statically indeterminate structure. We demonstrate the fundamentals of this behaviour by considering a simple (indeterminate) propped cantilever.

Propped cantilever plastic hinge 1 | EngineeringSkills.com

It can be shown that the bending moments at A and C are:

MA=3PL16(11)M_A = \frac{3PL}{16} \tag{11}
MC=5PL32(12)M_C =\frac{5PL}{32} \tag{12}

We’re going to set ourselves a simple task; to determine the value of the applied load, PP at which the structure collapses. In doing so we’ll need to take account of plastic hinge formation and moment redistribution.

We can start by stating some basic assumptions:

  • Let the span L=2mL=2 \:m
  • Let the yield moment My=12kNmM_y = 12\: kNm
  • Let the plastic moment Mp=15kNmM_p = 15 \:kNm

Let’s also assume the the beam must accommodate an applied working load P=25kNP = 25 \:kN.

Analysis for P=25 kN

Under normal working conditions, the moments at A and C are:

MA=3×25×216=9.4kNm(13)M_A = \frac{3\times 25\times 2}{16}=9.4 \:\:kNm \tag{13}
MC=5×25×232=7.8kNm(14)M_C = \frac{5\times 25 \times 2}{32} = 7.8 \:\:kNm \tag{14}

Both moments are less than the yield moment therefore the beam is safe under working conditions. So let’s push things further – remember we want to know the load at which the beam collapses.

Analysis for P=32 kN

Again determining the moments in the structure:

MA=12kNm=My(15)M_A = 12 \:\:kNm = M_y \tag{15}
MC=10kNm(16)M_C = 10\:\:kNm \tag{16}

Therefore the cross-section at location A has reached its yield limit; i.e. the extreme fibres in the cross-section at this location have reached their yield stress. In a linearly elastic analysis we would call this the maximum allowable load and end things here. But as we know, this beam can take more, so let’s push things further still.

Analysis for P=40 kN

In this case, the moments are:

MA=15kNm=MP(17)M_A = 15 \:\:kNm = M_P \tag{17}

At this point, a plastic hinge has formed at A. Any further loading will cause free rotation to occur at this location.

MC=12.5kNm>My(18)M_C = 12.5\:\:kNm >M_y \tag{18}

At C, the section has exceeded its elastic limit and yielding has started. At location C the section is undergoing elastoplastic bending.

Propped cantilever plastic hinge 2 | EngineeringSkills.com

It’s important to remember that the structure started out as indeterminate to 1 degree i.e. it had one redundancy. It has just lost this redundancy with the formation of a plastic hinge at A. The structure will not collapse yet, but if one more hinge forms, the structure will become a mechanism and will collapse.

Analysis for P=45 kN (Moment Redistribution!)

At this load, the moment at A remains 15 kNm. Remember, the section can’t withstand any more moment once a hinge has formed. Any further loading after hinge formation simply causes rotation about the plastic hinge.

So if any additional load is applied, we must have moment redistribution in the beam. This is inevitable; the load is continuing to increase, but a hinge has formed at A allowing the section to rotate at that location. If the structure is to avoid collapse, the increased moment within the structure must be redistributed from the support at A and into the span, further increasing the moment at C.

The formation of a plastic hinge at A means that the structure now behaves as a simply supported beam. In other words, all of the load applied between P=40 kN and P=45 kN is resisted by the structure as if it were a simply supported beam.

Propped cantilever plastic hinge 3 | EngineeringSkills.com

Now determining the moment at C:

Mc=12.5+2.5=15kNm=MP(19)M_c = 12.5 + 2.5 = 15 \:\:kNm = M_P \tag{19}

We can see that at this value of applied load, a hinge will form at C. The structure now becomes a mechanism and collapses.

Propped cantilever plastic hinge 4 | EngineeringSkills.com

The load at which collapse occurs is 45 kN while the load at which yielding first occurs is 32 kN. So this structure could withstand approximately 40% more load before collapse than an elastic analysis would suggest.

This ability to safely sustain loading after plastic hinge formation is often utilised in structural design. Large span steel portal frames are often designed to develop plastic hinges. Care must obviously be taken to ensure sufficient redundancy remains to ensure collapse does not occur.

Plastic Collapse Analysis

We can define a load factor λi\lambda_i as the collapse load for mechanism ii, divided by the structures’ working load,

λi=collapseloadformechanismiworkingload(20)\lambda_i = \frac{collapse\:\:load \:\:for \:\:mechanism \:\:i}{working \:\:load} \tag{20}

For a structure with multiple collapse mechanisms, the collapse load factor is the smallest value of the load factor – the load factor at which the structure will collapse. For our example above, the collapse load factor is,

λC=4525=1.8(21)\lambda_C = \frac{45}{25} = 1.8 \tag{21}

We might explore plastic collapse analysis in another tutorial series. But for now you should have a good understanding of the material and cross-section behaviour that leads to collapse of ductile structures and the additional capacity that exists, after material yielding but before plastic hinge formation.

In case you missed it, the first tutorial in this 2 part series can be found here.

getting-started
Dr Seán Carroll
BEng (Hons), MSc, PhD, CEng MIEI, FHEA
Hi, I’m Seán, the founder of EngineeringSkills.com (formerly DegreeTutors.com). I hope you found this tutorial helpful. After spending 10 years as a university lecturer in structural engineering, I started this site to help more people understand engineering and get as much enjoyment from studying it as I do. Feel free to get in touch or follow me on any of the social accounts.

Dr Seán Carroll's latest courses.

Analytical Modelling of Plate and Shell Structures: Part 1 - Plates

Analytical Modelling of Plate and Shell Structures: Part 1 - Plates

A practical guide to the analysis of circular and rectangular plates under load, from first principles.

Fundamentals of Reinforced Concrete Design to Eurocode 2

Fundamentals of Reinforced Concrete Design to Eurocode 2

An introduction to ultimate limit state design for bending and shear with optional calculation automation using Python.

Modelling and Analysis of Non-linear Cablenet Structures using Python and Blender

Modelling and Analysis of Non-linear Cablenet Structures using Python and Blender

Learn how to combine parametric modelling, exploratory form-finding and iterative analysis techniques to simulate 3D tensile structures.

Non-linear Finite Element Analysis of 2D Catenary & Cable Structures using Python

Non-linear Finite Element Analysis of 2D Catenary & Cable Structures using Python

Build an iterative solution toolbox to analyse structures that exhibit geometric non-linearity due to large deflections.


Do you have some knowledge or expertise you'd like to share with the EngineeringSkills community?
Check out our guest writer programme - we pay for every article we publish.