三维分布源法中面元流体力的积分.pdfVIP

第 卷第 期 船舶力学 6 4 Vol . 6 No . 4 2002 年 月 8 Journal of Ship Mechanics Aug . 2002 文章编号： ( ) 1007 - 7294 2002 04 - 0024 - 09 三维分布源法中面元流体力的积分 宋竞正，段文洋，邹元杰，秦洪德 ( 哈尔滨工程大学船舶工程学院，黑龙江 哈尔滨 150001) 摘要：三维分布源法广泛应用于各种海洋结构物的波浪载荷计算。其中格林函数及其空间导数沿面元的积分和 面元上流体力的计算是十分重要的。本文首先阐述了基于四边形面元几何特性的二重高斯型求积公式。接着讨 论了如何应用这一公式计算三维分布源法中面元中心诱导速度和面元流体力积分，给出了相应的算例。最后还 导出了适合计算任意形状四边形面元上波浪主干扰力的解析表达式，并利用数值计算结果进行了验证。 关键词：高斯型求积；格林函数；波浪主干扰力( F - K 力) 中图分类号： 文献标识码： U661. 32 A 三维分布源法是计算各种大型浮体和海洋结构物波浪载荷的一种常用方法。在线性势流理论范 畴内，浮动结构物在入射波浪作用下的运动响应及波浪载荷计算可以分解为绕射问题和辐射问题。三 ( ) ( ) 维分布源法利用频域格林函数 脉动点源 ，将物体周围的三维扰动流场 绕射 辐射 速度势的描述归 / 结为物体湿表面上的分布源的积分，然后根据物面条件建立关于未知源强的积分方程。如果将物体湿 ( ) 表面离散为一系列的四边形或三角形的单元 称为面元 ，并假定分布源强在各自面元上为常数，则积 分方程便转化为代数方程组。通过求解线性方程组可以得到分布源强，进而确定扰动流场。而作用于 ( ) 物体表面上的波浪力，包括波浪主干扰力 限Frouce - Kriloff 力 、绕射力以及与运动响应有关的静水 恢复力和辐射力，将得自相应水动压力沿物体表面的积分。如何正确地计算格林函数及其空间偏导数 的面元积分，是三维分布源法中的关键所在。此外，尽管Frouce - Kriloff 力的计算是一个常规问题，但 由于它作为波浪干扰力的主要部分，其数值精确与否直接关系到浮体的运动响应，因此寻求解析表达 式计算面元波浪主干扰力是很有意义的。 ( ) 本文将在面元几何特性计算的基础上，给出适合任意形状凸四边形 包括退化为三角形 面元积 分的二重高斯型求积公式，并讨论如何应用这一求积公式来计算三维有限水深格林函数及其偏导数 的面元积分。二重高斯型求积公式还可应用于面元波浪力的计算。本文导出了面元波浪主干扰力的解 析表达式，且该公式易于在计算机程序中实现。 ! 四边形面元的几何特性 设空间四边形 取其中心 为原点，建立局部坐标系 ，则该四边形在坐标平面 A A A A