1. Introduction
The quantitative assessment of the safety of an existing building requires, as is known, the implementation of a numerical model that is able to describe the geometric and mechanical characteristics of the building. For a reinforced concrete (r.c.) building, according to the Italian technical codes [
1,
2], ASCE 41-17 [
3], and the Eurocodes [
4], this requires, in particular, knowledge of the construction details (e.g., the quantity and arrangement of reinforcements) and the strength of the materials (concrete and steel).
The formulation of a numerical model, in all its phases, is inevitably affected by epistemic uncertainties, e.g., those due to the possible alternatives in the choice of the representative model of the building [
5,
6,
7] or of the actions to be considered, and by random uncertainties, e.g., those due to the intrinsic variability of the mechanical parameters throughout the structure, which can therefore be described only from a probabilistic point of view.
This is the consequence of the variability of these properties within the building, which may be due to the different characteristics of the materials used, corresponding to different days during which the concrete was poured, different casting conditions, etc.
It is worth noting that, in practical investigations, modern survey technologies can ensure the exact description of the geometrical features of the investigated building, while the design documentation, if available, may provide data regarding the dimensions and positions of the steel bars inside each structural element. Of course, the acquisition of a given level of knowledge of the investigated building depends on the number of tests that are conducted to assess the mechanical properties of the materials.
For this reason, the codes usually suggest a minimum number of tests to be performed, but even if the number of such tests and the experimental technique are the same, a huge number of combinations of investigated elements and of points within the same element may be chosen, depending on the technician’s criteria and/or for logistic reasons.
This is an intrinsic source of uncertainty, even neglecting other aspects that can have effects on the results of experimental investigations, such as the environmental conditions in which the tests are performed, the instruments utilised, and several other elements depending on the chosen technique [
8,
9,
10].
Therefore, for the same knowledge level, a huge number of different experimental results may be obtained, and, consequently, for each different combination of experimental data, this may correspond to a different value of the assessed safety level of the building with respect to seismic events.
For random uncertainties, usually, codes suggest discrete knowledge levels with the aim of describing the availability of information on the existing building examined. Based on the acquired knowledge level, two possible approaches are suggested by standards: the amplification of the force demand or a reduction in the materials’ resistance. For example, FEMA-356 [
11] prescribes the introduction of a knowledge factor, dependent on the acquired knowledge level, which amplifies the force demand.
However, for the random uncertainty of mechanical parameters—which is the focus of this paper—the Italian technical codes [
1,
2], ASCE 41-17 [
3], and the Eurocodes [
4] introduce a confidence factor (CF) which reduces the mean resistance of materials as provided by on-site tests, with the CF depending on the knowledge level (KL). A larger value of CF is suggested by such codes for a lower number—and therefore a lower representativeness—of the tests carried out. In detail, the codes define three knowledge levels based on the number of acquired data: KL1 is limited knowledge, KL2 is extended knowledge, and KL3 is full knowledge. It is worth noting that different values for the confidence factor correspond to the different knowledge levels. According to the Italian technical code [
2] and the Eurocodes [
4], the confidence factor related to limited knowledge is CF
KL1 = 1.35, while the confidence factor corresponding to extended knowledge is CF
KL2 = 1.20 and the confidence factor corresponding to full knowledge is CF
KL3 = 1.00.
The possibility of taking into account the many uncertainties that arise in the procedure for the seismic safety assessment of existing buildings with a single confidence factor applied to material strengths has been, and is still, widely debated in the scientific literature [
12,
13,
14,
15,
16]. In detail, starting from the approach proposed in [
12], the present research explores the effects of uncertainties related to the choice of on-site test locations on the assessment of seismic safety indexes.
Specifically, this paper investigates, through the case study of a typical reinforced concrete structure of the 1960s built in Italy, the probability that a safety assessment based on a pre-defined knowledge level and the value consequently assigned to the confidence factor can actually lead to the technician in charge obtaining a precautionary estimate of the building’s capacity to resist seismic actions.
To this end, the actual level of seismic safety of the building—evaluated by assuming perfect knowledge (never available in practice) of the mechanical properties of all the structural elements [
12]—is compared with the safety level evaluated under the assumption that the tests were carried out, and thus that the mechanical properties were estimated, only for a small number of structural elements. In this case, as is usual, the safety index was assessed by introducing the elements’ average strengths divided by a CF whose value was fixed according to the level of knowledge corresponding to the number of tests.
For a given knowledge level, the probability distribution of the estimated safety level was evaluated in the paper by means of a Monte Carlo approach, which allows phenomena with significant uncertainties to be modelled by considering randomly generated scenarios and by evaluating the outcomes of the investigated model for each of them. In this way, the mean value of an output, its distribution, and its minimum or maximum values can be estimated. This method is used widely to evaluate the impact of risk and uncertainties in the prediction of outcomes. In the case examined here, the Monte Carlo method was applied to simulate a large set of alternative selections of the structural elements to be tested, with the same number of such elements randomly picked out, and by estimating the safety index of the building for each scenario.
Moreover, to take into account that in real cases perfect knowledge of structural characteristics can never be attained, nine alternative realisations of the reference model—with equivalent statistical properties—were considered, and therefore a total number of 9 × 103 analyses were carried out for each one of the three knowledge levels prescribed by the codes.
This procedure was intended to evaluate the effect on safety estimation uncertainty of different possible choices of structural element locations that,
ceteris paribus, can be made by the technician in charge (denoted as the analyst in the following, according to [
12]), even when the number of such tests is assumed to be constant.
In this way, the paper deals with an intrinsic uncertainty in safety evaluation, starting from the observation that no strategy can assure that the different positions of the tests chosen by different analysts—even when they perform the same number of tests—will lead to unique “best” values of the average strengths of concrete and reinforcing steel. Thus, the paper explores how the unavoidable uncertainty of these values affects the evaluation of the safety index obtained by a given analyst (i.e., what the probability is for a given analyst of underestimating or overestimating the “true”—but never really knowable—safety index). The rest of the paper is organised in four parts:
Section 2 describes the case study and the main assumptions of the performed analyses,
Section 3 describes the numerical model implemented to assess building safety levels,
Section 4 reports the statistical analyses and discusses the effects of the confidence factor values on the accuracy of the assessed vulnerability index, and
Section 5 presents the conclusions.
3. Finite Element Model of the Building
A 3D finite element model of the building was implemented with beam elements and floors acting as rigid diaphragms, assuming the columns to be fixed at the base (
Figure 1a).
In order to represent the possible variation in the mechanical properties, corresponding, for example, to different conditions of the preparation, casting, or curing of the concrete or to different batches of steel, the value of the concrete compressive strength (fc) was considered to be in the range of 9 MPa to 29 MPa, while for the steel strength (fs), the range of variability was considered to be between 300 MPa and 500 MPa.
On this basis, each floor of the building was divided into five parts (see
Figure 1b). Each part was considered representative of a different casting day or batch of steel, and, for each part and for each floor, the values of the concrete and steel strengths were randomly chosen in the mentioned ranges, as summarised in
Table 1 and
Table 2, respectively.
The aforementioned procedure aimed to fix the “real” distribution of strengths in the examined structure and, on that basis, to evaluate the “true” seismic safety index of the “real” structure, which represents the final scope of seismic vulnerability assessment.
For the “real” structure, the randomly extracted concrete strength values were characterised by a mean value of 19.6 MPa, a standard deviation equal to 6.1 MPa, and a coefficient of variation equal to 0.3, while the chosen steel strengths had a mean value of 396 MPa, a standard deviation of 63 MPa, and a coefficient of variation equal to 0.16. As it is reasonable to expect, the steel strength was characterised, therefore, by a variability lower than that of the concrete.
These values were assigned to the beams and columns of each of the different regions (
Figure 1b) of the reference model of the structure, from now on denoted as M
0 (i.e., model 0) to distinguish it from alternative realisations (
Section 4) of the same “true” model, and represent the actual mechanical strengths of these elements. The perfect knowledge of these—never attainable in practice, but here introduced as a hypothesis—allowed evaluation of the “true” level of seismic safety of the building once all the other geometrical and loading details were also known, as assumed in the following, in order to focus attention on the random variation of mechanical properties.
The dynamic analyses of the finite element model, obtained by introducing the values in
Table 1 and
Table 2, were performed considering the spectral acceleration corresponding to the Life-Safety Limit State and assuming a behaviour factor (q) equal to 2, and the seismic safety index (α) of the model was evaluated as the ratio between the peak ground acceleration capacity (PGA
C) and the peak ground acceleration demand (PGA
D).
It is worth noting that the assessed seismic safety index was evaluated by considering “exact” strength values for each of the structural elements and not a mean value constant for the whole structure.
The obtained seismic safety index for the reference model was 0.315, and the analysis highlighted that the safety of the structure was mostly conditioned by the combined axial compression and bending moment capacity. Of course, the structural element capacity that conditions the safety of buildings in other cases studies may be different from the one considered here. However, the results discussed in
Section 4 highlight some issues that are also expected to be observed in other possible structural typologies and/or different mechanisms of limit state attainments.
The subsequent step of the numerical investigation described here aimed at evaluating the impact that, during the acquisition of a fixed knowledge level, the alternative locations of the in situ tests selected by different analysts can have on the estimated safety level, all other conditions being equal. Different locations may in fact correspond to different experimentally acquired strengths and, consequently, to different average values of the strengths adopted by the technician in charge of the structural analyses on the basis of the available data. The obtained mean values and the corresponding confidence factors were then utilised to assess the seismic safety indexes.
To compare the evaluated indexes, it was assumed that the attainment of the Life-Safety Limit State is always conditioned by the combined axial compression and bending moment capacity of the columns. Without affecting the more general validity of the proposed approach, which will be further explored in the following part of this research, such an assumption allows, by means of the interpolating surfaces introduced in
Section 4, a reduction in the computational effort required.
4. Effects of Uncertainties Related to Test Locations
According to the Italian technical codes [
1,
2], the seismic vulnerability assessment of an existing structure always requires the execution of in situ tests on the structural elements. Their number depends on the choice of the knowledge level, which in turn is related to the confidence factor (CF) (see
Section 1) introduced to pass from the average values of test results to the design values of materials’ strengths. All other factors being equal (the geometry of the building, the quantity and arrangement of the reinforcements, etc.), as assumed in this paper, the design strength of concrete and steel, on the other hand, affect the safety level of the building provided by the structural analysis.
However, even when the number of tests can be assumed as fixed (e.g., according to the indications of technical codes), the technician in charge may make alternative choices regarding the position of the tests. Therefore, for the same knowledge level, different analysts may estimate different mechanical properties for the same building.
As a consequence, even if the numerical model developed to assess the response of the investigated structure and the confidence factor are the same with respect to any other characteristic of the building, different analysts may estimate different values of the seismic safety index because of the variation in the mean values of the concrete and steel strengths that each analyst assumed to be representative of the existing building on the basis of the available in situ tests. Mean values of concrete and steel strength can therefore be viewed as samples for a statistical analysis of the seismic safety index, as discussed below.
The procedure here proposed is based on a Monte Carlo simulation to select the structural elements whose strengths are assumed to be measured in situ, with the aim of quantifying the uncertainty in the evaluated safety index as the location of the tests varies while keeping their number fixed.
In detail, the probability distribution of the safety indexes evaluated by 100 different analysts was obtained, i.e., considering the mean values of concrete and steel strengths corresponding to 100 different combinations of test locations.
To reduce the computational effort, however, for each given knowledge level, only 15 numerical FEM models were analysed, each one corresponding to a different combination of the mean strengths (
and
) virtually obtained by different analysts. As is usual in engineering practice, the obtained mean values were considered to be constant throughout the whole structure during the assessment of the safety index, which is different from the procedure (
Section 3) for the evaluation of the actual index, where the knowledge of the local strength of each element was assumed and introduced in the model. The evaluated safety indexes for each
,
pair are listed in
Table 3 for knowledge level KL2 and a confidence factor (CF
KL2) equal to 1.2 [
2]. The safety indexes (
) corresponding to other combinations of
and
may be obtained by fitting the obtained value by means of a polynomial law (see
Figure 2). Similarly,
Table 4 shows the safety indexes (
) corresponding to knowledge level KL1 and CF
KL1 = 1.35.
At first the analyses with knowledge level KL2 are described below, assuming that no other information (e.g., results of tests during the building process) was available on the strength of materials. Under this hypothesis and with reference to the columns, i.e., the elements by which the safety level is conditioned in this case study (see
Section 2), the Italian technical codes [
1,
2] suggest extended in situ testing that involves the testing of two samples of steel bars per storey and two concrete cores for each 300 m
2 of floor, and then six cores per storey for the investigated building (about 800 m
2 plan). For the five-floor building considered here, therefore, the strengths (
and
) to be considered in the safety evaluation were obtained for each analyst as the average of the strengths of 30 cores and 10 steel bars extracted by columns.
A Monte Carlo approach was adopted to define 100 different cases. For each case, the procedure simulated the extraction of six values of concrete strength and two values of steel strength for each floor, corresponding to a random subset of the columns, and, based on such values, the mean values of the concrete compressive strength (
) and steel strength (
) were evaluated. The obtained values are plotted in
Figure 3.
By means of the interpolating surface in
Figure 2, the seismic safety index was then evaluated for each combination of strength values in
Figure 3, and their cumulative distribution is shown in
Figure 4. Moreover, in
Figure 4, the cumulative distribution of the Gaussian curve which fits the evaluated safety indexes is plotted.
The red broken line in
Figure 4 corresponds to the safety index evaluated by means of the reference model (M
0) (
Section 2), which assumed perfect knowledge (never attainable in practice) of the effective strength of each structural element. The safety index of the reference model can therefore be assumed as an ideal target, whose comparison with the indexes evaluated by different analysts may allow estimation of the reliability of the safety level attainable when a finite (and usually small) number of in situ tests are available, as in real cases.
Figure 4 highlights that for the randomly selected locations of in situ tests corresponding to the average values of
and
in
Figure 3, the probability of overestimating the seismic safety is negligible once the confidence factor CF
KL2 = 1.2 is adopted in the analysis, as suggested by the Eurocodes [
4] and the Italian codes [
2] for the level of knowledge KL2.
To check with a wider set of analysts the probability of obtaining an unsafe evaluation of the safety index, the aforementioned procedure was repeated another nine times for different sets of 100 randomly selected values of concrete and steel strength, leading to another nine cumulative distribution functions of the safety index, which are shown in
Figure 5 together with the curve in
Figure 4. A total number of 10
3 combinations were therefore analysed to obtain the data plotted in
Figure 5.
This confirms that for the reference model (M
0) corresponding to the concrete and steel resistance in
Table 1 and
Table 2, this CF
KL2 value is well calibrated because it effectively guarantees, in almost all cases, a slightly precautionary estimation of the structural capacity.
Due to the inhomogeneity of the mechanical properties in
Table 1 and
Table 2, however, the capacity of the reference model (M
0) considered up to this point—and therefore its safety risk (α
0 = 0.315)—was conditioned by the effective distribution of such characteristics in the different regions of each storey, which is never completely knowable in practice.
To evaluate the impact of such a distribution on the capacity of the model assumed as a reference to describe the real structure, therefore, eight alternative models (M1–M8) were considered to be representative of eight different realisations of the building described so far. These were obtained by randomly re-distributing the steel and the concrete properties of the elements between the different 25 regions of the building and thus obtaining different combinations of , without changing the parameters (i.e., the mean values and standard deviations) of their values. The procedure performed for the M0 model was also carried out for each of the eight new models, leading to 10 cumulative distribution functions for each model, i.e., a total number of 9 × 103 analyses.
For each model,
Figure 6 compares the cumulative distribution functions of the safety indexes evaluated by 10 different sets of analysts with the actual safety index of the model.
Figure 7a shows, for each model and for each set of analysts, the probability of underestimation (
p*) of the actual safety index (blue circles) and the average μ
p* (red crosses)—also reported in
Table 5—as a function of the actual safety index of each model. Although for one of the models (model M
7) the mean probability of overestimating the actual safety index was about 50%, the overall mean value of the probability of underestimation (
p*) for the nine models was 92% (
Table 5), thus confirming that the confidence factor CF
KL2 allows a cautionary evaluation of the seismic safety index.
Based on the statistical distribution of the safety indexes of the nine models (M
0–M
8) and denoting their average values and standard deviations μ
α and σ
α, respectively,
Figure 7a also reports, as a vertical green line, the average value (μ
α = 0.313), and it can be verified that, in this case, the probability of underestimation of the safety index by the analysts is almost equal to 100%.
Figure 7a also shows with a green broken line the value of the safety index μ
α − σ
α = 0.299 with an 84% probability of exceedance, and in this case the probability of underestimation of the safety index by the analysts is equal to 87%.
The same results can be shown by plotting the probability (
pexc) of exceeding the actual safety index, as reported in
Figure 7b.
It is worth noting that for a major part of the analysed models (i.e., seven models out of the total number of nine), almost all the analysts underestimated the building’s actual structural capacity with an average probability of at least 95% and that, considering all of the nine models, the mean value of the probability of underestimation was equal to 92%. Thus, adopting a confidence factor CFKL2 = 1.2, the probability of overestimating the real structural capacity of the sample case discussed here is limited to about 8%. This value can be considered acceptable when compared with other sources of uncertainties which are usually involved in structural capacity assessment, such as those related to the choice of software, the method adopted for the analysis (linear or non-linear), and those strictly connected to the experimental tests.
To investigate the influence of the number of in situ tests, the procedure was repeated by assuming that three concrete cores and one steel bar were extracted at each of the five storeys, with a total number of 15 values of concrete strength and 5 values of steel strength randomly selected in a subset of the columns for each of the 100 simulated samples.
In this case, corresponding to the knowledge level KL1 of the Italian technical code, a confidence factor CF
KL1 = 1.35 was suggested [
2]. As shown in
Figure 8, which reports the cumulative distribution functions of the seismic safety indexes obtained by means of an interpolating surface based on
Table 4 (confidence factor CF
KL1 = 1.35), the probability of an unsafe evaluation becomes negligible, as also shown in
Figure 9, which was obtained by means of the same procedure adopted to evaluate
Figure 7.
This confirms that for the lowest level of knowledge allowed by the codes, the value of CFKL1 = 1.35 is well calibrated in the case study described here.
The above-described procedure was also carried out to investigate the case of the full-knowledge level, i.e., KL3, when the confidence factor is equal to 1.0 [
2], and, in this case, the seismic safety index of the sample building discussed here was always overestimated.
Although this aspect deserves to be further explored in the future by examining different structures, what seems to emerge clearly from these circumstances is that the choice of the confidence factor cannot always and only be related to the number of available tests and that it should instead explicitly take into account their higher or lower homogeneity. Accordingly, it should eventually be assumed to be equal to 1 only when the dispersion of experimental results is particularly low, unlike the case reported here.
For the building considered here, it was in fact observed that, even if the nine alternative (and statistically equivalent) models M
0–M
8 are characterised by the same mean values of concrete and steel strength, i.e.,
fc* = 19.64 MPa and
fs* = 396 MPa, their safety indexes vary with respect to their mean values (
μα) in an interval ±1.46 times their standard deviations
σα. Moreover, the safety indexes of models M
0–M
8 (probably because they are conditioned by the structural elements with the lowest resistance) were generally lower than the “ideal” safety index (
α* = 0.337), corresponding to the hypothesis of absolute homogeneity of the mechanical characteristics of concrete and steel, i.e., equal to f
c* and f
s* for each structural element, similarly to the results of [
12]. As a consequence, assuming the average values of experimental tests as design strengths, as when CF = 1 is assigned, can lead to the overestimation of safety indexes in cases of significant material inhomogeneity, as in those of models M
0–M
8.
This suggests that, in such cases, the confidence factor should be precautionarily assumed to be greater than 1, independently of the number of in situ tests, in order to take into account that the structural safety of a building can be highly conditioned by the structural elements with the lowest capacities. The suggestion of adopting in such cases a confidence factor greater that 1 is in accordance with results in [
12].
Future research will therefore be devoted to applying the above-described procedure to different probability distributions of mechanical properties, with the aim of exploring how and to what extent the calibration of the confidence factor should include the dispersion of experimental results, an issue that is only qualitatively considered in the Italian codes [
2].
Extending the results discussed here, which were obtained by adopting the linear analysis with behaviour factor, the future steps of the research will also consider the effect of the introduction of non-linear analyses (for example, push-over analyses, etc.) for the assessment of the seismic safety indexes of buildings with mechanical parameter distributions similar to those examined here, with the aim of verifying the influence of their mean values and their variations on safety indexes.
5. Conclusions
In this paper, the case study of a framed reinforced concrete building designed to resist only gravitational and static loads, representative of a widespread typology of structures built in the 1960s in Italy, is examined. Indeed, buildings of this kind usually stand in seismic-prone zones, according to the current codes, and this requires improvements to their structural and building performances to guarantee acceptable levels of safety and habitability conditions. The actuality of this topic is also confirmed by the attention devoted by building codes to the analysis of existing buildings.
These codes prescribe, among other procedures, the execution of experimental tests, the numbers of which depend on the desired knowledge level.
Although this topic is still debated in the scientific literature [
12,
13,
14,
15,
16], the Italian technical codes [
1,
2] and the Eurocodes [
4] prescribe the introduction of a confidence factor, based on the acquired knowledge level, in the assessment of seismic safety indexes, which serves to reduce the mean values of experimentally evaluated mechanical strengths of structural materials.
It is worth noting that the positions in which the tests are performed depend on several factors, and, clearly, different analysts may choose different locations; consequently, they may assess different values of the mechanical properties and obtain different safety indexes.
In the present paper, starting from the approach in [
12], the case study of a residential r.c. building built in the 1960s in Italy is discussed with the aim of investigating the implications of the different choices that the technician in charge (the “analyst”) can make regarding the structural elements to be tested for obtaining a prescribed level of knowledge.
The case study of an r.c. building was examined, assuming, initially, perfect knowledge of the material mechanical properties of each structural element and, on this basis, evaluating the “actual” safety index.
Then, numerical analyses were performed by considering the three knowledge levels prescribed by the building technical codes and, for each level, numerically simulating the acquisition of experimental strengths.
In detail, a Monte Carlo approach was adopted to consider the effects of the different locations of the tests. A numerical simulation of 103 possible combinations of concrete and steel mechanical properties, as virtually assessed for each knowledge level suggested by the codes, was carried out, and the corresponding seismic safety indexes were evaluated by introducing the confidence factor suggested by technical codes for each knowledge level.
To evaluate the capacity of the model assumed as a reference to describe the real structure, the same procedure was applied to nine different realisations of the same building, and thus a total number of 9x103 combinations were analysed for each knowledge level.
The analyses performed for the sample case show that the confidence factors corresponding to KL1 and KL2 [
2] could effectively provide a precautionary assessment of the structural capacity of the building. With the level of knowledge KL2, in fact, for the analysed models (representative of nine alternative realisations of the reference model with equivalent statistical properties), almost all the analysts underestimated the actual building’s structural capacity with an average probability of 92%. Thus, adopting a confidence factor CF
KL2 = 1.2, the probability of overestimating the real structural capacity of the sample case discussed here was limited to about 8%, while for the level of knowledge KL1 (CF
KL1 = 1.35), the probability of an unsafe evaluation became negligible.
For the full-knowledge level, the safety index was instead overestimated when a value of CFKL3 = 1 was assumed, thus suggesting that the calibration of this parameter should not only depend on the number of available in situ tests but must consider their dispersion.
As is widespread in engineering practice, in the present paper a linear numerical analysis was performed. However, further studies must be performed to extend the findings to the assessment of safety indexes by means of non-linear approaches.