(57)【要約】【課題】走査モアレ法やサンプリングモアレ法は静的な状態での変位ひずみ分布計測では多くの利点を有するが、動的な変位ひずみ分布計測をする場合、走査型レーザ顕微鏡を使用すると高分解能な画像を得るのに時間がかかり、あるいは面内変形の測定では試料ステージまたは顕微鏡の走査線を９０°回転する必要があり、その間の試料の変形の影響を受けると言う問題があって、迅速かつ正確に構造材料の二次元の変位ひずみ分布を測定できなかった。【解決手段】走査線数を少なくして、変位とひずみの測定感度を低下させることなく、より高速にモアレ画像を取得した。また交差格子を試料格子に用いて、ｘとｙの２方向の走査点間隔を適宜配置して、試料の回転をしないで取得した一枚のモアレ縞画像から両方向の一次元走査モアレ縞に分離し、高速モアレ法と組み合わせて迅速かつ正確に構造材料の二次元の変位ひずみ分布を測定できるようになった。【選択図】図１
2017146202.tif 000002
【特許請求の範囲】
【請求項１】
縦横の走査点からなり基本走査線間隔Ｔ0の整数倍の間隔Ｔ(整数は１以上の自然数)に走査し得る走査線を有する走査型顕微鏡を備え、
試料台に載置されその表面に一方向にピッチｐ(ｐはＴ0に近接しているとする)の規則格子が生成された試料の表面を前記走査型顕微鏡の走査線により走査して前記試料表面のモアレ縞画像を取得し、
当該取得したモアレ縞画像の走査点の輝度情報から前記試料の変位ひずみ量を計算し、
前記計算した試料の変位ひずみ量を表示する、
走査モアレ方法による変位ひずみ分布測定システムであって、
前記走査型顕微鏡の走査点間隔Ｔを前記ピッチｐの２以上の整数倍またはその近傍とし前記規則格子と前記走査点が小さな不一致またはずれが生じるようにしたことを特徴とする高速走査モアレ方法による高速変位ひずみ分布測定システム。
【請求項２】
前記規則格子ピッチｐの相対的ひずみ量ε(ｎ)は前記一方向の変位ｕ(ｎ)の一次微分から得たことを特徴とする請求項１に記載の高速変位ひずみ分布測定システム。
ただし、前記一方向をｙ方向とし、ｍは前記モアレ縞の次数、ｕ(ｎ)は前記基本走査線に関してｙ方向におけるモアレ縞の変位とする。
【数１７】
2017146202.tif 000019
【数１８】
2017146202.tif 000020
【請求項３】
前記ひずみ量εactualを変形前後の前記規則格子の幾何学的関係から得たことを特徴とする請求項２に記載の高速変位ひずみ分布測定システム。
ただしｐ0はｙ方向における変形前の試料の規則格子ピッチとする。
【数１９】
2017146202.tif 000021
【請求項４】
縦横の走査点からなり基本走査線間隔Ｔ0の整数倍の間隔Ｔ(整数は１以上の自然数)に走査し得る走査線を有する走査型顕微鏡を備え、
試料台に載置されその表面に直交する二方向にピッチｐ(ｐはＴ0に近接しているとする)の規則交差格子が生成された試料の表面を前記走査型顕微鏡の走査線により走査して前記試料表面の二次元モアレ縞画像を取得し、
当該取得した二次元モアレ縞画像の走査点の輝度情報から前記試料変位ひずみ量を計算し、
前記計算した試料の変位ひずみを計測して表示する、
走査モアレ方法による変位ひずみ分布測定システムであって、
前記規則交差格子と前記基本間隔Ｔ0の走査点とを近接して配置し前記規則交差格子と前記走査点が小さな不一致またはずれが生じるようにしたことを特徴とする二次元走査モアレ方法による変位ひずみ分布測定システム。
【請求項５】
前記規則交差格子の相対的ひずみ量εi(ｉ＝ｘ，ｙ)は前記二次元モアレ縞画像に画像処理(例えば、フーリエ変換やローパスフィルタ等)を施して直交する二方向の一次元走査モアレ縞に分離した後に前記規則交差格子の二方向における基本走査線に相対的な変位ｕi(ｉ＝ｘ，ｙ)の一次微分から求めた事を特徴とする請求項４に記載の変位ひずみ分布測定システム。
ただし、前記二方向はｘとｙ方向とする。
【数２０】
2017146202.tif 000022
【数２１】
2017146202.tif 000023
【請求項６】
前記二方向の規則交差格子の前記ひずみ量は、変形前後の前記規則交差格子ピッチの幾何学的関係(関数)から求めたことを特徴とする請求項５に記載の変位ひずみ分布測定システム。
ただし変形前のｘ方向またはｙ方向の試料格子ピッチをｐi_0(ｉ＝ｘ，ｙ)とする。
【数２２】
2017146202.tif 000024
【請求項７】
縦横に基本走査線間隔Ｔ0の整数倍の間隔Ｔ(整数は１以上の自然数)に走査し得る走査点からなる走査線を有する走査型顕微鏡を備え、
試料台に載置されその表面に直交する二方向にピッチｐ(ｐはＴ0に近接しているとする)の規則交差格子が生成された試料の表面を前記走査型顕微鏡の走査線により走査して前記試料表面の二次元モアレ縞画像を取得し、
当該取得した二次元モアレ縞画像の走査点の輝度情報から前記試料変位ひずみ量を計算し、
前記計算した試料の変位ひずみ量を表示する、
走査モアレ方法による変位ひずみ分布測定システムであって、
前記走査型顕微鏡の直交する二方向の走査線間隔Ｔx、Ｔyを前記ピッチｐの整数倍(整数は１以上の自然数)またはその近傍となるようにし、かつ、
少なくとも前記二方向の走査線間隔の一つは前記ピッチｐの２以上の整数倍またはその近傍とし、
前記規則交差格子と前記走査線間隔Ｔx、Ｔyの走査点とを近接して配置し前記規則交差格子と前記走査点が小さな不一致またはずれが生じるようにしたことを特徴とする高速二次元走査モアレ方法による高速変位ひずみ分布測定システム。
【請求項８】
前記規則交差格子の相対的ひずみ量εi(ｎi)(ｉ＝ｘ，ｙ)は前記二次元モアレ縞画像に画像処理(例えば、フーリエ変換やローパスフィルタ等)を施して直交する二方向の一次元走査モアレ縞に分離した後に前記規則交差格子の二方向における走査点に相対的な変位ｕi(ｎi)(ｉ＝ｘ，ｙ)の一次微分から求めた事を特徴とする請求項７に記載の変位ひずみ分布システム。
ただし、前記二方向はｘとｙ方向とし、ｍiはｉ方向(ｉ＝ｘ，ｙ)における前記モアレ縞の次数、ｕi(ｎi)は前記基本走査線に関してｉ方向におけるモアレ縞の変位とする。
【数２３】
2017146202.tif 000025
【数２４】
2017146202.tif 000026
【請求項９】
前記二方向の規則交差格子の前記ひずみ量は、変形前後の前記規則交差格子ピッチの幾何学的関係(関数)から求めたことを特徴とする請求項８に記載の変位ひずみ分布測定システム。
ただし変形前のｘ方向またはｙ方向の試料格子ピッチをｐi_0(ｉ＝ｘ，ｙ)とする。
【数２５】
2017146202.tif 000027
【請求項１０】
少なくとも対物レンズと二方向に整列した撮像素子からなるＣＣＤ・ＣＭＯＳイメージセンサを備えた撮影装置を用いてその表面に一方向にピッチＰの規則格子が生成された測定物の表面を撮影してモアレ縞画像を取得し、
当該取得したモアレ縞画像の撮像素子に記憶された輝度情報から前記測定物の変位ひずみ量を計算し、
前記計算した測定物の変位ひずみ量を表示する、
サンプリングモアレ方法による変位ひずみ分布測定システムであって、
前記一方向の前記ピッチＰの前記イメージセンサにおけるピッチｐの整数倍(整数は１以上の自然数)またはその近傍にある前記一方向の撮像素子列に前記モアレ縞画像を記憶し、
前記記憶された輝度情報から前記測定物の変位ひずみ量を計算したことを特徴とする高速サンプリングモアレ方法による高速変位ひずみ分布測定システム。
【請求項１１】
少なくとも対物レンズと二方向に整列した撮像素子からなるＣＣＤ・ＣＭＯＳイメージセンサを備えた撮影装置を用いてその表面に前記二方向にピッチＰの規則交差格子が生成された測定物の表面を撮影して二次元モアレ縞画像を取得し、
当該取得した二次元モアレ縞画像の撮像素子に記憶された輝度情報から前記測定物の変位ひずみ量を計算し、
前記計算した測定物の変位ひずみ量を表示する、
サンプリングモアレ方法による変位ひずみ分布測定システムであって、
前記二方向の前記ピッチＰの前記イメージセンサにおけるピッチｐの整数倍(整数は１以上の自然数)またはその近傍にある撮像素子列であって、かつ、
少なくとも前記一方向の撮像素子列は前記ピッチｐの２以上の整数倍またはその近傍にある、
前記二方向の撮像素子列に前記モアレ縞画像を記録し、
前記記録された輝度情報から前記測定物の変位ひずみ量を計算したことを特徴とする高速二次元サンプリングモアレ方法による高速変位ひずみ分布測定システム。
【請求項１２】
前記二方向は直交していることを特徴とする請求項１０又は請求項１１のいずれか１項に記載する高速変位ひずみ分布測定システム。
【請求項１３】
請求項１乃至請求項１２のいずれか１項に記載する前記取得した前記モアレ縞画像または前記二次元モアレ縞画像を当該処理する事を特徴とするプログラムおよび当該プログラムを収納した記録媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、モアレ法による高速かつ高精度で構造材料の二次元変位ひずみ分布を測定に関し、主に実験力学や非破壊検査の光学測定技術の分野に属する。
【背景技術】
【０００２】
構造材料の変位ひずみ分布を測定する技術は、ナノマイクロ材料の機械的特性評価から大型インフラ構造物の健全性評価まで幅広く利用されている。
【０００３】
これまでに走査モアレ縞を利用したマイクロ/ナノスケールでの変位ひずみ計測では、レーザ走査顕微鏡(Laser Scanning Microscope: ＬＳＭ)(非特許文献１)、走査電子顕微鏡(Scanning Electron Microscope: ＳＥＭ)(特許文献１)、原子間力顕微鏡(Atomic Force Microscope: ＡＦＭ)、透過型電子顕微鏡(Transmission Electron Microscope: ＴＥＭ)などが用いられてきた。
【０００４】
またミリ/メートリスケールでの変位ひずみ計測では、電荷結合素子(Charge-couple Device: ＣＣＤ)、相補性金属酸化膜半導体(Complementary Metal Oxide Semiconductor: ＣＭＯＳ)などのイメージセンサを利用したサンプリングモアレ法が用いられている(特許文献２)。
【０００５】
これらの手法では、構造材料表面に自然に存在するまたは人工的に付けた繰り返し模様を、各種走査型顕微鏡(例えば、ＬＳＭ、ＳＥＭ、ＡＦＭ、ＴＥＭ等)を用いた走査モアレ法では模様ピッチと近い間隔で走査(スキャニング)、または各種イメージセンサ(例えば、ＣＣＤ、ＣＭＯＳ等)を用いたサンプリングモアレ法では画像処理で撮影画像上での模様ピッチと近い画素間隔で間引き処理(ダウンサンウリング)して得られるモアレ縞の変化から定量的に構造材料表面の変位ひずみ分布を評価できる。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特許第５１５２６１５号公報(微小な規則的パターンの配列の乱れ検出測定方法)
【特許文献２】特許第４８３１７０３号公報(物体の変位測定方法)
【非特許文献】
【０００７】
【非特許文献１】Q.H. Wang, T. Hiroshi, S. Kishimoto, Y. Tanaka, and Y. Kagawa, Moire techniques based on memory function of laser scanning microscope for deformation measurement at micron/submicron scales, International Journal of Automation Technology, Vol.9, No.5, pp.494-501(2015).
【発明の概要】
【発明が解決しようとする課題】
【０００８】
従来の走査モアレ法やサンプリングモアレ法は静的な状態での変位ひずみ分布計測では多くの利点を有するが、動的な変位ひずみ分布計測の際に以下のようにいくつかの問題点が挙げられる。
【０００９】
１．走査モアレ法では、レーザ光源を走査して画像を取得するため、高分解能な画像を得るのに時間がかかる。
モアレ画像取得に時間がかかると、変位ひずみ分布評価に必要なモアレパターンが試料クリープの影響を受けやすく、測定誤差が発生しやすい。
【００１０】
２．明確なモアレ縞が生成するためには、走査間隔に近いピッチがある試料格子のみを使用することができるので、使用可能な試料格子の範囲が狭い。
特に、一般的にレーザ走査型顕微鏡(ＬＳＭ)装置で用いられる対物レンズの倍率は通常５×、１０×、２０×、５０×、１００×、１５０×といった具合に不連続であり、ＬＳＭモアレ法を様々なアプリケーションへ応用する際の障害となる。
【００１１】
３．材料性能評価のために不可欠であり、ｘ方向及びｙ方向の両方の面内変形の測定には、従来の走査モアレ法で、試料ステージまたは顕微鏡の走査線を９０°回転する必要がある。
回転前後のモアレ縞を得るために、顕微鏡は二回で走査する必要がある。
回転プロセスや二回走査が不便であり、試料の位置誤差をもたらすことになる。
この場合、時間依存現象を捉えるための動的変形測定を適用することは困難である。
【００１２】
上述の問題を踏まえて、迅速かつ正確に構造材料の二次元の変位ひずみ分布を測定する方法が求められる。
また、実用上では使用可能な試料格子ピッチの範囲を拡大できる手法が望まれている。
【課題を解決するための手段】
【００１３】
図１に、公知のモアレ法と、そのモアレ法を応用した、迅速かつ正確に構造材料の二次元変位ひずみ分布を実現するための本発明の三つの新変形測定方法を提案する。
第１の方法は、図１(ｂ)右上に示す高速走査モアレ法であり、第２の方法は図１(ｃ)左下に示す二次元走査モアレ法であり、第三の方法は図１(ｄ)右下に示した高速二次元走査モアレ法である。
なお、ここでは、走査モアレは回転ミラーなどで機械的にレーザ光源等を走査して得られる各走査点での輝度情報を取得するのに対して、サンプリングモアレ法はＣＣＤ面に配置された二次元画素配列を飛び飛び(間引き)の画素での明るさ情報をサンプリング点として取り出しているという点では異なるものの、離散的にサンプリング点での輝度情報を取得している点は本質的に同じである。
よって本説明で走査モアレ法と言う時には、図面を含めて、走査モアレ法とサンプリングモアレ法に共通な特徴が含まれるものとする。
【００１４】
これらの手法によって、変位ひずみ分布を評価するうえで必要となるモアレ画像の取得速度の向上と二次元の変形測定精度の向上が同時に達成することができる。
【００１５】
１．高速走査モアレ法
高速縞走査とは従前の変形計測の縞走査方法に要する時間よりも走査時間の短縮を意味する。
より詳しくは従前の半分以下の走査線数で同等の精度の変形計測ができる。
【００１６】
最初に、試料表面上に規則格子を作製する。
その後、試料台上に設置した試料をレーザ走査型顕微鏡(ＬＳＭ)や走査型電子顕微鏡(ＳＥＭ)または原子間力顕微鏡(ＡＦＭ)などの走査型顕微鏡より観察する。
【００１７】
次に、走査線の間隔を、試料格子ピッチの整数(≧２)倍に近くなるように調整する。
走査線(試料格子ピッチの整数倍)と試料格子間の干渉から、高速走査モアレパターンが生成することができる(図２)。
この場合、従来の走査モアレ法に比べて走査線数が少ないため、より高速でモアレ画像を取得できる。
【００１８】
最後に、高速走査モアレ法より得られるモアレ縞画像から、全視野で試料の変位とひずみ分布を測定することができる。
【００１９】
１．１  高速走査モアレ縞の形成原理
初めに従前の走査モアレ縞の形成原理を紹介する。
従前のモアレ縞は試料の格子と走査線(参照格子)間の干渉から発生する。
試料格子のピッチと走査線の間隔がそれぞれｙ方向にｐとＴ0として定義されている場合、試料格子と走査線の位相表現はそれぞれ２πｙ／ｐ、２πｙ／Ｔ0と表せる。
y方向は走査線の方向に対して垂直である。
【００２０】
従前の走査モアレ縞が発生する理由は、走査線と試料の格子間の小さな不一致またはずれであり、Ｔ0はｐに近接している。
１次モアレ縞の位相は次式で表現できる。
【００２１】
【数１】
2017146202.tif 000003
【００２２】
ただしφhは高次周波数の位相(２πｙ／ｐ、２πｙ／Ｔ0、２πｙ(１／ｐ＋１／Ｔ0))としφ0はその初期位相とする。
場合によっては、走査線間隔が試料格子ピッチに近接している条件を満足するのは容易でない。
【００２３】
この条件では、図２に示したように、走査線間隔が試料格子ピッチの整数倍に近接していれば、すなわちＴ(ｎ)がｎｐ(ｎは正の整数で２以上)に近接していれば、高速走査モアレ縞が現われる。
高速走査は走査解像度を落として走査線間隔を大きく走査する事で実現される。
高速走査線間隔は次の式で定義される。
【００２４】
【数２】
2017146202.tif 000004
【００２５】
ただしＴ0は試料格子ピッチに近接した基本走査線間隔とする。
そうすると高速走査モアレ縞は次式で表される。
【００２６】
【数３】
2017146202.tif 000005
【００２７】
ただしφh (n)は高次周波数の位相であって(２πｙ／ｐ、２πｙ／Ｔ0、２πｙ(１／ｐ＋１／Ｔ0)、２πｙ(１／ｐ−１／Ｔ0を含む)、φ0(n)はその初期位相とする。 
【００２８】
１．２  高速走査モアレ法の変形計測原理
モアレ縞はおおよそ周期的分布状態であるので高速走査モアレ縞のモアレ間隔は次式で表される。
【００２９】
【数４】
2017146202.tif 000006
【００３０】
ただし、ｄ(ｎ)はｙ方向のモアレ間隔であり、φh(ｎ)とφ0(n)は高次周波数の位相と高速走査モアレ縞の初期位相とする。
数式(３)と(４)のモアレ位相の基本周波数要素を比較すれば、次式に示す試料格子ピッチと高速走査線間隔と高速走査モアレ間隔の相関を得ることができる。
【００３１】
【数５】
2017146202.tif 000007
【００３２】
数式(２)を数式(５)に代入するとモアレ間隔はｎに依存せず一定であることがわかる。
よって次の２つの数式が得られる。
【００３３】
【数６】
2017146202.tif 000008
【００３４】
【数７】
2017146202.tif 000009
【００３５】
低走査線解像度のモアレ縞は基本走査線間隔の関数であり相対変位の等高線である。
モアレ縞次数をｍとして、基本走査線に関してｙ方向におけるモアレ縞の変位は次式で表される。
【００３６】
【数８】
2017146202.tif 000010
【００３７】
Ｔ0間隔の基本走査線に相対的なピッチｐの試料格子のひずみはｙ方向の相対変位の一次微分から得ることができる。
【００３８】
【数９】
2017146202.tif 000011
【００３９】
ただし方向変数ｙの漸増の半分はひずみを計算する場合において半微分ステップと呼ぶ。
試料格子の実際のひずみは試料格子ピッチ(または基本走査間隔)の変形前後の幾何学的関係から得られる。
実際のひずみは次式によって与えられる。
【００４０】
【数１０】
2017146202.tif 000012
【００４１】
ただし、ｐ0はｙ方向における変形前の試料格子ピッチとする。
【００４２】
１．３  高速走査モアレ法による変形計測感度
高速走査モアレ法ではＴ0間隔の基本走査線に相対的なピッチｐの試料格子のひずみは式、ε(ｎ)＝(ｐ−Ｔ0)／Ｔ0で表される。
この数式と数式(７)とを比較すれば、Ｔ(ｎ)＜ｎｐの場合には式、ε(ｎ)＝Ｔ0／(ｄ−Ｔ0)で得られ、Ｔ(ｎ)＜ｎｐの場合には式、ε(ｎ)＝−Ｔ0／(ｄ＋Ｔ0)で得られる。
モアレ間隔は通常ｙ方向には基本走査線間隔よりはるかに大きいことから相対ひずみは式、|ε(ｎ)|≒Ｔ0／ｄ＝Ｔ(ｎ)／(ｎｄ)で簡素化される。
【００４３】
次に提案法の変形計測精度について理論的に述べる。
数式(８)から提案法の変位感度は縞模様につきＴ(ｎ)／ｎであり、従来の走査モアレ法の縞模様あたりの変位感度と同じである。
ひずみ感度は測定可能な最小のひずみ、すなわち全視野で２個のモアレ縞のみ観測される場合のひずみ、ｄ＝Ｎ(ｎ)Ｔ(ｎ)と定義される。ただしＮ(ｎ)は走査線数とする。
【００４４】
簡素化された相対ひずみ式に代入すれば高速モアレ法のひずみ感度は式、Ｔ(ｎ)／[ｎ Ｎ(ｎ)Ｔ(ｎ)]＝１／[ｎ Ｎ(ｎ)]となる。
Ｔ0がｐに近接している場合の基本走査線数をＮ0とすれば全視野ではｄ＝Ｎ0Ｔ0となるので、ひずみ感度は、Ｔ(ｎ)／[ｎＮ0Ｔ0]＝１／Ｎ0となって、結局従来のモアレ法の感度と同じになることがわかる。
【００４５】
言い換えれば提案法の変位とひずみ計測感度は従来の走査モアレ法と同等である。
公知の位相シフト法と組み合わせることにより変形計測精度を改善する事ができる。
【００４６】
本発明で得られる高速走査モアレ縞は、従来の走査モアレ縞と同じモアレ間隔を有している。
よって高速走査モアレ法は、変位とひずみの測定感度を低下させることなく、高い速度でモアレ画像を取得できることから、従来の走査モアレ法より、試料のクリープを受けにくい。
また変形計測のためにより多くのピッチの試料格子に適用することができる。
【００４７】
２．二次元走査モアレ法
二次元走査モアレ法は、単一の二次元走査モアレ画像から、二方向の変形が同時に測定できることを示している。
【００４８】
最初に、試料表面上に規則的な交差格子を作製する。その後で、試料を走査型顕微鏡の試料台上に載置する。
【００４９】
次に、走査顕微鏡の走査線は多くの走査点で構成しているので、垂直方向と水平方向共に走査点の間隔が等しく、かつ、試料格子ピッチに近接するように走査解像度を調整する。
走査点と交差格子間との干渉から、二次元走査モアレ縞が生成する(図３)。
【００５０】
最後に、二次元走査モアレ法により、全視野で試料の二方向の変位とひずみ分布を同時に測定することができる。
【００５１】
２．１  二次元走査モアレ縞の生成原理
走査顕微鏡の走査プロセスでは並行走査線が知られている。
実際のところ、図３に示すように、走査型顕微鏡では走査線は多くの走査点から構成されている。
【００５２】
試料格子が交差格子である場合、交差格子と走査点との干渉により二次元走査モアレ縞が発生する。ｘ方向の走査点間隔(Ｔx)とｙ方向の走査点間隔(Ｔy)は等しいからＴx＝Ｔy＝Ｔ0となる。
【００５３】
２．２  二次元走査モアレ法による変形計測原理
図４に示すように、二次元走査モアレ縞は２つの従来の一次元走査モアレ縞を二方向に合成したものと考えてよい。
フーリエ変換と逆フーリエ変換を用いて図４に示すように二次元走査モアレ縞をｘ方向とｙ方向の一次元走査モアレ縞に分離することができる。
交差格子の二方向における走査点に相対的な変位は次式で表すことができる。
【００５４】
【数１１】
2017146202.tif 000013
【００５５】
ただし、図４に示すように、モアレ縞次数をｍ、ｘ方向とｙ方向における指数(図の左上を基点とする)を、走査点間隔をＴとする。
走査点に対する試料格子のひずみ、すなわち、二方向の相対ひずみは、それぞれ、二方向での相対変位の一次微分から得ることができる。
【００５６】
【数１２】
2017146202.tif 000014
【００５７】
ただし、方向変数ｉの漸増の半分はひずみを計算する場合において半微分ステップと呼ぶ。
試料格子の実際のひずみは、変形前後の試料格子ピッチや走査点間隔との幾何学的関係(関数)から求めることができる。
試料格子の実際の二方向ひずみは次式による。
【００５８】
【数１３】
2017146202.tif 000015
【００５９】
ただし変形前のｘ方向またはｙ方向の試料格子ピッチをｐi_0とする。
相対ひずみが最初に計測されるのは試料格子を観測することなく走査モアレが低倍率で記録されるので全視野が顕著に拡大できることによる。
より小さな格子を用いれば、より高精細な変位・ひずみ分布が得られる。
【００６０】
せん断ひずみは次式により求められるεxy＝∂ｕx／∂ｙ＋∂ｕy／∂ｘ。試料の実際のせん断ひずみは走査点に対する相対せん断ひずみとなる。
本発明は主にｘ方向とｙ方向の通常ひずみの測定に焦点をあてている。
【００６１】
二次元走査モアレ縞を用いて二方向での変形を測定するために、例えば、フーリエ変換または各種ローパスフィルタの画像処理を使用して、二方向に２つの一次元モアレ縞に分離することができる。
【００６２】
二次元走査モアレ法は、試料や走査方向を回転させることなく、単一の二次元モアレ縞を使用して、二次元の変位とひずみ分布を測定することができる。
従って、この方法は、二次元変形計測のために、簡単な操作、節約の時間、および高精度の利点を有する。
【００６３】
３．高速二次元走査モアレ法
高速二次元走査モアレ法は、高速走査モアレ法と二次元走査モアレ法を融合する方法である。最初に、試料表面上に規則的な交差格子を作製する。
その後で、試料を走査型顕微鏡の試料台上に載置する。
【００６４】
次に、走査解像度を垂直方向の走査点の間隔が試料格子ピッチの整数(≧２)倍数に近くに、水平方向の走査点の間隔が試料格子ピッチの整数(≧１)倍数に近くなるように調整する。
高速走査点と試料交差格子間の干渉は、高速二次元走査モアレ縞を生成する。
【００６５】
最後に、高速二次元走査モアレ法による、全視野で試料の二方向の変位とひずみ分布を同時に迅速に測定することができる。
【００６６】
図５は変形計測における高速二次元走査モアレ縞とモアレ像処理の生成原理を表している。
高速二次元走査モアレ縞は高速走査点と交差格子間の干渉で生成する。
ｙ方向の走査点間隔は基本走査点間隔の２倍以上、かつ、ｘ方向の間隔は基本走査点間隔と等しいまたはその倍数、Ｔi(ｎi)＝ｎi Ｔ0、ｎy≧２, ｎx≧１,である。
ただしＴ0は試料格子ピッチに近接する基本走査点間隔とし、ｉはｘ方向またはｙ方向とする。
基本走査点に相対する二方向の試料格子変位は次式で表される。
【００６７】
【数１４】
2017146202.tif 000016
【００６８】
ただし、図５に示すように、ｉ方向におけるモアレ縞次数をｍi、走査点間隔をＴとする。
走査点に対する試料格子のひずみ、すなわち、二方向の相対ひずみは、それぞれ、二方向での相対変位の一次微分から得ることができる。
【００６９】
【数１５】
2017146202.tif 000017
【００７０】
試料格子の実際のひずみは、変形前後の試料格子ピッチや走査点間隔との幾何学的関係(関数)から求めることができる。
試料格子の実際の二方向ひずみは次式による。
【００７１】
【数１６】
2017146202.tif 000018
【００７２】
高速二次元走査モアレ縞が高速で記録され、二次元走査モアレ縞のような二方向において同じモアレ間隔を有する。
【００７３】
高速二次元走査モアレ法は、高速走査モアレ法と二次元走査モアレ法の両方の利点を有している。
【００７４】
図６に、本発明で述べたこれらの走査モアレ法による変位ひずみ分布の解析フローチャートを示す。
解析手順は３つの走査モアレ法に共通であり、（ステップ１）まず試料表面に格子を作製し、（ステップ２）いずれかのモアレ画像を取得する。
【００７５】
（ステップ３）次に二次元のモアレ縞が含まれる場合は、（ステップ４）二次元走査モアレ縞をｘ方向とｙ方向の一次元走査モアレ縞に分離する前処理を行ってから、（ステップ５）適宜一次元または二次元の変位解析処理を行い、最後にひずみ解析処理を行う（ステップ６）。
【００７６】
図７にレーザ顕微鏡を用いて金属材料の引張試験の実験装置の一例を示す。
これらの方法は、様々な材料や構造の広視野や全視野での非破壊ミクロン/ナノスケールの変形を測定することができる。
【００７７】
４．ＣＣＤカメラ又はＣＭＯＳカメラへの拡張
本発明で述べたこれらの走査モアレ法は、走査線または走査点から、デジタル格子または、電荷結合素子(ＣＣＤ)カメラ又は相補型金属酸化膜半導体(ＣＭＯＳ)カメラの画素に変更することにより、ミリメートル/メートルスケールの変形を測定するために拡張することができる。
図８に示す提案した高速ＣＣＤ又はＣＭＯＳカメラを使用して、高速サンプリングモアレ法よりミリメートル／メートルスケールの変形を迅速に測定することができる。
【発明の効果】
【００７８】
本発明は迅速に正確に簡単に全視野での二次元変位ひずみ分布を測定でき、以下の効果を有する。
効果１：高速走査モアレ法と高速二次元走査モアレ法で、モアレ画像取得速度が高く、使用可能な試料格子ピッチの範囲が広い。
高速変形測定は、計測感度を低下させることなく、変位ひずみ分布を計測できる。
試料クリープの影響を受けにくく、動的測定に適している。
効果２：二次元走査モアレ法と高速二次元走査モアレ法で、二次元変形測定精度は高い。
試料または走査方向の回転動作および２回の走査を必要としないため、回転による誤差がなく、二方向の変形を同時に測定することができる。
また実験操作は簡単で、簡易に動的計測に利用できる。
効果３：開発したこれらの３つの走査モアレ法で、広視野での変位とひずみ分布を非破壊かつ定量的に測定することができる。
格子そのものを観察するではなく、拡大現象であるモアレ縞から変位ひずみ分布を評価しているため、試料格子ピッチの千倍以上の視野を観察できる。
【図面の簡単な説明】
【００７９】
【図１】提案した３つの変形測定方法(高速、二次元、高速二次元走査モアレ法/サンプリングモアレ法)の説明図である。
【図２】(ａ)従来の走査/サンプリングモアレ縞、及び(ｂ)高速走査/サンプリングモアレ縞の形成原理を表す図である。
【図３】二次元走査/サンプリングモアレ縞の生成原理を表す図である。
【図４】二次元走査/サンプリングモアレ縞と二方向に２つの一次元モアレ縞の関係を表す図である。
【図５】高速二次元走査/サンプリングモアレ縞の生成原理、および二方向に２つの一次元モアレ縞との関係を表す図である。
【図６】提案した変形測定方法のフローチャートである。
【図７】提案した変形測定方法の実験装置を撮影した写真である。
【図８】提案した高速ＣＣＤ又はＣＭＯＳカメラ装置の概略図である。図８右上は高速一次元記録モード、図８右下は高速二次元記録モードを表す図である。
【図９】高速走査/サンプリングモアレ縞と従来の走査/サンプリングモアレ縞のシミュレーション結果を表す図である。
【図１０】走査モアレ縞を生成するために使用する、シリコン上にナノインプリントリソグラフィによって作製したマイクロン交差格子を撮影した画像である。
【図１１】レーザ走査型顕微鏡で撮影した、シリコン上に生成した高速走査モアレ縞及び従来の走査モアレ縞を比較する図である。
【図１２】(ａ)従来の走査モアレ法から測定した、(ｂ)高速走査モアレ法から測定した、基本的な走査線に対する相対変形の結果を表す図である。
【図１３】従来の走査モアレ法と提案した高速走査モアレ法について測定した相対ひずみと測定した試料格子ピッチに関する比較検証図である。
【図１４】レーザ走査型顕微鏡で撮影した、シリコン上に生成した二次元走査モアレ縞と高速二次元走査モアレ縞の画像である。
【図１５】レーザ走査顕微鏡で撮影した二次元走査モアレ縞の処理前後の画像。(ａ)は図１２(ａ)中の任意の正方形で標識した二次元走査モアレパターン。(ｂ)と(ｃ)は、各ｘ方向とｙ方向における、二次元走査モアレパターンをフーリエ変換と逆フーリエ変換後でのモアレパターンである。
【図１６】レーザ走査顕微鏡で撮影した高速二次元走査モアレ縞の処理前後の画像。(ａ)は図１２(ｂ)中の正方形で標識した高速二次元走査モアレパターン。(ｂ)と(ｃ)は、各ｘ方向とｙ方向における、高速二次元走査モアレパターンをフーリエ変換と逆フーリエ変換後でのモアレパターンである。
【図１７】(ａ)と(ｂ)は、各提案した二次元走査モアレ法と高速二次元走査モアレ法について測定した相対ひずみと測定した試料格子ピッチに関する比較検証図である。
【図１８】(ａ)と(ｂ)と(ｃ)は、それぞれ提案した高速走査、二次元走査、及び高速二次元走査モアレ法から、計算した試料格子ピッチの相対的な差に関する比較図である。
【発明を実施するための形態】
【００８０】
以下に実施例を示す。
【実施例１】
【００８１】
１．シミュレーションと実験から高速走査モアレ縞の生成
本実施例では、本発明による高速走査モアレ縞の生成方法、及び従来の走査モアレ縞との関係を示す。
【００８２】
図９に走査モアレ法とサンプリングモアレ法に共通して、一次元のモアレ縞の生成に関するシミュレーション結果を示す。
図９は、異なる走査解像度または走査間隔でダウンサンプリングした場合に得られるモアレ縞を示している。
【００８３】
図９(ａ)は試料格子に相当する１３画素のピッチを有する余弦型平行格子である。
これに対して、図９(ｂ)に示すように基本的な走査線間隔Ｔ0を１２画素として、格子試料とＴ0＝１２画素を有する走査格子の重ね合わせから、従来の走査モアレ縞を生成できる。 
【００８４】
次に、図９(ｃ)にＴ(１)＝Ｔ0＝１２画素のサンプリング間隔(走査間隔に該当する)でのダウンサンプリングと輝度補間より得られるサンプリングモアレ縞は図９(ｂ)と同様に従来のモアレ縞を生成することができる。
従って、以降はサンプリングモアレを試料格子と走査線の干渉より得られる重畳モアレと走査モアレの代わりに、異なる走査解像度の場合の比較として使用した。
【００８５】
図９(ｃ)と同様にＴ(２)＝２４画素, Ｔ(３)＝３６画素, Ｔ(４)＝４８画素の異なるサンプリング間隔でダウンサンプリング(縦方向の間引き処理)と輝度補間の画像処理より得られるモアレ縞を、それぞれ図９(ｄ)、図９(ｅ)、図９(ｆ)に示す。
【００８６】
この結果からｎの値(n＝1, 2, 3, 4)に依らず、生成されるモアレ縞のピッチ間隔は変わらないことを確認できる。
すなわち、走査線間隔Ｔ(ｎ)は基本的な走査線間隔Ｔ0のｎ倍であれば、モアレ縞の間隔は変わらないことを意味する。
【００８７】
ただし、図９(ａ)〜図９(ｆ)中に示すモアレ縞画像のグレースケールの範囲は−１から１まで同じであり、得られるモアレ縞画像のコントラストはｎの増加に伴ってわずかに減少する。
しかし、シミュレーション及び実際の実験では、画像のヒストグラムを調整することで、モアレ縞画像のコントラストを高くすることができる。
本シミュレーション結果から、従来の走査モアレ法で得られるモアレ縞の間隔と、高速走査モアレ法で得られるモアレ縞の間隔が変わらないことを確認できる。
【００８８】
レーザ走査型顕微鏡での実験からも、高速走査モアレ縞と従来のモアレ縞の関係を示す観察結果が得られる。
図１０にＵＶナノインプリントリソグラフィよりシリコンウェハ上に作製した0.79μmのピッチを有する交差格子を示す。
【００８９】
作製した試料格子を有するシリコンウェハを、次にレーザ走査型顕微鏡(レーザーテックハイブリッドL3-URW)の試料台上に載置した。
レーザ走査顕微鏡を用いた顕微鏡画像を撮影する際に、ｙ方向に走査線数を1024、512及び256ラインに変更することができる。
【００９０】
20倍の倍率の対物レンズで、ｙ方向に走査線数を512に設定した場合、走査線間隔は試料格子ピッチ(ｐ＝0.79μm)に近いＴ(１)＝Ｔ0＝0.726μmに相当する。
従って、試料格子ピッチ0.79μmと走査線間隔0.726μmの干渉から、図９(ａ)に示すように従来の走査モアレ縞画像を生成できる。
【００９１】
同じ対物レンズを用いて、走査線数を256ラインに設定したときに、走査線間隔はＴ(２)＝２Ｔ0＝1.452μmに該当する。
この場合、試料格子ピッチ0.79μmの倍である1.58μmと走査線間隔1.452μmの干渉から、図１１(ｂ)に示すように低解像度走査モアレ縞画像が得られる。
【００９２】
ここで図１１(ｃ)と図１１(ｄ)に正方形で囲んだ同じ領域での走査モアレ縞の拡大画像を示す。
図１１(ｃ)と図１１(ｄ)から、高速走査モアレ縞と従来の走査モアレ縞のピッチ間隔は同じであり、モアレ縞の輝度値分布もほぼ同じであることを確認できる。
【００９３】
このことから、高速走査モアレ縞を用いて従来の走査モアレ法と同様に材料の変位ひずみ分布を算出できることを示唆している。
【実施例２】
【００９４】
２．高速走査モアレ法の変形測定と検証
この実施例では、従来の走査モアレ法と比較することで、本発明である第１手法の高速走査モアレ法の変形測定精度を確認する。
【００９５】
式(８)及び式(９)による、図１１(ｃ)と図１１(ｄ)にレーザ走査モアレ縞から、基本的な走査線(Ｔ0＝0.726μm)に対して、格子試料の変形分布を算出した。
【００９６】
ここでは２つのモアレ縞の間の変位を、１方向の放物線補間(２次補間)法から得た。
半微分ステップの値が10である場合、図１１(ｃ)に示すＴ(１)＝0.726 μmで得られたモアレ縞画像から算出した変位分布とひずみ分布をそれぞれ図１２(ａ)と図１２(ｂ)にプロットした。
【００９７】
また図１１(ｄ)に示すＴ(２)＝1.452μmで得られたモアレ縞画像から算出した変位分布とひずみ分布をそれぞれ図１２(ｃ)と図１２(ｄ)にプロットした。
基本的な走査線(Ｔ0＝0.726μm)に対して、異なる半微分ステップにおける、図１３(ａ)に示すＴ(１)＝0.726μmから得られた平均ひずみ値は0.08992であり、Ｔ(２)＝1.452μmから得られた平均ひずみ値は0.08981である。
【００９８】
従来の走査モアレ法(Ｔ(１)＝Ｔ0＝0.726 μm)の結果と比較すると、Ｔ(２)＝1.452 μmのときの高速走査モアレ法からのひずみの相対差はわずか0.12％である。
これらの変位分布とひずみ分布間のわずかな差異は、主に画像撮影時に発生した振動と外部ノイズに起因するものと考えられる。
【００９９】
図１３(ｂ)に示すように、高速走査モアレ法と従来の走査モアレ法を用いて得られたひずみから算出した試料格子ピッチの差は0.2nm以下である。
【０１００】
これらの結果から、本発明である第１手法の高速走査モアレ法は高速でありながら、従来の走査モアレ法と同程度の計測精度で変位ひずみ分布を算出できることを実証し、本手法の有効性を示した。
【実施例３】
【０１０１】
３．二次元(以下、２Ｄとも言う)走査モアレ縞と高速２Ｄ走査モアレ縞の生成
この実施例では、レーザ走査型顕微鏡で二次元走査モアレ縞と高速二次元走査モアレ縞の生成結果を示している。
【０１０２】
二次元走査モアレ縞を観察するためには、交差格子は試料格子になることが不可欠である。
先に実施例１と同様に、ＵＶナノインプリントリソグラフィによって、シリコンウェハ上に0.79μmのピッチを有する交差格子(図１０と同じ)を作製した。
このとき使用したナノインプリントのレジストはPAK01であった。
【０１０３】
次に交差試料格子をレーザ走査型顕微鏡で観察した。
ｙ方向の走査点数を512に設定した。
対物レンズの倍率は20倍であるときに、ｙ方向の走査点間隔は試料格子ピッチ(T<P)より少し小さい0.726μmであった。
ｘ方向の走査点間隔は、ｙ方向の走査点の間隔に等しい。
【０１０４】
試料格子の主方向の一方が水平方向とほぼ平行であったとき、試料格子とレーザ走査点の間の干渉から、図１４(ａ)に示すような二次元レーザ走査モアレ縞を生成できる。
【０１０５】
ここで同じく20倍率の対物レンズを用いて走査点数をｙ方向に256ラインに設定した場合、ｙ方向の走査点間隔は試料格子ピッチの２倍に近い1.452μmであった。
ｘ方向の走査点間隔はまだ0.726μmのままである。
この場合、走査速度は、従来の走査速度より２倍高い。
【０１０６】
試料格子の主方向の一方が水平方向とほぼ平行であるとき、試料格子とレーザ高速走査点の間の干渉から、図１４(ｂ)に示すような高速二次元レーザ走査モアレ縞を得ることができる。
そのモアレ縞のピッチ間隔は従来の走査モアレ法と変わらない。
【実施例４】
【０１０７】
３．２Ｄ走査モアレ法と高速２Ｄ走査モアレ法の変形測定と検証
この実施例では、従来の走査モアレ法と比較して、本発明である第２手法の二次元走査モアレ法と、本発明である第３手法の高速二次元走査モアレ法の変形測定精度を確認する。
二次元走査モアレ法による相対変形を測定するため、図１４(ａ)中の正方形領域を評価した。
【０１０８】
図１４(ａ)に示すように、二次元走査モアレ縞の評価領域の面積は217.7×217.7 μm2であった。
フーリエ変換と逆フーリエ変換の画像処理を用いて適切な周波数分離処理を施すことで、図１５(ａ)の二次元モアレパターンをそれぞれ、ｘ方向の一次元モアレ縞画像(図１５(ｂ))とｙ方向の一次元モアレ縞画像(図１５(ｃ))に分離することができる。
【０１０９】
また図１６(ａ)に、高速二次元走査モアレ法による相対変形を測定するため、図１４(ｂ)中の正方形領域を示す。
なお、図１５(ａ)と図１６(ａ)の領域は同じである。
フーリエ変換と逆フーリエ変換の画像処理を用いて適切な周波数分離処理を施すことで、図１６(ａ)の二次元モアレパターンをそれぞれ、ｘ方向の一次元モアレ縞画像(図１６(ｂ))とｙ方向の一次元モアレ縞画像(図１６(ｃ))に分離することができる。
【０１１０】
式(１１)によれば、ｘ方向及びｙ方向の走査点に対して試料格子の変位が得られる。
式(１２)を用いて、これら２つの方向の走査点に対して試料格子のひずみを算出できる。
図１７に二次元走査モアレ法と高速二次元走査モアレ法より算出した相対ひずみの平均値を示す。
【０１１１】
この２つのモアレ法から測定したｘ方向又はｙ方向の相対ひずみはほぼ同じであることがわかる。
ｘ方向の相対ひずみは0.091であり、ｙ方向の相対ひずみは0.090である。
両方が従来の走査モアレ法から得られたひずみ値である0.090とよく一致している。
【０１１２】
二次元走査モアレ法と高速二次元走査モアレ法の変形測定精度も、算出した試料格子ピッチから検証することができる。
【０１１３】
図１７(ｂ)に示すように、二次元走査モアレ法と高速二次元走査モアレ法から算出した試料格子ピッチがほぼ同じである。
ここで測定した試料格子ピッチは、ｘ方向では0.792μmであり、ｙ方向では0.791μmである。
ｘ方向とｙ方向にピッチの差はわずか1nmのみであることから高精度な測定が実現できているといえる。
【０１１４】
図１８に、従来の走査モアレ法より測定した試料格子ピッチ0.7904μmに対して、二次元走査モアレ法と高速二次元走査モアレ法より測定した試料格子ピッチの相対差をプロットした。
図１８に示すとおり、算出した試料格子ピッチの相対差はｘ方向では0.2％であり、ｙ方向では0.05％である。
【０１１５】
これらの結果から本発明である第２手法の二次元走査モアレ法、と本発明である第３手法の高速二次元走査モアレ法は、従来の走査モアレ法と同程度の計測精度で変位・ひずみ分布を算出できることを実証し、本手法の有効性を示した。
【０１１６】
さらに、二次元走査モアレ法は、試料や走査方向を回転させることなく、単一の二次元モアレ縞を使用して、二次元の変位とひずみ分布を測定することができる。
高速二次元走査モアレ法は、高速と簡便な二次元変形測定の両方の利点を有する。
【実施例５】
【０１１７】
４．高速ＣＣＤ又はＣＭＯＳカメラへの拡張
図８は提案した高速ＣＣＤ又はＣＭＯＳカメラ装置の概略図である
【０１１８】
図８右上は高速一次元モアレ縞画像記録モード、図８右下は高速二次元モアレ縞画像記録モードを表す図である。
ＣＣＤ／ＣＭＯＳ検出装置(イメージセンサ)を用い、縦横に整列した基本撮像素子間隔の縦または横、または縦横の各整数倍の撮像素子列に記録し、記録された輝度情報を取り出して、実施例１乃至実施例４に示したように、モアレ縞画像処理すればよい。
【０１１９】
また、全撮像素子に記録する公知の記録モードでは、測定物の表面に生成されたピッチＰの規則格子または規則交差格子のイメージセンサにおけるピッチｐに従って、実施例１乃至実施例４に示したように、基本撮像素子間隔の縦または横、または縦横の各整数倍の撮像素子列に記録し、記録された輝度情報を取り出してアレ縞画像処理すればよい。
こうすると、画像処理の高速化を図る事ができる。
【０１２０】
また以上の実施例の画像処理は、パーソナルコンピュータにおいてプログラムを作成しておよび画像処理ツールを用いて行ったが、当該処理手段を備える走査型顕微鏡装置上で行ってもよいし、また、そのプログラム等を分割して走査型顕微鏡装置とパーソナルコンピュータで協同して行ってもよい。
【０１２１】
本発明の画像の演算処理は、パーソナルコンピュータ上のＣ/Ｃ++プログラムにより実行して計測結果を求めたが、プログラムの種類は本発明を実施する計測実行環境により適宜変更することができる。
【０１２２】
またプログラムは内蔵式、組み込み式(Imbedded)、読み込み式、あるいはダウンロード方式を問わない。
【産業上の利用可能性】
【０１２３】
本発明である高速二次元走査モアレ法は、航空宇宙、自動車、電子パッケージ、バイオ医薬などの産業分野における材料の全視野変位ひずみ測定に有用である。
好適には、強化プラスチック、金属複合材料、セラミック複合材料、積層材料、フィルム/基板構造、パッケージ構造、などの金属、ポリマー、セラミック、半導体、複合材料及び構造を含む各種材料の変形を解析するのに適している。
特に以下の３つの材料評価に有効である。
１．亀裂発生、亀裂伝播および損傷挙動の評価
２．界面最適設計のための界面接合性能の評価
３．機械特性、残留応力および構造物の健全性評価
【符号の説明】
【０１２４】
１  試料表面の格子
２，３，４，５  参照格子
６  走査線・サンプリング列
７  公知の走査モアレ縞・公知のサンプリングモアレ縞
８  高速走査線・高速サンプリング列
９  高速走査モアレ縞・高速サンプリングモアレ縞
１０  走査点・サンプリング点
１１  主方向がｘ方向の走査線・サンプリング列
１２  主方向がｙ方向の走査線・サンプリング列
１３  顕微鏡(走査型レーザ顕微鏡、走査型電子顕微鏡、走査型プローブ顕微鏡など)
１４  機械荷重(引張、圧縮、３点曲げ)
１５  並進・回転荷重
１６  熱負荷
１７  電気負荷
１８  ＣＣＤ／ＣＭＯＳ検出装置(イメージセンサ)
１９  対物レンズ
２０  公知の記録モード
２１  高速一次元モアレ縞記録モード
２２  高速二次元モアレ縞記録モード
【図１】
2017146202.tif 000028
【図２】
2017146202.tif 000029
【図３】
2017146202.tif 000030
【図４】
2017146202.tif 000031
【図５】
2017146202.tif 000032
【図６】
2017146202.tif 000033
【図７】
2017146202.tif 000034
【図８】
2017146202.tif 000035
【図９】
2017146202.tif 000036
【図１０】
2017146202.tif 000037
【図１１】
2017146202.tif 000038
【図１２】
2017146202.tif 000039
【図１３】
2017146202.tif 000040
【図１４】
2017146202.tif 000041
【図１５】
2017146202.tif 000042
【図１６】
2017146202.tif 000043
【図１７】
2017146202.tif 000044
【図１８】
2017146202.tif 000045