Elastostatic response of a pile embedded in a transversely isotropic half‐space under transverse loading

In the framework of elastostatics, a mathematical treatment is presented for the boundary value problem of the interaction of a flexible cylindrical pile embedded in a transversely isotropic half‐space under transverse loadings. Taking the pile region as a stiffened subdomain of an extended half‐space, the formulation of the interaction problem is reduced to a Fredholm integral equation of the second kind. The necessary set of Green's functions for the transversely isotropic half‐space is obtained by means of a method of potentials. The resulting Green's functions are incorporated into a numerical procedure for the solution of the integral equation. The theoretical response of the pile is presented in terms of bending moment, displacement and slope profiles for a variety of transversely isotropic materials so that the effect of different anisotropy parameters can be meaningfully discussed. Copyright © 2013 John Wiley & Sons, Ltd.     

Axisymmetric interaction of a rigid disc with a transversely isotropic half‐space

A theoretical formulation is presented for the determination of the interaction of a vertically loaded disc embedded in a transversely isotropic half‐space. By means of a complete representation using a displacement potential, it is shown that the governing equations of motion for this class of problems can be uncoupled into a fourth‐order partial differential equation. With the aid of Hankel transforms, a relaxed treatment of the mixed‐boundary value problem is formulated as dual integral equations, which, in turn, are reduced to a Fredholm equation of the second kind. In addition to furnishing a unified view of existing solutions for zero and infinite embedments, the present treatment reveals a severe boundary‐layer phenomenon, which is apt to be of interest to this class of problems in general. The present solutions are analytically in exact agreement with the existing solutions for a half‐space with isotropic material properties. To confirm the accuracy of the numerical evaluation of the integrals involved, numerical results are included for cases of different degrees of the material anisotropy and compared with existing solutions. Further numerical examples are also presented to elucidate the influence of the degree of the material anisotropy on the response. Copyright © 2009 John Wiley & Sons, Ltd.     

Interaction of a cylindrical thin‐walled pile with an inhomogeneous isotropic elastic medium

This paper presents a rigorous analysis for the static interaction of a cylindrical thin‐walled pile with an inhomogeneous isotropic elastic half‐space under vertical, horizontal, and torsional forces individually applied at the top of pile. The inhomogeneity is specified with the exponential variation of shear modulus along depth of the embedding medium, and the Poisson's ratio is assumed to be constant. By means of a set of Green's functions for pile and soil medium and satisfying the compatibility conditions between the 2 interacting media, the formulation is reduced to coupled Fredholm integral equations. Using the adaptive‐gradient elements, capable of capturing the singular stress transfer at both ends of the pile, a numerical procedure is developed and utilized for evaluating the relevant integral equations and studying the inhomogeneity effect on the soil‐pile interaction responses. The analysis results have been validated for different soil‐pile modulus ratios under axial load and for a Poisson's ratio of 0.3 under lateral load. The procedure does not consider the nonlinear behavior of the soil medium or plastic yielding in the pile section, and the impact of the unreliable results for the case of high Poisson's ratio is not examined.     

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.     

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.     

Load–Deflection response of transversely isotropic piles under lateral loads

In general, pile materials are assumed to be isotropic during the analysis of the load–deflection response of piles under lateral loads. However, commonly used materials such as reinforced concrete and timber as well as potentially promising new pile materials such as fiber reinforced polymers are typically transversely isotropic materials. Experimental studies have shown that transversely isotropic materials have a high ratio of section longitudinal modulus to the section in‐plane shear modulus (E_zz/G_xz) compared to the value for isotropic materials. The high modulus ratio leads to a more significant shear deformation effect in beam bending. To account for the shear deformation effect, the Timoshenko Beam Theory has been adopted in deriving the solutions for the load–deflection response of transversely isotropic piles under lateral loads instead of the Classical (Euler–Bernoulli) Beam Theory. The load–deflection responses depend on the shear effect coefficient, the lateral soil resistance, the embedment ratio, and the boundary conditions. The deflection of the pile, if the shear deformation effect is considered, is always larger than if it is neglected. Copyright © 2000 John Wiley & Sons, Ltd.     

Displacement ring load Green's functions for saturated porous transversely isotropic tri‐material full‐space

Based on the Biot's poroelastic theory and using scalar potential functions both the ring load and point load displacement Green's functions for a transversely isotropic saturated porous full‐space composed of an upper half‐space, a finite thickness middle layer and a lower half‐space is analytically presented for the first time. It is assumed that each region consists of a different transversely isotropic material. The equations of poroelastodymanics in terms of the solid displacements and the pore fluid pressure are uncoupled with the help of two scalar potential functions, so that the governing equations for the potential functions are either a second order wave equation or a repeated wave‐heat transfer equation of sixth order. With the aid of Fourier expansion with respect to circumferential direction and Hankel integral transforms with respect to the radial direction in cylindrical coordinate system, the response is determined in the form of line integrals in the real space, followed by theorem of inverse Hankel integral transforms. The solutions degenerate to a single phase elastic material, and the results are compared with previous studies, where an excellent agreement may be observed with the results provided in the literature. Some examples of displacement Green's functions are finally given to illustrate the solution. Copyright © 2016 John Wiley & Sons, Ltd.     

Elastic solutions for stresses in a transversely isotropic half‐space subjected to three‐dimensional buried parabolic rectangular loads

In many areas of engineering practice, applied loads are not uniformly distributed but often concentrated towards the centre of a foundation. Thus, loads are more realistically depicted as distributed as linearly varying or as parabola of revolution. Solutions for stresses in a transversely isotropic half‐space caused by concave and convex parabolic loads that act on a rectangle have not been derived. This work proposes analytical solutions for stresses in a transversely isotropic half‐space, induced by three‐dimensional, buried, linearly varying/uniform/parabolic rectangular loads. Load types include an upwardly and a downwardly linearly varying load, a uniform load, a concave and a convex parabolic load, all distributed over a rectangular area. These solutions are obtained by integrating the point load solutions in a Cartesian co‐ordinate system for a transversely isotropic half‐space. The buried depth, the dimensions of the loaded area, the type and degree of material anisotropy and the loading type for transversely isotropic half‐spaces influence the proposed solutions. An illustrative example is presented to elucidate the effect of the dimensions of the loaded area, the type and degree of rock anisotropy, and the type of loading on the vertical stress in the isotropic/transversely isotropic rocks subjected to a linearly varying/uniform/parabolic rectangular load. Copyright © 2002 John Wiley & Sons, Ltd.     

Dynamic responses of a pile embedded in a layered poroelastic half‐space to harmonic lateral loads

A single pile embedded in a layered poroelastic half‐space subjected to a harmonic lateral load is investigated in this study. Based on Biot's theory, the frequency domain fundamental solution for a horizontal circular patch load applied in the layered poroelastic half‐space is derived via the transmission and reflection matrices method. Utilizing Muki and Sternberg's method, the second kind of Fredholm integral equation describing the dynamic interaction between the layered half‐space and the pile subjected to a top harmonic lateral load is constructed. The proposed methodology is validated by comparing results of this paper with some existing results. Numerical results show that for a two‐layered half‐space, the thickness of the upper softer layer has pronounced influences on the dynamic response of the pile and the half‐space. For a three‐layered half‐space, the presence of a softer middle layer in the layered half‐space will enhance the compliance for the pile significantly, while a stiffer middle layer will diminish the dynamic compliance of the pile considerably. Copyright © 2009 John Wiley & Sons, Ltd.     

Tensionless–frictionless interaction of flexible annular foundation with a transversely isotropic multi‐layered half‐space

A transversely isotropic multi‐layered half‐space, with axis of material symmetry perpendicular to the free surface, supports a flexible either annular or solid circle foundation. The contact area of the foundation and the half‐space is considered to be both frictionless and tensionless. The foundation is assumed to be affected by a vertical static axisymmetric load. Detailed analysis of the interaction of these two systems with different thickness of layers is the target of this paper. With the use of ring load Green's functions for both the foundation and the continuum half‐space, an integral equation accompanied with some inequalities is introduced to model the complex BVP. With the incorporation of ring‐shape FEM, we are capable of capturing both regular and singular solution smoothly. The validity of the combination of the analytical and numerical method is proved with comparing the results of this paper with a number of benchmark cases of both linear and nonlinear interaction of circular and annular foundation with half‐space. Some new illustrations are presented to portray the aspect of the anisotropy and layering of the half‐space. Copyright © 2014 John Wiley & Sons, Ltd.     

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.     

Poroelastic model for pile–soil interaction in a half‐space porous medium due to seismic waves

In this paper, frequency domain dynamic response of a pile embedded in a half‐space porous medium and subjected to P, SV seismic waves is investigated. According to the fictitious pile methodology, the problem is decomposed into an extended poroelastic half‐space and a fictitious pile. The extended porous half‐space is described by Biot's theory, while the fictitious pile is treated as a bar and a beam and described by the conventional 1‐D structure vibration theory. Using the Hankel transformation method, the fundamental solutions for a half‐space porous medium subjected to a vertical or a horizontal circular patch load are established. Based on the obtained fundamental solutions and free wave fields, the second kind of Fredholm integral equations describing the vertical and the horizontal interaction between the pile and the poroelastic half‐space are established. Solution of the integral equations yields the dynamic response of the pile to plane P, SV waves. Numerical results show the parameters of the porous medium, the pile and incident waves have direct influences on the dynamic response of the pile–half‐space system. Significant differences between conventional single‐phase elastic model and the poroelastic model for the surrounding medium of the pile are found. Copyright © 2007 John Wiley & Sons, Ltd.     

横观各向同性饱和地基竖向振动的衰减特性

基于提出的横观各向同性饱和多孔介质Biot波动方程的一般解,研究了饱和半空间地基在竖向点源简谐激振荷载作用下地表振动的衰减特征,分析了激振频率以及横观各向同性饱和土介质的各向异性参数和孔隙渗透系数对地表振动特征的影响。计算结果表明,低频和高频激振时,地表位移衰减特性存在明显差异;在饱和土的各向异性参数中,纵向和水平方向动态渗透系数比值和刚度系数比值对地表位移衰减影响最大,这也说明采用各向同性饱和介质的动力学模型不能准确地描述具有明显各向异性特性的饱和土地基的动力特性。     

交通荷载下横观各向同性土的Biot固结分析

在竖向荷载和切向荷载共同作用下,求解了具有下卧基岩的横观各向同性土的Biot固结问题。首先,基于Biot固结理论得到控制方程,同时引入3个状态变量;然后,对时间t进行Laplace变换和对半径r进行Hankel变换。利用Laplace-Hankel联合积分变换求解状态方程,得到交通荷载作用下横观各向同性土骨架位移、孔隙水压力等的一般积分形式解。结合算例验证了该方法的正确性,并分析了切向荷载和竖向荷载共同作用下对横观各向同性土层沉降和孔压的影响,可以发现切向荷载和竖向荷载共同作用所产生的前期沉降明显小于竖向荷载单独作用的情况。为了更好地模拟实际情况,竖向荷载还采用正弦循环荷载。在循环荷载作用下,沉降变化周期滞后于荷载周期。     

Dynamic response to fluid extraction from a poroelastic half‐space

In this study, the dynamic response of a poroelastic half‐space to a point fluid sink is investigated using Biot's dynamic theory of poroelasticity. Based on Biot's theory, the governing field equations are re‐formulated in frequency domain with solid displacement and pore pressure. In a cylindrical coordinate system, a method of displacement potentials for axisymmetric displacement field is proposed to decouple the Biot's field equations to three scalar Helmholtz equations, and then the general solution to axisymmetric problems are obtained. The full‐space fundamental singular solution for a point sink is also derived using potential methods. The mirror‐image method is finally applied to construct the fundamental solution for a point sink buried in a poroelastic half‐space. Furthermore, a numerical study is conducted for a rock, that is, Berea sandstone, as a representative example. Copyright © 2013 John Wiley & Sons, Ltd.     

Rocking vibrations of a rigid embedded foundation in a poroelastic half‐space

This paper presents an analysis of the rocking vibrations of a rigid cylindrical foundation embedded in poroelastic soil. The foundation is subjected to time‐harmonic rocking excitation and is perfectly bonded to the surrounding soil. The soil underlying the foundation base is represented by a homogeneous poroelastic half‐space, whereas the soil along the side of the foundation is modeled as an independent poroelastic stratum composed of a series of infinitesimally thin layers. The behavior of the soil is governed by Biot's poroelastodynamic theory. The contact surface between the foundation base and the poroelastic soil is assumed to be smooth and either fully permeable or impermeable. The dynamic interaction problem is solved by employing a simplified analytical method. Some numerical results for the nondimensional rocking dynamic impedance and nondimensional angular displacement amplitude of the foundation are presented to show the effect of nondimensional frequency of excitation, poroelastic material parameters, hydraulic boundary condition, depth ratio and mass ratio of the foundation. Copyright © 2009 John Wiley & Sons, Ltd.     

Time‐dependent behaviour of a rigid foundation on a transversely isotropic soil layer

An analytical solution is presented in this paper to study the time‐dependent settlement behaviour of a rigid foundation resting on a transversely isotropic saturated soil layer. The governing equations for a transversely isotropic saturated soil, within Biot's poroelasticity framework, are solved by means of Laplace and Hankel transforms. The problem is subsequently formulated in the Laplace transform domain in terms of a set of dual integral equations that are further reduced to a Fredholm integral equation of the second kind and solved numerically. The developed analytical solution is validated via comparison with the existing analytical solution for an isotropic saturated soil case, and adopted as a benchmark to examine the sensitivities of the mesh refinement and the locations of truncation boundaries in the finite element simulations using ABAQUS. Particular attention is paid to the influences of the degree of soil anisotropy, boundary drainage condition, and the soil layer thickness on the consolidation settlement and contact stress of the rigid foundation. Copyright © 2009 John Wiley & Sons, Ltd.     

横观各向同性地基中管桩扭转振动响应解析解

考虑地基沉积过程中产生的竖向和水平向力学性质的差异,对横观各向同性地基中管桩扭转振动频域响应进行了理论研究。基于横观各向同性材料的本构关系以及桩-土耦合扭转振动,建立了桩土系统定解问题,通过Laplace变换和分离变量法求得了桩周土和桩芯土扭转振动位移形式解。通过桩-土接触面的连续条件,求得了管桩扭转频域响应解析解,并得到了桩顶复动刚度和速度导纳的表达式。将所得解退化到横观各向同性地基中实心桩解以及均匀地基中管桩解,并与已有文献进行了对比,验证了解的合理性。通过数值算例,分析了桩周土和桩芯土的横观各向同性力学参数对桩顶扭转复刚度及速度导纳的影响。     

The axisymmetric consolidation of a semi‐infinite transversely isotropic saturated soil

This paper presents an exact analytical solution to fully coupled axisymmetric consolidation of a semi‐infinite, transversely isotropic saturated soil subjected to a uniform circular loading at the ground surface. The analysis is under the framework of Biot's general theory of consolidation. First, the governing equations of consolidation are transformed into a set of equivalent partial differential equations with the introduction of two auxiliary variables. These partial differential equations are then solved using Hankel–Laplace integral transforms. Once solutions in the transformed domain have been obtained, the actual solutions in the physical domain for displacements and stress components of the solid matrix, pore‐water pressure and fluid discharge can be finally obtained by direct numerical inversion. The accuracy of the numerical solutions developed is confirmed by comparison with an existing exact solution for an isotropic and saturated soil that is a special case of the more general problem addressed. Numerical analyses are also presented to investigate the influence of the degree of material anisotropy on the consolidation settlement. Copyright © 2005 John Wiley & Sons, Ltd.     

Impedances of rigid cylindrical foundations embedded in transversely isotropic soils

A complete formulation and implementation for assessment of the response to dynamic loads of cylindrical rigid structures embedded in transversely isotropic elastic half‐spaces is presented. The analysis is performed in the frequency domain and the steady‐state structure response is obtained. The method is based on a non‐singular version of the indirect boundary element method which uses influence functions, instead of Green's functions, as fundamental solutions. These influence functions are the response of an elastic half‐space to distributed, internally applied loads. The proposed method imposes full bonding contact between the foundation and the surrounding soil. Numerical results for displacement (vertical and horizontal) and rotation (twisting and rocking) impedances, showing the influence of the soil anisotropy, are presented. Results for the soil–structure interface tractions and for the displacement field throughout the half‐space are also shown. Copyright © 2006 John Wiley & Sons, Ltd.