喬印虎，張春燕，汪濤，陳杰平，緱瑞賓

(安徽科技學(xué)院機(jī)械工程學(xué)院，安徽鳳陽 233100)

發(fā)動機(jī)缸體加工熱分析方法研究

喬印虎，張春燕，汪濤，陳杰平，緱瑞賓

(安徽科技學(xué)院機(jī)械工程學(xué)院，安徽鳳陽 233100)

對發(fā)動機(jī)缸體的數(shù)控加工工藝進(jìn)行了規(guī)劃。為提高加工質(zhì)量，以其中銑削頂平面為例，用有限元法研究加工過程中切削熱對工件變形的影響，進(jìn)行了有限元仿真，根據(jù)仿真結(jié)果改變加工工藝參數(shù)。

發(fā)動機(jī)缸體;加工工藝;加工精度;有限元熱分析

2 工藝規(guī)程設(shè)計與精銑頂平面分析

2.1 工藝路線

依據(jù)粗基準(zhǔn)面的選擇原則和精基準(zhǔn)面的選擇原則對加工階段進(jìn)行劃分，包括粗加工階段、半精加工階段、精加工階段、光整加工階段，工藝路線方案如下:

2.2 精銑頂平面計算

圖1 加工參數(shù)設(shè)置

3 加工過程的熱分析

3.1 熱力學(xué)基本理論

加工過程主要的傳熱方式為熱傳導(dǎo)，用公式(1)表示，其次是熱對流，用公式 (2)表示。熱傳導(dǎo)是兩個接觸物體之間的能量交換或者一個物體內(nèi)部由于溫度梯度而引起的內(nèi)部能量交換。熱對流是物體與周圍介質(zhì)之間發(fā)生的熱交換，周圍介質(zhì) (如空氣或水)的流動帶走物體的熱量，產(chǎn)生熱交換，對流作為表面邊界條件施加［6］。

式中:q為熱流率，k為熱傳導(dǎo)率，A為熱傳導(dǎo)面積，ΔT/L為熱傳導(dǎo)方向的熱梯度。

式中:q為熱流率，h為對流換熱系數(shù)，A為對流表面面積，ΔT=Ts－T∞，Ts為對流表面溫度，T∞為環(huán)境流體溫度。

對體積為V、邊界為S的連續(xù)介質(zhì)，其能量守恒方程為［7］:

式中:vi是速度場，U為給定能量;Q為給定體積熱流;bi為給定體積力;Pi為單位面積上的邊界力;H為邊界S上的單位面積的熱流強(qiáng)度。對體積為V、密度為ρ的連續(xù)介質(zhì)，可建立積分力平衡方程。即:

引入柯西應(yīng)力分量σij，壓力可用柯西應(yīng)力表示為:

式中:ni表示表面S的單位法線方向。將力平衡方程引入能量守恒方程 (3)，可得到熱－機(jī)耦合的能量守恒方程。

依據(jù)虛功原理，建立結(jié)構(gòu)位移ui的方程式:

假設(shè)忽略慣性項的影響，則式 (7)右端第二項可去掉。并且假設(shè)物體的能量方程和力平衡方程都是建立在當(dāng)前的構(gòu)型上，則可用弱耦合的增量非線性有限元法處理熱－機(jī)耦合問題。

3.2 熱分析有限元方法

3.2.1 熱分析塊體單元

在ALGOR軟件中，將缸體按有限元方法劃分網(wǎng)格，如圖2所示，形成8節(jié)點單元并檢驗?zāi)Ｐ?，每個節(jié)點有一個溫度自由度，考慮與溫度相關(guān)的各向同性、各向異性材料性能，定義包含中節(jié)點，包含中節(jié)點的六面體單元為20節(jié)點六面體單元，進(jìn)入Result環(huán)境確定適當(dāng)?shù)墓?jié)點號。

圖2 待加工缸體有限元模型

3.2.2 熱流率計算方法

選擇與“Heat Flow Calculation”相對應(yīng)的下拉式菜單中的選項;“Projected at Centroid”:零件的熱通量由節(jié)點溫度通過傅里葉法計算得到; “Nonlinear Based on BC”，承受對流或輻射載荷的零件外表面的熱通量，通過對流或輻射邊界條件的輸入?yún)?shù)以及節(jié)點溫度得到，對內(nèi)表面的熱流率沒有影響;“Linear Based on BC”:與“Nonlinear Based on BC”計算方法類似，只是承受輻射載荷作用表面的熱通量將被線性化，為得到輻射或?qū)α鬟吔缟蠝?zhǔn)確的熱通量輸出，可以利用“Nonlinear Based on BC”或“Linear Based on BC”選項。

3.2.3 分析步驟

銑刀參數(shù)見表1。

表1 銑刀參數(shù)

打開模型。在 Pro/E中建模，將模型保存為*.igs格式。啟動ALGOR，分析類型選擇穩(wěn)態(tài)熱傳導(dǎo)，打開缸體模型:(1)劃分網(wǎng)格;(2)修改單位系統(tǒng)，Unit Systm選擇自定義，修改 Enrgy單位為Btu，確定;(3)定義材料。選擇材料Iron，Gray Cast ASTM A-48 Grade 25;(4)指定溫度邊界，指定溫度，輸入Magnitude=1 278;(5)施加對流界面。選擇所有暴露在空氣中的面，右擊鼠標(biāo)，選擇添加→表面對流載荷，點擊Calculate，將通過流體參數(shù)計算對流換熱系數(shù)，在“總體”面板中輸入流速17.6，在“流體屬性”面板中，輸入如下參數(shù):質(zhì)量: 1.167×10－7;動力黏滯系數(shù):2.56×10－9;導(dǎo)熱系數(shù):3.5×10－7;比熱:92.8。再設(shè)置對流換熱系數(shù)為1.829×10－6。設(shè)置約束條件、施加載荷并設(shè)置求解參數(shù)，保存并由Perform Analysis求解。

3.3 結(jié)果分析

通過圖3、4可以看出:粗加工與精加工時，熱應(yīng)力結(jié)果不同。除了可以輸出熱應(yīng)力結(jié)果外，還能輸出切削溫度、變形量的值。正如經(jīng)典理論所描述，最高溫度出現(xiàn)在距離刀尖一定距離的位置。利用ALGOR后處理技術(shù)得到粗、精加工的應(yīng)力分布，精加工時熱應(yīng)力較小，這與實際加工是相符的，證明了該有限元模型的有效性和合理性。

圖3 粗加工時熱應(yīng)力

圖4 精加工時熱應(yīng)力

4 結(jié)論

采用有限元方法對發(fā)動機(jī)缸體切削過程進(jìn)行了模擬，并分析了切削時應(yīng)力的分布情況，所得結(jié)果與傳統(tǒng)切削理論相符。其模擬結(jié)果可以指導(dǎo)切削條件的合理確定、刀具幾何參數(shù)及切削用量的合理選擇，從而提高刀具壽命、降低生產(chǎn)成本。金屬切削加工過程涉及到彈、塑性變形、斷裂、劇烈摩擦、相變以及刀具磨損等，這些現(xiàn)象是相互影響、相互關(guān)聯(lián)的。盡管切削加工的有限元模擬與實際相比還有不足，但是該技術(shù)已成為今后發(fā)展的主要方向。

參考文獻(xiàn):

【1】王良.我國航空發(fā)動機(jī)制造技術(shù)的現(xiàn)狀與挑戰(zhàn)［J］.航空制造技術(shù)，2008，25(S1):32－37.

【2】唐志濤，劉戰(zhàn)強(qiáng)，艾興，等.金屬切削加工熱彈塑性大變形有限元理論及關(guān)鍵技術(shù)研究［J］.中國機(jī)械工程，2007，18(6):746－751.

【3】陸郁，廖最巧.應(yīng)用有限元分析優(yōu)化缸體曲軸孔精鏜加工［J］.裝備制造技術(shù)，2011(10):141－144.

【4】馮勇，汪木蘭，王保升.高速切削熱及溫度預(yù)測研究進(jìn)展［J］.機(jī)械設(shè)計與制造，2012(5):261－263.

【5】李洪.機(jī)械加工工藝手冊［M］.北京:北京出版社1999.

【6】寇曉東.ALGOR結(jié)構(gòu)分析高級教程［M］.北京:清華大學(xué)出版社，2008.

【7】李目，陳蔚芳.銑削加工中工件變形仿真預(yù)測方法研究［J］.機(jī)械制造，2010，48(1):51－54.

Study on Thermal Analysis Method for Engine Block Machining

QIAO Yinhu，ZHANG Chunyan，WANG Tao，CHEN Jieping，GOU Ruibin
(Anhui Science and Technology University，F(xiàn)engyang Anhui 233100，China)

The CNC machining process of the engine block was designed.Taking the milling top plane as example，the finite element simulation was used to study the thermal deformation of the engine block.Then the processing parameters were changed to reduce the thermal impact on the workpiece to improve processing quality.

Engine block;Machining process;Machining precision;Finite element thermal analysis

［5］可以確定:底面的加工方案為銑平面:粗銑 (IT9)，粗糙度為6.3 μm，一般不淬硬的平面，精銑的粗糙度可以較小。

TH061

A

1001－3881(2014)8－009－3

10.3969/j.issn.1001－3881.2014.08.004

2013－02－20

國家自然科學(xué)基金項目 (61164012);安徽科技學(xué)院重點建設(shè)學(xué)科 (AKXK20102－5);安徽科技學(xué)院青年基金(ZRC2013337);安徽省優(yōu)秀青年人才基金重點項目 (2013SQRL062ZD)

喬印虎 (1979—)，男，在讀博士，講師，主要研究方向為數(shù)控與檢測、汽車、風(fēng)電葉片。E－mail:qyh7926@ 163.com。

發(fā)動機(jī)作為汽車的五大部件之一，其制造在汽車工業(yè)中具有舉足輕重的地位。同樣，在新一代航空發(fā)動機(jī)性能的提高中，制造技術(shù)與材料的貢獻(xiàn)率為50%～70%;在發(fā)動機(jī)減重的貢獻(xiàn)率中，制造技術(shù)和材料的貢獻(xiàn)率占70%～80%，這充分表明材料與工藝是制約新一代航空發(fā)動機(jī)發(fā)展的關(guān)鍵技術(shù)［1］。切削加工過程是一個復(fù)雜的強(qiáng)熱力耦合的動態(tài)物理過程，借助數(shù)值仿真技術(shù)，特別是有限元法，可有效預(yù)測熱力耦合作用下的切屑形態(tài)，工件、刀具以及切屑中的非均勻應(yīng)力場、應(yīng)變場和溫度場，以及切削力、刀具磨損，工件加工表層殘余應(yīng)力等一些物理量的分布情況，進(jìn)行工藝參數(shù)優(yōu)化。唐志濤等［2］導(dǎo)出了熱彈塑性大變形耦合控制方程，建立了正交切削加工鋁合金7050T7451有限元模型，對切屑形態(tài)、切削力、切削溫度以及應(yīng)力場和應(yīng)變場等物理量的分布進(jìn)行了有效預(yù)測。陸郁等人［3］介紹了有限元軟件 Advantedge在研究機(jī)加工工藝參數(shù)對切削性能影響方面的運(yùn)用，通過比較仿真計算結(jié)果，從中得出各種切削加工性能指標(biāo)，從而優(yōu)化加工參數(shù)，并通過實例——缸體曲軸孔的精鏜加工，闡述了仿真方法。數(shù)值仿真技術(shù)已成為國內(nèi)外眾多學(xué)者研究切削熱的主要方法，理論不斷發(fā)展深入［4］。

1 發(fā)動機(jī)缸體零件分析

有3組加工:底面、側(cè)面3個孔，頂面的4個孔以及左視圖上的2個孔。以底面為主要加工的表面，有底面的銑加工;另一組加工是側(cè)面的曲軸孔;以頂面為主要加工的表面，有頂面的銑加工;以前端面為主要加工表面，有斜面的銑加工;以后端面為主要加工表面，有端面的銑加工;以曲軸聯(lián)接面為主要加工表面，有端面的銑加工。