skip to main content


Title: An earthquake-source-based metric for seismic fragility analysis
The seismic fragility of a system is the probability that the system enters a damage state under seismic ground motions with specified characteristics. Plots of the seismic fragilities with respect to scalar ground motion intensity measures are called fragility curves. Recent studies show that fragility curves may not be satisfactory measures for structural seismic performance, since scalar intensity measures cannot comprehensively characterize site seismicity. The limitations of traditional seismic intensity measures, e.g., peak ground acceleration or pseudo-spectral acceleration, are shown and discussed in detail. A bivariate vector with coordinates moment magnitude m and source-to-site distance r is proposed as an alternative seismic intensity measure. Implicitly, fragility surfaces in the (m, r)-space could be used as graphical representations of seismic fragility. Unlike fragility curves, which are functions of scalar intensity measures, fragility surfaces are characterized by two earthquake-hazard parameters, (m, r). The calculation of fragility surfaces may be computationally expensive for complex systems. Thus, as solutions to this issue, a bi-variate log-normal parametric model and an efficient calculation method, based on stochastic-reduced-order models, for fragility surfaces are proposed.  more » « less
Award ID(s):
1639669
NSF-PAR ID:
10060257
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Bulletin of Earthquake Engineering
ISSN:
1570-761X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Loss of operation or devastating damage to buildings and industrial structures, as well as equipment housed in them, has been observed due to earthquake-induced vibrations. A common source of operational downtime is due to the performance reduction of vital equipment, which are sensitive to the total transmitted acceleration. A well-known method of protecting such equipment is seismic isolation of the equipment itself (or a group of equipment), as opposed to the entire structure due to the lower cost of implementation. The first objective of this dissertation is assessing a rolling isolation system (RIS) based on existing design guidelines for telecommunications equipment. A discrepancy is observed between the required response spectrum (RRS) and the one and only accelerogram recommended in the guideline. Several filters are developed to generate synthetic accelerograms that are compatible with the RRS. The generated accelerograms are used for probabilistic assessment of a RIS that is acceptable per the guideline. This assessment reveals large failure probability due to displacement demands in excess of the displacement capacity of the RIS. When the displacement demands on an isolation system are in excess of its capacity, impacts result in spikes in transmitted acceleration. Therefore, the second objective of this dissertation is to design impact prevention/mitigation mechanisms. A dual-mode system is proposed where the behavior changes when the displacement exceeds a predefined threshold. A new piecewise optimal control approach is developed and applied to find the best possible mechanism for the region beyond the threshold. By utilizing the designed curves obtained from the proposed optimal control procedure, a Kelvin-Voigt device is tuned for illustrative purposes. On the other hand, the preference for protecting equipment decreases as the earthquake intensity increases. In extreme seismic loading, the response mitigation of the primary structure (i.e., life safety and collapse prevention) is of greater concern than protecting isolated equipment. Therefore, the third objective of this dissertation is to develop an innovative dual-mode system that can behave as equipment isolation under low to moderate seismic loading and passively transition to behave as a vibration absorber for the primary structure under extreme seismic loading. To reduce the computational cost of simulating a large linear elastic structure with nonlinear attachments (i.e., equipment isolation with cubic hardening nonlinearity), a reduced order modeling method is introduced that can capture the behavior of such nonlinear coupled systems. The method is applied to study the feasibility of dual-mode vibration isolation/absorber. To this end, nonlinear transmissibility curves for the roof displacement and isolated mass total acceleration are developed from the steady-state responses of dual-mode systems using the harmonic balanced method. The final objective of this dissertation is to extend the reduced order modeling method developed for linear elastic structure with nonlinear attachment to inelastic structures (without attachments). The new inelastic model condensation (IMC) method uses the modal properties of the full structural model (in the elastic range) to construct a linear reduced order model in conjunction with a hysteresis model to capture the hysteretic inter-story restoring forces. The parameters of these hysteretic forces are easily tuned, in order to fit the inelastic behavior of the condensed structure to that of the full model under a variety of simple loading scenarios. The fidelity of structural models condensed in this way is demonstrated via simulation for different ground motion intensities on three different building structures with various heights. The simplicity, accuracy, and efficiency of this approach could significantly alleviate the computational burden of performance-based earthquake engineering. 
    more » « less
  2. null (Ed.)
    We present a database and analyze ground motions recorded during three events that occurred as part of the July 2019 Ridgecrest earthquake sequence: a moment magnitude (M) 6.5 foreshock on a left‐lateral cross fault in the Salt Wells Valley fault zone, an M 5.5 foreshock in the Paxton Ranch fault zone, and the M 7.1 mainshock, also occurring in the Paxton Ranch fault zone. We collected and uniformly processed 1483 three‐component recordings from an array of 824 sensors spanning 10 seismographic networks. We developed site metadata using available data and multiple models for the time‐averaged shear‐wave velocity in the upper 30 m (⁠VS30⁠) and for basin depth terms. We processed ground motions using Next Generation Attenuation (NGA) procedures and computed intensity measures including spectral acceleration at a number of oscillator periods and inelastic response spectra. We compared elastic and inelastic response spectra to seismic design spectra in building codes to evaluate the damage potential of the ground motions at spatially distributed sites. Residuals of the observed spectral accelerations relative to the NGA‐West2 ground‐motion models (GMMs) show good average agreement between observations and model predictions (event terms between about −0.3 and 0.5 for peak ground acceleration to 5 s). The average attenuation with distance is also well captured by the empirical NGA‐West2 GMMs, although azimuthal variations in attenuation were observed that are not captured by the GMMs. An analysis considering directivity and fault‐slip heterogeneity for the M 7.1 event demonstrates that the dispersion in the near‐source ground‐motion residuals can be reduced. 
    more » « less
  3. Summary

    Ground motion selection is a crucial step in probabilistic seismic performance assessment of structural systems. Particularly, identifying ground motion records compatible with a specific hazard level has been the major focus of past studies. We propose a multivariate return period (MRP)‐based record selection methodology to further improve the multivariate hazard consistency over a vector of intensity measures (IMs). Unlike the traditional univariate return period anchored on a scalar IM, MRP generalizes the return period concept by accommodating the joint rate of exceedance of a vector of IMs, thereby providing more holistic seismic hazard characterization. By leveraging MRP in linking the seismic hazard to a vector IMs, the proposed MRP‐based ground motion selection methodology for the first time offers a mathematically rigorous yet practical solution for multivariate hazard consistency in ground motion selection. The merit of the MRP‐based ground motion selection is demonstrated by comparing the resulting target spectra and seismic demand estimates for several case‐study structures with other state‐of‐the‐art ground motion selection alternatives. From the results, the MRP‐based ground motion selection employing Kendall's distribution function turns out to be a promising alternative, offering favorable new features including (a) moderate target spectra intensity and moderately low target spectra standard deviation; (b) superior convergence with the increase of the conditioning IM dimension and ability to approximate higher‐dimensional multivariate hazard consistency with lower dimensional conditioning IMs; and (c) capability to realistically capture the multimodal spectral shape owing to the incorporation of multivariate Gaussian mixture distribution in generating target spectra.

     
    more » « less
  4. Seismic isolation systems for buildings are generally selected to achieve higher seismic performance objectives, such as continued operation or immediate occupancy following a design earthquake event. However, recent large scale tests have suggested that these objectives may be compromised if the shaking includes large vertical acceleration components that are damaging to the nonstructural components and contents. Some research has been conducted to develop three dimensional isolation systems that can isolate the structure from both the horizontal and vertical components of ground motion. In several cases, systems have been proposed without much justification of the target design parameters. Rocking has been noted as a potential concern for structures with 3D isolation systems, and complex systems have been proposed to control the rocking. In this study, the fundamental dynamic response of structures with 3D isolation systems is explored. Target horizontal and vertical spectra for a representative strong motion site were developed based on NEHRP recommendations, and horizontal and vertical ground motions were selected that best fit the target spectra when the same amplitude scale factor was applied to all three motion components. Using a simple model of a rigid block resting on linear isolation bearings, the following aspects are evaluated for a wide range of horizontal and vertical isolation periods: response modes and severity of rocking, horizontal and vertical displacement demands in the isolation bearings, and attenuation of both horizontal and vertical accelerations in the structure relative to the ground acceleration. Preliminary results point to a number of useful observations. For example, rocking appears to be an issue only if the horizontal and vertical isolation periods are closely spaced. Helical spring isolation systems that have been applied to a few structures have this characteristic. However, if the horizontal isolation period is large relative to the vertical isolation period, troublesome rocking can be avoided. In addition, other researchers have proposed systems with vertical isolation periods on the order of 2 seconds, which require large displacement and damping capacity. However, preliminary results suggest that vertical isolation periods as low as 0.5 seconds will be effective in attenuating the vertical acceleration. Limiting the vertical isolation period will make design of a 3D isolation system more feasible with respect to vertical displacement capacity and avoiding rocking. 
    more » « less
  5. null (Ed.)
    Rapidly growing societal needs in urban areas are increasing the demand for tall buildings with complex structural systems. Many of these buildings are located in areas characterized by high seismicity. Quantifying the seismic resilience of these buildings requires comprehensive fragility assessment that integrates iterative nonlinear dynamic analysis (NDA). Under these circumstances, traditional finite element (FE) analysis may become impractical due to its high computational cost. Soft-computing methods can be applied in the domain of NDA to reduce the computational cost of seismic fragility analysis. This study presents a framework that employs nonlinear autoregressive neural networks with exogenous input (NARX) in fragility analysis of multi-story buildings. The framework uses structural health monitoring data to calibrate a nonlinear FE model. The model is employed to generate the training dataset for NARX neural networks with ground acceleration and displacement time histories as the input and output of the network, respectively. The trained NARX networks are then used to perform incremental dynamic analysis (IDA) for a suite of ground motions. Fragility analysis is next conducted based on the results of the IDA obtained from the trained NARX network. The framework is illustrated on a twelve-story reinforced concrete building located at Oklahoma State University, Stillwater campus. 
    more » « less