一种漂浮式风力机流体动力学耦合分析方法及装置与流程

1.本技术涉及流体动力学仿真模拟技术领域，更具体的说，是涉及一种漂浮式风力机流体动力学耦合分析方法及装置。


背景技术：

2.海上风能相对于陆地风能，风速更高、风能稳定性更好，开发海上风能可以有效的缓解我国沿海城市的用电压力。随着海上风电朝着深远海域发展，海上漂浮式风力机成为解决深远海风能利用的主要发电设备。
3.由于海上漂浮式风力机受到风、浪、流联合作用，且风的强非线性变化导致海上漂浮式风力机的空气动力学载荷存在较大的非线性变化，在控制漂浮式风力机的运动时，需要进行时域范围下的流体动力学耦合分析，才能保证运动控制的合理性和安全性。目前，针对海上漂浮式风力机的流体动力学时域耦合分析方法，通常采用黏性计算流体力学(computational fluiddynamics,cfd)方法，该方法计算精度高。
4.然而，传统cfd方法在进行仿真计算时，需要围绕漂浮式风力机几何模型构建数量庞大的网格节点，且对每个网格节点需要求解纳维-斯托克斯方程，耗费时间巨大，无法实时进行流体动力学耦合分析。
5.如何实现对漂浮式风力机进行高效流体动力学耦合分析，是需要关注的问题。


技术实现要素：

6.鉴于上述问题，提出了本技术以便提供一种漂浮式风力机流体动力学耦合分析方法及装置，以提高漂浮式风力机的流体动力学耦合分析的效率。
7.为了实现上述目的，现提出具体方案如下：
8.一种漂浮式风力机流体动力学耦合分析方法，包括：
9.分析漂浮式风力机的浮式基础，确定并存储所述漂浮式风力机在预设时长的当前时间步下的水动力载荷，以及确定并存储静水回复力矩阵；
10.分析所述漂浮式风力机中缆绳的当前状态信息，确定并存储所述当前时间步下所述漂浮式风力机的系泊力载荷；
11.基于所述漂浮式风力机在所述当前时间步下的机体参数信息和/或所处于水体的水体参数信息，确定并存储所述漂浮式风力机在所述当前时间步下的静态浮力载荷、黏性力载荷；
12.获取所述漂浮式风力机的叶片气动参数，以及所述漂浮式风力机的各气动节点的坐标和运动信息；
13.基于所述叶片气动参数与所述各气动节点的坐标和运动信息，通过叶素动量理论及其修正算法，计算并存储所述漂浮式风力机在当前时间步下的气动载荷；
14.基于所述漂浮式风力机的当前位置，融合所述气动载荷、所述系泊力载荷、所述水动力载荷、所述静态浮力载荷、所述黏性力载荷以及所述静水回复力矩阵，确定、更新并存
储所述漂浮式风力机在下一时间步的坐标及运动信息，并将下一时间步作为当前时间步，返回执行分析漂浮式风力机的浮式基础，确定并存储所述漂浮式风力机在预设时长的当前时间步下的水动力载荷的步骤，直到达到预设时间步为止。
15.可选的，所述分析漂浮式风力机的浮式基础，确定所述漂浮式风力机在预设时长的当前时间步下的水动力载荷，包括：
16.创建漂浮式风力机的浮式基础的几何模型，对所述几何模型进行网格划分，得到网格划分后的几何模型；
17.对所述网格划分后的几何模型进行水动力学分析，得到附加质量矩阵、阻尼系数和幅值响应算子；
18.基于所述阻尼系数，确定所述漂浮式风力机在预设时长的当前时间步下的阻尼系数函数；
19.利用jonswap波浪谱，创建模拟所述漂浮式风力机当前航行时的波浪函数，并确定所述波浪函数的幅值、相位；
20.根据所述波浪函数的幅值、相位以及所述幅值响应算子，确定所述漂浮式风力机的波浪激励力载荷，并将所述波浪激励力载荷、所述阻尼系数函数以及所述附加质量矩阵组成水动力载荷。
21.可选的，所述幅值响应算子包括响应幅值和响应相位；
22.根据所述波浪函数的幅值、相位以及所述幅值响应算子，确定所述漂浮式风力机的波浪激励力载荷，包括：
23.利用如下函数，确定所述漂浮式风力机的波浪激励力载荷：
[0024][0025]
其中，|h(ωj)|为所述响应幅值，δ(ωj)为所述响应相位，aj为所述波浪函数的幅值，t为所述预设时长，ωjt εj为所述波浪函数的相位，n为所述波浪函数的预设划分次数。
[0026]
可选的，分析所述漂浮式风力机中缆绳的当前状态信息，确定所述当前时间步下所述漂浮式风力机的系泊力载荷，包括：
[0027]
分析所述漂浮式风力机中缆绳的当前状态信息，得到所述缆绳的长度，所述缆绳在水平方向上和竖直方向上的长度，所述缆绳与水中泥沙之间的摩擦系数，所述缆绳的弹性模量、截面积，所述缆绳在水中的重力信息；
[0028]
利用如下方程组，得到所述缆绳的水平分力信息和竖直分力信息：
[0029]
[0030][0031]
其中，xf表示所述缆绳在水平方向上的长度，zf表示所述缆绳在竖直方向上的长度，ω表示所述缆绳在水中的重力，l表示所述缆绳的长度，cb表示所述摩擦系数，e表示所述缆绳的弹性模量，a表示所述缆绳的截面积，hf表示所述缆绳的水平分力信息，vf表示所述缆绳的竖直分力信息；
[0032]
根据所述水平分力信息和竖直分力信息，确定所述漂浮式风力机的系泊力载荷信息。
[0033]
可选的，所述机体参数信息包括所述漂浮式风力机稳定时在所述水体下的部分的体积；
[0034]
所述水体参数信息包括水体密度；
[0035]
基于所述漂浮式风力机在所述当前时间步下的机体参数信息和/或所处于水体的水体参数信息，确定所述漂浮式风力机在所述当前时间步下的静态浮力载荷，包括：
[0036]
根据所述体积以及所述水体密度，确定所述漂浮式风力机在所述水体上稳定时的静态浮力载荷。
[0037]
可选的，所述水体参数信息包括水体密度，所述机体参数信息包括所述漂浮式风力机的机体长轴方向上的横截面积；
[0038]
基于所述漂浮式风力机在所述当前时间步下的机体参数信息和/或所处于水体的水体参数信息，确定所述漂浮式风力机在所述当前时间步下的黏性力载荷，包括：
[0039]
基于所述波浪引起的水质点的变化，确定所述水质点与所述漂浮式风力机的机体短轴方向上的相对速度；
[0040]
根据所述水体的密度、所述横截面积、所述相对速度、预设的第一系数，确定所述漂浮式风力机在所述当前时间步下的黏性力载荷。
[0041]
可选的，所述叶片气动参数包括转子叶片的升力、阻力和扭矩系数；
[0042]
基于所述叶片气动参数与所述各气动节点的坐标和运动信息，通过叶素动量理论及其修正算法，计算并存储所述漂浮式风力机在当前时间步下的气动载荷，包括：
[0043]
根据所述各气动节点的坐标和运动信息，确定所述漂浮式风力机中包含各气动节点的塔柱所受到的作用于塔架表面的法向力；
[0044]
根据所述转子叶片的升力、阻力和扭矩系数，确定所述转子叶片的风力信息；
[0045]
根据所述转子叶片的风力信息和所述法向力，确定所述漂浮式风力机在当前时间步下的气动载荷。
[0046]
可选的，所述水动力载荷包括波浪激励力载荷、阻尼系数函数以及附加质量矩阵；
[0047]
基于所述漂浮式风力机的当前位置，融合所述气动载荷、所述系泊力载荷、所述水动力载荷、所述静态浮力载荷、所述黏性力载荷以及所述静水回复力矩阵，确定所述漂浮式风力机在下一时间步的坐标及运动信息，包括：
[0048]
利用如下微分方程，得到所述漂浮式风力机在所述预设时长后的下一时间步的位移量、速度和加速度：
[0049][0050]
其中，[m]表示基于所述漂浮式风力机的机体质量所得到的质量矩阵， [a]
inf
表示所述附加质量矩阵，k(t)表示以所述预设时长为自变量的阻尼系数函数，[c]表示所述静水回复力矩阵，f
waves
表示所述波浪激励力载荷对应的波浪激励力，f
mooring
表示所述系泊力载荷对应的系泊力，f
drag
表示所述黏性力载荷对应的黏性力，f
wind
表示所述气动载荷对应的风力，f
static
表示所述静态浮力载荷对应的静水浮力，t表示所述预设时间，表示加速度，表示速度，q表示位移量，所述位移量、速度和加速度组成运动信息；
[0051]
基于所述漂浮式风力机的当前位置以及所述预设时长后的下一时间步的位移量，确定所述漂浮式风力机在下一时间步的坐标。
[0052]
可选的，在所述返回执行获取所述漂浮式风力机的叶片气动参数，以及所述漂浮式风力机的各气动节点的坐标和运动信息的步骤，直到达到预设时间步为止之后，还包括：
[0053]
汇总所述漂浮式风力机在各时间步下的各载荷、坐标及运动信息，以供对所述漂浮式风力机进行性能预报。
[0054]
一种漂浮式风力机流体动力学耦合分析装置，包括：
[0055]
浮式基础分析单元，用于分析漂浮式风力机的浮式基础，确定并存储所述漂浮式风力机在预设时长的当前时间步下的水动力载荷，以及确定并存储静水回复力矩阵；
[0056]
系泊力载荷确定单元，用于分析所述漂浮式风力机中缆绳的当前状态信息，确定并存储所述当前时间步下所述漂浮式风力机的系泊力载荷；
[0057]
双载荷确定单元，用于基于所述漂浮式风力机在所述当前时间步下的机体参数信息和/或所处于水体的水体参数信息，确定并存储所述漂浮式风力机在所述当前时间步下的静态浮力载荷、黏性力载荷；
[0058]
气动参数获取单元，用于获取所述漂浮式风力机的叶片气动参数，以及所述漂浮式风力机的各气动节点的坐标和运动信息；
[0059]
气动载荷计算单元，用于基于所述叶片气动参数与所述各气动节点的坐标和运动信息，通过叶素动量理论及其修正算法，计算并存储所述漂浮式风力机在当前时间步下的气动载荷；
[0060]
运动信息确定单元，用于基于所述漂浮式风力机的当前位置，融合所述气动载荷、所述系泊力载荷、所述水动力载荷、所述静态浮力载荷、所述黏性力载荷以及所述静水回复力矩阵，确定、更新并存储所述漂浮式风力机在下一时间步的坐标及运动信息，并将下一时间步作为当前时间步，返回执行所述浮式基础分析单元用于分析漂浮式风力机的浮式基础，确定并存储所述漂浮式风力机在预设时长的当前时间步下的水动力载荷的处理步骤，直到达到预设时间步为止。
[0061]
借由上述技术方案，本技术通过分析漂浮式风力机的浮式基础，确定并存储所述漂浮式风力机在预设时长的当前时间步下的水动力载荷以及静水回复力矩阵，分析所述漂浮式风力机中缆绳的当前状态信息，确定并存储所述当前时间步下所述漂浮式风力机的系泊力载荷，基于所述漂浮式风力机在所述当前时间步下的机体参数信息和/或所处于水体的水体参数信息，确定并存储所述漂浮式风力机在所述当前时间步下的静态浮力载荷、黏
性力载荷，获取所述漂浮式风力机的叶片气动参数，以及所述漂浮式风力机的各气动节点的坐标和运动信息，基于所述叶片气动参数与所述各气动节点的坐标和运动信息，通过叶素动量理论及其修正算法，计算并存储所述漂浮式风力机在当前时间步下的气动载荷，基于所述漂浮式风力机的当前位置，融合所述气动载荷、所述系泊力载荷、所述水动力载荷、所述静态浮力载荷、所述黏性力载荷以及所述静水回复力矩阵，确定、更新并存储所述漂浮式风力机在下一时间步的坐标及运动信息，并将下一时间步作为当前时间步，返回执行获取所述漂浮式风力机的叶片气动参数，以及所述漂浮式风力机的各气动节点的坐标和运动信息的步骤，直到达到预设时间步为止。由此可见，通过分析漂浮式风力机的当前时间步下的气动载荷、系泊力载荷、水动力载荷、静态浮力载荷、黏性力载荷以及静水回复力矩阵，并将各参数耦合进行分析，从而得到下一时间步的运动信息及位置，避免了对漂浮式风力机的几何模型构建庞大的网格节点，进行复杂的流场变化分析，提高了对漂浮式风力机流体动力学耦合分析的效率。
附图说明
[0062]
通过阅读下文优选实施方式的详细描述，各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的，而并不认为是对本技术的限制。而且在整个附图中，用相同的参考符号表示相同的部件。在附图中：
[0063]
图1为本技术实施例提供的一种漂浮式风力机流体动力学耦合分析方法的流程图；
[0064]
图2为本技术实施例提供的一种漂浮式风力机流体动力学耦合分析装置的结构示意图；
[0065]
图3为本技术实施例提供的一种漂浮式风力机流体动力学耦合分析设备的结构示意图。
具体实施方式
[0066]
下面将结合本技术实施例中的附图，对本技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本技术一部分实施例，而不是全部的实施例。基于本技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本技术保护的范围。
[0067]
本技术方案可以基于具备数据处理能力的终端实现，该终端可以是电脑、服务器、云端等。本技术方案中的参数计算方法可以基于modelica语言实现。
[0068]
接下来，结合图1所述，本技术的漂浮式风力机流体动力学耦合分析方法可以包括如下步骤：
[0069]
步骤s110、分析漂浮式风力机的浮式基础，确定并存储所述漂浮式风力机在预设时长的当前时间步下的水动力载荷，以及确定并存储静水回复力矩阵。
[0070]
具体的，所述漂浮式风力机的浮式基础可以是漂浮式风力机在水面以下的部分，分析漂浮式风力机的浮式基础可以是对漂浮式风力机在水面以下的部分建立模型。所述预设时长可以表示漂浮式风力机模型从起始位置运动至目标位置的最短时间，预设时长可以自定义。
[0071]
可以理解的是，由于需要对漂浮式风力机进行流体动力学耦合分析时融合多个载荷，因此需要对所分析得到的水动力载荷以及静水回复力矩阵预先存储，以便从内存里提取。
[0072]
步骤s120、分析所述漂浮式风力机中缆绳的当前状态信息，确定并存储所述当前时间步下所述漂浮式风力机的系泊力载荷。
[0073]
具体的，缆绳的当前状态信息可以是缆绳在水平、竖直方向上的长度，缆绳在水中的重力，缆绳的长度，缆绳与水中泥沙之间的摩擦系数，缆绳的弹性模量、截面积等。
[0074]
可以理解的是，由于需要对漂浮式风力机进行流体动力学耦合分析时融合多个载荷，因此需要对所分析得到的系泊力载荷预先存储，以便从内存里提取。
[0075]
步骤s130、基于所述漂浮式风力机在所述当前时间步下的机体参数信息和/或所处于水体的水体参数信息，确定并存储所述漂浮式风力机在所述当前时间步下的静态浮力载荷、黏性力载荷。
[0076]
可以理解的是，由于需要对漂浮式风力机进行流体动力学耦合分析时融合多个载荷，因此需要对所分析得到的静态浮力载荷、黏性力载荷预先存储，以便从内存里提取。
[0077]
步骤s140、获取所述漂浮式风力机的叶片气动参数，以及所述漂浮式风力机的各气动节点的坐标和运动信息。
[0078]
其中，所述坐标可以表示漂浮式风力机的各气动节点的具体位置，运动信息可以表示漂浮式风力机的各气动节点的移动速度和加速度。
[0079]
具体的，对于首个时间步的各气动节点的坐标和运动信息可以初始化处理，初始化处理可以是随机生成的，首个时间步的后续时间步的各气动节点的坐标和运动信息可以是基于，对漂浮式风力机流体动力学耦合分析得到的坐标和运动信息得到的。
[0080]
步骤s150、基于所述叶片气动参数与所述各气动节点的坐标和运动信息，通过叶素动量理论及其修正算法，计算并存储所述漂浮式风力机在当前时间步下的气动载荷。
[0081]
可以理解的是，由于需要对漂浮式风力机进行流体动力学耦合分析时融合多个载荷，因此需要对所分析得到的气动载荷预先存储，以便从内存里提取。
[0082]
步骤s160、基于所述漂浮式风力机的当前位置，融合各载荷以及所述静水回复力矩阵，确定、更新并存储所述漂浮式风力机在下一时间步的坐标及运动信息，并将下一时间步作为当前时间步。
[0083]
具体的，确定漂浮式风力机在下一时间步的坐标及运动信息的方式可以是基于多体动力学模型中的运动控制方程得到的。
[0084]
可以理解的是，由于需要分析多个时间步下的漂浮式风力机的坐标及运动信息，因此需要更新并存储每一时间步下的坐标及运动信息，以便分析后续时间步下漂浮式风力机的各载荷。
[0085]
步骤s170、判断是否达到预设时间步，若是，则结束，若否，则返回执行步骤s110。
[0086]
具体的，预设时间步可以表示需要对漂浮式风力机流体动力学耦合分析的最少次数下的时间步，预设时间步可以自定义。
[0087]
可以理解的是，若当前时间步达到预设时间步，可以表示已经对漂浮式风力机流体动力学耦合分析完毕，此时得到各时间步下的漂浮式风力机的各载荷、坐标以及运动信息，结束对漂浮式风力机流体动力学耦合分析。若当前时间步未达预设时间步，可以表示仍
未完成对漂浮式风力机流体动力学耦合分析，则返回执行步骤s110。
[0088]
本实施例提供的漂浮式风力机的运动控制方法，分析漂浮式风力机的浮式基础，确定并存储所述漂浮式风力机在预设时长的当前时间步下的水动力载荷以及静水回复力矩阵，分析所述漂浮式风力机中缆绳的当前状态信息，确定并存储所述当前时间步下所述漂浮式风力机的系泊力载荷，基于所述漂浮式风力机在所述当前时间步下的机体参数信息和/或所处于水体的水体参数信息，确定并存储所述漂浮式风力机在所述当前时间步下的静态浮力载荷、黏性力载荷，获取所述漂浮式风力机的叶片气动参数，以及所述漂浮式风力机的各气动节点的坐标和运动信息，基于所述叶片气动参数与所述各气动节点的坐标和运动信息，通过叶素动量理论及其修正算法，计算并存储所述漂浮式风力机在当前时间步下的气动载荷，基于所述漂浮式风力机的当前位置，融合所述气动载荷、所述系泊力载荷、所述水动力载荷、所述静态浮力载荷、所述黏性力载荷以及所述静水回复力矩阵，确定、更新并存储所述漂浮式风力机在下一时间步的坐标及运动信息，并将下一时间步作为当前时间步，返回执行获取所述漂浮式风力机的叶片气动参数，以及所述漂浮式风力机的各气动节点的坐标和运动信息的步骤，直到达到预设时间步为止。由此可见，通过分析漂浮式风力机的当前时间步下的气动载荷、系泊力载荷、水动力载荷、静态浮力载荷、黏性力载荷以及静水回复力矩阵，并将各参数耦合进行分析，从而得到下一时间步的运动信息及位置，避免了对漂浮式风力机的几何模型构建庞大的网格节点，进行复杂的流场变化分析，提高了对漂浮式风力机流体动力学耦合分析的效率。
[0089]
本技术的一些实施例中，对上述步骤s110、分析漂浮式风力机的浮式基础，确定所述漂浮式风力机在预设时长的当前时间步下的水动力载荷的过程进行介绍，该过程可以包括：
[0090]
s1、创建漂浮式风力机的浮式基础的几何模型，对所述几何模型进行网格划分，得到网格划分后的几何模型。
[0091]
具体的，漂浮式风力机的浮式基础的几何模型可以是通过modelica中的ansys workbench中hydrodynamic diffraction组件创建的。
[0092]
s2、对所述网格划分后的几何模型进行水动力学分析，得到附加质量矩阵、阻尼系数和幅值响应算子。
[0093]
具体的，幅值响应算子可以由响应幅值和响应相位组成，其中，响应幅值可以表示为|h(ωj)|，响应相位可以表示为δ(ωj)。
[0094]
s3、基于所述阻尼系数，确定所述漂浮式风力机在预设时长的当前时间步下的阻尼系数函数。
[0095]
具体的，阻尼系数函数可以表示为：
[0096][0097]
其中，b
ij
(ω)可以表示阻尼系数，t可以表示预设时间，下标i、j为自由度方向。
[0098]
s4、利用jonswap波浪谱，创建模拟所述漂浮式风力机当前航行时的波浪函数，并确定所述波浪函数的幅值、相位。
[0099]
具体的，所述波浪函数可以表示为：
[0100][0101]
其中，aj可以表示波浪函数的幅值，t可以表示所述预设时长，ωjt εj可以表示波浪函数的相位，n为所述波浪函数划分的次数，表示波浪函数的模拟精度，n可以自定义。
[0102]
s5、根据所述波浪函数的幅值、相位以及所述幅值响应算子，确定所述漂浮式风力机的波浪激励力载荷，并将所述波浪激励力载荷、所述阻尼系数函数以及所述附加质量矩阵组成水动力载荷。
[0103]
具体的，波浪激励力载荷可以表示为：
[0104][0105]
本技术提供的实施例中，能够通过基于漂浮式风力机的浮式基础创建模型，并分析得到波浪激励力载荷、阻尼系数函数以及附加质量矩阵，从而确定漂浮式风力机在当前时间步下的水动力载荷。
[0106]
本技术的一些实施例中，对上述步骤s120、分析所述漂浮式风力机中缆绳的当前状态信息，确定并存储所述当前时间步下所述漂浮式风力机的系泊力载荷的过程进行介绍，该过程可以包括：
[0107]
s1、分析所述漂浮式风力机中缆绳的当前状态信息，得到所述缆绳的长度，所述缆绳在水平方向上和竖直方向上的长度，所述缆绳与水中泥沙之间的摩擦系数，所述缆绳的弹性模量、截面积，所述缆绳在水中的重力信息。
[0108]
s2、确定所述缆绳的水平分力信息和竖直分力信息。
[0109]
具体的，可以利用如下方程组，得到所述缆绳的水平分力信息和竖直分力信息：
[0110][0111][0112]
其中，xf表示所述缆绳在水平方向上的长度，zf表示所述缆绳在竖直方向上的长度，ω表示所述缆绳在水中的重力，l表示所述缆绳的长度，cb表示所述摩擦系数，e表示所述缆绳的弹性模量，a表示所述缆绳的截面积，hf表示所述缆绳的水平分力信息，vf表示所述缆绳的竖直分力信息。
[0113]
s3、根据所述水平分力信息和竖直分力信息，确定所述漂浮式风力机的系泊力载荷信息。
[0114]
可以理解的是，漂浮式风力机包含缆绳，因此作用在缆绳上的水平分力信息对应的水平分力，以及竖直分力信息对应的竖直分力，两者的合力可以作为漂浮式风力机的系
泊力载荷。
[0115]
本技术提供的实施例中，能够通过分析漂浮式风力机中缆绳的当前状态信息，得到缆绳的长度，缆绳在水平方向上和竖直方向上的长度，缆绳与水中泥沙之间的摩擦系数，缆绳的弹性模量、截面积，缆绳在水中的重力信息，从而确定漂浮式风力机在当前时间步下的系泊力载荷。
[0116]
本技术的一些实施例中，当上述实施例提到的机体参数信息为漂浮式风力机稳定时在所述水体下的部分的体积，以及，水体参数信息为水体密度时，对上述实施例提到的基于所述漂浮式风力机在所述当前时间步下的机体参数信息和/或所处于水体的水体参数信息，确定并存储所述漂浮式风力机在所述当前时间步下的静态浮力载荷的过程进行介绍，该过程可以包括：
[0117]
s1、根据所述体积以及所述水体密度，确定所述漂浮式风力机在所述水体上稳定时的静态浮力载荷。
[0118]
具体的，可以预先在漂浮式风力机航行在水体上且稳定时，测量漂浮式风力机在水体下的部分的体积，以及水体密度，确定所述漂浮式风力机此刻的静态浮力载荷信息。也可以预先在漂浮式风力机航行在水体上且稳定时，测量漂浮式风力机在水体下的部分的体积，以及水体密度并存储与内存，当需要计算静浮力时，调用所述体积，以及水体密度并计算所述静态浮力载荷。
[0119]
本技术提供的实施例中，能够通过漂浮式风力机稳定时在所述水体下的部分的体积，以及水体密度，确定漂浮式风力机在当前时间步下的静浮力载荷。
[0120]
本技术的一些实施例中，当上述实施例提到的机体参数信息为漂浮式风力机的机体长轴方向上的横截面积，以及，水体参数信息为水体密度时，对上述实施例提到的基于所述漂浮式风力机在所述当前时间步下的机体参数信息和/或所处于水体的水体参数信息，确定并存储所述漂浮式风力机在所述当前时间步下的黏性力载荷的过程进行介绍，该过程可以包括：
[0121]
s1、基于所述波浪引起的水质点的变化，确定所述水质点与所述漂浮式风力机的机体短轴方向上的相对速度。
[0122]
可以理解的是，当漂浮式风力机航行在平稳的水体上时，水质点位于漂浮式风力机的长轴上，当存在波浪时会使水质点偏离漂浮式风力机的长轴，那么此时的偏离的瞬时速度可以为所述水质点与所述漂浮式风力机的机体短轴方向上的相对速度。
[0123]
s2、根据所述水体的密度、所述横截面积、所述相对速度、预设的第一系数，确定所述漂浮式风力机在所述当前时间步下的黏性力载荷。
[0124]
具体的，黏性力载荷可以表示为：
[0125][0126]
其中，cd可以表示预设的第一系数，预设的第一系数可以表示拖拽力系数，ρ可以表示水体密度，a可以表示漂浮式风力机的机体长轴方向上的横截面积，μ
x
可以表示水质点与漂浮式风力机的机体短轴方向上的相对速度。
[0127]
本技术提供的实施例中，能够通过水体密度，漂浮式风力机的机体长轴方向上的横截面积，以及水质点与漂浮式风力机的机体短轴方向上的相对速度，确定漂浮式风力机
在当前时间步下的黏性力载荷。
[0128]
本技术的一些实施例中，当上述实施例提到的叶片气动参数为转子叶片的升力、阻力和扭矩系数时，对上述步骤s150、基于所述叶片气动参数与所述各气动节点的坐标和运动信息，通过叶素动量理论及其修正算法，计算并存储所述漂浮式风力机在当前时间步下的气动载荷的过程进行介绍，该过程可以包括：
[0129]
s1、根据所述各气动节点的坐标和运动信息，确定所述漂浮式风力机中包含各气动节点的塔柱所受到的作用于塔架表面的法向力。
[0130]
s2、根据所述转子叶片的升力、阻力和扭矩系数，确定所述转子叶片的风力信息。
[0131]
具体的，可以基于转子叶片的升力和阻力计算出合力，将扭矩系数与合力相乘得到转子叶片的风力。
[0132]
s3、根据所述转子叶片的风力信息和所述法向力，确定所述漂浮式风力机在当前时间步下的气动载荷。
[0133]
可以理解的是，所述风力信息对应的风力作用于转子叶片，所述法向力作用于塔柱，漂浮式风力机可以包含转子叶片和塔柱，因此所述风力与法向力的合力可以作为漂浮式风力机的风力载荷信息。
[0134]
本技术提供的实施例中，能够通过漂浮式风力机的各气动节点的坐标和运动信息确定塔柱所受到的作用于塔架表面的法向力，并根据所述转子叶片的升力、阻力和扭矩系数，确定所述转子叶片的风力信息，从而确定漂浮式风力机在当前时间步下的风力载荷。
[0135]
本技术的一些实施例中，当上述实施例提到的水动力载荷为波浪激励力载荷、阻尼系数函数以及附加质量矩阵时，对上述实施例提到的基于所述漂浮式风力机的当前位置，融合所述气动载荷、所述系泊力载荷、所述水动力载荷、所述静态浮力载荷、所述黏性力载荷以及所述静水回复力矩阵，确定所述漂浮式风力机在下一时间步的坐标及运动信息的过程进行介绍，该过程可以包括：
[0136]
s1、确定所述漂浮式风力机在所述预设时长后的下一时间步的位移量、速度和加速度。
[0137]
具体的，可以利用如下微分方程，得到所述漂浮式风力机在所述预设时长后的下一时间步的位移量、速度和加速度：
[0138][0139]
其中，[m]表示基于所述漂浮式风力机的机体质量所得到的质量矩阵， [a]
inf
表示所述附加质量矩阵，k(t)表示以所述预设时长为自变量的阻尼系数函数，[c]表示所述静水回复力矩阵，f
waves
表示所述波浪激励力载荷对应的波浪激励力，f
mooring
表示所述系泊力载荷对应的系泊力，f
drag
表示所述黏性力载荷对应的黏性力，f
wind
表示所述气动载荷对应的风力，f
static
表示所述静态浮力载荷对应的静水浮力，t表示所述预设时间，表示加速度，表示速度，q表示位移量，所述位移量、速度和加速度组成运动信息。
[0140]
s2、基于所述漂浮式风力机的当前位置以及所述预设时长后的下一时间步的位移量，确定所述漂浮式风力机在下一时间步的坐标。
[0141]
具体的，所述位移量可以包含方向和距离，基于所述漂浮式风力机的当前位置以
及需要移动的方向和距离，可以确定目的位置的坐标。
[0142]
本技术提供的实施例中，能够通过漂浮式风力机当前时间步下的波浪激励力载荷、阻尼系数函数、附加质量矩阵、气动载荷、系泊力载荷、静态浮力载荷、黏性力载荷以及静水回复力矩阵，计算得到漂浮式风力机在下一时间步的坐标和运动信息。
[0143]
考虑到对漂浮式风力机流体动力学耦合分析所存储的数据进行处理，本技术的一些实施例中，在上述实施例提到的在当前时间步达到预设时间步之后，还可以包括：汇总所述漂浮式风力机在各时间步下的各载荷、坐标及运动信息，得到对漂浮式风力机在预设时间步下的所有流体动力学耦合分析的所有结果，以供对所述漂浮式风力机进行性能预报。
[0144]
下面对本技术实施例提供实现的漂浮式风力机流体动力学耦合分析装置进行描述，下文描述实现的漂浮式风力机流体动力学耦合分析装置与上文描述实现的漂浮式风力机流体动力学耦合分析方法可相互对应参照。
[0145]
参见图2，图2为本技术实施例公开的一种实现的漂浮式风力机流体动力学耦合分析装置的结构示意图。
[0146]
如图2所示，该装置可以包括：
[0147]
浮式基础分析单元11，用于分析漂浮式风力机的浮式基础，确定并存储所述漂浮式风力机在预设时长的当前时间步下的水动力载荷，以及确定并存储静水回复力矩阵；
[0148]
系泊力载荷确定单元12，用于分析所述漂浮式风力机中缆绳的当前状态信息，确定并存储所述当前时间步下所述漂浮式风力机的系泊力载荷；
[0149]
双载荷确定单元13，用于基于所述漂浮式风力机在所述当前时间步下的机体参数信息和/或所处于水体的水体参数信息，确定并存储所述漂浮式风力机在所述当前时间步下的静态浮力载荷、黏性力载荷；
[0150]
气动参数获取单元14，用于获取所述漂浮式风力机的叶片气动参数，以及所述漂浮式风力机的各气动节点的坐标和运动信息；
[0151]
气动载荷计算单元15，用于基于所述叶片气动参数与所述各气动节点的坐标和运动信息，通过叶素动量理论及其修正算法，计算并存储所述漂浮式风力机在当前时间步下的气动载荷；
[0152]
运动信息确定单元16，用于基于所述漂浮式风力机的当前位置，融合所述气动载荷、所述系泊力载荷、所述水动力载荷、所述静态浮力载荷、所述黏性力载荷以及所述静水回复力矩阵，确定、更新并存储所述漂浮式风力机在下一时间步的坐标及运动信息，并将下一时间步作为当前时间步，返回执行所述浮式基础分析单元11用于分析漂浮式风力机的浮式基础，确定并存储所述漂浮式风力机在预设时长的当前时间步下的水动力载荷的步骤，直到达到预设时间步为止。
[0153]
本技术实施例提供的漂浮式风力机流体动力学耦合分析装置可应用于漂浮式风力机流体动力学耦合分析设备，如终端：手机、电脑等。可选的，图3 示出了漂浮式风力机流体动力学耦合分析的硬件结构框图，参照图3，漂浮式风力机流体动力学耦合分析设备的硬件结构可以包括：至少一个处理器1，至少一个通信接口2，至少一个存储器3和至少一个通信总线4；
[0154]
在本技术实施例中，处理器1、通信接口2、存储器3、通信总线4的数量为至少一个，且处理器1、通信接口2、存储器3通过通信总线4完成相互间的通信；
[0155]
处理器1可能是一个中央处理器cpu，或者是特定集成电路asic (application specific integrated circuit)，或者是被配置成实施本发明实施例的一个或多个集成电路等；
[0156]
存储器3可能包含高速ram存储器，也可能还包括非易失性存储器 (non-volatile memory)等，例如至少一个磁盘存储器；
[0157]
其中，存储器存储有程序，处理器可调用存储器存储的程序，所述程序用于：
[0158]
分析漂浮式风力机的浮式基础，确定并存储所述漂浮式风力机在预设时长的当前时间步下的水动力载荷，以及确定并存储静水回复力矩阵；
[0159]
分析所述漂浮式风力机中缆绳的当前状态信息，确定并存储所述当前时间步下所述漂浮式风力机的系泊力载荷；
[0160]
基于所述漂浮式风力机在所述当前时间步下的机体参数信息和/或所处于水体的水体参数信息，确定并存储所述漂浮式风力机在所述当前时间步下的静态浮力载荷、黏性力载荷；
[0161]
获取所述漂浮式风力机的叶片气动参数，以及所述漂浮式风力机的各气动节点的坐标和运动信息；
[0162]
基于所述叶片气动参数与所述各气动节点的坐标和运动信息，通过叶素动量理论及其修正算法，计算并存储所述漂浮式风力机在当前时间步下的气动载荷；
[0163]
基于所述漂浮式风力机的当前位置，融合所述气动载荷、所述系泊力载荷、所述水动力载荷、所述静态浮力载荷、所述黏性力载荷以及所述静水回复力矩阵，确定、更新并存储所述漂浮式风力机在下一时间步的坐标及运动信息，并将下一时间步作为当前时间步，返回执行分析漂浮式风力机的浮式基础，确定并存储所述漂浮式风力机在预设时长的当前时间步下的水动力载荷的步骤，直到达到预设时间步为止。
[0164]
可选的，所述程序的细化功能和扩展功能可参照上文描述。
[0165]
本技术实施例还提供一种存储介质，该存储介质可存储有适于处理器执行的程序，所述程序用于：
[0166]
分析漂浮式风力机的浮式基础，确定并存储所述漂浮式风力机在预设时长的当前时间步下的水动力载荷，以及确定并存储静水回复力矩阵；
[0167]
分析所述漂浮式风力机中缆绳的当前状态信息，确定并存储所述当前时间步下所述漂浮式风力机的系泊力载荷；
[0168]
基于所述漂浮式风力机在所述当前时间步下的机体参数信息和/或所处于水体的水体参数信息，确定并存储所述漂浮式风力机在所述当前时间步下的静态浮力载荷、黏性力载荷；
[0169]
获取所述漂浮式风力机的叶片气动参数，以及所述漂浮式风力机的各气动节点的坐标和运动信息；
[0170]
基于所述叶片气动参数与所述各气动节点的坐标和运动信息，通过叶素动量理论及其修正算法，计算并存储所述漂浮式风力机在当前时间步下的气动载荷；
[0171]
基于所述漂浮式风力机的当前位置，融合所述气动载荷、所述系泊力载荷、所述水动力载荷、所述静态浮力载荷、所述黏性力载荷以及所述静水回复力矩阵，确定、更新并存储所述漂浮式风力机在下一时间步的坐标及运动信息，并将下一时间步作为当前时间步，
返回执行分析漂浮式风力机的浮式基础，确定并存储所述漂浮式风力机在预设时长的当前时间步下的水动力载荷的步骤，直到达到预设时间步为止。
[0172]
可选的，所述程序的细化功能和扩展功能可参照上文描述。
[0173]
最后，还需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个
……”
限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
[0174]
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间可以根据需要进行组合，且相同相似部分互相参见即可。
[0175]
对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本技术。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本技术的精神或范围的情况下，在其它实施例中实现。因此，本技术将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。