CN109376427B - 考虑池壁效应的船舶波浪增阻的三维数值方法 - Google Patents
考虑池壁效应的船舶波浪增阻的三维数值方法 Download PDFInfo
- Publication number
- CN109376427B CN109376427B CN201811236877.5A CN201811236877A CN109376427B CN 109376427 B CN109376427 B CN 109376427B CN 201811236877 A CN201811236877 A CN 201811236877A CN 109376427 B CN109376427 B CN 109376427B
- Authority
- CN
- China
- Prior art keywords
- ship
- calculating
- motion
- wave
- equation
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- 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
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Aviation & Aerospace Engineering (AREA)
- Automation & Control Theory (AREA)
- Navigation (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供的是一种考虑池壁效应的船舶波浪增阻的三维数值方法,包括:读取网格文件,进行船舶静水力计算;计算泰勒展开边界元方法所需的边界积分方程所涉及的影响系数矩阵;叠模速度势及其空间一、二阶导数和Mj项求解;计算池壁镜像源产生的影响系数矩阵;积分格式自由面条件时域步进自由面离散网格中心点处的速度势;泰勒展开边界元法直接时域扰动波浪力计算;船舶大幅运动预报方程建模,求解运动方程;船舶波浪增阻计算;依据船舶运动响应RAO和增阻RAO,进行不规则波中船舶大幅运动和波浪增阻谱分析计算。利用本发明的方法能够预报计及池壁效应的船舶运动RAO,波浪增阻RAO,以及不规则波中船舶运动和波浪增阻谱分析结果。
Description
技术领域
本发明涉及的是一种数值模拟方法,特别涉及一种利用三维泰勒展开边界元方法计算考虑池壁效应的低航速肥大型船舶波浪增阻的数值方法。
背景技术
船舶波浪增阻值是评估船舶CO2排放等级的关键指标,在新船设计开发过程中为满足IMO国际海事组织的船舶CO2排放指数,以及运力过剩的情况下。船舶设计人员选择降低设计航速,从而降低主机功率。以便满足船舶CO2排放指数。对于油船、散货船等肥大型船舶主机储备功率不足,会导致恶劣海况下船舶安全性能降低。因此IMO国际海事组织正在讨论制定船舶最小装机功率评估导则。导则核心即如何评估恶劣海况下船舶增阻对航向稳定性的影响。
船舶设计人员依然在拖曳水池中开展低航速船舶波浪增阻试验研究。低航速船舶会在船前方形成扰动波,因此船模与池壁之间存在干扰,会对最终的波浪增阻测量有影响,因此需要通过数值模拟提前分析共振频率点。
船舶运动预报数值模拟中需要考虑定常叠模势对非定常速度势的影响。因此自由面条件实施对数值预报精度也有影响。目前常用的自由面条件包括NK自由面条件和DB自由面条件。因为DB自由面条件设计叠模势二阶导数的计算,因此要比NK自由面条件要复杂。而尖角边界处速度势二阶导数精确求解是数值模拟的难点。IMO临时导则规定8级海况。因此需要考虑船舶在恶劣海况中大幅运动对波浪增阻的影响,这是数值预报精度的关键点。
发明内容
本发明的目的在于提供一种能够为物理试验工况选择提供准则的考虑池壁效应的船舶波浪增阻的三维数值方法。
本发明的目的是这样实现的:
步骤1,读取网格文件,利用网格信息进行船舶静水力计算;
步骤2,计算泰勒展开边界元方法所需的边界积分方程所涉及的影响系数矩阵;
步骤3,叠模速度势及其空间一阶、二阶导数和Mj项求解;
步骤4,计算池壁镜像源产生的影响系数矩阵;与步骤2类似,区别是源点空间位置不同,其空间位置取决于池壁宽度;
步骤5,积分格式自由面条件时域步进自由面离散网格中心点处的速度势,并且在远场利用阻尼区,抑制扰动波浪返回流场,影响数值结果;
步骤6,泰勒展开边界元法直接时域扰动波浪力计算;
步骤7,船舶大幅运动预报方程建模,采用四阶龙格库塔方法步进求解运动方程;
步骤8,船舶波浪增阻计算;
步骤9,依据船舶运动响应RAO和增阻RAO,进行不规则波中船舶大幅运动和波浪增阻谱分析计算。
本发明提供了一种考虑池壁效应的肥大型船舶(油船、散货船)波浪增阻数值方法,寻找共振频率点,为物理试验工况选择提供准则。
本发明利用泰勒展开边界元方法精确求解定常和非定常速度势的二阶导数。能够精确预报船舶波浪增阻值。利用镜像格林函数法构造边界积分方程,模拟池壁对船舶运动及波浪增阻的影响。
利用本发明提出的方法能够预报计及池壁效应的船舶运动RAO,波浪增阻RAO,以及不规则波中船舶运动和波浪增阻谱分析结果。
附图说明
图1是本发明的流程图。
图2是船舶与镜像虚体间的坐标变换。
具体实施方式
下面举例对本发明做更详细的描述。
1)读取网格文件自动提取船舶水线信息(包括首尾驻点空间坐标,水线拟合曲线函数)。结合流场匹配边界智慧参数和水线信息,自动生成符合边界元方法的水面离散网格。并基于船体离散网格,计算船舶排水体积,浮心,漂心,惯性矩,湿表面积等静水力参数。检查静水力参数数值结果与物理船舶参数的误差,以此检验船舶网格质量。
2)本发明利用泰勒展开边界元求解各速度势成分及其空间一阶、二阶导数。泰勒展开核心思想是基于格林第三公式形成的边界积分方程进行数值离散求解的方法。对于三维问题,将浮体湿表面离散为若干四边形或三角形单元,在每一单元上,取单元节点坐标均值为中心,在面元中点对偶极强度作泰勒展开并保留一阶导数项,对源强作泰勒展开只保留一阶导数项。并引入场点的切向一阶导数来封闭方程组,从而构成了关于偶极强度、偶强的一阶导数为未知数,源强为已知变量的线性代数方程组。其中偶强切向一阶导数的影响系数包含主值(归一化后为正/负二分之一)。上述操作方法产生的利用边界单元求解边界积分方程数值解的方法称为泰勒展开边界元方法。对于任意面元i可得到如下简化的一阶泰勒展开边界元方法的离散方程组,i=1,2,…,N,
上式各矩阵中元素表达式:
该方法可同时求解速度势及两个相互正交的切向方向导数。在利用物面法向不可穿透条件,即构成了当地局部坐标系的速度场。可实现速度场在局部坐标下及大地坐标系下的转换。数值结果证明,该方法可明显改善流域边界拐角处的切向诱导速度的计算精度。引入辅助函数再次利用泰勒展开边界元方法求解φ的空间,即速度势空间二阶导数。因此该方法涉及到若干影响系数矩阵计算。
3)总速度势可分解为定常速度势、非定常入射势、辐射速度势和绕射速度势。即:
定常速度势又可分解为来流速度势和叠模速度势,即:
Φb=-Ux+Φ (3)
非定常势定解问题的物面条件涉及定常势影响,即mj项。对于DB假设,其边值定解问题为:
基于DB线性假设mj项为:
4)为满足池壁处边界条件,需对简单格林函数修正得到镜像格林函数。该主要方法为:将无限个简单格林函数关于池壁的镜像累加得到满足池壁条件的格林函数,镜像格林函数G可表示为:
式中:为满足简单格林函数;下标i表示第i个镜像虚体对本体的格林函数,当i=0时,则表示本体对本体影响的格林函数。沿着池壁方向,船舶本体上场点、源点对应在镜像虚体上的坐标可表示为pn(x,yn,z)、qn(ξ,ηn,ζ),如图2所示。镜像虚体上的源点对船舶本体上场点存在影响,其关于坐标变换可表示如下:
故而,镜像虚体上的源点可统一表示为:
ηn=nb+(-1)nη
同样地,船舶上的场点对应的镜像点pn可表示为:
yn=nb+(-1)ny (8)
利用上述镜像虚体和本体间场点和源点坐标转换关系,得到相应的镜像格林函数,并计算镜像格林函数的影响系数。
5)本专利程序采用时域直接求解非定常扰动速度势,辐射速度势和绕射速度势一起求解。非定常扰动速度势定解问题如式(9),利用泰勒展开边界元法时域步进求解该初边值问题。
本发明利用阻尼区法消除远场扰动波浪,防止扰动波浪反射至船体附近,导致数值模拟结果失真。引入阻尼区后,动力学和运动学自由面条件如下式所示:
式中μ1=3μ0(l-l0)2/L3,l是距离船舶重心纵向位置的辐射距离,阻尼区起始边缘位于l=l0,μ0和L分别为阻尼区阻尼强度及阻尼区长度。η是波面升高。本发明运用的是综合自由面条件,将上式中的波面升高消除,在转化成积分格式自由面条件,即为本发明最终所用的自由面条件。
计及阻尼区效应后,定解问题中自由面条件右端项F的表达式表示如下:
采用积分格式自由面条件步进自由面上各离散单元中心点处的速度势。以任意函数f(t)为例,阐述积分格式自由面条件的核心思想:即对被积函数作时间二次积分。
同理对自由面条件作时间二次积分得:
利用梯形法计算积分格式自由面条件。从而实现自由面条件时间步进。一旦扰动速度势求解完后,利用伯努利方程在平均湿表面上积分,可得到扰动波浪载荷。
6)在获得定常速度势,非定常扰动势后利用伯努利方程,可得到船体湿表面各离散网格中心点处的压强,对压强近场压力积分后可得到作用于船舶上的力和力矩。
7)根据船舶大幅运动预报方程建模,采用运动方程可进行船舶在顶浪中非线性运动评估。依据牛顿第二定律,船舶垂荡,横摇,纵摇三自由度耦合运动方程为:
利用四阶龙格库塔方法步进求解运动方程。
8)当速度势及其空间一阶、二阶导数计算完成,船舶六自由度运动时历信号计算完成后,利用近场压力积分公式,保留二阶压力载荷项,如下式所示即为船舶波浪增阻时历表达式。对该式取平均值即可得到船舶波浪增阻值。
9)得到船舶六自由度运动RAO和波浪增阻RAO后,利用谱分析方法能够得到各级海况下船舶运动和波浪增阻的统计值。ITTC双参数谱如下式所示:
其中,T1为谱心周期;H1/3为有义波高,ω为圆频率,S(ω)为海浪谱密度。将六自由度运动RAO值带入谱分析公式,便得到船舶在实际海况中的运动幅值。
Claims (2)
1.一种考虑池壁效应的船舶波浪增阻的三维数值方法,其特征是包括如下步骤:
步骤1,读取网格文件,利用网格信息进行船舶静水力计算;
步骤2,计算泰勒展开边界元方法所需的边界积分方程所涉及的影响系数矩阵;
所述计算泰勒展开边界元方法所需的边界积分方程所涉及的影响系数矩阵具体包括:将浮体湿表面离散为若干四边形或三角形单元,在每一单元上,取单元节点坐标均值为中心,在面元中点对偶极强度作泰勒展开并保留一阶导数项,对源强作泰勒展开只保留一阶导数项,引入场点的切向一阶导数来封闭方程组,构成关于偶极强度、偶强的一阶导数为未知数,源强为已知变量的线性代数方程组;
对于任意面元i可得到如下简化的一阶泰勒展开边界元方法的离散方程组,i=1,2,…,N,
上式各矩阵中元素表达式:
步骤3,叠模速度势及其空间一阶、二阶导数和Mj项求解;
总速度势可分解为定常速度势、非定常入射势、辐射速度势和绕射速度势;即:
定常速度势又可分解为来流速度势和叠模速度势,即:
Φb=-Ux+Φ
非定常势定解问题的物面条件涉及定常势影响,即mj项;对于DB假设,其边值定解问题为:
基于DB线性假设mj项为:
步骤4,计算池壁镜像源产生的影响系数矩阵;
将无限个简单格林函数关于池壁的镜像累加得到满足池壁条件的格林函数,镜像格林函数G表示为:
式中:为满足简单格林函数;下标i表示第i个镜像虚体对本体的格林函数,沿着池壁方向,船舶本体上场点、源点对应在镜像虚体上的坐标可表示为pn(x,yn,z)、qn(ξ,ηn,ζ),镜像虚体上的源点对船舶本体上场点存在影响,其关于坐标变换可表示如下:
镜像虚体上的源点统一表示为:
ηn=nb+(-1)nη
船舶上的场点对应的镜像点pn表示为:
yn=nb+(-1)ny
利用上述镜像虚体和本体间场点和源点坐标转换关系,得到相应的镜像格林函数,并计算镜像格林函数的影响系数;
步骤5,积分格式自由面条件时域步进自由面离散网格中心点处的速度势,并且在远场利用阻尼区,抑制扰动波浪返回流场,影响数值结果;
非定常扰动速度势定解问题如下:
引入阻尼区后,动力学和运动学自由面条件如下式所示:
式中μ1=3μ0(l-l0)2/L3,l是距离船舶重心纵向位置的辐射距离,阻尼区起始边缘位于l=l0,μ0和L分别为阻尼区阻尼强度及阻尼区长度;η是波面升高;
计及阻尼区效应后,定解问题中自由面条件右端项F的表达式表示如下:
采用积分格式自由面条件步进自由面上各离散单元中心点处的速度势;利用梯形法计算积分格式自由面条件;实现自由面条件时间步进;一旦扰动速度势求解完后,利用伯努利方程在平均湿表面上积分,得到扰动波浪载荷;
步骤6,泰勒展开边界元法直接时域扰动波浪力计算;
在获得定常速度势,非定常扰动势后利用伯努利方程,得到船体湿表面各离散网格中心点处的压强,对压强近场压力积分后得到作用于船舶上的力和力矩:
步骤7,船舶大幅运动预报方程建模,采用四阶龙格库塔方法步进求解运动方程;
船舶垂荡,横摇,纵摇三自由度耦合运动方程为:
利用四阶龙格库塔方法步进求解运动方程;
步骤8,船舶波浪增阻计算;
船舶波浪增阻时历表达式为:
步骤9,依据船舶运动响应RAO和增阻RAO,进行不规则波中船舶大幅运动和波浪增阻谱分析计算;
ITTC双参数谱如下式所示:
其中,T1为谱心周期;H1/3为有义波高,ω为圆频率,S(ω)为海浪谱密度;将六自由度运动RAO值带入谱分析公式,便得到船舶在实际海况中的运动幅值:
2.根据权利要求1所述的考虑池壁效应的船舶波浪增阻的三维数值方法,其特征是:所述的网格文件为船舶水线信息,所述船舶水线信息包括首尾驻点空间坐标、水线拟合曲线函数;所述利用网格信息进行船舶静水力计算具体包括:结合流场匹配边界智慧参数和水线信息,自动生成符合边界元方法的水面离散网格,基于船体离散网格计算船舶静水力参数,所述静水力参数包括排水体积、浮心、漂心、惯性矩、湿表面积,检查静水力参数数值结果与物理船舶参数的误差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811236877.5A CN109376427B (zh) | 2018-10-23 | 2018-10-23 | 考虑池壁效应的船舶波浪增阻的三维数值方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811236877.5A CN109376427B (zh) | 2018-10-23 | 2018-10-23 | 考虑池壁效应的船舶波浪增阻的三维数值方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109376427A CN109376427A (zh) | 2019-02-22 |
CN109376427B true CN109376427B (zh) | 2022-07-12 |
Family
ID=65401345
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811236877.5A Active CN109376427B (zh) | 2018-10-23 | 2018-10-23 | 考虑池壁效应的船舶波浪增阻的三维数值方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109376427B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110717216B (zh) * | 2019-08-30 | 2023-05-23 | 南京航空航天大学 | 不规则波下带柔性气囊直升机横摇响应预报方法 |
CN112052522B (zh) * | 2020-09-23 | 2022-09-27 | 哈尔滨工程大学 | 一种基于疲劳强度的船体结构优化简化计算方法 |
CN112182983B (zh) * | 2020-11-09 | 2023-07-25 | 中国船舶科学研究中心 | 计及海底地形及波浪影响的浮体水弹性响应分析方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103984793A (zh) * | 2014-03-21 | 2014-08-13 | 上海交通大学 | 考虑液舱晃荡影响的flng运动预报方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100430290C (zh) * | 2002-04-06 | 2008-11-05 | 程根清 | 船舶及筏的滚动式航行方法 |
CN108549616B (zh) * | 2018-02-28 | 2021-11-23 | 哈尔滨工程大学 | 基于g-n波浪模型的船舶在恶劣海况中大幅运动的预报方法 |
-
2018
- 2018-10-23 CN CN201811236877.5A patent/CN109376427B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103984793A (zh) * | 2014-03-21 | 2014-08-13 | 上海交通大学 | 考虑液舱晃荡影响的flng运动预报方法 |
Non-Patent Citations (1)
Title |
---|
三维移动脉动源格林函数的对称性及其应用;肖汶斌等;《国防科技大学学报》;20150428(第02期);138-145 * |
Also Published As
Publication number | Publication date |
---|---|
CN109376427A (zh) | 2019-02-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109325309B (zh) | 船舶大幅横摇运动的三维数值模拟方法 | |
CN109472065B (zh) | 预报恶劣海况下低速肥大型船舶波浪增阻的三维数值方法 | |
CN109376427B (zh) | 考虑池壁效应的船舶波浪增阻的三维数值方法 | |
CN109446634B (zh) | 基于泰勒展开边界元方法的船舶运动预报方法 | |
CN109344531B (zh) | 预报多浮体结构物波漂载荷的三维频域数值方法 | |
Papanikolaou et al. | Investigation into the capsizing of damaged Ro-Ro passenger ships in waves | |
CN110457795A (zh) | 基于神经网络的高层建筑物风致位移状态评估方法 | |
Xu et al. | A method for the prediction of extreme ship responses using design-event theory and computational fluid dynamics | |
Lee et al. | A study on ship performance in waves using a RANS solver, part 2: Comparison of added resistance performance in various regular and irregular waves | |
Guedes Soares et al. | Analysis of design wave loads on an FPSO accounting for abnormal waves | |
CN110083988A (zh) | 一种船舶水下辐射噪声评估方法 | |
Duan et al. | A fast time domain method for predicting of motion and excessive acceleration of a shallow draft ship in beam waves | |
CN110046459B (zh) | 一种半潜式海洋平台总体方案水下辐射噪声评估方法 | |
Wang et al. | Experimental and numerical studies of the wave-induced responses of a river-to-sea ship | |
Kim et al. | A study on prediction of whipping effect of very large container ship considering multiple sea states | |
CN106528924B (zh) | 一种应用于侧壁式气垫船的湿甲板砰击预报方法 | |
Cakici et al. | Effects of hull form parameters on seakeeping for YTU gulet series with cruiser stern | |
Lungu | Unsteady numerical simulation of the behavior of a ship moving in head sea | |
Priyanto et al. | Probability distribution of wave run up and dynamic response on a large volume semi-submersible | |
Niotis et al. | Verification of damage ship survivability with computational fluid dynamics | |
Kim et al. | Slamming impact design loads on large high speed naval craft | |
Duan et al. | Prediction of nonlinear acceleration response of a large container ship and the validation of excessive acceleration failure mode | |
Castro et al. | Calibration of a Clustering Algorithm to Improve the Nearshore Wave Prediction System | |
CN115048806B (zh) | 风浪数值模拟方法 | |
Ganjgani et al. | Investigation of Buoy Size and Location on Hydrodynamic Response and Mooring Tension of the DeepCwind Floating Wind Turbine |
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 |