CN109325309B - 船舶大幅横摇运动的三维数值模拟方法 - Google Patents

船舶大幅横摇运动的三维数值模拟方法 Download PDF

Info

Publication number
CN109325309B
CN109325309B CN201811236527.9A CN201811236527A CN109325309B CN 109325309 B CN109325309 B CN 109325309B CN 201811236527 A CN201811236527 A CN 201811236527A CN 109325309 B CN109325309 B CN 109325309B
Authority
CN
China
Prior art keywords
ship
wave
motion
calculating
amplitude
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
Application number
CN201811236527.9A
Other languages
English (en)
Other versions
CN109325309A (zh
Inventor
陈纪康
段文洋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Engineering University
Original Assignee
Harbin Engineering University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN201811236527.9A priority Critical patent/CN109325309B/zh
Publication of CN109325309A publication Critical patent/CN109325309A/zh
Application granted granted Critical
Publication of CN109325309B publication Critical patent/CN109325309B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/13Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Architecture (AREA)
  • Civil Engineering (AREA)
  • Structural Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)

Abstract

本发明提供的是一种船舶大幅横摇运动的三维数值模拟方法,包括:读取网格文件,进行船舶静水力计算;计算泰勒展开边界元方法所需的边界积分方程所涉及的影响系数矩阵;叠模速度势及其空间一、二阶导数和Mj项求解;时域格林函数及其空间法向导数求解;计算横摇阻尼系数;不规则波分解,线性叠加成入射波时历;泰勒展开边界元法直接时域扰动波浪力计算;入射波浪力与静水回复力计算;大幅运动预报方程建模,采用四阶龙格库塔方法步进求解运动方程,进行船舶在顶浪或者斜浪中非线性运动评估;进行不规则波中船舶大幅运动数值模拟和特征统计。利用本发明的方法能够预报集装箱船迎浪规则波中大幅运动及不规则波中船舶大幅运动数值模拟和特征统计。

Description

船舶大幅横摇运动的三维数值模拟方法
技术领域
本发明涉及的是一种模拟方法,特别涉及一种利用三维泰勒展开边界元方法模拟船舶大幅横摇运动的方法。
背景技术
首尾变化剧烈的船舶,在特定的海浪工况下容易发生参数横摇现象,即大幅横摇运动。特别是集装箱船极容易发生该现象。船舶一旦发生参数横摇现象则会对货物带来损失,对船舶安全带来影响。因此船舶参数横摇现象发生概率统计是水动力学研究的热点。
船舶运动预报数值模拟中需要考虑定常叠模势对非定常速度势的影响。因此自由面条件实施对数值预报精度也有影响。目前常用的自由面条件包括NK自由面条件和DB自由面条件。因为DB自由面条件设计叠模势二阶导数的计算,因此要比NK自由面条件要复杂。而尖角边界处速度势二阶导数精确求解是数值模拟的难点。
集装箱船舶尾部均为大外飘船型,利用时域格林函数方法求解存在数值发散问题。因此对于大外飘船型运动时域预报均采用简单格林函数法。对于远方辐射条件实施有阻尼区法和时域格林函数法。相比于阻尼区法,时域格林函数方法大大减少自由面网格数目,能够提高数值模拟的效率。
发明内容
本发明的目的在于提供一种能够为船舶参数横摇现象发生概率统计提供基础的船舶大幅横摇运动的三维数值模拟方法。
本发明的目的是这样实现的:
步骤1,读取网格文件,利用网格信息进行船舶静水力计算;
步骤2,计算泰勒展开边界元方法所需的边界积分方程所涉及的影响系数矩阵;
步骤3,叠模速度势及其空间一阶、二阶导数和Mj项求解;
步骤4,时域格林函数及其空间法向导数求解;
步骤5,计算横摇阻尼系数;
步骤6,不规则波分解,线性叠加成入射波时历;
步骤7,泰勒展开边界元法直接时域扰动波浪力计算;
步骤8,基于弱散射理论入射波浪力(Froude-Krylov力)与静水回复力计算;
步骤9,船舶大幅运动预报方程建模,采用四阶龙格库塔方法步进求解运动方程,进行船舶在顶浪或者斜浪中非线性运动评估;
步骤10,依据船舶大幅运动响应RAO,进行不规则波中船舶大幅运动数值模拟和特征统计。
集装箱船在船首尾部,舭部均存在尖角。因此精确求解尖角边界处叠模势,非定常扰动速度势切向诱导速度是数值模拟的关键。本发明是一种面向船舶参数横摇直接时域模拟的泰勒展开边界元方法,利用泰勒展开边界元方法能够精确求解任意流域边界的二阶导数。弱散射理论可以考虑瞬时波面下入射波浪力和恢复力,考虑船舶大幅运动。结合这两点本发明提出了面向船舶参数横摇直接时域模拟的泰勒展开边界元方法。
利用本发明提出的方法能够预报集装箱船迎浪规则波中大幅运动RAO及不规则波中船舶大幅运动数值模拟和特征统计。
附图说明
图1是本发明的流程图。
具体实施方式
下面举例对本发明做更详细的描述。
1)读取网格文件自动提取船舶水线信息(包括首尾驻点空间坐标,水线拟合曲线函数)。结合流场匹配边界智慧参数和水线信息,自动生成符合边界元方法的水面和匹配面离散网格。并基于船体离散网格,计算船舶排水体积,浮心,漂心,惯性矩,湿表面积等静水力参数。检查静水力参数数值结果与物理船舶参数的误差,以此检验船舶网格质量。
2)本发明利用泰勒展开边界元求解各速度势成分及其空间一阶、二阶导数。泰勒展开核心思想是基于格林第三公式形成的边界积分方程进行数值离散求解的方法。对于三维问题,将浮体湿表面离散为若干四边形或三角形单元,在每一单元上,取单元节点坐标均值为中心,在面元中点对偶极强度作泰勒展开并保留一阶导数项,对源强作泰勒展开只保留一阶导数项。并引入场点的切向一阶导数来封闭方程组,从而构成了关于偶极强度、偶强的一阶导数为未知数,源强为已知变量的线性代数方程组。其中偶强切向一阶导数的影响系数包含主值(归一化后为正/负二分之一)。上述操作方法产生的利用边界单元求解边界积分方程数值解的方法称为泰勒展开边界元方法。对于任意面元i可得到如下简化的一阶泰勒展开边界元方法的离散方程组,i=1,2,…,N,
Figure BDA0001838287860000021
上式各矩阵中元素表达式:
Figure BDA0001838287860000022
Figure BDA0001838287860000031
Figure BDA0001838287860000032
式中:上角标i和j表示面元编号。以矩阵中某一元素表达式为例做一解释:如:
Figure BDA0001838287860000033
该方法可同时求解速度势及两个相互正交的切向方向导数。在利用物面法向不可穿透条件,即构成了当地局部坐标系的速度场。可实现速度场在局部坐标下及大地坐标系下的转换。数值结果证明,该方法可明显改善流域边界拐角处的切向诱导速度的计算精度。引入辅助函数
Figure BDA00018382878600000311
再次利用泰勒展开边界元方法求解φ的空间,即速度势空间二阶导数。因此该方法涉及到若干影响系数矩阵计算。
3)总速度势可分解为定常速度势、非定常入射势、辐射速度势和绕射速度势。即:
Figure BDA0001838287860000034
定常速度势又可分解为来流速度势和叠模速度势,即:
Φb=-Ux+Φ (3)
非定常势定解问题的物面条件涉及定常势影响,即mj项。对于DB假设,其边值定解问题为:
Figure BDA0001838287860000035
基于DB线性假设mj项为:
Figure BDA0001838287860000036
式中:
Figure BDA0001838287860000037
为物面各点处的位移。
Figure BDA0001838287860000038
为船体平动位移,
Figure BDA0001838287860000039
为船体转动位移,
Figure BDA00018382878600000310
为船体湿表面上各点位置矢径。
4)本专利利用匹配法外传扰动波浪,在匹配面上利用时域格林函数构造匹配条件。时域格林函数为:
Figure BDA0001838287860000041
5)横摇阻尼计算
利用Ikeda方法计算船舶横摇阻尼系数,注意此时得到的总阻尼为无因次等效线性阻尼系数。该方法中,按如下方式无因次化:
Figure BDA0001838287860000042
其中ρ为水密度,
Figure BDA0001838287860000043
为排水体积,B为型宽,等效线性阻尼系数可表示为:
Figure BDA0001838287860000044
其中B441为线性阻尼系数,B443为立方阻尼系数,η4为横摇幅值,ωη为横摇固有频率。
为了更准确的求解非线性阻尼系数,本发明分别求出不同横摇幅值下的等效线性阻尼系数,然后根据上式进行最小二乘拟合,即可求得以上非线性阻尼系数。
6)将波浪谱进行傅里叶变换,得到对应波浪谱的子波参数;本方法将输入的波浪时历利用ITTC双参数谱的波谱公式通过有义波高和周期表示如下:
Figure BDA0001838287860000045
其中,T1为谱心周期;H1/3为有义波高,ω为圆频率,S(ω)为海浪谱密度。
利用等分能量法,将谱能量分为num份,各组成波的波幅如下:
Figure 1
其中ω1~ω2为截断频率,利用线性叠加原理即可得到入射波浪的时历表达式:
Figure BDA0001838287860000047
其中εi为各组成波的相位,而且
Figure BDA0001838287860000048
7)本发明采用时域直接求解非定常扰动速度势,辐射速度势和绕射速度势一起求解。非定常扰动速度势定解问题如式(9),利用泰勒展开边界元法时域步进求解该初边值问题。
Figure BDA0001838287860000051
式中自由面条件右端项F的表达式表示如下:
Figure BDA0001838287860000052
采用积分格式自由面条件步进自由面上各离散单元中心点处的速度势。以任意函数f(t) 为例,阐述积分格式自由面条件的核心思想:即对被积函数作时间二次积分。
Figure BDA0001838287860000053
同理对自由面条件作时间二次积分得:
Figure BDA0001838287860000054
利用梯形法计算积分格式自由面条件。从而实现自由面条件时间步进。一旦扰动速度势求解完后,利用伯努利方程在平均湿表面上积分,可得到扰动波浪载荷。
8)瞬时波面以下入射波浪力(Froude-Krylov力)与静水回复力计算;由伯努利方程,可知静水面以下入射波压力为:
Figure BDA0001838287860000055
静水压力为:
Ps=-ρgz
由于对于入射波力和静水回复力的计算积分至瞬时湿表面,故存在入射波波面在静水以上的情况,因此本发明采用Wheeler法评估静水面以上入射波压力的贡献:
PIS=-ρgz+ρgζ(t)ek(z-ζ(t)) (z≤ζ(t)) (16)
此方法保证了波浪压力随深度增加而呈指数衰减。
为了计算瞬时湿表面下静水恢复力和入射波浪力,先将面网格沿船体表面划分至甲板,采用四边形和三角形网格。在运动方程计算得到瞬时船体运动姿态后,通过转换关系将原始正浮状态的网格信息转换为瞬时姿态下的网格信息。在每个时间步判断船体上的湿表面网格在入射波面以下,以上还是与入射波面相截。瞬时网格中,水面以下的部分即为瞬时湿表面。
在获得了瞬时入射波面与船体湿表面交界面以下的瞬时湿表面后,将入射波和静水压力 PIS在瞬时船体湿表面积分,就获得了船舶受到的非线性入射波力和静水恢复力,表示如下:
Figure BDA0001838287860000061
Figure BDA0001838287860000062
在上式中,入射波和静水恢复力矩相对船舶重心的随体坐标系取矩。其原因是后续的船舶运动转动方程是在相对于船舶重心的随体坐标系下建立起来的。
9)根据船舶大幅运动预报方程建模,采用运动方程可进行船舶在顶浪中非线性运动评估。依据牛顿第二定律,船舶垂荡,横摇,纵摇三自由度耦合运动方程为:
Figure BDA0001838287860000063
利用四阶龙格库塔方法步进求解运动方程。
10)计算船舶大幅运动响应RAO,进行不规则波中船舶大幅运动数值模拟和特征统计。
根据船舶运动方程时域模拟,获得船舶运动响应时历,由此利用修正的周期图平均法(也称为Welch方法)获得船舶航行中遭遇的运动响应谱Si(ω)和海浪谱Sζ(ω)
Figure BDA0001838287860000064
其中L为将时历分的段数,Ii(ω)为加了海明窗后每一段的修正周期图。
假设运动响应谱和海浪谱之间满足线性传递函数,按照谱分析基本理论,可获得船舶运动响应RAO,如下式所示。
Figure BDA0001838287860000065
上式中Si(ω)代表由某一自由度响应曲线得到的频谱函数,Sζ(ω)代表波面的频谱函数。
利用逐波分析理论对模拟的运动响应时历有义值、运动极值和振荡统计进行分析。

Claims (6)

1.一种船舶大幅横摇运动的三维数值模拟方法,其特征是包括如下步骤:
步骤1,读取网格文件,利用网格信息进行船舶静水力计算;
提取船舶水线信息,结合流场匹配边界智慧参数和水线信息,生成符合边界元方法的水面和匹配面离散网格,并基于船体离散网格,所述的船舶水线信息包括首尾驻点空间坐标、水线拟合曲线函数;所述的利用网格信息进行船舶静水力计算具体包括:计算船舶静水力参数,所述的静水力参数包括排水体积、浮心、漂心、惯性矩以及湿表面积;检查静水力参数数值结果与物理船舶参数的误差;
步骤2,计算泰勒展开边界元方法所需的边界积分方程所涉及的影响系数矩阵;
将浮体湿表面离散为若干四边形或三角形单元,在每一单元上,取单元节点坐标均值为中心,在面元中点对偶极强度作泰勒展开并保留一阶导数项,对源强作泰勒展开只保留一阶导数项,引入场点的切向一阶导数来封闭方程组,构成关于偶极强度、偶强的一阶导数为未知数,源强为已知变量的线性代数方程组;
步骤3,叠模速度势及其空间一阶、二阶导数和Mj项求解;
步骤4,时域格林函数及其空间法向导数求解;
步骤5,计算横摇阻尼系数;
步骤6,不规则波分解,线性叠加成入射波时历;
步骤7,泰勒展开边界元法直接时域扰动波浪力计算;
步骤8,基于弱散射理论入射波浪力与静水回复力计算;
步骤9,船舶大幅运动预报方程建模,采用四阶龙格库塔方法步进求解运动方程,进行船舶在顶浪或者斜浪中非线性运动评估;
步骤10,依据船舶大幅运动响应RAO,进行不规则波中船舶大幅运动数值模拟和特征统计。
2.根据权利要求1所述的船舶大幅横摇运动的三维数值模拟方法,其特征是所述的时域格林函数及其空间法向导数求解是利用匹配法外传扰动波浪,在匹配面上利用时域格林函数构造匹配条件。
3.根据权利要求1所述的船舶大幅横摇运动的三维数值模拟方法,其特征是所述的计算横摇阻尼系数具体包括:利用Ikeda方法计算船舶横摇阻尼系数,按如下方式无因次化,
Figure FDA0003885040050000011
其中ρ为水密度,
Figure FDA0003885040050000012
为排水体积,B为型宽,等效线性阻尼系数表示为:
Figure FDA0003885040050000021
其中B441为线性阻尼系数,B443为立方阻尼系数,η4为横摇幅值,ωη为横摇固有频率,
分别求出不同横摇幅值下的等效线性阻尼系数,然后根据上式进行最小二乘拟合。
4.根据权利要求1所述的船舶大幅横摇运动的三维数值模拟方法,其特征是所述的不规则波分解,线性叠加成入射波时历具体包括:
将波浪谱进行傅里叶变换,得到对应波浪谱的子波参数,将输入的波浪时历利用ITTC双参数谱的波谱公式通过有义波高和周期表示如下:
Figure FDA0003885040050000022
其中,T1为谱心周期;H1/3为有义波高,ω为圆频率,S(ω)为海浪谱密度,
利用等分能量法,将谱能量分为num份,各组成波的波幅如下:
Figure FDA0003885040050000023
其中ω1~ω2为截断频率,利用线性叠加原理得到入射波浪的时历表达式:
Figure FDA0003885040050000024
其中εi为各组成波的相位,而且
Figure FDA0003885040050000025
5.根据权利要求1所述的船舶大幅横摇运动的三维数值模拟方法,其特征是所述的泰勒展开边界元法直接时域扰动波浪力计算具体包括:采用时域直接求解非定常扰动速度势、辐射速度势和绕射速度势一起求解。
6.根据权利要求1所述的船舶大幅横摇运动的三维数值模拟方法,其特征是所述的基于弱散射理论入射波浪力与静水回复力计算具体包括:
采用Wheeler法评估静水面以上入射波压力的贡献,将面网格沿船体表面划分至甲板,采用四边形和三角形网格,在运动方程计算得到瞬时船体运动姿态后,通过转换关系将原始正浮状态的网格信息转换为瞬时姿态下的网格信息,在每个时间步判断船体上的湿表面网格在入射波面以下、以上还是与入射波面相截,瞬时网格中,水面以下的部分即为瞬时湿表面;
在获得了瞬时入射波面与船体湿表面交界面以下的瞬时湿表面后,将入射波和静水压力在瞬时船体湿表面积分,得到船舶受到的非线性入射波力和静水恢复力。
CN201811236527.9A 2018-10-23 2018-10-23 船舶大幅横摇运动的三维数值模拟方法 Active CN109325309B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811236527.9A CN109325309B (zh) 2018-10-23 2018-10-23 船舶大幅横摇运动的三维数值模拟方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811236527.9A CN109325309B (zh) 2018-10-23 2018-10-23 船舶大幅横摇运动的三维数值模拟方法

Publications (2)

Publication Number Publication Date
CN109325309A CN109325309A (zh) 2019-02-12
CN109325309B true CN109325309B (zh) 2022-12-13

Family

ID=65263159

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811236527.9A Active CN109325309B (zh) 2018-10-23 2018-10-23 船舶大幅横摇运动的三维数值模拟方法

Country Status (1)

Country Link
CN (1) CN109325309B (zh)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109902427B (zh) * 2019-03-12 2022-07-15 哈尔滨工程大学 一种船舶横摇阻尼粘流模拟流域自动化网格划分方法
CN110083882B (zh) * 2019-04-04 2022-09-23 河海大学 一种晃动水槽对波浪形态影响的模拟方法
CN110115840B (zh) * 2019-05-16 2021-10-29 腾讯科技(深圳)有限公司 虚拟场景中的物体移动控制方法、装置、终端及存储介质
CN110717216B (zh) * 2019-08-30 2023-05-23 南京航空航天大学 不规则波下带柔性气囊直升机横摇响应预报方法
CN110728023B (zh) * 2019-09-06 2023-11-21 久瓴(江苏)数字智能科技有限公司 老虎窗构件节点生成方法、装置、计算机设备和存储介质
CN110765583B (zh) * 2019-09-28 2021-06-22 成都理工大学 基于三次Bezier的分段曲线拟合方法
CN111708081B (zh) * 2020-05-29 2022-04-15 成都理工大学 考虑衰减频散的深度域地震记录合成方法
CN112307419B (zh) * 2020-10-19 2022-05-24 华南理工大学 一种船舶垂向运动响应的快速预报方法
CN113111603B (zh) * 2021-04-07 2022-07-15 哈尔滨工程大学 一种双浮体平台波浪激励力及运动响应预报方法
CN114925458B (zh) * 2022-06-06 2023-05-23 中国船舶科学研究中心 一种平整冰中船舶回转数值模拟方法
CN115495919B (zh) * 2022-09-30 2023-05-26 中国船舶科学研究中心 一种基于格子玻尔兹曼的时域缓坡方程的数值求解方法
CN115526052B (zh) * 2022-10-08 2023-04-28 南京航空航天大学 基于等分频率-能量-方向角策略不规则波浪谱解析方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010031399A1 (en) * 2008-09-19 2010-03-25 Decision3 Sp/F System for dynamically optimizing the operation of a ship
CN107600352A (zh) * 2017-08-22 2018-01-19 哈尔滨工程大学 一种基于船舶水动力在线预报的船舶减摇控制系统
CN108549616A (zh) * 2018-02-28 2018-09-18 哈尔滨工程大学 基于g-n波浪模型的船舶在恶劣海况中大幅运动的预报方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010031399A1 (en) * 2008-09-19 2010-03-25 Decision3 Sp/F System for dynamically optimizing the operation of a ship
CN107600352A (zh) * 2017-08-22 2018-01-19 哈尔滨工程大学 一种基于船舶水动力在线预报的船舶减摇控制系统
CN108549616A (zh) * 2018-02-28 2018-09-18 哈尔滨工程大学 基于g-n波浪模型的船舶在恶劣海况中大幅运动的预报方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Side Wall Effects on the Hydrodynamics of a Floating Body by Image Green Function Based on TEBEM;Jikang Chen 等;《Journal of Marine Science and Application》;20171207;第17卷(第3期);第353-361页 *
基于泰勒展开边界元法和镜像格林函数的浮体池壁效应水动力计算;王隶加,等;《第十四届全国水动力学学术会议暨第二十八届全国水动力学研讨会》;20170808;第113-119页 *

Also Published As

Publication number Publication date
CN109325309A (zh) 2019-02-12

Similar Documents

Publication Publication Date Title
CN109325309B (zh) 船舶大幅横摇运动的三维数值模拟方法
CN109472065B (zh) 预报恶劣海况下低速肥大型船舶波浪增阻的三维数值方法
CN108549616B (zh) 基于g-n波浪模型的船舶在恶劣海况中大幅运动的预报方法
CN109446634B (zh) 基于泰勒展开边界元方法的船舶运动预报方法
Zhang et al. Study on prediction methods and characteristics of ship underwater radiated noise within full frequency
Song et al. Modelling of water wave interaction with multiple cylinders of arbitrary shape
Houtani et al. Generation of a spatially periodic directional wave field in a rectangular wave basin based on higher-order spectral simulation
CN109344531A (zh) 预报多浮体结构物波漂载荷的三维频域数值方法
CN109376427B (zh) 考虑池壁效应的船舶波浪增阻的三维数值方法
Kim et al. Fatigue life prediction of ship and offshore structures under wide-banded non-Gaussian random loadings: Part II: extension to wide-banded non-Gaussian random processes
Yang et al. The impact of coastal reclamation on tidal and storm surge level in Sanmen Bay, China
Ning et al. Experimental and numerical investigation of the hydrodynamic characteristics of submerged breakwaters in waves
CN103778282A (zh) 一种基于最大横流累积频率的航道有效宽度设计方法
CN106528924B (zh) 一种应用于侧壁式气垫船的湿甲板砰击预报方法
Nokob et al. Diffraction and radiation loads on open cylinders of thin and arbitrary shapes
Wei et al. Integrated approach for estimating extreme hydrodynamic loads on elevated pile cap foundation using environmental contour of simulated typhoon wave, current, and surge conditions
Yim et al. Nonlinear ocean wave models and laboratory simulation of high seastates and rogue waves
Ghafari et al. Potential-based boundary element method to calculate the hydrodynamic drift force on the floating cylinder
Du et al. A combined first-and second-order theory for the deckwetness prediction of sandglass-type floating body in irregular head waves
Priyanto et al. Probability distribution of wave run up and dynamic response on a large volume semi-submersible
Giorgi et al. NLFK4ALL: An open-source demonstration toolbox for computationally efficient nonlinear Froude-Krylov force calculations
Street et al. Modeling the wave climate in the Baltic sea
Duan et al. PTEBEM for wave drift forces based on hydrodynamic pressure integration
Wang A second order random wave model for predicting the power performances of a wave energy converter
Chen et al. Side Wall Effects on the Hydrodynamics of a Floating Body by Image Green Function Based on TEBEM: Keynote Contribution for the International Workshop on Wave Loads and Motions of Ships and Offshore Structures, Harbin, China, 5-7 November, 2017

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