机械振动测试系统的设计【含CAD图纸、说明书】
机械振动测试系统的设计1 开题依据1.1 选题背景用试验方法测量机械的振动量(如位移、速度和加速度等)和系统特征参数(如固有频率、阻尼、振型等),以及振动环境的模拟等,都属于振动测试。研究机械振动时通常采用理论分析和测试两种手段。通过测试可验证理论分析计算的正确性,提供所需的修正依据。20世纪80年代以来,振动测试仪器有了显著的进步,如传递函数分析仪、实时频率分析仪和快速傅里叶分析仪的相继应用,并与电子计算机相结合,为振动测试和测试结果的分析处理提供了方便的条件,从而也进一步推动了振动理论的研究和发展。系统的振动特性也可以应用激光全息照相法拍下实物或模型在振动时的全息照片,根据全息照片中的干涉条纹图案来分析。1.2 本课题意义 对于机械振动测试系统目前大多数的厂家还不具备先进的技术,装置大多从国外引进,从而设备成本较高。为了满足生产需要而又不需要投资太多的资金,我们需要不断学习国外的技术并不断地创新,研制出自己的产品。因此,我选择这个课题,对我来讲,意义非常深刻,也是一次不错的学习及自我能力提高的机会。2 文献综述2.1调研针对这次机械振动测试的题目,我特别到工厂进行调研。我了解到的机械振动测试的主要原理,方法和基本的过程,机械振动的研究可归结为机械系统的激励、响应和振动特性三个方面的问题。统称为机械阻抗曲线。在已知其中两个方面的情况下可求第三方面的问题。测得的机械阻抗数据通常以幅频特性曲线和相频特性曲线、实部和虚部频率特性曲线或幅相频率特性曲线(Nyquist图)3种形式表达,与之相对应,振动测试的基本内容包括:已知激励和系统的振动特性情况下求响应,即振动量的测量;已知激励和响应的情况下求系统的振动特性,即系统特征参数的测定,也称参数识别;已知系统的振动特性和响应的情况下求激励,即环境预测,即用时域描述的振动量时,这种测试称为振动环境模拟试验。 3.振动测试方案3.1 振动测试方案1的提出和组成原理图振动系统动态测试仪器的设计过程是:首先按图1所示的流程图布置好振动系统动态测试仪器的硬件系统;然后,在个人计算机上开发相应的应用软件;最后,振动系统动态测试与模态识别的全过程,都是在个人计算机上通过软面板对话框和适当的控制按钮来完成。信号发生器被测机械振动系统输入输出个人计算机系统及动态特性测试分析软 件X(n)振动传感器及信号调理电路Y1(n)被测振动系统的模态参数和数学模型的输出激振器振动传感器及信号调理电路数据采集卡DAQ,包含: 多路开关 放大器 采样保持器S/H A/D转换器YN(n) 图1 振动系统动态测试仪器原理图3.2 激振方式、激振器和振动传感器的布置激振方式可分为正弦稳态扫频和宽频带激振两大类。通常采用“正弦稳态扫频激振”,其优点是激振能量集中,信噪比高,对于线性系统的动态特性测试具有很高的测量精度。首先,根据测试分析要求及振动系统的机械结构特点,确定激振点的位置和测量点的数量N及各测量点的位置。对于空间结构系统,往往需要测量1个点的3个方向,这时应选用三轴振动传感器。安装激振器时,必须注意激振力的正确施加。在安装测振传感器时,应保证传感器能正确感受被测体的振动。小型振动系统必须选用微型传感器,以避免传感器附加质量影响振动系统动态特性。3.3 振动系统动态测试应用软件的总体研究 根据对前面振动系统动态测试虚拟仪器原理图的分析可知,振动系统动态测试应用软件的设计流程分为如下步骤: (1)计算机发出指令控制“信号发生器”产生激励信号,再通过“激振器”施加于被测振动系统。激励信号的数字序列为X(n)。 (2) 等待被测振动系统稳定一定时间,然后计算机再发出指令同步采集数据序列X(n)和Y1(n),Y2(n), ,YN(n) 。采样时间间隔或采样频率应符合采样定理及频率分辨率的要求。 (3)利用按本原理设计的专用软件分析处理X(n)和Y1(n),Y2(n),YN(n), 就可得出被测振动系统的模态参数和相应数学模型。显然,整套应用软件应包括两大部分,即操作控制部分和数据处理部分。由于目前计算机控制技术和虚拟仪器技术的发展,已足以解决此应用软件中的操作控制部分,且鉴于篇幅有限,所以笔者对应用软件的操作控制部分不再作进一步讨论,下面重点研究应用软件的数据处理部分。3.4 数据处理专用软件的设计原理3.4.1 滤波方法 由于数据序列在采集过程中难免会受到各种干扰, 所以必须采取有效的滤波算法对采集得到的数据序列Y1(n),Y2(n),YN(n)进行修正。通常应采用相关滤波算法,因为相关滤波可以有效地把真正由X(n)通过振动系统产生的信号从Y1(n),Y2(n),YN(n)中分离出来。将经过滤波算法处理后的数据序列记为YN1(n),YN2(n),YNN(n) ,n=1,2, ,D 。D为采样点数量。3.4.2 频率响应函数的确定 将确定性信号看成随机信号,采用离散数字随机信号的相关谱分析原理。 XK(k) = FFT X(n) YK(k) = FFT YN(n) SX(k) = XK(k)2 /D SY(k) = YK(k)2 /D SXY(k) = XK*(k)YK(k) /D XK*(K)为XK(k)的共轭复数。 为减少数据处理误差,要对SX(k),SY(k),SXY(k)采取平滑处理算法。分段时,相邻两段重叠50的效果最佳,即将记录段X(n),YN(n)分成多个样本,分别进行上述运算后再求平均值。VXY2(k) = SXY(k)2 / SX(k) SY(k), 相干系数VXY越接近1效果越好。H(k) =SXY(k) / SX(k), 由变量代换关系“k/(DTS) =f(TS为采样时间间隔)”可得到测试出的频率响应函数H(f) 。 对每一个数据序列Y1(n),Y2(n),YN(n)均用上述方法处理,就可以得到测试出来的N个频率响应函数H1(f),H2(f),HN(f) 。3.4.3 拟合求解模态参数机械振动系统模态参数识别,又称为曲线拟合,即采用最小二乘法将测试所得的频率响应值与系统模型值进行曲线拟合。优化准则就是使实测的频率响应值与理论数学模型对应值的总均方误差E极小。最小二乘曲线拟合过程中存在的有关问题,如方程组的病态问题,可根据数值分析理论,采用正交多项式的方法加以解决。通常机械振动系统的“理论数学模型”可用有理分式形式的传递函数来表示,即 H(s) = N(s)/D(s) = a0sm+a1sm-1+am-1s+am/ sn+b1sn-1+bn-1s+bn上式中,mn,n=2q,q为系统的自由度数。一般选定一个具体q值,令m=2q-1,然后循环拟合,直到满足E控制精度Eg。 对于每一个频率响应函数,设由D个采样点得到的D个测试频率响应值记为HC(f),对应的D个理论模型值记为HL(f),则总方差为E=(k=1,D)| HC(f)HL(f)|2。这是非线性优化问题,为使用方便和保证精度,可转化为线性优化问题。令s=jw=j2f,误差e=N(f) D(f)HC(f),则总方差转化为E=|N(f) D(f)HC(f)|2。此时,总方差函数为待求系数的线性函数,最小二乘解已变成线性优化问题。首先通过最小二乘拟合识别出有理分式的系数ai(i=1,2,m)和bi(i=1,2,n),然后由特征方程D(s)0解得极点sr及sr*(r=1,2,q),由留数公式ArH(s)(s-sr)| S=Sr”求出留数Ar及Ar* 。则传递函数可表示为H(s)r=1 , N Ar/(s-sr) + Ar*/(s-sr*)最后以sr= -r wr +j wr (1- r2)1/2,sr为共轭复数,求得各阶模态固有频率wr和阻尼比r 。对于线性系统,模态固有频率和阻尼比是系统的总体固有参数,理论上应该不随测量点而改变。依次对N个测量点的频率响应函数进行曲线拟合,可得到N组模态固有频率和阻尼比,再求出平均值作为系统的模态固有频率和阻尼比,进而,再求出各阶振型r 。3.4.4 其它输出振动系统动态测试所得的系统固有的各阶模态参数和用模态参数表示的各测量点,对应于输入点的数学模型。振动系统动态测试数据处理的全过程,可以用目前流行的程序设计语言VC+或Java2开发成专用的数据处理程序软件包,然后通过相关程序语言接口方法嵌入到用LabVIEW或LabWindows平台开发的振动系统动态测试虚拟仪器中。采用本虚拟仪器,可以很容易通过计算机多次重复振动系统动态测试的整个过程,进一步提高振动系统动态特性的测试精度。利用动态测试得到的数学模型,通过Matlab6.1仿真环境(或自己编制程序)可求出激振点处各种输入信号各测量点处的输出响应信号。3.5振动测试方案2的组成根据课题的需要,以及对实际情况的考察和比较,我们设计出系统的方案,如下图所示: 91073扫频信号源数采分析仪计算机系统及分析软件打印机或绘图仪电机调压器124568111213141516 系统方案设计图 振动测试与控制实验台由弹性体系统(包括简支梁、悬臂梁、薄壁圆板、单自由度系统、二自由度系统、多自由度系统模型)配以主动隔振、被动隔振用的空气阻尼减震器、单式动力吸振器、复式动力吸振器等组成。是完成振动与振动控制等近30个实验的试验平台 1、激振系统 激振系统包括: DH1301正弦扫频信号源 JZ-1型接触式激振器 JZF-1型非接触式激振器 偏心电动机、调压器 力锤(包括测力传感器) 2、测振系统DH5922动态采集分析仪 MT-3T型磁电式振动速度传感器 DH103压电式加速度传感器 WD302电涡流位移传感器 测力传感器3、动态采集分析系统 信号调理器 数据采集仪 计算机系统(或笔记本电脑) 控制与基本分析软件 模态分析软件系统软硬件配置合理、完善,能完成几乎全部与振动相关的试验,可完成各种模态试验(小到梁、板、构件,大到房屋、桥梁等),可完成频率耦合、结构修改、灵敏度分析等动力学试验;4 研究路线和技术方法振动量的测量测量机械系统某些选定点上的振幅(位移、速度和加速度)、频率、相位、振动的时间历程和频谱等。这种测量通常在机械系统的工作状态下进行,以了解其实际振动状况。也可将其转换为数字量后分析处理。对某些精密和大型机械设备的振动监控和诊断所作的测量也属这种性质。 对模拟量可直接分析处理,振动量测量按振动信号和转换方式可分为电测法、光测法和机械测振法,其中以电测法应用最为广泛。图1为一个较完整的振动量电测系统。测振传感器(拾振器)将机械振动量转换为与它成比例的电量。 测试结果的分析和数据处理 测试结果所获得的原始数据有两种表现形式:一种是模拟量,常用的测振传感器有发电型(如压电式、电动式和磁电式等)和电参数变化型(如电感式、电容式、电阻式和涡流式等)两类。以保证输出电信号的幅值和相位均不失真。不同类型的传感器需要配接不同类型的中间测量变换装置(图2)。最基本的要求是:在其工作频率范围内的幅频特性平坦、相频特性呈线性关系,中间测量变换装置对传感器输出的电信号进行前置变换(电阻抗变换)、微积分运算、放大、调制和解调等,以便驱动后接的分析或显示、记录设备。分析设备完成对信号的频率分析。显示、记录设备给出振动信号(经过分析的或未经过分析的)的波形,并用数字或模拟方式指示出测量结果,应校准温度的影响等。以便于储存、分析信号和进行数据处理。 系统特征参数的测定主要是应用机械阻抗测试技术,以获得机械阻抗数据(有时亦称频率响应数据),从而得到系统的特征参数如固有频率、阻尼、刚度、质量和振型等;还可通过模态分析求取系统在各阶模态下的特征参数,既模态参数。这一测试过程称为模态参数识别。这种测定通常在机械系统的非工作状态或模型试验情况下进行,以求全面了解其动态响应特性。尤其是测振传感器必须进行定期校准和定度。若在工作状态下进行,则常称为在线识别。在机械阻抗测试技术中,施加的激励有简谐、瞬态和随机3种类型,试验根据不同试验对象按相应的试验规范进行,故机械阻抗测试也相应地分3类。 振动环境模拟试验研究或考核试验对象在强度、寿命和功能方面的抗振性。简谐激励又有单点激励和多点激励两种形式。这种模拟试验分为周期性振动试验、随机振动试验和冲击试验 3种。周期性振动试验一般采用耐共振、耐扫频和耐预定频率试验 3种形式。在进行振动寿命试验时,为了缩短试验时间常采用提高振动量级的办法,即幅频特性;测定振动响应与激励间的相位差随激励频率变化的关系,即强化试验。提高的程度,即强化系数, 简谐激励机械阻抗的测试 以简谐力作为激励并保持其幅值恒定,应根据试件的振动响应特性和疲劳强度分析来考虑。试验根据不同试验对象按相应的试验规范进行,并用模拟振动试验机来实现。 测试系统的校准和定度为了保证测试结果的可靠性和测试精度,对所使用的仪器,尤其是测振传感器必须进行定期校准和定度。以求全面了解其动态响应特性。在进行重要的或特殊的试验前,这种测定通常在机械系统的非工作状态或模型试验情况下进行,常直接对整套测试系统进行现场校准和定度。测试系统最基本的校准项目包括灵敏度、频率响应和线性度。既模态参数。此外,从而得到系统的特征参数如固有频率、阻尼、刚度、质量和振型等;还可通过模态分析求取系统在各阶模态下的特征参数,根据需要还可进行某些特殊的校准,如所测振级变化范围大时,应校准动态线性范围;高温下测试时,以便于储存、分析信号和进行数据处理。应校准温度的影响等。测振传感器的校准在测试系统的校准中具有特别重要的意义。校准方法主要有两种:一是绝对校准法,分析设备完成对信号的频率分析。二是比较校准法。无论是对测振传感器或对组成测试系统的仪器和对整个测试系统的校准,最基本的要求是:在其工作频率范围内的幅频特性平坦、相频特性呈线性关系,不同类型的传感器需要配接不同类型的中间测量变换装置(图2)。以保证输出电信号的幅值和相位均不失真。 测试结果的分析和数据处理测试结果所获得的原始数据有两种表现形式:一种是模拟量,如电压和电流等;一种是数字量。对不同的数据形式,其中以电测法应用最为广泛。分析处理方法也不相同。 对模拟量可直接分析处理,也可将其转换为数字量后分析处理。前者设备较简单,后者精度和速度较高。 当模拟量是振动的时间历程,即用时域描述的振动量时,即环境预测,分析处理的主要内容就是进行各种频谱分析,以了解测试对象在频率域内的振动特性。即系统特征参数的测定,常用的频谱分析仪有恒定百分比带宽式、恒定带宽式、采用压缩时间历程的实时分析仪和具有并联滤波器的实时分析仪等。当模拟量是频率响应时,可归结为对机械阻抗数据的分析处理。测得的机械阻抗数据通常以幅频特性曲线和相频特性曲线、实部和虚部频率特性曲线或幅相频率特性曲线(Nyquist图)3种形式表达,统称为机械阻抗曲线。因此分析处理的主要内容是:根据机械阻抗曲线,通过模态分析,系统的振动特性也可以应用激光全息照相法拍下实物或模型在振动时的全息照片,识别测试对象在选定频率范围内的各阶模态参数和建立它的数字模型。 模拟量的数字分析处理是将测得的振动模拟量信号,经过模-数转换器变为相应的数字量,并与电子计算机相结合,然后输入数据处理机进行各种必要的分析。当测试结果直接以数字量表示时,则可利用软件在电子计算机上分析处理。20世纪80年代以来,4 毕业设计内容4.1 原理图的设计 测振传感器(拾振器)将机械振动量转换为与它成比例的电量。常用的测振传感器有发电型(如压电式、电动式和磁电式等)和电参数变化型(如电感式、电容式、电阻式和涡流式等)两类。不同类型的传感器需要配接不同类型的中间测量变换装置中间测量变换装置对传感器输出的电信号进行前置变换(电阻抗变换)、微积分运算、放大、调制和解调等,以便驱动后接的分析或显示、记录设备。分析设备完成对信号的频率分析。显示、记录设备给出振动信号(经过分析的或未经过分析的)的波形,并用数字或模拟方式指示出测量结果,以便于储存、分析信号和进行数据处理。4.2 机械结构设计4.3 控制设计硬件实现硬件实现的关键是要设计一个能够快速跟踪信号频率变化的倍频电路。利用锁相环的输出信号动态跟踪输入信号频率的这一特点,将VCO的输出信号经一分频器后送入相位比较器,则锁定时VCO输出信号的频率将等于K倍的输入信号的频率,即Kfi,适当地选择分频系数K就可以得到所需要的倍频信号的输出。本文所选的锁相环芯片为CD4046,使用中无需解调信号输出。比较器PC2用于比较两输入波形前沿或后沿相位的上跳或下跳(本电路采用上跳边沿触发方式)。如图所示,该电路主要由输入信号、隔直电容、平移、整形、锁相环组成。首先取键相位信号CIR作为测频输入,通过10F电容滤去太高的直流电压,送入MAX903比较放大器调节平移,输出的脉冲由74HCT14进一步整形后作为锁相环的输入信号;由锁相环电路实现倍频。为了实现FFT分析,在一个周期内采集2 n个数据,本设计中n取为6,也即是在一个周期内采集64个数据。分频器的输出信号CLK_64由CPLD产生。锁相环的输出信号CLK64实现AD的采样控制。本电路可以保证在有键相位的情况下,信号采集的第一起点在以键相位为基准,旋转36064点的位置上(因系统每周期采样64点)。软件实现图是本系统通过CPLD实现整周期采样的波形仿真图。实验中AD选用MAX1292,FIF0选用IDT7202。根据MAX1292的采样时序和FIFO的写时序,在CPLD中采用MOORE有限状态机来实现对时序的控制。RESET为系统的复位信号,CLK64为转速信号的64倍频;GLOBAL_CLK为CPLD的时钟信号;FF为FIFO的满标志位;AD_DATA为写入AD中的控制字节;AD INT、AD_WR、AD_RD、AD_HBEN、AD_CS为AD的相关控制信号:AD_START为启动AD采集信号,由MCU发出。AD_CHANNEL为AD的通道选择控制字,硬件电路上用两位开关来实现单通道、双通道及四通道的采集。图中为四通道采集。EF、HF、FIFOINT为FIF0单元控制模块。CLK 64、CLK 1000为分频器模块,用来实现整周期采样和自适应抗混叠滤波。实现AD整周期采样的过程描述如下:首先,CPLD在接收到CLK64信号的上升沿后,开始进行转换。MAX1292按照时序依次将四个通道的模拟量转换为数字量。其中,在每个通道的数据转换过程中,都要判断AD INT信号是否为低,若为低则表示转换完成,可以进行读操作。此时,再去判断FIFO的FF标志,若为1则表明FIFO未满,可以将数据写入FIFO中。当四个通道都完成如上操作后,再继续判断CLK64信号,循环执行,完成在一个转速信号周期内对各路信号采集64次。本装置中,被测的最大转速信号为1kHz,GLOBAL_CLK的信号频率为7MHz。毕业设计进度二月份:开题报告三月份:英文翻译1万字(10页)1-2个星期 完成绪论,第一章 拿出一个方案四月份:(理论设计) 第二章 方案比较 具体设计 建立数学模型五月份:第三章 理论 实验仿真 设计图纸六月初:中英文摘要 打印图纸,整理论文参考文献1 黄文虎,等设备故障诊断原理、技术及应用M北京:科学出版社,1994 2 赵俊超集成电路设计VHDL教程M北京:希望电子出版社,20023 李方泽,刘馥清,王正工程振动测试与分析M北京:高等教育出版社,19924 胡时岳、朱继梅:机械振动与冲击测试技术,科学出版社,北京,19835 闻邦椿 顾家柳 等.高等转子动力学理论、技术与应用M.北京:机械工业出版社,2000.1-76.6 陈德新 徐岚.灰度相关法在流动可视化图像解析中的应用J.华北水利水电学院学报,2000,21(1):47-50,55.7 万军,蒋世祥,蔡智勇; 旋转机械振动信号的小波包分解及故障检测 J;汽轮机技术; 2002年02期8 江志农,李艳妮; 旋转机械轴心轨迹特征提取技术研究 J;振动、测试与诊断; 2007年02期; 7-10+779 郭月强; 振动信号的测试与分析及其软件系统的开发 D;北京工业大学; 2002年10 张令弥. 振动测试与动态分析. 航空工业出版社,1992. 11 黄长艺, 严普强. 机械工程测试技术基础, 机械工业出版社, 1999. 12 姚天任, 江太辉. 数字信号处理. 华中科技大学出版社,2000. 目 录1 概述2 系统方案设计 2.1振动测试理论 2.2 振动测试技术 2.3 测试问题的提出 2.4 问题的解决和研究内容 2.5 方案1的设计 2.6 方案2的设计3 振动测试系统设计 3.1 机械系统设计 3.2 实验装置的选型 3.3 振动测试的激励4 实验结果及其分析 4.1 实验操作步骤 4.2 实验结果参考文献振动系统动态测试虚拟仪器的研究摘要:通过把振动系统动态分析理论与数字信号处理方法相结合,提出并系统地论述了振动系统动态测试虚拟仪器的开发原理。1 引言 随着数字计算机仿真技术的发展,虚拟仪器技术的应用越来越广泛,所谓虚拟仪器就是在通用计算机上增加一组软件和相关硬件,使使用者在操纵这台计算机时,就像是在操作一台他自己设计的专用传统电子仪器。在虚拟仪器系统中,硬件仅仅是为了解决信号的输入输出,而软件才是整个仪器系统的关键。任何一个使用者都可以通过修改软件的方法,很方便地调整仪器系统的功能与规模。虚拟仪器技术给用户一个充分发挥自己才能及想象力的空间,用户(而不是仪器厂家)可以根据自己的需求随心所欲地设计自己的仪器系统。因此有软件就是仪器的说法。 本文正是根据数字计算机仿真技术和虚拟仪器技术的思想,从简单实用的角度出发,把振动系统的动态测试理论和数字信号处理方法相结合,提出并系统地论述了振动系统动态测试虚拟仪器的开发原理。2 动态测试的基本概念 工程中大多数实际系统均用线性系统来描述,而线性系统的输入输出之间存在着简单的因果关系。因此,可以通过对被测系统输入输出物理量的测量和分析来确定系统的动态特性,这就是动态测试。动态特性的数学模型有多种形式:时域中常用的有微分方程、差分方程和状态方程;复数域中有传递函数、结构图;频域中有频率响应函数等。对机械振动系统进行动态特性辨识的主要目的是得到系统的频域特性或传递函数,并进一步获得机械振动力学系统的模态参数。频率响应函数用频率响应数据和曲线来表示,可以把频率响应函数看成是描述线性系统的非参数模型。对于机械振动系统,其动态特性常用“固有频率、阻尼比和振型”等所谓的模态参数来描述。在频率响应函数测试的基础上,可以通过参数识别的方法,即试验模态分析来建立机械振动系统的模态参数。 动态测试技术的主要内容就是对被测振动系统进行激励,通过振动测试、数据采集和信号分析,由输入和输出确定机械振动系统的动态特性。可见,动态测试分为两大部分:频率响应函数的测绘和模态参数的识别。机械振动系统动态特性测试与分析的重点是,如何利用测试分析所得的频率响应函数曲线,通过模态识别的方法,求解出表征机械振动系统动态特性的模态参数。模态识别可分为单模态识别和多模态识别,单模态识别是建立在小阻尼、弱耦合假设的基础上。对于模态密集的复杂结构或大阻尼情况,则需采用多模态分析。随着数字计算机技术的深入发展,现代模态参数识别大都采用多模态识别。3 振动系统动态测试虚拟仪器的设计3.1 虚拟仪器的提出和组成原理图虚拟仪器的优点在于能够共享计算机辅助测试系统的硬件资源和输入输出接口软件,用户只需设计有关的应用软件便可以实现相应的测试分析功能。这样不仅节约成本,而且还具有“实用性强,准确度高,效价比良,灵活性好,全自动化”的特点。基于数字计算机仿真技术和虚拟仪器技术的快速发展,通过把机械振动系统动态测试技术和数字信号处理方法相结合,提出了一种构建“振动系统动态测试虚拟仪器”的原理和方法。振动系统动态测试虚拟仪器的设计过程是:首先按流程图布置好振动系统动态测试虚拟仪器的硬件系统;然后,在个人计算机上开发相应的应用软件;最后,振动系统动态测试与模态识别的全过程,都是在个人计算机上通过软面板对话框和适当的控制按钮来完成。3.2 激振方式、激振器和振动传感器的布置激振方式可分为正弦稳态扫频和宽频带激振两大类。通常采用“正弦稳态扫频激振”,其优点是激振能量集中,信噪比高,对于线性系统的动态特性测试具有很高的测量精度。首先,根据测试分析要求及振动系统的机械结构特点,确定激振点的位置和测量点的数量N及各测量点的位置。对于空间结构系统,往往需要测量1个点的3个方向,这时应选用三轴振动传感器。安装激振器时,必须注意激振力的正确施加。在安装测振传感器时,应保证传感器能正确感受被测体的振动。小型振动系统必须选用微型传感器,以避免传感器附加质量影响振动系统动态特性。3.3 振动系统动态测试应用软件的总体研究 根据对前面振动系统动态测试虚拟仪器原理图的分析可知,振动系统动态测试应用软件的设计流程分为如下步骤: (1)计算机发出指令控制“信号发生器”产生激励信号,再通过“激振器”施加于被测振动系统。激励信号的数字序列为X(n)。 (2) 等待被测振动系统稳定一定时间,然后计算机再发出指令同步采集数据序列X(n)和Y1(n),Y2(n), ,YN(n) 。采样时间间隔或采样频率应符合采样定理及频率分辨率的要求。 (3)利用按本原理设计的专用软件分析处理X(n)和Y1(n),Y2(n),YN(n), 就可得出被测振动系统的模态参数和相应数学模型。显然,整套应用软件应包括两大部分,即操作控制部分和数据处理部分。由于目前计算机控制技术和虚拟仪器技术的发展,已足以解决此应用软件中的操作控制部分,且鉴于篇幅有限,所以笔者对应用软件的操作控制部分不再作进一步讨论,下面重点研究应用软件的数据处理部分。3.4 数据处理专用软件的设计原理3.4.1 滤波方法 由于数据序列在采集过程中难免会受到各种干扰, 所以必须采取有效的滤波算法对采集得到的数据序列Y1(n),Y2(n),YN(n)进行修正。通常应采用相关滤波算法,因为相关滤波可以有效地把真正由X(n)通过振动系统产生的信号从Y1(n),Y2(n),YN(n)中分离出来。将经过滤波算法处理后的数据序列记为YN1(n),YN2(n),YNN(n) ,n=1,2, ,D 。D为采样点数量。3.4.2 频率响应函数的确定 将确定性信号看成随机信号,采用离散数字随机信号的相关谱分析原理。 XK(k) = FFT X(n) YK(k) = FFT YN(n) SX(k) = XK(k)2 /D SY(k) = YK(k)2 /D SXY(k) = XK*(k)YK(k) /D XK*(K)为XK(k)的共轭复数。 为减少数据处理误差,要对SX(k),SY(k),SXY(k)采取平滑处理算法。分段时,相邻两段重叠50的效果最佳,即将记录段X(n),YN(n)分成多个样本,分别进行上述运算后再求平均值。VXY2(k) = SXY(k)2 / SX(k) SY(k), 相干系数VXY越接近1效果越好。H(k) =SXY(k) / SX(k), 由变量代换关系“k/(DTS) =f(TS为采样时间间隔)”可得到测试出的频率响应函数H(f) 。 对每一个数据序列Y1(n),Y2(n),YN(n)均用上述方法处理,就可以得到测试出来的N个频率响应函数H1(f),H2(f),HN(f) 。3.4.3 拟合求解模态参数机械振动系统模态参数识别,又称为曲线拟合,即采用最小二乘法将测试所得的频率响应值与系统模型值进行曲线拟合。优化准则就是使实测的频率响应值与理论数学模型对应值的总均方误差E极小。最小二乘曲线拟合过程中存在的有关问题,如方程组的病态问题,可根据数值分析理论,采用正交多项式的方法加以解决。通常机械振动系统的“理论数学模型”可用有理分式形式的传递函数来表示,即 H(s) = N(s)/D(s) = a0sm+a1sm-1+am-1s+am/ sn+b1sn-1+bn-1s+bn上式中,mn,n=2q,q为系统的自由度数。一般选定一个具体q值,令m=2q-1,然后循环拟合,直到满足E控制精度Eg。 对于每一个频率响应函数,设由D个采样点得到的D个测试频率响应值记为HC(f),对应的D个理论模型值记为HL(f),则总方差为E=(k=1,D)| HC(f)HL(f)|2。这是非线性优化问题,为使用方便和保证精度,可转化为线性优化问题。令s=jw=j2f,误差e=N(f) D(f)HC(f),则总方差转化为E=|N(f) D(f)HC(f)|2。此时,总方差函数为待求系数的线性函数,最小二乘解已变成线性优化问题。首先通过最小二乘拟合识别出有理分式的系数ai(i=1,2,m)和bi(i=1,2,n),然后由特征方程D(s)0解得极点sr及sr*(r=1,2,q),由留数公式ArH(s)(s-sr)| S=Sr”求出留数Ar及Ar* 。则传递函数可表示为H(s)r=1 , N Ar/(s-sr) + Ar*/(s-sr*)最后以sr= -r wr +j wr (1- r2)1/2,sr为共轭复数,求得各阶模态固有频率wr和阻尼比r 。对于线性系统,模态固有频率和阻尼比是系统的总体固有参数,理论上应该不随测量点而改变。依次对N个测量点的频率响应函数进行曲线拟合,可得到N组模态固有频率和阻尼比,再求出平均值作为系统的模态固有频率和阻尼比,进而,再求出各阶振型r 。3.4.4 其它输出振动系统动态测试所得的系统固有的各阶模态参数和用模态参数表示的各测量点,对应于输入点的数学模型。振动系统动态测试数据处理的全过程,可以用目前流行的程序设计语言VC+或Java2开发成专用的数据处理程序软件包,然后通过相关程序语言接口方法嵌入到用LabVIEW或LabWindows平台开发的振动系统动态测试虚拟仪器中。采用本虚拟仪器,可以很容易通过计算机多次重复振动系统动态测试的整个过程,进一步提高振动系统动态特性的测试精度。利用动态测试得到的数学模型,通过Matlab6.1仿真环境(或自己编制程序)可求出激振点处各种输入信号各测量点处的输出响应信号。4 结论经过实践检验,振动系统动态测试虚拟仪器在具体运用中具有“实用性强、准确度高、效价比良、重复性好及全自动化等优点。可以将数据采集功能和数据处理功能的软硬件集成在一起,生产出智能化振动系统动态测试虚拟仪器,也可以把数据处理软件单独形成专用的振动系统动态测试数据处理软件包,无论采取何种方式都具有广阔的应用前景。Virtual Instrument Research for Dynamic Testing of Vibration SystemAbstract: Through combining dynamic analyzing theory of vibration system with digital singles process method, developing principle of virtual instrument for dynamic test of vibration system is proposed and discussed.1. INTRODUCTIONWith the development of digital computer simulation technology, virtual instrument technology applications become wider, and the so-called virtual instrument is on the increase in the general computer a group of related hardware and software, so that users in the manipulation of this computer, as in the operation of a his design for the traditional electronic devices. In the virtual instrument system, the hardware simply to solve the input and output signals, and software is the key to the whole instrument. Any user can modify the software through the method, it is very easy to adjust the function of the system and apparatus scale. Virtual instrument technology to users a fully display their talent and imagination of space, users (rather than equipment manufacturers) can be arbitrary in accordance with their own needs and design their own instrument system. It is software is the apparatus argument. In this paper, in accordance with the digital computer simulation and virtual instrument technology ideas, from simple practical point of view, the vibration theory of the dynamic test and digital signal processing method, and proposed a systematic exposition vibration system dynamic testing Virtual Instrument development principle.2. The basic concept of dynamic testingMost of the actual system works were used to describe the linear system, the input-output linear system exists between the simple causal relationship. Therefore, it can be measured by the input-output system physical measurement and analysis to determine the dynamic characteristics, and this is dynamic testing. Dynamic characteristics of the mathematical model takes many forms: the time domain is commonly used differential equations, differential equations and equation of state; a complex domain transfer function, structure diagram; frequency domain, such as a frequency response function. Mechanical vibration system of the dynamic characteristics is the main purpose of identification by the system or frequency domain characteristics of the transfer function, and further mechanical vibration mechanical system modal parameters. Frequency response function with data and the frequency response curve, frequency response function can be described as a non-linear system parameter model. The mechanical vibration system, the dynamic characteristics of common natural frequency, damping and mode shape than the so-called modal parameters to describe. In the frequency response function test on the basis of parameters can be identified, that is, experimental modal analysis to establish the mechanical vibration system modal parameters. Dynamic testing technology is the main content of the measured vibration of the incentive system, through vibration testing, data acquisition and signal analysis, input and output determined by the mechanical vibration system dynamic characteristics. Thus, dynamic test is divided into two parts: the frequency response function mapping and identification of the modal parameters. Dynamic mechanical vibration system of testing and analysis of the key is how to use the proceeds of testing frequency response function curves through modal identification method for the characterization of mechanical vibration system dynamic characteristics of modal parameters. Modal Recognition state recognition can be divided into single-mode and multi-modal identification, identification of single-mode state is built on small damping loose on the basis of assumptions. The modal-intensive complex structure or damping, required a multi-modal analysis. With the in-depth development of the digital computer technology, the modern modal parameters to identify most of a multi-modal identification.3 Virtual Instrument Desire for Dynamic Testing of Vibration System3.1 Virtual Instrument made and composed schematicVirtual instrument has the advantage of being able to share computer-aided test the systems hardware resources and input-output interface software, users design of the application software will be able to achieve the corresponding test analysis functions. This will not only save costs, but also had a practical, accurate and high-titer, flexibility, and the fully automated features. Based on digital computer simulation and virtual instrument technology and the rapid development of the mechanical vibration system through dynamic testing technology, and digital signal processing method, the proposed construction of a dynamic vibration test system virtual machines the principles and methods. Dynamic vibration testing virtual instrument design process are as follows: First flow chart layout by good dynamic vibration test system hardware Virtual Instrument System; Then, in the corresponding development of personal computer software applications; Finally, the dynamic vibration testing and modal identification the entire process is in the personal computer through the soft panel dialog and appropriate control button.3.2 Vibration, vibration and vibration sensor layoutSinusoidal excitation mode can be divided into static and broadband excitation Swept two categories. Usually adopts a sine steady Sweep exciting, the merits of excitation energy is concentrated, high signal to noise ratio for linear dynamic characteristics of the test system with high precision. First, in accordance with requirements analysis and testing of mechanical vibration system structure, determine the location and vibration measurement points N and the number of the location of measurement points. The spatial structure system, often requires a measurement point three directions, then should use three-axis vibration sensors. Exciter installed, it should be noted excitation of the right to impose. Vibration sensors in the installation, should ensure that the sensor can accurately measured feelings of the vibration. Small vibration system to select micro-sensors, in order to avoid affecting the quality vibration sensors additional dynamic characteristics of the system.3.3 Dynamic vibration tests on the overall application softwareUnder the previous system dynamic vibration testing Virtual Instrument schematic analysis, we can see the dynamic vibration test application software design process is divided into the following steps: (1) directive issued computer control signal generator, produced encouraging signals, and through the Vibrator imposed on the measured vibration system. Encouraging signals in the digital sequence X (n). (2) the measured vibration system is stable for a certain period of time, and then issued instructions to the computer data acquisition synchronization sequence X (n) and Y1 (n), Y2 (n), YN (n). Sampling frequency or time interval should be consistent with the sampling theorem and frequency resolution requirements. (3) make use of the principle of design analysis software for X (n) and Y1 (n), Y2 (n), YN (n), measured vibration system can arrive at the modal parameters and the corresponding mathematical model . Obviously, application software package should include two parts, namely, operational control and data-processing part of some. As computer control technology and the development of virtual instrument technology has been sufficient to resolve this application in the operational control of light and space is limited, so the author on the application of software is no longer part of the operational control for further discussion below focuses on the application software data-processing part.3.4 Data processing software for the design principle3.4.1 filtering method The sequence data in the acquisition process will inevitably be all kinds of interference, it is necessary to take effective filtering algorithm to the data acquisition sequence Y1 (n), Y2 (n), YN (n) be amended. Usually associated filtering algorithm should be used, because the relevant filter can effectively from the real X (n) through the vibration signal from the system Y1 (n), Y2 (n), YN (n) separated. After filtering algorithm will be dealt with after the sequence data recorded as YN1 (n), YN2 (n), YNN (n), n = 1,2, D. D is the number of sampling points. 3.4.2 determine the frequency response function Signal will be uncertainty as random signal, a number of discrete random signals in the spectrum analysis of the relevant principles. XK (k) = FFT X (n) YK (k) = FFT YN (n) SX (k) = | XK (k) | 2 / D SY (k) = | YK (k) | 2 / D SXY (k) = XK * (k) YK (k) / D XK * (K) XK (k) Conjugate complex. In order to reduce data processing error, the SX (k), SY (k), SXY (k) to smoothing algorithm. Subparagraph, the two overlap adjacent 50 best results, the forthcoming record of the X (n), YN (n) is divided into a number of samples, respectively, to carry out such operations after seeking average. VXY2 (k) = | SXY (k) | 2 / SX (k) SY (k), the more coherent coefficient close to 1 VXY better results. H (k) = SXY (k) / SX (k), variable substitution k / (DTS) = f (TS for the sampling time interval) available to test the frequency response function H (f). For each data series Y1 (n), Y2 (n), YN (n) are dealing with the above method, it can be tested by the frequency response function of N H1 (f), H2 (f), HN (f).3.4.3 Fitting Solution modal parameters Mechanical vibration system modal parameter identification, also known as curve fitting, that is, using the least square method will test the frequency response of the system model and the value curve fitting. Optimization criterion is the value of the measured frequency response corresponds with the theoretical mathematical model of the total value of the minimum mean square error E. Least-squares fitting in the course of the relevant issues, such as the sick equations, numerical analysis according to theory, orthogonal polynomials method used to solve the problem. Usually mechanical vibration systems theory mathematical model can be rational fraction of the transfer function to indicate that H (s) = N (s) / D (s) = + a0sm a1sm-1 + + + am-1s am / sn-1 + + b1sn-1s + + bn bn - On, m n, n = 2q, q freedom of the system. Q general selected a specific value of the m = 2q-1, and then cycle fitting, until the E control accuracy Eg. For each frequency response function, a sampling from D-D tests by the frequency response of mind for HC (f), the corresponding value of D-theoretical model in mind for the HL (f), the total variance for E = (k = 1, D) | HC (f)-HL (f) | 2. This is the non-linear optimization problem, as easy to use and guarantee accuracy, and can be transformed to a linear optimization problem. Order s = = j2 f jw error e = N (f) - D (f) HC (f), the total variance into E = | N (f) - D (f) HC (f) | 2. At this time, the total variance function question for coefficient of linear functions, linear least squares solutions have become optimization problem. First identified by the least square fitting rational fraction of the coefficient ai (i = 1, 2, m) and the bi (i = 1, 2, n), and then by the characteristic equation D (s) = 0 Solutions in the pole sr and sr * (r = 1,2, q), by the residue formula Ar = H (s) (s-sr) | S = Sr sought stay of Ar and Ar *. will transfer function can be expressed as H (s) = r = 1, N Ar / (s-sr) + Ar * / (s-sr *) Finally, sr = - r wr + j wr (1 - r2) 1 / 2, sr for complex conjugate, and get the band Modal wr natural frequency and damping ratio r. For linear systems, modal natural frequency and damping ratio is inherent in the overall system parameters, in theory, should not change with the measurement point. N followed by the measurement point to the frequency response curve fitting function, can be inherent in group N modal frequencies and damping ratio, then calculated the average as a modal system natural frequency and damping ratio, in turn, obtained the order again r mode.3.4.4 Other Output System Dynamic vibration test systems inherent in the band modal parameters and the use of modal parameters of the measurement point, corresponding to the input point of the mathematical model. Dynamic vibration test data-processing system in the whole process can be used currently popular programming language VC + + or Java2 developed into a dedicated data processing package, and then through the relevant procedures language interface method using LabVIEW embedded development platform or LabWindows Dynamic vibration testing virtual instruments. Using the virtual instrument, can be easily repeated by the computer system dynamic vibration testing of the whole process, and further enhance the dynamic characteristics of vibration test system accuracy. The use of dynamic testing mathematical model through Matlab6.1 simulation environment (or their own development process) can be derived exciting point of the input signal measurement point response to the output signal. 4 Conclusion After the test of practice, vibration, dynamic testing system in the specific application of virtual instruments in a practical, accurate and high-titer, good repeatability, and the entire automation, etc. can function data acquisition and data processing hardware and software integration together to produce intelligent system dynamic vibration testing virtual instrument, can also form a separate data processing software for the system dynamic vibration test data-processing package, no matter in what ways have broad prospects.
收藏