Creating Fragility Functions for Performance-Based Earthquake Engineering:
Summary (4 min read)
BACKGROUND AND OBJECTIVES
- This paper summarizes such a standard developed for ATC-58.
- See Porter et al. (2006) for more detail, examples, commentary, and alternative approaches.
- Method A is not applicable when one knows the maximum EDP to which each specimen was subjected, but not the value of EDP at which specimens actually failed.
- The methods proposed here are no substitute for understanding the processes that lead to damage, but are intended to help practitioners and scholars create fragility functions from damage data.
DOCUMENTATION REQUIREMENTS
- Four requirements are proposed for documenting fragility functions: 1. Description of specimens.
- Indicate whether EDP is the value at which damage occurred (Method A data) or the maximum to which each specimen was subjected (Methods B, C, and U).
- Define damage measures (DMs) quantitatively in terms of repairs required.
DAMAGE STATE PROBABILITY
- Fdm edp denotes the fragility function for damage state dm, defined as the probability that the component reaches or exceeds damage state dm, given a particular EDP value (Equation 1), and idealized by a lognormal distribution (Equation 2): Fdm edp P DM dm EDP = edp 1.
- Fdm edp = ln edp/xm 2 where denotes the standard normal cumulative distribution function (e.g., normsdist in Microsoft Excel), xm denotes the median value of the distribution, and denotes the logarithmic standard deviation.
- The authors use the lognormal because it fits a variety of structural component failure data well (e.g., Beck et al.
- Both xm and are established for each component type and damage state using methods presented later.
- The probability that the component is in damage state dm, given EDP=edp, is given by.
METHOD A, ACTUAL EDP: ALL SPECIMENS FAILED AT OBSERVED EDP
- These are the most informative data for creating fragility functions.
- They are most common where DM can be associated with a point on the observed force-deformation behavior of a component, such as a yield point.
- Alternatively, specimens are subjected to increasing levels of EDP.
- The test is interrupted after each level of EDP is imposed, and the specimen examined for damage.
M number of specimens tested to failure i index of specimens, i 1,2 , . . .M ri EDP at which damage was observed to occur in specimen i.
- One tests the resulting fragility function using the Lilliefors goodness-of-fit test (presented below).
- If it passes at the 5% significance level, the fragility function is acceptable.
- Example 1. Aslani (2005) provides a table of peak transient drift ratios at which 43 specimens of pre-1976 reinforced concrete slab-column connections experienced cracking of no more than 0.3 mm width, repaired by applying a surface coating.
- The data are repeated in Table 2 with original specimen numbers.
- The lognormal distribution with these parameters passes the Lilliefors goodness-of-fit test at the 5% significance level.
METHOD B, BOUNDING EDP: SOME SPECIMENS FAILED, PEAK EDP KNOWN
- Here, the data include the maximum EDP to which each of M specimens was subjected, and knowledge of whether the specimen exceeded the damage state of interest.
- Data must not be biased by damage state, i.e., specimens must not be selected because they experienced damage.
- The data are grouped into bins by ranges of EDP, where each bin has approximately the same number of specimens in it.
- These serve as independent data points of failure probability and EDP.
- The following approach converts Equation 2 to a linear regression problem by taking the inverse Gaussian cumulative distribution function of each side and fitting a line ŷ=sx+c to the data (e.g., see “probability paper” in Ang and Tang 1975).
N number of EDP bins
- Porter et al. (2006) presents an alternative approach using a leastsquares fit to the binary failure data, i.e., to the pairs of EDP and a binary (0,1) failure indicator.
- The alternative approach avoids errors associated with bin-average EDPs.
- Consider the damage statistics in Figure 2, which depicts motor control centers (MCCs) observed after various earthquakes in 45 facilities.
- Crosshatched boxes represent MCCs that experienced noticeable earthquake effect such as shifting but that remained operable.
- For each bin, the values of xj− x̄ and yj− ȳ are calculated as shown.
METHOD C, CAPABLE EDP: NO SPECIMENS FAILED, EDPS ARE KNOWN
- It addresses the best case for this type of data, i.e., many specimens, none of which had apparent distress, and several of which were subjected to EDP near the maximum value.
- The specimens in this bin without apparent distress are assigned 0% subjective failure probability, 10% for specimens with distress not suggestive of imminent failure, and 50% for specimens with distress suggestive of imminent failure.
- ANCO Engineers, Inc. (1983) performed shake-table tests on ceiling systems with various lateral restraints.
- Peak diaphragm acceleration (PDA) from nine of these tests is recorded in Table 5.
- Failure required replacement of damaged grid and tiles.
METHOD D, DERIVED FRAGILITY FUNCTIONS
- The capacity of some components can be calculated by modeling the component as a structural system, and determining the EDP (e.g., acceleration or shear deformation) that would cause the system to reach dm.
- Other components may be amenable to fault tree analysis; e.g., see Vesely et al. (1981).
- Let r denote the calculated capacity of the component to resist damage state dm, including consideration of any anchorage or bracing.
METHOD E, EXPERT OPINION
- There are several methods for eliciting expert opinion, from ad hoc to structured processes involving multiple experts, self-judgment of expertise, and iteration to examine major discrepancies between experts.
- The method (introduced for the first time here) employs Spetzler and von Holstein (1972) for probability encoding and Dalkey et al. (1970) for expert qualification, with some useful simplifications.
- To use Method E, select experts with professional experience in the design or postearthquake observation of the component.
- Representative images should be offered to the experts and recorded.
- If an expert refuses to provide estimates or limits them to certain conditions, either narrow the component definition accordingly and iterate, or ignore that expert’s response and analyze the remaining ones.
N number of experts providing judgment about a value i index of experts, i 1,2 , . . .N
- If the results of the survey produce 0.4, and this low value of cannot be justified, use the judged xl to anchor the fragility function, apply =0.4, and calculate the resulting value of xm.
- Kennedy and Short (1994) show that by establishing the EDP at which the component has 10% failure probability, the overall reliability of the component is insensitive to , hence the value of directly encoding experts’ judgment of this value in particular.
- Stone cladding on the exterior of retail buildings may fall in earthquakes.
- Consider 2-in. x 6-in. x 1-3/16-in. stone veneer adhered to a concrete masonry unit substrate with thin-bed mortar (liquid latex mixed with Portland cement, 100% coverage).
- Create a fragility function for the probability that any given stone would fall from the building (posing a life-safety threat) and require replacement, as a function of the peak transient drift ratio of the story on which the stone is applied.
METHOD U, UPDATING A FRAGILITY FUNCTION WITH NEW DATA
- Here, the data are a pre-existing fragility function and M specimens with known damage state and maximum EDP.
- It is not necessary that any of the specimens failed.
- The method uses Bayes’ Theorem (e.g., Ang and Tang 1975) to revise xm and of an existing fragility function with new observations of M specimens whose EDP and damage state have been observed.
- For those familiar with Bayesian updating, the prior probability distribution of xm is taken as lognormal with median equal to the xm value in the pre-existing fragility function, and logarithmic standard deviation taken as 0.707 of the pre-existing fragility function, consistent with a compound lognormal fragility function and r= u=0.707 .
- Their joint distribution is approximated by five discrete points (xmj, j), each with probability-like weight wj (where j=1,2 , . . .5).
ASSESSING FRAGILITY FUNCTION QUALITY
- The previous section provided mathematical procedures for developing fragility functions.
- Issues associated with the quality of those fragility functions are now addressed, particularly the treatment of competing EDPs, goodness-of-fit testing, dealing with fragility functions that cross, and how to assign an overall quality level to a fragility function.
CONSIDERING COMPETING EDPS
- One may be uncertain which is the best EDP to use.
- In such a case, create fragility functions for each alternative and choose the fragility function with the lowest .
- See Porter et al. (2006) for choosing between EDPs with differing COV.
GOODNESS OF FIT
- A goodness-of-fit test checks that an assumed distribution adequately fits the data.
- It is a special case of the Kolmogorov-Smirnov (K-S) test, applicable when the parameters of the distribution are estimated from the same data as are being compared with the distribution, as is the case here.
FRAGILITY FUNCTIONS THAT CROSS
- Some components have two or more fragility functions.
- Any two lognormal fragility functions i and j with medians xmj xmi and logarithmic standard deviations i j cross if: edp exp jln xmi − iln xmj j − i : i j edp exp jln xmi − iln xmj j − i : i j 25.
- This produces a negative probability of being in damage state i under Equation 3b.
- Figure 6a illustrates the point: F2 has a higher than F1, and F3 has a lower than F2.
- Two methods are proposed to deal with the problem.
ASSIGNING A SINGLE QUALITY LEVEL TO A FRAGILITY FUNCTION
- Fragility functions come from data with varying quantity and quality.
- It is based solely on the authors’ judgment.
- The analyst should report the quality of fragility functions used with any loss estimate.
CONCLUSIONS
- Six methods for creating fragility functions were presented, including three new ones: one for dealing with cases where no failure has been observed, another for situations where one must rely on expert opinion, and a third for updating an existing fragility function with new damage observations.
- The procedures are under consideration as a standard for ATC-58, a technology-transfer project by the Applied Technology Council to bring PEER’s performance-based earthquake engineering methodology to practice.
- The procedures are intended for engineering professionals who will eventually use PBEE.
- Little unfamiliar math is involved, and no calculus.
- A larger document, Porter et al. (2006), presents these procedures with more commentary, some alternative approaches, and more sample problems.
Did you find this useful? Give us your feedback
Citations
896 citations
Cites background or methods from "Creating Fragility Functions for Pe..."
...…of damage, static structural analyses, or judgment (e.g., Kennedy and Ravindra 1984; Kim and Shinozuka 2004; Calvi et al. 2006; Villaverde 2007; Porter et al. 2007; Shafei et al. 2011), but here the focus is on so-called analytical fragility functions developed from dynamic structural analysis....
[...]
...…lognormally distributed; this is a common assumption that has been confirmed as reasonable in a number of cases (e.g., Ibarra and Krawinkler 2005, Porter et al. 2007, Bradley and Dhakal 2008, Ghafory-Ashtiany et al. 2011, Eads et al. 2013); however, it is not required and alternate assumptions…...
[...]
...This fragility-fitting approach has been used widely and is denoted “Method A” by Porter et al. (2007)....
[...]
...Another alternative is fitting the function using “Method B” of Porter et al. (2007), which transforms the observed fractions of collapse so that linear regression can be used to estimate the fragility function parameters....
[...]
776 citations
257 citations
241 citations
209 citations
Cites methods from "Creating Fragility Functions for Pe..."
...This method was first applied in earthquake engineering studies for seismic damage assessment (i.e., Porter et al. 2007) and has now been adapted to tsunami damage probability estimation....
[...]
References
3,923 citations
2,679 citations
"Creating Fragility Functions for Pe..." refers methods in this paper
...The following approach converts Equation 2 to a linear regression problem by taking the inverse Gaussian cumulative distribution function of each side and fitting a line ŷ=sx+c to the data (e.g., see “probability paper” in Ang and Tang 1975)....
[...]
...From the basic definitions of xm and e.g., Ang and Tang 1975 , xm = exp 1M i=1 M ln ri = 1M − 1 i=1 M ln ri/xm 2 4 One tests the resulting fragility function using the Lilliefors goodness-of-fit test (presented below)....
[...]
[...]
1,266 citations
1,122 citations
768 citations
"Creating Fragility Functions for Pe..." refers methods in this paper
...The following approach converts Equation 2 to a linear regression problem by taking the inverse Gaussian cumulative distribution function of each side and fitting a line ŷ=sx+c to the data (e.g., see “probability paper” in Ang and Tang 1975)....
[...]
..., Ang and Tang 1975) to revise xm and of an existing fragility function with new observations of M specimens whose EDP and damage state have been observed. Some explanation may be useful to readers unfamiliar with Bayesian updating. It is recognized here that xm and are themselves uncertain, and can be assigned probability distributions. The distributions are revised based on how likely it is that the observed damage would have occurred for various possible values of xm and . For those familiar with Bayesian updating, the prior probability distribution of xm is taken as lognormal with median equal to the xm value in the pre-existing fragility function, and logarithmic standard deviation taken as 0.707 of the pre-existing fragility function, consistent with a compound lognormal fragility function and r= u=0.707 . The prior of is taken as normal with expected value equal to the of the pre-existing fragility function, and coefficient of variation (COV) of 0.21. This COV is selected because it provides for 98% probability that is within the bounds of 0.5 and 1.5 times the prior , which agrees with the observed range for of 0.2 to 0.6. The distributions of xm and are assumed to be independent. Their joint distribution is approximated by five discrete points (xmj, j), each with probability-like weight wj (where j=1,2 , . . .5). Using a method described in Julier (2002), the values of xmj, j, and wj are chosen so that the first five moments of the discrete joint distribution match those of the continuous joint distribution....
[...]
...From the basic definitions of xm and e.g., Ang and Tang 1975 , xm = exp 1M i=1 M ln ri = 1M − 1 i=1 M ln ri/xm 2 4 One tests the resulting fragility function using the Lilliefors goodness-of-fit test (presented below)....
[...]
Related Papers (5)
Frequently Asked Questions (13)
Q2. what is the fragility function for a stone?
Create a fragility function for the probability that any given stone would fall from the building (posing a life-safety threat) and require replacement, as a function of the peak transient drift ratio of the story on which the stone is applied.
Q3. What is the value of encoding experts’ judgment of the value of a component?
Kennedy and Short (1994) show that by establishing the EDP at which the component has 10% failure probability, the overall reliability of the component is insensitive to , hence the value of directly encoding experts’ judgment of this value in particular.
Q4. how to create a fragility function from method c data?
To create a fragility function from Method-C data, letri EDP experienced by specimen i i=1,2 , . . .M rmax=maxi ri rd minimum EDP experienced by any specimen with distress ra the smaller of rd and 0.7rmax MA number of specimens without apparent distress and with ri ra MB number of specimens at any ri with distress not suggestive of imminent failure MC number of specimens at any ri with distress suggestive of imminent failure rm=rmax if MB+MC 0=0.5· rmax+ra otherwise S subjective failure probability at rmS = 0.5MC + 0.1MB / MA + MB + MC 14Use Table 4 to determine Fdm rm and Equation 15 to determine and xm.= 0.4z = −1 Fdm rmxm = rmexp − z 15Example.
Q5. What is the prior probability distribution of xm?
For those familiar with Bayesian updating, the prior probability distribution of xm is taken as lognormal with median equal to the xm value in the pre-existing fragility function, and logarithmic standard deviation taken as 0.707 of the pre-existing fragility function, consistent with a compound lognormal fragility function and r= u=0.707 .
Q6. What is the subjective failure probability for a bin of specimens without distress?
The specimens in this bin without apparent distress are assigned 0% subjective failure probability, 10% for specimens with distress not suggestive of imminent failure, and 50% for specimens with distress suggestive of imminent failure.
Q7. What is the only possible expression for calculating fragility functions?
No calculus is required, and the only possibly unfamiliar expression is the Gaussian distribution, typically available in spreadsheet software.
Q8. How many times does an expert say that the EDP is a given?
0.4xm = 1.67xl 18 Regarding Equation 18, it is common for experts to express overconfidence in an uncertain variable, such as the EDP at which damage will occur.
Q9. What is the purpose of this paper?
This paper introduces a set of procedures for creating fragility functions from various kinds of data: (A) actual EDP at which each specimen failed; (B) bounding EDP, in which some specimens failed and one knows the EDP to which each specimen was subjected; (C) capable EDP, where specimen EDPs are known but no specimens failed; (D) derived, where fragility functions are produced analytically; (E) expert opinion; and (U) updating, in which one improves an existing fragility function using new observations.
Q10. What is the probability that the component is in damage state?
The probability that the component is in damage state dm, given EDP=edp, is given by=Fdm edp − Fdm+1 edp 1 dm N=Fdm edp dm = N 3where N denotes the number of possible damage states for the component, in addition to the undamaged state, and dm=0 denotes the undamaged state.
Q11. How do you calculate the fragility function of a ceiling?
The capacity of some components can be calculated by modeling the component as a structural system, and determining the EDP (e.g., acceleration or shear deformation) that would cause the system to reach dm.
Q12. What is the simplest way to determine the fragility of a stone?
The method uses Bayes’ Theorem (e.g., Ang and Tang 1975) to revise xm and of an existing fragility function with new observations of M specimens whose EDP and damage state have been observed.
Q13. What is the method for eliciting expert opinion on uncertain quantities?
To properly elicit expert opinion on uncertain quantities requires attention to clear definitions, biases, assumptions, and expert qualifications.