Dynamic stress intensity factors around two rectangular cracks in an infinite elastic plate under impact load

A dynamic problem for two equal rectangular cracks in an infinite elastic plate is considered. The two cracks are placed perpendicular to the plane surfaces of the plate. An incoming shock tensile stress is returned by the cracks. In the Laplace transform domain, the boundary conditions at the two sides of the plate are satisfied using the Fourier transform technique. The mixed boundary conditions are reduced to dual integral equations. Crack displacement is expanded in a series of functions which are zero outside of the cracks. The unknown coefficients in the series are determined by the Schmidt method. The stress intensity factors are defined in the Laplace transform domain and these are inverted using a numerical method.     

Fracture analysis of a weak-discontinuous interface in a symmetrical functionally gradient composite strip loaded by anti-plane impact

The anti-plane impact fracture analysis was performed for a weak-discontinuous interface in a symmetrical functionally gradient composite strip. A new bi-parameter exponential function was introduced to simulate the continuous variation of material properties. Using Laplace and Fourier integral transforms, we reduced the problem to a dual integral equation and obtained asymptotic analytical solution of crack-tip stress field. Based on the numerical solution of the second kind of Fredholm integral equation transformed from the dual integral equation, the effects of the two non-homogeneity parameters on DSIF were discussed. It was indicated that the relative stiffness of the interface and the general stiffness of the whole structure are two important factors affecting the impact fracture behavior of the weak-discontinuous interface. The greater the relative stiffness of the interface is, the higher the value of the dynamic stress intensity factor will be. The greater the general stiffness of the whole structure is, the shorter the time for DSIF to arrive at the peak value and then to stabilize to the steady one. If the general stiffness of the whole structure is great enough, there will be an oscillation between the peak and steady values of DSIF.     

Dynamic stress intensity factors of multiple cracks in a functionally graded orthotropic half-plane

In the present paper dynamic stress intensity factor and strain energy density factor of multiple cracks in the functionally graded orthotropic half-plane under time-harmonic loading are investigated. By utilizing the Fourier transformation technique the stress fields are obtained for a functionally graded orthotropic half-plane containing a Volterra screw dislocation. The variations of the material properties are assumed to be exponential forms which the equilibrium has an analytical solution. The dislocation solution is utilized to formulate integral equation for the half-plane weakened by multiple smooth cracks under anti-plane deformation. The integral equations are of Cauchy singular type at the location of dislocation which are solved numerically to obtain the dislocation density on the faces of the cracks. The dislocation densities are employed to determined stress intensity factor and strain energy density factors (SEDFs) for multiple smooth cracks under anti-plane deformation. Numerical examples are provided to show the effects of material properties and the crack configuration on the dynamic stress intensity factors and SEDFs of the functionally graded orthotropic half-plane with multiple curved cracks.     

Scattering of SHPlease define abbreviation SH waves by arc-shaped interface cracks between a cylindrical magneto-electro-elastic inclusion and matrix: near fields

Summary In this paper, the scattering of SH waves by a magneto-electro-elastic cylindrical inclusion partially debonded from its surrounding magneto-electro-elastic material is investigated by using the wavefunction expansion method and a singular integral equation technique. The debonding regions are modeled as multiple arc-shaped interface cracks with non-contacting faces. The magneto-electric impermeable boundary conditions are adopted. By expressing the scattered fields as wavefunction expansions with unknown coefficients, the mixed boundary-value problem is firstly reduced to a set of simultaneous dual-series equations. Then, dislocation density functions are introduced as unknowns to transform these dual-series equations to Cauchy singular integral equations of the first type,which can be numerically solved easily. The solution is valid for arbitrary number and size of the arc-shaped interface cracks. Finally, numerical results of the dynamic stress intensity factors are presented for the cases of one debond. The effects of incident direction, crack configuration and various material parameters on the dynamic stress intensity factors are discussed. The solution of this problem is expected to have applications in the investigation of dynamic fracture properties of magneto-electro-elastic materials with cracks.The work was supported by the National Natural Science Fund of China (Project No. 19772029) and the Research Fund for Doctors of Hebei Province, China (Project No. B2001213).     

New stress intensity factor solutions of a periodic array of cracks in residual stress field

Many important applications of crack mechanics involve self-equilibrating residual or thermal stress fields. For these types of problems, the traditional fracture mechanics approach based on the superposition principle has ignored the effect of crack surface contact when the crack-tip propagates into the residual compressive region. Contact between the crack faces and the wedging action are responsible for subsequent crack-tip reopening, which often leads to a much larger mode I stress intensity factor. In this study, an analytical approach is used to study the effect of crack face contact for a period array of collinear cracks embedded in several typical residual stress fields. It is found that the nonlinear contact between crack surfaces dominates the cracking behavior in residual/thermal stress fields, which is responsible for crack coalescence.     

A study on stress intensity factors and singular stress fields of three-dimensional interface crack

By using the finite-part integral concepts and limit technique, the hypersingular integrodifferential equations of three-dimensional (3D) planar interface crack were obtained; then the dominant-part analysis of 2D hypersingular integral was further used to investigate the stress fields near the crack front theoretically, and the accurate formulae were obtained for the singular stress fields and the complex stress intensity factors. After that, a numerical method is proposed to solve the hypersingular integrodifferential equations of 3D planar interface crack, and the problem of elliptical planar crack is then considered to show the application of the method. The numerical results obtained are satisfactory. Project supported by the Foundation of Solid Mechanics Open Research Laboratory of State Education Commission at Tongji University and the National Natural Science Foundation.     

Dynamic behavior of two collinear anti-plane shear cracks in a functionally graded layer bonded to dissimilar half planes

In recent years, the functionally graded materials (FGMs) have been widely applied in extremely high temperate environment. In this paper, the dynamic behavior of two collinear cracks in FGM layer bonded to dissimilar half planes under anti-plane shear waves is studied by the Schmidt method. By using the Fourier transform technique, the present problem can be solved with a dual integral equation. These equations are solved using the Schmidt method. The present method is used to illustrate the fundamental behavior of the interacting cracks in FGMs under dynamic loading. Furthermore, the effects of the geometry of the interacting cracks, the shear stress wave velocity of the materials and the frequency of the incident wave on the Dynamic Stress Intensity Factor are investigated.     

Dynamic stress intensity factors of two 3D rectangular cracks in a transversely isotropic elastic material under a time-harmonic elastic P-wave

The dynamic stress intensity factors (DSIFs) of two 3D rectangular cracks in a transversely isotropic elastic material under an incident harmonic stress wave are investigated by generalized Almansi’s theorem and the Schmidt method in the present paper. Using 2D Fourier transform and defining the jumps of displacement components across the crack surface as the unknown functions, three pairs of dual integral equations are derived. To solve the dual integral equations, the jumps of the displacement components across the crack surfaces are expanded in a series of Jacobi polynomials. Numerical examples are provided to show the effects of the geometric shape of the rectangular crack, the characteristics of the harmonic wave and the distance between two rectangular cracks on the DSIFs of the transversely isotropic elastic material.     

Interaction of elastic waves with a periodic array of collinear inplane cracks

A boundary integral equation method is applied to the study of the interaction of plane elastic waves with a periodic array of collinear inplane cracks. Numerical results are presented for the dynamic stress intensity factors. The effects of the wave type, wave frequency, wave incidence angle, and crack spacing on the dynamic stress intensity factors are analyzed in detail. The project supported by the Committee of Science and Technology of Shanghai and Tongji University     

Analysis of mode Ⅲ crack perpendicular to the interface between two dissimilar strips

The present work is concerned with the problem of mode Ⅲ crack perpendicular to the interface of a bi-strip composite. One of these strips is made of a functionally graded material and the other of an isotropic material, which contains an edge crack perpendicular to and terminating at the interface. Fourier transforms and asymptotic analysis are employed to reduce the problem to a singular integral equation which is numerically solved using Gauss-Chebyshev quadrature formulae. Furthermore, a parametric study is carried out to investigate the effects of elastic and geometric characteristics of the composite on the values of stress intensity factor.     

3D dynamic stress intensity factors at three rectangular cracks in an infinite elastic medium subjected to a time-harmonic stress wave

Summary Dynamic stresses around three coplanar cracks in an infinite elastic medium are determined in the paper. Two of the cracks are equal, rectangular and symmetrically situated on either side of the centrally located rectangular crack. Time-harmonic normal traction acts on each surface of the three cracks. To solve the problem, two kind of solutions are superposed: one is a solution for a rectangular crack in an infinite elastic medium, and the other one is that for two rectangular cracks in an infinite elastic medium. The unknown coefficients in the combined solution are determined by applying the boundary conditions at the surfaces of the cracks. Finally, stress intensity factors are calculated numerically for several crack configurations. Received 14 July 1998; accepted for publication 2 December 1998     

Dynamic stress intensity factors of two collinear cracks in orthotropic medium subjected to time-harmonic disturbance

Dynamic stresses are obtained for an infinite orthotropic medium weakened by two collinear cracks. Time-harmonic elastic waves are interrupted at normal incidence by the line cracks. Fourier transform is applied reducing the problem to solving a pair of dual integral equations. Solution method involves expanding the crack surface displacement in a series of functions that vanish along a collinear line outside the cracks. The unknown coefficients in the series are evaluated by using the Schmidt method. Dynamic stress intensity factors are computed and displayed graphically for an orthotropic medium that corresponds to the elastic properties of boron-epoxy composite.     

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.     

Dynamic behavior of two collinear interface cracks between two dissimilar functionally graded piezoelectric/piezomagnetic material strips

The dynamic interaction of two collinear interface cracks between two dissimilar functionally graded piezoelectric/piezomagnetic material strips subjected to the anti-plane shear harmonic stress waves was investigated. By using the Fourier transform, the problem can be solved with the help of a pair of triple integral equations in which the unknown variable is jump of displacement across the crack surfaces. These equations are solved using the Schmidt method. Numerical examples are provided to show the effect of the functionally graded parameter, the circular frequency of the incident waves and the thickness of the strip upon stress, electric displacement and magnetic flux intensity factors of cracks.     

Stress intensity factors for curved and kinked cracks in plane extension

A singular integral equation containing the crack opening displacement (COD) is developed for solving plane elasticity problems. The crack may contain any number of kinks at different intervals and orientations, such as a saw-tooth shape. Cracks in the form of a sine wave can also be treated. The crack tip stress intensity factors are evaluated for a variety of crack shapes and the results are displayed graphically. The distance between the crack tips is found to be a dominant factor on the crack tip stress intensity while the angle between the tangent to the crack tip and load direction determines the proportion of Mode I and II stress intensity factors.     

2D dynamic analysis of cracks and interface cracks in piezoelectric composites using the SBFEM

The dynamic stress and electric displacement intensity factors of impermeable cracks in homogeneous piezoelectric materials and interface cracks in piezoelectric bimaterials are evaluated by extending the scaled boundary finite element method (SBFEM). In this method, a piezoelectric plate is divided into polygons. Each polygon is treated as a scaled boundary finite element subdomain. Only the boundaries of the subdomains need to be discretized with line elements. The dynamic properties of a subdomain are represented by the high order stiffness and mass matrices obtained from a continued fraction solution, which is able to represent the high frequency response with only 3–4 terms per wavelength. The semi-analytical solutions model singular stress and electric displacement fields in the vicinity of crack tips accurately and efficiently. The dynamic stress and electric displacement intensity factors are evaluated directly from the scaled boundary finite element solutions. No asymptotic solution, local mesh refinement or other special treatments around a crack tip are required. Numerical examples are presented to verify the proposed technique with the analytical solutions and the results from the literature. The present results highlight the accuracy, simplicity and efficiency of the proposed technique.     

Finite-part integral and boundary element method for interaction between two parallel planar cracks in three-dimensional finite bodies

By using the Somigliana representation and the concepts of finite-part integrals, a set of hypersingular integral equations of the interaction between two parallel planar cracks in a three-dimensional finite body subjected to arbitrary loads is derived, and then its numerical method is proposed by the finite-part integral method combined with the boundary element method. According to the analytic theory of hypersingular integral equations, the square root models of displacement discontinuities in the elements near the crack front are applied, and thus the computational precision is raised. Based on this, the stress intensity factors can be directly calculated. Finally, the stress intensity factors of several typical interaction problems are calculated.     

A crack perpendicular to and terminating at a bimaterial interface

Using dislocation simulation approach, the basic equation for a finite crack perpendicular to and terminating at a bimaterial interface is formulated. A novel expansion method is proposed for solving the problem. The complete solution to the problem, including the explicit formulae for theT stresses ahead of the crack tip and the stress intensity factors are presented. The stress field characteristics are analysed in detail. It is found that normal stresses {ie27-1} and {ie27-2} ahead of the crack tip, are characterised byQ fields if the crack is within a stiff material and the parameters |p _T | and |q _T | are very small, whereQ is a generalised stress intensity factor for a crack normal to and terminating at the interface. If the crack is within a weak material, the normal stresses {ie27-3} and {ie27-4} are dominated by theQ field plusT stress. This work was supported by the Swedish Research Council for Engineering Sciences.     

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     

基体裂纹与界面刚性线的弹性干涉

研究两种材料界面上的刚性线与其它任意位置处直线裂纹弹性干涉的反平面问题。基于界面上刚性线与任意位置处螺型位错干涉的基本解,运用连续位错密度模型法将问题转化为奇异积分方程。用半开型积分法求解奇异积分方程,得到位错密度函数的离散值,计算裂纹尖端处的应力强度因子。算例说明该方法可用于工程实际问题。