Stress singularities for three-dimensional corners using the boundary integral equation method

The boundary integral equation method is developed to study three-dimensional asymptotic singular stress fields at vertices of a pyramidal notch or inclusion in an isotropic elastic space. Two-dimensional boundary integral equations are used for the infinite body with pyramidal notches and inclusions when either stresses or displacements are specified on its surface. Applying the Mellin integral transformation reduces the problem to one-dimensional singular integral equations over a closed, piece-wise smooth line. Using quadrature formulas for regular and singular integrals with Hilbert and logarithmic kernels, these integral equations are reduced to a homogeneous system of linear algebraic equations. Setting its determinant to zero provides a characteristic equation for the determination of the stress singularity power. Numerical results are obtained and compared against known eigenvalues from the literature for an infinite region with a conical notch or inclusion, for a Fichera vertex, and for a half-space with a wedge-shaped notch or inclusion.     

Singular integral equations in elastic dielectrics

Galerkin representations for the displacement vector, polarization vector and the potential field are obtained by elementary matrix inversions of the equations of equilibrium. Matrices of fundamental solutions of an infinite elastic dielectric continuum subjected to a concentrated body force, an electric force, and a charge density, are constructed. Theorems are proved on the discontinuity of double layer potentials and R, M, M operators of single layer potentials. By means of these theorems, the solution of the two basic boundary value problems has been reduced to the solution of a system of seven singular integral equations.     

Analysis of the Singularity for the Vertical Contact Problem of Line Crack-Inclusion

ＩｎｔｒｏｄｕｃｔｉｏｎＵｐｔｏｎｏｗ ,ｔｈｅｔｅｃｈｎｉｃａｌｌｉｔｅｒａｔｕｒｅｏｎｓｅｐａｒａｔｅｃｒａｃｋｓ,ｖｏｉｄｓ,ｉｎｃｌｕｓｉｏｎｓａｎｄｔｈｅｉｎｔｅｒａｃｔｉｏｎｓｂｅｔｗｅｅｎｃｒａｃｋｓａｎｄｉｎｃｌｕｓｉｏｎｓｈａｖｅｂｅｅｎｑｕｉｔｅｅｘｔｅｎｓｉｖｅ.Ｈｏｗｅｖｅｒ,ｔｈｅｃｏｎｔａｃｔｐｒｏｂｌｅｍｓｏｆｃｒａｃｋ＿ｉｎｃｌｕｓｉｏｎｄｏｎｏｔｓｅｅｍｔｏｂｅａｓｗｉｄｅｌｙｓｔｕｄｉｅｄ .Ｔｈｉｓｐａｐｅｒｃａｎｂｅｒｅｇａｒｄｅｄａｓｔｈｅｆｕｒｔｈｅｒｒｅｓｅ…     

Equilibrium of layers weakened by loaded cuts

In this paper, a symmetric boundary value problem of the stress analysis for an equilibrium of layer with end-supports covered by diaphragms and weakened by several loaded stress raisers is investigated. The given boundary value problem is reduced to an infinite system of singular integral equations of the second kind. The expressions for stress components in an elastic layer weakened by stress raisers are presented. Based on the developed analytical procedure, extensive numerical investigations have been conducted. The results of these investigations are illustrated graphically exposing some novel qualitative and quantitative knowledge about stress concentration in the layer depending on some geometric parameters of stress raisers and Poisson’s ratio of a layer material.     

Green's function for thermopiezoelectric plates with holes of various shapes

Summary Thermoelectroelastic problems for holes of various shapes embedded in an infinite matrix are considered in this paper. Based on Lekhnitskii's formalism, the technique of conformal mapping and the exact electric boundary conditions on the hole boundary, the thermoelectroelastic Green's function has been obtained analytically in terms of a complex potential. As an application of the proposed function, the problem of an infinite plate containing a crack and a hole is analysed. A system of singular integral equations for the unknown temperature discontinuity and the discontinuity of elastic displacement and electric potential (EDEP) defined on crack faces is developed and solved numerically. Numerical results for stress and electric displacement (SED) intensity factors of the crack-hole system are presented to illustrate the application of the proposed formulation. Received 7 October 1998; accepted for publication 26 January 1999     

Bending of layers weakened by cuts

In the present paper, a method proposed by one of the authors is extended to a class of skew-symmetric elastic problems for the stress analysis of a layer supported by sliding fixed supports and weakened by several stress raisers. The corresponding boundary value problem is reduced to an infinite system of one-dimensional singular integral equations of the second kind. The expressions for the stress components in an elastic layer weakened by stress raisers are presented. Based on the developed analytical algorithm, extensive numerical investigations have been conducted. The results of these investigations are illustrated graphically exposing some novel qualitative and quantitative knowledge about stress concentration in the layer depending on some geometric parameters of stress raisers and Poisson’s ratio of a layer material.     

垂直磁压电材料界面三维裂纹的超奇异积分法

应用有限部积分概念和广义位移基本解，垂直于磁压电双材料界面三维复合型裂纹问题被转 化为求解一组以裂纹表面广义位移间断为未知函数的超奇异积分方程问题. 进而，通过主部 分析法精确地求得裂纹尖端光滑点附近的奇性应力场解析表达式. 然后，通过将裂纹表面 位移间断未知函数表达为位移间断基本密度函数与多项式之积，使用有限部积分法对超奇异 积分方程组建立了数值方法. 最后，通过典型算例计算，讨论了广义应力强度因子的变化规 律.     

Hypersingular integral equation method for a three-dimensional crack in anisotropic electro-magneto-elastic bimaterials

Using Green’s functions, the extended general displacement solutions of a three-dimensional crack problem in anisotropic electro-magneto-elastic (EME) bimaterials under extended loads are analyzed by the boundary element method. Then, the crack problem is reduced to solving a set of hypersingular integral equations (HIE) coupled with boundary integral equations. The singularity of the extended displacement discontinuities around the crack front terminating at the interface is analyzed by the main-part analysis method of HIE, and the exact analytical solutions of the extended singular stresses and extended stress intensity factors (SIFs) near the crack front in anisotropic EME bimaterials are given. Also, the numerical method of the HIE for a rectangular crack subjected to extended loads is put forward with the extended crack opening dislocation approximated by the product of basic density functions and polynomials. At last, numerical solutions of the extended SIFs of some examples are obtained.     

The Problem of an External Circular Crack Under Asymmetric Loadings

ＩｎｔｒｏｄｕｃｔｉｏｎＷｈｅｎｔｈｅｂｏｕｎｄａｒｙｉｎｔｅｇｒａｌｅｑｕａｔｉｏｎｍｅｔｈｏｄｉｓａｐｐｌｉｅｄｔｏｃｒａｃｋａｎａｌｙｓｉｓ,ｏｎｌｙｎｕｍｅｒｉｃａｌｓｏｌｕｔｉｏｎｓｃａｎｂｅｏｂｔａｉｎｅｄ ,ｓｕｃｈａｓ:ｔｈｅｔｙｐｉｃａｌｗｏｒｋｓｏｆＳｎｙｄｅｒａｎｄＣｒｕｓｅ[1],Ｃｒｏｕｃｈ[2 ],Ｂｌａｎｄｆｏｒｄｅｔａｌ.[3],Ｐｏｒｔｅｌａｅｔａｌ.[4 ],Ｂｕｉ[5 ],Ｗｅａｖｅｒ[6 ]ａｎｄＷＡＮＧｅｔａｌ.[7- 9].Ｉｔｉｓｖｅｒｙｄｉｆｆｉｃｕｌｔｔｏａｐｐｌｙｔｈｅｂｏｕｎｄａｒ…     

横观各向同性材料的三维断裂力学问题

从三维横观各向同性材料弹性力学理论出发, 使用Hadamard有限部积分概念, 导出了三维状态下单位位移间断(位错)集度的基 本解. 在此基础上, 进一步运用极限理论, 将任意载荷作用下, 三维无限大横观各向 同性材料弹性体中, 含有一个位于弹性对称面内的任意形状的片状裂纹问题, 归结为求 解一组超奇异积分方程的问题. 通过二维超奇异积分的主部分析方法, 精确地求得了裂纹前沿光滑点附近的应力奇异指数和奇异应力场, 从而找到了以裂纹表面位移间断表示的应力强度因子表达式及裂纹局部扩展所提供 的能量释放率. 作为以上理论的实际应用，最后给出了一个圆形片状裂纹问题 的精确解例和一个正方形片状裂纹问题的数值解例. 对受轴对称法向均布载荷作用下圆形片状裂纹问题, 讨论了超奇异积分方程的精确求解方法, 并获得了位移间断和应力强度因子的封闭解, 此结果与现有理论解完全一致.     

Transverse cracking of orthotropic composite laminates: a fracture mechanics approach

This research is concerned with the fracture mechanics of a laminated composite medium, which contains a central layer sandwiched by two outer layers. There is a periodic array of cracks in the central layer along the central axis of the medium. Fourier transform is used to reduce the problem to the solution of a system of dual integral equations, which are solved by the singular integral equation technique. Rigorous fracture mechanics analysis, which exactly satisfies all boundary conditions of the problem, is conducted. Numerical solutions for the crack tip field and the stress in the medium are obtained for various values such as crack length, crack spacing and layer thickness. Results are also given for the reduction of the equivalent Young’s modulus of the laminate due to multiple cracking. The cases of axial extension and residual temperature change of the composite medium are accounted for.     

三维间断位移法及强奇异和超奇异积分的处理方法

从积分方程Somigliana等式出发,导出三维状态下单位位错集度的基本解.在此基础上,建立了边界积分方程,并给出了其离散形式.对强奇异和超奇异积分,采用了Hadamard定义的有限部分积分来处理.最后,给出了计算裂纹应力强度因子的算例,并与解析解进行了比较,证实了该方法的有效性.     

Cracked elastic layer under compressive mechanical loads

We consider boundary value problem in which an elastic layer containing a finite length crack is under compressive loading. The crack is parallel to the layer surfaces and the contact between crack surfaces are either frictionless or with adhesive friction or Coulomb friction.Based on fourier integral transformation techniques the solution of the formulated problems is reduced to the solution of a singular integral equation, then, using Chebyshev’s orthogonal polynomials, to an infinite system of linear algebraic equations. The regularity of these equations is established. The expressions for stress and displacement components in the elastic layer are presented. Based on the developed analytical algorithm, extensive numerical investigations have been conducted.The results of these investigations are illustrated graphically, exposing some novel qualitative and quantitative knowledge about the stress field in the cracked layer and their dependence on geometric and applied loading parameters. It can be seen from this study that the crack tip stress field has a mode II type singularity.     

同震断层三维裂纹扩展波动时域超奇异积分法

应用波动时域超奇异积分法将P波、S波和磁电热弹多场耦合作用下同震断层任意形状三维裂纹扩展问题转化为求解以广义位移间断率为未知函数的超奇异积分方程组问题;定义了广义应力强度因子,得到裂纹前沿广义奇异应力增量解析表达式;应用波动时域有限部积分概念及体积力法,为超奇异积分方程组建立了数值求解方法,编制了FORTRAN程序,以三维矩形裂纹扩展问题为例,通过典型算例,研究了广义应力强度因子随裂纹位置变化规律;分析了同震断层裂纹扩展中力、磁、电场辐射规律.      

Interaction of general plane P wave and cylindrical inclusion partially debonded from its viscoelastic matrix

The interaction of a general plane P wave and an elastic cylindrical inclusion of infinite length partially debonded from its surrounding viscoelastic matrix of infinite extension is investigated. The debonded region is modeled as an arc-shaped interface crack between inclusion and matrix with non-contacting faces. With wave functions expansion and singular integral equation technique, the interaction problem is reduced to a set of simultaneous singular integral equations of crack dislocation density function. By analysis of the fundamental solution of the singular integral equation, it is found that dynamic stress field at the crack tip is oscillatory singular, which is related to the frequency of incident wave. The singular integral equations are solved numerically, and the crack open displacement and dynamic stress intensity factor are evaluated for various incident angles and frequencies. The project supported by the National Natural Science Foundation of China (19872002) and Climbing Foundation of Northern Jiaotong University     

导数场边界积分方程分析近边界应力分布

研究二维弹性力学问题边界积分方程,通过分部积分变换消除了常规导数边界积分方程中的超奇异积分,获得仅含强奇异积分的应力自然边界积分方程.对于近边界应力的计算,进一步运用正则化算法解析计算其中的几乎强奇异积分.较常规边界元法相比,应力自然边界积分方程可以求解离边界更加接近的内点应力值.算例证明了文中方法的可应用性和有效性.     

On the dynamic behaviour of a piezoelectric laminate with multiple interfacial collinear cracks

This paper investigates the dynamic behaviour of a piezoelectric laminate containing multiple interfacial collinear cracks subjected to steady-state electro-mechanical loads. Both the permeable and impermeable boundary conditions are examined and discussed. Based on the use of integral transform techniques, the problem is reduced to a set of singular integral equations, which can be solved using Chebyshev polynomial expansions. Numerical results are provided to show the effect of the geometry of interacting collinear cracks, the applied electric fields, the electric boundary conditions along the crack faces and the loading frequency on the resulting dynamic stress intensity and electric displacement intensity factors.     

Application of hypersingular integral equation method to three-dimensional crack in electromagnetothermoelastic multiphase composites

A three-dimensional crack problem in electromagnetothermoelastic multiphase composites (EMTE-MCs) under extended loads is investigated in this paper. Using Green’s functions, the extended general displacement solutions are obtained by the boundary element method. This crack problem is reduced to solving a set of hypersingular integral equations coupled with boundary integral equations, in which the unknown functions are the extended displacement discontinuities. Then, the behavior of the extended displacement discontinuities around the crack front terminating at the interface is analyzed by the main-part analysis method of hypersingular integral equations. Analytical solutions for the extended singular stresses, the extended stress intensity factors (SIFs) and the extended energy release rate near the crack front in EMTE-MCs are provided. Also, a numerical method of the hypersingular integral equations for a rectangular crack subjected to extended loads is put forward with the extended displacement discontinuities approximated by the product of basic density functions and polynomials. In addition, distributions of extended SIFs varying with the shape of the crack are presented. The results show that the present method accurately yields smooth variations of extended SIFs along the crack front.     

Boundary integral equations of unique solutions in elasticity

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｂｏｕｎｄａｒｙｅｌｅｍｅｎｔｍｅｔｈｏｄ（ＢＥＭ）ｐｒｏｖｉｄｅｓａｎａｔｔｒａｃｔｉｖｅａｌｔｅｒｎａｔｉｖｅｆｏｒｔｈｅａｎａｌｙｓｉｓｏｆｅｎｇｉｎｅｅｒｉｎｇｐｒｏｂｌｅｍｓ．Ｉｔｓｍａｉｎａｄｖａｎｔａｇｅｓａｒｅｅｃｏｎｏｍｉｃａｌａｎｄｐａｒｔｉｃｕｌａｒｌｙｃｏｎｖｅｎｉｅｎｔｆｏｒｕｎｂｏｕｎｄｅｄｄｏｍａｉｎａｎｄｓｔｒｅｓｓｃｏｎｃｅｎｔｒａｔｉｏｎｐｒｏｂｌｅｍｓ．Ｔｈｅｂｏｕｎｄａｒｙｉｎｔｅｇｒａｌｅｑｕａｔｉｏｎ（ＢＩＥ）ｉｓｔｈｅ…     

Analysis of three-dimensional edge cracks under tensile loading

Three-dimensional edge cracks are analyzed using the Self-Similar Crack Expansion (SSCE) method with a boundary integral equation technique. The boundary integral equations for surface cracks in a half space are presented based on a half space Green's function (Mindlin, 1936). By using the SSCE method, the stress intensity factors are determined by the crack-opening displacement over the crack surface. In discrete boundary integral equations, the regular and singular integrals on the crack surface elements are evaluated by an analytical method, and the closed form expressions of the integrals are given for subsurface cracks and edge crakcs. This globally numerical and locally analytical method improves the solution accuracy and computational effort. Numerical results for edge cracks under tensile loading with various geometries, such as rectangular cracks, elliptical cracks, and semi-circular cracks, are presented using the SSCE method. Results for stress intensity factors of those surface breaking cracks are in good agreement with other numerical and analytical solutions.