本發明屬于內燃機燃油系統，涉及一種船用噴射共軌泵的柱塞副動態磨損量計算方法，涉及一種船用噴射共軌泵中柱塞與套筒磨損量的計算方法，尤其涉及一種基于油膜壓力迭代與微凸體接觸的柱塞副動態磨損量計算方法。
背景技術：
1、高壓共軌泵是船舶發動機燃料供給系統的核心單元，其性能直接影響發動機的動力輸出、燃油經濟性與排放水平。在高壓共軌泵中，柱塞與套筒構成的柱塞副是最為關鍵的精密摩擦副之一。柱塞副在工作過程中承受極高的周期性壓力（通常可達160mpa以上），并在微米級配合間隙內進行高速往復運動，其摩擦學狀態與磨損進程直接決定了共軌泵的容積效率、壓力保持能力與使用壽命。
2、隨著全球航運業對清潔能源需求的日益增長，甲醇作為一種低碳燃料，在船用發動機中的應用逐漸廣泛。然而，甲醇燃料的引入也為發動機關鍵部件帶來了新的技術挑戰，尤其是燃料供給系統中精密偶件的磨損問題。甲醇燃料的物化特性與傳統船用燃料存在顯著差異，其親水性強，易與潤滑油發生互溶，導致潤滑油出現乳化現象。乳化會顯著降低潤滑油的黏度和油膜強度，使其潤滑性能大幅下降。在采用自潤滑設計的柱塞副中，由于缺乏持續外部潤滑油補給，潤滑條件更為苛刻。甲醇引起的潤滑油乳化進一步加劇了邊界潤滑或干摩擦狀態，加速材料磨損。
3、柱塞副磨損研究在此背景下顯得尤為重要。柱塞與套筒之間的微米級配合間隙若磨損增大，將導致燃油泄漏增加、噴射壓力下降、霧化質量惡化，進而引起發動機功率降低、能耗上升和排放超標。長期運行中，柱塞副磨損還可能引發突發性故障，造成停航維修，帶來巨大的經濟損失。因此，準確預測和評估柱塞的磨損量，對實現狀態監測、優化維護周期、提高系統可靠性和壽命具有重要意義。
4、然而，現有技術針對船用高壓共軌泵柱塞副磨損量的計算方法研究仍存在一些局限：
5、1）現有磨損模型多基于穩態假設，未充分考慮柱塞在復雜運動過程中的微運動及其對油膜狀態與接觸載荷的動態影響，導致磨損預測與實際情況存在偏差；
6、2）缺乏對表面形貌影響的系統建模，多數方法未將柱塞與套筒的表面粗糙度納入磨損計算中，無法準確反映混合潤滑狀態下微凸體接觸對磨損的貢獻；
7、3）磨損計算與多工況耦合能力不足，傳統研究多分析單一因素對摩擦學特性的作用，忽略了微運動、表面形貌的耦合作用，現有方法往往未整合凸輪轉速、入口壓力、表面形貌等多因素對磨損過程的綜合影響，難以實現變工況下的動態磨損量評估。
8、因此，基于現有設計方法的不足，針對船用噴射共軌泵的柱塞副磨損量計算，本領域的研究人員需要發明新的計算方法，旨在克服上述缺陷。這些方法應能更精準地描述復雜曲面接觸狀態，幫助優化設計和制造過程，從而提高裝備的可靠性、穩定性及其使用壽命。
技術實現思路
1、有鑒于此，本發明為了解決現有磨損模型未整合微運動、表面形貌等多因素對磨損過程的綜合影響，難以實現變工況下的動態磨損量評估的問題，提供一種船用噴射共軌泵的柱塞副動態磨損量計算方法。
2、為達到上述目的，本發明提供如下技術方案：
3、一種船用噴射共軌泵的柱塞副動態磨損量計算方法，包括以下步驟：
4、s1、建立針對復雜服役工況的動力學模型
5、輸入柱塞副的運動初始條件，建立柱塞副動力學模型，并計算柱塞在運動過程中所受的各力和力矩大小；
6、s2、基于偏心距與傾斜角度的油膜壓力與油膜厚度迭代計算
7、輸入柱塞相對于柱塞套的初始偏心距和初始傾斜角度；
8、基于所述初始偏心距和初始傾斜角度，計算初始油膜厚度分布；
9、基于所述初始油膜厚度分布，計算對應的油膜壓力分布；
10、判斷由所述油膜壓力分布計算得到的油膜承載力與柱塞所受側向載荷之間的平衡誤差是否滿足收斂要求；
11、若否，則修正所述偏心距和傾斜角度，并基于修正后的值再次執行計算油膜厚度分布的步驟，進行迭代計算；
12、若是，則輸出當前迭代下的最終油膜厚度分布和油膜壓力分布；
13、s3、計入表面粗糙度的混合潤滑模型修正與迭代
14、計算柱塞與柱塞套接觸表面的綜合粗糙度；
15、基于所述綜合粗糙度，對雷諾模型進行修正，以表征混合潤滑狀態；
16、使用修正后的雷諾模型，重新計算油膜壓力分布；
17、判斷重新計算后，油膜承載力與側向載荷的平衡誤差是否滿足收斂要求；
18、若否，則返回步驟s2修正偏心距和傾斜角度的步驟，并基于修正后的值再次執行計算油膜厚度分布的步驟，進行再次迭代計算；
19、若是，則輸出修正模型下的最終油膜厚度分布和油膜壓力分布；
20、s4、柱塞副動態磨損量計算
21、基于步驟s3輸出的最終油膜壓力分布與油膜厚度分布，計算微凸體接觸壓力；
22、結合archard磨損模型，基于微凸體接觸壓力和相對滑移距離，計算柱塞副的磨損量。
23、本發明的有益效果在于：
24、1、本發明所公開的船用噴射共軌泵的柱塞副動態磨損量計算方法，通過建立柱塞副動力學模型，求解偏心距和傾斜角度隨凸輪轉角的動態方程（采用龍格-庫塔法），首次將柱塞的瞬態微運動量化。在此基礎上，將微運動參數耦合到油膜厚度與壓力的迭代求解中，使油膜狀態（承載力）能夠實時響應運動與受力變化。這種“動力學模型-油膜狀態”的閉環耦合，從根本上改變了將接觸條件視為靜態的傳統思路，使磨損計算能夠模擬真實工作循環中的動態行為，顯著提高了預測的時域準確性和工況真實性。
25、2、本發明所公開的船用噴射共軌泵的柱塞副動態磨損量計算方法，系統性地引入了表面形貌參數，通過公式?計算接觸副的綜合粗糙度，并創新性地使用膜厚比作為潤滑狀態判據，在判定為混合潤滑后，采用平均流量模型，引入壓力流量因子、剪切流量因子和接觸流量因子對雷諾方程進行修正，使油膜壓力計算包含了粗糙度對潤滑油流動的影響。基于greenwood-tripp微凸體接觸理論計算混合潤滑區的實際接觸壓力，將宏觀的油膜潤滑與微觀的粗糙峰接觸統一到一個模型中，使磨損計算機制從“光滑表面假設”躍升至“真實表面接觸”，能更真實地反映混合潤滑下微凸體相互作用導致的磨損。
26、3、本發明所公開的船用噴射共軌泵的柱塞副動態磨損量計算方法，綜合集成了結構參數（柱塞直徑、長度等）、運動與工況參數（凸輪轉速、油膜入口壓力）、界面特性（表面粗糙度）以及求解得到的動態變量（偏心距、傾斜角、油膜厚度）。方法通過迭代算法自然納入了這些因素的相互影響（例如，凸輪轉速變化會改變滑移速度與動力，進而影響微運動和油膜狀態）。因此，該方法能夠通過調整輸入參數，直接模擬不同工作條件（如不同轉速、不同燃料壓力、不同加工表面質量）下的磨損進程，實現變工況的動態磨損量預測。這為從“離線設計評估”轉向“在線狀態監測與壽命預測”提供了可量化、可調參的全面理論模型和計算依據。
27、本發明的其他優點、目標和特征在某種程度上將在隨后的說明書中進行闡述，并且在某種程度上，基于對下文的考察研究對本領域技術人員而言將是顯而易見的，或者可以從本發明的實踐中得到教導。本發明的目標和其他優點可以通過下面的說明書來實現和獲得。