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

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

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

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

非均匀复合材料的动态热弹性断裂力学分析

对非均匀复合材料的动态热弹性断裂力学问题进行了研究,假设材料参数沿厚度方向为变化的,沿该方向将复合材料划分为许多单层,取每一单层材料参数为常数,应用Fourier变换法,在Laplace域内推导出了控制问题的奇异积分方程组,给出了热应力强度因子的表达式,然后利用Laplace数值反演,得出了裂纹尖端的动态应力强度因子．本文的方法具有以下特点：（1）多个垂直于厚度方向的裂纹,（2）材料可以为正交各向异性：（3）考虑了惯性效应．作为算例,研究了带有两个裂纹的功能梯度结构,分析了材料参数的变化对应力强度因子的影响．     

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

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

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

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

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

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

正交各向异性功能梯度材料反平面裂纹尖端应力场

采用积分变换-对偶积分方程方法,研究了正交各向异性功能梯度材料反平面裂纹问题,文中假定材料沿两个主轴方向的剪切模量成比例按双参数梯度模型变化,通过求解对偶积分程并考虑变形Bessel函数的渐特性,推导出了裂纹尖端应力场,最后考察了材料非均匀性及正交性对应力强度因子的影响。     

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

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

圆柱形界面相裂纹在扭转载荷作用下的动态断裂

本文研究了位于界面相中的圆柱形界面裂纹的扭转冲击问题.采用Laplace、Fourier变换和位错密度函数将混合边值问题转化为求解Cauchy核奇异积分方程,利用Laplace数值反演技术计算了动态应力强度因子.讨论了材料特性和结构的几何尺寸对动态应力强度因子的影响.结果表明,随着界面相厚度的增加,无量纲化的动态应力强度因子减小.当裂纹靠近剪切弹性模量大的材料时,无量纲化的动态应力强度因子增大,反之减小.界面相两侧不同的材料组合对裂尖动态应力强度因子的影响是随着剪切弹性模量和质量密度的比值的增加而减小.界面相中裂纹长度对裂尖动态应力强度因子的影响比其他因素的影响大.     

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

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

功能梯度条共线Griffith裂纹反平面剪切冲击

研究正交各向异性功能梯度条中多个共线Griffith裂纹的反平面剪切冲击问题．材料两个方向的剪切模量假定按比例同时以特定的梯度变化．采用Laplace和Fourier变换及引进位错密度函数将问题化为求解Cauchy奇异积方程，进而化为代数方程数值求解．考查材料非均匀性、正交性和功能梯度条高度对裂尖动态断裂特性的影响．动应力强度因子的数值结果显示：增加剪切模量的梯度和(或)增加垂直于裂纹面方向的剪切模量，可以抑制动应力强度因子的幅度；若功能梯度条较薄，增大条形域的高度也可抑制裂纹扩展．     

Torsional impact of transversely isotropic solid with functionally graded shear moduli and a penny-shaped crack

The torsional impact response of a penny-shaped crack in an unbounded transversely isotropic solid is considered. The shear moduli are assumed to be functionally graded such that the mathematics is tractable. Laplace transform and Hankel transform are used to reduce the problem to solving a Fredholm integral equation. The crack tip stress fields are obtained. Investigated are the influence of material nonhomogeneity and orthotropy on the dynamic stress intensity factor. The peak value of the dynamic stress intensity factor can be suppressed by increasing the shear moduli's gradient and/or increasing the shear modulus in a direction perpendicular to the crack surface.     

Dynamic response for functionally graded materials with penny-shaped cracks

This paper provides a method for studying the penny-shaped cracks configuration in functionally graded material(FGM) structures subjected to dynamic or steady loading. It is assumed that the FGMs are transversely isotropic and all the material properties only depend on the axial coordinatez. In the analysis, the elastic region is treated as a number of layers. The material properties are taken to be constants for each layer. By utilizing the Laplace transform and Hankel transform technique, the general solutions for the layers are derived. The dual integral equations are then obtained by introducing the mechanical boundary and layer interface conditions via the flexibility/stiffness matrix approach. The stress intensity factors are computed by solving dual integral equations numerically in Laplace transform domain. The solution in time domain is obtained by utilizing numerical Laplace inverse. The main advantage of the present model is its ability for treating multiple crack configurations in FGMs with arbitrarily distributed and continuously varied material properties by dividing the FGMs into a number of layers with the properties of each layer slightly different from one another. This work was supported by Failure Mechanics Laboratory of State Education Commission and the Post-doctor Research Fund of China.     

Pre-kinking of a moving crack in a magnetoelectroelastic material under in-plane loading

A constant moving crack in a magnetoelectroelastic material under in-plane mechanical, electric and magnetic loading is studied for impermeable crack surface boundary conditions. Fourier transform is employed to reduce the mixed boundary value problem of the crack to dual integral equations, which are solved exactly. Steady-state asymptotic fields near the crack tip are obtained in closed form and the corresponding field intensity factors are expressed explicitly. The crack speed influences the singular field distribution around the crack tip and the effects of electric and magnetic loading on the crack tip fields are discussed. The crack kinking phenomena is investigated using the maximum hoop stress intensity factor criterion. The magnitude of the maximum hoop stress intensity factor tends to increase as the crack speed increases.     

TORSIONAL IMPACT RESPONSE OF A PENNY-SHAPED CRACK IN A FUNCTIONAL GRADED STRIP

The torsional impact response of a penny-shaped crack in a nonhomogeneous strip is considered. The shear modulus is assumed to be functionally graded such that the mathematics is tractable. Laplace and Hankel transforms were used to reduce the problem to solving a Fredholm integral equation. The crack tip stress field is obtained by considering the asymptotic behavior of Bessel function. Explicit expressions of both the dynamic stress intensity factor and the energy density factor were derived. And it is shown that, as crack driving force, they are equivalent for the present crack problem. Investigated are the effects of material nonhomogeneity and strip‘s highness on the dynamic fracture behavior.Numerical results reveal that the peak of the dynamic stress intensity factor can be suppressed by increasing the nonhomogeneity parameter of the shear modulus, and that the dynamic behavior varies little with the adjusting of the strip‘ s highness.     

Mixed-mode fracture analysis of orthotropic functionally graded materials under mechanical and thermal loads

Mixed-mode fracture problems of orthotropic functionally graded materials (FGMs) are examined under mechanical and thermal loading conditions. In the case of mechanical loading, an embedded crack in an orthotropic FGM layer is considered. The crack is assumed to be loaded by arbitrary normal and shear tractions that are applied to its surfaces. An analytical solution based on the singular integral equations and a numerical approach based on the enriched finite elements are developed to evaluate the mixed-mode stress intensity factors and the energy release rate under the given mechanical loading conditions. The use of this dual approach methodology allowed the verifications of both methods leading to a highly accurate numerical predictive capability to assess the effects of material orthotropy and nonhomogeneity constants on the crack tip parameters. In the case of thermal loading, the response of periodic cracks in an orthotropic FGM layer subjected to transient thermal stresses is examined by means of the developed enriched finite element method. The results presented for the thermally loaded layer illustrate the influences of the material property gradation profiles and crack periodicity on the transient fracture mechanics parameters.     

Transient dynamic response of a coated piezoelectric strip with a vertical crack

In this study, the dynamic response of a coated piezoelectric strip containing a crack vertical to the interfaces under normal impact load is considered. Based on the superposition principle and the integral transform techniques, the solution in the Laplace transformed plane is obtained in terms of a singular integral equation. The order of stress singularity around the tip of the terminated crack is also obtained. The singular integral equation is solved by using the Gauss–Jacobi integration formula, and the numerical Laplace inversion is then carried out to obtain the resulting dynamic stress and electric displacement intensities. The effects of the material properties and the geometric parameters on the dynamic stress intensity factor and the dynamic energy density factors are shown graphically.     

Impact failure prediction of Mode III crack in orthotropic functionally graded strip

Mode III impact of a crack in an orthotropic functionally graded strip is investigated. The shear moduli in two directions of the material are assumed to vary proportionately with gradient. Laplace transform and Fourier cosine transform are used to reduce the problem to solving a Fredholm integral equation. The crack tip stress field is obtained by considering the asymptotic behavior of Bessel function. Energy density factor criterion is applied to obtain the maximum of minimum energy density and direction of crack initiation. Numerical results are given graphically. The effects of orthotropy, nonhomogeneity and height of the strip on the energy density factor are discussed.