張宗熠，高殿榮，許森浩，孫亞楠，張 琰

（1.燕山大學(xué) 機械工程學(xué)院，河北 秦皇島 066004；2.中國航空工業(yè)集團公司西安飛行自動控制研究所，西安 710000）

機械裝備的現(xiàn)代化發(fā)展要求液壓系統(tǒng)具有更高的機電液一體化程度［1］。近些年來，在液壓系統(tǒng)動力源中，通過將電動機、聯(lián)軸器和液壓泵進行高度集成，形成了一種高度集成化的液壓動力單元——電機泵。電機泵具有結(jié)構(gòu)緊湊、功率密度大、噪聲低、效率高、控制特性好等一系列優(yōu)點，展現(xiàn)了良好的發(fā)展?jié)摿蛻?yīng)用前景［2-5］。溫升是電機泵在設(shè)計時需要考慮的主要性能指標(biāo)之一，對于單位時間發(fā)熱量大的高速電機泵來說尤為重要。過高的溫升加速了定子繞組絕緣層的老化，導(dǎo)致永磁材料退磁，降低了工作介質(zhì)的黏度，嚴(yán)重影響了高速電機泵的工作性能和工作可靠性，甚至直接引起高速電機泵的故障［6］。因此，對高速電機泵的溫升和溫度分布特性進行準(zhǔn)確計算是保證電機泵正常工作的基礎(chǔ)，同時為高速電機泵的優(yōu)化設(shè)計提供了必要的參考。

電機和變壓器等電氣元件和電機泵具有相似的發(fā)熱機理和傳熱結(jié)構(gòu)，目前關(guān)于電機和變壓器的溫度場計算所采用的數(shù)值計算方法主要有等效熱網(wǎng)絡(luò)法和有限元法［7-12］。等效熱網(wǎng)絡(luò)法是應(yīng)用圖論原理，通過網(wǎng)絡(luò)的拓?fù)浣Y(jié)構(gòu)進行熱場分析的一種方法。在對高速電機泵的溫升進行計算時，采用有限元計算方法可以對電機泵的整體溫度場進行計算，不僅可以獲得關(guān)鍵組件的平均溫升，還可以獲得最大溫度和過熱點的位置，且具有較高的計算精度。張鳳閣等［13］分別采用熱網(wǎng)絡(luò)法和有限元法對一種新型磁障轉(zhuǎn)子無刷雙饋電機的溫度分布情況進行了計算，并通過實驗結(jié)果驗證了兩種數(shù)值計算方法的正確性。黃毅等［14］對大功率高壓開關(guān)電源變壓器的損耗分析及散熱控制進行了研究，計算了不同負(fù)載條件下的溫升情況及散熱所需風(fēng)量，通過ANSYS軟件對溫度場進行分析。SHINGO［15］采用有限元法對由鐵損和銅損引起的開關(guān)磁阻電機的穩(wěn)態(tài)溫升進行了研究。研究表明：與測量值相比，考慮開關(guān)磁阻電機中內(nèi)部空氣狀況可以更精確地計算電機溫升。

有關(guān)學(xué)者對電機泵的損耗及溫度分析也做了一定的研究。冀宏等［16］對一種非浸油結(jié)構(gòu)的電機葉片泵內(nèi)的溫度場進行了數(shù)值模擬，分析了主要部件的溫度分布情況。許丹丹［17］對一種液壓電機葉片泵的溫度場進行了計算。計算結(jié)果表明：相比風(fēng)冷方式，油冷方式能夠顯著降低電機泵的溫升。綜上所述，國內(nèi)外專家學(xué)者對不同電機和電機泵的溫度場進行了一定研究，但是未有針對高速柱塞式液壓電機泵溫度場的研究。

為此，本文首先分析了高速電機泵發(fā)熱及傳熱機制，建立了高速電機泵溫度場的流固耦合傳熱計算方法，應(yīng)用此方法對一種高速電機泵的溫度場進行了數(shù)值計算，分析了側(cè)流道面積、側(cè)流道個數(shù)、側(cè)流道與泵內(nèi)表面距離分布的影響溫度場的影響，并通過正交試驗法對不同的流道形式進行了優(yōu)化設(shè)計，計算方法和結(jié)果對類似結(jié)構(gòu)的高速電機泵的溫升計算具有重要的參考意義。

1 數(shù)值計算方法及計算模型

1.1 發(fā)熱和傳熱機制分析

高速電機泵在運行過程中，由于電磁能量轉(zhuǎn)化產(chǎn)生的損耗和摩擦損耗會在系統(tǒng)中以熱量的形式產(chǎn)生。發(fā)熱損耗主要包含了電磁損耗、各個摩擦副的摩擦損耗和油液黏性摩擦損耗。這些熱量一方面被工作介質(zhì)吸收，帶入液壓系統(tǒng)當(dāng)中；另一方面通過電機泵殼體與周圍空氣通過對流換熱的方式傳送出去。在電機泵內(nèi)部還會以熱傳導(dǎo)的形式在各部件之間進行熱量傳遞，由于發(fā)熱熱源和散熱條件的雙重作用，電機泵在穩(wěn)態(tài)運行時達(dá)到一個穩(wěn)態(tài)的熱平衡狀態(tài)。因此，高速電機泵在工作工程中的溫度傳遞過程是一個涉及電磁場、流體場和溫度場的多物理場耦合問題，高速電機泵整機生熱與散熱路徑如圖1所示。

1.2 數(shù)值計算方法

為了求解高速電機泵的穩(wěn)態(tài)溫升和溫度分布特性，首先需要建立涉及電磁場、流體場和溫度場的數(shù)值算法。電磁場計算的目的在求解電機泵運行當(dāng)中電磁能量轉(zhuǎn)化所產(chǎn)生的熱量，忽略整個計算場中流體場和溫度場對電磁損耗的影響，則電磁場與流體場和溫度場之間只有單向的耦合作用。因此，可先對電磁場進行求解，再將損耗作為熱源加載于流體場和溫度場之上。

在流體場和溫度場的求解問題上，同時涉及到流體計算域和固體計算域。流體計算域和固體計算域之間存在著熱量的交換。流體計算域的流體流動及傳熱過程受物理基本定律的支配，這些物理基本定律包括質(zhì)量守恒定律、動量守恒定律和能量守恒定律。針對不同的物理學(xué)定律，可以對其用不同的偏微分方程進行描述。

穩(wěn)定狀態(tài)下，不可壓縮流體質(zhì)量守恒方程和動量守恒方程可分別表示為：

式中：u、v和w分別為x、y和z方向的速度分量；p為流體壓力；ρ為流體密度；μ為流體動力黏度；Su、Sv和Sw為廣義源項。

能量守恒方程可以表示為：

式中：T為流體溫度；k為導(dǎo)熱系數(shù)；cp為比熱容；ST為能量耗散項。

在流體能量守恒方程中，如果將對流項和廣義源項中的速度場量去掉，則流體能量守恒方程可退化為固體導(dǎo)熱微分方程。因此，對于流體域和固體域，方程具有相同的形式，可應(yīng)用同一網(wǎng)格系統(tǒng)和求解器進行計算。

1.3 冷卻結(jié)構(gòu)及計算域模型

高速電機泵整機模型如圖2所示。取消了柱塞泵的殼體，將柱塞泵直接安裝在電機轉(zhuǎn)子內(nèi)，柱塞缸體與電機轉(zhuǎn)子采用花鍵聯(lián)接。采用油道內(nèi)循環(huán)的布置方式，當(dāng)油液由入口進入電機泵內(nèi)腔時，先經(jīng)過進油口后進入電機泵前端蓋腔內(nèi)后，油液可通過電機泵外殼部分的流道、定子槽孔流道以及定子與轉(zhuǎn)子之間的油隙這3部分進入電機泵后端蓋腔內(nèi)，在經(jīng)過核心部分的流道進入泵的吸油口，在柱塞泵內(nèi)隨柱塞腔的旋轉(zhuǎn)后到出油口后排到電機泵外。采用此種結(jié)構(gòu)可將溫度較低的液壓油當(dāng)作冷卻介質(zhì)充分地在電機泵內(nèi)循環(huán)流動，從而能更高效地降低電磁場產(chǎn)生的熱量、泵內(nèi)軸承摩擦產(chǎn)生的熱量以及柱塞泵摩擦做功所產(chǎn)生的熱量，從高速電機泵中提取的流場計算模型如圖3所示。

1.4 邊界條件及計算假設(shè)

入口采用速度入口邊界條件，入口流量與液壓泵的額定流量相同，出口采用壓力出口，出口壓力為0 MPa。外表面采用對流換熱邊界條件，流體域與固體域的交界面采用耦合邊界條件。

由于高速電機泵內(nèi)結(jié)構(gòu)較為復(fù)雜，為了合理簡化求解過程，做出以下基本假設(shè)：由于液壓泵組件內(nèi)元件眾多，暫忽略液壓泵組內(nèi)流場情況，只對電機定轉(zhuǎn)子部分的溫度場進行計算。

忽略電機泵在裝配過程中出現(xiàn)的縫隙結(jié)構(gòu)，將電機泵看作一個整體進行計算分析。本文中暫不涉及柱塞泵內(nèi)流場的分析，故此部分將不進行具體的計算。

2 計算結(jié)果分析

2.1 溫度場分析

首先使用初始參數(shù)進行數(shù)值計算，計算參數(shù)及樣式如表1所示。

表1 流道計算參數(shù)及樣式

將磁場有限元計算出的熱量的分布數(shù)據(jù)導(dǎo)入Fluent軟件，在進行條件設(shè)定后對其進行求解。

如圖4所示，油液進入電機泵腔內(nèi)沿著流道向泵的吸油口流入，在側(cè)流道處將定子繞組所產(chǎn)生的熱帶走。在半徑方向上，隨半徑增大溫度場在越過側(cè)面流道后溫度有明顯下降，約為8℃，這說明在電機泵殼增加側(cè)流道是有效的；軸向上，泵內(nèi)后腔體內(nèi)油液的溫度高于前端蓋腔內(nèi)的油液溫度，溫差為5℃；定子槽的溫度場分布偏向后腔體部分，這是由于在油液由前端流至后端的過程中油液的溫度會逐步增加，前端溫度較低、后端溫度較高，這樣會使軸向上定子繞組槽所散發(fā)出的溫度分布場不對稱。

2.2 流道參數(shù)對溫度場分布的影響分析

為了優(yōu)化電機泵的散熱減少由溫度導(dǎo)致的熱應(yīng)力與熱變形對整體結(jié)構(gòu)外形以及內(nèi)部配合精度的影響，在其結(jié)構(gòu)設(shè)計上：一是通過合理的流道結(jié)構(gòu)參數(shù)設(shè)計使其最高溫度下降；二是通過不同的結(jié)構(gòu)形式使內(nèi)部的溫度場的分布更加均勻，減小溫度集中區(qū)域的密度。通過側(cè)流道截面積、側(cè)流道個數(shù)、流道距內(nèi)壁面距離這3個因素來分析（各包含3個水平值）溫度場分布以及溫度值的變化。如表2所示，將采用如下幾種組合形式來對其分析。

表2 計算流道樣式參數(shù)

2.2.1 側(cè)流道面積變化對溫度場分布的影響

其他參數(shù)保持為初始計算參數(shù)，只改變側(cè)流道截面積，分別為2 904、2 934、2 964 mm23個數(shù)值來觀察溫度場的變化，結(jié)果如圖5所示。

結(jié)合圖4、5中3種不同溫度分布場可以發(fā)現(xiàn)：電機泵前腔體內(nèi)的溫度比后腔體內(nèi)的溫度平均高3～5℃，隨截面積增大，最高溫度不是線性變化；在沿軸向上，側(cè)流道截面積增大溫度的分布區(qū)域更加均勻，且溫度變化梯度也越穩(wěn)定。通過這一現(xiàn)象可以得出：隨著冷卻流道截面積的增加，它能使電機泵產(chǎn)生的熱量更多地接觸到冷卻流道，并且在降溫的過程中溫度場分布更加均勻；在不改變其他條件時，由于電機泵入口的流量是一定的，增加側(cè)流道的截面積相當(dāng)于減小單位時間內(nèi)通過單個流道油液的流量，這種現(xiàn)象會使油液在單位時間內(nèi)“運送”走的熱量減??；側(cè)流道截面積增大會使側(cè)流道更多地與高溫部分接觸。這一方面會使油液的整體溫度更高；另一方面也會使泵殼上的最高溫度略微降低。所以增加側(cè)流道截面積可以帶走更多的熱量，也會使溫度場的分布更加均勻。截面積增大會導(dǎo)致最高溫度有所降低，同時也會使泵內(nèi)的平均溫度達(dá)到最低，隨著截面積的增大，降溫效果更顯著。

由圖6可以看出：軸向溫度分布上，在第1、2種水平下高溫區(qū)域分布是關(guān)于中間部分對稱的，在第3種水平下是相對靠向后腔位置的；從兩端溫度數(shù)值上來看，在第1、2種水平下前端部分的溫度是高于后端的溫度，第3種水平下前端溫度較低，后端溫度較高。這表明：運用第1、2種水平的結(jié)構(gòu)設(shè)計會導(dǎo)致油液在側(cè)流道所能帶走的熱量減少，這樣會使已經(jīng)流出側(cè)流道的部分油液在進入后腔內(nèi)時，油液溫度升高得不夠明顯，而在第3種結(jié)構(gòu)設(shè)計下冷卻油液則能夠更加充分地吸收熱量。綜合來看，側(cè)流道截面積的選擇與流體的運動狀態(tài)密切相關(guān)，截面積增大使得接觸到壁面的油液更多，在這種因素下相對較低的流速會使散熱效果明顯。

2.2.2 側(cè)流道個數(shù)對溫度場分布的影響

保持其他各項參數(shù)不變，側(cè)流道的個數(shù)分別設(shè)置為6、9、12，對3種不同的情況進行模擬，得出如下的結(jié)果：如圖4與圖7，隨側(cè)流道個數(shù)增大泵內(nèi)整體最高溫度增大。當(dāng)流道數(shù)為9時，溫度高于63℃的范圍處于中心處且有向泵殼表面擴散的趨勢；當(dāng)流道數(shù)為12時，這種趨勢反而減少。這主要因為：在保證側(cè)流道截面積相同時改變流道個數(shù)相當(dāng)于改變單個側(cè)流道的截面積。在流道數(shù)為12時，相當(dāng)于用更多的側(cè)流道阻擋熱量的發(fā)散，增加了單位面積上冷卻油道的占比。這樣會導(dǎo)致油液在通過截面積小的孔時以更快的速度流出，這樣油液與溫度較高的壁面換熱的時間減少，反而降低了換熱的效率。

如圖8所示，在3種水平下泵殼外表面的最高溫度基本都保持在相同的數(shù)值。不同的是，在流道個數(shù)分別為6、9時，溫度分布是基本對稱的，但是在側(cè)流道數(shù)為9時，兩端的溫度略微高于側(cè)流道個數(shù)為6時的狀態(tài)，這說明在側(cè)流道數(shù)為9時的散熱的效果反而略差，但在這兩種水平下，泵殼后端的溫度均比前端的溫度要低，這也說明在這兩種水平下，側(cè)流道的換熱效果對整體的溫度分布而言有一定的效果。但當(dāng)流道數(shù)為12時，溫度最高的區(qū)域的分布是明顯向后偏移的，且后端的溫度高于前端的溫度，并且此時前端的溫度為3種水平的最低，后端的溫度為3種水平的最高，這種現(xiàn)象說明了在周向上流道數(shù)為12時，散熱的效果主要體現(xiàn)在了最高溫度的分布上。與改變測流道截面積的大小相比，側(cè)流道個數(shù)的影響主要反映在溫度場軸向上的分布，但對于降低整體的最高溫度的效果不顯著。這種現(xiàn)象說明：在此尺寸下有限地增加側(cè)流道的個數(shù)能使整體的溫度分布朝著有利的方向改善，但對控制最高溫度值無明顯效果。

2.2.3 側(cè)流道與泵內(nèi)表面距離對溫度場分布的影響

在前1小結(jié)中發(fā)現(xiàn)：當(dāng)側(cè)流道個數(shù)為12時效果優(yōu)于流道數(shù)為9的狀態(tài)，這說明在周向上增大單位面積上冷卻流道的占比是有效的，在半徑方向上設(shè)置不同的流道距離核心的距離，可以使溫度較高的泵殼以不同的溫度來接觸冷卻油道。因此，分別設(shè)計側(cè)流道底端距泵內(nèi)壁面7、10、13 mm觀察溫度場的分布狀態(tài)，如圖9、10所示。

結(jié)合圖4與圖10，溫度場的最高溫度并不總是隨著側(cè)流道距離內(nèi)壁面的增大有規(guī)律性變化。與水平2相比水平1、3的最高溫度略低，3種水平下其溫度分布場均偏向于后腔部分。這表明：適當(dāng)改變側(cè)流道距離內(nèi)表面的距離，會使溫度場分布向后腔偏移；但當(dāng)側(cè)流道與內(nèi)表面的距離減小時，流道在截面的形狀會顯得細(xì)長，并且在對比3種水平下的數(shù)值結(jié)果時發(fā)現(xiàn)：水平2的最高溫度能保持在63.5℃附近，但水平1的最高溫度維持在60℃附近，水平3的最高溫度在61℃附近。這種現(xiàn)象表明了側(cè)流道距離內(nèi)壁面這項因素是不太顯著的。但此項只是在微小尺寸上（＜13 mm）分析得出的結(jié)論，在對整體結(jié)構(gòu)繼續(xù)擴大后發(fā)現(xiàn)：當(dāng)徑向尺寸變化較大時（＞15 mm），其結(jié)果就很顯著，因受限于電機泵整體的尺寸限制，泵殼的直徑不能無限增大，故在此條件下，側(cè)流道內(nèi)的溫度場的分布情況并不理想。

觀察圖10可以發(fā)現(xiàn)：不同側(cè)流道距內(nèi)表面的距離對于泵殼外表面溫度的影響不明顯。造成這種現(xiàn)象的原因如下：在小幅度改變不同的距離，對于流道在周向上的分布情況影響幾乎可以忽略，故對于泵殼表面最高溫度的影響收效甚微。

3 流道參數(shù)的優(yōu)化設(shè)計

通過設(shè)計更為全面的側(cè)流道的幾何模型，應(yīng)用正交試驗的方法來進一步研究各流道參數(shù)對溫度場的影響，并探究其中規(guī)律。

3.1 側(cè)流道參數(shù)的正交試驗設(shè)計

在設(shè)計正交試驗時首先對要解決的問題拆分成不同的因素與水平的組合。為了更全面地考慮側(cè)流道幾何參數(shù)對溫度場分布的影響，設(shè)計了側(cè)流道截面形狀A(yù)、側(cè)流道軸向樣式B、側(cè)流道截面面積C、側(cè)流道個數(shù)D、側(cè)流道距內(nèi)壁面的距離E以及側(cè)流道進出口直徑F 6個因素，參見表3?？紤]在徑向上單位面積內(nèi)冷卻流道的占有率對溫度場分布的影響，設(shè)計了橢圓形、扇形和梯形3種不同的側(cè)流道的截面形狀，并在標(biāo)準(zhǔn)的基礎(chǔ)模型上分別減少與增加10%的側(cè)流道在相同的半徑向上的面積占比，流道模型的截面形狀如圖11所示。螺旋式流道具有強化對流換熱的效果，設(shè)計了在軸向的前端與后端有一定的角度差的流道，旋轉(zhuǎn)角度分別為0°、15°與30°。目的在于在保證相同截面積的前提下通過變相增加過流的長度從而能控制通過側(cè)流道的時間來觀察對溫度場分布的影響，流道模型軸向如圖12所示。

3.2 正交試驗結(jié)果分析

結(jié)合以上因素與水平的個數(shù)選擇6因素3水平的正交表，但現(xiàn)有正交表系中并沒有這種，根據(jù)原則選取L18（37）正交表。對各不同形狀參數(shù)的模型進行幾何建模劃分網(wǎng)格后運用Fluent軟件進行有限元分析。最終結(jié)果如表4所示，選擇的指標(biāo)為最高溫度，這項參數(shù)主要反映出冷卻流道的效果是否明顯。

基于以上情況對側(cè)流道幾何參數(shù)對泵殼內(nèi)溫度場的分布影響情況的因素以及所對應(yīng)的水平表如表3所示。

在對以上正交試驗的結(jié)果分析完成后就可以對其進行極差分析，極差能直接反映出結(jié)果的離散程度，也能反映出個水平間的重要程度。

表3 側(cè)流道幾何參數(shù)的因素與水平

表4 側(cè)流道幾何參數(shù)的實驗方案及其計算結(jié)果

在表5中K1×A的單元格中的數(shù)據(jù)為：在因素側(cè)流道截面形狀都為水平1情況下的6次試驗的結(jié)果之和，而K1×A單元格為上述值的6次試驗的平均值，其余單元格內(nèi)的數(shù)值計算方法均類比。而極差這一項為每個因素的均值中最大項與最小項的差值，其余類比計算。由于極差直接反映出本試驗中單項因素的權(quán)重，能明顯地反映出不同因素間的重要程度。根據(jù)此可以輕易判斷出每個因素的權(quán)重即F＞C＞D＞E＞A＞B，但在極差值中，由圖13可看出因素A、C、D、E、F這5項明顯比B項重要。

表5 極差分析結(jié)果

前面幾組情況反映出油液在流道流經(jīng)時的長度會影響溫度場的分布，但從此次正交試驗的結(jié)果來看，其影響并不強烈。故在此基礎(chǔ)上又增設(shè)了幾種狀況，分別為：在原始模型中側(cè)流道在軸向長度上按0.2倍依次增加4次，因要保持軸向長度相同，故采取改變進出口的螺旋角來變相改變軸向長度。通過一系列的有限元分析，最終發(fā)現(xiàn)：在側(cè)流道長度上相差20%以上（為保持軸向長度相同，通過改變螺旋角度來變相增加過流長度）時，且側(cè)流道的基礎(chǔ)的長度大于0.4 m時側(cè)流道的軸向樣式這項因素有明顯的效果，但由于結(jié)構(gòu)參數(shù)的限制，尺寸不能無限改變，故此項系數(shù)在這種結(jié)構(gòu)下效果并不顯著。除此之外，由于因素C、F都能直接對流體在側(cè)流道內(nèi)的流速產(chǎn)生強烈的影響，這也驗證了此前的結(jié)論。

根據(jù)正交試驗的結(jié)果，選取最終的最優(yōu)組合，再次進行幾何建模與有限元分析，最終結(jié)果顯示：在最優(yōu)組合下，泵內(nèi)最高溫度維持在61℃附近，且平均溫度降低至43℃附近，泵殼外表面的溫度場明顯偏向后端部分，且外表面的最高溫度降低至52℃附近。

4 結(jié)論

1）無論在何種參數(shù)下，泵殼上溫度場的分布基本都是關(guān)于中間部分對稱分布的，溫度高的部分主要在中間部分側(cè)流道的內(nèi)側(cè)，溫度低的部分集中在兩端，由于受到液壓油的影響，降溫的效果明顯。

2）泵殼上溫度分布主要受側(cè)流道內(nèi)油液流動狀態(tài)影響較大，當(dāng)其流速變化較大時（＞7%），溫度場的分布會有明顯的變化，主要體現(xiàn)在：有更低的最高溫度，整體溫度場分布向后腔體偏移，后腔體內(nèi)油液溫度會明顯高于前腔體的溫度，這樣能有力地說明側(cè)流道對于整體溫度場的控制是有作用的。

3）側(cè)流道進出口直徑、側(cè)流道截面積與側(cè)流道的個數(shù)對于溫度場的分布影響最大，這3個因素本質(zhì)上也是改變側(cè)流道內(nèi)油液的流速，這也證明擁有更高的流速對于內(nèi)部溫度場的控制有積極的作用，但在實際情況下油液流速與流量不能無限增大，只能在保證流速的前提下盡量優(yōu)化內(nèi)部結(jié)構(gòu)，使油液更加充分地接觸到溫度高的區(qū)域。