JP2021148537A - 計測方法、計測装置、計測システム及び計測プログラム - Google Patents

Abstract

【課題】移動体の静止時と移動時の荷重の差を考慮して構造物のたわみ波形を算出可能な計測方法を提供すること。【解決手段】第１観測点情報を用いて、移動体の先頭の部位が第１の観測点を通過した先頭時刻から各部位が前記第１の観測点を通過した時刻までの時間と、前記先頭時刻から各部位の前記第１の観測点への作用に対する応答である第１の物理量の総和を所定の配分比で配分する時刻までの時間とに基づいて、前記第１の物理量を補正する補正係数を算出するステップと、前記第１観測点情報及び第２観測点情報と、所定の係数と、前記補正係数と、構造物のたわみの近似式とに基づいて、前記複数の部位による前記構造物のたわみ波形を算出するステップと、前記構造物のたわみ波形を加算して、前記移動体による前記構造物のたわみ波形を算出するステップと、を含む、計測方法。【選択図】図２６

Description

本発明は、計測方法、計測装置、計測システム及び計測プログラムに関する。

特許文献１には、橋梁の維持管理をする上で、橋梁を通過する大型車両の車軸重量が、橋梁の損傷を予測するために重要な情報であって、この軸重測定ため、橋梁の主桁に設置したひずみ計から車両通過時のひずみ値を連続測定し、軸重を算出する手法Weight In Motionが提案されており、橋梁の主桁に配置されたひずみ計で計測したひずみ波形に基づいて、橋梁を通過する車両の車重を計測する橋梁通過車両監視システムが記載されている。詳細には、橋梁通過車両監視システムは、ひずみ計を配置して、ひずみ計で計測したひずみ波形から車軸の通過タイミングを検出して車両の軸間比率を算出し、算出した軸間比率と軸間距離データベースに登録された軸間距離から算出される軸間比率とを比較し、車両の軸間距離、車速および車種を特定する。また、橋梁通過車両監視システムは、車軸の通過タイミングに合わせて、基準軸重ひずみ波形を時間軸上に配置したひずみ波形を生成し、基準軸重ひずみ波形とひずみ計で計測したひずみ波形とを比較して各軸の軸重を算出する。そして、橋梁通過車両監視システムは、各軸の軸重を合計することにより車重を算出する。

特開２００９−２３７８０５号公報

車両の走行時における各車軸の荷重は、重心からの距離モーメントの影響を受けるため、車両の静止時の各車軸の荷重とは大きく異なる場合がある。特許文献１に記載のシステムでは、ひずみ波形と軸間距離データベースを用いることにより橋梁の変位を計測せずに車両の車重を計測することはできるが、車両等の移動体の静止時と移動時の荷重の差を考慮して、橋梁等の構造物を移動する移動体による構造物のたわみ波形を算出することはできない。

本発明に係る計測方法の一態様は、
移動体が構造物を移動する第１方向に沿って並ぶ前記構造物の第１の観測点及び第２の観測点のうちの前記第１の観測点を観測する観測装置による観測情報に基づいて、前記移動体の複数の部位が前記第１の観測点をそれぞれ通過した時刻及び前記複数の部位のそれぞれの前記第１の観測点への作用に対する応答である第１の物理量を含む第１観測点情報を取得する第１観測点情報取得ステップと、
前記第２の観測点を観測する観測装置による観測情報に基づいて、前記複数の部位が前記第２の観測点をそれぞれ通過した時刻及び前記複数の部位のそれぞれの前記第２の観測点への作用に対する応答である第２の物理量を含む第２観測点情報を取得する第２観測点情報取得ステップと、
前記第１観測点情報を用いて、前記複数の部位のうちの先頭の部位が前記第１の観測点を通過した時刻である第１の先頭時刻から前記複数の部位が前記第１の観測点をそれぞれ通過した時刻までの時間と、前記第１の先頭時刻から前記第１の物理量の総和を所定の配分比で配分する時刻である第１の基準時刻までの時間とを算出する第１時間算出ステップと、
前記第１の先頭時刻から前記複数の部位が前記第１の観測点をそれぞれ通過した時刻までの時間と、前記第１の先頭時刻から前記第１の基準時刻までの時間とに基づいて、前記第１の物理量を補正する補正係数を算出する補正係数算出ステップと、
前記第１観測点情報及び前記第２観測点情報と、所定の係数と、前記補正係数と、前記構造物のたわみの近似式とに基づいて、前記複数の部位のそれぞれによる前記構造物のたわみ波形を算出するたわみ波形算出ステップと、
前記たわみ波形算出ステップで算出した前記複数の部位のそれぞれによる前記構造物のたわみ波形を加算して、前記移動体による前記構造物のたわみ波形を算出する移動体たわみ波形算出ステップと、を含む。

前記計測方法の一態様は、
前記第２観測点情報を用いて、前記複数の部位のうちの前記先頭の部位が前記第２の観測点を通過した時刻である第２の先頭時刻から前記複数の部位が前記第２の観測点をそれぞれ通過した時刻までの時間と、前記第２の先頭時刻から前記第２の物理量の総和を前記所定の配分比で配分する時刻である第２の基準時刻までの時間とを算出する第２時間算出ステップをさらに含み、
前記補正係数算出ステップでは、
前記第１の先頭時刻から前記複数の部位が前記第１の観測点をそれぞれ通過した時刻までの時間と、前記第１の先頭時刻から前記第１の基準時刻までの時間と、前記第２の先頭時刻から前記複数の部位が前記第２の観測点をそれぞれ通過した時刻までの時間と、前記第２の先頭時刻から前記第２の基準時刻までの時間とに基づいて、前記第１の物理量及び前記第２の物理量を補正する前記補正係数を算出してもよい。

前記計測方法の一態様において、
前記配分比をＲ：１−Ｒとしたとき、Ｒは０．４以上０．６以下であってもよい。

前記計測方法の一態様において、
前記第１の観測点は前記構造物の第１の端部に設定され、
前記第２の観測点は前記構造物の前記第１の端部とは異なる第２の端部に設定されてもよい。

前記計測方法の一態様において、
前記構造物は、橋梁の上部構造であって、
前記上部構造は、隣り合う橋台と橋脚、隣り合う２つの橋台、又は、隣り合う２つの橋脚のいずれか１つに渡された構造であり、
前記上部構造の両端部は、前記隣り合う橋台と橋脚の位置、前記隣り合う２つの橋台の位置、又は、前記隣り合う２つの橋脚の位置にあり、
前記橋梁は、道路橋又は鉄道橋であってもよい。

前記計測方法の一態様において、
前記移動体は、鉄道車両、自動車、路面電車、建設車両、又は軍用車両であり、
前記複数の部位のそれぞれは車軸又は車輪であってもよい。

前記計測方法の一態様において、
前記構造物のたわみの近似式は、前記構造物の構造モデルに基づく式であってもよい。

前記計測方法の一態様において、
前記構造モデルは、両端を支持した単純梁であってもよい。

前記計測方法の一態様において、
前記構造物のたわみの近似式は、前記第１の観測点と前記第２の観測点との中央位置におけるたわみの最大振幅で規格化された式であってもよい。

前記計測方法の一態様において、
前記第１の観測点を観測する前記観測装置及び前記第２の観測点を観測する前記観測装置は、加速度センサーであってもよい。

前記計測方法の一態様において、
前記第１の観測点を観測する前記観測装置及び前記第２の観測点を観測する前記観測装置は、衝撃センサー、マイクロホン、歪計又はロードセルであってもよい。

前記計測方法の一態様において、
前記構造物は、ＢＷＩＭ（Bridge Weigh in Motion）が機能する構造であってもよい。

本発明に係る計測装置の一態様は、
移動体が構造物を移動する第１方向に沿って並ぶ前記構造物の第１の観測点及び第２の観測点のうちの前記第１の観測点を観測する観測装置による観測情報に基づいて、前記移動体の複数の部位が前記第１の観測点をそれぞれ通過した時刻及び前記複数の部位のそれぞれの前記第１の観測点への作用に対する応答である第１の物理量を含む第１観測点情報を取得する第１観測点情報取得部と、
前記第２の観測点を観測する観測装置による観測情報に基づいて、前記複数の部位が前記第２の観測点をそれぞれ通過した時刻及び前記複数の部位のそれぞれの前記第２の観測点への作用に対する応答である第２の物理量を含む第２観測点情報を取得する第２観測点情報取得部と、
前記第１観測点情報を用いて、前記複数の部位のうちの先頭の部位が前記第１の観測点を通過した時刻である第１の先頭時刻から前記複数の部位が前記第１の観測点をそれぞれ通過した時刻までの時間と、前記第１の先頭時刻から前記第１の物理量の総和を所定の配分比で配分する時刻である第１の基準時刻までの時間とを算出する第１時間算出部と、
前記第１の先頭時刻から前記複数の部位が前記第１の観測点をそれぞれ通過した時刻までの時間と、前記第１の先頭時刻から前記第１の基準時刻までの時間とに基づいて、前記第１の物理量を補正する補正係数を算出する補正係数算出部と、
前記第１観測点情報及び前記第２観測点情報と、所定の係数と、前記補正係数と、前記構造物のたわみの近似式とに基づいて、前記複数の部位のそれぞれによる前記構造物のたわみ波形を算出するたわみ波形算出部と、
前記たわみ波形算出部が算出した前記複数の部位のそれぞれによる前記構造物のたわみ波形を加算して、前記移動体による前記構造物のたわみ波形を算出する移動体たわみ波形算出部と、を含む。

本発明に係る計測システムの一態様は、
前記計測装置の一態様と、
前記第１の観測点を観測する前記観測装置と、
前記第２の観測点を観測する前記観測装置と、を備える。

本発明に係る計測プログラムの一態様は、
移動体が構造物を移動する第１方向に沿って並ぶ前記構造物の第１の観測点及び第２の観測点のうちの前記第１の観測点を観測する観測装置による観測情報に基づいて、前記移動体の複数の部位が前記第１の観測点をそれぞれ通過した時刻及び前記複数の部位のそれぞれの前記第１の観測点への作用に対する応答である第１の物理量を含む第１観測点情報を取得する第１観測点情報取得ステップと、
前記第２の観測点を観測する観測装置による観測情報に基づいて、前記複数の部位が前記第２の観測点をそれぞれ通過した時刻及び前記複数の部位のそれぞれの前記第２の観測点への作用に対する応答である第２の物理量を含む第２観測点情報を取得する第２観測点情報取得ステップと、
前記第１観測点情報を用いて、前記複数の部位のうちの先頭の部位が前記第１の観測点を通過した時刻である第１の先頭時刻から前記複数の部位が前記第１の観測点をそれぞれ通過した時刻までの時間と、前記第１の先頭時刻から前記第１の物理量の総和を所定の配分比で配分する時刻である第１の基準時刻までの時間とを算出する第１時間算出ステップと、
前記第１の先頭時刻から前記複数の部位が前記第１の観測点をそれぞれ通過した時刻までの時間と、前記第１の先頭時刻から前記第１の基準時刻までの時間とに基づいて、前記第１の物理量を補正する補正係数を算出する補正係数算出ステップと、
前記第１観測点情報及び前記第２観測点情報と、所定の係数と、前記補正係数と、前記構造物のたわみの近似式とに基づいて、前記複数の部位のそれぞれによる前記構造物のたわみ波形を算出するたわみ波形算出ステップと、
前記たわみ波形算出ステップで算出した前記複数の部位のそれぞれによる前記構造物のたわみ波形を加算して、前記移動体による前記構造物のたわみ波形を算出する移動体たわみ波形算出ステップと、をコンピューターに実行させる。

計測システムの構成例を示す図。 センサー及び観測点の配置例を示す図。 センサー及び観測点の配置例を示す図。 加速度センサーが検出する加速度の説明図。 車軸情報の一例を示す図。 センサー及び観測点の配置例を示す図。 センサー及び観測点の配置例を示す図。 観測点に対して検出される加速度の一例を示す図。 図８の各時刻の加速度振幅を加速度強度に変換した図。 図９の加速度強度を所定の閾値で２値化した図。 図１０に対して退出時刻のパターンをスライドさせた図。 橋梁の上部構造の構造モデルの説明図。 橋梁の上部構造の構造モデルの説明図。 規格化たわみ量の波形の一例を示す図。 規格化たわみ量モデルの一例を示す図。 補正されない場合の各車軸による橋梁のたわみ波形の一例を示す図。 ４軸車、３軸車、２軸車の重心位置及び車軸の位置の一例を示す図。 時間Δｔ_ｋの一例を示す図。 規格化時間と加算規格化衝撃パワーとの相関の一例を示す図。 規格化時間とシフト加算規格化衝撃パワーとの相関の一例を示す図。 基準時間と基準時刻との関係を示す図。 補正された場合の各車軸による橋梁のたわみ波形の一例を示す図。 車両による橋梁のたわみ波形の一例を示す図。 荷重試験の結果をプロットした図。 計測変位とたわみ波形から算出される変位の一例を示す図。 第１実施形態の計測方法の手順の一例を示すフローチャート図。 第１実施形態における第１時間算出ステップの手順の一例を示すフローチャート図。 第１実施形態における補正係数算出ステップの手順の一例を示すフローチャート図。 第１実施形態における計測装置の構成例を示す図。 第２実施形態の計測方法の手順の一例を示すフローチャート図。 第２実施形態における第１時間算出ステップの手順の一例を示すフローチャート図。 第２時間算出ステップの手順の一例を示すフローチャート図。 第２実施形態における補正係数算出ステップの手順の一例を示すフローチャート図。 第２実施形態における計測装置の構成例を示す図。

以下、本発明の好適な実施形態について図面を用いて詳細に説明する。なお、以下に説明する実施の形態は、特許請求の範囲に記載された本発明の内容を不当に限定するものではない。また以下で説明される構成の全てが本発明の必須構成要件であるとは限らない。

１．第１実施形態
１−１．計測システム
以下では、構造物が橋梁の上部構造であり、移動体が車両である場合を例に挙げ、本実施形態の計測方法を実現するための計測システムについて説明する。本実施形態に係る橋梁を通過する車両は、鉄道車両、自動車、路面電車、建設車両、又は軍用車両等の重量が大きく、ＢＷＩＭ（Bridge Weigh in Motion）で計測可能な車両である。ＢＷＩＭは、橋梁を「はかり」に見立て、橋梁の変形を計測することにより、橋梁を通行する車両の重量、軸数などを測定する技術である。変形やひずみなどの応答から通行する車両の重量を解析可能な橋梁の上部構造は、ＢＷＩＭが機能する構造物であり、橋梁の上部構造への作用と応答の間の物理的なプロセスを応用するＢＷＩＭシステムが通行する車両の重量の計測を可能にする。

図１は、本実施形態に係る計測システムの一例を示す図である。図１に示すように、本実施形態に係る計測システム１０は、計測装置１と、橋梁５の上部構造７に設けられる少なくとも１つのセンサー２１と、少なくとも１つのセンサー２２と、を有している。また、計測システム１０は、サーバー２を有してもよい。

橋梁５は上部構造７と下部構造８からなり、上部構造７は、床板Ｆ、主桁Ｇ、不図示の横桁等からなる橋床７ａと、支承７ｂと、を含む。下部構造８は、橋脚８ａと、橋台８ｂと、を含む。上部構造７は、隣り合う橋台８ｂと橋脚８ａ、隣り合う２つの橋台８ｂ、又は、隣り合う２つの橋脚８ａのいずれか１つに渡された構造である。上部構造７の両端部は、隣り合う橋台８ｂと橋脚８ａの位置、隣り合う２つの橋台８ｂの位置、又は、隣り合う２つの橋脚８ａの位置にある。

計測装置１と各センサー２１，２２とは、例えば、不図示のケーブルで接続され、ＣＡＮ（Controller Area Network）等の通信ネットワークを介して通信を行う。あるいは、計測装置１と各センサー２１，２２とは、無線ネットワークを介して通信を行ってもよい。

例えば、各センサー２１は、移動体である車両６の上部構造７への進入による衝撃を表すデータを出力し、各センサー２２は、車両６の上部構造７からの退出による衝撃を表すデータを出力する。本実施形態では、各センサー２１，２２は加速度センサーであり、例えば、水晶加速度センサーであってもよいし、ＭＥＭＳ（Micro Electro Mechanical Systems）加速度センサーであってもよい。

本実施形態では、各センサー２１は上部構造７の長手方向の第１の端部に設置され、各センサー２２は上部構造７の長手方向の第１の端部とは異なる第２の端部に設置されている。

各センサー２１は、車両６の上部構造７への進入の際に生じる上部構造７の加速度を検出し、各センサー２２は、車両６の上部構造７からの退出の際に生じる上部構造７の加速度を検出する。すなわち、本実施形態では、各センサー２１は、車両６の上部構造７への進入を検知する加速度センサーであり、各センサー２２は、車両６の上部構造７からの退出を検知する加速度センサーである。

計測装置１は、各センサー２１，２２から出力される加速度データに基づいて、車両６の走行による上部構造７の撓みの変位を算出する。

計測装置１とサーバー２とは、例えば、携帯電話の無線ネットワーク及びインターネット等の通信ネットワーク４を介して、通信を行うことができる。計測装置１は、車両６が上部構造７を走行した時刻や車両６の走行による上部構造７の変位等の情報をサーバー２に送信する。サーバー２は、当該情報を不図示の記憶装置に記憶し、例えば、当該情報に基づいて過積載の車両の監視や上部構造７の異常判定等の処理を行ってもよい。

なお、本実施形態では、橋梁５は、道路橋であり、例えば、鋼橋や桁橋、ＲＣ（Reinforced-Concrete）橋等である。

図２及び図３は、各センサー２１，２２の上部構造７への設置例を示す図である。なお、図２は、上部構造７をその上方から見た図であり、図３は、図２をＡ−Ａ線又はＢ−Ｂ線で切断した断面図である。

図２及び図３に示すように、上部構造７は、移動体である車両６が移動し得る第１〜第Ｎの経路としてのＮ個のレーンＬ_１〜Ｌ_Ｎ、及びＫ個の主桁Ｇ_１〜Ｇ_Ｋを有している。ここで、Ｎ，Ｋはそれぞれ１以上の整数である。なお、図２及び図３の例では、主桁Ｇ_１〜Ｇ_Ｋの各位置がレーンＬ_１〜Ｌ_Ｎの各境界の位置と一致しており、Ｎ＝Ｋ−１であるが、主桁Ｇ_１〜Ｇ_Ｋの各位置がレーンＬ_１〜Ｌ_Ｎの各境界の位置と一致している必要はなく、Ｎ≠Ｋ−１であってもよい。

図２及び図３の例では、上部構造７の長手方向の第１の端部ＥＡ１において、主桁Ｇ_１〜Ｇ_Ｋ−１のそれぞれにセンサー２１が設けられ、上部構造７の長手方向の第２の端部ＥＡ２において、主桁Ｇ_１〜Ｇ_Ｋ−１のそれぞれにセンサー２２が設けられている。図２及び図３の例では、Ｎ＝Ｋ−１であり、主桁Ｇ_Ｋにセンサー２１，２２が設けられていないが、主桁Ｇ_Ｋにセンサー２１，２２が設けられ、主桁Ｇ_１〜Ｇ_Ｋ−１のいずれか１つにセンサー２１，２２が設けられていなくてもよい。あるいは、Ｎ＝Ｋであり、主桁Ｇ_１〜Ｇ_Ｋのそれぞれにセンサー２１，２２が設けられていてもよい。

なお、各センサー２１，２２を上部構造７の床板Ｆに設けると、走行車両によって破壊するおそれがあり、また橋床７ａの局部的な変形により測定精度が影響を受けるおそれがあるため、図２及び図３の例では、各センサー２１，２２は上部構造７の主桁Ｇ_１〜Ｇ_Ｋ−１に設けられている。

本実施形態では、Ｎ個のセンサー２１に対応付けてＮ個の観測点Ｐ_１〜Ｐ_Ｎがそれぞれ設定されている。観測点Ｐ_１〜Ｐ_Ｎは、車両６が上部構造７を移動する第１方向と交差する第２方向に沿って並ぶ上部構造７のＮ個の観測点である。図２及び図３の例では、１以上Ｎ以下の各整数ｊに対して、観測点Ｐ_ｊは、第１の端部ＥＡ１において、主桁Ｇ_ｊに設けられたセンサー２１の鉛直上方向にある床板Ｆの表面の位置に設定されている。すなわち、主桁Ｇ_ｊに設けられたセンサー２１は、観測点Ｐ_ｊを観測する観測装置である。観測点Ｐ_ｊを観測するセンサー２１は、車両６の走行により観測点Ｐ_ｊに生じる加速度を検出可能な位置に設けられていればよいが、観測点Ｐ_ｊに近い位置に設けられることが望ましい。このように、観測点Ｐ_１〜Ｐ_ＮはＮ個のセンサー２１と１対１の関係にある。

また、本実施形態では、Ｎ個のセンサー２２に対応付けてＮ個の観測点Ｑ_１〜Ｑ_Ｎがそれぞれ設定されている。観測点Ｑ_１〜Ｑ_Ｎは、車両６が上部構造７を移動する第１方向と交差する第３方向に沿って並ぶ上部構造７のＮ個の観測点である。図２及び図３の例では、１以上Ｎ以下の各整数ｊに対して、観測点Ｑ_ｊは、第２の端部ＥＡ２において、主桁Ｇ_ｊに設けられたセンサー２２の鉛直上方向にある床板Ｆの表面の位置に設定されている。すなわち、主桁Ｇ_ｊに設けられたセンサー２２は、観測点Ｑ_ｊを観測する観測装置である。観測点Ｑ_ｊを観測するセンサー２２は、車両６の走行により観測点Ｑ_ｊに生じる加速度を検出可能な位置に設けられていればよいが、観測点Ｑ_ｊに近い位置に設けられることが望ましい。このように、観測点Ｑ_１〜Ｑ_ＮはＮ個のセンサー２２と１対１の関係にある。

本実施形態では、Ｎ個の観測点Ｐ_１〜Ｐ_Ｎは、それぞれレーンＬ_１〜Ｌ_Ｎに対応付けられている。同様に、Ｎ個の観測点Ｑ_１〜Ｑ_Ｎは、それぞれレーンＬ_１〜Ｌ_Ｎに対応付けられている。レーンＬ_ｊに対応付けて設定される観測点Ｐ_ｊ及び観測点Ｑ_ｊは、車両６が上部構造７を移動する第１方向に沿って並んでいる。図２及び図３の例では、第１方向は、上部構造７のレーンＬ_１〜Ｌ_Ｎに沿うＸ方向、すなわち、上部構造７の長手方向である。また、第２方向及び第３方向は、車両６が走行する上部構造７の面内においてＸ方向と直交するＹ方向、すなわち、上部構造７の幅員方向である。ただし、レーンＬ_１〜Ｌ_Ｎがそれぞれ曲線状である場合等は、第２方向と第３方向は一致しなくてもよい。また、第２方向及び第３方向は、第１方向と直交していなくてもよく、例えば、上部構造７の車両６が進入する側の端から観測点Ｐ_１〜Ｐ_Ｎまでの距離や、上部構造７の車両６が退出する側の端から観測点Ｑ_１〜Ｑ_Ｎまでの距離が異なっていてもよい。なお、１以上Ｎ以下の各整数ｊに対して、観測点Ｐ_ｊは「第１の観測点」の一例であり、観測点Ｑ_ｊは「第２の観測点」の一例である。

なお、センサー２１，２２の数及び設置位置は、図２及び図３に示した例には限定されず種々の変形実施が可能である。

計測装置１は、各センサー２１，２２から出力される加速度データに基づいて、第１方向であるＸ方向と、第２方向及び第３方向であるＹ方向とそれぞれ交差する第４方向の加速度を取得する。観測点Ｐ_１〜Ｐ_Ｎ，Ｑ_１〜Ｑ_Ｎは、衝撃によりＸ方向及びＹ方向と直交する方向に変位するので、計測装置１は、衝撃の大きさを正確に算出するために、Ｘ方向及びＹ方向と直交する第４方向、すなわち、床板Ｆの法線方向の加速度を取得するのが望ましい。

図４は、センサー２１，２２が検出する加速度を説明する図である。センサー２１，２２は、互いに直交する３軸の各軸方向に生じる加速度を検出する加速度センサーである。

車両６の上部構造７への進入により観測点Ｐ_１〜Ｐ_Ｎに加わる衝撃を検出するために、各センサー２１は、３つの検出軸であるｘ軸、ｙ軸、ｚ軸のうち、１軸が第１方向及び第２方向と交差する方向となるように設置される。同様に、車両６の上部構造７からの退出により観測点Ｑ_１〜Ｑ_Ｎに加わる衝撃を検出するために、各センサー２２は、３つの検出軸であるｘ軸、ｙ軸、ｚ軸のうち、１軸が第１方向及び第３方向と交差する方向となるように設置される。図２及び図３の例では、第１方向はＸ方向であり、第２方向及び第３方向はＹ方向であるから、各センサー２１，２２は、１軸がＸ方向及びＹ方向と交差する方向となるように設置される。観測点Ｐ_１〜Ｐ_Ｎ，Ｑ_１〜Ｑ_Ｎは、衝撃によりＸ方向及びＹ方向と直交する方向に変位するので、衝撃の大きさを正確に検出するために、理想的には、各センサー２１，２２は、１軸をＸ方向及びＹ方向と直交する方向、すなわち、床板Ｆの法線方向に合わせて設置される。

ただし、各センサー２１，２２を上部構造７に設置する場合、設置場所が傾いている場合もある。計測装置１は、各センサー２１，２２の３つの検出軸の１軸が、床板Ｆの法線方向に合わせて設置されなくても、概ね法線方向に向いていることで誤差は小さく無視できる。また、計測装置１は、各センサー２１，２２の３つの検出軸の１軸が、床板Ｆの法線方向に合わせて設置されなくても、ｘ軸、ｙ軸、ｚ軸の加速度を合成した３軸合成加速度によって、各センサー２１，２２の傾斜による検出誤差の補正を行うことができる。また、各センサー２１，２２は、少なくとも鉛直方向にほぼ平行な方向に生ずる加速度、あるいは、床板Ｆの法線方向の加速度を検出する１軸加速度センサーであってもよい。

以下、計測装置１が実行する本実施形態の計測方法の詳細について説明する。

１−２．車軸情報の生成
本実施形態では、計測装置１は、観測装置としてのＮ個のセンサー２１による観測情報である加速度データに基づいて、移動体である車両６の複数の部位が観測点Ｐ_ｊをそれぞれ通過した時刻及び当該複数の部位のそれぞれの観測点Ｐ_ｊへの作用に対する応答である第１の物理量を含む第１観測点情報を取得する。同様に、本実施形態では、計測装置１は、観測装置としてのＮ個のセンサー２２による観測情報である加速度データに基づいて、車両６の複数の部位が観測点Ｑ_ｊをそれぞれ通過した時刻及び当該複数の部位のそれぞれの観測点Ｑ_ｊへの作用に対する応答である第２の物理量を含む第２観測点情報を取得する。ここで、ｊは１以上Ｎ以下の各整数である。

本実施形態では、車両６が備える複数の車軸又は車輪による荷重が上部構造７に印加されるものと考え、第１観測点情報及び第２観測点情報を取得する対象となる複数の部位のそれぞれは、車軸又は車輪である。以下、本実施形態では、複数の部位のそれぞれは車軸であるものとする。

また、本実施形態では、加速度センサーである各センサー２１は、複数の車軸のそれぞれの観測点Ｐ_ｊへの作用による加速度を検出する。同様に、加速度センサーである各センサー２２は、複数の車軸のそれぞれの観測点Ｑ_ｊへの作用による加速度を検出する。

本実施形態では、図２に示したように、観測点Ｐ_１〜Ｐ_Ｎは第１の端部ＥＡ１に設定され、観測点Ｑ_１〜Ｑ_Ｎは第２の端部ＥＡ２に設定される。したがって、車両６の複数の車軸が観測点Ｐ_ｊをそれぞれ通過した時刻を、各車軸の上部構造７への進入時刻、より詳細にはレーンＬ_ｊへの進入時刻とみなすことができる。また、車両６の複数の車軸が観測点Ｑ_ｊをそれぞれ通過した時刻を、各車軸の上部構造７からの退出時刻、より詳細にはレーンＬ_ｊからの退出時刻とみなすことができる。

したがって、本実施形態では、第１観測点情報は、車両６の各車軸のレーンＬ_ｊへの進入時刻及び各車軸がレーンＬ_ｊに進入する時の作用に対する応答である第１の物理量としての加速度強度を含む。また、第２観測点情報は、車両６の各車軸のレーンＬ_ｊからの退出時刻及び各車軸がレーンＬ_ｊから退出する時の作用に対する応答である第２の物理量としての加速度強度を含む。

さらに、車両６の各車軸の進入と退出は対応するから、第１観測点情報及び第２観測点情報を層別することができ、第１観測点情報、第２観測点情報及びこれらの層別情報を含めて車軸情報と呼ぶことにする。

すなわち、車軸情報は、第１観測点情報及び第２観測点情報に加えて、車軸毎の、レーンＬ_ｊへの進入時刻、進入時の加速度強度、レーンＬ_ｊからの退出時刻及び退出時の加速度強度の対応情報や、車両６と車軸毎の当該対応情報との対応情報を含む。したがって、車軸情報により、上部構造７を通過した車両６毎に、各車軸が通過したレーンＬ_ｊや、観測点Ｐ_ｊ，Ｑ_ｊを通過した時刻及び通過時の加速度強度が特定される。

図５に車軸情報の一例を示す。図５の例では、１列目〜４列目の情報は車両番号が１の車両６に関する情報である。１列目の情報は車軸番号が１である先頭の車軸に関する情報であり、２列目の情報は車軸番号が２である２番目の車軸に関する情報であり、３列目の情報は車軸番号が３である３番目の車軸に関する情報であり、４列目の情報は車軸番号が４である４番目の車軸に関する情報である。例えば、１列目の対応情報は、車両番号が１の車両６の車軸番号が１である先頭の車軸について、レーンＬ_２への進入時刻がｔｉ１１であり、進入時の加速度強度がｐａｉ１１であり、レーンＬ_２からの退出時刻がｔｏ１１であり、退出時の加速度強度がｐａｏ１１であることを示している。

また、５列目〜６列目の情報は車両番号が２の車両６に関する情報である。５列目の情報は車軸番号が１である先頭の車軸に関する対応情報であり、６列目の情報は車軸番号が２である２番目の車軸に関する対応情報である。例えば、５列目の対応情報は、車両番号が２の車両６の車軸番号が１である先頭の車軸について、レーンＬ_１への進入時刻がｔｉ２１であり、進入時の加速度強度がｐａｉ２１であり、レーンＬ_１からの退出時刻がｔｏ２１であり、退出時の加速度強度がｐａｏ２１であることを示している。

また、７列目〜８列目の情報は車両番号が３の車両６に関する情報である。７列目の情報は車軸番号が１である先頭の車軸に関する対応情報であり、８列目の情報は車軸番号が２である２番目の車軸に関する対応情報である。例えば、７列目の対応情報は、車両番号が３の車両６の車軸番号が１である先頭の車軸について、レーンＬ_１への進入時刻がｔｉ３１であり、進入時の加速度強度がｐａｉ３１であり、レーンＬ_１からの退出時刻がｔｏ３１であり、退出時の加速度強度がｐａｏ３１であることを示している。

一例として、図６及び図７に、Ｎ＝２の場合の各センサー２１，２２及び観測点Ｐ_１，Ｐ_２，Ｑ_１，Ｑ_２の配置例を示し、図６及び図７に示す配置例の場合に、計測装置１が車軸情報を生成する手順について説明する。

図６は、上部構造７をその上方から見た図であり、図７は、図６をＡ−Ａ線又はＢ−Ｂ線で切断した断面図である。図６及び図７の例では、２個のセンサー２１が、上部構造７の第１の端部ＥＡ１において主桁Ｇ_１，Ｇ_３にそれぞれ設けられ、２個のセンサー２２が、上部構造７の第２の端部ＥＡ２において主桁Ｇ_１，Ｇ_３にそれぞれ設けられている。また、レーンＬ_１に対応する観測点Ｐ_１，Ｑ_１がそれぞれ主桁Ｇ_１に設けられたセンサー２１，２２の鉛直上方向にある床板Ｆの表面の位置に設定され、レーンＬ_２に対応する観測点Ｐ_２，Ｑ_２がそれぞれ主桁Ｇ_３に設けられたセンサー２１，２２の鉛直上方向にある床板Ｆの表面の位置に設定されている。主桁Ｇ_１に設けられたセンサー２１は観測点Ｐ_１を観測し、主桁Ｇ_３に設けられたセンサー２１は観測点Ｐ_２を観測する。また、主桁Ｇ_１に設けられたセンサー２２は観測点Ｑ_１を観測し、主桁Ｇ_３に設けられたセンサー２２は観測点Ｑ_２を観測する。計測装置１は、車軸情報を生成するために、各センサー２１，２２が検出した各時刻の加速度を振幅に変換し、加速度強度を取得する。

図８は、４軸の車両６がレーンＬ_２を走行した場合に観測点Ｐ_１，Ｐ_２，Ｑ_１，Ｑ_２に対して検出される加速度の一例を示す図である。また、図９は、図８の各時刻の加速度振幅を加速度強度に変換した図である。図８及び図９の例では、車両６がレーンＬ_２を走行しているので、車両６の４つの車軸が観測点Ｐ_２，Ｑ_２をそれぞれ通過する時刻において大きな加速度強度が取得されている。４つの車軸が観測点Ｐ_２をそれぞれ通過する時刻において取得される加速度強度は第１観測点情報に含まれる。また、４つの車軸が観測点Ｑ_２をそれぞれ通過する時刻において取得される加速度強度は第２観測点情報に含まれる。

そして、計測装置１は、取得した加速度強度が所定の閾値を超えた時刻を先頭の車軸から順に各車軸が観測点Ｐ_２，Ｑ_２を通過した時刻、すなわち、各車軸のレーンＬ_２への進入時刻及びレーンＬ_２からの退出時刻として取得する。

図１０は、図９の加速度強度を所定の閾値で２値化した図である。図１０の例では、４つの車軸のそれぞれのレーンＬ_２への進入時刻とレーンＬ_２からの退出時刻が取得されている。４つの車軸のそれぞれのレーンＬ_２への進入時刻は第１観測点情報に含まれる。また、４つの車軸のそれぞれのレーンＬ_２からの退出時刻は第２観測点情報に含まれる。

さらに、計測装置１は、４つの車軸のそれぞれのレーンＬ_２への進入時刻のパターン１と、４つの車軸のそれぞれのレーンＬ_２からの退出時刻のパターン２とを比較し、当該２つのパターンが同一の車両６の通過によるものが否かを判定する。４つの車軸の間隔は変化しないので、車両６が上部構造７を走行する速度が一定であれば、パターン１，２は一致することになる。例えば、計測装置１は、先頭の車軸の進入時刻と退出時刻を一致させるように、パターン１，２のいずれかの時刻をスライドさせ、２番目〜４番目の各車軸の進入時刻と退出時刻の差が所定の閾値以下の場合はパターン１，２が同一の車両６の通過によるものと判定し、当該差が所定の閾値よりも大きい場合はパターン１，２が２台の車両６の通過によるものと判定する。なお、計測装置１は、２台の車両６が同じ速度で１つのレーンを連行する場合において先の車両６の複数の車軸と後の車両６の複数の車軸をすべて１台の車両６の車軸と誤判定しないように、２つの連続する車軸の進入時刻又は退出時刻の間隔が規定以上の時間差である場合には、当該２つの車軸の進入時刻及び退出時刻を２つの車両６に分ければよい。

図１１は、図１０に対して、先頭の車軸の進入時刻と退出時刻を一致させるように、４つの車軸のそれぞれのレーンＬ_２からの退出時刻のパターン２をスライドさせた図である。なお、図１１は、図１０に対して横軸方向が拡大されている。図１１の例では、４つの車軸のそれぞれのレーンＬ_２への進入時刻のパターン１と、４つの車軸のそれぞれのレーンＬ_２からの退出時刻のパターン２とがほぼ一致しており、パターン１，２が同一の車両６の通過によるものと判定される。

そして、計測装置１は、図１０に示したレーンＬ_２への４つの進入時刻、図９に示した観測点Ｐ_２の４つの加速度強度のピーク値、図１０に示したレーンＬ_２からの４つの退出時刻、及び、図９に示した観測点Ｑ_２の４つの加速度強度のピーク値を、先頭から順に対応づけることにより、先頭の車軸の対応情報、２番目の車軸の対応情報、３番目の車軸の対応情報及び４番目の車軸の対応情報を取得する。さらに、計測装置１は、レーンＬ_２を走行した車両６と４つの車軸の対応情報とを対応付けた対応情報を取得する。これらの情報は、第１観測点情報及び第２観測点情報とともに、車軸情報に含まれる。

計測装置１は、車軸情報により、上部構造７のレーンＬ_ｊを通過した任意の車両６に対して、当該車両６の各車軸の観測点Ｐ_ｊへの進入時刻、各車軸による観測点Ｐ_ｊの加速度強度、各車軸の観測点Ｑ_ｊからの退出時刻、及び、各車軸による観測点Ｑ_ｊの加速度強度を特定することができる。

１−３．たわみ波形の生成
本実施形態では、橋梁５の上部構造７において、床板Ｆと主桁Ｇ_１〜Ｇ_Ｋなどで構成される橋床７ａが１つ或いは複数の連続配置される構成として考え、計測装置１は、１つの橋床７ａの変位を長手方向の中央部における変位として算出する。上部構造７に印加される荷重は上部構造７の一端から他端へ移動する。この時、荷重の上部構造７上の位置と荷重量を用いて、上部構造７の中央部の変位であるたわみ量を表すことができる。本実施形態では、車両６の車軸が上部構造７上を移動するときのたわみ変形を、１点荷重の梁上の移動によるたわみ量の軌跡として表すために、図１２に示す構造モデルを考え、当該構造モデルにおいて、中央部におけるたわみ量を算出する。図１２において、Ｐは荷重である。ａは、車両６が進入する側の上部構造７の端からの荷重位置である。ｂは、車両６が退出する側の上部構造７の端からの荷重位置である。ｌは、上部構造７の両端の間の距離である。図１２に示す構造モデルは、両端を支点とする両端を支持した単純梁である。

図１２に示す構造モデルにおいて、車両６が進入する側の上部構造７の端の位置をゼロとしてたわみ量の観測位置をｘとしたとき、単純梁の曲げモーメントＭは式（１）で表される。

式（１）において、関数Ｈ_ａは式（２）のように定義される。

式（１）を変形し、式（３）が得られる。

一方、曲げモーメントＭは式（４）で表される。式（４）において、θは角度であり、Ｉは二次モーメントであり、Ｅはヤング率である。

式（４）を式（３）に代入し、式（５）が得られる。

式（５）を観測位置ｘについて積分する式（６）を計算し、式（７）が得られる。式（７）において、Ｃ_１は積分定数である。

さらに、式（７）を観測位置ｘについて積分する式（８）を計算し、式（９）が得られる。式（９）において、Ｃ_２は積分定数である。

式（９）において、θｘはたわみ量を表し、θｘをたわみ量ｗに置き換えて式（１０）が得られる。

図１２より、ｂ＝ｌ−ａなので、式（１０）は式（１１）のように変形される。

ｘ＝０でたわみ量ｗ＝０として、ｘ≦ａよりＨ_ａ＝０であるから、式（１１）にｘ＝ｗ＝Ｈ_ａ＝０を代入して整理すると、式（１２）が得られる。

また、ｘ＝ｌでたわみ量ｗ＝０として、ｘ＞ａよりＨ_ａ＝１であるから、式（１１）にｘ＝ｌ，ｗ＝０，Ｈ_ａ＝１を代入して整理すると、式（１３）が得られる。

式（１３）にｂ＝ｌ−ａを代入し、式（１４）が得られる。

式（１０）に式（１２）の積分定数Ｃ_１及び式（１３）の積分定数Ｃ_２を代入し、式（１５）が得られる。

式（１５）を変形し、荷重Ｐが位置ａに印加された時の観測位置ｘにおけるたわみ量ｗは、式（１６）で表される。

図１３に、たわみ量の観測位置ｘが単純梁の中央位置に固定されている条件で、すなわちｘ＝ｌ／２の時に、荷重Ｐが単純梁の一端から他端へ移動する様子を示す。

荷重位置ａが観測位置ｘ＝ｌ／２よりも左側にある時、ｘ＞ａよりＨ_ａ＝１であるから、式（１６）にｘ＝ｌ／２，Ｈ_ａ＝１を代入し、式（１７）が得られる。

式（１７）にｌ＝ａ＋ｂを代入して整理すると式（１８）が得られる。

式（１８）にａ＋ｂ＝ｌを代入し、荷重Ｐの位置が中央の観測位置ｘ＝ｌ／２よりも左側にある場合の観測位置ｘのたわみ量ｗ_Ｌは式（１９）のようになる。

一方、荷重位置ａが観測位置ｘ＝ｌ／２よりも右側にある時、ｘ≦ａよりＨ_ａ＝０であるから、式（１６）にｘ＝ｌ／２，Ｈ_ａ＝０を代入し、式（２０）が得られる。

式（２０）にｌ＝ａ＋ｂを代入して整理すると、荷重Ｐの位置が中央の観測位置ｘ＝ｌ／２よりも右側にある場合の観測位置ｘのたわみ量ｗ_Ｒは式（２１）のようになる。

また、荷重位置ａが観測位置ｘ＝ｌ／２と同じである時、ｘ≦ａよりＨ_ａ＝０であるから、式（１６）にＨ_ａ＝０，ａ＝ｂ＝ｌ／２を代入して整理すると、式（２２）が得られる。

さらに、式（２２）にａ＝ｌ／２を代入すると、荷重Ｐの位置が中央の観測位置と同じである場合の観測位置ｘのたわみ量ｗは式（２３）のようになる。

両端支点の単純梁では、荷重Ｐが中央にある場合が最大たわみ変位となるので、式（２３）より、最大たわみ量ｗ_ｍａｘは式（２４）で表される。

荷重Ｐの位置が中央の観測位置ｘ＝ｌ／２よりも左側にある場合の観測位置ｘのたわみ量ｗ_Ｌを最大たわみ量ｗ_ｍａｘで除算して最大たわみ量ｗ_ｍａｘで規格化すると、式（１９）及び式（２４）より、式（２５）が得られる。

式（２５）においてａ／ｌ＝ｒと置くと式（２６）が得られる。

一方、荷重Ｐの位置が中央の観測位置ｘ＝ｌ／２よりも右側にある場合の観測位置ｘのたわみ量ｗ_Ｒを最大たわみ量ｗ_ｍａｘで除算して最大たわみ量ｗ_ｍａｘで規格化すると、式（２１）及び式（２４）より、式（２７）が得られる。

ここで、ａ／ｌ＝ｒ，ａ＋ｂ＝ｌより、ｂ＝ｌ×（１−ｒ）であるから、式（２７）にｂ＝ｌ×（１−ｒ）を代入し、式（２８）が得られる。

式（２５）、式（２７）とまとめて、単純梁上を荷重Ｐが移動する際に中央部で観測される最大たわみ量で規格化された規格化たわみ量ｗ_ｓｔｄは、式（２９）で表される。

式（２９）において、ｒ＝ａ／ｌ，１−ｒ＝ｂ／ｌは、単純梁の支点間の距離ｌに対する荷重Ｐの位置の比を示し、式（３０）に示すように変数Ｒを定義する。

式（３０）を用いて式（２９）は式（３１）に置き換えられる。

式（３０）及び式（３１）は、単純梁の中央に観測位置がある場合、荷重Ｐの位置が中央よりも右側と左側でたわみ量が対称となることを示している。

図１４に、観測位置ｘ＝ｌ／２の場合の規格化たわみ量ｗ_ｓｔｄの波形の一例を示す。図１４において、横軸は荷重Ｐの位置であり、縦軸は規格化たわみ量ｗ_ｓｔｄである。図１４の例では、単純梁の支点間の距離ｌ＝１である。

前述の車軸情報に含まれるのは、車両６の各車軸のレーンＬ_ｊへの進入時刻及びレーンＬ_ｊからの退出時刻、すなわち、車両６が上部構造７の両端の位置をそれぞれ通過した時刻であるので、上部構造７の両端の位置を車軸の進入時刻及び退出時刻に対応させて、荷重位置ａ，ｂを時間に置き換える。ただし、車両６の速度はほぼ一定で、位置と時刻はほぼ比例するものとする。

上部構造７の左端の荷重位置を進入時刻ｔ_ｉと対応させ、上部構造７の右端の荷重位置を退出時刻ｔ_ｏと対応させると、左端からの荷重位置ａは進入時刻ｔ_ｉからの経過時刻ｔ_ｐに置き換える。経過時刻ｔ_ｐは式（３２）で表される。

また、支点間の距離ｌは、進入時刻ｔ_ｉから退出時刻ｔ_ｏまでの時間ｔ_ｓに置き換えられる。時間ｔ_ｓは式（３３）で表される。

車両６の速度は一定なので、荷重位置ａが上部構造７の中央にある時刻ｔ_ｃは式（３４）で表される。

以上のように位置を時間に置き換えて、荷重Ｐの位置は、式（３５）及び式（３６）のようになる。

式（３５）及び式（３６）を式（２９）に代入し、時間に置き換えた規格化たわみ量ｗ_ｓｔｄは、式（３７）で表される。

あるいは、式（３０）及び式（３１）より、変数Ｒを時間に置き換えて、最大振幅で規格化された規格化たわみ量ｗ_ｓｔｄは、式（３８）で表される。

時間経過と規格化たわみ量の関係付けを観測データとして扱うことを考慮して、規格化たわみ量ｗ_ｓｔｄを、両端支点の単純梁上の単一集中荷重の移動による梁中央の観測位置の規格化たわみ量モデルｗ_ｓｔｄ（ｔ）に置き換えて、式（３８）は式（３９）のようになる。式（３９）は、構造物である上部構造７のたわみの近似式であり、上部構造７の構造モデルに基づく式である。具体的には、式（３９）は、車両６が移動するレーンＬ_ｊにおける観測点Ｐ_ｊと観測点Ｑ_ｊとの中央位置におけるたわみの最大振幅で規格化され、最大値が１となる式である。

この規格化たわみ量モデルｗ_ｓｔｄ（ｔ）に必要な時間情報は、前述の車軸情報から得られる。規格化たわみ量モデルｗ_ｓｔｄ（ｔ）は、上部構造７の中央位置において最大たわみ量ｗ_ｍａｘとなるので、式（４０）が得られる。

また、前述の式（２３）に示されるたわみ量ｗは、荷重Ｐの位置が中央の観測位置と同じである場合の観測位置ｘ＝ｌ／２のたわみ量であり、最大たわみ量ｗ_ｍａｘと一致するので、式（４１）が得られる。

図１５に、規格化たわみ量モデルｗ_ｓｔｄ（ｔ）の一例を示す。図１５の例では、進入時刻ｔ_ｉ＝４、退出時刻ｔ_ｏ＝６であり、時刻ｔ_ｃ＝（ｔ_ｉ＋ｔ_ｏ）／２＝５において、規格化たわみ量モデルｗ_ｓｔｄ（ｔ）は、上部構造７の中央位置において最大たわみ量ｗ_ｍａｘ＝１となっている。

仮に、構造物である上部構造７は、ＢＷＩＭ（Bridge Weigh in Motion）が機能すると仮定し、両端を支点とした単純梁に近似した変形をすると考える。また、移動体である車両６は上部構造７の一方の端部から上部構造７を概ね一定速度で通過し、もう一方の端部に移動するのであるから、上部構造７の中間部と上部構造７の端部は同じ荷重の作用を受けるので、観測される上部構造７の変位は、車軸情報から得られる車軸の加速度強度ａ_ｐに近似的に比例すると考えることができる。

比例係数を、車軸情報から得られる車軸の加速度強度ａ_ｐと所定の係数ｐとの積であるものとして、式（４２）により、各車軸による上部構造７のたわみ波形Ｈ（ｔ）が得られる。なお、加速度強度ａ_ｐは、車軸情報に含まれる進入時の加速度強度でもよいし、退出時の加速度強度でもよいし、進入時の加速度強度と退出時の加速度強度との平均値等の統計値でもよい。

式（３９）を式（４２）に代入し、たわみ波形Ｈ（ｔ）は式（４３）で表される。

これまでは、上部構造７に単一荷重Ｐが印加されるものとしていたが、車両６が走行するレーンＬ_ｊには車両６の各車軸による荷重が印加されるので、式（４３）は式（４４）のようにたわみ波形Ｈ_ｊｋ（ｔ）に置き換えられる。式（４４）において、ｋは車軸番号を表す整数であり、ｊはレーン番号を表す整数である。式（４４）に示すように、たわみ波形Ｈ_ｊｋ（ｔ）は、所定の係数ｐと加速度強度ａ_ｐｊｋとの積に比例する。

図１６に、レーンＬ_ｊを走行する車両６に含まれる各車軸による上部構造７のたわみ波形の一例を示す。図１６の例では、車両６は４軸車両であり、４つのたわみ波形Ｈ_ｊ１（ｔ），Ｈ_ｊ２（ｔ），Ｈ_ｊ３（ｔ），Ｈ_ｊ４（ｔ）が示されている。図１６の例では、先頭と２番目の車軸による荷重は相対的に小さく、３番目と４番目の車軸による荷重は相対的に大きいため、たわみ波形Ｈ_ｊ１（ｔ），Ｈ_ｊ２（ｔ）の最大振幅は相対的に小さく、たわみ波形Ｈ_ｊ３（ｔ），Ｈ_ｊ４（ｔ）の最大振幅は相対的に大きくなっている。

１−４．補正係数の算出
図１７に、４軸車、３軸車、２軸車の重心位置及び車軸の位置の一例を示す。図１７の例では、４軸車は、２番目の車軸と３番目の車軸の間に重心がある。また、３軸車は、先頭の車軸と２番目の車軸の間に重心がある。また、２軸車は、先頭の車軸と２番目の車軸の間に重心がある。

車両６の各車軸の加速度強度ａ_ｐｊｋは、車両６の重心から各車軸までの距離モーメントの影響を受けるため、静止時の各車軸の荷重によって生じる加速度強度とは大きく異なる場合がある。その場合、式（４４）によって算出されるたわみ波形Ｈ_ｊｋ（ｔ）と実際に計測される変位波形との近似性が低下する。そこで、本実施形態では、計測装置１は、第１観測点情報に基づいて、式（４４）のたわみ波形Ｈ_ｊｋ（ｔ）を、車両６の重心から各車軸までの距離モーメントの影響を考慮したたわみ波形に置き換える。

車両６の複数の車軸のうちの車軸番号が１である先頭の車軸が観測点Ｐ_ｊを通過した時刻である先頭時刻をｔ_ｊ１とすると、計測装置１は、式（４５）のように、先頭時刻ｔ_ｊ１から車軸番号がｋである車軸が観測点Ｐ_ｊを通過した時刻ｔ_ｊｋまでの時間Δｔ_ｊｋを算出する。ｋは１以上ｌａｓｔ以下の各整数である。なお、先頭時刻ｔ_ｊ１は「第１の先頭時刻」の一例である。

図１８に、式（４５）によって算出される時間Δｔ_ｊ２〜Δｔ_{ｊｌａｓｔ}の一例を示す。時間Δｔ_ｊ１はゼロであるため、図示が省略されている。なお、図１８は、４軸車、３軸車、２軸車がそれぞれレーンＬ_ｊを走行したときの各車軸による観測点Ｐ_ｊの加速度強度を２値化した図である。

次に、計測装置１は、先頭時刻ｔ_ｊ１から最後の車軸が観測点Ｐ_ｊを通過した時刻までの時間をΔｔ_{ｊｌａｓｔ}として、式（４６）により、式（４５）で表される時間Δｔ_ｊｋを時間Δｔ_{ｊｌａｓｔ}で規格化した規格化時間Δｔ_{ｊｓｔｄｋ}を算出する。

次に、計測装置１は、式（４７）により、車軸番号が１である先頭の車軸から車軸番号がｋである車軸までの各車軸による観測点Ｐ_ｊの加速度強度ａ_ｐｊ１〜ａ_ｐｊｋの積算値を、車軸番号が１である先頭の車軸から車軸番号がｌａｓｔである最後の車軸までの各車軸による観測点Ｐ_ｊの加速度強度ａ_ｐｊ１〜ａ_{ｐｊｌａｓｔ}の総和で規格化した加算規格化衝撃パワー_ａｄｄ｛ａ_{ｐｊ＿ｓｔｄｋ}｝を算出する。

図１９に、式（４６）によって算出される規格化時間Δｔ_{ｊｓｔｄｋ}と、式（４７）によって算出される加算規格化衝撃パワー_ａｄｄ｛ａ_{ｐｊ＿ｓｔｄｋ}｝との相関の一例を示す。図１９において、横軸は規格化時間Δｔ_{ｊｓｔｄｋ}であり、縦軸は加算規格化衝撃パワー_ａｄｄ｛ａ_{ｐｊ＿ｓｔｄｋ}｝である。図１９の例では、車両６は４軸車であり、４つの点がプロットされている。破線は２つの点を結ぶ直線である。図１９に示すように、加算規格化衝撃パワー_ａｄｄ｛ａ_{ｐｊ＿ｓｔｄｋ}｝は、規格化時間Δｔ_{ｊｓｔｄｋ}に対して単調増加する。

次に、車両６の荷重を、車両６の重心よりも前の荷重と重心よりも後の荷重に配分する所定の配分比をＲ：１−Ｒ＝_ｉｍａｇＲ_Ｆ：（１−_ｉｍａｇＲ_Ｆ）とし、計測装置１は、式（４８）のように、加算規格化衝撃パワー_ａｄｄ｛ａ_{ｐｊ＿ｓｔｄｋ}｝から_ｉｍａｇＲ_Ｆを減算したシフト加算規格化衝撃パワー_{ｓｈｉｆｔ＿ａｄｄ}｛ａ_{ｐｊ＿ｓｔｄｋ}｝を算出する。

この配分比_ｉｍａｇＲ_Ｆ：（１−_ｉｍａｇＲ_Ｆ）は、計測できない未知数であるので、例えば、一般的な車両の操縦安定性から推奨される重心前後の荷重の配分比が想定される。例えば、この配分比を０．４：０．６とすると、_ｉｍａｇＲ_Ｆは０．４である。なお、一般的な車両の操縦安定性から想定される_ｉｍａｇＲ_Ｆは０．４以上０．６以下である。

図２０に、図１９に示される規格化時間Δｔ_{ｊｓｔｄｋ}と加算規格化衝撃パワー_ａｄｄ｛ａ_{ｐｊ＿ｓｔｄｋ}｝から_ｉｍａｇＲ_Ｆ＝０．４を減算したシフト加算規格化衝撃パワー_{ｓｈｉｆｔ＿ａｄｄ}｛ａ_{ｐｊ＿ｓｔｄｋ}｝との相関を示す。図２０において、横軸は規格化時間Δｔ_{ｊｓｔｄｋ}であり、縦軸はシフト加算規格化衝撃パワー_{ｓｈｉｆｔ＿ａｄｄ}｛ａ_{ｐｊ＿ｓｔｄｋ}｝である。

次に、計測装置１は、規格化時間Δｔ_{ｊｓｔｄｋ}とシフト加算規格化衝撃パワー_{ｓｈｉｆｔ＿ａｄｄ}｛ａ_{ｐｊ＿ｓｔｄｋ}｝との相関から、シフト加算規格化衝撃パワー_{ｓｈｉｆｔ＿ａｄｄ}｛ａ_{ｐｊ＿ｓｔｄｋ}｝がゼロとなる基準時間Δｔ_ｊＧＯＣを内挿する。

次に、計測装置１は、車軸番号が１〜ｓの車軸に対するシフト加算規格化衝撃パワー_{ｓｈｉｆｔ＿ａｄｄ}｛ａ_{ｐｊ＿ｓｔｄｓ}｝が負の値であり、車軸番号がｓ＋１〜ｌａｓｔの車軸に対するシフト加算規格化衝撃パワー_{ｓｈｉｆｔ＿ａｄｄ}｛ａ_{ｐｊ＿ｓｔｄｓ}｝が正の値である車軸番号ｓを求める。

式（４９）のように、車軸番号がｓの車軸に対して、規格化時間Δｔ_{ｊｓｔｄｓ}と負の値であるシフト加算規格化衝撃パワー_{ｓｈｉｆｔ＿ａｄｄ}｛ａ_{ｐｊ＿ｓｔｄｓ}｝を示す座標を（ｘ_ｓ，ｙ_ｓ）とする。

同様に、式（５０）のように、車軸番号がｓ＋１の車軸に対して、規格化時間Δｔ_{ｊｓｔｄｓ＋１}と正の値であるシフト加算規格化衝撃パワー_{ｓｈｉｆｔ＿ａｄｄ}｛ａ_{ｐｊ＿ｓｔｄｓ＋１}｝を示す座標を（ｘ_ｓ＋１，ｙ_ｓ＋１）とする。

座標（ｘ_ｓ，ｙ_ｓ）の点と座標（ｘ_ｓ＋１，ｙ_ｓ＋１）の点を通る直線は式（５１）で与えられる。

計測装置１は、式（５２）により、式（５１）で表される直線とｙ＝０との交点のｘ座標を算出し、算出したｘ座標を基準時間Δｔ_ｊＧＯＣとする。

ここで、車軸番号が１である先頭の車軸から車軸番号がｌａｓｔである最後の車軸までの加速度強度ａ_ｐｊ１〜ａ_{ｐｊｌａｓｔ}の総和を所定の配分比_ｉｍａｇＲ_Ｆ：（１−_ｉｍａｇＲ_Ｆ）で配分する時刻を基準時刻ｔ_ｊＧＯＣとすると、基準時間Δｔ_ｊＧＯＣは、先頭の車軸が観測点Ｐ_ｊを通過した時刻から基準時刻ｔ_ｊＧＯＣまでの時間である。基準時刻ｔ_ｊＧＯＣは、車両６の重心が観測点Ｐ_ｊを通過したと推定される時刻である。

図２１に、図２０に示される規格化時間Δｔ_{ｊｓｔｄｋ}とシフト加算規格化衝撃パワー_{ｓｈｉｆｔ＿ａｄｄ}｛ａ_{ｐｊ＿ｓｔｄｋ}｝との相関に対して算出される基準時間Δｔ_ｊＧＯＣと基準時刻ｔ_ｊＧＯＣとの関係を示す。本実施形態では、先頭の車軸が観測点Ｐ_ｊを通過した時刻をゼロとしているので、基準時間Δｔ_ｊＧＯＣと基準時刻ｔ_ｊＧＯＣとは等しい値である。なお、基準時刻ｔ_ｊＧＯＣは「第１の基準時刻」の一例である。

次に、計測装置１は、式（５３）のように、各車軸の規格化時間Δｔ_{ｊｓｔｄｋ}と基準時間Δｔ_ｊＧＯＣとの時間差である規格化時間差Δｔ_{ｊｓｔｄ−ＧＯＣｋ}を算出する。

次に、計測装置１は、式（５４）により、シフト加算規格化衝撃パワー_{ｓｈｉｆｔ＿ａｄｄ}｛ａ_{ｐｊ＿ｓｔｄｋ}｝が負の値である車軸番号が１〜ｓの各車軸の規格化時間差Δｔ_{ｊｓｔｄ−ＧＯＣｋ}と、車軸番号が１〜ｓの車軸の規格化時間差Δｔ_{ｊｓｔｄ−ＧＯＣ１}〜Δｔ_{ｊｓｔｄ−ＧＯＣｓ}の和との比を算出し、補正係数Ｒ_ｊｋとする。ｋは１以上ｓ以下の各整数である。

同様に、計測装置１は、式（５５）により、シフト加算規格化衝撃パワー_{ｓｈｉｆｔ＿ａｄｄ}｛ａ_{ｐｊ＿ｓｔｄｋ}｝が正の値である車軸番号がｓ＋１〜ｌａｓｔの各車軸の規格化時間差Δｔ_{ｊｓｔｄ−ＧＯＣｋ}と、車軸番号がｓ＋１〜ｌａｓｔの車軸の規格化時間差Δｔ_{ｊｓｔｄ−ＧＯＣｓ＋１}〜Δｔ_{ｊｓｔｄ−ＧＯＣｌａｓｔ}の和との比を算出し、補正係数Ｒ_ｊｋとする。ｋはｓ＋１からｌａｓｔ以下の各整数である。

式（５４）及び式（５５）により、車軸番号が１〜ｌａｓｔの車軸の補正係数Ｒ_ｊ１〜Ｒ_{ｊｌａｓｔ}が得られる。補正係数Ｒ_ｊｋは、距離モーメントの影響を低減させるように加速度強度ａ_ｐｊｋを補正する係数であり、たわみ波形Ｈ_ｊｋ（ｔ）を算出するための前述の式（４４）は、補正係数Ｒ_ｊｋを用いて式（５６）に置き換えられる。

計測装置１は、式（５６）により、車軸番号がｋである車軸による上部構造７のたわみ波形Ｈ_ｊｋ（ｔ）を算出する。

図２２に、式（５６）によって算出される各車軸による上部構造７のたわみ波形の一例を示す。図２２には、図１６に示した４つのたわみ波形がそれぞれ補正されたたわみ波形Ｈ_ｊ１（ｔ），Ｈ_ｊ２（ｔ），Ｈ_ｊ３（ｔ），Ｈ_ｊ４（ｔ）が示されている。図２２の例を図１６の例と比較すると、たわみ波形Ｈ_ｊ１（ｔ），Ｈ_ｊ２（ｔ）の最大振幅は増加し、たわみ波形Ｈ_ｊ３（ｔ），Ｈ_ｊ４（ｔ）の最大振幅は減少している。すなわち、補正係数Ｒ_ｊ１〜Ｒ_{ｊｌａｓｔ}によって、加速度強度ａ_ｐｊ１〜ａ_{ｐｊｌａｓｔ}に対する距離モーメントの影響が補正されている。

式（５７）に示すように、レーンＬ_ｊを走行する車両６による上部構造７のたわみ波形ＣＰ_ｊｍ（ｔ）は、各車軸による上部構造７のたわみ波形Ｈ_ｊｋ（ｔ）を加算して得られる。式（５７）において、ｍは車両番号を表す整数であり、ｋは車軸番号を表す整数であり、ｊはレーン番号を表す整数である。

計測装置１は、式（５７）により、レーンＬ_ｊを走行する車両番号がｍである車両６による上部構造７のたわみ波形ＣＰ_ｊｍ（ｔ）を算出する。

図２３に、図２２に示した４つのたわみ波形Ｈ_ｊ１（ｔ），Ｈ_ｊ２（ｔ），Ｈ_ｊ３（ｔ），Ｈ_ｊ４（ｔ）を加算して得られる、車両番号がｍの車両６による上部構造７のたわみ波形ＣＰ_ｊｍ（ｔ）を実線で示す。また、図２３には、図１６に示した４つのたわみ波形Ｈ_ｊ１（ｔ），Ｈ_ｊ２（ｔ），Ｈ_ｊ３（ｔ），Ｈ_ｊ４（ｔ）を加算して得られるたわみ波形ＣＰ_ｊｍ（ｔ）も破線で示されている。図２３に示すように、補正係数Ｒ_ｊ１〜Ｒ_{ｊｌａｓｔ}によって、たわみ波形ＣＰ_ｊｍ（ｔ）の最大振幅や最大振幅となる時刻が補正されている。

１−５．変位の算出
式（５７）に示されるたわみ波形ＣＰ_ｊｍ（ｔ）と、観測した変位ＣＵ（ｔ）の相関を多項式で近似する。例えば、式（５８）のように、変位ＣＵ（ｔ）をたわみ波形ＣＰ_ｊｍ（ｔ）の１次式で近似する。式（５８）において、ｓは１次係数であり、ｉは０次係数である。

１次係数ｓ及び０次係数ｉは、例えば、複数の試験車両による荷重試験によって算出される。例えば、レーンＬ_ｊの中央に変位計を設置し、複数の試験車両をそれぞれ単独でレーンＬ_ｊを走行させ、計測装置１は、車軸情報を生成するとともに、変位計によって計測された変位を取得する。そして、計測装置１は、計測された変位の最大値をＣＵ_ｍａｘとし、車軸情報から式（５７）によって得られるたわみ波形ＣＰ_ｊｍ（ｔ）の最大値をＣＰ_{ｊｍ−ｍａｘ}としてグラフにプロットし、近似直線の１次係数ｓ_ｃｕ及び０次係数ｉ_ｃｕを求める。

図２４は、６台の試験車両による荷重試験の結果をプロットした図である。図２４において、横軸は試験車両による上部構造７のたわみ波形の最大値ＣＰ_{ｊｍ−ｍａｘ}であり、縦軸は計測変位の最大値ＣＵ_ｍａｘである。図２４において、６つの点は直線状に配置されており、６つの点に対する近似直線が点線で示されている。図２４の例では、この近似直線の１次係数ｓ_ｃｕは３０８４．４３５９４４であり、０次係数ｉ_ｃｕは０．２２９１８０１７４である。

計測装置１は、１次係数ｓ_ｃｕ及び０次係数ｉ_ｃｕと、未知の車両６の車軸情報から式（５７）によって得られるたわみ波形ＣＰ_ｊｍ（ｔ）とを用いて、式（５９）によりレーンＬ_ｊの中央の変位ＣＵ_ｅｓｔ（ｔ）を算出する。

図２５に、計測変位ＣＵ（ｔ）と、式（５７）によって得られる補正されたたわみ波形ＣＰ_ｊｍ（ｔ）から算出される変位ＣＵ_ｅｓｔ（ｔ）の一例を示す。図２５において、実線は計測変位ＣＵ（ｔ）を示し、破線はたわみ波形ＣＰ_ｊｍ（ｔ）から算出される変位ＣＵ_ｅｓｔ（ｔ）を示す。なお、図２５には、式（４４）によって算出される各車軸による上部構造７のたわみ波形Ｈ_ｊｋ（ｔ）を加算して得られるたわみ波形ＣＰ_ｊｍ（ｔ）、すなわち、補正されていないたわみ波形ＣＰ_ｊｍ（ｔ）から算出される変位ＣＵ_ｅｓｔ（ｔ）も一点鎖線で示されている。図２５に示すように、補正されていないたわみ波形ＣＰ_ｊｍ（ｔ）から算出される変位ＣＵ_ｅｓｔ（ｔ）は計測変位ＣＵ（ｔ）との差が大きいが、補正されたたわみ波形ＣＰ_ｊｍ（ｔ）から算出される変位ＣＵ_ｅｓｔ（ｔ）は計測変位ＣＵ（ｔ）との差が非常に小さい。したがって、計測装置１は、レーンＬ_ｊの中央の変位を計測することなく、式（５９）により、レーンＬ_ｊを走行した未知の車両６による中央の変位を精度よく算出することができる。

式（５９）において、０次係数ｉ_ｃｕは小さい値であり、式（５９）にｉ_ｃｕ＝０を代入し、式（５６）及び式（５７）より、式（６０）が得られる。

式（６０）より、所定の係数ｐと１次係数ｓ_ｃｕは入れ替えることができるので、所定の係数ｐは１次係数ｓ_ｃｕと同様の機能を有する係数である。すなわち、所定の係数ｐは、観測点Ｐ_ｊと観測点Ｑ_ｊとの間の上部構造７の部位のたわみと、上部構造７の当該部位の変位との相関を近似する関数の係数である。

１−６．計測方法
図２６は、第１実施形態の計測方法の手順の一例を示すフローチャート図である。本実施形態では、計測装置１が図２６に示す手順を実行する。

図２６に示すように、まず、計測装置１は、観測点Ｐ_１〜Ｐ_Ｎを観測するＮ個のセンサー２１による観測情報に基づいて、車両６の複数の車軸が観測点Ｐ_１〜Ｐ_Ｎのいずれかを通過した時刻及び複数の車軸のそれぞれの観測点Ｐ_１〜Ｐ_Ｎのいずれかへの作用に対する応答である第１の物理量としての加速度強度を含む第１観測点情報を取得する（ステップＳ１）。前述の通り、Ｎ個のセンサー２１はそれぞれ加速度センサーであり、Ｎ個のセンサー２１による観測情報は、観測点Ｐ_１〜Ｐ_Ｎに生じた加速度の検出情報である。計測装置１は、Ｎ個のセンサー２１がそれぞれ検出した加速度に基づいて第１観測点情報を取得する。このステップＳ１は、第１観測点情報取得ステップである。

次に、計測装置１は、観測点Ｑ_１〜Ｑ_Ｎを観測するＮ個のセンサー２２による観測情報に基づいて、車両６の複数の車軸が観測点Ｑ_１〜Ｑ_Ｎのいずれかを通過した時刻及び複数の車軸のそれぞれの観測点Ｑ_１〜Ｑ_Ｎのいずれかへの作用に対する応答である第２の物理量としての加速度強度を含む第２観測点情報を取得する（ステップＳ２）。前述の通り、Ｎ個のセンサー２２はそれぞれ加速度センサーであり、Ｎ個のセンサー２２による観測情報は、観測点Ｑ_１〜Ｑ_Ｎに生じた加速度の検出情報である。計測装置１は、Ｎ個のセンサー２２がそれぞれ検出した加速度に基づいて、第２観測点情報を取得する。このステップＳ２は、第２観測点情報取得ステップである。

次に、計測装置１は、ステップＳ１で取得した第１観測点情報を用いて、複数の車軸のうちの先頭の車軸が観測点Ｐ_ｊを通過した先頭時刻ｔ_ｊ１から複数の車軸が観測点Ｐ_ｊをそれぞれ通過した時刻ｔ_ｊ１〜ｔ_{ｊｌａｓｔ}までの時間である規格化時間Δｔ_{ｊｓｔｄ１}〜Δｔ_{ｊｓｔｄｌａｓｔ}と、先頭時刻ｔ_ｊ１から複数の車軸による加速度強度ａ_ｐｊ１〜ａ_{ｐｊｌａｓｔ}の総和を所定の配分比_ｉｍａｇＲ_Ｆ：（１−_ｉｍａｇＲ_Ｆ）で配分する時刻である基準時刻ｔ_ｊＧＯＣまでの時間である基準時間Δｔ_ｊＧＯＣとを算出する（ステップＳ３）。具体的には、計測装置１は、前述の式（４５）〜式（５２）により、規格化時間Δｔ_{ｊｓｔｄ１}〜Δｔ_{ｊｓｔｄｌａｓｔ}と基準時間Δｔ_ｊＧＯＣとを算出する。このステップＳ３は、第１時間算出ステップである。

次に、計測装置１は、ステップＳ３で算出した規格化時間Δｔ_{ｊｓｔｄ１}〜Δｔ_{ｊｓｔｄｌａｓｔ}及び基準時間Δｔ_ｊＧＯＣに基づいて、複数の車軸による加速度強度ａ_ｐｊ１〜ａ_{ｐｊｌａｓｔ}を補正する補正係数Ｒ_ｊ１〜Ｒ_{ｊｌａｓｔ}を算出する（ステップＳ４）。具体的には、計測装置１は、前述の式（５３）〜式（５５）により、補正係数Ｒ_ｊ１〜Ｒ_{ｊｌａｓｔ}を算出する。このステップＳ４は、補正係数算出ステップである。

次に、計測装置１は、ステップＳ１で取得した第１観測点情報及びステップＳ２で取得した第２観測点情報と、所定の係数ｐと、ステップＳ４で算出した補正係数Ｒ_ｊｋと、上部構造７のたわみの近似式とに基づいて、複数の車軸のそれぞれによる上部構造７のたわみ波形Ｈ_ｊｋ（ｔ）を算出する（ステップＳ５）。具体的には、計測装置１は、前述の式（５６）により、各レーンＬ_ｊを走行した車両６の各車軸による上部構造７のたわみ波形Ｈ_ｊｋ（ｔ）を算出する。このステップＳ５は、たわみ波形算出ステップである。

次に、計測装置１は、前述の式（５７）により、ステップＳ５で算出した車両６の複数の車軸のそれぞれによる上部構造７のたわみ波形Ｈ_ｊｋ（ｔ）を加算して、車両６による上部構造７のたわみ波形ＣＰ_ｊｍ（ｔ）を算出する（ステップＳ６）。このステップＳ６は、移動体たわみ波形算出ステップである。

次に、計測装置１は、前述の式（５９）により、ステップＳ６で算出した車両６による上部構造７のたわみ波形ＣＰ_ｊｍ（ｔ）に基づいて、上部構造７の変位ＣＵ_ｅｓｔ（ｔ）を算出する（ステップＳ７）。このステップＳ７は、変位算出ステップである。

次に、計測装置１は、ステップＳ７で算出した上部構造７の変位ＣＵ_ｅｓｔ（ｔ）をサーバー２に出力する（ステップＳ８）。このステップＳ８は、出力ステップである。

計測装置１は、計測を終了するまで（ステップＳ９のＮ）、ステップＳ１〜Ｓ８の処理を繰り返し行う。

図２７は、図２６のステップＳ３である第１時間算出ステップの手順の一例を示すフローチャート図である。

図２７に示すように、まず、計測装置１は、整数ｊを１に設定し（ステップＳ３０）、第１観測点情報を用いて、レーンＬ_ｊを走行した車両６の有無を判定する（ステップＳ３１）。

そして、計測装置１は、ステップＳ３１においてレーンＬ_ｊを走行した車両６があると判定した場合は（ステップＳ３２のＹ）、第１観測点情報を用いて、前述の式（４５）により、時間Δｔ_ｊ１〜Δｔ_{ｊｌａｓｔ}を算出する（ステップＳ３３）。すなわち、計測装置１は、車両６の先頭の車軸が観測点Ｐ_ｊを通過した先頭時刻ｔ_ｊ１から複数の車軸が観測点Ｐ_ｊをそれぞれ通過した時刻ｔ_ｊ１〜ｔ_{ｊｌａｓｔ}までの時間Δｔ_ｊ１〜Δｔ_{ｊｌａｓｔ}を算出する。

次に、計測装置１は、前述の式（４６）により、ステップＳ３３で算出した時間Δｔ_ｊ１〜Δｔ_{ｊｌａｓｔ}を時間Δｔ_{ｊｌａｓｔ}で除算して規格化時間Δｔ_{ｊｓｔｄ１}〜Δｔ_{ｊｓｔｄｌａｓｔ}を算出する（ステップＳ３４）。

次に、計測装置１は、前述の式（４７）により、第１観測点情報の加速度強度ａ_ｐｊ１〜ａ_ｐｊｋの積算値を、加速度強度ａ_ｐｊ１〜ａ_{ｐｊｌａｓｔ}の総和で除算し、加算規格化衝撃パワー_ａｄｄ｛ａ_{ｐｊ＿ｓｔｄ１}｝〜_ａｄｄ｛ａ_{ｐｊ＿ｓｔｄｌａｓｔ}｝を算出する（ステップＳ３５）。

次に、計測装置１は、前述の式（４８）により、ステップＳ３５で算出した加算規格化衝撃パワー_ａｄｄ｛ａ_{ｐｊ＿ｓｔｄ１}｝〜_ａｄｄ｛ａ_{ｐｊ＿ｓｔｄｌａｓｔ}｝から_ｉｍａｇＲ_Ｆを減算してシフト加算規格化衝撃パワー_{ｓｈｉｆｔ＿ａｄｄ}｛ａ_{ｐｊ＿ｓｔｄ１}｝〜_{ｓｈｉｆｔ＿ａｄｄ}｛ａ_{ｐｊ＿ｓｔｄｌａｓｔ}｝を算出する（ステップＳ３６）。

次に、計測装置１は、前述の式（４９）〜式（５２）により、座標（Δｔ_{ｊｓｔｄｓ}，_{ｓｈｉｆｔ＿ａｄｄ}｛ａ_{ｐｊ＿ｓｔｄｓ}｝（＜０））と座標（Δｔ_{ｊｓｔｄｓ＋１}，_{ｓｈｉｆｔ＿ａｄｄ}｛ａ_{ｐｊ＿ｓｔｄｓ＋１}｝（＞０））とを通直線とｙ＝０との交点のｘ座標を基準時間Δｔ_ｊＧＯＣとする（ステップＳ３７）。

計測装置１は、ステップＳ３１においてレーンＬ_ｊを走行した車両６がないと判定した場合は（ステップＳ３２のＮ）、ステップＳ３３〜Ｓ３７の処理を行わない。

計測装置１は、整数ｊがＮでない場合は（ステップＳ３８のＮ）、整数ｊに１を加算し（ステップＳ３９）、ステップＳ３１〜Ｓ３７の処理を繰り返し行う。

そして、整数ｊがＮになると（ステップＳ３８のＹ）、計測装置１は第１時間算出ステップの処理を終了する。

図２８は、図２６のステップＳ４である補正係数算出ステップの手順の一例を示すフローチャート図である。

図２８に示すように、まず、計測装置１は、整数ｊを１に設定し（ステップＳ４１）、図２７のステップＳ３１においてレーンＬ_ｊを走行した車両６があると判定した場合は（ステップＳ４２のＹ）、前述の式（５３）により、図２７のステップＳ３４で算出した規格化時間Δｔ_{ｊｓｔｄ１}〜Δｔ_{ｊｓｔｄｌａｓｔ}と図２７のステップＳ３７で算出した基準時間Δｔ_ｊＧＯＣとの差を算出し、規格化時間差Δｔ_{ｊｓｔｄ−ＧＯＣ１}〜Δｔ_{ｊｓｔｄ−ＧＯＣｌａｓｔ}とする（ステップＳ４３）。

次に、計測装置１は、前述の式（５４）により、ステップＳ４３で算出した規格化時間差Δｔ_{ｊｓｔｄ−ＧＯＣ１}〜Δｔ_{ｊｓｔｄ−ＧＯＣｓ}を規格化時間差Δｔ_{ｊｓｔｄ−ＧＯＣ１}〜Δｔ_{ｊｓｔｄ−ＧＯＣｓ}の和で除算し、補正係数Ｒ_ｊ１〜Ｒ_ｊｓを算出する（ステップＳ４４）。

次に、計測装置１は、前述の式（５５）により、ステップＳ４３で算出した規格化時間差Δｔ_{ｊｓｔｄ−ＧＯＣｓ＋１}〜Δｔ_{ｊｓｔｄ−ＧＯＣｌａｓｔ}を規格化時間差Δｔ_{ｊｓｔｄ−ＧＯＣｓ＋１}〜Δｔ_{ｊｓｔｄ−ＧＯＣｌａｓｔ}の和で除算し、補正係数Ｒ_ｊｓ＋１〜Ｒ_{ｊｌａｓｔ}を算出する（ステップＳ４５）。

計測装置１は、図２７のステップＳ３１においてレーンＬ_ｊを走行した車両６がないと判定した場合は（ステップＳ４２のＮ）、ステップＳ４３〜Ｓ４５の処理を行わない。

計測装置１は、整数ｊがＮでない場合は（ステップＳ４６のＮ）、整数ｊに１を加算し（ステップＳ４７）、ステップＳ４２〜Ｓ４５の処理を繰り返し行う。

そして、整数ｊがＮになると（ステップＳ４６のＹ）、計測装置１は補正係数算出ステップの処理を終了する。

１−７．計測装置の構成
図２９は、第１実施形態における計測装置１の構成例を示す図である。図２９に示すように、計測装置１は、制御部１１０と、第１通信部１２０と、記憶部１３０と、第２通信部１４０と、操作部１５０と、を有している。

制御部１１０は、上部構造７に設置された各センサー２１，２２から出力される加速度データに基づいて、車両６が上部構造７を走行した時刻や上部構造７の変位等を算出する。

第１通信部１２０は、各センサー２１，２２から、加速度データを受信する。各センサー２１，２２から出力される加速度データは、例えば、デジタル信号である。第１通信部１２０は、各センサー２１，２２から受信した加速度データを制御部１１０に出力する。

記憶部１３０は、制御部１１０が計算処理や制御処理を行うためのプログラムやデータ等を記憶するメモリーである。また、記憶部１３０は、制御部１１０が所定のアプリケーション機能を実現するためのプログラムやデータ等を記憶している。記憶部１３０は、例えば、ＲＯＭ（Read Only Memory）やフラッシュＲＯＭ、ＲＡＭ（Random Access Memory）等の各種ＩＣ（Integrated Circuit）メモリーやハードディスク、メモリーカードなどの記録媒体等により構成される。

記憶部１３０は、コンピューターにより読み取り可能な装置や媒体である不揮発性の情報記憶装置を含み、各種のプログラムやデータ等は当該情報記憶装置に記憶されていてもよい。情報記憶装置は、光ディスクＤＶＤ、ＣＤ等の光ディスク、ハードディスクドライブ、或いはカード型メモリーやＲＯＭ等の各種のメモリー等であってもよい。また、制御部１１０が通信ネットワーク４を介して各種のプログラムやデータ等を受信して記憶部１３０に記憶させてもよい。

第２通信部１４０は、通信ネットワーク４を介して、制御部１１０の計算結果等の情報をサーバー２に送信する。

操作部１５０は、ユーザーからの操作データを取得し、制御部１１０に送信する処理を行う。

制御部１１０は、第１観測点情報取得部１１１と、第２観測点情報取得部１１２と、第１時間算出部１１３と、補正係数算出部１１４と、たわみ波形算出部１１５と、移動体たわみ波形算出部１１６と、変位算出部１１７と、係数値算出部１１８と、出力処理部１１９と、を備えている。

第１観測点情報取得部１１１は、観測点Ｐ_１〜Ｐ_Ｎを観測するＮ個のセンサー２１による観測情報に基づいて、車両６の複数の車軸が観測点Ｐ_１〜Ｐ_Ｎのいずれかを通過した時刻及び複数の車軸のそれぞれの観測点Ｐ_１〜Ｐ_Ｎのいずれかへの作用に対する応答である第１の物理量としての加速度強度を含む第１観測点情報を取得する処理を行う。すなわち、第１観測点情報取得部１１１は、図２６における第１観測点情報取得ステップの処理を行う。第１観測点情報取得部１１１が取得した第１観測点情報は、記憶部１３０に記憶される。

第２観測点情報取得部１１２は、観測点Ｑ_１〜Ｑ_Ｎを観測するＮ個のセンサー２２による観測情報に基づいて、車両６の複数の車軸が観測点Ｑ_１〜Ｑ_Ｎのいずれかを通過した時刻及び複数の車軸のそれぞれの観測点Ｑ_１〜Ｑ_Ｎのいずれかへの作用に対する応答である第２の物理量としての加速度強度を含む第２観測点情報を取得する処理を行う。すなわち、第２観測点情報取得部１１２は、図２６における第２観測点情報取得ステップの処理を行う。第２観測点情報取得部１１２が取得した第２観測点情報は、記憶部１３０に記憶される。

第１時間算出部１１３は、第１観測点情報取得部１１１が取得した第１観測点情報を用いて、複数の車軸のうちの先頭の車軸が観測点Ｐ_ｊを通過した先頭時刻ｔ_ｊ１から複数の車軸が観測点Ｐ_ｊをそれぞれ通過した時刻ｔ_ｊ１〜ｔ_{ｊｌａｓｔ}までの時間である規格化時間Δｔ_{ｊｓｔｄ１}〜Δｔ_{ｊｓｔｄｌａｓｔ}と、先頭時刻ｔ_ｊ１から複数の車軸による加速度強度ａ_ｐｊ１〜ａ_{ｐｊｌａｓｔ}の総和を所定の配分比_ｉｍａｇＲ_Ｆ：（１−_ｉｍａｇＲ_Ｆ）で配分する時刻である基準時刻ｔ_ｊＧＯＣまでの時間である基準時間Δｔ_ｊＧＯＣとを算出する処理を行う。すなわち、第１時間算出部１１３は、図２６における第１時間算出ステップの処理を行う。第１時間算出部１１３が算出した規格化時間Δｔ_{ｊｓｔｄ１}〜Δｔ_{ｊｓｔｄｌａｓｔ}及び基準時間Δｔ_ｊＧＯＣは、記憶部１３０に記憶される。

補正係数算出部１１４は、第１時間算出部１１３が算出した規格化時間Δｔ_{ｊｓｔｄ１}〜Δｔ_{ｊｓｔｄｌａｓｔ}及び基準時間Δｔ_ｊＧＯＣに基づいて、複数の車軸による加速度強度ａ_ｐｊ１〜ａ_{ｐｊｌａｓｔ}を補正する補正係数Ｒ_ｊ１〜Ｒ_{ｊｌａｓｔ}を算出する処理を行う。すなわち、補正係数算出部１１４は、図２６における補正係数算出ステップの処理を行う。補正係数算出部１１４が算出した補正係数Ｒ_ｊ１〜Ｒ_{ｊｌａｓｔ}は、記憶部１３０に記憶される。

たわみ波形算出部１１５は、第１観測点情報取得部１１１が取得した第１観測点情報及び第２観測点情報取得部１１２が取得した第２観測点情報と、所定の係数ｐと、補正係数算出部１１４が算出した補正係数Ｒ_ｊｋと、上部構造７のたわみの近似式とに基づいて、複数の車軸のそれぞれによる上部構造７のたわみ波形Ｈ_ｊｋ（ｔ）を算出する処理を行う。すなわち、たわみ波形算出部１１５は、図２６におけるたわみ波形算出ステップの処理を行う。たわみ波形算出部１１５が算出したたわみ波形Ｈ_ｊｋ（ｔ）は、記憶部１３０に記憶される。また、所定の係数ｐ及び上部構造７のたわみの近似式は、あらかじめ記憶部１３０に記憶されている。

移動体たわみ波形算出部１１６は、たわみ波形算出部１１５が算出した車両６の複数の車軸のそれぞれによる上部構造７のたわみ波形Ｈ_ｊｋ（ｔ）を加算して、車両６による上部構造７のたわみ波形ＣＰ_ｊｍ（ｔ）を算出する処理を行う。すなわち、移動体たわみ波形算出部１１６は、図２６における移動体たわみ波形算出ステップの処理を行う。移動体たわみ波形算出部１１６が算出したたわみ波形ＣＰ_ｊｍ（ｔ）は、記憶部１３０に記憶される。

変位算出部１１７は、移動体たわみ波形算出部１１６が算出した車両６による上部構造７のたわみ波形ＣＰ_ｊｍ（ｔ）に基づいて、上部構造７の変位ＣＵ_ｅｓｔ（ｔ）を算出する処理を行う。すなわち、変位算出部１１７は、図２６における変位算出ステップの処理を行う。変位算出部１１７が算出した変位ＣＵ_ｅｓｔ（ｔ）は、記憶部１３０に記憶される。

係数値算出部１１８は、複数の試験車両がそれぞれ上部構造７を単独で走行した時に、不図示の変位計が計測した上部構造７の変位の最大値ＣＵ_ｍａｘと移動体たわみ波形算出部１１６が算出したたわみ波形ＣＰ_ｊｍ（ｔ）の最大値ＣＰ_{ｊｍ−ｍａｘ}との相関を近似する近似直線を求め、前述の式（５９）の１次係数ｓ_ｃｕの値及び０次係数ｉ_ｃｕの値を算出する処理を行う。係数値算出部１１８が算出した１次係数ｓ_ｃｕの値及び０次係数ｉ_ｃｕの値は、記憶部１３０に記憶される。

出力処理部１１９は、変位算出部１１７が算出した上部構造７の変位ＣＵ_ｅｓｔ（ｔ）を、第２通信部１４０を介してサーバー２に出力する処理を行う。すなわち、出力処理部１１９は、図２６における出力ステップの処理を行う。

例えば、制御部１１０は、操作部１５０からの操作データに基づいて、未知の車両６が上部構造７を走行した時刻や上部構造７の変位等を算出する第１モードと、１次係数ｓ_ｃｕの値及び０次係数ｉ_ｃｕの値を算出する第２モードとを切り替える。例えば、Ｎ個のセンサー２１及びＮ個のセンサー２２が上部構造７に設置された後、制御部１１０が第２モードに設定された状態で複数の試験車両による荷重試験が行われ、荷重試験の終了後、制御部１１０は第１モードに設定される。

本実施形態では、制御部１１０は、記憶部１３０に記憶された各種のプログラムを実行するプロセッサーであり、記憶部１３０に記憶された計測プログラム１３１を実行することにより、第１観測点情報取得部１１１、第２観測点情報取得部１１２、第１時間算出部１１３、補正係数算出部１１４、たわみ波形算出部１１５、移動体たわみ波形算出部１１６、変位算出部１１７、係数値算出部１１８、出力処理部１１９の各機能を実現する。換言すれば、計測プログラム１３１は、図２６に示したフローチャートの各手順を、コンピューターである計測装置１に実行させるプログラムである。

プロセッサーは、例えば各部の機能が個別のハードウェアで実現されてもよいし、或いは各部の機能が一体のハードウェアで実現されてもよい。例えば、プロセッサーはハードウェアを含み、そのハードウェアは、デジタル信号を処理する回路及びアナログ信号を処理する回路の少なくとも一方を含むことができる。プロセッサーは、ＣＰＵ（Central Processing Unit）、ＧＰＵ（Graphics Processing Unit）、或いはＤＳＰ（Digital Signal Processor）等であってもよい。ただし、制御部１１０は、ＡＳＩＣ（Application Specific Integrated Circuit）などのカスタムＩＣ（Integrated Circuit）として構成され、各部の機能を実現してもよいし、ＣＰＵとＡＳＩＣとによって各部の機能を実現してもよい。

なお、制御部１１０は係数値算出部１１８を含まなくてもよい。例えば、サーバー２あるいは他の装置が、１次係数ｓ_ｃｕの値及び０次係数ｉ_ｃｕの値を算出する処理を行い、これらの値を計測装置１の記憶部１３０に記憶させてもよい。

１−８．作用効果
以上に説明した第１実施形態の計測方法では、計測装置１は、観測点Ｐ_１〜Ｐ_Ｎを観測するＮ個のセンサー２１による観測情報に基づいて、車両６の各車軸について観測点Ｐ_１〜Ｐ_Ｎのいずれかを通過した時刻および加速度強度を含む第１観測点情報を取得する。また、計測装置１は、観測点Ｑ_１〜Ｑ_Ｎを観測するＮ個のセンサー２２による観測情報に基づいて、車両６の各車軸について観測点Ｑ_１〜Ｑ_Ｎのいずれかを通過した時刻および加速度強度を含む第２観測点情報を取得する。また、計測装置１は、第１観測点情報を用いて、式（４５）〜式（５２）により、先頭の車軸が観測点Ｐ_ｊを通過した先頭時刻ｔ_ｊ１から複数の車軸が観測点Ｐ_ｊをそれぞれ通過した時刻ｔ_ｊ１〜ｔ_{ｊｌａｓｔ}までの時間である規格化時間Δｔ_{ｊｓｔｄ１}〜Δｔ_{ｊｓｔｄｌａｓｔ}と、先頭時刻ｔ_ｊ１から複数の車軸による加速度強度ａ_ｐｊ１〜ａ_{ｐｊｌａｓｔ}の総和を所定の配分比_ｉｍａｇＲ_Ｆ：（１−_ｉｍａｇＲ_Ｆ）で配分する基準時刻ｔ_ｊＧＯＣまでの基準時間Δｔ_ｊＧＯＣとを算出する。また、計測装置１は、式（５３）〜式（５５）により、規格化時間Δｔ_{ｊｓｔｄ１}〜Δｔ_{ｊｓｔｄｌａｓｔ}及び基準時間Δｔ_ｊＧＯＣに基づいて、複数の車軸による加速度強度ａ_ｐｊ１〜ａ_{ｐｊｌａｓｔ}を補正する補正係数Ｒ_ｊ１〜Ｒ_{ｊｌａｓｔ}を算出する。そして、計測装置１は、第１観測点情報及び第２観測点情報と、所定の係数ｐと、補正係数Ｒ_ｊ１〜Ｒ_{ｊｌａｓｔ}と、上部構造７のたわみの近似式（３９）とに基づいて、式（５６）により各車軸による上部構造７のたわみ波形Ｈ_ｊｋ（ｔ）を算出し、たわみ波形Ｈ_ｊｋ（ｔ）を加算して車両６による上部構造７のたわみ波形ＣＰ_ｊｍ（ｔ）を算出する。したがって、第１実施形態の計測方法によれば、計測装置１は、構造物である上部構造７の変位を計測せずに、上部構造７を移動する移動体である車両６による上部構造７のたわみ波形を算出することができる。

また、第１実施形態の計測方法によれば、補正係数Ｒ_ｊ１〜Ｒ_{ｊｌａｓｔ}により、各車軸の荷重に対する重心からの距離モーメントの影響が小さくなるように、各車軸による加速度強度ａ_ｐｊ１〜ａ_{ｐｊｌａｓｔ}が補正されるので、車両６の静止時と移動時の荷重の差を考慮して、上部構造７のたわみ波形の算出精度が向上する。

また、第１実施形態の計測方法では、計測装置１は、相関式（５９）により、車両６による上部構造７のたわみ波形ＣＰ_ｊｍ（ｔ）からレーンＬ_ｊの変位ＣＵ_ｅｓｔ（ｔ）を算出することができる。したがって、第１実施形態の計測方法によれば、計測装置１は、上部構造７の変位を計測することなく、上部構造７の変位を推定することができる。

また、第１実施形態の計測方法によれば、上部構造７の変位を計測するための観測点を設定する必要がないので、観測装置による観測点の数が低減され、計測システムが簡易になるため、計測に要するコストが低減される。

また、第１実施形態の計測方法によれば、観測点Ｐ_１〜Ｐ_Ｎ，Ｑ_１〜Ｑ_Ｎは上部構造７の両端部に設定され、上部構造７の中央部等に観測点が設定されないので、計測システム１０の構築やメンテナンスが容易になり、計測に要するコストが低減される。

また、第１実施形態の計測方法によれば、計測装置１は、上部構造７を通過する車両６の車軸重量による上部構造７の変形であるたわみ波形を算出することができるため、上部構造７の損傷を予測するための橋梁５の維持管理のために十分な情報を提供することができる。

２．第２実施形態
第１実施形態の計測方法では、計測装置１は、第２観測点情報を用いずに第１観測点情報を用いて補正係数Ｒ_ｊ１〜Ｒ_{ｊｌａｓｔ}を算出するが、第２実施形態の計測方法では、計測装置１は、第１観測点情報及び第２観測点情報を用いて補正係数Ｒ_ｊ１〜Ｒ_{ｊｌａｓｔ}を算出する。以下、第２実施形態について、第１実施形態と同様の構成要素には同じ符号を付して第１実施形態と重複する説明を省略又は簡略し、主に第１実施形態と異なる内容について説明する。

本実施形態では、車両６の複数の車軸のうちの車軸番号が１である先頭の車軸が観測点Ｐ_ｊを通過した時刻である第１の先頭時刻をｔ_ｊ１＿１とすると、計測装置１は、式（４５）と同様の式（６１）のように、第１の先頭時刻ｔ_ｊ１＿１から車軸番号がｋである車軸が観測点Ｐ_ｊを通過した時刻ｔ_ｊｋ＿１までの時間Δｔ_ｊｋ＿１を算出する。ｋは１以上ｌａｓｔ以下の各整数である。

次に、計測装置１は、第１の先頭時刻ｔ_ｊ１＿１から最後の車軸が観測点Ｐ_ｊを通過した時刻までの時間をΔｔ_{ｊｌａｓｔ＿１}として、式（４６）と同様の式（６２）により、式（６１）で表される時間Δｔ_ｊｋ＿１を時間Δｔ_{ｊｌａｓｔ＿１}で規格化した第１の規格化時間Δｔ_{ｊｓｔｄｋ＿１}を算出する。

次に、計測装置１は、式（４７）と同様の式（６３）により、車軸番号が１である先頭の車軸から車軸番号がｋである車軸までの各車軸による観測点Ｐ_ｊの加速度強度ａ_{ｐｊ１＿１}〜ａ_{ｐｊｋ＿１}の積算値を、車軸番号が１である先頭の車軸から車軸番号がｌａｓｔである最後の車軸までの各車軸による観測点Ｐ_ｊの加速度強度ａ_{ｐｊ１＿１}〜ａ_{ｐｊｌａｓｔ＿１}の総和で規格化した第１の加算規格化衝撃パワー_ａｄｄ｛ａ_{ｐｊ＿ｓｔｄｋ＿１}｝を算出する。

次に、車両６の重心よりも前の荷重と重心よりも後の荷重に配分する所定の配分比をＲ：１−Ｒ＝_ｉｍａｇＲ_Ｆ：（１−_ｉｍａｇＲ_Ｆ）とし、計測装置１は、式（４８）と同様の式（６４）のように、第１の加算規格化衝撃パワー_ａｄｄ｛ａ_{ｐｊ＿ｓｔｄｋ＿１}｝から_ｉｍａｇＲ_Ｆを減算した第１のシフト加算規格化衝撃パワー_{ｓｈｉｆｔ＿ａｄｄ}｛ａ_{ｐｊ＿ｓｔｄｋ＿１}｝を算出する。

次に、計測装置１は、車軸番号が１〜ｓの車軸に対する第１のシフト加算規格化衝撃パワー_{ｓｈｉｆｔ＿ａｄｄ}｛ａ_{ｐｊ＿ｓｔｄｓ＿１}｝が負の値であり、車軸番号がｓ＋１〜ｌａｓｔの車軸に対する第１のシフト加算規格化衝撃パワー_{ｓｈｉｆｔ＿ａｄｄ}｛ａ_{ｐｊ＿ｓｔｄｓ＿１}｝が正の値である車軸番号ｓを求める。

式（４９）と同様、式（６５）のように、車軸番号がｓの車軸に対して、第１の規格化時間Δｔ_{ｊｓｔｄｓ＿１}と負の値である第１のシフト加算規格化衝撃パワー_{ｓｈｉｆｔ＿ａｄｄ}｛ａ_{ｐｊ＿ｓｔｄｓ＿１}｝を示す座標を（ｘ_ｓ＿１，ｙ_ｓ＿１）とする。

同様に、式（５０）と同様、式（６６）のように、車軸番号がｓ＋１の車軸に対して、第１の規格化時間Δｔ_{ｊｓｔｄｓ＋１＿１}と正の値である第１のシフト加算規格化衝撃パワー_{ｓｈｉｆｔ＿ａｄｄ}｛ａ_{ｐｊ＿ｓｔｄｓ＋１＿１}｝を示す座標を（ｘ_{ｓ＋１＿１}，ｙ_{ｓ＋１＿１}）とする。

座標（ｘ_ｓ＿１，ｙ_ｓ＿１）の点と座標（ｘ_{ｓ＋１＿１}，ｙ_{ｓ＋１＿１}）の点を通る直線は式（５１）と同様の式（６７）で与えられる。

計測装置１は、式（５２）と同様の式（６８）により、式（６７）で表される直線とｙ＝０との交点のｘ座標を算出し、算出したｘ座標を第１の基準時間Δｔ_{ｊＧＯＣ＿１}とする。

ここで、車軸番号が１である先頭の車軸から車軸番号がｌａｓｔである最後の車軸までの各車軸による観測点Ｐ_ｊの加速度強度ａ_{ｐｊ１＿１}〜ａ_{ｐｊｌａｓｔ＿１}の総和を所定の配分比_ｉｍａｇＲ_Ｆ：（１−_ｉｍａｇＲ_Ｆ）で配分する時刻を第１の基準時刻ｔ_{ｊＧＯＣ＿１}とすると、第１の基準時間Δｔ_{ｊＧＯＣ＿１}は、先頭の車軸が観測点Ｐ_ｊを通過した時刻から第１の基準時刻ｔ_{ｊＧＯＣ＿１}までの時間である。第１の基準時刻ｔ_{ｊＧＯＣ＿１}は、車両６の重心が観測点Ｐ_ｊを通過したと推定される時刻である。

さらに、車両６の複数の車軸のうちの車軸番号が１である先頭の車軸が観測点Ｑ_ｊを通過した時刻である第２の先頭時刻をｔ_ｊ１＿２とすると、計測装置１は、式（６１）と同様の式（６９）のように、第２の先頭時刻ｔ_ｊ１＿２から車軸番号がｋである車軸が観測点Ｑ_ｊを通過した時刻ｔ_ｊｋ＿２までの時間Δｔ_ｊｋ＿２を算出する。

次に、計測装置１は、第２の先頭時刻ｔ_ｊ１＿２から最後の車軸が観測点Ｑ_ｊを通過した時刻までの時間をΔｔ_{ｊｌａｓｔ＿２}として、式（６２）と同様の式（７０）により、式（６９）で表される時間Δｔ_ｊｋ＿２を時間Δｔ_{ｊｌａｓｔ＿２}で規格化した第２の規格化時間Δｔ_{ｊｓｔｄｋ＿２}を算出する。

次に、計測装置１は、式（６３）と同様の式（７１）により、車軸番号が１である先頭の車軸から車軸番号がｋである車軸までの各車軸による観測点Ｑ_ｊの加速度強度ａ_{ｐｊ１＿２}〜ａ_{ｐｊｋ＿２}の積算値を、車軸番号が１である先頭の車軸から車軸番号がｌａｓｔである最後の車軸までの各車軸による観測点Ｑ_ｊの加速度強度ａ_{ｐｊ１＿２}〜ａ_{ｐｊｌａｓｔ＿２}の総和で規格化した第２の加算規格化衝撃パワー_ａｄｄ｛ａ_{ｐｊ＿ｓｔｄｋ＿２}｝を算出する。

次に、計測装置１は、式（６４）と同様の式（７２）のように、第２の加算規格化衝撃パワー_ａｄｄ｛ａ_{ｐｊ＿ｓｔｄｋ＿２}｝から_ｉｍａｇＲ_Ｆを減算した第２のシフト加算規格化衝撃パワー_{ｓｈｉｆｔ＿ａｄｄ}｛ａ_{ｐｊ＿ｓｔｄｋ＿２}｝を算出する。

次に、計測装置１は、車軸番号が１〜ｓの車軸に対する第２のシフト加算規格化衝撃パワー_{ｓｈｉｆｔ＿ａｄｄ}｛ａ_{ｐｊ＿ｓｔｄｓ＿２}｝が負の値であり、車軸番号がｓ＋１〜ｌａｓｔの車軸に対する第２のシフト加算規格化衝撃パワー_{ｓｈｉｆｔ＿ａｄｄ}｛ａ_{ｐｊ＿ｓｔｄｓ＿２}｝が正の値である車軸番号ｓを求める。

式（６５）と同様、式（７３）のように、車軸番号がｓの車軸に対して、第２の規格化時間Δｔ_{ｊｓｔｄｓ＿２}と負の値である第２のシフト加算規格化衝撃パワー_{ｓｈｉｆｔ＿ａｄｄ}｛ａ_{ｐｊ＿ｓｔｄｓ＿２}｝を示す座標を（ｘ_ｓ＿２，ｙ_ｓ＿２）とする。

同様に、式（６６）と同様、式（７４）のように、車軸番号がｓ＋１の車軸に対して、第２の規格化時間Δｔ_{ｊｓｔｄｓ＋１＿２}と正の値である第２のシフト加算規格化衝撃パワー_{ｓｈｉｆｔ＿ａｄｄ}｛ａ_{ｐｊ＿ｓｔｄｓ＋１＿２}｝を示す座標を（ｘ_{ｓ＋１＿２}，ｙ_{ｓ＋１＿２}）とする。

座標（ｘ_ｓ＿２，ｙ_ｓ＿２）の点と座標（ｘ_{ｓ＋１＿２}，ｙ_{ｓ＋１＿２}）の点を通る直線は式（６７）と同様の式（７５）で与えられる。

計測装置１は、式（６８）と同様の式（７６）により、式（７５）で表される直線とｙ＝０との交点のｘ座標を算出し、算出したｘ座標を第２の基準時間Δｔ_{ｊＧＯＣ＿２}とする。

ここで、車軸番号が１である先頭の車軸から車軸番号がｌａｓｔである最後の車軸までの各車軸による観測点Ｑ_ｊの加速度強度ａ_{ｐｊ１＿２}〜ａ_{ｐｊｌａｓｔ＿２}の総和を所定の配分比_ｉｍａｇＲ_Ｆ：（１−_ｉｍａｇＲ_Ｆ）で配分する時刻を第２の基準時刻ｔ_{ｊＧＯＣ＿２}とすると、第２の基準時間Δｔ_{ｊＧＯＣ＿２}は、先頭の車軸が観測点Ｑ_ｊを通過した時刻から第２の基準時刻ｔ_{ｊＧＯＣ＿２}までの時間である。第２の基準時刻ｔ_{ｊＧＯＣ＿２}は、車両６の重心が観測点Ｑ_ｊを通過したと推定される時刻である。

次に、計測装置１は、式（７７）により、各車軸の第１の規格化時間Δｔ_{ｊｓｔｄｋ＿１}と第２の規格化時間Δｔ_{ｊｓｔｄｋ＿２}の平均を算出し、規格化時間Δｔ_{ｊｓｔｄｋ}とする。

また、計測装置１は、式（７８）により、第１の基準時間Δｔ_{ｊＧＯＣ＿１}と第２の基準時間Δｔ_{ｊＧＯＣ＿２}の平均を算出し、基準時間Δｔ_ｊＧＯＣとする。

そして、計測装置１は、前述の式（５３）によって規格化時間差Δｔ_{ｊｓｔｄ−ＧＯＣｋ}を算出し、前述の式（５４）及び式（５５）によって車軸番号が１〜ｌａｓｔの車軸の補正係数Ｒ_ｊ１〜Ｒ_{ｊｌａｓｔ}を算出する。

図３０は、第２実施形態の計測方法の手順の一例を示すフローチャート図である。本実施形態では、計測装置１が図３０に示す手順を実行する。

図３０に示すように、まず、計測装置１は、図２６のステップＳ１と同様、第１観測点情報取得ステップの処理を行う（ステップＳ１０１）。

次に、計測装置１は、図２６のステップＳ２と同様、第２観測点情報取得ステップの処理を行う（ステップＳ１０２）。

次に、計測装置１は、ステップＳ１０１で取得した第１観測点情報を用いて、複数の車軸のうちの先頭の車軸が観測点Ｐ_ｊを通過した第１の先頭時刻ｔ_ｊ１＿１から複数の車軸が観測点Ｐ_ｊをそれぞれ通過した時刻ｔ_ｊ１＿１〜ｔ_{ｊｌａｓｔ＿１}までの時間である第１の規格化時間Δｔ_{ｊｓｔｄ１＿１}〜Δｔ_{ｊｓｔｄｌａｓｔ＿１}と、第１の先頭時刻ｔ_ｊ１＿１から複数の車軸による加速度強度ａ_{ｐｊ１＿１}〜ａ_{ｐｊｌａｓｔ＿１}の総和を所定の配分比_ｉｍａｇＲ_Ｆ：（１−_ｉｍａｇＲ_Ｆ）で配分する時刻である第１の基準時刻ｔ_{ｊＧＯＣ＿１}までの時間である第１の基準時間Δｔ_{ｊＧＯＣ＿１}とを算出する（ステップＳ１０３）。具体的には、計測装置１は、前述の式（６１）〜式（６８）により、第１の規格化時間Δｔ_{ｊｓｔｄ１＿１}〜Δｔ_{ｊｓｔｄｌａｓｔ＿１}と第１の基準時間Δｔ_{ｊＧＯＣ＿１}とを算出する。このステップＳ１０３は、第１時間算出ステップである。

次に、計測装置１は、ステップＳ１０２で取得した第２観測点情報を用いて、複数の車軸のうちの先頭の車軸が観測点Ｑ_ｊを通過した第２の先頭時刻ｔ_ｊ１＿２から複数の車軸が観測点Ｑ_ｊをそれぞれ通過した時刻ｔ_ｊ１＿２〜ｔ_{ｊｌａｓｔ＿２}までの時間である第２の規格化時間Δｔ_{ｊｓｔｄ１＿２}〜Δｔ_{ｊｓｔｄｌａｓｔ＿２}と、第２の先頭時刻ｔ_ｊ１＿２から複数の車軸による加速度強度ａ_{ｐｊ１＿２}〜ａ_{ｐｊｌａｓｔ＿２}の総和を所定の配分比_ｉｍａｇＲ_Ｆ：（１−_ｉｍａｇＲ_Ｆ）で配分する時刻である第２の基準時刻ｔ_{ｊＧＯＣ＿２}までの時間である第２の基準時間Δｔ_{ｊＧＯＣ＿２}とを算出する（ステップＳ１０４）。具体的には、計測装置１は、前述の式（６９）〜式（７６）により、第２の規格化時間Δｔ_{ｊｓｔｄ１＿２}〜Δｔ_{ｊｓｔｄｌａｓｔ＿２}と第２の基準時間Δｔ_{ｊＧＯＣ＿２}とを算出する。このステップＳ１０４は、第２時間算出ステップである。

次に、計測装置１は、ステップＳ１０３で算出した第１の規格化時間Δｔ_{ｊｓｔｄ１＿１}〜Δｔ_{ｊｓｔｄｌａｓｔ＿１}及び第１の基準時間Δｔ_{ｊＧＯＣ＿１}と、ステップＳ１０４で算出した第２の規格化時間Δｔ_{ｊｓｔｄ１＿２}〜Δｔ_{ｊｓｔｄｌａｓｔ＿２}及び第２の基準時間Δｔ_{ｊＧＯＣ＿２}とに基づいて、複数の車軸による観測点Ｐ_ｊの加速度強度ａ_{ｐｊ１＿１}〜ａ_{ｐｊｌａｓｔ＿１}及び観測点Ｑ_ｊのａ_{ｐｊ１＿２}〜ａ_{ｐｊｌａｓｔ＿２}を補正する補正係数Ｒ_ｊ１〜Ｒ_{ｊｌａｓｔ}を算出する（ステップＳ１０５）。具体的には、計測装置１は、前述の式（７７）及び式（７８）により、規格化時間Δｔ_{ｊｓｔｄｋ}及び基準時間Δｔ_ｊＧＯＣを算出し、式（５３）〜式（５５）により、補正係数Ｒ_ｊ１〜Ｒ_{ｊｌａｓｔ}を算出する。このステップＳ１０５は、補正係数算出ステップである。

次に、計測装置１は、図２６のステップＳ５と同様、たわみ波形算出ステップの処理を行う（ステップＳ１０６）。

次に、計測装置１は、図２６のステップＳ６と同様、移動体たわみ波形算出ステップの処理を行う（ステップＳ１０７）。

次に、計測装置１は、図２６のステップＳ７と同様、変位算出ステップの処理を行う（ステップＳ１０８）。

次に、計測装置１は、図２６のステップＳ８と同様、出力ステップの処理を行う（ステップＳ１０９）。

計測装置１は、計測を終了するまで（ステップＳ１１０のＮ）、ステップＳ１０１〜Ｓ１０９の処理を繰り返し行う。

図３１は、図３０のステップＳ３である第１時間算出ステップの手順の一例を示すフローチャート図である。

図３１に示すように、まず、計測装置１は、整数ｊを１に設定し（ステップＳ１３０）、第１観測点情報を用いて、レーンＬ_ｊを走行した車両６の有無を判定する（ステップＳ１３１）。

そして、計測装置１は、ステップＳ１３１においてレーンＬ_ｊを走行した車両６があると判定した場合は（ステップＳ１３２のＹ）、第１観測点情報を用いて、前述の式（６１）により、時間Δｔ_ｊ１＿１〜Δｔ_{ｊｌａｓｔ＿１}を算出する（ステップＳ１３３）。すなわち、計測装置１は、車両６の先頭の車軸が観測点Ｐ_ｊを通過した第１の先頭時刻ｔ_ｊ１＿１から複数の車軸が観測点Ｐ_ｊをそれぞれ通過した時刻ｔ_ｊ１＿１〜ｔ_{ｊｌａｓｔ＿１}までの時間Δｔ_ｊ１＿１〜Δｔ_{ｊｌａｓｔ＿１}を算出する。

次に、計測装置１は、前述の式（６２）により、ステップＳ１３３で算出した時間Δｔ_ｊ１＿１〜Δｔ_{ｊｌａｓｔ＿１}を時間Δｔ_{ｊｌａｓｔ＿１}で除算して第１の規格化時間Δｔ_{ｊｓｔｄ１＿１}〜Δｔ_{ｊｓｔｄｌａｓｔ＿１}を算出する（ステップＳ１３４）。

次に、計測装置１は、前述の式（６３）により、第１観測点情報の加速度強度ａ_{ｐｊ１＿１}〜ａ_{ｐｊｋ＿１}の積算値を、加速度強度ａ_{ｐｊ１＿１}〜ａ_{ｐｊｌａｓｔ＿１}の総和で除算し、第１の加算規格化衝撃パワー_ａｄｄ｛ａ_{ｐｊ＿ｓｔｄ１＿１}｝〜_ａｄｄ｛ａ_{ｐｊ＿ｓｔｄｌａｓｔ＿１}｝を算出する（ステップＳ１３５）。

次に、計測装置１は、前述の式（６４）により、ステップＳ１３５で算出した第１の加算規格化衝撃パワー_ａｄｄ｛ａ_{ｐｊ＿ｓｔｄ１＿１}｝〜_ａｄｄ｛ａ_{ｐｊ＿ｓｔｄｌａｓｔ＿１}｝から_ｉｍａｇＲ_Ｆを減算して第１のシフト加算規格化衝撃パワー_{ｓｈｉｆｔ＿ａｄｄ}｛ａ_{ｐｊ＿ｓｔｄ１＿１}｝〜_{ｓｈｉｆｔ＿ａｄｄ}｛ａ_{ｐｊ＿ｓｔｄｌａｓｔ＿１}｝を算出する（ステップＳ１３６）。

次に、計測装置１は、前述の式（６５）〜式（６８）により、座標（Δｔ_{ｊｓｔｄｓ＿１}，_{ｓｈｉｆｔ＿ａｄｄ}｛ａ_{ｐｊ＿ｓｔｄｓ＿１}｝（＜０））と座標（Δｔ_{ｊｓｔｄｓ＋１＿１}，_{ｓｈｉｆｔ＿ａｄｄ}｛ａ_{ｐｊ＿ｓｔｄｓ＋１＿１}｝（＞０））とを通直線とｙ＝０との交点のｘ座標を第１の基準時間Δｔ_{ｊＧＯＣ＿１}とする（ステップＳ１３７）。

計測装置１は、ステップＳ１３１においてレーンＬ_ｊを走行した車両６がないと判定した場合は（ステップＳ１３２のＮ）、ステップＳ１３３〜Ｓ１３７の処理を行わない。

計測装置１は、整数ｊがＮでない場合は（ステップＳ１３８のＮ）、整数ｊに１を加算し（ステップＳ１３９）、ステップＳ１３１〜Ｓ１３７の処理を繰り返し行う。

そして、整数ｊがＮになると（ステップＳ１３８のＹ）、計測装置１は第１時間算出ステップの処理を終了する。

図３２は、図３０のステップＳ４である第２時間算出ステップの手順の一例を示すフローチャート図である。

図３２に示すように、まず、計測装置１は、整数ｊを１に設定し（ステップＳ１４１）、図３１のステップＳ１３１においてレーンＬ_ｊを走行した車両６があると判定した場合は（ステップＳ１４２のＹ）、第２観測点情報を用いて、前述の式（６９）により、時間Δｔ_ｊ１＿２〜Δｔ_{ｊｌａｓｔ＿２}を算出する（ステップＳ１４３）。すなわち、計測装置１は、車両６の先頭の車軸が観測点Ｑ_ｊを通過した第２の先頭時刻ｔ_ｊ１＿２から複数の車軸が観測点Ｑ_ｊをそれぞれ通過した時刻ｔ_ｊ１＿２〜ｔ_{ｊｌａｓｔ＿２}までの時間Δｔ_ｊ１＿２〜Δｔ_{ｊｌａｓｔ＿２}を算出する。

次に、計測装置１は、前述の式（７０）により、ステップＳ１４３で算出した時間Δｔ_ｊ１＿２〜Δｔ_{ｊｌａｓｔ＿２}を時間Δｔ_{ｊｌａｓｔ＿２}で除算して第２の規格化時間Δｔ_{ｊｓｔｄ１＿２}〜Δｔ_{ｊｓｔｄｌａｓｔ＿２}を算出する（ステップＳ１４４）。

次に、計測装置１は、前述の式（７１）により、第２観測点情報の加速度強度ａ_{ｐｊ１＿２}〜ａ_{ｐｊｋ＿２}の積算値を、加速度強度ａ_{ｐｊ１＿２}〜ａ_{ｐｊｌａｓｔ＿２}の総和で除算し、第２の加算規格化衝撃パワー_ａｄｄ｛ａ_{ｐｊ＿ｓｔｄ１＿２}｝〜_ａｄｄ｛ａ_{ｐｊ＿ｓｔｄｌａｓｔ＿２}｝を算出する（ステップＳ１４５）。

次に、計測装置１は、前述の式（７２）により、ステップＳ１４５で算出した第２の加算規格化衝撃パワー_ａｄｄ｛ａ_{ｐｊ＿ｓｔｄ１＿２}｝〜_ａｄｄ｛ａ_{ｐｊ＿ｓｔｄｌａｓｔ＿２}｝から_ｉｍａｇＲ_Ｆを減算して第２のシフト加算規格化衝撃パワー_{ｓｈｉｆｔ＿ａｄｄ}｛ａ_{ｐｊ＿ｓｔｄ１＿２}｝〜_{ｓｈｉｆｔ＿ａｄｄ}｛ａ_{ｐｊ＿ｓｔｄｌａｓｔ＿２}｝を算出する（ステップＳ１４６）。

次に、計測装置１は、前述の式（７３）〜式（７６）により、座標（Δｔ_{ｊｓｔｄｓ＿２}，_{ｓｈｉｆｔ＿ａｄｄ}｛ａ_{ｐｊ＿ｓｔｄｓ＿２}｝（＜０））と座標（Δｔ_{ｊｓｔｄｓ＋１＿２}，_{ｓｈｉｆｔ＿ａｄｄ}｛ａ_{ｐｊ＿ｓｔｄｓ＋１＿２}｝（＞０））とを通直線とｙ＝０との交点のｘ座標を第２の基準時間Δｔ_{ｊＧＯＣ＿２}とする（ステップＳ１４７）。

計測装置１は、図３１のステップＳ１３１においてレーンＬ_ｊを走行した車両６がないと判定した場合は（ステップＳ１４２のＮ）、ステップＳ１４３〜Ｓ１４７の処理を行わない。

計測装置１は、整数ｊがＮでない場合は（ステップＳ１４８のＮ）、整数ｊに１を加算し（ステップＳ１４９）、ステップＳ１４１〜Ｓ１４７の処理を繰り返し行う。

そして、整数ｊがＮになると（ステップＳ１４８のＹ）、計測装置１は第２時間算出ステップの処理を終了する。

図３３は、図３０のステップＳ５である補正係数算出ステップの手順の一例を示すフローチャート図である。

図３３に示すように、まず、計測装置１は、整数ｊを１に設定し（ステップＳ１５１）、図３１のステップＳ１３１においてレーンＬ_ｊを走行した車両６があると判定した場合は（ステップＳ１５２のＹ）、前述の式（７７）により、図３１のステップＳ１３４で算出した第１の規格化時間Δｔ_{ｊｓｔｄ１＿１}〜Δｔ_{ｊｓｔｄｌａｓｔ＿１}と図３２のステップＳ１４４で算出した第２の規格化時間Δｔ_{ｊｓｔｄ１＿２}〜Δｔ_{ｊｓｔｄｌａｓｔ＿２}との平均値をそれぞれ算出し、規格化時間Δｔ_{ｊｓｔｄ１}〜Δｔ_{ｊｓｔｄｌａｓｔ}とする（ステップＳ１５３）。

次に、計測装置１は、前述の式（７８）により、図３１のステップＳ１３７で算出した第１の基準時間Δｔ_{ｊＧＯＣ＿１}と図３２のステップＳ１４７で算出した第２の基準時間Δｔ_{ｊＧＯＣ＿２}との平均値を算出し、基準時間Δｔ_ｊＧＯＣとする（ステップＳ１５４）。

次に、計測装置１は、前述の式（５３）により、ステップＳ１５３で算出した規格化時間Δｔ_{ｊｓｔｄ１}〜Δｔ_{ｊｓｔｄｌａｓｔ}とステップＳ１５４で算出した基準時間Δｔ_ｊＧＯＣとの差を算出し、規格化時間差Δｔ_{ｊｓｔｄ−ＧＯＣ１}〜Δｔ_{ｊｓｔｄ−ＧＯＣｌａｓｔ}とする（ステップＳ１５５）。

次に、計測装置１は、前述の式（５４）により、ステップＳ１５５で算出した規格化時間差Δｔ_{ｊｓｔｄ−ＧＯＣ１}〜Δｔ_{ｊｓｔｄ−ＧＯＣｓ}を規格化時間差Δｔ_{ｊｓｔｄ−ＧＯＣ１}〜Δｔ_{ｊｓｔｄ−ＧＯＣｓ}の和で除算し、補正係数Ｒ_ｊ１〜Ｒ_ｊｓを算出する（ステップＳ１５６）。

次に、計測装置１は、前述の式（５５）により、ステップＳ１５５で算出した規格化時間差Δｔ_{ｊｓｔｄ−ＧＯＣｓ＋１}〜Δｔ_{ｊｓｔｄ−ＧＯＣｌａｓｔ}を規格化時間差Δｔ_{ｊｓｔｄ−ＧＯＣｓ＋１}〜Δｔ_{ｊｓｔｄ−ＧＯＣｌａｓｔ}の和で除算し、補正係数Ｒ_ｊｓ＋１〜Ｒ_{ｊｌａｓｔ}を算出する（ステップＳ１５７）。

計測装置１は、図３１のステップＳ１３１においてレーンＬ_ｊを走行した車両６がないと判定した場合は（ステップＳ１５２のＮ）、ステップＳ１５３〜Ｓ１５７の処理を行わない。

計測装置１は、整数ｊがＮでない場合は（ステップＳ１５８のＮ）、整数ｊに１を加算し（ステップＳ１５９）、ステップＳ１５２〜Ｓ１５７の処理を繰り返し行う。

そして、整数ｊがＮになると（ステップＳ１５８のＹ）、計測装置１は補正係数算出ステップの処理を終了する。

図３４は、第２実施形態における計測装置１の構成例を示す図である。図３４に示すように、第１実施形態と同様、計測装置１は、制御部１１０と、第１通信部１２０と、記憶部１３０と、第２通信部１４０と、操作部１５０と、を有している。

第１通信部１２０、記憶部１３０、第２通信部１４０及び操作部１５０がそれぞれ行う処理は、第１実施形態と同様であるため、その説明を省略する。

制御部１１０は、上部構造７に設置された各センサー２１，２２から出力される加速度データに基づいて、車両６が上部構造７を走行した時刻や車両６による荷重等を算出する。

制御部１１０は、第１観測点情報取得部１１１と、第２観測点情報取得部１１２と、第１時間算出部１１３と、補正係数算出部１１４と、たわみ波形算出部１１５と、移動体たわみ波形算出部１１６と、変位算出部１１７と、係数値算出部１１８と、出力処理部１１９と、第２時間算出部１６１と、を備えている。

第１観測点情報取得部１１１、第２観測点情報取得部１１２、たわみ波形算出部１１５、移動体たわみ波形算出部１１６、変位算出部１１７、係数値算出部１１８及び出力処理部１１９がそれぞれ行う処理は、第１実施形態と同様であるため、その説明を省略する。

第１時間算出部１１３は、第１観測点情報取得部１１１が取得した第１観測点情報を用いて、複数の車軸のうちの先頭の車軸が観測点Ｐ_ｊを通過した第１の先頭時刻ｔ_ｊ１＿１から複数の車軸が観測点Ｐ_ｊをそれぞれ通過した時刻ｔ_ｊ１＿１〜ｔ_{ｊｌａｓｔ＿１}までの時間である第１の規格化時間Δｔ_{ｊｓｔｄ１＿１}〜Δｔ_{ｊｓｔｄｌａｓｔ＿１}と、第１の先頭時刻ｔ_ｊ１＿１から複数の車軸による観測点Ｐ_ｊの加速度強度ａ_{ｐｊ１＿１}〜ａ_{ｐｊｌａｓｔ＿１}の総和を所定の配分比_ｉｍａｇＲ_Ｆ：（１−_ｉｍａｇＲ_Ｆ）で配分する時刻である第１の基準時刻ｔ_{ｊＧＯＣ＿１}までの時間である第１の基準時間Δｔ_{ｊＧＯＣ＿１}とを算出する処理を行う。すなわち、第１時間算出部１１３は、図３０における第１時間算出ステップの処理を行う。第１時間算出部１１３が算出した第１の規格化時間Δｔ_{ｊｓｔｄ１＿１}〜Δｔ_{ｊｓｔｄｌａｓｔ＿１}及び第１の基準時間Δｔ_{ｊＧＯＣ＿１}は、記憶部１３０に記憶される。

第２時間算出部１６１は、第２観測点情報取得部１１２が取得した第２観測点情報を用いて、複数の車軸のうちの先頭の車軸が観測点Ｑ_ｊを通過した第２の先頭時刻ｔ_ｊ１＿２から複数の車軸が観測点Ｑ_ｊをそれぞれ通過した時刻ｔ_ｊ１＿２〜ｔ_{ｊｌａｓｔ＿２}までの時間である第２の規格化時間Δｔ_{ｊｓｔｄ１＿２}〜Δｔ_{ｊｓｔｄｌａｓｔ＿２}と、第２の先頭時刻ｔ_ｊ１＿２から複数の車軸による観測点Ｑ_ｊの加速度強度ａ_{ｐｊ１＿２}〜ａ_{ｐｊｌａｓｔ＿２}の総和を所定の配分比_ｉｍａｇＲ_Ｆ：（１−_ｉｍａｇＲ_Ｆ）で配分する時刻である第２の基準時刻ｔ_{ｊＧＯＣ＿２}までの時間である第２の基準時間Δｔ_{ｊＧＯＣ＿２}とを算出する処理を行う。すなわち、第２時間算出部１６１は、図３０における第２時間算出ステップの処理を行う。第２時間算出部１６１が算出した第２の規格化時間Δｔ_{ｊｓｔｄ１＿２}〜Δｔ_{ｊｓｔｄｌａｓｔ＿２}及び第２の基準時間Δｔ_{ｊＧＯＣ＿２}は、記憶部１３０に記憶される。

補正係数算出部１１４は、第１時間算出部１１３が算出した第１の規格化時間Δｔ_{ｊｓｔｄ１＿１}〜Δｔ_{ｊｓｔｄｌａｓｔ＿１}及び第１の基準時間Δｔ_{ｊＧＯＣ＿１}と、第２時間算出部１６１が算出した第２の規格化時間Δｔ_{ｊｓｔｄ１＿２}〜Δｔ_{ｊｓｔｄｌａｓｔ＿２}及び第２の基準時間Δｔ_{ｊＧＯＣ＿２}とに基づいて、複数の車軸による観測点Ｐ_ｊの加速度強度ａ_{ｐｊ１＿１}〜ａ_{ｐｊｌａｓｔ＿１}及び観測点Ｑ_ｊの加速度強度ａ_{ｐｊ１＿２}〜ａ_{ｐｊｌａｓｔ＿２}を補正する補正係数Ｒ_ｊ１〜Ｒ_{ｊｌａｓｔ}を算出する処理を行う。すなわち、補正係数算出部１１４は、図３０における補正係数算出ステップの処理を行う。補正係数算出部１１４が算出した補正係数Ｒ_ｊ１〜Ｒ_{ｊｌａｓｔ}は、記憶部１３０に記憶される。

例えば、制御部１１０は、操作部１５０からの操作データに基づいて、未知の車両６が上部構造７を走行した時刻や車両６による荷重等を算出する第１モードと、１次係数ｓ_ｃｐの値及び０次係数ｉ_ｃｐの値を算出する第２モードとを切り替える。例えば、Ｎ個のセンサー２１及びＮ個のセンサー２２が上部構造７に設置された後、制御部１１０が第２モードに設定された状態で複数の試験車両による荷重試験が行われ、荷重試験の終了後、制御部１１０は第１モードに設定される。

第１実施形態と同様、制御部１１０は、記憶部１３０に記憶された各種のプログラムを実行するプロセッサーであり、記憶部１３０に記憶された計測プログラム１３１を実行することにより、第１観測点情報取得部１１１、第２観測点情報取得部１１２、第１時間算出部１１３、補正係数算出部１１４、たわみ波形算出部１１５、移動体たわみ波形算出部１１６、変位算出部１１７、係数値算出部１１８、出力処理部１１９、第２時間算出部１６１の各機能を実現する。換言すれば、計測プログラム１３１は、図３０に示したフローチャートの各手順を、コンピューターである計測装置１に実行させるプログラムである。ただし、制御部１１０は、ＡＳＩＣなどのカスタムＩＣとして構成され、各部の機能を実現してもよいし、ＣＰＵとＡＳＩＣとによって各部の機能を実現してもよい。

なお、制御部１１０は係数値算出部１１８を含まなくてもよい。例えば、サーバー２あるいは他の装置が、１次係数ｓ_ｃｕの値及び０次係数ｉ_ｃｕの値を算出する処理を行い、これらの値を計測装置１の記憶部１３０に記憶させてもよい。

以上に説明した第２実施形態の計測方法では、計測装置１は、第１観測点情報を用いて、式（６１）〜式（６８）により、第１の規格化時間Δｔ_{ｊｓｔｄ１＿１}〜Δｔ_{ｊｓｔｄｌａｓｔ＿１}及び第１の基準時間Δｔ_{ｊＧＯＣ＿１}を算出する。また、計測装置１は、第２観測点情報を用いて、式（６９）〜式（７６）により、第２の規格化時間Δｔ_{ｊｓｔｄ１＿２}〜Δｔ_{ｊｓｔｄｌａｓｔ＿２}及び第２の基準時間Δｔ_{ｊＧＯＣ＿２}を算出する。また、計測装置１は、式（７７）、式（７８）、式（５３）〜式（５５）により、第１の規格化時間Δｔ_{ｊｓｔｄ１＿１}〜Δｔ_{ｊｓｔｄｌａｓｔ＿１}及び第１の基準時間Δｔ_{ｊＧＯＣ＿１}と、第２の規格化時間Δｔ_{ｊｓｔｄ１＿２}〜Δｔ_{ｊｓｔｄｌａｓｔ＿２}及び第２の基準時間Δｔ_{ｊＧＯＣ＿２}とに基づいて、補正係数Ｒ_ｊ１〜Ｒ_{ｊｌａｓｔ}を算出する。そして、計測装置１は、第１観測点情報及び第２観測点情報と、所定の係数ｐと、補正係数Ｒ_ｊ１〜Ｒ_{ｊｌａｓｔ}と、上部構造７のたわみの近似式（３９）とに基づいて、式（５６）により各車軸による上部構造７のたわみ波形Ｈ_ｊｋ（ｔ）を算出し、たわみ波形Ｈ_ｊｋ（ｔ）を加算して車両６による上部構造７のたわみ波形ＣＰ_ｊｍ（ｔ）を算出する。したがって、第２実施形態の計測方法によれば、計測装置１は、構造物である上部構造７の変位を計測せずに、上部構造７を移動する移動体である車両６による上部構造７のたわみ波形を算出することができる。

また、第２実施形態の計測方法によれば、補正係数Ｒ_ｊ１〜Ｒ_{ｊｌａｓｔ}により、各車軸の荷重に対する重心からの距離モーメントの影響が小さくなるように、各車軸による加速度強度ａ_ｐｊ１〜ａ_{ｐｊｌａｓｔ}が補正されるので、車両６の静止時と移動時の荷重の差を考慮して、上部構造７のたわみ波形の算出精度が向上する。

さらに、第２実施形態の計測方法では、計測装置１は、式（７７）により、第１の規格化時間Δｔ_{ｊｓｔｄ１＿１}〜Δｔ_{ｊｓｔｄｌａｓｔ＿１}と第２の規格化時間Δｔ_{ｊｓｔｄ１＿２}〜Δｔ_{ｊｓｔｄｌａｓｔ＿２}の平均値を算出して規格化時間Δｔ_{ｊｓｔｄｋ}とし、式（７８）により、第１の基準時間Δｔ_{ｊＧＯＣ＿１}と第２の基準時間Δｔ_{ｊＧＯＣ＿２}の平均値を算出して基準時間Δｔ_ｊＧＯＣとする。そして、計測装置１は、規格化時間Δｔ_{ｊｓｔｄｋ}と基準時間Δｔ_ｊＧＯＣとを用いて、式（５３）〜式（５５）により、補正係数Ｒ_ｊ１〜Ｒ_{ｊｌａｓｔ}を算出する。したがって、第２実施形態の計測方法によれば、平均化によって規格化時間Δｔ_{ｊｓｔｄｋ}や基準時間Δｔ_ｊＧＯＣに含まれるランダムノイズ成分が低減されるので、補正係数Ｒ_ｊ１〜Ｒ_{ｊｌａｓｔ}の精度が向上し、上部構造７のたわみ波形の算出精度がさらに向上する。

また、第２実施形態の計測方法によれば、第１実施形態の計測方法と同様、計測装置１は、上部構造７の変位を計測することなく、上部構造７の変位を推定することができ、計測に要するコストが低減される。

３．変形例
本発明は本実施形態に限定されず、本発明の要旨の範囲内で種々の変形実施が可能である。

上記の第１実施形態では、計測装置１は、車両６の各車軸が観測点Ｐ_ｊを通過した時刻及び各車軸の観測点Ｐ_ｊに対する加速度強度を含む観測点情報を第１観測点情報として取得し、車両６の各車軸が観測点Ｑ_ｊを通過した時刻及び各車軸の観測点Ｑ_ｊに対する加速度強度を含む観測点情報を第２観測点情報として取得する。これに対して、計測装置１は、車両６の各車軸が観測点Ｑ_ｊを通過した時刻及び各車軸の観測点Ｑ_ｊに対する加速度強度を含む観測点情報を第１観測点情報として取得し、車両６の各車軸が観測点Ｐ_ｊを通過した時刻及び各車軸の観測点Ｐ_ｊに対する加速度強度を含む観測点情報を第２観測点情報として取得してもよい。そして、計測装置１は、この第１観測点情報を用いて、規格化時間Δｔ_{ｊｓｔｄ１}〜Δｔ_{ｊｓｔｄｌａｓｔ}及び基準時間Δｔ_ｊＧＯＣを算出し、規格化時間Δｔ_{ｊｓｔｄ１}〜Δｔ_{ｊｓｔｄｌａｓｔ}及び基準時間Δｔ_ｊＧＯＣに基づいて、補正係数Ｒ_ｊ１〜Ｒ_{ｊｌａｓｔ}を算出してもよい。

また、上記の各実施形態では、計測装置１は、式（５９）により変位ＣＵ_ｅｓｔ（ｔ）を算出しているが、所定の相関式により、変位ＣＵ_ｅｓｔ（ｔ）を荷重に変換してもよい。例えば、上部構造７の観測位置Ｒ_ｋの荷重ＣＷ_ｋ（ｔ）と変位ｘ_ｋ（ｔ）の関係は式（７９）ように表される。ここで、荷重ＣＷ_ｋ（ｔ）はＢＷＩＭにおいて変位波形が対応する荷重波形とする。式（７９）の１次係数Ｓｃ_ｋｋ及び０次係数Ｉｃ_ｋは、複数の試験車両による荷重試験によって得られる。

式（７９）において、Ｉｃ_ｋが十分に小さいものとすると、式（８０）が得られる。

式（８０）において、変位ｘ_ｋ（ｔ）を変位ＣＵ_ｅｓｔ（ｔ）に置き換えて、荷重ＣＷ_ｋ（ｔ）と変位ＣＵ_ｅｓｔ（ｔ）の相関式は、式（８１）のようになる。計測装置１は、相関式（８１）により、変位ＣＵ_ｅｓｔ（ｔ）を荷重ＣＷ_ｋ（ｔ）に変換することができる。

また、上記の各実施形態では観測点Ｐ_１〜Ｐ_Ｎを観測する観測装置及び観測点Ｑ_１〜Ｑ_Ｎを観測する観測装置は、それぞれ加速度センサーであるが、これに限られず、例えば、衝撃センサー、マイクロホン、歪計又はロードセルであってもよい。観測装置と観測点が１対１に対応している必要はなく、１つの観測装置が観測点Ｐ_１〜Ｐ_Ｎ，Ｑ_１〜Ｑ_Ｎの一部又は全部を観測してもよい。

衝撃センサーは、車両６の各車軸の観測点Ｐ_１〜Ｐ_Ｎ，Ｑ_１〜Ｑ_Ｎへの作用に対する応答として衝撃加速度を検出する。計測装置１は、観測点Ｐ_１〜Ｐ_Ｎに対する衝撃加速度に基づいて第１観測点情報を取得し、観測点Ｑ_１〜Ｑ_Ｎに対する衝撃加速度に基づいて第２観測点情報を取得する。マイクロホンは、車両６の各車軸の観測点Ｐ_１〜Ｐ_Ｎ，Ｑ_１〜Ｑ_Ｎへの作用に対する応答として音響を検出する。計測装置１は、観測点Ｐ_１〜Ｐ_Ｎに対する音響に基づいて第１観測点情報を取得し、観測点Ｑ_１〜Ｑ_Ｎに対する音響に基づいて第２観測点情報を取得する。歪計、ロードセルは、車両６の各車軸の観測点Ｐ_１〜Ｐ_Ｎ，Ｑ_１〜Ｑ_Ｎへの作用に対する応答として応力変化を検出する。計測装置１は、観測点Ｐ_１〜Ｐ_Ｎの応力変化に基づいて第１観測点情報を取得し、観測点Ｑ_１〜Ｑ_Ｎの応力変化に基づいて第２観測点情報を取得する。

また、上記の各実施形態では、車両６がレーンＬ_１〜Ｌ_Ｎを走行する方向はすべて同じであるが、レーンＬ_１〜Ｌ_Ｎのうちの少なくとも１つのレーンと他のレーンとは車両６の走行方向が異なっていてもよい。例えば、車両６がレーンＬ_１では観測点Ｐ_１から観測点Ｑ_１に向かう方向に走行し、レーンＬ_２では観測点Ｑ_２から観測点Ｐ_２に向かう方向に走行してもよい。この場合、計測装置１は、観測点Ｐ_１を観測するセンサー２１から出力される加速度データに基づいて、車両６のレーンＬ_１への進入時刻を取得し、観測点Ｑ_１を観測するセンサー２２から出力される加速度データに基づいて、車両６のレーンＬ_１からの退出時刻を取得する。また、計測装置１は、観測点Ｑ_２を観測するセンサー２２から出力される加速度データに基づいて、車両６のレーンＬ_２への進入時刻を取得し、観測点Ｐ_２を観測するセンサー２１から出力される加速度データに基づいて、車両６のレーンＬ_２からの退出時刻を取得する。

また、上記の各実施形態では、各センサー２１，２２は、それぞれ上部構造７の主桁Ｇに設けられているが、上部構造７の表面や内部、床板Ｆの下面、橋脚８ａ等に設けられていてもよい。また、上記の各実施形態では、橋梁５として道路橋を例に挙げたが、これに限られず、例えば、橋梁５は鉄道橋であってもよい。また、上記の各実施形態では、構造物として橋梁の上部構造を例に挙げたが、これに限られず、構造物は移動体の移動によって変形するものであればよい。

上述した実施形態および変形例は一例であって、これらに限定されるわけではない。例えば、各実施形態および各変形例を適宜組み合わせることも可能である。

本発明は、実施の形態で説明した構成と実質的に同一の構成、例えば、機能、方法及び結果が同一の構成、あるいは目的及び効果が同一の構成を含む。また、本発明は、実施の形態で説明した構成の本質的でない部分を置き換えた構成を含む。また、本発明は、実施の形態で説明した構成と同一の作用効果を奏する構成又は同一の目的を達成することができる構成を含む。また、本発明は、実施の形態で説明した構成に公知技術を付加した構成を含む。

１…計測装置、２…サーバー、４…通信ネットワーク、５…橋梁、６…車両、７…上部構造、７ａ…橋床、７ｂ…支承、Ｇ…主桁、Ｆ…床板、８…下部構造、８ａ…橋脚、８ｂ…橋台、１０…計測システム、２１…センサー、２２…センサー、１１０…制御部、１１１…第１観測点情報取得部、１１２…第２観測点情報取得部、１１３…第１時間算出部、１１４…補正係数算出部、１１５…たわみ波形算出部、１１６…移動体たわみ波形算出部、１１７…変位算出部、１１８…係数値算出部、１１９…出力処理部、１２０…第１通信部、１３０…記憶部、１３１…計測プログラム、１４０…第２通信部、１５０…操作部、１６１…第２時間算出部

Claims (15)

1. 移動体が構造物を移動する第１方向に沿って並ぶ前記構造物の第１の観測点及び第２の観測点のうちの前記第１の観測点を観測する観測装置による観測情報に基づいて、前記移動体の複数の部位が前記第１の観測点をそれぞれ通過した時刻及び前記複数の部位のそれぞれの前記第１の観測点への作用に対する応答である第１の物理量を含む第１観測点情報を取得する第１観測点情報取得ステップと、
   前記第２の観測点を観測する観測装置による観測情報に基づいて、前記複数の部位が前記第２の観測点をそれぞれ通過した時刻及び前記複数の部位のそれぞれの前記第２の観測点への作用に対する応答である第２の物理量を含む第２観測点情報を取得する第２観測点情報取得ステップと、
   前記第１観測点情報を用いて、前記複数の部位のうちの先頭の部位が前記第１の観測点を通過した時刻である第１の先頭時刻から前記複数の部位が前記第１の観測点をそれぞれ通過した時刻までの時間と、前記第１の先頭時刻から前記第１の物理量の総和を所定の配分比で配分する時刻である第１の基準時刻までの時間とを算出する第１時間算出ステップと、
   前記第１の先頭時刻から前記複数の部位が前記第１の観測点をそれぞれ通過した時刻までの時間と、前記第１の先頭時刻から前記第１の基準時刻までの時間とに基づいて、前記第１の物理量を補正する補正係数を算出する補正係数算出ステップと、
   前記第１観測点情報及び前記第２観測点情報と、所定の係数と、前記補正係数と、前記構造物のたわみの近似式とに基づいて、前記複数の部位のそれぞれによる前記構造物のたわみ波形を算出するたわみ波形算出ステップと、
   前記たわみ波形算出ステップで算出した前記複数の部位のそれぞれによる前記構造物のたわみ波形を加算して、前記移動体による前記構造物のたわみ波形を算出する移動体たわみ波形算出ステップと、を含む、計測方法。
2. 前記第２観測点情報を用いて、前記複数の部位のうちの前記先頭の部位が前記第２の観測点を通過した時刻である第２の先頭時刻から前記複数の部位が前記第２の観測点をそれぞれ通過した時刻までの時間と、前記第２の先頭時刻から前記第２の物理量の総和を前記所定の配分比で配分する時刻である第２の基準時刻までの時間とを算出する第２時間算出ステップをさらに含み、
   前記補正係数算出ステップでは、
   前記第１の先頭時刻から前記複数の部位が前記第１の観測点をそれぞれ通過した時刻までの時間と、前記第１の先頭時刻から前記第１の基準時刻までの時間と、前記第２の先頭時刻から前記複数の部位が前記第２の観測点をそれぞれ通過した時刻までの時間と、前記第２の先頭時刻から前記第２の基準時刻までの時間とに基づいて、前記第１の物理量及び前記第２の物理量を補正する前記補正係数を算出する、請求項１に記載の計測方法。
3. 前記配分比をＲ：１−Ｒとしたとき、Ｒは０．４以上０．６以下である、請求項１又は２のいずれか一項に記載の計測方法。
4. 前記第１の観測点は前記構造物の第１の端部に設定され、
   前記第２の観測点は前記構造物の前記第１の端部とは異なる第２の端部に設定される、請求項１乃至３のいずれか一項に記載の計測方法。
5. 前記構造物は、橋梁の上部構造であって、
   前記上部構造は、隣り合う橋台と橋脚、隣り合う２つの橋台、又は、隣り合う２つの橋脚のいずれか１つに渡された構造であり、
   前記上部構造の両端部は、前記隣り合う橋台と橋脚の位置、前記隣り合う２つの橋台の位置、又は、前記隣り合う２つの橋脚の位置にあり、
   前記橋梁は、道路橋又は鉄道橋である、請求項１乃至４のいずれか一項に記載の計測方法。
6. 前記移動体は、鉄道車両、自動車、路面電車、建設車両、又は軍用車両であり、
   前記複数の部位のそれぞれは車軸又は車輪である、請求項１乃至５のいずれか一項に記載の計測方法。
7. 前記構造物のたわみの近似式は、前記構造物の構造モデルに基づく式である、請求項１乃至６のいずれか一項に記載の計測方法。
8. 前記構造モデルは、両端を支持した単純梁である、請求項７に記載の計測方法。
9. 前記構造物のたわみの近似式は、前記第１の観測点と前記第２の観測点との中央位置におけるたわみの最大振幅で規格化された式である、請求項１乃至８のいずれか一項に記載の計測方法。
10. 前記第１の観測点を観測する前記観測装置及び前記第２の観測点を観測する前記観測装置は、加速度センサーである、請求項１乃至９のいずれか一項に記載の計測方法。
11. 前記第１の観測点を観測する前記観測装置及び前記第２の観測点を観測する前記観測装置は、衝撃センサー、マイクロホン、歪計又はロードセルである、請求項１乃至９のいずれか一項に記載の計測方法。
12. 前記構造物は、ＢＷＩＭ（Bridge Weigh in Motion）が機能する構造である、請求項１乃至１１のいずれか一項に記載の計測方法。
13. 移動体が構造物を移動する第１方向に沿って並ぶ前記構造物の第１の観測点及び第２の観測点のうちの前記第１の観測点を観測する観測装置による観測情報に基づいて、前記移動体の複数の部位が前記第１の観測点をそれぞれ通過した時刻及び前記複数の部位のそれぞれの前記第１の観測点への作用に対する応答である第１の物理量を含む第１観測点情報を取得する第１観測点情報取得部と、
    前記第２の観測点を観測する観測装置による観測情報に基づいて、前記複数の部位が前記第２の観測点をそれぞれ通過した時刻及び前記複数の部位のそれぞれの前記第２の観測点への作用に対する応答である第２の物理量を含む第２観測点情報を取得する第２観測点情報取得部と、
    前記第１観測点情報を用いて、前記複数の部位のうちの先頭の部位が前記第１の観測点を通過した時刻である第１の先頭時刻から前記複数の部位が前記第１の観測点をそれぞれ通過した時刻までの時間と、前記第１の先頭時刻から前記第１の物理量の総和を所定の配分比で配分する時刻である第１の基準時刻までの時間とを算出する第１時間算出部と、
    前記第１の先頭時刻から前記複数の部位が前記第１の観測点をそれぞれ通過した時刻までの時間と、前記第１の先頭時刻から前記第１の基準時刻までの時間とに基づいて、前記第１の物理量を補正する補正係数を算出する補正係数算出部と、
    前記第１観測点情報及び前記第２観測点情報と、所定の係数と、前記補正係数と、前記構造物のたわみの近似式とに基づいて、前記複数の部位のそれぞれによる前記構造物のたわみ波形を算出するたわみ波形算出部と、
    前記たわみ波形算出部が算出した前記複数の部位のそれぞれによる前記構造物のたわみ波形を加算して、前記移動体による前記構造物のたわみ波形を算出する移動体たわみ波形算出部と、を含む、計測装置。
14. 請求項１３に記載の計測装置と、
    前記第１の観測点を観測する前記観測装置と、
    前記第２の観測点を観測する前記観測装置と、を備えた、計測システム。
15. 移動体が構造物を移動する第１方向に沿って並ぶ前記構造物の第１の観測点及び第２の観測点のうちの前記第１の観測点を観測する観測装置による観測情報に基づいて、前記移動体の複数の部位が前記第１の観測点をそれぞれ通過した時刻及び前記複数の部位のそれぞれの前記第１の観測点への作用に対する応答である第１の物理量を含む第１観測点情報を取得する第１観測点情報取得ステップと、
    前記第２の観測点を観測する観測装置による観測情報に基づいて、前記複数の部位が前記第２の観測点をそれぞれ通過した時刻及び前記複数の部位のそれぞれの前記第２の観測点への作用に対する応答である第２の物理量を含む第２観測点情報を取得する第２観測点情報取得ステップと、
    前記第１観測点情報を用いて、前記複数の部位のうちの先頭の部位が前記第１の観測点を通過した時刻である第１の先頭時刻から前記複数の部位が前記第１の観測点をそれぞれ通過した時刻までの時間と、前記第１の先頭時刻から前記第１の物理量の総和を所定の配分比で配分する時刻である第１の基準時刻までの時間とを算出する第１時間算出ステップと、
    前記第１の先頭時刻から前記複数の部位が前記第１の観測点をそれぞれ通過した時刻までの時間と、前記第１の先頭時刻から前記第１の基準時刻までの時間とに基づいて、前記第１の物理量を補正する補正係数を算出する補正係数算出ステップと、
    前記第１観測点情報及び前記第２観測点情報と、所定の係数と、前記補正係数と、前記構造物のたわみの近似式とに基づいて、前記複数の部位のそれぞれによる前記構造物のたわみ波形を算出するたわみ波形算出ステップと、
    前記たわみ波形算出ステップで算出した前記複数の部位のそれぞれによる前記構造物のたわみ波形を加算して、前記移動体による前記構造物のたわみ波形を算出する移動体たわみ波形算出ステップと、をコンピューターに実行させる、計測プログラム。

Images