形貌与微变形全场光学同时测量方法

蒋艳鹏, 吴思进, 杨连祥

蒋艳鹏, 吴思进, 杨连祥. 形貌与微变形全场光学同时测量方法[J]. 应用光学, 2017, 38(1): 67-71. DOI: 10.5768/JAO201738.0103002
引用本文: 蒋艳鹏, 吴思进, 杨连祥. 形貌与微变形全场光学同时测量方法[J]. 应用光学, 2017, 38(1): 67-71. DOI: 10.5768/JAO201738.0103002
Jiang Yanpeng, Wu Sijin, Yang Lianxiang. Simultaneous measurement of contour and micro-deformation using full-field optical methods[J]. Journal of Applied Optics, 2017, 38(1): 67-71. DOI: 10.5768/JAO201738.0103002
Citation: Jiang Yanpeng, Wu Sijin, Yang Lianxiang. Simultaneous measurement of contour and micro-deformation using full-field optical methods[J]. Journal of Applied Optics, 2017, 38(1): 67-71. DOI: 10.5768/JAO201738.0103002

形貌与微变形全场光学同时测量方法

基金项目: 

北京市教育委员会科技计划面上项目 KM201511232004

国家自然科学基金项目 51275054

详细信息
    作者简介:

    蒋艳鹏(1991-),女,山东烟台人,硕士研究生,主要从事光电检测方面的研究工作。E-mail:939086488@qq.com

    通讯作者:

    吴思进(1979-),男,博士,讲师,主要从事数字散斑干涉等光电检测技术研究工作,E-mail: swu@bistu.edu.cn

  • 中图分类号: O439;TN247

Simultaneous measurement of contour and micro-deformation using full-field optical methods

  • 摘要: 在实际工程应用中,对材料形貌和结构变形等参量的检测是必不可少的,而且往往需要进行多参量同时测量。针对该背景,采用数字散斑干涉与数字条纹投影相结合的测量方法, 设计了一种集成光路,通过在数字散斑干涉实验光路中引入一个投影设备,实现物体表面形貌和微变形的同时测量。所提出的方法具有全场非接触测量的优点,且测量光路简单、操作方便、效率高、可靠性强。该方法的形貌测量分辨率优于10 μm,形变测量分辨率优于30 nm。
    Abstract: In practical engineering applications, contour and micro-deformation of materials and structures are essential parameters to be tested, and furthermore, usually need to be tested simultaneously. To meet this requirement, a method that combines digital speckle pattern interferometry (DSPI) and digital fringe projection (DFP) is proposed. An integrated optical setup is built by embedding a digital projector in the DSPI setup. With the new setup, surface contours and deformations can be measured simultaneously. The proposed method has the advantages of full-field, non-contact measurement, simple optical arrangement, ease of operation, high efficiency and reliability. The resolution of contour measurement is up to 10 μm, while the deformation measurement resolution is up to 30 nm.
  • 在工程领域,单一参量的测量往往比较容易,有多种成熟的技术可以选择,而多参量的同时测量则较为困难。通常通过组合多种测量技术进行多参量的测量,但如此则无法满足多技术同时同位置测量的要求。而工程中多参量的测量非常必要,比如当测量材料和结构的力学特性时,通常需要同时测量其全场形貌和微变形的分布,如此可以更加全面地反映被测对象的特性,同时还可以提高检测效率。例如在航天领域,航天器高精度抛物面天线是航天器与地面沟通的唯一装置,主要用于传导天地之间的电信号,其型面精度和变形量直接关系到传导精度。因此,在制造过程中对形貌和变形量的全场测量非常关键。又比如用于抓取捕获对接舱的航天机械臂末端关节,其结构尺寸精度和变形量直接关系到能否成功抓取,因此需要对形貌和变形进行精确测量。在其他领域同样存在许多例子需要进行全场形貌和变形的同时测量。可见,全场形貌和变形的同时测量技术能够满足包括国家重大工程在内的众多工程的测试需求,具有重要的研究价值。光学测量方法容易实现全场测量,同时具有高精度和非接触等优点,是理想的全场形貌和变形测量方法。

    随着工业技术的飞速发展以及现代工艺生产水平的提高,物体的外形设计趋于多样化,物体表面形貌的测量备受广大研究者的关注。典型的形貌测量方法采用三坐标测量机和柔性关节臂等工具进行扫描测量,无法实现全场测量,测量精度虽高,但却存在测量速度慢、价格昂贵、对环境较敏感及易对物体造成损伤等缺点。形貌的全场测量主要采用光学测量方法,主要包括激光干涉法和数字条纹投影法(digital fringe projection, DFP)等。激光干涉法测量精度高,但主要针对光学器件,测量对象较少,且难以实现大面积物体的测量[1]。数字条纹投影法通过主动向被测物投射光栅条纹,根据所得到的图像分布并结合光学三角原理计算表面形貌,具有系统简单、成本低、容易实现等优点,广泛应用于物体形貌测量中[2-3],并在计算机视觉、生物医学、考古、逆向工程等实际应用中发挥了重要作用[4-5]

    物体表面微变形的测量方法有很多,包括莫尔条纹法、数字图像相关法和数字散斑干涉法(digital speckle pattern interferometry, DSPI)等等。莫尔条纹法是利用莫尔条纹的光学放大原理,实现表面变形的全场测量,方法简单,抗干扰能力强,但测试灵敏度低[6]。数字图像相关法通过图像相关运算跟踪物体表面的随机散斑,计算全场变形和应变,方法简单,但散斑不易制造且容易改变被测物本身的属性,测试灵敏度也偏低[7]。数字散斑干涉法基于光学干涉原理进行测量,利用变形引起的干涉相位变化进行测量,具有适用对象广,测量灵敏度高等特点,其三维变形测试分辨率达到数十纳米,是理想的微变形测试技术[8-10],在航空航天[11]、精密机械[12]以及其他精密测量领域都有重要的应用。

    本文设计一种集成光路,通过光路复用方法在一个光路中同时实现数字条纹投影和数字散斑干涉技术,利用一个相机采集物体表面图像,从而实现了物体表面形貌和变形的同时同位置测量。测试系统简单,并且可靠性较高。

    图 1是数字条纹投影系统的结构图,由计算机编码正弦条纹图案,投影仪对参考平面和物体分别投影正弦条纹图,由CCD相机采集条纹图并提取相位,参考相位与被测物体相位相减后得到表示物体高度的相位,物体相位分布与高度关系为

    图  1  数字条纹投影系统结构图
    Figure  1.  System configuration of DFP
    $$ h\left( {x, y} \right) = \frac{{L \cdot \Delta \varphi \left( {x, y} \right)}}{{2{\rm{ \mathit{ π} }} \cdot \frac{d}{P} + \Delta \varphi \left( {x, y} \right)}} $$ (1)

    式中:h(x, y)为物体的高度分布;Δφ(x, y)为物体的相位变化量;L为相机光心到参考平面的垂直距离;d为投影光源与相机光心之间的距离;P为投影条纹间距。

    数字条纹投影通过建立相位与物体高度信息之间的关系,从而求取物体表面形貌,因此相位求解是关键技术。

    求解相位的方法有很多,其中最常用的是时间相移技术,即在时间序列上采集多幅条纹图,通过在相邻2幅图像之间引入特定的相位变化来实现相移。根据相移量和算法的不同,该技术可分为三步相移、四步相移、五步相移等等,其中四步相移法是目前应用比较广泛且计算相对简便的方法。

    四步相移法通过控制投影光栅的移动,在相邻两幅图之间引入π/2(即P/4条纹移动量)的相移量,当条纹投影到参考平面时记录4幅条纹图,其强度分别为

    $$ \begin{array}{l} {I_{0i}}\left( {x, y} \right) = a\left( {x, y} \right) + b\left( {x, y} \right)\cos \left[ {{\varphi _0}\left( {x, y} \right) + \frac{{\rm{ \mathit{ π} }}}{2}\left( {i - 1} \right)} \right]\\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\left( {i = 1, 2, 3, 4} \right) \end{array} $$ (2)

    式中:a(x, y)是背景光强;b(x, y)是与对比度相关的系数;φ0(x, y)是参考相位。联立以上式子可得相位:

    $$ {\varphi _0}\left( {x, y} \right) = \arctan \frac{{{I_{04}} - {I_{02}}}}{{{I_{01}} - {I_{03}}}} $$ (3)

    同理将参考平面更换为被测物体,可得相位φ1(x, y)。

    最终得到相位变化量为

    $$ \Delta \varphi \left( {x, y} \right) = {\varphi _1}\left( {x, y} \right) - {\varphi _0}\left( {x, y} \right) $$ (4)

    由于在求解相位的过程中用到了反正切函数,得到的相位值是在[-π, π)之间的包裹相位,并非实际相位值,因此需要对其进行解包裹运算,得到连续相位才能用于物体形貌测量中。

    在数字散斑干涉测量过程中,激光器发出相干光经分光镜分成两束:一束作为物光经扩束镜扩束照亮物体表面后产生漫反射;另一束作为参考光与反射的物光在相机图像传感器表面发生干涉,由相机记录散斑干涉图。

    相机采集到的物体变形前的散斑图强度为

    $$ I = {I_0}\left( {I + \gamma \cos \theta } \right) $$ (5)

    式中:θ=ϕ0-ϕr为物光和参考光之间的相位差;${I_0} = a_o^2 + a_r^2 $为背景强度;$\gamma = \frac{{2{a_o}{a_r}}}{{a_o^2 + a_r^2}} $为对比度。

    物体变形后的光强为

    $$ I' = {I_0}\left[ {1 + \gamma \cos \left( {\theta + \mathit{\Delta }} \right)} \right] $$ (6)

    式中,Δ=(ϕo-ϕr)-(ϕ′o-ϕ′r)=ϕo-ϕ′o是物体变形后引起的相位变化。

    相位差与物体变形之间的关系为

    $$ \begin{array}{l} \mathit{\Delta } = \frac{{2{\rm{ \mathit{ π} }}}}{\mathit{\lambda }}d \cdot s = \frac{{2{\rm{ \mathit{ π} }}}}{\mathit{\lambda }}\left( {Au + B\upsilon + Cw} \right) = \\ \;\;\;\;\;\frac{{2{\rm{ \mathit{ π} }}}}{\mathit{\lambda }}\left[ {\left( {\sin \alpha } \right)u + \left( {1 + \cos \alpha } \right)w} \right] \end{array} $$ (7)

    式中:s=(A, B, C)T为灵敏度矢量;而d=(u, v, w)为物体的变形矢量;u, v为面内变形量;w为面外变形量;λ为激光波长;α为入射角。当α非常小时,可得:$ \mathit{\Delta } = \frac{{4{\rm{ \mathit{ π} }}}}{\mathit{\lambda }}w$。由此可以得到物体面外变形量为

    $$ w = \frac{\lambda }{{4{\rm{ \mathit{ π} }}}}\mathit{\Delta } $$ (8)

    根据(8)式可知,要想求得物体变形量,关键是要对相位进行求解。数字散斑干涉法采用与数字条纹投影类似的方法进行相位求解,不同的是,数字散斑干涉通过压电陶瓷驱动平面镜移动来改变参考光的光程,从而引入特定的相位变化产生四步相移;而数字条纹投影则是通过光栅条纹移动,即每次移动光栅周期的1/4来产生相移。

    数字散斑干涉测量光路分两种:面外变形测量和面内变形测量。本文讨论面外变形测量,其光路原理图如图 2所示,当激光光源、相机和被测物中心组成的平面在xoz平面时,其灵敏度矢量为(sinα, 0, 1+cosα)。

    图  2  数字散斑干涉面外变形测量原理图
    Figure  2.  Schematic diagram of out-of-plane deformation measurement using DSPI

    实验光路如图 3所示。由计算机生成4幅条纹图,由投影仪投影到物体表面,然后由相机记录这些条纹图;同时激光器发出的光经分光镜分成2束:物光和参考光,2束光发生干涉产生散斑干涉图,由相机采集,即利用一个相机分别采集四步相移条纹图和散斑干涉图,然后由计算机进行处理。实验用到的投影仪型号是EPSON EB-C301MN,分辨率为1 280 pixel×800 pixel,采用CCD相机,其分辨率为1 040 pixel×1 360 pixel,激光器中心波长为532 nm。被测物是由铜片制成的曲率半径为200 mm的曲面。

    图  3  实验测量光路
    Figure  3.  Experimental setup

    首先由计算机生成4幅正弦光栅条纹图。图 4(a)是相机采集到的投影在被测物体上的条纹图。图 4 (b)是通过计算得到的物体相位包裹图,再应用相位解包裹技术得到被测物体的连续相位图,如图 4(c)所示。将被测物连续相位减去参考平面连续相位即可得到对应于物体表面高度信息的相位差分布图,如图 5(a)所示。最后通过标定得到如图 5(b)所示的物体形貌图。

    图  4  被测物投影及相位图
    Figure  4.  Sinusoidal grating projection and phase distribution
    图  5  相位差分布及三维形貌图
    Figure  5.  Phase difference distribution and three-dimensional contour

    本实验中条纹间距P=32像素,L=100 mm, d=900 mm, 通过公式(1)可知,每移动一级条纹,高度变化145 μm,而四步相移法的测相分辨率至少可以达到π/10,因此该系统形貌测量分辨率可达7 μm,优于10 μm的测量分辨率。

    通过搭建实验光路,对被测物变形进行测量,实验过程中通过旋转螺旋测微丝杆可以使被测物产生一个微小变形,此时相机分别采集物体变形前后的4幅散斑图,通过四步相移技术解相位得到物体的面外变形相位图如图 6所示。相位图解包裹后得到连续相位,最后得到二维面外变形图和三维面外变形图,如图 7(a)(b)所示。从图 7中可以看出变形量分布在0.11 μm~4.44 μm之间。

    图  6  面外变形相位图
    Figure  6.  Phase map of out-of-plane deformation
    图  7  物体变形图
    Figure  7.  Out-of-plane deformation

    数字散斑干涉技术能够实现高精度的变形测量,根据公式(8)可知,其相位图中每一级条纹对应半个波长的变形量,即相位变化2π,变形量为半个波长。考虑到相移技术至少能够实现π/10的测相分辨率,因此传统上该技术能够实现λ/20的变形测量分辨率。本实验系统测量变形的光路与传统数字散斑干涉技术一致,其中激光波长为532 nm,其测量分辨率达到27 nm,优于30 nm的离面变形测量分辨率。

    本文将数字散斑干涉与数字条纹投影2种技术集成到一个光路中, 从而实现了物体形貌与微变形的同时测量,提高了在实际测量中的工作效率。在本文研究的基础上,未来可基于所测得的表面形貌数据,对微变形数据进行修正,减小数字散斑干涉技术测量非平面物体的测量误差。

  • 图  1   数字条纹投影系统结构图

    Figure  1.   System configuration of DFP

    图  2   数字散斑干涉面外变形测量原理图

    Figure  2.   Schematic diagram of out-of-plane deformation measurement using DSPI

    图  3   实验测量光路

    Figure  3.   Experimental setup

    图  4   被测物投影及相位图

    Figure  4.   Sinusoidal grating projection and phase distribution

    图  5   相位差分布及三维形貌图

    Figure  5.   Phase difference distribution and three-dimensional contour

    图  6   面外变形相位图

    Figure  6.   Phase map of out-of-plane deformation

    图  7   物体变形图

    Figure  7.   Out-of-plane deformation

  • [1] 师途, 杨甬英, 张磊, 等.非球面光学元件的面形检测技术[J].中国光学, 2014, 7(1):26-46. http://d.old.wanfangdata.com.cn/Periodical/zggxyyygxwz201401002

    Shi Tu, Yang Yongying, Zhang Lei, et al. Surface testing methods of aspheric optical elements[J]. Chinese Optics, 2014, 7(1):26-46. http://d.old.wanfangdata.com.cn/Periodical/zggxyyygxwz201401002

    [2] 尚忠义, 李伟仙, 董明利, 等.基于四步相移光栅投影的三维形貌测量系统[J].应用光学, 2015, 36(4):584-589. http://d.old.wanfangdata.com.cn/Periodical/yygx201504015

    Shang Zhongyi, Li Weixian, Dong Mingli, et al. 3D shape measurement system on fringe projection in 4-step phase shifting[J]. Journal of Applied Optics, 2015, 36(4):584-589. http://d.old.wanfangdata.com.cn/Periodical/yygx201504015

    [3]

    Lin H, Gao J, Mei Q, et al. Adaptive digital fringe projection technique for high dynamic range three-dimensional shape measurement[J]. Optics Express, 2016, 24(7):7703-7718. doi: 10.1364/OE.24.007703

    [4] 黄亚楠, 娄小平.基于多频外差原理的相位校正及匹配方法研究[J].应用光学, 2014, 35(2):237-241. http://d.old.wanfangdata.com.cn/Periodical/yygx201402012

    Huang Yanan, Lou Xiaoping. Phase correction and matching based on multi-frequency heterodyne method[J]. Journal of Applied Optics, 2014, 35(2):237-241. http://d.old.wanfangdata.com.cn/Periodical/yygx201402012

    [5]

    Almarazcabral C, Gonzalezbarbosa J, Villa J, et al. Fringe projection profilometry for panoramic 3D reconstruction[J]. Optics & Lasers in Engineering, 2016, 78:106-112. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=f8e1dd81ecfac2be4e8944caf98ef6a0

    [6]

    Huang Y, Huang B, Xu H, et al. Moiré fringe method of using warping deformation measurement of electronic components[J]. Proceedings of SPIE: The International Society for Optical Engineering, 2010, 7657:76571K-76571K-6. https://www.researchgate.net/publication/252946811_Moire_fringe_method_of_using_warping_deformation_measurement_of_electronic_components

    [7]

    Shao Xinxing, Dai Xiangjun, Chen Zhenning, et al. Real-time 3D digital image correlation method and its application in human pulse monitoring[J]. Applied Optics, 2016, 55(4):696-704. doi: 10.1364/AO.55.000696

    [8]

    Yang Lianxiang, Zhu Lianqing. Review of electronic speckle pattern interferometry(ESPI) for three dimensional displacement measurement[J]. Chinese Journal of Mechanical Engineering, 2014, 27(1):1-13. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=jxgcxb-e201401001

    [9]

    Wu Sijin, Zhu Lianqing, Pan Siyang, et al. Spatiotemporal three-dimensional phase unwrapping in digital speckle pattern interferometry[J]. Optics Letters, 2016, 41(5):1050-1053. doi: 10.1364/OL.41.001050

    [10]

    Liu Kai, Wu Sijin, Gao Xinya, et al. Simultaneous measurement of in-plane and out-of-plane deformations using dual-beam spatial-carrier digital speckle pattern interferometry[J]. Applied Mechanics & Materials, 2015, 782:316-325. http://cn.bing.com/academic/profile?id=a1556074d21ac829e2a5bd995056e179&encoded=0&v=paper_preview&mkt=zh-cn

    [11]

    Kumar M, Shakher C. Measurement of elastic and thermal properties of composite materials using digital speckle pattern interferometry[J]. Proceeding of SPIE, 2015, 9660(12):4498-510. http://cn.bing.com/academic/profile?id=da61863fa8af669935c5f6f2d5c39df3&encoded=0&v=paper_preview&mkt=zh-cn

    [12]

    Viotti M R, Albertazzi A, Bonomo D, et al. Radial in-plane digital speckle pattern interferometer combined with instrumented indentation[J]. Optics and Lasers in Engineering. 2015, 71:1-8. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=0f6496226dea8f899d45a78f1266441b

  • 期刊类型引用(6)

    1. 于淼,唐小军,杨耀东,吴思进. 基于二维剪切散斑干涉的无人机缺陷检测. 科技创新与应用. 2023(02): 5-10 . 百度学术
    2. 王永红,包凤卿,张肖,赵琪涵,陈维杰,闫佩正. 数字散斑干涉三维变形测量技术研究进展. 应用光学. 2020(04): 681-689 . 本站查看
    3. 王永红,王欢庆,吴双乐,赵琪涵,孙方圆. 基于机器视觉的材料力学性能全场测量方法. 测控技术. 2020(08): 35-46 . 百度学术
    4. 闫浩,隆军,刘驰越,潘淑媛,左超,蔡萍. 数字全息技术及散斑干涉技术在形变测量领域的发展及应用. 红外与激光工程. 2019(06): 154-166 . 百度学术
    5. 吴思进,杨靖,潘思阳,李伟仙,杨连祥. 数字散斑干涉术和时空三维相位解包裹用于非连续表面动态变形测量(英文). 光子学报. 2018(02): 127-135 . 百度学术
    6. 刘承运,吴思进,李伟仙,杨连祥. 应用于数字散斑干涉测量的相移器结构设计. 自动化仪表. 2018(07): 61-64 . 百度学术

    其他类型引用(9)

图(7)
计量
  • 文章访问数:  864
  • HTML全文浏览量:  294
  • PDF下载量:  98
  • 被引次数: 15
出版历程
  • 收稿日期:  2016-07-03
  • 修回日期:  2016-10-27
  • 刊出日期:  2016-12-31

目录

/

返回文章
返回