JP6583761B2 - 周期模様を利用した三次元形状・変位・ひずみ測定装置、方法およびそのプログラム - Google Patents
周期模様を利用した三次元形状・変位・ひずみ測定装置、方法およびそのプログラム Download PDFInfo
- Publication number
- JP6583761B2 JP6583761B2 JP2018541896A JP2018541896A JP6583761B2 JP 6583761 B2 JP6583761 B2 JP 6583761B2 JP 2018541896 A JP2018541896 A JP 2018541896A JP 2018541896 A JP2018541896 A JP 2018541896A JP 6583761 B2 JP6583761 B2 JP 6583761B2
- Authority
- JP
- Japan
- Prior art keywords
- moire
- phase
- lattice
- distribution
- displacement
- 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
- 238000000034 method Methods 0.000 title claims description 176
- 238000006073 displacement reaction Methods 0.000 title claims description 44
- 230000000737 periodic effect Effects 0.000 title description 2
- 238000009826 distribution Methods 0.000 claims description 71
- 238000005070 sampling Methods 0.000 claims description 45
- 238000004364 calculation method Methods 0.000 claims description 13
- 238000000691 measurement method Methods 0.000 claims description 8
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 claims 1
- 238000005259 measurement Methods 0.000 description 80
- 239000011295 pitch Substances 0.000 description 71
- 230000010363 phase shift Effects 0.000 description 34
- 238000004458 analytical method Methods 0.000 description 25
- 239000000463 material Substances 0.000 description 15
- 238000004088 simulation Methods 0.000 description 12
- 230000008569 process Effects 0.000 description 11
- 229910052782 aluminium Inorganic materials 0.000 description 9
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 9
- 230000008859 change Effects 0.000 description 8
- 238000010586 diagram Methods 0.000 description 8
- 238000005305 interferometry Methods 0.000 description 8
- 230000000694 effects Effects 0.000 description 7
- 238000012795 verification Methods 0.000 description 7
- 238000003384 imaging method Methods 0.000 description 5
- 238000011156 evaluation Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000007796 conventional method Methods 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 239000002131 composite material Substances 0.000 description 2
- 238000012790 confirmation Methods 0.000 description 2
- 230000032798 delamination Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000001066 destructive effect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 229910052751 metal Inorganic materials 0.000 description 2
- 239000002184 metal Substances 0.000 description 2
- 150000002739 metals Chemical class 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 229920000642 polymer Polymers 0.000 description 2
- 239000004065 semiconductor Substances 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 1
- 239000000919 ceramic Substances 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000004100 electronic packaging Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- -1 hybrid structures Substances 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000007769 metal material Substances 0.000 description 1
- 229910044991 metal oxide Inorganic materials 0.000 description 1
- 150000004706 metal oxides Chemical class 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000001127 nanoimprint lithography Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000003908 quality control method Methods 0.000 description 1
- 230000003252 repetitive effect Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000009864 tensile test Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000004154 testing of material Methods 0.000 description 1
- 239000010409 thin film Substances 0.000 description 1
- 230000005428 wave function Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/16—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge
- G01B11/161—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge by interferometric means
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/16—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/16—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge
- G01B11/165—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge by means of a grating deformed by the object
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/24—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
- G01B11/25—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures by projecting a pattern, e.g. one or more lines, moiré fringes on the object
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/24—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
- G01B11/25—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures by projecting a pattern, e.g. one or more lines, moiré fringes on the object
- G01B11/254—Projection of a pattern, viewing through a pattern, e.g. moiré
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J9/00—Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength
- G01J9/02—Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength by interferometric methods
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L1/00—Measuring force or stress, in general
- G01L1/24—Measuring force or stress, in general by measuring variations of optical properties of material when it is stressed, e.g. by photoelastic stress analysis using infrared, visible light, ultraviolet
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/41—Refractivity; Phase-affecting properties, e.g. optical path length
- G01N21/45—Refractivity; Phase-affecting properties, e.g. optical path length using interferometric methods; using Schlieren methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/50—Depth or shape recovery
- G06T7/521—Depth or shape recovery from laser ranging, e.g. using interferometry; from the projection of structured light
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J9/00—Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength
- G01J9/02—Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength by interferometric methods
- G01J2009/0234—Measurement of the fringe pattern
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/41—Refractivity; Phase-affecting properties, e.g. optical path length
- G01N21/45—Refractivity; Phase-affecting properties, e.g. optical path length using interferometric methods; using Schlieren methods
- G01N2021/456—Moire deflectometry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/84—Systems specially adapted for particular applications
- G01N21/88—Investigating the presence of flaws or contamination
- G01N21/8806—Specially adapted optical and illumination features
- G01N2021/8829—Shadow projection or structured background, e.g. for deflectometry
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Theoretical Computer Science (AREA)
- Optics & Photonics (AREA)
- Length Measuring Devices By Optical Means (AREA)
Description
高精度、広視野で三次元形状・変位・ひずみ分布測定のために、２段モアレ技術とそのプログラムを開発することを目的とする。これは、非破壊計測、実験力学と光学技術の分野に属する。
物体の形状や変形を計測することは、力学的特性、構造不安定性、亀裂伝播と残留応力の評価のため極めて重要な技術である。製品の小型化や集積化に伴い、ミリ、ミクロン、ナノスケールレベルでの変形計測は材料科学の分野で大きな注目を集めている。
現時点で物体の三次元形状計測やその微小変形計測に一般的に使用される手法は、光学的全視野計測法として、デジタル画像相関(Digital Image Correlation, DIC)法、電子スペックルパターン干渉法(Electronic Speckle Pattern Interferometry, ESPI)、幾何学的位相分析(Geometric Phase Analysis, GPA)、縞画像の位相解析方法(例えば、フーリエ変換法、位相シフト法、時空位相シフト法等)、とモアレ法がある。
DIC法は、高い測定精度を有するが測定中外乱ノイズの影響を受けやすいのが欠点である。ESPI は振動に非常に敏感である。GPAは試料の格子を直接観察することから走査モアレ法と比較して観察視野が狭い。フーリエ変換法は平面内の小規模な変形を分析する別の方法であるが、使用する格子を観察することから視野が小さくなる (視野の幅は格子ピッチの通常100倍未満)。
モアレ法は、変形の前後のモアレ縞の変化から変形分布を非破壊的測定する効果的な方法である。中でも、顕微鏡モアレ法、CCDまたはCMOSモアレ法(以降、CCDモアレ法と記す)、モアレ干渉法、サンプリングモアレ法、デジタルモアレ法、オーバーラップ(幾何学的)モアレ法は大きな注目を集め、様々な材料および構造の変形計測に適用されている。
顕微鏡モアレ法とモアレ干渉法は、顕微鏡を用いたミクロン/ナノスケール変形計測に適している。CCDモアレ法、サンプリングモアレ法、デジタルモアレ法、(幾何学的)オーバーラップモアレ法は格子画像を取得できる場合、ナノスケールからメートルスケールの変形を測定できる。
モアレ法は、金属、高分子、複合材料などの材料のひずみや橋梁、建物の構造の変形の測定に適用される。また格子画像のモアレ縞の位相を解析することで、変形計測の測定精度を大幅に向上させることができる。
縞格子(fringe or grating、以降単に縞と呼ぶ)画像の位相解析として、図１のように、複数枚の位相がシフトされた縞画像を用いて位相を解析する位相シフト法や、１枚の縞画像のみを用いて位相を解析するサンプリングモアレ法、と位相シフト法とサンプリングモアレ法を融合した時空位相シフト法がある。
位相シフト法と時空位相シフト法は複数枚の縞画像を撮影する必要があり、動的計測への適用ができない。サンプリングモアレ法は１枚の縞画像から位相を解析でき、動的計測への適用が可能であるが、ノイズが多い環境では、精度の良い計測ができないといった問題点がある。１枚の縞画像を撮影するだけで、高精度な三次元形状変形計測方法が望まれている。
Qinghua Wang, Hiroshi Tsuda and Huimin Xie, Developments and Applications of Moire Techniques for Deformation Measurement, Structure Characterization and Shape Analysis, Vol. 8, No. 3, pp. 188-207 (2015).
顕微鏡モアレ法、CCDモアレ法、モアレ干渉法、サンプリングモアレ法、デジタルモアレ法、(幾何学的)オーバーラップモアレ法を含む従来のモアレ法では、視野と測定精度のバランスをとることは困難である。
顕微鏡モアレ法、CCDモアレ法、モアレ干渉法、デジタルモアレ法、オーバーラップモアレ法において、主に使用される解析法は、モアレ縞(フリンジ)を中心とした方法である。変形分布の計算にモアレ縞の中心線のみが使用されるため、変形の測定精度は高くない。
特許文献１に記載している時空間の位相シフト法を用いる場合は、顕微鏡とCCDモアレ法での測定確度を改善することができる。ただし、位相がシフトした複数枚のモアレ縞や格子イメージを記録する操作は、時間がかかり、振動によって容易に影響を受ける動的測定には適していない。
サンプリングモアレ法では、空間位相シフト法を利用して測定精度を向上できるが、高精度と広視野を同時に実現できない。より小さい格子ピッチを有する試料は、高い測定精度を与える反面、視野が小さくなる。より大きい格子ピッチ(3画素以上)を有する試料は、広視野にできるが、変形を計算する空間位相シフト法を実行できない。
高い変形計測精度と広視野(図２)、または高精度と動的計測を同時に達成することが課題となる。
広視野と高精度、または高精度と動的計測を同時に達成するための革新的な変形計測技術としてこの特許で２段モアレ法を提案する。
コアとなるアイデアは、記録または生成された従来のモアレ縞を、位相シフトされた２段モアレ縞を生成するための格子として扱うことである。２段モアレ縞の位相は空間位相シフト法または時空位相シフト法などのアルゴリスムを利用してより正確に算出することができる。
測定対象物表面と基準面に投影された格子の位相差から形状を、物体表面に貼り付けた格子の変形前後の２段モアレ縞の位相差から変位とひずみを決定することができる。
その結果、物体の三次元形状と変形分布を広視野で正確に、もしくは高精度で動的に測定可能になる。
その結果、物体の三次元形状と変形分布を広視野で正確に、もしくは高精度で動的に測定可能になる。
２段モアレ縞の形成
顕微鏡の走査線、またはカメラのピクセル配列、または仮想格子(ダウンサンプリングピクセル配列)、またはデジタル格子は、参照格子とすることができる。
試料の格子と参照格子間のわずかな不一致またはずれがある場合、従来のモアレ縞が試料の格子と参照格子との干渉から出現する(図３(a))。
顕微鏡の走査線、またはカメラのピクセル配列、または仮想格子(ダウンサンプリングピクセル配列)、またはデジタル格子は、参照格子とすることができる。
試料の格子と参照格子間のわずかな不一致またはずれがある場合、従来のモアレ縞が試料の格子と参照格子との干渉から出現する(図３(a))。
従来のモアレ縞には、顕微鏡モアレ縞、CCDモアレ縞、サンプリングモアレ縞またはデジタルモアレ縞が含まれている。
これら従来のモアレ縞を格子として扱う。
これら従来のモアレ縞を格子として扱う。
従来のモアレ縞の間隔 Pm に近いピッチT(2)でダウンサンプリング処理を行い、ダウンサンプリングの開始点を1 ピクセルずつ(T(2)-1)回ずらして、それぞれのダウンサンプリング画像に対して輝度補間を行えば、図３(b)のように T(2)ステップの位相がシフトされた２段モアレ縞を得ることができる。
２段モアレ縞を生成する方法、およびその位相分布を算出する方法は3つある。以下、それぞれの方法を記す。
２段モアレ縞の位相分布の算出原理１(モアレ縞を格子として扱う場合)
図４に１つ目の２段モアレ縞の生成、および位相分布を算出する方法の原理を示す。
試料の格子のピッチと参照格子(顕微鏡走査線、またはカメラ撮像画素デジタル格子)のx方向(またはy方向、以降省略)のピッチをそれぞれPとTと定義する。
図４に１つ目の２段モアレ縞の生成、および位相分布を算出する方法の原理を示す。
試料の格子のピッチと参照格子(顕微鏡走査線、またはカメラ撮像画素デジタル格子)のx方向(またはy方向、以降省略)のピッチをそれぞれPとTと定義する。
試料格子および参照格子の位相表現はそれぞれφs=2πx/P(またはφs=2πy/P)とφr=2πx/T(またはφr=2πy/T)になる。
A(x,y)をモアレ縞の振幅、φ0を初期位相とし、B(x,y)には背景輝度と高次周波数(2πx/P、2πx(1/P+1/T)、など)が含まれるとして、従来のモアレ縞の輝度分布は式(1)で表せる。試料の格子ピッチPと参照格子またはダウンサンプリングピッチTのわずかな違いによって、ピッチ間隔Pm(1/Pm=1/P-1/T)のモアレ縞が生成される。
A(x,y)をモアレ縞の振幅、φ0を初期位相とし、B(x,y)には背景輝度と高次周波数(2πx/P、2πx(1/P+1/T)、など)が含まれるとして、従来のモアレ縞の輝度分布は式(1)で表せる。試料の格子ピッチPと参照格子またはダウンサンプリングピッチTのわずかな違いによって、ピッチ間隔Pm(1/Pm=1/P-1/T)のモアレ縞が生成される。
高次周波数を除外すれば、従来のモアレ縞の位相は式(2)で表せる。
形状変形計測精度を向上させるために、従来のモアレ縞はサンプリングモアレ法における空間的位相シフト法により解析される。
従来のモアレ縞を位相がシフトされた２段モアレ縞を生成する際に用いる「格子」として扱う。
ダウンサンプリングピッチをT(2) とすると、A(2)(x,y)を２段モアレ縞の振幅、kを位相シフト手順の番号とし、B(2)(x,y)には２段モアレ縞の背景輝度と高次周波数(2πx/P、2πx(1/P+1/T)、など)が含まれるとして、２段モアレ縞の輝度は式(３)で表現できる。
その結果、２段モアレ縞の位相は式(４)で表され、１次元の離散的フーリエ変換(Discrete Fourier transform: DFT)または高速フーリエ変換(Fast Fourier transform: FFT)より算出できる。
ただし、回転因子Wを式(５)として定義する。
ここで求めた２段モアレ縞の位相を利用して、元のモアレ縞φmを式(６)、元の格子の位相φsを式(７)より求めることができる。
このように、従来のモアレ縞を格子とみなし、空間的位相解析方法であるサンプリングモアレ法を適用することで、本来の広視野という特徴に加えて、高精度な解析も同時にできるようになる。
２段モアレ縞の位相分布の算出原理２(格子画像を２段階にダウンサンプリングする場合)
図５に２つ目の２段モアレ縞の位相分布を算出する方法の原理を示す。
ここでは、物体表面に投影された１次元の垂直方向の格子または試料表面に貼り付けた格子のピッチをPとする。このときに撮影される格子画像の輝度分布は式(8)で表される。
図５に２つ目の２段モアレ縞の位相分布を算出する方法の原理を示す。
ここでは、物体表面に投影された１次元の垂直方向の格子または試料表面に貼り付けた格子のピッチをPとする。このときに撮影される格子画像の輝度分布は式(8)で表される。
a(x,y)、b(x,y)、φ0はそれぞれ格子の振幅、背景輝度と初期位相である。
撮影された格子画像に対して、x方向にサンプリング(カメラ撮像画素デジタル格子)ピッチTでダウンサンプリングと輝度補間すると、格子およびサンプリング格子の位相表現はそれぞれφs=2πx/Pとφr=2πx/Tになる。
撮影された格子画像に対して、x方向にサンプリング(カメラ撮像画素デジタル格子)ピッチTでダウンサンプリングと輝度補間すると、格子およびサンプリング格子の位相表現はそれぞれφs=2πx/Pとφr=2πx/Tになる。
A(x,y)をモアレ縞の振幅、φ0を初期位相とし、B(x,y)には背景輝度と高次周波数(2πx/P、2πx(1/P+1/T)、など)が含まれるとして、サンプリングモアレ法によって得られる従来のモアレ縞の輝度は式(9)で表せる。
式(1)と異なり、ここでは同時に複数枚の位相がシフトされたモアレ縞が得られている。
式(1)と異なり、ここでは同時に複数枚の位相がシフトされたモアレ縞が得られている。
ここで得られた複数枚の位相がシフトされたモアレ縞を２段モアレ縞の生成に用いる「格子」として扱う。それぞれの位相シフトモアレ縞に対して、さらに２段階目のダウンサンプリングとして、ダウンサンプリングピッチをT(2)とすると、A(2)(x,y)を２段モアレ縞の振幅、kとtを位相シフト手順の番号とし、B(2)(x,y)は２段モアレ縞の背景輝度として、２段モアレ縞の輝度分布は式(10)で表現できる。
この場合、T×T(2)枚の位相がシフトされた２段モアレ縞が得られる。
この場合、T×T(2)枚の位相がシフトされた２段モアレ縞が得られる。
その結果、２段モアレ縞の位相は式(11)で表され，２次元離散的フーリエ変換または高速フーリエ変換より算出できる。
ここで求めた２段モアレ縞の位相を利用して、元のモアレ縞φmを式(6)、元の格子の位相φsを式(7)より求めることができる。
また、物体表面に投影された１次元の斜め方向の格子または試料表面に貼り付けた格子をデジタルカメラで撮影して得られる格子画像の輝度分布は式(12)で表される。
ここで、Px、Pyはそれぞれ格子のx方向とy方向のピッチ間隔である。
撮影された格子画像に対して、x方向とy方向にそれぞれサンプリング(カメラ撮像画素デジタル格子)ピッチTxとTyでダウンサンプリングと輝度補間すると、２次元サンプリングモアレ法によって得られる従来のモアレ縞の輝度分布は式(13)で表せる。
撮影された格子画像に対して、x方向とy方向にそれぞれサンプリング(カメラ撮像画素デジタル格子)ピッチTxとTyでダウンサンプリングと輝度補間すると、２次元サンプリングモアレ法によって得られる従来のモアレ縞の輝度分布は式(13)で表せる。
ここで得られた複数枚の位相がシフトされたモアレ縞を２段モアレ縞の生成に用いる「格子」として扱う。
それぞれの位相シフトモアレ縞に対して、さらに２段階目のダウンサンプリングとして、x方向とy方向のダウンサンプリングピッチをTx (2)、Ty (2)とすると、A(2)(x,y)を２段モアレ縞の振幅、kxとkyをそれぞれx方向とy方向の位相シフト手順の番号とし、B(2)(x,y)は２段モアレ縞の背景輝度として、複数枚の位相がシフトされた２段モアレ縞の輝度分布は式(14)で表現できる。
それぞれの位相シフトモアレ縞に対して、さらに２段階目のダウンサンプリングとして、x方向とy方向のダウンサンプリングピッチをTx (2)、Ty (2)とすると、A(2)(x,y)を２段モアレ縞の振幅、kxとkyをそれぞれx方向とy方向の位相シフト手順の番号とし、B(2)(x,y)は２段モアレ縞の背景輝度として、複数枚の位相がシフトされた２段モアレ縞の輝度分布は式(14)で表現できる。
これらの２段モアレ縞を用いて、２段モアレ縞の位相を式(15)に示すとおり、(多次元)離散的フーリエ変換または高速フーリエ変換より算出できる。
ここで求めた２段モアレ縞の位相を利用して、元のモアレ縞φmを式(6)、元の格子の位相φsを式(7)より求めることができる。
２段モアレ縞の位相分布の算出原理３(格子画像をダウンサンプリング後にモアレ縞のピッチを調整する方法)
図６に３つ目の２段モアレ縞のピッチ間隔を調整してから位相分布を算出する方法の原理を示す。
本発明である２段モアレ法は、１回目のダウンサンプリングより生成されるモアレ縞を格子とみなし、再度ダウンサンプリング処理を行うが、モアレ縞のピッチ間隔が疎な場合、２段モアレ縞のダウンサンプリングピッチが大きくなり、空間分解能が低下することがある。
本発明である２段モアレ法は、１回目のダウンサンプリングより生成されるモアレ縞を格子とみなし、再度ダウンサンプリング処理を行うが、モアレ縞のピッチ間隔が疎な場合、２段モアレ縞のダウンサンプリングピッチが大きくなり、空間分解能が低下することがある。
この場合、適切なピッチ間隔で２回目のダウンサンプリング処理をできるように、以下にモアレ縞の位相値を利用して、モアレ縞のピッチを調整できる方法を述べる。
従来のサンプリングモアレ法で、格子ピッチPに近いピッチTでダウンサンプリングと輝度補間処理を行い、得られた疎なモアレ縞の位相値に、式(16)に示す整数であるNを乗算した値を再度正弦波または余弦波関数より、密な(乗算型)モアレ縞を生成できる。
この方法で得られた調整後のモアレ縞の間隔に近い値のダウンサンプリングピッチで２段モアレ縞を生成し、これまでに述べた２段モアレ縞の位相算出方法に適用すれば、高精度な２段モアレ縞の位相解析が行える。
この方法において、調整係数であるNを変えることで、精度を優先するのかまたは空間分解能を優先するのかを選択的に決めることができるのが利点である。
表１にこれまでに述べた本発明である３つの手法の特徴と適用条件および従来の位相解析方法との比較をまとめた。記録される格子そのものまたはモアレ縞の間隔が密か粗であるかによって最適な方法を選ぶ。表２に３次元形状・面外変位計測における従来法と本発明の比較と、表３に面内変位・ひずみ計測における従来法と本発明の比較をまとめている。
従来法に比べて、高精度かつ高速、または高精度かつ広視野の形状変形計測ができるようになる。
従来法に比べて、高精度かつ高速、または高精度かつ広視野の形状変形計測ができるようになる。
２段モアレ縞による物体の三次元形状計測の原理
物体を測定する方法の１例として、例えば、平面物体である基準面を格子投影と撮影を行い、基準面の位相分布φrefとし、同様に基準面の前に設置した物体を格子投影と撮影を行い、物体の位相分布をφobjectとする。ここで、基準面と物体の位相差と物体の高さの関係は式(17)に示すように、近似的に比例関係にある。
物体を測定する方法の１例として、例えば、平面物体である基準面を格子投影と撮影を行い、基準面の位相分布φrefとし、同様に基準面の前に設置した物体を格子投影と撮影を行い、物体の位相分布をφobjectとする。ここで、基準面と物体の位相差と物体の高さの関係は式(17)に示すように、近似的に比例関係にある。
したがって、基準面を奥行き(高さ)方向にあらかじめ数カ所動かしながら、それぞれの２カ所の距離(高さ)とそのときの位相差の関係をキャリブレーションして、位相・高さ換算係数κを求めておけば、物体の高さを非接触で測定することができる。
２段モアレ縞による変形(変位・ひずみ分布)計測の原理
構造材料に外力(力や熱変化など)を負荷して変形させると、試料表面の格子のピッチが変形のためx方向にPからP'に変わる。
参照格子Tのピッチが不変の場合、式(2)と同様に、変形後の従来のモアレ縞は式(18)になる。
構造材料に外力(力や熱変化など)を負荷して変形させると、試料表面の格子のピッチが変形のためx方向にPからP'に変わる。
参照格子Tのピッチが不変の場合、式(2)と同様に、変形後の従来のモアレ縞は式(18)になる。
試料格子のピッチがPからP'に変わると、試料の試料格子の位相がφs=2πy/P'に変わるので、変形前後の試料の位相差は式(19)になる。
式(2)、(18)、(19)から、従来のモアレ縞の位相差は変形前後の試料格子の位相差に等しい。
変形に伴い、試料の格子のピッチが変形のためx方向にPからP'に変更し、参照格子T(2)のピッチが不変の場合、式(4)と同様に、変形後の２段モアレ縞は式(21)になる。
式(4)、(20)、(21)により、２段モアレ縞の位相差は、変形前後の従来のモアレ縞の位相差と変形前後の試料格子の位相差に等しいことがわかる。
uをx方向の試料格子の変位を表現すると仮定すると、変形した試料格子の位相はφs'(x)=2π(x-u)/Pとして表せる。
変形前の試料格子の位相φs(x)=2πx/pに基づいて、試料格子の位相差はΔφs(x)=φs’(x)-φs(x)=-2πu/pによって計算できる。その結果、変位は試料格子の位相差に比例する。
ここで、pは試料表面の格子の実寸ピッチである。
x方向の試料の変位は、式(23)によって決定される。
ここで、pは試料表面の格子の実寸ピッチである。
x方向の試料の変位は、式(23)によって決定される。
x方向の試料のひずみは、変位uxの一次微分から測定可能である。
同様にy方向の格子ピッチもpである場合、y方向の変位とひずみ測定は、式(25)と式(26)より算出できる。
(第１の実施形態：三次元形状計測)
以下、図面を参考しながら本発明の第１の実施形態について説明する。
以下、図面を参考しながら本発明の第１の実施形態について説明する。
図７は、本実施形態に係る物体の三次元形状測定システムにおける計算機１、撮像素子３と格子投影装置４、および測定対象物５と基準面(平面物体)６の配置例を示す図である。
図７に示す例では、三次元形状測定システムは、測定対象とする三次元物体の形状を測定する。
図７に示す例では、三次元形状測定システムは、測定対象とする三次元物体の形状を測定する。
格子投影の制御部１１から設定された(その基準面でピッチ寸法および高さ係数(その規則格子の位相と高さの変換係数)が知られた)格子模様が格子投影装置４より投影され、投影されたパターンを撮像素子３(各種カメラや顕微鏡)用いて縞格子画像の記録部１２より物体の形状に応じてゆがんだ格子模様を撮影する。
撮影された画像に対して、本発明の方法による２段モアレ縞の生成部と位相解析の演算部１３で２段モアレ縞画像の位相分布を算出し、三次元形状の演算部１４で位相分布からさらに高さ分布を算出する。
その計測結果をモニター表示部２で表示される。
その計測結果をモニター表示部２で表示される。
(第２の実施形態：変位・ひずみ計測)
図８は、本実施形態に係る材料の変位・ひずみ分布測定システムにおける計算機１、撮像素子３と測定対象とする構造材料７、微細格子８(１次元または２次元)の配置例を示す図である。
図８は、本実施形態に係る材料の変位・ひずみ分布測定システムにおける計算機１、撮像素子３と測定対象とする構造材料７、微細格子８(１次元または２次元)の配置例を示す図である。
図８に示す例では、変形測定システムは、測定対象とする構造材料表面に付与されそのピッチ寸法が知られた微細格子を撮像素子３より格子画像またはモアレ縞を記録する。
付与されとは、規則格子を材料表面に貼り付けるだけでなく、例えば金属材料で観察される原子配列画像に見られるように試料表面にあらかじめ存在する繰り返し模様であってもよい。
記録部１２で記録されたモアレ縞または格子画像から生成されたモアレ縞を、２段モアレ縞の生成部と位相解析の演算部１３で各処理し、さらに変位・ひずみの演算部１５で計測した変位・ひずみ分布結果がモニター２で表示される。
付与されとは、規則格子を材料表面に貼り付けるだけでなく、例えば金属材料で観察される原子配列画像に見られるように試料表面にあらかじめ存在する繰り返し模様であってもよい。
記録部１２で記録されたモアレ縞または格子画像から生成されたモアレ縞を、２段モアレ縞の生成部と位相解析の演算部１３で各処理し、さらに変位・ひずみの演算部１５で計測した変位・ひずみ分布結果がモニター２で表示される。
形状変形計測手順
図９に２段モアレ法による形状・変形計測のフローチャート示す。
図９に２段モアレ法による形状・変形計測のフローチャート示す。
(Ｓ１)三次元形状計測の場合は物体表面に格子パターンを投影する。
変位・ひずみ分布計測の場合は試料表面上に規則格子を作製する、または試料自体の周期構造を格子としてもよい。
変位・ひずみ分布計測の場合は試料表面上に規則格子を作製する、または試料自体の周期構造を格子としてもよい。
(Ｓ２)次に画像記録装置の下に試料格子を配置する。
画像記録装置はレーザー走査型顕微鏡(LSM)、走査型電子顕微鏡(SEM)、透過型電子顕微鏡(TEM)や原子間力顕微鏡(AFM)などの走査型顕微鏡、電荷結合素子(CCD)カメラまたは相補的金属酸化物半導体(CMOS)カメラなどのデジタルカメラ、またはモアレ干渉計であってよい。
次に顕微鏡の走査線間隔またはカメラのピクセルサイズを試料格子ピッチ(またはその整数(2以上の自然数)倍やその整数を分母とする分数倍)の近傍に調整する。
試料格子と顕微鏡走査線またはピクセル配列との間の干渉により広視野CCD/CMOSモアレ縞が画像記録装置に記録される。
画像記録装置はレーザー走査型顕微鏡(LSM)、走査型電子顕微鏡(SEM)、透過型電子顕微鏡(TEM)や原子間力顕微鏡(AFM)などの走査型顕微鏡、電荷結合素子(CCD)カメラまたは相補的金属酸化物半導体(CMOS)カメラなどのデジタルカメラ、またはモアレ干渉計であってよい。
次に顕微鏡の走査線間隔またはカメラのピクセルサイズを試料格子ピッチ(またはその整数(2以上の自然数)倍やその整数を分母とする分数倍)の近傍に調整する。
試料格子と顕微鏡走査線またはピクセル配列との間の干渉により広視野CCD/CMOSモアレ縞が画像記録装置に記録される。
または２つのコヒーレントレーザー干渉によって引き起こされる干渉モアレ縞が干渉計によって記録される。
顕微鏡モアレ、CCD/CMOSモアレ、干渉モアレは、モアレ縞を直接観察、記録する手法である。
顕微鏡モアレ、CCD/CMOSモアレ、干渉モアレは、モアレ縞を直接観察、記録する手法である。
これ以外に格子画像を撮影して画像処理によりモアレ縞(たとえばデジタルモアレ、サンプリングモアレ、オーバーラップモアレなど)を生成する方法がある。
これらの手法においては視野を拡大するために、撮影画像上の試料格子ピッチは2ピクセルまで小さくできる。
これらの手法においては視野を拡大するために、撮影画像上の試料格子ピッチは2ピクセルまで小さくできる。
サンプリングモアレ法やデジタルモアレ法から得られるモアレ縞はダウンサンプリングと輝度補間、または試料格子と試料格子のピッチ(または整数倍や分数倍数)に近いデジタル格子との干渉によって生成される。オーバーラップ(幾何学的)モアレ縞は、変形前後の試料格子間の重ね合わせによって生成される。
(Ｓ３)これら従来のモアレ縞を使用して、その間隔に近いピッチによりダウンサンプリング、輝度補間、そして空間位相シフト法を用いて、位相がシフトされた２段モアレ縞を生成する。
それから位相がシフトされた２段モアレ縞から２段モアレ縞の位相分布を求めることができる。
それから位相がシフトされた２段モアレ縞から２段モアレ縞の位相分布を求めることができる。
(Ｓ４)変形前後の２段モアレ縞の位相差から(Ｓ５)変位・ひずみ分布を広視野で正確に定量評価することができる。
この形状変形計測法を２段モアレ法と命名する。
２段モアレ法は、モアレ縞からデジタル格子を利用してその位相差が同一となるモアレ縞を生成し得る限り多段階的に拡張することが可能である。
よって２段モアレ法には本発明の課題に適合する限度において多段モアレ法が含まれる。
２段モアレ法は、モアレ縞からデジタル格子を利用してその位相差が同一となるモアレ縞を生成し得る限り多段階的に拡張することが可能である。
よって２段モアレ法には本発明の課題に適合する限度において多段モアレ法が含まれる。
また、一枚の変形後のモアレ縞のみから変形していない部分のモアレ縞の位相変化(位相勾配)と変形している部分のモアレ縞位相変化とを比較することで変形前のモアレ縞を推定できる場合には、一枚の変形後のモアレ縞からひずみを求めることもできる。
モアレ縞の位相変化を求める際には推定されたモアレ縞と変形後のモアレ縞とを変形前後のモアレ縞として２段モアレ法を適用してもよい。
ただしこの手法は、図７に示した第１の実施形態で説明した三次元形状計測において格子を斜めから物体に投影するため基準面におけるゆがんだ格子模様が撮影されるような場合にはそのゆがんだ格子画像から基準面での正規な格子画像(基準面であったなら得られたであろう)が推定される場合を除いて不向きである。
モアレ縞の位相変化を求める際には推定されたモアレ縞と変形後のモアレ縞とを変形前後のモアレ縞として２段モアレ法を適用してもよい。
ただしこの手法は、図７に示した第１の実施形態で説明した三次元形状計測において格子を斜めから物体に投影するため基準面におけるゆがんだ格子模様が撮影されるような場合にはそのゆがんだ格子画像から基準面での正規な格子画像(基準面であったなら得られたであろう)が推定される場合を除いて不向きである。
本発明のプログラムの流れ
次に本発明に係る形状変形測定処理について説明する。図１０は、本発明に係る形状変形測定処理を示すフローチャートである。図１０に示す処理は、ステップS101、S102、S103、S104とS105を有する。
図１０に示す処理では、ステップS101の後、順番にステップS105の処理を行う。
プログラム言語は上述の数式が演算可能であれば適宜選択できそのプログラムを実行する計算機は専用機・汎用機のいずれでもよい。
格子模様またはモアレ縞を記録した(S101)後、ピッチ寸法、高さ係数、モアレ縞を生成し得る整数等パラメータを入力して(S102)、２段モアレ法を用いて必要な演算行い(S103)、２段モアレ縞の位相分より形状・変位・ひずみ量を算出して(S104)、その結果を表示する(S105)。
次に本発明に係る形状変形測定処理について説明する。図１０は、本発明に係る形状変形測定処理を示すフローチャートである。図１０に示す処理は、ステップS101、S102、S103、S104とS105を有する。
図１０に示す処理では、ステップS101の後、順番にステップS105の処理を行う。
プログラム言語は上述の数式が演算可能であれば適宜選択できそのプログラムを実行する計算機は専用機・汎用機のいずれでもよい。
格子模様またはモアレ縞を記録した(S101)後、ピッチ寸法、高さ係数、モアレ縞を生成し得る整数等パラメータを入力して(S102)、２段モアレ法を用いて必要な演算行い(S103)、２段モアレ縞の位相分より形状・変位・ひずみ量を算出して(S104)、その結果を表示する(S105)。
本発明は、２段モアレ法による広視野で高い測定精度でまたは高精度で動的に全視野変位とひずみを測定できるため、以下の効果がある(図１１)。
効果１：提案手法は、従前の顕微鏡、CCD、デジタル、オーバーラップモアレ法とモアレ干渉法に使用される縞中央化法と比較して、位相解析により顕著に測定精度が向上する。ひずみ感度は、理論的には10-6 に達する。また本発明によれば、各種顕微鏡の走査歪みに起因したひずみ測定誤差を除去できる。
効果２：提案手法の解析領域は、従来のサンプリング/デジタル/オーバーラップモアレ法と比して、数十から数百倍大きくできる。
効果３：変形計測の記録に必要なのは単一のモアレパターンまたは単一格子イメージのみである。
提案手法は、時空間の位相シフト法を適用できない動的な材料試験に適している。
提案手法は、時空間の位相シフト法を適用できない動的な材料試験に適している。
効果４：この技法は、良好なノイズ特性があって、モアレパターンまたは格子イメージを記録できる場合、ナノスケールからメータースケールの非破壊変形計測に使用することができる。
シミュレーションによる任意の繰り返し模様の２段モアレ縞の生成原理の確認
ここでは、任意の２次元繰り返し模様を対象に、２段モアレ法を適用した場合の結果を示す。
図１２(a)は物体表面に存在するまたは付与した繰り返し模様(英語A)とする。
図１２(a)の模様に対して、x方向およびy方向に一定ピッチで１回目のダウンサンプリングによって得られるモアレ縞画像は図１２(b)である。
ここでは、任意の２次元繰り返し模様を対象に、２段モアレ法を適用した場合の結果を示す。
図１２(a)は物体表面に存在するまたは付与した繰り返し模様(英語A)とする。
図１２(a)の模様に対して、x方向およびy方向に一定ピッチで１回目のダウンサンプリングによって得られるモアレ縞画像は図１２(b)である。
このモアレ縞画像を格子とみなし、さらにx方向およびy方向に一定ピッチで２回目のダウンサンプリングによって得られるモアレ縞画像は図１２(c)である。
元の繰り返し模様であるAを拡大して、さらに拡大していることがわかる。
これを２段モアレと呼ぶ。
元の繰り返し模様であるAを拡大して、さらに拡大していることがわかる。
これを２段モアレと呼ぶ。
シミュレーションによる倍増型モアレ縞の生成原理の確認
本発明では、格子画像をダウンサンプリング処理して得られるモアレ縞の間隔は疎であることがある。
このような場合は、本発明である第３の方法により、適切なモアレ縞の間隔になるように、モアレ縞の間隔を密となるような倍増型モアレ縞を再構成することができる。
図１３にシミュレーションによる(倍増型)モアレ縞の生成結果を示す。
本発明では、格子画像をダウンサンプリング処理して得られるモアレ縞の間隔は疎であることがある。
このような場合は、本発明である第３の方法により、適切なモアレ縞の間隔になるように、モアレ縞の間隔を密となるような倍増型モアレ縞を再構成することができる。
図１３にシミュレーションによる(倍増型)モアレ縞の生成結果を示す。
図１３(a)はシミュレーションにより生成されたピッチ間隔が3.02画素の格子パターンである。
図１３(b)は図１３(a)の格子画像をダウンサンプリングピッチ２画素で間引き処理して得られた２枚の位相がシフトされたモアレ縞画像である。
格子ピッチP(3.02画素)とサンプリングピッチT(2画素)が大きく異なっているため、モアレ縞の生成ができていないことが確認できる。
格子ピッチP(3.02画素)とサンプリングピッチT(2画素)が大きく異なっているため、モアレ縞の生成ができていないことが確認できる。
図１３(c)は図１３(a)の格子画像をダウンサンプリングピッチ３画素で間引き処理して得られた３枚の位相がシフトされたモアレ縞画像である。
格子ピッチP(3.02画素)とサンプリングピッチT(３画素)が近いため、きれいなモアレ縞を生成できるが、モアレ縞の間隔が疎であることがわかる。
このまま２段モアレ法を適用すると、空間分解能が著しく低下するという問題点がある。
格子ピッチP(3.02画素)とサンプリングピッチT(３画素)が近いため、きれいなモアレ縞を生成できるが、モアレ縞の間隔が疎であることがわかる。
このまま２段モアレ法を適用すると、空間分解能が著しく低下するという問題点がある。
図１３(d)は図１３(a)の格子画像をダウンサンプリングピッチ４画素で間引き処理して得られた４枚の位相がシフトされた間引き数が４画素のときに生成されるモアレ縞画像である。
図１３(b)と同様に、格子ピッチP(3.02画素)とサンプリングピッチT(4画素)が大きく異なっているため、モアレ縞の生成ができていないことがわかる。
図１３(b)と同様に、格子ピッチP(3.02画素)とサンプリングピッチT(4画素)が大きく異なっているため、モアレ縞の生成ができていないことがわかる。
図１３(e)はサンプリングピッチが３画素のときに生成されたモアレ縞画像(図１３(c))を用いて算出したモアレ縞の位相分布画像である。
この位相分布画像を用いて、調整係数Nを１から２０まで変えた場合に再構成された(倍増型)モアレ縞画像の結果を図１３(f)に示す。
このように、調整係数Nを変えることで、モアレ縞の間隔を疎から密に調整できるようになる。
このように、調整係数Nを変えることで、モアレ縞の間隔を疎から密に調整できるようになる。
格子投影法による物体の三次元形状の非接触計測の実験検証
本発明で提案する２段モアレ法による物体の三次元形状の非接触計測の精度向上を実験よりその効果を確認した。
本発明で提案する２段モアレ法による物体の三次元形状の非接触計測の精度向上を実験よりその効果を確認した。
図１４は、図７に示すような３次元形状測定装置を用いて、格子投影装置より投影された格子をカメラで撮影した格子画像を２段モアレ法に適用した場合の解析結果例を示す。
ここでは基準面に一定周期(プロジェクタでは8画素周期)の正弦波を投影して、撮影した格子画像を図１４(a)に示す。
ここでは基準面に一定周期(プロジェクタでは8画素周期)の正弦波を投影して、撮影した格子画像を図１４(a)に示す。
図１４(a)に対して、従来のサンプリングモアレ法を適用して、横方向に25画素のダウンサンプリングと輝度補間をして、位相がシフトされた25枚のモアレ縞を用いて計算したモアレ縞の位相分布を図１４(b)に示す。
これらのモアレ縞画像を位相がシフトされた格子画像を見なし、さらに横方向に100画素のダウンサンプリングと輝度補間を行い、生成される100×25の2500枚の位相シフトされた２段モアレ縞画像を利用して、２次元離散的フーリエ変換より算出した２段モアレ縞の位相分布を図１４(c)に示す。
２段モアレ縞の位相分布にサンプリング点の位相を加える(式(7))ことで得られる元の投影格子の位相分布を図１４(d)に示す。
この実施例の場合、１枚の格子画像から2500枚の位相がシフトされたモアレ縞を生成でき、より高精度な位相解析が行えるようになる。
この実施例の場合、１枚の格子画像から2500枚の位相がシフトされたモアレ縞を生成でき、より高精度な位相解析が行えるようになる。
図１５は図１４(a)に示すような基準面の格子画像を従来技術である位相シフト法、サンプリングモアレ法、時空位相シフト法および本発明による解析結果の比較結果である。
図１５(a)は撮影された8枚の位相シフトされた格子画像を従来の位相シフト法より算出した位相分布を横方向に空間微分した(位相勾配)結果である。
左図は位相勾配分布図であり、右図は横中央１ラインの断面図である。
図１５(a)は撮影された8枚の位相シフトされた格子画像を従来の位相シフト法より算出した位相分布を横方向に空間微分した(位相勾配)結果である。
左図は位相勾配分布図であり、右図は横中央１ラインの断面図である。
同様に図１５(b)は撮影された１枚の格子画像を従来のサンプリングモアレ法(間引き数は２５画素)より算出した結果である。
図１５(c)は撮影された８枚の位相シフトされた格子画像を従来の時空位相シフト法(間引き数は25画素)より算出した結果である。
図１５(d)は撮影された１枚の格子画像を本発明(間引き数は１段目では25画素、２段目では100画素)より算出した結果である。
基準面の位相勾配はある傾きをもったなめらかな直線であることが正しい結果である。
基準面の位相勾配はある傾きをもったなめらかな直線であることが正しい結果である。
図１５の結果からわかるように、従来の位相シフト法は測定誤差が大きい。
サンプリングモアレ法は１枚の格子画像で位相算出できるが、測定誤差がある。
時空位相シフト法は時間と空間の両方の輝度情報を使用するため、位相シフト法とサンプリングモアレ法よりも精度が高いが、複数枚の画像が必要であるため、動的な計測に向かない。
サンプリングモアレ法は１枚の格子画像で位相算出できるが、測定誤差がある。
時空位相シフト法は時間と空間の両方の輝度情報を使用するため、位相シフト法とサンプリングモアレ法よりも精度が高いが、複数枚の画像が必要であるため、動的な計測に向かない。
一方、本発明によれば、１枚の格子画像のみを用いて、最も精度の良い結果が得られていることから、本発明の有効性が示された。
次に、本発明である第２の方法を利用した高精度で動的な３次元形状計測の実施例を示す。
図１６は、図７に示すような３次元形状測定装置を用いて、格子投影装置より投影された格子をカメラより撮影した基準面と物体(手)の格子画像を２段モアレ法に適用し、基準面と物体の位相差から物体の高さを求めた解析結果である。
図１６は、図７に示すような３次元形状測定装置を用いて、格子投影装置より投影された格子をカメラより撮影した基準面と物体(手)の格子画像を２段モアレ法に適用し、基準面と物体の位相差から物体の高さを求めた解析結果である。
図１６(a)は上から順にそれぞれ基準面と物体に格子ピッチ比が８と９の２種類の正弦波状格子パターンを投影して異なる角度に設置したカメラで撮影した画像である。
これらの画像に対して、それぞれ２段モアレ法より算出した位相分布を図１６(b)に示す。
ここでは、奥行きのある物体でも高さを測定できるように、格子ピッチ比が８と９の２種類の格子パターンを用いている。
これらの画像に対して、それぞれ２段モアレ法より算出した位相分布を図１６(b)に示す。
ここでは、奥行きのある物体でも高さを測定できるように、格子ピッチ比が８と９の２種類の格子パターンを用いている。
格子ピッチ比が８における基準面と物体の位相差を図１６(c)に示し、同様に格子ピッチ比が９における基準面と物体の位相差を図１６(d)に示す。
格子ピッチ比が８の格子の場合の解析条件は図１４と同じである。
格子ピッチ比が９の格子の場合の解析条件は１回目のダウンサンプリングピッチは29画素で、２回目のダウンサンプリングピッチは110画素である。
格子ピッチ比が８の格子の場合の解析条件は図１４と同じである。
格子ピッチ比が９の格子の場合の解析条件は１回目のダウンサンプリングピッチは29画素で、２回目のダウンサンプリングピッチは110画素である。
図１６(c)と図１６(d)の２種類の位相差結果を用いて、既存の複数ピッチの格子による位相接続のアルゴリスムに適用して、得られた物体(手)の位相接続後の位相分布が図１６(e)である。
図１６(e)にキャリブレーションで得られる高さ係数をかければ最終的に物体の高さ分布を得ることができる。
図１７は、従来の位相シフト法と、サンプリングモアレ法、および本発明である２段モアレ法でそれぞれ得られた手の高さ測定結果である。
図１７(a)は、従来の位相シフト法による測定結果で、格子ピッチ比が８の格子画像を2π/8ずつ位相がシフトされた８枚の画像と、格子ピッチ比が９の格子画像を2π/9ずつ位相がシフトされた９枚の画像の計１７枚の画像を用いて求めた手の形状データである。
左側は形状分布であり、右側は横１ラインのデータを示している。
左側は形状分布であり、右側は横１ラインのデータを示している。
図１７(b)は、従来のサンプリングモアレ法による測定結果で、格子ピッチ比８と９のそれぞれ１枚で計２枚の画像のみ用いて動的な形状計測を行えるが、測定結果にノイズが多いのが問題点である。
図１７(c)は、本発明である２段モアレ法で得られた結果である。
この実施例では、投影された格子パターンのSN(Signal-noise)比が極めて悪い条件での解析結果である。
したがって、従来の位相シフト法(図１７(a))や時空位相シフト法は多くの位相シフト格子画像を必要とするため、動的計測への適用が困難である。
一方でサンプリングモアレ法(図１７(b))は動的計測に向いているが、測定結果にノイズが多く、うまく形状を測定できない箇所が多く存在する。
この実施例では、投影された格子パターンのSN(Signal-noise)比が極めて悪い条件での解析結果である。
したがって、従来の位相シフト法(図１７(a))や時空位相シフト法は多くの位相シフト格子画像を必要とするため、動的計測への適用が困難である。
一方でサンプリングモアレ法(図１７(b))は動的計測に向いているが、測定結果にノイズが多く、うまく形状を測定できない箇所が多く存在する。
対して本発明(図１７(c))によれば、格子ピッチ比８と９のそれぞれ１枚で計２枚の画像のみ用いて測定誤差の少ない手の形状を算出できることを確認できた。
図１７の実験結果からわかるように、本発明は高精度でかつ動的な三次元形状計測に有効である。
ただし、２段モアレ法を用いる本手法は、広い空間の輝度情報を用いるため、急激に変化する物体の計測に向かないことに留意する必要がある。
図１７の実験結果からわかるように、本発明は高精度でかつ動的な三次元形状計測に有効である。
ただし、２段モアレ法を用いる本手法は、広い空間の輝度情報を用いるため、急激に変化する物体の計測に向かないことに留意する必要がある。
均一な変形計測において提案手法のシミュレーション検証
実施例４は提案する２段モアレ法をシミュレーションで実施し、均一な変形計測した場合の正確性を検証する例である。
実施例４は提案する２段モアレ法をシミュレーションで実施し、均一な変形計測した場合の正確性を検証する例である。
10ピクセルピッチの規則格子は線型に伸張することが可能で現在のひずみは100με(図１８(a))。規則格子の格子像の長さは10000ピクセルだった。
ここでいうマイクロンストレイン(με)はひずみ量の単位で10のマイナス6乗である。
ここでいうマイクロンストレイン(με)はひずみ量の単位で10のマイナス6乗である。
参照格子のピッチが９ピクセルの時に変形前後の公知のディジタルサンプリングモアレ縞が得られた(図１８(b))。
変形前後の公知のモアレ縞の間隔はいずれも約90ピクセルだった。
変形前後の公知のモアレ縞の間隔はいずれも約90ピクセルだった。
公知のモアレ縞を規則格子として扱い88ピクセルでダウンサンプリングを行い、空間位相シフト法と結合して変形前後の２段モアレ縞を求めた(図１８(c))。
本発明を用いて変形前後の２段モアレ縞の位相を計算した(図１８(d))。
次に同一の分布特性を有する位相の変化、または変位を測定した(図１８(e))。
次に同一の分布特性を有する位相の変化、または変位を測定した(図１８(e))。
最後にひずみの分布が決定されて、平均的ひずみ101μεが得られた。
プリセットひずみ100μεに比較すると計測誤差は僅か1μεであり、本測定方法の正確性が確認できた。
プリセットひずみ100μεに比較すると計測誤差は僅か1μεであり、本測定方法の正確性が確認できた。
不均一変形計測において提案手法のシミュレーション検証
本実施例は不均一な変位とひずみを測定された分布特性(図１９)が理論的特性と同一であることを示す例である。
１ピッチが8.1ピクセルの規則格子(図１９(a)左、変形前)にMATLAB(登録商標)関数peaks(x,y)から得た数値をその位相に加えて変形した(図１９(a)右、変形後)。
格子像は256×180ピクセルだった。
１ピッチが8.1ピクセルの規則格子(図１９(a)左、変形前)にMATLAB(登録商標)関数peaks(x,y)から得た数値をその位相に加えて変形した(図１９(a)右、変形後)。
格子像は256×180ピクセルだった。
参照格子を７ピクセルとすると変形前後の公知のモアレ縞が現れた(図１９(b))。
変形前後の公知のモアレ縞の間隔は46〜57ピクセルだった。
変形前後の公知のモアレ縞の間隔は46〜57ピクセルだった。
その公知のモアレ縞を41ピクセルでダウンサンプルして空間位相シフト法と結合して変形前後の位相シフトした２段モアレ縞を求めた(図１９(c))。
本発明を用いて変形前後のモアレ縞の位相を求めた(図１９(d))。
位相の変化(図１９(e))、変位とひずみ分布を測定した。
測定された変位とひずみ(図１９(f))の分布特性は理論的に求めた結果と合致した。
位相の変化(図１９(e))、変位とひずみ分布を測定した。
測定された変位とひずみ(図１９(f))の分布特性は理論的に求めた結果と合致した。
ランダムノイズのある傾斜格子の提案手法のシミュレーション検証
本実施例ではランダムノイズを有する傾斜格子の変位測定を説明する。
格子線の傾斜角は0.0086radで垂直方向の格子ピッチは2.0572ピクセル、格子像は1300×900ピクセルだった。
傾斜格子に格子の輝度振幅の10％のランダムノイズを追加し、50με、100με、500με、1000με、2000με、3000μεに相当する引張りひずみを負荷した場合の画像をシミュレーションにより生成した。
本実施例ではランダムノイズを有する傾斜格子の変位測定を説明する。
格子線の傾斜角は0.0086radで垂直方向の格子ピッチは2.0572ピクセル、格子像は1300×900ピクセルだった。
傾斜格子に格子の輝度振幅の10％のランダムノイズを追加し、50με、100με、500με、1000με、2000με、3000μεに相当する引張りひずみを負荷した場合の画像をシミュレーションにより生成した。
ひずみが0、および1000μεの場合のモアレ縞、および２段モアレ縞生成プロセスを図２０に示す。２段モアレ縞の位相差から、変位およびひずみを測定できる。
シミュレーションで与えたひずみを横軸にして、計測された平均ひずみの誤差を図２１(a)に示す。
図２１(b)は測定されたひずみの標準偏差を示している。
すべての相対誤差は1%未満であり、標準偏差は4×10−5未満であることから、高精度なひずみ計測ができていることが確認できた。
すべての相対誤差は1%未満であり、標準偏差は4×10−5未満であることから、高精度なひずみ計測ができていることが確認できた。
載荷したアルミニウムの変形の実験的検証
本実施例では提案する変形計測の妥当性と正確性を実験により示す。
アルミニウム試料は、レーザー顕微鏡(図２２)の下で 引張されていた。
アルミニウム試料の寸法は長さ・幅・厚みはそれぞれ27、6.3、0.5(mm)であった。
UVナノインプリントリソグラフィーによりアルミニウム試料の表面に3.0μmピッチの格子を作製した。
本実施例では提案する変形計測の妥当性と正確性を実験により示す。
アルミニウム試料は、レーザー顕微鏡(図２２)の下で 引張されていた。
アルミニウム試料の寸法は長さ・幅・厚みはそれぞれ27、6.3、0.5(mm)であった。
UVナノインプリントリソグラフィーによりアルミニウム試料の表面に3.0μmピッチの格子を作製した。
試料の反対側にひずみゲージを貼り付けたのち、顕微鏡の下で引張負荷を与えた。
対物レンズの倍率を５倍にして試料上の格子と走査線との干渉で生じたレーザー走査モアレ縞を観測した。
引張過程中の一連のレーザー走査モアレ縞を記録した。
対物レンズの倍率を５倍にして試料上の格子と走査線との干渉で生じたレーザー走査モアレ縞を観測した。
引張過程中の一連のレーザー走査モアレ縞を記録した。
ひずみゲージの値が3000μεの場合を例にして、引張負荷下のアルミニウム試料の変形計測プロセスを図２３に示す。
提案手法から評価された平均ひずみとひずみゲージから評価されたひずみを用いて得られたアルミニウム試料の応力-ひずみ曲線を、図２４に示す。
提案手法から評価された平均ひずみ値はひずみゲージ法のそれと良く一致し、提案手法の妥当性と測定精度を検証することができた。
提案手法から評価された平均ひずみ値はひずみゲージ法のそれと良く一致し、提案手法の妥当性と測定精度を検証することができた。
提案された二次モアレ法とそのソフトウェアは、ナノスケールからメートルスケールに至る様々な素材や構造の全視野変位・ひずみ測定に使用できる。
解析可能な対象には、航空宇宙、自動車、電子パッケージ、生物医学、軍事、材料製造などの産業分野における金属、高分子、セラミックス、半導体、複合材料、ハイブリッド構造、薄膜がある。
産業分野における代表的なアプリケーションとして、主に次の４つが挙げられる。
産業分野における代表的なアプリケーションとして、主に次の４つが挙げられる。
１)はく離、亀裂伝播と不安定性モードの評価
提案手法は全視野変形から応力集中領域を見つけることによってき裂の発生領域を予測することができる。
提案手法は全視野変形から応力集中領域を見つけることによってき裂の発生領域を予測することができる。
き裂伝播形態はき裂先端近傍の変形分布を測定することによって評価できる。
様々な機械的負荷、電気負荷、熱負荷、磁気負荷、およびそれらが結合した混合負荷による不安定性モードと層間はく離、座屈、き裂を含む損傷形態を定量的に解析できる。
様々な機械的負荷、電気負荷、熱負荷、磁気負荷、およびそれらが結合した混合負荷による不安定性モードと層間はく離、座屈、き裂を含む損傷形態を定量的に解析できる。
２)残留ひずみ・応力分布の評価、最適化設計指針の提供
材料内部の残留ひずみとその構造的安定性への影響を応力解放法(加熱、穿孔法、リングコア法など) と組み合わせることによって検出できる。
材料内部の残留ひずみとその構造的安定性への影響を応力解放法(加熱、穿孔法、リングコア法など) と組み合わせることによって検出できる。
界面近傍の変形分布計測は、破壊し易い領域を見つける助けになり、材料の強化・強靭化および材料および界面の最適設計の指針を提供する。
３)機械的特性評価と構造健全性監視
変形計測情報は、応力-ひずみ曲線やヤング率、ポアソン比、弾性限界、降伏強度、極限強度などの材料定数を決定するのに使用できる。
機械的、電気的負荷または熱負荷下がもたらす変位とひずみを評価することができ、構造健全性を監視する助けとなる。
変形計測情報は、応力-ひずみ曲線やヤング率、ポアソン比、弾性限界、降伏強度、極限強度などの材料定数を決定するのに使用できる。
機械的、電気的負荷または熱負荷下がもたらす変位とひずみを評価することができ、構造健全性を監視する助けとなる。
４)高さ深さ平坦性の評価と製造品質の制御
測定した３次元形状から、対象物の高さ、深さ、表面平坦性と面外変位を定量的に評価することができる。製造品質の制御、機械視覚、自動処理などの工業分野に有用である。
測定した３次元形状から、対象物の高さ、深さ、表面平坦性と面外変位を定量的に評価することができる。製造品質の制御、機械視覚、自動処理などの工業分野に有用である。
１ 計算機
２ モニター
３ 撮像素子(各種カメラや顕微鏡)
４ 格子投影装置
５ 測定対象とする三次元物体
６ 基準面(平面物体)
７ 測定対象とする構造材料
８ 微細格子(１次元または２次元)
１１ 格子投影の制御部
１２ 縞格子画像の記録部
１３ ２段モアレ縞の生成部と位相解析の演算部
１４ 三次元形状の演算部
１５ 変位・ひずみの演算部
２ モニター
３ 撮像素子(各種カメラや顕微鏡)
４ 格子投影装置
５ 測定対象とする三次元物体
６ 基準面(平面物体)
７ 測定対象とする構造材料
８ 微細格子(１次元または２次元)
１１ 格子投影の制御部
１２ 縞格子画像の記録部
１３ ２段モアレ縞の生成部と位相解析の演算部
１４ 三次元形状の演算部
１５ 変位・ひずみの演算部
Claims (8)
- 変位分布を測定する測定方法であって、
試料の表面に表れる格子を示す格子画像を所定の方向に第１のダウンサンプリングおよび第１の輝度補間を行って第１のモアレ縞を取得する第１のステップと、
前記第１のモアレ縞を前記方向にさらに第２のダウンサンプリングおよび第２の輝度の補間を行って第２のモアレ縞を生成する第２のステップと、
変位前後の前記第２のモアレ縞の位相分布の差から前記所定の方向または奥行き方向への変位分布を演算する第３のステップと、
を有する測定方法。 - 前記第１のステップは、
前記第１のダウンサンプリングにおける第１のダウンサンプリングピッチＴの逆数のｔ（ｔは、０以上Ｔ−１の整数）倍でそれぞれ前記方向に位相をシフトしたＴ枚の前記第１のモアレ縞を生成するステップを有し、
前記第２のステップは、
Ｔ枚の前記第１のモアレ縞のそれぞれについて、
前記第２のダウンサンプリングにおける第２のダウンサンプリングピッチＴ （２） の逆数のｋ（ｋは、０以上Ｔ （２） −１の整数）倍でそれぞれ前記方向に位相をシフトしたＴ （２） 枚の前記第２のモアレ縞を生成するステップと、
Ｔ・Ｔ （２） 枚の前記第２のモアレ縞をフーリエ変換して前記第２のモアレ縞の位相分布を算出するステップと、
を有する請求項１に記載の測定方法。 - 前記第１のモアレ縞の位相に調整係数Ｎ（Ｎは、２以上の自然数）を乗じて前記第１のモアレ縞のピッチを調整するステップを有し、
ピッチを調整した前記第１のモアレ縞を前記第２のステップに適用する
請求項１または請求項２に記載の測定方法。 - 前記第１のステップは、
前記格子画像を所定の第１のダウンサンプリングピッチで走査して前記第１のモアレ縞を記録する
請求項１または請求項３に記載の測定方法。 - 前記格子は前記試料に付与された格子であって、
前記変位分布を前記方向に微分してひずみ分布を算出する第４のステップ
を有する請求項１から請求項４のいずれか一項に記載の測定方法。 - 前記位相分布の差は、前記試料を設置する基準面に投影装置により投影された前記格子に基づく前記第２のモアレ縞の位相分布と、前記投影装置により前記試料に投影された前記格子に基づく前記第２のモアレ縞の位相分布との差であり、
前記第３のステップは、
前記奥行き方向への変位と位相差との関係を示す換算係数に基づいて、前記位相分布の差から前記試料の形状として前記奥行き方向への変位分布を算出するステップ
を有する請求項１から請求項３のいずれか一項に記載の測定方法。 - 試料の表面に表れる格子を示す格子画像を所定の方向に第１のダウンサンプリングおよび第１の輝度補間を行って第１のモアレ縞を記録する記録部と、
前記第１のモアレ縞を前記方向にさらに第２のダウンサンプリングおよび第２の輝度の補間を行って第２のモアレ縞を生成する生成部と、
変位前後の前記第２のモアレ縞の位相分布の差から前記所定の方向もしくは奥行き方向への変位分布を演算する演算部と、
を有する測定装置。 - 測定装置のコンピュータに、
試料の表面に表れる格子を示す格子画像を所定の方向に第１のダウンサンプリングおよび第１の輝度補間を行って第１のモアレ縞を取得する第１のステップと、
前記第１のモアレ縞を前記方向にさらに第２のダウンサンプリングおよび第２の輝度の補間を行って第２のモアレ縞を生成する第２のステップと、
変位前後の前記第２のモアレ縞の位相分布の差から前記所定の方向または奥行き方向への変位分布を演算する第３のステップと、
を実行させるためのプログラム。
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2016188217 | 2016-09-27 | ||
| JP2016188217 | 2016-09-27 | ||
| PCT/JP2017/020435 WO2018061321A1 (ja) | 2016-09-27 | 2017-06-01 | 周期模様を利用した三次元形状・変位・ひずみ測定装置、方法およびそのプログラム |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPWO2018061321A1 JPWO2018061321A1 (ja) | 2019-04-18 |
| JP6583761B2 true JP6583761B2 (ja) | 2019-10-02 |
Family
ID=61760268
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2018541896A Active JP6583761B2 (ja) | 2016-09-27 | 2017-06-01 | 周期模様を利用した三次元形状・変位・ひずみ測定装置、方法およびそのプログラム |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US10655954B2 (ja) |
| EP (1) | EP3489621A4 (ja) |
| JP (1) | JP6583761B2 (ja) |
| WO (1) | WO2018061321A1 (ja) |
Families Citing this family (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2018156702A1 (en) * | 2017-02-23 | 2018-08-30 | Nikon Corporation | Measurement of a change in a geometrical characteristic and/or position of a workpiece |
| JP6718160B2 (ja) * | 2017-02-23 | 2020-07-08 | 国立研究開発法人産業技術総合研究所 | 残留熱ひずみ測定方法、残留熱ひずみ測定装置、及びそのプログラム |
| JP6897398B2 (ja) * | 2017-07-31 | 2021-06-30 | セイコーエプソン株式会社 | 三次元形状測定装置、ロボットシステム、及び三次元形状測定方法 |
| JP6982336B2 (ja) * | 2017-10-24 | 2021-12-17 | 国立研究開発法人産業技術総合研究所 | 計測方法、コンピュータプログラム及び計測システム |
| US11226195B1 (en) * | 2018-06-27 | 2022-01-18 | United States Of America Represented By The Secretary Of The Air Force | Method and system for measuring strain in a 3D printed part |
| CN109307475B (zh) * | 2018-09-04 | 2020-01-14 | 西安交通大学 | 差频主动扫描式光栅位移传感器及测量方法 |
| CN109084698B (zh) * | 2018-10-25 | 2024-05-24 | 中国测绘科学研究院 | 一种用于工程建筑物形变监测的方法及*** |
| JP7276016B2 (ja) * | 2019-09-04 | 2023-05-18 | 株式会社Ihi | 非接触計測装置 |
| CN110987677B (zh) * | 2019-12-25 | 2021-08-10 | 上海交通大学 | 提高微形变加载精度的方法及装置 |
| CN111476788B (zh) * | 2020-04-27 | 2023-08-25 | 武汉精立电子技术有限公司 | 一种显示屏夹层缺陷检测方法及*** |
| CN112414304B (zh) * | 2020-11-18 | 2022-06-21 | 天津科技大学 | 基于激光光栅投射的焊后焊缝表面三维测量方法 |
| CN113048901B (zh) * | 2021-03-05 | 2023-06-27 | 中国建筑材料科学研究总院有限公司 | 基于光流算法测量显微镜下纳米尺度三维形变的方法 |
| CN114322841B (zh) * | 2021-12-06 | 2022-11-18 | 重庆邮电大学 | 一种投影光栅相移生成动态三维测量方法及*** |
| CN115307746B (zh) * | 2022-08-10 | 2023-08-18 | 苏州科技大学 | 一种单色光的干涉测量相位计算方法 |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH081371B2 (ja) * | 1985-11-14 | 1996-01-10 | 大阪大学長 | 実時間変位分布測定方法および装置 |
| JP4831703B2 (ja) * | 2008-04-23 | 2011-12-07 | 国立大学法人 和歌山大学 | 物体の変位測定方法 |
| US9389068B2 (en) | 2012-03-14 | 2016-07-12 | National Institute Of Advanced Industrial Science And Technology | Method and device for analysing phase distribution of fringe image using high-dimensional intensity information, and program for the same |
| JP6120459B2 (ja) * | 2013-07-18 | 2017-04-26 | 国立研究開発法人産業技術総合研究所 | 規則性模様による変位分布のための測定方法、装置およびそのプログラム |
| TWI471522B (zh) | 2013-07-25 | 2015-02-01 | Nat Univ Tsing Hua | 使用陰影雲紋法線上即時檢測表面形貌與平面外變形之方法與系統 |
| WO2015102546A1 (en) * | 2013-12-31 | 2015-07-09 | Aselsan Elektronik Sanayi Ve Ticaret Anonim Sirketi | A liquid cooled device enclosure |
| JP6323929B2 (ja) * | 2014-01-30 | 2018-05-16 | 国立大学法人 和歌山大学 | 計測装置及び橋梁検査方法 |
| JP6472675B2 (ja) * | 2015-02-05 | 2019-02-20 | 国立研究開発法人産業技術総合研究所 | 二次モアレ縞による顕微鏡走査ゆがみの影響を受けない変形測定方法 |
| JP6813162B2 (ja) | 2016-02-17 | 2021-01-13 | 国立研究開発法人産業技術総合研究所 | モアレ法による高速変位・ひずみ分布測定方法及び測定装置 |
-
2017
- 2017-06-01 JP JP2018541896A patent/JP6583761B2/ja active Active
- 2017-06-01 EP EP17855280.8A patent/EP3489621A4/en not_active Withdrawn
- 2017-06-01 US US16/328,309 patent/US10655954B2/en not_active Expired - Fee Related
- 2017-06-01 WO PCT/JP2017/020435 patent/WO2018061321A1/ja unknown
Also Published As
| Publication number | Publication date |
|---|---|
| JPWO2018061321A1 (ja) | 2019-04-18 |
| US10655954B2 (en) | 2020-05-19 |
| EP3489621A1 (en) | 2019-05-29 |
| EP3489621A4 (en) | 2020-04-08 |
| WO2018061321A1 (ja) | 2018-04-05 |
| US20190212130A1 (en) | 2019-07-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6583761B2 (ja) | 周期模様を利用した三次元形状・変位・ひずみ測定装置、方法およびそのプログラム | |
| JP5818218B2 (ja) | 高次元輝度情報を用いた縞画像の位相分布解析方法、装置およびそのプログラム | |
| JP4831703B2 (ja) | 物体の変位測定方法 | |
| JP6813162B2 (ja) | モアレ法による高速変位・ひずみ分布測定方法及び測定装置 | |
| JP6534147B2 (ja) | 単一カメラによる物体の変位と振動の測定方法、装置およびそのプログラム | |
| Sutton et al. | Metrology in a scanning electron microscope: theoretical developments and experimental validation | |
| JP6753592B2 (ja) | 変形測定方法、変形測定装置、及びそのプログラム | |
| JP6983402B2 (ja) | 汚れの影響を受けにくいひずみ分布測定方法とそのプログラム | |
| KR102521925B1 (ko) | 잔류 열변형 측정 방법, 잔류 열변형 측정 장치, 및 그 프로그램 | |
| Zhu et al. | Accurate 3D measurement system and calibration for speckle projection method | |
| WO2015008404A1 (ja) | 規則性模様による変位分布のための測定方法、装置およびそのプログラム | |
| JP5130513B2 (ja) | ３次元変位ひずみ計測方法及び装置 | |
| JP5466325B1 (ja) | 物体に取り付けた格子の画像から物体の物理量を測定する方法 | |
| Ri | Accurate and fast out-of-plane displacement measurement of flat objects using single-camera based on the sampling moiré method | |
| CN110631487B (zh) | 一种利用激光散斑自相关技术测量横向微位移的方法 | |
| JP6533914B2 (ja) | 計測方法、計測装置、計測プログラム及び計測プログラムを記録した、コンピュータ読み取り可能な記録媒体 | |
| JP6472675B2 (ja) | 二次モアレ縞による顕微鏡走査ゆがみの影響を受けない変形測定方法 | |
| Morimoto et al. | Theory and application of sampling moiré method | |
| JP2013185942A (ja) | 計測方法 | |
| Dowling | Three Dimensional Digital Image Correlation for Dynamic Measurements | |
| Fuqiang | Development of Digital Image Correlation Technique for Engineering Measurement | |
| Arikawa et al. | Pattern position correction for measuring large deformation in speckle interferometry | |
| Levin et al. | Measurement of nanomovements of an object from the optical phase image | |
| Zhu et al. | Measurement of three-dimensional deformation using speckle projection method | |
| Ri et al. | Hybrid Moiré Method for Three-dimensional Stress Analysis |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20181207 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20190219 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20190228 |
|
| 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: 20190806 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20190822 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 6583761 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |