Interfacial fracture analysis of bonded dissimilar strips with a functionally graded interlayer under antiplane deformation

This paper provides the solution to the problem of dissimilar, homogeneous semi-infinite strips bonded through a functionally graded interlayer and weakened by an embedded or edge interfacial crack. The bonded system is assumed to be under antiplane deformation, subjected to either traction-free or clamped boundary conditions along its bounding planes. Based on the Fourier integral transform, the problem is formulated in terms of a singular integral equation which has a simple Cauchy kernel for the embedded crack and a generalized Cauchy kernel for the edge crack. In the numerical results, the effects of geometric and material parameters of the bonded system on the crack-tip stress intensity factors are presented in order to quantify the interfacial fracture behavior in the presence of the graded interlayer.     

Thermoelastic problem of steady-state heat flows disturbed by a crack at an arbitrary angle to the graded interfacial zone in bonded materials

Plane thermoelasticity solutions are presented for the problem of a crack in bonded materials with a graded interfacial zone. The interfacial zone is treated as a nonhomogeneous interlayer having spatially varying thermoelastic moduli between dissimilar, homogeneous half-planes. The crack is assumed to exist in one of the half-planes at an arbitrary angle to the graded interfacial zone, disturbing uniform steady-state heat flows. The Fourier integral transform method is employed in conjunction with the coordinate transformations of field variables in the basic thermoelasticity equations. Formulation of the current nonisothermal crack problem lends itself to the derivation of two sets of Cauchy-type singular integral equations for heat conduction and thermal stress analyses. The heat-flux intensity factors and the thermal-stress intensity factors are defined and evaluated in order to quantify the singular characters of temperature gradients and thermal stresses, respectively, in the near-tip region. Numerical results include the variations of such crack-tip field intensity factors versus the crack orientation angle for various combinations of material and geometric parameters of the dissimilar media bonded through the thermoelastically graded interfacial zone. The dependence of the near-tip thermoelastic singular field on the degree of crack-surface partial insulation is also addressed.     

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.     

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.     

An axisymmetric problem of an embedded crack in a graded layer bonded to a homogeneous half-space

In an attempt to simulate non-uniform coating delamination, the elasto-static problem of a penny shaped axisymmetric crack embedded in a functionally graded coating bonded to a homogeneous substrate subjected to crack surface tractions is considered. The coating’s material gradient is parallel to the axisymmetric direction and is orthogonal to the crack plane. The graded coating is modeled as a non-homogeneous medium with an isotropic constitutive law. Using Hankel transform, the governing equations are converted into coupled singular integral equations, which are solved numerically to yield the crack tip stress intensity factors. The Finite Element Method was additionally used to model the crack problem. The main objective of this paper is to study the influence of the material non-homogeneity and the crack position on the stress intensity factors for the purpose of gaining better understanding on the behavior of graded coatings.     

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.     

Dynamic stress intensity factor of a crack perpendicular to the weak-discontinuous interface in a nonhomogeneous coating–substrate structure

A mechanical model was established for the antiplane dynamic fracture problem of a functionally graded coating–substrate structure with a coating crack perpendicular to the weak-discontinuous interface. The problem was reduced to a Cauchy singular integral equation by the methods of Laplace and Fourier integral transforms. Erdogan’s collocation method and the Laplace numerical inversion proposed by Miller and Guy were used to calculate the dynamic stress intensity factors. Three conclusions were drawn through parametric studies: (a) unlike the conclusion drawn for an interfacial crack, reducing the weak discontinuity of the interface will not necessarily decrease the dynamic stress intensity factor (DSIF) of the coating crack perpendicular to the interface; (b) increasing the stiffness of the substrate when that of the coating is fixed, or decreasing the stiffness of coating when that of the substrate is fixed, will be beneficial for the reduction of the DSIF of a coating crack perpendicular to the interface; and (c) the free surface has a greater influence on the DSIF than the interface does, and the effect of the interface on the DSIF is greater than that of the material stiffness in the crack-tip region.     

The dynamic response of an edge crack in a functionally graded orthotropic strip

The dynamic response of a functionally graded orthotropic strip with an edge crack perpendicular to the boundaries is studied. The material properties are assumed to vary continuously along the thickness direction. Laplace and Fourier transforms are applied to reduce the problem to a singular integral equation. Numerical results are presented to illustrate the influences of parameters such as the nonhomogeneity constant and geometry parameters on the dynamic stress intensity factors (SIFs).     

Antiplane crack at the interface of two bonded dissimilar graded piezoelectric materials

The problem of an antiplane crack situated in the interface of two bonded dissimilar graded piezoelectric half-spaces is considered under the permeable crack assumption. The mechanical and electrical properties of the half-spaces are considered for a class of functional forms for which the equilibrium equation has analytic solutions. By using an integral transform technique, the problem is reduced to dual integral equations which are transformed into a Fredholm integral equation by introducing an auxiliary function. The stress intensity factors are obtained in explicit form in terms of auxiliary functions. By solving the Fredholm integral equation numerically, the numerical results for stress intensity factors are obtained which have been displayed graphically to show the influence of the graded piezoelectric materials.     

功能梯度材料涂层平面裂纹分析

研究粘接于均质基底材料上功能梯度涂层平面裂纹问题. 假设功能梯度材料剪切模量的倒数为坐标的线性函数，而泊松比为常数. 采用Fourier变换和传递矩阵法将该混合边值问题化为奇异积分方程组，通过数值求解获得 应力强度因子. 考察了材料梯度变化形式、结构几何尺寸和材料梯度参数对裂纹应力强度因子的影响，发现 功能梯度材料涂层尺寸、裂纹长度以及材料梯度参数均对应力强度因子有显著影响.     

Cylindrical interface crack in a hollow layered functionally graded cylinder under static torsion

A hollow functionally graded composite cylinder under static torsion, which consists of an inner and outer elastic circular tube with a cylindrical interface crack, is studied in this work. By utilizing Fourier integral transform method, the mixed boundary value problem is reduced to a Cauchy singular integral equation, from which the numerical results of the stress intensity factor are obtained by the Lobatto–Chebyshev quadrature technique. Numerical results demonstrate the coupled effects of geometrical, physical, and functionally graded parameters on the interfacial fracture behavior.     

The Anti-plane Dynamic Fracture Analysis of Functionally Graded Materials with Arbitrary Spatial Variations of Material Properties

Anti-plane dynamic fracture analysis is presented for functionally graded materials (FGM) with arbitrary spatial variations of material properties. The FGM with the material properties varying continuously in an arbitrary manner is modeled as a multi-layered medium with the elastic modulus and mass density varying linearly in each sub-layer and continuous at the interfaces between two adjacent sub-layers. With this linearly inhomogeneous multi-layered model, the problem of a crack in a graded interfacial zone bonded to two homogeneous half-spaces or in a coating bonded to a homogeneous half-space subjected to the anti-plane shear impact load is investigated. Laplace and Fourier transforms and transfer matrix are applied to reduce the associated mixed boundary value problem to a Cauchy singular integral equation which is solved numerically in the Laplace transformed domain. The dynamic stress intensity factors (DSIF) are obtained by using the numerical technique of Laplace inversion.     

Dynamic investigation of a functionally graded layered structure with a crack crossing the interface

The dynamic response of a functionally graded layered structure with a crack crossing the interface is analyzed. The in-plane impact loading condition is considered. By using the Laplace and Fourier integral transforms, singular integral equation method and residue theory, the present problem is reduced to a singular integral equation in the Laplace transform domain. The influences of Young’s modulus ratio, thickness ratio, and crack length and location on the dynamic stress intensity factors (DSIFs) are investigated. Particularly, the DSIFs corresponding to different crack locations are shown in the case when the crack center moves from one layer to another layer through the interface. The peak and static values and overshoot characteristics of the DSIFs are analyzed. It is found that these values typically exhibit kinking behavior when the crack tips arrive at the interface. This study is different from previous other investigations in the following respects: (1) the dynamic response of a crack crossing the interface of a functionally graded structure is studied analytically, which has hardly been done in the past and (2) the present model can be reduced to some important problems, such as a functionally graded coating-substrate structure with a crack in the graded coating or homogeneous substrate or one intersecting the interface.     

Elliptical crack normal to functionally graded interface of bonded solids

This paper presents an analysis of an elliptical crack that is perpendicular to a functionally graded interfacial zone between two fully bonded solids. The functionally graded interfacial zone is treated as a non-homogeneous solid layer with its elastic modulus varying in the thickness direction. A generalized Kelvin solution based boundary element method is employed for the calculation of the stress intensity factors associated with the three-dimensional crack problem. The elliptical crack surface is subject to either uniform normal traction or uniform shear traction. The stress intensity factors are examined by taking into account the effects of the non-homogeneity parameter and thickness of the functionally graded interfacial zone, as well as the crack distance to the zone. The SIF values are further incorporated into the S-criterion for prediction of crack growth. The paper presents the most possible direction and location of the elliptical crack growth under an inclined tensile (or compressive) load. The paper further presents results of the critical external loads that would cause the elliptical crack to grow at the most possible location and along the most possible direction. The paper also examines the effects of external load direction and material and geometrical parameters on the critical loads.     

Fracture analysis of mode-II crack perpendicular to imperfect bimaterial interface

The problem of a mode-II crack close to and perpendicular to an imperfect interface of two bonded dissimilar materials is investigated.The imperfect interface is modelled by a linear spring with the vanishing thickness.The Fourier transform is used to solve the boundary-value problem and to derive a singular integral equation with the Cauchy kernel.The stress intensity factors near the left and right crack tips are evaluated by numerically solving the resulting equation.Several special cases of the mode-II crack problem with an imperfect interface are studied in detail.The effects of the interfacial imperfection on the stress intensity factors for a bimaterial system of aluminum and steel are shown graphically.The obtained observation reveals that the stress intensity factors are dependent on the interface parameters and vary between those with a fully debonded interface and those with a perfect interface.     

平行于功能梯度材料夹层的币型裂纹起裂条件

分析了功能梯度材料中币型裂纹的扩展问题．裂纹平行于无限域中功能梯度材料夹层,受有与裂纹面成任意角度的拉应力．假定功能梯度材料夹层与两个半无限域均匀介质完全粘合,其弹性模量沿厚度方向变化．采用基于层状材料广义Kelvin基本解的边界元方法分析裂纹问题,给出了均布正应力和剪应力作用下裂纹的应力强度因子、将应力强度因子耦合于应变能密度断裂判据,讨论了裂纹体在拉伸应力作用下的起裂条件．     

Mode III eccentric crack in a functionally graded piezoelectric strip

This paper studies the internal crack problem located within one functionally graded piezoelectric strip. One crack is normal to the edge of the strip and the material properties vary along the direction of crack length. Three different boundary conditions and both impermeable and permeable cases are discussed. The problem can be reduced to a system of singular integral equations and solved by using the Gauss–Chebyshev formulas. The results show that the edge boundary conditions and the nonhomogeneous parameter significantly control the magnitudes of stress and electric displacement intensity factors.     

Diffraction of torsional elastic waves by a peripheral edge crack around a spherical cavity

This paper deals with the problem of finding the stress distribution in the neighborhood of a peripheral edge crack in a spherical cavity. The crack is excited by a torsional standing wave.The problem is solved by using integral transforms and is reduced to the solution of a singular integral equation. The solution of this equation is obtained numerically by the method due to Erdogan, Gupta, and Cook, and the stress intensity factors are displayed graphically.     

Thermal shock analyses for a strip containing an infinite row of periodically distributed cracks normal to its edge

The transient thermal stress problem of a semi-infinite plate containing an infinite row of periodically distributed cracks normal to its edge is investigated in this paper. The elastic medium is assumed to be cooled suddenly on the crack-containing edge. By the superposition principle, the formulation leads to a mixed boundary value problem, with the negating tractions arisen from the thermal stresses for a crack-free semi-infinite plate. The resulting singular integral equation is solved numerically. The effects on the stress intensity factors due to the presence of periodically distributed cracks in a semi-infinite plate are illustrated. For both the edge crack and the embedded crack arrays, the stress intensity factors increase, due to the reduction of the shielding effect, as the stacking cracks are more separated. For the case of embedded crack array, one has the further conclusion that the stress intensity factors decline as the crack array shifts from the plate edge.     

Anti-plane dynamic hole–crack interaction in a functionally graded piezoelectric media

The anti-plane dynamic problem of a functionally graded piezoelectric plane containing a hole–crack system is treated by a non-hypersingular traction-based boundary integral equation method. The material parameters vary exponentially in the same manner in an arbitrary direction. The system is loaded by an incident SH-type wave, and impermeable boundary conditions are assumed. Using a frequency-dependent fundamental solution of the wave equation, the boundary value problem is transformed into a system of integro-differential equations along the boundary of the hole and on the crack line. Its numerical solution yields the dynamic stress intensity factors and stress concentration factors. A parametric study reveals their dependence on the hole–crack scenario and its geometry, characteristics of the dynamic load and magnitude and direction of material inhomogeneity.