路面文件的生成

第3章路面文件的生成

3.1 随机路面不平度的拟合理论[4][5][17][18][39]

大量的测量分析结果表明，路面不平度具有随机、平稳和各态历经的特性，可以用平稳随机过程理论来分析描述。通常把道路垂直纵断面与道路表面的交线作为路面不平度的样本，通过样本的数学特征――方差或功率谱密度函数来描述路面。均值为零时，方差可以反映路面不平度大小的总体情况；功率谱密度函数能够表示路面不平度能量在空间频域的分布，它说明了路面不平度或者说路面波的结构。当功率谱密度用坐标图表示时，坐标上功率谱密度曲线下的面积就是路面不平度方差。从功率谱密度函数不仅能了解路面不平度的结构，还能反映出路面的总体特征。因此，功率谱密度函数（PSD）是路面不平度的最重要数学特征。

文献[17、39]介绍了采用多种方法生成路面不平度的时域模型，如：滤波白噪声生成法（线性滤波法），基于有理函数PSD模型的离散时间随机序列生成法，根据随机信号的分解性质所推演的谐波叠加法（也称频谱表示法），以及基于幂函数功率谱的快速Fourier反变换生成法等。

白噪声激励模拟的基本思想是：将路面高程的随机波动抽象为满足一定条件的白噪声, 然后经一假定系统进行适当变换而拟合出路面随机不平度的时域模型。

离散时间随机序列生成法的基本思想是：基于Parkhilovsk ii 提出的另一种有理函数形式的功率谱密度表达式[17]，建立路面不平度时间离散化模拟的递推公式[17]。

谐波叠加法的基本思想是：随机正弦波（或其他谐波）叠加法采用以离散谱逼近目标随机过程的模型, 是另一种离散化数值模拟路面的方法。

基于幂函数功率谱的快速Fourier反变换生成法的基本思想是：由功率谱密度的离散采样构造出频谱, 然后对频谱进行Fourier逆变换得到时域模拟的轨道

不平顺激励函数。

但是，不管是标准道路谱还是实测道路谱，其PSD 是路面不平度的一个统计量。因此，对应于测量范围内某一种确定的路面不平度，其PSD 是唯一的；但对于给定的PSD ，其模拟设计的路面不平度并不唯一，也就是说频域模型和时域模型并非一对一的映射，因此从频域模型所得的路面不平度的时域模型只能看成是满足给定路谱的全部可能的路面不平度中的一个样本函数。

在本文所讨论的方法中, 谐波叠加法是目前使用较普遍的方法, 特别适合用于国标道路谱时域模型的生成。尽管谐波叠加法计算量相对较大，但该方法理论基础严密、算法简单直观、数学基础严密、是一种高保真的频域模型转换方法。因此，本文采用谐波叠加法来构建随机路面不平度的频域模型。

谐波叠加法拟合不平路面的原理是：设路面高程为平稳的、遍历的均值为0的Gaussian 过程，则可以用不同形式的三角级数进行模拟。以下以正弦波为例加以描述。随机正弦波（或其它谐波）叠加法采用以离散谱逼近目标随机过程的模型，是一种离散化数值模拟路面的方法。随机信号，比如随机风压或路表轮廓，可以通过离散Fourier 分析变换分解为一系列具有不同频率和幅值的正弦波。谱密度就等于由带宽划分的这些正弦波幅值的平方。

按照国标GB7031-86《车辆振动输入――路面平度表示方法》，路面不平度位移功率谱密度拟合表达式采用下式：

()()00w d d n G n G n n -⎛⎫= ⎪⎝⎭ n >0 （3－1）

式中：0n —参考频率，100.1m n -=；()0d G n —路面不平度系数，21m /m -；()d G n —路面不平度，21m /m -；w —频率指数，经验值w ＝2。

已知在空间频率12n n n <<内的路面位移谱密度为()d G n ，利用平稳随机过

程的平均功率的频谱展开性质，路面不平度的方差2z σ为

()212

n z d n G n dn σ=⎰ （3－2）

将区间（1n ，2n ）划分为n 个小区间，取每个小区间的中心频率

()1,2,mid i n i -=…,n 处的谱密度值()d mid i G n -代替()d G n 在整个小区间内的值，则

式（3－2）离散化后近似为

()2

1n

z d

m id i i i G n n σ-=≈⋅∆∑ （3－3）

对应每个小区间，现在要找到具有频率()1,2,m i d i n i n -=…,且其标准差为

()sin 2mid i i n x πθ-+ （3－4）

将对应于各个小区间的正弦波函数叠加起来，就得到频域路面随机位移输入

()

()1sin 2n

m id i i i q x n x πθ-==+∑

（3－5）

式中：θ—[0，2π]上均匀分布的随机数；x —频域路面的X 方向

自此，通过离散频域路面的X 方向的x 值，就可以得到空间频域下随机路面Y 方向的值。

理论是这样的，然后按照李军的《ADAMS 实例教程》132页的格式编制路面文件。上面的这些理论X 表示前进方向，Y 表示路面的不平度（高度），Z 表示路面的宽度。上面的公式只是X 和Y 的关系，也就是一个二维路谱。Z 的大小任意定无所谓。

有的理论是把这样的路谱横向并列放，就是好多笔平行放一样，这样形成的路谱只是伪三维路谱，并不是真正的三位路谱。真正的三维路谱是Z(x,y)＝f(x,y)这样的关系。

3.2 随机路面的生成

在ADAMS 中，不平的路面是由一系列三角形的平面单元组合成的一个三维表面。原理图如图3-1所示，数字1、2、3等等表示节点（Node ），这些节点的x 、y 坐标要满足一定的规律，z 坐标仅仅表示路面的宽度；由这些节点按一定的规律组成路面单元（Element ）；再在路面单元里设置静摩擦系数和动摩擦系数，就能

模拟真实的路面。

图3-1 随机不平路面原理图

而且，随机路谱的编制要满足轮胎的要求：第一，路谱的位置要处于轮胎的下方；第二，路谱向上的方向要指向轮胎所处的一侧；第三，路谱的大小要根据仿真的需要确定。

按照上述路面原理图和上一节的不平路面拟合理论，编制出随机路面VB程序（附录C），程序的流程图如图3-2所示，程序界面如图3-3所示。