陳貞屹，趙文川，張啟燦，漢 語(yǔ)，劉元坤

基于立體相位測(cè)量偏折術(shù)的預(yù)應(yīng)力薄鏡面形檢測(cè)

陳貞屹1,2，趙文川2，張啟燦1*，漢 語(yǔ)2，劉元坤1

1四川大學(xué)電子信息學(xué)院，四川 成都 610065；2中國(guó)科學(xué)院光電技術(shù)研究所，四川 成都 610209

應(yīng)力拋光技術(shù)通過(guò)在鏡面上施加預(yù)定載荷，將包括自由曲面在內(nèi)的非球面轉(zhuǎn)化為球面進(jìn)行加工，對(duì)加工鏡面的形變進(jìn)行精準(zhǔn)檢測(cè)是實(shí)現(xiàn)高精度應(yīng)力拋光的關(guān)鍵。利用立體相位測(cè)量偏折術(shù)對(duì)預(yù)應(yīng)力薄鏡進(jìn)行鏡面面形和形變檢測(cè)，獲得被測(cè)鏡表面的連續(xù)相位分布，結(jié)合表面法線唯一性與梯度分布積分，最終得到被測(cè)鏡的高度分布和面形。模擬了系統(tǒng)誤差成分，同時(shí)采用旋轉(zhuǎn)平均法對(duì)系統(tǒng)誤差進(jìn)行標(biāo)定去除，保證和提高了測(cè)量精度。對(duì)一塊口徑320 mm，球面半徑5200 mm的預(yù)應(yīng)力薄鏡面形和變形量進(jìn)行測(cè)量，靜態(tài)測(cè)量結(jié)果與三坐標(biāo)機(jī)測(cè)量結(jié)果對(duì)比，動(dòng)態(tài)應(yīng)變測(cè)量結(jié)果與有限元仿真結(jié)果對(duì)比，分別一致吻合，表明本文方法具備微米級(jí)的測(cè)量精度，相比于干涉儀和三坐標(biāo)機(jī)更適用于大面形變化的預(yù)應(yīng)力薄鏡檢測(cè)。

光學(xué)面形檢測(cè)；立體相位測(cè)量偏折術(shù)；系統(tǒng)誤差；預(yù)應(yīng)力拋光

1 引 言

拼接鏡面主動(dòng)光學(xué)技術(shù)在大口徑望遠(yuǎn)鏡研制過(guò)程中的成功應(yīng)用，幫助天文望遠(yuǎn)鏡突破了口徑限制[1]。然而隨著主鏡口徑的不斷增大，構(gòu)成主鏡所需的離軸子鏡數(shù)量也呈幾何量級(jí)增長(zhǎng)。為實(shí)現(xiàn)離軸子鏡的高效率加工，加州大學(xué)Nelson教授等[2-3]提出了應(yīng)力拋光技術(shù)，該技術(shù)能夠?qū)ㄗ杂汕嬖趦?nèi)的非球面轉(zhuǎn)化為球面進(jìn)行加工，與數(shù)控子孔徑加工技術(shù)相比，材料去除效率大幅提高[4]。

應(yīng)力拋光技術(shù)的關(guān)鍵之一是對(duì)鏡面變形進(jìn)行精準(zhǔn)測(cè)量，現(xiàn)今主要測(cè)量手段包括干涉儀和三坐標(biāo)機(jī)。通常，此類被測(cè)鏡球面曲率半徑較大，長(zhǎng)達(dá)數(shù)米，若使用干涉測(cè)量[5]，干涉儀需放在曲率中心處，測(cè)量距離也達(dá)到數(shù)米，空氣擾動(dòng)將大大影響測(cè)量結(jié)果；另外，預(yù)應(yīng)力薄鏡的變形量通常都在微米量級(jí)以上，超出了干涉儀的動(dòng)態(tài)范圍而無(wú)法測(cè)量。三坐標(biāo)機(jī)測(cè)量屬于接觸式測(cè)量，測(cè)頭需要與鏡面接觸，采樣密度低，測(cè)量速度慢，無(wú)法對(duì)面形進(jìn)行實(shí)時(shí)反饋。

相位測(cè)量偏折術(shù)(Phase measuring deflectometry，PMD)基于結(jié)構(gòu)光三維檢測(cè)技術(shù)，具有系統(tǒng)簡(jiǎn)單、成本低、動(dòng)態(tài)范圍大、靈敏度高、測(cè)量效率高的特點(diǎn)[6-7]，經(jīng)過(guò)高精度的系統(tǒng)標(biāo)定與誤差校準(zhǔn)，可以實(shí)現(xiàn)高精度的面形測(cè)量。同時(shí)，通過(guò)條紋的相對(duì)變形，可以很精確地檢測(cè)到面形的變化，也就是說(shuō)，該技術(shù)可以非常方便容易地測(cè)量面形的相對(duì)變化，非常適合預(yù)應(yīng)力薄鏡的檢測(cè)需求。本文研究將立體相位測(cè)量偏折術(shù)運(yùn)用到預(yù)應(yīng)力薄鏡的面形檢測(cè)中，同時(shí)結(jié)合旋轉(zhuǎn)平均法消除系統(tǒng)誤差影響，實(shí)現(xiàn)了對(duì)一塊口徑320 mm，球面半徑5200 mm的預(yù)應(yīng)力薄鏡初始狀態(tài)面形和面形變形量的快速測(cè)量，結(jié)果與三坐標(biāo)機(jī)以及有限元仿真理論影響力函數(shù)對(duì)比一致。

2 基本原理

測(cè)量系統(tǒng)如圖1所示，由相機(jī)、被測(cè)鏡、顯示屏組成。兩個(gè)相機(jī)通過(guò)被測(cè)鏡的反射，拍攝顯示屏上依次相移的水平和垂直方向正弦條紋，經(jīng)過(guò)相移條紋的對(duì)應(yīng)相位計(jì)算和相位展開(kāi)后，得到攜帶了被測(cè)鏡面形信息的連續(xù)相位分布。按照光線可逆原理，假設(shè)光線從相機(jī)發(fā)出，相機(jī)每個(gè)像素點(diǎn)發(fā)出的光線經(jīng)被測(cè)面反射后最終入射到顯示屏上。在針孔相機(jī)模型下，由連續(xù)相位分布可以得到相機(jī)每個(gè)像素點(diǎn)入射到顯示屏上的坐標(biāo)。

在立體相位測(cè)量偏折術(shù)中，Camera1任意像素點(diǎn)1出射到被測(cè)鏡上的物點(diǎn)坐標(biāo)可通過(guò)假設(shè)來(lái)確定。如圖2所示，1,2分別為Camera1和Camera2的相機(jī)光心，實(shí)際情況下，1對(duì)應(yīng)的物點(diǎn)位置應(yīng)在被測(cè)鏡表面的點(diǎn)，1發(fā)出的光線入射到顯示屏上的點(diǎn)為1，Camera2對(duì)點(diǎn)的成像點(diǎn)為2，2發(fā)出的光線入射到顯示屏上的點(diǎn)為2?，F(xiàn)假設(shè)1點(diǎn)對(duì)應(yīng)的物點(diǎn)位置在￠處，則Camera2的成像點(diǎn)由2變?yōu)椤?，經(jīng)過(guò)￠點(diǎn)的光線與顯示屏的交點(diǎn)由2變?yōu)椤?。由1，￠，1和￠2，￠，￠2的坐標(biāo)可分別算出Camera1，Camera2中￠點(diǎn)的表面法線1，2。根據(jù)法線唯一性原則[8-9]，如果假設(shè)的點(diǎn)￠在真實(shí)物點(diǎn)處，則1應(yīng)該與2相等，否則繼續(xù)假設(shè)物點(diǎn)位置，反復(fù)迭代前述過(guò)程，直到1=2為止，由此點(diǎn)的坐標(biāo)便可找到。

圖1 測(cè)量系統(tǒng)示意圖

圖2 雙目PMD原理示意圖

理論上對(duì)所有點(diǎn)都使用該法線唯一性就可求出整個(gè)被測(cè)鏡的面形和高度，但由于雙目視覺(jué)測(cè)量精度對(duì)噪聲比較敏感，測(cè)得精度有限，計(jì)算速度也非常慢[10]。本文將利用假設(shè)點(diǎn)的法線束來(lái)獲取梯度分布[10]，結(jié)合一個(gè)可靠點(diǎn)的準(zhǔn)確高度，通過(guò)Southwell模型區(qū)域波前重構(gòu)積分[11-13]得到初始面形分布，反復(fù)迭代此過(guò)程得到被測(cè)面真實(shí)面形[14]。由于選取的可靠點(diǎn)往往不是被測(cè)鏡中心點(diǎn)，導(dǎo)致積分的采樣面坐標(biāo)系與物體坐標(biāo)系間有一定旋轉(zhuǎn)，在得到恢復(fù)的矢高面形后，應(yīng)將面形旋轉(zhuǎn)回物體坐標(biāo)系。

從光強(qiáng)中提取的相位信息除了被測(cè)反射鏡的高度調(diào)制相位，還有系統(tǒng)誤差，因此最后恢復(fù)的反射鏡面形也包含了系統(tǒng)誤差。這些誤差來(lái)自于條紋屏平面性、相機(jī)噪聲、計(jì)算積分時(shí)誤差擴(kuò)散和系統(tǒng)標(biāo)定引入的誤差等[15]，其中系統(tǒng)標(biāo)定的誤差是主要影響因素[16]。如何去除這些系統(tǒng)誤差將取決于其組成成份。為了檢驗(yàn)和分析此系統(tǒng)中標(biāo)定不準(zhǔn)確帶來(lái)的系統(tǒng)誤差成份，根據(jù)標(biāo)定的系統(tǒng)位置參數(shù)建立系統(tǒng)模型[17]，并模擬了系統(tǒng)誤差。圖3為標(biāo)定誤差較大時(shí)的系統(tǒng)誤差分布，表1為系統(tǒng)誤差的Zernike多項(xiàng)式分解。表1的結(jié)果證明了PMD中系統(tǒng)誤差主要是低階非旋轉(zhuǎn)對(duì)稱部分[18]。因此，本文采用旋轉(zhuǎn)平均法[19-21]來(lái)標(biāo)定和去除系統(tǒng)誤差，有效保證和提升測(cè)量精度。

圖3 模擬系統(tǒng)誤差分布

3 實(shí) 驗(yàn)

實(shí)驗(yàn)裝置由2個(gè)CCD相機(jī)(PointGrey, GS3-U3-28S5M-C，分辨率為1920 pixels′1440 pixels，像素尺寸為4.54 μm)、一個(gè)80寸(1寸=2.54 cm)液晶顯示屏(Sharp, LCD-80X818A，分辨率為3840 pixels′2160 pixels)、被測(cè)彎月型薄鏡及其支撐裝置和計(jì)算機(jī)構(gòu)成，如圖4所示。所使用的2個(gè)CCD相機(jī)分上下放置在顯示器同一側(cè)，取下方相機(jī)為主相機(jī)。被測(cè)預(yù)應(yīng)力薄鏡系統(tǒng)如圖5(a)所示，口徑為320 mm，其球面半徑為5200 mm，被固定在一個(gè)轉(zhuǎn)臺(tái)上，可靈活旋轉(zhuǎn)。薄鏡后表面中心有一個(gè)固定支撐點(diǎn)，外沿均勻分布6個(gè)機(jī)電式力促動(dòng)器(分別記為1,2,3,4,5,6號(hào)電機(jī))，如圖5(b)所示。被測(cè)薄鏡位置距離顯示屏約2000 mm左右。

測(cè)量前首先進(jìn)行系統(tǒng)標(biāo)定，先標(biāo)定2個(gè)相機(jī)的內(nèi)參[22-23]并計(jì)算彼此相對(duì)位置[24]，再用標(biāo)準(zhǔn)靶面當(dāng)作反射鏡標(biāo)定屏幕與主相機(jī)的相對(duì)位置[24-25]。圖6是系統(tǒng)位置參數(shù)標(biāo)定結(jié)果示意圖，1號(hào)靶面即為標(biāo)定相機(jī)與屏幕關(guān)系的標(biāo)定靶，標(biāo)定雙目系統(tǒng)的靶面位置在圖中已省略。

圖4 測(cè)量裝置圖

表1 系統(tǒng)誤差的Zernike系數(shù)

圖5 被測(cè)預(yù)應(yīng)力薄鏡。(a) 薄鏡旋轉(zhuǎn)支撐結(jié)構(gòu)側(cè)視圖；(b) 薄鏡背面支撐點(diǎn)位置分布示意圖

圖6 系統(tǒng)位置參數(shù)標(biāo)定結(jié)果示意圖

3.1 初始狀態(tài)面形測(cè)量

測(cè)量中采用了時(shí)間相位展開(kāi)算法[26-27]進(jìn)行相位計(jì)算，拍攝變形條紋時(shí)在顯示器上分別用不同頻率、水平和豎直不同方向上的各3組條紋圖相移4次，拍攝耗時(shí)約40 s，圖7為上下相機(jī)拍攝的一幀圖像。使用了旋轉(zhuǎn)平均法進(jìn)行系統(tǒng)誤差標(biāo)定，每隔60°轉(zhuǎn)動(dòng)薄鏡一次，一共轉(zhuǎn)動(dòng)6次完成一周360°的測(cè)量，在每個(gè)旋轉(zhuǎn)位置上分別測(cè)量3次，求平均值作為該旋轉(zhuǎn)度數(shù)位置下的測(cè)量面形，6個(gè)度數(shù)下測(cè)量的矢高如圖8所示。

圖9為旋轉(zhuǎn)平均法求出的系統(tǒng)誤差，圖10為去系統(tǒng)誤差、平移、傾斜和旋轉(zhuǎn)對(duì)稱項(xiàng)后的鏡面面形，其PV=4.53 μm，RMS=0.754 μm。圖11為三坐標(biāo)機(jī)的測(cè)量數(shù)據(jù)除去平移、傾斜、旋轉(zhuǎn)對(duì)稱部分后的面形誤差分布(PV=7.879 μm, RMS=1.163 μm)。對(duì)比兩圖可以看出，本文方法的測(cè)量結(jié)果與三坐標(biāo)機(jī)檢測(cè)結(jié)果分布大體一致。需要說(shuō)明的是，在進(jìn)行三坐標(biāo)測(cè)量時(shí)，預(yù)應(yīng)力薄鏡鏡面豎直向上，而在本系統(tǒng)測(cè)量中，預(yù)應(yīng)力薄鏡鏡面是水平方向。

3.2 變形量測(cè)量

預(yù)應(yīng)力薄鏡背面的力促動(dòng)電機(jī)可對(duì)鏡體施加雙向載荷，除了每個(gè)電機(jī)獨(dú)立工作對(duì)鏡面產(chǎn)生變形以外，6個(gè)電機(jī)還可以同時(shí)施加不同的校正力，產(chǎn)生球差、像散、慧差、三葉草像差?，F(xiàn)將6個(gè)力促動(dòng)電機(jī)分別單獨(dú)施加5 N的力向外頂，通過(guò)對(duì)變形前后兩次面形測(cè)量結(jié)果相減，得到各電機(jī)單獨(dú)加力后的鏡面變形量。圖12分別是ANSYS仿真的1~6號(hào)促動(dòng)器施加5 N外力的理論影響力函數(shù)分布，圖13為實(shí)測(cè)結(jié)果。圖14為理論分析與實(shí)際測(cè)量的PV、RMS折線圖。

圖7 顯示豎直方向相移條紋時(shí)拍攝的條紋圖。(a) 下相機(jī)圖拍攝；(b) 上相機(jī)圖拍攝

圖8 6個(gè)旋轉(zhuǎn)角度下的薄鏡表面矢高測(cè)量結(jié)果。(a) 0°；(b) 60°；(c) 120°；(d) 180°；(e) 240°；(f) 300°

圖9 系統(tǒng)誤差分布

圖10 去除系統(tǒng)誤差后的鏡面面形

圖11 三坐標(biāo)機(jī)測(cè)量的薄鏡面形結(jié)果

從圖中可以看出，理論上單個(gè)電機(jī)施加5 N的力能造成鏡面約1.68 μm的變形量，與實(shí)際測(cè)量結(jié)果一致，幅度和分布也一致。各電機(jī)分布中心處的低點(diǎn)是因?yàn)殓R面后表面中心為固定支撐點(diǎn)時(shí)，邊緣促動(dòng)電機(jī)外頂造成的；加力電機(jī)另一邊的鏡面突起也是由于中心被固定后，另一邊被作用點(diǎn)一同帶起。

接著對(duì)預(yù)應(yīng)力薄鏡的像差變形進(jìn)行了測(cè)量，四種像差的理論目標(biāo)校正函數(shù)如圖15所示，對(duì)應(yīng)變形后的檢測(cè)結(jié)果如圖16所示，表2為理論目標(biāo)校正函數(shù)與實(shí)際變形檢測(cè)結(jié)果的PV、RMS對(duì)比。圖表展示的結(jié)果表明，4種像差對(duì)應(yīng)薄鏡變形的本文檢測(cè)結(jié)果和相應(yīng)的目標(biāo)校正函數(shù)是一致的。

4 結(jié) 論

本文充分利用偏折術(shù)測(cè)量速度快、動(dòng)態(tài)范圍大和靈敏度高的特點(diǎn)，采用立體相位測(cè)量偏折術(shù)，對(duì)預(yù)應(yīng)力薄鏡的初始面形和應(yīng)變量進(jìn)行了檢測(cè)。與三坐標(biāo)機(jī)檢測(cè)方法相比，該方法結(jié)構(gòu)簡(jiǎn)單，可以一次完成非接觸全場(chǎng)測(cè)量，可以對(duì)面形變形進(jìn)行快速反饋。與干涉方法相比，該方法有更大的測(cè)量動(dòng)態(tài)范圍，受環(huán)境振動(dòng)和氣流影響更小。這些特點(diǎn)使該方法非常適用于大口徑、大變形量的自由面應(yīng)力薄鏡檢測(cè)。文中模擬了系統(tǒng)誤差成分，并通過(guò)旋轉(zhuǎn)平均去除系統(tǒng)誤差，保證和提高了檢測(cè)精度。本文檢測(cè)結(jié)果與三坐標(biāo)機(jī)檢測(cè)結(jié)果的細(xì)微差別主要是被測(cè)鏡體放置方式的不同，導(dǎo)致重力影響的不同而引起的。實(shí)驗(yàn)表明，本文方法的檢測(cè)結(jié)果與三坐標(biāo)機(jī)測(cè)量結(jié)果，以及促動(dòng)器的理論目標(biāo)函數(shù)結(jié)果相一致。

圖12 各個(gè)促動(dòng)電機(jī)單獨(dú)施力的理論影響力函數(shù)。(a) 1號(hào)；(b) 2號(hào)；(c) 3號(hào)；(d) 4號(hào)；(e) 5號(hào)；(f) 6號(hào)

圖13 各個(gè)促動(dòng)電機(jī)單獨(dú)施力的變形量測(cè)量結(jié)果。(a) 1號(hào)；(b) 2號(hào)；(c) 3號(hào)；(d) 4號(hào)；(e) 5號(hào)；(f) 6號(hào)

圖14 電機(jī)理論變形與實(shí)測(cè)變形對(duì)比。(a) PV；(b) RMS

圖15 4種像差的理論目標(biāo)校正函數(shù)。(a) 球差；(b) 像散；(c) 慧差；(d) 三葉草像差

圖16 4種像差的變形量實(shí)測(cè)結(jié)果。(a) 球差；(b) 像散；(c) 慧差；(d) 三葉草像差

表2 像差理論值與實(shí)際測(cè)量值對(duì)比

本文拓展了相位測(cè)量偏折術(shù)的應(yīng)用領(lǐng)域，將其應(yīng)用于預(yù)應(yīng)力薄面形檢測(cè)中，下一步可對(duì)鏡面面形進(jìn)行實(shí)時(shí)監(jiān)測(cè)，對(duì)力促動(dòng)器的性能和效果進(jìn)行實(shí)時(shí)測(cè)量反饋，提高校正的動(dòng)態(tài)能力。另外，在后續(xù)工作中，可在標(biāo)定顯示屏與主相機(jī)位置關(guān)系時(shí)使用平面性更好的標(biāo)準(zhǔn)反射鏡，或者用三坐標(biāo)機(jī)先對(duì)顯示屏的平面性進(jìn)行檢測(cè)，消除非平面性引入的誤差，以進(jìn)一步提高測(cè)量精度。

[1] Li X N, Jiang Z B, Gong X F,Stressed mirror annular polishing for scale-down TMT primary segments[J]., 2016, 9912: 99120A.

[2] Lubliner J, Nelson J E. Stressed mirror polishing. 1: A technique for producing nonaxisymmetric mirrors[J]., 1980, 19(14): 2332–2340.

[3] Nelson J E, Gabor G, Hunt L K,Stressed mirror polishing. 2: fabrication of an off-axis section of a paraboloid[J]., 1980, 19(14): 2341–2352.

[4] Li X N, Zhang H Y, Cui X Q,Study on stressed mirror polishing with continoues polishing machine for large aperture off-axis aspheric mirror[J]., 2012, 53(2): 161–170.

李新南, 章海鷹, 崔向群, 等. 大口徑離軸非球面鏡預(yù)應(yīng)力環(huán)拋方法研究[J]. 天文學(xué)報(bào), 2012, 53(2): 161–170.

[5] Valente M, Lewis B, Melena N,Advanced surface metrology for meter-class optics[J]., 2013, 8838: 88380F.

[6] Huang L, Idir M, Zuo C,Review of phase measuring deflectometry[J]., 2018, 107: 247–257.

[7] Zhao W C, Zhou M, Liu H T,. The off-axis aspheric mirror testing based on the fringe reflection technique[J]., 2018, 45(7): 170663.

趙文川, 周敏, 劉海濤, 等. 離軸非球面的條紋反射檢測(cè)技術(shù)[J]. 光電工程, 2018, 45(7): 170663.

[8] Liu F M, Lin J R, Sun Y B,Calibration method for stereo phase measuring deflectometry system[J]., 2020, 57(5): 051202.

劉方明, 林嘉睿, 孫巖標(biāo), 等. 一種立體相位偏折測(cè)量系統(tǒng)標(biāo)定方法[J]. 激光與光電子學(xué)進(jìn)展, 2020, 57(5): 051202.

[9] Knauer M C, Kaminski J, H?usler G. Phase measuring deflectometry: a new approach to measure specular free-form surfaces[J]., 2004, 5457: 366–376.

[10] 毛心潔. 基于雙目立體視覺(jué)的鏡面物體三維輪廓測(cè)試技術(shù)[D]. 南京: 南京理工大學(xué), 2017: 48–49.

[11] Southwell W H. Wave-front estimation from wave-front slope measurements[J]., 1980, 70(8): 998–1006.

[12] Zhu Y J, Zhu L X, Zhong J P,Integration algorithm in three-dimensional plane Shape measuring by fringe reflection[J]., 2018, 47(11): 110–119.

朱勇建, 朱立新, 鐘建平, 等. 條紋反射法測(cè)量三維面形中的積分算法[J]. 光子學(xué)報(bào), 2018, 47(11): 110–119.

[13] Yi J Y. Study on mobile phone shell inside and outside surface quality inspection based on fringe projection and fringe reflection technologies[D]. Chengdu: University of Electronic Science and Technology of China, 2016: 41–43.

易京亞. 基于條紋投影和條紋反射的手機(jī)殼內(nèi)外表面質(zhì)量檢測(cè)方法[D]. 成都: 電子科技大學(xué), 2016: 41–43.

[14] Zhao W C, Graves L R, Huang R,Iterative surface construction for blind deflectometry[J]., 2016, 9684: 96843X.

[15] Zhao W C, Su X Y, Liu Y K,Testing an aspheric mirror based on phase measuring deflectometry[J]., 2010, 37(5): 1338–1341.

趙文川, 蘇顯渝, 劉元坤, 等. 基于相位偏折術(shù)的非球面鏡檢測(cè)方法[J]. 中國(guó)激光, 2010, 37(5): 1338–1341.

[16] Faber C, Olesch E, Krobot R,Deflectometry challenges interferometry: the competition gets tougher![J]., 2012, 8493: 84930R.

[17] Korsch D, Hunter W R.[M]. Academic Press, 1991.

[18] Su P, Khreishi M, Huang RPrecision aspheric optics testing with SCOTS: a deflectometry approach[]., 2013, 8788: 87881E.

[19] Song W H, Wu F, Hou X. Method to test rotationally asymmetric surface deviation with high accuracy[J]., 2012, 51(22): 5567–5572.

[20] Song W, Wu F, Hou X,. Absolute calibration of a spherical reference surface for a Fizeau interferometer with the shift-rotation method of iterative algorithm[J]., 2013, 52(3): 033601.

[21] Yang P, Wu F, Hou X. Simulation analysis of absolute measurement for rotationally asymmetric surface error[J]., 2011, 38(1): 93–97, 102.

楊鵬, 伍凡, 侯溪. 旋轉(zhuǎn)非對(duì)稱項(xiàng)面形誤差絕對(duì)檢測(cè)的仿真分析[J]. 光電工程, 2011, 38(1): 93–97, 102.

[22] Zhang Z. A flexible new technique for camera calibration[J]., 2000, 22(11): 1330–1334.

[23] Gai Q Y. Optimization of stereo matching in 3D reconstruction based on binocular vision[J]., 2018, 960: 012029.

[24] Liu Y K, Su X Y. Camera calibration with planar crossed fringe patterns[J]., 2012, 123(2): 171–175.

[25] Zhou T, Chen K, Wei H Y,. Improved system calibration for specular surface measurement by using reflections from a plane mirror[J]., 2016, 55(25): 7018–7028.

[26] Saldner H O, Huntley J M. Temporal phase unwrapping: application to surface profiling of discontinuous objects[J]., 1997, 36(13): 2770–2775.

[27] Zhao W J, Chen W J, Su X Y. The comparison of several time phase unwrapping methods[J]., 2016, 53(1): 110–117.

趙文靜, 陳文靜, 蘇顯渝. 幾種時(shí)間相位展開(kāi)方法的比較[J]. 四川大學(xué)學(xué)報(bào)(自然科學(xué)版), 2016, 53(1): 110–117.

Shape measurement of stressed mirror based on stereoscopic phase measuring deflectometry

Chen Zhenyi1,2, Zhao Wenchuan2, Zhang Qican1*, Han Yu2, Liu Yuankun1

1School of Electronics and Information Engineering, Sichuan University, Chengdu, Sichuan 610065, China;2Institute of Optics and Electronics, Chinese Academy of Sciences, Chengdu, Sichuan 610209, China

Measuring device setup

Overview:Stressed polishing technology was firstly proposed by Jerry E.Nelson in 1980, considered as an effective method of aspheric fabrication. According to the calculated results of elastic mechanics, this method exerts an external force on the mirror to form a deformation which opposite to the desired deformation from spherical surface to off-axis aspheric surface. Under the state of deformation, the mirror is polished into sphere surface. After removing the external force, the required off-axis aspheric surface can be obtained. Because aspheric fabrication is converted to spherical fabrication, tools with large diameter can be used and efficiency is greatly improved.

The key to achieve high precision of stressed polishing is to test the deformation of mirror with high precision. However, the main methods of surface measurement nowadays are interferometer and CMM. If interferometer is used, its dynamic range can only support the detection below micron deformation. If CMM is used, probe may scratch the mirror surface, and the detection tempo is very slow. Furthermore, interferometer and CMM are both expensive and complex equipments.

So, aimed at stressed polishing above micron deformation, stereoscopic phase measuring deflectometry was used to test its surface topography and deformation. It is low cost and convenient technique and only screen, camera, and computer were needed when implemented. More importantly, characteristics such as high dynamic range, full-field three-dimensional measurement and excellent performance in medium and high frequency were brought in, which are very suitable for the test of stressed mirror.

When measuring, firstly calculating the unwrapped phase distribution through CCD cameras, then calculating the height of a specific point on the measured surface using normal consistency constraint, and finally the full-field height distribution was obtained by Southwell gradient integral algorithm. To improve the measuring accuracy, composition of systematic errors were simulated, proved that it mainly includes low-order non-rotational symmetry items. According to simulating results, errors were calibrated and removed by-step averaging method to get a absolute surface topography.

In this paper, the absolute surface topography and the deformation of a stressed mirror with a diameter of 320 mm, spherical radius of 5200 mm were measured. The measuring results were consistent with the corresponding result of CMM and finite element simulation, indicating that this proposed method is on the level of micron in terms of accuracy and more suitable for the test of stressed mirror compared with interferometer and CMM.

Citation: Chen Z Y, Zhao W C, Zhang Q C,Shape measurement of stressed mirror based on stereoscopic phase measuring deflectometry[J]., 2020, 47(8): 190435

Shape measurement of stressed mirror based on stereoscopic phase measuring deflectometry

Chen Zhenyi1,2, Zhao Wenchuan2, Zhang Qican1*, Han Yu2, Liu Yuankun1

1School of Electronics and Information Engineering, Sichuan University, Chengdu, Sichuan 610065, China;2Institute of Optics and Electronics, Chinese Academy of Sciences, Chengdu, Sichuan 610209, China

Stressed polishing technology transforms aspheric fabrication into spherical fabrication by applying predetermined loads on the surface of the mirror. The key to achieve high precision of stressed polishing is to test the surface deformation with high precision. Stereoscopic phase measuring deflectometry was used to test the surface topography and the deformation of stressed mirror. After obtained unwrapped phase distribution, and combined with normal consistency constraint and gradient integral algorithm, the height distribution was finally obtained. Composition of systematic errors were simulated. Also, the errors were calibrated and removed by-step averaging method in this system, which improved the measuring precision. In this paper, the surface topography and the deformation of a stressed mirror with a diameter of 320 mm, spherical radius of 5200 mm were measured. The measuring results were consistent with the corresponding result of CMM and finite element simulation, indicating that this proposed method is on the level of micron in terms of accuracy and more suitable for the test of stressed mirror compared with interferometer and CMM.

optical testing; stereoscopic phase measuring deflectometry; systematic error; stressed mirror polishing

TN247；TH741

A

10.12086/oee.2020.190435

: Chen Z Y, Zhao W C, Zhang Q C,. Shape measurement of stressed mirror based on stereoscopic phase measuring deflectometry[J]., 2020,47(8): 190435

陳貞屹，趙文川，張啟燦，等. 基于立體相位測(cè)量偏折術(shù)的預(yù)應(yīng)力薄鏡面形檢測(cè)[J]. 光電工程，2020，47(8): 190435

Supported by NationalNaturalScienceFoundation of China (61675141, 61505216) and The Youth Innovation Promotion Association, Chinese Academy of Sciences

* E-mail: zqc@scu.edu.cn

2019-07-23；

2019-10-22

國(guó)家自然科學(xué)基金資助項(xiàng)目(61675141，61505216)；中國(guó)科學(xué)院青年創(chuàng)新促進(jìn)會(huì)資助

陳貞屹(1995-)，男，碩士研究生，主要從事光學(xué)元件三維形貌測(cè)量的研究。E-mail：13628045693@163.com

張啟燦(1974-)，男，博士，教授，主要從事三維傳感、動(dòng)態(tài)三維測(cè)量等研究。E-mail：zqc@scu.edu.cn