(完整word版)随机振动分析报告

Good luck!-Alex-dreamer(南理工)

一：目的：根据abaqus爱好者提出的PSD随机振动分析，提出功率谱如何定义及如何加载？如果功率谱是加速度的平方，如何加载？如果在输入点施加载荷功率谱如何定义？本文将给出详细的分析过程。

二：随机振动基本概念

1. 随机振动的输入量和输出量都是概率统计值，因此存在不确定性。输入量为PSD （功率谱密度）曲线，分为加速度、速度、位移或者力的PSD曲线；最常见的是加速度PSD，常用语BASE MOTION基础约束加载。

2. 随机振动的响应符合正态分布，PSD实际上是随机变量的能量分布，也就是在不同频率上的方差值，反映不同频率处的振动能量，PSD曲线所围成的面积是随机变量总响应的方差值；

3. RMS为随机变量的标准方差，将PSD曲线包络面积开平方即为RMS。

4. 随机振动输出的位移、应力、应变等值都是对应不同频率的方差值（即PSD值），量纲为x^2，当然也可以输出这些变量的均方根值（即RMS值）；abaqus6.10以上版本可以直接在场变量里面输出设置。见下文。

5. 如果是单个激励源，定义为非相关性分析，如是多个激励源，则需要定义相关性参数。因此出现type=uncorrelated。

三：模型简介：

1）该模型很简单，是hypermesh中一个双孔模型。

2）网格划分在hypermesh中完成，保证了雅克比>0.7以及网格其它质量的要求。网格与几何具有较高的吻合度。

3）方案1（对应connect模型）：在上方两个孔采用全约束方式，且加载的功率谱PSD密度是加速度功率谱，也就是说基于BASE基础约束，进行随机振动

PSD分析。结果分析底部孔处某节点的结果响应。

4）方案2（对应connect模型）：在底部圆孔施加载荷force类型的功率谱PSD，与前者不同的是，这个不是基础施加PSD，而上某输入位置施加PSD。

5）重点和难点见一下详细介绍。

四：建模和导入在此略去。

1）导入模型；

施加约束

Fig 1 导入网格模型

2）分析布定义以及设置

第一步先进行频率提取分析，这里提取前20阶模态。采用lanczos法。

第二步进行随机振动分析布设置，低频5Hz开始，高频2000Hz截至。

Damping设置：开始模态阶数：1，结束模态阶数20，阻尼0.001。

Fig 3 模态分析设置

Fig 4 随机振动参数设置

Fig 5 随机振动参数设置

3）结果输出设置

我们关心均方根应力输出，那么输出变量是勾选RS即可。其他变量输出读者自己选择。

4）历史变量输出。假如我们关心底部某节点的输出。

Fig 6 均方根输出设置

Fig 7 历史输出设置5）边界条件设置

对基础点进行约束，也就是我们施加PSD的位置。

Fig 8 历史输出设置

6）如何编辑关键字？Model-→Edit Keywords---→random_force-m *Elastic

207000., 0.3

** BOUNDARY CONDITIONS

*Modal Damping

** OUTPUT REQUESTS

** FIELD OUTPUT: F-Output-2

*Output, field

*Node Output

A, AR, U, UR, V, VR

** HISTORY OUTPUT: H-Output-1

*Output, history

*Node Output, nset=output

A1, A2, A3, AR, AR1, AR2, AR3, RA1

RA2, RA3, RAR1, RAR2, RAR3, RU1, RU2, RU3 RUR1, RUR2, RUR3, RV1, RV2, RV3, RVR1, RVR2 RVR3, U1, U2, U3, UR, UR1, UR2, UR3

V1, V2, V3, VR, VR1, VR2, VR3

*End Step

6）结果分析

Fig 9 应力输出

Fig 10 均方根应力输出

Fig 11加速度响应输出

Fig 12 输出点加速度响应输出

Fig 13 模态阶数与频率关系输出

方案二：这次采用指定点输入激励PSD，来研究模型响应。原节点处固定方式不变，关键字有一定差异，详见下文。

Tip：我这里输入定义了节点集，因此引入INPUT命名。

1）分析布同上；输出点我变了，读者可以随意定义。

2）边界条件同上；

3）关键字编辑（对应connect-force模型）：

*Elastic

** BOUNDARY CONDITIONS

*Modal Damping

** OUTPUT REQUESTS

结果分析：

其它输出这里就不详细分析了。

机械故障诊断案例分析

六、诊断实例例1：圆筒瓦油膜振荡故障的诊断某气体压缩机运行期间，状态一直不稳定，大部分时间振值较小，但蒸汽透平时常有短时强振发生，有时透平前后两端测点在一周内发生了20余次振动报警现象，时间长者达半小时，短者仅1min左右。图1-7是透平1＃轴承的频谱趋势，图1-8、图1-9分别是该测点振值较小时和强振时的时域波形和频谱图。经现场测试、数据分析，发现透平振动具有如下特点。图1-7 1*轴承的测点频谱变化趋势图1-8 测点振值较小时的波形与频谱

图1-9 测点强振时的波形和频谱 (1)正常时，机组各测点振动均以工频成分）幅值最大，同时存在着丰富的低次谐波成分，并有幅值较小但不稳定的（相当于×）成分存在，时域波形存在单边削顶现象，呈现动静件碰磨的特征。 (2)振动异常时，工频及其他低次谐波的幅值基本保持不变，但透平前后两端测点出现很大的×成分，其幅度大大超过了工频幅值，其能量占到通频能量的75%左右。 (3)分频成分随转速的改变而改变，与转速频率保持×左右的比例关系。 (4)将同一轴承两个方向的振动进行合成，得到提纯轴心轨迹。正常时，轴心轨迹稳定，强振时，轴心轨迹的重复性明显变差，说明机组在某些随机干扰因素的激励下，运行开始失稳。 (5)随着强振的发生，机组声响明显异常，有时油温也明显升高。诊断意见：根据现场了解到，压缩机第一临界转速为3362r/min，透平的第一临界转速为8243r/min，根据上述振动特点，判断故障原因为油膜涡动。根据机组运行情况，建议降低负荷和转速，在加强监测的情况下，维持运行等待检修机会处理。生产验证：机组一直平稳运行至当年大检修。检修中将轴瓦形式由原先的圆筒瓦更改为椭圆瓦后，以后运行一直正常。例2：催化气压机油膜振荡某压缩机组配置为汽轮机十齿轮箱＋压缩机，压缩机技术参数如下：工作转速：7500r/min出口压力：轴功率：1700kW 进口流量：220m3 /min 进口压力：转子第一临界转速：2960r/min 1986年7月，气压机在运行过程中轴振动突然报警，Bently 7200系列指示仪表打满量程，轴振动值和轴承座振动值明显增大，为确保安全，决定停机检查。

随机振动名词解释

"脉冲响应函数" 英文对照 impulse response function; "脉冲响应函数" 在学术文献中的解释 1、h(t)是在初始时刻作用以单位脉冲而使单自由度系统产生的响应,所以称为脉冲响应函数.1·1·2频率响应函数H(ω)=1k-ω2m+iωcH(ω)是角频率为ω的单位简谐激励所引起的结构稳态简谐响应的振幅,称为频率响应函数,也称为转换函数文献来源 2、Yεi,jtt+s作为时间间隔s的一个函数,度量了在其他变量不变的情况下Yi,t+s对Yj,t的一个脉冲的反应,因此称为脉冲响应函数文献来源 "频率响应函数" 英文对照 frequency response function; "频率响应函数" 在学术文献中的解释 1、频率响应函数是指系统输出信号与输入信号的比值随频率的变化关系它是衡量高速倾斜镜工作性能的一个重要指标.通过抑制谐振峰可以改善高速倾斜镜的使用性能文献来源 2、经傅利叶变换,得到频域内的导纳(一般用速度导纳来表示)表达式 Hv(ω)=v(ω)F(ω)=jω-ω2M+jωC+K(2)H(ω)又称为频率响应函数文献来源 3、y（t）＝A0eiωty（t）＝iωA0eiωt（6）将（6）代入（3）得A0eiωt（RCiω＋1）＝Ajeiωt（7）和A0Aj＝1RCiω＋1＝U（iω）（8）U（iω）称为频率响应函数文献来源 "传递函数" 英文对照 transfer function of; transfer function; transfer function - noise; "传递函数" 在学术文献中的解释 1、由于传递函数的定义是两个拉普拉斯变换之比,所以使用时必须准确知道传递函数的类型,即,是位移、速度,还是加速度传递函数,才能避免出错文献来源 2、而传递函数的定义是两个分量之比为两个传感器之间优势波的传递函数.它给我们的启发是任取两个已知传感器组成一个传递函数通过分析传递函数的特征可以判断两个分量的优势波和非优势波文献来源

随机信号分析实验报告

一、实验名称微弱信号的检测提取及分析方法二、实验目的 1．了解随机信号分析理论如何在实践中应用 2．了解随机信号自身的特性，包括均值、方差、相关函数、频谱及功率谱密度等 3．掌握随机信号的检测及分析方法三、实验原理 1．随机信号的分析方法在信号与系统中，我们把信号分为确知信号和随机信号。其中随机信号无确定的变化规律，需要用统计特新进行分析。这里我们引入随机过程的概念，所谓随机过程就是随机变量的集合，每个随机变量都是随机过程的一个取样序列。随机过程的统计特性一般采用随机过程的分布函数和概率密度来描述，他们能够对随机过程作完整的描述。但由于在实践中难以求得，在工程技术中，一般采用描述随机过程的主要平均统计特性的几个函数，包括均值、方差、相关函数、频谱及功率谱密度等来描述它们。本实验中算法都是一种估算法，条件是N要足够大。 2．微弱随机信号的检测及提取方法因为噪声总会影响信号检测的结果，所以信号检测是信号处理的重要内容之一，低信噪比下的信号检测是目前检测领域的热点，而强噪声背景下的微弱信号提取又是信号检测的难点。噪声主要来自于检测系统本身的电子电路和系统外空间高频电磁场干扰等，通常从以下两种不同途径来解决 ①降低系统的噪声，使被测信号功率大于噪声功率。 ②采用相关接受技术，可以保证在信号功率小于噪声功率的情况下，人能检测出信号。对微弱信号的检测与提取有很多方法，常用的方法有：自相关检测法、多重自相法、双谱估计理论及算法、时域方法、小波算法等。对微弱信号检测与提取有很多方法，本实验采用多重自相关法。多重自相关法是在传统自相关检测法的基础上，对信号的自相关函数再多次做自相关。即令：式中，是和的叠加；是和的叠加。对比两式，尽管两者信号的幅度和相位不同，但频率却没有变化。信号经过相关运算后增加了信噪比，但其改变程度是有限的，因而限制了检测微弱信号的能力。多重相关法将当作x（t），重复自相关函数检测方法步骤，自相关的次数越多，信噪比提高的越多，因此可检测出强噪声中的微弱信号。

梁的振动实验报告

梁的振动实验报告实验目的改变梁的边界条件，对比分析不同边界条件，梁的振动特性（频率、振型等）。对比理论计算结果与实际测量结果。正确理解边界条件对振动特性的影响。实验内容对悬臂梁、简支梁进行振动特性对比，利用锤击法测量系统模态及阻尼比等。实验原理 1、固有频率的测定悬臂梁作为连续体的固有振动，其固有频率为：，其一、二、三、四阶时，简支梁的固有频率为：其一、二、三、四阶时，其中E为材料的弹性模量，I为梁截面的最小惯性矩，ρ为材料密度，A为梁截面积，l为梁的长度。试件梁的结构尺寸：长L=610mm, 宽b=49mm, 厚度h=8.84mm. 材料参数： 45＃钢，弹性模量E＝210 (GPa), 密度＝7800 (Kg/m3) 横截面积：A＝4.33*10-4 (m2), 截面惯性矩：J＝=2.82*10-9(m4) 则梁的各阶固有频率即可计算出。 2、实验简图

图1 悬臂梁实验简图

图2简支梁实验简图实验仪器本次实验主要采用力锤、加速度传感器、YE6251数据采集仪、计算机等。图3和图4分别为悬臂梁和简支梁的实验装置图。图5为YE6251数据采集仪。

图3 悬臂梁实验装置图图4 简支梁实验简图图5 YE6251数据采集分析系统实验步骤 1:"在教学装置选择"中，选择结构类型为"悬臂梁"，如果选择等份数为17，将需要测量17个测点。 2:本试验可采用多点激励，单点响应的方式，如果是划分为17等份,请将拾振点放在第5点。 3:请将力锤的锤头换成尼龙头,并将力通道的低通滤波器设置为1KHz,将拾振的加速度通道的低通滤波器设置为2KHz。 4:用力锤对第1点激振，对应的激励为f1,响应为1，平均3次，对应的数据为第1批数据，以此类推，测量完全部测点。 5:选择"教学装置模态分析和振型动画显示"，调入测量数据进行分

随机信号分析课程设计报告

随机信号分析课程设计报告标准化工作室编码[XX968T-XX89628-XJ668-XT689N]

随机信号分析课程设计报告题目学院信息电子技术专业电子信息工程班级 15级1班学籍号 1 姓名朱李伟指导教师刘文科信息电子技术学院 2018年6月18日

实验二随机过程的模拟与数字特征一、实验目的 1. 学习利用MATLAB模拟产生随机过程的方法。 2. 熟悉和掌握特征估计的基本方法及其MATLAB 实现。二、实验原理 1.正态分布白噪声序列的产生 MATLAB提供了许多产生各种分布白噪声序列的函数，其中产生正态分布白噪声序列的函数为randn。函数：randn 用法：x = randn(m,n) 功能：产生 m×n 的标准正态分布随机数矩阵。如果要产生服从N (,) 分布的随机序列，则可以由标准正态随机序列产生。如果X ~ N(0,1)，则N (,)。 2.相关函数估计 MATLAB提供了函数 xcorr用于自相关函数的估计。函数：xcorr 用法：c= xcorr (x,y) c= xcorr (x) c= xcorr (x,y ,'opition') c= xcorr (x, ,'opition') 功能：xcorr(x,y) 计算X (n ) 与Y (n)的互相关，xcorr(x)计算X (n )的自相关。 option 选项可以设定为： 'biased' 有偏估计。 'unbiased' 无偏估计。 'coeff' m = 0 时的相关函数值归一化为1。 'none' 不做归一化处理。 3.功率谱估计对于平稳随机序列X(n)，如果它的相关函数满足那么它的功率谱定义为自相关函数R X(m)的傅里叶变换：功率谱表示随机信号频域的统计特性，有着重要的物理意义。我们实际所能得到的随机信号的长度总是有限的，用有限长度的信号所得的功率谱只是真实功率谱的估计，称为谱估计或谱分析。功率谱估计的方法有很多种，这里我们介绍基于傅里叶分析的两种通用谱估计方法。（1）自相关法

振动实验报告剖析

振动与控制系列实验姓名：李方立学号：201520000111 电子科技大学机械电子工程学院

实验1 简支梁强迫振动幅频特性和阻尼的测量一、实验目的 1、学会测量单自由度系统强迫振动的幅频特性曲线。 2、学会根据幅频特性曲线确定系统的固有频率f 0和阻尼比。二、实验装置框图图3.1表示实验装置的框图图3-1 实验装置框图 K C X 图3-2 单自由度系统力学模型三、实验原理单自由度系统的力学模型如图3-2所示。在正弦激振力的作用下系统作简谐强迫振动，设激振力F 的幅值B 、圆频率ωo(频率f=ω／2π)，系统的运动微分方程式为：扫频信号源动态分析仪计算机系统及分析软件打印机或绘图仪简支梁振动传感器激振器力传感器质量块 M

或 M F x dt dx dt x d M F x dt dx n dt x d F Kx dt dx C dt x d M /2/222 22 2 222=++=++=++ωξωω （3-1）式中：ω—系统固有圆频率 ω =K/M n ---衰减系数 2n=C/M ξ---相对阻尼系数 ξ=n/ω F ——激振力 )2sin(sin 0ft B t B F πω== 方程①的特解，即强迫振动为： ) 2sin()sin(0?π?ω-=-=f A A x （3-2）式中：A ——强迫振动振幅 ? --初相位 2 0222024)(/ωωωn M B A +-= （3-3）式(3-3)叫做系统的幅频特性。将式(3-3)所表示的振动幅值与激振频率的关系用图形表示，称为幅频特性曲线(如图3-3所示)： 3-2 单自由度系统力学模型 3-3 单自由度系统振动的幅频特性曲线图3-3中，Amax 为系统共振时的振幅；f 0为系统固有频率，1f 、2f 为半功率点频率。振幅为Amax 时的频率叫共振频率f 0。在有阻尼的情况下，共振频率为： 2 21ξ-=f f a (3-4) 当阻尼较小时，0f f a =故以固有频率0f 作为共振频率a f 。在小阻尼情况下可得 01 22f f f -= ξ (3-5) 1f 、2f 的确定如图3-3所示： M X C K

转动设备常见振动故障频谱特征案例分析

转动设备常见振动故障频谱特征及案例分析一、不平衡转子不平衡是由于转子部件质量偏心或转子部件出现缺损造成的故障，它是旋转机械最常见的故障。结构设计不合理，制造和安装误差，材质不均匀造成的质量偏心，以及转子运行过程中由于腐蚀、结垢、交变应力作用等造成的零部件局部损坏、脱落等，都会使转子在转动过程中受到旋转离心力的作用，发生异常振动。转子不平衡的主要振动特征： 1、振动方向以径向为主，悬臂式转子不平衡可能会表现出轴向振动； 2、波形为典型的正弦波； 3、振动频率为工频，水平与垂直方向振动的相位差接近90度。案例：某装置泵轴承箱靠联轴器侧振动烈度水平13.2 mm／s，垂直11.8mm ／s，轴向12.0 mm／s。各方向振动都为工频成分，水平、垂直波形为正弦波，水平振动频谱如图1所示，水平振动波形如图2所示。再对水平和垂直振动进行双通道相位差测量，显示相位差接近90度。诊断为不平衡故障，并且不平衡很可能出现在联轴器部位。

解体检查未见零部件的明显磨损，但联轴器经检测存在质量偏心，动平衡操作时对联轴器相应部位进行打磨校正后振动降至2.4 mm／s。二、不对中转子不对中包括轴系不对中和轴承不对中两种情况。轴系不对中是指转子联接后各转子的轴线不在同一条直线上。轴承不对中是指轴颈在轴承中偏斜，轴颈与轴承孔轴线相互不平行。通常所讲不对中多指轴系不对中。不对中的振动特征： 1、最大振动往往在不对中联轴器两侧的轴承上，振动值随负荷的增大而增高；

2、平行不对中主要引起径向振动，振动频率为2倍工频，同时也存在工频和多倍频，但以工频和2倍工频为主； 3、平行不对中在联轴节两端径向振动的相位差接近180度； 4、角度不对中时，轴向振动较大，振动频率为工频，联轴器两端轴向振动相位差接近180度。案例：某卧式高速泵振动达16.0 mm／s，由振动频谱图(图3)可以看出，50 Hz（电机工频）及其2倍频幅值显著，且2倍频振幅明显高于工频，初步判定为不对中故障。再测量泵轴承箱与电机轴承座对应部位的相位差，发现接近180度。解体检查发现联轴器有2根联接螺栓断裂，高速轴上部径向轴瓦有金属脱落现象，轴瓦间隙偏大；高速轴止推面磨损，推力瓦及惰性轴轴瓦的间隙偏大。检修更换高速轴轴瓦、惰性轴轴瓦及联轴器联接螺栓后，振动降到A区。三、松动机械存在松动时，极小的不平衡或不对中都会导致很大的振动。通常有三种类型的机械松动，第一种类型的松动是指机器的底座、台板和基础存在结构松动，或水泥灌浆不实以及结构或基础的变形，此类松动表现出的振动频谱主要为1x。第二种类型的松动主要是由于机器底座固定螺栓的松动或轴承座出现裂纹引起，其振动频谱除1X外，还存在相当大的2X分量，有时还激发出1／2X和3X振动

北理工随机信号分析实验报告

本科实验报告实验名称：随机信号分析实验

实验一随机序列的产生及数字特征估计一、实验目的 1、学习和掌握随机数的产生方法。 2、实现随机序列的数字特征估计。二、实验原理 1、随机数的产生随机数指的是各种不同分布随机变量的抽样序列（样本值序列）。进行随机信号仿真分析时，需要模拟产生各种分布的随机数。在计算机仿真时，通常利用数学方法产生随机数，这种随机数称为伪随机数。伪随机数是按照一定的计算公式产生的，这个公式称为随机数发生器。伪随机数本质上不是随机的，而且存在周期性，但是如果计算公式选择适当，所产生的数据看似随机的，与真正的随机数具有相近的统计特性，可以作为随机数使用。 (0，1)均匀分布随机数是最最基本、最简单的随机数。(0，1)均匀分布指的是在[0，1]区间上的均匀分布，即 U(0，1)。实际应用中有许多现成的随机数发生器可以用于产生(0，1)均匀分布随机数，通常采用的方法为线性同余法，公式如下： )(m od ,110N ky y y n n -= N y x n n /= 序列{}n x 为产生的(0，1)均匀分布随机数。下面给出了上式的3组常用参数： 1、10 N 10,k 7==，周期7 510≈?； 2、（IBM 随机数发生器）31 16 N 2,k 23,==+周期8 510≈?； 3、（ran0）31 5 N 21,k 7,=-=周期9 210≈?；由均匀分布随机数，可以利用反函数构造出任意分布的随机数。定理 1.1 若随机变量 X 具有连续分布函数F X (x)，而R 为(0,1)均匀分布随机变量，则有 )(1R F X x -= 由这一定理可知，分布函数为F X (x)的随机数可以由(0,1)均匀分布随机数按上式进行变

随机信号分析大作业

随机信号分析实验报告信息25班 2120502123 赵梦然

作业题三：利用Matlab 产生一个具有零均值、单位方差的的高斯白噪声随机序列X(n)，并通过一脉冲响应为 (0.8)(0)0 n n h n else =≥??? 的线性滤波器。 (1) 产生一个具有零均值、单位方差的的高斯白噪声随机序列X(n)，检验其一维概率密度函数是否与理论相符。 (2) 绘出输入输出信号的均值、方差、自相关函数及功率谱密度的图形，讨论输出信号服从何种分布。 (3) 试产生在[-1,+1]区间均匀分布的白噪声序列，并将其替换高斯白噪声通过上述系统。画出此时的输出图形，并观察讨论输出信号服从何种分布。作业要求 (1) 用MATLAB 编写程序。最终报告中附代码及实验结果截图。 (2) 实验报告中必须有对实验结果的分析讨论。提示： (1) 可直接使用matlab 中已有函数产生高斯白噪声随机序列。可使用hist 函数画出序列的直方图，并与标准高斯分布的概率密度函数做对比。 (2) 为便于卷积操作，当N 很大时，可近似认为h(N)=0。卷积使用matlab 自带的conv 函数。 (3) 分析均值、方差等时，均可使用matlab 现有函数。功率谱密度和自相关函数可通过傅里叶变换相互获得。傅里叶变换使用matlab 自带的fft 函数。 (4) 作图使用plot 函数。

一、作业分析：本题主要考察的是加性高斯白噪声相关问题，因此构造一个高斯白噪声十分重要，故在本题中使用randn函数随机生成一个个符合高斯分布的数据，并由此构成高斯白噪声；而且由于白噪声是无法完全表示的，故此根据噪声长度远大于信号长度时可视为高斯白噪声，构造了一个长度为2000的高斯白噪声来进行试验。二、作业解答：（1）matlab程序为： x-1000:1:1000; k=1*randn(1,length(x));% 生成零均值单位方差的高斯白噪声。 [f,xi]=ksdensity(x);%利用ksdensity函数估计样本的概率密度。 subplot(1,2,1); plot(x,k); subplot(1,2,2); plot(xi,f); 实验结果为：

机械振动实验报告分析

实验三：简谐振动幅值测量一、实验目的 1、了解振动位移、速度、加速度之间的关系。 2、学会用压电传感器测量简谐振动位移、速度、加速度幅值二、实验仪器安装示意图三、实验原理由简谐振动方程：)sin()(?ω-=t A t f 简谐振动信号基本参数包括：频率、幅值、和初始相位，幅值的测试主要有三个物理量，位移、速度和加速度，可采取相应的传感器来测量，也可通过积分和微分来测量，它们之间的关系如下：根据简谐振动方程，设振动位移、速度、加速度分别为x 、v 、a ，其幅值分别为X 、V 、A ： )sin(?ω-=t X x )cos()cos(?ω?ωω-=-==t V t X x v )sin()sin(2?ω?ωω-=--==t A t X x a 式中：ω——振动角频率 ?——初相位所以可以看出位移、速度和加速度幅值大小的关系是：X V A X V 2ωωω===,。振动信号的幅值可根据位移、速度、加速度的关系，用位移传感器或速度传感器、加速度传感器进行测量，还可采用具有微积分功能的放大器进行测量。在进行振动测量时，传感器通过换能器把加速度、速度、位移信号转换成电信号，经过放大器放大，然后通过AD 卡进行模数转换成数字信号，采集到的数字信号为电压变化量，通过软件在计算机上显示出来，这时读取的数值为电压值，通过标定值进行换算，就可计算出振动量的大

小。 DASP 通过示波调整好仪器的状态（如传感器档位、放大器增益、是否积分以及程控放大倍数等）后，要在DASP 参数设置表中输入各通道的工程单位和标定值。工程单位随传感器类型而定，或加速度单位，或速度单位，或位移单位等等。传感器灵敏度为K CH （PC/U ）（PC/U 表示每个工程单位输出多少PC 的电荷，如是力，而且参数表中工程单位设为牛顿N ，则此处为PC/N ；如是加速度，而且参数表中工程单位设为m/s 2 ，则此处为PC/m/s 2 ）； INV1601B 型振动教学试验仪输出增益为K E ；积分增益为K J （INV1601 型振动教学试验仪的一次积分和二次积分K J =1）； INV1601B 型振动教学试验仪的输出增益：加速度：K E = 10(mV/PC) 速度：K E = 1 位移：K E = 0.5 则DASP 参数设置表中的标定值K 为： )/(U mV K K K K J E CH ??= 四、实验步骤 1、安装仪器把激振器安装在支架上，将激振器和支架固定在实验台基座上，并保证激振器顶杆对简支梁有一定的预压力（不要露出激振杆上的红线标识），用专用连接线连接激振器和INV1601B 型振动教学试验放大仪的功放输出接口。把带磁座的加速度传感器放在简支梁的中部，输出信号接到 INV1601B a 加速度。 2、打开INV1601B 型振动教学试验仪的电源开关，开机进入DASP2006 标准版软件的主界面，选择单通道按钮。进入单通道示波状态进行波形示波。 3、在采样参数设置菜单下输入标定值K 和工程单位m/s 2 ,设置采样频率为4000Hz ，程控倍数1倍。 4、调节INV1601B 型振动教学试验仪频率旋钮到40Hz 左右，使梁产生共振。 5、在示波窗口中按数据列表进入数值统计和峰值列表窗口，读取当前振动的最大值。 6、改变档位v (mm /s )、d (mm )进行测试记录。 7、更换速度和电涡流传感器分别测量a (m /s 2 )、v (mm /s )、d (mm )。

(完整word版)随机振动分析报告

Alex-dreamer制作PSD：（可以相互传阅学习，但是鄙视那些拿着别人成果随意买卖！）PSD随机振动应用领域很广，比如雷达天线，飞机，桥梁，天平，地面，等等行业。虽然现在对这方面公开资料很少，但是我相信以后会越来越多，发展的越来越成熟。学术的浪潮总体是向前的，不会因为几个大牛保密自己的成果就会阻止我们对PSD研究，因此结合我的经验和爱好，我研究了一下两种PSD加载分析。我标价的原则是含金量大小和花费我的时间以及我的经验值，如果你觉得值，就买；不值就不要下了。因为我始终认为：士为知己者死，女为悦己者容。算是互相尊重。如果你得到这份资料，那就祝你好运！ Good luck!-Alex-dreamer(南理工) 一：目的：根据abaqus爱好者提出的PSD随机振动分析，提出功率谱如何定义及如何加载？如果功率谱是加速度的平方，如何加载？如果在输入点施加载荷功率谱如何定义？本文将给出详细的分析过程。二：随机振动基本概念 1. 随机振动的输入量和输出量都是概率统计值，因此存在不确定性。输入量为PSD （功率谱密度）曲线，分为加速度、速度、位移或者力的PSD曲线；最常见的是加速度PSD，常用语BASE MOTION基础约束加载。 2. 随机振动的响应符合正态分布，PSD实际上是随机变量的能量分布，也就是在不同频率上的方差值，反映不同频率处的振动能量，PSD曲线所围成的面积是随机变量总响应的方差值； 3. RMS为随机变量的标准方差，将PSD曲线包络面积开平方即为RMS。 4. 随机振动输出的位移、应力、应变等值都是对应不同频率的方差值（即PSD值），量纲为x^2，当然也可以输出这些变量的均方根值（即RMS值）；abaqus6.10以上版本可以直接在场变量里面输出设置。见下文。 5. 如果是单个激励源，定义为非相关性分析，如是多个激励源，则需要定义相关性参数。因此出现type=uncorrelated。三：模型简介： 1）该模型很简单，是hypermesh中一个双孔模型。 2）网格划分在hypermesh中完成，保证了雅克比>0.7以及网格其它质量的要求。网格与几何具有较高的吻合度。 3）方案1（对应connect模型）：在上方两个孔采用全约束方式，且加载的功率谱PSD密度是加速度功率谱，也就是说基于BASE基础约束，进行随机振动 PSD分析。结果分析底部孔处某节点的结果响应。 4）方案2（对应connect模型）：在底部圆孔施加载荷force类型的功率谱PSD，与前者不同的是，这个不是基础施加PSD，而上某输入位置施加PSD。

adams振动分析实例中文版

1.问题描述研究太阳能板展开前和卫星或火箭分离前卫星的运行。研究其发射振动环境及其对卫星各部件的影响。 2.待解决的问题在发射过程中，运载火箭给敏感部分航天器部件以高载荷。每个航天器部件和子系统必学设计成能够承受这些高载荷。这就会带来附加的质量，花费高、降低整体性能。更好的选择是设计运载火箭适配器（launch vehicle adapter）结构。这部分，将设计一个（launch vehicle adapter）的隔离mount，以在有效频率范围降低发射震动传到敏感部件的部分。关心的敏感部件在太阳能板上，对70-100HZ的输入很敏感，尤其是垂直于板方向的。三个bushings将launch vehicle adapter和火箭连接起来。Bushing的刚度和阻尼影响70-100HZ范围传递的震动载荷。所以设计问题如下：找到运载火箭适配器系统理想刚度和阻尼从而达到以下目的：传到航天器的垂直加速度不被放大； 70-100HZ传递的水平加速度最小。 3.将要学习的 Step1——build：在adams中已存在的模型上添加输入通道和振动执行器来时系统振动，添加输出通道测量响应。 Step2——test:定义输入范围并运行一个振动分析来获得自由和强迫振动响应。 Step3——review：对自由振动观察模态振型和瞬态响应，对强迫振动，观察整体响应动画，传递函数。 Step4——improve：在横向添加力并检查传递加速度，改变bushing的刚度阻尼并将结果作比较。添加频域测量供后续设计研究和优化使用。

3.1需创建的东西：振动执行器、输入通道、输出通道完全非线性模型打开模型在install dir/vibration/examples/tutorial satellite 文件夹下可将其复制到工作木录。加载Adams/vibration模块：Tools/ plugin Manager. 仿真卫星模型：仿真看其是否工作正常，仿真之前关掉重力，这个仿真太阳能板在太空中的位置。关掉重力：Settings——Gravity ；仿真：tool面板——simulation ，设置仿真时间是15s，步长为500；点击，将停在仿真后mode 返回最初的模型状态：点击，把重力打开，这时模型回到振动分析准确的发射状态。

随机信号实验报告

随机信号分析实验报告目录随机信号分析 (1) 实验报告 (1) 理想白噪声和带限白噪声的产生与测试 (2) 一、摘要 (2) 二、实验的背景与目的 (2) 背景： (2) 实验目的: (2) 三、实验原理 (3) 四、实验的设计与结果 (4) 实验设计: (4) 实验结果: (5) 五、实验结论 (12) 六、参考文献 (13) 七、附件 (13) 1

理想白噪声和带限白噪声的产生与测试一、摘要本文通过利用MATLAB软件仿真来对理想白噪声和带限白噪声进行研究。理想白噪声通过低通滤波器和带通滤波器分别得到低通带限白噪声和帯通带限白噪声。在仿真的过程中我们利用MATLAB工具箱中自带的一些函数来对理想白噪声和带限白噪声的均值、均方值、方差、功率谱密度、自相关函数、频谱以及概率密度进行研究，对对它们进行比较分析并讨论其物理意义。关键词：理想白噪声带限白噪声均值均方值方差功率谱密度自相关函数、频谱以及概率密度二、实验的背景与目的背景：在词典中噪声有两种定义：定义1：干扰人们休息、学习和工作的声音，引起人的心理和生理变化。定义2：不同频率、不同强度无规则地组合在一起的声音。如电噪声、机械噪声，可引伸为任何不希望有的干扰。第一种定义是人们在日常生活中可以感知的，从感性上很容易理解。而第二种定义则相对抽象一些，大部分应用于机械工程当中。在这一学期的好几门课程中我们都从不同的方面接触到噪声，如何的利用噪声，把噪声的危害减到最小是一个很热门的话题。为了加深对噪声的认识与了解，为后面的学习与工作做准备，我们对噪声进行了一些研究与测试。实验目的: 了解理想白噪声和带限白噪声的基本概念并能够区分它们，掌握用MATLAB 或c/c++软件仿真和分析理想白噪声和带限白噪声的方法，掌握理想白噪声和带限白噪声的性质。

随机振动试验报告

随机振动试验报告高等桥梁结构试验报告讲课老师: 张启伟(教授) 姓名: 史先飞学号: 1232627 试验报告 1 试验目的 1.过试验进一步加深对结构模态分析理论知识的理解; 2.熟悉随机振动试验常用仪器的性能与操作方法; 3.复习和巩固随机振动数据测量和分析中有关基本概念; 4.掌握通过多点激振、单点拾振的方法，利用DASP2005软件进行模态分析的基本操作步骤。

2 试验仪器和设备 1. ZJY-601振动与控制教学实验仪系统(ZJY-601A型振动教学实验仪、激励锤、YJ9-A型压电型加速度传感器等)。 2. DASP 16通道接口箱。 3. 装有“DASP2005智能数据采集和信号分析系统”软件的PC机。 4. 有关设备之间的联接电缆。 3 试验原理 3.1模态叠加原理 N自由度线性振动系统的运动微分方程是一组耦合的方程组: 引入模态矩阵Φ和模态坐标(广义坐标或主坐标)q，使X= Φq。如果阻尼矩阵能对角化，方程组即可解耦: 解耦后的第i个方程为: 可见，采用固有振型描述振动的模态坐标后，N自由度线性振动系统的振动响应可以表示为N阶模态响应的叠加。 3.2实模态理论实模态理论建立在无阻尼的假设基础上。在实模态理论中，模态频率就是系统的无阻 ,尼模态固有频率错误～未找到引用源。;而固有振型矩阵中的各元素都是实数，它们之间i 的相位差是0?或180?。系统在P点激励，l点测量的频响函数为:

K,,式中，称为频率比，，为模态固有频率。当，则: ,,,,,/,,,iiiiiMi 取频响函数矩阵的一列或一行，如第P列，就可确定振动系统的全部动力特性(模态参数)。 3.3伪实模态理论某些有阻尼振动系统有时会出现与实模态一样的实数振型，而非复数振型，但其模态 2,,,,,1固有频率为，具有这种性质的振动系统的模态称为伪实模态。伪实模态理diii 论仅适应于阻尼矩阵可解耦，即可采用固有振型矩阵正交化模态称为伪实模态。在伪实模态下，各测点的相位差都是0?或180?。伪实模态理论仅适应于阻尼矩阵可解耦，即可采用固有振型矩阵正交化的情况。一般情况下，阻尼矩阵对角化的充要条件为: 上式也是有阻尼振动系统方程解耦的充要条件。总之，H(ω)建立了模态参数与频响函数的关系。因此，利用实验测出的H(ω) 值，即可计算出系统的模态参数。根据频响函数的互易定理及模态理论，只需 H(ω)矩阵的一列(或一行)即可求出全部模态参数。

随机信号分析实验报告(基于MATLAB语言)

随机信号分析实验报告 ——基于MATLAB语言姓名： _ 班级： _ 学号：专业：

目录实验一随机序列的产生及数字特征估计 (2) 实验目的 (2) 实验原理 (2) 实验内容及实验结果 (3) 实验小结 (6) 实验二随机过程的模拟与数字特征 (7) 实验目的 (7) 实验原理 (7) 实验内容及实验结果 (8) 实验小结 (11) 实验三随机过程通过线性系统的分析 (12) 实验目的 (12) 实验原理 (12) 实验内容及实验结果 (13) 实验小结 (17) 实验四窄带随机过程的产生及其性能测试 (18) 实验目的 (18) 实验原理 (18) 实验内容及实验结果 (18) 实验小结 (23) 实验总结 (23)

实验一随机序列的产生及数字特征估计实验目的 1.学习和掌握随机数的产生方法。 2.实现随机序列的数字特征估计。实验原理 1.随机数的产生随机数指的是各种不同分布随机变量的抽样序列（样本值序列）。进行随机信号仿真分析时，需要模拟产生各种分布的随机数。在计算机仿真时，通常利用数学方法产生随机数，这种随机数称为伪随机数。伪随机数是按照一定的计算公式产生的，这个公式称为随机数发生器。伪随机数本质上不是随机的，而且存在周期性，但是如果计算公式选择适当，所产生的数据看似随机的，与真正的随机数具有相近的统计特性，可以作为随机数使用。 (0,1)均匀分布随机数是最最基本、最简单的随机数。(0,1)均匀分布指的是在[0,1]区间上的均匀分布， U(0,1)。即实际应用中有许多现成的随机数发生器可以用于产生(0,1)均匀分布随机数，通常采用的方法为线性同余法，公式如下： y0=1，y n=ky n(mod N) ? x n=y n N 序列{x n}为产生的(0,1)均匀分布随机数。定理1.1若随机变量X 具有连续分布函数F x(x)，而R 为(0,1)均匀分布随机变量，则有 X=F x?1(R) 2.MATLAB中产生随机序列的函数 (1)(0,1)均匀分布的随机序列函数：rand 用法：x = rand(m,n) 功能：产生m×n 的均匀分布随机数矩阵。 (2)正态分布的随机序列函数：randn 用法：x = randn(m,n) 功能：产生m×n 的标准正态分布随机数矩阵。如果要产生服从N(μ,σ2)分布的随机序列，则可以由标准正态随机序列产生。 (3)其他分布的随机序列分布函数分布函数二项分布binornd 指数分布exprnd 泊松分布poissrnd 正态分布normrnd 离散均匀分布unidrnd 瑞利分布raylrnd 均匀分布unifrnd X2分布chi2rnd 3.随机序列的数字特征估计对于遍历过程，可以通过随机序列的一条样本函数来获得该过程的统计特征。这里我们假定随机序列X(n)为遍历过程，样本函数为x(n)，其中n=0,1,2,……N-1。那么，

机械振动实验报告

《机械振动基础》实验报告（2015年春季学期）姓名学号班级专业机械设计制造及其自动化报告提交日期2015.05.07 哈尔滨工业大学

报告要求 1.实验报告统一用该模板撰写，必须包含以下内容： (1)实验名称 (2)实验器材 (3)实验原理 (4)实验过程 (5)实验结果及分析 (6)认识体会、意见与建议等 2.正文格式：四号字体，行距为1.25倍行距； 3.用A4纸单面打印；左侧装订； 4.报告需同时提交打印稿和电子文档进行存档，电子文档由班长收齐，统一发送至：liuyingxiang868@https://www.360docs.net/doc/3312109694.html,。 5.此页不得删除。评语：教师签名：年月日

实验一报告正文一、实验名称：机械振动的压电传感器测量及分析二、实验器材 1、机械振动综台实验装置(压电悬臂梁) 一套 2、激振器一套 3、加速度传感器一只 4、电荷放大器一台 5、信号发生器一台 6、示波器一台 7、电脑一台 8、NI9215数据采集测试软件一套 9、NI9215数据采集卡一套三、实验原理信号发生器发出简谐振动信号，经过功率放大器放大，将简谐激励信号施加到电磁激振器上，电磁激振器振动杆以简谐振动激励安装在激振器上的压电悬臂梁。压电悬臂梁弯曲产生电流显示在示波器上，可以观测悬臂梁的振动情况；另一方面，加速度传感器安装在电磁激振器振动杆上，将加速度传感器与电荷放大器连接，将电荷放大器与数据采集系统连接，并将数据采集系统连接到计算机（PC机）上，操作NI9215数据采集测试软件，得到机械系统的振动响应变化曲线，可以观测电磁激振器的振动信号，并与信号发生器的激励信号作对比。实验中的YD64-310型压电式加速度计测得的加速度信号由DHF-2型电荷放大器后转变为一个电压信号。电荷放大器的内部等效电路如图1所示。 q

振动分析实例

旋转机械诊断监测管理系统（TDM）在电厂的应用摘要：介绍了应用旋转机械诊断监测管理系统（TDM）的硬件及软件组成；深入分析了#4汽轮机组9瓦轴振异常的原因，获取包括转速、波德图、频谱、倍频的幅值和相位等故障特征数据，从而为专业的故障诊断人员提供数据及专业的图谱，协助机组诊断维护专家深入分析机组运行状态，并成功处理了9瓦的轴振异常。关键词：应用旋转机械诊断监测管理系统（TDM），组成，异常振动，分析，解决 The Application of the Turbine Diagnosis Management (TDM) on Shanxi Zhangshan Electric Power co., Ltd Li Gang He Xiao Ming Kou Delin (The College of Power and Mechanical Engineering Wuhan University Wuhan 430072) Abstract: Introduce the hardware and software of the Turbine Diagnosis Management (TDM). Analysis the reasons of #9 bearing’s abnormal vibration of unit 4.Receives the characteristic data of the speed, Bode diagram, frequency phase, mult-frequency’s value and phase.Offers the professional data ,charts to the experts. Helps the experts diagnosis deeply the status of the unit 4. And solve the problem successfully. Key words：Turbine Diagnosis Management (TDM), Composition, abnormal vibration, Analysis, solution 引言汽轮机轴系监测系统(TSI)可以对汽轮机轴系参数起到基本的监测和安全保护作用，但TSI 缺少对机组振动数据的深入挖掘，使得许多振动方面的问题停留在表面,如在机组冲转、在负荷变化，主、调汽阀门进行切换和单/顺阀切换等工况变化时振动的分析研究。而旋转机械诊断监测管理系统（TDM）则填补了此项功能。它的主要作用在于对机组运行过程中的数据进行深入分析，获取包括转速、振动波形，频谱、倍频的幅值和相位等故障特征数据，从而为专业的故障诊断人员提供波德图、频谱图、瀑布图、级联图、轴心轨迹等专业的数据及图谱，协助机组诊断维护专家深入分析机组轴系运行状态，解决机组在实际运行中遇到的问题。 1. TDM 的硬件及软件的组成漳山电厂采用北京英华达公司生产的EN8001旋转机械振动监测分析故障诊断专家系统EN8001系统是由硬件系统和软件系统组成，硬件系统主要由下位高速智能数据采集、信息处理、信息数据存储管理系统和服务器、上位机工程师站及附件构成，硬件系统采用积木式模块化的结构，配置灵活，上下位硬件系统通过工业以太网络集成。系统软件由三大部分构成：数据采集软件，数据库软件和分析诊断软件构成。数据采集软件负责数据采集，它能自动识别机组的运行状态，如开停机、升降速及正常或异常状态，并根据机组的状态进行数据采集。在稳定运行状态下，数据硬件采集系统以定时方式进行采集，而在升降速状态下则根据转速的变化进行采集。数据库软件负责数据的存储，它由升降速数据库、历史数据库及事件数据库等组成，它根据机组的不同状态把有关数据存到不同的数据库中，以便于后续分析。分析诊断软件主要用于对各种数据进行在线或离线分析，以判断机组的运行状态并能自动给出机组故障原因和处理 1

ABAQUS软件随机振动分析 final

ABAQUS软件随机振动分析在工程中，结构一般需要对它进行随机振动分析。典型的例子是：通过机床的振动响应分析进行机床的结构设计，通过对结构的地震响应分析。在电子产品设计中，ABAQUS软件不仅仅能对电子产品进行冲击、热场、加工等过程进行数值模拟，还可以对电子产品在随机振动下产品的响应性能做出很好预测，以优化产品设计。本例题就某电子产品在随机激励作用下的响应结构为例，采用如下图所示的简化模型，分析在特定随机激励（如图2）中，分析该结构的响应。图 1 某电子产品结构简化图图2 随机激励的谱分布载荷边界条件为：四个底座固支，并在分析过程中，受到随机激励。需要分析整个结构在运动过程中的响应。启动ABAQUS/CAE，在Start Session对话框中，选择Create Model Database按钮。

一导入模型由于IGES文件给的是实体模型，我们在计算中产用shell模型，所以我们需要通过 ABAQUS/CAE中对shell的编辑功能对模型进行修改。导入IGES文件成Shell格式。 1．在主菜单选择File -＞Import-＞Part, 进入Import Part对话框。选择相应的 IGES文件，点击按钮。 2．在弹出的Create Part From IGES File 对话框中，如下图，对话框的Topology选择Shell选项，Name选项填写random。二利用CAE编辑修改模型在主菜单选择Shape -＞Shell-＞Remove Face,用鼠标点击选择模型中的面，选上之后面会变红色，点击鼠标中键，就可以去掉该面。重复操作，得到下图模型。

随机信号分析资料报告习题

随机信号分析习题一 1. 设函数???≤>-=-0 , 0 ,1)(x x e x F x ，试证明)(x F 是某个随机变量ξ的分布函数。并求下列概率：)1(<ξP ，)21(≤≤ξP 。 2. 设),(Y X 的联合密度函数为 (), 0, 0 (,)0 , other x y XY e x y f x y -+?≥≥=? ?，求{}10,10<<<

8. 两个随机变量1X ，2X ，已知其联合概率密度为12(,)f x x ,求12X X +的概率密度？ 9. 设X 是零均值，单位方差的高斯随机变量，()y g x =如图，求()y g x =的概率密度 ()Y f y \ 10. 设随机变量W 和Z 是另两个随机变量X 和Y 的函数 22 2 W X Y Z X ?=+?=? 设X ，Y 是相互独立的高斯变量。求随机变量W 和Z 的联合概率密度函数。 11. 设随机变量W 和Z 是另两个随机变量X 和Y 的函数 2() W X Y Z X Y =+?? =+? 已知(,)XY f x y ，求联合概率密度函数(,)WZ f z ω。 12. 设随机变量X 为均匀分布，其概率密度1 ,()0X a x b f x b a ?≤≤? =-???，其它（1）求X 的特征函数,()X ?ω。（2）由()X ?ω，求[]E X 。 13. 用特征函数方法求两个数学期望为0，方差为1，互相独立的高斯随机变量1X 和2X 之和的概率密度。 14. 证明若n X 依均方收敛，即 l.i.m n n X X →∞ =，则n X 必依概率收敛于X 。 15. 设{}n X 和{}n Y (1,2, )n =为两个二阶矩实随机变量序列，X 和Y 为两个二阶矩实随机变量。若l.i.m n n X X →∞ =，l.i.m n n Y Y →∞ =，求证lim {}{}m n m n E X X E XY →∞→∞ =。