Simulation of Crack Propagation in Functionally Graded Materials Under Mixed-Mode and Non-Proportional Loading

Automatic simulation of crack propagation in homogeneous and functionally graded materials is performed by means of a remeshing algorithm in conjunction with the finite element method. The crack propagation is performed under mixed-mode and non-proportional loading. Each step of crack growth simulation consists of calculation of mixed-mode stress intensity factors by means of a novel formulation of the interaction integral method, determination of crack growth direction based on a specific fracture criterion, and local automatic remeshing along the crack path. The present approach requires a user-defined crack increment at the beginning of the simulation. Crack trajectories obtained by the present numerical simulation are compared with available experimental results.     

Finite element evaluation of mixed mode stress intensity factors in functionally graded materials

This paper is directed towards finite element computation of fracture parameters in functionally graded material (FGM) assemblages of arbitrary geometry with stationary cracks. Graded finite elements are developed where the elastic moduli are smooth functions of spatial co‐ordinates which are integrated into the element stiffness matrix. In particular, stress intensity factors for mode I and mixed‐mode two‐dimensional problems are evaluated and compared through three different approaches tailored for FGMs: path‐independent J^*_k‐integral, modified crack‐closure integral method, and displacement correlation technique. The accuracy of these methods is discussed based on comparison with available theoretical, experimental or numerical solutions. Copyright © 2001 John Wiley & Sons, Ltd.     

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.     

Three-dimensional static and dynamic analyses of the functionally graded cylinder bonded to the laminated plate under general loading

A graded finite element method code based on Rayleigh–Ritz energy formulation is developed and implemented to study the elastic behavior of a layered plate loaded by a solid isotropic cylinder and a functionally graded interlayer. The applied nonaxisymmetric loading to the inner cylinder induces a stress concentration in the flexible part of the joint. The effects of different thicknesses and power law exponents of functionally graded interlayer on the distribution of displacements and stresses are investigated, which verifies the ability of functionally graded material to control the stress and displacement waves. The time-dependent response of the structure is also obtained based on Newmark's time integration method.     

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.     

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.     

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.     

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.     

Nonlinear vibration and instability analysis of functionally graded CNT-reinforced cylindrical shells conveying viscous fluid resting on orthotropic Pasternak medium

In this study, nonlinear vibration and instability of embedded temperature-dependent cylindrical shell conveying viscous fluid resting on temperature-dependent orthotropic Pasternak medium are investigated. The equivalent material properties of nanocomposites are estimated using rule of mixture. Both cases of uniform distribution and functionally graded distribution patterns of reinforcements are considered. Based on orthotropic Mindlin shell theory, the governing equations are derived. Generalized differential quadrature method is applied for obtaining the frequency and critical fluid velocity of a system. The effects of different parameters, such as distribution type of single-walled carbon nanotubes (SWCNTs), volume fractions of SWCNTs, and Pasternak medium are discussed.     

Fracture problems,vibration, buckling,and bending analyses of functionally graded materials: A state-of-the-art review including smart FGMS

Composite materials fail under extreme working conditions, particularly at high temperature, due to delamination (separation of fibers from matrix). And therefore it is needed to switch over functionally graded materials (FGMs) which can sustain at high temperature conditions (250–2000°C). There is a need to analyze the fracture and fatigue characteristics of FGM structures and so through this review the emphasis is given on fracture analysis of FGM materials. It has been reported that a combination of extended finite element method and isogeometric analysis methodologies has been used for general mixed-mode crack propagation problems after the introduction of extended isogeometric analysis. Furthermore, recent computational advances have been in the form of multiscale simulations where the part of model is simulated by a finer modeling scale, which can represent details of the material behavior and the interacting effects of material constituents in the finest way. The review is also focused on new advances in analytical and numerical methods for the stress, vibration, and buckling analyses of FGMs. Emphasis has been primarily on to restrict 2D analysis with sorts of compromise in the accuracy of results. First shear deformation theory (FSDT) and third-order shear deformation theory have been extensively used among the various 2D plate theories. FSDT can help us in terms of getting reasonably accurate results with less computational afford. This paper also outlines review on carbon nanotubes (CNT) reinforced FGMs, functionally graded nanocomposites, functionally graded single-walled CNT, FG nanobeam as well as functionally graded piezoelectric materials. Future applications would be based on these smart materials which are supposed to serve us in adverse conditions. Of course, with rise and advent of promising nanotechnology and its potential impact on aerospace industry as well as on other areas, it becomes important to us to compile this review article.     

圆板状SiC/C功能梯度材料残余热应力特征有限元分析

功能梯度材料残余热应力的大小及分布对其性能有效发挥及长期稳定使用有着较大的负面影响,为了尽可能充分发挥材料性能,增加材料的使用寿命,需尽可能减小残余应力以及使其合理分布.本文采用ANSYS有限元分析软件对不同叠层工艺参数的等离子体第一壁候选材料--SiC/C功能梯度材料(FGM)的残余热应力进行了数值模拟,获得了使热应力有效缓和的较适宜的工艺参数,对实际研发制备目标材料也可提供一些理论参照.相关结果表明,适量增加梯度叠层数及中间梯度层厚度可逐步有效缓和残余热应力,同时,针对本文今后应用的仍以炭材料为主体的炭基陶瓷保护层复合SiC/C FGM而言,纯SiC层厚度应取较小值,而叠层成分分布指数应取0.8～1.0为宜.     

Cohesive fracture modeling of elastic-plastic crack growth in functionally graded materials

This work investigates elastic-plastic crack growth in ceramic/metal functionally graded materials (FGMs). The study employs a phenomenological, cohesive zone model proposed by the authors and simulates crack growth by the gradual degradation of cohesive surfaces ahead of the crack front. The cohesive zone model uses six material-dependent parameters (the cohesive energy densities and the peak cohesive tractions of the ceramic and metal phases, respectively, and two cohesive gradation parameters) to describe the constitutive response of the material in the cohesive zone. A volume fraction based, elastic-plastic model (extension of the original Tamura-Tomota-Ozawa model) describes the elastic-plastic response of the bulk background material. The numerical analyses are performed using WARP3D, a fracture mechanics research finite element code, which incorporates solid elements with graded elastic and plastic properties and interface-cohesive elements coupled with the functionally graded cohesive zone model. Numerical values of volume fractions for the constituents specified at nodes of the finite element model set the spatial gradation of material properties with isoparametric interpolations inside interface elements and background solid elements to define pointwise material property values. The paper describes applications of the cohesive zone model and the computational scheme to analyze crack growth in a single-edge notch bend, SE(B), specimen made of a TiB/Ti FGM. Cohesive parameters are calibrated using the experimentally measured load versus average crack extension (across the thickness) responses of both Ti metal and TiB/Ti FGM SE(B) specimens. The numerical results show that with the calibrated cohesive gradation parameters for the TiB/Ti system, the load to cause crack extension in the FGM is much smaller than that for the metal. However, the crack initiation load for the TiB/Ti FGM with reduced cohesive gradation parameters (which may be achieved under different manufacturing conditions) could compare to that for the metal. Crack growth responses vary strongly with values of the exponent describing the volume fraction profile for the metal. The investigation also shows significant crack tunneling in the Ti metal SE(B) specimen. For the TiB/Ti FGM system, however, crack tunneling is pronounced only for a metal-rich specimen with relatively smaller cohesive gradation parameter for the metal.     

Thermal fracture analysis of orthotropic functionally graded materials using an equivalent domain integral approach

A new computational method based on the equivalent domain integral (EDI) is developed for mode I fracture analysis of orthotropic functionally graded materials (FGMs) subjected to thermal stresses. By using the constitutive relations of plane orthotropic thermoelasticity, generalized definition of the J-integral is converted to an equivalent domain integral to calculate the thermal stress intensity factor. In the formulation of the EDI approach, all the required thermomechanical properties are assumed to have continuous spatial variations through the functionally graded medium. Developed methodology is integrated into a fracture mechanics research finite element code FRAC2D using graded finite elements that possess cubic interpolation. Steady-state and transient temperature distribution profiles in orthotropic FGMs are computed using the finite elements based heat transfer analysis software HEAT2D. EDI method is validated and domain independence is demonstrated by comparing the numerical results obtained using EDI to those calculated by an enriched finite element method and to those available in the literature. Single and periodic edge crack problems in orthotropic FGMs are examined in order to study the influences of principal thermal expansion coefficient and thermal conductivity components, relative crack length and crack periodicity on the thermal stress intensity factors. Numerical results show that among the three principal thermal expansion coefficient components, the in-plane component perpendicular to the crack axis has the most significant influence on the mode I stress intensity factor. Gradation profile of the thermal expansion coefficient parallel to the crack axis is shown to have no effect on the outcome of the fracture analysis.     

A complex variable approach for asymptotic Mode-III crack-tip fields in an anisotropic functionally graded material

This paper addresses asymptotic full crack-tip fields for an anti-plane (Mode-III) stationary crack in an anisotropic functionally graded material. A monoclinic material that has a material symmetry plane is considered. The complex variable approach and the asymptotic analysis are used to solve a perturbed Laplace equation resulting from material anisotropy and gradation. The out-of-plane displacement and stress solutions are provided for a crack in exponentially and linearly graded anisotropic materials by considering material gradation either parallel or perpendicular to the crack. The characteristics of the asymptotic solutions in an anisotropic functionally graded material are compared with those for anisotropic homogeneous and isotropic graded materials. Finally, engineering significance of the present work is discussed.     

New implementation of MLBIE method for heat conduction analysis in functionally graded materials

The meshless local boundary integral equation (MLBIE) method with an efficient technique to deal with the time variable are presented in this article to analyze the transient heat conduction in continuously nonhomogeneous functionally graded materials (FGMs). In space, the method is based on the local boundary integral equations and the moving least squares (MLS) approximation of the temperature and heat flux. In time, again the MLS approximates the equivalent Volterra integral equation derived from the heat conduction problem. It means that, the MLS is used for approximation in both time and space domains, and we avoid using the finite difference discretization or Laplace transform methods to overcome the time variable. Finally the method leads to a single generalized Sylvester equation rather than some (many) linear systems of equations. The method is computationally attractive, which is shown in couple of numerical examples for a finite strip and a hollow cylinder with an exponential spatial variation of material parameters.     

Numerical and experimental study of molten pool behaviour and defect formation in multi-material and functionally graded materials laser powder bed fusion

Laser powder bed fusion (LPBF) of multi-material and functionally graded materials (FGM) has attracted significant research interest due to its ability to fabricate components with superior performance compared with those manufactured with single powder material. However, the forming mechanisms of various defects remain unknown. In this paper, a DEM-CFD model was first established to obtain an in-depth understanding of this process. It was discovered that the defects including partially melted and un-melted Invar36 powder were embedded in the lower level of the powder layer; this was attributed to the low laser absorptivity, low melting point and high thermal conductivity of the Cu10Sn powder. Inter-layer defects were more likely to occur with an increased powder layer thickness. In addition, the scanned track width was found related to an equilibrium achieved among the thermal properties of the powder mixture. Process parameters were optimised to obtain FGM structures without defects in both horizontal and vertical directions. Invar36/Cu10Sn samples were fabricated with a multi-material LPBF system using different mixed powder contents and laser volumetric energy densities (VEDs). By increasing the VED, fewer defects were observed between the interface of two processed powder layers, which had a good agreement with the modelling results.     

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.     

Four-parameter functionally graded cracked plates of arbitrary shape: A GDQFEM solution for free vibrations

The dynamic behavior of moderately thick FGM plates with geometric discontinuities and arbitrarily curved boundaries is investigated. The Generalized Differential Quadrature Finite Element Method (GDQFEM) is proposed as a numerical approach. The irregular physical domain in Cartesian coordinates is transformed into a regular domain in natural coordinates. Several types of cracked FGM plates are investigated. It appears that GDQFEM is analogous to the well-known Finite Element Method (FEM). With reference to the proposed technique the governing FSDT equations are solved in their strong form and the connections between the elements are imposed with the inter-element compatibility conditions. The results show excellent agreement with other numerical solutions obtained by FEM.     

Numerical simulation of elastic plastic fatigue crack growth in functionally graded material using the extended finite element method

In the present work, the extended finite element method has been used to simulate the fatigue crack growth problems in functionally graded material in the presence of holes, inclusions, and minor cracks under plastic and plane stress conditions for both edge and center cracks. Both soft and hard inclusions have been implemented in the problems. The validity of linear elastic fracture mechanics theory is limited to the brittle materials. Therefore, the elastic plastic fracture mechanics theory needs to be utilized to characterize the plastic behavior of the material. A generalized Ramberg-Osgood material model has been used for modeling purposes.     

四边固支功能梯度板中波的传播

摘要：考虑剪切变形和转动惯性的影响,采用一阶剪切变形板理论和小应变的应变－位移关系,利用Hamilton原理推得运动控制方程,并应用特征值方程得到频散方程。给出了波在功能梯度板中传播的频散,相速度和群速度随波数变化的曲线,分析了不同的功能梯度材料指数对波传播的影响规律。