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     

Elastodynamic analysis of a crack at an arbitrary angle to the graded interfacial zone in bonded half-planes under antiplane shear impact

A solution is provided for the elastodynamic problem of a crack at an arbitrary angle to the graded interfacial zone in bonded media under the action of antiplane shear impact. The interfacial zone is modeled by a nonhomogeneous interlayer with the spatially varying shear modulus and mass density in terms of power functions between the two dissimilar, homogeneous half-planes. Based on the use of Laplace and Fourier integral transforms and the coordinate transformations of basic field variables, formulation of the transient crack problem is reduced to solving a Cauchy-type singular integral equation in the Laplace transform domain. The crack-tip response in the physical domain is recovered via the inverse Laplace transform and the values of dynamic mode III stress intensity factors are obtained as a function of time. A comprehensive parametric study is then presented of the effects of crack obliquity on the overshoot behavior of the transient crack-tip response, by plotting the peak values of the dynamic stress intensity factors versus the crack orientation angle for various material and geometric combinations of the bonded system.     

Transient crack growth under creep conditions

Transient crack growth in an elastic/power-law creeping material is investigated under antiplane shear loading and small-scale-yielding conditions. At time t = 0 the solid is suddenly loaded far from the crack by tractions that correspond to the elastic crack-tip stress distribution. At that time the crack begins to propagate at a constant velocity. The stress fields evolve in a complex manner as the crack propagates due to the competing effects of stress relaxation due to constrained creep and stress elevation due to the instantaneous elastic material response to crack growth. From detailed finite element calculations it is shown that these fields can be approximated by a simple matching of three asymptotic singular crack-tip solutions. A characteristic stress, distance and time are defined for this problem which provide a normalization that accounts for any crack velocity, loading and all material properties for a given creep exponent n. Results are presented for crack-tip stresses, strains, crack opening displacements and creep zones.     

An embedded crack in a functionally graded orthotropic coating bonded to a homogeneous substrate under a frictional Hertzian contact

In this paper, we consider the elasto-static problem of an embedded crack in a graded orthotropic coating bonded to a homogeneous substrate subject to statically applied normal and tangential surface loading. The crack direction is parallel to the free surface. The coating is graded in the thickness direction and is orthogonal to the crack direction. This coating is modelled as a non-homogeneous medium with an orthotropic stress–strain law. The equivalent crack surface stresses are first obtained and substituted in the plane elasticity equations. Using integral transforms, the governing equations are converted into singular integral equations which are solved numerically to yield the displacement field as well as the crack-tip stress intensity factors. This study presents a complete theoretical formulation for the problem in the static case. A numerical predictive capability for solving the singular integral equations and computing the crack-tip stress intensity factors is proposed. Since the loading is compressive, a previously developed crack-closure algorithm is applied to avoid interpenetration of the crack faces. The main objective of the paper is to investigate the effects of the material orthotropy and non-homogeneity of the graded coating on the crack-tip stress intensity factors, with and without using the crack-closure algorithm, for the purpose of gaining better understanding on the behavior and design of graded coatings.     

Penny-shaped cracks

The Somigliana formula is used to reduce an arbitrary elastic crack problem to a system of three integral equations for the components of displacement discontinuity. For the case of a penny shaped crack situated in an infinite isotropic medium with the crack faces subjected to arbitrary tractions, the integral equations are solved explicitly. In particular integral formulae are obtained for the stresses on the plane of the crack beyond the crack-tip, and hence for the stress intensity factors. The special case of uni-directional shear traction on the crack is examined.     

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

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

Dynamic crack-tip fields in rate-sensitive solids

The asymptotic stress and strain fields near the tip of a crack which propagates dynamically in a rate-sensitive solid are obtained under anti-plane shear and plane strain conditions. The problem is formulated within the context of a small-strain theory for a solid whose mechanical behavior under high strain rates is described by an elastic-viscoplastic constitutive relation. It is shown that, if the stresses are singular at the crack-tip, the viscoplastic relation is equivalent asymptotically to an elastic-non-linear viscous relation. Furthermore, for a certain range of the material parameter which characterizes the rate-sensitivity of the material, the elastic strain-rates near the propagating crack tip are shown to have the same asymptotic radial dependence near the propagating crack-tip as the inelastic strain-rates. This determines the order of the stress singularity uniquely. The governing equations for anti-plane shear and plane strain are then derived. The numerical results for the stress and strain fields are presented for anti-plane shear and plane strain. For the present model, the results suggest that under small-scale yielding conditions, there exists a minimum velocity for stable steady crack propagation. The implication that a terminal velocity for a running crack may exist is also discussed.     

Intergranular and transgranular crack growth at triple junction boundaries in ordered intermetallics

Coincident site-lattice (CSL) and random grain boundaries (GBs) effects on intergranular and transgranular crack propagation paths in ordered intermetallics that are subjected to high rates of strain are investigated. A three dimensional dislocation density based multiple slip crystalline formulation and computational scheme are used for a detailed understanding and accurate characterization of interrelated deformation and failure mechanisms that can occur due to the generation, trapping, interaction, and annihilation of mobile and immobile dislocation densities that are generally associated with finite strain high strain-rate plasticity in L1₂ ordered intermetallics. Results from this study indicate that intergranular crack growth is along the GBs, normal to the stress-axis, and is due to the dominance of normal stresses in the crack-tip region. Transgranular crack growth is along slip-planes, and is due to the dominance of shear stresses in the crack-tip region.     

A numerical approach for a crack problem by Gauss–Chebyshev quadrature

In this paper, a problem of a crack in an orthotropic strip is studied under plane strain conditions. It is assumed that normal displacements and shear stresses do not act on neither of the boundaries of the strip. Cauchy-type singular integral equation for the crack problem is derived by using the theory of plane elasticity and the Fourier transformation technique. A quadrature collocation approach is adopted for the numerical solutions of the singular integral equation. The effect of relative thickness and mechanical properties of strip on Mode I stress intensity factors (SIFs) are examined under different loading conditions. Some sample results are given for SIFs; also, material orthotropy and geometrical effects are discussed in detail.     

Orthotropic semi-infinite medium with a crack under thermal shock

A cracked orthotropic semi-infinite plate under thermal shock is investigated. The thermal stresses are generated due to sudden cooling of the boundary by ramp function temperature change. The superposition technique is used to solve the problem. The crack problem is formulated by applying the thermal stresses obtained from the uncracked plate with opposite sign to be the only external loads on the crack surfaces as the crack surface tractions. The Fourier transform technique is used to solve the problem leading to a singular equation of the Cauchy type. The singular integral equation is solved numerically using the expansion method. The influence of the material orthotropy on the stress intensity factors is shown by comparing the results obtained for different orthotropic materials and isotropic materials in the case of plane stress. The numerical results of the stress intensity factors are demonstrated as a function of time, crack length, location of the crack and the duration of the cooling rate.     

黏弹性体界面裂纹的冲击响应

研究两半无限大黏弹性体界面Griffith裂纹在反平面剪切突出载荷下,裂纹尖端动应力强度因子的时间响应,首先,运用积分变换方法将黏弹性混合黑社会问题化成变换域上的对偶积分方程,通过引入裂纹位错密度函数进一步化成Cauchy型奇异积分方程,运用分片连续函数法数值求解奇异积分方程,得到变换域内的动应力强度因子,再用Laplace积分变换数值反演方法,将变换域的解反演到时间域内,最终求得动应力强度因子的时间响应,并对黏弹性参数的影响进行分析。     

ANTIPLANE CIRCULAR INCLUSION WITH A CURVED CRACK CROSSING THE BOUNDARY

The weakly singular integral equation used to solve the problem of the curved crack crossing the boundary of the antiplane circular inclusion is presented. Using the principal part analysis method of the Cauchy type integral equation, the singular stress index at the intersection and the singular stress of angular regions near the intersection are obtained. By using the singular stress obtained, the stress intensity factor at the intersection is, defined. After the numerical solution of the integral equation, the stress intensity factors at the end points of the crack and intersection are obtainable. The research is supported by National Natural Science Foundation of China (No. 59879012) and is the project of Chinese Foundation of State Education Commission (No. 98024832).     

A Long Crack Penetrating a Circular Inhomogeneity

This paper presents the effects of elastic mismatch and crack-tip position on the stress intensity factors of a long crack penetrating a circular inhomogeneity. The analysis relies on closed-form solutions, derived using complex variable techniques, for the stresses and displacements produced by dislocations positioned inside and outside the inhomogeneity. Dislocation distributions are introduced to express the traction boundary condition along the crack surfaces as a system of singular integral equations, whose solution is obtained through a numerical procedure. It is shown that if the elastic mismatch is interpreted correctly, then the stress intensity factors of this micromechanical model are very good approximations to those computed using a Monte Carlo finite element model of a long crack in a polycrystalline plate with compliant grain boundaries.     

ANALYSIS OF A CRACK IN A FUNCTIONALLY GRADED STRIP WITH A POWER FORM SHEAR MODULUS

The plane strain problem of a crack in a functionally graded strip with a power form shear modulus is studied. The governing equation in terms of Airy＇s stress function is solved exactly by means of Fourier transform. The mixed boundary problem is then reduced to a system of singular integral equations and is solved numerically to obtain the stress intensity factor at crack-tip. The maximum circumferential stress criterion and the strain energy density criterion are both employed to predict the direction of crack initiation. Numerical examples are given to show the influence of the material gradation models and the crack sizes on the mode-I and mode-II stress intensity factors. The dependence of the critical kink-angle on the crack size is examined and it is found that the crack kink-angle decreases with the increase of the normalized crack length, indicating that a longer crack tends to follow the original crack-line while it is much easier for a shorter crack to deviate from the original crack-line.     

与两相材料界面接触的裂纹对SH波的散射

利用积分变换方法得出了两相材料中作用简谐集中力时的格林函数．根据所得的格林函数并利用Betti-Rayleigh互易定理得出了与界面接触裂纹的散射波场．裂纹的散射波场可分解为两部分，一部分为奇异的散射场，另一部分为有界的散射场．利用分解后的散射场，可得裂纹在SH波作用下的超奇异积分方程．根据裂纹散射场的奇异部分和Cauchy型奇异积分的性质得出了裂纹和界面接触点处的奇性应力指数和接触点角形域内的奇性应力．利用所得的奇性应力定义了裂纹和界面接触点处的动应力强度因子．对所得超奇异积分方程的数值求解可得裂纹端点和接解点处的应力强度因子。     

On the load transfer from elastic inclusions to an infinite elastic orthotropic plane with cuts

The present paper deals with the problem of load transfer from elastic inclusions to an infinite elastic orthotropic plane with cuts located on one of the principal orthotropy directions. The constitutive system of equations of this problem is derived under the assumption that the inclusions are in a uniaxial stress state. The obtained system consists of a singular integro-differential equation and a singular integral equation for the jumps of the tangential stresses acting on the inclusion shores and for the derivative of the the cut opening function. The behavior of solutions of the system of constitutive equations at the endpoints of the inclusions and cuts is studied, and the solution of this system is constructed by the numerical-analytic discrete singularity method.     

Analysis of rotating solids with cracks by the boundary element method

The use of the boundary element method for the solution of linear elastic fracture mechanics problems, without body forces, is quite extensive since the method is intrinsically well suited to the analysis of high stress gradients associated with crack problems. The crack-tip stresses for rotating bodies are similar to the stresses for stationary bodies and therefore all the advantages of the boundary element procedure can be encompassed in the extension of the method to the solution of rotating bodies with cracks. In the present analysis, the additional volume integral that arises from the treatment of inertial body forces is eliminated by using the well-known particular integral procedure. The matrix ill-conditioning that results from the need to model co-planar crack surfaces of non-symmetrical cracks is avoided by using the multi-region approach. The accuracy of the numerical solutions is improved by utilizing quarter-point elements with traction singular enhancement at the crack-tip. The procedure is then applied to the solution of arbitrary cracks in two- and three-dimensional rotating bodies.     

Spread of plasticity from stacked stress concentrations

A large elastic solid containing an infinite sequence of slitlike relaxed cracks with a constant distance of vertical separation is considered. The solid is deforming under plane strain shear conditions (mode II). The plastic relaxation around each of the cracks is represented by a suitable distribution of edge dislocations coplanar with the crack itself, the distribution being determined from a singular integral equation. This equation is solved numerically using an expansion of the non-singular part of the kernel in a series of Chebyshev polynomials. Solutions are obtained for the extent of spread of plasticity from each crack and for the associated dislocation distribution as a function of the crack spacing and the applied load. The results are applied to a brief discussion of the fracture process at stress concentrations using the crack opening displacement criterion.     

Eccentric crack in a piezoelectric strip bonded to half planes

We consider the problem of determining the singular stresses and electric fields in a piezoelectric ceramic strip containing an eccentric Griffith crack off the centre line bonded to two elastic half planes under anti-plane shear loading using the continuous crack-face condition. Fourier transforms are used to reduce the problem to the solution of two pairs of dual integral equations, which are then expressed to a Fredholm integral equation of the second kind. Numerical values on the stress intensity factor and energy release rate are obtained.     

Griffith crack moving in a piezoelectric strip

Summary  The dynamic problem of an impermeable crack of constant length 2a propagating along a piezoelectric ceramic strip is considered under the action of uniform anti-plane shear stress and uniform electric field. The integral transform technique is employed to reduce the mixed-boundary-value problem to a singular integral equation. For the case of a crack moving in the mid-plane, explicit analytic expressions for the electroelastic field and the field intensity factors are obtained, while for an eccentric crack moving along a piezoelectric strip, numerical results are determined via the Lobatto–Chebyshev collocation method for solving a resulting singular integral equation. The results reveal that the electric-displacement intensity factor is independent of the crack velocity, while other field intensity factors depend on the crack velocity when referred to the moving coordinate system. If the crack velocity vanishes, the present results reduce to those for a stationary crack in a piezoelectric strip. In contrast to the results for a stationary crack, applied stress gives rise to a singular electric field and applied electric field results in a singular stress for a moving crack in a piezoelectric strip. Received 14 August 2001; accepted for publication 24 September 2002 The author is indebted to the AAM Reviewers for their helpful suggestions for improving this paper. The work was supported by the National Natural Science Foundation of China under Grant 70272043.