The axial loading of a rigid circular anchor plate embedded in an elastic half-space

The present paper examines the axisymmetric problem related to the loading of a rigid circular anchor plate which is embedded in bonded contact with an isotropic elastic half-space. A Hankel transform development of the governing equations is used to reduce the associated mixed boundary value problem to a set of coupled Fredholm integral equations of the second kind. These equatons are solved in a numerical fashion to generate results of engineering interest. In particular, the results indicate the influence of the depth of embedment on the axial stiffness of the rigid anchor plate.     

The influence of a boundary fracture on the elastic stiffness of a deeply embedded anchor plate

The problem of the axial loading of a rigid disk-shaped anchor plate embedded in an isotropic elastic medium of infinite extent is examined. At the boundary of the disk anchor plate the elastic medium contains a cracked region of finite extent. The presence of the cracked region decreases the elastic stiffness of the anchor plate. The mathematical formulation of the problem is developed, and a numerical scheme is presented which can be used to solve the resulting coupled integral equations. The numerical technique is used to evaluate the results, which illustrate the manner in which the elastic stiffness of the anchor plate is influenced by the extent of cracking. Similar results are developed for the flaw shearing mode stress intensity factor at the external boundary of the cracked region.     

An indirect boundary integral equation method for poroelasticity

This paper presents an indirect boundary integral equation method for analysis of quasi-static, time-harmonic and transient boundary value problems related to infinite and semi-infinite poroelastic domains. The present analysis is based on Biot's theory for poroelastodynamics with fluid viscous dissipation. The solution to a given boundary value problem is reduced to the determination of intensities of forces and fluid sources applied on an auxiliary surface defined interior to the surface on which the boundary conditions are specified. A coupled set of integral equations is established to determine the intensities of forces and fluid sources applied on the auxiliary surface. The integral equations are solved numerically in the Laplace domain for quasi-static and transient problems, and in the frequency domain for time-harmonic excitations. The kernel functions of the integral equation correspond to appropriate Green's functions for a poroelastic full space or half-space. The convergence and numerical stability of the present scheme are established by considering a number of bench mark problems. The versatility of the present method is demonstrated by studying the quasi-static response of a rigid spheroidal anchor, and time-harmonic and transient response of a rigid semi-circular tunnel.     

Elastic lateral dynamic impedance functions for a rigid cylindrical shell type foundation

Elastic lateral dynamic impedance functions are defined as the ratio of the lateral dynamic force/moment to the corresponding lateral displacement/rotation at the top ending of a foundation at very small strains. Elastic lateral dynamic impedance functions have a defining influence on the natural frequencies of offshore wind turbines supported on cylindrical shell type foundations, such as suction caissons, bucket foundations, and monopiles. This paper considers the coupled horizontal and rocking vibration of a cylindrical shell type foundation embedded in a fully saturated poroelastic seabed in contact with a seawater half‐space. The formulation of the coupled seawater–shell–seabed vibration problem is simplified by treating the shell as a rigid one. The rigid shell vibration problem is approached by the integral equation method using ring‐load Green's functions for a layered seawater‐seabed half‐space. By considering the boundary conditions at the shell–soil interface, the shell vibration problem is reduced to Fredholm integral equations. Through an analysis of the corresponding Cauchy singular equations, the intrinsic singular characteristics of the problem are rendered explicit. With the singularities incorporated into the solution representation, an effective numerical method involving Gauss–Chebyshev method is developed for the governing Fredholm equations. Selected numerical results for the dynamic contact load distributions, displacements of the shell, and lateral dynamic impedance functions are examined for different shell length–radius ratio, poroelastic materials, and frequencies of excitation. Copyright © 2016 John Wiley & Sons, Ltd.     

Rigid cylinder in a transversely isotropic half‐space under lateral loads

A boundary integral equation method is presented for a rigid cylindrical pipe‐pile of finite length embedded in a transversely isotropic half‐space under lateral loads. In the framework of three‐dimensional elastostatics, the complicated soil‐structure interaction problem is shown to be reducible to three coupled Fredholm integral equations. Through an analysis of the associated Cauchy singular kernels, the intrinsic singular characteristics of the radial, angular, and vertical interfacial load transfers are rendered explicit. By means of a complicated numerical procedure, detailed results on the three‐dimensional load–transfer process are provided for benchmark comparison and practical applications. Copyright © 2011 John Wiley & Sons, Ltd.     

A theoretical estimate of the tilt of a surface foundation due to an inclined anchor load applied at the interior of the soil medium

This paper examines the problem of the interaction between a loaded rigid circular foundation located at the surface of an isotropic elastic halfspace and an inclined concentrated anchor load which is located at a finite depth along the axis of symmetry. Such inclined loads can be induced by, for example, anchor regions supporting earth retaining structures. The loaded rigid circular foundation resting in smooth contact with the elastic soil mass experiences a displacement and a tilt due to the action of the inclined anchor load. The magnitude of the rotational settlement is evaluated in exact closed form.     

Forced vertical and horizontal movements of a rectangular rigid foundation on a transversely isotropic half‐space

An analytical investigation of a half‐space containing transversely isotropic material under forced vertical and horizontal displacements applied on a rectangular rigid foundation is presented in this paper. With the goal of a rigorous solution to the shape‐ and rigidity‐ induced singular mixed boundary value problem, the formulation employs scalar potential representation, the Fourier expansion and the Hankel integral transforms method to obtain the surface arbitrary point‐load solution in cylindrical coordinate system. The obtained Green's functions are rewritten in rectangular coordinate system, allowing the response of the half‐space because of an arbitrary distributed load on a rectangular surface area be given in terms of a double integral. The numerical evaluations of stresses are done with the use of an element, which is singular at the edge and the corner of the rectangle. Upon the imposition of the rigidity displacement boundary condition for a rigid foundation and the use of a set of two‐dimensional adaptive‐gradient elements, which can capture the singular behavior in the contact stress effectively, a set of new numerical results are presented to illustrate the effect of transverse isotropy on the foundation response. Copyright © 2012 John Wiley & Sons, Ltd.     

Vertical dynamic response of a rigid foundation embedded in a poroelastic soil layer

A simplified analytical method is presented for the vertical dynamic analysis of a rigid, massive, cylindrical foundation embedded in a poroelastic soil layer. The foundation is subjected to a time‐harmonic vertical loading and is perfectly bonded to the surrounding soil in the vertical direction. The soil underlying the foundation base is represented by a single‐layered poroelastic soil based on rigid bedrock while the soil at the side of the foundation is modeled as an independent poroelastic layer composed of a series of infinitesimally thin layers. The behavior of the soil is governed by Biot's poroelastodynamic theory and its governing equations are solved by the use of Hankel integral transform. The contact surface between the foundation base and the soil is smooth and fully permeable. The dynamic interaction problem is solved following standard numerical procedures. The accuracy of the present solution is verified by comparisons with the well‐known solutions obtained from other approaches for both the elastodynamic interaction problem and poroelastodynamic interaction problem. Numerical results for the vertical dynamic impedance and response factor of the foundation are presented to demonstrate the influence of nondimensional frequency of excitation, soil layer thickness, poroelastic material parameters, depth ratio and mass ratio on the dynamic response of a rigid foundation embedded in a poroelastic soil layer. Copyright © 2008 John Wiley & Sons, Ltd.     

The eccentric loading of a rigid circular foundation embedded in an isotropic elastic medium

This paper examines the problem of the eccentric loading of a rigid circular disc-shaped foundation embedded in bonded contact with an istropic elastic medium of infinite extent. The solution of this problem is achieved by using a Hankel integral transform technique which reduces the problem to two sets of dual integral equations. These systems of dual integral equations represent the equations which govern the axisymmetric and asymmetric states of deformation induced by the loaded foundation. Explicit results are derived both for the displacement and rotation of the circular foundation and for the contact stress at the interface.     

The response of a deep rigid anchor due to undrained elastic deformations of the surrounding soil medium

The equations governing the undrained linear elastic behaviour of a saturated soil are formally similar to the equations governing slow of an incompressible Newtonian viscous fluid. This principle of equivalence can then be effectively employed to obtain the load-deflection reiationship for a deep rigid anchor with the shape of a solid of revolution which is embedded in bonded contact with an unbounded incompressible elastic medium. It is found that the load-deflection relationship for the deep rigid anchor can be directly recovered from the expression for the drag induced on an impermeable object with the same size and shape as the anchor, which is appropriately placed in a slow viscous flow region of uniform velocity.     

A contact problem for a poroelastic halfspace containing an embedded inextensible membrane

The paper examines the axisymmetric problem of the indentation of a poroelastic halfspace that is reinforced with an inextensible permeable/impermeable membrane located at a finite depth by a rigid indenter. The constitutive behavior of the poroelastic halfspace is described by the three-dimensional theory of poroelasticity proposed by M.A. Biot. The contact conditions between the indenter and the poroelastic halfspace are varied to accommodate both adhesive/frictionless contact and impermeable/permeable conditions. The formulation of the mixed boundary value problems uses the stress function approaches applicable to semi-infinite domains. Successive applications of Laplace and Hankel integral transforms are used to reduce the mixed boundary value problems to sets of coupled Fredholm integral equations of the second kind. These integral equations are solved using numerical approaches, applicable both for the solution of the systems of coupled equations and for Laplace transform inversion, to examine the time-dependent displacement of the rigid indenter. The analytical-numerical estimates for the time-dependent displacements of the rigid indenter are compared with results obtained using a finite element approach.     

The Reissner–Sagoci problem for a transversely isotropic half‐space

A transversely isotropic linear elastic half‐space, z?0, with the isotropy axis parallel to the z‐axis is considered. The purpose of the paper is to determine displacements and stresses fields in the interior of the half‐space when a rigid circular disk of radius a completely bonded to the surface of the half‐space is rotated through a constant angle θ₀. The region of the surface lying out with the circle r?a, is free from stresses. This problem is a type of Reissner–Sagoci mixed boundary value problems. Using cylindrical co‐ordinate system and applying Hankel integral transform in the radial direction, the problem may be changed to a system of dual integral equations. The solution of the dual integral equations is obtained by an approach analogous to Sneddon's (J. Appl. Phys. 1947; 18 :130–132), so that the circumferential displacement and stress fields inside the medium are obtained analytically. The same problem has already been approached by Hanson and Puja (J. Appl. Mech. 1997; 64 :692–694) by the use of integrating the point force potential functions. It is analytically proved that the present solution, although of a quite different form, is equivalent to that given by Hanson and Puja. To illustrate the solution, a few plots are provided. The displacements and the stresses in a soil deposit due to a rotationally symmetric force or boundary displacement may be obtained using the results of this paper. Copyright © 2006 John Wiley & Sons, Ltd.     

An analytical solution for stresses induced by a post-tensioned anchor in rocks containing two perpendicular joint sets

A new analytical approach for the analysis of stress around a post-tensioned anchor in rock with two perpendicular joint sets is presented. The solution considers nonlinear shear stresses developed along the anchor bond length as well as debonding at the tendon–grout interface. The bearing plate effects are also considered in the analyses. The following assumptions are made: (1) homogeneous and orthotropic elastic rock; (2) half-space rock mass with plane strain conditions; (3) trilinear bond-slip model for the behavior of the tendon–grout interface; (4) an elastic anchor. The employed methodology is to decompose the anchor problem into two problems of simpler loadings: stresses produced under a rigid bearing plate and stresses induced by interfacial shear stresses mobilized along the bond length. Based on the proposed solution, an illustrative example is given and the results show that a large compressive region is formed around the anchor free length. Tensile stress concentrations are also observed around the bond length and near the rock surface outside the bearing plate. A parametric study is also conducted to investigate the effects of joint properties on the induced stresses. The results show that as the number of joints intersecting the anchor axis increases, the magnitude of compressive stresses in the free zone decreases, while the size of compression zone around the anchor increases. To verify the results of the analytical approach, a comparison is made with numerical results obtained by using the finite element method. The analytical solutions compare very well with the results obtained by numerical method.     

Rocking vibration of a rigid strip footing on saturated soil

The dynamic response of a rigid strip footing lying on saturated soil is greatly affected by the pore pressure induced by a rocking moment. To consider the complex behavior of the soil under the rocking load, an analytical solution for a rigid strip foundation on saturated soil under a rocking moment is developed under the framework of Biot’s coupling theory. The boundary-value problem for the governing coupling equations for saturated soil is solved using a Fourier transform to yield a pair of dual integral equations. These dual integral equations are transformed into a set of linear equations using an infinite series of orthogonal Jacobi polynomials to yield the compliance functions. In addition, a parametric study has been carried out to examine the influence of: (1) the dimensionless frequency, (2) the dynamic permeability and (3) the Poisson’s ratio on saturated soil under a rocking rigid strip footing.     

Seismic uplift capacity of inclined strip anchors in sand using upper bound limit analysis

Pseudo-static approach is adopted in this paper to determine the seismic uplift capacity of an inclined strip anchor using upper bound limit analysis. Two different failure mechanisms are considered to obtain the magnitudes of unit weight component of uplift factor f_γE for different values of soil friction angle, interface friction of anchor plate, anchor inclination, embedment ratio and horizontal seismic acceleration coefficient. The failure mechanism 1 consists of a triangular and quadrilateral rigid blocks; whereas the failure mechanism 2 comprises a logarithmic spiral failure zone with varied focus, sandwiched between a triangular and quadrilateral rigid blocks. It is observed that the magnitude of uplift factor f_γE decreases significantly with the increase in seismic acceleration but increases with the increase in embedment ratio and roughness of the anchor surface. However, a mixed trend in the values of f_γE can be observed for different inclination of the anchor, which is clearly discussed in this paper. The results are compared with the existing values in the literature and the significance of the present methodology for designing the inclined strip anchor is discussed.     

A note on the consolidation settlement of a rigid circular foundation on a poroelastic halfspace

This paper uses Biot's poroelasticity approach to examine the consolidation behaviour of a rigid foundation with a frictionless base in contact with a poroelastic halfspace. The mathematical development of the mixed boundary value problem involves a set of dual integral equations in the Laplace transform domain which cannot be conveniently solved by employing conventional procedures. In this paper, a numerical solution is developed using a scheme where the contact normal stress is approximated by a discretized equivalent. The influence of limiting drainage boundary conditions at the entire surface of the halfspace on the degree of consolidation of the rigid circular foundation is investigated. The results obtained in this study are compared with the corresponding results given in the literature. Copyright © 2016 John Wiley & Sons, Ltd.     

隧道挤压型变形与支护特性研究

针对挤压型大变形问题,通过木寨岭隧道科研试验与相关资料分析,得出挤压型大变形具有变形量大、持续时间长的特点,一般两侧为变形优势部位。测试数据表明,刚性支护中围岩压力、拱架应力、锚杆轴力值均较大,其中拱架应力值可能超过材料屈服强度,边墙锚杆较拱部锚杆作用明显,围岩压力与变形增量呈反比关系;在挤压型变形隧道中允许有控制性地释放围岩变形,能够减小支护结构上的围岩压力,有利于隧道结构的长期稳定。文中阐述了可让式支护原理,对可让式支护结构预留变形量进行初步预测,数值计算表明,采用刚性支护将位移控制在较小的范围内需要足够大的刚度,而可让式支护具有能够允许围岩一定程度的变形同时减小支护结构受力的优势,可在铁路隧道中推广使用。     

基于3D-FSM的预紧力端锚数值模拟及其多核并行化研究

在对传统锚杆荷载传递机制分析的基础上,提出了一种考虑托板、锚杆与岩体相互作用的数值模拟方法：通过Kelvin基本解计算出锚杆集中力对围岩的影响,同时就锚杆自由段与相应岩体两端点的位移差相等建立位移方程,结合先前开发的3D-FSM数值模拟系统中的表面受力平衡方程进行联立求解,利用所得结果可以计算域内任意点的应力及位移变化,形成完整的预紧力端锚边界元数值模拟系统。通过与Flac3D系统模拟对比,验证了该系统的可靠性。为提高运算效率,对该系统进行基于OpenMP的多核并行化改进,给出了改进的基本思路和加速比对比图。由于边界元本身具有建模简单、计算区域大、计算精度高等优点,因此,这种模拟方法有很大的应用价值。     

Rocking vibration of a rigid circular disc in a transversely isotropic full‐space

In this paper, forced rocking vibration of a rigid circular disc placed in a transversely isotropic full‐space, where the axis of material symmetry of the full‐space is normal to the surface of the plate, is analytically investigated. Because of using the Fourier series and Hankel integral transforms, the mixed boundary‐value problem is transformed into two separate pairs of integral equations called dual integral equations. The dual integral equations involved in this paper are reduced to Fredholm integral equations of the second kind. With the aid of contour integration, the governing integral equation is numerically evaluated in the general dynamic case. The reduced static case of the dual integral equations is solved analytically and the vertical displacement, the contact pressure and the static impedance/compliance function are explicitly determined, and it is shown that the pressure in between the plate and the full‐space and the compliance function reduced for isotropic half‐space are identical to the previously published solutions. The dynamic contact pressure in between the disc and the space and also the related impedance function are numerically evaluated in general dynamic case and illustrated. It is shown that the singularity exists in the contact pressure at the edge of the disc is the same as the static case. To show the effect of material anisotropy, the numerical evaluations are given for some different transversely isotropic materials and compared. Copyright © 2010 John Wiley & Sons, Ltd.     

广义Gibson饱和地基中刚性圆板的扭转柔度系数

考虑到地基在竖直方向上的非均匀性,结合扭转振动的特点,建立了简谐扭转动荷载作用时剪切模量随深度增大的广义Gibson饱和地基的动力方程,通过积分变换求解了动力方程。考虑到半空间地基表面处应力自由、埋置圆板所在平面为混合边界和无穷深度处为波的辐射边界等条件,得到了广义Gibson饱和地基中刚性圆板扭转振动时的对偶积分方程,通过合适的变换转化成了一个第2类Fredholm积分方程,求解了相应的动力响应问题。对比静扭距作用时的荷载-位移关系,给出了动力柔度系数和扭转角位移幅值的表达式,并把所研究的问题进行退化且与前人成果进行了对比。数值研究表明：当基础的埋置深度小于5倍基础半径时,广义Gibson饱和地基中埋置基础的扭转振动存在明显的边界层现象,且埋置深度越小,边界层现象越明显。