Thermoelastic contact instability of a functionally graded layer and a homogeneous half-plane

The conductive heat transfer between two elastic bodies in the static contact can cause the system to be unstable due to the interaction between the thermoelastic distortion and pressure-dependent thermal contact resistance. This paper investigates the thermoelastic contact instability of a functionally graded material (FGM) layer and a homogeneous half-plane using the perturbation method. The FGM layer and half-plane are exposed to a uniform heat flux and are pressed together by a uniform pressure. The material properties of the FGM layer vary exponentially along the thickness direction. The characteristic equation governing the thermoelastic stability behavior is obtained to determine the stability boundary. The effects of the gradient index, layer thickness and material combination on the critical heat flux are discussed in detail through a parametric study. Results indicate that the thermoelastic stability behavior can be modified by adjusting the gradient index of the FGM layer.     

Theoretical modeling and analysis of thermal fracture of semi-infinite functionally graded materials with edge cracks

The present investigation is devoted to a problem of the interaction of two edge cracks inclined arbitrary to the boundary of a non-homogeneous half-plane, which is a functionally graded layer on a homogeneous substrate. The functionally graded properties vary exponentially in thickness direction. One cycle of cooling from sintering temperature is considered. An approach based on integral equations is used and a solution is obtained, then the stress intensity factors are calculated and direction of the initial crack propagation is evaluated by using the maximum circumferential stress criterion. Influence of geometrical and material (inhomogeneity) parameters on the fracture characteristics is investigated. This study can serve as a part of the modeling of the fracture process in FGM coatings under cyclic heating–cooling thermal loading.     

Dynamic stress intensity factors for two parallel interface cracks between a nonhomogeneous bonding layer and two dissimilar elastic half-planes subject to an impact load

Some composite materials are constructed of two dissimilar half-planes bonded by a nonhomogeneous elastic layer. In the present study, a crack is situated at the interface between the upper half-plane and the bonding layer of such a material, and another crack is located at the interface between the lower half-plane and the bonding layer. The material properties of the bonding layer vary continuously from those of the lower half-plane to those of the upper half-plane. Incoming shock stress waves impinge upon the two interface cracks normal to their surfaces. Fourier transformations were used to reduce the boundary conditions for the cracks to two pairs of dual integral equations in the Laplace domain. To solve these equations, the differences in the crack surface displacements were expanded in a series of functions that are zero-valued outside the cracks. The unknown coefficients in the series were solved using the Schmidt method so as to satisfy the conditions inside the cracks. The stress intensity factors were defined in the Laplace domain and were inverted numerically to physical space. Dynamic stress intensity factors were calculated numerically for selected crack configurations.     

Thermal stresses around a crack in the nonhomogeneous interfacial layer between two dissimilar elastic half-planes

Thermal stresses around a crack in the interfacial layer between two dissimilar elastic half-planes are solved. The surfaces of the crack are assumed to be insulated. The material constants of the interfacial layer are assumed to vary continuously from those of the upper half-plane to those of the lower half-plane. Uniform heat flows perpendicular the crack. Stress intensity factors are calculated numerically for several thicknesses of the interfacial layer.     

Thermoelastic instability of functionally graded materials with interaction of frictional heat and contact resistance

Both of the frictional heat and thermal contact resistance have a grave responsibility for the localized high temperature (hot spots) at the contact region, which is known as one of the most dangerous appearances in the brakes systems. In this paper, we study the thermoelastic instability (TEI) of a functionally graded material (FGM) half-plane sliding against a homogeneous half-plane at the in-plane direction. The interaction of the frictional heat and thermal contact resistance is taken into account in the TEI analysis. The material properties of the FGM half-plane are supposed to follow the exponential function along the thickness direction. The coupled TEI problem of FGMs is solved by using the perturbation method. The frictionally excited TEI of FGMs is also considered by neglecting the effect of the thermal contact resistance. The results show that the thermal contact resistance, sliding speed and gradient index have significant influence on the TEI. It is found that the variation of the gradient index of FGMs can increase the critical sliding speed and critical heat flux, and therefore improve the TEI of the sliding system.     

Boundary element analysis of crack problems in functionally graded materials

The present study examines the crack problems in a functionally graded material (FGM) whose upper and bottom surfaces are fully bonded with dissimilar homogeneous materials. A so-called generalized Kelvin solution based boundary element method is used in the numerical examination. The multi-region method and the eight-node traction-singular boundary elements are used for the crack evaluation. The layer discretization technique is utilized to approximate the depth material non-homogeneity of the FGM layer. The proposed method can deal with any depth variations in both the shear modulus and the Poisson ratio of the FGMs. Results of the present analysis are compared very well with the exact analytical solutions available in the literature, which demonstrates that the proposed method can accurately evaluate the stress intensity factors (SIFs) for cracks in FGMs. The paper further evaluates the effect of the functionally graded variations in the Poisson ratio on the stress intensity factors. The paper also assesses the elliptical cracks in the FGM system. The paper presents the influence of both the non-homogeneity and the thickness of the FGM layer on the three SIFs associated with the elliptical cracks.     

Mixed-mode fracture analysis of orthotropic functionally graded materials under mechanical and thermal loads

Mixed-mode fracture problems of orthotropic functionally graded materials (FGMs) are examined under mechanical and thermal loading conditions. In the case of mechanical loading, an embedded crack in an orthotropic FGM layer is considered. The crack is assumed to be loaded by arbitrary normal and shear tractions that are applied to its surfaces. An analytical solution based on the singular integral equations and a numerical approach based on the enriched finite elements are developed to evaluate the mixed-mode stress intensity factors and the energy release rate under the given mechanical loading conditions. The use of this dual approach methodology allowed the verifications of both methods leading to a highly accurate numerical predictive capability to assess the effects of material orthotropy and nonhomogeneity constants on the crack tip parameters. In the case of thermal loading, the response of periodic cracks in an orthotropic FGM layer subjected to transient thermal stresses is examined by means of the developed enriched finite element method. The results presented for the thermally loaded layer illustrate the influences of the material property gradation profiles and crack periodicity on the transient fracture mechanics parameters.     

Dynamics of a system comprising a pre-stressed orthotropic layer and pre-stressed orthotropic half-plane under the action of a moving load

This paper investigates the dynamic response to a moving load of a system comprising an initially stressed covering layer and initially stressed half-plane, within the scope of the piecewise-homogeneous body model utilizing three-dimensional linearized wave propagation theory in the initially stressed body. It was assumed that the materials of the layer and half-plane are anisotropic (orthotropic), and that the velocity of the line-located moving load is constant as it acts on the free face of the covering layer. The investigations were made for a two-dimensional problem (plane-strain state) under subsonic velocity of the moving load for complete and incomplete contact conditions. Corresponding numerical results were obtained for the stiffer layer and soft half-plane system in which the modulus of elasticity of the covering layer material (for the moving direction of the load) is greater than that of the half-plane material, which was assumed to be isotropic. Numerical results are presented and discussed for the critical velocity and stress distribution for various values of the problem parameters. In particular, it was established that, the critical velocity of the moving load is controlled mainly with a Rayleigh wave speed of a half-plane material and the initial stretching of the covering layer causes to increase these values.     

Dynamics of a system comprising an orthotropic layer and orthotropic half-plane under the action of an oscillating moving load

This paper investigates the dynamic response to a time-harmonic oscillating moving load of a system comprising a covering layer and half-plane, within the scope of the piecewise-homogeneous body model utilizing of the exact equations of the linear theory of elastodynamics. It is assumed that the materials of the layer and half-plane are anisotropic (orthotropic), and that the velocity of the line-located time-harmonic oscillating moving load is constant as it acts on the free face of the covering layer. Our investigations were carried out for a two-dimensional problem (plane-strain state) under subsonic velocity for a moving load in complete and incomplete contact conditions. The corresponding numerical results were obtained for the stiffer layer and soft half-plane system in which the modulus of elasticity of the covering layer material (for the moving direction of the load) is greater than that of the half-plane material. Numerical results are presented and discussed for the critical velocity, displacement and stress distribution for various values of the problem parameters. In particular, it is established that the critical velocity of the moving load is controlled mainly with a Rayleigh wave speed of a half-plane material and the existence of the oscillation of the moving load causes two types of critical velocity to appear: one of which is less, but the other one is greater than that attained for the case where the mentioned oscillation is absent.     

Fracture analysis of functionally graded coatings: plane deformation

A new multi-layered model for fracture analysis of functionally graded materials (FGMs) with arbitrarily varying elastic moduli under plane deformation has been developed. In this model, the FGM is divided into several sub-layers and in each sub-layer the shear modulus is assumed to be a linear function of the depth while the Poisson's ratio is assumed to be a constant. With this new model, an interface crack problem of a functionally graded coating bonded to a homogeneous half-plane under normal and shear loading is investigated. Employment of transfer matrix method and Fourier integral transform technique reduces the problem to a system of Cauchy singular integral equations which are solved numerically. Stress intensity factors of an interface crack are obtained for the cases of the shear modulus varying in an exponential manner and in a sinusoidal manner. Comparison of the present new model to other existing models shows that the new one is more efficient.     

Driving forces and kinking of an oblique crack in bonded nonhomogeneous materials

Summary  Plane elasticity solutions are presented for the problem of an oblique crack in two bonded media. The material model under consideration consists of a homogeneous half-plane with an arbitrarily oriented crack and a nonhomogeneous half-plane. The Fourier integral transform method is employed in conjunction with the coordinate transformations of field variables in the basic elasticity equations. Formulation of the crack problem results in having to solve a system of singular integral equations for arbitrary crack surface tractions. A crack perpendicular to or along the bonded interface between the homogeneous and nonhomogeneous constituents arises as a limiting case. In the numerical results, the values of mixed-mode stress intensity factors are provided for various combinations of relevant geometric and material parameters of the bonded media. Subsequently, the infinitesimal kinks from the tips of a main crack are presumed, with the corresponding local driving forces being evaluated in terms of the stress intensities of the main crack. The criterion of maximum energy release rate is applied with the aim of making some conjectures concerning the likelihood of kinking and the probable kink direction based on the approximation of local homogeneity and brittleness of the crack-tip behavior. Received 25 September 2001; accepted for publication 13 February 2002     

Electroelastic Analysis of an Internal Interface Crack in a Half-Plane Consisting of Two Bonded Dissimilar Piezoelectric Quarter-Planes

The problem of an interface crack in a half-plane consisting of two bonded dissimilar piezoelectric quarters is considered under antiplane shear and inplane electric loading. The problem is solved under the electrically permeable assumption for a crack. The integral transform technique is employed to reduce the problem to triple integral equations, which is further converted to a hypersingular integral equation for the crack sliding displacement. By solving the resulting equation analytically, the electroelastic field along the interface and the energy release rate are obtained in explicit form, respectively. Several examples are given to illustrate the influence of the material properties and the crack position on the energy release rate.     

Thermally induced fracture of a smart functionally graded composite structure

Discussed is the fracture behavior of a cracked smart actuator on a substrate under thermal load. The actuator is made of piezoelectric material with functionally graded material (FGM) properties. Integral transform method is used to reduce the problem to the solution of a set of singular integral equations and is solved numerically. This paper is completed by including graphical plots of the thermal flow, stress and electric displacement intensity factors around the crack for different crack positions and material gradients. Directions of crack initiation are also predicted by using the energy density criterion.     

Two-dimensional thermoelastic contact problem of functionally graded materials involving frictional heating

The two-dimensional thermoelastic sliding frictional contact of functionally graded material (FGM) coated half-plane under the plane strain deformation is investigated in this paper. A rigid punch is sliding over the surface of the FGM coating with a constant velocity. Frictional heating, with its value proportional to contact pressure, friction coefficient and sliding velocity, is generated at the interface between the punch and the FGM coating. The material properties of the coating vary exponentially along the thickness direction. In order to solve the heat conduction equation analytically, the homogeneous multi-layered model is adopted for treating the graded thermal diffusivity coefficient with other thermomechanical properties being kept as the given exponential forms. The transfer matrix method and Fourier integral transform technique are employed to convert the problem into a Cauchy singular integral equation which is then solved numerically to obtain the unknown contact pressure and the in-plane component of the surface stresses. The effects of the gradient index, Peclet number and friction coefficient on the thermoelastic contact characteristics are discussed in detail. Numerical results show that the distribution of the contact stress can be altered and therefore the thermoelastic contact damage can be modified by adjusting the gradient index, Peclet number and friction coefficient.     

Thermal fracture of a functionally graded/homogeneous bimaterial with system of cracks

The thermal fracture of a bimaterial consisting of a homogeneous material and a functionally graded material (FGM) with a system of internal cracks and an interface crack is investigated. The bimaterial is subjected to a heat flux. The thermal properties of FGM are assumed to be continues functions of the thickness coordinate, while the elastic properties are constants. The method of the solution is based on the singular integral equations. For a special case where the interface crack is much larger than the internal cracks in the FGM the asymptotic analytical solution of the problem is obtained as series in a small parameter (the ratio between sizes of the internal and interface crack) and the thermal stress intensity factors (TSIFs) are derived as functions of geometry of the problem and material characteristics. A parametric analysis of the effects of the location and orientation of the cracks and of the inhomogeneity parameter of FGM’s thermal conductivity on the TSIFs is performed. The results are applicable to such kinds FGMs as ceramic/ceramic FGMs, e.g., TiC/SiC, MoSi₂/Al₂O₃ and MoSi₂/SiC, and also some ceramic/metal FGMs.     

竖向裂缝对功能梯度板的热响应影响

采用条单元法分析了竖向裂缝对功能梯度板的热响应影响。功能梯度板的材料属性沿板厚方向连续线性变化,温度变化热源作用于带竖向裂缝的功能梯度板上。通过研究得出了裂缝距离热源不同位置、不同长度裂缝及不同频率温度变化热源等状态时,功能梯度板上表面的热响应情况。采用条单元法研究带竖向裂缝的功能梯度板的热响应问题,为带裂缝的功能梯度板的热研究提供了一定的指导。     

半平面多边缘裂纹反平面问题的奇异积分方程

利用复变函数和奇异积分方程方法,求解弹性范围内半平面多边缘裂纹的反平面问题．提出了满足半平面边界自由的由分布位错密度表示的单边缘裂纹的基本解,此基本解由主要部分和辅助部分组成．将半平面多边缘裂纹问题看作是许多单边缘裂纹问题的叠加,建立了一组Cauchy型奇异积分方程．然后,利用半开型积分法则求解该奇异积分方程,得到了裂纹端处的应力强度因子．最后,给出了几个数值算例．     

STRESS INTENSITY FACTOR OF AN ANTI-PLANE CRACK PARALLEL TO THE WEAK/MICRO-DISCONTINUOUS INTERFACE IN A BI-FGM COMPOSITE

The problem considered is a mode Ⅲ crack lying parallel to the interface of an exponential-type functional graded material （FGM） strip bonded to a linear-type FGM substrate with infinite thickness. By applying the Fourier integral transform, the problem was reduced as a Cauchy singular integral equation with an unknown dislocation density function. The collocation method based on Chebyshev polynomials proposed by Erdogan and Gupta was used to solve the singular integral equation numerically. With the numerical solution, the effects of the geometrical and physical parameters on the stress intensity factor （SIF） were analyzed and the following conclusions were drawn： （a） The region affected by the interface or free surface varies with the material rigidity, and higher material rigidity will lead to bigger affected region. （b） The SIF of the crack in the affected region and parallel to the micro-discontinuous interface is lower than those of the weak discontinuous cases. Reducing the weak-discontinuity of the interface will be beneficial to decrease the SIF of the interface-parallel crack in the region affected by the interface. （c） The effect of the free surface on SIF is more remarkable than that of the interface, and the latter is still more notable than that of the material rigidity. When the effects of the interface and free surface are fixed, increase of the material rigidity will enhance the value of SIF.     

Anti-plane shear of an arbitrary oriented crack in a functionally graded strip bonded with two dissimilar half-planes

An internal crack located within a functionally graded material (FGM) strip bonded with two dissimilar half-planes and under an anti-plane load is considered. The crack is oriented in an arbitrary direction. The material properties of strip are assumed to vary exponentially in the thickness direction and two half-planes are assumed to be isotropic. Governing differential equations are derived and to reduce the difficulty of the problem dealing with solution of a system of singular integral equations Fourier integral transform is employed. Semi closed form solution for the stress distribution in the medium is obtained and mode III stress intensity factor (SIF), at the crack tip is calculated and its validity was verified. Finally, the effects of nonhomogeneous material parameter and crack orientation on the stress intensity factor are studied.     

Effects of material gradients on transient dynamic mode-III stress intensity factors in a FGM

This paper presents a transient dynamic crack analysis for a functionally graded material (FGM) by using a hypersingular time-domain boundary integral equation method. The spatial variations of the material parameters of the FGM are described by an exponential law. A numerical solution procedure is developed for solving the hypersingular time-domain traction BIE. To avoid the use of time-dependent Green’s functions which are not available for general FGM, a convolution quadrature formula is adopted for approximating the temporal convolution, while a Galerkin method is applied for the spatial discretization of the hypersingular time-domain traction BIE. Numerical results for the transient dynamic stress intensity factors for a finite crack in an infinite and linear elastic FGM subjected to an impact anti-plane crack-face loading are presented and discussed. The effects of the material gradients of the FGM on the transient dynamic stress intensity factors and their dynamic overshoot over the corresponding static stress intensity factors are analyzed.