Effects of wave passage on the relevant dynamic properties of structures with flexible foundation

An evaluation of the wave passage effects on the relevant dynamic properties of structures with flexible foundation is presented. A simple soil–structure system similar to that used in practice to take into account the inertial interaction effects by the soil flexibility is studied. The kinematic interaction effects due to non‐vertically incident P, SV and Rayleigh waves are accounted for in this model. The effective period and damping of the system are obtained by establishing an equivalence between the interacting system excited by the foundation input motion and a replacement oscillator excited by the free‐field ground motion. In this way, the maximum structural response could be estimated from standard free‐field response spectra using the period and damping of the building modified by both the soil flexibility and the travelling wave effects. Also, an approximate solution for the travelling wave problem is examined over wide ranges of the main parameters involved. Numerical results are computed for a number of soil–structure systems to identify under which conditions the effects of wave passage are important. It comes out that these effects are generally negligible for the system period, but they may significantly change the system damping since the energy dissipation within the soil depends on both the wave radiation and the diffraction and scattering of the incident waves by the foundation. Copyright © 2001 John Wiley & Sons, Ltd.     

金银岛岩土台阵原位剪切波速剖面评估

土层的剪切波速是描述土动力学特性的重要参数之一。利用金银岛岩土台阵记录的8次浅源地方震的弱震动数据,使用解卷积的地震干涉测量法识别的剪切波走时,评估了该台阵两个水平方向的原位剪切波速剖面。结果表明:在44.2m深度以上,估计与实测的平均剪切波速剖面基本一致,而在44.2m到103.6m深度范围前者大于后者;本研究估计的平均剪切波速剖面比MEHTA等（2007）估计的平均剪切波速剖面更接近实测结果,在44.2m到103.6m深度范围,本研究估计的平均剪切波速与MEHTA等（2007）估计的平均剪切波速相近,二者均大于实测的平均剪切波速。     

Three-dimensional non-linear soil–building interaction analysis in the lakebed zone of Mexico city during the hypothetical Guerrero earthquake

The 1985 Michoacan earthquake (M=8·1) caused very severe damage to mid-rise buildings in the lakebed zone of Mexico City, which is approximately 400 km from the epicentre in the Pacific Ocean. In the present study, we perform a three-dimensional (3-D) non-linear soil–building interaction analysis for several types of low- to high-rise buildings during the hypothetical Guerrero earthquake, and try to understand the real cause of heavy damage to mid-rise buildings in the lakebed zone during the 1985 Michoacan earthquake. We make a reasonable estimation of the input earthquake motions and the local site effects. The non-linear soil-building interaction analysis explains the damage pattern observed during the 1985 earthquake, although other analyses do not. We realize that all the factors from the earthquake source to the building superstructure must be taken into account adequately. © 1998 John Wiley & Sons, Ltd.     

In-plane foundation-soil interaction for embedded circular foundations

Foundation soil interaction is studied using an analytical two-dimensional model, for circular foundations embedded in a homogeneous elastic half-space and for incident plane P- and SV- and for surface Rayleigh waves. The scattered waves are expanded in complete series of cyclindrical wave functions. A detailed analysis is presented of the foundation response to unit amplitude incident waves as a function of the type of incident waves and angle of incidence, the depth of the embedment and the foundation mass per unit length.It is shown that free-field translations and point rotation approximate well the foundation input motion only for very long incident waves. For shorter incident waves, those in general overestimate the foundation input motion. Neglecting the rotation of the foundation input motion (which is usually done in practice) may eliminate a major contribution to the base excitation of buildings and may cause nonconservative estimates of the forces in these buildings. Incident waves appear as ‘longer’ to a shallow foundation than to a deeper foundation. Therefore, deeper foundations are more effective in reflecting and scattering the short incident waves.     

Reconstruction of Tsunami Inland Propagation on December 26, 2004 in Banda Aceh, Indonesia, through Field Investigations

This paper presents the results from an extensive field data collection effort following the December 26, 2004 earthquake and tsunami in Banda Aceh, Sumatra. The data were collected under the auspices of TSUNARISQUE, a joint French-Indonesian program dedicated to tsunami research and hazard mitigation, which has been active since before the 2004 event. In total, data from three months of field investigations are presented, which detail important aspects of the tsunami inundation dynamics in Banda Aceh. These include measurements of runup, tsunami wave heights, flow depths, flow directions, event chronology and building damage patterns. The result is a series of detailed inundation maps of the northern and western coasts of Sumatra including Banda Aceh and Lhok Nga. Among the more important findings, we obtained consistent accounts that approximately ten separate waves affected the region after the earthquake; this indicates a high-frequency component of the tsunami wave energy in the extreme near-field. The largest tsunami wave heights were on the order of 35 m with a maximum runup height of 51 m. This value is the highest runup value measured in human history for a seismically generated tsunami. In addition, our field investigations show a significant discontinuity in the tsunami wave heights and flow depths along a line approximately 3 km inland, which the authors interpret to be the location of the collapse of the main tsunami bore caused by sudden energy dissipation. The propagating bore looked like a breaking wave from the landward side although it has distinct characteristics. Patterns of building damage are related to the location of the propagating bore with overall less damage to buildings beyond the line where the bore collapsed. This data set was built to be of use to the tsunami community for the purposes of calibrating and improving existing tsunami inundation models, especially in the analysis of extreme near-field events.     

Energy dissipation by nonlinear soil strains during soil–structure interaction excited by SH pulse

Three variants of a two-dimensional (2-D) model of a building supported by a rectangular, flexible foundation embedded in nonlinear soil are analyzed. The building, the foundation, and soil have different physical properties. The building is assumed to be linear, but the foundation and the soil can experience nonlinear deformations. It is shown that the work spent for the development of nonlinear strains in the soil can consume a significant part of the input wave energy, and thus less energy is available for excitation of the building. The results help explain why the damage, during the 1994 Northridge earthquake in California, to residential buildings in the areas that experienced large strains in the soil was absent or reduced.     

Implementation of viscoelastic mud-induced energy attenuation in the third-generation wave model,SWAN

The interaction of waves with fluid mud can dissipate the wave energy significantly over few wavelengths. In this study, the third-generation wave model, SWAN, was advanced to include attenuation of wave energy due to interaction with a viscoelastic fluid mud layer. The performances of implemented viscoelastic models were verified against an analytical solution and viscous formulations for simple one-dimensional propagation cases. Stationary and non-stationary test cases in the Surinam coast and the Atchafalaya Shelf showed that the inclusion of the mud-wave interaction term in the third-generation wave model enhances the model performance in real applications. A high value of mud viscosity (of the order of 0.1 m²/s) was required in both field cases to remedy model overestimation at high frequency ranges of the wave spectrum. The use of frequency-dependent mud viscosity value improved the performance of model, especially in the frequency range of 0.2–0.35 Hz in the wave spectrum. In addition, the mud-wave interaction might affect the high frequency part of the spectrum, and this part of the wave spectrum is also affected by energy transfer from wind to waves, even for the fetch lengths of the order of 10 km. It is shown that exclusion of the wind input term in such cases might result in different values for parameters of mud layer when inverse modeling procedure was employed. Unlike viscous models for wave-mud interaction, the inverse modeling results to a set of mud parameters with the same performance when the viscoelastic model is used. It provides an opportunity to select realistic mud parameters which are in more agreement with in situ measurements.     

Estimating kinematic interaction of raft foundations from earthquake records and its effects on structural response

A practical method for estimating kinematic interaction from earthquake records is presented. The kinematic interaction is characterized by a two-parameter model and these parameters can be estimated by using a frequency-domain systems identification method. The simple model can be used to model both wave passage effects and the effects of incoherent wave fields. Numerical simulation tests show that kinematic interaction parameters can be estimated to their best accuracy by using building base responses and the free-field excitation and can also be estimated by using building responses, base responses and the free-field excitation. The method was applied to two buildings with raft foundations and it was found that kinematic interaction was significant during earthquakes. Published theoretical models (wave passage effect) for vertically incident SH waves can be used to estimate the transfer functions up to 4–5 Hz and the models for horizontally propagating waves under-predict the estimated transfer functions by a significant amount at frequencies beyond about 1–2 Hz. Theoretical models for a massless rigid foundation under the excitation of an incoherent wave field predict the general trend of the estimated transfer function reasonably well over a large frequency range. The results of numerical examples show that the recorded response spectral attenuation of basement records at high frequencies with respect to the free-field is mainly caused by kinematic interaction, while the changes in storey shear and overturning moment in a structure due to soil flexibility are mainly the results of inertial interaction.     

Regional seismic responses of shallow basins incorporating site‐city interaction analyses on high‐rise building clusters

This study assesses the 3D amplification effects in shallow basins and quantifies the effects of site‐city interaction (SCI) on high‐rise buildings. A regional‐scale 3D spectral element simulation is conducted on the Tuen Mun‐Yuen Long basin, which contains multiple subbasins with heterogeneous and nonlinear soil profiles, while 3D city models with various building layouts are fully integrated into the basin model for our SCI study. We found a good correlation between spectral amplification factors and soil depths. Site response is significantly amplified at basin edges and centers due to surface waves generated at basin edges and the focusing effects stemming from 3D basin geometry. Transfer functions of 3D basins can be up to fourfold at fundamental frequencies as compared to 1D response, and further amplifications occur at high frequencies due to surface waves. In the SCI simulations, we observe wave trapping in the open space amid buildings resulting in energy concentration and up to twofold PGA amplifications. The wave trapping effect diminishes as the space between buildings increase beyond their range of influence (～100 m). The SCI analyses show that destructive kinetic energy in superstructures increases 28% in one horizontal direction but decreases 22% in the other. Our study concluded that, 1D site response analysis can significantly underestimate the seismic demand in shallow basins. Site‐city interaction of high‐rise buildings increases the short‐period spectra of ground motions, leading to an increase in their story accelerations by up to 50% and to a substantial decrease in the seismic safety of short structures in their vicinity.     

Hydrodynamic constraints and storm wave characteristics on a sub‐horizontal shore platform

Few studies of wave processes on shore platforms have addressed the hydrodynamic thresholds that control wave transformation and energy dissipation, especially under storm conditions. We present results of a field experiment conducted during a storm on a sub‐horizontal shore platform on the east coast of Auckland, New Zealand. Small (<0.5 m) locally generated waves typically occur at the field site, whereas during the experiment the offshore wave height reached 2.3 m. Our results illustrate the important control that platform morphology has on wave characteristics. At the seaward edge of the platform a scarp abruptly descends beneath low tide level. Wave height immediately seaward of the platform was controlled by the incident conditions, but near the cliff toe wave height on the platform was independent of incident conditions. Results show that a depth threshold at the seaward platform edge > 2.5 times the gravity wave height (0.05–0.33 Hz) is necessary for waves to propagate onto the platform without breaking. On the platform surface the wave height is a direct function of water depth, with limiting maximum wave height to water depth ratios of 0.55 and 0.78 at the centre of the platform and cliff toe, respectively. A relative ‘platform edge submergence’ (water depth/water height ratio) threshold of 1.1 is identified, below which infragravity (<0.05 Hz) wave energy dominates the platform energy spectra, and above which gravity waves are dominant. Infragravity wave height transformation across the platform is governed by the relative platform edge submergence. Finally, the paper describes the first observations of wave setup on a shore platform. During the peak of the storm, wave setup on the platform at low tide (0.21 m) is consistent with measurements from planar sandy beaches, but at higher tidal stages the ratio between incident wave height and maximum setup was lower than expected. Copyright © 2014 John Wiley & Sons, Ltd.     

Translation,torsion, and wave excitation of a building during soil–structure interaction excited by an earthquake SH pulse

A two-dimensional (2-D) model of a building supported by a rectangular, flexible foundation embedded in the soil is analyzed. The building, the foundation, and the soil have different physical properties. The building is assumed to be linear, but the soil and the foundation can experience nonlinear deformations. While the work spent for the development of nonlinear strains in the soil can consume a significant part of the input wave energy—and thus less energy is available for the excitation of the building—the nonlinear response in the soil and the foundation does not signficantly alter the nature of excitation of the base of the building. It is noted that the response of a building can be approximated by translation and torsion of the base for excitation by long, strong motion waves.     

SV波在饱和土体自由表面的反射

基于饱和两相介质弹性波动方程分析SV波在饱和土体自由表面的反射问题,引入波动方程的势函数解答,求解出二维问题中SV波入射情况下饱和土体自由场的位移、速度、加速度和应力响应。在饱和土体自由场响应解析解基础上,建立SV波入射下饱和土体自由场静、动力有限元模型。建模中考虑了如下几方面因素:(1)在不同分析步,对土体单元赋予不同材料本构。通过*model change命令进行单元生死设定,从而实现在初始应力场平衡的静力状态下采用DuncanChang本构模型,而地震波动输入时采用Davidenkov动力本构模型;(2)采用多孔介质黏弹性人工边界条件,在人工边界上分别施加固相和液相介质的弹簧和阻尼来模拟饱和土体中能量的传播;(3)将地震波转化为作用在人工边界上的等效地震荷载,施加到人工边界节点上;(4)土体单元采用4结点平面应变孔压单元(CPE4P)。有限元计算与解析解比较结果表明:SV波在垂直入射和掠入射时,竖向位移响应为零;在45°左右入射时,水平位移响应最大;60°左右入射时,竖向位移响应最大。这些结论与解析解吻合较好,本文模型为建立土-结构动力相互作用模型打下良好的基础。     

Two-dimensional translation,rocking, and waves in a building during soil-structure interaction excited by a plane earthquake SV-wave pulse

A two-dimensional (2-D) model of a building supported by a rectangular, flexible foundation embedded in the soil is analyzed for excitation by an incident plane SV-wave. The incidence is below the critical angle. The building is assumed to be anisotropic and linear while the soil and the foundation are assumed to be isotropic and can experience nonlinear deformations. In general the work spent for the development of nonlinear strains in the soil can consume a significant part of the input wave energy and thus less energy is available for the excitation of the building. We show that the energy distribution in the building depends on the nature of the incident wave and differs substantially between the cases of incident P- and SV-waves. However, for both excitation by a plane SV-wave pulse and excitation by a P-wave, we show that the nonlinear response in the soil and the foundation does not significantly change the nature of excitation of the base of the building. It is noted that the building response can be approximated by translation and rocking of the base only for excitation by long, strong motion waves.     

黄土丘陵河谷场地地震动特征研究

研究黄土丘陵河谷场地在地震作用下强地面运动特征的变化情况,可以揭示强震对该类场地上震害的触发机理。结合黄土高原的地貌特征,建立具有代表性的动力数值分析模型,通过输入不同幅值、频谱特性和持续时间的地震波,对起伏地形和覆盖黄土层共同影响下的黄土河谷场地进行地震反应分析。结果表明:黄土层和地形耦合作用控制了地表的PGA变化,使其趋于复杂,在同一输入波不同振幅作用下,与基岩河谷各测点相比,黄土覆盖河谷场地的地震动频谱幅值均有所增加,并且频谱主峰均向高频移动。在不同地震波输入下,场地不同部位的固有频率受地形高程和土层影响;而地震动大小和频谱幅值不仅与场地的基本频谱和地形起伏有关,也与输入地震波的频谱成分相关。输入波PGA与地震频谱特征都不变时,同一场地输出的地震频谱形状具有相似的特征,随着地震持时增长,能量向场地基本频率附近集中,从而可能导致场地上相应频率建筑物震动幅值增加,造成累积破坏。     

多层弹性半空间中的地震波（一）

为了了解地震震源和地球介质的性质,很有必要对地震波的辐射、传播和衰减问题作仔细的分析。作为一种近似,可以暂且忽略地球的曲率,把传播地震波的地球介质视为多层半空间。为简便起见,地震波的衰减问题另作考虑。这样,便需要研究多层、均匀、各向同性和完全弹性半空间中地震震源辐射的地震波传播问题。 用哈斯克尔（Haskell）矩阵法解多层介质中弹性波的传播问题是很方便的。如果     

基于反卷积干涉法的土层地震反应特性研究

美国阿拉斯加安克雷奇德兰尼公园的土层地震反应观测台阵建有6个井下观测点和1个地表观测点, 每个观测点设置一个三分量加速度传感器, 该台阵建成后记录了多次地震事件. 本文通过反卷积地震干涉法对这些地震记录进行土层反应分析, 根据该分析结果提取土层结构的等效剪切波速和阻尼比. 所有层位的地震记录对地表记录的反卷积波形均反映出在土层中传播的上、 下行波, 其上行与下行至每一层的时间差别明显. 根据每一层的到时差, 确定土层中的等效剪切波速和分层等效剪切波速, 该结果与现场土层等效剪切波速测试结果吻合较好; 根据等效剪切波速计算出的土层卓越频率, 其一致性亦较好; 根据上、 下行波的峰值分析确定该场地土层的等效阻尼比, 也与软土的阻尼比经验值相当. 这些参数为后续的土层反应模拟和土结相互作用研究奠定了良好的数据基础. 结果表明, 反卷积干涉法能够用于土层反应分析和土层地震反应特性的提取.      

土层剪切波速测试中的不确定性对场地地震动参数的影响分析——以Ⅲ类场地为例

本文选取山东地区2个Ⅲ类场地的工程地质勘探及土层剪切波速等资料,将土层厚度按5个深度段,每个分段给出了4个土层剪切波速的改变量,通过改变不同深度段土层剪切波速,建立了19种土层地震反应分析模型,分析了不同深度段,不同概率水平下土层剪切波速的变化对场地地震动参数的影响。研究表明,不同深度段土层剪切波速的变化对场地地震动参数的影响有差异。具体表现为,土层剪切波速的改变在1—10m、11—40m和地震输入界面处三个深度段对地震动加速度峰值影响较大;其中,41—70m和71—100m两个深度段剪切波速的改变对地震动加速度峰值影响小;在土层深度1—10m时,剪切波速降低,峰值变大,剪切波速的改变与峰值的改变呈负相关;在其它深度段,剪切波速降低,峰值变小,剪切波速的改变与峰值的改变呈正相关。剪切波速的改变在1—10m和11—40m两个深度段对地震加速度反应谱影响较大;在41—70m、71—100m和地震输入界面三个深度段对地震加速度反应谱影响很小。     

Threshold seismic energy and liquefaction distance limit during the 2008 Wenchuan earthquake

Estimating the possible region of liquefaction occurrence during a strong earthquake is highly valuable for economy loss estimation, reconnaissance efforts and site investigations after the event. This study identified and compiled a large amount of liquefaction case histories from the 2008 Wenchuan earthquake, China, to investigate the relationship between the attenuation of seismic wave energy and liquefaction distance limit during this earthquake. Firstly, we introduced the concept of energy absorption ratio, which is defined as the absorbed energy of soil divided by the imparted energy of seismic waves at a given site, and the relationship between the energy absorption ratio and the material damping ratio was established based on shear stress–strain loop of soil element and the seismic wave propagation process from the source to the site. Secondly, the threshold imparted seismic energy of liquefaction was obtained based on existing researches of absorbed energy required to trigger liquefaction of sandy soils and the ground motion attenuation characteristics of the 2008 Wenchuan earthquake, and the liquefaction distance limit of this earthquake was estimated according to the proposed magnitude–energy–distance relationship. Finally, the field liquefaction database of 209 sites of the 2008 Wenchuan earthquake was used to validate such an estimation, and the field observed threshold imparted seismic energy to cause liquefaction in recent major earthquakes worldwide was back-analyzed to check the predictability of the present method, and several possible mechanisms were discussed to explain the discrepancy between the field observations and the theoretical predictions. This study indicates that seismic energy attenuation and liquefaction distance limit are regional specific and earthquake dependent, and 382 J/m³ is the average level of threshold imparted seismic energy to cause liquefaction for loose saturated sandy soils, and the corresponding liquefaction distance limit is approximately 87.4 km in fault distance for a M_w?=?7.9 event in the Chengdu Plain. The proposed regional energy attenuation model and threshold imparted seismic energy may be considered as an approximate tool in evaluating the liquefaction hazard during potential earthquakes in this area.     

Soil characterization of Tınaztepe region (İzmir/Turkey) using surface wave methods and nakamura (HVSR) technique

To determine the shear wave velocity structure and predominant period features of T?naztepe in ?zmir, Turkey, where new building sites have been planned, active–passive surface wave methods and single-station microtremor measurements are used, as well as surface acquisition techniques, including the multichannel analysis of surface waves (MASW), refraction microtremor (ReMi), and the spatial autocorrelation method (SPAC), to pinpoint shallow and deep shear wave velocity. For engineering bedrock (V _s > 760 m/s) conditions at a depth of 30 m, an average seismic shear wave velocity in the upper 30 m of soil (AV_s30) is not only accepted as an important parameter for defining ground behavior during earthquakes, but a primary parameter in the geotechnical analysis for areas to be classified by V _s30 according to the National Earthquake Hazards Reduction Program (NEHRP). It is also determined that Z1.0, which represents a depth to V _s = 1000 m/s, is used for ground motion prediction and changed from 0 to 54 m. The sediment–engineering bedrock structure for T?naztepe that was obtained shows engineering bedrock no deeper than 30 m. When compared, the depth of engineering bedrock and dominant period map and geology are generally compatible.     

Blind deconvolution methodology for site-response evaluation exclusively from ground-surface seismic recordings

A novel blind deconvolution methodology for identification of the local site characteristics based on two seismograms recorded on the free surface of a sediment site is presented. The approach does not require recordings at depth nor at a nearby rock outcrop, and eliminates the need for any prior parameterization of source and site characteristics. It considers that the surface recordings are the result of the convolution of the ‘input motion at depth' with transfer functions (channels) representing the characteristics of the transmission path of the waves from the input location to each recording station. The input motion at depth is considered to be the common component in the seismograms (same input in a statistical sense). The channel characteristics are considered to be the part in the seismograms that is non-common, since the travel path of the waves from the input motion location at depth to each recording station is different, due to spatially variable site effects. By means of blind deconvolution, the algorithm eliminates what is common in the seismograms, namely the input motion at depth, and retains what is different, namely the transfer functions of the site from the input location to each recording station. It estimates the site response in both frequency and time domains, and identifies the duration of the site's transfer functions. The methodology is applied herein to synthetic data at realistic sites for performance validation. The blindly estimated results are in almost perfect agreement with the actual site characteristics, indicating that the approach is a promising new tool for seismic site-response identification from recorded data.