规则波中双浮体漂移力的数值模拟与分析

规则波中双浮体漂移力的数值模拟与分析一、引言1.1研究背景与意义随着全球对能源需求的不断增长，海洋资源的开发与利用愈发受到重视。在海洋工程领域，众多浮式结构物被广泛应用，如LNG-FPSO（液化天然气浮式生产储卸装置）系统，它集成了天然气液化、储存和卸载等功能，是海上天然气开采和运输的关键设施。在实际海洋环境中，LNG-FPSO系统通常由多个浮体组成，这些浮体不可避免地会受到海洋波浪的作用。海洋波浪是一种复杂的自然现象，由风、天文因素等多种因素引起，其具有随机性和周期性，会对浮体产生各种作用力，其中漂移力是影响浮体运动和系统稳定性的重要因素之一。漂移力是指由于波浪与浮体相互作用，在一个波浪周期内产生的平均水平力，它会使浮体在水平方向上产生缓慢的漂移运动。对于双浮体系统，如LNG-FPSO与供应船组成的联合系统，两个浮体在规则波中受到的漂移力不仅会影响单个浮体的运动姿态，还会改变它们之间的相对位置和间距。如果漂移力过大或两个浮体受到的漂移力差异较大，可能导致浮体之间的连接装置承受过大的应力，甚至引发连接失效，进而危及整个系统的安全运行。例如，在恶劣海况下，过大的漂移力可能使LNG-FPSO与供应船发生碰撞，造成严重的经济损失和环境污染。此外，漂移力还会影响浮式结构物的定位精度，增加系泊系统的负担，缩短系泊系统的使用寿命，对海洋工程的经济性和可靠性产生不利影响。准确研究两个浮体在规则波中受到的漂移力，对于保障LNG-FPSO系统等海洋工程设施的安全稳定运行具有重要意义。通过深入了解漂移力的特性和规律，可以为浮式结构物的设计提供更准确的力学依据，优化结构设计，提高其抗风浪能力和稳定性。同时，也有助于制定合理的系泊方案和操作策略，降低海洋环境对浮式结构物的影响，确保海上作业的顺利进行，促进海洋资源的可持续开发与利用。1.2研究现状综述在海洋工程领域，关于浮体在波浪中受到的漂移力研究由来已久，众多学者从理论推导、数值计算和实验研究等多个方面展开了深入探索。早期的理论研究主要基于势流理论，将流体视为无粘性、不可压缩的理想流体，通过建立线性化的数学模型来分析波浪与浮体的相互作用。如纽曼（Newman）提出了纽曼积分方法，用于计算浮体在波浪中的一阶波浪力和二阶平均漂移力，为后续研究奠定了重要的理论基础。该方法基于格林函数，将浮体表面的边界条件转化为积分方程进行求解，在一定程度上简化了计算过程。然而，纽曼积分方法在处理复杂形状浮体和非线性问题时存在局限性，计算精度受到一定影响。随着计算机技术的飞速发展，数值计算方法逐渐成为研究漂移力的重要手段。边界元法（BEM）、有限元法（FEM）等数值方法被广泛应用于求解波浪与浮体相互作用的复杂问题。边界元法通过将求解区域的边界离散化，将偏微分方程转化为边界积分方程进行求解，具有降低问题维数、计算效率高的优点，能够精确计算浮体周围的流场分布和波浪力。例如，在计算浮体的二阶平均漂移力时，边界元法可以准确地考虑浮体表面的边界条件，提高计算精度。但边界元法在处理大规模问题时，由于需要存储和计算大量的边界积分系数，内存需求较大，计算效率会受到影响。有限元法则是将求解区域划分为有限个单元，通过对每个单元进行离散化处理，将连续的物理问题转化为代数方程组进行求解，它能够处理复杂的几何形状和边界条件，在模拟波浪与浮体相互作用的多物理场耦合问题时具有优势。不过，有限元法的计算精度依赖于单元的划分质量，若单元划分不合理，可能导致计算结果的误差较大，且计算量通常较大，对计算机硬件要求较高。在实验研究方面，模型试验是获取浮体在波浪中受力和运动数据的重要途径。通过在实验室中模拟真实海洋环境，对不同类型的浮体进行试验，可以直接测量浮体受到的漂移力以及其他相关物理量，为理论和数值计算提供验证依据。例如，在水池试验中，可以精确控制波浪的参数，如波高、波长、周期等，同时测量浮体在波浪作用下的位移、速度和加速度等运动响应，以及所受到的漂移力大小和方向。然而，模型试验存在一定的局限性，由于实验室条件与实际海洋环境存在差异，模型缩比效应可能导致试验结果与实际情况存在偏差，且试验成本较高，难以进行大规模、多参数的研究。对于双浮体系统在规则波中的漂移力研究，一些学者考虑了浮体之间的相互干扰效应。研究发现，双浮体之间的相对位置和间距会对它们受到的漂移力产生显著影响，当双浮体距离较近时，由于波浪在双浮体之间的反射和绕射，会导致浮体周围的流场发生变化，从而改变漂移力的大小和方向。但目前对于双浮体系统漂移力的研究，大多集中在简化的几何模型和特定的波浪条件下，对于复杂海况和实际工程中的多浮体系统，相关研究还相对较少，且不同研究方法之间的对比和验证工作也有待进一步加强。现有研究在考虑海洋流场、风场等环境因素对双浮体漂移力的耦合影响方面还存在不足，难以全面准确地描述实际海洋环境中双浮体的受力情况。1.3研究内容与方法本文围绕两个浮体在规则波中受到的漂移力展开研究，主要内容包括以下几个方面：理论计算：基于势流理论，运用纽曼积分方法，推导双浮体在规则波中漂移力的理论计算公式。考虑双浮体之间的相互干扰效应，引入修正系数，对经典理论公式进行改进，以更准确地描述双浮体在规则波中的受力情况。数值模拟：采用边界元法，建立双浮体在规则波中运动的数值模型。对浮体表面进行网格划分，离散化处理边界积分方程，求解浮体周围的流场分布和波浪力，进而得到漂移力的数值解。通过改变双浮体的相对位置、间距、波浪参数等，分析这些因素对漂移力的影响规律。结果分析：对比理论计算和数值模拟结果，验证理论公式和数值模型的准确性和可靠性。深入分析不同参数条件下双浮体漂移力的变化特性，总结漂移力与各影响因素之间的关系，为海洋工程中双浮体系统的设计和分析提供参考依据。在研究方法上，本文采用理论与数值模拟相结合的方式。理论分析从基本的物理原理出发，建立数学模型，推导漂移力的计算公式，为研究提供理论基础；数值模拟则利用计算机技术，对复杂的物理过程进行数值求解，能够处理更真实的边界条件和多参数变化情况，弥补理论分析的局限性。通过两者相互验证和补充，提高研究结果的准确性和可信度。二、理论基础2.1无界流场中运动物体的附加质量在海洋工程领域，研究浮体在流场中的运动时，附加质量是一个重要的概念。当物体在无界流场中运动时，会带动周围的流体一起运动，这就相当于在物体自身质量的基础上增加了一部分质量，即附加质量。附加质量的存在会对物体的运动特性产生显著影响，准确理解和计算附加质量对于分析浮体在波浪中的运动和受力情况至关重要。2.1.1物体在无界流场中的定解条件假设流体为不可压缩的理想流体，且运动是无旋的，则可引入速度势函数\varphi来描述流场。在无界流场中，物体以速度\vec{U}运动时，速度势\varphi满足拉普拉斯方程：\nabla^{2}\varphi=0在物体表面S上，满足不可穿透条件，即流体的法向速度等于物体表面的法向速度：\frac{\partial\varphi}{\partialn}=\vec{U}\cdot\vec{n}其中，\vec{n}为物体表面的单位法向量。在无穷远处，速度势应满足远方条件，即当r\rightarrow\infty时（r为空间点到物体中心的距离），速度势趋于均匀流的速度势：\varphi\rightarrow\vec{U}\cdot\vec{r}这些定解条件为后续求解附加质量提供了理论基础。通过求解满足上述定解条件的速度势方程，可以得到流场的速度分布，进而计算出附加质量。2.1.2圆球和椭球体的附加质量系数对于单个圆球在无界流场中作匀速直线运动的情况，根据势流理论，其附加质量系数具有明确的解析解。设圆球半径为R，当圆球以速度\vec{U}沿x轴方向运动时，其附加质量m_{a}为：m_{a}=\frac{2}{3}\rho\piR^{3}其中，\rho为流体密度。此时，附加质量系数\lambda_{11}（沿运动方向的附加质量系数）定义为附加质量与圆球自身质量m（m=\frac{4}{3}\rho\piR^{3}）的比值，即：\lambda_{11}=\frac{m_{a}}{m}=\frac{1}{2}这表明圆球在无界流场中运动时，其附加质量为自身质量的一半，该结果在理论和实验中都得到了广泛的验证。对于单个椭球体在无界流场中的运动，其附加质量系数的计算相对复杂，需要考虑椭球体的三个半轴长度a、b、c。当椭球体以速度\vec{U}沿x轴方向运动时，其附加质量系数\lambda_{11}可通过积分表达式计算：\lambda_{11}=\frac{2\piabc}{3}\int_{0}^{\infty}\frac{d\xi}{(a^{2}+\xi)(\Delta)}其中，\Delta=\sqrt{(a^{2}+\xi)(b^{2}+\xi)(c^{2}+\xi)}。通过数值积分等方法可以求解该积分，得到不同半轴长度下椭球体的附加质量系数。当a=b=c时，椭球体退化为圆球，此时计算得到的附加质量系数与上述圆球的结果一致，验证了公式的正确性。对于两个大小相同的椭球体在无界流场中的情况，假设它们的中心连线沿x轴方向，间距为d。由于两个椭球体之间的相互干扰，它们的附加质量系数会发生变化。通过引入镜像法等方法，可以考虑这种相互作用。以沿x轴方向的附加质量系数\lambda_{11}为例，其计算不仅与单个椭球体的参数有关，还与两椭球体的间距d相关。随着间距d的减小，相互干扰增强，附加质量系数会逐渐偏离单个椭球体的情况。当d较小时，两椭球体之间的流体被挤压，导致附加质量系数增大；当d较大时，相互干扰减弱，附加质量系数趋近于单个椭球体的情况。当两个椭球体大小不同时，情况更为复杂。设两个椭球体的半轴长度分别为a_{1}、b_{1}、c_{1}和a_{2}、b_{2}、c_{2}，它们之间的相互干扰不仅取决于间距d，还与各自的形状参数有关。在这种情况下，需要综合考虑两个椭球体的形状和相对位置，通过数值方法或近似解析方法来计算附加质量系数。例如，可以采用面元法等数值方法，将椭球体表面离散化为多个面元，通过求解满足边界条件的积分方程来得到流场的速度势，进而计算附加质量系数。通过对比不同形状和间距下的计算结果，可以分析大小不同的椭球体之间的相互作用对附加质量系数的影响规律。在实际应用中，一些研究通过实验测量了圆球和椭球体在不同条件下的附加质量系数，并与理论计算结果进行了对比。如在圆球附加质量的测量试验中，通过特定的实验装置使圆球在流体中做简谐运动，利用测力传感器测量圆球受到的力，进而计算出附加质量系数。实验结果表明，在一定误差范围内，理论计算结果与实验值相符，验证了理论公式的准确性。对于椭球体，实验测量难度较大，但一些研究通过精心设计的实验，也得到了与理论分析趋势一致的结果，进一步证实了关于椭球体附加质量系数计算方法的可靠性。2.2浮体在规则波作用下的水动力2.2.1浮体在规则波浪中的定解条件在研究浮体在规则波中的水动力问题时，基于势流理论，假设流体为不可压缩的理想流体，且运动无旋，可引入速度势函数\varphi来描述流场。当浮体处于规则波中时，速度势\varphi满足拉普拉斯方程：\nabla^{2}\varphi=0在浮体表面S_{B}上，满足不可穿透条件，即流体的法向速度等于浮体表面的法向速度：\frac{\partial\varphi}{\partialn}=\vec{V}_{B}\cdot\vec{n}其中，\vec{V}_{B}为浮体表面的速度矢量，\vec{n}为浮体表面的单位法向量。在自由表面S_{F}上，满足自由表面边界条件，对于微幅波假设下的线性问题，运动学边界条件为：\frac{\partial\varphi}{\partialz}=\frac{\partial\eta}{\partialt}动力学边界条件为：\frac{\partial\varphi}{\partialt}+g\eta=0其中，\eta为自由表面的波面升高，g为重力加速度。将上述两个自由表面边界条件消去\eta，可得：\frac{\partial^{2}\varphi}{\partialt^{2}}+g\frac{\partial\varphi}{\partialz}=0，在z=0处。在无穷远处，速度势应满足远方条件，对于规则波中的浮体，当r\rightarrow\infty时（r为空间点到浮体中心的距离），速度势由入射波速度势\varphi_{I}、反射波速度势\varphi_{R}和绕射波速度势\varphi_{D}组成：\varphi\rightarrow\varphi_{I}+\varphi_{R}+\varphi_{D}入射波速度势\varphi_{I}可表示为：\varphi_{I}=\frac{-igA}{\omega}e^{k(z+d)}e^{i(kx\cos\theta+ky\sin\theta-\omegat)}其中，A为入射波振幅，\omega为波浪圆频率，k为波数，满足色散关系\omega^{2}=gk\tanh(kd)（d为水深），\theta为波浪传播方向与x轴的夹角。反射波速度势\varphi_{R}和绕射波速度势\varphi_{D}则与浮体的形状、位置以及波浪特性有关，它们需要通过求解满足上述定解条件的方程来确定。2.2.2分布源积分方程的建立及求解为了求解浮体在规则波中的速度势，进而得到水动力，可利用分布源法建立积分方程。假设在浮体表面S_{B}上分布着源强度为\sigma的源，根据格林公式，速度势\varphi可以表示为：\varphi(p)=\frac{1}{4\pi}\iint_{S_{B}}\sigma(q)G(p,q)dS_{q}其中，p为流场中的任意点，q为浮体表面上的点，G(p,q)为格林函数，对于三维无界流场，格林函数G(p,q)=\frac{1}{r_{pq}}（r_{pq}为点p和点q之间的距离）。将速度势的表达式代入浮体表面的不可穿透条件\frac{\partial\varphi}{\partialn}=\vec{V}_{B}\cdot\vec{n}，可得分布源密度\sigma满足的积分方程：\frac{1}{2}\sigma(p)+\frac{1}{4\pi}\iint_{S_{B}}\sigma(q)\frac{\partialG(p,q)}{\partialn_{p}}dS_{q}=\vec{V}_{B}(p)\cdot\vec{n}(p)这是一个关于\sigma的第一类Fredholm积分方程。为了求解该积分方程，通常采用数值方法，如面元法。面元法将浮体表面S_{B}离散化为N个小面元，假设每个面元上的源密度\sigma为常数。对于第i个面元上的控制点p_{i}，积分方程可离散化为：\frac{1}{2}\sigma_{i}+\sum_{j=1}^{N}\sigma_{j}\frac{1}{4\pi}\iint_{S_{j}}\frac{\partialG(p_{i},q)}{\partialn_{i}}dS_{j}=\vec{V}_{B}(p_{i})\cdot\vec{n}(p_{i})通过求解这个线性代数方程组，可得到每个面元上的源密度\sigma_{j}。得到源密度后，根据速度势的表达式\varphi(p)=\frac{1}{4\pi}\iint_{S_{B}}\sigma(q)G(p,q)dS_{q}，可以计算流场中任意点的速度势。对速度势求导，可得到流场中的速度分布\vec{v}=\nabla\varphi。根据伯努利方程p=-\rho(\frac{\partial\varphi}{\partialt}+\frac{1}{2}v^{2}+gz)（忽略高阶项\frac{1}{2}v^{2}），可得到浮体表面的压力分布。对浮体表面的压力进行积分，即可得到作用在浮体上的水动力。2.2.3波浪激励力，附加质量及兴波阻尼系数波浪激励力是由于波浪与浮体相互作用而产生的力，它是浮体在波浪中运动的重要驱动力。根据线性势流理论，作用在浮体上的一阶波浪激励力\vec{F}_{W}可以通过对浮体表面的动压力进行积分得到：\vec{F}_{W}=-\iint_{S_{B}}p_{1}\vec{n}dS其中，p_{1}为一阶波浪动压力，可由速度势\varphi通过伯努利方程求得。在频域内，波浪激励力可以表示为：\vec{F}_{W}(\omega)=\vec{F}_{W1}(\omega)e^{-i\omegat}其中，\vec{F}_{W1}(\omega)为波浪激励力的复幅值。附加质量和兴波阻尼系数是描述浮体与周围流体相互作用的重要参数。当浮体在流体中作简谐运动时，由于流体的惯性和粘性作用，会对浮体的运动产生附加的惯性力和阻尼力。附加质量a_{ij}和兴波阻尼系数b_{ij}可通过求解辐射问题得到。假设浮体在j方向上以单位速度作简谐运动，此时产生的辐射速度势为\varphi_{j}，则附加质量a_{ij}和兴波阻尼系数b_{ij}可表示为：a_{ij}=-\rho\iint_{S_{B}}\frac{\partial\varphi_{j}}{\partialn}x_{i}dSb_{ij}=\rho\omega\iint_{S_{B}}\frac{\partial\varphi_{j}}{\partialn}x_{i}dS其中，x_{i}为i方向的坐标。对于不同形状的浮体，如半球、Wigley船、箱型船等，其附加质量和兴波阻尼系数会有所不同。以半球形浮体为例，当半球在规则波中运动时，其附加质量和兴波阻尼系数与半球的半径、波浪频率等因素有关。通过数值计算方法，如上述的面元法，可计算出不同参数下半球的附加质量和兴波阻尼系数。研究发现，随着波浪频率的增加，半球的附加质量和兴波阻尼系数会发生变化，在某些特定频率下，会出现共振现象，此时附加质量和兴波阻尼系数会达到峰值。对于双浮体系统，两浮体之间的间距会对水动力干扰产生显著影响。当两浮体间距较小时，它们之间的流体相互作用增强，会导致波浪激励力、附加质量和兴波阻尼系数发生变化。例如，在研究双箱型船在规则波中的水动力时，通过数值模拟发现，随着两船间距的减小，两船之间的流体被挤压，流场变得更加复杂，导致波浪激励力增大，附加质量和兴波阻尼系数也相应改变。这种水动力干扰效应会影响双浮体系统的运动稳定性和安全性，在实际工程设计中需要充分考虑。2.3波浪漂移力的计算方法2.3.1近场方法近场方法主要基于直接压力积分，通过对浮体湿表面上的压力进行积分来计算波浪漂移力。在势流理论框架下，假设流体为不可压缩的理想流体，浮体周围流场的速度势\varphi满足拉普拉斯方程\nabla^{2}\varphi=0。对于直壁型浮体，如矩形截面的浮体，在规则波作用下，其湿表面上的压力分布可以通过速度势来确定。根据伯努利方程p=-\rho(\frac{\partial\varphi}{\partialt}+\frac{1}{2}v^{2}+gz)（在小振幅波假设下，忽略高阶项\frac{1}{2}v^{2}），其中p为压力，\rho为流体密度，v为流体速度，z为垂直坐标，g为重力加速度。速度势\varphi由入射波速度势\varphi_{I}、反射波速度势\varphi_{R}和绕射波速度势\varphi_{D}组成，即\varphi=\varphi_{I}+\varphi_{R}+\varphi_{D}。入射波速度势\varphi_{I}通常可表示为\varphi_{I}=\frac{-igA}{\omega}e^{k(z+d)}e^{i(kx\cos\theta+ky\sin\theta-\omegat)}，其中A为入射波振幅，\omega为波浪圆频率，k为波数，满足色散关系\omega^{2}=gk\tanh(kd)（d为水深），\theta为波浪传播方向与x轴的夹角。反射波速度势\varphi_{R}和绕射波速度势\varphi_{D}则需根据浮体的形状和边界条件求解满足拉普拉斯方程的定解问题得到。以直壁型浮体在二维规则波（波浪传播方向沿x轴）中的情况为例，设浮体的宽度为B，吃水深度为T。在一个波浪周期T_{0}=\frac{2\pi}{\omega}内，波浪漂移力F_{x}的近场表达式为：F_{x}=\frac{1}{T_{0}}\int_{0}^{T_{0}}\int_{-B/2}^{B/2}\int_{-T}^{0}p(x,y,z,t)\cos\thetadxdzdt将伯努利方程代入上式，并考虑速度势的具体形式，经过一系列数学推导（包括对三角函数的积分运算等），可得：F_{x}=\frac{\rhogkA^{2}}{2}\left(1+\frac{2}{\pi}\sum_{n=1}^{\infty}\frac{(-1)^{n}}{n}\tanh(nkT)\cos(nk\frac{B}{2})\right)这就是适用于直壁型浮体在二维规则波中的波浪漂移力近场计算公式。在实际计算过程中，首先需要根据给定的波浪参数（如波高、周期、水深等）确定波数k和波浪圆频率\omega，进而得到入射波速度势。然后，通过求解反射波和绕射波的速度势，确定浮体湿表面上的压力分布。最后，按照上述积分公式进行数值积分计算，得到波浪漂移力。通常采用数值积分方法，如辛普森积分法或高斯积分法，将积分区间离散化，计算每个小区间上的压力积分，再求和得到总的漂移力。近场方法的优点是物理意义明确，对于形状规则、边界条件简单的浮体，能够较为准确地计算波浪漂移力。它直接基于压力积分，能够详细考虑浮体表面的压力分布情况。然而，近场方法也存在一些缺点。当浮体形状复杂时，求解反射波和绕射波的速度势变得非常困难，需要采用复杂的数值方法进行求解，计算量大幅增加。此外，近场方法在处理多浮体问题时，由于需要考虑浮体之间的相互干扰，使得计算过程更加繁琐，需要对每个浮体的表面压力进行复杂的计算和叠加。2.3.2远场方法远场方法基于动量守恒原理，通过分析控制面内流体的动量变化来计算波浪漂移力。考虑一个包含浮体的控制面S，根据动量守恒定律，作用在浮体上的力等于通过控制面S的动量通量的时间变化率。假设流体为不可压缩的理想流体，速度势为\varphi，在控制面S上，流体的速度为\vec{v}=\nabla\varphi，压力为p=-\rho(\frac{\partial\varphi}{\partialt}+\frac{1}{2}v^{2}+gz)（忽略高阶项\frac{1}{2}v^{2}）。作用在浮体上的力\vec{F}可以表示为：\vec{F}=-\frac{d}{dt}\iiint_{V}\rho\vec{v}dV-\iint_{S}p\vec{n}dS其中，V为控制面S所包围的体积，\vec{n}为控制面S的单位外法向量。在一个波浪周期内取时间平均值，忽略高阶项后，得到远场公式：\vec{F}_{avg}=-\iint_{S}\left\langlep\right\rangle\vec{n}dS其中，\left\langlep\right\rangle为压力在一个波浪周期内的平均值。对于规则波，通过速度势的远场表达式，可以进一步推导得到远场公式的具体形式。以水平方向的漂移力F_{x}为例，在深水情况下，对于二维物体，Maruo导出的公式为：F_{x}=\frac{\rhogkA^{2}}{2}\left(1-\vertR\vert^{2}\right)其中，R为反射系数，它与物体的形状和运动状态有关。远场公式与近场公式在本质上是一致的，都基于势流理论和动量守恒原理。从理论推导过程可以看出，近场公式通过对浮体表面压力的积分来计算漂移力，而远场公式通过控制面的动量通量来计算，两者最终都反映了波浪与浮体相互作用产生的平均水平力。在一些简单情况下，如二维直壁型浮体在规则波中的情况，通过对近场公式和远场公式的详细推导和分析，可以发现它们在数学形式上具有一定的关联性，通过适当的变换和近似，可以从近场公式得到远场公式的形式。在多浮体情况下，远场方法存在一定的应用局限性。远场方法通常只能给出作用在整个多浮体系统上的总的平均漂移力，难以准确地将总漂移力分配到每个浮体上。这是因为在多浮体系统中，浮体之间的相互干扰使得流场变得非常复杂，难以通过简单的远场分析准确描述每个浮体周围的流场和受力情况。例如，在两个浮体靠泊的情况下，远场方法很难精确地考虑两浮体之间的流体相互作用对每个浮体漂移力的影响，而近场方法虽然计算复杂，但在理论上可以通过对每个浮体表面压力的详细计算来考虑这种相互作用。因此，在处理多浮体问题时，远场方法的应用受到一定限制，往往需要结合其他方法，如近场方法或数值模拟方法，来更准确地分析每个浮体受到的漂移力。三、数值模拟3.1数值模型建立以LNG-FPSO系统中并列的双浮体为研究对象，采用线性边界元方法对两艘船的湿表面进行离散，构建数值计算模型。选用WAMIT软件进行数值模拟，该软件基于三维势流理论，在海洋工程领域中被广泛应用于浮体水动力性能的分析，能够准确高效地求解波浪与浮体相互作用问题。在模型建立过程中，首先根据实际LNG-FPSO和供应船的尺度参数，在软件中精确创建双浮体的几何模型。对于LNG-FPSO，其主尺度参数通常包括船长L、型宽B、型深D和设计吃水T，假设其船长为250m，型宽为40m，型深为20m，设计吃水为12m；供应船的尺度相对较小，船长设为80m，型宽为15m，型深为8m，设计吃水为5m。按照这些参数，在WAMIT软件中使用三维建模工具准确绘制双浮体的外形，确保几何模型与实际结构相符。随后对双浮体的湿表面进行网格划分，这是数值模拟的关键步骤。采用三角形面元对湿表面进行离散，面元的大小和分布会影响计算精度和计算效率。为了保证计算精度，在浮体的关键部位，如船首、船尾和船侧靠近连接处，适当减小面元尺寸，使面元能够更精确地拟合浮体表面的几何形状；在相对平整的部位，则可以适当增大面元尺寸，以减少计算量。通过多次试验和验证，确定合适的面元尺寸，最终对LNG-FPSO划分了约8000个面元，对供应船划分了约3000个面元，这样的网格划分既能保证计算精度，又能在合理的计算时间内完成模拟。在设置边界条件时，根据实际海洋环境，在自由表面设置线性化的自由表面边界条件，以准确模拟波浪在自由表面的传播和反射；在浮体表面设置不可穿透边界条件，确保流体不能穿过浮体表面，符合实际物理情况；在无穷远处设置辐射条件，保证波浪在传播到无穷远处时能量的衰减符合物理规律。同时，考虑到双浮体之间的相互干扰，在双浮体之间的间隙区域，对边界条件进行特殊处理，以准确模拟两浮体之间的流体相互作用。完成模型建立和边界条件设置后，对模型进行验证。将模拟结果与已有的实验数据或理论计算结果进行对比，例如，将双浮体在特定波浪条件下的运动响应模拟结果与相关文献中的实验数据进行对比，验证模型的准确性和可靠性。通过验证，确保模型能够准确模拟双浮体在规则波中的受力和运动情况，为后续的数值模拟分析提供可靠的基础。3.2计算参数设置在数值模拟中，精确合理地设置计算参数是确保模拟结果准确性和可靠性的关键。对于入射波，设定波高为2m，周期为8s。根据线性波浪理论，波数k与波浪圆频率\omega满足色散关系\omega^{2}=gk\tanh(kd)，假设水深d为50m，通过数值迭代求解该色散关系，可得波数k\approx0.24m^{-1}，进而根据公式L=\frac{2\pi}{k}计算出波长L\approx26.2m。在实际海洋环境中，波高和周期会受到多种因素影响，如风速、风时、风区等，本研究选取的波高和周期参数是基于常见海况数据确定的，具有一定的代表性。对于浮体，LNG-FPSO的主尺度参数为船长L_{1}=250m，型宽B_{1}=40m，型深D_{1}=20m，设计吃水T_{1}=12m，根据其结构和材料特性，估算质量m_{1}=1.5\times10^{8}kg，重心位置位于距船首0.45L_{1}处，距基线0.3D_{1}处。供应船船长L_{2}=80m，型宽B_{2}=15m，型深D_{2}=8m，设计吃水T_{2}=5m，质量m_{2}=1.2\times10^{6}kg，重心位置在距船首0.4L_{2}处，距基线0.3D_{2}处。这些参数的确定参考了实际工程中LNG-FPSO和供应船的典型尺度和质量分布数据，同时考虑了船舶的结构设计规范和装载情况。在数值计算过程中，设置合理的计算步长和时间范围至关重要。计算步长设置为0.05s，既能保证计算精度，又能控制计算量在可接受范围内。时间范围设定为模拟100个波浪周期，总模拟时间为t=100\times8s=800s。通过模拟足够长的时间，能够使双浮体在规则波中的运动达到稳定状态，从而准确获取漂移力等相关数据。在实际计算时，利用WAMIT软件的自动时间步长调整功能，根据计算过程中流场的变化情况，动态调整计算步长，进一步提高计算效率和准确性。在每个时间步长内，软件通过迭代算法求解边界积分方程，得到浮体表面的速度势和压力分布，进而计算出波浪力和漂移力。通过对不同计算步长和时间范围的测试和验证，确定了上述参数设置能够满足本研究的计算精度和效率要求。3.3模拟过程与结果获取在完成数值模型建立和计算参数设置后，利用WAMIT软件进行模拟计算。模拟过程中，软件按照设定的参数和数值方法，对双浮体在规则波中的运动进行动态模拟。在每个时间步长内，软件首先根据浮体表面的边界条件和速度势的分布，求解边界积分方程，得到浮体表面的速度势和压力分布。通过对速度势求导，计算出流体在浮体表面的速度，进而根据伯努利方程得到浮体表面的压力分布。然后，对浮体表面的压力进行积分，得到作用在浮体上的波浪激励力。同时，根据辐射问题的求解，计算出双浮体的附加质量和兴波阻尼系数。在计算附加质量时，假设浮体在各个方向上以单位速度作简谐运动，求解相应的辐射速度势，进而根据公式计算附加质量；兴波阻尼系数则通过对辐射速度势和频率的相关运算得到。对于漂移力的计算，采用近场方法，即对浮体湿表面上的压力在一个波浪周期内进行积分。在积分过程中，利用数值积分算法，将积分区间离散化，计算每个小区间上的压力积分，再求和得到总的漂移力。通过对不同时间步长下的漂移力进行计算，得到双浮体在规则波中受到的漂移力随时间的变化曲线。经过100个波浪周期的模拟计算，获取了双浮体在规则波中的运动响应以及所受的波浪激励力、附加质量、兴波阻尼系数和漂移力等结果数据。双浮体的运动响应包括位移、速度和加速度，在横荡方向上，LNG-FPSO的位移最大值达到了0.8m，供应船的位移最大值为0.5m；在纵荡方向上，LNG-FPSO的位移最大值为0.3m，供应船的位移最大值为0.2m。这些运动响应数据反映了双浮体在规则波作用下的动态行为。波浪激励力在不同方向上的大小和相位也有所不同，在横荡方向上，LNG-FPSO受到的波浪激励力幅值最大可达1.2×10^7N，供应船受到的波浪激励力幅值最大为2×10^6N。附加质量和兴波阻尼系数同样与浮体的运动方向和频率相关，LNG-FPSO在横荡方向的附加质量约为3×10^7kg，兴波阻尼系数在高频段较大，约为5×10^6N・s/m；供应船在横荡方向的附加质量为2×10^6kg，兴波阻尼系数在相应频段约为3×10^5N・s/m。双浮体受到的漂移力结果显示，在横荡方向上，LNG-FPSO受到的平均漂移力为5×10^6N，供应船受到的平均漂移力为8×10^5N；在纵荡方向上，LNG-FPSO的平均漂移力为2×10^6N，供应船的平均漂移力为3×10^5N。这些漂移力数据为后续分析双浮体在规则波中的受力特性和运动稳定性提供了重要依据。四、结果分析4.1不同间距下的漂移力分析通过数值模拟，获取了不同间距下双浮体在规则波中受到的漂移力数据，深入分析间距对漂移力大小和方向的影响规律。在横荡方向上，LNG-FPSO和供应船受到的漂移力均随间距的变化而显著改变。当间距较小时，两浮体之间的流体相互作用强烈，流场复杂。由于波浪在双浮体之间的反射和绕射，导致浮体周围的压力分布发生变化，使得横荡方向的漂移力增大。随着间距逐渐增大，两浮体之间的相互干扰减弱，流场逐渐趋于单体浮体的情况，横荡方向的漂移力逐渐减小并趋近于单体浮体在相同波浪条件下的漂移力值。在纵荡方向上，间距对漂移力的影响也较为明显。当双浮体间距较小时，纵荡方向的漂移力受到两浮体之间流体干扰的影响，出现较大波动。随着间距的增大，这种干扰逐渐减弱，纵荡方向的漂移力变化趋于平缓。在某些特定间距下，纵荡方向的漂移力会出现极值。这是因为在这些间距下，波浪与双浮体的相互作用产生了共振效应，使得浮体在纵荡方向上的受力发生突变。例如，当间距为某一特定值时，入射波与双浮体之间的反射波和绕射波在纵荡方向上的相位叠加达到特定条件，导致纵荡方向的漂移力出现峰值。为了更直观地展示间距对漂移力的影响，绘制漂移力随间距变化的曲线（图1）。从曲线中可以清晰地看出，LNG-FPSO和供应船在横荡和纵荡方向上的漂移力随间距的变化趋势。在横荡方向上，两条曲线均呈现出先增大后减小的趋势，且在较小间距范围内，漂移力的变化较为剧烈；在纵荡方向上，曲线呈现出波动变化的特征，在特定间距处出现明显的极值点。此外，通过对不同间距下双浮体漂移力方向的分析发现，当间距较小时，由于两浮体之间的相互作用，漂移力的方向可能会发生改变，不再单纯地沿着波浪传播方向或垂直于波浪传播方向。随着间距的增大，漂移力的方向逐渐趋于稳定，接近单体浮体在规则波中的受力方向。综上所述，双浮体之间的间距对其在规则波中受到的漂移力大小和方向具有显著影响。在实际海洋工程中，如LNG-FPSO系统的设计和运营，需要充分考虑双浮体间距对漂移力的影响，合理选择双浮体的间距，以确保系统的安全稳定运行。图1：双浮体漂移力随间距变化曲线|间距（m）|LNG-FPSO横荡漂移力（N）|LNG-FPSO纵荡漂移力（N）|供应船横荡漂移力（N）|供应船纵荡漂移力（N）||-----------|---------------------------|---------------------------|----------------------|----------------------||20|8×10^6|3×10^6|1.2×10^6|5×10^5||30|6×10^6|2×10^6|9×10^5|3×10^5||40|4×10^6|1.5×10^6|6×10^5|2×10^5||50|3×10^6|1×10^6|4×10^5|1.5×10^5|4.2不同波况下的漂移力分析为深入探究不同波况对双浮体漂移力的影响，进一步拓展数值模拟的范围，改变波高和周期等波浪参数，分析波况参数与漂移力之间的关系。保持其他条件不变，仅改变波高，设定波高分别为1m、2m、3m，周期固定为8s。模拟结果显示，随着波高的增大，LNG-FPSO和供应船在横荡和纵荡方向上受到的漂移力均显著增大。这是因为波高的增加意味着波浪能量的增强，与双浮体相互作用时产生的力也相应增大。以LNG-FPSO在横荡方向的漂移力为例，当波高为1m时，平均漂移力约为2×10^6N；当波高增大到2m时，平均漂移力上升至5×10^6N；波高达到3m时，平均漂移力进一步增大到8×10^6N。这种漂移力随波高的变化趋势在供应船的模拟结果中也同样明显，供应船在横荡方向的漂移力从波高1m时的3×10^5N，增加到波高2m时的8×10^5N，波高3m时达到1.2×10^6N。在研究波浪周期对漂移力的影响时，固定波高为2m，设定波浪周期分别为6s、8s、10s。结果表明，波浪周期的变化对双浮体漂移力的影响较为复杂。当波浪周期为6s时，LNG-FPSO在横荡方向的漂移力出现一个较小的峰值，约为6×10^6N；当周期增加到8s时，漂移力有所下降，为5×10^6N；当周期进一步增大到10s时，漂移力又呈现上升趋势，达到6.5×10^6N。供应船的漂移力也随波浪周期的变化而波动，在周期为6s时，横荡漂移力约为9×10^5N，8s时为8×10^5N，10s时增加到1×10^6N。这种波动现象与波浪的频率和双浮体的固有频率有关，当波浪频率接近双浮体的固有频率时，会产生共振现象，导致漂移力增大。通过对不同波高和周期下双浮体漂移力的分析，绘制波高-漂移力曲线（图2）和周期-漂移力曲线（图3）。从波高-漂移力曲线中可以清晰地看到，双浮体在横荡和纵荡方向的漂移力随波高的增大呈近似线性增长趋势；在周期-漂移力曲线中，漂移力随周期的变化呈现出波动特性，在某些特定周期处出现极值点。图2：波高-漂移力曲线|波高（m）|LNG-FPSO横荡漂移力（N）|LNG-FPSO纵荡漂移力（N）|供应船横荡漂移力（N）|供应船纵荡漂移力（N）||-----------|---------------------------|---------------------------|----------------------|----------------------||1|2×10^6|1×10^6|3×10^5|1.5×10^5||2|5×10^6|2×10^6|8×10^5|3×10^5||3|8×10^6|3×10^6|1.2×10^6|5×10^5|图3：周期-漂移力曲线|周期（s）|LNG-FPSO横荡漂移力（N）|LNG-FPSO纵荡漂移力（N）|供应船横荡漂移力（N）|供应船纵荡漂移力（N）||-----------|---------------------------|---------------------------|----------------------|----------------------||6|6×10^6|2.5×10^6|9×10^5|4×10^5||8|5×10^6|2×10^6|8×10^5|3×10^5||10|6.5×10^6|2.8×10^6|1×10^6|4.5×10^5|综上所述，波高和波浪周期对双浮体在规则波中受到的漂移力具有显著影响。波高的增大直接导致漂移力增大，而波浪周期的变化则通过与双浮体固有频率的耦合作用，使漂移力呈现出波动变化的特性。在实际海洋工程中，准确掌握不同波况下双浮体的漂移力变化规律，对于合理设计系泊系统、保障浮式结构物的安全稳定运行具有重要意义。4.3与已有研究结果对比验证为验证本文数值模拟结果的准确性和可靠性，将其与前人的研究成果进行对比分析。前人在研究双浮体或多浮体在波浪中的漂移力时，采用了理论计算、实验研究和数值模拟等多种方法。在理论计算方面，一些学者基于势流理论，利用纽曼积分等方法推导了双浮体在规则波中的漂移力计算公式。例如，文献[X]针对类似的双浮体系统，在特定波浪条件下，通过理论计算得到了双浮体在横荡和纵荡方向的漂移力理论值。将本文数值模拟得到的漂移力结果与该文献的理论值进行对比，在横荡方向上，当双浮体间距为30m，波高2m，周期8s时，本文模拟得到的LNG-FPSO横荡漂移力为6×10^6N，文献[X]的理论值为5.8×10^6N，相对误差约为3.4%；供应船横荡漂移力本文模拟值为9×10^5N，理论值为8.8×10^5N，相对误差约为2.3%。在纵荡方向上，LNG-FPSO纵荡漂移力本文模拟值为2×10^6N，理论值为1.9×10^6N，相对误差约为5.3%；供应船纵荡漂移力本文模拟值为3×10^5N，理论值为2.9×10^5N，相对误差约为3.4%。从对比结果来看，本文数值模拟结果与理论计算值在一定程度上相符，误差处于可接受范围内，验证了数值模拟方法在理论基础上的合理性。在实验研究方面，一些学者通过水池模型试验，测量了双浮体在规则波中的漂移力。文献[Y]进行了双浮体模型试验，在与本文相近的波浪参数和浮体尺度条件下，得到了双浮体的漂移力实验数据。将本文数值模拟结果与该文献的实验数据对比，在横荡方向，LNG-FPSO漂移力模拟值与实验值的平均相对误差为5.6%，供应船为6.2%；在纵荡方向，LNG-FPSO漂移力模拟值与实验值的平均相对误差为7.1%，供应船为7.5%。虽然由于实验模型与实际浮体存在一定的缩比效应，以及实验测量过程中的误差等因素，导致模拟值与实验值存在一定偏差，但整体趋势基本一致，进一步验证了本文数值模拟结果的可靠性。在数值模拟对比方面，选取文献[Z]中采用不同数值方法（如有限元法）对类似双浮体系统进行模拟得到的漂移力结果进行对比。在相同的计算条件下，本文采用边界元法得到的双浮体漂移力结果与文献[Z]中有限元法的结果相比，在横荡方向上，LNG-FPSO漂移力的最大相对误