So Fistik
So Fistik
So Fistik
SOFiSTiK | 2018
VERiFiCATiON MANUAL
SOFiSTiK AG
This manual is protected by copyright laws. No part of it may be translated, copied or reproduced, in any form or by any means,
without written permission from SOFiSTiK AG. SOFiSTiK reserves the right to modify or to release new editions of this manual.
The manual and the program have been thoroughly checked for errors. However, SOFiSTiK does not claim that either one is
completely error free. Errors and omissions are corrected as soon as they are detected.
The user of the program is solely responsible for the applications. We strongly encourage the user to test the correctness of all
calculations at least by random sampling.
Front Cover
Project: New SOFiSTiK Office, Nuremberg | Contractor: WOLFF & MLLER, Stuttgart | Architecture: WABE-PLAN ARCHITEKTUR, Stuttgart |
Structural Engineer: Boll und Partner. Beratende Ingenieure VBI, Stuttgart | MEP: GM Planen + Beraten, Griesheim | Lead Architect: Gerhard P.
Wirth gpwirtharchitekten, Nuremberg | Vizualisation: Armin Dariz, BiMOTiON GmbH
Contents
Contents
Introduction 3
1 About this Manual 3
1.1 Layout and Organization of a Benchmark .............................................. 3
1.2 Finding the Benchmark of interest ....................................................... 3
1.3 Symbols .................................................................................... 5
2 Index by Categories 7
2.1 Design Code Benchmarks................................................................ 7
22 DCE-EN20: Design of a Steel I-section for Bending, Shear and Axial Force 141
22.1 Problem Description....................................................................... 141
22.2 Reference Solution ........................................................................ 142
22.3 Model and Results ........................................................................ 142
22.4 Design Process ............................................................................ 143
22.5 Conclusion ................................................................................. 146
22.6 Literature ................................................................................... 147
27 DCE-EN25: Shear between web and flanges of Hollow CS acc. DIN EN 1992-2 193
27.1 Problem Description....................................................................... 193
27.2 Reference Solution ........................................................................ 193
27.3 Model and Results ........................................................................ 194
27.4 Design Process ............................................................................ 197
27.5 Conclusion ................................................................................. 201
27.6 Literature ................................................................................... 202
28 DCE-EN27: Design of Quad Elements - Layer Design and Baumann Method 203
28.1 Problem Description....................................................................... 203
28.2 Reference Solution ........................................................................ 203
28.3 Model and Results ........................................................................ 204
28.4 Design Process ............................................................................ 208
28.5 Conclusion ................................................................................. 215
28.6 Literature ................................................................................... 215
29 DCE-EN28: Design of a T-section for bending and shear using Design Elements 217
29.1 Problem Description....................................................................... 217
29.2 Reference Solution ........................................................................ 217
29.3 Model and Results ........................................................................ 217
29.4 Design Process ............................................................................ 220
29.5 Conclusion ................................................................................. 222
29.6 Literature ................................................................................... 222
33 DCE-EN32: Crack width calculation of reinforced slab acc. DIN EN 1992-1-1 263
33.1 Problem Description....................................................................... 263
33.2 Reference Solution ........................................................................ 263
33.3 Model and Results ........................................................................ 264
33.4 Design Process ............................................................................ 267
33.4.1 Material .................................................................................. 267
33.4.2 Without direct calculation .............................................................. 267
33.4.3 With direct calculation .................................................................. 268
33.5 Conclusion ................................................................................. 271
33.6 Literature ................................................................................... 271
36 DCE-MC1: Creep and Shrinkage Calculation using the Model Code 2010 297
36.1 Problem Description....................................................................... 297
36.2 Reference Solution ........................................................................ 297
36.3 Model and Results ........................................................................ 297
36.4 Design Process ............................................................................ 299
36.5 Conclusion ................................................................................. 303
36.6 Literature ................................................................................... 303
37 DCE-MC2: Creep and Shrinkage Calculation using the Model Code 1990 305
37.1 Problem Description....................................................................... 305
37.2 Reference Solution ........................................................................ 305
37.3 Model and Results ........................................................................ 305
37.4 Design Process ............................................................................ 307
37.5 Conclusion ................................................................................. 310
37.6 Literature ................................................................................... 310
Introduction
The primary objective of this manual is to verify the capabilities of SOFiSTiK by means of nontrivial
problems which are bound to reference solutions.
To this end, this manual contains a compilation of a number of selected computational benchmarks, each
benchmark focusing on a specific (mechanical/design) topic. The obtained results from the SOFiSTiK
analysis are contrasted with corresponding reference solutions1 .
The tasks covered by SOFiSTiK, address a broad scope of engineering applications and it is therefore
not possible to validate all specific features with known reference solutions in terms of this Verification
Manual. An attempt has been made though, to include most significant features of the software with
respect to common problems of general static and dynamic analysis of structures.
• Problem Description
• Reference Solution
• Conclusion
• Literature
First, the problem description is given, where the target of the benchmark is stated, followed by the
reference solution, where usually a closed-form analytical solution is presented, when available. The
next section is the description of the model, where its properties, the loading configuration, the analysis
method and assumptions, further information on the finite element model, are presented in detail. Finally,
the results are discussed and evaluated with respect to the reference solution and a final conclusion for
the response of the software to the specific problem is drawn. Last but not least, the textbooks and
references used for the verification examples are listed, which are usually well known and come from
widely acclaimed engineering literature sources.
Overview
1 Where available, analytical solutions serve as reference. Where this is not feasible, numerical or empirical solutions are
referred to. In any case, the origin of the reference solution is explicitly stated.
Procedure(s): LSTP
Topic(s): SOIL
Module(s): TALPA
Input file(s): passive earth pressure.dat
As it can be seen, the available categories are the element type, the analysis type, the procedure, the
topics and the modules. For each category that is provided in the description table, a hyperlink is created,
linking each example to the global categories tables. In this manner, the user has a direct overview of
the attributes involved in each problem, and at the same time is able to browse by category through the
Verification Manual focusing only on the one of his interest. Table 1.1 provides an overview of all the
categories options that are available.
Categories Options
Continuum 3D
Continuum 2D (plane strain)
Continuum axisymmetric
Shell
FE beam 3D
Nonlinear FE beam 3D (AQB)
Element Type
Fiber beam 2D
Fiber beam 3D
Truss element
Cable element
Spring element
Damping element
Couplings
Geometrically nonlinear
Physically nonlinear
Static
Potential problem
Buckling analysis
Load stepping
Phi-C reduction
Soil related
Topic
Categories Options
Seismic
Fire design
AQB
AQUA
ASE
Module BDK
BEMESS
CSM
DYNA
SOFiLOAD
SOFiMSHC
STAR2
TALPA
TENDON
1.3 Symbols
For the purpose of this manual the following symbols and abbreviations apply.
SOF. SOFiSTiK
Ref. reference
Tol. tolerance
cs cross-section
sect. section
temp. temperature
homog. homogeneous
Be. benchmark
con. construction
SDOF single degree of freedom
2 Index by Categories
Subsequent tables show all Benchmarks included in this Verification Manual, indexed by category.
Eurocode for Design of Concrete DIN EN 1992-1-1 DCE-EN1, DCE-EN2, DCE-EN3, DCE-EN4,
Structures - General rules and DCE-EN5, DCE-EN6, DCE-EN7, DCE-EN8,
rules for building incl. German DCE-EN9, DCE-EN10, DCE-EN11, DCE-
National Annex EN12, DCE-EN17, DCE-EN18, DCE-EN19,
DCE-EN21, DCE-EN22, DCE-EN23, DCE-
EN25
, DCE-EN27
, DCE-EN28, DCE-EN31, DCE-EN32, DCE-
EN33
Eurocode for Design of Concrete DIN EN 1992-2 DCE-EN25
Structures - Concrete bridges
incl. German National Annex
Eurocode for Design of Steel DIN EN 1993-1-1 DCE-EN13, DCE-EN14, DCE-EN15, DCE-
Structures - General rules and EN20, DCE-EN24
rules for buildings incl. German , DCE-EN29, DCE-EN30, DCE-EN34
National Annex
CEB-FIP Model Code for MC 1990 DCE-MC2
Concrete Structures 1990
CEB-FIP Model Code for MC 2010 DCE-MC1
Concrete Structures 2010
Swiss Standard for Concrete SIA 262 DCE-SIA1
Structures
Design of Cross Sections and of AQB DCE-EN1, DCE-EN2, DCE-EN3, DCE-EN4, DCE-
Prestressed Concrete and EN5, DCE-EN6, DCE-EN7, DCE-EN8, DCE-EN9,
Composite Cross Sections DCE-EN10, DCE-EN11, DCE-EN12, DCE-EN13,
DCE-EN14, DCE-EN17, DCE-EN18, DCE-EN19,
DCE-EN20, DCE-EN22, DCE-EN23, DCE-EN25
, DCE-EN28, DCE-EN30, DCE-MC1, DCE-MC2
3.1 Objectives
3.1.1 About SOFiSTiK
SOFiSTiK finite element software has been continuously developed since 1981. It is currently used
by more than 10000 customers worldwide. SOFiSTiK is a multipurpose tool with extensive capabilities
which fall into a wide spectrum of engineering analyses such as static and dynamic structural analysis,
modal and buckling eigenvalue problems, nonlinearities and higher order effects, geotechnics and tunnel
analysis, heat transfer and fire analysis, as well as numerous types of other applications.
3.2 Organisation
3.2.1 Software Release Schedule
The SOFiSTiK software release schedule is characterized by a two-year major release cycle. The
first customer shipment (FCS) of a SOFiSTiK major release is preceded by an extensive BETA testing
period. In this phase - after having passed all internal test procedures (Section 3.2.4: Software Release
Procedure) - the new product is adopted for authentic engineering projects both by SOFiSTiK and by
selected customers. For a two-year transition period, subsequent major releases are fully supported in
parallel, as shown in Fig. 3.1.
FCS
Planned: Summer 2019
BETA SOFiSTiK 2020
The major release cycle is supplemented by a two-month service pack cycle. Service packs are quality
assured, which means they have passed both the continuous testing procedures and the functional tests
(Section 3.2.2: SQA Modules - Classification). They are available for download via the SOFiSTiK update
tool SONAR.
Software updates for the current version (service packs) include bug-fixes and minor new features only;
major new developments with increased potential regarding side- effects are reserved for major releases
with an obligatory extensive testing period.
Preventive provisions essentially concern the organization of the development process. They aim at
minimizing human errors by a high degree of automatism and by avoiding error-prone stress situations.
These provisions comprise:
• Strict phase differentiation: Prior to any software release (also for service packs), the development
phase is followed up by a consolidation phase . This phase is characterized by extensive functional
testing. No new features are implemented, only test feedback is incorporated. For major releases,
an additional BETA test phase is scheduled.
Analytic provisions provide for the actual testing of the software products. Continuous Testing directly
accompanies the development process: Automated and modular regression tests assure feedback at a
very early stage of the development (Section 3.3.3: Continuous Testing). Functional Testing is carried
out in particular during the consolidation phases. These tests essentially involve manual testing; they
focus on comprehensive workflow tests and product oriented semantic tests.
SQA
Development Process
Continuous Testing
Phase differentiation, build- and Functional Testing
Automated, modular and
publishing mechanism
continuous regression testing Focus: workflow tests, product
Organizational Provision oriented semantic tests
Instruments (analytic)
(preventive)
3.2.3 Responsibilities
The consistent implementation of quality assurance procedures is responsibly coordinated by the man-
aging board executive for products.
Product
Development
Management
Continuous Integration
Functional testing
Continuous Testing
Integrating customer
Implementing feedback
feedback
Corporate Activity
Adaption of processes
Definition of objectives
Coordinated by managing board
Release
Requirements
3.3 Instruments
3.3.1 CRM System
Each request from our customers is traced by means of a Customer Relation Management (CRM)
System assuring that no case will be lost. Detailed feedback to the customer is provided via this system.
Possible bug fixes or enhancements of the software are documented with version number and date
in corresponding log files. These log files are published via RSS-feed to our customers. In this way,
announcement of available software updates (service-pack or hotfix) is featured proactively. Moreover,
information is provided independent of and prior to the actual software update procedure.
Further sources of information are the electronic newsletter/ newsfeeds and the internet forum
(www.sofistik.de / www.sofistik.com).
Automated Automated
Continuous Testing Continuous Integration
procedure procedure
Development/ PM
Assessing feedback
Committing modifications
Continuous integration denotes the automated process, assuring that all executed and committed mod-
ifications of the program’s code basis are directly integrated via rebuild into the internal testing environ-
ment.
Upon completion of the integration, the continuous testing procedure is triggered automatically. This
procedure executes a standardized testing scenario using the newly updated software. Test results are
prepared in form of compact test protocols allowing for quick assessment.
The executed tests are so-called regression tests. Regression tests examine by means of associated
reference solutions wether the conducted modifications of the code basis cause undesired performance
in other already tested parts of the program.
Together, continuous integration and continuous testing form the basis for a quality control that directly
accompanies the development process. This way, possibly required corrections can be initiated promptly.
SOFiSTiK has successfully implemented this procedure. Currently, the continuous test database com-
prises more than 3000 tests.
It is the credo of SOFiSTiK that a high-quality product can only be created and maintained by highly qual-
ified personnel. Continuing education of the staff members is required by SOFiSTiK and it is supported
by an education program which involves both in- house trainings and provisions of external trainings on
a regular basis.
3.5 Disclaimer
Despite all efforts to achieve the highest possible degree of reliability, SOFiSTiK cannot assure that the
provided software is bug-free or that it will solve a particular problem in a way which is implied with the
opinion of the user in all details. Engineering skill is required when assessing the software results.
Overview
h d mEd
s1 s1
• Design stress-strain curves for concrete and reinforcement (Section 3.1.7, 3.2.7)
Ac Fc
d z
As Fs
εs
In singly reinforced beams and slabs, the conditions in the cross-section at the ultimate limit state, are
assumed to be as shown in Fig. 4.2. The design stress-strain diagram for reinforcing steel considered
in this example, consists of an inclined top branch, as presented in Fig. 4.3 and as defined in DIN EN
1992-1-1:2004 [1] (Section 3.2.7).
A
σ
ƒtk,c = 525N/ mm2
ƒyk
ƒyd = ƒyk / γs
B A Idealised
B Design
εd = 25 0 /00
Figure 4.3: Idealised and Design Stress-Strain Diagram for Reinforcing Steel
B 500A d = 17.0 cm
b = 1.0 m
Material:
Design Load:
MEds 25 · 10−3
μEds = = = 0.061
b · d2 · ƒcd 1.0 · 0.172 · 14.17
Tab. 9.1 [3]: General Chart for up to
ε = 25 · 10−3; ζ = 0.97 → σs1d = 456.52 MP C50/ 60 - Cross-section without com-
pression reinforcement
1 MEds
s1 = · + NEd = 3.328 cm2 / m
σs1d ζ·d
1
s1 = · (ω · b · d · ƒcd + NEd ) = 3.334 cm2 / m
σsd
MEds NEd
s1 = ks · + · κs = 3.333 cm2 / m
d σs1d
1 The tools used in the design process are based on steel stress-strain diagrams, as
defined in [1] 3.2.7:(2), Fig. 3.8, which can be seen in Fig. 4.3.
2 The sections mentioned in the margins refer to DIN EN 1992-1-1:2004 (German Na-
4.5 Conclusion
This example shows the calculation of the required reinforcement for a slab section under bending.
Various different reference solutions are employed in order to compare the SOFiSTiK results to. It has
been shown that the results are reproduced with excellent accuracy.
4.6 Literature
[1] DIN EN 1992-1-1/NA: Eurocode 2: Design of concrete structures, Part 1-1/NA: General rules and
rules for buildings - German version EN 1992-1-1:2005 (D), Nationaler Anhang Deutschland -
Stand Februar 2010. CEN. 2010.
[2] F. Fingerloos, J. Hegger, and K. Zilch. DIN EN 1992-1-1 Bemessung und Konstruktion von
Stahlbeton- und Spannbetontragwerken - Teil 1-1: Allgemeine Bemessungsregeln und Regeln
für den Hochbau. BVPI, DBV, ISB, VBI. Ernst & Sohn, Beuth, 2012.
[3] K. Holschemacher, T. Müller, and F. Lobisch. Bemessungshilfsmittel für Betonbauteile nach Eu-
rocode 2 Bauingenieure. 3rd. Ernst & Sohn, 2012.
[4] Beispiele zur Bemessung nach Eurocode 2 - Band 1: Hochbau. Ernst & Sohn. Deutschen Beton-
und Bautechnik-Verein E.V. 2011.
[5] R. S. Narayanan and A. W. Beeby. Designers’ Guide to EN 1992-1-1 and EN 1992-1-2 - Eurocode
2: Design of Concrete Structures. Thomas Telford, 2005.
Overview
As2
h d d2 MEd
As1
• Design stress-strain curves for concrete and reinforcement (Section 3.1.7, 3.2.7)
εc σcd
d2 Fs2d
σs2d
As2 εs2
Fc
d zs2
z
zs1
As1
d1 εs1 σs1d Fs1d
Figure 5.2: Stress and Strain Distributions in the Design of Doubly Reinforced Cross-sections
In doubly reinforced rectangular beams, the conditions in the cross-section at the ultimate limit state, are
assumed to be as shown in Fig. 5.2. The design stress-strain diagram for reinforcing steel considered
in this example, consists of an inclined top branch, as presented in Fig. 5.3 and as defined in DIN EN
1992-1-1:2004 [1] (Section 3.2.7).
A
σ ƒtk,c = 525N/ mm2
ƒyk
ƒyd = ƒyk / γs
B A Idealised
B Design
εd = 25 0 /00
Figure 5.3: Idealised and Design Stress-Strain Diagram for Reinforcing Steel
B 500A d = 35.0 cm
d2 = 5.0 cm
b = 25 cm
Material:
NEd = 0
MEds = MEd − NEd · zs1 = 135 kNm
from design chart for μEds,m = 0.296 and d2 / d = 0.143 : Tab. 9.1 [3]: General Chart for up to
εs1 = 4.30 · 10−3 ; εs2 = −2.35 · 10−3 ; ζ = z/ d = 0.813 C50/ 60 - Section with compression re-
inforcement
1 MEds,m ΔMEds
As1 = · + + NEd = 10.73 cm2
σs1d ζ·d d − d2
1 ΔMEds
As2 = · = 2.47 cm2
|σs2d | d − d2
1 The tools used in the design process are based on steel stress-strain diagrams, as
defined in [1] 3.2.7:(2), Fig. 3.8, which can be seen in Fig. 5.3.
2 The sections mentioned in the margins refer to DIN EN 1992-1-1:2004 (German Na-
ƒcd
As2 = · (ω2 · b · d) = 2.52 cm2
ƒyd
ρ1 = 1.021; ρ2 = 1.097
(interpolated values for d2 / d = 0.143 and ks1 = 2.740)
MEds NEd
As1 = ρ1 · ks1 · + = 10.79 cm2
d σs1d
MEds
As2 = ρ2 · ks2 · = 2.43 cm2
d
5.5 Conclusion
This example shows the calculation of the required reinforcement for a rectangular beam cross-section
under bending. Various different reference solutions are employed in order to compare the SOFiSTiK
results to. It has been shown that the results are reproduced with excellent accuracy.
5.6 Literature
[1] DIN EN 1992-1-1/NA: Eurocode 2: Design of concrete structures, Part 1-1/NA: General rules and
rules for buildings - German version EN 1992-1-1:2005 (D), Nationaler Anhang Deutschland -
Stand Februar 2010. CEN. 2010.
[2] F. Fingerloos, J. Hegger, and K. Zilch. DIN EN 1992-1-1 Bemessung und Konstruktion von
Stahlbeton- und Spannbetontragwerken - Teil 1-1: Allgemeine Bemessungsregeln und Regeln
für den Hochbau. BVPI, DBV, ISB, VBI. Ernst & Sohn, Beuth, 2012.
[3] K. Holschemacher, T. Müller, and F. Lobisch. Bemessungshilfsmittel für Betonbauteile nach Eu-
rocode 2 Bauingenieure. 3rd. Ernst & Sohn, 2012.
[4] Beispiele zur Bemessung nach Eurocode 2 - Band 1: Hochbau. Ernst & Sohn. Deutschen Beton-
und Bautechnik-Verein E.V. 2011.
[5] R. S. Narayanan and A. W. Beeby. Designers’ Guide to EN 1992-1-1 and EN 1992-1-2 - Eurocode
2: Design of Concrete Structures. Thomas Telford, 2005.
Overview
hƒ
MEd
h d
zs
As1 As1
d1
b
Figure 6.1: Problem Description
• Design stress-strain curves for concrete and reinforcement (Section 3.1.7, 3.2.7)
εc
zs
As1
εs1
Figure 6.2: Stress and Strain Distributions in the Design of T-beams
In doubly reinforced rectangular beams, the conditions in the cross-section at the ultimate limit state, are
assumed to be as shown in Fig. 6.2. The design stress-strain diagram for reinforcing steel considered
in this example, consists of an inclined top branch, as presented in Fig. 6.3 and as defined in DIN EN
1992-1-1:2004 [1] (Section 3.2.7).
A
σ ƒtk,c = 525N/ mm2
ƒyk
ƒyd = ƒyk / γs
B A Idealised
B Design
εd = 25 0 /00
Figure 6.3: Idealised and Design Stress-Strain Diagram for Reinforcing Steel
B 500A d = 60.0 cm
d1 = 5.0 cm
b = 30 cm
beƒ ƒ = 258 cm
hƒ = 18 cm
Material:
NEd = 0
MEds = MEd − NEd · zs1 = 425 kNm
d 60
Tab. 9.3 [3]: kd −Table for up to kd = p =p = 4.67
C50/ 60 - Rectangular section without MEds / b 425/ 2.58
compression reinforcement
Referring to the table for kd = 4.67 and after interpolation
defined in [1] 3.2.7:(2), Fig. 3.8, which can be seen in Fig. 6.3.
2 The sections mentioned in the margins refer to DIN EN 1992-1-1:2004 (German Na-
6.5 Conclusion
This example shows the calculation of the required reinforcement for a T-beam under bending. Two
different reference solutions are employed in order to compare the SOFiSTiK results to. It has been
shown that the results are reproduced with excellent accuracy.
6.6 Literature
[1] DIN EN 1992-1-1/NA: Eurocode 2: Design of concrete structures, Part 1-1/NA: General rules and
rules for buildings - German version EN 1992-1-1:2005 (D), Nationaler Anhang Deutschland -
Stand Februar 2010. CEN. 2010.
[2] F. Fingerloos, J. Hegger, and K. Zilch. DIN EN 1992-1-1 Bemessung und Konstruktion von
Stahlbeton- und Spannbetontragwerken - Teil 1-1: Allgemeine Bemessungsregeln und Regeln
für den Hochbau. BVPI, DBV, ISB, VBI. Ernst & Sohn, Beuth, 2012.
[3] K. Holschemacher, T. Müller, and F. Lobisch. Bemessungshilfsmittel für Betonbauteile nach Eu-
rocode 2 Bauingenieure. 3rd. Ernst & Sohn, 2012.
[4] Beispiele zur Bemessung nach Eurocode 2 - Band 1: Hochbau. Ernst & Sohn. Deutschen Beton-
und Bautechnik-Verein E.V. 2011.
[5] R. S. Narayanan and A. W. Beeby. Designers’ Guide to EN 1992-1-1 and EN 1992-1-2 - Eurocode
2: Design of Concrete Structures. Thomas Telford, 2005.
Overview
As,tot / 2 h d d2 MEd
d1 = d2 NEd
As,tot / 2
• Design stress-strain curves for concrete and reinforcement (Section 3.1.7, 3.2.7)
εc σcd
d2 σs2d Fs2d
As2 εs2 Fc
d zs2
z
zs1
As1
d1 εs1 σs1d Fs1d
Figure 7.2: Stress and Strain Distributions in the Design of Doubly Reinforced Cross-sections
In doubly reinforced rectangular beams, the conditions in the cross-section at the ultimate limit state, are
assumed to be as shown in Fig. 7.2. The design stress-strain diagram for reinforcing steel considered
in this example, consists of an horizontal top branch, as presented in Fig. 7.3 and as defined in DIN EN
1992-1-1:2004 [1] (Section 3.2.7).
A
σ ƒtk,c
ƒyk
ƒyd = ƒyk / γs
B A Idealised
B Design
Figure 7.3: Idealised and Design Stress-Strain Diagram for Reinforcing Steel
d1 = d2 = 5.0 cm
b = 30 cm
Material:
Design Load:
NEd = −1785 kN
MEd = 382 kNm
ed MEd 382
= = = 0.428 < 3.5
h NEd · h −1785 · 0.50
ωtot = 0.60
b·h
As,tot = ωtot · = 35.19 cm2
ƒyd / ƒcd
As,tot
As1 = As2 = = 17.6 cm2
2
1 The tools used in the design process are based on steel stress-strain diagrams, as
defined in [1] 3.2.7:(2), Fig. 3.8, which can be seen in Fig. 7.3.
2 The sections mentioned in the margins refer to DIN EN 1992-1-1:2004 (German Na-
7.5 Conclusion
This example shows the calculation of the required reinforcement for a rectangular beam cross-section
under bending with axial force. It has been shown that the results are reproduced with excellent accu-
racy.
7.6 Literature
[1] DIN EN 1992-1-1/NA: Eurocode 2: Design of concrete structures, Part 1-1/NA: General rules and
rules for buildings - German version EN 1992-1-1:2005 (D), Nationaler Anhang Deutschland -
Stand Februar 2010. CEN. 2010.
[2] F. Fingerloos, J. Hegger, and K. Zilch. DIN EN 1992-1-1 Bemessung und Konstruktion von
Stahlbeton- und Spannbetontragwerken - Teil 1-1: Allgemeine Bemessungsregeln und Regeln
für den Hochbau. BVPI, DBV, ISB, VBI. Ernst & Sohn, Beuth, 2012.
[3] K. Holschemacher, T. Müller, and F. Lobisch. Bemessungshilfsmittel für Betonbauteile nach Eu-
rocode 2 Bauingenieure. 3rd. Ernst & Sohn, 2012.
[4] Beispiele zur Bemessung nach Eurocode 2 - Band 1: Hochbau. Ernst & Sohn. Deutschen Beton-
und Bautechnik-Verein E.V. 2011.
Overview
MEdz
each As,tot / 4
d1 NEd
= h MEdy
d2
b
b1 = b2
• Design stress-strain curves for concrete and reinforcement (Section 3.1.7, 3.2.7)
εc σcd
d2 Fs2d
As2 εs2 σs2d
d zs2 Fc
z
zs1
As1
d1 εs1 σs1d Fs1d
Figure 8.2: Stress and Strain Distributions in the Design of Doubly Reinforced Cross-sections
The design stress-strain diagram for reinforcing steel considered in this example, consists of an inclined
top branch, as presented in Fig. 8.3 and as defined in DIN EN 1992-1-1:2004 [1] (Section 3.2.7).
A
σ
ƒtk,c
ƒyk
ƒyd = ƒyk / γs
B A Idealised
B Design
Figure 8.3: Idealised and Design Stress-Strain Diagram for Reinforcing Steel
b = 40 cm
Material:
Design Load:
NEd = −1600 kN
MEdy = 500 kNm
MEdz = 450 kNm
• ν = −0.4
b·h
As,tot = ωtot · = 113.1 cm2
ƒyd / ƒcd
1 The tools used in the design process are based on steel stress-strain diagrams, as
defined in [1] 3.2.7:(2), Fig. 3.8, which can be seen in Fig. 8.3.
2 The sections mentioned in the margins refer to DIN EN 1992-1-1:2004 (German Na-
8.5 Conclusion
This example shows the calculation of the required reinforcement for a rectangular beam cross-section
under double axial bending with compressive axial force. It has been shown that the results are repro-
duced with excellent accuracy.
8.6 Literature
[1] DIN EN 1992-1-1/NA: Eurocode 2: Design of concrete structures, Part 1-1/NA: General rules and
rules for buildings - German version EN 1992-1-1:2005 (D), Nationaler Anhang Deutschland -
Stand Februar 2010. CEN. 2010.
[2] F. Fingerloos, J. Hegger, and K. Zilch. DIN EN 1992-1-1 Bemessung und Konstruktion von
Stahlbeton- und Spannbetontragwerken - Teil 1-1: Allgemeine Bemessungsregeln und Regeln
für den Hochbau. BVPI, DBV, ISB, VBI. Ernst & Sohn, Beuth, 2012.
[3] K. Holschemacher, T. Müller, and F. Lobisch. Bemessungshilfsmittel für Betonbauteile nach Eu-
rocode 2 Bauingenieure. 3rd. Ernst & Sohn, 2012.
[4] Beispiele zur Bemessung nach Eurocode 2 - Band 1: Hochbau. Ernst & Sohn. Deutschen Beton-
und Bautechnik-Verein E.V. 2011.
Overview
As,tot / 2
h d z
VEd
As,tot / 2
• Design stress-strain curves for concrete and reinforcement (Section 3.1.7, 3.2.7)
Fcd
V
α N
d z θ
V V(cot θ − cot α) M
s Ftd
shear reinforcement
Figure 9.2: Shear Reinforced Members
The design stress-strain diagram for reinforcing steel considered in this example, consists of an inclined
top branch, as presented in Fig. 9.3 and as defined in DIN EN 1992-1-1:2004 [1] (Section 3.2.7).
A
σ
ƒtk,c
ƒyk
ƒyd = ƒyk / γs
B A Idealised
B Design
Figure 9.3: Idealised and Design Stress-Strain Diagram for Reinforcing Steel
B 500A b = 30 cm
d = 45.0 cm
VRd,cc = 0.5 · 0.48 · 30 1/ 3 · (1 − 0) · 0.3 · 0.384 (NDP) 6.2.3 (2): σcd = NEd / Ac
1.2 + 0
cot θ = = 1.60 (NDP) 6.2.3 (2): The angle θ should be
1 − 85.91 / 343.25 limited by Eq. 6.7DE
As,req / s = VEd / (ƒyd · z · cot θ) = 12.84 cm2 / m 6.2.3 (3): Eq. 6.8
ƒyd = ƒyk / γs = 435 MP
1.0 ≤ cot θ ≤ 2.5 → cot θ = 1.60 (choose for comparison) 6.2.3 (2): Eq. 6.7N
The angle θ should be limited by Eq.
As,req / s = VEd / (ƒyd · z · cot θ) = 12.18 cm2 / m 6.7N
1 The tools used in the design process are based on steel stress-strain diagrams, as
defined in [1] 3.2.7:(2), Fig. 3.8, which can be seen in Fig. 9.3.
2 The sections mentioned in the margins refer to DIN EN 1992-1-1:2004 (German Na-
wise specified.
9.5 Conclusion
This example shows the calculation of the required reinforcement for a rectangular cross-section under
shear force. It has been shown that the results are reproduced with excellent accuracy.
9.6 Literature
[1] DIN EN 1992-1-1/NA: Eurocode 2: Design of concrete structures, Part 1-1/NA: General rules and
rules for buildings - German version EN 1992-1-1:2005 (D), Nationaler Anhang Deutschland -
Stand Februar 2010. CEN. 2010.
[2] F. Fingerloos, J. Hegger, and K. Zilch. DIN EN 1992-1-1 Bemessung und Konstruktion von
Stahlbeton- und Spannbetontragwerken - Teil 1-1: Allgemeine Bemessungsregeln und Regeln
für den Hochbau. BVPI, DBV, ISB, VBI. Ernst & Sohn, Beuth, 2012.
[6] EN 1992-1-1: Eurocode 2: Design of concrete structures, Part 1-1: General rules and rules for
buildings. CEN. 2004.
Overview
beƒ ƒ
hƒ
VEd
h d
zs
As1 As1
d1
b
Figure 10.1: Problem Description
• Design stress-strain curves for concrete and reinforcement (Section 3.1.7, 3.2.7)
Fcd
V
α N
d z θ
V V(cot θ − cot α) M
s Ftd
shear reinforcement
Figure 10.2: Shear Reinforced Members
The design stress-strain diagram for reinforcing steel considered in this example, consists of an inclined
top branch, as presented in Fig. 10.3 and as defined in EN 1992-1-1:2004 [6] (Section 3.2.7).
A
σ ƒtk,c
ƒyk
ƒyd = ƒyk / γs
B A Idealised
B Design
Figure 10.3: Idealised and Design Stress-Strain Diagram for Reinforcing Steel
B 500A d = 53.0 cm
d1 = 7.0 cm
b = 30 cm
beƒ ƒ = 180 cm
hƒ = 15 cm
As1 = 15 cm2
The intermediate steps of calculating the required reinforcement are also validated in this example.
First we calculate the design value for the shear resistance VRd,c for members not requiring shear
reinforcement.
Checking the results in AQB, we can see that SOFiSTiK outputs also Vrd1,c = 62.52 kN.
Just to test this result, if we input a shear force of VEd = 62.51 kN just below the value for VRd,c , AQB
will not output any value for cot θ and the minimum reinforcement will be printed (M). If we now give a
value of VEd = 62.53 kN just larger than VRd,c , then AQB will start increasing cot θ and the minimum
reinforcement will be printed. If we continue increasing VEd , AQB will continue increasing cot θ until
it reaches the upper limit ofcot θ = 2.5 with using the minimum reinforcement. If now the minimum
reinforcement is exceeded, AQB starts calculating a value for the required reinforcement.
Another option to test this limit of VRd,c = 62.52 kN, would be to keep cot θ = 1.0 and now with
VEd = 62.53 kN, AQB calculates a value for the required reinforcement larger than the minimum
reinforcement. For the maximum value of the angle θ, hence cot θ = 1.0, the maximum value allowed
for VEd can be calculated as 755.57 kN. This can be found in AQB results as the Vrd2,c = 755.57 kN
for the case of cot θ = 1.0. Giving as an input a shear force just above this value VEd = 755.58 triggers
a warning ”Shear design not possible”.
Next step is the validation of VRd,m . When the design shear force VEd exceeds VRd,m then cot θ
must be decreased so that VEd = VRd,m . The reference result for VRd,m is 521.08 kN. Inputing a
value just below that, should give a cot θ = 2.5, whereas for a value just above should give cot θ < 2.5.
This can be verified easily in AQB output for VEd = 521.07 and 521.09 kN, respectively.
Also the minimum reinforcement is calculated exactly by AQB with a value of 2.63 cm2 / m.
SOF. Ref.
6.2.3 (2): Eq. 6.7N: The angle θ should 1.0 ≤ cot θ ≤ 2.5 → start with cot θ = 2.50
be limited by Eq. 6.7N
6.2.3 (3): Eq. 6.9 VRd,m = b · z · ν · ƒcd / (cot θ + tn θ)
ƒck
6.2.2 (6): Eq. 6.6N ν = 0.6 · 1 − = 0.528
250
Min. reinforcement:
1 The tools used in the design process are based on steel stress-strain diagrams, as
defined in [6] 3.2.7:(2), Fig. 3.8, which can be seen in Fig. 10.3.
2 The sections mentioned in the margins refer to EN 1992-1-1:2004 [6], unless other-
wise specified.
p
ρ,mn = 0.08 · ƒck / ƒyk = 0.08 · 30/ 500 = 0.0008763
p
9.2.2 (5): Eq. 9.5N
Required reinforcement:
VEd
As,req / s = 6.2.3 (3): Eq. 6.8
ƒyd · z · cot θ
α = 1.0
20
VRd,m = 1.0 · 0.3 · 0.477 · 0.60 · = 0.85859 MN
1.0 + 1.0
VRd,m = 858.59 kN
20
VRd,m = 0.3 · 0.477 · 0.528 · = 0.52108 MN
2.5 + 0.4
VRd,m = 521.08 kN ≥ VEd = 63 kN
0.063
As,req / s = · 1002 = 1.2151 cm2
434.78 · 0.477 · 2.5
• For VEd = 63 kN and cot θ = 1.0:
20
VRd,m = 0.3 · 0.477 · 0.528 · = 0.7555 MN
1.0 + 1.0
VRd,m = 755.5 kN ≥ VEd = 63 kN
0.521
As,req / s = · 1002 = 3.037 cm2
434.78 · 0.477 · 1.0
• For VEd = 756 kN and cot θ = 1.0:
Shear design not possible, because the cot θ value is fixed and can’t
be iterated.
20
VRd,m = 0.3 · 0.477 · 0.528 · = 0.52108 MN
2.5 + 0.4
VRd,m = 521.08 kN ≥ VEd = 450 kN
0.450
As,req / s = · 1002 = 8.679 cm2
434.78 · 0.477 · 2.5
10.5 Conclusion
This example shows the calculation of the required reinforcement for a T-beam under shear force. It has
been shown that the results are reproduced with excellent accuracy.
10.6 Literature
[6] EN 1992-1-1: Eurocode 2: Design of concrete structures, Part 1-1: General rules and rules for
buildings. CEN. 2004.
Overview
As,tot / 2
h d z
NEd
VEd
As,tot / 2
• Design stress-strain curves for concrete and reinforcement (Section 3.1.7, 3.2.7)
Fcd
V
α N
d z θ
V V(cot θ − cot α) M
s Ftd
shear reinforcement
Figure 11.2: Shear Reinforced Members
The design stress-strain diagram for reinforcing steel considered in this example, consists of an inclined
top branch, as presented in Fig. 11.3 and defined in DIN EN 1992-1-1:2004 [1].
A
σ
ƒyk
ƒyd = ƒyk / γs
B A Idealised
B Design
Figure 11.3: Idealised and Design Stress-Strain Diagram for Reinforcing Steel
d = 45.0 cm
cV, = 3.6 cm
SOF. Ref.
Material:
NEd = −500.0 kN
1/ 3
VRd,c = CRd,c · k · (100 · ρ1 · ƒck ) + 0.12 · σcp · b · d (NDP) 6.2.2 (1): Eq. 6.2a: Design value
for shear resistance VRd,c for members
with a minimum of not requiring design shear reinforce-
ment
1/ 2
νmn = (0.0525/ γc ) · k 3/ 2 · ƒck = 0.41249 (NDP) 6.2.2 (1): Eq. 6.3aDE: νmn for
d ≤ 600 mm
VRd,c,mn = 109.68 kN
σcp = NEd / Ac < 0.2 · ƒcd (NDP) 6.2.2 (1): Eq. 6.2
σcp > 0 for compression
σcp = −500 · 10−3 / 0.15 · 106 = −3.3333 N/ mm2 < 3.4
VRd,c = 0.1 · 1.6667 · (1.432 · 30) 1/ 3 + 0.112 · 3.3333 · 0.3 ·
0.45 = 132.71 kN
defined in [1] 3.2.7:(2), Fig. 3.8, which can be seen in Fig. 11.3.
2 The sections mentioned in the margins refer to DIN EN 1992-1-1:2004 (German Na-
σcd
1/ 3
(NDP) 6.2.3 (2): Eq. 6.7bDE VRd,cc = c · 0.48 · ƒck · 1 − 1.2 · b · z
c = 0.5 ƒcd
(NDP) 6.2.3 (2): σcd = NEd / Ac σcp = NEd / Ac
6.2.3 (3): Eq. 6.8 As,req / s = VEd / (ƒyd · z · cot θ) = 11.27 cm2 / m
ƒyd = ƒyk / γs = 435 MP
(NDP) 6.2.3 (3): Eq. 6.9 VRd,m = b · z · ν1 · ƒcd / (cot θ + tn θ)
Maximum shear force VRd,m occurs
for θ = 45◦ : cot θ = tn θ = 1 VRd,m = 0.3 · 0.384 · 0.75 · 17 / (1 + 1) = 734.4 kN
(NDP) ν1 = 0.75 · ν2 = 0.75, ν2 = 1
for ≤ C50/ 60
11.5 Conclusion
This example shows the calculation of the required reinforcement for a rectangular beam cross-section
under shear with compressive axial force. It has been shown that the results are reproduced with
excellent accuracy.
11.6 Literature
[1] DIN EN 1992-1-1/NA: Eurocode 2: Design of concrete structures, Part 1-1/NA: General rules and
rules for buildings - German version EN 1992-1-1:2005 (D), Nationaler Anhang Deutschland -
Stand Februar 2010. CEN. 2010.
[2] F. Fingerloos, J. Hegger, and K. Zilch. DIN EN 1992-1-1 Bemessung und Konstruktion von
Stahlbeton- und Spannbetontragwerken - Teil 1-1: Allgemeine Bemessungsregeln und Regeln
für den Hochbau. BVPI, DBV, ISB, VBI. Ernst & Sohn, Beuth, 2012.
Overview
As,tot / 2
h d
TEd
VEd
As,tot / 2
• Design stress-strain curves for concrete and reinforcement (Section 3.1.7, 3.2.7)
• Guidelines for shear (Section 6.2) and torsion design (Section 6.3)
z
Ak
cover centre-line
TEd
teƒ / 2 V
teƒ
The design stress-strain diagram for reinforcing steel considered in this example, consists of an inclined
top branch, as presented in Fig. 12.3 and as defined in DIN EN 1992-1-1:2004 [1] (Section 3.2.7).
A
σ ƒtk,c
ƒyk
ƒyd = ƒyk / γs
B A Idealised
B Design
Figure 12.3: Idealised and Design Stress-Strain Diagram for Reinforcing Steel
d = 65.0 cm
SOF. Ref.
Material:
(NDP) 3.1.6: (1)P, Eq. (3.15): cc = ƒcd = cc · ƒck / γc = 0.85 · 35/ 1.5 = 19.833 MP
0.85 considering long term effects
(NDP) 3.2.2: (3)P: yield strength ƒyk = ƒyk = 500 MP
500 MP
3.2.7: (2), Fig. 3.8 ƒyd = ƒyk / γs = 500/ 1.15 = 434.78 MP
Design Load:
VRd,c = 111.74 kN
4.5 · 35
175 · 1 + = 700 > 111.74
175 · 0.3
→ requirement of Eq. NA.6.31.2 is not met
1 The tools used in the design process are based on steel stress-strain diagrams, as
defined in [1] 3.2.7:(2), Fig. 3.8, which can be seen in Fig. 12.3.
2 The sections mentioned in the margins refer to DIN EN 1992-1-1:2004 (German Na-
Torsional reinforcement
teƒ ƒ ,1 = teƒ ƒ ,2 = 2 · 50 = 100 mm (so = s = ss = 50 mm) 6.3.1 (3): Solid sections may be mod-
elled by equivalent thin-walled sections
Ak = (h − s − so ) · b − teƒ ƒ ,1 = 100 mm
(Fig. 12.2)
Ak : area enclosed by the centre-line
Ak = (700 − 50 − 50) · (300 − 100) = 120000 mm2 = 0.12 m2
As,req / s = TEd · tn θ / (ƒyd · 2Ak ) (NDP) 6.3.2 (3): Eq. (NA.6.28.1)
As,req = TEd · k · cot θ / (ƒyd · 2Ak ) = 5.37 cm2 (T) 6.3.2 (3): Eq. 6.28
Shear reinforcement
As,req / s = VEd / (ƒyd · z · cot θ) = 6.89 cm2 / m 6.2.3 (3): Eq. 6.8
ƒyd = ƒyk / γs = 435 MP
12.5 Conclusion
This example shows the calculation of the required reinforcement for a rectangular beam cross-section
under shear and torsion. It has been shown that the results are reproduced with excellent accuracy.
12.6 Literature
[1] DIN EN 1992-1-1/NA: Eurocode 2: Design of concrete structures, Part 1-1/NA: General rules and
rules for buildings - German version EN 1992-1-1:2005 (D), Nationaler Anhang Deutschland -
Stand Februar 2010. CEN. 2010.
[2] F. Fingerloos, J. Hegger, and K. Zilch. DIN EN 1992-1-1 Bemessung und Konstruktion von
Stahlbeton- und Spannbetontragwerken - Teil 1-1: Allgemeine Bemessungsregeln und Regeln
für den Hochbau. BVPI, DBV, ISB, VBI. Ernst & Sohn, Beuth, 2012.
[4] Beispiele zur Bemessung nach Eurocode 2 - Band 1: Hochbau. Ernst & Sohn. Deutschen Beton-
und Bautechnik-Verein E.V. 2011.
Overview
beƒ ƒ
hƒ
VEd
h d
As1 As1
d1
• Design stress-strain curves for concrete and reinforcement (Section 3.1.7, 3.2.7)
Fd
Fd beƒ ƒ
sƒ Δ
θƒ
The design stress-strain diagram for reinforcing steel considered in this example, consists of an inclined
top branch, as presented in Fig. 13.3 and as defined in DIN EN 1992-1-1:2004 [1] (Section 3.2.7).
A
σ ƒtk,c = 525N/ mm2
ƒyk
ƒyd = ƒyk / γs
B A Idealised
B Design
εd = 25 0 /00
Figure 13.3: Idealised and Design Stress-Strain Diagram for Reinforcing Steel
175
1 1
15
1 1
60
67.5 40 67.5
d1 = 7.0 cm
0.00
SOFiSTiK AG * Bruckmannring 38 * 85764 Oberschleißheim Page 1
SOFiSTiK 2016-3 WINGRAF - GRAPHICS FOR FINITE ELEMENTS (V 18.03) 2016-01-15
Single span beam
The system with its loading as well as the moment and shear force are shown in Fig. 13.5. The reference
calculation steps [4] 2.00
are presented in the4.00next section and the6.00results are given in8.00Table 13.2.
SOFiSTiK AG - www.sofistik.de
0.00 10.00 m
All loads (in components), Loadcase 1 sum_PZ=1550.0 kN . , (1 cm 3D = unit) Beam line load (force) in global Z (Unit=160.2 M 1 : 41
Y X
kN/m ) (Max=155.0)
Z
155.0 155.0 155.0 155.0 155.0 155.0 155.0 155.0 155.0 155.0 0.00
0.00
1001 1002 1003 1004 1005 1006 1007 1008 1009 1010
697.5
697.5
1240
1240
1628
1628
SOFiSTiK AG - www.sofistik.de
1860
1938
All loads (in components), Loadcase 1 sum_PZ=1550.0 kN . , (1 cm 3D = unit) Beam line load (force) in global Z (Unit=160.2 M 1 : 41
Y X
kN/m ) (Max=155.0)
Z
697.5
1240
1628
1628
1860
1860
1938
Asƒ / sƒ [cm2 / m]
0.00 2.00 4.00 6.00 8.00 10.00 m
5.72 5.79
Beam Elements , Number of element (Max=1010) M 1 : 41
Y X
Beam Elements , Bending moment My, Loadcase 1 sum_PZ=1550.0 kN . , 1 cm 3D = 1000. kNm (Max=1938.)
Z
Material:
Design loads
ω ≈ 0.063, ζ ≈ 0.967 and ξ ≈ 0.086 (interpolated) Tab. 9.2 [3]: ω−Table for up to C50/ 60
without compression reinforcement
1
As1 = · (ω · b · d · ƒcd + NEd ) NEd = 0
σsd
1
As1 = · (0.063 · 1.75 · 0.68 · 14.17) · 1002 = 23.27 cm2
456.52
z = ζ · d = 0.967 · 0.68 m ≈ 65.7 cm
1 The tools used in the design process are based on steel stress-strain diagrams, as
defined in [1] 3.2.7:(2), Fig. 3.8, which can be seen in Fig. 13.3.
2 The sections mentioned in the margins refer to DIN EN 1992-1-1:2004 (German Na-
ΔFd
DIN EN 1992, 6.2.4 (3),Eq. 6.20: hƒ Ed =
is the thickness of flange at the junc- hƒ · Δ
tions, Δ is the length under consider-
ation, ΔFd is the change of the normal In our case Δ = 1.0 because the beam length is = 1.00 m.
force in the flange over the length Δ
Please note that AQB is outputting the results per length.
409.36
In AQB output Ed = τ − V Ed = = 0.272 kN/ m2 = 2.72 MP
15 · 100
Checking the maximum Rd,m value to prevent crushing of the
struts in the flange
DIN EN 1992, 6.2.4 (4),Eq. 6.22 Ed ≤ Rd,m = ν · ƒcd · snθƒ · cosθƒ
ν = ν1
b → hƒ , z → Δ, c = 0.5
σcd
1/ 3
VRd,cc = c · 0.48 · ƒck · 1 − 1.2 · · hƒ · Δ
ƒcd
0
VRd,cc = 0.5 · 0.48 · 251/ 3 · 1 − 1.2 · · 0.15 · 1.0
14.17
VRd,cc = 0.1052 MN = 105.26 kN
1.2
cotθ = = 1.619
1 − 105.26/ 409.36
1 1
tn θ = = = 0.619 → θ = 31.75◦
cotθ 1.619
Rd,m = 0.75 · 14.17 · sn31.75 · cos31.75 = 4.7502 MP
VRd,m = 712.53 kN
Rd,c = k · ƒctd
Ed · hƒ
sƒ = DIN EN 1992, 6.2.4 (4),Eq. 6.21
cotθƒ · ƒyd
2.72 · 0.15
sƒ = · 1002 = 5.79 cm2
1.619 · 434.78
13.5 Conclusion
This example is concerned with the calculation of the shear between web and flanges of T-sections. It
has been shown that the results are reproduced with good accuracy.
13.6 Literature
[1] DIN EN 1992-1-1/NA: Eurocode 2: Design of concrete structures, Part 1-1/NA: General rules and
rules for buildings - German version EN 1992-1-1:2005 (D), Nationaler Anhang Deutschland -
Stand Februar 2010. CEN. 2010.
[2] F. Fingerloos, J. Hegger, and K. Zilch. DIN EN 1992-1-1 Bemessung und Konstruktion von
Stahlbeton- und Spannbetontragwerken - Teil 1-1: Allgemeine Bemessungsregeln und Regeln
für den Hochbau. BVPI, DBV, ISB, VBI. Ernst & Sohn, Beuth, 2012.
[3] K. Holschemacher, T. Müller, and F. Lobisch. Bemessungshilfsmittel für Betonbauteile nach Eu-
rocode 2 Bauingenieure. 3rd. Ernst & Sohn, 2012.
[4] Beispiele zur Bemessung nach Eurocode 2 - Band 1: Hochbau. Ernst & Sohn. Deutschen Beton-
und Bautechnik-Verein E.V. 2011.
Overview
beƒ ƒ
hƒ
VEd
h d
Shear section
zs
As1 As1
d1
b
Figure 14.1: Problem Description
• Design stress-strain curves for concrete and reinforcement (Section 3.1.7, 3.2.3)
≤ 30◦ VEd b
h1 ≤ 10d
old concrete anchorage
The design stress-strain diagram for reinforcing steel considered in this example, consists of an inclined
top branch, as presented in Fig. 14.3 and as defined in DIN EN 1992-1-1:2004 [1] (Section 3.2.7).
A
σ ƒtk,c
ƒyk
ƒyd = ƒyk / γs
B A Idealised
B Design
Figure 14.3: Idealised and Design Stress-Strain Diagram for Reinforcing Steel
d1 = 7.0 cm
b = 40 cm , beƒ ƒ = 250 cm
zs = 95.56 cm
Material:
S = h · b · (zs − h / 2) = 0.18058 m3
0.8 · 0.18058
τ= = 2.1669 MP
0.16667 · 0.4
0.8 · 0.18058
T = = 0.86676 MN/ m = 866.76 kN/ m
0.16667
The shear section with a length of 0.4
m is split into two equal parts with b =
0.2 m T = 866.76 / 2 = 433.38 kN/ m
State :
As
ρ1 = = 0.0 → Rd,c = 0.0
b d
with a minimum of
defined in [1] 3.2.7:(2), Fig. 3.8, which can be seen in Fig. 14.3.
2 The sections mentioned in the margins refer to DIN EN 1992-1-1:2004 (German Na-
c = 0.50 and μ = 0.9 for indented surface 6.2.5 (2): c, μ: factors depending on the
roughness of the interface
ƒctd = αct · ƒctk;0.05 / γc (NDP) 3.1.6 (2)P: Eq. 3.16
αct = 0.85
ƒctd = 0.85 · 1.80 / 1.5 = 1.02 3.1.2 (3): Tab. 3.1 - Strength for con-
crete: ƒctk;0.05 = 1.8 MP:
s s
Rd = 0.5 · 1.02 + 0 + · 435 · (1.2 · 0.9 · 1 + 0) ρ= : area of reinforcement cross-
0.2 · 1.0 b ·
ing interface / area of joint
s
Rd = 0.51 + · 469.56 = 2.1669
0.2
s = 7.07 cm2 / m
Design Load:
From the calculated inner lever arms for the two states we get a ratio:
z
= 0.7664
z
The associated design shear flow VEd is:
Ed ≤ Rd
Tab. 9.2 [3]: ω−Table for up to C50/ 60 ω = 0.03971 and ξ = 0.9766 (interpolated)
without compression reinforcement
1
NEd = 0 As1 = · (ω · b · d · ƒcd + NEd ) = 39.44 cm2
σsd
Design Load:
The shear section with a length of 0.4 T = V / z = 800 / 1.19 = 672.268 kN/ m
m is split into two equal parts with b =
0.2 m T = 672.268 / 2 = 336.134 kN/ m
1/ 3
Rd,c = CRd,c · k · (100 · ρ1 · ƒck ) + 0.12 · σcp
14.5 Conclusion
This example shows the calculation of the required reinforcement for a T-section under shear at the
interface between concrete cast at different times. It has been shown that the results are reproduced
with excellent accuracy. Small deviations occur because AQUA calculates (by using FEM analysis) the
shear stresses more accurate compared to the reference example.
14.6 Literature
[1] DIN EN 1992-1-1/NA: Eurocode 2: Design of concrete structures, Part 1-1/NA: General rules and
rules for buildings - German version EN 1992-1-1:2005 (D), Nationaler Anhang Deutschland -
Stand Februar 2010. CEN. 2010.
[2] F. Fingerloos, J. Hegger, and K. Zilch. DIN EN 1992-1-1 Bemessung und Konstruktion von
Stahlbeton- und Spannbetontragwerken - Teil 1-1: Allgemeine Bemessungsregeln und Regeln
für den Hochbau. BVPI, DBV, ISB, VBI. Ernst & Sohn, Beuth, 2012.
[3] K. Holschemacher, T. Müller, and F. Lobisch. Bemessungshilfsmittel für Betonbauteile nach Eu-
rocode 2 Bauingenieure. 3rd. Ernst & Sohn, 2012.
[4] Beispiele zur Bemessung nach Eurocode 2 - Band 1: Hochbau. Ernst & Sohn. Deutschen Beton-
und Bautechnik-Verein E.V. 2011.
Overview
As0
h d z MEd
NEd
As
• Design stress-strain curves for concrete and reinforcement (Section 3.1.7, 3.2.7)
• Basic assumptions for calculation of crack widths (Section 7.3.2, 7.3.3, 7.3.4)
A
ϕ h−
c B
ε2 = 0
h d
A C
E D
hc,eƒ
B ε1 5(c + ϕ/ 2)
The design stress-strain diagram for reinforcing steel considered in this example, consists of an inclined
top branch, as defined in DIN EN 1992-1-1:2004 [1] (Section 3.2.7).
b = 30.0 cm
k = 0.3 mm
given reinforcement
given reinforcement
Material:
3.1.6: (1)P, Eq. (3.15): cc = 0.85 con- ƒcd = cc · ƒck / γc = 0.85 · 25/ 1.5 = 14.17 MP
sidering long term effects
3.2.7: (2), Fig. 3.8 ƒyd = ƒyk / γs = 500/ 1.15 = 434.78 MP
Design Load:
7.3.2 (2): Eq. 7.1: As,mn minimum area As,mn · σs = kc · k · ƒct,eƒ ƒ · Act
of reinforcing steel within tensile zone
(NDP) 7.3.2 (2): ƒct,eƒ ƒ mean value of ƒct,eƒ ƒ = ƒctm ≥ 3.0 MP
concrete tensile strength
When the cracking time can’t be placed ƒct,eƒ ƒ = 2.6 < 3.0 MP → ƒct,eƒ ƒ = 3.0 MP
with certainty in the first 28 days then
ƒct,eƒ ƒ ≥ 3.0 MP for normal concrete
Tab. 3.1: ƒctm = 2.6 MP for C 25
Act = b · heƒ ƒ = 0.3 · 0.5 = 0.15 m2
7.3.2 (2): Act is the area of concrete
within tensile zone, for pure bending of
a rectangular section its height is h/ 2 fct,eff
heƒ ƒ = h/ 2
Ç
σs = k · 3.48 · 106 / ϕ∗
s
= 207.846 MP (NDP) 7.3.3 Tab. 7.2DE (a)
where k = 0.3 m the prescribed
maximum crack width
As,req = 0.4 · 0.8 · 3.0 · 0.15 · 104 / 207.846 = 6.928cm2
As,mn · σs = kc · k · ƒct,eƒ ƒ · Act 7.3.2 (2): Eq. 7.1: As,mn minimum area
of reinforcing steel within tensile zone
ƒct,eƒ ƒ = ƒctm ≥ 3.0 MP
ƒct,eƒ ƒ = 2.6 < 3.0 MP → ƒct,eƒ ƒ = 3.0 MP Tab. 3.1: ƒctm = 2.6 MP for C 25
As,mn · σs = kc · k · ƒct,eƒ ƒ · Act 7.3.2 (2): Eq. 7.1: As,mn minimum area
of reinforcing steel within tensile zone
ƒct,eƒ ƒ = ƒctm ≥ 3.0 MP
ƒct,eƒ ƒ = 2.6 < 3.0 MP → ƒct,eƒ ƒ = 3.0 MP Tab. 3.1: ƒctm = 2.6 MP for C 25
3.0 · 50
⇒ heƒ ƒ = = 37.5cm
3.0 + 1.0
ϕs = ϕ∗
s
· ƒct,eƒ ƒ / 2.9 N/ mm2
ϕ∗
s
= 25 · 2.9 / 3.0 = 24.16667
Ç
(NDP) 7.3.3 Tab. 7.2DE (a) σs = k · 3.48 · 106 / ϕ∗
s
= 207.846 MP
where k = 0.3 m the prescribed
maximum crack width σc
7.3.2 (2): Eq. 7.2: kc for bending with kc = 0.4 · 1 − ≤1
axial force of rectangular sections k1 · (h/ h∗ ) · ƒct,eƒ ƒ
7.3.3: Control of cracking without direct Control of cracking without direct calculation
calculation
Examples calculated in this section are • Case V: NEd = 0.0
w.r.t. Table 7.2DE, here Table 7.3N is
not relevant
(NDP) 7.3.2 (2): ƒct,eƒ ƒ mean value of
ƒct,eƒ ƒ = ƒctm ≥ 3.0 MP
concrete tensile strength
Tab. 3.1: ƒctm = 2.6 MP for C 25
ƒct,eƒ ƒ = 2.6 < 3.0 MP → ƒct,eƒ ƒ = 3.0 MP
(NDP) 7.3.3: Eq. 7.7.1DE: The maxi-
σs · As ƒct,eƒ ƒ
mum bar diameters should be modified
ϕs = ϕ∗
s
· ≥ ϕ∗
s
·
for load action 4(h − d) · b · 2.9 2.9
As = 24.5 cm2 prescribed reinforce-
ment 264.06 · 24.50
ϕs = 25 mm = ϕ∗
s
· = ϕ∗
s
· 4.6476
4(100 − 96) · 30 · 2.9
→ ϕ∗
s
= 5.3791 mm
Ç
(NDP) 7.3.3 Tab. 7.2DE (a) σs = k · 3.48 · 106 / ϕ∗
s
= 440.53 MP
where k = 0.3 m the prescribed
maximum crack width σs = 264.06 < 440.53 MP
In case the usage factor becomes 1.0 then the stresses σs are equal,
as it can be seen in Case V below.
ƒct,eƒ ƒ 3.0
and ϕs = 25 mm > ϕ∗
s
· = 5.3791 · = 5.5646
2.9 2.9
also control the steel stress with respect to the calculated strains
or
(NDP) 7.3.3 Tab. 7.2DE (a) from Tab. 7.2DE and for σs = 264.04 ≈ 260.0 MP
σs · As ƒct,eƒ ƒ
(NDP) 7.3.3: Eq. 7.7.1DE: The maxi- → ϕ∗
s
= 15 mm → ϕs = ϕ∗
s
· ≥ ϕ∗
s
·
mum bar diameters should be modified 4(h − d) · b · 2.9 2.9
for load action
264.06 · 24.50
ϕs = 15 · = 69.69 mm > 25 mm
4(100 − 96) · 30 · 2.9
→ crack width control is passed with given reinforcement.
ƒct,eƒ ƒ = ƒctm ≥ 3.0 MP (NDP) 7.3.2 (2): ƒct,eƒ ƒ mean value of
concrete tensile strength
ƒct,eƒ ƒ = 2.6 < 3.0 MP → ƒct,eƒ ƒ = 3.0 MP Tab. 3.1: ƒctm = 2.6 MP for C 25
from Tab. 7.2DE and for σs = 509.15 ≈ 510.0 MP (NDP) 7.3.3 Tab. 7.2DE (a)
σs · As
→ ϕ∗
s
≈ 3.9 mm → ϕs = ϕ∗
s
· (NDP) 7.3.3: Eq. 7.7.1DE: The maxi-
4(h − d) · b · 2.9 mum bar diameters should be modified
for load action
509.15 · 12.0
ϕs = 3.9 · = 17.12 mm < 25 mm
4(100 − 96) · 30 · 2.9
→ crack width control is not passed with given reinforcement.
from Tab. 7.2DE and for σs = 436.43 MP (NDP) 7.3.3 Tab. 7.2DE (a)
σs · As
→ ϕ∗
s
≈ 5.6 mm → ϕs = ϕ∗
s
· 7.3.3: Eq. 7.7.1DE: The maximum bar
4(h − d) · b · 2.9 diameters should be modified for load
action
436.43 · 14.54
ϕs = 5.6 · = 25.45 mm ≥ 25 mm
4(100 − 96) · 30 · 2.9
→ crack width control passed with additional reinforcement.
436.36 · 14.54
ϕs = 25 mm = ϕ∗
s
·
4(100 − 96) · 30 · 2.9
→ ϕ∗
s
= 5.4849 mm
Ç
σs = k · 3.48 · 106 / ϕ∗
s
(NDP) 7.3.3 Tab. 7.2DE (a)
where k = 0.3 m the prescribed
p maximum crack width
σs = 0.3 · 3.48 · 106 / 5.4849 = 436.28 MP
Here we can notice that the stresses are equal leading to a usage
factor of 1.0
7.3.4 (1): Eq. 7.10: where A0p and ρp,eƒ ƒ = (As + ξ21 · A0p ) / Ac,eƒ ƒ
Ac,eƒ ƒ are defined in 7.3.2 (3)
7.3.2 (3): where A0p is the area of pre or A0p = 0.0 cm2
post-tensioned tendons within Ac,eƒ ƒ
ρp,eƒ ƒ = 24.50 / 390 = 0.06282
ƒct,eƒ ƒ
σs − kt · · (1 + αe · ρp,eƒ ƒ )
ρp,eƒ ƒ σs
7.3.4 (1): Eq. 7.9: the difference of the εsm − εcm = ≥ 0.6 ·
mean strain in the reinforcement and in Es Es
the concrete
2/ 3
(NDP) 7.3.2 (2): ƒct,eƒ ƒ mean value of ƒct,eƒ ƒ = ƒctm = 0.30 · ƒck = 2.565 ≈ 2.6 MP
concrete tensile strength, here no mini-
mum value of ƒct,eƒ ƒ ≥ 3.0 MP is set ƒct,eƒ ƒ = 2.6 < 3.0 MP → ƒct,eƒ ƒ = 3.0 MP
Tab. 3.1: ƒctm = 2.6 MP for C 25 or
2/ 3
ƒctm = 0.30 · ƒck = 2.565 MP 2.565
264.06 − 0.4 · · (1 + 6.354 · 0.06282)
0.06282
kt : factor dependent on the duration of εsm − εcm =
the load, kt = 0.4 for long term loading 200000
264.06
εsm − εcm = 1.2060 · 10−3 > 0.6 · = 0.79218 · 10−3
200000
ϕ σs · ϕ
(NDP) 7.3.4 (3): Eq. 7.11: sr,m is the sr,m = ≤
maximum crack spacing 3.6 · ρp,eƒ ƒ 3.6 · ƒct,eƒ ƒ
25
sr,m = = 110.545 mm
3.6 · 0.06282
σs · ϕ
sr,m ≤ = 611.25 mm
3.6 · ƒct,eƒ ƒ
7.3.4 (1): Eq. 7.8: k the crack width k = sr,m · (εsm − εcm ) = 110.545 · 1.2060 · 10−3
k = 0.133 < 0.3 mm → Check for crack width passed with given
reinforcements
Stress limitation
15.5 Conclusion
This example shows the calculation of crack widths. Various ways of reference calculations are demon-
strated, in order to compare the SOFiSTiK results to. It has been shown that the results are reproduced
with excellent accuracy.
15.6 Literature
[1] DIN EN 1992-1-1/NA: Eurocode 2: Design of concrete structures, Part 1-1/NA: General rules and
rules for buildings - German version EN 1992-1-1:2005 (D), Nationaler Anhang Deutschland -
Stand Februar 2010. CEN. 2010.
[2] F. Fingerloos, J. Hegger, and K. Zilch. DIN EN 1992-1-1 Bemessung und Konstruktion von
Stahlbeton- und Spannbetontragwerken - Teil 1-1: Allgemeine Bemessungsregeln und Regeln
für den Hochbau. BVPI, DBV, ISB, VBI. Ernst & Sohn, Beuth, 2012.
Overview
b
tƒ
r
t
Mz
Vy
h y
presented below and the results are given in Table 16.2. The utilisation level of allowable plastic forces
are calculated and compared to SOFiSTiK results.
b = 20.0 cm Mz = 20 kNm
h = 19.0 cm
tƒ = 1.0 cm, t = 0.65 cm
r = 1.8 cm
Where the shear force is present allowance should be made for its effect
on the moment resistance. This effect may be neglected, where the 6.2.8: Bending and shear
shear force is less than half the plastic shear resistance.
p
AV · ƒy / 3
Vp,Rd = 6.2.6 (2): Eq. 6.18: Vp,Rd the design
γM0 plastic shear resistance
AV the shear area w.r.t. y− and z−axis,
AVy = 2 · Aƒ nge (in the y−direction only contribution the two flanges ) respectively
AVy = 2 · tƒ · b = 2 · 1 · 20 = 40 cm2
AVz = 53.8 − 2 · 20 · 1 + (0.65 + 2 · 1.8) · 1 6.2.6 (3): The shear area AV may be
taken as follows for rolled I-sections with
h = h − 2 · tƒ = 17 cm load parallel to the web
AVz = 18.08 cm2 > 1 · 17 · 0.65 = 11.05 6.2.6 (3): η may be conservatively taken
equal to 1.0
p
40 · 23.5 / 3 6.1 (1): Partial factor γM0 = 1.00 is
Vp,Rd,y = = 542.70 kN recommended
1.00
p
18.08 · 23.5 / 3
Vp,Rd,z = = 245.30 kN
1.00
• Case :
Vp,Rd,y = 611.85 kN
Vp,Rd,z = 228.38 kN
we have:
wise specified.
VEd 200
= = 0.3268 < 0.5
Vp,Rd,y 611.85
• Case :
VEd 300
= = 0.490 < 0.5
Vp,Rd,y 611.85
Mz 20
= = 0.418
Mz,V,Rd 47.89
b) EC3 Tables:
Vp,Rd,y = 542.71 kN
Vp,Rd,z = 245.32 kN
• Case :
VEd 200
= = 0.368 < 0.5
Vp,Rd,y 542.71
• Case :
VEd 300
= = 0.552 > 0.5
Vp,Rd,y 542.71
6.2.8 (3): The reduced moment resis- → reduction of moment resistance due to shear
tance should be taken as the design re-
sistance of the cross-section, calculated
using a reduced yield strength (1 − ρ) ·
ƒy for the shear area
100 VERiFiCATiON MANUAL - Design Code Benchmarks | SOFiSTiK 2018
DCE-EN13: Design of a Steel I-section for Bending and Shear
2
2 · VEd
ρ= −1
Vp,Rd,y
2
2 · 300
= −1
542.71
= 0.011143
Wp,z · ƒy
Mz,V,Rd = (1 − ρ) ·
γM0
Mz 20
= = 0.4223
Mz,V,Rd 47.356
16.5 Conclusion
This example shows the calculation of the resistance of steel cross-section for bending and shear. It has
been shown that the results are reproduced with excellent accuracy.
y Vy Vy
Figure 16.2: Visualisation of areas and approach for calculating the Vp,Rd,y
Vz Vz
Figure 16.3: Visualisation of areas and approach for calculating the Vp,Rd,z
By using FEM analysis the Vp,Rd,y and Vp,Rd,z values are deviating from the reference values. The
reason behind this difference is shown in Fig. 16.2, 16.3, 16.4 and 16.5.
FEM analysis represents a real physical behaviour, because as shown in Fig. 16.2 and 16.3 the calcu-
lated area is taken approximately for the reference Vp,Rd,y and Vp,Rd,z values. The results by using
FEM analysis are more accurate in reality than the reference example.
SOFiSTiK AG * Bruckmannring 38 * 85764 Oberschleißheim SOFiSTiK AG * Bruckmannring 38 * 85764 Oberschleißheim Page 1 Page 1
SOFiSTiK 2018 WINGRAF - GRAPHICS FOR FINITE ELEMENTS (V 19.00) SOFiSTiK 2018 WINGRAF - GRAPHICS FOR FINITE ELEMENTS2016-06-03
(V 19.00) 2016-06-03
0.375
0.366
0.356
0.347
0.337
0.328
0.319
0.309
0.300
-0.05
0.291
-0.05
0.281
0.272
0.262
0.253
0.244
0.234
0.225
0.216
SOFiSTiK AG - www.sofistik.de
0.206
0.197
0.00
0.00
0.187
0.178
0.169
0.159
0.150
0.141
0.131
0.122
0.112
0.103
0.05
0.05
0.094
0.084
0.075
0.066
0.056
0.047
0.037
0.028
0.019
0.009
0.15 0.10 0.05 0.00 -0.05 0.15 -0.10 0.10 -0.15 0.05 m 0.00 -0.05 -0.10 m
0.000
Shear stress of the Cross section in Element, Loadcase 2 Vy = 1.0 , 1 cm 3D = 0.866 MPa (Max=0.377)
Shear stress of the Cross section from middle of element, LoadcaseM21Vy
: 1.23
Y X
Z
Figure 16.4: Vp,Rd,y - SOFiSTiK Results (FEM) Y X
Z
= 1.0 , from 1.6442e-04 to 0.375 step 0.0094 MPa M 1 : 1.23
SOFiSTiK AG * Bruckmannring 38 * 85764 Oberschleißheim SOFiSTiK AG * Bruckmannring 38 * 85764 Oberschleißheim Page 1 Page 1
SOFiSTiK 2018 WINGRAF - GRAPHICS FOR FINITE ELEMENTS (V 19.00) SOFiSTiK 2018 WINGRAF - GRAPHICS FOR FINITE ELEMENTS 2016-06-03
(V 19.00) 2016-06-03
0.895
0.870
0.848
0.825
0.803
0.781
0.759
0.736
0.714
-0.05
0.692
-0.05
0.669
0.647
0.625
0.602
0.580
0.558
0.535
0.513
SOFiSTiK AG - www.sofistik.de
0.491
0.469
0.00
0.00
0.446
0.424
0.402
0.379
0.357
0.335
0.312
0.290
0.268
0.245
0.05
0.05
0.223
0.201
0.178
0.156
0.134
0.112
0.089
0.067
0.045
0.022
0.15 0.10 0.05 0.00 0.15
-0.05 0.10
-0.10 0.05
-0.15 m 0.00 -0.05 -0.10 m
0.002
Y X
Z
Shear stress of the Cross section in Element, Loadcase 3 Vz = 1.0 , 1 cm 3D = 2.16 MPa
Figure 16.5: Vp,Rd,z - SOFiSTiK Results (FEM) Y X
Z
Shear stress of the Cross section from middle of element, Loadcase M
(Max=0.895) 3 Vz
1 : =1.23
1.0 , from 0.0024 to 0.895 step 0.0223 MPa M 1 : 1.23
16.6 Literature
[7] EN 1993-1-1: Eurocode 3: Design of concrete structures, Part 1-1: General rules and rules for
buildings. CEN. 2005.
[8] Schneider. Bautabellen für Ingenieure. 21th. Bundesanzeiger Verlag, 2014.
Overview
tƒ
r
t
h y
A diagrammatic representation of the four classes of section is given in Fig. 17.2, where a cross-section
is subjected to an increasing major axis bending moment until failure [9].
Bending Moment
Curvature
Figure 17.2: Idealized Moment Curvature Behaviour for Four Classes of Cross-sections
h = 462.0 mm
tƒ = 17.0 mm
t = 9.6 mm
r = 10.2 mm
ε = 235/ ƒy = 0.924
Æ
Tab. 5.2: Maximum width-to-thickness
ratios for compression parts
The role of cross-section classification is to identify the extent to which
the resistance and rotation capacity of cross-sections is limited by its
Tab. 5.5.1(1): Classification of cross- local buckling resistance.
section basis
• Bending:
Tab. 5.2 (sheet 1): Internal compres- For the web: c = h − 2tƒ − 2r
sion parts
c / t = (462 − 2 · 17 − 2 · 10.2) / 9.6 = 42.46
Class 1 cross-sections are those which can form a plastic hinge with
the rotation capacity required from plastic analysis without reduction of
5.5.2 (1): Classification the resistance.
• Compression:
Tab. 5.5.2 (6): A cross-section is clas- Overall classification for bending: Class 4
sified according to the highest least
favourable class Class 4 cross-sections are those in which local buckling will occur be-
fore the attainment of yield stress in one or more parts of the cross-
5.5.2 (1): Classification section.
1 The sections mentioned in the margins refer to EN 1993-1-1:2005 [7] unless other-
wise specified.
17.5 Conclusion
This example shows the classification of steel cross-sections for bending and compression. It has been
shown that the results are reproduced with excellent accuracy.
17.6 Literature
[7] EN 1993-1-1: Eurocode 3: Design of concrete structures, Part 1-1: General rules and rules for
buildings. CEN. 2005.
[9] Structural Eurocodes - Extracts from the structural Eurocodes for students of structural design.
BSI - British Standards Institution. 2007.
Overview
My
b
tƒ
r
t
h y
My
t = 9.6 mm
r = 10.2 mm
A = 94.48 cm2
z = 1046.5 cm4
T = 66.23 cm4
= 516297.12 cm6
SOF. Ref.
Mcr [kNm] 154.05 154.26
ε = 235/ ƒy = 0.924
Æ
Tab. 5.2: Maximum width-to-thickness
ratios for compression parts
Design Load:
MEd = 150kNm
Wp,y · ƒy
Mc,Rd = Mp,Rd,y = = 447.31 kNm 6.2.5 (2):Eq. 6.13: Bending resistance
γM0 Mc,Rd for Class 1 cross-section
p v 6.1 (1): γM0 = 1.00
π 2 E
π Ez GT u
Mcr = 1+
t
6.3.2.2 (2): Mcr , the elastic critical mo-
L GT L2 ment for ltb is based on gross cross sec-
tional properties
p v
3.142 · 108.42
3.14 · 2197.74 · 53.496 u Various empirical or approximate formu-
Mcr = 1+
t
lae exist for the determination of Mcr
8 53.496 · 82
2
h i
LT = 0.5 1 + αLT λLT − λLT,0 + βλLT 6.3.2.3 (1): LT for rolled sections in
bending. Recommended values:
λLT,0 = 0.4 (maximum value)
h / b = 462/ 154.4 = 2.99 > 2
β = 0.75 (minimum value)
6.3.2.3 (1): Table 6.5: Recommenda-
for rolled I-sections and h / b > 2 → buckling curve c tion for the selection of ltb curve for
cross-sections using Eq. 6.57
for buckling curve c → αLT = 0.49 6.3.2.2 (2): Table 6.5: Recommenda-
tion values for imperfection factors αLT
LT = 0.5 1 + 0.49 1.703 − 0.4 + 0.75 · 1.7032 = 1.907
for ltb curves
1
χLT = r = 0.321 6.3.2.3 (1): Eq. 6.57: χLT reduction fac-
2
LT + 2LT − β λLT tor for ltb
1
but χLT = 0.321 ≤ 1.0 and χLT = 0.321 ≤ 2
= 0.345
λLT
ƒy
Mb,Rd = χLT Wy = 143.587 kNm 6.3.2.1 (3): Eq. 6.55: Mb,Rd The design
γM1 buckling resistance moment of laterally
unrestrained beam
MEd For Class 1 section Wy = Wp,y
= 1.045 → Beam fails in LTB
Mb,Rd 6.1 (1): γM1 = 1.00 recommended
value
1 The sections mentioned in the margins refer to EN 1993-1-1:2005 [7] unless other-
wise specified.
18.5 Conclusion
This example shows the ckeck for lateral torsional buckling of steel members. It has been shown that
the results are reproduced with excellent accuracy.
18.6 Literature
[7] EN 1993-1-1: Eurocode 3: Design of concrete structures, Part 1-1: General rules and rules for
buildings. CEN. 2005.
[9] Structural Eurocodes - Extracts from the structural Eurocodes for students of structural design.
BSI - British Standards Institution. 2007.
Overview
h d My
Np
zp
Ap
• Stress-strain curves for concrete and prestressing steel (Section 3.1.7, 3.3.6)
NP / A NP · zp / W2 M/ W2
+ + =
zp
Ap
NP / A NP · zp / W1 M/ W1
In rectangular cs, which are prestressed and loaded, stress conditions are developed, as shown in Fig.
19.2, where the different contributions of the loadings can be seen. The design stress-strain diagrams
for prestressing steel is presented in Fig. 19.3, as defined in DIN EN 1992-1-1:2004 [1] (Section 3.3.6).
A
σp
ƒpk
ƒp0,1k
ƒpd = ƒp0,1k / γs
B A Idealised
B Design
εp
Figure 19.3: Idealised and Design Stress-Strain Diagram for Prestressing Steel
d = 95.0 cm
L = 20.0 m
1
0.916 0.00
3572.3 0.057
0.917
3574.2
0.917 0.057
1.00
3576.1 0.057
0.918
3578.1
0.918 0.114
2.00
3580.2 0.057
0.919
3582.6
0.920 0.171
3.00
3585.6 0.055
0.920
3588.9
0.223
Material Properties
0.921 4.00
3592.4 0.050
0.922
3596.1
0.923 0.271
5.00
3600.0 0.045
0.924
3604.2
0.926 0.313
6.00
3608.5 0.038
0.927
3613.1
0.928 0.347
7.00
3618.0 0.031
0.929
3623.0
0.931 0.374
8.00
3628.3 0.022
0.932
3633.8
0.933 0.390
9.00
3639.5 0.011
0.935
10001 1001 10002 1002 10003 1003 10004 1004 10005 1005 10006 1006 10007 1007 10008 1008 10009 1009 10010
3645.4
0.396
2
0.937 10.00
L
3651.6 0.000
0.938
3657.8
0.940 0.390
11.00
Ap = 28.5 cm2
3663.8 -0.011
0.941
3680.5 0.944
0.944 0.347
3682.2 13.00
3678.8 -0.031
Table 19.1: (continued)
0.942
3673.9
0.941 0.313
14.00
3669.2 -0.038
DCE-EN17: Stress Calculation at a Rectangular Prestressed Concrete CS
0.938
3656.6
0.937 0.223
16.00
-0.050
10011 1010 10012 1011 10013 1012 10014 1013 10015 1014 10016 1015 10017 1016 10018 1017 10019 1018 10020
3634.3
0.000
3
20.00
-0.057
Loading (at = 10 m)
0.944
0.932
117
DCE-EN17: Stress Calculation at a Rectangular Prestressed Concrete CS
0 −4.82 −4.82
My [kNm] −185.91 −185.91
V 1 −4.51 −4.51
My [kNm] −156.11 −156.11
Material:
19 wires with area of 150 mm2 each, giving a total of Ap = 28.5 cm2
Cross-section:
Ac = 1.0 · 1.0 = 1 m2
The force applied to a tendon, i.e. the force at the active end during
tensioning, should not exceed the following value 5.10.2.1 (1)P: Prestressing force during
tensioning - Maximum stressing force
Pm = Ap · σp,m 5.10.2.1 (1)P: Eq. 5.41: Pm maximum
stressing force
where σp,m = min 0.80ƒpk ; 0.90ƒp0,1k
(NDP) 5.10.2.1 (1)P: σp,m maximum
stress applied to the tendon
Pm = Ap · 0.80 · ƒpk = 28.5 · 10−4 · 0.80 · 1770 = 4035.6 kN
Pm0 () = Ap · σp,m0 () 5.10.3 (2): Eq. 5.43: Pm0 initial pre-
stress force at time t = t0
1 The tools used in the design process are based on steel stress-strain diagrams, as
defined in [1] 3.3.6: Fig. 3.10, which can be seen in Fig 19.3.
2 The sections mentioned in the margins refer to DIN EN 1992-1-1:2004 (German Na-
Load Actions:
→ g1 = γ · A = 25 · 1 = 25 kNm
DIN EN 1990/NA [10]: (NDP) A.1.3.1 Safety factors at ultimate limit state
(4): Tab. NA.A.1.2 (B): Partial factors
for actions
(NDP) 2.4.2.2 (1): Partial factors for pre-
Actions (unfavourable) Safety factor at final state
stress • permanent γG = 1.35
• prestress γP = 1.00
Mg = g1 · L2 / 8 = 1250 kNm
The concrete stresses may be deter- Calculation of stresses σc,b at = 10.0 m middle of the span:
mined for each construction stage un-
der the total quasi-permanent combina- Position of the tendon: z = 0, 3901 m
tion σc {G + Pm0 + ψ2 · Q}
In this Benchmark no variable load Q is
• Case : Prestress at construction stage section 0 (P cs0)
examined
−σc
Mp
Np
zp
Pm0,=10
+σc
Np = −3653.0 kN
• Case : Prestress and self-weight at con. stage sect. 0 (P+G cs0)
−σc
Mg + Mp
Np
zp
Pm0,=10
+σc
−3653.0 −185.91
σc,b = + = −4.82 MP
0.9926 0.1633
Np My
σc,b = + W1,cs1 cross-section moduli for con-
Ade W1,cs1 truction stage 1
−3653.0 −1406.11
σc,b = + = −11.76 MP
1.013 0.172
• Case V: Prestress and self-weight at con. stage sect. 1 (P+G cs1)
−3653.0 −156.11
σc,b = + = −4.51 MP
1.013 0.172
19.5 Conclusion
This example shows the calculation of the stresses, developed in the concrete cross-section due to
prestress and bending. It has been shown that the results are reproduced with excellent accuracy.
19.6 Literature
[1] DIN EN 1992-1-1/NA: Eurocode 2: Design of concrete structures, Part 1-1/NA: General rules and
rules for buildings - German version EN 1992-1-1:2005 (D), Nationaler Anhang Deutschland -
Stand Februar 2010. CEN. 2010.
[2] F. Fingerloos, J. Hegger, and K. Zilch. DIN EN 1992-1-1 Bemessung und Konstruktion von
Stahlbeton- und Spannbetontragwerken - Teil 1-1: Allgemeine Bemessungsregeln und Regeln
für den Hochbau. BVPI, DBV, ISB, VBI. Ernst & Sohn, Beuth, 2012.
[4] Beispiele zur Bemessung nach Eurocode 2 - Band 1: Hochbau. Ernst & Sohn. Deutschen Beton-
und Bautechnik-Verein E.V. 2011.
[10] DIN EN 1990/NA: Eurocode: Basis of structural design, Nationaler Anhang Deutschland DIN EN
1990/NA:2010-12. CEN. 2010.
Overview
h My
Np
zp
Ap
• Time dependent losses of prestress for pre- and post-tensioning (Section 5.10.6)
The time dependant losses may be calculated by considering the following two reductions of stress [1]:
• due to the reduction of strain, caused by the deformation of concrete due to creep and shrinkage,
under the permanent loads
• the reduction of stress in the steel due to the relaxation under tension.
In this Benchmark the stress loss due to creep and shrinkage will be examined.
Ap = 28.5 cm2
Material:
Ecm = 34077 N/ mm2 3.1.2: Tab. 3.1: Ecm , ƒck and ƒcm for
C 35/ 45
ƒck = 35 N/ mm2
ƒcm = 43 N/ mm2
19 wires with area of 150 mm2 each, giving a total of Ap = 28.5 cm2
Cross-section:
Ac = 1.0 · 1.0 = 1 m2
Load Actions:
Mg = g1 · L2 / 8 = 1250 kNm
−σc
Mg + Mp
Np
zp
Pm0,=10
+σc
−3653.0 −185.91
σc,QP stress in concrete σc,QP = + = −4.82 MP
0.9926 0.1633
3.1.4 (6): Eq. 3.8: εcs total shrinkage εcs = εcd + εc
strain
3.1.4 (6): Eq. 3.9: εcd drying shrinkage εcd (t) = βds (t, ts) · kh · εcd,0
strain
The development of the drying shrinkage strain in time is strongly de-
pends on βds (t, ts) factor. SOFiSTiK accounts not only for the age at
start of drying ts but also for the influence of the age of the prestressing
t0 . Therefore, the calculation of factor βds reads:
3.1.4 (6): Eq. 3.10: βds βds = βds (t, ts ) − βds (t0 , ts )
(t − ts ) (t0 − ts )
3.1.4 (6): h0 the notional size (mm) of βds = Ç − Ç
the cs h0 = 2Ac / = 500 mm (t − ts ) + 0.04 · h30 (t0 − ts ) + 0.04 · h30
(1000028 − 0) (28 − 0)
βds = p − p
(1000028 − 0) + 0.04 · 5003 (28 − 0) + 0.04 · 5003
βds = 0.99955 − 0.05892 = 0.94063
εc (t) = βs (t) · εc (∞) 3.1.4 (6): Eq. 3.11: εc autogenous
shrinkage strain
εc (∞) = 2.5 (ƒck − 10) · 10−6 = 2.5 (35 − 10) · 10−6 3.1.4 (6): Eq. 3.12: εc (∞)
βs = βs (t) − βs (t0 ) 3.1.4 (6): Eq. 3.13: βs
p p p p
βs = 1 − e−0.2· t − 1 − e−0.2· t0 = e−0.2· t0 − e−0.2· t
βs = 0.347
εcs = 1.668 · 10−4 + 2.169 · 10−5 = −18.85 · 10−5 negative sign to declare losses
16.8 p
β(ƒcm ) = p = 16.8/ 43 = 2.562 Annex B.1 (1): Eq. B.4: β(ƒcm ) factor
ƒcm for effect of concrete strength on creep
0.7
35
α1 = ≤ 1 = 0.8658 Annex B.1 (1): Eq. B.8c: α1 , α2 , α3
ƒcm coefficients to consider influence of con-
crete strength
0.2
35
α2 = ≤ 1 = 0.9597
ƒcm
0.5
35
α3 = ≤ 1 = 0.9022
ƒcm
1 − 80/ 100
ϕRH = 1+ p
3
· 0.8658 · 0.9597 = 1.1691
0.1 · 500
1
β(t0 ) = Annex B.1 (1): Eq. B.5: β(t0 ) factor
0.1 + t00.20 for effect of concrete age at loading on
creep
!α
9
t0 = t0,T · 1.2
+1 ≥ 0.5 Annex B.1 (2): Eq. B.9: t0,T tempera-
2 + t0,T ture adjusted age of concrete at loading
adjusted according to expression B.10
⇒ βc (t, t0 ) = 0.9996
Annex B.1 (3): The values of ϕ(t, t0 ) ϕ(t, t0 ) = 1.463 · 0.9996/ 1.05 = 1.393
given above should be associated with
the tangent modulus Ec According to EN, the creep value is related to the tangent Young’s
3.1.4 (2): The values of the creep coef-
modulus Ec , where Ec being defined as 1.05 · Ecm . To account for
ficient, ϕ(t, t0 ) is related to Ec , the tan-
gent modulus, which may be taken as this, SOFiSTiK adopts this scaling for the computed creep coefficient
1.05 · Ecm (in SOFiSTiK, all computations are consistently based on Ecm ).
Ep
εcs · Ep + 0.8Δσpr + ϕ(t, t0 ) · σc,QP
Ecm
ΔPc+s+r = Ap ·Δσp,c+s+r = Ap
Ep Ap Ac
1+ 1+ 2
zcp [1 + 0.8ϕ(t, t0 )]
5.10.6 (2): Eq. 5.46: ΔPc+s+r time de- Ecm Ac c
pendent losses of prestress
In this example only the losses due to creep and shrinkage are taken
into account, the reduction of stress due to relaxation (Δσpr ) is ignored.
20.5 Conclusion
This example shows the calculation of the time dependent losses due to creep and shrinkage. It has
been shown that the results are in very good agreement with the reference solution.
20.6 Literature
[1] DIN EN 1992-1-1/NA: Eurocode 2: Design of concrete structures, Part 1-1/NA: General rules and
rules for buildings - German version EN 1992-1-1:2005 (D), Nationaler Anhang Deutschland -
Stand Februar 2010. CEN. 2010.
[2] F. Fingerloos, J. Hegger, and K. Zilch. DIN EN 1992-1-1 Bemessung und Konstruktion von
Stahlbeton- und Spannbetontragwerken - Teil 1-1: Allgemeine Bemessungsregeln und Regeln
für den Hochbau. BVPI, DBV, ISB, VBI. Ernst & Sohn, Beuth, 2012.
[4] Beispiele zur Bemessung nach Eurocode 2 - Band 1: Hochbau. Ernst & Sohn. Deutschen Beton-
und Bautechnik-Verein E.V. 2011.
Overview
b
t t
t As2
C
SC
h
t As1
with respect to DIN EN 1992-1-1/NA [11] [2]. The results are given in Table 21.2.
t = 400.0 cm LC 912:
L = 20.0 m Vz = 660 kN, My = 4950 kNm, Mt = −50.20 kNm
ΔσRsk (N∗ )
[MP] 152.17 152.2
γs,ƒ t
STEP 1: Material
STEP 2: Cross-section
Minimun reinforcements:
Permanent: Loadcase 1
Variable: Loadcase 2, 3
• LC 911 G
MAXP + MY : 1.00 * G
• LC 912 G+2
MAXF + MY : 1.00 * G + ψ1 * LC 2
• LC 913 G+3
MAXF + MY : 1.00 * G + ψ1 * LC 3
1 The tools used in the design process are based on steel stress-strain diagrams, as
otherwise specified.
τT : shear stresses resulting from tor- τT = −1/ WT · Mt = −0.34484 · 0.189 = 0.065 · 10−3 MP
sion
τ = τQ + τT = 427.770·10−3 + 0.065·10−3 = 427.835·10−3 MP
τ
θ: angle of compression struts σ =
α: angle of shear reinforcement cot θ + cot α
α = 90o ⇒ sin α = 1.0, cot α = 0.0
A rather nasty problem is the evaluation of the shear. The DIN design
6.8.2(3): In the design of shear rein- code allows a simple solution based on a corrected value for the incli-
forcement the inclination of the com- nation of the compressive struts:
pressive struts θƒ t may be calculated
by Eq. 6.65 p
tn θƒ t = tn θ
6.8.2(3): Eq. 6.65: tn θƒ t
1740.048 · 10−3
σ = = 994.313 · 10−3 MP
1.75 + 0.0
ƒQ · τQ ƒT · τT
σs = + σs : steel stresses
(cot θƒ t + cot α) · sin α cot θƒ t
B0 · ƒr Bb · ƒr
ƒT = , and ƒQ = As/ ct = As / 2 = 4.11 cm2 / m
As/ ct As/ ct
B0
Accordingly, we calculate the stresses for the rest of the loadcases. For
each loadcase the stresses are calculated for two cases, for τT and for
−τT , in order to determine the most unfavorable case. The results are
presented in Table 21.4
From the table, the minimum and maximum value of the steel stress is
determined:
As the exact fatigue stress check is not available, the simplified methods
according to DIN EN 1992-1-1/NA (Sect. 6.8, Fatigue) are selected via
λ : Coeff. equiv. stress range shear the coefficients λs , λt , λ , λc .
links, here input as 1.0
The admissible sways of the damage equivalent stress range for the
6.8.2 (2): Eq. 6.64: η factor for effect of shear links are obtained, as follows:
different bond behaviour
As + Ap
η= , since Ap = 0 Δσs,eq (N∗ ) = λ · (σs,m − σs,mn ) = 1.0 · (391.77 − 314.79)
As + Ap ξ · ϕs / ϕp
Æ
(NDP) 6.8.7 (1): Eq. 6.76: ƒcd,ƒ t The design fatigue strength of concrete is determined by:
k1 = 1.0
3.1.2 (6): Eq. 3.2: βcc (t0 ) ƒck
p
βcc (t0 ) = es·(1− 28/ t = 1.0 ƒcd,ƒ t = k1 · βcc (t0 ) · ƒcd · 1 −
250
35
ƒcd,ƒ t = 1.0 · 1.0 · 19.83 · 1 − = 17.0567 MP
250
In the case of the compression struts of members subjected to shear,
the concrete strength ƒcd,ƒ t should be reduced by the strength reduc-
(NDP) 6.8.7 (3): Verification of concrete tion factor ν1 according to 6.2.3(3).
under compression or shear
(NDP) 6.2.3 (3): ν1 , ν2 ν2 = 1 for ≤ C50/ 60
ν1 = 0.75 · ν2 = 0.75
v
σcd,m,eq u σcd,mn,eq
+ 0.43 · 1− ≤1
t
ƒcd,ƒ t,red σcd,m,eq
v !
u σcd,mn,eq
σcd,m,eq ≤ ƒcd,ƒ t,red · 1.0 − 0.43 · 1−
t
σcd,m,eq
v
t 0.9933
σcd,m,eq ≤ 12.7925 · 1.0 − 0.43 · 1−
1.2362
Accordingly the above verification is done for the minimum and maxi-
mum nonlinear stresses of concrete, as calculated from AQB Fig. 21.3,
at the defined 0 TOP0 point of the cross-section.
v !
u σcd,mn,eq
σcd,m,eq ≤ ƒcd,ƒ t · 1.0 − 0.43 · 1 −
t
σcd,m,eq
v
t 5.69
σcd,m,eq ≤ 17.0567 · 1.0 − 0.43 · 1−
6.60
Stress limitation
-5.69
-4.30
-2.90
-1.47 -5.69
-0.02
-9.37
-9.37
-9.37
-9.37 -5.69
-9.37 -4.30
-9.37
-2.90
-1.47
-0.02
411.74
411.74
411.74
411.74
411.74
411.74
-6.60
-5.00
-3.37
-1.71 -6.60
-0.02
-10.97
-10.97
-10.97 -6.60
-10.97
-10.97 -5.00
-10.97
-3.37
-1.71
-0.02
479.20
479.20
479.20
479.20
479.20
479.20
21.5 Conclusion
This example shows the verification of a reinforced concrete beam to fatigue. It has been shown that
AQB follows the fatigue verification procedure, as proposed in DIN EN 1992-1-1/NA [11] [2]. The in-
significant deviation arises from the fact that the benchmark (reference) results have been calculated by
using the BEM analysis. By introducing the FEM analysis, AQUA calculates now the 1/ WVz , 1/ WVy
and 1/ WT values more accurate.
21.6 Literature
[2] F. Fingerloos, J. Hegger, and K. Zilch. DIN EN 1992-1-1 Bemessung und Konstruktion von
Stahlbeton- und Spannbetontragwerken - Teil 1-1: Allgemeine Bemessungsregeln und Regeln
für den Hochbau. BVPI, DBV, ISB, VBI. Ernst & Sohn, Beuth, 2012.
[11] DIN EN 1992-1-1/NA:2013-04: Eurocode 2: Design of concrete structures, Part 1-1/NA: Gen-
eral rules and rules for buildings - Nationaler Anhang Deutschland, Ersatz für DIN EN 1992-1-
1/NA:2011-01 und DIN EN 1992-1-1/NA Berichtigung 1:2012-06. DIN. 2013.
Overview
First it will be shown that a loading below the ultimate capacity will be covered by selected software even
if deviations in the applied formulas will affect the expected results slightly. Then it will be demonstrated
that an advanced analysis using additional features within the provisions of the design code allows a
better approach to the ultimate capacity.
b
tƒ
r
t
Vz
My
h y
h = 524 mm N = −5000 kN
tƒ = 40 mm, t = 21 mm
r = 27 mm
Material:
ƒy = 235 N/ mm2
Cross-sectional properties:
AV 99.5 A − 2 · b · tƒ
= = = 0.289 6.2.9.1 (5): =
A 344.3 A
p
129.498 · 235/ 3
Vp,Rd,z = = 1757 kN 6.2.6 (3)a: The shear area AVz may be
1.00 taken as follows for rolled I - sections
with load parallel to the web
Wp,y · ƒy
Mc,Rd = = 1667 kNm 6.2.5 (2), Eq. 6.13: The design resis-
γM0 tance for bending, for Class 1 cross-
section
Where the shear and axial force are present allowance should be made
for the effect of both shear force and axial force on resistance moment. 6.2.10 (1): Bending, shear and axial
force
Provided that the design value of the shear force VEd does not exceed
the 50% of the design plastic shear resistance Vp,Rd no reduction of
the resistances for bending and axial force need to be made.
VEd
= 0.797 > 0.5
Vp,Rd,z
→ shear resistance limit exceeded 6.2.10 (3): Bending shear and axial
force
Provided that the design value of the shear force VEd exceeds the 50%
of the design plastic shear resistance Vp,Rd,z , the design resistance
of the cross-section to combinations of moment and axial force should
be calculated using a reduced yield strength (1 − ρ) · ƒy for the shear
area. Therefore all following formulas will be modified for the shear force
allowance.
1 The sections mentioned in the margins refer to EN 1993-1-1:2005 [7] unless other-
wise specified.
2
2 · VEd
6.2.10 (3) the design resistance should ρ= −1 = 0.352
be calculated using a reduced yield Vp,Rd,z
strength (1 − ρ) · ƒy for the shear area
ρ · A2
Wp,y − · ƒy
4 · t
6.2.8 (5): Eq. 6.30: The reduced de- My,V,Rd = = 1581 kNm
sign plastic resistance moment allowing γM0
for shear force may alternatively be ob-
tained for I-sections with equal flanges My
and bending about major axis = 0.277
My,V,Rd
My
= 0.270 with no reduction
Mc,Rd
A · ƒy
6.2.4 (2), Eq. 6.10: The design resis- Np,Rd = = 8091 kN
tance for compression γM0
NV,Rd = 7018.6 kN
NEd
n = = 0.712
NV,Rd
NEd
= 0.618 with no reduction
Np,Rd
Where an axial force is present, allowance should be made for its effect
6.2.9.1 (1): Bending and axial force - on the plastic moment resistance.
Class 1 and 2 cross-sections
0.25 · NV,Rd
6.2.9.1 (4): For doubly symmetrical NEd > mn 0.5 · h · t · ƒy · (1 − ρ)
I-sections, allowance need not to be
made for the effect of the axial force on γM0
the plastic resistance moment about the (
y-y axis when the criteria are satisfied. 1754.7 kN
NEd = 5000 > mn
709.5 kN
NEd
n= = 0.712
NV,Rd
∗ = mn A 344.4
0.5
∗ = 0.289
1 − 0.712
MN,y,Rd = 1581 · = 531 kNm
1 − 0.5 · 0.289
My,Ed 450
= = 0.846
MNV,y,Rd 531
0.352 · 129.52
My,V,Rd = 7094 − · ƒy Using A = A = AVz insead of
4 · 2.1 three different values: 93.2 6= 99.5 6=
129.5
My,V,Rd = 1501.8 kNm Hit: The reduction by the plastic bend-
ing resistance of the web may be ob-
1 − 0.712 tained with a better precision.
MN,y,Rd = · 1501 = 532 kNm
1 − 0.5 · 0.376
MEd 450
= = 0.846
MNV,y,Rd 532
22.5 Conclusion
Interaction formulas are not expected to be true for all cases [12]. Thus deviations are unavoidable.
The first reason of the deviation is that the interaction formulas are not linear. Thus an utilisation factor
of 0.5 does not mean that the ultimate load is twice the current loading. If the utilisation factor of the
normal force or the shear force become 1.0 or larger then the formulas are not applicable any more. So
the utilisation formula should be rewritten from the original form of:
M 1 − 0.5 ·
· ≤ 1.00 (22.1)
My,V,Rd N
1−
NV,Rd
M N
· (1 − 0.5 · Vz ) + ≤ 1.00 (22.2)
My,V,Rd NV,rd
These are two different curves intersecting for the critical value of 1.0, but the singularity for the utilisation
N/ NV,Rd = 1 is avoided. The second formula yields a value of 0.956 with the EN reference data.
The simplified interaction according to equation 6.2 of EN 1993-1-1 yields an utilisation value of 1.04.
Program RUBSTAHL QST-I [12] yields an utilisation of 1.00.
The second reason for the deviations is that the software follows a more general principle and uses only
a single value for the shear capacity as indicated in the last row of table 22.2. Thus the software has to
select either the web only or the dog-bone-shape including parts of the flanges, as can be visualized in
Fig. 22.2.
AV AVz
These deviations become larger, if we try to find the true capacity of the section. It has to be noted,
that most interaction formulas for shear and normal stress are working with partial areas and are not
accounting for the true interaction of stresses. There are different possibilities in the literature and there
are some numerical procedures by Katz [13] or Osterrieder [14] allowing for more detailed evaluations.
The solution with the optimization process from Osterrieder [14] yields a true utilisation factor of 0.93 for
the given forces. Thus the true ultimate plastic force combination is: N = 5380 kN, Vz = 1505 kN,
My = 483.9 kNm. The software provides for these forces with the modified interaction formulas an
utilisation of 1.08. If one follows strictly the design code, much higher values of about 2.7 are obtained,
showing the strong non proportionality of those formulas.
The solution according to Katz [13] uses a non-linear analysis, thus we have finite strains and the hard-
ening effect of the steel strain law is used, as indicated in EN 1993-1-5:2006 Appendix C.6 [15]. Then it
is no problem to obtain even higher ultimate forces: N = 5800 kN, Vz = 1624 kN, My = 522 kNm,
giving a utilisation factor of 0.86 for the given loading. The normal and shear stress distribution in Figure
22.3 is obtained with a maximum compressive strain of −1.61%◦ .
-22
-19
5.7
1.3
19
22
0
7.3
6.6
5
7
.89
52
17
.98
0.3
61
6
-53
.53
1 8
0.2 4.9
12 21
4
3.8
8 5.5
13 23
0
1.5
3 7.3
14 24
AV
-94
6
6.4
1 4.6
.01
14 25
94
16
.24
23
13
8.5
3.9
2.4
6
5
Non linear Axial Stress, f = 1.080 Non linear Shear Stress, f = 1.080 Non linear von Mises stress, f = 1.080
22.6 Literature
[7] EN 1993-1-1: Eurocode 3: Design of concrete structures, Part 1-1: General rules and rules for
buildings. CEN. 2005.
[12] R. Kindmann and U. Krüger. Stahlbau, Teil 1: Grundlagen. Ernst & Sohn, 2015.
[13] C. Katz. “Fließzonentheorie mit Interaktion aller Stabschnittgrößen bei Stahltragwerken”. In:
Stahlbau (1997).
[14] P. Osterrieder. “Plastic Bending and Torsion of Open Thin-Walled Steel Members”. In: Proceedings
of EUROSTEEL 2005, 4th European Conference on Steel and Composite Structures (Maastricht,
2005).
[15] DIN EN 1993-1-5: Eurocode 3: Design of steel structures, Part 1-5: Plated structural elements -
German version EN 1993-1-1:2006 + AC:2009. CEN. 2010.
Overview
beƒ ƒ
hƒ
VEd
h d
As1 As1
d1
The time dependant displacements are calculated by multiplying the length of the beam with the creep
(εcc ) and shrinkage (εcs ) strain:
• the creep deformation of concrete is calculated according to (EN 1992-1-1, 3.1.4, Eq. 3.6)
• the total shrinkage strain is calculated according to (EN 1992-1-1, 3.1.4, Eq. 3.8)
tendon geometry is simplified as much as possible and modelled as a horizontal force, therefore tendons
are not subject of this benchmark. The beam consists of a T-Beam cs and is loaded with a horizontal
prestressing force from time t1 = 100 days to time t2 = 300 days. The self-weight is neglected. A
calculation of the creep and shrinkage is performed in the middle of the span with respect to EN 1992-
1-1:2004 [6]. The calculation steps are presented below and the results are given in Table 23.2 for the
calculation with CSM. For calculating the real creep and shrinkage (RCRE) an equivalent loading is
used, see Fig. 23.2 and Fig. 23.3.
Real loading
shrinkage
loading
creep
t [days]
ts t0 t1 = 100 t2 = 300 t∞ = 30 yr.
Real loading
Figure 23.3: Equivalent loading and displacement for real creep and shrinkage (RCRE)
3. Calculating the displacements (creep and shrinkage) at time before the loading is inactive (t2 ≈
300 and t2 < 300 days).
4. Calculating the displacements at time when the loading is inactive (t2 ≈ 300 and t2 > 300 days).
RH = 80 hƒ = 40 cm
b = 40 cm
L = 20.0 m
Material:
3.1.2: Tab. 3.1: Ecm , ƒck and ƒcm for Ecm = 34077 N/ mm2
C 35/ 45
ƒck = 35 N/ mm2
ƒcm = 43 N/ mm2
Load Actions:
Self weight per length is neglected: γ = 0 kN/ m (to simplify the exam-
ple as much as possible)
NEd = −900 kN
A = 280 · 40 + 60 · 80 = 16000 cm2
• Calculating creep:
• Calculating shrinkage:
3.1.4 (6): Eq. 3.8: εcs total shrinkage εcs = εcd + εc
strain
3.1.4 (6): Eq. 3.9: εcd drying shrinkage εcd (t) = βds (t, ts ) · kh · εcd,0
strain
1 The tools used in the design process are based on steel stress-strain diagrams, as
otherwise specified.
βds = βds (t, ts ) − βds (t0 , ts ) 3.1.4 (6): Eq. 3.10: βds
(t − ts ) (t0 − ts )
βds = Ç − Ç 3.1.4 (6): h0 the notional size (mm) of
(t − ts ) + 0.04 · h30 (t0 − ts ) + 0.04 · h30 the cs h0 = 2Ac / = 500 mm
(100 − 3) (7 − 3)
βds = p − p
(100 − 3) + 0.04 · 4003 (7 − 3) + 0.04 · 4003
βds = 0.232 − 0.01235 = 0.22026
43
εcd,0 = 0.85 (220 + 110 · 4) · exp −0.12 · · 10−6 · 0.7564
10
Annex B.2 (1): αds1 , αds1 coefficients
depending on type of cement.
εcd,0 = 2.533 · 10−4 For class N αds1 = 4, αds2 = 0.12
Drying shrinkage:
εc (t) = βs (t) · εc (∞) 3.1.4 (6): Eq. 3.11: εc autogenous
shrinkage strain
εc (∞) = 2.5 (ƒck − 10) · 10−6 = 2.5 · (35 − 10) · 10−6 3.1.4 (6): Eq. 3.12: εc (∞)
βs = βs (t) − βs (t0 ) 3.1.4 (6): Eq. 3.13: βs
p p p p
βs = 1 − e−0.2· t − 1 − e−0.2· t0 = e−0.2· t0 − e−0.2· t
βs = 0.4537
Autogenous shrinkage:
Total shrinkage:
ε absolute shrinkage strain εcs = −2.84 · 10−5 + (−4.04) · 10−5 = −6.881 · 10−5
negative sign to declare losses
Calculating displacement:
Δ1,cs = εcs · L/ 2
Δ1,cs = −0.6881 mm
• Calculating shrinkage:
3.1.4 (6): Eq. 3.8: εcs total shrinkage εcs = εcd + εc
strain
3.1.4 (6): Eq. 3.9: εcd drying shrinkage εcd (t) = βds (t, ts) · kh · εcd,0
strain
The development of the drying shrinkage strain in time is strongly de-
pends on βds (t, ts) factor. SOFiSTiK accounts not only for the age at
start of drying ts but also for the influence of the age of the prestressing
βds = βds (t, ts ) − βds (t0 , ts ) 3.1.4 (6): Eq. 3.10: βds
(t − ts ) (t0 − ts )
βds = Ç − Ç 3.1.4 (6): h0 the notional size (mm) of
(t − ts ) + 0.04 · h30 (t0 − ts ) + 0.04 · h30 the cs h0 = 2Ac / = 500 mm
(300 − 3) (100 − 3)
βds = p − p
(300 − 3) + 0.04 · 4003 (100 − 3) + 0.04 · 4003
βds = 0.2487
43
εcd,0 = 0.85 (220 + 110 · 4) · exp −0.12 · · 10−6 · 0.7564
10
Annex B.2 (1): αds1 , αds1 coefficients
depending on type of cement.
εcd,0 = 2.533 · 10−4 For class N αds1 = 4, αds2 = 0.12
Drying shrinkage:
εc (t) = βs (t) · εc (∞) 3.1.4 (6): Eq. 3.11: εc autogenous
shrinkage strain
εc (∞) = 2.5 · (ƒck − 10) · 10−6 = 2.5 (35 − 10) · 10−6 3.1.4 (6): Eq. 3.12: εc (∞)
βs = βs (t) − βs (t0 ) 3.1.4 (6): Eq. 3.13: βs
p p p p
βs = 1 − e−0.2· t − 1 − e−0.2· t0 = e−0.2· t0 − e−0.2· t
βs = 0.104
Autogenous shrinkage:
Total shrinkage:
Δ3,cs = εcs · L/ 2
Δ3,cs = −0.5218 mm
• Calculating creep:
⇒ βc (t, t0 ) = 0.4916
ϕ(t, t0 ) = ϕ0 · βc (t, t0 )
ϕeƒ ƒ (t, t0 ) = 0.57/ 1.05 = 0.5445 Annex B.1 (3): The values of ϕ(t, t0 )
given above should be associated with
According to EN, the creep value is related to the tangent Young’s the tangent modulus Ec
3.1.4 (2): The values of the creep coef-
modulus Ec , where Ec being defined as 1.05 · Ecm . To account for
ficient, ϕ(t, t0 ) is related to Ec , the tan-
this, SOFiSTiK adopts this scaling for the computed creep coefficient gent modulus, which may be taken as
(in SOFiSTiK, all computations are consistently based on Ecm ). 1.05 · Ecm
• Calculating shrinkage:
3.1.4 (6): Eq. 3.8: εcs total shrinkage εcs = εcd + εc
strain
3.1.4 (6): Eq. 3.9: εcd drying shrinkage εcd (t) = βds (t, ts) · kh · εcd,0
strain
The development of the drying shrinkage strain in time is strongly de-
pends on βds (t, ts) factor. SOFiSTiK accounts not only for the age at
start of drying ts but also for the influence of the age of the prestressing
t0 . Therefore, the calculation of factor βds reads:
3.1.4 (6): Eq. 3.10: βds βds = βds (t, ts ) − βds (t0 , ts )
(t − ts ) (t0 − ts )
3.1.4 (6): h0 the notional size (mm) of βds = Ç − Ç
the cs h0 = 2Ac / = 500 mm (t − ts ) + 0.04 · h30 (t0 − ts ) + 0.04 · h30
(11250 − 3) (300 − 3)
βds = p − p
(11250 − 3) + 0.04 · 4003 (300 − 3) + 0.04 · 4003
βds = 0.49097
43
εcd,0 = 0.85 (220 + 110 · 4) · exp −0.12 · · 10−6 · 0.7564
10
Annex B.2 (1): αds1 , αds1 coefficients
depending on type of cement.
For class N αds1 = 4, αds2 = 0.12 εcd,0 = 2.533 · 10−4
Drying shrinkage:
εc (t) = βs (t) · εc (∞) 3.1.4 (6): Eq. 3.11: εc autogenous
shrinkage strain
εc (∞) = 2.5 · (ƒck − 10) · 10−6 = 2.5 (35 − 10) · 10−6 3.1.4 (6): Eq. 3.12: εc (∞)
βs = βs (t) − βs (t0 ) 3.1.4 (6): Eq. 3.13: βs
p p p p
βs = 1 − e−0.2· t − 1 − e−0.2· t0 = e−0.2· t0 − e−0.2· t
βs = 0.03130
Autogenous shrinkage:
Total shrinkage:
Δ5,cs = εcs · L/ 2
Δ5,cs = −0.9212 mm
• Calculating creep:
16.8 p
Annex B.1 (1): Eq. B.4: β(ƒcm ) factor β(ƒcm ) = p = 16.8/ 43 = 2.562
for effect of concrete strength on creep ƒcm
0.7
35
Annex B.1 (1): Eq. B.8c: α1 , α2 , α3 α1 = = 0.8658 ≤ 1
coefficients to consider influence of con- ƒcm
crete strength
0.2
35
α2 = = 0.9597 ≤ 1
ƒcm
0.5
35
α3 = = 0.9022 ≤ 1
ƒcm
1 − 80/ 100
ϕRH = 1+ p
3
· 0.8658 · 0.9597 = 1.1852
0.1 · 400
1
Annex B.1 (1): Eq. B.5: β(t0 ) factor β(t0 ) =
for effect of concrete age at loading on 0.1 + t00.20
creep
!α
9
Annex B.1 (2): Eq. B.9: t0,T tempera- t0 = t0,T · 1.2
+1 ≥ 0.5
ture adjusted age of concrete at loading 2 + t0,T
adjusted according to expression B.10
e−(4000/ [273+T(Δt )]−13.65) · Δt
Pn
Annex B.1 (3): Eq. B.10: tT temper- tT = =1
ature adjusted concrete age which re-
places t in the corresponding equations
t0,T = 300 · e−(4000/ [273+20]−13.65) = 300 · 1.0 = 300.0
0
9
Annex B.1 (2): Eq. B.9: α a power ⇒ t0 = 300 · +1 = 300
which depends on type of cement 2 + 3001.2
For class N α = 0
1
β(t0 ) = = 0.3097
0.1 + 3000.20
With MEXT NO 1 EIGE VAL 0.3773 the exponent 0.3 can be modified to
0.3773.
0.3773
(t − t0 )
Annex B.1 (1): Eq. B.7: βc (t, t0 ) co- βc (t, t0 ) =
efficient to describe the development of (βH + t − t0 )
creep with time after loading
βH = 1.5 · 1 + (0.012 · RH)18 · h0 + 250 · α3 ≤ 1500 · α3
Annex B.1 (1): Eq. B.8: βH coefficient
depending on relative humidity and no-
βH = 1.5 · 1 + (0.012 · 80)18 · 400 + 250 · 0.9022
tional member size
⇒ βc (t, t0 ) = 0.9641
ϕ(t, t0 ) = ϕ0 · βc (t, t0 )
ϕeƒ ƒ (t, t0 ) = 0.91/ 1.05 = 0.8637 Annex B.1 (3): The values of ϕ(t, t0 )
given above should be associated with
According to EN, the creep value is related to the tangent Young’s the tangent modulus Ec
3.1.4 (2): The values of the creep coef-
modulus Ec , where Ec being defined as 1.05 · Ecm . To account for
ficient, ϕ(t, t0 ) is related to Ec , the tan-
this, SOFiSTiK adopts this scaling for the computed creep coefficient gent modulus, which may be taken as
(in SOFiSTiK, all computations are consistently based on Ecm ). 1.05 · Ecm
σc
εcc (t, t0 ) = ϕ(t, t0 ) ·
Ecs
−0.5625
εcc (t, t0 ) = 0.91 · = −1.4087 · 10−5
36334.29
Δ
ε= → Δ5,cc = εcc · L/ 2
Δ5,cc = −1.4087 · 10−5 · 10000 mm = −0.1408 mm
Δ4010 = 0 mm
Δ4015 = Δ1,cs
Δ4015 = −0.688 mm
Δ4020 = −0.8431 mm
Δ4025 = −1.45314 mm
Δ4030 = −1.29814 mm
Δ4035 ≈ −2.08 mm
23.5 Conclusion
This example shows the calculation of the time dependent displacements due to creep and shrinkage.
It has been shown that the results are in very good agreement with the reference solution.
23.6 Literature
[2] F. Fingerloos, J. Hegger, and K. Zilch. DIN EN 1992-1-1 Bemessung und Konstruktion von
Stahlbeton- und Spannbetontragwerken - Teil 1-1: Allgemeine Bemessungsregeln und Regeln
für den Hochbau. BVPI, DBV, ISB, VBI. Ernst & Sohn, Beuth, 2012.
[6] EN 1992-1-1: Eurocode 2: Design of concrete structures, Part 1-1: General rules and rules for
buildings. CEN. 2004.
[16] fib Model Code 2010. fib Model Code for Concrete Structures 2010. International Federation for
Structural Concrete (fib). 2010.
Overview
h My
Np
zp
Ap
• Annex D: Detailed calculation method for prestressing steel relaxation losses (Section D.1)
In this Benchmark the stress loss due to relaxation will be examined, creep and shrinkage losses are
neglected and disabled.
in benchmark 17. The beam consists of a rectangular cs and is prestressed and loaded with its own
weight. A calculation of the relaxation stress losses is performed in the middle of the span with respect
to EN 1992-1-1:2004 [6]. The calculation steps are presented below and the results are given in Table
24.2 for the calculation with AQB.
RH = 80 L = 20.0 m
Ap = 28.5 cm2
ΔPpr
[%] 1.12% 1.12%
P0
• Material:
Ecm = 34077 N/ mm2 3.1.2: Tab. 3.1: Ecm , ƒck and ƒcm for
C 35/ 45
ƒck = 35 N/ mm2
ƒcm = 43 N/ mm2
19 wires with area of 150 mm2 each, giving a total of Ap = 28.5 cm2
• Cross-section:
Ac = 1.0 · 1.0 = 1 m2
The force applied to a tendon, i.e. the force at the active end during
tensioning, should not exceed the following value
Pm = Ap · σp,m
1 The tools used in the design process are based on steel stress-strain diagrams, as
5.10.3 (2): Eq. 5.43: Pm0 initial pre- Pm0 () = Ap · σp,m0 ()
stress force at time t = t0
where σp,m0 () = min 0.75ƒpk ; 0.85ƒp0,1k
(NDP) 5.10.3 (2): σp,m0 () stress in the
tendon immediately after tensioning or
transfer Pm0 = Ap · 0.75 · ƒpk = 28.5 · 10−4 · 0.75 · 1770 = 3783.4 kN
P0 = 3653 kN
σp0 1281.755
μ= = = 0.724155 ≥ μmn = 0.55
ƒpk 1770
The relaxation loss may be obtained from the manufactures test certifi-
cates or defined as percentage ratio of the variation of the prestressing
stress over the initial prestressing stress, should be determined by ap-
EC2-1-1, 3.3.2 (7): Eq. 3.29 plying next expression:
0.75·(1−μ)
Δσpr t
= 0.66 · ρ1000 · e9.1·μ · · 10−5
σp 1000
0.75·(1−μ)
t
t = 1000 h → = 1.00
1000
The values for ρ1000 can be either assumed equal to 8% for Class 1,
EC2-1-1, 3.3.2 (6) 2.5% for Class 2, and 4% for Class 3, or taken from the certificate.
Δσpr
= 0.66 · ρ1000 · e9.1·0.724155 · 10−5 · 100 = 1.200616 %
σp
AQB reduces the initial stress according to DIN 1045-1: DIN 1045-1 8.7.3 (6)
Ep
εcs · Ep + Δσpr + · ϕ(t, t0 ) · σc,QP
Ecm
ΔPc+s+r = Ap DIN 1045-1 8.7.3 (6): ΔPc+s+r time de-
Ep Ap Ac
pendent losses of prestress
1+ · 1+ · zcp [1 + 0.8 · ϕ(t, t0 )]
2
Ecm Ac c
Creep and shrinkage is not taken into account therefore we have:
εcs · Ep = 0
Ep
· ϕ(t, t0 ) · σc,QP = 0
Ecm
1 + 0.8 · ϕ(t, t0 ) = 1
15.388
Δσp,c+s+r =
195000 28.5 · 10−4 0.9926
1+ · · 1+ · 0.3901 2
34077 0.9926 0.08214
15.388
Δσp,c+s+r = = 14.70 MP
1.046644
Now we have:
σp0 = 1281.755 − 0.3 · 14.70 = 1277.345 MP DIN 1045-1 8.7.3 (6)
Now with the reduced stress the μ value is calculated again (iteration):
σp0 1277.345
μ= = = 0.72166 ≥ μmn = 0.55
ƒpk 1770
Δσpr
= 0.66 · ρ1000 · e9.1·0.72166 · 10−5 · 100 = 1.1736 %
σp
Δσpr,c = (εk0 + εky · zcp − εkz · ycp ) · Ep εk0 , εky and εkz - tendon strain
ΔPpr
ΔPpr [%] =
P0
40.9545
ΔPpr [%] = = 0.01120 · 100% = 1.12 %
3653
24.5 Conclusion
This example shows the calculation of the time dependent losses due to relaxation. It has been shown
that the reference solution and the AQB solution are in very good agreement.
24.6 Literature
[6] EN 1992-1-1: Eurocode 2: Design of concrete structures, Part 1-1: General rules and rules for
buildings. CEN. 2004.
Overview
h My
Np
zp
Ap
• Annex D: Detailed calculation method for prestressing steel relaxation losses (Section D.1)
In this Benchmark the stress loss due to relaxation will be examined, creep and shrinkage losses are
neglected and disabled.
defined in Table 25.1. Further information about the tendon geometry and prestressing can be found
in benchmark 17. The beam consists of a rectangular cs and is prestressed and loaded with its own
weight. A calculation of the relaxation stress losses is performed in the middle of the span with respect
to DIN EN 1992-1-1:2004 [1] [2]. The calculation steps are presented below and the results are given in
Table 25.2 for the calculation with AQB.
RH = 80 L = 20.0 m
Ap = 28.5 cm2
ΔPpr
[%] 2.274% 2.274%
P0
• Material:
Ecm = 34077 N/ mm2 3.1.2: Tab. 3.1: Ecm , ƒck and ƒcm for
C 35/ 45
ƒck = 35 N/ mm2
ƒcm = 43 N/ mm2
19 wires with area of 150 mm2 each, giving a total of Ap = 28.5 cm2
• Cross-section:
Ac = 1.0 · 1.0 = 1 m2
The force applied to a tendon, i.e. the force at the active end during
tensioning, should not exceed the following value
Pm = Ap · σp,m
1 The tools used in the design process are based on steel stress-strain diagrams, as
5.10.3 (2): Eq. 5.43: Pm0 initial pre- Pm0 () = Ap · σp,m0 ()
stress force at time t = t0
where σp,m0 () = min 0.75ƒpk ; 0.85ƒp0,1k
(NDP) 5.10.3 (2): σp,m0 () stress in the
tendon immediately after tensioning or
transfer Pm0 = Ap · 0.75 · ƒpk = 28.5 · 10−4 · 0.75 · 1770 = 3783.4 kN
P0 = 3653 kN
σp0 1281.755
μ= = = 0.724155 ≥ μmn = 0.55
ƒpk 1770
The formulas in section DIN EN 1992-1-1, 3.3.2 (7) are not used
for the calculation, according to DiBt [17] for DIN EN 1992-1-1 [6] the
relaxation values are obtained from Table 25.3.
Please refer to Fig. 25.3 to see the differences between DIN EN 1992-
1-1 and EN 1992-1-1 (for t=500000 h).
0.45 <1%
0.50 -
0.55 - 1.0 1.2
0.60 - 1.2 2.5 2.8
0.65 - 1.3 2.0 4.5 5.0
0.70 - 1.0 2.0 3.0 6.5 7.0
0.75 - 1.2 2.5 3.0 4.5 9.0 10.0
25 %
CLASS R1 new (old R1)
Relaxation losses Δσpr / σp for t=500000 h
20 %
15 %
13.0 % abZ
CLASS R2 new (old R3)
10 %
9.0 %
old R2
6.5 %
5% 4.5 %
2.5 %
0%
0.50 0.55 0.60 0.65 0.70 0.75 0.80
μ = σp / ƒpk
Figure 25.3: The relaxation loss differences according to EC2 and abZ
(DiBt)
t = 1000 h
AQB reduces the initial stress according to DIN 1045-1: DIN 1045-1 8.7.3 (6)
Ep
εcs · Ep + Δσpr + · ϕ(t, t0 ) · σc,QP
Ecm
ΔPc+s+r = Ap DIN 1045-1 8.7.3 (6): ΔPc+s+r Time de-
Ep Ap Ac
pendent losses of prestress
1+ · 1+ 2
· zcp [1 + 0.8 · ϕ(t, t0 )]
Ecm Ac c
Creep and shrinkage is not taken into account therefore we have:
εcs · Ep = 0
Ep
· ϕ(t, t0 ) · σc,QP = 0
Ecm
1 + 0.8 · ϕ(t, t0 ) = 1
31.82
Δσp,c+s+r =
195000 28.5 · 10−4 0.9926
1+ · · 1+ · 0.39012
34077 0.9926 0.08214
31.82
Δσp,c+s+r = = 30.402 MP
1.046644
Now we have:
DIN 1045-1 8.7.3 (6) σp0 = 1281.755 − 0.3 · 30.4020 = 1272.631 MP
Now with the reduced stress the μ value is calculated again (iteration):
σp0 1272.631
μ= = = 0.719 ≥ μmn = 0.55
ƒpk 1770
εk0 , εky and εkz - tendon strain Δσpr,c = (εk0 + εky · zcp − εkz · ycp ) · Ep
25.5 Conclusion
This example shows the calculation of the time dependent losses due to relaxation. It has been shown
that the reference solution and the AQB solution are in very good agreement.
25.6 Literature
[1] DIN EN 1992-1-1/NA: Eurocode 2: Design of concrete structures, Part 1-1/NA: General rules and
rules for buildings - German version EN 1992-1-1:2005 (D), Nationaler Anhang Deutschland -
Stand Februar 2010. CEN. 2010.
[2] F. Fingerloos, J. Hegger, and K. Zilch. DIN EN 1992-1-1 Bemessung und Konstruktion von
Stahlbeton- und Spannbetontragwerken - Teil 1-1: Allgemeine Bemessungsregeln und Regeln
für den Hochbau. BVPI, DBV, ISB, VBI. Ernst & Sohn, Beuth, 2012.
[6] EN 1992-1-1: Eurocode 2: Design of concrete structures, Part 1-1: General rules and rules for
buildings. CEN. 2004.
[17] DiBt. Allgemeine bauaufsichtliche Zulassung. Deutsches Institut für Bautechnik. 2015.
Overview
NEd
My,Ed
qz
My,Ed
NEd
Mz,Ed
tƒ = 16.0 mm
t = 10.2 mm
A = 115.5 cm2
y = 48197 cm4
z = 2142 cm4
T = 88.57 cm4
λy 0.195 0.195
χy 1.0 1.0
λz 0.927 0.927
ϕz 1.054 1.054
χz 0.644 0.644
μy 1.000 1.000
μz 0.937 0.937
y 1.138 1.138
z 1.500 1.500
λ0 0.759 0.757
C1 1.194 1.2
MzE d = 25 kNm
NEd = 25 kNm
Tab. 5.5: Classification of cross-section The cross-section is classified as Class 1, as demonstrated in [18].
π 2 Ey
6.3.1.2 (1): Ncr is the elastic critical Ncr,y = = 71035 kN
force for the relevant buckling mode L2
v
uAƒ
y
λy = = 0.195
t
6.3.1.2 (1): λ non dimensional slender-
ness for class 1 cross-sections Ncr,y
2
h i
z = 0.5 1 + αz λz − 0.2 + λz
for rolled I-sections with h / b > 1.2 and buckling about z-z axis →
buckling curve b
6.3.1.2 (2): Table 6.1: Imperfection fac- for buckling curve b → αz = 0.34
tors for buckling curves
z = 1.054
1
6.3.1.2 (1): Eq. 6.49: χz reduction fac- χz = r = 0.644 ≤ 1.0
2
tor for buckling
z + 2z − λz
Auxiliary terms:
NEd
1−
Ncr,y
Annex A: Tab. A.1: Interaction factors μy = = 1.0
kj (6.3.3(4)), Auxiliary terms NEd
1 − χy
Ncr,y
1 The sections mentioned in the margins refer to EN 1993-1-1:2005 [7] unless other-
wise specified.
NEd
1−
Ncr,z
μz = = 0.937
NEd
1 − χz
Ncr,z
Wp,y
y = = 1.138 ≤ 1.5
We,y
Wp,z
z = = 1.568 > 1.5 → z = 1.5
We,z
The general formula for combined end moments and transverse loads
is used here.
π 2 Ey |δz |
NEd
Cmy,0 = 1 + −1 Annex A: Tab. A.2: Equivalent uniform
2
L My,Ed,rght Ncr,y moment factors Cm,0
δz = 3.33 mm
Cmy,0 = 1.001
NEd
Cmz,0 = 0.79 + 0.21ψz + 0.36(ψz − 0.33) = 0.771
Ncr,z
My,Ed,rght A
Annex A: Tab. A.1: εy for class 1 cross- εy = = 2.383
section NEd We,y
p
αLT εy
Cmy = Cmy,0 + 1 − Cmy,0 = 1.001
1 + αLT εy
p
αLT
CmLT = C2my v = 1.139 ≥ 1.0
NEd NEd
u
1− 1−
t
Ncr,z Ncr,T
6.3.2.2: Lateral torsional buckling The general case method is chosen here.
curves - General case
Mcr = 1068 kNm, determined by eigenvalue analysis
v v
uW
t p,y ƒy t 2194 · 10−6 · 235 · 106
u
6.3.2.2 (1): λLT non dimensional slen- λLT = = = 0.695
derness for lateral torsional buckling Mcr 1079 · 103
2
h i
LT = 0.5 1 + αLT λLT − 0.2 + λLT
6.3.2.2 (1): Table 6.4: Recommenda- for rolled I-sections and h / b > 2 → buckling curve b
tion for the selection of ltb curve for
cross-sections using Eq. 6.56 for buckling curve b → αLT = 0.34
6.3.2.2 (2): Table 6.4: Recommenda-
tion values for imperfection factors αLT
LT = 0.825
for ltb curves
1
6.3.2.2 (1): Eq. 6.56: χLT reduction fac- χLT = r
2
tor for ltb
LT + 2LT − λLT
1.6 1.6
2
Cyy = 1 + y − 1 · C2my λm · C2my λm
2− −
y y
NEd
· − bLT
Nc,Rk
γM1
We,y
Cyy = 0.981 ≥ = 0.879
Wp,y
2
λ0 My,Ed
cLT = 10 · αLT · 4
= 0.471
Mp,y,Rk
5 + λz C
my · χLT,mod
γM1
2
Cmy λm 2
NEd
Cyz = 1 + (z − 1)
2 − 14
− c LT
z5 Nc,Rk
γM1
v
u z We,z
Cyz = 0.862 ≥ 0.6 = 0.439
t
y Wp,z
λ0 My,Ed Mz,Ed
dLT = 2 · αLT · 4 Mp,y,Rk Mp,z,Rk
0.1 + λz C
my · χLT,mod Cmz
γM1 γM1
= 0.348
2
C2my λm NEd
Czy = 1 + y − 1
2−
14 − dLT
y5 Nc,Rk
γM1
t y We,y
v
Czy = 0.842 ≥ 0.6 = 0.459
z Wp,y
λ0 My,Ed
eLT = 1.7 · αLT · 4
= 0.721
Mp,y,Rk
0.1 + λz C
my · χLT,mod
γM1
1.6 1.6
2
Czz = 1 + (z − 1) 2− · C2mz λm − · C2mz λm − eLT
z z
NEd
·
Nc,Rk
γM1
We,z
Czz = 1.013 ≥ = 0.667
Wp,z
NEd C
mLT Cmy · My,Ed
+ μy
Nc,Rk NEd Mp,y,Rk
χLT,mod
χy 1− Cyy ·
γM1 Ncr,y γM1
v
u Cmz Mz,Ed
t z
+0.6
NEd Mp,z,Rk
y
1− Cyz
Ncr,z γM1
6.3.3 (4): Eq. 6.61: Members which are
subjected to combined bending and ax-
ial compression should satisfy:
→ Satisfactory
v
NEd u
t y CmLT Cmy · My,Ed
+ μz 0.6
Nc,Rk NEd Mp,y,Rk
z χLT,mod
χz 1− Czy
γM1 Ncr,y γM1
Cmz Mz,Ed
+
NEd Mp,z,Rk
1− Czz
Ncr,z γM1
6.3.3 (4): Eq. 6.62: Members which are
subjected to combined bending and ax-
ial compression should satisfy:
v
500 t 1.138 1.139
u
+ 0.937 +0.6
2715 1.500 0.821
0.644
1.0
= 0.868 ≤ 1
→Satisfactory
26.5 Conclusion
This example shows the ckeck for lateral torsional buckling of steel members. The small deviations
that occur in some results come from the fact that there are some small differences in the sectional
values between SOFiSTiK and the reference solution. Therefore, these deviations are of no interest for
the specific verification process. In conclusion, it has been shown that the results are reproduced with
excellent accuracy.
26.6 Literature
[7] EN 1993-1-1: Eurocode 3: Design of concrete structures, Part 1-1: General rules and rules for
buildings. CEN. 2005.
[18] ECCS Technical Committee 8 Stability. Rules for Member Stability in EN 1993-1-1, Background
documentation and design guidelines. Tech. rep. No 119. European Convention for Constructional
Steelwork (ECCS), 2006.
Overview
hh
h
d d1
bh
• Design stress-strain curves for concrete and reinforcement (Section 3.1.7, 3.2.7)
The design stress-strain diagram for reinforcing steel considered in this example, consists of an inclined
top branch, as presented in Fig. 27.2 and as defined in DIN EN 1992-1-1:2004 [1] (Section 3.2.7).
A
σ
ƒtk,c = 525N/ mm2
ƒyk
ƒyd = ƒyk / γs
B A Idealised
B Design
εd = 25 0 /00
Figure 27.2: Idealised and Design Stress-Strain Diagram for Reinforcing Steel
25 110 25
1 1
20
75
35
1 1
20
160
d1 = 5.0 cm
DCE-EN25: Shear between web and flanges of Hollow CS acc. DIN EN 1992-2
SOFiSTiK AG * Bruckmannring 38 * 85764 Oberschleißheim Page 1
SOFiSTiK 2016-3 WINGRAF - GRAPHICS FOR FINITE ELEMENTS (V 18.03) 2016-01-15
Single span 155.0
beam 155.0 155.0 155.0 155.0 155.0 155.0 155.0 155.0 155.0
The system with its loading as well as the moment and shear force are shown in Fig. 27.4-27.7. The
0.00
reference calculation steps are presented in the next section and the results are given in Table 27.2.
155.0 155.0 155.0 155.0 155.0 155.0 155.0 155.0 155.0 155.0
SOFiSTiK AG - www.sofistik.de SOFiSTiK AG - www.sofistik.de
0.00
All loads (in components), Loadcase 1 sum_PZ=1550.0 kN . , (1 cm 3D = unit) Beam line load (force) in global Z (Unit=160.2 M 1 : 41
Y X
kN/m ) (Max=155.0)
Z
Figure 27.4: Loaded Structure (Pg )
0.00
1001 1002Loadcase 1 sum_PZ=1550.0
All loads (in components), 1003 kN 1004
. , (1 cm 3D = unit)1005 1006
Beam line load (force) 1007
in global Z (Unit=160.2 1008 1009 1010 M 1 : 41
SOFiSTiK YAGX * Bruckmannring
kN/m 38 * 85764 Oberschleißheim
) (Max=155.0) Page 1
SOFiSTiK Z2016-3 WINGRAF - GRAPHICS FOR FINITE ELEMENTS (V 18.03) 2016-02-17
697.5
697.5
1240
1240
SINGLE SPAN BEAM
1628
1628
1860
1860
1938
0.00
1001 1002 1003 1004 1005 1006 1007 1008 1009 1010
697.5
697.5
Figure 27.5: Resulting Bending Moment My
1240
1240
2.00
-2.00
1628
1628
0.00 2.00 4.00 6.00 8.00 10.00 m
1860
1860
1938
-500.0
-400.0
-300.0
-200.0
-100.0
0.00 2.00 4.00 6.00 8.00 10.00 m
SOFiSTiK AG - www.sofistik.de
0.00
1001 1002 1003 1004 1005 1006 1007 1008 1009 1010
100.0
-2.00
2.00
775.0
620.0
465.0
310.0
Beam Elements , Torsional moment Mt, Loadcase 1 sum_PZ=1550.0 kN . , 1 cm 3D = 500.0 kNm (Min=-500.0) (Max=500.0) M 1 : 41
SOFiSTiK AG - www.sofistik.de
Y X
Number of element, Beam Elements(Max=1010), Nodal sequences
Z
0.00
1001 1002 1003 1004 1005 1006 1007 1008 1009 1010
-155.0
-310.0
-465.0
-620.0
-775.0
Y X
At beam 1001
Beam Elements , Shear force Vz, Loadcase 1 sum_PZ=1550.0 kN . , 1 cm 3D = 500.0 kN (Min=-775.0) (Max=775.0) SOF. Ref. M 1 : 41
Number of element, Beam Elements(Max=1010), Nodal sequences
Z
Material:
Design loads
ω ≈ 0.053973, ζ ≈ 0.9658, ξ = 0.07833 (interpolated) Tab. 9.2 [3]: ω−Table for up to C50/ 60
without compression reinforcement
1
As1 = · (ω · b · d · ƒcd + NEd ) NEd = 0
σsd
1
As1 = · (0.053973 · 1.6 · 0.7 · 17.0) · 1002 = 22.52 cm2
456.52
z = ζ · d = 0.9658 · 0.70 m ≈ 67.00 cm
1 The tools used in the design process are based on steel stress-strain diagrams, as
defined in [1] 3.2.7:(2), Fig. 3.8, which can be seen in Fig. 27.2.
2 The sections mentioned in the margins refer to DIN EN 1992-1-1:2004 (German Na-
MEd,=1.0 MEd,=0.0 hƒ · b
ΔFd = − ·
z z hƒ · bh
ΔFd
352.35
In AQB output Ed,V = τ − V Ed,V = = 0.176 kN/ m2 = 1.76 MP
20 · 100
Checking the maximum Rd,m value to prevent crushing of the
struts in the flange
DIN EN 1992, 6.2.4 (4),Eq. 6.22 Ed,V ≤ Rd,m = ν · ƒcd · snθƒ · cosθƒ
ν = ν1
b → hƒ , z → Δ, c = 0.5
σcd
1/ 3
VRd,cc = c · 0.48 · ƒck · 1 − 1.2 · · hƒ · Δ
ƒcd
0
VRd,cc = 0.5 · 0.48 · 301/ 3 · 1 − 1.2 · · 0.20 · 1.0
17.00
The internal forces must be in equilibrium, see Fig. 27.8!!! If the internal
forces are not in equilibrium state, then AQB will print erroneous results.
Therefore the internal forces must be taken from the system (calculated
by using e.g. ASE or STAR2).
z
σz τy = τy
x τz = τz
y τz τyz = τzy
τzy τz
τyz
σ σ
τy
σy τy
σz
1.2
cotθ = = 2.0807 ≥ 1.75 → 1.75
1 − 0.745/ 1.76
AQB is checking and iterating if the cotθ is less than (<) the maximum
cotθ defined in the norm. If yes, then this value will be taken by AQB.
1 1
tn θ = = = 0.5714 → θ = 29.74488◦
cotθ 1.75
Rd,m = 0.75 · 17.00 · sn29.7448 · cos29.7448 = 5.492 MP
VRd,m = 1098.461 kN
Rd,c = k · ƒctd
Ed · hƒ
DIN EN 1992, 6.2.4 (4),Eq. 6.21 sƒ ,V =
cotθƒ · ƒyd
1.76 · 0.20
sƒ ,V = · 1002 = 4.63 cm2
1.75 · 434.78
• Calculating the necessary torsional reinforcement:
HINT: Please note that in AQB the torsional stress IS NOT calculated
by using the formula:
Mt
τt =
2 · Ak · teƒ ƒ
τt = Mt · 1/ Wt
• Total reinforcement:
27.5 Conclusion
This example is concerned with the calculation of the shear between web and flanges of a Hollow CS. It
shows partially the work-flow how AQB calculates the shear between web and flanges.
Please note that it is very difficult to show all steps how AQB calculates internally the τ and σ stresses,
therefore some steps are skipped. The reference example is just an approximation that shows the results
by using hand-calculation. It has been shown that the results calculated by using hand-calculation and
the AQB module are reproduced with very good accuracy.
The τt and τ values between cross-section points are interpolated and calculated by using Finite
Element Method.
Figure 27.9: The τt stresses (vector) calculated by using Finite Element Method
Figure 27.10: The τt stresses (fill) calculated by using Finite Element Method
27.6 Literature
[1] DIN EN 1992-1-1/NA: Eurocode 2: Design of concrete structures, Part 1-1/NA: General rules and
rules for buildings - German version EN 1992-1-1:2005 (D), Nationaler Anhang Deutschland -
Stand Februar 2010. CEN. 2010.
[3] K. Holschemacher, T. Müller, and F. Lobisch. Bemessungshilfsmittel für Betonbauteile nach Eu-
rocode 2 Bauingenieure. 3rd. Ernst & Sohn, 2012.
[19] DIN EN 1992-2: Eurocode 2: Design of concrete structures - Part 2, Part 2: Concrete Bridges -
Design of Detailing Rules - German version EN 1992-2:2005 (D), Nationaler Anhang Deutschland
- 2012. CEN. 2012.
Overview
Slab height=19 cm
12 482 24 382 12
Field 1 Field 2
A B C
• Design stress-strain curves for concrete and reinforcement (Section 3.1.7, 3.2.7)
The verification of the BEMESS results will be examined. A complete and detailed hand-calculation of
the results is not possible because of described BEMESS-strategy, which should be here to exhaustive.
For this reason, some results (e.g. internal forces) will be taken as outputted and further used in the
hand-calculation.
The design stress-strain diagram for reinforcing steel considered in this example, consists of an inclined
top branch, as presented in Fig. 28.2 and as defined in EN 1992-1-1:2004 [6] (Section 3.2.7).
A
σ
ƒtk,c = 525N/ mm2
ƒyk
ƒyd = ƒyk / γs
B A Idealised
B Design
εd = 25 0 /00
Figure 28.2: Idealised and Design Stress-Strain Diagram for Reinforcing Steel
-37.2
A
B
33.6
C
19.0
The system with its loading are shown in Fig. 28.4-28.10. The reference calculation steps are presented
in the next section and the results are given in Table 28.4.
16
.1
16.1
Field 2
Field 1
16.1
16
.1
8.57
16
16.1
.1
Field 2
Field 1
16.1
16
8.
8.57
.1
57
Figure 28.6: Layer Design - Quad 10026, Loadcase 1001, st = 5.59 cm2 / m
Figure 28.7: Layer Design - Quad 10726, Loadcase 1002, st = 5.02 cm2 / m
Figure 28.8: Layer Design - Quad 80376, Loadcase 1003, st = 2.73 cm2 / m
Summary of results:
16.1
16
.1
8.
8.57
57
Field 2
Field 1
16
16.1
.1
8. 8.57
57
Material:
To make the example more simple, the slab will be designed only for
the maximum and minimum moment mEd (Field 1, Field 2 and over the
middle support). The reduction of the moment over the middle support
will be neglected in this example.
|mEd,B,Red |
μEds =
b · d2 · ƒcd
37.16 · 10−3
μEds =
1.0 · 0.162 · 13.33
μEds = 0.1089
From tables:
μ ω ξ ζ σsd [MPa]
0.108867 0.11575 0.14298 0.940522 452.69
1
s,req = ω · b · d · ƒcd ·
σsd
1
s,req = 0.11576 · 1.00 · 0.16 · 13.33 ·
452.70
s,req = 5.45 cm2 / m
defined in [1] 3.2.7:(2), Fig. 3.8, which can be seen in Fig. 28.2.
2 The sections mentioned in the margins refer to EN 1992-1-1:2004 [6], unless other-
wise specified.
|mEd,B,Red |
μEds =
b · d2 · ƒcd
33.69 · 10−3
μEds =
1.0 · 0.162 · 13.33
μEds = 0.0987
From tables:
μ ω ξ ζ σsd [MPa]
0.0987 0.10429 0.12882 0.9464 455.25
1
s,req = ω · b · d · ƒcd ·
σsd
1
s,req = 0.10429 · 1.00 · 0.16 · 13.33 ·
455.25
s,req = 4.88 cm2 / m
|mEd,B,Red |
μEds =
b · d2 · ƒcd
19.09 · 10−3
μEds =
1.0 · 0.162 · 13.33
μEds = 0.05593
From tables:
μ ω ξ ζ σsd [MPa]
0.05593 0.05774 0.08222 0.968533 456.52
1
s,req = ω · b · d · ƒcd ·
σsd
1
s,req = 0.05774 · 1.00 · 0.16 · 13.33 ·
456.52
s,req = 2.698 cm2 / m
For each reinforcement layer: With the use of internal forces in local
element direction, the internal forces in main direction and the accom-
payning angle is calculted with the following equation
m + myy r
m/ = ± 0.5 · (myy − m )2 + 4 · m2y
2
my
tn 2φ0 = 2 ·
m − myy
The internal lever arm is calculated separately for the internal forces in
main moment direction and for the internal forces in main normal forces
direction. This is done with the theory explained in the paper from Prof.
Dr. Ing. Ulrich P. Schmitz [20]. The program choose the unfavorable
lever arm from both results for the next analysis step lever arm z (This
is why for each layer of reinforcement two lever arms are calculated
within the program).
Calculating N :
n m
N = +
2 z
0.00 37.16
N = +
2 0.1514
N = 245.442 kN/ m
Calculating Ny :
nyy myy
Ny = +
2 z
0.00 7.43
Ny = +
2 0.1514
Ny = 49.075 kN/ m
Calculating Ny :
ny ny
Ny = +
2 z
0.00 0.00
Ny = +
2 0.1514
Ny = 0.00 kN/ m
The next step is the transformation of the panel forces N , Ny , Ny into
the main forces N and N :
N + Ny r
N = + 0.5 · (Ny − N )2 + 4 · N2y
2
245.44 + 49.07 Æ
N = + 0.5 · (49.07 − 245.4)2 + 4 · 0.002
2
N = 147.256 + 98.165 = 245.421
N + Ny r
N = − 0.5 · (Ny − N )2 + 4 · N2y
2
245.44 + 49.07 Æ
N = − 0.5 · (49.07 − 245.4)2 + 4 · 0.002
2
N = 147.256 − 98.165 = 49.091
Ny
tn 2φ0 = 2 ·
N − Ny
0.00
tn 2φ0 = 2 ·
245.442 − 49.07
tn 2φ0 = 0.00 → φ = 0
Required reinforcement:
N2
k=
N1
49.091
k= = 0.20
245.41
k ≥ tn(α + π/ 4) · tn α = 0
N − N
Z = N + · sin 2α · (1 − tn α)
2
245.421 − 49.091
Z = 245.421 + · sin 0 · (1 − tn 0)
2
Z = 245.421
Z
s =
σsd
245.421
s = = 5.375 cm2 / m
456.52
Calculating N :
n m
N = +
z z
0.00 33.70
N = +
0.1514 0.1514
N = 222.589 kN/ m
Calculating Ny :
nyy myy
Ny = +
z z
0.00 6.74
Ny = +
0.1514 0.1514
Ny = 44.51 kN/ m
Calculating Ny :
ny ny
Ny = +
z z
0.00 0.00
Ny = +
0.1514 0.1514
Ny = 0.00 kN/ m
The next step is the transformation of the panel forces N , Ny , Ny into
the main forces N and N :
N + Ny r
N = + 0.5 · (Ny − N )2 + 4 · N2y
2
222.589 + 44.51 Æ
N = + 0.5 · (44.51 − 222.589)2 + 4 · 0.002
2
N = 133.5495 + 89.0395 = 222.589
N + Ny r
N = − 0.5 · (Ny − N )2 + 4 · N2y
2
222.589 + 44.51 Æ
N = − 0.5 · (44.51 − 222.589)2 + 4 · 0.002
2
N = 133.5495 − 89.0395 = 44.51
Ny
tn 2φ0 = 2 ·
N − Ny
0.00
tn 2φ0 = 2 ·
222.589 − 44.51
tn 2φ0 = 0.00 → φ = 0
Required reinforcement:
N2
k=
N1
44.51
k= = 0.199
222.589
k ≥ tn(α + π/ 4) · tn α = 0
N − N
Z = N + · sin 2α · (1 − tn α)
2
222.589 − 44.51
Z = 222.589 + · sin 0 · (1 − tn 0)
2
Z = 222.589
Z
s =
σsd
222.589
s = = 4.875 cm2 / m
456.52
Field 2 (QUAD 80376):
Calculating N :
n m
N = +
z z
0.00 19.09
N = +
0.1514 0.1549
N = 123.24 kN/ m
Calculating Ny :
nyy myy
Ny = +
z z
0.00 3.86
Ny = +
0.1549 0.1549
Ny = 24.91 kN/ m
Calculating Ny :
ny ny
Ny = +
z z
0.00 0.00
Ny = +
0.1549 0.1549
Ny = 0.00 kN/ m
The next step is the transformation of the panel forces N , Ny , Ny into
the main forces N and N :
N + Ny r
N = + 0.5 · (Ny − N )2 + 4 · N2y
2
123.24 + 24.91 Æ
N = + 0.5 · (24.91 − 123.24)2 + 4 · 0.002
2
0.00
tn 2φ0 = 2 ·
123.24 − 24.91
tn 2φ0 = 0.00 → φ = 0
Required reinforcement:
N2
k=
N1
24.91
k= = 0.202
123.24
k ≥ tn(α + π/ 4) · tn α = 0
N − N
Z = N + · sin 2α · (1 − tn α)
2
123.24 − 24.91
Z = 123.24 + · sin 0 · (1 − tn 0)
2
Z = 123.24
Z
s =
σsd
123.24
s = = 2.70 cm2 / m
456.52
28.5 Conclusion
This example shows the calculation of the required reinforcement for a one-way slab under bending. It
has been shown that the results are reproduced with very good accuracy.
28.6 Literature
[1] DIN EN 1992-1-1/NA: Eurocode 2: Design of concrete structures, Part 1-1/NA: General rules and
rules for buildings - German version EN 1992-1-1:2005 (D), Nationaler Anhang Deutschland -
Stand Februar 2010. CEN. 2010.
[4] Beispiele zur Bemessung nach Eurocode 2 - Band 1: Hochbau. Ernst & Sohn. Deutschen Beton-
und Bautechnik-Verein E.V. 2011.
[6] EN 1992-1-1: Eurocode 2: Design of concrete structures, Part 1-1: General rules and rules for
buildings. CEN. 2004.
[20] Ulrich P. Schmitz. “Biegebemessung nach der neuen E-DIN 1045-1”. In: Bauinformatik JOURNAL
(1999), pp. 49–51.
Overview
beƒ ƒ
2.4 m
qd + gd
0 = 14.04 m
Figure 29.1: Sketch of the slab and longitudinal system of single t-beam member
sections, and superimposes these with the corresponding forces of the eccentric beam to result in the
design values of the bending moment MEd and shear force VEd .
beƒ ƒ
hƒ
Quad element
h Beam element
b0 b00
b0
Figure 29.2: Cross sectional view of the model out of quad and beam elements
The sectional properties and loadings of the member are given in table 29.1. The analysis gives results
of the individual quad and beam elements. DECREATOR is then used to integrate all forces at the
explicit design sections for use in design. Since, in this model, the centroid of the chosen cross section
lies on the design element, no transformation of forces from design element axis to the centroid takes
place during computation. However, this full cross section will be used for design. The final analysis and
design results are given in Table 29.2, the reference calculation procedures of which are presented in
the following section.
beƒ ƒ = 240 cm
b00
= 20.0 cm
b0 = 6.0 cm
b = 26.0 cm
Summarized in Table 2 are the results of the reference example in [4] (Ref.-2), those obtained based on
design tables (Ref.-1) as presented in the next section, and those from SOFiSTiK computations. For sim-
plification purposes the reference example uses rough approximate procedures in obtaining reinforce-
ment areas, while the hand calculation and the SOFiSTiK results come from iterative approximations
and are as close to the exact requirements as possible.
In this regard, in calculating the reinforcement area in Ref.-2, the full thickness of the flange was assumed
to roughly approximate the compressive zone. The moment arm is thereby reduced and hence it resulted
in a larger reinforcement area. In a similar manner, the shear reinforcements in Ref.-2 are larger since
the example selects a suitable reinforcement diameter and spacing in advance to compute the area,
Design of the T-beam for bending and shear are carried out using the T-
cross section, as shown in figure 29.3. References on the side margins
can be further reviewed for notations and explanations.
beff,1
b1 d
beƒ ƒ = beƒ ƒ , + b ≤ b
P
5.3.2.1: (3) Eq. (5.7) and Fig. (5.3) for
notations
5.3.2.1: (3) Eq. (5.7a) beƒ ƒ , = 0.2 · b + 0.1 · 0 ≤ 0.2 · 0
b = 0.5 · (beƒ ƒ − b )
Using design table for T-beams: μs - Table for T-beams (Tab. 9.4) [3]
MEds 1086.0
μEds = = = 0.0594
beƒ ƒ · d2 · ƒcd 2.4 · 0.622 · 19.83
Design section at: = 0.5 · 0.2 + 0.62 = 0.72 m Distance from center of elastomer sup-
port with thickness b=200 mm [4]
Design shear force: VRd,s = 272.0 kN
VRd,s
(ASW / S) =
z · ƒyd · cotθ
272.0
(ASW / S) = = 9.33 cm2 / m
0.558 · 434.78 · 1.2
29.5 Conclusion
It has been demonstrated in this example that, using Design Elements in SOFiSTiK, a structural member
modeled out of different finite element types can be designed with a unifying full cross section. The
critical internal forces were well approximated and reasonable design results were obtained.
29.6 Literature
[3] K. Holschemacher, T. Müller, and F. Lobisch. Bemessungshilfsmittel für Betonbauteile nach Eu-
rocode 2 Bauingenieure. 3rd. Ernst & Sohn, 2012.
[4] Beispiele zur Bemessung nach Eurocode 2 - Band 1: Hochbau. Ernst & Sohn. Deutschen Beton-
und Bautechnik-Verein E.V. 2011.
[6] EN 1992-1-1: Eurocode 2: Design of concrete structures, Part 1-1: General rules and rules for
buildings. CEN. 2004.
Overview
Z
XY N = 300 kN
q = 4 kN/ m
ƒy = 235 N/ mm2 c = 32 mm
ν = 0.3 r = 18 mm
G = 81000 N/ mm2 tƒ = 10 mm
y = 82.8 mm
z = 49.8 mm
y = 3690 cm4
z = 1340 cm4
t = 21 cm4
λy 1.029 1.029
λz 0.855 0.855
Nd = 300 kN
q = 4.0 kN/ m
Mz,d = 0
Buckling lengths:
Lcr,y = 8.00 m
Lcr,z = 4.00 m
ξ = 1.35
Characteristic values:
y = 8.28 cm
z = 4.98 cm
zp = −9.5 cm
Critical loading:
Ncr - elastic critical force for the relevant Ncr,z = π 2 · 21000 · 1340/ 4002 = 1736 kN
buckling mode based on the gross cross
sectional properties c2 = (108000 + 0.039 · 4002 · 21.0)/ 1340 = 178.4 cm2
p
Mcr - elastic critical moment for lateral- Mcr = 1.35 · 1736 · 178.4 + 0.25 · 9.52 + 0.5 · (−9.5) · 102
torsional buckling
Mcr = 220.9 kNm
1 The sections mentioned in the margins refer to DIN EN 1993-1-1:2005 [21] unless
otherwise specified.
-300.0
7.00
7.00
6.00
6.00
5.00
5.00
4.00
4.00
32
.0
3.00
3.00
2.00
2.00
1.00
1.00
-300.0
0.00
0.00
p
λLT = 100.9/ 220.9 = 0.676 λLT - non dimensional slenderness for
lateral torsional buckling
h/ b = 190/ 200 = 0.95 < 2.0 → Class b
1
χLT = r
2
LT + 2LT − β · λLT
2
h i
LT = 0.5 · 1 + αLT · (λLT − λLT,0 ) + β · λLT
χLT 1
Mb,Rd = Mp,y,Rd · ·
ƒ γM1
ƒ = 1 − 0.5 · (1 − kc ) · 1 − 2.0 · (λLT − 0.82 ) ƒ - the value f may be defined in National
Annex, see EC 3, §6.3.2.3(2)
ƒ = 0.972
0.88 1
Mb,Rd = 100.9 · ·
0.972 1.1
Mb,Rd = 83.045 kNm mb,Rd - design buckling resistance mo-
ment - EC 3, §6.3.2.1, Eq. 6.55
Equivalent uniform moment factors:
αh = Mh / Ms = 0.00
ψ = 1.00
αh = Mh / Ms = 0.75
ψ = 0.00
z y
y z
ψ · Mh ψ · Mh
ψ · Ms = 0.75 · Mh
Ms Mh
Mh
1.1
kyy = 0.95 · 1 + 1.029 − 0.2) · 300 ·
0.580 · 1264.3
= 1.304
1.1
≤ 0.95 · 1 + 0.8 · 300 · = 1.292
0.580 · 1264.3
≤ 1.80
30.5 Conclusion
This example shows the check for lateral torsional buckling of steel members. The small deviations that
occur in some results come from the fact that there are some small differences in the sectional values
and elastic critical loadings (Mcr , Ncr ). Therefore, these deviations are of no interest for the specific
verification process. In conclusion, it has been shown that the results are reproduced with excellent
accuracy.
30.6 Literature
[8] Schneider. Bautabellen für Ingenieure. 21th. Bundesanzeiger Verlag, 2014.
[21] DIN EN 1993-1-1:2005 Eurocode 3: Design of steel structures, Part 1-1: General rules and rules
for buildings - Deutsche Fassung EN 1993-1-1:2005 + AC:2009. CEN. 2010.
Overview
bƒ 1
tƒ 1
t
h
y T
tƒ 2
bƒ 2
This benchmark example is used to verify and compare the SOFiSTiK results with the ECCS reference
example [22].
5500
17 -5500
5.5
17
5.5
17
5.5
17
5.5
17
5.5
17
5.5
17
5.5
17
5.5
17
5.5
17
5.5
17 -5500
5.5
00 Loading
0.00 My =
-4.002.00N · e
N,y -2.00 4.00 0.00 6.00
N
2 00 0 00 2 00 4 00
ν = 0.3 tƒ2 = 20 mm
G = 81000 N/ mm2 t = 10 mm
z = 153200.0 cm4
NEd = 5500 kN
1. CROSS-SECTION RESISTANCE
bƒ 1 = bƒ 2 = b − 2 · t = 600 − 2 · 10 = 580 mm
h = h − tƒ 1 − tƒ 2 = 600 − 10 − 20 = 570 mm
bƒ 1 580
= = 58.0 > 42 · ϵ = 42 · 0.924 = 38.8 (Class 4)
tƒ 1 10
bƒ 2 580
= = 29.0 < 42 · ϵ = 38.8 (Class 3)
tƒ 2 20
Web (compression):
h 570
= = 57.0 > 42 · ϵ = 38.8 (Class 4) Cross-section classification, EN 1993-
t 10 1-1, Table 5.2
otherwise specified.
1 1 1
2
4
2
4
2
4
3
Figure 31.3: Effective area for uniform compression
tƒ 1 + tƒ 1 h + tƒ 2
Sy = b · tƒ 1 · h − + 2 · h · t ·
2 2
1 10 + 20
Sy = 600 · 10 · 600 − +2
29400 2
570 + 20
· 570 · 10 · 600 −
2
Sy 6.873 · 106
rt = = = 233.8 mm
A 29400.0
where:
Sy is the first moment of area of the gross cross section with respect to
the centroid of the lower flange (y-y axis),
rt is the distance from the centroid of the lower flange to the centroid
of the gross cross-section
ψ = 1.0 → kσ = 4.0
bƒ 1 580
Plate elements without longitudinal stiff- λp = p = p = 1.105
eners, EN 1993-1-5, 4.4 tƒ 1 · 28.4 · ϵ · kσ 10 · 28.4 · 0.924 · 4.0
p
λp = 1.105 > 0.5 + 0.085 − 0.055 · Ψ
p
= 0.5 + 0.085 − 0.055 · 1 = 0.673
ψ = 1.0 → kσ = 4.0
h 570
λp = p = p = 1.086 Plate elements without longitudinal stiff-
t · 28.4 · ϵ ·
kσ 10 · 28.4 · 0.924 · 4.0 eners, EN 1993-1-5, 4.4
p
λp = 1.086 > 0.5 + 0.085 − 0.055 · 1 = 0.673
Aeƒ ƒ = A − (ƒ · tƒ 1 + 2 · · t )
tƒ 1 + tƒ 2 10 + 20
rƒ = h − − rt = 600 − = 351.2 mm
2 2
tƒ 2
r = h + − rT − be1, −
2 2
10 151.3
r = 570 + − 233.8 − 209.3 − = 61.2 mm
2 2
2 · r · · t + rƒ · ƒ · tƒ 1
eN,y =
Aeƒ ƒ
where:
eN,y is the shift of centroid of the effective area relative to the centre
of gravity of the gross cross section determined assuming uniform
axial compression.
rTeƒ ƒ ,N is the distance from the centroid of the bottom flange to the
centroid of the effective cross-section under uniform compression.
1 1 1
3
Figure 31.4: Effective area for bending
Aeƒ ƒ = A − ƒ · tƒ 1
rƒ · ƒ · tƒ 1 351.2 · 159.5 · 10
ΔrT,m = = = 20.1 mm
Aeƒ ƒ 24778.1
ƒ · tƒ3
!
eƒ ƒ ,y
= y + Aeƒ ƒ · ΔrT,M
2 − + ƒ · tƒ · (rƒ + ΔrT,M )2
12
159.5 · 103
eƒ ƒ ,y
= 1.748 · 109 + 27804.9 · 20.12 − +
12
where:
eƒ is the effective second moment of area (cross section under pure
ƒ ,y
bending) with respect to y-y considering the effective width of the
upper flange.
The effective section moduli at the upper and lower edge of the girder’s
web, Weƒ
and Weƒ are, respectively:
ƒ ,y,1 ƒ ,y,2
eƒ ƒ ,y
Weƒ =
ƒ ,y,1 tƒ 2
h + − rTeƒ ƒ ,M
2
1.540 · 109
Weƒ ƒ ,y,1
= = 4.20 · 106 mm3
20
570 − − 213.6
3
eƒ ƒ ,y
Weƒ =
ƒ ,y,2 tƒ 2
rTeƒ ƒ ,M −
2
1.540 · 109
Weƒ ƒ ,y,2
= = 7.558 · 106 mm3
20
213.6 −
2
Web (bending):
σ2 My,Ed / Weƒ ƒ ,y,2
Weƒ ƒ ,y,1 4.20 · 106
Ψ = = = =
σ1
My,Ed / Weƒ ƒ ,y,1
Weƒ ƒ ,z,2
7.558 · 106
Ψ = −0.56 > −1
h 570 42 · ϵ 42 · 0.924
= = 57.0 > =
t 10 0.67 + 0.33 · Ψ 0.67 − 0.33 · 0.56
h
= 57.0 > 79.8 (Class 3)
t
The web is at least of Class 3.
The effective section modulus Weƒ ƒ ,y for the design resistance to uni-
form bending is defined as the smallest value of the effective section
moduli at the centroid of the upper and lower flange, Weƒ ƒ ,y, and
Weƒ ƒ ,y,2 , respectively:
eƒ ƒ ,y 1.540 · 109
Weƒ ƒ ,y,1 = =
tƒ 1 + tƒ 2 10 + 20
h + − rTeƒ ƒ ,M 570 − − 213.6
2 2
Weƒ ƒ ,y,1 = 4.144 · 106 mm3
eƒ ƒ ,y 1.540 · 109
Weƒ ƒ ,y,2 = = = 7.205 · 106 mm3
zTeƒ ƒ 213.6
2. STABILITY CHECK
NRk is calculated assuming that the NRk = A · ƒy = 24778.1 · 275 = 6813977.5 N = 6813.97 kN
cross-section is subject only to stresses
due to uniform axial compression, EN My,Rk = Wp,y · ƒy = 4.144 · 106 · 275 = 1139.6 · 10−6 Nm
1993-1-5, 4.3(4).
My,Rk is calculated assuming the cross
My,Rk = 1139.6 kNm
section is subject only to bending
stresses, EN 1993-1-5, 4.3(4)
STEP 2: Reduction coefficients due to flexural buckling, χy and χz
Lcr,y = β · L = 4.00 m
v
t Aeƒ ƒ
Lcr,y A
EN 1993-1-1, Eq. 6.51 λy = ·
y λ1
v v
t y t 174700.00
y = = = 24.38 cm
A 294.00
λ1 = 93.9 · ϵ
v
u 235 v t 235
ϵ= = = 0.9244
t
ƒy 275
v
t 247.78
400 294.00
EN 1993-1-1, Eq. 6.51 λy = ·
24.38 93.9 · 0.9244
λy = 0.173
ϕ = 0.51
1
χy = s ≤ 1.0 χ, ϕ, EN 1993-1-1, Eq. 6.49
2
ϕ+ ϕ2 − λy
1
χy = p
0.51 + 0.512 − 0.1732
χy = 1.01 ≤ 1.0 → χy = 1.0
Lcr,z = β · L = 4.00 m
v
t Aeƒ ƒ
Lcr,z A
λz = · EN 1993-1-1, Eq. 6.51
z λ1
v v
t z t 153200.00
z = = = 22.82 cm
A 294.00
v
t 247.78
400 294.00
λz = · EN 1993-1-1, Eq. 6.51
22.82 93.9 · 0.9244
λz = 0.185
ϕ = 0.51
1
χz = r ≤ 1.0 χ, ϕ, EN 1993-1-1, Eq. 6.49
2
ϕ+ ϕ2 − λz
1
χz = p
0.51 + 0.512 − 0.1732
χz = 1.01 ≤ 1.0 → χz = 1.0
Cmy = 0.6 + 0.4 · ψy ≥ 0.4 Interaction factors k,j for members not
susceptible to torsional deformations,
Annex B, Table B.1
NEd NEd
kyy = Cmy · 1 + 0.6 · λy · ≤ Cmy · 1 + 0.6 ·
χy · NRk χy · NRk
5500
kyy = 1.0 · 1 + 0.6 · 0.173 ·
1.0 · 6813.97
5500
≤ 1.0 · 1 + 0.6 ·
0.173 · 6813.97
Final expression
31.5 Conclusion
In the reference example, the effective area Aeƒ ƒ is determined assuming that the cross-section is sub-
jected only to stresses due to uniform axial compression (EN 1993-1-5, 4.3(3)) Ac,eƒ ƒ = ρ · Ac . The
effective section modulus Weƒ ƒ is determined assuming the cross-section is subject to only bending
stresses (EN 1993-1-5, 4.3(3)).
By using the NEFF SIG SMIN input it is possible to define only one effective cross-section for the de-
sign and stability check, therefore the effective section modulus is determined assuming that the cross-
section is subject only to stresses due uniform axial compression. The Aeƒ ƒ as well as Weƒ ƒ ,y and
Weƒ ƒ ,z values are calculated in SOFiSTiK for the effective cross-section as shown in Fig. 31.3. This
approach checks the MOST UNFAVOURABLE case where all plates are under compression.
By using the iterative method (EN 1993-1-5, Annex E) for calculating the effective cross-section proper-
ties, the effective CS properties will be calculated for the current stress state, so it gives more realistic
and economical results as shown in table 31.2. The iterative method can be used ONLY for the THIN-
WALLED cross-sections. In Fig. 31.5 you will find the comparison between ”SIG NEFF”, ”Iterative
approach” and the reference.
1 1
2
SOF. Iterative 4
SOF. SIG NEFF
REF.
2
4
3
Figure 31.5: Comparison of the b,eƒ ƒ values
31.6 Literature
[7] EN 1993-1-1: Eurocode 3: Design of concrete structures, Part 1-1: General rules and rules for
buildings. CEN. 2005.
[8] Schneider. Bautabellen für Ingenieure. 21th. Bundesanzeiger Verlag, 2014.
[21] DIN EN 1993-1-1:2005 Eurocode 3: Design of steel structures, Part 1-1: General rules and rules
for buildings - Deutsche Fassung EN 1993-1-1:2005 + AC:2009. CEN. 2010.
[22] D. Beg et al. Design of plated structures. Ernst & Sohn and ECCS, 2010.
[23] EN 1993-1-5: Eurocode 3: Design of steel structures, Part 1-5: Plated structural elements. CEN.
2006.
Overview
For the concrete, strength class C35/45 (ƒck = 35 MP, γc = 1.5) is assumed, for the reinforcing steel,
grade B500B (ƒyk = 500 MP, Es = 205 GP, γs = 1.15, ductility class B). The factored design load
accounting for self-weight, dead load and imposed load is ed = 14.67 kN/ m2 .
Column Slab
450/450 mm h=240 mm
3
eƒ ƒ , = 6.75 m
eƒ ƒ ,y = 6.75 m
1
C B A
Figure 32.1: Model
19
ot 1 1 ot
VEd
In the reference example the reaction forces are estimated by using contributive areas, therefore the
example has been splitted into three models to show the punching for
β 1.10 1.10
β 1.40 1.40
VEd,red [kN] − −
β 1.35 1.35
32.4.1 Material
3.1 Concrete Concrete 35/45
ƒyk
ƒyd = = 435.00 N/ mm2
γs
Es = 205000 MP
Ductility class: B
d = h − c, − 0.5 · ϕ
= 240 − 30 − 10
= 200 mm
1 The sections mentioned in the margins refer to DIN EN 1992-1-1 [1] unless otherwise
specified.
dy = d − ϕ
= 200 − 20
= 180 mm
0.5 ·
0.5 · y
0.5 · y
3
0.5 · y
y
0.5 · y
2
0.6 · y
y
0.4 · y
1
C B A
Figure 32.3: Load distribution - columns
Effective depth d
d + dy
d=
2
0.2 + 0.18
d=
2
d = 0.19 m
24
19
45
2d
Perimeter 0 and 1
0 = 4 · 0.45 = 1.80 m
1 = 4.19 m
β · VEd
6.4.3, Eq. 6.38 Ed =
· d
VEd = 809 kN
BEMESS is reducing the VEd value by dead load of the slab.
Where:
1.1 ≤ β ≤ 1.8
MEd
β=1+k· 1
VEd ·
W1
R
The W value is calculated acc. W = 0
|e|d
6.4.3, Eq. 6.39 and Eq. 6.40
0.809 · 1.10
Ed = = 1.118 MN/ m2
4.19 · 0.19
Shear resistance without punching reinforcement
0 1.80
for = = 9.5 > 4
d 0.19
0.18
Rd,c = · k · (100 · ρ · ƒck )1/ 3 + 0.1 · σcp ≥ mn + 0.10 · σcp NDP, 6.4.4(1), Eq. 6.47, slabs without
γc punching reinforcement
with:
Reinforcement ratio
31.42
ρ, = = 0.0157
100 · 20
31.42
ρ,y = = 0.0175
100 · 18
0.18
Rd,c = · 2.0 · (100 · 0.0166 · 35)1/ 3
1.5
Rd,c = 0.928 MN/ m2 > mn
ot − 0
ot =
2·π
6.05 − 1.80
=
2·π
= 0.67 m ⇒≈ 3.52 · d
with:
sr = 0.75 · d
(Ed − 0.75 · Rd,c ) · 1 · d
As =
d
1.5 · · ƒyd,eƒ
sr
(1.118 − 0.75 · 0.928) · 4.19 · 0.19
As = · 104
1
1.5 · · 297
0.75
As = 5.66 cm2
Detailing of reinforcement
The spacing of link legs around a perimeter should not exceed 1.5 · d
within the first control perimeter (2 · d from loaded area), and should not
exceed 2 · d for perimeters outside the first control perimeter where that
part of the perimeter is assumed to contribute to the shear capacity.
2.40
s1 = 2.40 m ⇒ min n = =9
1.5 · 0.19
3.29
s2 = 3.29 m ⇒ min n = = 12
1.5 · 0.19
Min. punching reinforcement:
p
0.08 ƒck
As,mn = · · sr · st (NCI), 9.4.3: Eq. (9.11DE)
1.5 ƒyk
p
35
As,mn = 0.05333 · · 0.75 · 1.5 · 1.92
500
As,mn = 0.26 cm2
24
19
45
2d
2d
Edge
Critical perimeter
1 = 2.54 m
For edge and corner columns the effective perimeter is reduced based
by using the Sector Method (See Fig. 32.6). The Sector Method de-
livers the effective perimeter of the punching round cut. The ratio
/ 0 is output in % in the result list.
BOUNDARY OF SLAB
2d
s
OPENING
Control perimeter
6·d
VEd = 319 kN
β · VEd
Ed = Ed , 6.4.3, Eq. 6.38
· d
β = 1.4 NDP, 6.4.3(6),
β = 1.4 for edge columns
1.40 · 0.319
Ed =
2.54 · 0.19
Ed = 0.925 MN/ m2
with
k = 2.0
20.11
ρ = = 0.01
100 · 20
31.42
ρy = = 0.0175
100 · 18
0.18
Rd,c = · 2.0 · (100 · 0.0132 · 35)1/ 3
1.5
Rd,c = 0.860 MN/ m2 > mn
< 0.925 MN/ m2
VEd
(NCI) for 6.4.5 (4): For ot → Rd,c ot = β ·
is used according to 6.2.2 (1), CRd,c = Rd,c · d
0.15/ γc and not 0.18/ γc
0.319
ot = 1.4 ·
0.15
0.860 · · 0.19
0.18
ot = 3.28 m
ot − o
ot =
π
3.28 − 3 · 0.45
= = 0.61 m →≈ 3.21 · d
π
Punching reinforcement is required until (3.21 − 1.5) · d = 1.71 · d
sr = 0.5 · d
1 · d
As = (Ed − 0.75 · R,dc ) ·
d
1.5 · · ƒyd,eƒ
sr
0.5 · y 0.5 ·
3
0.5 · y
y
0.5 · y
2
0.6 · y
y
0.4 · y
1
C B A
Figure 32.7: Load distribution - wall
d = 190 mm
b1 = b = 350 mm < 3 · d
1
= b = 350 mm
2
< 3 · d − 0.5 · b1 = 3 · 190 − 175 = 395 mm
B
225
350
350
ed = 14.67 kN/ m2
Load distribution:
VEd = ALE · ed
Critical perimeter:
k = 2.0
ρ = 0.014
≤ 0.02
ƒcd
≤ 0.50 · = 0.023
ƒyd
0.18
Rd,c = · 2.0 · (100 · 0.0140 · 35)1/ 3
1.5
Rd,c = 0.878 MN/ m2 > mn
< 1.208 MN/ m2
ot − o
ot =
π
3.69 − 3 · 0.35
= = 0.84 m →≈ 4.42 · d
π
Punching reinforcement is required until (4.42 − 1.5) · d = 2.92 · d
sr = 0.5 · d
1 · d
As = (Ed − 0.75 · R,dc ) ·
d
1.5 · · ƒyd,eƒ
sr
32.5 Conclusion
The program searches for the single support nodes (single columns, wall ends as well as wall corners),
and performs a punching check for these points. Nodes with less than 5 kN support reaction are not
considered! Because the focus of the verification example is punching, the value ρ is overtaken from
the verification example.
It has been shown that the results are reproduced with excellent accuracy.
32.6 Literature
[1] DIN EN 1992-1-1/NA: Eurocode 2: Design of concrete structures, Part 1-1/NA: General rules and
rules for buildings - German version EN 1992-1-1:2005 (D), Nationaler Anhang Deutschland -
Stand Februar 2010. CEN. 2010.
[4] Beispiele zur Bemessung nach Eurocode 2 - Band 1: Hochbau. Ernst & Sohn. Deutschen Beton-
und Bautechnik-Verein E.V. 2011.
Overview
My
A0s N
h
As
• Design stress-strain curves for concrete and reinforcement (Section 3.1.7, 3.2.7)
• Basic assumptions for calculation of crack widths (Section 7.3.2, 7.3.3, 7.3.4)
A
ϕ h−
c B
ε2 = 0
h d
A C
E D
hc,eƒ
B ε1 5(c + ϕ/ 2)
The design stress-strain diagram for reinforcing steel considered in this example, consists of an inclined
top branch, as defined in DIN EN 1992-1-1:2004 [1] (Section 3.2.7).
10.0
.0
10
.0
10
10.0
15.19
27.00
35.44
40.50
42.19
40.50
35.44
B 500A d = 21 cm
b = 1000.0 cm
ϕs0 = 10.0 mm
k,req = 0.3 mm
Table 33.2: Results - Crack width calculation without direct method (quad element 10226)
ϕs [mm] 10 10
ϕ∗
s
[mm] 11.337 11.337
Table 33.3: Results - Crack width calculation with direct method (quad element 10226)
33.4.1 Material
Concrete: γc = 1.50 (NDP) 2.4.2.4: (1), Tab. 2.1DE: Partial
factors for materials
Steel: γs = 1.15
ƒcd = cc · ƒck / γc = 0.85 · 25/ 1.5 = 14.17 MP 3.1.6: (1)P, Eq. (3.15): cc = 0.85 con-
sidering long term effects
ƒyk = 500 MP 3.2.2: (3)P: yield strength ƒyk =
500 MP
ƒyd = ƒyk / γs = 500/ 1.15 = 434.78 MP 3.2.7: (2), Fig. 3.8
Design Load:
NEd = 0.0
405.09 · 5.5
ϕs = 10 mm = ϕ∗
s
· = ϕ∗
s
· 0.475
4(24 − 21) · 100 · 2.9
ϕs = 10 = ϕ∗
s
· 0.640 ≥ ϕ∗
s
· 0.882
ϕs 10
→ ϕ∗
s
= = = 11.337 mm
0.882 0.882
Ç
σs,tb = k · 3.48 · 106 / ϕ∗
s
(NDP) 7.3.3 Tab. 7.2DE (a)
where k = 0.3 mm the prescribed
p maximum crack width
σs,tb = 0.3 · 3.48 · 106 / 11.337 = 303.448 MP
mEd mEd
s,req = ≈
σs,tb · z σs,tb · 0.9 · d
42.11 · 100
s,req = = 7.34 cm2 / m0
30.3448 · 0.9 · 21
42.11 · 10−3
σs = = 405.09 MP
5.5 · 10−4 · 0.9 · 0.21
Effective area Ac,eƒ ƒ of concrete in tension surrounding the reinforce-
ment :
ƒct,eƒ ƒ
σs − kt · · (1 + αe · ρp,eƒ ƒ )
ρp,eƒ ƒ σs
Crack width k , 7.3.4, Eq. 7.9 ϵsm − ϵcm = ≥ 0.6 ·
Es Es
7.3.4, kt = 0.4 for long term loading, kt = 0.4
kt = 0.6 for short term loading
0
As + ξ21 · Ap
7.3.4, Eq. 7.10 ρp,eƒ ƒ =
Ac,eƒ ƒ
Es 200000
αe = = = 6.354
Ecm 31476
A0p = 0
0
As + ξ21 · Ap As
7.3.4, Eq. 7.10 ρp,eƒ ƒ = =
Ac,eƒ ƒ Ac,eƒ ƒ
5.5
ρp,eƒ ƒ = = 0.00916
600
2.56
405.09 − 0.4 · · (1 + 6.35 · 0.00916)
0.00916
Crack width k , 7.3.4, Eq. 7.9 ϵsm − ϵcm =
200000
σs 405.09
· 0.6 = · 0.6 = 0.00122
Es 200000
ϵsm − ϵcm = 0.00143 ≥ 0.00122
ϕ σs · ϕ
sr,m = ≤ (NDP) 7.3.4 (3), Eq. 7.11
3.6 · ρp,eƒ ƒ 3.6 · ƒct,eƒ ƒ
10 405.09 · 10
sr,m = ≤
3.6 · 0.00917 3.6 · 2.56
sr,m = 303.03 mm ≤ 439 mm
The steps are now iterated and the reinforcement will be increased
until k = req = 0.3 mm.
σs = 337.58 MP
ƒct,eƒ ƒ
σs − kt · · (1 + αe · ρp,eƒ ƒ )
ρp,eƒ ƒ σs
ϵsm − ϵcm = ≥ 0.6 · Crack width k , 7.3.4, Eq. 7.9
Es Es
kt = 0.4 7.3.4, kt = 0.4 for long term loading,
kt = 0.6 for short term loading
0
As + ξ21 · Ap
ρp,eƒ ƒ = 7.3.4, Eq. 7.10
Ac,eƒ ƒ
Es 200000
αe = = = 6.354 7.3.4, αe is the ratio Es / Ecm
Ecm 31476
A0p = 0
0
As + ξ21 · Ap As
ρp,eƒ ƒ = =
Ac,eƒ ƒ Ac,eƒ ƒ
6.6
ρp,eƒ ƒ = = 0.0110
600
2.56
337.58 − 0.4 · · (1 + 6.35 · 0.0110)
0.0110
7.3.4, Eq. 7.9 ϵsm − ϵcm =
200000
σs 337.58
· 0.6 = · 0.6 = 0.00101
Es 200000
ϵsm − ϵcm = 0.0011899 ≥ 0.00101
ϕ σs · ϕ
(NDP) 7.3.4 (3), Eq. 7.11 sr,m = ≤
3.6 · ρp,eƒ ƒ 3.6 · ƒct,eƒ ƒ
10 337.58 · 10
sr,m = ≤
3.6 · 0.0110 3.6 · 2.56
sr,m = 252.52 mm ≤ 366.30 mm
33.5 Conclusion
This example shows the calculation of crack widths. Various ways of reference calculations are demon-
strated, in order to compare the SOFiSTiK results to. It has been shown that the results are reproduced
with excellent accuracy.
33.6 Literature
[1] DIN EN 1992-1-1/NA: Eurocode 2: Design of concrete structures, Part 1-1/NA: General rules and
rules for buildings - German version EN 1992-1-1:2005 (D), Nationaler Anhang Deutschland -
Stand Februar 2010. CEN. 2010.
[2] F. Fingerloos, J. Hegger, and K. Zilch. DIN EN 1992-1-1 Bemessung und Konstruktion von
Stahlbeton- und Spannbetontragwerken - Teil 1-1: Allgemeine Bemessungsregeln und Regeln
für den Hochbau. BVPI, DBV, ISB, VBI. Ernst & Sohn, Beuth, 2012.
Overview
The main goal of this benchmark is to verify and compare the SOFiSTiK results with the reference
example Beispiele zur Bemessung nach Eurocode 2 - Band 1: Hochbau [4, Example 10]. In this example
the Ultimate Limit State (ULS) of the prefabricated column is verified.
D1
1.9 m
A
Detail A β = 2.1
6.2 m
225 225
450
β = 2.1
The design stress-strain diagram for reinforcing steel considered in this example, consists of an inclined
top branch, as defined in DIN EN 1992-1-1:2004 [1] (Section 3.2.7).
1 See conclusion (Section 34.5) for more details why the results are different compared to the reference example.
Dead loading
Gk,1 = 400 kN
Gk, = 431 kN
P
Variable loading
Snow: Qk,s = 68 kN
34.4.3 Materials
Concrete
Class C 30/37
3.1.6: (1)P, Eq. (3.15): cc = 0.85 con- ƒcd = cc · ƒck / γc = 0.85 · 30/ 1.5 = 17.0 MP
sidering long term effects
1 The tools used in the design process are based on steel stress-strain diagrams, as
Steel
Class B 500B
γs = 1.15
0,y = βy · co
y,0 p
λy = = 13.0 · 12/ 0.45 = 100 λy , slenderness
y
0,z = βz · co
z,0 p
λz = = 6.2 · 12/ 0.40 = 54 λz , slenderness
z
Slenderness criterion for isolated member
NEd 0.684
n= = = 0.224 n is the relative normal force, (NDP)
Ac · ƒcd 0.40 · 0.45 · 17.0 5.8.3.1 (1)
y
38 38
200
z
400
200
225 225
450
16 16
λm = p = = 34 for n = 0.224 < 0.41
n 0.224
Because λy = 100 and λz are > λm = 34 second order effects
5.8.3.1(2) should be taken into account in both directions (y and z).
34.4.5 Imperfections
For isolated members, the effect of imperfection may be taken into ac-
count:
0
e is the eccentricity acc. 5.2(7), a) Eq. e = θ ·
(5.2) 2
θ is the inclination value acc. Eq. (5.1) θ = θ0 · αh · αm
1
θ0 is the basic inclination value, Eq. θ0 =
(5.1) 200
2 2
αh is the reduction factor for length of αh = p = p = 0.803 < 1.0
the height 6.2
0,y = 13.00 m
1
θ = 0.803 · 0.79 · = 0.00317 = 1/ 315
200
13000
e,z = = 21 mm
2 · 315
Bending about the z-axis:
0,z = 6.20 m
1
θ = 0.803 · 1.0 · = 0.00402 ≈ 1/ 250
200
6200
e,y = = 13 mm
2 · 250
K1 = 1 for λ > 35
ƒcd · Ac
n = 1 + ω = 1 +
ƒyd · As
e1 [mm] 21 21 21
1 20
e2 = K1 · · [mm] 199 199 199
r 10
etot = e0 + e + e2 [mm] 329 430 378
ƒcd
As,tot = ωtot · b · h · [cm2 ] 13.40 12.70 16.20
ƒyd
M0Eqp
ϕ(∞, t0 ), 3.1.4, Fig. 3.1 a) ϕeƒ = ϕ(∞, t0 ) ·
M0Ed
Kƒ φ = 1 + β · φeƒ ≥ 1.0
ƒck λ
β = 0.35 + − ≥0
200 150
300 100
β = 0.35 + − = −0.17 < 0
200 150
Kφ = 1 → Creep and shrinkage is neglected!
Buckling length:
0,z = 6.20 m
Slenderness:
λz = 54 < λm = 86
co 6.20
= = 15.50 > 2.5
h 0.40
The column should be investigated for the second order effects.
13 6200 2 · 13
= 1.14 · (1 − 2 · ) − 0.02 · ≤1−
400 400 400
= 0.756 ≤ 0.94
0.70 · 30
ƒcd,p = = 14.0 N/ mm2
1.5
NRd = 0.45 · 0.40 · 14.0 · 0.756 = 1.9 sMN > |NEd | = 0.684 MN
NEd = 431 kN
d = 0.41 m
0.15
VRd,c = · k · (100 · ρ · ƒck ) 1/ 3 + 0.12 · σcp · b · d >
γc
VRd,c,mn Shear design, 6.2.2, Eq. (6.2)
v
t 200
k =1+ ≤ 2.0
d
v
t 200
k =1+ = 1.70
410
As 8.04
ρ = = = 0.0049 < 0.02
b · d 40 · 41
ƒck = 30 N/ mm2
NEd 0.431
σcp = = = 2.39 MN/ m2 C30/37 → 0.2 · ƒcd = 3.4 MP
Ac 0.40 · 0.45
0.15
VRd,c = · 1.70 · (100 · 0.49 · 30)1/ 3 + 0.12 · 2.39 · 0.40 ·
1.5
0.41 >= 0.115 MN
34.5 Conclusion
This example shows the calculation of column design by using the norminal curvature method. Various
ways of reference calculations are demonstrated, in order to compare the reference results to SOFiSTiK.
The main reason for minor deviations (See Tab. 34.3 in Section 34.3) is that in the reference benchmark,
the 2nd order effects and the reinforcement are approximated by using analytical formulas. In SOFiSTiK,
the Kr value is iterated until Kr < 1 · 10−4 , if this condition is met, then the iteration is stopped. With a
hand calculation this is very difficult to achieve.
START
New As
As · ƒyd
n = 1 +
Ac · ƒcd
NEd
n=
Ac · ƒcd
n + n
Kr =
n − nb
No
|Kr,ne − Kr,od | < 1 · 10−4
Yes
END
2. Find the critical deflection e0 , which depends on the current loadcase combination. Consider
different end-moments according EN 1992-1-1:2004, 5.8.8.2 (2), Eq. 5.32
3. Start iteration of design moment MEd (Theory II. Order). The iteration will stop, when the coefficient
Kr achieved convergence.
4. Additional cross section design on different locations (without imperfection and theory II.O.)
34.6 Literature
[1] DIN EN 1992-1-1/NA: Eurocode 2: Design of concrete structures, Part 1-1/NA: General rules and
rules for buildings - German version EN 1992-1-1:2005 (D), Nationaler Anhang Deutschland -
Stand Februar 2010. CEN. 2010.
[2] F. Fingerloos, J. Hegger, and K. Zilch. DIN EN 1992-1-1 Bemessung und Konstruktion von
Stahlbeton- und Spannbetontragwerken - Teil 1-1: Allgemeine Bemessungsregeln und Regeln
für den Hochbau. BVPI, DBV, ISB, VBI. Ernst & Sohn, Beuth, 2012.
[4] Beispiele zur Bemessung nach Eurocode 2 - Band 1: Hochbau. Ernst & Sohn. Deutschen Beton-
und Bautechnik-Verein E.V. 2011.
Overview
The main goal of this benchmark is to verify and compare the SOFiSTiK results with the ECCS reference
example Beg et al. [22, Example 2.4-3].
• EBPlate (2007)
• EN 1993-1-5 rules
Figure 35.2: SOFiSTiK FEM Model - Steel plate with 2 horizontal stiffeners
hs
σc
b1
ts
b b2 = b1
b3 = b1
a = 1800 mm, b = 1800 mm, b1 = 600 mm, hs = 100 mm, ts = 10 mm.
S 355, ƒy = 355 N/ mm2 , ϵ = 0.81
• Calculating the critical plate buckling stress (EN 1993-1-5 [23], Annex A.2)
To calculate the critical plate buckling stress the loading, σc = 1.0 N/ mm2 is used. The critical elastic
stress will be calculated by multiplying the minimum eigenvalue with the unity stress σc = 1.0.
σcr,p = αcr · σc
ϵ = 0.81 b1 = 600 mm
hs = 100 mm
ts = 10 mm
t = 12 mm
N
σcr 275.793 274.170 268.721 2682 290
mm2
1 EBPlate, V2.01
2 The results were overtaken from the ECCs reference example Beg et al. [22, Example 2.4-3]
1.0923
1.0798
1.0731
1.02 1.02 1.02 1.03
93
1.0664
0.9
1.0598
1.0531
1.0465
1.0398
1.0331
1.0265
1.0198 4
1.0131 1.0
1.01 1.01 1.01 1.01
1.0065 1.07 1.01 1.07
0.9998 1.01
0.9931 1.08 1.08
0.9865
0.9798
0.9731
0.9665 0.991
1.06
0.9598 1.06 0.991
0.9465
0.9398
0.9331
0.9265 1.08 1.08
0.9198 1.01
0.9131 1.08 1.09
0.9065 1.01
04
1.
0.8998
0.8932
0.8865
0.8798
0.8732
0.8665
0.8598 1.01
0.8532
0.8465 1.03 0.993 1.03
0.8398
0.8332
0.8265
0.8257
Figure 35.3: SOFiSTiK FEM Model - Stress distribution along the plate for σc = 1.0 N/ mm2
35.4.1 Klöppel
σcr,p is given with the following equation:
σcr,p = kσ,p · σE
Where:
π2 · E · t2
σE =
12 · (1 − ν 2 ) · b2
π 2 · 210000 · 122
σE = = 8.436 N/ mm2
12 · (1 − 0.32 ) · 18002
kσ,p is the elastic critical plate buckling coefficient according to Klöppel.
1800
α= = = 1.0
b 1800
As bs · t 100 · 10
δ= = = = 0.05
b·t b·t 1800 · 12
(s + As · e2 ) · 12 · (1 − ν 2 )
γ=
b · t3
ts
b3s · + bs · ts · e2 · 12 · (1 − ν 2 )
12
γ=
b · t3
10
+
1003 · 100 · 10 · 502 · 12 · 1 − 0.32
12
γ= = 11.70
1800 · 123
Note that parameter α and δ above are not the same as in EN 1993-1-
5, Annex A.1, where the procedure for plates stiffened with more than
two stiffeners is given.
kσ = 32.5
36
σ1 100
γ, δ b/ 3
30 15
γ, δ b/ 3 80
b/ 3
10 ψ · σ1 = σ1 = · b 60
50
20 6 40 30
kσ 25 20
4 15
10
10 2 6
24
1 1
0 0
0
0.7 1.0 1.5 2.0 2.5 3.0 3.5 3.8
α
35.4.2 EBPlate
The usual procedure (calculation of buckling modes) for the calculation
of critical stresses is presented by using EBPlate.
b1 b2
As,1
s,1
b
b2 b3
As,2
s,2
b
Amped,1
mped,2
bmped,1 bmped,2
bmped
Stiffeners I and II
b1 / 2 ts b2 / 2
e2
e1
bs
Figure 35.7
595 + 590
As,1 = + 10 · 12 + 100 · 10 = 8230 mm2
2
b3s · ts (b1 + b2 ) · 0.5 + ts,1 · t 3
s,1 = + + bs · ts · e21 +
12 12
(b1 + b2 ) · 0.5 + ts,1 · t · e22
c = 2998 mm
π 2 · E · s,1
σcr,s =
As,1 · 2
E · t 3 · b · 2
+
4 · π 2 · (1 − ν 2 ) · As,1 · b21 · b22
In case of a stress gradient over the plate width, the critical plate buck-
ling stress should be properly interpolated from the position of the stiff-
ener to the most stressed edge of the plate. In this case no stress gra-
dient over the depth of the plate is present. Therefore, the critical plate
buckling stress is equal to the critical stress calculated for the buckling
of the stiffener on the elastic support:
σcr,p = σcr,p
= σcr,s = 322 N/ mm2
Lumped stiffener
mped = s
+ = 3.675 · 106 + 3.675 · 106 = 7.35 · 106 mm4
s
v
mped · b2mped,1 · b2mped,2
u
u
4
c = 4.33 ·
t
t 3 · bmped
v
6 2 2
4 7.35 · 10 · 900 · 900
u
c = 4.33 ·
t
123 · 1800
c = 4832 mm
π 2 · E · mped
σcr,mped =
Amped · 2
E · t 3 · bmped · 2
+
4 · π 2 · (1 − ν 2 ) · Amped · b2mped,1 · b2mped,2
Finally we have:
35.4.4 ABAQUS
The results from ABAQUS have be overtaken from the ECCS reference
example Beg et al. [22, Example 2.4-3], see Table 35.2 for more details.
35.5 Conclusion
The critical plate buckling stress was calculated by using:
• SOFiSTiK
• Klöppel diagrams
• EBPlate
• ABAQUS
• EN 1993-1-5 rules
The methods used in the calculation give very similar results. The advantage of SOFiSTiK compared to
Eurocode formulas and other tools (that are not using FEM) is that the stiffeners can be added customly.
The cases with variable height of plates can be analysed as well (haunches). In conclusion, it has been
shown that the SOFISTiK results are reproduced with excellent accuracy.
(a) 1st buckling mode: σcr,p = 275.79 N/ mm2 (b) 2nd buckling mode: σcr,p = 324.22 N/ mm2
35.6 Literature
[7] EN 1993-1-1: Eurocode 3: Design of concrete structures, Part 1-1: General rules and rules for
buildings. CEN. 2005.
[22] D. Beg et al. Design of plated structures. Ernst & Sohn and ECCS, 2010.
[23] EN 1993-1-5: Eurocode 3: Design of steel structures, Part 1-5: Plated structural elements. CEN.
2006.
Overview
h My
Np
zp
Ap
Ap = 28.5 cm2
Material:
εcs (t, ts ) = εcbs (t) + εcds (t, ts ) 5.1.9.4.4: Eq. 5.1-75: εcs (t, ts ) total
shrinkage or swelling strains
Calculating the basic shrinkage:
εcbs (t) = εcbs0 (ƒcm ) · βbs (t) 5.1.9.4.4: Eq. 5.1-76: εcbs (t) is the ba-
sic shrinkage
2.5
0.1 · ƒcm
ε(ƒcm ) = −αbs · · 10−6 5.1.9.4.4: Eq. 5.1-78: εcds0 (ƒcm ) is the
6 + 0.1 · ƒcm notional shrinkage coefficient
p p
βbs (t) = 1 − ep(−0.2 · 36500) − 1 − ep(−0.2 · 28)
βbs (t) = 1 − ep(−38.2099) − (1 − ep(−1.0583))
βbs (t) = 0.347
5.1.9.4.4: Eq. 5.1-80 εcds0 = [(220 + 110 · αds1 ) · ep(−αds2 · ƒcm )] · 10−6
See table 5.1-12 Coefficients (αds ) are depending on the type of cement.
RH 3
−1.55 · 1 −
, ƒ or 40 ≤ RH < 99% · βs1
5.1.9.4.4: Eq. 5.1-81 βRH = 100
ƒ or RH ≥ 99% · βs1
0.25,
0.1
35
βs1 = ≤ 1.0
ƒcm
35 0.1
5.1.9.4.4: Eq. 5.1-83 βs1 = = 0.9796 ≤ 1.0
43
99% · βs1 = 99 · 0.9796 = 96.98
80 3
βRH = −1.55 · 1 − = −0.7564
100
SOFiSTiK accounts not only for the age at start of drying ts but also
for the influence of the age of prestressing, so the time development
function reads:
v v
u t − ts u t0 − ts
βds (t − ts ) = −
t t
5.1.9.4.4: Eq. 5.1-82
0.035 · h2 + (t − ts ) 0.035 · h2 + (t0 − ts )
v v
u
36500 u
28
βds (t − ts ) =
t t
−
0.035 · 5002 + 36500 0.035 · 5002 + 28
βds (t − ts ) = 0.8981 − 0.05669 = 0.8416
The drying shrinkage is calculated:
εcds (t, ts ) = εcds0 (ƒcm ) · βRH · βds (t − ts )
n
X
4000
tT = Δt · ep 13.65 − 5.1.10.2: Eq. 5.1-85; tT the adjusted
=1
273 + T(Δt ) concrete age, during the effect of ele-
vated or reduced temperatures on the
n
4000
X maturity of concrete
t0,T = 28 · ep 13.65 − = 27.947 days
=1
273 + 20
RH 80
1− 1−
100 100
5.1.9.4.3(b): Eq. 5.1-69 β(RH) = v = v = 0.251
h 500
u u
3
t 3
t
0.1 · 0.1 ·
100 100
1 1
5.1.9.4.3: Eq.5.1-70 βdc (t0 ) = = = 0.4886
0.2
0.1 + t0,dj 0.1 + 27.947
γ(t0 )
t − t0
5.1.9.4.3(b): Eq. 5.1-71a; βdc (t, t0 ) the βdc (t, t0 ) =
development of drying creep with time βh + (t − t0 )
with:
1 1
5.1.9.4.3: Eq. 5.1-71b γ(t0 ) = == 0.3376
3.5 2.962
2+ p
t0,dj
v v
u 35
t 35
u
αƒcm = = = 0.9021
t
5.1.9.4.3: Eq. 5.1-71d
ƒcm 43
5.1.9.4.3: Eq. 5.1-71c βh = 1.5 · h + 250 · αƒcm ≤ 1500 · αƒcm
βh = 1.5 · 500 + 250 · 0.9021 = 975.548 ≤ 1352.15
0.3376
36500 − 28
5.1.9.4.3: Eq. 5.1-71a βdc (t, t0 ) = = 0, 9911
975.548 + (36500 − 28)
According to Model Code 2010 [16], the creep value is related to the
tangent Youngs modulus Ec , where Ec being defined as 1.05 · Ecm . To
account for this, SOFiSTiK adopts this scaling for the computed creep
coefficient (in SOFiSTiK, all computations are consistently based on
Ecm .
ϕ(t, t0 ) = 1.64/ 1.05 = 1.56
36.5 Conclusion
This example shows the calculation of the creep and shrinkage using fib Model Code 2010 [16]. It has
been shown that the results are in very good agreement with the reference solution.
36.6 Literature
[16] fib Model Code 2010. fib Model Code for Concrete Structures 2010. International Federation for
Structural Concrete (fib). 2010.
Overview
h My
Np
zp
Ap
Ap = 28.5 cm2
ϕ0 1.57 - 1.566
Material:
The total shrinkage or swelling strans εcs (t, ts ) may be calculated from
εcs (t, ts ) = εcs0 · βs (t − ts ) 2.1.6.4.4: Eq. 2.1-74; εcs (t, ts ) is the
total or swelling strain
Calculating the notional shrinkage:
3 3
RH 80
βsRH = 1 − =1− = 0.488 2.1.6.4.4: Eq.2.1-78; RH0 = 100 %
RH0 100
βRH = −1.55 · 0.488 = −0.7564
βs = βs (t − ts ) − βs (t0 − ts )
0.5 0.5
36500 28
βs = −
350 · 52 + 36500 350 · 52 + 28
βs = 0.8981 − 0.05647 = 0.8416
n
X
4000
t0,T = 28 · ep 13.65 − = 27.947 days
=1
273 + 20/ 1
α
9
2.1.6.4.3(c): Eq.2.1-71; the effect of t0,dj = t0,T · 1.2
+ 1 ≥ 0.5 dys
type of cement on the creep coefficient 2 + (t0,T / t1,T )
of concrete may be taken into account
by using the modified age at loading 0
9
t0,dj ; α = 0 for cement class N
t0,dj = 27.947 · +1 = 27.947 ≥ 0.5 dys
2 + 27.9471.2
1 1
2.1.6.4.3(b): Eq. 2.1-68 β(t0 ) = = = 0.48862
0.1 + (t0 / t1 )0.2 0.1 + (27.947/ 1)0.2
βH = 1359.702 ≤ 1500
0.3
(36500 − 28)/ 1
βc (t − t0 ) = = 0.989
1359.702 + (36500 − 28)/ 1
ϕ0 = 1.254 · 2.556 · 0.48862 = 1.566
The creep value is related to the tangent Youngs modulus, where the
tangent modulus being defined as 1.05 · Ecm . To account for this,
SOFiSTiK adopts this scaling for the computed creep coefficient (in
SOFiSTiK, all computations are consistently based on the secant mod-
ulus of elasticity).
1.5489
ϕ(t, t0 ) = = 1.47
1.05
37.5 Conclusion
This example shows the calculation of the creep and shrinkage using Model Code 1990 [24]. It has been
shown that the results are in very good agreement with the reference solution.
37.6 Literature
[24] CEB-FIP Model Code 1990. Model Code for for Concrete Structures 1990. Euro-International
Concrete Committee. 1991.
Overview
For the concrete, strength class C30/37 (ƒck = 30 MP, γc = 1.5) is assumed, for the reinforcing steel,
grade B500B (ƒyk = 500 MP, Es = 205 GP, γs = 1.15, ductility class B). The factored design load
accounting for self-weight, dead load and imposed load is qd = 15.6kN/ m2 .
0,39
26
2,38
280
column
0,37
26 slab
3
600 3
600
0,1326
C1 C4
0,13
2,42
2,16
Opening
0,13
2,67
2,8
560
534
483
26 0,13
5,73
C2 C5 0,9
180
1146
2,67
0,13
2,8
2,42
2,16
560
534
483
26 0,13
0,13
C3 C6
26 0,13
2574
,87 0,13 2,87
574 0,8
160 0,77
154
7,57
1513
C30/ 37 h = 26 cm qd = 15.6kN/ m2
B500B d1 = 4.0 cm
d = 21.0 cm
In the reference example the reaction forces are estimated by using contributive areas. The results are
given in Table 38.2.
1.681 m 1.681 m
ke 0.9 0.9
ψ 1.31 % 1.325 %
ψy 1.41 % 1.424 %
kr 1.03 1.02786
STEP 1: Material
ηƒ c · ηt · ƒck 1 · 1 · 30
The dimensioning value of the concrete ƒcd = = = 20 N/ mm2
compressive strength; 2.3.2.3; Eq. (2) γc 1.5
p
(ƒck )
p
0.3 · ηt · 0.3 · 1.0 · 30
The dimensioning value of shear stress τcd = = = 1.095 N/ mm2
limit; 2.3.2.4; Eq.(3) γc 1.5
Dm = 32mm
STEP 2: Reinforcement
Es = 205000 MP
Ductility class: B
d = h − d1
1.6
= 26 − 4.8 +
2
= 20.4 cm
Inner:
= 2 · + 2 · b + d · π
= 2 · 26 + 2 · 26 + 20.4 · 3.14
= 168.1 cm
Vd = 686.1 kN
1 The sections mentioned in the margins refer to SIA 262:2013 [25] unless otherwise
specified.
Bemess takes into account the min. value of ke . For inner columns ke = 0.9, SIA
262:2013, 4.3.6.2.5
1
ke = mn 0.9; 1 + (38.1)
e
1+
b
Md 5.46 kNm
e =
− Δe =
= 7.95 mm
Vd 686.1 kN
v v
t4 t4
u u
b = · Ac = · 0.2063 = 0.5129 m = 512.9 mm
π π
1
= mn 0.9,
e
1+
b
1
= mn
0.9,
7.95
1+
512.9
= mn (0.9, 0.98)
= 0.9
STEP 5: Rotations
The distances rs, and rs,y are calculated from the results of the flexural
analysis, one can obtain the distances between the center of the column
rs = 1.166 m
rsy = 1.248 m
bs = 1.5 ·
p
rs · rsy ≤ mn
p
= 1.5 · 1.248 · 1.166
= 1.8094 m
3/ 2
rsy ƒsd msdy
ψy = 1.2 · · ·
d Es mRd
3/ 2
1.248 434.78 105.85
= 1.2 · · ·
0.204 205000 112.306
= 1.424 %
1
kr =
0.45 + 0.18 · ψ · d · kg
1
=
0.45 + 0.18 · ψ · d · kg
1
=
0.45 + 0.18 · 0.01424 · 204 · 1
= 1.02786
In BEMESS there isn’t any option to set the Dm value and the Dm
is strictly defined:
VRd,c = kr · τcd · d ·
= 1.02786 · 1.095 · 0.204 · 1.5129
= 0.34736 MN = 347.36 KN
Punching strength with shear reinforcement: 4.3.6.5; Eq. (65) and Eq. (66)
Vd,s
snβ · As =
ke · σsd
343.05
=
0.9 · 43.478
= 8.766 cm2
Es · ψ ƒbd d
σsd = · 1+ · ≤ ƒsd
6 ƒsd ϕs
205000 · 0.01424 2.703 204
= · 1+ · ≤ 434.78
6 434.78 16
= 487.00 · (1 + 0.006216 · 12.75) ≤ 434.78
= 525.60 ≤ 434.78
= 434.78 N/ mm2
0.6861
ot = = 3.7170 m
1.02786 · 1.095 · 0.164
ot − 4 · bc 3.717 − 4 · 0.26
rot = = = 0.426 m
2·π 2·π
38.5 Conclusion
This example shows the calculation of punching of flat slabs and it has been shown that the results are
reproduced with excellent accuracy.
38.6 Literature
[16] fib Model Code 2010. fib Model Code for Concrete Structures 2010. International Federation for
Structural Concrete (fib). 2010.
[25] SIA 262:2013 Concrete Structures. Swiss Society of Engineers and Architects. 2013.
[26] S. Lips, A. Muttoni, and M.F. Ruiz. “Punching of flat slabs: Design Example”. In: fib Model Code
2010 (2011).
[27] A. Muttoni et al. Background to fib Model Code 2010 shear provisions - part II: punching shear.
Tech. rep. fib, 2013.