An interaction integral method for analysis of cracks in orthotropic functionally graded materials

This paper presents two new interaction integrals for calculating stress-intensity factors (SIFs) for a stationary crack in two-dimensional orthotropic functionally graded materials of arbitrary geometry. The method involves the finite element discretization, where the material properties are smooth functions of spatial co-ordinates and two newly developed interaction integrals for mixed-mode fracture analysis. These integrals can also be implemented in conjunction with other numerical methods, such as meshless method, boundary element method, and others. Three numerical examples including both mode-I and mixed-mode problems are presented to evaluate the accuracy of SIFs calculated by the proposed interaction integrals. Comparisons have been made between the SIFs predicted by the proposed interaction integrals and available reference solutions in the literature, generated either analytically or by finite element method using various other fracture integrals or analyses. An excellent agreement is obtained between the results of the proposed interaction integrals and the reference solutions. The authors would like to acknowledge the financial support of the U.S. National Science Foundation (NSF) under Award No. CMS-9900196. The NSF program director was Dr. Ken Chong.     

Thermal fracture analysis of a functionally graded orthotropic strip with a crack

This paper is concerned with the thermal fracture problem of a functionally graded orthotropic strip, where the crack is situated parallel to the free edges. All the material properties are assumed to be dependent only on the coordinate y (perpendicular to the crack surfaces). By using Fourier transform, the thermoelastic problem is reduced to those that involve a system of singular integral equations. Numerical results are presented to show the effects of the crack position and the material distribution on the thermal stress intensity factors.     

An accurate scheme for mixed‐mode fracture analysis of functionally graded materials using the interaction integral and micromechanics models

The interaction integral is a conservation integral that relies on two admissible mechanical states for evaluating mixed‐mode stress intensity factors (SIFs). The present paper extends this integral to functionally graded materials in which the material properties are determined by means of either continuum functions (e.g. exponentially graded materials) or micromechanics models (e.g. self‐consistent, Mori–Tanaka, or three‐phase model). In the latter case, there is no closed‐form expression for the material‐property variation, and thus several quantities, such as the explicit derivative of the strain energy density, need to be evaluated numerically (this leads to several implications in the numerical implementation). The SIFs are determined using conservation integrals involving known auxiliary solutions. The choice of such auxiliary fields and their implications on the solution procedure are discussed in detail. The computational implementation is done using the finite element method and thus the interaction energy contour integral is converted to an equivalent domain integral over a finite region surrounding the crack tip. Several examples are given which show that the proposed method is convenient, accurate, and computationally efficient. Copyright © 2003 John Wiley & Sons, Ltd.     

Interaction integrals for thermal fracture of functionally graded piezoelectric materials

This paper presents domain form of the interaction integrals based on three independent formulations for computation of the stress intensity factors and electric displacement intensity factor for cracks in functionally graded piezoelectric materials subjected to steady-state thermal loading. Each of the formulation differs in the way auxiliary fields are imposed in the evaluation of interaction integral and each of them results in a consistent form of the interaction integral in the sense that extra terms naturally appear in their derivation to compensate for the difference in the chosen crack tip asymptotic fields of homogeneous and functionally graded piezoelectric medium.     

Interaction integrals for thermal fracture of functionally graded materials

This paper addresses finite element evaluation of the non-singular T-stress and mixed-mode stress intensity factors in functionally graded materials (FGMs) under steady-state thermal loads by means of interaction integral. Interaction integral provides an accurate and efficient numerical framework in evaluating these fracture parameters in FGMs under thermal as well as mechanical loads. We use a non-equilibrium formulation and the corresponding auxiliary (secondary) fields tailored for FGMs. Graded finite elements have been developed to account for the spatial gradation of thermomechanical properties. This paper presents various numerical examples in which the accuracy of the present method is verified.     

Interaction integrals for fracture analysis of functionally graded magnetoelectroelastic materials

This paper presents the domain form of interaction integrals based on three independent formulations for computation of stress intensity factors, electric displacement intensity factors and magnetic induction intensity factors for cracks in functionally graded magnetoelectroelastic materials. Conservation integrals of J-type are derived based on the governing equations for magnetoelectroelastic media and the crack tip asymptotic fields of homogeneous magnetoelectroelastic medium as auxiliary fields. Each of the formulations differs in the way auxiliary fields are imposed in the evaluation of interaction integrals and each of them results in a consistent form of the interaction integral in the sense that extra terms naturally appear in their derivation to compensate for the difference in the chosen crack tip asymptotic fields of homogeneous and functionally graded magnetoelectroelastic medium. The additional terms play an important role of ensuring domain independence of the presented interaction integrals. Comparison of numerically evaluated intensity factors through the three consistent formulations with those obtained using displacement extrapolation method is presented by means of two examples.     

Crack analysis in unidirectionally and bidirectionally functionally graded materials

Mixed-mode crack analysis in unidirectionally and bidirectionally functionally graded materials is performed by using a boundary integral equation method. To make the analysis tractable, the Young's modulus of the functionally graded materials is assumed to be exponentially dependent on spatial variables, while the Poisson's ratio is assumed to be constant. The corresponding boundary value problem is formulated as a set of hypersingular traction boundary integral equations, which are solved numerically by using a Galerkin method. The present method is especially suited for straight cracks in infinite FGMs. Numerical results for the elastostatic stress intensity factors are presented and discussed. Special attention of the analysis is devoted to investigate the effects of the material gradients and the crack orientation on the elastostatic stress intensity factors.     

Scaled boundary polygons with application to fracture analysis of functionally graded materials

This study presents a framework for the development of polygon elements based on the scaled boundary FEM. Underpinning this study is the development of generalized scaled boundary shape functions valid for any n‐sided polygon. These shape functions are continuous inside each polygon and across adjacent polygons. For uncracked polygons, the shape functions are linearly complete. For cracked polygons, the shape functions reproduce the square‐root singularity and the higher‐order terms in the Williams eigenfunction expansion. This allows the singular stress field in the vicinity of the crack tip to be represented accurately. Using these shape functions, a novel‐scaled boundary polygon formulation that captures the heterogeneous material response observed in functionally graded materials is developed. The stiffness matrix in each polygon is derived from the principle of virtual work using the scaled boundary shape functions. The material heterogeneity is approximated in each polygon by a polynomial surface in scaled boundary coordinates. The intrinsic properties of the scaled boundary shape functions enable accurate computation of stress intensity factors in cracked functionally graded materials directly from their definitions. The new formulation is validated, and its salient features are demonstrated, using five numerical benchmarks. Copyright © 2014 John Wiley & Sons, Ltd.     

3D crack analysis in functionally graded materials

Elastostatic crack analysis in three-dimensional, continuously non-homogeneous, isotropic and linear elastic functionally graded materials and structures is presented in this paper. A boundary-domain-integral equation formulation is applied for this purpose, which uses the elastostatic fundamental solutions for homogeneous, isotropic and linear elastic materials and involves a domain-integral due to the material’s non-homogeneity. To avoid displacement gradients in the domain-integral, normalized displacements are introduced. The domain-integral is transformed into boundary-integrals over the global boundary of the cracked solids by using the radial integration method. A meshless scheme is developed, which requires only the conventional boundary discretization and additional interior nodes instead of interior cells or meshes. Numerical examples for three-dimensional crack problems in continuously non-homogeneous, isotropic and linear elastic FGMs are presented and discussed, to show the effects of the material gradation on the crack-opening-displacements and the stress intensity factors.     

Stress intensity factors for functionally graded solid cylinders

This paper presents the mode I stress intensity factors for functionally graded solid cylinders with an embedded penny-shaped crack or an external circumferential crack. The solid cylinders are assumed under remote uniform tension. The multiple isoparametric finite element method is used. Various types of functionally graded materials and different gradient compositions for each type are investigated. The results show that the material property distribution has a quite considerable influence on the stress intensity factors. The influence for embedded cracks is quite different from that for external cracks.     

T-stress in orthotropic functionally graded materials: Lekhnitskii and Stroh formalisms

A new interaction integral formulation is developed for evaluating the elastic T-stress for mixed-mode crack problems with arbitrarily oriented straight or curved cracks in orthotropic nonhomogeneous materials. The development includes both the Lekhnitskii and Stroh formalisms. The former is physical and relatively simple, and the latter is mathematically elegant. The gradation of orthotropic material properties is integrated into the element stiffness matrix using a “generalized isoparametric formulation” and (special) graded elements. The specific types of material gradation considered include exponential and hyperbolic-tangent functions, but micromechanics models can also be considered within the scope of the present formulation. This paper investigates several fracture problems to validate the proposed method and also provides numerical solutions, which can be used as benchmark results (e.g. investigation of fracture specimens). The accuracy of results is verified by comparison with analytical solutions.     

An extended element free Galerkin method for fracture analysis of functionally graded materials

An extended element free Galerkin method (XEFGM) has been adopted for fracture analysis of functionally graded materials (FGMs). Orthotropic enrichments functions are used along with the sub-triangle technique for enhancing the Gauss quadrature accuracy near the crack, and the incompatible interaction integral method is employed to calculate the stress intensity factors. Numerical simulations have proved that XEFGM provides more accurate results by less number of nodes (DOFs) in comparison with the unenriched EFGM and other conventional methods for several FGM problems with different crack locations and loadings. The results have been compared with the reference results, showing the reliability, stability, and efficiency of present XEFGM.

Received 9 June 2014 Accepted 17 September 2014.     

A continuum shape sensitivity method for fracture analysis of isotropic functionally graded materials

This paper presents a new continuum shape sensitivity method for calculating mixed-mode stress-intensity factors for a stationary crack in two-dimensional, linear- elastic, isotropic FGMs with arbitrary geometry. The method involves the material derivative concept taken from continuum mechanics, the mutual potential energy release rate, and direct differentiation. Since the governing variational equation is differentiated prior to discretization, resulting sensitivity equations are independent of approximate numerical techniques, such as the finite element method, boundary element method, mesh-free method, or others. The discrete form of the mutual potential energy release rate is simple and easy to calculate, as it only requires multiplication of displacement vectors and stiffness sensitivity matrices. By judiciously selecting the velocity field, the method only requires displacement response in a subdomain close to the crack tip, thus making the method computationally efficient. Seven finite-element based numerical examples, which comprise mode-I and mixed-mode deformations and/or single or multiple interacting cracks, are presented to evaluate the accuracy of the fracture parameters calculated by the proposed method. Comparisons have been made between stress-intensity factors predicted by the proposed method and available reference solutions in the literature, generated either analytically or numerically using various other fracture integrals or analyses. Excellent agreement is obtained between the results of the proposed method and previously obtained solutions. Therefore, shape sensitivity analysis provides an attractive alternative to fracture analysis of cracks in homogeneous and non-homogeneous materials.     

Cohesive modeling of dynamic fracture in functionally graded materials

The dynamic fracture of functionally graded materials (FGMs) is modeled using an explicit cohesive volumetric finite element scheme that incorporates spatially varying constitutive and failure properties. The cohesive element response is described by a rate-independent bilinear cohesive failure model between the cohesive traction acting along the cohesive zone and the associated crack opening displacement. A detailed convergence analysis is conducted to quantify the effect of the material gradient on the ability of the numerical scheme to capture elastodynamic wave propagation. To validate the numerical scheme, we simulate dynamic fracture experiments performed on model FGM compact tension specimens made of a polyester resin with varying amounts of plasticizer. The cohesive finite element scheme is then used in a parametric study of mode I dynamic failure of a Ti/TiB FGM, with special emphasis on the effect of the material gradient on the initiation, propagation and arrest of the crack.     

功能梯度材料薄球壳热屈曲问题

研究了由金属和陶瓷组成的功能梯度材料薄球壳热屈曲问题。用张量方法推导得到轴对称球壳稳定性方程。将热本构方程应用到球壳稳定性方程中,得到以位移表示的球壳热屈曲方程组。分别考虑均布外压和温度作用,采用伽辽金法计算分析简支球壳的热屈曲问题,给出薄球壳厚度、物性参数变化、内外表面温差变化引起的临界温度变化趋势和临界压力变化趋势。     

考虑尺度依赖的平面正交各向异性功能梯度微梁的自由振动分析

基于一种新修正偶应力理论建立了微尺度平面正交各向异性功能梯度梁的自由振动模型。模型中包含两个材料尺度参数,能够分别描述两个正交方向上不同程度的尺度效应。当梁的几何尺寸远大于材料尺度参数时,本文模型亦可自动退化为相应的传统宏观模型。基于哈密顿原理推导了运动控制方程并以简支梁的自由振动为例分析了几何尺寸、功能梯度变化指数等对尺度效应产生的影响。算例结果表明:采用本文模型所预测的梁自振频率总是大于传统理论的结果,即捕捉到了尺度效应。尺度效应会随着梁几何尺寸的增大而逐渐减弱并在几何尺寸远大于尺度参数时消失;高阶自振频率所体现出的尺度效应较低阶自振频率更加明显。此外,功能梯度变化指数对尺度效应也有一定的影响。     

Scattering of the harmonic anti-plane shear waves by a crack in functionally graded piezoelectric materials

In the present paper, the dynamic behavior of a Griffth crack in the functionally graded piezoelectric material (FGPM) is investigated. It is assumed that the elastic stiffness, piezoelectric constant, dielectric permittivity and mass density of the FGPM vary continuously as an exponential function, and that FGPM is under the anti-plane mechanical loading and in-plane electrical loading. By using the Fourier transform and defining the jumps of displacement and electric potential 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 and electric potential components across the crack surface are expanded in a series of Jacobi polynomial. Numerical examples are provided to show the effects of material properties on the stress and the electric displacement intensity factors.     

Simulation of crack paths in functionally graded materials

One of the main interests of fracture mechanics in functionally graded materials is the influence of such an inhomogeneity on crack propagation processes. Using the Griffith’ energy principle, the change of energy has to be calculated, if the crack starts to propagate. In homogeneous linear-elastic structures (asymptotically precise) formulas for the energy release rate are known, but a direct transfer of these methods to functionally graded materials can lead to very inaccurate results. Moreover, the influence of the inhomogeneity on the crack path cannot be seen. Here, a simple model for functionally graded materials is introduced. For this model, a formula for the change of potential energy is derived, giving detailed information on the effect of the gradation on crack propagation.     

Antiplane crack analysis of a functionally graded material by a BIEM

Elastostatic analysis of an antiplane crack in a functionally graded material (FGM) is performed by using a hypersingular boundary integral equation method (BIEM). An exponential law is applied to describe the spatial variation of the shear modulus of the FGM. A Galerkin method is applied for the numerical solution of the hypersingular traction BIE. Both unidirectional and bidirectional material gradations are investigated. Stress intensity factors for an infinite and linear elastic FGM containing a finite crack subjected to an antiplane crack-face loading are presented and discussed. The influences of the material gradients and the crack orientation on the stress intensity factors are analyzed.     

Finite element simulations of crack propagation in functionally graded materials under flexural loading

Cracks in stepped and continuously graded material specimens under flexural loading were investigated via finite element analysis. Calculation of mechanical energy release rates and propagation angles with crack-opening displacement correlation and the local symmetry (K_II = 0) criterion, respectively, provided results most efficiently and accurately, as compared with compliance and J-integral approaches and other deflection criteria. A routine was developed for automatic crack extension and remeshing, enabling simulation of incremental crack propagation. Effects of gradient profile and crack geometry on crack-tip stresses and crack propagation path are examined, and implications of these for optimal design of graded components against failure by fast fracture are discussed.