[go: up one dir, main page]

JP7739245B2 - 構造物評価システム、構造物評価装置、構造物評価方法及びコンピュータプログラム - Google Patents

構造物評価システム、構造物評価装置、構造物評価方法及びコンピュータプログラム

Info

Publication number
JP7739245B2
JP7739245B2 JP2022148052A JP2022148052A JP7739245B2 JP 7739245 B2 JP7739245 B2 JP 7739245B2 JP 2022148052 A JP2022148052 A JP 2022148052A JP 2022148052 A JP2022148052 A JP 2022148052A JP 7739245 B2 JP7739245 B2 JP 7739245B2
Authority
JP
Japan
Prior art keywords
temperature
unit
pavement
correction
elastic wave
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2022148052A
Other languages
English (en)
Other versions
JP2024043068A (ja
Inventor
隆 碓井
英文 高峯
一雄 渡部
哲也 釘宮
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Toshiba Corp
Original Assignee
Toshiba Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Toshiba Corp filed Critical Toshiba Corp
Priority to JP2022148052A priority Critical patent/JP7739245B2/ja
Priority to CN202310129747.6A priority patent/CN117723627A/zh
Priority to US18/173,235 priority patent/US20240094166A1/en
Priority to EP23158754.4A priority patent/EP4339578A1/en
Publication of JP2024043068A publication Critical patent/JP2024043068A/ja
Application granted granted Critical
Publication of JP7739245B2 publication Critical patent/JP7739245B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M5/00Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings
    • G01M5/0066Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings by exciting or detecting vibration or acceleration
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M5/00Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings
    • G01M5/0033Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings by determining damage, crack or wear
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/04Analysing solids
    • G01N29/045Analysing solids by imparting shocks to the workpiece and detecting the vibrations or the acoustic waves caused by the shocks
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/04Analysing solids
    • G01N29/06Visualisation of the interior, e.g. acoustic microscopy
    • G01N29/0654Imaging
    • G01N29/069Defect imaging, localisation and sizing using, e.g. time of flight diffraction [TOFD], synthetic aperture focusing technique [SAFT], Amplituden-Laufzeit-Ortskurven [ALOK] technique
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/04Analysing solids
    • G01N29/07Analysing solids by measuring propagation velocity or propagation time of acoustic waves
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/14Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object using acoustic emission techniques
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/22Details, e.g. general constructional or apparatus details
    • G01N29/32Arrangements for suppressing undesired influences, e.g. temperature or pressure variations, compensating for signal noise
    • G01N29/326Arrangements for suppressing undesired influences, e.g. temperature or pressure variations, compensating for signal noise compensating for temperature variations
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/44Processing the detected response signal, e.g. electronic circuits specially adapted therefor
    • G01N29/4409Processing the detected response signal, e.g. electronic circuits specially adapted therefor by comparison
    • G01N29/4427Processing the detected response signal, e.g. electronic circuits specially adapted therefor by comparison with stored values, e.g. threshold values
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S5/00Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/18Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using ultrasonic, sonic, or infrasonic waves
    • G01S5/22Position of source determined by co-ordinating a plurality of position lines defined by path-difference measurements
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00Indexing codes associated with group G01N29/00
    • G01N2291/02Indexing codes associated with the analysed material
    • G01N2291/023Solids
    • G01N2291/0232Glass, ceramics, concrete or stone
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00Indexing codes associated with group G01N29/00
    • G01N2291/02Indexing codes associated with the analysed material
    • G01N2291/025Change of phase or condition
    • G01N2291/0258Structural degradation, e.g. fatigue of composites, ageing of oils
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00Indexing codes associated with group G01N29/00
    • G01N2291/02Indexing codes associated with the analysed material
    • G01N2291/028Material parameters
    • G01N2291/02818Density, viscosity
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00Indexing codes associated with group G01N29/00
    • G01N2291/02Indexing codes associated with the analysed material
    • G01N2291/028Material parameters
    • G01N2291/02827Elastic parameters, strength or force
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00Indexing codes associated with group G01N29/00
    • G01N2291/26Scanned objects
    • G01N2291/263Surfaces

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Pathology (AREA)
  • Biochemistry (AREA)
  • Immunology (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Signal Processing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Description

本発明の実施形態は、構造物評価システム、構造物評価装置、構造物評価方法及びコンピュータプログラムに関する。
橋梁のような構造物の表面にセンサを設置することで、構造物内部で発生する弾性波を検出することができる。さらに、複数のセンサを構造物の表面に設置することで、各センサで検出した弾性波の到達時刻の差に基づいて、弾性波の発生源(以下「弾性波源」という)の位置を標定することができる。外部から構造物の表面に衝撃を与えた場合にも構造物内部で弾性波が発生する。このような場合にも各センサで検出した弾性波の到達時刻の差に基づいて弾性波源の位置を標定することができる。
構造物内部における弾性波の伝搬経路に損傷がある場合、弾性波の伝搬が妨げられる。構造物内部の損傷により弾性波の伝搬が妨げられると、一部のセンサで弾性波を検出することができない。その結果、弾性波源の標定結果の精度が低下してしまう。降雨時の雨滴による路面への衝突のように空間的に一様に付与される衝撃を構造物の表面に与えて対面に設置された各センサによって弾性波を検出した場合、内部に損傷を有する領域では弾性波源の密度が低下して観測される。このような特性を利用して、構造物の劣化状態(構造物内部の損傷有無)を評価することができる。
従来、構造物の劣化状態の評価精度を向上させるために種々の補正方法が提案されている。例えば、構造物の表面に与えられる衝撃が一様ではない場合であっても、衝撃に応じて定まる補正値を用いて、位置標定に基づく情報を補正する技術が提案されている。しかしながら、従来のいずれの手法においても温度に関しては特に言及されていない。そのため、季節や時間帯等の異なる温度環境下で計測した場合には、温度による変化を補正することができず、評価精度を向上させることができない場合があった。
国際公開第2017/217034号 国際公開第2022/014004号 特開2022-70711号公報
本発明が解決しようとする課題は、構造物の劣化状態の評価精度を向上させることができる構造物評価システム、構造物評価装置、構造物評価方法及びコンピュータプログラムを提供することである。
実施形態の構造物評価システムは、複数のセンサと、位置標定部と、補正部と、評価部とを持つ。複数のセンサは、構造物から発生した弾性波を検出する。位置標定部は、前記複数のセンサそれぞれによって検出された複数の弾性波に基づいて、前記複数の弾性波の発生源の位置を標定する。補正部は、前記構造物の温度に応じて定まる補正値を用いて、前記位置標定部による位置標定に基づく情報を補正する。評価部は、補正された前記情報に基づいて前記構造物の劣化状態を評価する。記構造物は、少なくとも車両が走行可能な舗装部と、前記舗装部を支持する床版部とで構成される。前記構造物の温度は、前記舗装部の温度に関連する値である。
第1の実施形態における構造物評価システムの構成を示す図。 第1の実施形態における信号処理部の機能構成例を示す図。 第1の実施形態における弾性波の計測方法を説明するための図。 第1の実施形態における音響インピーダンスと透過率との関係を示す図。 第1の実施形態における温度と弾性波の密度との間の関係性を説明するための図。 第1の実施形態における補正部が行う補正処理の説明図。 第1の実施形態における構造物評価システムによる劣化状態の評価処理の流れを示すシーケンス図。 第1の実施形態における効果を説明するための図。 第2の実施形態における補正部が行う補正処理の説明図。
以下、実施形態の構造物評価システム、構造物評価装置、構造物評価方法及びコンピュータプログラムを、図面を参照して説明する。
(概要)
発明者らの検討によって、温度によって弾性波源の密度の計測結果が変動することが分かった。そのため、従来の補正方法では、温度による変化を補正することができず、季節や時間帯等異なる温度環境下で計測した場合に、構造物の劣化状態の評価精度を向上させることできない場合があることが分かった。そこで、実施形態における構造物評価システムでは、弾性波の計測時における温度の情報に基づいて、弾性波源の位置標定に基づく情報を補正することで、計測時の温度による影響を踏まえて、構造物の劣化状態の評価精度を向上させる。ここで、位置標定に基づく情報は、弾性波源の位置標定を行う前に用いる情報(例えば、弾性波の計測閾値)、又は、標定結果を用いて得られる情報である。
(第1の実施形態)
図1は、第1の実施形態における構造物評価システム100の構成を示す図である。構造物評価システム100は、構造物50の健全性の評価に用いられる。以下の説明において、評価とは、ある基準に基づいて構造物50の健全性の度合い、すなわち構造物50の劣化状態を決定することを意味する。以下の実施形態における構造物50は、温度に対する応答の異なる少なくとも2つの材料で構成される。例えば、構造物50が、車両が走行する橋梁である場合、構造物50は走行可能な舗装部と、舗装部を支持する床版部とで構成される。舗装部51は例えばアスファルトであり、床版部52は例えばコンクリートである。
以下の説明では、構造物50が橋梁である場合を例に説明するが、構造物50は橋梁に限定される必要はない。構造物50は、温度に対する応答の異なる少なくとも2つの材料で構成される構造物であって、亀裂の発生または進展、あるいは外的衝撃(例えば雨、人工雨など)に伴い弾性波11が発生する構造物であればどのようなものであってもよい。なお、橋梁は、河川や渓谷等の上に架設される構造物に限らず、地面よりも上方に設けられる種々の構造物(例えば高速道路の高架橋)なども含む。
構造物50の劣化状態の評価に影響を及ぼす損傷としては、例えば亀裂、空洞、土砂化等の弾性波11の伝搬を妨害する構造物内部の損傷がある。ここで、亀裂には、縦方向の亀裂、横方向の亀裂及び斜め方向の亀裂等が含まれる。縦方向の亀裂とは、路面に対して垂直な方向に生じている亀裂である。横方向の亀裂とは、路面に対して水平な方向に生じている亀裂である。斜め方向の亀裂とは、路面に対して水平及び垂直以外の方向に生じている亀裂である。土砂化とは、主にアスファルトとコンクリート床版の境界部でコンクリートが土砂状に変化する劣化である。
構造物評価システム100は、複数のセンサ20-1~20-n、信号処理部30及び構造物評価装置40を備える。複数のセンサ20-1~20-nそれぞれと信号処理部30とは、有線により通信可能に接続される。信号処理部30と構造物評価装置40とは、有線又は無線により通信可能に接続される。なお、以下の説明では、センサ20-1~20-nを区別しない場合にはセンサ20と記載する。
図1に示すように、構造物50上を車両10が通過した際、車両10の走行部Wと路面との接触により、路面に対して荷重がかかる。荷重によるたわみにより、多数の弾性波11が構造物50内に発生する。構造物50の下面(例えば、床版部52)に設置された各センサ20は、構造物50内で発生した弾性波11を検出することができる。なお、ここでは、車両10により構造物50内で弾性波が発生する例を示しているが、構造物50内で弾性波を発生させる方法は他の方法であってもよい。例えば、外部から構造物50に衝撃を与えることにより弾性波を発生させてもよい。
センサ20は、圧電素子を有し、構造物50内部から発生する弾性波11を検出する。センサ20は、構造物50の面上における弾性波11を検出することが可能な位置に設置される。例えば、センサ20-1~20-nは、路面、側面及び底面のいずれかの面上に、車両走行軸方向と、車両走行軸直交方向にそれぞれ同一もしくは異なる間隔で離間して設置される。車両走行軸方向とは、車両が路面を走行する方向を表す。車両走行軸直交方向とは、車両走行軸方向に垂直な方向を表す。センサ20は、検出した弾性波11を電気信号に変換する。以下の説明では、センサ20が、構造物50の底面に設置されている場合を例に説明する。
センサ20には、例えば10kHz~1MHzの範囲に感度を有する圧電素子が用いられる。センサ20は、周波数範囲内に共振ピークをもつ共振型、共振を抑えた広帯域型等の種類があるが、センサ20の種類はいずれでもよい。センサ20が弾性波11を検出する方法は、電圧出力型、抵抗変化型及び静電容量型等があるが、いずれの検出方法でもよい。
センサ20に代えて加速度センサが用いられてもよい。この場合、加速度センサは、構造物50内部で発生する弾性波11を検出する。そして、加速度センサは、センサ20と同様の処理を行うことによって、検出した弾性波11を電気信号に変換する。
センサ20と信号処理部30との間には、例えば不図示の増幅器及びA/D変換器が設けられる。
増幅器は、センサ20から出力された電気信号を増幅する。増幅器は、増幅後の電気信号をA/D変換器に出力する。増幅器は、例えばA/D変換器において処理ができる程度に電気信号を増幅する。
A/D変換器は、増幅された電気信号を量子化してデジタル信号に変換する。A/D変換器は、デジタル信号を信号処理部30に出力する。
信号処理部30は、A/D変換器から出力されたデジタル信号を入力とする。信号処理部30は、入力したデジタル信号に対して信号処理を行う。信号処理部30が行う信号処理は、例えば、ノイズ除去、パラメータ抽出等である。信号処理部30は、信号処理後のデジタル信号を含む送信データを生成する。信号処理部30は、生成した送信データを構造物評価装置40に出力する。
信号処理部30は、アナログ回路又はデジタル回路を用いて構成される。デジタル回路は、例えばFPGA(Field Programmable Gate Array)やマイクロコンピュータにより実現される。デジタル回路は、専用のLSI(Large-Scale Integration)により実現されてもいい。また信号処理部30は、フラッシュメモリ等の不揮発メモリや、取り外し可能なメモリを搭載してもよい。
図2は、第1の実施形態における信号処理部30の機能構成例を示す図である。信号処理部30は、波形整形フィルタ301、ゲート生成回路302、到達時刻決定部303、特徴量抽出部304、送信データ生成部305、メモリ306及び出力部307を備える。
波形整形フィルタ301は、入力されたデジタル信号から所定の帯域外のノイズ成分を除去する。波形整形フィルタ301は、例えばデジタルバンドパスフィルタ(BPF)である。波形整形フィルタ301は、ノイズ成分除去後のデジタル信号(以下「ノイズ除去信号」という。)をゲート生成回路302及び特徴量抽出部304に出力する。
ゲート生成回路302は、波形整形フィルタ301から出力されたノイズ除去信号を入力とする。ゲート生成回路302は、入力したノイズ除去信号に基づいてゲート信号を生成する。ゲート信号は、ノイズ除去信号の波形が持続しているか否かを示す信号である。
ゲート生成回路302は、例えばエンベロープ検出器及びコンパレータにより実現される。エンベロープ検出器は、ノイズ除去信号のエンベロープを検出する。エンベロープは、例えばノイズ除去信号を二乗し、二乗した出力値に対して所定の処理(例えばローパスフィルタを用いた処理やヒルベルト変換)を行うことで抽出される。コンパレータは、ノイズ除去信号のエンベロープが所定の閾値以上であるか否かを判定する。
ゲート生成回路302は、ノイズ除去信号のエンベロープが所定の閾値以上となった場合、ノイズ除去信号の波形が持続していることを示す第1のゲート信号を到達時刻決定部303及び特徴量抽出部304に出力する。一方、ゲート生成回路302は、ノイズ除去信号のエンベロープが所定の閾値未満になった場合、ノイズ除去信号の波形が持続していないことを示す第2のゲート信号を到達時刻決定部303及び特徴量抽出部304に出力する。なお、ゲート生成回路302において、ノイズ除去信号の波形が持続しているか否かをエンベロープに基づいて判定する構成を示したが、ゲート生成回路302は、ノイズ除去信号そのものや、絶対値を適用した信号に対して処理を行ってもよい。このゲート生成に用いる閾値を計測閾値と記載する。
到達時刻決定部303は、不図示の水晶発振器などのクロック源から出力されるクロックと、ゲート生成回路302から出力されたゲート信号とを入力とする。到達時刻決定部303は、第1のゲート信号が入力されている間に入力されたクロックを用いて、弾性波到達時刻を決定する。到達時刻決定部303は、決定した弾性波到達時刻を時刻情報として送信データ生成部305に出力する。到達時刻決定部303は、第2のゲート信号が入力されている間には処理を行わない。到達時刻決定部303は、クロック源からの信号をもとに、電源投入時からの累積の時刻情報を生成する。具体的には、到達時刻決定部303は、クロックのエッジをカウントするカウンタとし、カウンタのレジスタの値を時刻情報とすればよい。カウンタのレジスタは所定のビット長を有するように決定される。
特徴量抽出部304は、波形整形フィルタ301から出力されたノイズ除去信号と、ゲート生成回路302から出力されたゲート信号とを入力とする。特徴量抽出部304は、第1のゲート信号が入力されている間に入力されたノイズ除去信号を用いて、ノイズ除去信号の特徴量を抽出する。特徴量抽出部304は、第2のゲート信号が入力されている間に処理を行わない。特徴量は、ノイズ除去信号の特徴を示す情報である。
特徴量は、例えば波形の振幅[mV]、波形の立ち上がり時間[usec]、ゲート信号の持続時間[usec]、ゼロクロスカウント数[times]、波形のエネルギー[arb.]、周波数[Hz]及びRMS(Root Mean Square:二乗平均平方根)値等である。特徴量抽出部304は、抽出した特徴量に関するパラメータを送信データ生成部305に出力する。特徴量抽出部304は、特徴量に関するパラメータを出力する際に、特徴量に関するパラメータにセンサIDを対応付ける。センサIDは、構造物50の健全性の評価対象となる領域(以下「評価領域」という。)に設置されているセンサ20を識別するための識別情報を表す。
波形の振幅は、例えばノイズ除去信号の中で最大振幅の値である。波形の立ち上がり時間は、例えばゲート信号の立ち上がり開始からノイズ除去信号が最大値に達するまでの時間T1である。ゲート信号の持続時間は、例えばゲート信号の立ち上がり開始から振幅が予め設定される値よりも小さくなるまでの時間である。ゼロクロスカウント数は、例えばゼロ値を通る基準線をノイズ除去信号が横切る回数である。
波形のエネルギーは、例えば各時点においてノイズ除去信号の振幅を二乗したものを時間積分した値である。なお、エネルギーの定義は、上記例に限定されず、例えば波形の包絡線を用いて近似されたものでもよい。周波数は、ノイズ除去信号の周波数である。RMS値は、例えば各時点においてノイズ除去信号の振幅を二乗して平方根により求めた値である。
送信データ生成部305は、センサIDと、時刻情報と、特徴量に関するパラメータとを入力とする。送信データ生成部305は、入力したセンサIDと、時刻情報と、特徴量に関するパラメータとを含む送信データを生成する。
メモリ306は、送信データを記憶する。メモリ306は、例えばデュアルポートRAM(Random Access Memory)である。
出力部307は、メモリ306に記憶されている送信データを構造物評価装置40に逐次出力する。
図1に戻って説明を続ける。構造物評価装置40は、通信部41、制御部42、記憶部43及び表示部44を備える。
通信部41は、信号処理部30から出力された送信データを受信する。
制御部42は、構造物評価装置40全体を制御する。制御部42は、CPU(Central Processing Unit)等のプロセッサやメモリを用いて構成される。制御部42は、プログラムを実行することによって、取得部421、イベント抽出部422、位置標定部423、分布生成部424、補正部425及び評価部426として機能する。
取得部421、イベント抽出部422、位置標定部423、分布生成部424、補正部425及び評価部426の機能部のうち一部または全部は、ASIC(Application Specific Integrated Circuit)やPLD(Programmable Logic Device)、FPGAなどのハードウェアによって実現されてもよいし、ソフトウェアとハードウェアとの協働によって実現されてもよい。プログラムは、コンピュータ読み取り可能な記録媒体に記録されてもよい。コンピュータ読み取り可能な記録媒体とは、例えばフレキシブルディスク、光磁気ディスク、ROM、CD-ROM等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置などの非一時的な記憶媒体である。プログラムは、電気通信回線を介して送信されてもよい。
取得部421、イベント抽出部422、位置標定部423、分布生成部424、補正部425及び評価部426の機能の一部は、予め構造物評価装置40に搭載されている必要はなく、追加のアプリケーションプログラムが構造物評価装置40にインストールされることで実現されてもよい。
取得部421は、各種情報を取得する。例えば、取得部421は、通信部41によって受信された送信データを取得する。例えば、取得部421は、構造物50の温度情報を取得する。構造物50の温度情報は、評価対象となる構造物50付近の温度情報である。評価対象となる構造物50付近の温度情報は、例えば評価対象となる構造物50付近の気温であってもよいし、日射などによる温度上昇の影響を直接受ける構造物50の影響部分(例えば、舗装部)の温度であってもよい。
取得部421は、構造物50の温度情報を取得することができればどのような方法で取得してもよい。例えば、取得部421は、構造物50付近に設置された温度センサによって計測された温度情報を取得してもよいし、天気予報等で提供される構造物50付近の温度情報を取得してもよいし、影響部分から温度情報を直接取得してもよい。取得部421が、構造物50付近に設置された温度センサによって計測された温度情報を取得する場合、又は、影響部分から温度情報を直接取得する場合には、構造物評価システム100は温度センサを備える。取得部421は、取得した送信データや温度情報を記憶部43に保存する。第1の実施形態において取得部421は、評価対象となる構造物50を走行した車両10に関する情報(以下「車両情報」という。)を所定期間分取得する。
車両情報は、例えば所定の期間に構造物50を走行した車両10の車種情報と、構造物50の車両10の通行量情報である。車両10の車種情報には、少なくとも車両10のトレッド幅の情報が含まれる。車両10の通行量情報は、所定の期間に構造物50をどの車種の車両10が何台通過したかを示す情報である。取得部421は、車両情報を、ユーザの入力によって取得してもよいし、交通情報を記憶しているサーバから取得してもよい。
イベント抽出部422は、記憶部43に記憶されている送信データの中から1イベントにおける送信データを抽出する。イベントとは、構造物50で起こった弾性波発生事象を表す。本実施形態における弾性波発生事象は、車両10による路面の通過である。1回のイベントが発生した場合、複数のセンサ20で略同時刻に弾性波11が検出されることになる。すなわち、記憶部43には、略同時刻に検出された弾性波11に関する送信データが記憶されていることになる。そこで、イベント抽出部422は、所定の時間窓を設け、到達時刻が時間窓の範囲内に存在する全ての送信データを1イベントにおける送信データとして抽出する。イベント抽出部422は、抽出した1イベントにおける送信データを位置標定部423に出力する。
時間窓の範囲Twは、対象とする構造物50における弾性波伝搬速度vと、最大のセンサ間隔dmaxを用いて、Tw≧dmax/vの範囲になるように決定してもよい。誤検出を避けるためには、Twをできるだけ小さい値に設定することが望ましいため、実質的にはTw=dmax/vとすることができる。弾性波伝搬速度vは、予め求められていてもよい。
位置標定部423は、センサ位置情報と、イベント抽出部422によって抽出された複数の送信データそれぞれに含まれるセンサID及び時刻情報とに基づいて弾性波源の位置を標定する。
センサ位置情報には、センサIDに対応付けてセンサ20の設置位置に関する情報が含まれる。センサ位置情報は、例えば緯度および経度、あるいは構造物50の基準となる位置からの水平方向および垂直方向の距離などのセンサ20の設置位置に関する情報を含む。位置標定部423は、センサ位置情報を予め保持している。センサ位置情報は、位置標定部423が弾性波源の位置標定を行う前であればどのタイミングで位置標定部423に記憶されてもよい。
センサ位置情報は、記憶部43に記憶されていてもよい。この場合、位置標定部423は、位置標定を行うタイミングで記憶部43からセンサ位置情報を取得する。弾性波源の位置の標定には、カルマンフィルタ、最小二乗法などが用いられてもよい。位置標定部423は、計測期間中に得られた弾性波源の位置情報を分布生成部424に出力する。
分布生成部424は、位置標定部423から出力された複数の弾性波源の位置情報を入力とする。分布生成部424は、入力した複数の弾性波源の位置情報を用いて、弾性波源分布を生成する。弾性波源分布は、弾性波源の位置が示された分布を表す。より具体的には、弾性波源分布は、横軸を通行方向の距離とし、縦軸を幅方向の距離として、評価対象となる構造物50を表した仮想的なデータ上において弾性波源の位置を示す点が示された分布である。分布生成部424は、弾性波源分布を用いて弾性波源密度分布を生成する。弾性波源密度分布は、弾性波源分布において予め定められたエリア毎に、各エリアに含まれる弾性波源の数に応じて求められる密度の値が示された分布を表す。
補正部425は、取得部421によって取得された構造物50の温度情報に応じて定まる補正値を用いて、位置標定部423による位置標定に基づく情報を補正する。構造物50の温度情報は、舗装部51の温度に関連する値であり、例えば気温、床版部52下面の温度、舗装部51を走行する車両又は人手で計測した路面温度(舗装部51の温度)、日射量、計測中の温度の統計値、舗装部51の温度と床版部52の温度の差、舗装部51の温度と床版部52の温度の比のいずれかである。舗装部51の温度は、気温又は日射量の少なくともいずれかにより推定された値であってもよい。床版部52の温度は、床版下面の温度である。
第1の実施形態において構造物50の温度情報に応じて定まる補正値とは、弾性波源密度分布の密度を補正するための補正値である。第1の実施形態において位置標定に基づく情報は、弾性波源密度分布である。
評価部426は、補正された情報に基づいて構造物50の劣化状態を評価する。具体的には、評価部426は、補正された弾性波源密度分布を用いて構造物50の劣化状態を評価する。例えば、評価部426は、弾性波源の密度が閾値以上の領域を健全な領域と評価し、弾性波源の密度が閾値未満の領域を損傷領域と評価する。弾性波源と、各センサ20それぞれまでの弾性波の伝搬経路を示す波線が通過した領域が評価領域となる。そこで、評価部426は、評価領域内において構造物の劣化状態を評価する。
記憶部43には、取得部421によって取得された送信データ、温度情報、車両情報及び補正テーブルが記憶される。補正テーブルは、補正部425が補正を行う際に使用するテーブルである。補正テーブルには、各温度に対するヤング率の値が対応付けて登録される。例えば、補正テーブルには、温度20℃の場合のヤング率E20の値や、温度25℃の場合のヤング率E25の値等が登録される。これにより、取得部421によって温度情報をヤング率に変換して処理に用いることができる。記憶部43は、磁気ハードディスク装置や半導体記憶装置などの記憶装置を用いて構成される。
表示部44は、評価部426の制御に従って評価結果を表示する。さらに、表示部44は、評価部426の制御に従って、弾性波源と、各センサ20までの弾性波の伝搬経路を示す波線群を例えば投影法により表示する。表示部44は、液晶ディスプレイ、有機EL(Electro Luminescence)ディスプレイ等の画像表示装置である。表示部44は、画像表示装置を構造物評価装置40に接続するためのインタフェースであってもよい。この場合、表示部44は、評価結果を表示するための映像信号を生成し、自身に接続されている画像表示装置に映像信号を出力する。
次に、補正部425による補正処理について具体的に説明する。まず補正部425による補正処理を説明するにあたり必要になる温度と弾性波源の密度との関係性について説明する。図3は、第1の実施形態における弾性波の計測方法を説明するための図である。図3では、弾性波の計測方法の一例として、路面を走行する車両の走行部W(例えば、タイヤ)が発生させる弾性波に基づいて、構造物50の内部損傷を評価することを考える。ここでは、構造物50が、車両が走行する道路である場合を例に説明する。なお、車両については特に限定はなく、一般の走行車でもよい。一般の走行車が発生する弾性波を利用することで、検査のために車両規制の必要がなく、利用者の利便性向上に寄与することになる。車両の走行部Wと路面とが接触することで、路面に対して荷重がかかる。これにより、構造物50内に弾性波が発生する。発生した弾性波は、構造物50内を伝搬して、路面とは異なる面(例えば、底面)に設置された各センサ20で検出される。
構造物50は、通常、舗装部51と床版部52で構成される。一般的な舗装部51の厚みは80mm程度、床版部52の厚みは180mm~230mm程度である。路面を速度Vで走行する車両の走行部Wと、舗装部51の面の相互作用により生じた弾性波(振幅Asrc)は床版部52下面へ向かって伝搬する。途中、弾性波は、舗装部51と床版部52の界面53でそれぞれの音響インピーダンスの違いにより一部が反射、一部が透過する。この時の透過率をLとする。さらに床版部52下面へ伝搬した弾性波はセンサ20に到達し、センサ20により検出される。検出される弾性波の振幅をAとする。床版部52下面には、複数のセンサ20が設置され、弾性波はそれぞれ異なる時間差を持って各センサ20に検出される。各センサによって検出された際の時間差に基づいて、弾性波源SRの位置が標定される。
次に、温度変化が生じた場合の弾性波の挙動について説明する。材料のヤング率は、温度依存性を有することが知られている。これは、温度Tが上がると、ポテンシャル特性に基づいて平均原子間距離rが増加することに起因している(熱膨張)。通常、舗装部51のアスファルトと、床版部52のコンクリートでは、直接日射の影響を受ける舗装部51の温度変化に伴う熱膨張の影響が支配的となる。熱膨張に伴い、ヤング率Eも変化することになる。例えば、原子間距離rが増加すると、ヤング率Eが下がるという関係がある。温度とヤング率の関係は、以下の式(1)のWachtman Equationで近似できることが知られている。温度とヤング率の関係は、参考文献1に具体的に示されている。
(参考文献1:横浜国立大学, “材料と温度と強度”, [online] ,[令和4年9月8日検索], インターネット<URL: http://www.materep.ynu.ac.jp/archives/info/lec9>)
式(1)における係数Bは材料固有の値である。式(1)により、温度が下がるとヤング率Eも下がるという関係が導かれる。通常の材料ではT=0Kでは傾きがゼロに近づくが、300K前後ではほぼ線形(準線形関係)とみなせるとされている。ヤング率Eと体積弾性率Kの間には、ポアソン比νを介して以下の式(2)に示す関係がある。
式(2)により、ヤング率Eが下がると体積弾性率Kも下がるという関係が導かれる。さらに音響インピーダンスZは、密度ρと体積弾性率Kで以下の式(3)のように定義される。
式(3)により、体積弾性率Kが下がると音響インピーダンスZも下がるという関係が導かれる。ここで、図4のように、音響インピーダンスの異なる2つの材料(材料M1,M2)を伝搬する弾性波を考える。図4において、Z<Zrefである。音響インピーダンスZの材料M1から、音響インピーダンスZrefの材料M2へ透過する弾性波の透過率Lは、以下の式(4)のように表すことができる。特に、Z<<Zrefの場合の透過率Lは、2Z/Zrefの近似式で表すことができる。
アスファルトとコンクリートの音響インピーダンスは通常、式(4)に示す近似式が適用できる条件となる。例えば、材料M1が舗装部51に相当し、材料M2が床版部52に相当する。舗装部51は、日射が直接当たるため、床版部52に比べて温度変化が大きい。したがって、透過率Lは、舗装部51に相当する音響インピーダンスZの変化が大きく作用する。そのため、式(4)により、音響インピーダンスZが下がると透過率Lも下がるという関係が導かれる。材料減衰を無視すると、センサ20に到達する弾性波の振幅Aは、以下の式(5)に示すように弾性波源の振幅Asrcに透過率Lを乗じた値となる。
式(5)により、透過率Lが下がると弾性波の振幅Aも下がるという関係が導かれる。センサ20で検出された弾性波の対数振幅AdBが所定の計測閾値を超過した場合に、ヒットとして認識される。このように、対数振幅AdBが所定の計測閾値を超えない場合には、ヒットとして認識されない。そのため、弾性波の振幅Aと弾性波のヒット数(AEヒット数)とは、弾性波の振幅Aが下がれば弾性波のヒット数が下がるという関係になる。さらに、対数振幅AdBとヒット数Naeの間には、経験的に以下の式(6)に示す関係があることが知られている。これは地震学におけるGutenberg-Richter則と類似する。Gutenberg-Richter則については、参考文献2に記載されている。
(参考文献2:B. Gutenberg and C. F. Richter, “Seismicity of the Earth and Associated Phenomena”, Princeton Univ. Press, 1949.)
ヒット数Naeのうち一部が複数のセンサ20で検出され、イベント数Nevとして認識される。この時のイベント率をPevで表す。イベントとは、構造物50で起こった弾性波発生事象を表す。本実施形態における弾性波発生事象は、例えば車両による路面の通過である。1回のイベントが発生した場合、複数のセンサ20で略同時刻に弾性波が検出される(ヒットが発生する)ことになる。そのため、イベント数Nevは、同時に複数のセンサ20で検出されたヒット数Naeに起因する。したがって、ヒット数Naeが下がればイベント数Nevも下がるという関係になる。車両の通過数をncar、計測エリアをSとすると、単位面積、車両1台あたりの弾性波源の密度Dは以下の式(7)のように表すことができる。
式(7)により、イベント数Nevが下がると弾性波源の密度Dも下がるという関係が導かれる。これらの関係を整理すると、舗装部51のヤング率Eと温度Tの間には、負の相関関係があり、舗装部51のヤング率Eと弾性波源の密度Dの間には、両対数軸上での線形関係がある。すなわち、以下の式(8)の関係が成り立つ。
従って、温度Tと弾性波源の密度Dには、負の相関関係があることが見いだされる。上述した依存関係をまとめると、図5の関係が導かれる。図5は、温度Tと弾性波の密度Dとの間の関係性を説明するための図である。図5に示すように、舗装部51の温度T→原子間距離→ヤング率→体積弾性率→音響インピーダンス→透過率→AE振幅→AEヒット数→AE標定数→AE密度のように影響が伝搬することが分かる。以下の説明では、原子間距離、ヤング率、体積弾性率、音響インピーダンス、透過率、AE振幅、AEヒット数及びAE標定数を、温度からAE密度までの関係性を示す中間パラメータとして説明する。中間パラメータは、舗装部51の温度に関連する値である。その結果、温度Tが変化すると弾性波の密度Dが変化するということが分かる。例えば、温度Tが上がれば弾性波の密度Dが下がり、温度Tが下がれば弾性波の密度Dが上がるということが分かる。
(第1の実施形態における補正方法について)
第1の実施形態における補正部425が行う補正処理について図6を用いて説明する。図6において、Eexは計測時の温度に対応するヤング率を表し、Dexは計測した弾性波源の密度(弾性波源密度分布内のあるエリアの密度)を表し、Erefは基準温度(例えば、20℃)に対応するヤング率を表し、Dcompは補正後の弾性波源の密度を表す。図6において、Eex,Dex,Eref,Dcompは対数で表される。さらに、図6において、補正直線60は、弾性波源の密度Dとヤング率Eとの関係に基づく補正に用いられる直線である。補正直線60は、予め複数の温度条件で求めた弾性波源の密度Dとヤング率Eの両対数プロットを作成することで得られる。なお、補正直線60は、予め複数の温度条件で求めた弾性波源の密度Dとヤング率Eに応じて所定の傾きを有する。
例えば、補正直線60は、以下のように求められる。事前に、異なる温度条件(例えば、20℃、25℃、30…等)で、弾性波源の密度Dとヤング率Eを求める。各温度条件で求められた弾性波源の密度Dとヤング率Eそれぞれについて対数をとる。これにより得られる各温度条件に基づく弾性波源の密度Dの対数とヤング率Eの対数との両対数プロットを作成することで補正直線60が得られる。このように得られた補正直線60は記憶部43に記憶されていてもよい。
補正部425は、弾性波の計測時の温度情報と、記憶部43に記憶されている補正テーブルと、補正直線60とに基づいて弾性波源密度分布を補正する。この際、補正部425は、弾性波源密度分布におけるエリア毎に密度を補正することで弾性波源密度分布を補正する。まず補正部425は、取得部421から温度情報を取得する。次に、補正部425は、補正テーブルを参照して、取得した温度情報で示される温度に対応するヤング率の値を取得する。例えば、温度が温度Texである場合、補正部425は、計測時の温度Texに対応するヤング率Eexを補正テーブルから取得する。
次に、補正部425は、弾性波源密度分布における1つエリア内の弾性波源の密度Dexを取得する。補正部425は、弾性波源の密度Dとヤング率Eの両対数プロット(横軸がヤング率Eの対数を表し、縦軸が弾性波源の密度Dの対数を表す)上において、ヤング率Eexと弾性波源の密度Dexとの交点に補正直線60を設定する。補正部425は、弾性波源の密度Dとヤング率Eの両対数プロット上に設定した補正直線60に従って、基準温度に対応するヤング率Erefから導かれる補正後の弾性波源の密度Dcompを補正値として取得する。そして、補正部425は、1つエリアの弾性波源の密度Dを補正値(例えば、弾性波源の密度Dcomp)になるように補正する。補正部425は、同様の処理をエリア毎に行う。これにより、弾性波源密度分布が補正される。
図7は、第1の実施形態における構造物評価システム100による劣化状態の評価処理の流れを示すシーケンス図である。図7の処理は、評価対象となる構造物50を車両10が走行したことに応じて実行される。
評価対象となる構造物50を車両10が走行すると、車両10の走行部が路面に接触する。これにより、構造物50内で弾性波11が発生する。複数のセンサ20それぞれは、構造物50内で発生した弾性波11を検出する(ステップS101)。複数のセンサ20それぞれは、検出した弾性波11を電気信号に変換して信号処理部30に出力する(ステップS102)。複数のセンサ20それぞれから出力された電気信号は、不図示の増幅器に増幅される。増幅後の電気信号は、A/D変換器によってデジタル信号に変換される。
信号処理部30は、A/D変換器から出力されたデジタル信号を入力する。信号処理部30の到達時刻決定部303は、各弾性波11の到達時刻を決定する(ステップS103)。具体的には、到達時刻決定部303は、第1のゲート信号が入力されている間に入力されたクロックを用いて弾性波到達時刻を決定する。到達時刻決定部303は、決定した弾性波到達時刻を時刻情報として送信データ生成部305に出力する。到達時刻決定部303は、この処理を、入力した全てのデジタル信号に対して行う。
信号処理部30の特徴量抽出部304は、第1のゲート信号が入力されている間に入力されたデジタル信号であるノイズ除去信号を用いて、ノイズ除去信号の特徴量を抽出する(ステップS104)。特徴量抽出部304は、抽出した特徴量に関するパラメータを送信データ生成部305に出力する。送信データ生成部305は、センサIDと、時刻情報と、特徴量に関するパラメータとを含む送信データを生成する(ステップS105)。出力部307は、送信データを構造物評価装置40に逐次出力する(ステップS106)。
構造物評価装置40の通信部41は、信号処理部30から出力された送信データを受信する。取得部421は、通信部41によって受信された送信データを取得する。取得部421は、取得した送信データを記憶部43に記録する(ステップS107)。イベント抽出部422は、記憶部43に記憶されている送信データの中から1イベントにおける送信データを抽出する。イベント抽出部422は、抽出した1イベントにおける送信データを位置標定部423及び分布生成部424に出力する。
位置標定部423は、イベント抽出部422から出力された送信データに含まれるセンサID及び時刻情報と、予め保持しているセンサ位置情報とに基づいて弾性波源の位置を標定する(ステップS108)。具体的には、まず位置標定部423は、複数のセンサ20それぞれへの弾性波11の到達時刻の差を算出する。次に、位置標定部423は、センサ位置情報と、到達時刻の差の情報とを用いて弾性波源の位置を標定する。
位置標定部423は、計測期間中にイベント抽出部422から1イベントの送信データが出力される度にステップS108の処理を実行する。これにより、位置標定部423は、複数の弾性波源の位置を標定する。位置標定部423は、複数の弾性波源の位置情報を分布生成部424に出力する。取得部421は、外部から温度情報を取得する(ステップS109)。取得部421は、取得した温度情報を補正部425に出力する。なお、取得部421が温度情報を取得するタイミングは、このタイミングに限らず、図7の処理が開始される前、又は、補正部425による補正が行われる前であればよい。
分布生成部424は、位置標定部423から出力された複数の弾性波源の位置情報を用いて弾性波源分布を生成する。具体的には、分布生成部424は、得られた複数の弾性波源の位置情報で示される弾性波源の位置を、仮想的なデータ上にプロットすることによって弾性波源分布を生成する。分布生成部424は、生成した弾性波源分布を用いて弾性波源密度分布を生成する(ステップS110)。具体的には、まず分布生成部424は、弾性波源分布を所定の区画に区切ることによって複数のエリアに分割する。次に、分布生成部424は、エリア毎に弾性波源の密度を算出する。そして、分布生成部424は、エリア毎に、算出したエリア毎の弾性波源の密度の値を割り当てることによって弾性波源密度分布を生成する。分布生成部424は、生成した弾性波源密度分布を補正部425に出力する。
補正部425は、分布生成部424から出力された弾性波源密度分布と、記憶部43に記憶されている補正テーブルと、取得部421から出力された温度情報とを入力する。補正部425は、入力した弾性波源密度分布と、補正テーブルと、温度情報とに基づいて弾性波源密度分布を補正する(ステップS111)。補正部425は、補正後の弾性波源密度分布を評価部426に出力する。
評価部426は、補正後の弾性波源密度分布を用いて、構造物の劣化状態を評価する。評価部426は、評価結果を表示部44に出力する。表示部44は、評価部426から出力された評価結果を表示する(ステップS112)。例えば、表示部44は、評価結果として、補正後の弾性波源密度分布を表示してもよいし、損傷領域とみられる領域を他の領域と表示態様を変えて表示してもよい。
図8は、第1の実施形態における効果を説明するための図である。図8の左図は、異なる温度環境下で測定した際の弾性波源の密度の評価結果を表し、図8の右図は、実施形態における手法に基づく温度補正結果を表す。温度補正前(図8の左図)の平均の弾性波源の密度は2.0(+0.4~-0.43)であるのに対して、温度補正を行った後(図8の右図)は平均2.2(+0.28~-0.28)となり、誤差範囲が32%減少したことが分かる。
以上のように構成された構造物評価システム100では、温度に対応するヤング率と弾性波源の密度とに基づいて弾性波源密度分布を補正する。これにより、温度による影響を加味した弾性波源密度分布を生成することができる。このように得られた弾性波源分布を用いて構造物の劣化状態を評価することで、構造物の劣化状態の評価精度を向上させることが可能になる。
(第2の実施形態)
第2の実施形態では、第1の実施形態よりも簡易な方法で補正を行う構成について説明する。第2の実施形態において、システム構成及び各装置の構成は第1の実施形態と同様である。第2の実施形態では、補正テーブルを用いない点、補正直線及び補正部425の処理が第1の実施形態と異なる。以下、相違点を中心に説明する。
補正部425は、取得部421によって取得された構造物50の温度情報に応じて定まる補正値を用いて、位置標定部423による位置標定に基づく情報を補正する。第2の実施形態において構造物50の温度情報に応じて定まる補正値とは、弾性波源密度分布の密度を補正するための補正値である。第2の実施形態において位置標定に基づく情報は、弾性波源密度分布である。
(第2の実施形態における補正方法について)
第2の実施形態における補正部425が行う補正処理について図9を用いて説明する。図9において、Texは計測時の温度を表し、Dexは計測した弾性波源の密度(あるエリアの密度)を表し、Trefは基準温度(例えば、20℃)を表し、Dcompは補正後の弾性波源の密度を表す。さらに、図9において、補正直線61は、弾性波源の密度Dと温度Tとの関係に基づく補正に用いられる直線である。補正直線61は、予め複数の温度条件で求めた弾性波源の密度Dと温度Tをプロットした時の傾きを線形近似により求めることで得られる直線である。なお、補正直線61は、予め複数の温度条件で求めた弾性波源の密度Dと温度Tに応じて所定の傾きを有する。補正直線61は記憶部43に記憶されていてもよい。
補正部425は、弾性波の計測時の温度情報と、補正直線61とに基づいて弾性波源密度分布を補正する。この際、補正部425は、弾性波源密度分布におけるエリア毎に密度を補正することで弾性波源密度分布を補正する。まず補正部425は、取得部421から温度情報を取得する。例えば、補正部425が取得した温度情報で示される温度がTexであるとする。次に、補正部425は、弾性波源密度分布における1つエリア内の弾性波源の密度Dexを取得する。補正部425は、弾性波源の密度Dと温度Tの関係を表すグラフ(横軸が温度Tを表し、縦軸が弾性波源の密度Dを表す)上において、温度Texと弾性波源の密度Dexとの交点に補正直線61を設定する。補正部425は、弾性波源の密度Dと温度Tの関係を表すグラフ上に設定した補正直線61に従って、基準温度Trefから導かれる補正後の弾性波源の密度Dcompを補正値として取得する。そして、補正部425は、1つエリアの弾性波源の密度Dを補正値(例えば、弾性波源の密度Dcomp)になるように補正する。補正部425は、同様の処理をエリア毎に行う。これにより、弾性波源密度分布が補正される。
このように、第2の実施形態における補正部425は、取得した温度情報で示される温度が高い場合には弾性波源の密度Dを所定の割合で高くなる方向に補正し、温度が低い場合には弾性波源の密度Dを所定の割合で低くなる方向に補正する。一般的に舗装部と気温(温度)とは比例関係にあることが知られている(例えば、参考文献3参照)。
(参考文献3:Li. T,“Influencing Parameters on Tire-Pavement Interaction Noise: Review, Experiments”, and Design Considerations. Designs 2018, 2, 38.)
参考文献1によれば、予め複数の温度条件で求めた、弾性波源の密度Dと気温との関係をプロットしたときの傾きを線形近似により求め、補正してもよいことが分かる。特に、気温が高い時には密度を所定の割合で高くなる方向に補正する、また気温が低い時には、密度を所定の割合で低くなる方向に補正するように補正値を求めることができる。
(第3の実施形態)
第3の実施形態では、中間パラメータに基づいて密度を補正する構成について説明する。第3の実施形態において、システム構成及び各装置の構成は第1の実施形態及び第2の実施形態と同様である。第3の実施形態では、補正テーブルを用いない点、補正直線及び補正部425の処理が第1の実施形態及び第2の実施形態と異なる。以下、相違点を中心に説明する。
第3の実施形態における補正部425は、中間パラメータ(例えば、原子間距離、ヤング率、体積弾性率、音響インピーダンス、透過率、AE振幅、AEヒット数、AE標定数)のうちいずれか1つと、弾性波源の密度Dとの関係に基づいて得られる補正値を用いて、弾性波源の密度Dを補正する。
(第3の実施形態における補正方法について)
第3の実施形態における補正部425が行う補正処理について説明する。ここでは、中間パラメータとして音響インピーダンスZを用いる場合について説明する。なお、他の中間パラメータを用いた場合においては、以下に示す音響インピーダンスを他の中間パラメータにおける用語に置き換えればよい。ここで、Zexは計測時の音響インピーダンスを表し、Dexは計測した弾性波源の密度(あるエリアの密度)を表し、Zrefは所定の温度における基準音響インピーダンス(例えば、20℃における音響インピーダンス)を表し、Dcompは補正後の弾性波源の密度を表すものとする。中間パラメータとして音響インピーダンスZを用いる場合の補正直線は、弾性波源の密度Dと音響インピーダンスZとの関係に基づく補正に用いられる直線である。補正直線は、予め複数の音響インピーダンス条件で求めた弾性波源の密度Dと音響インピーダンスZをプロットした時の傾きを線形近似により求めることで得られる直線である。なお、補正直線は、予め複数の音響インピーダンス条件で求めた弾性波源の密度Dと音響インピーダンスZに応じて所定の傾きを有する。補正直線は記憶部43に記憶されていてもよい。
補正部425は、弾性波の計測時の音響インピーダンス情報と、補正直線とに基づいて弾性波源密度分布を補正する。この際、補正部425は、弾性波源密度分布におけるエリア毎に密度を補正することで弾性波源密度分布を補正する。まず補正部425は、取得部421から音響インピーダンス情報を取得する。例えば、補正部425が取得した音響インピーダンス情報で示される音響インピーダンスがZexであるとする。次に、補正部425は、弾性波源密度分布における1つエリア内の弾性波源の密度Dexを取得する。補正部425は、弾性波源の密度Dと音響インピーダンスZの関係を表すグラフ(横軸が音響インピーダンスZを表し、縦軸が弾性波源の密度Dを表す)上において、音響インピーダンスZexと弾性波源の密度Dexとの交点に補正直線を設定する。補正部425は、弾性波源の密度Dと音響インピーダンスZの関係を表すグラフ上に設定した補正直線に従って、基準音響インピーダンスZrefから導かれる補正後の弾性波源の密度Dcompを補正値として取得する。そして、補正部425は、1つエリアの弾性波源の密度Dを補正値(例えば、弾性波源の密度Dcomp)になるように補正する。補正部425は、同様の処理をエリア毎に行う。これにより、弾性波源密度分布が補正される。
以上のように構成された第3の実施形態における構造物評価システム100では、中間パラメータを用いた場合においても第1の実施形態や第2の実施形態と同様に、温度による影響を加味した弾性波源密度分布を生成することができる。このように得られた弾性波源分布を用いて構造物の劣化状態を評価することで、構造物の劣化状態の評価精度を向上させることが可能になる。
(第4の実施形態)
第1の実施形態から第3の実施形態では、密度を補正する構成を示した。第4の実施形態では、ヒットを検出するための閾値を変更することによってヒット数を補正する構成について説明する。第4の実施形態において、システム構成は第1の実施形態から第3の実施形態と同様である。第4の実施形態では、補正部425が信号処理部30に備えられる点と、補正部425の処理が第1の実施形態から第3の実施形態と異なる。以下、相違点を中心に説明する。
第4の実施形態では、補正部425が、信号処理部30に備えられる。補正部425は、構造物50の温度情報に応じて定まる補正値を用いて、位置標定部423による位置標定に基づく情報を補正する。第4の実施形態において構造物50の温度情報に応じて定まる補正値とは、ゲート生成回路の閾値を補正するための補正値である。第4の実施形態において位置標定に基づく情報は、ゲート生成回路の閾値である。補正部425は、構造物評価装置40に備えられてもよい。その場合は、構造物評価装置40内の補正部425から、信号処理部30内のゲート生成回路302に対して、ゲート生成回路302の閾値に関する情報が送信されるように構成することができる。
(第4の実施形態における補正方法について)
補正部425は、ゲート生成回路302における閾値を変更することでヒット数を補正する。より具体的には、補正部425は、予め定められた温度毎と閾値との対応付けがなされた補正テーブルを参照し、ゲート生成回路302における閾値を補正する。閾値を下げると、ヒットが多く検出され、閾値を上げるとヒットの検出が少なくなる。ヒット数と密度は相関関係があるため、補正部425は、弾性波源の密度Dを下げる場合には閾値を上げ、弾性波源の密度Dを上げる場合には閾値を下げるように補正を行う。例えば、補正部425は、計測前において温度情報を外部から取得する。温度情報は、構造物評価装置40の取得部421から取得してもよいし、その他から取得してもよい。補正部425は、補正テーブルを参照し、取得した温度情報で示される温度に対応する閾値を補正値として取得する。補正部425は、取得した補正値となるようにゲート生成回路302の閾値を補正する。
構造物評価装置40では、分布生成部424は、生成した弾性波源密度分布を評価部426に出力する。評価部426は、分布生成部424から出力された弾性波源密度分布に基づいて構造物の劣化状態を評価する。
(各実施形態に共通の変形例1)
信号処理部30は、構造物評価装置40に備えられてもよい。
(各実施形態に共通の変形例2)
上記の各実施形態では、複数のセンサ20-1~20-nが、1台の信号処理部30に接続されている構成を示した。構造物評価システム100は、複数台の信号処理部30を備え、各センサ20それぞれが異なる信号処理部30に接続されてもよい。
(各実施形態に共通の変形例3)
構造物評価装置40が備える各機能部は、一部又は全てが他の装置に備えられてもよい。例えば、構造物評価装置40が備える表示部44が他の装置に備えられてもよい。このように構成される場合、構造物評価装置40は、評価結果を、表示部44を備える他の装置に送信する。表示部44を備える他の装置は、受信した評価結果を表示する。
以上説明した少なくともひとつの実施形態によれば、構造物50から発生した弾性波を検出する複数のセンサ20と、複数のセンサ20それぞれによって検出された複数の弾性波に基づいて、複数の弾性波の発生源の位置を標定する位置標定部423と、構造物の温度に応じて定まる補正値を用いて、位置標定部423による位置標定に基づく情報を補正する補正部425と、補正された情報に基づいて構造物の劣化状態を評価する評価部426とを持つことにより、構造物の劣化状態の評価精度を向上させることすることができる。
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
20、20-1~20-n…センサ,30…信号処理部,40…構造物評価装置,41…通信部,42…制御部,43…記憶部,44…表示部,421…取得部,422…イベント抽出部,423…位置標定部,424…分布生成部,425…補正部,426…評価部,301…波形整形フィルタ,302…ゲート生成回路,303…到達時刻決定部,304…特徴量抽出部,305…送信データ生成部,306…メモリ,307…出力部

Claims (16)

  1. 構造物から発生した弾性波を検出する複数のセンサと、
    前記複数のセンサそれぞれによって検出された複数の弾性波に基づいて、前記複数の弾性波の発生源の位置を標定する位置標定部と、
    前記構造物の温度に応じて定まる補正値を用いて、前記位置標定部による位置標定に基づく情報を補正する補正部と、
    補正された前記情報に基づいて前記構造物の劣化状態を評価する評価部と、
    を備え
    前記構造物は、少なくとも車両が走行可能な舗装部と、前記舗装部を支持する床版部とで構成され、
    前記構造物の温度は、前記舗装部の温度に関連する値である構造物評価システム。
  2. 前記構造物は、温度に対する応答の異なる少なくとも2つの材料で構成される構造を有する、請求項1に記載の構造物評価システム。
  3. 前記位置標定に基づく情報は、密度に関する情報を含み、
    前記補正部は、前記補正値を用いて、前記舗装部の温度が高い場合には密度が高くなる方向へ補正し、前記舗装部の温度が低い場合には密度が低くなる方向へ補正する、請求項1又は2に記載の構造物評価システム。
  4. 前記補正部は、各温度におけるヤング率の値が登録された補正テーブルを参照し、前記補正値を算出する、請求項に記載の構造物評価システム。
  5. 前記補正部は、前記構造物の温度が高い場合に、前記複数の弾性波を検出するための閾値を下げる方向へ補正し、前記構造物の温度が低い場合に前記閾値を上げる方向へ補正する、請求項1又は2に記載の構造物評価システム。
  6. 前記舗装部の温度に関連する値は、気温である、請求項1又は2に記載の構造物評価システム。
  7. 前記舗装部の温度に関連する値は、前記床版部の下面の温度である、請求項1又は2に記載の構造物評価システム。
  8. 前記舗装部の温度に関連する値は、前記舗装部を走行する車両で計測した路面温度である、請求項1又は2に記載の構造物評価システム。
  9. 前記舗装部の温度に関連する値は、日射量である、請求項1又は2に記載の構造物評価システム。
  10. 前記舗装部の温度に関連する値は、原子間距離、ヤング率、体積弾性率、音響インピーダンス、透過率、反射率、弾性波の振幅、弾性波のヒット数、弾性波の標定数のいずれかの値である、請求項1又は2に記載の構造物評価システム。
  11. 前記構造物の温度は、前記舗装部の温度と前記床版部の温度の差、もしくは、比のいずれかに関連する値である、請求項1又は2に記載の構造物評価システム。
  12. 前記舗装部の温度は、気温又は日射量のいずれかにより推定された値、あるいは、前記舗装部の表面の温度であり、前記床版部の温度は、床版下面の温度である、請求項11に記載の構造物評価システム。
  13. 前記舗装部の温度に関連する値は、計測中の温度の統計値である、請求項1又は2に記載の構造物評価システム。
  14. 構造物から発生した弾性波を検出する複数のセンサそれぞれによって検出された複数の弾性波に基づいて、前記複数の弾性波の発生源の位置を標定する位置標定部と、
    前記構造物の温度に応じて定まる補正値を用いて、前記位置標定部による位置標定に基づく情報を補正する補正部と、
    補正された前記情報に基づいて前記構造物の劣化状態を評価する評価部と、
    を備え
    前記構造物は、少なくとも車両が走行可能な舗装部と、前記舗装部を支持する床版部とで構成され、
    前記構造物の温度は、前記舗装部の温度に関連する値である構造物評価装置。
  15. 構造物から発生した弾性波を検出し、
    検出した複数の弾性波に基づいて、前記複数の弾性波の発生源の位置を標定し、
    前記構造物の温度に応じて定まる補正値を用いて、位置標定に基づく情報を補正し、
    補正された前記情報に基づいて前記構造物の劣化状態を評価し、
    前記構造物は、少なくとも車両が走行可能な舗装部と、前記舗装部を支持する床版部とで構成され、
    前記構造物の温度は、前記舗装部の温度に関連する値である構造物評価方法。
  16. コンピュータに、
    構造物から発生した弾性波を検出する複数のセンサそれぞれによって検出された複数の弾性波に基づいて、前記複数の弾性波の発生源の位置を標定する位置標定ステップと、
    前記構造物の温度に応じて定まる補正値を用いて、前記位置標定ステップにおける位置標定に基づく情報を補正する補正ステップと、
    補正された前記情報に基づいて前記構造物の劣化状態を評価する評価ステップと、
    を実行させ
    前記構造物は、少なくとも車両が走行可能な舗装部と、前記舗装部を支持する床版部とで構成され、
    前記構造物の温度は、前記舗装部の温度に関連する値であるためのコンピュータプログラム。
JP2022148052A 2022-09-16 2022-09-16 構造物評価システム、構造物評価装置、構造物評価方法及びコンピュータプログラム Active JP7739245B2 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
JP2022148052A JP7739245B2 (ja) 2022-09-16 2022-09-16 構造物評価システム、構造物評価装置、構造物評価方法及びコンピュータプログラム
CN202310129747.6A CN117723627A (zh) 2022-09-16 2023-02-17 构造物评价系统、装置、方法以及非暂时性记录介质
US18/173,235 US20240094166A1 (en) 2022-09-16 2023-02-23 Structure evaluation system, structure evaluation apparatus, structure evaluation method, and non-transitory computer readable recording medium
EP23158754.4A EP4339578A1 (en) 2022-09-16 2023-02-27 Structure evaluation system, structure evaluation apparatus, structure evaluation method, and computer program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2022148052A JP7739245B2 (ja) 2022-09-16 2022-09-16 構造物評価システム、構造物評価装置、構造物評価方法及びコンピュータプログラム

Publications (2)

Publication Number Publication Date
JP2024043068A JP2024043068A (ja) 2024-03-29
JP7739245B2 true JP7739245B2 (ja) 2025-09-16

Family

ID=85384350

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2022148052A Active JP7739245B2 (ja) 2022-09-16 2022-09-16 構造物評価システム、構造物評価装置、構造物評価方法及びコンピュータプログラム

Country Status (4)

Country Link
US (1) US20240094166A1 (ja)
EP (1) EP4339578A1 (ja)
JP (1) JP7739245B2 (ja)
CN (1) CN117723627A (ja)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001241045A (ja) 2000-03-01 2001-09-04 Natl Inst Of Advanced Industrial Science & Technology Meti 岩盤崩壊の検知方法及び岩盤内亀裂の安定性評価方法
JP2012163561A (ja) 2011-02-08 2012-08-30 Boeing Co:The 構造健全性監視システム
US20140123758A1 (en) 2011-05-20 2014-05-08 Hatch Ltd. Furnace structural integrity monitoring systems and methods
JP2021096266A (ja) 2018-09-13 2021-06-24 株式会社東芝 構造物評価システム、構造物評価装置及び構造物評価方法
WO2022014004A1 (ja) 2020-07-16 2022-01-20 株式会社東芝 構造物評価システム、構造物評価装置及び構造物評価方法
CN114152679A (zh) 2021-10-28 2022-03-08 航天材料及工艺研究所 一种超低温液体环境下的钛合金气瓶声发射二维平面定位方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5242758A (en) * 1975-09-30 1977-04-02 Shimadzu Corp Device for measuring the position of vibration sources and the frequen cy of generation thereof
PL1893972T3 (pl) * 2005-02-22 2017-10-31 Hatch Ltd Układy, sposoby i urządzenie do bezinwazyjnej i nieniszczącej kontroli pieców hutniczych oraz podobnych zbiorników
FR3060743B1 (fr) * 2016-12-15 2019-05-17 Commissariat A L'energie Atomique Et Aux Energies Alternatives Procede et systeme de controle de sante integre d'une structure mecanique par ondes elastiques diffuses
JP7160695B2 (ja) * 2019-01-09 2022-10-25 株式会社東芝 検出システム、検出方法、および情報処理装置
EP4049019A4 (en) * 2019-10-25 2024-03-20 Molex, LLC ULTRASONIC TRANSDUCER PLATE FOR CONDITION MONITORING OF A STRUCTURAL ELEMENT

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001241045A (ja) 2000-03-01 2001-09-04 Natl Inst Of Advanced Industrial Science & Technology Meti 岩盤崩壊の検知方法及び岩盤内亀裂の安定性評価方法
JP2012163561A (ja) 2011-02-08 2012-08-30 Boeing Co:The 構造健全性監視システム
US20140123758A1 (en) 2011-05-20 2014-05-08 Hatch Ltd. Furnace structural integrity monitoring systems and methods
JP2021096266A (ja) 2018-09-13 2021-06-24 株式会社東芝 構造物評価システム、構造物評価装置及び構造物評価方法
WO2022014004A1 (ja) 2020-07-16 2022-01-20 株式会社東芝 構造物評価システム、構造物評価装置及び構造物評価方法
CN114152679A (zh) 2021-10-28 2022-03-08 航天材料及工艺研究所 一种超低温液体环境下的钛合金气瓶声发射二维平面定位方法

Also Published As

Publication number Publication date
JP2024043068A (ja) 2024-03-29
CN117723627A (zh) 2024-03-19
EP4339578A1 (en) 2024-03-20
US20240094166A1 (en) 2024-03-21

Similar Documents

Publication Publication Date Title
JP6368040B2 (ja) 構造物評価システム、構造物評価装置及び構造物評価方法
JP7323698B2 (ja) 構造物評価システム、構造物評価装置及び構造物評価方法
US20180074019A1 (en) Structure evaluation system, structure evaluation apparatus, and structure evaluation method
CN111183358A (zh) 构造物评价系统、构造物评价装置以及构造物评价方法
JP6833770B2 (ja) 構造物評価システム、構造物評価装置及び構造物評価方法
US11519883B2 (en) Structure evaluation system, structure evaluation apparatus, and structure evaluation method
CN114487126B (zh) 车辆信息推测系统、车辆信息推测装置、车辆信息推测方法以及计算机程序
JP7739245B2 (ja) 構造物評価システム、構造物評価装置、構造物評価方法及びコンピュータプログラム
US20180074023A1 (en) Structure evaluation system, structure evaluation apparatus, and structure evaluation method
CN119246673A (zh) 构造物评价系统、构造物评价装置以及构造物评价方法
JP7765008B2 (ja) 構造物評価システム、構造物評価装置及び構造物評価方法
JP7508427B2 (ja) 構造物評価システム、構造物評価装置及び構造物評価方法
US20250027834A1 (en) Structure evaluation system, structure evaluation method, and structure evaluation apparatus
JP2025106735A (ja) 速度推定装置、評価システム及び速度推定方法
JP7757364B2 (ja) 波形信号処理システム、構造物評価システム及び波形信号処理方法
JP2024134786A (ja) 構造物評価システム及び構造物評価方法
JP2025141617A (ja) 損傷検出装置、損傷検出システム及び損傷検出方法

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20240906

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20250423

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20250527

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20250728

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20250805

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20250903

R150 Certificate of patent or registration of utility model

Ref document number: 7739245

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150