路面文件的生成

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

2、平度的最重要数学特征。文献17、39介绍了采用多种方法生成路面不平度的时域模型，如：滤波白噪声生成法（线性滤波法），基于有理函数PSD模型的离散时间随机序列生成法,根据随机信号的分解性质所推演的谐波叠加法（也称频谱表示法），以及基于幂函数功率谱的快速Fourier反变换生成法等。白噪声激励模拟的基本思想是：将路面高程的随机波动抽象为满足一定条件的白噪声,然后经一假定系统进行适当变换而拟合出路面随机不平度的时域模型。离散时间随机序列生成法的基本思想是：基于Parkhilovskii提出的另一种有理函数形式的功率谱密度表达式17，建立路面不平度时间离散化模拟的递推公式17。谐波叠加法的基本思想是：

3、随机正弦波（或其他谐波）叠加法采用以离散谱逼近目标随机过程的模型,是另一种离散化数值模拟路面的方法。基于幕函数功率谱的快速Fourier反变换生成法的基本思想是：由功率谱密度的离散采样构造出频谱,然后对频谱进行Fourier逆变换得到时域模拟的轨道不平顺激励函数。但是，不管是标准道路谱还是实测道路谱，其PSD是路面不平度的一个统计量。因此，对应于测量范围内某一种确定的路面不平度，其PSD是唯一的；但对于给定的PSD，其模拟设计的路面不平度并不唯一，也就是说频域模型和时域模型并非一对一的映射，因此从频域模型所得的路面不平度的时域模型只能看成是满足给定路谱的全部可能的路面不平度中的一个样本函数。在

4、本文所讨论的方法中,谐波叠加法是目前使用较普遍的方法,特别适合用于国标道路谱时域模型的生成。尽管谐波叠加法计算量相对较大，但该方法理论基础严密、算法简单直观、数学基础严密、是一种高保真的频域模型转换方法。因此，本文采用谐波叠加法来构建随机路面不平度的频域模型。谐波叠加法拟合不平路面的原理是：设路面高程为平稳的、遍历的均值为0的Gaussian过程，则可以用不同形式的三角级数进行模拟。以下以正弦波为例加以描述。随机正弦波(或其它谐波)叠加法采用以离散谱逼近目标随机过程的模型，是一种离散化数值模拟路面的方法。随机信号，比如随机风压或路表轮廓，可以通过离散Fourier分析变换分解为一系列具有不同频

5、率和幅值的正弦波。谱密度就等于由带宽划分的这些正弦波幅值的平方。按照国标GB7031-86车辆振动输入一一路面平度表示方法路面不平度位移功率谱密度拟合表达式采用下式：GnGn)=Gn)dd0n031)式中：n式中：n参考频率，n0.1m，i；00G(n)路面不平度系数，m2/m-1；d0G(n)路面不平度，m2/m-1；w频率指数，经验值w=2。d已知在空间频率nnn内的路面位移谱密度为G(n)，利用平稳随机过12d程的平均功率的频谱展开性质，路面不平度的方差2为z2Jn2G(n)dn(3-2)znd1将区间(n，n)划分为n个小区间，取每个小区间的中心频率12n(i1,2,)处的谱密度值G(

6、n)代替G(n)在整个小区间内的值，则mid，idmid，id式(32)离散化后近似为,2为G(n)Anzdmid-iii1对应每个小区间，现在要找到具有频率n(i1,2,mid-i33)n)且其标准差为33)n)且其标准差为J2G(n)J2G(n)Ansindmid-i将对应于各个小区间的正弦波函数叠加起来(2冗nx+mid-i34)0)i就得到频域路面随机位移输nmid-inmid-i35)式中：自此，q(x)KJ2G(n)Ansin(2冗dmidiii100，2n上均匀分布的随机数；X频域路面的X方向通过离散频域路面的X方向的x值，就可以得到空间频域下随机路面Y方向的值。理论是这样的，然

7、后按照李军的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坐标仅仅表示路面的宽度；由这些节点按一

8、定的规律组成路面单元(Element)；再在路面单元里设置静摩擦系数和动摩擦系数，就能模拟真实的路面。39628RoadmapGenerationHowtoforAdams12,路谱的位置要处于轮胎的下方；第二，路谱向上的方向要指向轮胎所处的一侧；第三，路谱的大小要根；模拟真实的路面。39628RoadmapGenerationHowtoforAdams12,路谱的位置要处于轮胎的下方；第二，路谱向上的方向要指向轮胎所处的一侧；第三，路谱的大小要根；图3-1随机不平路面原理图龙忻r.；：:i町曲闷鼻吗CJLk：护严咿.也1而且，随机路谱的编制要满足轮胎的要求：第/Road/WidthNode据

9、仿真的需要确定。按照上述路面原理图和上一节的不平路面拟合理论，编制出随机路面VB程序（附录C）,程序的流程图如图3-2所示，程序界面如图3-3所示。图3-2随机路面文件生成的流程框图图3-3随机不平路面的VB程序界面根据上面的理论公式，选用国标GB7031-86车辆振动输入一一路面平度表示方法的B级路面编制成频域.rdf格式的路谱文件（附录D）。B级随机路面的频域模型如图3-4所示。路面X方向路面X方向图3-4B级随机路面的空间谱3.3随机路面的验证为了和国标GB7031-86车辆振动输入一一路面平度表示方法规定的路面空间频域谱图作比较，把程序生成的B级随机路面的数据点再作谱分析，得到按照谐波

10、叠加法理论拟合出B级随机路面的空间频域谱图如图3-5所示。0.0080.0060.0040.002-o-r00.010.020.030.040.050.060.070.080.091/m0.0080.0060.0040.002-o-r00.010.020.030.040.050.060.070.080.091/m0.1路面空间频域图3-5B级随机路面的据点的空间频域谱图根据国标GB7031-86车辆振动输入一一路面平度表示方法提供B级路面的“各倍频中心频率n(m-i)”值和“各倍频中心频率的位移功率谱密度g(n)cdc(10-6m2/m-i)”值拟合出标准B级路面的功率谱谱图。如图3-6所示，

11、把标准B级路面的功率谱谱图和编制生成的B级路面的功率谱谱图加以对比，可以很明显的看出，按照上面“谐波叠加法来构建随机路面不平度的频域模型”理论建立的随机路面是很好的。U-LlEEtlEEII111IIII|g0.008-0.007-0.006-0.005-0.004-Li.nns-0.002-Li.nni-图3-6U-LlEEtlEEII111IIII|g0.008-0.007-0.006-0.005-0.004-Li.nns-0.002-Li.nni-图3-6标准B级路面的功率谱谱图（黑点）和编制生成的B级路面的功率谱谱图（红色曲0.009-线）加以对比3.4小结本章依照谐波叠加法理论，并按照ADAMS软件的编制不平路面文件的格式要求，编制出B级随机路面。由