CN108489602B - 一种基于模态测试技术的惯性参数识别方法 - Google Patents
一种基于模态测试技术的惯性参数识别方法 Download PDFInfo
- Publication number
- CN108489602B CN108489602B CN201810268963.8A CN201810268963A CN108489602B CN 108489602 B CN108489602 B CN 108489602B CN 201810268963 A CN201810268963 A CN 201810268963A CN 108489602 B CN108489602 B CN 108489602B
- Authority
- CN
- China
- Prior art keywords
- inertia
- mass
- parameter identification
- test
- matrix
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 55
- 238000012360 testing method Methods 0.000 title claims abstract description 46
- 238000005516 engineering process Methods 0.000 title abstract description 8
- 230000004044 response Effects 0.000 claims abstract description 16
- 239000011159 matrix material Substances 0.000 claims description 27
- 238000005316 response function Methods 0.000 claims description 13
- 230000009466 transformation Effects 0.000 claims description 11
- 230000005284 excitation Effects 0.000 claims description 10
- 230000001133 acceleration Effects 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000013095 identification testing Methods 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 4
- 230000003595 spectral effect Effects 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 238000012935 Averaging Methods 0.000 claims description 2
- 230000035945 sensitivity Effects 0.000 claims description 2
- 238000000691 measurement method Methods 0.000 claims 1
- 238000004458 analytical method Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000004913 activation Effects 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000000725 suspension Substances 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M1/00—Testing static or dynamic balance of machines or structures
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M1/00—Testing static or dynamic balance of machines or structures
- G01M1/10—Determining the moment of inertia
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M1/00—Testing static or dynamic balance of machines or structures
- G01M1/12—Static balancing; Determining position of centre of gravity
- G01M1/122—Determining position of centre of gravity
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明实施例公开了一种基于振动模态测试技术的结构惯性参数识别方法,设计振动测试与分析技术,能够在较高支撑刚度的边界条件和较大噪声的环境条件下对结构惯性参数进行快速而准确的识别。适用于工厂环境下较大噪声的测试工作。本发明的方法包括:基于结构频率响应特性的直接参数识别法和剩余惯量法的组合方法。
Description
技术领域
本发明涉及振动测试技术领域,特别是一种基于模态测试技术的惯性参数识别方法。
背景技术
质量、质心位置、转动惯量和惯性积等惯性参数是结构动态仿真、性能优化、有限元分析、悬置安装位置和角度选取等过程中所需的重要参数。基于振动测试技术的结构惯性参数识别方法逐渐成为新的研究热点,目前主要分为三种方法,包括直接参数识别法,剩余惯量法和模态模型法。剩余惯量法应用较多;直接参数识别法根据结构最高阶刚体模态和第一阶弹性模态之间的加速度频响函数谱线实部数据,不依赖其他已知量,直接求出结构质量、质心、转动惯量和惯性积等共10个参数,但该方法对噪声和弹性模态的干扰十分敏感,很难得到准确结果;模态模型法通过识别质量归一化刚体模态振型,依据质量矩阵关于刚体模态振型的正交性,求出10个惯性参数,但该方法需要一次锤击得到4-6阶刚体模态,这往往很难做到,且该方法惯性参数的识别精度取决于刚体模态振型的识别精度,而对振型的识别精度往往不如对频率的识别精度。此外,噪声和刚体振型的耦合也会对结果造成不利影响。剩余惯量法精度较好,对噪声和弹性模态干扰不太敏感,但需要提前已知质量参数,而这往往是待求的未知量,与测试的目的不符。此外,为了使测试点的分布反映结构几何轮廓,需要布置大量测试点,增大了测试工作量,达不到高效测试的目的。
发明内容
为了克服现有技术存在的问题,本发明实施例提供一种新的测试方法,能够满足在较高弹簧支撑刚度和含有弹性模态和噪声干扰的情况下,一次测试得到所有10个惯性参数结果,并保持较高精度的方法。
为达到上述目的,本发明的实施采用以下技术方案:
本发明提供一种用于识别结构惯性参数的测试方法,所述方法用于一种结构的惯性参数识别测试,所述测试由于支撑***具有较大刚度,并且测试环境包含较大噪声,所述方法包括:
针对所述结构第一阶弹性模态振型,找到第一阶弹性模态振型节点位置,在符合惯性参数识别测试选点原则的前提下,在节点或其附近位置布置测点;
针对所述的测点分布,进行锤击法测试,得到包含最高阶刚体模态和第一阶弹性模态的频响函数矩阵;
利用频响函数矩阵代入所述直接参数识别法计算非第一阶弹性模态主要振动方向的质量值;
利用所述质量值,使用所述剩余惯量法,求出剩余9个惯性参数。
采用上述技术方案与现有技术相比,其有益效果为:
本发明提供的一种基于振动测试技术的惯性参数识别方法,与现有的剩余惯量法相比,本实施例结合直接参数识别法,实现了在不依赖任何已知量的前提下,在具有较高支承刚度和较大噪声干扰的测试环境下对结构惯性参数准确识别的目的,并且减少了测试点的数量,大幅降低工作量,有效地避免了直接参数识别法和剩余惯量法单独使用所面临的缺陷。
附图说明
图1为本发明基于模态测试技术的惯性参数识别方法的流程图;
图2为某型列车转向架构架模态试验测点分布图;
图3为某型列车转向架构架第一阶弹性模态振型和节点位置图;
图4为某型列车转向架构架惯性参数识别测试测点分布图;
图5为某型列车转向架构架频响函数谱线及质量线选带情况。
具体实施方式
本发明实施例提供了一种适用于较高刚度弹簧支撑、含有较大噪声干扰环境下的结构惯性参数识别测试技术,能够在弹性模态和噪声影响的条件下,保持较高的识别精度,并且减少了测点数量,大大降低测试工作强度,具备一定的工程实践意义。
本发明方法用于实际工厂环境下的结构惯性参数识别,实际工厂环境下的结构惯性参数识别受到较大刚度的弹簧支撑***影响,导致刚体模态不够低,并且测试环境中存在较大噪声影响,为达到上述目的,本发明的实施采用如下步骤:
步骤一:校正力锤和传感器灵敏度;
步骤二:均匀布置测点,高频采样,使用锤击法测试,3-5次锤击平均,测试结构第一阶弹性模态振型及其节点位置;
步骤三:在符合惯性参数识别测试选点原则的前提下,在第一阶弹性模态节点处或其附近重新布置测点,低频采样,8-10次锤击平均,获取包含刚体模态和第一阶弹性模态的频响函数矩阵;
步骤四:将频响函数矩阵、测点坐标及输入输出方向等信息代入直接参数识别法中,选择一段相对平直的加速度频响函数谱线,计算质量参数,取非第一阶弹性模态主要振动方向的质量值作为质量识别值;
步骤五:将识别的质量值代入剩余惯量法中计算剩余9个参数。
步骤三中惯性参数识别测试选点原则是:1、响应测试点的分布要反映结构几何轮廓,均匀分布,不能都处于一条直线上,不能位于对称轴或转动节线上,响应点坐标转换矩阵条件数不宜超过20左右;2、激励自由度个数至少6个,每个方向至少1个,力线方向不能全部相同,且力线不能穿过质心,激振力坐标转换矩阵条件数不宜超过20左右。
质量数值包括X、Y、Z、三个方向的结果,并不是每个结果都可用,具体计算过程为:响应点坐标转换矩阵:
激振力坐标转换矩阵:
为防止某些参数因数值过小导致相对误差过大,可考虑将参考坐标系原点平移一段距离,以保证测试数据的准确性。
步骤四中直接参数识别法的计算原理是通过最小二乘法求解6×6阶质量矩阵(此处脚标中的m、n指m行、n列):
其中R为响应点坐标转换矩阵,W为激励点坐标转换矩阵,H为实测频响函数矩阵实部数据。质量矩阵M各元素分布:
步骤五剩余惯量法的计算原理是:
识别质心:
识别转动惯量和惯性积:
其中m为已知质量,Jxxo,Jyyo,Jzzo,Jxyo,Jxzo,Jyzo分别为三个转动惯量和三个惯性积,Mx,My,Mz为广义激励力矩。
本发明适用于具有较高弹簧刚度支撑结构在含较大噪声环境下进行的惯性参数快速测试,以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进,这些改进也应视为本发明的保护范围。
Claims (3)
1.一种基于振动模态测试的惯性参数识别方法,其特征在于,该方法针对实测结构的加速度频率响应函数矩阵,识别出第一阶弹性模态振型,找到模态节点位置;
利用所述弹性模态节点位置作为响应测试点,采用直接参数识别法求出测试点三个方向对应的质量参数,并选择第一阶弹性模态振动幅值较小的方向对应的质量数值作为识别结果;
基于所述质量数值通过剩余惯量法求解剩余9个惯性参数;
所述识别方法具体包括以下步骤:
步骤一:校正力锤和传感器灵敏度;
步骤二:均匀布置测点,高频采样,使用锤击法测试,3-5次锤击平均,测试结构第一阶弹性模态振型及其节点位置;
步骤三:在符合惯性参数识别测试选点原则的前提下,在第一阶弹性模态节点处或其附近重新布置测点,低频采样,8-10次锤击平均,获取包含刚体模态和第一阶弹性模态的频响函数矩阵;
步骤四:将频响函数矩阵、测点坐标及输入输出方向信息代入直接参数识别法中,选择平直的加速度频响函数谱线,计算质量参数,取非第一阶弹性模态主要振动方向的质量值作为质量识别值;
步骤五:将识别的质量值代入剩余惯量法中计算剩余9个参数;
步骤三中惯性参数识别测试选点原则是:
a、响应测试点的分布要反映结构几何轮廓,均匀分布,不都处于一条直线上,不位于对称轴或转动节线上,响应点坐标转换矩阵条件数不超过20;
b、激励自由度个数至少6个,每个方向至少1个,力线方向不全部相同,且力线不穿过质心,激振力坐标转换矩阵条件数不超过20;
所述质量数值包括X、Y、Z、三个方向的结果,具体计算过程为:
响应点坐标转换矩阵:
激振力坐标转换矩阵:
两端同时右乘f-1,并使用最小二乘法,得到:
其中R为响应点坐标转换矩阵,W为激励点坐标转换矩阵,H为实测频响函数矩阵实部数据;
所述质量矩阵M6×6表示为:
取弹性模态振幅较小的方向的m作为质量参数识别值。
2.根据权利要求1所述的一种基于振动模态测试的惯性参数识别方法,其特征在于,将所述质量参数识别值m作为剩余惯量法的必要已知量代入计算,得到质心,包括Xc,Yc,Zc,转动惯量包括Jxxo,Jyyo,Jzzo,和惯性积包括Jxyo,Jyzo,Jxzo。
3.根据权利要求1所述的一种基于振动模态测试的惯性参数识别方法,其特征在于,在设计响应点坐标转换矩阵和激振力坐标转换矩阵时,为防止某些参数因数值过小导致相对误差过大,将参考坐标系原点平移一段距离,以保证测试数据的准确性。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810268963.8A CN108489602B (zh) | 2018-03-29 | 2018-03-29 | 一种基于模态测试技术的惯性参数识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810268963.8A CN108489602B (zh) | 2018-03-29 | 2018-03-29 | 一种基于模态测试技术的惯性参数识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108489602A CN108489602A (zh) | 2018-09-04 |
CN108489602B true CN108489602B (zh) | 2020-07-07 |
Family
ID=63317309
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810268963.8A Active CN108489602B (zh) | 2018-03-29 | 2018-03-29 | 一种基于模态测试技术的惯性参数识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108489602B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111062143B (zh) * | 2019-12-30 | 2022-11-04 | 重庆长安汽车股份有限公司 | 一种汽车发动机惯性参数识别方法 |
CN112611511B (zh) * | 2020-11-30 | 2022-03-29 | 桂林电子科技大学 | 一种基于加速度频响函数的复杂构件惯性参数获取方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
BR8402546A (pt) * | 1984-05-28 | 1986-01-07 | Paulo Euzebio | Aperfeicoamento em trava inercial para encostos de bancos de veiculos |
JP3757839B2 (ja) * | 2001-10-10 | 2006-03-22 | 株式会社明電舎 | 供試体の慣性モーメント測定方法とその装置 |
JP2009156425A (ja) * | 2007-12-27 | 2009-07-16 | Aisin Aw Co Ltd | 自動変速機 |
CN100595084C (zh) * | 2008-05-28 | 2010-03-24 | 天津大学 | 大功率柴油机动力总成六点悬置***设计方法 |
CN104296924B (zh) * | 2014-11-04 | 2017-10-10 | 广西大学 | 动力总成惯性参数测试方法 |
-
2018
- 2018-03-29 CN CN201810268963.8A patent/CN108489602B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN108489602A (zh) | 2018-09-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106767917B (zh) | 一种斜置冗余惯导***标定误差模型建模方法 | |
CN109822574A (zh) | 一种工业机器人末端六维力传感器标定的方法 | |
CN102636140B (zh) | 测量空间六自由度运动的伸缩拉杆式并联装置 | |
CN107356387B (zh) | 一种模态试验中多传感器附加质量消除方法 | |
CN108168774B (zh) | 一种空间矢量力校准方法 | |
CN110375694B (zh) | 基于便携关节式坐标测量机圆光栅测角误差的自校准方法 | |
CN106468623B (zh) | 一种整车状态下动力总成悬置***刚体模态参数测试方法 | |
CN105241630A (zh) | 应用于激波风洞测力试验的脉冲型杆式应变天平 | |
CN108489602B (zh) | 一种基于模态测试技术的惯性参数识别方法 | |
CN102679925A (zh) | 机器人定位误差测量方法 | |
CN106767673B (zh) | 一种卫星高精度光学敏感载荷的指向测量方法 | |
CN104296924B (zh) | 动力总成惯性参数测试方法 | |
CN113111547A (zh) | 基于缩减基的频域有限元模型修正方法 | |
Xu et al. | Dynamic decoupling and compensating methods of multi-axis force sensors | |
CN110516350B (zh) | 一种基于各向异性加权的ers点误差修正方法 | |
CN104296897B (zh) | 基于星箭连接环应变测量的星箭六自由度界面力计算方法 | |
CN106871933B (zh) | 一种消除频响函数中多传感器附加质量的方法 | |
CN109342008B (zh) | 基于单应性矩阵的风洞试验模型迎角单相机视频测量方法 | |
Sun et al. | Design of a novel Six-axis force/torque sensor based on strain gauges by finite element method | |
CN109061723B (zh) | 一种隧洞岩爆孕育过程的微震源高精度定位方法及*** | |
Prato et al. | A reliable sampling method to reduce large sets of measurements: a case study on the calibration of digital 3-axis MEMS accelerometers | |
CN110221245B (zh) | 联合估计目标位置和非视距误差的鲁棒tdoa定位方法 | |
CN109443333B (zh) | 一种陀螺阵列反馈加权融合方法 | |
CN103177192A (zh) | 五孔探针角度特性曲线数据处理方法 | |
CN107563009B (zh) | 一种汽车转向***cae模态识别方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |