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.     

An integral equation solution for three‐dimensional heat extraction from planar fracture in hot dry rock

In the numerical simulation of heat extraction by circulating water in a fracture embedded in geothermal reservoir, the heat conduction in the reservoir is typically assumed to be one‐dimensional and perpendicular to the fracture in order to avoid the discretization of the three‐dimensional reservoir geometry. In this paper we demonstrate that by utilizing the integral equation formulation with a Green's function, the three‐dimensional heat flow in the reservoir can be modelled without the need of discretizing the reservoir. Numerical results show that the three‐dimensional heat conduction effect can significantly alter the prediction of heat extraction temperature and the reservoir life as compared to its one‐dimensional simplification. Copyright © 2003 John Wiley & Sons, Ltd.     

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.     

Integral equation solution of heat extraction‐induced thermal stress in enhanced geothermal reservoirs

During fluid injection in enhanced geothermal systems, thermo‐mechanical processes can play an important role. In fact, the phenomena of reservoir seismicity and the variation of injectivity with respect to injection water temperature can be attributed to the induced thermal stresses. In this paper, a three‐dimensional integral equation formulation is presented for calculating thermally induced stresses associated with the cooling of a fracture in a geothermal reservoir. By utilizing Green's function in the integral equation, the three‐dimensional heat flow and stresses in the reservoir are modelled without discretizing the reservoir. The formulation is implemented in a computer program for the solution of injection into an infinite fracture as well as for the injection/extraction in an arbitrarily shaped fracture. Copyright © 2005 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.     

A quasistatic analysis of a plate on consolidating layered soils by analytical layer‐element/finite element method coupling

In this paper, a coupling method between finite element and analytical layer‐elements is utilized to analyze the time‐dependent behavior of a plate of any shape and finite rigidity resting on layered saturated soils. Based on the integral transform techniques together with the aid of an order reduction method, an analytical layer‐element solution is derived from the governing equations for three‐dimensional Biot consolidation with respect to a Cartesian coordinate system and then extended to be the fundamental solution for the layered saturated soil under a point load. The Mindlin plate is modeled by eight‐noded isoparametric elements. The governing equations of the interaction between soil and plate in the Laplace‐Fourier transformed domain are deduced by referring to the coupling theory of FEM/BEM, and the final solution is obtained by applying numerical inversion. Numerical examples concerned with the time‐dependent response of a plate are performed to demonstrate the influence of soil and plate properties on the interaction process. Copyright © 2014 John Wiley & Sons, Ltd.     

Simplified analytical solution for point load acting behind a cantilever wall

The analysis of earth pressure resulting from an applied point load behind a retaining wall is to be carried out in three dimensions as the earth pressure distribution varies both in vertical and horizontal directions. In this research, the semi‐empirical expressions developed by Terzaghi [Trans. ASCE 1954; 119 ] for earth pressure due to point load are integrated and a correction factor is introduced to develop equivalent equations for two‐dimensional analyses. A new parameter referred to as the design width, which is the horizontal distance retained by the individual vertical retaining element is introduced. The used procedure and the resulted equations are tested and verified by adopting different design width values covering the practical range. The resulted equations together with the semi‐empirical expressions of earth pressure due to line load developed by Terzaghi [Trans. ASCE 1954; 119 ] can be easily used for two‐dimensional analysis of the effect of point load. Copyright © 2011 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.     

Frictional contact of laminated elastic half‐spaces allowing interface cavities. Part 2: Numerical results

The paper presents some numerical results of the problem on contact of stratified elastic half‐spaces possessing unevenness of their surface layers. The approach of the Jacobi polynomials is used to reduce the singular integral equation of the problem to a system of linear algebraic equations. On the basis of the system's solution the dependence of gaps' geometric characteristics on external load and friction is investigated for two types of initial boundary disturbances, namely a local symmetric recess and a periodic set of recesses at a flat boundary. Numerical results are presented in figures for a few elastic and geometric parameters of joint components. The analysis shows that elastic properties and geometric structure of a laminated half‐space have a significant influence on the shape of interface gaps and their location. Copyright © 2001 John Wiley & Sons, Ltd.     

矩形运动荷载作用下软土地基的三维振动研究

通过引入势函数,并利用Helmholtz原理和Fourier变换技术,研究了运动荷载作用下有限层厚软土地基的振动,考虑了矩形分布荷载作用下振动的三维特性,使得分析更符合工程实际,给出了运动荷载作用下饱和黏弹性地基三维振动的应力、位移和孔隙水压力响应的积分形式解答。利用Fourier数值逆变换进行算例分析,讨论了荷载速度对位移及孔隙水压力分布的影响。结果表明,位移幅值随荷载速度的增加而增大,荷载不同位置处孔隙水压力的分布有很大差异。     

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.     

Integral equation solution of heat extraction from a fracture in hot dry rock

In the study of heat extraction by circulating water in a fracture embedded in geothermal reservoir, the heat conduction in the reservoir is typically assumed to be one dimensional and perpendicular to the fracture. In this paper, we demonstrate that by an integral equation formulation utilizing Green's function, the multi‐dimensional heat flow in the reservoir can be modelled. In the resulting numerical solution system, the discretization of reservoir geometry is entirely eliminated, leading to a much more efficient scheme. The multi‐dimensional heat conduction effect as compared to its one‐dimensional simplification is studied. Copyright © 2001 John Wiley & Sons, Ltd.     

A simulation using a hybrid method for predicting fault zones ahead of a tunnel face

By monitoring and analysing three‐dimensional absolute displacement during tunnel excavation, it is possible to predict discontinuity zones ahead of tunnel face. This paper presents results of the three‐dimensional displacement obtained from three‐dimensional simulations using hybrid method for assessing the effect of discontinuity zones on monitoring data. The influence of fault zones and in situ initial stresses on the three‐dimensional tunnel movement is evaluated. Copyright © 2002 John Wiley & Sons, Ltd.     

Analysis of piles subjected to lateral soil movements using a three‐dimensional displacement approach

An analytical approach using the three‐dimensional displacement of a soil is investigated to provide analytical solutions of the horizontal response of a circular pile subjected to lateral soil movements in nonhomogeneous soil. The lateral stiffness coefficient of the pile shaft in nonhomogeneous soil is derived from the rocking stiffness coefficient that is obtained from the analytical solution, taking into account the three‐dimensional displacement represented in terms of scalar potentials in the elastic three‐dimensional analysis. The relationship between horizontal displacement, rotation, moment, and shear force of a pile subjected to lateral soil movements in nonhomogeneous soil is obtainable in the form of the recurrence equation. For the relationship between the lateral pressure and the horizontal displacement, it is assumed that the behavior is linear elastic up to lateral soil yield, and the lateral pressure is constant under the lateral soil yield. The interaction factors between piles subjected to both lateral load and moment are calculated, taking into account the lateral soil movement. The formulation of the lateral displacement and rotation of the pile base subjected to lateral loads in nonhomogeneous soils is presented by taking into account the Mindlin equation and the equivalent thickness for soil layers in the equivalent elastic method. For lateral movement, lateral pressure, bending moment, and interaction factors, there are small differences between results obtained from the 1‐D and the 3‐D displacement methods except a very flexible pile. Copyright © 2015 John Wiley & Sons, Ltd.     

Three‐dimensional finite element study of arching behavior in slope/drilled shafts system

Two‐dimensional slope stability analysis for a slope with a row of drilled shafts needs a mechanism to take into account the three‐dimensional effect of the soil arching due to the spaced drilled shafts on slope. To gain a better understanding of the arching mechanisms in a slope with evenly spaced drilled shafts socketed into a stable stratum (or a rock layer), the three‐dimensional finite element modelling technique was used for a comprehensive parametric study, where the nonlinear and plastic nature of the soil and the elastic behavior of the drilled shafts as well as the interface frictions were modelled. Various factors were varied in the parametric study to include (1) the rigidity of the drilled shafts as influenced by its diameter, modulus of elasticity, and total length; (2) shafts spacing and location on the slope; (3) the material properties of rock and the socket length of shaft; and (4) the soil movement and strength parameters. Evidences of soil arching and reduction in the stresses and displacements through the load transfer mechanisms due to the presence of the drilled shafts were elucidated through the finite element method (FEM) simulation results. Design charts based on regression analysis of FEM simulation results were created to obtain a numerical value of the load transfer factor for the arching mechanism provided by the drilled shafts on the slope. Observations of the arching behavior learned from the FEM simulations provide an insight into the behavior of drilled shafts stabilized slope. Copyright © 2009 John Wiley & Sons, Ltd.     

Physics‐based simulation of multiple interacting crack growth in brittle rocks driven by thermal cooling

Crack growth in hot brittle rocks, driven by thermal cooling, was simulated using a coupled two‐dimensional discrete element and heat transport model that explicitly includes the random initiation and subsequent propagation of interacting cracks. The model clearly predicts that a quasi‐hierarchical array of subparallel cracks, oriented along the direction of the temperature gradient, is formed under small to moderately large thermally generated strain load conditions. The simulation results also demonstrate that, after an initial transient, thermal cracks propagate in a stable fashion with a velocity that scales with ~ t^{− 1/2}. However, under large thermal strain loads, a more complicated geometry composed of cracks that curve and coalesce develops during the later stages of crack growth. Copyright © 2016 John Wiley & Sons, Ltd.     

Response evaluation for horizontally loaded fixed‐head pile groups using 3‐D non‐linear analysis

The response of laterally loaded pile foundations may be significantly important in the design of structures for such loads. A static horizontal pile load test is able to provide a load–deflection curve for a single free‐head pile, which significantly differs from that of a free‐ or fixed‐head pile group, depending on the particular group configuration. The aim of this paper is to evaluate the influence of the interaction between the piles of a group fixed in a rigid pile cap on both the lateral load capacity and the stiffness of the group. For this purpose, a parametric three‐dimensional non‐linear numerical analysis was carried out for different arrangements of pile groups. The response of the pile groups is compared to that of the single pile. The influence of the number of piles, the spacing and the deflection level to the group response is discussed. Furthermore, the contribution of the piles constituting the group to the total group resistance is examined. Finally, a relationship is proposed allowing a reasonable prediction of the response of fixed‐head pile groups at least for similar soil profile conditions. Copyright © 2005 John Wiley & Sons, Ltd.     

Glyph and hyperstreamline representation of stress and strain tensors and material constitutive response

Results of numerical analyses of boundary value problems in geomechanics include output of three‐dimensional stress and strain states. Two‐dimensional plots of stress–stress or stress–strain quantities, often used to represent such output, do not fully communicate the evolution of stress and strain states. This paper describes the use of glyphs and hyperstreamlines for the visual representation of three dimensional stress and strain tensors in geomechanics applications. Glyphs can be used to represent principal stress states as well as normal stresses at a point. The application of these glyphs is extended in this paper to represent strain states. The paper introduces a new glyph, called HWY glyph for the representation of shear tensor components. A load step‐based hyperstreamline is developed to show the evolution of a stress or strain tensor under a general state of loading. The evolution of stress–strain states from simulated laboratory tests and a general boundary value problem of a deep braced excavation are represented using these advanced visual techniques. These visual representations facilitate the understanding of complex multidimensional stress–strain soil constitutive relationships. The visual objects introduced in this paper can be applied to stress and strain tensors from general boundary value problems. Copyright © 2003 John Wiley & Sons, Ltd.     

Tensor visualizations in computational geomechanics

We present a novel technique for visualizing tensors in three dimensional (3D) space. Of particular interest is the visualization of stress tensors resulting from 3D numerical simulations in computational geomechanics. To this end we present three different approaches to visualizing tensors in 3D space, namely hedgehogs, hyperstreamlines and hyperstreamsurfaces. We also present a number of examples related to stress distributions in 3D solids subjected to single and load couples. In addition, we present stress visualizations resulting from single‐pile and pile‐group computations. The main objective of this work is to investigate various techniques for visualizing general Cartesian tensors of rank 2 and it's application to geomechanics problems. Copyright © 2002 John Wiley & Sons, Ltd.     

Dynamic displacement and stress solutions for viscoelastic half‐spaces subjected to harmonic concentrated loads using the Radon and Fourier transforms

In this article a numerical solution for a three‐dimensional isotropic, viscoelastic half‐space subjected to concentrated surface stress loadings is synthesized with the aid of the Radon and Fourier integral transforms. Dynamic displacement and stress fields are computed for points at the surface and inside the domain. The analysis is performed in the frequency domain. Viscoelastic effects are incorporated by means of the elastic–viscoelastic correspondence principle. The equations of motion are solved in the Radon–Fourier transformed domain. Inverse transformations to the physical domain are accomplished numerically. The scheme used to perform the numerical inverse transformations is addressed. The solution is validated by comparison with results available in the literature. A set of original dynamic displacement and stress solutions for points within the half‐space is presented. Copyright © 2009 John Wiley & Sons, Ltd.