CN110470422B - 一种基于拉索振动频谱分析的基频优化方法 - Google Patents
一种基于拉索振动频谱分析的基频优化方法 Download PDFInfo
- Publication number
- CN110470422B CN110470422B CN201910778668.1A CN201910778668A CN110470422B CN 110470422 B CN110470422 B CN 110470422B CN 201910778668 A CN201910778668 A CN 201910778668A CN 110470422 B CN110470422 B CN 110470422B
- Authority
- CN
- China
- Prior art keywords
- frequency
- fundamental frequency
- cable
- theoretical
- inhaul cable
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000004458 analytical method Methods 0.000 title claims abstract description 19
- 238000005457 optimization Methods 0.000 title claims abstract description 14
- 238000001845 vibrational spectrum Methods 0.000 title claims abstract description 13
- 238000001228 spectrum Methods 0.000 claims abstract description 22
- 238000010183 spectrum analysis Methods 0.000 claims abstract description 17
- 238000005259 measurement Methods 0.000 claims abstract description 8
- 238000004364 calculation method Methods 0.000 claims abstract description 7
- 238000005070 sampling Methods 0.000 claims description 5
- 230000003595 spectral effect Effects 0.000 claims description 3
- 238000000547 structure data Methods 0.000 abstract description 2
- 229910000831 Steel Inorganic materials 0.000 description 4
- 239000010959 steel Substances 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
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/04—Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes for measuring tension in flexible members, e.g. ropes, cables, wires, threads, belts or bands
-
- 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/04—Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes for measuring tension in flexible members, e.g. ropes, cables, wires, threads, belts or bands
- G01L5/042—Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes for measuring tension in flexible members, e.g. ropes, cables, wires, threads, belts or bands by measuring vibrational characteristics of the flexible member
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明公开了桥梁结构数据分析技术领域的一种基于拉索振动频谱分析的基频优化方法，旨在解决现有技术中拉索基频往往由于能量较低而难以识别，导致索力计算误差较大的技术问题。所述方法包括如下步骤：获取拉索的理论基频；采集拉索振动原始数据；对拉索振动原始数据进行功率谱分析，获取拉索振动功率谱频域数据；基于理论基频对拉索振动功率谱频域数据进行频谱分析，获取拉索的实测基频。
Description
技术领域
本发明涉及一种基于拉索振动频谱分析的基频优化方法，属于桥梁结构数据分析技术领域。
背景技术
拉索是桥梁结构的重要受力构件，拉索索力的大小，直接关系到桥梁整体结构的受力及变形，因而精确获取索力值对于桥梁安全评估具有重大意义。现有技术中，基频法是拉索频率法计算索力的常用方法，然而在实际测量中，由于各种噪声以及设备精度的影响，使得拉索的基频（即一阶频率）往往由于能量较低而难以识别，导致索力计算误差较大。
发明内容
针对现有技术的不足，本发明的目的在于提供一种基于拉索振动频谱分析的基频优化方法，以解决现有技术中拉索基频往往由于能量较低而难以识别，导致索力计算误差较大的技术问题。
为解决上述技术问题，本发明所采用的技术方案是：
一种基于拉索振动频谱分析的基频优化方法，包括如下步骤：
获取拉索的理论基频；
采集拉索振动原始数据；
对拉索振动原始数据进行功率谱分析，获取拉索振动功率谱频域数据；
基于理论基频对拉索振动功率谱频域数据进行频谱分析，获取拉索的实测基频。
进一步地，获取拉索的理论基频，包括：
检测获取成桥索力和索力系数；
根据成桥索力和索力系数，计算获取拉索的理论基频。
进一步地，根据成桥索力和索力系数，计算获取拉索的理论基频，包括：基于弦理论频率法计算获取拉索的理论基频，其计算公式如下：
式中，
为拉索的理论基频，
为成桥索力，K为索力系数。
进一步地，所述拉索振动原始数据的采集数量大于102×F，式中，F为拉索振动原始数据的采样频率。
进一步地，功率谱分析的频率范围为［0，F/2］Hz，式中，F为拉索振动原始数据的采样频率。
进一步地，基于理论基频对拉索振动功率谱频域数据进行频率分析，获取拉索的实测基频，包括：
提取理论基频所对应的预设阶次频率，所述预设阶次频率不少于两个；
基于预设阶次频率和理论基频在拉索振动功率谱频域数据中截取频率区间，所述频率区间不少于两个；
对每个频率区间分别进行拉索频谱分析，提取频率区间内功率谱密度峰值所对应的频率，作为峰值频率；
各峰值频率分别除以其所处频率区间对应预设阶次频率的阶次，构成基频数组；
提取基频数组中最接近理论基频的元素，作为拉索的实测基频。
进一步地，所述预设阶次频率为理论基频所对应的前M阶频率。
进一步地，M＝7。
进一步地，所述预设阶次频率，其表达式如下：
式中，
为预设阶次频率，
为理论基频；
基于预设阶次频率和理论基频在拉索振动功率谱频域数据中截取的频率区间，其表达式如下：
进一步地，所述基频数组，其表达式如下：
式中，
为基频数组，
为峰值频率；
所述实测基频，其表达式如下：
其中，
；
式中，
为拉索的实测基频，N为实测基频对应的阶次，
为理论基频，
为基频数组中元素与理论基频差值的绝对值所构成的数组。
与现有技术相比，本发明所达到的有益效果：本发明方法通过拉索理论基频划分实测拉索振动数据功率谱频域区间，计算拉索前7阶振动频率，然后将各阶频率除以对应阶次得到拉索基频数组，选择基频数组中最接近理论基频的数值作为实测基频，能够有效解决基频法索力计算过程中，一阶频率峰值难以识别及高阶频率约束条件较多所带来的误差，因而更为高效准确，同时有利于集成算法智能化获取海量拉索振动数据的基频信息，对于索力检测及监测具有重要意义，可以广泛应用于工程实践中。
附图说明
图1是本发明方法实施例中的拉索振动原始数据时程图；
图2是本发明方法实施例中的测点SL-18-01功率谱图。
具体实施方式
下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案，而不能以此来限制本发明的保护范围。
本发明具体实施方式提供了一种基于拉索振动频谱分析的基频优化方法，本实施例以郑州黄河公铁两用大桥2017年07月01日凌晨00:00:00-00:00:10共计10分钟内的拉索振动监测数据为例，解析本发明方法的各个步骤。
郑州黄河公铁两用大桥为石郑铁路客运专线及河南省规划的中原黄河公路大桥跨越黄河的共用桥梁，桥位距下游京珠高速公路黄河大桥约6km。规划的中原黄河公路大桥北起河南省新乡市境内原阳县原武镇阎庄，接现国道G107线，向南经原阳县原武镇，跨越黄河，与新建的G107辅道相接。大桥主桥全长1684m，共两联，第一联为120＋5x168＋120m的六塔连续钢桁结合梁斜拉桥。上层公路桥面宽32.5m，下层铁路桥面为双线客运专线，线间距7m。主桁为三角形桁式，横向三片桁布置，中桁垂直，边桁倾斜。钢桁梁上弦杆与混凝土桥面板结合形成公路结合桥面，下层铁路桥面为正交异性整体钢桥面板。主塔为钢箱结构，纵桥向为“人”字形布置，设置在中桁。公路混凝土桥面板纵横向均为全预应力结构。
拉索振动测点位于5#塔两侧边索及第三根索，选取边索测点SL-18-01监测数据，采样频率F为100Hz，采集时间为10min，在该采集时间内所采集到的数据量为
。
本发明方法包括如下步骤：
步骤一，获取拉索振动原始数据、成桥索力
及索力系数K。
所述拉索振动原始数据为通过边索测点SL-18-01监测获取的索力测试原始数据data，该拉索振动原始数据的时序图如图1所示。为确保精确获得实测基频，一般要求索力测试原始数据data的采集数量
，本实施例中的数据采集数量显然满足此要求。
步骤二，由成桥索力
及索力系数K计算拉索的理论基频。
基于弦理论频率法，拉索的理论基频为：
式中，
为拉索的理论基频。
步骤三，基于理论基频截取实测拉索振动功率谱频域数据的频率区间，并计算各个频域区间内功率谱密度峰值对应的频率。
首先，对索力测试原始数据data进行功率谱分析，获取拉索振动功率谱频域数据。功率谱分析的频率范围通常为［0，F/2］Hz。本实施例中，由于F＝100Hz，则功率谱分析的频率范围为［0，50］Hz，功率谱分析获取的测点SL-18-01功率谱如图2所示；
然后，提取理论基频所对应的前7阶频率，即：
式中，
为理论基频所对应的前7阶频率；
接着，利用前7阶频率在拉索振动功率谱频域数据中截取7段用于拉索频谱分析的频率区间，每段频率区间记为：
最后，对每个频率区间进行拉索频谱分析，提取频率区间内功率谱密度峰值所对应的频率，即峰值频率；基于频率区间和峰值频率建立频率分析表，具体如表1所示。所述峰值频率记为：
。
表
：频率分析表
本实施例中，之所以提取理论基频所对应的前7阶频率，是由于拉索前7阶频率相对稳定，能量集中，峰值明显，同时，拉索高阶频率的精度容易受到拉索刚度、约束条件、垂度等因素影响，从而达到准确获得拉索各阶频率，提高计算效率，保证拉索频率识别精度。
步骤四，7段频域区间中的峰值频率分别除以其所处频率区间对应频率阶次，构成基频数组，其表达式如下：
式中，
为基频数组。
步骤五，提取基频数组中最接近理论基频的元素，作为拉索的实测基频。
先将基频数组
中的元素分别与理论基频
做差，再取绝对值，记为：
，式中，
为基频数组中元素与理论基频差值的绝对值所构成的数组；实测基频的表达式，记为：，式中，
为拉索的实测基频，N为实测基频对应的阶次。可以看出，由理论基频与其所对应的第五阶频率，在拉索振动功率谱频域数据中截取的频率区间内功率谱密度峰值所对应的基频，与理论基频的差值为0，因而最为接近。则可以得出，拉索的实测基频为
Hz。
本发明方法通过拉索理论基频划分实测拉索振动数据功率谱频域区间，计算拉索前7阶振动频率，然后将各阶频率除以对应阶次得到拉索基频数组，选择基频数组中最接近理论基频的数值作为实测基频，该方法有效解决基频法索力计算过程中，一阶频率峰值难以识别及高阶频率约束条件较多所带来的误差，因而更为高效准确，同时有利于集成算法智能化获取海量拉索振动数据的基频信息，对于索力检测及监测具有重要意义，可以广泛应用于工程实践中。
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和变形，这些改进和变形也应视为本发明的保护范围。
Claims (7)
1.一种基于拉索振动频谱分析的基频优化方法，其特征是，包括如下步骤：
获取拉索的理论基频；
采集拉索振动原始数据；
对拉索振动原始数据进行功率谱分析，获取拉索振动功率谱频域数据；
基于理论基频对拉索振动功率谱频域数据进行频谱分析，获取拉索的实测基频；
基于理论基频对拉索振动功率谱频域数据进行频率分析，获取拉索的实测基频，包括：
提取理论基频所对应的预设阶次频率，所述预设阶次频率不少于两个；
基于预设阶次频率和理论基频在拉索振动功率谱频域数据中截取频率区间，所述频率区间不少于两个；
对每个频率区间分别进行拉索频谱分析，提取频率区间内功率谱密度峰值所对应的频率，作为峰值频率；
各峰值频率分别除以其所处频率区间对应预设阶次频率的阶次，构成基频数组；
提取基频数组中最接近理论基频的元素，作为拉索的实测基频；
所述预设阶次频率，其表达式如下：
式中，
为预设阶次频率，
为理论基频；
基于预设阶次频率和理论基频在拉索振动功率谱频域数据中截取的频率区间，其表达式如下：
所述基频数组，其表达式如下：
式中，
为基频数组，
为峰值频率；
所述实测基频，其表达式如下：
其中，
；
式中，
为拉索的实测基频，N为实测基频对应的阶次，
为理论基频，
为基频数组中元素与理论基频差值的绝对值所构成的数组。
2.根据权利要求1所述的基于拉索振动频谱分析的基频优化方法，其特征是，获取拉索的理论基频，包括：
检测获取成桥索力和索力系数；
根据成桥索力和索力系数，计算获取拉索的理论基频。
3.根据权利要求2所述的基于拉索振动频谱分析的基频优化方法，其特征是，根据成桥索力和索力系数，计算获取拉索的理论基频，包括：基于弦理论频率法计算获取拉索的理论基频，其计算公式如下：
式中，
为拉索的理论基频，
为成桥索力，K为索力系数。
4.根据权利要求1所述的基于拉索振动频谱分析的基频优化方法，其特征是，所述拉索振动原始数据的采集数量大于102×F，式中，F为拉索振动原始数据的采样频率。
5.根据权利要求1所述的基于拉索振动频谱分析的基频优化方法，其特征是，功率谱分析的频率范围为［0，F/2］Hz，式中，F为拉索振动原始数据的采样频率。
6.根据权利要求1所述的基于拉索振动频谱分析的基频优化方法，其特征是，所述预设阶次频率为理论基频所对应的前M阶频率。
7.根据权利要求6所述的基于拉索振动频谱分析的基频优化方法，其特征是，M＝7。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910778668.1A CN110470422B (zh) | 2019-08-22 | 2019-08-22 | 一种基于拉索振动频谱分析的基频优化方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910778668.1A CN110470422B (zh) | 2019-08-22 | 2019-08-22 | 一种基于拉索振动频谱分析的基频优化方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110470422A CN110470422A (zh) | 2019-11-19 |
| CN110470422B true CN110470422B (zh) | 2020-11-17 |
Family
ID=68512851
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910778668.1A Active CN110470422B (zh) | 2019-08-22 | 2019-08-22 | 一种基于拉索振动频谱分析的基频优化方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110470422B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114544062B (zh) * | 2022-01-26 | 2024-07-23 | 四川升拓检测技术股份有限公司 | 一种模拟测试悬索张拉力的系统和方法 |
| CN114459596B (zh) * | 2022-02-11 | 2022-08-02 | 中冶建筑研究总院有限公司 | 基于动力测试的拉索频率综合自验证分析方法 |
| CN114398712B (zh) * | 2022-03-25 | 2022-07-01 | 西南交通大学 | 斜拉索实时索力的计算方法、装置、设备及可读存储介质 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101201282A (zh) * | 2007-12-20 | 2008-06-18 | 宁波大学 | 一种用于斜拉桥索力检测的基频识别方法 |
| CN101368860A (zh) * | 2008-09-12 | 2009-02-18 | 江苏工业学院 | 频率法检测斜拉桥斜拉索索力中fft数据的修正方法 |
| CN101586997A (zh) * | 2009-06-26 | 2009-11-25 | 贵州师范大学 | 一种拉索振动基频的计算方法 |
| CN102519651A (zh) * | 2011-12-13 | 2012-06-27 | 清华大学 | 一种振动方法测试斜拉桥索力时确定拉索基本频率的方法 |
| RU2613484C2 (ru) * | 2015-08-17 | 2017-03-16 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Сибирский государственный университет путей сообщения" (СГУПС) | Способ определения усилия натяжения вантового элемента моста |
| CN106840497A (zh) * | 2016-12-31 | 2017-06-13 | 嘉兴市纳杰微电子技术有限公司 | 桥梁拉索拉力识别方法 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101504324B (zh) * | 2009-02-26 | 2010-09-15 | 南京工业大学 | 拉索运营状态智能判别方法及系统 |
| CN105698942B (zh) * | 2016-03-11 | 2019-06-21 | 中国科学技术大学 | 一种基于光子轨道角动量的高精度波长计 |
| CN106200542A (zh) * | 2016-08-11 | 2016-12-07 | 北京大唐恒通机械输送技术有限公司 | 干渣机远程在线监测分析系统及方法 |
| CN106886628B (zh) * | 2017-01-12 | 2019-08-13 | 西安电子科技大学 | 一种基于索张力不确定性的平面薄膜天线薄膜形状确定方法 |
-
2019
- 2019-08-22 CN CN201910778668.1A patent/CN110470422B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101201282A (zh) * | 2007-12-20 | 2008-06-18 | 宁波大学 | 一种用于斜拉桥索力检测的基频识别方法 |
| CN101368860A (zh) * | 2008-09-12 | 2009-02-18 | 江苏工业学院 | 频率法检测斜拉桥斜拉索索力中fft数据的修正方法 |
| CN101586997A (zh) * | 2009-06-26 | 2009-11-25 | 贵州师范大学 | 一种拉索振动基频的计算方法 |
| CN102519651A (zh) * | 2011-12-13 | 2012-06-27 | 清华大学 | 一种振动方法测试斜拉桥索力时确定拉索基本频率的方法 |
| RU2613484C2 (ru) * | 2015-08-17 | 2017-03-16 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Сибирский государственный университет путей сообщения" (СГУПС) | Способ определения усилия натяжения вантового элемента моста |
| CN106840497A (zh) * | 2016-12-31 | 2017-06-13 | 嘉兴市纳杰微电子技术有限公司 | 桥梁拉索拉力识别方法 |
Non-Patent Citations (2)
| Title |
|---|
| 基于振动频率法的两端固支拉索;甘泉等;《水利与建筑工程学报》;20180228;第16卷(第1期);59-65 * |
| 振动频率法测定索力原理及其影响因素分析;刘记海等;《山西建筑》;20110930;第37卷(第25期);171-172 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110470422A (zh) | 2019-11-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110470422B (zh) | 一种基于拉索振动频谱分析的基频优化方法 | |
| Fenerci et al. | Strong wind characteristics and dynamic response of a long-span suspension bridge during a storm | |
| CN101532917B (zh) | 一种桥梁承载能力快速荷载试验方法 | |
| WO2021036751A1 (zh) | 支座反力影响线曲率的连续梁损伤识别方法 | |
| CN102253311B (zh) | 一种基于实时线路参数测量的双端测距方法 | |
| CN103335589B (zh) | 无缝线路钢轨温度应力传感节点的标定方法 | |
| CN101762347B (zh) | 一种用半波法测量多跨钢拉索索力的方法 | |
| CN101865683A (zh) | 基于应变计的梁式桥位移动态测量系统及测量方法 | |
| CN209670980U (zh) | 一种隧道盾构管片现场监测系统 | |
| CN110285909B (zh) | 基于同步压缩变换的索承桥梁瞬时索力计算方法 | |
| CN107476204A (zh) | 一种高速铁路连续梁的线性控制方法 | |
| CN105716534B (zh) | 一种基于长标距光纤光栅传感器的系杆拱桥挠度识别方法 | |
| CN102506986A (zh) | 自立式铁塔及大跨越输电塔模态和振动测试系统及其方法 | |
| CN105320596A (zh) | 一种基于倾角仪的桥梁挠度测试方法及其系统 | |
| CN106482917A (zh) | 一种斜拉桥主梁动态挠度的检测方法 | |
| CN104167076B (zh) | 一种覆冰输电线路薄弱环节预警方法 | |
| CN109855770B (zh) | 一种基于钢轨横向加速度功率谱密度检测温度力的方法 | |
| CN108318126B (zh) | 一种基于钢轨竖向振动特性检测轨下垫板刚度的方法 | |
| CN110360984A (zh) | 一种地表沉降的大范围分布式监测系统及方法 | |
| CN108108568A (zh) | 一种桥梁结构在线安全监测中剔除温度对低频采样指标影响的实时方法 | |
| CN112883480B (zh) | 一种利用倾角与位移联合测试的桥梁荷载试验评估方法 | |
| CN102937492A (zh) | 预应力混凝土桥梁绝对应力监测方法 | |
| CN109992815B (zh) | 一种基于关联性雨流计数的支座累计位移计算方法 | |
| CN111982057A (zh) | 一种主桥钢梁合龙口间距监测方法 | |
| CN112781656A (zh) | 下穿高铁隧道施工时对相交既有隧道安全监测方法 |
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 | ||
| CB02 | Change of applicant information |
Address after: 210061 No. 8 Paneng Road, Jiangbei New District, Nanjing City, Jiangsu Province Applicant after: China Railway Bridge and Tunnel Technology Co.,Ltd. Address before: 210061 No. 8 Paneng Road, Jiangbei New District, Nanjing City, Jiangsu Province Applicant before: CHINA RAILWAY MAJOR BRIDGE (NANJING) BRIDGE AND RUNNEL INSPECT & RETROFIT Co.,Ltd. |
|
| CB02 | Change of applicant information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |