Calculating and Interpreting the Second Moment of Area

by Callum Wilson
|
Published: 29th July 2023
|
Tutorial

Welcome to this tutorial which is focussed on the second moment of area. By the end of the tutorial, you will be comfortable explaining what the second moment of area is, why it’s important to engineers, how to calculate it and how to interpret the values in the context of civil and structural engineering.

1.0 Why is the second moment of area important?

Structural engineers ultimately have two decisions to make when designing individual structural members:

Decision 1 – What is the geometric profile of the structural member – its cross-section shape should be used?

Decision 2 – What are the structural member’s material properties?

For all structural members, both the geometric profile AND the material properties play an important role in determining if the member is suitable to deal with the design effects i.e. bending moments, shear forces, axial forces or torsion. So, if we want to increase the structural capacity of a structural member, shall we improve the geometric profile, improve the material properties or both?

Stiffness is a structural member property which we can easily describe. We will consider two descriptions of member stiffness:

Axial Stiffness is the resistance of a member to axial deformation under the action of axial loading. We say that axial stiffness = \frac{EA}{L}

Bending Stiffness is the resistance of a member to bending deformation under the action of bending moments. We say that bending stiffness = \frac{EI}{L}

For our discussion, let us assume that the material we intend to use is steel. Steel has an Elastic Modulus, E = 210000 MPa.

So if we want to increase the axial and bending stiffness of our structural member, what shall we do? Firstly, let’s assume that the length of the member needs to remain unchanged in order to fulfil its function – so we can’t reduce L. Secondly, let’s assume that all grades of structural steel (S275, S355 and higher) all have E = 210000 \:MPa – so we can’t increase E.

What options does this leave us with? The only other way to increase axial and bending stiffness is to increase the area (A) and the second moment of area (I).

Increasing the area (A) and hence the axial stiffness is intuitive as we are familiar with calculating the cross-sectional area of different shapes. It doesn’t matter where within the cross-section we add the ‘extra area‘, the axial stiffness is increased in proportion to the increase in area – Simple

Increasing the second moment of area (I) and hence the bending stiffness is less intuitive and often less well understood. The rest of this tutorial will improve our intuition for and understanding of the second moment of area. This will give us the tools to make sensible recommendations to improve the bending stiffness of structural members.

The second moment of area is also referred to as the moment of inertia (just a fancy name!). We designate the second moment of area with the symbol I. I always has the dimensions of \text{Length}^4. We must be careful to ensure that we are always consistent with our input measurements when calculating the second moment of area. In this tutorial, we will exclusively use millimetres (mm).

2.0 What is the second moment of area?

So far, we have recognised that the second moment of area is the property which governs the bending stiffness of structural members (assuming a given length and elastic modulus).

The term second moment of area contains two terms which we are familiar with. We know that area describes the quantity of structural material in a given region and that a moment refers to the product of a physical quantity and a distance.

Putting these two terms together, we say that the second moment of area is a description of the distribution of area about a specified axis.

There are two further fundamental statements about the second moment of area which we should keep in mind throughout this tutorial:

1. The second moment of area is a geometric property. If we change the material type or grade, the second moment of area DOES NOT change.

2. As we’ve already mentioned, the second moment of area has dimensions of length to the fourth power, L^4.

3.0 Geometric Axes and Sign Conventions

It is important that axes are clearly and consistently labelled when calculating the second moment of area of a section. The following standard axis definitions are adopted here, which are derived from the right-hand axis system and right-hand screw rule conventions.

img1_second_moment_of_area | DegreeTutors.com
Fig 1. Axes convention

When we refer to bending actions we must always define the axis about which the bending occurs.

img2_second_moment_of_area | DegreeTutors.com
Fig 2. Bending moment convention.

4.0 Locating the Centroid

The precursor to calculating the second moment of area of any shape is to define the position of its centroid. The position of the centroid gives us the point through which the centroidal axes pass. The centroidal axis is the axis about which we determine the global second moments of area.

Let’s consider a simple shape – a rectangle. If we wanted to balance a 2D rectangle on a single spike, where should we place the spike? 

img3_second_moment_of_area | DegreeTutors.com
Fig 3. A simple 2D shape balanced on a spike.

As a hint, consider three different options – ‘A’, ‘B’ and ‘C’…..

img4_second_moment_of_area | DegreeTutors.com
Fig 4. Option A
img5_second_moment_of_area | DegreeTutors.com
Fig 5. Option B.
img6_second_moment_of_area | DegreeTutors.com
Fig 6. Option C.

Intuitively we would place the spike at Option B, right? In fact, what we have done is located the centroid of this simple shape….without any formula! But what if the shape is not a simple rectangle?

img7_second_moment_of_area | DegreeTutors.com
Fig 7. A compound shape.

Not so simple, right? To formalise this idea, we can introduce the following formula, which will allow us to define the centroid of any shape, simple or otherwise. The centroid is a single point in 2D space. Therefore we need a Y-coordinate and a Z-coordinate to define its position.

    \begin{align*}Y_c &= \frac{\sum Ay}{\sum A} \\\\Z_c &= \frac{\sum Az}{\sum A}\end{align*}

We will consider the compound section example above to understand how these equations are used.

First, we need to define our ‘reference origin’. We will measure the ‘y’ and the ‘z’ from this point. The position of the reference origin is arbitrary and can be selected based on user preference. To simplify notation, our reference origin will be placed at the bottom left corner of the overall geometry; this means that y and z are always positive.

Second, we need to split our complicated geometry into square/rectangular regions, being careful not to overlap regions (as this would mean that the same areas would be counted twice). For this example, we have three rectangular regions.

img8_second_moment_of_area | DegreeTutors.com
Fig 8. Splitting the compound shape into rectangular regions.

Ok, so we are now ready to start using our formulae.

We will begin by calculating Y_c, which is the Y-coordinate of the centroid of the overall geometry. To do this, we invoke the first formula Y_c = \frac{\sum Ay}{\sum A}. This formula tells us that we need to find the sum of the products of the area of each region (1/2/3) and its perpendicular distance from the region’s local centroid (C^1/C^2/C^3) back to the origin (parallel to the Y-axis).

Note: the local centroid of each region is found following the same logic used when placing our spike at ‘Option B’

    \begin{equation*}Y_c = \frac{(\textrm{Area 1 x Distance (Y) from Centroid of 1 to Origin)+(Area 2 x Distance (Y) from Centroid of 2 to Origin)+ (Area 3 x Distance (Y) from Centroid of 3 to Origin)}}{\textrm{(Area 1)+(Area 2)+(Area 3)}}\end{equation*}

    \begin{equation*}Y_c = \frac{(60\times10\times5)+(60\times10\times30)+(20\times10\times65)}{(60\times10)+(60\times10)+(20\times10)} = 24.29mm\end{equation*}

Next, we will calculate Z_c, which is the Z-coordinate of the centroid of the overall geometry. In a similar fashion, we invoke the second formula Z_c = \frac{\sum Az}{\sum A}. Here we are taking the area of each region (1/2/3) and multiplying it by the perpendicular distance from the region’s centroid back (C^1/C^2/C^3) to the origin (parallel to the Z-axis).

    \begin{equation*}Z_c = \frac{(\textrm{Area 1 x Distance (Z) from Centroid of 1 to Origin)+(Area 2 x Distance (Z) from Centroid of 2 to Origin)+ (Area 3 x Distance (Z) from Centroid of 3 to Origin)}}{\textrm{(Area 1)+(Area 2)+(Area 3)}}\end{equation*}

    \begin{equation*}Z_c = \frac{(60\times10\times30)+(60\times10\times65)+(20\times10\times70)}{(60\times10)+(60\times10)+(20\times10)} = 50.71mm\end{equation*}

Ok, so we have calculated the coordinates for the centroid of the overall geometry. If we wanted to balance this 2D geometry on a spike we would place the spike at [Y_c, Z_c].

Let’s draw the centroid onto our diagram.

Note: you might have noticed that the centroid of the shape does not lie within the shape itself… interesting right!

img9_second_moment_of_area | DegreeTutors.com
Fig 9. Centroid of the compound shape.

5.0 Second Moment of Area of Rectangles

For any two-dimensional geometry that lies within the Y-Z plane there are three different components which are required to describe the second moment of area:

(1) The second moment of area about the Y axis:

    \begin{equation*}I_{yy} = \int_{A}^{}z^2 \mathrm{d}A\end{equation*}

(2) The second moment of area about the Z axis:

    \begin{equation*}I_{zz} = \int_{A}^{}y^2 \mathrm{d}A\end{equation*}

(3) The product second moment of area:

    \begin{equation*}I_{yz} = \int_{A}^{} yz \mathrm{d}A\end{equation*}

Initially, these integral expressions may appear daunting, but once again, we will consider a simple example. This will reveal that in fact, these integrals are just a mathematical way to say ‘the distribution of area about a given axis’.

As civil and structural engineers, we will work primarily with sections that can be split up into square/rectangular regions (in the same way that we divided the geometry in §4.0). Therefore, we consider square/rectangular shapes in the following derivations.

So let us consider a simple rectangular geometry of height, h, and breadth, b. From §4.0, we know that the centroid of the geometry is at Y_c = \frac{b}{2}, Z_c = \frac{h}{2}. We place the origin of our Y-Z axes at the centroid.

img10_second_moment_of_area | DegreeTutors.com
Fig 10. A rectangular shape placed onto a y-z axes.

We will begin with I_{yy}:

From the above diagram, we can say that \mathrm{d}A = dzdy. \mathrm{d}A is simply the area of a very small region within the geometry. This allows us to re-write our formula for I_{yy}:

    \begin{equation*}I_{yy} = \int\int z^2 \mathrm{d}z\mathrm{d}y\end{equation*}

Considering the geometry of the section and the centroid, we can say that:

    \begin{equation*}I_{yy} = \int_{\frac{-b}{2}}^{\frac{b}{2}}\int_{\frac{-h}{2}}^{\frac{h}{2}} z^2 \mathrm{d}z\mathrm{d}y\end{equation*}

We evaluate the inner integral to give,

    \begin{equation*}I_{yy} = \int_{\frac{-b}{2}}^{\frac{b}{2}}\frac{1}{3}\frac{h^3}{4}\mathrm{d}y\end{equation}

And then the outer integral to give,

    \begin{equation*}I_{yy} = \frac{1}{3}\frac{h^3}{4}b = \frac{bh^3}{12}\end{equation*}

Moving onto I_{zz}:

    \begin{equation*}$I_{zz} = \int_{\frac{-b}{2}}^{\frac{b}{2}} \int_{\frac{-h}{2}}^{\frac{h}{2}} y^2 \mathrm{d}z\mathrm{d}y$\end{equation*}

We evaluate the inner integral to give,

    \begin{equation*}$I_{zz} = \int_{\frac{-b}{2}}^{\frac{b}{2}}h y^2 \mathrm{d}y $\end{equation*}

And then the outer integral to give,

    \begin{equation*}$I_{zz} = \frac{1}{3}\frac{b^3}{4}h = \frac{b^3h}{12}$\end{equation*}

And finally I_{yz}:

    \begin{equation*}$I_{yz} = \int_{\frac{-b}{2}}^{\frac{b}{2}}\int_{\frac{-h}{2}}^{\frac{h}{2}} yz \mathrm{d}z\mathrm{d}y$ = 0\end{equation}

Note: for a rectangle, I_{yz} is zero. This is an important result which we will discuss later.

6.0 Second Moment of Area of Simple Shapes

Whilst structural engineers will most commonly encounter sections comprising of multiple rectangular regions, occasionally, one might need to calculate the second moment of area of other common shapes.

The following standard results should be considered for reference:

Rectangle

img11_second_moment_of_area | DegreeTutors.com
Fig 11. Second moment of area of a rectangle.

    \begin{align*}I_{yy} &= \frac{bh^3}{12}\\\\I_{zz} &= \frac{b^3h}{12}\end{align*}

Circle

img12_second_moment_of_area | DegreeTutors.com
Fig 12. Second moment of area of a circle.

    \begin{equation*}I_{yy} = I_{zz} =\frac{\pi r^4}{4}\end{equation*}

Hollow circle or annulus

img13_second_moment_of_area | DegreeTutors.com
Fig 13. Second moment of area of a hollow circle or annulus.

    \begin{equation*}I_{yy} = I_{zz} =\frac{\pi (r_o^4-r_i^4)}{4}\end{equation*}

Triangle

img14_second_moment_of_area | DegreeTutors.com
Fig 14. Second moment of area of a triangle.

    \begin{align*}I_{yy} &= \frac{bh^3}{36}\\\\I_{zz}  &= \frac{b^3h}{36}\end{align*}

7.0 Parallel Axis Theorem

We now have a standard result for the second moment of area of a rectangle (and other simple geometries if we need them!) However, a solid rectangle is rarely used as a structural geometry so now we begin to think about the second moment of area of more complex geometries, and in particular those which are used in structural engineering applications.

So the question is, how can we use our understanding of the rectangle to define the second moment of area of more complex shapes?

The link is the parallel axis theorem. This theorem allows us to determine the second moment of area of a shape about an axis which is not a centroidal axis.

The parallel axis theorem states:

    \begin{equation*}I_{yy}^G = I_{yy}^L + Az^2\end{equation}

where I_{yy}^L is the local second moment of area about the region’s centroid, A is the area of the sub-geometry in question and z is the perpendicular distance between the global and local centroidal axes. Similarly,

    \begin{equation*}I_{zz}^G = I_{zz}^L + Ay^2\end{equation*}

where I_{zz}^L is the local second moment of area about the region’s centroid and the other symbols have the same meaning. And finally,

    \begin{equation*}I_{yz}^G = I_{yz}^L - Ayz\end{equation*}

where I_{yz}^L is the local product second moment of area about the region’s centroid, A is the area of the sub-geometry in question, y is the perpendicular distance between the global centroid and the local centroid and z is the perpendicular distance between the global centroid and the local centroid.

So, if we can split our complex overall geometry into a finite number of rectangular regions, we can use the parallel axis theorem to define I_{yy}^G, I_{zz}^G and I_{yz}^G i.e. the gross second moment of area about the entire geometries centroidal axis.

So, in general, we can say

I_{yy}^G = \sum I_{yy}^L + \sum Az^2

I_{zz}^G = \sum I_{zz}^L + \sum Ay^2

I_{yz}^G = \sum I_{yz}^L - \sum Ayz

The above formulae are simply saying that the gross second moment of area is the sum of the local second moments of area for each region in addition to a contribution coming from the position of the region relative to the overall geometries’ centroidal axes.

We now have all the necessary tools to complete three examples, so let’s make a start!

8.0 Example 1 – I Section

The first example is a 610 x 305 x 149 UKB section. So-called I-sections are typically used in applications which require resistance against large bending moments.

Note: when calculating the second moment of area by hand, we typically ignore the radius junction (root radius) between web and flanges. This simplification can be accepted as the impact on the overall second moment of area is minimal.

img15_second_moment_of_area | DegreeTutors.com
Fig 15. General arrangement of example 1.

Step 1: Calculate the Global Centroid

We start by calculating the centroid using the following formula,

    \begin{align*}y_c &= \frac{\sum Ay}{\sum A}\\\\z_c &= \frac{\sum Az}{\sum A}\end{align*}

Remember, we take our reference origin to be at the bottom left corner of the section.

img16_second_moment_of_area | DegreeTutors.com
Fig 16. Reference origin for example 1.

    \begin{equation*}Y_c = \frac{(\textrm{Area 1 x Distance (Y) from Centroid of 1 to Origin)+(Area 2 x Distance (Y) from Centroid of 2 to Origin)+ (Area 3 x Distance (Y) from Centroid of 3 to Origin)}}{\textrm{(Area 1)+(Area 2)+(Area 3)}}\end{equation*}

    \begin{equation*}Y_c = \frac{(304.8 \times 19.7 \times \frac{304.8}{2})+(304.8 \times 19.7 \times \frac{304.8}{2})+((612.4-2(19.7)) \times 11.8 \times \frac{304.8}{2})}{(304.8 \times 19.7)+(304.8 \times 19.7)+((612.4-2(19.7)) \times 11.8)} = 152.4mm\end{equation*}

    \begin{equation*}Z_c = \frac{(\textrm{Area 1 x Distance (Z) from Centroid of 1 to Origin)+(Area 2 x Distance (Z) from Centroid of 2 to Origin)+ (Area 3 x Distance (Z) from Centroid of 3 to Origin)}}{\textrm{(Area 1)+(Area 2)+(Area 3)}}\end{equation*}

    \begin{equation*}Z_c = \frac{(304.8\times19.7\times(612.4-\frac{19.7}{2})+(304.8\times19.7\times\frac{19.7}{2})+((612.4-2(19.7))\times11.8\times\frac{612.4}{2})}{(304.8\times19.7)+(304.8\times19.7)+((612.4-2(19.7))\times11.8)} = 306.2mm\end{equation*}

img17_second_moment_of_area | DegreeTutors.com
Fig 17. Centroid of example 1.

We should always reflect on whether our calculation results are realistic. So, is our calculated centroid in agreement with what we would expect?

Hint: We have a doubly-symmetric section so we would expect the centroid to be at the midpoint (\frac{b}{2}, \frac{h}{2})…which is exactly what we have shown using our general formula for the centroid.

Now we can proceed to determine I_{yy}^G, I_{zz}^G and I_{yz}^G

Step 2: Calculate I_{yy}^G

    \begin{equation*}$I_{yy}^G = \sum I_{yy}^L + \sum Az^2$\end{equation*}

    \begin{equation*}$I_{yy}^G = (I_{yy^1} + I_{yy^2} + I_{yy^3} + ((\textrm{Area 1} \times \Delta \textrm{z}^2)+(\textrm{Area 2} \times \Delta \textrm{z}^2) + (\textrm{Area 3} \times \Delta \textrm{z}^2))\end{equation*}

    \begin{align*}I_{yy}^G &= \left(\frac{304.8\times19.7^3}{12} + \frac{304.8\times19.7^3}{12} + \frac{11.8\times((612.4-2(19.7))^3}{12} \right) + \\&+ \left(304.8\times19.7 \times (-306.2+602.55)^2\right)  +\\&+ \left( 304.8\times19.7 \times (-306.2+0.5\times 19.7)^2\right) + \\&+ \left(11.8\times (612.4-2(19.7)) \times 0^2 \right)\end{align}

    \begin{equation*}I_{yy}^G = 1240066179 \textrm{ mm}^4\end{equation*}

Step 3: Calculate I_{zz}^G

    \begin{equation*}I_{zz}^G = \sum I_{zz}^L + \sum Ay^2\end{equation*}

    \begin{equation*}I_{zz}^G = (I_{zz^1} + I_{zz^2} + I_{zz^3} + ((\textrm{Area 1} \times \Delta \textrm{y}^2)+(\textrm{Area 2} \times \Delta \textrm{y}^2) + (\textrm{Area 3} \times \Delta \textrm{y}^2))\end{equation*}

    \begin{align*}I_{zz}^G &= \left(\frac{19.7\times304.8^3}{12} + \frac{19.7\times304.8^3}{12} + \frac{((612.4-2(19.7))\times 11.8^3}{12}\right) + \\&+ \left(304.8\times19.7 \times 0^2\right) + \\&+ \left(304.8\times19.7 \times 0^2\right) + \\&+ \left(11.8\times (612.4-2(19.7)) \times 0^2\right)\end{align*}

    \begin{equation*}I_{zz}^G = 93052101.09 \textrm{ mm}^4\end{equation*}

Step 4: Calculate I_{yz}^G

    \begin{equation*}I_{yz}^G = \sum I_{yz}^L - \sum Ayz\end{equation*}

    \begin{equation*}I_{yz}^G = (I_{yz^1} + I_{yz^2} + I_{yz^3} - (((\textrm{Area 1} \times \Delta \textrm{y} \times \Delta \textrm{z} )+(\textrm{Area 2} \times \Delta \textrm{y} \times \Delta \textrm{z}) + (\textrm{Area 3} \times \Delta \textrm{y}\times \Delta \textrm{z})))\end{equation*}

    \begin{align*}I_{yz}^G =& - (304.8\times19.7 \times (304.6 -0.5\times19.7) \times 0 ) +\\& - (304.8\times19.7 \times (0.5\times19.7-304.6) \times 0) +\\& - (11.8\times (612.4-2\times19.7) \times 0 \times 0)\end{align*}

    \begin{equation*}I_{yz}^G = 0\end{equation*}

So, to summarise the first example:

    \begin{align*}I_{yy}^G &= 1240066179 \textrm{ mm}^4 \approx 1.240 \times 10^9 \textrm{ mm}^4\\I_{zz}^G &= 93052101.09 \textrm{ mm}^4 \approx 9.305 \times 10^7 \textrm{ mm}^4\\I_{yz}^G &= 0\end{align*}

Let’s think about what these results mean. We have shown that depending on the orientation of the beam, there is a two-order-of-magnitude difference in the second moment of area (I_{yy}^G compared to I_{zz}^G). If there is a significant difference in the second moment of area, then, by inspection, the bending stiffness of a beam depends on the orientation in which it is used.

Clearly, therefore, if we want to maximise the bending stiffness, we should utilise the largest second moment of area of the section, I_{yy}^G, and therefore apply load to the beam parallel to the Z-axis.

9.0 Example 2 – Channel Section

The second example is a 150 x 90 x 10 unequal angle section. Unequal angles are often used as structural members under axial force but can also be used to resist small bending moments.

img18_second_moment_of_area | DegreeTutors.com
Fig 18. General arrangement of example 2.

Step 1: Calculate the Global Centroid

Remember, we take our reference origin to be at the bottom left corner of the section.

img19_second_moment_of_area | DegreeTutors.com
Fig 19. Reference origin for example 2.

    \begin{equation*}Y_c = \frac{(\textrm{Area 1 x Distance (Y) from Centroid of 1 to Origin)+(Area 2 x Distance (Y) from Centroid of 2 to Origin)}}{\textrm{(Area 1)+(Area 2)}}\end{equation*}

    \begin{equation*}Y_c = \frac{150\times 10 \times \frac{10}{2})+(80 \times 10 \times 50))}{(150 \times 10)+(80 \times 10)} = 20.65mm\end{equation*}

    \begin{equation*}Z_c = \frac{(\textrm{Area 1 x Distance (Z) from Centroid of 1 to Origin)+(Area 2 x Distance (Z) from Centroid of 2 to Origin)}}{\textrm{(Area 1)+(Area 2)}}\end{equation*}

    \begin{equation*}Z_c = \frac{(150\times 10 \times \frac{150}{2})+(80 \times 10 \times \frac{10}{2})}{(150 \times 10)+(80 \times 10)} = 50.65mm\end{equation*}

img20_second_moment_of_area | DegreeTutors.com
Fig 20. Centroid for example 2.

Is this what we would expect?

This example is less intuitive. We have a section which is asymmetric across both the Y and Z axes. It would have been difficult to determine the centroid without using our formulae.

Now we can proceed to determine I_{yy}^G, I_{zz}^G and I_{yz}^G.

Step 2: Calculate I_{yy}^G

    \begin{equation*}I_{yy}^G = \sum I_{yy}^L + \sum Az^2\end{equation*}

    \begin{equation*}I_{yy}^G = (I_{yy^1} + I_{yy^2} + ((\textrm{Area 1} \times \Delta \textrm{z}^2)+(\textrm{Area 2} \times \Delta \textrm{z}^2))\end{equation*}

    \begin{align*}I_{yy}^G &= \frac{10\times150^3}{12} + \frac{80\times10^3}{12} + \\& + (150\times 10 \times (-50.65+75)^2) + \\& + (80\times10 \times (-50.65+5)^2) \end{align*}

    \begin{equation*}I_{yy}^G = 5375688 \textrm{ mm}^4\end{equation*}

Step 3: Calculate I_{zz}^G

    \begin{equation*}I_{zz}^G = \sum I_{zz}^L + \sum Ay'^2\end{equation*}

    \begin{equation*}I_{zz}^G = (I_{zz^1} + I_{zz^2} + ((\textrm{Area 1} \times \Delta \textrm{y}^2)+(\textrm{Area 2} \times \Delta \textrm{y}^2))\end{equation*}

    \begin{align*}I_{zz}^G &= \frac{150\times10^3}{12} + \frac{10\times80^3}{12} + \\& + (150\times 10 \times (-20.65+5)^2) + \\& + (80\times10 \times (-20.65+50)^2) \end{align*}

    \begin{equation*}I_{zz}^G = 1495688\textrm{ mm}^4\end{equation*}

Step 4: Calculate I_{yz}^G

    \begin{equation*}I_{yz}^G = \sum I_{yz}^L - \sum Ayz\end{equation*}

    \begin{equation*}I_{yz}^G = I_{yz^1} + I_{yz^2} - ((\textrm{Area 1} \times \Delta \textrm{y} \times \Delta \textrm{z} )+(\textrm{Area 2} \times \Delta \textrm{y} \times \Delta \textrm{z}))\end{equation*}

    \begin{equation*}I_{yz}^G = - (150\times10 \times (-20.65+5)(-50.65+75))-(80\times10 \times (-20.65+50)(-50.65+5))\end{equation*}

    \begin{equation*}I_{yz}^G = 1643478\textrm{ mm}^4\end{equation*}

To summarise:

    \begin{align*}I_{yy}^G &= 5375688 \textrm{ mm}^4 \approx 5.376 \times 10^6 \textrm{ mm}^4\\I_{zz}^G &= 1495688\textrm{ mm}^4 \approx 1.496 \times 10^6 \textrm{ mm}^4\\I_{yz}^G &= 1643478\textrm{ mm}^4 \approx 1.643 \times 10^6 \textrm{ mm}^4\end{align*}

We have shown that depending on the orientation of the beam, there is a small difference in the second moment of area (I_{yy^G} compared to I_{zz^G}). The difference is much smaller than what we observed for the I-section. This makes sense as we have already mentioned that channel sections aren’t typically used within applications which require large bending capacity.

Interestingly, unlike the I-section, the product second moment of area (I_{yz^G}) does not equal 0. We will discuss the importance of this in the final sections of this article.

10.0 Example 3 – Z Section

The third example is a Z-section. Z-sections are commonly used on pitched roof structures to span between rafters and support roofing materials.

img21_second_moment_of_area | DegreeTutors.com
Fig 21. General arrangement of example 3.

Step 1: Calculate the Global Centroid

Remember, we take our origin to be at the bottom left corner of the section.

img22_second_moment_of_area | DegreeTutors.com
Fig 22. Reference origin for example 3.

    \begin{equation*}Z_c = \frac{(\textrm{Area 1 x Distance (Z) from Centroid of 1 to Origin)+(Area 2 x Distance (Z) from Centroid of 2 to Origin)+ (Area 3 x Distance (Z) from Centroid of 3 to Origin)}}{\textrm{(Area 1)+(Area 2)+(Area 3)}}\end{equation*}

    \begin{equation*}Z_c = \frac{(100x10x(200-5))+(100x10x5)+(180x10x100)}{(100x10)+(100x10)+(180x10)} = 100mm\end{equation*}

    \begin{equation*}Y_c = \frac{(\textrm{Area 1 x Distance (Y) from Centroid of 1 to Origin)+(Area 2 x Distance (Y) from Centroid of 2 to Origin)+ (Area 3 x Distance (Y) from Centroid of 3 to Origin)}}{\textrm{(Area 1)+(Area 2)+(Area 3)}}\end{equation*}

    \begin{equation*}Y_c = \frac{(100\times10\times50)+(100\times10\times140)+(180\times10\times95)}{(100\times10)+(100\times10)+(180\times10)} = 95mm\end{equation*}

img23_second_moment_of_area | DegreeTutors.com
Fig 23. Centroid for example 3.

Is this what we would expect? We have a symmetric section (symmetric about an inclined plane) and so we would expect the centroid to be at the origin of a rotated axis, which is exactly what we have shown using the formula.

Now we can proceed to determine I_{yy}^G, I_{zz}^G and I_{yz}^G.

Step 2: Calculate I_{yy}^G

    \begin{equation*}$I_{yy}^G = \sum I_{yy}^L + \sum Az^2$\end{equation*}

    \begin{equation*}I_{yy}^G = (I_{yy^1} + I_{yy^2} + I_{yy^3} + ((\textrm{Area 1} \times \Delta \textrm{z}^2)+(\textrm{Area 2} \times \Delta \textrm{z}^2) + (\textrm{Area 3} \times \Delta \textrm{z}^2))\end{equation*}

    \begin{align*}I_{yy}^G &= \left( \frac{100\times10^3}{12} + \frac{100\times10^3}{12} + \frac{10\times180^3}{12} \right) +\\&+ (100\times10 \times (-100+195)^2) +\\&+(100\times10 \times (-100+5)^2) +\\&+(180\times 10 \times (-100+100)^2)\end{align*}

    \begin{equation*}I_{yy}^G = 22926666 \textrm{ mm}^4\end{equation*}

Step 3: Calculate I_{zz}^G

    \begin{equation*}I_{zz}^G = \sum I_{zz}^L + \sum Ay^2\end{equation*}

    \begin{equation*}I_{zz}^G = (I_{zz^1} + I_{zz^2} + I_{zz^3} + ((\textrm{Area 1} \times \Delta \textrm{y}^2)+(\textrm{Area 2} \times \Delta \textrm{y}^2) + (\textrm{Area 3} \times \Delta \textrm{y}^2))\end{equation*}

    \begin{align*}I_{zz}^G &= \left(\frac{10\times100^3}{12} + \frac{10\times100^3}{12} + \frac{180\times10^3}{12}\right) + \\&+(100\times10 \times (-95+140)^2) +\\&+(100\times10 \times (-95+50)^2) +\\&+(180\times 10 \times (-95+95)^2)\end{align*}

    \begin{equation*}I_{zz}^G = 5731666 \textrm{ mm}^4\end{equation*}

Step 4: Calculate I_{yz}^G

    \begin{equation*}I_{yz}^G = \sum I_{yz}^L - \sum Ayz\end{equation*}

    \begin{equation*}I_{yz}^G = (I_{yz^1} + I_{yz^2} + I_{yz^3} - (((\textrm{Area 1} \times \Delta \textrm{y} \times \Delta \textrm{z} )+(\textrm{Area 2} \times \Delta \textrm{y} \times \Delta \textrm{z}) + (\textrm{Area 3} \times \Delta \textrm{y}\times \Delta \textrm{z})))\end{equation*}

    \begin{align*}I_{yz}^G &= -(100\times10 \times (-95+140)(-100+195)) + \\&- (100\times10 \times (-95+50)(-100+5)) + \\&-(180\times 10 \times (-95+95)(-100+100))\end{align*}

    \begin{equation*}I_{yz}^G = -8550000 \textrm{ mm}^4\end{equation*}

To summarise:

    \begin{align*}I_{yy}^G &= 22926666 \textrm{ mm}^4 \approx 2.293 \times 10^7 \textrm{ mm}^4\\I_{zz}^G &= 5731666 \textrm{ mm}^4 \approx 5.732 \times 10^6 \textrm{ mm}^4\\I_{yz}^G &= -8550000 \textrm{ mm}^4 \approx 8.55 \times 10^6 \textrm{ mm}^4\end{align*}

We have shown that depending on the orientation of the beam, there is a an order of magnitude difference in the second moment of area (I_{yy}^G compared to I_{zz}^G).

Interestingly, unlike the I-section the product second moment of area (I_{yz}^G) does not equal 0. We will develop our intuitive understanding of these results in the final section of this article.

11.0 Understanding the product second moment of area (I_{yz}^G)

At the start of this tutorial, we established an axes convention. In the three examples we have covered, the second moment of area has always been calculated with respect to the two axes (y and z). Let’s think about our results from examples 1, 2 and 3.

In example 1 (I-Section), we obtained I_{yz}^G = 0, in example 2 (Channel section), we obtained I_{yz}^G> 0 and in example 3 (Z section), we obtained I_{yz}^G < 0.

Thinking about the three examples we have covered, are the values of I_{yy}^G and I_{zz}^G the maximum possible values which we could obtain for that specific section?

It is not immediately obvious, but if we were to invoke a new axis system (u-v) which is slightly rotated compared to the y-z axes, could we obtain a higher value of I_{uu}^G and I_{vv}^G compared to our existing axes of I_{yy}^G and I_{zz}^G?

These new axes are called the principal bending axes because they represent the orientation of the section, which gives the maximum and minimum values for the second moment of area.

The magnitude of the product second moment of area (I_{yz}^G) indicates to us how ‘far away’ the assumed y-z axes are away from the principal bending axes (u-v). Let us revisit our three examples to understand how this works in practice.

In example 1 (I-Section), we obtained I_{yz}^G = 0; therefore, the y-z axes are coincident with the principal bending axes (u-v) and I_{yy}^G = I_{uu}^G and I_{zz}^G = I_{vv}^G.

In example 2 (Channel section), we obtained I_{yz}^G > 0 and therefore the y-z axes is not the principal bending axes and the u-v axes will yield I_{uu}^G > I_{yy}^G and I_{vv}^G < I_{zz}^G. The positive sign of I_{yz}^G indicates the direction which we need to rotate the y-z axes to obtain the principal bending axis (clockwise).

In example 3 (Z section), we obtained I_{yz}^G < 0, and therefore the y-z axes is not the principal bending axes and the u-v axes will yield I_{uu}^G > I_{yy}^G and I_{vv}^G < I_{zz}^G. The negative sign of I_{yz}^G indicates the direction which we need to rotate the y-z axes to obtain the principal bending axis (anti-clockwise).

12.0 Obtaining the orientation of the principal bending axes

As we have discussed, for sections where I_{yz}^G\neq 0, the principal bending axes (u-v) are not coincident with our assumed y-z axes.

It is important, therefore, that we can determine the orientation of the principal bending axes (u-v) so that we can identify the orientation of the structural member which will yield the largest bending stiffness (i.e. utilising the largest possible second moment of area).

We can do this using three simple formulae:

    \begin{equation*}I_{max} = I_{uu}^G = \frac{I_{yy}^G+I_{zz}^G}{2} + \sqrt{\frac{I_{yy}^G-I_{zz}^G}{2}^2+{I_{yz}^G}^2}\end{equation*}

    \begin{equation*}I_{min} = I_{vv}^G = \frac{I_{yy}^G+I_{zz}^G}{2} + \sqrt{\frac{I_{yy}^G-I_{zz}^G}{2}^2+{I_{yz}^G}^2}\end{equation*}

    \begin{equation*}\alpha= \tan^{-1}(\frac{{I_{yz}^G}}{I_{max}-I_{zz}^G})\end{equation*}

You may recognise these formula from a previous DegreeTutors tutorial on Mohr’s circle.  In fact, the same theory developed in that tutorial on stress analysis, maps directly across to this discussion of second moments of area. In particular the concept of principal stresses and principal axes is a direct analogue to the principal bending axis introduced here. 

Let’s revisit example 2 to understand these formula, recall:

    \begin{align*}I_{yy}^G &= 5375688 \textrm{ mm}^4 \approx 5.376 \times 10^6 \textrm{ mm}^4\\I_{zz}^G & = 1495688\textrm{ mm}^4 \approx 1.496 \times 10^6 \textrm{ mm}^4\\I_{yz}^G &= 1643478\textrm{ mm}^4\approx 1.643 \times 10^6 \textrm{ mm}^4\end{align*}

Evaluating I_{max}:

    \begin{align*}I_{max} &=\frac{5375688 + 1495688}{2} + \sqrt{\left(\frac{5375688 - 1495688}{2}\right)^2+1643478^2}\\&=5978249\textrm{ mm}^4\end{align*}

Similarly, I_{min} is,

    \begin{align*}I_{min} &= \frac{5375688 + 1495688}{2} - \sqrt{\left(\frac{5375688 - 1495688}{2}\right)^2+1643478^2}\\&=893126\textrm{ mm}^4\end{align*}

And finally, \alpha is:

    \begin{align*}\alpha &= \tan^{-1}(\frac{1643478}{5978249-1495688})\\&=20.13^\circ\end{align*}

We can re-draw example 2, showing the original axes system (y-z) and our new principal bending axes (u-v)

img24_second_moment_of_area | DegreeTutors.com
Fig 24. Superimposition of the principal (u-v) axes over the original y-z axes.

So, if we wanted to orientate the channel section in order to utilise the section’s maximum second moment of area, then we should apply transverse load parallel to the V axes (i.e. creating bending about the U axis).

13.0 Wrapping up

Within this tutorial, we have developed methods to calculate the centroid of a structural section and the three descriptions of the second moment of area for a section.

We have also developed a method for calculating the principal bending axes for a section which helped us recognise that certain sections are more suitable in certain applications than others. All structural sections have a specific orientation which will yield the maximum and minimum bending stiffness.

You might be asking yourself if all this calculation work is a necessary skill for a structural engineer? Well, for the majority of cases, the centroid and second moments of area of common structural sections are available from industry data tables (and these are always helpful for checking hand calculations!).

However, complex or bespoke sections will not be covered in data tables, and therefore it is essential that a structural engineer can calculate these fundamental properties. For example, steel plate girders are typically required when an off-the-shelf I-section isn’t suitable. A structural engineer would need to calculate the second moment of area for such a bespoke section.

Thats it for this tutorial. I hope you have found the content useful and informative – if you did, tell someone else who might also find it helpful!

Author

Callum Wilson

BSc (Hons), BEng (Hons)

Hi, I’m Callum, a contributor on DegreeTutors.com. I hope you have found this content interesting and useful. I am excited to be producing content for DegreeTutors.com as I enjoy the challenge of communicating technical content in a simple and accessible format. Feel free to get in touch or follow DegreeTutors on any of the social accounts.