潘玉晶，陳柏金，蘆光榮

(1.華中科技大學(xué) 材料成形與模具技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室，湖北 武漢 430074；2.蘭州蘭石能源裝備工程研究院有限公司，甘肅 蘭州 730050)

棒料彎曲形狀自動(dòng)測(cè)量研究

潘玉晶1，陳柏金1，蘆光榮2

(1.華中科技大學(xué) 材料成形與模具技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室，湖北 武漢 430074；2.蘭州蘭石能源裝備工程研究院有限公司，甘肅 蘭州 730050)

本文基于非接觸測(cè)量方式，提出采用二維線輪廓傳感器，全程自動(dòng)掃描棒料輪廓，利用最小二乘法對(duì)截面中心進(jìn)行擬合，由連續(xù)的截面中心確定棒料中心線，實(shí)現(xiàn)了棒料彎曲量計(jì)算以及彎曲形狀的確定算法，為棒料自動(dòng)校直工藝參數(shù)的計(jì)算奠定了基礎(chǔ)。實(shí)驗(yàn)研究證明，這種測(cè)量方案不僅掃描速度快，穩(wěn)定性好，而且對(duì)棒料的彎曲形狀測(cè)量準(zhǔn)確，滿足自動(dòng)校直工藝要求。

輪廓傳感器；棒料測(cè)量；最小二乘擬合；彎曲參數(shù)

在各種壓力機(jī)上生產(chǎn)的棒料常因加熱溫度、變形量、鍛造尺寸誤差等發(fā)生一定的彎曲或扭曲變形。據(jù)不完全統(tǒng)計(jì)，我國(guó)制造業(yè)每年生產(chǎn)數(shù)萬(wàn)根棒料零件，其中約有70％以上存在不同程度的彎曲或扭曲[1]，若這些棒料未得到有效校直將導(dǎo)致在后續(xù)加工過程中產(chǎn)生不均勻的切削量，造成較高的廢品率。棒料彎曲形狀總體可分為平面單弧度彎曲、多弧度彎曲和空間復(fù)雜彎曲[2]。對(duì)不同的彎曲形狀應(yīng)制定不同的校直工藝方案，故計(jì)算棒料彎曲參數(shù)、確定彎曲形狀對(duì)后續(xù)校直操作具有十分重要的意義。

目前，針對(duì)棒料零件的校直大多采用三點(diǎn)反彎方式，校直前對(duì)棒料彎曲形狀的確定，通常依據(jù)人工卡尺測(cè)量或接觸式傳感器測(cè)量[3]。前者依據(jù)工人經(jīng)驗(yàn)，憑肉眼判斷彎曲形態(tài)，確定校直位置以及校直下壓量，并進(jìn)行反復(fù)校正。不僅勞動(dòng)強(qiáng)度大，工作效率低，而且校直精度低，穩(wěn)定性差。接觸式傳感器測(cè)量，通過旋轉(zhuǎn)工件測(cè)量其輪廓尺寸，不僅旋轉(zhuǎn)裝置設(shè)計(jì)復(fù)雜，且在工件旋轉(zhuǎn)與傳感器接觸測(cè)量過程中，極易造成傳感器的損壞。本文結(jié)合企業(yè)生產(chǎn)實(shí)際情況，基于非接觸式測(cè)量，采用二維激光位移傳感器與校直機(jī)相結(jié)合，掃描并提取棒料輪廓特征數(shù)據(jù)，擬合截面輪廓中心，精確計(jì)算棒料彎曲參數(shù)，實(shí)現(xiàn)了棒料彎曲形狀的自動(dòng)檢測(cè)，并為全自動(dòng)校直奠定了堅(jiān)實(shí)的基礎(chǔ)。

1 棒料輪廓掃描

二維激光線輪廓傳感器是基于線結(jié)構(gòu)光視覺檢測(cè)方式的一種測(cè)量方法，掃描速度快，測(cè)量精度高，廣泛應(yīng)用于在線和現(xiàn)場(chǎng)實(shí)時(shí)測(cè)量等過程。本文采用Gocator2330傳感器，通過向棒料輪廓發(fā)射一束可見激光，在棒料表面形成一條激光線。傳感器內(nèi)部面陣掃描相機(jī)從另一角度感知棒料表面的反射激光，通過檢測(cè)反射角度從而確定目標(biāo)表面輪廓位置，經(jīng)內(nèi)部數(shù)據(jù)處理單元獲得棒料表面輪廓點(diǎn)集的位置信息。Gocator2330傳感器掃描的每條剖面輪廓包含的數(shù)據(jù)點(diǎn)數(shù)為1280個(gè)，豎直方向分辨率為0.006mm～0.014mm，水平方向分辨率為0.044mm～0.075mm，傳感器尺寸大小49mm×75mm×142mm。

傳感器安裝在校直機(jī)主機(jī)內(nèi)，校直工作臺(tái)與校直機(jī)主機(jī)采用一體化設(shè)計(jì)，如圖1所示。測(cè)量過程中棒料由左右托輥架裝夾固定，傳感器隨校直機(jī)一起勻速移動(dòng)對(duì)棒料表面進(jìn)行掃描測(cè)量，獲取棒料輪廓特征數(shù)據(jù)，并發(fā)送到控制程序分析處理，計(jì)算棒料相關(guān)彎曲參數(shù)。掃描示意圖如圖2所示。

2 測(cè)量數(shù)據(jù)處理

最小二乘法，通過最小化誤差的平方和尋找數(shù)據(jù)的最佳函數(shù)匹配。它用最簡(jiǎn)的方法求得一些不可知的真值，而令誤差平方之和最小[4]。傳感器掃描結(jié)束后，基于最小二乘法對(duì)掃描數(shù)據(jù)進(jìn)行分析處理，擬合各個(gè)截面輪廓，獲取截面中心，確定棒料軸線。

圖1 傳感器安裝固定

圖2 掃描示意圖

2.1 截面中心擬合

鍛造后棒料毛坯截面近似為圓，截面中心的確定即為圓心擬合的過程?；谟邢迾颖军c(diǎn)集進(jìn)行圓擬合的方法有很多，例如最小二乘法、加權(quán)平均法、Hough變換法等[5]。受棒料輪廓形狀、傳感器結(jié)構(gòu)以及現(xiàn)場(chǎng)測(cè)量環(huán)境等因素的影響，每次掃描所獲取的點(diǎn)十分密集但占整個(gè)圓周的比例不足1/2，且可能存在周圍物體的干擾，所以可選用最小二乘法擬合，不僅計(jì)算精度高，而且快速高效，符合在線檢測(cè)的實(shí)時(shí)性要求。

如圖2所示，針對(duì)某個(gè)截面樣本點(diǎn)集(Xi，Yi)(i=1……n)：

(1)確定初始參考圓心、半徑。截面點(diǎn)集中距離傳感器最近點(diǎn)M(Xj，Ymin)必為棒料輪廓點(diǎn)，在M點(diǎn)前后取一定數(shù)量的測(cè)量點(diǎn)，確保P～Q之間的點(diǎn)均為棒料輪廓點(diǎn)。然后利用最小二乘圓擬合算法[6]，確定初始圓心C0(X0，Y0)和初始半徑R0；

(2)剔除干擾點(diǎn)。遍歷截面所有樣本點(diǎn)，如果樣本點(diǎn)到初始圓心C0的距離d與初始半徑R0的差值大于5mm(閾值根據(jù)實(shí)際情況確定)，則將其剔出截面點(diǎn)集；

(3)擬合截面圓心和半徑。對(duì)剔除干擾點(diǎn)后的樣本點(diǎn)集進(jìn)行最小二乘擬合，得到截面圓心(Xr，Yr)，半徑R。

2.2 棒料軸線擬合

棒料軸線包括其實(shí)際彎曲軸線和棒料中心線。校直前由于棒料彎曲變形較為復(fù)雜，實(shí)際軸線為一條空間曲線。理想狀態(tài)下，棒料校直后軸線為標(biāo)準(zhǔn)的空間直線，與其中心線重合。對(duì)被測(cè)棒料的n個(gè)截面分別進(jìn)行截面中心擬合，獲得n個(gè)圓心坐標(biāo)。由于各截面掃描間距很小，依次連接各截面的擬合圓心所形成的空間折線，可近似看作是被測(cè)棒料的實(shí)際軸線。在實(shí)際生產(chǎn)操作中，對(duì)于端部存在翹曲變形的棒料，直接將其作為短材出廠或廢品處理。因此，可直接用棒料端部截面中心進(jìn)行中心線的擬合。對(duì)一系列離散數(shù)據(jù)點(diǎn)進(jìn)行直線擬合，常用的方法有最小二乘法、契比雪夫法等，其中最小二乘法最為常見和重要，根據(jù)已知數(shù)據(jù)坐標(biāo)計(jì)算到直線的偏差距離，滿足距離之和最小的即為所求直線方程[7]。空間直線的標(biāo)準(zhǔn)方程[8]可寫為：

整理可得

此時(shí)將空間直線方程轉(zhuǎn)換成了兩個(gè)平面方程，對(duì)兩個(gè)平面方程分別采用最小二乘數(shù)據(jù)擬合即可獲得最佳空間直線方程，即為棒料中心線。具體求解過程可參考文獻(xiàn)[7]。

3 棒料彎曲參數(shù)計(jì)算

3.1 截面中心偏移量計(jì)算

如圖3所示，棒料各截面中心與基準(zhǔn)中心存在一定的偏移，截面的偏移量為一個(gè)矢量，具有大小和方向。圖4中，假設(shè)點(diǎn)P為某截面的擬合圓心，AB為棒料中心線，過 P(Xp，Yp，Zp)作 AB 的垂線，垂足為 Q(Xq，Yq，Zq)。PQ 長(zhǎng)度即為該截面的偏移量：

圖3 棒料軸線示意圖

圖4 偏移角度計(jì)算示意圖

截面偏移方向用PQ與豎直面F的夾角?表示，下文簡(jiǎn)稱偏移角度。校直時(shí)按照偏移角度將彎曲截面中心旋轉(zhuǎn)至y軸正上方，根據(jù)旋轉(zhuǎn)方向，逆時(shí)針記為正，順時(shí)針記為負(fù)，旋轉(zhuǎn)角度范圍為-180°～180°。具體求解算法流程如下：

(1)求解棒料中心線AB的方向向量：任取AB上兩點(diǎn)可得其方向向量

(4)確定旋轉(zhuǎn)方向：由于向量之間的夾角為0～180°，故需進(jìn)一步計(jì)算與 Z 軸的夾角 β，從而判定 ?的正負(fù)；若 β∈(0，90°)，則 ?為正，逆時(shí)針旋轉(zhuǎn)，若 β∈(90°，180°)，，則 ?為負(fù)，順時(shí)針旋轉(zhuǎn)。

對(duì)所有截面重復(fù)上述步驟，即可得Di和?i，其中i=1…n，Di為第 i個(gè)截面的偏移量，?i為第 i個(gè)截面的偏移角度大小。

3.2 棒料彎曲形狀確定

棒料本身是一個(gè)剛體，其軸線在空間上是一條連續(xù)的曲線。當(dāng)采樣截面間距很小時(shí)，可直接遍歷各截面偏移量Di(i=1…i…j…n)，劃分棒料彎曲段，再確定各彎曲段內(nèi)最大彎曲量的截面位置。彎曲段劃分算法流程如圖5所示：若截面彎曲量從第i個(gè)截面開始呈遞增趨勢(shì)，然后在某截面開始遞減，直到第j個(gè)截面接近于零，且i-j截面之間的最大偏移量在需校直范圍內(nèi)，則i-j即為一個(gè)校直彎曲段，否則認(rèn)為是正常波動(dòng)，無(wú)需處理。由于采樣截面不一定是最大偏移量的截面，因此需對(duì)各彎曲段內(nèi)所有采樣截面的中心進(jìn)行曲線擬合，求其曲率最大的截面中心，即為該彎曲段最大彎曲量截面中心位置。

圖5 彎曲段劃分流程圖

劃分彎曲段后若某彎曲段內(nèi)所有截面偏移角度近似相等或在較小范圍內(nèi)波動(dòng)，說明該彎曲段中所有截面均朝著同一個(gè)方向彎曲，則可認(rèn)為該彎曲段為平面彎曲。若彎曲段內(nèi)各截面偏移角度在較大范圍內(nèi)波動(dòng)，則判定該段為空間扭曲。如果棒料只存在一個(gè)彎曲段，即為單弧度彎曲，存在多個(gè)彎曲段，即為多弧度彎曲。后續(xù)校直工序中則可將多弧度彎曲拆分成多個(gè)單弧度彎曲，空間扭曲則根據(jù)扭曲程度確定修正系數(shù)將其近似為平面彎曲。

4 實(shí)驗(yàn)及結(jié)論

如圖6所示，采用Gocator2330線輪廓傳感器對(duì)多批次棒料進(jìn)行掃描測(cè)量，擬合計(jì)算其直徑，最終誤差均在±1mm范圍內(nèi)。此外，在UG中建立如圖7所示的彎曲棒料模型，左端面中心為坐標(biāo)原點(diǎn)，MN段沿豎直面向下彎曲，偏移角度為-180°，PQ段為空間扭曲。然后在棒料軸線上等間距取500個(gè)點(diǎn)，記為截面中心點(diǎn)集，導(dǎo)出所有點(diǎn)三維坐標(biāo)。然后按照上述方法進(jìn)行擬合計(jì)算，得到彎曲段MN和PQ，其中MN段截面中心偏移角度均在-180°±2°范圍內(nèi)變化，PQ段截面中心偏移角度在較大范圍內(nèi)波動(dòng)，為空間彎曲。計(jì)算結(jié)果與UG模型參數(shù)相符。

本文采用二維線輪廓傳感器與校直機(jī)結(jié)合，基于最小二乘法擬合棒料截面中心，確定棒料軸線，計(jì)算截面偏移量和偏移角度，確定彎曲位置，依據(jù)彎曲段內(nèi)各截面偏移角度波動(dòng)范圍判斷彎曲形狀，實(shí)現(xiàn)了棒料自動(dòng)掃描測(cè)量及棒料彎曲參數(shù)的精確計(jì)算。設(shè)備自動(dòng)化程度高，檢測(cè)范圍大，運(yùn)動(dòng)控制精準(zhǔn)，測(cè)量精度高。實(shí)驗(yàn)表明，此方案計(jì)算結(jié)果與實(shí)際棒料彎曲情況一致，對(duì)后續(xù)校直工藝方案的制定具有重要的指導(dǎo)意義。

圖6 現(xiàn)場(chǎng)測(cè)量示意圖

圖7 彎曲棒料UG模型

[1] 翟 華.軸類零件校直工藝?yán)碚撗芯縖D].合肥：合肥工業(yè)大學(xué)，2003.

[2] 韓 江，陸青松，夏 鏈，等.軸類零件精密校直工藝?yán)碚摷皯?yīng)用研究[J].組合機(jī)床與自動(dòng)化加工技術(shù)，2007，(8)：74-77.

[3] 于 鵬.新型大鍛件尺寸測(cè)量系統(tǒng)的關(guān)鍵技術(shù)研究[D].上海：上海交通大學(xué)，2008.

[4] 賈小勇，徐傳勝，白 欣.最小二乘法的創(chuàng)立及其思想方法[J].西北大學(xué)學(xué)報(bào)(自然科學(xué)版)，2006，(3)：507-511.

[5] 李 淵，喬 兵，陸宇平.基于聚類方法的圓擬合算法[J].計(jì)算機(jī)應(yīng)用，2013，S2：206-208.

[6] 張應(yīng)芳，張愛國(guó)，蘇積聰，等.觀測(cè)圓上多個(gè)任意點(diǎn)及其圓心坐標(biāo)擬合[J].黑龍江工程學(xué)院學(xué)報(bào)，2014，(2)：27-31.

[7] 襲 楊.空間直線擬合的一種方法[J].齊齊哈爾：齊齊哈爾大學(xué)學(xué)報(bào)，2009，(2)：64-68.

[8] 呂林根，許子道.解析幾何[M].第四版.北京：高等教育出版社，2006.

Study on automatic measurement of bar bending shape

PAN Yujing,CHEN Baijin,LU Guangrong
(State Key Laboratory of Materials Processing and Die&Mould Technology,Huazhong University of Science and Technology,Wuhan 430074,Hubei China)

On the basis of non-contact measurement mode,the two-dimensional contour sensor has been adopted to automatically scan the bar profile.The section center has been fitted by use of least-square method.The bar centerline has been confirmed by the continuous section center.The calculation of the bar bending volume and the algorithm of bending profile have been achieved,which could lay foundation for calculation of automatic straighteningprocess parameters.The test study shows that this measure method has high scanning speed and good stability which can accurately measure the bar bending profile.It can satisfy the automatic straightening process demands.

Contour sensor;Bar measurement;Least-square fitting;Curvature parameter

TH273

A

10.16316/j.issn.1672-0121.2017.05.016

1672-0121(2017)05-0054-03

2017-03-12；

2017-04-25

潘玉晶(1991-)，女，碩士在讀，主攻自動(dòng)校直工藝研究。E-mail：1508162553@qq.com；陳柏金(1965-)，男，博士，教授，從事機(jī)電液一體化研究