随机振动分析基础课件

《随机振动分析基础课件》由会员分享，可在线阅读，更多相关《随机振动分析基础课件（96页珍藏版）》请在金锄头文库上搜索。

1、随机振动分析基础随机振动分析基础1 1 随机振动的特点随机振动的特点对某振动运动，其规律显示出相当的随机性对某振动运动，其规律显示出相当的随机性而不能用确定性的函数来表达，使得只能用而不能用确定性的函数来表达，使得只能用概率和统计的方法来描述，这种振动被称为概率和统计的方法来描述，这种振动被称为随机振动。随机振动。随机振动可以由系统构成参数本身有随机性随机振动可以由系统构成参数本身有随机性而导致，但在多数情况下主要由激振源的随而导致，但在多数情况下主要由激振源的随机性所引起。本节主要研究这后一种情况，机性所引起。本节主要研究这后一种情况，即确定性系统在随机激励下的振动响应。即确定性系统在随机激

2、励下的振动响应。 汽车方面的典型例子是路面的随机凹凸汽车方面的典型例子是路面的随机凹凸不平使行驶的汽车产生随机振动；不平使行驶的汽车产生随机振动；被切削工件表层软硬不均使车刀及刀架被切削工件表层软硬不均使车刀及刀架产生随机振动；产生随机振动；风对建筑结构的随机激励；风对建筑结构的随机激励；地震对结构的随机激励；地震对结构的随机激励；浪使船舶产生随机振动；浪使船舶产生随机振动；大气湍流使机翼产生随机振动等等。大气湍流使机翼产生随机振动等等。 下图为一随机振动的时间历程样本函数表下图为一随机振动的时间历程样本函数表示。所谓样本函数是指随机振动本身是以示。所谓样本函数是指随机振动本身是以时间时间t为

3、过程参变量的函数过程。为过程参变量的函数过程。随机振动时间历程样本函数随机振动时间历程样本函数 从随机性的物理性质出发，这样各不相从随机性的物理性质出发，这样各不相同的函数应有无穷多个，每一个只是一同的函数应有无穷多个，每一个只是一个样本，最后构成集合个样本，最后构成集合xi(t), i=1, 2, 3, ; t 0, )。取尽各种可能性的无穷多个样本函数的取尽各种可能性的无穷多个样本函数的集合称为样本函数空间。集合称为样本函数空间。取取t=tk时刻各样本函数瞬时值构成一个序列时刻各样本函数瞬时值构成一个序列X(s, tk)=xi(tk), i=1, 2, 3, ; tk0, )；s用于表记对

4、应不同的样本函数。每个用于表记对应不同的样本函数。每个xi(tk)是是tk时刻的瞬时振动幅值，称为是随时刻的瞬时振动幅值，称为是随机变量机变量X(s, tk)的在的在s=si的一个样本点；的一个样本点；所有样本点的集合所有样本点的集合Ssi就是就是随机变量随机变量X(s, tk)的样本空间；的样本空间； 随机变量随机变量X(s, tk)随样本点的不同随机地随样本点的不同随机地取不同的值，即取不同的值，即X(s, tk)是样本点是样本点sS的的函数。同时注意它也是过程参变量函数。同时注意它也是过程参变量tk0, )的函数。的函数。当随机变量蕴含的是样本点的函数的意当随机变量蕴含的是样本点的函数的

5、意义明显且希望强调它是过程参变量义明显且希望强调它是过程参变量t的函的函数时，简记此随机变量为数时，简记此随机变量为X(t)。 当当tk取不同值时，可得不同时刻的随机变取不同值时，可得不同时刻的随机变量量X(tk)；从原理上看，对于各样本函数是；从原理上看，对于各样本函数是时间的连续函数的随机振动，只有时间的连续函数的随机振动，只有t连续变连续变化为无穷多个时刻而得出无穷多组随机变化为无穷多个时刻而得出无穷多组随机变量量X(t)才能完整地描述一个随机振动。才能完整地描述一个随机振动。这样实际形成的是以时间为过程参数的一这样实际形成的是以时间为过程参数的一族随机变量，这样定义的随机变量族就被族随

6、机变量，这样定义的随机变量族就被称为随机过程。随机振动是一种典型的随称为随机过程。随机振动是一种典型的随机过程。另外，也可以选用其它参数为随机过程。另外，也可以选用其它参数为随机过程的过程参数。机过程的过程参数。 2 2 相关函数和功率谱密度函数相关函数和功率谱密度函数 1. 相关函数相关函数掌握随机变量的性质是通过了解它的概率结掌握随机变量的性质是通过了解它的概率结构，构，最自然最自然是通过其概率密度函数是通过其概率密度函数p(x) 或概或概率分布函数率分布函数P(x)。完整地掌握。完整地掌握p(x)或或P(x)通常通常比较困难，因此比较困难，因此常用的统计描述常用的统计描述是讨论随机是讨论

7、随机变量的各低阶矩数字特征，如数学期望（均变量的各低阶矩数字特征，如数学期望（均值），均方值和方差等。值），均方值和方差等。作为增加了过程参数作为增加了过程参数t的随机变量族的随机过的随机变量族的随机过程，可通过对随机变量数字特征程，可通过对随机变量数字特征( (矩函数矩函数) )对对过程参数的扩展定义来研究其统计特性。过程参数的扩展定义来研究其统计特性。对图对图1.5-11.5-1所示随机振动，取离散时刻所示随机振动，取离散时刻t1, ,t2, tn可得一族随机变量可得一族随机变量X1, ,X2, , Xn。这些随机变量的概率结构可由概率密度函数这些随机变量的概率结构可由概率密度函数及不同时

8、刻的随机变量间的联合概率密度函及不同时刻的随机变量间的联合概率密度函数表达为数表达为 p(x1, t1)，p(x2, t2)，p(x1, t1; x2, t2)，p(x2, t2; x3, t3)，(1.5-1)1.5-1)， p(x1, t1; x2, t2; x3, t3; ; xn, tn) 上上述述表表达达的的n维维概概率率密密度度函函数数能能够够近近似似描描述述原原连连续续的的随随机机过过程程的的统统计计特特性性，n越越大大近近似似程程度度越越高高，当当n趋趋于于无无穷穷大大时时，(1.5-1)就就完完全表达了该随机过程的统计特性。全表达了该随机过程的统计特性。类似于研究随机变量统计

9、性质时对类似于研究随机变量统计性质时对各阶矩各阶矩的定义，可定义随机过程的各阶矩函数如的定义，可定义随机过程的各阶矩函数如下：下： k=1, 2, 3, , nk, j =1, 2, 3, , n k, i, j =1, 2, 3, , n (1.5-2) 可以证明，可以证明，用矩函数或用概率密度函数用矩函数或用概率密度函数（或概率分布函数）来描述随机过程数学（或概率分布函数）来描述随机过程数学上是等价的上是等价的。理论上完整地确定一个随机过程，需要确理论上完整地确定一个随机过程，需要确定定所有各阶矩函数所有各阶矩函数，明显这对实际应用来，明显这对实际应用来说又是一个过分的苛求。因此，实践上特

10、说又是一个过分的苛求。因此，实践上特别强调运用低阶矩即别强调运用低阶矩即1 1，2 2阶矩函数。阶矩函数。 1阶矩函数称为均值函数，定义为阶矩函数称为均值函数，定义为 (1.5-3)2阶矩函数称为相关函数，定义为阶矩函数称为相关函数，定义为 (1.5-4)(1.5-4)针对的是一个随机过程，因而可更针对的是一个随机过程，因而可更细分地称为自相关函数，以双下标细分地称为自相关函数，以双下标xx代表；代表；如果研究的对象包括有两个随机过程如果研究的对象包括有两个随机过程X(t), Y(t)，可以类似地定义出互相关函数如，可以类似地定义出互相关函数如 (1.5-5)均值函数和相关函数虽然只是随机过程

11、的矩均值函数和相关函数虽然只是随机过程的矩函数表达系列中的两个低阶矩函数，但它们函数表达系列中的两个低阶矩函数，但它们却表征了随机过程许多重要统计特征。却表征了随机过程许多重要统计特征。特别对一类实际上很常见的高斯随机过程，特别对一类实际上很常见的高斯随机过程，其高阶矩函数可以由其高阶矩函数可以由1 1，2 2阶矩函数表示，因阶矩函数表示，因此，对高斯随机过程，均值函数和相关函数此，对高斯随机过程，均值函数和相关函数完全表征了它的概率结构。完全表征了它的概率结构。而对于非高斯过程，这两矩函数也代表了其而对于非高斯过程，这两矩函数也代表了其统计性质中非常重要的一大部分。统计性质中非常重要的一大部

12、分。 高斯随机过程，又称为正态随机过程，是高斯随机过程，又称为正态随机过程，是这样一种随机过程：它在任意时刻这样一种随机过程：它在任意时刻tk的状的状态都服从正态分布，即是高斯随机变量。态都服从正态分布，即是高斯随机变量。定义为：对于任意定义为：对于任意n，X(t)的的n个样本为个样本为X(t1)，X(t2)，X(tn)，记，记x=x1，x2，xnT，m mX=m mX(t1)，m m X(t2)，m m X(tn)T，X(t)的的n维联合概率密度函数为维联合概率密度函数为 其中其中C为协方差矩阵为协方差矩阵记记CX(ti, tj)=E (X(ti)m mX(ti)( X(tj)m mX(tj

13、)，i，j=1, 2, ，n E. .为数学期望运算，有为数学期望运算，有 高斯随机过程是最常见的随机过程之一。高斯随机过程是最常见的随机过程之一。因为根据因为根据中心极限定理中心极限定理，当某随机变量是有，当某随机变量是有大量相互独立的随机因素的综合影响而形成，大量相互独立的随机因素的综合影响而形成，且每一个别因素所起作用都相对微小时，该且每一个别因素所起作用都相对微小时，该随机变量往往近似服从正态分布。随机变量往往近似服从正态分布。实际中满足中心极限定理条件的随机现实际中满足中心极限定理条件的随机现象很多，当现象也是过程参数象很多，当现象也是过程参数t的函数时的函数时即形成高斯随机过程。即

14、形成高斯随机过程。高斯随机过程的高阶统计函数可以由高斯随机过程的高阶统计函数可以由1 1，2 2阶函数表示的性质以及线性变换不改变阶函数表示的性质以及线性变换不改变其高斯分布特性的性质使它在理论和实其高斯分布特性的性质使它在理论和实践上都有很大重要性。践上都有很大重要性。 在常见的均值函数为常数（或简单函数）在常见的均值函数为常数（或简单函数）的情况下，习惯通过线性位移变换将考的情况下，习惯通过线性位移变换将考虑的运动静态工作点移到零，这导致均虑的运动静态工作点移到零，这导致均值函数为零，称为零均值化。值函数为零，称为零均值化。因此最关心的系统动态特性的统计特征因此最关心的系统动态特性的统计特

15、征就主要由就主要由2 2阶矩函数阶矩函数相关函数代表，相关函数代表，这更加突出了相关函数在随机振动分析这更加突出了相关函数在随机振动分析中的作用。中的作用。 2. 2. 平稳随机过程平稳随机过程如果一随机过程与过程参数如果一随机过程与过程参数t的起点无关，的起点无关，该随机过程即称为平稳随机过程；该随机过程即称为平稳随机过程；如果随机过程与过程参数如果随机过程与过程参数t的起点有关，的起点有关，该随机过程则称为非平稳随机过程。该随机过程则称为非平稳随机过程。平稳随机过程的概率结构可表为平稳随机过程的概率结构可表为 p(x1, t1)p(x1, t1+t t)，p(x1, t1; x2, t2)

16、=p(x1, t1+t t ; x2, t2+t t)， ，p(x1, t1; x2, t2; ; xn, tn) =p(x1, t1+t t ; x2, t2+t t ; ; xn, tn+t t) (1.5-6)其中其中t为任意常数。为任意常数。如果如果(1.5-6)各式都成立，随机过程称为强各式都成立，随机过程称为强平稳；如果只有前平稳；如果只有前2式成立，则称为弱平稳。式成立，则称为弱平稳。一般情况下，判断是否强平稳非常困难，一般情况下，判断是否强平稳非常困难，所谓平稳指的是弱平稳，简称平稳。所谓平稳指的是弱平稳，简称平稳。对于高斯随机过程，因为高阶概率密度函对于高斯随机过程，因为高阶

17、概率密度函数可由数可由1，2阶概率密度函数表达，弱平稳阶概率密度函数表达，弱平稳即为强平稳。即为强平稳。 对平稳随机过程，对平稳随机过程，1 1阶概率密度函数与参阶概率密度函数与参数数t无关，可写为无关，可写为p(x)。由由(1.5-3)可知它可知它的均值函数与时间无关，记为的均值函数与时间无关，记为m mx。2 2阶概率密度函数只与时间差阶概率密度函数只与时间差t t=t2- - t1有关，有关，可写为可写为p(x1, x2, t t)。对于平稳随机过程的。对于平稳随机过程的2 2阶矩函数阶矩函数相关函数，有下列性质相关函数，有下列性质: : (1)Rxx(t1,t2)= Rxx(t t)，

18、 Rxy(t1,t2)=Rxy(t t)，t t=t2 - - t1。(2)Rxx(t t)=Rxx(- -t t)； 当当t t=0时，有时，有Rxx(0 0)=Ex(t)x(t)=s sxx2，其中，其中s sxx2为均方值且与时间无关。为均方值且与时间无关。(3)Rxy(t t)= Ryx(- -t t)。 （4）（5）（6）具体有具体有特别特别 （7）如果）如果Rxx(t t)是关于是关于t t的衰减函数，且的衰减函数，且均值函数均值函数m mx=0，有，有 如果如果m mx不等于零，则不等于零，则 （8 8）设设Z(t)=X(t)+Y(t)为为两两个个平平稳稳随随机机过过程程X(t)

19、和和Y(t)之和，有之和，有Rzz(t t)= Rxx(t t)+Ryy(t t)+Rxy(t t)+Ryx(t t)； 如果如果X(t)和和Y(t)不相关，就有不相关，就有Rzz(t t)= Rxx(t t)+Ryy(t t)。 3. 3. 各态历经性各态历经性 上述概率密度函数或矩函数都是在样本上述概率密度函数或矩函数都是在样本函数空间中定义的。即使计算低阶矩函函数空间中定义的。即使计算低阶矩函数也要求在样本函数空间中平均。数也要求在样本函数空间中平均。理论上，样本函数空间要求无穷多个样理论上，样本函数空间要求无穷多个样本函数才算完备。本函数才算完备。 实际上可以获取的试验样本函数的个数实

20、际上可以获取的试验样本函数的个数通常是非常有限的，这就为实际数字特通常是非常有限的，这就为实际数字特征量计算要求的样本平均带来了困难。征量计算要求的样本平均带来了困难。为解决这个问题，并通过对实际的随机为解决这个问题，并通过对实际的随机过程的大量观察分析，提出了对某些随过程的大量观察分析，提出了对某些随机过程有各态历经性的假设。机过程有各态历经性的假设。 所谓各态历经假设，是建筑在这样的物所谓各态历经假设，是建筑在这样的物理观察和总结上的：理观察和总结上的：如果某一随机过程的概率结构与时间如果某一随机过程的概率结构与时间t的的起点无关，对它的一个样本函数，如果起点无关，对它的一个样本函数，如果

21、演进时间无穷长，此样本函数可取得其演进时间无穷长，此样本函数可取得其样本函数空间所包含的所有概率可能性；样本函数空间所包含的所有概率可能性；即它可以经历系统所有的可能状态，这即它可以经历系统所有的可能状态，这随机过程就被称为有各态历经性。随机过程就被称为有各态历经性。 在此假设下，各数字特征计算所要求的在此假设下，各数字特征计算所要求的多个样本函数的样本平均就转化为只要多个样本函数的样本平均就转化为只要求一个样本函数的时域平均。求一个样本函数的时域平均。 建立数字特征量的时间平均表达。对考建立数字特征量的时间平均表达。对考察长度为察长度为T的一段样本函数的一段样本函数x( (t) )进行时间进

22、行时间平均，各阶矩函数的表达式为平均，各阶矩函数的表达式为 (1.5-7) (1.5-7)各式中的各式中的s标明它们都是随所用样本标明它们都是随所用样本函数不同而不同的随机变量，它们也与考函数不同而不同的随机变量，它们也与考察时间长度察时间长度T有关。有关。如果有如果有 (1.5-8)即即(1.5-7)中的各式当中的各式当T趋于无穷大时都趋于无穷大时都有极限存在，并且这些极限还与选用那有极限存在，并且这些极限还与选用那一个样本函数无关，则此随机过程一个样本函数无关，则此随机过程X( (t) )是是严格各态历经的；严格各态历经的；如果如果(1.5-8)中只有前两式成立，则称此中只有前两式成立，则

23、称此随机过程随机过程X( (t) )是是2 2阶各态历经。一般情况阶各态历经。一般情况下所讲的各态历经性，指的就是下所讲的各态历经性，指的就是2 2阶各态阶各态历经。历经。很明很明显，各，各态历经随机随机过程一定是平程一定是平稳随机随机过程；反之程；反之则不一定成立。不一定成立。 4. 4. 功率谱密度函数功率谱密度函数 设平稳随机过程设平稳随机过程X(t)的自相关函数为的自相关函数为Rxx(t t)，Rxx(t t)的傅立叶变换存在，记为的傅立叶变换存在，记为(1.5-9) Sxx(w w)称为称为X(t)的自功率谱密度函数，是的自功率谱密度函数，是w w的非负实偶函数，的非负实偶函数，w

24、w为圆频率。为圆频率。Sxx(w w)与与Rxx(t t)形成傅立叶变换对，即又有形成傅立叶变换对，即又有 要要(1.5-9)和和(1.5-10)的傅立叶变换对存的傅立叶变换对存在，对变换函数有一定数学性质的要求。在，对变换函数有一定数学性质的要求。根据傅立叶变换存在性的绝对可积条件，根据傅立叶变换存在性的绝对可积条件，可以证明，可以证明，(1.5-9)存在的要求是存在的要求是 (1.5-11)(1.5-10)(1.5-11)成立要求成立要求 ，这，这要求过程要求过程X(t)的均值的均值m mx=0，此对一般平稳，此对一般平稳随机过程成立（或可以通过零均值化成立）随机过程成立（或可以通过零均值

25、化成立）。其物理基础是：随着两时刻间隔其物理基础是：随着两时刻间隔t t趋于无趋于无穷大，两时刻对应的随机过程值之间的相穷大，两时刻对应的随机过程值之间的相关性应趋于零。关性应趋于零。 如均值不等于零，可通过线性变换如均值不等于零，可通过线性变换Y(t)=X(t)- -m mx来进行零均值化。来进行零均值化。可以对照的是，对一般平稳随机过程可以对照的是，对一般平稳随机过程X(t)，由于平稳性的要求，肯定没有由于平稳性的要求，肯定没有 ，因此不满足绝对可积条件使得，因此不满足绝对可积条件使得X(t)的傅立的傅立叶变换并不存在。叶变换并不存在。因此直接对因此直接对X(t)作傅立叶变换以讨论其频谱作

26、傅立叶变换以讨论其频谱是要非常小心的。是要非常小心的。 对于工程常用的以对于工程常用的以H Hz z计的频率变量计的频率变量f，因为，因为w w=2=2p pf，并考虑到，并考虑到e-jwtwt= = coswt - wt - jsinwt wt 以及以及自相关函数的偶函数性质自相关函数的偶函数性质, ,有有 (1.5-12)(1.5-13) 可以证明，有可以证明，有 Sxx(f)=2p pSxx(w w) (1.5-14)由于工程上负频率无意义，故分析中有由于工程上负频率无意义，故分析中有时要使用单边谱密度，其定义为时要使用单边谱密度，其定义为 (1.5-15)对自相关函数有对自相关函数有

27、(1.5-16)自功率谱密度函数的物理意义。对自功率谱密度函数的物理意义。对(1.5-13)令令t t = 0，有，有 可见，过程可见，过程X(t)的平均能量的平均能量s sxx2可由可由Sxx(f)在全频带在全频带(， ) (或单边谱或单边谱Gxx(f)在正频带在正频带(0， ) )上积分得到，由上积分得到，由此说明此说明Sxx(f)代表了随机振动过程代表了随机振动过程X(t)的的平均能量在频率域上的分布情况。平均能量在频率域上的分布情况。 考察平稳随机振动位移考察平稳随机振动位移X(t)功率谱密度函数功率谱密度函数与对速度和加速度的功率谱密度函数的关与对速度和加速度的功率谱密度函数的关系。

28、系。设对设对X(t)功率谱密度函数为功率谱密度函数为Sxx(w w)，则由自，则由自相关函数的微分性质，有相关函数的微分性质，有 (1.5-17)(1.5-18)对应于考虑两平稳随机过程对应于考虑两平稳随机过程X(t)和和Y(t)的的互相关函数，可以定义互功率谱密度函互相关函数，可以定义互功率谱密度函数为数为 (1.5-19)其逆傅立叶变换为其逆傅立叶变换为(1.5-20)(1.5-21)即即Syx(w w)是是Sxy(w w)的复共的复共轭。上。上标“*”代表代表复共复共轭。 可以证明可以证明 (1.5-22) 5. 5. 平稳随机过程的谱特性分类平稳随机过程的谱特性分类由于自功率谱密度函数

29、由于自功率谱密度函数Sxx(f)（或（或Gxx(f)）反映了平稳随机过程的平均能）反映了平稳随机过程的平均能量随频率分布的特性，可以根据它分类量随频率分布的特性，可以根据它分类一些典型平稳随机振动过程，以便把握一些典型平稳随机振动过程，以便把握它们的本质特征。它们的本质特征。1) 窄带平稳过程窄带平稳过程典型谱密度函数，相关函数和时域样本典型谱密度函数，相关函数和时域样本函数分别如下图函数分别如下图 (a)，(b)和（和（c）。谱带）。谱带宽相比于它的中心频率宽相比于它的中心频率w wc大概要小一个大概要小一个数量级。谱分量主要集中在中心频率附数量级。谱分量主要集中在中心频率附近，说明振动能量

30、主要在近，说明振动能量主要在w wc附近。附近。 图图 (a)窄带过程单边谱窄带过程单边谱 (b)相关函数相关函数Rx(t t) (c)时域样本函数时域样本函数x(t)一般随机振动经过缓冲系统后，响应往一般随机振动经过缓冲系统后，响应往往是窄带随机振动。往是窄带随机振动。例如汽车在凹凸不平的道路上行驶时，例如汽车在凹凸不平的道路上行驶时，车身的振动即为窄带随机振动。这是因车身的振动即为窄带随机振动。这是因为汽车具有缓冲系统，只有车身部件的为汽车具有缓冲系统，只有车身部件的固有频率附近频带的振动才能传至车身，固有频率附近频带的振动才能传至车身，使车身在比较窄的频带范围振动。使车身在比较窄的频带范

31、围振动。 观察窄带随机振动的时域历程曲线可以观察窄带随机振动的时域历程曲线可以发现，它的峰值变化是随机的，但却近发现，它的峰值变化是随机的，但却近似地具有周期性，好像是峰值随时间随似地具有周期性，好像是峰值随时间随机变化的正弦振动，所以有时也称它为机变化的正弦振动，所以有时也称它为准正弦振动。准正弦振动。 2) 2) 宽带平稳过程宽带平稳过程 宽带平稳过程的谱带宽分布于较宽的范宽带平稳过程的谱带宽分布于较宽的范围，带宽与中心频率相比是同数量级的，围，带宽与中心频率相比是同数量级的，或更大。其典型谱密度函数，相关函数或更大。其典型谱密度函数，相关函数和时域样本函数分别如下图。和时域样本函数分别如

32、下图。 图图 (a)宽宽带过程单边谱带过程单边谱 (b)相关函数相关函数Rx(t t) (c)时域样本函数时域样本函数x(t)3) 理想白噪声过程理想白噪声过程白噪声过程是宽带平稳过程的极限特例。白噪声过程是宽带平稳过程的极限特例。对此过程有对此过程有Sxx(f)=So(或或Gxx(f)=2So)，So为常数。谱密度在整个频率轴为常数。谱密度在整个频率轴(,)上上均匀分布，即频带宽达整个频率轴。均匀分布，即频带宽达整个频率轴。白噪声过程的自相关函数为白噪声过程的自相关函数为Rxx(t t)=2p pS0d d(t)，是集中在，是集中在t t=0处强度为处强度为2p pSo的脉冲函数。的脉冲函数

33、。 而而t t=0处等于无穷大，处等于无穷大，t t0 0处都等于零处都等于零的的d d函数性质说明白噪声过程函数性质说明白噪声过程“自己自己”与与“自己自己”的相关性为无穷大，任的相关性为无穷大，任何不是何不是“自己自己”的两点间的相关性为的两点间的相关性为零。零。由于能量被要求均匀分布于整个频率轴，由于能量被要求均匀分布于整个频率轴，白噪声过程的均方值为白噪声过程的均方值为s sxx2 2= Ex2(t)= Rxx(0)=。这要求其平均能量为无穷大，。这要求其平均能量为无穷大，这在物理上是不可实现的。因此理想白这在物理上是不可实现的。因此理想白噪声只有理论上的意义。噪声只有理论上的意义。在

34、某一有限频带内有常数谱密度在某一有限频带内有常数谱密度So的有的有限带宽白噪声过程物理上是可能实现；限带宽白噪声过程物理上是可能实现；如果该有限频带带宽如果该有限频带带宽“相对相对”较大，可较大，可以一定程度上视为理想白噪声过程。以一定程度上视为理想白噪声过程。 3 3 线性系统在平稳随机激励下的响应线性系统在平稳随机激励下的响应 1. 1. 单自由度系统单自由度系统 设系统是确定性的，激励为随机性的，设系统是确定性的，激励为随机性的，对应的随机运动微分方程可写为对应的随机运动微分方程可写为 （1.5-231.5-23） 设设Fr(t)为为平平稳稳随随机机激激励励力力过过程程，X即即为为随随机

35、机位位移移响响应应过过程程，z z和和w wn如如前前定定义义为为确确定的系统特征参数。定的系统特征参数。 假定初始时刻系统处于静止，即假定初始时刻系统处于静止，即 （1.5-241.5-24） 线性系统随机响应分析研究的是线性系统随机响应分析研究的是随机微分随机微分方程的均方意义下的解过程方程的均方意义下的解过程，要获得该解，要获得该解过程的概率结构及建立它们与激励随机过过程的概率结构及建立它们与激励随机过程的概率结构之间的联系。程的概率结构之间的联系。根据线性随机常微分方程解的构成理论，根据线性随机常微分方程解的构成理论，在线性情况下，对确定性系统，随机激励在线性情况下，对确定性系统，随机

36、激励的的均方解均方解与确定性激励的解有相同的形式。与确定性激励的解有相同的形式。由此，根据前零初始条件下响应的杜哈梅由此，根据前零初始条件下响应的杜哈梅积分表达式（积分表达式（1.2-281.2-28），有），有 h(t)为（为（1.2-251.2-25）定义的系统的单位脉冲）定义的系统的单位脉冲响应函数。考察响应过程的数字特征，如响应函数。考察响应过程的数字特征，如前所述，仍主要关心其前所述，仍主要关心其1 1，2 2阶矩函数。阶矩函数。 响响应均均值函数函数为 （1.5-251.5-25） 将（将（1.2-251.2-25）代入上式，有）代入上式，有 考虑稳态响应令考虑稳态响应令t，有，有

37、 系统（系统（1.5-231.5-23）的频率响应函数）的频率响应函数H(w w)为为（1.5-261.5-261.5-261.5-26） （1.5-271.5-271.5-271.5-27） （1.5-281.5-281.5-281.5-28） 为激励的均值，如果它为零，则响为激励的均值，如果它为零，则响应均值亦为零；否则，稳态响应的均值应均值亦为零；否则，稳态响应的均值将以将以H(0)为乘数放大；对瞬态情况，为乘数放大；对瞬态情况，（1.5-261.5-26）表明响应均值函数和时间有）表明响应均值函数和时间有关，因此是非平稳的。关，因此是非平稳的。 响应自相关函数为响应自相关函数为响应自相

38、关函数为响应自相关函数为 由（由（由（由（1.5-291.5-291.5-291.5-29）可）可）可）可见见，即使激励是平，即使激励是平，即使激励是平，即使激励是平稳稳随机随机随机随机过过程，程，程，程，响响响响应应一般也不是平一般也不是平一般也不是平一般也不是平稳稳随机随机随机随机过过程。程。程。程。这对计这对计算响算响算响算响应应的数字特征的数字特征的数字特征的数字特征带带来麻来麻来麻来麻烦烦，所幸的是可以通，所幸的是可以通，所幸的是可以通，所幸的是可以通过过考考考考虑虑稳态稳态情况使情况使情况使情况使问题问题相相相相对简对简化。化。化。化。 （1.5-291.5-291.5-291.5

39、-29） 由定义可知由定义可知将上代入相关函数（将上代入相关函数（1.5-291.5-29），有），有 （1.5-301.5-301.5-301.5-30） （1.5-311.5-311.5-311.5-31） 将（将（1.2-251.2-25）的）的h( (t) )代入（代入（1.5-311.5-31），所），所得再代入（得再代入（1.5-301.5-30），整理后考虑稳态），整理后考虑稳态情况令情况令 令令t t= =t2- -t1，对照定义，对照定义 （1.5-321.5-321.5-321.5-32） （1.5-331.5-331.5-331.5-33） （1.5-331.5-33）是

40、平稳随机激励与响应的功率）是平稳随机激励与响应的功率谱密度函数之间的关系，注意它只适用谱密度函数之间的关系，注意它只适用于稳态情况。于稳态情况。 以以类似的方法，可以推似的方法，可以推导出响出响应和激励和激励的自的自谱及互及互谱密度函数之密度函数之间的关系，有的关系，有 （1.5-341.5-341.5-341.5-34） （1.5-351.5-351.5-351.5-35） 2. 2. 多自由度系统多自由度系统 对于线性振动多自由度系统，随机运动微分方对于线性振动多自由度系统，随机运动微分方对于线性振动多自由度系统，随机运动微分方对于线性振动多自由度系统，随机运动微分方程扩展为矩阵方程，即程

41、扩展为矩阵方程，即程扩展为矩阵方程，即程扩展为矩阵方程，即 由于考虑的系统本身是确定性的，质量、阻尼由于考虑的系统本身是确定性的，质量、阻尼由于考虑的系统本身是确定性的，质量、阻尼由于考虑的系统本身是确定性的，质量、阻尼及刚度矩阵与前（及刚度矩阵与前（及刚度矩阵与前（及刚度矩阵与前（1.3-221.3-221.3-221.3-22）中定义的一样。）中定义的一样。）中定义的一样。）中定义的一样。（1.5-361.5-361.5-361.5-36） 和单自由度系统解法类似，随机响应和单自由度系统解法类似，随机响应向量的解法可以利用对确定性的线性向量的解法可以利用对确定性的线性多自由度系统的模态分析

42、法。多自由度系统的模态分析法。 设系系统本身本身满足利用足利用实模模态分析法的分析法的条件，令条件，令 （1.5-371.5-371.5-371.5-37） 为系统的正则化实模态矩阵；为系统的正则化实模态矩阵；Y为对应为对应的随机模态坐标响应向量。将（的随机模态坐标响应向量。将（1.5-1.5-3737）代入（）代入（1.5-361.5-36），左乘），左乘 ，由正，由正交性得解耦的随机模态响应坐标方程为交性得解耦的随机模态响应坐标方程为 （1.5-381.5-381.5-381.5-38） 假设系统初始条件为零，则各随机模态假设系统初始条件为零，则各随机模态坐坐 标响应的初始条件都为零。标响

43、应的初始条件都为零。 对杜哈梅杜哈梅积分，分， 可以可以证明它的一个等价形式明它的一个等价形式为 （1.5-391.5-391.5-391.5-39） （1.5-251.5-251.5-251.5-25） 在随机响应分析时，认为（在随机响应分析时，认为（1.5-391.5-39）与）与（1.5-251.5-25）是有差别的。差别在于）是有差别的。差别在于（1.5-391.5-39）的积分上下限扩展到正负无）的积分上下限扩展到正负无穷大是对应于系统的稳态响应；而穷大是对应于系统的稳态响应；而（1.5-251.5-25）积分上下限在）积分上下限在 0, ,t 之间而对之间而对应瞬态响应情况。用（应

44、瞬态响应情况。用（1.5-391.5-39）求取随）求取随机模态坐标响应机模态坐标响应Yi(t)的稳态解，为的稳态解，为 利利用用（1.5-401.5-40）可可求求取取各各Yi(t)的的均均值值函函数数，进进而而通通过过（1.5-371.5-37）求求得得各各Xi(t)的的均均值值函函数。数。（1.5-401.5-401.5-401.5-40） 计算对计算对Y=Y1(t) Y2(t) Yn(t)T 的相的相关函数矩阵。由定义有关函数矩阵。由定义有（1.5-411.5-411.5-411.5-41） 为为nn维对角矩阵维对角矩阵 对（对（1.5-411.5-41）两边进行傅立叶变换，有）两边进

45、行傅立叶变换，有 其中其中 HF F (w w)= F (hF F(t)。对对因为因为 Q(t)=F FT Fr(t)所以所以 （1.5-421.5-421.5-421.5-42） 其中其中（1.5-431.5-431.5-431.5-43） （1.5-441.5-441.5-441.5-44） 由类似的推导，对互谱矩阵，有由类似的推导，对互谱矩阵，有 （1.5-451.5-451.5-451.5-45） 可以利用复模态方法推导有非比例阻尼可以利用复模态方法推导有非比例阻尼的线性多自由度系统的随机激励功率谱的线性多自由度系统的随机激励功率谱与随机稳态响应功率谱之间的关系，它与随机稳态响应功率谱

46、之间的关系，它们也有（们也有（1.5-431.5-43）和（）和（1.5-451.5-45）形式。）形式。但如前所述，实模态分析法在大多数情但如前所述，实模态分析法在大多数情况下，已满足实际应用的需求。况下，已满足实际应用的需求。 例例 考虑如图的轿车模型，并增加考虑如图的轿车模型，并增加C1，C2以以代表轮胎和悬架的阻尼效应，设前后轮代表轮胎和悬架的阻尼效应，设前后轮与地面接触处的不平度激励与地面接触处的不平度激励X1(t)，X2(t)为零均值平稳随机过程，求随机响为零均值平稳随机过程，求随机响应应X，Q的稳态统计特征。的稳态统计特征。 轿车模型轿车模型 解：增加了阻尼效应的对随机基础激励的

47、解：增加了阻尼效应的对随机基础激励的运动微分方程为运动微分方程为 其中其中 对上两式两边作傅立叶变换，可整理出对上两式两边作傅立叶变换，可整理出频率响应函数矩阵频率响应函数矩阵H(w w)为为 各量定义为各量定义为 记记Z(w w)=X(w w) Q Q(w w)T， P(w w)=X1(w w) X2 2(w w)T，由上有由上有Z(w w)=H(w w)P(w w) 参照（参照（1.5-431.5-43），有），有 （a a a a） 设路面不平度激励设路面不平度激励X1(t)，X2(t)为零均为零均值平稳高斯随机过程，所以值平稳高斯随机过程，所以X1(t)与与X2(t)之间仅相差一时差之

48、间仅相差一时差t t0=l/v，v为车速，为车速，l为轮距，则路面激励的功率谱矩阵为为轮距，则路面激励的功率谱矩阵为 （b b b b） Szz(w w)的对角线元素分别为随机响应的对角线元素分别为随机响应X，Q Q的自功率谱密度函数。由上的自功率谱密度函数。由上（a a a a）式有式有 响应的互谱密度为响应的互谱密度为 考虑考虑(b)(b)式，整理后有式，整理后有 上两式说明，当两个激励只相差一个时差上两式说明，当两个激励只相差一个时差t t0时，系统两个自由度的响应相当于每个时，系统两个自由度的响应相当于每个自由度在一个激励下的响应；其频率响自由度在一个激励下的响应；其频率响应函数，分别

49、为应函数，分别为 和和 。 这是在考虑简化车辆动力学模型时，近这是在考虑简化车辆动力学模型时，近似地将车身的垂直自由度运动和俯仰自似地将车身的垂直自由度运动和俯仰自由度运动分别考虑的基础。由度运动分别考虑的基础。 本章结束本章结束习题题1 随机过程Xt的样本函数为： 式中a1，a2，w1，w2是常数，f1，f2为统计独立的在0，2p上均匀分布的随机变量，求自相关函数Rxx(t)。题2 某平稳随机过程的自相关函数为： 求其均值mx，方差，功率谱密度函数Sxx(f)和单边谱密度函数Gxx(f)。 题3 已知某振动系统的输入为力，输出为位移，系统位移响应的y(t)的自功率谱为： 求响应y(t)的自相关函数和均方值。 题4 系统示意如下图，设F1 (t)为均值为零的白噪声，其自功率谱密度函数为SFF(w)，求稳态情况下响应的自功率谱密度函数，互功率谱密度函数及各响应的均方值。 题5 如下图，系统由主系统（m1，k1）和副系统（m2，C2，k2）组成，设作用在m1上的F1(t)为零均值白噪声，试以响应y1(t)的均方值最小为条件确定副系统的m2，C2，k2。 题6 设线性系统随机运动方程为 其中： W(t)为 平 稳 白 噪 声 激 励 向 量 ， 有EW(t)=0，EW(t)WT(t+t)=Id(t)，I为单位矩阵，用实模态分析法求响应的相关函数矩阵RXX(t)。