Basic Reliability Formulation Handbook2
Basic Reliability Formulation Handbook2
HANDBOOK 2
RELIABILITY BACKGROUNDS
Probability density E(x), R(x)
Load effect E
0,06
Resistance R
0,04
dx
0,02
0,00
40
60
x x+dx
100
120
140
Random variable X
HANDBOOK 2
RELIABILITY BACKGROUNDS
PARTNERSHIP:
The Klokner Institute of the Czech Technical University in Prague (KI CTU),
convener, Prof. Milan Holick, Czech Republic
The Czech Chamber of Certified Engineers and Technicians Engaged in Construction
(CKAIT), Prof. Alois Materna, Czech Republic,
The Institute of Steel Constructions of the University of Technology Aachen (RWTH),
Prof. Gerhard Sedlacek, Germany
The Spanish Organisation for Scientific Research (IET), Spain, Dr. Angel Arteaga
The University of Pisa (UOP), Prof. Luca Sanpaolesi, Italy
The Netherlands Organisation for Applied Scientific Research (TNO), Prof. Ton
Vrouwenvelder, The Netherlands
The Institute for Metal Constructions (IMK), Dr. Igor Kovse, Slovenia
The Building Research Establishment (BRE), Prof. Haig Gulvanessian, The United
Kingdom
Prague 10.2005
Leonardo da Vinci Pilot Project CZ/02/B/F/PP-134007
II
II-1
II-1
II-1
II-1
II-1
II-1
II-3
II-6
II-8
II-10
II-11
II-11
4
IV
II-12
II-14
II-16
III-1
III-1
III-1
III-1
III-2
III-3
III-4
III-5
III-6
III-7
III-9
III-10
III-10
III-11
III-12
III-13
III-15
IV-1
IV-1
IV-1
IV-1
IV-2
IV-4
IV-4
IV-6
IV-8
IV-9
IV-9
IV-14
IV-15
ATTACHMENTS
1. EXCEL worksheet from the workbook dast.xls
2. Source file dast.c to the program dast.exe
3. Input file dast.i1 for the program dast.exe
4. Output file dast.o1 produced by dast.exe from dast.o1
5. Input file dast.i2 for the program dast.exe
6. Output file dast.o2 produced by dast.exe from dast.o2
7. EXCEL chart showing the data from dast.o2
IV-16
IV-17
IV-18
IV-21
IV-22
IV-23
IV-24
IV-24
V-1
INTRODUCTION
1.1
Background documents
1.2
General principles
2
GENERAL FRAMEWORK OF ASSESSMENT
2.1
Reasons for assessment
2.2
Common rules
2.3
General procedure
3
INVESTIGATION
3.1
Purpose
3.2
Statement
4
BASIC VARIABLES
4.1
General
4.2
Characteristic values
5
EVALUATION OF INSPECTION RESULTS
5.1
Updating in general
5.2
Updating of probability distribution
5.3
Updating of failure probability
5.4
Updating of characteristic and design values
6
STRUCTURAL ANALYSIS
7
VERIFICATION
8
ASSESSMENT IN THE CASE OF DAMAGE
9
FINAL REPORT AND DECISION
10
CONCLUDING REMARKS
REFERENCES
APPENDIX A - GENERAL FLOW OF ASSESSMENT OF EXISTING
STRUCTURES
APPENDIX B - PROBABILITY UPDATING
APPENDIX C - BAYESIAN METHOD FOR FRACTILE ESTIMATION
ATTACHMENTS
1. MATHCAD sheet Update.mcd
2. MATHCAD sheet BayesFract.mcd
VI
V-1
V-1
V-1
V-2
V-2
V-2
V-2
V-3
V-3
V-3
V-4
V-4
V-4
V-5
V-5
V-5
V-6
V-6
V-7
V-7
V-8
V-9
V-9
V-10
V-11
V-12
V-14
V-16
V-17
V-18
VI-1
VI-1
VI-1
VI-1
VI-2
VI-3
VI-4
VI-4
VI-5
VI-5
VI-8
VI-10
VI-10
VI-11
VI-11
VI-12
A-1
A-1
A-1
A-1
A-1
A-1
A-2
A-3
A-4
A-6
A-6
A-6
A-9
A-11
A-13
A-17
A-17
A-17
A-21
A-22
A-23
A-26
A-28
A-30
A-31
A-32
A-33
A-34
A-35
A-37
A-39
B-1
B-1
B-1
B-1
B-1
B-4
B-4
B-6
B-7
B-9
B-9
4.2
Basic variables
4.3
Tail sensitivity problem
5
AN EXAMPLE OF REINFORCED CONCRETE SLAB
5.1
General
5.2
A reinforced concrete slab
5.3
Design of a slab
5.4
Reliability consideration
6
ASSESMENT OF THE FAILURE PROBABILITY IN
GENERAL CASE
6.1
General
6.2
First and Second Order Reliability Methods (FORM and SORM)
7
SYSTEM RELIABILITY
7.1
General
7.2
Parallel systems
7.3
Series systems
8
CONCLUDING REMARKS
REFERENCES
B-15
B-15
B-15
B-19
B-19
B-20
B-21
B-21
B-21
ATTACHMENTS
1. MATHCAD sheet GammaRGQ.mcd
2. MATHCAD sheet Prindex.mcd
3. MATHCAD sheet RCBeam.mcd
4.M A TH CA D sheetRelRCB.m cd
5. EXCEL sheet RCBeam.xls
6. MATHEMATICA notebook Fit_distribution.nb
7. MATHEMATICA notebook FORM.nb
8. MATHEMATICA package Level2.m
9. MATHCAD sheet FORM2.mcd
10. MATHCAD sheet FORM7.mcd
11. EXCEL sheet FORM7.xls
12. MATLAB package FORM7.m
13. MATLAB function Lndens (x,mu,sigma,sk)
14. MATLAB function Lndist (x,mu,sigma,sk)
15. MATLAB function Ndens (x, mu,sigma)
16. MATLAB function "Ndinv (p)"
17. MATLAB function "Ndist (x, mu,sigma)"
B-22
B-24
B-26
B-27
B-29
B-32
B-34
B-40
B-41
B-44
B-45
B-47
B-48
B-50
B-52
B-53
B-54
B-55
B-10
B-11
B-11
B-11
B-12
B-12
B-14
C-1
C-1
C-1
C-1
C-2
C-4
C-5
C-5
C-6
C-7
C-8
5.2
Two variable actions
6
CONCLUDING REMARKS
REFERENCES
APPENDIX A - DIRECT COMPARISON OF LOAD EFFECTS
APPENDIX B - EFFECT OF THE RESISTANCE VARIABILITY
APPENDIX C NOTATION
ATTACHMENTS
1. MATLAB function generic(gR,wr,k)
2. MATLAB function "Action3(mr,sr,skr,Rd,k)
3. MATLAB function Lnpf(mr,sr,skr,me,se,ske)
4. MATLAB function Lndens(x,mu,sigma,sk)
5. MATLAB function Lndist(x,mu,sigma,sk)
6. MATLAB function Ndinv(p)
7. MATHCAD sheet Generic.mcd
8. MATHCAD sheet Load effect.mcd
C-8
C-8
C-10
C-13
C-13
C-15
C-19
C-20
C-21
C-23
C-24
C-29
C-31
C-32
C-33
C-34
C-39
FOREWORD
The Leonardo da Vinci Pilot Project CZ/02/B/F/PP-134007, Development of Skills
Facilitating Implementation of Structural Eurocodes addresses the urgent need to implement
the new system of European documents related to design of construction works and products.
These documents, called Eurocodes, are systematically based on recently developed Council
Directive 89/106/EEC The Construction Products Directive and its Interpretative Documents
ID1 and ID2. Implementation of Eurocodes in each Member State is a demanding task as each
country has its own long-term tradition in design and construction.
The project should enable an effective implementation and application of the new
methods for designing and verification of buildings and civil engineering works in all the
partner countries (CZ, DE, ES, IT, NL, SI, UK) and in other Member States. The need to
explain and effectively use the latest principles specified in European standards is apparent
from various enterprises, undertakings and public national authorities involved in construction
industry and also from universities and colleges. Training materials, manuals and software
programmes for education are urgently required.
The submitted Handbook 2 is one of 5 upcoming handbooks intended to provide
required manuals and software products for training, education and effective implementation
of Eurocodes:
Handbook 1: Basis of Structural Design
Handbook 2: Reliability Backgrounds
Handbook 3: Load Effect for Buildings
Handbook 4: Load Effect for Bridges
Handbook 5: Design of Buildings for Fire Situation
It is expected that the Handbooks will address the following intents in further
harmonisation of European construction industry:
- reliability improvement and unification of the process of design;
- development of the single market for products and for construction services;
- new opportunities for the trained primary target groups in the labour market.
The Handbook 2 is focused on the basis of structural reliability and risk engineering
related to Eurocodes. The following topics are treated in particular:
- basic concepts of structural reliability;
- elementary methods of the reliability theory;
- reliability differentiation and design working life;
- design assisted by testing;
- assessment of existing structures;
- basis of risk assessment.
Annex A to the Handbook 2 provides a review of "Basic Statistical Concepts and
Techniques" frequently used in the text. Annex B provides an extension of elementary
methods of structural reliability and annex C describes calibration procedures that may be
used for specification of reliability elements. The Handbook 2 is written in a user-friendly way
employing only basic mathematical tools. Attached software products accompanying a number
of examples enable applications of general rules in practice.
A wide range of potential users of the Handbooks and other training materials includes
practising engineers, designers, technicians, experts of public authorities, young people - high
school and university students. The target groups come from all territorial regions of the
partner countries. However, the dissemination of the project results is foreseen to be spread
into all Member States of CEN and other interested countries.
Prague 05/2005
10
Summary
Uncertainties affecting structural performance can never be entirely eliminated and
must be taken into account when designing any construction work. Various design methods
and operational techniques for verification of structural reliability have been developed and
worldwide accepted in the past. The most advanced operational method of partial factors is
based on probabilistic concepts of structural reliability and available experience. General
principles of structural reliability can be used to specify and further calibrate partial factors
and other reliability elements. Moreover, developed calculation procedures and convenient
software products can be used directly for verification of structural reliability using
probabilistic concepts and available experimental data.
1
INTRODUCTION
1.1
Background materials
Basic concepts of structural reliability are codified in a number of national standards,
in the new European document EN 1990 [1] and the International Standard ISO 2394 [2].
Additional information may be found in the background document developed by JCSS [3]
and in recently published handbook to EN 1990 [4]. Guidance for application of probabilistic
methods of structural reliability may be found in working materials provided by JCSS [5] and
in relevant literature listed in [4 and [5]. Elementary methods of the theory of reliability are
described in Chapter II and III in this Handbook 2.
1.2
General principles
General principles of structural reliability are described in both the international
documents EN 1990 [1] and ISO 2394 [2]. Basic requirements on structures are specified in
Section 2 of EN 1990 [1]: a structure shall be designed and executed in such a way that it will,
during its intended life, with appropriate degrees of reliability and in an economic way
-
sustain all actions and influences likely to occur during execution and use;
remain fit for the use for which it is required.
It should be noted that two aspects are explicitly mentioned: reliability and economy
(see also Handbook 1). However, in this Handbook 2 we shall be primarily concern with
reliability of structures, which include
-
structural resistance;
serviceability;
durability.
I-1
UNCERTAINTIES
2.1
Classification of uncertainties
It is well recognised that construction works are complicated technical systems
suffering from a number of significant uncertainties in all stages of execution and use.
Depending on the nature of a structure, environmental conditions and applied actions, various
types of uncertainties become more significant than the others. The following types of
uncertainties can be identified in general:
- natural randomness of actions, material properties and geometric data;
- statistical uncertainties due to a limited size of available data;
- uncertainties of the resistance and load effect models due to simplifications of actual
conditions;
- vagueness due to inaccurate definitions of performance requirements;
- gross errors in design, during execution and use;
- lack of knowledge concerning behaviour of new materials and actions in actual
conditions.
The order of the listed uncertainties corresponds approximately to the decreasing level
of current knowledge and available theoretical tools for their description and consideration in
design (see following sections). It should be emphasized that most of the above listed
uncertainties (randomness, statistical and model uncertainties) can never be eliminated
absolutely and must be taken into account when designing any construction work.
2.2
I-2
RELIABILITY
3.1
General
The term "reliability" is often used very vaguely and deserves some clarification.
Often the concept of reliability is conceived in an absolute (black and white) way the
structure either is or isnt reliable. In accordance with this approach the positive statement is
understood in the sense that a failure of the structure will never occur. This interpretation is
unfortunately an oversimplification. Although it may be unpleasant and for many people
perhaps unacceptable, the hypothetical area of absolute reliability for most structures (apart
from exceptional cases) simply does not exist. Generally speaking, any structure may fail
(although with a small or negligible probability) even when it is declared as reliable.
The interpretation of the complementary (negative) statement is usually understood
more correctly: failures are accepted as a part of the real world and the probability or
frequency of their occurrence is then discussed. In fact in the design it is necessary to admit a
certain small probability that a failure may occur within the intended life of the structure.
Otherwise designing of civil structures would not be possible at all. What is then the correct
interpretation of the keyword reliability and what sense does the generally used statement
the structure is reliable or safe have?
3.2
Definition of reliability
A number of definitions of the term reliability are used in literature and in national
and international documents. ISO 2394 [2] provides a definition of reliability, which is similar
to the approach of national standards used in some European countries: reliability is the
ability of a structure to comply with given requirements under specified conditions during the
intended life, for which it was designed. In quantitative sense reliability may be defined as the
complement of the probability of failure.
Note that the above definition of reliability includes four important elements:
- given (performance) requirements definition of the structural failure,
- time period assessment of the required service-life T,
- reliability level assessment of the probability of failure Pf,
- conditions of use limiting input uncertainties.
An accurate determination of performance requirements and thus an accurate
specification of the term failure are of uttermost importance. In many cases, when considering
the requirements for stability and collapse of a structure, the specification of the failure is not
I-3
Probability of failure
The most important term used above (and in the theory of structural reliability) is
evidently the probability of failure Pf. In order to defined Pf properly it is assumed that
structural behaviour may be described by a set of basic variables X = [X1, X2, ... , Xn]
characterizing actions, mechanical properties, geometrical data and model uncertainties.
Furthermore it is assumed that the limit state (ultimate, serviceability, durability or fatigue) of
a structure is defined by the limit state function (or the performance function), usually written
in an implicit form as
Z(X) = 0
(1)
The limit state function Z(X) should be defined in such a way that for a favourable (safe) state
of a structure the function is positive, Z(X) 0, and for a unfavourable state (failure) of the
structure the limit state function is negative, Z(X) < 0 (a more detailed explanation is given in
the following Chapters of this Handbook 2).
For most limit states (ultimate, serviceability, durability and fatigue) the probability of
failure can be expressed as
Pf = P{Z(X) < 0}
(2)
The failure probability Pf can be assessed if basic variables X = [X1, X2, ... , Xn] are
described by appropriate probabilistic (numerical or analytical) models. Assuming that the
basic variables X = [X1, X2, ... , Xn] are described by time independent joint probability density
function X(x) then the probability Pf can be determined using the integral
Pf =
X ( x )dx
(3)
Z( X ) < 0
More complicated procedures need to be used when some of the basic variables are
time-dependent. Some details concerning theoretical models for time-dependent quantities
(mainly actions) and their use for the structural reliability analysis are given in other Chapters
of this Handbook 2. However, in many cases the problem may be transformed to a timeindependent one, for example by considering in equation (2) or (3) a minimum of the function
Z(X) over the reference period T.
I-4
Reliability index
An equivalent term to the failure probability is the reliability index , formally defined
as a negative value of a standardized normal variable corresponding to the probability of
failure Pf. Thus, the following relationship may be considered as a definition
= U1 ( Pf )
(4)
Here u1 ( p f ) denotes the inverse standardised normal distribution function. At present the
reliability index defined by equation (4) is a commonly used measure of structural
reliability in several international documents [1], [2], [5].
It should be emphasized that the failure probability Pf and the reliability index
represent fully equivalent reliability measures with one to one mutual correspondence given
by equation (4) and numerically illustrated in Table 1.
Pf
Table 1. Relationship between the failure probability Pf and the reliability index .
101
102
103
104
105
106
107
1,3
2,3
3,1
3,7
4,2
4,7
5,2
In EN 1990 [1] and ISO 2394 [2] the basic recommendation concerning a required
reliability level is often formulated in terms of the reliability index related to a certain
design working life.
3.5
DESIGN TARGETS
4.1
I-5
Examples
Temporary structures
Replaceable structural parts, bearings, girders
Agricultural and similar structures
Building structures and common structures
Monumental building or civil structures, bridges
4.2
5.1
General
During their historical development the design methods have been closely linked to
the available empirical, experimental as well as theoretical knowledge of mechanics and the
theory of probability. The development of various empirical methods for structural design
gradually crystallized in the twentieth century in three generally used methods, which are, in
various modifications, still applied in standards for structural design until today: the
permissible stresses method, the global factor and partial factor methods. All these methods
are often discussed and sometimes reviewed or updated.
The following short review of historical development illustrates general formats of
above mentioned design methods, indicate relevant measures that are applied to take into
account various uncertainties of basic variables and to control resulting structural reliability.
In addition a short description of probabilistic methods of structural reliability and their role
in further development of design procedures is provided. Detailed description of probabilistic
methods of structural reliability is given in Chapter II, Chapter III and in Annex B of this
Handbook 2.
I-6
Permissible stresses
The first of the worldwide-accepted design methods for structural design is the method
of permissible stresses that is based on linear elasticity theory. The basic design condition of
this method can be written in the form
(5)
The coefficient k (greater than 1) is the only explicit measure supposed to take into account all
types of uncertainties (some implicit measures may be hidden). Moreover, only a local effect
(a stress) max is compared with the permissible stress per and, therefore, a local (elastic)
behaviour of a structure is used to guarantee its reliability. No proper way is provided for
treating geometric non-linearity, stress distribution and ductility of structural materials and
members. For that reasons the permissible stress method leads usually to conservative and
uneconomical design.
However, the main insufficiency of the permissible stress method is lack of possibility
to consider uncertainties of individual basic variables and computational models used to
assess load effects and structural resistances. Consequently, reliability level of structures
exposed to different actions and made of different material may be not only conservative
(uneconomical) but also considerably different.
5.3
(6)
Thus the calculated safety factor s must be greater than its specified value s0 (for example s0 =
1,9 is commonly required for bending resistance of reinforced concrete members). The global
safety factor method attempts to take into account realistic assumptions concerning structural
behaviour of members and their cross-sections, geometric non-linearity, stress distribution
and ductility; in particular through the resulting quantities of structural resistance R and action
effect E.
However, as in the case of the permissible stresses method the main insufficiency of
this method remains a lack of possibility to consider the uncertainties of particular basic
quantities and theoretical models. The probability of failure can, again, be controlled by one
explicit quantity only, by the global safety factor s. Obviously harmonisation of reliability
degree of different structural members made of different materials is limited.
5.4
(7)
where the design values of action effect Ed and structural resistance Rd are assessed
considering the design values of basic variables describing the actions Fd = F Fk, material
properties fd = fk /m, dimensions ad + a and model uncertainties d. The design values of
these quantities are determined (taking into account various uncertainties) using their
characteristic values (Fk, fk, ak, k), partial factors , reduction factors and other measures of
reliability [1, 2, 3, 4], Thus the whole system of partial factors and other reliability elements
I-7
Probabilistic methods
The probabilistic design methods introduced in the International Standard [2] are
based on a requirement that during the service life of a structure T the probability of failure Pf
does not exceed the design value pd or the reliability index is greater than its design value d
Pf Pd or > d
(8)
In EN 1990 [1] the basic recommended reliability index for ultimate limit states d =
3,8 corresponds to the design failure probability Pd = 7,2 10-5, for serviceability limit states
d = 1,5 corresponds to Pd = 6,7 10-2. These values are related to the design working life of
50 years that is considered for building structures and common structures. In general greater
- values should be used when a short reference period (one or five years) will be used for
verification of structural reliability.
It should be mentioned that probabilistic methods are not yet commonly used in
design praxis. However, the developed calculation procedures and software products (for
example [7, 8] and [10]) already enable the direct verification of structural reliability using
probabilistic concepts and available experimental data. Recently developed software product
CodeCal [10] is primarily intended for calibration of codes based on the partial factor method.
In Chapter II of this Handbook 2 numerical examples will be presented to illustrate the
methods discussed above.
6
In some cases there is a need to base the design on a combination of tests and
calculations, for instance if no adequate calculation model is available. The tests may vary
from wind tunnel tests to prototype testing of new structural materials, elements or
assemblies. Tests may also be carried out during or after execution to confirm the design
assumptions. The extreme example is a proof load. For design by testing the following types
of tests can be distinguished:
a) tests to establish directly the resistance for given loading conditions
b) tests to obtain specific material properties
c) tests to reduce model uncertainties in loads, load effects or resistance models
Test should be set up and evaluated in such a way that the usual required level of
reliability is achieved. The derivation of a characteristic or design value should take into
account the scatter of test data, statistical uncertainty associated with the number of tests and
prior statistical knowledge. If the response of the structure or structural member or the
I-8
CONCLUDING REMARKS
The basic concepts of the probabilistic theory of reliability are characterized by two
equivalent terms, the probability of failure Pf and the reliability index . Although they
provide limited information on the actual frequency of failures, they remain the most
important and commonly used measures of structural reliability. Using these measures the
theory of structural reliability may be effectively applied for further harmonisation of
reliability elements and for extensions of the general methodology for new, innovative
structures and materials.
Historical review of the design methods worldwide accepted for verification of
structural members indicates different approaches to considering uncertainties of basic
variables and computational models. The permissible stresses method proves to be rather
conservative (and uneconomical). The global safety factor and partial factor methods lead to
similar results. Obviously, the partial factor method, accepted in the recent EN documents,
represents the most advanced design format leading to a suitable reliability level that is
relatively close to the level recommended in EN 1990 ( = 3,8). The most important
advantage of the partial factor method is the possibility to take into account uncertainty of
individual basic variables by adjusting (calibrating) the relevant partial factors and other
reliability elements.
Various reliability measures (characteristic values, partial and reduction factors) in the
new structural design codes using the partial factor format are partly based on probabilistic
methods of structural reliability, partly (to a great extent) on past empirical experiences.
Obviously the past experience depends on local conditions concerning climatic actions and
traditionally used construction materials. These aspects may be considerably different in
different countries. That is why a number of reliability elements and parameters in the present
suite of European standards are open for national choice.
I-9
I - 10
1 := 0 , 1 .. 6
2 Probability
pn ( 1 , n ) := 1 ( 1 p1( 1) )
p1( 1) := pnorm( 1 , 0 , 1)
3 Reliability index
n ( 1 , n ) := qnorm1 ( 1 pnorm( 1 , 0 , 1) ) , 0 , 1
n
4 Numerical results
5 Graphical results
1 =
5 .10
pn ( 1 , 50)
p1( 1) =
pn ( 1 , 50) =
0.5
0.159
0.023
0.684
1.3510-3
0.065
3.16710-5
1.58210-3
2.86710-7
1.43310-5
9.86610 -10
4.93310-8
1 := 0 , 0.1.. 6
1 .10
nt := 3.8
pnt := p1( nt )
pnt = 7.235 10
pn ( 1 , 1)
pn ( 1 , 5)
pnt
6
n( 1 , 1)
n( 1 , 5)
n( 1 , 50)
nt
4
1
I - 11
Summary
Elementary methods of structural reliability are described considering a fundamental
case of two random variables when the limit state function is formulated as a difference
between the resulting structural resistance and load effect. The initial assumption of normal
distribution of both resulting variables is generalised to any type of probability distribution.
The described computational procedures are illustrated by a number of numerical examples,
which are supplemented by MATHCAD and EXCEL sheets. An extension of the elementary
methods of structural reliability is presented in Annex B.
1
INTRODUCTION
1.1
Background materials
Fundamental concepts and procedures of structural reliability are well described in a
number of national standards, in the new European document EN 1990 [1] and International
Standard ISO 2394 [2]. Additional information may be found in the background document
developed by JCSS [3] and in recently published handbook to EN 1990 [4]. Guidance on
application of the probabilistic methods of structural reliability may be found in publications
and working materials developed by JCSS [5] and in relevant literature listed in [4] and [5].
1.2
General principles
The theory of structural reliability considers all basic variables as random quantities
having appropriate types of probability distribution. Different types of distributions should be
considered for actions, material properties and geometrical data. In addition, model
uncertainties of actions and resistance models should be taken into account. Prior theoretical
models of basic variables and procedures for probabilistic analysis are indicated in JCSS
documents [5].
2
2.1
General
The fundamental task of the theory of structural reliability concerns a basic
requirement for the relation between the action effect E and the structural resistance R written
in the form of inequality
E<R
II-1
(1)
(2)
Equation (2) represents a fundamental form of the failure boundary called the limit
state (performance) function (see also Chapter I of this Handbook 2). It should be noted that
for some structural members and materials the assumption of sharp failure boundary might be
rather artificial and can be accepted as an approximation only. Such a case is indicated in the
following Example 1.
Example 1.
A steel rod indicated in Figure 1 has a tensile resistance R
= d2 fy / 4, where d denotes the diameter of the rod and fy the
yield point. The rod is loaded by a weight E = V, where
V denotes the volume and the bulk weight density of the load.
Thus the inequality (1) has the form
V < d2 fy / 4
The limit state function (2) can be then written as
d2 fy / 4 - V = 0
In this example, the limit state is defined as the state when
the stress in the rod reaches the yield point fy. This simplification
is accepted in many common cases, but (depending on a type of
E
structural steel) it may not correspond to the actual failure of the
rod. In particular when structural steel with significant ductility
and strain hardening is used, then a failure (rupture) will occur Figure 1. A rod.
when the stress reaches the ultimate strength of the steel, which is
a considerably greater than the yield point.
Attached MATHCAD sheets SteelRod.mcd, DesVRod.mcd may be used to make all
numerical calculations.
Both the variables E and R are generally random variables and the validity of
inequality (1) cannot be guaranteed absolutely, i.e. with the probability equal to 1 (the total
certainty). Therefore, it is necessary to accept the fact that the limit state described by
equation (2) may be exceeded and failure may occur with a certain small probability. The
essential objective of the reliability theory is to assess the probability of failure Pf and to find
the necessary conditions for its limited magnitude. For the simple condition in the form of
inequality (1), the probability of failure may be formally written as
Pf = P(E > R)
(3)
The random character of the action effect E and the resistance R, both expressed in
terms of a suitable variable (performance indicator) X (i.e. stress, force, bending moment,
deflection) is usually described by appropriate distribution function, i.e. by distribution
functions E(x), R(x) and by corresponding probability density functions E(x), R(x), where
x denotes a general point of the considered variable X used to express both the variables E and
R. Distributions of variables E and R further depend on appropriate parameters, e.g. on
II-2
Load effect E,
Gamma distribution,
E = 70, E = 7
Resistance R
log-normal distribution,
R = 100, R = 10
0.04
0.02
0.00
40
60
80
100
120
140
Random variable X
Figure 2. Action effect E and resistance R as random variables.
The desired conditions will certainly include the trivial inequality E < R (see Figure
2). Obviously, this requirement for mutual position of both distributions is not sufficient to
ensure specified failure probability Pf. The correct conditions should certainly include also
conditions for variances of both variables. This will be clarified by the following discussion of
fundamental cases of structural reliability.
2.2
First, consider a special case when one of the variables E and R, say the action effect
E, has a very low (negligible) variability comparing to the variability of resistance R. Then E
may be considered as non-random (deterministic) variable, i.e. such a variable that attains a
certain fixed value E0 (E = E0) in its every realization. This assumption may certainly be
considered as an approximation of some practical cases. One of these cases is the loaded steel
II-3
(4)
The value of distribution function R(E0) is usually assessed from tables for a
standardized random variable U, for which the value u0 corresponding to E0 is computed. It
follows from the general transformation formula
u0 = (e0 - R) / R
(5)
Random resistance R:
Log-nomal distribution,
R = 100, R = 10
0,04
0,02
0,00
-u0R
40
60
80
100
120
140
Random variable X
Figure 3. Deterministic effect of actions E and random resistance R.
The probability of failure is then given as
Pf = P(R < e0) = R(e0) = U(u0)
(6)
where U(u0) is the value of distribution function of a standardized random variable of the
appropriate distribution (e.g. normal or log-normal).
Note that the value -u0 is the distance of the fixed value E0 of action effect E from the
mean R of resistance R expressed in the units of standard deviation R. If the distribution of
resistance R is normal, then the defined distance is called the reliability index
= (R e0) / R
II-4
(7)
(8)
= U1 ( p f )
(9)
where u1 ( p f ) denotes the inverse standardised normal distribution function. At present the
reliability index defined by equation (9) is a commonly used measure of structural
reliability in several national and international documents (see also previous Chapter I of this
Handbook 2). Note, however, that the probability distribution of the resistance R may differ
from the normal distribution.
Example 2.
Consider that resistance R has the mean R = 100 (expressed in dimensionless units),
standard deviation R =10 (the coefficient of variation is VR = 0,10). For the deterministic
action effect it holds that e0 = 80 (see Figure 3). If R has normal distribution, then the
reliability index follows directly from equation (7)
= (100 - 80) / 10 = 2
and probability of failure follows from relation (8)
Pf = P(R < 80) = u(-2) = 0,023
where u(-2) is the value of the distribution function of the standardized normal distribution
for u = 2. However, if the distribution of R is not normal but lognormal with the lower limit
at zero (skewness R = 3VR + V R 3= 0,301 [9]), then it follows from equation (5)
u0 = (80 100) / 10 = 2
The probability of failure Pf is then given as
Pf = P(R < 80) = LN,U(2) = 0,014
where LN,U(-2) is the distribution function of the standardized random variable U with lognormal distribution having the lower bound at zero (the skewness = 0,301). The resulting
probabilities do not much differ but their values are rather high.
If the fixed value of the action effect decreases to e0 = 70, then for normal distribution
of resistance R the reliability index is = 3 and probability of failure is
Pf = P(R < 70) = u(-3) = 0,00135
If the distribution of resistance R is log-normal with the lower limit at zero, then
Pf = P(R < 70) = LN,U(-3) = 0,00021
The reliability index defined by equation (9) is then = U1 (0,00021) = 3,53 , i.e.
greater than the value 3, which holds if normal distribution of resistance R is assumed.
Obviously, when the load effect is only e0 = 70 the resulting failure probabilities are
remarkably lower than in the case when e0 = 80. Furthermore, the numerical example also
shows that the assumption concerning the type of distribution plays an important role and may
be, in some cases, decisive.
II-5
(10)
called the safety margin, has the normal distribution with parameters
Z = R - E
(11)
2
2
2 2
Z2 = R + E + 2 RE R E
(12)
where RE is the coefficient of correlation of R and E. It is often assumed that R and E are
mutually independent and RE = 0. Equation (3) for the probability of failure Pf can now be
modified to
Pf = P(E > R) = P(Z < 0) = Z(0)
(13)
and the whole problem is reduced to determining the distribution function Z(z) for z = 0,
which leads to the probabilities of the safety margin Z being negative. The distribution
function Z(0) is usually determined by transformation of the variable Z to standardised
random variable U. Using this equation, the value u0 corresponding to the value g = 0 is
u0 = (0 - Z) /Z = - Z /Z
(14)
(15)
The probability density function Z(z) of the safety margin Z is shown in Figure 4,
where the grey area under the curve Z(z) corresponds to the failure probability Pf.
Probability density Z(g)
0,04
0,03
1 Pf
0,02
u0 Z
0,01
Pf
0,00
10
10
20
30
Safety margin Z
Figure 4. Distribution of the safety margin Z.
II-6
40
50
= Z / Z =
R E
R2 + E2 + 2 RE R2 E2
(16)
If the quantities R and E are mutually independent, then the coefficient of correlation
RE vanishes (RE=0). Thus, the reliability index is the distance of the mean Z of the safety
margin Z to the origin, given in the units of the standard deviation Z.
Example 3.
Consider again the Example 2, in which the resistance R and the load effect E are
mutually independent random variables (RE=0) having normal distribution. The resistance R
has the mean R = 100, variance R = 10 (coefficient of variation is therefore only w = 0,10),
and the effect of actions E has the mean E = 80 and E = 8 (all expressed in dimensionless
units). It follows from equation (11) and (12) that
Z = 100 80 = 20
2
2Z = 102 + 82 = 12,81
As both the basic variables R and E have normal distributions, the reliability index
follows directly from equation (16)
= 20 / 12,81 = 1,56
and the probability of failure follows from relation (8)
Pf = P(Z < 0) = U(-1,56) = 0,059
If the variables E and R are not normal, then the distribution of the safety margin G is
not normal either and then the above-described procedure has to be modified. In a general
case, numerical integration or transformation of both variables into variables with normal
distribution can be used. The transformation into a normal distribution is primarily used in
software products.
An approximate simple procedure can be used for a first assessment of the failure
probability Pf. The safety margin Z may be approximated by a three-parameter lognormal
distribution. Assume that the distributions of E and R depend on the moment parameters E,
E, E, R, R and R. The mean and variance of the safety margin Z may be assessed from
the previous equations (11) and (12), which hold for variables with an arbitrary distribution.
Assuming mutual independence of E and R, the skewness Z of the safety Z may be estimated
using the approximate formula (see Annex A - Basic statistical concepts and techniques in
this Handbook 2)
Z =
R3 R E3 E
( R2 + R2 ) 3 / 2
(17)
Then it is assumed that the safety margin Z can be described with sufficient accuracy
by a log-normal distribution with determined moment parameters Z, Z and Z (equations
(11), (12) and (17)). It shows that this approximation offers satisfactory results if the
probability of failure is not too small.
II-7
Z = R - E = 100 50 = 50
G2 = 2R + 2E = 102 + 102 = 14,14
Z =
3
3
R R E E
2
R
+ 2E )
3/ 2
3
3
10 0,301 10 1,14
= 0,30
(102 + 102 )3 / 2
In the case of two random variables E and R having any distribution, the exact
determination of the failure probability Pf, defined by equation (3), may be obtained by
probability integration. Figure 5 is used to explain this procedure. Let the event A denote the
occurrence of the action effect E in the differential interval <x, x+dx>. Probability of the event
A is given as
P(A) = P(x < E < x+dx) = E(x) dx
(18)
Let us denote B as the event that resistance R occurs within the interval <-, x>.
Probability of the event B is [9] given as
P(B) = P(R < x) = R(x)
(19)
(20)
0,06
Resistance R
0,04
dx
0,02
0,00
40
60
x x+dx
100
120
140
Random variable X
Figure 5. Distribution of variables E and R.
The integration of the differential relationship (20) over the interval in which both
variables E and R occur simultaneously (generally the interval <-, >) leads to the relation
p f = R ( x ) E ( x ) dx
(21)
the integration of the relation (21) usually has to be carried out numerically or using the
simulation methods (e.g. direct Monte Carlo methods).
Attached MATHCAD sheet PrLnLn.mcd offers a simple programme that can be used
to evaluate the numerical integration of relation (21) assuming that both the variables E and R
can be described (at least approximately) by the general (three-parameter) lognormal
distribution.
Example 5.
The action effect E and the resistance R are described by a log-normal distribution
with the same parameters as in Example 4 (the Gumbel distribution for E was simply
substituted by the log-normal distribution having the same parameters). The approximate
solution in Example 4, based on log-normal distribution with the lower bound at zero, leads to
the probability of failure Pf = P(R < E) = LN,U(-3,54) = 0,00101. The numerical integration
according to relation (21) using the programme MATHCAD leads to a solution Pf = P(R < E)
= 0,000792, the programme VaP suggests a solution Pf = P(R < E) = 0,000707, which can be
considered as a very good approximation.
The probability of failure Pf assessed by the direct integration may be determined
using MATHCAD sheet PrLnLn.mcd for the given parameters of variables E and R (R =
II-9
0,5
-2
1,5
-1
Skew ness R
1
1,00E -04
1,00E -05
pf
1,00E -06
CONCLUDING REMARKS
Elementary methods of the structural reliability can be used to assess the reliability of
fundamental cases of two random variables when the limit state function is formulated as the
difference between the resulting structural resistance and load effect.
Basic principles of the reliability theory provide operational techniques that can be
used for estimating the partial factors of basic variables. The assessment of various reliability
measures in the new structural design codes is, however, partly based on historical and
empirical experiences. Obviously the past experience depends on local conditions including
climatic actions and traditionally used construction materials and, consequently, might be in
different countries considerably diverse. That is why number of reliability elements and
parameters in the present suite of European standards are open for national choice.
II-10
ATTACHMENTS
1. MATHCAD sheet SteelRod.mcd
Mathcad sheet SteelRod is intended to investigate an effect of the partial factor G on
reliability of a steel rod exposed to permanent load G.
2. MATHCAD sheet DesVRod.mcd
Mathcad sheet DesVRod is intended to investigate of sensitivity factor E and R and
design values Ed and Rd.
3. MATHCAD sheet PrLnLn.mcd
Mathcad sheet PrLnLn is intended for calculation of the failure probability Pf =
P{E>R}based on approximation of E and R by three parameter lognormal distribution.
II-11
Gk := 1
A ( G) :=
fk := 235
( Gk G)
Check:
fd
m := 1.10
fd :=
fk
m
3
A ( 1.35) = 6.32 10
280
vG := 0.1 G := vG G :=
XS := 1 XS := 0
XR := 1
f := fk vf := 0.08 f := vf f
235
XR := 0.00
vXR:=
XR
vXS:=
XR
XS
XS
R ( G) := f XR A ( G) E := G XS
vXR + vXR vf + vf
vE :=
R ( 1.35) = 1.77
2 Check:
vR = 0.08
vXS + vXS vG + vG
R := 3 vR + vR
E = 1
vE = 0.1
E := 2 vE
R ( G) := vR R ( G)
( R ( G) ) + ( E )
3
g ( G) :=
R ( 1.35) = 0.14
g ( 1.35) = 0.77
R R ( G) E E
g ( G)
E := vE E
g ( 1.35) = 0.17
g ( 1.35) = 0.09
g ( G)
g ( G)
Check:
Pf0 ( G) := pnorm( 0 ( G) , 0 , 1)
0 ( 1.35) = 4.44
C ( G) :=
g ( G) + 4 + g ( G) g ( G) + 4 g ( G)
1
ln 1 + C ( G)
mg ( G) := ln ( C ( G)
x0 ( G) := g ( G)
1
C ( G)
) + ln ( g ( G) ) ( 0.5) ln (1 + C ( G) 2)
g ( G)
Check:
II-12
x0 ( 1.35) = 4.85
1 ( 1.35) = 4.76
En ( x ) := dnorm( x , E , E)
E
k :=
E Eg ( x ) := dgamma( x , k)
2
E
:=
C ( G) :=
R ( G)
( R ( G) a ( G) )
m ( G) := ln ( R ( G) ) ln ( C ( G) ) ( 0.5) ln 1 + C ( G)
Probability lognormal distribution of R
aR ( G) := C ( G) + 3 C ( G)
s ( G) :=
ln 1 + C ( G)
Rln ( x , G) := plnorm [ ( x a ( G) ) , m ( G) , s ( G) ]
t := 3.8
Pfn ( G) :=
En ( x ) Rln ( x , G) dx
n ( G) := qnorm( Pfn ( G) , 0 , 1)
Pfg ( G) :=
Eg ( x ) Rln ( x , G) dx
g ( G) := qnorm( Pfg ( G) , 0 , 1)
7 Parametric study of G
Check:
0 ( G)
7
6
7
Reliability index
1 ( G)
g ( G)
n ( G)
t
1 ( 1.35) = 4.76
g ( 1.35) = 4.61
n ( 1.35) = 4.85
1.2
1.4
1.6
Note: Reliability assessment assuming normal distribution for E and R seems to be on a safe side
(leads to a lower bound for ], while assessment assuming three parameter distribution for the
reliability margin g seems to provide a more realistic estimate.
II-13
Gk := 1
A ( G) :=
fk := 235
( Gk G)
fk
m
3
Check:
fd
fd :=
m := 1.10
A ( 1.35) = 6.32 10
vG := 0.1 G := vG G :=
XS := 1 XS := 0
XR := 1
280
f := fk vf := 0.08 f := vf f
235
XR := 0.00
vXR:=
XR
vXS:=
XR
XS
XS
R ( G) := f XR A ( G) E := G XS
CoV: vR :=
vXR + vXR vf + vf
vE :=
R ( 1.35) = 1.77
2 Check:
vR = 0.08
vXS + vXS vG + vG
R := 3 vR + vR
E = 1
vE = 0.1
E := 2 vE
R ( G) := vR R ( G)
g ( G) :=
( R ( G) ) + ( E)
g ( G) :=
R ( 1.35) = 0.14
g ( 1.35) = 0.77
R R ( G) E E
g ( G)
E := vE E
g ( 1.35) = 0.17
g ( 1.35) = 0.09
E
g ( G)
R ( G) :=
R ( G)
g ( G)
E( G)
R( G)
0.8
0.6
0.4
1.2
1.4
1.6
Note that the sensitivity factor E is shown with the oposite sign (as a positive quantity).
II-14
:= 3.8
E0 := 0.7
R0 := 0.8
Ed ( G) := E E ( G) E
Rd ( G) := R ( G) R ( G) R ( G)
Ed0 ( G) := E E0 E
Rd0( G) := R ( G) R0 R ( G)
Ed ( 1.35) = 1.22
Ed0 ( 1.35) = 1.27
1.8
Ed ( G)
Rd ( G)
Rd ( 1.35) = 1.33
Rd0( 1.35) = 1.34
Rd0ln( 1.35) = 1.39
1.6
1.4
Ed0 ( G)
Rd0 ( G)
Rd0ln ( G)
1.2
1
0.8
1.2
1.4
1.6
Notes: 1) Figure shows that the partial factor G should be greater than about 1,25 otherwise the
design value of the load effect E d would be greater than the design value of the resistance Rd.
2) The design value of the resistance Rd determined assuming lognormal distribution with
the lower bound at zero is greater than Rd determined assuming the normal distribution.
II-15
E
E
E2 + 4 + E
CE( E) :=
3
x0E( E) := E
R2 + 4 + R
CR( R) :=
3
x0R( R) := R
x0E( 0) =
Check:
R + 4 R
CR( 0.3) =
if R 0
CR( R)
CE( 0.301) =
E 6 E otherwise
Distribution parameter CR
given by the skewness R:
Check:
E + 4 E
if E 0
CE( E)
R 6 R otherwise
x0R( 0.000) =
x1( E, R) :=
otherwise
R + 6 R if R 0
x0( 0.608, 1) =
x0R( R) otherwise
x1( 1 , 1) =
uuE( x, E) :=
uE( x) :=
ln uE( x) +
CE( E)
(x E)
+ ln CE( E) 1 + CE( E)
if E 0
sign ( E) ln 1 + CE( E)
uE( x) otherwise
E( x, E) :=
dnorm( uuE( x, E) , 0 , 1)
E uE( x) +
CE( E)
dnorm( uuE( x, E) , 0 , 1)
E
Standardised variable R:
ln 1 + CE( E)
uR ( x) :=
if CE( E) 0
E( 50, 0) =
E( 50, 0.0001) =
otherwise
( x R)
R
II-16
ln uR ( x) +
uuR ( x, R) :=
+ ln CR( R) 1 + CR( R) 2
CR( R)
if R 0
1
sign ( R) ln 1 + CR( R)
uR ( x) otherwise
Distribution function LN,X(x) = LN,U(u) = (uu):
R( x, R) := pnorm( uuR ( x, R) , 0 , 1)
12
uuR ( 50, 0) = 5
R( 0 , 0.3) = 1.237 10
pf ( E, R) :=
x1( E , R)
E( x, E) R( x, R) dx
x0( E , R)
0.01
1 .10
1 .10
1 .10
0.5
1.5
sE( E) := ln 1 + CE( E)
Distribution function of R:
pf ( E, R) :=
sR ( R) := ln 1 + CR( R)
Check:
4
E( x, E) R( x, R) dx
x0( E , R)
II-17
Summary
Basic reliability elements specified in current standards for structural design
commonly include failure probability related to a certain reference period T. Required
reliability level of buildings and other civil engineering works is usually specified by the
design (target) failure probability pd or by appropriate reliability index d corresponding to a
specified design working life Td (for example 50 years). In reliability verification the design
values d and Td are sometimes replaced with an alternative reliability index a derived from
the design values d and Td for a convenient reference period Ta (for example 1 year).
Submitted study clarifies relationships between the alternative elements a, Ta and
design values d, Td, and indicates relevant procedures for reliability verification when
alternative reference period Ta is considered. It is emphasised that verification based on a, Ta
should be distinguished from verification of temporary or auxiliary structures when the design
working life Td itself is short. Theoretical consideration and numerical examples show that
characteristic values and partial factors of basic variables describing material properties and
self-weight are significantly dependent on the relevant reference period.
1
INTRODUCTION
1.1
Background documents
Recent documents [1], national [2], [3] and international documents ([4] to [7])
provide general principles and guidance for application of probabilistic methods to structural
designs. The latest European document [5] and international standards [6] and [7] also
indicate a theoretical basis of the so called partial factor method and procedures for
determination of partial factors of material properties and actions using probabilistic
principles.
The basic reliability elements considered in these procedures include probability of
failure p (or equivalent reliability index ) corresponding to a certain reference period T used
in verification of structural reliability. The reference period T used in verification may or may
not coincide with the design working life Td, which is the time period during which a structure
is required to perform adequately. When the reference period used in reliability verification is
different from Td then it is called an alternative period and denoted in this study Ta.
1.2
General Principles
Basic probabilistic methods are used to analyse principles of reliability differentiation.
Similarly as in Chapter I in this Handbook two essentially different cases are distinguished in
the following:
- an alternative reference period Ta (for example 1 or 5 years), which is different from
the design working life Td (for example 50 years), is considered; this case is applicable when
probabilistic models related to the period Ta are more credible than those related to Td;
III - 1
The basic reliability measures include the probability of failure and reliability index as
introduced in Chapter I and II in this Handbook. The probability of structural failure Pf can be
generally defined as
Pf = P{Z(X) < 0}
(1)
The limit state (performance) function Z(X) is formulated in such a way that the
reliable (safe) domain of a vector of basic variables X = X1, X2, ... , Xn corresponds to the
inequality Z(X) > 0 while the failure domain to the complementary inequality Z(X) < 0. A
simple example of Z(X) describes the basic relationship between the resulting load effect E
and resistance R
Z(X) = Z = R E
(2)
The random variable Z in equation (2) is often called the reliability (safety) margin; its
mean Z, standard deviation Z and skewness Z may be derived from corresponding
characteristics of resulting variables R and E as indicated in Chapter II.
Instead of the failure probability Pf, the reliability index is frequently used in
reliability consideration as an equivalent quantity to Pf. The reliability index is related to the
failure probability Pf as already indicated in Chapter I
Pf = ()
(3)
(4)
(5)
where pd denotes specified design (target) failure probability corresponding to the target
reliability index d.
Conditions (4) or (5) have to be used by designers when probabilistic methods are
applied for verification of structural reliability. Indicative target values pd and d are declared
in some national standards (e.g. [2] and [3]) and recently also specified in international
documents (e.g. [4] to [7]) for various design conditions (limit states, failure consequences
and economic aspects).
3
Design working life Td is an assumed period of time for which a structure or part of it
is to be used for its intended purpose with anticipated maintenance but without major repair
being necessary. In recent documents of CEN [5] and ISO [6] indicative values of Td are
provided for five categories of structures as shown in Chapter I of this Handbook.
More detailed specification of structural categories and design working lives is
available in the ISO documents [6, 7]. In general the design working lives indicated in [2] are
greater (in some cases by 100 %) than those given in Chapter I. For example the design
working life for temporary structures indicated in [2] is 15 years, for agricultural structures 50
years, for apartment and office buildings 100 years, and for railway structures, dams, tunnels
and other underground engineering works 120 years.
Design failure probabilities pd are usually indicated in relation to the expected social
and economical consequences. Table 1 shows classification of target reliability levels
provided in EN 1990 [5]. Reliability indexes are given for two reference periods T (1 year
and 50 years) only, without any explicit link to the design working life Td. Similar -values as
in Table 1 are given in [3] for the ultimate limit states, for which, however, the design
working life Td = 80 years (for building structures) is considered.
It should be underlined that a couple of values (a and d) specified in Table 1 for
each reliability class (for 1 year and 50 years) correspond to the same reliability level.
Practical application of these values depends on the time period Ta considered in the
verification, which may be connected with available information concerning time variant
vector of basic variables X = X1, X2, ..., Xn. For example, if the reliability class 2 and 50 years
design working period is considered, then the reliability index d = 3,8 should be used in the
verification of structural reliability. The same reliability level corresponding to class 2 is
achieved when the time period Ta = 1 year and a = 4,7 is used. Thus, various reference
periods Ta, in general different from the design working life Td, may be used for achieving a
certain reliability level.
Table 1. Reliability classification in accordance with CEN [5]
Reliability Consequences for loss of
Reliability index
classes
human life, economical,
a for Ta= d for Td=
social and environmental
1 year
50 years
consequences
3 high
High
5,2
4,3
2 normal Medium
4,7
3,8
1 low
Low
4,2
III - 3
3,3
When the vector of basic variables X = X1, X2, ... , Xm is time variant, then failure
probability p is also time variant and should always be related to a certain reference period T,
which may be generally different from the design working life Td. Considering a structure of a
given reliability level, the design failure probability pd = pn related to a reference period Tn = n
T1 can be derived from the alternative probability pa = p1 corresponding to Ta = T1 (to simplify
notation note that previously used subscript "d" corresponds now to "n" and subscript "a" to
"1") using approximate relationship given in [6], [7]
Pn = 1 (1 P1)n
(6)
(7)
Here (.) denotes the distribution function of standardised normal distribution. Figure
1 shows variation of n with 1 for n = 5, 25, 50 and 100. Note that, if the reference period T1
is one year, then n indicates the number of years of the reference period Tn (n = Tn).
III - 4
5
4,5 3,8 n = 5
25
4
50
3,5
100
3
2,5
2
1,5
1
0,5
0
2,5
3
3,5
4,7
4
4,5
5,5
1
Figure 1. Variation of n with 1 for n = 5, 25, 50 and 100
Note that, if 1-year period would be used for specification of the target reliability level
of a structure, then Figure 1 provides information on the resulting failure probability
corresponding to a given working life Tn. For example, if the target reliability level is
specified by the reliability index 1 = 4,7 (corresponding to the probability p1 = 1,3 10-6),
then (as already mentioned) the reliability level of a structure having a working life, Tn = 50
years is characterised by 50 = 3,8. Similarly when a period Tn = 5 years is used, then 5 = 4,3
or when Tn = 100 years, then 100 = 3,6.
So, the reliability level of a structure can be specified using different time periods T,
which may not necessarily coincide with the design working life Td. This may be useful when
experimental data concerning time variant basic variables are available for a specific
reference period T (for example 1 or 5 years) that is different from the design working life Td.
In such a case, however, all the basic variables (including those that are time independent)
should be considered by appropriate design values related to the same reference period T.
The following simple example indicates the effect of using a reference period T different from
the design working life considering a resistance variable (strength) having lognormal
distribution.
5
III - 5
(8)
Rd = R exp (R a VR)
(9)
Taking into account equations (8) and (9) it follows that the partial factor is given as
(10)
Considering selected values of the coefficient of variation VR, Figure 2 shows the
partial factor R for lognormal distribution of R (equation (10)).
It follows from Figure 2 that when reliability of a structure is verified using a short
alternative reference period Ta (for example for example for Ta = 1 year when a = 4,7), the
partial factor R should generally be greater than in the case when the whole design working
life Td (for example for Td = 50 when d = 3,8) is considered. It may be noted that the partial
factor R of material property R increases with the increasing value of the reliability index a.
Similar conclusions can be expected for partial factors of other basic variables, in
particular for partial factors of permanent actions.
2,5
R
2
VR = 0,20
0,15
0,10
0,05
1,5
0,5
0
2 a
(11)
(12)
In equation (11) and (12) G denotes the mean, G the standard deviation, VG the coefficient
of variation and G = 0,7 the sensitivity factor of G. The partial factor G of G is defined as
[5], [6] a [7]
G = Gd / Gk
(13)
Taking into account equations (11) and (12) it follows from (13) that
G = (1 + 0,7 a VG)
(14)
Figure 3 shows variation of the partial factor G with the reliability index a for
selected values of the coefficient of variation VG = 0,05; 0,10; 0,15 and 0,20. Note that G =
1,35 (recommended in EN 1990 [5]) corresponds approximately to the reliability index a =
3,8 if the coefficient of variation is about 0,1 (the value in EN 1990 [5] was increased by 5%
to take into account model uncertainty).
2,5
G
2
VG = 0,20
0,15
0,10
0,05
1,5
0,5
0
2 a
Drafts of European documents for climatic actions due to temperature [8], snow [9]
and wind [10] indicate possible reduction of characteristic values Qk for temperature, snow
load and wind speed in case of shorter reference (return) period (for example 5 years) than 50
years considered in normal cases. Such a reduction may be applied in transient design
situations (for example during execution).
III - 7
(15)
(16)
1 Vs
1 Vs
(17)
where sk,50 is the characteristic snow load on the ground for 50-year return period and sk,n for
n-year return period, p denotes here the probability of sk,n being exceeded corresponding to n
years of return period and Vs is the coefficient of variation of annual maximum snow load,
(c) In accordance with EN 1991-1-4 [10] the basic wind speed vb,n having the return
period n years may be assessed using semi-empirical expressions
1 K ln( ln(1 p ))
vb,n = k vb,50, where k =
1 K ln( ln(0,98))
0 ,5
(18)
where vb,50 is the basic wind velocity for 50-year return period and vb,n for n-year return period
and p denotes here the probability of vb,n being exceeded corresponding to n years of return
period. The constant K in equation (18) follows from Gumbel distribution as K = Vv6/,
where Vv denotes coefficient of variation of annual wind speed. An approximate value K = 0,2
(which corresponds to the coefficient of variation Vv= 0,26) is used in the following
comparison of reduction coefficients k for considered climatic actions.
Table 3 shows reduction coefficients k for climatic actions (applied in a general
relationship Qk,n = k Qk,50) for selected return periods of n - years.
Table 3. Reduction coefficient k for climatic actions (Qk,n = k Qk,50) for different return
periods of n - years.
Return period
p
Reduction coefficient k for
of n-years
Tmax,n
Tmin,n
sn,n
vb,n
2 years
0,5
0,8
0,45
0,64
0,77
5 years
0,2
0,86
0,63
0,75
0,85
10 years
0,1
0,91
0,74
0,83
0,90
50 years
0,02
III - 8
EXAMPLES
Consider a steel structure having the design working life Td = 50 years, for which the
target failure probability is specified as pd = 7,2 10 5 (d = 3,8). Failure probability p for the
alternative reference period Ta = 1 year, which is considered in design due to data concerning
actions, will be lower than the target failure probability pd (p < pd and > d); from equation
(6):
pa = 1 - (1-7,2 10 5 )1/ 50 = 1,44 10 6
When the reference period Ta = T1 = 1 year is considered in design verification, then
the reliability index follows from equation (7) as
1 = -1(1,44 10 6) = 4,7
Reliability index 1 is greater than the target value d = 3,8 specified for the design working
life Td = 50 years.
Using equation (10) the partial safety factor R for Ta = T1 = 1 year assuming the
coefficient of variation VR = 0,08 (corresponding to the common variability of strength of
structural steel) the partial safety factor is given as (see also Figure 2)
III - 9
CONCLUDING REMARKS
(1) In present international documents the target values of failure are related to economic
aspects of safety measures and consequences of structural failure only vaguely, without
any explicit relation to various design working lives Td for different types of structures.
(2) When alternative failure probability pa is derived for a suitable reference period Ta from
the target failure probability pd and design working life Td, partial factors and
characteristic values of variable actions for pa and Ta should also be specified.
(3) For temporary structures, with a short design working life Td, the target failure probability
pd can be specified in accordance with the general principles of reliability differentiation;
reliability elements for basic variables should be derived for specified pd and Td.
(4) The partial factors derived for an alternative reference period Ta different from Td may
vary considerably from the values corresponding to the design working life Td depending
on Ta and distributions of relevant basic variables.
(5) The partial factor of self-weight G corresponding to an alternative reference period Ta
varies with a-values less significantly than the partial factor of material property R.
(6) Partial factors R derived for an alternative reference period Ta of one-year may be
considerably greater than R specified for the design working life Td.
(7) Following recommendations of Eurocodes, the characteristic value for climatic actions
due to snow corresponding to 5-year return (reference) period may be reduced to 75 % of
the characteristic values for 50-year return period, similarly the characteristic value of
wind speed may be reduced to 85 %, the maximum temperature to 86%, the minimum
temperature to 63%.
REFERENCES
[1]
[2]
III - 10
ATTACHMENTS
1. MATHCAD sheet GammaRG.mcd
MATHCAD sheet Gamma is intended for determination of the partial factor R of the
resistance R and the partial factor G of the permanent load G.
2. MATHCAD sheet PSI0.mcd
MATHCAD sheet PSI0 is intended for determination of the Combination factor 0
for accompanying action.
3. MATHCAD sheet PSI12.mcd
MATHCAD sheet PSI12 is intended for determination of the combination factor 12
for accompanying action.
III - 11
k := 1.65
d := 3.09
Sensitivity factors:
Prameters:
R := 0.8
:= 0 , 0.1.. 5
R( ) := R
E := 0.7
E( ) := E
V := 0 , 0.1.. 0.5
kn( V) := ( 1 k V)
Normal distribution
kln( V) :=
dn ( , V) := ( 1 R( ) V)
exp( k) ln 1 + V
dln( V) :=
1+ V
dln ( , V) :=
1+ V
exp( R( ) ) ln 1 + V
1+ V
Rn( , V) :=
GammaR
kn( V)
dn ( , V)
Rln( , V) :=
1.5
Rn( , 0.1)
kln( V)
dln( , V)
1.5
Rn( , 0.05)
Rln( , 0.05)
Rln( , 0.1)
Rn( , 0.15)
Rln( , 0.15)
0.5
0.5
Gn( , 0.05)
1.4
Gn( , 0.15)
Gn( 3.3 , V)
Gn( , 0.1)
Gn( 3.8 , V)
1.2
Gn( 4.3 , V)
exp( d ) ln 1 + V
0.2
0.4
V
III - 12
Check:
1 + 0.7 V
Approximation in EC 1990
0na( V, r) :=
1 + 0.7 V
qgamma( pnorm( c( r) , 0 , 1) , V )
r
pnorm( 0.7 , 0 , 1) , 0 , 1
c( r) := qnorm
c( 7) = 3.259
0g( V, 7)
0n( V, 7)
0.5
0na( V , 7)
0d( V, 7)
0.2
0.4
0.6
0.8
V=
0n( V, 10) =
0.1
0.727
0.747
0.664
0.2
0.548
0.581
0.474
0.3
0.423
0.465
0.352
0.4
0.33
0.378
0.268
0.5
0.258
0.312
0.205
III - 13
T := 0.5
V := 0 , 0.1.. 1
0.8
0g( 3.8 , V , N , 0.5 )
0g( 3.8 , V , N , 0.6 )0.6
0g( 3.8 , V , N , 0.7 )
0.4
0g( 3.8 , V , N , 1 )
0.2
0.2
0.4
0.6
0.8
Variation of Q with T
.
T := 0 , 0.1 .. 1
V := 0.3
0.8
0g( 3.3 , V , N , T )0.6
0g( 3.8 , V , N , T )
0g( 4.3 , V , N , T )0.4
0.2
0.2
0.4
0.6
0.8
III - 14
w := 0 , 0.1.. 1.1
:= 0 , 0.02.. 1.01
Rreliability index
:= 3.8
1 + qnorm( , 0 , 1) w
1 + qnorm( 0.98, 0 , 1) w
0.8
0.6
12g( w , 0.5)
12( w , 0.5)
0.4
0.2
.
0
0.2
0.4
0.6
0.8
The Gumbel distribution leads to a lower 12 factor than the normal distribution
III - 15
0.8
0.6
0.4
0.2
0.2
0.4
0.6
0.8
Examples of probability = 1 - / q where is fraction of the refence period (0.01 or 0.5) during which
Q1 and Q2 are exceeded, p probabilioty of Q being non zero
- imposed w=1.1,
short term 1: on 18 days a year = 1 - / q = 1- 0.01/0.05 ~ 0.8
~ 0.5
2: almost always on = 1 - / q = 1- 0.5/1
long term
= 1 - / q = 1- 0.01/0.009 ~ 0.1
= 1 - / q = 1- 0.5/0.009 ~ NA >> 0
2 ~ 0.0
1~ 0.2
= 1 - / q = 1- 0.01/0.014 ~ 0.3
=
1
/
q
=
10.5/0.014
~
NA
>>
0
2 0.0
12g( w , , q ) :=
q := 0 , 0.01.. 1
1 0.78 w 0.58 + ln ln 1
III - 16
0.8
12g( 1.1 , 0.01 , q)
0.6
12g( 0.5 , 0.01 , q)
12g( 1.1 , 0.5 , q)
12g( 0.5 , 0.5 , q)
0.4
0.2
0.2
0.4
0.6
0.8
III - 17
Summary
Under particular circumstances it may be favorable or necessary to carry out tests in
order to obtain certain design parameters. Typical parameters determined from the tests are
actions on the structure, resistance of the structure or structural component and material
properties. Tests can be performed also to calibrate parameters in the theoretical model of
resistance. The design value of the parameter is obtained from the test results as the estimated
value of a certain fractile of the parameter in question. The procedures are explained for the
determination of a single property and for the determination of a probabilistic model of
resistance.
1
INTRODUCTION
1.1
Background materials
The section 5.2 of the European standard EN 1990 [1] mentions briefly the most
general principles of the design assisted by testing and refers to Annex D of the same
standard, where the procedures are dealt with in detail. The international standard ISO 2394
[2] also explains the procedures of the design assisted by testing in its Annex D. The two
standards differ in a number of details regarding this subject. Some insight in the derivation of
the statistical formulations is given in the international standard ISO 12491 [3] and in
statistical literature (e.g. [5]).
1.2
General principles
Under particular circumstances it may be favorable or necessary to carry out the tests
in order to obtain certain design parameters. The examples of possible such circumstances
include:
calculation models are lacking or are inadequate;
large number of similar components will be used;
cases when the calculation model leads to very conservative results;
derivation of new design formulae;
confirmation of assumptions made in the design.
The unknown quantities which are evaluated as a result of the tests may be
actions on the structure (e.g. wind loads);
structural response under loading or accidental effect;
strength or stiffness of the structure or structural element.
The level of reliability of a structure designed by testing should be at least the same as
for structures designed only by calculation models.
IV - 1
(1)
D((aiXi))=ai2D(Xi)
(2)
where E(X) and D(X) represent the mean value (expectation) and the variance (dispersion) of
the random variable X. The statement about the mean value holds even if the variables XI are
not independent. Note also, that X is distributed arbitrarily.
We can look at the sample x1,x2,...,xn as a realization of the random variables
X1,X2,...,Xn. Then M=(Xi)/n is also a random variable whose realization is the average value
of the sample. If we put ai=1/n in the above equations, we get the expressions for the mean m
and standard deviation s of the sample average value as follows:
m = E(M) = E((Xi)/n) = (E(Xi))/n = n /n =
(3)
(4)
IV - 2
(5)
and also decreases with the number of samples. The ratio VM/VX is shown on the Figure 1.
1
0,9
0,8
VM /Vx
0,7
0,6
0,5
0,4
0,3
0,2
0,1
0
0
10
20
30
40
50
60
70
80
90
100
number of samples n
IV - 3
1
2
3
4
5
6
7
8
9
10
11
12
7,98
8
8
7,99
7,98
7,95
7,98
7,99
7,96
7,99
7,97
7,98
816
845
832
830
846
821
826
822
841
807
831
830
924
944
948
925
969
937
928
934
956
946
942
926
13
14
15
16
17
18
19
20
21
22
23
24
7,98
7,98
7,97
8
8
7,98
8
7,98
7,98
7,98
7,98
8
832
811
840
839
855
830
833
826
836
843
840
847
949
932
937
934
943
928
934
934
942
948
937
928
25
26
27
28
29
30
31
32
33
34
35
36
d
8
7,98
7,99
7,97
7,98
7,99
7,98
7,98
7,97
7,99
7,98
8
818
828
817
851
855
847
822
836
830
830
845
829
907
940
941
959
970
947
921
925
945
938
945
934
37
38
39
40
41
42
43
44
45
7,98
7,98
7,96
7,99
7,98
8
7,98
7,98
8
829
810
832
823
826
829
815
826
823
930
904
925
916
957
931
910
942
932
Vm/Vx
1
analytical
yield stress
tensile strength
diameter
0,8
0,6
0,4
0,2
0
0 5 10 15 20 25 30 35 40 45 50
number of samples n
2.1
General principles
This section gives expressions for deriving design values of the ultimate resistance or
serviceability parameters of a structure or a component and for deriving the design values of
material properties.
It is assumed that all variables follow normal or lognormal distribution and that there
is no prior knowledge about the value of the mean. Two cases are considered regarding the
IV - 4
(6)
(7)
As we have shown (equations (3) and (5)) the sample mean mx is a random variable
with the mean X and the standard deviation X/n. If we multiply the above inequality by -1,
add a value mx and divide by X/n , we get
P ((mx - x)n /x kn = u) =
(8)
where the expression on the left side of the inequality is standardized normal variable and u
is a fractile of the standardized normal distribution corresponding to the probability . The
value k is therefore:
k = u /n
(9)
Example 4.
We have the same situation as in the Example 3, but now the standard deviation x is
unknown. We are searching k such that we will have the probability that the true mean x of
the parameter X will be greater then the sample mean mx, according to equation:
Xk = mx k sx x
(10)
We make use of the following theorem: if all the terms are the same as in Example 3
and sx is the sample standard deviation:
sx = (xi - mx)2/(n-1)
(11)
(mx - x)n / sx
(12)
(13)
with t being the fractile of the t-distribution of degree n-1 corresponding to the probability .
IV - 5
(14)
(15)
This equation is equal to equation (6) with the coefficient of variation Vx given by:
Vx= x/mx
(16)
The m is the partial factor for the parameter X and it should be taken from the appropriate
Eurocode EN 1992 to EN 1998. The d is the design value of the conversion factor. This
factor covers the differences between the laboratory test conditions and conditions during the
actual use. The value kn is obtained from the prediction method of fractile estimation and is:
kn = -up (1/n + 1)1/2
(17)
10
1,72
20
1,68
30
1,67
1,64
1,92
1,76
1,73
1,64
The numbers in the Table are actually based on assumptions x known and x
unknown. In EN 1990 [1] these assumptions are replaced with assumption of knowledge of
the coefficient of variation VX. If the standard deviation x is known then Vx should be
computed from the equation (16). If x or Vx is unknown, then Vx is calculated from the
sample standard deviation sx (equation (11)) as
Vx = sx /mx
(18)
The cases above assume that the variable X is distributed normally. What if X is
distributed log-normally?
When the parameter X is distributed log-normally, we use the transformation ln X=Y
to obtain the variable Y which is distributed normally with N(Y,Y2) (see also Appendix A to
the handbook H1). The relationships between the mean and variance of both variables are
Y = ln(x2)/(x2+x2)
(19)
Y 2= ln(1+x2/x2)=ln(1+Vx2)
(20)
IV - 6
(21)
(22)
(23)
VY from the equation (18), kn from Table 2 and then we calculate Yk from equation (15) using
mY and VY instead of mx and Vx.
Finally we transform the computed characteristic value Yk = mY- k Y (or mY- k sY ) of
the variable Y to the characteristic value Xk of the original variable X:
Xk = exp(mY- k Y)
(24)
Xk = exp(mY- k sY)
(25)
We see that for both types of distribution the design value is greater for the case Vx.
unknown then for the case Vx. known. The reason for this is that we have assumed much
greater Vx.in the case Vx. known that it really is. If we had assumed Vx.=0,015, as obtained
from all the test results, then we would get:
Normal distribution, Vx. known:
m(d) = 0,8 942(1 1,800,015)/1,1=666,6 MPa.
Log-normal distribution, Vx. known, sY = (ln(1+0,0152))1/2=0,015
Yk= 6,85-1,80,015=6,82 ln(MPa)
m(d) = 0,8 exp(6,82)/1,1= 668,2 MPa.
2.3
(26)
In case this method is used, the relevant limit states and the required level or reliability should
be accounted for. The conversion factor d should cover all uncertainties not covered by the
test. The factor kd,n is obtained from the prediction method of fractile estimation with the
lower value of about 0,1 % (the probability p = 0,001). When the coefficient of variation Vx is
known, then parameter X is assumed to be distributed normally and value kd,n is
kd,n = -up (1/n + 1)1/2
(27)
with the 0,001 fractile of the standardized normal distribution up = -3,09 (the value 3,04 is
used in Eurocode 1990 [1]). When the coefficient of variation Vx is unknown, then a fractile tp
of the t-distribution with the degree of freedom n-1 corresponding to probability p=0,001 is
used instead of up. The factor kd,n is dependent on the number of samples n and is given in
Table 3:
Table 3. Values of kd,n for the direct assessment of the design value.
n
1
2
3
4
5
6
8
10
20
Vx
4,36 3,77 3,56 3,44
3,37 3,33 3,27 3,23 3,16
known
Vx
11,40 7,85 6,36 5,07 4,51 3,64
unknown
30
3,13
3,04
3,44
3,04
Example 6.
With the data from the example 5 calculate the design value of tensile strength using
the direct method.
We assume first the normal distribution of the tensile strength. For the case Vx.
known (Vx.=0,015), we have from the Table 3 kd,n=3,37 and the design value is
m(d) = 0,8 942(1 3,370,015)=715,5 MPa.
For the case Vx. unknown we have from the Table 2 kn=7,85 and the design value is
m(d) = 0,8 942(1 7,850,01973)=636,9 MPa
When log-normal distribution is assumed, for the case Vx. known (Vx.=0,015), we
have
Yk= 6,85-3,370,015=6,80 ln(MPa)
m(d) = 0,8 exp(6,80)= 717,9 MPa.
And for the case Vx. unknown
Yk= 6,85-7,850,01967= 6,70 ln(MPa)
IV - 8
We obtained higher values with the direct method of calculation. This is because we
used the same conversion factor d in both cases.
2.4
The procedures given in this section are intended for the calibration of resistance
models and for the derivation of design values from the tests undertaken to reduce
uncertainties in parameters of the resistance model.
Based on observations and theoretical considerations, a design model of the resistance
is developed. The statistical interpretation of the test results should then be used to validate
and adjust the model, until sufficient correlation between test and theoretical data is achieved.
As in the previous section two methods are considered: (a) by assessing the characteristic
value of resistance and (b) by directly assessing the design value of resistance. We start with
the method (a).
The following assumptions are made: the resistance function is a function of
statistically independent variables X=(X1,..,Xj), which are either normally or log-normally
IV - 9
(28)
The model should include all relevant basic variables Xi that affect the resistance. We then
compare the theoretical model with experimental results. Theoretical values rti are calculated
by substituting actual measured properties of the sample i in the theoretical model. This is to
be compared with measured resistance values rei . We plot the points (rti , rei) in a twodimensional diagram with rti on the abscissa and rti on the ordinate, as shown in Figure 4. If
the theoretical model is accurate, all the points should lie on the diagonal of the first quadrant.
Some scatter will always be present in realistic situations, but if any considerable deviation
from that line occurs, further investigation of the experimental procedures and theoretical
model should be done.
re
re = b rt
rt
Figure 4. The diagram of the experimental to theoretical resistance.
Next, represent the probabilistic model of the resistance as:
rp = b grt(X)
(29)
(30)
(31)
In absence of other information it is assumed that >0 and that it is distributed log-normally.
It means that =ln() is distributed normally. For each test value i we calculate
i= rei / (brti)
(32)
i= ln(i)
(33)
and
The estimated mean and variance s for are
= (i)/n
(34)
s2 = (i- )2 / (n-1)
(35)
IV - 10
(36)
(37)
which are considered to be independent and distributed normally, then by taking logarithm of
the above expression
ln(X) = ln(X1) + ln(X2) + ln(X3) ++ ln(Xj )
(38)
we have a sum of log-normally distributed variables. The variance of the sum of independent
variables is (see equation (2))
(39)
(40)
(41)
If the resistance function is more complex function, so that it cannot be expressed as a product
of basic variables, then Vrt is computed from the equation
Vrt2 = D[grt(X)]/grt2(Xm) =
g
1
( rt i )2
g (Xm)
X i
2
rt
(42)
(43)
(44)
Q = ln() = ln(1+V2)
(45)
Q = ln(r) = ln(1+Vr2)
(46)
rt = Qrt/Q
(47)
= Q /Q
(48)
The mean value of the resistance function is obtained from the theoretical model using
the mean values Xm of the basic variables:
rm = b grt(Xm)
IV - 11
(49)
(50)
The factor kn is taken from the Table 2 for the case "Vx unknown" and k is the value
of kn for large n (k =1,64). When only one variable is present in the resistance model, only
the term - kn Q from the equation (50) is taken into account.
When the design value of the resistance is assessed directly, rather that from the
characteristic value of resistance, then the procedure is the same with the only modification
that the values k and kn are replaced with the values k,d and kd,n, taken from the Table 3 for
the case " Vx unknown".
Example 7.
Consider the resistance model in the form F=Am /1000. As the basic variables we
take the section of a rod A and the tensile strength m . Let the mean and coefficient of
variation of the basic variables be E(A)=269,76 mm2, V(A)=0,00295, E(m)=936.5 MPa,
V(m)=0,01509. The theoretical Ft and experimental Fe data for various values of the basic
variables are given in the Table 4. Compute the characteristic value of the resistance at the
mean values of basic variables.
n
1
2
3
4
5
6
Fe
The diagram in Figure 5 shows the data from the Table 4 with the regression line
Ft=bFe. The value b is obtained from equation (30): b=0.9995.
300
290
280
270
260
250
240
230
220
210
200
200 210 220 230 240 250 260 270 280 290 300
Ft
Figure 6. The characteristic value of the resistance for various coefficients of variation
of the basic variables A and m .
If the design value is calculated directly, then k,d =3,04 and kn,d = 3,7223 by the
approximation formula, and the design value is
rd=252,5exp(-3,040,67850,01537-3,72230,73460,0166 -0,50.02262)=233,7 kN
The partial factor m for the model resistance in this case is:
(51)
Two software products for the calculation of the characteristic resistance, design
resistance and partial factors are given in the attachments. These products include the
following.
IV - 13
IV - 14
(A.1)
Xk = exp(lnX kn lnX)
(A.2)
or
Since the mean (X) of X can be expressed with the mean lnX and standard deviation lnX of
ln(X) by the relationship:
(X)=exp(lnX+lnX2/2)
(A.3)
(A.4)
(A.5)
lnX = (ln(Y))+(ln(Z))
(A.6)
(A.7)
If we now combine equations (A.7) and (A.4), then the characteristic value of X as a
combination of Y and Z is:
Xk=(X) exp(-kn Y lnY kn Z lnZ lnX2/2)
(A.8)
(A.9)
Finally, since we assume that there is no statistical uncertainty for theoretical resistance
function rt with respect to the number of samples n, we can substitute k for kn and we obtain
the equation (50).
IV - 15
ATTACHMENTS
1. EXCEL worksheet from the workbook dast.xls
2. Source file dast.c to the program dast.exe
3. Input file dast.i1 for the program dast.exe
4. Output file dast.o1 produced by dast.exe from dast.o1
5. Input file dast.i2 for the program dast.exe
6. Output file dast.o2 produced by dast.exe from dast.o2
7. EXCEL chart showing the data from dast.o2
IV - 16
IV - 17
IV - 18
variations
Vx[1]
Vx[2]
...
Vx[j]
table
Xa Xb Va1 Va2 Vb1 Vb2 na nb
calc
characteristic OR design OR partial
end
-----------------------------------------------------------------------------*/
read_data(s) char *s; { FILE *fd; char buf[501];
if ( (fd=fopen(s,"r"))==0 ) printf("can't open %s!",s), exit(1);
for ( ; fgets(buf,500,fd); ) {
LINE:
if ( buf[0]==';' || buf[0]=='\n' ) continue;
else if ( !strncmp("resistance",buf,10) ) {
for ( ; fgets(buf,500,fd); ) {
if ( n>=N_n ) printf("too many points!"), exit(1);
if ( sscanf(buf,"%lg %lg",Rtheor+n,Rexper+n)!=2 ) goto LINE;
n++;
}
}
else if ( !strncmp("variations",buf,10) ) {
for ( ; fgets(buf,500,fd); ) {
if ( nj>=N_j ) printf("too many variables!"), exit(1);
if ( sscanf(buf,"%lg",Vx+nj)!=1 ) goto LINE;
nj++;
}
}
else if ( !strncmp("table",buf,5) ) {
fgets(buf,500,fd);
if ( sscanf(buf,"%d %d %lg %lg %lg %lg %d %d",
&Xa,&Xb,&Va1,&Va2,&Vb1,&Vb2,&na,&nb)!=8 )
printf("error in 'table'!"), exit(1);
table=1;
}
else if ( !strncmp("calculate",buf,9) ) {
fgets(buf,500,fd);
if
( !strncmp("design",buf,6) ) calc=design;
else if ( !strncmp("partial",buf,7) ) calc=partial;
else if ( !strncmp("characteristic",buf,14) ) calc=characteristic;
else goto ERROR;
}
else if ( !strncmp("end",buf,3) ) break;
else ERROR: printf("error in file: '%s'",buf), exit(1);
}
fclose(fd);
if ( n<=0 || nj<=0 || nj>n ) printf("error in n/nj!"), exit(1);
if ( table && (Xa<=0 || Xb<=0 || Xa>nj || Xb>nj) )
printf("error in table!"), exit(1);
}
/*----------------------------------------------------------------------------Computes the factor of a value of model resistance according to EN 1990,
Anex D, 8.2.2 and 8.2.3. Depending on the value of calc (characteristic,
design, partial) returns the exponent from equation D.17a, D.21 and the ratio
of the two, respectfully.
If table!=0 prints intermediate and final calculatuions.
-----------------------------------------------------------------------------*/
IV - 19
IV - 20
IV - 21
IV - 22
IV - 23
0
1.033
1.08
1.154
1.235
1.322
1.414
1.51
1.611
1.716
1.826
1.939
0.05
1.08
1.109
1.173
1.249
1.333
1.423
1.519
1.619
1.724
1.833
1.946
0.1
1.154
1.173
1.221
1.288
1.366
1.452
1.545
1.644
1.747
1.855
1.967
0.15
1.235
1.249
1.288
1.346
1.417
1.498
1.587
1.683
1.785
1.891
2.002
0.2
1.322
1.333
1.366
1.417
1.482
1.558
1.643
1.736
1.835
1.94
2.05
0.25
1.414
1.423
1.452
1.498
1.558
1.63
1.711
1.801
1.898
2.001
2.109
0.3
1.51
1.519
1.545
1.587
1.643
1.711
1.79
1.877
1.972
2.073
2.18
0.35
1.611
1.619
1.644
1.683
1.736
1.801
1.877
1.962
2.054
2.154
2.26
2,6
2,4
2,2
2
1,8
1,6
1,4
0,4
1,2
0,4
0,3
0,5
Vx1
0,2
0,2
0,1
1
0
IV - 24
Vx2
0.4
1.716
1.724
1.747
1.785
1.835
1.898
1.972
2.054
2.146
2.244
2.349
0.45
1.826
1.833
1.855
1.891
1.94
2.001
2.073
2.154
2.244
2.341
2.446
0.5
1.939
1.946
1.967
2.002
2.05
2.109
2.18
2.26
2.349
2.446
2.549
Summary
The approach to assessment of an existing structure is in many aspects different from
that taken in designing the structure of a newly proposed building. The effects of the
construction process and subsequent life of the structure, during which it may have undergone
alteration, deterioration, misuse, and other changes to its as-built (as-designed) state, must be
taken into account. To assess existing structures general principles and rules of the Eurocode
EN 1990 Basis of structural design must be supplemented by specific procedures provided by
International Standards ISO.
1
INTRODUCTION
1.1
Background documents
Background documents related to assessment of existing structures are limited to few
national codes and three International Standards ISO 2394 [1], ISO 13822 [2] and ISO 12491
[3]. General principles and rules of the Eurocode EN 1990 Basis of structural design [4] must
be supplemented by specific procedures provided in the above mentioned International
Standards ISO [1,2,3] that are primarily used in this contribution. Additional information
concerning assessment of existing structures may be found in scientific papers and
publications, for example in the publications [5], [6] and [7].
1.2
General principles
Assessment of existing structures is becoming a more and more important and
frequent engineering task. Continued use of existing structures is of a great significance due to
environmental, economic and socio-political assets, growing larger every year. These aspects
are particularly relevant to buildings that always constitute a great social and economic value.
General principles of sustainable development regularly lead to the need for extension
of the life of a structure, in majority of practical cases in conjunction with severe economic
constraints. That is why assessment of existing structures often requires application of
sophisticated methods, as a rule beyond the scope of traditional design codes.
The approach to assessment of an existing structure is in many aspects different from
that taken in designing the structure of a newly proposed building. The effects of the
construction process and subsequent life of the structure, during which it may have undergone
alteration, deterioration, misuse, and other changes to its as-built (as-designed) state, must be
taken into account. However, even though the existing building may be investigated several
times, some uncertainty in behaviour of the basic variables shall always remain. Therefore,
similarly as in design of new structures, actual variation in the basic variables describing
actions, material properties, geometric data and model uncertainties are taken into account by
partial factors or other code provisions.
V-1
2.1
Common rules
Two common rules are usually accepted when assessing existing structures:
- Currently valid codes for verification of structural reliability should be considered,
historic codes valid in the period when the structure was designed should be used as
guidance documents only.
- Actual characteristics of structural materials, actions, geometric data and structural
behaviour should be considered, the original design documentation including drawings
should be used as guidance documents only.
The first rule should be applied in order to achieve similar reliability level as in case of
newly designed structures. The second principle should avoid negligence of any structural
condition that may affect actual reliability (in favourable or unfavourable way) of a given
structure.
Most of the current codes are developed assuming the concept of limit states in
conjunction with the partial factor method. In accordance with this method, which is mostly
considered here, basic variables are specified by characteristic or representative values. The
design values of the basic variables are determined on the basis of the characteristic
(representative) values and appropriate partial factors.
It follows from the second principle that a visual inspection of the assessed structure
should be made whenever possible. Practical experience shows that inspection of the site is
also useful to obtain a good feel for actual situation and state of the structure.
As a rule the assessment need not to be performed for those parts of the existing
structure that will not be affected by structural changes, rehabilitation, repair, change in use or
which are not obviously damaged or are not suspected of having insufficient reliability.
2.3
General procedure
In general, the assessment procedure consists of the following steps (see the flow chart
in Annex A to this Chapter):
- specification of the assessment objectives required by the client or authority;
- scenarios related to structural conditions and actions;
- preliminary assessment:
- study of available documentation;
- preliminary inspection;
- preliminary checks;
V-2
INVESTIGATION
3.1
Purpose
Investigation of an existing structure is intended to verify and update the knowledge
about the present condition (state) of a structure with respect to a number of aspects. Often,
the first impression of the structural condition will be based on visual qualitative
investigation. The description of possible damage of the structure may be presented in verbal
terms like: 'unknown, none, minor, moderate, severe, destructive'. Very often the decision
based on such an observation will be made by experts in purely intuitive way.
A better judgement of the structural condition can be made on the basis of
(subsequent) quantitative inspections. Typically, assessment of existing structures is a cyclic
process when the first inspection is supplemented by subsequent investigations. The purpose
of the subsequent investigations is to obtain a better feel for the actual structural condition
(particularly in the case of damage) and to verify information required for determination of
the characteristic and representative values of all basic variables. For all inspection
techniques, information on the probability of detecting damages if present, and the accuracy
of the results should be given.
3.2
Statement
Results of an investigation should be included in the statement that usually contains
the data describing
- actual state of the structure;
- types of structural materials and soils;
- observed damages;
- actions including environmental effects;
- available design documentation.
A proof loading is a special type of investigation. Based on such tests one may draw
conclusions with respect to:
- the bearing capacity of the tested member under the test load condition;
- other members;
- other load conditions;
- the behaviour of the system.
V-3
BASIC VARIABLES
4.1
General
In accordance with the above-mentioned general principles and rules, characteristic
and representative values of all basic variables shall be determined taking into account the
actual situation and state of the structure. Available design documentation is used as a
guidance material only. Actual state of the structure should be verified by its inspection to an
adequate extent. If appropriate, destructive or non-destructive inspections should be
performed and evaluated using statistical methods.
4.2
Characteristic values
For verification of the structural reliability using partial factor method, the
characteristic and representative values of basic variables shall be considered as follows:
(a) Dimensions of the structural elements shall be determined on the basis of
adequate measurements. However, when the original design documentation is
available and no changes in dimensions have taken place, the nominal dimensions
given in the documentation may be used in the analysis.
(b) Load characteristics shall be introduced with the values corresponding with the
actual situation verified by destructive or non-destructive inspections. When some
loads have been reduced or removed completely, the representative values can be
reduced or appropriate partial factors can be adjusted. When overloading has been
observed in the past it may be appropriate to increase adequately representative
values.
(c) Material properties shall be considered according to the actual state of the
structure verified by destructive or non-destructive inspections. When the original
design documentation is available and no serious deterioration, design errors or
construction errors are suspected, the characteristic values given in original design
may be used.
(d) Model uncertainties shall be considered in the same way as in design stage unless
previous structural behaviour (especially damage) indicates otherwise. In some
cases model factors, coefficients and other design assumptions may be established
from measurements on the existing structure (e.g. wind pressure coefficient,
effective width values, etc.).
Thus reliability verification of an existing structure should be backed up by inspection
of the structure including collection of appropriate data. Evaluation of prior information and
its updating using newly obtained measurements is one of the most important steps of the
assessment.
V-4
5.1
Updating in general
Using results of an investigation (qualitative inspection, calculations, quantitative
inspection, proof loading) the properties and reliability estimates of the structure may be
updated. Two different procedures can be distinguished:
(1) Updating of the structural failure probability.
(2) Updating of the probability distributions of basic variables.
Direct updating of the structural reliability (procedure (1)) can be formally carried out
using the following basic formula of probability theory:
P(F|I) =
P( F I )
P( I )
(1)
(2)
where fX(x|I) denotes the updated probability density function of X, fX(x) denotes the
probability density function of X before updating, X a basic variable or statistical parameter, I
inspection information, C normalising constant, and P(I|x) likelihood function. An illustration
of equation (2) is presented in Figure 1.
fX(x), fX(x|I)
X
prior xd
updated xd
V-5
P(F|I) =
(3)
( x | I )dx
g( x <0)
where fX(x|I) denotes the updated probability density function and g(x) < 0 denotes the failure
domain (g(x) being the limit state function). It should be proved that the probability P(F|I),
given the design values for its basic variables, does not exceed a specified target value.
5.4
xd = exp 0,5
(4)
(5)
where xd is the updated design value for X, updated mean value, probabilistic influence
coefficient, target reliability index, V updated coefficient of variation, and 2 = ln(1+V2).
The value of the target reliability index is discussed in ISO/CD 13822 [2], the values
of can be taken equal to those commonly used for new structures (0,7 for the dominating
load parameter, 0,8 for the dominating resistance parameter and 0,3 for non-dominating
variables according to ISO 2394 [1]).
Alternatively one might determine the characteristic value xk first and then calculate
the design value xd by applying the appropriate partial factor m:
xd = xk /m
(6)
For normal and lognormal random variable X the characteristic value xk then follows as
xk = (1 kV )
xk = exp k 0,5
(7)
(8)
where k = 1,64 (5% fractile of the standardised normal distribution) is usually used. It may be
helpful to consider both methods and to use the most conservative result.
V-6
m =
xk
k
=
xd
(9)
which follows from general relationship (6). All the symbols used in (9) are defined above (k
= 1,64 is usually used for the characteristic strength). Similar relationships between m and
may be derived for lognormal or other distributions.
6
STRUCTURAL ANALYSIS
Structural behaviour should be analysed using models that describe actual situation
and state of an existing structure. Generally the structure should be analysed for ultimate limit
states and serviceability limit states using basic variables and taking into account relevant
deterioration processes.
All basic variables describing actions, material properties, load and model
uncertainties should be considered as mentioned above. The uncertainty associated with the
validity and accuracy of the models should be considered during assessment, either by
adopting appropriate factors in deterministic verifications or by introducing probabilistic
model factors in reliability analysis.
When an existing structure is analysed, conversion factors reflecting the influence of
shape and size effect of specimens, temperature, moisture, duration-of-load effect, etc., should
be taken into account. The level of knowledge about the condition of components should be
also considered. This can be achieved by adjusting the assumed variability in either the load
carrying capacity of the components or the dimensions of their cross sections, depending on
the type of structure.
When deterioration of an existing structure is observed, the deterioration mechanisms
shall be identified and a deterioration model predicting the future performance of the structure
shall be determined on the basis of theoretical or experimental investigation, inspection, and
experience.
7
VERIFICATION
V-8
The final report on structural assessment and possible interim reports (if required)
should include clear conclusions with regard to the objective of the assessment based on
careful reliability assessment and cost of repair or upgrading. The report shall be concise and
clear. A recommended report format is indicated in Annex G to ISO 13822 [2].
If the reliability of an existing structure is sufficient, no action is required. If an
assessment shows that the reliability of a structure is insufficient, appropriate interventions
should be proposed. Temporary intervention may be recommended and proposed by the
engineer if required immediately. The engineer should indicate a preferred solution as a
logical follow-up to the whole assessment in every case.
It should be noted that the client in collaboration with the relevant authority should
make the final decision on possible interventions, based on engineering assessment and
recommendations. The engineer performing the assessment might have, however, the legal
duty to inform the relevant authority if the client does not respond in a reasonable time.
10
CONCLUDING REMARKS
V-9
V - 10
Sufficient reliability?
Yes
No
Intervention
Construction
- Repair
- Upgrading
- Demolition
Operation
- Monitoring
- Change in use
V - 11
Periodical inspection
Maintenance
(B.1)
For example consider a simply supported steel beam of the span L exposed to
permanent uniform load g and variable load q. The beam has the plastic section modulus W
and the steel strength fy.
Using the partial factor method the design condition Rd Sd > 0 between the design
value Rd of the resistance R and design value Sd of the load effect S may be written as
W fyk /m - (g gkL2/8 + q qkL2/8) > 0
(B.2)
where fyk denotes the characteristic strength, gk the characteristic (nominal) value of
permanent load g, qk the characteristic (nominal) value of variable load q, m partial factor of
the steel strength, g the partial factor of permanent load and q the partial factor of variable
load. By analogy to (B.2) the limit state function G(X) follows as
G(X) = R - S = W fy - (gL2/8 + qL2/8)
(B.3)
where all the basic variables are generally considered as random variables described by
appropriate probabilistic models.
To verify its reliability the beam has been investigated and a proof loading up to the
level qtest is carried out. It is assumed that gact is the actual value of the permanent load g. If
the beam resistance is sufficient the information I obtained is described as
I = {H > 0} = {W fy - (gact L2/8 + qtest L2/8) > 0 }
(B.4)
where fy is the actual steel strength, gact the actual permanent load assuming it has been
determined (using non-destructive methods) reasonably accurately.
To determine the updated probability of failure P(F| I) using equation (B.1) it is
necessary to assess the following two probabilities:
P(G(X) < 0 H > 0)=P(W fy -(gL2/8 + qL2/8) < 0 W fy-(gact L2/8 + qtest L2/8) > 0) (B.5)
P(H > 0) = P(W fy - (gact L2/8 + qtest L2/8) > 0)
(B.6)
Additional assumptions concerning the basic variables are needed. Having the results of (B.5)
and (B.6) the updated probability of failure P(G(X) < 0| H > 0) follows from (B.1).
Alternatively, considering results of the proof test, the probability density function
fR(r) of the beam resistance R = Wfy may be truncated below the proof load as indicated in
Figure B.1.
V - 12
fR(r)
updated resistance R
prior resistance R
r
resistance adequate to proof load
Figure B.1 Truncated effect of proof loading on structural resistance.
Obviously, the truncation of structural resistance R decreases the updated probability
of structural failure defined as
Pf = P(R S < 0)
and increase, therefore, the updated value of structural reliability.
V - 13
(B.7)
( s ) + n ( m )
2 2
where C is the normalising constant, (n') = 0 for n' = 0 and (n') = 1 otherwise. The prior
parameters m', s', n', ' are parameters asymptotically given as
1
s
E() = m', E() = s', V() =
, V ( ) =
(C.2)
2
m n
while the parameters n' and ' are independent and may be chosen arbitrarily (it does not hold
that ' = n' 1). In equation (C.2) E(.) denotes the expectation and V(.) the coefficient of
variation of the variable in brackets. Equations (C.2) may be used to make estimates for
unknown parameters n' and ' provided the values V() and V() are estimated using
experimental data or available experience.
The posterior distribution function "(,) of and is of the same type as the prior
distribution function, but with parameters m", s", n" and ", given as
n" = n' + n
" = ' + + ( n )
m"n"= n'm' + nm
"(s")2 + n"(m")2 = '(s')2 + n'(m')2 + s2 + nm2
(C.3)
where m and s are the sample mean and standard deviation, n is the size of the observed
sample and = n 1. The predictive value xp,pred of a fractile xp is then
x p , Bayes = m + t p s 1 + 1 / n
(C.4)
where tp is the fractile of the t-distribution (see Table C.1) with " degrees of freedom. If no
prior information is available, then n'= '= 0 and the characteristics m", n", s", " equal the
sample characteristics m, n, s, . Then equation (C.4) formally reduces to so called prediction
estimates of the fractile given as
x p ,pred = m + t p s 1 + 1 / n
(C.5)
where tp denotes again the fractile of the t-distribution (Table C.1) with degrees of freedom.
Furthermore, if the standard deviation is known (from the past experience), then = and
s shall be replaced by .
Example
A sample of n = 5 concrete strength measurements having the mean m = 29,2 MPa and standard
deviation s = 4,6 MPa is to be used to assess the characteristic value of the concrete strength fck
= xp, where p = 0,05. If no prior information is available, then n'= '= 0 and the
characteristics m", n", s", " equal the sample characteristics m, n, s, . The predictive value
of xp then follows from (C.5) as
V - 14
1
+ 1 = 15,5 MPa
6
where the value tp = 2.13 is taken from Table C.1 for 1 p = 0.95 and = 5 1 = 4. When
information from previous production is available the Bayesian approach can be effectively used.
Assume the following prior information
m = 30,1 MPa, V(m) = 0,50, s = 4,4 MPa, V(s) = 0,28
It follows from equation (C.2)
2
1 1
4,6 1
n =
6
<1 , =
2 0,282
30,1 0,50
The following characteristics are therefore considered: n' = 0 and ' = 6. Taking into account that
1
+ 1 = 15,5 MPa
6
where the value tp = 1.81 is taken from Table C.1 for 1 p = 0.95 and = 10.
In this example the resulting characteristic strength is greater (by about 10 %) than the
value obtained by prediction method without using prior information. Thus, when previous
information is available the Bayesian approach may improve (not always) the fractile estimate,
particularly in the case of a great variance of the variable. In any case, however, due caution
should be paid to the origin of the prior information with regard to the nature of considered
variable.
Table C.1 - Fractiles tp of the t-distribution with
1p
V - 15
degrees of freedom
1p
0,90 0,95 0,975
1,36 1,78 2,18
1,35 1,76 2,14
1,34 1,75 2,12
1,33 1,73 2,10
1,32 1,72 2,09
1,32 1,71 2,06
1,31 1,70 2,04
1,28 1,64 1,96
0,99
2,68
2,62
2,58
2,55
2,53
2,49
2,46
2,33
0,995
3,06
2,98
2,92
2,88
2,85
2,79
2,75
2,58
V - 16
where p' ~ P(Bi ) denotes prior probabilities, l~ P(A|B i ) likelihoods and p''~ P(Bi ) updated
(posterior) probabilities.
1 Reading data for apriory probabilties and likelihoods from file
Check vaue
p := DATA
Likelihoods
l := DATA
n := length ( p )
n=2
pp :=
( p l)
pp =
p l
0.889
0.111
p := ( 0.8 0.2 0 0 0 0 0 0 0 0 )
Likelihoods l i
l := ( 1. 0.5 1 1 1 1 1 1 1 1 )
pp :=
( p l)
p l
V - 17
0.889
0.111
pp = 4
pp = 0.889
0
pp = 0.111
1
m1 := 30.1
s1 := 4.4
Current observation:
m := 29.2
s := 4.6
n := 5
:= n 1
m1 V
n1( V) := floor
Size n1
n1( 0.5) = 0
2
2 V
Degrees of freedom 1
Updated size n2
Updated 2
s1
2( V , V) :=
1( V) := floor
1( 0.28) = 6
n2( V) := n + n1( V)
n2( 0.5) = 5
+ 1( V) 1 if n1( V) 1
2( 0.5, 0.28) = 10
+ 1( V) otherwise
4. Updated means and standard deviations
m2( V) :=
m n + m1 n1( V)
s2 ( V , V) :=
n2( V)
2( V , V)
p fractile V unknown
qt inverse Student's
distribution
p fractile V known
qnorm is inverse
normal distribution
ks ( V , V) := qt ( p , 2( V , V) ) 1 +
k( V) := qnorm( p , 0 , 1) 1 +
V - 18
n2( V)
n2( V)
p := 0.05
6. Fractile estimates
Standard deviation s2 unknown
xp( V , V) := m2( V) + ks ( V , V) s2 ( V , V)
24
22
xp( 0.5 , V)
20
xp( 0.1 , V)
18
16
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
xp ( 0.5) = 21.126
24
22
xp( V , 0.28)
20
xp( V , 0.8)
18
16
0
0.1
0.2
0.3
0.4
V
0.5
0.6
V - 19
0.7
0.8
n1( 0.5) = 0
n1( 0.1) = 2
20
15
n1( V)
n2( V)
10
2( V , 0.28)
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
s2 ( 0.5 , V)
s2 ( 0.1 , V)
4
0.1
0.2
0.3
0.4
V
V - 20
0.5
0.6
0.7
0.8
Summary
Traditional methods for designing of civil engineering structures and other engineering
systems are frequently based on the concept of target probability of failure. However, this
fundamental quantity is usually specified on the basis of comparative studies and past
experience only. Moreover, probabilistic design methods suffer from several deficiencies,
including lack of consideration for accidental and other hazard situations and their
consequences. Both of these extreme conditions are more and more frequently becoming
causes of serious failures and other adverse events. Available experience clearly indicates that
probabilistic design procedures may be efficiently supplemented by a risk analysis and
assessment, which can take into account various consequences of unfavourable events. It is
therefore anticipated that in addition to traditional probabilistic concepts the methods of
advanced engineering design will also commonly include criteria for acceptable risks.
1
INTRODUCTION
1.1
Background documents
Background documents of the risk analysis of civil engineering systems considered in
this contribution consist of a number of national and international documents [1] to [9]. It
should be noted that Eurocode EN 1990 [10] for design of civil structures is based on the
concept of the target probability of failure pd. However, it is well recognised that the reliability
of structures and other engineering systems suffers from a number of uncertainties that can
hardly be analysed and well described by probabilistic methods [11, 12]. Moreover, traditional
probabilistic concepts consider the significance of failure and other adverse events only very
vaguely [10]. That is why probabilistic methods are often supplemented by recently developing
methods of risk assessment [12]. In some countries, risk assessment even becomes compulsory
by law in the case of complex technical systems (power stations, tunnel routes).
1.2
General principles
General principles of the risk analysis and the common tools applied for investigating
civil engineering systems considered in this contribution follow the basic concepts presented in
documents [1] to [9]. The risk analysis is an important part of the risk assessment and the entire
risk management of a system as indicated in Figure 1 (adopted from [2]).
The risk analysis of a system consists of the use of all available information to estimate
the risk to individuals or populations, property or the environment, from identified hazards.
Risk assessment further includes risk evaluation (acceptance or treatment) as indicated in
Figure 1 (adopted from [2]). The whole procedure of the risk assessment is typically an
iterative process as indicated in Figure 2 (adopted from [9]). The first step in the risk analysis
involves the context (scope) definition related to the system and the subsequent identification
of hazards.
VI - 1
Risk management
Risk control
Risk assessment
Risk
analysis
Hazard
identification
Risk evaluation
Risk
estimation
Decisionmaking
Risk
acceptance
Monitoring
Risk
treatment
HAZARD IDENTIFICATION
A hazard is a set of circumstances, possibly occurring within a given system, with the
potential for causing events with undesirable consequences. For instance the hazard of a civil
engineering system may be a set of circumstances with the potential to an abnormal action (e.g.
fire, explosion) or environmental influence (flooding, tornado) and/or insufficient strength or
resistance or excessive deviation from intended dimensions. In the case of a chemical
substance, the hazard may be a set of circumstances likely to cause its exposure [2].
Hazard identification and modelling is a process to recognize the hazard and to define
its characteristics in time and space. In the case of civil engineering systems the hazards Hi may
be linked to various design situations of the building (as defined in [7]) including persistent,
transient and accidental design situation. As a rule Hi are mutually exclusive situations (e.g.
persistent and accidental design situations of a building). Then if the situation Hi occurs with
the probability P{Hi}, it holds P{Hi} = 1. If the situations Hi are not mutually exclusive, then
the analysis becomes more complicated.
VI - 2
Start
Hazard identification
Concequence analysis
Risk estimation
Risk assessment
Probabilty analysis
Risk analysis
Risk evaluation
No
Risk treatment
Acceptable risk?
Yes
Stop
Figure 2. Flowchart of iterative procedure for the risk assessment (adopted from [9]).
Note that in some documents (for example in the recent European document EN 1990
[10]) the hazard is defined as an event, while in risk analysis [2] it is usually considered as a
condition with the potential for causing event, thus as a synonym to danger.
3
Hazard scenario is a sequence of possible events for a given hazard leading to undesired
consequences. To identify what might go wrong with the system or its subsystem is the crucial
task to risk analysis. It requires detail examination and understanding of the system [6].
Nevertheless, a given system is often a part of a larger system. Consequently, modelling and
subsequent analysis of the system is a conditional analysis.
The modelling of relevant scenarios may be dependent on specific characteristics of the
system. For this reason a variety of techniques have been developed for the identification of
hazards (e.g. PHA HAZOP) and for the modelling of relevant scenarios (fault tree, event
tree/decision trees, causal networks). Detail description of these techniques is beyond the scope
of this contribution, may be however found in [6, 9] and other literature.
VI - 3
ESTIMATION OF PROBABILITIES
P{H }P{F | H }
i
(1)
Equation (1) can be used for the modification of the partial probabilities P{Hi}P{F|Hi}
(appropriate to the situations Hi) with the aim to comply with the design condition pF < pt,
where pt is a specified target probability of failure. The target value pt may be determined using
the probabilistic optimisation of an objective function describing, for example, the total cost.
The conditional probabilities P{F|Hi} must be determined by a detail probabilistic
analysis of the respective situations Hi under relevant scenarios. The traditional reliability
methods [8] assume that the failure F of the system can be well defined in the domain of the
vector of basic variables X. For example, it is assumed that a system failure may be defined by
the inequality g(x) < 0, where g(x) is the so called limit state function, where x is a realisation
of the vector X. Note that g(x) = 0 describes the boundary of the limit state, and the inequality
g(x) > 0 the safe state of a structure.
If the joint probability density fX(x|Hi) of basic variables X given situation Hi is known,
the conditional probability of failure P{F|Hi} can be then determined [6] using the integral
P{F|Hi} =
X
g ( x )<0
( x | H i )dx
(2)
It should be mentioned that the probability P{F|Hi} calculated using equation (2) suffers
generally from two essential deficiencies:
- uncertainty in the definition of the limit state function g(x),
- uncertainty in the theoretical model for the density function fX(x|Hi) of basic variables
X [8].
These deficiencies are most likely the causes of the observed discrepancy between the
determined probability pF and actual frequency of failures; this problem is particularly
disturbing in case of fire. Yet, the probability requirement pF < pt is generally accepted as a
basic criterion for design of structures.
In a risk analysis we need to know not only probability of the structural failure F but
probabilities of all events having unfavourable consequences. In general, the situations Hi may
cause a number of events Eij (e.g. excessive deformations, full development of the fire). The
required conditional probabilities P{Eij|Hi} must be estimated by a separate analysis using
various methods, for example the fault tree method or causal networks.
5
ESTIMATION OF CONSEQUENCES
ESTIMATION OF RISK
Risk is a measure of the danger that undesired events represent for humans,
environment or economic values. Risk is commonly expressed in the probability and
consequences of the undesired events. It is often estimated by the mathematical expectation of
the consequences of an undesired event. Then it is the product "probability consequences".
However, a more general interpretation of the risk involves probability and consequences in a
non-product form. This presentation is sometimes useful, particularly when a spectrum of
consequences, with each magnitude having its own probability of occurrence, is considered [2].
The estimation of risk is the process used to produce an estimate of a measure of risk.
As already stated above the risk estimation is based on the hazard identification and generally
contains the following steps: scope definition, frequency analysis, consequence analysis, and
their integration [2]. If there is one-to-one mapping between the consequences Cij,k and the
events Eij, then the risk component Rk related to the considered situations Hi is the sum
Rk =
ij , k
P{Eij | H i }P{H i }
(3)
ij
(4)
where R(x) denotes the degree of risk as a function of basic variables X, and fX(x) denotes joint
probability density function X.
7
LOGIC TREES
A number of different logic (decision) trees (fault tree, event tree, cause/consequence
chart) have been developed to analyse the risk of a system [11] to [13]. Applications of logic
trees significantly improve the completeness and clarity of the engineering work. The use of
this kind of tool is widespread in risk analysis and implies some important advantages.
Influences of the environment and of human activities can easily be considered simultaneously.
Logic trees can also enable the detection of the most effective countermeasures. Furthermore,
they can be easily understood by inexperienced persons and therefore can provide very
effective communication means between experts and public authorities.
The fault tree can be defined as a logical diagram for the representation of combinations
of influences that can lead to an undesired event. When establishing a fault tree, the undesired
VI - 5
OR
H
M1
M2
Sway
Vertical
Combined
AND
AND
AND
2
M2
M4
Sway
2b
M5
4
G+Q
M3
Vertical
M4 H
G+Q
M4
M3
M1
Combined
M5
VI - 6
Fire flashover
Collapse
Yes
0,01
Yes
0,05
Yes
No
0,99
0,00495
0,1
Yes
0,00000095
0,00001
No
0,95
No
0,09499905
0,99999
Design situation
Yes
0,000009
0,00001
No
0,9
No
0,899991
0,99999
Figure 4. Event tree describing the collapse of a structure under persistent and fire design
situation (all data are approximately related to a 50-year period of an administrative building
having the fire compartment area 250 m2 without sprinklers and with a protected steel
structure).
Logic trees may be supplemented by the consequences of events; graphical
representation of such a tree is called the cause/consequence-chart. The consequence chart
corresponds to an event tree with a suitable representation of expected consequences. For
example Figure 4 may include consequences linked to each failure probability (frequency per
year) of the structural collapse under given conditions. Then the tree may be used for the
cause/consequence or risk (utility) analysis.
The simplest form of the cause/consequence consideration is the so-called prior-analysis
of the risk (utility) when the basic statistical and probabilistic information is available prior to
any decision or activity. The prior analysis is an assessment of the risk associated with different
decisions; commonly used for comparing the risks corresponding to different decisions. The
posterior decision analysis differs from the prior analysis by considering possible changes in
the branching probabilities and/or the consequences due to risk reducing measures, risk
mitigating measures and the collection of additional information. The posterior decision
analysis may be used to evaluate different additional activities affecting the total risk.
Other important modification of logic trees is known as the pre-posterior decision
analysis. The aim of the pre-posterior decision analysis is to identify the optimal decisions with
regard to activities that may be performed in the future, e.g. planning of risk reducing activities
and/or the collection of new information. An important pre-requisite for the pre-posterior
decision analysis is the consideration of future actions that may be applied taking into account
the results of the planned activities.
VI - 7
BAYESIAN NETWORK
Another promising tool for the risk analysis seem to be Bayesian (believe) causal
networks [13,14]. A simple example of the causal network is shown in Figure 5. The network
containing only four chance nodes describes the structural failure under persistent and fire
design situation similarly as the event tree in Figure 4. Compared with the event tree shown in
Figure 4 the network in Figure 5 includes also the effect of sprinklers (node B). Note that the
directional arrows in Figure 5 indicate the causal links between interconnected chance nodes.
B - Sprinkler
D-Flashover
A-Situation
E-Collapse
Figure 5. The causal network describing the structural failure under persistent and fire design
situation.
The collapse of a structure depends on the probability of persistent and fire situation
and on the conditional probabilities of full development of fire, that depends on the ability of
sprinklers and on conditional probability of the structural collapse under the conditions given
by parents nodes (for example when fire is fully developed - fire flash over). Obviously the
causal network representation seems to be much more effective than the event tree version.
Moreover each node may have several states. Consequently, the input data are not indicated
directly in the graphical representation of the network but are given in the tables of conditional
probabilities.
The basic principle of probability calculation used in the Bayesian networks may be
illustrated considering the nodes A, B and D of the network in Figure 5. One child node D (Fire
flashover) is dependent on two parent nodes: A (Design situation) and B (Sprinklers). If the
parents' nodes A and B have the discrete states Ai and Bj, then the probability of the event Dk (a
particular state of the node D) is given by the formula
P(Dk) = P(Dk| AiBj)P(Ai) P(Bj)
(5)
Equation (5) represents the fundamental theoretical tool for analysing the Bayesian
network. The input data consist of the probabilities P(Ai) and P(Bj), and the conditional
probabilities P(Dk|AiBj). These extensive data are based on available statistical evidence,
probabilistic analysis or expert assessment (judgement) and are transparently summarised in
the tables of conditional probabilities.
Bayesian networks supplemented by decision and utility nodes called influence
diagrams [13,14] provide a powerful tool for the risk estimation. In fact the influence diagram
is a generalisation of the cause/consequence-chart discussed above. The main features of this
tool are illustrated by the example shown in Figure 6, which is an extension of the fundamental
task indicated in Figure 5. Figure 6 shows a simplified influence diagram developed recently
[15,16] for the risk analysis of buildings under persistent and fire design situation.
VI - 8
1-Situation
8-C8
6-Sprinklers
10-C10
7-Protection
17-C17
2-Sprinklers
4-Protection
3-Flashover
5-Collapse
15-Size
12-Smoke
16-Escape r.
14-Number
9-C9
11-C11
13-C13
Figure 6. The Bayesian network describing a structure under normal and fire design situations.
The chance nodes 1, 2, 3, 4, 5, 12 and 14 represent alternative random variables having
two or more states. The node 1-Situation describes the probability of fire start pfi,s = P(H2) and
the complementary probability 1 pfi,s of normal situation H1. The chance node 2-Sprinklers
describes the functioning of sprinklers provided that the decision (node 6) is positive; the
probability of the active state of the sprinklers given the fire start is assumed to be very high,
for example 0,999. The chance node 3-Flashover has two states: the design situation H3 (fire
design situation without flashover) and H4 (fire design situation with flashover when the fire is
fully developed).
When sprinklers are installed, the flashover in a compartment of 250 m2 has the positive
state with the conditional probability 0,002; if sprinklers are not installed then P{H4|H2} =
0,066 [15,16]. It is assumed that with the probabilities equal to the squares of the above values
the fire will flash over the whole building, thus the values 0,000004 and 0,0044 are considered
for the chance node 3. The chance node 4-Protection (introduced for formal computational
reasons) has identical states as the decision node 7-Protection. The chance node 5-Collapse
represents the structural failure that is described by the probability distribution linked to three
children nodes (1,3,4). This situation can hardly be modelled using a decision tree. Note that
the probability of collapse in the case of fire but not flashover may be smaller than in a
persistent situation, due to the lower imposed load.
VI - 9
DECISION-MAKING
The decision-making is generally based on the process of the risk acceptance and option
analysis (see Figure 1) that is sometimes referred to as the risk evaluation. The risk acceptance
is based on various criteria of risk that are reference points against which the results of the risk
analysis are to be assessed. The criteria are generally based on regulations, standards,
experience, and/or theoretical knowledge used as a basis for the decision about the acceptable
risk. Acceptance criteria and the criteria of risk may be sometimes distinguished [1]. Various
aspects may be considered, including cultural, social, psychological, economical and other
aspect [6], [17], [18] and [19]. Generally the acceptance criteria may be expressed verbally or
numerically [6].
Assuming for example that the acceptance limits Ck,d for the components Ck are
specified, then it is possible to design the structure on the basis of acceptable risks using the
criterion Ck < Ck,d, which may supplement the probability requirement pF < pt.
It should be noted that various levels of risk might be recognized, for example
acceptable risk, tolerable risk, and objective risk [6] (see the definitions of theses terms in [2]).
It is a remarkable fact that the public seems to be generally better prepared to accept certain
risks than to stand for specified probabilities of failure [17].
10
The consequences may generally include economic as well as social and environmental
costs [17,18,19]. An example is provided by the influence diagram shown in Figure 6 used to
assess the risk of a building due to fire. Thus, in order to compare all possible consequences it
is necessary to express all consequences in terms of a single unit. This seems to be an
extremely difficult task. One of the possible approaches is represented by the concept of the
Implied Cost of Averting a Fatality ICAF or Life Quality index LQI [19]. Table 1 shows values
of the cost ICAF for selected countries adopted from [19]. It appears that the cost ICAF may be
estimated to about 1 to 3 million of USD.
Table 1. The Implied Cost of Averting a Fatality ICAF(e)), financial data in PPP US$ (1999)
obtained from UN-HDR 2001, World Bank.
Country
US
Japan
Germany
UK
Czech Republic
Mexico
South Africa
Colombia
China
India
Nigeria
g- annual
income
e- life
time
34000
26000
25000
22000
8000
8800
9100
5900
3900
2400
800
VI - 10
77
81
77
77
75
72
55
70
70
63
47
2 wICAF(e)
working
[ 106]
part of e
0.15
2.6
0.15
2.1
0.125
1.9
0.125
1.7
0.15
0.6
0.15
0.6
0.15
0.5
0.15
0.4
0.15
0.3
0.15
0.1
0.18
0.04
e 1 w
(6)
ICAF (e) = g 1 (1 + ) e
e
where symbols g, e and w are defined in Table 1. However, the concept of the Implied
Cost of Averting a Fatality described by equation (6) is just one of possible approaches to the
complex problem of evaluating social consequences. At present further intensive investigation
is expected.
11
CONCLUDING REMARKS
VI - 11
Reliability index.
X(x) Probability density function of the vector of basic variables X.
N1 ( p F ) Inverse distribution function of a standardized normal variable.
VI - 12
Summary
Elementary concepts and techniques of the theory of probability and mathematical
statistics required for understanding of basic reliability methods are reviewed and illustrated
by a number of numerical examples. Computational procedures for determination of sample
characteristics, fractiles of common theoretical models and estimates for fractiles based on
small samples can be applied using the attached MATHCAD sheets.
1
INTRODUCTION
1.1
Background materials
Elementary concepts and techniques of the theory of probability and mathematical
statistics applicable to civil engineering are available in a number of standards [1 to 5],
background materials [6, 7, 8], software products [9, 10, 11] and books [12 to 24]. Additional
information may be found in the extensive literature listed in the books [12, 13] and others. In
particular, documents developed by JCSS [6, 7] and recently published handbook [8] are
closely related to the statistical techniques described in this text.
1.2
General principles
The theory of structural reliability is based on a general principle that all the basic
variables are considered as random variables having appropriate type of probability
distribution. Different types of distributions should be used for description of actions, material
properties and geometric data. Prior theoretical models of basic variables and procedures for
probabilistic analysis are indicated in JCSS documents. Sample characteristics are used as
estimates of population parameters. In addition the population fractiles must be often assessed
using small samples. MATHCAD sheets that supplement described computational procedures
can be effectively used in practical applications.
2
2.1
General
Actions, mechanical properties and geometric data are generally described by random
variables (mainly by continuous variables). A random variable X, (e.g. concrete strength), is
such a variable, which may take each of the values of a specified set of values (e.g. any value
from a given interval), with a known or estimated probability. As a rule, only a limited
number of observations, constituting a random sample x1, x2, x3,..., xn of size n taken from a
population, is available for a variable X. Population is a general statistical term used for the
totality of units under consideration, e.g. for all concrete produced under specified conditions
within a certain period of time. The aim of statistical methods is to make decisions concerning
A-1
Sample characteristics
A sample characteristic is a quantity used to describe the basic properties of a sample.
The three basic sample characteristics, which are most commonly used in practical
applications, are:
- the mean m representing the basic measure of central tendency;
- the variance s2 describing the basic measure of dispersion; and
- the coefficient of skewness giving the basic measure of asymmetry.
The sample mean m (an estimate of the population mean) is defined as the sum
m = ( xi) / n
(2.1)
with the summation being extended over all the n values of xi.
The sample variance s2 (an estimate of the population variance), is defined as:
s2 = ( (xi - m)2) / (n - 1)
(2.2)
the summation being again extended over all values xi. Sample standard deviation s is the
positive square root of the variance s2.
The sample coefficient of skewness (an estimate of the population skewness)
characterising asymmetry of the distribution is defined as
(2.3)
Thus, the coefficient of skewness is derived from the central moment of order 3
divided by s3. If the sample has more distant values to the right from the mean than to the left,
the distribution is said to be skewed to the right or to have a positive skewness. If the reverse
is true, it is said to be skewed to the left or to have a negative skewness.
In some cases two different samples may be taken from one population and their
combination is needed. If the original data are not available, then the characteristics of
combined sample may be determined using the characteristics of both samples. If the sample
sizes are n1, n2, the means m1, m2, standard deviations s1, s2 and skewnesses 1, 2, then the
combined sample of the size n = n1+n2 has the characteristics
m=
s2 =
n1m1 + n2 m2
n
(2.4)
s 3
n
n2
n3
(2.5)
The coefficient of variation v can be effectively used only if the mean m differs from
zero. When the mean is much less than the standard deviation, then the standard deviation
rather then the coefficient of variation should be considered as a measure of the dispersion.
The coefficient of variation v is often used as a measure of production quality; for concrete
A-2
Distribution function
Probability distribution is a term generally used for any function giving the probability
that a variable X belongs to a given set of values. The basic theoretical models used to
describe the probability distribution of a random variable may be obtained from a random
sample by increasing the sample size or by smoothing either the frequency distribution or the
cumulative frequency polygon.
An idealisation of a cumulative frequency polygon is the distribution function (x)
giving, for each value x, the probability that the variable X is less than or equal to x:
(x) = P (X x)
(2.6)
(2.7)
Note that Appendix 1 to this Chapter provides a review of selected theoretical models
of continuous random variables that are most frequently used in reliability analysis of civil
structures.
Example 2.1.
A continuous random variable, which may attain equally likely any point x within a
two-sided interval <a, b> (each point x has the same probability density (x)) is described by
a so-called uniform distribution shown in Figure 2.1.
1
(x) = (x-a)/(b-a)
Distribution function
x
a
(x) = 1/(b-a)
Probability density
function
x
A-3
(x)dx =
(x)dx = 1
(2.8)
Thus, the surface bounded by the horizontal axis x and the curve of the density function (x)
has the area equal to unity.
2.4
Population parameters
The population parameters are quantities used in describing the distribution of a
random variable, as estimated from one or more samples. As in the case of random samples,
three basic population parameters are commonly used in practical applications:
- the mean representing the basic measure of central tendency;
- the variance 2 as the basic measure of dispersion; and
- the coefficient of skewness .giving the degree of asymmetry.
The population mean , for a continuous variable X having the probability density
(x), is defined as
=
x (x)dx
(2.9)
the integral being extended over the interval of variation of the variable X. The population
variance 2, for a continuous variable X having the probability density function (x), is the
mean of the squared deviation of the variable from its mean:
2 = (x - )2 (x)dx
(2.10)
The population standard deviation is the positive square root of the population
variance 2.
The population coefficient of skewness, characterising asymmetry of the distribution,
is defined as
= (x - )3 (x)dx / 3
(2.11)
Another population parameter based on the fourth order moment is called kurtosis .
= (x - )4 (x)dx / 4 - 3
(2.12)
Note that for normal distribution (described in Section 3.1) the kurtosis defined by
equation (2.12) is zero. However, this parameter is used mainly in theoretical consideration.
Another important parameter of the population is the coefficient of variation V defined
similarly as the sample coefficient of variation
V=/
(2.13)
The same restriction on the practical use of V applies as in the case of samples.
Geometrically is actually the x coordinate of the centre of gravity of the area
bounded by the horizontal axis x and the curve of density function (x). Figure 2.2 shows an
example of probability density function of lognormal distribution illustrating the geometric
interpretation of the mean and standard deviation .
The measure of dispersion of a random variable X relative to the mean is given by
the central moment of the second order (moment of inertia) of the area, and standard deviation
A-4
is therefore the centroidal radius of gyration around the mean of the area bounded by the
(2.14)
In civil engineering the probabilities p = 0,001; 0,01; 0,05 and 0,10 are used most
frequently. The probability p is often written as a percentage (e.g. p = 0,1 %; 1 %; 5 %; 10
%). If this is done, then xp is called a percentile, for example the 5th percentile is used when p
= 5 %. If p=50 %, then xp is called the median. More details about the fractiles of continuous
variables are given in the following sections.
Probability density (x)
2.5
2.0
1.5
=1.0
1.0
0.5
0.0
0.4
0.6
0.8
1.0
1.2
1.4
1.6
x
Figure 2.2. Geometric illustration of the mean and standard deviation
1.8
Example 2.2.
Parameters of the uniform distribution from example 2.1 may be derived using
equations (2.9) to (2.13) as
A-5
3.1
Normal distribution
Most frequently used models of continuous random variables that are applied in
reliability analysis of civil structures are reviewed in Appendix 1 of this Chapter. From a
practical and theoretical point of view the most important type of distribution of a continuous
random variable is the normal (Laplace-Gauss) distribution. Symmetric normal distribution of
a variable X is defined on an unlimited interval - < x < (which can be undesirable in some
practical applications) and depends on two parameters only on the mean and on the
standard deviation . Symbolically it is often denoted as N(,).
The normal distribution is frequently used as a theoretical model of various types of
random variables describing some loads (self-weight), mechanical properties (strengths) and
geometrical properties (outer dimensions). It is convenient for symmetric random variable
with a relatively low variance (coefficient of variation V < 0,3). It fails when used for
asymmetric variables with great variance and skewness > 0,5.
The probability density function of a normal random variable X with a mean and
standard deviation is given by the exponential expression
1 x 2
1
(x) =
exp
2
2
(3.1)
(3.2)
The standardized random variable U has a zero mean and variance (standard
deviation) equal to one; symbolically it is often denoted as N(0, 1).
The probability density function of the standardized random variable U is then given
as a function of u
(u) =
1
2
exp u
2
2
(3.3)
Lognormal distribution
Generally one-sided limited asymmetric lognormal distribution is defined on a limited
interval x0 < x < or - < x < x0. Therefore it eliminates one of the undesirable properties of
the normal distribution. A lognormal distribution is generally dependent on three parameters.
Commonly the moment parameters are used: mean X, standard deviation X and skewness
x. If the skewness x is unknown or uncertain, the lower or upper bound x0 is used.
A-6
(3.4)
has a normal distribution. In this relation x0 denotes the lower or upper limit of distribution of
a variable X, which depends on skewness x. If the variable has a mean x and standard
deviation x, then the lower or upper limit can be expressed as
x0 = X - X /c
(3.5)
where the coefficient c is given by the value of skewness X according to the relation
X = c3 + 3c3
(3.6)
c = 2X + 4 + x
) (
13
2
X
+ 4 X
13
1 3
2
(3.7)
0.4
0.3
0.2
0.1
0.0
-3
-2
-1
0,5
1,0
1,5
2,0
u0= 1/c
-6,05
-3,10
-2,14
-1,68
A-7
1
ln u + + ln c 1 + c 2
c
u =
ln(1 + c 2 )
)
(3.8)
where (as above) u = (x-X)/X denotes the original standardised variable. The probability
density function LN,U(u) and the distribution function LN,U(u) = LN,X(x) of the lognormal
distribution are then given as
LN,U(u) =
(u )
1
u + ln(1 + c 2 )
c
(3.9)
(3.10)
where (u) and (u) denote the probability density and distribution function of the
standardised normal variable.
A special case is the popular lognormal distribution with a lower bound at zero (x0 =
0), which like the normal distribution, depends on two parameters only the mean X and the
standard deviation X (symbolically it is denoted LN(, )). In such a case it follows from
equations (3.5) that the coefficient c is equal to the coefficient of variation VX. It further
follows from equation (3.6) that the skewness X of the lognormal distribution with a lower
bound at zero is given by the value of the coefficient of variation VX as
X = 3V X + V 3X
(3.11)
Thus the lognormal distribution with the lower bound at zero (x0 = 0) always has a
positive skewness, which may have relatively high value (greater than 0,5); e.g. for the
coefficient of variation equal to 0,30 a coefficient of skewness Vx = 0,927 obtained from
relation (3.11). Applications of the lognormal distribution with the lower limit at zero (xo = 0)
can thus lead to unrealistic theoretical models (usually underestimating the occurrence of
negative and overestimating the occurrence of positive deviations from the mean), particularly
for higher values of coefficient of variation VX. Although the occurrence of negative values
can also be undesirable (unrealistic for most mechanical quantities), it is usually negligible
from a practical point of view.
Example 3.1.
Reinforcement cover layer of a reinforced concrete cross-section X has a mean =
25 mm and standard deviation = 10 mm. The probability density function (x) for a normal
distribution and for a lognormal distribution with a lower limit at zero is shown in Figure 3.2.
A-8
0.02
0.01
0.00
-10
10
20
30
40
50
60
70
80
90
Gamma distribution
Another popular type of one-side limited distribution is the type III Pearson
distribution. Its detailed description is e.g. in the book [13]. A special case of the type III
Pearson distribution with lower limit at zero is the gamma distribution. The probability
density function of this important distribution is dependent on two parameters only: on the
mean and standard deviation . To simplify the notation two auxiliary parameters and
k are often used
( x) = x
k
exp(x)
, = 2 ,k =
(k )
k 1
A-9
(3.12)
, =
, =
2
2
=
= 2V
(3.13)
The curve is bell shaped for k > 1, i.e. for skewness < 2 (in the inverse case it is a
decreasing function of x). For k , the gamma distribution approaches the normal
distribution with parameters and .
The gamma distribution is applied similarly as the lognormal distribution with lower
bound at zero. However, it varies from the lognormal distribution by its skewness, which is
equal to twice the coefficient of variation ( = 2V) and is thus lower than the skewness of
lognormal distribution, which is more than 50% higher (according to equation (3.11) it is
X = 3VX + VX3 ). That is the reason why the gamma distribution is more convenient for
describing some geometrical quantities and variable action that do not have a great skewness.
Probability density (x)
0.05
n = 157
m = 26,8
s = 11,1
v = 0,42
a = 0,40
0.04
Log-normal distribution
Gamma distribution
Normal distribution
Beta distribution
0.03
0.02
0.01
0.00
10
20
30
40
Concrete cover [mm]
50
60
70
Figure 3.3. Histogram and theoretical models for concrete cover of reinforcement
Example 3.2.
A sample of the size n = 157 experimental results of concrete cover of reinforcement
measurements has these characteristics: m = 26,8 mm, s = 11,1 mm and v = 0,42. It is a
relatively large sample, which can be used for the assessed skewness (furthermore a long-term
experience is available). A histogram of the obtained values and theoretical models of normal
distribution, lognormal distribution with origin at zero, gamma distribution and beta
distribution are shown in Figure 3.3, with help of which the appropriateness of the individual
models can be considered.
A - 10
Beta distribution
An interesting type of distribution is the so-called beta distribution (also called
Pearsons type I curve), which is defined on a both-side limited interval <a, b > (but this
interval can be arbitrarily extended and the distribution then approaches the normal
distribution). Generally it is dependent on four parameters and it is used mainly in those cases
when it is evident that the domain of the random variable is limited on both sides (some
actions and geometrical data, e.g. weight of a subway car, fire load intensity, concrete cover
of reinforcement in a reinforced concrete cross-section). The principal difficulty in practical
application is the need to estimate all the four parameters, for which credible data may not be
available.
The beta distribution is usually written in the form
( x) =
( x a )c 1 ( x b) d 1
B(c, d )(b a )c + d 1
(3.14)
c + d +1
cd
(3.15)
where g is an auxiliary parameter. From equations (3.15), relations for parameters c and d can
be derived
c=
a ( a )(b )
b ( a )(b )
1 , d =
1
2
ba
ba
(3.16)
(3.17)
2 g (d c)
3g 2 (2(c + d ) 2 + cd (c + d 6))
,=
3
(c + d + 2)
(c + d + 2)(c + d + 3)
(3.18)
=
=
Note that skewness and kurtosis are dependent only on the parameters c and d
(they are independent of the limits a and b). That is why the parameters c and d are called
shape parameters. In practical applications the distribution is used for c > 1 and d > 1
(otherwise the curve is J or U shaped), for c = d = 1 it becomes a uniform distribution, for c =
d = 2 it is the so-called parabolic distribution on the interval <a, b >. When c = d, the curve is
symmetric around the mean. When d , the curve becomes the type III Pearson
distribution (see Section 3.3). If c = d , it approaches the normal distribution. Depending
A - 11
2V
(3.19)
where the coefficient of variation V = / . For = 2V the curve becomes the type III
Pearson distribution (see Section 3.3). Therefore if the input parameters are the mean ,
standard deviation and skewness 2V, the beta distribution with the lower limit at zero (a
= 0) is fully described. The upper limit b of the beta distribution with the lower limit at zero
follows from the relation (3.15)
b=
(c + d )
c
(1 + V (2 + V ))
(2V )
(3.20)
In equation (3.20) the parameters c and d are substituted by the following expressions
c=
d=
(2V ) 2 (4 + 2 )
2V (V + 2) 2 (4 + 2 )
(2V ) 2 (4 + 2 ) 2 + V
2 (V + 2) 2 (4 + 2 ) 2V
(3.21)
(3.22)
d=
The upper bound of the distribution b follows from equation (3.20) that
b=
25 (4,407 + 12,927)
= 98,325
4,407
A - 12
0,04
0,03
0,02
0,01
0,00
-10
10
20
30
40
50
60
70
80
(3.23)
(3.24)
Both the parameters xmod, c of the Gumbel distribution can be assessed from the mean
and standard deviation
xmod = 0,577
c=
(3.25)
(3.26)
(3.27)
(3.28)
so the mean N and standard deviation N of maxima of populations that are N times greater
than the original population are
N = + ln N/c = + 0,78 ln N , N =
(3.29)
Thus the standard deviation N of the greater population is equal to the standard
deviation of the original population, N = , but the mean N is greater than the original value
by 0.78 lnN/c.
Example 3.4.
One-year maxima of wind pressure are described by Gumbel distribution with a mean
1 = 0,35 kN/m2, 1 = 0,06 kN/m2. The corresponding parameters of the 50-year maximum
value distribution, i.e. parameters 50 and 50, follow from equation (3.29)
A - 14
1(x)
50(x)
6
5
4
3
2
1
0
0.2
0.3
0.4
0.5
0.6
0.7
0.8
Wind pressure x
Figure 3.5. Distribution of maximum wind pressure over the periods of 1 year and 50 years.
The distribution function of type I minimal values distribution (Gumbel distribution of
minimum values) has the form
(x) = 1 - exp(-exp(-c(xmod - x)))
(3.30)
(3.31)
Both these parameters can be assessed from the mean and standard deviation
xmod = + 0,577
c=
(3.32)
(3.33)
The probability density function of the minimum values is symmetrical to the shape of
maximal values relative to mode xmod, as it is apparent from Figure 3.6.
A - 15
0,3
Distribution of the
minimum values
0,2
Distribution of the
maximum values
0,1
0,0
10
12
13
15
16
Variable x
Figure 3.6. The Gumbel distribution of the minimum and maximum values.
In a similar way the type II distribution, the so-called Frchet distribution, and type III
distribution, the so-called Weibull distribution, are defined. All the three types of distribution
complement each other with respect to the possible values of skewness . Each type covers a
certain area of skewnesses, as indicated in Figure 3.7.
type I
0
Distribution of the minimum values
type II type I
-1,14
1,14
type II
type III
A - 16
(3.34)
where X = [X1, X2, ... , Xn] denotes a vector of basic variables. Then the resulting variable Z is
a random variable and its characteristics may be derived from relevant characteristics of basic
variables X = [X1, X2, ... , Xn]. Usually three moment parameters, the mean , standard
deviation and skewness , are used for a first assessment of the resulting variable Z.
Experience shows that using derived moment parameters (, and ) three parameter
lognormal distribution provides satisfactory approximation of Z. However, the software VaP
[9] applies a more accurate approximation based on four moment parameters (, , and
kurtosis ).
Appendix 2 of this Chapter provides approximate expressions for fundamental
functions of two basic variables that can be used in assessment of failure probability in case of
small number of basic variables.
4
ESTIMATION OF FRACTILES
4.1
(4.1)
The same definition holds also for standardised random variable U (given by
transformation equation (3.2)) when in equation (4.1) U is substituted for X and up is
substituted for xp. Fractiles up of standardised random variables U are commonly available in
tables. Figure 4.1 illustrates the definition of the fractile described by equation (4.1) for
standardised random variable U; it shows distribution function (u), probability density
function (u), probability p (approximately equal to 0,05) and fractile up for normal
standardised distribution U.
In general fractile xp of the original random variable X may be calculated using tables
for up available for standardised random variables U with a relevant type of distribution. It
follows from transformation (3.2) that the fractile xp may be determined from the fractile of
the standardised random variable up (found in available tables) using relationship
xp = + up = (1 + up V)
(4.2)
where denotes the mean, the standard deviation and V the coefficient of variation
of the observed variable X.
If the probability p < 0,5, then the value xp is called the lower fractile, for p > 0,5 the
xp is called the upper fractile. Figure 4.2 shows the lower and upper fractiles up of a
A - 17
0.8
0.6
0.4
0.2
up
p
0.0
-3.0
-2.5
-2.0
-1.5
-1.0
-0.5
0.0
0.5
1.0
1.5
2.0
2.5
3.0
0.0
0.5
1.0
1.5
2.0
2.5
3.0
u
Probability density (u)
0.5
0.4
0.3
0.2
0.1
p
up
0.0
-3.0
-2.5
-2.0
-1.5
-1.0
-0.5
u
Figure 4.1. Definition of the fractile for the standardised random variable U.
The values up of the lower fractile of a standardized random variable U having normal
distribution for selected probabilities p are given in Table 4.1. Considering the symmetry of
the normal distribution, the values up of the upper fractile can be assessed from Table 4.1 by
substituting of p by 1 - p and by changing the sign of values up (from negative to positive).
A - 18
0,500
-5,199 -4,753 -4,265 -3,719 -3,091 -2,327 -1,645 -1,282 -0,841 0,000
0,3
0,2
0,1
1- p = 0,05
p = 0,05
u0,05 = -1,645
0,0
-3,5
-2,5
-1,5
u0,95 = 1,645
-0,5
0,5
2,5
1,5
3,5
1+V
(4.3)
where unorm,p is the fractile of a standardized random variable with normal distribution, is
the mean and V the coefficient of variation of the variable X. An approximation of relation
(4.3) is often applied in the form
xp exp (unorm,p V)
(4.4)
whose accuracy is fully satisfactory for the coefficient of variations V < 0,2, but it is
commonly used for greater V as well.
A - 19
0,50
0,80
0,90
0,95
0,24
0,77
0,97
1,89
1,28
1,42
1,49
0,20
0,81
1,04
1,21
1,45
1,65
1,77
0,15
0,84
1,13
1,34
1,68
1,99
2,19
0,08
0,85
1,21
1,49
1,98
2,46
2,81
0,00
0,84
1,28
1,65
2,33
3,09
3,72
0,80
1,32
1,77
2,70
3,86
4,94
0.74
1,32
1,85
3,03
4,70
6,40
0,68
1,29
1,89
3,31
5,51
7,97
0,61
1,24
1,89
3,52
6,24
9,52
-4
10
-3
10
0,01
0,05
0,10
Example 4.1.
Let us assess the fractile xp of a normal and lognormal distribution with lower limit at
zero for p = 0,001; 0,01; 0,05 and 0,10, if V = 0,3. We know that the lognormal distribution
with lower limit at zero has, in this case, a positive skewness = 0,927 (according to equation
(3.11)), which needs to be known for interpolation in Table 4.2. The resultant values xp are
given in the following table in the form of dimensionless coefficients xp/ (expressing the
ratio of the fractile to the mean), which were for normal and for lognormal distribution
assessed by different ways.
Table of coefficients xp/.
Probability p
Coefficient xp/ for
0,001
0,073
0,385
0,483
0,591
0,658
0,387
0,484
0,591
0,657
0,396
0,496
0,610
0,681
Table of coefficients xp/ shows the expected difference between the fractiles of
normal and of lognormal distributions. The lower fractile of normal distribution is
significantly lower than the corresponding fractile of lognormal distribution particularly for
small probabilities p. The table also shows that the approximate formula (4.4) provides
satisfactory results for computation of fractile of lognormal distribution (the error will
decrease with decreasing coefficient of variation V).
The fractile of gamma distribution can be calculated from the available tables for type
III Pearson distribution [12, 13]. To calculate the fractile of beta distribution, the available
tables of incomplete beta function may be used or it can be assessed by integration of
probability density function according to definition (4.1). However, when it is needed (and
A - 20
(4.5)
where mode xmod and parameter c are substituted by relations (3.25) and (3.26).
Example 4.2.
Let us determine the upper fractile of wind pressure from Example 3.4 described by
Gumbel distribution when probability p = 0,98 is considered. It is known from Example 3.5
that for the one-year maximum 1 = 0,35 kN/m2, 1 = 0,06 kN/m2. The fractile x0,98 for such
parameters follows from equation (4.5)
x0,98 = 0,35 (0,45 + 0,78 ln(-ln(0,98))) 0,06 = 0,51 kN/m2
The corresponding fractile of the maximum for a period of 50 years (as shown in
Example 3.4 that 50 = 0,53 kN/m2, 50 = 0,06 kN/m2) is
x0,98 = 0,53 (0,45 + 0,78 ln(-ln(0,98))) 0,06 = 0,69 kN/m2
Simple mathematical procedures, including the computation of fractile, are some of
the reasons of the wide popularity of Gumbel, distribution which is frequently used for
random variables describing climatic and other variable actions that are defined by maximal
values for a given period (e.g. during one year).
However, theoretical models are not always known in practical applications. In civil
engineering, the fractile of a random variable (e.g. strength of a new or unknown material) has
to be assessed from a limited sample, the size n of which may be very small (sometimes less
than 10). Furthermore, considered random variables may have a high variability (the
coefficient of variation is sometimes greater than 0,30). Assessment of the fractile of a
population from a sample is then a serious problem, which is in mathematical statistics solved
by various methods of estimation theory. In the following three basic methods are shortly
described: the coverage method, the prediction method and the Bayesian method for
estimation of the population fractile.
4.2
(4.6)
Thus, the estimator xp,cover is lower (on the safe side of the lower fractile) than the
unknown fractile xp with the probability (confidence) .
A - 21
(4.7)
If the standard deviation of the population is unknown, then the sample standard
deviation s is considered
xp,cover = m kp s
(4.8)
(4.9)
It can be shown that for growing n the estimator xp,pred defined in this way is
asymptotically approaching the unknown fractile xp. It can be also shown that the estimator
xp,pred corresponds approximately to the estimator obtained by the coverage method xp,cover for
confidence = 0,75 [4].
If the standard deviation of the population is known, then the lower p-fractile is
estimated by the value xp,pred according to the relation
xp,pred = m + up (1/n + 1)1/2
(4.10)
(4.11)
4.4
kp for = 0,75
1,64
tp(1/n+1)1/2
0
0
10
15
20
Figure 4.3. Coefficients kp and -tp(1/n + 1)1/2 for p = 0,05 and normal distribution of the
population (when is unknown).
Also the skewness (asymmetry) of the population may affect significantly the
estimator of the populations fractile. Tables 4.5 and 4.6 show the coefficients kp from
equation (4.8) for three value of the skewness = 1,0, 0,0 and 1,0, probability p = 0,05 and
confidence = 0,75 (Table 4.5) and = 0,95 (Table 4.6). Values of the coefficients from
Table 4.6 are shown in Figure 4.4.
Table 4.5. Coefficient kp from equation (4.8) for p = 0,05, = 0,75 and lognormal distribution
having skewness (when is not known).
Sample size n
Skewness
3
4
5
6
8
10
20
30
= 1,00
= 0,00
= +1,00
1,64
n
0
0
10
15
20
Figure 4.4. Coefficient kp for p = 0,05 and confidence = 0,95 (when is unknown).
A similar dependence on skewness may be observed in the case of the generalized
Students t-distribution for which the fractiles tp are given in Table 4.7. These values tp are
applied in the prediction method using formula (4.11) and further in the Bayes method. That
is why Table 4.7 gives directly the values of fractiles tp depending on the number of degrees
of freedom . Similarly as in Tables 4.5 and 4.6 the probability p = 0,05 and three skewnesses
= 1,0; 0,0 and 1,0 are considered.
Table 4.7. Coefficient tp from equation (4.11) for p = 0,05 and lognormal distribution with
skewness (when is unknown).
Coefficient tp for = n 1 degrees of freedom
Skewness
3
4
5
6
8
10
20
30
A - 25
1
+ 1 4,6 = 20,4 MPa
5
where the value tp = -1,74 is given in Table 4.7 for = 1,0 and = 5 - 1= 4. The resulting
strength is in this case by 10% greater than the value, which corresponds to the normal
distribution ( = 0).
4.5
A - 26
(4.12)
(4.13)
Both the unknown variables n' and ' may be assessed independently (generally ' n'
- 1), depending on previous experience with the degree of uncertainty of estimator of the
mean and standard deviation of the population.
The next step of the procedure applies the prediction method of fractile estimation.
The Bayes estimator xp,Bayes of the fractile is given by relationship similar to equation (4.11)
for prediction estimator, assuming that the standard deviation of the population is not
known
1/ 2
(4.14)
x p,Bayes = m' ' +t p (1 / n' ' +1) s' '
where t' p' = t'p' ( , p, ' ' ) is a fractile of the generalised Students t-distribution having an
appropriate skewness for '' degrees of freedom (that is generally different from the value
n'' 1).
If the Bayes method is applied for an assessment of material strength, the advantage
may be taken of the fact that the long-term variability is constant. Then the uncertainty of an
assessment of and the value v() are relatively small, variables ' assessed according to
equation (4.13) and '' assessed according to equation (4.12) are relatively high. This factor
may lead to a favourable decrease of the value t'p' and to augmentation of the estimator of the
lower fractile of xp according to equation (4.14). On the other hand, uncertainties in
assessment of the mean and the variable v() are usually great and previous information
may not affect significantly the resulting values n'' and m''.
If no previous information is available, then n' = ' = 0 and the resulting characteristics
m'', n'', s'', '' equal the sample characteristics m, n, s, . In this case the Bayes method is
reduced to the prediction method and equation (4.14) becomes equation (4.11); if is known
equation (4.10) is used. This particular form of the Bayes method, when no previous
information is available, is considered in Eurocode EN 1990 [1] and international standards
ISO [2, 3].
Example 4.4.
If previous experience was available for Example 4.3, the Bayes method could be
used. Suppose that the information is m = 30,1 MPa, V() = 0,50, s = 4,4 MPa, V() = 0,28.
It follows from equation (4.13) that
2
1
4,4 1
n' =
6
< 1 , ' =
2 0,282
30,1 0,50
Further on these values are thus considered: n' = 0 and ' = 6. Because = n 1 = 4, it
follows from equation (4.12)
n'' = 5, '' = 10, m'' = 29,2 MPa, s'' = 4,5 MPa.
From equation (4.14) the fractile estimate follows as
x p , Bayes = 29,2 1,81
1
+ 1 4,5 = 20,3 MPa
5
A - 27
1
+ 1 4,5 = 21,9 MPa
5
(4.15)
For design values xd of dominating variables it holds approximately that p = 0,001 (or
another value close to this one), i.e. it holds
P(X < xd) = 0,001
(4.16)
(4.17)
A - 28
1/2
up(1/n+1) , known 2,31 2,01 1,89 1,83 1,80 1,77 1,74 1,72 1,68 1,67 1,64
- 3,37 2,63 2,33 2,18 2,00 1,92 1,76 1,73 1,64
tp(1/n+1)1/2, unknown Table 4.9. Coefficients kn for a design value xd of a dominating variable, P(X < xd) = 0,001.
Sample size n
Coefficient
1
2
3
4
5
6
8
10 20 30
1/2
4,36
3,77
3,56
3,44
3,37
3,33
3,27
3,23
3,16
3,13
3,09
up(1/n+1) , known
- 11,4 7,85 6,36 5,07 4,51 3,64 3,44 3,09
tp(1/n+1)1/2, unknown Table 4.10. Coefficients kn for a design value xd of a non-dominating variable, P(X < xd) = 0,1.
Sample size n
Coefficient
1
2
3
4
5
6
8
10 20 30
1/2
up(1/n+1) , known 1,81 1,57 1,48 1,43 1,40 1,38 1,36 1,34 1,31 1,30 1,28
tp(1/n+1)1/2, unknown - 3,77 2,18 1,83 1,68 1,56 1,51 1,45 1,36 1,33 1,28
The assumption concerning knowledge of the standard deviation is replaced
(inaccurately) in the document by the assumption that the coefficient of variation V is known.
The original version of Table 4.9 [1] gives for the sample size of a wrong value of 3,04 for
the coefficients (correct is 3,09). Let us also note that when knowledge of the standard
deviation is assumed, Tables 4.8 to 4.10 give values of coefficients already for the sample
size n = 1. Application of these values is, however, associated with significant statistical
uncertainties and therefore a minimum sample size n = 3 is recommended here. Note, that
Table 4.10 (for 0,1 fractile) is included only in the prestandard ENV 1991-1 and not in the
final document EN 1990 [1].
Statistical methods for determining the characteristic and design values of resistance
variables are provided in Annex D Design assisted by testing of EN 1990 [1]. Relevant
basic variables describing structural resistance are described by lognormal distribution. The
whole procedure is described in detail in the Annex D. Attached MATHCAD sheet
Mod_est.mcd can be used to evaluate experimental data using the whole procedure. It is
provided with explanatory notes and needs no additional information.
In order to simplify computational procedure the assessment coefficients given in
Tables 4.8 and 4.9 are in the attached MATHCAD sheet Mod_est.mcd expressed using
built-in distribution function of normal and Student t- distribution. In accordance with the
principles of Annex D in [1] single variable and model representation of a resistance variable
R are distinguished. The results shown in the attached sheet indicates that both approaches
lead to similar results.
A - 29
A - 30
Domain
of X
Parameters
Mean
Standard
deviation
Skewness
1/(b a)
axb
a
b>a
(a + b)/2
(b a)/12
1 x 2
exp
2
2
x0 + c
3c+c3
3V+V3
k/
k/
2 /k
| x x 0 || c | 1 + c 2
exp ln
|x x 0| ln(1 + c 2 ) 2
x 1+V 2
exp ln
x ln(1 + V 2 ) 2
(2 ln(1 + c 2 ))
(2 ln(1 + V 2 ))
d 1
( x a ) (b x)
B(c, d )(b a ) c + d 1
c 1
x0 x < x0 = c
pro > 0,
c
<xx0
pro < 0
0x<
V = /
0x<
= /2
k = ( /)2
axb
a
b >a
c 1
d 1
0xb
d 1
( x) (b x)
B(c, d ) b c + d 1
A - 31
b >0
c 1
d 1
x < xmod =
0,5776 /
c = /(6)
2 g (d c)
(b a )
a+
,
,
cg
+
dg
c
+
d
+
2
(b a ) c
+
c + d +1 g = c + d +1
c+d
g=
cd
cd
b c
c+d
b
,
cg + dg
g=
xmod +
0,577/c
2 g (d c)
,
c+d +2
c + d +1 g = c + d +1
cd
cd
/(6c)
1,14
Function Z
aX+b
X2
1
X
*)
*)
aX + bY + c
The mean Z
Standard deviation Z
Skewness Z
aX + b
| a | X
X2 + X2
1 + V X2 V X3 X
X
aX + bY + c
2 X X2 + X X X
(V
2
X
2V X3 X
(a
2
X
2
X
+ b 2 Y2
X+Y
X + Y
2
X
+ Y2
XY
X Y
2
X
+ Y2
XY
X
Y
*)
X Y
*)
X (1 + VY2 VY3 Y )
Y
1/ 2
Z3
6 V X4 V X3 X
a 3 X3 X + b 3 Y3 Y
1/ 2
X3 Z3
1/ 2
Z3
1/ 2
X3 X + Y3 Y
Z3
1/ 2
X3 X Y3 Y
Z3
X Y (V X2 + VY2 + V X2 VY2 )1 / 2
X (V X2 + VY2 2 VY3 Y )
Y
1/ 2
A - 32
8 X3 X3 ( X + 3 V X )
X3 Y3 (V X3 X + VY3 Y + 6V X2 VY2 )
Z3
X3 (V X3 X VY3 Y + 6VY4 + 6V X2 VY2 )
Y3 Z3
Domain of
X
Rectangular
R(a,b)
Normal
N(,)
axb
a + p (b a)
+ up = (1+ up V)
m p
p from
Table 4.3
Lognormal,
general
LN(,,)
LN(,,x0)
x0 x <
1
pro > 0, 1
exp sign ( )u p ln(1 + c 2 )
2
c
1+ c
< x x0
pro < 0
+ x0
= x0 +
exp sign ( )u p ln(1 + c 2 )
2
1+ c
known
unknown
m kp s
kp from
Table 4.4
m+up(1/n+1)1/2
up from Table 4.1
m+tp(1/n+1)1/2 s
tp(1/n+1)1/2 from
Table 4.4
m kp s
m+up(1/n+1)1/2
m+tp(1/n+1)1/2 s
kp from 4.5
and 4.6
m p
m kp s
m+up(1/n+1)1/2
m+tp(1/n+1)1/2 s
p not given
kp from
Table 4.5
and 4.6
m p
p not given
+ up = (1+ up V)
Lognormal,
zero origin
LN(,)
0x<
1+V 2
exp u p ln(1 + V )
2
+ up = (1+ up V)
and up for lognormal distribution from Table 4.2
Gumbel
Gum(,)
x <
1
x mod ln( ln( p ))
c
(0,45 + 0,78 ln( ln( p )))
A - 33
ATTACHEMENTS
1. MATHCAD sheet DistFract.mcd Fractiles of basic types of distributions.
MATHCAD sheet DistFract.mcd is intended for determination of fractiles of selected
theoretical models.
2. MATHCAD sheet SampFract.mcd Estimation of sample fractile.
MATHCAD sheet SampFract.mcd is intended for determination of fractiles using
limited samples.
3. MATHCAD sheet Mod_est.mcd Estimation of models.
MATHCAD sheet Mod_est.mcd is intended for determination of fractiles using limited
samples taking into account model uncertainties.
A - 34
:= 1
V := 0.10
:= V
check:
:= 3.8 P := pnorm ( , 0 , 1)
P = 1.183 10
u ( p) := qnorm( p , 0 , 1)
ln ( p) :=
exp u ( p) ln 1 + V
1+ V
lna ( p) := exp( u ( p) V )
a := 1 , 0.5 .. 1
1
2 + 4 + a a2 + 4 a
a
C ( a) :=
1
ln 1 + C ( a)
x0 ( a) :=
1
C ( a)
) + ln ( ) ( 0.5) ln (1 + C ( a) 2)
Check:
2
exp sign ( a) u ( p) ln 1 + C ( a)
lng ( p , a) := 1
1
C ( a)
2
1 + C ( a)
A - 35
x0 ( 1) = 0.69
k :=
:=
gam ( p) :=
qgamma ( p , k )
Lower fractiles
V = 0.1
Check:
n ( 0.001) = 0.691
ln ( 0.001) = 0.731
lna ( 0.001) = 0.734
lng ( 0.001 , 1) = 0.801
lng ( 0.001 , 1) = 0.53
lng ( 0.001, 0.4) = 0.743
gam ( 0.001) = 0.719
n ( p )
0.9
ln ( p)
lng ( p , 1)
0.8
lng ( p , 1)
gam ( p )
gum ( p )
0.7
0.6
0.01
0.1
Upper fractiles
1.8
n ( p )
ln ( p)
1.6
lng ( p , 1)
lng ( p , 1)
1.4
gam ( p )
gum ( p )
1.2
0.9
0.92
0.94
0.96
p
A - 36
0.98
ks ( n ) := qt ( 0.95, n 1) 1 +
5% fractile V known
1
n
1
k( n ) := qnorm( 0.95, 0 , 1) 1 +
ds ( n ) := qt ( 0.999, n 1) 1 +
d( n ) := qnorm( 0.999, 0 , 1) 1 +
1
n
exp( ks ( n ) ) ln 1 + V
2
ds ( n , V) :=
1+ V
Single variable, V known: xk= k *x, xd= d*x
k( n , V) :=
exp( k( n ) ) ln 1 + V
2
d( n , V) :=
exp( ds ( n ) ) ln 1 + V
1+ V
1+ V
exp( d( n ) ) ln 1 + V
2
1+ V
rt := 0.707 := 0.707
2
( 2
)
2
2
dmod( n , V) := exp( 3.09rt + d( n ) + 0.5 V) V
3. Partial factor M
Check values:
Ms ( n , V) :=
M( n , V) :=
Model, V unknown
Model, V known
ks ( n , V)
ds ( n , V)
k( n , V)
d( n , V)
ksmod ( n , V)
Msmod ( n , V) :=
dsmod ( n , V)
kmod( n , V)
Mmod( n , V) :=
dmod( n , V)
A - 37
V := 0.01, 0.02.. .4
0.9
ks( 10 , V)
0.8
ds( 10 , V)
ksmod( 10 , V)
dsmod( 10 , V)
0.7
0.6
0.5
0.4
0.1
0.2
0.3
0.4
Figure 1. Characteristic and design values versus coefficient of variation V for n = 100
and weighting factors rt= 0,301, = 0,955
5. Graphs for M factors
1.8
Ms ( 10 , V)
M ( 10 , V)
1.6
Msmod ( 10 , V)
M mod( 10 , V)
1.4
1.2
0.1
0.2
0.3
0.4
Figure 2. Partial factor M versus coefficient of variation V for n=100 and weighting factors
rt= 0,301, = 0,955
A - 38
1. Experimental data
rt := ( 1 2 3 4 5 6 7 8 9 10 )
re := ( 1 3 4 4 5 5 6 8 10 9 )
Reading experimental data from the file "rdata.prn" located in the same directory
Check values:
1
re := DATA
nr := length ( rt)
nr = 21
The least square fit for y=a+bx b := slope ( rt , re) a := intercept ( rt , re)
The least square fit for y=bx, a = 0
b = 0.975
a = 0.032
b = 1.014
rt rt
1
re
b rt
mrt = 0.705
1.5
2. Check of experimental
data taken from the file
"rdata.prn"
and check of the least
square fit - slope b
0.5
:= ln( )
re
:=
b rt
re rt
b :=
mre = 0.719
0.5
1.5
rt
m := mean( )
Characteristics of :
s2 :=
nr var ( )
nr 1
Check values:
V := exp( s2) 1
V = 0.107
Vrt = 0.094
2
Vr := V + Vrt
Vr = 0.143
Standard deviations
Qrt := ln Vrt + 1
Weighting factors:
Q := ln V + 1
rt :=
Qrt
Q
:=
Q := ln Vr + 1
Q
Q = 0.142
rt = 0.663
A - 39
= 0.752
ks ( n ) := qt ( 0.95, n 1) 1 +
Check values:
k( n ) := qnorm( 0.95, 0 , 1) 1 +
n := 3 , 3.5.. 30
ks ( 10) = 1.923
1
n2
k( 10) = 1.745
1
ds ( 10) = 4.507
n
1
d( n ) := qnorm( 0.999, 0 , 1) 1 +
d( 10) = 3.278
n2
4. Characteristic and design values (relative values related to the mean
V := 0.0, 0.001.. .4
Note that the range variable V is generally used for the coefficient of variation of a single variable
V and for a model investigation Vr, these may be different, for example V = 0,12 and Vr = 0,142.
Single variable, V unknown: rk=ks*mr, rd=dsmr
ks ( n , V) :=
exp( ks ( n ) ) ln 1 + V
ds ( n , V) :=
1+ V
1+ V
exp( ds ( n ) ) ln 1 + V
exp( k( n ) ) ln 1 + V
d( n , V) :=
1+ V
exp( d( n) ) ln 1 + V
2
1+ V
Model, V unknown: rk= ks*rx, xd=dsmr, weighting factors rt and taken from the
above experimental data re and rt (for V < 0,4 approximately Q =~ V):
Check values:
ksmod ( nr , Vr) = 0.774
dsmod ( nr , Vr) = 0.608
Model, V known: rk=k *mx, rd= dmr, (for V < 0,4 approximately Q =~ V):
2
( 2
)
2
2
dmod( n , V) := exp( 3.09rt + d( n ) + 0.5 V) V
ks ( n , V)
ds ( n , V)
k( n , V)
M( n , V) :=
Model, V unknown
Msmod ( n , V) :=
Model, V known
Mmod( n , V) :=
d( n , V)
ksmod ( n , V)
dsmod ( n , V)
kmod( n , V)
dmod( n , V)
A - 40
6. The relative characteristic and design values, the model values estimated using
weighting factors determined above from experimental data given in the file "rdata.prn"
1
0.9
ks ( 10 , V)
ds ( 10 , V)
0.8
k( 10 , V)
d( 10 , V)
0.7
ksmod ( 10 , V)
dsmod ( 10 , V)
0.6
kmod( 10 , V)
dmod( 10 , V)
0.5
0.4
0.1
0.2
0.3
0.4
Figure 2. Variation of the characteristic and design values with coefficient of variation V for n = 10.
The char. value of X:
Example
rks( n , V) := b mrt ks ( n , V)
The model char. value of X: rksmod( n , V) := b mrt ksmod ( n , V) From data: rksmod( nr , Vr) = 0.553
7. Partial factor M, the model values estimated using weighting factors determined
above from experimental data given in the file "rdata.prn"
2
1.8
Ms ( 10 , V)
M ( 10 , V)
1.6
Msmod ( 10 , V)
M mod( 10 , V)
1.4
1.2
0.1
0.2
V
0.3
0.4
Figure 3. Variation of the partial factor M versus coefficient of variation V for n=10.
A - 41
ks ( 3 , 0.0) = 1
ks
Figure 4. Variation of the characteristic values ks with n and V
.
Ms ( 20, 0.1) = 1.208
Ms
Figure 5. Variation of the partial factor M with n and V .
A - 42
Summary
Using basic principles of the reliability theory described in Chapter II Elementary
methods of structural reliability I, the operational techniques for estimating partial factors of
basic variables are derived and applied to common permanent and variable loads. The
described computational procedures are illustrated by a number of numerical examples, which
are supplemented by MATHCAD, EXCEL and MATHEMATICA sheets.
1
INTRODUCTION
1.1
Background materials
Fundamental concepts and procedures of structural reliability are well described in a
number of national standards, in the new European document EN 1990 [1] and International
Standard ISO 2394 [2]. Additional information may be found in the background document
developed by JCSS [3] and in recently published handbook to EN 1990 [4]. Guidance on
application of the probabilistic methods of structural reliability may be found in publications
and working materials developed by JCSS [5] and in relevant literature listed in [4] and [5].
Guidance on structural systems and time dependent reliability may be found in [6] and [7].
1.2
General principles
The theory of structural reliability considers all basic variables as random quantities
having appropriate types of probability distribution. Different types of distributions should be
considered for actions, material properties and geometrical data. In addition, model
uncertainties of actions and resistance models should be taken into account. Prior theoretical
models of basic variables and procedures for probabilistic analysis are indicated in JCSS
documents [5].
This appendix is a direct extension of Chapter II "Elementary methods of structural
reliability" of the main text to which reference is frequently made.
2
DESIGN POINT
B-1
E / E
R / R
R
(1)
Ed = E EE
(2)
Here E and R denote the so-called sensitivity factors of variables E and R. The minus
signs in equations (1) and (2) are conventionally accepted in the documents CEN [1] and ISO
[2].
B-2
E = E / E + R
(3)
R = R/ E + R
(4)
(5)
2
2
E = E / E + R = 0,7
(6)
The validity of such an approximation is delimited by means of a condition for the ratio of the
standard deviations in the form
0,16 < E / R < 7,6
(7)
When this condition is not satisfied, then the sensitivity factor = 1,0 is
recommended to be used for the variable having a greater standard deviation. Let us remark
that this simplification is on the safe side as the sum of squared direction cosines should be
equal to one.
The design values Ed and Rd of the variables E and R are thus defined as fractiles of
normal distribution
P(E > Ed) = U(+E) = U(-0,7)
(8)
(9)
where U(u) denotes a standardized normal distribution. If = 3,8, then the design values ed
and rd are fractiles corresponding approximately to probabilities 0,999 and 0,001. Note that in
equation (9) the use of the symmetry of normal distribution is taken into account, i.e. of the
relationship 1 - U(+E) = u(-E).
When the load or resistance model contains several basic variables (other loads, more
materials, geometrical data), equations (8) and (9) holds only for the leading variables (the
most significant for the studied condition of reliability). For other (accompanying) variables
the requirements on design values are reduced and it holds
P(E > Ed) = U(+0,4E) = U(-0,28)
(10)
(11)
When = 3,8 the design values of accompanying (non-leading) variables are fractiles
approximately corresponding to probabilities 0,9 and 0,1.
Design values are the upper fractiles (for actions) or the lower fractiles (for
resistance), corresponding to certain probabilities of being exceeded (actions) or not reached
(resistance). For leading variables, the probabilities are given by the distribution function of
the normal standardized distribution for values u = +E and -R, in the case of non-leading
variables for reduced values u = +0,4E and 0,4R. These probabilities (for the lower
fractile approximately 0,001 for leading and 0,1 for accompanying variables) then serve to
determine the design values even for those variables, which do not have normal distribution.
Let us note that according to the general principles, it is necessary in the case of upper
fractiles (actions) to consider the complementary probabilities (close to the value 1).
B-3
Example 1.
The design values Ed and Rd of variables E and R from Example 4 will be assessed
assuming that the reliability index = 3,8, E = -0,7 and R = 0,8. According to equation (8),
it holds for E that
P(E > ed) = U(E) = U(-2,66) = 0,0039
The complementary probability is therefore 0,9961 and we obtain from equation
ed = - (0,45 + 0,78ln(-ln(p))) = 50 (0,45 + 0,78ln(-ln(0,9961)))10 = 88,75
Note that when the normal distribution is assumed, it is given
ed = + up = 50 + 2,66 10 = 76,6
According to equation (9), it holds for R
P(R < Rd) = U(-R) = U(-3,04) = 0,0012
For the lognormal distribution with the mean of 100 (units) and standard deviation of
10 (units) it follows
Rd exp(unorm,p V) = 100 exp(-3,040,10) = 73,79
For normal distribution we obtain
Rp = + up = 50 3,04 10 = 69,6
Obviously, it holds for the coordinates of the design point that ed > rd and the tie rod
does not satisfy the condition (1) of Chapter II (we know from the Example 4 of Chapter II
that is only 3,09). In order to satisfy the condition for a reliability index of 3,8, the
parameters of variables E and R would have to be modified.
The attached MATHCAD sheets StRod.mcd, DesVRod.mcd may be used to make all
numerical calculations.
3
PARTIAL FACTORS
3.1
Material properties
The above described reliability concepts may be used to assess the partial factors. The
attached MATHCAD sheets GammaRGQ.mcd cover all computational procedures described
below and may be used to make additional numerical calculations.
In accordance to EN 1990 [1] or ISO 2394 [2] the partial factor R of material
resistance R is defined as the fraction of its characteristic value Rk and the design value Rd,
thus
R = Rk / Rd
(12)
B-4
(13)
(14)
In equation (13) and (14) R denotes the mean, R the standard deviation, VR the coefficient of
variation and R = 0,8 the sensitivity factor of R.
Taking into account equations (13) and (14) it follows from (12) that the partial factor
R for a normal distribution of R can be assessed as
(15)
(16)
(17)
Taking into account equations (16) and (17) it follows from (12) that the partial factor R for a
lognormal distribution can be assessed as
(18)
Figures 2 and 3 show the variation of the partial factor R of the material property R
with the reliability index for selected values of the coefficient of variation VR given for
normal distribution by equation (15) (Figure 2) and lognormal distribution by equation (18)
(Figure 3).
2,5
2,0
VR = 0,20
0,15
0,10
0,05
1,5
1,0
0,5
0
Figure 2. Variation of R with for selected coefficients of variation VR = 0,05; 0,10; 0,15
and 0,20, and for normal distribution of R.
Generally the partial factor R increases with increasing . The increase of R is
considerably greater in the case of normal distribution (Figure 2) than in the case of lognormal
distribution (Figure 3). The effect of the type of distribution is particularly obvious for
B-5
R
2
VR = 0,20
0,15
0,10
0,05
1,5
0,5
0
Figure 3. Variation of R with for selected coefficients of variation VR = 0,05; 0,10; 0,15 and
0,20, and for lognormal distribution of R.
3.2
Permanent load
Consider a self-weight G having a normal distribution. Similarly as in the case of
material property when a reference period T instead of the design working life Td is used in
the reliability verification of a structure, then the design value of G should be determined for
T instead of Td. It is assumed that the characteristic value Gk of G is defined as the mean G
[1], [2] and [5]:
Gk = G
(19)
The design value Gd is given by equation (23) of Chapter II (see also documents [1], [2]) as
Gd = G G G = G + 0,7 G = G(1 + 0,7 VG)
(20)
In equation (20) G denotes the mean, G the standard deviation, VG the coefficient of
variation and G = 0,7 the sensitivity factor of G.
The partial factor G of G is given as [1], [2]
G = Gd / Gk
(21)
Taking into account equations (19) and (20) it follows from (21) that
G = (1 +0,7 VG)
(22)
Figure 4 shows the variation of the partial factor G with the reliability index for
selected values of the coefficient of variation VG = 0,05; 0,10; 0,15 and 0,20. Note that G =
1,35 (recommended in EN 1990 [1]) corresponds approximately to the reliability index =
3,8 if the coefficient of variation is about 0,1 (the value recommended in EN 1990 [1] was
further increased by 5% to take into account model uncertainty).
B-6
2,5
G
2
VG = 0,20
0,15
0,10
0,05
1,5
0,5
0
Figure 4. Variation of G with for selected coefficients of variation VG = 0,05; 0,10; 0,15 and
0,20, and for normal distribution of G.
It follows from Figures 2, 3 and 4 that less significant variation with -values should
be generally expected for the partial factor of self-weight G than for the partial factor of
material property R..
3.3
Variable load
A similar procedure as in the case of permanent load G can be used for estimation of
the partial factors Q for variable loads Q. Assuming Gumbel distribution the characteristic
value (0,98 fractile) is given as
Qk = Q (1 VQ (0,45 + 0,78 ln(ln(0,98))))
(23)
(24)
In equation (23) and (24) G denotes the mean, VQ the coefficient of variation of annual
extremes of Q and G = 0,7 the sensitivity factor of Q.
The partial factor Q of Q is given as [1], [2]
Q = Qd / Qk
(25)
Figure 5 shows the variation of Q with the coefficients of variation VQ for selected
values of assuming Gumbel distribution of Q. It appears that in case of a variable action Q
the reliability index has a significant effect on the partial factor Q. The following Figure 6
shows the variation of Q with the reliability index for selected values of the coefficients of
variation VQ assuming again Gumbel distribution of Q.
It follows from Figures 5 and 6 that for the reliability index = 3,8 and the coefficient
of variation VQ up to 0,5, the partial factor Q is less than 1,3. However the coefficient of
variation may be also greater than 0,5 and other distribution may be more adequate (see other
Chapters of this Handbook). That is why a conservative value Q = 1,5 is recommended in
EN 1990 [1].
B-7
1.5
= 4,3
1.4
1.3
= 3,8
1.2
= 3,3
1.1
0.9
0.1
0.2
0.3
0.4
0.5
VQ
Figure 5. Variation of Q with the coefficients of variation VQ for selected values of
assuming Gumbel distribution of Q.
2
VQ = 0,4
VQ = 0,3
1.5
VQ = 0,2
1
0.5
2.5
3.5
4.5
Figure 6. Variation of Q with the reliability index for selected values of the coefficients of
variation VQ assuming Gumbel distribution of Q.
B-8
4.1
General
In the Chapter II Elementary methods of structural reliability I of this handbook it
has been presented what is called the fundamental case of structural reliability. It is the case in
which the limit state function can be represented by only two random independent variables,
the effect of the actions and the resistance. This fundamental case represents a very interesting
case for introducing the reliability concepts due to its very intuitive reasoning, and that we are
used to this concepts of global action effect and resistance. Also, due to the fact that in two
dimensions only is easy render graphic simple representations.
But, unfortunately, only in a few cases the limit state function may be represented by
the fundamental case, reducing the structural reliability problem to a simple resistance versus
action effect formulation assuming their independents. In most cases, at least, a few more
variables will be needed. In general the resistance is a function of the material(s) properties
and the dimensions of the structure or element and the action effects depend on the various
applied loads and densities and dimensions of the structure. Even, not always the resistance
and action effects can be considered independent: for instance, some dimensions affect both
the actions and resistances. Even more, the action effects may depend on the response of the
structure as a whole (e.g., the dynamic actions).
The limit state function can be written as:
Z(X1, X2, X3, ...) = 0 ;
(26)
Z(X) = 0 ;
(27)
or in vector form:
where X = {X1, X2, X3, ...} is the vector of the random variables depending on time defining
the limit state function.
In this case also, Z(X) > 0 represents the safe region and Z(X) < 0 the unsafe region.
The probability of failure, in this case is obtained as:
Pf = P [Z(X) 0] =
...
fZ(x) dx;
(28)
Z ( X )< 0
where fZ(x) is the joint density probability distribution of the vector of variables X.
In the case of all the basic variables are independent, in a lot of cases we can adopt this
hypothesis as a good approximation, the joint density probability distribution of X is the
product of the marginal density probability distribution of each variable. Therefore, is
possible to put the equation (18) in the form:
Pf = P [Z(X) < 0] =
...
Z ( X )< 0
B-9
(29)
Basic variables
Basic variables are those variables needed to characterize and define the behaviour of
a structure, and its safety, relating to a determined limit state.
The designer has some degree of freedom in order to choose the basic variables: In
general, those used in the conventional design are chosen: dimensions, weights, loads,
material strengths, etc. In general, independent basic variables are considered as any
dependence adds complexity and it is difficult to define the degree of dependence. However,
some variables as the tension, compression strength and the modulus of elasticity of the
material are known to be dependent, but one can, generally, deal with them as they were
independent.
In order to analyse structural reliability it is necessary to characterize basic variables
statistically. That is, to define, at least, their distribution functions, parameters and correlation
matrix. The parameters of distributions can be estimated on the basis of data by means of the
usual statistical techniques: maximum likelihood, methods of moments, etc. The data should
be carefully investigated in order to exclude the outliers, to analyse trends, and so on. A
graphical representation of data and of the model adopted is generally useful.
Guidance concerning distribution functions and their parameters for the models of
common actions and resistances in structural reliability are given in reference [5].
Example 1:
Consider a basic variable and the following experimental data:
data={1.3,3.2,4.3,1.3,5.4,3.7,3.8,4.0,2.9,3.2,4.5,4.0,3.4,2.4,1.8,1.7,2.2,4.1,2.6,4.1,3.3,
3.5,3.7,2.4,2.8,2.5,3.,3.3,2.6,2.9,2.4,2.6,2.9,2.8,3.1,3.2,3.4,3.5,3.2,3.3,3.};
In the following the results obtained with the attached MATHEMATICA package
distribution-fit are used. The following statistical characteristics were obtained:
Number of data
41
Mean
3,105
Variance
0,700
Standard variation
0,837
Coefficient of variation
0,270
We try to fit the normal, lognormal and gamma distributions to data. The parameters
of the distributions are obtained by the method of moments. Assuming the following
intervals: {1.3-2.3, 2.3-2.7, 2.7-3.3, -3.7, 3.7-4.1, 4.1-6}, the Chi square test is performed.
In the figure are represented the histogram and the density distribution function (in
this case the normal distribution). The results obtained for each distribution are shown in the
following table:
normal
0.5
0.4
Distribution Parameters
Estimator
W
Normal
3,103- 0,847 0,612
Gamma
13,40 - 0,231 1,001
Lognormal 1,096 0,268 1,310
Confidence
Level %
96,0 - 99,6
91,0 98,6
86,0 97,1
0.3
0.2
0.1
2
The range of the confidence level is due to the different hypothesis. The moments are
known a priori, then the degrees of freedom are the number of intervals minus one; or
B - 10
.999
0.8
.99
.95
.9
0.6
0.4
.5
.1
.01
.001
0.2
5.1
General
Various design concepts mentioned above may be illustrated considering a simple
example of a reinforced concrete slab in an office building. The example shows how different
design methods (permissible stresses, global safety factor, partial factor method) treat
B - 11
Design of a slab
Consider first a simple drawing of the slab cross-section including simplified stress
distribution diagrams in the compressive concrete zone (rectangular and triangular) as
indicated in Figure 8.
fc
fc
0,8bxfc
0,5bxfc
0,8x
x
d M
As
Asfy
(30)
As fy(d 0,4 x) = M
(31)
The basic variables used in equations (30) and (31) are evident from Figure 8:
d denotes the effective depth, x the depth of the neutral axis, b the width of the slab
(considered as 1 m), As the area of the reinforcement, fc the concrete strength and fy the
reinforcement strength (yield point). The bending moment in equations (31) is given as
B - 12
M=
( g + q ) L2
8
(32)
(33)
Without going into technical details note, that equation (33) may be used
approximately for the global safety factor methods and partial safety factor method. Attached
EXCEL sheet RCBeam and MATHCAD sheet RCBeam may be applied to make necessary
calculations.
The classical permissible stresses method assumes a triangular compressive stress
block in the compressive concrete zone, which is also indicated in Figure 8, and linear stress
strain diagram. The reinforcement area As can be found using the attached MATHCAD sheet
RCBeam.mcd, which includes further calculation details.
Note that up to now all the basic variables have been considered as deterministic
quantities without considering any uncertainties that may potentially affect their actual values.
However, it is well recognised that some of the basic variables entering equation (31) might
have considerable scatter (particularly the load effect (bending moment) M, concrete strength
fc, and reinforcement strength fy). On the other hand the geometric data As, b and d seem to be
much less uncertain (almost fixed or deterministic).
To get a first estimate of the area As, one may take the mean (average) values of all the
basic variables involved. Intuitively it is clear that this might be not safe enough and instead
of the mean values somehow conservative (safe) values should be applied. Table 1 indicates
the mean values together with values of basic variables as used in design in accordance with
the rules of above-mentioned design methods.
Table 1. Input design (characteristic) values of loads and material strengths used in design
calculation using different design methods.
Design method
Basic variable
The mean
Permissible Global
Partial
value
stresses
safety factor factor
Permanent load g [kN/m2]
7
7
7
9,45
2
Imposed load q [kN/m ]
0,8
3
3
4,5
Concrete compressive strength fc [MPa] 30
5,5
20
13,3
Reinforcement tensile strength fy [MPa] 560
275
500
435
Table 1 clearly indicates the differences between considered design methods. For
example, the input values of the permanent load g used in design calculation in accordance
with the permissible stresses method and the global safety factor method are equal to the
mean value (7 kN/m2), while the design value in the partial factor method (9,45 kN/m2) is
obtained as a product of the characteristic value and the partial factor G = 1,35. The design
value of imposed load q used in design calculation in accordance with the permissible stresses
method and the global safety factor method are equal to the characteristic value (3 kN/m2),
while the design value (4,5 kN/m2) in the partial factor method is product of the characteristic
value and partial factor Q = 1,5.
B - 13
As [m2]
0,00038
0,00204
0,00082
0,00069
MR [kNm]
35,1
228,9
97,4
82,4
0
8.0
5,0
4,2
pf
0,5
0,4410-16
0,3210-7
0,1210-5
The attached EXCEL sheet RCBeam and MATHCAD sheet RCBeam may be used to
check the resulting data indicated in Table 2. These sheets also indicate further details of
applied computational procedures.
It follows from Table 2 that resulting reinforcement areas vary within a broad range
from As = 0,00038 m2 (the mean value estimate) up to As = 0,00204 m2 (the permissible stress
method). The most economic procedure appears to be provided by the partial factor methods
that leads to the lowest reinforcement area As = 0,00069 m2.
5.4
Reliability consideration
The rectangular stress distribution (indicated in Figure 8), assumed in design
calculation according to the global or the partial factor method, is accepted for reliability
investigation of the slab designed by any of the design methods illustrated in Table 2
(including the permissible stress method, which assumes a triangular stress diagram). Thus,
taking into account equilibrium conditions (30) and (31) the limit state function (16) may be
written as
Z ( X ) = As f y (d
As f y
2bf c
) M = As f y (d
As f y
2bf c
( g + q ) L2
8
(34)
Assuming the limit state function (34) the reliability indexes and the failure
probabilities pf can be determined using commercially available software VaP [9] or
COMREL [10] or the MATHEMATICA notebook FORM.nb. An approximate value of the
reliability index and the failure probability pf (with an acceptable accuracy) may be
obtained using the attached MATHCAD sheet RelRCB, which is based on a simple procedure
of numerical integration. The self-content MATHCAD sheet RelRCB is provided with
B - 14
6.1
General
There are several different procedures for assessing failure probability in a general
case of more variables:
Analytical: Only in a few, very simple cases it is possible to find the analytical correct
solution. Depends on the variables vector, all must be independent and normally
distributed, and on the limit state region, it must be defined by hyper-planes. It cannot
be considered as a general solution.
Numerical: It is an exact solution in the sense that we can get, in principle, all the
precision we need. The simple trapezoidal rule of integration gives, in general good
results if there are not too many variables (4 or 5). The complexity of integration
increases exponentially with the number of variables.
Monte Carlo methods: The Monte Carlo simulation techniques are bases on the
random sampling of the variables and carry a long number of artificial experiments.
The use of this approach is increasing notably nowadays, due to bigger power and
speed of the computers. The crude application of the method lead to the same
difficulties explained above. Variance reduction and importance methods are
employed to avoid those difficulties.
First and Second Order Reliability Methods (FORM and SORM): These approximate
methods give iterative algorithms that allow to obtain the reliability index, using a
linear or quadratic approximation to the limit state surface in the point of minimum
distance to the mean point of the variables.
6.2
(35)
Obtain the length of the minimum distance vector from the new origin point to the
tangent hyper-plane of the limit state surface, referred to the new variables, in the
intersection point of that vector with the limit state surface.
6 Obtain the design point, (X1d, X2d,...,Xnd) and the sensitivity coefficients, 1,
2,...,n, the unitary cosines of that vector.
Figure 1 represents the space of the fundamental case, with only two variables: the
action effect and the resistance. The origin of the vector corresponds to the mean point of
the variables X, i.e.: (X1, X2,...Xn). The minimum distance vector from this point to the limit
state surface is perpendicular to the hyper-plane (a straight line in the two variables case)
tangent to the limit state surface in the design point (xd1,xd2,...xdn). The sensitivity coefficients
represent the influence of the corresponding variable in the failure probabilities.
The following equations hold for the sensitivity coefficients
i 1;
i2 = 1
(36)
Usually these values i, are taken as positives if the correspond to a resistance variable
and negatives in the case of effects of action variables.
When the limit state surface is highly non-linear, the error of substituting the surface
by the tangent hyper-plane in that point can be important. In those cases a minor error is
obtained if the limit state surface is substituted by the tangent quadratic surface. That is:
taking the square term in the Taylor's series. In this case the method is called Second Order
Reliability Method or SORM.
As it is said, the reliability index and the probability of failure pf are related by the
formula:
Pf = (-);
=-1(1- Pf );
(37)
-1
where () stands for the standard normal distribution and () for its inverse.
Example 2
a) Design of a steel beam
Simply supported beam: IPE 240 S235
g
Span
Cross section area:
Resistance modulus,
Yield stress
L = 6,0 m
A = 39,12 10-4 m2
W = 3 243 10-6 m3
fy = 235 MPa
Actions:
Permanent load: gk = 7,0 kN/m
Variable load: qk = 3,0 kN/m
Limit State function
Z(X) = 1 Wfy - 2 (g + q) L2 /8
In the following table all the variables are described by the mean, standard deviation
and the distribution function, see also [5]:
B - 16
Symbol
L [m]
W [m3]
1 [-]
fy [MPa]
2 [-]
g [kN/m2]
q [kN/m2]
Mean
6
3 243 10-6
1
280
1
0.007
0 .0008
Standard deviation
0
0
0,1
19,6
0,2
0.007 0,1
0. 00080,6
Distribution function
Deterministic
Deterministic
Lognormal
Lognormal
Lognormal
Normal
Gumbel
The following results are obtained using the attached MATHEMATICA notebook
Form.nb and package Level II.m.
Case a: Results
Reliability index: = 3,82 ; Probability of failure: (-) = 6,6710-5
Variable
Sensitivity coefficient
1
-0,392
fy
-0,275
2
0,778
g
0,304
q
0,270
From the results can be seen that the reliability index is well adjusted, but also that the
limit state is very sensitive to the actions model uncertainty (influence coefficient 0,778).
Case b:
Imagine that we can study more deeply the actions and we arrive to the conclusion
that the mean and variances of the actions are the same as before, but now we have reduced
the uncertainty model and we have now a coefficient of variation of 10%.
Action effects model
2 [-]
0,1
Lognormal
fy
Sensitivity coefficient
-0,505
-0,354
0,505
0,332
0,503
We can see that the influence coefficients are more equilibrated, but now the
reliability index is quite high.
Case c:
Perhaps we can reduce the steel profile. We can try with a smaller profile:
New steel profile: IPE 220, W= 252 cm3
Now the results are the following:
Reliability index: = 3,74; Probability of failure: (-) = 9,2910-5
B - 17
fy
-0,367
0,524
0,383
0,411
Example 3
Consider the slab of the example in Section 5 (see Figure 8). From the equilibrium
consideration was obtained:
0,8 fc x b = As fy
As fy(d 0,4 x) = M
Using equilibrium conditions (30) and (31), and considering uncertainty of the models
of resistances, 1, and actions, 2, the following formula for the limit state function can be
derived:
As f y
g(X) = 1 A s f y d
2 fcb
2 ( g + q )L 2 / 8
The following table shows the variables that are considered in the study using the
MATHEMATICA notebook FORM.nb:
Variable
Symbol
Resistance model
Reinforcement tensile strength
Reinforcement area
Effective height
Concrete compressive strength
Load model
Permanent load
Imposed load
Span
Mean value
1 [-]
1
fy [MPa]
560
As [m2] 0,00069
d [m]
0,23
fc [MPa]
30
2 [-]
1
7
g [kN/m2]
2
q [kN/m ]
0,8
L [m]
6
Standard
deviation
0,1
30
0,0000345
0,01
5,5
0,2
0.7
0,48
-
Distribution
function
Lognormal
Lognormal
Normal
Normal
Lognormal
Lognormal
Normal
Gamma
Deterministic
Results:
Reliability index: = 3,56; Probability of failure: (-) = 1,8710-4
Variable
1
As
fy
d
fc
2
g
q
Sensitivity coefficient -0,383 -0,193 -0,201 -0,177 -0,018 0,761 0,300 0,274
From the results it is possible to draw the following conclusions: the reliability index
is a little bit lower; the coefficient of influence of the concrete resistance, fc, is almost zero, it
could be possible to consider it as deterministic without influencing the results; and the
uncertainty in the actions has the maximum influence.
Case b:
The new study is performed considering reduced uncertainty in the load model as
follows:
Load model
2 [-]
B - 18
0,1
Lognormal
1
-0.500
As
-0,259
fy
-0,263
d
-0,235
fc
-0,023
2
0,500
g
0,347
q
0,431
Case c:
Now the reliability index is a bit higher. Reducing the cross section of the reinforcing
steel bars:
Reinforcement area
Load model
As [m2]
2 [-]
0,00059
1
0,000590,05
0,1
Normal
Lognormal
Results:
Reliability index: = 3,86; Probability of failure: (-) = 5,6110-5
Variable
1
Sensitivity coefficient -0,504
As
-0,260
fy
-0,266
d
-0,234
fc
-0,020
2
0,504
g
0,368
q
0,401
From both the examples can be seen that reducing the uncertainty in the model of
actions, is possible to reduce the steel section without reducing (even increasing) the overall
reliability.
7
SYSTEM RELIABILITY
7.1
General
Even in the simplest case of one structural element: a beam or column, more than one
ultimate limit state function has to be considered: failures at positive or negative moments or
shear of the beam. In most cases, the structure has multiple elements what is called a
'structural system'.
The reliability of the system depends on the reliability of its elements: the effects of
actions in each element depend on applied loads, loads and resistances may not be
independent, there may be a correlation between the properties of the elements in different
parts of the structure. Furthermore, there are limit states for the structure as a whole like the
overall deflection or foundation settlement.
When all the different failure modes are identified, a 'fault-tree' or an 'event-tree of all
the failures modes can be established.
Example 4
Consider the simple portal frame of the following Figure 4 submitted to the horizontal
and vertical loads Q. Assuming a plastic behaviour, the frame has three possible modes of
collapse. In each path different plastic hinges will be formed:
a)
b)
c)
B - 19
Q
1
The failure in any path implies the failure of the structure and the event 'structural
failure' will be the union of the all n failure modes (n = 3, in the example case). And therefore
the probability of failure of the structure will be:
Pf = P(FS) = P(F1 F2 ... Fn)
(38)
where Fi is the event failure in mode i. For each mode m nodes or structural elements have to
fail. So the probability of failure of each mode, Fi, will be:
Pf(Fi) = P(F1i F2i ... Fmi)
(39)
where Fji is the event failure of the j-th element or node in the i-th failure mode. To reach the
collapse of the structure in the mode i, m elements or nodes must fail.
In general the structure is idealized as a parallel system, a series system or a
combination of both.
7.2
Parallel systems
In the parallel system are the elements so interconnected that the failure of one or
more of the elements does not mean the failure of the whole structure. Such a structure is
called a redundant structure. This redundancy can be active, when the redundant members can
activate before the limit state of any member was reached, or passive, when the redundant
members only act when a limit state is reached in some member.
It should be taken into account that any hyperstatic structure is not necessary a parallel
system: if the elements are brittle, the failure of any of them can lead to a new distribution of
stresses such that new failures are reached immediately.
The failure of a pure parallel system with m components is given by:
Pf sys = P
I
1
Fj = P I (Zj < 0)
(40)
where Fj is event of the failure of the j-th component with Zj limit state function. Thus:
P(Fj) =P(Zj < 0) (-)
(41)
(42)
ij ik ;
with j, k = 1, 2, ... , m
(43)
and ij is the sensitivity factor of the i-th variable in the j-th component limit state function.
B - 20
IF )]
(44)
(45)
7.3
Series systems
A series system is that the failure of any element leads to the collapse of the structure.
It is called a 'weakest link'. Any statically determined structure is a series system. In a similar
manner that was shown in the previous section, the failure probability in a pure series system
with m components is given by:
C = P
U
1
Fj
(46)
(47)
m
m
Max P( F j ) Pf sys Min P( F j ),1 .
1
1
(48)
These bounds are usually too wide. A method to find more precise bounds can be
found in [11].
8
CONCLUDING REMARKS
B - 21
12.
B - 23
k := 1.65
d := 3.09
Sensitivity factors:
R := 0.8
E := 0.7
Reduced values:
R( ) := R
E( ) := E
kn( V) := ( 1 k V)
Lognormal distribution
kln( V) :=
Check:
exp( k) ln 1 + V
1+ V
Gumbel distribution
dn ( , V) := ( 1 R( ) V)
Lognormal distribution
dln( , V) :=
exp( R( ) ) ln 1 + V
1+ V
dgum ( , V) := 1 V ( 0.45 + 0.78ln( ln( 1 pnorm( E( ) , 0 , 1) ) ) )
Gumbel distribution
kn( V)
Rln( , V) :=
dn ( , V)
kln( V)
1.5
1.5
Rn( , 0.05)
Rn( , 0.1)
dln ( , V)
Rln( , 0.05)
Rln( , 0.1)
Rn( , 0.15)
Rln( , 0.15)
0.5
0.5
B - 24
Gn( , 0.05)
1.4
Gn( 3.3 , V)
Gn( , 0.1)
Gn( , 0.15)
Gn( 3.8 , V)
1.2
Gn( 4.3 , V)
0.2
0.4
V
dgum( , V)
kgum( V)
Qgum( , V) :=
Partial factor Q:
2
Qgum( , 0.2)
Qgum( 3.3 , V)
Qgum( , 0.3)
Qgum( 3.8 , V)
Qgum( , 0.4)
Qgum( 4.3 , V)
1.4
1.2
0.2
0.4
V
MH, 21.2.2003
B - 25
B - 26
B - 27
B - 28
B - 29
B - 30
B - 31
6,00
7,00 gammaG=
3,00 gammaQ=
62,78
1,35
1,5
L
Rectangular cross-section
Dimensions
Factors
Concrete
Rebars
Reratio
=x/d<
Estimate
b= 1
d= 0,17
c= 1,5
cc= 1
s= 1,15
fcd=acc*fck/c= 13,3
fctm= 2,2
fck [MPa] = 20
fyd=fyk/= 434,8
fyk [MPa] = 500
min [%] =
0,130
max
0,45
x
z~
0,9 d
As~Md/(z*fyd)=
0,000944
A[m^2]= 0,000933
As
[%] =
0,55
>min ? PRAVDA
b
=x/d=
0,22
<max?
1,10
PRAVDA max[%] =
General Table
m = Md/
/(b*d^2*fcd)
0,00
0,05
0,10
0,15
0,20
0,25
0,30
0,35
0,40
=
1,8
Reratio
[%]
0,00000
0,15737
0,32376
0,50091
0,69124
0,89821
1,12714
1,38698
1,69521
As
f
= cd
f yd
bd
1,00
0,97
0,95
0,92
0,89
0,85
0,82
0,77
0,72
s
[%]
5,11
2,30
1,36
0,89
0,61
0,41
0,27
0,16
As
m^2
0,00000
0,00027
0,00055
0,00085
0,00118
0,00153
0,00192
0,00236
0,00288
c=0,35%
d
Md
[kNm]
0,0
19,3
38,5
57,8
77,1
96,3
115,6
134,9
154,1
f cd bd 2
1,6
1,4
1,2
1,0
0,8
0,6
0,4
0,2
0,0
0,00
0,05
0,10
0,15
0,20
0,25
M/(fcd*b*d^2)
B - 32
0,30
0,35
0,40
1.35
1.5
L
Rectangular cross-section
Dimensions
Factors
Concrete
Rebars
Reratio
=x/d<
Estimate
6.00
7.00 gammaG=
3.00 gammaQ=
62.78
b= 1
d= 0.17
c= 1.5
cc= 1
s= 1.15
acc*fck/c= 13.3
fck [MPa] = 20
fctm= 2.2
yd=fyk/= 434.8
fyk [MPa] = 500
min [%] =
0.130
max
0.45
x
z~
0,9 d
As~Md/(z*fyd)=
0.000944
A[m^2]= 0.000933
As
[%] =
0.55
>min ? PRAVDA
b
=x/d=
0.22
<max?
PRAVDA max[%] =
1.10
c=0,35%
d
General Table
m = Md/
/(b*d^2*fcd)
0.00
0.05
0.10
0.15
0.20
0.25
0.30
0.35
0.40
Reratio
[%]
0.00000
0.15737
0.32376
0.50091
0.69124
0.89821
1.12714
1.38698
1.69521
As
f
= cd
bd
f yd
1.8
=x/d
/d = 1 - 0,4*
0.00
0.06
0.13
0.20
0.28
0.37
0.46
0.57
0.69
1.00
0.97
0.95
0.92
0.89
0.85
0.82
0.77
0.72
s
[%]
5.11
2.30
1.36
0.89
0.61
0.41
0.27
0.16
As
m^2
0.00000
0.00027
0.00055
0.00085
0.00118
0.00153
0.00192
0.00236
0.00288
Md
[kNm]
0.0
19.3
38.5
57.8
77.1
96.3
115.6
134.9
154.1
f cd bd 2
1.6
1.4
1.2
1.0
0.8
0.6
0.4
0.2
0.0
0.00
0.05
0.10
0.15
0.20
0.25
M/(fcd*b*d^2)
B - 33
0.30
0.35
0.40
Data
Packages needed
$TextStyle
Data
data = {1.3, 4.3, 1.3, 5.4, 3.7, 3.8, 4., 2.9, 3.2, 4.5, 4., 3.4, 2.4,
1.8, 1.7, 2.2, 4.1, 2.6, 4.1, 3.3, 3.5, 3.7, 2.4, 2.8, 2.5, 3., 3.3, 2.6,
2.9, 2.4, 2.6, 2.9, 2.8, 3.1, 3.2, 3.4, 3.5, 3.2, 3.3, 3.};
distributions = 9"normal", "lognormal", "gamma"<;
numdistributions = LengthAdistributionsD;
Descriptivestatistic
ndata = LengthAdataD
mean = MeanAdataD
variance = VarianceAdataD
deviation = StandardDeviationAdataD
deviation
coefvariation =
mean
40
3.1025
0.718199
0.847466
0.273156
Distribution functions
Parameters of distribution
B - 34
!!!
deviation 6
par2 = NB F;
Quantiles
DoA
PrintAdistAiDD;
PrintA
NAQuantileAdistAiD, 0.5D, 4D, " ",
NAQuantileAdistAiD, 0.75D, 4D, " ",
NAQuantileAdistAiD, 0.9D, 4D, " ",
NAQuantileAdistAiD, 0.95D, 4D, " ",
NAQuantileAdistAiD, 0.99D, 4D, " ",
NAQuantileAdistAiD, 0.999D, 4D, " ",
NAQuantileAdistAiD, 0.9999D, 4DD;
PrintA" "D,
9i, 1, numdistributions<D
N o r m a l D i s t r i b u t i o nA3 . 1 0 2 5, 0 . 8 4 7 4 6 6D
3 . 1 02 5 3 . 6 7 4 1 1 4 . 1 88 5 7 4 . 4 9 6 4 6 5 . 0 7 4 5 . 7 2 1 3 7 6 . 2 54 2 4
L o g No r m a l D i s t r i b ut io nA1 . 0 9 6 2 3, 0 . 2 6 8 2 5 7D
2 . 9 92 8 5 3 . 5 8 6 4 5 4 . 2 20 7 5 4 . 6 52 8 1 5 . 58 6 1 2 6 . 8 5 6 5 4 8 . 1 1 6 33
B - 35
3 . 0 25 6 9 3 . 6 2 4 5 4 4 . 2 26 4 . 6 1 5 1 8 5 . 4 0 5 9 6 . 3 8 8 9 4 7 . 2 76 9 8
frecuenciesPiT
DoBAppendToBhystogr, 9intervalsPiT, =F;
total HintervalsPi + 1T . intervalsPiTM
frecuenciesPiT
AppendToBhystogr, 9intervalsPi + 1T, =F,
total HintervalsPi + 1T . intervalsPiTM
9i, numintervals<F
AppendToAhystogr, 9intervalsPnumintervals + 1T, 0.<D;
hystogram = ListPlotAhystogr, PlotJoined True, DisplayFunction IdentityD;
Degrees of freedom = 6
C o nf i d e n c e le v e l bet w e e n 9 5 . 9 0 6 2 a n d 9 9 . 5 7 8 2 %
B - 36
L o g N o r m a l D i s t r i b u t i o nA1 . 0 9 6 2 3, 0 . 2 6 8 2 5 7D
E s ti m a t o r W = 1 . 3 1 8 28
D e g r e e s o f f re e d o m = 6
C o n f i de n c e l e v e l b et w e e n 8 5 . 8 2 6 7 a n d 9 7 . 0 6 4 %
lognormal
0.5
0.4
0.3
0.2
0.1
2
GammaDistributionA13.4023, 0.23149D
Estimator W = 1.00414
Degrees of freedom = 6
gamma
0.5
0.4
0.3
0.2
0.1
2
Distribution functions
Distribution of the data points
datosord = SortAdataD;
i
probab = TableB , 9i, 1, ndata<F;
ndata + 1.
loglogprob = .LogA. LogAprobabDD;
B - 37
AspectRatio 1.2,
GridLines 9None, 90.95<<,
DisplayFunction $DisplayFunctionD;
B - 38
0.8
0.6
0.4
0.2
.999
.99
.95
.9
.5
.1
.01
.001
2
B - 39
FORM
This notebook compute the reliability index, failure probability and influence factors in level II, using the package
"Reliability`Level2". In this package those variables are determined through the algorithm "Normal Tail Approxima
tion" as is explained in the book of Madsen et al.: Methods of Structural Safety, pp. 94 and following.
The failure function of the limit state must be defined and, also, the independent basic variables given by a matrix
with a row for each variable with the kind of distribution function assumed, the mean and the standard deviation.
The following names of distribution functions are implemented:
"normal" or "gauss",
"lognormal" or "lognor",
"gamma",
"logistic",
"uniform",
"weibull", and
"extremes" or "gumbel".
The iterations defined in the loop of the algorithm is stopped when the differences between two consecutive values i
less than the error defined.
Needs["Reliability`Level2`"]
LImit State functio
Define the limit state function
gAz_D :=
zAA1DD HzAA2DD zAA3DDM HzAA4DD . HHzAA2DD zAA3DDM H2 zAA5DD bMMM .
zAA6DD HzAA7DD + zAA8DDM L2 8
Deterministic
variable
b = 1;
L = 6.;
Ramdomvariable
Define the variables matrix: distribution function, mean and standard deviation
m = {{"lognormal",
1,.1},
{"normal", 0.00082 ,0.00082*0.05},
{"lognormal",
560,560*0.054},
{"normal",0.22,0.01}
,
{"lognormal",
30,30*0.167},
{"lognormal",
1.2,.18},
{"normal", 0.007,0.007*.1},
{"gumbel", .0008,.0008*0.6
}
};
Failprob[m,g]
Obtain the reliability index beta ; the failure probabilities; and the sensitivity factor for each variable
b e t a = 3 . 8 0 0 0 4, P h iH .b e t aM = 7 . 2 3 3 6 7 1 0 .5 H 5 i t e r a t i o n sM
a l p h a = 9 .0 . 4 4 2 3 8 8, . 0 . 2 2 3 4 5 5, .0 . 2 3 0 8 9 8, . 0 . 2 1 6 8 0 7, .0 . 0 2 5 8 5 8 5, 0 . 6 6 1 5 3 5, 0 . 3 3 4 9 9 6, 0 . 3 2 1 7 4 7<
B - 40
(* Definition of Grad*)
Failprob[matrix_,g_]:=
Module[
{
n=First[Dimensions[matrix]],
znew=muiter=Transpose[matrix][[2]],
sigmaiter=Transpose[matrix][[3]],
distnor=NormalDistribution[0.,1.],
error= 10.^(-3),
(* error admited *)
betanterior=0,
ind=0,
z,aux,gradiente,gradpart, ceta,muzeta,sigmazeta,alpha,beta,prob },
z=Array[a,n];
gradiente=Grad[g[z],z] ;
Characterization[matrix]; (* Characterization of the variables *)
Do[
(* iterative cycle *)
ziter=znew;
partic=Take[{z[[1]]->ziter[[1]],z[[2]]->ziter[[2]],z[[3]]->ziter[[3]],
z[[4]]->ziter[[4]],z[[5]]->ziter[[5]],z[[6]]->ziter[[6]],
z[[7]]->ziter[[7]],z[[8]]->ziter[[8]],z[[9]]->ziter[[9]],
z[[10]]->ziter[[10]]},n];
(* iteration values *)
Do[
aux = invnormal[CDF[dist[i],ziter[[i]]]];
sigmaiter[[i]] = PDF[distnor,aux]/PDF[dist[i], ziter[[i]] ];
muiter[[i]] = ziter[[i]] - aux sigmaiter[[i]],
B - 41
B - 42
B - 43
E := 50
wE := 0.2
C( ) :=
x0( , , ) :=
2
+ 4+
3
C( )
E := 1.14
R := 3 wR + wR
R := wR R
Check:
2
+ 4
C( 0) = 0
14
if 0
x0( R , R , R) = 8.527 10
x0( E , E , E) = 22.522
6 otherwise
E := wE E
3
s ( ) :=
ln 1 + C( )
R( x) :=
n =
xR xR
xE xE
3.304
for i 1 .. n
2.904
Re
2.904
R( xR)
2.904
pf := pnorm ( , 0, 1)
Re xR Re qnorm( R( xR) , 0 , 1)
dnorm( qnorm( E( xE) , 0 , 1) , 0 , 1)
Ee
pf n =
E( xE)
4.7710-4
Ee xE Ee qnorm( E( xE) , 0, 1)
aR
aE
1.84310-3
Re Ee
1.84410-3
( Re2 + Ee2)
0.5
1.84410-3
Re
2
2
Re + Ee
3.4
0.5
Ee
3.2
xR Re + aR Re
2.8
xE Ee + aE Ee
3
n
B - 44
Distribution parameter C
given by the skewness :
+4+
C( ) :=
x0( , , ) :=
Distribution bound x0
( - 6 for zero ):
Check:
+ 4
C( 0) = 0
C( )
if 0
x0( 0 , 1 , 1) = 3.1038
6 otherwise
Standardised variable:
ln u ( x, , ) +
uu ( x, , , ) :=
+ ln C( ) 1 + C( ) 2
if 0
C( )
sign ( ) ln 1 + C( )
u ( x, , ) otherwise
( x, , , ) :=
u ( x, , ) +
C( )
2
ln( 1 + C( ) )
dnorm( uu ( x, , , ) , 0 , 1)
if 0
otherwise
( x, , , ) := pnorm( uu ( x, , , ) , 0 , 1)
Distribution function:
1
3
= 0 then the normal
distribution is used . Index := 4 a := 1
When X1 and X4 are
5
1
model uncertainties
1
then LN(,) is used.
6
7
15
1
0.05
10
1
:=
0.05
10
5
0.01
4
x0 = 8.8818 10
x0 = 8.3195 10
x := i
i
0
1.
100
10
:=
1
50
10
0.01
4
x :=
1
x0 = 0
5
4 4( 5 5
6 6
a (a x + a x )
1 2 2
3 3
0
0.15
0.301
0.301
:=
0.15
0.608
1.14
0
x0 = 3.7388
7 7)
a + a x a x + a x + a x
B - 45
i := 1 .. 7
x0 := x0 i , i , i
i
x0 = 8.3195 10
1
14
x0 = 8.5265 10
2
x0 = 0.05
7
x = 0.5455
1
187
4. FORM iteration procedure: Probability of failure pf is determined from the reliability index
Number of iterations
The initial guess
values of X
The value x2 is
calculated
from g (X)= 0
n :=
n := 1 .. 5
x
x
1
(5
a + a x a x + a x + a x
4 4
(2
a a x + a x
1
7 7
6 6
3 3
for j 1 .. n
for i 1 .. 7
( (i
) ) )
( x , i , i , i)
i
ei x ei qnorm ( ( x , i , i , i) , 0 , 1)
i
i
Equivalent
normal
distributions of
the basic
variables X.
x ei
i
u
i
Standardised variables
dnorm qnorm x , i , i , i , 0 , 1 , 0 , 1
ei
ei
(2
g a a x + a x e1
1
Derivatives of g(X):
3 3
g a a x e2
2
1 2 1
g a a x e3
3
1 3 1
(5
g a a x + a x e4
4
6 6
g a a x e5
5
4 5 4
g a a x e6
6
4 6 4
g a a x e7
7
Reliability index:
4 7 4
( g u )
Sensitivity factors:
(gg)
0.5
for i 1 .. 7
g
aa
i
Iteration of the
reliability index
Probability of
failure pf
(gg)
0.5
x ei aa ei
i
i
x
1
(5
a + a x a x + a x + a x
4 4
(2
6 6
a a x + a x
1
3 3
7 7
0
2.9022
2.9303
=
2.9292
2.929
2.9289
B - 46
pf := pnorm 5 , 0, 1
pf = 1.7006 10
A
i
0
1
2
3
4
5
6
7
0
1
2
3
4
5
6
7
B - 47
B - 48
B - 49
B - 50
B - 51
B - 52
B - 53
B - 54
B - 55
Summary
The basic European standard for design of buildings and other engineering works,
EN 1990 Basis of structural design, provides alternative design procedures, for which
national choice is allowed. One of the most important questions concerns three fundamental
combinations of actions for persistent and transient design situations in the Ultimate limit
states. Simple example of generic structural member shows, that the alternative load
combinations may lead to considerably different reliability levels. Probabilistic methods of
structural reliability theory are used to identify characteristic features of each combination and
to formulate general recommendations. It appears that further calibration studies concerning
structures made of different materials are needed during the examination period of EN 1990
in order to analyse all possible consequences of national choice.
1
INTRODUCTION
1.1
Background materials
Each part of Eurocodes, including basic document EN 1990 [1], contains a number of
the Nationally Determined Parameters (NDP) for which national choice is allowed. In
accordance with the Guidance paper L concerning the Construction Products Directive an
important two years period after date of availability of each Eurocode Part is allowed to fix
the NDPs. However it is expected that calibration will continue during the coexistence period,
which starts at the end of the National calibration period and lasts up to three years after the
national publication of the last Part of a Package.
1.2
General principles
Basic concepts of code calibration are mentioned in Annex C of EN 1990 [1], in the
International Standard ISO 2394 [2] and ISO 13822 [3]. Additional information may be found
in the background document developed by JCSS [4] and in recently published handbook [5]
to EN 1990 [1]. Guidance for application of probabilistic methods of structural reliability may
be also found in working materials provided by JCSS [6] and in relevant literature listed in [5]
and [6].
In general NDPs may be calibrated either by direct comparison or by probabilistic
methods. Results of both approaches are usually combined with judgement (as mentioned in
ISO 2394 [2]). In this study probabilistic approach is applied mainly, a direct comparison of
load effects is shortly described in Appendix A to this contribution. Note that for the
probabilistic calibration software products [7,8,9] can be used. In particular the programme
[9] is intended for calibration purposes. Special purpose MATLAB functions and
MATHCAD sheet attached to this Annex may be also used for calibration studies.
C-1
(1)
(2)
Ed = G Gk + Q Qk + W W Wk
(3)
The less favourable action effect from (2) and (3) should be considered. In equation
(3) is a reduction factor for unfavourable permanent actions G. Note that in equations (1) to
(3) + generally implies to be combined with.
C. In addition EN 1990 [1] allows further modification of alternative B, simplifying
equation (2) by considering permanent loads only, thus the load effect is then
Ed = G Gk
(4)
The less favourable action effect resulting from (3) and (4) is then considered. In
addition to the combinations A, B, C provided in EN 1990 [1] (for recommended values G =
1,35, Q = 1,5) an additional combination may be also considered in the analysis to illustrate
the sensitivity of the resulting reliability level to partial factors, and the possible effect of their
reduction.
If the leading action is wind W, then in equations (1) and (2) instead of reducing wind
action W by factor W, the imposed load Q should be reduced by the appropriate factor Q.
C-2
= (Qk+Wk)/(Gk+Qk+Wk), k = Wk/Qk
(5)
Note that a realistic range of is from 0,1 to 0,6. However in some cases the load ratio
may be very low if not zero (e.g. underground garages).
For a given design value of the load effect Ed the characteristic values of individual
actions Gk, Qk, Wk can be expressed using variables and k as follows
Gk =
Ed
(( Q ) Q + k ( W ) W )
( ) G +
(1 + k )(1 )
, Qk =
Gk
, Wk = k Qk
(1 + k )(1 - )
(6)
The factors , G and Q indicated in the first relationship of (6) in brackets are applied
in the same way (either yes or no) as in equations (1) to (4) for alternative combination rules
A, B and C.
For alternative A, equation (1) is valid in the whole range 0 1, whereas using
alternative B, equation (2) is valid in the interval 0 lim,B and equation (3) in the interval
lim,B 1. Correspondingly, for alternative C equation (4) is valid in the interval 0
lim,C and equation (3) in the interval lim,C 1. The limiting values lim,B and lim,C can be
derived from equations (2) to (5) as follows
lim,B=
G (1 )(1 + k )
G (1 )(1 + k ) + Q (a Q ) + W k (b W )
(7)
G (1 )(1 + k )
G (1 )(1 + k ) + Q a + W kb
(8)
lim,C=
Thus, the National Annexes should include the recommendation of one of the
alternatives indicated in EN 1990 [1] for a fundamental combination of actions in the
Ultimate limit states and partial factors G and Q for permanent and variable actions.
Considering a generic structural member it will be shown that the choice of these nationally
determined parameters may significantly affect the resulting reliability level. Partial and
C-3
Permanent G
Imposed Q
Climatic W
Combination factor
Reduction factor
1,0
0,7
0,6
0.85
-
1,35
1,5
1,5
In addition to the factors indicated in Table 1 other values will be used to make comparison of
Eurocode procedures with some national rules.
3
(9)
where R denotes the global resistance factor (commonly expected to be within the
range from 1 to 1,2). The significance of both values Rk and Rd is obvious from Figure 2,
where the random variable R is described by the probability density function R (R), and the
design value Rd is indicated as a particular value of R corresponding to a certain small
probability p of being violated.
R(R)
5%
Rd
Rk
Figure 2. Random variable R, the characteristic value Rk and design value Rd.
In design calculation of a structural member the design value Rd of the resistance R is
normally obtained by substituting design values Xdi for the random variables Xi, thus
Rd = R(Xdi)
C-4
(10)
Middle value
1,15
1,10
0,15
1,28
Range
1,0 1,3
1,0 1,20
0,10 0,25
1,10 1,40
Note 1. The coefficient of variation VR includes the variability of the model uncertainty assumed to have the
coefficient of variability 0,05.
Note 2. The values are different for BSI codes reflecting the fact that lower values of R are used for particular
materials.
A middle values for the global safety factor R = 1,15, for the coefficient of variation
VR = 0,15 and for the mean ratio = R/Rk = 1,28 are considered in the following example of
a code condition.
4
4.1
(11)
where factor R represents uncertainties of the resistance model R0(X) and factor E represents
uncertainties of the load effect model E0(X). The vector X denotes all the basic variables
entering the expressions for the resistance R(X) and the load effect E(X). Taking into account
general expressions (1) to (4) the load effect E0(X) may be written as
E(X) = E (G0 + Q0 + W0)
(12)
Considering the limit state function given by equation (11) and expression (12) giving
the load effect, it follows that basic variables R, G, Q, and W covering effects of model
uncertainties are defined as follows
R = R R0(X), G= E G0, Q= E Q0, W= E W0
(13)
Taking into account equation (13), the limit state function (11) may be written in a
simple form as
Z(X) = R (G + Q+ W)
C-5
(14)
(15)
It should be noted that normally (due to several reasons) the design resistance Rd(X) is greater
than the design load effect Ed(X), which may provide additional safety margin not considered
here.
Assuming a certain set of partial and combination factors , , and , the design
expression (15) can be used to specify the characteristic values Xk of each basic variable X.
The probabilistic characteristics (the mean, standard deviation) of each basic variable X can
be then related to its characteristic value Xk as indicated in Table 3.
Table 3. Probabilistic models of basic variables for time invariant reliability analysis using
Turkstra's rule (combination of 50-year maximum of leading action and an annual maximum
of accompanying action).
St.dev.
No. Category of
Name of basic
Sym. Dim- DistriMean
variables
variables
X ension bution
X
X
1
Actions
Permanent
G0
kN
N
Gk
0,1X
2
Imposed - 5 years
Q0 kN/m2
GU
0,2Qk
1,1X
GU
0,6Qk
2
Imposed - 50 y.
Q0 kN/m2
0,35X
3
Wind - 1 year
W0 kN/m2
GU
0,3Wk
0,5X
4
Wind - 50 year
W0 kN/m2
GU
0,7Wk
0,35X
5
Resistance
Resistance
R
kN/m2
LN Rk +1.65R 0,15X
6
Uncertainty
Uncertainty
LN
1
0.05
E
Probabilistic models indicated in Table 3 are based on data available in the
recommendation of JCSS [4,6] and literature [11,12,13,14]. As mentioned above the
probabilistic characteristics indicated in Table 3 represent just conventional models that might
be slightly conservative.
Note that the mean of a resistance R indicated in Table 3 in terms of the characteristic
value Rk and the standard deviation R may be assessed assuming a given coefficient of
variation VR using relationship
R = Rk exp(1,65 VR)
(16)
Under this assumption the mean resistance factor considered in Table 2 is given as
= R / Rk = exp(1,65 VR)
(17)
Considering the coefficient of variation VR = 0.15, the mean resistance factor becomes
= 1.28 as indicated in Table 2.
It should be emphasised that the probabilistic models of basic variables indicated in
Table 3 are primarily intended as "conventional models" in time invariant reliability analysis
C-6
Reliability measures
The probability of failure Pf is the basic reliability measure used in this study. It can be
expressed on the basis of a limit state (performance) function Z(X) defined in such a way that
a structure is considered to survive if Z(X)> 0 and to fail if Z(X) 0. An example of the
function Z(X) is given by equation (14). In a general case the failure probability Pf can be
determined using the integral
( X )dX
Pf = Prob(Z 0) =
(18)
g
g ( X ) 0
where g(X) denotes joint probability density distribution of the basic variable X,
which may not be , however, available.
Assume that both the resistance R(X) and the load effect E(X) represent a single
variable X used to analyse structural performance (e.g. axial force or bending moment that is
represented by R(X) and E(X)). Then the integration indicated in expression (18) may be
simplified and the probability Pf can then be expressed as:
Pf = Prob(Z(x) 0) = E ( x ) R ( x )dx
(19)
where E(x) denotes the probability density function of E(X), R(x) the distribution of R(X).
To use equation (19) both the probability density function E(Z) and the distribution function
R(x) must be known (at least in an approximate form). Simplified procedure based on
expression (19) is used in this study.
Note that there are commercially available software products (e.g. VaP, COMREL),
which can be used to determine the failure probability Pf in more complicated cases than
considered here (when expression (19) cannot be used). These software products were used in
this study to check results obtained by numerical integration based on expression (19).
In Annex C of EN 1990 an alternative measure of reliability is conventionally defined
by the reliability index , which is related to Pf as
Pf = ( )
(20)
where is the cumulative distribution function of the standardised normal distribution. The
relation between Pf and is indicated in Table 4.
Table 4. Relation between and Pf.
Pf
10-1
1,28
10-2
2,32
10-3
3,09
10-4
3,72
C-7
10-5
4,27
10-6
4,75
10-7
5,20
Table C2 of EN 1990 recommends for the ultimate limit state of buildings over a fifty
year design working life a target value of reliability index t = 3,8. If one year period is
considered in reliability verification, then t = 4,7. Both the equivalent reliability measures,
the failure probability Pf and the reliability index , are used in this study.
4.4
Sensitivity factors
Sensitivity factors of the First Order Reliability Methods (FORM) are normally used
[1,2] to calibrate design values of basic variables and partial safety factors. Considering the
limit state function Z(X) (reliability margin) given by equation (11), the sensitivity factors for
the four cumulative variables R, G, Q, W can be defined in terms of their standard deviations
R, G, Q and W as follows
R =
, G = G , Q = Q , W = W
g
g
g
g
(21)
g = R2 + G2 + Q2 + W2
(22)
5.1
C-8
A
B
C
B
B
A
R
W
Q
Figure 3. Variation of the reliability index , the failure probability Pf and the sensitivity
factors R, E, G, Q and W with the load ration for k = 0, for a generic cross-section
assuming R = 1,15 and the coefficient of variation VR = 0,15.
It follows from Figure 3 that for the assumed higher coefficient of variation VR = 0,15
only the combination rule A (i.e. expression (6.10) of EN 1990) [1] seems to be fully
acceptable ( > 3.8 and Pf.< 7,2310-5) in the interval 0 < < 0.8, however the reliability
level considerably varies with . In some cases the alternative A might lead to an uneconomic
design.
The alternative B (i.e. expression (6.10a) and (6.10b) of EN 1990) is acceptable in a
slightly shorter range of , 0 < < 0.7 than the variant A but provides obviously much more
uniform distribution of reliability level with . Obviously it would lead to a more economic
C-9
C - 10
A
B
B
C
C
B
B
A
G
E
Figure 4. Variation of the reliability index , the failure probability Pf and the sensitivity
factors R, E, G, Q and W with the load ratio for k = 0.75, for a generic cross section
assuming R = 1,15 and the coefficient of variation VR = 0,15.
It follows from Figure 4 that for the assumed coefficient of variation VR = 0,15 and the
consideration of two variable actions the reliability of the generic cross-section exposed to
two variable actions is considerably greater than the reliability of the same cross-section
exposed to one variable action only. This finding also indicates that the factor W may be
rather high. Note that the sensitivity factors R seems to be slightly greater than the values R
= 0,8 considered in EN 1990 [1] and E in absolute value is less than E = - 0,7 recommended
in [1]. This finding depends on assumed variability of basic variables.
C - 11
Figure 5. Variation of the reliability index with the load ratio and partial factor for
resistance R and k = 0 (i.e. imposed load Q is the only variable action), for the generic-cross
section assuming partial safety factors G = 1,35 and Q = 1,5, and the coefficient of variation
VR = 0,15.
It follows from Figure 5 that for the assumed variables the acceptable domain of the
load ratio and the coefficient of variation VR is limited by the contour line determined as an
intersection of the surface and the plain = 3,8 in Figure 5. Obviously with increasing R
reliability index increases, R = 1.15 would be satisfactory for most of the practical range of
the load ratio (for the load ratio .< 0,8).
C - 12
CONCLUDING REMARKS
The newly available EN 1990 provides alternative design procedures and parameters
that should be unambiguously specified in the National Annexes of Member States of CEN.
These alternative design procedures lead in some cases to significantly different reliability
levels. Preparation of National Annexes is therefore a complicated task for each Member
State. Furthermore, the Eurocode standards recognise the responsibility of the regulatory
authorities in each Member State and safeguard their right to determine values related to
regulatory safety matters at national level.
Simple examples of a generic structural member confirm the results of the earlier
studies that the reliability of structures, designed according to the alternative combination
rules provided in EN 1990 by expressions (6.10), (6.10a) and (6.10b), may vary considerably.
Expression (6.10) leads to the most reliable but in some cases uneconomical structures. Twin
expressions (6.10a) and (6.10b) provide a lower but comparatively most uniform reliability
level for all load ratios. Moreover, they seem to fully comply with EN recommendations
(reliability index 3,8 for a 50-year time period). The lowest reliability is obtained from the
third alternative, given by modified expression (6.10a) and expression (6.10b). This
alternative seems to lead to a rather low reliability level, particularly for structures exposed
mainly to a permanent load.
In order to make an unambiguous recommendation for National Annexes to EN 1990,
further investigations are urgently needed. Obviously more complicated structural elements,
made of various materials, should be analysed and compared. Such a calibration activity
should preferably be organised on an international level. The short-term objective of these
activities should be to develop the necessary background materials for preparation of the
National Annexes. The long-term objective should be to further harmonization of the
alternative design procedures considered during the next revision of the present generation of
Eurocodes.
It is expected that further calibration studies concerning structures made of different
materials will be needed during the examination period of EN 1990 (next few years) in order
to analyse all possible consequences of national choice.
REFERENCES
[1] EN 1990 Eurocode - Basis of structural design. CEN 2002.
[2] ISO 2394 General principles on reliability for structures, ISO 1998.
[3] ISO 13822. Basis for design of structures - Assessment of existing structures, ISO 2001.
[4] JCSS: Probabilistic model code. JCSS working materials, http://www.jcss.ethz.ch/, 2001.
[5] Gulvanessian, H. Calgaro, J.-A. Holick, M.: Designer's Guide to EN 1990, Eurocode:
Basis of Structural Design; Thomas Telford, London, 2002, ISBN: 07277 3011 8, 192 p.
[6] JCSS: Background documentation, Part 1 of EC 1 Basis of design, 1996.
[7] VaP, Variable Processor, version 1.6, ETH Zurich, 1997.
[8] COMREL, version 7.10, Reliability Consulting Programs, RCP MUNICH, 1999.
[9] CodeCal, Excel sheet developed by JCSS, http://www.jcss.ethz.ch/.
[10] Turkstra, C.J.: Application of Bayesian Decision Theory. Study No. 3 Structural
Reliability and Codified Design. Solid Mechanics Division, University of Waterloo,
Ontario, Canada, 1970.
[11] Holick M., Markov J.: Verification of load factors for concrete components by
reliability and optimization analysis: Background documents for implementing
C - 13
C - 14
E = Ed / (Gk + Qk + Wk)
(A.1)
where the design load effect Ed is given by one of equations (1) to (4) depending on the
combination rule considered (for example combination rules A, B or C). It follows from
equation (1) to (4) and (A.1) that in general deterministic E may be expressed as
(A.2)
where the factors in brackets (), (Q) and (W) are applied in accordance with the principles
of appropriate combination rule. For example assuming that Q is the leading variable load and
W is accompanying load, the combination rule A based on expression (6.10) of EN 1990 [1]
the global factor E follows from (1) and (A.2) as
E = (1 - ) G + (Q + k W W) / (1+ k)
(A.3)
Similarly the global factors E of other combination rules B and C may be obtained from
general expression (A.2). It follows from (2) that equation (A.2) becomes
E = (1 - ) G + (Q Q + k W W) / (1+ k)
(A.4)
E = (1 - ) G + (Q + k W W) / (1+ k)
(A.5)
E = (1 - ) G
(A.6)
C - 15
E
1.5
A
1.4
B
1.3
1.2
1.1
lim,C
1
lim,B
0.2
0.4
0.6
0.8
Figure A.1. The global load factor E for the combination rules A, B and C assuming k = 0.
E
1.5
1.4
A
1.3
B
C
1.2
C
1.1
lim,C
1
0.2
lim,B
0.4
0.6
0.8
Figure A.2. The global load factor E for the combination rules A, B and C assuming k = 0,75.
C - 16
(A.7)
The probability p is obviously dependent on the global load factor E determining the load
effect Ed as follows from equation (A.1). When probabilistic models of actions G, Q and W
are known, then for a given value of E the probability p may be determined. Let us remind
that in accordance with the principles of EN 1990 [1] (considering E = -0.7 and = 3.8) the
recommended value of the probability p given by equation (A.4) is
pE = (E ) = (- 0,7 3,8) = 0,004
(A.8)
E
p=0.001
1.5
A
B
1.4
p=0.002
1.3
p=0.004
C
1.2
C
1.1
lim,C
0
0.2
lim,B
0.4
0.6
0.8
Figure A.3. The global load factor E for the combination rules A, B and C assuming k = 0
and theoretical values of E corresponding to selected probabilities of E exceeding Ed
C - 17
1 p
Figure A.4. Variation of the theoretical and deterministic global factor E with the load ratio
and the probability p assuming the combination rule A.
Variation of the global factor E with the load ratio and the probability p clearly
indicates that the combination rule A, represented in Figure A.4 by a plane, is rather safe (and
perhaps uneconomic) substitution of the theoretical (probabilistic) values.
C - 18
(B.1)
where E = 0,8 3,8 = 3,04 as recommended in EN 1990 [1]. Note that for VR = 0,10
equation B.1 yields the partial factor R = 1,15.
Figure B.1 shows the variation of the partial factor R with the coefficient of variability
VR together with corresponding reliability index determined taking into account the partial
factor R as a function of VR given by equation (B.1). In Figure B.1 the combination rule A
and a generic cross section are considered only.
5
=0.3
=0.6
4.5
4
=0.0
3.5
3
2.5
2
1.5
1
wR
0
0.1
0.2
0.3
0.4
0.5
Figure B.1. Variation of the reliability index of a generic cross section with the coefficient
of variability VR for selected load ratios assuming the partial factor R as a function of VR.
Figure B.1 indicates that if the partial factor R is considered as a function of the
coefficient of variability VR, the effect of resistance variability is not essential. Considering a
realistic range of resistance variability 0,05 < VR < 0,25, differences in values seem to be
about 0,5. This finding justifies the concept of a generic cross section used in reliability
analysis of alternative load combinations.
C - 19
R
E
G
Q
W
()
G
Q
W
Q
W
C - 20
ATTACHMENTS
C - 21
7. MATHCAD sheet"Generic.mcd"
MATHCAD function "Generic" is intended for investigation of the combination rules
provided in EN 1990. A general structural member of the resistance r (including model
uncertainty) is considered.
8. MATHCAD sheet"Load effect.mcd"
MATHCAD sheet "LoadEffect" is intended for investigation of combination rules
provided in EN 1990 by expressions (6.10), (6.10a) and (6.10b) considering three loads: G, Q
and W. Turkstra's rule ( 50-year extremes of a leading and annual extremes of an
accompanying action) is applied.
C - 22
C - 23
C - 24
C - 25
C - 26
C - 27
C - 28
C - 29
C - 30
C - 31
C - 32
C - 33
Constants:
:= 0.85 Q := 0.70
W := 0.60
Normal distribution of G:
Gumbel distribution of Q:
k := 0.0
Q = mQ*Qk
Rd := 1
G = Gk
wG := 0.1
mQ := 0.6
wQ := 0.35
mW := 0.3 wW := 0.5
:= 1.0
w := 0.05
:= 1.0
Q := 1.14
W := 1.14
:= 3 w + w
3
3
w := 0.05
:= 3 w + w
Qk( , G, Q , W ) :=
Rd
G + ( Q + k W )
( 1 + k) ( 1 )
G( , G, Q , W )
( 1 + k) ( 1 )
Check:
Wk( , G, Q , W ) := k Qk( , G, Q , W )
G( , G, Q , W ) := wG G( , G, Q , W )
Gumbel distribution of Q:
Q( , G, Q , W ) := mQ Qk( , G, Q , W )
Gumbel distribution of W
W ( , G, Q , W ) := mW Wk( , G, Q , W ) W ( , G, Q , W ) := wW W ( , G, Q , W )
wG G( , G, Q , W ) + wQ Q( , G, Q , W ) + wW W ( , G, Q , W )
2
C - 34
E0( , G, Q , W )
E( , G, Q , W ) := E0( , G, Q , W )
E0( , G, Q , W )
wE( , G, Q , W ) :=
wE0( , G, Q , W ) + w + wE0( , G, Q , W ) w
E( , G, Q , W ) := E( , G, Q , W ) wE( , G, Q , W )
E0( , G, Q , W ) :=
Q( , G, Q , W ) Q + W ( , G, Q , W ) W
E0( , G, Q , W )
Skewness of E:
E( , G, Q , W ) :=
Parameter C:
1
2
2
E( , G, Q , W ) + 4 + E( , G, Q , W ) E( , G, Q , W ) + 4 E( , G, Q , W )
C( , G, Q , W ) :=
1
sE( , G, Q , W ) := ln 1 + C( , G, Q , W )
) x0( , G, Q, W) := E( , G, Q, W )
C( , G, Q , W )
E( , G, Q , W )
5 Resistance variables R = *R 0:
wR( wR0) :=
wR0 + w + wR0 w
Check:
C - 35
pf50( , m, wR0, G, Q , W ) :=
d( , G , Q , W)
t := 3.8 Auxiliary:
a := 0 , 0.05.. l + 0.05
( G, Q , W ) :=
wG = 0.1
k0 :=
1 Q
Check:
1 W
b = 0.6
G ( 1 ) ( 1 + k)
G ( 1 ) ( 1 + k)
G ( 1 ) ( 1 + k) + ( Q a + W k b )
0 := 3 , 3.5.. 5
mQ = 0.6
k0 = 0.75
G ( 1 ) ( 1 + k) + Q ( a Q) + W k ( b W )
a=1
am := 0 , 0.05.. la + 0.04
Check:
k=0
l = 0.31
la = 0.119
5.5
4.5
3.5
.
3
0.2
0.4
0.6
0.8
C - 36
1t( G, Q) := 3.8
2t( , m) := 3.8
1 , 1t
2 , 2t
C - 37
2( , wR0) := ln( , 1.15, wR0, 1.35, a 1.5, b 1.5) 2t( , wR0) := 3.8
2 , 2t
12 Reliability index versus m a wR: 4( m, wR0) := ln(0.4, m, wR0, 1.2, a 1.4, b 1.4) 4t(m, wR0) := 3.8
Figure 5: Reliability index
corresponding to expression
(6.10) versus parameter R
and the coefficient of variation
wR.
Parameters : k = 0
mQ = 0.6 wQ = 0.35 mW = 0.3 wW = 0.5
4 , 4t
C - 38
MATHCADsheet "LoadEfect"
Mathcad sheet "LoadEffect" is intended for investigation of combination rules
provided in EN 1990 by expressions (6.10), (6.10a) and (6.10b) considering three loads: G,
Q and W. Turkstra's rule ( 50 years extremes of a leading and annual extremes of an
accompnying action) is applied. MH, May 2003.
Design Load effect:
Ed
= () G*Gk + (Q)Q*Qk +(W
Stochastic model:
E=
* (G + Q + W)
Load effect E = * (G + Q + W) is discribed by athree parameter lognormal distribution
LN ( E, E, E), permanent load G by N(Gk, 0.1*Gk), 50 yearsimposed load Q by GUM(0.6
Qk, 0.35 Q), annual wind load W by GUM(0.3Wk, 0.5W), uncertainty by LN(1.0, 0.05).
Parameters: =( Qk+Wk)/(Gk+Qk+Wk), k=Wk/Qk, factorsG, Q,W, , Q, W.
1 Input data:
Ed := 1
:= 0.85 Q := 0.70
W := 0.60
k := 0.0
( 1 + k)
3 Probabilistic models
Permanent load G:
Normal distribution of G:
Variable load Q:
Gumbel distribution of Q:
Variable load W:
mW := 0.3 wW := 0.5
Model uncertainty
:= 1.0 w := 0.05
Q = mQ*Qk
G = Gk
wG := 0.1
mQ := 0.6
wQ := 0.35
Q := 1.14
W := 1.14
:= 3 w + w
4 Determination of the load variables G, Q and W for a given load effect Ed = Rd:
Characteristic values Gk= G, Qk and Wk determined assuming Ed=Rd. Thus G*Gk+ Q*Qk+ W*Wk=Rd
G( , G, Q , W ) :=
Qk( , G, Q , W ) :=
Ed
G +
( Q + k W )
( 1 + k) ( 1 )
G( , G, Q , W )
( 1 + k) ( 1 )
Check:
Wk( , G, Q , W ) := k Qk( , G, Q , W )
G( , G, Q , W ) := wG G( , G, Q , W )
Gumbel distribution of Q:
Q( , G, Q , W ) := mQ Qk( , G, Q , W )
Gumbel distribution of W
W ( , G, Q , W ) := mW Wk( , G, Q , W ) W ( , G, Q , W ) := wW W ( , G, Q , W )
wG G( , G, Q , W ) + wQ Q( , G, Q , W ) + wW W ( , G, Q , W )
2
C - 39
E0( , G, Q , W )
E( , G, Q , W ) := E0( , G, Q , W )
E0( , G, Q , W )
wE0( , G, Q , W ) + w + wE0( , G, Q , W ) w
E( , G, Q , W ) := E( , G, Q , W ) wE( , G, Q , W )
Q( , G, Q , W )
E0( , G, Q , W ) :=
Q + W ( , G, Q , W ) W
E0( , G, Q , W )
Skewness of E:
E( , G, Q , W ) :=
Parameter C:
1
E( , G, Q , W ) 2 + 4 + E( , G, Q , W ) E( , G, Q , W ) 2 + 4 E( , G, Q , W )
C( , G, Q , W ) :=
1
sE( , G, Q , W ) := ln 1 + C( , G, Q , W )
) x0( , G, Q, W) := E( , G, Q, W)
C( , G, Q , W )
E( , G, Q , W )
7 Theoretical value of the global factor for a given exceedance probability p of Ed:
p = P{ E > Ed }
Ed( p , , G, Q , W )
Ek( , G, Q , W )
C - 40
8 The global load factor versus ratio : limit for dominant action: k0 :=
a := if( k k0, 1 , Q)
Auxiliary quantities:
Limit value of for (6.10a) and (6.10b):
( G, Q , W ) :=
t := 3.8 Auxiliary:
a := 0 , 0.05.. l + 0.05
Turkstra's for 50 years:
1 W
a=1
k0 = 0.75
b = 0.6
G ( 1 ) ( 1 + k) + Q ( a Q) + W k ( b W )
la := a( 1.4, 1.6, 1.6)
l := ( 1.35, 1.5, 1.5)
0 := 1 , 1.1.. 1.3
mQ = 0.6
Check:
G ( 1 ) ( 1 + k)
G ( 1 ) ( 1 + k)
(
)
G 1 ( 1 + k) + ( Q a + W k b )
1 Q
c := 0 , 0.05.. la + 0.04
Check:
k=0
l = 0.31
la = 0.116
1.5
1.4
1.3
1.2
1.1
0.2
0.4
0.6
0.8
C - 41
1.5
1.4
1.2
0
0
1.1
0.2
0.4
0.6
, , , , a , b , c , l , la ,
0.8
.
G := 1.35 Q := 1.5
W := 1.5
( p , ) := p( p , , G, a Q , b W )
Parameters : k = 0
mQ = 0.6 wQ = 0.35
mW = 0.3 wW = 0.5
a=1
b = 0.6
C - 42