12th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12 Vancouver, Canada, July 12-15, 2015 1Towards Quantifying The Effect Of Aftershocks In Seismic Risk Assessment Fatemeh Jalayer Associate Professor, Dep. of Structures for Engineering and Architecture, University of Naples Federico II, Naples, Italy Hossein Ebrahimian Post-doctoral Researcher, Dep. of Structures for Engineering and Architecture, University of Naples Federico II, Naples, Italy Gaetano Manfredi Professor, Dep. of Structures for Engineering and Architecture, University of Naples Federico II, Naples, Italy ABSTRACT: The state of the art in seismic design and assessment of structures implicitly relies on structural redundancy in order to deal with the effect of the triggered sequence of aftershock events on a building. Calculating the time-dependent limit state exceedance probability for a structure considering both the main event and the triggered sequence of aftershocks is complicated both by the time-dependent rate of aftershock occurrence and also the cumulative damage caused by the sequence of events. Taking advantage of a methodology developed in previous works by the authors for post-mainshock risk assessment, the limit state probability due to a sequence of a mainshock and the triggered aftershocks is calculated herein. Moreover, a closed- and analytic-form approximation to the post-mainshock limit state probability given the mainshock magnitude is derived. This closed- and analytic formulation facilitates the estimation of the limit state probability by employing the standard tools in risk assessment such as the fragility curve for the intact structure. Applying the proposed methodology for a reinforced concrete (RC) concrete moment-resisting frame with infills subjected to the main event and the triggered sequence emphasizes the importance of taking into account the cumulative damage caused by the triggered sequence for the case-study model. Moreover, it is demonstrated that the proposed analytical approximation based on the fragility curve for the intact structure leads to a surprisingly close agreement with the best-estimate results obtained by considering the time- and event-dependent degradation in the structure. 1. INTRODUCTION Considerable efforts have been made in the past decade in order to take into account additional damage potential due to aftershocks in seismic vulnerability assessment. This issue is addressed by subjecting the structure to various mainshock-aftershock (MS-AS) sequences. These sequences are constructed by combining a MS record with another one or by repeating it to account for aftershock record (see e.g., Luco et al. 2004, Yeo and Cornell 2009, Goda and Taylor 2012, Raghunandan et al. 2014), or by using real MS-AS duos (e.g., Goda and Taylor 2012, Goda and Salami 2014). These procedures provide a sound estimation of post-earthquake functionality. On a different note, Jalayer et al. (2011a, b) proposed a methodology for time-dependent risk assessment in the post-mainshock environment by calculating the probability of exceeding prescribed limit states. Conditioned on the occurrence of a given MS event, this method aimed to consider the progressive damage 12th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12 Vancouver, Canada, July 12-15, 2015 2induced by the occurrence of a sequence of aftershocks (represented by a sequence of strong-motion events). Ebrahimian et al. (2014b) refined and implemented this methodology within a performance-based and adaptive framework in order to perform operational aftershock risk forecasting in the post-mainshock environment. This study strives to quantify and to compare the limit state exceedance probability due to a MS and the triggered sequence of aftershocks (MS+AS sequence) with respect to the risk calculated by only considering the mainshock event (MS). The procedure outlined herein for calculating the limit state probabilities considers explicitly both the time-dependent rate of occurrence for aftershocks and the cumulative damage caused by the seismic sequence (MS+ AS sequence). This non-trivial task is facilitated by deriving a simple closed-form and analytical approximation to the post-mainshock limit state exceedance probability. Based on the derived analytical formulation, a simple and approximate solution is presented later on by employing only the fragility curve for the intact structure. As a numerical example, time-dependent risk related to (MS+AS sequence) is calculated through both the outlined “best-estimate” and “approximate” procedures. Risk is evaluated in terms of the probability of exceeding the Near Collapse (NC) limit state for a typical RC frame building with infill panels located in L’Aquila, central Italy. 2. METHODOLOGY In this section, the “best-estimate” procedure for calculating the limit state exceedance probability due to a MS and the triggered aftershock sequence (MS+AS sequence) is presented. Risk is calculated herein in terms of the first-excursion probability of a desired (discrete) limit state. Since the triggered aftershocks following a given MS usually take place in a span in the order of weeks or months, it is assumed that no additional main event will take place in the considered time interval. 2.1. Time-dependent Risk Formulation Let =P(LS|I1) be the probability of the first-excursion of a desired limit state LS in time interval [0,t] given information I1. I1: The structure may be subjected to a MS event (neglecting the probability of having more than one MS event of interest) with magnitude Mms so that Ml≤Mms≤Mu (Ml and Mu are the site-specific lower- and upper-bound magnitudes) followed by a sequence of triggered aftershocks (AS sequence) with magnitudes Mas in the range of Ml,as≤Mas≤Mu (Ml,as is the lower bound magnitude for aftershocks). The dependence on I1 is dropped hereafter for the sake of brevity and readability. Nevertheless, unless otherwise specified, the probability terms derived in this work are all conditioned on I1. The probability can be further broken down as follows: LS is exceeded after the MS (ms), or it is exceeded during the AS sequence that is triggered by the MS given that it is not exceeded due to the MS (as): 1ms as ms (1) The probability term ms can be interpreted as the limit state probability associated with the MS event. Adopting the first-mode spectral acceleration Sa(T) as the ground motion intensity measure and assuming that the limit state exceedance can be described by a homogenous Poisson process with rate LS, 0 d ( ) ddmsms msLS ms msmssasasa sasa (2) where sams is the Sa(T) associated with the MS event; 0=P(LS|sams) is the structural fragility for the MS event (also known as the fragility of the intact structure); and ms is the site-specific hazard defined as the mean rate of exceeding Sa(T) (note that the rates LS, ms and time t should have consistent units; i.e., the product of the rate and time is dimensionless). Accordingly, the probability of exceeding the limit state ms in time interval [0,t] can be calculated as: 12th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12 Vancouver, Canada, July 12-15, 2015 3 1 expms LSt (3) The probability as is the limit state first-excursion probability associated with the AS sequence following the MS event given that the MS has not caused the limit state excursion. This term can be expanded based on Total Probability Theorem (Benjamin and Cornell 1970) over all possible MS wave forms xg and magnitudes m: ( , ),1| , | ( )d d1 | ,x mggmg g gas asNi ggm as iiP LS m x p x m p m x mN P LS m x (4) where P(LSas|m,xg) is the probability of LS first-excursion for an arbitrary AS event belonging to the AS sequence, given that the MS has a wave-form denoted by xg and magnitude m; p(xg|m) is the probability of MS wave-form given magnitude m; and p(m) is the probability of having MS with a magnitude m. Based on the working assumption (see e.g., Jalayer et al. 2012) that different wave-forms are equally likely to occur, Eq. (4) can be approximated by the average of P(LSas|mi,xig) values over a set of Ngm ground motion records (which have not caused LS excursion for the intact structure). To estimate P(LSas|mi,xig) in Eq. (4), the aftershock risk assessment procedure derived by Jalayer et al. (2011a) and later refined in Ebrahimian et al. (2014b) is used herein. This procedure leads to the estimation of time-dependent limit state first-excursion probabilities given that the MS magnitude and its wave-form is known, taking into account both the time-varying aftershock occurrence rate and also the cumulative damage caused by the sequence of aftershocks. Let Nas be the maximum number of AS events that may take place in the time interval [0,t]. The limit state probability P(LSas|mi,xig) can be expanded over all possible AS events that can take place in the time-interval [0,t] based on Total probability Theorem (Benjamin and Cornell, 1970): ,, ,1| ,| , , | ,asasas i i gNas i as as ii g i gnP LS m xP LS m x n P n m x (5) where P(LSas|mi,xig,nas) is the probability of exceeding the limit state LS for the first time given that exactly nas AS events take place andmi and xig. The term P(nas|mi,xig) is the conditional probability that exactly nas AS events take place given MS magnitude and its wave-form. This probability is defined herein by a non-homogenous Poisson probability distribution with a time-decaying rate based on the Modified Omori (MO) aftershock occurrence model (Utsu 1961) (see also Jalayer et al. 2011a and Ebrahimian et al. 2014a). Let as be the number of AS events with magnitude Ml,as≤Mas≤Mu taken place in time interval [0,t] followed by a MS with magnitude mi. Assuming that the earthquake occurrence is described by a non-homogeneous Poisson process with parameter as, the probability P(nas|mi,xig) can be calculated as (note that we have dropped the dependence on MS wave-form): | !as asnasas iaseP n m n (6) Note that the parameter as is both time- and MS magnitude-dependent (see Ebrahimian et al. 2014a for more details). The parameters of the MO model can be obtained from a generic territorial model. The probability P(LSas|mi,xig,nas) in Eq. (5) can be calculated by taking into account the set of mutually exclusive and collectively exhaustive (MECE) events that LS first-excursion happens at one and just one of the nas AS events (Ebrahimian et al. 2014b): 1,1 1| , , 1asn ki gas i as k ik iP LS m x n (7) where k denotes the probability of LS first-excursion due to the occurrence of the kth AS event in the sequence given that the limit state 12th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12 Vancouver, Canada, July 12-15, 2015 4has not exceeded in the previous (k-1) events. The probability term k can be expanded with respect to the first-mode spectral acceleration Sa(T) as follows (Ebrahimian et al. 2014b): d ( )1 ddasas ask k as asas assasasa sasa (8) where saas is Sa(T) associated with the kth aftershock event; k=P(LSas|saas) is an event-dependent fragility for the kth aftershock event (defined as the probability of exceeding the limit state LS due to the kth event given that it has not been exceeded due to the previous k-1 events; and as is the mean (in unit of time t) rate of exceeding Sa(T) equal to saas (see Ebrahimian et al. 2014a for the description of aftershock hazard calculation). 2.2. Closed-form approximation of the time-dependent limit state probabilities In this section, the term P(LSas|mi,xig) in Eq. (5) is approximated with a closed-form analytic expression. A preliminary version of this closed-form expression was proposed by Jalayer et al. (2014a). Assume that the set of probability terms {k|k=1:nas} is identical and equal to the time- and event-invariant function . Thus, P(LSas|mi,xig,nas) in Eq. (7) is calculated as the sum of a geometric series: 11( | )1 1 1, ,asasn kia inskgs aP LS m x n (9) Substituting P(LSas|mi,xig,nas) from Eq. (9) and P(nas|mi) from Eq. (6) into Eq. (5), and performing some simple algebraic manipulations: 1( )( | ) 1 1 !,as asasasasnN n aigas isn asP LS m x en (10) The second term in Eq. (10) is a filtered Poisson probability mass function with a binary filter having probability . Thus, Eq. (10) can be simplified as follows: 11( |!, )11as asasasasnN asigasniasP LS m xee n (11) Finally, the closed-form approximation to the time-dependent limit state probability P(LSas|mi,xig) is derived as (considering that the sum of the Poisson probability mass function terms is asymptotically equal to unity for large Nas values): ( | ) 1 e, xpigas i asP LS m x (12) The closed-form expression in Eq. (12) provides a simple analytic equation for calculating the time-dependent LS first-excursion probability due a sequence of aftershocks triggered by a MS of magnitude mi and wave-form xig. An extremely convenient proposal for can be LSms,which is calculated by dividing the limit state exceedance rate LS due to the MS from Eq. (2) by the seismicity rate of MS events Ml≤Mms≤Mu, denoted by ms, taking care that LS and ms should have the same temporal unit. 3. NUMERICAL APPLICATION The methodology described in Section 2 is applied herein in order to perform risk assessment for (MS+AS sequence) for a typical RC building with infill panels located in L’Aquila, central Italy. 3.1. Case Study Definitions The case-study structure is a shear building model of a 3-story 2-bay RC two-dimensional frame having one-bay infill panel. The first-mode period of the building is equal to 0.27sec. A comprehensive representation of the model with the attributed nonlinear behavior of columns and infill panel can be found in Ebrahimian et al. (2014b). Since the damage is mainly accumulated in the first story, the displacement of the first story is taken to be the engineering demand parameter (EDP) in this study. The nonlinear dynamic analyses are performed using OpenSees (http://opensees.berkeley.edu, ver. 12th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12 Vancouver, Canada, July 12-15, 2015 52.4.4) by adopting properly calibrated hysteretic models. Only one discrete limit state of Near Collapse (NC) is considered in this study, which is conservatively set to 10% drop in ultimate strength of the columns based on the pushover analysis of the structure (see Ebrahimian et al. 2014b for more details). The latitudinal and longitudinal coordinates of the site where the reference structure is located are [42.3450N, 13.4009E]. This site is located within the seismic zone 923 based on the ZS9 Italian Seismogenetic Zonation (Gruppo di Lavoro, 2004), in which the key seismicity parameters are: ms=0.14 (/year), Ml=4.76, Mu=7.06. The aftershock zone considered herein is the one presented in Ebrahimian et al. (2014a), which is located within the prescribed latitude/longitude cells of a 0.1 0.1 grid of area [42.10-42.70N, 13.15-13.65E,]. The aftershock magnitudes are within the range of 4.0≤Mas≤7.06. The occurrence of aftershocks are modeled by the MO model proposed by Jalayer et al (2011a). Based on a Bayesian updating framework, they have updated the Italian generic parameters of MO model based on the sequence of aftershocks following the L’Aquila earthquake of 6th April 2009. The risk calculations are performed for a time-interval of 100 days after the MS. 3.2. Seismic Hazard Assessment The site-specific seismic hazard for the main events is extracted from the website of INGV (Istituto Nazionale di Geofisica e Vulcanologia, Progetto INGV-DPC-S1, http://esse1.mi.ingv.it/) based on the Italian seismic zonation. It provides the site-specific hazard associated with major seismic events within the Italian territory. The aftershock hazard, calculated as the probability of exceeding a given spectral acceleration value given that a MS event has taken place, is obtained by integrating an adopted ground-motion prediction equation over all possible aftershock magnitudes and distances within the desired aftershock zone. An updated version of the ground-motion prediction relation proposed by Sabetta and Pugliese (1996) (more suitable for small-magnitude events) whose parameters are updated using a Bayesian approach according to the L’Aquila aftershock sequence is used herein (Ebrahimian et al. 2014a). A uniform seismicity pattern over the aftershock zone is considered. 3.3. Calculation of Event-dependent Fragilities The methodology for calculating the event-dependent fragilities denoted as {πk|k=1:nas} is comprehensively discussed in the previous works by the authors (Ebrahimian et al. 2014b). The structural performance variable denoted as YLS is calculated as following: 1max1k kk rLS kLS rD DY C D (13) where Dmax(k) is the maximum EDP due to the kth event; Dr(k-1) is the associated residual demand corresponding to the sequence of (k-1) events; CLS is the limit state capacity of the (intact) structure. This history-dependent parameter, which is the ratio of maximum demand increment to a reduced limit state capacity, reveals appropriate correlation with Sa(T) of the kth event within the sequence. For performing the Cloud Analysis (see Elefante et al. 2010 and Jalayer et al. 2014b), a set of 50 European (especially Italian) ground motion records, including mainshock ground-motions as well as aftershocks, are selected from the NGA-West2 database (Ancheta et al. 2014). The set of records covers a wide range of magnitudes from 4.5 up to 7.5, and closest distance to ruptured area up to 80 km. It is to note that accurate record selection can significantly affect the results (see e.g. Goda 2014); however, this issue is not a primary focus in this study. The set of event-dependent fragility curves {πk|k=1:13} are shown in Figure 1 together with the cloud regression associated with π0, which is the fragility of the intact structure. For simplicity of calculations when k>13, the event-dependent fragility πk=π13 is considered for calculation of k in Eq. (8). 12th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12 Vancouver, Canada, July 12-15, 2015 6 Figure 1: (a) Event-dependent fragilities πk, k=1:13, (b) Cloud regression leading to π0 3.4. Time-dependent Risk Calculation The time-dependent limit state probability due to (MS+AS sequence) =P(LS|I1) is calculated according to Eq. (1). The limit state exceedance probability in time [0,t] due to MS denoted by ms is calculated according to Eq. (3). The first-excursion limit state probability due to the AS sequence, as in Eq. (1), is calculated from Eq. (4) by integrating the post-mainshock limit state probability given MS magnitude and wave-form P(LSas|m,xg) over all possible MS magnitudes m and wave-forms xg. As mentioned before, this expression is approximated as the average of P(LSas|mi,xig) values for a suite of strong motion records which have not led to limit state excursion of the intact structure. P(LSas|mi,xig) is calculated according to Eq. (5). With reference to Eq. (5), P(nas|mi,xig) is estimated by the Poisson probability distribution from Eq. (6), and the limit state probability P(LSas|mi,xig,nas) (due to aftershocks) given the MS magnitude mi, wave-form xig and the number of aftershock events nas can be calculated from the recursive formulation in Eq. (7) as a function of {k|k=1:nas}. Accordingly, k is the probability of LS first-excursion due to kth aftershock event given that the limit state has not exceeded in the previous (k-1) events, and is estimated from Eq. (8). The probabilities ms (red dotted line), and the “best-estimate” for P (calculated through the procedure described in Section 2 and plotted as black solid lines) are all plotted in Figure 2. It can be observed that the consideration of the cumulative damage caused by the aftershocks for the studied case (building severely damaged by the MS), will significantly increase the estimated time-variant risk . Now, it is interesting to see how the derived analytical closed-form for P(LSas|mi,xig) in Eq. (12) estimates the effect of the cumulative damage caused by the aftershocks. The resulting “approximate” P is plotted in blue dashed line in Figure 2. It can be observed that, although obtained with significantly reduced analytic effort (almost the same analytical effort as the one needed for obtaining ms), the “approximate” time-dependent limit state probability P roughly manages to capture the effect of the cumulative damage caused by the AS sequence. It should be mentioned that if the MS damaged fragility is used in Eq. (8) (instead of ), the “approximate” P becomes almost identical to the “best-estimate” (it is not shown in Figure 2). The probability of exceeding the limit state in a reference time interval [t, t+t] (say e.g. t=1 day) can also be calculated as the difference of the exceedance probabilities in intervals [0,t+t] and [0,t] (see Jalayer et al. 2011a, b). The resulting daily limit state probabilities as a function of time t for t=1 day are plotted in Figure 3. Assuming equal number of events in 1 year between two Poisson 0 0.5 1 1.5 2 2.500.10.20.30.40.50.60.70.80.91Sa(T=0.27 sec) [g]ProbabilityEvent-dependent Fragilities fragility 0fragility 1fragility 2fragility 3fragility 4fragility 5fragility 6fragility 7fragility 8fragility 9fragility 10fragility 11fragility 12fragility 1310-3 10-2 10-1 100 10110-210-1100101YLS (0)Sa(T=0.27 sec) [g] cloud datacloud regression(a) (b) 12th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12 Vancouver, Canada, July 12-15, 2015 7processes with daily and annual unit times, we have calculated the equivalent admissible mean daily rate of exceedance equal to 5.6 ൈ 10-6 corresponding to a mean annual admissible rate of exceedance equal to 2ൈ10-3 (475 year return period). For the first 100 days elapsed after the MS, the daily limit state probability is decisively above the admissible threshold (in green), while it is below the green line considering only MS (red dotted line). A very close agreement between the “best-estimate” and the “approximate” results can be observed, as well. Figure 2: Comparison between the estimated risk by considering only MS and the MS+AS sequence Figure 3: Time-dependent daily risk 4. CONCLUSIONS The findings in this work can be summarized as follows: (1) A procedure for calculating the limit state exceedance probability due to (MS+AS sequence) is proposed. (2) This methodology explicitly takes into account the cumulative damage caused by the AS sequence for MDOF structures with degrading behavior. (3) A closed-form analytical expression is derived for calculating the post-mainshock limit state probability. The simplest application of this analytic formulation based on the fragility of intact structure leads to a much better agreement with the computationally extensive “best-estimate” results than expected. (4) For the particular building studied herein, considering the triggered AS sequence significantly affects the limit state probability in the immediate aftermath of the main event. The obtained results seem to be consistent with those obtained in the previous work by the authors (Ebrahimian et al. 2014b). (5) It is expected that the results are sensitive to the amount of damage sustained by the building due to the MS and the characteristics of the AS sequence. The proposed procedure is based on the strong working assumption that various MS events are equally likely to take place. It should also be reminded that the Cloud Analysis, employed herein for the calculation of event-dependent fragilities, is quite sensitive to the selection of records. 5. AKNOWLEDGEMENT This work was supported in part by project STRIT (Strumenti e Tecnologie per la gestione del Rischio delle Infrastrutture di Trasporto). This support is gratefully acknowledged. 6. REFERENCES Ancheta, T.D., et al. (2014). “NGA-West2 Database.” Earthquake Spectra, 30(3), 989-1005. 0 10 20 30 40 50 60 70 80 90 10010-610-510-410-310-210-1100t [day]ProbabilityProbability of exceeding the Limit State NC in time t PmsP, best-estimateP, approximate0 10 20 30 40 50 60 70 80 90 10010-710-610-510-410-310-210-1100t [day]ProbabilityDaily Probability of exceeding the Limit State NC in time t admissible levelconsidering only MSbest-estimateapproximate12th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12 Vancouver, Canada, July 12-15, 2015 8Benjamin, J.R., and Cornell, C.A. (1970). “Probability, statistics and decision for civil engineers.” McGraw Hill, 684pp. Ebrahimian, H., Jalayer, F., Asprone, D., Lombardi, A.M., Marzocchi, W., Prota, A., Manfredi, G. (2014a). “Adaptive Daily Forecasting of Seismic Aftershock Hazard.” Bulletin of Seismological Society of America, 104(1), 145-161. Ebrahimian, H., Jalayer, F., Asprone, D., Lombardi, A.M., Marzocchi, W., Prota, A., Manfredi, G. (2014b). “A performance-based framework for adaptive seismic aftershock risk assessment.” Earthquake Engineering and Structural Dynamics, 43(14), 2179-2197. Elefante, L., Jalayer, F., Iervolino, I., Manfredi, G. (2010) “Disaggregation-based response weighting scheme for seismic risk assessment of structures.” Soil Dynamics and Earthquake Engineering, 30(12), 1513-1527. Gruppo di Lavoro. (2004). Redazione della mappa di pericolosità sismica prevista dall’Ordinanza PCM 3274 del 20 marzo 2003. Rapporto Conclusivo per il Dipartimento della Protezione Civile, INGV, Milano-Roma, April 2004; 65 pp. + 5 appendixes (in Italian). Goda, K., and Taylor, C.A. (2012). “Effects of aftershocks on peak ductility demand due to strong ground motion records from shallow crustal earthquakes.” Earthquake Engineering and Structural Dynamics, 41(15), 2311-2330. Goda, K., and Salami, R. (2014). “Inelastic seismic demand estimation of wood-frame houses subjected to mainshock-aftershock sequences.” Bulletin of Earthquake Engineering, 12(2), 855-874. Goda, K. (2014). “Record selection for aftershock incremental dynamic analysis.” Earthquake Engineering and Structural Dynamics, DOI:10.1002/eqe.2513. Jalayer, F., Asprone, D., Prota, A., Manfredi, G.A. (2011a). “A decision support system for post-earthquake reliability assessment of structures subjected to aftershocks: an application to L’Aquila earthquake, 2009.” Bulletin of Earthquake Engineering, 9(4), 997-1014. Jalayer, F., Asprone, D., Prota, A., Manfredi, G. (2011b). “Multi-hazard upgrade decision making for critical infrastructure based on life-cycle cost criteria.” Earthquake Engineering and Structural Dynamics, 40(10), 1163-1179. Jalayer, F., Beck, J., Zareian, F. (2012). “Analyzing the sufficiency of alternative scalar and vector intensity measures of ground shaking based on information theory.” Journal of Engineering Mechanics, 138(3), 307-316. Jalayer, F., Ebrahimian, H., Asprone, D., Prota, A., Manfredi, G. (2014a). A closed-from approximation to adaptive structural fragility to aftershocks. Proceedings of the 10th National Conference in Earthquake Engineering, Earthquake Engineering Research Institute, Anchorage, AK. Jalayer, F., De Risi, R., Manfredi, G. (2014b). “Bayesian Cloud Analysis: efficient structural fragility assessment using linear regression.” Bulletin of Earthquake Engineering, DOI 10.1007/s10518-014-9692-z. Luco, N., Bazzurro, P., Cornell, C.A. (2004). Dynamic versus static computation of the residual capacity of mainshock-damaged building to withstand an aftershock. In: Proceedings 13th World Conference on Earthquake Engineering, Vancouver, B.C., Canada, August 1-6, Paper No. 2405. Raghunandan, M., Liel, A.B., and Luco, N. (2014) “Aftershock collapse vulnerability assessment of reinforced concrete frame structures.” Earthquake Engineering and Structural Dynamics, DOI:10.1002/eqe.2478. Sabetta, F., Pugliese, A. (1996). “Estimation of response spectra and simulation of nonstationary earthquake ground motions.” Bulletin of Seismological Society of America, 86(2), 337-352. Utsu, T. (1961). “A statistical study of the occurrence of aftershocks.” Geophysical Magazine, 30, 521-605. Yeo, G.L., and Cornell, C.A. (2009). “Equivalent constant rates for post-quake seismic decision making.” Structural Safety, 31(5), 443-447.
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- International Conference on Applications of Statistics and Probability in Civil Engineering (ICASP) (12th : 2015) /
- Towards quantifying the effect of aftershocks in seismic...
Open Collections
International Conference on Applications of Statistics and Probability in Civil Engineering (ICASP) (12th : 2015)
Towards quantifying the effect of aftershocks in seismic risk assessment Jalayer, Fatemeh; Ebrahimian, Hossein; Manfredi, Gaetano Jul 31, 2015
pdf
Page Metadata
Item Metadata
Title | Towards quantifying the effect of aftershocks in seismic risk assessment |
Creator |
Jalayer, Fatemeh Ebrahimian, Hossein Manfredi, Gaetano |
Contributor | International Conference on Applications of Statistics and Probability (12th : 2015 : Vancouver, B.C.) |
Date Issued | 2015-07 |
Description | The state of the art in seismic design and assessment of structures implicitly relies on structural redundancy in order to deal with the effect of the triggered sequence of aftershock events on a building. Calculating the time-dependent limit state exceedance probability for a structure considering both the main event and the triggered sequence of aftershocks is complicated both by the time-dependent rate of aftershock occurrence and also the cumulative damage caused by the sequence of events. Taking advantage of a methodology developed in previous works by the authors for post-mainshock risk assessment, the limit state probability due to a sequence of a mainshock and the triggered aftershocks is calculated herein. Moreover, a closed- and analytic-form approximation to the post-mainshock limit state probability given the mainshock magnitude is derived. This closed- and analytic formulation facilitates the estimation of the limit state probability by employing the standard tools in risk assessment such as the fragility curve for the intact structure. Applying the proposed methodology for a reinforced concrete (RC) concrete moment-resisting frame with infills subjected to the main event and the triggered sequence emphasizes the importance of taking into account the cumulative damage caused by the triggered sequence for the case-study model. Moreover, it is demonstrated that the proposed analytical approximation based on the fragility curve for the intact structure leads to a surprisingly close agreement with the best-estimate results obtained by considering the time- and event-dependent degradation in the structure. |
Genre |
Conference Paper |
Type |
Text |
Language | eng |
Notes | This collection contains the proceedings of ICASP12, the 12th International Conference on Applications of Statistics and Probability in Civil Engineering held in Vancouver, Canada on July 12-15, 2015. Abstracts were peer-reviewed and authors of accepted abstracts were invited to submit full papers. Also full papers were peer reviewed. The editor for this collection is Professor Terje Haukaas, Department of Civil Engineering, UBC Vancouver. |
Date Available | 2015-05-21 |
Provider | Vancouver : University of British Columbia Library |
Rights | Attribution-NonCommercial-NoDerivs 2.5 Canada |
DOI | 10.14288/1.0076186 |
URI | http://hdl.handle.net/2429/53285 |
Affiliation |
Non UBC |
Citation | Haukaas, T. (Ed.) (2015). Proceedings of the 12th International Conference on Applications of Statistics and Probability in Civil Engineering (ICASP12), Vancouver, Canada, July 12-15. |
Peer Review Status | Unreviewed |
Scholarly Level | Faculty |
Rights URI | http://creativecommons.org/licenses/by-nc-nd/2.5/ca/ |
Aggregated Source Repository | DSpace |
Download
- Media
- 53032-Paper_307_Jalayer.pdf [ 247.49kB ]
- Metadata
- JSON: 53032-1.0076186.json
- JSON-LD: 53032-1.0076186-ld.json
- RDF/XML (Pretty): 53032-1.0076186-rdf.xml
- RDF/JSON: 53032-1.0076186-rdf.json
- Turtle: 53032-1.0076186-turtle.txt
- N-Triples: 53032-1.0076186-rdf-ntriples.txt
- Original Record: 53032-1.0076186-source.json
- Full Text
- 53032-1.0076186-fulltext.txt
- Citation
- 53032-1.0076186.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.53032.1-0076186/manifest