JP2021147825A - Measurement method, measurement apparatus, measurement system, and measurement program - Google Patents
Measurement method, measurement apparatus, measurement system, and measurement program Download PDFInfo
- Publication number
- JP2021147825A JP2021147825A JP2020047304A JP2020047304A JP2021147825A JP 2021147825 A JP2021147825 A JP 2021147825A JP 2020047304 A JP2020047304 A JP 2020047304A JP 2020047304 A JP2020047304 A JP 2020047304A JP 2021147825 A JP2021147825 A JP 2021147825A
- Authority
- JP
- Japan
- Prior art keywords
- observation
- action
- physical quantity
- observation points
- observation point
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Abstract
Description
本発明は、計測方法、計測装置、計測システム及び計測プログラムに関する。 The present invention relates to a measuring method, a measuring device, a measuring system and a measuring program.
特許文献1には、橋梁の維持管理をする上で、橋梁を通過する大型車両の車軸重量が、橋梁の損傷を予測するために重要な情報であって、この軸重測定ため、橋梁の主桁に設置したひずみ計から車両通過時のひずみ値を連続測定し、軸重を算出する手法Weight In Motionが提案されており、橋梁の主桁に配置されたひずみ計で計測したひずみ波形に基づいて、橋梁を通過する車両の車重を計測する橋梁通過車両監視システムが記載されている。詳細には、橋梁通過車両監視システムは、ひずみ計を配置して、ひずみ計で計測したひずみ波形から車軸の通過タイミングを検出して車両の軸間比率を算出し、算出した軸間比率と軸間距離データベースに登録された軸間距離から算出される軸間比率とを比較し、車両の軸間距離、車速および車種を特定する。また、橋梁通過車両監視システムは、車軸の通過タイミングに合わせて、基準軸重ひずみ波形を時間軸上に配置したひずみ波形を生成し、基準軸重ひずみ波形とひずみ計で計測したひずみ波形とを比較して各軸の軸重を算出する。そして、橋梁通過車両監視システムは、各軸の軸重を合計することにより車重を算出する。
In
しかしながら、特許文献1に記載のシステムでは、車両の車重を計測することはできるが、システムの異常や車両等の移動体が移動する橋梁等の構造物の異常を判断するために利用可能な指標を算出することはできない。
However, although the system described in
本発明に係る計測方法の一態様は、
2以上の整数N及び前記整数Nよりも大きい整数Mに対して、移動体が構造物の第1〜第Nの経路のいずれかを移動する第1方向と交差する第2方向に沿って並ぶ前記構造物の第1〜第Mの観測点を観測する少なくとも1つの観測装置による観測情報に基づいて、前記第1〜第Mの観測点の物理量を取得する物理量取得ステップと、
1以上M以下の任意の整数i及び1以上M以下の任意の整数jに対して、前記第jの観測点の作用xjと、前記作用xjが前記第iの観測点に及ぼす作用との相関を示す関数をyijとしたとき、前記物理量取得ステップで取得した前記第iの観測点の前記物理量が、関数yi1〜yiMの値の和に等しいものとして、前記物理量取得ステップで取得した前記第1〜第Mの観測点の前記物理量のうちのN個の前記物理量を複数組選択し、前記複数組のそれぞれのN個の前記物理量に基づいて、前記第1〜第Mの観測点のうちの前記第1〜第Nの経路に対応付けられたN個の観測点の作用をそれぞれ複数算出する作用算出ステップと、
前記N個の観測点のそれぞれに対して、前記作用算出ステップで算出した複数の前記作用の統計値を算出する統計値算出ステップと、を含む。
One aspect of the measurement method according to the present invention is
For an integer N of 2 or more and an integer M larger than the integer N, the moving body is arranged along a second direction intersecting the first direction in which the moving body travels in any of the first to Nth paths of the structure. A physical quantity acquisition step of acquiring the physical quantity of the first to Mth observation points based on the observation information by at least one observation device for observing the first to Mth observation points of the structure.
For one or more M or less arbitrary integer i and 1 to M any integer j, the effect x j of the observation point of the first j, and action of the working x j is on the observation point of the i-th Assuming that the physical quantity of the i-th observation point acquired in the physical quantity acquisition step is equal to the sum of the values of the functions y i1 to y iM , when the function showing the correlation of is y ij, the physical quantity acquisition step A plurality of sets of N physical quantities among the acquired physical quantities of the first to M observation points are selected, and based on each of the N physical quantities of the plurality of sets, the first to Mth. An action calculation step for calculating a plurality of actions of each of the N observation points associated with the first to Nth paths among the observation points, and an action calculation step.
For each of the N observation points, a statistical value calculation step for calculating a plurality of statistical values of the action calculated in the action calculation step is included.
前記計測方法の一態様は、
前記統計値算出ステップで前記N個の観測点のそれぞれに対して算出した前記統計値に基づいて、異常を判断する異常判断ステップを含んでもよい。
One aspect of the measurement method is
An abnormality determination step for determining an abnormality may be included based on the statistical value calculated for each of the N observation points in the statistical value calculation step.
前記計測方法の一態様において、
前記統計値は、最尤値、分散又は中間値であってもよい。
In one aspect of the measurement method
The statistical value may be a maximum likelihood value, a variance, or an intermediate value.
前記計測方法の一態様は、
前記移動体とは異なる既知の移動体が単独で前記構造物を移動したときの前記第1〜第Mの観測点の物理量を取得し、前記第1〜第Mの観測点の前記物理量に基づいて、前記関数yijの係数の値を算出する係数値算出ステップを含んでもよい。
One aspect of the measurement method is
Obtain the physical quantity of the first to third observation points when a known moving body different from the moving body moves the structure alone, and based on the physical quantity of the first to M observation points. The coefficient value calculation step for calculating the coefficient value of the function y ij may be included.
前記計測方法の一態様において、
前記関数yijは、前記作用xjの多項式関数であってもよい。
In one aspect of the measurement method
The function y ij may be a polynomial function of the action x j.
前記計測方法の一態様において、
前記物理量取得ステップで取得する前記第1〜第Mの観測点の前記物理量は、変位、又は前記移動体による荷重であってもよい。
In one aspect of the measurement method
The physical quantity of the first to Mth observation points acquired in the physical quantity acquisition step may be a displacement or a load due to the moving body.
前記計測方法の一態様において、
前記観測装置は、加速度センサーであってもよい。
In one aspect of the measurement method
The observation device may be an acceleration sensor.
前記計測方法の一態様において、
前記観測装置は、接触式変位計、リング式変位計、レーザー変位計、感圧センサー、画像処理による変位計測機器又は光ファイバーによる変位計測機器であってもよい。
In one aspect of the measurement method
The observation device may be a contact type displacement meter, a ring type displacement meter, a laser displacement meter, a pressure sensitive sensor, a displacement measuring device by image processing, or a displacement measuring device by an optical fiber.
前記計測方法の一態様において、
前記物理量取得ステップで取得する前記第1〜第Mの観測点の前記物理量は、前記第1方向及び前記第2方向とそれぞれ交差する第3方向の前記物理量であってもよい。
In one aspect of the measurement method
The physical quantity of the first to Mth observation points acquired in the physical quantity acquisition step may be the physical quantity in the third direction intersecting the first direction and the second direction, respectively.
前記計測方法の一態様において、
前記移動体は、鉄道車両、自動車、路面電車、建設車両、又は軍用車両であってもよい。
In one aspect of the measurement method
The moving vehicle may be a railroad vehicle, an automobile, a tram, a construction vehicle, or a military vehicle.
前記計測方法の一態様において、
前記構造物は、橋梁の上部構造であって、
前記上部構造は、隣り合う橋台と橋脚、隣り合う2つの橋台、又は、隣り合う2つの橋脚のいずれか1つに渡された構造であり、
前記上部構造の両端部は、前記隣り合う橋台と橋脚の位置、前記隣り合う2つの橋台の位置、又は、前記隣り合う2つの橋脚の位置にあり、
前記橋梁は、道路橋又は鉄道橋であってもよい。
In one aspect of the measurement method
The structure is a superstructure of a bridge and
The superstructure is a structure passed to any one of adjacent abutments and piers, two adjacent abutments, or two adjacent piers.
Both ends of the superstructure are located at the positions of the adjacent abutments and piers, the positions of the two adjacent abutments, or the positions of the two adjacent piers.
The bridge may be a road bridge or a railway bridge.
前記計測方法の一態様において、
前記構造物は、BWIM(Bridge Weigh in Motion)が機能する構造であってもよい。
In one aspect of the measurement method
The structure may be a structure in which BWIM (Bridge Weigh in Motion) functions.
本発明に係る計測装置の一態様は、
2以上の整数N及び前記整数Nよりも大きい整数Mに対して、移動体が構造物の第1〜第Nの経路のいずれかを移動する第1方向と交差する第2方向に沿って並ぶ前記構造物の第1〜第Mの観測点を観測する少なくとも1つの観測装置による観測情報に基づいて、前記第1〜第Mの観測点の物理量を取得する物理量取得部と、
1以上M以下の任意の整数i及び1以上M以下の任意の整数jに対して、前記第jの観測点の作用xjと、前記作用xjが前記第iの観測点に及ぼす作用との相関を示す関数をyijとしたとき、前記物理量取得部が取得した前記第iの観測点の前記物理量が、関数yi1〜yiMの値の和に等しいものとして、前記物理量取得部が取得した前記第1〜第Mの観測点の前記物理量のうちのN個の前記物理量を複数組選択し、前記複数組のそれぞれのN個の前記物理量に基づいて、前記第1〜第Mの観測点のうちの前記第1〜第Nの経路に対応付けられたN個の観測点の作用をそれぞれ複数算出する作用算出部と、
前記N個の観測点のそれぞれに対して、前記作用算出部が算出した複数の前記作用の統計値を算出する統計値算出部と、を含む。
One aspect of the measuring device according to the present invention is
For an integer N of 2 or more and an integer M larger than the integer N, the moving body is arranged along a second direction intersecting the first direction in which the moving body travels in any of the first to Nth paths of the structure. A physical quantity acquisition unit that acquires the physical quantity of the first to Mth observation points based on the observation information by at least one observation device that observes the first to Mth observation points of the structure.
For one or more M or less arbitrary integer i and 1 to M any integer j, the effect x j of the observation point of the first j, and action of the working x j is on the observation point of the i-th when the function indicating the correlation of the y ij, the physical quantity of the observation point of the i-th said physical quantity acquisition unit has acquired, as equal to the sum of the values of the function y i1 ~y iM, said physical quantity obtaining unit A plurality of sets of N physical quantities among the acquired physical quantities of the first to M observation points are selected, and based on each of the N physical quantities of the plurality of sets, the first to Mth. An action calculation unit that calculates a plurality of actions of each of the N observation points associated with the first to Nth paths among the observation points, and an action calculation unit.
For each of the N observation points, a statistical value calculation unit for calculating a plurality of statistical values of the action calculated by the action calculation unit is included.
本発明に係る計測システムの一態様は、
前記計測装置の一態様と、
前記観測装置と、を備える。
One aspect of the measurement system according to the present invention is
One aspect of the measuring device and
The observation device is provided.
本発明に係る計測プログラムの一態様は、
2以上の整数N及び前記整数Nよりも大きい整数Mに対して、移動体が構造物の第1〜第Nの経路のいずれかを移動する第1方向と交差する第2方向に沿って並ぶ前記構造物の第1〜第Mの観測点を観測する少なくとも1つの観測装置による観測情報に基づいて、前記第1〜第Mの観測点の物理量を取得する物理量取得ステップと、
1以上M以下の任意の整数i及び1以上M以下の任意の整数jに対して、前記第jの観測点の作用xjと、前記作用xjが前記第iの観測点に及ぼす作用との相関を示す関数をyijとしたとき、前記物理量取得ステップで取得した前記第iの観測点の前記物理量が、関数yi1〜yiMの値の和に等しいものとして、前記物理量取得ステップで取得した前記第1〜第Mの観測点の前記物理量のうちのN個の前記物理量を複数組選択し、前記複数組のそれぞれのN個の前記物理量に基づいて、前記第1〜第Mの観測点のうちの前記第1〜第Nの経路に対応付けられたN個の観測点の作用をそれぞれ複数算出する作用算出ステップと、
前記N個の観測点のそれぞれに対して、前記作用算出ステップで算出した複数の前記作用の統計値を算出する統計値算出ステップと、をコンピューターに実行させる。
One aspect of the measurement program according to the present invention is
For an integer N of 2 or more and an integer M larger than the integer N, the moving body is arranged along a second direction intersecting the first direction in which the moving body travels in any of the first to Nth paths of the structure. A physical quantity acquisition step of acquiring the physical quantity of the first to Mth observation points based on the observation information by at least one observation device for observing the first to Mth observation points of the structure.
For one or more M or less arbitrary integer i and 1 to M any integer j, the effect x j of the observation point of the first j, and action of the working x j is on the observation point of the i-th Assuming that the physical quantity of the i-th observation point acquired in the physical quantity acquisition step is equal to the sum of the values of the functions y i1 to y iM , when the function showing the correlation of is y ij, the physical quantity acquisition step A plurality of sets of N physical quantities among the acquired physical quantities of the first to M observation points are selected, and based on each of the N physical quantities of the plurality of sets, the first to Mth. An action calculation step for calculating a plurality of actions of each of the N observation points associated with the first to Nth paths among the observation points, and an action calculation step.
For each of the N observation points, a computer is made to execute a statistical value calculation step for calculating a plurality of statistical values of the action calculated in the action calculation step.
以下、本発明の好適な実施形態について図面を用いて詳細に説明する。なお、以下に説明する実施の形態は、特許請求の範囲に記載された本発明の内容を不当に限定するものではない。また以下で説明される構成の全てが本発明の必須構成要件であるとは限らない。 Hereinafter, preferred embodiments of the present invention will be described in detail with reference to the drawings. The embodiments described below do not unreasonably limit the content of the present invention described in the claims. Moreover, not all of the configurations described below are essential constituent requirements of the present invention.
1.第1実施形態
1−1.計測システム
以下では、構造物が橋梁の上部構造であり、移動体が車両である場合を例に挙げ、本実施形態の計測方法を実現するための計測システムについて説明する。本実施形態に係る橋梁を通過する車両は、鉄道車両、自動車、路面電車、建設車両、又は軍用車両等の重量が大きく、BWIM(Bridge Weigh in Motion)で計測可能な車両である。BWIMは、橋梁を「はかり」に見立て、橋梁の変形を計測することにより、橋梁を通行する車両の重量、軸数などを測定する技術である。変形やひずみなどの応答から通行する車両の重量を解析可能な橋梁の上部構造は、BWIMが機能する構造物であり、橋梁の上部構造への作用と応答の間の物理的なプロセスを応用するBWIMシステムが通行する車両の重量の計測を可能にする。
1. 1. First Embodiment 1-1. Measurement system In the following, a measurement system for realizing the measurement method of the present embodiment will be described by taking as an example a case where the structure is the superstructure of a bridge and the moving body is a vehicle. The vehicle passing through the bridge according to the present embodiment is a vehicle having a large weight such as a railroad vehicle, an automobile, a tram, a construction vehicle, or a military vehicle, and can be measured by BWIM (Bridge Weigh in Motion). BWIM is a technology for measuring the weight, number of axes, etc. of a vehicle passing through a bridge by measuring the deformation of the bridge by using the bridge as a "scale". The bridge superstructure, which can analyze the weight of vehicles passing by from responses such as deformation and strain, is the structure on which BWIM functions and applies the physical process between the action and response of the bridge superstructure. The BWIM system enables the measurement of the weight of vehicles passing by.
図1は、本実施形態に係る計測システムの一例を示す図である。図1に示すように、本実施形態に係る計測システム10は、計測装置1と、橋梁5の上部構造7に設けられる複数のセンサー23と、を有している。また、計測システム10は、サーバー2を有してもよい。
FIG. 1 is a diagram showing an example of a measurement system according to the present embodiment. As shown in FIG. 1, the
橋梁5は上部構造7と下部構造8からなり、上部構造7は、床板F、主桁G、不図示の横桁等からなる橋床7aと、支承7bと、を含む。下部構造8は、橋脚8aと、橋台8bと、を含む。上部構造7は、隣り合う橋台8bと橋脚8a、隣り合う2つの橋台8b、又は、隣り合う2つの橋脚8aのいずれか1つに渡された構造である。上部構造7の両端部は、隣り合う橋台8bと橋脚8aの位置、隣り合う2つの橋台8bの位置、又は、隣り合う2つの橋脚8aの位置にある。
The
計測装置1と各センサー23とは、例えば、不図示のケーブルで接続され、CAN(Controller Area Network)等の通信ネットワークを介して通信を行う。あるいは、計測装置1と各センサー23とは、無線ネットワークを介して通信を行ってもよい。
The measuring
例えば、各センサー23は、移動体である車両6の移動による上部構造7の変位を算出するためのデータを出力する。本実施形態では、各センサー23は加速度センサーであり、例えば、水晶加速度センサーであってもよいし、MEMS(Micro Electro Mechanical Systems)加速度センサーであってもよい。
For example, each
本実施形態では、各センサー23は上部構造7の長手方向の中央部に設置されている。ただし、各センサー23は、上部構造7の変位を算出するための加速度を検出することができればよく、その設置位置は上部構造7の中央部に限定されない。
In the present embodiment, each
上部構造7の床板Fや主桁G等は、上部構造7を走行する車両6による荷重によって、垂直方向下方に撓む。各センサー23は、上部構造7を走行する車両6の荷重による床板Fや主桁Gの撓みの加速度を検出する。
The floor plate F, main girder G, and the like of the
計測装置1は、各センサー23から出力される加速度データに基づいて、車両6の走行による上部構造7の撓みの変位を算出する。また、計測装置1は、算出した変位から、上部構造7を走行する車両6による荷重を算出する。
The measuring
計測装置1とサーバー2とは、例えば、携帯電話の無線ネットワーク及びインターネット等の通信ネットワーク4を介して、通信を行うことができる。計測装置1は、車両6が上部構造7を走行した時刻や車両6の走行による上部構造7の変位等の情報をサーバー2に送信する。サーバー2は、当該情報を不図示の記憶装置に記憶し、例えば、当該情報に基づいて過積載の車両の監視や上部構造7の異常判定等の処理を行ってもよい。
The measuring
なお、本実施形態では、橋梁5は、道路橋であり、例えば、鋼橋や桁橋、RC(Reinforced-Concrete)橋等である。
In the present embodiment, the
図2及び図3は、各センサー23の上部構造7への設置例を示す図である。なお、図2は、上部構造7をその上方から見た図であり、図3は、図2をA−A線で切断した断面図である。
2 and 3 are views showing an example of installation of each
図2及び図3に示すように、上部構造7は、移動体である車両6が移動し得る第1〜第Nの経路としてのN個のレーンL1〜LN、及びK個の主桁G1〜GKを有している。ここで、Nは2以上の整数であり、Kは1以上の整数である。なお、図2及び図3の例では、主桁G1〜GKの各位置がレーンL1〜LNの各境界の位置と一致しており、N=K−1であるが、主桁G1〜GKの各位置がレーンL1〜LNの各境界の位置と一致している必要はなく、N≠K−1であってもよい。
As shown in FIGS. 2 and 3, the
図2及び図3の例では、上部構造7の長手方向の中央部CAにおいて、主桁G1〜GKのそれぞれにM個のセンサー23のそれぞれが設けられている。ここで、整数Mは整数Nよりも大きい。図2及び図3の例では、M=K=N+1であるが、K>M>Nであり、主桁G1〜GKの少なくとも1つにセンサー23が設けられていなくてもよい。
In the example of FIG. 2 and FIG. 3, in the longitudinal direction of the central portion CA of the
なお、各センサー23を上部構造7の床板Fに設けると、走行車両によって破壊するおそれがあり、また橋床7aの局部的な変形により測定精度が影響を受けるおそれがあるため、本実施形態では、図2及び図3の例では、各センサー23は上部構造7の主桁G1〜GKに設けられている。
If each
本実施形態では、M個のセンサー23に対応付けてM個の観測点R1〜RMがそれぞれ設定されている。観測点R1〜RMは、車両6が上部構造7のレーンL1〜LNのいずれかを移動する第1方向と交差する第2方向に沿って並ぶ上部構造7のM個の観測点である。図2及び図3の例では、1以上M以下の各整数jに対して、観測点Rjは、中央部CAにおいて、主桁Gjに設けられたセンサー23の鉛直上方向にある床板Fの表面の位置に設定されている。すなわち、主桁Gjに設けられたセンサー23は、観測点Rjを観測する観測装置である。観測点Rjを観測するセンサー23は、車両6の走行により観測点Rjに生じる加速度を検出可能な位置に設けられていればよいが、観測点Rjに近い位置に設けられることが望ましい。このように、観測点R1〜RMはM個のセンサー23と1対1の関係にある。
In the present embodiment, in association with the
本実施形態では、M個の観測点R1〜RMのうちのN個の観測点は、それぞれレーンL1〜LNに対応付けられている。また、M個の観測点R1〜RMのうちのM−N個の観測点は、レーンL1〜LNのいずれにも対応付けられない冗長な観測点である。図2及び図3の例では、第1方向は、上部構造7のレーンL1〜LNに沿うX方向、すなわち、上部構造7の長手方向である。また、第2方向は、車両6が走行する上部構造7の面内においてX方向と直交するY方向、すなわち、上部構造7の幅員方向である。ただし、第2方向は、第1方向と直交していなくてもよい。例えば、上部構造7の一方の端から観測点R1〜RMまでの距離が異なっていてもよい。なお、観測点R1〜RMは、それぞれ「第1の観測点」〜「第Mの観測点」の一例である。
In this embodiment, N pieces of the observation point of the
なお、センサー23の数及び設置位置は、図2及び図3に示した例には限定されず種々の変形実施が可能である。
The number and installation positions of the
計測装置1は、各センサー23から出力される加速度データに基づいて、第1方向であるX方向と第2方向であるY方向とそれぞれ交差する第3方向の加速度を取得する。観測点R1〜RMは、X方向及びY方向と直交する方向に撓むので、計測装置1は、撓みの加速度の大きさを正確に算出するために、X方向及びY方向と直交する第3方向、すなわち、床板Fの法線方向の加速度を取得するのが望ましい。
The measuring
図4は、センサー23が検出する加速度を説明する図である。センサー23は、互いに直交する3軸の各軸方向に生じる加速度を検出する加速度センサーである。
FIG. 4 is a diagram illustrating the acceleration detected by the
車両6の走行による観測点R1〜RMの撓みの加速度を検出するために、各センサー23は、3つの検出軸であるx軸、y軸、z軸のうち、1軸が第1方向及び第2方向と交差する方向となるように設置される。図2及び図3の例では、第1方向はX方向であり、第2方向はY方向であるから、各センサー23は、1軸がX方向及びY方向と交差する方向となるように設置される。観測点R1〜RMは、X方向及びY方向と直交する方向に撓むので、撓みの加速度を正確に検出するために、理想的には、各センサー23は、1軸をX方向及びY方向と直交する方向、すなわち、床板Fの法線方向に合わせて設置される。
To detect the acceleration of deflection of the observation point R 1 to R M by the running of the
ただし、各センサー23を上部構造7に設置する場合、設置場所が傾いている場合もある。計測装置1は、各センサー23の3つの検出軸の1軸が、床板Fの法線方向に合わせて設置されなくても、概ね法線方向に向いていることで誤差は小さく無視できる。また、計測装置1は、各センサー23の3つの検出軸の1軸が、床板Fの法線方向に合わせて設置されなくても、x軸、y軸、z軸の加速度を合成した3軸合成加速度によって、各センサー23の傾斜による検出誤差の補正を行うことができる。また、各センサー23は、少なくとも鉛直方向にほぼ平行な方向に生ずる加速度、あるいは、床板Fの法線方向の加速度を検出する1軸加速度センサーであってもよい。
However, when each
以下、計測装置1が実行する本実施形態の計測方法の詳細について説明する。
Hereinafter, the details of the measurement method of the present embodiment executed by the
1−2.作用の算出
車両6がレーンL1を走行した時、車両6による荷重により、レーンL1に対応する観測点、例えば観測点R1に作用x1が生じる。そのため、観測点R1は作用x1によって変位する。この時、車両6が単独で上部構造7を走行していれば、レーンL2〜LNを走行する車両による荷重によって観測点R2〜RMに生じる作用はゼロである。しかしながら、作用x1が観測点R2〜RMに及ぼす作用が生じるため、各観測点R2〜RMも変位することになる。そのため、観測点R1を観測するセンサー23だけでなく、観測点R2〜RMをそれぞれ観測するM−1個のセンサー23も、レーンL1を走行した車両によって生じる加速度を検出することになる。
1-2. When calculating the
一例として、図5及び図6に、N=2,M=3の場合の各センサー23及び観測点R1,R2,R3の配置例を示し、図7に、図5及び図6に示す配置例の場合に、各センサー23が検出する加速度の一例を示す。
As an example, in FIGS. 5 and 6, N = 2, each
図5は、上部構造7をその上方から見た図であり、図6は、図5をA−A線で切断した断面図である。図5及び図6の例では、3個のセンサー23が、上部構造7の中央部CAにおいて主桁G1,G2,G3にそれぞれ設けられている。また、レーンL1に対応する観測点R1が、主桁G1に設けられたセンサー23の鉛直上方向にある床板Fの表面の位置に設定され、レーンL2に対応する観測点R3が、主桁G3に設けられたセンサー23の鉛直上方向にある床板Fの表面の位置に設定されている。また、レーンL1,L2のいずれにも対応しない観測点R2が、主桁G2に設けられたセンサー23の鉛直上方向にある床板Fの表面の位置に設定されている。主桁G1に設けられたセンサー23は観測点R1を観測し、主桁G2に設けられたセンサー23は観測点R2を観測し、主桁G3に設けられたセンサー23は観測点R3を観測する。
FIG. 5 is a view of the
図7は、車両6がレーンL1を単独で走行した場合に各センサー23から出力される加速度データの一例を示す図である。なお、図7の各波形は、ピークを明瞭化するために各加速度データに対してフィルター処理を施した波形である。
Figure 7 is a diagram showing an example of acceleration data output from each
観測点R1を観測するセンサー23から出力される加速度データのピークPK1は、車両6の車軸の通過を示している。
Peak PK1 of the acceleration data output from the
観測点R1を観測するセンサー23から出力される加速度データのピークPK1は、車両6による観測点R1の作用x1に対応する。一方、観測点R2を観測するセンサー23から出力される加速度データのピークPK2は、車両6による観測点R1の作用x1が観測点R2に及ぼす作用に対応する。同様に、観測点R3を観測するセンサー23から出力される加速度データのピークPK3は、車両6による観測点R1の作用x1が観測点R3に及ぼす作用に対応する。車両6がレーンL1を単独で走行するので、ピークPK1は、ピークPK3よりも大きい。また、ピークPK2は、ピークPK1よりも若干小さい。
Peak of the acceleration data output from the
仮に、車両6がレーンL1を走行する時に他の車両6がレーンL2を並走した場合、観測点R1を観測するセンサー23は、レーンL1を走行する車両6による観測点R1の作用x1と、レーンL2を走行する他の車両6による観測点R2の作用x2が観測点R1に及ぼす作用との和に対応する加速度を検出することになる。同様に、観測点R3を観測するセンサー23は、レーンL2を走行する他の車両6による観測点R3の作用x3と、レーンL1を走行する車両6による観測点R1の作用x1が観測点R3に及ぼす作用との和に対応する加速度を検出することになる。また、観測点R2を観測するセンサー23は、レーンL1を走行する車両6による観測点R1の作用x1が観測点R2に及ぼす作用と、レーンL2を走行する他の車両6による観測点R3の作用x3が観測点R2に及ぼす作用との和に対応する加速度を検出することになる。そのため、2台の車両6が並走する場合は、図7の例に対して、ピークPK1,PK2,PK3はすべて大きくなるが、ピークPK1,PK2,PK3の大小関係が変わらない場合もある。しかも、車両6の重量によってピークPK1,PK2,PK3の大きさが変わるため、ピークPK1,PK2,PK3を単純に比較しただけでは、1台の車両6が単独で走行したのか2台の車両6が並走したのかを区別することができない。
Assuming that the
一方、図7の例のように車両6が単独で走行する場合、各センサー23が検出する加速度の波形は互いにピークが異なるものの近似している。すなわち、1以上M以下の任意の整数j,iに対して、車両6が観測点Rjに対応するレーンを走行した時の観測点Rjの作用xjと、作用xjが観測点Riに及ぼす作用との間には相関がある。そこで、本実施形態では、計測装置1は、この相関を利用して、作用x1〜xMを算出し、作用x1〜xMによる上部構造7の変位や、レーンL1〜LNを走行した車両6による荷重を算出する。
On the other hand, when the
まず、車両6が観測点Rjに対応するレーンを走行した時の観測点Rjの作用xjと、観測点Rjの作用xjが観測点Riに及ぼす作用との相関を示す関数yijを、式(1)のように定義する。j,iは、それぞれ、1以上M以下の任意の整数である。式(1)において、aijは1次係数であり、bijは0次係数である。式(1)に示すように、関数yijは、作用xjの多項式関数であり、具体的には1次多項式関数である。
First, the function shown with the action x j of the observation point R j when the
例えば、一例として、図5及び図6に示した配置例のように、N=2の場合を例に挙げる。車両6がレーンL1を走行した時、レーンL1の観測点Rjに作用x1が作用して観測された物理量とレーンL1の作用x1がレーンL2の観測点R2に作用して観測された物理量との相関は図8のように1次多項式関数で示される。
For example, as an example, as in the arrangement example shown in FIGS. 5 and 6, the case of N = 2 will be taken as an example. When the
一例として、観測点R1の作用x1と、作用x1が観測点Riに及ぼす作用との相関を示す関数yi1は式(2)のように定義される。 As an example, the working x 1 observation point R 1, acting x 1 functions y i1 showing a correlation between the action on the observation point R i is defined by the equation (2).
より具体的には、作用x1と、作用x1が観測点R1〜RMに及ぼす作用との相関を示す関数y11〜yM1はそれぞれ式(3)のようになる。 More specifically, the working x 1, so that each formula is a function y 11 ~y M1 showing a correlation between action acting x 1 is on the observation point R 1 ~R M (3).
次に、式(4)に示すように、観測点Riの変位giは、関数yi1〜yiMの値の和に等しいものとする。 Next, as shown in Equation (4), the displacement g i of the observation point R i is assumed to be equal to the sum of the values of the function y i1 ~y iM.
このとき、式(1)及び式(4)より、観測点R1〜RMの変位g1〜gMを要素とする変位ベクトルgは、式(5)のように表される。
At this time, from equation (1) and (4), the displacement vector g for the
式(5)において、ベクトルYの各要素ykは式(6)のように定義される。kは1以上M以下の任意の整数である。 In the equation (5), each element y k of the vector Y is defined as in the equation (6). k is an arbitrary integer of 1 or more and M or less.
実際に観測される観測点R1〜RMのそれぞれの変位u1〜uMを要素とする変位ベクトルuが、変位ベクトルgと等しいものとすると、式(7)が得られる。変位u1〜uMは、例えば、観測点R1〜RMに対応するM個のセンサー23が検出する加速度をそれぞれ2回積分することによって得られる。
Displacement vector u to actually observed by each of the displacement u 1 ~u M elements of the observation points R 1 to R M is, assuming equal displacement vector g, the formula (7) is obtained.
式(7)において、M個の観測点R1〜RMのうちのM−N個の観測点は、レーンL1〜LNのいずれにも対応付けられないので、車両6がレーンL1〜LNのいずれかを走行しても、当該M−N個の観測点には車両6の走行による作用は生じないものと考える。すなわち、当該M−N個の観測点の車両6による作用はすべてゼロであるものと考える。そうすると、レーンL1〜LNに対応付けられるN個の観測点の未知の作用を算出するために、式(7)で得られるM個の変位u1〜uMから任意のN個の変位を選択して連立方程式を立てることができる。M個の変位u1〜uMから選択可能なN個の変位の組み合わせ数Cは、式(8)で表される。
In the formula (7), M-N pieces of observation points of the
計測装置1は、複数組のN個の変位の連立方程式を選択し、各組の連立方程式を解いて、車両6がレーンL1〜LNをそれぞれ単独で走行した時のレーンL1〜LNに対応するN個の観測点の作用をそれぞれ複数算出する。
Measuring
一例として、図5及び図6に示した配置例のように、N=2,M=3の場合を例に挙げて、式(7)から作用x1,x3をそれぞれ複数算出する過程を詳細に説明する。M=3であるから、式(7)より式(9)が得られる。 As an example, as in the arrangement example shown in FIGS. 5 and 6, the process of calculating a plurality of actions x 1 and x 3 from the equation (7) is performed by taking the case of N = 2 and M = 3 as an example. This will be described in detail. Since M = 3, the equation (9) can be obtained from the equation (7).
図5及び図6の例では、観測点R1,R3はそれぞれレーンL1,L2に対応付けられ、観測点R2はレーンL1,L2のいずれにも対応付けられない。したがって、式(9)において、作用x1,x3が未知数であり、作用x2はゼロであるものと考え、式(6)よりy2=0であるから、式(10)が得られる。 In the examples of FIGS. 5 and 6, observation points R 1 and R 3 are associated with lanes L 1 and L 2 , respectively, and observation points R 2 are not associated with any of lanes L 1 and L 2. Therefore, in the equation (9), it is considered that the actions x 1 and x 3 are unknown and the action x 2 is zero, and from the equation (6), y 2 = 0, so that the equation (10) can be obtained. ..
あるいは、図5及び図6の例において、観測点R1,R2,R3の変位u1,u2,u3は、観測点R1の作用x1のみによる観測点R1,R2,R3の変位u1_1,u2_1,u3_1と、観測点R2の作用x2のみによる観測点R1,R2,R3の変位u1_2,u2_2,u3_2と、観測点R3の作用x3のみによる観測点R1,R2,R3の変位u1_3,u2_3,u3_3との和に等しいものと考え、式(10)を算出することもできる。
Alternatively, in the example of FIG. 5 and FIG. 6, the displacement u 1, u 2, u 3 observation points R 1,
具体的には、作用x1のみによる観測点R1,R2,R3の変位u1_1,u2_1,u3_1を要素とする変位ベクトルg1は式(11)で表される。 Specifically, the displacement vector g 1 having the displacements u 1-11 , u 2_1 , and u 3_1 of the observation points R 1 , R 2 , and R 3 due to the action x 1 as elements is expressed by Eq. (11).
また、作用x2はゼロであるものとし、作用x2のみによる観測点R1,R2,R3の変位u1_2,u2_2,u3_2を要素とする変位ベクトルg2は式(12)で表される。 Further, it is assumed that the action x 2 is zero, and the displacement vector g 2 having the displacements u 1-2 , u 2_2 , u 3_2 of the observation points R 1 , R 2 , and R 3 due to the action x 2 as elements is given by Eq. (12). It is represented by.
作用x3のみによる観測点R1,R2,R3の変位u1_3,u2_3,u3_3を要素とする変位ベクトルg3は式(13)で表される。
Acting x 3 observation point only by R 1,
式(11)、式(12)、式(13)より、変位ベクトルuは式(14)のようになり、式(14)を変形することで式(10)が得られる。 From the equations (11), (12), and (13), the displacement vector u becomes the equation (14), and the equation (10) can be obtained by modifying the equation (14).
式(10)において、2つの未知数である作用x1,x3を算出するための、3個の変位u1,u2,u3から選択可能な変位の組み合わせは、式(8)より3通りである。
In equation (10), the combination of displacements that can be selected from the three displacements u 1 , u 2 ,
選択される第1の組は変位u1,u2であり、作用x1,x3を未知数とする連立方程式は式(15)で与えられる。 The first set to be selected is the displacement u 1 , u 2 , and the simultaneous equations with the actions x 1 , x 3 as unknowns are given by Eq. (15).
式(15)を変形し、式(16)が得られる。 Equation (15) is modified to obtain equation (16).
式(16)より、式(17)のように作用x1,x3が求められる。 From the equation (16), the actions x 1 and x 3 can be obtained as in the equation (17).
選択される第2の組は変位u1,u3であり、作用x1,x3を未知数とする連立方程式は式(18)で与えられる。 The second set to be selected is the displacement u 1 , u 3 , and the simultaneous equations with the actions x 1 , x 3 as unknowns are given by Eq. (18).
式(18)を変形し、式(19)が得られる。 Equation (18) is modified to obtain equation (19).
式(19)より、式(20)のように作用x1,x3が求められる。 From the equation (19), the actions x 1 and x 3 can be obtained as in the equation (20).
選択される第3の組は変位u2,u3であり、作用x1,x3を未知数とする連立方程式は式(21)で与えられる。
Third set chosen is the
式(21)を変形し、式(22)が得られる。 Equation (21) is modified to obtain Equation (22).
式(22)より、式(23)のように作用x1,x3が求められる。 From the equation (22), the actions x 1 and x 3 can be obtained as in the equation (23).
変位u1,u2,u3はそれぞれが時刻tをパラメータとした数列であって、作用x1,x2,x3も同様の数列である。したがって、この3組の連立方程式を解いて式(17)、式(20)、式(23)で得られる作用x1,x3をまとめると、3組の作用x1は式(24)のようになり、3組の作用x3は式(25)のようになる。 The displacements u 1 , u 2 , and u 3 are sequences with the time t as a parameter, and the actions x 1 , x 2 , and x 3 are similar sequences. Therefore, when these three sets of simultaneous equations are solved and the actions x 1 and x 3 obtained by the equations (17), (20), and (23) are put together, the three sets of actions x 1 are of the equation (24). Then, the three sets of actions x 3 are as shown in equation (25).
図9に、図5及び図6に示した配置例の場合に、車両6が単独でレーンL1を走行した場合に観測される観測点R1,R2,R3の変位の一例を示す。また、図10に、図9に示す観測点R1,R2,R3の変位に対して式(17)、式(20)又は式(23)によって得られる作用x1,x3から計算される観測点R1,R3の変位の一例を示す。図9及び図10において、横軸は時間であり、縦軸は変位である。図9において、実線は観測点R1の変位を示し、一点鎖線は観測点R2の変位を示し、破線は観測点Riの変位を示す。また、図10において、実線は観測点R1の変位を示し、破線は観測点Riの変位を示す。
FIG. 9 shows an example of the displacement of the observation points R 1 , R 2 , and R 3 observed when the vehicle 6 travels alone in the lane L 1 in the case of the arrangement example shown in FIGS. 5 and 6. .. Further, in FIG. 10, it is calculated from the actions x 1 , x 3 obtained by the equations (17), (20) or (23) with respect to the displacements of the observation points R 1 , R 2 and R 3 shown in FIG. It is the an example of a displacement of the
図9に示すように、車両6が単独でレーンL1を走行するため、観測点R1の変位は観測点R2,R3の変位よりも大きい。また、作用x1が観測点R2,R3に及ぼす作用のために、観測点R2,R3の変位はゼロになっていない。そして、観測点R1の変位が最大となる時間と観測点R2,R3の変位が最大となる時間が一致している。これに対して、図10に示すように、作用x1から計算される観測点R1の変位は、図9に示す観測点R1の変位とほぼ同じであり、作用x3から計算される観測点R3の変位はゼロである。そして、図10に示す観測点R1の変位がともにピークを持つ一方、観測点R3の変位はピークを持たないので、車両6が単独でレーンL1を走行したことが特定されるとともに、作用x1による観測点R1の変位から車両6による荷重を算出することができる。
As shown in FIG. 9, since the
1−3.統計値の算出
本実施形態では、計測装置1は、レーンL1〜LNに対応するN個の観測点の各々に対して、式(7)から選択した複数組の連立方程式をそれぞれ解いて得られる複数の作用の統計値を算出する。統計値は、例えば、最尤値、標準偏差、分散又は中間値である。
1-3. In calculating this embodiment of the statistics, the measuring
一例として、作用x1についてn個の作用x1_1(t)〜x1_n(t)が得られた場合、n個の作用x1_1(t)〜x1_n(t)の平均による最尤値x1_avg(t)は式(26)で与えられる。 As an example, n pieces of action x 1_1 the operation x 1 (t) ~x If 1_n (t) is obtained, n-number of acts x 1_1 (t) ~x 1_n ( t) the maximum likelihood value x the average of 1_avg (t) is given by equation (26).
また、n個の作用x1_1(t)〜x1_n(t)の標準偏差x1_σ(t)は式(27)で与えられる。 Further, the standard deviation x 1_Shiguma of n working x 1_1 (t) ~x 1_n ( t) (t) is given by equation (27).
また、n個の作用x1_1(t)〜x1_n(t)の分散x1_σ 2(t)は式(28)で与えられる。 The dispersibility x 1_σ 2 of n working x 1_1 (t) ~x 1_n ( t) (t) is given by equation (28).
また、n個の作用x1_1(t)〜x1_n(t)の中間値x1_medi(t)は、nが奇数の場合は、作用x1_1(t)〜x1_n(t)を大きい順又は小さい順に並べて中央値として与えられる。また、n個の作用x1_1(t)〜x1_n(t)の中間値x1_medi(t)は、nが偶数の場合は、作用x1_1(t)〜x1_n(t)を大きい順又は小さい順に並べて中央の2値の平均値として与えられる。 The intermediate value x 1_Medi of n working x 1_1 (t) ~x 1_n ( t) (t) , if n is odd, descending order of acts x 1_1 (t) ~x 1_n ( t) or It is given as the median in ascending order. The intermediate value x 1_Medi of n working x 1_1 (t) ~x 1_n ( t) (t) , if n is an even number, descending order of acts x 1_1 (t) ~x 1_n ( t) or It is given as the average value of the two values in the center, arranged in ascending order.
図11に、図5及び図6に示した配置例の場合に、式(24)で算出される3つの作用x1_1(t),x1_2(t),x1_3(t)及び式(26)で算出される平均による最尤値x1_avg(t)の一例を示す。なお、図11では、作用及び最尤値が変位に換算して示されており、横軸は時間であり、縦軸は変位である。図11において、実線は作用x1_1(t)を示し、破線は作用x1_2(t)を示し、一点鎖線は作用x1_3(t)を示し、点線は平均による最尤値x1_avg(t)を示す。図11では、3つの作用x1_1(t),x1_2(t),x1_3(t)及び平均による最尤値x1_avg(t)の互いの差が小さく、式(28)で得られる分散x1_σ 2(t)は正常な範囲である。したがって、計測装置1は、式(26)で得られる平均による最尤値x1_avg(t)を観測点R1の作用として観測点R1の変位や荷重を算出することができる。
11, in the case of the exemplary arrangement shown in FIG. 5 and FIG. 6, three actions calculated by the equation (24) x 1_1 (t) , x 1_2 (t), x 1_3 (t) and (26 ) Is shown as an example of the maximum likelihood value x 1_avg (t) calculated by the average. In FIG. 11, the action and the maximum likelihood value are shown in terms of displacement, the horizontal axis is time, and the vertical axis is displacement. 11, the solid line shows the effect x 1_1 (t), the broken line shows the effect x 1_2 (t), one-dot chain line act x 1_3 indicates (t), the maximum likelihood value x 1_Avg dotted lines by mean (t) Is shown. In Figure 11, three working x 1_1 (t), x 1_2 (t), x 1_3 (t) and the difference between the mutual maximum likelihood value x 1_avg (t) is small due to the average, is obtained by Equation (28) distributed x 1_σ 2 (t) is the normal range. Accordingly, the measuring
1−4.異常判断
本実施形態では、計測装置1は、観測点R1〜RNのうちのレーンL1〜LNに対応付けられたN個の観測点のそれぞれに対して算出した統計値に基づいて、上部構造7の異常を判断する。具体的には、計測装置1は、各レーンLjに対応する観測点に対して、複数の作用の平均による最尤値から、複数の作用の中間値を減算する。そして、計測装置1は、各レーンLjに対応する観測点に対して、複数の作用の分散値が第1の閾値を超え、かつ、最尤値から中間値を減算した値が第2の閾値を超える場合、上部構造7に異常があると判断する。
1-4. The abnormality determination embodiment, the measuring
一例として、図5及び図6に示した配置例の場合に、計測装置1が統計値に基づいて上部構造7の異常を判断する過程について説明する。
As an example, in the case of the arrangement examples shown in FIGS. 5 and 6, the process in which the
式(29)により、式(26)で得られる観測点R1の3つの作用x1_1(t),x1_2(t),x1_3(t)の平均による最尤値x1_avg(t)から、作用x1_1(t),x1_2(t),x1_3(t)の中間値x1_medi(t)を減算した波形x1_avg−medi(t)が得られる。 The equation (29), three actions of observation points R 1 obtained by Equation (26) x 1_1 (t) , x 1_2 (t), the maximum likelihood value according to the average of x 1_3 (t) x 1_avg ( t) , The waveform x 1_avg-medi (t) obtained by subtracting the median value x 1_medi (t) of the actions x 1_1 (t), x 1_2 (t), and x 1_3 (t).
同様に、式(30)により、観測点R3の3つの作用x3_1(t),x3_2(t),x3_3(t)の平均による最尤値x3_avg(t)から、作用x3_1(t),x3_2(t),x3_3(t)の中間値x3_medi(t)を減算した波形x3_avg−medi(t)が得られる。 Similarly, the equation (30), the three actions of the observation point R 3 x 3_1 (t), x 3_2 (t), the maximum likelihood value according to the average of x 3_3 (t) x 3_avg ( t), acting x 3_1 (t), x 3_2 (t ), the waveform obtained by subtracting an intermediate value x 3_medi (t) of x 3_3 (t) x 3_avg- medi (t) is obtained.
そして、計測装置1は、観測点R1に対して、式(28)で得られる3つの作用x1_1(t),x1_2(t),x1_3(t)の分散x1_σ 2(t)が第1の閾値を超えるか否かを判定する。計測装置1は、分散x1_σ 2(t)が第1の閾値を超えている場合、波形x1_avg−medi(t)が第2の閾値を超えるか否かを判定する。計測装置1は、波形x1_avg−medi(t)が第2の閾値を超えた場合、計測システム10あるいは上部構造7が異常であると判断する。同様に、計測装置1は、観測点R3に対して、3つの作用x3_1(t),x3_2(t),x3_3(t)の分散x3_σ 2(t)が第1の閾値を超え、かつ、波形x1_avg−medi(t)が第2の閾値を超えた場合、計測システム10あるいは上部構造7が異常であると判断する。
The measuring
計測装置1は、計測システム10あるいは上部構造7が異常であると判断した場合、異常であることをサーバー2等に通知する。
When the measuring
なお、式(24)及び式(25)において0次係数b11〜b33が小さく、ゼロと仮定すると、式(24)及び式(25)は、それぞれ式(31)及び式(32)のようになる。 Assuming that the zero-order coefficients b 11 to b 33 are small in the equations (24) and (25) and are zero, the equations (24) and (25) are the equations (31) and (32), respectively. Will be.
式(31)及び式(32)より、各作用は、選択した2つの観測点の変位に1次係数a11〜a33から算出される率α,βをそれぞれ乗算した差分と等しい関係にある。2つの観測点の変位の差から作用を計算するので、1次係数a11〜a33は上部構造7の作用に対する応答を決定する構造強度に比例する量と見ることができる。この場合、各作用は、選択した2つの観測点の間の関係から導かれているので、各作用を求める複数の方程式は、上部構造7の複数の部位間の構造に関する係数で構成されていると言える。したがって、算出された複数の作用の平均による最尤値は、上部構造7の平均的強度に関係していて、複数の作用の分散は上部構造7の複数の部位間の強度のばらつきに関係していると考えることができる。このことから、計測装置1は、冗長な観測点を含む複数の観測点の観測情報に基づいて複数の作用を算出し、複数の作用を統計的に比較してばらつきが大きい場合は、計測システム10あるいは上部構造7に異常があると判断することで、計測システム10及び上部構造7の診断を行うことができる。
From equations (31) and (32), each action has a relationship equal to the difference obtained by multiplying the displacements of the two selected observation points by the rates α and β calculated from the linear coefficients a 11 to a 33, respectively. .. Since the action is calculated from the difference between the displacements of the two observation points, the first-order coefficients a 11 to a 33 can be seen as a quantity proportional to the structural strength that determines the response to the action of the
1−5.計測方法
図12は、第1実施形態の計測方法の手順の一例を示すフローチャート図である。本実施形態では、計測装置1が図12に示す手順を実行する。
1-5. Measurement Method FIG. 12 is a flowchart showing an example of the procedure of the measurement method of the first embodiment. In this embodiment, the measuring
図12に示すように、まず、計測装置1は、車両が単独で上部構造7を走行したときの観測点R1〜RMの変位u1〜uMを取得し、変位u1〜uMに基づいて、関数yijの係数aij,bijの値を算出する(ステップS1)。i,jは1以上M以下の任意の整数である。車両は、未知の移動体である車両6とは異なる既知の移動体である。既知の移動体とは、荷重、寸法、軸数等の情報がわかっている移動体であり、未知の移動体とは当該情報がわかっていない移動体である。このステップS1は係数値算出ステップである。
As shown in FIG. 12, first, the
次に、計測装置1は、観測点R1〜RMを観測するM個のセンサー23による観測情報に基づいて、車両6が上部構造7を移動するときの観測点R1〜RMの物理量として変位u1〜uMを取得する(ステップS2)。前述の通り、M個のセンサー23はそれぞれ加速度センサーであり、M個のセンサー23による観測情報は、観測点R1〜RMに生じた加速度の検出情報である。そして、この加速度は、第1方向であるX方向と第2方向であるY方向とそれぞれ交差する第3方向の加速度である。計測装置1は、M個のセンサー23がそれぞれ検出した第3方向の加速度を2回積分して前述の式(7)に含まれる変位ベクトルuを算出する。したがって、物理量取得ステップで計測装置1が取得する観測点R1〜RMの物理量としての変位u1〜uMは、X方向及びY方向とそれぞれ交差する第3方向の変位、例えば、X方向及びY方向とそれぞれ直交する第3方向の変位である。このステップS2は物理量取得ステップである。
Next, the measuring
次に、計測装置1は、観測点Riの変位uiが関数yi1〜yiMの値の和に等しいものとして、ステップS2で取得した変位u1〜uMのうちのN個の変位を複数組選択し、当該複数組のそれぞれのN個の変位に基づいて、観測点R1〜RMのうちのレーンL1〜LNに対応付けられたN個の観測点の作用をそれぞれ複数算出する(ステップS3)。iは1以上M以下の任意の整数である。具体的には、計測装置1は、前述の式(7)で得られるM個の変位u1〜uMからN個の変位を複数組選択し、当該複数組の変位の連立方程式を解いて、レーンL1〜LNに対応付けられたN個の観測点の作用をそれぞれ複数算出する。このステップS3は作用算出ステップである。
Next, the measuring
次に、計測装置1は、N個の観測点のそれぞれに対して、ステップS3で算出した複数の作用の統計値を算出する(ステップS4)。このステップS4は統計値算出ステップである。
Next, the measuring
次に、計測装置1は、N個の観測点のそれぞれに対してステップS4で算出した統計値に基づいて、異常を判断する(ステップS5)。このステップS5は異常判断ステップである。
Next, the measuring
次に、計測装置1は、N個の観測点のそれぞれに対して、ステップS5で算出した統計値に基づいて、変位を算出する(ステップS6)。例えば、レーンL1〜LNに対応付けられたN個の観測点のうちの観測点Rkに対して、ステップS5で算出した複数の作用の平均による最尤値や複数の作用の中間値を作用xkとして、前述の式(5)の右辺において、作用x1〜xNのうちの作用xkを除くすべての作用をゼロとして算出される変位gkを観測点Rkの変位とする。このステップS6は変位算出ステップである。
Next, the measuring
次に、計測装置1は、ステップS6で算出したN個の観測点の変位に基づいて、レーンL1〜LNのそれぞれを走行した車両6による荷重を算出する(ステップS7)。1以上N以下の各整数jに対して、レーンLjに対応付けられた観測点の変位とレーンLjを走行した車両6による荷重との間には相関があるので、あらかじめ車両による荷重試験において、この相関式の係数を算出しておく。計測装置1は、当該相関式に、レーンLjに対応付けられた観測点の変位を代入して、レーンLjを走行した車両6による荷重を算出することができる。このステップS7は荷重算出ステップである。
Next, the measuring
次に、計測装置1は、ステップS7で算出したレーンL1〜LNを走行した車両6の荷重をサーバー2に出力する(ステップS8)。このステップS8は、出力ステップである。
Next, the measuring
計測装置1は、計測を終了するまで(ステップS9のN)、ステップS2〜S8の処理を繰り返し行う。
The measuring
図13は、図12のステップS1である係数値算出ステップの手順の一例を示すフローチャート図である。 FIG. 13 is a flowchart showing an example of the procedure of the coefficient value calculation step which is step S1 of FIG.
図13に示すように、まず、計測装置1は、整数jを1に設定し(ステップS11)、観測点R1〜RMを観測するM個のセンサー23が、車両が単独でレーンLjを走行した時に検出した加速度を取得する(ステップS12)。
As shown in FIG. 13, first, the
次に、計測装置1は、ステップS12で取得した観測点Rjの加速度と観測点R1〜RMのそれぞれの加速度との相関を1次近似し、1次係数a1j〜aMjの値及び0次係数b1j〜bMjの値を算出する(ステップS13)。
Next, the measuring
計測装置1は、整数jがNでない場合は(ステップS14のN)、整数jに1を加算し(ステップS15)、ステップS11〜S13の処理を繰り返し行う。
When the integer j is not N (N in step S14), the measuring
そして、整数jがNになると(ステップS14のY)、計測装置1は、係数値算出ステップの処理を終了する。
Then, when the integer j becomes N (Y in step S14), the measuring
図14は、図12のステップS4である統計値算出ステップの手順の一例を示すフローチャート図である。 FIG. 14 is a flowchart showing an example of the procedure of the statistical value calculation step which is step S4 of FIG.
図14に示すように、まず、計測装置1は、整数iを1に設定し(ステップS41)、観測点Riが冗長な観測点ではない、すなわち、レーンL1〜LNのいずれかに対応付けられる観測点である場合(ステップS42のN)、観測点Riに対して、複数の作用の平均による最尤値xi_avg(t)を算出する(ステップS43)。
As shown in FIG. 14, the measuring
次に、計測装置1は、観測点Riに対して、複数の作用の分散xi_σ 2(t)を算出する(ステップS44)。 Next, the measuring device 1 calculates the variance x i_σ 2 (t) of a plurality of actions with respect to the observation point R i (step S44).
次に、計測装置1は、観測点Riに対して、複数の作用の中間値xi_medi(t)を算出する(ステップS45)。
Next, the measuring
計測装置1は、観測点Riが冗長な観測点である場合(ステップS42のY)、ステップS43〜S45の処理を行わない。
Measuring
そして、計測装置1は、整数iがMでない場合は(ステップS46のN)、整数iに1を加算し(ステップS47)、ステップS42〜S45の処理を繰り返し行い、整数iがMになると(ステップS46のY)、統計値算出ステップの処理を終了する。
Then, when the integer i is not M (N in step S46), the measuring
図15は、図12のステップS5である異常判断ステップの手順の一例を示すフローチャート図である。 FIG. 15 is a flowchart showing an example of the procedure of the abnormality determination step which is step S5 of FIG.
図15に示すように、まず、計測装置1は、整数iを1に設定し(ステップS50)、観測点Riが冗長な観測点ではない場合(ステップS51のN)、観測点Riに対して、図14のステップS43で算出した最尤値xi_avg(t)から、図14のステップS45で算出した中間値xi_medi(t)を減算した波形xi_avg−medi(t)を算出する(ステップS52)。
As shown in FIG. 15, the measuring
次に、計測装置1は、図14のステップS44で算出した分散xi_σ 2(t)が第1の閾値以下である場合(ステップS53のN)、最尤値xi_avg(t)を、観測点Riの作用xiとする(ステップS54)。 Next, the measuring device 1 observes the maximum likelihood value x i_avg (t) when the variance x i_σ 2 (t) calculated in step S44 of FIG. 14 is equal to or less than the first threshold value (N in step S53). and acting x i of the point R i (step S54).
また、計測装置1は、分散xi_σ 2(t)が第1の閾値よりも大きい場合(ステップS53のY)、ステップS52で算出した波形xi_avg−medi(t)が第2の閾値よりも大きい場合(ステップS55のY)、異常と判断し、サーバー2等に異常を通知する(ステップS56)。
Further, in the
そして、計測装置1は、中間値xi_medi(t)を、観測点Riの作用xiとする(ステップS57)。
Then, the measuring
計測装置1は、分散xi_σ 2(t)が第1の閾値以下の場合(ステップS53のN)、ステップS56の処理を行わずにステップS57の処理を行う。
When the variance x i_σ 2 (t) is equal to or less than the first threshold value (N in step S53), the measuring
また、計測装置1は、観測点Riが冗長な観測点である場合(ステップS51のY)、ステップS52〜S57の処理を行わない。
Further, the measuring
そして、計測装置1は、整数iがMでない場合は(ステップS58のN)、整数iに1を加算し(ステップS59)、ステップS51〜S57の処理を繰り返し行い、整数iがMになると(ステップS58のY)、異常判断ステップの処理を終了する。
Then, when the integer i is not M (N in step S58), the measuring
1−6.計測装置の構成
図16は、本実施形態における計測装置1の構成例を示す図である。図16に示すように、計測装置1は、制御部110と、第1通信部120と、記憶部130と、第2通信部140と、操作部150と、を有している。
1-6. Configuration of the Measuring Device FIG. 16 is a diagram showing a configuration example of the measuring
制御部110は、上部構造7に設置された各センサー23から出力される加速度データに基づいて、上部構造7の変位等を算出し、また、計測システム10あるいは上部構造7の異常を判断する。
The
第1通信部120は、各センサー23から、加速度データを受信する。各センサー23から出力される加速度データは、例えば、デジタル信号である。第1通信部120は、各センサー23から受信した加速度データを制御部110に出力する。
The
記憶部130は、制御部110が計算処理や制御処理を行うためのプログラムやデータ等を記憶するメモリーである。また、記憶部130は、制御部110が所定のアプリケーション機能を実現するためのプログラムやデータ等を記憶している。記憶部130は、例えば、ROM(Read Only Memory)やフラッシュROM、RAM(Random Access Memory)等の各種IC(Integrated Circuit)メモリーやハードディスク、メモリーカードなどの記録媒体等により構成される。
The
記憶部130は、コンピューターにより読み取り可能な装置や媒体である不揮発性の情報記憶装置を含み、各種のプログラムやデータ等は当該情報記憶装置に記憶されていてもよい。情報記憶装置は、光ディスクDVD、CD等の光ディスク、ハードディスクドライブ、或いはカード型メモリーやROM等の各種のメモリー等であってもよい。また、制御部110が通信ネットワーク4を介して各種のプログラムやデータ等を受信して記憶部130に記憶させてもよい。
The
第2通信部140は、通信ネットワーク4を介して、制御部110の計算結果等の情報や、計測システム10あるいは上部構造7が異常であることを示す情報をサーバー2に送信する。
The
操作部150は、ユーザーからの操作データを取得し、制御部110に送信する処理を行う。
The
制御部110は、物理量取得部111と、作用算出部112と、統計値算出部113と、異常判断部114と、変位算出部115と、荷重算出部116と、係数値算出部117と、出力処理部118と、を備えている。
The
物理量取得部111は、観測点R1〜RMを観測するM個のセンサー23による観測情報に基づいて、車両6が上部構造7を移動するときの観測点R1〜RMの物理量として変位u1〜uMを取得する。すなわち、物理量取得部111は、図12における物理量取得ステップの処理を行う。物理量取得部111が取得した変位u1〜uMは、記憶部130に記憶される。
Physical
作用算出部112は、1以上M以下の任意の整数i,jに対して、観測点Riの変位uiが関数yi1〜yiMの値の和に等しいものとして、物理量取得部111が取得した変位u1〜uMのうちのN個の変位を複数組選択し、当該複数組のそれぞれのN個の変位に基づいて、観測点R1〜RMのうちのレーンL1〜LNに対応付けられたN個の観測点の作用をそれぞれ複数算出する。すなわち、作用算出部112は、図12における作用算出ステップの処理を行う。作用算出部112が算出したN個の観測点のそれぞれの複数の作用は、記憶部130に記憶される。
統計値算出部113は、N個の観測点のそれぞれに対して、作用算出部112が算出した複数の作用の統計値を算出する。すなわち、統計値算出部113は、図12における統計値算出ステップの処理を行う。統計値算出部113が算出したN個の観測点のそれぞれの複数の作用の統計値は、記憶部130に記憶される。
The statistical
異常判断部114は、N個の観測点のそれぞれに対して統計値算出部113が算出した統計値に基づいて、異常を判断する。すなわち、異常判断部114は、図12における異常判断ステップの処理を行う。異常判断部114による判断結果の情報は、記憶部130に記憶される。
The
変位算出部115は、N個の観測点のそれぞれに対して、統計値算出部113が算出した統計値に基づいて、変位を算出する。すなわち、変位算出部115は、図13における変位算出ステップの処理を行う。変位算出部115が算出したN個の観測点のそれぞれの変位は、記憶部130に記憶される。
The
荷重算出部116は、変位算出部115が算出したN個の観測点の変位に基づいて、レーンL1〜LNのそれぞれを走行した車両6による荷重を算出する。すなわち、荷重算出部116は、図12における荷重算出ステップの処理を行う。荷重算出部116が算出したレーンL1〜LNのそれぞれを走行した車両6による荷重は、記憶部130に記憶される。
係数値算出部117は、1以上M以下の任意の整数i,jに対して、車両が単独で上部構造7を走行したときの観測点R1〜RMの変位u1〜uMを取得し、変位u1〜uMに基づいて、関数yijの係数aij,bijの値を算出する。すなわち、係数値算出部117は、図12における係数値算出ステップの処理を行う。係数値算出部117が算出した1次係数a11〜aMMの値及び0次係数b11〜bMMの値は、記憶部130に記憶される。
Coefficient
出力処理部118は、荷重算出部116が算出したレーンL1〜LNを走行した車両6の荷重を、第2通信部140を介してサーバー2に出力する処理を行う。すなわち、出力処理部118は、図12における出力ステップの処理を行う。
例えば、制御部110は、操作部150からの操作データに基づいて、未知の車両6による上部構造7の変位等の算出や計測システム10あるいは上部構造7の異常の判断を行う第1モードと、1次係数a11〜aMM及び0次係数b11〜bMMを算出する第2モードとを切り替える。例えば、N個のセンサー23が上部構造7に設置された後、制御部110が第2モードに設定された状態で複数の車両による荷重試験が行われ、荷重試験の終了後、制御部110は第1モードに設定される。
For example, the
本実施形態では、制御部110は、記憶部130に記憶された各種のプログラムを実行するプロセッサーであり、記憶部130に記憶された計測プログラム131を実行することにより、物理量取得部111、作用算出部112、統計値算出部113、異常判断部114、変位算出部115、荷重算出部116、係数値算出部117、出力処理部118の各機能を実現する。換言すれば、計測プログラム131は、図12に示したフローチャートの各手順を、コンピューターである計測装置1に実行させるプログラムである。
In the present embodiment, the
プロセッサーは、例えば各部の機能が個別のハードウェアで実現されてもよいし、或いは各部の機能が一体のハードウェアで実現されてもよい。例えば、プロセッサーはハードウェアを含み、そのハードウェアは、デジタル信号を処理する回路及びアナログ信号を処理する回路の少なくとも一方を含むことができる。プロセッサーは、CPU(Central Processing Unit)、GPU(Graphics Processing Unit)、或いはDSP(Digital Signal Processor)等であってもよい。ただし、制御部110は、ASIC(Application Specific Integrated Circuit)などのカスタムIC(Integrated Circuit)として構成され、各部の機能を実現してもよいし、CPUとASICとによって各部の機能を実現してもよい。
In the processor, for example, the functions of each part may be realized by individual hardware, or the functions of each part may be realized by integrated hardware. For example, a processor includes hardware, which hardware can include at least one of a circuit that processes a digital signal and a circuit that processes an analog signal. The processor may be a CPU (Central Processing Unit), a GPU (Graphics Processing Unit), a DSP (Digital Signal Processor), or the like. However, the
なお、制御部110は荷重算出部116を含まなくてもよい。また、制御部110は異常判断部114を含まなくてもよい。例えば、計測装置1が、レーンL1〜LNに対応付けられたN個の観測点に対してそれぞれ算出した複数の作用の統計値の情報をサーバー2に送信し、サーバー2が当該統計値に基づいて計測システム10あるいは上部構造7の異常を判断してもよい。また、制御部110は係数値算出部117を含まなくてもよい。例えば、サーバー2あるいは他の装置が、1次係数a11〜aMMの値及び0次係数b11〜bMMの値を算出する処理を行い、これらの値を計測装置1の記憶部130に記憶させてもよい。
The
1−7.作用効果
以上に説明した第1実施形態の計測方法では、計測装置1は、観測点R1〜RMを観測するM個のセンサー23による観測情報に基づいて、観測点R1〜RMの物理量として変位u1〜uMを取得する。そして、観測点Rjの作用xjと、作用xjが観測点Riに及ぼす作用との相関を示す関数をyijとしたとき、計測装置1は、変位uiが、関数yi1〜yiMの値の和に等しいものとして、式(7)で与えられる変位u1〜uMのうちのN個の変位を複数組選択し、当該複数組のそれぞれのN個の変位に基づいて、観測点R1〜RMのうちのレーンL1〜LNに対応付けられたN個の観測点の作用をそれぞれ複数算出する。これにより、計測装置1は、移動体である車両6が構造物である上部構造7の観測点Rjを通過する時の作用xjを他の作用から分離して複数算出することができる。例えば、複数の車両6が複数のレーンを並走する場合でも、計測装置1は、各レーンを走行する車両6による各観測点の作用が他の観測点に及ぼす作用の影響を排除して、各レーンを走行する車両6による各観測点の作用を複数算出することができる。
1-7. In the measurement method of the first embodiment described above advantageous effects, the measuring
さらに、計測装置1は、レーンL1〜LNに対応付けられたN個の観測点のそれぞれに対して、算出した複数の作用の統計値を算出する。計測システムや構造物である上部構造7が正常であれば、算出された複数の作用は同程度の大きさであるが、計測システム10あるいは上部構造7が異常であれば、算出された複数の作用の大きさがばらつくことになる。そして、複数の作用の大きさのばらつきは、複数の作用の統計値によって判断可能である。したがって、第1実施形態の計測方法によれば、計測装置1は、計測システム10の異常や上部構造7の異常を判断するために利用可能な指標としての統計値を算出することができる。
Furthermore, the measuring
また、第1実施形態の計測方法では、計測装置1は、レーンL1〜LNに対応付けられたN個の観測点のそれぞれに対して、算出した複数の作用の統計値に基づいて、変位をし、算出したN個の観測点の変位に基づいて、レーンL1〜LNのそれぞれを走行した車両6による荷重を算出する。したがって、第1実施形態の計測方法によれば、計測装置1は、他の作用から分離して算出された各観測点の複数の作用の統計値に基づいて、車両6の走行による各観測点の変位や荷重を精度良く算出することができる。例えば、複数の車両6が複数のレーンを並走する場合でも、計測装置1は、レーンLjを移動した車両6によるレーンLjの変位や荷重を精度良く算出することができる。この変位や荷重の情報により、例えば、計測装置1あるいはサーバー2は、過積載の車両の監視等の処理を精度良く行うことができる。
Further, in the measurement method of the first embodiment, the measuring
また、第1実施形態の計測方法によれば、計測装置1は、上部構造7を通過する車両6の車軸重量による上部構造7の変位や車両6の荷重を算出することができるため、上部構造7の損傷を予測するための橋梁5の維持管理のために十分な情報を提供することができる。
Further, according to the measurement method of the first embodiment, the measuring
2.第2実施形態
第1実施形態の計測方法では、物理量取得ステップにおいて、計測装置1は、観測点R1〜RNの物理量として観測点R1〜RNの変位u1〜uNを取得する。これに対して、第2実施形態の計測方法では、物理量取得ステップにおいて、計測装置1は、観測点R1〜RNの物理量として観測点R1〜RNの車両6による荷重w1〜wNを取得する。以下、第2実施形態について、第1実施形態と同様の構成要素には同じ符号を付して第1実施形態と重複する説明を省略又は簡略し、主に第1実施形態と異なる内容について説明する。
2. The measuring method of the second embodiment first embodiment, the physical quantity obtaining step, measuring
式(33)に示すように、観測点Riの車両6による荷重fiは、前述の関数yi1〜yiMの値の和に等しいものとする。
As shown in equation (33), the load f i by the
このとき、前述の式(1)及び式(33)より、観測点R1〜RMの車両6による荷重f1〜fMを要素とする荷重ベクトルfは、式(34)のように表される。
Table as in this case, from the above equations (1) and (33), load vector f to a
式(34)において、ベクトルYの各要素ykは式(35)のように定義される。kは1以上M以下の任意の整数である。 In equation (34), each element y k of the vector Y is defined as in equation (35). k is an arbitrary integer of 1 or more and M or less.
実際に観測される観測点R1〜RMのそれぞれの車両6による荷重w1〜wMを要素とする荷重ベクトルwが、荷重ベクトルfと等しいものとすると、式(36)が得られる。
Indeed load vector w to the
式(36)において、M個の観測点R1〜RMのうちのM−N個の観測点は、レーンL1〜LNのいずれにも対応付けられないので、車両6がレーンL1〜LNのいずれかを走行しても、当該M−N個の観測点には車両6の走行による作用は生じないものと考える。すなわち、当該M−N個の観測点の車両6による作用はすべてゼロであるものと考える。そうすると、レーンL1〜LNに対応付けられるN個の観測点の未知の作用を算出するために、式(16)で得られるM個の荷重w1〜wMから任意のN個の荷重を選択して連立方程式を立てることができる。M個の荷重w1〜wMから選択可能なN個の荷重の組み合わせ数Cは、前述の式(8)で表される。
In the formula (36), M-N pieces of observation points of the
計測装置1は、複数組のN個の荷重の連立方程式を選択し、各組の連立方程式を解いて、車両6がレーンL1〜LNをそれぞれ単独で走行した時のレーンL1〜LNに対応するN個の観測点の作用をそれぞれ複数算出する。
Measuring
一例として、前述の図5及び図6に示した配置例の場合、式(36)から選択される3組の連立方程式を解いて得られる作用x1,x3をまとめると、3組の作用x1は式(37)のようになり、3組の作用x3は式(38)のようになる。式(37)及び式(38)は、前述の式(24)及び式(25)のそれぞれに対して変位u1(t),u3(t)を荷重w1(t),w3(t)に置き換えた式である。 As an example, in the case of the arrangement example shown in FIGS. 5 and 6 described above, the actions x 1 and x 3 obtained by solving the three sets of simultaneous equations selected from the equation (36) can be summarized as three sets of actions. x 1 is as in equation (37), and the three sets of actions x 3 are as in equation (38). Equations (37) and (38) load displacements u 1 (t) and u 3 (t) with respect to the above equations (24) and (25), respectively, w 1 (t), w 3 ( It is an equation replaced with t).
図17は、第2実施形態の計測方法の手順の一例を示すフローチャート図である。本実施形態では、計測装置1が図17に示す手順を実行する。
FIG. 17 is a flowchart showing an example of the procedure of the measurement method of the second embodiment. In this embodiment, the measuring
図17に示すように、まず、計測装置1は、車両が単独で上部構造7を走行したときの観測点R1〜RMの車両による荷重w1〜wMを取得し、荷重w1〜wMに基づいて、関数yijの係数aij,bijの値を算出する(ステップS101)。i,jは1以上M以下の任意の整数である。車両は、未知の移動体である車両6とは異なる既知の移動体である。このステップS101は係数値算出ステップである。
As shown in FIG. 17, first, the
次に、計測装置1は、観測点R1〜RMを観測するM個のセンサー23による観測情報に基づいて、車両6が上部構造7を移動するときの観測点R1〜RMの物理量として車両6による荷重w1〜wMを取得する(ステップS102)。前述の通り、M個のセンサー23はそれぞれ加速度センサーであり、M個のセンサー23による観測情報は、観測点R1〜RMに生じた加速度の検出情報である。そして、この加速度は、第1方向であるX方向と第2方向であるY方向とそれぞれ交差する第3方向の加速度である。計測装置1は、M個のセンサー23がそれぞれ検出した第3方向の加速度に基づいて、式(36)に含まれる荷重ベクトルwを算出する。したがって、物理量取得ステップで計測装置1が取得する観測点R1〜RMの物理量としての荷重w1〜wMは、X方向及びY方向とそれぞれ交差する第3方向の荷重、例えば、X方向及びY方向とそれぞれ直交する第3方向の荷重である。このステップS102は物理量取得ステップである。
Next, the measuring
次に、計測装置1は、観測点Riの車両6による荷重wiが関数yi1〜yiMの値の和に等しいものとして、ステップS102で取得した荷重w1〜wMのうちのN個の荷重を複数組選択し、当該複数組のそれぞれのN個の荷重に基づいて、観測点R1〜RMのうちのレーンL1〜LNに対応付けられたN個の観測点の作用をそれぞれ複数算出する(ステップS103)。iは1以上M以下の任意の整数である。具体的には、計測装置1は、式(36)で得られるM個の荷重w1〜wMからN個の荷重を複数組選択し、当該複数組の荷重の連立方程式を解いて、レーンL1〜LNに対応付けられたN個の観測点の作用をそれぞれ複数算出する。このステップS103は作用算出ステップである。
Next, the measuring
次に、計測装置1は、N個の観測点のそれぞれに対して、ステップS103で算出した複数の作用の統計値を算出する(ステップS104)。統計値は、例えば、最尤値、標準偏差、分散又は中間値である。このステップS104は統計値算出ステップである。
Next, the measuring
次に、計測装置1は、N個の観測点のそれぞれに対してステップS104で算出した統計値に基づいて、異常を判断する(ステップS105)。このステップS105は異常判断ステップである。
Next, the measuring
次に、計測装置1は、N個の観測点のそれぞれに対して、ステップS105で算出した統計値に基づいて、車両6による荷重を算出する(ステップS106)。例えば、レーンL1〜LNに対応付けられたN個の観測点のうちの観測点Rkに対して、ステップS105で算出した複数の作用の平均による最尤値や複数の作用の中間値を作用xkとして、式(34)の右辺において、作用x1〜xNのうちの作用xkを除くすべての作用をゼロとして算出される荷重fkを観測点Rkの車両6による荷重とする。このステップS106は荷重算出ステップである。
Next, the measuring
次に、計測装置1は、ステップS106で算出したN個の観測点の車両6による荷重に基づいて、レーンL1〜LNのそれぞれの変位を算出する(ステップS107)。1以上N以下の各整数jに対して、レーンLjに対応付けられた観測点の車両6による荷重とレーンLjの変位との間には相関があるので、あらかじめ車両による荷重試験において、この相関式の係数を算出しておく。計測装置1は、当該相関式に、レーンLjに対応付けられた観測点の車両6による荷重を代入して、レーンLjの変位を算出することができる。このステップS107は変位算出ステップである。
Next, the measuring
次に、計測装置1は、ステップS107で算出したレーンL1〜LNの変位をサーバー2に出力する(ステップS108)。このステップS108は、出力ステップである。
Next, the measuring
計測装置1は、計測を終了するまで(ステップS109のN)、ステップS102〜S108の処理を繰り返し行う。
The measuring
第2実施形態における計測装置1の構成は、図16と同様であるため、その図示を省略する。第1実施形態と同様、計測装置1は、制御部110と、第1通信部120と、記憶部130と、第2通信部140と、操作部150と、を有している。
Since the configuration of the measuring
第1通信部120、記憶部130、第2通信部140及び操作部150がそれぞれ行う処理は、第1実施形態と同様であるため、その説明を省略する。
Since the processes performed by the
制御部110は、上部構造7に設置された各センサー23から出力される加速度データに基づいて、車両6による荷重等を算出し、また、計測システム10あるいは上部構造7の異常を判断する。第1実施形態と同様、制御部110は、物理量取得部111と、作用算出部112と、統計値算出部113と、異常判断部114と、変位算出部115と、荷重算出部116と、係数値算出部117と、出力処理部118と、を備えている。
The
物理量取得部111は、観測点R1〜RMを観測するM個のセンサー23による観測情報に基づいて、車両6が上部構造7を移動するときの観測点R1〜RMの物理量として車両6による荷重w1〜wMを取得する。すなわち、物理量取得部111は、図17における物理量取得ステップの処理を行う。物理量取得部111が取得した荷重w1〜wMは、記憶部130に記憶される。
Physical
作用算出部112は、1以上M以下の任意の整数i,jに対して、観測点Riの車両6による荷重wiが関数yi1〜yiMの値の和に等しいものとして、物理量取得部111が取得した荷重w1〜wMのうちのN個の荷重を複数組選択し、当該複数組のそれぞれのN個の荷重に基づいて、観測点R1〜RMのうちのレーンL1〜LNに対応付けられたN個の観測点の作用をそれぞれ複数算出する。すなわち、作用算出部112は、図17における作用算出ステップの処理を行う。作用算出部112が算出したN個の観測点のそれぞれの複数の作用は、記憶部130に記憶される。
The action calculation unit 112 acquires the physical quantity assuming that the load w i by the
統計値算出部113は、N個の観測点のそれぞれに対して、作用算出部112が算出した複数の作用の統計値を算出する。すなわち、統計値算出部113は、図17における統計値算出ステップの処理を行う。統計値算出部113が算出したN個の観測点のそれぞれの複数の作用の統計値は、記憶部130に記憶される。
The statistical
異常判断部114は、N個の観測点のそれぞれに対して統計値算出部113が算出した統計値に基づいて、異常を判断する。すなわち、異常判断部114は、図17における異常判断ステップの処理を行う。異常判断部114による判断結果の情報は、記憶部130に記憶される。
The
変位算出部115は、荷重算出部116が算出したN個の観測点の車両6による荷重に基づいて、レーンL1〜LNのそれぞれの変位を算出する。すなわち、変位算出部115は、図17における変位算出ステップの処理を行う。変位算出部115が算出したレーンL1〜LNのそれぞれの変位は、記憶部130に記憶される。
荷重算出部116は、N個の観測点のそれぞれに対して、統計値算出部113が算出した統計値に基づいて、車両6による荷重を算出する。すなわち、荷重算出部116は、図17における荷重算出ステップの処理を行う。荷重算出部116が算出したN個の観測点のそれぞれの車両6による荷重は、記憶部130に記憶される。
The
係数値算出部117は、1以上M以下の任意の整数i,jに対して、車両が単独で上部構造7を走行したときの観測点R1〜RMの車両による荷重w1〜wMを取得し、荷重w1〜wMに基づいて、関数yijの係数aij,bijの値を算出する。すなわち、係数値算出部117は、図17における係数値算出ステップの処理を行う。係数値算出部117が算出した1次係数a11〜aMMの値及び0次係数b11〜bMMの値は、記憶部130に記憶される。
Coefficient
出力処理部118は、変位算出部115が算出したレーンL1〜LNの変位を、第2通信部140を介してサーバー2に出力する処理を行う。すなわち、出力処理部118は、図17における出力ステップの処理を行う。
例えば、制御部110は、操作部150からの操作データに基づいて、未知の車両6による荷重等の算出や計測システム10あるいは上部構造7の異常の判断を行う第1モードと、1次係数a11〜aMM及び0次係数b11〜bMMを算出する第2モードとを切り替える。例えば、N個のセンサー23が上部構造7に設置された後、制御部110が第2モードに設定された状態で複数の車両による荷重試験が行われ、荷重試験の終了後、制御部110は第1モードに設定される。
For example, the
第1実施形態と同様、制御部110は、記憶部130に記憶された各種のプログラムを実行するプロセッサーであり、記憶部130に記憶された計測プログラム131を実行することにより、物理量取得部111、作用算出部112、統計値算出部113、異常判断部114、変位算出部115、荷重算出部116、係数値算出部117、出力処理部118の各機能を実現する。換言すれば、計測プログラム131は、図17に示したフローチャートの各手順を、コンピューターである計測装置1に実行させるプログラムである。ただし、制御部110は、ASICなどのカスタムICとして構成され、各部の機能を実現してもよいし、CPUとASICとによって各部の機能を実現してもよい。
Similar to the first embodiment, the
なお、制御部110は変位算出部115を含まなくてもよい。また、制御部110は異常判断部114を含まなくてもよい。例えば、計測装置1が、レーンL1〜LNに対応付けられたN個の観測点に対してそれぞれ算出した複数の作用の統計値の情報をサーバー2に送信し、サーバー2が当該統計値に基づいて計測システム10あるいは上部構造7の異常を判断してもよい。また、制御部110は係数値算出部117を含まなくてもよい。例えば、サーバー2あるいは他の装置が、1次係数a11〜aMMの値及び0次係数b11〜bMMの値を算出する処理を行い、これらの値を計測装置1の記憶部130に記憶させてもよい。
The
以上に説明した第2実施形態の計測方法では、計測装置1は、観測点R1〜RMを観測するM個のセンサー23による観測情報に基づいて、観測点R1〜RMの物理量として車両6による荷重w1〜wMを取得する。そして、観測点Rjの作用xjと、作用xjが観測点Riに及ぼす作用との相関を示す関数をyijとしたとき、計測装置1は、荷重wiが、関数yi1〜yiMの値の和に等しいものとして、式(36)で与えられる荷重w1〜wMのうちのN個の荷重を複数組選択し、当該複数組のそれぞれのN個の荷重に基づいて、観測点R1〜RMのうちのレーンL1〜LNに対応付けられたN個の観測点の作用をそれぞれ複数算出する。これにより、計測装置1は、移動体である車両6が構造物である上部構造7の観測点Rjを通過する時の作用xjを他の作用から分離して複数算出することができる。例えば、複数の車両6が複数のレーンを並走する場合でも、計測装置1は、各レーンを走行する車両6による各観測点の作用が他の観測点に及ぼす作用の影響を排除して、各レーンを走行する車両6による各観測点の作用を複数算出することができる。
The measuring method of the second embodiment described above, the measuring
さらに、計測装置1は、レーンL1〜LNに対応付けられたN個の観測点のそれぞれに対して、算出した複数の作用の統計値を算出する。計測システムや構造物である上部構造7が正常であれば、算出された複数の作用は同程度の大きさであるが、計測システム10あるいは上部構造7が異常であれば、算出された複数の作用の大きさがばらつくことになる。そして、複数の作用の大きさのばらつきは、複数の作用の統計値によって判断可能である。したがって、第2実施形態の計測方法によれば、計測装置1は、計測システム10の異常や上部構造7の異常を判断するために利用可能な指標としての統計値を算出することができる。
Furthermore, the measuring
また、第2実施形態の計測方法では、計測装置1は、レーンL1〜LNに対応付けられたN個の観測点のそれぞれに対して、算出した複数の作用の統計値に基づいて、車両6による荷重をし、算出したN個の観測点の車両6による荷重に基づいて、レーンL1〜LNのそれぞれの変位を算出する。したがって、第2実施形態の計測方法によれば、計測装置1は、他の作用から分離して算出された各観測点の複数の作用の統計値に基づいて、車両6の走行による各観測点の荷重や変位を精度良く算出することができる。例えば、複数の車両6が複数のレーンを並走する場合でも、計測装置1は、レーンLjを移動した車両6によるレーンLjの荷重や変位を精度良く算出することができる。この荷重や変位の情報により、例えば、計測装置1あるいはサーバー2は、過積載の車両の監視等の処理を精度良く行うことができる。
Further, in the measurement method of the second embodiment, the measuring
3.変形例
本発明は本実施形態に限定されず、本発明の要旨の範囲内で種々の変形実施が可能である。
3. 3. Modifications The present invention is not limited to the present embodiment, and various modifications can be carried out within the scope of the gist of the present invention.
上記の各実施形態では、車両6がレーンLjを走行した時の観測点Rjの作用xjと、観測点Rjの作用xjが観測点Riに及ぼす作用との相関を示す関数yijは、式(1)に示す1次多項式関数であるものとしたが、相関が直線でない場合は、関数yijは、式(39)に示すようなm次多項式関数であってもよい。
In the above embodiments, function indicating the action x j of the observation point R j when the
また、上記の各実施形態では、観測点R1〜RMを観測する観測装置は、それぞれ加速度センサーであるが、これに限られず、例えば、接触式変位計、リング式変位計、レーザー変位計、感圧センサー、画像処理による変位計測機器又は光ファイバーによる変位計測機器であってもよい。観測装置と観測点が1対1に対応している必要はなく、1つの観測装置が観測点R1〜RMの一部又は全部を観測してもよい。 In the embodiments described above, the observation apparatus for observing the observation point R 1 to R M is respectively the acceleration sensor is not limited to this, for example, the contact displacement sensor, a ring-type displacement meter, laser displacement meter , A pressure sensor, a displacement measuring device by image processing, or a displacement measuring device by an optical fiber. It is not necessary to the observation apparatus and observation point corresponds to one-to-one, may be a single observation device observes a part or the whole of the observation point R 1 to R M.
接触式変位計、リング式変位計、レーザー変位計、画像処理による変位計測機器、光ファイバーによる変位計測機器は、車両6の観測点R1〜RMの作用に対する応答として変位を計測する。計測装置1は、観測点R1〜RMの変位に基づいて、観測点R1〜RMの物理量としての変位又は車両6による荷重を算出する。感圧センサーは、車両6の観測点R1〜RMへの作用に対する応答として応力変化を検出する。計測装置1は、観測点R1〜RMの応力変化に基づいて、観測点R1〜RMの物理量としての変位又は車両6による荷重を算出する。
Contact displacement meter, the ring-type displacement meter, laser displacement meter, a displacement measuring device by image processing, a displacement measuring apparatus according to the optical fiber, measures the displacement in response to the action of observation points R 1 to R M of the
また、上記の各実施形態では、車両6がレーンL1〜LNを走行する方向はすべて同じであるが、レーンL1〜LNのうちの少なくとも1つのレーンと他のレーンとは車両6の走行方向が異なっていてもよい。
In the embodiments described above, all directions in which the
また、上記の各実施形態では、各センサー23は、それぞれ上部構造7の主桁Gに設けられているが、上部構造7の表面や内部、床板Fの下面、橋脚8a等に設けられていてもよい。また、上記の各実施形態では、橋梁5として道路橋を例に挙げたが、これに限られず、例えば、橋梁5は鉄道橋であってもよい。また、上記の各実施形態では、構造物として橋梁の上部構造を例に挙げたが、これに限られず、構造物は移動体の移動によって変形するものであればよい。
Further, in each of the above embodiments, each
上述した実施形態および変形例は一例であって、これらに限定されるわけではない。例えば、各実施形態および各変形例を適宜組み合わせることも可能である。 The above-described embodiments and modifications are merely examples, and the present invention is not limited thereto. For example, each embodiment and each modification can be combined as appropriate.
本発明は、実施の形態で説明した構成と実質的に同一の構成、例えば、機能、方法及び結果が同一の構成、あるいは目的及び効果が同一の構成を含む。また、本発明は、実施の形態で説明した構成の本質的でない部分を置き換えた構成を含む。また、本発明は、実施の形態で説明した構成と同一の作用効果を奏する構成又は同一の目的を達成することができる構成を含む。また、本発明は、実施の形態で説明した構成に公知技術を付加した構成を含む。 The present invention includes a configuration substantially the same as the configuration described in the embodiment, for example, a configuration having the same function, method and result, or a configuration having the same purpose and effect. The present invention also includes a configuration in which a non-essential part of the configuration described in the embodiment is replaced. The present invention also includes a configuration that exhibits the same effects as the configuration described in the embodiment or a configuration that can achieve the same object. Further, the present invention includes a configuration in which a known technique is added to the configuration described in the embodiment.
1…計測装置、2…サーバー、4…通信ネットワーク、5…橋梁、6…車両、7…上部構造、7a…橋床、7b…支承、G…主桁、F…床板、8…下部構造、8a…橋脚、8b…橋台、10…計測システム、23…センサー、110…制御部、111…物理量取得部、112…作用算出部、113…統計値算出部、114…異常判断部、115…変位算出部、116…荷重算出部、117…係数値算出部、118…出力処理部、120…第1通信部、130…記憶部、131…計測プログラム、140…第2通信部、150…操作部 1 ... Measuring device, 2 ... Server, 4 ... Communication network, 5 ... Bridge, 6 ... Vehicle, 7 ... Superstructure, 7a ... Abutment, 7b ... Bearing, G ... Main girder, F ... Floor board, 8 ... Substructure, 8a ... pier, 8b ... abutment, 10 ... measurement system, 23 ... sensor, 110 ... control unit, 111 ... physical quantity acquisition unit, 112 ... action calculation unit, 113 ... statistical value calculation unit, 114 ... abnormality judgment unit, 115 ... displacement Calculation unit, 116 ... Load calculation unit, 117 ... Coefficient value calculation unit, 118 ... Output processing unit, 120 ... First communication unit, 130 ... Storage unit, 131 ... Measurement program, 140 ... Second communication unit, 150 ... Operation unit
Claims (15)
1以上M以下の任意の整数i及び1以上M以下の任意の整数jに対して、前記第jの観測点の作用xjと、前記作用xjが前記第iの観測点に及ぼす作用との相関を示す関数をyijとしたとき、前記物理量取得ステップで取得した前記第iの観測点の前記物理量が、関数yi1〜yiMの値の和に等しいものとして、前記物理量取得ステップで取得した前記第1〜第Mの観測点の前記物理量のうちのN個の前記物理量を複数組選択し、前記複数組のそれぞれのN個の前記物理量に基づいて、前記第1〜第Mの観測点のうちの前記第1〜第Nの経路に対応付けられたN個の観測点の作用をそれぞれ複数算出する作用算出ステップと、
前記N個の観測点のそれぞれに対して、前記作用算出ステップで算出した複数の前記作用の統計値を算出する統計値算出ステップと、を含む、計測方法。 For an integer N of 2 or more and an integer M larger than the integer N, the moving body is arranged along a second direction intersecting the first direction in which the moving body travels in any of the first to Nth paths of the structure. A physical quantity acquisition step of acquiring the physical quantity of the first to Mth observation points based on the observation information by at least one observation device for observing the first to Mth observation points of the structure.
For one or more M or less arbitrary integer i and 1 to M any integer j, the effect x j of the observation point of the first j, and action of the working x j is on the observation point of the i-th Assuming that the physical quantity of the i-th observation point acquired in the physical quantity acquisition step is equal to the sum of the values of the functions y i1 to y iM , when the function showing the correlation of is y ij, the physical quantity acquisition step A plurality of sets of N physical quantities among the acquired physical quantities of the first to M observation points are selected, and based on each of the N physical quantities of the plurality of sets, the first to Mth. An action calculation step for calculating a plurality of actions of each of the N observation points associated with the first to Nth paths among the observation points, and an action calculation step.
A measurement method including a statistical value calculation step for calculating a plurality of statistical values of the action calculated in the action calculation step for each of the N observation points.
前記上部構造は、隣り合う橋台と橋脚、隣り合う2つの橋台、又は、隣り合う2つの橋脚のいずれか1つに渡された構造であり、
前記上部構造の両端部は、前記隣り合う橋台と橋脚の位置、前記隣り合う2つの橋台の位置、又は、前記隣り合う2つの橋脚の位置にあり、
前記橋梁は、道路橋又は鉄道橋である、請求項1乃至10のいずれか一項に記載の計測方法。 The structure is a superstructure of a bridge and
The superstructure is a structure passed to any one of adjacent abutments and piers, two adjacent abutments, or two adjacent piers.
Both ends of the superstructure are located at the positions of the adjacent abutments and piers, the positions of the two adjacent abutments, or the positions of the two adjacent piers.
The measurement method according to any one of claims 1 to 10, wherein the bridge is a road bridge or a railway bridge.
1以上M以下の任意の整数i及び1以上M以下の任意の整数jに対して、前記第jの観測点の作用xjと、前記作用xjが前記第iの観測点に及ぼす作用との相関を示す関数をyijとしたとき、前記物理量取得部が取得した前記第iの観測点の前記物理量が、関数yi1〜yiMの値の和に等しいものとして、前記物理量取得部が取得した前記第1〜第Mの観測点の前記物理量のうちのN個の前記物理量を複数組選択し、前記複数組のそれぞれのN個の前記物理量に基づいて、前記第1〜第Mの観測点のうちの前記第1〜第Nの経路に対応付けられたN個の観測点の作用をそれぞれ複数算出する作用算出部と、
前記N個の観測点のそれぞれに対して、前記作用算出部が算出した複数の前記作用の統計値を算出する統計値算出部と、を含む、計測装置。 For an integer N of 2 or more and an integer M larger than the integer N, the moving body is arranged along a second direction intersecting the first direction in which the moving body travels in any of the first to Nth paths of the structure. A physical quantity acquisition unit that acquires the physical quantity of the first to Mth observation points based on the observation information by at least one observation device that observes the first to Mth observation points of the structure.
For one or more M or less arbitrary integer i and 1 to M any integer j, the effect x j of the observation point of the first j, and action of the working x j is on the observation point of the i-th when the function indicating the correlation of the y ij, the physical quantity of the observation point of the i-th said physical quantity acquisition unit has acquired, as equal to the sum of the values of the function y i1 ~y iM, said physical quantity obtaining unit A plurality of sets of N physical quantities among the acquired physical quantities of the first to M observation points are selected, and based on each of the N physical quantities of the plurality of sets, the first to Mth. An action calculation unit that calculates a plurality of actions of each of the N observation points associated with the first to Nth paths among the observation points, and an action calculation unit.
A measuring device including a statistical value calculation unit for calculating a plurality of statistical values of the action calculated by the action calculation unit for each of the N observation points.
前記観測装置と、を備えた、計測システム。 The measuring device according to claim 13 and
A measurement system including the observation device.
1以上M以下の任意の整数i及び1以上M以下の任意の整数jに対して、前記第jの観測点の作用xjと、前記作用xjが前記第iの観測点に及ぼす作用との相関を示す関数をyijとしたとき、前記物理量取得ステップで取得した前記第iの観測点の前記物理量が、関数yi1〜yiMの値の和に等しいものとして、前記物理量取得ステップで取得した前記第1〜第Mの観測点の前記物理量のうちのN個の前記物理量を複数組選択し、前記複数組のそれぞれのN個の前記物理量に基づいて、前記第1〜第Mの観測点のうちの前記第1〜第Nの経路に対応付けられたN個の観測点の作用をそれぞれ複数算出する作用算出ステップと、
前記N個の観測点のそれぞれに対して、前記作用算出ステップで算出した複数の前記作用の統計値を算出する統計値算出ステップと、をコンピューターに実行させる、計測プログラム。 For an integer N of 2 or more and an integer M larger than the integer N, the moving body is arranged along a second direction intersecting the first direction in which the moving body travels in any of the first to Nth paths of the structure. A physical quantity acquisition step of acquiring the physical quantity of the first to Mth observation points based on the observation information by at least one observation device for observing the first to Mth observation points of the structure.
For one or more M or less arbitrary integer i and 1 to M any integer j, the effect x j of the observation point of the first j, and action of the working x j is on the observation point of the i-th Assuming that the physical quantity of the i-th observation point acquired in the physical quantity acquisition step is equal to the sum of the values of the functions y i1 to y iM , when the function showing the correlation of is y ij, the physical quantity acquisition step A plurality of sets of N physical quantities among the acquired physical quantities of the first to M observation points are selected, and based on each of the N physical quantities of the plurality of sets, the first to Mth. An action calculation step for calculating a plurality of actions of each of the N observation points associated with the first to Nth paths among the observation points, and an action calculation step.
A measurement program for causing a computer to execute a statistical value calculation step for calculating a plurality of statistical values of the action calculated in the action calculation step for each of the N observation points.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2020047304A JP2021147825A (en) | 2020-03-18 | 2020-03-18 | Measurement method, measurement apparatus, measurement system, and measurement program |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2020047304A JP2021147825A (en) | 2020-03-18 | 2020-03-18 | Measurement method, measurement apparatus, measurement system, and measurement program |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2021147825A true JP2021147825A (en) | 2021-09-27 |
Family
ID=77847844
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2020047304A Pending JP2021147825A (en) | 2020-03-18 | 2020-03-18 | Measurement method, measurement apparatus, measurement system, and measurement program |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2021147825A (en) |
-
2020
- 2020-03-18 JP JP2020047304A patent/JP2021147825A/en active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP7375637B2 (en) | Measurement method, measurement device, measurement system and measurement program | |
US11881102B2 (en) | Measurement method, measurement device, measurement system, and measurement program | |
US11761811B2 (en) | Measurement method, measurement device, measurement system, and measurement program | |
US11713993B2 (en) | Bridge displacement measurement method | |
US11714021B2 (en) | Measurement method, measurement device, measurement system, and measurement program | |
CN111521247A (en) | Measuring device, measuring system and measuring method | |
US11408761B2 (en) | Measurement method, measurement device, measurement system, and measurement program | |
JP2021147825A (en) | Measurement method, measurement apparatus, measurement system, and measurement program | |
JP2021147820A (en) | Measurement method, measurement apparatus, measurement system, and measurement program | |
JP6050415B2 (en) | Road surface smoothness evaluation method | |
JP2021148526A (en) | Measuring method, measuring device, measuring system, and measuring program | |
JP2021147818A (en) | Measurement method, measurement apparatus, measurement system, and measurement program | |
JP2021148525A (en) | Measuring method, measuring device, measuring system, and measuring program | |
JP2021147826A (en) | Measurement method, measurement apparatus, measurement system, and measurement program | |
JP2021148527A (en) | Measuring method, measuring device, measuring system, and measuring program | |
JP2021148538A (en) | Measuring method, measuring device, measuring system, and measuring program | |
US20230341289A1 (en) | Measurement Method, Measurement Device, Measurement System, And Non-Transitory Computer-Readable Storage Medium Storing Measurement Program | |
US11921012B2 (en) | Abnormality determination for bridge superstructure based on acceleration data | |
JP2021148535A (en) | Measuring method, measuring device, measuring system, and measuring program | |
US20230175918A1 (en) | Measurement Method, Measurement Device, Measurement System, And Measurement Program | |
US20230003575A1 (en) | Measurement Method, Measurement Device, Measurement System, And Measurement Program | |
US20230019808A1 (en) | Measurement Method, Measurement Device, Measurement System, And Measurement Program | |
JP2023084856A (en) | Measurement method, measurement device, measurement system, and measurement program |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
RD07 | Notification of extinguishment of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7427 Effective date: 20210719 |
|
RD04 | Notification of resignation of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7424 Effective date: 20210914 |
|
RD03 | Notification of appointment of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7423 Effective date: 20211101 |