梯度复合材料裂纹扩展路径和起裂载荷的有限元分析

为了模拟功能梯度材料(FGM)在工程应用中可能会出现的断裂问题并计算相应的开裂载荷,通过编写用户自定义UEL子程序将梯度扩展单元嵌入到ABAQUS软件中模拟功能梯度材料的物理场,并编写交互能量积分后处理子程序计算裂纹尖端的混合模式应力强度因子(SIF),采用最大周向应力准则编写子程序计算裂纹的偏转角,并模拟了裂纹扩展路径,计算了裂纹的起裂载荷。讨论了材料梯度参数对裂纹扩展路径以及起裂载荷的影响规律。通过与均匀材料的对比,验证了功能梯度材料断裂性能的优越性。研究表明:外载平行于梯度方向时,垂直梯度方向的初始裂纹朝着等效弹性模量小的方向扩展,且偏转角在梯度指数线性时出现峰值,并随着组分弹性模量比的增加而变大;当外载和初始裂纹均平行于梯度方向时,材料等效弹性模量和断裂韧性的增加或者梯度指数的减小都导致起裂载荷变大。     

The collinear cracks in an inhomogeneous medium subjected to transient load

Summary A laminate model is employed to solve the elastodynamic problem of a collinear crack in an inhomogeneous material. The inhomogeneous material is treated as a series of thinner layer. The Laplace and Fourier transforms are used to reduce the problem to a set of singular integral equations that is solved numerically. Numerical results of two collinear cracks in a functionally graded material strip are obtained to show the influence of material inhomogeneity and crack position on crack tip field intensities.     

Curved crack propagation in homogeneous and graded materials

Deflection and deviation of cracks commonly occurs because of asymmetry in crack‐tip stresses in both homogeneous materials and functionally graded materials (FGMs); yet the analysis of curved cracks has been limited to simple crack shapes, otherwise the analysis would involve extensive levels of computation. The present study investigates the approximation of curved cracks with simplified shapes. A simple analytical model justifying the use of crack‐shape approximations, developed in an earlier study on stationary curved cracks in homogeneous materials, is outlined. Then, the approach is applied to propagating cracks in both homogeneous and graded material structures. Results are presented from finite element (FE) simulations of crack propagation using exact and simplified crack shapes. The use of an approximated crack shape can provide basic estimates for crack propagation path and critical load. However, systematic divergence can occur between predictions for exact and approximated crack shapes, particularly in inhomogeneous material configurations, and so the development of solutions for non‐straight cracks in FGMs would be expedient.     

Some generalized anti-plane strain problems for an inhomogeneous elastic halfspace

This paper develops a generalized formulation of the theory of anti-plane deformations of a linear elastic solid, the Lamé constants of which depend on a single spatial variable. The generalized theory is applied to the study of a generalized screw dislocation, a semi-infinite crack and a finite crack located in a halfspace region where the elastic inhomogeneity is depth dependent. The problems examined in this paper are relevant to the modelling of cracks located in inhomogeneous geological materials and the study of surface defects associated with functionally graded materials.     

功能梯度压电双材料板中厚度-扭曲波的传播

分析了厚度-扭曲波在无限大功能梯度压电双材料板中的传播性能,板的上表面和下表面是机械自由和电学开路的,材料常数在厚度方向按指数规律变化。首先推导了满足控制方程和边界条件的电弹场,然后利用界面条件得到了厚度-扭曲波传播所应满足的关系。通过算例表明了材料梯度变化对厚度-扭曲波传播性能的影响,结果对功能梯度压电材料在声波器件中应用有参考价值。     

A Review of Dynamic Fracture Studies in Functionally Graded Materials

Abstract:  This article presents a review of dynamic fracture studies on functionally graded materials. A brief literature review on the fracture mechanics of graded materials is presented first. This is followed by a discussion on the higher-order asymptotic analysis of the transient elastic field surrounding the tip of a dynamically growing crack in a functionally graded material. A comprehensive experimental study of dynamic crack growth in model functionally graded material using the optical method of reflection photoelasticity and high-speed photography is then presented. The results are analysed to establish a generalised relationship between the crack velocity and the dynamic mode-I stress intensity factor (SIF). This relationship is found to be unique and is distinctly different from that previously established for the matrix material (polyester). Finally, an innovative experimental procedure is used to demonstrate the necessity of employing a fully transient stress-field representation in the analysis of optical data for an accurate prediction of the dynamic SIF history.     

Precise computation and error control of stress intensity factors and certain integral characteristics in anisotropic inhomogeneous materials

The numerical simulation of quasi-static crack propagation is closely related to the computation of characteristics such as stress intensity factors or energy release rates. In this work, ideas are proposed, how such quantities can be calculated precisely in linear elastic, anisotropic and inhomogeneous plane structures. Stress intensity factors and other local characteristics can be evaluated in terms of functionals depending on solutions of certain elasticity problems. The approach used here to calculate these functional values precisely with Galerkin finite elements is a dual-weighted-residual method for adaptive mesh refinement and a posteriori error control. Especially in structures under mixed-mode loadings contact of crack faces can occur. A numerical realization of mutual non-penetration conditions of inequality type for the crack faces and the effect of such constraints on stress intensity factors is shown. Numerical results are presented for anisotropic and functionally graded materials.     

Fracture behaviors for a functionally graded superconducting cylinder with a cylindrical crack

In this study, a double exponential model is proposed to investigate the cylindrical crack problem for a functionally graded superconducting cylinder. The stress intensity factors (SIFs) are analytically obtained by transforming the corresponding crack problem into dual integral equations. The effects of applied magnetic field, model parameters, and crack configuration on the SIFs are analyzed. Some important phenomena are observed. Among others, both decreasing the graded index of Young's modulus and increasing the introduced nondimensional exponent parameter in the critical current model can inhibit crack propagation. This study should be useful for the application of superconducting devices.     

Wave propagation and transient response of a FGM plate under a point impact load based on higher-order shear deformation theory

In this paper, the wave propagation and transient response of an infinite functionally graded plate under a point impact load are presented. The effective material properties of functionally graded materials (FGMs) for the plate are assumed to vary continuously through the plate thickness and be distributed according to a volume fraction power law along the plate thickness. Based on the higher-order shear deformation theory and considering the effect of the rotary inertia, the governing equations of the wave propagation in the functionally graded plate are derived by using the Hamilton’s principle. The analytic dispersion relation of the functionally graded plate is obtained by means of integral transforms and a complete discussion of dispersion for the functionally graded plate is given. Then, using the dispersion relation and integral transforms, exact integral solutions for the functionally graded plate under a point impact load are obtained. The transient response curves of the functionally graded plates are plotted and the influence of volume fraction distributions on transient response of functionally graded plates is analyzed. Finally, the solutions of the higher-order shear deformation theory and the first-order shear deformation theory are studied.     

Crack Tip Shielding or Anti-shielding due to Smooth and Discontinuous Material Inhomogeneities

This paper describes a theoretical model and related computational methods for examining the influence of inhomogeneous material properties on the crack driving force in elastic and elastic-plastic materials. Following the configurational forces approach, the crack tip shielding or anti-shielding due to smooth (e.g. graded layer) and discontinuous (e.g. bimaterial interface) distributions in material properties are derived. Computational post-processing methods are described to evaluate these inhomogeneity effects. The utility of the theoretical model and computational methods is demonstrated by examining a bimaterial interface perpendicular to a crack in elastic and elastic-plastic compact tension specimens.     

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.     

Periodically distributed parallel cracks in a functionally graded piezoelectric (FGP) strip bonded to a FGP substrate under static electromechanical load

In this paper, the Fourier integral transform–singular integral equation method is presented for the problem of a periodic array of cracks in a functionally graded piezoelectric strip bonded to a different functionally graded piezoelectric material. The properties of two materials, such as elastic modulus, piezoelectric constant and dielectric constant, are assumed in exponential forms and vary along the crack direction. The crack surface condition is assumed to be electrically impermeable or permeable. The mixed boundary value problem is reduced to a singular integral equation over crack by applying the Fourier transform and the singular integral equation is solved numerically by using the Lobatto–Chebyshev integration technique. The analytic expressions of the stress intensity factors and the electric displacement intensity factors are derived. The effects of the loading parameter λ, material constants and the geometry parameters on the stress intensity factor, the energy release ratio and the energy density factor are studied.     

Estimation of the Crack Length After Thermal Shock in FGM Strip

Thermal shock in a strip with a functionally graded Al₂O₃/Al interpenetrating network microstructure is modelled. The problem of thermal shock induced crack propagation is analyzed. A strong influence of residual stresses and the type of composition gradient on the crack length is predicted. The type of property gradient required in order to achieve an optimal thermal shock resistance is discussed.     

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.     

Fracture mechanics analysis models for functionally graded materials with arbitrarily distributed properties (Modes II and III problems)

In Wang et al. (2002), the mode I crack problem has been analyzed for a functionally graded material strip with arbitrarily distributed properties. For the completeness of the research, this paper provides an analysis method for the mode II in plane and mode III anti-plane elastic crack problems for a functionally graded material (FGM) strip containing a crack along the gradient direction. In the analysis, the gradient region is treated as a number of homogeneous layers stacked along the gradient direction. The analytical models developed in this paper for the modes II and III problems, together with the model developed in Wang et al. (2002) for the mode I crack problem complete the analysis of the mixed-mode crack problem for functionally graded materials with a crack parallel the gradient direction. Those analytical models allow the material properties be any continuous functions along the thickness direction of the strip.     

An analytical approach for fracture and fatigue in functionally graded materials

The problem of brittle crack propagation and fatigue crack growth in functionally graded materials (FGMs) is addressed. The proposed analytical approach can be used to estimate the variation of the stress-intensity factor as a function of the crack length in FGMs. Furthermore, according to the Paris’ law, the fatigue life and the crack-tip velocity of crack propagation can be predicted in the case of fatigue crack growth. A comparison with numerical results obtained according to the Finite Element method will show the effectiveness of the proposed approach. Detailed examples are provided in the case of three-point bending beam problems with either a FGM interlayer, or a FGM external coating. A comparison is presented between two types of grading in the elastic modulus: a continuous linear variation in the FGM layer and a discrete approximation with a multi-layered beam and a constant Young’s modulus in each layer.     

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.     

A thermomechanical fracture modeling and simulation for functionally graded solids using a residual-strain formulation

This paper addresses mixed-mode crack growth in two-dimensional functionally graded solids under thermomechanical loads, and investigates the effect of mechanical and thermal loads as well as the T-stress on their crack growth behavior. A novel residual strain-based formulation in the interaction integral method is developed and used for the accurate evaluation of mixed-mode stress intensity factors and/or the T-stress. Simulation of mixed-mode crack propagation in functionally graded materials including solid oxide fuel cells under thermomechanical loads is performed by means of the finite element method and the generalized interaction integrals in conjunction with a remeshing algorithm. An iterative procedure is used for crack growth simulation including the calculation of mixed-mode stress intensity factors and/or the T-stress by means of the generalized interaction integral method, determination of crack growth direction and crack initiation condition based on selected fracture criteria, and local automatic remeshing along the crack path. The present approach employs a user-defined crack increment at the beginning of the simulation. Crack trajectories and fracture parameters obtained by the present simulation for thermomechanical loads are assessed for some numerical examples in comparison with those for mechanical loads.     

Crack growth resistance behavior of a functionally graded material: computational studies

This paper describes crack growth resistance simulation in a ceramic/metal functionally graded material (FGM) using a cohesive zone ahead of the crack front. The plasticity in the background (bulk) material follows J₂ flow theory with the flow properties determined by a volume fraction based, elastic-plastic model (extension of the original Tamura-Tomota-Ozawa model). A phenomenological, cohesive zone model with 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) describes the constitutive response of the cohesive zone. Crack growth occurs when the complete separation of the cohesive surfaces takes place. The crack growth resistance of the FGM is characterized by a rising J-integral with crack extension (averaged over the specimen thickness) computed using a domain integral (DI) formulation. The 3-D 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. The paper describes applications of the cohesive zone model and the DI method to compute the J resistance curves for both single-edge notch bend, SE(B), and single-edge notch tension, SE(T), specimens having properties of a TiB/Ti FGM. The numerical results show that the TiB/Ti FGM exhibits significant crack growth resistance behavior when the crack grows from the ceramic-rich region into the metal-rich region. Under these conditions, the J-integral is generally higher than the cohesive energy density at the crack tip even when the background material response remains linearly elastic, which contrasts with the case for homogeneous materials wherein the J-integral equals the cohesive energy density for a quasi-statically growing crack.