Simplified subgrade model for three-dimensional soil-foundation interaction analysis

A simple mechanical model is presented for the three-dimensional dynamic soil-structure interaction analysis of surface foundations. The model is made of one-dimensional vertical beams with distributed mass and horizontal springs which interconnect the two adjacent beams. Its parameters are rather uniquely related with the soil properties alone and thus are minimally dependent on the loading condition and the foundation conditions like geometry, flexibility and size. Formulations are provided to determine the model parameters from the soil properties. Solving the governing equations of this model, expressions for the subgrade behavior in response to the applied load and soil-foundation interaction are developed in analytical forms for various cases. The dynamic and static response of three-dimensional surface foundations are computed by these expressions. It is verified that the model is well capable of reproducing the three-dimensional soil-structure interaction behavior.     

Dynamic response analysis of submerged soil by thin layer element method

A thin layer element method is formulated to compute the dynamic response of submerged soil. The formulation is based on Biot's equation describing the dynamic behavior of fluid-saturated elasto-porous medium. The dynamic response of submerged soil is computed for various cases by using the developed formulation. The effects of submerged conditions are examined for submerged soil deposits with a water level at and above the ground surface. It is found that both submerged conditions and water body above the ground surface can considerably affect the dynamic response of soil deposits.     

Simplified mechanical subgrade model for dynamic response analysis of shallow foundations

Baranov proposed an approach for the treatment of the subgrade in the dynamic response analysis of a rigid foundation which is partially embedded. The approach is based on a Winkler approximation together with wave equation solutions developed for simplified conditions. The paper extends and modifies Baranov's approach so as to be applicable to surface foundation problems, by introducing an internal coupling mechanism. Closed form expressions are presented for the steady state harmonic responses of a massless rigid footing and an elastic beam on the surface of the presented subgrade model. These expressions can accommodate the inhomogeneity of the subgrade in a straightforward manner. The model is found to perform well in reproducing the continuous subgrade medium behaviour. The dynamic responses of structures are often highly dependent on the subgrade modelling and therefore great caution is needed in modelling the subgrade.     

Dynamic analysis of stacked rigid blocks

The dynamic behavior of a structural model of two stacked rigid blocks subjected to ground excitation is examined. Assuming no sliding, the rocking response of the system standing free on a rigid foundation is investigated. The derivation of the equations of motion accounts for the consecutive transition from one pattern of motion to another, each being governed by a set of highly nonlinear differential equations. The system behavior is described in terms of four possible patterns of response and impact between either the two blocks or the base block and the ground. The equations governing the rocking response of the system to horizontal and vertical ground accelerations are derived for each pattern, and an impact model is developed by conservation of angular momentum considerations. Numerical results are obtained by developing an ad hoc computational scheme that is capable of determining the response of the system under an arbitrary base excitation. This feature is demonstrated by using accelerograms from the Northridge, CA, 1994, earthquake. It is hoped that the two-blocks model used herein can facilitate the development of more sophisticated multi-block structural models.     

Fluid–structure interaction analysis of 3-D rectangular tanks by a variationally coupled BEM–FEM and comparison with test results

A variationally coupled BEM–FEM is developed which can be used to analyse dynamic response, including free-surface sloshing motion, of 3-D rectangular liquid storage tanks subjected to horizontal ground excitation. The tank structure is modelled by the finite element method and the fluid region by the indirect boundary element method. By minimizing a single Lagrange function defined for the entire system, the governing equation with symmetric coefficient matrices is obtained. To verify the newly developed method, the analysis results are compared with the shaking-table test data of a 3-D rectangular tank model and with the solutions by the direct BEM–FEM. Analytical studies are conducted on the dynamic behaviour of 3-D rectangular tanks using the method developed. In particular, the characteristics of the sloshing response, the effect of the rigidity of adjacent walls on the dynamic response of the tanks and the orthogonal effects are investigated. © 1998 John Wiley & Sons, Ltd.     

Dynamic analysis of flexible rectangular fluid containers subjected to horizontal ground motion

In this paper, an analytical method is proposed to determine the dynamic response of 3‐D rectangular liquid storage tanks with four flexible walls, subjected to horizontal seismic ground motion. Fluid–structure interaction effects on the dynamic responses of partially filled fluid containers, incorporating wall flexibility, are accounted for in evaluating impulsive pressure. The velocity potential in which boundary conditions are satisfied is solved by the method of separation of variables using the principle of superposition. The impulsive pressure distribution is then computed. Solutions based on 3‐D modeling of the rectangular containers are obtained by applying the Rayleigh–Ritz method using the vibration modes of flexible plates with suitable boundary conditions. Trigonometrical functions that satisfy boundary conditions of the storage tank such that the flexibility of the wall is thoroughly considered are used to define the admissible vibration modes. The analysis is then performed in the time domain. Moreover, an analytical procedure is developed for deriving a simple formula that evaluates convective pressure and surface displacements in a similar rigid tank. The variation of dynamic response characteristics with respect to different tank parameters is investigated. A mechanical model, which takes into account the deformability of the tank wall, is developed. The parameters of such a model can be obtained from developed charts, and the maximum seismic loading can be predicted by means of a response spectrum characterizing the design earthquake. Accordingly, a simplified but sufficiently accurate design procedure is developed to improve code formulas for the seismic design of liquid storage tanks. Copyright © 2013 John Wiley & Sons, Ltd.     

An analytical model of a deformable cantilever structure rocking on a rigid surface: experimental validation

This paper describes an experimental program to examine the dynamic response of deformable cantilevers rocking on a rigid surface. The primary goal of the tests is to verify and validate a dynamic rocking model that describes the behavior of these structures. The benchmark response data was obtained from shaking‐table tests on deformable rocking specimens with different natural vibration frequencies and different aspect ratios excited by analytical pulses and recorded ground motions. The responses computed using the model are found to be in good agreement with the benchmark test results. Widely used impact, restitution and damping assumptions are revisited based on the experiment results and the analytical model findings. Copyright © 2015 John Wiley & Sons, Ltd.     

A simplified model for lateral response of large diameter caisson foundations—Linear elastic formulation

The transient response of large embedded foundation elements of length-to-diameter aspect ratio D/B=2–6 is characterized by a complex stress distribution at the pier–soil interface that cannot be adequately represented by means of existing models for shallow foundations or flexible piles. On the other hand, while three-dimensional (3D) numerical solutions are feasible, they are infrequently employed in practice due to their associated cost and effort. Prompted by the scarcity of simplified models for design in current practice, we here develop an analytical model that accounts for the multitude of soil resistance mechanisms mobilized at their base and circumference, while retaining the advantages of simplified methodologies for the design of non-critical facilities. The characteristics of soil resistance mechanisms and corresponding complex spring functions are developed on the basis of finite element simulations, by equating the stiffness matrix terms and/or overall numerically computed response to the analytical expressions derived by means of the proposed Winkler model. Sensitivity analyses are performed for the optimization of the truncated numerical domain size, the optimal finite element size and the far-field dynamic boundary conditions to avoid spurious wave reflections. Numerical simulations of the transient system response to vertically propagating shear waves are next successfully compared to the analytically predicted response. Finally, the applicability of the method is assessed for soil profiles with depth-varying properties. The formulation of frequency-dependent complex spring functions including material damping is also described, while extension of the methodology to account for nonlinear soil behavior and soil–foundation interface separation is described in the conclusion and is being currently investigated.     

Bridge–embankment interaction under transverse ground excitation

Seismic performance and dynamic response of bridge–embankments during strong or moderate ground excitations are investigated through finite element (FE) modelling and detailed dynamic analysis. Previous research studies have established that bridge–embankments exhibit increasingly flexible performance under high‐shear deformation levels and that soil displacements at bridge abutment supports may be significant particularly in the transverse direction. The 2D equation of motion is solved for the embankment, in order to evaluate the dynamic characteristics and to describe explicitly the seismic performance and dynamic response under transverse excitations accounting for soil nonlinearities, soil–structure interaction and imposed boundary conditions (BCs). Using the proposed model, equivalent elastic analysis was performed so as to evaluate the dynamic response of approach embankments while accounting for soil–structure interaction. The analytical procedures were applied in the case of a well‐documented bridge with monolithic supports (Painter Street Overcrossing, PSO) which had been instrumented and embankment participation was identified from its response records after the 1971 San Fernando earthquake. The dynamic characteristics and dynamic response of the PSO embankments were evaluated for alternative BCs accounting for soil–structure interaction. Explicit expressions for the evaluation of the critical embankment length L_c are provided in order to quantify soil contribution to the overall bridge system under strong intensity ground excitations. The dynamic response of the entire bridge system (deck–abutments–embankments) was also evaluated through simplified models that considered soil–structure interaction. Results obtained from this analysis are correlated with those of detailed 3D FE models and field data with good agreement. Copyright © 2007 John Wiley & Sons, Ltd.     

Separation and sliding between soil and structure during strong ground motion

This paper presents an approach to the problem of separation and sliding between soil and structure in the finite element analysis of dynamic soil-structure interaction problems. Joint elements are arranged along the contact surface between soil and structure and they have a property such that tensile forces are not transmitted between the planes representing structure and soil in the finite element analysis. The dynamic properties governing the sliding are determined by the Mohr-Coulomb failure law determined from the cohesion and the friction angle between soil and structure. The proposed method is applied to (i) a model of a reactor building resting on the free surface of layered ground and (ii) a buried foundation structure. The numerical computations reveal the following results: that the translation is dominant in the motion of the structure when sliding is taking place between soil and structure, and that the rocking is dominant in the rest of the response. The amplitude of the response during sliding is increased on any one point of the structure and decreased on any one point of the ground compared with that of the fixed condition at the interface. In the case of the buried structure, it is observed in the computed results that the structure and soil move in the opposite direction along the vertical contact surface and are separated from each other in the near surface region during the strong phase of the excitation.     

土层平均剪切波速确定及其可靠性验证——基于钻井台阵强震动记录

震害资料显示,场地条件对地震动特性以及工程结构破坏程度影响显著。为减少因场地效应而造成的经济损失和社会影响,在进行场地地震反应分析时,需最大限度地减小因场地土层模型参数的不确定性引起的地震动评估偏差,为工程结构地震反应分析选取并生成适当的地震动输入。随着强震动观测技术的逐渐发展,大量可靠的钻井台阵记录为地震过程中场地观测点的动力反应提供了直接数据。以美国加州地区La Cienega钻井台阵强震动观测数据为基础,利用互相关函数,对不同强度地震作用下场地土层的平均剪切波速进行分析,并在此基础上,以Cyclic 1D为模拟平台,建立一维自由场地地震反应有限元分析模型。分析结果表明：通过钻井台阵地震动观测数据识别,得到场地平均剪切波速,能够反映该场地的动力特性,数值模拟计算结果和台阵地震动记录基本吻合,可为数值模型参数选取提供依据。     

土石坝地震响应的非线性剪切条模型与对比分析

对一维剪切条计算模型进行改进,提出了土石坝非线性地震反应的简化计算方法。首先将坝体沿坝高离散为一系列的具有不同剪切刚度与阻尼比等参数特性的层状体系,建立了各层的振动控制方程及其边值条件,进而采用数学物理方程方法进行了求解,确定了体系的振动特性,并根据振型叠加原理和Duhamel积分确定了坝体地震反应的线弹性解。采用等价线性化方法考虑坝料的动力非线性性质,通过对线弹性地震响应的反复迭代计算,使得各层土的模量和阻尼比与其相应的剪应变水平相协调,确定出与非线性坝体系统相等效的线性解答,并将所得到的地震响应作为非线性地震响应的近似解。最后,以均质坝和心墙坝作为算例进行了具体的数值计算,将所得结果与有限元数值解进行对比分析,论证了所提方法的适用性和合理性。     

高阶晃动振型对LNG储罐地震响应的影响

为了研究高阶晃动振型对LNG储罐地震响应的影响,考虑高阶晃动振型,建立LNG储罐的简化力学模型,推导LNG储罐的运动控制方程,给出了LNG储罐的基底剪力、倾覆弯矩和储罐内液体晃动波高的表达式。以某16×10⁴ m³ LNG储罐为例,采用大型通用有限元分析软件ADINA System对其进行有限元模型分析,验证其修正模型的有效性,结果表明:高阶晃动振型对基底剪力和倾覆弯矩几乎无影响,但对晃动波高影响显著,尤其是长周期地震动作用下,并且考虑高阶晃动振型的晃动波高存在延时效应。提出的简化力学模型修正公式与有限元分析结果吻合较好,可以准确地预测LNG储罐地震响应。     

Evaluation of site‐dependent constant‐damage design spectra for reinforced concrete structures

An investigation on the validity of the conventional design approach known as constant displacement ductility is carried out. The hysteretic behaviour described by the Modified Takeda model is taken to represent the characteristics of reinforced concrete structural systems. The results presented in the form of seismic damage spectra indicate that the conventional design approach may not be valid because cumulative damage is excessively high. The inelastic design spectra based on the constant‐damage concept are proposed in terms of simplified expressions. The expressions are derived from constant‐damage design spectra computed by non‐linear response analysis for SDOF systems subjected to ground motions recorded on rock sites, alluvium deposits, and soft‐soil sites. The proposed expressions, which are dependent on the local soil conditions, are functions of target seismic damage, displacement ductility ratio and period of vibration. The seismic damage of structures that have been designed based on this new design approach is also checked by a design‐and‐evaluation approach. The results are found to be satisfactory. Copyright © 2004 John Wiley & Sons, Ltd.     

Nonlinear 3D finite element analysis of soil–pile–structure interaction system subjected to horizontal earthquake excitation

The dynamic response of a seismic soil–pile–structure interaction (SSPSI) system is investigated in this paper by conducting nonlinear 3D finite element numerical simulations. Nonlinear behaviors such as non-reflecting boundary condition and soil–pile–structure interaction modeled by the penalty method have been taken into account. An equivalent linear model developed from the ground response analysis and the modified Drucker–Prager model are separately used for soil ground. A comparison of the two models shows that the equivalent linear soil model results in an underestimated acceleration response of the structure under this ground shaking and the soil behavior should be considered as a fully-nonlinear constitutive model in the design process of the SSPSI system. It was also observed that the dynamic response of the system is greatly affected by the nonlinearity of soil–pile interface and is not sensitive to the dilation angle of the soil. Furthermore, the effect of the presence of pile foundations on SSPSI response is also analyzed and discussed.     

On the seismic behavior of cylindrical base-isolated liquid storage tanks excited by long-period ground motions

A common effective method to reduce the seismic response of liquid storage tanks is to isolate them at base using base-isolation systems. It has been observed that in many earthquakes, the foregoing systems significantly affect on the whole system response reduction. However, in exceptional cases of excitation by long-period shaking, the base-isolation systems could have adverse effects. Such earthquakes could cause tank damage due to excessive liquid sloshing. Therefore, the numerical seismic response of liquid storage tanks isolated by bilinear hysteretic bearing elements is investigated under long-period ground motions in this research. For this purpose, finite shell elements for the tank structure and boundary elements for the liquid region are employed. Subsequently, fluid–structure equations of motion are coupled with governing equation of base-isolation system, to represent the whole system behavior. The governing equations of motion of the whole system are solved by an iterative and step-by-step algorithm to evaluate the response of the whole system to the horizontal component of three ground motions. The variations of seismic shear forces, liquid sloshing heights, and tank wall radial displacements are plotted under various system parameters such as the tank geometry aspect ratio (height to radius), and the flexibility of the isolation system, to critically examine the effects of various system parameters on the effectiveness of the base-isolation systems against long-period ground motions. From these analyses, it may be concluded that with the installation of this type of base-isolation system in liquid tanks, the dynamic response of tanks during seismic ground motions can be considerably reduced. Moreover, in the special case of long-period ground motions, the seismic response of base-isolated tanks may be controlled by the isolation system only at particular conditions of slender and broad tanks. For the case of medium tanks, remarkable attentions would be required to be devoted to the design of base-isolation systems expected to experience long-period ground motions.     

Optimal damper placement for building structures including surface ground amplification

A systematic method for optimal added damper placement in building structures is developed, taking into account the response amplification due to the surface ground. Non-linear amplification of the surface ground is described by an equivalent linear model. Hysteretic damping of the surface ground and radiational damping into the semi-infinite visco-elastic ground are included in the model. An original steepest direction search algorithm is applied to the interaction model. Closed-form expressions of the inverse of the coefficient matrix (tri-diagonal matrix) enable one to compute the transfer function and its derivative with respect to design variables very efficiently. It is shown that the ratio of the fundamental natural period of the structure to that of the surface ground is a key parameter for characterizing the optimal damper placement. Several examples for different soil conditions are presented to demonstrate the effectiveness and validity of the present method.     

Dynamic behaviour and seismic response of El Infiernillo dam

El Infiernillo, a 145 m high rockfill dam in Mexico built in a narrow V-shaped canyon, was subjected to eight major earthquakes since its construction. In this study, the dynamic dam response is analysed using (i) the recorded November 15, 1975 bedrock-crest acceleration histories and (ii) the results of a 1970 full-scale test conducted by UNAM, in which eight upstream–downstream and longitudinal resonant frequencies and configurations were observed and documented. These observed and seismically induced dynamic responses are compared herein to predictions of two different numerical models of El Infiernillo dam; a newly developed simplified three-dimensional (3D) model, and a 3D finite element model. The dynamic dam response characteristics are assessed, and performance of the employed numerical models is evaluated. It is found that (i) higher modes of vibration had participated significantly in the recorded seismic response and (ii) upstream–downstream response is well represented by the two numerical models employed. Using the simplified model, the September 19, 1985 earthquake non-linear response is computed and shown to compare satisfactorily with its recorded counterpart in the upstream-downstream and vertical directions. The largest computed dynamic accelerations, stresses and strains are found to occur within the upper third of the dam body.     

Boundary differential equation method: simplified dynamic soil stiffnesses for embedded rigid foundations

With a simplified model and Galerkin's weighted residual procedure, two simple differential equations of dynamic behavior of a bounded rectangular medium are established along the boundaries in the x- and y-direction in the medium. Solutions of these equations yield closed form expressions of soil stiffnesses for various cases of a partially embedded rigid foundation, including the stiffnesses per depth of foundation with rectangular base area and the stifnesses of strip foundation. The developed procedure provides the definition of the weight functions, which are used in Galerkin's method for weighted residual. In addition to these weight functions, their conjugators are also suitable for weight functions. When the soil depth is finite, the original weight functions fail to produce physically meaningful results in some frequency range but the conjugators do not fail at any frequencies. The developed equations to compute soil stiffnesses for embedded foundations are simple yet capable of calculating the responses close to those computed by the much more elaborated finite element method.     

四种地下结构抗震设计简化分析方法对比

在地下结构抗震设计简化分析方法中,强制反应位移法将土层变形施加在有限元模型侧边界模拟地震作用,反应加速度法将土层加速度施加到整个有限元模型上模拟地震作用,此外还有仅将土层加速度施加到土层模型上模拟地震作用的方法。上述方法均规避了反应位移法中关于弹簧刚度的取值问题,提高了计算效率。本文以1个双跨箱形结构为例,用动力时程分析的计算结果作为校核,分析了强制反应位移法、反应加速度法和仅将土层加速度施加到土体中的简化分析方法在不同侧边距条件下的计算精度,再结合常用的反应位移法,对比分析了4种简化分析方法的误差。分析结果表明:使用强制反应位移法时,侧边距取为1倍结构宽度导致的误差最小,反应加速度法和仅在土体施加加速度的简化方法对侧边距取值不敏感,反应位移法在角点造成的误差最大。