Computations of Seismic Passive Resistance and Uplift Capacity of Horizontal Strip Anchors in Sand

In this paper, the limit equilibrium method is used to compute seismic passive earth pressure coefficients and the vertical uplift capacity of horizontal strip anchors in presence of both horizontal and vertical pseudo-static earthquake forces. By considering a simple planar failure surface, distribution of soil reaction is obtained through the use of Kötter’s equation. Presence of pseudo-static seismic forces induces a considerable reduction in the seismic passive earth pressure coefficients. The reduction in seismic passive earth pressure coefficients increases with increase in magnitude of the earthquake accelerations in both horizontal and vertical directions and with increase in wall friction angle. The vertical uplift capacity of horizontal strip anchor is obtained for various values of soil friction angle, embedment ratio and seismic acceleration coefficients in both horizontal and vertical directions by using rigorous computational optimization. Proper justification for selected value of wall friction angle is established. Results are presented in the form of non-dimensional breakout factor for anchor. A significant reduction in breakout factor is observed in presence of both the seismic acceleration coefficients whereas breakout factor increases with increase in soil friction angle and embedment ratio even under the seismic condition. Angles of failure planes keep changing with change in seismic acceleration coefficients and failure zone shifts towards the critical direction of seismic acceleration coefficients. Present results are compared and found in good agreement with some specific available results in literature.     

Seismic Passive Earth Thrust on Retaining Walls with Cohesive Backfills Using Pseudo-Dynamic Approach

In this paper, the pseudo-dynamic approach is used to estimate seismic passive earth thrust on retaining walls with cohesive-frictional backfills. The time-dependent pseudo-dynamic approach considers the influence of dynamic parameters such as the velocity of primary and shear waves, the period of lateral shaking, and the phase and amplitude variations of horizontal and vertical earthquake accelerations with depth. The failure plane behind the wall is assumed to be planar. The analysis is based on the equilibrium of forces which act within the failure wedge. The obtained results show that the backfill cohesion increases both the seismic passive earth thrust and the failure plane inclination angle with the horizontal plane. It is also observed that both horizontal and vertical seismic accelerations have decreasing effect on seismic passive earth thrust as well as failure plane inclination angle. The results of present pseudo-dynamic analysis propose a lower solution for seismic passive earth thrust compared to earlier pseudo-static solution available in the literature.     

考虑饱和土的地铁车站结构非线性地震反应研究

饱和土广泛存在于自然界之中,将土体视为固－液二相介质,基于biot饱和土介质动力方程的u-p形式研究了反映饱和土无约束域能量辐射效应影响的等效黏弹性人工边界单元,基于地震波转化为作用于人工边界节点上等效荷载的方法实现了波动输入。采用算例验证了等效人工单元的精度。建立饱和土-地铁车站非线性有限元整体模型,分析了考虑双相介质饱和土的波动传播,对比了饱和土及单相土在地震作用下地铁车站的地震反应,并与振动台试验结果进行了比较,研究表明：地震作用下采用考虑多孔介质的特性分析地基土与单相固体介质相比更加合理。通过对车站结构关键构件的加速度、位移和内力等进行分析,研究了饱和土－地下结构相互作用体系的地震反应。     

近断层强地震动预测中的有限断层震源模型

提出了近断层强地震动预测中建立活断层上设定地震有限断层震源模型的方法和步骤.首先,根据地震地质和地震活动性调查以及地球物理勘探等资料,确定活断层的空间方位和滑动类型; 然后,根据地震定标律确定活断层的宏观震源参数; 第三,将高强体模型与k平方滑动模型相结合,产生断层破裂面上的混合滑动分布.在此基础上,预测了与1994年Northridge地震断层类型、矩震级(M_w6.7)基本一致的设定地震的有限断层震源模型.最后,将预测的有限断层震源模型与基于地震学的、使用动力学拐角频率的地震动随机合成方法相结合,预测了1994年Northridge地震近断层12个基岩台站的加速度时程,并和实际记录进行了对比.结果表明,用上述方法和步骤建立的有限断层震源模型是可行、实用的.      

Near-field ground motion during earthquake preparation process

Wave number discretization method is applied to study the near-field of seismic sources embedded in a cracked elastic solid. Near-field solutions are obtained for horizontal and vertical line forces. Effects of modifications in cracks of focal region on ground motion, in the near-field, are studied numerically for different

? values of crack density

? saturation of cracks

? width of cracks, and

? regimes of connection between cracks.

An earthquake process is assumed to be going through five major stages. These stages represent continuous accumulation of stress, interconnections between cracks leading to eventual failure and drainage of fluid from cracks after the major shock. Variations in the velocity ratio of waves noted from accelerograms verify the process of preparation of an earthquake.P wave contribution to vertical acceleration is negligible when the source is a vertical line force andS wave contributes only a little to horizontal acceleration when the source is a horizontal line force.     

A simplified approach for studying the influence of vertical accelerations on seismic response of slopes

ABSTRACT

A simplified approach is presented for estimating permanent displacements in slopes as a result of both vertical and horizontal seismic accelerations. A study of 52 earthquake records showed that the time difference between maximum horizontal and vertical accelerations varied between 0 and 10.3 s. The approach is illustrated for an earth dam embankment by analysing the effects of five of the above earthquake records. The approach combines a pseudo-static slope stability analysis for estimation of the critical (or yield) horizontal-vertical acceleration combinations, and a Newmark type displacement analysis. Guidelines are presented for conservative choice of soil strength parameters of saturated clays for use in the stability analysis. While permanent displacements of up to 40 cm were predicted without considering the vertical acceleration component, no additional displacement above 3.5 cm resulted when this component was included. The predicted additional displacement was consistently less than 10%, and in 50% of the analyses, vertical acceleration led to smaller predicted displacements. The simple approach may be applied in analysis for any slope using real earthquake records. Using existing, empirical expressions for permanent displacement, based only on horizontal accelerations, the effect of the vertical accelerations may be conservatively estimated by increasing the displacement by 10%.     

桩板式抗滑挡墙地震响应的振动台试验研究

汶川地震路基震害调查表明,在顺层或堆积体边坡中的桩板式抗滑挡墙具有良好的抗震性能。为了更好地了解该结构的抗震性能和优化抗震设计方法,以大型振动台模型试验为手段对其进行研究。为明确地震作用下桩板式抗滑挡墙的地震响应特性,试验采用缩尺的卧龙台站实测地震波对模型激励。试验结果揭示了土压力沿桩身分布规律、桩体位移和边坡岩土体加速度的地震响应特征。研究表明,地震土压力沿桩身呈非线性分布,竖向地震荷载对水平加速度有放大效果。所以,双向加载时的地震土压力比水平单向加载时大,但二者差距在地震基本烈度VII、VIII度区域不显著。滑坡推力、滑床对桩的土体抗力和桩身位移均与输入地震动峰值加速度成正比,即随着地震动峰值加速度的增加,加速度放大比增大;滑动面材料剪切强度折减,滑坡推力、土体抗力和抗身位移均增大,且增大速率加快。此外,结合试验成果,建议了桩板式抗滑挡墙设计时地震综合影响系数Cz的合理取值,对应地震基本烈度VII、VIII、IX度区分别为0.2、0.35、0.4。试验结果有助于揭示该结构抗震机制,也为其抗震设计提供了可靠依据。     

Deformation of a layered half-space due to a very long tensile fault

The problem of the coseismic deformation of an earth model consisting of an elastic layer of uniform thickness overlying an elastic half-space due to a very long tensile fault in the layer is solved analytically. Integral expressions for the surface displacements are obtained for a vertical tensile fault and a horizontal tensile fault. The integrals involved are evaluated approximately by using Sneddon’s method of replacing the integrand by a finite sum of exponential terms. Detailed numerical results showing the variation of the displacements with epicentral distance for various source locations in the layer are presented graphically. The displacement field in the layered half-space is compared with the corresponding field in a uniform half-space to demonstrate the effect of the internal boundary. Relaxed rigidity method is used for computing the postseismic deformation of an earth model consisting of an elastic layer of uniform thickness overlying a viscoelastic half-space.     

Laboratory tests and numerical simulations on the impact of subsurface barriers to saltwater intrusion

The dams are built to supply water needs of people, to produce electricity and to supply irrigation water in agricultural activities since ancient times. They are very important because of their contribution to energy production. The construction of dams is very difficult and costly. It is vital that they are resistant to all kinds of effects, since they are built to prevent flooding in the region. The most important and dangerous of these effects are earthquake forces. It is very important that the dams are designed to be durable. In this study, the Darideresi-II Dam reservoir located in Isparta was examined in terms of the behavior under various earthquake loads. To reflect the behavior of the dam in a more realistic way, the finite element method is used. While developing this model, the material specifications and boundary conditions are taken into account. The dam reservoir was modeled using ANSYS program, and its behavior under different earthquake accelerations was investigated. The deformation and stress forces under earthquake accelerations are taken into account. By examining these results, it is interpreted how the Darideresi-II Dam reservoir will behave during earthquake.     

A detailed investigation of block dynamic sliding by the discontinuous deformation analysis

An extensive examination of the discontinuous deformation analysis (DDA) in block dynamic sliding modeling is carried out in this paper. Theoretical solutions for a single block sliding on an arbitrarily inclined plane by applying the horizontal/vertical seismic loadings to the sliding block as acceleration time histories or to the base as constraint displacement time histories are derived. As compared with the theoretical solutions, for a single block sliding, the DDA predicts the sliding displacements and block interaction forces accurately under various base incline angles and friction angles under both the harmonic loadings and a real seismic loading. The vertical seismic component may influence the block sliding displacements to different extent, and the DDA can capture these phenomena successfully and give accurate results. For the calculation of the single block relative sliding, both the theoretical and the DDA solutions indicate that applying the seismic accelerations as constraint displacement time histories (derived by integrating the seismic accelerations twice) to the base is equivalent to applying the seismic accelerations as volume forces to the sliding block in the opposite directions. The DDA modeling also demonstrates that this conclusion still stands for the case of multi‐block sliding. Copyright © 2012 John Wiley & Sons, Ltd.     

Seismic response of pile foundations in liquefiable soil: parametric study

The performance of pile foundations in liquefiable soil subjected to earthquake loading is a very complex process. The strength and stiffness of the soil decrease due to the increase in pore pressure. The pile can be seriously destroyed by the soil liquefaction during strong earthquakes. This paper presents the response of vertical piles in liquefiable soil under seismic loads. A finite difference model, known as fast Lagrangian analysis of continua, is used to study the pile behavior considering a nonlinear constitutive model for soil liquefaction and pile?Csoil interaction. The maximum lateral displacement and maximum pile bending moment are obtained for different pile diameters, earthquake predominant frequencies, Arias intensities, and peak accelerations. It is found that the maximum lateral displacement and the maximum pile bending moment increase when the predominant earthquake frequency value decreases for a given peak acceleration value.     

基于包络的汶川大地震高频地震波辐射区域反演及近场加速度合成

合理地估计汶川破坏区域的地震动有助于地震灾害的研究.通过利用芦山地震记录建立的加速度包络衰减关系和汶川地震近场30个台站的加速度包络,基于线源模型,采用差分进化方法反演了汶川地震断层面上高频 (＞1 Hz) 辐射区域分布.结果表明:断层面上高频辐射分布很不均匀,辐射较强的区域主要位于:(1) 产生较大地表破裂的映秀、北川和南坝区域;(2) 映秀和北川等凹凸体的周边区域,包括震中东北侧60~90 km区域、北川和南坝东北侧30 km处;(3) 断层破裂停止的东北端约30 km长的区域.其中,破裂贯穿到地表的映秀、北川和南坝是低频和高频辐射都很强的区域.对于无观测记录场点,选择其临近且场地条件类似的台站加速度提取平稳随机过程,结合高频辐射分布和衰减关系得到的包络,合成了加速度时程,可为汶川地震结构震害分析提供地震动输入.      

Equivalent-Linear Seismic Ground Response Analysis of Some Typical Sites in Mumbai

Any earthquake event is associated with a rupture mechanism at the source, propagation of seismic waves through underlying rock and finally these waves travel through the soil layers to the particular site of interest. The bedrock motion is significantly modified at the ground surface due to the presence of local soil layers above the bedrock beneath the site of interest. The estimation of the amplifications in ground response due to the local soil sites is a complex problem to the designers and the problem is more important for mega cities like Mumbai in India, where huge population may get affected due to devastations of earthquake. In the present study, the effect of local soil sites in modifying ground response is studied by performing one dimensional equivalent-linear ground response analysis for some of the typical Mumbai soil sites. Field borelog data of some typical sites in Mumbai city viz. Mangalwadi site, Walkeswar site, BJ Marg near Pandhari Chawl site are considered in this study. The ground responses are observed for range of input motions and the results are presented in terms of surface acceleration time history, ratio of shear stress to vertical effective stress versus time, acceleration response spectrum, Fourier amplitude ratio versus frequency etc. The typical amplifications of ground accelerations considering four strong ground motions with wide variation of low to high MHA, frequency contents and durations are obtained. Results show that MHA, bracketed duration, frequency content have significant effects on the amplification of seismic accelerations for typical 2001 Bhuj motion. The peak ground acceleration amplification factors are found to be about 2.50 for Mangalwadi site, 2.60 for Walkeswar site and 3.45 for BJ Marg site using 2001 Bhuj input motion. The response spectrum along various soil layers are obtained which will be useful for designers for earthquake resistant design of geotechnical structures in Mumbai for similar sites in the absence of site specific data.     

Numerical identification of soil-structure interaction pressures

A numerical approach is presented for the ‘identification’ (or back calculation) of the earth pressure acting on embedded or retaining structures. The procedure is applicable to structures of any shape and requires a set of in situ measurements that may include displacements of points on the structure, values of concentrated forces, values of distributed loads at some locations, etc. Possible limiting values of the unknown loads, non-linear structural behaviour, varying accuracies of the input data are accounted for in the problem formulation. Depending upon the type of problem, the solution is reached by means of the unconstrained or constrained minimization of a suitably defined error function. As an example, the proposed approach is applied to the identification of the earth pressure acting on some typical geotechnical engineering structures.     

Transient response of reservoir–dam–soil systems to earthquakes

A direct time-domain numerical procedure is proposed to analyse the transient dynamic response of two-dimensional reservoir–dam–soil systems. The reservoir extends to infinity and the dam is supported by an unbounded soil. The structure with either linear or non-linear material properties is modelled by the Finite Element Method (FEM). The soil is assumed to be an elastic, isotropic and homogeneous half-space represented by a boundary condition in the form of generalized impedance determined by the transient Lamb's solution due to a uniformly distributed traction imposed on the free surface, Guan and Novak.¹ Moreover, a technique is developed to include the influence of the reservoir on the dam in terms of nodal accelerations along their interface at different time steps. The advantages of the proposed procedure are obvious. For example, it avoids any additional discretization of the boundaries except the soil–dam interface, and the influence matrix of the fluid is obtained explicitly using shape functions defined at the upstream face of the dam without the finite analysis of the reservoir so that it works very efficiently. Numerical results for a system consisting of reservoir, elastic dam and foundation subjected to the San Fernando, 1971 earthquake ground motion are presented.     

Ground motion estimation at Guwahati city for an Mw 8.1 earthquake in the Shillong plateau

In this paper, the ground motion at Guwahati city for an 8.1 magnitude earthquake on Oldham fault in the Shillong plateau has been estimated by stochastic finite-fault simulation method. The corresponding acceleration time histories on rock level at several sites in the epicentral region have been computed. These results are validated by comparing them with the estimates obtained from Medvedev–Sponheuer–Karnik (MSK) intensity observations of 1897 Shillong earthquake. Using the local soil parameters, the simulated rock level acceleration time history at Guwahati city is further amplified up to the ground surface by nonlinear site response analysis. The results obtained are presented in the form of peak ground acceleration (PGA) contour map. The maximum amplification for PGA over Guwahati city is as high as 2.5. Based on the simulated PGA, the liquefaction susceptibility at several locations in the city has been estimated. The results are presented in the form of contours of factor of safety against liquefaction at different depths below the ground surface. It is observed that over a large part of the Guwahati city, the factor of safety against liquefaction is less than one, indicating that the city is highly vulnerable to liquefaction in the event of this earthquake. The contour maps obtained can be used in identifying vulnerable areas and disaster mitigation.     

土石坝拟静力抗震稳定分析的强度折减有限元法

基于拟静力抗震设计概念,提出利用强度折减有限单元法分析土石坝的抗震稳定性,给出了两种确定地震惯性力的方法：（1）依据《水工建筑物抗震设计规范》[1],并结合有关土石坝动态分布系数计算了沿坝高分布的地震惯性力;（2）直接利用土石坝有限元地震动力反应分析得到的单元节点加速度反应,依据建议的方法确定坝体各单元节点的地震惯性力。将上述计算确定的地震惯性力与其他形式的外荷载共同作用到土石坝上,采用强度折减有限元法确定土石坝坝体的拟静力抗震安全系数。对于稳定渗流期,水位降落期等不同工况,或需要考虑振动孔隙水压力作用的饱和无黏性土填筑坝等不同计算条件,给出了使用折减强度有限元法分析坝体抗震稳定性的实现途径和方法。研究表明,有限元法对边界条件、复杂断面条件和材料分区及荷载组合均具有较强的适应能力,因此,使用有限元法分析土石坝抗震稳定性具有显著的优越性。     

盾构扩挖地铁车站结构地震反应特性振动台试验

开展了近远场地震动作用下盾构扩挖地铁车站结构的振动台试验,分析了砂土模型地基的水平位移、地表变形、加速度、土压力反应及模型结构的加速度、应变等。结果表明：模型地基-结构体系的地震响应对中低频成分发育的地震波反应更为强烈;强震作用下地铁车站结构具有明显的空间效应,地下结构的存在将会改变模型地基表面变形的分布模式。小震时模型结构中柱的加速度反应自下而上逐渐增加,而大震时其反应规律变成先增大后减小;车站结构中板的加速度反应最大、底板次之、顶板最小;小震时,同等深度处模型结构的加速度反应与模型地基土的加速度反应大小相当,侧墙的动土压力自下而上逐渐增大;大震时,模型结构的加速度反应明显大于同深度处模型地基土的加速度反应,动土压力的最大值发生在扩挖隧道的拱肩和中间部位。基于震后模型结构的宏观现象和拉应变幅值,给出了砂土地基中盾构扩挖车站结构的地震损伤演化机制。     

Acceleration of near-field scattering from an inhomogeneous spherical shell

The three dimensional scattering of near-field, from a point source, is studied for acceleration in the time domain. The perturbation method is applied to define the acceleration for the first order scattering from a weak inhomogeneity in a homogeneous surrounding. A body force, arising from the interaction between the primary waves and the inhomogeneity, acts as the source generating the scattered motion. The acceleration of scattered waves is related to the velocity and density fluctuations of the inhomogeneity. No restrictions are placed on the inhomogeneity size or locations of the source and receiver. Decoupling of scattered motion enables the identification of different phases. Integral expressions are derived for the scattering acceleration due to the incidence of near-field wave (from an impulsive point force) at a radially inhomogeneous volume element. These integrals are solved further for scattering from an inhomogeneous spherical shell. The accelerations for back scattering are obtained as a special case. These accelerations are simple analytically solvable expressions in closed form. Only spherical asymmetry ofP wave velocity inhomogeneity can affect the scatteredS acceleration. ScatteredP acceleration is affected by the gradient ofS wave velocity inhomogeneity. The back scattering of near-field from a spherical shell, is independent of radial inhomogeneity ofP wave velocity. Inhomogeneity with smoothly perturbedS wave velocity does not back-scatter any acceleration. Accelerations are computed numerically for scattering from a part of inhomogeneous spherical shell. Hypothetical models are considered to study the effects of the distances of spherical shell from source, receiver, its thickness and its position relative to the direction of impulsive force.     

Modeling of strong motion generation area of the Uttarkashi earthquake using modified semiempirical approach

The semiempirical approach based on envelope summation method given by Midorikawa (Tectonophysics 218:287–295, 1993) has been modified in this paper for modeling of strong motion generation areas (SMGAs). Horizontal components of strong ground motion have been simulated using modifications in the semiempirical approach given by Joshi et al. (Nat Hazard 71:587–609, 2014). Various modifications in the technique account for finite rupture source, layering of earth, componentwise division of energy and frequency-dependent radiation pattern. In this paper, SMGAs of the Uttarkashi earthquake have been modeled. Two different isolated wave packets in the recorded accelerogram have been identified from recorded ground motion, which accounts for two different SMGAs in the entire rupture plane. The approximate locations of SMGAs within the rupture plane were estimated using spatio-temporal variation of 77 aftershocks. Source parameters of each SMGA were calculated from theoretical and observed source displacement spectra computed from two different wave packets in the record. The final model of rupture plane responsible for the Uttarkashi earthquake consists of two SMGAs, and the same has been used to simulate horizontal components of acceleration records at different station using modified semiempirical technique. Comparison of the observed and simulated acceleration records in terms of root mean square error confirms the suitability of the final source model for the Uttarkashi earthquake.