Non‐stationary seismic response of tanks with soil interaction by wavelets

A wavelet‐based random vibration theory has been developed for the non‐stationary seismic response of liquid storage tanks including soil interaction. The ground motion process has been characterized via estimates of statistical functionals of wavelet coefficients obtained from a single time history of ground accelerations. The tank–liquid–soil system has been modelled as a two‐degree‐of‐freedom (2‐DOF) system. The wavelet domain equations have been formulated and the wavelet coefficients of the required response state are obtained by solving two linear simultaneous algebraic equations. The explicit expression for the instantaneous power spectral density function (PSDF) in terms of the functionals of the input wavelet coefficients has been obtained. The moments of this PSDF are used to estimate the expected pseudo‐spectral acceleration (PSA) response of the tank. Parametric variations are carried out to study the effects of tank height, foundation natural frequency, shear wave velocity of soil and ratio of the mass of tank (including liquid) to the mass of foundation on the PSA responses of tanks. Copyright © 2001 John Wiley & Sons, Ltd.     

Wavelet‐based characterization of design ground motions

With the recent emergence of wavelet‐based procedures for stochastic analyses of linear and non‐linear structural systems subjected to earthquake ground motions, it has become necessary that seismic ground motion processes are characterized through statistical functionals of wavelet coefficients. While direct characterization in terms of earthquake and site parameters may have to wait for a few more years due to the complexity of the problem, this study attempts such characterization through commonly available Fourier and response spectra for design earthquake motions. Two approaches have been proposed for obtaining the spectrum‐compatible wavelet functionals, one for input Fourier spectrum and another for input response spectrum, such that the total number of input data points are 30–35% of those required for a time‐history analysis. The proposed methods provide for simulating ‘desired non‐stationary characteristics’ consistent with those in a recorded accelerogram. Numerical studies have been performed to illustrate the proposed approaches. Further, the wavelet functionals compatible with a USNRC spectrum in the case of 35 recorded motions of similar strong motion durations have been used to obtain the strength reduction factor spectra for elasto‐plastic oscillators and to show that about ±20% variation may be assumed from mean to 5 and 95% confidence levels due to uncertainty in the non‐stationary characteristics of the ground motion process. Copyright © 2002 John Wiley & Sons, Ltd.     

Identification of modal parameters of non‐classically damped linear structures under multi‐component earthquake loading

A method for parametric system identification of classically damped linear system in frequency domain is adopted and extended for non‐classically damped linear systems subjected up to six components of earthquake ground motions. This method is able to work in multi‐input/multi‐output (MIMO) case. The response of a two‐degree‐of‐freedom model with non‐classical damping, excited by one‐component earthquake ground motion, is simulated and used to verify the proposed system identification method in the single‐input/multi‐output case. Also, the records of a 10 storey real building during the Northridge earthquake is used to verify the proposed system identification method in the MIMO case. In this case, at first, a single‐input/multi‐output assumption is considered for the system and modal parameters are identified, then other components of earthquake ground motions are added, respectively, and the modal parameters are identified again. This procedure is repeated until all four components of earthquake ground motions which are measured at the base level of the building are included in the identification process. The results of identification of real building show that consideration of non‐classical damping and inclusion of the multi‐components effect of earthquake ground motions can improve the least‐squares match between the finite Fourier transforms of recorded and calculated acceleration responses. Copyright © 2006 John Wiley & Sons, Ltd.     

Simulations of non‐stationary frequency content and its importance to seismic assessment of structures

To verify the importance of the non‐stationary frequency characteristic of seismic ground motion, a joint time–frequency analysis technique of time signals, called chirplet‐based signal approximation, is developed to extract the non‐stationary frequency information from the recorded data. The chirplet‐based signal approximation is clear in concept, similar to Fourier Transform in mathematical expressions but with different base functions. Case studies show that the chirplet‐based signal approximation can represent the joint time–frequency variation of seismic ground motion quite well. Both the random models of uniform modulating process and evolutionary process are employed to generate artificial seismic waves. The joint time–frequency modulating function in the random model of evolutionary process is determined by chirplet‐based signal approximation. Finally, non‐linear response analysis of a SODF system and a frame structure is performed based on the generated artificial seismic waves. The results show that the non‐stationary frequency characteristic of seismic ground motion can significantly change the non‐linear response characteristics of structures, particularly when a structure goes into collapse phase under seismic action. It is concluded that non‐stationary frequency characteristic of seismic ground motion should be considered for the assessment of seismic capacity of structures. Copyright © 2002 John Wiley & Sons, Ltd.     

Estimate of input energy for elasto‐plastic SDOF systems during earthquakes based on discrete wavelet coefficients

The response of an elasto‐plastic single degree of freedom (SDOF) system to ground motion is estimated based on wavelet coefficients calculated by discrete wavelet transform. Wavelet coefficients represent both the time and frequency characteristics of input ground motion, and thus can be considered to be directly related to the dynamic response of a non‐linear system. This relationship between the energy input into an elastic SDOF system and wavelet coefficients is derived based on the assumption that wavelets deliver energy to the structure instantaneously and the quantity of energy is constant regardless of yielding. These assumptions are shown to be valid when the natural period of the system is in the predominant period range of the wavelet, the most common scenario for real structures, through dynamic response analysis of a single wavelet. The wavelet‐based estimation of elastic and plastic energy transferred by earthquake ground motion is thus shown to be in good agreement with the dynamic response analysis when the natural period is in the predominant range of the input. Copyright © 2005 John Wiley & Sons, Ltd.     

Free‐vibration analysis of a three‐dimensional soil–structure system

A procedure which involves a non‐linear eigenvalue problem and is based on the substructure method is proposed for the free‐vibration analysis of a soil–structure system. In this procedure, the structure is modelled by the standard finite element method, while the unbounded soil is modelled by the scaled boundary finite element method. The fundamental frequency, and the corresponding radiation damping ratio as well as the modal shape are obtained by using inverse iteration. The free vibration of a dam–foundation system, a hemispherical cavity and a hemispherical deposit are analysed in detail. The numerical results are compared with available results and are also verified by the Fourier transform of the impulsive response calculated in the time domain by the three‐dimensional soil–structure–wave interaction analysis procedure proposed in our previous paper. The fundamental frequency obtained by the present procedure is very close to that obtained by Touhei and Ohmachi, but the damping ratio and the imaginary part of modal shape are significantly different due to the different definition of damping ratio. This study shows that although the classical mode‐superposition method is not applicable to a soil–structure system due to the frequency dependence of the radiation damping, it is still of interest in earthquake engineering to evaluate the fundamental frequency and the corresponding radiation damping ratio of the soil–structure system. Copyright © 2001 John Wiley & Sons, Ltd.     

Wavelet‐based non‐stationary seismic rocking response of flexibly supported tanks

The non‐stationary rocking response of liquid storage tanks under seismic base excitations including soil interaction has been developed based on the wavelet domain random vibration theory. The ground motion has been characterized through statistical functionals of wavelet coefficients of the ground acceleration history. The tank–liquid–foundation system is modelled as a multi‐degree‐of‐freedom (MDOF) system with both lateral and rocking motions of vibration of the foundation. The impulsive and convective modes of vibration of the liquid in the tank have been considered. The wavelet domain coupled dynamic equations are formulated and then solved to get the expressions of instantaneous power spectral density function (PSDF) in terms of functionals of input wavelet coefficients. The moments of the instantaneous PSDF are used to obtain the stochastic responses of the tank in the form of coefficients of hydrodynamic pressure, base shear and overturning base moment for the largest expected peak responses. Parametric variations are carried out to study the effects of various governing parameters like height of liquid in the tank, height–radius ratio of the tank, ratio of total liquid mass to mass of foundation, and shear wave velocity in the soil medium, on the responses of the tank. Copyright © 2003 John Wiley & Sons, Ltd.     

Non-stationary seismic response of MDOF systems by wavelet transform

A wavelet-based formulation has been presented in this paper for the stochastic analysis of a linear multi-degree-of-freedom (MDOF) classically damped system subjected to earthquake ground motion. The ground motion has been modelled as a non-stationary process (both in amplitude and frequency) using wavelets. Closed-form expressions of the moments of the instantaneous Power Spectral Density Function (PSDF) of the response have been derived and used to predict the statistics of the response peak of any desired order. For illustration of the formulation, an example torsionally coupled multistoried building has been considered along with the twenty synthetically generated time-histories corresponding to an example ground motion process. © 1997 John Wiley & Sons, Ltd.     

Asymmetric one‐storey elastic systems with non‐linear viscous and viscoelastic dampers: Earthquake response

Investigated are earthquake responses of one‐way symmetric‐plan, one‐storey systems with non‐linear fluid viscous dampers (FVDs) attached in series to a linear brace (i.e. Chevron or inverted V‐shape braces).Thus, the non‐linear damper is viscous when the brace is considered rigid or viscoelastic (VE) when the brace is flexible. The energy dissipation capacity of a non‐linear FVD is characterized by an amplitude‐dependent damping ratio for an energy‐equivalent linear FVD, which is determined assuming the damper undergoes harmonic motion. Although this formulation is shown to be advantageous for single‐degree‐of‐freedom (SDF) systems, it is difficult to extend its application to multi‐degree‐of‐freedom (MDF) systems for two reasons: (1) the assumption that dampers undergo harmonic motion in parameterizing the non‐linear damper is not valid for its earthquake‐induced motion of an MDF system; and (2) ensuring simultaneous convergence of all unknown amplitudes of dampers is difficult in an iterative solution of the non‐linear system. To date, these limitations have precluded the parametric study of the dynamics of MDF systems with non‐linear viscous or VE dampers. However, they are overcome in this investigation using concepts of modal analysis because the system is weakly non‐linear due to supplemental damping. It is found that structural response is only weakly affected by damper non‐linearity and is increased by a small amount due to bracing flexibility. Thus, the effectiveness of supplemental damping in reducing structural responses and its dependence on the planwise distribution of non‐linear VE dampers were found to be similar to that of linear FVDs documented elsewhere. As expected, non‐linear viscous and VE dampers achieve essentially the same reduction in response but with much smaller damper force compared to linear dampers. Finally, the findings in this investigation indicate that the earthquake response of the asymmetric systems with non‐linear viscous or VE dampers can be estimated with sufficient accuracy for design applications by analysing the same asymmetric systems with all non‐linear dampers replaced by energy‐equivalent linear viscous dampers. Copyright © 2003 John Wiley & Sons, Ltd.     

Stochastic model for simulation of near‐fault ground motions

A parameterized stochastic model of near‐fault ground motion in two orthogonal horizontal directions is developed. The major characteristics of recorded near‐fault ground motions are represented. These include near‐fault effects of directivity and fling step; temporal and spectral non‐stationarity; intensity, duration, and frequency content characteristics; directionality of components; and the natural variability of ground motions. Not all near‐fault ground motions contain a forward directivity pulse, even when the conditions for such a pulse are favorable. The proposed model accounts for both pulse‐like and non‐pulse‐like cases. The model is fitted to recorded near‐fault ground motions by matching important characteristics, thus generating an ‘observed’ set of model parameters for different earthquake source and site characteristics. A method to generate and post‐process synthetic motions for specified model parameters is also presented. Synthetic ground motion time series are generated using fitted parameter values. They are compared with corresponding recorded motions to validate the proposed model and simulation procedure. The use of synthetic motions in addition to or in place of recorded motions is desirable in performance‐based earthquake engineering applications, particularly when recorded motions are scarce or when they are unavailable for a specified design scenario. Copyright © 2016 John Wiley & Sons, Ltd.     

A method of calculating the non‐linear seismic response of a building braced with viscoelastic dampers

Buildings are continually subject to dynamic loads, such as wind load, seismic ground motion, and even the load from internal utility machines. The recent trend of constructing more flexible high‐rise buildings underscores the importance of including viscoelastic dampers in building designs. Viscoelastic dampers are used to control the dynamic response of a building. If the seismic design is based only on the linear response spectrum, considerable error may occur when calculating the seismic response of a building; rubber viscoelastic dampers show non‐linear hysteretic damping that is quite different from viscous damping. This study generated a non‐linear response spectrum using a non‐linear oscillator model to simulate a building with viscoelastic dampers installed. The parameters used in the non‐linear damper model were obtained experimentally from dynamic loading tests. The results show that viscoelastic dampers effectively reduce the seismic displacement response of a structure, but transmit more seismic force to the structure, which essentially increases its seismic acceleration response. Copyright © 2003 John Wiley & Sons, Ltd.     

Asymmetric one‐storey elastic systems with non‐linear viscous and viscoelastic dampers: Simplified analysis and supplemental damping system design

Investigated is the accuracy in estimating the response of asymmetric one‐storey systems with non‐linear viscoelastic (VE) dampers by analysing the corresponding linear viscous system wherein all non‐linear VE dampers are replaced by their energy‐equivalent linear viscous dampers. The response of the corresponding linear viscous system is determined by response history analysis (RHA) and by response spectrum analysis (RSA) extended for non‐classically damped systems. The flexible and stiff edge deformations and plan rotation of the corresponding linear viscous system determined by the extended RSA procedure is shown to be sufficiently accurate for design applications with errors generally between 10 and 20%. Although similar accuracy is also shown for the ‘pseudo‐velocity’ of non‐linear VE dampers, the peak force of the non‐linear VE damper cannot be estimated directly from the peak damper force of the corresponding linear viscous system. A simple correction for damper force is proposed and shown to be accurate (with errors not exceeding 15%). For practical applications, an iterative linear analysis procedure is developed for determining the amplitude‐ and frequency‐dependent supplemental damping properties of the corresponding linear viscous system and for estimating the response of asymmetric one‐storey systems with non‐linear VE dampers from the earthquake design (or response) spectrum. Finally, a procedure is developed for designing non‐linear supplemental damping systems that satisfy given design criteria for a given design spectrum. Copyright © 2003 John Wiley & Sons, Ltd.     

Response of SDOF system to non‐stationary earthquake excitation

Earthquake excitation is often modelled by non‐stationary random process (i.e. uniformly modulated broad‐band excitation) for analysis of structural safety subjected to seismic hazards. In this paper, the non‐stationary response of a single‐degree‐of‐freedom (SDOF) system to non‐stationary earthquake motion is investigated for different shapes of modulating functions. The evolutionary power‐spectral density function (PSDF) of the displacement of the SDOF system is obtained using the time‐varying frequency response function and the PSDF of the earthquake excitation. The close form expressions for time‐varying frequency response function are derived for different shapes of the modulating functions. In order to study the effects of the shape of the modulating function, a comparison of the non‐stationary earthquake response of the SDOF system is also made for different modulating functions having the same energy content. Copyright © 2004 John Wiley & Sons, Ltd.     

Seismic response of a cable‐stayed bridge deck under multi‐component non‐stationary random ground motion

A Markov method of analysis is presented for obtaining the seismic response of cable‐stayed bridges to non‐stationary random ground motion. A uniformly modulated non‐stationary model of the random ground motion is assumed which is specified by the evolutionary r.m.s. ground acceleration. Both vertical and horizontal components of the motion are considered to act simultaneously at the bridge supports. The analysis duly takes into account the angle of incidence of the earthquake, the spatial correlation of ground motion and the quasi‐static excitation. A cable‐stayed bridge is analysed under a set of parametric variations in order to study the non‐stationary response of the bridge. The results of the numerical study indicate that (i) frequency domain spectral analysis with peak r.m.s. acceleration as input could provide more r.m.s. response than the peak r.m.s. response obtained by the non‐stationary analysis; (ii) the longitudinal component of the ground motion significantly influences the vertical vibration of the bridge; and (iii) the angle of incidence of the earthquake has considerable influence on the deck response. Copyright © 2003 John Wiley & Sons, Ltd.     

Estimating non‐linear site response by single period approximation

Site effects characterize the filtering mechanisms within the soil sedimentary layers overlying bedrock. In regions of high seismicity such as California where strong motion records are relatively abundant, site coefficients can be developed by regression of recorded ground shaking parameters. In regions of low‐to‐moderate seismicity or of high seismicity but with a paucity of recorded strong motion data, such empirical models cannot be obtained in the same way. This study describes the theoretical development of a simple, rational manual procedure to calculate site coefficients, based on a single period approximation (SPA), and to construct displacement response spectra (RSD) for soil sites. The proposed simplified model, which takes into account the non‐linear behaviour of soil that is dependent on the level of shaking, impedance contrast at the soil–bedrock interface and the plasticity of soil material, has been verified by comparison with results obtained from non‐linear shear wave analyses and data recorded during the 1994 Northridge earthquake. The proposed model is believed to be a convenient tool for calculating non‐linear site responses and constructing site‐specific response spectra, which has the potential of being incorporated into code provisions. Copyright © 2006 John Wiley & Sons, Ltd.     

Operator‐splitting method for real‐time substructure testing

It has been shown that the operator‐splitting method (OSM) provides explicit and unconditionally stable solutions for quasi‐static pseudo‐dynamic substructure testing. However, the OSM provides only an explicit target displacement but not an explicit target velocity, so that it is essentially an implicit method for real‐time substructure testing (RST) when the velocity‐dependent restoring force is considered. This paper proposes a target velocity formulation based on the forward difference of the predicted displacements so as to render the OSM explicit for RST. The stability and accuracy of the resulting OSM‐RST algorithm are investigated. It is shown that the OSM‐RST is unconditionally stable so long as the non‐linear stiffness and damping are of the softening type (i.e. the tangent stiffness and damping never exceed the initial values). The stability of the OSM‐RST for structures with infinite tangent damping coefficient or stiffness is also proved, and the stability of the method for MDOF structures with a non‐classical damping matrix is demonstrated by an energy criterion. The effects of actuator delay and compensation are analysed based on the bilinear approximation of the actuator step response. Experiments on damped SDOF and MDOF structures verify that the stability of the OSM‐RST is preserved when the experimental substructure generates velocity‐dependent reaction forces, whereas the stability of real‐time substructure tests based on the central difference method is worsened by the damping of the specimen. Copyright © 2005 John Wiley & Sons, Ltd.     

Estimating the seismic displacement of friction pendulum isolators based on non‐linear response history analysis

A procedure for developing equations that estimate the isolator displacement due to strong ground motion is applied to buildings isolated with the friction pendulum system. The resulting design equations, based on rigorous non‐linear analysis, offer an alternative to the iterative equivalent‐linear methods used by current U.S. building codes. The governing equations of the system are reduced to a form such that the median normalized displacement of the system due to an ensemble of ground motions is found to depend on only the isolation period—a function of the curvature of the isolator—and the friction force at incipient slip normalized by peak ground velocity. The normalization is effective in minimizing the dispersion of the normalized displacement for an ensemble of ground motions, implying that the median normalized displacement is a reliable estimate of response. The design equations reflect the significant (20 to 38%) increase in displacement when the excitation includes two lateral components of ground motion instead of just one component. Equivalent‐linear methods are shown to underestimate by up to 30% the exact median displacement determined by non‐linear response history analysis for one component of ground motion, and building codes include at most a 4.4% increase for a second component. Copyright © 2003 John Wiley & Sons, Ltd.     

Semi‐empirical model for site effects on acceleration time histories at soft‐soil sites. Part 2: calibration

A previously developed simplified model of ground motion amplification is applied to the simulation of acceleration time histories at several soft‐soil sites in the Valley of Mexico, on the basis of the corresponding records on firm ground. The main objective is to assess the ability of the model to reproduce characteristics such as effective duration, frequency content and instantaneous intensity. The model is based on the identification of a number of parameters that characterize the complex firm‐ground to soft‐soil transfer function, and on the adjustment of these parameters in order to account for non‐linear soil behavior. Once the adjusted model parameters are introduced, the statistical properties of the simulated and the recorded ground motions agree reasonably well. For the sites and for the seismic events considered in this study, it is concluded that non‐linear soil behavior may have a significant effect on the amplification of ground motion. The non‐linear soil behavior significantly affects the effective ground motion duration for the components with the higher intensities, but it does not have any noticeable influence on the lengthening of the dominant ground period. Copyright © 2004 John Wiley & Sons, Ltd.     

A nonlinear data‐driven model for synthetic generation of annual streamflows

A hybrid model that blends two non‐linear data‐driven models, i.e. an artificial neural network (ANN) and a moving block bootstrap (MBB), is proposed for modelling annual streamflows of rivers that exhibit complex dependence. In the proposed model, the annual streamflows are modelled initially using a radial basis function ANN model. The residuals extracted from the neural network model are resampled using the non‐parametric resampling technique MBB to obtain innovations, which are then added back to the ANN‐modelled flows to generate synthetic replicates. The model has been applied to three annual streamflow records with variable record length, selected from different geographic regions, namely Africa, USA and former USSR. The performance of the proposed ANN‐based non‐linear hybrid model has been compared with that of the linear parametric hybrid model. The results from the case studies indicate that the proposed ANN‐based hybrid model (ANNHM) is able to reproduce the skewness present in the streamflows better compared to the linear parametric‐based hybrid model (LPHM), owing to the effective capturing of the non‐linearities. Moreover, the ANNHM, being a completely data‐driven model, reproduces the features of the marginal distribution more closely than the LPHM, but offers less smoothing and no extrapolation value. It is observed that even though the preservation of the linear dependence structure by the ANNHM is inferior to the LPHM, the effective blending of the two non‐linear models helps the ANNHM to predict the drought and the storage characteristics efficiently. Copyright © 2007 John Wiley & Sons, Ltd.     

Influence of time‐varying frequency content in earthquake ground motions on seismic response of linear elastic systems

Earthquake ground motion records are nonstationary in both amplitude and frequency content. However, the latter nonstationarity is typically neglected mainly for the sake of mathematical simplicity. To study the stochastic effects of the time‐varying frequency content of earthquake ground motions on the seismic response of structural systems, a pair of closely related stochastic ground motion models is adopted here. The first model (referred to as ground motion model I) corresponds to a fully nonstationary stochastic earthquake ground motion model previously developed by the authors. The second model (referred to as ground motion model II) is nonstationary in amplitude only and is derived from the first model. Ground motion models I and II have the same mean‐square function and global frequency content but different features of time variation in the frequency content, in that no time variation of the frequency content exists in ground motion model II. New explicit closed‐form solutions are derived for the response of linear elastic SDOF and MDOF systems subjected to stochastic ground motion model II. New analytical solutions for the evolutionary cross‐correlation and cross‐PSD functions between the ground motion input and the structural response are also derived for linear systems subjected to ground motion model I. Comparative analytical results are presented to quantify the effects of the time‐varying frequency content of earthquake ground motions on the structural response of linear elastic systems. It is found that the time‐varying frequency content in the seismic input can have significant effects on the stochastic properties of system response. Copyright © 2016 John Wiley & Sons, Ltd.