ホーム > 特許ランキング > 日本電信電話株式会社
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】
(24)【登録日】2022-11-09
(45)【発行日】2022-11-17
(54)【発明の名称】クラック検出装置とその方法
(51)【国際特許分類】
G01N 27/02 20060101AFI20221110BHJP
【ＦＩ】
G01N27/02 Z
【請求項の数】
2
(21)【出願番号】P 2018214708
(22)【出願日】2018-11-15
(65)【公開番号】P2020085459
(43)【公開日】2020-06-04
【審査請求日】2021-02-18
(73)【特許権者】
【識別番号】000004226
【氏名又は名称】日本電信電話株式会社
(74)【代理人】
【識別番号】100083806
【弁理士】
【氏名又は名称】三好 秀和
(74)【代理人】
【識別番号】100129230
【弁理士】
【氏名又は名称】工藤 理恵
(72)【発明者】
【氏名】美濃谷 直志
(72)【発明者】
【氏名】松永 恵里
(72)【発明者】
【氏名】中村 昌人
(72)【発明者】
【氏名】津田 昌幸
【審査官】田中 洋介
(56)【参考文献】
【文献】特開平０９－２４３５８４（ＪＰ，Ａ）
【文献】特開２０１１－１７４８２３（ＪＰ，Ａ）
【文献】特開２０１４－１９０７６１（ＪＰ，Ａ）
【文献】特開平１０－２９３１０７（ＪＰ，Ａ）
【文献】中国特許出願公開第１０３９９５０３６（ＣＮ，Ａ）
【文献】松永恵理，インフラ劣化検知を目的とした塗装式センサの基礎検討，第６５回材料と環境討論会講演集，2018年10月12日，pp.349-350
(58)【調査した分野】(Int.Cl.，ＤＢ名)
Ｇ０１Ｎ ２７／００－２７／２４
ＪＳＴＰｌｕｓ／ＪＳＴ７５８０（ＪＤｒｅａｍＩＩＩ）
(57)【特許請求の範囲】
【請求項1】
構造物のクラックを検出するクラック検出装置であって、前記構造物に貼り付けられる伝導体－絶縁体－伝導体の三層構造のセンサ部と、
所定の周波数の範囲を掃引して前記センサ部のインピーダンスが最大及び最小になる複数の周波数を取得する周波数特性取得部と、
複数の前記周波数の不均一性からクラックの有無を判定するクラック有無判定部と、
クラックの位置と周波数のずれの方向の関係を記録したクラック位置テーブルと、
前記クラック有無判定部でクラックがあると判定された場合に、前記周波数特性取得部で取得した複数の前記周波数の内の２つの周波数の複数の組み合わせと、各組み合わせにおける２つの周波数の差の符号と、に対応するクラックの位置を、前記クラック位置テーブルを参照して検出するクラック位置検出部と、を備え、
前記クラック有無判定部は、
前記周波数特性取得部で取得した複数の前記周波数を、低い周波数から１から順番に付番し、該付番した番号でそれぞれの周波数を除算した正規化周波数を複数算出し、複数の前記正規化周波数が一定の場合はクラックが無いと判定し、複数の前記正規化周波数が一定でない場合はクラックありと判定する
ことを特徴とするクラック検出装置。
【請求項2】
構造物のクラックを検出するクラック検出装置が行うクラック検出方法であって、前記クラック検出装置は、
前記構造物に貼り付けられる伝導体－絶縁体－伝導体の三層構造のセンサ部を備え、
所定の周波数の範囲を掃引して前記センサ部のインピーダンスが最大及び最小になる複数の周波数を取得する周波数特性取得ステップと、
複数の前記周波数の不均一性からクラックの有無を判定するクラック有無判定ステップと、
前記クラック有無判定ステップでクラックがあると判定された場合に、前記周波数特性取得ステップで取得した複数の前記周波数の内の２つの周波数の複数の組み合わせと、各組み合わせにおける２つの周波数の差の符号と、に対応するクラックの位置を、クラックの位置と周波数のずれの方向の関係を記録したクラック位置テーブルを参照して検出するクラック位置検出ステップと、を行い、
前記クラック有無判定ステップは、
前記周波数特性取得ステップで取得した複数の前記周波数を、低い周波数から１から順番に付番し、該付番した番号でそれぞれの周波数を除算した正規化周波数を複数算出し、複数の前記正規化周波数が一定の場合はクラックが無いと判定し、複数の前記正規化周波数が一定でない場合はクラックありと判定する
ことを特徴とするクラック検出方法。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、構造物のクラックを検出するクラック検出装置とその方法に関する。
【背景技術】
【0002】
構造物のクラックを検出するためには、構造物の在る現場に行き、目視、及び機器を用いて点検する方法がある（例えば非特許文献１）。この方法では、作業者を派遣するコストが生じると共に常時監視することが不可能である。
【0003】
構造物のクラックを常時検出するためには、例えば非特許文献２に開示されたセンサを構造物に貼り付けて監視する方法がある。センサは、例えば、伝導体－絶縁体－伝導体の三層構造であり、そのインピーダンスの変化でクラックの有無を検出する。
【先行技術文献】
【非特許文献】
【0004】
【文献】コンクリート構造物のクラック検知ツール〔平成30年11月6日検索〕、インターネット（URL: http://www.kurabo.co.jp/chem/kk_crack_sensor.html）
【文献】原田隆郎、他２名、「PVDF フィルムセンサによるひび割れ検知に関する研究」、構造工学論文集Vol.59A(2013 年3 月)pp. 47-55.
【発明の概要】
【発明が解決しようとする課題】
【0005】
しかしながら、センサの形状は、予め決まっているため種々の構造の様々な形状に適用するのが難しいという課題がある。また、センサの端子に水等の誘電体が付着すると、センサのインピーダンスが変化する。この水等の誘電体が付着することで生じるインピーダンスの変化及び経年による変化と、クラックが生じたことによるインピーダンスの変化を見分けるのが困難であるという課題がある。
【0006】
本発明は、この課題に鑑みてなされたものであり、センサの形状に限定されず、クラックによるインピーダンスの変化を容易に見分けることが可能なクラック検出装置とその方法を提供することを目的とする。
【課題を解決するための手段】
【0007】
本発明の一態様に係るクラック検出装置は、構造物のクラックを検出するクラック検出装置であって、前記構造物に貼り付けられる伝導体－絶縁体－伝導体の三層構造のセンサ部と、所定の周波数の範囲を掃引して前記センサ部のインピーダンスが最大及び最小になる複数の周波数を取得する周波数特性取得部と、複数の周波数の不均一性からクラックの有無を判定するクラック有無判定部と、クラックの位置と周波数のずれの方向の関係を記録したクラック位置テーブルと、前記クラック有無判定部でクラックがあると判定された場合に、前記周波数特性取得部で取得した２つの周波数の差を取って符号を求め、該符号で前記クラック位置テーブルを参照してクラックの位置を検出するクラック位置検出部と
を備えることを要旨とする。
を備えることを要旨とする。
【0008】
また、本発明の一態様に係るクラック検出方法は、クラックを検出する対象の構造物に貼り付けられるセンサ部を含むクラック検出装置が実行するクラック検出方法であって、所定の周波数の範囲を掃引して前記センサ部のインピーダンスが最大及び最小になる複数の周波数を取得する周波数特性取得ステップと、複数の周波数の不均一性からクラックの有無を判定するクラック有無判定ステップと、前記クラック有無判定ステップでクラックがあると判定された場合に、前記周波数特性取得ステップで取得した２つの周波数の差を取って符号を求め、該符号でクラックの位置と周波数のずれの方向の関係を記録したクラック位置テーブルを参照してクラックの位置を検出するクラック位置検出ステップとを行うことを要旨とする。
【発明の効果】
【0009】
本発明によれば、センサの形状に限定されず、クラックによるインピーダンスの変化を容易に見分けることが可能なクラック検出装置とその方法を提供することができる。
【図面の簡単な説明】
【0010】
【図1】本発明の実施形態に係るクラック検出装置の機能構成例を示すブロック図である。
【図2】クラックがある場合のセンサ部の平面を模式的に示す図である。
【図3】クラックがある場合のセンサ部を伝送線で表した場合の回路モデルを示す図である。
【図4】センサ部のインピーダンスの周波数特性の例を模式的に示す図である。
【図5】クラックの位置を表すパラメータの例を示す図である。
【発明を実施するための形態】
【0011】
以下、本発明の実施形態について図面を用いて説明する。複数の図面中同一のものに
は同じ参照符号を付し、説明は繰り返さない。
は同じ参照符号を付し、説明は繰り返さない。
【0012】
（センサ部）
本発明の実施形態に用いるセンサ部について説明する。センサ部は、伝導体－絶縁体－伝導体の三層構造を、柔らかいフィルムで挟んだ構造である。
本発明の実施形態に用いるセンサ部について説明する。センサ部は、伝導体－絶縁体－伝導体の三層構造を、柔らかいフィルムで挟んだ構造である。
【0013】
センサ部は、クラックの有無を検出する対象の構造物に貼り付けて使用される。構造物に貼り付けられたセンサ部の構造物側の伝導体を一方の端子、外側の伝導体を他方の端子とみなすと、センサ部は一方の端子対AB及び他方の端子対CDで表される伝送路の等価回路で表現することができる。
【0014】
伝送路の一方の端子対ABからみたインピーダンスZtは以下の式で表せる。
【0015】
【数1】
【0016】
上式でZoは伝送線の特性インピーダンス、γは伝播定数（複素数）、Ltは伝送線の長さを表す。ここで一般的には伝搬定数γに含まれるαを減衰定数（実数）を無視して、位相定数（実数）βとｊ=（－１）0.5を用いて、γ＝ｊβとして式（１）を展開すると以下の式を得る。
【0017】
【数2】
【0018】
伝送線の単位長あたりのインダクタンスをLo、容量をCo、円周率をπ、周波数をω／２πとするとβ＝ω（Lo Co）０．５の関係が分布定数回路の理論で知られている。式（２）において、分母と分子に含まれるβLt＝ω（Lo Co）０．５Ltを変数とする三角関数により、インピーダンスは周期的にピークとなる周波数特性となる。インピーダンスがピークとなる周期が伝送線の長さLtに関係するのは、端子対CDが電磁波の反射点となって端子対ABと端子対CDの間で定在波がたつためである。構造物にクラックが生じた時では、クラックが電磁波の反射点となり定在波のたつ周波数が変わるため、インピーダンスがピークとなる周波数が変化する。
【0019】
検知対象である構造物が屋外に設置された場合、雨滴等によりセンサ部の端に水が付着することがある。また、構造物が屋内にある場合でも結露などによりセンサ部の端に水が付着することがある。この場合のセンサ部は、伝送線の他方の端子対CDに水の付着による寄生容量Cpが接続された回路モデルで表される。この時の端子対ABからみたインピーダンスはγ＝ｊβとして以下の式で表される。
【0020】
【数3】
【0021】
式（３）において分母と分子に含まれる三角関数の変数はβLt＋ωCｐZoとなり、インピーダンスがピークとなる定在波のたつ周波数は寄生容量Cpにも依存する。式（２）で表されるCpがないときの定在波のたつ周波数と比較して、Cpがあるときの定在波のたつ周波数は変化する。
【0022】
また、絶縁体の誘電率の温度・湿度等での変化や経年での変化でも定在波のたつ周波数は変化する。
【0023】
このように、構造物にクラックが生じた場合だけでなく、センサ部の端に水等の誘電体が付着した場合や誘電体の温度・湿度や経年による変化でも定在波のたつ周波数が変化する。このため、センサ部の端に水等の誘電体が付着した場合および誘電体の温度・湿度や経年による変化の場合と、クラックが生じた場合を見分ける必要がある。
【0024】
本発明の実施形態に係るクラック検出装置は、構造物にクラックが生じたことによるセンサ部のインピーダンスの変化とそれ以外によるインピーダンスの変化の見分けを容易にし、センサ部の形状による影響を受けない構成である。
【0025】
（クラック検出装置）
図１は、本発明の実施形態に係るクラック検出装置の機能構成例を示すブロック図である。図１に示すクラック検出装置１００は、センサ部１０、周波数特性取得部２０、クラック有無判定部３０、クラック位置テーブル４０、及びクラック位置検出部５０を備える。
図１は、本発明の実施形態に係るクラック検出装置の機能構成例を示すブロック図である。図１に示すクラック検出装置１００は、センサ部１０、周波数特性取得部２０、クラック有無判定部３０、クラック位置テーブル４０、及びクラック位置検出部５０を備える。
【0026】
センサ部１０は、構造物２００に貼り付けられる伝導体－絶縁体－伝導体の三層構造である。センサ部１０は、例えば、ＰＶＤＦフィルムで構成することができる。
【0027】
周波数特性取得部２０は、所定の周波数の範囲を掃引してセンサ部１０のインピーダンスが最大及び最小になる複数の周波数を取得する。周波数特性取得部２０は、例えば、インピーダンスアナライザ等の一般的な計測器で構成することもできる。又は、ＦＦＴアナライザで構成してもよい。ＦＦＴアナライザでセンサ部１０の周波数特性を測定し、複数の直列共振周波数と並列共振周波数を測定するようにしてもよい。具体例については後述する。
【0028】
クラック有無判定部３０は、周波数特性取得部２０で取得した複数の周波数の不均一性からクラックの有無を判定する。センサ部１０にクラックが無い場合は、周波数特性取得部２０で取得した複数の周波数のそれぞれの間隔は一定である。センサ部１０にクラックが生じた場合は、複数の周波数のそれぞれの間隔にバラツキが生じる。よって、複数の周波数の不均一性からクラックの有無を判定することができる。また、複数の周波数の不均一性からクラックの有無を判定するので、センサ部１０の大きさが限定されない。つまり、任意の大きさのセンサ部１０を用いることができる。詳しくは後述する。
【0029】
クラック位置テーブル４０は、センサ部１０上のクラックの位置と周波数のずれの方向の関係を記録したテーブルである。クラック位置テーブル４０は、予め設けられる。具体例は後述する。
【0030】
クラック位置検出部５０は、クラック有無判定部３０でクラックがあると判定された場合に、周波数特性取得部２０で取得した２つの周波数の差を取って符号を求め、該符号でクラック位置テーブル４０を参照してクラックの位置を検出する。クラックが生じた位置によって、センサ部１０のインピーダンス特性は変化するので、２つの周波数の差を取り、予め設定された基準と比較することでクラックの位置を特定することができる。詳しくは後述する。
【0031】
このように本実施形態に係るクラック検出装置１００は、構造物２００のクラックを検出するクラック検出装置であって、構造物２００に貼り付けられる伝導体－絶縁体－伝導体の三層構造のセンサ部１０と、所定の周波数の範囲を掃引してセンサ部１０のインピーダンスが最大及び最小になる複数の周波数を取得する周波数特性取得部２０と、複数の周波数の不均一性からクラックの有無を判定するクラック有無判定部３０と、クラックの位置と周波数のずれの方向の関係を記録したクラック位置テーブル４０と、クラック有無判定部３０でクラックがあると判定された場合に、周波数特性取得部２０で取得した２つの周波数の差を取って符号を求め、該符号でクラック位置テーブル４０を参照してクラックの位置を検出するクラック位置検出部５０とを備える。
【0032】
これにより、センサの形状に限定されず、クラックによるインピーダンスの変化を容易に見分けることが可能になる。
【0033】
以下でクラックを検知する動作の原理を詳しく説明する。図２にクラックがある場合のセンサ部１０の模式図を示し、図３にクラックがある場合の回路モデルを示す。図２の模式図では上部の伝導体の左端からLt1で右端からLt2の箇所に間隙がδlのクラックが生じている状態を示している。図３の回路モデルにおいてCcはクラックの間隙の部分に存在する浮遊容量を表し、L1とC1図２においてつながっている部分を伝送線とみなした場合の単位長さあたりのインダクタンスと容量を表している。長さLt1とLt2の伝送線の特性インピーダンスをZo（実数）、伝播定数をγ（複素数）としている。
【0034】
端子AとB間の電位差をVab、端子２と２‘間の電位差をV2、端子３と３’間の電位差をV3、端子CとD間の電位差をVcdとする。また、端子Aにおいて伝送線に沿って部分に流れる電流をIab、端子２において伝送線部分の電流をI2、端子３において伝送線部分の電流をI3、端子Cにおいて伝送線に沿って流れる電流をIcdとする。分布定数回路の理論に従い、端子Bには逆方向にIabが流れ、同様に端子２‘には逆方向にI2、端子３’には逆方向にI3、端子Dには逆方向にIcdが流れるとする。本回路モデルでは以下の回路方程式が成立する。
【0035】
【数4】
【0036】
ここで、K1とK2は長さLt1の伝送線の境界条件によって決定される定数であり、K3とK4は長さLt2の伝送線の境界条件によって決定される定数である。また上式でeは自然対数の底（ネイピア数）を表す。
【0037】
端子CDが開放であることからIcd=0となり、式（１３）からK3とK4の関係は以下の式となる。
【0038】
【数5】
【0039】
式（１４）を式（１０）に代入すると以下の式が得られる。
【0040】
【数6】
【0041】
同様に式（１４）を式（１１）に代入すると以下の式が得られる。
【0042】
【数7】
【0043】
式（１６）で式（１５）を割るとV3とI3の関係式が以下のようになる。
【0044】
【数8】
【0045】
式（９）からV2はV3とI3を用いて以下の式で表される。
【0046】
【数9】
【0047】
式（８）からI2はV3とI3を用いて以下の式で表される。
【0048】
【数10】
【0049】
式（１８）に代入する。
【0050】
【数11】
【0051】
ここで式（２０）から式（２２）のパラメータを用いると式（１８）と式（１９）はそれぞれ式（２３）と式（２４）のように表される。
【0052】
【数12】
【0053】
式（６）と式（７）からK1とK2は式（２５）と式（２６）で表される。
【0054】
【数13】
【0055】
ここでLp=L1+L2、Ln=-L1+L2すると、L1=(Lp-Ln)/2、L2=(Lp+Ln)/2となり、式（２５）と式（２６）は以下の式で表される。
【0056】
【数14】
【0057】
式（１７）、式（２７）、式（２８）を式（４）に代入するとVabに関して以下の式が得られる。
【0058】
【数15】
【0059】
ここで、以下のパラメータを導入すると、式（２９）は式（３２）のように展開できる。
【0060】
【数16】
【0061】
式（１７）、式（２７）、式（２６）を式（５）に代入するとＩabに関して以下の式が得られる。
【0062】
【数17】
【0063】
ここで式（３０）と式（３１）を用いて展開するとIabの表式は以下のようになる。
【0064】
【数18】
【0065】
端子ＡＢからみたインピーダンスZab=Vab/Iabは式（３２）と式（３３）から以下の式で表される。
【0066】
【数19】
【0067】
式（２）と式（３）の式展開と場合と同様に、βを位相定数（実数）、ｊを（－１）0.5として、γ＝ｊβとすると、Zabの分子は以下のように展開できる。
【0068】
【数20】
【0069】
Zabの分母は以下のように展開できる。
【0070】
【数21】
【0071】
クラックが入った場合では、伝送線のインピーダンスZabの式の分母と分子に含まれる三角関数の変数は、{β(Lt1+Lt2)+θz}とβ(Lt1-Lt2)の２種となる。クラックが入っておらず伝送線の端の付着物（寄生容量Cp）がない場合ではインピーダンスの式に含まれる三角関数の変数は式（A-2）よりβLtの１種であり、クラックが入っておらず伝送線の端の付着物（寄生容量Cp）がある場合のインピーダンスの式に含まれる三角関数の変数は式（A-3）より(βLt＋ωCｐZo)の１種である。この考察からクラックの有無によりインピーダンスの式に含まれる三角関数の変数の数が変化することが分かる。ここでLtは伝送線の長さを表し、Lt=Lt1+Lt2+δlの関係がある。
【0072】
クラックが無い場合でのインピーダンスのピークの周期は、β＝ω（Lo Co）０．５の関係とインピーダンスの式に含まれる三角関数の変数から考えて、Cpが無い場合で（Lo Co）０．５Lt、Cpがある場合で（Lo Co）０．５Lt＋CｐZoと周波数に依存せず一定である。
【0073】
これに対し、クラックがある場合において、２種の変数のうちβ(Lt1-Lt2)の成分に関しては、インピーダンスがピークとなる周期が（Lo Co）０．５(Lt1-Lt2)となる。{β(Lt1+Lt2)+θz}の成分に関しては、以下のように考察できる。
【0074】
以下の式で表されるパラメータθzはLp=L1+L2、式（２０）、式（２１）を用いて式（３８）のように表される。
【0075】
【数22】
【0076】
構造物に生じたクラックがまだ初期段階でクラック間隙δlが非常に短い場合には以下の近似が成立する。
【0077】
【数23】
【0078】
さらに周波数ω/2πが低い領域では式（３９）は以下の式（４０）となり、θzは周波数に正比例する。
【0079】
【数24】
【0080】
従って、構造物に生じたクラックがまだ初期段階でクラック間隙δlが非常に短く周波数が低い領域における{β(Lt1+Lt2)+θz}の成分の周期は{（Lo Co）０．５(Lt1+Lt2)+(L1+C1 Zo2)δl/(2Zo) }となる。
【0081】
また、クラックが初期段階でのLp、Ln、|Az|と以下の式で与えられる。
【0082】
【数25】
【0083】
式（４２）、式（４３）よりクラックがまだ初期段階で間隙δlが非常に短く周波数が低い領域ではLnと|Az|は周波数に対し一定となる。クラックが生じている個所の伝導体の幅は、クラックのない伝導体の幅と比較して狭いためLo<L1、Co>C1となる。式（４２）よりLnは負となる。
【0084】
次に、クラックがまだ初期段階で間隙δlが非常に短く周波数が低い領域でアドミタンスとインピーダンスがピークとなる周波数について考察する。アドミタンスはインピーダンスの逆数のため、アドミタンスのピークとなる周波数はZabの分子である式（３５）がゼロとなる周波数である。式（４１）から式（４３）を代入すると以下の式を得る。
【0085】
【数26】
【0086】
クラックがまだ初期段階で間隙δlが非常に短く周波数が低い領域ではZo>>ω(L1-LoC1/Co)δl/2となり式（４４）の第１項が支配的となるが、Zabの分子がゼロ近傍を議論することから式（４４）では第２項も無視せずに考慮することとした。
【0087】
支配的な式（４４）の第1項がゼロとなる周波数の近傍に式（４４）がゼロとなる周波数が存在すると考えられる。式（４４）の第1項のcos関数の大きさがゼロとなるとき第1項がゼロとなる。円周率をπとすると、cos関数の変数がπ/2の奇数倍となるときcos関数の大きさはゼロとなる。cos関数の性質によりｃos関数の大きさがゼロになる近傍では式（４４）の第１項のcos関数は以下の式で近似できる。
【0088】
【数27】
【0089】
式（４５）と式（４６）とβ=ω(LoCo)0.5を用いて、式（４４）がゼロとなるすなわちアドミタンスのピークとなる周波数ω/2πは以下の条件を満たす。
【0090】
【数28】
【0091】
式（４７）と式（４８）をまとめると以下の式となる。
【0092】
【数29】
【0093】
インピーダンスのピークとなる周波数はZabの分母である式（３６）がゼロとなる周波数である。式（４１）から式（４３）を代入すると以下の式を得る。
【0094】
【数30】
【0095】
クラックがまだ初期段階で間隙δlが非常に短く周波数が低い領域ではZo>>ω(L1-LoC1/Co)δl/2となり式（５０）の第１項が支配的となるが、Zabの分母がゼロ近傍を議論することから式（５０）では第２項も無視せずに考慮することとした。
【0096】
支配的な式（５０）の第1項がゼロとなる周波数の近傍に式（５０）がゼロとなる周波数が存在すると考えられる。式（５０）の第1項のsin関数の大きさがゼロとなるとき第1項がゼロとなる。Sin関数の変数がπの整数倍となるときsin関数の大きさはゼロとなる。sin関数の性質によりsin関数の大きさがゼロになる近傍では式（５０）の第１項のsin関数は以下の式で近似できる。
【0097】
【数31】
【0098】
式（５１）と式（５２）とβ=ω(LoCo)0.5を用いて、式（４４）がゼロとなるすなわちアドミタンスのピークとなる周波数ω/2πは以下の条件を満たす。
【0099】
【数32】
【0100】
式（５３）と式（５４）をまとめると以下の式となる。
【0101】
【数33】
【0102】
さらに式（４９）と式（５５）をまとめると以下の式を得る。
【0103】
【数34】
【0104】
式（５６）において周波数fm=ωm/2πを用いた。
【0105】
図４にインピーダンスの周波数特性の例を示す。図４の周波数f1やf3のインピーダンスの極小値はアドミタンスのピークを表す。アドミタンスとインピーダンスのピークの周波数を直流側から番号を付与すると式（５６）中のmと同じになる。
【0106】
クラックの無いδl=0の場合では式（２）や式（３）について議論したとおりアドミタンスとインピーダンスのピークの周波数は周期的であるため、隣り合うアドミタンスとインピーダンスのピークの周波数の差は一定である。このため、アドミタンスとインピーダンスのピークの周波数を直流側から付与した番号mでその周波数を規格化した周波数は一定となる。式（５６）においてもδl=0とすると右辺に周波数に依存する右辺の第３項がゼロとなるため、mで規格化した周波数fm/mは一定となる。
【0107】
クラックのある場合では、fm/mは一定ではなくなる。このため、アドミタンスとインピーダンスのピークとなる周波数を数点取得し、mで規格化した周波数の一致・不一致でクラックの有無が判定できる。
【0108】
次にクラックの位置について考察する。クラックの位置をLt1とLt2のいずれかで表すことができるが、以下の式で表される全体の長さLt1+Lt2で規格化したパラメータを導入しても表すことができる。
【0109】
【数35】
【0110】
図５にパラメータLocとトランスデューサ上の位置の関係を示す。アドミタンスとインピーダンスのピークとなる周波数を計測中に全体の長さLt1+Lt2は変化しないため式（５７）の分母は一定である。トランスデューサの中心ではLt1=Lt2となりLoc=0となる。LocはLt1と正比例するためクラックの位置が図５上の右側になる（Lt1が大きくなる）とLocは0より大きくなる。また、LocはLt2に-1を乗じた数に比例するためクラックの位置が図５上の左側になる（Lt2が大きくなる）とLocは0より小さくなる。
【0111】
式（５６）の右辺第３項のcos関数の変数にLt1-Lt2があるため、第３項のみを抽出すればパラメータLocとの関係式が得られると考えられる。クラックがあるためfm/mは一定ではないが、f1の近傍にあると考え、さらにf1はクラックが無い時にアドミタンスがピークになる最低の周波数fbの近傍にあると考えて以下の式でfmを表す。
【0112】
【数36】
【0113】
fbは1/4波長で定在波がたったときの周波数のため、トランスデューサの全体の長さLt1+Lt2と電磁波の速度(LoCo)-0.5とに以下の関係が成立する。
【0114】
【数37】
【0115】
式（５８）を用いて式（５６）右辺の第３項のcos関数を展開すると以下の式を得る。
【0116】
【数38】
【0117】
式（６０）の展開では2πδf(LoCo)0.5 (Lt1-Lt2)<<1とした。式（６０）に式（５９）を代入すると以下の式を得る。
【0118】
【数39】
【0119】
式（６１）において第１項がゼロより十分大きくδf/fb<<1とすれば式（５６）右辺の第３項のcos関数を以下の式で近似できる。
【0120】
【数40】
【0121】
式（６２）を用いて式（５６）は以下の式となる。
【0122】
【数41】
【0123】
任意のma番目の周波数とmb番目の周波数についてmで規格化した周波数の逆数の差をとるとLocに依存する第３項を取り出すことができ以下の式を得る。
【0124】
【数42】
【0125】
式（６４）の右辺は周波数がなくても導出可能であり、maとmbは予め決定できる数である。また、任意の長さLt1+Lt2のトランスデューサに対してLocの範囲は-1<Loc<1であるため、BLを除けばどのような長さのトランスデューサ対しても式（６４）の右辺は予め導出できる。前述の式（４２）で説明したとおりBL>0であるため、右辺の正負の符号は左辺と同じである。式（６４）の左辺のfmaとfmbは、アドミタンスとインピーダンスの周波数依存性の測定により得られる。このため、予めLocとma、mbの組合せに対する式（６５）で表される第１の評価関数の結果をクラック位置テーブルとして保持しておき、アドミタンスとインピーダンスの周波数依存性の測定結果から得られた式（６６）で表される第２の評価関数Evmの符号、ma、mbの組合せとクラック位置テーブル４０からクラックの位置であるLocの領域を逆引きすることができる。
【0126】
【数43】
【0127】
式（６５）と式（６６）でsgnは符号関数である。
【0128】
アドミタンスとインピーダンスの個数を4とした場合のクラック位置テーブル４０の例を図６（ａ）に示す。例えば、はじめにma=3とmb=1でのEvmの符号を評価しEvmが正の時では、次にma=4とmb=3でのEvmを評価する。ma=4とmb=3でのEvmが正の場合では図６（ｂ）に示すトランスデューサ上の（１）の領域にクラックがあることとなる。
【0129】
上述のように、予めLocとma、mbの組合せに対する式（６４）で表される評価関数の結果をクラック位置テーブル４０として保持しておき、アドミタンスとインピーダンスの周波数依存性の測定結果から得られた式（６５）で表されるEvmの符号、ma、mbの組合せとクラック位置テーブル４０からクラックの位置であるLocの領域を逆引きすることができる。
【0130】
周波数特性取得部２０の具体例について説明する。
【0131】
周波数特性取得部２０は、所定の周波数の範囲を掃引してセンサ部１０のインピーダンスが最大及び最小になる複数の周波数を取得する。図７に周波数特性取得部２０の周波数可変信号源２１、負荷部２２とセンサ部１０の回路モデルを示す。周波数可変信号源２１から出力された電圧信号は負荷部２２（抵抗Rs）を介してセンサ部１０に印加される。センサ部１０のインピーダンスが大きい時では端子A-B間の電圧Vaも大きいため、Vaがピークとなる周波数はインピーダンスのピークとなる周波数と一致する。センサ部１０のアドミタンスが大きい時では電流Isが大きくなる。電流Isが大きくなると負荷部２２での電圧降下が大きくなるため、Vs-Vaを測定すればアドミタンスの大きさがわかる。従ってVs-Vaがピークとなる周波数とアドミタンスがピークとなる周波数は一致する。
【0132】
図８に周波数特性取得部２０の構成例を示す。端子Ａからの信号は直交復調部２４と差動演算部２３に入力される。直交復調部２４では周波数可変信号源２１から入力される基準信号の同相成分と直交成分の振幅に端子Ａからの信号を分離し周波数特性算出部２６に出力する。差動演算部２３では負荷部２２からの信号と端子Ａからの信号の差をとって直交復調部２５に出力する。直交復調部２５では基準信号の同相成分と直交成分の振幅に差動演算部２３からの入力信号を分離し周波数特性算出部２６に出力する。周波数特性算出部２６では直交復調部２４の出力信号からインピーダンスのデータを生成し、直交復調部２５の出力信号からアドミタンスのデータを生成して周波数特性データとしてクラック有無判定部３０に送信する。
【0133】
所定の周波数範囲の周波数に対応した周波数特性データ取得したのち、アドミタンスとインピーダンスがピークとなる周波数を抽出するとともに低周波側から周波数に対応した番号を付与する。
【0134】
この後、アドミタンスとインピーダンスのピークの周波数を番号で規格化した複数の周波数が所定の不均一さ未満であればクラックが無いと判定する。規格化した複数の周波数が所定の不均一さ以上である場合には、２種の番号とピークの周波数を用いて式（６６）で表される第２の評価関数で評価した結果とクラック位置テーブル４０を照合してクラックの位置を決定する。クラックの有無の判定結果とクラックの位置の情報は、通信部（図示せず）を介して外部へ送信するようにしても良い。
【0135】
以上説明したように本実施形態に係るクラック検出装置１００は、構造物２００のクラックを検出するクラック検出装置であって、構造物２００に貼り付けられる伝導体－絶縁体－伝導体の三層構造のセンサ部１０と、所定の周波数の範囲を掃引してセンサ部１０のインピーダンスが最大及び最小になる複数の周波数を取得する周波数特性取得部２０と、複数の周波数の不均一性からクラックの有無を判定するクラック有無判定部３０と、クラックの位置と周波数のずれの方向の関係を記録したクラック位置テーブル４０と、クラック有無判定部３０でクラックがあると判定された場合に、周波数特性取得部２０で取得した２つの周波数の差を取って符号を求め、該符号で前記クラック位置テーブル４０を参照してクラックの位置を検出するクラック位置検出部５０とを備える。
【0136】
また、クラック有無判定部３０は、インピーダンスが最大及び最小になる複数の周波数を低い周波数から１から順番に付番し、該付番した番号でそれぞれの周波数を正規化した正規化周波数を求め、複数の正規化周波数のそれぞれが同じであればクラックが無いと判定し、複数の正規化周波数が得られればクラックありと判定する。この判定方法によれば、周波数を番号で規格化した複数の周波数が所定の不均一さでクラックの有無を判定するので、ゼンサ部１０のセンサの形状の影響を受けない。つまり、センサはどのような形状で有っても構わない。
【0137】
このように本実施形態に係るクラック検出装置１００は、クラックによるインピーダンスの変化を容易に見分けることを可能にする。なお、本発明は、上記の実施形態に限定されるものではなく、特許請求の範囲に記載された要旨の範囲内で変形が可能である。
【符号の説明】
【0138】
１０：センサ部
２０：周波数特性取得部
２１：周波数可変信号源
２２：負荷部
２３：差動演算部
２４，２５：直交復調部
２６：周波数特性算出部
３０：クラック有無判定部
４０：クラック位置テーブル
５０：クラック位置検出部
２０：周波数特性取得部
２１：周波数可変信号源
２２：負荷部
２３：差動演算部
２４，２５：直交復調部
２６：周波数特性算出部
３０：クラック有無判定部
４０：クラック位置テーブル
５０：クラック位置検出部
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】