某差速器齿轮的动态接触仿真与疲劳分析

为分析差速器齿轮的疲劳寿命,运用动态有限元与试验相结合的方法,研究了锥齿轮在啮合过程中的应力分布以及疲劳强度。首先基于CATIA软件对差速器的半轴齿轮、行星齿轮进行了参数化建模,并采用动态有限元法模拟了齿轮副在最大扭矩工况下的动态接触;再根据有限元强度分析结果,将最大接触应力作为静载输入,运用疲劳分析软件对齿轮副的接触疲劳性能进行计算;最后进行了差速器齿轮副的台架试验,并将仿真结果与台架试验进行了对比。结果表明:齿面最大接触应力产生于节圆附近,齿轮间的最大接触应力为1309MPa;半轴齿轮在90%存活率下的疲劳寿命为3.394×106;仿真结果与台架试验具有较好的一致性,齿轮满足疲劳寿命要求。将动态有限元和疲劳寿命分析方法相结合可以有效预测差速器齿轮疲劳寿命。     

基于矩法的应力疲劳寿命可靠性曲线分析方法

通过对Weibull分布、正态分布、指数分布等常用概率模型函数的中值疲劳寿命和高可靠度疲劳寿命关系的分析,建立了应力疲劳寿命可靠性通用表达式.指出疲劳寿命服从各常用分布的可靠性公式都可看作应力疲劳寿命通用公式的特例.基于概率模型分布函数的矩法理论,提出了适用于5种常用概率模型分布的应力疲劳寿命可靠性曲线分析方法.应用基于矩法的应力疲劳寿命可靠性曲线分析方法,获得了42CrMo硬齿面齿轮轮齿弯曲疲劳试验寿命数据的中值疲劳寿命和高可靠度疲劳寿命曲线方程.     

一种FRP累积损伤模型及其在结构疲劳寿命估算中的应用

推荐了一种应变损伤累积模型，能够考虑单向板面内多轴应力和平均应力的影响。只需要单向板在确定应力比下的若干典型疲劳试验结果，就可以预测相同材料体系多向层压结构在不同应力比的循环载荷下的疲劳寿命，有助于降低试验成本和工作量。研究了适用于多向层压结构剩余强度估算和疲劳寿命预测的步骤和程序。针对碳纤维／树脂基T300／QY8911复合材料，试验测定了三组典型单轴循环应力（[0]16拉－拉、[90]16拉－拉和[0/90]4S剪－剪）下的S－N曲线。以此为输入，预测四种多向铺层板在各种拉－拉循环应力下的疲劳寿命，寿命预测结果和相应的试验结果吻合良好。采用了保持计算和试验的载荷／强度比相对等值的方法来近似抵消层合效应对疲劳寿命的影响。强调了进一步发展能够定量估计层间应力影响与分层扩展过程的疲劳损伤模型的重要性。     

齿轮参数对三维裂纹应力强度因子影响的研究

建立了齿轮三维精确三齿模型,结合线弹性疲劳力学和边界元方法,运用专业断裂分析软件FRANC3D对齿轮齿根裂纹进行了数值分析,模拟求解了齿根椭圆形裂纹的裂尖三种类型的应力强度因子,并且探讨了三种应力强度因子随载荷、裂纹长度、模数、齿数、变位系数、裂纹角度等齿轮参数的变化规律。结果表明:I型、II型应力强度因子的数值在裂纹前缘呈近似于抛物线分布,而III型应力强度因子的数值在裂纹前缘呈近似于正弦曲线分布;应力强度因子随齿面载荷的变化呈线性规律;模数对裂纹两端的II型应力强度因子、裂纹中间的I型和III型应力强度因子影响较小,基本不变;齿数对应力强度因子影响较小,其在不同齿数下的最大差值仅为187 N.m-3/2;变位系数对其影响基本也呈线性变化;裂纹发生角度为60°时,I型应力强度因子比其它角度下的值都要大,因此应尽量避免产生此角度下的裂纹。本文研究为齿轮的断裂分析和寿命预测奠定了一定的基础。     

基于修正的局部应力应变法估算连接件疲劳寿命

飞机机体结构疲劳裂纹的萌生与扩展几乎都发生在紧固件连接处,预测结构连接件的疲劳寿命具有重要意义.为了准确计算连接件疲劳寿命,提出了一种修正的局部应力应变法,该方法首先采用应力严重系数法和修正Neuber法分析连接件高应力区的应力、应变,然后利用Manson-Coffin方程计算疲劳寿命,此外考虑到在中、高周疲劳里表面加工和尺寸因素的影响是不能忽略的,对应变-寿命曲线的弹性段作了修正以使其同时适用于高、低周疲劳寿命的估算.以搭接件为算例进行疲劳寿命计算,并使用MTS Landmark试验机对搭接件进行疲劳试验.结果表明,采用此方法估算的连接件疲劳寿命与试验结果相比误差小于16％,证明了此方法的有效性.     

表面裂纹疲劳扩展形状演化预测

应用新的应力强度因子经验公式预测了拉伸和纯弯曲载荷下疲劳扩展过程中表面裂纹形状的演化,并与Newman和Raju的经验公式以及实验结果进行了比较。结果表明:虽然对应力强度因子估算有一定差别,但两者对裂纹形状演化的预测有很好的一致性,预测结果与实验结果也比较吻合,从而表明所提出的经验公式能够用于预测疲劳扩展过程中表面裂纹形状演化。     

残余应力下厚壁筒表面裂纹的应力强度因子计算

本文首先介绍了边界元法计算裂纹尖端应力强度因子的基本理论,接着利用边界元法计算了在残余应力下不同厚壁筒内表面椭圆裂纹的应力强度因子,研究了其大不随椭圆裂纹不同而变化的规律,为厚壁筒结构的设计,制造以及疲劳寿命分析提供了许多有价值的参考资料。     

扩展有限元法在结构件疲劳寿命估算中的应用

为了高效、准确估算飞机典型结构件的疲劳扩展寿命,基于通用有限元分析软件ABAQUS的扩展有限元功能计算应力强度因子,采用NASGRO方程计算疲劳裂纹扩展速率。为了避免大量重复提交运算的繁琐过程,应用PYTHON语言开发了在ABAQUS疲劳裂纹扩展分析中反复调用扩展有限元结果和NASGRO方程的子程序,在指定的疲劳裂纹扩展路径上实现了对飞机连接件疲劳裂纹扩展速率与寿命的计算。结果表明,模拟的疲劳裂纹稳定扩展阶段的寿命曲线与实验结果吻合较好,疲劳寿命误差约为10%。     

考虑平均应力的同频拉扭多轴高周疲劳寿命评价方法

目前基于临界平面理论的高周疲劳寿命预测模型, 大都充分考虑了法向平均应力对材料疲劳寿命的影响, 但是没有有效反映剪切平均应力对疲劳寿命的影响. 通过分析7075-T651铝合金的试验数据发现, 与法向拉平均应力类似, 剪切平均应力同样对材料的疲劳寿命产生不利影响. 因此, 如果寿命预测模型中忽略剪切平均应力的影响, 存在明显剪切平均应力加载工况下, 预测寿命可能偏于危险. 由此, 本文定义具有较大法向应力的最大剪应力范围平面为临界平面, 建立了一个能够同时反映法向和剪切平均应力影响的高周疲劳寿命预测模型, 并给出了模型中材料常数的确定方法. 新模型首先将基于应变的Fatemi-Socie准则, 推广到材料的高周疲劳寿命预测, 给出了Fatemi-Socie准则的应力表述形式. 然后, 引入剪切和法向Walker因子, 反映剪切和法向平均应力对材料疲劳寿命的影响. 剪切和法向Walker因子的取值都介于0和1之间, 不同取值反映了材料对剪切和法向平均应力敏感程度的不同. 新模型适用于范围内的金属塑性材料. 利用5种材料在12种存在平均应力加载工况下的试验数据, 对所建模型进行了试验验证, 结果表明预测结果与试验结果吻合良好, 绝大多数寿命预测结果分布在3倍误差带以内.      

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

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

A MULTIAXIAL HIGH CYCLE FATIGUE LIFE PREDICTION MODEL CONSIDERING THE EFFECT OF MEAN STRESS FOR TENSION-TORSION LOADINGS WITH SAME FREQUENCY 1)

目前基于临界平面理论的高周疲劳寿命预测模型, 大都充分考虑了法向平均应力对材料疲劳寿命的影响, 但是没有有效反映剪切平均应力对疲劳寿命的影响. 通过分析7075-T651铝合金的试验数据发现, 与法向拉平均应力类似, 剪切平均应力同样对材料的疲劳寿命产生不利影响. 因此, 如果寿命预测模型中忽略剪切平均应力的影响, 存在明显剪切平均应力加载工况下, 预测寿命可能偏于危险. 由此, 本文定义具有较大法向应力的最大剪应力范围平面为临界平面, 建立了一个能够同时反映法向和剪切平均应力影响的高周疲劳寿命预测模型, 并给出了模型中材料常数的确定方法. 新模型首先将基于应变的Fatemi-Socie准则, 推广到材料的高周疲劳寿命预测, 给出了Fatemi-Socie准则的应力表述形式. 然后, 引入剪切和法向Walker因子, 反映剪切和法向平均应力对材料疲劳寿命的影响. 剪切和法向Walker因子的取值都介于0和1之间, 不同取值反映了材料对剪切和法向平均应力敏感程度的不同. 新模型适用于范围内的金属塑性材料. 利用5种材料在12种存在平均应力加载工况下的试验数据, 对所建模型进行了试验验证, 结果表明预测结果与试验结果吻合良好, 绝大多数寿命预测结果分布在3倍误差带以内.     

Experimental investigation of dynamic punch tests on isotropic and composite materials

An experimental investigation is conducted on the two-dimensional punch problem for isotropic materials and unidirectional fiber-reinforced composite materials under quasi-static and impact loading. Singular stresses are generated in the specimen near the punch corners, and the stress intensity factorK _Iis introduced to describe the singular stress field. Laser interferometry is used to measure in-plane stresses (transmission mode) and out-of-plane displacements (reflection mode) and then estimate the stress intensity factor. In the dynamic case, a high-speed photography technique was employed to capture the transient response of the specimen and measureK(t) just after the impact. In all the cases, a good agreement between the measurements ofK and theoretical predictions was found.     

Application of effective stress intensity factor crack closure model to evaluate train load sequence effects on crack growth rates

An effective stress intensity factor crack closure model is studied in relation to results of laboratory spectrum crack growth tests on compact tension specimens (CTS) fabricated from rail steel. Comparison of model predictions with test results for crack growth life is effected by means of an analysis of a center-cracked tension panel (CCT) subjected to an equivalent stress spectrum. The trends of the model predictions and test results agree as to the effect of changing cycle order in the spectrum, but the actual effect on crack growth life in the laboratory tests is found to be much stronger than the effect predicted by the crack closure model using the effective stress intensity factor.     

Photoelastic analysis of the singular stress field in a bimaterial wedge

This paper presents an experimental investigation of the singular stress field near the vertex of a bimaterial wedge using a digital photoelastic technique. Special attention is given to the casting of bimaterial wedge specimens and analysis technique for extracting stress intensity factors from photoelastic samples. Different bimaterial wedge specimens are made of two different photoelastic materials bonded through a special casting procedure and loaded in simple tension. A new multiple-parameter method is developed to obtain the stress intensity factor reliably from the isochromatic fringe patterns and the series representation of the stress field at the vertex of the wedge. Experimental results are compared with finite element predictions, and good agreement is observed.     

Birefringence measurements on polymer melts in an axisymmetric flow cell

 The stress-optical rule relates birefringence to stress. Consequently, measurement of flow birefringence provides a non-intrusive technique of measuring stresses in complex flows. In this investigation we explore the use of an axisymmetric geometry to create a uniaxial elongational flow in polymer melts. In axisymmetric flows both birefringence and orientation angle change continuously along the path of the propagating light. The cumulative influence of the material's optical properties along the light's integrated path makes determination of local birefringence in the melt impossible. One can nevertheless use birefringence measurements to compare with predictions from computer simulations as a means of evaluating the constitutive equations for the stress. More specifically, in this investigation we compare the light intensity transmitted through the experimental set-up vs entry position, with the theoretically calculated transmitted intensity distribution as a means of comparing experiment and simulation. The main complication in our experiments is the use of a flow cell that necessarily consists of materials of different refractive indices. This introduces refraction and reflection effects that must be modeled before experimental results can be correctly interpreted. We describe how these effects are taken into account and test the accuracy of predictions against experiments. In addition, the high temperatures required to investigate polymer melts mean that a further complication is introduced by thermal stresses present in the flow cell glass. We describe how these thermal-stresses are also incorporated in the simulations. Finally, we present some preliminary results and evaluate the success of the overall method. Received: 2 April 2001 Accepted: 27 August 2001     

Stress-strain state of an inclined elliptical defect in a plate under biaxial loading

A mathematical model for the stress-strain state of a plate with an inclined elliptical defect under biaxial loading is considered. Exact formulas for stresses in polar coordinates, displacements, principal stresses, maximum shear stress, and stress intensity in the case of a plane stress state of the plate were obtained by the Kolosov-Muskhelishvili method. Simulation results are compared with experimental data obtained by holographic interferometry.     

Photoelastic validation of a theoretical solution for an edge cracked disk

In this paper, the theoretical solution developed by Vaughan and Wu for the stress analysis of a circular disk with a radial edge crack extending to its center is validated by photoelasticity. The photoelastic results include the fullfield isochromatics as well as measurements of the maximum shear stress at a number of test points. Additionally, the experimental stress intensity factor is extracted from the photoelastic data by Irwin's two-parameter method. Good agreement is observed when the theoretical stress field and stress intensity factor are compared with the experimental results. It is concluded that the Vaughan-Wu solution can be confidently applied in fracture mechanics analyses. The application of such a solution to the stress analysis of two-dimensional bodies with complex geometries subjected to complex loading is also noted.     

Geometric functions of stress intensity factor solutions for spot welds in lap-shear specimens

In this paper, the stress intensity factor solutions for spot welds in lap-shear specimens are investigated by finite element analyses. Three-dimensional finite element models are developed for lap-shear specimens to obtain accurate stress intensity factor solutions. In contrast to the existing investigations of the stress intensity factor solutions based on the finite element analyses, various ratios of the sheet thickness, the half specimen width, the overlap length, and the specimen length to the nugget radius are considered in this investigation. The computational results confirm the functional dependence on the nugget radius and sheet thickness of the stress intensity factor solutions of [Zhang, S., 1997. Stress intensities at spot welds. International Journal of Fracture 88, 167–185; Zhang, S., 1999. Approximate stress intensity factors and notch stresses for common spot-welded specimens. Welding Journal 78, 173s–179s]. The computational results provide some geometric functions in terms of the normalized specimen width, the normalized overlap length, and the normalized specimen length to the stress intensity factor solutions of [Zhang, S., 1997. Stress intensities at spot welds. International Journal of Fracture 88, 167–185; Zhang, S., 1999. Approximate stress intensity factors and notch stresses for common spot-welded specimens. Welding Journal 78, 173s–179s] for lap-shear specimens. The computational results also indicate that when the spacing between spot welds decreases, the mode I stress intensity factor solution at the critical locations increases and the mode mixture of the stress intensity factors changes consequently. Finally, based on the analytical and computational results, the dimensions of lap-shear specimens and the corresponding approximate stress intensity factor solutions are suggested.