劉怡然，李麗君，杜月云

（1.商丘職業(yè)技術(shù)學(xué)院交通學(xué)院，河南 商丘 476000；（2.山東理工大學(xué)交通與車(chē)輛工程學(xué)院，山東 淄博 255049）

1 引言

隨著時(shí)代的發(fā)展，汽車(chē)已經(jīng)成為必不可少的家庭用品，但是隨著使用時(shí)間增加，汽車(chē)齒輪會(huì)出現(xiàn)不同程度的磨損，尤其是后橋螺旋錐齒輪的磨損。這是由于汽車(chē)在行駛過(guò)程中，后橋存在耦合豎向振動(dòng)，從而導(dǎo)致后橋螺旋錐齒輪出現(xiàn)磨損現(xiàn)象，如果不進(jìn)行及時(shí)的檢測(cè)，則會(huì)降低汽車(chē)行駛的安全性。研究有效的齒輪磨損檢測(cè)方法，提高齒輪磨損的檢測(cè)精度，提高汽車(chē)行駛的安全性，為人們的交通出行提供安全保障［1-2］。

文獻(xiàn)［3］提出基于小波包能量和調(diào)制信號(hào)雙譜邊帶估計(jì)的齒輪磨損檢測(cè)方法，采用小波包變換將齒輪磨損信號(hào)分解為多個(gè)頻段，然后通過(guò)調(diào)制信號(hào)雙譜邊帶估計(jì)方法從磨損信號(hào)中提取齒面磨損特征，最后通過(guò)檢測(cè)指標(biāo)完成齒輪磨損檢測(cè)。但是，該方法未考慮噪聲的干擾問(wèn)題，導(dǎo)致齒輪磨損檢測(cè)精度不足。文獻(xiàn)［4］提出基于逆向工程的齒輪磨損檢測(cè)方法，該方法以齒輪磨損檢測(cè)原理為基礎(chǔ)，構(gòu)建齒輪磨損區(qū)域的轉(zhuǎn)換矩陣，并且對(duì)磨損區(qū)域進(jìn)行誤差補(bǔ)償。結(jié)合補(bǔ)償點(diǎn)構(gòu)建NURBS曲面擬合方程，完成齒輪磨損的檢測(cè)。但是，該方法的磨損曲面擬合結(jié)果會(huì)出現(xiàn)未磨損區(qū)域，導(dǎo)致檢測(cè)效果不足。文獻(xiàn)［5］提出基于改進(jìn)Mask Scoring RCNN的齒輪磨損檢測(cè)方法，該方法將殘差網(wǎng)絡(luò)與特征金字塔網(wǎng)絡(luò)作為主干網(wǎng)絡(luò)，通過(guò)語(yǔ)義信息與細(xì)節(jié)特征融合進(jìn)行齒輪磨損特征的提取。對(duì)局部特征層進(jìn)行歸一化處理后，計(jì)算二值化掩碼圖形信息，完成齒輪磨損的檢測(cè)。但是該方法未考慮數(shù)據(jù)的去噪問(wèn)題，導(dǎo)致樣本數(shù)據(jù)中的噪聲過(guò)大，降低了齒輪磨損的檢測(cè)精度。

為了解決上述齒輪磨損檢測(cè)方法存在的問(wèn)題，提出Delaunay三角剖分的汽車(chē)螺旋錐齒輪磨損檢測(cè)方法。采用Delaunay三角剖分對(duì)齒面的離散數(shù)據(jù)相鄰掃描線進(jìn)行處理，完成齒面非特征離散數(shù)據(jù)分塊。采用Hermite 插值算法對(duì)劃分的Delaunay 三角網(wǎng)格進(jìn)行遍歷檢測(cè)，完成齒輪磨損的檢測(cè)。

2 齒面非特征離散數(shù)據(jù)分塊

利用NURBS（Non-Uniform Rational B-Splines）曲面構(gòu)建汽車(chē)螺旋錐齒輪模型［6］，通過(guò)Delaunay三角剖分進(jìn)行齒輪表面數(shù)據(jù)的非特征分塊處理。

Delaunay三角形網(wǎng)中任一三角形的外接圓范圍內(nèi)不會(huì)有其它點(diǎn)存在，并且以最近的點(diǎn)形成三角形，不論從齒面何處開(kāi)始處理，都會(huì)得到一致性的結(jié)果。將Delaunay三角剖分應(yīng)用于齒面非特征數(shù)據(jù)的分塊，可以精準(zhǔn)描述汽車(chē)螺旋錐齒輪表面的區(qū)域，有助于提高齒輪磨損的檢測(cè)精度。

利用張量積形式完成齒輪的NURBS曲面的構(gòu)建，表達(dá)式具體如下：

式中：i、b—曲面方向；A(i，b)—NURBS曲面函數(shù)；q—方向i的中心點(diǎn)個(gè)數(shù)；m—方向b的中心點(diǎn)個(gè)數(shù)；Bok—曲面中心點(diǎn)；No，l—方向i的樣條基函數(shù)；Nk，l—方向b的樣條基函數(shù)；ξok—權(quán)因子。

在通過(guò)掃描式方法測(cè)量出齒輪表面的離散數(shù)據(jù)后，利用Delaunay三角剖分原理處理表面的離散數(shù)據(jù)。設(shè)αi表示掃描線，則有：i=1，2，3，…，n，αl，αl+1表示兩條鄰近掃描線，αl(k)，k=1，2，3，…，y；αl+1(K)，K=1，2，3，…，y表示測(cè)試點(diǎn)，y表示測(cè)試點(diǎn)個(gè)數(shù)。

連接兩條相鄰掃描線上的起始點(diǎn)和終止點(diǎn)，再連接一條掃描線上各測(cè)量點(diǎn)對(duì)應(yīng)在另一條線上的最近測(cè)量點(diǎn)，將兩掃描線之間的空間區(qū)域被劃分為若干個(gè)區(qū)域，如圖1（a）所示。檢測(cè)出三角形及四邊形空間區(qū)域，并依據(jù)最小內(nèi)角最大準(zhǔn)則剖分四邊形空間區(qū)域，最終完成Delaunay三角剖分，如圖1（b）所示。

圖1 相鄰掃描線的Delaunay三角剖分Fig.1 Delaunay Triangulation of Adjacent Scan Lines

通過(guò)Delaunay三角剖分對(duì)齒面非特征離散數(shù)據(jù)進(jìn)行分塊處理，選定NURBS齒面模型上任意一三角形為中心三角形CenTri，計(jì)算全部中心三角形值，其公式如下：

式中：n=1，2，3，εn—兩個(gè)三角形之間的夾角；Y—中心三角形區(qū)域值。

選取區(qū)域值最小的中心三角形作為起始三角形，如圖2（a）所示。選擇與起始三角形夾角最小的相鄰三角形作為中心三角形，如滿足式（2），則排除起始中心三角，計(jì)算周邊三角形與起始三角形的夾角度數(shù)，取度數(shù)最小的三角形為中心三角形，否則，繼續(xù)重新選取起始三角形。依照上述操作，遍歷全部三角形完成齒輪非特征離散數(shù)據(jù)分塊，如圖2（b）所示。

圖2 齒輪模型的非特征離散數(shù)據(jù)分塊Fig.2 Non-Feature Discrete Data Blocking for Gear Models

通過(guò)三角區(qū)域插值算法計(jì)算出齒面數(shù)據(jù)，實(shí)現(xiàn)齒輪表面區(qū)域的全面描述。

3 Hermite插值算法的齒輪磨損檢測(cè)

為了準(zhǔn)確的檢測(cè)出齒輪的磨損區(qū)域，以上述齒面非特征離散數(shù)據(jù)分塊結(jié)果為基礎(chǔ)，采用Hermite插值算法中的局部均值分解算法進(jìn)行齒面磨損的檢測(cè)。具體的檢測(cè)步驟如下：

（1）設(shè)定原始離散數(shù)據(jù)為c(y)，將原始強(qiáng)噪離散數(shù)據(jù)c(y)變?yōu)樾?shù)數(shù)據(jù)［7-9］。

（2）在級(jí)聯(lián)雙穩(wěn)隨機(jī)共振（CBSR，cascadedbistable stochastic resonance）系統(tǒng)中對(duì)步驟（1）中得到的小參數(shù)數(shù)據(jù)進(jìn)行去噪處理［10］，得到去噪后的離散數(shù)據(jù)cn(y)。

（3）對(duì)去噪處理后的離散數(shù)據(jù)極值點(diǎn)ml實(shí)行延展處理，以得到新的序列cl。

（4）從新序列的一端開(kāi)始，在其上選擇最大點(diǎn)或最小點(diǎn)，通過(guò)Hermite插值算法生成兩條包絡(luò)線。

（5）計(jì)算局部均值函數(shù)，其式如下：

式中：q(y)—局部均值函數(shù)；Rup(y)—上包絡(luò)線；Rdown(y)—下包絡(luò)線。

（6）計(jì)算包絡(luò)估計(jì)函數(shù)，其公式如下：

式中：s(y)—包絡(luò)估計(jì)函數(shù)。

（7）將采用LMD分析方法從原始齒面數(shù)據(jù)中分離出包絡(luò)數(shù)據(jù)，并且計(jì)算包絡(luò)數(shù)據(jù)構(gòu)建PF（product function）分量幅值，其公式如下：

式中：fo(y)—PF分量幅值，fo(y)的取值范圍為［0，1］，當(dāng)式（5）的計(jì)算結(jié)果不在取值范圍內(nèi)則判定該區(qū)域內(nèi)出現(xiàn)磨損。

按照上述步驟無(wú)限次處理至所有齒輪的三角網(wǎng)格，直到全部的三角網(wǎng)格完成遍歷檢測(cè)。

4 實(shí)驗(yàn)分析

為了驗(yàn)證Delaunay三角剖分的汽車(chē)螺旋錐齒輪磨損檢測(cè)方法的有效性，進(jìn)行對(duì)比測(cè)試實(shí)驗(yàn)。

此次研究用選擇桑塔納2015款手動(dòng)汽車(chē)的螺旋錐齒輪作為研究對(duì)象，選擇同款齒輪20 個(gè)進(jìn)行研究，該螺旋錐齒輪工作平穩(wěn)、傳動(dòng)比大、傳遞動(dòng)力大、結(jié)構(gòu)緊湊，可以滿足此次實(shí)驗(yàn)的需求。

汽車(chē)的螺旋錐齒輪的材料屬性，如表1所示。此次實(shí)驗(yàn)使用的設(shè)備，如圖3所示。汽車(chē)螺旋錐齒輪的磨損情況，如圖4所示。磨損齒輪的Delaunay三角剖分結(jié)果，如圖5所示。

表1 汽車(chē)螺旋錐齒輪材料屬性Tab.1 Material Properties of Automotive Spiral Bevel Gears

圖3 實(shí)驗(yàn)設(shè)備Fig.3 Experimental Equipment

圖4 螺旋錐齒輪磨損情況Fig.4 Wear of Spiral Bevel Gears

圖5 磨損齒輪的Delaunay三角剖分結(jié)果Fig.5 Delaunay Triangulation Results of Worn Gears

以上述實(shí)驗(yàn)設(shè)備與材料屬性數(shù)據(jù)為基礎(chǔ)，為了得到更加精準(zhǔn)的實(shí)驗(yàn)結(jié)果，對(duì)實(shí)驗(yàn)流程進(jìn)行規(guī)劃。具體的實(shí)驗(yàn)流程為：

（1）保持齒輪磨損實(shí)驗(yàn)設(shè)備的持續(xù)運(yùn)轉(zhuǎn)，運(yùn)轉(zhuǎn)時(shí)間為200h，采集200h內(nèi)齒輪的磨損數(shù)據(jù)，數(shù)據(jù)采集的時(shí)間間隔為5min，一小時(shí)內(nèi)可收集12組磨損數(shù)據(jù)。

（2）設(shè)定以齒輪磨損深度、齒輪磨損率、磨損區(qū)域面積檢測(cè)精度為此次實(shí)驗(yàn)的測(cè)試指標(biāo)。

（3）選擇文獻(xiàn)［4］提出的基于逆向工程的檢測(cè)方法、文獻(xiàn)［5］提出的基于改進(jìn)Mask Scoring R-CNN的檢測(cè)方法作為此次實(shí)驗(yàn)的對(duì)比方法。

（4）按照設(shè)定的實(shí)驗(yàn)指標(biāo)，采用三種不同方法進(jìn)行實(shí)驗(yàn)驗(yàn)證。

4.1 齒輪磨損深度

齒輪磨損主要為主動(dòng)齒輪與從動(dòng)齒輪之間的磨損，齒輪的磨損深度會(huì)隨著運(yùn)行時(shí)間的增加發(fā)生變化。主動(dòng)齒輪與從動(dòng)齒輪的磨損深度結(jié)果，如圖6所示。

圖6 齒輪磨損深度結(jié)果Fig.6 Results of Gear Wear Depth

圖6中，磨損深度0為初始表面狀態(tài)，并且齒輪磨損程度較高的位置為齒頂與齒根。將三種方法的磨損深度檢測(cè)結(jié)果與實(shí)際結(jié)果進(jìn)行對(duì)比，發(fā)現(xiàn)所研究方法的磨損深度結(jié)果與實(shí)際測(cè)量基本一致，而文獻(xiàn)［4］與文獻(xiàn)［5］方法的磨損深度值與實(shí)際測(cè)量結(jié)果相差較大。因此，說(shuō)明所研究方法能夠?qū)X輪的磨損深度進(jìn)行有效的分析。

4.2 齒輪磨損率

齒輪磨損率是指單位時(shí)間內(nèi)齒輪材料的磨損量，可以對(duì)齒輪的磨損情況進(jìn)行有效的體現(xiàn)，因此同樣將三種方法的磨損率檢測(cè)結(jié)果與實(shí)際磨損率測(cè)算結(jié)果進(jìn)行對(duì)比分析。磨損率結(jié)果，如圖7所示。

圖7 齒輪磨損率結(jié)果Fig.7 Results of Gear Wear Rate

從圖7所示的齒輪磨損率對(duì)比結(jié)果中可以看出，在對(duì)20個(gè)齒輪進(jìn)行研究的過(guò)程中，所研究方法的磨損率與實(shí)際磨損率始終保持一致，最大誤差不超過(guò)0.01μm，而文獻(xiàn)［4］方法與文獻(xiàn)［5］方法的磨損率則與實(shí)際磨損率相差較大。因此，說(shuō)明所研究方法能夠獲得更加精準(zhǔn)的齒輪磨損率結(jié)果，能夠提高齒輪磨損的檢測(cè)效果。

4.3 齒輪磨損區(qū)域面積檢測(cè)精度

為了進(jìn)一步驗(yàn)證所研究方法的齒輪磨損檢測(cè)性能，以齒輪磨損區(qū)域面積檢測(cè)精度為指標(biāo)，同樣將所研究方法與文獻(xiàn)［4］方法、文獻(xiàn)［5］方法進(jìn)行對(duì)比驗(yàn)證。每次實(shí)驗(yàn)對(duì)20個(gè)齒輪的磨損區(qū)域面積進(jìn)行檢測(cè)，共進(jìn)行10次實(shí)驗(yàn)，每次實(shí)驗(yàn)的磨損區(qū)域面積檢測(cè)均值結(jié)果，如表2所示。

表2 齒輪磨損區(qū)域面積檢測(cè)精度Tab.2 Inspection Accuracy of Gear Wear Area

從表2所示的齒輪磨損區(qū)域面積檢測(cè)精度對(duì)比結(jié)果中可以看出，對(duì)比文獻(xiàn)［4］與文獻(xiàn)［5］方法，所研究方法的磨損區(qū)域面積檢測(cè)精度最高達(dá)到了98.7%，而兩種對(duì)比方法的磨損區(qū)域面積檢測(cè)精度最高值分別為79.4%、84.3%。因此，說(shuō)明所研究方法能夠?qū)X輪磨損區(qū)域的面積進(jìn)行精準(zhǔn)地檢測(cè)。

5 結(jié)束語(yǔ)

汽車(chē)螺旋錐齒輪磨損的檢測(cè)有利于汽車(chē)駕駛的安全性，對(duì)汽車(chē)的安全穩(wěn)定運(yùn)行具有十分重要的現(xiàn)實(shí)意義，因此提出Delaunay三角剖分的汽車(chē)螺旋錐齒輪磨損檢測(cè)方法。對(duì)齒面相鄰掃描線進(jìn)行Delaunay三角剖分，完成齒面非特征離散數(shù)據(jù)的分塊。根據(jù)分塊結(jié)果，采用Hermite插值計(jì)算PF分量復(fù)制，遍歷全部的三角網(wǎng)格完成磨損區(qū)域的檢測(cè)。

實(shí)驗(yàn)數(shù)據(jù)表明，所提出方法的磨損深度與磨損率結(jié)果與實(shí)際值基本一致，并且能夠?qū)δp區(qū)域面積進(jìn)行精準(zhǔn)地檢測(cè)，最高檢測(cè)精度達(dá)到98.7%。有效解決了傳統(tǒng)方法中存在的問(wèn)題，為汽車(chē)螺旋錐齒輪磨損檢測(cè)技術(shù)注入了新鮮血液。