实验模态分析的基本原理与主要过程
基本概念
模态分析
模态分析是研究结构动力特性一种方法,一般应用在工程振动领域。
其中,模态是指机械结构的固有振动特性,每一个模态都有特定的固有频率、阻尼比和模态振型。分析这些模态参数的过程称为模态分析。
按计算方法,模态分析可分为: * 计算模态分析:由有限元计算的方法取得待识别参数 * 试验模态分析:通过试验将采集的系统输入与输出信号经过参数识别获得待识别参数
通常,模态分析都是指试验模态分析。
实验模态分析
实验模态分析主要是通过模态试验,测量系统的振动响应信号或同时测量系统的激励信号和响应信号,从测量信号中,识别描述系统动力特征的有关参数。
试验模态分析结果可以验证和修改结构数学模型、综合试验分析模型、校核动态计算分析结果的有效性、分析结构产品性能、检查结构中的薄弱部位、鉴定结构总体设计的合理性。
待识别参数
- 物理参数:质量矩阵、刚度矩阵、阻尼矩阵等
- 模态参数:固有频率、衰减系数、模态矢量、模态刚度、模态阻尼
模态分析系统的组成
- 激振系统:使得系统产生稳态、瞬态或随机振动
- 稳态正弦激振法
- 瞬态激振法
- 快速正弦扫描激振法
- 脉冲锤击激振法
- 随机激振法
- 纯随机激振法
- 伪随机激振法
- 周期随机激振法
- 测量系统:用传感器测量实验对象的各个主要部位上的位移、速度或加速度信号,然后将这些信号与激振信号一起记录并保存
- 分析系统:将保存的激励信号和响应信号经过数模转换,采样并输入数字式分析仪或计算机中,用硬件或软件系统识别振动系统的模态参数
实验模态分析的基本原理
基本原理
本部分以多自由度动力学系统为例,阐明了实验模态分析的基本原理: 对于多自由度动力学常微分方程 \(\rightarrow\) 对其拉氏变换 \(\rightarrow\) 得到传递函数 \(\rightarrow\) 通过振型的正交性 \(\rightarrow\) 发现传递函数中含有模态参数的信息 \(\rightarrow\) 因此可以通过实验可以得到传递函数 \(\rightarrow\) 再从传递函数中识别得到振动系统有关参数 \(\rightarrow\) 最后进一步计算得到待求模态参数
传递函数与模态参数的关系
对于多自由度动力学系统,其动力学方程为:
\[ \textbf{M}\bar{\ddot{x} }+\textbf{C}\bar{ \dot{x} }+\textbf{K}\bar{x}=\bar{f} \]
对其进行拉氏变换:
\[ \textbf{M}\left[ s^2 \bar{X}(s) - s\bar{x}(0) - \bar{\dot{x} }(0) \right] + \textbf{C} \left[ s\bar{X}(s) - \bar{x}(0) \right] + \textbf{K} \bar{X}(s) = \bar{F}(s) \]
初始条件全为0时:
\[ \left[ s^2 \textbf{M} + s \textbf{C} + \textbf{K} \right] \bar{X}(s) = \bar{F}(s) \]
根据传递函数的定义,系统的传递函数为:
\[ \textbf{H}(s)=\frac{\bar{X}(s)}{\bar{F}(s)} = \frac{ 1}{ s^2 \textbf{M} + s \textbf{C} + \textbf{K} } \]
机械系统中传递函数的定义:零初始条件下线性系统响应(即输出)量的拉普拉斯变换与激励(即输入)量的拉普拉斯变换之比。
利用振型的正交性,当模态振型矩阵为$ [] $时,模态质量矩阵、模态阻尼矩阵、模态刚度矩阵为:
\[ \textbf{M}_r = [\varphi]^\mathrm{T}\textbf{M}[\varphi] = \mathrm{diag}(m_1 \, \cdots \, m_r \, \cdots \, m_N) \]
\[ \textbf{C}_r = [\varphi]^\mathrm{T}\textbf{C}[\varphi] = \mathrm{diag}(c_1 \, \cdots \, c_r \, \cdots \, c_N) \]
\[ \textbf{K}_r = [\varphi]^\mathrm{T}\textbf{K}[\varphi] = \mathrm{diag}(k_1 \, \cdots \, k_r \, \cdots \, k_N) \]
将上式代入传递函数,可得到:
\[ \textbf{H}(s)= [\varphi] \left[ s^2 \textbf{M}_r + s \textbf{C}_r + \textbf{K}_r \right]^\mathrm{-1}[\varphi]^\mathrm{T} \]
将传递函数矩阵展开可以发现,它的每一行或每一列都包含有所有的模态信息:模态质量矩阵、模态阻尼矩阵、模态刚度矩阵中所有元素以及所有模态矢量。
实验、传递矩阵和模态参数的关系
若固定在一点测量振动响应信号,而不断改变激励信号的作用点,这样就测量出了传递函数的一行;若固定在一点进行激励,而在不同点进行振动响应信号测量,这样就测量出了传递函数的一列。
从上一节的分析中可知,传递函数矩阵的每一行或每一列都包含有所有的模态信息。因此,只要通过实验测得传递函数矩阵的一行或一列再加以分析计算,就可以得到全部的模态信息。
从传递函数中求出模态参数的方法,被称为参数识别方法。一般通用的参数识别方法有两种:图解法和曲线拟合法。
图解法是利用频响函数曲线直接进行模态识别的方法。其优点是简便、直观,缺点是精度差、不能排除外界干扰。
曲线拟合法是用一条连续曲线去拟合一组离散的实验数据,然后利用拟合曲线识别有关参数的方法。该方法常利用图解法所识别的参数作为初始值进行迭代优化计算,直到满足精度为止,从而可以从传递函数中识别得到振动系统有关参数,最后进一步计算得到待求模态参数。
实验模态分析的主要过程
根据实验模态分析的基本原理,现将实验模态分析的主要过程总结如下: 1. 确定实验模型,将实验结构支撑起来(边界条件的确定) 1. 进行模态试验,利用上述激振方法进行激励实验结构(一般用锤击法),记录原点和各测点的激励、响应时间历程曲线 1. 对各测点的时间历程数据进行数字处理,利用传递函数与模态参数的关系和时/频域转换技术求出各测点的传递函数,并组成传递函数矩阵 1. 利用参数识别方法进行参数识别,求得模态参数 1. 进行模态动画展示