(57)【特許請求の範囲】
【請求項１】
被照射体の表面に音波を照射し、その表面の複数の測定箇所において振動速度を測定して、前記被照射体の強度の分布を判定する音波を用いた強度推定方法であって、
音波発信源から音波を照射し、前記被照射体の表面を振動させて、各測定箇所における振動速度を測定する工程と、
得られた振動速度を測定結果に基づいて、各測定箇所における周波数と振動速度との関係を表す振幅スペクトル（Ｓｆ）を求めた後、さらに周波数と振動エネルギー（ＰＳＤ）との関係を求め、得られた周波数と振動エネルギー（ＰＳＤ）との関係について特定範囲で積分値を求め、その積分値が最も小さい値であった測定箇所を最小ＰＳＤ部とする工程と、
下記式（Ｉ）に基づいて、前記最小ＰＳＤ部に対するその他の測定箇所の振動エネルギーの比である振動エネルギー比（ＶＥＲ（１））を求める工程と、
【数１】
0006684074.tif 000008
（式（Ｉ）中のＰＳＤminは、前記最小ＰＳＤ部における、周波数と振動エネルギー（ＰＳＤ）との関係についての積分値（積分範囲：ｆ1〜ｆ2）を意味し、式（Ｉ）中のＰＳＤxは、前記最小ＰＳＤ部以外の測定箇所における、周波数と振動エネルギー（ＰＳＤ）との関係についての積分値（積分範囲：ｆ1〜ｆ2）を意味する。）
求められた振動エネルギー比（ＶＥＲ（１））の値が、第１の閾値以下であった測定箇所を強度が強い箇所であり、第１の閾値以上であった測定箇所を強度が弱い箇所である、と判定する工程と、
を備える強度推定方法。
【請求項２】
さらに、下記式（ＩＩ）に基づいて、前記振幅スペクトル（Ｓｆ）からスペクトルエントロピー（Ｈ）を求める工程と、
【数２】
0006684074.tif 000009
（式（ＩＩ）においてＳｆは振幅スペクトルを意味する。）
前記求められた振動エネルギー比の値が第１の閾値以下であり、かつ、求められたスペクトルエントロピーの値が第２の閾値以上である箇所を、強度が強い箇所と判定し、一方、前記求められた振動エネルギー比の値が第１の閾値以上であり、かつ、求められたスペクトルエントロピーの値が第２の閾値以下である箇所を、強度が弱い箇所と判定する工程と、
を備える、請求項１に記載の強度推定方法。
【請求項３】
被照射体の表面に音波を照射し、その表面の複数の測定箇所において振動速度を測定して、前記被照射体の強度の分布を判定する音波を用いた強度推定方法であって、
音波発信源から音波を照射し、前記被照射体の表面を振動させて、各測定箇所における振動速度を測定する工程と、
得られた振動速度を測定結果に基づいて、各測定箇所における周波数と振動速度との関係を表す振幅スペクトル（Ｓｆ）を求めた後、さらに周波数と振動エネルギー（ＰＳＤ）との関係を求め、得られた周波数と振動エネルギー（ＰＳＤ）との関係について特定範囲で積分値を求め、その積分値が最も小さい値であった測定箇所を最小ＰＳＤ部とする工程と、
下記式（Ｉ）に基づいて、前記最小ＰＳＤ部に対するその他の測定箇所の振動エネルギーの比である振動エネルギー比（ＶＥＲ（１））を求める工程と、
【数３】
0006684074.tif 000010
（式（Ｉ）中のＰＳＤminは、前記最小ＰＳＤ部における、周波数と振動エネルギー（ＰＳＤ）との関係についての積分値（積分範囲：ｆ1〜ｆ2）を意味し、式（Ｉ）中のＰＳＤxは、前記最小ＰＳＤ部以外の測定箇所における、周波数と振動エネルギー（ＰＳＤ）との関係についての積分値（積分範囲：ｆ1〜ｆ2）を意味する。）
下記式（ＩＩ）に基づいて、前記振幅スペクトル（Ｓｆ）からスペクトルエントロピー（Ｈ）を求める工程と、
【数４】
0006684074.tif 000011
（式（ＩＩ）においてＳｆは振幅スペクトルを意味する。）
求められたスペクトルエントロピー（Ｈ）の値が、第２の閾値以下であった測定箇所を強度が弱い箇所であり、第２の閾値以上であった測定箇所を強度が強い箇所である、と判定する工程と、
を備える強度推定方法。
【請求項４】
被照射体の表面に音波を照射し、その表面の複数の測定箇所において振動速度を測定して、前記被照射体の強度の分布を判定する音波を用いた強度推定システムであって、
前記被照射体の表面を振動させ得る音波を発生させる音響発信源と、
前記被照射体の表面の振動速度を測定する計測器と、
得られた振動速度の測定結果を用いて、被照射体の強度の分布を特定する解析装置とを有し、
請求項１〜３のいずれか１に記載の強度推定方法を行うことができる、強度推定システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、音波を用いた強度推定方法および強度推定システムに関する。
【背景技術】
【０００２】
コンクリート強度を推定する手法としては、反発強度法（シュミットハンマー法）が一般的には用いられている。この手法は、コンクリートに打撃を与え、返ってきた衝撃の強さを基に強度を推定する手法である。
【０００３】
しかし、これらの方法では検査対象物に接触して計測する必要があるために、足場を必要とするような高所での適用性に問題があり、遠距離非接触での強度推定法の開発が長い間望まれていた。
【０００４】
一方、空中放射音波を用いる非接触音響検査法として、特許文献１による方法や、本発明者らによる非特許文献１による方法も知られている。
【０００５】
これらの非接触音響検査法は、遠方から、音波あるいは超音波により非検査対象物を加振し、印加された対象物の振動を、遠方において光学的に検出する検査手法である。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開平８−２４８００６号公報
【非特許文献】
【０００７】
【非特許文献１】杉本恒美，歌川紀之，他２名，「SLDVとLRADを用いた非破壊検査のための非接触音響映像法に関する研究」，日本音響学会2012年春季研究発表会講演論文集  1-5-7, March 2012
【発明の概要】
【発明が解決しようとする課題】
【０００８】
本発明者は鋭意研究開発を推し進め、欠陥検査にのみ用いられていた、空中放射音波を用いた非接触音響探査法を、健全部の評価（強度等）にも用いることができることを見出した。
【０００９】
本発明は、非接触で被照射体の健全部の評価をすることができる、音波を用いた強度推定方法および強度推定システムを提供することを目的とする。
【課題を解決するための手段】
【００１０】
本発明者は鋭意検討し、上記課題を解決する方法を見出し、本発明を完成させた。
本発明は次の（ｉ）〜（ｉｖ）である。
（ｉ）被照射体の表面に音波を照射し、その表面の複数の測定箇所において振動速度を測定して、前記被照射体の強度の分布を判定する音波を用いた強度推定方法であって、
音波発信源から音波を照射し、前記被照射体の表面を振動させて、各測定箇所における振動速度を測定する工程と、
得られた振動速度を測定結果に基づいて、各測定箇所における周波数と振動速度との関係を表す振幅スペクトル（Ｓｆ）を求めた後、さらに周波数と振動エネルギー（ＰＳＤ）との関係を求め、得られた周波数と振動エネルギー（ＰＳＤ）との関係について特定範囲で積分値を求め、その積分値が最も小さい値であった測定箇所を最小ＰＳＤ部とする工程と、
下記式（Ｉ）に基づいて、前記最小ＰＳＤ部に対するその他の測定箇所の振動エネルギーの比である振動エネルギー比（ＶＥＲ（１））を求める工程と、
【数１】
0006684074.tif 000002
（式（Ｉ）中のＰＳＤminは、前記最小ＰＳＤ部における、周波数と振動エネルギー（ＰＳＤ）との関係についての積分値（積分範囲：ｆ1〜ｆ2）を意味し、式（Ｉ）中のＰＳＤxは、前記最小ＰＳＤ部以外の測定箇所における、周波数と振動エネルギー（ＰＳＤ）との関係についての積分値（積分範囲：ｆ1〜ｆ2）を意味する。）
求められた振動エネルギー比（ＶＥＲ（１））の値が、第１の閾値以下であった測定箇所を強度が強い箇所であり、第１の閾値以上であった測定箇所を強度が弱い箇所である、と判定する工程と、
を備える強度推定方法。
（ｉｉ）さらに、下記式（ＩＩ）に基づいて、前記振幅スペクトル（Ｓｆ）からスペクトルエントロピー（Ｈ）を求める工程と、
【数２】
0006684074.tif 000003
（式（ＩＩ）においてＳｆは振幅スペクトルを意味する。）
前記求められた振動エネルギー比の値が第１の閾値以下であり、かつ、求められたスペクトルエントロピーの値が第２の閾値以上である箇所を、強度が強い箇所と判定し、一方、前記求められた振動エネルギー比の値が第１の閾値以上であり、かつ、求められたスペクトルエントロピーの値が第２の閾値以下である箇所を、強度が弱い箇所と判定する工程と、
を備える、上記（ｉ）に記載の強度推定方法。
（ｉｉｉ）被照射体の表面に音波を照射し、その表面の複数の測定箇所において振動速度を測定して、前記被照射体の強度の分布を判定する音波を用いた強度推定方法であって、
音波発信源から音波を照射し、前記被照射体の表面を振動させて、各測定箇所における振動速度を測定する工程と、
得られた振動速度を測定結果に基づいて、各測定箇所における周波数と振動速度との関係を表す振幅スペクトル（Ｓｆ）を求めた後、さらに周波数と振動エネルギー（ＰＳＤ）との関係を求め、得られた周波数と振動エネルギー（ＰＳＤ）との関係について特定範囲で積分値を求め、その積分値が最も小さい値であった測定箇所を最小ＰＳＤ部とする工程と、
下記式（Ｉ）に基づいて、前記最小ＰＳＤ部に対するその他の測定箇所の振動エネルギーの比である振動エネルギー比（ＶＥＲ（１））を求める工程と、
【数３】
0006684074.tif 000004
（式（Ｉ）中のＰＳＤminは、前記最小ＰＳＤ部における、周波数と振動エネルギー（ＰＳＤ）との関係についての積分値（積分範囲：ｆ1〜ｆ2）を意味し、式（Ｉ）中のＰＳＤxは、前記最小ＰＳＤ部以外の測定箇所における、周波数と振動エネルギー（ＰＳＤ）との関係についての積分値（積分範囲：ｆ1〜ｆ2）を意味する。）
下記式（ＩＩ）に基づいて、前記振幅スペクトル（Ｓｆ）からスペクトルエントロピー（Ｈ）を求める工程と、
【数４】
0006684074.tif 000005
（式（ＩＩ）においてＳｆは振幅スペクトルを意味する。）
求められたスペクトルエントロピー（Ｈ）の値が、第２の閾値以下であった測定箇所を強度が弱い箇所であり、第２の閾値以上であった測定箇所を強度が強い箇所である、と判定する工程と、
を備える強度推定方法。
（ｉｖ）被照射体の表面に音波を照射し、その表面の複数の測定箇所において振動速度を測定して、前記被照射体の強度の分布を判定する音波を用いた強度推定システムであって、
前記被照射体の表面を振動させ得る音波を発生させる音響発信源と、
前記被照射体の表面の振動速度を測定する計測器と、
得られた振動速度の測定結果を用いて、被照射体の強度の分布を特定する解析装置とを有し、
上記（ｉ）〜（ｉｉｉ）のいずれか１に記載の強度推定方法を行うことができる、強度推定システム。
【発明の効果】
【００１１】
本発明によれば、非接触で被照射体の健全部の評価をすることができる、音波を用いた強度推定方法および強度推定システムを提供することができる。
【図面の簡単な説明】
【００１２】
【図１】本発明の強度推定システムの好適態様を示す概略図である。
【図２】特定の情報処理の代表例を示すフローチャートである。
【図３】コンクリート供試体の外観を示す図である。
【図４】実施例で用いた強度推定システムの写真である。
【図５】コンクリート供試体の拡大図である。
【図６】振動エネルギー比とエントロピースペクトルの計算結果である。
【発明を実施するための形態】
【００１３】
本発明について説明する。
本発明は音波を用いた強度推定方法および強度推定システムである。本発明の強度推定方法および本発明の強度推定システムによれば、被照射体の強度の分布を正確に把握することができる。被照射体としては、例えば、コンクリート構造物、地面（土、砂、石、アスファルト等）、木、液体、人体が挙げられる。
【００１４】
本発明の強度推定方法は、本発明の強度推定システムによって実現することが好ましい。本発明の強度推定システムは、被照射体の表面に音波を照射し、その表面の複数の測定箇所において振動速度を測定して、被照射体の強度の分布を判定する強度推定システムであって、前記被照射体の表面を振動させ得る音波を発生させる音響発信源と、前記被照射体の表面の振動速度を測定する計測器と、得られた振動速度の測定結果を用いて、被照射体の強度の分布を特定する解析装置とを有し、前記解析装置によって特定の情報処理を行うことができる。本発明の強度推定システムとして、具体的には、例えば図１に示す装置が挙げられる。
【００１５】
（強度推定システム）
図１は、被照射体１の表面を振動させ得る音波を発生させる音響発信源１１と、被照射体１の表面の振動速度を測定する計測器１３と、被照射体１の強度の分布を特定するために用いる解析装置１５１を含むコンピュータ１５とを有する装置１０を示す概略図である。図１に示す装置１０は、さらに、任意波形発生装置１７およびアンプ１９を有しており、加えて、コンピュータ１５は制御装置１５２および表示部１５３を含んでおり、制御装置１５２によって任意波形発生装置１７を制御して、所望の周波数の音波を音響発信源１１から発生することができる。任意波形発生装置１７が発生するトリガ信号に制御装置１５２を同期させて計測することもできる。表示部１５３には、後に説明する振動速度分布図等を表示することができる。表示部とはディスプレイ画面等を意味する。
【００１６】
図１に示す本発明の強度推定システム（装置１０）において、音響発信源１１はフラットスピーカである。本発明の強度推定システムにおいて音響発信源の数やスピーカの角度等は特に限定されない。
【００１７】
音響発信源はフラットスピーカの他、パラメトリックスピーカも好ましく用いることができ、また、具体的に、アメリカンテクノロジー社製のＬＲＡＤ（登録商標）を好ましく用いることができる。また、ラウドスピーカを用いることもできるが、この場合は、音響発信源と被照射体との距離を比較的近くする。その他に用いることができる音響発信源としては、パルスレーザ、高圧ガスガン、衝撃波管が挙げられる。
【００１８】
また、音響発信源から被照射体へ照射される音波は、所望の周波数（ω）に調整することができ、かつ、被照射体の表面をその振動速度が計測器によって測定できる程度に、表面に平行方向ではない方向（好ましくは、表面に垂直方向）へ振動させることができる音波であればよく、空気中で振動振幅が減衰し難い可聴帯域の音波（音響波）が好ましい。なお、超音波は用い難い。超音波は空気中で振動振幅の減衰が大きいからである。
また、被照射体の共振周波数帯が不明な場合には、音響発信源から被照射体へ照射される音波は、ホワイトノイズであることが好ましい。全ての周波数を含んでいるからである。
【００１９】
音響発信源から被照射体へ音波を照射することで、被照射体の表面に９０ｄＢ以上の音圧を発生させることが好ましく、１００ｄＢ程度の音圧を発生させることがより好ましい。
【００２０】
図１に示す本発明の強度推定システム（装置１０）において、計測器１３はレーザドップラー振動計であることが好ましく、レーザ１３１を被照射体１に照射して、その表面の振動速度を測定することができる。得られた振動速度のデータは解析装置１５１で解析するために用いられる。
なお、本発明の強度推定システムにおいて計測器は、被照射体の表面の振動速度を非接触で測定できるものであれば特に限定されず、例えばレーザ変位計を用いることができ、レーザドップラー振動計であることが好ましい。被照射体と計測器とが比較的離れていても、被照射体の表面の振動を正確に測定することができるからである。
また、１度に１点の振動計測が可能なシングルレーザタイプのレーザ振動計を用いることは可能であるが、スキャニングレーザタイプのレーザ振動計を用いることが好ましい。スキャニング振動計であるレーザドップラー振動計としては、具体的に、ポリテックジャパン社製のＰＳＶ４００−Ｈ４が挙げられる。このレーザドップラー振動計は解析装置の一部および制御装置を含むものである。
【００２１】
図１に示す本発明の強度推定システム（装置１０）において、解析装置１５１は、被照射体１における強度の分布を特定するための特定の情報処理を行うことができるものであれば特に限定されない。この特定の情報処理は本発明の強度推定方法が備えるものであり、後に詳細に説明する。
【００２２】
図１に示す本発明の強度推定システム（装置１０）において、任意波形発生装置１７は、制御装置１５２の指令によって所望の周波数の音波を音響発信源１１から発生させることができる装置である。例えば、ノイズ波やバースト波を発生可能な市販のファンクションジェネレータ等を用いることができる。送信する音波の波形は通常この任意波形発生装置により制御することができる。通常は簡単のために手動で制御するが、解析装置側から制御するようにシステムを構成することも可能である。任意波形発生装置１７が発生するトリガ信号に制御装置１５２を同期させて計測することもできる。
また、アンプ１９は特に限定されず、例えば、市販オーディオアンプ等を用いることができる。
【００２３】
本発明の強度推定システムは、本発明の強度推定方法を実施することができる構成を備えている。
【００２４】
（強度推定方法）
次に、本発明の強度推定方法について説明する。
本発明の強度推定方法は、被照射体の表面に音波を照射し、その表面の複数の測定箇所において振動速度を測定して、被照射体の強度の分布を判定する音波を用いた強度推定方法である。
そして、音波発信源から音波を照射し、前記被照射体の表面を振動させて、各測定箇所における振動速度を測定する工程を備える。
この工程は、例えば前述の本発明の強度推定システムを用い、音響発信源から被照射体へホワイトノイズを照射し、レーザドップラー振動計などの計測器を用いて、その表面の各測定箇所における振動速度を測定して行うことができる。振動速度は測定時間（時刻）との関係として得ることができる。
【００２５】
なお、ここでの計測結果についてゲート処理を施して、目的信号のみを抽出することが好ましい。ゲート処理とは、時間的、周波数的に計測したい信号を取り出す処理であり、本発明者が既に行った特願２０１２−２５８８８８号に記載の処理が例示される。
【００２６】
本発明の強度推定方法は、さらに、以下に説明する特定の情報処理を行う各工程を備える。
【００２７】
本発明の強度推定方法は、特定の情報処理を行う工程として、前工程によって得られた振動速度を測定結果に基づいて、各測定箇所における周波数と振動速度との関係を表す振幅スペクトル（Ｓｆ）を求めた後、さらに周波数と振動エネルギー（ＰＳＤ）との関係を求め、得られた周波数と振動エネルギー（ＰＳＤ）との関係について特定範囲で積分値を求め、その積分値が最も小さい値であった測定箇所を最小ＰＳＤ部とする工程を備える。
【００２８】
この工程では、前工程によって得られた振動速度の測定結果をフーリエ変換して振幅スペクトル（Ｓｆ）を求め、さらに周波数と振動エネルギー（ＰＳＤ）との関係を求める。
振動エネルギー（ＰＳＤ）は、振動速度の２乗に比例する値である。
【００２９】
また、レーザドップラー振動計などの計測器に共振周波数が存在する場合および多重反射による影響が無視できない場合には、一度、フーリエ変換後、周波数領域で指定範囲内のエネルギーを計算することが好ましい。
【００３０】
次に、得られた周波数と振動エネルギー（ＰＳＤ）との関係について、特定範囲で積分値を求める。
積分する範囲は特に限定されないが、計測器の共振周波数を含まない範囲で積分することが好ましい。例えばレーザドップラー振動計等の共振周波数のノイズが存在する場合、そのノイズを含まない範囲を積分する範囲としてもよい。具体的には、例えばレーザドップラー振動計等の共振周波数のノイズが１ｋＨｚ以下に存在する場合、１３５０Ｈｚ〜８０００Ｈｚの範囲を積分する範囲としてもよい。
【００３１】
このようにして、各測定箇所について、上記の積分値が得られる。そして、その積分値が最も小さい値であった測定箇所を、最小ＰＳＤ部とする。
【００３２】
次に、本発明の強度推定方法は、特定の情報処理を行う工程として、前記最小ＰＳＤ部に対するその他の測定箇所の振動エネルギーの比である振動エネルギー比（ＶＥＲ（１））を求める。
【００３３】
振動エネルギー比（ＶＥＲ（１））は、前記最小ＰＳＤ部に対するその他の測定箇所の振動エネルギーの比であり、次の式（Ｉ）から算出される。
【００３４】
【数５】
0006684074.tif 000006
【００３５】
すなわち、式（Ｉ）における分母（式（Ｉ）中の「ＰＳＤmin」）が、前記最小ＰＳＤ部における、周波数と振動エネルギー（ＰＳＤ）との関係についての積分値（積分範囲：ｆ1〜ｆ2）である。そして、式（Ｉ）における分子（式（Ｉ）中の「ＰＳＤx」）が、前記最小ＰＳＤ部以外の測定箇所における、周波数と振動エネルギー（ＰＳＤ）との関係についての積分値（積分範囲：ｆ1〜ｆ2）である。
【００３６】
そして、求められた振動エネルギー比（ＶＥＲ（１））の値が、第１の閾値以下であった測定箇所を強度が強い（硬い）箇所であり、その値が第１の閾値以上であった測定箇所を強度が弱い（柔らかい）箇所である、と判定する。
【００３７】
第１の閾値は、被照射体の材質等によって概ね決まる値である。例えば、非照射体がコンクリートである場合、コンクリートの配合や相対的な評価にもよるが、第１の閾値は、１．２〜１．５程度である。
【００３８】
次に、本発明の強度推定方法は、特定の情報処理を行う工程として、前記振幅スペクトル（Ｓｆ）からスペクトルエントロピー（Ｈ）を求める。
【００３９】
各測定箇所におけるスペクトルエントロピー（Ｈ）は、次の式（ＩＩ）より算出する。
【００４０】
【数６】
0006684074.tif 000007
【００４１】
ここで、Ｓｆは振幅スペクトルである。また、スペクトルエントロピー（Ｈ）の周波数計算範囲は特に限定されないが、前述の、周波数と振動エネルギー（ＰＳＤ）との関係について積分値を求める範囲と同一とすることが好ましい。
【００４２】
そして、求められたスペクトルエントロピー（Ｈ）の値が、第２の閾値以下であった測定箇所を強度が弱い（柔らかい）箇所であり、その値が第２の閾値以上であった測定箇所を強度が強い（硬い）箇所である、と判定する。
【００４３】
第２の閾値は、被照射体の材質等によって概ね決まる値である。例えば、非照射体がコンクリートである場合、コンクリートの配合や相対的な評価にもよるが、第２の閾値は、１２．５程度である。
【００４４】
強度の判定については、上述したように、振動エネルギー比の値のみで判断してもよく、スペクトルエントロピーの値のみで判断してもよい。さらに、振動エネルギー比とスペクトルエントロピーの値を複合的に判定してもよい。すなわち、振動エネルギー比の値が第１の閾値以下であり、かつ、スペクトルエントロピーの値が第２の閾値以上である箇所を、強度が強い（硬い）箇所と判定し、一方、振動エネルギー比の値が第１の閾値以上であり、かつ、スペクトルエントロピーの値が第２の閾値以下である箇所を、強度が弱い（柔らかい）箇所と判定してもよい。
【００４５】
上記のような一連の処理についての代表例を図２にフローチャートとして表す。
図２に示すように、最初に、各計測点結果についてゲート処理（時間、周波数）を施し、目的信号のみを抽出する。次に振動エネルギーを計算する（レーザドップラー振動計に共振周波数が存在する場合には一度ＦＦＴ（フーリエ変換）後周波数領域で指定範囲内のエネルギーを計算することが好ましい。）。次に計測点中の最も低い値を示したエネルギーを健全部の基準とし、振動エネルギー比を計算する。
【００４６】
次に、前記振幅スペクトル（Ｓｆ）からスペクトルエントロピー（Ｈ）を求める。
【００４７】
次に、振動エネルギー比（ＶＥＲ（１））の値が第１の閾値以下であった測定箇所を強度が強い（硬い）箇所であり、その値が閾値以上であった測定箇所を強度が弱い（柔らかい）箇所である、と判定する。また、スペクトルエントロピー（Ｈ）の値が第２の閾値以下であった測定箇所を強度が弱い（柔らかい）箇所であり、その値が第２の閾値以上であった測定箇所を強度が強い（硬い）箇所である、と判定する。
【００４８】
このとき、振動エネルギー比（ＶＥＲ（１））の値が第１の閾値以下であり、かつ、スペクトルエントロピー（Ｈ）の値が第２の閾値以上である箇所を強度が強い（硬い）箇所であると判定してもよい。また、振動エネルギー比（ＶＥＲ（１））の値が第１の閾値以上であり、かつ、スペクトルエントロピー（Ｈ）の値が第２の閾値以下である箇所を強度が弱い（柔らかい）箇所であると判定してもよい。
【００４９】
このようにして、被照射体の健全部についての強度の分布を判定する。
具体的な方法は、後述する実施例において説明する。
【００５０】
次に、本発明者が、上記のような振動エネルギー比及びスペクトルエントロピーを用いた強度推定方法に関する本発明に至った経緯について記す。
【００５１】
非接触音響探査法では、空中放射音波によりコンクリート壁面を振動させ、この振動をレーザ計測により、振動速度スペクトルとして計測する。欠陥部と比べると健全部はほとんど振動していないことから、従来は健全部で計測される振動速度スペクトルの物理的な意味については考慮してこなかった。しかしながら、同じ健全部でも振動速度スペクトルにはばらつきがあり、特に強度が明瞭に異なる場合にはそのスペクトルの分布にも差が出ることが考えられる。例えば、硬い部分と柔らかい部分が同じコンクリート壁面に存在したとすると、柔らかい部分は硬い部分に比べると振動しやすいため、振動エネルギー（比）として比較した場合には、柔らかい部分の振動エネルギーは高くなることが想定される。また、このことはスペクトルの分布にも影響し、硬い部分は柔らかい部分に比べると平坦なスペクトル分布特性、すなわち白色度の高い特性を持っていることが想定される。振動エネルギー（比）とスペクトルエントロピーは元来、計測時の受光漏れに起因する計測不良点の識別のために考案されたものであるが、実はこのような対象物の強度推定にも用いることができる。
【実施例】
【００５２】
初めに、強度が異なる４種類の円筒状コンクリートを内包した供試体を用意した。円筒状コンクリートの半径は、約５０ｍｍである。この供試体の寸法図を図３に示す。
【００５３】
次に、図４に示すように供試体を配置した後、長距離音響放射装置（ＬＲＡＤ）から供試体へ音波を照射して、供試体表面に１００ｄＢの音圧を発生させた。そして、スキャニング振動計（ＳＬＤＶ）を用いて、その表面の複数個所について振動速度を測定した。測定位置は、図５に示すとおりである。強度が異なる円筒状コンクリートを、左から２、５、８、１１として、測定した。その他の１、３、４、６、７、９、１０、１２は、周囲のコンクリートの測定である。なお、これら４つの円筒状コンクリートの圧縮強度は、測定箇所２、５、８、１１について、それぞれ３０、４０、６０、８０Ｎ／ｍ2である。
【００５４】
この際の測定条件は、以下のとおりである。
計測周波数範囲：５００−７１００Ｈｚ
音源：ＬＲＡＤ（Ｌｏｎｇ  Ｒａｎｇｅ  Ａｃｏｕｓｔｉｃ  Ｄｅｖｃｅ）
ＬＲＡＤから測定面までの距離：２．５４９ｍ
ＳＬＤＶから測定面までの距離：３．１６５ｍ
音波：トーンバースト波（変調周波数２００Ｈｚ，ｉｎｔｅｒｖａｌ＝５０ｍｓ）
測定点での加算平均回数：５回
測定面での最大音圧は：１００ｄＢ程度
【００５５】
次に、各測定箇所（１２か所）の各々における振動速度の測定結果（振動速度と時間との関係）をフーリエ変換し、測定箇所ごとの振幅スペクトル（Ｓｆ）を求めた。
【００５６】
次に、得られた振幅スペクトル（Ｓｆ）に基づいて、周波数と振動エネルギー（ＰＳＤ）との関係を求めた。振動エネルギー（ＰＳＤ）は、振動速度の２乗に比例する値である。
【００５７】
そして、得られた周波数と振動エネルギー（ＰＳＤ）との関係において、ｆ1＝１３５０Ｈｚからｆ2＝８１９２Ｈｚまでの範囲で積分値（ＰＳＤ積分値）を求めた。
その結果、１番のＰＳＤ積分値が最も小さくなったため、１番を「最小ＰＳＤ部」と認定した。そして、他の測定箇所（１１か所）について、１番と対比する振動エネルギー比（ＶＥＲ（１））を求めた。すなわち、前述の式（Ｉ）における分母を１番のものとし、他の各測定箇所におけるものを分子として、各々における振動エネルギー比（ＶＥＲ（１））を求めた。
【００５８】
次に、各測定箇所における前述の式（ＩＩ）に基づいてスペクトルエントロピーを計算した。なお、スペクトルエントロピー（Ｈ）の周波数計算範囲は１３５０〜８１９２Ｈｚとした。
【００５９】
このようにして得られた、１〜１２番までの測定箇所の振動エネルギー比とエントロピースペクトルの計算結果を図６に示す。
振動エネルギー比の第１の閾値を１．４程度とし、エントロピースペクトルの第２の閾値を１２．５程度とすると、測定箇所２及び５は、振動エネルギー比が第１の閾値以上であり、エントロピースペクトルが第２の閾値以下であるので、強度が弱い柔らかいグループに分類される。一方、測定箇所８及び１１は、振動エネルギー比が第１の閾値以下であり、エントロピースペクトルが第２の閾値以上であるので、強度が強い硬いグループに分類される。
【００６０】
このように、非接触音響探査法の欠陥検出アルゴリズムを健全部に適用することで、コンクリートの強度推定に用いることができることが明らかになった。従来、健全部でのデータのばらつきには物理的な意味があるとは考えていなかったが、今回、初めて意味があることが示された。なお、本手法はコンクリート以外のものでも使用できるため、今後さらに適用範囲が広がっていくことが予想される。
【符号の説明】
【００６１】
１    被照射体
１０  装置
１１  音響発信源
１３  計測器
１３１  レーザ
１５  コンピュータ
１５１  解析装置
１５２  制御装置
１５３  表示部
１７  任意波形発生装置
１９  アンプ
【図１】
0006684074.tif 000012
【図２】
0006684074.tif 000013
【図３】
0006684074.tif 000014
【図４】
0006684074.tif 000015
【図５】
0006684074.tif 000016
【図６】
0006684074.tif 000017