船舶力学

船舶有航速时三维时域格林函数计算方法研究

PDF下载

鲁江 ¹ , 张楠 ¹ , 张新曙 ² , 顾民 ¹

船舶力学 | 流体力学 2024,28(1): 20-35

收起

船舶力学 | 流体力学 2024, 28(1): 20-35

船舶有航速时三维时域格林函数计算方法研究

全屏

鲁江¹, 张楠¹, 张新曙², 顾民¹

作者信息

^1.中国船舶科学研究中心，江苏　无锡　214082

^2.上海交通大学　海洋工程国家重点实验室，上海　200240

鲁江（1980-），男，博士，研究员，通讯作者，E-mail：lujiang1980@aliyun.com

张楠（1977-），男，博士，研究员。

通讯作者:

通讯作者，E-mail：lujiang1980@aliyun.com

Method of three-dimensional time-domain Green function for ships at forward speed

Jiang LU¹, Nan ZHANG¹, Xin-shu ZHANG², Min GU¹

Affiliations

^1.China Ship Scientific Research Center, Wuxi 214082, China

^2.State Key Laboratory of Ocean Engineering, Shanghai Jiao Tong University, Shanghai 200240, China

出版时间: 2024-01-20 doi: 10.3969/j.issn.1007-7294.2024.01.003

文章导航

摘要

收起

船舶工业CAE软件势流求解器的开发对势流计算方法的可靠性提出要求，本文针对船舶有航速时深水三维时域格林函数计算方法开展研究。首先引入三维时域自由面格林函数，其瞬时项的求解采用Hess & Smith方法，自由面记忆项的求解采用美国密歇根大学Beck团队的方法。其次，分区域推导出三维时域自由面格林函数记忆项及其导数的编程所用的数学展开表达式。然后，采用间接法模型求解分布源强和速度势，并采用压力直接积分法求出时域辐射力。最后，针对Wigley I船型开展了辐射力计算，并与公开的试验及计算结果进行对比，验证本文方法的可靠性。

关键词

三维时域格林函数 / 辐射力 / 船舶CAE自主软件

Abstract

收起

Developing a reliable potential flow solver is necessary for the ship CAE software. The method on the three-dimensional time-domain Green function of ships with forward speed in infinite depth was studied. Firstly, a three-dimensional transient free-surface Green function was introduced, and its Rankine part was calculated by Hess & Smith’s method, and its free-surface memory part was calculated by the method of Beck team from the University of Michigan, US. Secondly, the memory part of the Green function and its derivatives were obtained in the mathematical expressions for making the program in five regions. Thirdly, the source method was used to calculate the source and the radiation potential, and then the radiation force was obtained by integrating the pressure around the floating body. Finally, the radiation force of the Wigley I ship was calculated, and compared with the published experimental and numerical results. The reliability of the method and the code in this paper is confirmed.

Key words

three-dimensional transient free-surface Green function / radiation force / ship CAE software

引用本文

鲁江, 张楠, 张新曙, 顾民. 船舶有航速时三维时域格林函数计算方法研究. 船舶力学, 2024 , 28 (1) : 20 -35 . DOI: 10.3969/j.issn.1007-7294.2024.01.003

Jiang LU, Nan ZHANG, Xin-shu ZHANG, Min GU. Method of three-dimensional time-domain Green function for ships at forward speed[J]. Journal of Ship Mechanics, 2024 , 28 (1) : 20 -35 . DOI: 10.3969/j.issn.1007-7294.2024.01.003

正文

收起

0　引言

收起

船舶CAE软件是国产船舶行业工业软件体系的重要组成部分，对基于三维时域面元法的势流求解器开发提出了需求。在上个世纪80年代和90年代，基于边界元法的三维时域面元法得到了快速发展，但局限于计算机计算速度，主要停留在算法验证阶段，工程应用上仍以频域为主。2000年以来，随着计算机软硬件的快速发展，三维时域面元法的工程应用得到快速发展。2020年以来，船舶工业界把基于势流求解器和粘流求解器的国产船舶流体CAE软件开发提上了议程。三维时域面元法是直接在时域内建立初边值问题进行求解，在处理瞬态问题（抨击、冲击作用）、非线性失稳大幅运动、波浪中机动操纵等问题时具有频域方法不可替代的优势，但三维时域边界积分方程，不仅包含面积分和线积分，而且还包含整个历程的时间卷积计算，计算量巨大。由于三维时域面元法的工程应用优点以及近代计算机的快速发展，基于三维时域面元法的势流求解器开发重新引起船舶界的重视。

Finkelstain（1957）^[1]提出了线性自由面下的三维时域格林函数；Wehausen和Laitone（1960）^[5]给出了有航速深水三维时域格林函数积分形式；Cummins（1962）^[2]和Ogilvie（1964）^[3]首先讨论了非定常运动问题的时域直接求解方法；Wehausen（1967）^[4]详细讨论了零航速时非定常运动问题，给出了零航速三维时域问题的Haskind关系。目前，针对该格林函数主要有三种求解方法，介绍如下。

美国密歇根大学Beck团队的Liapis（1986）^[6]、King（1987）^[7]、Magee（1991）^[8]研究了船舶有航速时的三维时域格林函数法，提出一种三维时域格林函数计算方法，根据μ和β值的范围将时域格林函数分成（μ <0.7， β <10）、（μ <0.7， β ≥10）、（0.7≤μ <0.98）、（0.98≤μ ≤1， β <10）、（0.98≤μ ≤1， β ≥10）5个区域，每个区域用不同的级数展开和渐进展开进行计算。

Newman（1985，1990）^[9-10]提出满足10E-6和10E-10计算精度的两套三维时域格林函数计算方法，Newman团队的Bingham（1994）^[11]采用三维时域面元法进一步验证了该三维时域格林函数方法解决Newman-Kelvin线性运动的有效性。Lin和Yue（1990）^[12]基于Newman（1985）^[9]的三维时域格林函数计算方法，根据μ和β值的范围将时域格林函数分成5个区域，每个区域用不同的级数展开、渐进展开和剩余区间Chebyshev正交逼近进行计算，依据该格林函数的计算方法，Lin和Yue（1999）等^[13]进一步提出一种满足线性自由面条件的数值方法，将流场分为内场和外场，以内场Rankine源方法为核心，外场匹配时域格林函数来满足辐射条件，即混合源法。以此为基础，美国海军开发了船舶大幅运动评估软件LAMP，分为LAMP1线性版本、LAMP2弱非线性版本和LAMP4全非线性版本（Shin，2003）^[14]。

黄德波（1992）^[15]导出了三维时域格林函数及其导数的简化计算公式，对于β>8+1.515μ采用渐进展开，对于β<8+1.515μ使用级数展开，提出造表插值方案和相应的节点划分办法，三维时域格林函数法开始在国内得到学者的广泛应用；周正全（1993）等^[16]提出了反向辐射势表示绕射力的关系式，可认为是Wehausen（1967）^[4]的无航速时域Haskind关系的发展，称之为有航速时的Haskind关系，其三维时域格林函数的计算采用黄德波（1992）^[15]方法；王大云（1996）^[17]提出了一种三维水弹性时域分析方法，其三维时域格林函数的计算采用黄德波（1992）^[15]方法；卜淑霞等（2018）^[18]和储纪龙（2019）等^[19]采用三维时域混合源法研究了参数横摇现象，其三维时域格林函数的计算采用Newman（1985）^[9]和Lin和Yue（1990）^[12]的计算方法；陈纪康、段文洋等（2021）^[20]采用三维时域混合源法构建边界积分方程预报船舶运动，其外域三维时域格林函数及其导数采用黄德波（1992）^[15]的计算方法，并利用泰勒展开边界元法求解该边界积分方程，解决了有航速定解问题中涉及的非光滑边界处切向诱导速度精度低的问题。

Clement（1998）^[21]提出了三维时域格林函数及其导数的四阶微分方程方法；Duan和Dai（2001）^[22]、申亮和朱仁传等（2007）^[23]分别对三维时域格林函数及其导数的四阶微分方程方法做了推导和部分改进；Li、Ren等（2015）^[24]提出了三维时域格林函数及其导数的四阶微分方程精细积分方法，对Clement （1998）^[21]的四阶微分方程法做出了显著改进；杨鹏（2016）^[25]给出了一种三维时域水弹性运动方程的数值求解方法，其三维时域格林函数的计算采用Clement（1998）^[21]的四阶微分方程法；周文俊等（2021）^[26]提出了基于垂向积分形式时域格林函数的多域高阶面元法，其三维时域格林函数的计算采用Clement （1998）^[21]的四阶微分方程法。

Newman（1990）^[10]指出Beck团队对三维时域格林函数计算方法做出了实质性改进，但国内公开应用的很少，本文依据美国密歇根大学Beck团队的三维时域自由面格林函数理论方法，推导出可直接用于程序开发的三维时域自由面格林函数记忆项完整数学展开表达式，推导给出可直接用于程序开发的三维时域自由面格林函数记忆项导数的完整数学展开表达式，给出可编程应用的参数取值和求解说明，给出了本文计算结果和公开的试验结果、计算结果以及作者使用的加强切片法计算结果的对比验证，使得三维时域势流求解器的工程开发简单化，可促进波浪中操纵性、快速性和稳性等学科的发展。

1　数学模型

收起

空间固定坐标系O₀-X₀Y₀Z₀原点位于静水面，X₀轴为波浪传播方向，Y₀轴指向左为正，Z₀轴垂直水面向上为正。参考坐标系o-xyz随着船体以恒定速度U₀沿着x正方向前进，初始时刻和空间坐标系重合，不随船体转动而转动。随船坐标系o-x′y′z′固定于船体上，原点和参考坐标相同，船体正浮平衡时，与参考坐标系重合，随船体的转动而转动。坐标系如图1所示。设入射波浪向角为β，船舶航向角为χ，则β =-χ。

1.1　辐射力时域直接积分法

在势流理论里总的速度势Φ_T可写成如下表达式：

(1)

式中，-U₀x+Φ₀（x，y，z）为定常兴波速度势，ϕ₀（x，y，z，t）为入射波速度势，ϕ_i（x，y，z，t）为浮体运动产生的辐射速度势，i =1，2，⋯，6分别对应纵荡、横荡、垂荡、横摇、纵摇和首摇，ϕ₇（x，y，z，t）为绕射速度势，t为时间。

根据Bernoulli方程，作用于浮体湿表面上的流体压力可以表示为

(2)

式中，p（x，y，z，t）为t时刻坐标（x，y，z）处的压力，S_B（t）为t时刻浮体湿表面，ρ为液体密度，g为重力加速度。

根据压力和力/力矩的关系，忽略二阶速度势小量，由船体i模态运动引起的作用在j自由度方向上辐射力/力矩可以表示为

(3)

式中，j =1，2，⋯，6分别对应纵荡、横荡、垂荡、横摇、纵摇和首摇6个方向。

为了消除船体物面速度势空间导数的计算，许多学者采用Tuck^[27]理论假定。公式（3）的第二项是关于速度势的前进速度影响项，适用于分部积分法，该积分项在船首到船尾整船积分为零，公式（3）可以写成下面表达式：

(4)

式中，n_j为浮体湿表面上在j方向矢量，n₁是小量，n₂、n₃在x方向变化很慢，

。

式（4）离散成可编程的数学表达式：

(5)

式中，M是浮体湿表面上面元的个数，A_p是第p个面元的面积，[n_j]_p是第p个面元上的n_j，[m_j]_p是第p个面元上的m_j。

在势流理论里，辐射速度势ϕ_i （x，y，z，t）物面边界条件满足如下表达式，本文计算采用船体静水中湿表面。

(6)

式中，x_i （t），

分别是t时刻浮体i方向模态运动的位移和速度。

时域辐射力求解的关键是求出辐射速度势ϕ_i，下述章节将论述辐射速度势的求解方法。

1.2　三维时域自由面格林函数法

在势流理论里，时域求解非定常速度势的方法主要有时域自由面格林函数法和时域Rankine源法。自由面格林函数法满足线性自由面条件以及除物面条件外的其他所有边界条件，仅需在浮体湿表面上分布奇点，数值离散化的方程只含有浮体湿表面及水线上的未知量。

Wehausen和Laitone（1960）^[5]给出了深水有航速三维时域格林函数积分形式，Liapis（1986）^[6]和King（1987）^[7]采用脉冲函数法。三维时域格林函数积分形式可写成如下表达式：

(7)

式中，格林函数G（P，Q；t-τ）表示场点P（x，y，z）在t时刻由于源点Q（ξ，η，ζ）在时间τ时刻的脉冲扰动而产生的速度势；δ（t）为脉冲函数，

；H（t）为阶跃函数，

；J₀（）为第一类零阶Bessel函数，G₀（P，Q）和

分别表示三维时域格林函数中瞬时项和自由面记忆项，R表示源点Q（ξ，η，ζ）和场点P（x，y，z）之间的水平距离。

1.3　三维时域自由面格林函数及其导数计算方法

三维时域自由面格林函数的瞬时项采用Hess&Smith方法，参见文献[28]，不再论述。下面主要描述本文中三维时域自由面格林函数的记忆项具体采用的计算方法及推导给出的编程数学表达式。

对三维时域自由面格林函数记忆项表达式作如下无因次处理（King，1987）^[7]：

(8)

其中：

(8a)

(8b)

(8c)

(8d)

是求解三维时域自由面格林函数记忆项的关键，对公式（8）等量替换可得如下表达式：

(9)

三维时域自由面格林函数记忆项对空间导数的无因次表达式如下，和（King，1987）^[7]一致：

(10)

根据公式（8b），推导出μ对x，y，z的导数如下：

(11)

根据公式（8c），推导出β对x，y，z的导数如下：

(12)

根据公式（10），推导出

、

后，乘以对应的

，

即可求出三维时域自由面格林函数记忆项的空间导数。

参考文献[7]，本文时域格林函数记忆项求解分为五个区域，第一、二区域的分界线β为8.8，第三、四区域的分界线为0.97，并推导了可直接用于程序开发的三维时域格林函数记忆项及其导数的数学展开表达式。

第一区域，0≤μ<0.7，0≤β<8.8，时域格林函数记忆项采用Legendre多项式的级数展开，文献[6]给出了公式的前4项，从软件开发角度，推导给出可编程的数学表达式如下：

(13)

根据公式（13），可进一步推导出该区域对应的三维时域格林函数记忆项导数的数学表达式为

(14)

n的取值以计算结果收敛到计算精度为准，文献[7]没有给出具体值，n大于80时已满足计算精度10E-10的要求。

第二区域，0≤μ <0.7，β ≥8.8，采用渐进展开，格林函数公式分成两部分来论述。

(15)

文献[7]在第一项中给出了前3个表达式，从软件开发角度，推导给出该区域第一项可编程的数学表达式为

(16)

n的取值以计算结果收敛为准，文献[7]没有给出具体值，n大于100时已满足计算精度10E-6的要求。

该区域格林函数第二项表达式（King，1987）^[7]如下：

(17)

式中，θ =sin^-1 （μ）。

根据公式（16），从软件开发角度，推导给出该区域对应的三维时域格林函数记忆项第一项导数的数学表达式为

(18)

根据公式（17），从软件开发角度，推导给出该区域格林函数记忆项公式第二项第一部分导数的可编程数学表达式如下，其他部分导数推导原理类同。

(19)

第三区域，0.7≤μ <0.97， β ≥0，采用递推的Filon积分，格林函数公式分成两部分来论述。

(20)

该区域格林函数的第二项

及其导数和第二区域

相同，此处不再给出。文献[7]中

误写为2，需要注意。下面主要论述该区域格林函数第一部分。

根据文献[7]，推导出该区域格林函数第一部分的可编程数学表达式如下：

(21)

(21a)

公式中n的取值范围根据积分上限k_max来确定，k_max的确定是一个关键点。其中，

，

，δ =γh。

通过零阶贝塞尔函数J₀（）可求出f （k），J₀ （k²）表达式如下。n的取值以计算结果收敛为准，文献[7]没有给出具体值，n大于100时已满足计算精度10E-6的要求。

(22)

进一步通过公式

来确定k_max。本文设定区间k =0~30，步长为0.05，找出积分上限k_max。

文献[7]定义（γh）_min为小值，但没有给出具体值。本文（γh）_min的值设为0.35，Δk=h=0.05，n的值为k_max/（2h）取整数。

根据公式（21），推导给出该区域格林函数记忆项的第一项导数的可编程数学表达式：

(23)

(24)

(25)

(26)

第四区域，0.97≤μ <1， β <10，采用贝塞尔函数^[7]展开。

(27)

式中，a₀ =1， a₁ =-2.249 999 7， a₂ =1.265 620 8， a₃ =-0.316 386 6， a₄ =0.044 447 9， a₅ =0.000 210 0，

，

，w（）为误差函数。

按照公式（27），推导出循环公式（28），求出每一个I_m的实部和虚部，即可求出Im{I_{4n + 2}}

(28)

公式（27）中n最大值为5，公式（28）中对应的m最大值为22。

根据公式（27）和公式（28），推导给出该区域格林函数记忆项的导数的可编程数学表达式：

(29)

当μ ≈1时，格林函数记忆项及其导数表达式如下：

(30)

第五区域，0.97≤μ <1， β ≥10，采用误差函数的渐进展开（King，1987）^[7]，推导给出可编程的数学表达式：

(31)

根据公式（31），推导给出该区域格林函数记忆项公式导数的可编程数学表达式为

(32)

(33)

n的取值以计算结果收敛为准，文献[7]没有给出具体值，n大于30时已满足计算精度10E-6的要求。

当μ ≈1时，格林函数记忆项及其导数表达式如下：

(34)

1.4　间接法求解源密度速度势

同时在物面上布置源和偶极子的方法称为直接法，仅布置源的方法称为间接法。间接法可以通过边界积分方程得到船体湿表面的切向速度，更适合求解考虑瞬时湿表面的非线性大幅运动问题。文献[8]的公式（2.11）、（2.12）基于势流理论给出了如下求解分布源密度σ（P，t）的边界积分方程和速度势ϕ（P，t）表达式：

(35)

注：文献[8]的公式（2.11）水线积分前的符号误写为“-”，其离散公式（3.6）是正确的。

(36)

根据公式（36），本文用到的速度势在x方向的导数公式如下，y，z方向的导数类似：

(37)

式中：S_B（t）为t时刻浮体湿表面；S_B（τ）为τ时刻浮体湿表面；Γ（τ）为τ时刻水面与船体湿表面的交线；

为浮体湿表面面元的单位法向矢量，方向为指向浮体湿表面、离开流体域；

是自由面上单位法向量，方向为指向Γ（τ）、离开流体域；V_N是Γ（τ）在

方向的法向速度，V_N和V_n的关系为

，V_N可以近似取线元临近面元的法向速度在该线元切向速度的投影。

公式（35）~（36）离散写成可编程的数学表达式如下：

(38)

(39)

式中：M为浮体湿表面上面元个数，q为浮体湿表面上第q个面元，p为浮体湿表面上第p个面元；L为水线上的线元个数，l为水线上第l个线元；t_N表示当前的时间步N对应的时间，t_n表示开始计及记忆效应的第n步对应的时间，Δt为时间步长。

2　数值计算方法验证

收起

2.1　三维时域自由面格林函数记忆项及其导数计算及验证

图2给出了无量纲化格林函数记忆项及其导数值。图3中“Present”表示本文计算结果，“Yang，2016”表示文献Yang（2016）^[25]的结果，后续图中表示方法类似。本文计算结果和Yang（2016）^[25]采用Clement（1998）的微分方程法计算的结果吻合。

从图中可以看出，在μ接近0时（场点P和源点Q靠近水线），随时间β的增大，格林函数及其导数的振荡特性非常明显，能够影响时域波浪力求解的稳定性，是三维时域格林函数计算方法需要不断改进的原因之一。随着μ变大，格林函数及其导数值随时间β的增大而减小，最后趋于零，μ越大，其趋于零的速度越快。

2.2　Wigley I船型辐射力计算及验证

Wigley Ι船型型值参见文献[29]，本文船长设为9.81 m，船宽设为0.981 m，吃水设为0.613 125 m。ω为摇荡频率，L为船长，Δ为排水体积。图4和图5中“Present”表示本文计算结果，“Exp_Journee_1992”表示文献[29]的试验结果，“Bingham_1994”表示文献[11]的计算结果，“Estrip”表示本文采用加强切片法计算的结果。

通过三维时域格林函数计算Wigley I船型的垂荡、纵摇时域辐射力，根据时域辐射力曲线的振幅及其和波浪的相对相位，整理对应附加质量和阻尼系数（Present），并和公开的试验结果（Journée，1992）^[29]、计算结果（Bingham，1994）^[11]以及作者采用加强切片法（Kashiwagi et al，2010）^[30]计算的结果（EStrip）进行对比。本文网格数设为240，时间步长设为0.05 s，垂荡方向时域辐射力计算采用公式（3），纵摇方向时域辐射力计算采用公式（4）。图4给出了强制垂荡运动引起的垂向附加质量、阻尼系数以及在纵摇方向的附加质量和阻尼系数，图5给出了强制纵摇运动引起的纵摇方向的附加质量、阻尼系数以及垂荡方向附加质量和阻尼系数。本文计算结果和（Journée，1992）^[29]试验结果、（Bingham，1994）^[11]计算结果相比，整体吻合较好，但偏差还是存在。和（Journée，1992）^[29]试验结果出现偏差的原因主要在于真实流体是黏性的，强制摇荡试验时还有自由面影响，而计算是基于理想的理论假定，且只考虑平均湿表面。和（Bingham，1994）^[11]计算结果出现偏差的原因主要是计算资源的限制，网格、时间步长取值不一致，水线的计算处理不完全一致。作者同时采用加强切片法（EStrip）计算了水动力系数，和切片法计算结果相比，三维时域方法整体上优于切片理论的计算结果。后续采用三维时域混合源法进一步验证三维时域理论的有效性。

2.3　Wigley I船型强制运动产生的时域辐射力

通过三维时域格林函数进一步计算Wigley I船型在不同航速时的时域辐射力，图6给出了强制垂荡/纵摇振幅为0.05 m/0.05 rad，频率为3 rad/s，航速Fn=0.0、0.3、0.6时在垂荡/纵摇方向产生的时域力/力矩。尽管计算是基于静水湿表面，但力的求解公式（4）和物面条件公式（6）都含有和航速有关的m项，同时时域格林函数中场点和源点的距离和航速有关。航速增大后，该方法仍可以计算出稳定的垂荡和纵摇方向的时域辐射力/力矩，只是力/力矩振幅增大，相位变化微小。根据时域辐射力/力矩曲线的振幅及其和波浪的相位差，可以整理出对应的附加质量和阻尼系数。

3　结论

收起

本文依据美国密歇根大学Beck团队的三维时域自由面格林函数理论公式，推导出可直接用于程序开发的三维时域自由面格林函数记忆项及其导数的数学展开表达式，通过Wigley I船型验证了本文方法的可靠性，得出如下结论：

（1）本文采用的方法能够准确得出三维时域自由面格林函数记忆项及其导数值，重现格林函数及其导数的振荡特性。

（2）本文采用的三维时域自由面格林函数法能够可靠地求解Wigley I船型的垂荡、纵摇及其耦合方向的辐射力。

（3）本文推导给出了三维时域自由面格林函数记忆项及其导数的编程数学展开表达式，有助于从事此方面研究的人员进行编程应用，实用性强。

本文方法和代码可用于船舶工业CAE软件势流求解器，实船工程应用上需要在规避三维时域格林函数的振荡特性以及外飘船型、尾部浅吃水的速度势发散方面做深入研究。

基金

收起

工信部高技术船舶项目(2017（614）)

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

[1]

Finkelstein

A B

. The initial value problem for transient water waves[J]. Communications on Pure and Applied Mathematics, 1957, 10 (4): 511-522.

[2]

Cummins

W E

. The impulse response function and ship motions[J]. Schiffstechnik, 1962, 9: 101-109.

[3]

Ogilvie

T F

. Recent progress toward the understanding and prediction of ship motions[C]//Proceeding of 5^th Symposium on Naval Hydrodynamics, ONR, Washington, D.C., 1964.

[4]

Wehausen

J V

. Initial-value problem for the motion in an undulating sea of a body with fixed equilibrium position[J]. Journal of Engineering Mathematics, 1967, 1: 1-19.

[5]

Wehausen

J V

, Laitone

E V

. Surface waves[M]. Berlin: Handbuch der Physic, Springer-Verlag, 1960.

[6]

Liapis

S J

. Time-domain analysis of ship motions[D]. The University of Michigan, 1986.

[7]

King

. Time-domain analysis of wave exciting forces on ships and bodies[D]. The University of Michigan, 1987.

[8]

Magee

A R

. Large-amplitude ship motions in the time domain[D]. The University of Michigan, 1991.

[9]

Newman

J N

. The evaluation of free surface Green functions[C]//4th International Conference on Numerical Ship Hydrodynamics, Washington D C, 1985.

[10]

Newman

J N

. The approximation of free surface Green functions[C]//Proceeding of the Fritz Ursell Retirement Meeting, Manchester, England, 1990: 107-135.

[11]

Bingham

H B

. Simulating ship motions in the time domain[D]. Massachusetts Institute of Technology, 1994.

[12]

Lin

W M

, Yue

D K P

. Numerical solution for large-amplitude ship motions in time domain[C]//18th Symposium on Navel Hydrodynamics, Ann Arbor, Michigan, 1990.

[13]

Lin

W M

, Zhang

, Weems

, Yue

D K P

. A mixed source formulation for nonlinear ship-motion and wave-load simulations[C]//7th International Conference on Numerical Ship Hydrodynamics, Paris, France, 1999.

[14]

Shin

Y S

, Belenky

V L

, Lin

W M

, et al. Nonlinear time domain simulation technology for seakeeping and wave-load analysis for modern ship design[R]. ABS Technolgy Papeers, 2003.

[15]

黄德波，时域Green函数及其导数的数值计算[J]. 中国造船，1992，119(4): 16-25.

Huang

Debo

. Approximation of time-domain free surface function and its spatial derivatives[J]. Shipbuilding of China, 1992, 119(4): 16-25. (in Chinese)

[16]

周正全，张亮，戴遗山，船舶在波浪中航行时绕射问题的线性时域解[J]. 中国造船，1993，3(122): 1-25.

Zhou

Zhengquan

, Zhang

Liang

, Dai

Yishan

. Time-domain analysis of wave diffraction for ship motions with forward speed[J]. Shipbuilding of China, 1993, 3(122): 1-25. (in Chinese)

[17]

王大云. 三维船舶水弹性力学的线性时域分析方法[D]. 无锡：中国船舶科学研究中心，1996. (in Chinese)

Wang

Dayun

. Three-dimensional hydroelastic analysis of ships in time domain[D]. Wuxi: China Ship Scientific Research Center, 1996. (in Chinese)

[18]

卜淑霞，鲁江，顾民，等. 基于三维时域混合源法的顶浪不规则波参数横摇研究[J]. 船舶力学，2018，22(8): 926-934.

Shuxia

, Lu

Jiang

, Gu

Min

, et al. A study on parametric roll in irregular head waves by a 3D time domain hybrid panel method[J]. Journal of Ship Mechanics, 2018, 22(8): 926-934. (in Chinese)

[19]

储纪龙，顾民，鲁江，卜淑霞. 规则波中参数横摇直接稳性评估方法研究[J]. 水动力学研究与进展，2019，34(1): 39-44.

Chu

Jilong

, Gu

Min

, Lu

Jiang

, Bu

Shuxia

. Study on the direct stability assessment method of parametric roll in regular waves[J]. Journal of Hydrodynamics, 2019, 34(1): 39-44. (in Chinese)

[20]

Chen

J K

, Duan

W Y

, Ma

, Li

J D

. Time domain TEBEM method of ship motion in waves with forward speed by using impulse response function formulation[J]. Ocean Engineering, 2021, 227: 108617.

[21]

Clement

A H

. An ordinary differential equation for the Green function of time-domain free-surface hydrodynamics[J]. Journal of Engineering Mathematics, 1998, 33: 201-217.

[22]

Duan

W Y

, Dai

Y S

. New derivation of ordinary different equation for time-domain free-surface Green functions[J]. China Ocean Engineering, 2001, 15(4): 499-507.

[23]

申亮，朱仁传，缪国平，刘应中. 深水时域格林函数的实用数值计算[J]. 水动力学研究与进展（A辑），2007，22(3): 380-386.

Shen

Liang

, Zhu

Renchuan

, Miao

Guoping

, Liu

Yingzhong

. A practical numerical method for deep water time-domain Green function[J]. Journal of Hydrodynamics, 2007, 22(3): 380-386. (in Chinese)

[24]

Z F

, Ren

H L

, Tong

X W

, Li

. A precise computation method of transient free surface Green function[J]. Ocean Engineering, 2015(105): 318-326.

[25]

杨鹏. 船舶三维时域非线性水弹性响应研究[D]. 无锡：中国船舶科学研究中心，2016.

Yang

Peng

. 3D nonlinear hydroelastic response study of ships in time domain[D]. Wuxi: China Ship Scientific Research Center, 2016. (in Chinese)

[26]

周文俊，彭营豪，陆屿，陈曦. 基于垂向积分形式时域格林函数的多域高阶面元法及其在船舶波阻增加计算中的应用[J]. 中国造船，2021，62(4): 15-34.

Zhou

W J

, Peng

Y H

, Lu

, Chen

. Multi-domain higher-order boundary element method based on vertical integral time domain Green’s function and its application in calculation of added wave resistance[J]. Shipbuilding of China, 2021, 62(4):15-34. (in Chinese)

[27]

Ogilvie

T F

, Tuck

E O

. A rational strip theory for ship motions, Part 1[R]. The University of Michigan, 1969.

[28]

Hess

J L

, Smith

A M O

. Calculation of nonlifting potential flow about arbitrary three-dimensional bodies[J]. Journal of Ship Research, 1964, 8(2): 22-24.

[29]

Journée

J M J

. Experiments and calculations on 4 Wigley hull forms in head waves[R]. Delft University of Technology, 1992.

[30]

Kashiwagi

, Ikeda

, Sasagawa

. Effect of forward speed of a ship on added resistance in waves[J]. International Journal of Offshore and Polar Engineering, 2010, 20(2): 1-8.

2024年第28卷第1期

PDF下载

引用本文

BibTeX

文章信息

doi: 10.3969/j.issn.1007-7294.2024.01.003

接收时间：2023-07-22
首发时间：2026-03-21
出版时间：2024-01-20

补充材料

相关文章

文章信息

作者

出版历史

收稿日期：2023-07-22

基金

工信部高技术船舶项目(2017（614）)

作者信息

^1.中国船舶科学研究中心，江苏　无锡　214082

^2.上海交通大学　海洋工程国家重点实验室，上海　200240

通讯作者:

通讯作者，E-mail：lujiang1980@aliyun.com

参考文献

分享链接

https://castjournals.cast.org.cn/joweb/cblx/CN/10.3969/j.issn.1007-7294.2024.01.003

分享至

全文二维码

扫描看全文

引用本文

BibTeX

本文的引用情况

2种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT