CN115238409A - 一种风机多源模型耦合协同仿真方法、平台、介质及设备 - Google Patents

一种风机多源模型耦合协同仿真方法、平台、介质及设备 Download PDF

Info

Publication number
CN115238409A
CN115238409A CN202210787705.7A CN202210787705A CN115238409A CN 115238409 A CN115238409 A CN 115238409A CN 202210787705 A CN202210787705 A CN 202210787705A CN 115238409 A CN115238409 A CN 115238409A
Authority
CN
China
Prior art keywords
model
fan
establishing
file
function
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Withdrawn
Application number
CN202210787705.7A
Other languages
English (en)
Inventor
胡广
黄蓉
黄国燕
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
MingYang Smart Energy Group Co Ltd
Original Assignee
MingYang Smart Energy Group Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by MingYang Smart Energy Group Co Ltd filed Critical MingYang Smart Energy Group Co Ltd
Priority to CN202210787705.7A priority Critical patent/CN115238409A/zh
Publication of CN115238409A publication Critical patent/CN115238409A/zh
Priority to CN202310725246.4A priority patent/CN116910929A/zh
Withdrawn legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational design
    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/02—CAD in a network environment, e.g. collaborative CAD or distributed simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/06—Wind turbines or wind farms
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Physics (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Computational Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Automation & Control Theory (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种风机多源模型耦合协同仿真方法、平台、介质及设备，该方法首先通过在SolidWorks软件中建立风机的三维实体模型、在MATLAB软件中建立用于解析风能的转换过程的风机数学模型、在Bladed软件中建立风机运行参数模型以及在Simulink软件中搭建风机控制器模型后，再通过多源模型资源接口方法对三维实体模型、风机数学模型、风机运行参数模型和风机控制器模型进行封装，设计通用模型资源实体组件及组件描述文件，从而建立数据共享机制，最终进行各模型的耦合仿真并将结果体现在风机的各部件实体上；本发明可实时呈现风机各部件的真实运动状态，大幅度降低计算资源，节约时间成本，快速分析风机载荷变化情况与外部激励的关系，为风机控制器的调试提供技术支持。

Description

一种风机多源模型耦合协同仿真方法、平台、介质及设备
技术领域
本发明涉及风机仿真建模的技术领域，尤其是指一种风机多源模型耦合协同仿真方法、平台、介质及设备。
背景技术
风电平价时代来临，风机大型化是大势所趋，需兼顾风机控制策略、风能波动、电网扰动等对机组各部件的相互作用及动态特性。风机本身的气动模型与大部分其他模型直接和间接耦合，譬如，非扭转机械结构的振动形成一个相对风速并与自然风速叠加，同时气动模型输出的轴向推力、机舱点头力矩等构成非扭转机械模型的振动力源，如此形成两模型之间的气–弹耦合，气动模型输出的转矩和发电机电磁转矩共同作用于发电机轴系，风工况分布不均导致的气动力矩波动，以及电网故障激发的电磁暂态均可能引起发电机轴系发生扭振，由此产生的机电耦合等情况；鉴于风机涉及到空气动力学、多体动力学、液压、电气等多个专业领域，复杂程度高，为在产品开发前期全方位发现其潜在不足，对风机进行多物理联合仿真的需求日益明显。
发明内容
本发明的第一目的在于为解决现有技术中的不足，提供了一种风机多源模型耦合协同仿真方法，将SolidWorks、Matlab、Simulink和Bladed模型互通，建立数据共享机制，进行各模型的耦合仿真并将结果体现在风机各部件实体上。
本发明的第二目的在于提供一种风机多源模型耦合协同仿真平台。
本发明的第三目的在于提供一种非暂时性计算机可读介质。
本发明的第四目的在于提供一种计算设备。
本发明的第一目的通过下述技术方案实现：一种风机多源模型耦合协同仿真方法，该方法首先通过在SolidWorks软件中建立风机的三维实体模型、在MATLAB软件中建立用于解析风能的转换过程的风机数学模型、在Bladed软件中建立风机运行参数模型以及在Simulink软件中搭建风机控制器模型后，再通过多源模型资源接口方法对三维实体模型、风机数学模型、风机运行参数模型和风机控制器模型进行封装，设计通用模型资源实体组件及组件描述文件，从而建立数据共享机制，最终进行各模型的耦合仿真并将结果体现在风机的各部件实体上。
进一步，该方法包括以下步骤：
S1、建立三维实体模型、风机数学模型、风机运行参数模型和风机控制器模型；
S2、采用Matlab/C++混合编程方式，以Matlab Engine作为控制接口，实现对三维实体模型、风机数学模型、风机运行参数模型和风机控制器模型的控制；
S3、利用Matlab库函数实现与Matlab模型间的信息交互，设计通用模型资源实体组件.DLL及组件描述文件.xml，调用风机运行参数模型中的风机运行参数，以共享内存的方式实现通用模型资源实体组件.DLL与风机控制器模型间的信息交互，风机运行参数模型与外部软件的数据交互采用API方式，Simulink得到风机的风力荷载后加载至三维实体模型上，使风机的运动情况实时可视化。
进一步，在步骤S1中，具体执行以下操作：
在SolidWorks软件中建立风机的三维实体模型，在SolidWorks软件的控制动力学模块设置风机各部件的材料属性，在各部件间施加固定副、旋转副，创建接触力约束，在叶片上创建风力作为驱动；
在MATLAB软件中依次建立风机的气动、传动链、塔架和发电机的数学模型，解析风能的转换过程；气动数学模型不考虑风剪切在叶片上造成载荷不均衡现象，叶片上的气动力、气动力矩相等的基础上，调用预设的翼型气动数据库，沿叶片展向进行离散线积分，求取整支叶片产生的气动力，其中所述气动翼型数据库为翼型在不同攻角下，从切入至切出这一风速段内的推力系数和扭矩系数；传动链模型简化为叶轮—传动链和传动链—发电机两振系，均为质量—弹簧—阻尼系统，建立运动微分方程：
Figure BDA0003732159420000031
其中，k为弹簧刚度，c为阻尼器阻尼，Q为外力矩，θ为角位移，是相对值，J为惯性质量矩；
塔架模型视为质量—弹簧—阻尼系统，其变量关系满足：
Figure BDA0003732159420000032
其中，k为弹簧刚度，c为阻尼器阻尼，F为外力，x为位移，是相对于地面的相对值，M为惯性质量矩；
发电机模型为发电机和变频器组成的变速驱动装置，用一阶迟滞模型反应该变频器的转矩响应：
Figure BDA0003732159420000033
其中，Qd为所请求的扭矩，Qg为气隙扭矩，τe为一阶迟滞的时间常数，S表示该传递函数是原函数经过拉普拉斯变换后的表达式；
在Simulink软件中使用Logic and Bit Operations模块、Lookup Tables模块、Math Operations模块以及Port&Subsystems模块，实现对机组的启机、并网、停机、偏航对风、机组监测与保护、故障诊断与智能维护的控制逻辑，最终搭建出控制器模型。
进一步，在步骤S3中，具体执行以下操作：
1)设定共享内存区域，Bladed进程和Simulink进程能够对该共享内存区域的地址同时访问，能同时对共享内存区域中的数据进行读写，首先利用CreateFileMapping()函数对Bladed计算结果文件或Simulink指令文件进行存取操作，在该函数中对Bladed计算结果文件或Simulink指令文件进行操作等同于直接对共享内存数据进行操作；在成功创建共享文件后，CreateFileMapping()函数产生一个file-mapping核心对象，file-mapping核心对象记录了内存数据与共享文件数据的映射关系，再利用MapViewFile()函数，产生一个指向可用内存的指针，指向被映射的视图开头；
2)使用共享内存，在完成共享内存的定义与创建并建立了映射关系后，利用file-mapping核心对象使得其它需要访问共享内存的进程都能通过名称找到共享内存，并将其映射到各进程的地址空间中，构建成视图；再利用OpenFileMapping()函数，向访问共享内存的进程返回一个句柄，代表一个file-mapping核心对象，同时进程通过调用OpenFileMapping()函数获得一个指向共享内存的指针；
3)清理共享内存，在完成对共享内存的读写操作后，调用UnmapViewOfFile()函数，将原本由调用MapViewOfView()函数而获得的指针交出，再调用CloseHandle()函数，把file-mapping核心对象的句柄交出；
4)利用MATLAB软件对描述有共享内存接口调用信息的程序进行编译，生成MEX文件，利用S-Function将编译通过的MEX文件转换为被风机控制器模型应用的控件资源，使用mdlInitializeSizes()函数确定风机控制器模型的输入输出口的个数，在MEX文件中使用mdlInitializeSampleTimes()函数对S-Function函数的采样周期作出定义，最终采用内存共享技术来完成Bladed进程和Simulink进程之间的通信，确保数据的传输可靠。
本发明的第二目的通过下述技术方案实现：一种风机多源模型耦合协同仿真平台，包括：
基于Matlab的风机数学建模模组，用于建立风机的气动、传动链、塔架和发电机的数学模型；
基于Bladed的风机运行参数建模模组，用于对风机建立运行参数模型；
基于SolidWorks的风机可视化模组，用于对风机建立三维实体模型；
基于Simulink的风机控制器建模模组，用于对机组的启机、并网、停机、偏航对风、机组监测与保护、故障诊断与智能维护的控制逻辑；
通用模型资源实体模组，用于建立数据共享机制，最终进行各模型的耦合仿真。
本发明的第三目的通过下述技术方案实现：一种存储有指令的非暂时性计算机可读介质，当所述指令由处理器执行时，执行根据上述的风机多源模型耦合协同仿真方法的步骤。
本发明的第四目的通过下述技术方案实现：一种计算设备，包括处理器以及用于存储处理器可执行程序的存储器，所述处理器执行存储器存储的程序时，实现上述的风机多源模型耦合协同仿真方法。
本发明与现有技术相比，具有如下优点与有益效果：
1、通过该联合仿真平台可实时呈现风机各部件的真实运动状态，丰富了其它风机仿真软件的后处理功能，大幅度降低计算资源，节约时间成本，可快速分析风机载荷变化情况与外部激励两者的关系，为项目早期风机控制器的调试提供技术支持；
2、本发明提出Simulink与Bladed通过共享内存实现联合仿真方法的优点是无需对Matlab软件本身进行二次开发，在软件现有功能基础上提出一个通用的联合仿真技术，实现Simulink与Bladed的联合仿真数据交互和联合仿真过程中的时序控制。
附图说明
图1为各模型的数据交互示意图。
图2为SolidWorks风电机组拓扑图。
图3为简易风机本体模型框架图。
图4为Simulink控制策略模型逻辑流程图。
具体实施方式
下面结合具体实施例对本发明作进一步说明。
实施例1
参见图1至图4所示，为本实施例所提供的风机多源模型耦合协同仿真方法，该方法首先通过在SolidWorks软件中建立风机的三维实体模型、在MATLAB软件中建立用于解析风能的转换过程的风机数学模型、在Bladed软件中建立风机运行参数模型以及在Simulink软件中搭建风机控制器模型后，再通过多源模型资源接口方法对三维实体模型、风机数学模型、风机运行参数模型和风机控制器模型进行封装，设计通用模型资源实体组件及组件描述文件，从而建立数据共享机制，最终进行各模型的耦合仿真并将结果体现在风机的各部件实体上，包括以下步骤：
S1、建立三维实体模型、风机数学模型、风机运行参数模型和风机控制器模型，具体执行以下操作：
在SolidWorks软件中建立风机的三维实体模型，在SolidWorks软件的控制动力学模块设置风机各部件的材料属性，在各部件间施加固定副、旋转副，创建接触力约束，在叶片上创建风力作为驱动；
参见图2所示，风电机组传动链主要由叶片1、叶片2、叶片3、轮毂、主轴承、主轴、齿轮箱、刹车片、高速轴、联轴器、发电机、弹性支撑等部件组成，为确定传动链的动力学行为，需设定：1)各轴承、部件及弹性支撑等的定位、定形尺寸；2)轮毂等柔性部件的力学性能参数；3)各部件的质量、质心及转动惯量；4)联轴器、弹性支撑件的刚度、阻尼系数；5)轴承6个方向的刚度及阻尼系数；6)发电机的功率外特性曲线；7)对主轴-行星架、行星齿轮、一号齿轮-太阳轮、三号齿轮-四号齿轮，以及二号齿轮-输出轴创建连杆，以每个活动连杆的旋转中心创建旋转副，按照传动比设置相应的齿轮副，设置齿轮接触模式为实体接触；8)将SolidWorks的运动模型通过协同仿真生成一个基于Simulink的各类数据处理模组，将此模块添加到Simulink中。
在MATLAB软件中依次建立风机的气动、传动链、塔架和发电机的数学模型，解析风能的转换过程；1)确定单叶片的推力系数CT和单叶片扭矩系数CQ，计算出单叶片推力T为：
Figure BDA0003732159420000071
式中ρ为空气密度，V0为风速，A为叶轮扫风面积。
单叶片扭矩M为：
Figure BDA0003732159420000072
机舱处的力系点头弯矩M_tilt，偏航弯M_yaw，前后推力F_faf为：
M_tilt＝2/3*R*(T1*sin(Azimuth)+T2*sin(Azimuth+2*pi/3)+T3*sin(Azimuth+4*pi/3)
M_yaw＝2/3*R*(T1*cos(Azimuth)+T2*cos(Azimuth+2*pi/3)+T3*cos(Azimuth+4*pi/3)
F_faf＝T1+T2+T3
式中R为叶轮半径，Azimuth为方位角，T1、T2、T3分别为叶片1、2、3的推力。
叶轮承受扭矩Ta，压力F_sd为：
Ta＝M1+M2+M3
F_sd＝-2/R*(M1*sin(Azimuth)+M2*sin(Azimuth+2*pi/3)+M3*sin(Azimuth+4*pi/3)
式中M1、M2、M3分别为叶片1、2、3的扭矩。
2)传动链模型由叶轮—传动链和传动链—发电机两振系组成。
叶轮—传动链振系：
Ta＝TRotorIner-Tsh
Tsh＝c×(ω1-ω2)+k×(θ1-θ2)
传动链—发电机振系：
Tgen＝TGenreatorIner-Tsh&Loss
Tsh&Loss＝Tsh×Efficieny/GearboxRatio
上式中Ta为气动转矩，Tsh为传动力矩，GearboxRatio为齿轮箱传动比，Efficiency为传递效率，ω1为叶轮端角速度，ω2为发电机端转换至叶轮端的角速度，θ1为叶轮端角位移，θ2为发电机端转换至叶轮端的角位移。Tgen为发电机电磁转矩，TRotorIner、TGeneratorIner分别为叶轮和发电机的惯性力矩。
由TRotorIner和TGeneratorIner即可推导出叶轮的角位移θ_Rotor、角速度ω_Rotor，发电机的角位移θ_Gener、角速度ω_Gener以及风轮侧功率MecPower和发电机电磁功率ElePower，如下关系式：
Figure BDA0003732159420000081
Figure BDA0003732159420000082
Figure BDA0003732159420000083
Figure BDA0003732159420000084
MecPower＝Ta×ω_Rotor
ElePower＝Tgen×ω_Gener
3)将塔架模型视为质量—弹簧—阻尼系统，其运动规律满足：
前后方向：
Figure BDA0003732159420000091
左右方向：
Figure BDA0003732159420000092
式中k为弹簧刚度，c为阻尼器阻尼，F_faf为前后推力，F_sd为侧向推力，x为位移，是对大地的相对值，M为惯性质量矩，气动力臂为2H/3，H为塔架高度。
4)发电机模型为发电机和变频器组成的变速驱动装置，由一阶迟滞模型反应该响应，如下：
Figure BDA0003732159420000093
式中，Qd为所请求的扭矩，Qg为气隙扭矩，τe为一阶迟滞的时间常数，S表示该传递函数是原函数经过拉普拉斯变换后的表达式；
在Simulink软件中使用Logic and Bit Operations模块、Lookup Tables模块、Math Operations模块以及Port&Subsystems模块，实现对机组的启机、并网、停机、偏航对风、机组监测与保护、故障诊断与智能维护的控制逻辑，最终搭建出控制器模型。
S2、采用Matlab/C++混合编程方式，以Matlab Engine作为控制接口，实现对三维实体模型、风机数学模型、风机运行参数模型和风机控制器模型的控制；
S3、利用Matlab库函数实现与Matlab模型间的信息交互，设计通用模型资源实体组件.DLL及组件描述文件.xml，调用风机运行参数模型中的风机运行参数，以共享内存的方式实现通用模型资源实体组件.DLL与风机控制器模型间的信息交互，风机运行参数模型与外部软件的数据交互采用API方式，Simulink得到风机的风力荷载后加载至三维实体模型上，使风机的运动情况实时可视化，具体执行以下操作：
1)设定共享内存区域，Bladed进程和Simulink进程能够对该共享内存区域的地址同时访问，能同时对共享内存区域中的数据进行读写，首先利用CreateFileMapping()函数对Bladed计算结果文件或Simulink指令文件进行存取操作，在该函数中对Bladed计算结果文件或Simulink指令文件进行操作等同于直接对共享内存数据进行操作；在成功创建共享文件后，CreateFileMapping()函数产生一个file-mapping核心对象，file-mapping核心对象记录了内存数据与共享文件数据的映射关系，再利用MapViewFile()函数，产生一个指向可用内存的指针，指向被映射的视图开头；其中，Bladed计算结果文件包括风电机组的主轴、发电机转速，功率、风速、气动信息，Simulink指令文件包括Simulink搭建的控制器模型输出的变桨、转矩、偏航、解缆指令；
2)使用共享内存，在完成共享内存的定义与创建并建立了映射关系后，利用file-mapping核心对象使得其它需要访问共享内存的进程都能通过名称找到共享内存，并将其映射到各进程的地址空间中，构建成视图；再利用OpenFileMapping()函数，向访问共享内存的进程返回一个句柄，代表一个file-mapping核心对象，同时进程通过调用OpenFileMapping()函数获得一个指向共享内存的指针；
3)清理共享内存，在完成对共享内存的读写操作后，调用UnmapViewOfFile()函数，将原本由调用MapViewOfView()函数而获得的指针交出，再调用CloseHandle()函数，把file-mapping核心对象的句柄交出；
4)利用MATLAB软件对描述有共享内存接口调用信息的程序进行编译，生成MEX文件，利用S-Function将编译通过的MEX文件转换为可被风机控制器模型应用的控件资源，使用mdlInitializeSizes()函数确定风机控制器模型的输入输出口的个数，在MEX文件中使用mdlInitializeSampleTimes()函数对S-Function函数的采样周期作出定义，最终采用内存共享技术来完成Bladed进程和Simulink进程之间的通信，确保数据的传输可靠。
实施例2
本实施例公开了风机多源模型耦合协同仿真平台，包括：
基于Matlab的风机数学建模模组，用于建立风机的气动、传动链、塔架和发电机的数学模型；
基于Bladed的风机运行参数建模模组，用于对风机建立运行参数模型；
基于SolidWorks的风机可视化模组，用于对风机建立三维实体模型；
基于Simulink的风机控制器建模模组，用于对机组的启机、并网、停机、偏航对风、机组监测与保护、故障诊断与智能维护的控制逻辑；
通用模型资源实体模组，用于建立数据共享机制，最终进行各模型的耦合仿真。
实施例3
本实施例公开了一种存储有指令的非暂时性计算机可读介质，当所述指令由处理器执行时，执行根据实施例1所述的风机多源模型耦合协同仿真方法的步骤。
本实施例中的非暂时性计算机可读介质可以是磁盘、光盘、计算机存储器、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、U盘、移动硬盘等介质。
实施例4
本实施例公开了一种计算设备，包括处理器以及用于存储处理器可执行程序的存储器，所述处理器执行存储器存储的程序时，实现实施例1所述的风机多源模型耦合协同仿真方法。
本实施例中所述的计算设备可以是台式电脑、笔记本电脑、智能手机、PDA手持终端、平板电脑、可编程逻辑控制器(PLC，Programmable Logic Controller)、或其它具有处理器功能的终端设备。
以上所述之实施例只为本发明之较佳实施例，并非以此限制本发明的实施范围，故凡依本发明之形状、原理所作的变化，均应涵盖在本发明的保护范围内。

Claims (7)

1.一种风机多源模型耦合协同仿真方法，其特征在于：该方法首先通过在SolidWorks软件中建立风机的三维实体模型、在MATLAB软件中建立用于解析风能的转换过程的风机数学模型、在Bladed软件中建立风机运行参数模型以及在Simulink软件中搭建风机控制器模型后，再通过多源模型资源接口方法对三维实体模型、风机数学模型、风机运行参数模型和风机控制器模型进行封装，设计通用模型资源实体组件及组件描述文件，从而建立数据共享机制，最终进行各模型的耦合仿真并将结果体现在风机的各部件实体上。
2.根据权利要求1所述的一种风机多源模型耦合协同仿真方法，其特征在于，包括以下步骤：
S1、建立三维实体模型、风机数学模型、风机运行参数模型和风机控制器模型；
S2、采用Matlab/C++混合编程方式，以Matlab Engine作为控制接口，实现对三维实体模型、风机数学模型、风机运行参数模型和风机控制器模型的控制；
S3、利用Matlab库函数实现与Matlab模型间的信息交互，设计通用模型资源实体组件.DLL及组件描述文件.xml，调用风机运行参数模型中的风机运行参数，以共享内存的方式实现通用模型资源实体组件.DLL与风机控制器模型间的信息交互，风机运行参数模型与外部软件的数据交互采用API方式，Simulink得到风机的风力荷载后加载至三维实体模型上，使风机的运动情况实时可视化。
3.根据权利要求2所述的一种风机多源模型耦合协同仿真方法，其特征在于，在步骤S1中，具体执行以下操作：
在SolidWorks软件中建立风机的三维实体模型，在SolidWorks软件的控制动力学模块设置风机各部件的材料属性，在各部件间施加固定副、旋转副，创建接触力约束，在叶片上创建风力作为驱动；
在MATLAB软件中依次建立风机的气动、传动链、塔架和发电机的数学模型，解析风能的转换过程；气动数学模型不考虑风剪切在叶片上造成载荷不均衡现象，叶片上的气动力、气动力矩相等的基础上，调用预设的翼型气动数据库，沿叶片展向进行离散线积分，求取整支叶片产生的气动力，其中所述气动翼型数据库为翼型在不同攻角下，从切入至切出这一风速段内的推力系数和扭矩系数；传动链模型简化为叶轮—传动链和传动链—发电机两振系，均为质量—弹簧—阻尼系统，建立运动微分方程：
Figure FDA0003732159410000021
其中，k为弹簧刚度，c为阻尼器阻尼，Q为外力矩，θ为角位移，是相对值，J为惯性质量矩；
塔架模型视为质量—弹簧—阻尼系统，其变量关系满足：
Figure FDA0003732159410000022
其中，k为弹簧刚度，c为阻尼器阻尼，F为外力，x为位移，是相对于地面的相对值，M为惯性质量矩；
发电机模型为发电机和变频器组成的变速驱动装置，用一阶迟滞模型反应该变频器的转矩响应：
Figure FDA0003732159410000023
其中，Qd为所请求的扭矩，Qg为气隙扭矩，τe为一阶迟滞的时间常数，S表示该传递函数是原函数经过拉普拉斯变换后的表达式；
在Simulink软件中使用Logic and Bit Operations模块、Lookup Tables模块、MathOperations模块以及Port&Subsystems模块，实现对机组的启机、并网、停机、偏航对风、机组监测与保护、故障诊断与智能维护的控制逻辑，最终搭建出控制器模型。
4.根据权利要求2所述的一种风机多源模型耦合协同仿真方法，其特征在于，在步骤S3中，具体执行以下操作：
1)设定共享内存区域，Bladed进程和Simulink进程能够对该共享内存区域的地址同时访问，能同时对共享内存区域中的数据进行读写，首先利用CreateFileMapping()函数对Bladed计算结果文件或Simulink指令文件进行存取操作，在该函数中对Bladed计算结果文件或Simulink指令文件进行操作等同于直接对共享内存数据进行操作；在成功创建共享文件后，CreateFileMapping()函数产生一个file-mapping核心对象，file-mapping核心对象记录了内存数据与共享文件数据的映射关系，再利用MapViewFile()函数，产生一个指向可用内存的指针，指向被映射的视图开头；
2)使用共享内存，在完成共享内存的定义与创建并建立了映射关系后，利用file-mapping核心对象使得其它需要访问共享内存的进程都能通过名称找到共享内存，并将其映射到各进程的地址空间中，构建成视图；再利用OpenFileMapping()函数，向访问共享内存的进程返回一个句柄，代表一个file-mapping核心对象，同时进程通过调用OpenFileMapping()函数获得一个指向共享内存的指针；
3)清理共享内存，在完成对共享内存的读写操作后，调用UnmapViewOfFile()函数，将原本由调用MapViewOfView()函数而获得的指针交出，再调用CloseHandle()函数，把file-mapping核心对象的句柄交出；
4)利用MATLAB软件对描述有共享内存接口调用信息的程序进行编译，生成MEX文件，利用S-Function将编译通过的MEX文件转换为被风机控制器模型应用的控件资源，使用mdlInitializeSizes()函数确定风机控制器模型的输入输出口的个数，在MEX文件中使用mdlInitializeSampleTimes()函数对S-Function函数的采样周期作出定义，最终采用内存共享技术来完成Bladed进程和Simulink进程之间的通信，确保数据的传输可靠。
5.一种风机多源模型耦合协同仿真平台，其特征在于，包括：
基于Matlab的风机数学建模模组，用于建立风机的气动、传动链、塔架和发电机的数学模型；
基于Bladed的风机运行参数建模模组，用于对风机建立运行参数模型；
基于SolidWorks的风机可视化模组，用于对风机建立三维实体模型；
基于Simulink的风机控制器建模模组，用于对机组的启机、并网、停机、偏航对风、机组监测与保护、故障诊断与智能维护的控制逻辑；
通用模型资源实体模组，用于建立数据共享机制，最终进行各模型的耦合仿真。
6.一种存储有指令的非暂时性计算机可读介质，其特征在于，当所述指令由处理器执行时，执行根据权利要求1至4任意一项所述的风机多源模型耦合协同仿真方法的步骤。
7.一种计算设备，包括处理器以及用于存储处理器可执行程序的存储器，其特征在于，所述处理器执行存储器存储的程序时，实现权利要求1至4任意一项所述的风机多源模型耦合协同仿真方法。
CN202210787705.7A 2022-07-06 2022-07-06 一种风机多源模型耦合协同仿真方法、平台、介质及设备 Withdrawn CN115238409A (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202210787705.7A CN115238409A (zh) 2022-07-06 2022-07-06 一种风机多源模型耦合协同仿真方法、平台、介质及设备
CN202310725246.4A CN116910929A (zh) 2022-07-06 2023-06-19 一种风机多源模型耦合协同仿真方法、平台、介质及设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210787705.7A CN115238409A (zh) 2022-07-06 2022-07-06 一种风机多源模型耦合协同仿真方法、平台、介质及设备

Publications (1)

Publication Number Publication Date
CN115238409A true CN115238409A (zh) 2022-10-25

Family

ID=83671161

Family Applications (2)

Application Number Title Priority Date Filing Date
CN202210787705.7A Withdrawn CN115238409A (zh) 2022-07-06 2022-07-06 一种风机多源模型耦合协同仿真方法、平台、介质及设备
CN202310725246.4A Pending CN116910929A (zh) 2022-07-06 2023-06-19 一种风机多源模型耦合协同仿真方法、平台、介质及设备

Family Applications After (1)

Application Number Title Priority Date Filing Date
CN202310725246.4A Pending CN116910929A (zh) 2022-07-06 2023-06-19 一种风机多源模型耦合协同仿真方法、平台、介质及设备

Country Status (1)

Country Link
CN (2) CN115238409A (zh)

Also Published As

Publication number Publication date
CN116910929A (zh) 2023-10-20

Similar Documents

Publication Publication Date Title
He et al. Multi-body dynamics modeling and TMD optimization based on the improved AFSA for floating wind turbines
Jin et al. Optimization of tuned mass damper parameters for floating wind turbines by using the artificial fish swarm algorithm
Khaouch et al. Mechatronic modeling of a 750 kW fixed-speed wind energy conversion system using the Bond Graph Approach
CN108256210B (zh) 一种地震作用下的海上风机整体耦合分析方法
Abdel Hafeez et al. Flutter limit investigation for a horizontal axis wind turbine blade
Heinz Partitioned fluid-structure interaction for full rotor computations using CFD
Jonkman et al. Functional requirements for the WEIS toolset to enable controls co-design of floating offshore wind turbines
TW201943951A (zh) 用於設計風力機基座之方法及系統
Xie et al. Structural vibration control for the offshore floating wind turbine including drivetrain dynamics analysis
CN117313526A (zh) 一种漂浮式海上风力机的耦合协同仿真计算方法
CN114281029A (zh) 风力发电机的数字孪生系统和方法
Marten et al. Validation and comparison of a newly developed aeroelastic design code for VAWT
CN115238409A (zh) 一种风机多源模型耦合协同仿真方法、平台、介质及设备
CN116484671A (zh) 风力机叶片全耦合动力响应精细化分析方法
CN110807278A (zh) 一种齿轮系统的三维实体单元建模方法
Pao et al. Control co-design of wind turbines
Solcia et al. Efficient multirate simulation of complex multibody systems based on free software
CN111639428B (zh) 带襟翼的风力发电机的动力学仿真方法、设备及介质
El Beshbichi et al. Modelica‐AeroDyn: Development, benchmark, and application of a comprehensive object‐oriented tool for dynamic analysis of non‐conventional horizontal‐axis floating wind turbines
Bir Structural dynamics verification of rotorcraft comprehensive analysis system (RCAS)
Saleh et al. Computational Design Scheme for Wind Turbine Drive‐Train Based on Lagrange Multipliers
CN106960077A (zh) 一种风机齿轮箱动力学模型的建模与仿真方法
CN111159929A (zh) 复合材料叶片结构动力学数值模拟方法
CN205450761U (zh) 一种用于风电机组电控系统测试的硬件在环实验装置
Nejadkhaki et al. Integrative modeling platform for design and control of an adaptive wind turbine blade

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
WW01 Invention patent application withdrawn after publication

Application publication date: 20221025

WW01 Invention patent application withdrawn after publication