張德海，白代萍，楊勇，王良文，郭成

(1.鄭州輕工業(yè)學(xué)院機電工程學(xué)院，河南鄭州450002;2.西安交通大學(xué)機械工程學(xué)院，陜西西安710049)

冰箱門外殼多數(shù)屬于大曲率半徑?jīng)_壓件，回彈是其成形過程中存在的主要問題之一。而出于美觀化和實用化的考慮，VCM 鋼板和噴涂鋼板得到大量采用，其高塑性或小彈性模量使得回彈問題進一步突出，從成形工藝上進行控制是目前常用的方法，如對壓邊力、拉延筋力等的調(diào)整［1－2］?；跀?shù)值模擬的回彈補償算法，一般須借助有限元軟件，利用設(shè)計好的模面與毛坯建立相應(yīng)的數(shù)值模型，并進行成形回彈分析，再將回彈誤差量反向補償?shù)侥＞呔W(wǎng)格上。這類方法具有很好的通用性，但沒有考慮到制造工程多品種小批量的具體特點［3］。

三維光學(xué)檢測技術(shù)作為一種新興的非接觸測量方法，隨著計算機技術(shù)和光學(xué)技術(shù)的成熟而得到了迅猛的發(fā)展。該檢測方法對于大型或超大型工件、小型或超小型工件，尤其是具有大量A 級曲面或A 級曲線的不規(guī)則產(chǎn)品、非標(biāo)產(chǎn)品等，具有無可比擬的科學(xué)優(yōu)勢，并呈現(xiàn)出了強大的生命力，日益應(yīng)用到產(chǎn)品幾何尺寸檢測、逆向工程、材料力學(xué)性能測試和在線質(zhì)量檢測等方面。

作者對三維光學(xué)檢測技術(shù)進行了深入研究，并把該方法率先應(yīng)用到冰箱門外殼的成形回彈測試上，通過對實際零件測試獲得點云數(shù)模，與設(shè)計的CAD 數(shù)模進行基準(zhǔn)對齊后，根據(jù)實際需要可以在任意位置作截面線，獲得兩個數(shù)模之間的差值，即回彈值，包含方向和角度兩個方面。該回彈值可以對材料的成形性能進行分析，是光學(xué)檢測技術(shù)在家電行業(yè)尤其在薄板成形方面的積極嘗試和努力，應(yīng)用的成功將極大推動光學(xué)檢測技術(shù)的發(fā)展。

1 三維光學(xué)檢測技術(shù)

1.1 攝影測量技術(shù)和計算機雙目視覺技術(shù)相結(jié)合的測量方法

基于攝影測量技術(shù)，結(jié)合圖1，放置編碼點和非編碼點在測量對象周圍，沿圓周等分為8 個攝站［18］，每個攝站相機自身旋轉(zhuǎn)0°和90°各拍攝1 張照片，在測量對象正上方拍攝4 張照片。照片輸入攝影測量系統(tǒng)A，設(shè)置相機型號規(guī)格、編碼點類型、標(biāo)尺長度、溫度等參數(shù)，計算出非編碼點和編碼點的三維空間坐標(biāo)。利用輸出功能導(dǎo)出全局非編碼點空間坐標(biāo)矩陣θ(A)，將θ(A)導(dǎo)入以計算機雙目視覺為基礎(chǔ)的面掃描系統(tǒng)B，B 系統(tǒng)掃描一次獲得的點云矩陣要包含不少于5 個非編碼標(biāo)志點的點云矩陣θ(SUBA)，這里θ(SUBA)?θ(A)。

面掃描系統(tǒng)B 計算θ(SUBA)中非編碼點矩陣中的幾何信息，對比θ(A)中的非編碼點矩陣，逐一搜索匹配，直到找出θ(A)和θ(SUBA)對應(yīng)的同名非編碼點。根據(jù)同名非編碼點在攝影測量坐標(biāo)系和面掃描坐標(biāo)系的坐標(biāo)信息，應(yīng)用最小二乘算法獲得轉(zhuǎn)換矩陣ρ(A)，面掃描獲得的表面點坐標(biāo)矩陣和ρ(A)相乘，實現(xiàn)了掃描的所有點云坐標(biāo)矩陣轉(zhuǎn)換，并統(tǒng)一歸并到攝影測量的坐標(biāo)系中。多次拍攝獲得的點云矩陣?yán)萌址蔷幋a標(biāo)志的自動匹配，實現(xiàn)了拍攝點云的自動拼接，利用逆向后處理軟件Geomagic Studio 等獲得最終的單副點云。

通過點云獲取的CAD 數(shù)?？梢宰鳛楣ぜ貜椇蟮男螤?，原有的成形模具即設(shè)計的模具CAD 數(shù)模，兩者導(dǎo)入Geomagic Qualify 逆向軟件，經(jīng)“3－2－1”等操作對齊以后，根據(jù)工程需要可以在任意點作任意方向的截面，獲得兩條截線，截線上相應(yīng)關(guān)鍵點的差值就是各點的回彈量，包含回彈大小和角度兩個方面。

圖1 攝影測量相機攝站布局

1.2 攝影測量技術(shù)

利用立體像對中兩張像片的內(nèi)、外方位元素和像點坐標(biāo)計算對應(yīng)點的三維坐標(biāo)的方法，稱為立體像對的空間前方交會。相機在兩相鄰攝站S1、S2分別拍攝一張像片，構(gòu)成立體像對。任意一點A 在左、右像片的構(gòu)像分別為a1和a2。為了確定像點與物體點的數(shù)學(xué)關(guān)系，建立物體攝影測量坐標(biāo)系D－X1PY1PZ1P，Z1P軸。過左攝站S1建立與物體攝影測量坐標(biāo)系平行的像空間輔助坐標(biāo)系S1－X1Y1Z1，過右攝站S2建立與D－X1PY1PZ1P平行的像空間輔助坐標(biāo)系S2－X2Y2Z2。

設(shè)待測點A 在物體攝影測量坐標(biāo)系D－X1PY1PZ1P中的坐標(biāo)為(XA，YA，ZA)，對應(yīng)像點a1、a2在各自的像空間坐標(biāo)系中的坐標(biāo)為(x1，y1，－ f)和(x2，y2，－f)，在像空間輔助坐標(biāo)系中的坐標(biāo)分別為(X1，Y1，Z1)和(X2，Y2，Z2)。若已知兩張像片的外方位元素，就可以由像點的像空間坐標(biāo)計算出該點的像空間輔助坐標(biāo)，即

式中:R1、R2為由已知的外方位角元素算得的左、右像片的旋轉(zhuǎn)矩陣。

右攝站S2在S1－X1Y1Z1中的坐標(biāo)，即攝影基線B 的3 個分量BX、BY、BZ可由外方位直線元素計算得到:

（1）關(guān)于水流控制與增設(shè)臨時性的水利施工設(shè)施。從水流管理以及施工導(dǎo)流操作的視角來看，對于水流控制通常來講需要關(guān)注如下要點：首先需要妥善排除基坑內(nèi)的積水，并且做好截留堤的全面修建；其次應(yīng)當(dāng)妥善規(guī)劃下游以及上游的施工圍堰，確保能夠有序搭建并且拆除各區(qū)域的導(dǎo)流洞以及固有的圍堰設(shè)施；再次應(yīng)當(dāng)保障下游水流供應(yīng)并且確保順利完成下閘蓄水，此外還要全面做好注漿與封鎖導(dǎo)流洞的有關(guān)操作。

因左、右像空間輔助坐標(biāo)系與物體攝影測量坐標(biāo)系相互平行，且攝站點、像點、物體點三點共線，則由圖2 可得:

式中:N1、N2為左、右像點的點投影系數(shù)。

一般情況下，不同的點有不同的點投影系數(shù)值。根據(jù)式(3)可以得到用前方交會法計算出的物體點坐標(biāo)的公式，即，

圖2 空間前方交會

式(4)中N1、N2仍然未知，為此，結(jié)合式(2)有，

由式(5)中的一、三式聯(lián)立求解得，

式(4)和式(6)就是利用立體像對，在已知像片外方位元素的前提下，由像點坐標(biāo)計算待測點空間坐標(biāo)的前方交會公式。根據(jù)空間后方交會法分別計算出左、右兩張像片的外方位元素(XS1，YS1，ZS1，φ1，ω1，κ1)和(XS2，YS2，ZS2，φ2，ω2，κ2)。

按式(7)分別計算左、右像片的方向余弦，組成左、右像片各自的旋轉(zhuǎn)矩陣R1和R2;結(jié)合左、右像片的外方位直線元素，按式(2)計算攝影基線分量BX、BY、BZ;結(jié)合式(1)逐點計算各像點的像空間輔助坐標(biāo)(X1，Y1，Z1)和(X2，Y2，Z2);結(jié)合式(6)和式(4)可以逐點計算出點投影系數(shù)和待測物體的幾何信息點坐標(biāo)。

1.3 計算機雙目視覺

計算機雙目視覺測量是建立在計算機視覺理論上的三維測量方法。它利用兩臺位置相對固定的攝像機，從不同角度同時獲取同一景物的兩幅圖像，通過計算空間點在兩幅圖像中的像差來獲得其三維坐標(biāo)值。

假設(shè)投影柵光強正弦分布通過光學(xué)透鏡，光柵被投影至物體表面(圖3)，由兩個CCD 攝像頭組成成像系統(tǒng)，在CCD 耙面記錄到物體表面及光柵的影子。由于光強規(guī)則分布的光柵被物體表面調(diào)制而畸變，此畸變了的光柵帶有物體表面形狀信息。借助CCD 耙面上的點和物面、參考面的幾何關(guān)系，直接獲取空間的位置，然后重構(gòu)物體的三維表面形狀。處理像面上“點”的位置反映的光學(xué)圖像的相位即是整個測量的自動化過程。

圖3 編碼光和雙目視覺相結(jié)合的三維光學(xué)測量原理

2 實驗

以國內(nèi)某著名家電公司的BC-48 冰箱門外殼為研究對象，厚度0.6 mm。成形后噴塑保證外觀。應(yīng)用以上所述的檢測方法，可以獲得冰箱門外殼成形后的幾何形狀，其截面形狀和設(shè)計的模具CAD 數(shù)模進行對齊后對比，獲得兩者不同部位、不同方向的三維回彈量數(shù)值。

實驗步驟如下:

(1)攝影測量獲得物體的編碼點和非編碼點，如圖4所示。

(2)計算機雙目視覺技術(shù)獲得回彈后工件點云。

(3)在Geomagic Qualify 中將獲得的數(shù)模和原有的模具CAD 數(shù)模對齊;

(4)2D 比較，采用截線方法，在Geomagic Qualify 中垂直于x 方向等間距做平行截面，間距80 mm，選取6 個截面(圖5)。

圖4 照片導(dǎo)入XJTUDP 計算編碼點和非編碼點

圖5 等間距80 mm 做6 條平行截面

3 討論和分析

對于零件的一些重要截面，對回彈有嚴(yán)格要求。這就需要得到更加精確的回彈量，以便作出合適的回彈補償，保證截面的最終形狀滿足精度的要求。2D比較功能可以對零件的一些重要截面進行回彈檢測。

圖6 為第一個截面獲得的冰箱門外殼回彈數(shù)據(jù)，在Geomagic Qualify 中截取的截面上，實際表現(xiàn)為兩條截線。選取回彈敏感的測試點，沿著z 方向做回彈長度方向的回彈數(shù)值。對于角部回彈，選取兩條截線的切線，取其切線的夾角獲取。該組數(shù)值來自截面1上不同部位關(guān)鍵點的回彈量，讀取后獲得相應(yīng)的數(shù)值如圖7。操作時首先選擇所需要的邊界，可從3 個方向比較出邊界的偏差，分別為法向、切向和3D 偏差。其分別對應(yīng)于CAD 模型邊界處的法線方向偏差、CAD 模型邊界處的切線方向偏差和沿著指定邊界CAD 模型和點云數(shù)據(jù)的最小偏差。這里的回彈長度方向數(shù)值來源于邊界法向，回彈角度數(shù)值來源于邊界切線的延長線夾角。

圖6 CAD 數(shù)模和點云對齊后的回彈數(shù)值顯示(以第一個截面為例)

由圖7 可以看出，回彈最大值D1為7.346 mm 位于中間圓弧處，最小值D9為1.426 mm 位于下圓角R16處沿法線方向的回彈?；貜椊侵械淖畲笾礑13為13.791°位于上部圓角小圓弧和中間大圓弧相結(jié)線處，最小值D4為1.739°位于截面中間大圓弧和下部圓角接近部位。說明該方法可以實現(xiàn)回彈數(shù)值大小和回彈角度的檢測。

圖7 冰箱門外殼截面1 關(guān)鍵點回彈數(shù)據(jù)(1 到7僅僅代表回彈數(shù)值大小的降序排列)

4 結(jié)論

(1)采用攝影測量和計算機雙目視覺技術(shù)可以用來獲取冰箱門外殼的回彈量檢測，包含回彈大小和角度兩個方面;

(2)點云數(shù)模和CAD 數(shù)模對齊后做截線，獲得是二維回彈值，通過在不同方位做截線，可以獲得三維的回彈值;

(3)截面一中回彈最大值為7.346 mm 位于中間圓弧處，最小值1.426 mm 位于下圓角R16 處沿法線方向的回彈?；貜椊侵械淖畲笾禐?3.791°位于上部圓角小圓弧和中間大圓弧相結(jié)線處，最小值1.739°位于截面中間大圓弧和下部圓角接近部位;

(4)文中只給出了截線1 若干關(guān)鍵點的回彈數(shù)值，應(yīng)用提出的檢測方法進行了一定的分析，實際應(yīng)用時用戶可以根據(jù)工程需要適當(dāng)選取。通過和CAD軟件結(jié)合，獲得的三維數(shù)值為沖壓件回彈前后的形狀偏差，如何建立補償回彈變形的模具型面方法，并得到補償和控制回彈的合理沖壓模具型面，是光學(xué)檢測方法能否正確合理運用的根本。

【1】WANG H，LI G Y，ZHONG Z H.Optimization of Sheet Metal Forming Processes by Adaptive Response Surface Based on Intelligent Sampling Method［J］.Journal of Materials Processing Technology，2008，197(1/2/3):77－88.

【2】OHATA T，NAKAMURA Y，KATAYAMA T，et al.Development of Optimum Process Design System for Sheet Fabrication Using Response Surface Method［J］.Journal of Materials Technology，2003，143/144:667－672.

【3】楊偉俊，李東升，李小強，等.基于快速回彈補償?shù)南鹌つ乙簤撼尚文Ｃ嬖O(shè)計方法［J］.機械工程學(xué)報，2011，47(12):67－73.

【4】YOON H J，SONG K Y，KIM J S，et al.Longitudinal Strain Monitoring of Rail Using Adistributed Fiber Sensor Based on Brillouin Optical Correlation Domain Analysis［J］.NDT＆ E International，2011，44:637－644.

【5】TAKEDA M，MUTOH K.Fourier Transform Profilometry for the Automatic Measurement of 3-D Object Shapes［J］.Appl.Optics，1983，22:3977－3982.

【6】KHETAN R P.PhD Dissertation on Projection Moire［D］.New York，State University，1975，USA.

【7】GEISSLER P，DIERIG T，MALLOT H A.Three Dimensional Algorithms［M］.Computer Vision and Applications.A Guide for Students and Practitioners，Chapter 11(B.Ja¨hne and H.Haussecker，Eds).San Diego，USA:Academic Press，2000:397－438.

【8】Digital Measuring Microscopes，MM400/800 Series.Nikon Corporation，Japan.［EBOL］.http://www.nikon.com/products/instruments/lineup/measuring _ instruments/index.htm.

【9】Optical Imaging Profiler PLu2300.Sensofar－ Tech SL.，Terrassa，Spain.［EBOL］.http://www.sensofar.com/products/products_plu2300.html.

【10】劉建偉，蔣志強，劉元朋，等.輸電塔架承載變形的三維光學(xué)測量［J］.光學(xué)精密工程，2012，20(5):942－948.

【11】張陽，臧順來，郭翔，等.基于數(shù)字散斑應(yīng)變測量法的薄板各向異性力學(xué)性能研究［J］.材料工程，2012(4):6－11.

【12】牛小兵，林玉池，趙美，等.光柵投影三維測量的原理及關(guān)鍵技術(shù)分析［J］.儀器儀表學(xué)報，2001，22(z4):203－205.

【13】顧元亮，李新軍，呂曉東，等.標(biāo)志點拼接技術(shù)在大型薄壁件型面非接觸測量中的應(yīng)用［J］.塑性工程學(xué)報，2004，11(6):26－28.

【14】馬揚飚，鐘約先，李仁舉，等.注塑件產(chǎn)品的三維無接觸全場檢測［J］.塑性工程學(xué)報，2006，13(1):94－97.

【15】皮敏捷，徐科.基于多條激光線的鋼板表面缺陷三維檢測方法［J］.機電產(chǎn)品開發(fā)與創(chuàng)新，2007，20(2):123－125.

【16】嚴(yán)子深，段潤保.依據(jù)零件三維檢測數(shù)據(jù)的機械加工誤差診斷［J］.河北理工學(xué)院學(xué)報，2005，5(2):40－43.

【17】孔多，賀成柱，于忠武.基于攝影測量技術(shù)的板式換熱器靜態(tài)變形試驗研究［J］.機械研究與應(yīng)用，2012(1):37－40.

【18】張德海，梁晉，唐正宗，等.基于近景攝影測量和三維光學(xué)測量的大幅面測量新方法［J］.中國機械工程，2009，20(7):817－822.