改进的虚裂纹闭合技术及其在复合材料脱层分析中的应用

提出一种改进的虚裂纹闭合技术（ＶＣＣＴ）对轮胎带束（一咱复合材料结构）不同脱层长度下的应变能释放率进行了计算。将强度裂纹闭合的位移作了一种可变约束施加在变形后的裂纹尖端,约束反力作了裂纹尖端力计算,应变能释放率由裂纹尖端力和约束位移的乘积得到。该方法的优点是简便易行,在物理理解上直观,实施起来方便。本文就用该技术对带束脱层进行了断裂力学分析,给出数值结果,并对带束的脱层寿命进行了评价。     

扭力轴三维裂纹扩展寿命仿真研究

对疲劳载荷谱作用下三维表面裂纹,采用双重边界元理论求解裂纹前沿的应力应变场,运用Forman理论、最小应变能密度法和Elber模型,计算裂纹前沿各点的扩展长度、扩展方向和应力强度因子等特征量.根据增量步下裂纹几何形状的改变,对裂纹面进行网格重划分和迭代计算,模拟三维裂纹的扩展,预测裂纹扩展寿命.扭力轴表面裂纹扩展的仿真结果表明该方法合理可行.     

压力容器高应变区应变疲劳裂纹扩展试验研究

针对压力容器几何结构不连续高应变区断裂和疲劳破坏特点,本文通过对应变循环条件下疲劳裂纹扩展影响因素和控制参量的理论分析,提出将△J作为应变循环条件下疲劳裂纹扩展速率的控制参量,并在EPRI弹塑性断裂分析工程方法基础上,提出了△J的工程计算方法,并通过模拟应变疲劳裂纹扩展的SCT试样试验,验证了△J工程计算方法的可靠性,根据SCT试样应变疲劳裂纹扩展试验研究结果,在考虑了压力容器高应变区几何形状和材质变化对应变疲劳裂纹扩展参量影响后,提出了一种压力容器高应变区应变疲劳裂纹扩展速率工程估算方法。     

裂纹转子分叉、混沌行为研究中的映射延拓综合法

介绍一种能全面计算周期激励下非线性系统周期响应，拟周期响应和混沌响应的新算法－映射延拓综合法，它可方便地确定拟周期响应和混沌响应对应的系统参数区间。应用此算法对具有非线性刚度的裂纹转子系统裂纹扩展故障特征问题进行了研究，得到了以裂纹深度为分叉参数的系统稳态响应分叉解图。     

三点弯曲试样动态应力强度因子计算研究

利用Hopkinson压杆对三点弯曲试样进行冲击加载,采集了垂直裂纹面距裂尖2mm和与裂纹面成60°距裂尖5mm处的应变信号。根据裂尖附近测试的应变信号计算试样的动态应力强度因子,并与有限元计算结果进行比较,结果表明由于裂尖有一段疲劳裂纹区,通过裂尖附近应变信号来计算动态应力强度因子时,如果裂尖位置确定不准及粘贴应变片位置不够准确对计算结果将带来很大影响。因此利用应变片法计算动态应力强度因子时,为了获得更准确的计算结果,在实验后应对试件裂纹面进行分析测量,重新确定裂尖位置,必要时需对应变片至裂尖距离进行修正后再计算动态应力强度因子值。     

复合材料刚度性能表征的协同损伤力学模型

针对复合材料层合板的弥散型损伤,提出一个刚度性能表征的协同损伤力学模型. 该模型兼顾了微观物理损伤响应和宏观材料刚度性能表征. 从微观角度,建立细观RVE 模型求解裂纹表面张开位移和滑开位移,以此定义损伤张量,并在宏观上通过对材料应变和损伤表面位移进行均匀化处理,建立单向板或层合板的损伤刚度矩阵和损伤张量之间的联系. 以基体裂纹为例,详细分析并建立了横向裂纹和纵向裂纹的损伤本构. 计算了[±θ/90₄]_S 铺层层合板中基体横向裂纹对刚度性能的影响,结果表明该方法能够准确地预测复合材料层合板由损伤导致的刚度性能衰减.      

基于Erdogan基本解边界元法计算应力强度因子

引入含裂纹问题基本解(Erdogan基本解)，提出了基于Erdogan基本解的样条虚边界 元法，并阐述了该法在实施过程中的特点与具体做法. 采用该方法详细分析了若干 典型裂纹问题，全面考察了方法的计算精度和收敛情况，以及在求解复杂裂纹问题方面 的能力. 结果显示，该方法具有精度高、收敛快、计算能力强等优点，是裂纹问题分析中 一种具有竞争力的通用计算方法.     

数字图像相关法测量局域变形场中形函数和模板尺寸的影响

数字图像相关法测量均匀变形场已被普遍接受, 其测量结果可与应变片测量结果比较. 然而, 在工程测量中, 针对局域变形场(应变高度集中, 如波特文-勒夏特利埃带、试件缺口附近和裂纹尖端等), 应变片受限于其尺寸, 其测量结果是接触面内的平均应变值. 此时, 采用数字图像相关法能够测量这些局域变形场. 但形函数和模板尺寸等计算参数对计算结果影响很大, 这也导致使用者很难判断计算结果的可靠性. 论文通过对合金拉伸实验获得的不同应变梯度的波特文-勒夏特利埃带和模拟生成的带的计算分析, 发掘了形函数和模板尺寸作用于计算结果的深层机制, 证明了二阶形函数比一阶形函数更适用于高度非均匀的局域变形场. 提出了在局域应变场测量中, 当一阶和二阶形函数计算结果的相对误差小于10% 时, 二阶形函数的结果是可靠的判据.      

一种基于应变梯度监测Ⅰ型裂纹的方法

基于线弹性断裂力学中I型裂纹的欧文解答,解析推导了在单向拉伸作用下无限大平板中I型裂纹尖端应变梯度场,建立了应变梯度与裂纹扩展之间的关联;基于挠曲电效应建立了电极化强度与应变梯度之间的力电耦合关系,提出了一种利用应变梯度传感器监测I型裂纹的方法,获知裂纹尖端坐标和裂纹扩展长度.本研究拟为应用应变梯度传感器对工程结构中裂纹扩展的实时监测提供初步的理论依据及方法.挠曲电感应技术在结构健康监测领域前景广阔.     

结构可靠性分析的支持向量机响应面法

针对隐式极限状态可靠性分析问题,提出了一种支持向量机响应面法,该方法采用了与经典响应面法类似的迭代思想,由支持向量回归机替代经典响应面法中的固定多项式函数来构建响应面,并结合一次二阶矩法形成迭代过程.在此基础上,还对训练样本提出了一种改进的选取方法,从而进一步提高方法的效率.文中将所提方法与多种经典可靠性分析方法的计算结果对比分析,改进的支持向量机响应面法精度较高,调用结构分析程序的次数最少.     

A new approach to on-line vibration monitoring and fault diagnosis

This paper presents a finite element method by deriving the stiffness matrix of cracked element from the stress intensity factor. The method can be used to evaluate the strain response of the cantilever beam with a single edge crack. By using the method of identification in time domain, the modal parameters are identified and at the same time the effect of crack closure is taken into consideration. Analysis has also been made of the sensitivity of the strain to the faults. Test results show that the strain response is more sensitive to faults than the displacement response. Then, a distinguishing function is introduced, which is more effective than any other existing ones in failure detection. Finally a new approach to on-line vibration monitoring and fault diagnosis is given.     

A boundary element modeling of fatigue crack growth in a plane elastic plate

This paper presents an extension of a boundary element method to fatigue growth analysis of mixed-mode cracked plane elastic bodies. The method consists of the non-singular displacement discontinuity element presented by Crouch and Starfield and the crack-tip displacement discontinuity element due to the author. In the boundary element implementation the left or the right crack-tip element is placed locally at the corresponding left or right crack tip on top of non-singular displacement discontinuity elements that cover the entire crack surface and the other boundaries. Crack growth is simulated with an incremental crack extension analysis based on the modified maximum strain energy density criterion. In numerical simulation, for each increment of crack extension, remeshing of existing boundaries is not required because of an intrinsic feature of the boundary element method. Crack growth is simulated by adding new boundary elements on the incremental crack extension to the previous crack boundaries. At the same time, the element characters of some related elements are adjusted according to the manner in which the boundary element method is implemented. Some numerical results of fatigue growth in a plane elastic plate with a center-inclined crack under uniaxial cyclic loading are given.     

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

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

Fatigue crack propagation life analysis of solder joints under thermal cyclic loading

Fatigue crack propagation life analysis of solder joints under thermal cyclic loadings was investigated by the strain energy release rate method using finite element analysis. A relationship between the crack-growth rate and the strain energy release rate was derived. Finite element simulations were carried out to investigate the effect of crack growth along the interface of solder and lead in a solder joint assembly. The crack propagation life of the solder joint with an interface crack was predicted from the derived relationship between crack growth rate and values of the strain energy release rate. It was found that crack propagation life is much higher than the crack initiation life.     

A new finite element method for strain gradient theories and applications to fracture analyses

A new compatible finite element method for strain gradient theories is presented. In the new finite element method, pure displacement derivatives are taken as the fundamental variables. The new numerical method is successfully used to analyze the simple strain gradient problems – the fundamental fracture problems. Through comparing the numerical solutions with the existed exact solutions, the effectiveness of the new finite element method is tested and confirmed. Additionally, an application of the Zienkiewicz–Taylor C1 finite element method to the strain gradient problem is discussed. By using the new finite element method, plane-strain mode I and mode II crack tip fields are calculated based on a constitutive law which is a simple generalization of the conventional J₂ deformation plasticity theory to include strain gradient effects. Three new constitutive parameters enter to characterize the scale over which strain gradient effects become important. During the analysis the general compressible version of Fleck–Hutchinson strain gradient plasticity is adopted. Crack tip solutions, the traction distributions along the plane ahead of the crack tip are calculated. The solutions display the considerable elevation of traction within the zone near the crack tip.     

A viscoplastic study of crack-tip deformation and crack growth in a nickel-based superalloy at elevated temperature

Viscoplastic crack-tip deformation behaviour in a nickel-based superalloy at elevated temperature has been studied for both stationary and growing cracks in a compact tension (CT) specimen using the finite element method. The material behaviour was described by a unified viscoplastic constitutive model with non-linear kinematic and isotropic hardening rules, and implemented in the finite element software ABAQUS via a user-defined material subroutine (UMAT). Finite element analyses for stationary cracks showed distinctive strain ratchetting behaviour near the crack tip at selected load ratios, leading to progressive accumulation of tensile strain normal to the crack-growth plane. Results also showed that low frequencies and superimposed hold periods at peak loads significantly enhanced strain accumulation at crack tip. Finite element simulation of crack growth was carried out under a constant ΔK-controlled loading condition, again ratchetting was observed ahead of the crack tip, similar to that for stationary cracks.A crack-growth criterion based on strain accumulation is proposed where a crack is assumed to grow when the accumulated strain ahead of the crack tip reaches a critical value over a characteristic distance. The criterion has been utilized in the prediction of crack-growth rates in a CT specimen at selected loading ranges, frequencies and dwell periods, and the predictions were compared with the experimental results.     

准脆性材料中强不连续问题的数值方法若干应用及其研究进展

综述了模拟准脆性材料开裂过程的数值计算方法的研究进展和工程应用,比较了表征强不连续问题的显式非连续模型和隐式非连续模型的优缺点.结合混凝土粘结裂纹, 重点讨论了嵌入非连续模型,扩展有限元方法和富集有限元技术等非连续方法的构造特征和本质区别.从各种富集方法的理论完备性考察,以假定发展应变为基础的嵌入非连续方法虽然可以解决混凝土开裂过程中的应力锁死,满足内部边界的静力平衡条件以及反映开裂后的位移不连续问题,但嵌入非连续所采用的富集函数在开裂单元中并不能满足协调条件,使非连续两侧的应变不独立. 其局限性是由于富集自由度在单元的水平上引入,而以单位分解为基础的扩展有限元和富集有限元的富集函数以节点自由度的方式引入,除具有嵌入非连续的优点, 还可以有效消除嵌入非连续引起裂纹两侧应变的相互影响.文中同时指出了网格重构技术,弥散裂纹模型的局限性以及扩展有限元和富集有限元技术在构造方式上的细微差别.对于节点自由度方式引入的富集函数, 其操作困难性在文中也作了说明.      

一种新的计算各向异性材料裂纹尖端应力强度因子杂交元法

利用有限元特征分析法研究了平面各向异性材料裂纹端部的奇性应力指数以及应力场和位移场的角分布函数,以此构造了一个新的裂纹尖端单元。文中利用该单元建立了研究裂纹尖端奇性场的杂交应力模型,并结合Hellinger-Reissner变分原理导出应力杂交元方程,建立了求解平面各向异性材料裂纹尖端问题的杂交元计算模型。与四节点单元相结合,由此提出了一种新的求解应力强度因子的杂交元法。最后给出了在平面应力和平面应变下求解裂纹尖端奇性场的算例。算例表明,本文所述方法不仅精度高,而且适应性强。     

黏弹性界面断裂分析的增量“加料”有限元法

提出了一种适用于黏弹性界面裂纹问题的增量“加料” 有限元方法. 利用弹性界面裂纹尖端位移场的解答,通过对应原理和拉普拉斯逆变换近似方法,得到了黏弹性界面裂纹的尖端位移场. 用该位移场构造了黏弹性界面裂纹“加料” 单元和过渡单元位移模式,推导了增量“加料” 有限元方程,求解有限元方程可获得应力强度因子和应变能释放率等断裂参量. 建立了典型黏弹性界面裂纹平面问题“加料” 有限元模型,计算结果表明,对于弹性/黏弹性界面裂纹和黏弹性/黏弹性界面裂纹,该方法都能得到相当精确地断裂参量,并能很好地反映蠕变和松弛特性,可推广应用于黏弹性界面断裂问题的计算分析.      

Local crack tip strain behavior during fatigue crack initiation and propagation in polycarbonate

To investigate whether the two fatigue processes of crack initiation and propagation can be combined, the change of local notch root strain and its history, as well as the change of local crack tip strain and the local strain history, of a fatigued element ahead of the propagating crack tip up to failure in a polycarbonate subjected to low-cycle fatigue tests are measured by the fine grid method. As a result, the existence of a unified local strain field in which the two fatigue processes can be substantially combined is experimentally confirmed. Therefore, the local crack tip strain may be examined by a simpler, one-parameter approach for fatigue life estimation.