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.     

Isotropic/orthotropic two-layered strips containing cracks terminating at and going through an interface

Crack problems for isotropic/orthotropic two-layered strips have been investigated. A system of two singular integral equations can be derived by using Fourier integral transformation and boundary conditions of crack problems. After stress singularities at crack tips or other special points are determined for internal and edge cracks, and for cracks terminating at and going through the interface, the system of singular integral equations is solved numerically by Gauss-Jacobi or Gauss-Chebyshev integration formulas for stress intensity factors at the tips and other singular points of cracks. Finally, possible crack growth behavior for cracks approaching and going through the interface is discussed.     

Dynamic stress intensity factors for two parallel interface cracks between a nonhomogeneous bonding layer and two dissimilar elastic half-planes subject to an impact load

Some composite materials are constructed of two dissimilar half-planes bonded by a nonhomogeneous elastic layer. In the present study, a crack is situated at the interface between the upper half-plane and the bonding layer of such a material, and another crack is located at the interface between the lower half-plane and the bonding layer. The material properties of the bonding layer vary continuously from those of the lower half-plane to those of the upper half-plane. Incoming shock stress waves impinge upon the two interface cracks normal to their surfaces. Fourier transformations were used to reduce the boundary conditions for the cracks to two pairs of dual integral equations in the Laplace domain. To solve these equations, the differences in the crack surface displacements were expanded in a series of functions that are zero-valued outside the cracks. The unknown coefficients in the series were solved using the Schmidt method so as to satisfy the conditions inside the cracks. The stress intensity factors were defined in the Laplace domain and were inverted numerically to physical space. Dynamic stress intensity factors were calculated numerically for selected crack configurations.     

3D dynamic stress intensity factors around two parallel square cracks in an infinite elastic medium under impact load

Summary  Transient stresses around two parallel cracks in an infinite elastic medium are investigated in the present paper. The shape of the cracks is assumed to be square. Incoming shock stress waves impinge upon the two cracks normal to tzheir surfaces. The mixed boundary value equations with respect to stresses and displacements are reduced to two sets of dual integral equations in the Laplace transform domain using the Fourier transform technique. These equations are solved by expanding the differences in the crack surface displacements in a double series of a function that is equal to zero outside the cracks. Unknown coefficients in the series are calculated using the Schmidt method. Stress intensity factors defined in the Laplace transform domain are inverted numerically to the physical space. Numerical calculations are carried out for transient dynamic stress intensity factors under the assumption that the shape of the upper crack is identical to that of the lower crack. Received 2 February 2000; accepted for publication 10 May 2000     

Study into the Interaction of Cracks in an Elastic Half-Space under a Shock Load by Means of Boundary Integral Equations

The effect of a shock load on the interaction of circular cracks in an elastic half-space is studied. In the space of Fourier time transforms, the problem is reduced to a system of two-dimensional boundary integral equations in the form of the Helmholtz potential with unknown densities characterizing the discontinuities in the displacements of the opposite crack faces. Discrete analogs of those equations are constructed. As an example, two cracks are considered whose faces are under the action of shock tensile loads varying in time as the Heaviside function. The time dependences of the dynamic stress intensity factors are obtained. Their dependence on the relative position of the cracks in the half-space is analyzed.     

International conference on role of fracture mechanics in modern technology : Fukuoka, Japan, June 2–6, 1986

An infinite elastic plane containing two straight cracks of arbitrary length and location is analyzed within the framework of elastostatics. The mathematical formulation is based on the stress solution for a single crack and leads to a system of singular integral equations that govern the crack surface displacement densities. The solution series in terms of the reciprocal of the crack centre distance is not suitable for cracks that are spaced too closely. It is shown by way of examples that the method of asymptotic solution is convenient for developing approximation expressions of the stress and displacement field with certain characteristics. The formulas for the stress intensity factors and crack opening are given for the case of a constant tensile load. Graphical results are given for the variations of the stress intensity factors with parameters depending on the relative positions of the cracks.     

Asymptotic elastic solution for two straight cracks of arbitrary length and location

An infinite elastic plane containing two straight cracks of arbitrary length and location is analyzed within the framework of elastostatics. The mathematical formulation is based on the stress solution for a single crack and leads to a system of singular integral equations that govern the crack surface displacement densities. The solution series in terms of the reciprocal of the crack centre distance is not suitable for cracks that are spaced too closely. It is shown by way of examples that the method of asymptotic solution is convenient for developing approximation expressions of the stress and displacement field with certain characteristics. The formulas for the stress intensity factors and crack opening are given for the case of a constant tensile load. Graphical results are given for the variations of the stress intensity factors with parameters depending on the relative positions of the cracks.     

Dynamic interaction of multiple arbitrarily oriented planar cracks in a piezoelectric space: A semi-analytic solution

The problem of determining the electro-elastic fields around arbitrarily oriented planar cracks in an infinite piezoelectric space is considered. The cracks which are acted upon by a transient load are either electrically impermeable or permeable. A semi-analytic method based on the theory of exponential Fourier transformation is proposed for solving the problem in the Laplace transform domain. The Laplace transforms of the jumps in the displacements and electric potential across opposite crack faces are determined by solving a system of hypersingular integral equations. Once these displacement and electric potential jumps are obtained, the displacements and electric potential and other physical quantities of interest, such as the crack tip stress and electric displacement intensity factors, can be computed with the help of a suitable algorithm for inverting Laplace transforms. The stress and electric displacement intensity factors are computed for some specific cases of the problem.     

Mixed mode axisymmetric annular cracks in a finite layer: Off angle crack initiation

The solutions of axisymmetric Volterra type climb and glide edge dislocations are obtained in a layer by means of the Hankel transforms. Utilizing the same procedure, Green’s function solution is obtained for a layer under self-equilibration normal ring traction. The distributed dislocation technique is used to construct integral equations for a system of co-axial annular cracks where the layer is under axisymmetric normal loads. These equations are solved numerically to obtain dislocation density on the cracks surfaces. The results are employed to determine stress intensity factors for annular and penny-shaped cracks and the interaction between two co-axial penny-shaped cracks is studied. Moreover, the stress intensity factors of the interacting cracks are determined such that they can be further used in conjunction with strain energy density (SED) failure criterion to obtain the possible direction of crack initiation that may not be apparent under mixed mode conditions.     

Retardation of a crack with connections between the faces using an induced thermoelastic stress field

This paper considers local temperature variations near the tip of a crack in the presence of regions in which the crack faces interact. It is assumed that these regions are adjacent to the crack tip and are comparable in size to the crack size. The problem of local temperature variations consists of delay or retardation of crack growth. For a crack with connections between the crack faces subjected to external tensile loads, an induced thermoelastic stress field, and the stresses at the connections preventing crack opening, the boundary-value problem of the equilibrium of the crack reduces to a system of nonlinear singular integrodifferential equations with a Cauchy kernel. The normal and tangential stresses at the connections are found by solving this system of equations. The stress intensity factors are calculated. The energy characteristics of cracks with tip regions are considered. The limiting equilibrium condition for cracks with tip regions is formulated using the criterion of limiting stretching of the connections.Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 1, pp. 133–143, January–February, 2005     

界面下主微裂纹间干涉的研究

涉及两相正交各向异性体界面干涉问题的研究,多裂纹问题被分解为只含单裂纹的子问题,利用位错理论和裂面应力自由条件,列出一组可数值求解位错密度函数的奇异积分方程,从耐 注得应力强度因子。     

Two mode III internal cracks located within two bonded functionally graded piezoelectric half planes respectively

This paper studies the mode III crack problem of two bonded functionally graded piezoelectric half planes which contain a crack respectively. These two cracks are located normal to the interface. All the material properties are assumed to vary along the direction of the crack line. A system of singular integral equations for electrically impermeable and permeable cracks is derived and solved numerically by using the Gauss–Chebyshev integration formula. The influence of the nonhomogeneous parameters and the dependence of the crack interactions on the stress and electric displacement intensity factors are investigated.     

On quadruple integral equations related to a certain crack problem

An exact solution of a four part mixed boundary value problem representing a three colinear crack system connected with specified crack opening displacements between the cracks is obtained. The three cracks thus become one with pressure and/or opening displacement prescribed on the crack face. From considerations of dual symmetry and a formulation based on Papkovich-Neuber harmonic functions, the boundary value problem is reduced to solving a quadruple set of integral equations. An exact solution of these equations is derived using a modified finite Hilbert transform technique. The closed form results for the stress distributions and the crack-tip stress intensity factors are presented. Limiting cases of the solution yield results which agree with well known solutions.     

Interaction of two offset interfacial cracks in bonded dissimilar media with a functionally graded interlayer: Antiplane deformation

This paper is concerned with the problem of bonded dissimilar, homogeneous media with a functionally graded interlayer, weakened by two offset interfacial cracks under antiplane deformation. Based on the Fourier integral transform method, formulation of the crack problem is reduced to a system of Cauchy-type singular integral equations. The mode III stress intensity factors are defined and evaluated in terms of the solution to the integral equations. Numerical results include the variations of stress intensity factors versus offset distance between the two cracks for various combinations of material and other geometric parameters of the bonded system, addressing the interaction of the two neighboring interfacial cracks spaced apart by the graded interlayer.     

Opening-Function Simulation of the Three-Dimensional Nonstationary Interaction of Cracks in an Elastic Body

Integral relations between three-dimensional dynamic displacements (stresses) in an infinite elastic body with arbitrarily located plane cracks and discontinuities in the displacements of the opposite crack faces are presented. The influence of opening cracks on each other is considered in the problem on crack faces loaded by pulse forces. This problem is reduced to a system of boundary integral equations of the wave-potential type in a time domain. The dynamic mode I stress intensity factors are determined for two coplanar elliptic cracks under forces in the form of the Heaviside function     

Multiple cracks in thermoelectroelastic bimaterials

The formulation for thermal stress and electric displacement in an infinite thermopiezoelectric plate with an interface and multiple cracks is presented. Using Green's function approach and the principle of superposition, a system of singular integral equations for the unknown temperature discontinuity defined on each crack face is developed and solved numerically. The formulation can then be used to calculate some fracture parameters such as the stress–electric displacement and strain energy density factor. The direction of crack growth for many cracks in thermopiezoelectric bimaterials is predicted by way of the strain energy density theory. Numerical results for stress–electric displacement factors and crack growth direction at a particular crack tip in two crack system of bimaterials are presented to illustrate the application of the proposed formulation.     

Transient response of a piezoelectric ceramic with two coplanar cracks under electromechanical impact

The transient response of two coplanar cracks in a piezoelectric ceramic under antiplane mechanical and inplane electric impacting loads is investigated in the present paper. Laplace and Fourier transforms are used to reduce the mixed boundary value problems to Cauchy-type singular integral equations in Laplace transform domain, which are solved numerically. The dynamic stress and electric displacement factors are obtained as the functions of time and geometry parameters. The present study shows that the presence of the dynamic electric field will impede or enhance the propagation of the crack in piezoelectric ceramics at different stages of the dynamic electromechanical load. Moreover, the electromechanical response is greatly affected by the ratio of the space of the cracks and the crack length.     

Isolated crack in three-dimensional piezoelectric solid. Part II: Stress intensity factors for circular crack

This is part II of the work concerned with finding the stress intensity factors for a circular crack in a solid with piezoelectric behavior. The method of solution involves reducing the problem to a system of hypersingular integral equations by application of the unit concentrated displacement discontinuity and the unit concentrated electric potential discontinuity derived in part I [1]. The near crack border elastic displacement, electric potential, stress and electric displacement are obtained. Stress and electric displacement intensity factors can be expressed in terms of the displacement and the potential discontinuity on the crack surface. Analogy is established between the boundary integral equations for arbitrary shaped cracks in a piezoelectric and elastic medium such that once the stress intensity factors in the piezoelectric medium can be determined directly from that of the elastic medium. Results for the penny-shaped crack are obtained as an example.     

压电体中孔边Ⅲ型界面裂纹的动应力强度因子

采用Green函数法研究界面上含圆孔边界径向有限长度裂纹的两半无限压电材料对SH波的散射和裂纹尖端动应力强度因子问题.首先构造出具有半圆型凹陷半空间的位移Green函数和电场Green函数,然后采用裂纹"切割"方法构造孔边裂纹,并根据契合思想和界面上的连接条件建立起求解问题的定解积分方程.最后作为算例,给出了孔边界面裂纹尖端动应力强度因子的计算结果图并进行了讨论.     

Stress intensity factors for anisotropic layered composites: Part II—isolated double periodic cracks

The formulation in Part I (Theoret. Appl. Fracture Mech. 17, 205–219 (1992)) of this work for collinear cracks in alternate layers of an anisotropic laminate is extended to a system where each of the cracked layer contains a periodic array of parallel cracks. These cracks are also collinear such that they are periodic in two mutually perpendicular directions. Finite Fourier transform is applied at discrete points for one of the space variables reducing the problem to a singular integral equation. The unknown is expressible in terms of the crack opening displacement as in Part I. Displayed graphically are the normalized stress intensity factor, effective stiffness of the laminate, and the interlayer stresses. The local stress intensity factor for the double array crack system is always less than that for a single isolated crack depending on the periodicity ratio. The interlayer stresses directly ahead of the crack are elevated, the intensity of which increases with decreasing distance between the crack tip and interface. Increase in the thickness of the adjoining layer tends to decrease the interlayer stresses nearest to the crack tip, a result that is to be expected.