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

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

含裂纹功能梯度Euler-Bernoulli梁和Timoshenko梁的屈曲载荷计算与分析

含裂纹构件的屈曲载荷是结构是否安全的判定准则之一, 其计算与分析也是结构健康监测和安全评价中关注的重要问题。基于Euler-Bernoulli梁理论和Timoshenko梁理论, 建立了一种求解含裂纹功能梯度材料梁的屈曲载荷计算方法。首先裂纹导致的构件截面转角不连续性由转动弹簧模型进行模拟, 再根据功能梯度材料Euler-Bernoulli梁和Timoshenko梁的屈曲控制方程及其闭合解, 由传递矩阵法建立了求解含裂纹功能梯度材料梁在多种边界条件下屈曲载荷的循环递推公式和特征行列式, 使问题通过降阶的方法得到快速准确的解答。数值算例研究了剪切变形、 裂纹的不同数目及位置、 材料参数变化、 长细比和不同边界约束条件等对含裂纹功能梯度材料梁屈曲载荷的影响。结果表明该方法可以简单、 方便和准确地计算不同数目裂纹和任意边界条件下功能梯度材料梁的屈曲问题。      

含任意方向裂纹功能梯度材料的应力分析研究

功能梯度材料是在航空航天领域的需求背景下发展起来的,但由于生产技术及工作环境等方面的原因,功能梯度材料内部常常产生各种形式的裂纹并最终导致材料破坏,因此研究含任意方向裂纹功能梯度材料的断裂问题具有重要意义。以含有任意方向裂纹的功能梯度材料为对象,运用积分变换方法,给出了相应材料平面问题的位移场的形式解。通过引入辅助函数并利用相关条件,可将问题转化为求解一组带有Cauchy核的奇异积分方程,继而采用Lobatto-Chebyshev方法对奇异积分方程进行数值求解。最后分析了裂纹方向、材料非均匀指数、载荷条件对混合型应力强度因子的影响。      

一阶剪切理论下功能梯度梁与均匀梁静态解之间的相似关系

基于一阶剪切理论,研究了功能梯度材料Timoshenko 梁的静态弯曲解与对应的均匀材料梁的解的线性转换关系。通过比较功能梯度材料梁和均匀材料梁的无量纲控制方程,发现了它们弯曲解的线性相关性。在给定材料弹性模量沿横向非均匀变化规律后,可将功能梯度材料Timoshenko 梁在静载荷作用下的弯曲变形解用相同尺寸、相同载荷以及相同边界条件下的均匀材料Timoshenko 梁的弯曲变形解线性表示。这样,可将非均匀Timoshenko 梁弯曲问题的求解转化为对应的均匀材料Timoshenko 梁弯曲问题的求解和转换系数的计算,从而使得求解过程得以简化。     

梯度复合材料热传导分析的梯度单元法

给出了一种适用于梯度复合材料热传导分析的梯度单元, 采用细观力学方法描述材料变化的热物理属性, 通过线性插值和高阶插值温度场分别给出了4节点和8节点梯度单元随空间位置变化的热传导刚度矩阵。推导了在温度梯度载荷和热流密度载荷作用下, 矩形梯度板的稳态温度场和热通量场精确解。基于该精确解对比了连续梯度模型和传统的离散梯度模型的热传导有限元计算结果, 验证了梯度单元的有效性, 并讨论了相关参数对梯度单元的影响。结果表明, 梯度单元和均匀单元得到的温度场基本一致; 当热载荷垂直于材料梯度方向时, 梯度单元能够给出更加精确的局部热通量场; 当热载荷平行于材料梯度方向时, 4节点梯度单元性能恶化, 8节点梯度单元和均匀单元的计算结果与精确解吻合很好。     

在反平面变形下功能梯度压电板条的应力强度因子分析

利用Fourier变换技术将混合边值问题化为对偶积分方程,求解对偶积分方程得到应力强度因子以及电位移强度因子的表达式。最后通过数值计算讨论了材料参数、载荷条件以及裂纹的几何参数等对功能梯度压电材料中裂纹尖端应力强度因子的影响。     

非均匀热载荷作用下功能梯度梁的非线性静态响应

由于功能梯度材料结构沿厚度方向的非均匀材料特性,使得夹紧和简支条件的功能梯度梁有着相当不同的行为特征。该文给出了热载荷作用下,功能梯度梁非线性静态响应的精确解。基于非线性经典梁理论和物理中面的概念导出了功能梯度梁的非线性控制方程。将两个方程化简为一个四阶积分-微分方程。对于两端夹紧的功能梯度梁,其方程和相应的边界条件构成微分特征值问题;但对于两端简支的功能梯度梁,由于非齐次边界条件,将不会得到一个特征值问题。导致了夹紧与简支的功能梯度梁有着完全不同的行为特征。直接求解该积分-微分方程,得到了梁过屈曲和弯曲变形的闭合形式解。利用这个解可以分析梁的屈曲、过屈曲和非线性弯曲等非线性变形现象。最后,利用数值结果研究了材料梯度性质和热载荷对功能梯度梁非线性静态响应的影响。     

多铁性非均匀空心层合柱圆弧界面的反平面断裂分析

讨论反平面载荷作用下多铁性非均匀空心层合柱的圆弧界面裂纹问题,层合柱由梯度铁电层和梯度铁磁层粘接而成,界面处存在圆弧型裂纹。采用分离变量和Cauchy核奇异积分方程方法求解该断裂问题。通过讨论断裂参数的数值解,分析了梯度非均匀参数、几何与材料参数变动等对应力强度因子的影响。     

考虑热对材料参数影响的FGM梁热后屈曲特性研究

功能梯度梁热后屈曲特性研究对于推进功能梯度材料在航天器热防护设计中的应用有着重要意义。基于经典梁的几何非线性理论和物理中面的概念,建立了热载荷作用下功能梯度梁的运动微分方程,通过化简得到一个仅关于挠度的四阶微分-积分方程,并与固支边界条件构成特征值问题,分析研究了功能梯度梁的热后屈曲及在此基础上的振动问题。首先证明了通过哈密顿原理推导的运动方程是轴线可伸长理论的近似形式。接着考虑热对材料物性参数的影响,并对梁的长细比、功能梯度指数和温度比作了详细分析,研究了这些参数对热后屈曲路径和后屈曲振动的影响规律。结果表明：只有在长细比较大时才可以不考虑温度对材料物性参数的影响,否则误差较大;长细比、功能梯度指数和温度比的增大会增大无量纲热屈曲载荷,同时使屈曲路径和频率-载荷曲线向热载荷增大的方向移动。     

倾斜应力作用下垂直于功能梯度材料夹层的币型裂纹扩展问题

分析了功能梯度材料中币型裂纹扩展问题。该裂纹体受有与裂纹面成任意角度的张应力或压应力,裂纹垂直于无限域中功能梯度材料夹层。假定非均匀介质的功能梯度材料夹层与两个半无限域完全结合,其弹性模量沿厚度方向变化。利用已发表的裂纹应力强度因子数据和线弹性断裂力学的叠加原理,将应力强度因子耦合于最小应变能密度因子断裂判据,讨论了裂纹扩展的临界荷载;并讨论了荷载方向和材料性质对临界荷载的影响。     

Fracture toughness of the + 45° / – 45° interface of a laminate composite

Experiments are carried out to determine the delamination toughness for a crack along the interface between two transversely isotropic materials. The material chosen for study consists of carbon fibers embedded within an epoxy matrix. A crack is introduced between two layers of this material, with fibers in the upper layer along the + 45°-direction and those in the lower layer along the − 45°-direction both with respect to the crack plane. The Brazilian disk specimen is employed in the testing. To calibrate the specimens, stress intensity factors are obtained which result from the applied load, as well as residual curing stresses. It may be noted that all three modes are coupled, leading to a three-dimensional problem. The finite element method and a mechanical M-integral are employed to determine the stress intensity factors arising from the applied load. For the residual stresses, a three-dimensional conservative thermal M-integral is presented for stress intensity factor determination. The stress intensity factors found for the applied load and residual stresses are superposed to obtain a local energy release rate, together with two phase angles. From the load at fracture, the critical interface energy release rate or interface toughness as a function of phase angles ψ and ϕ is determined. Results are compared to a fracture criterion.     

The residual stress intensity factors for cold‐worked cracked holes: a technical note

Cold‐working of riveted holes reduces the stress intensity factor associated with cracks that may develop at the hole boundary, by creating a compressive residual stress field. The residual stress field is determined using the finite‐element method and the reduction of the stress intensity factor for different values of the interference is evaluated with the weight function method, in the case of an infinite plate made from an elastic–perfectly plastic material, and having a hole with two symmetrical cracks. Once the weight function of the structure is known, further calculation of the stress intensity factors for different loadings such as a remote uniform stress, or a point load that simulates the action of the rivet can be performed without difficulty.     

The Interaction between a Piezoelectric Screw Dislocation Dipole and an Elliptic Blunt Crack in Elliptical Inhomogeneity

The interaction problem of a piezoelectric screw dislocation dipole with a confocal elliptic blunt crack in elliptical inhomogeneity subjected to remote anti-plane stress field and in-plane electric field is investigated by using the complex method of elasticity. The exact closed-form solutions of a series of quantities, such as singularity stress field, image force and image torque acting on the center of screw dislocation dipole, stress intensity factor and electric displacement intensity factor of crack tip, energy release rate, and generalized strain energy density are obtained. Then the influence laws of remote load, the dip angle of dislocation dipole, the size of blunt crack, and the material constants on the quantities are analyzed. The numerical results show that the image force, image torque, stress intensity factor, and electric displacement intensity factor make periodic variation as the dip angle of dislocation dipole; the energy release rate of crack tip is negative when subjected to pure electric field, however, it can be positive or negative when subjected to the combined action of mechanical field and electric field; the sharp crack is not easy to expand in some combined action of mechanical field and electric field.     

A wedge crack in an anisotropic material under antiplane shear

A crack emanating from the apex of an infinite wedge in an anisotropic material under antiplane shear is investigated. An isotropic wedge crack subjected to concentrated forces is first solved by using the conformal mapping technique. The solution of an anisotropic wedge crack is obtained from that of the transformed isotropic wedge crack based on a linear transformation method. Expressions for the stress intensity factor for the anisotropic wedge crack with both concentrated and distributed loads are derived. The stress intensity factors are numerically calculated for generally orthotropic wedge cracks with various crack and wedge angles as well as anisotropic parameters.     

基于应力强度因子评估WCP形状对WCP/Fe复合材料热疲劳裂纹扩展行为的影响

应力强度因子在断裂力学中广泛应用于预测由远程载荷或残余应力引起的裂缝尖端附近应力状态。本文基于平面应力条件下应力强度因子建立WC_P形状与其尖端应力之间的规律,利用有限元分析软件对含不同形状WC_P的WC_P/Fe复合材料的热应力进行模拟仿真,研究WC_P形状对WC_P/Fe复合材料热疲劳裂纹扩展行为的影响。研究结果表明,WC_P的形状显著影响应力强度因子,进一步影响WC_P/Fe复合材料的热疲劳裂纹扩展行为。含球状和不规则状WC_P的WC_P/Fe复合材料的极限抗压强度分别约为460 MPa和370 MPa。含不规则状WC_P的WC_P/Fe复合材料因应力集中而容易产生脆性开裂现象。通过热震实验进行验证,发现实验结果与模拟仿真结果相近,说明有限元法的准确性,同时为WC_P/Fe复合材料的热疲劳裂纹扩展行为研究提供科学依据和理论基础。      

Experimental Investigations on Deformation and Fracture Behavior of Glass Sphere Filled Epoxy Functionally Graded Materials

In this paper, deformation and fracture behavior of glass sphere filled epoxy functionally graded materials (FGM) are numerically evaluated and experimentally studied. The fabrication of the FGM is described in detail, and the spatial gradation of elastic modulus and the microscopic structure in FGM are measured and analyzed. The deformation and fracture characterization of the FGM specimen with a crack oriented along the direction of the elastic gradient under three point bend are studied by the experimental and the finite element method. The influences of crack location at both the stiff and the compliant sides of the FGM specimen on crack initiation, deformation field and stress intensity factor are analyzed. The results are: (a) The neutral-axis in the FGM specimen under three-point-bending will shift toward the stiffer side; (b) The initial fracture load increases with the increase of elastic modulus at the crack tip; (c) The elastic gradients shield a crack on the compliant side and lower the stress intensity factor when compared to the one with crack on the stiff side. These results will be useful for better design and reliable evaluation of FGM.     

The treatment of thermal and residual stresses in fracture assessments

The behaviour of a crack under combined mechanical and thermal/residual loads is described in terms of a reference stress or equivalent mechanical loading. The reference stress can be readily established from a knowledge of the stress intensity factor for the thermal/residual load, the magnitude of the mechanical load, and the material stress-strain curve. The result enables the J-integral to be evaluated and is also used to define a procedure for the inclusion of thermal/residual stresses in the CEGB defect assessment method (R6). In particular, a detailed procedure is given for treating thermal/residual stresses in recent revisions of R6 which cater for structures made of materials which strongly work-harden.     

压电材料中反平面渗透裂纹的解

用复变函数的保角映射法,采用可渗透边界条件,研究了含裂纹的无限大压电材料在平面内电场和反平面荷载作用下的耦合场,得到了精确的解和场强度因子以及能量释放率。结果表明,电场强度在裂尖没有奇异性,应变、应力、电位移具有1/2阶的奇异性,能量释放率总是正的。     

FATIGUE CRACK GROWTH UNDER MIXED-MODE I AND II LOADING

Abstract— Mixed-mode fatigue crack growth has been studied using four point bend specimens under asymmetric loads. A detailed finite element analysis provides the stress intensity factors for curved cracks under different mixed-mode load conditions. Both fatigue crack growth direction and crack growth rate are studied. The maximum tangential stress and the minimum strain energy density criteria were found to provide satisfactory predictions of the crack growth directions. An effective stress intensity factor was used to correlate the fatigue crack growth rates successfully. It is found that the use of mode I fatigue crack growth rate properties results in a conservative crack growth rate prediction for mixed-mode load conditions.     

Determination of approximate point load weight functions for embedded elliptical cracks

This paper presents the application of weight function method for the calculation of stress intensity factors in embedded elliptical cracks under complex two-dimensional loading conditions. A new general mathematical form of point load weight function is proposed based on the properties of weight functions and the available weight functions for two-dimensional cracks. The existence of this general weight function form has simplified the determination of point load weight functions significantly. For an embedded elliptical crack of any aspect ratio, the unknown parameters in the general form can be determined from one reference stress intensity factor solution. This method was used to derive the weight functions for embedded elliptical cracks in an infinite body and in a semi-infinite body. The derived weight functions are then validated against available stress intensity factor solutions for several linear and non-linear stress distributions. The derived weight functions are particularly useful for the fatigue crack growth analysis of planer embedded cracks subjected to fluctuating non-linear stress fields resulting from surface treatment (shot peening), stress concentration or welding (residual stress).