第七章随机振动的响应分析

1、第七章随机振动的响应分析第七章 随机振动的响应分析7-1 单输入单输出的线性系统 7-2 多输入多输出的线性系统本章讨论机械或结构系统在随机激励作用下，激励激励系统系统响应响应三者之间的关系。系统有线性与非线性之分。大量工程问题，线性模型可得到逼真的结果。本课程只讨论线性系统问题。随机激励分两类：参数激励与非参数激励参数激励：系统本身的某些参数(如质量、刚度、阻尼等)随时间随机地变化而引起振动。非参数激励即由外界施加的激励。非参数激励又分为平稳的和非平稳的两类。本章研究常参数线性系统对平稳随机激励的响应当系统的激励(输入)是平稳过程时，由于常参数的假设，系统的响应(输出)也一定是平稳的。对于一

2、个常参数线性系统，它往往可能在不同位置上同时受到激励，即有多个输入；其响应也可能有很多个，而且不同位置处的响应也不同。对于线性系统来说，多输入与多输出问题可以在单输入与单输出问题的基础上应用叠加原理得到解决对于确定性振动，激励与响应之间的关系，一般用微分方程来描述，方程的非齐次项是确定的，初始条件也是确定的，因此响应也是确定的。在随机振动中，一般激励与响应都必须用概率统计的方法来描述。在激励与系统特性已知的情况下，只能求出响应的一些统计特征，如期望（均值）、相关函数、功率谱密度、均方值等。7-1 7-1 单输入单输出的线性系统单输入单输出的线性系统假定常参数线性系统只受到一个输入x(t)的作用

3、，其相应的响应(输出)为y(t)，如图所示。常参数线性振动系统常参数线性振动系统y(t)Output (response)输出（响应）输出（响应）x(t)Input (excitation)输入（激励）输入（激励）本章研究输入、输出和系统动态特性三者之间的关系，以及计算响应（输出）的统计特征的方法设x(t)是平稳随机过程X(t)的一个样本函数则系统输出y(t) 是另一平稳随机过程Y(t)的一个样本函数设系统的脉冲响应函数h(t)，则频率响应函数是H()。常参数线性振动系统常参数线性振动系统y(t)Output (response)输出（响应）输出（响应）x(t)Input (excitatio

4、n)输入（激励）输入（激励）一、响应的均值对于输入的一个样本函数，由卷积积分公式，可得输出的一个样本函数( )() ( )dy tx th( )() ( )dy tx th设想对于输入中的每个样本函数，都可按上式写出其对应的输出的样本函数。于是，可得到输出的集合平均为: ( )() ( )dE Y tEX thdhtXEtYE)()()(xtXEtXE)()(dhtYEx)()(dehHj)()(0dhH)()0( 输入与输出均值的关系式为： ( )= (0)YXE Y tHdhtYEx)()(dhH)()0(H(0)是一个常数，它表示输入X(t)与输出Y(t)中，频率=0这一成分（即直流分量

5、）之间的传递关系。)()()0(txtyH直流分量上式表明，当输入是平稳过程时，输出的均值与输入的均值只差一个乘子H(0)。若输入的均值为零，则输出的均值也一定为零。此结论可以推广到多输入与多输出的情形。 ( )=(0)YXE Y tH二、响应的自相关函数输出过程Y(t)的自相关函数定义为: ( ) ()E Y t Y t上式为输出的自相关函数之间的关系式。该式说明，对于常参数线性系统，若激励是平稳随机过程，则响应的自相关函数与自然时间无关，也一定是平稳的随机过程。则响应的自相关函数可表示为:122112( ) ( ) ()=() ()()ddYXRE Y t Y thhR 对输出的自相关函数

6、作傅立叶变换，便得到响应的自功率谱密度SY()为三、响应的自功率谱密度函数三、响应的自功率谱密度函数变换积分次序，并重新排列jj122112( )( ) ( ) ()()d dd YYXSRedehhR 1221jj1122j2121ddY()XS ()h()edh()eR ()e() 令=-1+2，由维纳辛钦关系式知，最后一个积分就是激励X(t)的自谱密度： j( )( )dXXSRe第二个积分就是脉冲响应函数h(2)的傅立叶变换，即频率响应函数H()。2j22( )()dHhe1221jj1122j2121ddY()XS ()h()edh()eR ()e() 前两个积分的不同在于指数中的正

7、负号的差别。经处理后得随机输入与输出的自谱密度关系式:上式是随机振动理论中一个极其重要的公式，指出了输入、输出与系统动态特性三者之间的关系。 1221jj1122j2121ddY()XS ()h()edh()eR ()e() 111)()()(dehHHj2( )( )( )( )( )( )YXXSHHSHS上式表明，若已知系统的增益因子|H()|和输入的自谱密度SX()，则可确定输出的自谱密度SY()。事实上，若已知SX() 、|H()|和SY() 三者中的任意两个，就可以确定第三个。此外，响应的自谱密度是与系统的相位因子无关的。2( )( )( )( )( )( )YXXSHHSHS已知

8、响应的自谱密度SY()，则可计算出响应的均方值EY2: 21(0)( )d2YYE YRS四、响应的均方值四、响应的均方值将随机输入与输出的自谱密度关系式代入上式2221( )( )d2YXE YHS注意：当均值为零时，均方值就等于方差。22YY在输入为理想白噪声的情况下，由于输入的自谱密度对于所有的频率都是常数，则响应的均方值公式可得到简化:220( ) d2YSH只要计算出如下的广义积分I值，便可求得响应的均方值：2( ) dIH2221( )( )d2YXE YHS由互相关函数的定义，可得激励与响应之间的互相关函数:( )( ) ()( )( )()d ( ) ( )()d( )()dX

9、YXRE X t Y tE X thX thE X t X thR五、激励与响应的互相关函数五、激励与响应的互相关函数常参数线性系统在受到平稳随机输入时，激励与响应之间的互相关函数正好等于脉冲响应函数与输入自相关函数的卷积( )() ( )dy tx th( )( )()dXYXRhR对互相关函数表达式作傅立叶变换，便可得到激励与响应之间的互功率谱密度。六、激励与响应的互谱密度六、激励与响应的互谱密度jjj-( )( )d =( ) ()d ( )()d d XYXYXSRehRdehRe 与前面计算响应自谱相似的方法，将上式改写为：jj()( )( )d()d() XYXSheRe ( )(

10、 )()dXYXRhR上式第一个积分是频率响应函数H()，第二个积分就是激励X(t)的自谱密度SX()( )( )( ) XYXSHS上式表明：输入与输出之间的互谱密度等于系统的频率响应函数与输入自谱密度函数的乘积。通过该式可完整地确定系统的频率特性H()。jj()( )( )d()d() XYXSheRe 由于H()是复数，它可表示为:( )( )j ( )HAB则互谱密度可以表示为:( ) ( )j ( )( )XYXSABS由于SX()是实偶函数，则互谱密度函数可表示为：( )( )( )( )( )arctg( )XYXXYBSHSA上式表明：互谱密度的幅值等于系统的增益因子与 输入自

11、谱的乘积互谱密度的幅角又等于系统的相位因子。( )( )( )( )( )arctg( )XYXXYBSHSA七、相干函数七、相干函数22( )( )( )( )XYXYXYSSS系统输入与输出的谱相干函数（又称凝聚函数）可通过下式来定义：SX()和SY()皆为实函数，故相干函数必为实函数。可以证明，对于所有频率，相干函数满足以下不等式：20( )1XY当输入与输出互不相关时，有RXY()=0，从而互谱密度SXY()=0，于是由定义知相干函数也等于零。对于线性系统，存在下列关系2( )( )( )( )( )( )YXXYXSHSSHS在线性系统的假设下，输入输出线性相关，有22222( )(

12、 )( )( )1( )( )( )( )( )XYXXYXYXXSHSSSSHS输入输出互不相关时，相干函数的值等于0；输入输出线性相关时，相干函数等于1。相干函数的值在0与1之间。如果相干函数值大于零但小于1，为以下三种情况之一 联系输入X(t)和输出Y(t)的系统是非线性的(2) 测量中有外界噪声干扰(3) 输出Y(t)是输入X(t)和其它输入的综合输出。j0( )tx tx e如图所示的单输入线性系统，假定只在输出测量中混有噪声，则实测得到的输出Z(t)是真实输出Y(t)与噪声干扰N(t)之和。( )( )( )Z tY tN tN(t)Z(t)H()X(t)Y(t)只讨论一种存在噪声

13、干扰的情况：假定X(t)与N(t)皆是均值为零的平稳随机过程，且N(t)与X(t)和Y(t)都是不相关的，则有：0)()()()(YNNYXNNXRRRR)()()()()()()()()()(XYXNXYXZRRRtNtXtYtXEtZtXER输入与实测输出之间的互相关函数：故有:)()(XYXZSS实测输出的自相关函数)()()()()()()()()()()()()(NYNNYYNYZRRRRRRtNtYtNtYEtZtZER实测输出的自谱密度)()()(NYZSSS输入X(t)与实测输出Z(t)的谱相干函数:( )( )( )( )( )XZXYZYNSSSSS22( )( )( )(

14、 )XZXZXZSSS得到222Z2( )( )( )1( )( )( )( )XXXXNHSSHSS将互谱密度与自谱密度式代入上式，并结合下式2( )( )( )( )( )( )XYXYXSHSSHS上式表明：在有噪声干扰的情况下，输入与实测输出的谱相干函数将小于1。因此，对于线性系统，可借助相干函数值来判断干扰影响的大小。)()()()()()()(XXZXXYXZSSHSHSS此外：虽输出中含有干扰，但通过实测信号的互谱密度以及输入信号的自谱密度可以精确的获得系统的频响特征222Z2( )( )( )1( )( )( )( )XXXXNHSSHSS7.2 7.2 多输入多输出的线性系统

15、多输入多输出的线性系统考虑某一具有m个输入Xi(t) (i=1,2, m)和n个输出Yk(t) (k=1,2, n)的常参数系统，假定每个输入Xi(t)都是平稳的随机过程。( ), 1,2, kihtim在系统有m个输入Xi(t)的情况下，对应于每一个输出Yk(t)，有m个脉冲响应函数：对于n个输出，则共有nm个脉冲响应函数，脉冲响应以矩阵形式可表示为：111212122212( )( )( )( )( )( )( )( )( )( )mmnnnmhthththththth thththt矩阵中各元素均加以两个脚标第一个脚标k(k=1,2, n)表示k处的响应(输出)；第二个脚标i(i=1,2

16、, m)表示i处的激励 (输入)频率响应函数是脉冲响应函数的傅立叶变换，因此，图示系统共有nm个频率响应函数，频率响应矩阵H()可表示为111212122212( )( )( )( )( )( )( )( )( )( )mmnnnmHHHHHHHHHH已知系统的激励和动态特性，便可确定系统响应的各个统计特征。m个输入和n个输出 可以表示为1122( )( )( )( )( ) ( ) ( )( )mnX tY tXtY tX tY tXtY t一、响应的均值一、响应的均值对于线性系统，每一个输出Yk(t) (k=1,2, n)都可以由对应于每个独立输入的响应Yki(t) (i=1,2, m)叠

17、加而成，如图所示。121( )( )( )( )( )mkkkkmkiiY tYtYtYtYt假定各个输入Xi(t)都是平稳的随机过程，则有iXXEi常数对应于每个独立输入(第i个输入)的响应Yki(t)的期望为(0)( )dkiiikiYXkiXkiE YHh11(0) (1,2, )kimmYkkiXkiiiE YE YHknk处总响应的期望为写成矩阵形式：XYH)0(YnYYY21XmXXX21111212122212(0)(0)(0)(0)(0)(0)(0)(0)(0)(0)mmnnnmHHHHHHHHHH系统的每一个输出Yk(t)是对应于各个独立输入的响应的叠加：11111( )(

18、)( )()d (1,2, )mmkkikiiiiY tYthX tkn二、响应的相关矩阵系统的响应以矩阵形式的形式可以表示为下式111( )( )()dY thX t其中Y(t)与X(t)是给出的列阵，而h(1)表示脉冲响应矩阵。( ):1( ):( ):1Y tnh tnmX tm111( )( )()dY thX t12( )( )( )( )mX tXtX tXt111212122212( )( )( )( )( )( )( )( )( )( )mmnnnmhthththththth thththt12( )( )( ) ( )nY tY tY tY t11111( )( )( )()

19、d (1,2, )mmkkikiiiiY tYthX tkn222()()()dY thX t对上式进行转置得222()()()dTTTYtXthn个输出的自相关与互相关函数构成一个nn阶的输出相关矩阵RY()1 11 212 12 2212( )( )( )( )( )( )( )( )( )( )nnnnn nYYYYYYY YY YY YYY YY YY YRRRRRRRRRR():1Y tn():1TYtn)()()()(21tYtYtYtYn)()()()(21tYtYtYtYnT上式也可以下列形式表示写成111222( ) ( )() ( ) ()d()()dTYTTRE Y t

20、YtEhX tXth若将该式中两个积分的乘积改写成二重积分，并交换求平均与积分的次序，可得 ()(1)nmm(1)()mmn112221121221( )() ()()()dd ()()()ddTTYTXRhE X tXthhRh111222( ) ( )() ( ) ()d()()dTYTTRE Y t YtEhX tXth该表达式给出了多输入与多输出系统的输出相关矩阵与输入相关矩阵之间的关系式，其中，输入相关矩阵为:( )( )()TXRE X t Xt系统的n个输出的自谱与互谱构成了一个nn阶的输出功率矩阵SY()。1 11 212 12 2212( )( )( )( )( )( )(

21、)( )( )( )nnnnn nYYYYYYY YY YY YYY YY YY YSSSSSSSSSS三、响应的谱矩阵三、响应的谱矩阵其中矩阵元素为输出功率矩阵也可以表示成下列表达式的形式：j( )( )d ( ,1,2, )k lk lY YY YSRek lnj( )( )dYYSRe将输出相关矩阵表达式代入上式，可得多输入与多输出系统的输出功率谱矩阵与输入功率谱矩阵间关系式1122j112221jj()j11121222( )()()()dd d ()d()d()()d ( )( )( ) TYXTXTXShRheheReheHSH 112221( )()()()ddTYXRhRh()

22、()()n mm mm n其中H()是系统的频率响应矩阵，HT()是H()的转置矩阵，H()的共轭矩阵为1j11( )()( )dHHhe而输入功率谱矩阵为j( )( )dXXSRe1122j112221jj()j11121222( )()()()dd d ()d()d()()d ( )( )( ) TYXTXTXShRheheReheHSH )()()()(TXYHSHS系统的m个激励和n个响应构成一个mn阶互相关矩阵RXY() 1 11 212 12 2212( )( )( )( )( )( )( )( )( )( )nnmmm nX YX YX YX YX YX YXYX YX YX Y

23、RRRRRRRRRR四、激励与响应的互相关矩阵四、激励与响应的互相关矩阵矩阵的元素为1( )( )()( )( )()d ( 1,2,; 1,2, ) j kmX YjkjkiiiRE Xt Y tE XthX tjmkn矩阵的元素为激励与响应的互相关矩阵也可表示为： dhRdhtXtXEdhtXtXEtYtXERTXTTTTTXY)()()()()()()()()()()(系统的m个激励和n个响应构成一个mn阶的互谱矩阵SXY() 1 11 212 12 2212( )( )( )( )( )( )( )( )( )( )nnmmm nX YX YX YX YX YX YXYX YX YX YSSSSSSSSSS五、激励与响应的互谱函数五、激励与响应的互谱函数其中，矩阵元素的定义为( )( )d (1,2,; 1,2, )j ki kjX YX YSRejmkn激励与响应的互谱矩阵也可表示为jj()j( )()( )d d ()d()( )d( )( )TXYXTTXXSRheReheSH 类似地，有下列表达式成立：)()()(1TYYXHSS)()()()()()(1XTTXSHHHSH)()()(TXXYHSS)()()(*TXXYHSS)()()(*XTXYSHS)()(*TXYYXSSY为输入，X为输出Y X ，H