Updated 15 April 2024
Reading time: 41 mins

Steel Column Design using Eurocode 3 - A Complete Guide

Learn how to design steel columns subject to compression and biaxial bending.
[object Object]
by Callum Wilson
Free PDF - Download the PDF version of this step-by-step-guide

Free PDF - Download the PDF version of this step-by-step-guide

Steel Design to Eurocode 3 - Part 2 of 4

Welcome to this, the second part in our series exploring structural steel design to Eurocode 3. The tutorial is focussed on steel column design.

In the previous tutorial on steel beam design, we introduced the Eurocode 3 document and in particular, Eurocode 3 Part 1-1: General rules and rules for buildings. We won't repeat this here, so refer back to part 1 now if you haven’t read it yet.

This is another quite long and detailed tutorial, so again, my advice is to tackle it in multiple sittings! It’s broken up as follows:

✅ Section 1: Eurocode 3 and Compression Members

In this first section, we’ll discuss how Eurocode 3 classifies compression members based on the actions that they must resist. This is critical as different classifications attract different design checks.

✅ Section 2: Steel Column Design Workflow

In section 2 we explore the column design workflow, using our what, where, why, how model. We’ll explain each step in the design workflow making sure to spend time understanding the motivations behind each design check as well as where in the code the design check is specified. This is quite a lengthy and detailed section but it will serve as a helpful reference, particularly as we move onto our worked example in the next section.

✅ Section 3: Steel Column Design - Worked Example

In this section, we work through a complete column design example. The design example has been constructed to attract as many of the design checks as possible in order to demonstrate the full breadth of the design process.

Essentially, we’ll be tackling a complex bi-axial bending design scenario. Once you’ve completed this example, 90% of practical column designs will look quick and easy by comparison!

✅ Section 4: Initial Sizing of Steel Columns

Luckily not all column designs are as involved as the example we considered in the previous section. So, it will be helpful to understand how the process is simplified for more straightforward design scenarios. We’ll cover this in section 4, giving you a more commonly used design flow for more typical column designs.

✅ Section 5: Wrapping up and Coming up next

In the final section, we briefly reflect on what we’ve achieved and briefly tee-up the next instalment of this steel design tutorial series!

1.0 Eurocode 3 and Compression Members

As we know, Eurocode 3 categorises design checks based on the internal forces developed in the member. Steel members which primarily resist axial compression are called columns. They represent one of the major structural components used across the world.

One of the most important ideas to grasp when considering the design of columns to Eurocode 3 is that although columns do predominantly resist axial compression, they may also resist bending moments - to a greater or lesser extent, depending on their configuration within the wider structure.

Eurocode 3 sets out a series of design checks which validate both the cross-sectional resistance and buckling resistance of compression members. In this sense, the design of columns and beams to Eurocode 3 is a similar process.

As is true for most Eurocode designs, we consider both the ultimate limit state (ULS) and serviceability limit states (SLS). However, for columns, we are almost exclusively concerned with the ULS case. Why is this? By the end of the article, you will be able to answer this question!

ℹ️ Note

A steel column is one type of compression member. Other types of compression members do exist, for example, struts within braced frames or the top chord of a sagging truss. This article will only cover the design of steel columns but the majority of the design checks can be applied to other, more general forms of compression member.

1.1 Types of Column

Recall that we reached the conclusion that steel columns are a form of compression member which predominantly resists axial compressive forces. Prior to undertaking any Eurocode 3 design work, it’s critical that we recognise three distinct types of column behaviour. Recognising which behaviour our column most closely represents will be extremely useful when it comes to completing a Eurocode 3-compliant design.

Note that this article is concerned with the design of structural steel columns. We will not discuss the theoretical background to column buckling. This is discussed in our column buckling series (part 1, part 2 and part 3).

1.1.1 Simple columns

As the name would suggest, simple columns represent the simplest form of column behaviour. So what makes a column simple?

In fact a simple column is any column which forms part of a simple structure. A simple structure is composed of members connected by pinned joints and which resists horizontal loads by way of bracing elements rather than via frame action.

The definition of a simple column is as follows…

'a simple column is a compression member which is designed for the prevailing and predominant axial force but also for a nominal bending moment which is induced in the column due to the eccentricity of the transferred force between beam and column at connection points.'

Importantly, the eccentricity, ee, of the vertical load is typically no larger than 100mm100\:mm and therefore the design bending moment acting upon the column is typically not large.

1.1.2 Uniaxial beam-columns

Simple structures are extremely popular in modern construction as they permit efficient design and fast construction. In many cases, bracing elements can provide horizontal stability, making simple construction a much-favoured choice.

However, there are structures where bracing elements are not suitable. For example, consider a typical portal-framed structure. When lateral loads are applied in the plane of the portal frames, bending moments that develop at the column-rafter connections are the mechanism by which lateral stability is achieved. Bracing is typically provided to resist lateral loads which act out of plane. We discuss lateral stability further in this tutorial.

From this observation, we can deduce that in response to in-plane loading, the portal columns will experience axial force and uniaxial bending i.e. bending about an axis perpendicular to the span of the rafters (typically the major principal axis for the column section).

This leads us to the notion of a 'beam-column'. There is no set magnitude of bending moment at which a simple column becomes a beam-column; this is left to the Engineer's judgment.

Put simply, if your column is experiencing significant bending moments about a single axis in addition to axial force, then it should be classified as a uniaxial beam-column.

1.1.3 Biaxial beam-columns

As you might have guessed, a biaxial beam-column is a further progression from the uniaxial beam-column. A biaxial beam-column experiences axial force and significant bending moments about both of its principal axes.

Remember: for both the uniaxial and biaxial beam-column cases, a suitable connection must be detailed between beam and column which permits the transfer of bending moments between these elements. We won’t be covering connection design in this tutorial.

Initially, it can be confusing to categorise a column into one of the above three categories. This tutorial will provide you with all of the tools needed to assess each type of column. Fortunately, if you wrongly classify a simple column as a bi-axial beam-column, the design checks will actually result in a more efficient design than if you had jumped straight into the simplifications afforded by a classification as a simple column.

2.0 Steel Column Design Workflow

The following section will identify and explain the design steps which are necessary for steel column design to Eurocode 3. Just like in the previous steel beam design tutorial, the purpose of this section is to highlight

  • what the checks are,
  • how to conduct the checks
  • and most importantly, why we are conducting these checks.

In the next section, we will apply all of these design checks to steel column design example.

📖 Have you got a copy of Eurocode 3?

The following discussion assumes that you have access to a copy of Eurocode 3. We'll reference various clauses - these, unless otherwise stated, refer to Eurocode 3. Unfortunately, it's not possible to reproduce the various extracts from Eurocode 3 within this article.

If you need to purchase a copy of Eurocode 3, you can do so from the British Standards Institute (BSI). A good alternative to purchasing the complete Eurocode is to purchase Structural Eurocodes: Extracts from the Structural Eurocodes for Students of Structural Design which contains the relevant sections of Eurocode 3 as well as other Eurocodes.

Structural Eurocodes: Extracts from the Structural Eurocodes for Students of Structural Design | EngineeringSkills.com

Structural Eurocodes: Extracts from the Structural Eurocodes for Students of Structural Design published by the British Standards Institute.

To help structure the discussion and aid your understanding, we’ll adopt a What, Why, Where and How template:

  • What - A summary of the design check.

  • Where - The relevant clauses within Eurocode 3.

  • Why - Why are we conducting this check?

  • How - A summary of how to conduct the check.

Preliminary step: General description for members under combined bending and axial compression

A steel compression member, depending on whether it falls into the simple column, uniaxial beam-column or biaxial beam-column classification, is subject to a different set of design checks within Eurocode 3. The design checks are split into cross-section resistance checks and member buckling checks.

Helpfully for us, Eurocode 3 provides interaction formulae which are used to assess both cross-sectional resistance and buckling resistance under the combined actions of axial compression and bending.

That said, understanding the nuances of these interaction checks is not so straightforward so we are deliberately starting this design workflow by considering the overarching interaction formulae. It is important that we recognise the 'big picture' before jumping into the detailed design steps.

Note: steel columns are typically class 1 or class 2 sections - (think about how bulky they are!). The following design workflow assumes a class 1 or class 2 section.

What?

For cross-sectional resistance:

[My,EdMN,y,Rd]α+[Mz,EdMN,z,Rd]β1\Big[\frac{M_{y,Ed}}{M_{N,y,Rd}}\Big]^\alpha + \Big[\frac{M_{z,Ed}}{M_{N,z,Rd}}\Big]^\beta \leq 1

where:

MN,RdM_{N,Rd} is the design plastic moment of resistance when reduced due to the axial force NEdN_{Ed}.

For member buckling

For any steel member subject to combined compression and bending, the following two expressions must be satisfied to satisfy member buckling criteria.

NEdχyNRkγM1+kyyMy,Ed+ΔMy,EdχLTMy,RkγM1+kyzMz,Ed+ΔMz,EdMz,RkγM11.0\frac{N_{Ed}}{\chi_y\frac{N_{Rk}}{\gamma_{M1}}} + k_{yy}\frac{M_{y,Ed}+\Delta M_{y,Ed}}{\chi_{LT}\frac{M_{y,Rk}}{\gamma_{M1}}} + k_{yz}\frac{M_{z,Ed}+\Delta M_{z,Ed}}{\frac{M_{z,Rk}}{\gamma_{M1}}} \leq 1.0
NEdχzNRkγM1+kzyMy,Ed+ΔMy,EdχLTMy,RkγM1+kzzMz,Ed+ΔMz,EdMz,RkγM11.0\frac{N_{Ed}}{\chi_z\frac{N_{Rk}}{\gamma_{M1}}} + k_{zy}\frac{M_{y,Ed}+\Delta M_{y,Ed}}{\chi_{LT}\frac{M_{y,Rk}}{\gamma_{M1}}} + k_{zz}\frac{M_{z,Ed}+\Delta M_{z,Ed}}{\frac{M_{z,Rk}}{\gamma_{M1}}} \leq 1.0

where:

NEd,My,Ed and Mz,EdN_{Ed}, M_{y,Ed} \textrm{ and } M_{z,Ed} are the design values of the compression force and the maximum design moments about the respective axes.

ΔMy,Ed and ΔMz,Ed\Delta M_{y,Ed} \textrm{ and } \Delta M_{z,Ed} are the moments induced in the section which occur due to the shift of the centroidal axis for class 4 sections.

χy and χz\chi_y \textrm{ and } \chi_z are the reduction factors which are required to describe flexural buckling (column buckling).

χLT\chi_{LT} is the reduction factor which is required to assess lateral torsional buckling.

kyy,kyz,kzy and kzzk_{yy}, k_{yz}, k_{zy} \textrm{ and } k_{zz} are interaction factors - more on these later!

Where?

  • For cross-sectional resistance: EN 1993-1-1: 2005 Clause 6.2.9
  • For member buckling: EN 1993-1-1: 2005 Clause 6.3.3

Why?

For cross-sectional resistance

If we think about an elastic stress distribution across a steel section subject to combined axial compression and bending we will recognise that the maximum compressive stress will occur at the extreme fibre of the section where the compressive stresses due to the axial force and the compressive stresses due to the bending moments are combined.

We know from elastic theory that first yield develops when the stress at any point reaches the yield strength of the material, fyf_y. Put simply, first yield can be reached sooner in sections resisting combined axial compression and bending and we must therefore account for this during design.

💡 Remember:

Earlier, we said we are considering class 1 or class 2 sections, which can, of course, develop full plastic capacity. This means that the first yielding (as described above) does not constitute section failure.

For member buckling

Looking at the two interaction formulae above we can easily see the relevance of the expressions by breaking down the constituent parts.

Without delving into the intricacies of the interaction factors kijk_{ij} just yet, we can say that a section under combined axial compression and bending has adequate buckling capacity if:

  • the ratio of design axial force to the flexural buckling resistance 1\leq 1, and

  • the ratio of design bending moment about the major bending axis to the lateral torsional buckling resistance 1\leq 1, and

  • the ratio of design bending moment about the minor bending axis to the moment of resistance about the minor bending axis 1\leq 1, and

  • the sum of the three above ratios 1\leq 1

Note: we repeat the above checks for the flexural buckling resistance in both axes.

How?

Throughout the rest of this design workflow, we are going to break down each design step. In essence, we need to obtain the values of all of the constants and variables in the above expressions...so let’s get started!

Step 1: Cross-section classification

What?

In step 1, we need to classify the geometry of the member cross-section. In Eurocode 3, four classifications are defined. The classification is based on the ability of a member to develop its full plastic section capacity. If you have completed the previous tutorial in this series, then these descriptions will be familiar:

  • Class 1 cross-sections are those that can form a plastic hinge with the rotation capacity required by plastic analysis without reducing the resistance.

  • Class 2 cross-sections are those which can develop their plastic moment resistance, but have limited rotation capacity because of local buckling.

  • Class 3 cross-sections are those in which the stress in the extreme compression fibre of the steel member, assuming an elastic distribution of stresses, can reach the yield strength, but local buckling is liable to prevent development of the plastic moment resistance.

  • Class 4 cross-sections are those in which local buckling will occur before the yield stress is reached in one or more parts of the cross-section.

Where?

  • For the definition of the classes: EN 1993-1-1: 2005 Clause 5.5.2
  • For the classification of internal compression parts: EN 1993-1-1: 2005 Table 5.2 (sheet 1 of 3)
  • For classification of outstand flanges: EN 1993-1-1: 2005 Table 5.2 (sheet 2 of 3)

Why?

Classifying our cross-section is the method chosen within Eurocode 3 to simplify the analysis of different types of members.

Cross-sections of class 3 and class 4 are susceptible to local buckling failure prior to yielding, whereas cross-sections of class 1 and class 2 are most likely to fail due to the lack of cross-section resistance, as opposed to local buckling.

This simple cross-section classification step allows us to predict the member’s behaviour, and complete the relevant design checks based on the thickness of its constituent parts.

How?

  1. Using the diagrams at the top of Table 5.2 (sheet 2 of 3), identify the outstand flange compression parts. For an I-section, this will be the flanges of the member.

  2. Calculate the value of cc, making sure to account for the thickness of surrounding elements and the root radius of any joints.

  3. From the section geometry drawing, determine the thickness of the flanges. This is typically tft_f for an I-section.

  4. Now calculate the value of the quotient ct\frac{c}{t} for the critical outstand flange.

  5. Next, determine the value of ϵ\epsilon by evaluating the quotient 235fy\sqrt{\frac{235}{f_y}} where fyf_y is the yield strength determined in step 1.

  6. We are considering a 'part' in pure compression and hence we need to use the left-most column of Table 5.2 (sheet 2 of 3).

  7. Now classify the internal compression part of the section by checking if ct\frac{c}{t} is less than 9ϵ9\epsilon (Class 1), 10ϵ10\epsilon (Class 2) or 14ϵ14\epsilon (Class 3).

  8. Now we move on to classify the internal compression parts. For an I-section this is typically the web of the section.

  9. Calculate the value of cc, making sure to account for the thickness of surrounding elements and the root radius of any joints.

  10. From the section geometry drawing, determine the thickness of the internal compression parts. Again, this is typically twt_w for an I-section.

  11. Now calculate the value of the quotient ct\frac{c}{t} for the internal compression member.

  12. We are considering a 'part' in combined bending and compression and hence we need to use the right-most column of Table 5.2 (sheet 1 of 3).

  13. To select our limiting values for part classification, we need to calculate α\alpha, compression proportion of the web.

  14. We define α\alpha as:

α=(h2(tf+r)+12NEdtwfy)cw\alpha = \frac{\Big(\frac{h}{2}- (t_f + r) + \frac{1}{2}\frac{N_{Ed}}{t_wf_y}\Big)}{c_w}
  1. Using our value of α\alpha we can compare our value for ct\frac{c}{t} with the limiting values.

  2. Once you have classified both the internal compression parts (the web) and the outstand flanges (the compression flange) then we are ready to classify the section.

  3. The overall classification of the section is the lowest of the classifications of the individual parts. For example, if the internal compression part classification is class 3 and the outstand flange classification is class 1, then the section is class 3.

  4. You might have guessed this, but if the cross-section parts do not meet the criteria for class 3 then they are class 4. Class 4 steel members are subject to an entire Eurocode 3 document - Eurocode 3: Design of steel structures - Part 1-5: General rules - Plated structural elements. For the time being, we will assume that our section is either class 1, class 2 or class 3.

ℹ️ Note

If you are designing a compression member which does not experience any bending moment, then during classification we would use the table columns from Table 5.2 which correspond to pure compression only. Other than this slight difference, the method is the same.

Step 2: Calculate NRdN_{Rd}

What?

Prior to calculating the cross-sectional resistance under the combined action of bending moment and axial force, it is mandatory to first calculate NRdN_{Rd}. NRdN_{Rd} is the design resistance of the section against axial compressive force - initially as a standalone action.

We don't completely ignore bending moments when calculating NRdN_{Rd} since we have assumed a bending stress distribution within the internal compression part (the web) during cross-section classification.

For class 1, 2 or 3 cross-sections, NRdN_{Rd} is given by the following:

NRd=AfyγM0N_{Rd} = \frac{Af_y}{\gamma_{M0}}

where AA is the area of the cross-section.

Where?

EN 1993-1-1: 2005 Clause 6.2.4

Why?

This is perhaps the simplest of all the Eurocode 3 design checks. We calculate the cross-section resistance against axial compression by multiplying the cross-sectional area of the section by the yield strength of the steel. Simple.

How?

This simple check can be completed by following the simple expression shown above. The expression is a little more complicated for class 4 cross-sections but for the time being, we won't concern ourselves with those.

Step 3: Combined cross-section resistance - conservative method

What?

There are a few different approaches which can be taken to assess the cross-sectional resistance under combined bending and compression. We will first consider the conservative approach.

This check simply states that the linear summation of the utilisation ratios of the stress resultants corresponding to the combination of design actions NEd,My,EdN_{Ed}, M_{y,Ed} and Mz,EdM_{z,Ed} must be less than unity. The expression we refer to is as follows:

NEdNRd+My,EdMy,Rd+Mz,EdMz,Rd1\frac{N_{Ed}}{N_{Rd}}+\frac{M_{y,Ed}}{M_{y,Rd}}+\frac{M_{z,Ed}}{M_{z,Rd}} \leq 1

where:

NRd,My,RdN_{Rd}, M_{y,Rd} and Mz,RdM_{z,Rd} are the design values of the resistance depending on cross-sectional classification and including any reduction that may be caused by shear effects.

Where?

EN 1993-1-1: 2005 Clause 6.2.1 (7) and Expression 6.2

Why?

This check does exactly what it says on the tin. We are checking that the utilisation ratios for the three leading internal design actions (axial force, major bending and minor bending) are less than unity and that the sum is less than unity. As stated, this is a conservative approach - after all, think about summing three numbers to be less than 1.

How?

  • For My,RdM_{y,Rd} and Mz,RdM_{z,Rd}, refer to the previous tutorial in this series.
  • For NRdN_{Rd}, refer to step 3 in this design workflow.

Step 4: Combined cross-section resistance - alternative method

An alternative method (to that covered in step 3) is provided within Eurocode 3. Typically this alternative method results in a more efficient design.

What?

This design check offers a more rigorous approach to assess the cross-sectional resistance against combined bending moments and axial force.

Where?

  • Design basis: EN 1993-1-1: 2005 Clause 6.2.9.1
  • Interaction criteria for uniaxial bending: EN 1993-1-1: 2005 Clause 6.2.9.1 Expression 6.31
  • Interaction criteria for biaxial bending: EN 1993-1-1: 2005 Clause 6.2.9.1 Expression 6.41

Why?

This alternative method is used for the same reason as the method presented in step 3 i.e. ensuring that there is adequate cross-sectional resistance against the combined bending moment and axial compression acting on the section.

How?

For My,RdM_{y,Rd} and Mz,RdM_{z,Rd}, refer to the previous tutorial. Once these values are defined, proceed to determine the reduced bending resistance (due to axial force), MN,y,RdM_{N,y,Rd} and MN,z,RdM_{N,z,Rd} using EN 1993-1-1: 2005 Clause 6.2.9.1

For NRdN_{Rd}, refer to step 2 in this design workflow.

Step 5: Compression member buckling resistance

What?

Now we move into the more challenging aspect of the Eurocode 3 design checks for beam-columns - member buckling resistance. Once again, we begin with the buckling resistance of the pure compression member. In this design check, we will not be considering the design bending moments - except of course that we maintain the cross-section classification that we obtained in step 1.

Where?

EN 1993-1-1: 2005 Clause 6.3.1.1

Why?

In step 3, we calculated the cross-sectional resistance of a pure compression member. However, as we discussed at the start of this tutorial, cross-sectional resistance is just one part of steel member design to Eurocode 3. We must also consider failure due to member buckling - which can occur at an axial compression force which is much less than the cross-sectional resistance, NRdN_{Rd} (again, refer to our tutorial series on theoretical column buckling behaviour or Euler buckling, referenced above).

How?

This is a multi-step design process, so we will break down the compression member buckling resistance check into multiple smaller steps. However, the final expression which we need to evaluate is as follows (for a class 1, 2 or 3 cross-section).

Nb,Rd=χAfyγM1N_{b,Rd} =\frac{\chi Af_y}{\gamma_{M1}}

where χ\chi is the reduction factor for the relevant buckling mode.

Step 5.1 Calculate the non-dimensional slenderness λˉ\bar{\lambda}

The non-dimensional slenderness is an indicator of the tendency of a compression member to fail due to flexural buckling. Column sections have a major and a minor buckling axis. As a result, we define two non-dimensional slenderness values, λyˉ\bar{\lambda_y} and λzˉ\bar{\lambda_z}.

In general, for a class 1, 2 or 3 cross-section, λˉ\bar{\lambda} is given by:

λˉ=AfyNcr\bar{\lambda} = \frac{Af_y}{N_{cr}}

where NcrN_{cr} is the elastic critical force for the relevant buckling mode based on the gross cross-sectional properties. An alternative way to express λˉ\bar{\lambda} is as follows.

λˉ=AfyNcr=Lcri1λ1\bar{\lambda} = \frac{Af_y}{N_{cr}} = \frac{L_{cr}}{i}\frac{1}{\lambda_1}

where:

LcrL_{cr} is the critical buckling length i.e. the distance between points of full flexural buckling restraint (for a given axis) and ii is the radius of gyration - taken about the same axis as LcrL_{cr} and λ1\lambda_1 is given by,

λ1=93.9235fy\lambda_{1} = 93.9\sqrt{\frac{235}{f_y}}

So, at the end of this process, we will have two values, λyˉ\bar{\lambda_y} and λzˉ\bar{\lambda_z}.

EN 1993-1-1: 2005 Clause 6.3.1.2

Step 5.2 Select the buckling curve and find α\alpha

Now we use Table 6.2 to determine the correct buckling curve. These buckling curves allow us to determine the imperfection factor aa which we must account for during our flexural buckling check.

Once we know the correct buckling curve, we use Table 6.1 to determine α\alpha.

EN 1993-1-1: 2005 Table 6.1 and Table 6.2

Step 5.3 Calculate the reduction factor χ\chi

Now we are at the point where we can calculate the reduction factor χ\chi. All of the sub-steps we have completed in step 5 have been leading us towards calculating this reduction factor.

First, we calculate Φ\Phi as follows:

Φ=0.5[1+α(λˉLT0.2)+λˉ2]\Phi = 0.5\Big[1+\alpha\Big(\bar{\lambda}_{LT}-0.2\Big)+\bar{\lambda}^2\Big]

And using Φ\Phi we can now calculate χ\chi as follows:

χ=1Φ+Φ2λˉ21.0\chi = \frac{1}{\Phi+\sqrt{\Phi^2-\bar{\lambda}^2}} \leq 1.0

Remember that we must complete this calculation for both axes of the section. So by the end of this step you should have two values, χy\chi_y and χz\chi_z.

EN 1993-1-1: 2005 Clause 6.3.1.2 Expression 6.49

Step 5.4 Calculate the buckling resistance for the beam Nb,RdN_{b,Rd}

We have now calculated our reduction factors χy\chi_y and χz\chi_z which leaves us to just determine the buckling resistance of the compression member. At this point, we need to choose the critical value for χ\chi i.e. the largest value, since this will produce the lowest buckling resistance.

Nb,Rd=χAfyγM1N_{b,Rd} = \chi A\frac{f_y}{\gamma_{M1}}

where γM1\gamma_{M1} is taken to be 1.001.00 in accordance with 6.1 NOTE 2B

EN 1993-1-1: 2005 Clause 6.3.1.1 Expression 6.47

Step 6: Combined buckling resistance

What?

We are now coming towards the end of the Eurocode 3 design checks for a compression member under combined bending and axial force. As we stated at the start of this design workflow, we need to validate the buckling resistance of the member under the action of the combined internal actions.

Recall the expressions:

NEdχyNRkγM1+kyyMy,Ed+ΔMy,EdχLTMy,RkγM1+kyzMz,Ed+ΔMz,EdMz,RkγM11.0\frac{N_{Ed}}{\chi_y\frac{N_{Rk}}{\gamma_{M1}}} + k_{yy}\frac{M_{y,Ed}+\Delta M_{y,Ed}}{\chi_{LT}\frac{M_{y,Rk}}{\gamma_{M1}}} + k_{yz}\frac{M_{z,Ed}+\Delta M_{z,Ed}}{\frac{M_{z,Rk}}{\gamma_{M1}}} \leq 1.0
NEdχzNRkγM1+kzyMy,Ed+ΔMy,EdχLTMy,RkγM1+kzzMz,Ed+ΔMz,EdMz,RkγM11.0\frac{N_{Ed}}{\chi_z\frac{N_{Rk}}{\gamma_{M1}}} + k_{zy}\frac{M_{y,Ed}+\Delta M_{y,Ed}}{\chi_{LT}\frac{M_{y,Rk}}{\gamma_{M1}}} + k_{zz}\frac{M_{z,Ed}+\Delta M_{z,Ed}}{\frac{M_{z,Rk}}{\gamma_{M1}}} \leq 1.0

Which we can simplify as a result of our section being class 1, 2 or 3 down to:

NEdχyNRkγM1+kyyMy,EdχLTMy,RkγM1+kyzMz,EdMz,RkγM11.0\frac{N_{Ed}}{\chi_y\frac{N_{Rk}}{\gamma_{M1}}} + k_{yy}\frac{M_{y,Ed}}{\chi_{LT}\frac{M_{y,Rk}}{\gamma_{M1}}} + k_{yz}\frac{M_{z,Ed}}{\frac{M_{z,Rk}}{\gamma_{M1}}} \leq 1.0
NEdχzNRkγM1+kzyMy,EdχLTMy,RkγM1+kzzMz,EdMz,RkγM11.0\frac{N_{Ed}}{\chi_z\frac{N_{Rk}}{\gamma_{M1}}} + k_{zy}\frac{M_{y,Ed}}{\chi_{LT}\frac{M_{y,Rk}}{\gamma_{M1}}} + k_{zz}\frac{M_{z,Ed}}{\frac{M_{z,Rk}}{\gamma_{M1}}} \leq 1.0

Note: The final component of the two expressions deals with the minor axis bending resistance for the member. You may have noted that on the denominator of the term, we do not apply the reduction factor χLT\chi_{LT}. This is because members do not exhibit lateral torsional buckling due to bending about the minor axis.

Where?

  • For the general expression: EN 1993-1-1: 2005 Clause 6.3.3 Expression 6.61 and 6.62
  • _For the kijk_{ij} factors_: EN 1993-1-1: 2005 Annex A or Annex B - depending on which method is chosen. In this tutorial we will choose Annex A method.

Why?

At this point, it should be clear that the buckling resistance of a member will typically govern the design. In the previous tutorial on steel beam design, we calculated the buckling resistance of members subject to bending moments alone and in step 4 of this article we calculated the buckling resistance of a compression member subject to axial compression only. It will come as no surprise then that the final step of the process is to check the buckling resistance under the combination of both bending and compression.

How?

We are going to use the more exact approach covered in Annex A of Eurocode 3 as it leads to more efficient designs, which in turn, leads to material and cost savings! It should be noted that the simplified approach covered in Annex B is more conservative (and hence still safe to use).

Once again, we are going to break down this Eurocode 3 design check into smaller and more manageable steps.

Step 6.1 Calculate the factor bLTb_{LT}

bLT=0.5aLTλˉ02My,EdχLTMpl,y,RdMz,EdMpl,z,Rdb_{LT} = 0.5a_{LT}\bar{\lambda}_0^2\frac{M_{y,Ed}}{\chi_{LT}M_{pl,y,Rd}}\frac{M_{z,Ed}}{M_{pl,z,Rd}}

where

aLT=1ITIy0a_{LT} = 1 - \frac{I_T}{I_y} \geq 0

where,

ITI_T is the St.Venant torsion constant (available in the section tables), IyI_y is the second moment of area about the major yyy-y axis and λˉ0\bar{\lambda}_0 is the non-dimensional slenderness for lateral-torsional buckling due to uniform bending moment i.e. when ψy=1\psi_y = 1

Note: head back to step 8.3 and 8.4 in the previous tutorial and calculate λˉ\bar{\lambda} for the case when ψy=1\psi_y = 1. Once calculated, take λˉ=λˉ0\bar{\lambda} = \bar{\lambda}_0.

EN 1993-1-1: 2005 Annex A Table A.1

Step 6.2 Calculate the factor cLTc_{LT}

cLT=10aLTλˉ025+λzˉ4My,EdCmyχLTMpl,y,Rdc_{LT} = 10a_{LT}\frac{\bar{\lambda}_0^2}{5+\bar{\lambda_z}^4}\frac{M_{y,Ed}}{C_{my}\chi_{LT}M_{pl,y,Rd}}

where

Cmy=Cmy,0+(1Cmy,0)ϵyaLT1+ϵyaLTC_{my} = C_{my,0} + \big(1-C_{my,0}\Big)\frac{\sqrt{\epsilon_y}a_{LT}}{1+\sqrt{\epsilon_y}a_{LT}}

using:

(for class 1, 2 and 3 cross sections)

ϵy=My,EdNEdAWel,y\epsilon_y = \frac{M_{y,Ed}}{N_{Ed}}\frac{A}{W_{el,y}}

and

Cmi,0C_{mi,0} is obtained from Table A.2 using knowledge of the bending moment diagram.

Note:

If

λˉ00.2C1(1NEdNcr,z)(1NEdNcr,TF)4\bar{\lambda}_0 \leq 0.2\sqrt{C_1}\sqrt[4]{\Big(1-\frac{N_{Ed}}{N_{cr,z}}\Big)\Big(1-\frac{N_{Ed}}{N_{cr,TF}}\Big)}

then

Cmy=Cmy,0C_{my} = C_{my,0}

EN 1993-1-1: 2005 Annex A Table A.1

Step 6.3 Calculate the factor dLTd_{LT}

dLT=2aLTλˉ00.1+λz4ˉMy,EdCmyχLTMpl,y,RdMz,EdCmzMpl,z,Rdd_{LT} = 2a_{LT}\frac{\bar{\lambda}_0}{0.1+\bar{\lambda_z^4}}\frac{M_{y,Ed}}{C_{my}\chi_{LT}M_{pl,y,Rd}}\frac{M_{z,Ed}}{C_{mz}M_{pl,z,Rd}}

where

Cmz=Cmz,0C_{mz} = C_{mz,0}

with Cmz,0C_{mz,0} coming from Annex A Table A.2 using knowledge of the bending moment diagram for the given axis.

EN 1993-1-1: 2005 Annex A Table A.1

Step 6.4 Calculate the factor eLTe_{LT}

eLT=1.7aLTλˉ00.1+λz4ˉMy,EdCmyχLTMpl,y,Rde_{LT} = 1.7a_{LT}\frac{\bar{\lambda}_0}{0.1+\bar{\lambda_z^4}}\frac{M_{y,Ed}}{C_{my}\chi_{LT}M_{pl,y,Rd}}

EN 1993-1-1: 2005 Annex A Table A.1

Step 6.5 Calculate the factors Cyy,Cyz,Czy,CzzC_{yy}, C_{yz}, C_{zy}, C_{zz}

Cyy=1+(wy1)[(21.6wyCmy2λmaxˉ1.6wyCmy2λmaxˉ2)nplbLT]Wel,yWpl,yC_{yy} = 1+\Big(w_y-1\Big)\Big[\Big(2-\frac{1.6}{w_y}C_{my}^2\bar{\lambda_{max}}-\frac{1.6}{w_y}C_{my}^2\bar{\lambda_{max}}^2\Big)n_{pl}-b_{LT}\Big]\geq \frac{W_{el,y}}{W_{pl,y}}
Cyz=1+(wz1)[(214Cmz2λmaxˉ2wz5)nplcLT]0.6wzwyWel,zWpl,zC_{yz} = 1+\Big(w_z-1\Big)\Big[\Big(2-14\frac{C_{mz}^2\bar{\lambda_{max}}^2}{w_z^5}\Big)n_{pl}-c_{LT}\Big]\geq 0.6\sqrt{\frac{w_z}{w_y}}\frac{W_{el,z}}{W_{pl,z}}
Czy=1+(wy1)[(214Cmy2λmaxˉ2wy5)npldLT]0.6wywzWel,yWpl,yC_{zy} = 1+\Big(w_y-1\Big)\Big[\Big(2-14\frac{C_{my}^2\bar{\lambda_{max}}^2}{w_y^5}\Big)n_{pl}-d_{LT}\Big]\geq 0.6\sqrt{\frac{w_y}{w_z}}\frac{W_{el,y}}{W_{pl,y}}
Czz=1+(wz1)[(21.6wzCmz2λmaxˉ1.6wzCmz2λmaxˉ2)npleLT]Wel,zWpl,zC_{zz} = 1+\Big(w_z-1\Big)\Big[\Big(2-\frac{1.6}{w_z}C_{mz}^2\bar{\lambda_{max}}-\frac{1.6}{w_z}C_{mz}^2\bar{\lambda_{max}}^2\Big)n_{pl}-e_{LT}\Big]\geq \frac{W_{el,z}}{W_{pl,z}}

where

npl=NEdNRk/γ1n_{pl} = \frac{N_{Ed}}{N_{Rk}/\gamma_1}
wy=Wpl,yWel,y1.5w_y = \frac{W_{pl,y}}{W_{el,y}}\leq 1.5
wz=Wpl,zWel,z1.5w_z = \frac{W_{pl,z}}{W_{el,z}}\leq 1.5
λmaxˉ=max(λyˉ,λzˉ)\bar{\lambda_{max}} = \textrm{max}\Big(\bar{\lambda_y},\bar{\lambda_z})

EN 1993-1-1: 2005 Annex A Table A.1

Step 6.6 Calculate the interaction factors Kyy,Kyz,Kzy,KzzK_{yy}, K_{yz}, K_{zy}, K_{zz}

We will make the assumption that our section is class 1 or class 2:

kyy=CmyCmLTμy1NEdNcr,y1Cyyk_{yy} = C_{my}C_{mLT}\frac{\mu_y}{1-\frac{N_{Ed}}{N_{cr,y}}}\frac{1}{C_{yy}}
kyz=Cmzμy1NEdNcr,z1Cyz0.6wzwyk_{yz} = C_{mz}\frac{\mu_y}{1-\frac{N_{Ed}}{N_{cr,z}}}\frac{1}{C_{yz}}0.6\sqrt{\frac{w_z}{w_y}}

where

μy=1NEdNcr,y1χyNEdNcr,y\mu_y = \frac{1-\frac{N_{Ed}}{N_{cr,y}}}{1-\chi_y\frac{N_{Ed}}{N_{cr,y}}}

and

kzy=CmyCmLTμz1NEdNcr,y1Czy0.6wywzk_{zy} = C_{my}C_{mLT}\frac{\mu_z}{1-\frac{N_{Ed}}{N_{cr,y}}}\frac{1}{C_{zy}}0.6\sqrt{\frac{w_y}{w_z}}
kzz=Cmzμz1NEdNcr,z1Czzk_{zz} = C_{mz}\frac{\mu_z}{1-\frac{N_{Ed}}{N_{cr,z}}}\frac{1}{C_{zz}}

where

μz=1NEdNcr,z1χzNEdNcr,z\mu_z = \frac{1-\frac{N_{Ed}}{N_{cr,z}}}{1-\chi_z\frac{N_{Ed}}{N_{cr,z}}}

Note:

CmLT=1.0 ifλˉ00.2C1(1NEdNcr,z)(1NEdNcr,TF)4C_{mLT} = 1.0\textrm{ if} \bar{\lambda}_0 \leq 0.2\sqrt{C_1}\sqrt[4]{\Big(1-\frac{N_{Ed}}{N_{cr,z}}\Big)\Big(1-\frac{N_{Ed}}{N_{cr,TF}}\Big)}
CmLT=Cmy2aLT(1NEdNcr,z)(1NEdNcr,T) ifλˉ0>0.2C1(1NEdNcr,z)(1NEdNcr,TF)4C_{mLT} = C_{my}^2\frac{a_{LT}}{\sqrt{\Big(1-\frac{N_{Ed}}{N_{cr,z}}\Big)\Big(1-\frac{N_{Ed}}{N_{cr,T}}\Big)}}\textrm{ if} \bar{\lambda}_0 > 0.2\sqrt{C_1}\sqrt[4]{\Big(1-\frac{N_{Ed}}{N_{cr,z}}\Big)\Big(1-\frac{N_{Ed}}{N_{cr,TF}}\Big)}

EN 1993-1-1: 2005 Annex A Table A.1

Step 6.7 Evaluate the combined buckling resistance inequalities

All that remains is to plug in all of the calculated values which we have determined in the preceding design steps. If we can show that the value of both inequalities is less than unity then our member has sufficient buckling resistance against the combined bending moments and axial compression.

NEdχyNRkγM1+kyyMy,EdχLTMy,RkγM1+kyzMz,EdMz,RkγM11.0\frac{N_{Ed}}{\chi_y\frac{N_{Rk}}{\gamma_{M1}}} + k_{yy}\frac{M_{y,Ed}}{\chi_{LT}\frac{M_{y,Rk}}{\gamma_{M1}}} + k_{yz}\frac{M_{z,Ed}}{\frac{M_{z,Rk}}{\gamma_{M1}}} \leq 1.0
NEdχzNRkγM1+kzyMy,EdχLTMy,RkγM1+kzzMz,EdMz,RkγM11.0\frac{N_{Ed}}{\chi_z\frac{N_{Rk}}{\gamma_{M1}}} + k_{zy}\frac{M_{y,Ed}}{\chi_{LT}\frac{M_{y,Rk}}{\gamma_{M1}}} + k_{zz}\frac{M_{z,Ed}}{\frac{M_{z,Rk}}{\gamma_{M1}}} \leq 1.0

EN 1993-1-1: 2005 Clause 6.3.3 Expression 6.61 and 6.62

3.0 Steel Column Design - Worked Example

We are now ready to implement our Eurocode 3 design workflow. The worked example deliberately includes as many of the Eurocode 3 column design checks as possible.

Mercifully, most columns do not experience large bending moments, and consequently, many of the design steps will be omitted for simpler design scenarios. Recognising this, section 4.0 will cover a set of simple rules and design simplifications that can be used for the sizing of typical, simple columns.

In each section of the worked example, look out for the references back to the design workflow introduced in the previous section i.e. step 1, step 2, step 3 etc.

3.1 Setting up the problem

In our example problem, we are going to consider a compression member which is subject to bi-axial bending. The characteristics of the problem are as summarised below. Try to think about how these characteristics will inform our selection of Eurocode 3 design checks. The problem is presented with a horizontal member but imagine this member orientated vertically in a structure.

  • The member is subject to an axial force parallel to the xx-axis - in a structure this axial force would result from gravity loads.

  • The member is subject to a concentrated load acting in the major principal plane i.e. parallel to the zz-axis.

  • The member is subject to a uniformly distributed load (UDL) acting in the minor principal plane i.e. parallel to the yy-axis.

  • Lateral deflection and twisting are prevented at both ends of the member and at the midspan point. Other than these points, the member is free to deflect laterally in both planes and to rotate.

3.2 Loading

3.2.1 Permanent Actions (GG)

  • Uniformly distributed load, W1W_1 acting in the positive Y-direction = 2.3kN/m2.3\:kN/m
  • Concentrated Load at Point A in negative Z-direction, P1=5.5kNP_1 = 5.5\:kN
  • Concentrated Load at Point B in negative X-direction, P2=40kNP_2 = 40\:kN

3.2.2 Variable Actions (QQ)

  • Concentrated Load at Point A in negative Z-direction, P1=7.5kNP_1 = 7.5\:kN
  • Concentrated Load at Point B in negative X-direction, P2=80kNP_2 = 80\:kN

3.3 Load Cases

3.3.1 Partial factors for actions

  • Partial factor for permanent actions, γG\gamma_G = 1.35
  • Partial factor for variable actions, γG\gamma_G = 1.50

3.3.2 Combination

γGj,supGj,sup+γQ,1Q1\gamma_{Gj,sup}G_{j,sup} +\gamma_{Q,1}Q_1

3.4 Analysis

Now, we progress to the analysis stage, which will allow us to determine the internal actions (against which we will conduct our member design checks). You can use whichever method of analysis you like - we won’t dwell on the method of analysis here.

3.4.1 Applied loads

UDL (self-weight): W=γGGg=1.35×2.3=3.11W = \gamma_{G}G_g = 1.35 \times 2.3 = 3.11kN/m

Concentrated Load at A: P1=γGG1+γQQ1=(1.35×5.5)+(1.5×7.5)=18.68P_1 = \gamma_{G}G_1 + \gamma_Q Q_1 = \Big(1.35\times5.5\Big) + \Big(1.5 \times 7.5\Big) = 18.68kN

Concentrated Load at B: P2=γGG2+γQQ2=(1.35×40)+(1.5×80)=174P_2 = \gamma_{G}G_2 + \gamma_Q Q_2 = \Big(1.35\times40\Big) + \Big(1.5 \times 80\Big) = 174kN

General arrangement of worked example | EngineeringSkills.com

Fig 1. General arrangement of worked example column.

3.4.2 Axial Force

Clearly, the column is subject to an axial compression force of 174kN174\:kN.

3.4.3 Bending moments

As we are considering bi-axial bending, it is necessary to look at bending moment diagrams for both the major and minor axes.

Major axis For the major axis, a single point load at mid-span gives us a simple bending moment diagram with peak value of 42kNm42\:kNm at mid-span.

Major axis bending moment diagram | EngineeringSkills.com

Fig 2. Major axis bending moment diagram.

Minor axis For the minor axis, an additional support is present, which results in a two-span continuous member. The peak moment is a hogging moment located over the central support with a value of 7.87kNm7.87\:kNm.

Minor axis bending moment diagram | EngineeringSkills.com

Fig 3. Minor axis bending moment diagram.

3.4.4 Shear Forces

The shear force diagrams are included in this analysis section but we will not be referring to them any further within this tutorial. For beam-columns, the shear forces very rarely govern. In fact, EN 1993-1-1: 2005 Clause 6.2.10 tells us that unless VEd0.5Vpl,RdV_{Ed} \geq 0.5V_{pl,Rd} then no reduction to axial or bending resistance. You will find that this clause is met in most cases.

Major axis shear force diagram | EngineeringSkills.com

Fig 4. Major axis shear force diagram.

Minor axis shear force diagram | EngineeringSkills.com

Fig 5. Minor axis shear force diagram.

3.4.5 Deflection

We obtain a separate deflection plot for the major and minor axes. For vertical members such as columns, the lateral deflection is less commonly reported.

However, in Eurocode 3 beam-column design checks, the deflection is used within buckling resistance checks. For that reason, these deflection plots are generated using the ULS combination of loads, rather than the SLS combination.

Major axis deflection plot (ULS) | EngineeringSkills.com

Fig 6. Major axis deflection plot (ULS).

Minor axis deflection plot (ULS) | EngineeringSkills.com

Fig 7. Minor axis deflection plot (ULS).

3.5 Initial section choice

For this worked example, we will try an S355 305x165x46 UB. We’ll discuss some useful tips on the initial sizing of columns and beam-columns at the end of this tutorial. This section is chosen to be quite conservative so that we have some confidence that all design checks will be clearly passed.

From the section properties tables we acknowledge the following:

Section Properties
  • Depth, h=306.6mmh = 306.6 \:mm
  • Width, b=165.7mmb =165.7 \:mm
  • Web thickness, tw=6.7mmt_w =6.7 \:mm
  • Flange thickness, tf=11.8mmt_f=11.8 \:mm
  • Root radius, r=8.9mmr=8.9 \:mm
  • Depth between fillets, d=265.2mmd=265.2 \:mm
  • Second moment of area about y-y axis, Iyy=9900cm4I_{yy} =9900\: cm^4
  • Second moment of area about z-z axis, Izz=896cm4I_{zz} =896 \:cm^4
  • Warping constant, Iw=0.195cm6I_w=0.195 \:cm^6
  • Torsion constant, IT=22.2cm4I_T=22.2 \:cm^4
  • Radius of gyration about y-y axis, iy=13cmi_y=13 \:cm
  • Radius of gyration about z-z axis, iz=3.908cmi_z=3.908 \:cm
  • Plastic modulus y-y axis, Wpl,y=720cm3W_{pl,y} =720\: cm^3
  • Plastic modulus z-z axis, Wpl,z=166cm3W_{pl,z} =166\: cm^3
  • Elastic modulus y-y axis, Wel,y=646cm3W_{el,y} =646\: cm^3
  • Elastic modulus z-z axis, Wel,z=108cm3W_{el,z} =108\: cm^3
  • Area, A=58.7cm2A=58.7 \:cm^2
  • Modulus of elasticity, E=210000N/mm2E=210000\: N/mm^2
  • Shear Modulus, G=81000N/mm2G=81000 \:N/mm^2

Note: the units for each of the above properties are those which are typically used in industry. You will find that using these units is more convenient and avoids using extremely large or small numbers!

3.6 Cross-section classification

Step 1: Determine the yield strength, fyf_y

  • Our section is S355 steel
  • We have 1616mm t40\leq t \leq 40 mm

fy=355\therefore f_y = 355 N/mm2^2

EN 1993-1-1: 2005 Table 3.1

Classify the cross-section - Outstand of compression flange

cf=btw2r2=165.76.72(8.9)2=70.6 mmc_f = \frac{b-t_w-2r}{2} = \frac{165.7-6.7-2(8.9)}{2} = 70.6\textrm{ mm}
cftf=70.611.8=5.98\frac{c_f}{t_f} = \frac{70.6}{11.8} = 5.98
ϵ=235355=0.81\epsilon = \sqrt{\frac{235}{355}} = 0.81

The limiting value for a Class 1 section is 9ϵ=7.299\epsilon = 7.29. Evaluating the limiting value:

5.987.295.98 \leq 7.29

\therefore The flange in compression is Class 1.

EN 1993-1-1: 2005 Table 5.2 (sheet 2 of 3)

Classify the cross-section - Web (part subject to bending and compression)

cw=d=265.2 mmc_w = d = 265.2\textrm{ mm}
cwtw=265.26.7=39.58\frac{c_w}{t_w} = \frac{265.2}{6.7} = 39.58
ϵ=235355=0.81\epsilon = \sqrt{\frac{235}{355}} = 0.81

Now given that our part is in combined bending and compression, we first have to calculate α\alpha

α=(h2(tf+r)+12NEdtwfy)cw=(306.62(11.8+8.9)+12174×1036.7×355)265.2=0.63\alpha = \frac{\Big(\frac{h}{2}- (t_f + r) + \frac{1}{2}\frac{N_{Ed}}{t_wf_y}\Big)}{c_w} = \frac{\Big(\frac{306.6}{2}- (11.8 + 8.9) + \frac{1}{2}\frac{174 \times10^3}{6.7\times355}\Big)}{265.2} = 0.63

0.630.50.63 \geq 0.5

\therefore The limiting value for a class 1 internal compression part is:

cwtw396ϵ13α1\frac{c_w}{t_w} \leq \frac{396\epsilon}{13\alpha -1}

Evaluating the limiting value:

39.5844.61=396(0.81)13(0.63)139.58 \leq 44.61 = \frac{396(0.81)}{13(0.63) -1}

\therefore The section of web in combined bending and compression is Class 1.

EN 1993-1-1: 2005 Table 5.2 (sheet 1 of 3)

Overall classification Both the flange and the web are Class 1. Therefore the overall section classification is Class 1.

3.7 ULS Checks

3.7.1 Cross-section resistance

We have a case of bi-axial bending and hence the following criterion must be used.

[My,EdMN,y,Rd]α+[Mz,EdMN,z,Rd]β1\Big[\frac{M_{y,Ed}}{M_{N,y,Rd}}\Big]^\alpha + \Big[\frac{M_{z,Ed}}{M_{N,z,Rd}}\Big]^\beta \leq 1

where for an I/H section, α=2\alpha = 2 and β=5n\beta =5n where n=NEdNpl,Rdn = \frac{N_{Ed}}{N_{pl,Rd}}

So in reality, we have to calculate three terms before we can evaluate the criteria for cross-sectional resistance; Npl,RdN_{pl,Rd}, MN,y,RdM_{N,y,Rd} and MN,z,RdM_{N,z,Rd}.

Step 2 Calculate NRdN_{Rd}

For a class 1 section:

Nc,Rd=Npl,Rd=AfyγM0=58.7×102×3551=2083.85kNN_{c,Rd} = N_{pl,Rd} = \frac{Af_y}{\gamma_{M0}} = \frac{58.7\times10^2\times355}{1} = 2083.85\textrm{kN}

EN 1993-1-1: 2005 Clause 6.2.4

Step 4 Alternative method for cross-section resistance

To make our design checks simpler, we enlist two specific clauses from Eurocode 3.

EN 1993-1-1: 2005 Clause 6.2.9.1 (4) tells us that:

for doubly symmetrical I- and H-sections or other flanges sections, allowance need not be made for the effect of the axial force on the plastic resistance moment about the y-y axis when both the following criteria are satisfied:

NEd0.25Npl,RdN_{Ed} \leq 0.25N_{pl,Rd}

and

NEd0.5hwtwfyγM0N_{Ed} \leq \frac{0.5h_wt_wf_y}{\gamma_{M0}}

so we proceed to test these two criterion:

174520.96=0.25×2083.85kN174 \leq 520.96 = 0.25\times 2083.85\textrm{kN}

and

174315.39=0.5×265.2×6.7×3551174 \leq 315.39 = \frac{0.5\times 265.2 \times 6.7 \times 355}{1}

So both criteria are met, which means that we can say the following for the plastic moment resistance about the y-y axis.

Mn,y,Rd=Mpl,y,Rd=Wpl,yfyγM0=720×103×3551=255.6kN.mM_{n,y,Rd} = M_{pl,y,Rd} = \frac{W_{pl,y}f_y}{\gamma_{M0}} = \frac{720 \times 10^3\times355}{1} = 255.6\textrm{kN.m}

EN 1993-1-1: 2005 Clause 6.2.9.1 (4) tells us that:

for doubly symmetrical I- and H-sections, allowance need not be made for the effect of the axial force on the plastic resistance moment about the z-z axis when:

NEdhwtwfyγM0N_{Ed} \leq \frac{h_wt_wf_y}{\gamma_{M0}}

so we proceed to test this criterion:

174630.70=265.2×6.7×3551174 \leq 630.70 = \frac{265.2 \times 6.7 \times 355}{1}

The criterion is met, which means that we can say the following for the plastic moment resistance about the z-z axis.

Mn,z,Rd=Mpl,z,Rd=Wpl,zfyγM0=166×103×3551=58.93kN.mM_{n,z,Rd} = M_{pl,z,Rd} = \frac{W_{pl,z}f_y}{\gamma_{M0}} = \frac{166\times 10^3\times355}{1} = 58.93\textrm{kN.m}

Now we can re-visit the governing criterion:

[My,EdMN,y,Rd]α+[Mz,EdMN,z,Rd]β1\Big[\frac{M_{y,Ed}}{M_{N,y,Rd}}\Big]^\alpha + \Big[\frac{M_{z,Ed}}{M_{N,z,Rd}}\Big]^\beta \leq 1
[42.03255.6]2+[7.8758.93]0.42=0.461\Big[\frac{42.03}{255.6}\Big]^2 + \Big[\frac{7.87}{58.93}\Big]^{0.42} = 0.46 \leq 1

Therefore the cross-section is adequate when we consider the combined actions of axial compression and bi-axial bending.

3.7.2 Buckling resistance

Step 5: Compression member buckling resistance

Before we progress with determining the combined buckling resistance of the member, we first need to determine the buckling resistance of the member, assuming it sees only axial compression.

Step 5.1 Calculate the non-dimensional slenderness λˉ\bar{\lambda}

We define two non-dimensional slenderness values, λyˉ\bar{\lambda_y} and λzˉ\bar{\lambda_z}.

λyˉ=Lcr,yiy1λ1=9000130176.4=0.91\bar{\lambda_y} = \frac{L_{cr,y}}{i_y}\frac{1}{\lambda_1} = \frac{9000}{130}\frac{1}{76.4} = 0.91
λzˉ=Lcr,ziz1λ1=450039176.4=1.51\bar{\lambda_z} = \frac{L_{cr,z}}{i_z}\frac{1}{\lambda_1} = \frac{4500}{39}\frac{1}{76.4} = 1.51

where:

λ1=93.9235fy=76.4\lambda_{1} = 93.9\sqrt{\frac{235}{f_y}} = 76.4

EN 1993-1-1: 2005 Clause 6.3.1.2

Step 5.2 Select the buckling curve and find αLT\alpha_{LT}

For h/b<1.2h/b < 1.2 and tf40t_f \leq 40mm, we consider the following:

  • For buckling about the yyy-y axis, use buckling curve a, so take α=0.21\alpha = 0.21
  • For buckling about the zzz-z axis, use buckling curve b, so take α=0.34\alpha = 0.34

EN 1993-1-1: 2005 Table 6.1 and Table 6.2

Step 5.3 Calculate the reduction factor χ\chi

Now we are at the point where we can calculate the reduction factor χ\chi. Recall that we need to calculate this for buckling about both axes.

First, we calculate Φ\Phi as follows:

Φy=0.5[1+αy(λyˉ0.2)+λyˉ2]=0.5[1+0.21(0.910.2)+0.912]=0.99\Phi_y = 0.5\Big[1+\alpha_y\Big(\bar{\lambda_y}-0.2\Big)+\bar{\lambda_y}^2\Big] = 0.5\Big[1+0.21\Big(0.91-0.2\Big)+0.91^2\Big] = 0.99
Φz=0.5[1+α(λzˉ0.2)+λzˉ2]=0.5[1+0.34(1.510.2)+1.512]=1.86\Phi_z = 0.5\Big[1+\alpha\Big(\bar{\lambda_z}-0.2\Big)+\bar{\lambda_z}^2\Big] = 0.5\Big[1+0.34\Big(1.51-0.2\Big)+1.51^2\Big] = 1.86

And using Φ\Phi we can now calculate χ\chi as follows:

χy=1Φy+Φy2λyˉ2=10.99+0.9920.912=0.721.0\chi_y = \frac{1}{\Phi_y+\sqrt{\Phi_y^2-\bar{\lambda_y}^2}} = \frac{1}{0.99+\sqrt{0.99^2-{0.91^2}}} = 0.72 \leq 1.0
χz=1Φz+Φz2λzˉ2=11.86+1.8621.512=0.341.0\chi_z = \frac{1}{\Phi_z+\sqrt{\Phi_z^2-\bar{\lambda_z}^2}} = \frac{1}{1.86+\sqrt{1.86^2-{1.51^2}}} = 0.34 \leq 1.0

EN 1993-1-1: 2005 Clause 6.3.1.2 Expression 6.49

Step 5.4 Calculate the buckling resistance for the beam Nb,RdN_{b,Rd}
Nb,y,Rd=χyAfyγM1=0.72×58.7×102×3551.0=1500kNN_{b,y,Rd} = \chi_y A\frac{f_y}{\gamma_{M1}} = 0.72 \times 58.7 \times 10^2\times \frac{355}{1.0} = 1500\:kN
Nb,z,Rd=χzAfyγM1=0.34×58.7×102×3551.0=708kNN_{b,z,Rd} = \chi_z A\frac{f_y}{\gamma_{M1}} = 0.34 \times 58.7 \times 10^2\times \frac{355}{1.0} = 708\:kN

where γM1\gamma_{M1} is taken to be 1.001.00 in accordance with 6.1 NOTE 2B

So we can say that the buckling resistance, Nb,RdN_{b,Rd}, is 708kN>174kN708\:kN > 174\:kN therefore, for the case of axial compression only, the member has adequate buckling resistance.

EN 1993-1-1: 2005 Clause 6.3.1.1 Expression 6.47

Step 6: Combined buckling resistance

Step 6.1 Calculate the factor bLTb_{LT}
bLT=0.5aLTλˉ02My,EdχLTMpl,y,RdMz,EdMpl,z,Rdb_{LT} = 0.5a_{LT}\bar{\lambda}_0^2\frac{M_{y,Ed}}{\chi_{LT}M_{pl,y,Rd}}\frac{M_{z,Ed}}{M_{pl,z,Rd}}

We break this expression down and calculate the individual terms one by one.

aLT=1ITIy=122.29000=0.9970a_{LT} = 1 - \frac{I_T}{I_y} = 1 - \frac{22.2}{9000} = 0.997 \geq 0
My,Ed=42.03kNmM_{y,Ed} = 42.03 \:kNm
Mz,Ed=7.87kNmM_{z,Ed} = 7.87 \:kNm
Mpl,y,Rd=255.6kNmM_{pl,y,Rd} = 255.6 \:kNm
Mpl,z,Rd=58.93kNmM_{pl,z,Rd} = 58.93 \:kNm

To calculate λˉ0\bar{\lambda}_0 and χLT\chi_{LT} we need to utilise the design checks we covered in our steel beam design tutorial. If you need to refresh yourself on these design checks, review that tutorial first. In this example, we’ll move swiftly through these results.

λˉ0=WyfyMcr\bar{\lambda}_0 = \sqrt{\frac{W_yf_y}{M_{cr}}}
Mcr,0=1×π2×210000×896×104(4500)2(0.195×1012896×104)+((4500)2×81000×22.2×104π2×210000×896×104)×106M_{cr,0} = 1 \times \frac{\pi^{2}\times210000\times 896 \times 10^4}{\Big(4500\Big)^2}\sqrt{\Big(\frac{0.195\times10^{12}}{896 \times 10^4}\Big)+\Big(\frac{(4500)^2\times81000\times22.2\times10^4}{\pi^2\times210000\times896 \times 10^4}\Big)} \times 10^{-6}
Mcr,0=186.53 kN.mM_{cr,0} = 186.53\textrm{ kN.m}
λˉ02=720×103×355186.53×106=1.17\bar{\lambda}_0^2 = \sqrt{\frac{720\times10^3\times355}{186.53\times10^6}} = 1.17
Mcr,LT=1.348×π2×210000×896×104(4500)2(0.195×1012896×104)+((4500)2×81000×22.2×104π2×210000×896×104)×106M_{cr,LT} = 1.348 \times \frac{\pi^{2}\times210000\times 896 \times 10^4}{\Big(4500\Big)^2}\sqrt{\Big(\frac{0.195\times10^{12}}{896 \times 10^4}\Big)+\Big(\frac{(4500)^2\times81000\times22.2\times10^4}{\pi^2\times210000\times896 \times 10^4}\Big)} \times 10^{-6}
Mcr,LT=251.44 kN.mM_{cr,LT} = 251.44\textrm{ kN.m}
λˉLT2=720×103×355251.44×106=1.008\bar{\lambda}_{LT}^2 = \sqrt{\frac{720\times10^3\times355}{251.44\times10^6}} = 1.008

We have a section with h/b<2h/b <2 and tf40t_f \leq 40mm and so choose LTB buckling curve a - therefore αLT=0.21\alpha_{LT} = 0.21

ΦLT,1=0.5[1+αLT(λˉLT0.2)+λˉLT2]=0.5[1+0.21(1.0080.2)+1.0082]=1.09\Phi_{LT,1} = 0.5\Big[1+\alpha_{LT}\Big(\bar{\lambda}_{LT}-0.2\Big)+\bar{\lambda}_{LT}^2\Big] = 0.5\Big[1+0.21\Big(1.008-0.2\Big)+1.008^2\Big] = 1.09
χLT=1ΦLT+ΦLT2λˉLT2=11.09+1.0921.0082=0.66\chi_{LT} = \frac{1}{\Phi_{LT}+\sqrt{\Phi_{LT}^2-\bar{\lambda}_{LT}^2}} = \frac{1}{1.09+\sqrt{1.09^2-1.008^2}} = 0.66

We can now come back and calculate bLTb_{LT}

bLT=0.5aLTλˉ02My,EdχLTMpl,y,RdMz,EdMpl,z,Rd=0.5×0.997×1.172×42.030.66×255.6×7.8758.93=0.023\begin{align*} b_{LT} &= 0.5a_{LT}\bar{\lambda}_0^2\frac{M_{y,Ed}}{\chi_{LT}M_{pl,y,Rd}}\frac{M_{z,Ed}}{M_{pl,z,Rd}} \\\\ &= 0.5\times0.997\times1.17^2\times \frac{42.03}{0.66\times255.6}\times \frac{7.87}{58.93} \\\\ &= 0.023 \end{align*}
Step 6.2 Calculate the factor cLTc_{LT}
cLT=10aLTλˉ025+λzˉ4My,EdCmyχLTMpl,y,Rdc_{LT} = 10a_{LT}\frac{\bar{\lambda}_0^2}{5+\bar{\lambda_z}^4}\frac{M_{y,Ed}}{C_{my}\chi_{LT}M_{pl,y,Rd}}

Once again, we break this expression down and calculate the individual terms one by one.

aLT=0.997a_{LT} = 0.997
λˉ0=1.17\bar{\lambda}_0 = 1.17
λzˉ=1.51\bar{\lambda_z} = 1.51
χLT=0.66\chi_{LT} = 0.66
My,Ed=42.03kN.mM_{y,Ed} = 42.03 \textrm{kN.m}
Mpl,y,Rd=255.6kN.mM_{pl,y,Rd} = 255.6 \textrm{kN.m}

Now we need to determine CmyC_{my}. The first step in this is to evaluate the inequality as follows:

λˉ00.2C1(1NEdNcr,z)(1NEdNcr,TF)4=0.21.348(1174917)(1174917)4=0.19\bar{\lambda}_0 \geq 0.2\sqrt{C_1}\sqrt[4]{\Big(1-\frac{N_{Ed}}{N_{cr,z}}\Big)\Big(1-\frac{N_{Ed}}{N_{cr,TF}}\Big)} = 0.2\sqrt{1.348}\sqrt[4]{\Big(1-\frac{174}{917}\Big)\Big(1-\frac{174}{917}\Big)} = 0.19

therefore

Cmy=Cmy,0+(1Cmy,0)ϵyaLT1+ϵyaLT=1.00348+(11.00348)×2.19×0.9971+2.19×0.997=1.0014C_{my} = C_{my,0} + \big(1-C_{my,0}\Big)\frac{\sqrt{\epsilon_y}a_{LT}}{1+\sqrt{\epsilon_y}a_{LT}} = 1.00348 + \big(1-1.00348\Big)\times \frac{\sqrt{2.19}\times0.997}{1+\sqrt{2.19}\times0.997} = 1.0014

using:

ϵy=My,EdNEdAWel,y=42.03174×58.7×104646×106=2.19\epsilon_y = \frac{M_{y,Ed}}{N_{Ed}}\frac{A}{W_{el,y}} = \frac{42.03}{174}\times\frac{58.7\times10^{-4}}{646\times10^{-6}} = 2.19

and returning back to the original expression for cLTc_{LT}

cLT=10aLTλˉ025+λzˉ4My,EdCmyχLTMpl,y,Rd=10×0.997×1.1725+1.514×42.031.0014×0.66×255.6=0.333c_{LT} = 10a_{LT}\frac{\bar{\lambda}_0^2}{5+\bar{\lambda_z}^4}\frac{M_{y,Ed}}{C_{my}\chi_{LT}M_{pl,y,Rd}} = 10\times0.997\times\frac{1.17^2}{5+{1.51}^4}\times\frac{42.03}{1.0014\times0.66\times255.6} = 0.333
Step 6.3 Calculate the factor dLTd_{LT}
dLT=2aLTλˉ00.1+λz4ˉMy,EdCmyχLTMpl,y,RdMz,EdCmzMpl,z,Rdd_{LT} = 2a_{LT}\frac{\bar{\lambda}_0}{0.1+\bar{\lambda_z^4}}\frac{M_{y,Ed}}{C_{my}\chi_{LT}M_{pl,y,Rd}}\frac{M_{z,Ed}}{C_{mz}M_{pl,z,Rd}}

Again, let us state what we know first.

aLT=0.997a_{LT} = 0.997
λˉ0=1.17\bar{\lambda}_0 = 1.17
λzˉ=1.51\bar{\lambda_z} = 1.51
χLT=0.66\chi_{LT} = 0.66
My,Ed=42.03kN.mM_{y,Ed} = 42.03 \textrm{kN.m}
Mz,Ed=7.87kN.mM_{z,Ed} = 7.87 \textrm{kN.m}
Mpl,y,Rd=255.6kN.mM_{pl,y,Rd} = 255.6 \textrm{kN.m}
Mpl,z,Rd=58.93kN.mM_{pl,z,Rd} = 58.93 \textrm{kN.m}
Cmy=1.0014C_{my} = 1.0014

So all that remains is to determine CmzC_{mz}

where

Cmz=Cmz,0C_{mz} = C_{mz,0}

and

Cmz,0=1+(π2EIzδxL2Mi,Ed(x)1)×NEdNcr,i=1+(π2×210000×896×1043.11×45004185×210000×896×10445002×3.11×4.528×1061)×174708=0.86\begin{align*} C_{mz,0} &= 1+\Big(\frac{\pi^2EI_z|\delta_x|}{L^2|M_{i,Ed(x)}|}-1\Big)\times\frac{N_{Ed}}{N_{cr,i}} \\\\ &= 1+\Big(\frac{\pi^2\times 210000 \times 896 \times 10^4|\frac{3.11\times4500^4}{185\times210000\times896\times10^4}|}{4500^2\times|\frac{3.11\times4.5^2}{8}|\times10^6|}-1\Big)\times\frac{174}{708} = 0.86 \end{align*}

Now revisiting the original expression for dLTd_{LT}

dLT=2aLTλˉ00.1+λz4ˉMy,EdCmyχLTMpl,y,RdMz,EdCmzMpl,z,Rd=2×0.997×1.170.1+1.514×42.031.0014×0.66×255.6×7.870.86×58.93=0.017\begin{align*} d_{LT} &= 2a_{LT}\frac{\bar{\lambda}_0}{0.1+\bar{\lambda_z^4}}\frac{M_{y,Ed}}{C_{my}\chi_{LT}M_{pl,y,Rd}}\frac{M_{z,Ed}}{C_{mz}M_{pl,z,Rd}} \\\\ &= 2\times 0.997\times \frac{1.17}{0.1+1.51^4}\times \frac{42.03}{1.0014\times0.66\times 255.6}\times\frac{7.87}{0.86\times 58.93} = 0.017 \end{align*}
Step 6.4 Calculate the factor eLTe_{LT}
eLT=1.7aLTλˉ00.1+λz4ˉMy,EdCmyχLTMpl,y,Rd=1.7×0.997×1.170.1+1.514×42.031.0014×0.66×255.6=0.093\begin{align*} e_{LT} &= 1.7a_{LT}\frac{\bar{\lambda}_0}{0.1+\bar{\lambda_z^4}}\frac{M_{y,Ed}}{C_{my}\chi_{LT}M_{pl,y,Rd}} \\\\ &= 1.7\times 0.997 \times \frac{1.17}{0.1+1.51^4}\times\frac{42.03}{1.0014\times 0.66 \times 255.6} = 0.093 \end{align*}
Step 6.5 Calculate the factors Cyy,Cyz,Czy,CzzC_{yy}, C_{yz}, C_{zy}, C_{zz}

Evaluating CyyC_{yy}

Cyy=1+(wy1)[(21.6wyCmy2λmaxˉ1.6wyCmy2λmaxˉ2)nplbLT]Wel,yWpl,yC_{yy} = 1+\Big(w_y-1\Big)\Big[\Big(2-\frac{1.6}{w_y}C_{my}^2\bar{\lambda_{max}}-\frac{1.6}{w_y}C_{my}^2\bar{\lambda_{max}}^2\Big)n_{pl}-b_{LT}\Big] \geq \frac{W_{el,y}}{W_{pl,y}}
Cyy=1+(1.111)[(21.61.11×1.00142×1.511.61.11×1.00142×1.512)0.0830.023]=0.970.90\begin{align*} C_{yy} &= 1+\Big(1.11-1\Big)\Big[\Big(2-\frac{1.6}{1.11}\times1.0014^2\times 1.51 -\frac{1.6}{1.11}\times1.0014^2\times 1.51^2\Big)0.083-0.023\Big] \\\\ &= 0.97 \geq 0.90 \end{align*}

Evaluating CyzC_{yz}

Cyz=1+(wz1)[(214Cmz2λmaxˉ2wz5)nplcLT]0.6wzwyWel,zWpl,zC_{yz} = 1+\Big(w_z-1\Big)\Big[\Big(2-14\frac{C_{mz}^2\bar{\lambda_{max}}^2}{w_z^5}\Big)n_{pl}-c_{LT}\Big]\geq 0.6\sqrt{\frac{w_z}{w_y}}\frac{W_{el,z}}{W_{pl,z}}
Cyz=1+(wz1)[(214Cmz2λmaxˉ2wz5)nplcLT]=0.790.454\begin{align*} C_{yz} &= 1+\Big(w_z-1\Big)\Big[\Big(2-14\frac{C_{mz}^2\bar{\lambda_{max}}^2}{w_z^5}\Big)n_{pl}-c_{LT}\Big] \\\\ &= 0.79 \geq 0.454 \end{align*}

Evaluating CzyC_{zy}

Czy=1+(wy1)[(214Cmy2λmaxˉ2wy5)npldLT]0.6wywzWel,yWpl,yC_{zy} = 1+\Big(w_y-1\Big)\Big[\Big(2-14\frac{C_{my}^2\bar{\lambda_{max}}^2}{w_y^5}\Big)n_{pl}-d_{LT}\Big]\geq 0.6\sqrt{\frac{w_y}{w_z}}\frac{W_{el,y}}{W_{pl,y}}
Czy=1+(wy1)[(214Cmy2λmaxˉ2wy5)npldLT]=0.8450.463\begin{align*} C_{zy} &= 1+\Big(w_y-1\Big)\Big[\Big(2-14\frac{C_{my}^2\bar{\lambda_{max}}^2}{w_y^5}\Big)n_{pl}-d_{LT}\Big]\\\\ &= 0.845 \geq 0.463 \end{align*}

Evaluating CzzC_{zz}

Czz=1+(wz1)[(21.6wzCmz2λmaxˉ1.6wzCmz2λmaxˉ2)npleLT]Wel,zWpl,zC_{zz} = 1+\Big(w_z-1\Big)\Big[\Big(2-\frac{1.6}{w_z}C_{mz}^2\bar{\lambda_{max}}-\frac{1.6}{w_z}C_{mz}^2\bar{\lambda_{max}}^2\Big)n_{pl}-e_{LT}\Big]\geq \frac{W_{el,z}}{W_{pl,z}}
Czz=1+(wz1)[(21.6wzCmz2λmaxˉ1.6wzCmz2λmaxˉ2)npleLT]=0.9120.651\begin{align*} C_{zz} &= 1+\Big(w_z-1\Big)\Big[\Big(2-\frac{1.6}{w_z}C_{mz}^2\bar{\lambda_{max}}-\frac{1.6}{w_z}C_{mz}^2\bar{\lambda_{max}}^2\Big)n_{pl}-e_{LT}\Big] \\\\ &= 0.912 \geq 0.651 \end{align*}

where

npl=NEdNRk/γ1=1742083.85/1.0=0.083n_{pl} = \frac{N_{Ed}}{N_{Rk}/\gamma_1} = \frac{174}{2083.85/1.0} = 0.083
wy=Wpl,yWel,y=720646=1.111.5w_y = \frac{W_{pl,y}}{W_{el,y}} = \frac{720}{646} = 1.11 \leq 1.5
wz=Wpl,zWel,z=166108=1.541.5wz=1.5w_z = \frac{W_{pl,z}}{W_{el,z}} = \frac{166}{108} = 1.54 \geq 1.5 \therefore w_z = 1.5
λˉmax=max(λyˉ,λzˉ)=1.51\bar{\lambda}_{max} = \textrm{max}\Big(\bar{\lambda_y},\bar{\lambda_z}) = 1.51
Step 6.6 Calculate the interaction factors Kyy,Kyz,Kzy,KzzK_{yy}, K_{yz}, K_{zy}, K_{zz}

Noting that our section is class 1:

kyy=CmyCmLTμy1NEdNcr,y1Cyy=1.0014×1.23×0.9811742533×10.97=1.336k_{yy} = C_{my}C_{mLT}\frac{\mu_y}{1-\frac{N_{Ed}}{N_{cr,y}}}\frac{1}{C_{yy}} = 1.0014 \times 1.23 \times \frac{0.98}{1-\frac{174}{2533}}\times\frac{1}{0.97} = 1.336
kyz=Cmzμy1NEdNcr,z1Cyz0.6wzwy=0.86×0.981174917×10.79×0.61.51.11=0.918k_{yz} = C_{mz}\frac{\mu_y}{1-\frac{N_{Ed}}{N_{cr,z}}}\frac{1}{C_{yz}}0.6\sqrt{\frac{w_z}{w_y}} = 0.86\times\frac{0.98}{1-\frac{174}{917}}\times \frac{1}{0.79}\times0.6\sqrt{\frac{1.5}{1.11}} = 0.918

where

μy=1NEdNcr,y1χyNEdNcr,y=1174253310.72×1742533=0.98\mu_y = \frac{1-\frac{N_{Ed}}{N_{cr,y}}}{1-\chi_y\frac{N_{Ed}}{N_{cr,y}}} = \frac{1-\frac{174}{2533}}{1-0.72\times\frac{174}{2533}} = 0.98
CmLT=Cmy2aLT(1NEdNcr,z)(1NEdNcr,T)=1.00142×0.997(1174917)(1174917)=1.23C_{mLT} = C_{my}^2\frac{a_{LT}}{\sqrt{\Big(1-\frac{N_{Ed}}{N_{cr,z}}\Big)\Big(1-\frac{N_{Ed}}{N_{cr,T}}\Big)}} = 1.0014^2\times \frac{0.997}{\sqrt{\Big(1-\frac{174}{917}\Big)\Big(1-\frac{174}{917}\Big)}} = 1.23

and

kzy=CmyCmLTμz1NEdNcr,y1Czy0.6wywz=1.0014×1.23×0.8711742533×10.845×0.61.111.5=0.703k_{zy} = C_{my}C_{mLT}\frac{\mu_z}{1-\frac{N_{Ed}}{N_{cr,y}}}\frac{1}{C_{zy}}0.6\sqrt{\frac{w_y}{w_z}} = 1.0014\times 1.23 \times \frac{0.87}{1-\frac{174}{2533}}\times \frac{1}{0.845}\times 0.6\sqrt{\frac{1.11}{1.5}} = 0.703
kzz=Cmzμz1NEdNcr,z1Czz=0.86×0.871174917×10.912=1.013k_{zz} = C_{mz}\frac{\mu_z}{1-\frac{N_{Ed}}{N_{cr,z}}}\frac{1}{C_{zz}} = 0.86\times \frac{0.87}{1-\frac{174}{917}}\times\frac{1}{0.912} = 1.013

where

μz=1NEdNcr,z1χzNEdNcr,z=117491710.34×174917=0.87\mu_z = \frac{1-\frac{N_{Ed}}{N_{cr,z}}}{1-\chi_z\frac{N_{Ed}}{N_{cr,z}}} =\frac{1-\frac{174}{917}}{1-0.34\times\frac{174}{917}} = 0.87
Step 6.7 Evaluate the combined buckling resistance inequalities

All that remains is to plug in all of the calculated values which we have determined in the preceding steps. If we can show that the value of both inequalities is less than unity then our member has sufficient buckling resistance against the combined bending moments and axial compression.

NEdχyNRkγM1+kyyMy,EdχLTMy,RkγM1+kyzMz,EdMz,RkγM1\frac{N_{Ed}}{\chi_y\frac{N_{Rk}}{\gamma_{M1}}} + k_{yy}\frac{M_{y,Ed}}{\chi_{LT}\frac{M_{y,Rk}}{\gamma_{M1}}} + k_{yz}\frac{M_{z,Ed}}{\frac{M_{z,Rk}}{\gamma_{M1}}}
(1740.72×2083.851.0)+(1.336×42.030.66×255.61.0)+(0.918×7.8758.931)=0.571.0\Big(\frac{174}{0.72\times\frac{2083.85}{1.0}}\Big) + \Big(1.336\times\frac{42.03}{0.66\times\frac{255.6}{1.0}}\Big) + \Big(0.918\times\frac{7.87}{\frac{58.93}{1}}\Big) = 0.57 \leq 1.0
NEdχzNRkγM1+kzyMy,EdχLTMy,RkγM1+kzzMz,EdMz,RkγM1\frac{N_{Ed}}{\chi_z\frac{N_{Rk}}{\gamma_{M1}}} + k_{zy}\frac{M_{y,Ed}}{\chi_{LT}\frac{M_{y,Rk}}{\gamma_{M1}}} + k_{zz}\frac{M_{z,Ed}}{\frac{M_{z,Rk}}{\gamma_{M1}}}
(1740.34×2083.851)+(0.703×42.030.66×255.61)+(1.013×7.8758.931)=0.561.0\Big(\frac{174}{0.34\times\frac{2083.85}{1}}\Big) + \Big(0.703\times\frac{42.03}{0.66\times\frac{255.6}{1}}\Big) + \Big(1.013\times\frac{7.87}{\frac{58.93}{1}}\Big) = 0.56 \leq 1.0

Since both inequalities have been satisfied, we have shown that the two governing criteria are met and therefore the proposed section has adequate buckling resistance against the combined action of axial compression and bi-axial bending. In fact, the utilisation is around 60% and therefore if we wanted to, the section size could be further optimised.

4.0 Initial sizing of compression members

The previous design example demonstrates the potential complexity involved in designing compression members to Eurocode 3. The chosen design example aimed to be as broad as possible to allow a full breadth of design checks to be conducted.

In reality, most columns are simple, as explained at the start of this tutorial. As a result, we can greatly simplify the design checks involved.

4.1 Simplified approach for simple columns

We can simplify the governing buckling criteria for simple columns down to the following expression

NEdχminAfy+My,EdχLTWyfy+1.5Mz,EdWzfy1.0\frac{N_{Ed}}{\chi_{min}Af_y} + \frac{M_{y,Ed}}{\chi_{LT}W_yf_y} + 1.5\frac{M_{z,Ed}}{W_zf_y} \leq 1.0

Which I am sure you will agree is much simpler than the two buckling criteria we used in our worked example.

A further simplification can be afforded by taking λˉLT=0.9λˉmin\bar{\lambda}_{LT} = 0.9\bar{\lambda}_{min}. So, in reality, a simple column can be designed in two steps:

  1. Check the combined resistance using the following:
NEdNRd+My,EdMy,Rd+Mz,EdMz,Rd1.0\frac{N_{Ed}}{N_{Rd}}+\frac{M_{y,Ed}}{M_{y,Rd}}+\frac{M_{z,Ed}}{M_{z,Rd}} \leq 1.0
  1. Check the combined buckling resistance using the following:
NEdχminAfy+My,EdχLTWyfy+1.5Mz,EdWzfy1.0\frac{N_{Ed}}{\chi_{min}Af_y} + \frac{M_{y,Ed}}{\chi_{LT}W_yf_y} + 1.5\frac{M_{z,Ed}}{W_zf_y} \leq 1.0

If we step back and look at these expressions, the only term which we need to perform multiple calculation steps on is χmin\chi_{min} - and we know from the worked example that this is one of the simpler Eurocode 3 calculations.

To make your life even easier, you can make an initial assumption that χmin\chi_{min} = 0.7, which allows a section to be chosen on the basis of area alone, using the following expression.

Amin=NEdγM1χminfyA_{min} = \frac{N_{Ed}\gamma_{M1}}{\chi_{min}f_y}

Once you've picked your section, then the above two criterion can be checked reasonably quickly.

5.0 Wrapping up and Coming up next

I hope you have found this tutorial on the design of steel compression members to Eurocode 3 helpful. This tutorial is intended to give you the skills and confidence to use Eurocode 3 Part 1-1 to design steel compression members, including those which must resist bi-axial bending in addition to axial compression. You should now have a good understanding of the different design checks which we conduct and the reason why we conduct them.

In our next tutorial in this series, we focus our attention on the design of steel trusses. As we know, Eurocode 3 does not have specific sections for this type of application and instead categorises members based on the internal forces which they must resist...can you think which internal forces we will be considering for steel members within truss structures?

We will answer all of these questions, and quite a few more, in the next instalment.

getting-started
Callum Wilson
BSc (Hons), BEng (Hons), MSc
Hi, I'm Callum, a regular contributor on EngineeringSkills.com. I hope you have found this content interesting and useful. I am excited to be producing content for EngineeringSkills.com as I enjoy the challenge of communicating technical content in a simple and accessible format. Feel free to get in touch or follow EngineeringSkills on any of the social accounts.

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.