路面不平度数值模拟研究报告

路面不平度的数值模拟研究

［摘要］

在汽车设计开发过程中，常需要预测、研究汽车零部件在时域内振动响应，于是在系统参数已知的情况下，需要即需有公路路面的随机不平度数据。本文研究了一种公路路面不平度的数值模拟新方法，即直接对已知路面不平度的功率谱密度经过一系列处理获得路面的不平度值，研究表明所得路面不平度数据的功率谱密度与所要求的准确一致，并且这种方法简洁实用、便于操作。

关键词：功率谱密度；路面不平度；傅立叶变换；采样

1、引言

汽车以一定的速度行驶时，路面的随机不平度通过轮胎、悬架等传递到车身上，并通过座椅将振动传递到人体。当把汽车近似为线性系统处理时，得到了路面不平度功率谱以及车辆系统的频响函数，就可以求出各响应物理量的功率谱，从而可分析车辆振动系统参数对各响应物理量的影响和评价平顺性。然而，汽车振动系统中包括许多非线性元件，如轮胎（有可能离地＞、渐变刚度悬架、液力减振器、橡胶减振块及悬架的干摩擦阻尼等。为获得更准确的结果，特别是在进行振动幅度较大的汽车可靠性等研究时，需采用非线性振动模型⑴。对于非线性系统，线性系统中熟知的叠加原理不再成立，不能直接采用频域方法进行研究，只能在时域中进行研究。另外，最近主动、半主动控制悬架的研究已经了人们充分重视，控制系统的反馈信号是时域信号，所以在进行控制策略研究时，也只能在时域中进行。对于这两类问题，所需的路面激励是时域或空间域信号，而非频域信号。获得路面随机不平度的方法有两种，一种是实验测试，一种是将路面不平度的功率谱密度变换为空间域激励

函数，近年来受到了广泛重视［1-4］。

1984年国际标准化组织在文件ISO/TC108/SC2N67中提出了路面不平度的功率谱密度表达式模型和分等方法。1986年，中国学者在进行了大量研究的基础上，也提出了类似的表达式和分等方法，制订了相应的国家标准，即GB7031- 86《车辆振动输入一

路面平度表示方法》。对于路面不平度空间域（或时域＞内的问题，各国学者进行了大量研究，早期的研究方法有谐波叠加法（或称三角级数合成法＞，该方法的基本思想是将路面不平度表示成大量具有随机相位的正弦或余弦之和。三角级数合成模型适用于模拟具有任意形状的谱密度的平稳随机过程，而且所得结果的样本是连续的，但该模型涉及大量三角函数运算⑵，计算效率低。除了谐波叠加法外，还有积分单位白噪声、滤波器整形白噪声的方法［5］以及利用ARMA模型的方法［2］等。文献⑹利用功率谱密度的逆变换对铁路轨道的随机过程进行了研究。在此文献的基础上，本文对利用GB7031-

86建议的公路路面功率谱密度的拟和表达式进行研究，获得分布所需频率范围内的离散功率谱密度数据，通过计算、分析获得路面不平度的离散傅立叶变换，进而通过傅立叶逆变换得到路面不平度值。通过上述整个过程以及算例进行研究，可知这种方法概念清楚、简单易行，并且利用这种方法得到的路面不平度的功率谱密度可以达到与所需的功率谱密度准确一致。

2由时域信号得到功率谱密度函数

只有了解了如何由时域信号得到其功率谱密度的过程，才能正确地根据所要求的功率谱密度得到时域信号，对于象路面不平度这样的空间域数据也是如此。设

------ 1 是一个各态历经的平稳随机过程，显然它不能满足绝对可积条件：

［X I ，所以回不存在傅立叶变换，为此引入一个辅助函数

⑺＜截尾函数）上：

<1)

满足绝对可积条件，存在傅立叶变换，即

是各态历经的平稳随机过程，于是其均方值为

根据＜2）

式，上式可进一步表达为

<3)

＜3）式中，上1是I的共轭复数。

当亠时，，所以

<4)

由于自相关函数与功率谱密度函数构成傅立叶变换对, 根据均方值与自相关函数之间的关系，可得

<5)

比较以上两式［8］，得到

上式的定义域为_______ I

只在.M I有非零值，若现在将F的•点左移—构成函数F，贝y有_ =［，即实际上

取|出| ,则

可见

L1

所以

1

显然若e仅在［0, T］上取得非零值，则有

以上定义、分析的功率谱密度为双边谱密度，即对的正负值均有定义。在工程实际中，由于厂|无意义，所以常根据的偶函数性质，把负频率范围的谱密度折算到正频率范围内，从而得到单边谱密度丨，即

所以只需求解的情况，便可得到单边谱密度.T ,于是

.二

Y || <6)

若实际中对回采样的时间间隔、采样点数和总采样时间分别为

1 —■，则其离散傅立叶变换的相临两点的频率差，即频率分辨率［9］为

____ | ，则＜6）式相应的离散形式为

进一步可简写为

<7)

以上建立了离散的时间信号与离散形式的功率谱密度之间的关系，按照＜7）式便可得到功率谱密度I 。但是若采样的时间间隔选取的不合理将产生频率

混叠效应［9］，采样总时间不合理将影响频率分辨率，从而产生较大误差

。若一需要的最小频率和最大频率分别为一，根据采样定理＜Shannon Theorem）可知，若一及」不发生频率混叠，应该有

LrJ

另外，频率间隔，即频率分辨率为

最小频率应满足

其中的离散傅立叶变换。

<8

是

I ＜10

由＜8）、＜9）、＜10）式可知总采样时间和采样点数________ I 、应满足

并且还需要

丨X】

对具有N个数据的离散的时域信号____________ I 而言，其离散傅立叶变换也是N 个数据，相临两个数据对应的频率差为，或丨的第.个点＜回）对

应的频率为卜］，而不发生频率混叠时而的最大频率为匚’，可见声的取值应为：=0, 1, 2,……，—。

以上讨论了如何由离散的时域信号得到其功率谱密度丨的过程，若_____________ I 不是时域信号而是路面不平度值，其功率谱密度的求解方法与上述

完全相同，只不过是上述的时间T对应距离L,采样时间间隔对应采样的距离间隔，时间频率应对应空间频率。

3由路面不平度的功率谱密度得到路面不平度

大量的实验测量表明，路面的不平度是具有零均值的、各态历经的平稳Gaus

s随机过程［2、5］。并且，通常用功率谱密度来描述路面的统计特性，路面的不平度

的垂直位移功率谱密度可用下式来拟合［10］

式中：一空间频率，它是波长泊勺倒数，表示每M长度中包含波的周期数，单位为」「参考空间频率，=0.1」