随着可靠性试验和环境试验的要求以及计算机技术的发展，振动试验从一开始的简单正弦定频试验到正弦扫频试验，然后再发展到随机振动试验，一步步地提高着振动试验的精度。
但是，最近的随机振动试验，尤其在运输振动试验中，不仅没能再现实际负载环境条件，而且即使试验合格还会出破损事故。也就是说，对试验精度的要求越来越严格。因此，急切需要对现在的随机振动试验进行改良，提供一种更接近实际负载情况和满足高精度试验要求的试验方法。于是，非高斯型随机振动试验方法便被逐渐地提上了应用的日程，随着技术的发展，现如今较好的振动控制仪基本上都能对应此种试验方法。
先来谈谈现在应用最广的随机振动试验（高斯型），是一种含有各种频率波形成份的不规则振动。通过设定试验条件的加速度密度谱（PSD），然后通过振动控制仪按照此PSD产生随机波形进行振动。如下图1便是一般的运输试验规格PSD曲线（JIS Z 0232）。
图1 加速度密度谱示例
随机振动试验，加速度的概率密度函数分布为正态分布（高斯分布，高斯型命名的来源），一般用下式表示，
式中，x为加速度，m为加速度的平均值，f（x）为概率密度分布，σ为标准偏差。图2为其概率密度分布曲线，图中可以看出随机振动试验的加速度符合高斯分布，超过3σ的大加速度出现的概率极低，仅有0.3%。
图2 高斯分布
接着再来看实际情况下发生的振动，以运输试验为例。在实际运输中，由于路面凹凸等的影响，加速度不满足高斯分布的情况很多。比如图3显示的是某货车架上的振动波形，图4是其振动加速度的概率分布情况。图4中可以看出实际分布不满足高斯分布，超过3σ的大加速度出现的概率有1.5%，相当于高斯分布的5倍。
图3 某货车架实际振动波形
图4 波形（图3）对应的概率密度曲线
第三，如何来定量非高斯型振动试验。还是以上面实际运输时的振动为例，发现货车架上实际的加速度分布不是高斯分布的场合比较多。通过学者的研究，提出了偏度S（Skewness）和峭度K（Kurtosis）两个概念来定量化非高斯型随机振动试验，定义式子如下，概率论学的好的，是不是感觉和概率论里面的3阶中心距和四阶中心矩的概念差不多？
式中，x为加速度，m为加速度平均值，N为数据点数，σ为标准偏差。
来看看偏度的含义。高斯型平稳随机振动的概率密度函数曲线是相对于均值m对称分布的，其偏度S为零。若平稳随机振动时间历程x（t）的偏度S不为零，则表明概率密度函数曲线关于m不对称。当S大于零的场合，概率密度函数曲线的峰值偏左，如图5（a），；当S小于零的场合，曲线的峰值偏右，如图5（b），其峰值对应于平均值m。
图5 偏度示意图
再来看看峭度的含义。峭度的取值范围为1～∞。高斯型平稳随机振动的峭度K等于3。与具有相同均值和均方值的高斯随机振动相比，当K小于3时，非高斯型的概率密度曲线的峰值较低，整个峰较宽，但尾部的量值较低；对于均匀分布信号，K为1.8。当K大于3时，非高斯型的概率密度曲线的峰值较高，整个峰较窄，但尾部的量值较高，极大值出现的概率较大，具体如图6所示。
图6 不同峭度的概率密度函数曲线（平均值为零，方差为1）
有了上面的概念，实际中的波形和概率密度曲线又是什么样的呢？看图7，高斯型随机平稳振动，即现今最普及的随机振动试验，K=3，S=0。
图8中，K值偏大（远大于3），就会出现较多的极大值。
图9中，S值偏大（远大于零），加速度会向正方向偏移分布较多。当然，如果S远小于零，加速度会向负方向偏移分布较多。
结束语：
通过以上的介绍，相信会对非高斯型随机振动试验有一定的认识。进入21世纪后，针对偏度和峭度的量化算法，通过计算机硬件产生非高斯型振动波形的研究如火如荼，并取得了很可喜的成果，已经开始在振动试验中得到了应用。有资源的工程师，建议自己编辑一下试验条件，运行一下非高斯型随机振动试验，尤其是通过耳朵可以比较一下两者之间的运行声音区别。
非高斯随机振动试验已在运输试验中得到广泛应用。目前，其在产品的疲劳评价试验中尚未推广开来，很多学者对于里面的评价方法等细节还在研究中，一旦得到广泛应用，相信对振动试验的精度产生飞跃性的影响，让我们拭目以待。
本文内容大部分来自日文原稿：
1 面向振动试验精度向上的技术课题—随机振动从高斯型向非高斯型，细山 亮、中岛 隆胜，2009年3月31日
参考资料：
1 《振动环境试验》，国防工业出版社，王帅、夏益霖、荣克林 著，P21非高斯平稳随机振动。
联系客服
