電機(jī)NVH仿真做過(guò)幾輪之后我最大的感受是這活兒不是單個(gè)軟件能搞定的也不是跑一次模態(tài)、出一張瀑布圖就叫“仿真完成”。真正的價(jià)值在于把結(jié)構(gòu)動(dòng)力學(xué)和電磁學(xué)串成一條完整的鏈路從“這個(gè)電機(jī)固有頻率在哪”一路追到“提速時(shí)哪個(gè)轉(zhuǎn)速會(huì)吵、會(huì)振”最后落在瀑布圖上讓問(wèn)題一眼可見。這篇文章就以“電機(jī)NVH仿真從模態(tài)分析到瀑布圖生成的探索”為主線把我實(shí)際走通的流程、每個(gè)環(huán)節(jié)的原理、以及那些文檔里不會(huì)寫的坑完整梳理一遍。適合正在做電機(jī)結(jié)構(gòu)設(shè)計(jì)、NVH仿真或者剛接觸振動(dòng)噪聲分析的朋友參考哪怕你是從零開始按這個(gè)思路也能搭建起自己的仿真路徑。1. 項(xiàng)目概述與整體思路拆解1.1 電機(jī)NVH到底在解決什么問(wèn)題電機(jī)的NVH本質(zhì)上是三個(gè)字母的事Noise噪聲、Vibration振動(dòng)、Harshness聲振粗糙度。但落到工程上真正讓人頭疼的是兩個(gè)東西一個(gè)是“電磁嘯叫”也就是大家常說(shuō)的高頻電嘯聲另一個(gè)是“結(jié)構(gòu)共振”某個(gè)轉(zhuǎn)速區(qū)間突然振得厲害、聲音刺耳。為什么電機(jī)特別容易出NVH問(wèn)題因?yàn)樗约壕褪且粋€(gè)激勵(lì)源。轉(zhuǎn)子轉(zhuǎn)動(dòng)時(shí)氣隙磁場(chǎng)會(huì)產(chǎn)生隨時(shí)間和空間變化的電磁力波這個(gè)力波作用在定子齒上把定子鐵芯、繞組、機(jī)殼整個(gè)結(jié)構(gòu)都“推”起來(lái)。電機(jī)結(jié)構(gòu)本身又是一個(gè)彈性系統(tǒng)有自己的固有頻率和振型。當(dāng)電磁力波的頻率和空間分布恰好跟結(jié)構(gòu)模態(tài)“對(duì)上”的時(shí)候就會(huì)出現(xiàn)共振振動(dòng)被放大噪聲一下子飆起來(lái)。所以電機(jī)NVH仿真說(shuō)白了就是回答三個(gè)問(wèn)題這個(gè)電機(jī)結(jié)構(gòu)有哪些固有頻率和振型電磁激勵(lì)產(chǎn)生的力波頻率和空間階次是多少兩者什么時(shí)候會(huì)碰上碰撞后振動(dòng)和噪聲有多大模態(tài)分析解決第一個(gè)問(wèn)題電磁力分析和響應(yīng)計(jì)算解決第二第三個(gè)問(wèn)題瀑布圖則是把“轉(zhuǎn)速-頻率-幅值”三者的關(guān)系可視化讓人一眼看出共振點(diǎn)在哪兒。1.2 從模態(tài)到瀑布圖的完整技術(shù)鏈路我做這個(gè)項(xiàng)目的整體鏈路是這樣的第一步結(jié)構(gòu)模態(tài)分析。把定子鐵芯、繞組、機(jī)殼、端蓋裝配體建模求解自由或約束狀態(tài)下的固有頻率和模態(tài)振型重點(diǎn)是關(guān)注低階模態(tài)和節(jié)徑數(shù)。第二步電磁力分析。用電磁場(chǎng)有限元跑一個(gè)轉(zhuǎn)速范圍或者幾個(gè)典型工況提取作用在定子齒表面的徑向電磁力波再對(duì)力波做空間和時(shí)間上的傅里葉分解得到力波的“階次-頻率-幅值”譜。第三步共振風(fēng)險(xiǎn)預(yù)判。把模態(tài)頻率和電磁力頻率放在一起對(duì)比疊成Campbell圖看有沒有交點(diǎn)、交點(diǎn)落在哪個(gè)轉(zhuǎn)速區(qū)間。第四步振動(dòng)響應(yīng)計(jì)算。把電磁力作為激勵(lì)加載到結(jié)構(gòu)模型上做諧響應(yīng)分析得到殼體表面在不同轉(zhuǎn)速下的振動(dòng)速度或加速度。第五步瀑布圖生成。把每個(gè)轉(zhuǎn)速下響應(yīng)信號(hào)的頻譜做出來(lái)按轉(zhuǎn)速堆疊成三維譜圖。從瀑布圖上可以清晰地看到階次線、共振帶以及它們交匯形成的“熱點(diǎn)”。這套鏈路的難度不在單點(diǎn)而在銜接。模態(tài)和電磁分別都是成熟學(xué)科但中間的數(shù)據(jù)映射、坐標(biāo)系對(duì)齊、單位統(tǒng)一、邊界條件匹配才是真正決定結(jié)果準(zhǔn)不準(zhǔn)的地方。1.3 為什么這條鏈路值得完整走一遍很多人做NVH仿真喜歡“缺哪補(bǔ)哪”共振了就跑個(gè)模態(tài)噪聲大了就算個(gè)電磁力。但這樣往往解決不了根本問(wèn)題。模態(tài)頻率差了你沒算電磁力就不知道是激勵(lì)問(wèn)題還是結(jié)構(gòu)問(wèn)題電磁力算準(zhǔn)了模態(tài)模型不對(duì)響應(yīng)結(jié)果照樣偏得離譜。把整條鏈路走完最直接的好處是能定位問(wèn)題源頭。舉個(gè)例子瀑布圖上如果某個(gè)轉(zhuǎn)速下出現(xiàn)了明顯的斜線加亮帶斜線代表階次激勵(lì)亮帶代表某階模態(tài)被激起。你順著這條斜線找到對(duì)應(yīng)的電磁力階次再對(duì)比亮帶頻率對(duì)應(yīng)哪階模態(tài)就能馬上判斷是應(yīng)該改極槽配合來(lái)削弱激勵(lì)還是應(yīng)該加筋加厚來(lái)把模態(tài)頻率移開。這種“對(duì)癥下藥”的效率遠(yuǎn)高于盲試。另外完整的NVH仿真模型也是后續(xù)做多目標(biāo)優(yōu)化的基礎(chǔ)?，F(xiàn)在電機(jī)設(shè)計(jì)越來(lái)越卷效率、轉(zhuǎn)矩、成本、NVH要一起權(quán)衡。沒有一條可信的NVH快速仿真鏈路優(yōu)化就無(wú)從談起。2. 模態(tài)分析邊界條件決定一切2.1 模態(tài)分析基礎(chǔ)固有頻率不是“玄學(xué)”模態(tài)分析的物理本質(zhì)是求解結(jié)構(gòu)在無(wú)外力作用下的自由振動(dòng)方程[ [M]{\ddot{x}} [K]{x} 0 ]求解得到的特征值對(duì)應(yīng)系統(tǒng)的固有頻率特征向量對(duì)應(yīng)振型。換句話說(shuō)固有頻率取決于結(jié)構(gòu)的質(zhì)量矩陣[M]和剛度矩陣[K]也就是材料屬性、幾何形狀和約束方式共同決定的結(jié)果。對(duì)于電機(jī)這種旋轉(zhuǎn)電機(jī)結(jié)構(gòu)最需要關(guān)注的是定子鐵芯和機(jī)殼的周向模態(tài)。因?yàn)殡姶帕Σㄖ饕趶较蜃饔米钊菀准て鸬囊彩嵌ㄗ拥膹较蚰B(tài)。這類模態(tài)可以用節(jié)徑數(shù)來(lái)描述0節(jié)徑是“呼吸模態(tài)”整個(gè)圓環(huán)均勻脹縮1節(jié)徑是“擺動(dòng)模態(tài)”整體像剛體一樣平移2節(jié)徑是“橢圓模態(tài)”圓變成橢圓3節(jié)徑是“三角形模態(tài)”以此類推。2節(jié)徑模態(tài)和0節(jié)徑模態(tài)在電機(jī)里最常被激發(fā)。定子鐵芯的2節(jié)徑固有頻率通常落在幾百赫茲到幾千赫茲的區(qū)間正好和電磁力波的頻率范圍重合所以共振風(fēng)險(xiǎn)很高。2.2 電機(jī)定子核心模態(tài)與節(jié)徑概念我們把定子鐵芯簡(jiǎn)化為一個(gè)圓環(huán)來(lái)看。圓環(huán)的周向模態(tài)和節(jié)徑數(shù)n的關(guān)系很直接n0時(shí)整個(gè)環(huán)在做均勻的徑向呼吸n2時(shí)一圈上有兩個(gè)對(duì)稱的“波峰波谷”形變是橢圓n3時(shí)是三角形形變。節(jié)徑越多模態(tài)頻率通常越高因?yàn)楦呖臻g階次的形變需要更多彎曲能量。但電機(jī)定子不是裸圓環(huán)上面還有齒槽、繞組、機(jī)殼。齒槽的存在會(huì)顯著增加局部剛度繞組則主要增加質(zhì)量機(jī)殼對(duì)整體剛度的貢獻(xiàn)則取決于材料和連接方式。所有這些都會(huì)讓模態(tài)頻率偏移所以仿真建模時(shí)必須把這些細(xì)節(jié)考慮進(jìn)去否則算出來(lái)的頻率可能偏差很大。這里我踩過(guò)一個(gè)坑第一次建模時(shí)把繞組簡(jiǎn)化成純附加質(zhì)量忽略了銅線、絕緣和浸漆形成的剛度貢獻(xiàn)結(jié)果2節(jié)徑模態(tài)算出來(lái)比測(cè)試值低了差不多10%。后來(lái)把繞組用等效各向異性材料處理——軸向和徑向剛度不同密度按填充率折算——才把誤差壓到3%以內(nèi)。所以繞組不是簡(jiǎn)單的“增重”它對(duì)定子剛度的影響在某些模態(tài)下非常顯著。2.3 邊界條件設(shè)置自由模態(tài)vs約束模態(tài)模態(tài)分析的邊界條件直接決定了結(jié)果的含義。做電機(jī)NVH復(fù)現(xiàn)共振問(wèn)題時(shí)我通常第一步先做自由-自由模態(tài)分析也就是說(shuō)機(jī)殼和端蓋不施加任何約束只模擬結(jié)構(gòu)“懸空”狀態(tài)下的固有特性。有人會(huì)問(wèn)電機(jī)明明安裝在車架上或有掛腳約束為什么要算自由模態(tài)原因有兩個(gè)。第一自由模態(tài)是結(jié)構(gòu)自身的固有屬性方便和實(shí)驗(yàn)對(duì)標(biāo)——實(shí)驗(yàn)室里用軟繩懸掛電機(jī)做錘擊法測(cè)模態(tài)測(cè)出來(lái)的就是自由模態(tài)。第二實(shí)際安裝約束非常復(fù)雜螺栓扭矩、接觸面、安裝支架剛度都難以精確模擬如果一開始就把這些不確定性引進(jìn)來(lái)反而說(shuō)不清頻率偏差到底來(lái)自哪。等自由模態(tài)對(duì)標(biāo)完了再根據(jù)實(shí)際安裝方式補(bǔ)充約束模態(tài)或裝配體模態(tài)。這里的原則是從簡(jiǎn)單到復(fù)雜逐步加約束每一步都確認(rèn)模型行為合理再進(jìn)入下一步。直接一上來(lái)就做全約束裝配體分析一旦結(jié)果不對(duì)排查成本會(huì)非常高。2.4 模態(tài)實(shí)驗(yàn)對(duì)標(biāo)與相關(guān)性仿真模型不是“算完就完了”必須跟實(shí)驗(yàn)對(duì)標(biāo)。模態(tài)實(shí)驗(yàn)通常用錘擊法或激振器法在定子或機(jī)殼表面布置加速度傳感器通過(guò)力錘激勵(lì)測(cè)得頻響函數(shù)FRF再用模態(tài)擬合軟件提取模態(tài)頻率、阻尼和振型。對(duì)標(biāo)時(shí)不僅要看頻率數(shù)值還要看振型形狀。這里用到一個(gè)很常用的指標(biāo)——模態(tài)置信準(zhǔn)則MAC值用來(lái)定量評(píng)估兩個(gè)模態(tài)振型的相關(guān)性。MAC值在0到1之間1表示完全相關(guān)。實(shí)際操作中仿真和測(cè)試的同階模態(tài)MAC值最好在0.8以上低于0.6基本可以斷定振型對(duì)不上需要回頭檢查建?；蜻吔鐥l件。我做過(guò)一次對(duì)標(biāo)頻率誤差只有2%但4階振型的MAC值只有0.5。后來(lái)發(fā)現(xiàn)是機(jī)殼和定子鐵芯之間的過(guò)盈配合在模型里被簡(jiǎn)化成了“綁定”實(shí)際接觸剛度不足導(dǎo)致高階振型形狀明顯不同。改成帶摩擦的接觸之后MAC值立刻升到了0.9以上。這就是我反復(fù)強(qiáng)調(diào)的頻率對(duì)得上不代表模態(tài)模型就對(duì)了振型相關(guān)性同樣關(guān)鍵。3. 電磁激勵(lì)與共振匹配NVH的“因”3.1 電磁力波的來(lái)源麥克斯韋應(yīng)力與徑向力波電機(jī)NVH的激勵(lì)源本質(zhì)上是氣隙磁場(chǎng)產(chǎn)生的麥克斯韋應(yīng)力。氣隙磁場(chǎng)的徑向分量(B_r)在定子齒表面上產(chǎn)生的徑向應(yīng)力可以表示為[ \sigma_r \frac{B_r^2}{2\mu_0} ]這個(gè)式子看著簡(jiǎn)單難點(diǎn)在于(B_r)本身是時(shí)間和空間的復(fù)雜函數(shù)。轉(zhuǎn)子旋轉(zhuǎn)時(shí)氣隙磁場(chǎng)中存在大量空間諧波和時(shí)間諧波不同諧波相互作用會(huì)產(chǎn)生一系列頻率、階次各不相同的徑向力波。力波有兩個(gè)關(guān)鍵屬性空間階次沿圓周分布幾個(gè)波峰波谷和時(shí)間頻率每秒波動(dòng)多少次?？臻g階次決定了它容易激起哪一階結(jié)構(gòu)模態(tài)時(shí)間頻率決定了它是否落在結(jié)構(gòu)固有頻率附近。為什么我們要特別關(guān)心低階力波因?yàn)榻Y(jié)構(gòu)的模態(tài)剛度通常隨空間階次增加而變大低階模態(tài)更容易被激勵(lì)出大變形。所以電磁設(shè)計(jì)時(shí)大家普遍追求“最小的非零力波階次盡量高”這樣才能避開低階結(jié)構(gòu)模態(tài)。3.2 力波階次計(jì)算極槽配合的底層邏輯力波的空間階次跟極數(shù)和槽數(shù)直接相關(guān)。對(duì)于永磁同步電機(jī)常見的徑向力波階次可以用這個(gè)式子描述[ r kN_s \pm 2p \quad (k 0, 1, 2, \ldots) ]其中(N_s)是定子槽數(shù)(p)是極對(duì)數(shù)(2p)是極數(shù)。實(shí)際中最小的非零力波階次通常就是(\min(kN_s \pm 2p))。舉個(gè)例子一臺(tái)8極48槽整數(shù)槽電機(jī)(2p8)(N_s48)。當(dāng)(k0)時(shí)(r8)所以最小力波階次是8階。8階力波對(duì)應(yīng)的空間波長(zhǎng)很短不容易激起低階結(jié)構(gòu)模態(tài)這也是整數(shù)槽繞組電機(jī)通常NVH表現(xiàn)較好的原因之一。再看一臺(tái)什么情況會(huì)出問(wèn)題。有些分?jǐn)?shù)槽集中繞組電機(jī)比如8極9槽雖然繞組系數(shù)和槽利用率有優(yōu)勢(shì)但力波階次會(huì)低得多。用同樣的公式(N_s9)(2p8)取(k1)時(shí)(r1)這意味著存在1階力波也就是在圓周上只有一個(gè)波峰的分布。這種低階力波很容易激發(fā)定子的1節(jié)徑模態(tài)中低頻電磁噪聲和振動(dòng)就很突出。這就是為什么分?jǐn)?shù)槽電機(jī)的NVH設(shè)計(jì)難度更大的底層原因。3.3 Campbell圖與共振頻率對(duì)齊判斷有了模態(tài)頻率也有了電磁力的頻率和階次下一步就是把它們放到一張圖上對(duì)照。Campbell圖就是干這個(gè)用的。橫軸是轉(zhuǎn)速縱軸是頻率圖上畫若干條從原點(diǎn)出發(fā)的斜線每條代表電磁激勵(lì)的某個(gè)階次激勵(lì)頻率(f \text{order} \times \text{轉(zhuǎn)速RPM} / 60)。同時(shí)把結(jié)構(gòu)模態(tài)頻率畫成水平線或近似水平線斜線穿過(guò)水平線的地方就是潛在的共振點(diǎn)。但不要見到交叉就說(shuō)“這里會(huì)共振”。還要滿足兩個(gè)條件第一頻率重合。這個(gè)比較容易理解激勵(lì)頻率等于或接近固有頻率頻率差最好在5%以內(nèi)否則共振放大效應(yīng)有限。第二空間階次匹配。電磁力波的空間階次必須和模態(tài)的節(jié)徑數(shù)一致或接近才能有效輸入能量。比如2階力波激發(fā)2節(jié)徑模態(tài)效率最高但如果模態(tài)是4節(jié)徑2階力波就很難把它激勵(lì)起來(lái)。階次差得越大耦合效率衰減越快。所以Campbell圖上即使看到很多交叉點(diǎn)真正需要關(guān)注的往往是那些“階次與節(jié)徑匹配”的交叉點(diǎn)。這也是為什么純用“頻率是否重合”來(lái)判斷共振經(jīng)常會(huì)判斷失誤的原因。3.4 電磁力提取與FFT處理從電磁場(chǎng)有限元里提取力波通常會(huì)遍歷一個(gè)電周期內(nèi)的多個(gè)時(shí)間步得到定子齒表面各節(jié)點(diǎn)的徑向電磁力隨時(shí)間的變化。然后對(duì)一個(gè)空間周期做空間FFT得到各階次力波的幅值對(duì)時(shí)間序列做時(shí)間FFT得到各頻率成分的幅值。這里有個(gè)實(shí)用的建議提取力的時(shí)候不要把每個(gè)節(jié)點(diǎn)的力都導(dǎo)出來(lái)那樣數(shù)據(jù)量巨大且沒意義。更好的做法是把每個(gè)定子齒頂面的力按面積積分得到“每齒的合力”再沿圓周展開成空間序列做FFT。這樣既保留了力波的空間分布特征又大幅壓縮了數(shù)據(jù)量。時(shí)間FFT時(shí)要注意采樣點(diǎn)數(shù)和采樣時(shí)長(zhǎng)。電磁力信號(hào)里既包含恒定的直流分量也包含脈動(dòng)分量。直流分量不產(chǎn)生振動(dòng)處理時(shí)可以直接濾掉脈動(dòng)的關(guān)鍵頻率通常是電頻率的2倍及其整數(shù)倍因?yàn)榱φ扔诖琶芷椒酱琶苊孔兓粋€(gè)周期力會(huì)變化兩個(gè)周期。4. 振動(dòng)響應(yīng)計(jì)算與瀑布圖生成的完整實(shí)操4.1 從電磁力到振動(dòng)響應(yīng)的傳遞路徑模態(tài)分析告訴我們結(jié)構(gòu)“有哪些頻率容易振”電磁力分析告訴我們“有什么頻率在激勵(lì)”而振動(dòng)響應(yīng)計(jì)算就是把兩者真正“接”起來(lái)算出實(shí)際振動(dòng)有多大。在頻域里做這事最直觀的方法是諧響應(yīng)分析。把每階電磁力波作為簡(jiǎn)諧激勵(lì)頻率設(shè)為力波的時(shí)間頻率空間分布按力波階次映射到定子齒表面然后求解結(jié)構(gòu)在這個(gè)頻率下的穩(wěn)態(tài)振動(dòng)響應(yīng)。把所有關(guān)注的頻率都算一遍就得到結(jié)構(gòu)表面在不同頻率下的振動(dòng)速度或加速度。這個(gè)過(guò)程的精度很大程度上取決于電磁力和結(jié)構(gòu)模型之間的數(shù)據(jù)傳遞。電磁分析里力作用在定子齒的齒頂面上但結(jié)構(gòu)網(wǎng)格和電磁網(wǎng)格通常不是一個(gè)體系需要做插值映射。映射時(shí)要注意保護(hù)總力和合力矩不然局部力的偏差會(huì)導(dǎo)致響應(yīng)結(jié)果失真。響應(yīng)輸出位置也要多想一步。定子鐵芯本身振動(dòng)但客戶關(guān)心的通常是機(jī)殼外表面特別是安裝腳、吊耳這些直接傳遞振動(dòng)的路徑。所以模型里至少要保留機(jī)殼和端蓋不能只單獨(dú)分析定子鐵芯。4.2 瀑布圖的數(shù)據(jù)來(lái)源時(shí)變轉(zhuǎn)速下的頻譜切片瀑布圖的橫軸是頻率縱軸是轉(zhuǎn)速或時(shí)間顏色表示振動(dòng)幅值或聲壓級(jí)。每一“行”數(shù)據(jù)本質(zhì)上是一個(gè)轉(zhuǎn)速下振動(dòng)信號(hào)的頻譜。所以做瀑布圖的第一步就是要有“轉(zhuǎn)速?gòu)男〉酱髵哌^(guò)去”的振動(dòng)響應(yīng)數(shù)據(jù)。這個(gè)數(shù)據(jù)來(lái)源可以是仿真也可以是試驗(yàn)。如果是仿真做法是在一個(gè)轉(zhuǎn)速范圍內(nèi)選若干轉(zhuǎn)速點(diǎn)比如500到6000轉(zhuǎn)/分每隔100轉(zhuǎn)取一個(gè)點(diǎn)。每個(gè)點(diǎn)都做一次諧響應(yīng)分析得到一個(gè)頻譜把所有點(diǎn)的頻譜按轉(zhuǎn)速堆疊就得到了瀑布圖。如果是試驗(yàn)則通常用加速度傳感器測(cè)殼體振動(dòng)同時(shí)用電控信號(hào)記錄轉(zhuǎn)速脈沖在一個(gè)勻加速或勻減速工況下連續(xù)采集數(shù)據(jù)然后按時(shí)間窗切片做FFT。每個(gè)時(shí)間窗對(duì)應(yīng)一個(gè)轉(zhuǎn)速把每個(gè)窗口的頻譜按轉(zhuǎn)速排列也能得到瀑布圖。這里有個(gè)細(xì)節(jié)值得注意瀑布圖里每個(gè)頻譜切片的質(zhì)量取決于窗函數(shù)的設(shè)置。如果窗函數(shù)太短頻率分辨率太低如果太長(zhǎng)窗口內(nèi)轉(zhuǎn)速變化太大頻譜就會(huì)被“抹糊”。通常我建議窗口內(nèi)的轉(zhuǎn)速變化控制在50轉(zhuǎn)/分以內(nèi)譜線數(shù)在800到1600線之間兼顧頻率分辨率和時(shí)域分辨率。下面是仿真生成瀑布圖的簡(jiǎn)化流程以Python為例供參考import numpy as np from scipy.signal import spectrogram speeds np.arange(500, 6001, 100) # 轉(zhuǎn)速序列500~6000 rpm步長(zhǎng)100 freq_axis np.fft.rfftfreq(nfft, d1/fs) # 頻率軸 waterfall [] # 每個(gè)轉(zhuǎn)速存一份頻譜 for rpm in speeds: # 假設(shè) response_spectrum(rpm) 返回該轉(zhuǎn)速下振動(dòng)信號(hào)的幅值譜 spec response_spectrum(rpm, freq_axis) waterfall.append(spec) waterfall_db 20 * np.log10(np.array(waterfall) / 1e-6 1e-12) # 轉(zhuǎn)為dB # 最后用 matplotlib 的 pcolormesh / pcolor 畫二維色圖即可當(dāng)然實(shí)際工程里我們很少自己在Python里造輪子生成瀑布圖一般直接用LMS Test.Lab、Simcenter、或者M(jìn)ATLAB的NVH工具箱直接出圖。但理解這個(gè)數(shù)據(jù)處理的邏輯很重要否則你會(huì)搞不明白為什么同一條瀑布圖不同人設(shè)置不同結(jié)果差異竟然這么大。4.3 瀑布圖怎么看階次線與共振帶識(shí)別瀑布圖出來(lái)以后怎么讀才是關(guān)鍵。我一般按三步來(lái)看。第一步找階次線。測(cè)得的振動(dòng)信號(hào)中那些“從原點(diǎn)出發(fā)、隨著轉(zhuǎn)速升高而頻率線性升高”的亮線就是階次激勵(lì)線。階次1對(duì)應(yīng)轉(zhuǎn)子每轉(zhuǎn)一圈振動(dòng)一次階次2就是轉(zhuǎn)一圈振動(dòng)兩次。電磁激勵(lì)的話永磁同步電機(jī)常見的電磁階次跟極數(shù)和供電頻率有關(guān)比如在瀑布圖上會(huì)出現(xiàn)2f電頻率相關(guān)的斜線換算到機(jī)械階次往往就是極數(shù)或極數(shù)倍數(shù)。第二步找共振帶。如果某階結(jié)構(gòu)模態(tài)被激勵(lì)就會(huì)在某個(gè)固定頻率位置出現(xiàn)一條“豎帶”或“水平亮帶”——不管轉(zhuǎn)速怎么變這條帶的頻率基本不動(dòng)因?yàn)楣逃蓄l率不隨轉(zhuǎn)速變化。它的亮度會(huì)沿著轉(zhuǎn)速方向變化轉(zhuǎn)速接近共振轉(zhuǎn)速時(shí)最亮。第三步找“熱點(diǎn)”。階次線穿過(guò)共振帶的區(qū)域就是共振最劇烈的轉(zhuǎn)速-頻率點(diǎn)。在這個(gè)點(diǎn)激勵(lì)頻率和固有頻率重合且激勵(lì)階次能有效激發(fā)對(duì)應(yīng)模態(tài)振動(dòng)幅值會(huì)出現(xiàn)一個(gè)明顯峰值。瀑布圖上的熱點(diǎn)往往對(duì)應(yīng)實(shí)際駕駛中“某個(gè)轉(zhuǎn)速下突然聽到嘯叫”的現(xiàn)象。我見過(guò)很多剛開始接觸NVH的人拿到瀑布圖就找最亮的點(diǎn)把所有亮點(diǎn)都當(dāng)成共振。實(shí)際上有些亮線只是正常電磁階次激勵(lì)幅值雖高但遠(yuǎn)低于結(jié)構(gòu)模態(tài)頻率并不會(huì)產(chǎn)生劇烈共振。所以看瀑布圖時(shí)必須結(jié)合模態(tài)頻率和階次信息一起判斷不能“看圖說(shuō)話”。4.4 從瀑布圖反推優(yōu)化方向瀑布圖不止是用來(lái)“看問(wèn)題”它還能直接告訴你“怎么改”。如果熱點(diǎn)出現(xiàn)在某個(gè)階次線和某條共振帶的交匯處你可以判斷這條共振帶對(duì)應(yīng)哪階模態(tài)這個(gè)階次線來(lái)自哪個(gè)電磁力諧波然后針對(duì)性地處理。舉個(gè)例子我遇到過(guò)一次8極48槽電機(jī)在3000轉(zhuǎn)/分左右出現(xiàn)明顯嘯叫。瀑布圖上看熱點(diǎn)頻率在2500Hz左右對(duì)應(yīng)的是12節(jié)徑模態(tài)階次線是一條斜率較高的線計(jì)算下來(lái)對(duì)應(yīng)電磁力的某個(gè)高次諧波。后來(lái)把注意力放在那對(duì)極槽配合的諧波優(yōu)化上通過(guò)調(diào)整磁鋼形狀減小這個(gè)諧波幅值熱點(diǎn)顯著減弱嘯叫消失了。如果問(wèn)題是模態(tài)頻率落在常用轉(zhuǎn)速區(qū)間內(nèi)那就要改結(jié)構(gòu)。比如在機(jī)殼上加軸向筋條、增加壁厚、調(diào)整端蓋連接剛度把模態(tài)頻率移到激勵(lì)較弱的頻段。這比削弱電磁激勵(lì)更直接但會(huì)增加重量和成本所以通常放在電磁優(yōu)化之后考慮。5. 常見問(wèn)題與排查技巧實(shí)錄5.1 模態(tài)仿真結(jié)果“偏剛”或“偏柔”怎么調(diào)這是模態(tài)分析里最常見的困惑。頻率算出來(lái)比測(cè)試高說(shuō)明模型偏剛反之則是偏柔。按我的經(jīng)驗(yàn)問(wèn)題通常出在幾個(gè)地方。第一個(gè)是繞組建模。繞組對(duì)整個(gè)定子模態(tài)頻率的影響非常大。建模太簡(jiǎn)單比如只當(dāng)質(zhì)量塊頻率會(huì)偏低如果浸漆、銅線、絕緣紙全都考慮進(jìn)去等效剛度又會(huì)偏高因?yàn)閷?shí)際繞組各匝之間并非完全剛性約束。建議用正交各向異性材料來(lái)等效徑向彈性模量取一個(gè)較低值周向取較高值密度按實(shí)測(cè)填充率折算。第二個(gè)是接觸連接。鐵芯疊片之間、鐵芯和機(jī)殼之間、機(jī)殼和端蓋之間都存在接觸剛度。如果全都設(shè)成綁定的剛性連接頻率會(huì)偏高如果設(shè)成無(wú)摩擦接觸又可能偏低。通常的做法是先綁定額定主要件再在機(jī)殼-鐵芯接觸部位用帶摩擦接觸并施加過(guò)盈量用模態(tài)測(cè)試數(shù)據(jù)來(lái)標(biāo)定接觸參數(shù)。第三個(gè)是阻尼的影響。模態(tài)頻率本身對(duì)阻尼不敏感但阻尼對(duì)共振幅值影響巨大。如果諧響應(yīng)分析的結(jié)果跟實(shí)測(cè)幅值差很多先檢查阻尼比設(shè)置結(jié)構(gòu)阻尼通常取0.5%~2%不要默認(rèn)取0。5.2 瀑布圖“糊成一團(tuán)”的原因瀑布圖做出來(lái)如果模糊、頻率軸糊成一片最常見的原因是窗函數(shù)和譜線數(shù)不匹配。頻率分辨率(\Delta f)和FFT窗口時(shí)長(zhǎng)(T)的關(guān)系是(\Delta f 1/T)。如果你用1秒的窗口頻率分辨率是1Hz想看5Hz以內(nèi)的模態(tài)細(xì)節(jié)就看不清楚。如果窗口太短分辨率不足很多相鄰頻率就會(huì)混在一起。但窗口拉長(zhǎng)又有風(fēng)險(xiǎn)尤其是在掃頻工況下窗口內(nèi)轉(zhuǎn)速變化太大相當(dāng)于把好幾個(gè)轉(zhuǎn)速狀態(tài)混在一起平均頻譜特征被模糊。解決辦法是提高采樣率、增大FFT點(diǎn)數(shù)、使用漢寧窗降低頻譜泄漏同時(shí)控制窗口內(nèi)轉(zhuǎn)速漂移不超過(guò)總測(cè)試范圍的2%~3%。另外顏色標(biāo)尺的設(shè)置也會(huì)影響判斷。默認(rèn)的顏色映射可能把小幅值差異掩蓋了建議手動(dòng)調(diào)整色標(biāo)下限把本底噪聲壓下去讓階次線和共振帶更容易辨認(rèn)。5.3 計(jì)算資源與收斂性問(wèn)題完整跑一輪“電磁力諧響應(yīng)瀑布圖”計(jì)算量不小。電磁分析幾分鐘到幾十分鐘不等諧響應(yīng)分析如果掃頻點(diǎn)很多每個(gè)工作頻率都要解一次大規(guī)模方程時(shí)間會(huì)成倍增加。這里分享一個(gè)省時(shí)間的技巧不需要每個(gè)轉(zhuǎn)速點(diǎn)都重算電磁力。如果電磁力只跟負(fù)載和電流角有關(guān)而對(duì)轉(zhuǎn)速本身不敏感忽略轉(zhuǎn)速引起的渦流和磁飽和變化那可以用幾個(gè)典型負(fù)載點(diǎn)算電磁力再換算到不同轉(zhuǎn)速下的頻率直接復(fù)用這些力數(shù)據(jù)。實(shí)測(cè)下來(lái)這個(gè)方法在工程精度范圍內(nèi)完全可行能省掉大半計(jì)算時(shí)間。諧響應(yīng)計(jì)算也可以用模態(tài)疊加法而不是直接法。模態(tài)疊加法只保留最重要的幾十階模態(tài)計(jì)算速度比直接法快一個(gè)數(shù)量級(jí)。前提是你已經(jīng)做過(guò)模態(tài)分析并且確認(rèn)關(guān)注的頻率范圍內(nèi)模態(tài)截?cái)嗾`差可控。通常保留0到最高關(guān)注頻率2倍范圍內(nèi)的模態(tài)就夠用了。5.4 仿真與測(cè)試對(duì)不齊時(shí)的排查思路對(duì)不齊的時(shí)候我的第一反應(yīng)不是懷疑求解器而是回頭檢查輸入條件。整理的排查順序是這樣的先查頻率單位。Hz和rad/s搞混的問(wèn)題我見過(guò)很多次一下就能差6.28倍。再查材料參數(shù)。銅繞組密度、硅鋼片彈性模量、機(jī)殼鋁合金的阻尼比這些參數(shù)來(lái)源不同廠家可能不一樣必須在項(xiàng)目一開始就確認(rèn)清楚不要隨便從材料庫(kù)拉一個(gè)“默認(rèn)鋁”就用。接著查邊界條件。實(shí)驗(yàn)里電機(jī)是懸置的還是剛性安裝的螺栓扭矩多大有沒有加墊片這些細(xì)節(jié)對(duì)模態(tài)頻率的影響可能遠(yuǎn)大于材料參數(shù)。最后查激勵(lì)加載。電磁力的幅值、相位、作用位置是否正確映射到了結(jié)構(gòu)網(wǎng)格上。尤其是力作用在齒頂面和齒根面的差別對(duì)響應(yīng)結(jié)果影響很大卻容易被忽略。6. 后續(xù)擴(kuò)展與個(gè)人體會(huì)6.1 把NVH仿真流程做成自動(dòng)化流水線流程跑通之后最有價(jià)值的改進(jìn)就是自動(dòng)化。模態(tài)分析、電磁力提取、諧響應(yīng)計(jì)算、瀑布圖生成每一步都有固定的輸入和輸出非常適合用腳本串起來(lái)。我現(xiàn)在會(huì)在電磁仿真軟件里把力數(shù)據(jù)按標(biāo)準(zhǔn)格式導(dǎo)出來(lái)然后用批處理腳本驅(qū)動(dòng)結(jié)構(gòu)求解器自動(dòng)完成諧響應(yīng)分析最后把結(jié)果統(tǒng)一回收到數(shù)據(jù)后處理腳本里集中出圖。一個(gè)轉(zhuǎn)速范圍幾十個(gè)點(diǎn)全自動(dòng)跑完只需要一個(gè)晚上白天就能拿到瀑布圖來(lái)分析。這在項(xiàng)目早期做多方案對(duì)比時(shí)太有用了能讓你在兩天內(nèi)篩掉大批不合理的結(jié)構(gòu)方案。要自動(dòng)化一定要注重命名規(guī)范和數(shù)據(jù)格式的統(tǒng)一。每一個(gè)仿真實(shí)例都要能追溯到它的電磁模型版本、結(jié)構(gòu)模型版本和求解設(shè)置。不然跑了十個(gè)方案最后都分不清哪張圖是哪個(gè)方案出來(lái)的自動(dòng)化就成了災(zāi)難。6.2 模態(tài)分析在更廣領(lǐng)域的遷移價(jià)值做電機(jī)NVH繞不開模態(tài)分析但模態(tài)分析的應(yīng)用范圍遠(yuǎn)超電機(jī)。龍門機(jī)床的模態(tài)分析關(guān)注的是整機(jī)在切削力作用下的動(dòng)態(tài)剛度通過(guò)錘擊或激振得到各階模態(tài)判斷切削顫振風(fēng)險(xiǎn)方法和電機(jī)定子模態(tài)分析一脈相承。結(jié)構(gòu)動(dòng)力學(xué)的基本邏輯是通用的換的是激勵(lì)源、邊界條件和關(guān)注頻段。其實(shí)“多模態(tài)”的思維也不只屬于模態(tài)分析。比如在信號(hào)處理里時(shí)頻分析同樣強(qiáng)調(diào)從多個(gè)維度看數(shù)據(jù)在不同領(lǐng)域間遷移方法論往往能碰撞出新解法。這也是我鼓勵(lì)剛?cè)胄械娜讼仍鷮?shí)做通一個(gè)完整流程的原因——流程背后的底層思維才是可以帶去所有后續(xù)項(xiàng)目里的資產(chǎn)。6.3 幾點(diǎn)個(gè)人經(jīng)驗(yàn)供參考最后說(shuō)幾句實(shí)在話。第一NVH仿真永遠(yuǎn)不要追求“一次算準(zhǔn)”工程上更可行的目標(biāo)是“趨勢(shì)正確、誤差可控”。頻率誤差3%以內(nèi)、熱點(diǎn)位置能對(duì)上、階次激勵(lì)能分清就已經(jīng)有很強(qiáng)的工程決策價(jià)值了。第二仿真的盡頭是實(shí)驗(yàn)。不管模型做得多精細(xì)最后都要用測(cè)試來(lái)驗(yàn)收。所以我始終建議在做仿真的同時(shí)提前想好后續(xù)實(shí)驗(yàn)怎么做測(cè)點(diǎn)放哪里、用什么激勵(lì)方式、關(guān)注什么頻段這些實(shí)驗(yàn)方案如果能跟仿真模型的設(shè)計(jì)一起考慮整個(gè)項(xiàng)目的效率會(huì)高非常多。第三多留一點(diǎn)時(shí)間給數(shù)據(jù)整理和文檔記錄。NVH項(xiàng)目的數(shù)據(jù)量非常大而且涉及多學(xué)科如果沒有良好的記錄習(xí)慣三個(gè)月后回頭查問(wèn)題光回憶就夠你折騰好幾天。