疲劳分析的数值计算方法及ANSYS疲劳分析实例解读

第十四章疲劳分析的数值计算方法

及实例

第一节引言

零件或构件由于交变载荷的反复作用，在它所承受的交变应力尚未达到静强度设计的许用应力情况下就会在零件或构件的局部位置产生疲劳裂纹并扩展、最后突然断裂。这种现象称为疲劳破坏。疲劳裂纹的形成和扩展具有很大的隐蔽性而在疲劳断裂时又具有瞬发性，因此疲劳破坏往往会造成极大的经济损失和灾难性后果。金属的疲劳破坏形式和机理不同与静载破坏，所以零件疲劳强度的设计计算不能为经典的静强度设计计算所替代，属于动强度设计。随着机车车辆向高速、大功率和轻量化方向的迅速发展，其疲劳强度及其可靠性的要求也越来越高。近几年随着我国铁路的不断提速，机车、车辆和道轨等铁路设施的疲劳断裂事故不断发生，越来越引起人们的重视。疲劳强度设计及其研究正在成为我国高速机车车辆设计制造中的一项不可缺少的和重要的工作。

金属疲劳的研究已有近150年的历史，有相当多的学者和工程技术人员进行了大量的研究，得到了许多关于金属疲劳损伤和断裂的理论及有关经验技术。但是由于疲劳破坏的影响因素多而复杂并且这些因素互相影响又与构件的实际情况密切相关，使得其应用性成果尚远远不能满足工程设计和生产应用的需要。据统计，至今有约90%的机械零部件的断裂破坏仍然是由直接于疲劳或者间接疲劳而引起的。因此，在21世纪的今天，尤其是在高速和大功率化的新产品的开发制造中，其疲劳强度或疲劳寿命的设计十分重要，并且往往需要同时进行相应的试验研究和试验验证。

疲劳断裂是因为在零件或构件表层上的高应力或强度比较低弱的部位区域产生疲劳裂纹，并进一步扩展而造成的。这些危险部位小到几个毫米甚至几十个微米的范围，零件或构件的几何缺口根部、表面缺陷、切削刀痕、碰磕伤痕及材料的内部缺陷等往往是这种危险部位。因此，提高构件疲劳强度的基本途径主要有两种。一种是机械设计的方法，主要有优化或改善缺口形状，改进加工工艺工程和质量等手段将危险点的峰值应力降下来；另一种是材料冶金的方法，即用热处理手段将危险点局部区域的疲劳强度提高，或者是提高冶金质量来减少金属基体中的非金属夹杂等材料缺陷等局部薄弱区域。在解决实际工程问题时，往往需要结合运用以上两种方法进行疲劳强度设计和研究。合理地利用各种提高疲劳强度的手段，可以有效地提高构件的疲劳强度或延长其疲劳寿命，并起到轻量化的作用。

关于疲劳问题的研究，基本上可分为疲劳裂纹的形成和扩展机理、规律方面的基础性研究和疲劳强度设计以及提高疲劳强度的有效途径等应用性研究。应用性研究虽然借鉴了基础性研究的成果，但因为它需要考虑更多的实际影响因素，所以它的研究更为复杂和困难。因此相比之下关于疲劳寿命的预测和疲劳强度设计等应用性研究要少得多，远远落后于实际工程的需求。过去，疲劳强度设计和寿命预测的研究主要是以试验为基础进行的。随着计算机应用技术和有限元数值计算理论及其应用的迅速发展，现在又兴起了基于大量试验数据的疲

劳强度设计和寿命计算的有限元数值计算方法，有力地推动着零部件疲劳强度设计的研究及应用的发展。

当前，业已发展形成了专用的疲劳分析软件，如MSC/FA TIGUE 等。此外许多著名的有限元分析软件也嵌套有功能较为齐全的疲劳强度计算模块，如MARC ，ANSYS ，以及I-Deas 中的CAE 等。这些软件疲劳强度计算模块的细节虽然不尽相同，但是其基本思路与算法大都相似。

本章将阐述疲劳强度设计的基本概念、疲劳强度的主要影响因素、疲劳强度设计的有关理论、基本设计方法等。最后，还将对疲劳强度分析的实例进行介绍。

第二节 疲劳载荷类型与S-N 曲线：

一、疲劳载荷的类型与基本术语

使零件或构件发生疲劳破坏的动载荷称为疲劳载荷，可分为为两类。一类是其大小和正负方向随时间周期性地变化的交变载荷，另一类是大小和正负方向随时间随机变化的随机载荷。交变载荷又称为循环载荷，是最为简单和基本的疲劳载荷形式。所研究结构部位因交变载荷引起的应力称为交变应力。

图14-1（a ）是一个典型的交变应力-时间的变化历程。图中循环应力的大小和正负方（拉压）向随着时间的变化而作周期性的变化。一个周期的应力变化过程称为一个应力循环。应力循环特点可用循环中的最大应力σmax 、最小应力σmin 和周期T （或频率f=1/T ）来描述。因为最大应力和最小应力的绝对值相等而正负号相反，故称这种交变载荷为对称循环应力。典型的循环载荷如圆轴类杆件的旋转弯曲、轴向拉压和平板零件的双向弯曲等，都可以在零件的表面或内部产生这样的交变应力。另外，轴类零件的双向扭转也可以产生类似的交变应力。

图14-1（a ） 对称循环交变载荷

图14-1（b ） 不对称循环交变载荷 在疲劳载荷的描述中经常使用应力幅σa 和应力范围△σ（也称为应力振幅、应力幅度）的概念，定义如下。 （14-1）

（14-2） a σσσσ2=-=∆m in m ax 2

min max a σσσ-=t σ

σ

m a x

σ

m i n △σ σa t

σm σm a x σm i n

σa

σ 0

应力幅σa 反映了交变应力在一个应力循环中变化大小的程度，它是使金属构件发生疲劳破坏的根本原因。

当研究的部位除承受有动载荷外，还有静载分量荷时，动静载荷的共同作用下的应力-时间变化曲线如图14-1（b ）所示。此时的载荷时间-变化曲线相当于把图14-1（a ）的对称循环应力曲线向上平移一个了静应力分量。这种的循环载荷称为不对称循环载荷，并用最小应力与最大应力的比值R 来描述循环应力的不对称程度，R 称为应力比，有时又称为不对称系数，即

（14-4）

由定义可知，当R =-1时的循环应力即为对称循环应力，当R ≠0时统称不对称循环应力。其中，R =0时为拉伸脉动应力，R =-∞时为压缩脉动循环。

循环应力中的静载分量通常称为平均应力，用σm 表示，可由下式求出。

（14-5） 静载分量或平均应力对构件的疲劳强度有一定的影响。压缩平均应力往往提高构件的疲劳强度，而拉伸平均应力往往降低构件的疲劳强度。因此，在疲劳强度和疲劳寿命的研究中，给定一个循环应力水平时，需要同时给出应力幅σa 和应力比R 、或者同时给出最大应力σmax 和平均应力σm ，也有时直接给出最大应力σmax 和最小应力σmin 来表示循环应力水平。

由以上各式可知，在应力幅、平均应力、应力比、最大应力和最小应力的参数中，只要已知其中的两个便可求出其它。如当已知σa 、R 时，其它参数便可由下式得到。 （14-6）

或者已知或σa 、σm 时，

（14-7）

二、材料的S-N 曲线与基本术语

一般情况下，材料所承受的循环载荷的应力幅越小，到发生疲劳破断时所经历的应力循环次数越长。S-N 曲线就是材料所承受的应力幅水平与该应力幅下发生疲劳破坏时所经历的应力循环次数的关系曲线。S-N 曲线一般是使用标准试样进行疲劳试验获得的。如图14-2所示，纵坐标表示试样承受的应力幅，有时也表示为最大应力，但二者一般都用σ表示；横坐标表示应力循环次数，常用N f 表示。为使用方便，在双对数坐标系下S-N 曲线被近似简化成两条直线。但也有很多情况只对横坐标取对数，此时也常把S-N 曲线近似简化成两条直线。

S-N 曲线中的水平直线部分对应的应力水平就是材料的疲劳极限，其原意为材料经受无

max min σσ=R a m a a R R R R R

σσσσσσ-+=-=-=1121212m in m ax 2

min max m σσσ+=a m a m σσσσσσ-=+=m in m ax