CN111079286B - 一种层冰工况船舶推力减额分数数值估计方法 - Google Patents
一种层冰工况船舶推力减额分数数值估计方法 Download PDFInfo
- Publication number
- CN111079286B CN111079286B CN201911294887.9A CN201911294887A CN111079286B CN 111079286 B CN111079286 B CN 111079286B CN 201911294887 A CN201911294887 A CN 201911294887A CN 111079286 B CN111079286 B CN 111079286B
- Authority
- CN
- China
- Prior art keywords
- ice
- thrust
- self
- ship
- model
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供的是一种层冰工况船舶推力减额分数数值估计方法。一、基于目标船舶开展裸船体的水阻力数值模拟,采集不同航速下的裸船体的水阻力。二、基于目标船舶开展船后自航数值模拟,采集不同航速、转速下的水阻力、螺旋桨有效推力以及扭矩。通过推力减额定义求解出船舶水中的推力减额。三、通过带附着冰的目标船自航模型开展船后自航模拟,采集不同航速、转速下的水阻力、螺旋桨有效推力以及扭矩。基于有效推力的差额来完成冰中的推力减额,求得冰水中的推力减额分数。该方法考虑了冰对伴流场的影响以及冰对螺旋桨的阻塞作用,忽视了碎冰与桨直接接触的影响,可以有效解决目前采用常规船舶计算推力减额方法误差较大的问题。
Description
技术领域
本发明涉及的是一种及冰区船舶性能预报方法,具体地说是一种层冰下船舶推力减额估计方法。
背景技术
层冰工况是极地船舶航行主要工况之一。自然条件形成的层冰是极地航行中常见的情形。极地船舶遭遇层冰时,会产生较大的冰载荷,使得极地船舶的推进性能的研究不同于常规船舶的研究方法。
在常规船舶拖曳试验中,船舶的牵引力与裸船体阻力相平衡。而在常规船舶自航试验中,船舶的总推力与裸船体阻力、船桨相互作用下阻力增值相平衡。其中阻力增值的物理起源是由于螺旋桨转动导致其前后的流体加速流动,导致边界层的剪切速率增加,船体后部的压力降低(前进模式)。类似于常规船舶自航试验,在覆冰水域航行,船舶自航条件可表示为:船舶的总推力与裸船体水阻力、裸船体冰阻力、以及船体、螺旋桨、水流和冰之间的相互作用引起的船体增阻。相比于常规船舶。冰的存在导致了额外的冰-船、冰-桨、以及冰-水之间的相互作用。其中冰-船的相互作用导致的是冰阻力的大幅度增加;冰-桨的作用表现为大块船底碎冰对螺旋桨的阻塞作用以及小块船底碎冰与桨的直接接触。而冰-水的相互作用表现为附着在船体底部的碎冰随着船体不断向前,流向船艉,导致原本不受螺旋桨影响区域的伴流场改变,从而影响推进性能。
针对层冰区极地船舶航行性能的研究,主要依赖于模型试验以及早期的实船数据,限于现阶段的数值手段,想要采用数值模拟方法实现层冰条件下对推力减额的估算十分困难,原因在于直接实现层冰条件下对推力减额的估算涉及船-冰-桨-流多相耦合,作用机理十分复杂。有研究者会直接采用常规船舶的推力减额计算方法来完成极地船舶的推力减额估算方法,误差较大。
发明内容
本发明的目的在于提供一种精度高的层冰工况船舶推力减额分数数值估计方法。
本发明的目的是这样实现的:
步骤1,输入目标船舶裸船体实尺度模型并进行网格划分、流体域设置,建立裸船体数值计算模型;设置不同航速的工况,进行数值计算,得到基于裸船体的水阻力;
步骤2,输入目标船舶自航实尺度模型并采用与步骤1相同的网格大小进行划分、相同的流体域进行设置,建立自航数值计算模型;设置不同航速、转速的工况,进行数值计算,得到自航模型的水阻力、推力以及扭矩;
步骤3,输入目标船舶带附着冰的自航实尺度模型并采用与步骤一相同的网格大小进行划分、相同的流体域进行设置,建立带附着冰的自航数值计算模型;设置不同航速、转速的工况,进行模型计算,得到层冰工况下典型时刻的带有附着冰自航模型的水阻力、推力以及扭矩;
步骤4,依据裸船体的水阻力以及自航模型的水阻力按照推力减额定义求解水中的推力减额;
步骤5,依据带有附着冰的自航模拟结果得到的有效推力以及不带附着冰的自航模拟结果得到的有效推力的差额,将该部分差额作为冰中的推力减额;
步骤6,利用冰中的推力减额和水中的推力减额得到层冰工况下不同航速、转速的推力减额分数。
为了克服采用数值方法求解层冰工况下推力减额分数困难问题,以及采用常规船舶推力减额误差较大的问题。本发明提出了一种层冰工况下考虑附着冰对伴流场的作用以及冰块对推进器的阻塞作用的船舶推力减额分数数值估算方法。为了估算极地船舶的推进减额分数,本发明是一种考虑冰对伴流场的作用、以及冰对桨的阻塞作用的层冰工况下推力减额估算方法。将冰中推力减额增值和水中推力减额增值单独求解,再将其叠加,实现层冰工况下的极地船舶推力减额的估算工作,该方法对于估算极地船舶推进性能、以及估算船舶破冰能力具有重要的意义。
本发明的特点包括:
将冰水介质下船舶的推力减额复杂问题分解为简单的冰中推力减额计算问题以及水中推力减额的计算问题,实现将船舶冰区自航模拟的复杂问题转化成带有层冰特征的常规船舶自航模拟问题。
估算层冰下船舶推力减额,将冰和水分离,采用裸船体模型、自航模型、以及带附着冰的自航模型来计算估算层冰下的船舶推力减额系数。
可结合模型层冰试验工况下,典型船-冰位置状态图片,来反演船舶推力减额数值估算工作。同时,可设定等间距的船底附着冰自航模型完成船舶推力减额估算工作。
与现有方法相比,本发明的有益效果是:通过带附着冰和不带附着冰的自航模型,考虑了冰的对于伴流场的影响以及冰对桨的阻塞作用,最终形成跟贴近于实际层冰工况下的推力减额数值估算方法。
附图说明
图1是本发明流程图。
图2a-图2c是裸船体实尺度模型图,其中:图2a为主视图;图2b为仰视图;图2c为左视图。
图3a-图3c是目标船自航实尺度模型(不带附着冰)图,其中:图3a为主视图;图3b为仰视图;图3c为左视图。
图4a-图4c是目标船自航实尺度模型(带附着冰)图,其中:图4a为主视图;图4b为仰视图;图4c为左视图。
具体实施方式
下面举例对本发明做更详细的描述。
结合图1,本发明的层冰工况船舶推力减额分数数值估计方法的具体过程是:
步骤1,提供目标船舶裸船体实尺度模型,并对模型进行网格划分、流体域设置,建立裸船体数值计算模型。设置不同航速的计算工况,开展数值计算,得到基于裸船体的水阻力。
步骤2,提供目标船舶自航实尺度模型,并对自航模型采用与步骤一相同的网格大小进行划分、相同的流体域进行设置,建立自航数值计算模型。设置不同航速、转速的计算工况,开展数值计算,得到自航模型的水阻力、推力以及扭矩。
步骤3,提供目标船舶带附着冰的自航实尺度模型,并对带附着冰的自航模型采用与步骤一相同的网格大小进行划分、相同的流体域进行设置,建立带附着冰的自航数值计算模型。设置不同航速、转速的计算工况,开展模型计算,得到层冰工况下典型时刻的带有附着冰自航模型的水阻力、推力以及扭矩。
步骤4,依据裸船体的水阻力以及自航模型的水阻力按照推力减额定义求解水中的推力减额。
步骤5,依据带有附着冰的自航模拟结果得到的有效推力以及不带附着冰的自航模拟结果得到的有效推力的差额,该部分差额产生的原因在于:附着冰对伴流场的影响、以及冰块对螺旋桨的阻塞作用。将该部分差额作为冰中的推力减额。
步骤6,利用冰中的推力减额和水中的推力减额可求得层冰工况下不同航速、转速的推力减额分数。
本发明可以采用多种商用软件实现本发明提出的推力减额估算方法,本实施例利用STAR-CCM商用软件模拟提出的层冰工况下推力减额数值估算方法。本实例带附着冰的自航模型,不限于图中形式,可依据船舶实际情况建立相应的带附着冰的模型。
一种层冰工况船舶推力减额分数数值估算方法流程图如图1所示。图2是裸船体实尺度模型。图3是目标船自航实尺度模型(不带附着冰)。图4是目标船自航实尺度模型(带附着冰)。
Claims (1)
1.一种层冰工况船舶推力减额分数数值估计方法,其特征是:
步骤1,输入目标船舶裸船体实尺度模型并进行网格划分和流体域设置,建立裸船体数值计算模型;设置不同航速的工况,进行数值计算,得到基于裸船体的水阻力;
步骤2,输入目标船舶自航实尺度模型并采用与步骤1相同的网格大小进行划分和相同的流体域进行设置,建立自航数值计算模型;设置不同航速和转速的工况,进行数值计算,得到自航模型的水阻力、推力以及扭矩;
步骤3,输入目标船舶带附着冰的自航实尺度模型并采用与步骤一相同的网格大小进行划分和相同的流体域进行设置,建立带附着冰的自航数值计算模型;设置不同航速和转速的工况,进行模型计算,得到层冰工况下典型时刻的带有附着冰自航模型的水阻力、推力以及扭矩;
步骤4,依据裸船体的水阻力以及自航模型的水阻力按照推力减额定义求解水中的推力减额;
步骤5,依据带有附着冰的自航模拟结果得到的有效推力以及不带附着冰的自航模拟结果得到的有效推力的差额,将该部分差额作为冰中的推力减额;
步骤6,利用冰中的推力减额和水中的推力减额得到层冰工况下不同航速和转速的推力减额分数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911294887.9A CN111079286B (zh) | 2019-12-16 | 2019-12-16 | 一种层冰工况船舶推力减额分数数值估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911294887.9A CN111079286B (zh) | 2019-12-16 | 2019-12-16 | 一种层冰工况船舶推力减额分数数值估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111079286A CN111079286A (zh) | 2020-04-28 |
CN111079286B true CN111079286B (zh) | 2022-06-17 |
Family
ID=70315022
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911294887.9A Active CN111079286B (zh) | 2019-12-16 | 2019-12-16 | 一种层冰工况船舶推力减额分数数值估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111079286B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111797464B (zh) * | 2020-06-05 | 2022-07-01 | 长江三峡通航管理局 | 升船机船舶-机械牵引系统设计集成系统 |
CN111783321B (zh) * | 2020-08-05 | 2021-05-11 | 中国水利水电科学研究院 | 一种缺乏资料区河渠冬季冰情发展过程的模拟方法 |
CN112034705B (zh) * | 2020-08-05 | 2022-05-03 | 智慧航海(青岛)科技有限公司 | 一种船舶动力定位时的推进器及舵控制方法 |
CN112379591B (zh) * | 2020-10-22 | 2022-05-03 | 智慧航海(青岛)科技有限公司 | 一种考虑推进器性能的推力分配优化方法 |
CN112231938A (zh) * | 2020-11-23 | 2021-01-15 | 中国船舶科学研究中心 | 一种基于数值模拟的船舶连续破冰阻力预报方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102298326A (zh) * | 2011-06-30 | 2011-12-28 | 哈尔滨工程大学 | 欠驱动auv自适应轨迹跟踪控制装置及控制方法 |
CN107902042A (zh) * | 2017-10-27 | 2018-04-13 | 华中科技大学 | 一种基于Simulink模型的船舶多轴电力推进系统 |
CN108846207A (zh) * | 2018-06-15 | 2018-11-20 | 哈尔滨工程大学 | 一种基于非线性模型辨识的船舶运动预报装置及其预报方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009021522A1 (en) * | 2007-08-14 | 2009-02-19 | Propeller Control Aps | Efficiency optimizing propeller speed control for ships |
-
2019
- 2019-12-16 CN CN201911294887.9A patent/CN111079286B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102298326A (zh) * | 2011-06-30 | 2011-12-28 | 哈尔滨工程大学 | 欠驱动auv自适应轨迹跟踪控制装置及控制方法 |
CN107902042A (zh) * | 2017-10-27 | 2018-04-13 | 华中科技大学 | 一种基于Simulink模型的船舶多轴电力推进系统 |
CN108846207A (zh) * | 2018-06-15 | 2018-11-20 | 哈尔滨工程大学 | 一种基于非线性模型辨识的船舶运动预报装置及其预报方法 |
Non-Patent Citations (4)
Title |
---|
Analysis of waterjet-hull interaction and its impact on the propulsion performance of a four-waterjet-propelled ship;Jie Gonga等;《Ocean Engineering》;20190401;第211-222页 * |
Evaluation of selected state-of-the-art methods for ship transit simulation in various ice conditions based on full-scale measurement;Fang Li等;《Cold Regions Science and Technology》;20180315;第151卷;第94-108页 * |
三桨船水动力性能试验研究;覃冲宇;《江西建材》;20150413(第161期);第206-207页 * |
碎冰条件下冰区船阻力性能试验研究;郭春雨等;《哈尔滨工程大学学报》;20160127;第37卷(第4期);第481-486页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111079286A (zh) | 2020-04-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111079286B (zh) | 一种层冰工况船舶推力减额分数数值估计方法 | |
Choi et al. | Resistance and propulsion characteristics of various commercial ships based on CFD results | |
Phillips et al. | Evaluation of manoeuvring coefficients of a self-propelled ship using a blade element momentum propeller model coupled to a Reynolds averaged Navier Stokes flow solver | |
Ponkratov et al. | Validation of ship scale CFD self-propulsion simulation by the direct comparison with sea trials results | |
KR101863746B1 (ko) | 운항선의 표준 운항 상태 선속-동력 해석 시스템 | |
Palmer et al. | Modelling tunnel thrusters for autonomous underwater vehicles | |
KR101859671B1 (ko) | 운항선의 표준운항상태 선속-동력 해석 방법 | |
Choi et al. | A numerical study of scale effects on performance of a tractor type podded propeller | |
KR102434772B1 (ko) | 선박의 시마진 도출 시스템 및 방법, 동 방법을 컴퓨터에서 실행하기 위한 컴퓨터 프로그램이 기록된, 컴퓨터 판독 가능한 기록 매체 | |
Zou et al. | Confined water effects on the viscous flow around a tanker with propeller and rudder | |
Nakashima et al. | A basic study for propulsive performance prediction of a cascade of wing sails considering their aerodynamic interaction | |
Wu et al. | Modeling and simulation of automatic berthing based on bow and stern thruster assist for unmanned surface vehicle | |
Bensow | Simulation of unsteady propeller blade loads using OpenFOAM | |
Gornicz et al. | Estimating berthing performance of twin shaft ship by means of CFD | |
Göksu et al. | Performance prediction Of DTMB 4119 and highly skewed model propellers using star-Ccm+ Cfd solver | |
Tang et al. | Performance simulation of a dual-direction ship in level ice | |
De Barros et al. | Progress towards a method for predicting AUV derivatives | |
Duman et al. | URANS approach in hull-propeller-rudder interaction of a surface combatant at high speed | |
Zhang et al. | Hydrodynamic performance of a self-propelled KCS at angle of drift including rudder forces | |
Sigmund et al. | Numerical prediction of the propulsion characteristics of ships in waves | |
CN111814410B (zh) | 一种细长体出水空泡与碎冰相互作用的预报方法 | |
CN102508960B (zh) | 变吃水作业状态下耙吸挖泥船动力定位的滤波器设计方法 | |
Yamaguchi et al. | A study on shape optimization for an underwater vehicle based on numerical simulation | |
Park et al. | CFD Analysis of marine propeller-hub vortex control device interaction | |
CN114444207B (zh) | 一种基于实尺度数值模拟的船舶航速预报方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |