模态空间—时域,频域与模态空间之间的联系..

模态空间—时域,频域与模态空间之间的联系..
模态空间—时域,频域与模态空间之间的联系..

模态空间—时域,频域与模态空间之间的联系

这是一个常见的问题。但这个问题会涉及到很多不同的方面,因此我们先从一个简单的示意图入手,并且尽量避免涉及到过多的数学知识。我们通过这幅图来讨论时域、频域、模态空间和物理空间所有这些不同的方面。示意图中要讨论的内容有很多,所以我们先将他们分开逐一讨论,最后再将其整合到一起。你应该还记得我们先前的讨论,你问我模态分析是怎么一回事(“模态分析的简要解释”),这对下面的讨论很有帮助。

首先,我们考虑一个简单的悬臂梁,并假定对梁的端部施加一个单位脉冲激励。梁端部的响应是系统所有阶模态响应的叠加显示(如黑色时域曲线所示);可以注意到结构在几个不同频率上都有响应。通过对时域信号进行傅立叶变换,可以将梁端部的时域响应变换到频域。这个过程伴随有相当多的数学知识,但它对于我们来说是一个常见的变换。通常我们将这个由时域信号变换过来的频域曲线称为频响函数,简称为FRF(如黑色频域曲线所示);注意这条曲线上的多个峰,其对应于系统的多阶固有频率。

在进一步讨论时域和频域之前,让我们先讨论一下图中左上部的物理模型。我们知道悬臂梁具有多个振动固有频率。在每个固有频率上,结构具有特定的变形形式,我们称之为模态振型。对于该悬臂梁,我们看到,第一阶弯曲模态如蓝色所示,第二阶弯曲模态如红色所示,第三阶弯曲模态如绿色所示。当然,还有其他更高阶模态没有显示出来,并且此处我们只讨论前三阶模态,但这也很容易适用到更高阶模态。

我们也可以利用图中右上部的解析集中质量模型或者有限元模型(黑色所示)来计算这个实物梁。此模型通常用方程组来表示,方程组中的各点或自由度(DOF)之间相互影响或者耦合。这意味着,如果你上拉模型的某一个自由度,其他自由度也会受到影响并产生位移。这种耦合意味着,如果想确定系统响应,方程组会非常复杂。当描述系统的方程数目越来越多时,方程组会更加复杂。我们通常用矩阵将所有的运动方程组合在一起,以描述系统响应,如下所示:

其中[M]、[C]、[K]分别为质量、阻尼和刚度矩阵,对应的还有加速度、速度和位移以及施加于系统上的激振力。通常质量矩阵为对角阵;阻尼矩阵和刚度矩阵为对称阵,具有非对角元素,以描述系统的不同方程或自由度之间的耦合程度。矩阵的大小依赖于我们用以描述系统的方程数。从数学上讲,我们求得所谓的特征值,并利用模态变换方程,将这些耦合的方程组转换为一组解耦的单自由度系统,在新坐标系统中,它由模态质量、模态阻尼和模态刚度等对角矩阵来描述,这个新坐标系统称为模态空间,

表示下:

因此我们可以看出,利用模态变换方程,从物理空间到模态空间的变换是一个将复杂的耦合物理方程组转换为一组简单的解耦的单自由度系统的过程。并且,在图中我们看到,这个解析模型可以分解为一组单自由度系统,其中描述第一阶模态的单自由度如蓝色所示,第二阶模态如红色所示,第三阶模态如绿色所示。模态空间允许我们方便地利用简单的单自由度系统来描述系统。

现在我们回过头来讨论时域和频域响应,如黑色曲线所示。我们知道,系统的响应可以由各阶模态对其的贡献得到。黑色所示的总响应是由各模型响应结果求和得到的,各模型如蓝色第一阶,红色第二阶,和绿色第三阶模态所示。不论我是在时域还是在频域来描述系统,这都是正确的。各域是等同的,只是数据呈现的视角不同而已。这与货币非常相似——当我从一个国家到另一个国家时,每个国家的货币看起来都不一样,但其实都是一回事。所以,我们可以看出,时域总响应是由各阶模态时域响应即蓝色第一阶、红色第二阶和绿色第三阶模态的时域响应的贡献而来。我们也可以看出,总频响是由各阶模态

频响即蓝色第一阶、红色第二阶和绿色第三阶模态的频响贡献而来。(此处只显示了频响函数的幅值部分;这个函数实际上是复数形式,需要用频响的幅值和相位或者实部和虚部两部分来正确表示)。

我们可以将解析模型分解为一组单自由度系统,从而能够确定每个单自由度系统的频响函数,如蓝色第一阶、红色第二阶、和绿色第三阶模态所示。也能够确定每个单自由度系统的时域响应,对于脉冲输入引起的单自由度系统的响应,可通过闭式解求得或者对每个单自由度系统的频响函数进行简单的傅立叶逆变换得到。也能够测量单位脉冲输入引起的梁端部响应,并且对系统每阶模态的响应进行滤波,这样我们可以观察系统每阶模态的响应,如蓝色第一阶、红色第二阶、和绿色第三阶模态示。(当然,这里我对很多理论进行了简化处理,使得我们能够理解这些概念。)

至此我们已经将示意图中的各个部分进行了单独分析,我认为我们应该已经很清楚地认识到时域、频域、模态空间和物理空间之间真的没有区别。每个域仅仅是为了呈现或观测数据的一种方便的手段。但是有时,在一个域中比另一个域更容易看清事物。例如,时域总响应不能清楚地表明有多少阶模态对梁的响应有贡献。然而频域中的总频响函数却能很清楚地表明有多少阶模态被激发,以及各阶模态的频率。所以我们经常从一个域变换到另外一个域,只不过是因为数据更易于解释。

虽然介绍的内容并不多,但我希望这个简图及解释能帮助你有更深入的理解。好好思考一下。如果您有关于模态分析的任何其它问题,欢迎垂询。

模态空间系列(一)

你能为我解释模态分析吗?

嗯…说来有点话长,但下面的解释人人都可理解

请我用简单的概念来解释模态分析,以便任何人都可以理解它,你不是第一个人。简言之,模态分析是一种方法,籍此,可以根据结构的频率、阻尼和振型等固有属性-其动态特性-来描述结构。这真够拗口的,那我们来解释这是什么意思。不钻技术牛角尖,我经常用一个简单平板的振动模态来解释模态分析。对于刚接触振动及模态分析的工程师们来讲,这种解释向来有益。

考虑一个自由支撑平板,施加常力于平板一角。我们通常从静态的意义上来看待一个力,它在平板内引起某种静态变形。但这里我要做的是施加一个按正弦方式变化的力,振荡频率固定的常力。我们将改变振荡频率,但不改变力的峰值-仅是力的振荡频率改变。另在平板一角安装一加速度计来测量激励引起的平板应。

如果现在测量平板响应,注意到当改变输入力的振荡频率时,响应幅值也发生变化。频率升高过程中,不同时刻点上,幅值有增也有减。这好像很奇怪,因为我们施加常力于系统,响应幅值却随输入力的振荡速率而变化。但这确确实实发生了—当施加的力的振荡速率越来越接近于系统固有频率(或共振频率)时,响应增大,当振荡速率为系统固有频率时,响应达到最大值。想想看,这真令人惊奇,因为我每时每刻都施加了相同幅值的力-仅是振荡速率改变而已!

这个时域数据提供了非常有用的信息。但是如果采集到时域数据,并利用快速傅立叶变换将它变换到频域,则可以求得所谓的频响函数。现在有几点要关注:系统共振频率处,这个函数上有峰值。输入激励的振荡速率等于峰值频率的位置,

观察到了时域最大响应。

现在如果将时域波形跟频响图形叠加在一起,会注意到时域波形达到最大值时的振荡频率与频响函数最大峰处的频率相一致。所以,既可以利用时域波形来确定幅值达到最大值处的频率,也可以用频响函数来确定固有频率何处发生。显然,用频响函数更容易求。

结构为何具有这些固有属性,你感到大为惊奇。对了,在这些固有频率处,变形形式也大为不同,依赖于激振力用哪个频率。

好了,我们来看一看,在每个固有频率处,结构上的变形形式是怎样的。在平板上均布45个加速度计,测量不同激振频率的平板响应幅值。如果在每个频率处驻留-每次一个固有频率-可以观察结构上的变形形式。图中显示了按某一阶系统

固有频率激励时,得到的变形形式。在第一阶固有频率驻留时,平板具有第一阶弯曲变形形式,如蓝色所示。在第二阶固有频率驻留时,平板具有第一阶扭转变形形式,如红色所示。在第三、四阶固有频率驻留时,第二阶弯曲和第二阶扭转变形形式如绿色和紫红色所示。这些变形形式称为结构的模态振型。(从纯粹数学角度讲,这不完全正确。但事实上,此处简单讨论起见,这些变形形式非常接近于模态振型)

你看,我

们设计的所有结构都具有这些固有频率和模态振型。从本质上讲,这些特性依赖于结构的质量和刚度,它决定了固有频率和模态振型何处存在。作为设计工程师,需要确定这些频率,并且需要知道当力激励结构时,它们是如何影响结构响应的。理解模态振型和受激结构如何振动,将有助于设计工程师设计出更优的结构。然而模态分析的内容很多,这只是一个非常简单的解释。

现在我们能够更好地理解模态分析是什么—它研究结构的固有特性。利用固有频率和模态振型(依赖于结构的质量和刚度分布)帮助设计噪声和振动方面应用的结构系统。我们利用模态分析来帮助设计所有类型的结构,包括汽车、飞机、航天器、计算机、网球拍、高尔夫球杆…等等。

我希望这个非常简明的介绍有助于解释模态分析是什么。我曾经利用上面的例子向我的妈妈解释模态分析,并且我认为她第一次真正了解我是从事什么工作的了。自此以后,听到妈妈给她的朋友们解释模态分析,使用了不少非常接近于模态分析的词,其中最接近的一次是她称之为面条分析(noodle analysis)…好了,

有点扯远了!

平板前四阶模态的振型动画:

ITD模态参数识别matlab修改版

%ITD法识别模态参数 clear clc close all hidden format long %% txt文件下输入 fni=input('ITD法模态参数识别-输入数据文件名:','s'); fid=fopen(fni,'r'); mn=fscanf(fid,'%d',1); %模态阶数 %定义输入实测数据类型 %ig=1时域数据如冲击响应、自由振动、互相关函数、随机减量法处理结果 %ig=2频域数据如频响函数实部和虚部数据 ig=fscanf(fid,'%f',1); %ig=1时,f为采样频率sf,ig=2时,f为频率间隔df f=fscanf(fid,'%f',1); fno=fscanf(fid,'%s',1); %输出数据文件名 b=fscanf(fid,'%f',[ig,inf]); %实测时域或频域数据 status=fclose(fid); %% clc; clear all; format long [FileName,PathName] = uigetfile('*.mat', 'Select the Mat-files of time signal'); %窗口读文件,并获取包含路径的文件名 if isequal(FileName,0) disp('User cancel the selection'); %如果取消选择则显示提示 return; else FULLFILE=fullfile(PathName,FileName); Signal_str= sprintf('User selected signal file: %s',FULLFILE); disp(Signal_str); Struct=load(FULLFILE); end c=fieldnames(Struct); %得到一个元胞数组,包含Struct中各个域名(倘若有多个的话) b=getfield(Struct,c{1}); %获取c{1}对应的域中的内容 b=b(3601:9600); %% %ig=1时域数据如冲击响应、自由振动、互相关函数、随机减量法处理结果 %ig=2频域数据如频响函数实部和虚部数据 ig=input('数据类型ig='); f=input('采样频率f=');%指定采样频率 mn=input('计算模态阶数mn=');%指定计算模态阶数

模态参数识别的单模态法,模态参数识别的导纳圆法

一.模态参数识别的单模态法 常见的单模态识别有三种方法:直接读数法(分量分析法)、最小二乘圆拟合法和差分法。 所谓单模态识别法,是指一次只识别一阶模态的模态参数,所用数据为该阶模态共振频率附近的频响函数值。待识别的这阶模态称为主导模态,余模态称为剩余模态,剩余模态的影响可以全部忽略或简化处理。 1. 直接读数法(分量分析法) 1)基本公式 所谓分量分析法就是讲频响函数分成实部分量和虚部分量来进行分析。 N 自由度结构系统结构,p 点激励l 点响应的实模态频响函数可表示如下: 2222222111 ()(1)(1)N r r lp r er r r r r g H j K g g ωωωω=??--=+??-+-+?? ∑ (1.1) 其中r er lr pr K K φφ= ,为第二阶等效刚度 /r r ωωω= g 2r r r ζω= ,为第r 阶模态结构阻尼比 当ω趋近于某阶模态的固有频率时,该模态起主导作用,称为主导模态或者主模 态。 在主模态附近,其他模态影响较小。若模态密度不是很大,各阶模态比较远离,其余模态的频响函数值在该模态附近很小,且曲线比较平坦,即几乎不随频率而变化,因此其余模态的影响可以用一个复常数来表示,第r 阶模态附近可用剩余模态表示成: 222222211 ()()(1)(1)R I r r lp C C er r r r r g H j H H K g g ωωωω??-= -++??-+-+?? (1.2) ()lp H ω的实部和虚部可分别表示如下: 222211 ()(1)R R r lp C er r r H H K g ωωω??-= +??-+?? (1.3) 2221 ()(1)I I r lp C er r r g H H K g ωω??-= +??-+?? (1.4)

HHT方法在结构模态参数时域识别中的应用

[收稿日期]2008209212  [作者简介]任春(19802),男,2002年大学毕业,硕士,工程师,现主要从事桥梁工程设计与结构工程方面的研究工作。 HHT 方法在结构模态参数时域识别中的应用 任 春 (林同炎?李国豪土建工程咨询有限公司,上海200092) 张继承 (长江大学城市建设学院,湖北荆州434023) 罗奇峰 (上海防灾救灾研究所,同济大学,上海200092) [摘要]简略介绍了Hilbert 2Huang 变换(H H T )中EMD 法和Hilbert 变换及其结合随机减量技术 (RD T )在结构模态参数识别中的应用,并且对22DOF 动力模型的非平稳响应的模态参数时域识别进行了 探讨。通过将采用RD T 方法处理得到的时程曲线与采用RD T +H H T 方法处理得到的时程曲线相比较, 阐释了H H T 与RD T 方法结合对结构固有模态参数时域识别具有相当精度的机理。最后,将识别结果与 传统的频域法识别作了比较,得出了一些有意义的结论。 [关键词]H H T ;随机减量;模态参数;非平稳响应 [中图分类号]TU31113[文献标识码]A [文章编号]167321409(2008)042N115204 结构模态参数识别的方法有频域识别和时域识别,目前一些较成熟的时域法有ITD 法、STD 法、复指数法、ERA 法和随机减量法(RD T )等[1]。时域法不要求给出结构的具体激振信息,但一般对结构反应数据都有特殊要求。比如随机减量技术(Random Decrement Technique ,RD T ),它要求结构反应时程是平稳的、且要识别的各阶模态的间隔较大[2]。Hilbert 2Huang 变换(简称H H T )是分析非线性和非平稳数据的强有力工具,国内外已有学者用H H T 和RD T 相结合的方法识别结构模态参数[3~7],在时域法模态参数识别方面取得了一定的成果。但研究局限于结构激振形式趋于理想或简单化,也未对H H T 结合RD T 方法时域识别结构固有模态参数的有效性进行一定深度的机理阐述。为此,笔者利用H H T 结合RD T 的模态参数识别方法,对简单有限元模型在强非线性、非平稳性地震波下的响应进行了比较精确的模态参数时域识别,进一步肯定了H H T 与RD T 相结合的结构模态参数时域识别对于非平稳响应的有效性,并且通过处理得到的自由衰减时程曲线的数值比较,从本质上解释了H H T 与RD T 结合的时域识别方法的识别结果具有相当精度的机理。 1 模态参数识别中的HHT 方法 Hilbert 2Huang 变换(简称H H T )是由美籍华人Norden E 1Huang 于1998提出的一种数据处理方法,该方法因能有效处理非线性、非平稳数据而见长[8]。 111 EMD 法(H uang 变换) EMD (Empirical Mode Decompo sition )分解法认为任何信号都由一些不同尺度的振动模态构成,这些振动模态既可以是线性的也可以是非线性的,并据此将信号分解为若干固有模态函数(Int rinsic Mode Function ,IM F )的和。每个IMF 根据信号自身相临极值点间的时间间隔来定义和区分,并通过一称为筛选(Shifting Process )的步骤来完成分解。文献[8,9]给出了EMD 分解方法的详细步骤。112 Hilbert 变换 Hilbert 变换是研究线性及非线性动力系统特性的一种实用方法,利用信号的Hilbert 变换可以得到瞬时幅值及瞬时频率变化曲线,从中识别出系统的动力特性参数。文献[10,11]对于受冲击荷载的SDOF 系统自由振动响应的无阻尼固有圆频率ω0、 无阻尼固有频率f 0和阻尼比ξ的识别原理进行了详尽地阐述。在模态识别过程中,Hilbert 变换应用于由RD T 处理得到的关于某个模态的自由衰减时程,结? 511?长江大学学报(自然科学版)  2008年12月第5卷第4期:理工Journal of Yangtze U niversity (N at Sci Edit) Dec 12008,Vol 15No 14:Sci &Eng

环境振动下模态参数识别方法综述.

环境振动下模态参数识别方法综述 摘要:模态分析是研究结构动力特性的一种近代方法,是系统识别方法在工程振动领域中的应用。环境振动是一种天然的激励方式,环境振动下结构模态参数识别就是直接利用自然环境激励,仅根据系统的响应进行模态参数识别的方法。与传统模态识别方法相比,具有显著的优点。本文主要是做了环境振动下模态识别方法的一个综述报告。 关键词:环境振动模态识别综述 Abstract: The modal analysis is the study of structural dynamic characteristics of a modern method that is vibration system identification methods in engineering applications in the field. Ambient vibration is a natural way of incentives, under ambient vibration modal parameter identification is the direct use of the natural environment, incentives, based only on the response of the system for modal parameter identification method. With the traditional modal identification methods, has significant advantages. This paper is a summary report of the environmental vibration modal identification method. Keywords: Ambient vibration ;modal parameters ;Review 随着我国交通运输事业的发展,各种形式的大、中型桥梁不断涌现,由于大型桥梁结构具有结构尺大、造型复杂、不易人工激励、容易受到环境影响、自振频率较低等特点,传统模态参数识别技术在应用上的局限性越来越突出。传统的振动试验采用重振动器或落锤激励桥梁,需要投入大量人力和试验设备,激励成本增高,难度大,而且对于桥梁这样的大型复杂结构,激励(输入)往往很难测得,也不适合长期监测的实验模态分析。 环境振动是指振幅很小的环境地面运动。系由天然的和(或)人为的原因所造成,例如风、海浪、交通干扰或机械振动等,受激结构的振幅较小,但响应涵盖频率丰富。系统或者结构的模态参数包括:模态频率、模态阻尼、模态振型等。模态参数识别是系统识别的一部分,通过模态参数的识别可以了解系统或结构的动力学特性,这些动力特性可以作为结构有限元模型修正、故障诊断、结构实时监测的评定标准和基础。环境振动下的模态参数识别就是利用自然环境激励,根据结构的动

模态参数识别方法的比较研究

模态参数识别方法的比较研究 发表时间:2017-09-07T14:07:39.937Z 来源:《防护工程》2017年第9期作者:安鹏强[导读] 本文将频域法、时域法和整体识别法识别模态参数的应用范围、存在的优缺点进行对比、分析和说明。 航天长征化学工程股份有限公司兰州分公司甘肃兰州 730050 摘要:本文将频域法、时域法和整体识别法识别模态参数的应用范围、存在的优缺点进行对比、分析和说明,对模态参数识别的研究方向具有指导意义。 关键词:模态参数识别;频域法;时域法;整体识别法 引言 多自由度线性振动系统的微分方程可以表达为[1]: [M]{x ?(t)}+[C]{x ?(t)}+[K]{x(t)}={f(t)} 通过将试验采集的系统输入与输出信号用于参数识别的方法中,进而对系统的模态质量、模态阻尼、模态刚度、模态固有频率及模态振型进行识别,这一过程称为结构的模态参数识别。本文将对模态参数识别的频域法、时域法及整体识别法三者的应用范围、存在的优缺点进行对比、分析和说明。 1频域法 模态参数识别的频域法是结合傅里叶变换理论[1]形成的,这种方法是从实测数据的频响函数曲线上对测试结构的模态参数进行估计。图解法[1]是最早的频域模态参数识别方法,随之,又陆续发展了导纳圆拟合法[2]、最小二乘迭代法[2]、有理式多项式法[2]等多种频域模态参数识别方法。 频域法的优点是直观、简便,噪声影响小,模态定阶问题易于解决。频域法识别模态参数的思路是首先借助实测频响函数曲线对模态参数进行粗略的估计,进而将初步观测的模态估计值作为一些频域识别法的最初输入值,通过反复的迭代获取最终的模态参数。频域识别方法对于实测频响函数的分布容易控制,其输人数据是主观人为的。频域中参数识别方法识别结果的精准度,取决于测试试验中获得的频响函数质量的好坏。判断实测频响函数的质量,就要看其曲线的光滑[2]和曲线的饱满程度[2],曲线越光滑越饱满的实测频响函数,用其进行参数识别时,识别精度越高。 2时域法 模态参数识别的时域法的研究与应用比频域法晚,时域法可以克服频域法的一些缺陷。时域模态参数识别的技术优点在于无需获得激励力即可进行参数的识别[3-7]。对于一些大型的工程结构如大坝、桥梁等,获取激励荷载不太容易,但容易测得他们在风、地脉动等环境激励下的响应数据,把这些响应数据用于时域中一些参数识别的方法上,即可对测试结构的模态参数进行识别。 时域法的优点不仅在于其无需激励设备、减少测试费用而且可以避免由信号截断而造成对识别精度的影响,并且可实现对大型工程结构的在线参数识别,真实地反映结构的动力特性。但是由于响应信号中含有大量的噪声,这会使得所识别的模态中含有虚假模态。目前,对于如何剔除噪声模态、优化识别过程中的一些参数问题、以及怎样更稳定、可靠地进行模态定阶等成为时域法研究中的重要课题。目前常用的判定模态真假的方法是稳定图方法[8],该方法的基本思想在于不同阶次的系统模型会对虚假模态的影响比较大,在稳定图中出现次数最多的模态可认为是系统的真实模态。 3整体识别法 结构模态参数识别的单输入单输出类型是针对单个响应点的数据进行相应的计算,从而得到该测点对应的模态频率、阻尼比和振型系数等动力参数,但是对于有多个测点的试验,若要用单输入单输出类型的识别方法对多自由度结构进行参数识别,则需要对各个测点单独计算来识别各个测点对应的模态参数,通过对各个测点分别计算处理,得到每一个测点数据所识别的模态参数,然后求取所有测点响应识别的算术平均值来作为整体结构最终的识别结果。理论上讲,用每个测点数据识别的结果应该是一样的,但实际测试实验中,因测试实验中测点布置位置的不同、测试中其他因素及识别方法上的不完善会使得各个测点的识别结果不同、识别精度不同及错误的识别结果等现象。因此,对于多测点的测试试验,用单输入单输出类型的识别方法进行参数识别不仅会因多次重复导致计算工作量复杂累赘而且识别结果的正确性及精度无法保证。 整体识别的方法避免了单输入单输出类型的一些不足之处。该方法通过将结构上的所有测点的实测数据同时进行识别计算,所识别得到的结果作为结构整体的模态参数,每阶模态的固有频率和阻尼比是唯一的,减小了随机误差,提高了识别进度,并且使得计算工作量大大减少。 4三种识别方法的比较分析 (1)频域内的模态参数识别方法方便、快捷,但在实际运用中人为的主观选择性对识别结果的影响较大; (2)基于环境激励的时域模态参数的识别方法具有测试试验的花费较少、测试相对安全,并且识别精度较高。因此,基于环境激励的时域模态参数的识别方法已成为科研工作者研究的热点问题。 (3)对于多测点的测试试验,用频域和时域的单输入单输出类型识别模态参数不仅会因多次重复导致计算工作量复杂累赘而且识别结果的正确性及精度无法保证。整体识别法将所有测点的数据同时进行处理计算,得到结构的整体识别结果。整体识别方法通过对所有测点数据同时进行识别计算,减小了随机误差,提高了识别进度,使得计算工作量大大减少。 (4)对比时域和频域识别方法对虚假模态的剔除,可以看出,频域中的剔除虚假模态主要依据模态频率在频幅曲线图上会出现峰值的原理,利用该峰值处的幅值角是否为0°或180°来剔除虚假模态;相对频域剔除虚假模态的方法来说,时域中的剔除虚假模态的方法有定量的精度判别指标。总体看来,时域识别方法无法判别是否已将系统的所有模态进行识别且对于阻尼比的确定还有待研究。参考文献 [1] 曹树谦,张德文,萧龙翔. 振动结构模态分析-理论、实验与应用[M]. 天津大学出版社,2001. [2] 王济,胡晓. Matlab在振动信号处理中的应用[M]. 水利水电出版社,2006.

ITD模态参数识别matlab修改版

%ITDxx识别模态参数 clear clc close all hidden format long %% txt文件下输入 fni=input('ITD法模态参数识别-输入数据文件名:','s'); fid=fopen(fni,'r'); mn=fscanf(fid,'%d',1);%模态阶数 %定义输入实测数据类型 %ig=1时域数据如冲击响应、自由振动、互相关函数、随机减量法处理结果%ig=2频域数据如频响函数实部和虚部数据 ig=fscanf(fid,'%f',1); %ig=1时,f为采样频率sf,ig=2时,f为频率间隔df f=fscanf(fid,'%f',1); fno=fscanf(fid,'%s',1);%输出数据文件名 b=fscanf(fid,'%f',[ig,inf]);%实测时域或频域数据 status=fclose(fid); %% clc; clear all;

format long [FileName,PathName] = uigetfile('*.mat', 'Select the Mat-files of time signal'); %窗口读文件,并获取包含路径的文件名 if isequal(FileName,0) disp('User cancel the selection');%如果取消选择则显示提示 return; else FULLFILE=fullfile(PathName,FileName); Signal_str= sprintf('User selected signal file:%s',FULLFILE); disp(Signal_str); Struct=load(FULLFILE); end c=fieldnames(Struct);%得到一个元胞数组,包含Struct中各个域名(倘若有多个的话) b=getfield(Struct,c{1}); %获取c{1}对应的域中的内容 b=b(3601:9600); %% %ig=1时域数据如冲击响应、自由振动、互相关函数、随机减量法处理结果 %ig=2频域数据如频响函数实部和虚部数据 ig=input('数据类型ig='); f=input('采样频率f=');%指定采样频率 mn=input('计算模态阶数mn=');%指定计算模态阶数

基于响应信号的结构模态参数提取方法

第36卷 第7期2008年 7月 华 中 科 技 大 学 学 报(自然科学版) J.Huazhong Univ.of Sci.&Tech.(Natural Science Edition )Vol.36No.7 J ul. 2008 收稿日期:2007209203. 作者简介:毛宽民(19642),男,副教授,E 2mail :kmmao4645@https://www.360docs.net/doc/4613651105.html,. 基金项目:国家重点基础研究发展计划资助项目(2005CB7244101);国家高技术研究发展计划资助项目 (2006AA04Z407). 基于响应信号的结构模态参数提取方法 毛宽民a ,b 李 斌a (华中科技大学a 机械科学与工程学院;b 数字制造与装备国家重点试验室,湖北武汉430074) 摘要:基于现有实验模态分析技术,提出了以一个响应信号作为参考信号,并且只利用响应信号提取结构模态参数的方法.以一个自由的钢梁为实验对象,通过与传统的用传递函数矩阵进行模态参数识别的实验模态分析法的识别结果比较,验证了所提出方法的有效性:固有频率识别精度和模态阻尼比的识别精度较高,误差分别不超过0.5%和18%;振型有一定的误差,但是总体趋势是一致的,能够反映结构的振动形态.该方法特别适合于用力锤或激振器无法激振的大型重型结构,如大型机床等设备,也适合于那些不宜用外力激振的设备,如高精密机床等. 关 键 词:模态识别;频率响应函数;振动;运行模态分析 中图分类号:T H113.1 文献标识码:A 文章编号:167124512(2008)0720077203 R esponse signals 2based structural modal parameter identif ication M ao Kuanmi n a ,b  L i B i n a (a College of Mechanical Science and Engineering ;b The State Key Laboratory of Digital Manufacturing Equipment and Technology ,Huazhong University of Science and Technology ,Wuhan 430074,China ) Abstract :On t he basis of present research o n experimental modal analysis technology ,a new met hod is p ut forward ,which uses one response signal as a reference signal and only use response signals to identify st ruct ural modal parameters.This met hod is especially applicable for big and heavy st ruct ure which can not be excited by hammer or exciter ,like big machine tool equip ment ;and it is also right for t he equip ment s like high 2precision machine tool ,which can not be excited by external forces.This paper does experiment on a free steel beam and testify t he effectiveness of t his met hod by comparing it wit h t he t raditional experimental modal analysis met hod ,which uses a t ransform f unction mat rix to identify modal parameters.The frequency identification precision is very high ,t he error is less t han 0.5%;modal damping ratio n identification p recision is very high ,t he error is less t han 18%;t he model shape is generally t he same wit h a certain difference but is able to reflect t he vibration state.K ey w ords :modal identification ;frequency response f unctio n ;vibration ;operation modal analysis 结构的动态性能主要是由结构的模态参数决定的,结构模态参数提取方法,主要是实验模态分析技术,已经发展得相当丰富[1].这些技术的基本思路是通过实验,在知道结构的激励和响应的情况下,通过频率响应函数(频域法)[2]或脉冲响应函数(时域法)[3]提取结构的模态参数.利用结构的振动响应信号提取结构的模态参数,已经得到 了实验模态分析领域研究人员的普遍关注,提出了许多相应的方法,包括ODS (运行变形形状)和OMA (运行模态分析)[4~8].但这些新方法无一例 外的是在假设结构激励为稳态白噪声激励条件下.显然实际情况并非如此. 本文在现有实验模态分析技术的研究基础上,提出了以一个响应信号作为参考信号,并且只

模态参数辨识的频域方法

模态参数辨识的频域方法 吕毅宁 目录 模态参数辨识的频域方法 (1) 单点输入单点输出(SISO) (1) 图解法............................................................................................................ 1 频域多参考点模态参数辨识(MIMO ) ............................................................ 2 频域模态测试和参数辨识的可控性和可观性. (5) 单点输入单点输出(SISO) 图解法 1) 峰值检测 半功率点 )(2 1 )()(21r j H j H j H ωωω= = (1) r r ωωωξ21 2-= (2) 2) 模态检测 () ir r jr r r r ij r jr ir r r r r jr ir r r ij Q A Q j j Q j H ψσψσσψψωσωψψω-= -= -= +-= ) ()( (3) 式中,r Q 是模态比例换算因子。 在上式中,() r ij A 是模态质量r m 和模态刚度r k 的函数,又由下面的关系 2r r r m k ω= (4) 联立即可求得模态质量和模态刚度。 3) 圆拟合法 固有频率

max ==ω ωωd ds r r (5) 振型 r er I ij g k H 1 -= (6) jr ir r er k k ??= (7) er k 是等效模态刚度,r r r k g η= 是等效结构阻尼。 ()r ij r I ij ir r r jr R g k )(2==-H ?? (8) 模态阻尼 r g ) 1(2tan 211 ωα-= (9) r g ) 1(2tan 222 -= ωα (10) 2 tan 2 tan 22 1 12ωωω+-= r r g (11) 模态刚度 由 r er r I ij g k H 1 )1(-= =ω (12) 可得 r r I ij er g H k )1(1 =-= ω (13) 模态质量 2 r r r k m ω= (14) 其他方法,如正交多项式曲线拟合法,非线性优化辨识方法。 频域多参考点模态参数辨识(MIMO ) 一个N 自由度粘性阻尼线性系统,对它施加P 个激励力,在N 个点上进行响应

模态分析与参数识别

模态分析方法在发动机曲轴上的应用研究 xx (xx大学 xxxxxxxx学院 , 山西太原 030051) 摘要:综述模态分析在研究结构动力特性中的应用,介绍模态分析的两大方法:数值模态分析与试验模态分析。并着重介绍目前的研究热点一一工作模态分析。通过发动机曲轴的模态分析这一具体的实例,综述了运行模态分析国内外研究现状,指出了其关键技术、存在问题以及研究发展方向。 关键词:模态分析数值模态试验模态工作模态 Abstract :Sums up methods of model analysis applied on the research of configuration dynamic;al characteristio. It introduces two methods of model analysis: numerical value model analysis and experimentation model analysis. Then it stresses the hotspot-working model analysis.Some key techniques, unsolved problems and research directions of OMA were also discussed. Key words:Model analysis Numerical value model analysis Experimentation model analysis Working model analysis 1、引言 1.1模态分析的基本概念 物体按照某一阶固有频率振动时,物体上各个点偏离平衡位置的位移是满足一定的比例关系的,可以用一个向量表示,这个就称之为模态。模态这个概念一般是在振动领域所用,你可以初步的理解为振动状态,我们都知道每个物体都具有自己的固有频率,在外力的激励作用下,物体会表现出不同的振动特性。 一阶模态是外力的激励频率与物体固有频率相等的时候出现的,此时物体的振动形态叫做一阶振型或主振型;二阶模态是外力的激励频率是物体固有频率的两倍时候出现,此时的振动外形叫做二阶振型,以依次类推。

各种模态分析方法总结与比较

各种模态分析方法总结与比较 一、模态分析 模态分析是计算或试验分析固有频率、阻尼比和模态振型这些模态参数的过程。 模态分析的理论经典定义:将线性定常系统振动微分方程组中的物理坐标变换为模态坐标,使方程组解耦,成为一组以模态坐标及模态参数描述的独立方程,以便求出系统的模态参数。坐标变换的变换矩阵为模态矩阵,其每列为模态振型。 模态分析是研究结构动力特性一种近代方法,是系统辨别方法在工程振动领域中的应用。模态是机械结构的固有振动特性,每一个模态具有特定的固有频率、阻尼比和模态振型。这些模态参数可以由计算或试验分析取得,这样一个计算或试验分析过程称为模态分析。这个分析过程如果是由有限元计算的方法取得的,则称为计算模记分析;如果通过试验将采集的系统输入与输出信号经过参数识别获得模态参数,称为试验模态分析。通常,模态分析都是指试验模态分析。振动模态是弹性结构的固有的、整体的特性。如果通过模态分析方法搞清楚了结构物在某一易受影响的频率围各阶主要模态的特性,就可能预言结构在此频段在外部或部各种振源作用下实际振动响应。因此,模态分析是结构动态设计及设备的故障诊断的重要方法。 模态分析最终目标是在识别出系统的模态参数,为结构系统的振动特性分析、振动故障诊断和预报以及结构动力特性的优化设计提供依据。 二、各模态分析方法的总结

(一)单自由度法 一般来说,一个系统的动态响应是它的若干阶模态振型的叠加。但是如果假定在给定的频带只有一个模态是重要的,那么该模态的参数可以单独确定。以这个假定为根据的模态参数识别方法叫做单自由度(SDOF)法n1。在给定的频带围,结构的动态特性的时域表达表示近似为: ()[]}{}{T R R t r Q e t h r ψψλ= 2-1 而频域表示则近似为: ()[]}}{ {()[]2ωλωψψωLR UR j Q j h r t r r r -+-= 2-2 单自由度系统是一种很快速的方法,几乎不需要什么计算时间和计算机存。 这种单自由度的假定只有当系统的各阶模态能够很好解耦时才是正确的。然而实际情况通常并不是这样的,所以就需要用包含若干模态的模型对测得的数据进行近似,同时识别这些参数的模态,就是所谓的多自由度(MDOF)法。 单自由度算法运算速度很快,几乎不需要什么计算和计算机存,因此在当前小型二通道或四通道傅立叶分析仪中,都把这种方法做成置选项。然而随着计算机的发展,存不断扩大,计算速度越来越快,在大多数实际应用中,单自由度方法已经让位给更加复杂的多自由度方法。 1、峰值检测 峰值检测是一种单自由度方法,它是频域中的模态模型为根据对系统极点进行局部估计(固有频率和阻尼)。峰值检测方法基于这样的事实:在固有频率附近,频响函数通过自己的极值,此时其实部为零(同相部分最

模态参数识别频域法

振动模态分析理论与应用 模态参数识别频域法 当系统阻尼为比例阻尼或小阻尼时,阻尼矩阵经模态坐标变换后可以对角化,模态参数为实数,频响函数可按实模态展开。若在p 点激励,在l 点测量,则频响函数可表示为对于粘性阻尼有 ∑ 1 2 ωω ξ2ωω1 )ω(N i i i i lp lp j D H =+= 对于结构阻尼有 ∑ 1 2ωω 1 )ω(N i i i lp lp jg D H =+= 以上两式即为实模态参数识别的基本公式 6.1 实模态识别图解法 6.1.1 共振法 这是一种经典的模态分析方法,其基本思想是:当激励频率在系统某阶固有频率r ω附近时, 该阶模态导纳便起主导作用,其余各阶模态导纳的影响可忽略不计。即 )ω(≈)ω(lpr lp H H 此时,整个系统等效于一个单自由度系统。利用幅频特性和相频特性,便可确定系统的模态参数(参看图6-1)。 在待测结构上选择l 个测试点,求其中某点P 对所有各点的位移导纳。点数l 一般应等于或大于拟选的模态数N (自由度数)。则p 点对任意点l 的位移导纳可作如下处理: 当激振频率在r 阶固有频率附近时有 () () 2 22 2∞ 1 2 ωωξ4ωω1≈ ωω ξ2ωω1 )ω(∑ ++==r r i r lp i i i i i lp lp j D j D H 因此,测得的幅频曲线)ω(lp H 的第r 个峰值位置(共振频率点),便可近似确定r 阶固有频率r ω。由r ω两侧半功率带宽,可以确定r 阶模态阻尼比)ω2/Δω(ξr r =。由r ω处位移

有 ()r r lp r lp D H ξ2)ω(= 所以 ()()r lp r r lp H D ωξ2= 由因为 ()r pr lr r lp k D φ φ= 故在令pr φ的值等于1(振型中各元素具有确定的比例,其绝对值可认为地指定,不妨取第r 阶振型第p 个元素pr φ的值等于1)时,由原点导纳曲线的峰值可得r 阶模态刚度为 ) ω(ξ21 r pp r r H k = 此外,当r ωω=时,l 个导纳的幅值分别为 r r pr r r p k H ξ2φφ|)ω(|11= r r pr r r p k H ξ2φφ|)ω(|22= r r pr lr r lp k H ξ2φφ|)ω(|= 写成矩阵形式 = lr r r r r pr r lp r p r p k H H H φφφξ2φ| )ω(|| )ω(||)ω(|2121 因此,第r 阶振型为 {}±±±==| )ω(||)ω(|| )ω(|φφ φφ2121r lp r p r p lr r r r H H H 为表示振型的几何形状,上试中各导纳幅值应考虑其相位,可用正负号表示同相或反相,对 于实模态,其振型向量的各分量都是实数,且只有大小和正负之差。因此,系统作固有振动时,各坐标点同时达到极值,同时通过平衡位置。用共振法确定模态参数,方法简单直观。但由于忽略了相邻模态的影响,识别出的模态精度不高,特别是识别振型和阻尼时,可能引起较大的误差。另外当各阶模态耦合较密时可能识别不出单个模态。因此这种方法一般只用于对模态的初步分析。 6.1.2分量分析法 分量分析法的思想是利用导纳的实频和虚频特性识别出系统的模态参数。其优点是能考虑其余模态的影响。

实验模态分析与参数识别报告

2015 年春季学期研究生课程考核 (读书报告、研究报告) 考核科目:实验模态分析 学生所在院(系): 学生所在学科: 学生姓名: 学号: 学生类别: 考核结果阅卷人

实验模态分析与参数识别报告 模态分析可分为实验模态分析与工作模态分析等。模态分析的最终目标是识别出系统的模态参数,为结构系统的振动分析、振动故障诊断和预报、结构动力特性的优化设计提供依据。模态是机械结构的固有振动特性,每一个模态具有特定的固有频率、阻尼比和模态振型。这些模态参数可以由计算或试验分析取得,这样一个计算或试验分析过程称为模态分析。 1、模态分析原理 模态分析的过程是将线性时不变系统振动微分方程组中的物理坐标变换为模态坐标,使方程组解耦,成为一组以模态坐标及模态参数描述的独立方程,坐标变换的变换矩阵为振型矩阵,其每列即为各阶振型。 []{}[]{}[]{}{}()M X C X K X F t ++= (1) 其中:[]M —质量矩阵,[]K —刚度矩阵,[]C —粘性阻尼矩阵,{}()F t —激励力的列阵。 振动模态是弹性结构固有的、整体的特性。如果通过模态分析方法搞清楚了结构物在某一易受影响的频率范围内,各阶主要模态的特性,就可能预知结构在此频段内,在外部或内部各种振源作用下实际振动响应,而且一旦通过模态分析知道模态参数并给予验证,就可以把这些参数用于设计过程,优化系统动态特性,或者研究把该结构连接到其他结构上时所产生的影响。 方程(1)经傅氏变换,可得频域内的振动方程: [][][]{}{}2()()()M j C K X F w w w w -++= (2) 对应于固有频率1ω的固有振型或模态向量以幅值最大点为参考点的表达式为:{}{}11max 1()()X X w w w =。它们亦即简谐自由振动的主振型,满足以下关系式: [][]{}2()0i K M w j -= (3) 此代数方程组的系数行列式等于零,即为特征方程式;[]M ,[]K 为实数对称矩阵,[] M 正定,[]K 为非负定,其特征值20ω和对应的特征向量为实数。 主振型矩阵[]{}{}{}1 2,,,,n j j j j 轾=臌为实模态矩阵。根据振型的正交性: [][][][]1T M M j j =,[][][][]1T K K j j =;系统阻尼为比例阻尼时, [][][][]1T C C j j =。

由传递函数转换成状态空间模型(1)

由传递函数转换成状态空间模型——方法多!!! SISO 线性定常系统 高阶微分方程化为状态空间表达式 SISO ()()()()()()m n u b u b u b y a y a y a y m m m n n n n ≥+++=++++--- 1102211 )(2 211110n n n n m m m a s a s a s b s b s b s G +++++++=--- 假设1+=m n 外部描述 ←—实现问题:有了内部结构—→模拟系统 内部描述 SISO ???+=+=du cx y bu Ax x 实现问题解决有多种方法,方法不同时结果不同。 一、 直接分解法 因为 1 0111 11()()()() ()()()() 1m m m m n n n n Y s Z s Z s Y s U s Z s U s Z s b s b s b s b s a s a s a ----?=? =?++++++++ ???++++=++++=----) ()()() ()()(11 11110s Z a s a s a s s U s Z b s b s b s b s Y n n n n m m m m 对上式取拉氏反变换,则 ???++++=++++=----z a z a z a z u z b z b z b z b y n n n n m m m m 1) 1(1)(1)1(1)(0 按下列规律选择状态变量,即设)1(21,,,-===n n z x z x z x ,于是有 ?????? ?+----===-u x a x a x a x x x x x n n n n 12113 221

模态参数辨识方法——综述

模态参数辨识方法综述 摘要:本文对模态分析和模态参数识别进行了综述,对当前识别方法的原理、识别精度及适用条件进行阐述和比较,提出环境激励下模态参数识别方法需解决的关键问题及模态分析在缺陷检测和结构优化中作用。 关键词:模态分析模态参数识别模态分析与缺陷检测结构工作模态 0引言 模态分析是将线性时不变系统振动微分方程组中的物理坐标变换为模态坐标,使方程组解耦,成为一组以模态坐标及模态参数描述的独立方程,坐标变换的变换矩阵为振型矩阵,其每列即为各阶振型。模态是机械结构的固有振动特性,每一个模态具有特定的固有频率、阻尼比和模态振型。这些模态参数可以由计算或试验分析取得,这样一个计算或试验分析过程称为模态分析。振动模态是弹性结构固有的、整体的特性。如果通过模态分析方法搞清楚了结构物在某一易受影响的频率范围内,各阶主要模态的特性,就可能预知结构在此频段内,在外部或内部各种振源作用下实际振动响应,而且一旦通过模态分析知道模态参数并给予验证,就可以把这些参数用于(重)设计过程,优化系统动态特性,或者研究把该结构连接到其他结构上时所产生的影响。模态分析的最终目标是识别出系统的模态参数,为结构系统的振动分析、振动故障诊断和预报、结构动力特性的优化设计提供依据。 解析模态分析可用有限元计算实现,而实验模态分析则是对结构进行可测可控的动力学激励,由激振力和响应的信号求得系统的频响函数矩阵,再在频域或转到时域采用多种识别方法求出模态参数,得到结构固有的动态特性,这些特性包括固有频率、振型和阻尼比等。有限元法是当前分析机械结构模态的主要方法,很多学者研究了单裂缝和多裂缝缺陷对不同结构动态特性的影响,但这些研究仅局限于出现缺陷结构的当前状态,考虑到缺陷在机械结构使用过程中的扩展,提出了模态分析与缺陷扩展理论相结合的方法分析缺陷的发展趋势,便于机械结构剩余寿命的评估,使已达到设计寿命的结构在失效前仍然发挥其功能,节约了经济成本。 一般模态识别方法是基于实验室条件下的频率响应函数进行的参数识别方法,它要求同时测得结构上的激励和响应信号。但是,在许多工程实际应用中,工作条件和实验室条件相差很大,对一些大型结构无法施加激励或施加激励费用很昂贵,因此要求识别结构在工作条件下的模态参数。工作模态参数识别方法与传统模态参数识别方法相比有如下特点:一、仅

实验八MATLAB状态空间分析

实验八 线性系统的状态空间分析 §8.1 用MATLAB 分析状态空间模型 1、状态空间模型的输入 线性定常系统状态空间模型 x Ax Bu y Cx Du =+=+ 将各系数矩阵按常规矩阵形式描述。 [][][]11 121120 10 1;;;n n n nn n n A a a a a a a B b b b C c c c D d ==== 在MATLAB 里,用函数SS()来建立状态空间模型 (,,,)sys ss A B C D = 例8.1 已知某系统微分方程 22d d 375d d y y y u t t ++= 求该系统的状态空间模型。 解:将上述微分方程写成状态空间形式 0173A ??=??--??,01B ??=???? []50C =,0D = 调用MATLAB 函数SS(),执行如下程序 % MATLAB Program example 6.1.m A=[0 1;-7 -3]; B=[0;1]; C=[5 0]; D=0; sys=ss(A,B,C,D) 运行后得到如下结果 a = x1 x2 x1 0 1

x2 -7 -3

b = u1 x1 0 x2 1 c = x1 x2 y1 5 0 d = u1 y1 0 Continuous-time model. 2、状态空间模型与传递函数模型转换 状态空间模型用sys 表示,传递函数模型用G 表示。 G=tf(sys) sys=ss(G) 状态空间表达式向传递函数形式的转换 G=tf(sys) Or [num,den]=ss2tf(A,B,C,D) 多项式模型参数 [num,den]=ss2tf(A,B,C,D,iu) [z,p,k]=ss2zp(A,B,C,D,iu) 零、极点模型参数 iu 用于指定变换所需的输入量,iu 默认为单输入情况。 传递函数向状态空间表达式形式的转换 sys=ss(G) or [A,B,C,D]=tf2ss(num,den) [A,B,C,D]=zp2ss(z,p,k) 例 8.2 11122211220.560.050.03 1.140.2500.1101001x x u x x u y x y x -??????????=+??????????-????????????????=??????? ????? 试用矩阵组[a ,b ,c ,d]表示系统,并求出传递函数。 % MATLAB Program example 6.2.m

相关文档
最新文档