随机振动功率谱密度值
功率谱密度

振动台在使用中经常运用的公式1、 求推力(F )的公式F=(m 0+m 1+m 2+ ……)A …………………………公式(1) 式中:F —推力(激振力)(N )m 0—振动台运动部分有效质量(kg ) m 1—辅助台面质量(kg )m 2—试件(包括夹具、安装螺钉)质量(kg )A — 试验加速度(m/s 2)2、 加速度(A )、速度(V )、位移(D )三个振动参数的互换运算公式 2.1 A=ωv ……………………………………………………公式(2) 式中:A —试验加速度(m/s 2)V —试验速度(m/s ) ω=2πf (角速度) 其中f 为试验频率(Hz )2.2 V=ωD ×10-3………………………………………………公式(3) 式中:V 和ω与“2.1”中同义D —位移(mm 0-p )单峰值2.3 A=ω2D ×10-3 ………………………………………………公式(4) 式中:A 、D 和ω与“2.1”,“2.2”中同义 公式(4)亦可简化为:A=D f ⨯2502式中:A 和D 与“2.3”中同义,但A 的单位为g1g=9.8m/s 2所以: A ≈D f ⨯252,这时A 的单位为m/s 2 定振级扫频试验平滑交越点频率的计算公式 3.1 加速度与速度平滑交越点频率的计算公式f A-V =VA28.6 ………………………………………公式(5)式中:f A-V —加速度与速度平滑交越点频率(Hz )(A 和V 与前面同义)。
3.2 速度与位移平滑交越点频率的计算公式DV f DV 28.6103⨯=- …………………………………公式(6) 式中:D V f -—加速度与速度平滑交越点频率(Hz )(V 和D 与前面同义)。
3.3 加速度与位移平滑交越点频率的计算公式f A-D =DA ⨯⨯23)2(10π ……………………………………公式(7) 式中:f A-D — 加速度与位移平滑交越点频率(Hz ),(A 和D 与前面同义)。
随机振动(振动频谱)计算(Random Vibration)

Random Vibration1. 定义1.1 功率谱密度当波的频谱密度乘以一个适当的系数后将得到每单位频率波携带的功率,这被称为信号的功率谱密度(power spectral density, PSD)。
功率谱密度谱是一种概率统计方法,是对随机变量均方值的量度。
1.2 均方根均方根(RMS)是指将N项的平方和除于N后,开平方的结果。
均方根值也是有效值,如对于220交流电,示波器显示的有效值或均方根值为220V。
2. 加速度功率谱密度2.1 单位加速度单位:m/s^2或g加速度功率谱密度单位:(m/s^2)^2/Hz或g^2/HzHz单位为:1/s,所以加速度功率谱密度单位也可写为:m^2/s^32.2功率谱密度函数功率谱密度函数曲线的纵坐标是(g²/Hz)。
功率谱曲线下的面积就是随机加速度的总方差(g²):σ²= ∫Φ(f)df其中:Φ(f)........功率谱密度函数σ ............. 均方根加速度3. 计算示例随机振动100-2000HZ,功率谱密度为0.01g^2/Hz,则其加速度峰值计算如下:σ²=0.01*(2000-100)=19σ=4.36g峰值加速度不大于3倍均方根加速度:13.08g4、SAE J 1455 随机振动要求4.1功率谱图4.1.1 Vertical axis4.1.2 Transverse axis4.1.3 Longitudinal axis4.2 Vertical axis加速度计算功率谱曲线下的面积:σ²=(40-5)0.016+0.5*(500-40)*0.016=4.24σ=2.06g峰值加速度不大于3倍均方根加速度:6.18g5. FGE随机振动要求5.1功率谱图5.2 要求在工作状态,振动频率范围:10Hz-1000Hz,振动方向:X、Y、Z三轴,试验时间:每轴各8h,加速度均方根为33.9 m/s²(3.46g)。
随机振动分析

程序支持多个PSD基础激励,但是不考虑其关联性,也就 是程序不支持计算不同PSD激励的关联性。
3.随机振动分析步骤
(4)计算结果 程序支持三个方向的位移,速度和加速度; 因为每个方向的计算结果是统计结果,因此不 能使用一般的方法进行合并。
如果需要输出应力和应变,可用的应力结果只有名义应变和应力, 剪切应变和应力,等效应力。
4.工程实例:电路板的随机振动计算
1.随机振动分析简介
什么是随机振动分析
– 基于概率的谱分析. – 典型应用如火箭发射时结构承受的载荷谱,每次发射的谱不同,但统 计规律相同.
1.随机振动分析简介
• 和确定性谱分析不同,随机振动不能用瞬态动力学分析代 替. • 应用基于概率的功率谱密度分析,分析载荷作用过程中的 统计规律
什么是PSD?
3.随机振动分析步骤
(2)分析设置
Analysis Settings > Output Controls (1)默认情况下,位移,速度和加速度响应是输出的; (2)为了不输出速度或加速度响应,可以将输出选项设置 为No。
3.随机振动分析步骤
(3)载荷和支撑条件
1)支撑条件必须在模态分析中进行设置; 2)PSD分析中只支持PSD基础激励,包括 -PSD加速度 -PSD G加速度 -PSD速度 -PSD位移
• PSD是激励和响应的方差随频率的变化。 – PSD曲线围成的面积是响应的方差. – PSD的单位是 方差/Hz (如加速度功率谱的单位是 G2/Hz). – PSD可以是位移、速度、加速度、力或压力.
2.随机振动分析理论
(1)随机振动激励分布规律 因为随机振动激励被假设为服从高斯正态分布,因此没有计算发生 概率为100%的结构响应。 在实际工程中,分布式激励更加普遍; 此外,高sigma激励发生的概率很低;
随机振动功率谱密度

701z010203040506070800.0020.0040.0060.0080.010.0120.0140.016频率(Hz)功率谱密度功率谱密度函数图(汉宁窗)1020304050607080-65-60-55-50-45-40-35-30-25-20-15频率(Hz)功率谱密度(d B )功率谱密度函数图(汉宁窗)经过matlab 频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.1378m/s2(70km/h,z 方向,第一次试验,前排)0.1378010203040506070800.511.522.5-3频率(Hz)功率谱密度频率加权后功率谱密度函数图(汉宁窗)701y010203040506070801234567-3频率(Hz)功率谱密度功率谱密度函数图(汉宁窗)1020304050607080-70-65-60-55-50-45-40-35-30-25-20频率(Hz)功率谱密度(d B )功率谱密度函数图(汉宁窗)经过matlab 频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.0164m/s2(70km/h,y 方向,第一次试验,前排)010203040506070800.511.522.53-5频率(Hz)功率谱密度频率加权后功率谱密度函数图(汉宁窗)701x010203040506070800.20.40.60.811.21.41.61.8-3频率(Hz)功率谱密度功率谱密度函数图(汉宁窗)01020304050607080-70-65-60-55-50-45-40-35-30-25频率(Hz)功率谱密度(d B )功率谱密度函数图(汉宁窗)经过matlab 频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.0070m/s2(70km/h,x 方向,第一次试验,前排)010203040506070800.511.522.533.5-6频率(Hz)功率谱密度频率加权后功率谱密度函数图(汉宁窗)702经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.0078m/s2(70km/h,x方向,第2次试验,前排)经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.0193m/s2(70km/h,y方向,第2次试验,前排)经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.1393m/s2(70km/h,z方向,第2次试验,前排)703经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.0035m/s2(70km/h,x方向,第1次试验,后排)经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.0038m/s2(70km/h,y方向,第1次试验,后排)经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.027m/s2(70km/h,z方向,第1次试验,后排)经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.0052m/s2(70km/h,x方向,第2次试验,后排)经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.0072m/s2(70km/h,y方向,第2次试验,后排)经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.0296m/s2(70km/h,z方向,第2次试验,后排)经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.0062m/s2(60km/h,x方向,第1次试验,前排)经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.0144m/s2(60km/h,y方向,第1次试验,前排)经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.1216m/s2(60km/h,z方向,第1次试验,前排)经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.0079m/s2(60km/h,x方向,第2次试验,前排)经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.0174m/s2(60km/h,y方向,第2次试验,前排)经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.1172m/s2(60km/h,z方向,第2次试验,前排)经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.0032m/s2(60km/h,x方向,第1次试验,后排)经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.0063m/s2(60km/h,y方向,第1次试验,后排)经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.0267m/s2(60km/h,z方向,第1次试验,后排)经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.0029m/s2(60km/h,x方向,第2次试验,后排)经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.0058m/s2(60km/h,y方向,第2次试验,后排)经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.0231m/s2(60km/h,z方向,第2次试验,后排)经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.0058m/s2(50km/h,x方向,第1次试验,前排)经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.0170m/s2(50km/h,y方向,第1次试验,前排)经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.1186m/s2(50km/h,z方向,第1次试验,前排)502(可疑数据)经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.0122m/s2(50km/h,x方向,第2次试验,前排)经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.0302m/s2(50km/h,y方向,第2次试验,前排)经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.1820m/s2(50km/h,z方向,第2次试验,前排)经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.0030m/s2(50km/h,x方向,第1次试验,后排)经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.0049m/s2(50km/h,y方向,第1次试验,后排)经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.0220m/s2(50km/h,z方向,第1次试验,后排)经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.0025m/s2(50km/h,x方向,第2次试验,后排)经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.0042m/s2(50km/h,y方向,第2次试验,后排)经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.0214m/s2(50km/h,z方向,第2次试验,后排)经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.0081m/s2(40km/h,x方向,第1次试验,前排)经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.0102m/s2(40km/h,y方向,第1次试验,前排)经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.1141m/s2(40km/h,z方向,第1次试验,前排)402(可疑数据)经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.0092m/s2(40km/h,x方向,第2次试验,前排) 最大值=0.0468比值=5.2011<9经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.0088m/s2(40km/h,y方向,第2次试验,前排) 最大值=0.0489比值=5.4457<9经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.1254m/s2(40km/h,z方向,第2次试验,前排)403经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.0029m/s2(40km/h,x方向,第1次试验,后排)经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.0045m/s2(40km/h,y方向,第1次试验,后排)经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.0182m/s2(40km/h,z方向,第1次试验,后排)404经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.0020m/s2(40km/h,x方向,第2次试验,后排)经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.0031m/s2(40km/h,y方向,第2次试验,后排)经过matlab频率加权法,利用功率谱密度函数计算得到加权加速度均方根值0.0162m/s2(40km/h,z方向,第2次试验,后排)。
[物理]随机振动--第7章-功率谱密度
![[物理]随机振动--第7章-功率谱密度](https://img.taocdn.com/s3/m/4a8ea9cdcf2f0066f5335a8102d276a2002960f0.png)
第7章功率谱密度函数7.1 自相关的物理意义及其傅里叶变换7.2 自功率谱密度函数及其性质7.3 窄带随机过程与宽带随机过程7.4 互功率谱密度函数及其性质7.5 共相谱、正交谱和相干函数7.1 自相关的物理意义及其傅里叶变换自相关函数的物理意义可以表达现在的波形与时间坐标平移后的波形之间的相似程度表达随机过程两个不同截口处的两个随机变量之间的相关程度自相关函数与原始信号具有相同的周期(频率)、衰减率(阻尼)动态特性可用来检测随机过程中是否含有周期成分,或者其信号特征自功率谱计算的依据自相关函数既包含了一个随机过程间隔时间的相关程度和依赖性,同时也包含了能量大小的信息。
不过要注意,相关性再也不是象相关系数那样能够用-1到1这样的数来表示相关大小了自相关函数的性质1:⑴自相关函数是偶函数x x R E X t X t E X t X t R自相关函数的性质2:⑵周期平稳过程的自相关函数也是周期函数,其周期与过程的周期相同。
x x R T E X t X t T E X t X t R自相关函数的性质3:⑶τ=0时的自相关函数就是均方值20x x xR E X t X t R E X t X t⑷如果随机过程不是周期过程,则:22222222 01000x x xxx x xxx xxxx x C R R R R时,随机变量与它自身是完全相关的时,两个随机变量之间将不再相关前提:不是周期函数若,则 2lim x xR自相关函数的性质4:⑸自相关函数是一个有界函数22222 110xxx x xxxR R一般τ越大,则两时刻的随机变量X(t1)和X(t1+τ)之间的相关性愈差。
τ↑,Rx(τ)↓。
自相关函数的性质5:一、自功率谱密度函数二、互功率谱密度函数自相关函数Rx(τ)描述“平均功率”随时差τ的变化→“平均功率”的时间结构。
功率谱密度S x (f):描述“平均功率”在频域(谱域)的分布→频率结构。
功率谱密度

功率谱密度谱是一种概率统计方法,是对随机变量均方值的量度。
一般用于随机振动分析,连续瞬态响应只能通过概率分布函数进行描述,即出现某水平响应所对应的概率。
功率谱密度是结构在随机动态载荷激励下响应的统计结果,是一条功率谱密度值—频率值的关系曲线,其中功率谱密度可以是位移功率谱密度、速度功率谱密度、加速度功率谱密度、力功率谱密度等形式。
数学上,功率谱密度值—频率值的关系曲线下的面积就是方差,即响应标准偏差的平方值。
谱是个很不严格的东西,常常指信号的Fourier变换,是一个时间平均(time average)概念功率谱的概念是针对功率有限信号的(能量有限信号可用能量谱分析),所表现的是单位频带内信号功率随频率的变换情况。
保留频谱的幅度信息,但是丢掉了相位信息,所以频谱不同的信号其功率谱是可能相同的。
有两个重要区别:1。
功率谱是随机过程的统计平均概念,平稳随机过程的功率谱是一个确定函数;而频谱是随机过程样本的Fourier变换,对于一个随机过程而言,频谱也是一个“随机过程”。
(随机的频域序列)2。
功率概念和幅度概念的差别。
此外,只能对宽平稳的各态历经的二阶矩过程谈功率谱,其存在性取决于二阶局是否存在并且二阶矩的Fourier变换收敛;而频谱的存在性仅仅取决于该随机过程的该样本的Fourier变换是否收敛。
热心网友回答提问者对于答案的评价:谢谢解答。
频谱分析(也称频率分析),是对动态信号在频率域内进行分析,分析的结果是以频率为坐标的各种物理量的谱线和曲线,可得到各种幅值以频率为变量的频谱函数F(ω)。
频谱分析中可求得幅值谱、相位谱、功率谱和各种谱密度等等。
频谱分析过程较为复杂,它是以傅里叶级数和傅里叶积分为基础的。
功率谱是个什么概念?它有单位吗?随机信号是时域无限信号,不具备可积分条件,因此不能直接进行傅氏变换。
一般用具有统计特性的功率谱来作为谱分析的依据。
功率谱与自相关函数是一个傅氏变换对。
功率谱具有单位频率的平均功率量纲。
随机振动-试验人员必须了解的参数及设置

随机振动-试验人员必须了解的参数及设置江苏省电子信息产品质量监督检验研究院谢杰一.简述近年来,随机振动试验在我院所有振动试验中的比例越来越高,原因有三:1、科学进步,此类设备的软件大量普及,一般只需在原来的电磁振动台加上一套控制软件及配套设备就可实行。
2、企业随着国际标准的大量采用,许多振动试验都采用随机振动。
3、随机振动相对传统的正弦振动有着无法比拟的优点,它能模拟各种实际运输条件下可能遇到的振动情况,如模拟公路运输,模拟铁路运输,模拟海运运输等等。
本文主要介绍对于试验人员来说必须了解的随机振动参数及设置要求。
二.随机振动数据上图是某一随机振动试验后的试验数据,对于试验人员来说,必须了解其中的一些参数含义。
曲线中,横坐标是频率,纵坐标是PSD,一般简称为频谱曲线。
PSD:Power spectrum density 功率谱密度PSD单位有二种:g2/Hz,(m2/Hz)2/Hz,二者之间换算:1 g2/Hz=96(m2/Hz)2/Hz PSD是随机振动中的重要参数,可理解为每频率单位中所含振动能量的大小,其值越大,相对应的频率段振幅值会变大,在试验中提高最低频率的PSD 值可明显感觉到振幅增大。
频谱曲线的特点:1、它是对数坐标,主要是为了表述画线方便。
2、它有一条平线或多条平线及斜线组成,平线和斜线之间首尾相连组成。
3、试验条件中,PSD值不变的是平线,用+dB/oct表示向上的斜线,用- dB/oct 表示向下的斜线。
如-3 dB/oct 表示每增加一倍频率,PSD值下降一半。
频谱曲线中,中间一条是设定曲线,上面二条和下面二条是设备的保护及中断线,附加在中间设定值上的变化曲线是振动台实际控制曲线。
三.频率的选择频率是随机振动的另一个重要参数,其单位是Hz,频率的选择一般与实践使用范围有关。
例如:海运试验条件频率较低,一般从1~100Hz,而且低频PSD 值较大,随机振动的感觉像乘海轮,振幅大,频率低。
中国公路和京沪铁路运输、欧洲和日本实际采集随机振动功率谱密度曲线及数据

0.181
5Hz~200Hz的Grms
0.167
注1:低频率进行试验时,位移峰峰值可能超过25.4mm;
注2:5Hz~200Hz适用于小型样品的振动试验。
0.000 028 03:36:00
30 0.253 0.233
频率 100H0z
等级 3 0.006 0.01 0.01
/ / 0.0 003 / 0.0 003 0.00 015 / 0.00 015 01:12:00 10 0.415 0.395
距离 L,km
时间,h:min
L≤200
0:15
200<L≤500
0:30
500<L≤1 000
1:00
1 000<L≤1 500
1:30
1 500<L≤2 000
2:00
2 000<L≤2 500
2:30
2500<L
3:00
GB/T 4857.23—××××
11
1.06
12:00
1.00
注:推荐每个运输轴向至少振动30分钟。
D.2 ISO 13355:2016 实际采集日本 20 吨、半载、弹簧减震汽车的随机振动。
D.2.1 图D.2是日本实际采集的功率谱密度曲线。
0.1
功率谱密度 g^2/Hz
0.01
0.001
0.0001
0.00001
图D.2 日本实际采集随机振动PSD曲线
gbt48572310d13基于12小时运输试验数据缩短试验时间和提高试验等级放大倍数的关系见表d2表d2试验时间与psd放大倍数试验时间psd等级放大倍数0
附录D (资料性附录) 欧洲和日本实际采集随机振动功率谱密度曲线及数据
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
模态分析主要目的:
获得结构的固有频Biblioteka ,可避免共振现象的发生当外界激励力的频率等于振动系统的固有频率时, 系统发生共振现象。
模态是结构所固有的一种特性,它与本身的材料 特性,约束方式,结构形状有关,而与输入的加载形 式等激励形式无关。
单自由度系统频响函数分析 粘性阻尼系统
•强迫振动方程及其解
m x c x kx f
3.2.随机振动
离散系统回归方程:
y回 kx b
偏差:
y (kx b)
偏差平方期望:
E(2) E[( y kx b) 2 ]
求偏导:
E (2 ) 0 k
E (2 ) 0 b
偏差平方差期望值最小
2 2 2 2 [E(2 )]min y k 2 x y (1 xy )
X (iw)e iwt d ]dt
T0
lim
1 1 X (iw)[ 2T0 2
T0
lim
1 X (iw) X (i )d 2T0
1 2 X (iw) d 2T0
T0
lim
sx ( )dd
Sx(w)为功率谱密度
•解的形式(s为复数)及拉氏变换:
x Xe
st
..
.
(ms2 cs k ) x(s) f (s)
自由振动
m x c x kx 0
ms cs k 0
2
..
.
s1, 2 0 j 0 1 2
实部:衰减因子,反映系统阻尼 虚部:有阻尼系统的固有频率
1 S x ( ) 2
-
R ( )e i d
随机响应 输入自相关系数
Gin ( ) Rx ( )
经傅里叶变换,时域转换为频域
1 Gin ( ) S x ( ) 1 Gin ( ) S x ( ) 2 2
- -
i R ( )e d R ( )e i d
T0
T0 2 T0 2
1 X (-i ) X (i )eiw d 2T0
2
- T0
lim
X (i ) 2T0
e i d
-
S ( x)e i d
自相关函数是功率谱密度 函数的傅里叶变换
Rx ( )
-
S ( )ei d
由傅里叶变换可知:
阻尼比 范围(0-1) 内为欠阻尼 无阻尼固有频率
c /(2 km)
0 k / m
传递函数
H ( s)
x( s ) 1 2 f ( s) m s cs k
频率响应函数
s jw
H ( ) 1 m 2 jwc k
无阻尼时
H ( )
1 m 2 k
以电阻为例,单位电阻在T0时间,做的功
E( x )
2
T0 2 T0 2
x 2 (t ) d t
当T0较长时,平均功率为:
lim
T0
T0
1 T0 1 T0
T0 2 T0 2
x 2 (t )dt
lim
T0 2 T0 2
1 x(t )dt[ 2
X (i )eiwt d ]dt
系统传递函数 频率响应函数
H ( )
Gout () H () Gin ()
x(t ) x(t )d
iwt t t x(t )[ x(i )eiwt eii d ]dt d
T T0 0
1 2T0 0
T T0 0 2 2 T T0 0 2 2
lim lim
T0
T0 2 T0 2
T0 1 iwt x(i )ei [ 2 x ( t ) e dt]d T0 2T0 2
一、振动系统
定则系统 振动系统 随机系统 定则系统:
参量的变化规律可以用时间确定函数表述的振动系统
随机系统:
无法用时间的确定函数描述,而只能用有关统计性描述的 振动系统
二、随机振动
激励 随机激励(激励不能用确定描述函数表示) 定则激励
随机振动:受随机激励的系统称为随机振动,即使定则系 统 ,受随机激励也是随机振动 。
两时间段相隔
则:Rxy ( ) E[ x(t ) y(t )]
Rx ( ) E[ x(t ) y(t )]
样本函数时间的平均
Rxy ( ) lim Rx ( ) lim
T0
1 T0 1 T0
T0 2 T0 2
x(t ) y(t )d x(t ) x(t )d
如齿轮轴虽然是定则系统,但受随机激励是随机振动
随机激励 随机响应
齿轮轴
三、动力学分析
模态分析 动力学分析
谐响应分析
随机振动分析
3.1、模态分析
模态分析经典定义:将线性时不变系统 振动微分方程组中的物理坐标变换为模 态坐标,使方程组解耦,成为一组以模 态坐标及模态参数描述的独立方程,坐 标变换的变换矩阵为振型矩阵,其每列 即为各阶振型。
工程中把随机过程处理成零平均值
xy (t1 , t2 )
x (t1 , t2 )
E ( x1 y 2 )
x y
1
2
E ( x1 x2 )
x x
1
2
分子成为互相关函数和自相关函数:
Rxy (t1, t2 ) E( x1 y1 ) E[ x(t1 ) y(t2 )]
Rx (t1 , t2 ) E( x1 x2 ) E[ x(t1 ) x(t2 )]
式中
xy
E{[ x E ( x)][y E ( y)]}
x y
两个随机变量可以看做两个不同随机过程取出的
xy (t1 , t2 )
E{[ x E ( x)][y E ( y)]}
x y
式中:x1 x(t1 ), y2 y(t2 ) 若两个随机过程从同一个随机过程取 出 E{[ x1 E ( x1 )][x2 E ( x2 )]} x (t1 , t2 ) x1 x2
傅里叶函数的时移性质
1 X (t ) 2 π
x(i )e it d
1 X (t ) 2 π
x (i )e i ( t ) d
Rx ( ) lim x lim
T T0 0
1 T0 0
T T0 0 2 2 T T0 0 2 2
T0
T0 2 T0 2
功率谱密度(PSD)(Power Spectral Density) 单位频带上的功率,当波的频谱密度乘以一个适 当的系数后将得到每单位频率波携带的功率。
由傅里叶变换可知:
X (i ) x(t )e
it
dt
1 it X (t ) x(i )e d 2 π