(Eng) Advanced Concept Training Steel 2011.0 - Nemetschek Scia
(Eng) Advanced Concept Training Steel 2011.0 - Nemetschek Scia
(Eng) Advanced Concept Training Steel 2011.0 - Nemetschek Scia
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
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:
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.
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
sheet welded Iw
L Angle Section
from library
U Channel Section
3
Steel
from library
T T Section
from library
sheet welded Tw
Z Z Section
from library
RS Rectangular Sections
from library
O Solid Tube
from library
from library
pairs 2CFUo
pairs 2CFUc
4
Advanced Concept Training
pairs 2CFCo
pairs 2CFCc
pairs 2CFLT
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]).
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
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:
8
Advanced Concept Training
Section classification
The classification is executed following EN 1993-1-1, art. 5.5.
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.
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
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 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)
$
/% 59.57
=
= = 1.05
28.4'( 28.4 ∙ 1.00 ∙ √4.0
1.05 − 0.22
= = 0,75
(1.05)
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.
The user can choose to perform an elastic check. This corresponds with a class 3 check.
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
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.
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
Consider bracing B234 (for load case 3DWind1). In this load case the bracing in under tension:
17
Steel
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
The classification on this position is Class 1 (column under compression and bending).
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,<: = CP,<: =
?@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
The classification on this position is Class 1 (column under compression and bending).
19
Steel
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.
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<:
21
Steel
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.
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:
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:
On the two beams a torsional moment (8 kNm) has been inputted and a small line load (-0.1 kN/m).
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.
C,9: C],9:
r s
o q +o q ≤1
Cp,,<: Cp,],<:
- 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:
If yz{ ≥ }, ~ y,,{
(1 − )
2 `9:
Where
= − 1
`P,<:
And
=a ∙b6> /√cd
`P,<: = ?@A
2 `9:
= − 1
`P,
,<:
Class 3 cross-sections
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
In absence of shear force, for Class 4 cross-sections the maximum longitudinal stress should satisfy
the criterion:
6>
,9: ≤ ?
@A
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.
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
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:
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:
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).
89:
≤1
8,<:
Where
∙=∙6>
8,<: =
?@
- For class 1, 2 or 3 cross sections
∙=MNN ∙6>
8,<: =
?@
- For class 4 cross sections
=
/
²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)
²²
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 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
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”:
33
Steel
' = (
( # # O)( O O /)
// # O)( # // O)
' = +4
²
4/ + ²/
=
(/+ ) + 8/
I
=
¡
C
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
²∙/""""p/²∙/Z¦""""" §
8 = = = 37522008 = 3752.2 '8
²9l
²² (.c¦) (O""")²
o
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
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
For a VARH element, Scia Engineer will use another calculation for the buckling length.
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
t² ¡
8 =
'²²
2. We suppose that
α 1
>
ξ -1 2
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
Example:
VARH.esa
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.
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:
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
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
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:
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:
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
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=
46
Advanced Concept Training
1 ² I
= (¯¡; + )
,
4® H" ©
1 ² I
= (¯¡; + )
,
4® H" ©
1
− b
= ob + ,
d + ,
d − 4u , ,
q
,
2u , ,
²
=
©
,
° ±
H
"
u =1−° ±
H"
47
Steel
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
λ
= ª
C
Mcr Critical bending moment (see next paragraph for this calculation)
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:
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:
/
¶´µ
And
,¾: =
f may be defined in the National Annex.
50
Advanced Concept Training
' = 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:
The default method for the calculation of the C-factors has been used. This is also indicated in the
preview of the check results:
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:
Where
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:
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.
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”:
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.
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.
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Å
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.
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:
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:
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
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”:
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
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 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
60
Advanced Concept Training
Example: Diaphragm.esa
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
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”.
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:
63
Steel
Unless second order analysis is carried out, members which are subjected to combined bending and
axial compression should satisfy:
Ø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.
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
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
In Scia Engineer:
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
67
Steel
Stiffners
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
69
Steel
This check can be executed following the EC-ENV and EC-EN 11993-1-1
art.6.4
(1) 2I
(2) 2Uo
(3) 2Uc
Consider member B1
On this beam, only a compression force of 500kN is applied
Properties of the batten:
70
Advanced Concept Training
H0 = = = 110,833BB
P /.#c^ /"ß §
∙=Í ∙ZZ"
Mu à:
= 0 if Lambda > 150
= 2−
Å:Å
¦#
if 75 < Lambda < 150
71
Steel
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
Å (//#")
With:
AAAA¤
²∙f g∙(/.#c^ /"ß § )
8 = = = 15064 '8
²∙9∙l566 TT
² (OZ"")²
1− − = 0,9466 > 0
#""p #""p
/#"ZOp O¦"" p
= = 3,32'8
∙Y9: ∙O^Z"p
P OZ""
VEd
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:
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
Nch,Ed = NChord = 170,78 kN (see check of Check of chord as beam in field between battens)
Buckling length = a
= = 32,411
ÇîPw® 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
(ÇîPw® 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:
Check of batten:
t thickness of batten
b width of batten
= = = 9.09'8
h9:∙Å c,cp∙//#"
Í"∙ /"∙
T
73
Steel
= = 16,9CDE
Yò∙/ \#O"""p∙"/
l Z.//∙/"è§
Sigma
= = 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
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.
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
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:
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.
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”.
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.
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”.
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.
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
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:
84
Advanced Concept Training
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
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
Increase sway
effects with: No Members Yes
1 with e0
1
1−
α cr
5.2.2(3)c 5.2.2(7)b
Section Check
86
Advanced Concept Training
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.
1
õ= ∙t ∙t
200 Í
tÍ = ≤ tÍ ≤ 1,0
√Í c
but
1
t = ª0,5 °1 + ±
B
The global imperfection can be inputted in Scia Engineer in the non linear combinations as explained
by the example below.
87
Steel
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:
These inclination functions are entered through “Main -> Library -> Structure, Analysis -> Initial
deformations”:
88
Advanced Concept Training
89
Steel
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:
̅ > 0,5
= ∙6>
•
pij
90
Advanced Concept Training
̅ =
= ∙6>
p
And
> 0,5
= ∙6> = ∙6>
Thus: p pij
> 0,25
/ /
p pij
Bow imperfection e0
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.
In Scia Engineer the user can choose between three options for the bow imperfection in the non
linear combinations:
91
Steel
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
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.
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.
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
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
This can be verified with Euler’s formula using the member length as the buckling length:
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.
Calculation of e0
= = 0,945
/
v",db
",#ú/rb d û f",#ú/rb d ûg vb
v",db 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Ï
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
ø" ,Å
1
-03
at mid section = 1,44E
mm²
øw; = K" ∙ 9∙l ∙ ø
p
> ∙ఎ,T¥ೣ
"
ø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
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.
Vertical deformation
The maximum displacement in the Z direction is 59.6 mm
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.
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:
99
Steel
Deformation in the beginning of the beam = 5,967mm and at the end = 11,863mm.
0,979B
å_Pî5 Pw5;",\¦\ = 5,967BB + ∗ (11,863BB − 5,967BB) = 9,504BB
1,632B
ÊK© å_ = = 1/2236
",¦c
/Zc
100
Advanced Concept Training
No the maximal relative deformation is taken as the deformation minus the deformation of the fixed
node.
ÊK© å_ = = 1/139
//,¦
/Zc
Around the y-y axis (in the z-direction), only the first node has been supported:
Length = 15,075 m
= 1,12
///¦^
//""
Check =
101
Steel
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.
Θ® = 154 √t + 20 [°C]
§
Smoldering fire:
During 20 minutes followed by the standard ISO 834 curve
102
Advanced Concept Training
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:
ઢࣂࢇ,࢚ = ࢙ࢎ
/y
öࢇ ࣋ࢇ
ࢎሶࢋ࢚ ઢ࢚
ࣅ /y ࣂࢍ,࢚ vࣂࢇ,࢚
ઢࣂࢇ,࢚ = { ∙ (ࣘ/) ઢ࢚ − °ࢋ} − ± ઢࣂࢍ,࢚
ࣘ
öࢇ ࣋ࢇ
Where:
For I-sections under nominal fire actions: 'ÂÍ = 0.9 [4 /`] /[4 /`]
ℎሶ 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 /`
103
Steel
And the netto heat flux can be calculated according EN 1991-1-2 article 3.1:
ℎሶ 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:
= 1,1 10v
P
P
Specific heat
/¦^"
For 735°I ≤ ߠÅ < 900°I: LÅ = 545 + ܬ/'݃ܭ
ఏ¥షçô
105
Steel
Thermal conductivity
Where:
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.
106
Advanced Concept Training
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.
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:
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:
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.
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:
= 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
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:
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
[1] Eurocode 3
Design of steel structures
Part 1 - 1 : General rules and rules for buildings
ENV 1993-1-1:1992, 1992
[3] R. Maquoi
ELEMENTS DE CONSTRUCTIONS METALLIQUE
Ulg , Faculté des Sciences Appliquées, 1988
[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
[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
117
Steel
[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
[18] H. Rubin,
Interaktionsbeziehungen für doppeltsymmetrische I- und Kasten-Querschnitte bei zweiachsiger
Biegung und Normalkraft
Der Stahlbau 5/1978, 6/1978
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
[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
[33] Eurocode 3
Design of steel structures
Part 1 - 1 : General rules and rules for buildings
EN 1993-1-1:2005
119