發(fā)動機(jī)機(jī)油濾清器濾層強度分析及優(yōu)化

劉攀勇, 陳振雷, 商 晴, 柳超杰, 林枝強

發(fā)動機(jī)機(jī)油濾清器濾層強度分析及優(yōu)化

劉攀勇, 陳振雷*, 商 晴, 柳超杰, 林枝強

（寧波大學(xué) 海運學(xué)院, 浙江 寧波 315832）

針對機(jī)油濾清器工作工況下進(jìn)出口壓差、機(jī)油濾層強度及導(dǎo)流樁高度等問題, 通過試驗測試與仿真相結(jié)合, 對濾清器初步設(shè)計進(jìn)行了評估及優(yōu)化, 以確保濾清器在工作工況下進(jìn)出口壓降及濾層強度能滿足要求. 首先進(jìn)行濾層性能試驗, 得到濾層的慣性阻力系數(shù)和黏性阻力系數(shù); 再通過濾層多孔介質(zhì)CFD分析, 對濾清器進(jìn)出口壓降進(jìn)行分析計算. 結(jié)果表明: 在-18℃、25℃和70℃的工況下, 進(jìn)出口壓降都小于10kPa, 滿足相關(guān)要求. 針對濾層的最大主應(yīng)力超過其抗拉強度的問題, 通過CAE仿真分析, 優(yōu)化濾層與導(dǎo)流樁間隙, 將濾層最大主應(yīng)力由110.1MPa降至36.99MPa, 小于其抗拉強度42.8MPa.

濾清器; 濾層; 黏性阻力系數(shù); 慣性阻力系數(shù); 強度分析

發(fā)動機(jī)的機(jī)油過濾器能優(yōu)先濾除對發(fā)動機(jī)磨損最有害的金屬屑和無機(jī)雜質(zhì)等物質(zhì), 從而能有效降低發(fā)動機(jī)摩擦副的磨損, 延長發(fā)動機(jī)使用壽命, 同時還能較好地改善發(fā)動機(jī)機(jī)油品質(zhì), 達(dá)到延長發(fā)動機(jī)壽命的作用[1-3]. 所以機(jī)油濾清器必須具有濾清能力強, 流通阻力小, 使用壽命長等特點[4]. 但影響濾清器過濾性能的主要是濾層, 如果濾紙所受主應(yīng)力大于其抗拉強度, 則會破壞濾紙, 導(dǎo)致濾層喪失過濾能力. 更為嚴(yán)重的是, 材質(zhì)較差或損壞的濾紙產(chǎn)生大量紙屑, 造成發(fā)動機(jī)潤滑系統(tǒng)油路堵塞或零部件異常, 導(dǎo)致發(fā)動機(jī)損壞或報廢.

因此大量學(xué)者針對濾清器過濾問題做了許多研究. Kaukopaasi等[5]在研究中證實, 濾清器多層復(fù)合濾紙在過濾效率和使用壽命上明顯優(yōu)于單層濾紙, 但是復(fù)合多少層可以達(dá)到最好效果還有待研究. 梁云等[6]用纖維梯度結(jié)構(gòu)設(shè)計理念制備了具有梯度過濾功能的濾紙, 但是要想過濾精度高, 孔徑就要變小, 阻力就會變大, 反之精度就會降低. 李亨等[7]綜合對比了多孔介質(zhì)-純流體耦合流動區(qū)域的數(shù)值模擬方法及現(xiàn)狀, 認(rèn)為多孔介質(zhì)模型可以用于多孔介質(zhì)中流體流動的數(shù)值模擬, 且此方法后續(xù)一直延用在濾層的CFD計算中.

本文針對過濾器進(jìn)出口壓降及濾紙強度不足的問題, 在完成濾紙試驗測試的基礎(chǔ)上, 通過CFD計算出在-18℃、25℃和70℃工況下的進(jìn)出口壓降以及濾紙所受最大主應(yīng)力, 最終通過優(yōu)化分析并完善了濾清器的下殼體以及改變?yōu)V層與導(dǎo)流樁間距, 在保證進(jìn)出口油壓的同時, 減小了濾紙受到的最大主應(yīng)力.

1 理論分析

1.1 數(shù)學(xué)模型

機(jī)油濾清器在工作流量范圍內(nèi)一般采用湍流模型, 在計算時采用Reynolds時均方程法對模型進(jìn)行計算[8]. 流體運動滿足質(zhì)量守恒定律和動量守恒定律, 基本瞬時控制方程組由質(zhì)量守恒方程和動量守恒方程組成[9].

質(zhì)量守恒方程:

動量守恒方程:

流體流過濾層(多孔介質(zhì))時遵守Darcy定律[10], 當(dāng)流經(jīng)均勻的多孔介質(zhì)時, 動量方程的源項由黏性損失項和慣性損失項組成, 計算公式如下[11]:

將上式在濾層厚度方向進(jìn)行積分, 可得:

1.2 邊界條件

2 模型建立

2.1 網(wǎng)格模型建立

使用Hypermesh軟件對機(jī)油濾清器進(jìn)行四面體網(wǎng)格劃分[12], 最小網(wǎng)格尺寸為0.1mm, 最大網(wǎng)格尺寸為5.12mm, 網(wǎng)格總數(shù)量為226萬, 濾層的厚度為1.2mm. 初始設(shè)計中, 濾層和導(dǎo)流樁之間的距離為0.4mm, 濾清器下殼體內(nèi)有8個長方體的導(dǎo)流樁, 當(dāng)濾層受力時, 導(dǎo)流樁對濾層起到支撐作用, 濾清器模型如圖1所示.

2.2 各工況邊界條件

當(dāng)濾清器處于工作狀態(tài), 流經(jīng)的機(jī)油油溫不同時, 機(jī)油的各項參數(shù)也有所不同. 3種不同溫度下濾清器中機(jī)油的各項相關(guān)參數(shù)見表1.

表1 各溫度下的機(jī)油參數(shù)

按設(shè)計要求, 過濾的機(jī)油在不同溫度條件下,流經(jīng)濾層時對濾清器進(jìn)出口油壓有一定壓差要求, 同時也對濾層所受最大主應(yīng)力有一定的要求. 設(shè)計的具體工況條件見表2.

表2 各工況邊界條件

由表2可知, 在-18℃且進(jìn)口流量為4L?min-1時, 進(jìn)出口壓降要小于15kPa; 在25℃且進(jìn)口流量為15L?min-1時, 進(jìn)出口壓降要小于10kPa; 在70℃且進(jìn)口流量為20L?min-1時, 進(jìn)出口壓降要小于15kPa. 并且需要滿足在25℃溫度下, 當(dāng)濾層上下壓差為85kPa負(fù)壓時, 濾層所受最大主應(yīng)力不能超過其本身的抗拉強度42.8MPa.

3 濾層性能測試

3.1 試驗原理

CFD計算進(jìn)出口油壓和濾層表面壓差時, 計算軟件對于多孔介質(zhì)所需要的輸入?yún)?shù)中, 缺失濾層的慣性阻力系數(shù)和黏性阻力系數(shù). 為獲得準(zhǔn)確的流體力學(xué)分析參數(shù), 需進(jìn)行濾層性能試驗[13], 以得到對應(yīng)濾層的慣性阻力系數(shù)和黏性阻力系數(shù).

根據(jù)式(4), 在已知機(jī)油黏性系數(shù)、密度以及濾層厚度的條件下, 若要得到濾層的黏性阻力系數(shù)和慣性阻力系數(shù), 需要試驗得到迎面風(fēng)速和濾層兩側(cè)的壓降.

試驗在規(guī)定壓差條件下, 測試一定時間內(nèi)垂直通過試樣給定面積的氣流流量, 計算出透氣率, 再將透氣率換算成迎面風(fēng)速, 即1L?(m2s)-1為10-3m?s-1. 具體測試數(shù)據(jù)見表3. 表3中的風(fēng)速值為每組試驗壓降下5次試驗風(fēng)速的平均值, 以確保數(shù)據(jù)的可靠性.

表3 風(fēng)速與壓降測試數(shù)據(jù)

3.2 試驗數(shù)據(jù)處理

根據(jù)上述表3數(shù)據(jù), 通過數(shù)據(jù)擬合可以得到以下壓降與風(fēng)速的表達(dá)式:

4 結(jié)果及優(yōu)化分析

4.1 結(jié)果分析

通過濾層的性能試驗, 得到了相對準(zhǔn)確的黏性阻力系數(shù)和慣性阻力系數(shù). 通過Fluent軟件計算得到初始設(shè)計下表2中3個工況的進(jìn)出口壓降分別為3.2、3.9和4.6kPa. 以-18℃的工況為例, 進(jìn)出口壓降云圖如圖3所示, 且3個工況下的進(jìn)出口壓降全部小于10kPa, 符合要求.

圖3 -18℃工況下進(jìn)出口壓降圖

但在25℃溫度下, 濾層上下壓差為85kPa時, 濾層的最大主應(yīng)力云圖如圖4所示. 由圖可見, 在濾層與導(dǎo)流樁接觸部位的最大主應(yīng)力為110.1MPa, 超過濾層抗拉強度42.8MPa, 不符合要求, 需要對其進(jìn)行優(yōu)化分析.

圖4 濾層最大主應(yīng)力云圖

圖5為濾層位移云圖, 由圖可知, 濾層的最大位移量為1.565mm, 但在初始設(shè)計中濾層和導(dǎo)流樁的間隙僅0.4mm. 所以在工作狀態(tài)時, 導(dǎo)流樁對濾層出現(xiàn)過度支撐, 導(dǎo)致濾層產(chǎn)生高達(dá)110.18MPa的應(yīng)力.

圖5 濾層位移云圖

4.2 優(yōu)化分析

導(dǎo)流樁在機(jī)油濾清器工作的過程中起著承擔(dān)濾層部分應(yīng)力的作用. 如果導(dǎo)流樁數(shù)量較少, 與濾層的間隙過大, 導(dǎo)流樁將不能對濾層產(chǎn)生有效支撐, 從而導(dǎo)致濾層在其邊緣產(chǎn)生過大應(yīng)力. 如果導(dǎo)流樁與濾層的間隙過小, 則會使得導(dǎo)流樁對濾層產(chǎn)生過度支撐, 導(dǎo)致濾層接觸局部產(chǎn)生過大應(yīng)力. 因此, 需要對導(dǎo)流樁的形狀及其與濾層的間隙進(jìn)行優(yōu)化, 使其滿足項目目標(biāo).

針對圖4導(dǎo)流樁局部產(chǎn)生應(yīng)力過大的情況, 對導(dǎo)流樁進(jìn)行優(yōu)化. 圖6為優(yōu)化后的下殼體模型, 初始設(shè)計的8個長方體導(dǎo)流樁改進(jìn)為23個表面帶有過渡圓弧狀的導(dǎo)流樁. 23個導(dǎo)流樁均勻分布為6行, 每行3~4個, 出油口兩端各分布1個, 均勻分散來自濾層的應(yīng)力. 同時導(dǎo)流樁的布設(shè)應(yīng)滿足基本導(dǎo)流作用, 避免因?qū)Я鳂斗植歼^密影響出油效率.

圖6 優(yōu)化后下殼體模型

采用Abaqus軟件[14], 進(jìn)行濾層應(yīng)力最優(yōu)設(shè)計計算分析. 分析表明, 逐步增加濾層與導(dǎo)流樁間隙將導(dǎo)致濾層中部與導(dǎo)流樁接觸的部位應(yīng)力減少, 同時濾層邊緣部位應(yīng)力增加. 間隙量達(dá)到2mm, 且濾層上下壓差為85kPa條件時, 濾層邊緣最大主應(yīng)力為36.99MPa, 濾層中部與導(dǎo)流樁接觸應(yīng)力為14.04MPa, 均小于濾層材料抗拉強度42.8MPa的要求, 符合條件. 圖7為優(yōu)化后的對應(yīng)應(yīng)力云圖. 圖8為優(yōu)化后濾層的位移云圖, 濾層的最大位移量為2.03mm, 此時濾層和導(dǎo)流樁剛好接觸, 最大應(yīng)力出現(xiàn)在濾層邊緣部位.

最后通過Fluent軟件計算得到表2的3個工況下的進(jìn)出口壓降分別為4.8、4.5和5.0kPa, 也滿足要求. 圖9為25℃工況下優(yōu)化后濾清器進(jìn)出口壓降云圖.

5 結(jié)語

針對過濾器進(jìn)出口壓降及濾層強度不足的問題, 在通過濾層性能測試獲得計算參數(shù)的基礎(chǔ)上, 進(jìn)行流場與強度仿真, 并通過優(yōu)化分析, 在不增加成本的條件下, 滿足了進(jìn)出口壓降及濾層強度要求, 結(jié)論如下:

圖7 優(yōu)化后濾層最大主應(yīng)力云圖

圖8 優(yōu)化后濾層位移云圖

圖9 優(yōu)化后25℃工況下進(jìn)出口壓降圖

(2)通過濾層多孔介質(zhì)CFD分析, 對濾清器進(jìn)出口壓降進(jìn)行分析計算. 結(jié)果表明在-18℃、25℃和70℃的3個工況下, 進(jìn)出口壓降都小于10kPa, 滿足要求.

(3)針對濾層最大主應(yīng)力超過其抗拉強度的問題, 通過Abaqus仿真分析, 優(yōu)化了濾層與導(dǎo)流樁間隙, 將濾層最大主應(yīng)力由110.1MPa降至36.99 MPa, 小于其抗拉強度42.8MPa.

[1] 韋唐凌, 鐘康海, 韋尚勇, 等. 基于CFD技術(shù)的離心式機(jī)油濾清器結(jié)構(gòu)優(yōu)化[J]. 內(nèi)燃機(jī)與配件, 2019(11):4-8.

[2] 趙海蓉, 龔濤, 王佳, 等. 機(jī)油濾清器流體仿真與結(jié)構(gòu)優(yōu)化[J]. 過濾與分離, 2017, 27(2):38-41.

[3] Alie A, Darwito P A. Improve of engine oil lifetime by using additional filter: A case study at PT Vale Indonesia Tbk[C]. AIP Conference Proceedings, Surabaya, Indonesia, 2019.

[4] 朱坤. 詳解機(jī)油濾清器[J]. 汽車維護(hù)與修理, 2015(2): 88-89.

[5] Kaukopaasi J G, Shah N B. Multi layered sheet: The media of the future for automotive filtration applications [J]. TAPPI Journal, 2002(1):127-132.

[6] 梁云, 胡健, 周雪松. 汽車用梯度結(jié)構(gòu)過濾紙的研制[J]. 中國造紙, 2005, 24(8):18-21.

[7] 李亨, 張錫文, 何楓. 論多孔介質(zhì)中流體流動問題的數(shù)值模擬方法[J]. 石油大學(xué)學(xué)報(自然科學(xué)版), 2000, 24 (5):111-116.

[8] 李佳, 劉震濤, 劉忠民, 等. 空氣濾清器流動過程仿真與試驗分析[J]. 浙江大學(xué)學(xué)報(工學(xué)版), 2012, 46(2): 327-332.

[9] 劉連明. 汽車空氣濾清器流場分析與結(jié)構(gòu)優(yōu)化[D]. 南昌: 江西農(nóng)業(yè)大學(xué), 2018.

[10] 張惠, 王國麗. 濾層仿真參數(shù)的確定[J]. 內(nèi)燃機(jī)與配件, 2010(9):18-20.

[11] Rybakov Y P, Semenova N V, Adam G, et al. Generalized Darcy’s law in filtration theory[J]. EPJ Web of Conferences, 2018, 173:2017-2020.

[12] 李偉, 王珂. 基于CFD的濾材仿真參數(shù)研究[J]. 汽車零部件, 2011(10):85-87.

[13] 曾鋒. 基于HyperMesh軟件的客車CAE前處理二次開發(fā)[J]. 客車技術(shù)與研究, 2019, 41(4):25-27.

[14] GB/T 14295-2019. 空氣過濾器[S].

[15] Chen D R, Pui D Y H, Liu B Y H. Optimization of pleated filter designs using a finite-element numerical model[J]. Aerosol Science and Technology, 1995, 23(4): 579-590.

Oil filter layer stress analysis and optimization

LIU Panyong, CHEN Zhenlei*, SHANG Qing, LIU Chaojie, LIN Zhiqiang

( Faculty of Maritime and Transportation, Ningbo University, Ningbo 315832, China )

Aiming at the problems concerning pressure difference between inlet and outlet, strength of oil filter layer and diversion pile height under the working condition of oil filter, the preliminary design of filter is evaluated and optimized by the combination of experiment and simulation, so that the pressure drop of inlet and outlet and the strength of filter layer are ensured to meet the requirements. First, the filter layer performance test is performed to obtain the inertial resistance coefficient and the viscous resistance coefficient of the filter layer. Then the pressure drop at the inlet and outlet of the filter is analyzed and calculated by CFD (Computational Fluid Dynamics) analysis of the porous media of the filter layer. The results show that, under three working conditions of -18℃, 25℃ and 70℃, the inlet and outlet pressure drop is always less than 10kPa, which meets the requirements. In view of the problem that the maximum principal stress of the filter layer exceeds its tensile strength, CAE simulation analysis is used to optimize the gap between the filter layer and the diversion piles, making the maximum principal stress of the filter layer being reduced from 110.1MPa to 36.99MPa, which is less than its tensile strength of 42.8MPa.

filter; filter layer; viscous drag coefficient; inertial drag coefficient; stress analysis

TK432

A

1001-5132（2021）01-0045-05

2019?12?05.

寧波大學(xué)學(xué)報（理工版）網(wǎng)址: http://journallg.nbu.edu.cn/

國家自然科學(xué)基金（11972203）; 寧波大學(xué)人才工程項目（理）.

劉攀勇（1995－）, 男, 安徽巢湖人, 在讀碩士研究生, 主要研究方向: 流體與動力機(jī)械結(jié)構(gòu)仿真. E-mail: 1070409780@qq.com

陳振雷（1963－）, 男, 浙江寧波人, 博士/教授, 主要研究方向: 高端制造仿真分析. E-mail: chenzhenlei@nbu.edu.cn

（責(zé)任編輯 章踐立）