Axial soil-pile interaction in a transversely isotropic half-space

A rigorous integral equation formulation is presented for the axisymmetric load-transfer analysis of a thin-walled pile embedded in a transversely isotropic half-space under axial load. By virtue of a set of ring-load Green’s functions for the pile and one for the half-space, the problem is shown to be reducible to a pair of Fredholm integral equations. Through a mathematical analysis of an auxiliary pair of Cauchy integral equations, the inherent singularities of the contact stress distributions are rendered explicit. With a direct incorporation of the singular nature of the resultant load transfers, the numerical solution of the integral equations is shown to be possible by an interpolation of regular functions only. Typical results for various material and geometrical conditions are presented, including a comparison with past classical solutions, to illustrate the effects of transverse anisotropy on the load-transfer process.     

Free vibration of four-parameter functionally graded spherical and parabolic shells of revolution with arbitrary boundary conditions

The objective of this work is to present a Haar Wavelet Discretization (HWD) method-based solution approach for the free vibration analysis of functionally graded (FG) spherical and parabolic shells of revolution with arbitrary boundary conditions. The first-order shear deformation theory is adopted to account for the transverse shear effect and rotary inertia of the shell structures. Haar wavelet and their integral and Fourier series are selected as the basis functions for the variables and their derivatives in the meridional and circumferential directions, respectively. The constants appearing in the integrating process are determined by boundary conditions, and thus the equations of motion as well as the boundary condition equations are transformed into a set of algebraic equations. The proposed approach directly deals with nodal values and does not require special formula for evaluating system matrices. Also, the convenience of the approach is shown in handling general boundary conditions. Numerical examples are given for the free vibrations of FG shells with different combinations of classical and elastic boundary conditions. Effects of spring stiffness values and the material power-law distributions on the natural frequencies of shells are also discussed. Some new results for the considered shell structures are presented, which may serve as benchmark solutions.     

Transverse impact of a ball on a sphere with allowance for waves in the target

The transverse impact of a solid projectile on an elastic spherical shell with a pivoting contour support has been studied. Inside the contact zone, the projectile-target interaction is described by a solution of the standard system of equations. Outside the contact zone, the points of the shell are displaced and the shell is deformed due to propagation of a nonstationary wave front. A solution in this region is constructed using ray series with variable coefficients representing jumps of the time derivatives of the unknown functions on the wave surface of strong discontinuity. These coefficients are determined to within arbitrary constants using momentless equations of motion of the shell points. The constants are determined by matching two solutions at the contact zone boundary. Using the obtained analytical expressions and plotted dependences for the contact force and dynamic inflection, it is possible to judge on the influence of the shell structure design on the dynamic characteristics of impact interaction.     

Boundary integral equation method for conductive cracks in two and three-dimensional transversely isotropic piezoelectric media

Using the fundamental solutions and the Somigliana identity of piezoelectric medium, the boundary integral equations are obtained for a conductive planar crack of arbitrary shape in three-dimensional transversely isotropic piezoelectric medium. The singular behaviors near the crack edge are studied by boundary integral equation approach, and the intensity factors are derived in terms of the displacement discontinuity and the electric displacement boundary value sum near the crack edge on crack faces. The boundary integral equations for two dimensional crack problems are deduced as a special case of infinite strip planar crack. Based on the analogy of the obtained boundary integral equations and those for cracks in conventional isotropic elastic material and for contact problem of half-space under the action of a rigid punch, an analysis method is proposed. As an example, the solution to conductive Griffith crack is derived.     

On the method of functional equations and the performance of desingularized boundary element methods

A correspondence is made between the reciprocal relation for linear elliptic partial differential equations and the Riesz integral representation. The former relates the boundary distributions and appropriate normal fluxes of two arbitrary solutions, and the latter expresses a continuous linear functional in terms of an integral involving a representing function. When sufficient regularity conditions are met, the representing function is identified with the unknown boundary distribution. In principle, the representing function may be expressed in terms of the images of a complete set of orthonormal basis functions with known normal fluxes, as suggested by Kupradze [Kupradze VD. On the approximate solution of problems in mathematical physics. Russ Math Surv 1967; 22: 59–107]; in practice, the representing function is computed by solving integral equations using boundary element methods. The basic procedure involves expressing the representing function in terms of finite-element or other basis functions, and requiring the satisfaction of the reciprocal relationship with a suitable set of test functions such as Green's functions and their dipoles. When the singular points are placed at the boundary, we obtain the standard boundary integral equation method. When the singular points are placed outside the domain of solution, we obtain a system of functional equations and associated class of desingularized boundary integral methods. When sufficient regularity conditions are met and the test functions comprise a complete set, then in the limit of infinite discretization the numerical solution converges to the unknown boundary distribution. An overview of formulations is presented with reference to Laplace's equation in two dimensions. Numerical experimentation shows that, in general, the solution obtained by desingularized methods becomes increasingly less accurate as the singular points of Green's functions move farther away from the boundary, but the loss of accuracy is significant only when the exact solution shows pronounced variations. Exceptions occur when the integral equation does not have a unique solution. In contrast, and in agreement with previous findings, the condition number of the linear system increases rapidly with the distance of the singular points from the boundary, to the extent that a dependable solution cannot be obtained when the singularities are located even a moderate distance away from the boundary. The desingularized formulation based on Green's function dipoles is superior in accuracy and reliability to the one that uses Green functions. The implementation of the method to the equations of elastostatics and Stokes flow are also discussed.     

Integral equation method for soulution of boundary value problems of structural mechanics,Part I. Ordinary differential equations

This paper presents a method of numerical solution of boundary value problems governed by a set of ordinary differential equations. The highest derivative is chosen as the unknown function. Governing equations are transformed into a set of integral equations. The kernels of integral equations turn out to be influence functions for deflection and/or bending moments of a corresponding beam, and can, therefore, be computed using well-known methods of structural analysis. Finally, the unknown function (highest derivative) is approximated in a defined manner and the solution is obtained through numerical integration.     

Implementation of thermoelastic forces in boundary element analysis of elastic contact and fracture mechanics problems

This paper presents a pseudo-body-force approach multi-domain boundary integral equation method for the analysis of thermoelastic and body-force type elastic contact and fracture mechanics problems. Using this approach only the boundaries of the bodies involved have to be discretized. The transformation of the domain integrals due to body-force and pseudo-force to their equivalent boundary integrals are shown. Also, it is shown that by employing the initial strain approach the same set of equivalent boundary integrals would be obtained. Isoparametric quadratic elements are employed to represent the geometries and the functions. This two-dimensional BEM thermoelastic implementation can be found very simple and can be applied to both harmonic and nonharmonic temperature distributions. The accuracy is asserted by applying it to several thermoelastic fracture mechanics and contact problems.     

Explicit evaluation of hypersingular boundary integral equations for acoustic sensitivity analysis based on direct differentiation method

This paper presents a new set of boundary integral equations for three dimensional acoustic shape sensitivity analysis based on the direct differentiation method. A linear combination of the derived equations is used to avoid the fictitious eigenfrequency problem associated with the conventional boundary integral equation method when solving exterior acoustic problems. The strongly singular and hypersingular boundary integrals contained in the equations are evaluated as the Cauchy principal values and Hadamard finite parts for constant element discretization without using any regularization technique in this study. The present boundary integral equations are more efficient to use than the usual ones based on any other singularity subtraction technique and can be applied to the fast multipole boundary element method more readily and efficiently. The effectiveness and accuracy of the present equations are demonstrated through some numerical examples.     

Boundary integral equations and the boundary element method for fracture mechanics analysis in 2D anisotropic thermoelasticity

This paper develops the Somigliana type boundary integral equations for fracture of anisotropic thermoelastic solids using the Stroh formalism and the theory of analytic functions. In the absence of body forces and internal heat sources, obtained integral equations contain only curvilinear integrals over the solid’s boundary and crack faces. Thus, the volume integration is eliminated and also there is no need to evaluate integrals over the contours in the mapped temperature domain as it was done before. In addition to finite solids, the case of an infinite anisotropic medium with a remote thermal load is also studied. The dual boundary element method for fracture of anisotropic thermoelastic solids is developed based on the obtained boundary integral equations. Presented numerical examples show the validity and efficiency of the obtained equations in the analysis of both finite and infinite solids with cracks.     

Buckling analysis of shear deformable shallow shells by the boundary element method

In this work a boundary element (BE) formulation for buckling problem of shear deformable shallow shells is presented. A set of five boundary integral equations are obtained by coupling two-dimensional plane stress elasticity with shear deformable plate bending (Reissner). The domain integrals appearing in the formulation (due to the curvature and due to the domain load) are transferred into equivalent boundary integrals. The BE formulation is presented as an eigenvalue problem, to provide direct evaluation of critical load factors and buckling modes. Several examples are presented. The BE results for a cylindrical shallow shell with different curvatures are compared with other numerical solutions and good agreements are obtained.     

An elastic contact problem for a neo-hookean half-space indented by a annular rigid stamp of arbitrary profiles

The present paper deals with an elastic contact problem for a neo-Hookean half-space indented by a annular rigid stamp of arbitrary profile, which constitutes a three part mixed boundary value problem. The problem is reduced into triple integral equations which are further reduced into an infinite set of simultaneous equations. Numerical results are illustrated graphically for the distributions of the normal displacement and stress in the case of parabolic concave punch for different values of a/b at $\text{z = 0}$     

Levi functions for linear elliptic systems with variable coefficients including shell equations

This paper gives an algorithm to construct Levi functions of arbitrary degree for elliptic systems of linear partial differential equations with variable (real-analytic) coefficients. Further, an indirect method is described to transform elliptic boundary value problems into a system of integral equations. This method is applied to the shell equations in the non-shallow case. (In the shallow case the shell equations have constant coefficients.) Some questions of discretization are discussed and numerical results are presented.     

离散多层圆筒在热冲击载荷下的弹性动力响应

离散多层圆筒由薄内筒和倾角错绕的钢带层组成,具有制造简便、成本低等优点。预测筒体在热冲击载荷下的热应力对强度设计和安全操作具有重要的应用价值。该文首次研究了离散多层圆筒在热冲击载荷作用下的热弹性动态响应。将内筒和钢带层的径向位移分别分解为满足给定应力边界条件的准静态解和满足初始条件的动态解,准静态解通过齐次线性方法确定,热弹性动态解通过有限Hankel积分变换和Laplace变换确定。根据内外层界面处位移连续条件,得到层间压力关于时间的第二类Volterra积分方程,利用Hermit二次三项式插值方法可求得该层间应力。最后将离散多层圆筒的热弹性动力响应与单层厚壁圆筒的响应进行了比较,并分析了钢带缠绕倾角和材料参数对热弹性动力响应的影响。     

Complex variable boundary integral method for linear viscoelasticity: Part I—basic formulations

The basic formulations (direct and indirect) of the complex variable boundary integral method for linear viscoelasticity are presented. Complex variable temporal integral equations for the formulations are obtained for viscoelastic solids whose behavior in shear is governed by a Boltzmann model while the bulk behavior is purely elastic. The functions involved in the integral equations are the time-dependent complex boundary tractions and displacements for the direct approach and the unknown time-dependent complex density functions for the indirect approaches. The temporal integral equations give the displacements and stresses at a point inside a viscoelastic region in terms of time convolution and space integrals over the boundary of this region. The equations are valid for the boundaries of arbitrary shapes provided that these boundaries are sufficiently smooth. Complex variable temporal boundary equations are obtained by taking the inner point to the boundary. Numerical treatment of spatial and time convolution integrals involved in the boundary equations is discussed.     

Forced torsional oscillations of multilayered solids

A new approach is proposed for obtaining the dynamic elastic response of a multilayered elastic solid caused by a forced torsional oscillation inside the solid. The elastodynamic Green’s function of the center of rotation and a point load method are used to solve the problem. The solution of the center of rotation for multilayered solids is obtained by solving a set of simultaneous linear algebraic equations using the boundary conditions for the singularity and for the layer interfaces. The solution of the forced torsional oscillation is formulated by integrating the Green’s function over the contact area with unknown surface traction. The dual integral equations of the unknown surface traction are established by considering the boundary conditions on the contact surface of the multilayered solid, which can be converted into a Fredholm integral equation of the second kind and solved numerically.     

Transient behavior of a few dies on an elastic half-space

Summary An asymptotic approach to dynamic interaction between a few distant dies and an elastic half-space is proposed. The transient motion of the dies under low-frequency vertical load is under consideration. The explicit expression for the fundamental singular solution of Lamb's problem is used to derive the boundary integral equation of contact. Then this equation is asymptotically simplified and solved numerically in combination with equations of motion of the dies.Equations obtained in the asymptotic limit describe both the die-medium dynamic interaction and the interaction between dies through the elastic medium. These equations take into account the energy dissipation phenomenon associated with energy transfer deep into the medium by outgoing elastic waves, of so called geometrical damping.Equations proposed are asymptotically correct within the corresponding range of parameters, as such improving the state-of-the-art.     

A comprehensive study on the 2D boundary element method for anisotropic thermoelectroelastic solids with cracks and thin inhomogeneities

This paper develops Somigliana type boundary integral equations for 2D thermoelectroelasticity of anisotropic solids with cracks and thin inclusions. Two approaches for obtaining of these equations are proposed, which validate each other. Derived boundary integral equations contain domain integrals only if the body forces or distributed heat sources are present, which is advantageous comparing to the existing ones. Closed-form expressions are obtained for all kernels. A model of a thin pyroelectric inclusion is obtained, which can be also used for the analysis of solids with impermeable, permeable and semi-permeable cracks, and cracks with an imperfect thermal contact of their faces. The paper considers both finite and infinite solids. In the latter case it is proved, that in contrast with the anisotropic thermoelasticity, the uniform heat flux can produce nonzero stress and electric displacement in the unnotched pyroelectric medium due to the tertiary pyroelectric effect. Obtained boundary integral equations and inclusion models are introduced into the computational algorithm of the boundary element method. The numerical analysis of sample and new problems proved the validity of the developed approach, and allowed to obtain some new results.     

Analysis of dynamic interaction between an inclusion and a nearby moving crack by BEM

In this paper, the dynamic interaction between an inclusion and a nearby moving crack embedded in an elastic medium is studied by the boundary element method (BEM). To deal with this problem, the multi-region technique and two kinds of time-domain boundary integral equations (BIEs) are introduced. The system is divided into two parts along the interface between the inclusion and the matrix medium. Each part is linear, elastic, homogeneous and isotropic. The non-hypersingular traction boundary integral equation is applied on the crack surfaces; while the traditional displacement boundary integral equation is used on the interface and external boundaries. In the numerical solution procedure, square root shape functions are adopted as to describe the proper asymptotic behavior in the vicinity of the crack-tips. The crack growth is modeled by adding new elements of constant length to the moving crack tip, which is controlled by the fracture criterion based on the maximum circumferential stress. In each time step, the direction and the speed of the crack advance are evaluated. The numerical results of the crack growth path, speed, dynamic stress intensity factors (DSIFs) and dynamic interface tractions for various material combinations and geometries are presented. The effect of the inclusion on the moving crack is discussed.     

加层电磁弹性材料界面裂纹瞬态响应

研究加层电磁弹性材料界面裂纹在反平面剪切冲击载荷和面内电磁冲击载荷作用下的动态响应问题。假设裂纹面是电磁不导通的。采用Laplace变换、Fourier变换和位错密度函数将混合边值问题转化为求解Laplace域内Cauchy奇异积分方程。讨论了磁冲击载荷、电冲击载荷、材料参数及加层厚度对能量释放率的影响。该问题的解有助于分析含裂纹电磁弹性材料的动态断裂特性。