三维多浮体水动力学问题中波浪交互理论的参数应用特性探究分析

  • 格式:pdf
  • 大小:1.08 MB
  • 文档页数:8

第28卷第7期 

2012年7月 科技通报 

BULI TIN OF SCIENCE AND TECHN0L0GY Vo1.28 No.7 July 2012 

三维多浮体水动力学问题中波浪交互理论的 

参数应用特性探究分析 

史琪琪 ,柏木正 ,杨建民 

(1.上海交通大学船舶海洋与建筑工程学院,上海200240;2 ,崔进举 

日本大阪大学,日本大阪567—0058) 

摘要:波浪交互理论是解决j维多浮体水动力问题的一种有效方法,该方法能正确处理相邻浮体造 

成的波浪干涉效应,尤其适用于浮体数目较多的场合。文中运用波浪交互理论,对不同间距下一个箱型 

多浮体系统的水动力问题进行求解,研究其中的关键参数一消逝波模数和傅立叶级数展开项数的变 

化,对浮体表面压力分布和波浪力等水动力结果所产生的影响,为波浪交互理论在三维多浮体水动力 

问题上的进一步深入应用提供依据。 

关键词:波浪交互理论;消逝波模数;傅立叶级数展开项;多浮体;水动力 

中图分类号:035 文献标识码:A 文章编号:1001—71 19(2012)07—0001—07 

Research on Application Characteristics of Parameters in 

Wave-interaction Theory when Solving Hydrodynamic 

Problems of 3一D Multiple Floating Bodies 

SHI Qiqi ,Masashi Kashiwagi 2,YANG Jianmin ,CUI j均u l 

(1.School of Naval Architecture and Civil Engineering,Shanghai Jiao Tong Unive ̄ity,Shanghai 200240,China; 

2.Department of Naval Architecture and Ocean Engineering,Osaka University,Osaka 567—0058,Japan) 

Abstract:he wave interaction theory is an effective way tO solve hydrodynamic problems for 3-D muhiple floating bodies,by 

correedy dealing with the wave interferential effects.Thus,it is especiMly suitable to be adopted in the case when the number 

of floating bodies becomes large.In this paper,the wave interaction theory is applied to solve the hydrodynamic problem of 

four identical box-shaped floating bodies as the separation distance between bodies varies,and the influence of two critical 

Dmmeters,the numbers of evanescent wave modes and of Fourier series in the —direction,on the results,such as the 

pressure distribution on the body surface and the wave exciting force,is studied.All of these will help tlS go further in solving 

3-D hydrodynamic problems of multiple floating bodies by the wave interaction theory. Key words:volatility;wave interaction theory,evanescent wave modes,numbers of Fourier series,multiple floating bodies, 

hydrodynamics 

0 引言 

在处理三维多浮体系统的水动力学问题时, 由于各浮体之间的相互作用对结果影响显著,分 

析时必须对此着重考虑。波浪交互理论是解决三 

维多浮体系统水动力问题的一种数值方法,该方 

收稿日期:2011一Ol一25 基金项目:国家自然科学基金资助项目(50879045) 作者简介:(1986一),女,浙江宁波人,上海交通大学硕士研究生,主要从事海洋工程水动力学方面的研究, 

Email:graceshqq@yahoo.com.cn。 通讯作者:杨建民,教授;Email:jmyang@sjtu_edu.en。

 2 第28卷 

法先考虑单个浮体的水动力效应,通过对速度势 

的变换进一步考虑由其他浮体产生的波浪干涉 

效应,进而求解整体的速度势问题。该方法计算 

效率高,对于浮体数目较大的问题尤为适用,能 

有效改善高阶边界元法等方法中计算时间随浮 

体个数增加而急剧增大的不足。 

波浪交互理论由Ohkusu[ ]首次运用到水波问 

题上,之后Kagemoto、Yue及其他人对此理论进 

行了拓展和延伸闭;该理论在解决多浮体问题时 

十分有效,Murai等还曾将该理论用于求解VLFS 

的水动力实例问题中[31。 

本文运用波浪交互理论,对4个对称分布的 

相同箱式浮体组成的系统的水动力问题进行求 

解,研究不同的浮体间距情况下,该理论中的两 

个关键性参数,消逝波模数n和0方向上傅立叶 

级数展开项数m的变化,对多浮体表面压力分布 

和波浪力等水动力结果产生的影响,探索其变化 

规律,为波浪交互理论在三维多浮体水动力问题 

上的进一步深入应用提供一定的依据。 

1 原理概述 

根据势流理论,流体假设为无旋、无粘和不 

可压缩,速度势满足拉普拉斯方程。流场中的总 

速度势可以分解为人射波速度势、散射势和辐射 

势 。其中, 和 s合称绕射势 。。 

, J= ,£ , J+ , J (1) 

本文中只考虑绕射势问题,由于假定结构物 

在平衡位置周围作微幅的简谐振荡,可将绕射势 

转化为定常问题求解。 

瞩,t)=-Re[ { + J}e ] (2) 

其中,g为重力加速度, 为人射波幅值, 

c J= J+ 为定常绕射势。 

定常绕射势必须满足Laplace方程、自由面、 

物面、和水底条件,运用格林定理分析上述问题, 

可以得到如下积分方程。 

cf ,Pj+qbo(Q)OG(P ̄Q)_d5=咖 ,Pj (3) 

其中,尸= , 和Q= y ,z 分别为流场中 

的点和浮体表面上的积分点,c 为立体角,G 

为满足自由表面条件和辐射条件的自由表面格 林函数。 

1.1高阶边界元法 

通过高阶边界元法可求解上述绕射势的积 

分方程,该方法将浮体表面离散成多个四边形单 

元,浮体表面坐标和对应速度势均可由对应的二 

阶形状函数表示。由迭代计算,可将边界积分方 

程离散成速度势的数值方程组。浮体自由表面内 

部设置若干附加的场点,以消除不规则频率的影 

响,采用最小二乘法求解相应结点上的绕射势l 5_。 

一旦得到绕射势,浮体各表面上的无因次波浪力 

可通过积分求得,进而获得系统的整体波浪力。 

1.2波浪交互理论 

由波浪交互理论求解绕射势积分方程,入射 

波可由局部柱坐标系( )下的单元入射波系 

列表示,满足拉普拉斯方程和各边界条件[61。 

相应于单元入射波,绕射势 D= 亦满 

足原积分方程。其中 为每个单元入射波对应 

的散射势。 

当流场中的点在远离浮体表面的外部区域 

fr>rI 1时,Bessel函数的附加定理成立。代人由格 

林定理推得的速度势中,可得流场中单元入射波 

对应的散射势, 

{ }_ } (4) 

其中,[ j-{ 枷, )为散射势的系数矩阵, 

单元散射势向量{ l定义如下, 

{ }:z0 (5)J-/rr ̄ i j=… 胎 e 【)) 

其中,m=O,±1,±2, 为0方向上傅立叶级数 

的展开项数; =1,2,…,∞为消逝波模数。 

不考虑水波相互作用时的入射波速度势 ,, 

可由第i个浮体柱坐标系下的 ,表示, 

,=={ IT{ } (6) 

对于第i个浮体而言,入射波分为两部分:一部分 

是由上式给出的入射波,另一部分是由周围其他 

浮体产生的散射势。 

Ne qb { { ̄iIt+∑{AJsf { 

=( lr+∑{ASsIT[Ti,])t@,/ 

其中,%为浮体个数,{A }为第 个浮体散 

射势表达式中的未知系数向量,[ 为{ }和{ } 

之间的坐标转换矩阵。 

因此,第i个浮体的散射势可以类似表示为,

 第7期 史琪琪等.三维多浮体水动力学问题中波浪交互理论的参数应用特性探究分析 3 

({ Y+∑{ }.I ) 帆} (8) 

芦 根据定义,此处的散射势与(7)式定义的 

A }应保持相等,因此可以得到联立的线性 

方程组。 

N { }一 ∑[珊{A }_吲 

for i=1,2,…,% (9) 

求解(9)式,可以得到系数{/4 },进而得到幅 

N8 值向量({ai} +∑/aJs} 'tTA)。 

j=l ̄/#i 则浮体表面的无因次压力分布,可由整体绕 

射势 表达如下, 

咖 咖 

% 

=({ }_r+∑IAJs}T[T ̄J)( }+ }(10) 

≠1 % 

=({ai}.r+∑IAis}T[TjJ){0 ̄ot 

芦 其中,{ }表示第i个浮体上相应于单元入 

射波,不考虑波浪相互作用的绕射势。 

作用在各表面上的波浪力积分可表示为, 

m { } dS 

其中, 为对应积分面, 为波浪力模态( = 

1,2,3分别表示纵向、横向和垂向), 为浮体表 

面法向量的k方向分量。 

将各个浮体上的波浪力积分叠加,可以得到 

多体系统的整体波浪力。 

在运用波浪交互理论分析浮体问的相互干 

涉影响时,该方法中的消逝波模数n和0方向上 

傅立叶级数展开项数m必须是有限的,因此它们 

可分别表示为n=1,2,…,Ⅳ和m=0,±1,±2,…,± 。 

即消逝波模数的总项数为Ⅳ,而0方向上傅立叶 

级数展开的总项数为2 1。总的说来,当各个浮 

体之间相距较远,且浮体为柱状时, 可以取较 

小的数值;而当水深相对较浅,浮体间距比较大 

时,Ⅳ可取较小的数值,这是由于消逝波衰减而 

不产生波浪干涉效应。 

在保证结果精确度的前提下,当消逝波模数 

和0方向上傅立叶级数展开项数可取较小值时, 

波浪交互理论的计算效率将得到很大的改善。因 

此适当地选取这两个参数的值,将是保证波浪交 互理论准确有效应用的关键所在。 

2 研究对象 

本文以四个对称布置的相同方型浮体组为 

例,在不同的浮体间距下,分别研究了波浪交互 

理论中消逝波模数n的变化,以及0方向 傅立 

叶级数展开项数m的变化,对三维多浮体水动力 

问题分析结果的影响规律,为该理论的进一步深 

入应用提供依据。模型布置如下图所示。 

Y 

图1计算模型布置 

Fig.1 Arrangement of the computed models