双轴载荷作用下源于椭圆孔的分支裂纹的一种边界元分析

利用一种边界元方法来研究双轴载荷作用下无限大板中源于椭圆孔的分支裂纹．该边界元方法由Crouch与Starfied建立的常位移不连续单元和笔者提出的裂尖位移不连续单元构成．在该边界元方法的实施过程中,左、右裂尖位移不连续单元分别置于裂纹的左、右裂尖处,而常位移不连续单元则分布于除了裂尖位移不连续单元占据的位置之外的整个裂纹面及其它边界,文中算例说明本数值方法对计算平面弹性裂纹的应力强度因子是非常有效的。该文对双轴载荷作用下无限大板中源于椭圆孔的分支裂纹的数值结果进一步证实本数值方法对计算复杂裂纹的应力强度因子的有效性,同时该数值结果可以揭示双轴载荷及裂纹体几何对应力强度因子的影响。     

An internal crack problem for an infinite elastic layer

The elastostatic plane problem of an infinite elastic layer with an internal crack is considered. The elastic layer is subjected to two different loadings, (a) the elastic layer is loaded by a symmetric transverse pair of compressive concentrated forces P/2, (b) it is loaded by a symmetric transverse pair of tensile concentrated forces P/2. The crack is opened by an uniform internal pressure p ₀ along its surface and located halfway between and parallel to the surfaces of the elastic layer. It is assumed that the effect of the gravity force is neglected. Using an appropriate integral transform technique, the mixed boundary value problem is reduced to a singular integral equation. The singular integral equation is solved numerically by making use of an appropriate Gauss–Chebyshev integration formula and the stress-intensity factors and the crack opening displacements are determined according to two different loading cases for various dimensionless quantities.     

弹性半平面中边界裂纹研究的新方法

讨论了载荷作用在裂纹面上的弹性半平面边界裂纹问题.研究以线弹性断裂力学为基础,采用复变函数方法以及Riemann-Hilbert(R-H)边值问题的一般理论,将问题分拆为含有限裂纹的全平面问题与无裂纹的半平面问题的叠加,计算得到裂纹尖端的应力强度因子.与文献结果比较,该方法具有精度高的优点.     

The stress intensity factors near the tips of two collinear cracks in composite under in-plane shear

In this paper, the stress-intensity factors for two collinear cracks in a composite bonded by an isotropic and an anisotropic half-plane were calculated. The cracks are paralell to the interface, and the crack surfaces are loaded by uniform shear stresses. By using Fourier transform, the mixed boundary value problem is reduced to a set of singular integral equations. For solving the integral equations, the crack surface displacements are expanded in triangular series and the unknown coefficients in the series are determined by the Schmidt method. The stress intensity factors for the cracks in the boron-fibre plastics and aluminium joined composite and in carbon-fibre reinforced plastics were calculated numerically.     

Equilibrium of an internal transverse crack in a semiinfinite elastic body with thin coating

Static elasticity problems for a half-plane and a strip weakened by a rectilinear transverse crack are studied. In each case, the upper boundary of the body is reinforced by a flexible patch. Various versions of conditions on the lower boundary are considered in the case of the strip. The crack is maintained in the open state by distributed normal forces. The method of generalized integral transforms reduces solving the problem for the equations of equilibriumto solving a singular integral equation of the first kind with the Cauchy kernel with respect to the derivative of the crack opening function. The solutions of the integral equation are constructed by the small parameter and collocation methods for various combinations of the geometric and physical parameters of the problem, and the structure of the solutions is analyzed. The values of the stress intensity factor (SIF) near the crack vertex are obtained.     

Determining the size and location of transverse cracks in beams

The spectral element method, which is very suitable for solving inverse dynamic problems, is combined with a stochastic genetic algorithm to give a scheme that can locate and size cracks in structural components. The mechanical model is based on an approach that separates the global structural dynamics from the local crack-tip zone dominated by singular stresses. The global model, consisting of connected waveguides, describes the structural dynamics using spectral elements. The local model, describing the crack-tip behavior, is based on the relation between the stress-intensity factor and the stored strain energy representing the crack region. The results are demonstrated with experimental data from an aluminum beam with a transverse crack.     

Multiple interacting cracks in an orthotropic layer

The stress fields in an orthotropic layer containing climb and glide edge dislocations are obtained by means of the complex Fourier transform. Stress analysis in the intact layer under in-plane point loads is also carried out. These solutions are employed to derive integral equations for the layers weakened by several interacting cracks subject to in-plane deformation. The integral equations are of Cauchy singular type. These equations are solved numerically for the density of dislocations on a crack surface. The dislocation densities are utilized to derive stress intensity factor for cracks. Several examples are solved and the interaction between the two cracks is investigated.     

Calculation of the scattering of elastic waves from a penny-shaped crack by the method of optimal truncation

The method of optimal truncation (MOOT) [1, 2], a least-squares boundary-residual method [3] for solving scattering problems, is applied to the plane circular crack. An equatorially cloven spherical inclusion is used to model the crack. Numerical advantages of this model are discussed and demonstrated. Results are given for cross sections for longitudinal waves incident on the crack at arbitrary angles. Both clear cracks and fluid-filled cracks are considered. A refinement of the method which would allow accurate calculation of dynamic stress-intensity factors is developed.     

Stress analysis near the tips of a transverse crack in an elastic semi-strip

The plane elastic problem for a semi-strip with a transverse crack is investigated. The initial problem is reduced to a one-dimensional continuous problem by use of an integral transformation method with a generalized scheme. The one-dimensional problem is first formulated as a vector boundary problem, and then reduced to a system of three singular integral equations(SIEs). The system is solved by use of an orthogonal polynomial method and a special generalized method. The contribution of this work is the consideration of kernel fixed singularities in solving the system. The crack length and its location relative to the semi-strip's lateral sides are investigated to simplify the problem's statement. This simplification reduces the initial problem to a system of two SIEs.     

Some transient coupled thermoelastic crack problems

This paper is concerned with determining the elastodynamic response of a plane strain medium containing a central crack deformed by the action of suddenly applied thermal and/or mechanical disturbances when the assumptions of the general theory of coupled thermoelasticity are assumed. Integral transform solution is employed to reduce the governing equations into integral equations of Fredholm type. A numerical inversion technique is used to compute the dynamic stress-intensity factors when the faces of the crack are subjected to constant heat flux and/or mechanical loading. Attention is focused on the overshoot in the stress-intensity factor and its time interval for non-stationary temperature fields, and to what degree it is influenced by the mutual dependence of the temperature and displacement fields inherent in the coupled theory of thermoelasticity.     

Thermoelastic stress intensity factors for a cracked layer between two different anisotropic solids

Steady-state anisotropic thermoelasticity equations are used to obtain the stress intensity factors for a cracked layer sandwiched between two different anisotropic elastic solids. The anisotropy is assumed to arise from discrete fibers whose orientation could alter with reference to the crack edges. A generalized plane deformation prevails in the dissimilar media domain with a line of discontinuity disturbing a uniform heat flow. The flexibility/stiffness matrix approach is used such that the crack problem reduces to solving two sets of singular integral equations. Numerical values of the crack tip stress-intensity factors are obtained for various crack size, crack location, crack surface insulation, fiber volume fraction and orientation angles. The results are displayed graphically.     

Determination of stress intensity factors in half-plane containing several moving cracks

The dynamic stress intensity factors in a half-plane weakened by several finite moving cracks are investigated by employing the Fourier complex transformation. Stress analysis is performed in a half-plane containing a single dislocation and without dislocation. An exact solution in a closed form to the stress fields and displacement is ob- tained. The Galilean transformation is used to transform between coordinates connected to the cracks. The stress components are of the Cauchy singular kind at the location of dislocation and the point of application of the the influence of crack length and crack running force. Numerical examples demonstrate velocity on the stress intensity factor.     

曲线裂纹和反平面圆形夹杂相交问题

建立了和反平面圆夹杂界面相交的曲线裂纹的弱奇异积分方程,利用Cauchy型奇异积分方程主部分析方法研究了穿过反平面圆夹杂界面的曲线裂纹在交点处的奇性应力指数以及交点处角形域内的奇性应力,并根据奇性应力定义了交点处的应力强度因子。通过对弱奇异积分方程的数值求解,可得裂纹端点和交点处的应力强度因子。     

Analysis of a Griffith crack at the interface of two piezoelectric materials under anti-plane loading

The main objective of this work is the contribution to the study of the piezoelectric structures which contain preexisting defect (crack). For that, we consider a Griffith crack located at the interface of two piezoelectric materials in a semi-infinite plane structure. The structure is subjected to an anti-plane shearing combined with an in-plane electric displacement. Using integral Fourier transforms, the equations of piezoelectricity are converted analytically to a system of singular integral equations. The singular integral equations are further reduced to a system of algebraic equations and solved numerically by using Chebyshev polynomials. The stress intensity factor and the electric displacement intensity factor are calculated and used for the determination of the energy release rate which will be taken as fracture criterion. At the end, numerical results are presented for various parameters of the problem; they are also presented for an infinite plane structure.     

Variation of the stress intensity factor along the front of a 3-D rectangular crack subjected to mixed-mode load

Summary  The singular integral equation method is applied to the calculation of the stress intensity factor at the front of a rectangular crack subjected to mixed-mode load. The stress field induced by a body force doublet is used as a fundamental solution. The problem is formulated as a system of integral equations with r ⁻³-singularities. In solving the integral equations, unknown functions of body-force densities are approximated by the product of polynomial and fundamental densities. The fundamental densities are chosen to express two-dimensional cracks in an infinite body for the limiting cases of the aspect ratio of the rectangle. The present method yields rapidly converging numerical results and satisfies boundary conditions all over the crack boundary. A smooth distribution of the stress intensity factor along the crack front is presented for various crack shapes and different Poisson's ratio. Received 5 March 2002; accepted for publication 2 July 2002     

Scattering of monochromatic elastic waves on a planar crack of arbitrary shape

Scattering of monochromatic elastic waves on an isolated planar crack of arbitrary shape is considered. The 2D-integral equation for the crack opening vector is discretized by Gaussian approximating functions. For such functions, the elements of the matrix of the discretized problem have forms of standard one-dimensional integrals that can be tabulated. For regular grids of approximating nodes, the matrix of the discretized problem has the Toeplitz structure, and the corresponding matrix–vector products can be calculated by the fast Fourier transform technique. The latter strongly accelerates the process of iterative solution of the discretized problem. Examples of calculations of crack opening vectors, dynamic stress-intensity factors, and differential cross-sections of circular (penny-shaped) and non-circular cracks for various incident wave fields are presented. For a penny-shaped crack and longitudinal incident waves normal to the crack plane, an efficient semi-analytical method of the solution of the scattering problem is developed. The results of both methods are compared in a wide frequency region of the incident field.     

Transient response of a piezoelectric ceramic strip with an eccentric crack under electromechanical impacts

The transient response of a piezoelectric strip with an eccentric crack normal to the strip boundaries under applied electromechanical impacts is considered. By using the Laplace transform, the mixed initial-boundary-value problem is reduced to triple series equations, then to a singular integral equation of the first kind by introducing an auxiliary function. The Lobatto–Chebyshev collocation technique is adopted to solve numerically the resulting singular integral equation. Dynamic field intensity factors and energy release rate are obtained for both a permeable crack and an impermeable crack. The effects of the crack position and the material properties on the dynamic stress intensity factor are examined and numerical results are presented graphically.     

The scattering of general SH plane wave by interface crack between two dissimilar viscoelastic bodies

The scattering of general SH plane wave by an interface crack between two dissimilar viscoelastic bodies is studied and the dynamic stress intensity factor at the crack-tip is computed. The scattering problem can be decomposed into two problems: one is the reflection and refraction problem of general SH plane waves at perfect interface (with no crack); another is the scattering problem due to the existence of crack. For the first problem, the viscoelastic wave equation, displacement and stress continuity conditions across the interface are used to obtain the shear stress distribution at the interface. For the second problem, the integral transformation method is used to reduce the scattering problem into dual integral equations. Then, the dual integral equations are transformed into the Cauchy singular integral equation of first kind by introduction of the crack dislocation density function. Finally, the singular integral equation is solved by Kurtz's piecewise continuous function method. As a consequence, the crack opening displacement and dynamic stress intensity factor are obtained. At the end of the paper, a numerical example is given. The effects of incident angle, incident frequency and viscoelastic material parameters are analyzed. It is found that there is a frequency region for viscoelastic material within which the viscoelastic effects cannot be ignored. This work was supported by the National Natural Science Foundation of China (No.19772064) and by the project of CAS KJ 951-1-20     

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.     

Stress analysis for an infinite strip weakned by periodic cracks

Stress analysis for an infinite stripcracks were assumed in a horizontal position,weakened by periodic cracks is studied. The and the strip was applied by tension “p“ in y-direction. The boundary value problem can be reduced into a complex mixed one. It is found that the EEVM ( eigenfunction expansion variational method) is efficient to solve the problem. The stress intensity factor at the crack tip and the T-stress were evaluated. From the deformation response under tension the cracked strip can be equivalent to an orthotropic strip without cracks. The elastic properties in the equivalent orthotropic strip were also investigated. Finally, numerical examples and results were given.