Extended MSRS rule for seismic analysis of bridges subjected to differential support motions

This paper introduces a generalized formulation of the multiple support response spectrum (MSRS) method (Earthquake Engng Struct. Dyn. 1992; 21 :713–740) and extends it by accounting for the quasi‐static contributions of truncated modes. The generalized formulation allows consideration of response quantities that involve support degrees of freedom (DOF). This situation arises for many response quantities of interest when rotational DOF are condensed out. New cross‐correlation coefficients are introduced in the extended rule and a parametric study is performed to gain insight and identify cases of ground motion spatial variability in which these terms are significant. An efficient computer implementation of the extended MSRS method is described and used for comprehensive analysis of two real bridge models with vastly different structural characteristics. The specified input is in accordance with standards used in engineering practice. The effects of differential support motions, including the influence of spatially varying soil conditions, on the pseudo‐static and dynamic components and the total response are examined and the improvement achieved with the extended MSRS method is assessed. Copyright © 2011 John Wiley & Sons, Ltd.     

Least‐squares reverse‐time migration in a matrix‐based formulation

This paper describes least‐squares reverse‐time migration. The method provides the exact adjoint operator pair for solving the linear inverse problem, thereby enhancing the convergence of gradient‐based iterative linear inversion methods. In this formulation, modified source wavelets are used to correct the source signature imprint in the predicted data. Moreover, a roughness constraint is applied to stabilise the inversion and reduce high‐wavenumber artefacts. It is also shown that least‐squares migration implicitly applies a deconvolution imaging condition. Three numerical experiments illustrate that this method is able to produce seismic reflectivity images with higher resolution, more accurate amplitudes, and fewer artefacts than conventional reverse‐time migration. The methodology is currently feasible in 2‐D and can naturally be extended to 3‐D when computational resources become more powerful.     

Approximate analysis methods for asymmetric plan base‐isolated buildings

An approximate method for linear analysis of asymmetric‐plan, multistorey buildings is specialized for a single‐storey, base‐isolated structure. To find the mode shapes of the torsionally coupled system, the Rayleigh–Ritz procedure is applied using the torsionally uncoupled modes as Ritz vectors. This approach reduces to analysis of two single‐storey systems, each with vibration properties and eccentricities (labelled ‘effective eccentricities’) similar to corresponding properties of the isolation system or the fixed‐base structure. With certain assumptions, the vibration properties of the coupled system can be expressed explicitly in terms of these single‐storey system properties. Three different methods are developed: the first is a direct application of the Rayleigh–Ritz procedure; the second and third use simplifications for the effective eccentricities, assuming a relatively stiff superstructure. The accuracy of these proposed methods and the rigid structure method in determining responses are assessed for a range of system parameters including eccentricity and structure flexibility. For a subset of systems with equal isolation and structural eccentricities, two of the methods are exact and the third is sufficiently accurate; all three are preferred to the rigid structure method. For systems with zero isolation eccentricity, however, all approximate methods considered are inconsistent and should be applied with caution, only to systems with small structural eccentricities or stiff structures. Copyright © 2001 John Wiley & Sons, Ltd.     

Seismic response of multi‐supported structures by proper orthogonal decomposition

The seismic analysis of structures is usually carried out considering the ground motion as fully‐correlated in space and determining the structural response by pseudo‐deterministic methods such as the response spectrum technique. Actually, the partial correlation of the seismic acceleration may influence heavily the behaviour of spatially extended structures, such as bridges, viaducts or pipelines. In order to take its partial correlation into account, the seismic ground motion is schematized as a stochastic process dependent on time and on space; the hypotheses of stationarity and homogeneity are used to obtain simple and general results. The influence of the partial correlation of the seismic ground motion on the structural response is investigated by introducing suitable Equivalent Spectra. The acceleration of the support‐points of the structure is represented by the Proper Orthogonal Decomposition (POD), defining the modes of the earthquake. The method is formulated for any kind of multi‐degree‐of‐freedom system and is applied, as a case study, to an ideal single‐storey multi‐supported frame with an axially rigid beam. In the case of two supports, the POD decouples the pseudo‐static and the dynamic contributions to the structural response. This property is preserved for structural systems with many supports, where only the lower modes of the earthquake, usually the first two POD modes, are responsible for the structural response. Copyright © 2003 John Wiley & Sons, Ltd.     

LMS‐based structural health monitoring of a non‐linear rocking structure

A structure's health or level of damage can be monitored by identifying changes in structural or modal parameters. This research directly identifies changes in structural stiffness due to modelling error or damage for a post‐tensioned pre‐cast reinforced concrete frame building with rocking beam column connections and added damping and stiffness (ADAS) elements. A structural health monitoring (SHM) method based on adaptive least mean squares (LMS) filtering theory is presented that identifies changes from a simple baseline model of the structure. This method is able to track changes in the stiffness matrix, identifying when the building is (1) rocking, (2) moving in a hybrid rocking–elastic regime, or (3) responding linearly. Results are compared for two different LMS‐based SHM methods using an L ₂ error norm metric. In addition, two baseline models of the structure, one using tangential stiffness and the second a more accurate bi‐linear stiffness model, are employed. The impact of baseline model complexity is then delineated. The LMS‐based methods are able to track the non‐linearity of the system to within 15% using this metric, with the error due primarily to filter convergence rates as the structural response changes regimes while undergoing the El Centro ground motion record. The use of a bi‐linear baseline model for the SHM problem is shown to result in error metrics that are at least 50% lower than those for the tangential baseline model. Errors of 5–15% with this L ₂ error norm are fairly stringent compared to the greater than 2 × changes in stiffness undergone by the structure, however, in practice the usefulness of the results is dependent on the resolution required by the user. The impact of sampling rate is shown to be negligible over the range of 200–1000Hz, along with the choice of LMS‐based SHM method. The choice of baseline model and its level of knowledge about the actual structure is seen to be the dominant factor in achieving good results. The methods presented require 2.8–14.0 Mcycles of computation and therefore could easily be implemented in real time. Copyright © 2005 John Wiley & Sons, Ltd.     

Large‐scale real‐time hybrid simulation involving multiple experimental substructures and adaptive actuator delay compensation

Real‐time hybrid simulation provides a viable method to experimentally evaluate the performance of structural systems subjected to earthquakes. The structural system is divided into substructures, where part of the system is modeled by experimental substructures, whereas the remaining part is modeled analytically. The displacements in a real‐time hybrid simulation are imposed by servo‐hydraulic actuators to the experimental substructures. Actuator delay compensation has been shown by numerous researchers to vitally achieve reliable real‐time hybrid simulation results. Several studies have been performed on servo‐hydraulic actuator delay compensation involving single experimental substructure with single actuator. Research on real‐time hybrid simulation involving multiple experimental substructures, however, is limited. The effect of actuator delay during a real‐time hybrid simulation with multiple experimental substructures presents challenges. The restoring forces from experimental substructures may be coupled to two or more degrees of freedom (DOF) of the structural system, and the delay in each actuator must be adequately compensated. This paper first presents a stability analysis of actuator delay for real‐time hybrid simulation of a multiple‐DOF linear elastic structure to illustrate the effect of coupled DOFs on the stability of the simulation. An adaptive compensation method then proposed for the stable and accurate control of multiple actuators for a real‐time hybrid simulation. Real‐time hybrid simulation of a two‐story four‐bay steel moment‐resisting frame with large‐scale magneto‐rheological dampers in passive‐on mode subjected to the design basis earthquake is used to experimentally demonstrate the effectiveness of the compensation method in minimizing actuator delay in multiple experimental substructures. Copyright © 2011 John Wiley & Sons, Ltd.     

Visco‐acoustic prestack reverse‐time migration based on the time‐space domain adaptive high‐order finite‐difference method

It is important to include the viscous effect in seismic numerical modelling and seismic migration due to the ubiquitous viscosity in an actual subsurface medium. Prestack reverse‐time migration (RTM) is currently one of the most accurate methods for seismic imaging. One of the key steps of RTM is wavefield forward and backward extrapolation and how to solve the wave equation fast and accurately is the essence of this process. In this paper, we apply the time‐space domain dispersion‐relation‐based finite‐difference (FD) method for visco‐acoustic wave numerical modelling. Dispersion analysis and numerical modelling results demonstrate that the time‐space domain FD method has great accuracy and can effectively suppress numerical dispersion. Also, we use the time‐space domain FD method to solve the visco‐acoustic wave equation in wavefield extrapolation of RTM and apply the source‐normalized cross‐correlation imaging condition in migration. Improved imaging has been obtained in both synthetic and real data tests. The migration result of the visco‐acoustic wave RTM is clearer and more accurate than that of acoustic wave RTM. In addition, in the process of wavefield forward and backward extrapolation, we adopt adaptive variable‐length spatial operators to compute spatial derivatives to significantly decrease computing costs without reducing the accuracy of the numerical solution.     

An efficient algorithm for identifying pulse‐like ground motions based on significant velocity half‐cycles

Ground motions with strong velocity pulses are of particular interest to structural earthquake engineers because they have the potential to impose extreme seismic demands on structures. Accurate classification of records is essential in several earthquake engineering fields where pulse‐like ground motions should be distinguished from nonpulse‐like records, such as probabilistic seismic hazard analysis and seismic risk assessment of structures. This study proposes an effective method to identify pulse‐like ground motions having single, multiple, or irregular pulses. To effectively characterize the intrinsic pulse‐like features, the concept of an energy‐based significant velocity half‐cycle, which is visually identifiable, is first presented. Ground motions are classified into 6 categories according to the number of significant half‐cycles in the velocity time series. The pulse energy ratio is used as an indicator for quantitative identification, and then the energy threshold values for each type of ground motions are determined. Comprehensive comparisons of the proposed approach with 4 benchmark identification methods are conducted, and the results indicate that the methodology presented in this study can more accurately and efficiently distinguish pulse‐like and nonpulse‐like ground motions. Also presented are some insights into the reasons why many pulse‐like ground motions are not detected successfully by each of the benchmark methods.     

Studies on seismic reduction of story‐incresed buildings with friction layer and energy‐dissipated devices

A new type of energy‐dissipated structural system for existing buildings with story‐increased frames is presented and investigated in this paper. In this system the sliding‐friction layer between the lowest increased floor of the outer frame structure and the roof of the original building is applied, and energy‐dissipated dampers are used for the connections between the columns of the outer frame and each floor of the original building. A shaking table test is performed on the model of the system and the simplified structural model of this system is given. The theory of the non‐classical damping approach is introduced to the calculation analyses and compared with test results. The results show that friction and energy‐dissipated devices are very effective in reducing the seismic response and dissipating the input energy of the model structure. Finally, the design scheme and dynamic time‐history analyses of an existing engineering project are investigated to illustrate the application and advantages of the given method. Copyright © 2003 John Wiley & Sons, Ltd.     

Fast and accurate three‐dimensional controlled source electromagnetic modelling†

We present a fast approximate method for three‐dimensional low frequency controlled source electro‐magnetic modeling. We apply the method to a synthetic model in a typical marine controlled source electromagnetic scenario, where conductivity and permittivity are different from the known background medium. For 3D configurations, fast computational methods are relevant for both forward and inverse modelling studies. Since this problem involves a large number of unknowns, it has to be solved efficiently to obtain results in a timely manner, without compromising accuracy. For this reason, the Born approximation, extended Born approximation and iterative extended Born approximation are implemented and compared with the full solution of the conjugate gradient fast Fourier transformation method. These methods are based on an electric field domain integral equation formulation. It is shown here how well the iterative extended Born approximation method performs in terms of both accuracy and speed with different configurations and different source positions. The improved accuracy comes at virtually no additional computational cost. With the help of this method, it is now possible to perform sensitivity analysis using 3D modelling in a timely manner, which is vital for controlled source electromagnetic applications. For forward modeling the solution at the sea‐bottom is of interest, because that is where the receivers are usually located. For inverse modeling, the accuracy of the solution in the target zone is important to obtain reasonably accurate conductivity values from the inversion using this approximate solution method. Our modelling studies show that the iterative extended Born approximation method is fast and accurate for both forward and inverse modelling. Sensitivity analysis as a function of the source position and different reservoir sizes validate the accuracy of the iterative extended Born approximation.     

Application of wavelet theory to identify yielding in seismic response of bi‐linear structures

The time–frequency and the time‐scale analysis methods are used in this paper to identify the dynamic characteristics of non‐linear seismic response of structural systems with single degree of freedom (SDOF) and multiple degrees of freedom (MDOF). Based on the floor acceleration response time histories of bi‐linear SDOF and MDOF structures, the current study compares the results of system identification using the short‐time Fourier transform (STFT), continuous wavelet transform (CWT) and discrete wavelet transform (DWT) methods. The aim is to identify the frequency variations and the time at on‐set of yielding and unloading of a bi‐linear structural system during seismic response. The results demonstrate that the CWT method is better than the STFT method in both time and frequency resolutions, and that the DWT method is the best at detecting the time at on‐set of yielding and unloading. Combining the results of CWT and DWT methods therefore provides accurate information of both frequency variations and yielding time in non‐linear seismic response. To alleviate the problems associated with noise‐contaminated signals, e.g. seismic response data recorded on site, the study suggests that low‐pass filtering be carried out before applying the DWT method to decompose the signals into multiple levels of details. Copyright © 2001 John Wiley & Sons, Ltd.     

Evaluation of a recently proposed record selection and scaling procedure for low‐rise to mid‐rise reinforced concrete buildings and its use for probabilistic risk assessment studies

This paper evaluates a recent record selection and scaling procedure of the authors that can determine the probabilistic structural response of buildings behaving either in the elastic or post‐elastic range. This feature marks a significant strength on the procedure as the probabilistic structural response distribution conveys important information on probability‐based damage assessment. The paper presents case studies that show the utilization of the proposed record selection and scaling procedure as a tool for the estimation of damage states and derivation of site‐specific and region‐specific fragility functions. The method can be used to describe exceedance probabilities of damage limits under a certain target hazard level with known annual exceedance rate (via probabilistic seismic hazard assessment). Thus, the resulting fragility models can relate the seismicity of the region (or a site) with the resulting building performance in a more accurate manner. Under this context, this simple and computationally efficient record selection and scaling procedure can be benefitted significantly by probability‐based risk assessment methods that have started to be considered as indispensable for developing robust earthquake loss models. Copyright © 2013 John Wiley & Sons, Ltd.     

Applicability of pushover methods for the seismic analysis of single‐column bent viaducts

An overview of the applicability of a typical single‐mode pushover method (the N2 method) and two typical multi‐mode pushover methods (the modal pushover analysis (MPA) and incremental response spectrum analysis (IRSA) methods) for the analysis of single column bent viaducts in the transverse direction is presented. Previous research, which was limited to relatively short viaducts supported by few columns, has been extended to longer viaducts with more bents. The single‐mode N2 method is accurate enough for bridges where the effective modal mass of the fundamental mode is at least 80% of the total mass. The applicability of this method depends on (a) the ratio of the stiffness of the superstructure to that of the bents and (b) the strength of the bents. In short bridges with few columns, the accuracy of the N2 method increases as the seismic intensity increases, whereas in long viaducts (e.g. viaducts with lengths greater than 500 m) the method is in general less effective. In the case of the analyzed moderately irregular long viaducts, which are common in construction design practice, the MPA method performed well. For the analysis of bridges where the modes change significantly, depending on the seismic intensity, the IRSA method is in principle more appropriate, unless a viaduct is torsionally sensitive. In such cases, all simplified methods should be used with care. Copyright © 2008 John Wiley & Sons, Ltd.     

A comparison of snowmelt‐derived streamflow from temperature‐index and modified‐temperature‐index snow models

Reliable estimation of the volume and timing of snowmelt runoff is vital for water supply and flood forecasting in snow‐dominated regions. Snowmelt is often simulated using temperature‐index (TI) models due to their applicability in data‐sparse environments. Previous research has shown that a modified‐TI model, which uses a radiation‐derived proxy temperature instead of air temperature as its surrogate for available energy, can produce more accurate snow‐covered area (SCA) maps than a traditional TI model. However, it is unclear whether the improved SCA maps are associated with improved snow water equivalent (SWE) estimation across the watershed or improved snowmelt‐derived streamflow simulation. This paper evaluates whether a modified‐TI model produces better streamflow estimates than a TI model when they are used within a fully distributed hydrologic model. It further evaluates the performance of the two models when they are calibrated using either point SWE measurements or SCA maps. The Senator Beck Basin in Colorado is used as the study site because its surface is largely bedrock, which reduces the role of infiltration and emphasizes the role of the SWE pattern on streamflow generation. Streamflow is simulated using both models for 6 years. The modified‐TI model produces more accurate streamflow estimates (including flow volume and peak flow rate) than the TI model, likely because the modified‐TI model better reproduces the SWE pattern across the watershed. Both models also produce better performance when calibrated with SCA maps instead of point SWE data, likely because the SCA maps better constrain the space‐time pattern of SWE.     

Alternative non‐linear demand estimation methods for probability‐based seismic assessments

Alternative non‐linear dynamic analysis procedures, using real ground motion records, can be used to make probability‐based seismic assessments. These procedures can be used both to obtain parameter estimates for specific probabilistic assessment criteria such as demand and capacity factored design and also to make direct probabilistic performance assessments using numerical methods. Multiple‐stripe analysis is a non‐linear dynamic analysis method that can be used for performance‐based assessments for a wide range of ground motion intensities and multiple performance objectives from onset of damage through global collapse. Alternatively, the amount of analysis effort needed in the performance assessments can be reduced by performing the structural analyses and estimating the main parameters in the region of ground motion intensity levels of interest. In particular, single‐stripe and double‐stripe analysis can provide local probabilistic demand assessments using minimal number of structural analyses (around 20 to 40). As a case study, the displacement‐based seismic performance of an older reinforced concrete frame structure, which is known to have suffered shear failure in its columns during the 1994 Northridge Earthquake, is evaluated. Copyright © 2008 John Wiley & Sons, Ltd.     

An automated cross‐validation method to assess seismic time‐to‐depth conversion accuracy: a case study on the Cooper and Eromanga basins,Australia

Selecting a seismic time‐to‐depth conversion method can be a subjective choice that is made by geophysicists, and is particularly difficult if the accuracy of these methods is unknown. This study presents an automated statistical approach for assessing seismic time‐to‐depth conversion accuracy by integrating the cross‐validation method with four commonly used seismic time‐to‐depth conversion methods. To showcase this automated approach, we use a regional dataset from the Cooper and Eromanga basins, Australia, consisting of 13 three‐dimensional (3D) seismic surveys, 73 two‐way‐time surface grids and 729 wells. Approximately 10,000 error values (predicted depth vs. measured well depth) and associated variables were calculated. The average velocity method was the most accurate overall (7.6 m mean error); however, the most accurate method and the expected error changed by several metres depending on the combination and value of the most significant variables. Cluster analysis tested the significance of the associated variables to find that the seismic survey location (potentially related to local geology (i.e. sedimentology, structural geology, cementation, pore pressure, etc.), processing workflow, or seismic vintage), formation (potentially associated with reduced signal‐to‐noise with increasing depth or the changes in lithology), distance to the nearest well control, and the spatial location of the predicted well relative to the existing well data envelope had the largest impact on accuracy. Importantly, the effect of these significant variables on accuracy were found to be more important than choosing between the four methods, highlighting the importance of better understanding seismic time‐to‐depth conversions, which can be achieved by applying this automated cross‐validation method.     

Evaluation and generalization of temperature‐based methods for calculating evaporation

Seven temperature‐based equations, each representing a typical form, were evaluated and compared for determining evaporation at two climatological stations (Rawson Lake and Atikokan) in north‐western Ontario, Canada. The comparison was first made using the original constant values involved in each equation, and then using the recalibrated constant values. The results show that when the original constant values were used, larger biases existed for most of the equations for both stations. When recalibrated constant values were substituted for the original constant values, six of the seven equations improved for both stations. Using locally calibrated parameter values, all seven equations worked well for determining mean seasonal evaporation values. For monthly evaporation values, the modified Blaney–Criddle method produced least error for all months for both stations, followed by the Hargreaves and Thornthwaite methods. The Linacre, Kharrufa and Hamon methods showed a significant bias in September for both stations. With properly determined constant values, the modified Blaney–Criddle, the Hargreaves and Thornthwaite methods can be recommended for estimating evaporation in the study region, as far as temperature‐based methods are concerned. Copyright © 2001 John Wiley & Sons, Ltd.     

Piecewise‐continuous velocity inversion

We suggest a new method to determine the piecewise‐continuous vertical distribution of instantaneous velocities within sediment layers, using different order time‐domain effective velocities on their top and bottom points. We demonstrate our method using a synthetic model that consists of different compacted sediment layers characterized by monotonously increasing velocity, combined with hard rock layers, such as salt or basalt, characterized by constant fast velocities, and low velocity layers, such as gas pockets. We first show that, by using only the root‐mean‐square velocities and the corresponding vertical travel times (computed from the original instantaneous velocity in depth) as input for a Dix‐type inversion, many different vertical distributions of the instantaneous velocities can be obtained (inverted). Some geological constraints, such as limiting the values of the inverted vertical velocity gradients, should be applied in order to obtain more geologically plausible velocity profiles. In order to limit the non‐uniqueness of the inverted velocities, additional information should be added. We have derived three different inversion solutions that yield the correct instantaneous velocity, avoiding any a priori geological constraints. The additional data at the interface points contain either the average velocities (or depths) or the fourth‐order average velocities, or both. Practically, average velocities can be obtained from nearby wells, whereas the fourth‐order average velocity can be estimated from the quartic moveout term during velocity analysis. Along with the three different types of input, we consider two types of vertical velocity models within each interval: distribution with a constant velocity gradient and an exponential asymptotically bounded velocity model, which is in particular important for modelling thick layers. It has been shown that, in the case of thin intervals, both models lead to similar results. The method allows us to establish the instantaneous velocities at the top and bottom interfaces, where the velocity profile inside the intervals is given by either the linear or the exponential asymptotically bounded velocity models. Since the velocity parameters of each interval are independently inverted, discontinuities of the instantaneous velocity at the interfaces occur naturally. The improved accuracy of the inverted instantaneous velocities is particularly important for accurate time‐to‐depth conversion.     

Residual structural capacity evaluation of steel moment‐resisting frames with dynamic‐strain‐based model updating method

This paper presents a method for evaluating the residual structural capacity of earthquake‐affected steel structures. The method first quantifies the damage severity of a beam by computing the dynamic‐strain‐based damage index. Next, the model used to analyze the structure is updated based on the damage index, to reflect the observed damage conditions. The residual structural capacity is then estimated in terms of changes in stiffness and strength, which can be applied by structural engineers, via a nonlinear static analysis of the updated model. The main contributions of this paper are in performance evaluation of the dynamic‐strain‐based damage index for seismically induced damage using a newly developed substructure testing environment, consideration of various damage patterns in composite beams, and extension of a local damage evaluation technique to a residual capacity estimation procedure by incorporating the model‐updating technique. In laboratory testing, the specimens were damaged quasi‐statically, and vibration tests were conducted as the damage proceeded. First, a bare steel beam–column connection was tested, and then a similar one with a floor slab was used for a more realistic case. The estimated residual structural capacities for these specimens were compared with the static test results. The results verified that the proposed method can provide fine estimates of the stiffness and strength deteriorations within 10% for the specimen without the floor slab and within 30% for that with the floor slab. Copyright © 2017 John Wiley & Sons, Ltd.     

A modal pushover analysis procedure to estimate seismic demands for unsymmetric‐plan buildings

An Erratum has been published for this article in Earthquake Engng. Struct. Dyn. 2004; 33:1429. Based on structural dynamics theory, the modal pushover analysis (MPA) procedure retains the conceptual simplicity of current procedures with invariant force distribution, now common in structural engineering practice. The MPA procedure for estimating seismic demands is extended to unsymmetric‐plan buildings. In the MPA procedure, the seismic demand due to individual terms in the modal expansion of the effective earthquake forces is determined by non‐linear static analysis using the inertia force distribution for each mode, which for unsymmetric buildings includes two lateral forces and torque at each floor level. These ‘modal’ demands due to the first few terms of the modal expansion are then combined by the CQC rule to obtain an estimate of the total seismic demand for inelastic systems. When applied to elastic systems, the MPA procedure is equivalent to standard response spectrum analysis (RSA). The MPA estimates of seismic demand for torsionally‐stiff and torsionally‐flexible unsymmetric systems are shown to be similarly accurate as they are for the symmetric building; however, the results deteriorate for a torsionally‐similarly‐stiff unsymmetric‐plan system and the ground motion considered because (a) elastic modes are strongly coupled, and (b) roof displacement is underestimated by the CQC modal combination rule (which would also limit accuracy of RSA for linearly elastic systems). Copyright © 2004 John Wiley & Sons, Ltd.