CN112833837B - 一种基于加速度传感器估计桥梁位移的方法 - Google Patents
一种基于加速度传感器估计桥梁位移的方法 Download PDFInfo
- Publication number
- CN112833837B CN112833837B CN202110083874.8A CN202110083874A CN112833837B CN 112833837 B CN112833837 B CN 112833837B CN 202110083874 A CN202110083874 A CN 202110083874A CN 112833837 B CN112833837 B CN 112833837B
- Authority
- CN
- China
- Prior art keywords
- displacement
- acceleration
- pseudo
- static
- gravity
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
- G01B21/02—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring length, width, or thickness
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- General Engineering & Computer Science (AREA)
- Operations Research (AREA)
- Computing Systems (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
本发明属于铁路桥梁检测领域,具体涉及一种基于加速度传感器估计桥梁位移的方法。考虑桥梁裂缝影响,以多自由度模型考虑桥墩的损伤,采用无参考位移估计方法,具体包括以下步骤:(1):应用有限脉冲响应滤波器来避免不确定边界条件下加速度二重积分增加的位移漂移,由滤波后的加速度近似得到高频零均值动态位移;(2):利用桥墩因重力引起的弯曲角度和重力在加速度轴上的投影来估算拟静力位移;(3):将步骤(1)得到的高频零均值动态位移和步骤(2)得到的拟静力位移叠加得到总估计位移。通过本发明,在信息不受限制的情况下,可估计损伤桩的无参考总位移,且具有较高的精度,实现铁路桥梁损伤构件横向位移的测量。
Description
技术领域
本发明属于铁路桥梁检测领域,具体涉及一种基于加速度传感器估计桥梁位移的方法。
背景技术
近年来,随着铁路运输量不断增加,因为桥梁横向振幅偏大,构件产生横向位移导致运输事故的情况越俩越多。
铁路桥梁在列车穿越事件下最大位移的变化是作为铁路业主和管理人员桥梁健康状况的重要指标。铁路工程师需要观察和监控桥梁的位移情况来决定优先修复位移较大的桥梁。然而,用传感器现场测量桥梁位移是很复杂的,其工作量大准确度不够高。
目前已有利用无参考位移估计法进行横向位移测量的研究,如论文《基于悬臂式结构的无参考位移估计方法及验证》,但是其是基于二自由度模型的研究计算,未考虑桥梁裂缝的影响,准确度不够高。
发明内容
本发明的目的在于提供一种基于加速度传感器估计桥梁位移的方法。
实现本发明目的的技术解决方案为:一种基于加速度传感器估计桥梁位移的方法,考虑桥梁裂缝影响,以多自由度模型考虑桥墩的损伤,采用无参考位移估计方法,具体包括以下步骤:
步骤(1):应用有限脉冲响应滤波器来避免不确定边界条件下加速度二重积分增加的位移漂移,由滤波后的加速度近似得到高频零均值动态位移;
步骤(2):利用桥墩因重力引起的弯曲角度和重力在加速度轴上的投影来估算拟静力位移;
步骤(3):将步骤(1)得到的高频零均值动态位移和步骤(2)得到的拟静力位移叠加得到总估计位移。
进一步的,分别在桥墩顶部,损伤点裂缝顶部,损伤点裂缝下部和桥墩底部处设置一个加速度计,测量得到的加速度用于步骤(1)高频零均值动态位移和步骤(2)拟静力位移的求取。
进一步的,所述步骤(1)中高频零均值动态位移的计算方法如下:
La为(2k+1)阶对角加权矩阵,
Lc是加权矩阵,维数为(2k+1)×(2k+3),
L代表对角线加权矩阵,定义为L=La×Lc,因此:
λ代表最佳的正则化因子,
λ=46.81N-1.95 (5)
其中N为周期窗口对应的采样点数。
进一步的,步骤(2)中,静力位移估算方法如下:
步骤(21):利用桥墩因重力引起的弯曲角度和重力在加速度轴上的投影来估算拟静力位移,用一个加速度计来测量物体的横向加速度并计算弯曲角度,传感器旋转角为θ,重力加速度g投影到x'轴得到的值为Ax,Ax等于加速度计测量的加速度,如式(6)所示:
Ax是重力加速度在x'轴上投影得到的值;R为所需最小分辨率;Q为倾角的测量范围;M为测量角度的最小极限;
步骤(22):集成一个额外的传感器来测量重力加速度g在y'轴上的投影,同样,桥墩因重力引起的弯曲角度,即传感器旋转角为θ,将测量到的y'轴加速度转换为倾角,如式(7)所示;通过融合x'和y'轴上的加速度数据,可以计算出更准确的倾斜角;
Ay是重力加速度g在y'轴上的投影;
采用简单移动平均(SMA)滤波器来消除高频零均值动态响应:
其中θp是伪静态的角度;θi代表了估计的地球旋转加速度;i表示第i次步长;α代表平均点,j表示从0开始取至α-1的第j个点数;
步骤(23):通过步骤(21)和(22)所得伪静态倾斜角推导出伪静态位移,推导公式如下:
其中Δp代表伪静态位移;θ1和θ4分别代表桩顶部和底部的伪静态旋转,Ln代表损伤点以下的桩的长度,Lm代表损伤点以上的桩的长度,θ3表示在多自由度模型中损伤点底部产生的伪静态旋转;θ2表示在多自由度模型中损伤点顶部产生的伪静态旋转。
进一步的,步骤(3)所述总位移估计方法如下:
总估计位移Δt是估计的高频零均值动态位移和伪静态位移的叠加:
Δt=Δd+Δp (10)。
本发明与现有技术相比,其显著优点在于:
(1)本发明提供的铁路桥梁横向位移测量方法基于多自由度模型进行分析计算,考虑了桥梁损伤部位裂缝对高频动态位移和低频伪静态位移估算的影响以及地基条件的影响对,对原有的横向位移测量方法进行了改进,准确度更高。
(2)本发明提供的铁路桥梁横向位移测量方法分析了此类桥梁在列车作用下的横向位移机理和损伤部位裂缝对桥梁整体损伤的影响,可进一步对列车运行舒适度,安全性等进行分析。
附图说明
图1为本发明铁路桥梁伪静态倾斜角计算相关参数参考图;
图2为本发明铁路桥梁多自由度横向位移测量实际应用示意图;
附图标记说明:
1-加速度计Ⅰ,2-加速度计II,3-加速度计III,4-加速度计Ⅳ,5-损伤点裂缝。
具体实施方式
下面结合附图对本发明作进一步详细描述。
一种基于加速度传感器估计桥梁位移的方法,该方法具体为,考虑到桥梁损伤部位裂缝影响,以一种多自由度模型为参考考虑桥墩的横向位移,采用无参考位移估计方法,通过有限脉冲响应滤波器计算铁路桥梁的高频零均值动态位移,在多自由度模型下利用桥墩因重力引起的弯曲角度和重力在加速度轴上的投影来计算低频伪静态位移,将两者叠加得出铁路桥梁损伤构件总横向位移。具体实施时,如图二所示,将加速度计Ⅰ1、加速度计II2、加速度计III3、加速度计Ⅳ4分别由上至下设置于桥墩顶部,损伤点裂缝5顶部,损伤点裂缝5下部,桥墩底部。通过下述方法,求出高频零均值动态位移和低频伪静态位移;
应用有限脉冲响应滤波器来避免不确定边界条件下加速度二重积分增加的位移漂移,由滤波后的加速度近似得到高频零均值动态位移。公式如下所示:
La为(2k+1)阶对角加权矩阵,
Lc是加权矩阵,维数为(2k+1)×(2k+3),
L代表对角线加权矩阵,定义为L=La×Lc,因此:
λ代表最佳的正则化因子,
λ=46.81N-1.95 (5)
其中N为周期窗口对应的采样点数。
静力位移估算方法如下:
①利用桥墩因重力引起的弯曲角度和重力在加速度轴上的投影来估算拟静力位移。用一个加速度计来测量物体的横向加速度并计算弯曲角度。如图一所示,传感器旋转角为θ,重力加速度g投影到x'轴得到的值为Ax,Ax等于加速度计测量的加速度,如式(6)所示:
Ax是重力加速度在x'轴上投影得到的值;R为所需最小分辨率;Q为倾角的测量范围;M为测量角度的最小极限。
②为了进一步提高测量精度,集成一个额外的传感器来测量重力加速度g在y'轴上的投影。同样,如图1所示,桥墩因重力引起的弯曲角度,即传感器旋转角为θ,将测量到的y'轴加速度转换为倾角,如式(7)所示。通过融合x'和y'轴上的加速度数据,可以估计出更准确的倾斜角。
Ay是重力加速度g在y'轴上的投影。
采用简单移动平均(SMA)滤波器来消除高频零均值动态响应:
其中θp是伪静态的角度;θi代表了估计的地球旋转加速度;i表示第i次步长;α代表平均点,j表示从0开始取至α-1的第j个点数。
③如图2所示,通过步骤①②所得伪静态倾斜角可推导出桥墩顶部,损伤点裂缝5顶部,损伤点裂缝5下部,桥墩底部的伪静态位移。推导公式如下:
其中Δp代表伪静态位移;θ1和θ4分别代表桩顶部和底部的伪静态旋转。Ln代表损伤点以下的桩的长度,Lm代表损伤点以上的桩的长度。θ3表示在多自由度模型中损伤点底部产生的伪静态旋转;θ2表示在多自由度模型中损伤点顶部产生的伪静态旋转。
将动态和伪静态估计位移叠加得到总位移估计。
进一步的,所述步骤(3)中所述总位移估计方法如下:
总位移估计总估计位移Δt是动态和伪静态估计位移的叠加:
Δt=Δd+Δp (10)
式(1)为高频动态位移分量的公式;
式(9)为低频伪静态位移分量的公式;
式(10)为叠加高频动态位移和低频伪静态位移的公式。
本发明基于铁路桥梁在列车荷载作用下产生的损伤和横向位移为工程背景,提供一种基于加速度传感器估计桥梁位移的方法,提出一种考虑桥梁裂缝损伤影响的多自由度模型,分析此类桥梁在列车荷载作用下和裂缝损伤影响下的横向振动机理,进一步铁路列车运行的舒适度、安全性进行综合分析,为我国既有铁路桥梁改造工程提供技术支持。
Claims (4)
1.一种基于加速度传感器估计桥梁位移的方法,其特征在于,考虑桥梁裂缝影响,以多自由度模型考虑桥墩的损伤,采用无参考位移估计方法,具体包括以下步骤:
步骤(1):应用有限脉冲响应滤波器来避免不确定边界条件下加速度二重积分增加的位移漂移,由滤波后的加速度近似得到高频零均值动态位移;
步骤(2):利用桥墩因重力引起的弯曲角度和重力在加速度轴上的投影来估算拟静力位移;
步骤(2)中,静力位移估算方法如下:
步骤(21):利用桥墩因重力引起的弯曲角度和重力在加速度轴上的投影来估算拟静力位移,用一个加速度计来测量物体的横向加速度并计算弯曲角度,传感器旋转角为θ,重力加速度g投影到x'轴得到的值为Ax,Ax等于加速度计测量的加速度,如式(6)所示:
Ax是重力加速度在x'轴上投影得到的值;R为所需最小分辨率;Q为倾角的测量范围;M为测量角度的最小极限;
步骤(22):集成一个额外的传感器来测量重力加速度g在y'轴上的投影,同样,桥墩因重力引起的弯曲角度,即传感器旋转角为θ,将测量到的y'轴加速度转换为倾角,如式(7)所示;通过融合x'和y'轴上的加速度数据,可以计算出更准确的倾斜角;
Ay是重力加速度g在y'轴上的投影;
采用简单移动平均(SMA)滤波器来消除高频零均值动态响应:
其中θp是伪静态的角度;θi代表了估计的地球旋转加速度;i表示第i次步长;α代表平均点,j表示从0开始取至α-1的第j个点数;
步骤(23):通过步骤(21)和(22)所得伪静态倾斜角推导出伪静态位移,推导公式如下:
其中Δp代表伪静态位移;θ1和θ4分别代表桩顶部和底部的伪静态旋转,Ln代表损伤点以下的桩的长度,Lm代表损伤点以上的桩的长度,θ3表示在多自由度模型中损伤点底部产生的伪静态旋转;θ2表示在多自由度模型中损伤点顶部产生的伪静态旋转;
步骤(3):将步骤(1)得到的高频零均值动态位移和步骤(2)得到的拟静力位移叠加得到总估计位移。
2.根据权利要求1所述的方法,其特征在于,分别在桥墩顶部,损伤点裂缝顶部,损伤点裂缝下部和桥墩底部处设置一个加速度计,测量得到的加速度用于步骤(1)高频零均值动态位移和步骤(2)拟静力位移的求取。
4.根据权利要求3所述的方法,其特征在于,步骤(3)总估计位移方法如下:
总估计位移Δt是估计的高频零均值动态位移和伪静态位移的叠加:
Δt=Δd+Δp (10)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110083874.8A CN112833837B (zh) | 2021-01-21 | 2021-01-21 | 一种基于加速度传感器估计桥梁位移的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110083874.8A CN112833837B (zh) | 2021-01-21 | 2021-01-21 | 一种基于加速度传感器估计桥梁位移的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112833837A CN112833837A (zh) | 2021-05-25 |
CN112833837B true CN112833837B (zh) | 2022-08-05 |
Family
ID=75929351
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110083874.8A Active CN112833837B (zh) | 2021-01-21 | 2021-01-21 | 一种基于加速度传感器估计桥梁位移的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112833837B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114322911B (zh) * | 2021-12-31 | 2023-08-04 | 重庆大学 | 一种联合卡尔曼滤波的桥梁路面平整度间接精准识别方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001033477A (ja) * | 1999-07-26 | 2001-02-09 | Life Tec Kenkyusho:Kk | 力学量センサ及びこの力学量センサを用いた変位測定装置、変位測定方法、並びに防振装置、防振方法 |
WO2002016926A1 (en) * | 2000-08-23 | 2002-02-28 | Mecon Limited | Monitoring structures |
CN106248335A (zh) * | 2016-07-28 | 2016-12-21 | 浙江工业大学 | 非通航孔桥墩防撞的监测***及实时诊断方法 |
CN110243329A (zh) * | 2019-06-25 | 2019-09-17 | 扬州大学 | 一种轻型墩铁路桥梁横向位移测量方法 |
CN110261132A (zh) * | 2019-06-21 | 2019-09-20 | 扬州大学 | 一种测量车桥***动力响应的模拟***及方法 |
JP6801216B2 (ja) * | 2016-04-20 | 2020-12-16 | セイコーエプソン株式会社 | 損傷検出システム、損傷検出装置、損傷検出方法及びプログラム |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102006005258B4 (de) * | 2006-02-02 | 2011-01-20 | Litef Gmbh | Verfahren zur Ermittlung von Belastungen/Schäden einer mechanischen Struktur |
-
2021
- 2021-01-21 CN CN202110083874.8A patent/CN112833837B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001033477A (ja) * | 1999-07-26 | 2001-02-09 | Life Tec Kenkyusho:Kk | 力学量センサ及びこの力学量センサを用いた変位測定装置、変位測定方法、並びに防振装置、防振方法 |
WO2002016926A1 (en) * | 2000-08-23 | 2002-02-28 | Mecon Limited | Monitoring structures |
JP6801216B2 (ja) * | 2016-04-20 | 2020-12-16 | セイコーエプソン株式会社 | 損傷検出システム、損傷検出装置、損傷検出方法及びプログラム |
CN106248335A (zh) * | 2016-07-28 | 2016-12-21 | 浙江工业大学 | 非通航孔桥墩防撞的监测***及实时诊断方法 |
CN110261132A (zh) * | 2019-06-21 | 2019-09-20 | 扬州大学 | 一种测量车桥***动力响应的模拟***及方法 |
CN110243329A (zh) * | 2019-06-25 | 2019-09-17 | 扬州大学 | 一种轻型墩铁路桥梁横向位移测量方法 |
Non-Patent Citations (1)
Title |
---|
桥梁智能检测技术研究与应用;钟继卫等;《桥梁建设》;20191022;第49卷(第51期);第1-6页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112833837A (zh) | 2021-05-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Chen et al. | Camera-based vibration measurement of the world war I memorial bridge in Portsmouth, New Hampshire | |
Moreu et al. | Reference-free displacements for condition assessment of timber railroad bridges | |
CN112378507B (zh) | 一种基于运动补偿的计算机视觉结构振动监测方法 | |
EP2510154A1 (en) | Method for determining the stress free temperature of the rail and/or the track resistance | |
CN109614674B (zh) | 一种高速铁路无砟轨道层间离缝脱空动态检测方法 | |
CN114964447B (zh) | 测量方法、测量装置、测量***以及存储介质 | |
CN111024421A (zh) | 基于弹性约束支承梁转角影响线的梁桥损伤检测方法 | |
Sekiya et al. | Visualization system for bridge deformations under live load based on multipoint simultaneous measurements of displacement and rotational response using MEMS sensors | |
Wu et al. | Accurate structural displacement monitoring by data fusion of a consumer-grade camera and accelerometers | |
Chen et al. | Stress influence line identification of long suspension bridges installed with structural health monitoring systems | |
CN112833837B (zh) | 一种基于加速度传感器估计桥梁位移的方法 | |
CN111337212A (zh) | 未知状态下基于转角的简支梁最大挠度测定方法 | |
Nagayama et al. | A numerical study on bridge deflection estimation using multi-channel acceleration measurement | |
CN116698316A (zh) | 基于车辆振动信号的桥梁振型识别方法 | |
CN111353238A (zh) | 一种基于车辆传感的桥墩冲刷深度识别方法 | |
Oats et al. | Digital image correlation advances in structural evaluation applications: a review | |
CN110658543B (zh) | 一种基于非接触式测量的高速铁路轨道几何参数检测方法 | |
Ozdagli et al. | Measuring total transverse reference-free displacements for condition assessment of timber railroad bridges: Experimental validation | |
Murray | Dynamic monitoring of rail and bridge displacements using digital image correlation | |
CN113415308A (zh) | 钢轨波磨检测方法及装置 | |
CN116678576A (zh) | 基于车辆振动信号的桥梁振型修正*** | |
CN116296180A (zh) | 基于双轴车接触响应空间位置关系的桥梁阻尼比识别方法 | |
Eick | Structural health monitoring of Inland navigation infrastructure | |
Choi et al. | 3D scanning technique for obtaining road surface and its applications | |
Greisen | Measurement, simulation, and analysis of the mechanical response of railroad track |
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 |