skip to main content

Attention:

The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 10:00 PM ET on Friday, December 8 until 2:00 AM ET on Saturday, December 9 due to maintenance. We apologize for the inconvenience.


Title: Surrogate modeling of structural seismic response using probabilistic learning on manifolds
Abstract

Nonlinear response history analysis (NLRHA) is generally considered to be a reliable and robust method to assess the seismic performance of buildings under strong ground motions. While NLRHA is fairly straightforward to evaluate individual structures for a select set of ground motions at a specific building site, it becomes less practical for performing large numbers of analyses to evaluate either (1) multiple models of alternative design realizations with a site‐specific set of ground motions, or (2) individual archetype building models at multiple sites with multiple sets of ground motions. In this regard, surrogate models offer an alternative to running repeated NLRHAs for variable design realizations or ground motions. In this paper, a recently developed surrogate modeling technique, called probabilistic learning on manifolds (PLoM), is presented to estimate structural seismic response. Essentially, the PLoM method provides an efficient stochastic model to develop mappings between random variables, which can then be used to efficiently estimate the structural responses for systems with variations in design/modeling parameters or ground motion characteristics. The PLoM algorithm is introduced and then used in two case studies of 12‐story buildings for estimating probability distributions of structural responses. The first example focuses on the mapping between variable design parameters of a multidegree‐of‐freedom analysis model and its peak story drift and acceleration responses. The second example applies the PLoM technique to estimate structural responses for variations in site‐specific ground motion characteristics. In both examples, training data sets are generated for orthogonal input parameter grids, and test data sets are developed for input parameters with prescribed statistical distributions. Validation studies are performed to examine the accuracy and efficiency of the PLoM models. Overall, both examples show good agreement between the PLoM model estimates and verification data sets. Moreover, in contrast to other common surrogate modeling techniques, the PLoM model is able to preserve correlation structure between peak responses. Parametric studies are conducted to understand the influence of different PLoM tuning parameters on its prediction accuracy.

 
more » « less
Award ID(s):
2131111
NSF-PAR ID:
10418905
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Earthquake Engineering & Structural Dynamics
Volume:
52
Issue:
8
ISSN:
0098-8847
Page Range / eLocation ID:
p. 2407-2428
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. : In order to evaluate urban earthquake resilience, reliable structural modeling is needed. However, detailed modeling of a large number of structures and carrying out time history analyses for sets of ground motions are not practical at an urban scale. Reduced-order surrogate models can expedite numerical simulations while maintaining necessary engineering accuracy. Neural networks have been shown to be a powerful tool for developing surrogate models, which often outperform classical surrogate models in terms of scalability of complex models. Training a reliable deep learning model, however, requires an immense amount of data that contain a rich input-output relationship, which typically cannot be satisfied in practical applications. In this paper, we propose model-informed symbolic neural networks (MiSNN) that can discover the underlying closed-form formulations (differential equations) for a reduced-order surrogate model. The MiSNN will be trained on datasets obtained from dynamic analyses of detailed reinforced concrete special moment frames designed for San Francisco, California, subject to a series of selected ground motions. Training the MiSNN is equivalent to finding the solution to a sparse optimization problem, which is solved by the Adam optimizer. The earthquake ground acceleration and story displacement, velocity, and acceleration time histories will be used to train 1) an integrated SNN, which takes displacement and velocity states and outputs the absolute acceleration response of the structure; and 2) a distributed SNN, which distills the underlying equation of motion for each story. The results show that the MiSNN can reduce computational cost while maintaining high prediction accuracy of building responses. 
    more » « less
  3. null (Ed.)
    Nonstructural components within mission-critical facilities such as hospitals and telecommunication facilities are vital to a community's resilience when subjected to a seismic event. Building contents like medical and computer equipment are critical for the response and recovery process following an earthquake. A solution to protecting these systems from seismic hazards is base isolation. Base isolation systems are designed to decouple an entire building structure from destructive ground motions. For other buildings not fitted with base isolation, a practical and economical solution to protect vital building contents from earthquake-induced floor motion is to isolate individual equipment using, for example, rolling-type isolation systems (RISs). RISs are a relatively new innovation for protecting equipment. These systems function as a pendulum-like mechanism to convert horizontal motion into vertical motion. An accompanying change in potential energy creates a restoring force related to the slope of the rolling surface. This study seeks to evaluate the seismic hazard mitigation performance of RISs, as well as propose and test a novel double RIS. A physics-based mathematical model was developed for a single RIS via Lagrange's equation adhering to the kinetic constraint of rolling without slipping. The mathematical model for the single RIS was used to predict the response and characteristics of these systems. A physical model was fabricated with additive manufacturing and tested against multiple earthquakes on a shake table. The system featured a single-degree-of-freedom (SDOF) structure to represent a piece of equipment. The results showed that the RIS effectively reduced accelerations felt by the SDOF compared to a fixed-base SDOF system. The single RIS experienced the most substantial accelerations from the Mendocino record, which contains low-frequency content in the range of the RIS's natural period (1-2 seconds). Earthquakes with these long-period components have the potential to cause impacts within the isolation bearing that would degrade its performance. To accommodate large displacements, a double RIS is proposed. The double RIS has twice the displacement capacity of a single RIS without increasing the size of the bearing components. The mathematical model for the single RIS was extended to the double RIS following a similar procedure. Two approaches were used to evaluate the double RIS's performance: stochastic and deterministic. The stochastic response of the double RIS under stationary white noise excitation was evaluated for relevant system parameters, namely mass ratio and tuning frequency. Both broadband and filtered (Kanai-Tajimi) white noise excitation were considered. The response variances of the double RIS were normalized by a baseline single RIS for a comparative study, from which design parameter maps were drawn. A deterministic analysis was conducted to further evaluate the double RIS in the case of nonstationary excitation. The telecommunication equipment qualification waveform, VERTEQ-II, was used for these numerical simulations. Peak transient responses were compared to the single RIS responses, and optimal design regions were determined. General design guidelines based on the stochastic and deterministic analyses are given. The results aim to provide a framework usable in the preliminary design stage of a double RIS to mitigate seismic responses. 
    more » « less
  4. 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
  5. Summary

    Robust estimation of collapse risk should consider the uncertainty in modeling of structures as well as variability in earthquake ground motions. In this paper, we illustrate incorporation of the uncertainty in structural model parameters in nonlinear dynamic analyses to probabilistically assess story drifts and collapse risk of buildings. Monte Carlo simulations with Latin hypercube sampling are performed on ductile and non‐ductile reinforced concrete building archetypes to quantify the influence of modeling uncertainties and how it is affected by the ductility and collapse modes of the structures. Inclusion of modeling uncertainty is shown to increase the mean annual frequency of collapse by approximately 1.8 times, as compared with analyses neglecting modeling uncertainty, for a high‐seismic site. Modeling uncertainty has a smaller effect on drift demands at levels usually considered in building codes; for the same buildings, modeling uncertainty increases the mean annual frequency of exceeding story drift ratios of 0.03 by 1.2 times. A novel method is introduced to relate drift demands at maximum considered earthquake intensities to collapse safety through a joint distribution of deformation demand and capacity. This framework enables linking seismic performance goals specified in building codes to drift limits and other acceptance criteria. The distributions of drift demand at maximum considered earthquake and capacity of selected archetype structures enable comparisons with the proposed seismic criteria for the next edition (2016) of ASCE 7. Subject to the scope of our study, the proposed drift limits are found to be unconservative, relative to the target collapse safety in ASCE 7. Copyright © 2016 John Wiley & Sons, Ltd.

     
    more » « less