基于高阶有限元法的均匀流中物体二阶波辐射问 - 辐射研究与辐射工艺学报杂志社投稿_期刊论文发表|版面费|电话|编辑部- 辐射研究与辐射工艺学报

一、稿件要求： 1、稿件内容应该是与某一计算机类具体产品紧密相关的新闻评论、购买体验、性能详析等文章。要求稿件论点中立，论述详实，能够对读者的购买起到指导作用。文章体裁不限，字数不限。 2、稿件建议采用纯文本格式(*.txt)。如果是文本文件，请注明插图位置。插图应清晰可辨，可保存为*.jpg、*.gif格式。如使用word等编辑的文本，建议不要将图片直接嵌在word文件中，而将插图另存，并注明插图位置。 3、如果用电子邮件投稿，最好压缩后发送。 4、请使用中文的标点符号。例如句号为。而不是.。 5、来稿请注明作者署名(真实姓名、笔名)、详细地址、邮编、联系电话、E-mail地址等，以便联系。 6、我们保留对稿件的增删权。 7、我们对有一稿多投、剽窃或抄袭行为者，将保留追究由此引起的法律、经济责任的权利。二、投稿方式： 1、请使用电子邮件方式投递稿件。 2、编译的稿件，请注明出处并附带原文。 3、请按稿件内容投递到相关编辑信箱三、稿件著作权： 1、投稿人保证其向我方所投之作品是其本人或与他人合作创作之成果，或对所投作品拥有合法的著作权，无第三人对其作品提出可成立之权利主张。 2、投稿人保证向我方所投之稿件，尚未在任何媒体上发表。 3、投稿人保证其作品不含有违反宪法、法律及损害社会公共利益之内容。 4、投稿人向我方所投之作品不得同时向第三方投送，即不允许一稿多投。若投稿人有违反该款约定的行为，则我方有权不向投稿人支付报酬。但我方在收到投稿人所投作品10日内未作出采用通知的除外。 5、投稿人授予我方享有作品专有使用权的方式包括但不限于：通过网络向公众传播、复制、摘编、表演、播放、展览、发行、摄制电影、电视、录像制品、录制录音制品、制作数字化制品、改编、翻译、注释、编辑，以及出版、许可其他媒体、网站及单位转载、摘编、播放、录制、翻译、注释、编辑、改编、摄制。 6、投稿人委托我方声明，未经我方许可，任何网站、媒体、组织不得转载、摘编其作品。

基于高阶有限元法的均匀流中物体二阶波辐射问

作者:

关键词:

摘要：

波-流-体相互作用是船舶与海洋工程中最具挑战性的课题之一．流的存在会显著改变波浪参数，如波高、频率等．船舶和钻井平台等海洋结构物在波-流联合作用下,其受力和运动发生显著变化,甚至会导致结构物的破坏从而影响其安全性．因此，对波-流-体相互作用研究在海洋工程领域备受关注．

文献[1]中基于时域线性理论,研究了位于均匀流中半圆柱的二维辐射和绕射问题,分析了均匀流对漂浮物体处的波浪爬高及作用在物体上的水动力的影响．随后,将此方法拓展到了三维垂直圆柱绕射问题[2]．文献[3]中也作了类似的研究，研究表明:当入射波传播方向和水流方向相同时,上游处的波浪爬高会显著增大,同时水流对作用在结构物上的水动力载荷有明显的影响,特别是对二阶平均力的影响尤为显著．文献[4-5]中采用时域二阶理论,计算了三维垂直圆柱处的波浪爬高及其所受水动力．国内学者也作了类似的工作[6-7]．目前完全非线性理论的波-流和结构物相互作用方法已经建立并得到应用[8-13]．相比较线性或二阶理论而言,基于非线性理论的数值计算需要更多的计算资源．

以上波-流-体相互作用研究主要是关于绕射问题,而对于波-流-体辐射问题的研究则很少,主要工作基于线性理论，如文献[1]．文中基于时域二阶理论,采用高阶有限元法研究均匀流中半圆柱的辐射问题,将文献[1]中的位于均匀流中的线性辐射问题扩展到二阶,计算自由表面波和半圆柱所受到的波浪力,分析水流对一、二阶波及作用在半圆柱上的一、二阶水动力的影响,同时讨论波和水动力在不同流速下的非线性特性．

1 控制方程

定义坐标系oxy(图1),x轴和静水面重合,y轴垂直向上,在无波动时水深为h．Sf为瞬时自由表面,Sb为半圆柱表面(物面),其单位外方向矢量为N=(Nx,Ny)．水底表面Sbot位置为y=-h,Sc为流体域?的人工截断边界．假设流体为均匀、不可压缩和无粘性的理想流体,且流体的运动为无旋．半圆柱运动时流场速度势φ满足Laplace方程:

若流场内存在大小为U0并沿x轴正方向的稳态流,则流场总的速度势为:

图1坐标系

Φ在流场内也满足拉普拉斯方程以及边界条件:

以及在Sc上满足辐射条件．方程中t、g和η分别为时间、重力加速度和波高,下标表示对其的导数．将方程式(2)代入式(3)～(6)并将φ和η进行摄动，展开到二阶:

式中：ε为小量(一般可取为波陡)；φ(0)为由水流作用引起的扰动势,而由其造成的波高忽略不计．φ(k)、η(k)(k=1,2)分别表示第k阶速度势和波高．同时,将水平位移X和角位移Θ也展开到二阶,在求解给定简谐运动的辐射问题时X和Θ只取到一阶,即X=εX(1),Θ=εΘ(1)．将方程式(7～8)代入方程式(1～6)可得各阶(k=0, 1, 2)定解问题:

式中，?(0)是由自由表面平均位置人工边界Sc、物面平均位置和水底部Sbot构成的流域．n=(nx,ny)为平均物面上的外法线方向．在方程中,表达式分别为:

φ(1)|2-

g(0)=-U0nx

在0阶时,η(0)=0以及此外,由于在时域内求解,需要提供初始条件．作用在物体上水动力为：

式中，p为压力,由Bernoulli方程有：

将式(2)代入式(15)再代入式(14)中并进行摄动，展开得:

F=F(0)+F(1)+F(2)+

式中，F(0)为由流引起的力,F(1)、F(2)和分别为一阶力、二阶力和二阶平均力[5]．

2 数值方法与计算分析

混合边值问题式(9～13)通过有限元法求解:

式中：K为有限元系数矩阵；F(k)为右端向量,其元素表达式分别为:

式中，Ni为单元插值函数，文中采用12节点等参单元．由于该类单元形函数为三次插值,对其二次求导后仍保持连续,故可以方便而精确地计算出二阶导数，如方程式(10、11、12)中的φxx、φyy和φxy．

由于在时域内计算,每一时刻自由面上的波高和速度势必须知道,文中采用四阶Adams-Bashforth格式计算(t+Δt)时步自由表面波高和速度势．同时,采用人工阻尼区(图1)来吸收反射波[14],保证外传波条件．

对均匀流中半径为a的半圆柱绕射问题进行计算, 水深h=10a．半圆柱做水平或垂直运动：

式中：A为运动幅值；ω为运动频率；t为时间．

先考虑付氏数Fn=0.064和0.128时的情况,付氏数定义为计算时取A/a=0.1．文献[1]中基于线性理论,用常数面元法研究半圆柱辐射问题，计算了离半圆柱无限远处的线性波(一阶波)．文中计算距离半圆柱中心5a处即x=±5a处线性波随无因次波数Ka的变化关系，波数K和ω的关系为ω2=KgtanhKh．图2为两种不同付氏数，文献[1]和文中计算的一阶波幅值|η(1)/A|结果的对比．无因次波数Ka在0.4～2.0，其符合程度总体较好．

文章来源：《辐射研究与辐射工艺学报》网址: http://www.fsyjyfsgyxb.cn/qikandaodu/2021/0331/683.html