功能梯度材料的黏弹性断裂问题

功能梯度材料(FGM)是一种不同于传统复合材料的新型工程复合材料 [1], 国内外关于FGM的断裂力学方面的研究发展非常迅速. 关于FGM静态裂纹问题,学者们研究了不同类型裂纹尖端场的应力强度因子 [2-5], 探讨了有限长裂纹在不用载荷作用下的传播等问题. 而关于动态裂纹问题,也已经取得很大成就 [6-9]. FGM一个很重要的应用是高温结构材料,在强大的热环境中,很多材料都呈现出黏弹性. 因此,研究FGM的黏弹性断裂力学非常具有实际价值.对此,众多研究 [10-14]提出不同的分析模型,并在不同受载条件,通过理论计算,分析了黏弹性裂纹尖端场的力学 行为.本文考查了功能梯度材料板条中界面裂纹垂直于梯度方向时的黏弹性断裂问题,首先利用有限元法求解线弹性功能梯度材料板条的裂纹尖端场,然后根据黏弹性的对应性原理,求解出黏弹性功能梯度材料板条裂纹问题的应力场强度因子.      

平行于功能梯度材料夹层的币型裂纹起裂条件

分析了功能梯度材料中币型裂纹的扩展问题．裂纹平行于无限域中功能梯度材料夹层,受有与裂纹面成任意角度的拉应力．假定功能梯度材料夹层与两个半无限域均匀介质完全粘合,其弹性模量沿厚度方向变化．采用基于层状材料广义Kelvin基本解的边界元方法分析裂纹问题,给出了均布正应力和剪应力作用下裂纹的应力强度因子、将应力强度因子耦合于应变能密度断裂判据,讨论了裂纹体在拉伸应力作用下的起裂条件．     

A micromechanical model for predicting the fracture toughness of functionally graded foams

In this paper, finite element method based micromechanical analysis is used to understand the fracture behavior of functionally graded foams. The finite element analysis uses a micro-mechanical model in conjunction with a macro-mechanical model in order to relate the stress intensity factor to the stresses in the struts of the foam. The stress intensity factor at the crack tip of the macro-mechanical model can be evaluated using either the J-contour integral or the stresses in the singularity-dominated zone. The fracture toughness is evaluated for various crack positions and length within the functionally graded foam. Then the relationship between the fracture toughness of the graded foam and the local density at the crack tip is studied. Convergence tests for both macro-mechanical and micro-mechanical model analysis were conducted in order to maintain adequate accuracy with reasonable computational time. Fracture toughness of homogenous foams and functionally graded foams for various cases are presented as a function of relative density. This study indicates that the fracture toughness of functionally graded foams mainly depends on the relative density at the crack-tip.     

Nonlocal theory solution of two collinear cracks in the functionally graded materials

In this paper, the interaction of two collinear cracks in functionally graded materials subjected to a uniform anti-plane shear loading is investigated by means of nonlocal theory. The traditional concepts of the nonlocal theory are extended to solve the fracture problem of functionally graded materials. To make the analysis tractable, it is assumed that the shear modulus varies exponentially with the coordinate vertical to the crack. By use of the Fourier transform, the problem can be solved with the help of a pair of triple integral equations, in which the unknown variable is the displacement on the crack surfaces. To solve the triple integral equations, the displacement on the crack surfaces is expanded in a series of Jacobi polynomials. Unlike the classical elasticity solutions, it is found that no stress singularity is present near the crack tips. The nonlocal elastic solutions yield a finite hoop stress at the crack tip, thus allowing us to use the maximum stress as a fracture criterion in functionally graded materials. The magnitude of the finite stress field depends on the crack length, the distance between two cracks, the parameter describing the functionally graded materials and the lattice parameter of the materials.     

Mode III fracture problem of a cracked FGPM surface layer bonded to a cracked FGPM substrate

This work deals with the mode III fracture problem of a cracked functionally graded piezoelectric surface layer bonded to a cracked functionally graded piezoelectric substrate. The cracks are normal to the interface and the electro-elastic material properties are assumed to be varied along the crack direction. Potential and flux types of boundary condition are assigned on the edge of the surface layer. The problem under the assumptions of impermeable and permeable cracks can be formulated to the standard singular integral equations, which are solved by using the Gauss–Chebyshev technique. The effects of the boundary conditions, the material properties and crack interaction on the stress and electric displacement intensity factors are discussed.     

黏弹性功能梯度材料裂纹问题的有限元方法

针对组分材料体积含量任意分布的黏弹性功能梯度材料裂纹问题建立有限元分析途径. 通过Laplace变换,将黏弹性问题转化到象空间中求解,基于反映材料非均匀的梯度单元和裂纹尖端奇异特性的奇异单元计算象空间中的位移、应力和应变场,应用虚拟裂纹闭合方法得到应变能释放率,分别由应力和应变能释放率确定应力强度因子. 给出这些断裂参量在物理空间和象空间之间的对应关系,由数值逆变换求出其在物理空间的相应值. 文中分析两端均匀受拉的黏弹性边裂纹板条,首先针对松弛模量表示为空间函数和时间函数乘积的特殊梯度材料进行计算,结合对应原理验证方法的有效性. 然后分析组分材料体积含量具有任意梯度分布的情形,由Mori-Tanaka方法预测象空间中的等效松弛模量. 计算结果表明,蠕变加载条件下,应变能释放率随时间增加,其增大程度与黏弹性组分材料体积含量相关. 由于梯度材料的非均匀黏弹性性质,产生应力重新分布,导致应力强度因子随时间变化,其变化范围与组分材料的体积含量分布方式有关.     

Stress Concentration Factor in Functionally Graded Plates with Circular Holes Subjected to Anti-plane Shear Loading

Stress concentration factors due to the presence of geometrical discontinuities (circular holes) in functionally graded plates are derived. The material property inhomogeneity is assumed to be in the radial direction originating at the center of the plate. Variable separable closed-form solutions are obtained for the stresses and displacements in functionally graded plates (without and with holes) subjected to anti-plane shear loading. The stresses in functionally graded plates without a hole are not homogeneous as it is in the case of homogeneous plates. Either a stress concentration (more than the applied stress) or dilution (less than the applied stress) occurs depending on whether the modulus increases (hardening graded material) or decreases (softening graded material) away from the center of the graded plate without a hole. A novel definition of the stress concentration factor due to the geometrical discontinuity in functionally graded plates is derived. The effect of the circular hole in functionally graded plates is to magnify (compared to homogeneous plates) the stress concentration when the modulus decreases away from the center of the hole (softening material). Beneficial reduction of the stress concentration factor is achieved in hardening functionally graded materials.     

一种特殊梯度分布的功能梯度圆柱壳的二维分析

功能梯度材料是一种新型材料,其结构分析已成为当今力学研究的热点。本文对一种特殊梯度分布的功能梯度材料圆柱壳进行了二维精确分析。从弹性力学平面应变问题的基本方程出发,引入应力函数,导出功能梯度材料圆柱壳受静载作用下的控制微分方程。假设材料的杨氏模量沿半径方向呈幂函数分布,泊松比为常数,利用分离变量法,导出了简支边界情况下功能梯度圆柱壳的精确解。通过算例分析了不同梯度变化时,功能梯度圆柱壳内的应力和位移变化规律。计算结果表明不同梯度分布的圆柱壳结构中的应力、位移沿厚度方向的变化规律是不同的,有时甚至差别很大。因此对于材料性质梯度变化的功能梯度材料圆柱壳,必须针对其自身特点,建立相应的理论分析模型。     

Modeling method for a crack problem of functionally graded materials with arbitrary properties—piecewise-exponential model

A new model, piecewise-exponential model (PE model), is developed to investigate the crack problem of the functionally graded materials (FGMs) with arbitrary properties. In the PE model, the functionally graded material is divided into some nonhomogeneous layers along the gradient direction of the properties, with each layer’s properties varying exponentially. By this way, the real material properties can be approached by a series of exponential functions. Since the real material properties are used on both surfaces of each nonhomogeneous layer, the nature of continuously varying properties of FGMs can be approached accurately. The influences of the local nonhomogeneity on the crack-tip fields can be fully considered. By using the new model, the fracture problem of a functionally graded strip with arbitrary properties and a crack vertical to the free surfaces is studied. The integral transform method, the theory of residues and the theory of singular integral equation are applied. Some representative samples with different kinds of nonhomogeneous properties are analyzed and the corresponding stress intensity factors (SIFs) are presented. It is shown that the PE mode is effective for investigating the crack problems of the FGMs with arbitrary properties.     

Mode-3 spontaneous crack propagation in unsymmetric functionally graded materials

The effects of the inhomogeneous material property variation on the mode-3 crack propagation characteristics have been analyzed. The spectral form of the elastodynamic boundary integral equations are derived for a functionally graded material, which is used as the numerical tool in our analysis. The material property gradient is assumed to vary unsymmetrically in the direction normal to the fracture plane. A parametric study has been performed by systematically varying the inhomogeneity length scale. In comparison with a homogeneous material, an unsymmetric functionally graded material offers lesser fracture resistance. The fracture resistance progressively decreases with increase in inhomogeneity, quantified by the increase in crack sliding displacement jumps, crack tip velocities and accelerations. The material property inhomogeneity affects only the transient crack propagation velocities, while the quasi-steady-state velocity remains unaltered.     

功能梯度材料裂纹能量释放率

 提出一个评价功能梯度材料裂纹扩展能量释放率的简单方法， 用双悬臂梁柔度法推导了典型梯度材料内裂纹沿梯度方向扩展的能量 释放率解析表达式，用算例分析了裂纹扩展方向不同能量释放率的差 别.     

Thermal loading in multi-layered and/or functionally graded materials: Residual stress field,delamination, fatigue and related size effects

The stress field and fracture propagation due to thermal loading in multi-layered and/or functionally graded composite materials are extensively analysed. Regarding fracture, we have focused the attention on delamination between the layers due to brittle or fatigue thermally induced crack propagations. The statically indeterminate stress analysis is solved coupling equilibrium, compatibility and constitutive equations. Fracture analysis is based on the classical Griffith’s criterion rewritten for composite structures under thermal loading. As an example, a two-layer prismatic structure is considered, each layer being composed by a different functionally graded material. The solution is particularized for the case of a linear grading. The size and shape effects are discussed and an optimization procedure is proposed. A numerical application of the findings to hard metal and diamond based cutters concludes the paper.     

功能梯度板条断裂分析

现存文献关于功能梯度材料断裂问题的研究大都假设材料性质为坐标的指数函数或幂函数,而对其它函数形式较少采用。本文假设功能梯度材料剪切模量为坐标的双曲函数,而泊松比为常量,研究功能梯度板条的混合型裂纹问题。利用Fourier积分变换技术将混合边值问题转化为一对奇异积分方程,通过数值求解奇异积分方程获得含裂纹功能梯度板条分别在剪切和法向载荷作用下的I型和Ⅱ型应力强度因子,并讨论了材料的非均匀性和裂纹相对尺寸对裂纹尖端应力强度因子的影响。     

Elliptical crack normal to functionally graded interface of bonded solids

This paper presents an analysis of an elliptical crack that is perpendicular to a functionally graded interfacial zone between two fully bonded solids. The functionally graded interfacial zone is treated as a non-homogeneous solid layer with its elastic modulus varying in the thickness direction. A generalized Kelvin solution based boundary element method is employed for the calculation of the stress intensity factors associated with the three-dimensional crack problem. The elliptical crack surface is subject to either uniform normal traction or uniform shear traction. The stress intensity factors are examined by taking into account the effects of the non-homogeneity parameter and thickness of the functionally graded interfacial zone, as well as the crack distance to the zone. The SIF values are further incorporated into the S-criterion for prediction of crack growth. The paper presents the most possible direction and location of the elliptical crack growth under an inclined tensile (or compressive) load. The paper further presents results of the critical external loads that would cause the elliptical crack to grow at the most possible location and along the most possible direction. The paper also examines the effects of external load direction and material and geometrical parameters on the critical loads.     

Transient thermal stress analysis of a functionally graded magneto-electro-thermoelastic strip due to nonuniform surface heating

This article is concerned with the theoretical analysis of the functionally graded magneto-electro-thermoelastic strip due to unsteady and nonuniform surface heating in the width direction. We analyze the transient thermal stress problem for a functionally graded strip constructed of the anisotropic and linear magneto-electro-thermoelastic materials using a laminated composite mode as one of theoretical approximation. The transient two-dimensional temperature is analyzed by the methods of Laplace and finite sine transformations. We obtain the solution for the simply supported and functionally graded magneto-electro-thermoelastic strip under a plane strain state. As an illustration, we carried out numerical calculations for a functionally graded strip composed of piezoelectric BaTiO₃ and magnetostrictive CoFe₂O₄, and examined the behaviors in the transient state for temperature change, stress, electric potential and magnetic potential distributions. Furthermore, the effects of the nonhomogeneity of material on the stresses, electric potential, and magnetic potential are investigated.     

功能梯度材料的平面断裂力学分析

针对材料参数在厚度方向可能按任意连续变化的梯度材料，给出了一个新的分层模型，利用该模型求解了面内加载下梯度界面层和涂层中的界面裂纹问题，借助Fburier积分技术和传递矩阵方法，将该问题化为一个Cauchy型奇异积分方程，通过数值求解，得到感兴趣的应力强度因子，对不同形式的杨氏模量和泊松比，计算了界面裂纹应力强度因子，结果表明泊松比的变化形式对应力强度因子影响不大，可当作常数处理，而杨氏模量的影响则很大。     

Basic solutions of multiple parallel symmetric mode-III cracks in functionally graded piezoelectric/piezomagnetic material plane

The Schmidt method is adopted to investigate the fracture problem of multiple parallel symmetric and permeable finite length mode-III cracks in a functionally graded piezoelectric/piezomagnetic material plane. This problem is formulated into dual integral equations, in which the unknown variables are the displacement jumps across the crack surfaces. In order to obtain the dual integral equations, the displacement jumps across the crack surfaces are directly expanded as a series of Jacobi polynomials. The results show that the stress, the electric displacement, and the magnetic flux intensity factors of cracks depend on the crack length, the functionally graded parameter, and the distance among the multiple parallel cracks. The crack shielding effect is also obviously presented in a functionally graded piezoelectric/piezomagnetic material plane with mul- tiple parallel symmetric mode-III cracks.     

Dynamic fracture behaviors of cracks in a functionally graded magneto-electro-elastic plate

In this paper the dynamic anti-plane problem for a functionally graded magneto-electro-elastic plate containing an internal or an edge crack parallel to the graded direction is investigated. The crack is assumed to be magneto-electrically impermeable. Integral transforms and dislocation density functions are employed to reduce the problem to Cauchy singular integral equations. Field intensity factors and energy release rate are derived, analyzed and partially calculated numerically. The effects of material graded index, loading combination parameter (including size and direction) and geometry criterion of the plate on the dynamic energy release rate are shown graphically. Numerical results indicate that increasing the graded index can all retard the crack extension, and that both the applied magnetic field loadings and electric field loadings play a dominant role in the dynamic fracture behaviors of crack tips.     

Stress intensity factors for three-dimensional cracks in functionally graded materials using enriched finite elements

A three-dimensional enriched finite-element methodology is presented to compute stress intensity factors for three-dimensional cracks contained in functionally graded materials (FGMs). A general-purpose 3D finite-element based fracture analysis program, FRAC3D, is enhanced to include this capability. First, using available solutions from the literature, comparisons have been made in terms of stresses under different loading conditions, such as uniform tensile, bending and thermal loads. Mesh refinement studies are also performed. The fracture solutions are obtained for edge cracks in an FGM strip and surface cracks in a finite-thickness FGM plate and compared with existing solutions in the literature. Further analyses are performed to study the behavior of stress intensity factor near the free surface where crack front terminates. It is shown that three-dimensional enriched finite elements provide accurate and efficient fracture solutions for three-dimensional cracks contained in functionally graded materials.     

Cylindrical interface crack in a hollow layered functionally graded cylinder under static torsion

A hollow functionally graded composite cylinder under static torsion, which consists of an inner and outer elastic circular tube with a cylindrical interface crack, is studied in this work. By utilizing Fourier integral transform method, the mixed boundary value problem is reduced to a Cauchy singular integral equation, from which the numerical results of the stress intensity factor are obtained by the Lobatto–Chebyshev quadrature technique. Numerical results demonstrate the coupled effects of geometrical, physical, and functionally graded parameters on the interfacial fracture behavior.