Quasi-static deformation due to two-dimensional seismic sources embedded in an elastic half-space in welded contact with a poroelastic half-space

The Biot linearized theory of fluid saturated porous materials is used to study the plane strain deformation of a two-phase medium consisting of a homogeneous, isotropic, poroelastic half-space in welded contact with a homogeneous, isotropic, perfectly elastic half-space caused by a two-dimensional source in the elastic half-space. The integral expressions for the displacements and stresses in the two half-spaces in welded contact are obtained from the corresponding expressions for an unbounded elastic medium by applying suitable boundary conditions at the interface. The case of a long dip-slip fault is discussed in detail. The integrals for this source are solved analytically for two limiting cases: (i) undrained conditions in the high frequency limit, and (ii) steady state drained conditions as the frequency approaches zero. It has been verified that the solution for the drained case (ω → 0) coincides with the known elastic solution. The drained and undrained displacements and stresses are compared graphically. Diffusion of the pore pressure with time is also studied.     

Fuzzy finite element analysis of a foundation on an elastic soil medium

One of the important factors that lead to errors in settlement predicitions is the degree of precision in obtaining the soil parameters. Most mathematical methods for reliability modelling offered to date in the area of geomechanics are based on the classical probabilistic approach, in which soil properties are treated as random variables. In this paper, a model based on the theory of fuzzy sets is presented to take account of the uncertainty in the soil behaviour. This proposed method considers the elastic modulus and Poisson's ratio as two fuzzy numbers in the elastic matrix. An example is given to show the possibility distributions of the displacements and the stresses at some locations in the soil medium. By means of a fuzzy inference scheme, the total possibility distributions or total membership functions of the finite element results may be obtained by considering the estimated error resulting from the mesh discretization.     

Influence functions for stress and displacement discontinuity elements in an anisotropic medium

Influence functions, that permit us to determine stresses and displacements at an arbitrary point in an infinite, homogeneous, linear elastic, anisotropic medium due to different three-dimensional (3-D) stress or displacement discontinuities distributed on infinite, flat, band-type elements, are presented. Any straight-line segment on the band, which is perpendicular to its infinite side, has the same distribution of the discontinuities. Along with the functions, their Taylor series approximations are also provided. The last can be useful to analyse stresses and displacements at points distant from the elements. The functions allow us to avoid procedures of numerical integration in the Indirect Boundary Element Method and/or the Displacement Discontinuity Method computer codes that are able to solve complete plane-strain problems with 3-D boundary conditions for an elastic, anisotropic medium. © 1997 by John Wiley & Sons, Ltd.     

Wave-induced displacements in seafloor sands

The stability and deformation of the seafloor in response to storm wave loading is an important consideration in the design of pipelines, anchors, and offshore structures such as gravity-and pile-supported platforms. Two quite different approaches have been proposed for investigating the effects of storm waves on seafloor sands. These approaches derive from quite different considerations of soil-water interaction. In one of these approaches (exemplified by Yamamoto and Madsen) the seafloor is visualized to be a porous elastic or Coulomb-damped porous nonlinear medium. In this approach, the response of the seafloor, i.e. the displacements, stresses and porewater pressures, are computed for a given wave. The response is transient and is, therefore, related to the largest wave in the storm wave group. On the other hand, in the approach proposed by Seed and Rahman, only the development of residual porewater pressures (for estimating liquefaction potential) due to the storm wave group are considered. The effects of transient porewater pressures and stresses on the seafloor sands are ignored. The residual porewater pressures are computed using properties based on initial effective stresses through a porewater generation-dissipation model. The application of these techniques to problems of offshore structures such as pile-supported platforms has led to considerable uncertainty and confusion. In this paper, both of these approaches are combined to investigate seafloor response to a storm wave group. For typical storm wave conditions that exist in the North Sea, it appears that the inclusion of damping, inertia and anisotropic permeabilities in the study is not important relative to ocean sands. Parametric studies show that, for a given wave, the thickness and the stiffness properties of the soil deposit dictate the response of the deposit. The displacements in the deposit are, therefore, affected by the generation of residual porewater pressures.     

Displacement and stress distributions under a uniform inclined rectangular load on a cross‐anisotropic geomaterial

In practical engineering, an applied rectangular area load is not often horizontally or vertically distributed but is frequently inclined at a certain angle with respect to the horizontal and vertical axes. Thus, the solutions of displacements and stresses due to such a load are essential to the design of foundations. This article yields the analytical solutions of displacements and stresses subjected to a uniform rectangular load that inclines with respect to the horizontal and vertical axes, resting on the surface of a cross‐anisotropic geomaterial. The planes of cross‐anisotropy are assumed to be parallel to the horizontal ground surface. The procedures to derive the solutions can be integrated the modified point load solutions, which are represented by several displacement and stresses elementary functions. Then, upon integrations, the displacement and stress integral functions resulting from a uniform inclined rectangular load for (1) the displacements at any depth, (2) the surface displacements, (3) the average displacements in a given layer, (4) the stresses at any depth, and (5) the average stresses in a given layer are yielded. The proposed solutions are clear and concise, and they can be employed to construct a series of calculation charts. In addition, the present solutions clarify the load inclinations, the dimensions of a loaded rectangle, and the analyzed depths, and the type and degree of geomaterial anisotropy profoundly affect the displacements and stresses in a cross‐anisotropic medium. Parametric results show that the load inclination factor should be considered when an inclined rectangular load uniformly distributed on the cross‐anisotropic material. Copyright © 2008 John Wiley & Sons, Ltd.     

水布垭特大断面导流洞围岩稳定性分析

湖北清江水布垭导流洞采用矿山法施工,为确保施工安全,以三维有限元对该隧洞处于断层破碎带的进口渐变段的施工过程进行了计算模拟,并根据计算结果确定了合理的施工方案.根据地下结构设计理论中的地层结构法,选用连续体计算模型,对围岩材料和混凝土材料的本构关系分别作理想弹塑性应力-应变关系和线弹性本构关系的简化,并采用摩尔-库仑屈服准则,分析了初期支护结构的受力特点、围岩塑性区的深度和范围以及洞周围岩位移的发展规律.结果显示,围岩塑性区深度较大,而且初期支护结构承受的最大拉、压应力均略大于混凝土的强度极限,指出该区段初期支护结构在采用钢筋格栅喷混凝土的基础上,应加强长、短系统锚杆进行围岩加固.     

Scattering of elastic wave by a partially sealed cylindrical shell embedded in poroelastic medium

Scattering of an elastic wave by a cylindrical shell embedded in poroelastic medium is investigated theoretically with the assumption that the shell material is also a porous elastic medium. The porous medium is modellized via Biot's theory and the scattering by cylindrical shell is expressed by the definition of scattering matrix. The normal mode expansion technique is employed for analysing the scattering field, and the asymptotic solutions of displacements, stresses and pore pressure are derived. Two limiting cases‐scattering by a poroelastic cylinder in Biot medium and a elastic cylindrical shell in elastic medium are obtained from the general solutions. The dispersion curves of displacement amplitude at the interface of shell and medium is compared with the case of elastic shell. The scattering amplitude associated with the fast, slow and transverse waves are identified by numerical simulation. Furthermore, the influence of the poroelastic property of shell material on scattering amplitude is analysed. Copyright © 2005 John Wiley & Sons, Ltd.     

非线性接触力学模型在地基-基础相互作用弹性分析中的应用

针对多体相互作用体系的非连续变形分析问题和接触问题,采用Mohr-Coulomb屈服准则和关联流动法则以及接触界面上的非线性应力分布模式,考虑接触界面特性提出了非线性接触力元模型,以结点位移和界面相互接触应力同时作为独立未知变量,建立了离散系统的总体控制方程.进而,通过数值求解能够直接确定变形体内的应力与变形、界面上的接触应力与离散体的位移与运动.将这种以接触力元为基础的多体系统分析方法具体应用于基础与地基相互作用分析,通过数值计算与分析探讨了地基与基础的相对刚度、荷载大小及其偏心距、地基与基础间界面力学参数对接触界面的应力分布和地基变形的影响,所得结果为工程中考虑基础与地基相互作用影响的设计与分析提供了参考依据.     

Static deformation of an orthotropic multilayered elastic half-space by two-dimensional surface loads

The transfer matrix approach is used to solve the problem of static deformation of an orthotropic multilayered elastic half-space by two-dimensional surface loads. The general problem is decoupled into two independent problems. The antiplane strain problem and the plane strain problem are considered in detail. Integral expressions for displacements and stresses at any point of the medium due to a normal line load and a shear line load, acting parallel to a symmetry axis, are obtained. In the case of a uniform half-space, closed form analytic expressions for displacements and stresses are derived. The procedure developed is quite easy and convenient for numerical computations.     

广蓄电站厂房的三维位移反分析

本文采用参数同步变化法对广蓄电站厂房进行全三维的位移反分析,以确定厂房区的地应力和岩体的材料参数。参数同步变化法基于测点埋设时和测试时的几何组构分别进行地应力场和材料参数同步变化的三维计算,取两次计算的位移差值作为优化目标函数的相关变量,与实测值进行比较。这种方法去除了人为的歧意性,适用于各种复杂的地质条件。     

Boundary element analysis of transversely isotropic bi-material halfspaces with inclined planes of isotropy and interfaces

In this paper, a single-region boundary element method (BEM) is presented for the analysis of transversely isotropic bi-material halfspaces with arbitrarily inclined planes of isotropy and material interfaces. The proposed BEM uses the fundamental solution of a transversely isotropic bi-material fullspace and five boundary element techniques. Infinite boundary elements are introduced to consider the far-fields of a transversely isotropic bi-material halfspace. The effective integration methods are proposed for dealing with various integrals in the discretized boundary integral equation. The stresses at internal points are obtained using the coordinate transformation of kernel functions, and the stresses on the boundary surface are calculated using an improved traction recovered method. Numerical verifications of displacements and stresses for a benchmark problem are conducted, and excellent agreement with previously published results is obtained. Numerical examples are presented to illustrate the influence of non-horizontal or horizontal planes of isotropy in bi-material halfspaces on the displacements and stresses induced by the tractions on the horizontal boundary surface. Results reveal that the elastic fields vary clearly with the dip angle of the isotropic plane and the stresses across the bi-material interface are closely related to the ratios of the elastic parameters of the bi-material.     

Analytical stress and displacement due to twin tunneling in an elastic semi‐infinite ground subjected to surcharge loads

The construction of twin tunnels at shallow depth has become increasingly common in urban areas. In general, twin tunnels are usually near each other, in which the interaction between tunnels is too significant to be ignored on their stability. The equivalent arbitrarily distributed loads imposed on ground surface were considered in this study, and a new analytical approach was provided to efficiently predict the elastic stresses and displacements around the twin tunnels. The interaction between 2 tunnels of different radii with various arrangements was taken into account in the analysis. We used the Schwartz alternating method in this study to reduce the twin‐tunnel problem to a series of problems where only 1 tunnel was contained in half‐plane. The convergent and highly accurate analytical solutions were achieved by superposing the solutions of the reduced single‐tunnel problems. The analytical solutions were then verified by the good agreement between analytical and numerical results. Furthermore, by the comparison on initial plastic zone and surface settlement between analytical solution and numerical/measured results of elastoplastic cases, it was proven that the analytical solution can accurately predict the initial plastic zone and its propagation direction and can qualitatively provide the reliable ground settlements. A parametric study was finally performed to investigate the influence of locations of surcharge load and the tunnel arrangement on the ground stresses and displacements. The new solution proposed in this study provides an insight into the interaction of shallow twin tunnels under surcharge loads, and it can be used as an alternative approach for the preliminary design of future shallow tunnels excavated in rock or medium/stiff clay.     

Strain‐softening analysis of a spherical cavity

This paper studies the excavation of a spherical cavity subjected to hydrostatic initial stresses in the infinite homogeneous and isotropic rock mass with strain‐softening Mohr–Coulomb (M‐C) and Hoek–Brown (H‐B) behaviors. Numerical solutions of the spherical cavity are obtained and the application to determining stress–strain curve of strain‐softening M‐C and H‐B rock mass is studied. A closed‐form solution for the elastic–brittle–plastic medium is introduced first, and then a numerical procedure that simplifies the strain‐softening process into a series of brittle–plastic ones is presented. The approach is validated against the facts that the strain‐softening process evolves into a brittle–plastic one when the softening slope is very steep, whereas it evolves into an elasto‐plastic one when the softening slope approaches zero. Numerical solutions for the prediction of displacements and stresses around the spherical cavity in the strain‐softening M‐C and H‐B rock mass are presented. On the basis of the analysis of the spherical cavity in strain‐softening rock mass, the stress–strain relationship at an infinitesimal cube around the cavity is obtained and discussed with different evolution laws for the strength parameters considered. Copyright © 2011 John Wiley & Sons, Ltd.     

Block displacement fields in the Altai-Sayan region and effective rheologic parameters of the Earth’s crust

The paper is focused on recent displacement rates in the Altai-Sayan region, obtained by hydroleveling, leveling, and satellite geodesy. Effective elastic moduli and viscosity parameters of the crust are used in the modeling of coseismic and tectonic processes. The elastic moduli are determined from measurements of periodic vertical displacements during seasonal loadings of the Sayano-Shushenskaya hydropower plant. We present the results of the modeling of coseismic displacements during the earthquakes of 10 February 2011 (M = 6.1) and 27 December 2011 (M = 6.7) in Tuva and West Sayan. The results of GPS determinations for postseismic displacements in the Chuya earthquake zone (Gorny Altai, 27 September 2003, M = 7.5) are analyzed; models for the geologic medium are selected; and its effective viscosity is estimated. The tectonic component of the recent crustal displacements in the Altai-Sayan region is defined.     

2-D deformation analysis of a half-space due to a long dip-slip fault at finite depth

Closed form analytical expressions of stresses and displacements at any field point due to a very long dip-slip fault of finite width buried in a homogeneous, isotropic elastic half-space, are presented. Airy stress function is used to derive the expressions of stresses and displacements which depend on the dip angle and depth of the upper edge of the fault. The effect of dip angle and depth of the upper edge of the fault on stresses and displacements is studied numerically and the results obtained are presented graphically. Contour maps for stresses and displacements are also presented. The results of Rani and Singh (1992b) and Freund and Barnett (1976) have been reproduced.     

半空间饱和土中圆形衬砌对弹性压缩波的散射

借助于Biot 波动理论和弹性波的传播理论,采用复变函数和多级坐标法,对半空间饱和土中圆形衬砌结构对弹性稳态压缩波的散射问题进行求解和分析。利用一个半径很大的圆弧来逼近半空间直边界,将待解问题转化为稳态弹性压缩波在一个大圆孔和一个弹性衬砌结构的散射问题。通过引入势函数,将饱和土的Biot波动方程和衬砌的弹性波动方程解耦成Helmholtz 方程,借助复变函数级数展开便可以预先写出该组Helmholtz方程的通解。然后,通过引用复变量,把饱和土和衬砌结构中的应力、位移及孔压用设定的势函数表示出来,再利用半空间饱和土和衬砌结构的连续性条件和近似直边界的圆弧边界和衬砌内边界的边界条件求解出该组势函数的特解。最后,利用势函数的特解,得到饱和土中的位移,应力和孔压及衬砌结构的位移和应力;变换不同的参数求解衬砌结构内外边界的动应力和孔压的集中系数,通过对算例结果的分析得出一系列有益的结论。     

Analytical study of ground responses induced by the excavation of quasirectangular tunnels at shallow depths

The construction of quasirectangular tunnels at shallow depths is becoming increasingly common in urban areas to efficiently utilize underground space and reduce the need for backfilling. To clarify the mechanical mechanism of the stresses and displacements around the tunnels, this study proposes analytical solutions that precisely account for quasirectangular tunnel shapes, the ground surface, the tunnel depth, and the ground's elastic/viscoelastic properties. The Schwarz alternating method combined with complex variable theory is employed to derive the elastic solution, and convergent and highly accurate solutions are obtained by superposing the solutions in the alternating iterations. Based on the solution and the extended corresponding principle for the viscoelastic problem, the time-dependent analytical solutions for the displacement are obtained for the ground assuming any viscoelastic model. The analytical solutions agree well with the finite element method (FEM) numerical results for models that are completely consistent, and qualitatively agree with field data. Furthermore, based on the stress solution combined with the Mohr-Coulomb failure criterion, the predicted initial plastic zone and propagation directions around the tunnels are qualitatively consistent with those determined by the limit analysis. A parametric study is performed to investigate the influences of the rectangular/quasirectangular tunnel shape, burial depth, and supporting pressure on the ground stresses and displacements.     

爆炸荷载作用下深埋圆形隧洞的饱和黏弹性土-衬砌系统动力响应

假定土骨架服从标准线性固体黏弹性本构关系,研究了深埋圆形隧洞的饱和黏弹性土-弹性衬砌耦合系统在轴对称爆炸作用下的瞬态动力响应。首先,基于饱和土的Biot模型和衬砌的弹性理论,通过引入势函数和Laplace变换,利用弹性衬砌和饱和黏弹性土界面处的连续性条件以及边界条件,得到饱和黏弹性土体和弹性衬砌位移、应力和孔隙水压力等在Laplace变换域中的解析解。其次,利用Laplace数值Crump逆变换得到耦合系统在时间域的动力响应,数值分析了不同土体模型下土体-衬砌耦合系统的径向位移和环向应力以及土体孔隙水压力等。结果表明：对不同土体模型的土体-衬砌耦合系统,其在爆炸载荷作用下的动力响应性态基本一致,但动力响应的振动周期和幅值等具有明显的差异。同时,对于饱和黏弹性土-弹性衬砌系统,土体黏性参数对土体径向位移和孔隙水压力有明显的影响,但对土体环向应力影响较小。     

Elastic solutions for a transversely isotropic half-space subjected to a point load

We rederive and present the complete closed-form solutions of the displacements and stresses subjected to a point load in a transversely isotropic elastic half-space. The half-space is bounded by a horizontal surface, and the plane of transverse isotropy of the medium is parallel to the horizontal surface. The solutions are obtained by superposing the solutions of two infinite spaces, one acting a point load in its interior and the other being free loading. The Fourier and Hankel transforms in a cylindrical co-ordinate system are employed for deriving the analytical solutions. These solutions are identical with the Mindlin and Boussinesq solutions if the half-space is homogeneous, linear elastic, and isotropic. Also, the Lekhnitskii solution for a transversely isotropic half-space subjected to a vertical point load on its horizontal surface is one of these solutions. Furthermore, an illustrative example is given to show the effect of degree of rock anisotropy on the vertical surface displacement and vertical stress that are induced by a single vertical concentrated force acting on the surface. The results indicate that the displacement and stress accounted for rock anisotropy are quite different for the displacement and stress calculated from isotropic solutions. © 1998 John Wiley & Sons, Ltd.     

弹性矩形板下横观各向同性多层地基分析

利用弹性矩形板与多层地基表面的竖向位移协调条件与光滑接触条件,由横观各向同性多层地基应力与位移非耦合的传递矩阵解,推导出弹性矩形板下竖向应力和位移的解析解。在此基础上,编制了相应的程序,并进行了数值计算。计算结果表明：矩形板刚度对板底竖向位移及板中心下的竖向应力有着较为显著的影响;板底竖向位移及板中心下的竖向应力随着板刚度的增加而减小,相同荷载作用下横观各向同性地基与均匀各向同性地基模型的计算结果差异较大,实际工程中很有必要采用更符合土体性质的横观各向同性地基模型