JP5629481B2 - 損傷診断システム - Google Patents
損傷診断システム Download PDFInfo
- Publication number
- JP5629481B2 JP5629481B2 JP2010058784A JP2010058784A JP5629481B2 JP 5629481 B2 JP5629481 B2 JP 5629481B2 JP 2010058784 A JP2010058784 A JP 2010058784A JP 2010058784 A JP2010058784 A JP 2010058784A JP 5629481 B2 JP5629481 B2 JP 5629481B2
- Authority
- JP
- Japan
- Prior art keywords
- mode
- frequency
- propagation
- subject
- damage
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating 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/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
- G01N29/4409—Processing the detected response signal, e.g. electronic circuits specially adapted therefor by comparison
- G01N29/4436—Processing the detected response signal, e.g. electronic circuits specially adapted therefor by comparison with a reference signal
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating 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/04—Analysing solids
- G01N29/07—Analysing solids by measuring propagation velocity or propagation time of acoustic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating 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/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
- G01N29/46—Processing the detected response signal, e.g. electronic circuits specially adapted therefor by spectral analysis, e.g. Fourier analysis or wavelet analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/01—Indexing codes associated with the measuring variable
- G01N2291/011—Velocity or travel time
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/028—Material parameters
- G01N2291/0289—Internal structure, e.g. defects, grain size, texture
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/04—Wave modes and trajectories
- G01N2291/042—Wave modes
- G01N2291/0427—Flexural waves, plate waves, e.g. Lamb waves, tuning fork, cantilever
Landscapes
- Physics & Mathematics (AREA)
- Immunology (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Biochemistry (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Acoustics & Sound (AREA)
- Mathematical Physics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Investigating Materials By The Use Of Optical Means Adapted For Particular Applications (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Description
本発明は、ラム波を利用した損傷診断システムに関する。
例えば、航空機等の機体のような、素材に対する強度と軽量化の双方が要求される分野においては、このような要求に応えるために、CFRP等の複合材料の大幅な適用化が不可欠である。そして、その複合材料からなる構造の高信頼性の確保とより効率的な設計のために、損傷探知技術（ヘルスモニタリング技術）が注目を集めている。
このような複合材料の損傷、欠陥等の探知を行う装置として、特許文献１，２には、ＦＢＧ（Fiber Bragg Grating）光ファイバセンサを用いた損傷探知装置が記載されている。光ファィバは、昨今では細径化（例えば、直径52[μm]）が進んでおり、構造物に埋め込んでも、当該構造物の強度の低下をあまり生じないため、その設置に関して自由度が高いという利点を備えている。
このような複合材料の損傷、欠陥等の探知を行う装置として、特許文献１，２には、ＦＢＧ（Fiber Bragg Grating）光ファイバセンサを用いた損傷探知装置が記載されている。光ファィバは、昨今では細径化（例えば、直径52[μm]）が進んでおり、構造物に埋め込んでも、当該構造物の強度の低下をあまり生じないため、その設置に関して自由度が高いという利点を備えている。
特許文献１，２記載の発明によれば、構造用複合材料の所定個所に固定配置されたピエゾ素子と、ピエゾ素子に信号を伝達する導線と、ピエゾ素子との間に構造用複合材料を構成する複合材料を挟んで固定配置されコア部に所定の波長光を反射するグレーティング部を有する光ファイバセンサと、コア部に光照射を行う光源と、グレーティング部からの反射光の特性を検出する特性検出手段とを用い、ピエゾ素子により加振し特性検出手段の出力の変化から損傷を探知する。また、特性検出手段としては、グレーティング部からの反射光の周波数特性を検出するスペクトラムアナライザが用いられる。
さらに、特許文献１記載の発明には、予め取得した正常な構造用複合材料による検出データとの比較を行うか、他の方法として、スペクトラムアナライザの検出する周波数分布による特定周波数の非振動時からの変動値に閾値を設定し、それ以下の場合は損傷ありと判定しても良い旨記載されている（段落0032）。
特許文献２には、スペクトラムアナライザに２つの光学フィルタが設けられ、スペクトラムアナライザが反射光を前記２つの光学フィルタを通して演算処理装置に出力することで反射光の波長振動信号を高感度に検出すること、演算処理装置が得られた波長振動信号に基づき被検体の損傷の規模に相当する値（DI値）を算出することが提案された。
損傷探知技術の一手法として、ラム（Lamb）波と呼ばれる形態の超音波を送受振し、受振波形の変化から損傷発生を診断する研究が行われている。ラム波とは薄板を伝播する超音波であり、減衰が比較的小さく長距離に渡って板構造を伝播するため、損傷探知に適した超音波の伝播形態である。
また、ラム波には多重モード性と速度分散性（周波数依存性）という二つの特徴があり、板厚と周波数に依存して、速度の異なる複数のモードが存在する。この複雑な特徴のため、従来では、ラム波の特定の周波数の情報のみを利用して、損傷検知が行われてきた。
また、ラム波には多重モード性と速度分散性（周波数依存性）という二つの特徴があり、板厚と周波数に依存して、速度の異なる複数のモードが存在する。この複雑な特徴のため、従来では、ラム波の特定の周波数の情報のみを利用して、損傷検知が行われてきた。
本発明は、このラム波の分散性という特徴を利用して、広帯域にわたるモード分散性を計測可能とし、従来よりも損傷検知に役立つより多くの情報が得られ、剥離長さ定量評価を可能とし、高精度で高信頼性の高い損傷の探知、診断を可能にする損傷診断システムを提供することを課題とする。
以上の課題を解決するための請求項１記載の発明は、被検体に広帯域ラム波の超音波振動を加える加振装置と、
前記被検体から伝達される広帯域ラム波を検知する振動検知センサと、
前記加振装置の発振動作を制御し、前記振動検知センサの出力値を処理して計測結果を出力する処理装置とを備え、
前記処理装置は、
前記加振装置による加振時に得られる前記振動検知センサの出力値の出力値を変換して周波数及び伝播時間の２次元に展開した伝播強度分布データを得て、ラム波の非対称モードを同定し、同定したラム波の非対称モードの中から選ばれる１つ又は２つ以上の特定の非対称モードについて、周波数ごとに伝播強度の最大値が生じる伝播時間を算出し、当該最大値が生じる周波数及び伝播時間の関係を特定し、
特定した前記周波数及び伝播時間の関係と、記憶手段に蓄積されている、損傷状態が既知の被検体において予め特定した前記周波数及び伝播時間の関係と対比して、当該被検体に対する損傷規模の推定演算を行って表示する損傷診断システムである。
前記被検体から伝達される広帯域ラム波を検知する振動検知センサと、
前記加振装置の発振動作を制御し、前記振動検知センサの出力値を処理して計測結果を出力する処理装置とを備え、
前記処理装置は、
前記加振装置による加振時に得られる前記振動検知センサの出力値の出力値を変換して周波数及び伝播時間の２次元に展開した伝播強度分布データを得て、ラム波の非対称モードを同定し、同定したラム波の非対称モードの中から選ばれる１つ又は２つ以上の特定の非対称モードについて、周波数ごとに伝播強度の最大値が生じる伝播時間を算出し、当該最大値が生じる周波数及び伝播時間の関係を特定し、
特定した前記周波数及び伝播時間の関係と、記憶手段に蓄積されている、損傷状態が既知の被検体において予め特定した前記周波数及び伝播時間の関係と対比して、当該被検体に対する損傷規模の推定演算を行って表示する損傷診断システムである。
請求項２記載の発明は、前記加振装置および前記振動検知センサは、前記被検体の表裏の同位置にそれぞれ設置され、
前記処理装置は、前記表裏の前記加振装置に逆位相の波を発振させることにより非対称に発振させ、表裏の前記振動検知センサの出力値同士を加算して算出したデータを変換して周波数及び伝播時間の２次元に展開した前記伝播強度分布データを得る請求項１に記載の損傷診断システムである。
前記処理装置は、前記表裏の前記加振装置に逆位相の波を発振させることにより非対称に発振させ、表裏の前記振動検知センサの出力値同士を加算して算出したデータを変換して周波数及び伝播時間の２次元に展開した前記伝播強度分布データを得る請求項１に記載の損傷診断システムである。
請求項３記載の発明は、前記処理装置は、
前記伝播強度分布データの最大値が生じる周波数及び伝播時間の関係を特定する代わりに、
前記周波数ごとに伝播強度の最大値が生じる伝播時間を算出し特定した周波数及び伝播時間の関係を示すデータを、横軸が周波数、縦軸が伝播時間で表示されるグラフ上にプロットし、ある周波数範囲における当該データの傾きを変化率として算出し、
算出した前記変化率と、記憶手段に蓄積されている、損傷状態が既知の被検体において予め算出した前記変化率と対比して、当該被検体に対する損傷規模の推定演算を行って表示する請求項１又は請求項２に記載の損傷診断システムである。
前記伝播強度分布データの最大値が生じる周波数及び伝播時間の関係を特定する代わりに、
前記周波数ごとに伝播強度の最大値が生じる伝播時間を算出し特定した周波数及び伝播時間の関係を示すデータを、横軸が周波数、縦軸が伝播時間で表示されるグラフ上にプロットし、ある周波数範囲における当該データの傾きを変化率として算出し、
算出した前記変化率と、記憶手段に蓄積されている、損傷状態が既知の被検体において予め算出した前記変化率と対比して、当該被検体に対する損傷規模の推定演算を行って表示する請求項１又は請求項２に記載の損傷診断システムである。
請求項４記載の発明は、前記処理装置は、
前記変化率を算出する代わりに、
前記伝播強度分布データの最大値における伝播時間の、前記被検体に損傷がない場合に対する減少量を算出し、
算出した前記減少量に基づいて、当該被検体に対する損傷規模の推定演算を行って表示する請求項３に記載の損傷診断システムである。
前記変化率を算出する代わりに、
前記伝播強度分布データの最大値における伝播時間の、前記被検体に損傷がない場合に対する減少量を算出し、
算出した前記減少量に基づいて、当該被検体に対する損傷規模の推定演算を行って表示する請求項３に記載の損傷診断システムである。
請求項５記載の発明は、被検体に広帯域ラム波の超音波振動を加える加振装置と、
前記被検体から伝達される広帯域ラム波を検知する振動検知センサと、
前記加振装置の発振動作を制御し、前記振動検知センサの出力値を処理して計測結果を出力する処理装置とを備え、
前記処理装置は、
前記加振装置による加振時に得られる前記振動検知センサの出力値の出力値を変換して周波数及び伝播時間の２次元に展開した伝播強度分布データを得て、ラム波の対称モードを同定し、同定したラム波の対称モードの中から選ばれる１つ又は２つ以上の特定の対称モードについて、周波数ごとに伝播強度の最大値が生じる伝播時間を算出し、当該最大値が生じる周波数及び伝播時間の関係を特定し、
特定した前記周波数及び伝播時間の関係と、記憶手段に蓄積されている、損傷状態が既知の被検体において予め特定した前記周波数及び伝播時間の関係と対比して、当該被検体に対する損傷規模の推定演算を行って表示する損傷診断システムである。
前記被検体から伝達される広帯域ラム波を検知する振動検知センサと、
前記加振装置の発振動作を制御し、前記振動検知センサの出力値を処理して計測結果を出力する処理装置とを備え、
前記処理装置は、
前記加振装置による加振時に得られる前記振動検知センサの出力値の出力値を変換して周波数及び伝播時間の２次元に展開した伝播強度分布データを得て、ラム波の対称モードを同定し、同定したラム波の対称モードの中から選ばれる１つ又は２つ以上の特定の対称モードについて、周波数ごとに伝播強度の最大値が生じる伝播時間を算出し、当該最大値が生じる周波数及び伝播時間の関係を特定し、
特定した前記周波数及び伝播時間の関係と、記憶手段に蓄積されている、損傷状態が既知の被検体において予め特定した前記周波数及び伝播時間の関係と対比して、当該被検体に対する損傷規模の推定演算を行って表示する損傷診断システムである。
請求項６記載の発明は、前記加振装置および前記振動検知センサは、前記被検体の表裏の同位置にそれぞれ設置され、
前記処理装置は、前記表裏の前記加振装置に同位相の波を発振させることにより対称に発振させ、表裏の前記振動検知センサの出力値同士を加算して算出したデータを変換して周波数及び伝播時間の２次元に展開した伝播強度分布データを得る請求項５に記載の損傷診断システムである。
前記処理装置は、前記表裏の前記加振装置に同位相の波を発振させることにより対称に発振させ、表裏の前記振動検知センサの出力値同士を加算して算出したデータを変換して周波数及び伝播時間の２次元に展開した伝播強度分布データを得る請求項５に記載の損傷診断システムである。
請求項７記載の発明は、前記処理装置は、
前記伝播強度分布データの最大値が生じる周波数及び伝播時間の関係を特定する代わりに、
前記伝播強度分布データの最大値における伝播時間の、前記被検体に損傷がない場合に対する増加量を算出し、
算出した前記増加量に基づいて、当該被検体に対する損傷規模の推定演算を行って表示する請求項５又は請求項６に記載の損傷診断システムである。
前記伝播強度分布データの最大値が生じる周波数及び伝播時間の関係を特定する代わりに、
前記伝播強度分布データの最大値における伝播時間の、前記被検体に損傷がない場合に対する増加量を算出し、
算出した前記増加量に基づいて、当該被検体に対する損傷規模の推定演算を行って表示する請求項５又は請求項６に記載の損傷診断システムである。
請求項８記載の発明は、前記処理装置は、前記推定演算を行って得られた前記被検体に対する損傷規模を表示する代わりに、特定した非対称モードまたは対称モードにおける前記周波数及び伝播時間の関係、算出した前記伝播時間の変化率、算出した前記減少量、または算出した前記増加量を表示する請求項１から請求項７のいずれか一項に記載の損傷診断システムである。
ラム波は、板状の振動伝播媒体の板厚中心に対して対称的な変位を持つ対称モード (symmetric modes; Ｓモード)と非対称な変位を持つ非対称モード(a symmetric modes; Ａモード)があり、基本波対称モード（Ｓ０）と基本波非対称モード（Ａ０）のそれぞれに対し高次なｎ次モード（Ｓｎ，Ａｎ）が無数に存在し、波動が複雑になる。
本発明者らの研究において、広帯域ラム波を発振、受振することにより、対称/非対称モードを分離する手法が構成され、これを利用して各モードについて分析した結果、Ｓ１モードが層間剥離部においてＳ０モード及びＡ１モードに変換されて伝播し、層間剥離部を通過すると再びＳ１モードに戻って伝播することが分かった。
また、Ａ１モードが層間剥離部においてＡ１モードより伝播速度の速いＳ０モードに変換されて伝播し、層間剥離部を通過すると再びＡ１モードに戻って伝播することが分かった。
すなわち、層間剥離部における速度変化が到達時間の変化となり、層間剥離部の振動伝播方向の長さにより、各モードの到達時間はそれぞれ独自の変化を示すことがわかった。
したがって、本発明によれば、周波数及び伝播時間の２次元に展開した伝播強度分布データを得、特定のモードについて前記データから損傷の影響による本モードの到達時間の変化を顕示する所定の特徴値（損傷規模の指標）を得ることで、損傷の有無や規模を診断することが可能となるという効果がある。
本発明者らの研究において、広帯域ラム波を発振、受振することにより、対称/非対称モードを分離する手法が構成され、これを利用して各モードについて分析した結果、Ｓ１モードが層間剥離部においてＳ０モード及びＡ１モードに変換されて伝播し、層間剥離部を通過すると再びＳ１モードに戻って伝播することが分かった。
また、Ａ１モードが層間剥離部においてＡ１モードより伝播速度の速いＳ０モードに変換されて伝播し、層間剥離部を通過すると再びＡ１モードに戻って伝播することが分かった。
すなわち、層間剥離部における速度変化が到達時間の変化となり、層間剥離部の振動伝播方向の長さにより、各モードの到達時間はそれぞれ独自の変化を示すことがわかった。
したがって、本発明によれば、周波数及び伝播時間の２次元に展開した伝播強度分布データを得、特定のモードについて前記データから損傷の影響による本モードの到達時間の変化を顕示する所定の特徴値（損傷規模の指標）を得ることで、損傷の有無や規模を診断することが可能となるという効果がある。
以下に本発明の一実施形態につき図面を参照して説明する。以下は本発明の一実施形態であって本発明を限定するものではない。
〔基本構成〕
まず、本実施形態の損傷探知システムの基本構成につき説明する。
図１は、構造用複合材料Ｚの損傷探知を行う損傷探知システム１０の概略構成図である。本実施形態では、構造用複合材料を被検体とする。
本実施形態では、被検体にラム波の超音波振動を加える加振装置として、ＭＦＣ（Macro Fiber Composite）アクチュエータを用いる。ＭＦＣアクチュエータは、極細の角柱状の圧電セラミックスを一方向に並べてエポキシ樹脂に埋め、その上下面に電極を接着したものであり、面内の一方向に比較的大きなひずみを発生させることが出来る。その特性から超音波発振素子としても使用可能であることが知られている。加振装置としてピエゾ素子などの他の発振アクチュエータを適用してもよい。
まず、本実施形態の損傷探知システムの基本構成につき説明する。
図１は、構造用複合材料Ｚの損傷探知を行う損傷探知システム１０の概略構成図である。本実施形態では、構造用複合材料を被検体とする。
本実施形態では、被検体にラム波の超音波振動を加える加振装置として、ＭＦＣ（Macro Fiber Composite）アクチュエータを用いる。ＭＦＣアクチュエータは、極細の角柱状の圧電セラミックスを一方向に並べてエポキシ樹脂に埋め、その上下面に電極を接着したものであり、面内の一方向に比較的大きなひずみを発生させることが出来る。その特性から超音波発振素子としても使用可能であることが知られている。加振装置としてピエゾ素子などの他の発振アクチュエータを適用してもよい。
図１に示すように、本実施形態の損傷探知システム１０は、構造用複合材料Ｚの損傷探知を行うべき箇所の近傍において構造用複合材料Ｚの面部に貼着されるＭＦＣアクチュエータ２１と、構造用複合材料Ｚの損傷探知を行うべき箇所の近傍に設置される振動検知センサとしての光ファイバセンサ３０と、ＭＦＣアクチュエータ２１の制御装置４１と、光ファイバセンサ３０から得られる反射光の波長特性を検出するスペクトラムアナライザ４２と、スペクトラムアナライザ４２の出力値を演算処理する演算処理装置５０とを備えている。なお、スペクトラムアナライザ４２の電源装置４３を図示した。
ＭＦＣアクチュエータ２１は、外部から駆動電圧を印加されると、その面内の一方向に比較的大きなひずみを発生させる。かかる性質を利用して、制御装置４１は、ＭＦＣアクチュエータ２１に対して駆動用の電圧を印加し、構造用複合材料Ｚに瞬間的な振動を加える。
光ファイバセンサ３０は、ＦＢＧ（Fiber Bragg Grating）光ファイバセンサであり、図２(a)の概略構成図に示すように、コア部３２に所定の波長光を反射するグレーティング部３３を有する光ファイバ３４から構成されている。
光ファイバ３４は、その一端部においてスペクトラムアナライザ４２に接続されており、当該スペクトラムアナライザ４２が有する光源により、所定範囲の波長帯域を網羅する照射光がコア部３２に入射される。このスペクトラムアナライザ４２から入射する光は、コア部３２を伝搬してグレーティング部３３でその一部の波長光のみが反射される。
光ファイバ３４は、その一端部においてスペクトラムアナライザ４２に接続されており、当該スペクトラムアナライザ４２が有する光源により、所定範囲の波長帯域を網羅する照射光がコア部３２に入射される。このスペクトラムアナライザ４２から入射する光は、コア部３２を伝搬してグレーティング部３３でその一部の波長光のみが反射される。
図２(b)は、コア部３２の光進行方向における屈折率変化を示す線図であり、図中の範
囲Ｌがグレーティング部３３における屈折率を示している。
かかる図に示すように、グレーティング部３３は、コア部３２の屈折率を一定の周期で変化するように構成されている。グレーティング部３３は、かかる屈折率の変化している境界部分で特定の波長の光のみ選択的に反射する。このグレーティング部３３に振動によりひずみ等の外乱が加えられると格子間隔の変化(伸縮)に伴って反射光の波長が変化する。
ここで、ＦＢＧ光ファイバセンサの反射光の波長変化ΔλBは、コアの実効屈折率をn、グレーティング間隔をΛ、ポッケルス係数をP11,P12、ポアソン比をν、印加歪をε、ファ
イバ材の温度係数をα、温度変化をΔTとすると次式（１）で表される(Alan D. Kersey, Fiber Grating Sensors”JOURNAL OF LIGHTWAVE TECHNOLOGY, Vol. 15, No. 8, 1997)。
囲Ｌがグレーティング部３３における屈折率を示している。
かかる図に示すように、グレーティング部３３は、コア部３２の屈折率を一定の周期で変化するように構成されている。グレーティング部３３は、かかる屈折率の変化している境界部分で特定の波長の光のみ選択的に反射する。このグレーティング部３３に振動によりひずみ等の外乱が加えられると格子間隔の変化(伸縮)に伴って反射光の波長が変化する。
ここで、ＦＢＧ光ファイバセンサの反射光の波長変化ΔλBは、コアの実効屈折率をn、グレーティング間隔をΛ、ポッケルス係数をP11,P12、ポアソン比をν、印加歪をε、ファ
イバ材の温度係数をα、温度変化をΔTとすると次式（１）で表される(Alan D. Kersey, Fiber Grating Sensors”JOURNAL OF LIGHTWAVE TECHNOLOGY, Vol. 15, No. 8, 1997)。
したがって、グレーティング部３３に振動を生じると、グレーティング部３３の歪み量εに変化を生じ、その結果、歪み量εに応じて反射光の波長が変動することとなる。振動がその振動源から良好に伝達されれば、グレーティング部３３は大きく歪みを生じ、波長の変化量ΔλBは大きく変動することとなるし、振動がその振動源から良好に伝達されな
い場合は、グレーティング部３３は小さく歪みを生じ、波長の変化量ΔλBは小さく変動
することとなる。
い場合は、グレーティング部３３は小さく歪みを生じ、波長の変化量ΔλBは小さく変動
することとなる。
ＭＦＣは繊維状の圧電素子に軸方向の垂直ひずみを発生させ、ＦＢＧは繊維状の光ファイバに発生した軸方向ひずみを検出している。これらの素子は共振周波数を持たずに広い周波数特性を有し、その挙動は強い指向性を有するため伝播経路が明確である。この二つの特徴によって、本計測システムでは、広帯域なラム波を指向性をもって伝播させることが可能となっている。また、ＦＢＧとＭＦＣは共に小型・軽量であり、柔軟で破断ひずみも高いため、積層板と一体化可能であり、大ひずみ下においても破断せずに高い信頼性がある、といった構造ヘルスモニタリングに適した特徴も有する。
図３(a)に光ファイバセンサとこれに接続したスペクトラムアナライザ４２の構成例を
示す。図３(a)に示すように、スペクトラムアナライザ４２は、光源６１と、光サーキュ
レータ６２と、ＡＷＧモジュール６３と、光電変換器６０とを備える。本構成例では、反射波長の異なる４つの光ファイバセンサ３０ａ〜ｄが直列に設けられた光ファイバ３４をスペクトラムアナライザ４２に接続する。最少構成としては、光ファイバセンサ３０は一つでよい。
示す。図３(a)に示すように、スペクトラムアナライザ４２は、光源６１と、光サーキュ
レータ６２と、ＡＷＧモジュール６３と、光電変換器６０とを備える。本構成例では、反射波長の異なる４つの光ファイバセンサ３０ａ〜ｄが直列に設けられた光ファイバ３４をスペクトラムアナライザ４２に接続する。最少構成としては、光ファイバセンサ３０は一つでよい。
光源６１は、光ファイバセンサ３０ａ〜ｄの反射波長の振動域をすべて含む広帯域の光源である。ラム波により光ファイバセンサ３０ａ〜ｄの反射波長が振動しても、常に完全な反射光が得られるようにするためである。
光サーキュレータ６２は、光源６１からの光を光ファイバセンサ３０ａ〜ｄ側へ進行させ、返ってきた光ファイバセンサ３０a〜ｄからの反射光を光ファイバ６９へ導出する。光ファイバ６９に導出された反射光はＡＷＧモジュール６３の入力ポートＰ０に導入される。
光サーキュレータ６２は、光源６１からの光を光ファイバセンサ３０ａ〜ｄ側へ進行させ、返ってきた光ファイバセンサ３０a〜ｄからの反射光を光ファイバ６９へ導出する。光ファイバ６９に導出された反射光はＡＷＧモジュール６３の入力ポートＰ０に導入される。
ＡＷＧモジュール６３は、ＡＷＧ基板６４を有する。ＡＷＧ基板６４には、光導波路技術によりガラス基板上にモノリシック集積された光波回路が形成されている。ＡＷＧ基板６４上の光波回路は、入出力スラブ導波路６５，６６とアレイ導波路６７と、出力導波路６８とを有し、入力ポートＰ０に並列接続された通過域の異なる８つの光学フィルタを構成している。ＡＷＧ基板６４上の光波回路は、波長多重された入力光を分配して８つの光学フィルタに通すことより波長分離し、８つの出力ポートＰ１〜Ｐ８にパラレル出力する。但し、実用時の出力ポートは８つに限定されない。
８つの出力ポートＰ１〜Ｐ８に対応する各光学フィルタの通過域を図３(b)のスペクト
ル図に示す。例えば図３(b)で、中心波長λ２の反射波長を有する光学ファイバセンサ３
０ｂの反射光入力分布７０が通過域７１に重なる部分に相当する反射光をある一つの光学フィルタが通過させ出力ポートＰ３へ出力すると共に、これに並行して通過域７２に重なる部分に相当する反射光を他の光学フィルタが通過させ出力ポートＰ４へ出力する。同様にして反射中心波長λ１の光学ファイバセンサ３０ａには出力ポートＰ１,Ｐ２が、反射中心波長λ３の光学ファイバセンサ３０ｃには出力ポートＰ５,Ｐ６が、反射中心波長λ４の光学ファイバセンサ３０ｄには出力ポートＰ７,Ｐ８がそれぞれ対応し、同様の原理で波長分離が可能である。上述したように最少構成としては、光ファイバセンサ３０は一つでよく、この場合、光学フィルタは２つで足りる。
ル図に示す。例えば図３(b)で、中心波長λ２の反射波長を有する光学ファイバセンサ３
０ｂの反射光入力分布７０が通過域７１に重なる部分に相当する反射光をある一つの光学フィルタが通過させ出力ポートＰ３へ出力すると共に、これに並行して通過域７２に重なる部分に相当する反射光を他の光学フィルタが通過させ出力ポートＰ４へ出力する。同様にして反射中心波長λ１の光学ファイバセンサ３０ａには出力ポートＰ１,Ｐ２が、反射中心波長λ３の光学ファイバセンサ３０ｃには出力ポートＰ５,Ｐ６が、反射中心波長λ４の光学ファイバセンサ３０ｄには出力ポートＰ７,Ｐ８がそれぞれ対応し、同様の原理で波長分離が可能である。上述したように最少構成としては、光ファイバセンサ３０は一つでよく、この場合、光学フィルタは２つで足りる。
代表して、１つの光ファイバセンサ３０からの反射光に対する処理内容を図４を参照して説明する。
図４(b)に示すように、光ファイバセンサ３０からの反射光の入力分布７３Ｔが現れる
。ＭＦＣアクチュエータ２１による加振時には、ＭＦＣアクチュエータ２１を発振源とするラム波が構造用複合材料Ｚを伝搬し、光ファイバセンサ３０は構造用複合材料Ｚから伝達されるラム波に応じて出力する反射光の波長を振動させる。この波長の振動を図示すると図４(a)の入力波７３Ｗとなる。
この波長の振動により、図４(b)に示す反射光入力分布７３Ｔは、上位、下位に交互に
シフトして微小に振動し、波長の値は増減を繰り返す。
このような波長振動において、図中７３Ｃは、反射光入力分布７３Ｔの中心波長の振動中心である。一方の光学フィルタの通過域７５Ｔの中心波長７５Ｃは振動中心７３Ｃの上域に固定されている。また、他方の光学フィルタの通過域７４Ｔの中心波長７４Ｃは振動中心７３Ｃの下域に固定されている。
また、中心波長７５Ｃ及び中心波長７４Ｃは、振動中心７３Ｃから反射光の波長振動の振幅以上に離れた位置に固定されている。
さらに、反射光入力分布７３Ｔの静止時において、上位の通過域７５Ｔの下位側のスロープ７５Ｔ−１は、反射光入力分布７３Ｔの上位側のスロープ７３Ｔ−１に交わり、上位の通過域７５Ｔと反射光入力分布７３Ｔとは波長振動の振幅以上の幅で重なる。
同様に、反射光入力分布７３Ｔの静止時において、下位の通過域７４Ｔの上位側のスロープ７４Ｔ−１は、反射光入力分布７３Ｔの下位側のスロープ７３Ｔ−２に交わり、下位の通過域７４Ｔと反射光入力分布７３Ｔとは波長振動の振幅以上の幅で重なる。
反射光入力分布７３Ｔに対し以上の位置関係に通過域７５Ｔ及び通過域７４Ｔを固定することにより、反射光の波長振動を高感度に検出することができる。
上位の光学フィルタは反射光入力分布７３Ｔが通過域７５Ｔに重なる部分に相当する反射光を通過させて出力する。同様に、下位の光学フィルタは反射光入力分布７３Ｔが通過域７４Ｔに重なる部分に相当する反射光を通過させて出力する。
図４(b)に示すように、光ファイバセンサ３０からの反射光の入力分布７３Ｔが現れる
。ＭＦＣアクチュエータ２１による加振時には、ＭＦＣアクチュエータ２１を発振源とするラム波が構造用複合材料Ｚを伝搬し、光ファイバセンサ３０は構造用複合材料Ｚから伝達されるラム波に応じて出力する反射光の波長を振動させる。この波長の振動を図示すると図４(a)の入力波７３Ｗとなる。
この波長の振動により、図４(b)に示す反射光入力分布７３Ｔは、上位、下位に交互に
シフトして微小に振動し、波長の値は増減を繰り返す。
このような波長振動において、図中７３Ｃは、反射光入力分布７３Ｔの中心波長の振動中心である。一方の光学フィルタの通過域７５Ｔの中心波長７５Ｃは振動中心７３Ｃの上域に固定されている。また、他方の光学フィルタの通過域７４Ｔの中心波長７４Ｃは振動中心７３Ｃの下域に固定されている。
また、中心波長７５Ｃ及び中心波長７４Ｃは、振動中心７３Ｃから反射光の波長振動の振幅以上に離れた位置に固定されている。
さらに、反射光入力分布７３Ｔの静止時において、上位の通過域７５Ｔの下位側のスロープ７５Ｔ−１は、反射光入力分布７３Ｔの上位側のスロープ７３Ｔ−１に交わり、上位の通過域７５Ｔと反射光入力分布７３Ｔとは波長振動の振幅以上の幅で重なる。
同様に、反射光入力分布７３Ｔの静止時において、下位の通過域７４Ｔの上位側のスロープ７４Ｔ−１は、反射光入力分布７３Ｔの下位側のスロープ７３Ｔ−２に交わり、下位の通過域７４Ｔと反射光入力分布７３Ｔとは波長振動の振幅以上の幅で重なる。
反射光入力分布７３Ｔに対し以上の位置関係に通過域７５Ｔ及び通過域７４Ｔを固定することにより、反射光の波長振動を高感度に検出することができる。
上位の光学フィルタは反射光入力分布７３Ｔが通過域７５Ｔに重なる部分に相当する反射光を通過させて出力する。同様に、下位の光学フィルタは反射光入力分布７３Ｔが通過域７４Ｔに重なる部分に相当する反射光を通過させて出力する。
したがって、反射光の波長の値が増加して反射光入力分布７３Ｔが上位にシフトすると、通過域７５Ｔを有する上位の光学フィルタの出力値は増加し、通過域７４Ｔを有する下位の光学フィルタの出力値は減少する。逆に、反射光の波長の値が減少して反射光入力分布７３Ｔが下位にシフトすると、通過域７５Ｔを有する上位の光学フィルタの出力値は減少し、通過域７４Ｔを有する下位の光学フィルタの出力値は増加する。
そのため、反射光の中心波長の変化が図４(a)に示す入力波７３Ｗにより振動する時、
通過域７５Ｔを有する上位の光学フィルタの出力値は、図４(c)に示す出力波７５Ｗを生
成し、通過域７４Ｔを有する下位の光学フィルタの出力値は、図４(c)に示す出力波７４
Ｗを生成する。図４(c)に示すように、出力波７４Ｗと出力波７５Ｗは逆位相の波動とな
る。
そのため、反射光の中心波長の変化が図４(a)に示す入力波７３Ｗにより振動する時、
通過域７５Ｔを有する上位の光学フィルタの出力値は、図４(c)に示す出力波７５Ｗを生
成し、通過域７４Ｔを有する下位の光学フィルタの出力値は、図４(c)に示す出力波７４
Ｗを生成する。図４(c)に示すように、出力波７４Ｗと出力波７５Ｗは逆位相の波動とな
る。
以上の原理により、図３に示すスペクトラムアナライザ４２は、加振時に８つの出力ポートＰ１〜Ｐ８それぞれに光波を出力し、これらを光電変換器６０により電気信号に変化して外部出力する。スペクトラムアナライザ４２の出力は、図示しないインターフェースを介してＡ／Ｄ変換されて演算処理装置５０に入力される。
演算処理装置５０は、図５に示すように、プログラムに従い演算処理を行うＣＰＵ５１と、各種の処理及び制御を行うためのプログラムを記憶するＲＯＭ５２と、各種の処理において一時的にデータ等を格納する作業領域となるＲＡＭ５３と、制御装置４１とデータの送受を図るインターフェース５４と、スペクトラムアナライザ４２からのデータを入力するインターフェース５５と、演算結果の表示データを表示装置５６に適正なフォーマットの画像信号に変換して表示装置５６に出力する画像出力インターフェース５７と、上記各構成間での各種指令又はデータの伝送を行うデータバス５８とを備えている。
損傷探知システム１０は、損傷探知の対象となる構造用複合材料Ｚに設置したＭＦＣアクチュエータ２１により構造用複合材料Ｚに対して振動を加えると共に、光ファイバセンサ３０により検出される振動波の伝搬状態から、光ファイバセンサ３０の近傍に損傷が発生しているか否かを探知する。そのために、演算処理装置５０は、ＲＯＭ５２に記憶された各種のプログラムをＲＡＭ５３を用いてＣＰＵ５１が処理することで以下に説明するような各種の機能を実行する。
ＣＰＵ５１は、ＲＯＭ５２に格納されたプログラムに従い、ＭＦＣアクチュエータ２１について駆動用電圧を印加するように制御装置４１の動作制御を行う。ＭＦＣアクチュエータ２１が複数ある場合は、選択するＭＦＣアクチュエータ２１としては、いずれのものでも良いが、例えば、振動の発生源とした場合に、光ファイバセンサ３０のグレーティング部３３までの間に構造用複合材料Ｚの損傷を生じやすい部分が存在するようなＭＦＣアクチュエータを選択することが望ましい。
ＣＰＵ５１は、ＲＯＭ５２に格納されたプログラムに従い、駆動用電圧を印加し、ＭＦＣアクチュエータ２１による加振中の一定期間にスペクトラムアナライザ４２から８並列で出力される出力波データを取得しＲＡＭ５３に記憶する処理を行う。
ＣＰＵ５１は、制御指令を発してＭＦＣアクチュエータ２１によりラム波の超音波振動を構造用複合材料Ｚに加振し、その加振動時に得られる光学フィルタの出力波７４Ｗと出力波７５Ｗの差動信号ｆ（ｔ）を数値化して得る。例えば、図７（ａ）に示す差動信号ｆ（ｔ）が得られる。
そしてＣＰＵ５１は、次式（２）によりｆ（ｔ）データにウェーブレット変換を施す。これによりｆ（ｔ）データは、周波数及び伝播時間の２次元に展開された伝播強度分布データに変換される。本データが、ラム波の光ファイバセンサ３０への伝播強度分布に相当し、図示すると例えば図７（ｃ）のごとくとなる。
そしてＣＰＵ５１は、次式（２）によりｆ（ｔ）データにウェーブレット変換を施す。これによりｆ（ｔ）データは、周波数及び伝播時間の２次元に展開された伝播強度分布データに変換される。本データが、ラム波の光ファイバセンサ３０への伝播強度分布に相当し、図示すると例えば図７（ｃ）のごとくとなる。
〔損傷探知動作〕
以上説明した基本構成を用い、さらに図１９又は図２１に示すようにＭＦＣアクチュエータ２１，２１及び光ファイバセンサ３０，３０をそれぞれ構造用複合材料Ｚの表裏同位置に設置して、以下の損傷探知動作を実行する。
以上説明した基本構成を用い、さらに図１９又は図２１に示すようにＭＦＣアクチュエータ２１，２１及び光ファイバセンサ３０，３０をそれぞれ構造用複合材料Ｚの表裏同位置に設置して、以下の損傷探知動作を実行する。
ＣＰＵ５１は、表裏のＭＦＣアクチュエータ２１，２１に同位相の波を発振させることにより、対称モードのみを構造用複合材料Ｚに加振し、上述したようにｆ（ｔ）データにウェーブレット変換を施して、例えば図１１(a)に示すごとくＳモードのみの周波数及び伝播時間の２次元展開データを得る。その後、ＣＰＵ５１は、図８に示す理論分散曲線に基づき、Ｓ０モード、Ｓ１モード、Ｓ２モード等のモードを特定して、特定モードの周波数ごとのウェーブレット係数最大値が生じる伝播時間を算出する。例えば、特定モードをＳ０モード及びＳ１モードとする場合は、図２３に示すごとくウェーブレット係数最大値が生じる周波数，伝播時間の関係が特定される。これが２次元展開データから抽出した一つの特徴値であって一つの計測結果である。
ＣＰＵ５１は、これを図２３に示すごとく表示装置５６に表示する。ＣＰＵ５１は、損傷状態が未知の被検体に対する計測結果と同一で損傷状態が既知の構造体に対する計測結果とともに表示する。これを検査者が参照し、それらの比較により損傷の有無や規模を推定することができる。
又は、Ｓ０モード及びＳ１モードの伝播時間は、層間剥離長さの増加により増加するので、ＣＰＵ５１は、例えば図２６に示すように損傷の無い場合に対する伝播時間の増加量の形態で表示する。これを検査者が参照し、損傷の有無や規模を推定することができる。
またＣＰＵ５１がさらに進んで、ＲＯＭ５２等に多数蓄積した損傷状態が既知の構造体に対する計測結果と、損傷状態が未知の被検体に対する計測結果とに基づき、当該被検体に対する損傷規模の推定演算を行い、その結果を表示装置５６に表示する形態ととってもよい。
なお、対称モード（Ｓモード）のデータを取得するために、上記方法に代えて、表裏の光ファイバセンサ３０，３０の出力値同士を加算することにより非対称モードを相殺して対称モードを強調した２次元展開データ（伝播強度分布データ）を得てもよい。
ＣＰＵ５１は、これを図２３に示すごとく表示装置５６に表示する。ＣＰＵ５１は、損傷状態が未知の被検体に対する計測結果と同一で損傷状態が既知の構造体に対する計測結果とともに表示する。これを検査者が参照し、それらの比較により損傷の有無や規模を推定することができる。
又は、Ｓ０モード及びＳ１モードの伝播時間は、層間剥離長さの増加により増加するので、ＣＰＵ５１は、例えば図２６に示すように損傷の無い場合に対する伝播時間の増加量の形態で表示する。これを検査者が参照し、損傷の有無や規模を推定することができる。
またＣＰＵ５１がさらに進んで、ＲＯＭ５２等に多数蓄積した損傷状態が既知の構造体に対する計測結果と、損傷状態が未知の被検体に対する計測結果とに基づき、当該被検体に対する損傷規模の推定演算を行い、その結果を表示装置５６に表示する形態ととってもよい。
なお、対称モード（Ｓモード）のデータを取得するために、上記方法に代えて、表裏の光ファイバセンサ３０，３０の出力値同士を加算することにより非対称モードを相殺して対称モードを強調した２次元展開データ（伝播強度分布データ）を得てもよい。
またＣＰＵ５１は、表裏のＭＦＣアクチュエータ２１，２１に逆位相の波を発振させることにより、非対称モードのみを構造用複合材料Ｚに加振し、上述したようにｆ（ｔ）データにウェーブレット変換を施して、例えば図１１(ｂ)に示すごとくＡモードのみの周波数及び伝播時間の２次元展開データを得る。その後、ＣＰＵ５１は、図８に示す理論分散曲線に基づき、Ａ０モード、Ａ１モード等のモードを特定して、特定モードの周波数ごとのウェーブレット係数最大値が生じる伝播時間を算出する。例えば、特定モードをＡ１モードとする場合は、図２２に示すごとくウェーブレット係数最大値が生じる周波数，伝播時間の関係が特定される。これが２次元展開データから抽出した一つの特徴値であって一つの計測結果である。
ＣＰＵ５１は、これを図２２に示すごとく表示装置５６に表示する。ＣＰＵ５１は、損傷状態が未知の被検体に対する計測結果と同一で損傷状態が既知の構造体に対する計測結果とともに表示する。これを検査者が参照し、それらの比較により損傷の有無や規模を推定することができる。
又は、Ａ１モードの伝播時間は損傷部でのＡ１モードより伝播速度の速いＳ０モードへの変換により減少するので、ＣＰＵ５１は、例えば図２５に示すように損傷の無い場合に対する伝播時間の減少量の形態で表示する。これを検査者が参照し、損傷の有無や規模を推定することができる。
またＣＰＵ５１は、Ａ１モードに関し、周波数対伝播時間の変化率を算出する。この変化率は、例えば、図２２の２５０〜４００（ｋＨｚ）の範囲の各試験片の計測データ群に近似直線を計算し、この近似直線の傾きが相当する。これも２次元展開データから抽出した一つの特徴値であって一つの計測結果である。ＣＰＵ５１は、この変化率（傾き）を例えば図２４に示すように数値やグラフで表示する。ここでも、ＣＰＵ５１は、損傷状態が未知の被検体に対する計測結果と同一で損傷状態が既知の構造体に対する計測結果とともに表示する。これを検査者が参照し、それらの比較により損傷の有無や規模を推定することができる。
またＣＰＵ５１がさらに進んで、ＲＯＭ５２等に多数蓄積した損傷状態が既知の構造体に対する計測結果と、損傷状態が未知の被検体に対する計測結果とに基づき、当該被検体に対する損傷規模の推定演算を行い、その結果を表示装置５６に表示する形態をとってもよい。推定演算の基礎データは、以上のＳ０モード及びＳ１モードの伝播時間の増加量、Ａ１モードの伝播時間の減少量、及び変化率（傾き）のうちいずれか１又は２以上とされる。
なお、非対称モード（Ａモード）のデータを取得するために、上記方法に代えて、表裏の光ファイバセンサ３０，３０の出力値同士を減算することにより対称モードを相殺して非対称モードを強調した２次元展開データ（伝播強度分布データ）を得てもよい。
ＣＰＵ５１は、これを図２２に示すごとく表示装置５６に表示する。ＣＰＵ５１は、損傷状態が未知の被検体に対する計測結果と同一で損傷状態が既知の構造体に対する計測結果とともに表示する。これを検査者が参照し、それらの比較により損傷の有無や規模を推定することができる。
又は、Ａ１モードの伝播時間は損傷部でのＡ１モードより伝播速度の速いＳ０モードへの変換により減少するので、ＣＰＵ５１は、例えば図２５に示すように損傷の無い場合に対する伝播時間の減少量の形態で表示する。これを検査者が参照し、損傷の有無や規模を推定することができる。
またＣＰＵ５１は、Ａ１モードに関し、周波数対伝播時間の変化率を算出する。この変化率は、例えば、図２２の２５０〜４００（ｋＨｚ）の範囲の各試験片の計測データ群に近似直線を計算し、この近似直線の傾きが相当する。これも２次元展開データから抽出した一つの特徴値であって一つの計測結果である。ＣＰＵ５１は、この変化率（傾き）を例えば図２４に示すように数値やグラフで表示する。ここでも、ＣＰＵ５１は、損傷状態が未知の被検体に対する計測結果と同一で損傷状態が既知の構造体に対する計測結果とともに表示する。これを検査者が参照し、それらの比較により損傷の有無や規模を推定することができる。
またＣＰＵ５１がさらに進んで、ＲＯＭ５２等に多数蓄積した損傷状態が既知の構造体に対する計測結果と、損傷状態が未知の被検体に対する計測結果とに基づき、当該被検体に対する損傷規模の推定演算を行い、その結果を表示装置５６に表示する形態をとってもよい。推定演算の基礎データは、以上のＳ０モード及びＳ１モードの伝播時間の増加量、Ａ１モードの伝播時間の減少量、及び変化率（傾き）のうちいずれか１又は２以上とされる。
なお、非対称モード（Ａモード）のデータを取得するために、上記方法に代えて、表裏の光ファイバセンサ３０，３０の出力値同士を減算することにより対称モードを相殺して非対称モードを強調した２次元展開データ（伝播強度分布データ）を得てもよい。
以上の実施形態にあっては、２つの光学フィルタの出力値の差動信号をウェーブレット変換の基礎データとしたが、本発明はこれに限定されず、１つの光学フィルタの出力値をウェーブレット変換の基礎データとしてもよい。
また以上の実施形態にあっては、特定モードのウェーブレット係数の最大ピーク値を算出したが、取得ラム波の特定モードの比較に適したパラメータであればいずれのパラメータの値であってもよい。
また以上の実施形態にあっては、光ファイバセンサの検出値を周波数及び伝播時間の２次元に展開する変換としてウェーブレット変換を適用したが、本発明はこれに限定されず、短時間フーリエ変換等他の変換手法を適用してもよい。
〔検証実験・解析〕
次に、本発明の理論説明及び本発明の実施に当たっての参考とするために、以下の実験・解析による検証内容を開示する。
次に、本発明の理論説明及び本発明の実施に当たっての参考とするために、以下の実験・解析による検証内容を開示する。
〔１．モード同定法（モード分離法）〕
まず始めに、擬似等方性CFRP積層板（T700S/2500、東レ（株）、[45/0/-45/90]3s、厚さ3.4mm）を対象として計測を行った。ＭＦＣ（M-2814-P2, Smart material Co., Ltd.）は長さ6mm、幅14mm、厚さ0.3mm、ＦＢＧセンサ（フジクラ（株））はセンサ長1.5mm、ポリイミド被覆付きで直径150μmである。両者を100mm離してCFRP積層板表面に貼り付け、計測を行った。貼りつけにはシアノアクリレート系接着剤であるアロンアルファ（コニシ（株））を用いた。ＭＦＣへの入力信号は、図６に示すようなfc=400kHzのサイン波1周期分にハミング窓を掛けた広帯域な信号を用いた。そしてＭＦＣで発振され、積層板を伝播してＦＢＧセンサで受振されたラム波の受振波形は、ノイズを除去するために32768回の計測波形の平均化を行った。その後、得られた受振波形に対して信号解析を行い、受振波に含まれるモード分散性を時間−周波数領域で表わした。信号解析には、窓関数として複素モルレ関数を用い、1次元複素連続ウェーブレット解析を行った。ＦＢＧセンサで受振した受振波形とそのフーリエスペクトル、およびウェーブレット変換結果を図７に示す。この結果より、特定の周波数で大きなピークが現れずに広帯域に渡る成分を受振していることが確認できる。また、ウェーブレット変換結果より、速度や周波数の異なる複数のモードが観察され、モード分散性があることが分かる。次に、この受振波形に現れる各モードを同定するために理論分散曲線を導出する。
まず始めに、擬似等方性CFRP積層板（T700S/2500、東レ（株）、[45/0/-45/90]3s、厚さ3.4mm）を対象として計測を行った。ＭＦＣ（M-2814-P2, Smart material Co., Ltd.）は長さ6mm、幅14mm、厚さ0.3mm、ＦＢＧセンサ（フジクラ（株））はセンサ長1.5mm、ポリイミド被覆付きで直径150μmである。両者を100mm離してCFRP積層板表面に貼り付け、計測を行った。貼りつけにはシアノアクリレート系接着剤であるアロンアルファ（コニシ（株））を用いた。ＭＦＣへの入力信号は、図６に示すようなfc=400kHzのサイン波1周期分にハミング窓を掛けた広帯域な信号を用いた。そしてＭＦＣで発振され、積層板を伝播してＦＢＧセンサで受振されたラム波の受振波形は、ノイズを除去するために32768回の計測波形の平均化を行った。その後、得られた受振波形に対して信号解析を行い、受振波に含まれるモード分散性を時間−周波数領域で表わした。信号解析には、窓関数として複素モルレ関数を用い、1次元複素連続ウェーブレット解析を行った。ＦＢＧセンサで受振した受振波形とそのフーリエスペクトル、およびウェーブレット変換結果を図７に示す。この結果より、特定の周波数で大きなピークが現れずに広帯域に渡る成分を受振していることが確認できる。また、ウェーブレット変換結果より、速度や周波数の異なる複数のモードが観察され、モード分散性があることが分かる。次に、この受振波形に現れる各モードを同定するために理論分散曲線を導出する。
上記実験と同様3.4mmのCFRP積層板における、伝播距離100mm地点での到達時間の理論分散曲線を導出したもの図８である。また、この分散曲線上で高次のモードの到達時間が急激に遅くなり、無限大に発散する周波数をカットオフ周波数と言う。この理論分散曲線と先ほどの受振波形のウェーブレット変換結果を比較すると、モード分散性は両者でよく一致していることが分かる。しかし、複数のモードが重なり合う300kHz以上の周波数領域においてはモードの同定を行うことは困難であり、正確なモード同定を行うためには、
これらの重なり合うモードを分離する必要がある。
そこで、これらのモードの分離手法として、ＭＦＣとＦＢＧセンサを積層板上下面の同じ位置にそれぞれ貼り付ける方法を用いた。図９に示すように、ＭＦＣを上下表面の同じ位置に貼り付け、お互いに同位相の波を発振すると、対称モードのみを発振することが出来る。反対に逆位相の波を発振すれば、非対称モードのみを発振することが出来る。
また、図１０に示すように、ＦＢＧセンサを上下表面の同じ位置に貼り付け、板の上下面での受振波形の和を取れば対称モードを、差を取れば非対称モードをそれぞれ分離することが出来る。
この二つの手法を用いて、受振ラム波に含まれるＳ（対称）モードとA（非対称）モードを分離し、それぞれウェーブレット変換を行った後に、先ほどの理論分散曲線と比較しモード同定を行った結果を図１１に示す。結果、モード分離をおこなうことによって複数モードの重なりが解消され、正確にモード同定をおこなうことが出来た。そして、受振波形にはＡ０,Ｓ０,Ａ１,Ｓ1,Ｓ２モードが存在することが確認された。
以上の結果より、上記のモード分離手法を用いれば、受振ラム波に含まれる各モードの同定が可能である。
これらの重なり合うモードを分離する必要がある。
そこで、これらのモードの分離手法として、ＭＦＣとＦＢＧセンサを積層板上下面の同じ位置にそれぞれ貼り付ける方法を用いた。図９に示すように、ＭＦＣを上下表面の同じ位置に貼り付け、お互いに同位相の波を発振すると、対称モードのみを発振することが出来る。反対に逆位相の波を発振すれば、非対称モードのみを発振することが出来る。
また、図１０に示すように、ＦＢＧセンサを上下表面の同じ位置に貼り付け、板の上下面での受振波形の和を取れば対称モードを、差を取れば非対称モードをそれぞれ分離することが出来る。
この二つの手法を用いて、受振ラム波に含まれるＳ（対称）モードとA（非対称）モードを分離し、それぞれウェーブレット変換を行った後に、先ほどの理論分散曲線と比較しモード同定を行った結果を図１１に示す。結果、モード分離をおこなうことによって複数モードの重なりが解消され、正確にモード同定をおこなうことが出来た。そして、受振波形にはＡ０,Ｓ０,Ａ１,Ｓ1,Ｓ２モードが存在することが確認された。
以上の結果より、上記のモード分離手法を用いれば、受振ラム波に含まれる各モードの同定が可能である。
〔２．特定モードの伝播時間の変化原因（層間剥離部におけるモード変換挙動）〕
前節において各モードの同定を行い、計測結果に含まれるモード分散性を理解することが出来るようになった。次に、このモード分散性の変化から生じるモード変換挙動を実験
と解析によって解明する。
前節において各モードの同定を行い、計測結果に含まれるモード分散性を理解することが出来るようになった。次に、このモード分散性の変化から生じるモード変換挙動を実験
と解析によって解明する。
（１）剥離部の板厚変化によるモード変換
ラム波の伝播速度は周波数と板厚の積に依存するため、板厚が変化するとラム波のモード分散性も変化する。よって、図１２に示すように、積層板内部に層間剥離が発生すると、剥離部における伝播経路の板厚は健全部よりも減少するため、健全部と剥離部ではモード分散性が異なる。このモード分散性の変化によって、健全部を伝播してきたラム波にはモード変換が生じ、健全部とは異なるモード形態で剥離部を伝播すると考えられる。
例えば、板厚3.4mmの積層板では周波数300kHzのラム波の伝播形態としてＡ０,Ｓ０,Ａ１モードの三つが存在するが、積層板の中央に層間剥離が発生し、剥離部での板厚が1.7mmに変化すると、Ａ０,Ｓ０モードの二つしか伝播形態が存在しない。
よって、健全部をＡ１モードとして伝播してきたラム波は剥離部でモード変換し、Ａ０,Ｓ０モードとして伝播していくこととなる。しかし、どちらのモードとして剥離部を伝播するかは理論分散曲線からは求められない。そこで、実際の層間剥離部で生じるモード変換挙動の解明を、実験と有限要素解析から行う。
ラム波の伝播速度は周波数と板厚の積に依存するため、板厚が変化するとラム波のモード分散性も変化する。よって、図１２に示すように、積層板内部に層間剥離が発生すると、剥離部における伝播経路の板厚は健全部よりも減少するため、健全部と剥離部ではモード分散性が異なる。このモード分散性の変化によって、健全部を伝播してきたラム波にはモード変換が生じ、健全部とは異なるモード形態で剥離部を伝播すると考えられる。
例えば、板厚3.4mmの積層板では周波数300kHzのラム波の伝播形態としてＡ０,Ｓ０,Ａ１モードの三つが存在するが、積層板の中央に層間剥離が発生し、剥離部での板厚が1.7mmに変化すると、Ａ０,Ｓ０モードの二つしか伝播形態が存在しない。
よって、健全部をＡ１モードとして伝播してきたラム波は剥離部でモード変換し、Ａ０,Ｓ０モードとして伝播していくこととなる。しかし、どちらのモードとして剥離部を伝播するかは理論分散曲線からは求められない。そこで、実際の層間剥離部で生じるモード変換挙動の解明を、実験と有限要素解析から行う。
（２）実験
実際の層間剥離の開始端部で生じるモード変換挙動の解明を行うために、擬似等方性CFRP積層板（T700S/2500、東レ（株）、[45/0/-45/90]3s、厚さ3.4mm）の板厚方向の中央に層間剥離が存在する場合を模擬した。なお、モード分離によって受振ラム波のモード同定を行うため、剥離部のモード分散性を計測するためには、模擬した層間剥離の内表面にＦＢＧセンサを貼り付ける必要がある。そこで、厚さ1.7mmのCFRP積層板を二枚用意し、一方にＦＢＧセンサを接着したのち、ＦＢＧセンサを接着した面が内側に来るように、板の片端から長さ60mm範囲をエポキシ系接着剤Araldite Standard（ハンツマン・アドバンスト・マテリアル（株））で接着した。二枚のCFRP積層板は、積層構成[45/0/-45/90]3sを模擬するために、積層構成[45/0/-45/90]3で作製し、接着面で対称になるように接着した。試験片の寸法を図１３に示す。また、板の幅は90mmである。
ＭＦＣ（M-2814-P2）は長さ6mm、幅14mm、厚さ0.3mmのものを用い、積層板の上下表面に一つずつ貼り付けた。また、ＭＦＣ先端から距離30mm,板厚3.4mmの地点（健全部）、および距離70mm,板厚1.7mmの地点（剥離部）の二点でＦＢＧセンサを積層板の上下表面に貼り付けラム波を受振する。実験に用いたＦＢＧセンサ（フジクラ（株））はセンサ長1.5mm、ポリイミド被覆付きで直径150μm。素子の貼り付けにはアロンアルファ（コニシ（株））を用いた。入力信号はfc=400kHzのサイン波1周期分にハミング窓を掛けた信号であり、受振波形はノイズを除去するため、32768回の計測波形の平均化を行う。上下面のＭＦＣを用いて、Ｓ（対称）モードのみの発振を行うことによって求めたＳモードのモード変換挙動を図１４に、Ａ（非対称）モードのみの発振を行うことによって求めたＡモードのモード変換挙動を図１５にそれぞれ示す。
実際の層間剥離の開始端部で生じるモード変換挙動の解明を行うために、擬似等方性CFRP積層板（T700S/2500、東レ（株）、[45/0/-45/90]3s、厚さ3.4mm）の板厚方向の中央に層間剥離が存在する場合を模擬した。なお、モード分離によって受振ラム波のモード同定を行うため、剥離部のモード分散性を計測するためには、模擬した層間剥離の内表面にＦＢＧセンサを貼り付ける必要がある。そこで、厚さ1.7mmのCFRP積層板を二枚用意し、一方にＦＢＧセンサを接着したのち、ＦＢＧセンサを接着した面が内側に来るように、板の片端から長さ60mm範囲をエポキシ系接着剤Araldite Standard（ハンツマン・アドバンスト・マテリアル（株））で接着した。二枚のCFRP積層板は、積層構成[45/0/-45/90]3sを模擬するために、積層構成[45/0/-45/90]3で作製し、接着面で対称になるように接着した。試験片の寸法を図１３に示す。また、板の幅は90mmである。
ＭＦＣ（M-2814-P2）は長さ6mm、幅14mm、厚さ0.3mmのものを用い、積層板の上下表面に一つずつ貼り付けた。また、ＭＦＣ先端から距離30mm,板厚3.4mmの地点（健全部）、および距離70mm,板厚1.7mmの地点（剥離部）の二点でＦＢＧセンサを積層板の上下表面に貼り付けラム波を受振する。実験に用いたＦＢＧセンサ（フジクラ（株））はセンサ長1.5mm、ポリイミド被覆付きで直径150μm。素子の貼り付けにはアロンアルファ（コニシ（株））を用いた。入力信号はfc=400kHzのサイン波1周期分にハミング窓を掛けた信号であり、受振波形はノイズを除去するため、32768回の計測波形の平均化を行う。上下面のＭＦＣを用いて、Ｓ（対称）モードのみの発振を行うことによって求めたＳモードのモード変換挙動を図１４に、Ａ（非対称）モードのみの発振を行うことによって求めたＡモードのモード変換挙動を図１５にそれぞれ示す。
図１４の結果より、Ｓモードのみの発振を行った場合、健全部ではＳ０モードとＳ１モードのみが観察される。図１２（ｃ）に示す板厚1.7mmの場合の理論分散曲線より、剥離部ではＳ１モードは800kHz以上にしか存在しないことが分かっているため、健全部で観察されたＳ１モードは剥離部においてモード変換し他のモードとして伝播していると考えられる。そこで、剥離部に存在するモードを観察すると、Ｓ０モードとＡ１モードの二つが観察される。この結果より、層間剥離の開始部では、「Ｓ１モード→Ｓ０モード,Ａ１モード」のモード変換が生じていることが確認された。
図１５の結果より、Ａモードのみの発振を行った場合、健全部ではＡ０モードとＡ１モードのみが観察される。図１２（ｃ）に示す板厚1.7mmの場合の理論分散曲線より、剥離部ではＡ１モードは500kHz以上にしか存在しないことが分かっているため、健全部で観察されたＡ１モードは剥離部においてモード変換し他のモードとして伝播していると考えられる。そこで、剥離部に存在するモードを観察すると、Ｓ０モードとＡ１モードの二つが観察される（Ａ０モードは到達時間が遅いためモード変換の対象としては考慮しない）。この剥離部で観察されたＡ１モードはモード変換ではなく、500kHz以上のＡ１モードがそのまま伝播したものと考えられる。よって、この結果より、層間剥離の開始部では、「Ａ１モード→Ｓ０モード」のモード変換が生じていることが確認された。
図１５の結果より、Ａモードのみの発振を行った場合、健全部ではＡ０モードとＡ１モードのみが観察される。図１２（ｃ）に示す板厚1.7mmの場合の理論分散曲線より、剥離部ではＡ１モードは500kHz以上にしか存在しないことが分かっているため、健全部で観察されたＡ１モードは剥離部においてモード変換し他のモードとして伝播していると考えられる。そこで、剥離部に存在するモードを観察すると、Ｓ０モードとＡ１モードの二つが観察される（Ａ０モードは到達時間が遅いためモード変換の対象としては考慮しない）。この剥離部で観察されたＡ１モードはモード変換ではなく、500kHz以上のＡ１モードがそのまま伝播したものと考えられる。よって、この結果より、層間剥離の開始部では、「Ａ１モード→Ｓ０モード」のモード変換が生じていることが確認された。
（３）有限要素解析による検証
上記（２）の実験から求めたモード変換挙動を検証するために、二次元の有限要素解析を行った。図１６に有限要素モデルとその寸法を示す。モデルの構築および解析にはLS-DYNA 971を使用した。解析モデルの要素は2Dshell要素（平面ひずみ）を用いた。Meshサイズは、波長の短い高周波も計算できるように十分小さく0.125mmとした。ＭＦＣとCFRP積層板の接触部はNode結合を行い、ＭＦＣへの入力波形は実験と同様、周波数400kHzのサイン波1周期分にハミング窓を掛けたものを用いた。また、LS-DYNAでは圧電効果を計算できないため、圧電効果をＭＦＣの伸縮方向の熱膨張率として付与して模擬した。以上のような条件のもと、図１６に示す三つの受振点（健全部：伝播距離20mm,板厚3.4mm、剥離部：伝播距離60mm,板厚1.7mm、健全部（剥離通過後）：伝播距離100mm,板厚3.4mm）において、x方向ひずみの時間履歴を計算し受振波形とした。上下面のＭＦＣを用いて、Ｓ（対称）モードのみの発振を行うことによって求めたＳモードのモード変換挙動を図１７に、Ａ（非対称）モードのみの発振を行うことによって求めたＡモードのモード変換挙動を図１８にそれぞれ示す。
上記（２）の実験から求めたモード変換挙動を検証するために、二次元の有限要素解析を行った。図１６に有限要素モデルとその寸法を示す。モデルの構築および解析にはLS-DYNA 971を使用した。解析モデルの要素は2Dshell要素（平面ひずみ）を用いた。Meshサイズは、波長の短い高周波も計算できるように十分小さく0.125mmとした。ＭＦＣとCFRP積層板の接触部はNode結合を行い、ＭＦＣへの入力波形は実験と同様、周波数400kHzのサイン波1周期分にハミング窓を掛けたものを用いた。また、LS-DYNAでは圧電効果を計算できないため、圧電効果をＭＦＣの伸縮方向の熱膨張率として付与して模擬した。以上のような条件のもと、図１６に示す三つの受振点（健全部：伝播距離20mm,板厚3.4mm、剥離部：伝播距離60mm,板厚1.7mm、健全部（剥離通過後）：伝播距離100mm,板厚3.4mm）において、x方向ひずみの時間履歴を計算し受振波形とした。上下面のＭＦＣを用いて、Ｓ（対称）モードのみの発振を行うことによって求めたＳモードのモード変換挙動を図１７に、Ａ（非対称）モードのみの発振を行うことによって求めたＡモードのモード変換挙動を図１８にそれぞれ示す。
図１７の結果より、層間剥離の開始部では、実験と同様「Ｓ１モード→Ｓ０モード,Ａ１モード」のモード変換が生じていることが確認された。また、層間剥離通過後の健全部では、剥離通過前の健全部と同様の分散性が観察され、Ｓ０モードとＳ１モードが観察された。このＳ１モードは剥離部を伝播して健全部に戻る際に、剥離部のＳ０モード,Ａ１モードが、健全部において再びＳ１モードにモード変換したものと考えられる。
次に、図１８の結果より、層間剥離の開始部では、実験と同様「Ａ１モード→Ｓ０モード」のモード変換が生じていることが確認された。また、層間剥離通過後の健全部では、
剥離通過前の健全部と同様の分散性が観察され、Ａ０モードとＡ１モードが観察された。
このＡ１モードは剥離部を伝播して健全部に戻る際に、剥離部のＳ０モードが健全部において再びＡ１モードにモード変換したものと考えられる。
以上より、実験によって求めたモード変換挙動の妥当性を示し、さらに層間剥離が終了する地点でのモード変換挙動に関しても解明した。結果、層間剥離部を通過する際に、以下の二つのモード変換挙動が存在することが確認された。
・「Ｓ１モード→Ｓ０モード,Ａ１モード→Ｓ１モード」
・ 「Ａ１モード→Ｓ０モード→Ａ１モード」
次に、図１８の結果より、層間剥離の開始部では、実験と同様「Ａ１モード→Ｓ０モード」のモード変換が生じていることが確認された。また、層間剥離通過後の健全部では、
剥離通過前の健全部と同様の分散性が観察され、Ａ０モードとＡ１モードが観察された。
このＡ１モードは剥離部を伝播して健全部に戻る際に、剥離部のＳ０モードが健全部において再びＡ１モードにモード変換したものと考えられる。
以上より、実験によって求めたモード変換挙動の妥当性を示し、さらに層間剥離が終了する地点でのモード変換挙動に関しても解明した。結果、層間剥離部を通過する際に、以下の二つのモード変換挙動が存在することが確認された。
・「Ｓ１モード→Ｓ０モード,Ａ１モード→Ｓ１モード」
・ 「Ａ１モード→Ｓ０モード→Ａ１モード」
このようなモード変換挙動により、剥離部を伝播する時のモードと、健全時に伝播する時のモードは異なる。例えば「Ａ１モード→Ｓ０モード→Ａ１モード」のモード変換では、図１９に示すように、健全時はすべての伝播経路をＡ１モードとして伝播するのに対し、剥離が発生するとその領域をＳ０モードとして伝播する。このとき、健全部（板厚3.4mm）のＡ１モードと剥離部のＳ０モード（板厚1.7mm）の伝播速度は図２０に示すように異なっており、Ｓ０モードの方がＡ１モードよりも速度が速い。よって、ＦＢＧセンサで受振されるＡ１モードの到達時間は、剥離が発生した場合の方が健全時よりも早くなる。この到達時間の変化は、健全部を伝播するモードと剥離部を伝播するモードの伝播速度の差と、剥離の長さに応じて生じる。よって、この差を指標とすることで、層間剥離の検知と、その剥離長さの定量評価が可能となる。
〔３．人工層間剥離検知による検証実験と解析〕
（１）検証実験
剥離部を通過した後の波に、実際に到達時間の変化が生じているかどうかを実験から検証し、本件発明の有効性を示す。
そこで、CFRP擬似等方性積層板の成形中に、板厚方向中央の隣り合う90°層の間に、厚さ50μmのテフロン（登録商標）フィルムを二枚重ねて埋め込むことで、板厚方向の中央に、剥離長さL=20,40,60mmの層間剥離を人工的に導入した３種類の積層板を作製した。そして、それらの剥離領域を通過するように広帯域なラム波を伝播させ、受振波形を計測した。その実験構成を図２１に示す。ＭＦＣとＦＢＧセンサはそれぞれ板の上下表面の同じ場所に接着し、前節でモード変換挙動を解明した場合と同様にモード分離を行う。
この実験構成を用いて、Ｌ=20,40,60mmの人工剥離を導入した積層板を用いて人工層間剥離の検知実験を行った。それらの結果と、健全な積層板（L=0mm）を対象にして計測を行った場合の結果を比較し、到達時間の変化について評価する。
そのために、受振波形をウェーブレット変換した後、各周波数のウェーブレット係数の最大値を抽出した。層間剥離が存在する場合の、このウェーブレット係数最大値の時刻の健全状態からの変化量が、到達時間の変化量に対応する。
上下面ＭＦＣを用いてＡモードを発振した場合に、ＦＢＧセンサで計測されるＡモードの「周波数ごとのウェーブレット係数最大値が現れる時間」を、Ｌ=0,20,40,60mmの場合について、到達時間に比較的大きな変化が生じている200〜700kHzのＡ１モードに関して、図２２に示す。
次に、上下面ＭＦＣを用いてＳモードを発振した場合に、ＦＢＧセンサで計測されるＳモードの「周波数ごとのウェーブレット係数最大値が現れる時間」を、Ｌ=0,20,40,60mmの場合について、到達時間に変化が生じている400〜600kHzのＳ０,Ｓ１モードに関して、図２３に示す。
図２２より、剥離の長さが長くなるにつれてＡ１モードの到達時間が早くなる様子が観察される。また、200〜500kHzの周波数領域のモード分散性の傾きが変化する様子も観察
される。これは、図２０に示されるように、Ａ１モードのカットオフ周波数に近くなる（周波数が低くなる）ほどＳ０モードとＡ１モードの伝播速度の差が大きくなるからである。
また、図２３からは、剥離の長さが長くなるにつれてＳ０,Ｓ１モードの到達時間が遅くなる様子が観察される。
以上の結果より、層間剥離が発生すると確かに到達時間に変化が生じており、本発明が有効であることが示された。
（１）検証実験
剥離部を通過した後の波に、実際に到達時間の変化が生じているかどうかを実験から検証し、本件発明の有効性を示す。
そこで、CFRP擬似等方性積層板の成形中に、板厚方向中央の隣り合う90°層の間に、厚さ50μmのテフロン（登録商標）フィルムを二枚重ねて埋め込むことで、板厚方向の中央に、剥離長さL=20,40,60mmの層間剥離を人工的に導入した３種類の積層板を作製した。そして、それらの剥離領域を通過するように広帯域なラム波を伝播させ、受振波形を計測した。その実験構成を図２１に示す。ＭＦＣとＦＢＧセンサはそれぞれ板の上下表面の同じ場所に接着し、前節でモード変換挙動を解明した場合と同様にモード分離を行う。
この実験構成を用いて、Ｌ=20,40,60mmの人工剥離を導入した積層板を用いて人工層間剥離の検知実験を行った。それらの結果と、健全な積層板（L=0mm）を対象にして計測を行った場合の結果を比較し、到達時間の変化について評価する。
そのために、受振波形をウェーブレット変換した後、各周波数のウェーブレット係数の最大値を抽出した。層間剥離が存在する場合の、このウェーブレット係数最大値の時刻の健全状態からの変化量が、到達時間の変化量に対応する。
上下面ＭＦＣを用いてＡモードを発振した場合に、ＦＢＧセンサで計測されるＡモードの「周波数ごとのウェーブレット係数最大値が現れる時間」を、Ｌ=0,20,40,60mmの場合について、到達時間に比較的大きな変化が生じている200〜700kHzのＡ１モードに関して、図２２に示す。
次に、上下面ＭＦＣを用いてＳモードを発振した場合に、ＦＢＧセンサで計測されるＳモードの「周波数ごとのウェーブレット係数最大値が現れる時間」を、Ｌ=0,20,40,60mmの場合について、到達時間に変化が生じている400〜600kHzのＳ０,Ｓ１モードに関して、図２３に示す。
図２２より、剥離の長さが長くなるにつれてＡ１モードの到達時間が早くなる様子が観察される。また、200〜500kHzの周波数領域のモード分散性の傾きが変化する様子も観察
される。これは、図２０に示されるように、Ａ１モードのカットオフ周波数に近くなる（周波数が低くなる）ほどＳ０モードとＡ１モードの伝播速度の差が大きくなるからである。
また、図２３からは、剥離の長さが長くなるにつれてＳ０,Ｓ１モードの到達時間が遅くなる様子が観察される。
以上の結果より、層間剥離が発生すると確かに到達時間に変化が生じており、本発明が有効であることが示された。
（２）有限要素法による検証
上記２．（３）で用いた二次元有限要素解析モデルの層間剥離長さをＬ=20, 40, 60mmで変化させ、上記実験と同様の実験構成で解析を行った。その後、実験と同様に、振幅最大値をＡ１モードとＳ０,Ｓ１モードに関して求め、到達時間の変化を観察すると、実験（図２２,図２３）の結果とよく一致し、同様の変化が観察された。
上記２．（３）で用いた二次元有限要素解析モデルの層間剥離長さをＬ=20, 40, 60mmで変化させ、上記実験と同様の実験構成で解析を行った。その後、実験と同様に、振幅最大値をＡ１モードとＳ０,Ｓ１モードに関して求め、到達時間の変化を観察すると、実験（図２２,図２３）の結果とよく一致し、同様の変化が観察された。
（３）剥離長さの定量評価
さらに、上記（１）の実験結果と上記（２）の解析結果から観察された到達時間の変化やモード分散性の傾きを指標として、剥離長さの定量評価の可能性を示す。
到達時間に変化が生じている周波数領域の振幅最大値のプロットから線形近似直線を算出し、その近似直線を用いて以下の指標を求めた。指標は実験結果と解析結果からそれぞれ求めた。各指標を剥離長さごとにプロットした結果を図２４,図２５,図２６に示す。
図２４は、250〜400kHzのＡ１モードのモード分散性の傾きを、図２５は300kHzにおけるＡ１モードの到達時間の減少量を、図２６は400kHz（解析は350kHz）におけるＳ０,Ｓ１モードの到達時間の増加量を示す。
図２４の結果より、剥離長さが大きくなるにつれ、モード分散性の傾きが小さくなっていく様子が観察される。また、図２５,図２６の結果より、剥離長さが大きくなるにつれ、Ａ１モードの到達時間の減少量と、Ｓ０,Ｓ１モードの到達時間の増加量が大きくなっていく様子が観察される。しかも、これらの指標は剥離長さにほぼ比例して変化している。よって、これらの指標を用いることで剥離長さの定量評価が可能である。
さらに、上記（１）の実験結果と上記（２）の解析結果から観察された到達時間の変化やモード分散性の傾きを指標として、剥離長さの定量評価の可能性を示す。
到達時間に変化が生じている周波数領域の振幅最大値のプロットから線形近似直線を算出し、その近似直線を用いて以下の指標を求めた。指標は実験結果と解析結果からそれぞれ求めた。各指標を剥離長さごとにプロットした結果を図２４,図２５,図２６に示す。
図２４は、250〜400kHzのＡ１モードのモード分散性の傾きを、図２５は300kHzにおけるＡ１モードの到達時間の減少量を、図２６は400kHz（解析は350kHz）におけるＳ０,Ｓ１モードの到達時間の増加量を示す。
図２４の結果より、剥離長さが大きくなるにつれ、モード分散性の傾きが小さくなっていく様子が観察される。また、図２５,図２６の結果より、剥離長さが大きくなるにつれ、Ａ１モードの到達時間の減少量と、Ｓ０,Ｓ１モードの到達時間の増加量が大きくなっていく様子が観察される。しかも、これらの指標は剥離長さにほぼ比例して変化している。よって、これらの指標を用いることで剥離長さの定量評価が可能である。
〔４．結論〕
以上のように、まず始めに、広帯域超音波送受振システムで計測される広帯域ラム波の各モードの同定を行った。そのための手法として、対称/非対称モードを分離する手法を提案し、この分離手法を用いればモード同定が可能なことを示した。
次に、層間剥離部におけるモード変換挙動を実験と解析から解明し、「Ｓ１モード→Ｓ０モード,Ａ１モード→Ｓ１モード」、「Ａ１モード→Ｓ０モード→Ａ１モード」の二つのモード変換挙動の存在を確認した。
その後、モード変換によるラム波の速度変化を利用した本発明に係る剥離検知手法の有効性を実験と解析から検証した。結果、剥離部における速度変化が到達時間の変化として観察されることを確認した。
最後に、Ａ１モードのモード分散性の傾きと、Ａ１モードの到達時間の減少量、Ｓ０,Ｓ１モードの到達時間の増加量の三つの指標のそれぞれにより、剥離長さの定量評価が可能であることを示した。
以上のように、まず始めに、広帯域超音波送受振システムで計測される広帯域ラム波の各モードの同定を行った。そのための手法として、対称/非対称モードを分離する手法を提案し、この分離手法を用いればモード同定が可能なことを示した。
次に、層間剥離部におけるモード変換挙動を実験と解析から解明し、「Ｓ１モード→Ｓ０モード,Ａ１モード→Ｓ１モード」、「Ａ１モード→Ｓ０モード→Ａ１モード」の二つのモード変換挙動の存在を確認した。
その後、モード変換によるラム波の速度変化を利用した本発明に係る剥離検知手法の有効性を実験と解析から検証した。結果、剥離部における速度変化が到達時間の変化として観察されることを確認した。
最後に、Ａ１モードのモード分散性の傾きと、Ａ１モードの到達時間の減少量、Ｓ０,Ｓ１モードの到達時間の増加量の三つの指標のそれぞれにより、剥離長さの定量評価が可能であることを示した。
１０ 損傷探知システム
２１ ＭＦＣアクチュエータ
３０ 光ファイバセンサ
３２ コア部
３３ グレーティング部
４２ スペクトラムアナライザ
５０ 演算処理装置
６３ ＡＷＧモジュール
６４ ＡＷＧ基板
６５，６６ 入出力スラブ導波路
６７ アレイ導波路
６８ 出力導波路
Ｐ０ 入力ポート
Ｐ１〜Ｐ８ 出力ポート
７０，７３Ｔ 反射光入力分布
７１，７２，７４Ｔ，７５Ｔ 通過域
７３Ｃ 振動中心
１０１ 剥離
１０２ ラム波
１０３ モード変換挙動
１０４ リフレクション
１０５ 接着層
２１ ＭＦＣアクチュエータ
３０ 光ファイバセンサ
３２ コア部
３３ グレーティング部
４２ スペクトラムアナライザ
５０ 演算処理装置
６３ ＡＷＧモジュール
６４ ＡＷＧ基板
６５，６６ 入出力スラブ導波路
６７ アレイ導波路
６８ 出力導波路
Ｐ０ 入力ポート
Ｐ１〜Ｐ８ 出力ポート
７０，７３Ｔ 反射光入力分布
７１，７２，７４Ｔ，７５Ｔ 通過域
７３Ｃ 振動中心
１０１ 剥離
１０２ ラム波
１０３ モード変換挙動
１０４ リフレクション
１０５ 接着層
Claims (8)
- 被検体に広帯域ラム波の超音波振動を加える加振装置と、
前記被検体から伝達される広帯域ラム波を検知する振動検知センサと、
前記加振装置の発振動作を制御し、前記振動検知センサの出力値を処理して計測結果を出力する処理装置とを備え、
前記処理装置は、
前記加振装置による加振時に得られる前記振動検知センサの出力値の出力値を変換して周波数及び伝播時間の２次元に展開した伝播強度分布データを得て、ラム波の非対称モードを同定し、同定したラム波の非対称モードの中から選ばれる１つ又は２つ以上の特定の非対称モードについて、周波数ごとに伝播強度の最大値が生じる伝播時間を算出し、当該最大値が生じる周波数及び伝播時間の関係を特定し、
特定した前記周波数及び伝播時間の関係と、記憶手段に蓄積されている、損傷状態が既知の被検体において予め特定した前記周波数及び伝播時間の関係と対比して、当該被検体に対する損傷規模の推定演算を行って表示することを特徴とする損傷診断システム。 - 前記加振装置および前記振動検知センサは、前記被検体の表裏の同位置にそれぞれ設置され、
前記処理装置は、前記表裏の前記加振装置に逆位相の波を発振させることにより非対称に発振させ、表裏の前記振動検知センサの出力値同士を加算して算出したデータを変換して周波数及び伝播時間の２次元に展開した前記伝播強度分布データを得ることを特徴とする請求項１に記載の損傷診断システム。 - 前記処理装置は、
前記伝播強度分布データの最大値が生じる周波数及び伝播時間の関係を特定する代わりに、
前記周波数ごとに伝播強度の最大値が生じる伝播時間を算出し特定した周波数及び伝播時間の関係を示すデータを、横軸が周波数、縦軸が伝播時間で表示されるグラフ上にプロットし、ある周波数範囲における当該データの傾きを変化率として算出し、
算出した前記変化率と、記憶手段に蓄積されている、損傷状態が既知の被検体において予め算出した前記変化率と対比して、当該被検体に対する損傷規模の推定演算を行って表示することを特徴とする請求項１又は請求項２に記載の損傷診断システム。 - 前記処理装置は、
前記変化率を算出する代わりに、
前記伝播強度分布データの最大値における伝播時間の、前記被検体に損傷がない場合に対する減少量を算出し、
算出した前記減少量に基づいて、当該被検体に対する損傷規模の推定演算を行って表示することを特徴とする請求項３に記載の損傷診断システム。 - 被検体に広帯域ラム波の超音波振動を加える加振装置と、
前記被検体から伝達される広帯域ラム波を検知する振動検知センサと、
前記加振装置の発振動作を制御し、前記振動検知センサの出力値を処理して計測結果を出力する処理装置とを備え、
前記処理装置は、
前記加振装置による加振時に得られる前記振動検知センサの出力値の出力値を変換して周波数及び伝播時間の２次元に展開した伝播強度分布データを得て、ラム波の対称モードを同定し、同定したラム波の対称モードの中から選ばれる１つ又は２つ以上の特定の対称モードについて、周波数ごとに伝播強度の最大値が生じる伝播時間を算出し、当該最大値が生じる周波数及び伝播時間の関係を特定し、
特定した前記周波数及び伝播時間の関係と、記憶手段に蓄積されている、損傷状態が既知の被検体において予め特定した前記周波数及び伝播時間の関係と対比して、当該被検体に対する損傷規模の推定演算を行って表示することを特徴とする損傷診断システム。 - 前記加振装置および前記振動検知センサは、前記被検体の表裏の同位置にそれぞれ設置され、
前記処理装置は、前記表裏の前記加振装置に同位相の波を発振させることにより対称に発振させ、表裏の前記振動検知センサの出力値同士を加算して算出したデータを変換して周波数及び伝播時間の２次元に展開した伝播強度分布データを得ることを特徴とする請求項５に記載の損傷診断システム。 - 前記処理装置は、
前記伝播強度分布データの最大値が生じる周波数及び伝播時間の関係を特定する代わりに、
前記伝播強度分布データの最大値における伝播時間の、前記被検体に損傷がない場合に対する増加量を算出し、
算出した前記増加量に基づいて、当該被検体に対する損傷規模の推定演算を行って表示することを特徴とする請求項５又は請求項６に記載の損傷診断システム。 - 前記処理装置は、前記推定演算を行って得られた前記被検体に対する損傷規模を表示する代わりに、特定した非対称モードまたは対称モードにおける前記周波数及び伝播時間の関係、算出した前記伝播時間の変化率、算出した前記減少量、または算出した前記増加量を表示することを特徴とする請求項１から請求項７のいずれか一項に記載の損傷診断システム。
Priority Applications (5)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2010058784A JP5629481B2 (ja) | 2010-03-16 | 2010-03-16 | 損傷診断システム |
| EP11157867.0A EP2369334B1 (en) | 2010-03-16 | 2011-03-11 | System and method for damage diagnosis |
| US13/047,517 US20110231112A1 (en) | 2010-03-16 | 2011-03-14 | System and Method for Damage Diagnosis |
| CN201110063840.9A CN102192954B (zh) | 2010-03-16 | 2011-03-16 | 损伤诊断系统 |
| US14/335,182 US10054568B2 (en) | 2010-03-16 | 2014-07-18 | System and method for damage diagnosis |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2010058784A JP5629481B2 (ja) | 2010-03-16 | 2010-03-16 | 損傷診断システム |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2011191230A JP2011191230A (ja) | 2011-09-29 |
| JP5629481B2 true JP5629481B2 (ja) | 2014-11-19 |
Family
ID=44010056
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2010058784A Active JP5629481B2 (ja) | 2010-03-16 | 2010-03-16 | 損傷診断システム |
Country Status (4)
| Country | Link |
|---|---|
| US (2) | US20110231112A1 (ja) |
| EP (1) | EP2369334B1 (ja) |
| JP (1) | JP5629481B2 (ja) |
| CN (1) | CN102192954B (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US12146855B2 (en) | 2019-01-15 | 2024-11-19 | Yamaha Hatsudoki Kabushiki Kaisha | Internal state detector and vehicle |
Families Citing this family (37)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8499632B1 (en) * | 2010-08-23 | 2013-08-06 | The Boeing Company | Characterizing anomalies in a laminate structure |
| US8788218B2 (en) | 2011-01-21 | 2014-07-22 | The United States Of America As Represented By The Secretary Of The Navy | Event detection system having multiple sensor systems in cooperation with an impact detection system |
| DE102012202189B3 (de) * | 2012-02-14 | 2013-02-07 | Sonotec Ultraschallsensorik Halle Gmbh | Verfahren und Vorrichtung zur Wavelet basierten Laufzeitmessung von Ultraschallimpulsen |
| GB2521661A (en) | 2013-12-27 | 2015-07-01 | Xsens As | Apparatus and method for measuring flow |
| CN104007175B (zh) * | 2014-05-09 | 2017-04-05 | 华南理工大学 | 一种悬臂柔性梁多裂缝损伤识别装置及方法 |
| CN103954697B (zh) * | 2014-05-20 | 2016-09-28 | 南京信息工程大学 | 一种基于分数阶微分的兰姆波去噪方法 |
| CN103971012B (zh) * | 2014-05-20 | 2017-01-18 | 南京信息工程大学 | 一种基于赛利斯模型和分数阶微分的兰姆波信号消噪方法 |
| US9791312B2 (en) | 2014-06-26 | 2017-10-17 | Baker Hughes Incorporated | OFDR system for localized vibration detection |
| JP6777301B2 (ja) * | 2015-03-10 | 2020-10-28 | 日本電気株式会社 | 検査装置及び検査方法 |
| JP6317708B2 (ja) * | 2015-06-12 | 2018-04-25 | 株式会社Ｓｕｂａｒｕ | 超音波探傷システム、超音波探傷方法及び航空機構造体 |
| CN105044212A (zh) * | 2015-06-30 | 2015-11-11 | 清华大学 | 一种多模态超声Lamb波复杂缺陷层析成像方法 |
| JP6329188B2 (ja) * | 2016-02-29 | 2018-05-23 | 株式会社Ｓｕｂａｒｕ | 超音波検査システム、超音波検査方法及び航空機構造体 |
| JP6109431B1 (ja) * | 2016-03-01 | 2017-04-05 | 三菱電機株式会社 | 超音波測定装置及び超音波測定方法 |
| EP3441755B1 (en) * | 2016-04-06 | 2025-07-02 | Subaru Corporation | Aircraft structural object and ultrasonic inspection method |
| CN108896935B (zh) * | 2016-04-08 | 2020-11-13 | 北京信息科技大学 | 一种利用超声脉冲诱发光栅变形测量磁场的方法 |
| KR101880194B1 (ko) * | 2016-05-30 | 2018-07-20 | 한국과학기술원 | 구조물의 진단 방법 및 진단 시스템 |
| CN106290566B (zh) * | 2016-08-04 | 2019-03-08 | 南京航空航天大学 | 非接触式评估复合材料结构的疲劳特性的方法及系统 |
| GB2545044B (en) | 2016-08-11 | 2018-02-07 | Guided Ultrasonics Ltd | Determining a thickness of a region of wall-or plate-like structure |
| EP3517445B1 (en) * | 2016-09-26 | 2022-08-24 | Subaru Corporation | Damage detection system and damage detection method |
| CN106441383B (zh) * | 2016-10-10 | 2019-08-23 | 山东大学 | 一种基于光纤光栅传感的主动激励式压滤机滤板健康检测系统 |
| CN108508097B (zh) * | 2017-02-28 | 2021-03-02 | 香港理工大学 | 一种基于光纤超声导波技术的铁轨裂纹监测系统 |
| CN110892261B (zh) | 2017-07-27 | 2023-05-30 | 株式会社斯巴鲁 | 超声波检查系统的制造方法 |
| CN108008021B (zh) * | 2018-01-09 | 2023-10-13 | 东莞理工学院 | 用于钢轨探伤的超声导波斜探头及其探伤方法 |
| KR102813799B1 (ko) * | 2018-05-04 | 2025-05-29 | 더 보잉 컴파니 | 음향 방출 시험을 위한 방법 및 장치 |
| CN110497425A (zh) * | 2018-05-17 | 2019-11-26 | 张家港市光武智能科技有限公司 | 一种基于兰姆波与瑞利波的机器人皮肤 |
| JP7131967B2 (ja) | 2018-05-30 | 2022-09-06 | 株式会社Subaru | 光検査システム、光検査方法及び航空機構造体 |
| JP6748675B2 (ja) * | 2018-06-29 | 2020-09-02 | 株式会社Ｓｕｂａｒｕ | 航空機管理システム |
| CN109298076B (zh) * | 2018-11-13 | 2021-05-04 | 中国核动力研究设计院 | 一种基于Lamb波的主动式阀门内漏损伤检测系统及方法 |
| CN110485296A (zh) * | 2019-08-15 | 2019-11-22 | 武汉理工大学 | 基于宏纤维复合材料正交异性钢桥面板疲劳振动减振方法 |
| CN111175379B (zh) * | 2020-01-21 | 2022-08-02 | 南京邮电大学 | 一种基于同步压缩小波变换的Lamb波板结构健康监测方法 |
| KR102323836B1 (ko) * | 2020-03-19 | 2021-11-08 | 부경대학교 산학협력단 | 객체 손상 검사 장치 및 이를 이용한 객체 손상 검사 방법 |
| CN111551630B (zh) * | 2020-04-23 | 2021-07-13 | 上海大学 | 一种基于空间-波数滤波器的损伤无波速定位方法 |
| CN111983032B (zh) * | 2020-08-25 | 2023-10-10 | 中国科学院电工研究所 | 一种光学元件损伤的在线监测方法及系统 |
| DE102020123975A1 (de) | 2020-09-15 | 2022-03-17 | Technische Universität Dresden, Körperschaft des öffentlichen Rechts | Vorrichtung und Messverfahren zur ortsaufgelösten in-situ Erfassung von Flüssigwasser auf einer Trägerstruktur sowie deren Verwendung in einem Brennstoffzellenstapel und in Plattenwärmeübertrager-Anwendungen |
| CN114428119B (zh) * | 2022-01-20 | 2023-11-21 | 重庆大学 | 一种用于反演各向异性特征的复合材料弹性常数的方法 |
| US20230334646A1 (en) * | 2022-04-14 | 2023-10-19 | Antech Systems, Inc. | System and Method for Mode Identification and Extraction in Guided Wave Signals and nondestructive methods for automatically sensing, locating, and/or characterizing problems in structures |
| CN115307774B (zh) * | 2022-08-10 | 2024-04-26 | 山东大学 | 一种基于超声导波的物体表面温度场检测方法及系统 |
Family Cites Families (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS59173753A (ja) * | 1983-03-23 | 1984-10-01 | Mitsubishi Electric Corp | 試料板体の性状を判定するための板波送受信装置 |
| IT1262407B (it) * | 1993-09-06 | 1996-06-19 | Finmeccanica Spa | Strumentazione utilizzante componenti in ottica integrata per la diagnostica di parti con sensori a fibra ottica inclusi o fissati sulla superficie. |
| US5814730A (en) * | 1996-06-10 | 1998-09-29 | Institute Of Paper Science And Technology And Georgia Institute Of Technology | Material characteristic testing method and apparatus using interferometry to detect ultrasonic signals in a web |
| JP3957276B2 (ja) * | 2001-05-14 | 2007-08-15 | 富士フイルム株式会社 | 超音波受信装置 |
| US7117742B2 (en) * | 2003-09-22 | 2006-10-10 | Advanced Structure Monitoring, Inc. | Sensors and systems for structural health monitoring |
| JP4377642B2 (ja) | 2003-09-26 | 2009-12-02 | 富士重工業株式会社 | 構造用複合材料の損傷探知システム |
| JP4470655B2 (ja) * | 2004-09-01 | 2010-06-02 | Ｊｆｅスチール株式会社 | 超音波によるスポット溶接部の評価方法及び装置 |
| US20090326834A1 (en) * | 2004-11-12 | 2009-12-31 | Sundaresan Mannur J | Systems, methods and computer program products for characterizing structural events |
| JP4742675B2 (ja) * | 2005-05-20 | 2011-08-10 | 三菱化学株式会社 | 筒状体内面付着層の厚さ測定方法 |
| JP2007023271A (ja) | 2005-06-16 | 2007-02-01 | Toyobo Co Ltd | ポリエステル重合触媒およびこれを用いて製造されたポリエステル並びにポリエステルの製造方法 |
| JP4804957B2 (ja) | 2006-02-27 | 2011-11-02 | 富士重工業株式会社 | 損傷探知システム |
| US8176786B2 (en) * | 2006-06-30 | 2012-05-15 | Carnegie Mellon University | Methods, apparatuses, and systems for damage detection |
| FR2905208B1 (fr) * | 2006-08-28 | 2008-12-19 | St Microelectronics Sa | Filtre a resonateurs a ondes de lamb couples. |
| US8419643B2 (en) * | 2006-09-21 | 2013-04-16 | Artann Laboratories Inc. | Ultrasonic method and apparatus for assessment of bone |
| US7809513B2 (en) * | 2007-04-16 | 2010-10-05 | Acellent Technologies, Inc. | Environmental change compensation in a structural health monitoring system |
| US7963165B2 (en) * | 2007-09-25 | 2011-06-21 | Los Alamos National Security, Llc | Non-contact feature detection using ultrasonic Lamb waves |
| CA2711894C (en) * | 2008-01-11 | 2016-10-18 | Pii Limited | Pipeline inspection apparatus and method using two different ultrasound wavemodes |
| JP4553047B2 (ja) * | 2008-03-12 | 2010-09-29 | セイコーエプソン株式会社 | ラム波型共振子及び発振器 |
-
2010
- 2010-03-16 JP JP2010058784A patent/JP5629481B2/ja active Active
-
2011
- 2011-03-11 EP EP11157867.0A patent/EP2369334B1/en active Active
- 2011-03-14 US US13/047,517 patent/US20110231112A1/en not_active Abandoned
- 2011-03-16 CN CN201110063840.9A patent/CN102192954B/zh active Active
-
2014
- 2014-07-18 US US14/335,182 patent/US10054568B2/en active Active
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US12146855B2 (en) | 2019-01-15 | 2024-11-19 | Yamaha Hatsudoki Kabushiki Kaisha | Internal state detector and vehicle |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2011191230A (ja) | 2011-09-29 |
| CN102192954B (zh) | 2016-07-20 |
| US20110231112A1 (en) | 2011-09-22 |
| US10054568B2 (en) | 2018-08-21 |
| EP2369334B1 (en) | 2023-06-07 |
| US20140330528A1 (en) | 2014-11-06 |
| CN102192954A (zh) | 2011-09-21 |
| EP2369334A1 (en) | 2011-09-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP5629481B2 (ja) | 損傷診断システム | |
| US8922789B2 (en) | Displacement measuring device and displacement measuring method | |
| US7470888B2 (en) | Damage detection system | |
| EP3112857B1 (en) | System comprising an aircraft structural object attached to an ultrasonic test system, and ultrasonic test method | |
| JP6235508B2 (ja) | 超音波探傷システム、超音波探傷方法及び航空機部品の製造方法 | |
| JP5586009B2 (ja) | 振動検出システム、該システムを用いた装置及び振動検出方法 | |
| JP5008182B2 (ja) | 衝撃探知システム | |
| Wang et al. | Phase-shifted fiber Bragg grating sensing network and its ultrasonic sensing application | |
| CN105277271A (zh) | 一种超声振动的相移光纤光栅传感探测系统及其应用 | |
| Soman et al. | Guided waves based damage localization based on mode filtering using fiber Bragg grating sensors | |
| JP2008139170A (ja) | 衝撃探知システム | |
| CN110554091B (zh) | 光学检查系统、光学检查方法以及航空器结构体 | |
| JP4214483B2 (ja) | Ｆｂｇ超音波センサの被検体への取り付け構造及び取り付け方法 | |
| Li et al. | Doppler effect-based fiber-optic sensor and its application in ultrasonic detection | |
| JP4565093B2 (ja) | 可動式ｆｂｇ超音波センサ | |
| Peters | Remote bonding of fiber Bragg grating sensors for Lamb wave measurements: a review | |
| Wee | Enhancement of ultrasonic lamb wave transfer to fiber bragg grating sensors through remote bonding | |
| Soman et al. | Optimal Design of a Sensor Network for Guided Wave-Based Structural Health Monitoring Using Acoustically Coupled Optical Fibers | |
| CN113447566A (zh) | 结构健全性诊断系统、诊断方法以及航空器结构体 | |
| Gopalakrishnan et al. | Computational techniques for damage detection, classification and quantification | |
| Lin et al. | Piezo-Optical Active Sensing With PWAS And FBG Sensors For Structural Health Monitoring | |
| Wee et al. | Sensitivity of contact-free fiber Bragg grating sensor to ultrasonic Lamb wave | |
| Okabe et al. | Delamination detection in fiber metal laminates using the mode conversion of Lamb waves | |
| Parali | Guided Wave Structural Health Monitoring Using Remote Excitation Through Thin Wire Waveguides | |
| Wandowski et al. | Study of Elastic Wave Sensing and Mode Separation Based on Piezoelectric and Fibre Bragg Grating Sensors |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20121203 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20130711 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20130723 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20130924 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20140415 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20140521 |
|
| 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: 20140924 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20141006 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 5629481 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |