Ductility demands of tall buildings subjected to base rocking induced by Rayleigh waves

Performance-based earthquake engineering requires accurate estimation of structural response associated with different damage states because of strong ground motion. In recent work (Meza-Fajardo and Papageorgiou, 2018, EESD), we demonstrated that a significant contribution to the response of elastic soil-structure systems for high-rise buildings is attributed to base rocking associated with Rayleigh waves. The present paper presents results of a study investigating the effects of Rayleigh waves on the response of soil-structure systems with nonlinear behavior at the level of the superstructure. By introducing a rigid-elastic rotational spring at the base of the building, we take into account the stiffness reduction due to damage to the lateral load-resisting system at its root, and with it, increased displacement demands. Considering different levels of ductility and post-yield stiffness, we investigate the impact of rocking because of Rayleigh waves on maximum and residual interstory drift ratios. Our results indicate that rocking due to surface waves should be an important consideration for design and evaluation of tall buildings, as inelastic action elongates their effective natural period, and consequently, they are more prone to be damaged by resonance and excitation of extended duration because of Rayleigh waves.     

Response of tall buildings to base rocking induced by Rayleigh waves

By now, it is well known that long‐period surface waves can induce resonant response in high‐rise buildings, in particular those located in sedimentary basins. Rayleigh wave passage has been reported to induce rocking motion at the base of the buildings which can increase displacement demands significantly. However, the building behavior to base rocking has not been extensively studied because commercially available instruments do not record rotational components of ground motion, and thus, rocking time histories have not been available to the analysts. In a recent study, we proposed an effective method for estimating the rocking associated with Rayleigh waves, which takes into account their frequency‐dependent phase velocities. In the present work, we select a number of recorded seismic motions which include surface waves on sedimentary basins from recent well‐recorded earthquake events. Then, we proceed to identify and extract the recorded surface waves by using the technique mentioned above. Using realistic soil‐structure analytical models that have been proposed in the published literature for high‐rise buildings, we study their response to Rayleigh waves as they respond to both translational and rocking motions. Of particular interest is to compare the response of such structures with and without the presence of rotational motions due to surface waves. Using the roof displacement and the building interstory drift as response quantities, our results indicate that demands are controlled by rotational (rocking) motions associated with Rayleigh waves.     

Response of a circular foundation on a uniform half-space to elastic waves

An alternative technique to obtain the dynamic response of a massless rigid circular foundation resting on a uniform elastic half-space when subjected to harmonic plane waves is presented. The technique relies on the use of an integral representation involving the free-field ground motion and the contact tractions obtained in the course of calculating the dynamic force–displacement relationship of the foundation for external forces. Tables listing the translational and rotational components of the response of the foundation for non-vertically incident SH, P, SV and Rayleigh waves are presented.     

A brief theory and computing of seismic ground rotations for structural analyses

The seismic ground rotations are important with respect to spatial structural models, which are sensitive to the wave propagation. The rotational ground motion can lead to significant increasing of structural response, instability and unusual damages of buildings. Currently, the seismic analyses often take into account the rocking and torsion motions separately using artificial accelerograms. We present an exact analytical method, proposed by Nazarov [15] for computing of three rotational accelerograms simultaneously from given translational records. The method is based on spectral representation in the form of Fourier amplitude spectra of seismic waves, corresponding to the given three-component translational accelerogram. The composition, directions and properties of seismic waves are previously determined in the form of a generalized wave model of ground motion. It is supposed that seismic ground motion can be composed by superposition of P, SV, SH- and surface waves. As an example, the dynamic response analysis of 25-story building is presented. Here recorded (low-frequency) and artificial (high-frequency) accelerograms were used; each of them includes three translational and three rotational components. In this structural analysis, we have clarified primarily conditions under which rotational ground motion should be taken into account. Next, we have calculated three rotational components of seismic ground motion. Then they were taken as additional seismic loads components for further seismic analysis of the building. Note, soil–structure interaction (SSI) is not considered in this study. For computing, we use the special software for structural analyses and accelerogram processing (FEA Software STARK ES and Odyssey software, Eurosoft Co., Russia). It was developed and is used in engineering practice in the Central Research Institute of Building Constructions (TsNIISK, Moscow, Russia).     

Structural response for six correlated earthquake components

The paper examines the effect on the structural response of the inevitable correlation which exists between the six earthquake components acting along a set of structural axes. The rotational components are expressed in terms of the spatial derivatives of the translational components. For the calculation of response, modal analysis is employed so that ground response spectra can also be used as seismic input. A methodology is developed to obtain the maximum mean square response which can occur in a structure, irrespective of its orientation with respect to the impinging seismic waves. The application of this methodology for the calculation of design response is advocated, especially for asymmetric structures. For the assumed model of seismic wave motion, the numerical results show a significant contribution to the response from the rotational components. This contribution is, however, expected to be reduced by structural foundation averaging and interaction effects. Further studies with more complete models of seismic wave motions, and their interaction with structural foundations, are thus warranted for a realistic evaluation and characterization of the rotational inputs for design purposes.     

地震面波产生的地震动转动分量研究

本文利用弹性波动理论对地面转动分量，即瑞利（Rayleigh)波和乐夫（Love)波产生的转动分量进行了研究，给出了相应的计算公式和计算方法，特别注意到面波的散射效应对转动分量的影响，并将这一特性引入到转动分量的求取中，使问题的解决更切合于实际，最后选取实际地震记录，利用得到的公式计算出地震面波产生的转动分量。     

Structural response for stochastic kinematic interaction

The influence of stochastic kinematic interaction (SKI) on structural response is investigated in this paper. The SKI is evaluated through a computational model based on the boundary element method (BEM) formulated in the frequency domain. The singular integrals required in the computation of BEM are evaluated in a closed form. It is assumed that the foundation input motion (FIM) is the result of the superposition of many plane, stationary, correlated stochastic SH‐, P‐ and SV‐waves travelling within a homogeneous viscoelastic soil at different angles. The results obtained indicate that the effect of SKI on the foundation response is qualitatively similar to that of wave passage. Both effects involve a reduction of translational components of the response at intermediate and high frequencies and creation of a rotational response component at intermediate frequencies, which decreases at high frequencies. While, it is found that the SKI decreases the maximum response of structures built on embedded rigid strip foundations excited by SH‐ and P‐waves, it increases the maximum response for SV‐waves, except when the natural frequency of the structure is less than 0.5 Hz and for short structures excited by shallowly incident SV‐waves. Copyright © 2001 John Wiley & Sons, Ltd.     

我国在井-含水层系统对地震波同震响应方面的研究进展

总结了井-含水层系统对地震波同震响应的研究历史和现状,对我国最近几年发表的关于此类研究的文献进行了分类总结,得到了目前国内关于井-含水层系统对地震波同震响应研究的六个方向的研究进展,在文章的最后列举了目前井-含水层系统对地震波同震响应方面的研究成果以及今后的研究方向。     

Seismic rotational displacement of gravity walls by pseudo-dynamic method: Passive case

Prediction of the seismic rotational displacements of retaining wall under passive condition is an important aspect of design in earthquake prone region. In this paper, the pseudo-dynamic method is used to compute the rotational displacements of rigid retaining wall supporting cohesionless backfill under seismic loading for the passive earth pressure condition. The proposed method considers time, phase difference and effect of amplification in shear and primary waves propagating through both the backfill and the retaining wall. The influence of ground motion characteristics on rotational displacement of the wall is evaluated. Also the effects of variation of parameters like wall friction angle, soil friction angle, amplification factor, shear wave velocity, primary wave velocity, period of lateral shaking, horizontal and vertical seismic accelerations on the rotational displacements are studied. The rotational displacement of the wall increases substantially with increase in amplification of both shear and primary waves, time of input motion, period of lateral shaking and decreases with increase in soil friction angle, wall friction angle. The rotational displacements of the wall also increase when the effect of wall inertia is taken into account. Results are provided in graphical form.     

Response spectra for torsion,rocking and rigid foundations

Simple procedures are proposed for computing response spectra for torsional and rocking input ground motions assuming horizontally travelling waves of constant shape. It is shown that harmonic relationships exist between the rotational spectra and the corresponding translational spectra, and that SV rather than PSV is the correct basis for deriving the angular displacement, velocity and acceleration response spectra. An approximation enabling the use of the standard tripartite logarithmic response spectra is discussed. Simple expressions for ‘accidental’ eccentricity and rocking input effects are presented. Also proposed are multipliers to spectral ordinates to account for the filtering effects of rigid base mats resting on Winkler type foundations. For wave transit times shorter than half the natural period of the structure, these multipliers can be approximated by the frequency dependent averaging coefficients given in the literature, which are dependent, however, on the response, rather than the input, frequency.     

北京地区上空OH转动温度的季节性变化

在北京东北方向的兴隆天文台,自主搭建的大气辐射观测仪器对OH夜气辉从2011年12月开始进行观测.利用高分辨率的OH(8-3)带的振转光谱计算了转动温度,并与TIMED/SABER探测的温度进行了比较.观测表明,两年(2012—2013)的OH(8-3)带转动温度平均值为203.0±11.2K,有明显的季节变化,冬季高,夏季低,温差可达60K.与SABER观测温度的季节变化一致.对日平均的转动温度进行年振荡和半年振荡分量的拟合分析表明,年振荡强度(10.8K)远大于半年振荡(2.7K).研究还发现,不同夜晚转动温度变化形态差别很大,既有很强的潮汐控制的波动,又有相对短周期的波动.     

Simplification of strong ground motion considering inelastic responses of structures

Simplification of strong ground motions to 1 cycle sine waves was investigated from the elastic and inelastic earthquake response analyses and response analyses under sine wave input using single‐degree‐of‐freedom systems. Strong ground motions could be simplified to 1 cycle sine waves if large plastic deformations, with ductility factor more than 2, were assumed. This is because the approximate maximum responses from input sine waves are determined by the initial response cycle, due to period elongation and plastic energy dissipation of the systems. A sine wave whose acceleration amplitude is the peak ground acceleration (PGA) and whose period is that of an equivalent 1 cycle sine wave is proposed. The period of an equivalent sine wave is easily obtained from the elastic response acceleration spectrum of a seismic record. This means that the inelastic responses are approximately determined by the PGA and an equivalent 1 cycle sine wave period. Therefore, an equivalent 1 cycle sine wave period provides a single index to express the frequency characteristics of a strong ground motion. Copyright © 2000 John Wiley & Sons, Ltd.     

A model for nonlinear wave interactions

Phenomena associated with weakly nonlinear waves are encountered in optics, plasma physics, and geophysics. Specific applications in oceanography include ocean surface waves, internal gravity waves, Rossby waves, and interactions between these wave systems. Techniques for studying these phenomena include direct numerical integration, Boltzmann-like transport equations, and Langevin transient response functions. A simple model, called thetest wave model, is described and used to illustrate comparison of these methods.     

六分量弹性波场数值模拟与分析

地震波传播过程中,质点的振动不仅包括三个独立的平移部分,还包括三个独立的旋转部分.本文基于一阶速度-应力弹性波方程,采用分裂完全匹配层(SPML)的吸收边界条件,推导了时间导数二阶精度和空间导数高阶精度的交错网格有限差分格式的弹性波速度与应力各分量计算公式,模拟了各向同性介质中均匀模型和层状模型下的六分量波场,并对二维各向同性层状模型下的三个分量地震记录做高分辨率线性拉东变换得到各自的频散能谱.数值模拟分析结果表明:(1)旋转分量的能量要比平动分量弱的多;(2)在平动分量上,面波能量强,频率低,反射P波能量较强,反射S波能量稍弱;在旋转分量上,反射P波能量很弱,S波能量强;(3)与平动分量相比,旋转分量的频散能谱效果更好,能看到基阶和完整的高阶面波,即旋转分量能反映更多的地下介质信息.     

Symmetropy of earthquake patterns: asymmetry and rotation in a disordered seismic source

We demonstrate that the idea of symmetropy can be used for quantification of earthquake patterns. The symmetropy can be considered as a measure of asymmetry. A pattern is richer in asymmetry when the symmetropy is smaller. The specific results of its applications are obtained as follows. In a discrete model of a seismic source with self-organized criticality, the spatial patterns of earthquakes during critical states and sub-critical states are distinguished by the behaviour of the symmetropy: sub-critical patterns show that the symmetropy is approximately a constant but this has various values during critical states. The critical patterns show asymmetric property without any asymmetric force from the outside and without asymmetric intracellular rule. We show that the emergence of asymmetric patterns is a generic feature of dynamic ruptures in our model. Such a generic asymmetry results from the model which is an inherently discrete system consisting of finite-sized cells. These cells may represent geometrical disordered fault zones. We further discuss rotational motions that generate seismic rotational waves. In micromorphic continuum theory, such rotations are attributed to dynamic ruptures in disordered systems. We note that the concept of disorder in this theory is expressed by a set of finite-sized microstructures and is consistent with the concept of disorder modelled in the present study. Thus, we suggest that the spatially asymmetric patterns of earthquakes might be related to the rotational motions, because both come from dynamic ruptures in a discrete fault zone without a well-defined continuum limit.     

Hydrodynamic coefficients of a monolithic circular offshore structure

Expressions for both the rectilinear and rotational inertial and damping coefficients for a circular monolithic tower of uniform radius are derived. The analysis matches the fluid velocity, derived from potential theory, with the structural velocity in sway. That is, the motions of the tower are assumed to be in a vertical plane. The analysis is then applied to a tower composed of (lumped-mass) elements, where the expressions for the added-mass and damping coefficients are shown to be functions of wave number. The added-mass is shown to be a product of two wave systems: a travelling wave system, which is responsible for the radiation damping, and a standing wave system, called the evanescent system, which is attached to the structure. The added-mass of the evanescent system is negative for small wave numbers, while that of the travelling waves is positive. The negative sign simply means that the inertial force of the evanescent waves is 180° out of phase with that of the travelling system. Furthermore, it is shown that the contributions of the two wave systems to the total added-mass of the structure counteract each other, resulting in a total added-mass which varies gradually with the wave number. Finally, the analysis is applied to an experiment, and results of the analysis and the experiment are found to agree rather well.     

Accidental torsion in buildings due to base rotational excitation

This investigation is concerned with accidental torsion in buildings resulting from rotational excitation (about a vertical axis) of the building foundations as a result of spatially non-uniform ground motions. Because of this accidental torsion, the displacements and deformations in the structural elements of the building are likely to increase. This increase in response is evaluated using actual base rotational excitations derived from ground motions recorded at the base of 30 buildings during recent California earthquakes. Accidental torsion has the effect of increasing the building displacements, in the mean, by less than 5 per cent for systems that are torsionally stiff or have lateral vibration periods longer than half a second. On the other hand, short period (less than half a second) and torsionally flexible systems may experience significant increases in response due to accidental torsion. Since the dependence between this increase in response and the system parameters is complex, two simplified methods are developed for conveniently estimating this effect of accidental torsion. They are the ‘accidental eccentricity’ and the ‘response spectrum’ method. The computed accidental eccentricities are much smaller than the typical code values, 0.05bb or 0.1b, except for buildings with very long plan dimensions (b ≥ 50 m). Alternatively, by using the response spectrum method the increase in response can be estimated by computing the peak response to each base motion independently and combining the peak values using the SRSS rule.     

Flexible support of a pile embedded in unsaturated soil under Rayleigh waves

Under the action of Rayleigh waves, pile head is easy to rotate with a concrete pile cap, and pure fixed-head condition is rarely achieved, which is a common phenomenon for it usually occurs on the precast piles with insufficient anchorage. In addition, the propagation characteristics of Rayleigh wave have been changed significantly due to the existence of capillary pressure and the coupling between phases in unsaturated soil, which significantly affects the pile-soil interaction. In order to study the above problems, a coupled vibration model of unsaturated soil–pile system subjected to Rayleigh waves is established on the basis that the pile cap is equivalent to a rigid mass block. Meanwhile, the soil constitution is simplified to linear-elastic and small deformations are assumed to occur during the vibration phase of soil–pile system. Then, the horizontal dynamic response of a homogeneous free-field unsaturated soil caused by propagating Rayleigh waves is obtained by using operator decomposition theory and variable separation method. The dynamic equilibrium equation of a pile is established by using the dynamic Winkler model and the Timoshenko beam theory, and the analytical solutions of the horizontal displacement, rotation angle, bending moment and shear force of pile body are derived according to the boundary conditions of flexible constraint of pile top. Based on the present solutions, the rationality of the proposed model is verified by comparing with the previous research results. Through parametric study, the influence of rotational stiffness and yield bending moment of pile top on the horizontal dynamic characteristics of Rayleigh waves induced pile is investigated in detailed. The analysis results can be utilized for the seismic design of pile foundation under Rayleigh waves.     

The role of strong motion rotations in the response of structures near earthquake faults

Early studies of earthquake strong motion assumed linear materials and small deformations. It was observed that under favorable conditions (long waves), the accompanying rotational motions are usually small, and so their effects could be neglected. In 1932, when Biot opted for the vibrational method of solution of the dynamic response problems [Trifunac MD. 75th anniversary of the response spectrum method—a historical review. Soil Dyn Earthquake Eng 2008 [in press].] in his formulation of the response spectrum concept, his choice of the discrete mathematical models of buildings further led to the conditions that did not explicitly require consideration of the rotations [Trifunac MD. Buildings as sources of rotational waves, Chapter I.5. In: Teisseyre R, Nagahama H, Majewski E, editors. Physics of asymmetric continua: extreme and fracture processes. Heidelberg, Germany: Springer; 2008 [in press].]. The engineering profession was not prepared in the 1930s and 1940s for Biot's new theory and first had to learn the basic dynamics of structures before it could question the wisdom and consequences of the vibrational versus the wave-propagation approaches to the solution. Also, there were too many other concerns, often caused by the modeling simplifications, that pushed the studies of the rotational motion further down to the low levels of priority. Even today, 40 years after the arrival of digital computers and the emergence of powerful numerical computational capabilities, which uncovered unexpectedly large families of chaotic solutions accompanying large deformations, as well as nonlinear response [Trifunac MD. Nonlinear problems in earthquake engineering. In: Springer's encyclopedia of complexity and system science, 2008 [in press] [94].], most researchers continue to ignore the role of rotations. Had Biot chosen the wave-propagation approach for the solution of the earthquake engineering problems in 1932, the “progress” might have been faster. The wave representation can be differentiated with respect to a space coordinate, giving the rotations at a point directly. In contrast, the lumped-mass models in the vibrational approach do not make this possible, and the closest one can come to considering rotations is in terms of average, per-floor rotation, or drift.     

On plane strain problems in magneto-thermo-visco-elasticity

Summary Plane strain problems on magneto-thermo-visco-elastic interactions in a parallel union of the Kelvin and Maxwell bodies are investigated using the basic equations of electrodynamics and thermo-visco-elasticity. Assuming that the applied magnetic field is transverse to the plane of deformation and that the material is a perfect conductor of electricity, it is seen that the heat sources and the potential part of the body forces produce longitudinal waves only and the rotational part of the body forces gives rise to transverse waves only. The effect of deformation on magnetic permeability is equivalent to an anisotropic rescaling of the primary magnetic field. The effect of the applied magnetic field on waves produced by a plane heat source is equivalent to increasing the value of the material constants which results in an increase in the speed of the waves.