CN110580391A - 一种柔性结构的基频模态测量方法 - Google Patents
一种柔性结构的基频模态测量方法 Download PDFInfo
- Publication number
- CN110580391A CN110580391A CN201910837491.8A CN201910837491A CN110580391A CN 110580391 A CN110580391 A CN 110580391A CN 201910837491 A CN201910837491 A CN 201910837491A CN 110580391 A CN110580391 A CN 110580391A
- Authority
- CN
- China
- Prior art keywords
- mode
- modal
- transient response
- subscript
- fundamental frequency
- 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.)
- Granted
Links
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种柔性结构的基频模态测量方法,通过有限元仿真建模分析得到结构的前三阶模态振型并按最大值进行归一化,然后根据模态叠加法,利用前三阶振型构造一个新的振型,再给结构的每个结点施加新振型时每个结点的位移改变量,通过瞬态响应分析获得结构的瞬态响应数据,并由瞬态响应数据通过特征系统实现算法辨识出结构的基频模态参数。本发明只需要结构的响应信号,能够有效地辨识出大型柔性结构的基频模态参数。
Description
技术领域
本发明涉及模态测量技术领域,尤其是一种柔性结构的基频模态测量方法。
背景技术
结构模态参数识别是动力学的核心内容之一,可以为动响应预示、故障诊断、安全性评估、结构优化设计等提供参考。模态参数识别理论体系已经较为成熟,广泛应用于航空、航天、土木、机械等工程领域。
由于大型柔性结构具有低频、模态密集的特征,在激振和模态辨识方面较为困难,所以针对大型柔性结构的基频模态测量成为亟待解决的技术问题。
发明内容
本发明所要解决的技术问题在于,提供一种柔性结构的基频模态测量方法,能够有效地辨识出大型柔性结构的基频模态。
为解决上述技术问题,本发明提供一种柔性结构的基频模态测量方法,包括如下步骤:
(1)通过有限元仿真软件建模分析计算得到柔性结构的频率f、振型S这两个模态参数,并记录前三阶模态振型ΔxJ i,再作归一化处理后对应结点的形变量其中ΔxJ i和中的下标i表示结点号,上标J表示模态阶次;
(2)构造新的叠加振型ΔXi,其中下标i表示结点号;
(3)给结构每个点xi施加由(2)得到的每个结点对应的新振型改变量ΔXi,其中下标i表示结点号;
(4)通过有限元仿真软件的瞬态响应分析计算初始位移下结构的瞬态响应函数Y(i);
(5)将(4)得到的瞬态响应数据Y(i),由特征系统实现算法辨识出结构的基频模态参数。
优选的,将步骤(1)中得到的前三阶原始模态振型[Δx1 i,Δx2 i,Δx3 i]按最大值进行归一化处理,即得到归一化模态振型其中坐标ΔxJ i和中的下标i表示结点号,上标J表示模态阶次。
优选的,将步骤(2)中得到的归一化模态振型构造新振型,即其中fi为第i阶对应的频率,中的下标i表示结点号,上标J表示模态阶次。
优选的,将步骤(3)中得到的新振型ΔXi施加到结构中对应的结点xi,其中下标i表示结点号;将步骤(3)中的施加初始振型后的结构通过瞬态响应分析,得到结构的每个结点xi对应的瞬态响应函数Y(i),其中下标i表示结点号。
优选的,将步骤(4)中得到的瞬态响应函数Y(i)利用特征系统实现算法辨识结构的模态参数,其中下标i表示结点号。
优选的,步骤(5)中,将步骤(4)得到的瞬态响应数据Y(i),由特征系统实现算法辨识出结构的基频模态参数具体为:
利用瞬态响应数据Y(i)构造Hankel矩阵,即:
对Hrs(0)进行奇异值分解Hrs(0)=UVWT;
由奇异值分解确定的阶次获得系统的最小实现:
对矩阵A进行特征值分解,求得特征值矩阵G,然后求出特征向量矩阵
根据矩阵A的特征值与系统特征值λr之间的关系确定模态频率ωr:
本发明的有益效果为:本发明通过给结构施加初始位移的方法,得到结构的瞬态响应函数Y(i),并由瞬态响应函数通过特征系统实现算法得到结构的基频模态参数;该方法能够有效的辨识出结构的基频模态参数,并且无需激励信号,为大型柔性结构的基频模态辨识提供方法。
附图说明
图1为本发明的方法流程示意图。
图2为本发明的悬臂梁结点位置示意图。
图3为本发明结点2的位移瞬态响应曲线示意图。
图4为本发明ERA算法振型与仿真振型MAC对比值示意图。
具体实施方式
如图1所示,一种柔性结构的基频模态测量方法,包括如下步骤:
(1)通过有限元仿真软件建模分析计算得到柔性结构的频率f、振型S这两个模态参数,并记录前三阶模态振型ΔxJ i,再作归一化处理后对应结点的形变量其中ΔxJ i和中的下标i表示结点号,上标J表示模态阶次;
(2)构造新的叠加振型ΔXi,其中下标i表示结点号;
(3)给结构每个点xi施加由(2)得到的每个结点对应的新振型改变量ΔXi,其中下标i表示结点号;
(4)通过有限元仿真软件的瞬态响应分析计算初始位移下结构的瞬态响应函数Y(i);
(5)将(4)得到的瞬态响应数据Y(i),由特征系统实现算法辨识出结构的基频模态参数。
(1)建模分析:通过相关参数,建立完整的结构有限元模型并进行模态分析;
表1结构材料参数表
长度/mm | 1200 |
弹性模量(MPa) | 2×10<sup>5</sup> |
密度(g/cm3) | 7.8 |
泊松比 | 0.3 |
通过仿真分析计算得出悬臂梁的前三节模态并按最大值作归一化处理。表2为前三阶模态各结点对应的振型改变量。
表2归一化振型表
频率/变形 | 结点2 | 结点3 | 结点4 | 结点5 | 结点6 | 结点7 |
0.841Hz | -0.045 | -0.165 | -0.338 | -0.545 | -0.770 | -1.000 |
5.116Hz | 0.231 | 0.611 | 0.758 | 0.482 | -0.168 | -1.000 |
13.950Hz | 0.570 | 0.897 | 0.135 | -0.731 | -0.416 | 1.000 |
(2)构造并施加初始位移:利用模态叠加法构造新的振型,并给悬臂梁施加新振型时每个结点对应的位移改变量,结构阻尼系数设为0.9,采样时间为0.005s,即采样频率为200Hz,共采样20000个步长,通过瞬态响应计算分析得到结点2~结点7的位移响应信号,如图2和图3所示。
(3)模态辨识;将瞬态响应的位移数据导出,利用特征系统实现算法,得到如表3所示的结果:
表3瞬态响应位移数据误差分析表
模态阶次 | ERA频率/Hz | 仿真频率/Hz | 频率误差/% |
1 | 0.852 | 0.841 | 1.31 |
2 | 5.322 | 5.116 | 4.03 |
3 | 14.575 | 13.950 | 4.48 |
如图4所示,分析可知,前三阶频率ERA识别结果与仿真结果吻合较好,1阶频率误差为1.31%,2阶频率误差为4.03%,3阶频率误差为4.48%。
根据模态置信度理论,利用本发明ERA算法得到的结果与仿真得到的振型结果的MAC值对比,从图中可以看出,ERA算法得到的试验模态振型与仿真得到的模态振型匹配度较高,模态置信度能达到0.95以上。
Claims (6)
1.一种柔性结构的基频模态测量方法,其特征在于,包括如下步骤:
(1)通过有限元仿真软件建模分析计算得到柔性结构的频率f、振型S这两个模态参数,并记录前三阶模态振型ΔxJ i,再作归一化处理后对应结点的形变量ΔxJ i’,其中ΔxJ i和ΔxJ i’中的下标i表示结点号,上标J表示模态阶次;
(2)构造新的叠加振型ΔXi,其中下标i表示结点号;
(3)给结构每个点xi施加由(2)得到的每个结点对应的新振型改变量ΔXi,其中下标i表示结点号;
(4)通过有限元仿真软件的瞬态响应分析计算初始位移下结构的瞬态响应函数Y(i);
(5)将(4)得到的瞬态响应数据Y(i),由特征系统实现算法辨识出结构的基频模态参数。
2.如权利要求1所述的柔性结构的基频模态测量方法,其特征在于,将步骤(1)中得到的前三阶原始模态振型[Δx1 i,Δx2 i,Δx3 i]按最大值进行归一化处理,即得到归一化模态振型[Δx1 i’,Δx2 i’,Δx3 i’],其中坐标ΔxJ i和ΔxJ i’中的下标i表示结点号,上标J表示模态阶次。
3.如权利要求1所述的柔性结构的基频模态测量方法,其特征在于,将步骤(2)中得到的归一化模态振型构造新振型,即其中fi为第i阶对应的频率,ΔxJ i’中的下标i表示结点号,上标J表示模态阶次。
4.如权利要求1所述的柔性结构的基频模态测量方法,其特征在于,将步骤(3)中得到的新振型ΔXi施加到结构中对应的结点xi,其中下标i表示结点号;将步骤(3)中的施加初始振型后的结构通过瞬态响应分析,得到结构的每个结点xi对应的瞬态响应函数Y(i),其中下标i表示结点号。
5.如权利要求1所述的柔性结构的基频模态测量方法,其特征在于,将步骤(4)中得到的瞬态响应函数Y(i)利用特征系统实现算法辨识结构的模态参数,其中下标i表示结点号。
6.如权利要求1所述的柔性结构的基频模态测量方法,其特征在于,步骤(5)中,将步骤(4)得到的瞬态响应数据Y(i),由特征系统实现算法辨识出结构的基频模态参数具体为:
利用瞬态响应数据Y(i)构造Hankel矩阵,即:
对Hrs(0)进行奇异值分解Hrs(0)=UVWT;
由奇异值分解确定的阶次获得系统的最小实现:
对矩阵A进行特征值分解,求得特征值矩阵G,然后求出特征向量矩阵
根据矩阵A的特征值与系统特征值λr之间的关系确定模态频率ωr:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910837491.8A CN110580391B (zh) | 2019-09-05 | 2019-09-05 | 一种柔性结构的基频模态测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910837491.8A CN110580391B (zh) | 2019-09-05 | 2019-09-05 | 一种柔性结构的基频模态测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110580391A true CN110580391A (zh) | 2019-12-17 |
CN110580391B CN110580391B (zh) | 2023-08-11 |
Family
ID=68812333
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910837491.8A Active CN110580391B (zh) | 2019-09-05 | 2019-09-05 | 一种柔性结构的基频模态测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110580391B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113406385A (zh) * | 2021-06-17 | 2021-09-17 | 哈尔滨工业大学 | 一种基于时域空间的周期信号基频确定方法 |
CN113935172A (zh) * | 2021-10-18 | 2022-01-14 | 西南交通大学 | 一种端接频变负载的传输线系统瞬态响应分析方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100031972A1 (en) * | 2008-05-28 | 2010-02-11 | Fbs, Inc. | Ultrasonic vibration system and method for removing/avoiding unwanted build-up on structures |
CN106484984A (zh) * | 2016-09-29 | 2017-03-08 | 北京大学 | 星载柔性附件热致微振动响应仿真分析平台 |
-
2019
- 2019-09-05 CN CN201910837491.8A patent/CN110580391B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100031972A1 (en) * | 2008-05-28 | 2010-02-11 | Fbs, Inc. | Ultrasonic vibration system and method for removing/avoiding unwanted build-up on structures |
CN106484984A (zh) * | 2016-09-29 | 2017-03-08 | 北京大学 | 星载柔性附件热致微振动响应仿真分析平台 |
Non-Patent Citations (2)
Title |
---|
任帅: "矩形薄板在周期性荷载作用下的动力响应研究", 《硕士电子期刊》 * |
李五欠: "基于振型叠加法的柔性结构形态重构研究", 《硕士电子期刊》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113406385A (zh) * | 2021-06-17 | 2021-09-17 | 哈尔滨工业大学 | 一种基于时域空间的周期信号基频确定方法 |
CN113406385B (zh) * | 2021-06-17 | 2022-01-21 | 哈尔滨工业大学 | 一种基于时域空间的周期信号基频确定方法 |
CN113935172A (zh) * | 2021-10-18 | 2022-01-14 | 西南交通大学 | 一种端接频变负载的传输线系统瞬态响应分析方法 |
CN113935172B (zh) * | 2021-10-18 | 2024-05-17 | 西南交通大学 | 一种端接频变负载的传输线系统瞬态响应分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110580391B (zh) | 2023-08-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2019011026A1 (zh) | 基于聚类分析的复合材料结构有限元模型修正方法 | |
CN108959794B (zh) | 一种基于深度学习的结构频响动力学模型修正方法 | |
CN110580391A (zh) | 一种柔性结构的基频模态测量方法 | |
CN108108559B (zh) | 一种基于子结构的结构响应获取方法及灵敏度获取方法 | |
CN109359333B (zh) | 一种包含多尺度形貌特征的体模型构建方法 | |
CN108460204B (zh) | 一种通过扬声器振动部件的受力和位移反推其材料动态力学参数的方法 | |
CN107908864A (zh) | 一种基于特征融合的复杂设备剩余使用寿命预测方法 | |
Diskin et al. | Grid convergence for turbulent flows | |
CN108021533A (zh) | 一种基于广义坐标系求解任意色散材料电磁特性的方法 | |
Liu et al. | Radial point interpolation based finite difference method for mechanics problems | |
CN103823406A (zh) | 一种基于模态质量分布矩阵的数控机床敏感环节辨识方法 | |
Abdollahzadeh et al. | Experimental and numerical investigation on large deformation reconstruction of thin laminated composite structures using inverse finite element method | |
Xu et al. | A ‘FE-Meshfree’TRIA3 element based on partition of unity for linear and geometry nonlinear analyses | |
Matsumoto et al. | Explicit evaluation of hypersingular boundary integral equation for 3-D Helmholtz equation discretized with constant triangular element | |
CN110083988B (zh) | 一种船舶水下辐射噪声评估方法 | |
CN104914167A (zh) | 基于序贯蒙特卡洛算法的声发射源定位方法 | |
CN109885896B (zh) | 一种基于复变差分灵敏度的非线性结构有限元模型修正方法 | |
CN111159934A (zh) | 一种铝基复合材料动力学模拟方法 | |
CN103984331B (zh) | 一种复杂产品装配质量的误差敏感性分析与控制方法 | |
CN113919078A (zh) | 一种尺寸效应下耦合应力梯度的缺口结构疲劳分析方法 | |
CN113704942A (zh) | 供水管网未知节点压力估计方法、装置及计算机系统 | |
Huang et al. | Residual stress identification in thin plates based on modal data and sensitivity analysis | |
CN110795820B (zh) | 工程结构裂纹问题求解方法以及装置 | |
CN105404740B (zh) | 一种基于联合指标与有限信息的修正灵敏度结构损伤识别方法 | |
Lee et al. | Two scale meshfree method for the adaptivity of 3-D stress concentration problems |
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 |