吳后吉 施 維 梁興雨
(1-廣東省智能制造研究所 廣東 廣州 510070 2-天津大學內燃機燃燒學國家重點實驗室)
活塞環與氣缸套是發動機最為關鍵的摩擦副之一,有研究表明[1],內燃機在工作過程中由氣缸套-活塞環摩擦副引起的摩擦損失占整個發動機摩擦損失的20%~50%,所以降低氣缸套-活塞環摩擦副的摩擦損失,對發動機動力性、經濟性具有重要意義。由于缸套-活塞環摩擦副運動時周圍的工作條件非常惡劣,所涉及到的因素非常繁雜,摩擦潤滑狀態的改變頻繁,每加入一個非線性影響因素進行多模型潤滑耦合分析都已經十分困難,但多模型潤滑耦合分析仍然是當今缸套-活塞環潤滑摩擦系統發展的主流方向。
本文針對不同發動機工作參數,諸如發動機轉速、負荷以及缸內壓力變化等影響因素為主要研究對象,再結合Reynolds 方程、微凸體接觸模型、以及表面粗糙度和潤滑油變粘度等多個非線性影響因素導入基本的缸套-活塞環流體潤滑模型中,建立流體潤滑耦合仿真模型,之后使用編程軟件結合Reynolds 邊界條件編制流體潤滑程序,對氣缸套-活塞環摩擦副進行流體潤滑多模型仿真耦合計算,在此基礎上,對不同發動機轉速,負荷,缸內壓力變化對活塞環潤滑性能進行探討,分析得出發動機工作參數對活塞環潤滑性能的影響,找到改善活塞環潤滑狀態與降低摩擦力所帶來的功率損耗的具體措施。
活塞環-缸套流體潤滑示意圖如圖1 所示。

圖1 活塞環-缸套流體潤滑示意圖
采用Patir 和Cheng[2-3]提出的平均雷諾方程進行計算:

式中:p 為平均油膜壓力;φx和φy分別為在x 和y 方向上的壓力流量因子;φs是剪切流量因子;h 為名義油膜厚度;μ 為潤滑劑動力粘度;U 為活塞運動速度,σ 為兩粗糙表面綜合粗糙度;hˉt為實際油膜厚度期望值采用Reynolds 粘度-溫度-壓力公式[4]計算。
變動力粘度:

式中:μ0為常溫下的動力粘度,壓強單位為帕斯卡,T為絕對溫度。
選活塞環與缸內最大距離處作為油膜起始點,仿真計算采用Reynolds 邊界條件,應滿足:

式中:xl(-b/2 <xl<b/2)為油膜破裂處的位置,可由求解Reynolds 方程時自動確定其位置。
由于考慮到粗糙度的存在,當膜厚比小于4 時,油膜比較薄,活塞環-氣缸套摩擦副處于混合潤滑狀態,2 表面微凸體發生接觸,在活塞環徑向受力平衡時,應考慮微凸體接觸力。假定微凸體高度服從高斯分布,曲率半徑為定值,可采用Greenwood[5]等提出的粗糙表面的接觸理論模型:

式中:對于工程表面來講,可取ηβσ=0.04,σ/β=10-3;Wa為微凸體接觸力;Ac為實際接觸面積;A 為名義接觸面積;E 為2 表面綜合彈性模量,F2.5(H)和F2(H)可根據參考文獻[6-8]得到。
使用SOR[9-10]法迭代求解Reynolds 方程,得到油膜厚度及油膜壓力分布后,可用下式計算摩擦力:

式中:τ0為材料剪切應力系數,可取τ0=2x10-6,α0為表面峰元剪切強度變化率,可取α0=0.1
考慮到發動機在進氣行程中點附近,這個時刻的活塞環相對速度較大,此時處于完全流體潤滑狀態,摩擦力低,摩擦磨損很小,工作條件較佳?;钊倪\動速度接近最大值,潤滑效果好,擠壓項可以忽略不計。參照這種思路,我們可以得到如下求解過程:首先選定進氣沖程中點附近的某一時刻,先估取一個最小油膜厚度值,給定初始油膜壓力及粘度、溫度等參數,經過修改最小油膜厚度迭代求解得到修正的最小油膜厚度和油膜壓力,然后接著進行下個時間步長(曲軸轉角)的計算,直到完成整個工作循環的計算,修正初始最小油膜厚度,直至整個循環始末最小膜厚值小于某個值時,程序停止。如圖2 所示。

圖2 計算流程圖
利用上述流體潤滑模型對某發動機進行仿真計算,該發動機活塞環組由三道環組成,分別是兩道氣環和一道油環。主要參數為:連桿長度L 為0.174 m,曲柄半徑R=0.075 5 m,活塞環半徑r=0.05m 缸套粗糙度σ=1.6 μm,活塞環表面粗糙度σ=0.8 μm。
圖3 表示的是在轉速為1 200 rpm、1 600 rpm、2 000rpm 的情況下發動機一個工作循環的最小油膜厚度曲線。橫坐標表示曲軸轉角,以進氣沖程階段開始計算,即曲軸轉角為零度時開始計算,縱坐標表示最小油膜厚度值。從圖3 中可以得出,當轉速增加時,在發動機每個沖程中,每個沖程的中間時刻的最小油膜厚度最大,因為此時活塞環刮油輪廓面與氣缸壁之間所產生的油膜擠壓變化率比較小,活塞在每個沖程階段的運動速度較大,潤滑狀態較佳,故而此時的最小油膜厚度為最大值。每個沖程階段中間時刻的最大的最小油膜厚度越大,發動機的最小油膜厚度值也就越大。在點火上止點時刻(在圖3 中曲軸轉角300°~400°這個范圍中)。

圖3 速度-最小油膜厚度關系圖
這個時刻的最小油膜厚度值是發動機整個工作循環中所產生的最小油膜厚度值當中的最小值,隨著轉速的增加,此時的最小油膜厚度值有變化,但不是很明顯,這是由于此時的最小油膜厚度值已經非常小,潤滑狀態為混合潤滑狀態,潤滑油膜難以形成或只能形成一層很薄的油膜,這時轉速的變化對這個時刻的最小油膜影響不是主要影響因素,潤滑狀態的轉化,發動機工作環境的惡劣性都會影響著油膜厚度的形成。最小油膜厚度的變化與活塞環徑向所受的外載有直接的關系,轉速的增加,造成活塞每次往返運動所產生的慣性力增加,最終結果導致活塞環徑向所受外力變化大,滿足所受合外力平衡條件決定著最小油膜厚度值的變化規律。
如圖4 所示,當轉速增加時,流體摩擦力也隨著增大,這是因為提高發動機轉速,潤滑油的內摩擦力與摩擦熱流也增大,潤滑油溫度升高,造成異常磨損。這也就說明了并非轉速越高越好,轉速的提高,雖然對加大潤滑油膜的厚度比較有利,但隨之帶來的由潤滑油產生的摩擦力也會增加,造成了摩擦功的損耗加大,降低發動機產生的熱能利用率。

圖4 速度-摩擦力關系圖
如圖5 所示,在轉速為1 600rpm 時為前提條件下,表示的是發動機在整個工作循環中不同負荷條件下缸內壓力隨曲軸轉角變化的關系圖。在轉速,粘度等參數不變的情況下,選取發動機負荷為40%、70%、100%3 個工作狀態進行討論。由圖5 分析,負荷的改變,對接近點火上止點附近的缸內壓力影響較大,而在其余曲軸轉角下的缸壓則沒有太大的影響,并且隨著負荷的增大,接近點火上止點附近的缸內壓力也隨之增大。

圖5 不同負荷的缸內壓力圖
圖6 表示的是發動機在整個工作循環中不同負荷條件下最小油膜厚度隨曲軸轉角變化的關系圖。圖中,給定一個發動機處于一定的負荷值,在進氣沖程階段,活塞運行到行程中間時刻的最小油膜厚度比其他時刻的最小油膜厚度值大,在點火上止點附近的最小油膜厚度值最小,且3 條不同負荷時的曲線在點火上止點附近的最小油膜厚度值相差不大,這是因為,負荷的改變,對缸內壓力的影響非常小,只影響這缸壓曲線的峰值處,在其余的情況40%的負荷、70%的負荷與100%的負荷的缸壓曲線幾乎重合。

圖6 負荷-最小油膜厚度關系圖
圖7 表示的是不同負荷下摩擦力隨曲柄轉角變化的關系圖,從圖中可以明顯地看到,改變發動機工作負荷對摩擦力幾乎沒有影響。

圖7 負荷-摩擦力關系圖
使用增壓技術對發動機進行增壓,利用增壓器將空氣或可燃混合氣進行預增壓,再送入氣缸進行反復燃燒,以達到提高發動機功率和改善發動機的經濟性。但是,發動機在采用廢氣渦輪增壓技術后,工作中產生的最高爆發壓力和平均溫度將大幅度提高,從而使發動機的機械性能、潤滑性能都會受到影響。所以,我們在進行理論計算時,考慮到增壓后的發動機進行分析研究,分析增壓后對潤滑性能的影響是很有必要的。在這里只分析活塞環在運動過程中的潤滑狀態進行分析。
通過在轉速為1 600 rpm 時為前提條件下,圖8表示了加載渦輪增壓器和沒有加載渦輪增壓器時發動機整個工作循環的最小油膜厚度曲線圖,橫坐標表示的是曲軸轉角,從進氣沖程時刻開始計算,縱坐標表示的是最小油膜厚度值。

圖8 缸內增壓-最小油膜厚度關系圖
增壓后的最高爆發壓力和溫度有所增大,在發動機一個工作循環中,增壓前和增壓后的發動機,除了發動機最高爆發壓力有所改變外,其余工作狀態下的缸內燃燒壓力變化很小,計算得到的最小油膜厚度值幾乎相等,而在發動機最高爆發壓力處,也就是接近點火上止點處,增壓后的最小油膜厚度值比正常情況下的最小油膜厚度值略小,原因是增壓后對發動機主要結構部件、材料、工藝等方面產生了更高的要求,而在沒有改變這些發動機零部件的情況下,增壓只能使發動機承受更高缸內混合氣壓縮和燃燒氣體的溫度和壓力以及燃燒室受熱零件的熱負荷,造成此時缸內工作環境有所惡化。而在圖9 中計算活塞環所受摩擦力的關系圖中,增壓后與增壓前活塞環所受摩擦力幾乎沒有發生變化。

圖9 缸內增壓-摩擦力關系圖
1)在缸套-活塞環摩擦副潤滑摩擦分析中,雖然每加入一個非線性影響因素進行多模型潤滑耦合分析已經十分困難,但多模型潤滑耦合分析仍然是當今缸套-活塞環潤滑摩擦系統發展的主流方向。
2)發動機的負荷變化和加載增壓器后的缸內壓力的變化對活塞環潤滑性能沒有多大影響。
3)不同轉速狀態下,當轉速增加時,最小油膜厚度增大,有利于提高潤滑性能;當轉速增加時,摩擦力也隨著增大,這是因為提高發動機轉速,潤滑油的內摩擦力與摩擦熱流也增大,潤滑油溫度升高,造成異常磨損。這也就說明了并非轉速越高越好,應當有個最佳值。轉速的提高,雖然對加大潤滑油膜的厚度比較有利,但隨之帶來的由潤滑油產生的摩擦力也會增加,造成了摩擦功熱的損耗加大,降低發動機產生的熱能的利用率。