CN117490823A - 一种基于视觉测量的新型振动提取与振型观测方法 - Google Patents
一种基于视觉测量的新型振动提取与振型观测方法 Download PDFInfo
- Publication number
- CN117490823A CN117490823A CN202311446779.5A CN202311446779A CN117490823A CN 117490823 A CN117490823 A CN 117490823A CN 202311446779 A CN202311446779 A CN 202311446779A CN 117490823 A CN117490823 A CN 117490823A
- Authority
- CN
- China
- Prior art keywords
- vibration
- edge
- modal
- order
- time sequence
- 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
- 238000005259 measurement Methods 0.000 title claims abstract description 20
- 238000000605 extraction Methods 0.000 title claims abstract description 19
- 238000000034 method Methods 0.000 title claims abstract description 18
- 239000011159 matrix material Substances 0.000 claims abstract description 9
- 230000000007 visual effect Effects 0.000 claims abstract description 6
- 238000005070 sampling Methods 0.000 claims description 9
- 238000000354 decomposition reaction Methods 0.000 claims description 6
- 230000004931 aggregating effect Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 239000002131 composite material Substances 0.000 claims description 3
- 238000013016 damping Methods 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 abstract description 5
- 238000012360 testing method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H9/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明涉及结构振动分析技术领域,公开了种基于视觉测量的新型振动提取与振型观测方法,包括如下步骤:步骤一、搭建测量平台,利用高速相机获取被测结构的运动视频数据;步骤二、定位边缘区域,获取运动视频中结构边缘所有位置的时序像素切片,组成时序矩阵;步骤三、基于偏差提取算法,提取结构边缘时序矩阵上亚像素精度的振动信号。本发明提出的基于视觉测量的新型振动提取与振型观测方法,首先无需在被测结构上布置接触式传感器,降低了测量难度和测量成本;其次通过视觉测量振动信号,提高了空间分辨率;最后利用偏差提取算法,仅关注结构边缘的时序扰动,提升了结构振动分析与振型观测的效率。
Description
技术领域
本发明涉及结构振动分析技术领域,尤其涉及一种基于视觉测量的新型振动提取与振型观测方法。
背景技术
振动现象广泛存在于各种工程应用中,如机械、车辆和建筑领域。运行中的工程结构和设备会产生大量振动信号,这些信号中包含着结构内在特性和状况信息,对于反映系统状态及其变化规律至关重要。随着科学技术的进步,非接触测量在振动测量领域变得越来越重要,成为了一个研究的关键方向。
相对于接触式传感器而言,基于视觉的非接触式测量方法具有多项优势,包括高精度、普适性、多点测量能力以及长期稳定的工作性能,并且能够提供更高的空间分辨率,使得基于相机的设备成为振动信号全场测量的理想选择,这为实际应用和研究领域带来了重要的价值和意义。
传统的接触式传感器在实验模态试验中通常能提供高精度和可靠的数据。然而,这些传感器需要在被测结构表面进行安装,复杂耗时且容易产生负载效应。
基于视觉的设备具有灵活的结构和更高的空间分辨率,使其更适合进行远程安装和全场测试。对于工程应用而言,考虑全场的振动提取过程相当耗时,结构振型观测流程复杂繁琐,难以实现快速、高精度振动信号提取与振型观测。
发明内容
为解决背景技术中所提出的技术问题,本发明提供一种基于视觉测量的新型振动提取与振型观测方法,利用结构边缘时序扰动提取边缘位置的振动信号,进而进行时间、空间上的解耦,实现对结构的振型进行观测的目的,在提升效率的同时降低了模态分析的难度与成本。
本发明采用以下技术方案实现:一种基于视觉测量的新型振动提取与振型观测方法,包括如下步骤:
步骤一、搭建测量平台,利用高速相机获取被测结构的运动视频数据;
步骤二、定位边缘区域,获取运动视频中结构边缘所有位置的时序像素切片,组成时序矩阵;
步骤三、基于偏差提取算法,提取结构边缘时序矩阵上亚像素精度的振动信号;
步骤四、利用变分模态分解算法对振动信号进行解耦,获取结构的各阶模态响应;
步骤五、根据振动理论,系统振动的总响应可以表示为对应权值的模态振型的线性叠加,根据解耦后的振动信号和各阶模态响应,获取结构各阶模态的空间权值;
步骤六、根据解耦后的振动信号和结构各阶模态的空间权值,可实现结构的振型观测。
具体地,所述步骤三中,按如下流程进行:
假设合成图像为I(x,y)且沿着x轴方向有一条平整的直线边缘,I(x,y)在y方向上做单自由度运动,组成I(x,y)的像素共有N列,提取时序上所有图像帧中第n列xn的像素灰度值作为采样图像且/>时序上的运动信号/>仅与t坐标值相关,因此的边缘轮廓可以用来表征/>定义图像t时刻xn处的像素采样为/>
若时序上无运动则像素采样与t坐标无关,即/>若时序上存在运动/>则有
采样图像的共同边缘轮廓/>是通过聚合时序上所有可用的像素采样来计算的,又因为运动信号/>为一微小值,所以可以用边缘轮廓的平均值来计算假设图像噪声/>在每个像素独立分布,则/>可以表示为
对使用一阶泰勒展开:
则像素在t坐标处的平均值为:
其中μf为在t方向上的均值,Dt为t方向上的图像帧数
由于μf同样为一微小值,因此再次使用泰勒展开有:
至此,在t坐标处的平均值就可以通过参数μf进行近似;
在实际的计算中,我们使用近似表示/>则运动信号可以通过计算与实际观测值的最小平方误差实现
此时运动信号可以近似表达为:
其中,为结构在时序上所有图像帧中第n列xn的像素灰度值,/>为标准边缘轮廓的运动信号,/>为/>的一维数值梯度,/>为结构振动信号;
由此可得结构边缘的振动信号,此时提取的振动信号是耦合的。
具体地,所述步骤四中,利用变分模态分解算法对某一非节点位置的边缘振动信号进行分解,得到解耦后的分量,获取结构的各阶模态响应qi(t),并对固有频率和阻尼比等模态参数进行估计。
具体地,根据模态叠加原理,结构各个边缘位置的振动可以表示为各阶模态响应的线性组合
p(t)为结构振动信号,k为被激发的最高阶模态阶数,φi为结构第i阶模态空间权值系数,qi(t)为结构第i阶的模态响应,αi是衰减系数
因此结构的边缘振动信号和模态响应的关系可以表示为:
实际应用中,利用解耦后的振动信号以及各阶模态响应,根据上述公式可以计算出结构各阶模态的空间权值φi。
相比现有技术,本发明的有益效果在于:
本发明提出的基于视觉测量的新型振动提取与振型观测方法,首先无需在被测结构上布置接触式传感器,降低了测量难度和测量成本;其次通过视觉测量振动信号,提高了空间分辨率;最后利用偏差提取算法,仅关注结构边缘的时序扰动,提升了结构振动分析与振型观测的效率。
附图说明
图1为本发明的悬臂梁结构试验现场图;
图2为本发明的结构边缘时序矩阵图;
图3为本发明的耦合振动信号图;
图4为本发明的解耦后前三阶模态响应及频谱图;
图5为本发明的可视化结构振型;
具体实施方式
下面,结合附图以及具体实施方式,对本发明做进一步描述,需要说明的是,在不相冲突的前提下,以下描述的各实施例之间或各技术特征之间可以任意组合形成新的实施例。
实施例1:
参照图1-5,本方案提出的一种基于视觉测量的新型振动提取与振型观测方法,包括如下步骤:
步骤一、搭建测量平台,利用高速相机获取被测结构的运动视频数据;如图1所示,即悬臂梁结构试验现场图。在实际应用中,被测结构件可以为梁、柱、剪力墙、桁架、车身等结构件。可以理解的是,本发明实施例并不对结构件的具体形式进行限定,任何用于受力或者起到支撑、连接作用的物理构件均可以作为结构件。
步骤二、定位边缘区域,获取运动视频中结构边缘所有位置的时序像素切片,组成时序矩阵;
步骤三、基于偏差提取算法,提取结构边缘时序矩阵上亚像素精度的振动信号;
具体地,所述步骤三中,按如下流程进行:
假设合成图像为I(x,y)且沿着x轴方向有一条平整的直线边缘,I(x,y)在y方向上做单自由度运动,组成I(x,y)的像素共有N列,提取时序上所有图像帧中第n列xn的像素灰度值作为采样图像且/>时序上的运动信号/>仅与t坐标值相关,因此的边缘轮廓可以用来表征/>定义图像t时刻xn处的像素采样为/>
若时序上无运动则像素采样与t坐标无关,即/>若时序上存在运动/>则有
采样图像的共同边缘轮廓/>是通过聚合时序上所有可用的像素采样来计算的,又因为运动信号/>为一微小值,所以可以用边缘轮廓的平均值来计算假设图像噪声/>在每个像素独立分布,则/>可以表示为
对使用一阶泰勒展开:
则像素在t坐标处的平均值为:
其中μf为在t方向上的均值,Dt为t方向上的图像帧数
由于μf同样为一微小值,因此再次使用泰勒展开有:
至此,在t坐标处的平均值就可以通过参数μf进行近似;
在实际的计算中,我们使用近似表示/>则运动信号可以通过计算与实际观测值的最小平方误差实现
此时运动信号fxn(t)可以近似表达为:
其中,为结构在时序上所有图像帧中第n列xn的像素灰度值,/>为标准边缘轮廓的运动信号,/>为/>的一维数值梯度,/>为结构振动信号;
由此可得结构边缘的振动信号,此时提取的振动信号是耦合的,如图3所示
步骤四、利用变分模态分解算法对振动信号进行解耦,获取结构的各阶模态响应;具体地,所述步骤四中,利用变分模态分解算法对某一非节点位置的边缘振动信号进行分解,得到解耦后的分量,获取结构的各阶模态响应qi(t),并对固有频率和阻尼比等模态参数进行估计。qi(t)如图4所示。
步骤五、根据振动理论,系统振动的总响应可以表示为对应权值的模态振型的线性叠加,根据解耦后的振动信号和各阶模态响应,获取结构各阶模态的空间权值;
具体地,根据模态叠加原理,结构各个边缘位置的振动可以表示为各阶模态响应的线性组合
p(t)为结构振动信号,k为被激发的最高阶模态阶数,φi为结构第i阶模态空间权值系数,qi(t)为结构第i阶的模态响应,αi是衰减系数
因此结构的边缘振动信号和模态响应的关系可以表示为:
实际应用中,利用解耦后的振动信号以及各阶模态响应,根据上述公式可以计算出结构各阶模态的空间权值φi。
步骤六、根据解耦后的振动信号和结构各阶模态的空间权值,可实现结构的振型观测。如图5所示。
上述实施方式仅为本发明的优选实施方式,不能以此来限定本发明保护的范围,本领域的技术人员在本发明的基础上所做的任何非实质性的变化及替换均属于本发明所要求保护的范围。
Claims (4)
1.一种基于视觉测量的新型振动提取与振型观测方法,其特征在于,包括如下步骤:
步骤一、搭建测量平台,利用高速相机获取被测结构的运动视频数据;
步骤二、定位边缘区域,获取运动视频中结构边缘所有位置的时序像素切片,组成时序矩阵;
步骤三、基于偏差提取算法,提取结构边缘时序矩阵上亚像素精度的振动信号;
步骤四、利用变分模态分解算法对振动信号进行解耦,获取结构的各阶模态响应;
步骤五、根据振动理论,系统振动的总响应可以表示为对应权值的模态振型的线性叠加,根据解耦后的振动信号和各阶模态响应,获取结构各阶模态的空间权值;
步骤六、根据解耦后的振动信号和结构各阶模态的空间权值,可实现结构的振型观测。
2.如权利要求1所述的一种基于视觉测量的新型振动提取与振型观测方法,其特征在于,所述步骤三中,具体地:
假设合成图像为I(x,y)且沿着x轴方向有一条平整的直线边缘,I(x,y)在y方向上做单自由度运动,组成I(x,y)的像素共有N列,提取时序上所有图像帧中第n列xn的像素灰度值作为采样图像且/>时序上的运动信号/>仅与t坐标值相关,因此/>的边缘轮廓可以用来表征/>定义图像t时刻xn处的像素采样为/>
若时序上无运动则像素采样与t坐标无关,即/>若时序上存在运动/>则有
采样图像的共同边缘轮廓/>是通过聚合时序上所有可用的像素采样来计算的,又因为运动信号/>为一微小值,所以可以用边缘轮廓的平均值来计算假设图像噪声/>在每个像素独立分布,则/>可以表示为
对使用一阶泰勒展开:
则像素在t坐标处的平均值为:
其中μf为在t方向上的均值,Dt为t方向上的图像帧数
由于μf同样为一微小值,因此再次使用泰勒展开有:
至此,在t坐标处的平均值就可以通过参数μf进行近似;
在实际的计算中,我们使用近似表示/>则运动信号可以通过计算与实际观测值的最小平方误差实现
此时运动信号可以近似表达为:
其中,为结构在时序上所有图像帧中第n列xn的像素灰度值,/>为标准边缘轮廓的运动信号,/>为/>的一维数值梯度,/>为结构振动信号;
由此可得结构边缘的振动信号,此时提取的振动信号是耦合的。
3.如权利要求1所述的一种基于视觉测量的新型振动提取与振型观测方法,其特征在于,所述步骤四中,利用变分模态分解算法对某一非节点位置的边缘振动信号进行分解,得到解耦后的分量,获取结构的各阶模态响应qi(t),并对固有频率和阻尼比等模态参数进行估计。
4.如权利要求1所述的,其特征在于,所述步骤五中,根据模态叠加原理,结构各个边缘位置的振动可以表示为各阶模态响应的线性组合
p(t)为结构振动信号,k为被激发的最高阶模态阶数,φi为结构第i阶模态空间权值系数,qi(t)为结构第i阶的模态响应,αi是衰减系数
因此结构的边缘振动信号和模态响应的关系可以表示为:
实际应用中,利用解耦后的振动信号以及各阶模态响应,根据上述公式可以计算出结构各阶模态的空间权值φi。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311446779.5A CN117490823B (zh) | 2023-11-02 | 2023-11-02 | 一种基于视觉测量的振动提取与振型观测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311446779.5A CN117490823B (zh) | 2023-11-02 | 2023-11-02 | 一种基于视觉测量的振动提取与振型观测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117490823A true CN117490823A (zh) | 2024-02-02 |
CN117490823B CN117490823B (zh) | 2024-07-05 |
Family
ID=89675728
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311446779.5A Active CN117490823B (zh) | 2023-11-02 | 2023-11-02 | 一种基于视觉测量的振动提取与振型观测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117490823B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020168589A1 (zh) * | 2019-02-19 | 2020-08-27 | 大连理工大学 | 一种自动检测高速铁路桥梁自由振动响应以识别模态的方法 |
CN115046717A (zh) * | 2022-06-08 | 2022-09-13 | 安徽农业大学 | 一种利用傅里叶变换轮廓术优化的结构振型可视化方法、装置及系统 |
GB202215178D0 (en) * | 2022-03-03 | 2022-11-30 | Univ Hefei Technology | Computer vision-based dynamic bridge shape recognition method |
CN115797335A (zh) * | 2023-01-31 | 2023-03-14 | 武汉地震工程研究院有限公司 | 用于桥梁振动测量的欧拉运动放大效果评估及优化方法 |
CN116012760A (zh) * | 2023-01-10 | 2023-04-25 | 国网安徽省电力有限公司马鞍山供电公司 | 一种基于欧拉-拉格朗日混合框架的结构振型可视化方法 |
-
2023
- 2023-11-02 CN CN202311446779.5A patent/CN117490823B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020168589A1 (zh) * | 2019-02-19 | 2020-08-27 | 大连理工大学 | 一种自动检测高速铁路桥梁自由振动响应以识别模态的方法 |
GB202215178D0 (en) * | 2022-03-03 | 2022-11-30 | Univ Hefei Technology | Computer vision-based dynamic bridge shape recognition method |
CN115046717A (zh) * | 2022-06-08 | 2022-09-13 | 安徽农业大学 | 一种利用傅里叶变换轮廓术优化的结构振型可视化方法、装置及系统 |
CN116012760A (zh) * | 2023-01-10 | 2023-04-25 | 国网安徽省电力有限公司马鞍山供电公司 | 一种基于欧拉-拉格朗日混合框架的结构振型可视化方法 |
CN115797335A (zh) * | 2023-01-31 | 2023-03-14 | 武汉地震工程研究院有限公司 | 用于桥梁振动测量的欧拉运动放大效果评估及优化方法 |
Non-Patent Citations (2)
Title |
---|
王昊: "基于视觉测量及振动理论的桥梁结构损伤识别研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》, 15 May 2021 (2021-05-15), pages 1 - 79 * |
胡朝辉;袁向荣;刘敏;: "振动试验视频图像测试技术", 噪声与振动控制, no. 03, 18 June 2011 (2011-06-18), pages 162 - 165 * |
Also Published As
Publication number | Publication date |
---|---|
CN117490823B (zh) | 2024-07-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4660779B2 (ja) | 移動装置の位置誤差評価方法およびその評価結果に基づく移動精度向上方法 | |
CN109282785B (zh) | 一种基于应变监测的弹性支承板的变形监测方法 | |
GB2394543A (en) | Positional measurement of a feature within an image | |
CN111174961A (zh) | 一种基于模态分析的索力光学测量方法及其测量系统 | |
Bornert et al. | Digital image correlation (chapter 6) | |
Cabaleiro et al. | A case study of measurements of deformations due to different loads in pieces less than 1 m from lidar data | |
Mangini et al. | Tag recognition: A new methodology for the structural monitoring of cultural heritage | |
CN117490823B (zh) | 一种基于视觉测量的振动提取与振型观测方法 | |
CN108983702B (zh) | 基于计算机显微视觉切片扫描技术的显微视觉系统的显微视场数字化扩展方法及系统 | |
CN112525085A (zh) | 基于三角形网格化图像技术的节点位移与应变测量方法 | |
Sali et al. | Calculation of plates on elastic foundation by the generalized equations of finite difference method | |
JP2019203724A (ja) | 建物基礎の管理方法と建物基礎群の管理方法 | |
Banks et al. | Detection of non-symmetrical damage in smart plate-like structures | |
Echerfaoui et al. | Laser interferometer based measurement for positioning error compensation in cartesian multi-axis systems | |
Hohmann et al. | Digital image correlation (dic): an advanced nondestructive testing method for life extension of nuclear power plants | |
CN113063356B (zh) | 测量材料定向形变离面位移的方法 | |
CN118095016B (zh) | 基于逆元法与有限元刚度矩阵构建的载荷识别方法 | |
Huntley et al. | Large-scale full-field metrology using projected fringes: some challenges and solutions | |
Le Blanc et al. | Image correlation applied to single crystal plasticity experiments and comparison to strain gage data | |
CN112380595B (zh) | 超高层结构的温度变形预测模型建立方法及预测方法 | |
CN113656959B (zh) | 一种基于光学测量的含节理岩石动态断裂模式的分类方法 | |
CN118463807B (zh) | 大理石平台运动位置温度误差标定及补偿方法、系统 | |
CN115096529B (zh) | 一种桥梁动挠度分布式测量装置及测量方法 | |
CN114510762B (zh) | 基于时间序列模型系数灵敏度的结构损伤识别方法及系统 | |
Chen et al. | Microdeformation measurement via single-element extended digital image correlation |
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 |