劉瀟瀟，張國政，董艷琪，魏廣超

(1.哈爾濱理工大學(xué)榮成學(xué)院，山東 威海 264300;2.內(nèi)蒙古工業(yè)大學(xué)(金川校區(qū))，內(nèi)蒙古 呼和浩特 010000)

0 引 言

水下潛器在工作過程中，導(dǎo)管螺旋槳的推力特性成為了其能否完成預(yù)定認(rèn)為的關(guān)鍵。但是對近海潛器的螺旋槳推力特性研究較少，基于此，我們對近海潛器導(dǎo)管螺旋槳的推力特性進行分析計算研究。

首先我們建立了控制方程和雷諾應(yīng)力湍流模型，在得到速度場的前提下，利用Fluent中的半隱式求解器，求出壓力場，進而求得近海海域?qū)Ч苈菪龢耐屏μ匦裕瑸榻：Ｓ蚶脤?dǎo)管螺旋槳推動的潛器控制設(shè)計奠定基礎(chǔ)。

1 選擇湍流模型

螺旋槳工作時，其運動狀況是一種高度復(fù)雜的運動，其中各運動參量隨時間和空間的變化而變化，形成湍流運動。在研究導(dǎo)管螺旋槳的推力特性時，需要首先建立控制方程湍流模型[1]，分析流場中螺旋槳的數(shù)值參數(shù)。

1.1 流動控制方程

由于導(dǎo)管螺旋槳工作于近海海域，因此假定流場中的流體為不可壓縮恒定流體，所以只需考慮質(zhì)量守恒方程和動量守恒方程：

1.2 雷諾應(yīng)力湍流模型

為了簡化算法，選擇雷諾應(yīng)力模型，將N-方程中的速度瞬變量分解為平均量和波動量：

2 控制方程的離散

在上述中，建立的控制方程為一組非線性偏微分方程，只知道流場的速度場，我們無法對控制方程進行求解。所以，需要對控制方程進行離散處理，進而對控制方程進行間接求解。

利用有限元體積法，對計算域進行離散化處理，得到離散方程組，利用插值法，間接得到方程的解。

對于有限元體積法的離散處理[2]，動量方程在笛卡爾坐標(biāo)下表達式為：

式中：φ=1，Γφ為φ的擴散系數(shù)；sφ為廣義源項；Fi為流通量和擴建通量的總統(tǒng)量。

3 構(gòu)造幾何模型

為了對導(dǎo)管螺旋槳在不同進給速度、不同轉(zhuǎn)速條件下進行數(shù)值仿真，需要構(gòu)建導(dǎo)管螺旋槳進行三維建模。

3.1 獲取螺旋槳三維參數(shù)坐標(biāo)

根據(jù)二維工程圖，利用導(dǎo)管螺旋槳切面的葉背值和葉面值進行坐標(biāo)變換[3]，將工程圖中導(dǎo)管螺旋槳的二維參數(shù)轉(zhuǎn)化為曲面型三維參數(shù)坐標(biāo)值，見表 1。

表1 螺旋槳參數(shù)

3.2 導(dǎo)管螺旋槳三維模型

利用Pro/E生成螺旋槳的三維實體模型，將三維模型導(dǎo)入GAMBIT中，利用GAMBIT建模功能完成導(dǎo)管螺旋槳槳轂和導(dǎo)管的三維幾何建模，最后利用GAMBIT中的布爾運算功能將葉片、漿轂以及導(dǎo)管的合并，從而形成導(dǎo)管螺旋槳的三維模型，見圖1。

圖1 導(dǎo)管螺旋槳三維圖 圖2 計算域

4 確定計算域及邊界條件

文中，采用的是Fluent有限元分析中的半陰式求解器，需要設(shè)定計算域(見圖2)，劃分網(wǎng)格和設(shè)定邊界條件。

4.1 設(shè)定計算域

螺旋槳定位為民用推動動力設(shè)備，假設(shè)螺旋槳置于無限大的淺水流域，為了模擬這一流場且減少邊界的影響，把計算域規(guī)定成一個直徑為4D圓柱體，螺旋槳槳軸中心置于離入口4D、離出口5D的圓面圓心處。

4.2 劃分網(wǎng)格

首先在線上布點，生成面網(wǎng)格，最后生成體網(wǎng)格。在靠近螺旋槳的方向上按照1:2對點進行加密。網(wǎng)格類型上，在螺旋槳的周圍區(qū)域設(shè)置為Tet/Hybrid單元TGrid類型，其他的網(wǎng)格采用四面體網(wǎng)格，利用Cooper生成體網(wǎng)格。

4.3 確定邊界條件

在進口邊界設(shè)置為速度進口邊界條件，設(shè)置進口速度為v；出口邊界定義為質(zhì)量流出出口，出口的靜壓力為0;圓柱體側(cè)面設(shè)置為壁面邊界條件，采用單通道計算域的兩側(cè)的矩形切面設(shè)置為周期性邊界；流場內(nèi)部的分割面均設(shè)為interior邊界條件。

圖3 網(wǎng)格劃分圖

5 流場數(shù)值求解

目前，只知道壓力場，對于螺旋槳的數(shù)值模擬計算，我們需要得到速度場。因此利用壓力耦合方程組的半隱式求法，根據(jù)壓力場，見圖4～6，求解出離散的動量方程，求得速度場。進而利用基于壓力的隱士求解器，得到流速云圖，由圖7可知，在螺旋槳槳葉的不同位置，流速并不相同。

圖4 槳葉正面壓力 圖5 槳葉正面壓力

圖6 動壓等壓圖 圖7 流速云圖

6 螺旋槳推力特性計算

在得到了速度場的前提下，接著需要得到近海海域?qū)Ч苈菪龢耐屏μ匦浴?/p>

首先采用雷諾壓力湍流模型通過計算得到的推力及扭矩，根據(jù)推算出推力系數(shù)、扭矩系數(shù)和敞水效率。計算結(jié)果如表2所列。

表2 導(dǎo)管螺旋槳推力表

7 結(jié) 語

導(dǎo)管螺旋槳的工作環(huán)境為不可壓縮恒定流體，建立了質(zhì)量守恒方程和動量方程，建立了雷諾應(yīng)力模型。接著對控制方程緊張了離散處理，構(gòu)造了幾何模型，設(shè)定了計算域和邊界條件，最后，利用壓力耦合方程組的半隱式求法對控制方程進行求解，得到速度場，進而得到近海潛器導(dǎo)管螺旋槳的推力特性。由表2可知，隨著電機額定功率增加，前向推力隨著增加，但是當(dāng)功率為500～600時，后向推力隨功率的增大而減少。此次研究得到了海海域?qū)Ч苈菪龢耐屏μ匦?，為近海海域潛器控制設(shè)計提供了參考。