CN114777910A - Cable multimode vortex-induced vibration monitoring method - Google Patents
Cable multimode vortex-induced vibration monitoring method Download PDFInfo
- Publication number
- CN114777910A CN114777910A CN202210342564.8A CN202210342564A CN114777910A CN 114777910 A CN114777910 A CN 114777910A CN 202210342564 A CN202210342564 A CN 202210342564A CN 114777910 A CN114777910 A CN 114777910A
- Authority
- CN
- China
- Prior art keywords
- vortex
- induced vibration
- modal
- cable
- time
- 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
- 238000000034 method Methods 0.000 title claims abstract description 19
- 238000012544 monitoring process Methods 0.000 title claims abstract description 17
- 238000001228 spectrum Methods 0.000 claims abstract description 28
- 230000001133 acceleration Effects 0.000 claims abstract description 21
- 238000000513 principal component analysis Methods 0.000 claims abstract description 8
- 230000003595 spectral effect Effects 0.000 claims abstract description 3
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000010183 spectrum analysis Methods 0.000 claims description 3
- 230000002159 abnormal effect Effects 0.000 abstract description 5
- 238000010586 diagram Methods 0.000 description 2
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 229910052802 copper Inorganic materials 0.000 description 1
- 239000010949 copper Substances 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000004044 response Effects 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
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
The invention discloses a cable multi-mode vortex-induced vibration monitoring method which comprises the steps of obtaining a cable acceleration power spectrum, extracting a power spectrum and extracting a multi-order frequency fiAnd corresponding spectral amplitude AiI takes 1 to N; calculating the energy concentration coefficient Aj/A1J is 2 to N; when energy concentration coefficient A2/A1And A3/A1While satisfying 0.1>A2/A1>0.01 and 0.05>A3/A1>0.005, mark as multi-modal vortex-induced vibration onset time; when the multi-modal vortex-induced vibration cannot be simultaneously detected, marking the time as the ending time of the multi-modal vortex-induced vibration; for N-1 energy concentration coefficients A in the duration of multi-modal vortex-induced vibrationj/A1And (3) performing principal component analysis, selecting an energy concentration coefficient with the principal component contribution degree exceeding 95% and ranking at the front, recording the number of the energy concentration coefficients with the accumulated contribution degree exceeding 95% as n, wherein n is the number of the modes participating in the multi-mode vortex-induced vibration, and forming the modes with the number of n, namely the abnormal vibration mode of the stay cable. The invention can realize the time for the vortex-induced vibration of the inhaul cable,And recognizing vibration multimodality.
Description
Technical Field
The invention relates to a bridge structure performance monitoring method, in particular to a cable multi-mode vortex-induced vibration monitoring method.
Background
The stay cable is a key bearing component of the cable-stayed bridge, and the main beam and the vehicle load are all transmitted to the cable tower through the stay cable, so the normal service of the stay cable has important significance on the safety of the bridge. The service state of the bridge is complex, the bridge is continuously subjected to environmental loads such as wind and sunshine, and under the action of specific wind working conditions, the inhaul cable generates vortex-induced vibration, so that the inhaul cable is not only subjected to fatigue damage, but also can drive the bridge to generate abnormal vibration. Along with the increase of the length of the stay cable of the cable-stayed bridge, the vortex-induced vibration of the stay cable is obviously different from the single-mode vortex-induced vibration of the main beam and is expressed as multi-mode vortex-induced vibration.
In order to sense the vibration characteristic of the cable, an acceleration sensor is usually installed on the cable, a vibration frequency power spectrum is calculated according to acceleration time-course data, and the power spectrum multi-peak data represents the vibration characteristic of the cable. However, the prior art cannot solve the problem of multi-mode superposition in cable vortex excitation, so that the abnormal vibration state of the cable is difficult to perceive by using cable power spectrum data.
Disclosure of Invention
The purpose of the invention is as follows: the invention aims to provide a monitoring method capable of realizing cable vortex-induced vibration multi-mode recognition.
The technical scheme is as follows: the invention discloses a cable multi-modal vortex-induced vibration monitoring method, which comprises the following steps of:
(1) obtaining a power spectrum of the acceleration of the inhaul cable, and extracting the 1 st order to the Nth order frequency f of the power spectrumiAnd corresponding power spectral amplitudeAi，i＝1，2，3，......，N；
(2) Calculating the energy concentration coefficient Aj/A1，j＝2，3，......，N；
(3) When energy concentration coefficient A2/A1And A3/A1While satisfying 0.1 > A2/A1> 0.01 and 0.05 > A3/A1When the frequency is more than 0.005, marking as the starting time of the multi-modal vortex-induced vibration; when the two can not be met simultaneously, marking as the ending time of the multi-modal vortex-induced vibration; obtaining the multi-modal vortex-induced vibration duration time according to the vibration starting time and the vibration ending time;
(4) for N-1 energy concentration coefficients A in the duration of multi-modal vortex-induced vibrationj/A1And performing principal component analysis, arranging the principal component components from high to low according to the sizes of the principal component components, and when the principal component components accumulated from high to low exceed a set value, adding the 1 st order frequency, namely the vortex-induced vibration multi-order frequency mode, to the modal frequency value corresponding to the numerator of the energy concentration coefficient corresponding to the accumulated components, wherein the number n of the accumulated components is the multi-order modal number.
According to the method, through the vibration acceleration power spectrum of the stay cable, the starting time and the ending time of the multi-modal vortex-induced vibration of the stay cable are determined by utilizing the identification characteristic value and the threshold value, the specific gravity of the energy concentration coefficients participating in different modes is further analyzed, and the mode quantity is determined through principal component analysis, so that the abnormal vibration sensing of the stay cable is realized, reliable information is provided for a bridge management engineer to know the working state of a bridge in real time, and the warning and disaster assessment can be performed in time.
Furthermore, in the step (1), the measured cable acceleration time sequence data is segmented according to a set time interval, and power spectrum analysis calculation is carried out on each segment of acceleration data through a PSD (phase-sensitive detector), LPSD (low power digital) SD (secure digital) digital or DPSD (differential Power digital) algorithm to obtain a cable acceleration power spectrum. The discrete inhaul cable acceleration information is converted into power spectrum data, the power spectrum data are identified in the frequency domain range, negative effects caused by wild values, drift values and other data defects can be eliminated, and the monitoring precision is improved.
Further, the set time interval is determined according to the comprehensive calculation power and the real-time performance, and may be 1min, 10min or longer, and is preferably 1 min.
Further, the value of N is not more than 10. The multi-modal participation number of the inhaul cable vortex-induced vibration is generally not more than 10, so the method selects the first 10 main modes for analysis, namely N is 10.
Further, in the step (4), the set value is less than 1, preferably 95%, which can ensure that the monitoring method has strong robustness and accuracy.
Has the beneficial effects that: compared with the prior art, the invention has the following remarkable advantages: the method comprises the steps of monitoring the vibration acceleration power spectrum of the stay cable, determining the starting time and the ending time of the multi-mode vortex-induced vibration of the stay cable, evaluating the proportion of energy concentration coefficients corresponding to different modes by using principal component analysis, and selecting multi-order mode frequencies of the vortex-induced vibration, so that the number of the modes is determined.
Drawings
FIG. 1 is a block flow diagram of the present invention;
FIG. 2 is an acceleration power spectrum of the inhaul cable in an interval of 1 min;
FIG. 3 is energy concentration coefficient timing data;
FIG. 4 is a timing chart of event determination;
FIG. 5 shows the energy concentration coefficient contribution obtained by principal component analysis.
Detailed Description
The present invention will be described in further detail with reference to the accompanying drawings.
As shown in fig. 1, a method for monitoring multi-modal vortex-induced vibration of a cable comprises the following steps:
(1) acquiring the acceleration response of the inhaul cable, and calculating the acceleration power spectrum of the inhaul cable by taking 1min as a time interval;
(2) extracting the 1 st order frequency f of the power spectrum1And its corresponding power spectrum amplitude A1Extracting the second order frequency f of the power spectrum2And its corresponding power spectrum amplitude A2And so on to obtain the first ten frequencies f of the power spectrum1～f10And its corresponding amplitude A1～A10；
(3) Taking 1min as a time interval, calculating 9 energy setsMiddle coefficient A2/A1、A3/A1、A4/A1、A5/A1、A6/A1、A7/A1、A8/A1、A9/A1、A10/A1；
(4) When energy concentration coefficient A2/A1And A3/A1Simultaneously satisfies 0.1 > A2/A1> 0.01 and 0.05 > A3/A1When the frequency is more than 0.005, marking the starting time of the multi-mode vortex-induced vibration;
(5) when energy concentration coefficient A2/A1And A3/A1Cannot satisfy 0.1 > A simultaneously2/A1> 0.01 and 0.05 > A3/A1When the frequency is more than 0.005, marking the ending time of the multi-mode vortex-induced vibration;
(6) calculating the duration of the multi-modal vortex-induced vibration according to the starting time and the ending time of the multi-modal vortex-induced vibration;
(7) for 9 energy concentration coefficients A collected in the duration of multi-modal vortex-induced vibration2/A1、A3/A1、A4/A1、A5/A1、A6/A1、A7/A1、A8/A1、A9/A1、A10/A1And (3) performing principal component analysis, arranging the energy concentration coefficients from high to low according to the sizes of corresponding principal component components, accumulating the principal component components from high to low, and when the accumulated component exceeds 95%, adding 1 st order frequency to a modal frequency value corresponding to a numerator of the coefficient corresponding to the accumulated component to obtain a vortex-induced vibration multi-order frequency mode, wherein the number n of the accumulated component is the number of the multi-order modes.
The following description will be given to the concrete implementation process of the present invention by taking the cable acceleration of the dual-purpose cable-stayed bridge of holy road and railway as an example.
(1) Bridge structures health monitoring system that dual-purpose cable-stay bridge of copper tomb highway installed additional can collect the acceleration data of cable, and these acceleration data are data drive's basis. Taking 1min as a time interval, performing segmented processing on all acquired acceleration time-course data, and performing power spectrum analysis calculation on each segment of data, as shown in a power spectrum schematic diagram obtained by calculation within 1min in fig. 2.
(2) As shown in FIG. 2, the 1 st order frequency f of the power spectrum is extracted1And its corresponding power spectrum amplitude A1Extracting the second order frequency f of the power spectrum2And its corresponding power spectrum amplitude A2And so on to obtain the first ten frequencies f of the power spectrum1～f10And its corresponding amplitude A1～A10。
(3) Taking 1min as a time interval, calculating 9 energy concentration coefficients A2/A1、A3/A1、A4/A1、A5/A1、A6/A1、A7/A1、A8/A1、A9/A1、A10/A1FIG. 3 shows the energy concentration coefficient A2/A1Time course curve of (c).
(4) When energy concentration coefficient A2/A1And A3/A1Simultaneously satisfies 0.1 > A2/A1> 0.01 and 0.05 > A3/A1And when the time is more than 0.005, marking the starting time of the multi-modal vortex-induced vibration. As in fig. 4, the start moment is marked on the map at time 25 and the start token is marked.
(5) When energy concentration coefficient A2/A1And A3/A1Cannot satisfy 0.1 > A simultaneously2/A1> 0.01 and 0.05 > A3/A1And when the time is more than 0.005, marking the ending time of the multi-mode vortex-induced vibration, and further calculating the duration time of the multi-mode vortex-induced vibration. Referring to fig. 4, the end of the vortex-induced vibration is marked at the 100 th time, and the duration of the vibration is 75 min.
(6) For 9 energy concentration coefficients A collected in the multi-modal vortex-induced vibration duration2/A1、A3/A1、A4/A1、A5/A1、A6/A1、A7/A1、A8/A1、A9/A1、A10/A1Principal component analysis is performed, the number n of principal component components required for the principal component accumulation to exceed 95% is calculated, and the concentration coefficients are recorded. As shown in FIG. 5, A2/A1、A3/A1And A4/A1If the cumulative principal component of (1) exceeds 95%, n is 3.
(6) As shown in table 1, the number n of modes and the first-order principal vibration frequency, i.e., the abnormal vibration mode of the cable, are formed, and the number n of the accumulated principal components is the number of the combined modes participating in the multi-modal vortex-induced vibration, and the analysis result is output.
TABLE 1 finally outputted vortex-induced vibration mode Table
| Combined vortex induced |
2 order/3 order/4 order |
| Number of combined |
3 |
Claims (5)
1. A cable multi-mode vortex-induced vibration monitoring method is characterized by comprising the following steps: the method comprises the following steps:
(1) acquiring a power spectrum of the acceleration of the inhaul cable, and extracting the frequencies f from the 1 st order to the Nth order of the power spectrumiAnd corresponding power spectral amplitude Ai，i＝1，2，3，......，N；
(2) Calculating the energy concentration coefficient Aj/A1，j＝2，3，......，N；
(3) When energy concentration coefficient A2/A1And A3/A1While satisfying 0.1 > A2/A1> 0.01 and 0.05 > A3/A1When the frequency is more than 0.005, marking as the starting time of the multi-modal vortex-induced vibration; can not be full at the same timeWhen the time is sufficient, marking the time as the ending time of the multi-modal vortex-induced vibration; obtaining the multi-modal vortex-induced vibration duration according to the vibration starting time and the vibration ending time;
(4) for N-1 energy concentration coefficients A in the duration of multi-modal vortex-induced vibrationj/A1And performing principal component analysis, arranging the principal component components from high to low according to the sizes of the principal component components, and when the principal component components accumulated from high to low exceed a set value, adding the 1 st order frequency, namely the vortex-induced vibration multi-order frequency mode, to the modal frequency value corresponding to the numerator of the energy concentration coefficient corresponding to the accumulated components, wherein the number n of the accumulated components is the multi-order modal number.
2. The inhaul cable multi-modal vortex-induced vibration monitoring method according to claim 1, wherein: in the step (1), the measured acceleration time sequence data of the inhaul cable is segmented according to a set time interval, and power spectrum analysis calculation is carried out on each segment of acceleration data through a PSD (phase-sensitive detector), LPSD (low power digital) or DPSD (differential pulse-width digital) algorithm to obtain an acceleration power spectrum of the inhaul cable.
3. The inhaul cable multimodal vortex-induced vibration monitoring method according to claim 2, wherein: and the set time interval is determined according to the comprehensive calculation power and the real-time performance.
4. The inhaul cable multi-modal vortex-induced vibration monitoring method according to claim 1, wherein: the value of N is not more than 10.
5. The inhaul cable multi-modal vortex-induced vibration monitoring method according to claim 1, wherein: in the step (4), the set value is less than 1.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210342564.8A CN114777910A (en) | 2022-04-02 | 2022-04-02 | Cable multimode vortex-induced vibration monitoring method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210342564.8A CN114777910A (en) | 2022-04-02 | 2022-04-02 | Cable multimode vortex-induced vibration monitoring method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114777910A true CN114777910A (en) | 2022-07-22 |
Family
ID=82426787
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210342564.8A Pending CN114777910A (en) | 2022-04-02 | 2022-04-02 | Cable multimode vortex-induced vibration monitoring method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114777910A (en) |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101586997A (en) * | 2009-06-26 | 2009-11-25 | 贵州师范大学 | Method for calculating guy cable vibrating base frequency |
| JP2012242363A (en) * | 2011-05-24 | 2012-12-10 | East Japan Railway Co | Bridge pier soundness evaluation method using fine motion in strong wind |
| JP2015230206A (en) * | 2014-06-04 | 2015-12-21 | 公益財団法人鉄道総合技術研究所 | State monitoring system, information processing unit, state monitoring method, program, and recording medium |
| CN110688963A (en) * | 2019-09-30 | 2020-01-14 | 哈尔滨工业大学 | Clustering algorithm-based large-span bridge vortex-induced vibration automatic identification method |
| CN111832099A (en) * | 2020-05-28 | 2020-10-27 | 东南大学 | Bridge structure damage identification method based on FBG and finite element model correction |
| CN112362274A (en) * | 2020-10-10 | 2021-02-12 | 中铁大桥局集团有限公司 | Method and system for monitoring, early warning and evaluating vortex-induced vibration of long-span bridge in operation period |
| CN112818455A (en) * | 2021-02-22 | 2021-05-18 | 深圳市市政设计研究院有限公司 | Bridge structure response monitoring method and system |
| CN113343541A (en) * | 2021-07-08 | 2021-09-03 | 石家庄铁道大学 | Vortex-induced vibration early warning method, device and terminal for long and large span bridge |
| CN114169768A (en) * | 2021-12-09 | 2022-03-11 | 西南交通大学 | Method and system for evaluating comfort of traffic flow passing bridge based on bridge vortex-induced vibration |
-
2022
- 2022-04-02 CN CN202210342564.8A patent/CN114777910A/en active Pending
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101586997A (en) * | 2009-06-26 | 2009-11-25 | 贵州师范大学 | Method for calculating guy cable vibrating base frequency |
| JP2012242363A (en) * | 2011-05-24 | 2012-12-10 | East Japan Railway Co | Bridge pier soundness evaluation method using fine motion in strong wind |
| JP2015230206A (en) * | 2014-06-04 | 2015-12-21 | 公益財団法人鉄道総合技術研究所 | State monitoring system, information processing unit, state monitoring method, program, and recording medium |
| CN110688963A (en) * | 2019-09-30 | 2020-01-14 | 哈尔滨工业大学 | Clustering algorithm-based large-span bridge vortex-induced vibration automatic identification method |
| CN111832099A (en) * | 2020-05-28 | 2020-10-27 | 东南大学 | Bridge structure damage identification method based on FBG and finite element model correction |
| CN112362274A (en) * | 2020-10-10 | 2021-02-12 | 中铁大桥局集团有限公司 | Method and system for monitoring, early warning and evaluating vortex-induced vibration of long-span bridge in operation period |
| CN112818455A (en) * | 2021-02-22 | 2021-05-18 | 深圳市市政设计研究院有限公司 | Bridge structure response monitoring method and system |
| CN113343541A (en) * | 2021-07-08 | 2021-09-03 | 石家庄铁道大学 | Vortex-induced vibration early warning method, device and terminal for long and large span bridge |
| CN114169768A (en) * | 2021-12-09 | 2022-03-11 | 西南交通大学 | Method and system for evaluating comfort of traffic flow passing bridge based on bridge vortex-induced vibration |
Non-Patent Citations (3)
| Title |
|---|
| 张勇等: "大跨度悬索桥涡激振动动态监控预测", 振动与冲击, no. 08, pages 148 - 155 * |
| 曹素功等: "基于监测大数据的涡激振动影响因素分析", 桥梁建设, no. 001, pages 59 - 64 * |
| 辛大波等: "大跨度桥梁主梁风雨致涡激振动试验研究", 哈尔滨工程大学学报, no. 09, pages 70 - 74 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| WO2020244288A1 (en) | Method and apparatus for evaluating truck driving behaviour based on gps trajectory data | |
| CN105241660B (en) | High-speed rail large bridge performance test methods based on health monitoring data | |
| WO2017202139A1 (en) | Bridge damage identification method based on long-gauge-length strain influence envelope | |
| CN108515984B (en) | Wheel damage detection method and device | |
| CN108170925A (en) | A kind of bridge damnification method for early warning based on arma modeling | |
| CN104331595B (en) | The mobile principal component correlation coefficient process of bridge damnification early warning | |
| CN102175768A (en) | Method and device for detecting defects and failures of high-speed rail based on vibration signals | |
| Cheung et al. | The application of statistical pattern recognition methods for damage detection to field data | |
| CN105923014B (en) | A kind of track transition Amplitude Estimation method based on evidential reasoning rule | |
| CN101266190A (en) | Normal traffic flow stayd-cable bridge mode parametric measurement device and method | |
| JP6874858B2 (en) | Damage diagnostic equipment, damage diagnostic methods, and damage diagnostic programs | |
| CN106895900B (en) | Method based on the bridge dynamic strain identification bicycle number of axle | |
| CN105045983A (en) | Axle ageing analysis method of high speed train on the basis of axle temperature data | |
| CN104833535A (en) | Railway vehicle tire tread scratch detection method | |
| JP2015102329A (en) | Train information inferring method, and soundness evaluating method for bridges | |
| CN114997249B (en) | Monitoring method and system for bridge expansion joint device | |
| CN111400793A (en) | Movable intelligent rapid monitoring and bearing capacity evaluation algorithm for viaduct bridge | |
| WO2017107790A1 (en) | Method and apparatus for predicting road conditions using big data | |
| CN112926384A (en) | Automatic modal identification method based on power spectrum transfer ratio and support vector machine | |
| CN114777910A (en) | Cable multimode vortex-induced vibration monitoring method | |
| CN107782548A (en) | One kind is based on to track vehicle parts detecting system | |
| CN114323512B (en) | Heavy-load vehicle identification method and system | |
| CN114935450A (en) | Failure monitoring and alarming method for large-span cable-stayed bridge inhaul cable damper | |
| CN115290133A (en) | Method and system for monitoring track structure at joint of light rail platform | |
| CN109711075B (en) | Steel girder bridge life and reliability analysis method based on sudden load nonlinear theory |
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 |