劉玉慧中航工業南京輕型航空動力有限公司
Abaqus有限元軟件在活塞發動機氣缸體強度計算中的應用
劉玉慧中航工業南京輕型航空動力有限公司
氣缸體是發動機的基礎零件和骨架,由于其結構和工作受力情況復雜,用常規的方法無法對其進行準確的受力分析。針對某型活塞發動機在試驗時氣缸體斷裂的事故,采用有限元數值計算方法,應用Abaqus有限元軟件,對其在工作載荷下的靜強度進行了分析,并根據計算結果對氣缸體進行了結構優化,優化后的方案得到了臺架試驗的驗證。
氣缸體;有限元;強度
氣缸體是發動機的基礎零件和骨架,通過它把發動機的曲柄連桿機構(包括活塞、連桿、曲軸、飛輪等零件)和配氣機構(包括缸蓋、凸輪軸等)以及供油、潤滑、冷卻等機構連接成一個整體。氣缸體工作在高載荷、磨損劇烈的條件下,承受較大的壓力,受力復雜,作用在缸蓋上的燃氣壓力通過缸蓋螺栓傳遞到缸體上,作用在活塞上的氣缸壓力產生了活塞向缸控右側的側壓力,此力以曲軸中心為支點,產生氣缸體的側翻力,曲軸的曲柄、連桿軸頸、連桿大頭的離心力作用對各主軸承孔蓋形成了彎矩,氣缸體承受彎矩。氣缸體結構復雜、壁厚差懸殊,易出現應力集中現象,因此,缸體的可靠性是人們在發動機設計開發階段要首先考慮的問題。
某型活塞發動機在開發試驗階段,氣缸體出現了斷裂事故。通過限元法對氣缸體進行了詳細分析,為缸體的優化設計提供了理論依據,并通過了臺架試驗的考核。有限元方法不僅節約了開發成本,還縮短了發動機的開發周期。
建立準確的數學模型是確保有限元數值模擬準確性的首要條件,包括有限元網格模型,邊界和載荷條件,所建的模型應能如實反映結構的幾何形狀、材料特性、邊界條件及承載情況等。特別是對復雜的工程問題作精細分析時,計算模型的建立至關重要。
2.1 結構及材料參數
采用軟件UG完成氣缸體三維實體模型。氣缸體材料為ZL101A,熱處理狀態T6,抗拉強度為275MPa;氣缸套材料為硼鑄鐵。
2.2 有限元模型
為了保證有限元計算的準確性,僅對計算精度影響較小的部分圓角進行適當簡化,采用10節點四面體實體單元劃分網格,對斷裂處的結構進行了局部細化,劃分網格后的單元總數為353212個,節點總數為670984個。
2.3 邊界條件及載荷
計算時考慮氣缸套與氣缸體之間的過盈裝配,過盈量為0.17mm;考慮溫度載荷,氣缸套均溫160℃,氣缸體均溫120℃;考慮力學載荷,取氣缸體爆發壓力最大及活塞對缸體側推力最大(曲軸轉動10°)時作為工作載荷,其中最大爆發壓力加載在缸體與缸蓋之間連接螺栓位置處;側推力加載在活塞裙部對應的缸體位置處,計算時不考慮缸蓋與氣缸體密封墊片的作用。計算載荷見表1。
約束氣缸體與缸蓋接觸面的法向位移,約束螺栓孔的徑向和軸向位移。

表1 計算載荷
2.4 變形及應力計算結果
氣缸體在計算載荷下最大位移量為0.43mm。氣缸體在計算載荷下的應力云圖見圖1,最大應力為828MPa,最大應力處由于受到螺栓孔約束的影響而產生較大的應力,此時該處的應力-應變關系已經不是線性關系,但由于有限元程序在計算時將結構變形當成彈性變形,因此該處的應力值不能反映真實情況,計算應力值大于真實值。不考慮螺栓孔位置處的應力,最大應力值為217MPa,最大應力位置為氣缸體試驗斷裂處。

圖1 氣缸體應力分布圖
氣缸體靜強度安全系數計算公式為:

式中:σb為材料的抗拉強度極限,σmax為最大當量應力。參考航空零部件安全性的一般要求,強度儲備系數n>1.5。根據式(1)可以得到n=1.26<1.5,表明氣缸體靜強度不滿足要求。
從應力分布可以看出,氣缸體在第一和第二散熱片之間(對應螺栓孔位置)出現應力集中現象,且加強筋在靠近散熱片位置應力偏大。根據計算分析結果,對氣缸體結構進行改進,在第一和第二散熱片之間局部增厚,并將加強筋加寬增厚。
對優化后的結構做同樣載荷及邊界條件下的靜強度有限元分析,氣缸體在第一和第二散熱片之間區域的最大應力由217MPa下降至130MPa,滿足靜強度要求,且改進后的結構相對原結構在加強筋及散熱片區域平均應力更低。
有限元法是一種先進的數值計算方法,通過有限元法計算分析了某型活塞發動機氣缸體的變形及應力情況,并對其進行了靜強度校核。實踐證明,此分析結果與實際情況比較吻合。
根據計算結果確定了氣缸體結構強度薄弱區,對其結構進行了優化,計算結果表明,優化后的結構應力值下降了40%,滿足靜強度要求。優化后的氣缸體裝機后在后續試驗及小批量生產中,均未出現斷裂故障。
[1]林海英,龍向陽.某型增壓航空活塞發動機仿真研究.航空動力學報,2009年06期
[2]馬斌,邵微,陳秉智.柴油發動機活塞靜強度分析.計算機輔助工程,2014年02期
劉玉慧(1982-),女,湖南省邵陽人,碩士,畢業院校;中國航空院,職稱:工程師,現就職于:中航工業南京輕型航空動力有限公司,研究方向:發動機結構設計及強度分析。