CN114777986A - A method and system for detecting looseness of transmission tower bolts based on correlation coefficient - Google Patents
A method and system for detecting looseness of transmission tower bolts based on correlation coefficient Download PDFInfo
- Publication number
- CN114777986A CN114777986A CN202210238911.2A CN202210238911A CN114777986A CN 114777986 A CN114777986 A CN 114777986A CN 202210238911 A CN202210238911 A CN 202210238911A CN 114777986 A CN114777986 A CN 114777986A
- Authority
- CN
- China
- Prior art keywords
- tower
- measured
- correlation coefficient
- signal
- modal
- 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.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L5/00—Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes
- G01L5/24—Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes for determining value of torque or twisting moment for tightening a nut or other member which is similarly stressed
- G01L5/246—Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes for determining value of torque or twisting moment for tightening a nut or other member which is similarly stressed using acoustic waves
Landscapes
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
Description
技术领域technical field
本发明涉及螺栓松动检测技术领域,更具体的说是涉及一种基于相关系数的输电铁塔螺栓松动检测方法。The invention relates to the technical field of bolt loosening detection, in particular to a method for detecting bolt loosening of a transmission tower based on a correlation coefficient.
背景技术Background technique
输电铁塔工作环境复杂,常年受到外界载荷反复作用,使得螺栓连接的预紧力下降甚至出现松动,从而导致输电铁塔的整体刚度下降,严重者可能出现输电铁塔倒塌的事故。因此研究一种快速有效的方法应用于检测输电铁塔的螺栓连接结构的松动情况势在必行。The working environment of the transmission tower is complex, and it is repeatedly acted by external loads all year round, which reduces the pre-tightening force of the bolted connection or even loosens, which leads to the decrease of the overall rigidity of the transmission tower. In severe cases, the transmission tower may collapse. Therefore, it is imperative to study a fast and effective method to detect the looseness of the bolted connection structure of the transmission tower.
针对螺栓连接结构松动问题,许多学者做了大量的仿真和实验研究。在检测方法上,周靖等将机器视觉原理用于螺栓松动旋转角度测量,进而判断螺栓状态,该方法拥有良好的精度,但输电铁塔的高度和螺栓分布情况会影响其视觉信号的获取;李路兵将超声波无损检测方法用于螺栓探伤,该方法具有抗干扰能力强、分辨率高等特点,但该方法对实验仪器的要求比较苛刻;江文强等通过观察螺栓在横向载荷下的应变过程来判断螺栓松动的过程,该方法可以精准地对螺栓轴向力进行测量从而判断螺栓的状态,但是在工程实际中,螺栓松动是长期的过程,利用应变片对其进行实时监测不太现实。Many scholars have done a lot of simulation and experimental research on the loosening of bolted structures. In the detection method, Zhou Jing et al. applied the principle of machine vision to the measurement of the bolt loosening rotation angle, and then judged the bolt state. This method has good accuracy, but the height of the transmission tower and the bolt distribution will affect the acquisition of its visual signal; Li Lubing; The ultrasonic non-destructive testing method is used for bolt flaw detection. This method has the characteristics of strong anti-interference ability and high resolution, but this method has strict requirements on experimental instruments; Jiang Wenqiang et al. This method can accurately measure the axial force of the bolt to judge the state of the bolt, but in engineering practice, the loosening of the bolt is a long-term process, and it is not realistic to use strain gauges to monitor it in real time.
目前采用最为广泛的方法为基于振动响应的检测法,但是由于现场测试中会受到各种噪声影响,因此在信号处理上,周文强等将经验模态分解(EMD)用于螺栓连接结构振动信号处理上,但EMD算法存在许多的问题,比如效率低下、模式混叠和端点效应等缺陷;李启月等将EEMD算法应用于爆破振动信号降噪上,虽然EEMD分解很好地解决了EMD算法地模态混叠地问题,但其会影响降噪效果。At present, the most widely used method is the detection method based on vibration response. However, due to the influence of various noises in the field test, Zhou Wenqiang et al. applied empirical mode decomposition (EMD) to the vibration signal of bolted connection structures in signal processing. In terms of processing, the EMD algorithm has many problems, such as low efficiency, mode aliasing and end-point effects. Li Qiyue et al. applied the EEMD algorithm to the noise reduction of blasting vibration signals. state aliasing problem, but it will affect the noise reduction effect.
因此,如何提供一种基于相关系数的能够有效降低噪声影响便于更精确检测输电铁塔螺栓松动的方法及系统是本领域技术人员亟需解决的问题。Therefore, how to provide a method and system based on the correlation coefficient that can effectively reduce the influence of noise and facilitate more accurate detection of the looseness of the transmission tower bolts is an urgent problem to be solved by those skilled in the art.
发明内容SUMMARY OF THE INVENTION
有鉴于此,本发明提供了一种基于相关系数的输电铁塔螺栓松动检测方法及系统以解决背景技术中提到的问题。In view of this, the present invention provides a method and system for detecting the looseness of a transmission tower bolt based on a correlation coefficient to solve the problems mentioned in the background art.
为了实现上述目的,本发明采用如下技术方案:In order to achieve the above object, the present invention adopts the following technical solutions:
一种基于相关系数的输电铁塔螺栓松动检测方法,从与输电铁塔螺栓连接的塔材中任选两个塔材作为塔材1和塔材2,包括以下步骤:A method for detecting looseness of transmission tower bolts based on correlation coefficient, selecting two tower materials as
S1.分别在塔材1和塔材2上各选取待测点,采集各个待测点的振动信号,每个塔材选取的待测点数量至少为1;S1. Select points to be measured on
S2.获取各个所述待测点的所述振动信号的时域图像和频域图像;S2. obtain the time domain image and frequency domain image of the vibration signal of each described point to be measured;
S3.求解所述塔材1的各个所述待测点分别与所述塔材2的各个待测点在时域和频域上的相关系数;S3. Solve the correlation coefficients between the points to be measured of the
S4.根据所述相关系数的变化检测输电铁塔螺栓松动情况,当所述相关系数的增大比例大于设定阈值时,则输电铁塔螺栓松动。S4. Detect the looseness of the bolts of the transmission tower according to the change of the correlation coefficient. When the increase ratio of the correlation coefficient is greater than the set threshold, the bolts of the transmission tower are loose.
优选的,S2中获取各个所述待测点的所述振动信号的时域图像和频域图像过程中还包括对所述振动信号进行降噪处理,所述降噪处理的具体内容包括:Preferably, the process of acquiring the time domain image and the frequency domain image of the vibration signal of each of the points to be measured in S2 also includes performing noise reduction processing on the vibration signal, and the specific content of the noise reduction processing includes:
S21.对各个所述振动信号进行CEEMDAN分解,分别得到K个模态分量;S21. carry out CEEMDAN decomposition to each described vibration signal, obtain K modal components respectively;
S22.计算各个所述模态分量的方差贡献率,并根据所述方差贡献率对各个所述模态分量进行筛选,剔除噪声分量;S22. Calculate the variance contribution rate of each of the modal components, and screen each of the modal components according to the variance contribution rate to eliminate noise components;
S23.将剩余的所述模态分量进行重组,并利用自相关函数进一步降噪处理,获得所述去噪信号。S23. Recombining the remaining modal components, and further denoising by using an autocorrelation function to obtain the denoised signal.
优选的,S21的具体内容包括:Preferably, the specific content of S21 includes:
S211.将正负高斯白噪声信号加入到所述振动信号中得到重构信号,然后将所述重构信号进行EMD分解,得到第一阶模态分量C1:S211. Add positive and negative Gaussian white noise signals to the vibration signal to obtain a reconstructed signal, and then perform EMD decomposition on the reconstructed signal to obtain the first-order modal component C 1 :
S212.求解CEEMDAN分解的第一个模态分量:S212. Solve the first modal component of the CEEMDAN decomposition:
S213.计算第1个残差:S213. Calculate the first residual:
S214.令y(t)=r1(t),重复S211、S212和S213,得到第2个模态分量和第2个残差分量:S214. Let y(t)=r 1 (t), repeat S211, S212 and S213 to obtain the second modal component and the second residual component:
S215.重复S211、S212、S213和S214,当得到的残差分量为单调函数时,分解停止,算法结束,设此时得到的模态分量个数为K,则所述振动信号y(t)分解为:S215. Repeat S211, S212, S213 and S214, when the obtained residual component is a monotonic function, the decomposition is stopped and the algorithm ends, and the number of modal components obtained at this time is set to K, then the vibration signal y(t) Decomposed into:
其中,y(t)为所述振动信号,y(t)+(-1)q∈υj(t)为所述重构信号,q=1,2,Ei(·)为EMD分解的第i个模态分量,为CEEMDAN分解的第k模态分量,νj是符合正态分布的高斯白噪声,j=1,2,...,N为添加白噪声的次数,ε为白噪声的标准表。Among them, y(t) is the vibration signal, y(t)+(-1) q ∈υ j (t) is the reconstructed signal, q=1, 2, E i (·) is the EMD decomposed signal the i-th modal component, is the kth modal component of CEEMDAN decomposition, ν j is Gaussian white noise conforming to the normal distribution, j=1, 2,..., N is the number of times of adding white noise, and ε is the standard table of white noise.
优选的,S22的具体内容包括:Preferably, the specific content of S22 includes:
计算各模态分量的方差D(k)从而得到各模态分量的方差贡献率e(k):Calculate the variance D(k) of each modal component to obtain the variance contribution rate e(k) of each modal component:
优选的,S3中Preferably, in S3
其中,ρxy为所述塔材1的各个所述待测点分别与所述塔材2的各个待测点的相关系数,μx和μy分别为所述塔材1和所述塔材2待测点的所述振动信号的数学期望,σx和σy分别为所述塔材1和所述塔材2待测点的所述振动信号的标准差。Wherein, ρ xy is the correlation coefficient between each point to be measured of the
一种基于相关系数的输电铁塔螺栓松动检测系统,包括信号采集模块、图像获取模块、相关系数求解模块和判定模块;A correlation coefficient-based bolt loosening detection system for a transmission tower includes a signal acquisition module, an image acquisition module, a correlation coefficient solution module and a determination module;
所述信号采集模块,用于分别从与输电铁塔螺栓连接的塔材1和塔材2各选取的待测点采集各个待测点的振动信号,每个塔材选取的待测点数量至少为1;The signal acquisition module is used to collect the vibration signal of each point to be measured from the points to be measured selected from the
所述图像获取模块,用于获取各个所述待测点的所述振动信号的时域图像和频域图像;The image acquisition module is used to acquire the time domain image and the frequency domain image of the vibration signal of each of the to-be-measured points;
所述相关系数求解模块,用于求解所述塔材1的各个所述待测点分别与所述塔材2的各个待测点在时域和频域上的相关系数;The correlation coefficient solving module is used to solve the correlation coefficients between the points to be measured of the
所述判定模块,用于根据所述相关系数的变化检测输电铁塔螺栓松动情况,当所述相关系数的变化比例大于设定阈值时,则输电铁塔螺栓松动。The determining module is used for detecting the looseness of the bolts of the transmission tower according to the change of the correlation coefficient, and when the change ratio of the correlation coefficient is greater than the set threshold, the bolts of the transmission tower are loose.
优选的,所述图像获取模块包括信号分解单元、模态分量筛选单元和重组降噪单元;Preferably, the image acquisition module includes a signal decomposition unit, a modal component screening unit and a recombination noise reduction unit;
所述信号分解单元,用于对各个所述振动信号进行CEEMDAN分解分别得到K个模态分量;The signal decomposition unit is used to perform CEEMDAN decomposition on each of the vibration signals to obtain K modal components respectively;
所述模态分量筛选单元,用于计算各个所述模态分量的方差贡献率,并根据所述方差贡献率对各个所述模态分量进行筛选,剔除噪声分量;The modal component screening unit is configured to calculate the variance contribution rate of each of the modal components, and screen each of the modal components according to the variance contribution rate to eliminate noise components;
所述重组降噪单元,用于将剩余的所述模态分量进行重组,并利用自相关函数进一步降噪处理,获得所述去噪信号。The recombination denoising unit is used for recombining the remaining modal components, and further denoising by using an autocorrelation function to obtain the denoising signal.
优选的,所述信号分解单元包括EMD分解子单元、CEEMDAN分解子单元和残差计算子单元;Preferably, the signal decomposition unit includes an EMD decomposition subunit, a CEEMDAN decomposition subunit and a residual calculation subunit;
所述EMD分解子单元,用于将正负高斯白噪声信号加入到所述振动信号中得到重构信号,并将所述重构信号进行EMD分解获得EMD模态分量;The EMD decomposition subunit is used to add positive and negative Gaussian white noise signals to the vibration signal to obtain a reconstructed signal, and perform EMD decomposition on the reconstructed signal to obtain an EMD modal component;
所述CEEMDAN分解子单元,用于求解CEEMDAN分解的模态分量;The CEEMDAN decomposition subunit is used to solve the modal component of the CEEMDAN decomposition;
所述残差计算子单元,用于计算残差。The residual calculation subunit is used to calculate the residual.
优选的,所述模态分量筛选单元包括方差计算子单元、方差贡献率计算子单元和剔除噪声子单元;Preferably, the modal component screening unit includes a variance calculation subunit, a variance contribution rate calculation subunit, and a noise elimination subunit;
所述方差计算子单元,用于计算各模态分量的方差;The variance calculation subunit is used to calculate the variance of each modal component;
所述方差贡献率计算子单元,用于根据所述方差计算各模态分量的方差贡献率;The variance contribution rate calculation subunit is used to calculate the variance contribution rate of each modal component according to the variance;
所述剔除噪声子单元,用于根据所述方差贡献率对各个所述模态分量进行筛选,剔除噪声分量。The noise-removing subunit is configured to screen each of the modal components according to the variance contribution rate, and remove the noise component.
优选的,所述的一种基于相关系数的输电铁塔螺栓松动检测系统,还包括预警模块,当所述判定模块判定输电铁塔螺栓松动时进行报警。Preferably, the system for detecting the looseness of the bolts of the transmission tower based on the correlation coefficient further includes an early warning module, which will give an alarm when the determination module determines that the bolts of the transmission tower are loose.
经由上述的技术方案可知,与现有技术相比,本发明公开提供了一种基于相关系数的输电铁塔螺栓松动检测方法及系统,通过对现场测得振动信号进行CEEMDAN-自相关函数降噪处理,使得振动信号共振频率对应的峰值得到很好的保留,非共振信息得到抑制,便于消除噪声的影响;通过相关系数分析法来检测输电铁塔螺栓松动情况,操作简单,快速有效,判定更准确。It can be seen from the above technical solutions that, compared with the prior art, the present invention provides a method and system for detecting the looseness of a transmission tower bolt based on a correlation coefficient. , so that the peak value corresponding to the resonance frequency of the vibration signal is well preserved, the non-resonant information is suppressed, and it is convenient to eliminate the influence of noise; through the correlation coefficient analysis method to detect the looseness of the transmission tower bolts, the operation is simple, fast and effective, and the determination is more accurate.
附图说明Description of drawings
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图获得其他的附图。In order to explain the embodiments of the present invention or the technical solutions in the prior art more clearly, the following briefly introduces the accompanying drawings that need to be used in the description of the embodiments or the prior art. Obviously, the accompanying drawings in the following description are only It is an embodiment of the present invention. For those of ordinary skill in the art, other drawings can also be obtained according to the provided drawings without creative work.
图1附图为本发明提供的输电铁塔螺栓松动检测方法示意图;Fig. 1 accompanying drawing is a schematic diagram of a method for detecting the looseness of a transmission tower bolt provided by the present invention;
图2附图为本发明提供的输电铁塔待测螺栓及待测点示意图;The accompanying drawing of Fig. 2 is a schematic diagram of a bolt to be measured and a point to be measured of a transmission tower provided by the present invention;
图3附图为本发明提供的CEEMDAN-自相关函数降噪的流程图;Fig. 3 accompanying drawing is the flow chart of CEEMDAN-autocorrelation function noise reduction provided by the present invention;
图4附图为本发明提供的实施例测点3的原始振动信号示意图;4 is a schematic diagram of the original vibration signal of the
图5附图为本发明提供的实施例测点3的原始信号和降噪信号对比图;The accompanying drawing of FIG. 5 is a comparison diagram of the original signal and the noise reduction signal of the
图6附图为本发明提供的实施例测点3的CEEMDAN分解图;Fig. 6 accompanying drawing is the CEEMDAN exploded view of
图7附图为本发明提供的实施例螺栓紧固时测点1-4降噪后振动时域图;Fig. 7 accompanying drawing is the vibration time domain diagram after noise reduction of measuring points 1-4 when bolts are tightened according to the embodiment of the present invention;
图8附图为本发明提供的实施例螺栓松动时测点1-4降噪后振动时域图;The accompanying drawing of FIG. 8 is a time domain diagram of vibration after noise reduction of measuring points 1-4 when the bolts are loosened according to the embodiment of the present invention;
图9附图为本发明提供的实施例螺栓紧固时测点1-4降噪后振动频域图;Fig. 9 accompanying drawing is the vibration frequency domain diagram after noise reduction of measuring points 1-4 when bolts are tightened according to the embodiment of the present invention;
图10附图为本发明提供的实施例螺栓松动时测点1-4降噪后振动频域图。FIG. 10 is a diagram of the vibration frequency domain after noise reduction of measuring points 1-4 when the bolts are loosened according to the embodiment of the present invention.
具体实施方式Detailed ways
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only a part of the embodiments of the present invention, but not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts shall fall within the protection scope of the present invention.
本发明实施例公开了一种基于相关系数的输电铁塔螺栓松动检测方法,从与输电铁塔螺栓连接的塔材中任选两个塔材作为塔材1和塔材2,如图1,包括以下步骤:The embodiment of the present invention discloses a method for detecting looseness of transmission tower bolts based on correlation coefficient. Two tower materials are selected from the tower materials connected with the transmission tower bolts as
S1.分别在塔材1和塔材2上各选取待测点,采集各个待测点的振动信号,每个塔材选取的待测点数量至少为1;S1. Select points to be measured on
S2.获取各个待测点的振动信号的时域图像和频域图像;S2. Obtain the time domain image and frequency domain image of the vibration signal of each point to be measured;
S3.求解塔材1的各个待测点分别与塔材2的各个待测点在时域和频域上的相关系数;S3. Solve the correlation coefficient between each point to be measured of
S4.根据相关系数的变化检测输电铁塔螺栓松动情况,当所述相关系数的增大比例大于设定阈值时,则输电铁塔螺栓松动。S4. Detect the looseness of the bolts of the transmission tower according to the change of the correlation coefficient. When the increase ratio of the correlation coefficient is greater than the set threshold, the bolts of the transmission tower are loose.
在本实施例中,选取灵敏度为4.44pC/N,增益为0.500mV/Pv的激振锤作为激振源对输电铁塔施加脉冲激励,在输电铁塔螺栓的两个连接板上分别选取两个测点,采用灵敏度为5.072mv/g的压电加速度传感器,设置信号采集仪的采样频率为5kHz。In this embodiment, an excitation hammer with a sensitivity of 4.44pC/N and a gain of 0.500mV/Pv is selected as the excitation source to apply pulse excitation to the transmission tower. point, use a piezoelectric acceleration sensor with a sensitivity of 5.072mv/g, and set the sampling frequency of the signal acquisition instrument to 5kHz.
在本实施例中,对铁塔待测螺栓进行紧固,实施S1-S4步骤,求解的相关系数作为初始相关系数,后期通过S1-S4测得的待测点的相关系数与初始相关系数对比,增大的比例大于设定阈值时,则输电铁塔螺栓松动。In this embodiment, the bolts to be measured of the iron tower are fastened, the steps S1-S4 are implemented, and the solved correlation coefficient is used as the initial correlation coefficient, and the correlation coefficient of the point to be measured measured by S1-S4 in the later stage is compared with the initial correlation coefficient, When the increased ratio is greater than the set threshold, the bolts of the transmission tower are loose.
为了进一步实施上述技术方案,如图2,S2中获取各个待测点的振动信号的时域图像和频域图像过程中还包括对振动信号进行降噪处理,降噪处理的具体内容包括:In order to further implement the above technical solution, as shown in Figure 2, the process of acquiring the time domain image and frequency domain image of the vibration signal of each point to be measured in S2 also includes noise reduction processing on the vibration signal, and the specific content of the noise reduction processing includes:
S21.对各个振动信号进行CEEMDAN分解,分别得到K个模态分量;S21. Perform CEEMDAN decomposition on each vibration signal to obtain K modal components respectively;
S22.计算各个模态分量的方差贡献率,并根据方差贡献率对各个模态分量进行筛选,剔除噪声分量;S22. Calculate the variance contribution rate of each modal component, and screen each modal component according to the variance contribution rate to eliminate noise components;
S23.将剩余的模态分量进行重组,并利用自相关函数进一步降噪处理,获得去噪信号。S23. Recombining the remaining modal components, and using the autocorrelation function for further noise reduction processing to obtain a denoised signal.
在实际应用中,CEEMDAN分解是EMD分解算法的改进算法,有效解决了EMD分解使模态混叠问题以及EEMD和CEEMD分解会残留白噪声的问题。In practical applications, CEEMDAN decomposition is an improved algorithm of EMD decomposition algorithm, which effectively solves the problem of mode aliasing caused by EMD decomposition and the problem of residual white noise caused by EEMD and CEEMD decomposition.
在本实施例中,将方差贡献率小于5%的模态分量剔除。In this embodiment, modal components whose variance contribution rate is less than 5% are eliminated.
为了进一步实施上述技术方案,S21的具体内容包括:In order to further implement the above technical solution, the specific content of S21 includes:
S211.将正负高斯白噪声信号加入到振动信号中得到重构信号,然后将重构信号进行EMD分解,得到第一阶模态分量C1:S211. Add positive and negative Gaussian white noise signals to the vibration signal to obtain a reconstructed signal, and then perform EMD decomposition on the reconstructed signal to obtain the first-order modal component C 1 :
S212.求解CEEMDAN分解的第一个模态分量:S212. Solve the first modal component of the CEEMDAN decomposition:
S213.计算第1个残差:S213. Calculate the first residual:
为了进一步实施上述技术方案,S21的具体内容还包括:In order to further implement the above technical solution, the specific content of S21 also includes:
S214.令y(t)=r1(t),重复S211、S212和S213,得到第2个模态分量和第2个残差分量:S214. Let y(t)=r 1 (t), repeat S211, S212 and S213 to obtain the second modal component and the second residual component:
S215.重复S211、S212、S213和S214,当得到的残差分量为单调函数时,分解停止,算法结束,设此时得到的模态分量个数为K,则振动信号y(t)分解为:S215. Repeat S211, S212, S213 and S214. When the obtained residual component is a monotonic function, the decomposition is stopped and the algorithm ends. Assuming that the number of modal components obtained at this time is K, the vibration signal y(t) is decomposed into :
其中,y(t)为振动信号,y(t)+(-1)q∈υj(t)为重构信号,q=1,2,Ei(·)为EMD分解的第i个模态分量,为CEEMDAN分解的第k模态分量,νj是符合正态分布的高斯白噪声,j=1,2,...,N为添加白噪声的次数,ε为白噪声的标准表。Among them, y(t) is the vibration signal, y(t)+(-1) q ∈υ j (t) is the reconstructed signal, q=1, 2, E i ( ) is the ith mode of the EMD decomposition state component, is the kth modal component of CEEMDAN decomposition, ν j is Gaussian white noise conforming to the normal distribution, j=1, 2,..., N is the number of times of adding white noise, and ε is the standard table of white noise.
在本实施例中,以测点3为例,首先利用EEMD分解中各参量的取值范围,采用试错法确定CEEMDAN分解所需的基本参量:设置正负高斯白噪声标准差为0.2,噪声加入的次数为500,允许的最大迭代次数为5000,最终将3测点的原始振动信号如图3,分解为11个模态参量(IMF1、...、IMF11)如图5。In this embodiment, taking
为了进一步实施上述技术方案,S22的具体内容包括:In order to further implement the above technical solution, the specific content of S22 includes:
计算各模态分量的方差D(k)从而得到各模态分量的方差贡献率e(k):Calculate the variance D(k) of each modal component to obtain the variance contribution rate e(k) of each modal component:
在本实施例中,根据上述方法计算测点3原始振动信号各模态分量的方差贡献率,如下表所示:In this embodiment, the variance contribution rate of each modal component of the original vibration signal of measuring
可以看出,11个分量中IMF3-IMF11的方差贡献率都在5%以下,说明这9组分量对原信号而言为噪声信号的干扰,所以可以直接将9组分量进行剔除,然后将剩余的IMF1和IMF2这两组分量进行重组,利用自相关函数进行相应的降噪处理,得到最后的去噪信号。It can be seen that the variance contribution rates of IMF3-IMF11 in the 11 components are all below 5%, indicating that these 9 components are the interference of noise signals to the original signal, so the 9 components can be directly eliminated, and then the remaining components The two components of IMF1 and IMF2 are recombined, and the autocorrelation function is used to perform corresponding noise reduction processing to obtain the final denoised signal.
按照上述测点3的降噪过程依次对各个测点进行降噪处理,得到螺栓松动前后各个测点滤波降噪处理后的振动时域图和频谱图,如图6-9。Perform noise reduction processing on each measurement point in turn according to the noise reduction process of
为了进一步实施上述技术方案,S3中In order to further implement the above technical solution, in S3
其中,ρxy为塔材1的各个待测点分别与塔材2的各个待测点的相关系数,μx和μy分别为塔材1和塔材2待测点的振动信号的数学期望,σx和σy分别为塔材1和塔材2待测点的振动信号的标准差。Wherein, ρ xy is the correlation coefficient of each point to be measured of
在本实施例中,塔材1和塔材2待测点的振动信号为经过降噪处理后的振动信号。In this embodiment, the vibration signals of the
在本实施例中,各测点之间的相关系数绝对值计算结果如下表所示:In this embodiment, the calculation result of the absolute value of the correlation coefficient between each measuring point is shown in the following table:
在本实施例中,当螺栓连接结构发生松动后,各个测点振动信号之间无论在时域上还是在频域上其相关系数升高,关联性增强。In this embodiment, when the bolted connection structure is loosened, the correlation coefficient between the vibration signals of each measuring point increases, whether in the time domain or in the frequency domain, and the correlation is enhanced.
在螺栓紧固时,塔材1和塔材2相当于一个整体,在脉冲激振的作用下,各个测点发生自激振动,测点振动信号之间的关联性不强;在螺栓发生完全松动后,塔材2不再和塔材1是一个整体,当用激振锤对塔材1施加冲击信号时,塔材1发生振动,碰撞塔材2,使得塔材2一起振动,而塔材2发生的振动是由塔材1引起的强迫振动,因此测点振动信号之间的关联性增强。When the bolts are tightened,
一种基于相关系数的输电铁塔螺栓松动检测系统,包括信号采集模块、图像获取模块、相关系数求解模块和判定模块;A correlation coefficient-based bolt loosening detection system for a transmission tower includes a signal acquisition module, an image acquisition module, a correlation coefficient solution module and a determination module;
信号采集模块,用于分别从与输电铁塔螺栓连接的塔材1和塔材2各选取的待测点采集各个待测点的振动信号,每个塔材选取的待测点数量至少为1;The signal acquisition module is used to collect the vibration signal of each to-be-measured point from the to-be-measured points selected from the
图像获取模块,用于获取各个待测点的振动信号的时域图像和频域图像;The image acquisition module is used to acquire the time domain image and frequency domain image of the vibration signal of each point to be measured;
相关系数求解模块,用于求解塔材1的各个待测点分别与塔材2的各个待测点在时域和频域上的相关系数;The correlation coefficient solving module is used to solve the correlation coefficient of each point to be measured of
判定模块,用于根据相关系数的变化检测输电铁塔螺栓松动情况,当所述相关系数的变化比例大于设定阈值时,则输电铁塔螺栓松动。The determination module is used to detect the looseness of the bolts of the transmission tower according to the change of the correlation coefficient. When the change ratio of the correlation coefficient is greater than the set threshold, the bolts of the transmission tower are loose.
为了进一步实施上述技术方案,图像获取模块包括信号分解单元、模态分量筛选单元和重组降噪单元;In order to further implement the above technical solution, the image acquisition module includes a signal decomposition unit, a modal component screening unit and a recombination noise reduction unit;
信号分解单元,用于对各个振动信号进行CEEMDAN分解分别得到K个模态分量;The signal decomposition unit is used to perform CEEMDAN decomposition on each vibration signal to obtain K modal components respectively;
模态分量筛选单元,用于计算各个模态分量的方差贡献率,并根据方差贡献率对各个模态分量进行筛选,剔除噪声分量;The modal component screening unit is used to calculate the variance contribution rate of each modal component, and screen each modal component according to the variance contribution rate to eliminate noise components;
重组降噪单元,用于将剩余的模态分量进行重组,并利用自相关函数进一步降噪处理,获得去噪信号。The recombination denoising unit is used for recombining the remaining modal components, and further denoising by using the autocorrelation function to obtain denoising signals.
为了进一步实施上述技术方案,信号分解单元包括EMD分解子单元、CEEMDAN分解子单元和残差计算子单元;In order to further implement the above technical solution, the signal decomposition unit includes an EMD decomposition subunit, a CEEMDAN decomposition subunit and a residual calculation subunit;
EMD分解子单元,用于将正负高斯白噪声信号加入到振动信号中得到重构信号,并将重构信号进行EMD分解获得EMD模态分量;The EMD decomposition sub-unit is used to add positive and negative Gaussian white noise signals to the vibration signal to obtain a reconstructed signal, and perform EMD decomposition on the reconstructed signal to obtain EMD modal components;
CEEMDAN分解子单元,用于求解CEEMDAN分解的模态分量;The CEEMDAN decomposition subunit is used to solve the modal components of the CEEMDAN decomposition;
残差计算子单元,用于计算残差。Residual calculation subunit for calculating residuals.
为了进一步实施上述技术方案,模态分量筛选单元包括方差计算子单元、方差贡献率计算子单元和剔除噪声子单元;In order to further implement the above technical solution, the modal component screening unit includes a variance calculation subunit, a variance contribution rate calculation subunit and a noise elimination subunit;
方差计算子单元,用于计算各模态分量的方差;The variance calculation subunit is used to calculate the variance of each modal component;
方差贡献率计算子单元,用于根据方差计算各模态分量的方差贡献率;The variance contribution rate calculation subunit is used to calculate the variance contribution rate of each modal component according to the variance;
剔除噪声子单元,用于根据方差贡献率对各个模态分量进行筛选,剔除噪声分量。The noise elimination subunit is used to filter each modal component according to the variance contribution rate, and eliminate the noise component.
为了进一步实施上述技术方案,一种基于相关系数的输电铁塔螺栓松动检测系统,还包括预警模块,当判定模块判定输电铁塔螺栓松动时进行报警。In order to further implement the above technical solution, a system for detecting looseness of transmission tower bolts based on a correlation coefficient further includes an early warning module, which gives an alarm when the determination module determines that the bolts of the transmission tower are loose.
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的装置而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。The various embodiments in this specification are described in a progressive manner, and each embodiment focuses on the differences from other embodiments, and the same and similar parts between the various embodiments can be referred to each other. For the device disclosed in the embodiment, since it corresponds to the method disclosed in the embodiment, the description is relatively simple, and the relevant part can be referred to the description of the method.
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。The above description of the disclosed embodiments enables any person skilled in the art to make or use the present invention. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be implemented in other embodiments without departing from the spirit or scope of the invention. Thus, the present invention is not intended to be limited to the embodiments shown herein, but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.
Claims (10)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210238911.2A CN114777986A (en) | 2022-03-11 | 2022-03-11 | A method and system for detecting looseness of transmission tower bolts based on correlation coefficient |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210238911.2A CN114777986A (en) | 2022-03-11 | 2022-03-11 | A method and system for detecting looseness of transmission tower bolts based on correlation coefficient |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114777986A true CN114777986A (en) | 2022-07-22 |
Family
ID=82424106
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210238911.2A Pending CN114777986A (en) | 2022-03-11 | 2022-03-11 | A method and system for detecting looseness of transmission tower bolts based on correlation coefficient |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114777986A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117782403A (en) * | 2024-02-27 | 2024-03-29 | 北京谛声科技有限责任公司 | Loose bolt positioning method, device and medium based on separation network |
-
2022
- 2022-03-11 CN CN202210238911.2A patent/CN114777986A/en active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117782403A (en) * | 2024-02-27 | 2024-03-29 | 北京谛声科技有限责任公司 | Loose bolt positioning method, device and medium based on separation network |
CN117782403B (en) * | 2024-02-27 | 2024-05-10 | 北京谛声科技有限责任公司 | Loose bolt positioning method, device and medium based on separation network |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Rivola et al. | Bispectral analysis of the bilinear oscillator with application to the detection of fatigue cracks | |
US8219361B2 (en) | Vibration analysis | |
TWI449883B (en) | Method for analyzing structure safety | |
RU2470280C2 (en) | Method for detection and automatic identification of rolling bearing damage | |
US5646350A (en) | Monitoring slow speed machinery using integrator and selective correction of frequency spectrum | |
US20140090472A1 (en) | Device and method for evaluation of vibrations | |
CN106768967B (en) | A kind of nondestructive detection method and system for looseness of flange fastening bolts | |
CN112461358B (en) | Bridge modal parameter identification method based on instantaneous frequency of vehicle-bridge system | |
CN110568081B (en) | Ancient building damage testing device and method | |
CN114777986A (en) | A method and system for detecting looseness of transmission tower bolts based on correlation coefficient | |
CN107748209A (en) | A kind of Structure Damage Identification | |
CN111780866B (en) | Diesel engine test base installation state natural frequency testing method and device | |
WO2020202567A1 (en) | Vibro-acoustic analysis method and device and equipment-abnormality-location estimation method and device | |
JP2006226716A (en) | Damage detection method and system of structure | |
CN113108973A (en) | Ultrasonic bolt pretension test system and test method | |
Goswami et al. | Efficient vibration signal denoising techniques for effective condition monitoring of gearbox | |
CN107192442A (en) | A kind of non-contact electric eddy sensor construction mode testing method | |
CN113358342A (en) | Bolt monitoring system and method for wind generating set | |
CN111486943B (en) | Industrial robot vibration state analysis method, system and device and readable storage medium | |
CN114084764A (en) | Elevator transportation quality detection method and detection system | |
CN101526417B (en) | Method for diagnosing structural damage by musical tone rule | |
JPH09269268A (en) | Bolt loosening detection method by strain waveform analysis | |
CN106629395B (en) | A kind of amplitude detection method of high pedestal jib crane Vibration Signal in Frequency Domain | |
JP2002348949A (en) | Dynamic seismic resistance assessment system for building | |
JPH01219653A (en) | Diagnostic apparatus of abnormality of gear device |
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 |