CN116972798A - 基于有限车辆响应的桥梁路面不平整度识别方法 - Google Patents
基于有限车辆响应的桥梁路面不平整度识别方法 Download PDFInfo
- Publication number
- CN116972798A CN116972798A CN202310949230.1A CN202310949230A CN116972798A CN 116972798 A CN116972798 A CN 116972798A CN 202310949230 A CN202310949230 A CN 202310949230A CN 116972798 A CN116972798 A CN 116972798A
- Authority
- CN
- China
- Prior art keywords
- vehicle
- matrix
- bridge
- representing
- response
- 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
- 230000004044 response Effects 0.000 title claims abstract description 57
- 238000000034 method Methods 0.000 title claims abstract description 56
- 239000013598 vector Substances 0.000 claims abstract description 56
- 230000001133 acceleration Effects 0.000 claims abstract description 24
- 238000006073 displacement reaction Methods 0.000 claims abstract description 24
- 230000008878 coupling Effects 0.000 claims abstract description 17
- 238000010168 coupling process Methods 0.000 claims abstract description 17
- 238000005859 coupling reaction Methods 0.000 claims abstract description 17
- 230000003746 surface roughness Effects 0.000 claims abstract description 10
- 238000001914 filtration Methods 0.000 claims abstract description 5
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 4
- 239000011159 matrix material Substances 0.000 claims description 95
- 238000005259 measurement Methods 0.000 claims description 11
- 238000012546 transfer Methods 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 230000005284 excitation Effects 0.000 claims description 4
- 238000013016 damping Methods 0.000 claims description 3
- 230000003993 interaction Effects 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 description 9
- 230000036541 health Effects 0.000 description 6
- 238000012544 monitoring process Methods 0.000 description 6
- 238000012545 processing Methods 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000007405 data analysis Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 238000013500 data storage Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 230000008520 organization Effects 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
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/30—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 roughness or irregularity of surfaces
-
- E—FIXED CONSTRUCTIONS
- E01—CONSTRUCTION OF ROADS, RAILWAYS, OR BRIDGES
- E01C—CONSTRUCTION OF, OR SURFACES FOR, ROADS, SPORTS GROUNDS, OR THE LIKE; MACHINES OR AUXILIARY TOOLS FOR CONSTRUCTION OR REPAIR
- E01C23/00—Auxiliary devices or arrangements for constructing, repairing, reconditioning, or taking-up road or like surfaces
- E01C23/01—Devices or auxiliary means for setting-out or checking the configuration of new surfacing, e.g. templates, screed or reference line supports; Applications of apparatus for measuring, indicating, or recording the surface configuration of existing surfacing, e.g. profilographs
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Civil Engineering (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Structural Engineering (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Architecture (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
本发明公开了一种基于有限车辆响应的桥梁路面不平整度识别方法,包括如下步骤:步骤一:将车载传感器布置在测量车车轴的中央位置或正上方的车厢上;步骤二:采用牵引车引导测量车匀速驶过待测桥梁,采集测量车的竖向响应;步骤三:构建车辆和桥梁的平衡方程,根据定义的状态向量得到车桥耦合系统的状态空间方程,并将车桥耦合系统的状态空间方程进行离散化表示;步骤四:以加速度响应作为观测值,将车辆加速度进行离散化表示;步骤五:以用k时刻状态向量表示k+2时刻的状态向量,得到k+2时刻的车辆位移;步骤六:考虑实际中存在的观测噪声,得到改进的系统观测方程;步骤七:根据扩展卡尔曼滤波算法,求解得到路面不平整度。
Description
技术领域
本发明属于桥梁健康监测技术领域,具体的为一种基于有限车辆响应的桥梁路面不平整度识别方法。
背景技术
目前,针对桥梁结构的健康监测技术主要分为直接测量和间接测量技术。直接测量技术通常是在桥梁结构表面或者关键部位安装、布置一系列的传感器来监测、记录桥梁在运营过程中的振动响应。间接测量主要包括图像技术、无人机技术和GPS等等。其中,车辆扫描法作为一种新型的桥梁状态间接评估技术,由于其具有高效、可移动性强以及更好的经济性等优势,近来年在中小型桥梁健康监测领域受到广泛的关注。
车辆扫描法主要的思想是:通过装备传感器的移动车辆来间接感知桥梁的振动,进而评估桥梁的动力特性,为桥梁的健康监测提供基础资料。然而,路面不平整度是车辆扫描法应用环节不可忽略的因素,它会显著削弱车辆响应谱中桥梁频率的可见性。因此,从车辆响应中辨识路面不平整度是车辆扫描法推广应用中的关键一环。近年来,基于车辆响应的桥梁粗糙度辨识技术应运而生,该技术主要涉及较为复杂的车辆模型,在实际应用过程中需要投入一定数量的传感器,即需要在车辆不同位置布置多个传感器,并结合先进的信号处理技术来识别路面不平整度。然而,由于不同传感器记录的桥梁响应受外界干扰,难以保持一致和稳定的状态,导致研究成果难以用于实际。此外,大量的传感器投入会提高桥梁健康监测的成本,极大地削弱了车辆扫描技术在经济性和操作便捷性等方面的优势。
发明内容
有鉴于此,本发明的目的在于提供一种基于有限车辆响应的桥梁路面不平整度识别方法,仅需采用单个车载传感器,即可识别得到桥梁路面不平整度。
为达到上述目的,本发明提供如下技术方案:
一种基于有限车辆响应的桥梁路面不平整度识别方法,包括如下步骤:
步骤一:将车载传感器布置在测量车车轴的中央位置或布置在位于测量车车轴中央正上方的车厢上;
步骤二:采用牵引车引导测量车匀速驶过待测桥梁,在测量车在待测桥梁上行驶过程中,信号采集系统记录测量车的竖向响应;
步骤三:构建车辆和桥梁的平衡方程,根据定义的状态向量Z得到车桥耦合系统的状态空间方程,并将车桥耦合系统的状态空间方程进行离散化表示;
步骤四:以加速度响应作为观测值,将车辆加速度进行离散化表示;
步骤五:以用k时刻状态向量Zk表示k+2时刻的状态向量Zk+2,得到k+2时刻的车辆位移yv,k+2;
步骤六:考虑实际中存在的观测噪声vk+1,得到改进的系统观测方程;
步骤七:根据扩展卡尔曼滤波算法,求解得到路面不平整度
进一步,所述车载传感器采用位移传感器,所述步骤二中,信号采集系统记录测量车的竖向位移响应。
进一步,所述步骤三中,车辆和桥梁的平衡方程为:
其中:Mb,Cb和Kb分别为桥梁的质量矩阵、阻尼矩阵和刚度矩阵;mv为车辆的质量;kv为车辆刚度;yv为车辆的竖向位移;为车辆竖向加速度;qb为桥梁节点的位移向量,维度为1×n,n为桥梁的自由度数目;/>为桥梁节点的速度向量;/>为桥梁节点的加速度向量;rc为待识别的路面不平整度;H为包含Hamilton插值函数的行向量;F为车辆与桥梁之间的相互作用力;且:
H1×n=[0,…0,N,0,…0]
F=HT[mvg+kv(yv-rc)-kvHqb]
其中:N为Hamilton插值函数;g为重力加速度。
进一步,所述步骤三中,状态向量Z定义为:
其中:表示车辆竖向速度;
车桥耦合系统的状态空间方程为:
其中:表示状态向量Z对时间t的一阶导数;t表示离散时间;A和B分别为状态方程的状态矩阵和输入矩阵;F为已知荷载矩阵;且:
车桥耦合系统的状态空间方程离散化表示为:
Zk+1=AkZk+Bkrc,k+Fk+wk
Ak=exp(Adt)≈(I+AΔt)
Bk=[exp(Adt)-I]A-1B≈dtB
Fk=[exp(Adt)-I]A-1F
其中:I为单位矩阵;w表示系统噪音。
进一步,所述步骤四中,车辆加速度为:
将车辆加速度离散化表示为:
其中:L表示输出矩阵;E表示传递矩阵。
进一步,所述步骤五中,k+2时刻的状态向量Zk+2为:
Zk+2=Ak+1AkZk+Ak+1Bkrc,k+Ak+1Fk+Bk+1rc,k+1+Fk+1+wk
Ak+1Ak=[U1,k+1U2,k+1U3,k+1U4,k+1]T
Ak+1Bk=[V1,k+1V2,k+1V3,k+1V4,k+1]T
k+2时刻的车辆位移yv,k+2为:
yv,k+2=U3,k+1Zk+V3,k+1rc,k
其中:U1,k+1表示Ak+1Ak矩阵第一行对应的子矩阵;U2,k+1表示Ak+1Ak矩阵第二至n+1行对应的子矩阵;U3,k+1表示Ak+1Ak矩阵第n+2行对应的子矩阵;U4,k+1表示Ak+1Ak矩阵第n+3至2n+2行对应的子矩阵;V1,k+1表示Ak+1Bk矩阵第一行对应的子矩阵;V2,k+1表示Ak+1Bk矩阵第二至n+1行对应的子矩阵;V3,k+1表示Ak+1Bk矩阵第n+2行对应的子矩阵;V4,k+1表示Ak+1Bk矩阵第n+3至2n+2行对应的子矩阵。
进一步,所述步骤六中,改进的系统观测方程为;
Yk+1=CkZk+1+Dkrc,k+1+vk+1其中:Yk+1表示改进的观测向量;vk+1表示观测噪声;Ck为改进的输出矩阵;Dk为改进的传递矩阵。
进一步,所述步骤七中,求解得到路面不平整度的步骤为:
71)估计状态初始值设置:
其中:Z0和r0分别为状态向量和未知输入的初始值;和/>分别为状态向量和未知输入的初始估计值;Pz,0|0为状态向量协方差初始值;E[·]表示数学期望;
72)时间更新阶段:
其中:和/>分别表示上一迭代步和当前迭代步计算的状态向量;Pz,k+1|k表示状态预测值/>的误差协方差矩阵;Pz,k|k表示状态预测值/>的误差协方差矩阵;Qk+1为k+1时刻噪声向量wk+1的协方差矩阵;
73)计算卡尔曼增益Kz,k+1:
其中:Ck+1|k表示由上一步迭代计算得到的输出矩阵;Rk+1表示观测噪声vk+1的协方差矩阵;
74)估计未知输入:
其中:Sk+1表示未知激励rk+1的误差协方差矩阵;Dk+1|k表示由上一步迭代计算得到的传递矩阵;Yk+1表示改进的观测向量;I表示单位矩阵;
75)测量更新阶段:
其中:和/>分别为状态向量Z和未知输入r的估计值。
本发明的有益效果在于:
本发明基于有限车辆响应的桥梁路面不平整度识别方法,首先,将路面不平整度作为车桥耦合系统的未知激励,并采用位移传感器采集测量车过桥时的竖向响应;然后,定义车桥耦合系统的状态向量和观测向量,分别建立相应的状态方程和离散方程,得到车桥耦合状态空间方程的系统矩阵;最后,定义状态初始值,联合卡尔曼滤波技术以及获得的系统矩阵重构路面不平整度。综上,本发明基于有限车辆响应的桥梁路面不平整度识别方法,仅需采用单个车载传感器,即可识别得到桥梁路面不平整度;由于车载传感器设为一个并布置在测量车车轴的中央位置或布置在位于测量车车轴中央正上方的车厢上,不仅具有经济性和操作便捷性等方面的优点,而且能够避免测量得到的竖向响应受外界干扰,从而能够保持测量的一致性和稳定性。
附图说明
为了使本发明的目的、技术方案和有益效果更加清楚,本发明提供如下附图进行说明:
图1为本发明基于有限车辆响应的桥梁路面不平整度识别方法实施例的流程图;
图2为路面不平整度辨识系统实施例的原理图;
图3为测量车系统的结构示意图;
图4为被测桥梁的数学模型;
图5为不同车辆行驶速度下识别得到的桥面不平整度;
图6为不同环境噪声下识别得到的桥面不平整度;
图7为将试验车置于路面识别得到的路面不平整度。
具体实施方式
下面结合附图和具体实施例对本发明作进一步说明,以使本领域的技术人员可以更好的理解本发明并能予以实施,但所举实施例不作为对本发明的限定。
如图1所示,本实施例基于有限车辆响应的桥梁路面不平整度识别方法,包括如下步骤。
步骤一:将车载传感器15布置在测量车10车轴14的中央位置或布置在位于测量车10车轴14中央正上方的车厢16上。本实施例中,车载传感器15采用位移传感器。
步骤二:采用牵引车11引导测量车10匀速驶过待测桥梁,在测量车10在待测桥梁上行驶过程中,信号采集系统记录测量车的竖向响应。本实施例中,信号采集系统记录测量车的竖向位移响应。
步骤三:构建车辆和桥梁的平衡方程,根据定义的状态向量Z得到车桥耦合系统的状态空间方程,并将车桥耦合系统的状态空间方程进行离散化表示。
车辆和桥梁的平衡方程为:
其中:Mb,Cb和Kb分别为桥梁的质量矩阵、阻尼矩阵和刚度矩阵;mv为车辆的质量;kv为车辆刚度;yv为车辆的竖向位移;为车辆竖向加速度;qb为桥梁节点的位移向量,维度为1×n,n为桥梁的自由度数目;/>为桥梁节点的速度向量;/>为桥梁节点的加速度向量;rc为待识别的路面不平整度;H为包含Hamilton插值函数的行向量,除了与车辆作用相对应的梁单元外,其他项均为零;F为车辆与桥梁之间的相互作用力;且:
H1×n=[0,…0,N,0,…0]
F=HT[mvg+kv(yv-rc)-kvHqb]
其中:N为Hamilton插值函数;g为重力加速度;
状态向量Z定义为:
其中:表示车辆竖向速度;
得到车桥耦合系统的状态空间方程为:
其中:表示状态向量Z对时间Z的一阶导数;t表示离散时间;A和B分别为状态方程的状态矩阵和输入矩阵;F为已知荷载矩阵;且:
考虑实际观测到的序列是离散的,且存在系统噪声w。因此,将状态空间方程进行离散化表示。本实施例中,车桥耦合系统的状态空间方程离散化表示为:
Zk+1=AkZk+Bkrc,k+Fk+wk
Ak=exp(Adt)≈(I+AΔt)
Bk=[exp(Adt)-I]A-1B≈dtB
Fk=[exp(Ad)-I]A-1F
其中:I为单位矩阵;w表示系统噪音。
步骤四:以加速度响应作为观测值,将车辆加速度进行离散化表示。
车辆的加速度可以通过测量的车辆位移信号求导获得,且对于车辆扫描法而言,能够测量的仅为车辆响应。具体的,车辆加速度为:
将车辆加速度离散化表示为:
其中:L表示输出矩阵;E表示传递矩阵。
步骤五:以用k时刻状态向量Zk表示k+2时刻的状态向量Zk+2,得到k+2时刻的车辆位移yv,k+2。
k+2时刻的状态向量Zk+2为:
Zk+2=Ak+1AkZk+Ak+1Bkrc,k+Ak+1Fk+Bk+1rc,k+1+Fk+1+wk
Ak+1Ak=[U1,k+1U2,k+1U3,k+1U4,k+1]T
Ak+1Bk=[V1,k+1V2,k+1V3,k+1V4,k+1]T
k+2时刻的车辆位移yv,k+2为:
yv,k+2=U3,k+1Zk+V3,k+1rc,k
其中:U1,k+1表示Ak+1Ak矩阵第一行对应的子矩阵;U2,k+1表示Ak+1Ak矩阵第二至n+1行对应的子矩阵;U3,k+1表示Ak+1Ak矩阵第n+2行对应的子矩阵;U4,k+1表示Ak+1Ak矩阵第n+3至2n+2行对应的子矩阵;V1,k+1表示Ak+1Bk矩阵第一行对应的子矩阵;V2,k+1表示Ak+1Bk矩阵第二至n+1行对应的子矩阵;V3,k+1表示Ak+1Bk矩阵第n+2行对应的子矩阵;V4,k+1表示Ak+1Bk矩阵第n+3至2n+2行对应的子矩阵。
步骤六:考虑实际中存在的观测噪声vk+1,得到改进的系统观测方程。
令k=k+1,改进的系统观测方程为:
Yk+1=CkZk+1+Dkrc,k+1+vk+1
其中::Yk+1表示改进的观测向量;vk+1表示观测噪声;Ck为改进的输出矩阵;Dk为改进的传递矩阵。
车辆位移响应yv,k+3可以直接由布置在车辆上的位移传感器记录得到,车辆加速度响应可以通过已获得的位移响应对时间dt两次微分得到。同理,若车辆配置的加速度传感器,则通过对车辆加速度响应对时间dt两次积分得到车辆位移响应。
步骤七:根据扩展卡尔曼滤波算法,求解得到路面不平整度
求解得到路面不平整度的步骤为:
71)估计状态初始值设置:
其中:Z0和r0分别为状态向量和未知输入的初始值;和/>分别为状态向量和未知输入的初始估计值;Pz,0|0为状态向量协方差初始值;E[·]表示数学期望;
72)时间更新阶段:
其中:下标k|k表示该数值由同步迭代计算得到,k+1|k表示该数值由上一步迭代得到。具体的,和/>分别表示上一迭代步和当前迭代步计算的状态向量;Pz,k+1|k表示状态预测值/>的误差协方差矩阵;Pz,k|k表示状态预测值/>的误差协方差矩阵;Qk+1表示k+1时刻系统噪声向量wk+1的协方差矩阵;
73)计算卡尔曼增益Kz,k+1:
其中:Ck+1|k表示由上一步迭代计算得到的输出矩阵;Rk+1表示观测噪声vk+1的协方差矩阵;
74)估计未知输入:
其中:Sk+1表示未知激励rk+1的误差协方差矩阵;Dk+1|k表示由上一步迭代计算得到的传递矩阵;Yk+1表示改进的观测向量;I表示单位矩阵;
75)测量更新阶段:
其中:和/>分别为状态向量Z和未知输入r的估计值。
如图2所示,为适用于本发明基于有限车辆响应的桥梁路面不平整度识别方法的路面不平整度辨识系统的原理图。本实施例的路面不平整度辨识系统,包括现场测量系统、数据分析处理平台和数据输出及显示终端。
本实施例中,现场测量系统包括测量车系统、数据采集模块、数据转换模块、数据通信模块和数据存储模块。本实施例的测量车系统用于实时测量竖向响应数据,数据采集模块采集传感器测量得到的竖向响应数据,数据转换模块将采集得到的竖向响应数据进行数据转换,数据通信模块将转换后的竖向响应数据传输至数据存储模块进行存储。如图3所示,本实施例的测量车系统包括测量车10和牵引车11,牵引车11用于引导测量车10移动,且测量车10可绕一转轴12相对于牵引车11转动、并可沿一滑轴13相对于牵引车11移动,转轴12和滑轴13相互垂直。本实施例中,转轴12与牵引车11转动配合,转轴12上设有滑套,滑轴13与滑套滑动配合,滑轴13与测量车10固定连接。测量车10包括车厢16,车厢16上安装有一根车轴14,车轴14分别与转轴12和滑轴13垂直,滑轴13固定安装在车厢16上。本实施例的车轴14中央位置处布置有一颗用于采集测量车10在被测路面上行驶过程中的竖向响应的车载传感器15。当然,在其他实施方式中,也可以将车载传感器15布置在位于车轴14中央位置正上方的车厢16上。
本实施例中,数据分析处理平台用于根据采集得到的竖向响应数据辨识路面不平整度。本实施例中,车载传感器15采用位移传感器,数据采集模块采集位移传感器测量得到的竖向位移响应。
本实施例中,数据输出及显示终端用于实时输出和显示所述数据分析处理平台的计算结果。
下面结合具体实例对本发明基于有限车辆响应的桥梁路面不平整度识别方法进行说明。
在数值验证中,采用图4所示的被测桥梁的数学模型。图4中桥梁的参数设置如下:桥梁长度L=25m,截面尺寸为A=3.2m2,桥梁密度为ρ=4800kg/m3,弹性模量E=2.75×1010N/m2。测量车参数设置如下:车辆刚度kv=200kN/m,车体质量mv=14,000kg,移动速度v=2m/s。采用国际标准组织ISO(International Standard Organization)8608所定义的功率谱密度函数PSD(Power Spectral Density)法进行模拟。ISO8608标准将路面平整度划分为A-E五个不同等级,其中A级代表路面情况最好,E级代表路面最粗糙。本实施例中,路面不平整度等级考虑为C级。
为验证本实施例基于有限车辆响应的桥梁路面不平整度识别方法的鲁棒性,开展数值研究验证了不同车辆行驶速度,不同环境噪声干扰下桥面不平整度的辨识效果。此外,为验证本实施例基于有限车辆响应的桥梁路面不平整度识别方法的通用性,在数值研究中特意模拟试验车置于路面的情形。上述三种工况下的粗糙度案例如下:
工况一:不同车辆行驶车速下利用本实施例基于有限车辆响应的桥梁路面不平整度识别方法识别的桥面不平整度,如图5所示;
工况二:不同环境噪声下利用本实施例基于有限车辆响应的桥梁路面不平整度识别方法识别的桥面不平整度,如图6所示;
工况三:不同应用场景利用本实施例基于有限车辆响应的桥梁路面不平整度识别方法识别的路面不平整度,如图7所示。
从三种工况的数值验证结果可以看出,计算得到的粗糙度与理论值无论在时域还是频域都表现较好的吻合度。此外,本实施例基于有限车辆响应的桥梁路面不平整度识别方法具有广泛的适用性,不仅适用于桥梁路面不平整度辨识,同样适用于其他路面的粗糙度辨识。本实施例路面不平整度辨识系统的粗糙度辨识方法辨识过程高效,辨识流程简单,辨识结果准确,能够为量大面广的桥梁结构桥面不平整度辨识提供新的技术支撑,服务桥梁健康监测,运营管理与维护。
此外,本实施例可以很好地拓展致其他领域的荷载辨识问题,例如车辆行车舒适度评估。
以上所述实施例仅是为充分说明本发明而所举的较佳的实施例,本发明的保护范围不限于此。本技术领域的技术人员在本发明基础上所作的等同替代或变换,均在本发明的保护范围之内。本发明的保护范围以权利要求书为准。
Claims (8)
1.一种基于有限车辆响应的桥梁路面不平整度识别方法,其特征在于:包括如下步骤:
步骤一:将车载传感器布置在测量车车轴的中央位置或布置在位于测量车车轴中央正上方的车厢上;
步骤二:采用牵引车引导测量车匀速驶过待测桥梁,在测量车在待测桥梁上行驶过程中,信号采集系统记录测量车的竖向响应;
步骤三:构建车辆和桥梁的平衡方程,根据定义的状态向量Z得到车桥耦合系统的状态空间方程,并将车桥耦合系统的状态空间方程进行离散化表示;
步骤四:以加速度响应作为观测值,将车辆加速度进行离散化表示;
步骤五:以用k时刻状态向量Zk表示k+2时刻的状态向量Zk+2,得到k+2时刻的车辆位移yv,k+2;
步骤六:考虑实际中存在的观测噪声vk+1,得到改进的系统观测方程;
步骤七:根据扩展卡尔曼滤波算法,求解得到路面不平整度
2.根据权利要求1所述基于有限车辆响应的桥梁路面不平整度识别方法,其特征在于:所述车载传感器采用位移传感器,所述步骤二中,信号采集系统记录测量车的竖向位移响应。
3.根据权利要求2所述基于有限车辆响应的桥梁路面不平整度识别方法,其特征在于:所述步骤三中,车辆和桥梁的平衡方程为:
其中:Mb,Cb和Kb分别为桥梁的质量矩阵、阻尼矩阵和刚度矩阵;mv为车辆的质量;kv为车辆刚度;yv为车辆的竖向位移;为车辆竖向加速度;qb为桥梁节点的位移向量,维度为1×n,n为桥梁的自由度数目;/>为桥梁节点的速度向量;/>为桥梁节点的加速度向量;rc为待识别的路面不平整度;H为包含Hamilton插值函数的行向量;F为车辆与桥梁之间的相互作用力;且:
H1×n=[0,…0,N,0,…0]
F=HT[mvg+kv(yv-rc)-kvHqb]
其中:N为Hamilton插值函数;g为重力加速度。
4.根据权利要求3所述基于有限车辆响应的桥梁路面不平整度识别方法,其特征在于:所述步骤三中,状态向量Z定义为:
其中:表示车辆竖向速度;
车桥耦合系统的状态空间方程为:
其中:表示状态向量Z对时间t的一阶导数;t表示离散时间;A和B分别为状态方程的状态矩阵和输入矩阵;F为已知荷载矩阵;且:
车桥耦合系统的状态空间方程离散化表示为:
Zk+1=AkZk+Bkrc,k+Fk+wk
Ak=exp(Adt)≈(I+AΔt)
Bk=[exp(Adt)-I]A-1B≈dtB
Fk=[exp(Adt)-I]A-1F
其中:I为单位矩阵;w表示系统噪音。
5.根据权利要求4所述基于有限车辆响应的桥梁路面不平整度识别方法,其特征在于:所述步骤四中,车辆加速度为:
将车辆加速度离散化表示为:
其中:L表示输出矩阵;E表示传递矩阵。
6.根据权利要求5所述基于有限车辆响应的桥梁路面不平整度识别方法,其特征在于:所述步骤五中,k+2时刻的状态向量Zk+2为:
Zk+2=Ak+1AkZk+Ak+1Bkrc,k+Ak+1Fk+Bk+1rc,k+1+Fk+1+wk
Ak+1Ak=[U1,k+1U2,k+1U3,k+1U4,k+1]T
Ak+1Bk=[V1,k+1V2,k+1V3,k+1V4,k+1]T
k+2时刻的车辆位移yv,k+2为:
yv,k+2=U3,k+1Zk+V3,k+1rc,k
其中:U1,k+1表示Ak+1Ak矩阵第一行对应的子矩阵;U2,k+1表示Ak+1Ak矩阵第二至n+1行对应的子矩阵;U3,k+1表示Ak+1Ak矩阵第n+2行对应的子矩阵;U4,k+1表示Ak+1Ak矩阵第n+3至2n+2行对应的子矩阵;V1,k+1表示Ak+1Bk矩阵第一行对应的子矩阵;V2,k+1表示Ak+1Bk矩阵第二至n+1行对应的子矩阵;V3,k+1表示Ak+1Bk矩阵第n+2行对应的子矩阵;V4,k+1表示Ak+1Bk矩阵第n+3至2n+2行对应的子矩阵。
7.根据权利要求6所述基于有限车辆响应的桥梁路面不平整度识别方法,其特征在于:所述步骤六中,改进的系统观测方程为;
其中:Yk+1表示改进的观测向量;vk+1表示观测噪声;Ck为改进的输出矩阵;Dk为改进的传递矩阵。
8.根据权利要求7所述基于有限车辆响应的桥梁路面不平整度识别方法,其特征在于:所述步骤七中,求解得到路面不平整度的步骤为:
71)估计状态初始值设置:
其中:Z0和r0分别为状态向量和未知输入的初始值;和/>分别为状态向量和未知输入的初始估计值;Pz,0|0为状态向量协方差初始值;E[·]表示数学期望;
72)时间更新阶段:
其中:和/>分别表示当前迭代步和上一迭代步计算的状态向量;Pz,k+1|k表示状态预测值/>的误差协方差矩阵;Pz,k|k表示状态预测值/>的误差协方差矩阵;Qk+1表示k+1时刻系统噪声向量wk+1的协方差矩阵;
73)计算卡尔曼增益Kz,k+1:
其中:Ck+1|k表示由上一步迭代计算得到的输出矩阵;Rk+1表示观测噪声vk+1的协方差矩阵;
74)估计未知输入:
其中:Sk+1表示未知激励rk+1的误差协方差矩阵;Dk+1|k表示由上一步迭代计算得到的传递矩阵;Yk+1表示改进的观测向量;I表示单位矩阵;
75)测量更新阶段:
其中:和/>分别为状态向量Z和未知输入r的估计值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310949230.1A CN116972798A (zh) | 2023-07-31 | 2023-07-31 | 基于有限车辆响应的桥梁路面不平整度识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310949230.1A CN116972798A (zh) | 2023-07-31 | 2023-07-31 | 基于有限车辆响应的桥梁路面不平整度识别方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116972798A true CN116972798A (zh) | 2023-10-31 |
Family
ID=88470892
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310949230.1A Pending CN116972798A (zh) | 2023-07-31 | 2023-07-31 | 基于有限车辆响应的桥梁路面不平整度识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116972798A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117454718A (zh) * | 2023-12-05 | 2024-01-26 | 重庆大学 | 一种考虑任意边界条件下识别桥梁频率的车辆扫描方法 |
-
2023
- 2023-07-31 CN CN202310949230.1A patent/CN116972798A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117454718A (zh) * | 2023-12-05 | 2024-01-26 | 重庆大学 | 一种考虑任意边界条件下识别桥梁频率的车辆扫描方法 |
CN117454718B (zh) * | 2023-12-05 | 2024-05-31 | 重庆大学 | 一种考虑任意边界条件下识别桥梁频率的车辆扫描方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Imine et al. | Road profile input estimation in vehicle dynamics simulation | |
Keenahan et al. | The use of a dynamic truck–trailer drive-by system to monitor bridge damping | |
CN110704801B (zh) | 桥梁集群结构运营安全智能监测与快速检测成套方法 | |
CN105000033B (zh) | 一种轨道几何不平顺的检测评估系统 | |
CN109855823A (zh) | 一种利用测试车进行桥梁结构损伤识别的方法 | |
Yang et al. | Using a single-DOF test vehicle to simultaneously retrieve the first few frequencies and damping ratios of the bridge | |
Li et al. | A two-step drive-by bridge damage detection using dual Kalman filter | |
CN107588915A (zh) | 一种桥梁影响线识别方法及系统 | |
CN101266190A (zh) | 正常交通流下斜拉桥模态参数测定装置及方法 | |
CN116972798A (zh) | 基于有限车辆响应的桥梁路面不平整度识别方法 | |
CN113353085B (zh) | 一种基于卡尔曼滤波理论的路面不平度识别方法 | |
CN110543706A (zh) | 一种基于车辆刹车作用的在役桥梁支座损伤诊断方法 | |
CN108228994B (zh) | 在越野道路随机激励下车辆及工作装置的应力计算方法 | |
Yang et al. | Closely spaced modes of bridges estimated by a hybrid time–frequency method using a multi-sensor scanning vehicle: Theory and practice | |
CN113432815B (zh) | 一种基于测量车振动响应的桥面响应重构方法 | |
CN116678576A (zh) | 基于车辆振动信号的桥梁振型修正系统 | |
Li et al. | Indirect Frequency Identification of Footbridges with Pedestrians Using the Contact-Point Response of Shared Scooters | |
CN114322911B (zh) | 一种联合卡尔曼滤波的桥梁路面平整度间接精准识别方法 | |
CN116698316A (zh) | 基于车辆振动信号的桥梁振型识别方法 | |
CN116972797A (zh) | 路面不平整度辨识系统 | |
CN113420362B (zh) | 一种适用于时变中小跨径桥梁的损伤识别方法及装置 | |
CN116558745A (zh) | 一种移动与非移动车辆共同作用下的桥梁频率提取方法 | |
CN116296180A (zh) | 基于双轴车接触响应空间位置关系的桥梁阻尼比识别方法 | |
Ngwangwa | Calculation of road profiles by reversing the solution of the vertical ride dynamics forward problem | |
WO2020261959A1 (ja) | 推定装置、推定方法、およびプログラム |
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 |