CN113494951A - 计测方法、计测装置、计测系统以及存储介质 - Google Patents

计测方法、计测装置、计测系统以及存储介质 Download PDF

Info

Publication number
CN113494951A
CN113494951A CN202110282521.0A CN202110282521A CN113494951A CN 113494951 A CN113494951 A CN 113494951A CN 202110282521 A CN202110282521 A CN 202110282521A CN 113494951 A CN113494951 A CN 113494951A
Authority
CN
China
Prior art keywords
observation point
observation
deflection
vehicle
measurement
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.)
Granted
Application number
CN202110282521.0A
Other languages
English (en)
Other versions
CN113494951B (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.)
Seiko Epson Corp
Original Assignee
Seiko Epson Corp
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 Seiko Epson Corp filed Critical Seiko Epson Corp
Publication of CN113494951A publication Critical patent/CN113494951A/zh
Application granted granted Critical
Publication of CN113494951B publication Critical patent/CN113494951B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H1/00—Measuring characteristics of vibrations in solids by using direct conduction to the detector
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01G—WEIGHING
    • G01G19/00—Weighing apparatus or methods adapted for special purposes not provided for in the preceding groups
    • G01G19/02—Weighing apparatus or methods adapted for special purposes not provided for in the preceding groups for weighing wheeled or rolling bodies, e.g. vehicles
    • G01G19/03—Weighing apparatus or methods adapted for special purposes not provided for in the preceding groups for weighing wheeled or rolling bodies, e.g. vehicles for weighing during motion
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01G—WEIGHING
    • G01G3/00—Weighing apparatus characterised by the use of elastically-deformable members, e.g. spring balances
    • G01G3/12—Weighing apparatus characterised by the use of elastically-deformable members, e.g. spring balances wherein the weighing element is in the form of a solid body stressed by pressure or tension during weighing
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M5/00—Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings
    • G01M5/0008—Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings of bridges
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01G—WEIGHING
    • G01G19/00—Weighing apparatus or methods adapted for special purposes not provided for in the preceding groups
    • G01G19/02—Weighing apparatus or methods adapted for special purposes not provided for in the preceding groups for weighing wheeled or rolling bodies, e.g. vehicles
    • G01G19/022—Weighing apparatus or methods adapted for special purposes not provided for in the preceding groups for weighing wheeled or rolling bodies, e.g. vehicles for weighing wheeled or rolling bodies in motion
    • G01G19/024—Weighing apparatus or methods adapted for special purposes not provided for in the preceding groups for weighing wheeled or rolling bodies, e.g. vehicles for weighing wheeled or rolling bodies in motion using electrical weight-sensitive devices

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Length Measuring Devices With Unspecified Measuring Means (AREA)
  • Bridges Or Land Bridges (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明公开了计测方法、计测装置、计测系统以及存储介质。计测方法，包括：基于观测装置的观测信息，获取包括移动体的多个部位通过了结构物的第一观测点的时刻以及对于所述多个部位对第一观测点的作用的响应即物理量的第一观测点信息的步骤；基于所述观测信息，获取包括所述多个部位通过了第二观测点的时刻以及对于所述多个部位对第二观测点的作用的响应即物理量的第二观测点信息的步骤；基于第一观测点信息及第二观测点信息、规定的系数、以及结构物的挠曲的近似式，计算所述多个部位所引起的结构物的挠曲波形的步骤；以及将所述多个部位所引起的结构物的挠曲波形相加，计算移动体所引起的结构物的挠曲波形的步骤。

Description

计测方法、计测装置、计测系统以及存储介质
技术领域
本发明涉及计测方法、计测装置、计测系统以及存储介质。
背景技术
在专利文献1中记载了以下桥梁通过车辆监视系统：在进行桥梁的维护管理方面，通过桥梁的大型车辆的车轴重量是用于预测桥梁损伤的重要信息，并且为了测定该轴重，提出了由设置于桥梁的主梁的应变仪连续测定车辆通过时的应变值并计算轴重的方法(Weight In Motion：动态称重)，基于用配置于桥梁的主梁的应变仪计测出的应变波形来计测通过桥梁的车辆的车重。
详细地，在桥梁通过车辆监视系统中，配置应变仪，根据用应变仪计测出的应变波形来检测车轴的通过定时并计算车辆的轴间比率，对算出的轴间比率与根据轴间距离数据库中已登记的轴间距离算出的轴间比率进行比较，确定车辆的轴间距离、车速以及车型。另外，桥梁通过车辆监视系统配合车轴的通过定时而生成将基准轴重应变波形配置在时间轴上的应变波形，对基准轴重应变波形与用应变仪计测出的应变波形进行比较，来计算各轴的轴重。并且，桥梁通过车辆监视系统通过对各轴的轴重进行求和，从而计算车重。
专利文献1：日本特开2009-237805号公报
发明内容
但是，在专利文献1所记载的系统中，虽然能够通过使用应变波形和轴间距离数据库而不用计测桥梁的位移来计测车辆的车重，但是无法基于用应变仪计测出的应变波形来计算在桥梁等结构物上移动的车辆等移动体所引起的结构物的挠曲波形。
本发明的计测方法的一方式包括：第一观测点信息获取步骤，基于对沿着移动体在结构物上移动的第一方向排列的所述结构物的第一观测点和第二观测点中的所述第一观测点进行观测的第一观测装置的观测信息，获取包括所述移动体的多个部位分别通过了所述第一观测点的时刻以及对于所述多个部位分别对所述第一观测点的作用的响应即物理量的第一观测点信息；第二观测点信息获取步骤，基于观测所述第二观测点的第二观测装置的观测信息，获取包括所述多个部位分别通过了所述第二观测点的时刻以及对于所述多个部位分别对所述第二观测点的作用的响应即物理量的第二观测点信息；挠曲波形计算步骤，基于所述第一观测点信息及所述第二观测点信息、规定的系数以及所述结构物的挠曲的近似式，计算所述多个部位各自所引起的所述结构物的挠曲波形；以及移动体挠曲波形计算步骤，将在所述挠曲波形计算步骤中算出的所述多个部位各自所引起的所述结构物的挠曲波形相加，计算所述移动体所引起的所述结构物的挠曲波形。
在所述计测方法的一方式中也可以是，所述第一观测点被设定在所述结构物的第一端部，所述第二观测点被设定在所述结构物的与所述第一端部不同的第二端部。
在所述计测方法的一方式中也可以是，所述移动体是铁路车辆、汽车、路面电车、建筑车辆或者军用车辆，所述多个部位分别是车轴或者车轮。
在所述计测方法的一方式中也可以是，所述结构物的挠曲的近似式是以所述第一观测点与所述第二观测点的中央位置处的挠曲的最大振幅进行归一化后的式子。
在所述计测方法的一方式中也可以是，所述结构物的挠曲的近似式是基于所述结构物的结构模型的式子。
在所述计测方法的一方式中也可以是，所述结构模型是将两端进行支承而成的简支梁。
在所述计测方法的一方式中也可以是，所述结构物的挠曲的近似式是正弦波的半波长的波形的式子。
在所述计测方法的一方式中也可以是，所述结构物是桥梁的上部结构，所述上部结构是架设于相邻的桥台和桥墩、相邻的两个桥台、或者相邻的两个桥墩中的任意一方的结构，所述上部结构的两端部位于所述相邻的桥台和桥墩的位置、所述相邻的两个桥台的位置、或者所述相邻的两个桥墩的位置，所述桥梁是公路桥或者铁路桥。
在所述计测方法的一方式中也可以是，所述规定的系数是使所述第一观测点与所述第二观测点之间的所述结构物的部位的挠曲与所述结构物的所述部位的位移或者对该部位的载荷的关联近似的函数的系数。
在所述计测方法的一方式中也可以是，观测所述第一观测点的所述第一观测装置和观测所述第二观测点的所述第二观测装置是加速度传感器。
在所述计测方法的一方式中也可以是，观测所述第一观测点的所述第一观测装置和观测所述第二观测点的所述第二观测装置是冲击传感器、麦克风、应变仪或者测压元件。
在所述计测方法的一方式中也可以是，所述结构物是BWIM(Bridge Weigh inMotion：桥梁动态称重)系统发挥功能的结构。
在本发明的计测装置的一方式中，包括：第一观测点信息获取部，基于对沿着移动体在结构物上移动的第一方向排列的所述结构物的第一观测点和第二观测点中的所述第一观测点进行观测的第一观测装置的观测信息，获取包括所述移动体的多个部位分别通过了所述第一观测点的时刻以及对于所述多个部位分别对所述第一观测点的作用的响应即物理量的第一观测点信息；第二观测点信息获取部，基于观测所述第二观测点的第二观测装置的观测信息，获取包括所述多个部位分别通过了所述第二观测点的时刻以及对于所述多个部位分别对所述第二观测点的作用的响应即物理量的第二观测点信息；挠曲波形计算部，基于所述第一观测点信息及所述第二观测点信息、规定的系数以及所述结构物的挠曲的近似式，计算所述多个部位各自所引起的所述结构物的挠曲波形；以及移动体挠曲波形计算部，将所述挠曲波形计算部算出的所述多个部位各自所引起的所述结构物的挠曲波形相加，计算所述移动体所引起的所述结构物的挠曲波形。
在本发明的计测系统的一方式中，具备：所述计测装置的一方式；所述第一观测装置；以及所述第二观测装置。
在本发明的存储介质的一方式中，存储使计算机执行以下步骤的计测程序：第一观测点信息获取步骤，基于对沿着移动体在结构物上移动的第一方向排列的所述结构物的第一观测点和第二观测点中的所述第一观测点进行观测的第一观测装置的观测信息，获取包括所述移动体的多个部位分别通过了所述第一观测点的时刻以及对于所述多个部位分别对所述第一观测点的作用的响应即物理量的第一观测点信息；第二观测点信息获取步骤，基于观测所述第二观测点的第二观测装置的观测信息，获取包括所述多个部位分别通过了所述第二观测点的时刻以及对于所述多个部位分别对所述第二观测点的作用的响应即物理量的第二观测点信息；挠曲波形计算步骤，基于所述第一观测点信息及所述第二观测点信息、规定的系数以及所述结构物的挠曲的近似式，计算所述多个部位各自所引起的所述结构物的挠曲波形；以及移动体挠曲波形计算步骤，将在所述挠曲波形计算步骤中算出的所述多个部位各自所引起的所述结构物的挠曲波形相加，计算所述移动体所引起的所述结构物的挠曲波形。
附图说明
图1是示出计测系统的构成例的图。
图2是示出传感器和观测点的配置例的图。
图3是示出传感器和观测点的配置例的图。
图4是加速度传感器所检测的加速度的说明图。
图5是示出车轴信息的一例的图。
图6是示出传感器和观测点的配置例的图。
图7是示出传感器和观测点的配置例的图。
图8是示出针对观测点检测的加速度的一例的图。
图9是将图8的各时刻的加速度振幅转换为加速度强度的图。
图10是将图9的加速度强度用规定阈值进行了二值化的图。
图11是使退出时刻的图案相对于图10滑动了的图。
图12是桥梁的上部结构的结构模型的说明图。
图13是桥梁的上部结构的结构模型的说明图。
图14是示出归一化挠曲量的波形的一例的图。
图15是示出归一化挠曲量模型的一例的图。
图16是示出各车轴所引起的桥梁的挠曲波形的一例的图。
图17是示出车辆所引起的桥梁的挠曲波形的一例的图。
图18是绘制了载荷试验的结果的图。
图19是示出根据计测位移和挠曲波形算出的位移的一例的图。
图20是示出第一实施方式的计测方法的顺序的一例的流程图。
图21是示出挠曲波形计算步骤的顺序的一例的流程图。
图22是示出第一实施方式的计测装置的构成例的图。
图23是绘制了载荷试验的结果的图。
图24是示出根据基于车辆的载荷和挠曲波形算出的载荷的一例的图。
图25是示出第二实施方式的计测方法的顺序的一例的流程图。
图26是示出第二实施方式的计测装置的构成例的图。
图27是示出第三实施方式的归一化挠曲量的波形的一例的图。
附图标记说明
1…计测装置，2…服务器，4…通信网络，5…桥梁，6…车辆，7…上部结构，7a…桥面，7b…支承体，G…主梁，F…地板，8…下部结构，8a…桥墩，8b…桥台，10…计测系统，21…传感器，22…传感器，110…控制部，111…第一观测点信息获取部，112…第二观测点信息获取部，113…挠曲波形计算部，114…移动体挠曲波形计算部，115…位移计算部，116…系数值计算部，117…输出处理部，118…载荷计算部，120…第一通信部，130…存储部，131…计测程序，140…第二通信部，150…操作部。
具体实施方式
以下，使用附图详细说明本发明优选的实施方式。此外，以下说明的实施方式并非对权利要求书所记载的本发明的内容进行不当限定。另外，以下说明的全部构成不限于是本发明的必要构成要件。
1.第一实施方式
1-1.计测系统
以下列举结构物是桥梁的上部结构且移动体是车辆的情况为例来说明用于实现本实施方式的计测方法的计测系统。本实施方式的通过桥梁的车辆是铁路车辆、汽车、路面电车、建筑车辆、或者军用车辆等重量大且能够用BWIM(Bridge Weigh in Motion：桥梁动态称重)来计测的车辆。BWIM是将桥梁比作“秤”来计测桥梁的变形从而测定在桥梁上通行的车辆的重量、轴数等的技术。能够根据变形或应变等响应来分析通行车辆的重量的桥梁的上部结构为BWIM系统发挥功能的结构物，应用对桥梁的上部结构的作用与响应之间的物理过程的BWIM系统能够计测通行车辆的重量。
图1是示出本实施方式的计测系统的一例的图。如图1所示，本实施方式的计测系统10具有计测装置1、设置于桥梁5的上部结构7的至少一个传感器21、以及至少一个传感器22。另外，计测系统10也可以具有服务器2。
桥梁5包括上部结构7和下部结构8，上部结构7包括：由地板F、主梁G、未图示的横梁等形成的桥面7a、以及支承体7b。下部结构8包括桥墩8a和桥台8b。上部结构7是架设于相邻的桥台8b和桥墩8a、相邻的两个桥台8b、或者相邻的两个桥墩8a中的任意一者的结构。上部结构7的两端部位于相邻的桥台8b和桥墩8a的位置、相邻的两个桥台8b的位置、或者相邻的两个桥墩8a的位置。
计测装置1和各传感器21、22例如由未图示的电缆连接，经由CAN(ControllerArea Network：控制器局域网)等通信网络进行通信。或者，计测装置1与各传感器21、22可以经由无线网络进行通信。
例如，各传感器21输出表示移动体即车辆6向上部结构7进入所引起的冲击的数据，各传感器22输出表示车辆6从上部结构7退出所引起的冲击的数据。在本实施方式中，各传感器21、22是加速度传感器，例如，可以是水晶加速度传感器，也可以是MEMS(MicroElectro Mechanical Systems：微机电系统)加速度传感器。
在本实施方式中，各传感器21设置于上部结构7的长边方向的第一端部，各传感器22设置于上部结构7的长边方向的与第一端部不同的第二端部。
各传感器21检测当车辆6向上部结构7进入时产生的上部结构7的加速度，各传感器22检测当车辆6从上部结构7退出时产生的上部结构7的加速度。即，在本实施方式中，各传感器21是检测车辆6向上部结构7进入的加速度传感器，各传感器22是检测车辆6从上部结构7退出的加速度传感器。
计测装置1基于从各传感器21、22输出的加速度数据，计算车辆6行驶所引起的上部结构7的挠曲的位移。
计测装置1与服务器2例如能够经由便携电话的无线网络和互联网等通信网络4进行通信。计测装置1将车辆6在上部结构7上行驶的时刻、车辆6的行驶所引起的上部结构7的位移等信息发送到服务器2。服务器2可以将该信息存储到未图示的存储装置，例如，基于该信息进行超载车辆的监视、上部结构7的异常判定等处理。
此外，在本实施方式中，桥梁5是公路桥，例如是钢桥、梁桥、RC(Reinforced-Concrete：钢筋混凝土)桥等。
图2和图3是示出各传感器21、22设置于上部结构7的设置例的图。此外，图2是从上方观察上部结构7时的图，图3是将图2用A-A线或者B-B线截断后的剖视图。
如图2和图3所示，上部结构7具有移动体即车辆6能够移动的作为第一～第N路径的N个车道L1～LN以及K个主梁G1～GK。在此，N、K分别是1以上的整数。此外，在图2和图3的例子中，各主梁G1～GK的位置与各车道L1～LN的边界的位置一致，是N＝K-1，但各主梁G1～GK的位置无需与各车道L1～LN的边界的位置一致，也可以是N≠K-1。
在图2和图3的例子中，在上部结构7的长边方向的第一端部EA1中，在主梁G1～GK-1分别设置有传感器21，在上部结构7的长边方向的第二端部EA2中，在主梁G1～GK-1分别设置有传感器22。在图2和图3的例子中，是N＝K-1，在主梁GK未设置有传感器21、22，但也可以是，在主梁GK设置传感器21、22，而在主梁G1～GK-1中的任意一方不设置传感器21、22。或者也可以是N＝K，在主梁G1～GK分别设置传感器21、22。
此外，当将各传感器21、22设置于上部结构7的地板F时，有可能会被行驶车辆破坏，另外，有可能由于桥面7a的局部变形而使测定精度受到影响，因此在图2和图3的例子中，各传感器21、22设置于上部结构7的主梁G1～GK-1。
在本实施方式中，与N个传感器21相对应地分别设定有N个观测点P1～PN。观测点P1～PN是沿着与车辆6在上部结构7上移动的第一方向交叉的第二方向排列的上部结构7的N个观测点。在图2和图3的例子中，对于1以上且N以下的各整数j，观测点Pj在第一端部EA1中被设定于位于在主梁Gj上设置的传感器21的铅直上方的地板F的表面的位置。即，设置于主梁Gj的传感器21是观测观测点Pj的观测装置。观测观测点Pj的传感器21只要设置于能够检测由于车辆6的行驶而在观测点Pj产生的加速度的位置即可，但优选设置于离观测点Pj近的位置。这样，观测点P1～PN与N个传感器21为一对一的关系。
另外，在本实施方式中，与N个传感器22相对应地分别设定有N个观测点Q1～QN。观测点Q1～QN是沿着与车辆6在上部结构7上移动的第一方向交叉的第三方向排列的上部结构7的N个观测点。在图2和图3的例子中，对于1以上且N以下的各整数j，观测点Qj在第二端部EA2中被设定于位于在主梁Gj上设置的传感器22的铅直上方的地板F的表面的位置。即，设置于主梁Gj的传感器22是观测观测点Qj的观测装置。观测观测点Qj的传感器22只要设置于能够检测由于车辆6行驶而在观测点Qj产生的加速度的位置即可，但优选设置于离观测点Qj近的位置。这样，观测点Q1～QN与N个传感器22为一对一的关系。
在本实施方式中，N个观测点P1～PN分别与车道L1～LN相对应。同样地，N个观测点Q1～QN分别与车道L1～LN相对应。与车道Lj相对应地设定的观测点Pj和观测点Qj沿着车辆6在上部结构7上移动的第一方向排列。在图2和图3的例子中，第一方向是沿向上部结构7的车道L1～LN的X方向、即上部结构7的长边方向。另外，第二方向和第三方向是在车辆6所行驶的上部结构7的面内与X方向正交的Y方向、即上部结构7的宽度方向。不过，在车道L1～LN分别是曲线状等情况下，第二方向与第三方向可以不一致。另外，第二方向和第三方向也可以不与第一方向正交，例如，上部结构7的从车辆6进入的一侧的端部到观测点P1～PN的距离、上部结构7的从车辆6退出的一侧的端部到观测点Q1～QN的距离可以不同。此外，对于1以上且N以下的各整数j，观测点Pj是“第一观测点”的一例，观测点Qj是“第二观测点”的一例。
此外，传感器21、22的数量和设置位置不限于图2和图3所示的例子，能够实施各种变形。
计测装置1基于从各传感器21、22输出的加速度数据，获取与作为第一方向的X方向和作为第二方向及第三方向的Y方向分别交叉的第四方向的加速度。观测点P1～PN、Q1～QN通过冲击而向与X方向及Y方向正交的方向位移，因此，计测装置1为了准确地计算冲击的大小，优选获取与X方向及Y方向正交的第四方向、即地板F的法线方向的加速度。
图4是说明传感器21、22所检测的加速度的图。传感器21、22是检测在彼此正交的三个轴的各轴向上产生的加速度的加速度传感器。
为了检测由于车辆6向上部结构7进入而施加于观测点P1～PN的冲击，各传感器21被设置成作为三个检测轴的x轴、y轴、z轴中的一个轴的方向为与第一方向及第二方向交叉的方向。同样地，为了检测由于车辆6从上部结构7退出而施加于观测点Q1～QN的冲击，各传感器22被设置成作为三个检测轴的x轴、y轴、z轴中的一个轴的方向为与第一方向及第三方向交叉的方向。在图2和图3的例子中，第一方向是X方向，第二方向及第三方向是Y方向，因此，各传感器21、22被设置成一个轴的方向为与X方向及Y方向交叉的方向。观测点P1～PN、Q1～QN由于冲击而向与X方向及Y方向正交的方向位移，因此，为了准确地检测冲击的大小，理想的是，各传感器21、22的一个轴被设置成其方向是与X方向及Y方向正交的方向、即与地板F的法线方向一致。
不过，在将各传感器21、22设置于上部结构7时，也有设置场所为倾斜的情况。即使各传感器21、22的三个检测轴中的一个轴没有设置成与地板F的法线方向一致，由于计测装置1大致朝向法线方向，所以误差也小，能够忽略。另外，即使各传感器21、22的三个检测轴中的一个轴没有设置成与地板F的法线方向一致，计测装置1通过将x轴、y轴、z轴的加速度进行了合成的三轴合成加速度，也能够校正各传感器21、22的倾斜所引起的检测误差。另外，各传感器21、22也可以是检测至少在与铅直方向大致平行的方向上产生的加速度或者地板F的法线方向的加速度的单轴加速度传感器。
以下，对计测装置1所执行的本实施方式的计测方法的详细内容进行说明。
1-2.车轴信息的生成
在本实施方式中，计测装置1基于作为观测装置的N个传感器21的观测信息即加速度数据，获取包括移动体即车辆6的多个部位分别通过了观测点Pj的时刻以及对于该多个部位分别对观测点Pj的作用的响应即物理量的第一观测点信息。同样地，在本实施方式中，计测装置1基于作为观测装置的N个传感器22的观测信息即加速度数据，获取包括车辆6的多个部位分别通过了观测点Qj的时刻以及对于该多个部位分别对观测点Qj的作用的响应即物理量的第二观测点信息。其中，j是1以上且N以下的各整数。
在本实施方式中，认为基于车辆6所具备的多个车轴或者车轮的载荷施加于上部结构7，成为获取第一观测点信息及第二观测点信息的对象的多个部位分别是车轴或者车轮。以下，在本实施方式中，多个部位分别为车轴。
另外，在本实施方式中，作为加速度传感器的各传感器21检测多个车轴各自对观测点Pj的作用所引起的加速度。同样地，作为加速度传感器的各传感器22检测多个车轴各自对观测点Qj的作用所引起的加速度。
在本实施方式中，如图2所示，观测点P1～PN被设定于第一端部EA1，观测点Q1～QN被设定于第二端部EA2。因而，能够将车辆6的多个车轴分别通过了观测点Pj的时刻视为各车轴向上部结构7进入的时刻、更详细地向车道Lj进入的时刻。另外，能够将车辆6的多个车轴分别通过了观测点Qj的时刻视为各车轴从上部结构7退出的时刻、更详细地从车道Lj退出的时刻。
因而，在本实施方式中，第一观测点信息包括车辆6的各车轴向车道Lj进入的时刻以及作为对于各车轴进入车道Lj时的作用的响应即物理量的加速度强度。另外，第二观测点信息包括车辆6的各车轴从车道Lj退出的时刻以及作为对于各车轴从车道Lj退出时的作用的响应即物理量的加速度强度。
进而，车辆6的各车轴的进入与退出对应，因此能够将第一观测点信息和第二观测点信息分层，将包括第一观测点信息、第二观测点信息以及这些分层信息在内的信息称为车轴信息。
即，车轴信息除了第一观测点信息和第二观测点信息之外还包括每一车轴的、向车道Lj进入的时刻、进入时的加速度强度、从车道Lj退出的时刻以及退出时的加速度强度的对应信息、车辆6与每一车轴的该对应信息的对应信息。因而，通过车轴信息，按通过了上部结构7的每一车辆6来确定各车轴通过了的车道Lj、通过了观测点Pj、Qj的时刻以及通过时的加速度强度。
在图5中示出车轴信息的一例。在图5的例子中，第1行～第4行的信息是关于车辆编号为1的车辆6的信息。第1行的信息是关于车轴编号为1的前车轴的信息，第2行的信息是关于车轴编号为2的第2个车轴的信息，第3行的信息是关于车轴编号为3的第3个车轴的信息，第4行的信息是关于车轴编号为4的第4个车轴的信息。例如，对于车辆编号为1的车辆6的车轴编号为1的前车轴，第1行的对应信息示出了向车道L2进入的时刻是ti11，进入时的加速度强度是pai11，从车道L2退出的时刻是to11，退出时的加速度强度是pao11。
另外，第5行～第6行的信息是关于车辆编号为2的车辆6的信息。第5行的信息是关于车轴编号为1的前车轴的对应信息，第6行的信息是关于车轴编号为2的第2个车轴的对应信息。例如，对于车辆编号为2的车辆6的车轴编号为1的前车轴，第5行的对应信息示出了向车道L1进入的时刻是ti21，进入时的加速度强度是pai21，从车道L1退出的时刻是to21，退出时的加速度强度是pao21。
另外，第7行～第8行的信息是关于车辆编号为3的车辆6的信息。第7行的信息是关于车轴编号为1的前车轴的对应信息，第8行的信息是关于车轴编号为2的第2个车轴的对应信息。例如，对于车辆编号为3的车辆6的车轴编号为1的前车轴，第7行的对应信息示出了向车道L1进入的时刻是ti31，进入时的加速度强度是pai31，从车道L1退出的时刻是to31，退出时的加速度强度是pao31。
作为一例，在图6和图7中示出N＝2的情况下的各传感器21、22以及观测点P1、P2、Q1、Q2的配置例，在图6和图7所示的配置例的情况下，对计测装置1生成车轴信息的顺序进行说明
图6是从上部结构7的上方观察上部结构7时的图，图7是将图6用A-A线或者B-B线截断后的剖视图。在图6和图7的例子中，两个传感器21在上部结构7的第一端部EA1分别设置于主梁G1、G3，两个传感器22在上部结构7的第二端部EA2分别设置于主梁G1、G3。另外，与车道L1对应的观测点P1、Q1分别设定于位于在主梁G1上设置的传感器21、22的铅直上方的地板F的表面的位置，与车道L2对应的观测点P2、Q2分别设定于位于在主梁G3上设置的传感器21、22的铅直上方的地板F的表面的位置。设置于主梁G1的传感器21观测观测点P1，设置于主梁G3的传感器21观测观测点P2。另外，设置于主梁G1的传感器22观测观测点Q1，设置于主梁G3的传感器22观测观测点Q2。计测装置1为了生成车轴信息，而将各传感器21、22检测出的各时刻的加速度转换为振幅，从而获取加速度强度。
图8是示出4轴车辆6在车道L2上行驶的情况下针对观测点P1、P2、Q1、Q2检测的加速度的一例的图。另外，图9是将图8的各时刻的加速度振幅转换为加速度强度的图。在图8和图9的例子中，车辆6正在车道L2上行驶，因此，在车辆6的4个车轴分别通过观测点P2、Q2的时刻获取了较大的加速度强度。在4个车轴分别通过观测点P2的时刻所获取的加速度强度包括在第一观测点信息中。另外，在4个车轴分别通过观测点Q2的时刻所获取的加速度强度包括在第二观测点信息中。
并且，计测装置1将所获取的加速度强度超过规定阈值的时刻作为各车轴从前车轴起按顺序通过了观测点P2、Q2的时刻、即作为各车轴向车道L2进入的时刻和从车道L2退出的时刻来获取。
图10是将图9的加速度强度用规定阈值进行二值化后的图。在图10的例子中，获取了4个车轴分别向车道L2进入的时刻和从车道L2退出的时刻。4个车轴分别向车道L2进入的时刻包括在第一观测点信息中。另外，4个车轴分别从车道L2退出的时刻包括在第二观测点信息中。
进而，计测装置1对4个车轴分别向车道L2进入的时刻的图案1和4个车轴分别从车道L2退出的时刻的图案2进行比较，判定该两个图案是否是基于同一车辆6通过而形成的图案。4个车轴的间隔不发生变化，因此若车辆6在上部结构7上行驶的速度固定，则图案1、2一致。例如，为了使前车轴的进入时刻与退出时刻一致，计测装置1使图案1、2中的任意一方的时刻滑动，在第2个～第4个各车轴的进入时刻与退出时刻之差是规定阈值以下的情况下，判定为图案1、2是基于同一车辆6通过而形成的图案，在该差大于规定阈值的情况下，判定为图案1、2是基于2辆车辆6通过而形成的图案。此外，计测装置1在2辆车辆6以相同的速度在一个车道上连行的情况下，为了不将在先的车辆6的多个车轴和在后的车辆6的多个车轴误判定为1辆车辆6的车轴，在两个连续的车轴的进入时刻或者退出时刻的间隔是规定以上的时间差的情况下，将该两个车轴的进入时刻和退出时刻按两个车辆6分开即可。
图11是为了使前车轴的进入时刻与退出时刻一致而使4个车轴分别从车道L2退出的时刻的图案2相对于图10滑动后的图。此外，图11相对于图10放大了横轴方向。在图11的例子中，4个车轴分别向车道L2进入的时刻的图案1与4个车轴分别从车道L2退出的时刻的图案2大致一致，判定为图案1、2是基于同一车辆6通过而形成的图案。
并且，计测装置1将图10所示的4个向车道L2进入的时刻、图9所示的观测点P2的4个加速度强度的峰值、图10所示的4个从车道L2退出的时刻、以及图9所示的观测点Q2的4个加速度强度的峰值从前车轴起按顺序建立对应关系，从而获取前车轴的对应信息、第2个车轴的对应信息、第3个车轴的对应信息以及第4个车轴的对应信息。进而，计测装置1获取将在车道L2上行驶的车辆6与4个车轴的对应信息建立了对应关系的对应信息。这些信息与第一观测点信息及第二观测点信息一起包括在车轴信息中。
计测装置1能够通过车轴信息而针对在上部结构7的车道Lj上通过了的任意的车辆6来确定该车辆6的各车轴向观测点Pj进入的时刻、基于各车轴的观测点Pj的加速度强度、各车轴从观测点Qj退出的时刻、以及基于各车轴的观测点Qj的加速度强度。
1-3.挠曲波形的生成
在本实施方式中，在桥梁5的上部结构7中，考虑是连续配置1个或多个由地板F和主梁G1～GK等构成的桥面7a的构成，计测装置1计算一个桥面7a的位移作为长边方向的中央部的位移。施加于上部结构7的载荷从上部结构7的一端向另一端移动。此时，能够使用载荷在上部结构7上的位置和载荷量来表示上部结构7的中央部的位移即挠曲量。在本实施方式中，为了将车辆6的车轴在上部结构7上移动时的挠曲变形表示为1点载荷在梁上的移动所引起的挠曲量的轨迹，而考虑图12所示的结构模型，在该结构模型中，计算中央部的挠曲量。在图12中，P是载荷。a是距离车辆6进入侧的上部结构7的端部的载荷位置。b是距离车辆6退出侧的上部结构7的端部的载荷位置。l是上部结构7的两端之间的距离。图12所示的结构模型是以两端为支点支撑着两端的简支梁。
在图12所示的结构模型中，当以车辆6进入侧的上部结构7的端部的位置为零点并将挠曲量的观测位置设为x时，简支梁的弯矩M用式(1)表示。
[式1]
Figure BDA0002979166730000151
在式(1)中，函数Ha由式(2)定义。
[式2]
Figure BDA0002979166730000161
将式(1)变形并得到式(3)。
[式3]
Figure BDA0002979166730000162
另一方面，弯矩M用式(4)表示。在式(4)中，θ是角度，I是二次力矩，E是杨氏模量。
[式4]
Figure BDA0002979166730000163
将式(4)代入式(3)，得到式(5)。
[式5]
Figure BDA0002979166730000164
计算将式(5)对于观测位置x进行积分的式(6)，得到式(7)。在式(7)中，C1是积分常量。
[式6]
Figure BDA0002979166730000165
[式7]
Figure BDA0002979166730000166
进而，计算将式(7)对于观测位置x进行积分的式(8)，得到式(9)。在式(9)中，C2是积分常量。
[式8]
Figure BDA0002979166730000167
[式9]
Figure BDA0002979166730000171
在式(9)中，θx表示挠曲量，将θx置换为挠曲量w而得到式(10)。
[式10]
Figure BDA0002979166730000172
根据图12，由于b＝l-a，因此，式(10)可以变形为式(11)。
[式11]
Figure BDA0002979166730000173
设为x＝0且挠曲量w＝0，根据x≤a，Ha＝0，因此当将x＝w＝Ha＝0代入式(11)并整理时，得到式(12)。
[式12]
C2＝0…(12)
另外，设为x＝l且挠曲量w＝0，根据x＞a，Ha＝1，因此当将x＝l、w＝0、Ha＝1代入式(11)并整理时，得到式(13)。
[式13]
Figure BDA0002979166730000174
将b＝l-a代入式(13)，得到式(14)。
[式14]
Figure BDA0002979166730000175
将式(12)的积分常量C1及式(13)的积分常量C2代入式(10)，得到式(15)。
[式15]
Figure BDA0002979166730000176
将式(15)变形，载荷P施加于位置a时的观测位置x的挠曲量w用式(16)表示。
[式16]
Figure BDA0002979166730000181
在图13中，示出在挠曲量的观测位置x固定于简支梁的中央位置的条件下、即当x＝l/2时，载荷P从简支梁的一端向另一端移动的样子。
当载荷位置a位于比观测位置x＝l/2更靠左侧时，根据x＞a，Ha＝1，因此将x＝l/2、Ha＝1代入式(16)，得到式(17)。
[式17]
Figure BDA0002979166730000182
当将l＝a+b代入式(17)并整理时，得到式(18)。
[式18]
Figure BDA0002979166730000183
将a+b＝l代入式(18)，载荷P的位置位于比中央的观测位置x＝l/2更靠左侧的情况下的观测位置x的挠曲量wL如式(19)所示。
[式19]
Figure BDA0002979166730000184
另一方面，当载荷位置a位于比观测位置x＝l/2更靠右侧时，根据x≤a，Ha＝0，因此将x＝l/2、Ha＝0代入式(16)，得到式(20)。
[式20]
Figure BDA0002979166730000185
当将l＝a+b代入式(20)并整理时，载荷P的位置位于比中央的观测位置x＝l/2更靠右侧的情况下的观测位置x的挠曲量wR如式(21)所示。
[式21]
Figure BDA0002979166730000191
另外，当载荷位置a与观测位置x＝l/2相同时，根据x≤a，Ha＝0，因此当将Ha＝0、a＝b＝l/2代入式(16)并整理时，得到式(22)。
[式22]
Figure BDA0002979166730000192
进而，当将a＝l/2代入式(22)时，载荷P的位置位于与中央的观测位置相同的情况下的观测位置x的挠曲量w如式(23)所示。
[式23]
Figure BDA0002979166730000193
关于两端支点的简支梁，在载荷P位于中央的情况下为最大挠曲位移，因此根据式(23)，最大挠曲量wmax用式(24)表示。
[式24]
Figure BDA0002979166730000194
当载荷P的位置位于比中央的观测位置x＝l/2更靠左侧的情况下的观测位置x的挠曲量wL除以最大挠曲量wmax而以最大挠曲量wmax进行归一化时，根据式(19)和式(24)，得到式(25)。
[式25]
Figure BDA0002979166730000195
在式(25)中，当设置为a/l＝r时，得到式(26)。
[式26]
Figure BDA0002979166730000196
另一方面，当载荷P的位置位于比中央的观测位置x＝l/2更靠右侧的情况下的观测位置x的挠曲量wR除以最大挠曲量wmax而以最大挠曲量wmax进行归一化时，根据式(21)及式(24)，得到式(27)。
[式27]
Figure BDA0002979166730000201
其中，根据a/l＝r及a+b＝l，b＝l×(1-r)，因此将b＝l×(1-r)代入式(27)，得到式(28)。
[式28]
Figure BDA0002979166730000202
总结式(25)、式(27)，以载荷P在简支梁上移动时在中央部观测的最大挠曲量进行归一化后的归一化挠曲量wstd用式(29)来表示。
[式29]
Figure BDA0002979166730000203
在式(29)中，r＝a/l、1-r＝b/l表示载荷P的位置相对于简支梁的支点间的距离l的比例，如式(30)来定义变量R。
[式30]
Figure BDA0002979166730000204
使用式(30)，将式(29)置换为式(31)。
[式31]
wstd＝3R-4R3…(31)
式(30)和式(31)示出在观测位置位于简支梁的中央的情况下，在载荷P的位置比中央更靠右侧和更靠左侧的位置处挠曲量为对称。
在图14中，示出观测位置x＝l/2的情况下的归一化挠曲量wstd的波形的一例。在图14中，横轴是载荷P的位置，纵轴是归一化挠曲量wstd。在图14的例子中，简支梁的支点间的距离l＝1。
前述的车轴信息所包括的是，车辆6的各车轴向车道Lj进入的时刻和从车道Lj退出的时刻、即车辆6分别通过了上部结构7的两端的位置的时刻，因此使上部结构7的两端的位置与车轴的进入时刻及退出时刻对应而将载荷位置a、b置换为时间。不过，设为车辆6的速度大致固定，位置与时刻大致成比例。
当使上部结构7的左端的载荷位置与进入时刻ti对应，使上部结构7的右端的载荷位置与退出时刻to对应时，距离左端的载荷位置a被置换为自进入时刻ti起的经过时刻tp。经过时刻tp用式(32)表示。
[式32]
tp＝t-ti…(32)
另外，支点间的距离l被置换为从进入时刻ti到退出时刻to的时间ts。时间ts用式(33)表示。
[式33]
ts＝to-ti…(33)
车辆6的速度是固定的，因此，载荷位置a位于上部结构7的中央的时刻tc用式(34)表示。
[式34]
Figure BDA0002979166730000211
如上所示，将位置置换为时间，载荷P的位置成为式(35)和式(36)。
[式35]
Figure BDA0002979166730000212
[式36]
Figure BDA0002979166730000221
将式(35)及式(36)代入式(29)，被置换为时间的归一化挠曲量wstd用式(37)表示。
[式37]
Figure BDA0002979166730000222
或者，根据式(30)及式(31)，将变量R置换为时间，以最大振幅进行归一化后的归一化挠曲量wstd用式(38)表示。
[式38]
Figure BDA0002979166730000223
考虑将时间经过与归一化挠曲量的关系作为观测数据进行处理，将归一化挠曲量wstd置换为两端支点的简支梁上的单一集中载荷的移动所引起的梁中央的观测位置的归一化挠曲量模型wstd(t)，式(38)成为式(39)。式(39)是作为结构物的上部结构7的挠曲的近似式，是基于上部结构7的结构模型的式子。具体地，式(39)是以车辆6所移动的车道Lj中的观测点Pj与观测点Qj的中央位置处的挠曲的最大振幅被归一化、且最大值成为1的式子。
[式39]
Figure BDA0002979166730000231
该归一化挠曲量模型wstd(t)所需的时间信息是从前述的车轴信息得到的。关于归一化挠曲量模型wstd(t)，在上部结构7的中央位置处成为最大挠曲量wmax，因此得到式(40)。
[式40]
Figure BDA0002979166730000232
另外，前述的式(23)所示的挠曲量w是载荷P的位置与中央的观测位置相同的情况下的观测位置x＝l/2的挠曲量，与最大挠曲量wmax一致，因此得到式(41)。
[式41]
Figure BDA0002979166730000233
在图15中，示出归一化挠曲量模型wstd(t)的一例。在图15的例子中，进入时刻ti＝4，退出时刻to＝6，在时刻tc＝(ti+to)/2＝5时，归一化挠曲量模型wstd(t)在上部结构7的中央位置处成为最大挠曲量wmax＝1。
假设对于作为结构物的上部结构7而言，BWIM(Bridge Weigh in Motion)系统发挥功能，认为进行与以两端为支点的简支梁近似的变形。另外，由于作为移动体的车辆6从上部结构7的一个端部以大致恒定速度通过而向另一个端部移动，所以上部结构7的中间部和上部结构7的端部受到相同的载荷的作用，因此能够认为所观测的上部结构7的位移与根据车轴信息得到的车轴的加速度强度ap近似地成比例。
将比例系数设为根据车轴信息得到的车轴的加速度强度ap与规定的系数p的乘积，通过式(42)，得到各车轴所引起的上部结构7的挠曲波形H(t)。此外，加速度强度ap可以是车轴信息所包括的进入时的加速度强度，也可以是退出时的加速度强度，还可以是进入时的加速度强度与退出时的加速度强度的平均值等统计值。
[式42]
H(t)＝papwstd(t)…(4Z)
将式(39)代入式(42)，挠曲波形H(t)用式(43)表示。
[式43]
Figure BDA0002979166730000241
至此，设为对上部结构7施加单一载荷P，但由于对车辆6所行驶的车道Lj施加基于车辆6的各车轴的载荷，因此，式(43)可以如式(44)所示置换为挠曲波形Hjk(t)。在式(44)中，k是表示车轴编号的整数，j是表示车道编号的整数。如式(44)所示，挠曲波形Hjk(t)与规定的系数p与加速度强度apjk的乘积成比例。
[式44]
Figure BDA0002979166730000242
在图16中，示出在车道Lj上行驶的车辆6所包括的各车轴所引起的上部结构7的挠曲波形的一例。在图16的例子中，车辆6是4轴车辆，示出了4个挠曲波形Hj1(t)、Hj2(t)、Hj3(t)、Hj4(t)。在图16的例子中，基于前车轴和第2个车轴的载荷相对小，基于第3个和第4个车轴的载荷相对大，因此，挠曲波形Hj1(t)、Hj2(t)的最大振幅相对小，挠曲波形Hj3(t)、Hj4(t)的最大振幅相对大。
如式(45)所示，在车道Lj上行驶的车辆6所引起的上部结构7的挠曲波形CPjm(t)是将各车轴所引起的上部结构7的挠曲波形Hjk(t)相加而得到的。在式(45)中，m是表示车辆编号的整数，k是表示车轴编号的整数，j是表示车道编号的整数。
[式45]
Figure BDA0002979166730000251
在图17中，示出将图16所示的4个挠曲波形Hj1(t)、Hj2(t)、Hj3(t)、Hj4(t)相加而得到的、车辆编号为m的车辆6所引起的上部结构7的挠曲波形CPjm(t)。
1-4.位移的计算
将式(45)所示的挠曲波形CPjm(t)与观测到的位移CU(t)的关联用多项式进行近似。例如，如式(46)所示，将位移CU(t)用挠曲波形CPjm(t)的1次式进行近似。在式(46)中，s是1次系数，i是0次系数。
[式46]
Figure BDA0002979166730000252
1次系数s和0次系数i例如是通过基于多个车辆的载荷试验而算出的。例如，在车道Lj的中央设置位移仪，使多个车辆分别单独地在车道Lj上行驶，计测装置1生成车轴信息，并且获取由位移仪计测到的位移。然后，计测装置1将计测到的位移的最大值设为CUmax，将根据车轴信息并通过式(45)得到的挠曲波形CPjm(t)的最大值设为CPjm-max绘制为坐标图，求出近似直线的1次系数scu和0次系数icu。
图18是绘制出基于6辆车辆的载荷试验的结果的图。在图18中，横轴是车辆所引起的上部结构7的挠曲波形的最大值CPjm-max，纵轴是计测位移的最大值CUmax。在图18中，6个点按直线状配置，相对于6个点的近似直线用虚线示出。在图18的例子中，该近似直线的1次系数scu是3084.435944，0次系数icu是0.229180174。
计测装置1使用1次系数scu和0次系数icu、以及根据未知的车辆6的车轴信息并通过式(45)得到的挠曲波形CPjm(t)，通过式(47)计算车道Lj的中央的位移CUest(t)。
[式47]
CUest(t)＝scuCPjm(t)+icu…(47)
在图19中，示出根据计测位移CU(t)和挠曲波形CPjm(t)算出的位移CUest(t)的一例。在图19中，实线示出计测位移CU(t)，虚线示出根据挠曲波形CPjm(t)算出的位移CUest(t)。如图19所示，根据挠曲波形CPjm(t)算出的位移CUest(t)近似于计测位移CU(t)。因而，计测装置1能够不用计测车道Lj的中央的位移而通过式(47)算出在车道Lj上行驶的未知的车辆6所引起的中央的位移。
在式(47)中，0次系数icu是小的值，将icu＝0代入式(47)，从式(44)及式(45)得到式(48)。
[式48]
Figure BDA0002979166730000261
根据式(48)，1次系数scu能够与规定的系数p进行替换，因此，规定的系数p是具有与1次系数scu同样的功能的系数。即，规定的系数p是使观测点Pj与观测点Qj之间的上部结构7的部位的挠曲与上部结构7的该部位的位移的关联近似的函数的系数。
1-5.计测方法
图20是示出第一实施方式的计测方法的顺序的一例的流程图。在本实施方式中，计测装置1执行图20所示的顺序。
如图20所示，首先，计测装置1基于观测观测点P1～PN的N个传感器21的观测信息，获取包括车辆6的多个车轴通过了观测点P1～PN中的任意一个的时刻以及作为对于多个车轴分别对观测点P1～PN中的任意一个的作用的响应即物理量的加速度强度的第一观测点信息(步骤S1)。如前所述，N个传感器21分别是加速度传感器，N个传感器21的观测信息是产生于观测点P1～PN的加速度的检测信息。计测装置1基于N个传感器21分别检测出的加速度，获取第一观测点信息。该步骤S1是第一观测点信息获取步骤。
接着，计测装置1基于观测观测点Q1～QN的N个传感器22的观测信息，获取包括车辆6的多个车轴通过了观测点Q1～QN中的任意一个的时刻以及作为对于多个车轴分别对观测点Q1～QN中的任意一个的作用的响应即物理量的加速度强度的第二观测点信息(步骤S2)。如前所述，N个传感器22分别是加速度传感器，N个传感器22的观测信息是产生于观测点Q1～QN的加速度的检测信息。计测装置1基于N个传感器22分别检测出的加速度，获取第二观测点信息。该步骤S2是第二观测点信息获取步骤。
接着，计测装置1基于在步骤S1中获取的第一观测点信息和在步骤S2中获取的第二观测点信息、规定的系数p、以及上部结构7的挠曲的近似式，计算多个车轴各自所引起的上部结构7的挠曲波形Hjk(t)(步骤S3)。具体地，计测装置1使用第一观测点信息和第二观测点信息生成前述的车轴信息，使用车轴信息和规定的系数p，通过前述的式(44)计算在各车道Lj上行驶的车辆6的各车轴所引起的上部结构7的挠曲波形Hjk(t)。该步骤S3是挠曲波形计算步骤。
接着，计测装置1通过前述的式(45)，将在步骤S3中算出的车辆6的多个车轴各自所引起的上部结构7的挠曲波形Hjk(t)相加，计算车辆6所引起的上部结构7的挠曲波形CPjm(t)(步骤S4)。该步骤S4是移动体挠曲波形计算步骤。
接着，计测装置1通过前述的式(47)，基于在步骤S4中算出的车辆6所引起的上部结构7的挠曲波形CPjm(t)，计算上部结构7的位移CUest(t)(步骤S5)。该步骤S5是位移计算步骤。
接着，计测装置1将在步骤S5中算出的上部结构7的位移CUest(t)输出到服务器2(步骤S6)。该步骤S6是输出步骤。
计测装置1反复进行步骤S1～S6的处理，直至将计测结束为止(步骤S7的N)。
图21是示出作为图20的步骤S3的挠曲波形计算步骤的顺序的一例的流程图。
如图21所示，首先，计测装置1将整数j设定为1(步骤S30)，使用第一观测点信息及第二观测点信息，对各车轴向车道Lj进入的时刻的图案1与各车轴从车道Lj退出的时刻的图案2进行比较(步骤S31)。
然后，计测装置1在图案1所包括的各车轴的进入时刻与图案2所包括的各车轴的退出时刻之差为阈值以下的情况下(步骤S32的Y)，将图案1所包括的各车轴的进入时刻及加速度强度和图案2所包括的各车轴的退出时刻及加速度强度与1辆车辆6建立对应关系地生成车轴信息(步骤S33)。
另外，计测装置1在图案1所包括的各车轴的进入时刻与图案2所包括的各车轴的退出时刻之差大于阈值的情况下(步骤S32的N)，不进行步骤S33的处理。
计测装置1在整数j不是N的情况下(步骤S34的N)，对整数j加上1(步骤S35)，反复进行步骤S31～S33的处理。
并且，当整数j成为N时(步骤S34的Y)，计测装置1将整数j设定为1(步骤S36)，使用在步骤S33中生成的车轴信息和规定的系数p，对于在车道Lj上行驶的各车辆6，计算各车轴所引起的上部结构7的挠曲波形Hjk(t)(步骤S37)。
计测装置1在整数j不是N的情况下(步骤S38的N)，对整数j加上1(步骤S39)，反复进行步骤S37的处理。
并且，当整数j成为N时(步骤S38的Y)，计测装置1结束挠曲波形计算步骤的处理。
1-6.计测装置的构成
图22是示出第一实施方式的计测装置1的构成例的图。如图22所示，计测装置1具有控制部110、第一通信部120、存储部130、第二通信部140以及操作部150。
控制部110基于从设置于上部结构7的各传感器21、22输出的加速度数据，计算车辆6在上部结构7上行驶的时刻、上部结构7的位移等。
第一通信部120从各传感器21、22接收加速度数据。从各传感器21、22输出的加速度数据例如是数字信号。第一通信部120将从各传感器21、22接收到的加速度数据输出到控制部110。
存储部130是存储控制部110用于进行计算处理、控制处理的程序、数据等的存储器。另外，存储部130存储控制部110用于实现规定的应用功能的程序、数据等。存储部130例如由ROM(Read Only Memory：只读存储器)、闪存ROM、RAM(Random Access Memory：随机存取存储器)等各种IC(Integrated Circuit：集成电路)存储器、硬盘、存储卡等记录介质等构成。
存储部130可以包括作为能够通过计算机读取的装置、介质的非易失性信息存储装置，各种程序、数据等存储到该信息存储装置。信息存储装置可以是光盘DVD、CD等光盘、硬盘驱动器或者卡式存储器、ROM等各种存储器等。另外，控制部110也可以经由通信网络4接收各种程序、数据等并将其存储到存储部130。
第二通信部140将控制部110的计算结果等信息经由通信网络4发送到服务器2。
操作部150进行获取来自用户的操作数据并将其发送到控制部110的处理。
控制部110具备第一观测点信息获取部111、第二观测点信息获取部112、挠曲波形计算部113、移动体挠曲波形计算部114、位移计算部115、系数值计算部116以及输出处理部117。
第一观测点信息获取部111进行基于观测观测点P1～PN的N个传感器21的观测信息来获取包括车辆6的多个车轴通过了观测点P1～PN中的任意一个的时刻以及作为对于多个车轴分别对观测点P1～PN中的任意一个的作用的响应即物理量的加速度强度的第一观测点信息的处理。即，第一观测点信息获取部111进行图20的第一观测点信息获取步骤的处理。第一观测点信息获取部111所获取的第一观测点信息被存储到存储部130中。
第二观测点信息获取部112进行基于观测观测点Q1～QN的N个传感器22的观测信息来获取包括车辆6的多个车轴通过了观测点Q1～QN中的任意一个的时刻以及作为对于多个车轴分别对观测点Q1～QN中的任意一个的作用的响应即物理量的加速度强度的第二观测点信息的处理。即，第二观测点信息获取部112进行图20的第二观测点信息获取步骤的处理。第二观测点信息获取部112所获取的第二观测点信息被存储到存储部130中。
挠曲波形计算部113进行基于第一观测点信息获取部111所获取的第一观测点信息及第二观测点信息获取部112所获取的第二观测点信息、规定的系数p、以及基于上部结构7的结构模型的上部结构7的挠曲的近似式来计算多个车轴各自所引起的上部结构7的挠曲波形Hjk(t)的处理。即，挠曲波形计算部113进行图20的挠曲波形计算步骤的处理。挠曲波形计算部113所算出的挠曲波形Hjk(t)被存储到存储部130中。另外，规定的系数p和上部结构7的挠曲的近似式被预先存储到存储部130中。
移动体挠曲波形计算部114进行将挠曲波形计算部113所算出的车辆6的多个车轴各自所引起的上部结构7的挠曲波形Hjk(t)相加从而计算车辆6所引起的上部结构7的挠曲波形CPjm(t)的处理。即，移动体挠曲波形计算部114进行图20的移动体挠曲波形计算步骤的处理。移动体挠曲波形计算部114所算出的挠曲波形CPjm(t)被存储到存储部130中。
位移计算部115进行基于移动体挠曲波形计算部114所算出的车辆6所引起的上部结构7的挠曲波形CPjm(t)来计算上部结构7的位移CUest(t)的处理。即，位移计算部115进行图20的位移计算步骤的处理。位移计算部115所算出的位移CUest(t)被存储到存储部130中。
系数值计算部116进行以下处理：当多个车辆分别在上部结构7上单独行驶时，求出使未图示的位移仪计测到的上部结构7的位移的最大值CUmax与移动体挠曲波形计算部114算出的挠曲波形CPjm(t)的最大值CPjm-max的关联近似的近似直线，计算前述的式(47)的1次系数scu的值和0次系数icu的值。系数值计算部116算出的1次系数scu的值和0次系数icu的值被存储到存储部130中。
输出处理部117进行将位移计算部115算出的上部结构7的位移CUest(t)经由第二通信部140输出到服务器2的处理。即，输出处理部117进行图20的输出步骤的处理。
例如，控制部110基于来自操作部150的操作数据，切换计算未知的车辆6在上部结构7上行驶的时刻、上部结构7的位移等的第一模式及计算1次系数scu的值和0次系数icu的值的第二模式。例如，在N个传感器21和N个传感器22设置于上部结构7后，在控制部110被设定为第二模式的状态下进行基于多个车辆的载荷试验，在载荷试验结束后，控制部110被设定为第一模式。
在本实施方式中，控制部110是执行存储部130中存储的各种程序的处理器，执行在存储部130中存储的计测程序131，从而实现第一观测点信息获取部111、第二观测点信息获取部112、挠曲波形计算部113、移动体挠曲波形计算部114、位移计算部115、系数值计算部116、输出处理部117的各功能。换言之，计测程序131是使作为计算机的计测装置1执行图20所示的流程图的各顺序的程序。
处理器例如可以通过独立的硬件来实现各部的功能，或者也可以通过一体的硬件来实现各部的功能。例如，处理器可以包括硬件，该硬件可以包括处理数字信号的电路和处理模拟信号的电路中的至少一方。处理器可以是CPU(Central Processing Unit：中央处理单元)、GPU(Graphics Processing Unit：图形处理单元)或者DSP(Digital SignalProcessor：数字信号处理器)等。不过，控制部110可以构成为ASIC(Application SpecificIntegrated Circuit：专用集成电路)等特制IC(Integrated Circuit：集成电路)来实现各部的功能，也可以通过CPU和ASIC来实现各部的功能。
此外，控制部110也可以不包括系数值计算部116。例如，也可以是，服务器2或者其他装置进行计算1次系数scu的值和0次系数icu的值的处理，将这些值存储到计测装置1的存储部130中。
1-7.作用效果
在以上说明的第一实施方式的计测方法中，计测装置1基于观测观测点P1～PN的N个传感器21的观测信息，获取包括车辆6的各车轴通过了观测点P1～PN中的任意一个的时刻和加速度强度的第一观测点信息。另外，计测装置1基于观测观测点Q1～QN的N个传感器22的观测信息，获取包括车辆6的各车轴通过了观测点Q1～QN中的任意一个的时刻和加速度强度的第二观测点信息。然后，计测装置1基于第一观测点信息及第二观测点信息、规定的系数p、基于上部结构7的结构模型的上部结构7的挠曲的近似式(39)，通过式(44)计算各车轴所引起的上部结构7的挠曲波形Hjk(t)，将挠曲波形Hjk(t)相加来计算车辆6所引起的上部结构7的挠曲波形CPjm(t)。因而，根据第一实施方式的计测方法，计测装置1能够不用计测作为结构物的上部结构7的位移来算出作为在上部结构7上移动的移动体的车辆6所引起的上部结构7的挠曲波形。
另外，在第一实施方式的计测方法中，计测装置1能够通过相关式(47)，根据车辆6所引起的上部结构7的挠曲波形CPjm(t)计算车道Lj的位移CUest(t)。因而，根据第一实施方式的计测方法，计测装置1能够不用计测上部结构7的位移来推定上部结构7的位移。
另外，根据第一实施方式的计测方法，无需设定用于计测上部结构7的位移的观测点，因此减少观测装置的观测点的数量，计测系统变得简单，因此，能够减少计测所需的成本。
另外，根据第一实施方式的计测方法，观测点P1～PN、Q1～QN被设定于上部结构7的两端部，在上部结构7的中央部等不设定观测点，因此，计测系统10的构建、维护变得容易，能够减少计测所需的成本。
另外，根据第一实施方式的计测方法，计测装置1能够算出通过上部结构7的车辆6的车轴重量所引起的上部结构7的变形即挠曲波形，因此能够为用于预测上部结构7的损伤的桥梁5的维持管理提供足够的信息。
2.第二实施方式
在第一实施方式的计测方法中，计测装置1基于车辆6所引起的上部结构7的挠曲波形CPjm(t)计算上部结构7的位移CUest(t)，而在第二实施方式的计测方法中，基于车辆6所引起的上部结构7的挠曲波形CPjm(t)，计算基于车辆6的载荷。以下，对第二实施方式的与第一实施方式同样的构成要素标注相同的附图标记而省略或者简化与第一实施方式重复的说明，主要对与第一实施方式不同的内容进行说明。
在本实施方式中，将前述的式(45)所示的挠曲波形CPjm(t)与基于车辆的虚拟载荷CW(t)的关联用多项式进行近似。例如，如式(49)所示，将载荷CW(t)用挠曲波形CPjm(t)的1次式进行近似。在式(49)中，s是1次系数，i是0次系数。
[式49]
Figure BDA0002979166730000331
1次系数s和0次系数i例如是通过基于多个车辆的载荷试验来算出的。例如，在载荷试验前分别计测出多个车辆的重量。使多个车辆分别单独地在车道Lj上行驶，计测装置1生成车轴信息，并且获取各车辆的重量。然后，计测装置1以各车辆的重量为载荷CWtest并以根据车轴信息且通过式(45)得到的挠曲波形CPjm(t)的最大值为CPjm-max而绘制坐标图，求出近似直线的1次系数scp及0次系数icp。
图23是描绘出基于5辆车辆的载荷试验的结果的图。在图23中，横轴是车辆所引起的上部结构7的挠曲波形的最大值CPjm-max，纵轴是车辆的重量CWtest。在图23中，5个点按直线状配置，用虚线示出相对于5个点的近似直线。在图23的例子中，该近似直线的1次系数scp是-13714.92504，0次系数icp是-0.880179651。
计测装置1使用1次系数scp和0次系数icp、以及根据未知的车辆6的车轴信息并通过式(45)得到的挠曲波形CPjm(t)，通过式(50)计算基于车辆6的载荷CWest(t)。
[式50]
CWest(t)＝scpCPjm(t)+icp…(50)
在图24中，示出根据基于车辆的载荷CW(t)和挠曲波形CPjm(t)算出的虚拟载荷CWest(t)的一例。在图24中，实线表示基于车辆的载荷CW(t)，虚线表示根据挠曲波形CPjm(t)算出的载荷CWest(t)。如图24所示，根据挠曲波形CPjm(t)算出的载荷CWest(t)与基于车辆的载荷CW(t)近似。因而，计测装置1能够不用计测车道Lj的中央的位移而通过式(50)算出基于在车道Lj上行驶的未知的车辆6的载荷。
在式(50)中，0次系数icp是较小的值，将icp＝0代入式(50)，通过前述的式(44)和式(45)得到式(51)。
[式51]
Figure BDA0002979166730000341
根据式(51)，1次系数scp能够与规定的系数p进行替换，因此，规定的系数p是具有与1次系数scp同样的功能的系数。即，规定的系数p是使观测点Pj与观测点Qj之间的上部结构7的部位的挠曲与车辆6对该部位的载荷的关联近似的函数的系数。
图25是示出第二实施方式的计测方法的顺序的一例的流程图。在本实施方式中，计测装置1执行图25所示的顺序。在图25中，对进行与图20同样的处理的步骤标注相同的附图标记，省略其说明。
如图25所示，与第一实施方式同样地，计测装置1进行步骤S1～S4的处理。
接着，计测装置1通过式(50)，基于在步骤S4中算出的车辆6所引起的上部结构7的挠曲波形CPjm(t)，计算基于车辆6的虚拟载荷CWest(t)(步骤S8)。该步骤S8是载荷计算步骤。
接着，计测装置1将在步骤S8中算出的车辆6的载荷CWest(t)输出到服务器2(步骤S9)。该步骤S9是输出步骤。
计测装置1反复进行步骤S1～S9的处理，直至结束计测为止(步骤S10的N)。
图26是示出第二实施方式的计测装置1的构成例的图。在图26中，针对与图22同样的构成要素标注相同的附图标记。如图26所示，与第一实施方式同样地，计测装置1具有控制部110、第一通信部120、存储部130、第二通信部140以及操作部150。
第一通信部120、存储部130、第二通信部140以及操作部150分别进行的处理与第一实施方式相同，因此，省略其说明。
控制部110基于从设置于上部结构7的各传感器21、22输出的加速度数据，计算车辆6在上部结构7上行驶的时刻、基于车辆6的载荷等。
控制部110具备第一观测点信息获取部111、第二观测点信息获取部112、挠曲波形计算部113、移动体挠曲波形计算部114、载荷计算部118、系数值计算部116以及输出处理部117。
第一观测点信息获取部111、第二观测点信息获取部112、挠曲波形计算部113以及移动体挠曲波形计算部114分别进行的处理与第一实施方式相同，因此，省略其说明。
载荷计算部118进行基于移动体挠曲波形计算部114算出的车辆6所引起的上部结构7的挠曲波形CPjm(t)来计算基于车辆6的载荷CWest(t)的处理。即，载荷计算部118进行图25的载荷计算步骤的处理。载荷计算部118算出的载荷CWest(t)被存储到存储部130中。
系数值计算部116进行以下处理：当多个车辆分别在上部结构7上单独行驶时，求出使车辆的重量CWtest与移动体挠曲波形计算部114算出的挠曲波形CPjm(t)的最大值CPjm-max的关联近似的近似直线，计算式(50)的1次系数scp的值和0次系数icp的值。系数值计算部116算出的1次系数scp的值和0次系数icp的值被存储到存储部130中。
输出处理部117进行将载荷计算部118算出的车辆6的载荷CWest(t)经由第二通信部140输出到服务器2的处理。即，输出处理部117进行图25的输出步骤的处理。
例如，控制部110基于来自操作部150的操作数据，切换计算未知的车辆6在上部结构7上行驶的时刻、基于车辆6的载荷等的第一模式及计算1次系数scp的值和0次系数icp的值的第二模式。例如，在N个传感器21和N个传感器22设置于上部结构7后，在控制部110被设定为第二模式的状态下进行基于多个车辆的载荷试验，在载荷试验结束后，控制部110被设定为第一模式。
与第一实施方式同样地，控制部110是执行存储部130中存储的各种程序的处理器，执行存储部130中存储的计测程序131，从而实现第一观测点信息获取部111、第二观测点信息获取部112、挠曲波形计算部113、移动体挠曲波形计算部114、载荷计算部118、系数值计算部116、输出处理部117的各功能。换言之，计测程序131是使作为计算机的计测装置1执行图25所示的流程图的各顺序的程序。不过，控制部110也可以构成为ASIC等特制IC来实现各部的功能，也可以通过CPU和ASIC来实现各部的功能。
此外，控制部110也可以不包括系数值计算部116。例如，也可以是，服务器2或者其他装置进行计算1次系数scp的值和0次系数icp的值的处理，将这些值存储到计测装置1的存储部130中。
根据以上说明的第二实施方式的计测方法，与第一实施方式的计测方法同样地，计测装置1能够不用计测作为结构物的上部结构7的位移而算出作为在上部结构7上移动的移动体的车辆6所引起的上部结构7的挠曲波形。
另外，在第二实施方式的计测方法中，计测装置1能够通过相关式(50)并根据车辆6所引起的上部结构7的挠曲波形CPjm(t)计算基于车辆6的载荷CWest(t)。因而，根据第二实施方式的计测方法，计测装置1能够不用计测上部结构7的位移来推定由车辆6施加于上部结构7的载荷。
另外，根据第二实施方式的计测方法，能够与第一实施方式的计测方法同样地减少计测所需的成本。
3.第三实施方式
在第一实施方式和第二实施方式的计测方法中，将上部结构7的挠曲的近似式设为基于上部结构7的结构模型的式子，因此如前述的式(29)所示，归一化挠曲量wstd在载荷位置a比l/2小的区间和比l/2大的区间式子是不同的。与此相对，在第三实施方式的计测方法中，将上部结构7的挠曲的近似式用正弦波的半波长的式子进行近似，使得在载荷位置a比l/2小的区间和比l/2大的区间式子相同。以下，对第三实施方式的与第一实施方式或者第二实施方式同样的构成要素标注相同的附图标记，省略或者简化与第一实施方式或者第二实施方式重复的说明，主要对与第一实施方式和第二实施方式不同的内容进行说明。
在本实施方式中，用式(52)表示归一化挠曲量wstd。
[式52]
Figure BDA0002979166730000381
在式(52)中，载荷位置a是0≤a≤l，因此通过式(52)，归一化挠曲量wstd用正弦半波进行近似。
在图27中，用实线示出在观测位置x＝l/2的情况下通过式(52)算出的归一化挠曲量wstd的波形的一例。在图27中，横轴是载荷P的位置，纵轴是归一化挠曲量wstd。在图27的例子中，简支梁的支点间的距离l＝10。在图27中，用虚线示出通过前述的式(29)算出的归一化挠曲量wstd的波形。
如图27所示，通过式(52)算出的归一化挠曲量wstd的波形与通过式(29)算出的归一化挠曲量wstd的波形近似，能够将式(29)置换为式(52)。通过该置换，前述的式(39)被置换为式(53)。式(53)是以使最大值成为1的方式被归一化的式子。另外，式(53)是作为结构物的上部结构7的挠曲的近似式，是正弦波的半波长的波形的式子。
[式53]
Figure BDA0002979166730000382
根据式(53)，前述的式(44)置换为式(54)。
[式54]
Figure BDA0002979166730000383
挠曲波形计算部113通过式(54)计算在车道Lj上行驶的车辆6的各车轴所引起的上部结构7的挠曲波形Hjk(t)。
移动体挠曲波形计算部114通过前述的式(45)，将挠曲波形计算部113算出的挠曲波形Hjk(t)相加，计算车辆编号为m的车辆6所引起的上部结构7的挠曲波形CPjm(t)。
位移计算部115使用1次系数scu和0次系数icu、以及移动体挠曲波形计算部114算出的挠曲波形CPjm(t)，通过式(47)计算车道Lj的中央的位移CUest(t)。
或者，载荷计算部118使用1次系数scp和0次系数icp、以及移动体挠曲波形计算部114算出的挠曲波形CPjm(t)，通过式(50)计算基于车辆6的载荷CWest(t)。
在以上说明的第三实施方式中，计测装置1基于第一观测点信息及第二观测点信息、规定的系数p、以及正弦波的半波长的波形的式子即上部结构7的挠曲的近似式(53)，通过式(54)计算各车轴所引起的上部结构7的挠曲波形Hjk(t)，将挠曲波形Hjk(t)相加并计算车辆6所引起的上部结构7的挠曲波形CPjm(t)。因而，根据第三实施方式的计测方法，计测装置1能够不用计测作为结构物的上部结构7的位移而算出作为在上部结构7上移动的移动体的车辆6所引起的上部结构7的挠曲波形。
另外，在第三实施方式的计测方法中，与第一实施方式的计测方法或者第二实施方式的计测方法同样地，计测装置1能够不用计测上部结构7的位移而推定上部结构7的位移或者由车辆6施加于上部结构7的载荷。
另外，根据第三实施方式的计测方法，与第一实施方式的计测方法及第二实施方式的计测方法同样地能够减少计测所需的成本。
4.变形例
本发明不限于本实施方式，能够在本发明的主旨的范围内实施各种变形。
在上述的各实施方式中，计测装置1算出了上部结构7的位移和载荷中的任意一方，但也可以算出上部结构7的位移和载荷双方。例如，计测装置1可以通过前述的式(47)算出位移CUest(t)，通过前述的式(50)算出载荷CWest(t)。
另外，也可以是，计测装置1通过式(47)算出位移CUest(t)，通过规定的相关式将位移CUest(t)转换为载荷。例如，上部结构7的观测位置Rk的载荷CWk(t)与位移xk(t)的关系可以如式(55)那样示出。其中，载荷CWk(t)在BWIM中被设为位移波形所对应的载荷波形。式(55)的1次系数Sckk和0次系数Ick可以通过基于多个车辆的载荷试验而得到。
[式55]
CWk(t)＝Sckk·xk(t)+Ick…(55)
在式(55)中，若Ick被设为足够小的值，则可以得到式(56)。
[式56]
CWk(t)＝Sckk·xk(t)…(56)
在式(56)中，将位移xk(t)置换为位移CUest(t)，载荷CWk(t)与位移CUest(t)的相关式如式(57)所示。计测装置1能够通过相关式(57)将位移CUest(t)转换为载荷CWk(t)。
[式57]
CWk(t)＝Sckk·CUest(t)…(57)
另外，计测装置1也可以通过式(50)算出载荷CWest(t)，通过规定的相关式将载荷CWest(t)转换为位移。在式(56)中，当将载荷CWk(t)置换为载荷CWest(t)并整理时，位移xk(t)与载荷CWest(t)的相关式如式(58)所示。计测装置1能够通过相关式(58)将载荷CWk(t)转换为位移xk(t)。
[式58]
Figure BDA0002979166730000401
另外，在上述的各实施方式中，观测观测点P1～PN的观测装置和观测观测点Q1～QN的观测装置分别是加速度传感器，但是并不局限于此，例如，也可以是冲击传感器、麦克风、应变仪或者测压元件。观测装置与观测点无需一对一地对应，也可以是一个观测装置观测观测点P1～PN、Q1～QN的一部分或者全部。
冲击传感器检测冲击加速度作为对于车辆6的各车轴对观测点P1～PN、Q1～QN的作用的响应。计测装置1基于针对观测点P1～PN的冲击加速度获取第一观测点信息，基于针对观测点Q1～QN的冲击加速度获取第二观测点信息。麦克风检测声响作为对于车辆6的各车轴对观测点P1～PN、Q1～QN的作用的响应。计测装置1基于针对观测点P1～PN的声响获取第一观测点信息，基于针对观测点Q1～QN的声响获取第二观测点信息。应变仪、测压元件检测应力变化作为对于车辆6的各车轴对观测点P1～PN、Q1～QN的作用的响应。计测装置1基于观测点P1～PN的应力变化获取第一观测点信息，基于观测点Q1～QN的应力变化获取第二观测点信息。
另外，在上述的各实施方式中，车辆6在车道L1～LN上行驶的方向全部相同，但也可以是车道L1～LN中的至少一个车道与其他车道的车辆6的行驶方向不同。例如，车辆6可以在车道L1上向从观测点P1去往观测点Q1的方向行驶，在车道L2上向从观测点Q2去往观测点P2的方向行驶。在这种情况下，计测装置1基于从观测观测点P1的传感器21输出的加速度数据，获取车辆6向车道L1进入的时刻，基于从观测观测点Q1的传感器22输出的加速度数据，获取车辆6从车道L1退出的时刻。另外，计测装置1基于从观测观测点Q2的传感器22输出的加速度数据，获取车辆6向车道L2进入的时刻，基于从观测观测点P2的传感器21输出的加速度数据，获取车辆6从车道L2退出的时刻。
另外，在上述的各实施方式中，各传感器21、22分别设置于上部结构7的主梁G，但也可以设置于上部结构7的表面或内部、地板F的下表面、桥墩8a等。另外，在上述的各实施方式中，作为桥梁5，列举了公路桥的例子，但是并不局限于此，例如，桥梁5也可以是铁路桥。另外，在上述的各实施方式中，作为结构物，列举了桥梁的上部结构的例子，但是并不局限于此，结构物只要是通过移动体的移动而变形的结构物即可。
上述的实施方式和变形例是一例，并不限于这些内容。例如，也能够将各实施方式及各变形例适当组合。
本发明包括与在实施方式中说明的构成实质上相同的构成、例如功能、方法以及结果相同的构成、或者目的及效果相同的构成。另外，本发明包括将在实施方式中说明的构成的非本质性部分置换后的构成。另外，本发明包括与在实施方式中说明的构成能够起到同一作用效果的构成或者能够实现同一目的的构成。另外，本发明包括对在实施方式中说明的构成附加有公知技术的构成。

Claims (15)

1.一种计测方法，其特征在于，包括：
第一观测点信息获取步骤，基于对沿着移动体在结构物上移动的第一方向排列的所述结构物的第一观测点和第二观测点中的所述第一观测点进行观测的第一观测装置的观测信息，获取包括所述移动体的多个部位分别通过了所述第一观测点的时刻以及对于所述多个部位分别对所述第一观测点的作用的响应即物理量的第一观测点信息；
第二观测点信息获取步骤，基于观测所述第二观测点的第二观测装置的观测信息，获取包括所述多个部位分别通过了所述第二观测点的时刻以及对于所述多个部位分别对所述第二观测点的作用的响应即物理量的第二观测点信息；
挠曲波形计算步骤，基于所述第一观测点信息及所述第二观测点信息、规定的系数以及所述结构物的挠曲的近似式，计算所述多个部位各自所引起的所述结构物的挠曲波形；以及
移动体挠曲波形计算步骤，将在所述挠曲波形计算步骤中算出的所述多个部位各自所引起的所述结构物的挠曲波形相加，计算所述移动体所引起的所述结构物的挠曲波形。
2.根据权利要求1所述的计测方法，其特征在于，
所述第一观测点被设定在所述结构物的第一端部，
所述第二观测点被设定在所述结构物的与所述第一端部不同的第二端部。
3.根据权利要求1或2所述的计测方法，其特征在于，
所述移动体是铁路车辆、汽车、路面电车、建筑车辆或者军用车辆，
所述多个部位分别是车轴或者车轮。
4.根据权利要求1所述的计测方法，其特征在于，
所述结构物的挠曲的近似式是以所述第一观测点与所述第二观测点的中央位置处的挠曲的最大振幅进行归一化后的式子。
5.根据权利要求1所述的计测方法，其特征在于，
所述结构物的挠曲的近似式是基于所述结构物的结构模型的式子。
6.根据权利要求5所述的计测方法，其特征在于，
所述结构模型是将两端进行支承而成的简支梁。
7.根据权利要求1所述的计测方法，其特征在于，
所述结构物的挠曲的近似式是正弦波的半波长的波形的式子。
8.根据权利要求1所述的计测方法，其特征在于，
所述结构物是桥梁的上部结构，
所述上部结构是架设于相邻的桥台和桥墩、相邻的两个桥台或者相邻的两个桥墩中的任意一方的结构，
所述上部结构的两端部位于所述相邻的桥台和桥墩的位置、所述相邻的两个桥台的位置或者所述相邻的两个桥墩的位置，
所述桥梁是公路桥或者铁路桥。
9.根据权利要求1所述的计测方法，其特征在于，
所述规定的系数是使所述第一观测点与所述第二观测点之间的所述结构物的部位的挠曲与所述结构物的所述部位的位移或者对该部位的载荷的关联近似的函数的系数。
10.根据权利要求1所述的计测方法，其特征在于，
观测所述第一观测点的所述第一观测装置和观测所述第二观测点的所述第二观测装置是加速度传感器。
11.根据权利要求1所述的计测方法，其特征在于，
观测所述第一观测点的所述第一观测装置和观测所述第二观测点的所述第二观测装置是冲击传感器、麦克风、应变仪或者测压元件。
12.根据权利要求1所述的计测方法，其特征在于，
所述结构物为桥梁动态称重系统发挥功能的结构。
13.一种计测装置，其特征在于，包括：
第一观测点信息获取部，基于对沿着移动体在结构物上移动的第一方向排列的所述结构物的第一观测点和第二观测点中的所述第一观测点进行观测的第一观测装置的观测信息，获取包括所述移动体的多个部位分别通过了所述第一观测点的时刻以及对于所述多个部位分别对所述第一观测点的作用的响应即物理量的第一观测点信息；
第二观测点信息获取部，基于观测所述第二观测点的第二观测装置的观测信息，获取包括所述多个部位分别通过了所述第二观测点的时刻以及对于所述多个部位分别对所述第二观测点的作用的响应即物理量的第二观测点信息；
挠曲波形计算部，基于所述第一观测点信息及所述第二观测点信息、规定的系数以及所述结构物的挠曲的近似式，计算所述多个部位各自所引起的所述结构物的挠曲波形；以及
移动体挠曲波形计算部，将所述挠曲波形计算部算出的所述多个部位各自所引起的所述结构物的挠曲波形相加，计算所述移动体所引起的所述结构物的挠曲波形。
14.一种计测系统，其特征在于，具备：
权利要求13所述的计测装置；
所述第一观测装置；以及
所述第二观测装置。
15.一种存储介质，其特征在于，存储使计算机执行以下步骤的计测程序：
第一观测点信息获取步骤，基于对沿着移动体在结构物上移动的第一方向排列的所述结构物的第一观测点和第二观测点中的所述第一观测点进行观测的第一观测装置的观测信息，获取包括所述移动体的多个部位分别通过了所述第一观测点的时刻以及对于所述多个部位分别对所述第一观测点的作用的响应即物理量的第一观测点信息；
第二观测点信息获取步骤，基于观测所述第二观测点的第二观测装置的观测信息，获取包括所述多个部位分别通过了所述第二观测点的时刻以及对于所述多个部位分别对所述第二观测点的作用的响应即物理量的第二观测点信息；
挠曲波形计算步骤，基于所述第一观测点信息及所述第二观测点信息、规定的系数以及所述结构物的挠曲的近似式，计算所述多个部位各自所引起的所述结构物的挠曲波形；以及
移动体挠曲波形计算步骤，将在所述挠曲波形计算步骤中算出的所述多个部位各自所引起的所述结构物的挠曲波形相加，计算所述移动体所引起的所述结构物的挠曲波形。
CN202110282521.0A 2020-03-18 2021-03-16 计测方法、计测装置、计测系统以及存储介质 Active CN113494951B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2020-047306 2020-03-18
JP2020047306A JP7400566B2 (ja) 2020-03-18 2020-03-18 計測方法、計測装置、計測システム及び計測プログラム

Publications (2)

Publication Number Publication Date
CN113494951A true CN113494951A (zh) 2021-10-12
CN113494951B CN113494951B (zh) 2023-07-21

Family

ID=77747829

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110282521.0A Active CN113494951B (zh) 2020-03-18 2021-03-16 计测方法、计测装置、计测系统以及存储介质

Country Status (3)

Country Link
US (1) US11761811B2 (zh)
JP (1) JP7400566B2 (zh)
CN (1) CN113494951B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7375637B2 (ja) 2020-03-18 2023-11-08 セイコーエプソン株式会社 計測方法、計測装置、計測システム及び計測プログラム
JP7396139B2 (ja) 2020-03-18 2023-12-12 セイコーエプソン株式会社 計測方法、計測装置、計測システム及び計測プログラム
JP2021147819A (ja) 2020-03-18 2021-09-27 セイコーエプソン株式会社 計測方法、計測装置、計測システム及び計測プログラム
JP7447586B2 (ja) * 2020-03-18 2024-03-12 セイコーエプソン株式会社 計測方法、計測装置、計測システム及び計測プログラム
JP2021148537A (ja) 2020-03-18 2021-09-27 セイコーエプソン株式会社 計測方法、計測装置、計測システム及び計測プログラム
JP7400566B2 (ja) 2020-03-18 2023-12-19 セイコーエプソン株式会社 計測方法、計測装置、計測システム及び計測プログラム

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102998081A (zh) * 2012-12-17 2013-03-27 黑龙江省博凯科技开发有限公司 运用多套捷联惯性系统进行桥梁监测的方法
US20170038307A1 (en) * 2014-04-16 2017-02-09 Nec Corporation Information processing device, information processing method and medium
CN107588915A (zh) * 2017-10-18 2018-01-16 厦门大学 一种桥梁影响线识别方法及系统
CN107883915A (zh) * 2016-09-30 2018-04-06 上海数久信息科技有限公司 一种桥梁动挠度检测方法及装置
CN108829998A (zh) * 2018-06-25 2018-11-16 石家庄铁道大学 一种桥梁挠度的检测方法、检测装置及终端设备
CN109992827A (zh) * 2019-02-20 2019-07-09 深圳高速工程顾问有限公司 桥梁结构预警方法、装置、计算机设备和存储介质

Family Cites Families (40)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4958306A (en) * 1988-01-06 1990-09-18 Pacific Northwest Research & Development, Inc. Pavement inspection apparatus
US20060137914A1 (en) 2004-12-23 2006-06-29 Osmos S.A. Method and device for measuring the weight applied to the ground by at least one axle
US7668692B2 (en) 2005-10-11 2010-02-23 Tatom Frank B Method for weighing vehicles crossing a bridge
JP5164100B2 (ja) 2008-03-26 2013-03-13 株式会社エヌ・ティ・ティ・データ 橋梁通過車両監視システム、橋梁通過車両監視方法、およびコンピュータプログラム
KR20090126362A (ko) * 2008-06-04 2009-12-09 한국과학기술원 균열 측정 방법 및 측정 장치
US8990032B2 (en) 2010-12-30 2015-03-24 Sensys Networks, Inc. In-pavement wireless vibration sensor nodes, networks and systems
KR101163764B1 (ko) * 2012-03-29 2012-07-16 디비정보통신 주식회사 고정식 차량 부정 계중 방지장치
AT513258B1 (de) 2012-12-13 2014-03-15 Univ Wien Verfahren zum Vermessen eines fahrenden Fahrzeugs
US10459615B2 (en) 2014-12-11 2019-10-29 Rdi Technologies, Inc. Apparatus and method for analyzing periodic motions in machinery
JP6531984B2 (ja) 2015-06-12 2019-06-19 学校法人五島育英会 加速度記録を用いた変位応答算出法
JP2017020796A (ja) 2015-07-07 2017-01-26 公益財団法人鉄道総合技術研究所 橋梁動的応答解析方法
JP6555662B2 (ja) 2015-09-15 2019-08-07 学校法人五島育英会 計測装置、計測方法、プログラム、および計測システム
US20170097278A1 (en) 2015-10-02 2017-04-06 Seiko Epson Corporation Measurement instrument, measurement method, measurement system, and program
US10198640B2 (en) * 2015-10-02 2019-02-05 Seiko Epson Corporation Measuring device, measuring system, measuring method, and program
JP6627501B2 (ja) 2015-12-28 2020-01-08 セイコーエプソン株式会社 計測装置、減衰特性算出方法、プログラム、および計測システム
JP6604200B2 (ja) 2015-12-28 2019-11-13 セイコーエプソン株式会社 加速度センサー、計測システム、および計測装置
JP6587194B2 (ja) * 2016-02-15 2019-10-09 パナソニックＩｐマネジメント株式会社 応力分布測定方法及び応力分布測定システム
US11327475B2 (en) * 2016-05-09 2022-05-10 Strong Force Iot Portfolio 2016, Llc Methods and systems for intelligent collection and analysis of vehicle data
NL2016794B1 (en) 2016-05-18 2017-11-30 Heijmans N V Method for determining the structural integrity of an infrastructural element
US11237546B2 (en) * 2016-06-15 2022-02-01 Strong Force loT Portfolio 2016, LLC Method and system of modifying a data collection trajectory for vehicles
CN105953996B (zh) 2016-06-30 2017-04-19 东南大学 一种基于冲击振动的桥梁检测评估方法与设备
JP6809691B2 (ja) 2016-07-26 2021-01-06 学校法人五島育英会 計測装置、計測方法、および計測システム
JP6543863B2 (ja) 2016-08-24 2019-07-17 公益財団法人鉄道総合技術研究所 鉄道橋の構造性能調査方法
JP6863383B2 (ja) * 2016-08-31 2021-04-21 日本電気株式会社 欠陥検出装置、欠陥検出方法、及びコンピュータ読み取り可能記録媒体
JP6775213B2 (ja) 2016-10-19 2020-10-28 学校法人五島育英会 計測装置、計測システム、プログラム及び計測方法
KR101693759B1 (ko) 2016-11-29 2017-01-09 한국건설기술연구원 로드셀이 구비된 신축이음장치를 이용한 교량 안전진단 장치 및 그 방법
JP6883768B2 (ja) 2017-02-28 2021-06-09 パナソニックＩｐマネジメント株式会社 モニタリングシステム
EP3929551A1 (en) 2017-04-26 2021-12-29 Panasonic Intellectual Property Management Co., Ltd. Stress measurement device, stress measurement system, and stress measurement method
JP6993143B2 (ja) 2017-09-07 2022-01-13 東海旅客鉄道株式会社 鉄道橋梁のたわみ測定装置
US20200363287A1 (en) 2017-11-16 2020-11-19 Nec Corporation Damage diagnosing device, damage diagnosing method, and recording medium having damage diagnosing program stored thereon
JP2020125913A (ja) * 2019-02-01 2020-08-20 セイコーエプソン株式会社 計測装置、計測システム及び計測方法
US11867592B2 (en) 2019-03-01 2024-01-09 Nec Corporation Structure displacement amount measurement apparatus
JPWO2020194539A1 (zh) * 2019-03-26 2020-10-01
JP6894656B1 (ja) 2019-07-26 2021-06-30 株式会社Ｔｔｅｓ 構造物の変位を推定するためのシステム、方法、及びプログラム
JP2021148537A (ja) 2020-03-18 2021-09-27 セイコーエプソン株式会社 計測方法、計測装置、計測システム及び計測プログラム
JP7375637B2 (ja) * 2020-03-18 2023-11-08 セイコーエプソン株式会社 計測方法、計測装置、計測システム及び計測プログラム
JP7400566B2 (ja) 2020-03-18 2023-12-19 セイコーエプソン株式会社 計測方法、計測装置、計測システム及び計測プログラム
JP2021147819A (ja) 2020-03-18 2021-09-27 セイコーエプソン株式会社 計測方法、計測装置、計測システム及び計測プログラム
JP7396139B2 (ja) 2020-03-18 2023-12-12 セイコーエプソン株式会社 計測方法、計測装置、計測システム及び計測プログラム
JP2022131027A (ja) * 2021-02-26 2022-09-07 セイコーエプソン株式会社 計測方法、計測装置、計測システム及び計測プログラム

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102998081A (zh) * 2012-12-17 2013-03-27 黑龙江省博凯科技开发有限公司 运用多套捷联惯性系统进行桥梁监测的方法
US20170038307A1 (en) * 2014-04-16 2017-02-09 Nec Corporation Information processing device, information processing method and medium
CN107883915A (zh) * 2016-09-30 2018-04-06 上海数久信息科技有限公司 一种桥梁动挠度检测方法及装置
CN107588915A (zh) * 2017-10-18 2018-01-16 厦门大学 一种桥梁影响线识别方法及系统
CN108829998A (zh) * 2018-06-25 2018-11-16 石家庄铁道大学 一种桥梁挠度的检测方法、检测装置及终端设备
CN109992827A (zh) * 2019-02-20 2019-07-09 深圳高速工程顾问有限公司 桥梁结构预警方法、装置、计算机设备和存储介质

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
袁向荣: "桥梁动挠度计算的加速度积分方法", 《广州大学学报(自然科学版)》 *

Also Published As

Publication number Publication date
CN113494951B (zh) 2023-07-21
JP2021147827A (ja) 2021-09-27
JP7400566B2 (ja) 2023-12-19
US11761811B2 (en) 2023-09-19
US20210293604A1 (en) 2021-09-23

Similar Documents

Publication Publication Date Title
CN113494951B (zh) 计测方法、计测装置、计测系统以及存储介质
CN113494950B (zh) 测量方法、测量装置、测量系统以及存储介质
US11714021B2 (en) Measurement method, measurement device, measurement system, and measurement program
US11881102B2 (en) Measurement method, measurement device, measurement system, and measurement program
US11713993B2 (en) Bridge displacement measurement method
US11408761B2 (en) Measurement method, measurement device, measurement system, and measurement program
JP2021148526A (ja) 計測方法、計測装置、計測システム及び計測プログラム
JP2021147820A (ja) 計測方法、計測装置、計測システム及び計測プログラム
JP2021148527A (ja) 計測方法、計測装置、計測システム及び計測プログラム
JP2021148535A (ja) 計測方法、計測装置、計測システム及び計測プログラム
JP2021148538A (ja) 計測方法、計測装置、計測システム及び計測プログラム
JP2021147818A (ja) 計測方法、計測装置、計測システム及び計測プログラム
US20230341289A1 (en) Measurement Method, Measurement Device, Measurement System, And Non-Transitory Computer-Readable Storage Medium Storing Measurement Program
CN115541150A (zh) 测量方法、测量装置、测量系统及存储介质
JP2021148525A (ja) 計測方法、計測装置、計測システム及び計測プログラム
US20230003575A1 (en) Measurement Method, Measurement Device, Measurement System, And Measurement Program
US20230341255A1 (en) Measurement Method, Measurement Device, Measurement System, And Non-Transitory Computer-Readable Storage Medium Storing Measurement Program
JP2021147826A (ja) 計測方法、計測装置、計測システム及び計測プログラム
CN115541152A (zh) 测量方法、测量装置、测量系统及存储介质
JP2021147825A (ja) 計測方法、計測装置、計測システム及び計測プログラム
JP2023084856A (ja) 計測方法、計測装置、計測システム及び計測プログラム
JP2023161194A (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
GR01 Patent grant
GR01 Patent grant