A practical method to transform frequency dependent impedance to time domain

In order to perform time history earthquake response analyses with consideration to both the dynamic soil–structure interaction and the non‐linear behaviour of the structure, it is important to transform the soil impedance in the frequency domain to the impulse response in the time domain. In this paper, a new transform method with high practicality is proposed. First, the formulation of the proposed transform method is described. Next, the validity of the method is examined using an example problem whose impulse response is analytically obtained. Then, the impedance of the rigid foundation on 2‐layered soil is transformed to the time domain, and the characteristics of the impulse response are investigated. Finally, time history earthquake response analyses of a structure on the soil using the obtained impulse response are carried out. The validity and the efficiency of the proposed method are confirmed through these investigations. Copyright © 2005 John Wiley & Sons, Ltd.     

Soil-structure interaction experiment using impulse response

A new experimental method using a finite soil model with no special treatement on its boundaries is employed for soil-structure intration problems to simulate the semi-infinitenesss of the actual soil medium. The present method utlizies the characteristics of transient response to an impluse load to obtain the impedance functions and effective input motions for surface and embedded foundations. This technique is applicable to a linear elastic system whose impulse response decreases to a small enough value before observing the reflected waves. The experimentally obtained impedance functions and effective input motions are compared with those obtained by the direct boundary integral equation method and the hybrid approach. Good agreement between the xperimental and analytical results validates the present method as well as the accuracy of the numerical tools.     

基岩弹性刚度对土层地震反应的影响

将基岩上均匀、各向同性土层的地震反应,简化为置于弹性支座上的一维剪切梁模型进行分析。将地震激励假定为白噪声谱,在随机边界激励下,主要探讨了土层与基岩2种介质间的波阻抗比、波速比、土层厚度和阻尼特性对土层地震反应的影响。计算结果表明,对于一定的土层厚度,在一定阻尼比条件下,土层和基岩的阻抗比小到一定程度时,可以将基岩假定为刚性约束,而误差可以控制在一定的范围内。     

A numerical comparison of time and frequency‐domain marine electromagnetic methods for hydrocarbon exploration in shallow water

In shallow water the frequency domain controlled source electromagnetic method is subject to airwave saturation that strongly limits the sensitivity to resistive hydrocarbon targets at depth. It has been suggested that time‐domain CSEM may offer an improved sensitivity and resolution of these deep targets in the presence of the airwave. In order to examine and test these claims, this work presents a side‐by‐side investigation of both methods with a main focus on practical considerations, and how these effect the resolution of a hydrocarbon reservoir. Synthetic noisy data for both time‐domain and frequency domain methods are simulated using a realistic frequency dependent noise model and frequency dependent scaling for representative source waveforms. The synthetic data studied here include the frequency domain response from a compact broadband waveform, the time‐domain step‐response from a low‐frequency square wave and the time‐domain impulse response obtained from pseudo‐random binary sequences. These data are used in a systematic resolution study of each method as a function of water‐depth, relative noise and stacking length. The results indicate that the broadband frequency domain data have the best resolution for a given stacking time, whereas the time‐domain data require prohibitively longer stacking times to achieve similar resolution.     

Embedded foundation in layered soil under dynamic excitations

The critical step in the substructure approach for the soil–structure interaction (SSI) problem is to determine the impedance functions (dynamic-stiffness coefficients) of the foundations. In the present study, a computational tool is developed to determine the impedance functions of foundation in layered soil medium. Cone frustums are used to model the foundation soil system. Cone frustums are developed based on wave propagation principles and force-equilibrium approach. The model is validated for its ability to represent the embedded foundation in layered medium by comparing the results with the rigorous analysis results. Various degrees of freedom, such as, horizontal, vertical and rocking are considered for this study.     

A method for nonlinear dynamic response analysis of semi-infinite foundation soil

This paper presents development of a special finite difference method for the nonlinear dynamic response analysis of semi-infinite foundation soil. Semi-infinite domain is mapped into the finite domain using special mapping. For the region of engineering interest, mapping is isometrical, and for far field, shrink mapping which transforms an infinite interval into a finite interval is adopted. Using linear and nonlinear constitutive models, the responses of semi-infinite foundation soil are computed using a proposed method with a small mesh model and an extensive mesh model. Surface loadings or incident earthquake waves are applied to the models in the computations. Good agreements were obtained among the theoretical and computed results of the two models and the effectiveness of the proposed method was demonstrated.     

Stability of the time‐domain analysis method including a frequency‐dependent soil–foundation system

A number of methods have been proposed that utilize the time‐domain transformations of frequency‐dependent dynamic impedance functions to perform a time‐history analysis. Though these methods have been available in literature for a number of years, the methods exhibit stability issues depending on how the model parameters are calibrated. In this study, a novel method is proposed with which the stability of a numerical integration scheme combined with time‐domain representation of a frequency‐dependent dynamic impedance function can be evaluated. The method is verified with three independent recursive parameter models. The proposed method is expected to be a useful tool in evaluating the potential stability issue of a time‐domain analysis before running a full‐fledged nonlinear time‐domain analysis of a soil–structure system in which the dynamic impedance of a soil–foundation system is represented with a recursive parameter model. Copyright © 2015 John Wiley & Sons, Ltd.     

Improved methods to transform frequency‐dependent complex stiffness to time domain

A method to transform the frequency‐dependent complex stiffness to the impulse response in the time domain was proposed in the previous paper. However, there is a problem in that the accuracy and the convergence of the transformed impulse response are not good in some cases. Moreover, the hysteretic damping was not considered in the previous study although it is essential for practical purposes. In this paper, transform method improvements are proposed. First, the accuracy and the convergence are improved by taking the concept of virtual mass into account. Then, a more improved method for transforming the complex stiffness with large hysteretic damping to the time domain is proposed using the least square method. It is well known that the rigorous transform of the hysteretic damping is impossible because it is non‐causal. So this method is thought to be an approximate causalization process. Copyright © 2006 John Wiley & Sons, Ltd.     

Systematic assessment of irregular building–soil interaction using efficient modal analysis

In this paper an efficient methodology applying modal analysis is developed to assess systematically the combined soil–structure interaction and torsional coupling effects on asymmetric buildings. This method is implemented in the frequency domain to accurately incorporate the frequency‐dependent foundation impedance functions. For extensively extracting the soil–structure interaction effects, a diagonal transfer matrix in the modal space is derived. A comprehensive investigation of asymmetric building–soil interaction can then be conveniently conducted by examining various types of response quantities. Results of parametric study show that the increasing height‐to‐base ratio of a structure generally amplifies its translational and torsional responses. Moreover, both the translational and torsional responses are reduced for the case where the two resonant frequencies are well separated and this reduction is enhanced with the decreasing values of the relative soil stiffness and the height‐to‐base ratio. The most noteworthy phenomenon may be the fact that the SSI effects can enlarge the translational response if the structure is slender and the two resonant frequencies are very close. Copyright © 2001 John Wiley & Sons, Ltd.     

Dynamic stiffness and kinematic response of single piles in inhomogeneous soil

The effect of soil inhomogeneity on dynamic stiffness and kinematic response of single flexural elastic piles to vertically-propagating seismic SH waves is explored. A generalized parabolic function is employed to describe the variable shear wave propagation velocity in the inhomogeneous stratum. A layered soil with piece-wise homogeneous properties is introduced to approximate the continuous inhomogeneity in the realm of a Beam-on-Dynamic-Winkler-Foundation model. The problem is treated numerically by means of a layer transfer-matrix (Haskell–Thompson) formulation, and validated using available theoretical solutions and finite-element analyses. The role of salient model parameters such as pile-head fixity conditions, pile-to-soil stiffness ratio, surface-to-base shear wave velocity ratio and rate of inhomogeneity is elucidated. A new normalization scheme for inertial and kinematic response of such systems is presented based on an average Winkler wavenumber. With reference to long piles in moderately inhomogeneous soils, results indicate that: (a) kinematic pile response is essentially governed by a single dimensionless frequency parameter accounting for pile-to-soil stiffness ratio, pile slenderness and soil inhomogeneity and (b) definition of a characteristic pile wavelength allows an approximate estimation of pile elastodynamic response for preliminary design or analysis. Issues related to domain discretization and Winkler moduli are discussed.     

A numerical model for the dynamic analysis of inclined pile groups

The paper presents a numerical model for the dynamic analysis of pile groups with inclined piles in horizontally layered soil deposits. Piles are modelled with Euler–Bernoulli beams, while the soil is supposed to be constituted by independent infinite viscoelastic horizontal layers. The pile–soil–pile interaction as well as the hysteretic and geometric damping is taken into account by means of two‐dimensional elastodynamic Green's functions. Piles cap is considered by introducing a rigid constraint; the condensation of the problem permits a consistent derivation of both the dynamic impedance matrix of the soil–foundation system and the foundation input motion. These quantities are those used to perform inertial soil–structure interaction analyses in the framework of the substructure approach. Furthermore, the model allows evaluating the kinematic stress resultants in piles resulting from waves propagating in the soil deposit, taking into account the pile–soil–pile interactions. The model validation is carried out by performing accuracy analyses and comparing results in terms of dynamic impedance functions, kinematic response parameters and pile stress resultants, with those furnished by 3D refined finite element models. To this purpose, classical elastodynamic solutions are adopted to define the soil–pile interaction problem. The model results in low computational demands without significant loss of precision, compared with more rigorous approaches or refined finite element models. Copyright © 2015 John Wiley & Sons, Ltd.     

水平向不均匀土对箱基动力阻抗的影响

采用薄层元素法和有限单元法,建立了地基-箱形基础动力相互作用的三维分析模型,该模型可考虑箱形基础周围土的不均匀性影响.利用该模型,分析了水平方向不均匀土对箱形基础的水平和摇摆阻抗的影响.分析了水平方向不均匀土的厚度、剪切弹性模量以及材料阻尼比的影响.分析结果表明:水平方向不均匀土使得箱形基础的各种阻抗均降低,特别是阻抗的虚部.水平方向不均匀土对摇摆阻抗的影响较大.水平方向不均匀土的剪切弹性模量的影响较大而其材料阻尼的影响较小.随着水平方向不均匀土的厚度的增加和基础埋深的增加,其影响增大.     

On the performance of lumped parameter models with gyro‐mass elements for the impedance function of a pile‐group supporting a single‐degree‐of‐freedom system

Lumped parameter models with a so called “gyro‐mass” element (GLPMs) have been proposed recently in response to a strong demand for efficiently and accurately representing frequency‐dependent impedance functions of soil–foundation systems. Although GLPMs are considered to be powerful tools for practical applications in earthquake engineering, some problems remain. For instance, although GLPMs show fairly close agreement with the target impedance functions, the accuracy of the transfer functions and the time‐histories of dynamic responses in structural systems comprising GLPMs have never been verified. Furthermore, no assessment has been performed on how much difference appears in the accuracy of dynamic responses obtained from GLPMs and those from conventional Kelvin–Voigt models comprising a spring and a dashpot arranged in parallel with various frequency‐independent constants. Therefore, in this paper, these problems are examined using an example of 2×4 pile groups embedded in a layered soil medium, supporting a single‐degree‐of‐freedom system subjected to ground motions. The results suggest that GLPMs are a new option for highly accurate computations in evaluating the dynamic response of structural systems comprising typical pile groups, rather than conventional Kelvin–Voigt models. Copyright © 2011 John Wiley & Sons, Ltd.     

Combined paraxial-consistent boundary conditions finite element model for simulating wave propagation in elastic half-space media

In this paper, a finite element model of a soil island is coupled to both a consistent transmitting boundary and a paraxial boundary, which are then used to model the propagation of waves in semi-infinite elastic layered media. The formulation is carried out in the frequency domain while assuming plane strain conditions. It is known that a discrete model of this type, while providing excellent results for a wide range of physical parameters in the context of a half-space problem, may deteriorate rapidly at low frequencies of excitation. This is so because at low frequencies the various waves in the model eventually attain characteristic wavelengths which exceed the distance of the bottom boundary, which then causes that boundary to fail. Also, the paraxial boundaries themselves break down at very low frequencies. In this paper, this difficulty is overcome and the model׳s performance is improved upon dramatically by incorporating an artificial buffer layer sandwiched between the bottom of the soil medium and the underlying elastic half-space. Applications dealing with rigid foundations resting on homogenous or layered half-space media are shown to exhibit significant improvement. Following extensive simulations, clear guidelines are provided on the performance of the coupled model and an interpretation is given on the engineering significance of the findings. Finally, clear recommendations are provided for the practical use of the proposed modelling strategy.     

Stability and identification for rational approximation of frequency response function of unbounded soil

Exact representation of unbounded soil contains the single output–single input relationship between force and displacement in the physical or transformed space. This relationship is a global convolution integral in the time domain. Rational approximation to its frequency response function (frequency‐domain convolution kernel) in the frequency domain, which is then realized into the time domain as a lumped‐parameter model or recursive formula, is an effective method to obtain the temporally local representation of unbounded soil. Stability and identification for the rational approximation are studied in this paper. A necessary and sufficient stability condition is presented based on the stability theory of linear system. A parameter identification method is further developed by directly solving a nonlinear least‐squares fitting problem using the hybrid genetic‐simplex optimization algorithm, in which the proposed stability condition as constraint is enforced by the penalty function method. The stability is thus guaranteed a priori. The infrequent and undesirable resonance phenomenon in stable system is also discussed. The proposed stability condition and identification method are verified by several dynamic soil–structure‐interaction examples. Copyright © 2009 John Wiley & Sons, Ltd.     

Direct time-domain soil profile reconstruction for one-dimensional semi-infinite domains

We discuss the inverse medium problem associated with the reconstruction of the heterogeneous material profile of a semi-infinite (layered) soil medium, directly in the time domain, based on the complete waveform response of the medium to interrogating waves. To tackle the inversion process, we use a partial-differential-equation-constrained optimization approach, supplemented with a time-dependent regularization scheme. We introduce an absorbing boundary to truncate the semi-infinite extent of the physical domain, and propose two schemes to refine the reconstructed profiles: the first is based on iteratively re-positioning the truncation boundary until convergence, and the second is based on optimizing the observation period, so as to exclude records with information beyond the truncation boundary. We present numerical results that attest to the efficacy of the proposed schemes in reconstructing sharp profiles of semi-infinite soil domains using both noise-free and noisy data, while in the presence of absorbing boundaries.     

Receiver function method in reflection seismology

The receiver function method was originally developed to analyse earthquake data recorded by multicomponent (3C) sensors and consists in deconvolving the horizontal component by the vertical component. The deconvolution process removes travel path effects from the source to the base of the target as well as the earthquake source signature. In addition, it provides the possibility of separating the emergent P and PS waves based on adaptive subtraction between recorded components if plane waves of constant ray parameters are considered. The resulting receiver function signal is the local PS wave's impulse response generated at impedance contrasts below the 3C receiver.We propose to adapt this technique to the wide‐angle multi‐component reflection acquisition geometry. We focus on the simplest case of land data reflection acquisition. Our adapted version of the receiver function approach consists in a multi‐step procedure that first removes the P wavefield recorded on the horizontal component and next removes the source signature. The separation step is performed in the τ?p domain while the source designature can be achieved in either the τ?p or the t?x domain. Our technique does not require any a priori knowledge of the subsurface. The resulting receiver function is a pure PS‐wave reflectivity response, which can be used for amplitude versus slowness or offset analysis. Stack of the receiver function leads to a high‐quality S wave image.     

Dynamic plate-soil interaction - finite and infinite, flexible and rigid plates on homogeneous, layered or Winkler soil

An integral method to calculate the solution of a homogeneous or layered soil due to a harmonic point load is described. An infinite plate at the surface of the soil can be introduced in this integration in wavenumber domain, too. Finite structures on the soil are calculated by a combined finite element and boundary element method, which makes use of the point load solution of the soil. The compliance functions for a vertical point load and some vibration modes are calculated for realistic parameters of the plate and the soil and for a wide range of frequencies. The influence of the stiffness of the soil and the foundation is investigated, showing that the soil mainly affects the low-frequent response whereas the structural properties are more important at higher frequencies. A rigid approximation of flexible plates is only found at low frequencies, if the elastic length is used as the radius of a rigid disk. At higher frequencies, a characteristic behaviour of the flexible plate of approximately is observed, what is in clear contrast to the compliance of rigid foundations. A plate on a visco-elastic support (Winkler soil) shows similar displacements as a plate on a homogeneous half-space, but the maximal stresses between the plate and the soil are considerably smaller which is found to be more realistic for a plate on a layered soil. For practical applications, a normalized diagram and some explicit formulas of the exact and the approximate solutions of an infinite plate on a homogeneous half-space are given, which is a useful model to represent the soil-structure interaction of flexible foundations.     

Response of soil media due to impulse loads and isolation using trenches

This paper examines the effectiveness of open and filled trenches for screening Rayleigh waves due to impulse loads for 3D problems. The effects of the geometric parameters of trench reducing ground displacements are studied in detail. The computational scheme involves the use of finite element method to discretize a large computational domain and the resulting integration dynamic equilibrium equations are integrated to get the response using Newmark's β method. Soil has been idealized as linear, isotropic continuum.     

Vibration isolation by trenches in continuously nonhomogeneous soil by the BEM

Vibration isolation of structures from ground-transmitted waves by open trenches in isotropic, linearly elastic or viscoelastic soil with a shear modulus varying continuously with depth is numerically studied. Both an exponential and a linear shear modulus variation with depth are used in this work. Waves produced by the harmonic motion of a rigid surface machine foundation are considered. The problem is solved by the frequency domain boundary element method employing the Green's function of Kausel-Peek-Hull for a thin layered half-space. Thus only the trench perimeter and the soil-foundation interface need essentially to be discretized. The proposed methodology is first tested for accuracy by solving two Rayleigh wave propagation problems in nonhomogeneous soil with known analytical solutions and/or for which experimental results are available. Then the method is applied to vibration isolation problems and the effect of the inhomogeneity on the wave screening effectiveness of trenches is studied.