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

第七章随机振动的响应分析第七章随机振动的响应分析§7-1单输入单输出的线性系统§7-2多输入多输出的线性系统本章讨论机械或结构系统在随机激励作用下，激励—系统—响应三者之间的关系。系统有线性与非线性之分。大量工程问题，线性模型可得到逼真的结果。本课程只讨论线性系统问题。随机激励分两类：参数激励与非参数激励参数激励：系统本身的某些参数(如质量、刚度、阻尼等)随时间随机地变化而引起振动。非参数激励即由外界施加的激励。非参数激励又分为平稳的和非平稳的两类。本章研究常参数线性系统对平稳随机激励的响应当系统的激励(输入)是平稳过程时，由于常参数的假设，系统的响应(输出)也一定是平稳的。对于一个常参数线性系统，它往往可能在不同位置上同时受到激励，即有多个输入；其响应也可能有很多个，而且不同位置处的响应也不同。对于线性系统来说，多输入与多输出问题可以在单输入与单输出问题的基础上应用叠加原理得到解决对于确定性振动，激励与响应之间的关系，一般用微分方程来描述，方程的非齐次项是确定的，初始条件也是确定的，因此响应也是确定的。在随机振动中，一般激励与响应都必须用概率统计的方法来描述。在激励与系统特性已知的情况下，只能求出响应的一些统计特征，如期望（均值）、相关函数、功率谱密度、均方值等。7-1单输入单输出的线性系统假定常参数线性系统只受到一个输入x(t)的作用，其相应的响应(输出)为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(excitation)输入（激励）一、响应的均值对于输入的一个样本函数，由卷积积分公式，可得输出的一个样本函数()()()dytxth()()()dytxth设想对于输入中的每个样本函数，都可按上式写出其对应的输出的样本函数。于是，可得到输出的集合平均为:[()]()()dEYtEXthdhtXEtYE)()]([)]([xtXEtXE)]([)]([dhtYEx)()]([dehHj)()(0dhH)()0(输入与输出均值的关系式为：[()]=(0)YXEYtHdhtYEx)()]([dhH)()0(H(0)是一个常数，它表示输入X(t)与输出Y(t)中，频率ω=0这一成分（即直流分量）之间的传递关系。)()()0(txtyH直流分量上式表明，当输入是平稳过程时，输出的均值与输入的均值只差一个乘子H(0)。若输入的均值为零，则输出的均值也一定为零。此结论可以推广到多输入与多输出的情形。[()]=(0)YXEYtH二、响应的自相关函数输出过程Y(t)的自相关函数定义为:[()()]EYtYt上式为输出的自相关函数之间的关系式。该式说明，对于常参数线性系统，若激励是平稳随机过程，则响应的自相关函数与自然时间无关，也一定是平稳的随机过程。则响应的自相关函数可表示为:122112()[()()]=()()()]ddYXREYtYthhR对输出的自相关函数作傅立叶变换，便得到响应的自功率谱密度SY(ω)为三、响应的自功率谱密度函数变换积分次序，并重新排列jj122112()()()()()]dddYYXSRedehhR1221jj1122j2121ddY()XS()h()edh()eR()e()令ξ=τ-θ1+θ2，由维纳—辛钦关系式知，最后一个积分就是激励X(t)的自谱密度：j()()dXXSRe第二个积分就是脉冲响应函数h(θ2)的傅立叶变换，即频率响应函数H(ω)。2j22()()dHhe1221jj1122j2121ddY()XS()h()edh()eR()e()前两个积分的不同在于指数中的正负号的差别。经处理后得随机输入与输出的自谱密度关系式:上式是随机振动理论中一个极其重要的公式，指出了输入、输出与系统动态特性三者之间的关系。1221jj1122j2121ddY()XS()h()edh()eR()e()111)()()(dehHHj2()()()()()()YXXSHHSHS上式表明，若已知系统的增益因子|H(ω)|和输入的自谱密度SX(ω)，则可确定输出的自谱密度SY(ω)。事实上，若已知SX(ω)、|H(ω)|和SY(ω)三者中的任意两个，就可以确定第三个。此外，响应的自谱密度是与系统的相位因子无关的。2()()()()()()YXXSHHSHS已知响应的自谱密度SY(ω)，则可计算出响应的均方值E[Y2]:21[](0)()d2πYYEYRS四、响应的均方值将随机输入与输出的自谱密度关系式代入上式2221[]()()d2πYXEYHS注意：当均值为零时，均方值就等于方差。22YY在输入为理想白噪声的情况下，由于输入的自谱密度对于所有的频率都是常数，则响应的均方值公式可得到简化:220()d2πYSH只要计算出如下的广义积分I值，便可求得响应的均方值：2()dIH2221[]()()d2πYXEYHS由互相关函数的定义，可得激励与响应之间的互相关函数:()[()()]()()()d()[()()]d()()dXYXREXtYtEXthXthEXtXthR五、激励与响应的互相关函数常参数线性系统在受到平稳随机输入时，激励与响应之间的互相关函数正好等于脉冲响应函数与输入自相关函数的卷积()()()dytxth()()()dXYXRhR对互相关函数表达式作傅立叶变换，便可得到激励与响应之间的互功率谱密度。六、激励与响应的互谱密度jjj-()()d=[()()]d()()ddXYXYXSRehRdehRe与前面计算响应自谱相似的方法，将上式改写为：jj()()()d()d()XYXSheRe()()()dXYXRhR上式第一个积分是频率响应函数H(ω)，第二个积分就是激励X(t)的自谱密度SX(ω)()()()XYXSHS上式表明：输入与输出之间的互谱密度等于系统的频率响应函数与输入自谱密度函数的乘积。通过该式可完整地确定系统的频率特性H(ω)。jj()()()d()d()XYXSheRe由于H(ω)是复数，它可表示为:()()j()HAB则互谱密度可以表示为:()[()j()]()XYXSABS由于SX(ω)是实偶函数，则互谱密度函数可表示为：()()()()()arctg()XYXXYBSHSA上式表明：互谱密度的幅值等于系统的增益因子与输入自谱的乘积互谱密度的幅角又等于系统的相位因子。()()()()()arctg()XYXXYBSHSA七、相干函数22()()()()XYXYXYSSS系统输入与输出的谱相干函数（又称凝聚函数）可通过下式来定义：SX(ω)和SY(ω)皆为实函数，故相干函数必为实函数。可以证明，对于所有频率ω，相干函数满足以下不等式：20()1XY当输入与输出互不相关时，有RXY(τ)=0，从而互谱密度SXY(ω)=0，于是由定义知相干函数也等于零。对于线性系统，存在下列关系2()()()()()()YXXYXSHSSHS在线性系统的假设下，输入输出线性相关，有22222()()()()1()()()()()XYXXYXYXXSHSSSSHS输入输出互不相关时，相干函数的值等于0；输入输出线性相关时，相干函数等于1。相干函数的值在0与1之间。如果相干函数值大于零但小于1，为以下三种情况之一(1)联系输入X(t)和输出Y(t)的系统是非线性的(2)测量中有外界噪声干扰(3)输出Y(t)是输入X(t)和其它输入的综合输出。j0()txtxe如图所示的单输入线性系统，假定只在输出测量中混有噪声，则实测得到的输出Z(t)是真实输出Y(t)与噪声干扰N(t)之和。()()()ZtYtNtN(t)Z(t)H(ω)X(t)Y(t)只讨论一种存在噪声干扰的情况：假定X(t)与N(t)皆是均值为零的平稳随机过程，且N(t)与X(t)和Y(t)都是不相关的，则有：0)()()()(YNNYXNNXRRRR)()()()]()()()([)]()([)(XYXNXYXZRRRtNtXtYtXEtZtXER输入与实测输出之间的互相关函数：故有:)()(XYXZSS实测输出的自相关函数)()()()()()()}]()()}{()([{)]()([)(NYNNYYNYZRRRRRRtNtYtNtYEtZtZER实测输出的自谱密度)()()(NYZSSS输入X(t)与实测输出Z(t)的谱相干函数:()()()()()XZXYZYNSSSSS22()()()()XZXZXZSSS得到222Z2()()()1()(()()())XXXXNHSSHSS将互谱密度与自谱密度式代入上式，并结合下式2()()()()()()XYXYXSHSSHS上式表明：在有噪声干扰的情况下，输入与实测输出的谱相干函数将小于1。因此，对于线性系统，可借助相干函数值来判断干扰影响的大小。)()()()()()()(XXZXXYXZSSHSHSS此外：虽输出中含有干扰，但通过实测信号的互谱密度以及输入信号的自谱密度可以精确的获得系统的频响特征222Z2()()()1()(()()())XXXXNHSSHSS7.2多输入多输出的线性系统考虑某一具有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个脉冲响应函数：X1(t)X2(t)Xm(t)Y1(t)Y2(t)Yn(t)······对于n个输出，则共有n×m个脉冲响应函数，脉冲响应以矩阵形式可表示为：111212122212()()()()()()()()()()mmnnnmhthththththththththt矩阵中各元素均加以两个脚标第一个脚标k(k=1,2,…,n)表示k处的响应(输出)；第二个脚标i(i=1,2