基于代理模型的海上风机支撑结构优化设计方法及系统

摘要

本发明公开基于代理模型的海上风机支撑结构优化设计方法及系统及系统，方法基于Kriging代理模型建立不同容量大小海上风电机组的支撑结构设计参数与优化设计的结构性能响应之间的关系，同时对塔架和单桩基础进行优化设计以找到整体支撑结构质量最轻的优化设计；采用本发明所述方法能够避免海上风机支撑结构优化设计时陷入局部最优解，导致找到的只是塔架或者单桩最轻设计的问题；采用本发明所用的方法和流程，在载荷计算后同时对塔架和基础结构进行校核和优化设计，能找到全局最优设计，使得结构整体性更强，刚度分布更均匀，受力更合理，最终达到整体支撑结构重量最轻，基于代理模型进行优化设计减少需要的计算机模拟结构分析次数，提高优化效率。

说明书

技术领域

本发明属于海上风力发电机组支撑结构设计技术领域，具体涉及基于代理模型的海上风机支撑结构优化设计方法及系统。

背景技术

受益于风电的技术进步和规模扩大，风电机组价格、风电开发投资成本及运行维护成本呈现不断下降趋势。从风电机组价格来看，海上风机支撑结构包括塔架和基础两部分：风机塔架成本约占海上风电项目投资成本8％左右，海上风机基础主要包括单桩、导管架、高桩承台等不同基础形式，一般占海上风电项目投资成本的14％左右，即整体支撑结构成本在总建设成本中占比约22％左右。因此，降低海上风电支撑结构成本能够有效降低海上风电的平准化度电成本。

当前国内海上风电项目在招投标时通常采用分步迭代设计方法，一般由风机厂家给出塔架设计并担保塔架工程量，评标过程中会对塔架重量进行排名和评分；标后详细设计阶段时风机厂家和设计院依次分别对塔架和基础进行优化设计。在此流程下，风机厂家会尽可能给出塔架最轻的局部最优设计方案，而塔架最轻的设计往往不是整体支撑结构最轻的全局最优设计方案。目前已有的研究往往只对整体支撑结构的一部分(塔架或者基础)进行优化设计，因此得到的结果具有一定局限性。

在设计海上风机支撑结构时都包含了三部分：载荷计算、塔架设计和基础设计。

1)载荷计算

海上风电支撑结构受到风、浪、流等多种环境载荷的联合作用。行业内大部分风机厂家采用GH-Bladed进行一体化建模和载荷计算。

一体化建模包含环境条件输入和整体支撑结构模型搭建两方面。其中，环境条件包括风资源参数、海洋水文参数、工程地质参数及其他特殊工况(海冰、地震、台风等)；整体支撑结构模型包括机头、塔架、泥面以上结构和基础(也可统称为基础结构)。

载荷计算中考虑了风浪异向作用，根据IEC规范需要考虑正常发电、紧急停机、开机、正常停机、空转、维修等多种工况，根据风浪联合分布可分成多达20000多个工况。

2)塔架设计

塔架设计中，需要进行塔架主体及局部结构的极限强度、屈曲强度和疲劳强度校核。极限强度校核包括塔架筒体、塔架法兰、门洞及海缆孔、锚栓笼等局部结构的校核；屈曲强度校核包括塔架筒体和门洞海缆孔等结构的校核；疲劳强度校核包括塔架筒体焊缝、法兰连接螺栓、门框及海缆孔、顶法兰、锚栓笼等结构的校核。

3)基础设计

在基础结构主体设计中，主要包含极限海况下的强度承载力分析、正常服役工况分析、船撞分析、地震工况分析等。荷载组合中考虑可能出现的最不利水位下的波浪、海流与风机运行载荷的极端组合作用。疲劳强度分析利用S-N曲线与Miner线性累计损伤理论进行疲劳计算。分别计算各个管节点在疲劳荷载作用下累积损伤程度，利用累积损伤程度评估结构的抗疲劳设计安全性。

当前国内风电行业大都采用的是分步迭代设计方法。图1给出了海上单桩基础整体支撑结构的示意图。如图1所示，整体支撑结构以设计交界面为分界，交界面以上为塔架，交界面以下为基础结构。图2给出了分步迭代设计方法的流程。首先由设计院提供项目的环境输入；风机厂家根据环境输入给出塔架和基础的初始构型并进行整体建模和载荷计算，在得到最优的塔架后把设计交界面处的载荷、塔架构型和频率要求提给设计院；接着设计院在给定载荷和塔架构型的前提下对基础结构进行校核和优化设计，并满足风机厂家给出的频率要求；最后风机厂家在得到优化的基础结构后判断是否收敛，如果满足则迭代结束，如不满足则重新建模并进行载荷计算。在这里收敛准则包含两类：一类是根据规范对塔架和基础校核的设计准则；另一类是本轮和上轮得到的优化设计的质量和频率差别是否在1％以内。

在这里，需要说明的是，在当前国内采用海上风机支撑结构设计方法时，在确定了初步构型以后(塔架和单桩基础的直径)一般需要2-4轮迭代才能收敛，每一轮都需要进行载荷计算、塔架和基础的设计优化。如果要进一步优化塔架和单桩的直径以找到整体支撑结构质量最轻的设计，将十分耗费计算时间，以致于影响项目进度。因此，在实际工程项目中为了尽快提供塔架和单桩基础的施工图，往往没有足够的时间用于进行优化，并且在此过程中，塔架和基础的设计和优化是依次先后进行的，是两个独立的设计域，目标都是找到各自设计域中的最优设计。因此在国内实际的项目中，最终得到的设计往往是塔架最轻的局部最优设计，而不是整体支撑结构最轻的全局最优设计。

如果我们要采用基于灵敏度的优化算法来求解海上风机支撑结构的最优设计，只能通过有限差分法求得结构性能对设计变量的灵敏度，工作量十分大，无法满足实际工期的要求。

发明内容

为了解决现有技术中存在的问题，本发明提供基于代理模型的海上风机支撑结构优化设计方法及系统，同时对塔架和基础结构进行优化设计，站在寻找全局最优设计的角度对海上风机支撑结构进行设计，降低海上风电平准化度电成本以及设计周期。

为了实现上述目的，本发明采用的技术方案是：基于代理模型的海上风机支撑结构优化设计方法及系统，包括以下步骤：

基于确定风机支撑结构的设计变量，采用抽样方法在设计空间内生成初始样本集合；

对所述初始样本集合中的每一个样本点进行数值模拟，得到相应的结构性能响应结果；

选择代理模型，基于所述初始样本集合和结构性能响应结果，采用所对应的拟合或插值方法建立根据初始样本集合和结构性能响应结果之间的映射关系；

使用代理模型进行结构优化设计，得到塔架和基础结构的初步设计结果，在进行结构优化设计时，在载荷计算后同时对塔架和基础结构进行校核和优化设计；

对所得初步设计结果采用数值仿真进行模拟，确认其可行性和最优性，并校验是否满足收敛准则；如果不满足收敛准则，则选择合适的加点准则添加样本点，更新代理模型，采用更新后的代理模型进行优化设计，直至满足收敛准则，则优化迭代结束。

所述设计变量为塔底直径和单桩直径；结构性能响应为塔架质量、单桩质量以及整体支撑结构频率。

采用中心复合试验方法生成初始样本集合。

结构性能响应为塔架质量、单桩质量以及整体支撑结构频率。

采用Kriging代理模型。

所述收敛准则为：

Δ

Δ

其中，Δ

采用代理模型优化过程中的优化列式为：

find:塔底和单桩基础的直径、壁厚、焊缝高度

minimum:整体支撑结构质量

subject to:

SRF

SRF

UC

Damage

单桩基础的设计抗压承载力≤单桩基础的允许抗压承载力 ⑤

单桩基础的设计变形值≤单桩基础的许用变形值 ⑥

其中：SRF

①塔架最轻的设计不是整体支撑结构最轻的设计；

②整体支撑结构最轻的设计所对应的塔架直径较大：3-5MW级别机组为6.0m，6-8MW机组为7.0m，且对应的单桩直径都稍比塔底直径大0.5m-1.0m范围内，即3-5MW级别机组单桩直径为6.5-7.0m，6-8MW机组单桩直径为7.5-8.0m。

本发明的另一个方面是提供基于代理模型的海上风机支撑结构优化设计系统，包括初始样本集合生成模块、结构性能响应获取模块、映射关系构建模块、优化模块、迭代校验模块；

初始样本集合生成模块基于风机支撑结构的设计变量，采用抽样方法在设计空间内生成初始样本集合；

结构性能响应获取模块用于对所述初始样本集合中的每一个样本点进行数值模拟，得到相应的结构性能响应结果；

映射关系构建模块用于选择代理模型，基于所述初始样本集合和结构性能响应结果，采用所对应的拟合或插值方法建立根据初始样本集合和结构性能响应结果之间的映射关系；

优化模块使用代理模型进行结构优化设计，得到塔架和基础结构的初步设计结果，在进行结构优化设计时，在载荷计算后同时对塔架和基础结构进行校核和优化设计；

迭代校验模块对所得初步设计结果采用数值仿真进行模拟，确认其可行性和最优性，并校验是否满足收敛准则；如果不满足收敛准则，则选择合适的加点准则添加样本点，更新代理模型，采用更新后的代理模型进行优化设计，直至满足收敛准则，则优化迭代结束。

本发明还提供一种计算机设备，包括一个或多个处理器以及存储器，存储器用于存储计算机可执行程序，处理器从存储器中读取部分或全部所述计算机可执行程序并执行，处理器执行部分或全部计算可执行程序时能实现本发明所述基于代理模型的海上风机支撑结构优化设计方法。

与现有技术相比，本发明至少具有以下有益效果：

采用本发明所述方法能够避免海上风机支撑结构优化设计时陷入局部最优解，导致找到的只是塔架或者单桩最轻设计的问题；采用本发明所用的方法和流程，在载荷计算后同时对塔架和基础结构进行校核和优化设计，能找到全局最优设计，使得结构整体性更强，刚度分布更均匀，受力更合理，最终达到整体支撑结构重量最轻，降低海上风机支撑结构造价的目的；基于代理模型进行优化设计减少需要的计算机模拟结构分析次数，提高优化效率。

附图说明

下面结合附图对本发明的示例性实施例进行详细描述，本发明的以上和其它特点及优点将变得更加清楚，附图中：

图1是海上支撑结构示意图。

图2是海上风机支撑结构分步迭代设计方法流程。

图3是代理模型优化算法流程。

图4是海上风机支撑结构整体优化设计方法流程。

图5是3-5MW级海上风电机组塔底直径及单桩直径初始样本点。

图6是6-8MW级海上风电机组塔底直径及单桩直径初始样本点。

图7是3-5MW级海上风电机组塔架质量代理模型。

图8是3-5MW级海上风电机组整体支撑结构质量代理模型。

图9是6-8MW级海上风电机组塔架质量代理模型。

图10是6-8MW级海上风电机组整体支撑结构质量代理模型。

具体实施方式

下面结合图3-图10和具体实施例对本发明作进一步详细说明。

参考图3中，代理模型优化算法的一般流程图，包括以下步骤：

步骤1，通过工程项目经验获取海上支撑结构设计特征并筛选设计变量，在本发明中设计变量为为塔底直径和单桩直径；

步骤2，采用抽样方法在设计空间内生成初始样本集合，本发明采用中心复合试验方法；

步骤3，对所述初始样本集合中的每一个样本点进行数值模拟分析，得到相应的结构性能响应，本发明中结构性能响应为塔架质量、单桩质量以及整体支撑结构频率。

步骤4，选择代理模型(建议但不限于Kriging代理模型)，采用所对应的合适的拟合或插值方法根据步骤2和步骤3中的输入和输出，建立它们之间的映射关系。

步骤5，使用代理模型进行结构优化设计，得到塔架和基础结构的初步设计结果，在进行结构优化设计时，在载荷计算后同时对塔架和基础结构进行校核和优化设计；

步骤6，对步骤5所得初步设计结果采用数值仿真进行模拟分析，确认其可行性和最优性，并校验是否满足收敛准则；如果不满足收敛准则，则选择合适的加点准则添加样本点，更新代理模型，返回步骤5。如果满足收敛准则，则优化迭代结束；这里的收敛准则指：本轮和上轮得到的优化设计的响应值及塔架质量、单桩质量以及整体支撑结构频率等物理量差别是否在1％以内。

在如上所述的步骤2对样本点进行数值模拟的方法，本发明提出的一种基于代理模型的海上风机支撑结构整体优化设计方法；参考图4所示支撑结构整体优化设计方法的流程，海上风电支撑结构整体优化设计方法与分步迭代法设计方法最大的不同是在载荷计算后同时对塔架和基础结构进行校核和优化设计，在整个设计域中寻找整体支撑结构最轻的全局最优设计。风机厂家和设计院在得到载荷后，同时对塔架和基础结构进行优化设计，站在寻找全局最优设计的角度进行设计。当前海上支撑结构的材料的采购费基本只考虑质量因素，通过支撑结构整体优化设计方法可以降低整体支撑结构的质量，因此达到降低海上风电平准化度电成本的目的。优化列式如下：

find:塔底和单桩基础的直径、壁厚、焊缝高度

minimum:整体支撑结构质量

subject to:

SRF

SRF

UC

Damage

单桩基础的设计抗压承载力≤单桩基础的允许抗压承载力 ⑤

单桩基础的设计变形值≤单桩基础的许用变形值 ⑥

其中：SRF

通过本发明提出的基于代理模型的海上风机支撑结构优化设计方法，可以得到的最优支撑结构：

①塔架最轻的设计不是整体支撑结构最轻的设计；

②整体支撑结构最轻的设计所对应的塔架直径较大(3-5MW级别机组为6.0m，6-8MW机组为7.0m)，且对应的单桩直径都稍比塔底直径大0.5m-1.0m范围内，即3-5MW级别机组为6.5-7.0m，6-8MW机组为7.5-8.0m。

以3-5MW和6-8MW级海上风电机组支撑结构优化设计为例：

步骤1，选择塔底直径和单桩直径为设计变量。

步骤2，选用中心复合试验方法生成初始样本点(如图5是3-5MW级海上风电机组塔底直径及单桩直径初始样本点；图6是6-8MW级海上风电机组塔底直径及单桩直径初始样本点)。

步骤3，采用海上风机支撑结构整体优化得到相应的结构性能响应(塔架质量、单桩质量、整体支撑结构频率等)。

步骤4，选择Kriging代理模型，采用所对应的合适的拟合或插值方法根据第2)，3)步中的输入和输出，建立它们之间的映射关系的代理模型(如图7和图8是3-5MW级海上风电机组塔架质量代理模型和整体支撑结构质量代理模型；图9和图10是6-8MW级海上风电机组塔底直径及单桩直径初始样本点)。

步骤5，使用代理模型代替原模型进行结构优化设计。

步骤6，对步骤5求得的优化设计采用原模型进行校核，确认其可行性和最优性，并校验收敛准则。如果不满足，则选择合适的加点准则添加样本点，更新代理模型，返回步骤5。如果满足，则优化迭代结束。

这里的收敛准则指：本轮和上轮得到的优化设计的响应值(塔架质量、单桩质量、整体支撑结构频率等物理量)差别是否在1％以内。

优化列式如下：

find:塔底直径，单桩直径

minimum:整体支撑结构质量

subject to:

SRF

SRF

UC

Damage

单桩基础的设计抗压承载力≤单桩基础的允许抗压承载力。

单桩基础的设计变形值≤单桩基础的许用变形值。

表1 3-5MW海上风机支撑结构计算结果

表2 6-8MW海上风机支撑结构计算结果

表1和表2分别给出了基于代理模型优化算法给出的3-5MW和6-8MW海上风机支撑结构计算结果，其中，对于3-5MW机组，塔架最轻的设计为D塔架＝5.5m，D单桩＝6.0m；整体支撑结构最轻的设计为D塔架＝6.0m，D桩＝6.5m。

对于6-8MW机组，塔架最轻的设计为D塔架＝6.5m，D单桩＝8.5m；整体支撑结构最轻的设计为D塔架＝7.0m，D桩＝7.5m。

本发明还提供一种计算机设备，包括但不限于一个或多个处理器以及存储器，存储器用于存储计算机可执行程序，处理器从存储器中读取部分或全部所述计算机可执行程序并执行，处理器执行部分或全部计算可执行程序时能实现本发明所述基于代理模型的海上风机支撑结构优化设计方法的部分步骤或所有步骤。

所述识别以太坊上庞氏骗局的的设备可以是笔记本电脑、平板电脑、桌面型计算机或工作站。

处理器可以是中央处理器(CPU)、数字信号处理器(DSP)、专用集成电路(ASIC)或现成可编程门阵列(FPGA)。

对于本发明所述存储器，可以是笔记本电脑、平板电脑、桌面型计算机、手机或工作站的内部存储单元，如内存、硬盘；也可以采用外部存储单元，如移动硬盘、闪存卡。