CN103983412B - 用于振动有限元模型修正航空电子设备工作模态测量方法 - Google Patents
用于振动有限元模型修正航空电子设备工作模态测量方法 Download PDFInfo
- Publication number
- CN103983412B CN103983412B CN201410238286.7A CN201410238286A CN103983412B CN 103983412 B CN103983412 B CN 103983412B CN 201410238286 A CN201410238286 A CN 201410238286A CN 103983412 B CN103983412 B CN 103983412B
- Authority
- CN
- China
- Prior art keywords
- avionic device
- response
- test point
- processing device
- operation mode
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 230000004044 response Effects 0.000 claims abstract description 60
- 238000012360 testing method Methods 0.000 claims abstract description 60
- 238000012545 processing Methods 0.000 claims abstract description 26
- 230000003595 spectral effect Effects 0.000 claims abstract description 16
- 230000001133 acceleration Effects 0.000 claims abstract description 13
- 238000004458 analytical method Methods 0.000 claims description 11
- 238000004422 calculation algorithm Methods 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000004088 simulation Methods 0.000 claims description 7
- 238000013016 damping Methods 0.000 claims description 3
- 238000007689 inspection Methods 0.000 claims description 3
- 239000007767 bonding agent Substances 0.000 claims description 2
- 230000008569 process Effects 0.000 abstract description 4
- 238000002474 experimental method Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 230000005284 excitation Effects 0.000 description 3
- 238000012544 monitoring process Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000036541 health Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012731 temporal analysis Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 238000000700 time series analysis Methods 0.000 description 1
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
一种用于振动有限元模型修正的航空电子设备工作模态测量方法,步骤如下:一、布置响应测试点,贴装压电式加速度传感器:在其电路板上布置n个响应测试点,每个响应测试点处贴装压电式加速度传感器;二:安装航空电子设备和测试设备,搭建成工作模态测试平台;三:采集该航空电子设备在随机振动条件下的时域响应数据;四:利用模态辨识软件进行数据处理,获取航空电子设备的工作模态参数;通过以上步骤,构建了一个用于振动有限元模型修正的航空电子设备工作模态测量系统,完成从随机振动信号产生到时域响应数据采集,从时域响应数据样本到互谱函数生成,再到航空电子设备工作模态参数识别的全过程。该方法实现在线模态识别,节省测试时间和费用。
Description
技术领域
本发明提供一种用于振动有限元模型修正的航空电子设备工作模态测量方法,该方法主要研究的是在输入激励信息未知的情况下,仅通过振动试验输出的时域响应数据求解生成包含模态信息的互谱函数,来辨识航空电子设备的工作模态,用于修正可靠性仿真试验中的振动有限元模型,同时也可以作为航空电子设备振动设计改进的依据,属于基于故障物理的可靠性评估技术领域。
背景技术
模态分析是通过计算或试验来获取模态参数的过程,是一种研究结构动力特性的方法,是系统辨别方法在工程振动领域中的应用。其中,通过试验,利用采集的系统输入与输出信号获取模态参数,称为试验模态分析;仅利用采集的系统输出信号来进行模态参数辨别,则称为工作模态分析。可见,试验模态分析方法需要利用激励和响应的完整信息进行参数识别,但对于某些实际工程结构,很难甚至无法获得输入激励的完整信息。而工作模态分析仅需测试振动响应数据,且这些数据直接来源于结构实际所经受的振动工作环境,因而识别结果更符合实际情况和边界条件,这就解决了传统分析方法中输入不可测的问题。
工作模态分析方法主要包括随机减量法、多参考点复指数法、特征系统实现算法、峰值拾取法、时间序列分析法、NExT(自然激励技术法)等,目前已用于航天、桥梁、建筑、机械等领域。航空电子设备的振动仿真有限元计算主要包括模态计算和随机振动响应计算。由于建模过程中对复杂结构的简化,有限元模型与实际设备之间存在一定的差距,为了保证有限元模型的准确性,需要进行模态测试以及随机响应测试来修正有限元模型。目前工程上首先采用的锤击法测量模态,而后将航空电子设备固定在振动试验台上施加随机振动来获得响应数据。这种两步试验测试的方法需要多种设备,试验步骤繁琐,工作量较大。
本专利提出一种用于振动有限元模型修正的航空电子设备工作模态测量方法,可以在设备工作状态下,例如随机振动情况下,获得时域响应数据,利用多项式算法辨识模态并得到航空电子设备的随机响应。获得的模态和随机响应可用于振动有限元模型的修正,也可以用于振动设计改进以及振动状态的在线监控。
发明内容
1、目的:本发明的目的在于提供一种用于振动有限元模型修正的航空电子设备工作模态测量方法,该方法利用现有的软硬件设备搭建成工作模态测试平台,仅依赖振动试验时域响应数据,生成包含航空电子设备模态信息的互谱函数,利用多项式算法来辨识航空电子设备的工作模态,用于改进设计,校核航空电子设备可靠性仿真试验的振动有限元模型,修正有限元计算结果,为实现工作模态和振动响应的同时校核提供可能。
2、技术方案:本发明是通过以下技术方案实现的,首先安装航空电子设备和测试设备,搭建成工作模态测试平台,将压电式加速度传感器贴装在航空电子设备的印制电路板上,航空电子设备通过卡具固定在振动试验台上,通过线缆实现振动控制系统与振动试验台、压电式加速度传感器与数据采集设备的连接;接着进行试验,采集航空电子设备在振动激励下产生的时域响应数据,利用振动控制系统产生随机振动信号,经由功率放大器来驱动振动试验台,利用数据采集设备输出航空电子设备的时域响应数据;最后在模态辨识软件中进行数据处理,获取航空电子设备的工作模态参数。
本发明一种用于振动有限元模型修正的航空电子设备工作模态测量方法,其具体步骤如下:
步骤一:布置响应测试点,贴装压电式加速度传感器:在航空电子设备的电路板上布置n个响应测试点,该n个范围为5≤n≤30;每个响应测试点处贴装一个压电式加速度传感器;
步骤二:安装航空电子设备和测试设备,搭建成工作模态测试平台,其具体做法是:
a)将航空电子设备通过卡具固定到振动试验台上;
b)用线缆将振动控制系统连接到振动试验台上;
c)用线缆将压电式加速度传感器连接到数据采集设备上;
步骤三:采集该航空电子设备在随机振动条件下的时域响应数据:由振动控制系统产生随机振动信号,该随机振动信号的谱型与该航空电子设备可靠性仿真试验随机振动谱型一致,加载方向为航空电子设备印制电路板法线方向,利用数据采集设备输出航空电子设备的时域响应数据。
步骤四:利用模态辨识软件进行数据处理,获取航空电子设备的工作模态参数。其获取步骤如下:
a)利用模态辨识软件的绘图模块建立航空电子设备的印制电路板模型,在该模型上设置响应测试点,这里的响应测试点与试验设备上的响应测试点一致,并对这些响应测试点依次进行编号;
b)选择一个响应测试点为参考响应点,利用谱计算模块计算生成互谱函数;
c)采用峰值函数法确定模态阶数;选择多项式算法识别各阶模态和阻尼,并拟合计算留数,检查之后保存振型;通过模态置信准则(MAC)来判断识别结果的准确程度,最后保存以上计算结果。
通过以上四个步骤,构建了一个用于振动有限元模型修正的航空电子设备工作模态测量系统,该系统利用现有的软件和硬件设备,完成从随机振动信号产生到时域响应数据采集,从时域响应数据样本到互谱函数生成,再到航空电子设备工作模态参数识别的全过程。
其中,在步骤2中所述的“振动控制系统”,是选用现有的振动控制系统,如LDS、VENZO等。
其中,在步骤2中所述的“数据采集设备”,是选用现有的数据采集设备,如DataPhysics、BBM、SOMAT等。
其中,在步骤4中所述的“模态辨识软件”,是选用现有的模态辨识软件,如ME’scope、Vib’SYS等。
3、优点及功效:本发明一种用于振动有限元模型修正的航空电子设备工作模态测量方法,具有以下优点:
a.从随机振动响应试验中就可以获取航空电子设备的工作模态信息,避免了目前航空电子设备振动有限元模型修正时模态和随机振动响应分别进行试验测量的方式,从而节省了测试时间和费用,加快了航空电子设备的研制流程。
b.利用实时响应数据进行航空电子设备工作模态的辨识,使模态分析由传统的主要针对静止的设备扩展到用于现场运行的结构,不仅可以实现在线模态识别,还能更加准确地反映航空电子设备的实际动态特性,可用于航空电子设备的在线健康监测和损伤诊断。
附图说明
图1是本发明的方法流程图。
图2是航空电子设备工作模态测试装置示意图。
图3是航空电子设备印制电路板模型图。
图4是互谱曲线图。
图5是电源模块第一阶模态动画图。
图6是电源模块第二阶模态动画图。
图7是电源模块第三阶模态动画图。
图中符号、代号说明如下:
图3、5、6、7中的等表示响应测试点的编号。
具体实施方式
下面将结合附图和实例对本发明做进一步的详细说明。
以下实例针对的航空电子设备是某航空参数处理设备,该设备主要负责采集信息并进行信号处理,由箱体、电路板、母板,继电器盒和航插组成,其中主要有5种模块,分别为电源模块、数字信号模块、模拟信号模块、CPU模块以及功率转换模块。实例按照如图1所示的流程对该航空参数处理设备整机进行工作模态辨识。
本发明一种用于振动有限元模型修正的航空电子设备工作模态测量方法,如图1所示,具体步骤如下:
步骤一:在航空参数处理设备上布置响应测试点,贴装压电式加速度传感器。由于分析对象是结构较为简单的印制电路板,这里采用传统的均分布置法布置了16个响应测试点,并用硬性粘接剂在每个响应测试点处贴装一个型号为Endevco7703A-50的压电式加速度传感器。
步骤二:安装航空参数处理设备和测试设备,搭建成工作模态测试平台,如图2所示,其具体做法是:
a)将该航空参数处理设备通过卡具固定到振动试验台上;
b)用线缆将LDS振动控制系统连接到振动试验台上;
c)用线缆将压电式加速度传感器连接到Data Physics730数据采集设备上。
步骤三:加载随机振动信号,采集该航空参数处理设备在随机振动条件下的时域响应数据。由LDS振动控制系统产生随机振动信号,该随机振动信号的谱型与该航空参数处理设备可靠性仿真试验随机振动谱型一致,加载方向为该航空参数处理设备印制电路板法线方向。利用Data Physics730数据采集设备输出该航空参数处理设备的实时振动响应数据,格式为.uff。
步骤四:利用ME’scope模态辨识软件进行数据处理,获取该航空参数处理设备的工作模态参数。其获取步骤如下:
a)利用ME’scope模态辨识软件的绘图模块建立该航空参数处理设备的印制电路板模型,在该模型上设置响应测试点,这里的响应测试点与试验设备上的响应测试点一致,并对这些响应测试点依次进行编号,如图3所示。
b)选择一个响应测试点为参考响应点,利用谱计算模块计算生成互谱函数,其中,某响应测试点与其参考点的互谱曲线如图4所示。
c)利用互谱函数,结合多项式算法,确定该航空参数处理设备的工作模态参数,并进行检查和保存。采用峰值函数法中的幅值或虚部估计,确定模态阶数;选择多项式算法分频段识别各阶模态和阻尼;选择多项式算法拟合计算留数,检查之后保存振型,振型表见表1;通过模态置信准则(MAC)来判断识别结果的准确程度,如表2所示;最后检查和保存以上计算结果。还可根据需要,关联测试数据与结构模型,将各阶振型进行动画显示,其中,电源模块的前三阶振型动画截图见图5~图7。
表1航空参数处理设备工作模态振型表
表2模态置信准则(MAC)表
一阶振型 | 二阶振型 | 三阶振型 |
一阶振型 | 1 | 0.609 | 0.673 |
二阶振型 | 0.609 | 1 | |
三阶振型 | 0.673 | 0.494 | 1 |
本发明提出了一种用于振动有限元模型修正的航空电子设备工作模态测量方法,该方法利用现有的软硬件设备搭建成工作模态测试平台,从随机振动响应试验中就可以获取航空电子设备的工作模态信息,避免了目前航空电子设备振动有限元模型修正时模态和随机振动响应分别进行试验测量的方式,也无需对输入激励进行测试,不仅使结果更符合实际情况和边界条件,还节省了大量测试时间和费用。而且,测试获得的实时响应数据可以更加准确地反映航空电子设备的动态特性,从而用于航空电子设备的在线健康监测和损伤诊断。该方法为基于故障物理的可靠性仿真分析方法在工程中更好的推广应用提供了技术基础。
Claims (1)
1.一种用于振动有限元模型修正的航空电子设备工作模态测量方法,其特征在于:其具体步骤如下:
步骤一:在航空参数处理设备上布置响应测试点,贴装压电式加速度传感器;分析对象是印制电路板,采用传统的均分布置法布置了16个响应测试点,并用硬性粘接剂在每个响应测试点处贴装一个型号为Endevco 7703A-50的压电式加速度传感器;
步骤二:安装航空参数处理设备和测试设备,搭建成工作模态测试平台,其具体做法是:
a)将该航空参数处理设备通过卡具固定到振动试验台上;
b)用线缆将LDS振动控制系统连接到振动试验台上;
c)用线缆将压电式加速度传感器连接到Data Physics 730数据采集设备上;
步骤三:加载随机振动信号,采集该航空参数处理设备在随机振动条件下的时域响应数据;由LDS振动控制系统产生随机振动信号,该随机振动信号的谱型与该航空参数处理设备可靠性仿真试验随机振动谱型一致,加载方向为该航空参数处理设备印制电路板法线方向;利用Data Physics 730数据采集设备输出该航空参数处理设备的实时振动响应数据,格式为.uff;
步骤四:利用ME’scope模态辨识软件进行数据处理,获取该航空参数处理设备的工作模态参数;其获取步骤如下:
a)利用ME’scope模态辨识软件的绘图模块建立该航空参数处理设备的印制电路板模型,在该模型上设置响应测试点,响应测试点与试验设备上的响应测试点一致,并对这些响应测试点依次进行编号;
b)选择一个响应测试点为参考响应点,利用谱计算模块计算生成互谱函数,
c)利用互谱函数,结合多项式算法,确定该航空参数处理设备的工作模态参数,并进行检查和保存;采用峰值函数法中的幅值或虚部估计,确定模态阶数;选择多项式算法分频段识别各阶模态和阻尼;选择多项式算法拟合计算留数,检查之后保存振型;通过模态置信准则MAC来判断识别结果的准确程度,最后检查和保存以上计算结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410238286.7A CN103983412B (zh) | 2014-05-30 | 2014-05-30 | 用于振动有限元模型修正航空电子设备工作模态测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410238286.7A CN103983412B (zh) | 2014-05-30 | 2014-05-30 | 用于振动有限元模型修正航空电子设备工作模态测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103983412A CN103983412A (zh) | 2014-08-13 |
CN103983412B true CN103983412B (zh) | 2017-06-06 |
Family
ID=51275472
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410238286.7A Expired - Fee Related CN103983412B (zh) | 2014-05-30 | 2014-05-30 | 用于振动有限元模型修正航空电子设备工作模态测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103983412B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104807534B (zh) * | 2015-05-21 | 2018-04-13 | 华北电力大学(保定) | 基于在线振动数据的设备固有振动模式自学习识别方法 |
CN105547625A (zh) * | 2016-01-15 | 2016-05-04 | 华北电力科学研究院有限责任公司 | 一种瓷套模态测试方法、装置及系统 |
CN106482827B (zh) * | 2016-11-11 | 2018-12-25 | 北京航空航天大学 | 基于互功率谱函数模态参数识别的电子产品振动dlp方法 |
CN106525226B (zh) * | 2016-11-28 | 2020-07-10 | 株洲中车时代电气股份有限公司 | 一种基于现场振动载荷识别的评估方法及系统 |
CN107145620B (zh) * | 2017-03-14 | 2019-08-06 | 浙江大学 | 一种基于随机减量技术的结构动力特性识别方法 |
CN106959210B (zh) * | 2017-03-22 | 2019-06-11 | 国网江苏省电力公司电力科学研究院 | 一种用于敞开式隔离开关的分合状态检测方法及装置 |
CN107131984B (zh) * | 2017-05-31 | 2019-05-24 | 航天东方红卫星有限公司 | 一种随机振动试验测力方法 |
CN108573084B (zh) * | 2017-12-21 | 2022-01-04 | 蒙上阳 | 环境振动试验方法及系统 |
CN111896201A (zh) * | 2019-05-05 | 2020-11-06 | 天津市房屋质量安全鉴定检测中心有限公司 | 一种建筑振动模态测试方法 |
CN110163134B (zh) * | 2019-05-10 | 2022-03-29 | 大连理工大学 | 一种基于分频段加权最小二乘的结构损伤区域识别方法 |
CN111879400B (zh) * | 2020-07-01 | 2022-07-05 | 北京控制工程研究所 | 一种航天器机电产品的组件模态参数测量系统及方法 |
CN114254533B (zh) * | 2021-12-02 | 2022-11-01 | 中国兵器科学研究院宁波分院 | 考核疲劳振动对产品组部件固定角度影响和预测的方法 |
CN116484512B (zh) * | 2023-06-22 | 2023-09-01 | 西北工业大学 | 一种航空发动机盘鼓类转子预紧状态的辨识方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3242260B2 (ja) * | 1994-05-11 | 2001-12-25 | 株式会社日立製作所 | 構造物の振動試験装置、構造物の振動試験方法、並びに構造物 |
CN102012304B (zh) * | 2010-09-21 | 2012-06-27 | 上海大学 | 一种多频率合成的振动试验装置 |
JP2012233814A (ja) * | 2011-05-06 | 2012-11-29 | Microsignal Kk | 振動試験制御装置 |
CN103245474B (zh) * | 2013-05-21 | 2015-07-15 | 中国航空动力机械研究所 | 弹性支承器刚度测量装置及测量方法 |
-
2014
- 2014-05-30 CN CN201410238286.7A patent/CN103983412B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN103983412A (zh) | 2014-08-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103983412B (zh) | 用于振动有限元模型修正航空电子设备工作模态测量方法 | |
CN106525226B (zh) | 一种基于现场振动载荷识别的评估方法及系统 | |
CN105022860B (zh) | Pcb焊点随机振动疲劳寿命预测方法与系统 | |
CN104268335B (zh) | 微组装组件振动疲劳寿命预测方法和系统 | |
CN104792450B (zh) | 一种设备对安装基座激励力的间接测量方法 | |
CN104239646A (zh) | 微组装组件振动疲劳寿命预测仿真模型验证方法和系统 | |
CN107421763B (zh) | 一种设备故障检测方法及装置 | |
CN202033101U (zh) | 一种振动及噪声测量仪 | |
CN105138770B (zh) | 基于间接可靠性特征量的航天产品可靠性仿真评估方法 | |
CN103064010B (zh) | 基于希尔伯特-黄变换的模拟电路故障元件参数估计方法 | |
CN107590321A (zh) | 一种管道振动响应的整体测量方法 | |
CN107525849A (zh) | 一种基于光纤光栅的单输入单输出试验模态分析系统及方法 | |
CN102538950A (zh) | 一种发动机零部件固有频率的声学测试方法 | |
JP2018087799A (ja) | 地盤特性を自動校正する現地型地震早期警報システム及び関連方法 | |
CN104462862A (zh) | 一种基于三次b样条尺度函数的机械结构动载荷识别方法 | |
CN107036808A (zh) | 基于支持向量机概率估计的风电机组齿轮箱复合故障诊断方法 | |
CN104180823A (zh) | 一种温度补偿方法及装置 | |
CN102680214B (zh) | 波纹管隔振性能试验装置与试验方法 | |
CN206208920U (zh) | 一种振动冲击测试中等效部件传感模块 | |
CN102998367B (zh) | 一种基于虚拟派生结构的损伤识别方法 | |
CN202547906U (zh) | 波纹管隔振性能试验装置 | |
Abdullahi et al. | Accelerometer Based Structural Health Monitoring System on the Go: Developing Monitoring Systems with NI LabVIEW. | |
CN106969893B (zh) | 非接触式构件刚度检测设备及方法 | |
CN116147758A (zh) | 一种便携式实船机械噪声快速测试评估装置及评估方法 | |
CN116155143A (zh) | 马达参数的获取方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170606 |