CN101769722B - 一种外差式时间序列散斑干涉测量物体变形的方法 - Google Patents
一种外差式时间序列散斑干涉测量物体变形的方法 Download PDFInfo
- Publication number
- CN101769722B CN101769722B CN2010101012525A CN201010101252A CN101769722B CN 101769722 B CN101769722 B CN 101769722B CN 2010101012525 A CN2010101012525 A CN 2010101012525A CN 201010101252 A CN201010101252 A CN 201010101252A CN 101769722 B CN101769722 B CN 101769722B
- Authority
- CN
- China
- Prior art keywords
- light
- deformation
- speckle
- heterodyne
- place
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000005305 interferometry Methods 0.000 title claims abstract description 17
- 230000002123 temporal effect Effects 0.000 title claims abstract description 9
- 230000010287 polarization Effects 0.000 claims abstract description 27
- GQYHUHYESMUTHG-UHFFFAOYSA-N lithium niobate Chemical compound [Li+].[O-][Nb](=O)=O GQYHUHYESMUTHG-UHFFFAOYSA-N 0.000 claims abstract description 5
- 238000001914 filtration Methods 0.000 claims description 5
- 238000012360 testing method Methods 0.000 claims description 5
- 238000002834 transmittance Methods 0.000 claims description 5
- 230000002452 interceptive effect Effects 0.000 claims description 2
- 238000001514 detection method Methods 0.000 abstract description 5
- XEEYBQQBJWHFJM-UHFFFAOYSA-N iron Substances [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 23
- 229910052742 iron Inorganic materials 0.000 description 12
- 238000005516 engineering process Methods 0.000 description 10
- 230000008569 process Effects 0.000 description 6
- 238000006073 displacement reaction Methods 0.000 description 5
- 238000005259 measurement Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 101100117236 Drosophila melanogaster speck gene Proteins 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000005210 holographic interferometry Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 230000003746 surface roughness Effects 0.000 description 1
- 238000011426 transformation method Methods 0.000 description 1
Images
Landscapes
- Length Measuring Devices By Optical Means (AREA)
- Instruments For Measurement Of Length By Optical Means (AREA)
Abstract
本发明公开了属于激光散斑干涉测量领域的一种外差式时间序列散斑干涉测量物体变形的方法。该方法结合了时间序列散斑干涉测量和外差式干涉测量方法,直接给出时变场的相位函数,并引入频率差,确保形变测量精度。本发明简化了物体变形检测系统,并通过使用铌酸锂偏振控制器更好地适应了高速CCD的存储速度。
Description
技术领域
本发明涉及激光散斑干涉测量领域,尤其涉及一种外差式时间序列散斑干涉测量物体变形的方法。
背景技术
当激光照射在具有漫反射性质的物体表面时,根据惠更斯原理,物体表面每一点都可以看成一个点光源,从物体表面反射的光在空间相干叠加,就会在整个空间发生干涉,形成随机分布的亮斑和暗斑,称为激光散斑。散斑现象普遍存在于光学成像过程中,随着物体的变形或运动而变化。最初人们主要研究的是如何减少散斑带来的影响,但在研究过程中发现,物体变形前后,散斑有一定的规律,且携带了物体表面变形或位移的信息,逐渐产生了许多应用。20世纪70年代初,激光散斑干涉测量方法得到了发展,它除了具有全息干涉测量方法的非接触式、可以遥感、直观、能给出全场情况等一系列优点外,还具有光路简单,对试件表面要求不高,对实验条件要求较低,计算方便等特点。激光散斑干涉测量技术的用途广泛,除了测量物体的位移、应变外,还可以用于无损探伤、物体表面粗糙度测量、振动测量等方面。
电子散斑干涉(ESPI)是在本世纪初就已经被广泛应用于漫射体表面位移或变形监测的现代光学计量方法。在以往的计量中,它们主要是以强度相关条纹的形式来表征所需要测量的位移或者变形量的,但是由于复杂的环境影响和条纹本身的复杂性与随机性,很难形成统一的条纹处理模式与定量解释方式。为了解决散斑计量中的定量检测问题,在ESPI技术中引入了相移技术(PSESPI)。该方法可以方便地对变形场进行定量检测,但需要同一变形状态下的多幅干涉图进行相位求解。同时相移装置的引入也增加了实验系统的复杂性和对环境噪声的敏感性,因此PSESPI不适合时变环境下的变形检测。载波电子散斑干涉技术(CESPI)是在ESPI装置中引入载波条纹,利用傅里叶变换技术,就可以解调出全场相位,因此它可用于无法使用PSESPI的场合,如连续变形或运动场等。
从传统的散斑法到ESPI、PSESPI和CESPI,位移场的测量都是基于两个状态,即变形前后的状态,都没有涉及到时间参数。随着大容量计算机、高速CCD等设备的出现,使得在一定时间间隔内记录数字序列图像变成了现实。90年代末,德国的C.Joenathan等人提出了一种新的测量物体变形的技术,称为时间序列散斑干涉技术(TSPI)。TSPI是将整个物体变形的全过程记录下来,然后用傅里叶变换技术提取出物体变形的信息,观察物体上每一个点随时间的变化。但是存在的一个问题是,测量精度不高;一旦待测物体上个别点的变形小于5λ,TSPI难以精确测量出这些点的变形。为了解决这个问题,本世纪初C.Joenathan等人又提出了引入恒定的频率差的方法,即外差干涉测量。C.Joenathan使用的产生频差的方法是:激光器射出的线偏振光经过空间滤波后,入射到一个1/2波片和一个1/4波片组成的频移器,其中1/2波片通过外界机械制动,以一定的角频率转动,1/4波片固定不动。当线偏振光经过旋转的1/2波片和固定的1/4波片后,就会变成含有两个频率的正交线偏振光,两束光具有一定的频率差。但是此方法也存在一定的缺点:1、旋转1/2波片属于机械式旋转,稳定性不高,且机械旋转产生的振动对整个光学系统有一定的影响;2、旋转波片时会产生一定的偏摆,使得旋转中的1/2波片所在的平面无法保证在旋转中能够垂直于入射光,对测量过程中的相位会有一定的影响;3、C.Joenathan采用的方法是激光器出射的光先过空间滤波器再过频移器,由于激光过空间滤波器后光束会扩大,而1/2波片和1/4波片由于制造技术的限制,使得通光孔径较小,阻挡了扩束后的激光,使得测量面积减小。
发明内容
本发明的目的是针对背景技术中所描述的外差干涉测量技术存在的问题,提出了一种外差式时间序列散斑干涉测量物体变形的方法。
其特征在于,包括以下步骤:
步骤一:选择He-Ne激光器,输出频率是ω的线偏振光;
步骤二:线偏振光通过偏振控制器后变成含有两个频率的正交线偏振光;
步骤三:上述含有两个频率的正交线偏振光经过空间滤波器进行空间滤波和扩束;
步骤四:扩束后的含有两个频率的正交线偏振光透射过普通分光镜,入射偏振分光棱镜并被分开成两束;
其中一束光被偏振分光棱镜反射向上通过检偏器入射平面反射镜,被反射后再次通过检偏器,在偏振分光棱镜处被反射到普通分光镜处,作为参考光;
另一束光透射过偏振分光棱镜后,入射到待测物体上,发生散射,带有物体变形信息的散射光反射回偏振分光棱镜处并透过它入射到普通分光镜处,作为测量光;
步骤五:所述测量光和参考光被普通分光镜反射,通过检偏器后干涉,形成物体变形散斑图;该散斑图成像于高速CCD上;
步骤六:对高速CCD获得的散斑图进行傅里叶变换,获得调制频率图和频谱图,其中,通过调制频率图确定物体各个部位变形的方向;通过频谱图获得物体基于时间的变形量。
所述偏振控制器为铌酸锂偏振控制器。
本发明结合了时间序列散斑干涉测量和外差式干涉测量方法的优点,简化了检测系统,很大程度地提高了测量整个物体变形的精度;通过使用铌酸锂偏振控制器更好地适应了高速CCD的存储速度。
附图说明
图1:外差式时间序列散斑干涉测量物体变形的方法的光路图。
具体实施方式
下面结合附图,对优选实施例作详细说明。应该强调的是,下述说明仅仅是示例性的,而不是为了限制本发明的范围及其应用。
如图1所示:光源是He-Ne激光器1,输出频率是ω的线偏振光,耦合进入铌酸锂偏振控制器2中,该偏振控制器通过外加电压得到含有两个频率的正交线偏振光,且通过改变电压改变外差调制频率。含有两个频率的正交线偏振光的琼斯矢量为:
其中,A是振幅。含有两个频率的正交线偏振光通过空间滤波器3进行空间滤波和扩束,然后透射过普通分光镜4,入射到偏振分光棱镜10。在偏振分光棱镜处,含有两个频率的正交线偏振光被分开,一束光被偏振分光棱镜反射向上通过检偏器11入射平面反射镜12,被反射后再次通过检偏器11,在偏振分光棱镜处被反射回到普通分光镜处,称为参考光;另一束光透射过偏振分光棱镜后,入射到待测物体—铁板13上,发生散射,带有物体变形信息的散射光反射回偏振分光棱镜处并透过它入射到普通分光镜处,称为测量光。
测量光和参考光被普通分光镜反射,通过检偏器5后干涉,形成散斑图。远心成像系统由第一凸透镜6、夹缝7和第二凸透镜8组成,作用是将散斑图成像于高速CCD 9上,高速CCD可将铁板变形的过程记录下来,散斑图强度函数可表示为:
I(x,y,t)=I0(x,y){1+Vcos[Φ0(x,y)+4ω0t]}
其中,I0(x,y)是干涉场平均强度,V是对比调制能见度,Φ0(x,y)是随机相位,4ω0是前面产生的外差调制频率,用来调制散斑强度。当待测铁板没有变形的时候,外差调制频率4ω0可在高速CCD上观察到。当铁板发生变形后,上式变为:
其中,Δz(x,y,t)是铁板变形函数,当铁板变形时向着偏振分光棱镜靠近时,该值为正;当铁板变形是远离偏振分光棱镜时,该值为负。
时间序列散斑干涉测量就是利用傅里叶变换方法将铁板变形函数提取出来,进而直观的观察铁板的变形。
在一定时间t内的铁板变形函数Δz产生的中值调制频率为:
铁板上不同的点在变形中产生的变形对应不同的中值调制频率,其值将会与外差调制频率相加或者相减,这取决于该点变形的方向。
铁板进行连续的运动后,引起散斑强度的时间调制变化,这个过程被高速CCD采集到大量的散斑图。对高速CCD获得的物体变形散斑强度图进行傅里叶变换,处理数据可获得调制频率图和频谱图。从获得的调制频率的大小可以得出待测物体各个部位变形的方向,调制频率增大则说明该点向着偏振分光镜移动;反之,调制频率减小则说明该点远离偏振分光棱镜。对于频谱图,其形式为:
F{I(x,y)}=A(f·y)+Q(f+f0·y)+Q*(f-f0·y) (2)
上式为(1)式的傅里叶变换,对其进行滤波,得到正一级或负一级分量,这里选取正一级分量Q(f+f0·y)。接着对其进行逆傅里叶变换:
F-1{Q}=q(x,y)
其中相位函数Φ(x,y)就包含在q(x,y)中,又:
Φ(x,y)=tan-1{Im[q(x,y)]/Re[q(x,y)]} (3)
其中,Re[q(x,y)]是q(x,y)的实部,Im[q(x,y)]是q(x,y)的虚部,经过(3)式即可得到相位函数:
最后对获得的相位函数Φ(x,y)进行展开,即可获得变形函数Δz(x,y,t)。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。
Claims (2)
1.一种外差式时间序列散斑干涉测量物体变形的方法,其特征在于,包括以下步骤:
步骤一:选择He-Ne激光器,输出频率是ω的线偏振光;
步骤二:线偏振光通过偏振控制器后变成含有两个频率的正交线偏振光;
步骤三:上述含有两个频率的正交线偏振光经过空间滤波器进行空间滤波和扩束;
步骤四:扩束后的含有两个频率的正交线偏振光透射过普通分光镜,入射偏振分光棱镜并被分开为两束;
其中一束光被偏振分光棱镜反射向上通过检偏器入射平面反射镜,被反射后再次通过检偏器,在偏振分光棱镜处被反射到普通分光镜处,作为参考光;
另一束光透射过偏振分光棱镜后,入射到待测物体上,发生散射,带有物体变形信息的散射光反射回偏振分光棱镜处并透过它入射到普通分光镜处,作为测量光;
步骤五:所述测量光和参考光被普通分光镜反射,通过检偏器后干涉,形成物体变形散斑图;该散斑图成像于高速CCD上;
步骤六:对高速CCD获得的散斑强度图进行傅里叶变换,获得调制频率图和频谱图,其中,通过调制频率图确定物体各个部位变形的方向;通过频谱图获得物体基于时间的变形量。
2.根据权利要求1所述的一种外差式时间序列散斑干涉测量物体变形的方法,其特征在于,所述偏振控制器为铌酸锂偏振控制器。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101012525A CN101769722B (zh) | 2010-01-26 | 2010-01-26 | 一种外差式时间序列散斑干涉测量物体变形的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101012525A CN101769722B (zh) | 2010-01-26 | 2010-01-26 | 一种外差式时间序列散斑干涉测量物体变形的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101769722A CN101769722A (zh) | 2010-07-07 |
CN101769722B true CN101769722B (zh) | 2011-08-17 |
Family
ID=42502725
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010101012525A Expired - Fee Related CN101769722B (zh) | 2010-01-26 | 2010-01-26 | 一种外差式时间序列散斑干涉测量物体变形的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101769722B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102221342B (zh) * | 2011-04-02 | 2013-04-24 | 北京交通大学 | 一种时域多波长外差散斑干涉测量物体变形的方法 |
CN102506716B (zh) * | 2011-10-24 | 2013-09-11 | 河南科技大学 | 面内位移和离面位移同时测量的激光散斑测量装置和方法 |
CN103144036B (zh) * | 2012-12-31 | 2015-07-22 | 南京理工大学 | 环抛阶段磨削量在线监测装置及磨削量在线监测方法 |
CN103148787B (zh) * | 2013-03-05 | 2016-12-07 | 北京航空航天大学 | 一种古斯汉欣位移传感测量方法及系统 |
FR3022346B1 (fr) * | 2014-06-16 | 2022-10-07 | Commissariat Energie Atomique | Dispositif et procede de caracterisation d'un faisceau de lumiere |
CN104949940B (zh) * | 2015-06-19 | 2018-03-06 | 苏州大学 | 一种测量散射物体散射函数实部和虚部的装置和方法 |
WO2018072143A1 (zh) * | 2016-10-19 | 2018-04-26 | 北京交通大学 | 一种动态激光多普勒成像全场测量方法 |
CN107167087A (zh) * | 2017-05-12 | 2017-09-15 | 天津大学 | 一种基于经验小波变换的物体变形测量方法 |
CN109039449A (zh) * | 2018-08-22 | 2018-12-18 | 安徽问天量子科技股份有限公司 | 一种光纤线缆覆冰的检测方法及其检测装置 |
CN110864640A (zh) | 2018-08-28 | 2020-03-06 | 合肥京东方显示技术有限公司 | 光学系统及利用感光相机测量物体应变的方法 |
CN110631487B (zh) * | 2019-11-08 | 2020-12-11 | 济南大学 | 一种利用激光散斑自相关技术测量横向微位移的方法 |
CN111721492A (zh) * | 2020-06-30 | 2020-09-29 | 武汉大学 | 电子产品抗冲击性能测试装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1617173A2 (en) * | 2004-07-13 | 2006-01-18 | Mitutoyo Corporation | System and method for improving accuracy in a speckle-based image correlation displacement sensor |
CN1932433A (zh) * | 2006-09-29 | 2007-03-21 | 山东师范大学 | 载频电子散斑位移场的分离方法 |
CN1971207A (zh) * | 2006-11-28 | 2007-05-30 | 山东师范大学 | 大剪切载频电子散斑干涉位移场分离方法 |
CN101520306A (zh) * | 2009-03-30 | 2009-09-02 | 哈尔滨工业大学 | 基于空间载波的干涉共焦测量装置与方法 |
-
2010
- 2010-01-26 CN CN2010101012525A patent/CN101769722B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1617173A2 (en) * | 2004-07-13 | 2006-01-18 | Mitutoyo Corporation | System and method for improving accuracy in a speckle-based image correlation displacement sensor |
CN1932433A (zh) * | 2006-09-29 | 2007-03-21 | 山东师范大学 | 载频电子散斑位移场的分离方法 |
CN1971207A (zh) * | 2006-11-28 | 2007-05-30 | 山东师范大学 | 大剪切载频电子散斑干涉位移场分离方法 |
CN101520306A (zh) * | 2009-03-30 | 2009-09-02 | 哈尔滨工业大学 | 基于空间载波的干涉共焦测量装置与方法 |
Non-Patent Citations (1)
Title |
---|
C.Joenathan etc..Speckle interferometry with temporal phase evaluation for measuring large-object deformation.《APPLIED OPTICS》.1998,第37卷(第13期),2608-2614. * |
Also Published As
Publication number | Publication date |
---|---|
CN101769722A (zh) | 2010-07-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101769722B (zh) | 一种外差式时间序列散斑干涉测量物体变形的方法 | |
CN102221342B (zh) | 一种时域多波长外差散斑干涉测量物体变形的方法 | |
CN102788558B (zh) | 散斑相关和散斑干涉相结合的三维变形测量系统及方法 | |
CN102322807B (zh) | 一种物体动态三维变形的实时测量方法 | |
CN104296698A (zh) | 一种超高精度的光学表面平整度测量方法 | |
CN102865811B (zh) | 基于正交双光栅的同步移相共光路干涉显微检测装置及检测方法 | |
US20180128591A1 (en) | Noise reduction techniques, fractional bi-spectrum and fractional cross-correlation, and applications | |
CN103245285A (zh) | 一种反射式点衍射载波同步移相干涉检测装置及检测方法 | |
Tiziani et al. | From speckle pattern photography to digital holographic interferometry | |
CN102230827B (zh) | 马赫-曾德点衍射干涉仪及激光复振幅重建方法 | |
EP3845857A1 (en) | Surface shape measurement device and surface shape measurement method | |
CN103615993B (zh) | 基于离轴显微干涉术的微结构测试系统及方法 | |
CN105571517A (zh) | 一种适用于光纤端面检测的改进型相干峰解调方法 | |
CN106247950A (zh) | 基于广义相移数字全息的微小位移测量方法 | |
Wang et al. | All-fiber differential interferometer for nanometric displacement measurement | |
Ruiz et al. | Depth-resolved imaging and displacement measurement techniques viewed as linear filtering operations | |
CN105674875B (zh) | 一种全视场低频外差点衍射干涉仪 | |
Morris et al. | Noise reduction in dynamic interferometry measurements | |
CN105571516A (zh) | 一种全视场低频外差干涉仪 | |
CN100529656C (zh) | 利用电子散斑干涉载频调制技术测量物体形貌的方法 | |
CN103322912B (zh) | 一种反射式点衍射离轴同步移相干涉检测装置与检测方法 | |
CN104280139A (zh) | 一种动态相位测量装置及方法 | |
CN106091974A (zh) | 一种物体形变测量仪器、方法和设备 | |
El-Morsy | A novel algorithm based on sub-fringe integration method for direct two-dimensional unwrapping phase reconstruction from the intensity of one-shot two-beam interference fringes | |
Singh et al. | Nanoscale measurement of in-plane and out-of-plane displacements of microscopic object by sensor fusion |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20110817 |
|
CF01 | Termination of patent right due to non-payment of annual fee |