周長(zhǎng)峰,張家豪,劉立邦,劉新磊,王希波
(山東交通學(xué)院 汽車工程學(xué)院,山東 濟(jì)南 250357)
作為汽車制動(dòng)系統(tǒng)的核心部件,制動(dòng)盤的振動(dòng)、噪聲特性對(duì)整車的噪聲、振動(dòng)和平順(noise vibration harshness,NVH)性能有重大影響。通過模態(tài)分析的手段,對(duì)制動(dòng)盤進(jìn)行模態(tài)特性研究,可以為解決車輛制動(dòng)時(shí)的振動(dòng)和噪聲問題奠定基礎(chǔ)[1]。
陳曉鵬[2]通過改變測(cè)點(diǎn)數(shù)量對(duì)制動(dòng)盤模態(tài)振型的差別進(jìn)行了描述;褚志剛等[3]和陳國(guó)強(qiáng)等[4]通過試驗(yàn)與仿真結(jié)果的對(duì)比,驗(yàn)證了有限元模型的準(zhǔn)確性;文獻(xiàn)[5-7]通過對(duì)制動(dòng)盤的試驗(yàn)?zāi)B(tài)分析,得到了模態(tài)振型與噪聲的關(guān)系;文獻(xiàn)[8-10]通過仿真求解制動(dòng)盤模態(tài),并對(duì)其進(jìn)行結(jié)構(gòu)優(yōu)化;焦標(biāo)強(qiáng)等[11]和譚祥軍等[12]通過試驗(yàn)的手段對(duì)制動(dòng)盤進(jìn)行模態(tài)分析,研究了參數(shù)辨識(shí)過程中的影響因素。目前的研究大多是單純的仿真或是通過仿真與試驗(yàn)結(jié)果對(duì)比來驗(yàn)證有限元模型的準(zhǔn)確性,對(duì)試驗(yàn)結(jié)果的準(zhǔn)確性判斷并未過多研究。
本文對(duì)某汽車制動(dòng)盤進(jìn)行計(jì)算模態(tài)分析,得出制動(dòng)盤前五階模態(tài)特性,并通過頻率響應(yīng)函數(shù)(frame repetition frequency,F(xiàn)RF)和模態(tài)判定準(zhǔn)則( modal assurance criteria,MAC)對(duì)試驗(yàn)?zāi)B(tài)結(jié)果進(jìn)行置信度研究[12],最后對(duì)比仿真與試驗(yàn)結(jié)果,驗(yàn)證所建模型的可靠性。本文所建立的模型可對(duì)制動(dòng)盤使用狀態(tài)下的振動(dòng)特性進(jìn)行預(yù)測(cè)。
本文分析的制動(dòng)盤外徑0.276 m,凸臺(tái)內(nèi)徑0.070 m,高度0.047 m,構(gòu)建制動(dòng)盤CAD模型后,將其導(dǎo)入Hyperworks軟件,經(jīng)過幾何清理,去除細(xì)小邊,生成有限元模型如圖1所示。由于制動(dòng)盤頂部壁厚較小,因此劃分體網(wǎng)格大小約為2 mm,共508 710個(gè)單元,116 106個(gè)節(jié)點(diǎn)。制動(dòng)盤的材料為HT250,其參數(shù)選用文獻(xiàn)[3]中的材料參數(shù),E=117 000 MPa,μ=0.27,ρ=7×10-9t/mm3。

圖1 制動(dòng)盤有限元模型Fig.1 Finite element model of brake disc
制動(dòng)盤自由狀態(tài)下各階模態(tài)如圖2所示。自由狀態(tài)下,制動(dòng)盤前五階彈性模態(tài)頻率分別為1 126.9、1 129.2、2 618.8、2 619.3、2 624.4 Hz。

圖2 自由狀態(tài)下制動(dòng)盤各階模態(tài)仿真結(jié)果Fig.2 Modal simulation results of brake disc in the free state
由圖2可以看出,制動(dòng)盤的模態(tài)振型可分為制動(dòng)摩擦面以及摩擦面與凸臺(tái)組合的變形,圖2(a)、2(b)的振型為盤邊的四階彎曲,二者為重根模態(tài),除振型之間存在一定的角度差,其模態(tài)頻率、模態(tài)振型基本相同。圖2(c)、2(d)為盤邊的六階彎曲,二者也是重根模態(tài),區(qū)別在于其節(jié)徑與x軸之間的夾角不同。圖2(e)為制動(dòng)盤盤體的呼吸振型,該振型存在一個(gè)節(jié)圓。從上述結(jié)果可以看出,前四階模態(tài)為摩擦面模態(tài),第五階為摩擦面與凸臺(tái)的組合模態(tài),仿真表明更高階的摩擦面模態(tài)還包括盤邊八階彎曲、十階彎曲等,更高階的摩擦面與凸臺(tái)的組合模態(tài)包括凸臺(tái)的伸縮變形等,因此在制動(dòng)盤更高階的模態(tài)頻率分析中,摩擦面模態(tài)、摩擦面與凸臺(tái)的組合模態(tài)一定會(huì)交替出現(xiàn),而且也一定會(huì)出現(xiàn)更多的重根模態(tài)。
模態(tài)測(cè)試通過計(jì)算被測(cè)結(jié)構(gòu)的頻率響應(yīng)函數(shù),得到結(jié)構(gòu)的模態(tài)參數(shù)。通過模態(tài)測(cè)試,可以評(píng)價(jià)結(jié)構(gòu)的動(dòng)態(tài)特性,對(duì)結(jié)構(gòu)的振動(dòng)故障進(jìn)行診斷和預(yù)報(bào),判斷振動(dòng)發(fā)生的根本原因,驗(yàn)證有限元模型的準(zhǔn)確性。在模態(tài)測(cè)試前,首先要制定制動(dòng)盤的測(cè)量方案。
以制動(dòng)盤中心軸線為z軸,制動(dòng)盤的幾何形狀關(guān)于x、y坐標(biāo)軸對(duì)稱,其屬于典型的回轉(zhuǎn)體零件,因此在試驗(yàn)時(shí)可選擇圓柱坐標(biāo)系進(jìn)行建模,采用半徑r、角度θ、高度z這三個(gè)坐標(biāo)描述測(cè)試點(diǎn),制動(dòng)盤x、y軸方向的剛度,遠(yuǎn)大于z軸方向的剛度,前幾階振動(dòng)主要沿z軸產(chǎn)生,因此只采集各測(cè)點(diǎn)z方向的加速度即可。
制動(dòng)盤整體質(zhì)量較小,粘貼傳感器時(shí),應(yīng)注意減小附加質(zhì)量的影響,可采用多點(diǎn)激勵(lì)、多點(diǎn)拾振的方法進(jìn)行,激勵(lì)為z向,加速度傳感器的方向與圓柱坐標(biāo)系三個(gè)方向一致。由前述仿真結(jié)果可知,若完整描述盤邊四階彎曲,需要在盤邊測(cè)8個(gè)點(diǎn),而完整描述盤邊六階彎曲,需要在盤邊測(cè)12個(gè)點(diǎn),但上述點(diǎn)的角度、方位是不同的。因此為保證試驗(yàn)時(shí)得到完整的上述振型,取上述兩種振型測(cè)點(diǎn)數(shù)的最小公倍數(shù),即在整個(gè)圓周上布置24個(gè)測(cè)點(diǎn),一共布置5層,共120個(gè)激勵(lì)點(diǎn),測(cè)點(diǎn)的編號(hào)為逆時(shí)針編號(hào),由外向里,編號(hào)逐漸增加,如圖3所示。

圖3 制動(dòng)盤測(cè)點(diǎn)布置圖Fig.3 Measuring points graph of brake disc
試驗(yàn)時(shí)將制動(dòng)盤用彈性繩懸掛起來,如圖4所示。因?yàn)橹苿?dòng)盤固有頻率較高,利用lms impact模塊進(jìn)行測(cè)量,頻率范圍選為0~3000 Hz。試驗(yàn)時(shí)4個(gè)加速度傳感器沿周向均布,以防破壞結(jié)構(gòu)的對(duì)稱性。試驗(yàn)采用傳感器固定,移動(dòng)力錘的方式進(jìn)行,每敲擊一個(gè)測(cè)點(diǎn),采集一組數(shù)據(jù)。

圖4 制動(dòng)盤自由模態(tài)試驗(yàn)懸掛圖Fig.4 Free modes experimental graph of brake disc
利用lms polymax軟件對(duì)測(cè)試數(shù)據(jù)處理,得到的穩(wěn)態(tài)圖如圖5所示。

圖5 lms軟件分析所得穩(wěn)態(tài)圖Fig.5 Stabilization graph obtained by lms software analysis
從圖5中可以看出,制動(dòng)盤傳遞函數(shù)固有頻率非常明顯,選取穩(wěn)定點(diǎn)較多的位置,就可以確定制動(dòng)盤的固有頻率,進(jìn)一步計(jì)算得到比較明顯的模態(tài)振型,如圖6所示。由圖6及表1中的數(shù)據(jù),對(duì)比仿真得到的模態(tài)振型,仿真和試驗(yàn)所得模態(tài)振型相似、模態(tài)頻率接近,模態(tài)頻率的誤差小于3%,說明所建制動(dòng)盤有限元模型準(zhǔn)確度較高,該模型可用于約束狀態(tài)下制動(dòng)盤模態(tài)特性的預(yù)測(cè)。

圖6 自由狀態(tài)下制動(dòng)盤試驗(yàn)?zāi)B(tài)Fig.6 Experimental modes of brake disc in the free state

表1 制動(dòng)盤仿真模態(tài)與試驗(yàn)?zāi)B(tài)結(jié)果比較
模態(tài)驗(yàn)證就是對(duì)模態(tài)參數(shù)估計(jì)結(jié)果的正確性進(jìn)行檢驗(yàn)。除了經(jīng)驗(yàn)判斷之外,還有許多手段進(jìn)行驗(yàn)證。第一級(jí)驗(yàn)證方法是把試驗(yàn)得到的FRF與從模態(tài)參數(shù)識(shí)別過程中綜合得出的FRF進(jìn)行比較;第二級(jí)驗(yàn)證是用模態(tài)置信準(zhǔn)則(MAC)等數(shù)學(xué)工具來檢驗(yàn)估計(jì)的模態(tài)質(zhì)量。
試驗(yàn)?zāi)B(tài)分析對(duì)測(cè)量所得FRF進(jìn)行曲線擬合,能夠得到結(jié)構(gòu)的模態(tài)參數(shù),根據(jù)擬合確定的模態(tài)參數(shù)又可以得到綜合的FRF。在模態(tài)驗(yàn)證階段可對(duì)比試驗(yàn)FRF與綜合的FRF,得出二者的相關(guān)性與誤差,二者的一致性主要受分析帶寬、殘余項(xiàng)和噪聲等因素的影響。在制動(dòng)盤上隨機(jī)選擇兩個(gè)測(cè)點(diǎn),上述兩點(diǎn)的試驗(yàn)FRF與綜合FRF對(duì)比如圖7所示。結(jié)果顯示二者在分析帶寬內(nèi)一致性較好,驗(yàn)證了試驗(yàn)所得模態(tài)參數(shù)的正確性。

圖7 制動(dòng)盤測(cè)點(diǎn)試驗(yàn)FRF與綜合FRF對(duì)比Fig.7 Comparison of measured FRF and comprehensive FRF of brake disc measurement points
MAC是振型向量之間的點(diǎn)積,用于評(píng)價(jià)兩個(gè)模態(tài)振型向量幾何上的相關(guān)性,如果MAC值接近0,表明兩個(gè)振型向量之間是正交的;若MAC值接近1,說明這兩個(gè)振型向量彼此平行[13]。MAC矩陣中對(duì)角元素是各階振型向量與自身的點(diǎn)積,因此MAC值為1,非對(duì)角元素為不同振型向量之間的點(diǎn)積,MAC值接近于0。制動(dòng)盤振型MAC矩陣如圖8所示,從圖中可以看出MAC矩陣中非對(duì)角元素值很小,表明試驗(yàn)測(cè)試得到的制動(dòng)盤前五階模態(tài)就是制動(dòng)盤本身的固有特性。

圖8 制動(dòng)盤振型MAC矩陣Fig.8 Vibration shape MAC matrix of brake disc
(1)測(cè)量并建立了制動(dòng)盤有限元模型,仿真得到制動(dòng)盤自由狀態(tài)下前五階彈性模態(tài),其頻率分布范圍為1000~3000 Hz。在這五階模態(tài)中,有兩組重根模態(tài)。
(2)模態(tài)測(cè)試前,可以根據(jù)仿真模態(tài)所關(guān)心的振型測(cè)點(diǎn)的最小公倍數(shù),確定制動(dòng)盤每周所布置的測(cè)點(diǎn)數(shù)。宜采用圓柱坐標(biāo)系,多點(diǎn)激勵(lì)、多點(diǎn)拾振的方法對(duì)制動(dòng)盤等回轉(zhuǎn)體零部件進(jìn)行模態(tài)試驗(yàn)。通過試驗(yàn)?zāi)軠?zhǔn)確地測(cè)試得到自由狀態(tài)下制動(dòng)盤前五階的彈性模態(tài),包括摩擦面模態(tài),摩擦面與凸臺(tái)的組合模態(tài),其振型分別為制動(dòng)盤盤邊的四階彎曲、六階彎曲和制動(dòng)盤的呼吸振型。
(3)由試驗(yàn)?zāi)B(tài)參數(shù)得到的綜合FRF與試驗(yàn)FRF相吻合;模態(tài)置信準(zhǔn)則矩陣中非對(duì)角元素值很小。這些分析結(jié)果表明,試驗(yàn)測(cè)試得到的制動(dòng)盤前五階模態(tài)就是制動(dòng)盤本身的固有特性。試驗(yàn)結(jié)果與仿真結(jié)果相吻合,表明了所建仿真模型的準(zhǔn)確性。
從分析結(jié)果也可以看出,在制動(dòng)盤更高階的模態(tài)頻率分析中,摩擦面模態(tài)、摩擦面與凸臺(tái)的組合模態(tài)一定會(huì)交替出現(xiàn),而且也一定會(huì)出現(xiàn)更多的重根模態(tài)。