(Eng) Advanced Concept Training Steel 2011.0 - Nemetschek Scia

Download as pdf or txt
Download as pdf or txt
You are on page 1of 123

Advanced Concept Training

Steel Code Check


Steel

All information in this document is subject to modification without prior notice. No part or this manual
may be reproduced, stored in a database or retrieval system or published, in any form or in any way,
electronically, mechanically, by print, photo print, microfilm or any other means without prior written
permission from the publisher. Scia is not responsible for any direct or indirect damage because of
imperfections in the documentation and/or the software.
© Copyright 2011 Nemetschek Scia nv. All rights reserved.

2
Table of contents

Contents
Introduction................................................................................................................................... 1
Cross-sections.............................................................................................................................. 2
Materials ........................................................................................................................................ 6
Section classification ................................................................................................................... 9
ULS Section Check..................................................................................................................... 17
Partial Safety factors....................................................................................................................... 17
Tension ............................................................................................................................................. 17
Compression ................................................................................................................................... 18
Bending moment ............................................................................................................................. 19
Shear................................................................................................................................................. 20
Torsion ............................................................................................................................................. 21
Combined check: Bending, shear and axial force ....................................................................... 26
ULS Stability Check.................................................................................................................... 29
Classification ................................................................................................................................... 29
Flexural buckling check.................................................................................................................. 30
Buckling factors ............................................................................................................................ 32
Buckling length ............................................................................................................................. 41
Buckling factors/lengths: manual input ........................................................................................ 44
Flexural buckling check in Scia Engineer .................................................................................... 45
Torsional Buckling .......................................................................................................................... 46
Lateral Torsional Buckling ............................................................................................................. 48
General ......................................................................................................................................... 48
Calculation of Mcr.......................................................................................................................... 53
LTB Restraints.............................................................................................................................. 56
Diaphragms .................................................................................................................................. 58
Lateral Torsional Buckling using LTBII – Not in Concept Edition ................................................ 62
Compression and bending check .................................................................................................. 64
Shear buckling check – EN 1993-1-5 ............................................................................................. 66
General ......................................................................................................................................... 66
Stiffners ........................................................................................................................................ 68
ULS Check for Battened compression members .................................................................... 70
ULS Check for Built-in beams ................................................................................................... 75
Optimisation................................................................................................................................ 78
Cross section optimisation ............................................................................................................ 78
Overall optimisation ........................................................................................................................ 82
nd
2 order calculation and imperfections .................................................................................. 85
Overview .......................................................................................................................................... 85
Alpha critical – Not in concept edition .......................................................................................... 87
Global frame imperfection ϕ........................................................................................................... 87
Bow imperfection ............................................................................................................................ 90
NEd > 25% Ncr ............................................................................................................................... 90
Bow imperfection e0 ..................................................................................................................... 91
Buckling shape as imperfection - ηcr – Not in concept edition .................................................. 93
SLS Check ................................................................................................................................... 98
Nodal displacement ........................................................................................................................ 98
Relative deformation ....................................................................................................................... 98
Fire resistance Check - Not in concept edition .................................................................... 102
General ........................................................................................................................................... 102
Temperature time curves ........................................................................................................... 102
Steel temperature ....................................................................................................................... 103
Steel properties .......................................................................................................................... 104

3
Steel

Fire resistance properties in Scia Engineer ............................................................................... 106


Resistance domain........................................................................................................................ 107
Principle...................................................................................................................................... 108
Example in Scia Engineer .......................................................................................................... 108
Temperature/time domain ............................................................................................................ 112
Principle...................................................................................................................................... 112
Example in Scia Engineer .......................................................................................................... 113
References and literature ........................................................................................................ 117

4
Advanced Concept Training

Introduction
This course will explain the calculation of steel following the EN 1993-1-1: General rules and rules for
buildings and EN 1993-1-2: Design of Steel Structures: Structural fire design.

Most of the options in the course can be calculated/checked in Scia Engineer with the Concept
edition.

For some supplementary checks an extra module (or edition) is required, but this will always be
indicated in those paragraphs.

1
Steel

Cross-sections
Scia Engineer will use the axes y-y and z-z respectively for the major and minor principal axes of the
cross section.

If the principal axes not coincide with the y-y and z-z axes following the EN 1993-1-1, also those axes
are indicated:

EN 1993-1-1 Scia Engineer EN 1993-1-1 Scia Engineer

In the steel code check, the strong axis will always be considered as the y-y axis. So when inputting a
profile in which the local y-axis does not corresponds with the strong axis, the axes are switched in the
steel code check and Scia Engineer will give a message about the axes.

Example: Local Axes.esa

The strong axis of this cross section in the local z-axis:

2
Advanced Concept Training

This beam is only loaded by a line load in the z-direction, resulting in a moment My and a shear force
Vz on the beam.

When looking at the steel code check, the axes are switched and the strong axis is taken as local y-
axis:

This is indicated in steel code check and also the internal forces are switched following the new local
axes.

The following sections are classified during the steel code checks

So as I profile only the indicated profiles are recognized. The other I-forms are not recognized as an I-
profile, and e.g. a classification cannot be executed for those forms.

I Symmetric I shape
from library

thin walled geometric

sheet welded Iw

RHS Rectangular Hollow Section


from library

thin walled geometric

CHS Circular Hollow Section


from library

thin walled geometric

L Angle Section
from library

thin walled geometric

U Channel Section

3
Steel

from library

thin walled geometric

T T Section
from library

thin walled geometric

sheet welded Tw

PPL Asymmetric I shape


from library

thin walled geometric

sheet welded Iwn

Z Z Section
from library

thin walled geometric

RS Rectangular Sections
from library

thin walled geometric

O Solid Tube
from library

thin walled geometric

Σ Cold formed sections

from library

pairs 2CFUo

pairs 2CFUc

4
Advanced Concept Training

pairs 2CFCo

pairs 2CFCc

pairs 2CFLT

General Cross Section: thin walled

NUM Numerical section


Numerical section

COM all others

The standard profile sections (from library) have fixed sections properties and dimensions, which have
to be present in the profile library. The section properties and the required dimension properties are
described in the Theoretical Background (Ref.[32]).

Example: Cross Section.esa


- 3 cross sections:
o HEA300 from library
o HEA300 inputted as a general cross section (Imported from dwg)
o HEA 300 inputted as a general shape
- Only the first profile is recognized as a symmetric I-shape
- The second and third profile are not recognized as a symmetric I-shape
- The first profile will be classified as an I-profile and a plastic check will be performed.
- The second and third profile cannot be classified, so an elastic check will be performed. This will
results in different checks:

Remark:
The general cross section cannot be imported in Scia Engineer with the concept edition. For
this option the Professional or Expert edition is necessary or module esa.07.

5
Steel

Materials

For standard steel grades, the yield strength fy and tensile strength fu are defined according to the
thickness of the element.

The standard steel grades as defined in Table 3.1 of the EN 1993-1-1 are:

6
Advanced Concept Training

Those materials are included in Scia Engineer:

The Histar and Fritenar steel grades have been implemented according to Arcelor.

With the option “Thickness range” the influence of the thickness on the yield strength fy and tensile
strength fu are defined. When inputting a new user defined material in Scia Engineer, also this
thickness range for this material can be inputted:

Example: NA_Material_Strength_Application.esa

In this example a material of the National Annex of Belgium has been inputted manually: S 275 J2.
This material has a lot of thickness range strength reductions:

7
Steel

And for each beam the correct yield strength has been taken into account according to the inputted
table above:

CS1 CS2 CS3 CS4

8
Advanced Concept Training

Section classification
The classification is executed following EN 1993-1-1, art. 5.5.

Four classes of cross sections are defined:


• Class 1 (EC3, NEN) or PL-PL (DIN) section
Cross sections which can form a plastic hinge with the rotation capacity required for plastic
analysis

• Class 2 (EC3, NEN) or EL-PL (DIN) section


Cross sections which can develop their plastic moment resistance, but have limited rotation
capacity.

• Class 3 (EC3, NEN) or EL-EL (DIN) section


Cross sections in which the calculated stress in the extreme compression fibre of the steel
member can reach its yield strength, but local buckling is liable to prevent development of the
plastic moment resistance.

• Class 4 (EC3, NEN) or Slender section


Cross sections in which it is necessary to make explicit allowance for the effects of local buckling
when determining their moment resistance or compression resistance.

This classification depends on the proportions of each of its compression elements.


For each intermediary section, the classification is determined and the related section check is
performed. The classification can change for each intermediary point.

For each load case/combination, the critical section classification over the member is used to perform
the stability check. So, the stability section classification can change for each load case/combination.
However, for non-prismatic sections, the stability section classification is determined for each
intermediary section.

9
Steel

The classification check in Scia Engineer will be executed following tables 5.2 of the EN 1993-1-1:

10
Advanced Concept Training

11
Steel

If this classification results in a class 4 profile, the effective cross section will be calculated according to
EN 1993-1-5.

For each load case and combination, the most critical effective area properties are saved:
• Aeff is the effective area of the cross section when subject to uniform compression.
• W eff is the effective section modulus of the cross-section when subject only to moment about the
relevant axis.
• eN is the shift of the relevant centroidal axis when the cross section is subject to uniform
compression.

With these critical properties, the steel code check is performed.

Example: Industrial Hall.esa


In this example the classification will be done for an IPE750x137 profile, resulting in a class 4 cross
section. And afterwards the calculation of the effective shape following EN 1993-1-5 will be given.
Consider column B28

The classification has been executed at the bottom of the column (position = 0.00m)

On this position a hinged support has been inputted, so on this position the column is not subjected to
bending.

12
Advanced Concept Training

Table 5.2 (sheet 1): Table 5.2 (sheet 2):

c = H – 2 * thickness flange – 2* radius c = B/2 – thickness web/2 – radius


c = 753mm – 2x17mm – 2x17mm = 685mm c = 263mm/2 – 11,5mm/2 – 17 = 109
t = 11,5 mm t = 17 mm
c/t = 59,57 c/t = 6,40

With  = 235/ = 1 With  = 235/ = 1


Maximum ratio class 1: 33 Maximum ratio class 1: 9
Maximum ratio class 2: 38 Maximum ratio class 2: 10
Maximum ratio class 3: 42 Maximum ratio class 3: 14

c/t = 59.57 > 42 => Class 4 c/t = 6.40 < 9 => Class 1

In Scia Engineer:

13
Steel

This cross section has a classification: class 4, so effective properties has to be calculated.
Those properties are also given in Scia Engineer in the preview of the steel code check, just below the
title “STABILITY CHECK”.

The calculation of the sectional area Aeff is given in below.

In this section there is a uniform compression force over the web:

The calculation of the effective section will be performed following EN 1993-1-5:2006, Tabel 4.1 and
EN 1993-1-5:2006/AC:2009 article 9). In Table 4.1 the uniform compression situation will be used in
this example:

14
Advanced Concept Training


≤ 0,5 + 0.085 − 0.055 = 0,673 =1

> 0,5 + 0.085 − 0.055 = 0,673 − 0.22)/
 = ( 

$ /%

= [ / ]",# =
28.4'(

ε = 1,00 (S235)

kσ = 4,0 (Tabel 4.1 of EN 1993-1-5:2006)

$ /% 59.57
=
= = 1.05
28.4'( 28.4 ∙ 1.00 ∙ √4.0

1.05 − 0.22
= = 0,75
(1.05)

b./ = b. = 0.5 ∙ b.00 = 0,5 ∙ ρ ∙ b = 0,5 ∙ 0,75 ∙ 685 mm = 256,88 mm

4566 = [263 ∙ 17 + 256,88 ∙ 11,5 + 17 ∙ 11,5] ∙ 2 = 15152 mm²

In this calculation the roundings at the corners between the flanges and the web are not taken into
account. Therefore the result in Scia Engineer will be a bit higher:

Also the other properties of this effective cross section can be calculated.

15
Steel

The calculated classification in Scia Engineer can be overruled by 2 settings in the Steel member
data:

Section classification

The user can choose between a classification calculation by program. Or can overwrite this and
choose for class 1, 2 or 3. Since classification 4 is not described for all cross section in the
Eurocode, this option can not be chosen.

Elastic Check Only

The user can choose to perform an elastic check. This corresponds with a class 3 check.

Classification in Scia Engineer

Following classes are possible per profile type in Scia Engineer:

Section Classification Class 1 Class2 Class3 Class 4 Warping Warping


check Class
1&2 check Class 3

I x x x x x x x
RHS x x x x x
CHS x x x x
L x x
U x x x x
T x x
PPL x x x x
RS x x
Z =3 x
O =3 x

Σ x x x x

NUM =3 x

COM =3 x

And in the chapter “Cross sections” an explanation is given which sections are recognized under which
form type (I, RHS, L, CHS, ...).

16
Advanced Concept Training

ULS Section Check

In this chapter, first the Partial safety factors are explained and afterwards a short explanation of all the
section checks is given.

The section check can be found in Scia Engineer under “Steel -> Beams -> ULS check -> Check”. In
this menu the user can choose to look at the “Brief” output and or the “Detailed output”.

With the brief output the results are shown in one line.
With the detailed output, the results of all the unity checks are shown including a reference to the used
formula in the EN 1993-1-1 for each check.
In below the Section Check will be explained and the detailed results are shown.

Partial Safety factors


The partial safety factors are taken from EN 1993-1-1 art. 6.1.

And the following safety factors are taken into account:

γM0 = 1,00 Resistance of cross-section


γM1 = 1,00 Resistance of members to instability accessed by member checks
γM0 = 1,25 Resistance of cross-section in tension to fracture

Those factors can also be found in the National Annex of EN 1993-1-1 in Scia Engineer:

Tension
The tension check will be executed following EN 1993-1-1 art. 6.2.3.

89:
≤1
8;,<:
=∙6>
8;,<: = ? : the design plastic resistance of the cross-section
@A

Example: Industrial Hall.esa

Consider bracing B234 (for load case 3DWind1). In this load case the bracing in under tension:

17
Steel

1,7663 ∙ 10 BB ∙ 235CDE


8;,<: = = 41,51'8
1.0
32,32'8
FGH% IℎKL' = = 0,78
41,51'8

Compression
The compression check will be executed following EN 1993-1-1 art. 6.2.4.

89:
≤1
8,<:

Where
=∙6>
- 8,<: = ? For class 1, 2 or 3 cross-sections
@A
=MNN ∙6>
8,<: =
?@A
- For class 4 cross-sections

Example: Industrial Hall.esa

Consider column B28 (for combination CO1-ULS).

The classification on this position is Class 1 (column under compression and bending).

1,7500 ∙ 10O BB ∙ 235CDE


8,<: = = 4112,5 '8
1.0
161.00'8
FGH% IℎKL' = = 0,04
4112.5'8

18
Advanced Concept Training

Bending moment
The bending moment check for My and Mz will be executed following EN 1993-1-1 art. 6.2.5.

C9:
≤1
C,<:

Where

QRS ∙6>
C,<: = C P,<: =
?@A
- For class 1 or 2 cross-sections

QMS,TUV ∙6>
- C,<: = C5P,<: = ?@A
For class 3 cross-sections

QMNN ∙6>
C,<: =
?@A
- For class 4 cross-sections

Example: Industrial Hall.esa

Consider column B28 (for combination CO1-ULS).

The classification on this position is Class 1 (column under compression and bending).

19
Steel

W P, ∙  4.865 ∙ 10Z BB³ ∙ 235CDE


C,,<: = = = 1,14328 ∙ 10\ 8BB = 1143,28 '8B
XY" 1.0

W P,] ∙  6.141 ∙ 10# BB³ ∙ 235CDE


C],,<: = = = 1,4431 ∙ 10^ 8BB = 144,31 '8B
XY" 1.0

300.30 '8B
FGH% IℎKL' C = = 0,26
1143,28 '8B

0.13 '8B
FGH% IℎKL' C_ = = 0,00
144,31 '8B

Shear
The shear check for Vy and Vz will be executed following EN 1993-1-1 art. 6.2.6.

`9:
≤1
`,<:

20
Advanced Concept Training

For plastic design Vc,Rd the absence of torsion, is the design plastic shear resistance Vpl,Rd:
=a ∙b6> /√cd
` P,<: =
?@A

With:
Av: the shear area.
The formula for Av depends on the cross-section (see EN 1993-1-1 article 6.2.6(3)).

For elastic design Vc,Rd is the design elastic shear resistance. The following criterion for a critical point
of the cross-section may be used unless the buckling verification in section 5 of EN 1993-1-5 applies:
e9:
≤1
 ∙ f3 ∙ XC0 g

e9: =
hij ∙k
l ;
Where

Where the shear force is combined with a torsional moment, the plastic resistance Vpl,Rd should be
reduced as specified in the next paragraph.

Example: Industrial Hall.esa

Consider column B28 (for combination CO1-ULS).

The classification on this position is Class 1 (column under compression and bending).

Torsion
EN 1993-1-1 article 6.2.7.

m9:
≤1
m<:

Where TRd is the design torsional resistance of the cross-section.

21
Steel

m9: = m;,9: + mn,9:


TEd should be considered as the sum of two internal effects:

Where Tt,Ed is the internal St. Venant torsion


Tw,Ed is the internal warping torsion

Scia Engineer will take into account the St. Venant torsion automatically. If the user would like to
calculate also with warping torsion, this option should be activated in the “Buckling and relative
lengths” of this beam:

The check box Warping check should be activated and for the beginning and end of the beam, the
user should indicate if this end is Fixed or Free for warping.

Example: Warping.esa

As quoted in the EN 1993-1-1, art. 6.2.7(7), for closed hollow sections the effects of torsional warping
can be neglected and in the case of a member with an open cross section, such as I of H, it may be
assumed that the effects of St.Venant torsion can be neglected.

In this example 4 beams are inputted:


- IPE 180 – No Warping has been activated
- SHS 180/180/10.0 – No Warping has been activated
- IPE 180 – Warping has been activated in the buckling data
- SHS 180/180/10.0 –Warping has been activated in the buckling data

On the IPE 180 beams a torsional moment of 1 kNm has been inputted.
On the SHS 180/180/10.0 a torsional moment of 30 kNm has been inputted.

22
Advanced Concept Training

For the SHS profiles, there is no difference between the option “Warping” activated or not. This was
expected as warping can be neglected for closed hollow section. Both beams give the following
torsional check:

For IPE profiles, the warping torsion cannot be neglected. So when this is activated, it has a big
influence on the resistance value for torsion:

Torsion check IPE 180 without warping activated:

Torsion check IPE 180 with warping activated:

23
Steel

In this case the difference between the option “warping” activated or not, give a big difference in the
unity check.

The torsional properties of the cross sections are not calculated correctly in a lot of cases.
The next example give an explanation about this calculation.

Example: Torsion.esa

In this example two beams are inputted, with the same properties, only the torsional properties of the
cross sections are different:

24
Advanced Concept Training

CS1 CS2
The option “FEM analysis” is not activated: The option “FEM analysis” is activated:

So the value for It and Iw will be calculated with a


FEM analysis:

So the value for It will be calculated with a


general formula, which does only give an
estimation of the It, but not a correct value. The
value for Iw will not be calculated and Iw will
6
taken as 0.00 mm .

On the two beams a torsional moment (8 kNm) has been inputted and a small line load (-0.1 kN/m).

This result in the following Section Check:

25
Steel

The difference in the results is in the torsion check. For the first beam no torsion check can be
executed, since no results for torsional stiffnesses can be found.

For beam B2 (CS2) the torsion check can be executed:

Combined check: Bending, shear and axial force


This check will be executed following EN 1993-1-1 articles 6.2.8 – 6.2.9 -6.2.10.

For bi-axial bending the following criterion may be used:

C,9: C],9:
r s
o q +o q ≤1
Cp,,<: Cp,],<:

in which α and β are defined as follows:

- I and H sections:
α = 2; β = 5n but β > 1
- Circular hollow sections:
α = 2; β = 2
- Rectangular hollow sections:

t = u =
/,ZZ
but α = β < 6
/v/,/cw²

G=p
pij
Where
RS,xj

26
Advanced Concept Training

The values MN,y,Rd and MN,z,Rd depends on the moment resistance, reduced with a factor depending on
“n”, the check of the normal force:

Cp,,<: = C P,,<: (1 − G)/(1 − 0,5En )

Cp,],<: = C P,],<: (1 − G)/b1 − 0,5E6 d

And Mpl,Rd depends on the yield strength fy.

If yz{ ≥ }, ~ y€,,‚{

The yield strength will be reduced with a factor ρ.

(1 − )

2 `9:
Where

=ƒ − 1„
` P,<:

And
=a ∙b6> /√cd
` P,<: = ?@A

When torsion is present ρ should be obtained from

2 `9:

=ƒ − 1„
` P,…,<:

Class 3 cross-sections

EN 1993-1-1 article 6.2.9.2.

In absence of shear force, for Class 4 cross-sections the maximum longitudinal stress should satisfy
the criterion:

6>
†,9: ≤ ?
@A

Class 4 cross-sections

EN 1993-1-1 article 6.2.9.3.

In absence of shear force, for Class 4 cross-sections the maximum longitudinal stress should satisfy
the criterion:

6>
†,9: ≤ ?
@A

The following criterion should be met:

89: C,9: + 89: ∙ Kp C],9: + 89: ∙ Kp]


+ + ≤1
  
4566 ∙ W566,,‡ˆw ∙ W566,],‡ˆw ∙
XY" XY" XY"

Where Aeff is the effective area of the cross-section when subjected to uniform compression
W eff, min is the effective section modulus of the cross-section when subjected only to moment
about the relevant axis

27
Steel

eN is the shift of the relevant centroidal axis when the cross-section is subjected to
compression only.

Example: Industrial hall.esa

Consider column B28 (for combination CO1-ULS).

 n = 0,04

 Unity check for shear force is smaller than 0,5, thus no reduction of the yield strength for the
combined check.

28
Advanced Concept Training

ULS Stability Check

The stability check can be found in Scia Engineer under “Steel -> Beams -> ULS check -> Check”. In
this menu the user can choose to look at the “Brief” output and or the “Detailed output”.

With the brief output the results are shown in one line.
With the detailed output, the results of all the unity checks are shown including a reference to the used
formula in the EN 1993-1-1 for each check.
In below the Section Check will be explained and the detailed results are shown.

Not all stability checks are inputted in Scia Engineer for all cross sections. In the table below an
overview of the possible checks for each section is given:

Section Buckling LTB Shear buckling


I x x x
RHS x x
CHS x x
L x (*2)
U x (*2) x
T x (*2)
PPL x x x
RS x (*2)
Z x (*2)
O x (*2)
Σ (*1) (*2) x

NUM (*1) (*2)


COM (*1) (*2)

(*1) Buckling curve are introduced manually by user


(*2) General formula for Mcr – see chapter ‘Lateral Torsional Buckling’

And in the chapter “Cross sections” an explanation is given which sections are recognized under which
form type (I, RHS, L, CHS, ...).

Classification
There is a difference on the section of classification in the Section check and the Stability check:

Classification in the section check:


Here the classification is done for each section on the member and afterwards the Section check will
be executed with the classification and the internal forces on this section.

Classification in the Stability Check


In the stability check the highest classification along the member is used for the stability check on all
the sections.

So when having a beam with the highest check on 2m from the beginning of the beam:
- The section check will take into account the internal forces on 2m from the beginning of the beam
and perform a classification with those internal forces.

29
Steel

- The stability check will perform a stability check with the internal forces on 2m of the beginning of
the beam, but not with the classification on this place. Perhaps this section has the highest
moment in comparison with a compression force. The classification can result here in Class one.
But in the beginning of the beam, the bending moment is equal to zero, so perhaps the
classification here will be Class 4. The stability check will thus take the internal forces of section
2m on the member, but take into account the classification Class 4 (from 0m on the member).

Flexural buckling check


The flexural buckling check will be executed following EN 1993-1-1 art. 6.3.1.

89:
≤1
8‰,<:

Where
Š∙=∙6>
8‰,<: =
?@‹
- For class 1, 2 or 3 cross sections

Š∙=MNN ∙6>
8‰,<: =
?@‹
- For class 4 cross sections

The reduction factor χ will be calculated as follows:

Œ=
/
‘
Ž²v
but χ < 1,0

Φ 0,5 “1 + tb ̅ − 0,2d + ̅ •
with

̅
=∙6>
–p For class 1, 2 or 3 cross-sections
—˜
=MNN ∙6>
–
p—˜
For class 4 cross-sections

8 =
™²9l
Ncr Critical normal force (Euler force)

š²›²

α Imperfection depending on the buckling curve:

The selection of the buckling curve for a cross section is done with EN 1993-1-1, Table 6.2

30
Advanced Concept Training

In Scia Engineer the calculation of the buckling coefficient ky or kz can be done automatically or can be
inputted by the user.

In the next paragraph the calculation of the Buckling factors is explained.

Example: Buckling Curves.esa

In Table 6.2 general buckling curves are given for the most common profiles. It seems that those
curves can be changed in Scia Engineer, but the program will always take the buckling curves as
defined in the EN 1993-1-1 into account.

But for general sections, or sections not available in this table 6.2 the user can choose a buckling curve
for this cross section.

31
Steel

In this example 3 beams are inputted with 3 different cross sections:

- B1: CS1 – IPE 180:

o buckling curve a for y-y according to code

o Stability check uses curve a, according to code

- B2: CS2 – IPE 180

o buckling curve a for y-y according to code

o Manual input by user: curve c

o Stability check uses curve a, according to code

- B3: CS3 – IPE180 + T50:

o non-standard section: no buckling curve according to code

o Manual input by user: curve c

o Stability check uses curve c, according to manual input

Buckling factors

General method

By default, for the calculation of the buckling ratios, two approximate formulas are used: one formula
for a non sway structure (resulting in a buckling factor smaller (or equal) than 1) and one formula for a
sway structure (resulting in a buckling factor higher (of equal) than 1.

So it is important for this method that the user chooses the correct option sway or non sway for the two
local directions:
yy: buckling around the local y-axis (so deformation in the direction of the local z-axis)
zz: buckling around the local z-axis (so deformation in the direction of the local y-axis)

32
Advanced Concept Training

The option “sway” or “non sway” can be chosen in the menu “Steel -> Beams -> Steel Setup” for the
whole structure:

It the option is checked (as here for buckling around the y-axis), it is indicated that the construction is
not braced enough and can be considered as “sway”.

This can also be changed for each beam separately, using the buckling data. This can be found under
“Steel -> Beams -> Member Check data -> Member buckling data”:

For Sway yy and Sway zz there are 3 options:


- According to Steel>Beams>Setup: The same option will be taken as in the steel setup as shown
above.
- Yes: this beam will be considered as “sway” in this direction.
- No: this beam will be considered as “non sway” in this direction.

33
Steel

The following formulas are used for the buckling ratios:

• for a non sway structure :

' = (œ ‹
(œ œ‘ Ž#œ‹ Ž#œ‘ ŽO)(œ‹ œ‘ ŽOœ‹ ŽOœ‘ Ž/)
‹ œ‘ Ž//œ‹ Ž#œ‘ ŽO)(œ‹ œ‘ Ž#œ‹ Ž//œ‘ ŽO)

• for a sway structure :

' = – +4
™²
œ‹ †

with k the buckling factor


L the system length
E the modulus of Young
I the moment of inertia
Ci the stiffness in node i
Mi the moment in node i
φi the rotation in node i

4/  + ž²/
=
ž  (/+  ) + 8/ 

Iˆ Ÿ
ˆ =
 ¡

Iˆ =
¢ˆ

The values for Mi and φi are approximately determined by the internal forces and the deformations,
calculated by load cases which generate deformation forms, having an affinity with the buckling form.
So when performing a linear calculation, in the background 2 additional load cases are calculated, just
to calculated the buckling factors for the elements.

This calculation is automatically done when calculating the construction linearly. So when calculating
non linear, the user should also perform a linear calculation otherwise no buckling factors are
calculated and no steel code check can be performed.

The following load cases are considered in the linear calculation for the calculation of the buckling
factors:

- load case 1:
o on the beams, the local distributed loads qy=1 N/m and qz=-100 N/m are used
o on the columns the global distributed loads Qx =10000 N/m and Qy =10000 N/m are
used.

- load case 2:
o on the beams, the local distributed loads qy=-1 N/m and qz=-100 N/m are used
o on the columns the global distributed loads Qx =-10000 N/m and Qy=-10000 N/m are
used.

The used approach gives good results for frame structures with perpendicular rigid or semi-rigid
beam connections. For other cases, the user has to evaluate the presented bucking ratios.

34
Advanced Concept Training

Example: Buckling Factor.esa

Consider Column B1:


- L = 4000mm
- Set as sway
- In node N1 : My = 0 kNm => C2 = ρ2 = 0.0
- This node N1 defines ρ2 because ρ2 is always the smallest of the two.
- In node N2 for Loadcase LC1:
o My1 = 79833 kNm
o φ1 = fiy = 1530.6 mrad
o C1 = My1/ φ1 = 79833 kNm / 1530.6 mrad = 52.158 kNm/mrad
10
= 5,2158 x 10 Nmm/rad
o E = 210 000 N/mm²
4
o Iy = 162700000 mm
‹A‹A ¤TT
#./#^∙ ∙O"""‡‡
= = 6.106
£U › ˜¥j
ρ1 = ¤
9l /"""" ∙/Z¦""""" ‡‡ §
o
TT‘

Oœ‹ œ‘ Ž™‘¨‹ O∙Z./"Z∙"."Ž™‘ ∙Z./"Z


= = = 1.0
™‘ (œ‹ Žœ‘ )Ž^œ‹ œ‘ ™‘ (Z./"ZŽ".")Ž^∙Z./"Z∙".""
o

' = – + 4 = 1.0– + 4 = 2.37


™‘ ™‘
œ‹ † Z.//∙/.""
o

™²∙/""""p/‡‡²∙/Z¦""""" ‡‡§
8 = = = 37522008 = 3752.2 '8
™²9l
š²›² (.c¦)‘ (O""")²
o

Those values can also be found in Scia Engineer:

Under “Steel -> Beams -> Steel slenderness” the buckling length ky will be found:

This value can also be found in the stability check through “Steel -> Beams -> ULS Check -> Check”
under the buckling parameters. And here also the critical normal force Ncr can be found:

35
Steel

Calculation of the buckling factors for crossing diagonals

As the previous default method is only valid for perpendicular connection, this can not be used for the
calculation of the buckling factors for diagonals.

In the DIN 1880 Teil 2, Table 15 a method is given for the calculation of the buckling factor for crossing
diagonals. In Scia Engineer this option is also implemented. In this method the buckling length sK is
calculated in function of the load distribution in the element and sK is not a purely geometrical data. The

This method is only applicable for 2 diagonals with a hinged or rigid connection in the middle. To use
this functionality in Scia Engineer, the user has to connect the two diagonals with the option “Cross-
link”, which can be found in the menu Structure:

When connecting two beams with this option, in the buckling and relative lengths properties of the two
beams the option X-diagonals will be checked automatically:

And Scia Engineer will use the method from the DIN to calculate the buckling factors for the diagonals.

36
Advanced Concept Training

NOTE
When using cross-links, this option is automatically activated. The user must verify if this is wanted or
not. For example, when modelling purlins and rafters using cross-links, the option crossing diagonals
may not be activated.

With:
sK buckling length
L member length
l1 length of supporting diagonal
I moment of inertia (in the buckling plane) of the member

37
Steel

I1 moment of inertia (in the buckling plane) of the supporting diagonal


N compression force in member
N1 compression force in supporting diagonal
Z tension force in supporting diagonal
E elastic modulus

Calculation of the buckling length for a VARH-element

For a VARH element, Scia Engineer will use another calculation for the buckling length.

A VARH element is defined as follows:


The member has the properties of a symmetric I section where only the height is linear variable along
the member. The system length for buckling around the local yy axis (strong axis), is equal to the
member length.
For a VARH element (form node i to node j), we can define

L beam length
Ii, Ij moment of inertia at end i and j
Ai, sectional area at end i and j
Aj
E modulus of Young
Ncr critical Euler force
Ri , beam stiffness at end i and j
Rj

The stiffness R and R' is given by:


M
R=
Φ
L
R i′ = R i
EIi
L
R j′ = R j
EIi
Ij
ξ=
Ii

The critical Euler force is given by:

t² ¡ˆ
8 =
'²Ÿ²

To calculate α, the next steps are followed:

1. Calculate L, Ii, Ij, Ri, Rj, R'i, R'j, ξ

2. We suppose that
α 1
>
ξ -1 2

3. Calculate a, b, c and d as follows

38
Advanced Concept Training

α2 1
β= -
(ξ - 1 ) 2 4

1 1
a= [1 + (ξ - 1)( - β cotg(β ln ξ)]
α 2 2
1 β ξ -1
b = c = 2 [1 - ]
α sin (β ln ξ) ξ
1 (ξ - 1) 1
d= [1 + ( + β cotg(β ln ξ)]
α 2 ξ 2

4. For a beam in non-sway system, we solve


1 + a R i′ + d R j′ + (ad - bc) R i′ R j′ = 0

For a beam in sway system, we solve


R i′ (1 - a α2) + R j′ (1 - d α2) - α2 + R i′ R j′ (a - b - c + d - α2 (ad - bc)) = 0

5. When a solution is found, we check if


α 1
>
ξ -1 2

6. If not, then recalculate a,b,c en d as follows :


1 1
(ξ - 1)(( - β) ξβ - ( + β) ξ-β)
1 2 2
a = 2 [1 + ]
α ξβ - ξβ
1 2β(ξ - 1)
b = c = 2 [1 - ]
α ξ (ξβ - ξ-β)
1 1
(ξ - 1)(( - β) ξ-β - ( + β) ξβ)
1 2 2
d = 2 [1 + β -β
]
α ξ(ξ - ξ )

and resolve the proper equation of 4.

Example:

VARH.esa

350 mm 300 mm 750 mm 800 mm

Consider column B1
The VARH element going from 800mm to 300mm
The VARH is internally divided into a number of prismatic members. In this case, in 5 parts, but this
can be changed in the Mesh setup:

39
Steel

- in node N1:
o hi = 350 mm
-4 4
o Ii = 2.403 10 m
o M=0 → Ri, Ri’ = 0.0
- in node N1:
o fiy = 1512 mrad
o My = 79818 kNm
o E = 210000 N/mm²
o hj = 750 mm
9 4
o Ii = 1.363 x 10 mm
10
o Rj = 5.28 x 10
o R’j = 4.2
o ξ = 2.38
o α = 1.75
 8 = 9626.16 '8

40
Advanced Concept Training

Buckling length
In the previous paragraph the general calculation of the buckling factors has been explained for all type
of elements. With this buckling factor the buckling length of the beam will be calculated as follows:

© = ' ∙ Ÿ

With:
l the buckling length
k the buckling factor
L the system length

The system length of a beam is defined by the property “Buckling and relative lengths” of the beam.

Example: Buckling Length.esa

Consider column B4. This column has 3 nodes: N5, N10 and N6. The local X direction goes from the
bottom to the top of the column, so internal for Scia Engineer N5 is the first node of this column and N6
the last one.

In the properties window of this column the user can choose for “Buckling and relative lengths”. With
this option, the system length of the beam can be inputted. Default the following option will appear:

Explanation of the System Lengths:


- The first node (according to the local x-axis) is node N5, the last one is node N6.
- yy direction:
o This means around the local y-axis. So the column will deform in the z-direction.
o Around the y-axis, node N5 is supported. In node N10 no beam can be found in the local z-
direction in this point, thus column B3 is not supported around the y-axis in node N10. In
node N6 a horizontal beam in the local z-direction can be found and the column will be
supported around the local y-axis (yy) in node N6. This is indicated with the triangles in this
window:
 Supported in node N5
 Not supported in node N10
 Supported in node N6

41
Steel

- zz directions:
o This means around the local z-axis. So the column will deform in the y-direction.
o Around the z-axis, node N5 is supported. In node N10 a horizontal beam in the local y-
direction can be found and the column will be supported around the local z-axis (zz) in node
N10. Also in node N6 a horizontal beam in the local y-direction can be found and the column
will also be supported around the local z-axis (zz) in node N6. This is indicated with the
triangles in this window:
 Supported in node N5
 Supported in node N10
 Supported in node N6

- The system length will be taken as follows:


o Around the y-axis: the length between node N5 and N6: so 3m
o Around the z-axis: the length between node N5 and N10 for the first part of the beam (1,8m)
and the length between N10 and N6 for the second part of the beam: so 1.2m.
o This can also be found in the menu “Steel -> Beams -> Steel slenderness”:

o In this window the user can easily check the system length (Ly and Lz), the buckling factors
(ky and kz) and the buckling length (ly = ky x Ly and lz = kz x Lz).

- When clicking on Edit in the buckling data window, the triangles indicating the system lengths of
each part of a beam can be changed in the second tab “Buckling data”:

- Compare beam B3 and beams (B2+B30) with each other: they should have exactly the same system
lengths. The only difference between those columns is that beam B3 was inputted as a beam of 3m
and beams (B2+B30) are divided in two parts. Scia Engineer will consider those two beams also as
one buckling system:
o When the local axis are exactly in the same direction (so in this case the local x-axis is in the
same direction and the angle between the beams is exactly 180°).
o If no hinge has been inputted between the two beams.

42
Advanced Concept Training

- There is also a graphical representation of the system lengths:


o Select beam B2 (or B30)
o Go to the actions menu and click on Graphical input of system length:

o And now the user can choose between the yy and the zz axis:

- Now take a look at beam B13. The system lengths are the following (as expected, because there are
horizontal beams in the two directions on each node):

When looking at the rendered view, it will be clear that beam B16 is too weak to have an influence
on the system length of beam B13. In Scia Engineer there is possibility to exclude a beam from a
buckling system.
So select beam B16 and go in the Actions menu to Buckling data:

43
Steel

In this menu it is possible to indicate that beam B16 is a secondary beam and should not be taken
into account in the system lengths:

When looking now at member B13 and changing the Buckling and relative lengths back to “defaults”
the member B16 will not be included in the system lengths:

Buckling factors/lengths: manual input


The principles of the buckling factor and buckling length are explained in the previous paragraphs.
Those default settings can be changed in the menu “Buckling and relative lengths” or in the steel menu
with the “member buckling data”. This last one has a higher priority, so Scia Engineer will first look if
there are “member buckling data” on a beam. If yes, the program will take into account those
properties, if no, the program will look at the properties in “buckling and relative lengths” of the
properties window of the member.
In the menu “Steel -> Beams -> Member Check data -> Member buckling data” the user can change
the “Buckling ky, kz coefficients or buckling lengths” on “user input”.

44
Advanced Concept Training

And now the user can choose for ky and kz if those coefficients are calculated (default in Scia
Engineer) or input the buckling factor or even directly the buckling length:

And those properties can be placed on one or more members.

Flexural buckling check in Scia Engineer


Once all the buckling factors and system lengths have been inputted correctly, the Flexural buckling
check can be executed in Scia Engineer.

Example: Industrial hall.esa

Consider column B28:

The classification of beam B28 is class 4, so an effective cross section has been calculated:

Scia Engineer will first show the Buckling parameters of this beam:

45
Steel

In below the results of the yy direction are explained.

This direction has been set on “sway” in the steel Setup and with this option the Buckling factor k is
calculated.
The length of the column is 6.900m. So the buckling length = 2.60 x 6.900 m = 17.9 m

With this buckling length the Critical Euler load Ncr can be calculated. Afterwards the slenderness and
the Relative slenderness Lambda can be calculated with the critical Euler load.

An IPE750 profile has a buckling curve a, resulting in an Imperfection factor alpha = 0.21.

With those properties the reduction factor χ will be calculated, which will filled in the following formula:

Š∙=MNN ∙6>
8‰,<: = ?@‹

This result in a buckling resistance Nb,Rd = 3207.07 kN for flexural buckling around the local y-axis.

Exactly the same principle can be repeated for flexural buckling around the local z-axis.
This will result in a lower buckling resistance: Nb,Rd = 1754,66 kN.

The lowest buckling resistance will be used in the flexural buckling check:

Torsional Buckling
The check on torsional buckling can be important for profiles in which the position of the shear centre is
not the same as the centre of gravity of this section.

The design buckling resistance Nb,Rd for torsional or torsional-flexural buckling (according to EN 1993-
1-1) shall be obtained using buckling curve b, and with relative slenderness given by:

‰
̅ = ª u=


 = min ( ,… ,  ,…­ )

46
Advanced Concept Training

1 ž² I‡
 =  (¯¡; + )

4® H" © …

1 ž² I‡
 = (¯¡; + )

4® H" © …

H" = H + H] + "

1
− –b

 = ob +  ,… d +  ,… d − 4u ,  ,… q
,…­
2u , ,

ž² 
 =
© 
,
° ±
H
" 
u =1−° ±
H"

with βA the ratio Aeff/A


fyb the basic yield strength
σcr the critical stress
σcr,T the elastic critical stress for torsional buckling
σcr,TF the elastic critical stress for torsional-flexural buckling
G the shear modulus
E the modulus of elasticity
IT the torsion constant of the gross section
CM the warping constant
iy the radius of gyration about yy-axis
iz the radius of gyration about zz-axis
lT the buckling length of the member for torsional buckling (= Lyz)
y0 the position of the shear centre
ly the buckling length for flexural buckling about the yy-axis

47
Steel

Example: Flexural Torsional buckling.esa


Look at the steel code check for member B1:

Lateral Torsional Buckling


General

General case
The flexural buckling check will be executed following EN 1993-1-1 art. 6.3.2.

C9:
≤1
C‰,<:


Where
C‰,<: = Œ›… W
XY/

W = W P, for class 1 or 2 cross-sections


W = W5P, for class 3 cross-sections
W = W566 for class 4 cross-sections

The reduction factor χLT will be calculated as follows:



Œ›… =
/
but χ < 1,0
²³ Ž–‘ ‘
´µ v¶´µ

Φ›… = 0,5“1 + α›… bλ ›… − 0,2d + λ ¹º •


with

W 
λ ›… = ª
C

Mcr Critical bending moment (see next paragraph for this calculation)

α Imperfection depending on the buckling curves:

48
Advanced Concept Training

With the following recommended buckling curves for lateral torsional buckling:

›…
For slendernesses λ ≤ λ ›…,» of for YYij ≤ 됟m,¼ ² lateral torsional buckling effects may be ignored.
—˜

Following EN 1993-1-1 article 6.3.2.3: λ ›…,» = 0.4 but this can be adapted in the national annex of a
country and also in Scia Engineer:

Lateral torsional buckling for rolled sections or equivalent welded sections

In the EN 1993-1-1 a distinction is made between lateral torsional buckling for general cases and for
rolled sections or equivalent welded sections.
This distinction can also be chosen in Scia Engineer through “Steel -> Beams -> Steel setup” this can
be chosen for all the beams:

49
Steel

With this option the reduction factor for Lateral Torsional buckling is calculated a bit differently:

1
Œ›… =
Φ›… + –Φ¹º

− βλ ¹º

Œ›… ≤ 1,0
Œ›… ≤ ‘
But:
/
¶´µ

Φ›… = 0,5“1 + α›… bλ ›… − λ ›…," d + βλ ¹º •

And

λ ›…,» = 0,4 (maximum value)

u= 0,75 (minimum value)

The reduction factor Œ›… may be modified as follows:

Œ›…
Œ›…,‡¾: =

f may be defined in the National Annex.

The following minimum values are recommended:

 = 1 − 0,5(1 − ' )[1 − 2,0bλ ›… − 0,8d ]




But f < 1,0

50
Advanced Concept Training

With kc by default taken from the next table:

But alternatively kc can also be calculated from the factor C1:

' = 1/I/

In Scia Engineer the user can choose between the standard method or the calculation of kc in function
of C1 (by default kc will be taken from Table 6.6) in “Steel -> Beams -> Steel setup”:

51
Steel

And also the values for beta and LambdaLT,0 can be adapted in Scia Engineer:

Lateral Torsional Buckling Check in Scia Engineer


This LTB check will also be executed in Scia Engineer.

Example: Industrial hall.esa

Consider column B6:

The default method for the calculation of the C-factors has been used. This is also indicated in the
preview of the check results:

First the Lateral Torsional Buckling Check will be shown:

And afterwards the parameters for the calculation of Mcr will be shown.

Below this “Mcr Parameters”-window the calculation method for the C parameters will be indicated.

52
Advanced Concept Training

Calculation of Mcr
For I sections (symmetric and asymmetric), and Rectangular Hollow Sections (RHS), the elastic
critical moment for Lateral Torsional Buckling Mcr is calculated by the following formula:

ž   ¡] '  ¡n (' ∙ Ÿ) ¯¡;


¿ª° ±

C = I/ + + bI _® − Ic _À d − “I _® − Ic _À •Á
(' ∙ Ÿ)² 'n ¡] ž  ¡]


Where

E the Youn modulus of elasticity (E = 210000 N/mm² for steel)


G the shear modulus (G = 80770 N/mm² for steel)
'∙Ÿ the lateral torsional buckling length of the beam between points which have
lateral restraint (= lLTB).
Iw the warping constant
It the torsional constant
Iz the moment of inertia about the minor axis
_® the distance between the point of load application and the shear center
'n A factor which refers to end warping. Unless special provision for warming fixity
is made, kw should be taken as 1,0.
_À _À = _Â − 0,5 Ã= ( ² + _²)
]
Ä4
l>

For doubly symmetric cross-sections: _À = 0


_Â the coordinate of the shear center

C1, C2 and C3 are factors depending on the loading and end restraint conditions.

In Scia Engineer are different methods implemented for the calculation of those C1 and C2 factors.
Those methods are explained further in this chapter.

Haunched sections (I+Ivar, Iw+Plvar, Iw+Iwvar, Iw+Ivar, I+Iwvar) and composed rail sections (Iw+rail,
Iwn+rail, I+rail, I+2PL+rail, I+PL+rail, I+2L+rail, I+Ud+rail) are considered as equivalent asymmetric I
sections.

For the other supported sections as defined above, the elastic critical moment for Lateral Torsional
Buckling Mcr is given by:

ž   ¡] ¡n (' ∙ Ÿ) ¯¡;


C = ª +
(' ∙ Ÿ)² ¡] ž   ¡]

with E the modulus of elasticity


G the shear modulus
'∙Ÿ the lateral torsional buckling length of the beam between points which
have lateral restraint (= lLTB).
Iw the warping constant
It the torsional constant
Iz the moment of inertia about the minor axis

53
Steel

LTBII
It is also possible to calculate Mcr with a more precise calculation, a second order Lateral Torsional
Buckling calculation. This will be explained further in this chapter.

General – calculation of C1, C2 and C3 factors


The coefficients C1, C2 and C3 can be calculated in Scia Engineer according to three different
methods:
o ENV 1993-1-1 Annex F
o ECCS 119/Galea
o Lopez, Young, Serna

By default the method according to ECCS 119/Galea is applied. The following paragraphs give more
information on these methods.

The user can choose between those 3 methods in “Steek -> Beams -> Steel setup”:

ENV 1993-1-1 Annex F

When this setting is chosen, the moment factors are determined according to ENV 1993-1-1 Annex
F Ref.[5].

For determining the moment factors (EN 1993-1-1: C1, C2 and C3) for lateral torsional buckling (LTB),
we use the standard tables.
The current moment distribution is compared with some standard moment distributions. This standard
moment distributions are moment lines generated by a distributed q load, a nodal F load, or where the
moment line reach a maximum at the start or at the end of the beam.

The standard moment distributions which is closest to the current moment distribution, is taken for the
calculation of the factors C1, C2 and C3.

54
Advanced Concept Training

ECCS 119/Galea

When this setting is chosen, the moment factors are determined according to ECCS 119 Annex B
Ref.[34].
The figures given in this reference for C1 and C2 in case of combined loading originate from Ref.[28]
which in fact also gives the tabulated values of those figures as well as an extended range.

The actual moment distribution is compared with several standard moment distributions. These
standard moment distributions are moment lines generated by a distributed q load, a nodal F load, or
where the moment line is maximum at the start or at the end of the beam.

The standard moment distribution which is closest to the actual moment distribution, is taken for the
calculation of the factors C1 and C2.

Galea gives results only for C1 and C2 factors.


C3 is taken from ECCS 119 Annex B tables 63 and 64. The C3 is determined based on the case of
which the C1 value most closely matches the table value.

Lopez, Yong, Serna

When this method is chosen, the moment factors are determined according to Lopez, Yong, Serna
Ref[35].

When using this method the coefficient C2 and C3 are set to zero.

The coefficient C1 is calculated as follows:


b1 − √'d b1 − √'d
ª√'4/ + o 4 q + 4
2 2
I/ =
4/

Where:

' = '/ '

C‡Å†

+ t/ C/ + t C + tc Cc + tO CO + t# C#
4/ =
(1 + t/ + t + tc + tO + t# )C‡Å†

C/ + 2C + 3Cc + 4CO + 5C#


4 = Æ Æ
9C‡Å†

t/ = 1 − '

'/c
t = 5
'

1 1
tc = 5 ° + ±
'/ '

'c
tO = 5
'/

t/ = 1 − '/

55
Steel

With:
k1 Taken equal to kw
k2 Taken equal to kw
M1, M2, M3, M4, M5 The moments My determined on the buckling system in the given sections as
shown on the above figure. These moments are determined by dividing the beam
into 10 parts (11 sections) and interpolating between these sections.
Mmax The maximal moment My along the LTB system.

This method is only supported in case both k and kw equal 0.50 or 1.00.

Comparison of the 3 calculation methods

In below an example in Scia Engineer in which the three methods are calculated:

Example: Cfactors.esa

In the steel setup the chosen calculation method for the C factors has been changed. In below an
overview of the results for those factors for the three methods:

ENV 1993-1-1 Annex F ECCS 119/Galea Lopez, Yong, Serna

Mcr = 1576.03 kNm Mcr = 1118.50 kNm Mcr = 842.64 kNm

There can be a big difference between the three calculation models.


In the method following “Lopez, Yong, Serna” the values for C2 and C3 are always taken equal to zero.
When comparing the C1 factors, the method following “ECCS 119/Galea” and “Lopez, Yong, Serna”
are approximately the same (1.26 and 1.20 respectively), but the C1 factor following the “ENV 1993-1-
1 Annex F” results in total different value: 1.89.

LTB Restraints

56
Advanced Concept Training

In Scia Engineer Lateral Torsional Buckling restraints can be inputted. Those restraints will change the
Lateral Torsional Buckling Length, used for the calculation of Mcr:

ž   ¡] '  ¡n (Ÿ›…Ç ) ¯¡;


¿ª° ±

C = I/ + + bI _® − Ic _À d − “I _® − Ic _À •Á
(Ÿ›…Ç )² 'n ¡] ž  ¡]


Where

'∙Ÿ the lateral torsional buckling length of the beam between points which have
lateral restraint (= lLTB).

This length will be taken as the distance between two LTB restraints.

Fixed LTB restraints are defined on top flange or on bottom flange. The LTB lengths for the
compressed flange are taken as distance between these restraints. The LTB moments factors are
calculated between these restraints.

The restraints can be inputted via “Steel -> Beams -> Member Check data -> LTB Restraints”.

And only the restraints on the compressed side are taken into account.

Example: LTB_Restraints.esa
In this example the same beam has been inputted three times:
Those are the results for Load Case LC1:

57
Steel

B1 B2

λ ›…,» = 0,4
And for beam B3, 6 LTB restraints are inputted, so no LTB calculation will be executed:
If λ ›… ≤ λ ›…,» or for ij ≤ λ ›…,» ² lateral torsional buckling effects may be ignored. With:
Y
Y—˜

So for beam B3 the following LTB check is displayed:

When looking at load case LC2, the top side of the beam will be under tension, so Scia Engineer
will not take into account the effects of the LTB Restraints:

Diaphragms
When diaphragms (steel sheeting) are used, the torsional constant It is adapted for symmetric or
asymmetric I sections, channel sections, Z sections, cold formed U, C and Z sections.

58
Advanced Concept Training

A diaphragm can be inputted in Scia Engineer via “Steel -> Beams -> Member Check data ->
Diaphragms”:

The settings for the diaphragm are:

Diaphragms LIB With this option the user can choose between diaphragms of
the library of Scia Engineer or input his own diaphragm.
Position z A diaphragm is only taken into account at the
compression side of the beam for the LTB check.
k The value of coefficient k depends on the number of spans of
the diaphragm:
k = 2 for 1 or 2 spans,
k = 4 for 3 or more spans.
Diaphragm position The position of the diaphragm may be either positive or
negative.
Positive means that the diaphragm is assembled in a way so
that the width is greater at the top side.
Negative means that the diaphragm is assembled in a way so
that the width is greater at the bottom side.
Bolt position Bolts may be located either at the top or bottom side of the
diaphragm.

59
Steel

Bold pitch Bolts may be either:


- in every rib (i.e. "br"),
- in each second rib (i.e. "2 br").
Frame distance The distance of frames
Length The length of the diaphragm (shear field.)

The diaphragm will only have an influence on the torsional constant It and will taken into account in the
calculation for Mcr for the Lateral Torsional Buckling check.

The torsional constant It is adapted with the stiffness of the diaphragms:

The value for It equals the previous value for It (thus It of the beam) plus a supplementary stiffness
calculated with the values of the diaphragm:

©²
¡;,ˆ: = ¡; + ÈÉÊℎ IË
ž²¯

= + +
Where:
/ / / /
̾ Í £Î £Î@,Ï £Î =,š £ÎÐ,Ï

  ¡566
IËY,š = '
Ñ
$Å 
IË=,š = I/"" Ò Ó H $Å ≤ 125
100


IË=,š = 1.25 I/"" Ò Ó H 125 < $Å < 200
100

3 ∙   ∙ ¡Â
IËÕ,š ≈
(ℎ − %)

ѳ
¡Â =
12

with l the LTB length


G the shear modulus
vorh the actual rotational stiffness of diaphragm

CθM, the rotational stiffness of the diaphragm
k
CθA,k the rotational stiffness of the connection between the diaphragm and
the beam
CθP,k the rotational stiffness due to the distortion of the beam
k numerical coefficient
= 2 for single or two spans of the diaphragm
= 4 for 3 or more spans of the diaphragm
EIeff bending stiffness of per unit width of the diaphragm
s spacing of the beam
ba the width of the beam flange (in mm)
C100 rotation coefficient - see table
h beam height
t thickness beam flange
s thickness beam web

60
Advanced Concept Training

Example: Diaphragm.esa

Consider member B1:


o LTB Length = 7 m
o C1 = 1.13, C2= 0.45, C3 = 0.53
o Mcr = 282 kNm

Consider member B2
4
o It,id = 10880887 mm
o LTB Length = 7 m
o C1 = 1.13, C2= 0.45, C3 = 0.53
o Mcr = 332 kNm

The results for the calculation of the diaphragm are shown in the preview of the steel check under the
Lateral Torsional Buckling check:

61
Steel

Lateral Torsional Buckling using LTBII – Not in Concept Edition


This option has been inputted in Scia Engineer as a separate module and it is not included in an
edition. The necessary module for this option is esasd.14.

For a detailed Lateral Torsional Buckling analysis, a link was made to the Friedrich + Lochner LTBII
calculation program.
nd
LTBII is the abbreviation of “Lateral Torsional Buckling with 2 order calculation”.

The Frilo LTBII solver can be used in 2 separate ways:

1. Calculation of Mcr through eigenvalue solution


nd
2. 2 Order calculation including torsional and warping effects

For both methods, the member under consideration is sent to the Frilo LTBII solver and the respective
results are sent back to Scia Engineer.

The single element is taken out of the structure and considered as a single beam, with:
o Appropriate end conditions for torsion and warping
o End and begin forces
o Loadings
o Intermediate restraints (diaphragms, LTB restraints)
nd
To use this option in Scia Engineer the functionality “LTB 2 order” should be activated.
And afterward with the option “Steel -> Beams -> Member check data -> LTB II member data” the user
can input LTB-data on a beam.

Example: LTBII.esa
In this example the same beam will be calculated with LTBII data on it and without in three
configurations:

When looking at the check for LC1 there is not a big difference between the beams with or without the
LTBII data:

62
Advanced Concept Training

Only for the calculation with the diaphragm the results for the calculation with LTBII data are favourable.
In Scia Engineer the calculation of a diaphragm is only calculation in one direction.
In Frilo the rigidity of the diaphragm against rotation and translation is taken into account, which results
in a higher stiffness of the diaphragm and thus a better unity check.

When looking at the results for Load Case LC2, the direction of the inputted line load is changed. The
upper side of the beams is under tension for this load case. Scia Engineer does take into account LTB
restraints or a diaphragm on the tension side. But Frilo can take this also into account. The restraints
and the diaphragm will have a lower influence than at the compression side, but there will be an
augmentation of the rigidity:

Also in the detailed output this calculation is indicated:

63
Steel

Compression and bending check


The compression and bending check for a member will be executed following EN 1993-1-1 art. 6.3.2.

Unless second order analysis is carried out, members which are subjected to combined bending and
axial compression should satisfy:

8 C,9: + ∆C,9: C],9: + ∆C],9:


Œ 8<š + ' + '] ≤1
9:
C,<š C],<š
ŒŸm X XY/
XY/ Y/

8 C,9: + ∆C,9: C],9: + ∆C],9:


Œ_8<š + '] + ']] ≤1
9:
C,<š C],<š
XY/
ŒŸm X XY/
Y/

Øz{ , ÙÚ,z{ , ÙÛ,z{ are the design values of the compression force and the maximum moments about
the y-y and z-z axis along the member respectively.

ÜÚ and ÜÛ are the reduction factors due to flexural buckling

Ü݁ is the reduction factor due to lateral torsional buckling

The interaction factors kyy, kyz, kzy and kzz have been derived from EN 1993-1-1 Annex A (alternative
method1) or from Annex B (alternative method 2).

The choice between Interaction Method 1 or 2 can be made in Scia Engineer in the National Annex
parameters:

64
Advanced Concept Training

Example: Industrial Hall.esa

Consider column B28 (for combination CO1-ULS).

In Scia Engineer first all calculated formulas, in this example following EN 1993 Annex 1, are given:

And afterwards the check is given as two times a sum of three values as in formulas (6.61) and (6.62)
of the EN 1993-1-1:

65
Steel

89: C,9: + ∆C,9: C],9: + ∆C],9:


+ ' + '] ≤1
Œ 8<š C C],<š
ŒŸm X,<š XY/
XY/ Y/

8 C,9: + ∆C,9: C],9: + ∆C],9:


Œ_8<š + '] + ']] ≤1
9:
C,<š C],<š
XY/
ŒŸm X XY/
Y/

In Scia Engineer:

Shear buckling check – EN 1993-1-5


General
The shear buckling check will check if the web of the cross section can locally buckle.
This check is not included in EN 1993-1-1, but in 2N 1993-1-5: Design of steel structures – Part 1-5:
Plated structural elements.

First it will be checked if the slenderness is lesser than a certain value, because for webs with a small
slenderness, this check should not be executed.

For plates with hw/t greater than 72ε/η for an unstiffened web should be checked for resistance to
shear.

With:

Þ=–
c#
6> [p/‡‡²]

η will be defined in the National Annex. The value η = 1.20 is recommended for steel grades
up to and including S460. For higher steel grades η = 1.00 is recommended.

If the slenderness is higher than the minimum one, the shear buckling check is executed following EN
1993-1-5, formula (5.10) (this is a check on the shear force) and formula (7.1) (which is the check on
the interaction between shear force, bending moment and axial force). Both formulas are checked in
Scia Engineer.

66
Advanced Concept Training

Example: Industrial Hall.esa

Consider Beam B50

This is an IPE 160 profile: hw = 160 – 2 x 7.40 = 145.2 mm. And t = 5 mm


 hw / t = 29.04

This should be checked with the value: 72ε/η = 72 x 1.00 / 1.2 = 60


 hw / t < 72ε/η
 The shear buckling check does not have to be executed.

This is also indicated in Scia Engineer:

Consider Beam B28 (for combination CO1-ULS):

This is an IPE750x137 profile: hw = 753 – 2 x 17 = 719 mm. And t = 11.5 mm


 hw / t = 62.52

This should be checked with the value: 72ε/η = 72 x 1.00 / 1.2 = 60


 hw / t > 72ε/η
 The shear buckling check has to be executed.

67
Steel

Stiffners

In EN 1993-1-5, Figure 5.1 the definitions of stiffeners on a member are given:

Those stiffeners will influence the total length for shear buckling. Also in Scia Engineer those stiffeners
can be inputted on a beam and will only have influence on the shear buckling check.

Those stiffeners can be inputted via “Steel -> Beams -> Member check data -> Stiffeners”.

The stiffeners define the field dimensions (a,d) which are only relevant for the shear buckling check.
When no stiffeners are defined, the value for ‘a’ is taken equal to the member length.

Example: Stiffeners.esa

In this example 2 identical beams are inputted. B1 without stiffeners and B2 with stiffeners:

Let’s now compare the shear buckling check with and without stiffeners:

68
Advanced Concept Training

Without stiffeners With stiffeners


Shear buckling check: Shear buckling check:

69
Steel

ULS Check for Battened compression members

This check can be executed following the EC-ENV and EC-EN 11993-1-1
art.6.4

The following section pairs are supported as battened compression member :

(1) 2I

(2) 2Uo

(3) 2Uc

Two links (battens) are used.

The following additional checks are performed:


- buckling resistance check around weak axis of single chord with Nf,Sd
- section check of single chord, using internal forces
- section check of single batten, using the internal forces
ho

Example – manual calculation of a battened compression member


EC_EN_Battened_Compression_Members_I.esa

Consider member B1
On this beam, only a compression force of 500kN is applied
Properties of the batten:

70
Advanced Concept Training

Properties of the IPE 330 profile

Values in Scia Engineer:

Manual calculation of those values:


l length = Member length – (distance from begin) – (distance from end)
= 5000mm – 200mm – 200mm
= 4600mm
a distance between battens = l / (number of divisions) = 4600mm / 4 = 1150 mm
h0 distance between centroids of chords = 210mm
Ich Iz of the IPE330 profile
Ach area if one I-profile, so A of IPE 330
Lambda = l / i0 = 4600mm / 110,833 mm = 41.5
with
©/ = 0.5 ∙ ℎ0 ∙ 4Lℎ + 2 ∙ ¡Lℎ
= 0.5 ∙ (210BB) ∙ (6260BB ) + 2 ∙ (7881000BBO ) = 1.538  10^ BBO

H0 = – =– = 110,833BB
P‹ /.#c^ † /"ß ‡‡§
∙=Í ∙ZZ"‡‡‘

Mu à:
= 0 if Lambda > 150
= 2−
›Å‡‰:Å
¦#
if 75 < Lambda < 150

= 1.0 if Lambda < 75

71
Steel

And Lambda = 41.5 < 75 in this example => à = 1.00


leff 0.5 ∙ ℎ0 ∙ 4Lℎ + 2 ∙ à ∙ ¡Lℎ
= 0.5 ∙ (210BB)² ∙ (6260BB ) + 2 ∙ (1.00) ∙ (7881000BBO ) = 1.538  10^ BBO
∙™‘ ∙9∙lÍ
= $å% æÈ ≤
O∙9∙lÍ
‘∙â—ã ãA ő
ő ∙f/Ž ∙ g
Sv
V∙âä ¥

With:
n = number of planes of battens: n=2
Ib = (tickness batten) x (width of batten)³ / 12 = (7 mm) x (220 mm)³ / 12
6 4
= 6,211 x 10 mm
‘‹AAAA¤
O∙f g∙(¦^^/"""‡‡§ )
= = 2.438  10O '8
O∙9∙lÍ TT‘
‘∙â—ã ãA ‘∙(çßߋAAATT§ ) ‘‹ATT
ő ∙f/Ž ∙ g (//#"‡‡)‘ ∙°/Ž ∙ ±
V∙âä ¥ ‘∙(è,‘‹‹ é ‹Aè êê§ ) ‹‹ëATT

And
‘‹AAAA¤
∙™‘ ∙9∙lÍ ∙™‘ ∙f g∙(¦^^/"""‡‡§ )
= TT‘
= 2.470  10O '8
ő (//#")‘

 æÈ = 2.438  10O '8

e0 = l/500 = 4600mm/500 = 9.2 mm


MEd,I Mz in the cross section under consideration = 0 kNm
MEd will be calculated as follows:

If f1 − − g>0 => C Ä =


p9: p9: p9:∙5"ŽY9:_l
p kÌ ¤ij ¤ij
/v v
¤—˜ ía
6
Otherwise: = > Med = 1 x 10 kNm

With:
‘‹AAAA¤
™²∙f g∙(/.#c^ † /"ß ‡‡§ )
8 = = = 15064 '8
™²∙9∙l566 TT‘
›² (OZ""‡‡)²

1− − = 0,9466 > 0
#""šp #""šp
/#"ZOšp O¦"" šp

 C Ä = = = 4860'8BB = 4,860'8B


p9:∙5"ŽY9:_l #""šp∙\.‡‡Ž"šp‡
¤ij ¤ij ".\OZZ
/v v
¤—˜ ía

= = 3,32'8
™∙Y9: ™∙O^Z"šp‡‡
P OZ""‡‡
VEd

Check of chord as beam in field between battens

The check of the chord will be executed with a section check following the EN 1993-1-1 for the chord
profile with the following internal forces on one batten:

NG = NChord = 0.5 ∙ 8 Ä + = 0.5 ∙ 500'8 + = 270,78 '8


Y5:∙Í"∙=Í (O^Z"šp‡‡)∙(/"‡‡)∙(ZZ"‡‡²)
∙P566 ∙/.#c^ † /"ß ‡‡§

VG = VEd /2 = 3,32kN / 2 = 1,66kN


= = 0.95'8B
h5:∙Å c,cšp∙//#"‡‡
O O
MG =

72
Advanced Concept Training

For an I-section a classification is made, for an U-section always an elastic check is done.
This section check will result in a unity check of 0,19

Buckling check of chord

Nch,Ed = NChord = 170,78 kN (see check of Check of chord as beam in field between battens)
Buckling length = a

= = 32,411
ÇîšPˆw® P5w®;Í //#"‡‡
§
– â—ã –çßߋAAATT
Slenderness =
ï—ã è‘èATT²

‘çë¤


ZZ"‡‡‘ ∙
Relative slenderness Lambda – = 0.373
=Í∙6 TT‘
p ð /c#/""" p

With:
‘‹AAAA¤
™²∙f g∙(¦^^/"""‡‡§ )
8LÊ] = = = 12351 '8
™²∙9∙lÍ TT‘
(ÇîšPˆw® P5w®;Í)² (//#"‡‡)²

Buckling curve = b, this is the buckling curve for the IPE 330 around the z-axis
Imperfection alpha for buckling curve b = 0.34
Reduction factor Chi can be calculated as explained in the buckling check.

= = 0.17
pÍ9: ¦",¦^šp
p‰<: /Z/^, šp
Unity check =

With:

8$, ñÄ = ÊKÄåL%HÉG EL%ÉÊ IℎH = 0.94 = 1618,21 '8


=Í ∙6 ZZ"‡‡² ∙¦#p/‡‡²
?@‹ /.""

Check of batten:

t thickness of batten
b width of batten
= = = 9.09'8
h9:∙Å c,cšp∙//#"‡‡
Í"∙ /"‡‡∙
T

73
Steel

MG = 0.954 kNm (see previous check)

= = 16,9CDE
Yò∙‰/ \#O"""p‡‡∙"‡‡/
l‰ Z.//∙/"臇§
Sigma

(calculation of Ib see calculation of Sv at the properties)

= = 0.06
kˆ®‡Å /Z,\YÕÅ
Unity check sigma: = N> ‘çë¤/TT²
ó@A ‹.AA

= = 5.9CDE
… \"\"p
‰∙; "‡‡∙¦‡‡
Tau =

= = 0.04
…Åî #,\YÕÅ
N>/√ô ‘çë¤/TT²/√ô
Unity check Tau: =
ó@A ‹.AA

74
Advanced Concept Training

ULS Check for Built-in beams

The built-in beams are of the types SFB (Slim Floor Beam), IFB (Integrated Floor Beam) or THQ (Top
Hat Beams).
The relevant section classification according to the proper code is performed. Classification for class 1,
2, 3 and 4 is supported. For the ULS section check, the local plate bending is taken into account for the
plastic moment capacity and the bending stresses in the section. Out-of-balance loading is checked.

This out-of-balance check has been inputted by a factor “Q”.

When inputting a line load on such a section, the option “Bottom flange” can be checked and
afterwards a Q factor (an out-of-balance factor) can be inputted in the properties window of this line
load. This Q factor indicates the distribution of the force:

If Q = 0, both sides are loaded equally, Q=1: one side is loaded the other side isn’t.

The position of the load is indicated in the cross section library by the value sw:

75
Steel

Example: Built-in beams.esa

In this example a beam with an IFB profile has been inputted with a line load of 100kN/m and an out-of-
balance factor Q=0.8:

In the steel code check, the out-of-balance check is given:

76
Advanced Concept Training

77
Steel

Optimisation
In Scia Engineer there are two ways of performing an optimization:
1. Cross section optimisation: Optimisation of a chosen cross section.
2. Overall optimisation: Optimisation of one or more (or all) cross sections at the same time.

Both options are explained in the next two paragraphs.

Cross section optimisation


With this option the user can optimise a cross section for the steel code check. This will be explained in
the example below.

Example: Optimisation.esa

The optimisation can be found via “Steel -> Beams -> ULS Checks -> Check”. First the steel code
check has to be calculated, so first click on “Refresh”.

When clicking now on “Autodesign”, Scia Engineer will give the following message:

And in the properties window the option for “Filter” has been changed into “Cross-section”. Now choose
here for the profile IPE 140 and click again on “Refresh” and afterwards on “Autodesign”.

Now the “Autodesign of the cross-section” window will open:

78
Advanced Concept Training

First you can fill in the maximal check. Normally this is a unity check of 1.00.
The Maximum unity check displayed below is the maximum unity check for the IPE140, found in this
project.

When clicking on “Search for optimal”, Scia Engineer will suggest the smallest IPE profile which will
resist the calculated internal forces.

In this example an IPE 200 will be found. When clicking on “OK” now, Scia Engineer will replace this
IPE140 profile by an IPE200 profile automatically.

Now the steel code check of this IPE200 profile is executed with the internal forces calculated with the
properties of the IPE140 profile. For example the self weight will not be taken into account correctly. So
the project should be recalculated, before accepting those new unity check results.

It is also possible to use a list of sections.


With this list it is possible to indicate which sections can be used or not. For example the section
IPE140A can be filtered out the possibilities for the autodesign.

Recalculate the project and let’s make a list for the columns. This can be done through “Library ->
Structure, analysis, Cross-section list” and choose for “Rolled cross-section one type”.

Now select the profiles that should be added to this list:

When performing an Autodesign in the steel menu now for the columns HEA220, this list can be
chosen:

79
Steel

For the moment this profile has a maximum unity check of 0.368, so Scia Engineer will search for the
smallest profile from this “LIST1” which will pass the unity check

Close this example without saving! It will be used in the next paragraph again!

In the next example the different options of the cross-section optimisation are explained.

Example: Optimisation2.esa

Go here to “Steel -> Beams -> ULS checks -> Check” and go to the Autodesign menu for this beam.

This profile has a lot of parameters, so the user can choose which parameters should be updated:

80
Advanced Concept Training

Also the option “Advanced Autodesign” will appear. With this option several parameters can be
optimised as the same time and restriction for each parameter can be given. So input the following
options:

So the thicknesses thb and tha are not autodesigned. All other options can be adapted by Scia
Engineer. When clicking on “Search for an optimal”, the following profile is suggested by Scia
Engineer:

81
Steel

With the option “Set value” a value for a certain parameter can be set. So select the parameter “Ba”
and click on “Set value”. Now a value of 500mm can be inputted. And the Maximum unity check for this
profile will be adapted automatically.

Overall optimisation
It is also possible to perform an overall optimisation in Scia Engineer. With this option, one or more
profiles can be optimised at the same time. Afterwards the calculation will be restarted and the internal
forces are recalculated with the new cross-sections, followed again with a new optimisation.

This iterative process can


o Or stop because all profiles does not have to be autodesigned and the same profile was
found as in the previous step
o Or stop because the maximum number of iteration steps has been reached if this is inputted
by the user.

It is advised to input a number of iteration steps, otherwise this optimisation process can become a
loop and will stop after 99 iteration steps. This will cost a lot of calculation time.

The principle of the overall optimisation process is explained by the following example.

Example: Optimisation.esa

If this example is still open from the previous chapter, please close it and reopen it without saving.

Calculate the project and go to “Main -> Calculation, Mesh -> Autodesign”.
Click here on Add item and choose for “Cross-section steel check” and add all cross-sections in this
Autodesign process.

The combination, for which the optimisation calculation has to be calculated, can be chosen.

82
Advanced Concept Training

Also for each profile, the user can indicate if this profile has to be autodesigned only in the “Up”
direction (so only become bigger) or in the ‘Up and down” direction (and also can become a smaller
profile as result). With this last option, there is a chance that the iterative process will become a loop.

Now Close this window. And click in the next window on Opimisation Routine:

Now the user can choose for a maximum number of iteration steps:
o Determine automatically: the iteration process will stop if all the profiles are autodesigned
and no different result will be found in a certain iteration step. So no maximum of number of
iteration steps will be inputted.
o Limit number of iteration: maximum of number of iteration steps.

83
Steel

Input 4 iteration steps as limit number and click on “Start”.

Now the iteration process will start. The iteration process will stop after no difference will be found or
after 4 iteration steps. In this example Scia Engineer will stop after 3 iteration steps, because all
profiles remain the same after the third step.

After this process an information window about this iteration will be displayed.

ITERATION STEP 1:

Changed profiles:
o CS1: HEA220 - > HEA200A
o CS2: IPE140 -> IPE200
o CS3: HFLeq45x45x5 -> HFLeq80x80x10
o CS4: IPE100 -> IPE140

ITERATION STEP 2:

Changed profiles:
o CS1: HEA220A - > HEA160
o CS3: HFLeq80x80x10 -> HFLeq90x90x9
o CS4: IPE140 -> IPE140A

ITERATION STEP 3:

No profile is adapted. So the iteration process is stopped.

84
Advanced Concept Training

2nd order calculation and imperfections

Overview
Global analysis aims at determining the distribution of the internal forces and moments and the
corresponding displacements in a structure subjected to a specified loading.
The first important distinction that can be made between the methods of analysis is the one that
separates elastic and plastic methods. Plastic analysis is subjected to some restrictions.
Another important distinction is between the methods, which make allowance for, and those, which
neglect the effects of the actual, displaced configuration of the structure. They are referred to
respectively as second-order theory and first-order theory based methods.

The second-order theory can be adopted in all cases, while first-order theory may be used only when
the displacement effects on the structural behavior are negligible.

The second-order effects are made up of a local or member second-order effects, referred to as the P-
δ effect, and a global second-order effect, referred to as the P-∆ effect.

On the next page an overview of the global analysis following the EN 1993-1-1, chapter 5, will be
given:
o All the rules in this overview are given in the EN 1993-1-1 art. 5. For each step the rule will
be indicated. The first rule (αcr > 10) will be explained in EN 1993-1-1 art. 5.2.1(3).
o In this overview 3 paths are defined:
 Path 1: In this path a first order calculation will be executed
 Path 2: In this path a second order calculation will be executed with global (and
bow) imperfections.
 Path 3: In this path a second order calculation will be executed with the buckling
shape of the construction as imperfection.
o The calculation will become more precise when choosing for a higher path.
o The lower paths will result in a faster calculation, because a first order calculation can be
executed without iterations, but this first-order theory may be used only when the
displacement effects on the structural behavior are negligible.
o In the next paragraphs the rules in this overview will be explained.

85
Steel

Structural Frame Stability

5.2.1(3)

αcr ≥ 10
Yes No
5.2.2(3)b 5.2.2(3)a 5.3.2(11)

Global Imperfection φ

5.3.2(6)

No Yes
αcr ≥ 3
5.2.2(5)
NEd > 25% Ncr ηcr
(member)

e0 if e0 in all
required members

1a 1b 2a 2b 2c 3

1st Order Analysis 2nd Order Analysis

Increase sway
effects with: No Members Yes
1 with e0
1
1−
α cr
5.2.2(3)c 5.2.2(7)b

lb based on a global lb taken equal to L


buckling mode

Stability Check in plane


5.2.2(7)a

Stability Check out of plane + LTB Check

Section Check

With: ηcr Elastic critical buckling mode.


L Member system length
lb Buckling Length

86
Advanced Concept Training

Alpha critical – Not in concept edition


The calculation of alpha critical is done by a stability calculation in Scia Engineer. For this calculation a
Professional or an Expert edition is necessary. The stability calculation has been inputted in module
esas.13.
st
According to the EN 1993-1-1, 1 Order analysis may be used for a structure, if the increase of the
relevant internal forces or moments or any other change of structural behaviour caused by
deformations can be neglected. This condition may be assumed to be fulfilled, if the following criterion

t = ≥ 10 for elastic analysis


­—˜
is satisfied:

­ij

With: αcr: the factor by which the design loading has to be increased
to cause elastic instability in a global mode.
FEd: the design loading on the structure.
Fcr: the elastic critical buckling load for global instability,
based on initial elastic stiffnesses.
nd
If αcr has a value lower then 10, a 2 Order calculation needs to be executed. Depending on the type
of analysis, both Global and Local imperfections need to be considered.
nd
EN1993-1-1 prescribes that 2 Order effects and imperfections may be accounted for both by the
global analysis or partially by the global analysis and partially through individual stability checks of
members.

Global frame imperfection ϕ


The global frame imperfection will be inputted for the whole structure by an imperfection value φ. This
value can be calculated with the following formula (EN 1993-1-1 art. 5.3.2(3)a):

1
õ= ∙t ∙t
200 Í ‡
tÍ = ≤ tÍ ≤ 1,0
 
√Í c
but

1
t‡ = ª0,5 °1 + ±
B

With: h The height of the structure in meters


m The number of columns in a row including only those columns which carry a vertical
load NEd not less than 50% of the average value of the vertical load per column in
the plane considered.

The global imperfection can be inputted in Scia Engineer in the non linear combinations as explained
by the example below.

87
Steel

Example: Industrial Hall.esa

In this example, two global imperfections functions will be used: one according to global X-direction
and one according to the global Y direction. It is not necessary to combine both imperfections in the
same combination.
In Scia Engineer the user can choose between four options for the global imperfection in the non
linear combinations:

o None: No global imperfection will be taken into account


o Simple inclination: the user will input a deformation in the global X- or/and Y-direction as
the imperfection in mm per meter height in the global Z-direction [mm/m].
o Inclination functions: input a user inclination function or the inclination function of the EN
1993-1-1 (this option is explained below).
o Deformation from loadcase: with this option the user can choose for a calculated loadcase
and the deformation of that loadcase will be used as initial global imperfection.

In this example the option “inclination functions” has been chosen.

These inclination functions are entered through “Main -> Library -> Structure, Analysis -> Initial
deformations”:

88
Advanced Concept Training

The Type is chosen as “According to code”, with a standard imperfection of 1/200.


The height of the construction is 8.4m for both inclination functions.
There are 6 columns in the X-direction, but in the middle only 2 columns are inputted. Because a long
part of structure only has 2 columns in the X-direction, in this example the number of columns in this
direction has been inputted as “2”.
There are 11 columns in the Y-direction. But the columns at the end are smaller than the middle ones.
So in this example it is decided to input “9” columns in the Y-direction.
The inclination function for the x-direction (Def_X) in Scia Engineer is displayed below:

In this example 5 non linear combinations are inputted:


o 1,35 x Self Weight + 1,35 x Self Weight Cladding + 1,5 x Maintenance
+ 0,75 x Snow + 0,9 x 3DWind10
o 1,35 x Self Weight + 1,35 x Self Weight Cladding + 0,75 x Snow + 1,5 x 3DWind13
o 1,35 x Self Weight + 1,35 x Self Weight Cladding + 0,75 x Snow + 1,5 x 3DWind14
o 1,35 x Self Weight + 1,35 x Self Weight Cladding + 0,75 x Snow + 1,5 x 3DWind15
o 1,35 x Self Weight + 1,35 x Self Weight Cladding + 0,75 x Snow + 1,5 x 3DWind16

89
Steel

All combinations are entered four times:


o NC1-NC5: inclination in the positive X- direction
o NC6-NC10: inclination in the negative X- direction
o NC11-NC15: inclination in the positive Y- direction
o NC16-NC20: inclination in the negative Y- direction

Bow imperfection
NEd > 25% Ncr
The relative initial local bow imperfections of members for flexural buckling is given as the value:
e0/L.

This bow imperfection does not have to be applied on each member as given in EN 1993-1-1 art.
5.3.2(6):

The bow imperfection has to be applied when the normal force NEd in a member is higher than 25%
of the member’s critical load Ncr:
When performing the global analysis for determining end forces and end moments to be used in
member checks according to 6.3 local imperfections may be neglected. However for frames sensitive
to second order effects local bow imperfections of members additionally to global sway imperfections
should be introduced in the structural analysis of the frame of each compressed member where the
following conditions are met:

• At least one moment resistant joint at one member end

̅ > 0,5 –
= ∙6>

pij

90
Advanced Concept Training

̅ = –
= ∙6>
p—˜
And

> 0,5 –
= ∙6> = ∙6>
Thus: –p pij
—˜

> 0,25
/ /
p—˜ pij

Øz{ > 0,25 Øö÷

Bow imperfection e0

The initial bow imperfection is given by:

Where L is the member length.

As explained in the previous chapter, the bow imperfection has to be


applied when the normal force NEd in a member is higher than 25% of
the member’s critical buckling load Ncr. If NEd < 25% Ncr, the user can
choose to apply this bow imperfection or not.

The buckling curve used for calculation of the imperfection is the curve inputted in the cross-section
library. For standard sections, the curve according to the code is automatically used, for non-
standard cross sections (as general cross sections) the user needs to input the buckling curve
manually.

Example: Industrial hall.esa

In Scia Engineer the user can choose between three options for the bow imperfection in the non
linear combinations:

91
Steel

o None: No bow imperfection will be taken into account


o Simple curvature: the user will input a deformation “f” or “1/f” for all the members, where
“1/f” corresponds with the value “e0/L” as explained earlier.
o According to buckling data: with this option the user can choose for a local imperfection as
defined in the buckling data, so for each member following his own buckling curve (this
option is explained below).

In this example the option “According to buckling data” has been chosen.
The buckling data can be inputted via “Steel -> Beams -> Member Check data -> Member buckling
data”

92
Advanced Concept Training

The user can choose between 6 options:


o According to code – elastic: the elastic value following the buckling curve of the cross
section will be used.
o According to code – plastic: the plastic value following the buckling curve of the cross
section will be used.
o According to code – elastic – only if required: the elastic value following the buckling
curve of the cross section will be used if Ned > 25% Ncr.
o According to code – plastic – only if required: the plastic value following the buckling
curve of the cross section will be used if Ned > 25% Ncr
o Manually input of bow imperfection: the user can input manually a value for this bow
imperfection.
o No bow imperfection: no bow imperfection is taken into account for the member.

In this example on all beams the bow imperfection “According to code – elastic” has been inputted.
Except for the diagonals no bow imperfection has been inputted in this example.

Buckling shape as imperfection - ηcr – Not in concept edition

The calculation of alpha critical is done by a stability calculation in Scia Engineer. For this calculation a
Professional or an Expert edition is necessary. The stability calculation has been inputted in module
esas.13.

As an alternative the global and bow imperfection can be replaced by the buckling shape as
imperfection (path 3 from the global diagram).

To input geometrical imperfections, the functionality Nonlinearity > Initial deformations and
curvature and Stability must be activated.

93
Steel

So first the stability calculation (linear or non linear) is calculated. Afterwards in the non linear
combinations, the user can choose the buckling shape he wants to take into account. So first he can
choose for the stability combination and just below for the calculated mode. If for example the first and
second mode give a negative value for the alpha critical, it could be possible that the user wants to use
rd
the 3 calculated buckling mode:

The last option that has to be inputted is a value for the maximum deformation. This is the deformation
of the node which has the biggest deformation of the construction. Scia Engineer will recalculate all the
displacement of the other nodes using this maximal deformation.

Since the buckling shape is dimensionless, EN 1993-1-1 gives the formula to calculate the amplitude
ηinit of the imperfection. In below an example of this calculation is given for a simple example.

Example: Buckling Shape.esa

In this example, the procedure to calculate a buckling shape is illustrated for a column.

The column has a cross-section of type IPE 300, is fabricated from S235 and has the following relevant
properties:

E = 210.000 N/mm²
fy = 235 N/mm²
γM1 = 1.00
L = 5000 mm
A = 5380 mm²
4
Iy = 83560000mm
3
W pl,y = 628400 mm

Calculation of the buckling shape

First a Stability calculation is done using a load of 1kN. This way, the elastic critical buckling load Ncr
is obtained.
In order to obtain precise results, the Number of 1D elements is set to 10. This can be done in “Setup
-> Mesh”. In addition, the Shear Force Deformation is neglected so the result can be checked by a
manual calculation.

94
Advanced Concept Training

The stability calculation gives the following result:

This can be verified with Euler’s formula using the member length as the buckling length:

ž² ¡ ž² ∙ 2100008/BB² ∙ 8,3560 ∙ 10¦ BBO


8 = = = 6927,51 '8
©² (5000BB)²

The following picture shows the mesh nodes of the column and the corresponding buckling shape:

th
Using for example an Excel worksheet, the buckling shape can be approximated by a 4 grade
polynomial.

95
Steel

A polynomial has the advantage that the second derivative can easily be calculated.

 ø = −2,11 ∙ 10v// ∙  O + 2,11 ∙ 10v¦ ∙  c − 6,95 ∙ 10v# ∙   − 2,29 ∙ 


 ø" ,‡Å† = −2,53 ∙ 10v/" ∙   + 1,27 ∙ 10vZ ∙  − 1,39 ∙ 10vO

Calculation of e0

8<š =  ∙ 4 = 235 ∙ 5380 BB = 1264300 8


p
‡‡‘

C<š =  ∙ W P = 235 ∙ 628400 BBc = 147674000 8BB (L©EÑÑ2)


p
‡‡‘

̅ = 8<š /8 = 12643008/68852808 = 0,43

t = 0,21 for buckling curve a

Œ= = 0,945
/
‘ ‘ ‘
v",dŽb
",#ú/Žrb d û Ž–f",#ú/Žrb d‘ ûg vb
v",dŽb  d‘

These intermediate results can be verified through SCIA ENGINEER when performing a Steel Code
Check on the column:

 d‘
ü∙bý
/v
 K" = t ∙ b ̅ − 0,2d ∙ ∙ = t ∙ b ̅ − 0,2d ∙
YxÏ ó@‹ YxÏ
‘
pxÏ d
/vŠ∙b pxÏ

 K" = 0,21 ∙ (0,43 − 0,2) ∙ = ~, ~þ ࢓࢓


/O¦Z¦O""" p‡‡
/ZOc"" p

The required parameters have now been calculated so in the final step the amplitude of the
imperfection can be determined.

Calculation of ηinit
The mid section of the column is decisive ⇒ x = 2500

ø at mid section = -3504

ø" ,‡Å†
1
-03
at mid section = 1,44E
mm²
 øˆwˆ; = K" ∙ 9∙l ∙ ø
p—˜
> ∙ఎ—˜,T¥ೣ
"

 øˆwˆ; = 5,57 BB ∙ ‘‹AAAA¤ ∙ 3504


Z^^#^" p
∙^c#Z""""‡‡§ ∙/,OO∙/"షô //‡‡²
TT‘

 øˆwˆ; = ~, ૜૝ ࢓࢓

96
Advanced Concept Training

This value can now be inputted as amplitude of the buckling shape for imperfection.

97
Steel

SLS Check

Nodal displacement
The nodal displacement defines the maximum global deflections in the vertical and horizontal
directions.
The following values are controlled in the example below:
o Limit for horizontal deflections δ is h/150
o Limit for vertical deflection δmax is L/200

Example: Industrial Hall.esa

Look at “Results -> Displacement of nodes” and look at the combination CO2 – SLS.

Horizontal deformation
The maximum displacement in the X direction is 25.3mm on a height of 6.9m.
And in the Y-direction 27.4mm on a height of 8.1m.

Limit for Limit for horizontal deflection δ is h/150


 6900/150 = 46 mm → 25.3mm < 46 mm → OK
 8100/150 = 54 mm → 27.4mm < 54mm → OK

Vertical deformation
The maximum displacement in the Z direction is 59.6 mm

Limit for Limit for vertical deflection δ is L/200


 60000/2000 = 300 mm → 59.6 mm < 300 mm → OK

Relative deformation
For each beam type, limiting values for the relative deflections are set, using the menu “Steel -> Beams
-> Setup -> Relative deformations”:

98
Advanced Concept Training

With the option “Steel -> Beams -> SLS Checks - Relative deformation”, the relative deformations can
be checked. The relative deformations are given as absolute value, relative value related to the span,
or as unity check related to the limit for the relative value to the span.

The span is defined in the menu “Buckling and relative lengths”.

Two options for this span are possible and explained below:

2 nodes supported

When two nodes are supported in this menu as shown below, the deformation is taken as the
maximum deflection of the beam in comparison with a line connecting the two end nodes:

In below an example on this principle:

99
Steel

Calculation for this relative deformation:

Deformation in the beginning of the beam = 5,967mm and at the end = 11,863mm.

Maximum deformation is at 0,979m from the beginning at the beam.


So this punt has already displaced 9,5mm (see blue line at the picture):

0,979B
å_‰Pî5 Pˆw5;",\¦\‡ = 5,967BB + ∗ (11,863BB − 5,967BB) = 9,504BB
1,632B

And uz, relative = 10,2mm – 9,5mm = 0,7mm

And suppose the length of this beam = 1632mm:

 ÊK© å_ = = 1/2236
",¦c‡‡
/Zc‡‡

1 node supported and the other free

This is the case if one end node is free:

In below an example on this principle:

100
Advanced Concept Training

No the maximal relative deformation is taken as the deformation minus the deformation of the fixed
node.

And uz, relative = 12,057mm – 0,355mm = 11,7mm

And suppose the length of this beam = 1632mm:

 ÊK© å_ = = 1/139
//,¦‡‡
/Zc‡‡

Now this principle is shown at the example of the industrial hall.

Example: Industrial hall

Consider beam B10.

Around the y-y axis (in the z-direction), only the first node has been supported:

Length = 15,075 m

Deformation uz = 84.6 mm (and 0 mm at the beginning of the beam)

uzrelative = 84.6mm – 0.0 mm = 84.6 mm

uzrelative / Length = 84.6mm / 15075mm = 1/178

= 1,12
///¦^
//""
Check =

101
Steel

Fire resistance Check - Not in concept edition

For the fire resistance calculation a Professional or an Expert edition is necessary. The fire resistance
check has been inputted in module esasd.05.xx (esasd.05.01 for the EC-EN).

General

The fire resistance check in Scia Engineer has been inputted following the EN 1993-1-2 – simple
calculation model.

For a selected temperature curve, the temperature in the material after a required period is calculated.
And with this material temperature, the material characteristics are adapted.

The required levels of fire safety depend on factors such as:


- type of occupancy
- height and size of the building
- effectiveness of fire brigade action
- active measures such as vents and sprinklers

Temperature time curves


The user can choose between 4 nominal temperature-time curves in Scia Engineer:

Standard temperature-time curve: Θ® = 20 + 345 ©É݃/" (8 % + 1) [°C]


α = 25 W/m²K

External fire curve: Θ® = 660 (1 − 0,687 K v",c ; − 0,313 K vc,^ ; ) + 20 [°C]


α = 25 W/m²K

Hydrocarbon curve : Θ® = 1080 (1 − 0,325 K v",/Z¦ ; − 0,675 K v,# ; ) + 20 [°C]


α = 50 W/m²K

Θ® = 154 √t + 20 [°C]
§
Smoldering fire:
During 20 minutes followed by the standard ISO 834 curve

102
Advanced Concept Training

This can be set in “Steel -> Beams -> Steel setup”:

This is the temperature of the air during the time.

Steel temperature
Afterwards the steel temperature will be calculated after a certain time with the following formulas. This
steel temperature will be assumed as a uniform temperature in the whole section:

Unprotected steel member:

ઢࣂࢇ,࢚ = ࢑࢙ࢎ
࡭࢓ /y
öࢇ ࣋ࢇ
ࢎሶ࢔ࢋ࢚ ઢ࢚

Protected steel member:

ࣅ ࡭ /y ࣂࢍ,࢚ vࣂࢇ,࢚
ઢࣂࢇ,࢚ = { ∙ (૚Žࣘ/૜) ઢ࢚ − °ࢋ૚} − ૚± ઢࣂࢍ,࢚

 öࢇ ࣋ࢇ

Where:

'ÂÍ Correction factor for shadow effect

For I-sections under nominal fire actions: 'ÂÍ = 0.9 [4‡ /`]‰ /[4‡ /`]

All other cases: 'ÂÍ = [4‡ /`]‰ /[4‡ /`]

4‡ /` Section factor for unprotected steel members [1/m]


In table 4.2 (EN 1993-1-2) some section factors are calculated for unprotected steel
members.

[4‡ /`]‰ the box value for the section factor

ℎሶ w5; the design value on the net heat flux per unit area [W/m²]
This value should be obtained from EN 1991-1-2 with εf = 1,0 and εm = 0,7.
L‫݌ ݌‬
¢ =
LE E p
d ∙ 4 /`

4 /` section factor – see also table 4.3. En 1993-1-2

103
Steel

And the netto heat flux can be calculated according EN 1991-1-2 article 3.1:

ℎሶ w5; = ℎሶ w5;, + ℎሶ w5;,

ℎሶ w5;,
= t bΘ® − Θ‡ d [W/m²]
= heat transfer by convection

ℎሶ w5;,
= Φ ∙ ‡ ∙ 6 ∙  ∙ [(Θ + 273)O − (Θ‡ + 273)O ]
= heat transfer by radiation
[W/m²]

With
αc Coefficient of heat transfer by convection [W/m²K]
Θ® The gas temperature in the vicinity of the fire exposed member [°C]
This temperature may be adopted as nominal temperature-time curves as
given below
Θ‡ The surface temperature of the member [°C]
Stephan Boltzmann constant (=5,67 ∙ 10v^ W/m²K )
4
σ
εm Surface emissivity of fire = 0,7 (EN 1993-1-2)
εf Emissivity of fire = 1
Φ Configuration factor - Φ = 1,0. A lower value may be chosen to take account of
so called position and shadow effects (calculation is given in EN 1991-1-2 –
Annex G).
Θ Is the effective radiation temperature of the fire environment [°C]
In case of fully fire engulfed members, the radiation temperature Θ may be
represented by the gas temperature Θ® around that member.

All the parameters of the previous formulas can be adapted in the Steel setup:

Steel properties
Most of the steel properties will change by a different temperature, so once the steel temperature is
known, the steel properties can be calculated. In below the properties for carbon steel from EN 1993-1-
2, art.3. The properties for stainless steel can be found in EN 1993-1-2: 2005: Annex C.

Reduction factors

Effective yield strength, relative to yield strength at 20°C: ',ఏ = ,ఏ /
Proportional limit, relative to yield strength at 20°C: 'Õ,ఏ = Õ,ఏ /
Slope of linear elastic range, relative to slope at 20°C: '9,ఏ =  Å,ఏ / Å

104
Advanced Concept Training

Thermal elongation

The relative thermal elongation of steel Δ©/© should be determined from the following:

For 20°I ≤ ߠÅ < 750°I:

= 1,2  10v# ߠÅ + 0,4  10v^ ߠÅ −


୼P
P
2,416  10vO
For 750°I ≤ ߠÅ < 860°I:

= 1,1  10v
୼P
P

For 860°I ≤ ߠÅ < 1200°I:

= 2  10v# ߠÅ − 6,2  10vc


୼P
P

Specific heat

For 20°I ≤ ߠÅ < 600°I: LÅ = 425 + 7,73  10v/ ߠÅ −


1,69  10vc ߠÅ + 2,22  10vZ ߠÅc ‫ܬ‬/'݃‫ܭ‬

For 600°I ≤ ߠÅ < 735°I: LÅ = 666 + ‫ܬ‬/'݃‫ܭ‬


/c""
¦c^vఏ¥

/¦^"
For 735°I ≤ ߠÅ < 900°I: LÅ = 545 + ‫ܬ‬/'݃‫ܭ‬
ఏ¥షçô‹

For 900°I ≤ ߠÅ < 1200°I: LÅ = 650 ‫ܬ‬/'݃‫ܭ‬

105
Steel

Thermal conductivity

For 20°I ≤ ߠÅ < 800°I:


Å = 54 − 3,33  10v ߠÅ W/B‫ܭ‬

For 800°I ≤ ߠÅ < 1200°I:


Å = 27,3 W/B‫ܭ‬

Fire resistance properties in Scia Engineer


In Scia Engineer the user can input fire resistance properties on a steel member by “Steel -> Beams ->
Member Check data -> Fire resistance”:

Where:

Time resistance: specifies the required fire resistance

Buckling ratio ky: When this ration is greater than 0,0, the ‘normal’ defined buckling ratio ky
will be overruled by this value.

Buckling ratio kz: When this ration is greater than 0,0, the ‘normal’ defined buckling ratio kz
will be overruled by this value.

Fire exposure: The section may be exposed to fire on all or only three sides.

Covered flange: When a section is exposed to fire on only three sides, the covered flange
can be chosen here.

Protection: Defines whether the section is protected or not.

106
Advanced Concept Training

Insulation LIB: Here the insulation properties can be inputted:

Thickness [mm]: Thickness of the insulation.

k2: This parameters is the adaptation factor for non-uniform temperature


distribution along the beam. This factor κ2 is used for the check on the
design moment resistance Mfi,θ,Rd.

Temperature-Time curve: The user can choose for another temperature time as defined in the steel
setup or can even input his own temperature-time curve here.

Resistance domain
In Scia Engineer 2 calculation methods are implemented:
- Resistance domain
- Temperature time domain.

The choice between those analysis types can be made in “Steel -> Beams -> Steel setup”:

In this chapter the fire resistance check following the Resistance domain will be explained and in the
next paragraph the Temperature/time domain.

107
Steel

Principle
The user will choose the used temperature time curve and will input a required fire resistance time.
After this time the temperature of the gas and afterwards of the steel will be calculated.

With this steel temperature the reduced properties will be calculated and a fire resistance check
according to EN 1993-1-2, art. 4 will be executed with those adapted steel properties. This check will
result in a unity check, which is the fire resistance check for the resistance domain.

Example in Scia Engineer


This principle is explained with an example in Scia Engineer.

Example: Industrial hall.esa

Consider member B28:

Following fire resistance properties are inputted on this column:

So the fire resistance will be checked after 30 minutes (= 1800 seconds) with a fibre board protection
and the buckling factors are taken equal as the buckling factors of the normal steel code check.

The fire resistance check is executed via “Steel -> Beams -> ULS check -> Check - fire resistance” for
this column and for combination “CO3 - Fire”, resulting in a unity check of 0.61:

108
Advanced Concept Training

When looking at the default output, this calculation is given by Scia Engineer:

First the partial safety factors are given:

Afterwards the material properties (not adapted by the temperature) are given:

And the fire resistance properties as inputted in Scia Engineer. Here is also indicated that the fire
resistance check has been executed after 30 minutes of fire.

And now a graph is shown with the gas temperature (in this example follow the ISA 834 curve), the
steel temperature calculated with a protection and the reduction of the yield strength.

109
Steel

And after all the unity check is shown with the reduced properties:

110
Advanced Concept Training

111
Steel

So in this example the Compression and bending check will result in a unity check op 0.61.

Temperature/time domain
In Scia Engineer 2 calculation methods are implemented:
- Resistance domain
- Temperature time domain.

The choice between those two analysis types can be made in “Steel -> Beams -> Steel setup”. In this
chapter the fire resistance check following the Temperature/time domain will be explained.

Principle
The user will choose the used temperature time curve and will input a required fire resistance time.
After this time the temperature of the gas and afterwards of the steel will be calculated.

Now also the critical steel temperature will be calculated. And the Fire resistance check according to
the temperature/time domain will be the ratio of the real steel temperature after a chosen time and the
critical steel temperature.

This critical steel temperature θa,cr can be calculated with a simple formula of in an iterative way and
this choice can be made in “Steel -> Beams -> Steel setup”:

112
Advanced Concept Training

When the Iterative process is not checked, the critical steel temperature is calculated with a simple
formula:

ߠÅ, = 39,19 ©G Ò − 1Ó + 482


/
",\Z¦O ఓAô,ßôô

Where µ0 is the degree of utilization. This means the unity check following EN 1993-1-2 at time = 0sec,
so without augmentation of the temperature.

NOTE
This simple formula is only valid if now stability phenomena or deformation criteria have to be taken
into account.

When the Iterative process option is checked, the critical steel temperature will be calculated with an
iterative process. So first an estimation of this critical temperature will be chosen and the unity check
following EN 1993-1-2 will be executed, if this check is lower than one, a higher critical temperature is
chosen or when this check is higher than one, a lower temperature is chosen. Now this unity check is
recalculated just until the moment this unity check gives a result for this critical steel temperature
between 0.99 and 1.01.

This is a more accurate procedure to calculate the critical temperature and this method is also valid if
stability phenomena or deformation criteria have to be taken into account.

Example in Scia Engineer


This principle is explained with an example in Scia Engineer.

Example: Industrial hall.esa

Consider member B28:

And change in “Steel -> Beams -> Steel setup” the analysis type to “Temperature/time domain” without
activating the option “iterative process”.

The fire resistance check is executed via “Steel -> Beams -> ULS check -> Check - fire resistance” for
this column and for combination “CO3 - Fire”, resulting in a unity check of 0.90:

113
Steel

When looking at the default output, this calculation is given by Scia Engineer:

So the check following EN 1993-1-2 will be given at t = 0min, so 20°C, without any reduction of the
steel properties. This check will result in a low unity check, which is the degree of utilisation µ0:

In this example again the Compression and bending check will result in the highest unity check = 0.26
= µ0.

This value is used in the simple formula for the critical steel temperature:

ߠÅ, = 39,19 ©G Ò − 1Ó + 482


/
",\Z¦O ఓAô,ßôô

ߠÅ, = 39,19 ©G ú − 1û + 482 = 683 °I


/
",\Z¦O (".Z)ô.ßôô

And the steel temperature after 30 minutes is 614.72°C.

= 0.90
Z/O.¦°£
Z^c.Zc°£
The unity check is:

All those values are also given in the overview table in the preview of the fire resistance check:

114
Advanced Concept Training

And also the Fire resistance time is given in this table:


This member van resist fire for 37.68 minutes.

NOTE
As stated before, this simple calculation method can only be used if no stability phenomena have to be
taken into account. In this example the stability causes the highest unity check and is thus taken into
account, so this method is not correct and should not be used!

Now go to “Steel -> Beams -> Steel setup” and check the option “iterative process”.

The fire resistance check is executed via “Steel -> Beams -> ULS check -> Check - fire resistance” for
this column and for combination “CO3 - Fire”, resulting in a unity check of 0.99:

When looking at the default output, this calculation is given by Scia Engineer:

So the check following EN 1993-1-2 will be given at θa,cr =622.7°C. This temperature is calculated
iterative resulting in a unity check following the EN 1993-1-2 equal to 1:

115
Steel

And indeed the highest unity check will be equal to 1.00 for this critical temperature.

= 0.99
Z/O.¦°£
Z.Z^°£
So for this case the unity check is:

All those values are also given in the overview table in the preview of the fire resistance check:

And also the Fire resistance time is given in this table:


This member van resist fire for 30.78 minutes.

NOTE
As stated before, this simple calculation method can only be used if no stability phenomena have to be
taken into account. In this example the critical steel temperature was calculated using the iterative
process, so this method is correct and can be used as an alternative for the resistance domain.

116
Advanced Concept Training

References and literature

[1] Eurocode 3
Design of steel structures
Part 1 - 1 : General rules and rules for buildings
ENV 1993-1-1:1992, 1992

[2] Essentials of Eurocode 3


Design Manual for Steel Structures in Building
ECCS - N° 65, 1991

[3] R. Maquoi
ELEMENTS DE CONSTRUCTIONS METALLIQUE
Ulg , Faculté des Sciences Appliquées, 1988

[4] ENV 1993-1-3:1996


Eurocode 3 : Design of steel structures
Part 1-3 : General rules
Supplementary rules for cold formed thin gauge members and sheeting
CEN 1996

[5] Eurocode 3
Design of steel structures
Part 1 - 1/ A1 : General rules and rules for buildings
ENV 1993-1-1:1992/A1, 1994

[6] Eurocode 3
Design of steel structures
Part 1 - 2 : General rules - Structural fire design
ENV 1993-1-2:2005, 2005

[7] Model Code on Fire Engineering


ECCS - N° 111
May 2001

[8] Eurocode 1
Basis of design and actions on structures
Part 2-2 : Actions on structures - Actions on structures exposed to fire
ENV 1991-2-2:1995

[9] DIN 18800 Teil 1


Stahlbauten
Bemessung und Konstruktion
DK 693.814.014.2, November 1990

[10] DIN 18800 Teil 2


Stahlbauten
Stabilitätsfälle, Knicken von Stäben und Stabwerken
DK 693.814.074.5, November 1990

[11] DIN 18800 Teil 3


Stahlbauten
Stabilitätsfälle, Plattenbeulen
DK 693.814.073.1, November 1990

[12] G. Hünersen, E. Fritzsche


Stahlbau in Beispielen

117
Steel

Berechnungspraxis nach DIN 18 800 Teil 1 bis Teil 3 (11.90)


Werner-Verlag, Düsseldorf 1991

[13] E. Kahlmeyer
Stahlbau nach DIN 18 800 (11.90)
Werner-Verlag, Düsseldorf

[14] Beuth-Kommentare
Stahlbauten
Erläuterungen zu DIN 18 800 Teil 1 bis Teil 4, 1.Auflage
Beuth Verlag, Berlin-Köln 1993

[15] Stahlbau Kalender 1999


DSTV
Ernst & Sohn, 1999

[16] Roik, Carl, Lindner


Biegetorsionsprobleme gerader dünnwandiger Stäbe
Verlag von Wilhelm Ernst & Sohn
1972

[17] DASt-Richtlinie 016


Bemessung und konstruktive Gestaltung von Tragwerken aus dünnwandigen kaltgeformted
Bauteilen
Stahlbau-Verlagsgesellschaft - 1992

[18] H. Rubin,
Interaktionsbeziehungen für doppeltsymmetrische I- und Kasten-Querschnitte bei zweiachsiger
Biegung und Normalkraft
Der Stahlbau 5/1978, 6/1978

[19] Stahl im Hochbau


14. Auflage, Band I / Teil 2
1986, Verlag Stahleisen mbH, Düsseldorf

[20] Staalconstructies TGB 1990


Basiseisen en basisrekenregels voor overwegend statisch belaste constructies
NEN 6770, december 1991

[21] Staalconstructies TGB 1990


Stabiliteit
NEN 6771, december 1991

[22] Staalconstructies TGB 1990


Stabiliteit
NEN 6771, januari 2000

[23] NEN 6072


Rekenkundige bepaling van de brandwerendheid van bouwdelen
Staalconstructies
December 1991

[24] NEN 6072/A2 - Wijzigingsblad


Rekenkundige bepaling van de brandwerendheid van bouwdelen
Staalconstructies
December 2001

[25] NEN 6702


Belastingen en vervormingen TGB 1990
December 1991

[26] Handleiding moduul STACO VGI

118
Advanced Concept Training

Staalbouwkundig Genootschap
Staalcentrum Nederland
5684/82

[27] Newmark N.M. A simple approximate formula for effective end-fixity of columns
J.Aero.Sc. Vol.16 Feb.1949 pp.116

[28] Stabiliteit voor de staalconstructeur


uitgave Staalbouwkundig Genootschap

[29] Rapportnr. BI-87-20/63.4.3360


Controleregels voor lijnvormige constructie-elementen
IBBC Maart 1987

[30] Y. Galéa
Flambement des poteaux à inertie variable
Construction Métallique 1-1981

[31] I. Vayas,
Interaktion der plastischen Grenzschnittgrössen doppelsymmetrischer I-Querschnitte
Stahlbau 69 (2000), Heft 9

[32] Scia Engineer Steel Code Check


Theoretical Background
04/2011

[33] Eurocode 3
Design of steel structures
Part 1 - 1 : General rules and rules for buildings
EN 1993-1-1:2005

[34] Rules for Member Stability in EN 1993-1-1


Background documentation and design guidelines
ECCS - N° 119
2006

[35] Lateral-Torsional buckling of steel beams:


A general expression for the moment gradient factor.
A. López, D. J. Yong, M. A. Serna
Stability and Ductility of Steel Structures, 2006.

119

You might also like