你知道怎么用LTspice分析振動數(shù)據(jù)嗎？

數(shù)字技術(shù)的進步?jīng)]有絲毫放緩的跡象，已滲透到我們生活的方方面面。為機器提供智能并非奧威爾式的反烏托邦；由于自動化反饋環(huán)路有助于減少直接維護時間，因此可提高工廠自動化的效率。

工業(yè)4.0描述了將大數(shù)據(jù)的優(yōu)勢帶入工廠車間的概念。裝有傳感器的機器可監(jiān)控自身的性能并相互通信，從而共同分擔整個工作載荷，同時向后臺提供重要的診斷信息，而且無論在同一座建筑物里還是在不同的大陸都可以實現(xiàn)。

對ADI產(chǎn)品進行的一項快速調(diào)查顯示，ADI主要致力于為工業(yè)物聯(lián)網(wǎng) (IIOT)提供解決方案，即從傳感器到云的各種穩(wěn)定可靠的高性能信號鏈組件。

在工業(yè)自動化中的一個應用領(lǐng)域就是狀態(tài)監(jiān)控(CbM),通過仔細校準機器的標稱工作特性，然后使用本地傳感器密切監(jiān)控機器本身的狀態(tài)。偏離標稱信號的狀態(tài)即表示機器需要維護。因此，配備狀態(tài)監(jiān)控系統(tǒng)的機器可根據(jù)實際需要進行維護，而不是相對隨意地安排維修計劃。

確定電機運行狀態(tài)的一種比較好的方法是檢查其振動特征。ADI的MEMS技術(shù)可用于持續(xù)監(jiān)控電機的振動特征，并與已知無故障電機的特征比較，由此判斷電機的運行狀況。事實上，每種電機故障都有其自己的獨特諧波特征。通過查看振動模式的諧波成分，可以檢測軸承、內(nèi)環(huán)和外環(huán)，甚至齒輪箱齒中的故障。

為了產(chǎn)生用于在LTspice中進行傅里葉分析的數(shù)據(jù)，將三個ADXL1002加速度計連接到電機，以測量側(cè)向、垂向和縱向（分別為X、Y和Z）振動。

將振動數(shù)據(jù)下載并保存到Microsoft Excel電子表格中。在500 kSPS速率下進行數(shù)據(jù)采樣，通過一秒振動數(shù)據(jù)得到三列Microsoft Excel數(shù)據(jù)，每列數(shù)據(jù)長500,000行。

現(xiàn)在可檢查此數(shù)據(jù)的諧波成分，以確定電機的運行狀況。傅里葉分析是從波形中提取分量頻譜的數(shù)學過程。純正弦波的頻譜中僅包含一個頻率，稱為基波頻率。如果正弦波失真，將出現(xiàn)除基波頻率之外的其他頻率。通過分析電機振動模式的頻譜，可精確地診斷其運行狀況。

由于能夠執(zhí)行傅里葉分析的硬件和軟件通常價格很高，所以這里我們介紹一種可以對MEMS數(shù)據(jù)進行傅里葉分析的方法，基本上無需任何成本。

LTspice是一款功能強大、可免費使用的電路仿真器，它可以使用從狀態(tài)監(jiān)控系統(tǒng)的MEMS傳感器中獲取的振動數(shù)據(jù)，通過傅里葉分析繪制任何波形的頻譜。

通過數(shù)據(jù)格式，LTspice能夠生成傅里葉分析圖，其中每個振動數(shù)據(jù)點都與其相應的時間戳配對。

使用Microsoft Excel將數(shù)據(jù)轉(zhuǎn)換成這種格式相對比較容易。過程如下。

首先，將數(shù)據(jù)列分成Excel文件中的三個工作表，命名為X、Y和Z。

在數(shù)據(jù)左側(cè)插入一列——此列為每個數(shù)據(jù)值的時間戳。

由于在一秒內(nèi)提取了500,000個數(shù)據(jù)樣本，每個數(shù)據(jù)點間隔2 μs。因此，在新列的第一個單元格中，輸入2E-6

代表2μs處的第一個時間戳。

填充其余時間戳列數(shù)值的最簡單方式是使用Series命令。在Microsoft Excel的搜索框中，鍵入“Series”以顯示菜單選項。

從下拉菜單中選擇填充系列或模式(Fill Series or Pattern), 然后選擇系列…(Series…)。

此時出現(xiàn)對話框，選中列(Columns)和線性(Linear)單選按鈕。在步進值(Step value)中輸入2E-6，在停止值(Stop value)中輸入1。

點擊確定(OK)填充左列數(shù)據(jù)時間戳，從2 μs遞增到1秒。先填充前幾個值，然后將光標一直拖到數(shù)據(jù)范圍末尾的底部單元格，也可達到同樣的目的——但對于500,000行數(shù)據(jù)，需要拖得很長。

現(xiàn)在就得到LTspice可以處理的數(shù)據(jù)格式。

如果數(shù)據(jù)集很大，采樣間隔短，則Microsoft Excel可能會將時間戳四舍五入為不恰當?shù)男?shù)位數(shù)。如果出現(xiàn)這種情況，則突出顯示第一列，然后選擇然后選擇格式化(Format) > 格式化單元格(Format Cells)。

選擇合適的小數(shù)位數(shù)。

在填充時間戳列并擴展有效位數(shù)后，將每個工作表的兩列復制到記事本或其他文本編輯器文件中。

總共應該有三個文本文件，其中包含狀態(tài)監(jiān)控系統(tǒng)中X、Y和Z軸的振動數(shù)據(jù)。

現(xiàn)在，可將此數(shù)據(jù)直接讀入LTspice中。

在LTspice中構(gòu)建原理圖。在該設計中，有六個電壓源分別對應于故障和非故障的X，Y，Z軸的數(shù)據(jù)。這樣就可以對新電機的振動數(shù)據(jù)執(zhí)行傅里葉分析，并將分析結(jié)果與疑似故障電機數(shù)據(jù)的傅里葉分析進行比較。此方法的一大優(yōu)勢是新（非故障）電機的頻率圖可以疊加在疑似故障電機的頻率圖上，因此，性能差異一目了然。

.options plotwinsize=0 numdgt=15

去除了LTspice中的默認壓縮設置，有時會產(chǎn)生更清晰的結(jié)果。如果忽略此行，仿真運行速度會更快，但產(chǎn)生的結(jié)果可能不太精確。

完成原理圖后，右鍵單擊每個電壓源，選擇高級(Advanced)按鈕，選中PWL文件(PWL File)單選按鈕，然后輸入包含振動數(shù)據(jù)的相應文本文件的文件名，如圖12所示。這將創(chuàng)建一個分段線性電壓源，其中包含一系列電壓及其相應的時間實例。如果這些文本文件與LTspice文件存儲在同一目錄中，則操作會更簡單。

故障電機和非故障電機的仿真結(jié)果：該實驗在一臺轉(zhuǎn)速為587.3 rpm的電機上進行，電機的軸承出現(xiàn)故障，外環(huán)未對準，負載為12磅；同時還顯示了同一轉(zhuǎn)速下無故障電機的振動模式。顯然，與非故障電機相比，故障電機的振動特征幅度明顯更高。

突出顯示波形(Waveform)窗口，然后從菜單欄中選擇查看(View) > FFTT。這將基于瞬態(tài)數(shù)據(jù)計算FFT。

從數(shù)據(jù)可以看到，在35000V這樣如此高的失調(diào)電壓上，我們通過數(shù)字只能看到很小的變化。在LTspice中進行仿真時，這些數(shù)據(jù)會轉(zhuǎn)換成一個35,000 V的直流失調(diào)電壓，并在此失調(diào)電壓上還會疊加一個交流波形。

在傅里葉分析圖中，此失調(diào)電壓在頻譜位置的直流點上表現(xiàn)為很大的一個尖刺，因此，當LTspice自動縮放Y軸時，相關(guān)諧波比例極小。右鍵單擊X軸，指定高于直流電壓的頻率范圍，由此可忽略直流失調(diào)電壓——5 Hz至1 kHz應該足夠。

右鍵單擊Y軸，選擇線性(Linear)單選按鈕以查看諧波。

在圖形區(qū)單擊鼠標右鍵，可添加額外的繪圖窗格，即可將振動頻譜成分以X、Y和Z圖分別呈現(xiàn)。

可以清楚地看到電機的10 Hz旋轉(zhuǎn)頻率，以及60 Hz、142 Hz和172 Hz處存在明顯的諧波。雖然本文不會分析電機內(nèi)部的哪些組件導致了這些諧波，但毫無疑問，振動模式因電機磨損而改變。

ADI的MEMS加速度計系列能夠提供關(guān)鍵數(shù)據(jù)，進而在早期檢測出電機故障，但這只是解決方案的一半。必須通過傅里葉分析仔細研究這些數(shù)據(jù)。遺憾的是，能夠執(zhí)行傅里葉分析的設備或軟件通常很昂貴。而LTspice能夠免費精確分析CbM數(shù)據(jù)，從而實現(xiàn)早期檢測和診斷機器故障。