基于加速度信号的振动疲劳分析方法研究
作者：王宗乐房芳
摘要：本论文首先介绍了结构静态疲劳和振动疲劳的区别和基本原理，以及利用加速度信号进行振动疲劳分析的流程和方法。阐述了振动疲劳计算的理论和公式，说明汽车结构振动疲劳寿命预测的工作流程,文中最后对某支架的振动疲劳寿命进行了仿真。
主题字：随机振动；频响函数，振动疲劳
引言
疲劳耐久是汽车产品开发的一项重要指标，同时也是影响产品品牌价值和产品顾客感知的重要依据。疲劳耐久性指标的设计和开发又是一件既费时又耗钱的工作，传统意义上的耐久性开发工作主要包括三个过程：设计-验证-改进-再验证-批准。这种开发方法就是一种试错（trial error）的方法，即使取得好的效果，也会浪费产品开发时间，降低产品的市场竞争力。
计算机辅助工程（CAE）最大的优势就是提供了一种全新的思路和工具，可以在产品物理样机之前进行汽车产品疲劳耐久性的性能设计和分析，从而减少耐久性性能开发的时间和成本。通常来说，利用计算机辅助工程进行耐久性性能开发有两种方法：一种是时域方法，另一种是频域方法。所谓时域方法就是利用测试仪器，通常是应变片获取结构随着时间变化的载荷数据，然后利用线性叠加的方法获取确定的应力时间历程，根据结构模型的特征和材料的疲劳特性，利用miner 法则计算结构的时域疲劳寿命和分布。时域方法的优点是计算简单，方法成熟，能够在大结构疲劳问题上取得满意的结果；同时时域方法的最大缺点就是提前假设载荷信号是确定性、周期性的，并且忽略了结构对外界载荷的动态响应。频域方法对载荷信号的要求比较简单，只要获取相应的加速度信号，并且载荷和响应都是以功率谱密度（PSD）信号形式给出。与时域信号相比，信号容量更少，计算效率明显提高；频域方法同时考虑了路面作用于汽车结构的载荷随机性，并且兼顾了结构对外界载荷的动态响应，这也是频率方法计算汽车结构疲劳寿命优于时域方法的重要方面。
文章对频域内的疲劳耐久性分析方法进行了研究，通过分析频域疲劳方法的基本原理和基本方法，形成了针对汽车结构疲劳寿命研究的固定流程。进行频域疲劳耐久性分析目前在汽车行业应用较少，本文的作为一种尝试，共同研究探索适合汽车产品耐久性开发方法和思路。
1 基于加速度信号的振动疲劳分析流程和思路
基于加速度激励信号的频域振动疲劳分析流程和时域法疲劳寿命分析有着明显不同的分析流程和思路，频域分析兼顾了结构对外界载荷的动态响应，更加符合汽车结构的使用情况。经过对于实际项目的分析和不断摸索，我们总结出了利用频域方法开展振动疲劳的基本流程：
图1－1 基于加速度信号的频域疲劳寿命分析流程
2 振动疲劳寿命分析的理论研究
2．1 振动与疲劳的关系
美国材料试验协会（ASTM）在“疲劳试验及数据统计之有关术语的标准定义”（ASTM
E206-72）中定义“疲劳是在某点或某些点承受扰动应力，且在足够多的循环扰动作用之后形成裂纹或完全断裂的材料中所发生的局部的、永久结构变化的发展过程。”很明显结构疲劳与其受到的振动载荷密切相关。工程实际中，绝大多数结构或设备都处于不同程度的振动环境之中，它们直接或间接受到各种不同动态载荷的激励作用，在长时间的振动损伤累积作用下，最终发生疲劳断裂。汽车行驶时要经受发动机产生的振动以及地面不平、紧急刹车等引起的动态载荷，一些安装附件的支架结构容易受到振动信号的作用而发生短期疲劳。
结构的疲劳过程包括裂纹的形成、裂纹扩展和快速断裂三个阶段。由于结构裂纹改变了结构的动力特性，同样也就改变了结构动力特性对疲劳的作用。例如，趋向共振的结构裂纹扩展速率加快，结构疲劳破坏加剧；与此同时，由于疲劳裂纹的扩展，结构的自然频率发生变化，振动激励的频率偏离结构自然频率，结构远离共振状态，因此结构的裂纹扩展速率又会开始下降。
从以上分析，可以发现振动与疲劳密不可分，解决疲劳问题必须同时兼顾结构对外界激励的动态响应特征。以往简单的忽略结构对外界载荷动态特性的疲劳分析方法是不准确的，特别是对于一些受到振动信号作用的支架类零部件，更加需要根据振动激励信号，开展振动疲劳寿命分析和研究。
2．2 振动疲劳寿命计算理论
如同基于时域疲劳分析方法一样，振动疲劳寿命分析的关键一步就是在获取结构的动态响应应力信号PSD 以后，根据应力的PSD 信号，计算结构的疲劳寿命问题。根据上面的分析流程，对获取的信号进行预处理，得到可以通过程序进行疲劳计算的输入数据，然后利用MSC.nastran 软件开展结构频率相应分析，获取结构单位激励信号的响应，利用结构响应分布以及测试的激励信号，配合结构疲劳数据获取疲劳寿命数据。
3 结构频率响应分析
频率响应分析主要研究结构在外界稳态振动激励作用下的响应问题，解决频响分析的数值方法目前主要有2 种：直接方法和模态方法。相比较模态方法，计算简单，所需计算资源较少，这里主要介绍模态分析求解频响分析。
3．1 无阻尼模态频率响应分析
模态法主要是利用结构的模态振形对运动方程进行解耦，缩减矩阵，使得数值求解更有效率。
模态解决此方程的第一步就是利用模态振形把变量从物理坐标系转换到模态坐标系：
上式中[Φ]是表示全部的结构模态振形，也就是利用了结构的所有模态；实际上只是利用了部分振形参与计算，因此上面的公式也是个近似计算转换关系。
首先忽略结构阻尼，并且将坐标转换到模态坐标下面，运动方程两边同事除以e iωt，运动方
程变化为：
到现在运动方程已经完全在模态坐标下了，不过方程仍然是耦合的，无法进行简单的求解。为对上面的运动方程进行解耦，对运动方程左乘[φ]T，这样方程变形为：
其中：
解耦运动方程的最后一步就是利用模态向量的正交性，上面的模态（广义）质量矩阵、模态（广义）刚度矩阵都是对角矩阵。对角矩阵在偏离对角线上的矩阵元素都是零，正是这些非对角元素造成了方程耦合。至此运动方程完全解耦，成为一系列解耦的单自由度矩阵
其中：
mi 第i 阶模态质量
ki 第i 阶模态刚度
pi 第i 阶模态载荷
通过上面单自由度方程计算出了模态响应（模态坐标）ξ(ω)，物理响应利用前面的坐标转换得出：
