劉小龍,孟 婥,張玉井,孫以澤
(東華大學 機械工程學院,上海 201620)
碳纖維編織是由編織錠子放線、纖維傾斜交織,最終匯聚在芯模上形成預制件的過程[1]。編織過程極易發生摩擦導致纖維異常:如錠子放線過程中碳纖維與編織錠子上的導紗輪處產生摩擦,形成纖維纏繞狀態;纖維長時間磨損導致斷纖;碳纖維單絲斷裂纏繞在纖維筒上形成纖維環狀態。若無法及時發現編織過程中纖維異常狀態,將嚴重干擾正常的編織生產。Maidl等[2]的研究表明異常狀態的出現導致編織生產率降低并增加了錯誤原因分析和錯誤糾正時間的額外成本。Ebel等[3]的研究表明編織缺陷的出現,如纖維斷裂和纖維環,可導致機器停機時間高達總生產時間的26%,極大降低了編織機的生產效率。因此,生產時需及時去除異常纖維以保證編織效率。
綜上,大都研究都不適用于編織工藝,而視覺檢測提升了生產效率,相對人工有更高的檢測精度。因此,本文針對斷纖、纖維纏繞及纖維環狀態分別用NCC匹配算法、基于圖像分割的二次定位和幾何特征檢測算法及積分投影算法進行檢測,以期利用視覺檢測技術實現編織過程的纖維檢測,從而提高檢測效率及編織質量。
纖維狀態檢測系統的任務是檢測編織過程中有無纖維異常狀態,該系統主要包括相機、鏡頭、光源、傳感器等,系統構成如圖1所示。其中,相機芯片為CMOS類型且分辨率為1 440像素×1 080像素,光源為2個環形白色LED光源。相機1采集纖維環狀態和導紗輪1的纖維纏繞圖像,相機2采集斷纖狀態和導紗輪2的纖維纏繞圖像。
圖1 系統構成Fig.1 System Components
檢測系統算法流程分為預處理階段和狀態識別階段。系統整體檢測流程為:相機采集到當前位置圖像,經過預處理階段和狀態識別階段檢測纖維狀態,當檢測到纖維異常狀態時進行處理。預處理階段包括劃分感興趣區域和直方圖均衡化、中值濾波對圖像進行處理,預處理完成后針對產生斷纖、纖維纏繞及纖維環的位置分別采用NCC匹配算法、基于圖像分割的二次定位和幾何特征檢測算法及積分投影算法進行檢測。算法流程如圖2所示。
圖2 檢測系統算法流程Fig.2 Detection system algorithm flow
斷纖是編織過程中常見的一種異常現象,正常編織時由于張力的存在,導紗輪2在錠子的上部,當斷纖時張力為0,導紗輪2下墜到錠子底部,從而確定此纖維已經斷裂。纖維狀態特征如圖3所示。因此可通過導紗輪位置的變化作為特征間接表征斷纖,導紗輪2區域特征為白色,斷纖后的區域灰度變化明顯,采用NCC匹配算法可檢測斷纖狀態。
圖3 纖維狀態特征Fig.3 Fiber state characteristics.
2.1.1 基于NCC匹配算法的斷纖檢測
圖像匹配算法可以對不同時間采集到的多幅圖像進行不同程度的幾何變換,以達到各圖像中相同的區域在幾何上對準,該方法分為基于區域和基于特征兩大類[12]。本文選用的NCC匹配算法是基于區域的圖像配準[13],即利用像素點的灰度信息來確定待配準圖像之間的幾何變換關系。NCC匹配算法需要模板圖像和待配準圖像,若模板圖像T大小為m×m,待配準圖像S大小為n×n,則匹配過程如下:
①模板圖在待配準圖像上進行連續平移操作,每一次平移所覆蓋的子圖區域記為P(i,j),其中(i,j)為區域的左上角坐標位置;
②根據歸一化相關函數計算區域P與輸入的待配準圖像間的相關系數值;
③對待配準圖像進行自上而下、自左而右的遍歷搜索。歸一化相關系數N(x,y)的計算通常分為去均值(見式(1))和不去均值(見式(2))2種,通過式(1)或(2)可計算區域P與輸入的待配準圖像間的歸一化相關系數值,相關系數最大位置即為最佳匹配位置。本文采用式(1)進行計算。
(1)
(2)
式(1)和式(2)中:S是搜索圖像;T是模板圖像;z1和z2分別代表圖像匹配窗口內的像素平均值;歸一化相關系數N(x,y)的取值范圍為[-1,1],其值越大代表2幅圖像的相關性越大;N(x,y)=-1時,圖像之間完全不同;N(x,y)=1時,圖像之間完全相同。
2.1.2 檢測過程分析
模板圖I、待配準圖像J、歸一化相關系數最大時的匹配位置坐標(xmax,ymax)及歸一化相關系數值N如圖4所示。
圖4 匹配數據Fig.4 Match data
由圖4可知,根據NCC算法計算得到纖維正常與斷纖2種狀態下錠子位置的歸一化相關系數N,當歸一化相關系數N較大時,該位置的纖維為正常狀態;當歸一化相關系數N較小時,該位置的纖維為斷纖狀態。通過多次分析正常纖維狀態圖像的歸一化相關系數,得到該值波動結果如圖5所示。
圖5 波動范圍Fig.5 Fluctuation range
由圖5可知,統計分析得到的歸一化相關系數取值在區間[0.8,1]的范圍內波動,實驗設定當歸一化系數N>0.8時,為纖維正常狀態。
正常纖維在膠料的作用下會形成一致性較強的束,如圖6(a)所示,但隨著編織過程進行,纖維與錠子上的導紗輪不斷摩擦且在放線過程中纖維張力交替變化,使得纖維的集束性變差導致纖維產生如圖6(b)所示的纖維分離現象,當纖維產生一定程度的分離時會導致如圖6(c)所示的單絲斷裂現象,斷裂的單絲會纏繞在導紗輪上,纖維纏繞會加速正常纖維的磨損,從而加速纖維斷裂、影響編織成型效果。
圖6 纖維纏繞的形成過程Fig.6 The formation process of filament winding.(a)Normal fiber ;(b)Fiber separation;(c)Monofilament separation and fracture
針對導紗輪上纖維纏繞的檢測,采用粗、精定位共計二次定位確定導紗輪區域。二次定位過程如圖7所示,纖維纏繞檢測過程如圖8所示,具體為:
①在采集圖像上利用感興趣區域對錠子上存在導紗輪的區域進行粗定位,得到如圖8纖維纏繞檢測過程所示的粗定位區域;
②得到步驟①中的區域后,利用閾值分割、形態學運算和特征提取的方式進行精定位,提取得到如圖8所示精定位的導紗輪區域;
③填充干擾區域。
圖7 定位過程Fig.7 Positioning process
圖8 纖維纏繞檢測過程Fig.8 Filament winding detection process
通過上述步驟得到待識別的區域,后對待識別區域進行二值化操作,計算得到該區域中黑色像素區域面積在整體區域的面積占比如表1所示。
表1 黑色區域面積占比Tab.1 The proportion of black area
由表1可知,可以根據黑色區域面積占比判別纖維正常區域和纖維纏繞區域。對導紗輪的纖維纏繞設定面積占比高于15%時,當前區域有纖維纏繞,需要進行處理。
編織過程中纖維環的形成會使纖維張力反復變化,從而造成纖維進一步磨損,導致纖維偏離正常位置,從而造成編織缺陷。在采集的圖像上得到感興趣區域如圖9所示,可見纖維環的存在使該區域的灰度變化明顯。
圖9 感興趣區域Fig.9 Area of interest
積分投影法得到纖維區域后進行投影,對感興趣區域中的像素按行、列分別累加求和再求均值,即
(3)
(4)
式(3)和式(4)中:(x,y)為位置坐標;Q(x,y)代表位置像素;r為行數;c為列數;H(x)和V(y)分別代表第x行的平均像素和第y列的平均像素。
圖10 積分投影曲線Fig.10 Integral projection curve.(a)Horizontal projection;(b)Vertical projection
由式(3)(4)得如圖10所示的水平投影曲線和垂直投影曲線。由圖10可知,纖維正常和纖維環狀態在水平投影方向灰度變化區別較小,較難確定是否含有纖維環,在垂直投影方向上260~280列的灰度變化明顯,該波谷位置即為纖維環區域,由于無需確定纖維環具體位置,所以可根據垂直投影方向上的差值判別纖維正常區域和纖維纏繞區域,通過分析不同纖維環圖像垂直方向上投影的差值,得到如圖11所示的差值曲線。
圖11 差值曲線Fig.11 Difference curve
由圖11可知,A、B為高度差極值點,纖維環狀態在垂直方向上投影的差值在60以上,因此設定高度差大于60則存在纖維環。
圖像采集平臺如圖12所示。圖12(a)中為環形編織機,白色框為固定相機的支架,圖12(b)為左圖框的放大圖。
圖12 圖像采集平臺Fig.12 Platform of image acquisition(a)Field installation;(b)Local location
為驗證算法的可行性,對圖像采集平臺中采集到的圖片利用前述算法進行纖維狀態檢測,同時在不同光照(圖像最大閾值分別為210、230和255)、運動模糊、密度和為0.5和0.7的椒鹽噪聲圖像的條件下進行對比實驗,實驗所用的計算機為Win10系統的ipc-610整機(研華科技),編織機為環形編織機,其中實驗中的識別時間不包含增加干擾所用時間。編織機運行參數根據控制系統界面自主設定,界面參考如圖13所示。
圖13 控制系統界面Fig.13 Control system interface
斷纖的檢測流程圖如圖14所示。
圖14 檢測流程Fig.14 Detection process
不同環境下的識別準確率及識別時間統計如表2所示。斷纖檢測結果如圖15所示。
圖15 不同環境下的檢測結果Fig.15 Test results in different environments
表2 對比檢測結果Tab.2 Compare test results
由表2可知,在實驗的不同環境中得到的圖像結果均顯示,所用算法對于斷纖狀態的檢測識別準確率均在94%以上,識別時間在0.5 s左右,識別精度較好。
纖維纏繞檢測流程圖如圖16所示。
圖16 檢測流程Fig.16 Detection process
纖維纏繞檢測實驗結果如圖17所示。不同環境下的識別準確率及識別時間如表3所示。根據表3可以發現算法識別準確率在94%以上,識別時間在0.4 s~0.5 s,檢測結果和檢測時間能滿足基本的應用要求。
圖17 不同環境下的檢測結果Fig.17 Test results in different environments
表3 對比檢測結果Tab.3 Compare test results
纖維環狀態檢測流程圖如圖18所示。
圖18 檢測流程Fig.18 Detection process
纖維環檢測結果如圖19所示。不同環境下的識別準確率及識別時間如表4所示。
圖19 不同環境下的檢測結果Fig.19 Test results in different environments
表4 對比檢測結果Tab.4 Compare test results
由表4可知,所用算法的平均識別準確率為93.50%,精度較高,識別時間較短。
根據前述對不同纖維狀態檢測的實驗數據可知,纖維狀態檢測算法在不同環境下雖有所差異,但總體識別率較高。檢測時間較短。斷纖檢測所用算法的平均識別準確率為95.75%,平均識別時間為0.54 s;纖維纏繞狀態所用算法的平均識別準確率為95.25%,平均識別時間為0.47 s;纖維環狀態所用算法的平均識別準確率為93.50%,平均識別時間為0.48 s。由此可知,算法可以克服一定的外在環境干擾,相對于人工檢測,提高了對纖維狀態的檢測效率。同時結合算法,設計了如圖20所示的纖維狀態檢測系統界面。
圖20 檢測系統軟件界面Fig.20 Detection system software interface
通過利用機器視覺技術,設計了一種纖維狀態檢測系統,分別對3種編織過程纖維狀態進行了實驗研究,得出以下結論:
①針對編織過程中出現的斷纖狀態,將正常纖維狀態作為匹配模板并將采集圖像作為待匹配對象,利用NCC匹配算法得到不同狀態間的相關系數以完成檢測。實驗表明,在不同環境下該算法的平均識別準確率為95.75%。
②針對編織過程中出現的纖維纏繞狀態,采用基于圖像分割的二次定位提取區域并計算區域二值化后的黑色區域面積占比進行檢測。實驗表明,在不同環境下該算法的平均識別準確率為95.25%,檢測時間在0.5 s內。
③針對編織過程中出現的纖維環狀態,經過濾波后進行積分投影實現了對該狀態的檢測。實驗表明,在不同環境下該算法的平均識別準確率為93.50%。
基于機器視覺進行纖維狀態檢測為編織過程纖維狀態的自動化檢測提供了參考依據和方法。