1、PDF外文：http:/ 1 中文 7800 字 出处： Journal of System Design and Dynamics, 2009, 3(6): 918-931 采用线性预测模型对铁道车辆车体进行模态分析 Takahiro TOMIOKA， Tadao TAKIGAMI， Ken-Ichiro AIDA 车辆噪声与振动实验室，铁路技术研究所 2 8 38 Hikari cho， Kokubunji-shi， 东京， 185 8540 邮箱： tom
2、iokartri.or.jp 摘 要 本文介绍了一个采用线性预测模型对铁路车 辆进行的模态特性识别研究。由静止或运行试验获得的实际铁路车辆的输入 (激振力或轴箱加速度 )和输出 (车体加速度 )之间的关系可由一个 ARX(自回归线性预测模型 )表示，同时模态参数的提取过程也能被详细描述。一个合适的模型定阶 (即 ARX 模型中预测系数的阶 )应是从实际应用的角度考虑的。分析数据得出两个不同部件的平均估计误差的实现过程被提出，同时他们对决定模型定阶的有效性被评估。使用 ARX 模型得出MIMO(多输入 多输出 )的合适性也被描述。结果表明，利用所提出的方法，详细的模态特
3、性可以被成功地从静止、运行测试测得的数据中 确定。 关键词： 铁路，模态分析，线性预测模型，信号处理，弯曲振动 1.介绍 要提高铁路车辆的行驶质量，重要的是要抑制车体纵向弯曲振动。为抑制这种振动，第一步要做的是对车体频率、模态属性等振动特性进行识别。静止和运行的振动测试通常就是为达此目的而进行的。运行测试在车辆运行时通过一个实际的商业服务性的途径对车辆进行了频率特性的分析和行驶质量的评估。平稳振动测试则适用于确定车体的模态性能，这是因为输入 (激振力 )和输出 (响应加速度 )之间的关系是明确的。 由于进行铁道车辆测试的成本很高
4、，所以通过单一的测量试验同时 评估乘坐质量和模态性能是非常有效的。作者已经介绍了一种从静止测试 )2)(1( 中来评估运行质量的方法。本文主要介绍通过运行测试来评估车体模 2 态性能的技术。 在运行过程中的铁道车辆的输入 /输出关系是复杂和不稳定的。车辆受到多输入的作用，而激励条件会在很短的时间改变。因此运行测试难以确定车体频率和模态性能。 为了应对这一挑战，作者尝试运用类似 ARX(自回归线性预测模型 )的 LPM(线性预测模型 )来分析铁路车辆振动 )7()3( ，因为 LPM 模式对待短时间数据和多输入 多输出 (MIMO)
5、的问题更有效。然而，确定模型定阶 (即 ARX 模型中预测系数的阶 )仍是有问题的。本文介绍的对铁路车辆车体的模态识别采用 ARX 模型，而适当模型定阶的确定则是从实际使用的观点出发。 2. 引入线性预测模型 (LPM)的必要性 在运行过程中，由于赛道条件和运行速度总是在不断变化，铁路车辆的激励条件每时每刻都不相同。因此，车体中诱发弯曲的振动幅度也每时每刻都在变化。 当使用有足够的频率分辨率 F 的 FFT(快速傅立叶变换 )对这样的非定常的振动数据进行分析时，例如用 F = 0.1Hz 来计算加速度的 PSD(功率谱密度 )时
6、，因为缺乏平均我们难以得到可靠的结果。我们能使用足够的数据长度例如 60 秒来解决这个问题。 然而，对于一辆以 300 公里每小时的速度运行的列车，它能在60 秒内行驶 5 公里。这样的话，对于有着显著弯曲振动发生的某一指定部分的数据长度有限的分析， FFT 的方法就不再适合。 作者研究了利用 LPM 来分析铁路车辆振动特性的适用性，并表明了它是很有效的。 LPM 不仅能确定车体模态性能，同时还能确定车体的频率特性 ，例如运行时的加速度的 PSD )6)(5)(4)(3( 。通常，铁路车辆在八个轮子上运行，因此它在运行过程中受到八个垂直方向的激励。 LP
7、M 模型可以被很容易地扩展用以容纳多个输入，这是另一种比运用 FFT 来进行模态分析优越的地方。 3.运用 LPM 进行模态分析 本节概述了利用 LPM 对铁道车辆车体的模态分析，分析中将铁路车辆看成一个多输入多输出 (MIMO)系统。这里分析过程的叙述是基于以前发表的文献 )9)(8( 。 3 3.1 预测系数的计算 )8( 假设输入信号为 u(n)，输出信号为 y(n)，这些信号形成一个任意采样时间 t 下的离散数据序列。在这里， n 表示数据样本的数目。现在，我们利用
8、样品数为 m的过去的输入和输出数据乘以加权系数对数据样本数为 n的输出信号进行预测。如下： MmMmMmMm mnubmnyany1)(1)( )()()( ， (1) 表示预测误差。我们获得如下 方程： )()()()( )(1 )(1 )( nmnubmnyany MyMmMmMmMm ，
9、 (2) 这个等式表明了 M 阶 ARX(自回归线性预测模型 )中输入输出信号的关系。输入输出信号是矢量 u(n)=u1(n)， .， up(n)T 和 y(n)=y1(n)， .， yQ(n)T ， P和 Q 输入输出的级数， ， T 表示向量的转置。 )(Mma和 )(Mmb表示 Q Q 和 Q P 阶矩阵。 接下来我们写出反向形式的输入输出等式 : )()()()( )()(1)(1nmnycmnudnu MuMmMmMmMm ， (3) 这里md和mc表示 P P
10、和 P Q 阶的预测系数的矩阵。 )(Mu(n)表示 u(n)的预测误差。由于等式 (2)和 (3)是独立的，因此它们能被合并。我们得出以下等式。 MmMMm nmnxAnx1)()( )()()( ， (4) 这里 x(n)=uT (n) yT (n)T 是合并输入输出二矢量一系列的时间数据得到的，)()( nM = )(Mu T (n) )()( nTMy T 表示合并的预测误差， )(MmA 表示以下的包含预测系数的分块矩阵： )()()()()(MmMmMmMmMmabcdA 。 方程 (4)表示 x(n)可以表示为有 P + Q 模型独立变量的 AR(自回归 )模型， )(MmA能够运用现有的普通 AR 模型计算得出。在这项研究中，我们采用 Burg 法这被认为
