陳科良,邱 滿,杜民獻(xiàn),夏天軍,吳吉平,熊勇剛
(湖南工業(yè)大學(xué) 機(jī)械工程學(xué)院,湖南 株洲 412007)
以往研究表明,發(fā)動機(jī)振動是產(chǎn)生噪聲的主要原因,發(fā)動機(jī)的噪聲主要是其殼體表面的輻射噪聲[1-2]。數(shù)值分析與試驗(yàn)相結(jié)合的研究方法,已被廣泛應(yīng)用于結(jié)構(gòu)振動與輻射噪聲的預(yù)測。本文基于這種方法,首先對結(jié)構(gòu)的模態(tài)參數(shù)是否正確進(jìn)行判斷,并進(jìn)行有限元結(jié)果分析[3-4],以便得到殼體更準(zhǔn)確的模態(tài)參數(shù);再把模態(tài)分析的結(jié)果與模態(tài)試驗(yàn)結(jié)果相比較,以確定發(fā)動機(jī)殼體Pro/E模型建立的有效性以及試驗(yàn)方法的正確性[5]。
試驗(yàn)采用的動態(tài)信號分析儀QLVSA-3為重慶大學(xué)測試中心開發(fā)的QLV系列虛擬儀器,最高測量頻率選擇為3000Hz,儀器如圖1所示。
圖1 動態(tài)信號分析儀Fig.1 Dynamic signal analyzer
在發(fā)動機(jī)殼體表面布置24個測試點(diǎn),待發(fā)動機(jī)轉(zhuǎn)速穩(wěn)定后,依次對各測試點(diǎn)進(jìn)行振動測量,并在測試點(diǎn)上施加激勵力,同時測量其他測試點(diǎn)的響應(yīng)。經(jīng)過擬合、模態(tài)綜合和數(shù)據(jù)處理得到傳遞函數(shù)后,便能獲得反映被測構(gòu)件動態(tài)特性的各階固有頻率和振型等模態(tài)參數(shù)。測試點(diǎn)分布如圖2所示。
圖2 測試點(diǎn)分布Fig.2 Test points distribution
發(fā)動機(jī)殼體模態(tài)試驗(yàn)測得的前四階固有頻率如表1所示。將表1中的結(jié)果與自由模態(tài)分析結(jié)果作對比,以驗(yàn)證試驗(yàn)方法的正確性。
表1 試驗(yàn)測得固有頻率Table1 The natural frequency measured in experiment
在有限元仿真運(yùn)算中,首先是建立一個接近于實(shí)體的三維模型,再通過有限元前處理模塊建立合理的有限元模型,最后通過有限元分析軟件進(jìn)行分析求解。由于ANSYS軟件在建模時有一定的局限性,因此采用Pro/E進(jìn)行殼體建模。要建立與實(shí)物完全相同的模型難度較大,在建立三維模型時作一些修改:如原實(shí)物上的2個類似橢圓的部分均為不規(guī)則曲面,建模時簡化成平面[6];忽略實(shí)物邊上的細(xì)小部件。最后建立的三維模型如圖3所示。
圖3 殼體三維模型Fig.3 Shell three-dimensional modal
在Pro/E中建立三維實(shí)體模型并保存為IGES格式文件,導(dǎo)入ANSA中進(jìn)行前處理,建立有限元模型。為了能建立起有限元網(wǎng)格,在對有限元計(jì)算精度影響不大的情況下,對實(shí)體模型進(jìn)行簡化:殼體邊上的螺栓孔采用直孔處理;有限元前處理時忽略因倒圓角產(chǎn)生的尖角[7],對少許微小的倒圓角直接改成直角或把倒圓角改大,忽略一些局部細(xì)小結(jié)構(gòu)。建立有限元模型主要是為了分析有限元的動態(tài)響應(yīng),不需要進(jìn)行應(yīng)力分析,故不需考慮結(jié)構(gòu)的應(yīng)力,因此可忽略殼體結(jié)構(gòu)上一些細(xì)小的孔。
考慮到發(fā)動機(jī)殼體是典型的薄殼結(jié)構(gòu)[8],所以進(jìn)行網(wǎng)格劃分時,選擇的是20節(jié)點(diǎn)SOILD186單元,劃分網(wǎng)格的單元尺寸大小為4mm。根據(jù)以往的研究[9],為保證結(jié)果的精確性,此類殼體一般選用四面體單元計(jì)算,共49026個。最終得到合格的發(fā)動機(jī)殼體有限元模型,如圖4所示。發(fā)動機(jī)殼體為鋁合金材料,材料參數(shù)如表2所示。
表2 殼體材料參數(shù)Table2 Shell material parameters
自由模態(tài)分析時,不考慮約束的影響,只需計(jì)算出結(jié)構(gòu)本身的固有頻率[10]。利用 ANSYS 軟件提供的Block Lanczos法對4種單元尺寸的有限元模型的自由模態(tài)進(jìn)行求解,Block Lanczos 法可以提高計(jì)算效率和精度[11]。從Block Lanczos法求解的結(jié)果中提取發(fā)動機(jī)殼體模態(tài)分析的前四階自由模態(tài)固有頻率,如表3所示。
表3 前四階自由模態(tài)固有頻率Table2 Free mode natural frequencies for the first four orders
取模態(tài)試驗(yàn)得到的前四階固有頻率與自由模態(tài)得到的前四階固有頻率進(jìn)行對比,以驗(yàn)證有限元模型建立的正確性,如表4所示。
表4 試驗(yàn)?zāi)B(tài)固有頻率與自由模態(tài)固有頻率對比Table4 Natural frequencies comparison between experimental mode and free mode
由表中的數(shù)據(jù)可以看出,雖然試驗(yàn)數(shù)據(jù)有一定的誤差,但所測的前四階固有頻率與有限元自由模態(tài)計(jì)算的固有頻率基本吻合,這說明該有限元模型符合要求,可以進(jìn)行下一步分析。圖5為一至四階自由模態(tài)的位移等值線云圖。
圖5 自由模態(tài)的位移等值線云圖Fig.5 The displacement isoline nephogram for free mode
由圖5a可知,殼體右上角部分有沿高度方向向上的振動,對應(yīng)的左下角也有相似的振動,殼體表面上下2個類似橢圓的連接處產(chǎn)生了一定的彎曲。由圖5b可知,發(fā)動機(jī)殼體邊上有2個部位振動較強(qiáng)烈,振動部位有所轉(zhuǎn)移,但與一階模態(tài)時振動基本相似;其模態(tài)振型表現(xiàn)為彎曲模態(tài),左右兩側(cè)往上翹起,彎曲部位在殼體最外層小橢圓外側(cè)。由圖5c可知,三階模態(tài)振動情況比較復(fù)雜,殼體邊緣呈S形變形,振動部位較多。由圖5d可知,其振動部位相比一階和二階模態(tài)有所上移,殼體表面出現(xiàn)較強(qiáng)的振動,表現(xiàn)在殼體各曲面上和它們的交接處;彎曲部位在曲面的交接處。
采用與實(shí)際相符的約束條件對結(jié)構(gòu)進(jìn)行約束模態(tài)分析,不同約束條件下,結(jié)構(gòu)的固有頻率和振型都會發(fā)生變化。約束模態(tài)與實(shí)際情況越接近,對發(fā)動機(jī)殼體的振動噪聲研究越有意義。約束模態(tài)分析時,對殼體的11個螺孔和殼體邊界施加約束[12]。殼體的約束如圖6所示。
圖6 殼體的約束Fig.6 Shell constraint
在其他初始條件與自由模態(tài)分析相同的情況下,為了形成更精確的對比,本文采用模態(tài)分析的前八階固有頻率。約束模態(tài)與自由模態(tài)固有頻率的對比見表5和圖7。
表5 約束模態(tài)與自由模態(tài)固有頻率對比Table5 The natural frequencies comparison for constraint mode and free mode
圖7 約束模態(tài)和自由模態(tài)固有頻率對比Fig.7 The natural frequencies comparison for constraint mode and free mode
從表5和圖7可以看出,殼體在施加約束后固有頻率發(fā)生了較大變化,各階固有頻率均有較大提高。
圖8為發(fā)動機(jī)殼體約束模態(tài)分析得到的第一階和第八階位移等值線云圖。
圖8 約束模態(tài)的位移等值線云圖Fig.8 The displacement isoline nephogram for the constraint mode
對比約束模態(tài)和自由模態(tài)位移等值線云圖可知,第一階約束模態(tài)相比第一階自由模態(tài),因施加約束,殼體螺栓處振動變?nèi)酰饕駝訁^(qū)域上移,振動最強(qiáng)烈的地方在殼體3個曲面的連接部分表面和連接處;第八階約束模態(tài)相比第八階自由模態(tài),殼體振動情況發(fā)生改變,但主要振動區(qū)域還是大體相同,殼體上最外層類似橢圓的中心處振動依然強(qiáng)烈。以上說明,自由模態(tài)分析僅能反映殼體的振動形態(tài),而約束模態(tài)卻能夠在一定程度上改變其固有頻率和振動形態(tài)[13]。
1) 通過試驗(yàn)測試數(shù)據(jù)與模態(tài)分析數(shù)據(jù)的比較,驗(yàn)證了基于Pro/E-ANSA三維模型建立的可行性和試驗(yàn)方法的正確性。
2)對比約束模態(tài)和自由模態(tài),前者比后者對發(fā)動機(jī)振動特性的影響更明顯,殼體施加約束會大大提高其固有頻率,并且約束下的低階振型受到較大影響,而高階振型影響較小。約束模態(tài)下雖提高了殼體的固有頻率,改變了殼體的振型,但振型并沒有削弱。
3)相同工況下,殼體振動的大小取決于機(jī)體結(jié)構(gòu)本身的剛度。所以,提高殼體的剛度是最直接有效地減小殼體振動的措施,其主要方法有增加壁厚、加筋等。
[1] 夏天軍,陳科良,熊勇剛. 發(fā)動機(jī)殼體輻射噪聲試驗(yàn)研究[J]. 湖南工業(yè)大學(xué)學(xué)報(bào),2013,27(2):45-47.Xia Tianjun,Chen Keliang,Xiong Yonggang. Experimental Investigation on Engine Shell Radiation Noise[J]. Journal of Hunan University of Technology,2013,27(2):45-47.
[2] 楊 誠,周 科,陳 旭. 發(fā)動機(jī)殼體輻射噪聲預(yù)測[J].江蘇大學(xué)學(xué)報(bào):自然科學(xué)版,2010,31(4):393-396.Yang Cheng,Zhou Ke,Chen Xu. Prediction of Noise Radiated from Engine Shel[J]. Journal of Jiangsu University:Natural Science Edition,2010,31(4):393-396.
[3] 崔燕娟,汪志義,王 忠,等. 485柴油機(jī)機(jī)體有限元模擬及結(jié)構(gòu)改進(jìn)[J]. 小型內(nèi)燃機(jī)與摩托車,2009,38(5):61-64.Cui Yanjuan,Wang Zhiyi,Wang Zhong,et al. 485 Diesel Engine Finite Element Simulation and Improvement[J]. Small Internal Combustion Engine and Motorcycle,2009,38(5):61-64.
[4] 奚志朋,孫 平,靳永軍,等. 柴油機(jī)機(jī)體有限元結(jié)構(gòu)分析[J]. 拖拉機(jī)與農(nóng)用運(yùn)輸車,2006,33(1):24-26.Xi Zhipeng,Sun Ping,Jin Yongjun,et al. Finite Element Analysis (FEA) for Cylinder Block Structure of Diesel Engine[J]. Tractor & Farm Transporter,2006,33(1):24-26.
[5] 郭常立,張保成,趙丕歡,等. 柴油機(jī)整機(jī)結(jié)構(gòu)有限元模態(tài)分析[J]. 車用發(fā)動機(jī),2008(6):37-40.Guo Changli,Zhang Baocheng,Zhao Pihuan,et al. Diesel Engine Structure Finite Element Modal Analysis[J]. Vehicle Engine,2008(6):37-40.
[6] 王智興,樊文欣. 發(fā)動機(jī)噪聲測試環(huán)境的研究[J]. 車用發(fā)動機(jī),2005(2):49-51.Wang Zhixing,F(xiàn)an Wenxin. Enviromental Research of Engine Nosie Testing[J]. Vechicle Engine,2005(2):49-51.
[7] 沈慧軍,黃鼎友,殷承維. 發(fā)動機(jī)噪聲的聲強(qiáng)測試與分析[J]. 噪聲與振動控制,2008,28(4):132-135.Shen Huijun,Huang Dianyou,Yin Chengwei. Testing Sound Intensity and Analyzing Noise of Engine[J]. Noise and Vibration Control,2008,28(4):132-135.
[8] 張 磊,任海軍,艾曉玉,等. 發(fā)動機(jī)薄壁件結(jié)構(gòu)振動優(yōu)化[J]. 振動工程學(xué)報(bào),2010,23(6):687-697.Zhang Lei,Ren Haijun,Ai Xiaoyu,et al. Vibration Optimization of Engine Covers[J]. Journal of Vibration Engineering,2010,23(6): 687-697.
[9] 熊勇剛,李 順,夏天軍,等. 發(fā)動機(jī)殼體噪聲輻射控制研究[J]. 湖南工業(yè)大學(xué)學(xué)報(bào),2012,26(3):44-47.Xiong Yonggang,Li Shun,Xia Tianjun,et al. Study on Controlling Noises Radiated from Engine Shell[J]. Journal of Hunan University of Technology,2012,26(3):44-47.
[10]姜紹忠,閻文兵,陳 越,等. 柴油機(jī)表面輻射噪聲控制研究[J]. 機(jī)械與電子,2011(5):71-72.Jiang Shaozhong,Yan Wenbing,Chen Yue,et al. Radiation Noise Control of Diesel Engine[J]. Mechanical and Electronic,2011(5):71-72.
[11]李忠杰,王樹宗,練永慶. 某活塞發(fā)動機(jī)機(jī)體組合結(jié)構(gòu)模態(tài)分析[J]. 艦船科學(xué)技術(shù),2007,29(4):125-127.Li Zhongjie,Wang Shuzong,Lian Yongqing. Modal Analysis on Engine Block Combination Structure of a Piston Engine[J]. Ship Science and Technology,2007,29(4):125-127.
[12]金 巖,郝志勇. 車用汽油機(jī)噪聲源的聲強(qiáng)測量識別[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2006,37(2):14-16.Jin Yan,Hao Zhiyong. Investigation on Main Noise Sources of Vehicle Gasoline Engine by Sound Intensity Measurement[J]. Transactions of the Chinese Society for Agricultural Machinery,2006,37(2):14-16.
[13]王震武,胡志新,袁森柱,等. 摩托車噪聲源測量與降噪方法研究[J]. 傳感器與微系統(tǒng),2010,29(7):37-39.Wang Zhenwu,Hu Zhixin,Yuan Senzhu,et al. Noise Source Measurement and Study on Noise Reduction of Motorcycle[J]. Transducer and Microsystem Technologies,2010,29(7):37-39.