地下隧道轴向地震动土作用分析

以弹性基岩上覆层状场地中刚性衬砌隧道为模型,采用间接边界元方法求解衬砌隧道所受的沿轴向地震动土作用,通过参数分析揭示轴向动土作用的幅值大小、空间分布等基本规律。研究表明,土-隧道动力相互作用对地震动土作用的空间分布形式影响较小,但对地下隧道所受地震动土作用峰值大小具有显著影响,隧道主要位置点的地震动土作用峰值与隧道相应位置处自由场土层应力相比放大1.7~2.4倍。论文最后提出一个轴向地震动土作用的简化计算方法。     

Seismic analysis of underground tunnels by the 2.5D finite/infinite element approach

A procedure for the seismic analysis of underground tunnels using recorded free-field earthquakes based on the 2.5D finite/infinite element approach is presented. The near and far fields of the half space are modeled by finite and infinite elements, respectively. Using the 1D wave theory, the nodal force and displacement on the near-field boundary are computed for each spectral frequency of the earthquake. Then, equivalent seismic forces are computed for the near-field boundary for the earthquake spectrum. By assuming the soil-tunnel system to be uniform along the tunnel axis, the 2.5D approach can account for the wave transmission along the tunnel axis, which reduces to the 2D case for infinite transmission velocity. The horizontal and vertical components of the 1999 Chi-Chi Earthquake (TCU068) are adopted as the free-field motions in the numerical analysis. The maximal stresses and distribution patterns of the tunnel section under the P- and SV-waves are thoroughly studied by the 2.5D and 2D approaches, which should prove useful to the design of underground tunnels.     

Static equivalent method for the kinematic interaction analysis of single piles

This paper presents a static equivalent approach to estimate the maximum kinematic interaction effects on piles subjected to lateral seismic excitation. Closed-form expressions are reported for the evaluation of the maximum free-field soil movements and for the computation of maximum pile shear force and bending moments. Firstly, modal analysis, combined with a suitable damped response spectrum, is used to evaluate the maximum free-field response. Secondly, the pile is schematised as a Winkler's beam subjected to equivalent static forces defined according to soil vibration modal shapes and amplitude. The method may be applied by using response spectra suggested by National Standards or those obtained with accelerograms. The procedure proposed may be conveniently implemented in simple spreadsheets or in commercial finite element programs and easily used by practicing engineers. Method accuracy is demonstrated by comparing the results with those obtained with a more rigorous model. Good results may be achieved by considering only the first soil vibration mode making the procedure straightforward for practical design purposes.     

Evaluation of the seismic response of stone pagodas using centrifuge model tests

This study attempts to propose dynamic centrifuge model tests as a method of seismic risk assessment in order to discover how stone architectural heritages with masonry structures have endured seismic load, and whether there is any possibility of future earthquake damage. Dynamic centrifuge tests have been conducted for one fifteenth scale models of Seok-ga-tap and the five-storey stone pagoda of Jeongnimsa temple site, which are Korean representative stone pagodas. In order to make input motions of the earthquake simulator, site investigation and site-specific response analysis have been performed. The models of two stone pagodas, which have the same number of pieces with the real structures, have been produced and the dynamic centrifuge tests have been conducted for the model pagodas. Accelerometers were attached at different heights of the pagoda. The measured acceleration records and frequency responses were analysed during dynamic centrifuge test. Two real earthquake records, Hachinohe and Ofunato earthquakes and a sweeping signal with ranged frequency were utilised for input motions of dynamic centrifuge tests to evaluate the behaviour of the stone pagodas. For Seok-ga-tap models, it was observed that acceleration tends to be amplified with height. The third floor body shows at most 2.5 amplification of acceleration in comparison to the surface ground. The amplification was at a frequency of 3.83 Hz and it was considered as the natural frequency of the pagoda. For the five-storey stone pagoda, the seismic wave energy significantly reduced while it passed the first body floor, and then the peak acceleration was gradually amplified upwards. It was found that the pagodas did not collapse when the peak acceleration of ground surface was raised to 0.4 g. Given that the maximum design seismic acceleration specified in Korean seismic design guide is 0.22 g and the amplification ratio of peak acceleration in the supporting ground of the pagodas ranges from 1.45 to 1.74, it can be shown that the two pagodas are stable against 2400-year return period earthquake level, and have excellent seismic performance.     

Bounding surface plasticity model for the seismic liquefaction analysis of geostructures

This paper presents the constitutive relations and the simulative potential of a new plasticity model developed mainly for the seismic liquefaction analysis of geostructures. The model incorporates the framework of critical state soil mechanics, while it relies on bounding surface plasticity with a vanished elastic region to simulate the non-linear soil response. Key constitutive ingredients of the new model are: (a) the inter-dependence of the critical state, the bounding and the dilatancy (open cone) surfaces on the basis of the state parameter ψ, (b) a (Ramberg–Osgood type) non-linear hysteretic formulation for the “elastic” strain rate, (c) a discontinuously relocatable stress projection center related to the “last” load reversal point, which is used for mapping the current stress point on model surfaces and as a reference point for introducing non-linearity in the “elastic” strain rate and finally (d) an empirical index of the directional effect of sand fabric evolution during shearing, which scales the plastic modulus. In addition, the paper outlines the calibration procedure for the model constants, and exhibits its accuracy on the basis of a large number of laboratory element tests on Nevada sand. More importantly, the paper explores the potential of the new model by presenting simulations of the VELACS centrifuge tests of Models No 1 and 12, which refer to the free-field liquefaction response of Nevada sand and the seismic response of a rigid foundation on the same sand, respectively. These simulations show that the new model can be used successfully for the analysis of widely different boundary value problems involving earthquake soil liquefaction, with the same set of model constants calibrated on the basis of laboratory element tests.     

Seismic behaviors of earth-core and concrete-faced rock-fill dams by dynamic centrifuge tests

The investigation on the seismic behavior of dams becomes crucial but is limited to lack of experimental or field data. This paper aims to experimentally simulate two major dam types of earth-core rock-fill dam and concrete-faced rock-fill dam by dynamic centrifuge tests to investigate the seismic response of the dam. A series of staged centrifuge tests was performed by applying real earthquake records from 0.05 to 0.5g. The distributions of amplification ratio differed depending on the magnitude of earthquake loading and the zoning condition. The amplification ratio at the crest increased in the bedrock acceleration that exceeds 0.3g and strongly influenced by the loosening behavior of the upper part. The residual settlements and horizontal displacement at the dam crest were small. Shallow surface sliding was dominant failure. The maximum tensile stress on the face slab by dynamic loading occurred at a height of around 4/5 near the upstream water level. Finally, two-dimensional numerical simulations were performed in an effort to simulate the centrifuge models. The centrifuge tests and numerical analysis obtained mostly comparable results, thus confirming that centrifuge modeling reasonably simulates the seismic behavior of dams.     

Centrifuge modeling of seismic response of layered soft clay

Centrifuge modeling is a valuable tool used to study the response of geotechnical structures to infrequent or extreme events such as earthquakes. A series of centrifuge model tests was conducted at 80g using an electro-hydraulic earthquake simulator mounted on the C-CORE geotechnical centrifuge to study the dynamic response of soft soils and seismic soil–structure interaction (SSI). The acceleration records at different locations within the soil bed and at its surface along with the settlement records at the surface were used to analyze the soft soil seismic response. In addition, the records of acceleration at the surface of a foundation model partially embedded in the soil were used to investigate the seismic SSI. Centrifuge data was used to evaluate the variation of shear modulus and damping ratio with shear strain amplitude and confining pressure, and to assess their effects on site response. Site response analysis using the measured shear wave velocity, estimated modulus reduction and damping ratio as input parameters produced good agreement with the measured site response. A spectral analysis of the results showed that the stiffness of the soil deposits had a significant effect on the characteristics of the input motions and the overall behavior of the structure. The peak surface acceleration measured in the centrifuge was significantly amplified, especially for low amplitude base acceleration. The amplification of the earthquake shaking as well as the frequency of the response spectra decreased with increasing earthquake intensity. The results clearly demonstrate that the layering system has to be considered, and not just the average shear wave velocity, when evaluating the local site effects.     

场地地震反应分析中侧向人工边界距离选取及影响因素分析

结合均匀场地模型和分层均匀场地模型,采用数值试验的方法研究了场地地震反应分析中侧向人工边界距离的选取问题。分析了各主要因素对人工边界距离选取的影响,得到了如下结论:单侧人工边界距离与土层深度之比L/h是影响人工边界距离选取的主要参数;随着侧向人工边界距离的增加,有限元模型的计算误差并不一直减小,在局部范围内,会出现计算误差随着边界距离的增加而增大的现象;场地硬度、阻尼比和激振频率等因素不会影响计算误差随边界距离的变化规律,但采用自由场边界时受输入地震动频谱特性的影响较为明显;对水平激振下的一维场地模型来说,采用结合自由度边界时所需边界距离最短,自由场边界次之,采用粘性边界达到收敛时所需边界距离最长;无论哪种边界条件,当边界距离大于某一范围后,都能得到稳定的计算结果;相对于均质场地模型来说,实际的分层均匀场地建模时可适当减小侧向人工边界距离。     

Seismic response of retaining walls with cohesive backfill: Centrifuge model studies

Observations from recent earthquakes show that retaining structures with non-liquefiable backfills perform extremely well; in fact, damage or failures related to seismic earth pressures are rare. The seismic response of a 6-m-high braced basement and a 6-m free-standing cantilever wall retaining a compacted low plasticity clay was studied in a series of centrifuge tests. The models were built at a 1/36 scale and instrumented with accelerometers, strain gages and pressure sensors to monitor their response. The experimental data show that the seismic earth pressure on walls increases linearly with the free-field PGA and that the earth pressures increase approximately linearly with depth, where the resultant acts near 0.33 H above the footing as opposed to 0.5–0.6 H, which is suggested by most current design methods. The current data suggest that traditional limit equilibrium methods yield overly conservative earth pressures in areas with ground accelerations up to 0.4g.     

Dam-reservoir-massed foundation system and travelling wave along reservoir bottom

In the present paper, effects of non-uniform excitation due to spatially variation of seismic waves under the reservoir bottom on linear and nonlinear responses of arch dams are studied. Foundation is assumed to be massed and infinite elements are utilized to model semi-infinite medium via the far-end boundary of the foundation FE model. A continuum crack propagation model based on the smeared crack approach in 3D space is introduced. Reservoir's water is assumed compressible and the coupled system is solved using the staggered method. As a case study, Amir-Kabir double curvature arch dam in Iran is selected to investigate seismic behavior of the system. Two cases are analyzed in which wave travelling velocities are taken as 650 m/s and infinity. According to the results, non-uniform excitation leads to reduction in seismic response of the system and, in addition, frequency content of the response is different completely in comparison with the system under the uniform excitation.     

Significance of small strain damping and dilation parameters in numerical modeling of free-field lateral spreading centrifuge tests

Lateral spreads are complex dynamic phenomena that are challenging to represent numerically. In this paper numerical models are developed and calibrated using the displacement, acceleration, and pore water pressure time histories recorded in a free-field lateral spreading centrifuge test. The calibrated numerical model then is used to predict another free-field lateral spreading centrifuge test using the same soil profile but different input acceleration time history. The computed response shows good agreement with the centrifuge test measurements. This paper demonstrates that even in a large strain problem, such as lateral spreading, small strain damping plays an important role in numerical simulation results; it also shows the need to have pressure dependent dilation parameters in the employed soil constitutive model implemented in order to simultaneously reproduce measurements of pore water pressure, acceleration and lateral displacement.     

土-结构相互作用分析中几类地震输入方式的比较

地震作用下土-结构动力相互作用的整体有限元分析需要在人工边界处输入地震动。目前可能采用的地震输入方法包括黏弹性边界自由场输入方法、自由场应力方法、自由场位移方法以及侧边界自由方法。由于采用近似人工边界条件或者未完全考虑地震自由场,上述地震输入方法均为近似方法。本文以大开地铁车站二维有限元分析为例,根据规范建议的边界位置,研究了上述地震输入方法的精度,研究成果可为土-结构相互作用分析的合理地震输入提供一定参考。     

平面P-SV波入射时非均匀饱和土 自由场地的响应

基于Biot多孔介质波动模型,研究了非均匀饱和土层对平面P-SV波入射时的动力响应.考虑饱和土地基的物理力学特性沿厚度方向连续变化,利用亥姆霍兹矢量分解原理和动力刚度法,分析了平面入射P-SV波在非均匀饱和土层中的反射和透射,并给出了基岩表面和自由表面处反射系数和透射系数的计算表达式.基于理论推导结果,数值分析了平面SV波入射下非均匀饱和土自由场地的动力响应,其中假设饱和土地基的物理力学性质沿土层深度按幂律梯度变化.数值结果表明,平面SV波入射所引起的地面位移与基岩位移之比均随土层厚度和土体的非均匀程度、波的入射角和入射频率的增加而减小,且其竖向位移比的减小更为显著,厚土层对地震波的耗散作用尤为明显.      

基于UPFs的珊瑚岛礁场地地震反应分析方法研究

合理且高效地模拟珊瑚砂非线性动力特性、远场无限地基辐射阻尼以及海域岛礁动水压力的影响是进行珊瑚岛礁抗震安全分析的关键技术问题.以通用有限元软件ANSYS为研究工具,基于UPFs二次开发灵活性的特点,建立了适用于珊瑚砂地基条件下的岛礁场地地震反应分析时域计算模型.通过创建一种新的珊瑚砂等价线性单元描述岛礁场地的非线性动力...     

古田地震水口水电站重力坝强震观测记录的谱结构特征

通过对水口水电站重力坝强震反应台站在古田地震中获取的强震反应观测资料进行信噪比、反应谱和功率谱分析,得到如下结论:①大坝在0.7～15 Hz频率段的振动特性较为可信;②坝基和自由场输入地震动富于高频,峰值加速度反应谱存在较大差异;③坝基输入地震动存在差异性,建议今后此类大坝抗震设计时考虑多点地震动输入;④单个卓越频率携带的能量对反应谱影响不大,反应谱是和输入地震动总能量相关的;⑤坝体刚度较大,此次地震中还处于线弹性状态。初步了解了强震记录的地震动特性和大坝结构的抗震性能,对认识水库地震近场地震动特性和重力坝地震反应有一定的参考意义。     

Results and analysis of soil-structure interaction experiments performed in the centrifuge

In this paper a centrifuge model that is capable of realistically representing soil-structure systems subjected to earthquake-like excitation is used to create a data pool which demonstrates the influence of (i) the frequencies of the structure, (ii) the foundation embedment and (iii) the foundation shape on radiation damping and soil-structure interaction effects for a structure on a semi-infinite soil layer over bedrock. The centrifuge model used in this study was developed and validated by the authors in an earlier publication,¹ and employs an internal method of earthquake simulation, and the clay-like material, Duxseal, to absorb wave reflections at the boundary of the soil sample. The results of the experimental study are used to compute damping and stiffness values of a two-degree-of-freedom piecewise-linear numerical model of the soil-structure systems. The experimental parameter values are then compared to the values computed using classical text book formulae. The analysis demonstrates the value of the experimental data in validating and developing soil-structure interaction theory, and confirms the accuracy of classical text book formulae in the linear range.     

Effect of underground cavities on surface earthquake ground motion under SH wave propagation

A theoretical approach is presented to study the antiplane seismic response of underground structures, subjected to the incidence of both plane and cylindrical waves. The structure is assumed to be a circular inclusion embedded in a homogenous, isotropic and linear visco‐elastic halfspace. The inclusion may consist either of a cavity, with or without a ring‐shaped boundary, or it may be filled in with a linear‐elastic material, without loss of generality. The analytical solution is obtained using expansions of wave functions in terms of Bessel and Hankel functions, relying on the technique of images and the use of Graf's addition theorem to enforce the boundary conditions. The effects of underground cavities on surface earthquake ground motion are studied as a function of the size of the cavity, its embedment depth, the frequency content of the excitation, the incidence angle and the distance from the axis of symmetry of the cavity itself. A simple application of Rayleigh's method allows us to verify that the ground surface response is dominated by the fundamental vibration mode of the portion of soil between the cavity and ground surface itself, in the frequency range of interest for engineering purposes. A simple relationship to estimate the fundamental natural frequency as a function of the embedment depth of the cavity is given. Finally, amplification factors on response spectra are obtained, to provide a practical insight into the effect of an underground cavity on surface ground motion during real earthquakes. Copyright © 2009 John Wiley & Sons, Ltd.     

An indirect boundary element method applied to simulate the seismic response of alluvial valleys for incident P,S and Rayleigh waves

A boundary integral formulation is presented and applied to model the ground motion on alluvial valleys under incident P, S and Rayleigh waves. It is based on integral representations for the diffracted and the refracted elastic waves using single-layer boundary sources. This approach is called indirect BEM in the literature as the sources' strengths should be obtained as an intermediate step. Boundary conditions lead to a system of integral equations for boundary sources. A discretization scheme based on the numerical and analytical integration of exact Green's functions for displacements and tractions is used. Various examples are given for two-dimensional problems of diffraction of elastic waves by soft elastic inclusion models of alluvial deposits in an elastic half-space. Results are displayed in both frequency and time domains. These results show the significant influence of locally generated surface waves in seismic response and suggest approximations of practical interest. For shallow alluvial valleys the response and its resonant frequencies are controlled by a coupling mechanism that involves both the simple one-dimensional shear beam model and the propagation of surface waves.     

大型河谷场地地震动特性研究

采用有限差分方法,通过算例研究了大型河谷场地地震动特性分析中的人工边界的选取方法,对比分析了不同人工边界的选取对数值模拟结果的影响,确定了散射场地震反应分析输入边界的地震动输入方法,认为在进行有限差分动力计算时,模型两侧施加自由场边界的模拟效果要优于两侧施加粘性边界的模拟效果;同时,对FLAC计算软件进行二次开发,对2个地形差异较大的河谷场地,采用线性和非线性摩尔-库仑模型进行了地震反应对比分析,研究了河谷场地地震动幅值和频谱特性随地形变化的规律。模拟数据表明,河谷场地谷底处地震动基本无放大作用,地势凸起处放大作用则较为明显;当考虑土体非线性时,随着地震动强度的增加,放大作用逐渐减小;谱分析结果表明,地势凸起处受高频地震动的影响显著,而地势平坦的谷底则受低频地震动的影响显著。     

On the effective seismic input for non-linear soil-structure interaction systems

Two equivalent semi-discrete formulations are presented for the problem of the transient response of soil-structure interaction systems to seismic excitation, considering linear behaviour of the soil material and arbitrary non-linear structural properties. One formulation results in a direct method of analysis in which the motion in the structure and the entire soil medium, rendered finite by an artificial absorbing boundary, is determined simultaneously. The other represents a substructuring technique in which the structure and the soil are analysed separately. The forces induced in the discretized system by the incident seismic motion are obtained as part of the general formulation by using the free-field motion of the unaltered soil as the earthquake input. It is shown that these forces act within the soil region in the direct method, but only on the soil-structure interface in the substructure formulation. Both sets of forces, however, involve only the displacements and tractions acting on the fictitious surface in the unaltered (linear) soil which coincides with the soil-structure interface of the complete system. It is shown, further, that the free-field displacements alone define a minimal set of data for evaluating the seismic response of the structure, since the tractions and displacements on that surface are interrelated. In practice, the minimal set must be obtained by extrapolating the available information, as the free-field ground motion at a site is usually specified at a single reference point.