Method, system and device for identifying vortex-induced vibration of bridge

Abstract

The invention relates to the field of bridge structure performance monitoring, in particular to a method, system and device for recognizing vortex-induced vibration of a bridge, and the method comprises the steps: carrying out the mean removal and high-pass filtering of time history monitoring signal data # imgabs0 # of bridge vibration; determining a 0 phase point and a # imgabs2 phase point in the # imgabs1 #, and eliminating redundant data before the 0 phase point and redundant data after the # imgabs3 phase point to obtain a data sequence # imgabs4 # of a complete period; according to the # imgabs5 #, an amplitude sequence # imgabs7 # of an upper envelope line of the # imgabs6 # is obtained through calculation; the identification index # imgabs9 # of bridge vortex-induced vibration is obtained through calculation according to # imgabs8 #, # imgabs10 # is set to be a preset identification threshold value, when # imgabs11 # is larger than or equal to # imgabs12 #, vortex-induced vibration of the bridge is identified, the system adopts the method for identifying the bridge vortex-induced vibration, the device can execute the method for identifying the bridge vortex-induced vibration, and through the setting, the method for identifying the bridge vortex-induced vibration can be used for identifying the bridge vortex-induced vibration. The stability degree of the bridge vortex-induced vibration identification index is improved, and the calculation amount of the analysis process is reduced.

Description

Technical Field

[0001] The invention relates to the field of bridge structure performance monitoring, and in particular to a method, system and device for identifying vortex-induced vibration of a bridge. Background Art

[0002] During the operation of a bridge, vortex-induced vibrations may occur in the bridge due to the influence of the external environment. The vortex-induced vibrations of the bridge will damage the structure of the bridge, so it is necessary to monitor the vortex-induced vibrations of the bridge. Timely alarms are issued at the early stage of the vortex-induced vibration of the bridge and real-time monitoring of the entire process allows the maintenance unit to take emergency response measures in a timely manner; after the vortex-induced vibration ends, graded warnings are implemented to provide a basis for the maintenance unit to assess the safety of the bridge. These are the two main issues of concern in bridge vortex-induced vibration monitoring. Existing research methods for vortex-induced vibration identification mainly include various analysis methods such as energy factor method, Hilbert transform method, clustering algorithm and machine learning algorithm.

[0003] For example, the energy factor method uses the power spectrum of the bridge vibration acceleration signal and the principal component analysis method to evaluate the proportion of energy concentration coefficients corresponding to different modes, thereby determining the multi-modal vortex-induced vibration results of the main beam of the suspension bridge. This calculation method is relatively simple, but its disadvantage is that the index is also sensitive to non-vortex-induced signals, and the overall quantitative index is not conducive to alarm classification. When the quantitative index is set too high, the timeliness decreases; when it is too low, the false alarm rate increases significantly and the accuracy is low.

[0004] For another example, a machine learning algorithm trains a model by comparing the vortex-induced vibration characteristics of the target bridge with those of other bridges in the source domain to achieve intelligent identification of vortex-induced vibration. However, this method requires the acceleration amplitude, main frequency PSD amplitude, energy concentration coefficient R, and main frequency relative change D as feature inputs, and requires existing manually identified vortex-induced vibration signals as source domain samples. As a result, this method requires a large amount of calculation when identifying vortex-induced vibration of bridges, which reduces the efficiency of the identification work.

[0005] It can be seen that when the existing methods are used to identify bridge vortex-induced vibration, misjudgment may occur due to insufficient indicator stability, or the efficiency of identifying bridge vortex-induced vibration may be low due to large amount of calculation. Summary of the invention

[0006] The purpose of the present invention is to overcome the shortcomings of the prior art in which, when identifying vortex-induced vibration of a bridge, misjudgment may occur due to insufficient indicator stability, or the efficiency of identifying vortex-induced vibration of the bridge may be reduced due to large amount of calculation, and to provide a method, system and device for identifying vortex-induced vibration of the bridge.

[0007] In a first aspect, the present invention provides a method for identifying vortex-induced vibration of a bridge, comprising the following steps: S1: Time-course monitoring signal data of bridge vibration Perform averaging and high-pass filtering. , where i represents the time-course monitoring signal data The number of the data point in is a constant; S2: OK The 0 phase point and phase point, and remove the redundant data before the 0 phase point and Redundant data after the phase point to obtain a complete cycle data sequence , ,in, represents the number of complete cycles, A sequence of data representing a complete cycle The number of the data point in is a constant; S3: According to Calculated Amplitude sequence of the upper envelope , ; S4: According to The identification index of bridge vortex-induced vibration is calculated ,set up is the preset recognition threshold. Greater than or equal to When , it is identified that the bridge is experiencing vortex-induced vibration.

[0008] By monitoring the bridge vibration signal data over time De-averaging and high-pass filtering are performed to reduce the overall signal baseline offset caused by the DC component and the trend drift caused by low-frequency noise, thereby improving subsequent recognition accuracy.

[0009] This method is based on the time domain characteristics of bridge vortex-induced vibration, takes the signal envelope as the judgment basis, and realizes the monitoring and identification of bridge vortex-induced vibration events by picking up the envelope of the time domain signal and quantifying the identification indicators. Compared with the existing technology, this method has the advantages of high indicator stability, good indicator quantification effect, simple bridge vortex-induced vibration analysis process, small calculation amount, strong real-time performance, and visualization.

[0010] The step S2 comprises: S21: Establishing a rectangular coordinate system ,Will Each data point in is represented by a rectangular coordinate system and is obtained In the rectangular coordinate system The sequence of intersections of the axes , The expression is: ; S22: Settings The curve with a positive slope is The sequence of axis intersections is ,set up The curve with a negative slope is Sequence of axis intersections ,and ,in: The expression is: , The expression is: ; S23: Settings The sequence of all points on a curve with a positive slope is ,set up The sequence of all points on a curve with a negative slope is ,and , and Respectively expressed as: in, Representation sequence The minimum value in Representation sequence The maximum value in Representation sequence The minimum value in Representation sequence The maximum value in .

[0011] In step S3, the local maximum method, interpolation method, polynomial fitting method, Hilbert transform method, empirical mode decomposition method, or wavelet transform method are used to calculate Amplitude sequence of the upper envelope .

[0012] Calculated using the Hilbert transform method Amplitude sequence of the upper envelope , the step S3 comprises: S31: Yes Performing Hilbert transform, we get ,in, Representatives The sequence after Hilbert transformation, is a constant; S32: According to and Get the amplitude sequence of the upper envelope , It is expressed as: .

[0013] In step S4, the identification index pass The minimum value of The ratio of the average values ​​of is calculated. It is expressed as: ; in, express The minimum value in express The average value of .

[0014] Identify indicators pass The minimum value of This method improves the quantization effect of the stability of the signal amplitude, and the amount of calculation is small. Real-time calculation occupies less computing resources, which improves work efficiency.

[0015] The identification threshold Including the first recognition threshold , the second recognition threshold and the third recognition threshold ,and , when Greater than or equal to and Less than When , it is identified that the bridge is in the embryonic stage of vortex-induced vibration; when Greater than or equal to and Less than When , it is identified that the bridge is in the development stage of vortex-induced vibration; when Greater than or equal to When , it is identified that the bridge is in the stable stage of vortex-induced vibration.

[0016] Said The value of is 0.05-0.15, The value of is 0.45-0.55, The value of is 0.75-0.85. Preferably, The value of is 0.1, The value of is 0.5, The value of is 0.8.

[0017] Said The duration is 0.9-1.1min.

[0018] Since the vortex-induced vibration of the bridge is a development process, the data stability is higher within the time period of 0.9-1.1min. Therefore, selecting the data with a time period of 0.9-1.1min for identification can improve the accuracy of judging the vortex-induced vibration state of the bridge. Preferably, the The duration is 1min.

[0019] In a second aspect, the present invention provides a system for identifying vortex-induced vibration of a bridge, which adopts the method for identifying vortex-induced vibration of a bridge, including a bridge vibration signal processing module, a redundant data elimination module, an upper envelope calculation module, a bridge vortex-induced vibration identification index calculation module and a bridge vortex-induced vibration judgment module, wherein the bridge vibration signal processing module is used to perform de-averaging and high-pass filtering on the time-history monitoring signal data of the bridge vibration; The redundant data elimination module is used to eliminate redundant data to obtain complete cycle bridge vibration time history monitoring signal data; The upper envelope calculation module is used to calculate the upper envelope of the bridge vibration time history monitoring signal data of a complete cycle; The bridge vortex-induced vibration identification index calculation module is used to calculate the identification index of the bridge vortex-induced vibration; The bridge vortex-induced vibration judgment module is used to compare the bridge vortex-induced vibration identification index with a preset threshold value and obtain a judgment result.

[0020] Through this setting, the efficiency of identifying vortex-induced vibration of bridges is improved.

[0021] In a third aspect, the present invention provides a system device for identifying vortex-induced vibrations of bridges, comprising at least one processor and a memory communicatively connected to the at least one processor; the memory stores instructions executable by the at least one processor, and the instructions are executed by the at least one processor to enable the at least one processor to execute the method for identifying vortex-induced vibrations of bridges.

[0022] Through this arrangement, the convenience of using the method for identifying vortex-induced vibration of a bridge is improved.

[0023] Compared with the prior art, the present invention has the following beneficial effects: 1. The present invention provides a method for identifying vortex-induced vibration of a bridge. The method is based on the time domain characteristics of vortex-induced vibration of the bridge and takes the signal envelope as the judgment basis. Through this method, the envelope of the time domain signal is picked up and the identification index is quantified to realize the monitoring and identification of vortex-induced vibration events of the bridge. Compared with the prior art, this method has the advantages of high index stability, good index quantification effect, simple bridge vortex-induced vibration analysis process, small calculation amount, strong real-time performance, visualization, etc. 2. The present invention provides a system for identifying vortex-induced vibration of a bridge. Through this configuration, the efficiency of identifying vortex-induced vibration of a bridge is improved; 3. The present invention provides a system device for identifying vortex-induced vibration of a bridge. Through this configuration, the convenience of using the method for identifying vortex-induced vibration of a bridge is improved. BRIEF DESCRIPTION OF THE DRAWINGS

[0024] Figure 1 This is a flow chart of the method for identifying vortex-induced vibration of a bridge according to Example 1 provided by the present invention.

[0025] Figure 2 It is a waveform diagram of a vortex-induced vibration signal containing redundant data in Example 1 provided by the present invention.

[0026] Figure 3 This is a waveform diagram of a random vibration signal containing redundant data in Example 1 provided by the present invention.

[0027] Figure 4 It is the vortex-induced vibration waveform diagram and its upper envelope waveform diagram after removing redundant data in Example 1 provided by the present invention.

[0028] Figure 5 It is a random vibration waveform diagram and an upper envelope waveform diagram after removing redundant data in Example 1 provided by the present invention.

[0029] Figure 6 It is a distribution diagram of the vortex-induced vibration signal and its upper envelope in polar coordinates of Example 1 provided by the present invention.

[0030] Figure 7 It is a distribution diagram of the random vibration signal and its upper envelope in polar coordinates of Example 1 provided by the present invention.

[0031] Figure 8 It is a comparison diagram of the vibration signal waveform and the identification index of Example 1 provided by the present invention. DETAILED DESCRIPTION

[0032] The present invention is further described in detail below in conjunction with specific embodiments. However, this should not be understood as the scope of the above subject matter of the present invention being limited to the following embodiments, and all technologies realized based on the content of the present invention belong to the scope of the present invention.

[0033] Unless otherwise specified, in the description of the specific embodiments of the present invention, the terms indicating the orientation or position relationship such as "up", "down", "left", "right", "center", "inside", "outside", etc. are all expressions based on the orientation or position relationship shown in the drawings, or are the orientation or position relationship when the invented product / equipment / device is usually used. These terms of orientation or position relationship are only for the convenience of describing the scheme of the present invention or simplifying the description in the specific embodiments, so as to facilitate the technicians to quickly understand the scheme, and do not indicate or imply that a specific device / component / element must have a specific orientation, or be constructed and operated in a specific position relationship, and therefore cannot be understood as a limitation on the present invention.

[0034] In addition, if the terms "horizontal", "vertical", "overhanging", "parallel" and the like appear, it does not mean that the corresponding devices / components / elements are required to be absolutely horizontal or vertical or overhanging or parallel, but may be slightly tilted or have deviations. For example, "horizontal" only means that its direction is more horizontal than "vertical", and does not mean that the structure must be completely horizontal, but may be slightly tilted. Alternatively, it can be simplified to mean that the corresponding devices / components / elements are set in directions such as "horizontal", "vertical", "overhanging", "parallel", etc., and can have an error / deviation of ±10% relative to the corresponding direction setting, more preferably an error / deviation within ±8%, more preferably an error / deviation within ±6%, more preferably an error / deviation within ±5%, and more preferably an error / deviation within ±4%. As long as the corresponding device / component / element is within the error / deviation range, it can still achieve its role in the scheme of the present invention.

[0035] In addition, the expressions “first”, “second”, “third”, etc., which appear in the terms, are merely used to distinguish the description of the same or similar components and should not be understood as emphasizing or implying the relative importance of specific components.

[0036] In addition, in the description of the embodiments of the present invention, "several", "plurality" and "a number" represent at least 2. It can be any number such as 2, 3, 4, 5, 6, 7, 8, 9, and even more than 9.

[0037] In addition, in the description of the technical solution of the present invention, unless otherwise clearly specified / defined / restricted, the terms "set", "install", "connect", "connected", "provided with", "laid", and "arranged" should be understood in a broad sense, for example, it can be a fixed connection, a detachable connection, or an integral connection, and can be welding, riveting, bolting, threading, and other commonly used connection means in the field. This connection can be a mechanical connection, an electrical connection, or a communication connection; it can be a direct connection, or an indirect connection through an intermediate medium, and it can be the internal connection of two elements.

[0038] Example 1 like Figure 1-Figure 8 As shown, this embodiment provides a method for identifying vortex-induced vibration of a bridge, comprising the following steps: S1: Time-course monitoring signal data of bridge vibration De-averaging and high-pass filtering are performed. , where i represents the time-course monitoring signal data The number of the data point in is a constant; when the bridge undergoes vortex-induced vibration, the waveform of the time-history monitoring signal data is as follows Figure 2 As shown in the figure, when the bridge vibrates randomly, the waveform of the time-history monitoring signal data is as follows: Figure 3 shown.

[0039] S2: OK The 0 phase point and phase point, and remove the redundant data before the 0 phase point and Redundant data after the phase point to obtain a complete cycle data sequence , when the bridge is subjected to vortex-induced vibration, ,in, represents the number of complete cycles, A sequence of data representing a complete cycle The number of the data point in is a constant; S3: According to Calculated Amplitude sequence of the upper envelope , ;like Figure 4 As shown, Figure 2 The upper envelope waveform of the vortex-induced vibration waveform can be obtained by processing the time-course monitoring signal data in the above steps. Figure 6 for Figure 4 The distribution diagram of the waveform in polar coordinates; Figure 5 As shown, Figure 3 The upper envelope waveform of the random vibration waveform can be obtained by processing the time-course monitoring signal data in the above steps. Figure 7 for Figure 5 Distribution of the waveform in polar coordinates.

[0040] S4: According to The identification index of bridge vortex-induced vibration is calculated ,set up is the preset recognition threshold. Greater than or equal to When , it is identified that the bridge is experiencing vortex-induced vibration.

[0041] When the bridge vibrates, its vibration signal waveform and identification indicators The comparison chart is as follows Figure 8 As shown in the figure, it can be seen that by picking up the upper envelope of the time domain signal and quantifying the identification index, and comparing the upper envelope of the signal with the identification index, the process of monitoring and identifying the vortex-induced vibration event of the bridge can be more clearly displayed.

[0042] By monitoring the bridge vibration signal data over time De-averaging and high-pass filtering are performed to reduce the overall signal baseline offset caused by the DC component and the trend drift caused by low-frequency noise, thereby improving subsequent recognition accuracy.

[0043] This method is based on the time domain characteristics of bridge vortex-induced vibration, takes the signal upper envelope as the judgment basis, and realizes the monitoring and identification of bridge vortex-induced vibration events by picking up the time domain signal upper envelope and quantifying the identification index. Compared with the existing technology, this method has the advantages of high indicator stability, good indicator quantification effect, simple bridge vortex-induced vibration analysis process, small calculation amount, strong real-time performance, and visualization.

[0044] The step S2 comprises: S21: Establishing a rectangular coordinate system ,Will Each data point in is represented by a rectangular coordinate system and is obtained In the rectangular coordinate system The sequence of intersection points of the axes , The expression is: ; S22: Settings The curve with a positive slope is The sequence of axis intersections is ,set up The curve with a negative slope is Sequence of axis intersections ,and ,in: The expression is: , The expression is: ; S23: Settings The sequence of all points on a curve with a positive slope is ,set up The sequence of all points on a curve with a negative slope is ,and , and Respectively expressed as: in, Representation sequence The minimum value in Representation sequence The maximum value in Representation sequence The minimum value in Representation sequence The maximum value in .

[0045] In step S3, the local maximum method, interpolation method, polynomial fitting method, Hilbert transform method, empirical mode decomposition method, or wavelet transform method are used to calculate Amplitude sequence of the upper envelope .

[0046] Calculated using the Hilbert transform method Amplitude sequence of the upper envelope , the step S3 comprises: S31: Yes Performing Hilbert transform, we get ,in, Representatives The sequence after Hilbert transformation, is a constant; S32: According to and Get the amplitude sequence of the upper envelope , It is expressed as: .

[0047] In step S4, the identification index pass The minimum value of The ratio of the average values ​​of is calculated. It is expressed as: ; in, express The minimum value in express The average value of .

[0048] Identify indicators pass The minimum value of This method improves the quantization effect of the stability of the signal amplitude, and the amount of calculation is small. Real-time calculation occupies less computing resources, which improves work efficiency.

[0049] The identification threshold Including the first recognition threshold , the second recognition threshold and the third recognition threshold ,and , when Greater than or equal to and Less than When , it is identified that the bridge is in the embryonic stage of vortex-induced vibration; when Greater than or equal to and Less than When , it is identified that the bridge is in the development stage of vortex-induced vibration; when Greater than or equal to When , it is identified that the bridge is in the stable stage of vortex-induced vibration.

[0050] Said The value of is 0.1, The value of is 0.5, The value of is 0.8.

[0051] Said The duration is 1min.

[0052] Since the vortex-induced vibration of a bridge is a developing process, the data stability is higher within a time series of 1 minute. Therefore, selecting data with a time series of 1 minute for identification can improve the accuracy of judging the vortex-induced vibration state of the bridge.

[0053] Example 2 This embodiment provides a system for identifying vortex-induced vibration of a bridge, which adopts the method for identifying vortex-induced vibration of a bridge, including a bridge vibration signal processing module, a redundant data elimination module, an upper envelope calculation module, a bridge vortex-induced vibration identification index calculation module and a bridge vortex-induced vibration judgment module, wherein the bridge vibration signal processing module is used to perform mean removal and high-pass filtering on the time-history monitoring signal data of bridge vibration; The redundant data elimination module is used to eliminate redundant data to obtain complete cycle bridge vibration time history monitoring signal data; The upper envelope calculation module is used to calculate the upper envelope of the bridge vibration time history monitoring signal data of a complete cycle; The bridge vortex-induced vibration identification index calculation module is used to calculate the identification index of the bridge vortex-induced vibration; The bridge vortex-induced vibration judgment module is used to compare the bridge vortex-induced vibration identification index with a preset threshold value and obtain a judgment result.

[0054] Through this setting, the efficiency of identifying vortex-induced vibration of bridges is improved.

[0055] Example 3 This embodiment provides a system device for identifying vortex-induced vibration of a bridge, including at least one processor and a memory communicatively connected to the at least one processor; the memory stores instructions executable by the at least one processor, and the instructions are executed by the at least one processor so that the at least one processor can execute the method for identifying vortex-induced vibration of a bridge.

[0056] Through this arrangement, the convenience of using the method for identifying vortex-induced vibration of a bridge is improved.

[0057] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions and improvements made within the spirit and principles of the present invention should be included in the protection scope of the present invention.

Claims

1. A method for identifying vortex-induced vibration of a bridge, characterized in that: The following steps are involved: S1: Time-course monitoring signal data of bridge vibration De-averaging and high-pass filtering are performed. , where i represents the time-course monitoring signal data The number of the data point in is a constant; S2: OK The 0 phase point and phase point, and remove the redundant data before the 0 phase point and Redundant data after the phase point to obtain a complete cycle data sequence , ,in, represents the number of complete cycles, A sequence of data representing a complete cycle The number of the data point in is a constant; S3: According to Calculated Amplitude sequence of the upper envelope , ; S4: According to The identification index of bridge vortex-induced vibration is calculated ,set up is the preset recognition threshold. Greater than or equal to When , it is identified that the bridge is experiencing vortex-induced vibration.

2. A method for identifying vortex-induced vibration of a bridge according to claim 1, characterized in that: The step S2 comprises: S21: Establishing a rectangular coordinate system ,Will Each data point in is represented by a rectangular coordinate system and is obtained In the rectangular coordinate system The sequence of intersections of the axes , The expression is: ； S22: Settings The curve with a positive slope is The sequence of axis intersections is ,set up The curve with a negative slope is Sequence of axis intersections ,and ,in: The expression is: , The expression is: ; S23: Settings The sequence of all points on a curve with a positive slope is ,set up The sequence of all points on a curve with a negative slope is ,and , and Respectively expressed as: in, Representation sequence The minimum value in Representation sequence The maximum value in Representation sequence The minimum value in Representation sequence The maximum value in .

3. A method for identifying vortex-induced vibration of a bridge according to claim 1, characterized in that: In step S3, the local maximum method, interpolation method, polynomial fitting method, Hilbert transform method, empirical mode decomposition method, or wavelet transform method are used to calculate Amplitude sequence of the upper envelope .

4. A method for identifying vortex-induced vibration of a bridge according to claim 3, characterized in that: Calculated using the Hilbert transform method Amplitude sequence of the upper envelope , the step S3 comprises: S31: Yes Performing Hilbert transform, we get ,in, Representatives The sequence after Hilbert transformation, is a constant; S32: According to and Get the amplitude sequence of the upper envelope , It is expressed as: .

5. The method for identifying vortex-induced vibration of a bridge according to claim 1, characterized in that: In step S4, the identification index pass The minimum value of The ratio of the average values ​​of is calculated. It is expressed as: ； in, express The minimum value in express The average value of .

6. A method for identifying vortex-induced vibration of a bridge according to any one of claims 1 to 5, characterized in that: The identification threshold Including the first recognition threshold , the second recognition threshold and the third recognition threshold ,and , when Greater than or equal to and Less than When , it is identified that the bridge is in the embryonic stage of vortex-induced vibration; when Greater than or equal to and Less than When , it is identified that the bridge is in the development stage of vortex-induced vibration; when Greater than or equal to When , it is identified that the bridge is in the stable stage of vortex-induced vibration.

7. A method for identifying vortex-induced vibration of a bridge according to claim 6, characterized in that: Said The value of is 0.05-0.15, The value of is 0.45-0.55, The value is 0.75-0.

85.

8. The method for identifying vortex-induced vibration of a bridge according to claim 1, characterized in that: Said The duration is 0.9-1.1min.

9. A system for identifying vortex-induced vibration of a bridge, characterized in that: A method for identifying vortex-induced vibration of a bridge as claimed in any one of claims 1 to 8, comprising a bridge vibration signal processing module, a redundant data elimination module, an upper envelope calculation module, a bridge vortex-induced vibration identification index calculation module and a bridge vortex-induced vibration judgment module, wherein the bridge vibration signal processing module is used to perform de-averaging and high-pass filtering on the time-history monitoring signal data of the bridge vibration; The redundant data elimination module is used to eliminate redundant data to obtain complete cycle bridge vibration time history monitoring signal data; The upper envelope calculation module is used to calculate the upper envelope of the bridge vibration time history monitoring signal data of a complete cycle; The bridge vortex-induced vibration identification index calculation module is used to calculate the identification index of the bridge vortex-induced vibration; The bridge vortex-induced vibration judgment module is used to compare the bridge vortex-induced vibration identification index with a preset threshold value and obtain a judgment result.

10. A system device for identifying vortex-induced vibration of a bridge, characterized in that: comprising at least one processor, and a memory communicatively connected to the at least one processor; The memory stores instructions that can be executed by at least one processor, and the instructions are executed by at least one processor so that the at least one processor can execute the method for identifying vortex-induced vibration of a bridge as described in any one of claims 1-8.

Images