«上一篇
文章快速检索     高级检索
下一篇»
  哈尔滨工程大学学报  2021, Vol. 42 Issue (7): 967-974  DOI: 10.11990/jheu.202005057
0

引用本文  

龚逸纲, 卢东强, 浦俊. 三层流体中入射波作用下直立圆柱支撑圆板的水弹性响应[J]. 哈尔滨工程大学学报, 2021, 42(7): 967-974. DOI: 10.11990/jheu.202005057.
GONG Yigang, LU Dongqiang, PU Jun. Wave-induced hydroelastic response of a circular plate supported by a vertical bottom-mounted cylinder in a three-layer fluid[J]. Journal of Harbin Engineering University, 2021, 42(7): 967-974. DOI: 10.11990/jheu.202005057.

基金项目

国家自然科学基金面上项目(11872239)

通信作者

卢东强, E-mail: dqlu@shu.edu.cn

作者简介

龚逸纲, 男, 博士研究生;
卢东强, 男, 研究员, 博士生导师

文章历史

收稿日期:2020-05-30
网络出版日期:2021-04-26
三层流体中入射波作用下直立圆柱支撑圆板的水弹性响应
龚逸纲 1, 卢东强 1,2,3, 浦俊 1     
1. 上海大学 力学与工程科学学院, 上海 200072;
2. 上海市应用数学和力学研究所, 上海 200072;
3. 上海市能源工程力学重点实验室, 上海 200072
摘要:将海洋超大型浮式结构简化成弹性薄板, 固结于贯穿海洋全深度的直立圆柱以限制其刚体位移, 将海洋近似为有限深三层均匀流体。本文在理想、不可压流体、无旋流场和小振幅波的假设下, 解析研究了入射波浪与超大型浮式结构的水弹性相互作用, 求解波浪的散射和结构受力。将整个流场划分为开阔水域与板覆盖水域, 分别在2个区域对速度势函数进行展开, 利用开阔水域垂向本征函数和适合三层流体的新型内积定义式, 对匹配方程以及直立圆柱侧面边界条件做内积, 以简化方程并最终求解。分析了不同流体、弹性板参数与不同位置波入射时结构的水弹性响应。研究结果表明:高阶衰减模态对结果的贡献较小, 在力、弯矩上基本可以忽略不计。表面波模态与界面波模态入射相比, 引起的波幅在开阔水域各个界面上较大, 但相应模态的波数较小。在界面波模态入射时, 有能量从界面波模态转移到了表面波模态。
关键词超大型浮式结构物    三层流体    本征函数展开法    势流理论    水弹性    色散关系    内积关系    贝塞尔函数    
Wave-induced hydroelastic response of a circular plate supported by a vertical bottom-mounted cylinder in a three-layer fluid
GONG Yigang 1, LU Dongqiang 1,2,3, PU Jun 1     
1. School of Mechanics and Engineering Science, Shanghai University, Shanghai 200072, China;
2. Shanghai Institute of Applied Mathematics and Mechanics, Shanghai 200072, China;
3. Shanghai Key Laboratory of Mechanics in Energy Engineering, Shanghai 200072, China
Abstract: A very large floating structure is modeled as an elastic thin plate, which is fixed with a vertical bottom-mounted cylinder to restrain the displacement of the floating rigid body. The ocean is simplified as a three-layer fluid, and each layer of fluid is homogeneous. Under the assumptions of ideal, incompressible fluid, irrotational flow motion, and small amplitude wave, the hydroelastic interaction between incident waves and the very large floating structure is studied analytically to demonstrate the wave scattering and the forces on the structure. The whole flow field is divided into open-water and plate-covered regions, for which the velocity potential functions are expanded. The vertical eigenfunctions in the open-water region and a newly defined inner product for the three-layer fluid are used to perform an inner product operation on the matching equations and boundary conditions so as to simplify the equation and solve it finally. The hydroelastic responses with different parameters of the fluid, elastic plate, and different modes of incident waves are analyzed. The results show that the higher-order evanescent mode has a small contribution to the result, which can be ignored in the calculation for the force and bending moment. Compared with the incidence of the interface wave mode, the surface wave mode causes a larger amplitude at each interface in the open-water region and a smaller wavenumber of the corresponding mode. When the interface wave mode is incident, there is an energy transfer from the interface wave mode to the surface one.
Keywords: VLFS    three-layer fluid    method of matched eigenfunction expansions    potential flow    hydroelasticity    dispersion relation    inner product    Bessel function    

海洋超大型浮式结构物(very large floating structures, VLFS)作为一个远离陆地、稳固可靠的海上基地, 对存储、开采海洋能源以及军事防护都有重大意义, 引起了学者的广泛关注。超大型浮式结构的厚度远小于水平方向的尺寸, 使它的弯曲刚度变得较小, 在入射波浪和其他外加载荷作用下其水弹性响应十分重要。在理论研究中, 通常运用Kirchhoff-Love薄板理论进行分析[1-2]。对于二维问题, 则将其视作为Timoshenko和Woinwsky-Krieger弹性梁模型。Watanabe[3]、赵存宝等[4]提出了采用中厚度板理论, 可以避免采用Kirchhoff-Love板理论时应力合力不准确的问题。本文将VLFS简化为弹性薄板模型, 固结一根贯穿海洋深度的直立圆柱来固定它的刚体位移, 同时起到一定的支撑作用。在数学上, 外载荷与弹性薄板互相作用这一问题, 在不可压势流理论的框架下可以归结成为一个Laplace方程初边值问题, 常见的求解途径分为频域解法与时域解法。

频域分析适合于线性稳态问题, 主要有模态分析法与本征函数展开法2类。Wu等[5]在计算有限长弹性板在有限水深下的水动力响应时, 将板的刚体运动和弯曲变形用梁的自然模态展开, 然后解出流体相应的绕射势和辐射势, 再运用Bernoulli方程计算出流体压强, 最后根据梁方程计算出每个模态的振幅。Kashiwagi[6]将板的垂向位移用梁的自然模态展开, 将压力用B样条函数展开, 通过Galerkin方法得到压力表达式系数, 计算高效且准确性较高, 成为工程实际中广泛应用的一种方法。模态分析法的难点在于模态函数的选取, 不同模态函数的选取对应的计算时间差异很大。本征函数展开法则是在只考虑结构物水弹性响应而不考虑刚体位移的前提下, 将速度势函数通过分离变量法展开, 其关键就是求解展开的速度势函数中未知系数。Sahoo等[7]构造了内积关系式, 使得本征函数在新内积意义下具有正交性, 简化了计算。Xu等[8]验证了直接使用开阔水域的垂向本征函数可以得到更高的计算效率。

对于实际的密度分层海洋, 两层流体是常见的简化模型。Linton等[9]基于线性水波理论考虑了二层流体中波浪和水平放置圆柱的相互作用并运用多极展开法进行求解。Lin等[10]研究了有限深两层流体中固结在浮冰中的直立圆柱的绕射问题, 讨论了表面波模态和界面波模态入射时, 不同波模态对圆柱体的水平作用力和弯矩的影响。Lin等[11]认为在二层有限深流体中内波入射下, 弹性板的存在会使2种波模态之间的能量交换增加。更精细地看, 海洋是一种三层式的布局。表层密度变化十分缓慢, 随着深度不断增加, 在一个深度密度梯度突然增大, 出现一个密度跃层, 最后密度变化又趋于平缓。Rusås等[12]提出了一种基于积分方程的完全非线性数值格式, 研究了在三层流体中水平传播的界面对称孤立波。Meng等[13]探究了三层流体中漂浮的弹性薄板受到波浪引起的水弹性响应。浦俊[14]研究了三层流体中斜入射波作用下半无限板的水弹性响应, 得到了不同波模态入射时的临界入射角。对于弹性板—直立圆柱模型, Meng等[15]导出了多层流体漂浮多模块板的水弹性波色散关系。

本文将海洋简化为三层密度不同的均匀流体, 下层流体密度大于上层。同时在界面处密度发生突变, 并有界面波入射, 其中, 第2层流体是为了考虑密度跃层的影响而假定的一层流体。

1 数学模型建立

本文主要讨论三层流体中入射重力波与直立圆柱支撑弹性圆板复合结构的相互作用。如图 1所示, 采用三维柱坐标系, 取竖直向上的为z轴, 并将圆盘的中心设为坐标的原点。弹性圆板与直立圆柱的连接方式为固支, 直立圆柱被视为刚体。弹性圆板的厚度为d, 半径为R, 直立圆柱的半径为a。将z=0, -H1, -H2, -H3分别定义为平静水面、第1层流体与第2层流体的交界(上界面)、第2层流体与第3层流体的交界(下界面)以及水底。不同层流体之间有明显且稳定的界面。第1层流体(-H1z < 0)、第2层流体(-H2z<-H1)与第3层流体(-H2z<-H3)的密度分别为ρ1ρ2ρ3, 且ρ3>ρ2>ρ1, 每层流体的厚度为hm(m=1, 2, 3)。将整个水域划分为开阔水域(Rr<∞)和板覆盖区(ar < R)。在z=0处受到表面入射波的作用, 并在z=-H1以及z=-H2处分别受到上界面、下界面入射波的影响。

Download:
图 1 三层流体中重力波与圆形弹性板及直立圆柱相互作用 Fig. 1 Schematic diagram of gravity wave with circular elastic plate and vertical cylinder in three-layer fluid

本文基于势流理论, 假设流体理想且不可压, 运动无旋。令Φm(r, θ, z, t)(m=1, 2, 3)为第1、2、3层流体速度势函数, ξm(r, θ, t)(m=1, 2, 3)为自由水面、上下界面的铅垂位移。将时间自变量t从上述变量中分离出来, 得到:

$ \varPhi_{m}(r, \theta, z, t) =\operatorname{Re}\left[\phi_{m}(r, \theta, z) \mathrm{e}^{\mathrm{i} \omega t}\right] $ (1)
$ \xi_{m}(r, \theta, t) =\operatorname{Re}\left[\xi_{m}(r, \theta) \mathrm{e}^{\mathrm{i} \omega t}\right] $ (2)

速度势函数ϕm(r, θ, z)在整个流域内满足Laplace方程:

$ \left(\nabla^{2}+\frac{\partial^{2}}{\partial z^{2}}\right) \phi_{m}=0, a \leqslant r<\infty,-H_{3}<z<0 $ (3)

其中:

$ \nabla^{2}=\frac{\partial^{2}}{\partial r^{2}}+\frac{\partial}{r \partial r}+\frac{\partial^{2}}{r^{2} \partial \theta^{2}} $ (4)

一阶近似下的自由水面边界条件为:

$ -\omega^{2} \phi_{1}+g \frac{\partial \phi_{1}}{\partial z}=0, R \leqslant r<\infty, z=0 $ (5)

板覆盖区水面边界条件为:

$ \begin{gathered} -\rho_{1} \omega^{2} \phi_{1}+\left(D \nabla^{4}-M_{\mathrm{e}} \omega^{2}+\rho_{1} g\right)=\frac{\partial \phi_{1}}{\partial z}=0, \\ a \leqslant r<R, z=0 \end{gathered} $ (6)

式中:D为弹性圆板的抗弯刚度;Me=ρed为弹性圆板单位面积上的质量。

水底的非渗透条件为:

$ \frac{\partial \phi_{3}}{\partial z}=0, \quad z=-h_{3} $ (7)

在各层流体之间的界面处流体速度、压力满足连续性条件:

$ \frac{\partial \phi_{n}}{\partial z}=\frac{\partial \phi_{n+1}}{\partial z}=-\mathrm{i} \omega \xi_{n+1}, \quad z=-H_{n}, a \leqslant r<\infty $ (8)
$ \begin{gathered} \gamma_{n}=\left(\frac{\omega^{2}}{g} \phi_{n}-\frac{\partial \phi_{n}}{\partial z}\right)=\frac{\omega^{2}}{g} \phi_{n+1}-\frac{\partial \phi_{n+1}}{\partial z}, \\ z=-H_{n}, a \leqslant r<\infty \end{gathered} $ (9)

式中:γn=ρn/ρn+1; n=1, 2。

考虑弹性圆板的外缘是自由端, 即在r=R处, 弹性圆板的力与弯矩为零, 在r=a处, 直立圆柱与弹性圆板固支, 因此其位移与转角为0。

$ \left[\frac{\partial}{\partial r^{2}}+v\left(\frac{1}{r} \frac{\partial}{\partial r}+\frac{1}{r^{2}} \frac{\partial^{2}}{\partial \theta^{2}}\right)\right] \xi_{1}=0, r=R $ (10)
$ {\left[\frac{\partial}{\partial r^{2}} \nabla^{2}+\frac{1-v}{r^{2}}\left(\frac{\partial}{\partial r}-\frac{1}{r^{2}}\right) \frac{\partial^{2}}{\partial \theta^{2}}\right] \xi_{1}=0, r=R} $ (11)
$ \xi_{1}=0, \quad r=a $ (12)
$ \frac{\partial \xi_{1}}{\partial r}=0, \quad r=a $ (13)

在板覆盖区与开阔水域的交界处, 即r=R处, 速度势满足压力与速度连续的匹配条件为:

$ \left.\phi_{m}\right|_{r=R^{-}}=\left.\phi_{m}\right|_{r+R^{+}}, \quad m=1,2,3 $ (14)
$ \left.\frac{\partial \phi_{m}}{\partial r}\right|_{r=R^{-}}=\left.\frac{\partial \phi_{m}}{\partial r}\right|_{r=R^{+}}, \quad m=1,2,3 $ (15)

在直立圆柱的侧面, 流体与固体之间没有间隙, 因此不能有沿圆柱法向的速度分量:

$ \left.\frac{\partial \phi_{m}}{\partial r}\right|_{r=a^{+}}=0, \quad m=1,2,3 $ (16)
2 数学模型求解过程

色散关系满足方程[13]

$ D_{6} \omega^{6}-D_{4} \omega^{4}+D_{2} \omega^{2}-D_{0}=0 $ (17)

其中具体参数定义以及垂向本征函数表达式见文献[13]。

在三层流体的开阔水域, 对于给定频率的入射波会以3个不同的波数进行传播: 表面波模态(k01), 上界面波模态(k02)以及下界面波模态(k03), 同时会有无穷多个虚数波数的耗散波ikn(n=1, 2, 3, …)产生。而在板覆盖区, 除了与开阔水域中相同数量的波数$\tilde{k}_{0_{1}}、\tilde{k}_{0_{2}}、\tilde{k}_{0_{3}}、\mathrm{i} \tilde{k}_{n}$(n=1, 2, 3, …), 还会额外产生2个复根$\mathrm{i} \tilde{k}_{I}、\mathrm{i} \tilde{k}_{I I}$。不同区域对应的波数可以由不同区域的色散关系分别求解出来。

为简洁, 将速度势函数ϕm写成分段函数形式:

$ \phi=\left\{\begin{array}{ll} \phi_{1}, & -H_{1}<z<0 \\ \phi_{2}, & -H_{2} \leqslant z \leqslant-H_{1} \\ \phi_{3}, & -H_{3}<z<-H_{2} \end{array}\right. $ (18)

在整个流场中, 速度势函数可以分解为开阔水域速度势ϕO(r, θ, z)以及板覆盖区速度势$\widetilde{\phi}$(r, θ, z)。在开阔水域, 速度势可分解为入射势ϕI(r, θ, z)及绕射势ϕD(r, θ, z), ϕO=ϕI+ϕD。根据本征函数展开法, 各势函数可以展开为:

$ \phi^{\mathrm{I}}(r, \theta, z)=\mathrm{I}_{0_{m}} \sum\limits_{i} \mathrm{i}^{n} \mathrm{~J}_{n}\left(k_{0_{m}} r\right) Z_{0_{m}} $ (19)
$ \begin{aligned} \phi^{\mathrm{D}}(r, \theta, z)=& \sum\limits_{n}\left[\sum\limits_{m} A_{n, 0_{m}} \mathrm{H}_{n}\left(k_{0_{m}} r\right) Z_{0_{m}}+\right.\\ &\left.\sum\limits_{i} A_{n, i} \mathrm{~K}_{n}\left(k_{i} r\right) Z_{i}\right] \end{aligned} $ (20)
$ \begin{gathered} \widetilde{\phi}(r, \theta, z)= \\ \sum\limits_{n}\left\{\sum\limits_{m}\left[B_{n, 0_{m}} \mathrm{~J}_{n}\left(\tilde{k}_{0_{m}} r\right)+C_{n, 0_{m}} \mathrm{Y}_{n}\left(\tilde{k}_{0_{m}} r\right)\right] \widetilde{Z}_{0_{m}}+\right.\\ \left.\sum\limits_{j}\left[B_{n, j} \mathrm{I}_{n}\left(\tilde{k}_{j} r\right)+C_{n, j} \mathrm{~K}_{n}\left(\tilde{k}_{j} r\right)\right] \widetilde{Z}_{j}\right\} \end{gathered} $ (21)

其中:

$ \mathrm{I}_{0_{m}}=-\mathrm{i} \omega {\xi}_{m}\left[\frac{\partial V\left(k_{0_{m}}, z\right)}{\partial z}\right]_{z=-H_{m-1}}^{-1} $ (22)
$ \left\{\begin{array}{l} Z_{0_{m}}(z)=V\left(k_{0_{m}}, z\right) \\ Z_{i}(z)=V\left(\mathrm{i} k_{i}, z\right) \\ \widetilde{Z}_{0_{m}}(z)=V\left(\tilde{k}_{0_{m}}, z\right) \\ \widetilde{Z}_{j}(z)=V\left(\mathrm{i} \tilde{k}_{j}, z\right) \end{array}\right. $ (23)

式中: 求和号$\sum\limits_{i} \sum\limits_{j} \text { 与 } \sum\limits_{m}$的求和范围分别是$\sum\limits_{i=1}^{\infty}$$\sum\limits_{j=I, I I, 1}^{\infty} 与\sum\limits_{m=1}^{3}$; 特别地, $\sum\limits_{n} \text { 为 } \sum\limits_{n=-\infty}^{\infty} {\rm{e}}^{\mathrm{i}n \theta}$的缩写; Jn(·), Yn(·), Hn(·), In(·) 与Kn(·)分别为n阶第一类Bessel函数、第二类Bessel函数、第一类Hankel函数、第一类修正Bessel函数以及第二类修正Bessel函数; An, 0mAn, iBn, 0mBn, jCn, 0mCn, j为待定的复系数。

为了得到上述复系数, 将展开后的速度势函数代入弹性圆板与直立圆柱在z=0处的边界条件(10)~(13)。

根据垂向本征函数的正交性, 引入垂向本征函数的内积定义为:

$ \begin{aligned} &P_{l n}=<Z_{l}, Z_{n}>=\int_{-H_{3}}^{-H_{2}} Z_{l} Z_{n} \mathrm{~d} z+ \\ &\gamma_{2} \int_{-H_{2}}^{-H_{1}} Z_{l} Z_{n} \mathrm{~d} z+\gamma_{1} \gamma_{2} \int_{-H_{1}}^{0} Z_{l} Z_{n} \mathrm{~d} z \end{aligned} $ (24)

其中下标取值范围为l, n为01, 02, 03, 1, 2, 3, …。可以通过计算验证在ln时, Pln=0。首先将速度势函数代入匹配条件(14)、(15)以及直立圆柱侧面边界条件(16)。再将本征函数Zl(z)(l=01, 02, 03, 1, 2, 3, …)分别对上述匹配条件与直立圆柱侧面边界条件做内积:

$ \begin{gathered} \sum\limits_{n}\left[I_{l} i^{n} \mathrm{~J}_{n}\left(k_{l} R\right) P_{l l}+A_{n, l} \mathrm{H}_{n}\left(k_{l} R\right) P_{l l}\right]= \\ \sum\limits_{n}\left\{\sum\limits_{m}\left[B_{n, 0_{m}} \mathrm{~J}_{n}\left(\tilde{k}_{0_{m}} R\right)+C_{n, 0_{m}} \mathrm{Y}_{n}\left(\tilde{k}_{0_{m}} R\right)\right] Q_{0_{m} l}+\right. \\ \left.\sum\limits_{j}\left[B_{n, j} \mathrm{I}_{n}\left(\tilde{k}_{j} R\right)+C_{n, j} \mathrm{~K}_{n}\left(\tilde{k}_{j} R\right)\right] Q_{j l}\right\}, \\ l \text { 为 } 0_{1}, 0_{2}, 0_{3} \end{gathered} $ (25)
$ \begin{gathered} \sum\limits_{n} A_{n, i} \mathrm{~K}_{n}\left(k_{i} R\right) P_{i i}= \\ \sum\limits_{n}\left\{\sum\limits_{m}\left[B_{n, 0_{m}} \mathrm{~J}_{n}\left(\tilde{k}_{0_{m}} R\right)+C_{n, 0_{m}} \mathrm{Y}_{n}\left(\tilde{k}_{0_{m}} R\right)\right] Q_{0_{m} i}+\right. \\ \left.\sum\limits_{j}\left[B_{n, j} \mathrm{I}_{n}\left(\tilde{k}_{j} R\right)+C_{n, j} \mathrm{~K}_{n}\left(\tilde{k}_{j} R\right)\right] Q_{j i}\right\}, \\ l=i=1,2, \cdots \end{gathered} $ (26)
$ \begin{gathered} \sum\limits_{n}\left[I_{l} i^{n} \mathrm{~J}_{n}^{\prime}\left(k_{l} R\right) P_{l l}+A_{n, l} \mathrm{H}_{n}^{\prime}\left(k_{l} R\right) P_{l l}\right]=\\ \sum\limits_{n}\left\{\sum\limits_{m}\left[B_{n, 0_{m}} \mathrm{~J}_{n}^{\prime}\left(\tilde{k}_{0_{m}} R\right)+C_{n, 0_{m}} \mathrm{Y}_{n}^{\prime}\left(\tilde{k}_{0_{m}} R\right)\right] Q_{0_{m} l}+\right. \\ \left.\sum\limits_{j}\left[B_{n, j} \mathrm{I}_{n}^{\prime}\left(\tilde{k}_{j} R\right)+C_{n, j} K_{n}^{\prime}\left(\tilde{k}_{j} R\right)\right] Q_{j l}\right\}, \\ l \text { 为 } 0_{1}, 0_{2}, 0_{3} \end{gathered} $ (27)
$ \begin{gathered} \sum\limits_{n} A_{n, i} \mathrm{~K}_{n}^{\prime}\left(k_{i} R\right) P_{i i}= \\ \sum\limits_{n}\left\{\sum\limits_{m}\left[B_{n, 0_{m}} \mathrm{~J}_{n}^{\prime}\left(\tilde{k}_{0_{m}} R\right)+C_{n, 0_{m}} \mathrm{Y}_{n}^{\prime}\left(\tilde{k}_{0_{m}} R\right)\right] Q_{0_{m} i}+\right. \\ \left.\sum\limits_{j}\left[B_{n, j}\mathrm{~I}_{n}^{\prime}\left(\tilde{k}_{j} R\right)+C_{n, j} K_{n}^{\prime}\left(\tilde{k}_{j} R\right)\right] Q_{j i}\right\} ,\\ l=i=1,2, \cdots \end{gathered} $ (28)
$ \begin{gathered} \sum\limits_{n}\left\{\sum\limits_{m}\left[B_{n, 0_{m}} \mathrm{~J}_{n}^{\prime}\left(\tilde{k}_{0_{m}} a\right)+C_{n, 0_{m}} \mathrm{Y}_{n}^{\prime}\left(\tilde{k}_{0_{m}} a\right)\right] Q_{0_{m} l}+\right. \\ \left.\sum\limits_{j}\left[B_{n, j} \mathrm{I}_{n}^{\prime}\left(\tilde{k}_{j} a\right)+C_{n, j} \mathrm{~K}_{n}^{\prime}\left(\tilde{k}_{j} a\right)\right] Q_{j l}\right\}=0 ,\\ l \text { 为 } 0_{1}, 0_{2}, 0_{3} \end{gathered} $ (29)
$ \begin{gathered} \sum\limits_{n}\left\{\sum\limits_{m}\left[B_{n, 0_{m}} \mathrm{~J}_{n}^{\prime}\left(\tilde{k}_{0_{m}} a\right)+C_{n, 0_{m}} \mathrm{Y}_{n}^{\prime}\left(\tilde{k}_{0_{m}} a\right)\right] Q_{0_{m} i}+\right. \\ \left.\sum\limits_{j}\left[B_{n, j} \mathrm{I}_{n}^{\prime}\left(\tilde{k}_{j} a\right)+C_{n, j} \mathrm{~K}_{n}^{\prime}\left(\tilde{k}_{j} a\right)\right]_{j i}\right\}=0 ,\\ l=i=1,2, \cdots \end{gathered} $ (30)

对无穷累加项i, j做截断, 取i=1, 2, …, M; j=I, II, 1, 2, …, M以及n=-N, …, -1, 0, 1, …, NMN为正整数。根据上述边界条件及匹配条件得到一个(2N+1)(3M+13)维的方程组, 相应的未知系数为(2N+1)(3M+13)个, 方程封闭。通过求解该方程组, 求解出待定系数, 可以解出速度势函数表达式, 从而进一步得出其他物理量的值。

在得到了整个流场内势函数表达式之后, 可以进一步分析弹性圆板与直立圆柱这2个结构物上的受力情况。根据本文所采用的线性势流理论, 直立圆柱侧面(r=a)处的任意一点的动压可以根据线性化后的Bernoulli方程给出。

3 数值结果与讨论

为不失一般性, 取流场深度H3, 重力加速度g以及上层流体的密度ρ1做无量纲化。后文中出现的物理量均为无量纲化后的物理量。

采用上述方法进行数值运算。弹性圆板的参数数值为E=1.2×106, ν=0.3, ρe=0.9, d=0.01。流场中流体参数为γ1=γ2=0.9, h1=h2=0.2, h3=0.6。首先对求解的正确性进行检验并判断收敛性。正确性的判据选取为检验开阔水域与板覆盖区在交界处的波流能守恒。参考文献[16], 波流能εf可以表示为:

$ {\varepsilon _f} = \overline {\int_{ - {H_3}}^{{\zeta _1}} p u{\rm{d}}z} \cong - \rho \int_{ - {H_3}}^0 {\overline {\frac{{\partial \mathit{\Phi }}}{{\partial t}}\frac{{\partial \mathit{\Phi }}}{{\partial r}}} {\rm{d}}z} $ (31)

M=5, N=15, R=4, a=1, ω=0.5, d=0.01的算例中, 相对误差为:

$ \frac{|E-\widetilde{E}|}{\min [E, \widetilde{E}]} \approx 0.000\ 85 \% $ (32)

对于衰减模态项数的截断参数M, 当M≥5时, 衰减模态对结果的影响就可忽略不计, 本文后续计算均在M=7下进行。对于Bessel函数阶数截断参数N, 数值的收敛性随ω的变化而变化。在ω较小时, 如ω=0.5, N=10便已经收敛, 但在ω=2时, N=20才能达到收敛。本文后续计算全部取N=20。

图 2展示了在不同入射频率ω下上界面上竖直位移的幅值。在ω=0.1时, 波受到直立圆柱的干扰较小, 但当ω增大到1时, 直立圆柱的存在对后面的流场产生了很大的影响, 符合低频长波绕射能力强, 对下游的影响更大。同时, 随着频率的增加, 直立圆柱一周的竖向位移不断减小。

Download:
图 2 表面波模态入射时, 不同ω下上界面处的幅值 Fig. 2 The amplitude at the upper-layer interfaces of different ω due to the incident surface wave mode

图 3展示了在表面波模态入射时, 不同板厚d, 密度比γ和每层流体的厚度h下, 表面、上下界面θ=π处的竖直位移幅值, 其中选定ω=0.5, ζ1=0.01, ζ2=ζ3=0。当板厚d不断变大时, 弹性板的抗弯刚度不断增加, 因此板的挠度更小同时板覆盖区的波幅也相对较小, 在x=-1处满足弹性圆板与直立圆柱的固结条件, 同时因为抗弯刚度的增加, 使得开阔水域的波幅增大。当密度比γ变化时, 板的挠度并没有明显的变化, 但当密度比不断减小, 明显在自由水面上激发了新的模态, 波形变得更加粗糙, 但在上下界面模态波数则不断变小, 同时波幅有所增加。最后, 当第1层与第2层流体的厚度不断增加时, 表面波波形基本不变, 上下界面处的波幅则变得更小, 激发的模态波长也更长, 表明在表面波模态入射时, 上层流体越厚, 传到下层的能量也就越少。

Download:
图 3 表面波模态入射时, 不同d, γ, h下不同界面θ=π处的幅值 Fig. 3 The amplitude at the upper-layer interfaces of different d, γ, h due to the incident surface wave mode

图 4讨论了表面波模态入射时, 不同半径直立圆柱、弹性圆板的径向压力s(a, θ)。图 4(a)给出的是弹性圆板与直立圆柱半径基本相同的模型, 可以看出, 当直立圆柱的半径不断增大, 受到的径向压力就越小, 同时随着直立圆柱尺寸增加, 圆柱后方压力分布变得更加复杂。在a=1时可以与图 4(b)相较得出, 板的存在吸收了部分的能量。

Download:
图 4 表面波模态入射时, 不同圆柱半径a与弹性圆板半径R下圆柱的径向压力s(a, θ) Fig. 4 The radial pressure s(a, θ) of the cylinder under different radii a and R of the elastic circular plate due to incident surface wave mode

图 5为表面波模态入射时, 不同入射频率ω下圆柱收到的切向压力q(z)。因为在界面上势函数并不连续, 因此水平方向的压力会有突变。随着入射频率不断增加, 水平方向压力减小, 但当ω取很小的值, 如ω=0.1, 圆柱受到的压力更均匀, 因此水平方向的压力就变得很小。

Download:
图 5 表面波模态入射时, 不同ω下圆柱水平方向受到的压力q(z) Fig. 5 The horizontal pressure q(z) of the cylinder under different ω due to incident surface wave mode

图 6对比了同样波幅下表面、上下界面波模态分别入射时, 各个界面θ=π处的竖直位移幅值。只有在表面波模态入射时, 在表面上会引起较大的振幅, 而当上下界面波模态入射时, 在开阔水域会激起波数更大的新模态, 但波幅更小。在板覆盖水域则在上下界面处引起较大的波幅, 但板的挠度几乎为零。相比表面入射波, 上下界面入射波所产生的水平方向压力与径向压力都要小得多, 上表面波模态入射时, 在第2层流体内, 水平方向压力明显减小, 与下表面波模态入射时相反。

Download:
图 6 表面、上下界面波模态分别入射时, 各个界面θ=π处的幅值以及圆柱所受到的径向压力与水平方向的压力 Fig. 6 The amplitudes at each interface and the radial and horizontal pressure on the cylinder due to the incident wave modes of the surface and the upper and lower interfaces
4 结论

1) 衰减模态对于结果的贡献较小, 在力、弯矩上基本可以忽略不计, 在数值计算时并不需要考虑高阶项的影响。而Bessel函数的阶数对数值收敛性有影响, 求解高频短波对阶数要求更高。

2) 在低频长波入射时, 行进波的绕射能力强, 受到弹性圆板直立圆柱结构的影响很小; 当频率增加, 模型对行进波的影响就越大。

3) 表面波模态与界面波模态入射相比, 引起的波幅在开阔水域各个界面上较大, 但相应模态的波数越小。

4) 在界面波模态入射时, 有能量从界面波模态转移到了表面波模态。

参考文献
[1]
KASHIWAGI M. A time-domain mode-expansion method for calculating transient elastic responses of a pontoon-type VLFS[J]. Journal of marine science and technology, 2000, 5(2): 89-100. DOI:10.1007/PL00010631 (0)
[2]
KOROBKIN A A, KHABAKHPASHEVA T I. Regular wave impact onto an elastic plate[J]. Journal of engineering mathematics, 2006, 55(1/2/3/4): 127-150. DOI:10.1007/s10665-005-0191-8 (0)
[3]
WATANABE E, UTSUNOMIYA T, WANG C M, et al. Benchmark hydroelastic responses of a circular VLFS under wave action[J]. Engineering structures, 2006, 28(3): 423-430. DOI:10.1016/j.engstruct.2005.08.014 (0)
[4]
赵存宝, 翟林娜. 组合式超大型浮体结构在外载荷作用下的水弹性响应研究[C]//第八届全国流体力学学术会议论文摘要集. 兰州: 中国力学学会, 2014: 23.
ZHAO Cunbao, ZHAI Linna. Research on hydroelastic response of combined very large floating structure under external load[C]//Abstracts of the 8th National Conference on Fluid Mechanics. Lanzhou: Chinese Society of Mechanics, 2014: 23. (0)
[5]
WU Chong, WATANABE E, UTSUNOMIYA T. An eigenfunction expansion-matching method for analyzing the wave-induced responses of an elastic floating plate[J]. Applied ocean research, 1995, 17(5): 301-310. DOI:10.1016/0141-1187(95)00023-2 (0)
[6]
KASHIWAGI M. A B-spline Galerkin scheme for calculating the hydroelastic response of a very large floating structure in waves[J]. Journal of marine science and technology, 1998, 3(1): 37-49. DOI:10.1007/BF01239805 (0)
[7]
SAHOO T, YIP T L, CHWANG A T. Scattering of surface waves by a semi-infinite floating elastic plate[J]. Physics of fluids, 2001, 13(11): 3215-3222. DOI:10.1063/1.1408294 (0)
[8]
XU Feng, LU Dongqiang. An optimization of eigenfunction expansion method for the interaction of water waves with an elastic plate[J]. Journal of hydrodynamics, 2009, 21(4): 526-530. DOI:10.1016/S1001-6058(08)60180-8 (0)
[9]
LINTON C M, MCIVER M. The interaction of waves with horizontal cylinders in two-layer fluids[J]. Journal of fluid mechanics, 1995, 304: 213-229. DOI:10.1017/S002211209500440X (0)
[10]
LIN Q, LU D Q. Water wave diffraction by a bottom-mounted circular cylinder clamped to an elastic plate floating on a two-layer fluid[J]. European journal of mechanics-B/fluids, 2014, 44: 10-21. DOI:10.1016/j.euromechflu.2013.11.004 (0)
[11]
LIN Q, DU X G, KUANG J, et al. Forces and moments exerted by incident internal waves on a plate-cylinder structure in a two-layer fluid[J]. Meccanica, 2019, 54(10): 1545-1560. DOI:10.1007/s11012-019-01032-0 (0)
[12]
RUSÅS P O, GRUE J. Solitary waves and conjugate flows in a three-layer fluid[J]. European journal of mechanics-B/fluids, 2002, 21(2): 185-206. DOI:10.1016/S0997-7546(01)01163-3 (0)
[13]
MENG Qingrui, LU Dongqiang. Hydroelastic interaction between water waves and thin elastic plate floating on three-layer fluid[J]. Applied mathematics and mechanics, 2017, 38(4): 567-584. DOI:10.1007/s10483-017-2185-6 (0)
[14]
浦俊, 卢东强. 三层流体中斜入射波作用下半无限板的水弹性响应[J]. 力学学报, 2019, 51(6): 1614-1629.
PU Jun, LU Dongqiang. Hydroelastic response of a semi-infinite plate due to obliquely incident waves in a three-layer fluid[J]. Chinese journal of theoretical and applied mechanics, 2019, 51(6): 1614-1629. (0)
[15]
MENG Q R, LU D Q. Scattering of flexural-gravity waves by a group of elastic plates floating on a stratified fluid[J]. European journal of mechanics-B/fluids, 2018, 67: 329-340. DOI:10.1016/j.euromechflu.2017.09.010 (0)
[16]
李玉成, 滕斌. 波浪对海上建筑物的作用[M]. 2版. 北京: 海洋出版社, 2002.
LI Yucheng, TENG Bin. Wave action on maritime structures[M]. 2nd ed. Beijing: Ocean Press, 2002. (0)