A partially insulated interface crack between a graded orthotropic coating and a homogeneous orthotropic substrate under heat flux supply

In this paper, thermal response of an orthotropic functionally graded coating-substrate structure with a partially insulated interface crack under heat flux supply is considered. It is assumed that there exists thermal resistance to heat conduction through the crack region. The mixed boundary value problems are reduced to a system of singular integral equations and solved numerically. Higher order asymptotic terms for the singular integral kernels are considered to improve the accuracy and the convergence efficiency of the numerical integrals. Numerical results are presented to show the effects of the orthotropy parameters, thermo-elastic nonhomogeneity parameters, and dimensionless thermal resistance on the temperature distribution and the thermal stress intensity factors (SIFs).     

Thermal stress intensity factors for an interface crack in a functionally graded layered structures

The fracture behavior of a functionally graded layered structure (FGLS) with an interface crack under thermal loading is investigated. Considering new boundary conditions, it is assumed that interface crack is partly insulated, and the temperature drop across the crack surfaces is the result of the thermal resistance due to the heat conduction through the crack region. The problem is formulated in terms of a system of singular integral equations. Numerical results are presented to show the influence of the material nonhomogeneity parameters and the dimensionless thermal resistance on the thermal stress intensity factors (TSIFs).     

Thermal response of a partially insulated interface crack in a graded coating-substrate structure under thermo-mechanical disturbance: Energy release and density

Consider the thermal fracture problem of a functionally graded coating-substrate structure of finite thickness with a partially insulated interface crack subjected to thermal-mechanical supply. A new model is proposed that the heat conduction through the crack region occurs and the temperature drop across the crack surfaces is the result of the thermal resistance. For the first time, real fundamental solutions are derived for the fracture analysis of functionally graded materials. The complicated mixed boundary problems of equations of heat conduction and elasticity are converted analytically into singular integral equations, which are solved numerically. The asymptotic expressions with higher order terms for the singular integral kernels are considered to improve the accuracy and efficiency of the numerical integration. Explicit expressions of various failure modes including stress intensity factors, energy release rate and strain energy density, are provided. Numerical results are presented to illustrate the effects of non-homogeneity parameters and the dimensionless thermal resistance on the temperature distribution along the crack surfaces and extended crack line, the thermal stress intensity factors and minimum strain energy density.     

Buckling of an orthotropic graded coating with an embedded crack bonded to a homogeneous substrate

To simulate buckling of nonuniform coatings, we consider the problem of an embedded crack in a graded orthotropic coating bonded to a homogeneous substrate subjected to a compressive loading. The coating is graded in the thickness direction and the material gradient is orthogonal to the crack direction which is parallel with the free surface. The elastic properties of the material are assumed to vary continuously along the thickness direction. The principal directions of orthotropy are parallel and perpendicular to the crack orientation. The loading consists of a uniform compressive strain applied away from the crack region. The graded coating is modeled as a nonhomogeneous medium with an orthotropic stress–strain law. Using a nonlinear continuum theory and a suitable perturbation technique, the plane strain problem is reduced to an eigenvalue problem describing the onset of buckling. Using integral transforms, the resulting plane elasticity equations are converted analytically into singular integral equations which are solved numerically to yield the critical buckling strain. The Finite Element Method was additionally used to model the crack problem. The main objective of the paper is to study the influence of material nonhomogeneity on the buckling resistance of the graded layer for various crack positions, coating thicknesses and different orthotropic FGMs.     

正交各向异性功能梯度材料反平面裂纹尖端应力场

采用积分变换-对偶积分方程方法,研究了正交各向异性功能梯度材料反平面裂纹问题,文中假定材料沿两个主轴方向的剪切模量成比例按双参数梯度模型变化,通过求解对偶积分程并考虑变形Bessel函数的渐特性,推导出了裂纹尖端应力场,最后考察了材料非均匀性及正交性对应力强度因子的影响。     

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.     

基于近场动力学微分算子的含裂纹正交各向异性板热传导分析方法

采用近场动力学微分算子(Peridynamic Differential Operator, PDDO)理论建立正交各向异性板热传导的非局部模型。通过构造近场动力学函数,将边界条件和热传导方程由局部微分形式转化为非局部积分形式,引入Lagrange乘子,采用变分分析对含裂纹正交各向异性板温度及裂纹尖端的热通量分布进行求解。通过对比算例,验证了该模型具有较好的收敛性和有效性。分析了正交各向因子、材料铺设角、裂纹倾角及间距对裂纹尖端热通量的影响。结果表明,基于PDDO建立的含裂纹正交各向异性板热传导模型,考虑了热传导问题中的非局部性,能有效提高计算精度,预测含裂纹板中裂纹尖端出现的奇异性。     

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.     

An exact solution for the steady-state thermoelastic response of functionally graded orthotropic cylindrical shells

We analyze the steady-state response of a functionally graded thick cylindrical shell subjected to thermal and mechanical loads. The functionally graded shell is simply supported at the edges and it is assumed to have an arbitrary variation of material properties in the radial direction. The three-dimensional steady-state heat conduction and thermoelasticity equations, simplified to the case of generalized plane strain deformations in the axial direction, are solved analytically. Suitable temperature and displacement functions that identically satisfy the boundary conditions at the simply supported edges are used to reduce the thermoelastic equilibrium equations to a set of coupled ordinary differential equations with variable coefficients, which are then solved by the power series method. In the present formulation, the cylindrical shell is assumed to be made of an orthotropic material, although the analytical solution is also valid for isotropic materials. Results are presented for two-constituent isotropic and fiber-reinforced functionally graded shells that have a smooth variation of material volume fractions, and/or in-plane fiber orientations, through the radial direction. The cylindrical shells are also analyzed using the Flügge and the Donnell shell theories. Displacements and stresses from the shell theories are compared with the three-dimensional exact solution to delineate the effects of transverse shear deformation, shell thickness and angular span.     

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.     

Thermal stresses factors of orthotropic graded layers with a finite line bond

In the paper, the problem of a finite line bond between two orthotropic functionally graded strips under thermal loading is considered. Considering some new boundary conditions, it is assumed that the temperature drop across the finite line bond is the result of the thermal conductivity index controlling heat conduction through the bond region. Using Fourier transforms technique, the therm-elastic mixed boundary value problems are reduced to a system of singular integral equations which can be solved approximately by applying the Chebyshev polynomials. The numerical results for the temperature and displacement field as well as the thermal stress intensity factors (TSIFs) are presented. The influence of the thickness of the layers and the thermo-elastic nonhomogeneity parameters on the temperature distribution and the TSIFs is discussed in detail. These results can be expected to be used for the purpose of gaining better understanding of the thermo-mechanical behavior of layered structures.     

Micromechanics-based thermoelastic model for functionally graded particulate materials with particle interactions

Thermoelastic behavior of functionally graded particulate materials is investigated with a micromechanical approach. Based on a special representative volume element constructed to represent the graded microstructure of a macroscopic material point, the relation between the averaged strains of the particle and matrix phases is derived with pair-wise particle interactions, and a set of governing equations for the thermoelastic behavior of functionally graded materials is presented. The effective coefficient of thermal expansion at a material point is solved through the overall averaged strain of two phases induced by temperature change under the stress-free condition, and is shown to exhibit a weak anisotropy due to the particle interactions within the graded microstructure. When the material gradient is eliminated, the proposed model predicts the effective coefficient of thermal expansion for uniform composites as expected. If the particle interactions are disregarded, the proposed model recovers the Kerner model. The proposed semi-analytical scheme is consistent and general, and can handle any thermal loading variation. As examples, the thermal stress distributions of graded thermal barrier coatings are solved for two types of thermal loading: uniform temperature change and steady-state heat conduction in the gradation direction.     

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.     

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.     

Dynamic stress intensity factor for cracked functionally graded orthotropic medium under time-harmonic loading

Dynamic stress intensity factor for a Griffith crack in functionally graded orthotropic materials under time-harmonic loading is investigated in the present paper. By using the Fourier transform and defining the jumps of displacement components across the crack surface as the unknown functions, two pairs of dual integral equations are derived. To solve the dual integral equations, the jumps of the displacement components across the crack surface are expanded in a series of Jacobi polynomial. Numerical examples are provided to show the effects of material properties and the crack configuration on the dynamic stress intensity factors of the functionally graded orthotropic materials with a Griffith crack.     

功能梯度材料热传导问题的近场动力学模型

当材料中存在不连续性缺陷的时候,传统的基于连续介质理论的方法在研究热传导问题时存在诸多不便．我们基于近场动力学方法(Peridynamics)建立了功能梯度材料的热传导模型,并且研究了在温度荷载作用下功能梯度材料的温度场的变化表现．该文概述了PD方法应用于热传导问题时的详细理论基础,描述了其建模思路以及计算体系,给出了使用PD方法模拟结构承受温度荷载时的计算格式,讨论了不同梯度形式对功能梯度材料内热传导的影响．算例结果表明：通过PD模型所得到的温度场与解析解吻合较好,证明了PD方法在分析功能梯度材料热传导行为等问题时的可行性．     

DYNAMIC STRESS FIELD AROUND THE MODE Ⅲ CRACK TIP IN AN ORTHOTROPIC FUNCTIONALLY GRADED MATERIAL

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎｒｅｃｅｎｔｙｅａｒｓ,ｇｒｅａｔａｔｔｅｎｔｉｏｎｓｈａｖｅｂｅｅｎｐａｉｄｔｏｔｈｅｒｅｓｅａｒｃｈｏｆＦｕｎｃｔｉｏｎａｌｌｙＧｒａｄｅｄＭａｔｅｒｉａｌｓ(ＦＧＭ).Ｆｒｏｍｔｈｅｖｉｅｗｐｏｉｎｔｓｏｆａｐｐｌｉｅｄｍｅｃｈａｎｉｃｓ,ＦＧＭａｒｅｎｏｎ＿ｈｏｍｏｇｅｎｅｏｕｓｓｏｌｉｄｓ.Ｔｈｅｎｏｎ＿ｈｏｍｏｇｅｎｅｉｔｙｏｆＦＧＭｈａｓａｇｒｅａｔｉｎｆｌｕｅｎｃｅｏｎｔｈｅｉｒｍｅｃｈａｎｉｃａｌｂｅｈａｖｉｏｒ,ｅｓｐｅｃｉａｌｌｙｗｈｅｎｔｈｅｃｏｍｐ…     

Thermomechanical vibration analysis of a functionally graded shell with flowing fluid

This paper reports the results of an investigation into the vibration of functionally graded cylindrical shells with flowing fluid, embedded in an elastic medium, under mechanical and thermal loads. By considering rotary inertia, the first-order shear deformation theory (FSDT) and the fluid velocity potential, the dynamic equation of functionally graded cylindrical shells with flowing fluid is derived. Here, heat conduction equation along the thickness of the shell is applied to determine the temperature distribution and material properties are assumed to be graded distribution along the thickness direction according to a power-law in terms of the volume fractions of the constituents. The equations of eigenvalue problem are obtained by using a modal expansion method. In numerical examples, effects of material composition, thermal loading, static axial loading, flow velocity, medium stiffness and shell geometry parameters on the free vibration characteristics are described. The new features in this paper are helpful for the application and the design of functionally graded cylindrical shells containing fluid flow.     

Chaotic vibrations of an orthotropic FGM rectangular plate based on third-order shear deformation theory

In this paper, an analysis on the nonlinear dynamics and chaos of a simply supported orthotropic functionally graded material (FGM) rectangular plate in thermal environment and subjected to parametric and external excitations is presented. Heat conduction and temperature-dependent material properties are both taken into account. The material properties are graded in the thickness direction according to a simple power law distribution in terms of the volume fractions of the constituents. Based on the Reddy’s third-order share deformation plate theory, the governing equations of motion for the orthotropic FGM rectangular plate are derived by using the Hamilton’s principle. The Galerkin procedure is applied to the partial differential governing equations of motion to obtain a three-degree-of-freedom nonlinear system. The resonant case considered here is 1:2:4 internal resonance, principal parametric resonance-subharmonic resonance of order 1/2. Based on the averaged equation obtained by the method of multiple scales, the phase portrait, waveform and Poincare map are used to analyze the periodic and chaotic motions of the orthotropic FGM rectangular plate. It is found that the motions of the orthotropic FGM plate are chaotic under certain conditions.