基于子结构法的成层场地中沉积河谷的散射分析

土层场地对地震动的影响明显,研究成层场地局部复杂地形的波动场特性,对于抗震设防具有重要意义。基于土-结构相互作用理论,含有局部复杂地形的成层场地,可分解成广义结构（即近场复杂地形及其周围土体）和具有规则开挖边界的成层无限地基（即远场）。远场的格林函数可通过精细积分求解频域-波数域内的对偶波动方程获得,利用傅里叶逆变换得到频域-空间域内的柔度矩阵。近场可采用比例边界有限元进行模拟,通过高性能连分式的传递边界求解动刚度。自由场波动可通过传递矩阵法求得,将动刚度和自由场波动位移代入控制方程即可得到散射场的动力响应。数值算例验证了方法的准确性,并利用提出的方法讨论含有软夹层场地局部复杂场地的波动场特性。     

On numerical computation of impedance functions for rigid soil-structure interfaces embedded in heterogeneous half-spaces

We present fully-discrete procedures for computing the impedance functions of rigid massless soil-structure interfaces that are embedded in arbitrarily heterogeneous half-spaces. The finite element method (FEM) is used for obtaining the wave responses of (visco-)elastic half-spaces truncated by Perfectly Matched Layers (PMLs), which provide the wave absorbing boundary conditions. The devised FEM-PML approach is verified in both time and frequency domains by using various benchmark solutions. Requirements on the prescribed input excitations for obtaining accurate impedances in the time domain as well as the relative computational cost of time- and frequency domain solutions are investigated. Accuracy of the implemented PMLs in extracting the impedance functions is also examined in comparison to Lysmer–Kuhlemeyer dashpots; and it was found that this simplified boundary treatment is generally inadequate. The utility of the proposed method is demonstrated by extracting the impedance matrix of rectangular and circular voids embedded in a linearly stiffening half-space. Impedance functions for such complex soil-structure systems are shown to be highly coupled and frequency-dependent due to wave reflections and interference caused by the soil heterogeneity and interface geometry. Fully discrete approaches, such as those proposed herein, are necessary to adequately capture these effects.     

A practical and efficient numerical scheme for the analysis of steady state unconfined seepage flows

The scaled boundary finite‐element method (SBFEM), a novel semi‐analytical technique, is applied to the analysis of the confined and unconfined seepage flow. This method combines the advantages of the finite‐element method and the boundary element method. In this method, only the boundary of the domain is discretized; no fundamental solution is required, and singularity problems can be modeled rigorously. Anisotropic and nonhomogeneous materials satisfying similarity are modeled without additional efforts. In this paper, SBFE equations and solution procedures for the analysis of seepage flow are outlined. The accuracy of the proposed method in modeling singularity problems is demonstrated by analyzing seepage flow under a concrete dam with a cutoff at heel. As only the boundary is discretized, the variable mesh technique is advisable for modeling unconfined seepage analyses. The accuracy, effectiveness, and efficiency of the method are demonstrated by modeling several unconfined seepage flow problems. Copyright © 2011 John Wiley & Sons, Ltd.     

A nonlinear HHT-α method with elastic–plastic soil–structure interaction in a coupled SBFEM/FEM approach

This article presents a new method for the calculation of elastic–plastic building ground deformations and elastic–plastic building ground failure including wave propagation in the ground. The presented procedure is a hybrid method, based on several common calculation methods. Included is a nonlinear calculation with the finite element method (FEM), a nonlinear HHT-alpha method with full Newton–Raphson iteration and the scaled boundary finite element method (SBFEM). The presented method can be used as a tool for the accurate calculation of building ground deformations and the stability of the subsoil with included dynamic loading.     

局部复杂场地Rayleigh波传播特性分析的SBFEM模型研究

基于土-结构动力相互作用（SSI）的原理,将求解包含局部复杂地形的散射场问题,分解为规则开挖场地的动刚度、包含复杂构造散射体的动刚度和自由场波动输入3个问题的求解,并依据子结构控制方程,提出了一种高精度的比例边界有限（SBFEM）计算模型。通过计算半圆形河谷和半圆形沉积河谷两个算例,验证了SBFEM计算模型的正确性。较为系统地讨论了矩形峡谷、垂直断层和孤凸山体等常见复杂地形对Rayleigh波场的散射作用。结果表明：峡谷深度、软弱断层破碎带和孤凸山体对Rayleigh波都有一定阻挡作用;峡谷宽度、断层宽度和硬性断层填充物等因素对Rayleigh波散射作用影响不明显。SBFEM分析模型可应用于隔震沟、隔震帷幕设计,对大坝、桥梁等工程的选址提供有益的建议。     

ABAQUS动力无限元人工边界研究

针对动力场天然无限地基的数值模拟与地震波输入问题进行了一些有意义的研究,评述了现有动力计算常用无限元的优缺点,详细阐述了ABAQUS无限元理论体系框架,并加以改进,提出一种考虑外域地震动影响的ABAQUS动力无限元人工边界。采用等效边界力的叠加原理,对入射波和散射波分开处理,视入射波和散射波在边界上互不影响,将输入地震动转化为作用于有限元无限元交界面上的等效应力的方法来解决外源波的入射问题。算例验证结果表明：内源振动和固定边界会出现失真和扰动现象,同时该计算结果与黏弹性边界的计算结果对比可知,该方法对外行散射波的过滤作用优于黏弹性边界。因此,改进的ABAQUS动力无限元人工边界理论方法有效且具有一定的稳定性。     

Finite element modeling of wave propagation problems in multilayered soils resting on a rigid base

A new finite element scheme is proposed, in this paper, for solving two-dimensional wave propagation problems in multilayered soils resting on a rigid base. The multilayered soils are treated as multiple horizontal layers of lateral infinite extension in geometry. Since these horizontal layers can be truncated by two artificially truncated vertical boundaries, two high-order artificial boundary conditions are applied for propagating the incoming waves from the interior domain into the far field of the system. Both the semi-analytical method and the truncated boundary migration procedure are used to derive the high-order artificial boundary conditions, which are comprised of a physically meaningful dashpot and a generalized energy absorber. The main advantage of using the proposed finite element scheme is that the derived artificial boundary condition can be straightforwardly implemented in the finite element analysis, without violating the band/sparse structure of the conventional finite element equation. The related numerical examples have demonstrated that the proposed finite element scheme is of high accuracy in dealing with wave propagation problems in multiple horizontal layers.     

Comparison of transmitting boundaries in dynamic finite element analyses using explicit time integration

Explicit time integration schemes provide an efficient solution to non-linear dynamic finite element analyses of geotechnical problems especially when high frequency response is important. Such explicit time integration schemes require one of two distinct transmitting boundary formulations to overcome the problem of radiation damping. These are the superposition boundary approach, which involves the cancellation of the reflected waves by combining the solutions of two different boundary conditions and the viscous boundary approach, which involves the absorption of incident wave energy by frequency independent viscous dashpots. The theoretical justification of these two approaches and their means of implementation are reviewed. The solutions obtained using the two different boundary approaches to the problem of a rigid massless circular footing vibrating on an elastic half-space are compared with an independent theoretical solution. The performance of the boundaries for problems involving step loading is also examined and the implications for any loading pattern with a non-zero time average are discussed.     

Prediction of elastic settlement of rectangular raft foundation — a coupled FE–BE approach

A numerical method of analysis is proposed for computation of the elastic settlement of raft foundations using a FEM–BEM coupling technique. The structural model adopted for the raft is based on an isoparametric plate bending finite element and the raft–soil interface is idealized by boundary elements. Mindlin's half-space solution is used as a fundamental solution to find the soil flexibility matrix and consequently the soil stiffness matrix. Transformation of boundary element matrices are carried out to make it compatible for coupling with plate stiffness matrix obtained from the finite element method. This method is very efficient and attractive in the sense that it can be used for rafts of any geometry in terms of thickness as well as shape and loading. Depth of embedment of the raft can also be considered in the analysis. Copyright © 1999 John Wiley & Sons, Ltd.     

An approximate spring–dashpot artificial boundary for transient wave analysis of fluid-saturated porous media

Practical civil engineering problems are usually formulated in an infinite half-space domain, and a selected finite domain is required to analyze the dynamic responses of a fluid-saturated porous medium by the finite element method (FEM). Devising a method to deal with the boundaries of the finite domain is the key issue for this open system. In this paper, a two-dimensional spring–dashpot artificial boundary (SDAB) for transient analysis in a fluid-saturated porous media is developed. Based on Biot’s dynamic theory of fluid-saturated porous media, the normal and tangential boundary stress formulae are deduced for out-going cylindrical body waves. The boundary stress is proportional to displacement and velocity, thus continuously distributed dashpots and springs can be placed on the artificial boundaries in the normal and tangential directions to simulate the energy absorption of the infinite media outside of the finite domain for the interior distributed source problems. In this paper, the input seismic motion can be realized by applying an equivalent load on the SDAB for the seismic scattering problems of exterior distributed sources. Numerical examples are given and the analyzed results show that the SDAB and the method of wave motion input have good stability and acceptable accuracy.     

Coupled meshfree and fractal finite element method for unbounded problems

This paper presents a coupling technique for integrating the element-free Galerkin method (EFGM) with the fractal finite element method (FFEM) to analyze unbounded problems in the half-space. FFEM is adopted to model the far field of an unbounded domain and EFGM is used in the near field. In the transition region interface elements are employed. The shape functions of interface elements which comprise both the element-free Galerkin and the finite element shape functions, satisfy the consistency condition thus ensuring convergence of the proposed coupled EFGM–FFEM. The proposed method combines the best features of EFGM and FFEM, in the sense that no structured mesh or special enriched basis functions are necessary. The numerical results show that the proposed method performs extremely well converging rapidly to the analytical solution. Also a parametric study is carried out to examine the effects of the integration order, the similarity ratio, the weight function, the scaling parameter and the number of transformation terms, on the quality of the numerical solutions.     

Boundary element approach to the dynamic stiffness functions of circular foundations

A boundary element approach for time harmonic axisymmetric problems using the complete space point load fundamental solution is presented. The fundamental solution is integrated numerically along the azimuthal co-ordinate of each axisymmetric element. To increase the accuracy of the numerical integration a simple co-ordinate transformation is proposed. The approach is applied to the computation of the dynamic stiffness functions of rigid circular foundations on layered viscoelastic soils. Three different sites are considered: a uniform half-space, a soil layer on a half-space, and a soil consisting of four horizontal layers and a compliant half-space. The numerical results obtained by the proposed approach for surface circular foundations are very close to corresponding published results obtained by different procedures.     

地震荷载作用下大坂盾构隧道管片配筋校核

将盾构隧道管片视为分段曲梁、管片之间为铰接关系、管片与地层之间为非线性接触关系,边界采用无限元吸收地震波的反向与折射,地层服从Mohr-Coulomb弹塑性本构模型,地震波采用EI-Centro波,采用地层-结构法,应用隐式动力有限元分析地震荷载作用下管片轴力与弯矩分布情况。根据有限元计算结果,结合结构力学方法,校核管片配筋是否合理。与荷载结构法相比,其方法可模拟围岩与管片、管片与管片之间的相互关系,并可考虑锚杆等支护结构对管片受力的影响,因而,计算结果更为精确。研究为盾构隧道管片配筋优化设计,降低隧道支护成本,提供了一种管片配筋校核计算方法与设计依据。     

Wave scattering and diffraction of subsurface cavities in layered half-space for incident SV-P and SH waves

A novel approach is proposed to deal with the problem of wave scattering and diffraction of subsurface cavities embedded in stratified half-space. The subsurface cavity with complex surroundings is treated as a substructure. The continuity condition at the interface between the substructure and the far field of stratified half-space is maintained by applying a free-field approach. As the boundary of the free-field ground is regular, the construction of the dynamic matrices and the evaluation of the wave input on the interface become considerably easier. Based on the previous work with some improvement, a novel approach for evaluation of Green's functions in stratified half-space is presented. The wave equation is decoupled into the one for SV-P wave components and the other one for SH wave component. The precise integration technique ensures high accuracy of the solution of wave equations. The layer merging technique and the dual form equation make it possible to obtain Green's function in closed-form solution of matrix equations. Numerical examples validate accuracy and efficiency of the proposed approach.     

P波作用下跨断层隧道轴线地震响应分析

基于近场波动有限元方法并结合黏弹性人工边界条件,针对3D断层场地,通过求解等效二维场地地震响应,从而获得3D场地地震动输入的自由场响应,并将自由场响应转化为3D模型边界面上的等效节点力,从而建立含断层3D场地P波入射的倾斜输入方法。自由场算例验证所提方法具有较好精度,进而基于建立的输入方法,开展了跨断层隧道地震响应的数值模拟研究。数值模拟结果表明:在P波作用下,隧道跨断层部位处于拉、压、剪切的复杂受力状态,且断层处的隧道衬砌地震响应明显大于其他部位的地震响应;围岩的力学性质与断层的力学性质相差越大,断层处衬砌的地震响应放大越明显;断层处衬砌地震响应随断层深度的增加而增加。另外,跨断层隧道的地震响应受P波入射角度的影响较大,随P波入射角度的增加,断层处隧道衬砌的轴力、弯矩先增加后减小,而剪力具有逐渐减小的规律。     

地震分析中人工边界处理与地震动输入方法研究

基于柱面波波动方程,推导建立了适用于土-结构地震动力相互作用分析的地震动输入和人工边界的处理方法.其中,地震动的输入是通过在人工边界上施加等效节点力来实现的,等效节点力的大小与入射地震波波速成正比;而人工边界的处理方法使得人工边界条件不仅在时间上是局部的,而且在空间上也是局部的.这种处理方法简单、有效,物理意义清晰,且很容易在有限元法中实现,结合Newmark时间积分是无条件稳定的.为了验证方法的有效性和精度,给出了两个算例,分别用于检验人工边界条件的性能以及地震动输入方法的正确性.算例分析结果表明,所提出的方法是十分有效的.     

有限压缩层地基上筏基沉降的边界元分析

基于有限压缩层地基模型和Reissner板的边界积分方程,考虑柱荷载作用区域的形状和筏板的横向剪切变形效应,建立了有限压缩层地基上厚筏基础与地基相互作用分析的边界元方程和系统的数值方法。对于弹性半空间模型,其柔度方程可视为有限压缩层模型的特殊情况。把筏板作为自由边界条件处理,被剖分为一系列三角形或矩形网格,假设基底反力在网格内均匀分布,以便与现有的地基沉降计算模式相一致。计算表明,虽然基底反力在内部网格相接处不连续,但并不影响计算结果,反而消除了边界基底压力计算值过大的现象。将该方法与其他方法的计算结果进行比较,显示了该方法的有效性。计算结果表明,对于实际复杂的筏板基础,无需划分太多单元即可得到较高的计算精度     

波的传播及其与介质间断相互作用的数值模拟研究

本文用有限差分法模拟研究二维空间中波的传播、散射及其与介质间断的相互作用。文章第一部分运用所模拟的纵波波源讨论了半无限空间表面台阶地形或槽形裂隙所造成的波的散射。这些间断可模拟具有波长量级的悬崖或狭谷。第二部分用差分法与扰动法结合而成的混合法处理半无限空间表面有另一种物质构成的表面浅障碍时波的传播,以模拟山脊、堑沟等地形的影响。此混合法也被用于讨论半无限空间内含有水平夹层或孔隙的问题。所得研究结果对波的传播及其与介质间断的相互作用的理论研究和地震波破坏效应的分析和隔离等应用研究有一定的意义。     

水平-摇摆耦合振动波阻板主动隔振三维分析

为研究三维地基中波阻板（简称WIB）的隔振效果,建立了以薄层法基本解答作为格林函数的半解析边界单元法。对水平-摇摆耦合激振作用下WIB隔振系统的主动隔振进行了计算,详细分析了WIB参数对隔振效果的影响。研究表明,WIB可以有效地隔振;增加WIB的厚度、提高WIB的弹性模量、减小WIB 的埋深可以有效提高隔振效果,而增加WIB的直径可以显著降低WIB范围内地表位移振幅,但对WIB范围外的地表位移影响较为复杂。