多金屬結(jié)核概念車浮游體的外形設(shè)計及阻力特性分析

馮冠華，李智剛，馮迎賓，梁洪光，程陽銳，吳東華

(1. 沈陽理工大學(xué) 機(jī)械工程學(xué)院,遼寧 沈陽 110159；2. 中國科學(xué)院 沈陽自動化研究所 機(jī)器人學(xué)國家重點(diǎn)實(shí)驗(yàn)室，遼寧 沈陽 110016；3. 長沙礦治研究院有限責(zé)任公司 深海礦產(chǎn)資源開發(fā)利用技術(shù)國家重點(diǎn)實(shí)驗(yàn)室，湖南 長沙 410012)

多金屬結(jié)核概念車浮游體的外形設(shè)計及阻力特性分析

馮冠華1，李智剛*2，馮迎賓2，梁洪光2，程陽銳3，吳東華3

(1. 沈陽理工大學(xué) 機(jī)械工程學(xué)院,遼寧 沈陽 110159；2. 中國科學(xué)院 沈陽自動化研究所 機(jī)器人學(xué)國家重點(diǎn)實(shí)驗(yàn)室，遼寧 沈陽 110016；3. 長沙礦治研究院有限責(zé)任公司 深海礦產(chǎn)資源開發(fā)利用技術(shù)國家重點(diǎn)實(shí)驗(yàn)室，湖南 長沙 410012)

針對深海采礦實(shí)際需求，提出將集礦車上的液壓站與集礦車分離，設(shè)計成相對獨(dú)立、懸浮于集礦車前部上方的浮游體。設(shè)計了多金屬結(jié)核概念車浮游體(以下簡稱浮游體)外形，并利用三維設(shè)計軟件SolidWorks建立了其三維幾何模型。在此基礎(chǔ)上，對幾何模型進(jìn)行簡化處理，建立了浮游體流體動力學(xué)計算模型。利用流體動力學(xué)仿真分析軟件ANSYS CFX，采用雷諾時均算法(RANS算法)和兩方程的k-ε湍流模型，得到了不同工況的阻力數(shù)據(jù)。采用Matlab曲線擬合功能，研究了浮游體阻力特性，驗(yàn)證了外形設(shè)計方案的可行性。

深海采礦；浮游體外形；流體動力學(xué)仿真分析；阻力特性

0 引言

隨著陸地礦產(chǎn)資源的逐漸枯竭和人類對金屬資源需求的日益增加，海底礦產(chǎn)必將成為21世紀(jì)的重要資源[1]。但是這些海底礦產(chǎn)大都存在于深海，必須依靠深海采礦裝備——集礦車對其進(jìn)行開發(fā)，集礦車設(shè)計的好壞將直接影響深海采礦的成敗。在深海采礦過程中，集礦車要具有適應(yīng)深海惡劣環(huán)境和海底地貌變化、連續(xù)無故障作業(yè)、越障越溝、保證采集率和不給深海生態(tài)環(huán)境帶來過度破壞等能力[1-3]。

液壓站是集礦車的重要組成部分，其在整個車體中的具體位置如何布置，并不影響其功能的實(shí)現(xiàn)。但液壓站的體積和重量較大，會影響集礦車整體布局，增加整車重量，使集礦車結(jié)構(gòu)復(fù)雜化，拆卸維修困難，且增加整車能耗。另外，液壓站工作時會產(chǎn)生強(qiáng)烈震動，嚴(yán)重影響集礦車導(dǎo)航控制系統(tǒng)的正常工作。

本文正是綜合考慮以上因素，提出將液壓站從集礦車中分離出來，將液壓站設(shè)計成一個相對獨(dú)立于集礦車、懸浮于集礦車前部上方的浮游體。采用懸浮液壓站技術(shù)，一方面能縮小集礦車體積、減輕整車重量，便于集礦車整體優(yōu)化布局，提高其在海底稀軟底質(zhì)上的行走能力；另一方面浮游體作業(yè)時懸浮于集礦車的前部上方，能有效擴(kuò)展監(jiān)控范圍，保證集礦車作業(yè)過程的安全。本文借助計算流體動力學(xué)仿真分析軟件ANSYS CFX對浮游體的阻力特性進(jìn)行了數(shù)值仿真研究。

1 浮游體外形設(shè)計

多金屬結(jié)核概念車浮游體(以下簡稱浮游體)主要組成部分包括：水下電機(jī)、液壓泵、儲油艙、補(bǔ)償器、推進(jìn)器、電子艙、框架、攝像機(jī)及其他傳感器等。先根據(jù)集礦車對液壓系統(tǒng)的具體需求，確定液壓站尺寸、重量等參數(shù)。再根據(jù)對浮游體的作業(yè)深度、巡航速度、承載能力、總重量和最大尺寸范圍等參數(shù)的具體要求，設(shè)計浮游體的各部分外形并整體布局，確保滿足其所需要的運(yùn)動能力和穩(wěn)定性[4]。設(shè)計的多金屬結(jié)核概念車浮游體參數(shù)為：長度2 419 mm、寬度1 340 mm、高度1 487 mm、重量(空氣中)1 500 kg、排水量1 328 L和穩(wěn)心高11 mm。建立的多金屬結(jié)核概念車浮游體幾何模型如圖1所示。

圖1 多金屬結(jié)核概念車浮游體幾何模型Fig.1 Geometric model of floating body for polymetallic nodule concept vehicle

2 流體動力學(xué)理論分析

浮游體整體外形結(jié)構(gòu)，決定其工作時表面的流場。采用CFD(計算流體動力學(xué)仿真分析)軟件ANSYS CFX來研究浮游體阻力特性，其數(shù)值仿真計算效率和準(zhǔn)確度較高，可以滿足設(shè)計要求。計算流體動力學(xué)仿真軟件ANSYS CFX是在流動基本方程控制下對流動的數(shù)值模擬，其采用有限體積法對控制方程和計算區(qū)域進(jìn)行離散，添加邊界條件后求解。本文計算數(shù)學(xué)模型采用工程湍流數(shù)值計算中最為廣泛的雷諾時均算法(RANS算法)和兩方程的k-ε模型，其中k是湍流動能，ε是湍流耗散率[5]。

雷諾時均RANS方程組為：

(1)

兩方程的k-ε湍流模型為：

(2)

上式中湍流動能生成項(xiàng)為：

(3)

(4)

在兩方程的k-ε湍流模型中，湍流黏性系數(shù)定義為：

(5)

兩方程k-ε湍流模型系數(shù)為：

Cμ=0.09，Cε1=1.44，Cε2=1.92，σk=1.0，σε=1.3。

3 流體動力學(xué)計算分析

首先，在完成多金屬結(jié)核概念車浮游體整體結(jié)構(gòu)設(shè)計的基礎(chǔ)上，確定其流域范圍，同時利用三維設(shè)計軟件SolidWorks進(jìn)一步分析、簡化已建立的浮游體幾何模型，建立浮游體流體動力學(xué)計算模型。然后，將浮游體流體動力學(xué)計算模型導(dǎo)入到網(wǎng)格劃分軟件ANSYS ICEM CFD中生成網(wǎng)格，利用計算流體動力學(xué)軟件ANSYS CFX分析得到浮游體流體動力學(xué)模型。最后，進(jìn)一步計算得到浮游體阻力等數(shù)據(jù)，并采用Matlab曲線擬合功能，得到浮游體相應(yīng)的阻力系數(shù)Cd曲線，并研究其阻力特性。

3.1 網(wǎng)格生成

根據(jù)浮游體整體外形尺寸，取流域的長、寬、高為浮游體本體長、寬、高的5倍并取整，最終確定計算流域?yàn)殚L12 000 mm、寬7 000 mm、高7 000 mm的長方體域。利用SolidWorks對浮游體整體外形進(jìn)行適當(dāng)簡化，以浮游體幾何模型的重心為原點(diǎn)建立坐標(biāo)系，建立流體動力學(xué)計算模型，如圖2所示。根據(jù)浮游體已有的各個局部尺寸，同時考慮后續(xù)網(wǎng)格劃分計算時間和精度等，將流體動力學(xué)模型分為IN、OUT、WALL、BIGBODY、MIDBODY、SMALLBODY和RING七個Part，各個Part的意義如表1所示。

圖2 流體動力學(xué)計算模型Fig.2 Computational model of fluid dynamics

表1 各個Part的意義Tab.1 Meaning of parts

采用ANSYS ICEM CFD軟件的非結(jié)構(gòu)自動體網(wǎng)格生成技術(shù)，定義全局網(wǎng)格參數(shù)，即設(shè)置全局網(wǎng)格比例尺(Scale factor)為1.0，最大網(wǎng)格尺寸(Max element)為400.0。定義全局殼網(wǎng)格參數(shù)，即網(wǎng)格類型和網(wǎng)格生成方法分別為Quad Dominant和Patch Dependent。定義全局體網(wǎng)格參數(shù)，即網(wǎng)格類型和網(wǎng)格生成方法分別為Tetra/Mixed和Robust(Octree)。定義全局棱柱網(wǎng)格參數(shù)，即設(shè)置生長方式(Growth law )為exponential，初始高度(Initial height)為0.1，高度比(Height ratio)為1.2，邊界層數(shù)(Number of layers)為3層。為了保證計算精度，同時又減少網(wǎng)格規(guī)模，提高數(shù)值計算效率，不同Part上定義了不同的最大網(wǎng)格尺寸，如表2所示。

表2 各個Part的網(wǎng)格尺寸Tab.2 Mesh size for parts

在定義網(wǎng)格參數(shù)基礎(chǔ)上，生成網(wǎng)格，如圖3所示。網(wǎng)格質(zhì)量如圖4所示，網(wǎng)格質(zhì)量都大于0.2，滿足后續(xù)計算要求。

圖3 生成網(wǎng)格Fig.3 Generated mesh

圖4 網(wǎng)格質(zhì)量Fig.4 Mesh quality

3.2 流體動力學(xué)計算分析

流體動力學(xué)數(shù)值計算及后處理采用仿真軟件ANSYS CFX，在前處理中定義邊界條件時，保持浮游體本體在流域中的位置不變，而改變流域的朝向[6]，如圖5所示。

圖5 計算流域示意圖Fig.5 Sketch map of computational basin

流域入口按照速度邊界條件[7-8]，即按照攻角及漂角進(jìn)行取值：

(6)

其中，Vin是入口處來流速度，α是攻角，β是漂角。

流域出口按照壓力邊界條件，即取相對靜壓值為零。流域側(cè)壁按照自由滑移壁面邊界條件，而浮游體本體表面按照無滑移壁面邊界條件。

設(shè)定各個參數(shù)的收斂殘差值Residual Target為1.0E-4。設(shè)置最大迭代次數(shù)Max. Iterations為80，迭代將近60步后計算結(jié)果開始收斂，圖6所示為浮游體流體動力學(xué)計算殘差變化情況，可以判斷計算結(jié)果已經(jīng)收斂。最終計算得到了浮游體在每個攻角(分別為0°、±2°、±4°、±6°、±8°和±10°)下不同速度(1～6 kn)時的阻力及阻力矩數(shù)據(jù)，共66組數(shù)據(jù)。

圖6 計算殘差變化情況Fig.6 Variation of residual

3.3 阻力特性分析

在研究浮游體阻力時，分為主體部分和附屬體部分。浮游體主體部分阻力包括摩擦阻力、形狀阻力和興波阻力。通常把附屬體部分阻力的摩擦阻力計入到主體的摩擦阻力中，剩余的阻力作為一項(xiàng)，為附體阻力[9-10]。浮游體水下航行的總阻力為：

R=RF+RX+RW+RAP

(7)

其中：RF是摩擦阻力，RX是形狀阻力，RW是興波阻力，RAP是附體阻力。

本文所設(shè)計的浮游體在水下航行時，其工作深度超過三分之一總長，浮游體的興波阻力幾乎為零，故不考慮其興波阻力，則浮游體水下航行時的阻力主要由摩擦阻力、形狀阻力和附體阻力三部分組成。即：

R=RF+RX+RAP

(8)

上式中摩擦阻力和形狀阻力分別為：

(9)

(10)

其中：ρ是海水密度；V是浮游體速度；S是浮游體浸濕表面積；(CF+ΔCF)是摩擦阻力系數(shù)，CF是光滑平板的摩擦阻力系數(shù)，ΔCF是粗糙度附加增量；CE是形狀阻力系數(shù)。

將附屬體部分的阻力并入到形狀阻力中去，主要從摩擦阻力和形狀阻力計算分析浮游體阻力特性。此時浮游體水下航行的總阻力R為：

(11)

其中：Cd為阻力系數(shù)。

可得到阻力系數(shù)Cd定義如下：

(12)

利用通過仿真分析軟件ANSYS CFX得到的阻力數(shù)據(jù)，采用Matlab最小二乘法擬合，并結(jié)合上述阻力系數(shù)Cd的定義，得到了不同攻角工況下阻力系數(shù)Cd計算結(jié)果，并擬合出阻力系數(shù)曲線，如圖7所示。

圖7 阻力系數(shù)曲線Fig.7 Drag coefficient curve

通過不同攻角工況下阻力系數(shù)曲線圖可以看出，在正攻角工況下，阻力系數(shù)Cd隨攻角的增加而減少。而在負(fù)攻角工況下，阻力系數(shù)Cd隨攻角絕對值的增加而增加。盡管阻力系數(shù)Cd隨攻角的變化而增加或減少，但阻力系數(shù)變化幅度都很小，這一點(diǎn)從圖8和圖9的兩種攻角工況下浮游體表面壓力云圖得到了驗(yàn)證。浮游體表面壓力分布對其攻角的改變不敏感，說明本文設(shè)計的浮游體整體外形滿足設(shè)計需求，驗(yàn)證了外形設(shè)計方案的可行性。

圖8 0°攻角工況下浮游體表面壓力云圖Fig.8 Pressure contour of the surface of floating body at the angle of attack 0°

圖9 -6°攻角工況下浮游體表面壓力云圖Fig.9 Pressure contour of the surface of floating body at the angle of attack -6°

4 小結(jié)

本文針對深海采礦實(shí)際需求，提出了一種新的設(shè)計思路，將液壓站分離出來設(shè)計成獨(dú)立的浮游體，浮游體與集礦車協(xié)同工作，使深海集礦車獲得了更好的穩(wěn)定性及操縱性，提高采集率。運(yùn)用SolidWorks完成了多金屬結(jié)核概念車浮游體外形設(shè)計，建立了其流體動力學(xué)計算模型。借助ANSYS CFX軟件，對浮游體進(jìn)行了流體動力學(xué)仿真分析，模擬了其在攻角為0°、±2°、±4°、±6°、±8°、±10°和速度為1～6 kn工況下的水動力性能，并進(jìn)一步研究了浮游體的阻力特性，得知本文設(shè)計的浮游體表面壓力分布對其攻角的改變不敏感。驗(yàn)證了外形設(shè)計方案的可行性，為浮游體下一步控制策略研究提供了重要依據(jù)，并為未來深海集礦工程實(shí)際應(yīng)用打下了基礎(chǔ)。

致謝 感謝中國科學(xué)院沈陽自動化研究所及長沙礦治研究院責(zé)任公司給予項(xiàng)目經(jīng)費(fèi)資助，并在整個研究過程中提供實(shí)驗(yàn)室及實(shí)驗(yàn)設(shè)備等；感謝陳質(zhì)二博士及王亞興博士在整個仿真實(shí)驗(yàn)過程中給予指導(dǎo)及建議；感謝張洪彬碩士在論文選題、構(gòu)思及撰寫過程中給予指導(dǎo)及建議。

[1] DAI Yu, LIU Shao-jun. Researches on deep ocean mining robots: status and development[J]. Robot,2013,35(3):363-375.

戴瑜,劉少軍.深海采礦機(jī)器人研究：現(xiàn)狀與發(fā)展[J].機(jī)器人,2013,35(3):363-375.

[2] LIU Shao-jun, LIU Chang, DAI Yu. Status and progress on researches and developments of deep ocean mining equipments[J]. Journal of Mechanical Engineering,2014,50(2):8-18.

劉少軍,劉暢,戴瑜.深海采礦裝備研發(fā)的現(xiàn)狀與進(jìn)展[J].機(jī)械工程學(xué)報,2014,50(2):8-18.

[3] DAI Yu, LIU Shao-jun. Establishment of the dynamic model of the total deep ocean mining system and fast simulation of its integrated operation process[J]. Journal of Mechanical Engineering,2012,48(9):79-88.

戴瑜,劉少軍.深海采礦整體系統(tǒng)動力學(xué)建模及聯(lián)動開采作業(yè)過程仿真分析[J].機(jī)械工程學(xué)報,2012,48(9):79-88.

[4] MENG Ling-shuai, LIN Yang, ZHENG Rong, et al. Mechanical design and implementation of a modular autonomous underwater vehicle[J]. Robot,2016,38(4):395-401.

孟令帥,林揚(yáng),鄭榮,等.模塊化自主水下機(jī)器人的機(jī)械設(shè)計與實(shí)現(xiàn)[J].機(jī)器人,2016,38(4):395-401.

[5] ZHANG De-liang. A course in computational fluid dynamics[M]. Beijing: Higher Education Press,2010:19-36.

張德良.計算流體力學(xué)教程[M].北京：高等教育出版社,2010：19-36.

[6] HU Zhi-qiang, YI Rui-wen, LIN Yang, et al. Numerical calculation methods for hydrodynamics of unmanned underwater vehicles based on body-fixed coordinate frames[J]. Chinese Science Bulletin,2013,58(Suppl.Ⅱ):55-66.

胡志強(qiáng),衣瑞文,林揚(yáng),等.基于隨體坐標(biāo)系的水下機(jī)器人水動力數(shù)值計算方法[J].科學(xué)通報,2013,58(增刊Ⅱ):55-66.

[7] SOHANKAR A, NORBERG C, DAVIDSON L. Low-reynolds-number flow around a square cylinder at incidence: study of blockage, onset of vortex shedding and outlet boundary condition[J]. International Journal for Number Methods in Fluids,1998,26(1):39-56.

[8] WANG Mei-ting, QI Yong-feng, DAI Zhi-guang, et al. Research on shape and linear sailing resistance character of mini underwater robot[J]. Machinery Design & Manufacture,2013,12(1):35-38.

王妹婷,齊永鋒,戴志光,等.小型水下機(jī)器人外形及其直航阻力特性研究[J].機(jī)械設(shè)計與制造,2013,12(1):35-38.

[9] HAO Ying-ze，LIN Fan-cai，LIU Bai-shun. Sailed resistance analysis and drag reduction measures of submarine[C]//Proceedings of oversize ship maneuvering and ship security and management，2003:141-144.

郝英澤，林凡彩，劉百順.淺析潛艇水下航行阻力及減阻措施[C]//特大型船舶操縱和船舶安全與管理論文集，2003:141-144.

[10] NGUYEN B，HOPKIN D. Modeling autonomous underwater vehicle(AUV)operations in mine hunting[C]//Oceans Europe #Brest France，2005:533-538.

The contour design and analysis for drag characteristic of floating body for polymetallic nodule concept vehicle

FENG Guan-hua1, LI Zhi-gang*2, FENG Ying-bin2, LIANG Hong-guang2, CHENG Yang-rui3, WU Dong-hua3

(1.SchoolofMechanicalEngineering,ShenyangLigongUniversity,Shenyang110159,China; 2.StateKeyLaboratoryofRobotics,ShenyangInstituteofAutomation,ChineseAcademyofSciences,Shenyang110016,China; 3.StateKeyLaboratoryofExploitationandUtilizationofDeep-SeaMineralResource,ChangshaResearchInstituteofMiningandMetallurgyCo.Ltd.,Changsha410012,China)

Aiming at the practical demand of deep ocean mining, the hydraulic station on the nodule collector was proposed to separate from the nodule collector and hydraulic station was designed as a floating body, relatively independent and suspended above the front of the nodule collector. Contour of floating body for polymetallic nodule concept vehicle (referred to as floating body) was designed, and its 3D model was established adopting SolidWorks. The hydrodynamic and computational model of the floating body was established, after simplifying the 3D model. Reynolds average navier-stokes algorithm (RANS algorithm) and turbulence model were applied, and drag data under the condition of different working conditions was obtained adopting fluid dynamics simulation and analysis software ANSYS CFX. Drag characteristic of floating body was studied adopting curve fitting function of Matlab, and the feasibility of contour design was verified.

deep ocean mining; contour of floating body; fluid dynamics simulation and analysis; drag characteristic

10.3969/j.issn.1001-909X.2017.01.010.

2015-12-03

2017-01-09

國家高技術(shù)研究發(fā)展計劃(“863”計劃)項(xiàng)目資助(2012AA091201)

馮冠華(1991-)，男，河南淇縣人，主要從事水下機(jī)器人水動力分析以及載體技術(shù)研究。E-mail:fengguahua@sia.cn

*通訊作者：李智剛(1970-)，男，研究員，主要從事水下機(jī)器人控制系統(tǒng)研究、動力傳輸及保護(hù)技術(shù)研究、產(chǎn)品開發(fā)及工程化。E-mail:lzg@sia.cn

TP242.3

A

1001-909X(2017)01-0080-06

10.3969/j.issn.1001-909X.2017.01.010

馮冠華，李智剛，馮迎賓，等.多金屬結(jié)核概念車浮游體的外形設(shè)計及阻力特性分析[J].海洋學(xué)研究,2017,35(1):80-85,

FENG Guan-hua, LI Zhi-gang, FENG Ying-bin, et al. The contour design and analysis for drag characteristic of floating body for polymetallic nodule concept vehicle[J].Journal of Marine Sciences,2017,35(1):80-85, doi:10.3969/j.issn.1001-909X.2017.01.010.