陳曉++陳小兵++龔艷

摘要：機(jī)械化高效施藥是目前對(duì)噴藥機(jī)械的基本要求，所以風(fēng)送植保機(jī)械被普遍使用。風(fēng)機(jī)是其中的關(guān)鍵部件，利用PROE和ANSYS這2種軟件對(duì)強(qiáng)力風(fēng)送裝備中多翼型離心風(fēng)機(jī)葉輪進(jìn)行實(shí)體建模、靜力學(xué)分析和模態(tài)分析。通過靜力學(xué)分析，得到葉輪在工作時(shí)的等效云圖，分析得到葉輪的最大應(yīng)力出現(xiàn)在葉片與前盤等的焊接處；通過模態(tài)分析，得到葉輪的前八階模態(tài)的振型云圖，驗(yàn)證了葉輪工作的安全性。利用有限元分析方法，可以提高結(jié)果的準(zhǔn)確性，為優(yōu)化設(shè)計(jì)提供理論依據(jù)。

關(guān)鍵詞：離心風(fēng)機(jī)；葉輪；有限元分析

中圖分類號(hào)： S491文獻(xiàn)標(biāo)志碼： A

文章編號(hào)：1002-1302（2017）09-0180-04

風(fēng)送植保機(jī)械由于具有作業(yè)效率高、霧滴在靶標(biāo)物上沉積均勻、農(nóng)藥利用率高和流失少、對(duì)土壤和地下水等環(huán)境污染小等特點(diǎn)，在歐美發(fā)達(dá)國(guó)家得到廣泛應(yīng)用。噴霧技術(shù)是施藥作業(yè)的關(guān)鍵，最理想的施藥效果就是將農(nóng)藥霧滴均勻地噴灑在靶標(biāo)上，由于靶標(biāo)是立體的，霧滴很難穿透到達(dá)指定的地點(diǎn)，而風(fēng)送噴霧技術(shù)的應(yīng)用大大增加了霧滴穿透性，它是利用氣流的動(dòng)能把藥液霧滴吹送到指定的位置，大大地改善了藥液的霧化性能，增強(qiáng)了霧滴穿透性，提高了霧滴在植物中的沉積，成為了改善藥液分布的主要措施[1]。

離心風(fēng)機(jī)是組成風(fēng)送植保機(jī)械噴灑系統(tǒng)的關(guān)鍵部件。離心風(fēng)機(jī)工作時(shí)，葉輪通過磁電轉(zhuǎn)機(jī)子直接接在汽油機(jī)曲軸上，隨著發(fā)動(dòng)機(jī)高速轉(zhuǎn)動(dòng)，風(fēng)機(jī)產(chǎn)生的高速氣流分3部分流出，一部分氣流經(jīng)前半蝸的冷卻風(fēng)口流出用以冷卻發(fā)動(dòng)機(jī)；少部分氣流經(jīng)蝸殼上部的擋風(fēng)板進(jìn)入藥箱，在藥液表面形成正壓，藥液就在風(fēng)壓作用下經(jīng)輸液短管送到噴頭處；大部分高速氣流經(jīng)風(fēng)機(jī)出口進(jìn)入噴管，在噴頭喉口處流速最大，將流出的藥液[LL]沖擊成極細(xì)的霧滴，并隨氣流噴出[2]。

離心風(fēng)機(jī)部件主要包括風(fēng)機(jī)前殼、后殼和葉輪。葉輪是離心風(fēng)機(jī)的關(guān)鍵部件之一。其中葉片和前后盤的幾何參數(shù)及形狀都對(duì)離心風(fēng)機(jī)的性能有很大影響[2]。同樣，隨著植保機(jī)械技術(shù)的提高，對(duì)風(fēng)機(jī)強(qiáng)度方面的要求也越來越高。通過對(duì)葉輪進(jìn)行有限元分析可以提高葉輪設(shè)計(jì)的準(zhǔn)確性以及提高優(yōu)化葉輪的效率。

1葉輪的三維實(shí)體建模

葉輪是風(fēng)機(jī)的旋轉(zhuǎn)部件，是風(fēng)機(jī)傳遞給氣體能量的唯一元件，離心風(fēng)機(jī)的葉輪形狀、尺寸和數(shù)量，決定了它的工作能力和性能參數(shù)。離心風(fēng)機(jī)的葉輪一般由前盤、后盤、葉片和軸盤等組成[3]。

1.1葉輪的幾何參數(shù)

本研究主要是對(duì)意大利TIFONE公司的CANNON 80S型強(qiáng)力風(fēng)送噴霧機(jī)上的多翼型離心風(fēng)機(jī)進(jìn)行有限元分析。葉輪的幾何參數(shù)如表1所示。

1.2葉輪實(shí)體建模

葉輪的實(shí)體模型通過PROE軟件來實(shí)現(xiàn)，葉輪的輪盤和葉片是通過焊接固定在一起的，由于不考慮焊接對(duì)葉輪的影響，所以在建模時(shí)可以將二者看為一個(gè)整體。最終的實(shí)體模型如圖1所示。

完成建模后，將PROE默認(rèn)的PRT格式的文件轉(zhuǎn)化為可以導(dǎo)入到ANSYS Workbench中的IGS文件。IGS是CAD文件的一種通用格式，可以有效地避免不能成功導(dǎo)入而帶來的麻煩。

葉輪材料參數(shù)如表2所示。

2風(fēng)機(jī)葉輪的靜力分析

本研究的離心風(fēng)機(jī)應(yīng)用于大型植保機(jī)械上，它的功用是產(chǎn)生高速氣流，將藥液破碎霧化或?qū)⑺幏鄞瞪ⅲ⒅拖蜻h(yuǎn)方。風(fēng)機(jī)的主要部件包括前殼、后殼和葉輪。而葉輪是風(fēng)機(jī)的旋轉(zhuǎn)部件，是風(fēng)機(jī)傳遞給氣體能量的唯一元件，離心風(fēng)機(jī)的性能參數(shù)尤為重要[4]?？梢酝ㄟ^有限元分析來判斷風(fēng)機(jī)葉輪的強(qiáng)度。

2.1葉輪網(wǎng)格的劃分

ANSYS中網(wǎng)格的結(jié)構(gòu)和數(shù)量將會(huì)影響計(jì)算結(jié)果的精度和計(jì)算規(guī)模的大小。一般來講，網(wǎng)格數(shù)量增加，計(jì)算精度會(huì)有所提高，但同時(shí)計(jì)算規(guī)模也會(huì)增加，所以在確定網(wǎng)格數(shù)量時(shí)應(yīng)權(quán)衡2個(gè)因素綜合考慮。網(wǎng)格疏密是指在結(jié)構(gòu)不同部位采用大小不同的網(wǎng)格，在計(jì)算數(shù)據(jù)變化梯度較大的部位（如應(yīng)力集中處），為了較好地反映數(shù)據(jù)變化規(guī)律，需要采用比較密集的網(wǎng)格。而在計(jì)算數(shù)據(jù)變化梯度較小的部位，為減小模型規(guī)模，則應(yīng)劃分相對(duì)稀疏的網(wǎng)格。

在ANSYS Workbench中，主要有自動(dòng)劃分網(wǎng)格、四面體網(wǎng)格劃分、六面體網(wǎng)格劃分和掃掠法等。本研究采用自動(dòng)網(wǎng)格劃分的方法，對(duì)在PROE中進(jìn)行建模的三維葉輪模型進(jìn)行網(wǎng)格劃分，得到的有限元網(wǎng)格模型如圖2所示。

2.2施加工作載荷與約束

根據(jù)葉輪的實(shí)際工作狀況，葉輪隨風(fēng)機(jī)軸轉(zhuǎn)動(dòng)。葉輪沒有軸向移動(dòng)并且輪盤上的6個(gè)螺栓孔表面轉(zhuǎn)動(dòng)時(shí)也沒有周向位移。所以對(duì)葉輪靜力分析時(shí)施加載荷的約束條件為：6個(gè)螺栓孔表面沒有軸向和周向的位移，即定義螺栓孔表面的自由度為0。

葉輪在工作中主要受到離心力和氣動(dòng)力的影響，由于氣動(dòng)力對(duì)葉輪產(chǎn)生的影響與離心力相比要小得多，所以在研究時(shí)只考慮離心力對(duì)葉輪的影響[2]。離心力可以通過施加角速度來定義，其角速度為210 rad/s。

2.3靜力學(xué)分析結(jié)果

經(jīng)過ANSYS的計(jì)算求解后，得到了葉輪在工作時(shí)的等效應(yīng)力圖、等效應(yīng)變?cè)茍D和總變形量云圖。圖3所示為葉輪的等效應(yīng)力圖，葉輪的應(yīng)力主要集中于葉片與前盤等的焊接處，最大值為137.37 MPa，小于材料的屈服強(qiáng)度355 MPa，可以認(rèn)為葉輪滿足強(qiáng)度要求。葉輪在總體上所呈現(xiàn)的趨勢(shì)是離旋轉(zhuǎn)軸越遠(yuǎn)，應(yīng)力越大。而葉片與輪盤焊接的中部應(yīng)力更為集中，此外在6個(gè)螺栓孔的位置周圍也存在應(yīng)力集中現(xiàn)象。

圖4、圖5分別為葉輪的等效應(yīng)變?cè)茍D和總變形量云圖。從等效應(yīng)變?cè)茍D中可以看出，應(yīng)力集中的地方也是應(yīng)變較大的地方，應(yīng)變的最大值為6.148 8×10-4 mm。在總變形量云圖中，位移沿徑向方向成圓環(huán)行分布，位移最大的位置是在葉片的外側(cè)邊緣上，最大的位移為0.387 32 mm。在應(yīng)變?cè)茍D中，每個(gè)葉片的變形基本上是相同的，基本上呈現(xiàn)的是循環(huán)對(duì)稱的特點(diǎn)。在總變形量云圖中，可以發(fā)現(xiàn)葉輪輪盤的變形量比葉片的變形量小很多，葉片的變形出現(xiàn)了扭曲的現(xiàn)象。

通過葉輪的分布云圖可知，應(yīng)力最大處發(fā)生在葉片與輪盤的焊接處，可以改進(jìn)焊接工藝來減少應(yīng)力集中，以此來提高葉輪的強(qiáng)度，增加風(fēng)機(jī)的壽命。

3模態(tài)分析

模態(tài)分析主要是為了確定機(jī)械鼻尖的振動(dòng)特性，包括固有頻率和振型，為避免發(fā)生共振提供參考。

離心風(fēng)機(jī)在旋轉(zhuǎn)工作時(shí)，葉輪的激振源主要是離心力的激振頻率。葉輪的結(jié)構(gòu)相對(duì)復(fù)雜，直接計(jì)算本身的固有頻率比較困難。借助ANSYS對(duì)其進(jìn)行模態(tài)分析，可以得到相應(yīng)模態(tài)的振型，可以在風(fēng)機(jī)設(shè)計(jì)時(shí)為避免共振而提供參考的依據(jù)。

模態(tài)分析是研究結(jié)構(gòu)動(dòng)力特性的一種方法。模態(tài)是機(jī)械結(jié)構(gòu)的固有振動(dòng)特性，采用在靜力學(xué)分析中運(yùn)用的載荷和約束，對(duì)葉輪進(jìn)行預(yù)應(yīng)力模態(tài)分析。隨著模態(tài)分析階數(shù)的不斷增大，模態(tài)分析的結(jié)果誤差將會(huì)變大[5]。所以根據(jù)需要，本研究只進(jìn)行前八階的模態(tài)分析，結(jié)果如表3所示。

由表3可以看出，葉輪的一階和二階模態(tài)的振動(dòng)頻率非常接近，應(yīng)該是振動(dòng)方程的二重根。它們的振型也是相似的，如圖6、圖7所示，分別為葉輪的一階模態(tài)和二階模態(tài)振型云圖，它們?cè)诠舱駮r(shí)存在一個(gè)相位差，即共振時(shí)振幅變化的位置不同。在這兩階模態(tài)下，葉輪的振型可以看成是由2個(gè)近似對(duì)稱的部分組成，離對(duì)稱線越遠(yuǎn)，其位移量越大。圖中的位移數(shù)值只是葉輪在固有頻率下振動(dòng)時(shí)各點(diǎn)位移的比值，并不是實(shí)際的位移值。

圖8、圖9分別是葉輪的三階模態(tài)振型云圖和四階模態(tài)振型云圖，其徑向位移成圓環(huán)狀向外遞增，最大位移出現(xiàn)在葉輪邊緣和葉片外側(cè)邊緣。

圖10、圖11分別是葉輪的五階模態(tài)和六階模態(tài)振型云圖，它們的振動(dòng)頻率相近，應(yīng)該是振動(dòng)方程的兩重根。它們的振[CM（25]型也是相似的，只是在共振時(shí)存在一個(gè)相位差。輪盤的最大位移出現(xiàn)在輪盤邊緣的2個(gè)區(qū)域，這2個(gè)區(qū)域近似對(duì)稱，葉片的最大位移出現(xiàn)在其外側(cè)邊緣上。

圖12是葉輪的七階模態(tài)振型云圖，葉輪的振型可以看成是由2個(gè)近似對(duì)稱的部分組成，離對(duì)稱線越遠(yuǎn)，其位移量越大。

圖13是葉輪的八階模態(tài)振型云圖，葉輪的振型近似對(duì)稱分布，但是位移量最大的點(diǎn)出現(xiàn)在葉輪邊緣的一側(cè)，此時(shí)，葉輪輪盤也存在扭轉(zhuǎn)振動(dòng)。

綜合前八階的振型云圖可知，最大得到相對(duì)位移都出現(xiàn)在輪盤的外側(cè)，有效地避免了葉片發(fā)生破壞?？梢栽黾虞啽P的厚度來提高輪盤的剛度。

離心風(fēng)機(jī)葉輪的干擾頻率與葉輪的轉(zhuǎn)速有關(guān)，其關(guān)系如下：

式中：n為轉(zhuǎn)速，r/min；f為干擾頻率，Hz。本研究的葉輪轉(zhuǎn)速為2 000 r/min，由上述公式可得到葉輪的干擾頻率為 33.3 Hz。由前八階的模態(tài)分析可以得到葉輪工作時(shí)有效地避免了共振現(xiàn)象，保證其安全。

4結(jié)論

葉輪是風(fēng)機(jī)的關(guān)鍵部件，傳統(tǒng)的強(qiáng)度計(jì)算是將葉輪進(jìn)行相應(yīng)的轉(zhuǎn)化，采用相應(yīng)的公式進(jìn)行計(jì)算，計(jì)算結(jié)果與實(shí)際情況存在較大的誤差[6]。本研究采用PROE和ANSYS對(duì)植保機(jī)械上的離心風(fēng)機(jī)進(jìn)行了靜力學(xué)分析和模態(tài)分析，為葉輪的優(yōu)化設(shè)計(jì)提供了依據(jù)。運(yùn)用PROE和ANSYS對(duì)產(chǎn)品進(jìn)行優(yōu)化，可以減少強(qiáng)度和模態(tài)分析所需要的時(shí)間，降低優(yōu)化所需的成本，為優(yōu)化設(shè)計(jì)提供有效的數(shù)據(jù)。由靜力學(xué)分析和模態(tài)分析的結(jié)果可以得知該葉輪的結(jié)構(gòu)和材料滿足工作要求。葉片與輪盤焊接處的應(yīng)力比較集中，在進(jìn)行優(yōu)化設(shè)計(jì)時(shí)，可以考慮增加相應(yīng)的厚度來提高結(jié)果強(qiáng)度；葉片外側(cè)邊緣的位移較大，在優(yōu)化設(shè)計(jì)時(shí)應(yīng)考慮相應(yīng)情況進(jìn)而減少其位移。

通過ANSYS軟件靜力分析可以很直觀地了解葉輪最大應(yīng)力出現(xiàn)的位置，便于在設(shè)計(jì)過程中對(duì)葉輪進(jìn)行改進(jìn)優(yōu)化。利用ANSYS對(duì)葉輪進(jìn)行模態(tài)分析可以得到葉輪各階的固有頻率和相對(duì)應(yīng)的振型，為避免風(fēng)機(jī)發(fā)生共振提供了有力的理論依據(jù)。而采用傳統(tǒng)的計(jì)算方法，所得到的最大應(yīng)力只是平均應(yīng)力，不是實(shí)際最大應(yīng)力。所以采用傳統(tǒng)的計(jì)算方法所得到的結(jié)果通常比實(shí)際的結(jié)果小。雖然通過有限元分析方法進(jìn)行的輕度校核與實(shí)際也存在一定誤差，但比傳統(tǒng)計(jì)算方法所得的誤差小。在進(jìn)行葉輪設(shè)計(jì)的過程中，建議采用有限元法對(duì)風(fēng)機(jī)進(jìn)行強(qiáng)度校核，但風(fēng)機(jī)的最終優(yōu)化結(jié)果還需要經(jīng)過試驗(yàn)實(shí)踐的驗(yàn)證。

參考文獻(xiàn)：

[1]劉建，呂新民，黨革榮，等. 植保機(jī)械的研究現(xiàn)狀與發(fā)展趨勢(shì)[J]. 西北農(nóng)林科技大學(xué)學(xué)報(bào)（自然科學(xué)版），2003，31（增刊1）：202-204.

[2]張玉成，儀登利，馮殿義，等. 通風(fēng)機(jī)設(shè)計(jì)與選型[M]. 北京：化學(xué)工業(yè)出版社，2011.

[3]黃炎，趙滿全，黃鵬飛. 前置式雙圓盤割草機(jī)主軸有限元分析[J]. 農(nóng)機(jī)化研究，2013，35（8）：26-27.

[4]蔡兆麟，韓海燕，彭鑫. 旋轉(zhuǎn)葉輪三維應(yīng)力分布及其改善措施[J]. 華中科技大學(xué)學(xué)報(bào)（自然科學(xué)版），2003，31（12）：81-83.

[5]鄭磊，尹健，紀(jì)斌，等. 基于ANSYS的脫離滾筒模態(tài)分析[J]. 農(nóng)機(jī)化研究，2013，35（4）：48-49.

[6]李燕平. 離心風(fēng)機(jī)葉輪的有限元分析及優(yōu)化[D]. 楊凌：西北農(nóng)林科技大學(xué)，2012：23-24.

[7]凌桂龍，丁金濱，溫正. ANSYS Workbench 15.0從入門到精通[M]. 北京：清華大學(xué)出版社，2012：114-115.