張亮，紀(jì)勝軍，2，李樹(shù)珍，沈江潔，石磊

(1.河北科技師范學(xué)院機(jī)電工程學(xué)院，河北秦皇島 066004；2.山東華宇工學(xué)院機(jī)械工程學(xué)院，山東德州 253000；3.河北科技師范學(xué)院海洋資源與環(huán)境學(xué)院，河北秦皇島 066004)

0 前言

水動(dòng)力特性是水下機(jī)器人設(shè)計(jì)和運(yùn)動(dòng)控制的基礎(chǔ)[1]。由于水下機(jī)器人運(yùn)動(dòng)時(shí)周圍流場(chǎng)復(fù)雜，水動(dòng)力特性分析困難，需要通過(guò)試驗(yàn)獲取各項(xiàng)水動(dòng)力系數(shù)[2]。國(guó)內(nèi)的水下機(jī)器人測(cè)試系統(tǒng)一般較為大型[3]，如上海交通大學(xué)的多功能拖曳水池整體尺寸長(zhǎng)300 m、寬16 m、深7.5 m，最大拖曳速度10 m/s[4]；中國(guó)船舶重工集團(tuán)公司第七○二研究所擁有的操縱性旋臂水池，直徑48 m、深4.5 m，最大角速度1 rad/s，臂端速度12 m/s，可進(jìn)行船舶、潛水器等的操縱性理論及試驗(yàn)測(cè)試，預(yù)報(bào)各型水中運(yùn)動(dòng)體的操縱性能。而國(guó)外的測(cè)試系統(tǒng)更具有針對(duì)性，專業(yè)性也更高，如德國(guó)宇航院(DLR)研發(fā)的小水洞，試驗(yàn)段尺寸0.26 m×0.33 m，長(zhǎng)1.26 m，試驗(yàn)速度0.5 m/s[5]；ZAHARI M等為研究渦激振動(dòng)是由湍流旋渦脫落現(xiàn)象引起的，設(shè)計(jì)了一個(gè)為了教學(xué)目的的低成本水洞[6]。

本文作者針對(duì)低速水下機(jī)器人模型試驗(yàn)測(cè)取水動(dòng)力特性的需求，基于相似原理[7]，設(shè)計(jì)一套用于測(cè)量小型低速水下機(jī)器人水動(dòng)力系數(shù)的試驗(yàn)系統(tǒng)——室內(nèi)小型低速開(kāi)口式直流水洞。

1 水洞設(shè)計(jì)

針對(duì)小型水下機(jī)器人模型設(shè)計(jì)的室內(nèi)小型低速開(kāi)口式直流水洞如圖1所示，主要包括水洞整流段(過(guò)渡段二1、擴(kuò)散段2、蜂窩器3、穩(wěn)定段4、收縮段5、過(guò)渡段一6)、水箱7、閥8、水泵9、水管10、法蘭11。

圖1 水洞CAD模型

水洞大小主要是由試驗(yàn)段決定的，試驗(yàn)?zāi)Ｐ痛笮Q定了試驗(yàn)段的大小[8]。本文作者借鑒劉彥偉和劉瑩[9]對(duì)有邊界和無(wú)邊界場(chǎng)的模擬分析對(duì)比，得出阻塞比不得超過(guò)25%，否則結(jié)果沒(méi)有實(shí)際意義。借鑒英國(guó)NPL水洞的設(shè)計(jì)，試驗(yàn)段的長(zhǎng)度通常是模型長(zhǎng)度的3倍，以保證模型在更均勻的流場(chǎng)中[10]。

根據(jù)以上要求和模型的尺寸比例大小，以100 mm×100 mm×34 mm機(jī)器人模型為基準(zhǔn)，水洞試驗(yàn)段穩(wěn)流區(qū)域達(dá)到300 mm×200 mm×68 mm即可。

2 水洞動(dòng)力學(xué)仿真試驗(yàn)

2.1 水洞模擬條件設(shè)定

將SolidWorks建立的水洞三維模型導(dǎo)入Space claim軟件，仿真條件設(shè)定如下：

(1)求解域設(shè)置

在Fluent模塊中，檢查網(wǎng)格，保證沒(méi)有負(fù)體積網(wǎng)格，默認(rèn)國(guó)際單位制。

(2)物理模型設(shè)置

選擇標(biāo)準(zhǔn)κ-ε模型，設(shè)置計(jì)算域材料為水，選擇速度入口，過(guò)渡段二入口為模型模擬入口，D=80 mm，基于流量相等原則，入口速度為2.6 m/s。

(3)求解方案設(shè)置

選擇分離求解程序，Simple算法用于壓力和速度的耦合。添加用于收斂檢查的報(bào)告定義——?dú)埐钋€和阻力曲線，曲線收斂才說(shuō)明數(shù)值可信。選擇初始化，并確保初始化值不能大于10-6，最后進(jìn)行求解計(jì)算[11]。

2.2 水動(dòng)力學(xué)仿真結(jié)果分析

通過(guò)模擬計(jì)算，得到如圖2所示流體的運(yùn)動(dòng)狀態(tài)云圖。

圖2 水洞流體運(yùn)動(dòng)云圖

由圖2(a)可知，水流在入口處給定的速度2.6 m/s，水流經(jīng)過(guò)擴(kuò)散段發(fā)散撞入蜂窩網(wǎng)經(jīng)過(guò)水流方向和狀態(tài)的整合進(jìn)入穩(wěn)定段，水流變得有規(guī)律且穩(wěn)定流動(dòng)，但由于水流經(jīng)過(guò)蜂窩器的大面積摩擦，水流速度降低，在收縮段處水流速度急劇增大到4.64 m/s，經(jīng)過(guò)過(guò)渡段噴向開(kāi)口試驗(yàn)段，水流均勻穩(wěn)定流動(dòng)。觀察圖2(b)可以發(fā)現(xiàn)，試驗(yàn)中心速度相同，水流均勻流動(dòng)，但整體會(huì)有向周圍擴(kuò)散方向速度均勻減小的趨勢(shì)。

為研究水洞試驗(yàn)段是否符合試驗(yàn)要求的范圍，分別截取距離試驗(yàn)段入口0.7、0.9、1.1、1.3 m位置的中心速度云圖進(jìn)行分析。不同截面的速度云圖如圖3所示，選擇距離試驗(yàn)段入口0.7～1.3 m處作為試驗(yàn)范圍，滿足試驗(yàn)條件：(1)最大速度達(dá)到2 m/s；(2)滿足水流流速均勻性、穩(wěn)定性的條件；(3)滿足試驗(yàn)條件的試驗(yàn)區(qū)域大于試驗(yàn)區(qū)域最小限制范圍(300 mm×200 mm×68 mm)，滿足阻塞率不大于25%、長(zhǎng)度不小于3倍模型長(zhǎng)度的要求。

圖3 不同截面仿真速度云圖

3 水洞試驗(yàn)與標(biāo)定

3.1 動(dòng)力系統(tǒng)選型

流量Q用流體橫截面積A和流體速度v的乘積來(lái)表示

Q=Av

流體橫截面積A為過(guò)渡段圓形變方形的面積，圓的內(nèi)接正方形還可以增大水流速度，提高水流品質(zhì)[12]。圓直徑為80 mm，內(nèi)接正方形的邊長(zhǎng)取56 mm。經(jīng)模擬，此處速度取4.1 m/s，試驗(yàn)穩(wěn)流段方可達(dá)到試驗(yàn)速度要求，流量Q=46.3 m3/h。

根據(jù)水洞需要的流量和水頭損耗，選用U8000變頻水壓調(diào)節(jié)器及IRG80-125A型立式單級(jí)離心泵(性能參數(shù)見(jiàn)表1)構(gòu)成水洞動(dòng)力系統(tǒng)，如圖4所示，變頻水壓調(diào)節(jié)器壓力范圍為0.05～1.50 MPa，水泵出口處安裝不銹鋼止回閥，防止水倒流。

表1 立式單級(jí)離心泵性能

圖4 水洞動(dòng)力系統(tǒng)

3.2 水洞結(jié)構(gòu)及工作原理

水洞樣機(jī)如圖5所示，由水洞動(dòng)力系統(tǒng)1、控制面板2、水池3、閥4、管道5、過(guò)渡段二3.1、擴(kuò)散段3.2、穩(wěn)定段3.3、收縮段3.4、過(guò)渡段一3.5組成。水洞動(dòng)力系統(tǒng)1為水的流動(dòng)提供動(dòng)力，通過(guò)操作控制面板2控制動(dòng)力系統(tǒng)，通過(guò)改變壓力的大小來(lái)控制水流速度，水池3為開(kāi)口試驗(yàn)段，管道5連接水池和動(dòng)力系統(tǒng)，傳遞水流的作用；圖5(b)中各段為水流整流段，決定了水流流場(chǎng)質(zhì)量。

圖5 水洞樣機(jī)

3.3 水洞系統(tǒng)調(diào)試

通過(guò)改變水泵的水壓調(diào)節(jié)器的壓力控制電機(jī)轉(zhuǎn)速調(diào)節(jié)試驗(yàn)段水流流速。為保護(hù)水泵電機(jī)安全，減少葉片磨損，設(shè)置壓力0.5～1.0 MPa可調(diào)，測(cè)量位置為穩(wěn)流試驗(yàn)段中心，測(cè)量并記錄10個(gè)時(shí)均速度并取平均值。壓力與水洞流速對(duì)應(yīng)值見(jiàn)表2。

表2 壓力與水洞流速對(duì)應(yīng)關(guān)系

3.4 水洞流場(chǎng)校核

截面速度分布均勻度和流速穩(wěn)定性是表示流場(chǎng)質(zhì)量的兩個(gè)重要參數(shù)[13]。

(1)速度均勻性

速度均勻性表示某一時(shí)刻試驗(yàn)段斷面上速度分布的均勻性[14]。根據(jù)水洞流場(chǎng)模擬云圖，在試驗(yàn)段內(nèi)選取4個(gè)截面(距離入口0.7、0.9、1.1、1.3 m)，每個(gè)截面選取45個(gè)點(diǎn)，測(cè)試點(diǎn)位置如圖6所示，給定2 m/s的水速，將便攜式流速儀豎直置于測(cè)速點(diǎn)，單點(diǎn)測(cè)量各個(gè)截面點(diǎn)的時(shí)均速度，將數(shù)據(jù)繪制成云圖，如圖7所示。

圖6 內(nèi)測(cè)試點(diǎn)位置示意

圖7 不同截面實(shí)測(cè)速度云圖

分析試驗(yàn)得到的云圖可知：水流雖有從中心向四周擴(kuò)散的趨勢(shì)，但是中間速度均勻分布，可達(dá)到2 m/s的最大水流速度要求，距試驗(yàn)段入口0.7 m處截面內(nèi)260 mm×185 mm(x(-0.11, 0.15),y(-0.11, 0.075))矩形區(qū)域，1.1 m處截面內(nèi)200 mm×175 mm(x(-0.10, 0.10),y(-0.10, 0.075))矩形區(qū)域及1.3 m處截面內(nèi)200 mm×100 mm(x(-0.10, 0.10),y(-0.025, 0.075))矩形區(qū)域，滿足試驗(yàn)條件，所以滿足試驗(yàn)條件的試驗(yàn)區(qū)域大小為600 mm×200 mm×100 mm，位于距離試驗(yàn)段入口0.7～1.3 m處；或試驗(yàn)條件的試驗(yàn)區(qū)域大小為400 mm×200 mm×175 mm，位于距離試驗(yàn)段入口0.7～1.1 m處，全部大于最小限制范圍(300 mm×200 mm×68 mm)，滿足阻塞率不大于25%、長(zhǎng)度不小于3倍模型長(zhǎng)度的要求。計(jì)算得到此試驗(yàn)平臺(tái)的適用范圍為：模型長(zhǎng)不大于200 mm，模型來(lái)流方向的橫截面積不大于5 000 mm2。

(2)速度穩(wěn)定性

速度穩(wěn)定性是指在一時(shí)間段內(nèi)，每隔一定時(shí)間測(cè)量試驗(yàn)段內(nèi)相同位置的流速大致相同。在試驗(yàn)段中心，給定不同速度，每隔10 min測(cè)量一組試驗(yàn)數(shù)據(jù)，共測(cè)量6組。由表3可以看出，時(shí)均速度最大偏差在±0.2 m/s以內(nèi)，符合速度穩(wěn)定性要求。

表3 試驗(yàn)段速度穩(wěn)定性測(cè)量數(shù)據(jù) 單位：m·s-1

4 結(jié)論

基于相似原理設(shè)計(jì)一種針對(duì)低速水下機(jī)器人的水動(dòng)力試驗(yàn)系統(tǒng)。經(jīng)仿真模擬及試驗(yàn)分析，最終確定了該試驗(yàn)系統(tǒng)適用的模型尺度范圍和試驗(yàn)位置。此試驗(yàn)系統(tǒng)適用于速度不大于2 m/s，模型長(zhǎng)不大于200 mm，模型來(lái)流方向的橫截面積不大于5 000 mm2的水動(dòng)力試驗(yàn)，為低速水下機(jī)器人的水動(dòng)力試驗(yàn)提供了一種方案。