JP2020085459A - クラック検出装置とその方法 - Google Patents

クラック検出装置とその方法 Download PDF

Info

Publication number
JP2020085459A
JP2020085459A JP2018214708A JP2018214708A JP2020085459A JP 2020085459 A JP2020085459 A JP 2020085459A JP 2018214708 A JP2018214708 A JP 2018214708A JP 2018214708 A JP2018214708 A JP 2018214708A JP 2020085459 A JP2020085459 A JP 2020085459A
Authority
JP
Japan
Prior art keywords
crack
frequencies
frequency
unit
impedance
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.)
Granted
Application number
JP2018214708A
Other languages
English (en)
Other versions
JP7174236B2 (ja
Inventor
美濃谷 直志
Naoshi Minoya
直志 美濃谷
恵里 松永
Eri Matsunaga
恵里 松永
昌人 中村
Masato Nakamura
昌人 中村
昌幸 津田
Masayuki Tsuda
昌幸 津田
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.)
Nippon Telegraph and Telephone Corp
Original Assignee
Nippon Telegraph and Telephone 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 Nippon Telegraph and Telephone Corp filed Critical Nippon Telegraph and Telephone Corp
Priority to JP2018214708A priority Critical patent/JP7174236B2/ja
Priority to US17/285,431 priority patent/US11454607B2/en
Priority to PCT/JP2019/043146 priority patent/WO2020100641A1/ja
Publication of JP2020085459A publication Critical patent/JP2020085459A/ja
Application granted granted Critical
Publication of JP7174236B2 publication Critical patent/JP7174236B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N27/00Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
    • G01N27/02Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance
    • G01N27/028Circuits therefor
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N27/00Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
    • G01N27/02Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance
    • G01N27/04Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance
    • G01N27/20Investigating the presence of flaws
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N27/00Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
    • G01N27/72Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating magnetic variables
    • G01N27/82Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating magnetic variables for investigating the presence of flaws
    • G01N27/90Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating magnetic variables for investigating the presence of flaws using eddy currents
    • G01N27/9046Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating magnetic variables for investigating the presence of flaws using eddy currents by analysing electrical signals
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N27/00Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
    • G01N27/72Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating magnetic variables
    • G01N27/82Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating magnetic variables for investigating the presence of flaws
    • G01N27/90Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating magnetic variables for investigating the presence of flaws using eddy currents
    • G01N27/9073Recording measured data

Landscapes

  • Chemical & Material Sciences (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Electrochemistry (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)

Abstract

【課題】センサの形状に限定されず、クラックによるインピーダンスの変化を容易に見分けることが可能なクラック検出装置を提供する。【解決手段】構造物200に貼り付けられる伝導体−絶縁体−伝導体の三層構造のセンサ部10と、所定の周波数の範囲を掃引して前記センサ部のインピーダンスが最大及び最小になる複数の周波数を取得する周波数特性取得部20と、複数の周波数の不均一性からクラックの有無を判定するクラック有無判定部30と、クラックの位置と周波数のずれの方向の関係を記録したクラック位置テーブル40と、クラック有無判定部30でクラックがあると判定された場合に、周波数特性取得部20で取得した2つの周波数の差を取って符号を求め、該符号でクラック位置テーブル40を参照してクラックの位置を検出するクラック位置検出部50とを備える。【選択図】図1

Description

本発明は、構造物のクラックを検出するクラック検出装置とその方法に関する。
構造物のクラックを検出するためには、構造物の在る現場に行き、目視、及び機器を用いて点検する方法がある(例えば非特許文献1)。この方法では、作業者を派遣するコストが生じると共に常時監視することが不可能である。
構造物のクラックを常時検出するためには、例えば非特許文献2に開示されたセンサを構造物に貼り付けて監視する方法がある。センサは、例えば、伝導体−絶縁体−伝導体の三層構造であり、そのインピーダンスの変化でクラックの有無を検出する。
コンクリート構造物のクラック検知ツール〔平成30年11月6日検索〕、インターネット(URL: http://www.kurabo.co.jp/chem/kk_crack_sensor.html) 原田隆郎、他2名、「PVDF フィルムセンサによるひび割れ検知に関する研究」、構造工学論文集Vol.59A(2013 年3 月)pp. 47-55.
しかしながら、センサの形状は、予め決まっているため種々の構造の様々な形状に適用するのが難しいという課題がある。また、センサの端子に水等の誘電体が付着すると、センサのインピーダンスが変化する。この水等の誘電体が付着することで生じるインピーダンスの変化及び経年による変化と、クラックが生じたことによるインピーダンスの変化を見分けるのが困難であるという課題がある。
本発明は、この課題に鑑みてなされたものであり、センサの形状に限定されず、クラックによるインピーダンスの変化を容易に見分けることが可能なクラック検出装置とその方法を提供することを目的とする。
本発明の一態様に係るクラック検出装置は、構造物のクラックを検出するクラック検出装置であって、前記構造物に貼り付けられる伝導体−絶縁体−伝導体の三層構造のセンサ部と、所定の周波数の範囲を掃引して前記センサ部のインピーダンスが最大及び最小になる複数の周波数を取得する周波数特性取得部と、複数の周波数の不均一性からクラックの有無を判定するクラック有無判定部と、クラックの位置と周波数のずれの方向の関係を記録したクラック位置テーブルと、前記クラック有無判定部でクラックがあると判定された場合に、前記周波数特性取得部で取得した2つの周波数の差を取って符号を求め、該符号で前記クラック位置テーブルを参照してクラックの位置を検出するクラック位置検出部と
を備えることを要旨とする。
また、本発明の一態様に係るクラック検出方法は、クラックを検出する対象の構造物に貼り付けられるセンサ部を含むクラック検出装置が実行するクラック検出方法であって、所定の周波数の範囲を掃引して前記センサ部のインピーダンスが最大及び最小になる複数の周波数を取得する周波数特性取得ステップと、複数の周波数の不均一性からクラックの有無を判定するクラック有無判定ステップと、前記クラック有無判定ステップでクラックがあると判定された場合に、前記周波数特性取得ステップで取得した2つの周波数の差を取って符号を求め、該符号でクラックの位置と周波数のずれの方向の関係を記録したクラック位置テーブルを参照してクラックの位置を検出するクラック位置検出ステップとを行うことを要旨とする。
本発明によれば、センサの形状に限定されず、クラックによるインピーダンスの変化を容易に見分けることが可能なクラック検出装置とその方法を提供することができる。
本発明の実施形態に係るクラック検出装置の機能構成例を示すブロック図である。 クラックがある場合のセンサ部の平面を模式的に示す図である。 クラックがある場合のセンサ部を伝送線で表した場合の回路モデルを示す図である。 センサ部のインピーダンスの周波数特性の例を模式的に示す図である。 クラックの位置を表すパラメータの例を示す図である。 図1に示すクラック位置テーブルの例とクラックの位置を示す図である。 図1に示す周波数特性取得部の概念的な構成を示すブロック図である。 図7に示す周波数特性取得部の具体的な機能構成例を示すブロック図である。
以下、本発明の実施形態について図面を用いて説明する。複数の図面中同一のものに
は同じ参照符号を付し、説明は繰り返さない。
(センサ部)
本発明の実施形態に用いるセンサ部について説明する。センサ部は、伝導体−絶縁体−伝導体の三層構造を、柔らかいフィルムで挟んだ構造である。
センサ部は、クラックの有無を検出する対象の構造物に貼り付けて使用される。構造物に貼り付けられたセンサ部の構造物側の伝導体を一方の端子、外側の伝導体を他方の端子とみなすと、センサ部は一方の端子対AB及び他方の端子対CDで表される伝送路の等価回路で表現することができる。
伝送路の一方の端子対ABからみたインピーダンスZtは以下の式で表せる。
Figure 2020085459
上式でZoは伝送線の特性インピーダンス、γは伝播定数(複素数)、Ltは伝送線の長さを表す。ここで一般的には伝搬定数γに含まれるαを減衰定数(実数)を無視して、位相定数(実数)βとj=(−1)0.5を用いて、γ=jβとして式(1)を展開すると以下の式を得る。
Figure 2020085459
伝送線の単位長あたりのインダクタンスをLo、容量をCo、円周率をπ、周波数をω/2πとするとβ=ω(Lo Co)0.5の関係が分布定数回路の理論で知られている。式(2)において、分母と分子に含まれるβLt=ω(Lo Co)0.5Ltを変数とする三角関数により、インピーダンスは周期的にピークとなる周波数特性となる。インピーダンスがピークとなる周期が伝送線の長さLtに関係するのは、端子対CDが電磁波の反射点となって端子対ABと端子対CDの間で定在波がたつためである。構造物にクラックが生じた時では、クラックが電磁波の反射点となり定在波のたつ周波数が変わるため、インピーダンスがピークとなる周波数が変化する。
検知対象である構造物が屋外に設置された場合、雨滴等によりセンサ部の端に水が付着することがある。また、構造物が屋内にある場合でも結露などによりセンサ部の端に水が付着することがある。この場合のセンサ部は、伝送線の他方の端子対CDに水の付着による寄生容量Cpが接続された回路モデルで表される。この時の端子対ABからみたインピーダンスはγ=jβとして以下の式で表される。
Figure 2020085459
式(3)において分母と分子に含まれる三角関数の変数はβLt+ωCpZoとなり、インピーダンスがピークとなる定在波のたつ周波数は寄生容量Cpにも依存する。式(2)で表されるCpがないときの定在波のたつ周波数と比較して、Cpがあるときの定在波のたつ周波数は変化する。
また、絶縁体の誘電率の温度・湿度等での変化や経年での変化でも定在波のたつ周波数は変化する。
このように、構造物にクラックが生じた場合だけでなく、センサ部の端に水等の誘電体が付着した場合や誘電体の温度・湿度や経年による変化でも定在波のたつ周波数が変化する。このため、センサ部の端に水等の誘電体が付着した場合および誘電体の温度・湿度や経年による変化の場合と、クラックが生じた場合を見分ける必要がある。
本発明の実施形態に係るクラック検出装置は、構造物にクラックが生じたことによるセンサ部のインピーダンスの変化とそれ以外によるインピーダンスの変化の見分けを容易にし、センサ部の形状による影響を受けない構成である。
(クラック検出装置)
図1は、本発明の実施形態に係るクラック検出装置の機能構成例を示すブロック図である。図1に示すクラック検出装置100は、センサ部10、周波数特性取得部20、クラック有無判定部30、クラック位置テーブル40、及びクラック位置検出部50を備える。
センサ部10は、構造物200に貼り付けられる伝導体−絶縁体−伝導体の三層構造である。センサ部10は、例えば、PVDFフィルムで構成することができる。
周波数特性取得部20は、所定の周波数の範囲を掃引してセンサ部10のインピーダンスが最大及び最小になる複数の周波数を取得する。周波数特性取得部20は、例えば、インピーダンスアナライザ等の一般的な計測器で構成することもできる。又は、FFTアナライザで構成してもよい。FFTアナライザでセンサ部10の周波数特性を測定し、複数の直列共振周波数と並列共振周波数を測定するようにしてもよい。具体例については後述する。
クラック有無判定部30は、周波数特性取得部20で取得した複数の周波数の不均一性からクラックの有無を判定する。センサ部10にクラックが無い場合は、周波数特性取得部20で取得した複数の周波数のそれぞれの間隔は一定である。センサ部10にクラックが生じた場合は、複数の周波数のそれぞれの間隔にバラツキが生じる。よって、複数の周波数の不均一性からクラックの有無を判定することができる。また、複数の周波数の不均一性からクラックの有無を判定するので、センサ部10の大きさが限定されない。つまり、任意の大きさのセンサ部10を用いることができる。詳しくは後述する。
クラック位置テーブル40は、センサ部10上のクラックの位置と周波数のずれの方向の関係を記録したテーブルである。クラック位置テーブル40は、予め設けられる。具体例は後述する。
クラック位置検出部50は、クラック有無判定部30でクラックがあると判定された場合に、周波数特性取得部20で取得した2つの周波数の差を取って符号を求め、該符号でクラック位置テーブル40を参照してクラックの位置を検出する。クラックが生じた位置によって、センサ部10のインピーダンス特性は変化するので、2つの周波数の差を取り、予め設定された基準と比較することでクラックの位置を特定することができる。詳しくは後述する。
このように本実施形態に係るクラック検出装置100は、構造物200のクラックを検出するクラック検出装置であって、構造物200に貼り付けられる伝導体−絶縁体−伝導体の三層構造のセンサ部10と、所定の周波数の範囲を掃引してセンサ部10のインピーダンスが最大及び最小になる複数の周波数を取得する周波数特性取得部20と、複数の周波数の不均一性からクラックの有無を判定するクラック有無判定部30と、クラックの位置と周波数のずれの方向の関係を記録したクラック位置テーブル40と、クラック有無判定部30でクラックがあると判定された場合に、周波数特性取得部20で取得した2つの周波数の差を取って符号を求め、該符号でクラック位置テーブル40を参照してクラックの位置を検出するクラック位置検出部50とを備える。
これにより、センサの形状に限定されず、クラックによるインピーダンスの変化を容易に見分けることが可能になる。
以下でクラックを検知する動作の原理を詳しく説明する。図2にクラックがある場合のセンサ部10の模式図を示し、図3にクラックがある場合の回路モデルを示す。図2の模式図では上部の伝導体の左端からLt1で右端からLt2の箇所に間隙がδlのクラックが生じている状態を示している。図3の回路モデルにおいてCcはクラックの間隙の部分に存在する浮遊容量を表し、L1とC1図2においてつながっている部分を伝送線とみなした場合の単位長さあたりのインダクタンスと容量を表している。長さLt1とLt2の伝送線の特性インピーダンスをZo(実数)、伝播定数をγ(複素数)としている。
端子AとB間の電位差をVab、端子2と2‘間の電位差をV2、端子3と3’間の電位差をV3、端子CとD間の電位差をVcdとする。また、端子Aにおいて伝送線に沿って部分に流れる電流をIab、端子2において伝送線部分の電流をI2、端子3において伝送線部分の電流をI2、端子Cにおいて伝送線に沿って流れる電流をIcdとする。分布定数回路の理論に従い、端子Bには逆方向にIabが流れ、同様に端子2‘には逆方向にI2、端子3’には逆方向にI3、端子Dには逆方向にIcdが流れるとする。本回路モデルでは以下の回路方程式が成立する。
Figure 2020085459
ここで、K1とK2は長さLt1の伝送線の境界条件によって決定される定数であり、K3とK4は長さLt2の伝送線の境界条件によって決定される定数である。また上式でeは自然対数の底(ネイピア数)を表す。
端子CDが開放であることからIcd=0となり、式(13)からK3とK4の関係は以下の式となる。
Figure 2020085459
式(14)を式(10)に代入すると以下の式が得られる。
Figure 2020085459
同様に式(14)を式(11)に代入すると以下の式が得られる。
Figure 2020085459
式(16)で式(15)を割るとV3とI3の関係式が以下のようになる。
Figure 2020085459
式(9)からV2はV3とI3を用いて以下の式で表される。
Figure 2020085459
式(8)からI2はV3とI3を用いて以下の式で表される。
Figure 2020085459
式(18)に代入する。
Figure 2020085459
ここで式(20)から式(22)のパラメータを用いると式(18)と式(19)はそれぞれ式(23)と式(24)のように表される。
Figure 2020085459
式(6)と式(7)からK1とK2は式(25)と式(26)で表される。
Figure 2020085459
ここでLp=L1+L2、Ln=-L1+L2すると、L1=(Lp-Ln)/2、L2=(Lp+Ln)/2となり、式(25)と式(26)は以下の式で表される。
Figure 2020085459
式(17)、式(27)、式(28)を式(4)に代入するとVabに関して以下の式が得られる。
Figure 2020085459
ここで、以下のパラメータを導入すると、式(29)は式(32)のように展開できる。
Figure 2020085459
式(17)、式(27)、式(26)を式(5)に代入するとIabに関して以下の式が得られる。
Figure 2020085459
ここで式(30)と式(31)を用いて展開するとIabの表式は以下のようになる。
Figure 2020085459
端子ABからみたインピーダンスZab=Vab/Iabは式(32)と式(33)から以下の式で表される。
Figure 2020085459
式(2)と式(3)の式展開と場合と同様に、βを位相定数(実数)、jを(−1)0.5として、γ=jβとすると、Zabの分子は以下のように展開できる。
Figure 2020085459
Zabの分母は以下のように展開できる。
Figure 2020085459
クラックが入った場合では、伝送線のインピーダンスZabの式の分母と分子に含まれる三角関数の変数は、{β(Lt1+Lt2)+θz}とβ(Lt1-Lt2)の2種となる。クラックが入っておらず伝送線の端の付着物(寄生容量Cp)がない場合ではインピーダンスの式に含まれる三角関数の変数は式(A-2)よりβLtの1種であり、クラックが入っておらず伝送線の端の付着物(寄生容量Cp)がある場合のインピーダンスの式に含まれる三角関数の変数は式(A-3)より(βLt+ωCpZo)の1種である。この考察からクラックの有無によりインピーダンスの式に含まれる三角関数の変数の数が変化することが分かる。ここでLtは伝送線の長さを表し、Lt=Lt1+Lt2+δlの関係がある。
クラックが無い場合でのインピーダンスのピークの周期は、β=ω(Lo Co)0.5の関係とインピーダンスの式に含まれる三角関数の変数から考えて、Cpが無い場合で(Lo Co)0.5Lt、Cpがある場合で(Lo Co)0.5Lt+CpZoと周波数に依存せず一定である。
これに対し、クラックがある場合において、2種の変数のうちβ(Lt1-Lt2)の成分に関しては、インピーダンスがピークとなる周期が(Lo Co)0.5(Lt1-Lt2)となる。{β(Lt1+Lt2)+θz}の成分に関しては、以下のように考察できる。
以下の式で表されるパラメータθzはLp=L1+L2、式(20)、式(21)を用いて式(38)のように表される。
Figure 2020085459
構造物に生じたクラックがまだ初期段階でクラック間隙δlが非常に短い場合には以下の近似が成立する。
Figure 2020085459
さらに周波数ω/2πが低い領域では式(39)は以下の式(40)となり、θzは周波数に正比例する。
Figure 2020085459
従って、構造物に生じたクラックがまだ初期段階でクラック間隙δlが非常に短く周波数が低い領域における{β(Lt1+Lt2)+θz}の成分の周期は{(Lo Co)0.5(Lt1+Lt2)+(L1+C1 Zo2)δl/(2Zo) }となる。
また、クラックが初期段階でのLp、Ln、|Az|と以下の式で与えられる。
Figure 2020085459
式(42)、式(43)よりクラックがまだ初期段階で間隙δlが非常に短く周波数が低い領域ではLnと|Az|は周波数に対し一定となる。クラックが生じている個所の伝導体の幅は、クラックのない伝導体の幅と比較して狭いためLo<L1、Co>C1となる。式(42)よりLnは負となる。
次に、クラックがまだ初期段階で間隙δlが非常に短く周波数が低い領域でアドミタンスとインピーダンスがピークとなる周波数について考察する。アドミタンスはインピーダンスの逆数のため、アドミタンスのピークとなる周波数はZabの分子である式(35)がゼロとなる周波数である。式(41)から式(43)を代入すると以下の式を得る。
Figure 2020085459
クラックがまだ初期段階で間隙δlが非常に短く周波数が低い領域ではZo>>ω(L1-LoC1/Co)δl/2となり式(44)の第1項が支配的となるが、Zabの分子がゼロ近傍を議論することから式(44)では第2項も無視せずに考慮することとした。
支配的な式(44)の第1項がゼロとなる周波数の近傍に式(44)がゼロとなる周波数が存在すると考えられる。式(44)の第1項のcos関数の大きさがゼロとなるとき第1項がゼロとなる。円周率をπとすると、cos関数の変数がπ/2の奇数倍となるときcos関数の大きさはゼロとなる。cos関数の性質によりcos関数の大きさがゼロになる近傍では式(44)の第1項のcos関数は以下の式で近似できる。
Figure 2020085459
式(45)と式(46)とβ=ω(LoCo)0.5を用いて、式(44)がゼロとなるすなわちアドミタンスのピークとなる周波数ω/2πは以下の条件を満たす。
Figure 2020085459
式(47)と式(48)をまとめると以下の式となる。
Figure 2020085459
インピーダンスのピークとなる周波数はZabの分母である式(36)がゼロとなる周波数である。式(41)から式(43)を代入すると以下の式を得る。
Figure 2020085459
クラックがまだ初期段階で間隙δlが非常に短く周波数が低い領域ではZo>>ω(L1-LoC1/Co)δl/2となり式(50)の第1項が支配的となるが、Zabの分母がゼロ近傍を議論することから式(50)では第2項も無視せずに考慮することとした。
支配的な式(50)の第1項がゼロとなる周波数の近傍に式(50)がゼロとなる周波数が存在すると考えられる。式(50)の第1項のsin関数の大きさがゼロとなるとき第1項がゼロとなる。Sin関数の変数がπの整数倍となるときsin関数の大きさはゼロとなる。sin関数の性質によりsin関数の大きさがゼロになる近傍では式(50)の第1項のsin関数は以下の式で近似できる。
Figure 2020085459
式(51)と式(52)とβ=ω(LoCo)0.5を用いて、式(44)がゼロとなるすなわちアドミタンスのピークとなる周波数ω/2πは以下の条件を満たす。
Figure 2020085459
式(53)と式(54)をまとめると以下の式となる。
Figure 2020085459
さらに式(49)と式(55)をまとめると以下の式を得る。
Figure 2020085459
式(56)において周波数fm=ωm/2πを用いた。
図4にインピーダンスの周波数特性の例を示す。図4の周波数f1やf3のインピーダンスの極小値はアドミタンスのピークを表す。アドミタンスとインピーダンスのピークの周波数を直流側から番号を付与すると式(56)中のmと同じになる。
クラックの無いδl=0の場合では式(2)や式(3)について議論したとおりアドミタンスとインピーダンスのピークの周波数は周期的であるため、隣り合うアドミタンスとインピーダンスのピークの周波数の差は一定である。このため、アドミタンスとインピーダンスのピークの周波数を直流側から付与した番号mでその周波数を規格化した周波数は一定となる。式(56)においてもδl=0とすると右辺に周波数に依存する右辺の第3項がゼロとなるため、mで規格化した周波数fm/mは一定となる。
クラックのある場合では、fm/mは一定ではなくなる。このため、アドミタンスとインピーダンスのピークとなる周波数を数点取得し、mで規格化した周波数の一致・不一致でクラックの有無が判定できる。
次にクラックの位置について考察する。クラックの位置をLt1とLt2のいずれかで表すことができるが、以下の式で表される全体の長さLt1+Lt2で規格化したパラメータを導入しても表すことができる。
Figure 2020085459
図5にパラメータLocとトランスデューサ上の位置の関係を示す。アドミタンスとインピーダンスのピークとなる周波数を計測中に全体の長さLt1+Lt2は変化しないため式(57)の分母は一定である。トランスデューサの中心ではLt1=Lt2となりLoc=0となる。LocはLt1と正比例するためクラックの位置が図5上の右側になる(Lt1が大きくなる)とLocは0より大きくなる。また、LocはLt2に-1を乗じた数に比例するためクラックの位置が図5上の左側になる(Lt2が大きくなる)とLocは0より小さくなる。
式(56)の右辺第3項のcos関数の変数にLt1-Lt2があるため、第3項のみを抽出すればパラメータLocとの関係式が得られると考えられる。クラックがあるためfm/mは一定ではないが、f1の近傍にあると考え、さらにf1はクラックが無い時にアドミタンスがピークになる最低の周波数fbの近傍にあると考えて以下の式でfmを表す。
Figure 2020085459
fbは1/4波長で定在波がたったときの周波数のため、トランスデューサの全体の長さLt1+Lt2と電磁波の速度(LoCo)-0.5とに以下の関係が成立する。
Figure 2020085459
式(58)を用いて式(56)右辺の第3項のcos関数を展開すると以下の式を得る。
Figure 2020085459
式(60)の展開では2πδf(LoCo)0.5 (Lt1-Lt2)<<1とした。式(60)に式(59)を代入すると以下の式を得る。
Figure 2020085459
式(61)において第1項がゼロより十分大きくδf/fb<<1とすれば式(56)右辺の第3項のcos関数を以下の式で近似できる。
Figure 2020085459
式(62)を用いて式(56)は以下の式となる。
Figure 2020085459
任意のma番目の周波数とmb番目の周波数についてmで規格化した周波数の逆数の差をとるとLocに依存する第3項を取り出すことができ以下の式を得る。
Figure 2020085459
式(64)の右辺は周波数がなくても導出可能であり、maとmbは予め決定できる数である。また、任意の長さLt1+Lt2のトランスデューサに対してLocの範囲は-1<Loc<1であるため、BLを除けばどのような長さのトランスデューサ対しても式(64)の右辺は予め導出できる。前述の式(42)で説明したとおりBL>0であるため、右辺の正負の符号は左辺と同じである。式(64)の左辺のfmaとfmbは、アドミタンスとインピーダンスの周波数依存性の測定により得られる。このため、予めLocとma、mbの組合せに対する式(65)で表される第1の評価関数の結果をクラック位置テーブルとして保持しておき、アドミタンスとインピーダンスの周波数依存性の測定結果から得られた式(66)で表される第2の評価関数Evmの符号、ma、mbの組合せとクラック位置テーブル40からクラックの位置であるLocの領域を逆引きすることができる。
Figure 2020085459
式(65)と式(66)でsgnは符号関数である。
アドミタンスとインピーダンスの個数を4とした場合のクラック位置テーブル40の例を図6(a)に示す。例えば、はじめにma=3とmb=1でのEvmの符号を評価しEvmが正の時では、次にma=4とmb=3でのEvmを評価する。ma=4とmb=3でのEvmが正の場合では図9(b)に示すトランスデューサ上の(1)の領域にクラックがあることとなる。
上述のように、予めLocとma、mbの組合せに対する式(64)で表される評価関数の結果をクラック位置テーブル40として保持しておき、アドミタンスとインピーダンスの周波数依存性の測定結果から得られた式(65)で表されるEvmの符号、ma、mbの組合せとクラック位置テーブル40からクラックの位置であるLocの領域を逆引きすることができる。
周波数特性取得部20の具体例について説明する。
周波数特性取得部20は、所定の周波数の範囲を掃引してセンサ部10のインピーダンスが最大及び最小になる複数の周波数を取得する。図7に周波数特性取得部20の周波数可変信号源21、負荷部22とセンサ部10の回路モデルを示す。周波数可変信号源21から出力された電圧信号は負荷部22(抵抗Rs)を介してセンサ部10に印加される。センサ部10のインピーダンスが大きい時では端子A-B間の電圧Vaも大きいため、Vaがピークとなる周波数はインピーダンスのピークとなる周波数と一致する。センサ部10のアドミタンスが大きい時では電流Isが大きくなる。電流Isが大きくなると負荷部22での電圧降下が大きくなるため、Vs-Vaを測定すればアドミタンスの大きさがわかる。従ってVs-Vaがピークとなる周波数とアドミタンスがピークとなる周波数は一致する。
図8に周波数特性取得部20の構成例を示す。端子Aからの信号は直交復調部24と差動演算部23に入力される。直交復調部24では周波数可変信号源21から入力される基準信号の同相成分と直交成分の振幅に端子Aからの信号を分離し周波数特性算出部26に出力する。差動演算部23では負荷部22からの信号と端子Aからの信号の差をとって直交復調部25に出力する。直交復調部25では基準信号の同相成分と直交成分の振幅に差動演算部23からの入力信号を分離し周波数特性算出部26に出力する。周波数特性算出部26では直交復調部24の出力信号からインピーダンスのデータを生成し、直交復調部25の出力信号からアドミタンスのデータを生成して周波数特性データとしてクラック有無判定部30に送信する。
所定の周波数範囲の周波数に対応した周波数特性データ取得したのち、アドミタンスとインピーダンスがピークとなる周波数を抽出するとともに低周波側から周波数に対応した番号を付与する。
この後、アドミタンスとインピーダンスのピークの周波数を番号で規格化した複数の周波数が所定の不均一さ未満であればクラックが無いと判定する。規格化した複数の周波数が所定の不均一さ以上である場合には、2種の番号とピークの周波数を用いて式(66)で表される第2の評価関数で評価した結果とクラック位置テーブル40を照合してクラックの位置を決定する。クラックの有無の判定結果とクラックの位置の情報は、通信部(図示せず)を介して外部へ送信するようにしても良い。
以上説明したように本実施形態に係るクラック検出装置100は、構造物200のクラックを検出するクラック検出装置であって、構造物200に貼り付けられる伝導体−絶縁体−伝導体の三層構造のセンサ部10と、所定の周波数の範囲を掃引してセンサ部10のインピーダンスが最大及び最小になる複数の周波数を取得する周波数特性取得部20と、複数の周波数の不均一性からクラックの有無を判定するクラック有無判定部30と、クラックの位置と周波数のずれの方向の関係を記録したクラック位置テーブル40と、クラック有無判定部30でクラックがあると判定された場合に、周波数特性取得部20で取得した2つの周波数の差を取って符号を求め、該符号で前記クラック位置テーブル40を参照してクラックの位置を検出するクラック位置検出部50とを備える。
また、クラック有無判定部30は、インピーダンスが最大及び最小になる複数の周波数を低い周波数から1から順番に付番し、該付番した番号でそれぞれの周波数を正規化した正規化周波数を求め、複数の正規化周波数のそれぞれが同じであればクラックが無いと判定し、複数の正規化周波数が得られればクラックありと判定する。この判定方法によれば、周波数を番号で規格化した複数の周波数が所定の不均一さでクラックの有無を判定するので、ゼンサ部10のセンサの形状の影響を受けない。つまり、センサはどのような形状で有っても構わない。
このように本実施形態に係るクラック検出装置100は、クラックによるインピーダンスの変化を容易に見分けることを可能にする。なお、本発明は、上記の実施形態に限定されるものではなく、特許請求の範囲に記載された要旨の範囲内で変形が可能である。
10:センサ部
20:周波数特性取得部
21:周波数可変信号源
22:負荷部
23:差動演算部
24,25:直交復調部
26:周波数特性算出部
30:クラック有無判定部
40:クラック位置テーブル
50:クラック位置検出部

Claims (3)

  1. 構造物のクラックを検出するクラック検出装置であって、
    前記構造物に貼り付けられる伝導体−絶縁体−伝導体の三層構造のセンサ部と、
    所定の周波数の範囲を掃引して前記センサ部のインピーダンスが最大及び最小になる複数の周波数を取得する周波数特性取得部と、
    複数の周波数の不均一性からクラックの有無を判定するクラック有無判定部と、
    クラックの位置と周波数のずれの方向の関係を記録したクラック位置テーブルと、
    前記クラック有無判定部でクラックがあると判定された場合に、前記周波数特性取得部で取得した2つの周波数の差を取って符号を求め、該符号で前記クラック位置テーブルを参照してクラックの位置を検出するクラック位置検出部と
    を備えることを特徴とするクラック検出装置。
  2. 前記クラック有無判定部は、
    インピーダンスが最大及び最小になる複数の周波数を低い周波数から1から順番に付番し、該付番した番号でそれぞれの周波数を正規化した正規化周波数を求め、複数の正規化周波数のそれぞれが同じであればクラックが無いと判定し、複数の正規化周波数が得られればクラックありと判定する
    ことを特徴とする請求項1に記載のクラック検出装置。
  3. クラックを検出する対象の構造物に貼り付けられるセンサ部を含むクラック検出装置が実行するクラック検出方法であって、
    所定の周波数の範囲を掃引して前記センサ部のインピーダンスが最大及び最小になる複数の周波数を取得する周波数特性取得ステップと、
    複数の周波数の不均一性からクラックの有無を判定するクラック有無判定ステップと、、
    前記クラック有無判定ステップでクラックがあると判定された場合に、前記周波数特性取得ステップで取得した2つの周波数の差を取って符号を求め、該符号でクラックの位置と周波数のずれの方向の関係を記録したクラック位置テーブルを参照してクラックの位置を検出するクラック位置検出ステップと
    を行うことを特徴とするクラック検出方法。

JP2018214708A 2018-11-15 2018-11-15 クラック検出装置とその方法 Active JP7174236B2 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2018214708A JP7174236B2 (ja) 2018-11-15 2018-11-15 クラック検出装置とその方法
US17/285,431 US11454607B2 (en) 2018-11-15 2019-11-01 Crack detection device and method therefor
PCT/JP2019/043146 WO2020100641A1 (ja) 2018-11-15 2019-11-01 クラック検出装置とその方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2018214708A JP7174236B2 (ja) 2018-11-15 2018-11-15 クラック検出装置とその方法

Publications (2)

Publication Number Publication Date
JP2020085459A true JP2020085459A (ja) 2020-06-04
JP7174236B2 JP7174236B2 (ja) 2022-11-17

Family

ID=70731537

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2018214708A Active JP7174236B2 (ja) 2018-11-15 2018-11-15 クラック検出装置とその方法

Country Status (3)

Country Link
US (1) US11454607B2 (ja)
JP (1) JP7174236B2 (ja)
WO (1) WO2020100641A1 (ja)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09243584A (ja) * 1996-03-07 1997-09-19 E Syst Inc 構造物の一区域の試験及び検査を行う装置及び方法
JPH10293107A (ja) * 1997-02-19 1998-11-04 Tokin Corp 積層セラミックコンデンサの内部欠陥検査方法
JP2011174823A (ja) * 2010-02-25 2011-09-08 Railway Technical Research Institute 亀裂監視装置及び亀裂監視方法
CN103995036A (zh) * 2014-06-09 2014-08-20 河南理工大学 一种利用电化学阻抗谱实时监测水泥基材料裂缝的方法
JP2014190761A (ja) * 2013-03-26 2014-10-06 Railway Technical Research Institute 亀裂監視材及び亀裂監視装置

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4644259A (en) * 1984-10-31 1987-02-17 North American Philips Corporation Nondestructive testing of multilayer ceramic capacitors
JP2004028907A (ja) * 2002-06-27 2004-01-29 Shiyoui Ko 構造物健全度評価装置及びその装置を用いた構造物の健全度評価方法
JP6882688B2 (ja) * 2018-05-14 2021-06-02 日本電信電話株式会社 クラック検知システム及びクラック検知方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09243584A (ja) * 1996-03-07 1997-09-19 E Syst Inc 構造物の一区域の試験及び検査を行う装置及び方法
JPH10293107A (ja) * 1997-02-19 1998-11-04 Tokin Corp 積層セラミックコンデンサの内部欠陥検査方法
JP2011174823A (ja) * 2010-02-25 2011-09-08 Railway Technical Research Institute 亀裂監視装置及び亀裂監視方法
JP2014190761A (ja) * 2013-03-26 2014-10-06 Railway Technical Research Institute 亀裂監視材及び亀裂監視装置
CN103995036A (zh) * 2014-06-09 2014-08-20 河南理工大学 一种利用电化学阻抗谱实时监测水泥基材料裂缝的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
松永恵理: "インフラ劣化検知を目的とした塗装式センサの基礎検討", 第65回材料と環境討論会講演集, JPN6022011581, 12 October 2018 (2018-10-12), pages 349 - 350, ISSN: 0004736772 *

Also Published As

Publication number Publication date
US11454607B2 (en) 2022-09-27
WO2020100641A1 (ja) 2020-05-22
US20210364459A1 (en) 2021-11-25
JP7174236B2 (ja) 2022-11-17

Similar Documents

Publication Publication Date Title
CN104090214B (zh) 一种电缆故障检测及老化分析方法
US7174255B2 (en) Self-processing integrated damage assessment sensor for structural health monitoring
US20130204555A1 (en) Method and Apparatus for Electrically Locating a Fault in a Cable
CN104808195B (zh) 用于通过分析处理阻抗包络线识别超声换能器的故障的方法和设备
US20100211339A1 (en) Diagnostic method for determining deformations in a transformer winding
CN1666109A (zh) 时-频域反射仪的装置及方法
JP2011511679A5 (ja)
Giaquinto et al. Criteria for automated estimation of time of flight in TDR analysis
US20170219454A1 (en) Leak identification in a mass transport network
CN111551790B (zh) 用于确定用于电力系统的绝缘的介电质等效电路图的元件的方法和设备
WO2020100641A1 (ja) クラック検出装置とその方法
Gelman et al. Fatigue crack diagnostics: A comparison of the use of the complex bicoherence and its magnitude
JP6882688B2 (ja) クラック検知システム及びクラック検知方法
EP3298419A1 (fr) Procede d&#39;analyse d&#39;un cable, impliquant un traitement d&#39;amplification de la signature d&#39;un defaut non franc
KR102061616B1 (ko) 비선형 스펙트럼 상관관계를 이용한 구조물 진단 방법 및 구조물 진단 시스템
Trentadue et al. Detecting damage through the processing of dynamic shapes measured by a PSD-triangular laser sensor
Nissilä et al. Diagnosing simultaneous faults using the local regularity of vibration signals
US11435309B2 (en) Method and apparatus for detecting ply inconsistencies within composite objects
JP4550464B2 (ja) 地絡点標定方法および装置
Giaquinto et al. Accuracy analysis in the estimation of ToF of TDR signals
KR101188161B1 (ko) 전자 장치의 측정 방법 및 측정 장치
Scandurra et al. Accurate QTF sensing approach by means of narrow band spectral estimation
CN116015432B (zh) 基于光感和遥感的光缆监控方法、装置、设备及存储介质
Cui et al. Fatigue crack localisation based on empirical mode decomposition and pre-selected entropy
US11624687B2 (en) Apparatus and method for detecting microcrack using orthogonality analysis of mode shape vector and principal plane in resonance point

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20210218

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20220329

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20220525

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: 20221004

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20221017

R150 Certificate of patent or registration of utility model

Ref document number: 7174236

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150