第1页共6页编号：时间：2021年x月x日书山有路勤为径，学海无涯苦作舟页码：第1页共6页采用工作模态识别法进行模态分析郑敏,申凡,史东锋,陈怀海.中国航空学报(英文版),2004,17(4):215–219摘要：对于一些大型工作结构,有时很难测得输入信号,只能单独利用实测响应数据进行工作模态识别。本文针对工程实际中的这一情况,首先利用响应间互相关函数同脉冲响应函数在表达形式上的相似性,推导出了多参考点时域工作模态复指数法,继而又从响应间的互相关函数入手,推导出了多参考点频域工作模态识别法。最后采用一飞机模型对这两种方法进行了试验验证,并将所得结果做了对比分析,结果表明,两种多参考点方法都能较有效地单独从实测响应数据提取结构的模态参数,且频域工作模态识别法比时域工作模态复指数法识别精度更高。关键词：模态识别;工作状态;系统识别;自然激励;环境激励工作条件下的航空发动机及核反应堆，飞行中的导弹和航空器，受风或车辆载荷激励的大型建筑输入力的测量是不可实现的。在实验室中通过人为激励进行模拟实验是非常困难的。从操作机构中提取模态参数能够极大地降低实验花费。因此，最近几年在单响应数据中模态鉴定系统得到越来越多的关注。已经有一些方法从单输出数据中估算模态参数，例如功率谱密度函数的峰值选取法，自回归滑动模型，随机子空间法，用易卜拉欣时间域技术随机减量处理，最大熵法，最小二乘曲线拟合技术等等。本文，将推导出在白噪声激励下测量反馈间互相关函数表达式。接着提出了一个多参考时域的的操作模型复指数法。1理论方面1.1相关函数及脉冲响应函数假设输入时白噪声，那么在输出n及输出p之间的互相关函数就可以获得第2页共6页第1页共6页编号：时间：2021年x月x日书山有路勤为径，学海无涯苦作舟页码：第2页共6页N是模态的个数，λr是第r个复特征值。是第r个模态矩阵的第n个成分。Ψr,Qpr是常数。众所周知，脉冲响应函数被描述如下：是有点l导致的n点的脉冲响应，是模态参与因数。因此，它表明互相关函数是负指数函数的总和并且与脉冲响应函数具有相同的形式。1.2工作模态复指数方法作为参考服务的N反应与P反应之间自动与互相关函数矩阵表示如下：表示N×2N模态矩阵，表示2N×2N复特征值矩阵，Q表示充满的2N×2N常数矩阵。让模型的阶位为q，下面的方程式可以通过在传统的聚参考复指数方法中用互相关函数矩阵代替脉冲响应函数矩阵。U=V×BU=[,第3页共6页第2页共6页编号：时间：2021年x月x日书山有路勤为径，学海无涯苦作舟页码：第3页共6页T表示矩阵或向量的转置。K是互相关函数的样本数。系数矩阵B（1），B（2），。。。B(q)能够从方程式4中获得。对于第r模态可获得下列方程式其中Qr第r个适量矩阵的Q，Zr=，△为样本区间。一旦知道系数矩阵EQ（1）可以改写成一个产生于特征值Zr和相应的左特征向量Qr的广义特征值问题。通过下式可以获得系统极点λr模态可以通过方程式3获得。1.3频域工作模态识别方法方程式3的拉普拉斯变换可以写为其中G（s）是相关函数矩阵R(T)的拉普拉斯变换，I是单位矩阵。一个新的矩阵方程就构成了G(s)是R(T)的一次导数的拉普拉斯变换；。众所周知，必定存在系统矩阵A满足下列关系:第4页共6页第3页共6页编号：时间：2021年x月x日书山有路勤为径，学海无涯苦作舟页码：第4页共6页将上述方程式自右乘O(s)可以获得下式将方程式8代入方程式10可以获得也可以为方程式12也可以表示为将方程式13代入方程式11，并使s=jw，可以获得下式表示测量响应的功频率谱密度函数矩阵。方程式14可以表述为对于在测量的频率范围内所有离散频率ω1,ω2,...,ωK都可以获得一组标量方程。这些方程也可以写为。从方程式16中可以获得：其中+表示广义逆。当测量的功率谱密度已知时可以推导出矩阵A。求解方程式9的特征值在于从Λ对角线中得到系统特征值λr。原有第5页共6页第4页共6页编号：时间：2021年x月x日书山有路勤为径，学海无涯苦作舟页码：第5页共6页振动系统的第r个模态仅仅是矩阵的第r个特征值2应用飞机模型由三个紧固在飞机机身上三个点的橡胶绳悬置在空中。向两个振动器中施加两个互不相关的白噪声信号。MVMAS23多振动测量和分析系统有64的通道能...
