CN114913688B - 一种交通连续流作用下桥梁耦合振动响应预警方法 - Google Patents

一种交通连续流作用下桥梁耦合振动响应预警方法 Download PDF

Info

Publication number
CN114913688B
CN114913688B CN202210547758.1A CN202210547758A CN114913688B CN 114913688 B CN114913688 B CN 114913688B CN 202210547758 A CN202210547758 A CN 202210547758A CN 114913688 B CN114913688 B CN 114913688B
Authority
CN
China
Prior art keywords
bridge
traffic flow
flow state
matrix
next moment
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
Application number
CN202210547758.1A
Other languages
English (en)
Other versions
CN114913688A (zh
Inventor
韩智强
李华腾
王晓军
赵洲
陈嘉毅
刘世忠
刘秀英
杨春霞
张驿丹
晋民杰
贾志绚
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Taiyuan University of Science and Technology
Original Assignee
Taiyuan University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Taiyuan University of Science and Technology filed Critical Taiyuan University of Science and Technology
Priority to CN202210547758.1A priority Critical patent/CN114913688B/zh
Publication of CN114913688A publication Critical patent/CN114913688A/zh
Application granted granted Critical
Publication of CN114913688B publication Critical patent/CN114913688B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G08—SIGNALLING
    • G08G—TRAFFIC CONTROL SYSTEMS
    • G08G1/00—Traffic control systems for road vehicles
    • G08G1/01—Detecting movement of traffic to be counted or controlled
    • G—PHYSICS
    • G08—SIGNALLING
    • G08B—SIGNALLING OR CALLING SYSTEMS; ORDER TELEGRAPHS; ALARM SYSTEMS
    • G08B31/00—Predictive alarm systems characterised by extrapolation or other computation using updated historic data
    • G—PHYSICS
    • G08—SIGNALLING
    • G08G—TRAFFIC CONTROL SYSTEMS
    • G08G1/00—Traffic control systems for road vehicles
    • G08G1/01—Detecting movement of traffic to be counted or controlled
    • G08G1/0104—Measuring and analyzing of parameters relative to traffic conditions
    • G08G1/0125—Traffic data processing
    • G—PHYSICS
    • G08—SIGNALLING
    • G08G—TRAFFIC CONTROL SYSTEMS
    • G08G1/00—Traffic control systems for road vehicles
    • G08G1/01—Detecting movement of traffic to be counted or controlled
    • G08G1/0104—Measuring and analyzing of parameters relative to traffic conditions
    • G08G1/0137—Measuring and analyzing of parameters relative to traffic conditions for specific applications
    • G—PHYSICS
    • G08—SIGNALLING
    • G08G—TRAFFIC CONTROL SYSTEMS
    • G08G1/00—Traffic control systems for road vehicles
    • G08G1/01—Detecting movement of traffic to be counted or controlled
    • G08G1/052—Detecting movement of traffic to be counted or controlled with provision for determining speed or overspeed
    • G—PHYSICS
    • G08—SIGNALLING
    • G08G—TRAFFIC CONTROL SYSTEMS
    • G08G1/00—Traffic control systems for road vehicles
    • G08G1/065—Traffic control systems for road vehicles by counting the vehicles in a section of the road or in a parking area, i.e. comparing incoming count with outgoing count
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • Emergency Management (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
  • Traffic Control Systems (AREA)

Abstract

本发明公开了一种交通连续流作用下桥梁耦合振动响应预警方法，包括步骤：一、初始化桥梁交通流；二、实时采集桥梁上交通流；三、构建当前交通流状态下的车桥耦合矩阵并求解对应的桥梁动力响应增量；四、获取当前交通流状态下桥梁的动力响应结果；五、采集下一时刻桥梁上交通流、构建下一时刻交通流状态下的车桥耦合矩阵并求解对应的桥梁动力响应增量；六、获取下一时刻交通流状态下桥梁的动力响应结果；七、计算下一时刻交通流状态下桥梁冲击系数；八、实现对交通连续流作用下桥梁耦合振动响应的预警监控。本发明通过计算随机连续流下桥梁动力响应，将所得冲击系数与规范值对比，提出预警指标，对桥梁结构损伤、性能评定、运营安全评估提供参考。

Description

一种交通连续流作用下桥梁耦合振动响应预警方法
技术领域
本发明属于交通连续流作用下桥梁耦合振动分析技术领域，具体涉及一种交通连续流作用下桥梁耦合振动响应预警方法。
背景技术
近年来，随着我国公路交通运输事业飞速发展，车辆类型呈现多样化，车辆数量、密度、载重不断增加，车辆行驶速度显著提升，公路桥梁作为公路建设中的咽喉，在设计上逐渐向轻型化、柔性化、大跨度化发展，这使得车桥耦合振动问题愈发明显，由此引发的安全事故更不可忽视。尤其在随机连续流作用下，交通荷载组成随机且复杂，桥上的车辆数量、车辆类型、车辆间距、行驶速度均具有很大的随机性，只研究既定荷载作用下的桥梁动力行为不够准确，与桥梁实际运营时所受随机连续流荷载影响存在较大差异。因此，现如今缺少一种针对随机车流下桥梁的耦合振动响应预警方法。
发明内容
本发明所要解决的技术问题在于针对上述现有技术中的不足，提供一种交通连续流作用下桥梁耦合振动响应预警方法，通过计算随机连续流下桥梁动力响应，将计算所得冲击系数与规范值对比，提出相关预警指标，对桥梁结构损伤、性能评定、运营安全评估提供一定参考，避免只通过既定荷载作用下的桥梁动力行为对桥梁耦合振动响应进行评估，便于推广使用。
为解决上述技术问题，本发明采用的技术方案是：一种交通连续流作用下桥梁耦合振动响应预警方法，其特征在于，该方法包括以下步骤：
步骤一、初始化桥梁交通流，以桥梁上无车辆为桥梁上交通流初始状态；
步骤二、实时采集桥梁上交通流，所述交通流包括车辆数量、车辆类型、车辆间距和车辆行驶速度；
步骤三、构建当前交通流状态下的车桥耦合矩阵并求解对应的桥梁动力响应增量
Figure BDA0003650170090000021
其中，
Figure BDA0003650170090000022
为当前交通流状态下车辆质量矩阵，MB为桥梁质量矩阵，
Figure BDA0003650170090000023
为当前交通流状态下车辆阻尼矩阵，CB为桥梁阻尼矩阵，
Figure BDA0003650170090000024
为当前交通流状态下车桥阻尼耦合矩阵；
Figure BDA0003650170090000025
为当前交通流状态下桥车阻尼耦合矩阵，
Figure BDA0003650170090000026
为当前交通流状态下车辆刚度矩阵，KB为桥梁刚度矩阵，
Figure BDA0003650170090000027
为当前交通流状态下车桥刚度耦合矩阵，
Figure BDA0003650170090000028
为当前交通流状态下桥车刚度耦合矩阵，
Figure BDA0003650170090000029
为当前交通流状态下车辆所受外力矩阵，
Figure BDA00036501700900000210
为当前交通流状态下桥梁所受外力矩阵；
Yi为当前交通流状态下车辆位移增量矩阵，
Figure BDA00036501700900000211
为当前交通流状态下车辆速度增量矩阵，
Figure BDA00036501700900000212
为当前交通流状态下车辆加速度增量矩阵；
Qi为当前交通流状态下桥梁位移增量矩阵，
Figure BDA00036501700900000213
为当前交通流状态下桥梁速度增量矩阵，
Figure BDA00036501700900000214
为当前交通流状态下桥梁加速度增量矩阵；
桥梁位移增量矩阵、桥梁速度增量矩阵和桥梁加速度增量矩阵构成桥梁动力响应增量；
步骤四、获取当前交通流状态下桥梁的动力响应结果：桥梁上交通流初始状态为0向量，累加当前交通流状态下桥梁动力响应增量，获取当前交通流状态下桥梁的动力响应结果；
步骤五、采集下一时刻桥梁上交通流、构建下一时刻交通流状态下的车桥耦合矩阵并求解对应的桥梁动力响应增量
Figure BDA00036501700900000215
其中，
Figure BDA00036501700900000216
为下一时刻交通流状态下车辆质量矩阵，
Figure BDA00036501700900000217
为下一时刻交通流状态下车辆阻尼矩阵，
Figure BDA00036501700900000218
为下一时刻交通流状态下车桥阻尼耦合矩阵；
Figure BDA00036501700900000219
为下一时刻交通流状态下桥车阻尼耦合矩阵，
Figure BDA0003650170090000031
为下一时刻交通流状态下车辆刚度矩阵，
Figure BDA0003650170090000032
为下一时刻交通流状态下车桥刚度耦合矩阵，
Figure BDA0003650170090000033
为下一时刻交通流状态下桥车刚度耦合矩阵，
Figure BDA0003650170090000034
为下一时刻交通流状态下车辆所受外力矩阵，
Figure BDA0003650170090000035
为下一时刻交通流状态下桥梁所受外力矩阵；
Yi+1为下一时刻交通流状态下车辆位移增量矩阵，
Figure BDA0003650170090000036
为下一时刻交通流状态下车辆速度增量矩阵，
Figure BDA0003650170090000037
为下一时刻交通流状态下车辆加速度增量矩阵；
Qi+1为下一时刻交通流状态下桥梁位移增量矩阵，
Figure BDA0003650170090000038
为下一时刻交通流状态下桥梁速度增量矩阵，
Figure BDA0003650170090000039
为下一时刻交通流状态下桥梁加速度增量矩阵；
步骤六、获取下一时刻交通流状态下桥梁的动力响应结果：
当下一时刻交通流状态下桥梁位移增量矩阵Qi+1和当前交通流状态下桥梁位移增量矩阵Qi的维度一致时，累加下一时刻交通流状态下桥梁的动力响应增量，获取下一时刻交通流状态下桥梁的动力响应结果；
当下一时刻交通流状态下桥梁位移增量矩阵Qi+1和当前交通流状态下桥梁位移增量矩阵Qi的维度不一致且存在新的车辆上桥时，对当前交通流状态下桥梁的动力响应结果对应的矩阵末位补0，再与下一时刻交通流状态下桥梁的动力响应增量进行累加，获取下一时刻交通流状态下桥梁的动力响应结果；
当下一时刻交通流状态下桥梁位移增量矩阵Qi+1和当前交通流状态下桥梁位移增量矩阵Qi的维度不一致且存在车辆下桥时，对当前交通流状态下桥梁的动力响应结果对应的矩阵首位补0，再与下一时刻交通流状态下桥梁的动力响应增量进行累加，获取下一时刻交通流状态下桥梁的动力响应结果；
步骤七、根据公式
Figure BDA00036501700900000310
计算下一时刻交通流状态下的桥梁冲击系数μ，其中，
Figure BDA00036501700900000311
为下一时刻交通流状态下桥梁位移响应结果最大值，
Figure BDA0003650170090000041
为下一时刻交通流状态下桥梁静态位移最大值；
当下一时刻交通流状态下的桥梁冲击系数μ大于桥梁规范冲击系数时，交通连续流作用下桥梁耦合振动响应发出预警信号；
当下一时刻交通流状态下的桥梁冲击系数μ不大于桥梁规范冲击系数时，交通连续流作用下桥梁耦合振动响应不发出预警信号；
步骤八、循环步骤五至步骤七，实现对交通连续流作用下桥梁耦合振动响应的预警监控。
上述的一种交通连续流作用下桥梁耦合振动响应预警方法，其特征在于：所述下一时刻交通流状态下桥梁静态位移最大值
Figure BDA0003650170090000044
通过在有限元软件中模拟下一时刻交通流对应的车辆在桥梁上低速通行时获取的桥梁位移最大值，所述低速为大于0且小于10km/h的车速。
上述的一种交通连续流作用下桥梁耦合振动响应预警方法，其特征在于：步骤七中，当桥梁的自振频率f≤1.5Hz时，桥梁规范冲击系数为0.05；
当桥梁的自振频率1.5Hz＜f≤14Hz时，桥梁规范冲击系数为0.1767Inf-0.0157；
当桥梁的自振频率f＞14Hz时，桥梁规范冲击系数为0.45。
上述的一种交通连续流作用下桥梁耦合振动响应预警方法，其特征在于：相邻的两个采集桥梁上交通流时刻间隔为1min～5min。
上述的一种交通连续流作用下桥梁耦合振动响应预警方法，其特征在于：步骤三和步骤五中，车辆所受外力矩阵和桥梁所受外力矩阵中对应的一组外力为相互作用力，即
Figure BDA0003650170090000042
其中，Fvn为桥上第n辆车辆所受外力，FBn为桥辆所受第n辆车辆带来的外力，kn为第n辆车辆的轮轴的车轮刚度，cn为第n辆车辆的轮轴的阻尼系数，Δn为相邻的两个采集桥梁上交通流时刻间隔内第n辆车辆的轮轴的车轮与桥梁间的相对位移，
Figure BDA0003650170090000043
为相邻的两个采集桥梁上交通流时刻间隔内第n辆车辆的轮轴的车轮与桥梁间的相对速度。
本发明与现有技术相比具有以下优点：
1、本发明通过构建当前交通流状态下的车桥耦合矩阵并求解对应的桥梁动力响应增量，获取当前交通流状态下桥梁的动力响应结果；再获取下一时刻交通流状态下桥梁的动力响应结果，利用变维度方式适用桥上车辆变化的动力响应结果计算，使用效果好，便于推广使用。
2、本发明通过计算随机连续流下桥梁动力响应，将计算所得冲击系数与规范值对比，提出相关预警指标，对桥梁结构损伤、性能评定、运营安全评估提供一定参考，可靠性高。
3、本发明方法步骤简单，迭代计算精度高，预警提示效果好，避免只通过既定荷载作用下的桥梁动力行为对桥梁耦合振动响应进行评估，便于推广使用。
综上所述，本发明通过计算随机连续流下桥梁动力响应，将计算所得冲击系数与规范值对比，提出相关预警指标，对桥梁结构损伤、性能评定、运营安全评估提供一定参考，避免只通过既定荷载作用下的桥梁动力行为对桥梁耦合振动响应进行评估，便于推广使用。
下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。
附图说明
图1为本发明的方法流程框图。
具体实施方式
如图1所示，本发明的一种交通连续流作用下桥梁耦合振动响应预警方法，包括以下步骤：
步骤一、初始化桥梁交通流，以桥梁上无车辆为桥梁上交通流初始状态；
步骤二、实时采集桥梁上交通流，所述交通流包括车辆数量、车辆类型、车辆间距和车辆行驶速度；
步骤三、构建当前交通流状态下的车桥耦合矩阵并求解对应的桥梁动力响应增量
Figure BDA0003650170090000061
其中，
Figure BDA0003650170090000062
为当前交通流状态下车辆质量矩阵，MB为桥梁质量矩阵，
Figure BDA0003650170090000063
为当前交通流状态下车辆阻尼矩阵，CB为桥梁阻尼矩阵，
Figure BDA0003650170090000064
为当前交通流状态下车桥阻尼耦合矩阵；
Figure BDA0003650170090000065
为当前交通流状态下桥车阻尼耦合矩阵，
Figure BDA0003650170090000066
为当前交通流状态下车辆刚度矩阵，KB为桥梁刚度矩阵，
Figure BDA0003650170090000067
为当前交通流状态下车桥刚度耦合矩阵，
Figure BDA0003650170090000068
为当前交通流状态下桥车刚度耦合矩阵，
Figure BDA0003650170090000069
为当前交通流状态下车辆所受外力矩阵，
Figure BDA00036501700900000610
为当前交通流状态下桥梁所受外力矩阵；
Yi为当前交通流状态下车辆位移增量矩阵，
Figure BDA00036501700900000611
为当前交通流状态下车辆速度增量矩阵，
Figure BDA00036501700900000612
为当前交通流状态下车辆加速度增量矩阵；
Qi为当前交通流状态下桥梁位移增量矩阵，
Figure BDA00036501700900000613
为当前交通流状态下桥梁速度增量矩阵，
Figure BDA00036501700900000614
为当前交通流状态下桥梁加速度增量矩阵；
桥梁位移增量矩阵、桥梁速度增量矩阵和桥梁加速度增量矩阵构成桥梁动力响应增量；
步骤四、获取当前交通流状态下桥梁的动力响应结果：桥梁上交通流初始状态为0向量，累加当前交通流状态下桥梁动力响应增量，获取当前交通流状态下桥梁的动力响应结果；
步骤五、采集下一时刻桥梁上交通流、构建下一时刻交通流状态下的车桥耦合矩阵并求解对应的桥梁动力响应增量
Figure BDA00036501700900000615
其中，
Figure BDA00036501700900000616
为下一时刻交通流状态下车辆质量矩阵，
Figure BDA00036501700900000617
为下一时刻交通流状态下车辆阻尼矩阵，
Figure BDA00036501700900000618
为下一时刻交通流状态下车桥阻尼耦合矩阵；
Figure BDA00036501700900000619
为下一时刻交通流状态下桥车阻尼耦合矩阵，
Figure BDA00036501700900000620
为下一时刻交通流状态下车辆刚度矩阵，
Figure BDA00036501700900000621
为下一时刻交通流状态下车桥刚度耦合矩阵，
Figure BDA00036501700900000622
为下一时刻交通流状态下桥车刚度耦合矩阵，
Figure BDA00036501700900000623
为下一时刻交通流状态下车辆所受外力矩阵，
Figure BDA00036501700900000624
为下一时刻交通流状态下桥梁所受外力矩阵；
Yi+1为下一时刻交通流状态下车辆位移增量矩阵，
Figure BDA0003650170090000071
为下一时刻交通流状态下车辆速度增量矩阵，
Figure BDA0003650170090000072
为下一时刻交通流状态下车辆加速度增量矩阵；
Qi+1为下一时刻交通流状态下桥梁位移增量矩阵，
Figure BDA0003650170090000073
为下一时刻交通流状态下桥梁速度增量矩阵，
Figure BDA0003650170090000074
为下一时刻交通流状态下桥梁加速度增量矩阵；
步骤六、获取下一时刻交通流状态下桥梁的动力响应结果：
当下一时刻交通流状态下桥梁位移增量矩阵Qi+1和当前交通流状态下桥梁位移增量矩阵Qi的维度一致时，累加下一时刻交通流状态下桥梁的动力响应增量，获取下一时刻交通流状态下桥梁的动力响应结果；
当下一时刻交通流状态下桥梁位移增量矩阵Qi+1和当前交通流状态下桥梁位移增量矩阵Qi的维度不一致且存在新的车辆上桥时，对当前交通流状态下桥梁的动力响应结果对应的矩阵末位补0，再与下一时刻交通流状态下桥梁的动力响应增量进行累加，获取下一时刻交通流状态下桥梁的动力响应结果；
当下一时刻交通流状态下桥梁位移增量矩阵Qi+1和当前交通流状态下桥梁位移增量矩阵Qi的维度不一致且存在车辆下桥时，对当前交通流状态下桥梁的动力响应结果对应的矩阵首位补0，再与下一时刻交通流状态下桥梁的动力响应增量进行累加，获取下一时刻交通流状态下桥梁的动力响应结果；
步骤七、根据公式
Figure BDA0003650170090000075
计算下一时刻交通流状态下的桥梁冲击系数μ，其中，
Figure BDA0003650170090000076
为下一时刻交通流状态下桥梁位移响应结果最大值，
Figure BDA0003650170090000077
为下一时刻交通流状态下桥梁静态位移最大值；
当下一时刻交通流状态下的桥梁冲击系数μ大于桥梁规范冲击系数时，交通连续流作用下桥梁耦合振动响应发出预警信号；
当下一时刻交通流状态下的桥梁冲击系数μ不大于桥梁规范冲击系数时，交通连续流作用下桥梁耦合振动响应不发出预警信号；
步骤八、循环步骤五至步骤七，实现对交通连续流作用下桥梁耦合振动响应的预警监控。
本实施例中，所述下一时刻交通流状态下桥梁静态位移最大值
Figure BDA0003650170090000083
通过在有限元软件中模拟下一时刻交通流对应的车辆在桥梁上低速通行时获取的桥梁位移最大值，所述低速为大于0且小于10km/h的车速。
本实施例中，步骤七中，当桥梁的自振频率f≤1.5Hz时，桥梁规范冲击系数为0.05；
当桥梁的自振频率1.5Hz＜f≤14Hz时，桥梁规范冲击系数为0.1767Inf-0.0157；
当桥梁的自振频率f＞14Hz时，桥梁规范冲击系数为0.45。
本实施例中，相邻的两个采集桥梁上交通流时刻间隔为1min～5min。
本实施例中，步骤三和步骤五中，车辆所受外力矩阵和桥梁所受外力矩阵中对应的一组外力为相互作用力，即
Figure BDA0003650170090000081
其中，Fvn为桥上第n辆车辆所受外力，FBn为桥辆所受第n辆车辆带来的外力，kn为第n辆车辆的轮轴的车轮刚度，cn为第n辆车辆的轮轴的阻尼系数，Δn为相邻的两个采集桥梁上交通流时刻间隔内第n辆车辆的轮轴的车轮与桥梁间的相对位移，
Figure BDA0003650170090000082
为相邻的两个采集桥梁上交通流时刻间隔内第n辆车辆的轮轴的车轮与桥梁间的相对速度。
本发明使用时，通过构建当前交通流状态下的车桥耦合矩阵并求解对应的桥梁动力响应增量，获取当前交通流状态下桥梁的动力响应结果；再获取下一时刻交通流状态下桥梁的动力响应结果，利用变维度方式适用桥上车辆变化的动力响应结果计算，使用效果好；通过计算随机连续流下桥梁动力响应，将计算所得冲击系数与规范值对比，提出相关预警指标，对桥梁结构损伤、性能评定、运营安全评估提供一定参考，可靠性高；迭代计算精度高，预警提示效果好，避免只通过既定荷载作用下的桥梁动力行为对桥梁耦合振动响应进行评估。
以上所述，仅是本发明的较佳实施例，并非对本发明作任何限制，凡是根据本发明技术实质对以上实施例所作的任何简单修改、变更以及等效结构变化，均仍属于本发明技术方案的保护范围内。

Claims (5)

1.一种交通连续流作用下桥梁耦合振动响应预警方法，其特征在于，该方法包括以下步骤：
步骤一、初始化桥梁交通流，以桥梁上无车辆为桥梁上交通流初始状态；
步骤二、实时采集桥梁上交通流，所述交通流包括车辆数量、车辆类型、车辆间距和车辆行驶速度；
步骤三、构建当前交通流状态下的车桥耦合矩阵并求解对应的桥梁动力响应增量
Figure FDA0003650170080000011
其中，
Figure FDA0003650170080000012
为当前交通流状态下车辆质量矩阵，MB为桥梁质量矩阵，
Figure FDA0003650170080000013
为当前交通流状态下车辆阻尼矩阵，CB为桥梁阻尼矩阵，
Figure FDA0003650170080000014
为当前交通流状态下车桥阻尼耦合矩阵；
Figure FDA0003650170080000015
为当前交通流状态下桥车阻尼耦合矩阵，
Figure FDA0003650170080000016
为当前交通流状态下车辆刚度矩阵，KB为桥梁刚度矩阵，
Figure FDA0003650170080000017
为当前交通流状态下车桥刚度耦合矩阵，
Figure FDA0003650170080000018
为当前交通流状态下桥车刚度耦合矩阵，
Figure FDA0003650170080000019
为当前交通流状态下车辆所受外力矩阵，
Figure FDA00036501700800000110
为当前交通流状态下桥梁所受外力矩阵；
Yi为当前交通流状态下车辆位移增量矩阵，
Figure FDA00036501700800000111
为当前交通流状态下车辆速度增量矩阵，
Figure FDA00036501700800000112
为当前交通流状态下车辆加速度增量矩阵；
Qi为当前交通流状态下桥梁位移增量矩阵，
Figure FDA00036501700800000113
为当前交通流状态下桥梁速度增量矩阵，
Figure FDA00036501700800000114
为当前交通流状态下桥梁加速度增量矩阵；
桥梁位移增量矩阵、桥梁速度增量矩阵和桥梁加速度增量矩阵构成桥梁动力响应增量；
步骤四、获取当前交通流状态下桥梁的动力响应结果：桥梁上交通流初始状态为0向量，累加当前交通流状态下桥梁动力响应增量，获取当前交通流状态下桥梁的动力响应结果；
步骤五、采集下一时刻桥梁上交通流、构建下一时刻交通流状态下的车桥耦合矩阵并求解对应的桥梁动力响应增量
Figure FDA0003650170080000021
其中，
Figure FDA0003650170080000022
为下一时刻交通流状态下车辆质量矩阵，
Figure FDA0003650170080000023
为下一时刻交通流状态下车辆阻尼矩阵，
Figure FDA0003650170080000024
为下一时刻交通流状态下车桥阻尼耦合矩阵；
Figure FDA0003650170080000025
为下一时刻交通流状态下桥车阻尼耦合矩阵，
Figure FDA0003650170080000026
为下一时刻交通流状态下车辆刚度矩阵，
Figure FDA0003650170080000027
为下一时刻交通流状态下车桥刚度耦合矩阵，
Figure FDA0003650170080000028
为下一时刻交通流状态下桥车刚度耦合矩阵，
Figure FDA0003650170080000029
为下一时刻交通流状态下车辆所受外力矩阵，
Figure FDA00036501700800000210
为下一时刻交通流状态下桥梁所受外力矩阵；
Yi+1为下一时刻交通流状态下车辆位移增量矩阵，
Figure FDA00036501700800000211
为下一时刻交通流状态下车辆速度增量矩阵，
Figure FDA00036501700800000212
为下一时刻交通流状态下车辆加速度增量矩阵；
Qi+1为下一时刻交通流状态下桥梁位移增量矩阵，
Figure FDA00036501700800000213
为下一时刻交通流状态下桥梁速度增量矩阵，
Figure FDA00036501700800000214
为下一时刻交通流状态下桥梁加速度增量矩阵；
步骤六、获取下一时刻交通流状态下桥梁的动力响应结果：
当下一时刻交通流状态下桥梁位移增量矩阵Qi+1和当前交通流状态下桥梁位移增量矩阵Qi的维度一致时，累加下一时刻交通流状态下桥梁的动力响应增量，获取下一时刻交通流状态下桥梁的动力响应结果；
当下一时刻交通流状态下桥梁位移增量矩阵Qi+1和当前交通流状态下桥梁位移增量矩阵Qi的维度不一致且存在新的车辆上桥时，对当前交通流状态下桥梁的动力响应结果对应的矩阵末位补0，再与下一时刻交通流状态下桥梁的动力响应增量进行累加，获取下一时刻交通流状态下桥梁的动力响应结果；
当下一时刻交通流状态下桥梁位移增量矩阵Qi+1和当前交通流状态下桥梁位移增量矩阵Qi的维度不一致且存在车辆下桥时，对当前交通流状态下桥梁的动力响应结果对应的矩阵首位补0，再与下一时刻交通流状态下桥梁的动力响应增量进行累加，获取下一时刻交通流状态下桥梁的动力响应结果；
步骤七、根据公式
Figure FDA0003650170080000031
计算下一时刻交通流状态下的桥梁冲击系数μ，其中，
Figure FDA0003650170080000032
为下一时刻交通流状态下桥梁位移响应结果最大值，
Figure FDA0003650170080000033
为下一时刻交通流状态下桥梁静态位移最大值；
当下一时刻交通流状态下的桥梁冲击系数μ大于桥梁规范冲击系数时，交通连续流作用下桥梁耦合振动响应发出预警信号；
当下一时刻交通流状态下的桥梁冲击系数μ不大于桥梁规范冲击系数时，交通连续流作用下桥梁耦合振动响应不发出预警信号；
步骤八、循环步骤五至步骤七，实现对交通连续流作用下桥梁耦合振动响应的预警监控。
2.按照权利要求1所述的一种交通连续流作用下桥梁耦合振动响应预警方法，其特征在于：所述下一时刻交通流状态下桥梁静态位移最大值
Figure FDA0003650170080000034
通过在有限元软件中模拟下一时刻交通流对应的车辆在桥梁上低速通行时获取的桥梁位移最大值，所述低速为大于0且小于10km/h的车速。
3.按照权利要求1所述的一种交通连续流作用下桥梁耦合振动响应预警方法，其特征在于：步骤七中，当桥梁的自振频率f≤1.5Hz时，桥梁规范冲击系数为0.05；
当桥梁的自振频率1.5Hz＜f≤14Hz时，桥梁规范冲击系数为0.1767Inf-0.0157；
当桥梁的自振频率f＞14Hz时，桥梁规范冲击系数为0.45。
4.按照权利要求1所述的一种交通连续流作用下桥梁耦合振动响应预警方法，其特征在于：相邻的两个采集桥梁上交通流时刻间隔为1min～5min。
5.按照权利要求4所述的一种交通连续流作用下桥梁耦合振动响应预警方法，其特征在于：步骤三和步骤五中，车辆所受外力矩阵和桥梁所受外力矩阵中对应的一组外力为相互作用力，即
Figure FDA0003650170080000041
其中，Fvn为桥上第n辆车辆所受外力，FBn为桥辆所受第n辆车辆带来的外力，kn为第n辆车辆的轮轴的车轮刚度，cn为第n辆车辆的轮轴的阻尼系数，Δn为相邻的两个采集桥梁上交通流时刻间隔内第n辆车辆的轮轴的车轮与桥梁间的相对位移，
Figure FDA0003650170080000042
为相邻的两个采集桥梁上交通流时刻间隔内第n辆车辆的轮轴的车轮与桥梁间的相对速度。
CN202210547758.1A 2022-05-18 2022-05-18 一种交通连续流作用下桥梁耦合振动响应预警方法 Active CN114913688B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210547758.1A CN114913688B (zh) 2022-05-18 2022-05-18 一种交通连续流作用下桥梁耦合振动响应预警方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210547758.1A CN114913688B (zh) 2022-05-18 2022-05-18 一种交通连续流作用下桥梁耦合振动响应预警方法

Publications (2)

Publication Number Publication Date
CN114913688A CN114913688A (zh) 2022-08-16
CN114913688B true CN114913688B (zh) 2023-02-14

Family

ID=82769056

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210547758.1A Active CN114913688B (zh) 2022-05-18 2022-05-18 一种交通连续流作用下桥梁耦合振动响应预警方法

Country Status (1)

Country Link
CN (1) CN114913688B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107357951A (zh) * 2017-05-22 2017-11-17 郑州大学 基于参数识别技术的桥梁冲击系数优化方法
CN110334371A (zh) * 2019-04-18 2019-10-15 朱思宇 一种基于有限元模型的车-桥耦合系统振动计算方法
CN110530589A (zh) * 2019-09-05 2019-12-03 哈尔滨工业大学 一种考虑连续车流耦合效应的梁式桥冲击系数检测方法
WO2021159569A1 (en) * 2020-02-10 2021-08-19 Dalian University Of Technology Method for engineering structural fiexibility identification considering the impact of noise
CN113392451A (zh) * 2021-06-09 2021-09-14 哈尔滨工业大学 基于车-桥梁耦合作用力修正的桥梁模型更新方法、系统、存储介质及设备

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107357951A (zh) * 2017-05-22 2017-11-17 郑州大学 基于参数识别技术的桥梁冲击系数优化方法
CN110334371A (zh) * 2019-04-18 2019-10-15 朱思宇 一种基于有限元模型的车-桥耦合系统振动计算方法
CN110530589A (zh) * 2019-09-05 2019-12-03 哈尔滨工业大学 一种考虑连续车流耦合效应的梁式桥冲击系数检测方法
WO2021159569A1 (en) * 2020-02-10 2021-08-19 Dalian University Of Technology Method for engineering structural fiexibility identification considering the impact of noise
CN113392451A (zh) * 2021-06-09 2021-09-14 哈尔滨工业大学 基于车-桥梁耦合作用力修正的桥梁模型更新方法、系统、存储介质及设备

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
基于主成分-逐步回归的大跨弯连续刚构桥冲击系数计算;韩智强等;《公路交通科技》;20220131;第39卷(第1期);全文 *
多点激励简支梁桥车桥耦合振动响应;韩智强等;《科学技术与工程》;20211231;第21卷(第18期);全文 *
大跨连续梁桥车桥耦合振动响应;韩智强等;《科学技术与工程》;20220430;第22卷(第11期);全文 *

Also Published As

Publication number Publication date
CN114913688A (zh) 2022-08-16

Similar Documents

Publication Publication Date Title
CN109612559B (zh) 基于分布式长标距光纤光栅传感器的桥梁式动态称重方法
CN101281117B (zh) 大跨度轨道交通桥梁损伤识别方法
CN104615888B (zh) 一种基于广义最小残差方法的桥梁移动车辆荷载识别方法
CN107389168A (zh) 一种用于桥梁的车辆移动载荷的识别方法
CN108573601B (zh) 一种基于wim数据的交通安全风险场构建方法
CN105923014B (zh) 一种基于证据推理规则的轨道高低不平顺幅值估计方法
CN104598931B (zh) 基于Kaczmarz代数迭代重建方法的桥梁移动车辆荷载识别方法
CN101739816A (zh) 交通车辆道路行驶安全分析方法
CN107657112A (zh) 一种轨道列车多级能量吸收的耐冲击优化设计方法
CN110398280A (zh) 一种测量车辆荷载的车载终端及测量、计算方法
CN112051048B (zh) 一种基于移动车辆动力作用的空心板桥铰缝快速评估方法
CN112765727B (zh) 一种轻卡驱动负荷耐久试验强化系数计算方法
CN111976731B (zh) 基于车辆频域响应的路面不平度识别方法
CN114913688B (zh) 一种交通连续流作用下桥梁耦合振动响应预警方法
CN106908128A (zh) 基于Landweber迭代法的桥面多轴移动荷载的识别方法
CN110530589B (zh) 一种考虑连续车流耦合效应的梁式桥冲击系数检测方法
CN114323510A (zh) 一种简支梁桥挠度快速检测方法
CN111666623B (zh) 双轴车辆缓慢通过连续梁桥跨中位移影响线快速提取方法
CN114023065A (zh) 一种利用视频解析数据智能诊断路口服务等级的算法
CN112326787B (zh) 一种基于专属测试车多点快速静态采集的梁式桥识别方法
CN111581723B (zh) 一种车辆缓慢通过简支梁桥跨中位移影响线快速提取方法
CN106768234A (zh) 基于经典行作用法的桥面多轴移动荷载的识别方法
CN110987499A (zh) 一种桥梁动载试验方法
CN102829968A (zh) 汽车波纹挠性节路谱采集方法及其采集系统
CN113190996B (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