基于proE的車床撥叉加工專用銑床夾具設(shè)計(jì)-含工藝【三維PROE】【銑斷4mm槽】【含6張圖紙及及檔全套】

喜歡就充值下載吧，資源目錄下展示的全都有，下載后全都有，所見即所得，歡迎充值購(gòu)買哦=喜歡就充值下載吧，資源目錄下展示的全都有，下載后全都有，所見即所得，歡迎充值購(gòu)買哦=喜歡就充值下載吧，資源目錄下展示的全都有，下載后全都有，所見即所得，歡迎充值購(gòu)買哦= 裝訂線湖南文理學(xué)院芙蓉學(xué)院工程與技術(shù)系機(jī)械加工工序卡片產(chǎn)品型號(hào)零件圖號(hào)工 序 號(hào)10產(chǎn)品名稱831007零件名稱撥叉工序名稱銑材料毛坯種類同時(shí)加工件數(shù)牌號(hào)硬度鑄件1HT200加工車間切 削 液設(shè) 備機(jī)加工名 稱型 號(hào)夾具名稱及 編 號(hào)輔助工具銑床X5140專用夾具更改內(nèi)容工步號(hào)工 步 內(nèi) 容刀 具量 具走刀長(zhǎng)度（mm）單邊余量（mm）走刀次數(shù)切削深度進(jìn)給量（nn/r或mm/min）主軸轉(zhuǎn)速（r/min）切削速度（m/min）基本時(shí)間（min）粗銑小頭孔的上端面高速鋼盤銑刀高度尺20.92001500.1設(shè)計(jì)者指導(dǎo)教師共頁(yè)第頁(yè)注：大而復(fù)雜的零件允許另用紙繪出工序圖，附在卡片后面。裝訂線湖南文理學(xué)院芙蓉學(xué)院工程與技術(shù)系機(jī)械加工工序卡片產(chǎn)品型號(hào)零件圖號(hào)工 序 號(hào)20產(chǎn)品名稱831007零件名稱撥叉工序名稱銑材料毛坯種類同時(shí)加工件數(shù)牌號(hào)硬度鑄件1HT200加工車間切 削 液設(shè) 備機(jī)加工名 稱型 號(hào)夾具名稱及 編 號(hào)輔助工具銑床X5140專用夾具更改內(nèi)容工步號(hào)工 步 內(nèi) 容刀 具量 具走刀長(zhǎng)度（mm）單邊余量（mm）走刀次數(shù)切削深度進(jìn)給量（nn/r或mm/min）主軸轉(zhuǎn)速（r/min）切削速度（m/min）基本時(shí)間（min）粗銑中間孔的上端面高速鋼盤銑刀高度尺20.92001500.1設(shè)計(jì)者指導(dǎo)教師共頁(yè)第頁(yè)裝訂線湖南文理學(xué)院芙蓉學(xué)院工程與技術(shù)系機(jī)械加工工序卡片產(chǎn)品型號(hào)零件圖號(hào)工 序 號(hào)0產(chǎn)品名稱831007零件名稱撥叉工序名稱銑材料毛坯種類同時(shí)加工件數(shù)牌號(hào)硬度鑄件1HT200加工車間切 削 液設(shè) 備機(jī)加工名 稱型 號(hào)夾具名稱及 編 號(hào)輔助工具銑床X5140專用夾具更改內(nèi)容工步號(hào)工 步 內(nèi) 容刀 具量 具走刀長(zhǎng)度（mm）單邊余量（mm）走刀次數(shù)切削深度進(jìn)給量（nn/r或mm/min）主軸轉(zhuǎn)速（r/min）切削速度（m/min）基本時(shí)間（min）粗銑中間孔的下端面高速鋼盤銑刀高度尺20.92001500.1設(shè)計(jì)者指導(dǎo)教師共頁(yè)第頁(yè)裝訂線湖南文理學(xué)院芙蓉學(xué)院工程與技術(shù)系機(jī)械加工工序卡片產(chǎn)品型號(hào)零件圖號(hào)工 序 號(hào)40產(chǎn)品名稱831007零件名稱撥叉工序名稱銑材料毛坯種類同時(shí)加工件數(shù)牌號(hào)硬度鑄件1HT200加工車間切 削 液設(shè) 備機(jī)加工名 稱型 號(hào)夾具名稱及 編 號(hào)輔助工具銑床X5140專用夾具更改內(nèi)容工步號(hào)工 步 內(nèi) 容刀 具量 具走刀長(zhǎng)度（mm）單邊余量（mm）走刀次數(shù)切削深度進(jìn)給量（nn/r或mm/min）主軸轉(zhuǎn)速（r/min）切削速度（m/min）基本時(shí)間（min）精銑小頭孔的上端面高速鋼盤銑刀高度尺20.92901330.83設(shè)計(jì)者指導(dǎo)教師共頁(yè)第4頁(yè)裝訂線湖南文理學(xué)院芙蓉學(xué)院工程與技術(shù)系機(jī)械加工工序卡片產(chǎn)品型號(hào)零件圖號(hào)工 序 號(hào)50產(chǎn)品名稱831007零件名稱撥叉工序名稱銑材料毛坯種類同時(shí)加工件數(shù)牌號(hào)硬度鑄件1HT200加工車間切 削 液設(shè) 備機(jī)加工名 稱型 號(hào)夾具名稱及 編 號(hào)輔助工具銑床X5140專用夾具更改內(nèi)容工步號(hào)工 步 內(nèi) 容刀 具量 具走刀長(zhǎng)度（mm）單邊余量（mm）走刀次數(shù)切削深度進(jìn)給量（nn/r或mm/min）主軸轉(zhuǎn)速（r/min）切削速度（m/min）基本時(shí)間（min）精銑中間孔的上端面高速鋼盤銑刀高度尺20.92901330.83設(shè)計(jì)者指導(dǎo)教師共頁(yè)第 5頁(yè)裝訂線湖南文理學(xué)院芙蓉學(xué)院工程與技術(shù)系機(jī)械加工工序卡片產(chǎn)品型號(hào)零件圖號(hào)工 序 號(hào)60產(chǎn)品名稱831007零件名稱撥叉工序名稱銑材料毛坯種類同時(shí)加工件數(shù)牌號(hào)硬度鑄件1HT200加工車間切 削 液設(shè) 備機(jī)加工名 稱型 號(hào)夾具名稱及 編 號(hào)輔助工具銑床X5140專用夾具更改內(nèi)容工步號(hào)工 步 內(nèi) 容刀 具量 具走刀長(zhǎng)度（mm）單邊余量（mm）走刀次數(shù)切削深度進(jìn)給量（nn/r或mm/min）主軸轉(zhuǎn)速（r/min）切削速度（m/min）基本時(shí)間（min）精銑中間孔的下端面高速鋼盤銑刀高度尺3.50.41120801.9設(shè)計(jì)者指導(dǎo)教師共頁(yè)第 6頁(yè)裝訂線湖南文理學(xué)院芙蓉學(xué)院工程與技術(shù)系機(jī)械加工工序卡片產(chǎn)品型號(hào)零件圖號(hào)工 序 號(hào)70產(chǎn)品名稱831007零件名稱撥叉工序名稱鉆材料毛坯種類同時(shí)加工件數(shù)牌號(hào)硬度鑄件1HT200加工車間切 削 液設(shè) 備機(jī)加工名 稱型 號(hào)夾具名稱及 編 號(hào)輔助工具鉆床Z525專用夾具更改內(nèi)容工步號(hào)工 步 內(nèi) 容刀 具量 具走刀長(zhǎng)度（mm）單邊余量（mm）走刀次數(shù)切削深度進(jìn)給量（nn/r或mm/min）主軸轉(zhuǎn)速（r/min）切削速度（m/min）基本時(shí)間（min）擴(kuò)兩小端頭孔高速鋼直柄麻花鉆高度尺20.92001000.098設(shè)計(jì)者指導(dǎo)教師共頁(yè)第 7頁(yè)裝訂線湖南文理學(xué)院芙蓉學(xué)院工程與技術(shù)系機(jī)械加工工藝過(guò)程卡片產(chǎn) 品 型 號(hào)撥叉831007零 件 圖 號(hào)產(chǎn) 品 名 稱撥叉零 件 名 稱撥叉材料牌號(hào)HT200毛坯種類鑄件毛坯尺寸每坯件數(shù)2每臺(tái)件數(shù)1工序號(hào)工 序 內(nèi) 容加 工車 間設(shè) 備夾 具刀 具量 具單件工進(jìn)(min)名 稱型 號(hào)名 稱規(guī) 格名 稱規(guī) 格名 稱規(guī) 格10粗銑小頭孔的上端面機(jī)銑床X5140專用夾具高速鋼盤銑刀高度尺20粗銑中間孔的上端面機(jī)銑床X5140專用夾具高速鋼盤銑刀高度尺30粗銑中間孔的下端面機(jī)銑床X5140專用夾具高速鋼盤銑刀高度尺40精銑小頭孔的上端面機(jī)銑床X5140專用夾具高速鋼盤銑刀高度尺50精銑中間孔的上端面機(jī)銑床X5140專用夾具高速鋼盤銑刀高度尺60精銑中間孔的下端面機(jī)銑床X5140專用夾具高速鋼盤銑刀高度尺70擴(kuò)兩小端頭孔機(jī)鉆床Z525專用夾具高速鋼麻花鉆頭高度尺設(shè)計(jì)者指導(dǎo)教師共 頁(yè)第 頁(yè) 機(jī)械原理基于局部平均分解的階次跟蹤分析及其在齒輪故障診斷中的應(yīng)用Junsheng Cheng, Kang Zhang, Yu Yang關(guān)鍵詞： 階次跟蹤分析 局部平均分解 解調(diào) 齒輪 故障診斷摘要：局部平均分解(LMD)是一種新的自適應(yīng)時(shí)頻分析方法,這種方法特別適合處理多分量的調(diào)幅信號(hào)和調(diào)頻(AM-FM)信號(hào)。通過(guò)使用LMD方法,可以將任何復(fù)雜的信號(hào)分解為一系列的產(chǎn)品功能PF分量(PFs),每個(gè)PF分量都是純調(diào)頻信號(hào)和包絡(luò)信號(hào)的乘積，且通過(guò)純調(diào)頻信號(hào)可以獲得具有物理意義的瞬時(shí)頻率。從理論上講,每個(gè)PF分量都是一個(gè)單分量的AM-FM信號(hào)。 因此,可以將LMD的過(guò)程看作是信號(hào)解調(diào)的過(guò)程。齒輪發(fā)生故障時(shí),振動(dòng)信號(hào)呈現(xiàn)明顯的AM-FM特征。因此，針對(duì)齒輪升降速過(guò)程中故障振動(dòng)信號(hào)為多分量的調(diào)制信號(hào)，以及故障特征頻率隨轉(zhuǎn)速變化的特點(diǎn)，提出了一種基于LMD和階次跟蹤分析的齒輪故障診斷方法。齒輪箱的故障診斷實(shí)驗(yàn)表明本文提出的方法能有效地提出齒輪故障診斷特征。1 引言齒輪傳動(dòng)是機(jī)械設(shè)備中常見的傳動(dòng)方式, 故對(duì)齒輪進(jìn)行故障診斷具有重要意義。齒輪故障診斷的關(guān)鍵一步是故障特征的提取。一方面，傳統(tǒng)的齒輪故障診斷方法的重點(diǎn)在一個(gè)固定的旋轉(zhuǎn)速度檢測(cè)振動(dòng)信號(hào)的頻譜分析。 而齒輪作為一種旋轉(zhuǎn)部件, 其升降速過(guò)程的振動(dòng)信號(hào)往往包含了豐富的狀態(tài)信息, 一些在平穩(wěn)運(yùn)行時(shí)不易反映的故障特征在升降速過(guò)程中可能會(huì)充分地表現(xiàn)出來(lái)1,此外，來(lái)自齒輪振動(dòng)信號(hào)的暫態(tài)過(guò)程中，速度依賴性總是顯示非平穩(wěn)特征。如果頻譜分析直接應(yīng)用于非平穩(wěn)振動(dòng)信號(hào)，混頻將不可避免的發(fā)生，這將對(duì)故障特征提取帶來(lái)不良影響。在以往的研究中，為了跟蹤技術(shù)，通常利用振動(dòng)信號(hào)中添加旋轉(zhuǎn)機(jī)械軸轉(zhuǎn)速信息，已經(jīng)成為一個(gè)在旋轉(zhuǎn)機(jī)械故障診斷2,3的重要途徑。從本質(zhì)上講，階次跟蹤分析技術(shù)可以在時(shí)域非平穩(wěn)信號(hào)轉(zhuǎn)換成角域靜止，可以突出的旋轉(zhuǎn)速度相關(guān)的振動(dòng)信息和抑制無(wú)關(guān)的信息。因此，階次跟蹤分析是在助跑過(guò)程中齒輪的故障特征提取和運(yùn)行了一個(gè)可取的方法 另一方面，當(dāng)發(fā)生故障的齒輪振動(dòng)信號(hào)，拿起在運(yùn)行和運(yùn)行過(guò)程中始終存在的振幅特性調(diào)制和頻率調(diào)制（AMFM）。為了提取齒輪故障振動(dòng)信號(hào)的調(diào)制特征，解調(diào)分析是最流行的方法之一 4,5 。然而，傳統(tǒng)的解調(diào)方法，如希爾伯特變換解調(diào)和傳統(tǒng)包絡(luò)分析有其自身的局限性 6 。這些缺點(diǎn)包括兩個(gè)方面：（1）在實(shí)踐中大多數(shù)的齒輪故障振動(dòng)信號(hào)都是多組分是調(diào)頻信號(hào)。這些信號(hào)，在傳統(tǒng)的解調(diào)方法，他們通常是通過(guò)帶通濾波器分解成單組分是調(diào)頻信號(hào)的解調(diào)，然后提取的頻率和振幅信息。然而，這兩個(gè)數(shù)載波頻率的載波頻率成分和幅值都難以在實(shí)踐中被確定，所以帶通濾波器的中心頻率的選擇具有主體性，將解調(diào)誤差和使它提取機(jī)械故障振動(dòng)信號(hào)的特征是無(wú)效的；（2）由于希爾伯特不可避免的窗口效應(yīng)變換，當(dāng)使用希爾伯特變換提取調(diào)制信息，目前的非瞬時(shí)響應(yīng)特性，即，在調(diào)制信號(hào)被解調(diào)以及打破中間部分的兩端會(huì)再次產(chǎn)生調(diào)制，使振幅指數(shù)衰減的方式得到的波動(dòng)，然后解調(diào)誤差將增加 7 。為了克服第一個(gè)缺點(diǎn)，一個(gè)合適的分解方法應(yīng)尋找獨(dú)立的多分量信號(hào)為多個(gè)單組分是調(diào)頻信號(hào)的包絡(luò)分析之前。由于EMD（經(jīng)驗(yàn)?zāi)B(tài)分解）自適應(yīng)復(fù)雜多分量信號(hào)分解為一系列固有模態(tài)函數(shù)（IMF）的瞬時(shí)頻率的物理意義 8,9 ，基于EMD的階比跟蹤方法已廣泛應(yīng)用于齒輪故障診斷 13 。然而，仍然存在許多不足之處 14 ，如在EMD的端點(diǎn)效應(yīng)和模態(tài)混 15 ，仍在進(jìn)行。此外，對(duì)原信號(hào)通過(guò)EMD分解，產(chǎn)生了由希爾伯特變換（上面提到的）缺點(diǎn)是不可避免的在IMF進(jìn)行希爾伯特變換的包絡(luò)分析。此外，有時(shí)無(wú)法解釋的負(fù)瞬態(tài)頻率時(shí)會(huì)出現(xiàn)瞬時(shí)頻率計(jì)算每個(gè)IMF進(jìn)行希爾伯特變換 16 局部均值分解（LMD）是一種新型的解調(diào)分析方法，特別適合于處理多組分的幅度調(diào)制和頻率調(diào)制（AM調(diào)頻）信號(hào) 16 。用LMD，任何復(fù)雜的信號(hào)可以分解成許多產(chǎn)品功能（PFS），每一種產(chǎn)品的包絡(luò)線信號(hào)（獲得直接由分解）的PF瞬時(shí)振幅可以得到一個(gè)純粹的頻率調(diào)制信號(hào)從一個(gè)良好定義的瞬時(shí)頻率可以計(jì)算。在本質(zhì)上，每個(gè)PF正是一種單組分我調(diào)頻信號(hào)。因此，LMD的程序可以，事實(shí)上，作為解調(diào)過(guò)程。調(diào)制信息可以通過(guò)頻譜分析的瞬時(shí)振幅（包絡(luò)信號(hào)，直接獲得通過(guò)分解）每個(gè)PF分量進(jìn)行希爾伯特變換，而不是由PF分量。因此，當(dāng)LMD和EMD方法分別應(yīng)用到解調(diào)分析，與EMD，LMD的突出優(yōu)點(diǎn)是避免希爾伯特變換。此外，LMD迭代過(guò)程中所采用的手段和當(dāng)?shù)氐姆炔黄交牡胤接肊MD的三次樣條的方法，這可能帶來(lái)的包絡(luò)的誤差和影響的精度瞬時(shí)頻率和振幅。此外，與EMD端點(diǎn)效應(yīng)相比并不明顯，因?yàn)樵贚MD方法更快的速度和算法的迭代次數(shù)更少 17 。 基于以上分析，階次跟蹤和解調(diào)技術(shù)，LMD最近的發(fā)展，科學(xué)相結(jié)合，并應(yīng)用于齒輪故障診斷過(guò)程中各軸速度。首先，訂單跟蹤技術(shù)被用于將從時(shí)間域的齒輪振動(dòng)信號(hào)角域。其次，分解角域重采樣信號(hào)的PF系列LMD，因此組件和相應(yīng)的瞬時(shí)振幅和瞬時(shí)頻率可以得到的。最后，進(jìn)行頻譜分析的故障信息含有顯性PF分量的瞬時(shí)幅值。從實(shí)驗(yàn)的振動(dòng)信號(hào)，表明該方法能有效地提取故障特征和分類準(zhǔn)確齒輪工作狀態(tài)的分析結(jié)果。 本文的組織如下。第2節(jié)是一個(gè)給定的LMD方法理論。在第3節(jié)中的齒輪故障診斷方法中，以技術(shù)和LMD跟蹤相結(jié)合的提出和實(shí)踐應(yīng)用表明，提出的方法。此外，LMD和基于EMD的比較也在第3節(jié)提到了基礎(chǔ)的方法。最后，我們得出了第4部分的結(jié)論。2 LMD 方法 LMD方法的本質(zhì)是通過(guò)迭代從原始信號(hào)中分離出純調(diào)頻信號(hào)和包絡(luò)信號(hào),然后將純調(diào)頻信號(hào)和包絡(luò)信號(hào)相乘便可以得到一個(gè)瞬時(shí)頻率具有物理意義的PF分量,循環(huán)處理直至所有的PF分量分離出來(lái)對(duì)任意信號(hào)x(t),其分解過(guò)程如16: ( 1) 確定原始信號(hào)第i個(gè)局部極值及其對(duì)應(yīng)的時(shí)刻,計(jì)算相鄰兩個(gè)局部極值和的平均值 (1) 將所有平均值點(diǎn)mi在其對(duì)應(yīng)的時(shí)間段,內(nèi)伸一線段,然后用滑動(dòng)平均法進(jìn)行0平滑處理,得到局均值m11(t) 。 ( 2) 采用局部極值點(diǎn)計(jì)算局部幅值 : =| -|/2 (2)將所有局部幅值點(diǎn)ai在其對(duì)應(yīng)的時(shí)間段,內(nèi)伸成一條線段,然后采用滑動(dòng)平均法進(jìn)行平滑處理,得到包估計(jì)函數(shù)a11(t) 。 ( 3) 將局部均值函數(shù)m11(t)從原始信號(hào)x(t)中分離來(lái), 即去掉一個(gè)低頻成分,得到 h11(t)=x(t)-m11(t) (3) ( 4)用h11(t)除以包絡(luò)估計(jì)函數(shù)A11( t)以對(duì)h11(t)進(jìn)行解調(diào),得到 s11(t)=h11(t)/A11(t) (4) 對(duì)s11( t)重復(fù)上述步驟便能得到s11(t)的包絡(luò)估計(jì)函數(shù)A12(t),若A12(t)不等于1,則s11( t)不是一個(gè)純調(diào)頻信號(hào)需要重復(fù)上述迭代過(guò)程n次,直至s1n(t)為一個(gè)純調(diào)頻信號(hào),即 s1n(t)的包絡(luò)估計(jì)函數(shù) A1(n+1)(t)=1,所以,有 (5) (6)為理論上, 迭代終止的條件 (7) 在實(shí)踐中,一種變體會(huì)提前確定。如果1a1(n + 1)(t)1 +and1s1n(t)1,然后迭代過(guò)程將停止( 5) 把迭代過(guò)程中產(chǎn)生的所有包絡(luò)估計(jì)函數(shù)相乘便可以得到包絡(luò)信號(hào)( 瞬時(shí)幅值函數(shù)) : （8）( 6) 將包絡(luò)信號(hào)A1(t)和純調(diào)頻信號(hào)s1n(t)相乘便可以得到原始信號(hào)的第一個(gè)PF分量: PF1(t）=a1(t)s1n(t) ( 9)PF1(t)包含了原始信號(hào)中頻率值最高的成分,是一個(gè)單分量的調(diào)幅-調(diào)頻信號(hào),PF1(t)的瞬時(shí)幅值就是包絡(luò)信號(hào)A1(t),PF1(t)的瞬時(shí)頻率f1(t)則可由純調(diào)頻信號(hào)s1n(t)求出,即: (10) ( 7)將第一個(gè)PF分量PF1(t)從原始信號(hào)x(t)中分離出來(lái), 得到一個(gè)新的信號(hào)u1(t),將u1( t)作為原始數(shù)據(jù)重復(fù)以上步驟,循環(huán)k次,直到 uk為一個(gè)單調(diào)函數(shù)為止,即: (11) 原始信號(hào)x(t)能夠被所有的PF分量和uk重構(gòu),即: (12) 產(chǎn)品功能p的數(shù)量在哪里.此外,相應(yīng)的完整的時(shí)頻分布可以通過(guò)組裝瞬時(shí)幅度和瞬時(shí)頻率的PF組件。3 基于階次跟蹤分析與 L M D 的齒輪故障診斷3.1 階次跟蹤分析 階次跟蹤分析首先根據(jù)參考軸的轉(zhuǎn)速信息對(duì)時(shí)域信號(hào)進(jìn)行等角度重采樣, 將時(shí)域非平穩(wěn)信號(hào)轉(zhuǎn)換為角域平穩(wěn)信號(hào), 再對(duì)角域平穩(wěn)信號(hào)進(jìn)行譜分析得到階次譜。階次跟蹤分析能夠提取信號(hào)中與參考軸轉(zhuǎn)速有關(guān)的信息, 同時(shí)抑制與轉(zhuǎn)速無(wú)關(guān)的信號(hào), 因此非常適合分析旋轉(zhuǎn)機(jī)械在變轉(zhuǎn)速過(guò)程下的振動(dòng)信號(hào)。實(shí)現(xiàn)階次跟蹤分析技術(shù)的關(guān)鍵在于, 如何實(shí)現(xiàn)被分析信號(hào)相對(duì)于參考軸的等角度重采樣, 即階次重采樣。常用的階次重采樣方法有硬件階次跟蹤法 6、計(jì)算階次跟蹤法 7和基于瞬時(shí)頻率估計(jì)的階次跟蹤法 8等。硬件階次跟蹤法直接通過(guò)專用的模擬設(shè)備實(shí)現(xiàn)信號(hào)的等角度重采樣,實(shí)時(shí)性好,但只適用于軸轉(zhuǎn)速較穩(wěn)定的情況,且成本很高;基于瞬時(shí)頻率估計(jì)的階次跟蹤法不需要專門的硬件設(shè)備,無(wú)需考慮硬件安裝問題,且成本較低, 但是不適用于分析多分量信號(hào),而實(shí)際工程信號(hào)大多為多分量信號(hào), 因此其實(shí)際應(yīng)用意義不大;COT法通過(guò)軟件的形式實(shí)現(xiàn)等角度重采樣,分析精度高, 對(duì)被分析的信號(hào)沒有特別的要求,并且無(wú)需特定的硬件, 因此是一種應(yīng)用廣泛的階次跟蹤分析方法。根據(jù)試驗(yàn)條件采用COT法實(shí)現(xiàn)信號(hào)的階次重采樣,其具體步驟如下:1. 對(duì)振動(dòng)信號(hào)和轉(zhuǎn)速信號(hào)分兩路同時(shí)進(jìn)行等時(shí)間間隔(間隔為$t)采樣,得到異步采樣信號(hào);2. 通過(guò)轉(zhuǎn)速信號(hào)計(jì)算等角度增量 $H 所對(duì)應(yīng)的時(shí)間序列ti ;3. 根據(jù)時(shí)間序列ti的值,對(duì)振動(dòng)信號(hào)進(jìn)行插值,求出其對(duì)應(yīng)的幅值,得到振動(dòng)信號(hào)的同步采樣信號(hào),即角域平穩(wěn)信號(hào);4.使用LMD分解平衡角重采樣信號(hào),因此sPF系列組件和相應(yīng)的瞬間振幅和瞬時(shí)頻率可以獲得5.光譜分析應(yīng)用于每個(gè)PF的瞬時(shí)振幅組件,然后我們有訂單譜3.2 齒輪故障診斷實(shí)例升降速過(guò)程中的齒輪故障振動(dòng)信號(hào)通常是多分量的調(diào)幅-調(diào)頻信號(hào),并且故障特征頻率會(huì)隨著轉(zhuǎn)速的變化而改變。針對(duì)升降速過(guò)程齒輪故障振動(dòng)信號(hào)的這些特點(diǎn), 提出了基于階次跟蹤分析和 LM D 的齒輪故障診斷方法。首先采用階次跟蹤分析將齒輪升降速過(guò)程的時(shí)域振動(dòng)信號(hào)轉(zhuǎn)換成角域平穩(wěn)信號(hào);然后對(duì)角域信號(hào)進(jìn)行LMD分解,得到一系列PF分量,以及各個(gè)PF分量的瞬時(shí)幅值和瞬時(shí)頻率; 最后對(duì)各個(gè)PF分量的瞬時(shí)幅值進(jìn)行頻譜分析,便可以有效地提取出齒輪故障特征。為了驗(yàn)證方法的正確性,在旋轉(zhuǎn)機(jī)械試驗(yàn)臺(tái)上進(jìn)行了齒輪正常和齒根裂紋兩種工況的試驗(yàn)。該系統(tǒng)中, 電機(jī)輸入軸齒輪齒數(shù)z1=55, 輸出軸齒輪齒數(shù)z2 = 75。在輸入軸齒輪齒根上加工出小槽,以模擬齒根紋故 障, 因此齒輪嚙合階次xm=55,故障特征階次xc=1。圖1和圖2所示分別為由轉(zhuǎn)速傳感器測(cè)得的輸入軸瞬時(shí)轉(zhuǎn)速n(t),以及由振動(dòng)傳感器測(cè)得的齒輪故障 振動(dòng)加速度a(t),其中采樣頻率為8192H z,采樣時(shí)間為20s從圖1可以看出,輸入軸轉(zhuǎn)速首先從150r/min逐漸加速至1410r/min, 然后再減速到820r/min,而加速度信號(hào)的幅值也隨著作出了相應(yīng)的變化。不失一般性,截取圖2中5 7s升速過(guò)程的信號(hào) a1(t)進(jìn)行分析。 圖 1 輸 入軸的瞬時(shí)轉(zhuǎn)速 n ( t )圖 2 齒輪故障振動(dòng)加速度信號(hào) a( t )值在秩序O=55和O=110相應(yīng)的齒輪嚙合秩序和雙。因此這意味著頻率混淆現(xiàn)象已經(jīng)在很大程度上消除。然而,為j1()仍然是一個(gè)多個(gè)組件MA-MF信號(hào)。因此,一邊頻帶反映故障特征頻率模糊。有效地提取故障特征,應(yīng)用LMD j - 1(),因此七PF組件和殘?jiān)梢缘玫綀D6所示,這意味著LMD解調(diào)的進(jìn)展。因此,它是可以提取齒輪故障特性,利用頻譜分析的瞬時(shí)振幅PF組件包含主要故障信息。通過(guò)分析,我們知道失敗的主要信息包括在第一個(gè)PF組件。因此,無(wú)花果。7和8給瞬時(shí)振幅a1()的第一個(gè)PF組件PF 1()和相應(yīng)的秩序光譜的a1(),很明顯,有不同的光譜峰值在第一順序(O = 1)對(duì)應(yīng)齒輪階次跟蹤功能,符合齒輪的實(shí)際工況。圖9和圖10顯示轉(zhuǎn)速信號(hào)的n(t)和振動(dòng)加速度信號(hào)的時(shí)域波形s(t)齒輪分別與破碎的牙齒,采樣率為8192 Hz和總樣品時(shí)間是20年代。斷齒故障引入輸入軸上的齒輪與激光切割槽的牙根。首先,一段信號(hào)s1(t)5 s-7年代為進(jìn)一步分析的進(jìn)步是攔截;其次,假設(shè)樣本點(diǎn)每旋轉(zhuǎn)400;第三,角域信號(hào)為j1()圖11所示可以通過(guò)執(zhí)行命令重采樣s1(t);第四,LMD適用于j-1();最后,相應(yīng)的秩序頻譜圖12所示的瞬時(shí)振幅首先PF組件PF 1()可以了,很明顯,有不同的光譜峰值(比在圖8)在第一順序(O = 1)階次跟蹤分析對(duì)應(yīng)于齒輪故障功能,符合齒輪的實(shí)際工況。同樣的,我們同樣可以做正常的齒輪。轉(zhuǎn)速信號(hào)n(t)和振動(dòng)的時(shí)域波形加速度信號(hào)s(t)的正常齒輪分別列在無(wú)花果。13和14,采樣率為8192 Hz和總樣品時(shí)間是20多歲。在上述相同的方法應(yīng)用于原始信號(hào)圖14所示,結(jié)果無(wú)花果所示。15和16。圖15顯示了角域j - 1()執(zhí)行順序重采樣后的信號(hào)部分(5s-7年代在籌備進(jìn)展)的原始信號(hào)。圖16顯示了相應(yīng)的瞬時(shí)振幅譜第一個(gè)PF組件,很難找到齒輪故障特征,也符合實(shí)際的工作狀態(tài)的裝備。目前,多組分的另一個(gè)競(jìng)爭(zhēng)解調(diào)方法AM-FM信號(hào),即經(jīng)驗(yàn)?zāi)Ｊ椒纸?EMD)存在,已經(jīng)被廣泛應(yīng)用于信號(hào)解調(diào)分析(7、22)。為了比較兩個(gè)EMD方法,取代LMD,我們能做的同樣使用EMD進(jìn)行重采樣信號(hào)無(wú)花果所示。圖4、11和15 圖 3 齒輪故障振動(dòng)加速度信號(hào)的頻譜圖 4 階次重采樣后的齒輪故障振動(dòng) 加速度信號(hào)圖5 j1()的階次譜分別,因此可以獲得一系列國(guó)際貨幣基金組織(IMF)組件。此外,相應(yīng)的瞬時(shí)振幅和國(guó)際貨幣基金組織每個(gè)組件的瞬時(shí)頻率可以通過(guò)希爾伯特變換計(jì)算。通過(guò)分析,我們知道,IMF主要特征信息包含在第一個(gè)組件。因此,只有應(yīng)用于瞬時(shí)頻譜分析第一個(gè)國(guó)際貨幣基金組織(IMF)組件的振幅。無(wú)花果。17日至19日給訂單頻譜對(duì)應(yīng)三種振動(dòng)信號(hào)的破解斷層、斷齒故障和正常的齒輪,分別,很明顯,訂單跟蹤分析基于EMD也可以提取齒輪故障特性,確定齒輪的工作狀態(tài)。盡管EMD和LMD都可以分解原始信號(hào)實(shí)際上,兩種方法之間的差異仍然存在。EMD方法比較,如第一節(jié)中所述,LMD有更多迭代次數(shù)少等優(yōu)點(diǎn),不明顯的效果和更少的瞬時(shí)頻率的虛假成分,可以使用更多的應(yīng)用在實(shí)踐中。 圖 6 角域信號(hào)j1( )的LMD分解結(jié)果圖 7 PF1()的瞬時(shí)幅值A(chǔ)1()圖 8 第1個(gè)PF分量的幅值譜圖 9 輸入軸的瞬時(shí)轉(zhuǎn)速 n(t)圖 1 0 正常齒輪的振動(dòng)加速度信號(hào) a(t)圖11 階次重采樣后的正常齒輪振動(dòng)加速度信號(hào)j1(）圖 12 第一個(gè)PF分量的幅值譜圖13 輸入軸轉(zhuǎn)速r(t)正常齒輪前和過(guò)程中圖圖14 齒輪的振動(dòng)加速度信號(hào)(t)在正常狀態(tài)圖15 相應(yīng)的振動(dòng)加速度信號(hào)為j1()角域通過(guò)應(yīng)用順序重采樣tos(t)圖14所示。圖17 第一個(gè)IMF分量的幅值譜 圖 18 第一個(gè)IMF分量的幅值譜3 結(jié)論 在齒輪故障診斷技術(shù)、階次跟蹤是一個(gè)著名的技術(shù),可用于故障檢測(cè)的旋轉(zhuǎn)機(jī)器采用振動(dòng)信號(hào)。針對(duì)齒輪故障振動(dòng)信號(hào)的調(diào)制特點(diǎn)在助跑和破敗的和缺點(diǎn)在齒輪經(jīng)常可以發(fā)相關(guān)軸轉(zhuǎn)速在瞬態(tài)過(guò)程中,階次跟蹤和技術(shù)LMD相結(jié)合用于齒輪故障診斷。從理論分析和實(shí)驗(yàn)結(jié)果以下幾點(diǎn)得出結(jié)論: ( 1) 在分析齒輪變轉(zhuǎn)速狀態(tài)下的振動(dòng)信號(hào)時(shí),轉(zhuǎn)速波動(dòng)會(huì)引起頻譜圖出現(xiàn)頻率混疊, 而階次跟蹤分析通過(guò)對(duì)信號(hào)進(jìn)行階次重采樣能夠在很大程度上消除頻率混疊, 使頻譜圖的譜線清晰可讀。 ( 2) 齒輪故障時(shí)的振動(dòng)信號(hào)為一多分量的調(diào)幅- 調(diào)頻信號(hào), 采用LMD方法能將其分解為若干個(gè)PF分量之和,同得到各個(gè)PF分量的瞬時(shí)幅值和瞬時(shí)頻率, 實(shí)現(xiàn)了原信號(hào)的解調(diào)。對(duì)含有齒輪故障特征的PF分量的瞬時(shí)幅值進(jìn)行頻譜分析, 能夠準(zhǔn)確地提取出齒輪故障特征信息。 圖19 階次的第一個(gè)國(guó)際貨幣基金組織(IMF)組件的正常使用EMD齒輪 ( 3) 對(duì)齒輪正常和齒根裂紋兩種工況的振動(dòng)信號(hào)進(jìn)行了分析,分析結(jié)果表明, 本文方法能夠準(zhǔn)確地反映出齒輪的實(shí)際工況。References1 S.K. Lee, P.R. White, Higher-order timefrequency analysis and its application to fault detection in rotating machinery, Mechanical Systems and Signal Processing 11 (1997) 637650.2 Mingsian Bai, Jiamin Huang, Minghong Hong, Fucheng Su, Fault diagnosis of rotating machinery using an intelligent order tracking system, Journal of Sound and Vibration 280 (2005) 699718.3 JianDa Wu, YuHsuan Wang, PengHsin Chiang, Mingsian R. Bai, A study of fault diagnosis in a scooter using adaptive order tracking technique and neural network, Expert Systems with Applications 36 (1) (2009) 4956.4 J. Ma, C.J. Li, Gear defect detection through model-based wideband demodulation of vibrations, Mechanical System and Signal Process 10 (5) (1996) 653665.5 R.B. Randall, J. Antoni, S. chobsaard, The relationship between spectral correlation and envelope analysis in the diagnostics of bearing faults and other cyclostationary machine signals, Mechanical Systems and Signal Processing 15 (5) (2001) 945962.6 He Lingsong, Li Weihua, Morlet wavelet and its application in enveloping, Journal of Vibration Engineering. 15 (1) (2002) 119122.7 Cheng Junsheng, Yu Dejie, Yang Yu, The application of energy operator demodulation approach based on EMD in machinery fault diagnosis, Mechanical Systems and Signal Processing 21 (2) (2007) 668677.8 N.E. Huang, Z. Shen, S.R. Long, The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis, Proceedings of the Royal Society of London Series 454 (1998) 903995.9 N.E. Huang, Z. Shen, S.R. Long, A new view of nonlinear water waves: the Hilbert spectrum, Annual Review of Fluid Mechanics 31 (1999) 417457.10 B.L. Eggers, P.S. Heyns, C.J. Stander, Using computed order tracking to detect gear condition aboard a dragline, Journal of the Southern AfricanInstitute of Mining and Metallurgy 107 (2007) 18.11 Q. Gao, C. Duan, H. Fan, Q. Meng, Rotating machine fault diagnosis using empirical mode decomposition, Mechanical Systems and Signal Processing 22 (2008) 10721081.12 F.J. Wu, L.S. Qu, Diagnosis of subharmonic faults of large rotating machinery based on EMD, Mechanical Systems and Signal Processing 23 (2009) 467475.13 K.S. Wang, P.S. Heyns, Application of computed order tracking, VoldKalman filtering and EMD in rotating machine vibration, Mechanical Systems and Signal Processing 25 (2011) 416430.14 Junsheng Cheng, Dejie Yu, Yu Yang, Application of support vector regression machines to the processing of end effects of HilbertHuang transform, Mechanical Systems and Signal Processing 21 (3) (2007) 11971211.15 Marcus Datig, Torsten Schlurmann, Performance and limitations of the HilbertHuang transformation (HHT) with an application to irregular water waves, Ocean Engineering 31 (14) (2004) 17831834.16 Jonathan S. Smith, The local mean decomposition and its application to EEG perception data, Journal of the Royal Society, Interface 2 (5) (2005) 443454.17 Junsheng Cheng, Yi Yang, Yu Yang A rotating machinery fault diagnosis method based on local mean decomposition, Digital Signal Processin 22 (2) (2012) 356366.18 K.M. Bossley, R.J. Mckendrick, Hybrid computed order tracking, Mechanical Systems and Signal Processing 13 (4) (1999) 627641.19 JianDa Wu, Mingsian R. Bai, Fu Cheng Su, Chin Wei Huang, An expert system for the diagnosis of faults in rotating machinery using adaptive order tracking algorithm, Expert Systems with Applications 36 (3) (2009) 54245431.20 Guo Yu, Qin Shuren, Tang Baoping, Ji Yuebo, Order tracking of rotating machinery based on instantaneous frequencies estimation, Chinese Journalof Mechanical Engineering. 39 (3) (2003) 3236.21 Yu Dejie, Yang Yu, Cheng Junsheng, Application of timefrequency entropy method based on HilbertHuang transform to gear fault diagnosis, Measurement 40 (2007) 823830.22 R.T. Rato, M.D. Ortigueira, A.G. Batista, On the HHT, its problems, and some solutions, Mechanical Systems and Signal Processing 22 (6) (2008) 13741394.An order tracking technique for the gear fault diagnosis using local meandecomposition methodJunsheng Cheng, Kang Zhang, Yu YangState Key Laboratory of Advanced Design and Manufacturing for Vehicle Body, Hunan University, Changsha, 410082, PR ChinaCollege of Mechanical and Vehicle Engineering, Hunan University, Changsha, 410082, PR Chinaa r t i c l ei n f oa b s t r a c tArticle history:Received 17 November 2010Received in revised form 13 December 2011Accepted 30 April 2012Available online 28 May 2012Local mean decomposition (LMD) is a new self-adaptive timefrequency analysis method,which is particularly suitable for the processing of multi-component amplitude-modulatedand frequency-modulated (AMFM) signals. By using LMD, any complicated signal can bedecomposed into a number of product functions (PFs), each of which is the product of anenvelope signal and a purely frequency modulated signal from which physically meaningfulinstantaneous frequencies can be obtained. Theoretically, each PF is exactly a mono-componentAMFM signal. Therefore, the procedure of LMD can be regarded as the process of demodulation.While fault occurs in gear, the vibration signals would exactly present AMFM characteristics.Therefore, targeting the modulation feature of gear fault vibration signal in run-ups and run-downs and the fact that fault characteristics found in gear vibration signal could often be relatedto revolution of the shaft in the transient process, a gear fault diagnosis method in which ordertracking technique and local mean decomposition is put forward. The analysis results from thepractical gearbox vibration signal demonstrate that the proposed algorithm is effective in gearfault feature extraction. 2012 Elsevier Ltd. All rights reserved.Keywords:Order tracking techniqueLocal mean decompositionDemodulationGearFault diagnosis1. IntroductionGears are the important and frequently encountered components in the rotating machines that find widespread industrialapplications. Therefore, the corresponding gear fault diagnosis has been the subject of extensive research.The key step of gear fault diagnosis is the extraction of fault feature. On the one hand, the conventional gear fault diagnosismethods focus on examining the frequency spectrum analysis of vibration signal at a fixed rotation speed. Unfortunately, theinformation obtained thus is only partial because some faults maybe do not respond significantly at the fixed operation speed.Since faults commonly found in gear could often be related to revolution of the shaft, more comprehensive information may beacquired by measuring the gear vibration signal in the process of run-up and run-down 1. In addition, vibration signals derivedfrom gear in the transient process that are speed-dependent always display non-stationary feature. If frequency spectrum analysisis directly applied to the non-stationary vibration signal, frequency mixing would occur inevitably, which will bring undesirableeffect to the fault feature extraction. In past research, order-tracking technique, which normally exploits a vibration signalsupplemented with information of shaft speed of rotating machinery, has become one of the significant approaches for faultdiagnosis in rotating machinery 2,3. Essentially, order-tracking technique can transform a non-stationary signal in time domaininto stationary one in angular domain, which can highlight the vibration information related to rotation speed and restrain theunrelated information. Therefore, order tracking is a desirable method to extract gear fault feature in the process of run-up andrun-down.Mechanism and Machine Theory 55 (2012) 6776 Corresponding author at: State Key Laboratory of Advanced Design and Manufacturing for Vehicle Body, Hunan University, Changsha, 410082, PR China.Tel.: +86 731 88664008; fax: +86 731 88711911.E-mail address: (J. Cheng).0094-114X/$ see front matter 2012 Elsevier Ltd. All rights reserved.doi:10.1016/j.mechmachtheory.2012.04.008Contents lists available at SciVerse ScienceDirectMechanism and Machine Theoryjournal homepage: the other hand, while faults occur in gears, the vibration signal picked up in run-up and run-down process always presentthe characteristics of amplitude-modulated and frequency-modulated (AMFM). In order to extract the modulation featureof gear fault vibration signals, demodulation analysis is one of the most popular methods 4,5. However, conventionaldemodulation approaches such as Hilbert transform demodulation and traditional envelope analysis have their own limitations6. These drawbacks include two aspects: (1) in practice most gear fault vibration signals are all multi-component AMFMsignals. For these signals, in conventional demodulation approaches, they are usually decomposed into single component AMFMsignals by band-pass filter and then demodulated to extract frequencies and amplitudes information. However, both the numberof the carrier frequency components and the magnitude of the carrier frequency are hard to be determined in practice, so theselection of central frequency of band-pass filter carries great subjectivity that would bring demodulation error and make itineffective to extract the characteristic of machinery fault vibration signal; (2) owing to the inevitable window effect of Hilberttransform, when Hilbert transform is used to extract the modulate information, the demodulation results present non-instantaneous response characteristic, that is, at the two ends of the modulated signal which has been demodulated as well as themiddle part with break would produce modulation again, which makes the amplitude get fluctuation in an exponentialattenuation way, and then the demodulation error would increase 7. In order to overcome the first drawback, an appropriatedecomposition method should be looked for to separate multi-component signal into a number of single component AMFMsignals before the envelope analysis. Since EMD (Empirical mode decomposition) could adaptively decompose a complicatedmulti-component signal into a sum of intrinsic mode functions (IMFs) whose instantaneous frequencies have physicalsignificance 8,9, order tracking method based on EMD has been widely used in the gear fault diagnosis 1013. However, therestill exist many deficiencies in EMD such as the end effects 14 and modes mixing 15 that are still underway. In addition, afterthe original signal is decomposed by EMD, the drawback produced by Hilbert transform (above mentioned) is inevitable whenIMF is performed envelope analysis by Hilbert transform. Moreover, sometimes the unexplainable negative instantaneousfrequency would appear when calculating instantaneous frequency by performing Hilbert transform to each IMF 16.Local mean decomposition (LMD) is a novel demodulation analysis method, which is particularly suitable for the processing ofmulti-component amplitude-modulated and frequency-modulated (AMFM) signals 16. By using LMD, any complicated signalcan be decomposed into a number of product functions (PFs), each of which is the product of an envelope signal (obtaineddirectly by the decomposition) from which instantaneous amplitude of the PF can be obtained and a purely frequency modulatedsignal from which a well-defined instantaneous frequency could be calculated. In essence, each PF is exactly a mono-componentAMFM signal. Therefore, the procedure of LMD could be, in fact, regarded as the process of demodulation. Modulationinformation can be extracted by performing spectrum analysis to the instantaneous amplitude (envelope signal, obtained directlyby the decomposition) of each PF component rather than by performing Hilbert transform to the PF components. Hence, whenLMD and EMD are applied to the demodulation analysis respectively, compared with EMD, the prominent advantage of LMD is toavoid the Hilbert transform. In addition, the LMD iteration process which uses smoothed local means and local magnitudes avoidsthe cubic spline approach used in EMD, which maybe bring the envelope errors and influence on the precision of theinstantaneous frequency and amplitude. Moreover, compared with EMD the end effect is not obvious in LMD approach because offaster algorithm speed and less iterative times 17.Based upon the above analysis, order-tracking analysis and the recent development of demodulation techniques, LMD, arecombined and applied to the gear fault diagnosis of various shaft speeds process. Firstly, order tracking technique is used totransform the gear vibration signals from time domain to angular domain. Secondly, decompose the re-sampling signal of angulardomain by LMD, thus s series PF components and corresponding instantaneous amplitudes and instantaneous frequencies can beobtained. Finally, spectrum analysis is carried out to the instantaneous amplitudes of the PF component containing dominant faultinformation. The analysis results from the experimental vibration signal show that the proposed method can extract fault featureof the gear effectively and classify working condition accurately.This paper is organized as follows. A theory of the LMD approach is given in Section 2. In Section 3 a gear fault diagnosisapproach in which order tracking technique and LMD are combined is put forward and the practice applications of proposedmethod are demonstrated. In addition, the comparison between LMD-based and EMD-based method is also given in Section 3.Finally, we offer the conclusion in Section 4.2. LMD analysis methodAs mentioned above, the nature of LMD is to demodulate AMFM signals. By using LMD a complicated signal can bedecomposed into a set of product functions, each of which is the product of an envelope signal and a purely frequency modulatedsignal. Furthermore, the completed timefrequency distribution of the original signal can be obtained. For any signal x(t), it can bedecomposed as follows 16:(1) Determine all local extrema niof the original signal x(t), and then the mean value miof two successive extrema niand ni+1can be calculated bymini ni121All mean value miof two successive extreme are connected by straight lines, and then local mean function m11(t)can be formed by using moving averaging to smooth the local means mi.68J. Cheng et al. / Mechanism and Machine Theory 55 (2012) 6776(2) A corresponding envelope estimate aiis given byainini1?22Similarly, the envelope estimate aiis smoothed in the same way and the corresponding envelope function a11(t) isformed.(3) The local mean function m11(t) is subtracted from the original signal x(t) and the resulting signal h11(t) is given byh11t x t m11t 3(4) h11(t) can be amplitude demodulated by dividing it by envelope function a11(t)s11t h11t =a11t 4Ideally, s11(t) is a purely frequency modulated signal, namely, the envelope function a12(t) of s11(t) should satisfya12(t)=1. If a12(t)1, then s11(t) is regarded as the original signal and the above procedure needs to be repeateduntil a purely frequency modulated signal s1n(t) that meets 1s1n(t)1 is derived. In other words, envelopefunction a1(n+1)(t) of the resulting s1n(t) should satisfy a1(n+1)(t)=1. Thereforeh11t x t m11t h12 s11t m12t h1nt s1 n1t m1nt 8:5in which,s11t h11t =a11t s12t h12t =a12t s1nt h1nt =a1nt 8:6where the objective is thatlimna1nt 17In practice, a variation can be determined in advance. If 1a1(n+1)(t)1+ and 1s1n(t)1, then iterativeprocess would be stopped.(5) Envelope signal a1(t), namely, instantaneous amplitude function, can be derived by multiplying together the successiveenvelope estimate functions that are acquired during the iterative process described above.a1t a11t a12t a1nt nq1a1qt 8where q is the times of the iterative process.(6) Multiplying envelope signal a1(t) by the purely frequency modulated signal s1n(t) the first product function PF1of theoriginal signal can be obtained.PF1t a1t s1nt 9PF1contains the highest frequency oscillations of the original signal. Meantime, it is a mono-component AMFMsignal, whose instantaneous amplitude is exactly the envelope signal a1(t) and instantaneous frequency is definedfrom the purely frequency modulated signal s1n(t) asf1t 12d arccos s1nt ?dt10(7) Subtract the first PF component PF1(t) from the original signal x(t) and we have a new signal u1(t), which becomes the neworiginal signaland the whole of the above procedure is repeated,i.e. up tok times,until ukbecomes monotonic functionu1t x t PF1t u2t u1t PF2t ukt uk1t PFkt 8:1169J. Cheng et al. / Mechanism and Machine Theory 55 (2012) 6776Thus, the original signal x(t) was decomposed into k-product and a monotonic function ukx t Xkp1PFpt ukt 12where p is the number of the product function.Furthermore, the corresponding complete timefrequency distribution could be obtained by assembling the instantaneousamplitude and instantaneous frequency of all PF components.3. The gear fault diagnosis method based on order tracking technique and LMD3.1. Order tracking analysis and the corresponding fault diagnosis methodOrder-tracking technique could transform a non-stationary signal in time domain into a stationary signal in angular domain byapplying equi-angular re-sampling to vibration signal with reference to shaft speed. Furthermore, order spectrum can be obtainedby using spectrum analysis to stationary signal in angular domain, thus the information related to rotation speed can behighlighted and the unrelated one could be restrained. Therefore, order-tracking is suitable for the vibration signal analysis ofrotation machine.There are three popular techniques for producing synchronously sampled data: a traditional hardware solution, computedorder tracking (COT) and order tracking based on estimation of instantaneous frequency 1820. The traditional hardwareapproach, which uses specialized hardware to dynamically adapt the sample rate, is only suitable for the case that rotating speedof shaft is relatively smooth, thus resulting to a high cost. The method of order tracking based on estimation of instantaneousfrequency has no need for specialized hardware and thus cost is relatively low, however, it has failed to analyze multiplecomponent signal. While in practice most gear fault vibration signals exactly present the characteristic of multi-component.Therefore, this technique has little practice significance. COT technique realized equi-angular re-sampling by software, thereforeit not only requires no specialized hardware, but also have no limitation for analysis signal that means it is more flexible and moreaccurate. Just for this reason, COT is introduced into the gear fault detection in this paper.The step of the gear fault diagnosis method based on order tracking technique and LMD can be listed as follows:(1) The vibration signals and a tachometer signal are asynchronously sampled, that is, they are sampled conventionally atequal time incrementst;(2) Calculate the time series ticorresponding to equi-angular increments by tachometer signals;(3) According to the time series ti, apply interpolation to the vibration signals, thus the synchronous sampling signal, namely,stationary signal in angular domain, can be obtained;(4) Use LMD to decompose the equi-angular re-sampling signal, thus s series PF components and corresponding instantaneousamplitudes and instantaneous frequencies can be acquired;(5) Apply spectrum analysis to the instantaneous amplitude of each PF component, and then we have the order spectrum.3.2. ApplicationSince the gear fault vibration signal in run-up and run-down process are always multiple component AMFM signals and faultfeature frequency would vary with rotation speed, the fault diagnosis method in which order tracking technique and LMD arecombined would be suitable for gear fault detection.To verify the effectiveness of the proposed method, the fault diagnosis method based on order tracking technique and LMDwas applied to the experimental gear vibration signals analysis. An experiment has been carried out on the rotating machinerytest rig that is used for modeling different gear faults 21. Here we consider three working conditions that are gear with normalcondition, with cracked tooth and with broken tooth. Standard gears with teeth number z=55 and z=75 are used on input andoutput shafts respectively, in which the crack fault is introduced into the gear on the input shaft by cutting slot with laser in theroot of tooth, and the width of the slot is 0.15 mm, as well as its depth is 0.3 mm. Therefore, the mesh order is xm=55 and thefault feature order is xc=1. Figs. 1 and 2 give the rotation speed signal r(t) picked up by a tachometer and vibration accelerationsignal s(t) of the gear with crack fault collected by a piezoelectric acceleration sensor respectively, in which the sample frequencyis 8192 Hz and total sample time is 20 s, and from which we know the speed of input shaft increased gradually from 150 rpm to1410 rpm, then decreased to 820 rpm. Meantime, the amplitude of vibration acceleration signal accordingly changed, from whicha section of signal s1(t) of 5 s7 s in the run-up progress is intercepted for further analysis. Fig. 3 gives the spectrum of s1(t) byapplying spectrum analysis directly to vibration signal. For the rotation speed changes with time, the frequency mixing arises.Therefore, it is impossible to find meshing frequency and fault feature frequency in Fig. 3. As a result, actual gear workingcondition cannot be identified. Replace direct spectrum analysis by the order tracking method. Firstly, assume sample point perrotation is 400, namely, the maximum analysis order is 200. Secondly, angular domain signal j1() shown in Fig. 4 can be obtainedby performing order re-sampling to s1(t), in which horizontal ordinate has changed from time to radian. Thirdly, thecorresponding order spectrum of j1() can be calculated that is illustrated in Fig. 5, from which we can find obvious spectral peak70J. Cheng et al. / Mechanism and Machine Theory 55 (2012) 6776values at order O=55 and O=110 corresponding to gear meshing order and the double. Thus it means that frequency aliasingphenomenon has been eliminated to a large degree. However, j1() is still a multiple component MAMF signal. Therefore, sidefrequency band reflecting fault feature frequency is indistinct. To extract fault characteristic effectively, apply LMD to j1(), thusseven PF components and a residue can be obtained shown in Fig. 6, which means LMD is a demodulation progress. Therefore, it ispossible to extract gear fault feature by utilizing spectrum analysis to the instantaneous amplitude of PF component containingdominant fault information. By analysis, we know that the main failure information is included in the first PF component.Therefore, Figs. 7 and 8 give instantaneous amplitude a1() of the first PF component PF1() and the corresponding orderspectrum of a1(), from which it is clear that there are distinct spectral peak value at the 1st order (O=1) corresponding to gearfault feature order xc, which accords with the actual working condition of the gear.Figs. 9 and 10 show the rotation speed signal n(t) and the time domain waveform of vibration acceleration signal s(t) of thegear with broken tooth respectively, in which the sample rate is 8192 Hz and total sample time is 20 s. The broken tooth fault isintroduced into the gear on the input shaft by cutting slot with laser in the root of tooth. Firstly, a section of signal s1(t) of 5 s7 sin the run-up progress is intercepted for further analysis; secondly, assume sample point per rotation is 400; thirdly, angulardomain signal j1() shown in Fig. 11 can be obtained by performing order re-sampling to s1(t); fourthly, apply LMD to j1();finally, the corresponding order spectrum shown in Fig. 12 of instantaneous amplitude of the first PF component PF1() can beacquired, from which it is clear that there are distinct spectral peak value (it is bigger than that in Fig. 8) at the 1st order (O=1)corresponding to gear fault feature order xc, which accords with the actual working condition of th