龐永杰，王慶云,2，李偉坡，王建

(1.哈爾濱工程大學(xué) 水下機(jī)器人技術(shù)重點實驗室,黑龍江 哈爾濱 150001；2.昆明船舶設(shè)備集團(tuán)有限公司,云南 昆明 650217)

螺旋槳及其運行對潛艇操縱性水動力影響的模型試驗研究

龐永杰1，王慶云1,2，李偉坡1，王建1

(1.哈爾濱工程大學(xué) 水下機(jī)器人技術(shù)重點實驗室,黑龍江 哈爾濱 150001；2.昆明船舶設(shè)備集團(tuán)有限公司,云南 昆明 650217)

潛艇模型；螺旋槳；水動力系數(shù)；操縱性試驗

在潛艇初步設(shè)計階段，其操縱性水動力系數(shù)的準(zhǔn)確獲取對操縱性的預(yù)報起關(guān)鍵性的作用。由于艇后螺旋槳的抽吸作用改變了艇體尾部流場分布和邊界層厚度，從而導(dǎo)致潛艇水動力性能產(chǎn)生變化。姚震球等[1]基于滑移網(wǎng)格對帶槳與不帶槳艇體阻力與艇尾流場進(jìn)行了數(shù)值模擬，PAN Yucun等[2]對無槳SUBOFF潛艇操縱性水動力系數(shù)進(jìn)行了數(shù)值模擬，Yingsan Wei等[3]對帶槳潛艇尾部流體噪聲進(jìn)行了模擬與分析，A.G.L. Holloway等[4]對不同外型潛艇主艇體完成了回轉(zhuǎn)試驗的數(shù)值模擬，ZHANG Nan等[5]對近水面潛艇直航槳艇干擾完成了數(shù)值模擬，A. Saeidinezhad等[6]對無槳潛艇完成了拘束模型攻角試驗，分析了周圍流場的變化機(jī)理，CAO Liushuai等[7]運用不同湍流模型對SUBOFF無槳潛艇回轉(zhuǎn)運動完成了數(shù)值模擬，S. Zaghi等[8]基于重疊網(wǎng)格對無槳潛艇操縱面完成了數(shù)值模擬，獲取了舵角水動力系數(shù)，WU Xiaocui等[9]對SUBOFF無槳潛艇直航、回轉(zhuǎn)等工況完成了數(shù)值模擬，結(jié)果與實驗值符合較好?？梢娭T多研究者均針對潛艇不帶槳情況獲取操縱性水動力系數(shù)，從而進(jìn)行潛艇操縱性能的預(yù)報，但是潛艇實際航行則是在帶槳情況下完成的，故無槳理想化的模型水動力系數(shù)的獲取將導(dǎo)致操縱性預(yù)報有所失真，從而影響潛艇的操縱性能。

故此，以某全附體潛艇縮比模型為研究對象，在循環(huán)水槽中利用 VPMM 完成帶槳和不帶槳潛艇模型斜航、純升沉、純俯仰、純橫蕩、純搖艏及舵角試驗，獲取該潛艇模型的操縱性線性水動力系數(shù)，并對比分析了變化的趨勢及原因。

1 試驗?zāi)Ｐ汀⒃O(shè)備及目的

1.1 模型及循環(huán)水槽

以某全附體潛艇為研究對象，該潛艇模型包括主艇體、圍殼、圍殼舵、十字型尾舵翼，主尺度及主要參數(shù)見表1。在舵翼前加裝輔助翼以保證水流光順的通過舵翼表面，有效提高舵效。圍殼舵、十字型尾舵翼斷面采用NACA0012翼型，模型外殼采用玻璃鋼材料制作；按照縮比模型試驗標(biāo)準(zhǔn)，完成槳模的繪制，并采用不銹鋼加工制作，基本參數(shù)見表2。

表1 潛艇模型基本參數(shù)

表2 螺旋槳基本參數(shù)

試驗在哈爾濱工程大學(xué)水下機(jī)器人技術(shù)重點實驗室循環(huán)水槽中完成，圖1為試驗中安裝完畢的帶槳全附體潛艇模型。循環(huán)水槽穩(wěn)流工作段尺寸：長7 m，寬1.7 m，水深1.5 m，即水槽橫截面積SA=1.7×1.5 m2；最大流速2 m/s，常用工作流速為0.3～1.8 m/s；流場不均勻度小于2%；流速不穩(wěn)定性小于2%；水平流向偏差小于±0.5°；湍流度0.7%～2.0%；自由液面駐波高度小于7 mm；從阻塞效應(yīng)來看，要求模型橫截面積與水槽工作段濕截面積SA之比小于3%，就本文而言，如果以主艇體橫截面積算，S/SA=2.86%，滿足試驗要求。

圖1 帶槳潛艇模型安裝示意圖Fig.1 Installation diagram of submarine model with propeller

1.2 模型配置與安裝及試驗流速選擇

模型加工以后，為完成操舵試驗，需在模型艏部和艉部加裝水密舵機(jī)及驅(qū)動螺旋槳旋轉(zhuǎn)的水密電機(jī)，然后進(jìn)行模型配置和浮態(tài)調(diào)整。按照《操縱性水動力系數(shù)試驗指導(dǎo)書》，為滿足模型在水中處于正浮狀態(tài)，同時減小非定常試驗過程中內(nèi)流對測力的影響，在模型空腔內(nèi)配備一定量與水密度相當(dāng)?shù)奈镔|(zhì)，并計算得到水動力模型的質(zhì)量、轉(zhuǎn)動慣量、重心位置，具體數(shù)值見表1；以浮心為坐標(biāo)原點，測力天平安裝要求測力中心點和浮心的縱向位置誤差應(yīng)該不大于1 mm，以減小試驗測力誤差。

流速選定：當(dāng)水動力模型Re超過臨界雷諾數(shù)Recr，模型的尺度效應(yīng)可以忽略。物體繞流運動中，邊界層內(nèi)流動狀態(tài)轉(zhuǎn)捩及流動分離不僅取決于來流Re數(shù)，而且依賴于來流湍流度、壁面粗糙度以及外部主流壓力梯度，因此試驗段來流湍流度ε對邊界層、阻力和升力都有很大影響，應(yīng)盡可能降低ε值，然而水槽的固有屬性使得水流湍流度ε一般為1.2%～2.0%[10]，其中湍流度ε對平板Recr影響的試驗結(jié)果表明：當(dāng)ε=1.2%～2.0%時，Recr≈(3.5～5.0)×105。根據(jù)《操縱性水動力系數(shù)試驗指導(dǎo)書》意見，試驗縮比模型長度為3.08 m,理論上在上述工作流速段試驗中Re數(shù)均超過Recr，但盡可能使Re較大，減小尺度效應(yīng)，同時避免流速過高，水流不穩(wěn)定，自由表面興波較大，水中氣泡含量過高等對試驗的影響，選取試驗流速為1.2 m/s。對應(yīng)Re約為3.6×106，滿足試驗要求。

1.3 試驗?zāi)康?/p>

為研究螺旋槳及其運行對潛艇操縱性水動力系數(shù)的影響，對無槳與帶槳全附體潛艇模型進(jìn)行拘束模型試驗，比較兩種模型下獲取的水動力系數(shù)的變化情況，為后期潛艇潛器操縱性預(yù)報仿真和CFD數(shù)值模擬奠定一定的基礎(chǔ)和提供參考依據(jù)。

2 基本理論與拘束模型試驗

2.1 潛艇操縱運動線性方程

潛艇在水下的運動可分為水平面運動和垂直面運動，在螺旋槳推動作用下，通過對潛艇的操縱面控制，以達(dá)到潛艇水下空間運動。采用根據(jù)國際拖曳水池會議(ITTC)和造船與輪機(jī)工程學(xué)會(SNAME)術(shù)語公報的體系推薦的潛艇坐標(biāo)系如圖2所示，慣性坐標(biāo)系為E-ξηζ，運動坐標(biāo)系為o-xyz。

不考慮橫滾和縱蕩，在線性假定下，可得潛艇在水平面和垂直面的四自由度操縱運動方程[11]：

水平面運動：

(1)

垂直面運動：

(2)

圖2 坐標(biāo)系示意圖Fig.2 The diagram of coordinate

由上述操縱運動分析，試驗內(nèi)容即可分為垂直面和水平面試驗，垂直面試驗包括沖角、純升沉、純俯仰、升降舵、圍殼舵試驗，為避免測力天平支桿的影響，模型以圍殼指向水槽槽底的方式安裝；水平面試驗包括漂角、純橫蕩、純搖艏、方向舵試驗[12]。

2.2 拘束模型定常試驗

2.2.1 斜航試驗

圖3 垂向力Z、俯仰力矩M隨垂向速度w的變化Fig.3 Vertical force Z,pitch moment M changes with vertical speed w

2.2.2 舵角試驗

表3 潛艇模型水動力系數(shù)

2.3 拘束模型簡諧運動試驗

簡諧運動試驗分為：純升沉、純俯仰、純橫蕩、純搖艏試驗。即驅(qū)動模型在來流中做純升沉或純俯仰運動，或?qū)⒛Ｐ蜐L轉(zhuǎn)90°安裝做純橫蕩或純搖艏運動。純俯仰和純搖艏運動時要求模型迎流方向與中心處軌跡曲線相切，來流速度U，中心處橫蕩速度v=Utanθ≈Uθ0sinωt，θ0為給定模型運動初始角位移，當(dāng)θ0較小時，一階近似為簡諧運動。此時后桿對前桿的滯后角ε滿足：

(3)

通過運動分解可把測得的拘束力(力矩)寫為方程(4)所示形式，用下標(biāo)in表示慣性力(同相分量)，用下標(biāo)out表示阻尼力(正交分量)，其余為常量：

(4)

圖4 垂向力Z、俯仰力矩M隨升降舵舵角δs變化Fig.4 Vertical force Z、pitching moment M changes with elevator δs

2.3.1 純升沉和純橫蕩試驗

垂直面純升沉運動線性方程：

(5)

模型純升沉運動參數(shù)：

(6)

測力天平與數(shù)據(jù)采集系統(tǒng)將拘束力測出并將其分解，從而求得純升沉水動力系數(shù)：

(7)

同理按上述方式試驗，可以得到純橫蕩水動力系數(shù)：

(8)

試驗得到對應(yīng)阻尼力、慣性力隨速度、加速度的變化趨勢，經(jīng)過線性擬合并結(jié)合表1中水動力模型參數(shù)代入式(7)、(8)可得到相應(yīng)水動力系數(shù)，無因次化處理得到的水動力系數(shù)匯總于表3中。以純升沉為例，測試結(jié)果如圖5、6所示。

圖5 純升沉慣性力Zin-aω2試驗結(jié)果Fig.5 Test results of pure heave Zin-aω2

圖6 純升沉慣性力矩Min-aω2試驗結(jié)果Fig.6 Test results of pure heave Min-aω2

2.3.2 純俯仰和純搖艏試驗

根據(jù)垂直面純俯仰運動線性方程為

(9)

由于俯仰角速度的平方為二階小量，忽略-mzGq2項。模型純俯仰運動參數(shù)為

(10)

測力天平與數(shù)據(jù)采集系統(tǒng)將拘束力測出并將其分解，從而求得純俯仰水動力系數(shù)：

(11)

同理按上述方式試驗，可以得到純搖艏水動力系數(shù)：

(12)

試驗得到對應(yīng)阻尼力、慣性力隨角速度、角加速度的變化趨勢，經(jīng)過線性擬合并結(jié)合表1中水動力模型參數(shù)代入式(11)、(12)可得到相應(yīng)水動力系數(shù)，無因次化處理得到的水動力系數(shù)匯總于表3中。以純俯仰為例，測試結(jié)果如圖7～10所示。

圖7 純俯仰阻尼力Zout-θ0ω試驗結(jié)果Fig.7 Test results of pure pitch Zout-θ0ω

圖8 純俯仰慣性阻尼力Zin-θω2試驗結(jié)果Fig.8 Test results of pure pitch Zin-θω2

圖9 純俯仰阻尼力矩Mout-θ0ω試驗結(jié)果Fig.9 Test results of pure pitch Mout-θ0ω

圖10 純俯仰慣性力矩Min-θω2試驗結(jié)果Fig.10 Test results of pure pitch Min-θω2

3 試驗結(jié)果與分析

按照上述試驗，所有結(jié)果匯總于表3，一共測試得到的22個線性類水動力系數(shù)中，對比無槳與帶槳值，不難發(fā)現(xiàn)有14個水動力系數(shù)呈減小趨勢，8個呈增加趨勢。

3.1 斜航

其次由于循環(huán)水槽本身尺寸的限制，導(dǎo)致斜航試驗時，螺旋槳尾流與循環(huán)水槽壁面和自由液面的干擾勢必給試驗測力帶來誤差，且沖角試驗的影響可能大于漂角試驗。因為本試驗是在主艇體截面積滿足S/SA=2.86%<3%條件下進(jìn)行的，而如果考慮圍殼等其他附體的影響，可有艇體最大橫截面積SM/SA=3.9%>3%，可見圍殼的存在產(chǎn)生了一定的阻塞效應(yīng)，且阻塞效應(yīng)主要發(fā)生在平行于圍殼的平面內(nèi)。在沖角試驗時，圍殼指向槽底，螺旋槳尾流場將與槽底或自由面耦合干擾加劇艇體周圍流場的變化，而漂角試驗時，模型繞軸x旋轉(zhuǎn)90°安裝，測力方向恰好正交于圍殼，故影響減小。

3.2 舵角試驗

3.3 簡諧運動試驗

4 結(jié)論

通過在循環(huán)水槽中，利用VPMM運動機(jī)構(gòu)對無槳與帶槳潛艇模型完成拘束模型試驗，獲取該潛艇模型的四自由度線性水動力系數(shù)，得到以下結(jié)論：

2)其余水動力系數(shù)變化基本保持在10%以內(nèi)；

3)由于水槽尺寸的限制，試驗過程中螺旋槳尾流場與水槽壁面和自由面的干擾，未免給水動力系數(shù)的測量帶來影響，有待進(jìn)一步改善；

4)此研究只是定性分析了水動力系數(shù)的變化原因，故也需CFD數(shù)值模擬相互對比論證并進(jìn)行流場變化物理機(jī)理的詳細(xì)分析與討論，進(jìn)一步研究尾部螺旋槳的運行對操縱性水動力系數(shù)的影響。

[1]姚震球, 高慧, 楊春蕾. 基于滑移網(wǎng)格的帶螺旋槳艇體尾流場數(shù)值分析方法[J]. 江蘇科技大學(xué)學(xué)報: 自然科學(xué)版, 2008, 22(2): 15-20. YAO Zhenqiu, GAO Hui, YANG Chunlei. Numerical simulation of interaction between submarine and propeller based on approach of sliding mesh[J]. Journal of Jiangsu university of science and technology: natural science edition, 2008, 22(2): 15-20.

[2]PAN Yucun, ZHANG Huaixin, ZHOU Qidou. Numerical Prediction of submarine hydrodynamic coefficients using CFD simulation[J]. Journal of hydrodynamics, series B, 2012, 24(6): 840-847.

[3]WEI Yingsan, WANG Yongsheng. Unsteady hydrodynamics of blade forces and acoustic responses of a model scaled submarine excited by propeller’s thrust and side-forces[J]. Journal of sound and vibration, 2013, 332(8): 2038-2056.

[4]ZHANG J T, MAXWELL J A, GERBER A G, et al. Simulation of the flow over axisymmetric submarine hulls in steady turning[J]. Ocean engineering, 2013, 57: 180-196.

[5]ZHANG Nan, ZHANG Shengli. Numerical simulation of hull/propeller interaction of submarine in submergence and near surface conditions[J]. Journal of hydrodynamics, ser. b, 2014, 26(1): 50-56.

[6]SAEIDINEZHAD A, DEHGHAN A A, MANSHADI M D. Experimental investigation of hydrodynamic characteristics of a submersible vehicle model with a non-axisymmetric nose in pitch maneuver[J]. Ocean engineering, 2015, 100: 26-34.

[7]CAO Liushuai, ZHU Jun, WAN Wenbin. Numerical Investigation of submarine hydrodynamics and flow field in steady turn[J]. China ocean engineering, 2016, 30(1): 57-68.

[8]ZAGHI S, DI MASCIO A, BROGLIA R, et al. Application of dynamic overlapping grids to the simulation of the flow around a fully-appended submarine[J]. Mathematics and computers in simulation, 2015, 116: 75-88.

[9]WU Xiaocui, WANG Yiwei, HUANG Chenguang, et al. An effective CFD approach for marine-vehicle maneuvering simulation based on the hybrid reference frames method[J]. Ocean engineering, 2015, 109: 83-92.

[10]高霓. 微小型水下潛器近自由液面操縱性預(yù)報[D]. 哈爾濱: 哈爾濱工程大學(xué), 2013: 68-78. GAO Ni. Maneuverability prediction of mini underwater vehicle near free surface[D]. Harbin: Harbin Engineering University, 2013: 68-78.

[11]丁正良, 胡瑞芝, 孫百超, 等. 水平型循環(huán)水槽的水動力性能[Z]. 哈爾濱: 哈爾濱船舶工程學(xué)院, 1985: 5-10. DING Zhengliang, HU Ruizhi, SUN Baichao, et al. Hydrodynamic performance of horizontal circulating water channel[Z]. Harbin: Harbin Shipbuilding Engineering University, 1985: 5-10.

[12]王慶云, 龐永杰, 李偉坡, 等. 系列舵翼潛艇水動力系數(shù)數(shù)值計算及試驗研究[J]. 艦船科學(xué)技術(shù), 2015, 37(11): 21-26. WANG Qingyun, PANG Yongjie, LI Weipo, et al. Numerical calculation and experimental study of hydrodynamic coefficients of submarine of a series of rudder and wing[J]. Ship science and technology, 2015, 37(11): 21-26.

Model test study of influence of propeller and its rotation on hydrodynamics of submarine maneuverability

PANG Yongjie1, WANG Qingyun1,2, LI Weipo1,WANG Jian1

(1.Science and Technology on Underwater Vehicle Laboratory, Harbin Engineering University, Harbin 150001, China；2.Kunming Shipbuilding Equipment Co., Ltd., Kunming 650217,China)

submarine model; propeller; hydrodynamic coefficient; maneuverability tests

2015-12-16.

時間：2016-12-21.

國家自然科學(xué)基金項目(2012T50331).

龐永杰(1955-), 男, 教授,博士生導(dǎo)師； 王慶云(1989-), 男, 碩士研究生.

王慶云, E-mail：wangqingyun@hrbeu.edu.cn.

10.11990/jheu.201512054

U661

A

1006-7043(2017)01-0109-06

龐永杰，王慶云，李偉坡，等. 螺旋槳及其運行對潛艇操縱性水動力影響的模型試驗研究[J]. 哈爾濱工程大學(xué)學(xué)報， 2017, 38(1): 110-115. PANG Yongjie, WANG Qingyun, LI Weipo，et al. Model test study of influence of propeller and its rotation on hydrodynamics of submarine maneuverability[J]. Journal of Harbin Engineering University, 2017, 38(1): 110-115.

網(wǎng)絡(luò)出版地址：http://www.cnki.net/kcms/detail/23.1390.u.20161221.1524.006.html