亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

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

        2017-04-21 03:33:24馮冠華李智剛馮迎賓梁洪光程陽銳吳東華
        海洋學(xué)研究 2017年1期
        關(guān)鍵詞:模型

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

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

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

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

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

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

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

        0 引言

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

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

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

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

        多金屬結(jié)核概念車浮游體(以下簡稱浮游體)主要組成部分包括:水下電機、液壓泵、儲油艙、補償器、推進器、電子艙、框架、攝像機及其他傳感器等。先根據(jù)集礦車對液壓系統(tǒng)的具體需求,確定液壓站尺寸、重量等參數(shù)。再根據(jù)對浮游體的作業(yè)深度、巡航速度、承載能力、總重量和最大尺寸范圍等參數(shù)的具體要求,設(shè)計浮游體的各部分外形并整體布局,確保滿足其所需要的運動能力和穩(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ù)值仿真計算效率和準確度較高,可以滿足設(shè)計要求。計算流體動力學(xué)仿真軟件ANSYS CFX是在流動基本方程控制下對流動的數(shù)值模擬,其采用有限體積法對控制方程和計算區(qū)域進行離散,添加邊界條件后求解。本文計算數(shù)學(xué)模型采用工程湍流數(shù)值計算中最為廣泛的雷諾時均算法(RANS算法)和兩方程的k-ε模型,其中k是湍流動能,ε是湍流耗散率[5]。

        雷諾時均RANS方程組為:

        (1)

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

        (2)

        上式中湍流動能生成項為:

        (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進一步分析、簡化已建立的浮游體幾何模型,建立浮游體流體動力學(xué)計算模型。然后,將浮游體流體動力學(xué)計算模型導(dǎo)入到網(wǎng)格劃分軟件ANSYS ICEM CFD中生成網(wǎng)格,利用計算流體動力學(xué)軟件ANSYS CFX分析得到浮游體流體動力學(xué)模型。最后,進一步計算得到浮游體阻力等數(shù)據(jù),并采用Matlab曲線擬合功能,得到浮游體相應(yīng)的阻力系數(shù)Cd曲線,并研究其阻力特性。

        3.1 網(wǎng)格生成

        根據(jù)浮游體整體外形尺寸,取流域的長、寬、高為浮游體本體長、寬、高的5倍并取整,最終確定計算流域為長12 000 mm、寬7 000 mm、高7 000 mm的長方體域。利用SolidWorks對浮游體整體外形進行適當簡化,以浮游體幾何模型的重心為原點建立坐標系,建立流體動力學(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],即按照攻角及漂角進行取值:

        (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 阻力特性分析

        在研究浮游體阻力時,分為主體部分和附屬體部分。浮游體主體部分阻力包括摩擦阻力、形狀阻力和興波阻力。通常把附屬體部分阻力的摩擦阻力計入到主體的摩擦阻力中,剩余的阻力作為一項,為附體阻力[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隨攻角的增加而減少。而在負攻角工況下,阻力系數(shù)Cd隨攻角絕對值的增加而增加。盡管阻力系數(shù)Cd隨攻角的變化而增加或減少,但阻力系數(shù)變化幅度都很小,這一點從圖8和圖9的兩種攻角工況下浮游體表面壓力云圖得到了驗證。浮游體表面壓力分布對其攻角的改變不敏感,說明本文設(shè)計的浮游體整體外形滿足設(shè)計需求,驗證了外形設(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è)計成獨立的浮游體,浮游體與集礦車協(xié)同工作,使深海集礦車獲得了更好的穩(wěn)定性及操縱性,提高采集率。運用SolidWorks完成了多金屬結(jié)核概念車浮游體外形設(shè)計,建立了其流體動力學(xué)計算模型。借助ANSYS CFX軟件,對浮游體進行了流體動力學(xué)仿真分析,模擬了其在攻角為0°、±2°、±4°、±6°、±8°、±10°和速度為1~6 kn工況下的水動力性能,并進一步研究了浮游體的阻力特性,得知本文設(shè)計的浮游體表面壓力分布對其攻角的改變不敏感。驗證了外形設(shè)計方案的可行性,為浮游體下一步控制策略研究提供了重要依據(jù),并為未來深海集礦工程實際應(yīng)用打下了基礎(chǔ)。

        致謝 感謝中國科學(xué)院沈陽自動化研究所及長沙礦治研究院責任公司給予項目經(jīng)費資助,并在整個研究過程中提供實驗室及實驗設(shè)備等;感謝陳質(zhì)二博士及王亞興博士在整個仿真實驗過程中給予指導(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.

        戴瑜,劉少軍.深海采礦機器人研究:現(xiàn)狀與發(fā)展[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].機械工程學(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].機械工程學(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.

        孟令帥,林揚,鄭榮,等.模塊化自主水下機器人的機械設(shè)計與實現(xiàn)[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.

        胡志強,衣瑞文,林揚,等.基于隨體坐標系的水下機器人水動力數(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].機械設(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”計劃)項目資助(2012AA091201)

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

        *通訊作者:李智剛(1970-),男,研究員,主要從事水下機器人控制系統(tǒng)研究、動力傳輸及保護技術(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.

        猜你喜歡
        模型
        一半模型
        一種去中心化的域名服務(wù)本地化模型
        適用于BDS-3 PPP的隨機模型
        提煉模型 突破難點
        函數(shù)模型及應(yīng)用
        p150Glued在帕金森病模型中的表達及分布
        函數(shù)模型及應(yīng)用
        重要模型『一線三等角』
        重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
        3D打印中的模型分割與打包
        少妇性荡欲视频| 日本刺激视频一区二区| 国产情侣一区二区三区| 人人妻人人澡人人爽人人精品97| 在线视频制服丝袜中文字幕| 中文字幕色视频在线播放| 青青草视频在线观看入口| 亚洲av色香蕉一区二区三区老师| 91丝袜美腿亚洲一区二区| 久久视频在线| 国产小毛片| 手机在线观看亚洲av| 人妻av有码中文字幕| 小宝极品内射国产在线| 精品无码久久久九九九AV| 亚洲成av人片在久久性色av| 日本va中文字幕亚洲久伊人| 国产av国片精品jk制服| 亚洲永久精品ww47| 无码一区二区三区人| 国产三级精品三级在线专区2| 美女脱了内裤张开腿让男人桶网站| 五十路熟女一区二区三区| av永久天堂一区二区三区蜜桃 | 国产一区二区波多野结衣| 麻豆AV免费网站| 国产一区二区三区精品成人爱| 免费观看全黄做爰大片| 亚洲av无码国产剧情| 国产精品久久久久亚洲| 国产不卡在线观看视频| 国产成人无码精品久久二区三区| 成人在线激情网| 国产欧美激情一区二区三区| 日本午夜艺术一区二区| 女人被狂躁c到高潮视频| 在线观看欧美精品| 亚洲大片一区二区三区四区| 午夜dv内射一区二区| 亚洲熟妇无码av另类vr影视| 亚洲一级无码AV毛片久久|