The indentation of a half-space of hexagonal elastic material by a circular punch of arbitrary end-profile

Summary The problem is considered of the indentation by a smooth rigid punch of a half-space composed of linear elastic material of hexagonal symmetry whose plane boundary is parallel to the basal planes. The case is considered in which the area of contact between the punch and the half-space is circular, the end of the punch with is in contact with the half-space having an arbitrary profile. An integral equation is formulated and solved for the boundary value of the normal displacement in the half-space, and an expression is derived for the distribution of pressure under the punch.     

Sliding frictional contact analysis of functionally graded piezoelectric layered half-plane

This paper investigates the sliding frictional contact problem of a layered half-plane made of functionally graded piezoelectric materials (FGPMs) in the plane strain state. It is assumed that the punch is a perfect electrical insulator with zero electric charge distribution, and the friction within the contact region is of Coulomb type. The electro-elastic properties of the FGPM layer vary exponentially along the thickness direction. The fundamental solutions for the applied concentrated linear forces perpendicular and parallel to the FGPM layer surface are obtained. Using the superposition theorem, the problem is reduced to a Cauchy singular integral equation which is then numerically solved to determine the contact tractions, contact region, maximum indentation depth, electrical potential and electromechanical fields. Numerical results show that both the material property gradient and the friction coefficient have significant influence on the contact performance of the FGPM layered half-plane.     

Boundary integral equation method to solve embedded planar crack problems under shear loading

The solution of three-dimensional planar cracks under shear loading are investigated by the boundary integral equation method. A system of two hypersingular integral equations of a three-dimensional elastic solid with an embedded planar crack are given. The solution of the boundary integral equations is succeeded taking into consideration an appropriate Gauss quadrature rule for finite part integrals which is suitable for the numerical treatment of any plane crack without a polygonal contour shape and permit the fast convergence for the results. The stress intensity factors at the crack front are calculated in the case of a circular and an elliptic crack and are compared with the analytical solution.     

Edge crack in an elastic layer resting on Winkler foundation

The paper deals with the stress analysis near a crack tip in an elastic layer resting on Winkler foundation. The edge crack is assumed to be normal to the lower boundary plane. The upper surface of the layer is loaded by given forces normal to the boundary. The considered problem is solved by using the method of Fourier transforms and dual integral equations, which are reduced to a Fredholm integral equation of the second kind. The stress intensity factor is given in the term of solution of the Fredholm integral equation and some numerical results are presented.     

On the indentation of an inhomogeneous anisotropic elastic material by multiple straight rigid punches

A generalised plane strain problem concerning the indentation of an inhomogeneous anisotropic elastic material by mutiple straight rigid punches is considered. The problem is reduced to a boundary integral equation with the stresses over the contact regions being represented in terms of Chebyschev polynomials. The boundary integral equation is solved numerically for some particular antiplane contact problems involving one or two contact regions and the stress intensity factors at the ends of the contact regions are calculated. The effect of anisotropy and inhomogeneity on the stress intensity factors is examined through the illustrative examples. The analysis is relevant for a class of geomechanics problems involving inhomogeneous materials.     

Contact between an elastically supported circular plate and a rigid indenter

The axially symmetric problem of a finite circular plate loaded at its center by a smooth, rigid punch is solved by superposing an infinite layer elasticity solution with a pure bending plate theory solution. The problem is reduced to dual integral equations, which are further reduced to a single Fredholm integral equation of the second kind. The Fredholm equation is numerically solved and the results are used to compute contact stresses under the indenter as well as the overall load-deflection behavior. The problem is formulated to model a partially fixed edge around the plate's perimeter, and calculations are carried out for the limiting cases of simple supports and complete fixity. Various ratios of plate diameter to plate thickness are studied, and the results are compared to both Hertzian contact theory and standard plate theory.     

A modified scaled boundary finite element method for three‐dimensional dynamic soil‐structure interaction in layered soil

This paper is devoted to the analysis of elastodynamic problems in 3D‐layered systems which are unbounded in the horizontal direction. For this purpose, a finite element model of the near field is coupled to a scaled boundary finite element model (SBFEM) of the far field. The SBFEM is originally based on describing the geometry of a half‐space or full‐space domain by scaling the geometry of the near field / far field interface using a radial coordinate. A modified form of the SBFEM for waves in a 2D layer is also available. None of these existing formulations can be used to describe a 3D‐layered medium. In this paper, a modified SBFEM for the analysis of 3D‐layered continua is derived. Based on the use of a scaling line instead of a scaling centre, a suitable scaled boundary transformation is proposed. The derivation of the corresponding scaled boundary finite element (SBFE) equations in displacement and stiffness is presented in detail. The latter is a nonlinear differential equation with respect to the radial coordinate, which has to be solved numerically for each excitation frequency considered in the analysis. Various numerical examples demonstrate the accuracy of the new method and its correct implementation. These include rigid circular and square foundations embedded in or resting on the surface of layered homogeneous or inhomogeneous 3D soil deposits over rigid bedrock. Hysteretic damping is assumed in some cases. The dynamic stiffness coefficients calculated using the proposed method are compared with analytical solutions or existing highly accurate numerical results. Copyright © 2011 John Wiley & Sons, Ltd.     

The mixed-type integral equations for transient elastodynamic plane crack problems

In the present paper, by use of the boundary integral equation method and the techniques of Green fundamental solution and singularity analysis, the dynamic infinite plane crack problem is investigated. For the first time, the problem is reduced to solving a system of mixed-typed integral equations in Laplace transform domain. The equations consist of ordinary boundary integral equations along the outer boundary and Cauchy singular integral equations along the crack line. The equations obtained are strictly proved to be equivalent with the dual integral equations obtained by Sih in the special case of dynamic Griffith crack problem. The mixed-type integral equations can be solved by combining the numerical method of singular integral equation with the ordinary boundary element method. Further use the numerical method for Laplace transform, several typical examples are calculated and their dynamic stress intensity factors are obtained. The results show that the method proposed is successful and can be used to solve more complicated problems.     

Solution of plane elasticity problems by the displacement discontinuity method. I. Infinite body solution

This paper is concerned with the development of a numerical procedure for solving complex boundary value problems in plane elastostatics. This procedure—the displacement discontinuity method—consists simply of placing N displacement discontinuities of unknown magnitude along the boundaries of the region to be analyzed, then setting up and solving a system of algebraic equations to find the discontinuity values that produce prescribed boundary tractions or displacements. The displacement discontinuity method is in some respects similar to integral equation or ‘influence function’ techniques, and contrasts with finite difference and finite element procedures in that approximations are made only on the boundary contours, and not in the field. The method is illustrated by comparing computed results with the analytical solutions of two boundary value problems: a circular disc subjected to diametral compression, and a circular hole in an infinite plate under a uniaxial stress field. In both cases the numerical results are in excellent agreement with the exact solutions.     

动水压力问题的一种边界元解法

本文将动水压力设定为有限个半平面解答的线性组合,它精确满足域内的方程,然后令解答在N 个边界点满足给定的边界条件求解设定的解答中的待定常数。该方法有简单,避免奇异积分、解设定常数的方程组状态好等优点。文中给出的算例说明了方法的正确性。     

Finite difference analysis for laminar flow heat transfer in concentric annuli with simultaneously developing hydrodynamic and thermal boundary layers

The results of a finite difference analysis are presented for the problem of incompressible laminar flow heat transfer in concentric annuli with simultaneously developing hydrodynamic and thermal boundary layers, the boundary conditions of one wall being isothermal and the other wall adiabatic. This corresponds to the fundamental solution of the third kind according to the four fundamental solutions classified by Reynolds, Lundberg and McCuen¹. Firstly, the hydrodynamic entry length problem, based on the boundary layer simplifications of the Navier–Stokes equations, was solved by means of an extension of the linearized finite difference scheme used previously by Bodia and Osterle² to solve a similar problem between parallel plates. The energy equation is then solved, using the velocity profiles previously obtained, by means of an implicit finite difference technique. The accuracy of the numerical solution was checked by comparing results for the annulus of radius ratio 0.25 with the avaiable solution of Shumway and McEligot³.     

Numerical solution of the Dirichlet initial boundary value problem for the heat equation in exterior 3-dimensional domains using integral equations

A numerical method for the Dirichlet initial boundary value problem for the heat equation in the exterior and unbounded region of a smooth closed simply connected 3-dimensional domain is proposed and investigated. This method is based on a combination of a Laguerre transformation with respect to the time variable and an integral equation approach in the spatial variables. Using the Laguerre transformation in time reduces the parabolic problem to a sequence of stationary elliptic problems which are solved by a boundary layer approach giving a sequence of boundary integral equations of the first kind to solve. Under the assumption that the boundary surface of the solution domain has a one-to-one mapping onto the unit sphere, these integral equations are transformed and rewritten over this sphere. The numerical discretisation and solution are obtained by a discrete projection method involving spherical harmonic functions. Numerical results are included.     

应力边界元法解平面热弹性问题

本文提出了求解平面热弹性问题的应力边界元法。利用应力法由平面热弹性问题的基本方程出发,简要地叙述了边界积分方程的建立,给出了位移单值条件。这种方法适用于应力边界值问题。作为数值计算例,计算了圆形区域和具有偏心圆孔的圆形区域的热应力,得到了满意的结果。应力边界元法也可应用于平板弯曲问题。     

A boundary integral approach for plane analysis of electrically semi-permeable planar cracks in a piezoelectric solid

A plane electro-elastostatic problem involving arbitrarily located planar stress free cracks which are electrically semi-permeable is considered. Through the use of the numerical Green's function for impermeable cracks, the problem is formulated in terms of boundary integral equations which are solved numerically by a boundary element procedure together with a predictor–corrector method. The crack tip stress and electric displacement intensity factors can be easily computed once the boundary integral equations are properly solved.     

Boundary integral equation method for linear porous-elasticity with applications to soil consolidation

For physical phenomena governed by the Biot model of porous-elasticity, a reciprocal relation, similar to the Betti's recoprocal theorem in elasticity, is constructed in Laplace transformed space. Integrating the reciprocal relation enables one to formulate boundary integral equations. The fundamental kernels for the integral equations are solved in closed forms for the case of isotropic material. Numerical implementation of two-dimensional problems includes finite element ideas of discretization and polynomial interpolation, and numerical inversion of a Laplace transform. Practical applications of the method are found in consolidation problems in soils which contain compressible as well as incompressible pore fluids. Also, as a numerical experiment, consolidation of partially saturated soil is simulated and interesting phenomena are observed. The currently developed boundary integral equation method (BIEM) for porous-elasticity may be viewed as an efficient and accurate alternative of existing finite element and finite difference methods. For linear consolidation problems, application of BIEM is always preferred to the other numerical methods whenever possible.     

Complex hypersingular integral equation for the piece-wise homogeneous half-plane with cracks

New complex hypersingular integral equation (CHSIE) is derived for the half-plane containing the inclusions (which can have the different elastic properties), holes, notches and cracks of the arbitrary shape. This equation is obtained by superposition of the equations for each homogeneous region in a half-plane. The last equations follow from the use of complex analogs of Somigliana's displacement and stress identities (SDI and SSI) and Melan's fundamental solution (FS) written in a complex form. The universal numerical algorithm suggested before for the analogous problem for a piece-wise homogeneous plane is extended on case of a half plane. The unknown functions are approximated by complex Lagrange polynomials of the arbitrary degree. The asymptotics for the displacement discontinuities (DD) at the crack tips are taken into account. Only two types of the boundary elements (straight segments and circular arcs) are used to approximate the boundaries. All the integrals involved in CHSIE are evaluated in a closed form. A wide range of elasticity problems for a half-plane with cracks, openings and inclusions are solved numerically.     

Numerical treatment of finite part integrals in 2-D boundary element analysis with application infracture mechanics

For the solution of problems in fracture mechanics by the boundary element method usually the subregion technique is employed to decouple the crack surfaces. In this paper a different procedure is presented. By using the displacement boundary integral equation on one side of the crack surface and the hypersingular traction boundary integral equation on the opposite side, one can renounce the subregion technique.An essential point when applying the traction boundary integral equation is the treatment of the thus arising hypersingular integrals. Two methods for their numerical computation are presented, both based on the finite part concept. One may either scale the integrals properly and use a specific quadrature rule, or one may apply the definition formula for finite part integrals and transform the resulting regular integrals into the usual element coordinate system afterwards. While the former method is restricted to linear or circular approximations of the boundary geometry, the latter one allows for arbitrary curved (e.g. isoparametric) elements. Two numerical examples are enclosed to demonstrate the accuracy of the two boundary integral equations technique compared with the subregion technique.     

The axisymmetric double contact problem for a frictionless elastic layer indented by an elastic cylinder

This paper is concerned with the smooth receding contact between an elastic layer and a half space when the layer is compressed by a frictionless semi-infinite elastic cylinder. Upon loading, the contact along the layer-subspace interface shrinks to a circular area, radius of which is unknown. The analysis leads to a system of singular integral equations of the second kind. The integral equations are solved numerically and the contact pressures, extent of contact and the stress intensity factor round the edge of the cylinder are calculated for various material combinations.     

The mode I crack problem for layered piezoelectric plates

The plane strain singular stress problem for piezoelectric composite plates having a central crack is considered. For the case of the crack which is normal to and ends at the interface between the piezoelectric plate and the elastic layer, the order of stress singularity around the tip of the crack is obtained. The Fourier transform technique is used to formulate the problem in terms of a singular integral equation. The singular integral equation is solved by using the Gaus–Jacobi integration formula. Numerical calculations are carried out, and the main results presented are the variation of the stress intensity factor as functions of the geometric parameters, the piezoelectric material properties and the electrical boundary conditions of the layered composites.     

A receding contact problem between a functionally graded layer and two homogeneous quarter planes

In this paper, the plane problem of a frictionless receding contact between an elastic functionally graded layer and two homogeneous quarter planes is considered when the graded layer is pressed against the quarter planes. The top of the layer is subjected to normal tractions over a finite segment. The graded layer is modeled as a non-homogeneous medium with a constant Poisson’s ratio and exponentially varying shear modules. The problem is converted into the solution of a Cauchy-type singular integral equation in which the contact pressure and the receding contact half-length are the unknowns using integral transforms. The singular integral equation is solved numerically using Gauss–Jacobi integration. The corresponding receding contact half-length that satisfies the global equilibrium condition is obtained using an iterative procedure. The effect of the material non-homogeneity parameter on the contact pressure and on the length of the receding contact is investigated.