梯度材料中矩形裂纹的对偶边界元方法分析

采用对偶边界元方法分析了梯度材料中的矩形裂纹. 该方法基于层状材料基本解,以非裂纹边界的位移和面力以及裂纹面的间断位移作为未知量. 位移边界积分方程的源点配置在非裂纹边界上,面力边界积分方程的源点配置在裂纹面上. 发展了边界积分方程中不同类型奇异积分的数值方法. 借助层状材料基本解,采用分层方法逼近梯度材料夹层沿厚度方向力学参数的变化. 与均匀介质中矩形裂纹的数值解对比,建议方法可以获得高精度的计算结果. 最后,分析了梯度材料中均匀张应力作用下矩形裂纹的应力强度因子,讨论了梯度材料非均匀参数、夹层厚度和裂纹与夹层之间相对位置对应力强度因子的影响.      

动态载荷下功能梯度复合材料的圆币形裂纹问题

研究了动态载荷下功能梯度材料中的圆币形裂纹问题．假设材料为横观各向同性,并且含有多个垂直于厚度方向的裂纹,材料参数沿轴向（与裂纹面垂直的方向）为变化的,沿该方向将材料划分为许多单层,各单层材料参数为常数,利用Hankel变换祛,在Laplace域内推导出了控制问题的对偶积分方程组．利用Laplace数值反演,得出了裂纹尖端的动态应力强度因子和能量释放率．研究了含两个裂纹的功能梯度接头结构,分析了材料非均匀性参数对应力强度因子和能量释放率的影响．     

功能梯度材料涂层平面裂纹分析

研究粘接于均质基底材料上功能梯度涂层平面裂纹问题. 假设功能梯度材料剪切模量的倒数为坐标的线性函数，而泊松比为常数. 采用Fourier变换和传递矩阵法将该混合边值问题化为奇异积分方程组，通过数值求解获得 应力强度因子. 考察了材料梯度变化形式、结构几何尺寸和材料梯度参数对裂纹应力强度因子的影响，发现 功能梯度材料涂层尺寸、裂纹长度以及材料梯度参数均对应力强度因子有显著影响.     

功能梯度材料裂纹尖端动态应力场

研究受反平面剪切作用的功能梯度材料动态裂纹问题,通过积分变换-对偶积分方程方法推出了裂纹尖端动态应力场,时间域内的动态应力强度因子由Laplace数值反演获得,研究结果表明功能梯度材料的梯度越大,相应的裂纹问题的动态应力强度因子值越低。     

功能梯度压电/压磁双材料的周期界面裂纹问题

研究反平面载荷作用下压电/压磁双材料的周期界面裂纹问题,压电/压磁双材料由有限厚度的功能梯度压电层和功能梯度压磁层粘结而成.为便于分析,假设压电层和压磁层的材料性质沿着裂纹的法线方向呈指数变化,基于分离变量和Hilbert核奇异积分方程方法,获得应力强度因子的数值解.数值算例讨论层厚、周期带长度、梯度参数以及材料参数变动等对应力强度因子的影响.结果发现层厚以及裂纹间距的增大会降低裂纹尖端应力强度因子,梯度参数的改变对应力强度因子也有显著的影响.材料参数变动的讨论发现弹性参数的变动对应力强度因子影响最大,其次为电参数,磁参数的变动对应力强度因子影响最小.     

功能梯度材料涂层平面运动裂纹分析

研究粘结于均匀材料基底上功能梯度材料涂层平面运动裂纹问题,假设功能梯度材料剪切模量和密度为坐标的指数函数,而泊松比为常数.采用Fourier变换和传递矩阵法将该混合边值问题转化为一对奇异积分方程,通过数值求解奇异积分方程组获得功能梯度材料涂层平面运动裂纹的应力强度因子.考察了结构几何尺寸、裂纹运动速度以及材料梯度参数对运动裂纹的应力强度因子的影响,发现材料梯度参数、结构几何尺寸、裂纹长度以及运动速度均对功能梯度材料动态断裂行为有显著影响.     

功能梯度板条断裂分析

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

功能梯度压电压磁材料粘结的Ⅲ型裂纹问题

利用奇异积分方程方法研究两个半无限大的功能梯度压电压磁材料粘结,在渗透和非 渗透边界条件下的III型裂纹问题. 首先通过积分变换构造出原问题的形式解，然 后利用边界条件通过积分变换与留数定理得到一组奇异积分方程， 最后利用Gauss-Chebyshev方法进行数值 求解，讨论材料参数、材料非均匀参数以及裂纹几何形状等对裂纹尖端应力 强度因子的影响. 从结果中可以看出，压电压磁复合材料中反平面问题的应力奇异性 形式与一般弹性材料中的反平面问题应力奇异形式相同，但材料梯度参数对功能梯度压电压 磁复合材料中的应力强度因子和电位移强度因子有很大的影响.     

任意梯度分布功能梯度涂层平面裂纹分析

提出可以分析任意梯度功能梯度材料的分层模型,并采用该模型研究功能梯度涂层平面裂纹问题.采用Fourier变换和传递矩阵法将该混合边值问题化为奇异积分方程组,通过数值求解获得应力强度因子.考察了分层模型的有效性,以及材料梯度变化形式、结构几何尺寸和材料梯度参数对裂纹应力强度因子的影响,发现结构几何尺寸、材料梯度变化形式、...     

功能梯度板条Ⅲ型运动裂纹问题研究

现存文献关于梯度材料断裂问题的研究大都是假设材料参数为坐标的指数函数或幂函数,而其它函数形式较少采用.本文假设功能梯度材料剪切模量和密度的倒数均为坐标的线性函数,而泊松比为常量,研究功能梯度板条的反平面运动裂纹问题.利用Fourier积分变换技术和传递矩阵法将混合边值问题化为一对奇异积分方程,通过数值求解奇异积分方程获得板条运动裂纹在反平面载荷作用下的动态应力强度因子,并讨论了裂纹运动速度、裂纹相对尺寸、以及材料非均匀性对动态应力强度因子的影响,结果证明梯度参数、裂纹速度和几何尺寸对材料动态断裂行为有显著影响.     

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.     

Stress intensity factors for an oblique edge crack in a coating/substrate system with a graded interfacial zone under antiplane shear

This paper investigates the edge crack problem for a coating/substrate system with a functionally graded interfacial zone under the condition of antiplane deformation. With the interfacial zone being modeled by a nonhomogeneous interlayer having the continuously varying shear modulus between the dissimilar, homogeneous phases of the coated medium, the coating is assumed to contain an edge crack at an arbitrary angle to the interfacial zone. The Fourier integral transform method is used in conjunction with the coordinate transformations of basic field variables. Formulation of the proposed crack problem is then reduced to solving a singular integral equation with a generalized Cauchy kernel. The mode III stress intensity factors are defined and evaluated in terms of the solution to the integral equation. In the numerical results, the values of the stress intensity factors are plotted, illustrating the effects of the crack orientation angle for various material and geometric combinations of the coating/substrate system with the graded interfacial zone.     

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.     

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

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

Effects of material gradients on transient dynamic mode-III stress intensity factors in a FGM

This paper presents a transient dynamic crack analysis for a functionally graded material (FGM) by using a hypersingular time-domain boundary integral equation method. The spatial variations of the material parameters of the FGM are described by an exponential law. A numerical solution procedure is developed for solving the hypersingular time-domain traction BIE. To avoid the use of time-dependent Green’s functions which are not available for general FGM, a convolution quadrature formula is adopted for approximating the temporal convolution, while a Galerkin method is applied for the spatial discretization of the hypersingular time-domain traction BIE. Numerical results for the transient dynamic stress intensity factors for a finite crack in an infinite and linear elastic FGM subjected to an impact anti-plane crack-face loading are presented and discussed. The effects of the material gradients of the FGM on the transient dynamic stress intensity factors and their dynamic overshoot over the corresponding static stress intensity factors are analyzed.     

Interfacial fracture analysis of bonded dissimilar strips with a functionally graded interlayer under antiplane deformation

This paper provides the solution to the problem of dissimilar, homogeneous semi-infinite strips bonded through a functionally graded interlayer and weakened by an embedded or edge interfacial crack. The bonded system is assumed to be under antiplane deformation, subjected to either traction-free or clamped boundary conditions along its bounding planes. Based on the Fourier integral transform, the problem is formulated in terms of a singular integral equation which has a simple Cauchy kernel for the embedded crack and a generalized Cauchy kernel for the edge crack. In the numerical results, the effects of geometric and material parameters of the bonded system on the crack-tip stress intensity factors are presented in order to quantify the interfacial fracture behavior in the presence of the graded interlayer.     

DYNAMIC BEHAVIOR OF SEVERAL CRACKS IN FUNCTIONALLY GRADED STRIP SUBJECTED TO ANTI-PLANE TIME-HARMONIC CONCENTRATED LOADS

This paper contains a theoretical formulations and solutions of multiple cracks sub- jected to an anti-plane time-harmonic point load in a functionally graded strip. The distributed dislocation technique is used to construct integral equations for a functionally graded material strip weakened by several cracks under anti-plane time-harmonic load. These 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 evaluate the stress intensity factor and strain energy density factors （SEDFs） for multiple cracks with differ- ent configurations. Numerical calculations are presented to show the effects of material properties and the crack configuration on the dynamic stress intensity factors and SEDFs of the functionally graded strip with multiple curved cracks.     

FRACTURE ANALYSIS OF A FUNCTIONALLY GRADED STRIP UNDER PLANE DEFORMATION

In this paper the plane elasticity problem for a functionally graded strip containing a crack is considered. It is assumed that the reciprocal of the shear modulus is a linear function of the thickness-coordinate, while the Possion's ratio keeps constant. By utilizing the Fourier transformation technique and the transfer matrix method, the mixed boundary problem is reduced to a system of singular integral equations that are solved numerically. The influences of the geometric parameters and the graded parameter on the stress intensity factors and the strain energy release rate are investigated. The numerical results show that the graded parameters, the thickness of the strip and the crack size have significant effects on the stress intensity factors and the strain energy release rate.     

Three-dimensional mixed-mode stress intensity factors for cracks in functionally graded materials using enriched finite elements

Three-dimensional enriched finite elements are used to compute mixed-mode stress intensity factors (SIFs) for three-dimensional cracks in elastic functionally graded materials (FGMs) that are subject to general mixed-mode loading and constraint conditions. The method, which advantageously does not require special mesh configuration/modifications and post-processing of finite element results, is an enhancement of previous developments applied so far on isotropic homogeneous and isotropic interface cracks. The spatial variation of FGM material properties is taken into account at the level of element integration points. To validate the developed method, two- and three-dimensional mixed-mode fracture problems are selected from the literature for comparison. Two-dimensional cases include: inclined central crack in a large FGM medium under uniform tensile strain and stress loadings, a slanted crack in a finite-size FGM plate under exponentially varying tensile stress loading and an edge crack in a finite-size plate under shear traction load. The three-dimensional example models a deflected surface crack in a finite-size FGM plate under uniform tensile stress loading. Comparisons between current results and those from analytical and other numerical methods yield good agreement. Thus, it is concluded that the developed three-dimensional enriched finite elements are capable of accurately computing mixed-mode fracture parameters for cracks in FGMs.