Critical response of structures to multicomponent earthquake excitation

This paper aims to develop an improved understanding of the critical response of structures to multicomponent seismic motion characterized by three uncorrelated components that are defined along its principal axes: two horizontal and the vertical component. An explicit formula, convenient for code applications, has been derived to calculate the critical value of structural response to the two principal horizontal components acting along any incident angle with respect to the structural axes, and the vertical component of ground motion. The critical response is defined as the largest value of response for all possible incident angles. The ratio r_cr/r_srss between the critical value of response and the SRSS response—corresponding to the principal components of ground acceleration applied along the structure axes—is shown to depend on three dimensionless parameters: the spectrum intensity ratio γ between the two principal components of horizontal ground motion characterized by design spectra A(T_n) and γA(T_n); the correlation coefficient α of responses r_x and r_y due to design spectrum A(T_n) applied in the x‐ and y‐directions, respectively; and β = r_y/r_x. It is demonstrated that the ratio r_cr/r_srss is bounded by 1 and . Thus the largest value of the ratio is , 1.26, 1.13 and 1.08 for γ = 0, 0.5, 0.75 and 0.85, respectively. This implies that the critical response never exceeds times the result of the SRSS analysis, and this ratio is about 1.13 for typical values of γ, say 0.75. The correlation coefficient α depends on the structural properties but is always bounded between −1 and 1. For a fixed value of γ, the ratio r_cr/r_srss is largest if β = 1 and α = ±1. The parametric variations presented for one‐storey buildings indicate that this condition can be satisfied by axial forces in columns of symmetric‐plan buildings or can be approximated by lateral displacements in resisting elements of unsymmetrical‐plan buildings. Copyright © 2000 John Wiley & Sons, Ltd.     

Seismic motion and response prediction alternatives

Statistical methods are available which predict the maximum response of simple oscillators given the peak acceleration (A_p), peak velocity (V_p) or peak displacement (D_p) of seismic ground motions. An alternative parameter, namely an ordinate (or ordinates) of the Fourier amplitude spectrum of ground motion acceleration, FS(f), may in fact be a preferred predictor of peak response, especially in a frequency range close to f. Other statistical methods (attenuation laws) use distance R and other parameters such as magnitude (M), Modified Mercalli epicentral Intensity (I_o) and Modified Mercalli site Intensity (MMI or I_s) to predict spectral velocity (S_v(f)), etc. In using such approaches, it is most desirable to know the total uncertainty in the predicted peak response of the system given the starting parameter values. An extensive strong motion data set is used to study these questions, The most direct prediction models are found to be preferable (have lower prediction dispersion) but data may not be available in all regions to permit their use.     

Prediction of spatially distributed seismic demands in specific structures: Structural response to loss estimation

A companion paper has investigated the effects of intensity measure (IM) selection in the prediction of spatially distributed response in a multi‐degree‐of‐freedom structure. This paper extends from structural response prediction to performance assessment metrics such as probability of structural collapse; probability of exceeding a specified level of demand or direct repair cost; and the distribution of direct repair loss for a given level of ground motion. In addition, a method is proposed to account for the effect of varying seismological properties of ground motions on seismic demand that does not require different ground motion records to be used for each intensity level. Results illustrate that the conventional IM, spectral displacement at the first mode, S_de(T₁), produces higher risk estimates than alternative velocity‐based IM's, namely spectrum intensity, SI, and peak ground velocity, PGV, because of its high uncertainty in ground motion prediction and poor efficiency in predicting peak acceleration demands. Copyright © 2009 John Wiley & Sons, Ltd.     

黄河阶地特征与地震动参数的相关性研究

在分析已有资料的基础上划分了兰州盆地与建设工程分布密切相关的T₀、T₁、T₂、T₃和T₄级黄河阶地,建立201个土层地震反应模型。通过一维等效线性化计算和反应谱分析,得出兰州盆地沉积阶地50年超越概率10%地表地震动参数,分析阶地高度和v_S≤500 m/s覆盖层厚度特征与地震动参数峰值加速度Am和加速度反应谱特征周期T_g的相关性。表明兰州盆地T₀~T₂阶地覆盖层厚度与50年超越概率10%A_m呈正相关,T₃及以上阶地覆盖层厚度对A_m增大有明显的减小作用。T_g值随T₀~T₃阶地覆盖层厚度的增加而变大,当覆盖层厚度进一步变大,T_g值不再同步增大,阶地覆盖层厚度对T_g的影响是有限的,阶地海拔高度与地表50年超越概率10%地震动参数没有关系。     

Hyporheic exchange mechanism driven by flood wave

To study the hyporheic exchange driven by a single peak flood-induced water level fluctuation (i.e., flood wave), a method combining numerical simulation with theoretical derivation was proposed based on the Inbuk Stream, Korea, where flooding occurs frequently. The hyporheic exchanges induced by different flood waves were investigated by varying amplitude (A), duration (T), wave type parameter (r), and rising duration (t_p), which were adopted from the real-time stream stage fluctuations. Additionally, the idea of constant upstream flood volume (CUFV) condition for flood waves was put forward, and the effects of “Botan” (T/A) and peak number (N) on hyporheic exchange were studied. The results showed that the hyporheic exchange flux (q) was controlled by the water level h (sine-type) and its change rate v (cosine-type), and was proportional to the polynomial of them q ∝ (ω∙ h + v), where ω is the angular frequency of the flood wave. Based on this mechanism, the influence principles on hyporheic exchanges of the typical flood wave parameters (A, T, r and t_p) as well as T/A and N under CUFV condition were clarified. The main characteristic variables of hyporheic exchange, which were maximum aquifer storage and residence time, were positively correlated. They also had positive relations to the integral of the flood wave over time, which increased when the wave became higher, wider, rounder and less skewed. However, when CUFV condition was imposed, the residence time was positively correlated with T/A, whereas the maximum aquifer storage was negatively correlated with T/A. With the increase in N, water exchanged more frequently and some water returned to the stream early, leading to the slight decrease in maximum aquifer storage and residence time. These findings enriched the theory of hyporheic exchange driven by surface water fluctuation and be of great significance to enhance pollutant degradation in the hyporheic zone downstream of reservoirs.     

Drought analysis of monthly hydrological sequences: a case study of Canadian rivers

Abstract

Two parameters of importance in hydrological droughts viz. the longest duration, L_T and the largest severity, S_T (in standardized form) over a desired return period, T years, have been analysed for monthly flow sequences of Canadian rivers. An important point in the analysis is that monthly sequences are non-stationary (periodic-stochastic) as against annual flows, which fulfil the conditions of stochastic stationarity. The parameters mean, μ, standard deviation, σ (or coefficient of variation), lag1 serial correlation, ρ, and skewness, γ (which is helpful in identifying the probability distribution function) of annual flow sequences, when used in the analytical relationships, are able to predict expected values of the longest duration, E(L_T ) in years and the largest standardized severity, E(S_T ). For monthly flow sequences, there are 12 sets of these parameters and thus the issue is how to involve these parameters to derive the estimates of E(L_T ) and E(S_T ). Moreover, the truncation level (i.e. the monthly mean value) varies from month to month. The analysis in this paper demonstrates that the drought analysis on an annual basis can be extended to monthly droughts simply by standardizing the flows for each month. Thus, the variable truncation levels corresponding to the mean monthly flows were transformed into one unified truncation level equal to zero. The runs of deficits in the standardized sequences are treated as drought episodes and thus the theory of runs forms an essential tool for analysis. Estimates of the above parameters (denoted as μ_av, σ_av, ρ_av, and γ_av) for use in the analytical relationships were obtained by averaging 12 monthly values for each parameter. The product- and L-moment ratio analyses indicated that the monthly flows in the Canadian rivers fit the gamma probability distribution reasonably well, which resulted in the satisfactory prediction of E(L_T ). However, the prediction of E(S_T ) tended to be more satisfactory with the assumption of a Markovian normal model and the relationship E(S_T ) ≈ E(L_T ) was observed to perform better.     

Spectral shape proxies and nonlinear structural response

In this paper, spectral-shape-based intensity measures (IMs) are discussed with respect to ordinary, pulse-like and narrow-band records. First, the analyses address the ability of these IMs to capture the peak and cumulative damage potentials of ground motions. Second, a new vector-valued ground motion IM based on the spectral acceleration at the first mode of the structure, Sa(T₁), and a parameter proxy for the spectral shape, namely N_p, is introduced. The vector 〈Sa,N_p〉 is compared to other state-of-the-art IMs in terms of estimation of the seismic response of nonlinear single degree of freedom systems, reinforced concrete and steel moment resisting frames. Results show that 〈Sa,N_p〉 may be especially useful to represent the ground motion potential in the case of records with peculiar spectral shape. Further, it is shown that 〈Sa,N_p〉 has the properties of efficiency, sufficiency and scaling robustness. Finally, a scalar ground motion IM based on Sa(T₁) and N_p is also discussed, and the possibility to compute the seismic hazard analysis for it is illustrated.     

Prediction of capillary air-liquid interfacial area vs. saturation function from relationship between capillary pressure and water saturation

This short communication investigates if the capillary air-liquid interfacial area vs. saturation relationship A_lv(S) can be predicted from the capillary pressure vs. saturation relationship S(h), using the theoretical sample scale model of Diamantopoulos and Durner (2013, 2015). We selected three published experimental datasets, where S(h) and A_lv(S) relationships had been measured for the same porous media. The sample scale model was fitted to the retention curve S(h) of each porous medium and then used to predict the air-liquid interfacial area A_lv(S). We also included in the analysis the thermodynamic models of Leverett (1941) and Grant and Gerhard (2007). For two sandy materials and especially for high saturation values, the model predicted the capillary A_lv(S) successfully, which was in one case given by a pore-network model simulation (Kibbey and Chen, 2012) and in the other case experimentally determined (Brusseau et al., 2006). For glass bead experiments, the contact angle needed to be fitted to properly describe the experimental A_lv(S) curve.     

Artificial ground motion compatible with specified peak velocity and target spectrum

In this paper, a method, which synthesizes the artificial ground motion compatible with the specified peak velocity as well as the target acceleration response spectrum, was proposed. In this method, firstly, an initial acceleration time history α8^（0） （t）, which satisfies the prescribed peak ground acceleration, the target spectral acceleration ST（ω, ζ）,and the specified intensity envelope, is generated by the traditional method that generates the requency domain; secondly,α8^（0） （t）is further modulated by superimposing narrow-band time histories upon it in the time domain to make its peak velocity, approach the target peak ground velocity, and at the same time to improve its fitting precision to the target spectrum. Numerical examples show that this algorithm boasts high calculation precisions.     

2021年青海玛多7.4级地震强震动记录及特征分析

2021年5月22日,青海省果洛州玛多县发生7.4级地震,中国强震动观测网络在主震中捕获16组强震动数据。对48条三分向加速度记录进行基线校正、滤波等常规处理,计算相应的地震动参数,发现位于断层破裂前向位置的63DAW台NS向记录的地震动速度波形具有长周期分量丰富的特征。分析6个典型台站的单自由度加速度反应谱,并与我国建筑抗震设计谱比较,分析此次地震的频谱特性。将实际观测到的PGA、PGV和S_a(T=0.1 s、T=1.0 s、T=2.0 s、T=5.0 s)与国内广泛使用的几种地震动预测模型对比,研究此次地震的影响场。通过分析S_a-S_d曲线,探讨此次地震靠近断层区域地面运动大位移与桥梁落梁震害间可能存在的关系。     

Validation of 3-D basin structure models for long-period ground motion simulation in the Osaka basin, western Japan

We studied the applicability of two types of existing three-dimensional (3-D) basin velocity structure models of the Osaka basin, western Japan for long-period ground motion simulations. We synthesized long-period (3–20 s) ground motions in the Osaka basin during a M6.5 earthquake that occurred near the hypothetical Tonankai earthquake source area, approximately 200 km from Osaka. The simulations were performed using a 3-D finite-difference method with nonuniform staggered grids using the two basin velocity structure models. To study the ground motion characteristics inside the basin, we evaluated the wave field inside the basin using the transfer functions derived from the synthetics at the basin and a reference rock site outside the basin. The synthetic waveforms at the basin site were obtained by a convolution of the calculated transfer function and the observed waveform at the reference rock site. First, we estimated the appropriate Q values for the sediment layers. Assuming that the Q value depends on the S wave velocity V _S and period T, it was set to Q = (1/3V _S)(T ₀/T) where V _S is in m/s and the reference period T ₀ is 3.0 s. Second, we compared the synthetics and the observations using waveforms and pseudovelocity response spectra, together with a comparison of the velocity structures of the two basin models. We also introduced a goodness-of-fit factor to the pseudovelocity response spectra as an objective index. The synthetics of both the models reproduced the observations reasonably well at most of the stations in the central part the basin. At some stations, however, especially where the bedrock depth varies sharply, there were noticeable discrepancies in the simulation results of the models, and the synthetics did not accurately reproduce the observation. Our results indicate that the superiority of one model over the other cannot be determined and that an improvement in the basin velocity structure models based on simulation studies is required, especially along the basin edges. We also conclude that our transfer function method can be used to examine the applicability of the basin velocity structure models for long-period ground motion simulations.     

On the magnetic flux linkage of an electrically-conducting fluid

Abstract

A method has recently been proposed for finding the radius r_c of the electrically-conducting fluid core of a planet of outer radius r_s from observations of the magnetic field B in the accessible region near or well above the surface of the planet (Hide, 1978). The method is based on the supposition that when the magnetic field is produced by hydromagnetic dynamo action in the core, implying that the magnetic Reynolds number R there is large, (a) fluctuations in B can occur everywhere on the comparatively short advective time-scale τ _A associated with fluid motions in the core and so can fluctuations in the quantity N, defined for any closed surface S as the total number of intersection of magnetic lines of force with S, provided that S lies well outside the core, but (b) at the surface of the core, where lines of magnetic force emerge from their region of origin, concomitant fluctuations in N are negligibly small, of the order of τ _A/τ_O where τ _O (= Rτ_A ) is the Ohmic decay time of the core.

A proof of this supposition follows directly from the general expression derived in the present paper showing that when S is a material surface the time rate of change of N is equal to minus twice the line integral of the current density divided by the electrical conductivity around all the lines on S where the magnetic field is tangential to S. This expression (which Palmer in an accompanying paper rederives and extends to the relativistic case using the mathematical formalism of Cartan’s exterior calculus) also provides a direct demonstration of the well-known result that although high electrical conductivity, sufficient to make R ? 1, is a necessary condition for hydromagnetic dynamo action, such action is impossible in a perfect conductor, when R→ ∞.     

Selection and scaling of earthquake records in assessment of structures in low-to-moderate seismicity zones

The behaviour of a building during an earthquake depends on many factors including the characteristics of the ground motion, the configuration and dynamic characteristics of the structure, the hysteretic behaviour of the structural components, the building's non-structural parts, and the local soil behaviour. When conducting seismic assessment studies, controlling the first factor involves selecting and scaling appropriate earthquake records. It is a common practice to scale the records to fit a target acceleration response S_a determined from a pre-determined spectrum or attenuation law. In doing so, two key-problems must be faced: (i) the value adopted for the “effective” period T₁,which lengthens as the structure enters the non-linear range, and (ii) the reliability of the reference value adopted for S_a(T₁). In this work, the influence of selecting appropriate values for these two parameters—T₁ and S_a(T₁)—on the response of the building is investigated from the study of two structures, of 4 and 8 stories. The study does not focus on evaluating a particular seismic code; rather, it uses a European attenuation model as the target spectrum. Non-linear dynamic response analyses are carried out to investigate the Collapse Prevention State of these structures subjected to records representative of a Performance-based Earthquake Engineering level characterized by 2% exceedence in 50 years with 2475 years of return time. The records are selected from a database of European earthquakes, and the failure of the structures is evaluated by statistical means. From the results of the analyses, it is proposed that in seismic assessment studies the reference period for scaling should be 1.1 times the fundamental period of the structure, and at least 30 records should be used to guarantee reliable results.     

Attenuation characteristics of coda waves and estimation of Q c values in eastern China

For short-period near-earthquake records in eastern China, from the empirical attenuation formula of coda ground motion amplitudeA with timeτ: lgA=G−2. 235 lgτ, using the single scattering theory modified with epicentral distance, we obtain the curve family of corrected coda amplitudeA _c(r,t), andω/2Q _c values for each time interval of coda. From this,Q _c(f,h) values, which correspond to each observational average frequency and sampling depth, are calculated. The results substantially agree with those observationalQ _c values in Yunnan, Beijing and central Asia.     

Scaling of ductility and damage‐based strength reduction factors for horizontal motions

The conventional approach of obtaining the inelastic response spectra for the aseismic design of structures involves the reduction of elastic spectra via response modification factors. A response modification factor is usually taken as a product of (i) strength factor, R_S, (ii) ductility factor, R_μ, and (iii) redundancy factor, R_R. Ductility factor, also known as strength reduction factor (SRF), is considered to primarily depend on the initial time period of the single‐degree‐of‐freedom (SDOF) oscillator and the displacement ductility demand ratio for the ground motion. This study proposes a preliminary scaling model for estimating the SRFs of horizontal ground motions in terms of earthquake magnitude, strong motion duration and predominant period of the ground motion, geological site conditions, and ductility demand ratio, with a given level of confidence. The earlier models have not considered the simultaneous dependence of the SRFs on various governing parameters. Since the ductility demand ratio is not a complete measure of the cumulative damage in the structure during the earthquake‐induced vibrations, the existing definition of the SRF is sought to be modified with the introduction of damage‐based SRF (in place of ductility‐based SRF). A parallel scaling model has been proposed for estimating the damage‐based SRFs. This model considers damage and ductility supply ratio as parameters instead of ductility demand ratio. Through a parametric study on ductility‐based SRFs, it has been shown that the hitherto assumed insensitivity of earthquake magnitude and strong motion duration may not be always justified and that the initial time period of the oscillator plays an important role in the dependence of SRF on these parameters. Further, the damage‐based SRFs are found to show similar parametric dependence as observed in the case of the ductility‐based SRFs. Copyright © 2000 John Wiley & Sons, Ltd.     

An inversion of site response and Lg attenuation using Lg waveform

Based on spectral ratio method, a joint inversion method was used to obtain parameters of Lg wave attenuation and site response. The inversion method allows simple and direct (two-parameter) determination of Lg wave attenua- tion and site response from sparse spectral data, which are not affected by radiation pattern factor and different response of same instrument after geometrical spreading. The method was used successfully for estimating site re- sponse of stations of Zhejiang Seismic Network and measuring Lg wave attenuation. The study is based on 20 earth- quakes occurred in northeast of Taiwan with magnitude MS5.0~6.7 and 960 seismic wave records from 16 stations in Zhejiang area from 2002 to 2005. The parameters of site response and Lg attenuation were calculated with a fre- quency interval of 0.2 Hz in the range of 0.5 Hz to 10 Hz. Lg wave attenuation coefficient corresponding to U-D, E-W and N-S components are γ ( f )=0.00175 f 0.43485, γ ( f )=0.00145f 0.48467 and γ ( f )=0.0021f 0.41241, respectively. It is found that the site response is component-independent. It is also found that the site response of QIY station is significant above the frequency of 1.5 Hz, and that the site response of NIB station is low for most frequency     

Near‐fault ground motions,and the response of elastic and inelastic single‐degree‐of‐freedom (SDOF) systems

In order to investigate the response of structures to near‐fault seismic excitations, the ground motion input should be properly characterized and parameterized in terms of simple, yet accurate and reliable, mathematical models whose input parameters have a clear physical interpretation and scale, to the extent possible, with earthquake magnitude. Such a mathematical model for the representation of the coherent (long‐period) ground motion components has been proposed by the authors in a previous study and is being exploited in this article for the investigation of the elastic and inelastic response of the single‐degree‐of‐freedom (SDOF) system to near‐fault seismic excitations. A parametric analysis of the dynamic response of the SDOF system as a function of the input parameters of the mathematical model is performed to gain insight regarding the near‐fault ground motion characteristics that significantly affect the elastic and inelastic structural performance. A parameter of the mathematical representation of near‐fault motions, referred to as ‘pulse duration’ (T_P), emerges as a key parameter of the problem under investigation. Specifically, T_P is employed to normalize the elastic and inelastic response spectra of actual near‐fault strong ground motion records. Such normalization makes feasible the specification of design spectra and reduction factors appropriate for near‐fault ground motions. The ‘pulse duration’ (T_P) is related to an important parameter of the rupture process referred to as ‘rise time’ (τ) which is controlled by the dimension of the sub‐events that compose the mainshock. Copyright © 2004 John Wiley & Sons, Ltd.     

Empirical evaluation of peak ground velocity and displacement as a function of elastic spectral ordinates for design

In the framework of the revision of Part 1 of Eurocode 8, this study aims at developing new empirical correlations to compute peak values of ground velocity (PGV) and displacement (PGD) as a function of elastic spectral ordinates for design. At variance with the expressions for PGV and PGD currently adopted in the Eurocode 8, based solely on the peak ground acceleration (PGA), in this paper reference is made to spectral ordinates of the short and intermediate period range, namely S_s, which is the constant acceleration spectral ordinate, and S₁, which is the spectral ordinate at 1 s. On the one hand, a relationship between PGV and the product (S_s?S₁) was found based on the regression analysis on a high‐quality strong‐motion dataset. On the other hand, the PGD was estimated by extrapolating to long periods the constant displacement branch of the elastic response spectrum, introducing a correlation between the corner period T_D and S₁. For this purpose, results of a long period probabilistic seismic hazard assessment study for Italy, encompassing low to high seismicity areas, were considered. Furthermore, verification of the proposed relationship against strong‐motion records was carried out, and differences justified in terms of the concept of uniform hazard spectrum.