CN110399683B - 基于频域幅值谱相似性滤波技术的桥梁冲击系数提取方法 - Google Patents
基于频域幅值谱相似性滤波技术的桥梁冲击系数提取方法 Download PDFInfo
- Publication number
- CN110399683B CN110399683B CN201910685464.3A CN201910685464A CN110399683B CN 110399683 B CN110399683 B CN 110399683B CN 201910685464 A CN201910685464 A CN 201910685464A CN 110399683 B CN110399683 B CN 110399683B
- Authority
- CN
- China
- Prior art keywords
- frequency
- bridge
- relative deviation
- low
- static
- 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
Images
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种基于频域幅值谱相似性滤波技术的桥梁冲击系数提取方法,针对准静态分量受测试条件干扰程度低、冲击系数检测不需要信号实时处理等工程特点,通过综合比对多条动位移时程曲线的频域幅值谱结果,可根据多条曲线幅值谱之间的相似性特点,自适应的确定滤波截止频率,有效解决了常用于动静分离的低通滤波法截止频率估计困难的问题,提高了实测冲击系数的计算准确性。本发明不增加检测试验的工作量,数据处理不依赖于人工判读,降低了试验人员主观因素对检测结果的影响,同时低通滤波截止频率的确定考虑到实测信号自身的特点,降低了由于滤波器设计偏差对试验结果的干扰,改善了目前试验因素影响冲击系数检测评定结果的技术现状。
Description
技术领域
本发明涉及一种冲击系数提取方法,具体涉及一种桥梁动力荷载试验检测中,基于桥梁动位移实测数据,用于桥梁冲击系数计算的试验数据处理方法。
背景技术
冲击系数是反映桥梁动力响应特性的重要指标,是桥梁设计中的重要参数。由冲击系数在各国规范中的发展历程可知,冲击系数的设计值是一个半经验半理论的取值,对冲击系数设计值的研究离不开现场试验数据的支持。同时,随着对桥梁检测与状态评估工作的重要性认识的逐渐提高,实测既有桥梁的冲击系数来评定桥梁的动力状态,已成为桥梁工作者的一项重要工作。因此,对冲击系数的实测检验方法进一步研究具有理论与实际的双重意义。
我国的桥梁设计规范(JTG D60-2015)中,冲击系数定义为:在汽车通过桥梁的效应时间曲线上,由最大静效应位置处,量测得到的最大动效应与最大静效应的比值。在工程应用中,冲击系数的测量多基于动位移测试数据,所以在时间-位移曲线上精确地获取最大动、静位移是计算实测冲击系数的关键。其中,最大动位移的获取较为容易,只需在曲线上量取位移最大值;但是,最大静位移是动位移曲线中的隐含信息,需要对实测动位移数据进行分析处理,实现动、静位移分量分离后才能够获得。
公路桥梁荷载试验规程(JTG/T J21-01-2015)中提供了一种通过直接的波峰-波谷平均法在位移时程曲线上计算最大静位移的方法。此方法对理想信号具有较好效果,也是目前普遍采用的实测冲击系数计算方法。但在实际操作中,由于噪声等其他不确定因素的干扰,实测时程曲线往往不存在完美的明显峰值和谷值,计算波峰、波谷的选择可能因人而异,导致无法做到对所有桥梁的稳定客观评价。
铁路桥梁检定规范(铁运函[2004]120号)中提供了一种通过车辆低速通过桥梁的跑车试验来获得准静态的试验方法。但此方法一方面需要进行额外的跑车试验,增加试验工作量;另一方面,Szurgott等人的研究表明,对于某些形式的桥梁,如跨径较大的混凝土梁桥,由于结构阻尼引起的变形滞后性,在相同车辆荷载作用下,高速跑车试验位移时程曲线中的准静态分量可能小于低速跑车试验获得的准静态位移曲线,导致以此方法测得的冲击系数偏小,甚至出现负值冲击系数,不符合实际情况。
除去上述两种规范推荐方法,研究者们也从不同角度对实测冲击系数的计算方法进行研究。王永平等人利用Taylor级数对动位移曲线中的静态趋势项进行逼近。李伟钊等人采用低通滤波法进行动静分离,根据桥梁的动响应数据确定滤波器截止频率。于明策同样采用低通滤波法进行动静分离,但采用了逐步试算法确定截止频率。
综合目前的研究成果,低通滤波法是改进实测冲击系数计算方法的一个可行途径,在公路桥梁荷载试验规程中也建议,在条件允许的情况下可考虑采用低通滤波的方法进行时程曲线中动静位移分离。但公路桥梁荷载试验规程中并未给出具体的滤波方法设计指标,而低通滤波法的滤波结果会由于滤波器类型和截止频率选择的影响,导致滤波结果的不确定性,目前的研究多集中于截止频率的确定方法,但尚未有被广泛接受的滤波器设计方法。
发明内容
本发明的目的是提供一种基于频域幅值谱相似性滤波技术的桥梁冲击系数提取方法,以实现实测动位移信号的动、静分量分离,用于桥梁实测冲击系数的计算。
本发明的目的是通过以下技术方案实现的:
一种基于频域幅值谱相似性滤波技术的桥梁冲击系数提取方法,包括如下步骤:
步骤一、使用同一加载车辆进行两次跑车试验,获得两条数据长度均为N的桥梁动位移曲线信号S1n和S2n,n=1,2,…,N;
步骤二、对步骤一获得的两条桥梁动位移曲线信号做FFT变换,得到两组动位移信号的频域幅值谱A1k和A2k,k=1,2,…,N,将两条跑车曲线的频域幅值谱进行对位比较,计算二者相同频率处的幅值相对偏差倍率Rk:
步骤三、设定一幅值相对偏差倍率阈值,在低频段由低向高对各频点的幅值相对偏差倍率进行判定,幅值相对偏差倍率小于幅值相对偏差倍率阈值的频点即认为由准静态位移控制,予以保留,偏差倍率大于等于幅值相对偏差倍率的频点全部置零,相当于进行理想矩形窗函数低通滤波;
步骤四、对经过低通置零处理的频谱进行FFT逆变换,得到准静态位移曲线S';
步骤五、在原始动位移曲线上量测得到最大动位移Ydmax,在动静分离后得到的准静态位移曲线上量测得到最大静位移Yjmax,则计算冲击系数实测值μ为:
本发明具有如下优点:
1、本发明从一个全新的角度,给出了滤波器截止频率的确定方法。
2、本发明针对准静态分量受测试条件干扰程度低、冲击系数检测不需要信号实时处理等工程特点,通过综合比对多条动位移时程曲线的频域幅值谱结果,可根据多条曲线幅值谱之间的相似性特点,自适应的确定滤波截止频率,有效地解决了常用于动静分离的低通滤波法截止频率估计困难的问题,提高了实测冲击系数的计算准确性。
3、相对于目前普遍采用的实测冲击系数计算方法,本发明不增加检测试验的工作量,数据处理不依赖于人工判读,降低了试验人员主观因素对检测结果的影响,同时低通滤波截止频率的确定考虑到实测信号自身的特点,相较于目前较多采用的低通滤波器设计方法,降低了由于滤波器设计偏差对试验结果的干扰,改善了目前试验因素影响冲击系数检测评定结果的技术现状。
附图说明
图1为计算模型;
图2为检测试验获得的简支梁跨中动位移信号,(a)样本1,(b)样本2;
图3为动位移频域幅值谱;
图4为动位移频域幅值偏差倍率;
图5为准静态位移曲线提取效果。
具体实施方式
下面对本发明的技术方案作进一步的说明,但并不局限于此,凡是对本发明技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围,均应涵盖在本发明的保护范围中。
在桥梁冲击系数的现场检测试验中,直接试验结果数据为桥梁的动位移曲线信号。其中,准静态位移趋势项包含了动位移信号的大部分低频成分,而且由于频率较低,一般的车辆振动、路面颠簸和仪器噪声都不会明显的干扰到此频段,测试结果的稳定性较高,所以若能从动位移曲线测试信号中分离出此低频成分,即可实现动位移曲线动、静位移分量的分离,并进一步完成冲击系数的实测计算。
动位移曲线中的准静态分量曲线形状只与冲击系数检测试验的试验车重量有关,其他的检测试验影响因素的改变不会对准静态分量产生实质性的影响。对于一座检测桥梁,只要加载车辆重量保持不变,经过多次跑车激振获得的各条动位移曲线,虽然可能由于随机因素的干扰表现为不同的波动形态,但其中所包含的准静态位移分量是不变的。另一方面,冲击系数检测试验对检测信号的处理不要求实时性,不要求在信号记录的同时完成滤波操作,所以可允许对检测中获得的多条动位移曲线进行比较分析,利用多条曲线之间的相似性关系来完成准静态位移趋势项的分离。
在此基础上,本发明提出了一种改进的低通滤波法用于准静态位移的提取,此方法通过信号频域幅值谱相似性特点确定低通滤波器的截止频率,具体方法如下:
①使用同一加载车辆进行两次跑车试验,获得两条数据长度均为N的桥梁动位移曲线信号S1n和S2n,(n=1,2,…,N)。
②对记录的动位移信号做FFT变换,得到两组动位移信号的频域幅值谱A1k和A2k,(k=1,2,…,N)。将两条跑车曲线的频域幅值谱进行对位比较,计算二者相同频率处的幅值相对偏差倍率:
③设定一幅值相对偏差倍率(阈值大小根据实际信号的幅值相对偏差倍率Rk确定,根据发明人目前的研究与试验经验,取值大小不应超过0.5)阈值,在低频段由低向高对各频点的幅值相对偏差倍率进行判定,幅值相对偏差倍率小于阈值的频点即认为由准静态位移控制,予以保留。由于准静态位移主要由低频信息控制,所以当偏差开始超过阈值时,此频率即作为低通滤波的截止频率,后续高频频点不再进行判定,全部置零,相当于进行理想矩形窗函数低通滤波。
④对经过低通置零处理的频谱进行FFT逆变换,得到准静态位移曲线S'。
⑤在原始动位移曲线上量测得到最大动位移Ydmax,在动静分离后得到的准静态位移曲线上量测得到最大静位移Yjmax,则可计算冲击系数实测值为:
以跨径L=30m简支梁桥为例,模拟计算30t重两轴卡车驶过桥梁的工况,计算模型如图1所示。对桥梁模型,截面抗弯刚度EI=1.9×1011N·m2;单位长度质量m=31000kg。对车辆模型,簧下质量m1f=1000kg,m1r=2000kg;轮胎弹性模量kft=4800kN·m,krt=9600kN·m;轮胎阻尼系数cft=6kN·s·m,crt=12kN·s·m;悬挂系统弹性模量kfs=1200kN·m,krs=2400kN·m;悬挂系统阻尼系数cfs=5kN·s·m,crs=10kN·s·m;车身质量m2=25000kg;车身转动惯量J2=1×105kg·m2;车轴至重心距离lf=3.7m,lr=1m。
模拟计算两次跑车试验测得的动位移曲线如图2所示,假定两次跑车试验车速相同,均为30km/h,但两次跑车的行车轨迹并不完全相同,计算采用不同的桥面不平整曲线,并计入随机白噪声干扰。
对此两组动位移信号进行FFT变换,其频域幅值谱如图3所示。为保证两组数据的频谱能够准确对位,在进行FFT变换前,应注意将两组信号截取为相同的数据长度。由于信号中低频与高频部分的幅值相差过大,为方便显示,图3中的坐标y轴采用了对数坐标。
计算两条动位移曲线幅值谱之间的幅值相对偏差倍率Rk,结果如图4所示,由于只有低频段对准静态位移的提取有效,为使图形更加清晰,图4中只给出低频部分的一段。
取幅值相对偏差倍率阈值为0.4,根据图4中结果,动位移频谱中频率0.6Hz以下频点的频率偏差小于阈值,可认为是准静态趋势项所控制的频点。将其余高频点全部置零,而后针对新得到的频谱进行FFT逆变换,得到低频滤波后的准静态位移逼近曲线,结果如图5所示。
由于曲线局部位置的高频信息未能保留,所以滤波后的逼近曲线相对于理论真实曲线在转折点和端点的细节上存在细微偏差,但总体上拟合情况良好,峰值位置处偏差为0.64%,可满足工程精度要求。一般在桥梁检测的跑车试验中,同一工况均需进行多次跑车以保证测得有效的数据,所以本发明的方法不增加实际试验的工作量,而且不需要预先对截止频率做出估计,避免了人为滤波器设计偏差对试验结果造成的影响。
由图5中分离得到的的准静态位移曲线可得到冲击系数试验的最大静位移Yjmax=8.67×10-4m。同样由图2中的两条动位移实测曲线可得样本1的最大动位移Ydmax1=10.23×10-4m,样本2的最大动位移Ydmax2=10.69×10-4m。于是样本1对应的冲击系数实测值μ1=0.18,样本2对应的冲击系数实测值μ2=0.23,由于二者所对应的桥面平整情况不相同,所以测得的冲击系数存在差异。
Claims (2)
1.一种基于频域幅值谱相似性滤波技术的桥梁冲击系数提取方法,其特征在于所述方法包括如下步骤:
步骤一、使用同一加载车辆进行两次跑车试验,获得两条数据长度均为N的桥梁动位移曲线信号S1n和S2n,n=1,2,…,N;
步骤二、对步骤一获得的两条桥梁动位移曲线信号做FFT变换,得到两组动位移信号的频域幅值谱A1k和A2k,k=1,2,…,N,将两条跑车曲线的频域幅值谱进行对位比较,计算二者相同频率处的幅值相对偏差倍率Rk,所述幅值相对偏差倍率Rk的计算公式如下:
步骤三、设定一幅值相对偏差倍率阈值,在低频段由低向高对各频点的幅值相对偏差倍率进行判定,幅值相对偏差倍率小于幅值相对偏差倍率阈值的频点即认为由准静态位移控制,予以保留,幅值相对偏差倍率大于等于幅值相对偏差倍率阈值的频点全部置零,相当于进行理想矩形窗函数低通滤波;
步骤四、对经过低通置零处理的频谱进行FFT逆变换,得到准静态位移曲线S';
步骤五、在步骤一获得的桥梁动位移曲线上量测得到最大动位移Ydmax,在步骤四动静分离后得到的准静态位移曲线S'上量测得到最大静位移Yjmax,则计算冲击系数实测值μ为:
2.根据权利要求1所述的基于频域幅值谱相似性滤波技术的桥梁冲击系数提取方法,其特征在于所述幅值相对偏差倍率阈值不超过0.5。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910685464.3A CN110399683B (zh) | 2019-07-27 | 2019-07-27 | 基于频域幅值谱相似性滤波技术的桥梁冲击系数提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910685464.3A CN110399683B (zh) | 2019-07-27 | 2019-07-27 | 基于频域幅值谱相似性滤波技术的桥梁冲击系数提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110399683A CN110399683A (zh) | 2019-11-01 |
CN110399683B true CN110399683B (zh) | 2020-05-05 |
Family
ID=68326269
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910685464.3A Active CN110399683B (zh) | 2019-07-27 | 2019-07-27 | 基于频域幅值谱相似性滤波技术的桥梁冲击系数提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110399683B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114459657B (zh) * | 2022-04-14 | 2022-07-01 | 西南交通大学 | 冲击荷载自动化识别方法、电子设备和存储介质 |
CN115344913B (zh) * | 2022-06-28 | 2024-04-19 | 中铁大桥勘测设计院集团有限公司 | 一种用于桥梁变形评价的弦测值计算方法 |
CN117392515B (zh) * | 2023-09-15 | 2024-04-30 | 中亿丰(苏州)城市建设发展股份有限公司 | 一种基于视觉传感的桥梁结构计量检测系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107228724A (zh) * | 2017-05-24 | 2017-10-03 | 中南大学 | 桥梁动力冲击系数提取方法 |
CN109388871A (zh) * | 2018-09-27 | 2019-02-26 | 成都理工大学 | 一种考虑结构参数不确定性的结构系统随机振动分析方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7668692B2 (en) * | 2005-10-11 | 2010-02-23 | Tatom Frank B | Method for weighing vehicles crossing a bridge |
CN202305156U (zh) * | 2011-11-03 | 2012-07-04 | 中铁西南科学研究院有限公司 | 桥梁动静位移测试装置 |
CN102620900B (zh) * | 2012-04-01 | 2014-06-04 | 长安大学 | 基于动载试验检测桥梁冲击系数的方法 |
CN104215421A (zh) * | 2014-09-11 | 2014-12-17 | 华北水利水电大学 | 一种快速测定桥梁冲击系数的方法 |
EP3187838B1 (en) * | 2015-12-29 | 2020-03-25 | Neostrain spolka z ograniczona odpowiedzialnoscia | System for vehicles weight preselection and evaluation of the technical state of road infrastructure |
CN109002673B (zh) * | 2018-10-08 | 2022-07-05 | 哈尔滨工业大学 | 一种基于车辆制动冲击作用的桥梁基础冲刷识别方法 |
-
2019
- 2019-07-27 CN CN201910685464.3A patent/CN110399683B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107228724A (zh) * | 2017-05-24 | 2017-10-03 | 中南大学 | 桥梁动力冲击系数提取方法 |
CN109388871A (zh) * | 2018-09-27 | 2019-02-26 | 成都理工大学 | 一种考虑结构参数不确定性的结构系统随机振动分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110399683A (zh) | 2019-11-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110399683B (zh) | 基于频域幅值谱相似性滤波技术的桥梁冲击系数提取方法 | |
CN105092467B (zh) | 一种高速铁路车轮踏面擦伤快速监测装置及方法 | |
CN108515984B (zh) | 一种车轮伤损检测方法及装置 | |
CN108804814B (zh) | 一种提取路面谱参数的方法 | |
Sun et al. | Modeling indirect statistics of surface roughness | |
CN104120644A (zh) | 一种基于重力加速度传感器的路面平整度检测方法 | |
CN100520317C (zh) | 固定式车辆动态称重系统 | |
CN104792937A (zh) | 一种基于车载重力加速度传感器的桥头跳车检测评价方法 | |
CN109829252B (zh) | 一种基于影响线识别的桥梁状况快速评级方法 | |
CN110619106A (zh) | 一种桥梁损伤定位方法及其定量方法 | |
CN108180983A (zh) | 自适应时频域混合滤波的应急救援车辆振动位移重构方法 | |
CN107655700A (zh) | 一种道路滑行阻力的确定方法及装置 | |
CN115758289B (zh) | 一种基于多任务学习神经网络的钢轨波磨识别方法 | |
CN113343928A (zh) | 变速路段高速铁路钢轨波磨检测方法及装置、计算机设备 | |
Velinsky et al. | Vehicle energy dissipation due to road roughness | |
CN109919511B (zh) | 既有铁路线形评估方法及系统 | |
CN113415308B (zh) | 钢轨波磨检测方法及装置 | |
Rasmussen et al. | Dynamic spring rate performance of rolling tires | |
CN116296180A (zh) | 基于双轴车接触响应空间位置关系的桥梁阻尼比识别方法 | |
CN116678576A (zh) | 基于车辆振动信号的桥梁振型修正系统 | |
CN113641951B (zh) | 一种基于车辆传感技术的桥梁振型识别方法 | |
CN108827509A (zh) | 一种电动汽车的道路滑行阻力的计算方法和装置 | |
CN105066959A (zh) | 路面纵断面的高程信息获取方法 | |
Klockiewicz et al. | Simulation study of the method of random kinematic road excitation’s reconstruction based on suspension dynamic responses with signal disruptions | |
CN107655708A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |