CN108007375B - 一种基于合成波长双光源剪切散斑干涉的三维变形测量方法 - Google Patents

一种基于合成波长双光源剪切散斑干涉的三维变形测量方法 Download PDF

Info

Publication number
CN108007375B
CN108007375B CN201711367603.5A CN201711367603A CN108007375B CN 108007375 B CN108007375 B CN 108007375B CN 201711367603 A CN201711367603 A CN 201711367603A CN 108007375 B CN108007375 B CN 108007375B
Authority
CN
China
Prior art keywords
phase
wavelength
plane displacement
plane
shearing
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.)
Expired - Fee Related
Application number
CN201711367603.5A
Other languages
English (en)
Other versions
CN108007375A (zh
Inventor
魏连锁
郭媛
胡现成
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.)
Qiqihar University
Original Assignee
Qiqihar University
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 Qiqihar University filed Critical Qiqihar University
Priority to CN201711367603.5A priority Critical patent/CN108007375B/zh
Publication of CN108007375A publication Critical patent/CN108007375A/zh
Application granted granted Critical
Publication of CN108007375B publication Critical patent/CN108007375B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/16Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge
    • G01B11/161Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge by interferometric means
    • G01B11/162Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge by interferometric means by speckle- or shearing interferometry

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

一种基于合成波长双光源剪切散斑干涉的三维变形测量方法,属于光学成像技术领域。所述方法如下:红色和绿色激光分别通过光分束器,经反光镜和光束扩展器照射在被测物体上,在其表面发生反射,反射光经凸透镜通过剪切镜进行光的干涉,在CCD相机中采集干涉图像;通过压电陶瓷控制器控制PZT相移量,红色激光经彩色相机R通道,绿色激光经彩色相机G通道,分别被CCD相机同时记录;通过压力阀调节真空箱中的压强,压强值通过压力表显示;对干涉图像进行傅里叶变换,推算合成波长相位条纹数与单波长条纹数关系;计算面内位移导数阵和离面位移导数阵,根据该测量系统的光学几何特性,进而得到物体三维变形阵,实现物体三维变形测量。

Description

一种基于合成波长双光源剪切散斑干涉的三维变形测量方法
技术领域
本发明属于光学成像技术领域,具体涉及一种基于合成波长双光源剪切散斑干涉的三维变形测量方法。
背景技术
作为主要测量方法的数字散斑干涉法和混合光测方法在三维变形测量中已经取得了一定的研究成果,但仍然存在以下几个方面的问题:
(1)测量光路设置复杂操作不便、测量范围小、对环境和硬件设备要求高,不能满足工程环境测量。
(2)当测量由缺陷引起的深度变形时,导致CCD上采集到的散斑干涉条纹图过于密集,出现欠采样问题,常用的数字图像处理方法失效,无法得到真实相位。
(3)在散斑干涉条纹图像处理上,相位去噪滤波、提取和解包裹算法还存在去噪效果不明显、丢失或破坏有用信息、相位提取算法复杂度高、大计算量解包运行时间长等。
发明内容
本发明的目的是为了解决三维变形测量中存在的精度低、抗干扰性差的问题,提供一种基于合成波长双光源剪切散斑干涉的三维变形测量方法。
为实现上述目的,本发明采取的技术方案如下:
一种基于合成波长双光源剪切散斑干涉的三维变形测量方法,所述方法步骤如下:
步骤一:两束不同波长的红色激光和绿色激光通过光分束器,并经反光镜和光束扩展器照射放置在真空箱内的被测物体上,在被测物体表面发生反射,反射光经凸透镜L1并在剪切干涉装置中通过剪切镜进行光的干涉,在CCD相机中采集干涉图像;
步骤二:通过PZTcontroller控制PZT相移量,红色激光经彩色相机R通道,绿色激光经彩色相机G通道,分别被CCD相机同时记录,R通道和G通道位于剪切干涉装置内部;
步骤三:通过压力阀调节真空箱中的压强,压强值通过压力表显示;由于步骤一采集的干涉图像中包含着红色光谱信息和绿色光谱信息,在工作波长λr的情况下对被测物体表面进行移相干涉测量,相应得到N个釆样点,进而得到N个相位数据φ1r,φ2r,…,φNr;然后,更换工作波长为λg,同样能够得到N个相位数φ1g,φ2g,…,φNg
步骤四:对干涉图像进行傅里叶变换,分离红色光谱信息和绿色光谱信息,由相移法提取包裹相为:红色波长相位变化量减去绿色波长相位变化量即得合成波长相位变化量,相位条纹数等于推得合成波长相位条纹数与单波长条纹数关系其中,为相位改变量,Nr为红色波长相位条纹数,Ns为合成波长相位条纹数,λr为红色波长,λs为合成波长。
如果剪切量是沿x轴向,则相位变化可表示为
式中,u,v,w分别是物体沿x,y,z轴向的位移分量,α、β分别为x-z面和y-z面的光照角,λ为波长。
当激光器和CCD相机置于x-z面,且剪切沿x轴向,则合成光束1采集的相位变化为
合成光束2采集的相位变化为
由此,可计算面内位移相位变化和离面位移相位变化;
(1)面内位移相位变化计算
沿x、y轴向的面内位移相位变化为
由上两式,可得物体面内位移相位变化的一阶导数
(2)离面位移相位变化计算
剪切沿x、y轴向的离面位移相位变化分别为:
由上两式,可得离面位移相位变化的一阶导数
当光照角α很小且接近于0时,则离面位移导数为
步骤五:面内位移导数阵Gin和离面位移导数阵Gout可分别表述为
式中,u、v、w分别是x、y、z轴向的位移分量,则空间三维变形阵G可写为
根据该测量系统的光学几何特性有:当合成光束1所采集的相位减去合成光束2所采集的相位,即为面内位移相位,通过面内位移相位变化可得物体面内位移的一阶导数;当合成光束1照射时采集的相位加上合成光束2照射时所采集的相位,即为离面位移相位,通过离面位移相位变化,能够得到物体离面位移的一阶导数,从而得到物体的三维变形阵,因此能够实现物体三维变形的测量。
本发明相对于现有技术的有益效果是:
(1)采用两个不同波长的红光合成光束1和绿光合成光束2照射被测物体,设计光路简单、非接触、全场测量、精度高、速度快的三维变形测量系统,解决三维变形测量中存在的精度低、抗干扰性差等问题,从而满足工程环境三维变形的测量要求,有望为物体三维变形测量提供一种新的光学研究方法。
(2)采用上下双光源照射,同步得到物体面内位移和离面位移,解决目前物体三维变形量无法同步获得的问题。合成等效波长光源的使用,使光束中心光斑能量增强,波长更长,测量范围更大,有助于降低频谱混叠程度,大幅度降低相位条纹密度,解决深变形导致干涉条纹过密,易使后续相位处理失败,无法获得三维变形的情况。
附图说明
图1为合成等效波长的上下双光源剪切散斑干涉测量原理图。
具体实施方式
下面结合附图和实施例对本发明的技术方案作进一步的说明,但并不局限于此,凡是对本发明技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围,均应涵盖在本发明的保护范围中。
具体实施方式一:如图1所示,本实施方式记载的是一种基于合成波长双光源剪切散斑干涉的三维变形测量方法,所述方法步骤如下:
步骤一:两束不同波长的红色激光和绿色激光通过光分束器,并经反光镜和光束扩展器照射放置在真空箱内的被测物体上,在被测物体表面发生反射,反射光经凸透镜L1并在剪切干涉装置中通过剪切镜进行光的干涉,在CCD相机中采集干涉图像;
步骤二:通过PZTcontroller(压电陶瓷控制器)控制PZT相移量,红色激光经彩色相机R通道,绿色激光经彩色相机G通道,分别被CCD相机同时记录,R通道和G通道位于剪切干涉装置内部;
步骤三:通过压力阀调节真空箱中的压强,压强值通过压力表显示;由于步骤一采集的干涉图像中包含着红色光谱信息和绿色光谱信息,在工作红光波长λr的情况下对被测物体表面进行移相干涉测量,相应得到N个釆样点,进而得到N个相位数据φ1r,φ2r,…,φNr;然后,更换绿波工作波长为λg,同样能够得到N个相位数φ1g,φ2g,…,φNg
步骤四:对干涉图像进行傅里叶变换,分离红色光谱信息和绿色光谱信息,由相移法提取包裹相为:红色波长相位变化量减去绿色波长相位变化量即得合成波长相位变化量,相位条纹数等于推得合成波长相位条纹数与单波长条纹数关系来有效降低干涉条纹的密集程度,避免由干涉条纹过于密集而导致欠采样的问题;
如果剪切量是沿x轴向,则相位变化可表示为
式中,u,v,w分别是物体沿x,y,z轴向的位移分量,α、β分别为x-z面和y-z面的光照角,λ为波长。
当激光器和CCD相机置于x-z面,且剪切沿x轴向,则合成光束1采集的相位变化为
合成光束2采集的相位变化为
由此,计算面内位移相位变化和离面位移相位变化;
(1)面内位移相位变化计算
沿x、y轴向的面内位移相位变化为
由上两式,得物体面内位移相位变化的一阶导数
(2)离面位移相位变化计算
剪切沿x、y轴向的离面位移相位变化分别为:
由上两式,得离面位移相位变化的一阶导数
当光照角α很小且接近于0时,则离面位移导数为
步骤五:面内位移导数阵和离面位移导数阵可分别表述为
式中,u、v、w分别是x、y、z轴向的位移分量,则空间三维变形阵可写为
根据该测量系统的光学几何特性有:当合成光束1所采集的相位减去合成光束2所采集的相位,即为面内位移相位,通过面内位移相位变化可得物体面内位移的一阶导数;当合成光束1照射时采集的相位加上合成光束2照射时所采集的相位,即为离面位移相位,通过离面位移相位变化,得到物体离面位移的一阶导数,从而能够得到物体的三维变形阵,因此能够实现物体三维变形的测量。
图像采集卡采集数幅干涉图样,经PC滤波和相位解包裹得面内位移和离面位移导数,由此实现了物体三维变形的测量。
在分析现有散斑干涉光路结构及成像原理的基础上,建立新的剪切散斑干涉法测量三维变形光路测量系统,为进一步研究满足工程环境三维变形的测量系统奠定基础。
物体三维变形可分为两部分,即面内变形(即面内位移)和离面变形(即离面位移)。双束合成等效波长的上下双光源剪切散斑干涉测量系统基本原理如图所示,双束合成等效波长激光器和CCD相机置于xoz面,剪切沿y轴向。
当合成等效波长激光束照在漫射体表面,携带被测物体表面信息的反射光经凸透镜L1,在剪切干涉装置中进行光的干涉,进而在CCD上相干成像。P[P(x,y,z)]是物体变形前测试面上的任意一点,物体变形后则为点P′[P(x+u,y+v,z+w)],由于变形前后照射在点P上的光程不等于照射在点P′上的光程,光程的改变引起了相位的变化,通过不同面上相位差的测量可得物体的三维变形。物体变形前,通过PZT产生n步相移,在不同相移下分别记录同一状态的n幅干涉图,经空域相移法可得物体变形前的相位,同理,可得物体变形后的相位,则物体变形后的相位减去物体变形前的相位即为物体变形引起的相位变化量。
与剪切干涉装置中心线所成角度α的激光束为合成光束1,与剪切干涉装置中心线所成角度-α的激光束为合成光束2(其中-π/2<α<π/2)。
实施例1
假设选定该测量系统中,二极管泵浦激光器发出的是绿光,波长λg=532nm,氦氖激光器发出的是红色光,波长为λr=632.8nm,由于激光有很好的相干性,两束激光相干后,光束中心光斑能量增强,可得高功率相干光束,使得波长更长、测量范围更大,有助于降低频谱混叠程度,其合成等效波长为λs
经计算可得合成波长,即λs=3.3398μm。

Claims (1)

1.一种基于合成波长双光源剪切散斑干涉的三维变形测量方法,其特征在于:所述方法步骤如下:
步骤一:两束不同波长的红色激光和绿色激光通过光分束器,并经反光镜和光束扩展器照射放置在真空箱内的被测物体上,在被测物体表面发生反射,反射光经凸透镜L1并在剪切干涉装置中通过剪切镜进行光的干涉,在CCD相机中采集干涉图像;
步骤二:通过PZT controller控制PZT相移量,红色激光经彩色相机R通道,绿色激光经彩色相机G通道,分别被CCD相机同时记录,R通道和G通道位于剪切干涉装置内部;
步骤三:通过压力阀调节真空箱中的压强,压强值通过压力表显示;由于步骤一采集的干涉图像中包含着红色光谱信息和绿色光谱信息,在工作红光波长λr的情况下对被测物体表面进行移相干涉测量,相应得到N个釆样点,进而得到N个相位数据φ1r,φ2r,…,φNr;然后,更换工作绿光波长为λg,同样能够得到N个相位数φ1g,φ2g,…,φNg
步骤四:对干涉图像进行傅里叶变换,分离红色光谱信息和绿色光谱信息,由相移法提取包裹相为:红色波长相位变化量减去绿色波长相位变化量即得合成波长相位变化量,相位条纹数等于推得合成波长相位条纹数与单波长条纹数关系Nr为红色波长相位条纹数,Ns为合成波长相位条纹数,λr为红色波长,λs为合成波长;
如果剪切量是沿x轴向,则相位变化可表示为
式中,u,v,w分别是物体沿x,y,z轴向的位移分量,α、β分别为x-z面和y-z面的光照角,λ为波长;
当激光器和CCD相机置于x-z面,且剪切沿x轴向,则合成光束1采集的相位变化为
合成光束2采集的相位变化为
由此,能够计算面内位移相位变化和离面位移相位变化;
(1)面内位移相位变化计算
沿x、y轴向的面内位移相位变化为
由上两式,得物体面内位移相位变化的一阶导数
(2)离面位移相位变化计算
剪切沿x、y轴向的离面位移相位变化分别为:
由上两式,得离面位移相位变化的一阶导数
当光照角α很小且接近于0时,则离面位移导数为
步骤五:面内位移导数阵和离面位移导数阵可分别表述为
式中,u、v、w分别是x、y、z轴向的位移分量,则空间三维变形阵可写为
根据该测量系统的光学几何特性有:当合成光束1所采集的相位减去合成光束2所采集的相位,即为面内位移相位,通过面内位移相位变化可得物体面内位移的一阶导数;当合成光束1照射时采集的相位加上合成光束2照射时所采集的相位,即为离面位移相位,通过离面位移相位变化,得物体离面位移的一阶导数,从而得物体的三维变形阵,因此可实现物体三维变形的三维测量。
CN201711367603.5A 2017-12-18 2017-12-18 一种基于合成波长双光源剪切散斑干涉的三维变形测量方法 Expired - Fee Related CN108007375B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711367603.5A CN108007375B (zh) 2017-12-18 2017-12-18 一种基于合成波长双光源剪切散斑干涉的三维变形测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711367603.5A CN108007375B (zh) 2017-12-18 2017-12-18 一种基于合成波长双光源剪切散斑干涉的三维变形测量方法

Publications (2)

Publication Number Publication Date
CN108007375A CN108007375A (zh) 2018-05-08
CN108007375B true CN108007375B (zh) 2019-09-24

Family

ID=62059642

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711367603.5A Expired - Fee Related CN108007375B (zh) 2017-12-18 2017-12-18 一种基于合成波长双光源剪切散斑干涉的三维变形测量方法

Country Status (1)

Country Link
CN (1) CN108007375B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109539981A (zh) * 2018-10-30 2019-03-29 成都飞机工业(集团)有限责任公司 一种高温高频振动耦合环境下试件相关图像获取方法
WO2020132960A1 (zh) * 2018-12-26 2020-07-02 合刃科技(深圳)有限公司 缺陷检测方法及缺陷检测系统
CN110068283A (zh) * 2019-04-25 2019-07-30 东华大学 一种应用于面内应变的数字散斑传感系统
FR3100881B1 (fr) 2019-09-13 2021-08-13 Ifp Energies Now Procédé de mesure de la déformation d’une conduite sous-marine par interférométrie
CN110779456A (zh) * 2019-11-08 2020-02-11 桂林电子科技大学 一种太赫兹波段超表面相移装置及其测量方法
CN111024001B (zh) * 2019-12-13 2022-01-07 北京信息科技大学 滚转角测量方法及装置
CN111412850B (zh) * 2020-03-23 2021-09-21 昆明理工大学 一种基于单摄像头的高温三维数字图像相关测量系统及方法
CN114440785B (zh) * 2022-02-08 2023-02-07 山东大学 基于光干涉原理的材料光致形变系数测量装置及方法

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0275982A1 (fr) * 1987-01-20 1988-07-27 Communaute Europeenne De L'energie Atomique (Euratom) Dispositif d'analyse interférométrique des déformations d'un corps
CN1932437A (zh) * 2006-09-29 2007-03-21 山东师范大学 利用空间载频电子散斑干涉测量物体三维变形的方法
CN1971207A (zh) * 2006-11-28 2007-05-30 山东师范大学 大剪切载频电子散斑干涉位移场分离方法
CN101349549A (zh) * 2008-09-11 2009-01-21 清华大学 一种高温云纹干涉变形测量系统
CN101871769A (zh) * 2010-06-02 2010-10-27 山东师范大学 三维变形场同时载频调制的电子散斑检测方法
WO2011144413A1 (de) * 2010-05-18 2011-11-24 Siemens Aktiengesellschaft Verfahren und system zum ermitteln von 3d-geometriedaten von objekten
CN103278272A (zh) * 2013-05-31 2013-09-04 哈尔滨工业大学 基于激光散斑技术的焊接接头应力集中测试的装置及应力集中测试方法
CN105136806A (zh) * 2015-07-24 2015-12-09 合肥工业大学 基于空间载波的双方向剪切散斑干涉系统及测量方法
CN206223105U (zh) * 2016-12-14 2017-06-06 盐城工学院 基于数字散斑干涉的变形及斜率同步测量装置

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0275982A1 (fr) * 1987-01-20 1988-07-27 Communaute Europeenne De L'energie Atomique (Euratom) Dispositif d'analyse interférométrique des déformations d'un corps
CN1932437A (zh) * 2006-09-29 2007-03-21 山东师范大学 利用空间载频电子散斑干涉测量物体三维变形的方法
CN1971207A (zh) * 2006-11-28 2007-05-30 山东师范大学 大剪切载频电子散斑干涉位移场分离方法
CN101349549A (zh) * 2008-09-11 2009-01-21 清华大学 一种高温云纹干涉变形测量系统
WO2011144413A1 (de) * 2010-05-18 2011-11-24 Siemens Aktiengesellschaft Verfahren und system zum ermitteln von 3d-geometriedaten von objekten
CN101871769A (zh) * 2010-06-02 2010-10-27 山东师范大学 三维变形场同时载频调制的电子散斑检测方法
CN103278272A (zh) * 2013-05-31 2013-09-04 哈尔滨工业大学 基于激光散斑技术的焊接接头应力集中测试的装置及应力集中测试方法
CN105136806A (zh) * 2015-07-24 2015-12-09 合肥工业大学 基于空间载波的双方向剪切散斑干涉系统及测量方法
CN206223105U (zh) * 2016-12-14 2017-06-06 盐城工学院 基于数字散斑干涉的变形及斜率同步测量装置

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
丁效红等.《一种数字剪切散斑检测离面位移的计算方法》.《激光杂志》.2014, *
张霖等.《激光电子散斑干涉测微小位移》.《物理实验》.1999, *
郭媛等.《双波长剪切散斑干涉法在复合材料缺陷检测中的应用》.《光子学报》.2015, *
郭媛等.《基于剪切散斑干涉技术的物体变形动态检测》.《应用光学》.2017, *

Also Published As

Publication number Publication date
CN108007375A (zh) 2018-05-08

Similar Documents

Publication Publication Date Title
CN108007375B (zh) 一种基于合成波长双光源剪切散斑干涉的三维变形测量方法
US9062959B2 (en) Wavelength scanning interferometer and method for aspheric surface measurement
CN105066908B (zh) 一种基于多波长和多偏振态的数字全息三维形貌检测装置
Zuo et al. Direct continuous phase demodulation in digital holography with use of the transport-of-intensity equation
JP7082137B2 (ja) スペクトル制御干渉法による曲率半径測定
GB2505106A (en) Method and arrangement for short coherence holography
CN104534979A (zh) 一种多波长相移显微成像系统及方法
EP2572157A1 (en) Method and device for high resolution full field interference microscopy
Kumar et al. Portable digital holographic camera featuring enhanced field of view and reduced exposure time
CN110989314A (zh) 一种基于菲涅尔双面镜连续太赫兹波自参考数字全息成像方法
CN105159044A (zh) 基于双波长数字全息技术的反射式显微成像装置
CN105277136B (zh) 基于双波长数字全息技术的透射式显微成像装置及其方法
CN117870534A (zh) 一种基于led照明的双波长动态全息显微成像系统与方法
Licznerski et al. Novel double path shearing interferometer in corneal topography measurements
CN113758440B (zh) 一种改变观察矢量的数字全息三维形貌测量方法及装置
Abdelsalam et al. Digital holographic shape measurement using Fizeau microscopy
CN208937181U (zh) 一种可消除oct共轭镜像的定相差双路线阵谱域oct装置
US20190195615A1 (en) Cascade fourier domain optical coherence tomography
Friedman et al. Hybrid reflective interferometric system combining wide-field and single-point phase measurements
Lyalikov High-sensitive interferometric control of the quality of diffractive elements
CN110360923B (zh) 一种被测面可旋转的相移干涉仪及测量方法
Zvyagin et al. Image reconstruction in full-field Fourier-domain optical coherence tomography
Itoh et al. Coherence-based 3-D and spectral imaging and laser-scanning microscopy
Körner et al. Multilevel optical sectioning based on digital holography with a femtosecond frequency comb laser
Tome et al. Phase-measuring interferometry: applications and techniques

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190924

Termination date: 20211218