Dynamic stress intensity factor due to concentrated loads on a propagating semi-infinite crack in orthotropic materials

The elastodynamic response of an infinite orthotropic material with a semi-infinite crack propagating at constant speed under the action of concentrated loads on the crack faces is examined. Solution for the stress intensity factor history around the crack tip is found for the loading modes I and II. Laplace and Fourier transforms along with the Wiener-Hopf technique are employed to solve the equations of motion. The asymptotic expression for the stress near the crack tip is analyzed which lead to a closed-form solution of the dynamic stress intensity factor. It is found that the stress intensity factor for the propagating crack is proportional to the stress intensity factor for a stationary crack by a factor similar to the universal function k(v) from the isotropic case. Results are presented for orthotropic materials as well as for the isotropic case.     

Dynamic stress field for torsional impact of a penny-shaped crack in a transversely isotropic functional graded strip

The torsional impact response of a penny-shaped crack in a transversely isotropic strip is considered. The shear moduli are assumed to be functionally graded such that the mathematics is tractable. Laplace and Hankel transforms are used to reduce the problem to solving a Fredholm integral equation. The crack tip stress field is obtained by considering the asymptotic behavior of Bessel function. Investigated are the effects of material nonhomogeneity and orthotropy and strip’s highness on the dynamic stress intensity factor. The peak of the dynamic stress intensity factor can be suppressed by increasing the shear moduli’s gradient and/or increasing the shear modulus in a direction perpendicular to the crack surface. The dynamic behavior varies little with the increasing of the strip’s highness.     

The anti-plane shear problem of debonding of two solids made of power law hardening materials

Debonding of two different solids made of power law hardening materials is studied for the case of anti-plane shear loading mode by using an interface crack model. The stresses and the stress intensity factor at the interface crack are determined analytically. Using these analytical results, the constitutive equations by Hencky–Ilyushin and the general equation of energy in the neighborhood of the crack tip, the adhesion energy for the loading mode under consideration is found analytically. It can be observed that for the particular case of two linearly elastic materials and a homogeneous linearly elastic material the solution found here is in excellent agreement with the solutions found in the literature.     

Boundary element analysis of three-dimensional crack problems in two joined transversely isotropic solids

The present paper presents a boundary element analysis of penny-shaped crack problems in two joined transversely isotropic solids. The boundary element analysis is carried out by incorporating the fundamental singular solution for a concentrated point load in a transversely isotropic bi-material solid of infinite space into the conventional displacement boundary integral equations. The conventional multi-region method is used to analyze the crack problems. The traction-singular elements are employed to capture the singularity around the crack front. The values of the stress intensity factors are obtained by using crack opening displacements. The numerical scheme results are verified with the closed-form solutions available in the literature for a penny-shaped crack parallel to the plane of the isotropy of a homogeneous and transversely isotropic solid of infinite extent. The new problem of a penny-shaped crack perpendicular to the interface of a transversely isotropic bi-material solid is then examined in detail. The crack surfaces are subject to the three normal tractions and the uniform shear traction. The associated stress intensity factor values are obtained and analyzed. The present results can be used for the prediction of the stability of composite structures and the hydraulic fracturing in deep rock strata and reservoir engineering.     

Thermoelastic and infrared-thermography methods for surface strains in cracked orthotropic composite materials

Two quantitative thermoelastic strain analysis (TSA) experimental methods are proposed to determine the surface strain fields in mechanically loaded orthotropic materials using the spatial distribution of temperature gradient measured from the surface. Cyclic loadings are applied to orthotropic composite specimens to achieve adiabatic conditions. The small change in surface temperatures that resulted from the change in the elastic strain energy is measured using a high sensitivity infrared (IR) camera that is synchronized with the applied loading. The first method is applied for layered orthotropic composites with a coat layer made of isotropic or in-plane transversely isotropic material. In this case, one material parameter (pre-calibrated from the surface) is required to map the strain invariant to the temperature gradients. The proposed method can be used together with Lekhnitskii’s elasticity solution to quantify the full strain field and determine mixed-mode stress intensity factors (SIFs) for crack tips in composite plates subjected to off-axis loading. The second method is formulated for orthotropic layers without a coat and it requires thermo-mechanical calibrations for two material parameters aligned with the material axes. The virtual crack closure technique (VCCT), Lekhnitskii’s and Savin’s elasticity solutions, and finite element (FE) analyses are used for demonstrations and validations of the second experimental method. The SIFs from the TSA methods are very sensitive to the uncertainty in the location of the crack tip and the unknown inelastic or damage zone size around the crack tip. The two experimental methods are effective in generating the strain fields around notched and other FRP composites.     

On the Mode I stress intensity factor for an external circular crack with fibre bridging

This paper presents a theoretical analysis of an external matrix crack located in a unidirectional fibre-reinforced elastic solid modelled as a transversely isotropic material. The presence of matrix cracking with fibre continuity introduces bridging action that has an influence on the stress intensity factors at the crack tip of the external crack. This paper presents a model for the bridged crack, where the fibre ligaments induce a constant displacement-dependent traction constraint over the external crack. This gives rise to a Fredholm integral equation of the second kind, which can be solved in an approximate fashion. We examine the specific problem where the bridged external circular crack is loaded by a doublet of concentrated forces. Numerical results are presented to illustrate the influence of the fibre–matrix modular ratio and the location of the loading on the bridged-crack opening mode stress intensity factor.     

Dynamic fracture mechanics study of an electrically impermeable mode III crack in a transversely isotropic piezoelectric material under pure electric load

The dynamic response of an electrically impermeable Mode III crack in a transversely isotropic piezoelectric material under pure electric load is investigated by treating the electric loading process as a transient impact load, which may be more appropriate to mimic the real service environment of piezoelectric materials. The stress intensity factor, the mechanical energy release rate, and the total energy release rate are derived and expressed as a function of time for a given applied electric load. The theoretical results indicate that a purely electric load can fracture the piezoelectric material if the stress intensity factor or the mechanical energy release rate is used as a failure criterion.     

Dielectric breakdown model for a conductive crack and electrode in piezoelectric materials

The strip dielectric breakdown (DB) model introduced by Zhang and Gao [T.Y. Zhang, C.F. Gao, Fracture behavior of piezoelectric materials, Thero. Appl. Fract. Mech. 41 (2004) 339–379] is used to study the generalized 2D problem of a conductive crack and an electrode in an infinite piezoelectric material. The energy release rate and stress intensity factors are derived based on the Stroh formalism, and then they are applied as failure criteria to predict the critical fracture loads. It is found that the DB strip may take the shielding effect on a conductive crack or electrode. For the case of an electrode, the local energy release rate and stress intensity factor become zero when DB happens ahead of the electrode tip. For the case of a mode-I conductive crack in a transversely isotropic piezoelectric solid, the results based on the DB model show that the critical stress intensity factor linearly increases as the applied electric field parallel to the poling direction increases, while it linearly decreases as the applied electric field anti-parallel to the poling direction increases. Finally, the upper and lower bounds of the actual critical fracture loads are proposed for a conductive crack in a piezoelectric material under combined mechanical–electrical loads.     

Stress concentration around a strongly oblate spheroidal cavity in a transversely isotropic medium

Expressions for the Eshelby tensor of a strongly oblate spheroidal region in a transversely isotropic medium are given explicitly. Based on the equivalent inclusion method, three dimensional stress concentration around the spherodal cavity subjected to remote uniform loading is analyzed and the associated stress concentration factor is determined. Analogously to two-dimensional blunted cracks, the so called stress rounding factor is introduced so that the connection between the crack tip stress and the stress intensity factor in linear fracture mechanics is established. Numerical values of the stress rounding factor for several representative cases of transversely isotropic symmetry are given.     

On higher order parameters in cracked composite plates under far‐field pure shear

This paper presents the analytical solution of the crack tip fields as well as the crack parameters in an infinitely large composite plate with a central crack subjected to pure shear loading. To this end, the complex variable method is employed to formulate an asymptotic solution for the crack tip fields in an anisotropic plane. Using a stress‐based definition of the crack tip modes of loading, only the mode II crack parameters are found to be non‐zero under pure shear load. Special focus is given to the determination of the higher order parameters of the crack tip asymptotic field, particularly the first non‐singular term, ie, the T‐stress. Unlike the isotropic materials, in which the T‐stress is zero under pure shear, it is found that the T‐stress is non‐zero for the case of anisotropic materials, being the only material‐dependent crack tip stress parameter. The veracity of our exact crack tip fields is assessed and verified through a comparison made with respect to the finite element (FE) solution. Finally, we demonstrate the significance of the T‐stress on stresses near the crack tip in composite plates under pure shear loads.     

Analysis of elliptical cracks perpendicular to the interface of two joined transversely isotropic solids

This paper presents a boundary element analysis of elliptical cracks in two joined transversely isotropic solids. The boundary element method is developed by incorporating the fundamental singular solution for a concentrated point load in a transversely isotropic bi-material solid of infinite space into the conventional displacement boundary integral equations. The multi-region method is used to analyze the crack problems. The traction-singular elements are employed to capture the singularity around the crack front. The values of stress intensity factors (SIFs) are obtained by using crack opening displacements. The results of the proposed method compare well with the existing exact solutions for an elliptical crack parallel to the isotropic plane of a transversely isotropic solid of infinite extent. Elliptical cracks perpendicular to the interface of transversely isotropic bi-material solids of either infinite extent or occupying a cubic region are further examined in detail. The crack surfaces are subject to the uniform normal tractions. The stress intensity factor values of the elliptical cracks of the two types are analyzed and compared. Numerical results have shown that the stress intensity factors are strongly affected by the anisotropy and the combination of the two joined solids.     

Elastodynamic analysis of the Finite punch and finite crack problems in orthotropic materials

The transient elastodynamic response of the finite punch and finite crack problems in orthotropic materials is examined. Solution for the stress intensity factor history around the punch corner and crack tip is found. Laplace and Fourier transforms together with the Wiener–Hopf technique are employed to solve the equations of motion in terms of displacements. A detailed analysis is made in the simplified case when a flat rigid punch indents an elastic orthotropic half-plane, the punch approaches with a constant velocity normally to the boundary of the half-plane. An asymptotic expression for the singular stress near the punch corner is analyzed leading to an explicit expression for the dynamic stress intensity factor which is valid for the time the dilatational wave takes to travel twice the punch width. In the crack problem, a finite crack is considered in an infinite orthotropic plane. The crack faces are loaded by impact uniform pressure in mode I. An expression for the dynamic stress intensity factor is found which is valid while the dilatational wave travels the crack length twice. Results for orthotropic materials are shown to converge to known solutions for isotropic materials derived independently.     

K‐dominance of static crack tip in functionally gradient materials

K‐dominance of static crack tip in functionally gradient materials (FGMs) with a crack oriented along the direction of the elastic gradient is studied through coherent gradient sensing (CGS), digital speckle correlation method (DSCM) and finite element method (FEM). In the direction of crack propagation, the shear modulus has a linear variation with constant mass density and Poisson's ratio. First, the CGS and DSCM governing equations related to the measurements and the elastic solutions at mode I crack in FGMs are obtained in terms of the stress intensity factor, material constants and graded index. Secondly, two kinds of FGMs specimens and one homogenous specimen are prepared to observe the influences of the property variation on the K‐dominance. Then, CGS and DSCM experiments using three‐point‐bending of FGMs and homogenous beams are performed. Thirdly, based on the results of the experiments, the stress intensity factors of three kinds of specimens are calculated by CGS and DSCM. Meanwhile, the stress intensity factors are obtained by FEM. Finally, comparing the results from CGS, DSCM and FEM, the K‐dominance of mode‐I static crack tip in FGMs is discussed in detail.     

A boundary element method for determining the effect of holes on the stress distribution around a crack

A boundary integral procedure is obtained for examining the effect of a finite number of holes on the crack tip stress intensity factors for a plane crack in a homogeneous anisotropic elastic material. Numerical results for specific examples involving a transversely isotropic material are given.     

An anisotropic layered material with a crack

Summary The problem of an anisotropic layered material which contains a plane crack in its interior is considered here. The problem is reduced to a set of Fredholm integral equations of the second kind which may be solved iteratively. Once these integral equations are solved, the crack tip stress intensity factors may be readily computed. Numerical results for some particular examples involving transversely isotropic materials are given here.     

Plane strain crack problems for elastic materials with variable properties

Distributions of stress, strain and displacement occurring at the tip of a crack in a material with properties dependent on the type of loading are investigated for the conditions of plane strain in both far-field tensile and shear loads. The causes of the dependence of material properties on the type of external forces are the various inhomogeneities such as microcracks, pores, inclusions or reinforcing components in a material. The behaviour of these inhomogeneities depends substantially on the conditions of loading or deformation. Hence, the deformation properties of a material are not fixed intrinsic material characteristics that are invariant to the loading conditions, but rather the macroproperties of such materials are stress-state-dependent ones, and this effect becomes more noticeable as the volume content of the inhomogeneities increases. The asymptotic solutions of crack problems are obtained on the basis of proposed stress-strain relations describing not only the stress-state dependence of material properties, but the interrelation between the characteristics of volume and shear deformation as well. In a non-uniform stress state the primary macrohomogeneous material becomes an heterogeneous one. The use of the stress function is not effective for the solution of plane strain crack problems for the materials under consideration. Therefore, an approach based on the corresponding representation for the strains is used. It is shown that the commonly used suppositions of the symmetry or anti-symmetry in the stress distribution relative to the crack plane can not be accepted, since they do not allow all the boundary conditions to be satisfied. The opening of the crack surfaces in the case of far shear field is observed. The influence of stress-state sensitivity of material properties on the values of the stress intensity factor is more significant for tensile crack than for the crack in far shear field.     

Plane strain crack problems for elastic materials with variable properties

Distributions of stress, strain and displacement occurring at the tip of a crack in a material with properties dependent on the type of loading are investigated for the conditions of plane strain in both far-field tensile and shear loads. The causes of the dependence of material properties on the type of external forces are the various inhomogeneities such as microcracks, pores, inclusions or reinforcing components in a material. The behaviour of these inhomogeneities depends substantially on the conditions of loading or deformation. Hence, the deformation properties of a material are not fixed intrinsic material characteristics that are invariant to the loading conditions, but rather the macroproperties of such materials are stress-state-dependent ones, and this effect becomes more noticeable as the volume content of the inhomogeneities increases. The asymptotic solutions of crack problems are obtained on the basis of proposed stress-strain relations describing not only the stress-state dependence of material properties, but the interrelation between the characteristics of volume and shear deformation as well. In a non-uniform stress state the primary macrohomogeneous material becomes an heterogeneous one. The use of the stress function is not effective for the solution of plane strain crack problems for the materials under consideration. Therefore, an approach based on the corresponding representation for the strains is used. It is shown that the commonly used suppositions of the symmetry or anti-symmetry in the stress distribution relative to the crack plane can not be accepted, since they do not allow all the boundary conditions to be satisfied. The opening of the crack surfaces in the case of far shear field is observed. The influence of stress-state sensitivity of material properties on the values of the stress intensity factor is more significant for tensile crack than for the crack in far shear field.     

An elastic-plastic finite element analysis of crack tip fields under biaxial loading conditions

A square plate containing a central crack and subjected to biaxial stresses has been studied by a finite element analysis. An elastic analysis shows that the crack opening displacement and stress of separation ahead of the crack tip are not affected by the mode of biaxial loading and therefore the stress intensity factor adequately describes the crack tip states in an elastic continuum.

An elastic-plastic analysis involving more than localized yielding at the crack tip provides different solutions of crack tip stress fields and crack face displacements for the different modes of biaxial loading.

The equi-biaxial loading mode causes the greatest separation stress but the smallest plastic shear ear and crack displacement. The shear loading system induces the maximum size of shear ear and crack displacement but the smallest value of crack tip separation stress.

     

Dynamic stress intensity factor of a functionally graded material under antiplane shear loading

Summary The dynamic response of a finite crack in an unbounded Functionally Graded Material (FGM) subjected to an antiplane shear loading is studied in this paper. The variation of the shear modulus of the functionally graded material is modeled by a quadratic increase along the direction perpendicular to the crack surface. The dynamic stress intensity factor is extracted from the asymptotic expansion of the stresses around the crack tip in the Laplace transform plane and obtained in the time domain by a numerical Laplace inversion technique. The influence of graded material property on the dynamic intensity factor is investigated. It is observed that the magnitude of dynamic stress intensity factor for a finite crack in such a functionally graded material is less than in the homogeneous material with a property identical to that of the FGM crack plane.     

T-stress based fracture model for cracks in isotropic materials

A modified version of the T-stress based fracture model, developed by Cotterell and Rice, is proposed. In this version an experimentally determined T_crit value is included in the model. The model is then applied to a branched crack in an isotropic material, and the direction of growth of the crack is predicted qualitatively. The branched crack problem is solved using the method of dislocations and a singular integral equation is obtained. The singular integral equation is solved using three different numerical techniques and their respective merits are discussed. The stress intensity factors and the T-stress in front of the branched crack tip are evaluated numerically. It is shown that the T-stress and the stress intensity factors are insensitive to the order of the singularity assumed at the reentrant wedge corner of the branched crack. For an uniaxial load and short kink length it is demonstrated that the kink will turn from its initial direction and realign with the main crack. However if the loading is biaxial then the direction of kink growth depends strongly on the applied transverse stress σ_xx. For a longer initial kink length the direction of kink growth depends on both the kink angle and loading.