CN1936531A - 风力发电机组风轮功率输出特性模拟试验方法 - Google Patents
风力发电机组风轮功率输出特性模拟试验方法 Download PDFInfo
- Publication number
- CN1936531A CN1936531A CN 200610048043 CN200610048043A CN1936531A CN 1936531 A CN1936531 A CN 1936531A CN 200610048043 CN200610048043 CN 200610048043 CN 200610048043 A CN200610048043 A CN 200610048043A CN 1936531 A CN1936531 A CN 1936531A
- Authority
- CN
- China
- Prior art keywords
- wind
- unit
- wheel
- formula
- wind speed
- 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.)
- Granted
Links
Images
Landscapes
- Wind Motors (AREA)
Abstract
一种风力发电机组风轮功率输出特性模拟试验方法属于风力发电机技术领域。本发明方法通过模拟试验装置进行试验,试验装置包括上位机、模拟风轮和试验对象电机,其中模拟风轮包括变频器和异步电机,变频器分别与上位机、异步电机和试验对象电机相连,异步电机和试验对象电机相连。该方法包括建立风的数学模型、确定风轮动力特性方程和计算传递函数,由上位机通过软件模块实现。本发明方法不必应用于现场,只需在实验室就可测试发电机和控制器的性能,提高了产品质量,加快了产品成熟周期并节约大量资金。
Description
技术领域
本发明属于风力发电机技术领域,特别涉及一种风力发电机组风轮功率输出特性模拟试验方法。
技术背景
目前,中国的并网型风力发电机组生产厂已有30多个,产品的定性和设计没有科学的试验测试手段,很难生产出高质量的风力发电机组与国外的风力发电机组制造厂商竞争。这样就制约了产品的迅速产业化和商业化,失去了占领技术制高点和占领市场的机会。目前在国外也仅仅只有小功率普通拖动发电系统和仿真软件,没有实际模拟风轮输出特性的设备和实现的理论方法。
发明内容
针对现有技术中存在的问题,本发明提供一种风力发电机组风轮功率输出特性模拟试验方法。本发明方法是对设计完成的风力发电机组产品进行产品定型和优化设计所必需的试验。本发明方法所测得的数据可作为技术设计的依据。
本发明方法所依赖的装置包括上位机1、模拟风轮2和试验对象电机5,其中模拟风轮包括变频器3和异步电机4,如图1所示,变频器3分别与上位机1、异步电机4和试验对象电机5相连,异步电机4和试验对象电机5相连。本发明方法的试验对象可为任何原型风力机风轮特性模拟输出,如用于变速恒频发电系统的风力机仿真拖动模拟,风力发电机组的出厂性能试验测试系统的模拟拖动等。
本发明方法包括建立风的数学模型和求出风轮动力特性方程。下面进行详细介绍:
一、建立风的数学模型
(1)基本风。基本风在风轮机正常运行过程中一直存在,它决定了风力发电机向系统输送额定功率的大小,基本上反映了风电场平均风速的变化,它可以由风电场测风所得的韦布尔分布参数确定:
K——韦布尔分布的形状参数;Г——伽马函数。
对于基本风模型,在实验室模拟时无需考虑实际风场风速的分布情况,只需考虑作用在风轮机上的平均风速是多少,所以基本风模型指定一个具体的平均风速值即可。
(2)阵风。在风速变化的过程中,描述风速突然变化的特性,一般用阵风来表示,可在基本风上叠加一阵风分量VWG。
式中:VWG——阵风风速,单位为m/s;t——时间;TG——周期,单位为s;
T1G——阵风启动时间,单位为s;maxG——阵风最大值,单位为m/s。
(3)渐变风。对风速的渐变特性可以用渐变风成份来模拟,可在基本风上叠加一个渐变风分量VWR。
式中:VWR——渐变风速,单位为m/s;maxR——渐变风最大值,单位为m/s;
T1R——起始时间,单位为s;T2R——终止时间,单位为s;
TR——保持时间,单位为s。
(4)随机噪声风速。为反映风速变化的随机特性,在风速模拟中可以用随机噪声风速成份来表示,可在基本风上叠加一个随机分量VWN。
式中:ω——角频率;ωi——初象角;_i——均匀分布的随机变量,_i=0~2π;
i——序列;KN——地表粗糙系数,一般取KN=0.004;
F——扰动范围,单位为m2;μ——相对高度的平均风速,单位为m/s。
综合上述四种风速成份,模拟实际作用在风轮机上的风速为
式中:VW——作用在风轮机上的风速,单位为m/s。
二、根据风轮不同的特性模拟要求,将风的模型特性带入风轮动力特性表达式中,求出其特性方程
(1)风能利用系数
风能利用系数反映风轮机捕获风能效率的高低,贝兹(Betz)极限证明了理想风轮的最大风能利用效率
实际应用当中达不到最大值。对于变桨距风轮机来说,影响风能利用系数的有两个主要因素:1、风轮机叶片的叶尖速比λ(风轮机叶片尖端的旋转线速度与风速的比值);2、叶片的桨距角β(叶片的翼弦与水平面的夹角),如图2所示。
升力表达公式如下:
式中:Fy——升力;cy——升力系数;S——叶片在气流方向的投影面积,单位为m2;
V——气流速度,单位为m/s;ρ——空气密度,单位为kg/m3;
对于流线型翼型中cy与攻角α的关系如图3所示。
风轮在启动后以某一速度旋转时叶片的受力情况有所改变,这时相对于叶片上叶尖的气流速度将是垂直于风轮旋转面的气流速度V与该叶片的旋转线速度ωr的矢量和ω,而且这时的攻角已经不是V与翼弦的夹角β,而是ω与翼弦的夹角α。如图4所示。
将风轮所受升力代入式(6)得:
升力系数用以下公式来近似:
cp=-1.6x2+1.9x+0.02 (x=β-arctanλ) (9)
(2)风轮机的机械输出功率
风轮的动力特性方程,通过风轮机旋转面的空气所蕴含的动能为P:
风轮机从风中捕获的机械能量是空气动能的一部分。风轮机的机械输出功率为Pw:
式中:ρ——空气密度,单位为kg/m3;D——风轮直径,单位为m;
A——风轮机的截面积,单位为m2;V——空气速度,单位为m/s;
Cp——风能利用系数。
(3)风轮机的机械转矩TW
式中:Ω——转子转速;N——发电机转动速度。
(4)风力机风轮的输出特性方程
风轮转速的输出状态方程为
式中:
W=V——风速;β——桨距角;Ω——转子转速;
Kd——传动链刚度;Irot——风轮转动惯量;
Igen——发电机转动惯量;Cd——阻尼系数;
三、求出传递函数
传递函数G:
G=C(sI-A)-1B (14)
式中:A、B、C——矩阵,s——拉氏变换因子。
当前状态的传递函数值,作为异步电机的控制输入,迅速调节实现实验模拟状态的输出结果。
本发明方法最终由上位机实现,上位机包括三个模块:界面显示模块、数值处理模块和上、下位机之间的通讯模块,通讯模块分别与界面显示模块、数值处理模块连接,如图5所示。上位机在监控界面上显示数据和模拟过程曲线,其过程包括以下步骤:
步骤一、输入命令,当输入命令为模式选取时进入步骤四;
步骤二、开始;
步骤三、初始化;
步骤四、试验模式选取,即选取风的模型;
步骤五、读取数据,数据包括风速、环境温度、风轮转矩、风轮转速、风轮输出功率、异步电机功率;
步骤六、计算风轮转速输出状态方程中矩阵A、B、C;
步骤七、由式(14)计算传递函数;
步骤八、显示模块显示输出;
步骤九、通讯模块进行通讯;
步骤十、D/A输出给变频器,对异步电机进行调整,返回步骤一。
其中显示模块显示输出过程包括以下步骤:
步骤一、输入用户命令;
步骤二、与风力发电机参数数据库连接;
步骤三、调用函数(连接数据库或从数据库读数或刷新界面或写数据);
步骤四、数据处理;
步骤五、显示。
通讯模块通讯过程包括以下步骤:
步骤一、读数据;
步骤二、判断是否为周期数据,是,进入步骤五,否则进入步骤三;
步骤三、发送非周期数据读取命令;
步骤四、读数据;
步骤五、初始化;
步骤六、切换到Sercos2段配置运行参数;
步骤七、运行参数读取;
步骤八、切换到Sercos4段周期通讯;
步骤九、周期通讯;
步骤十、周期运行。
本发明方法不必应用于现场,只需在实验室就可测试发电机和控制器的性能,提高了产品质量,加快了产品成熟周期并节约大量资金。
附图说明
图1为本发明方法所依赖的装置结构示意图;
图2为桨距角示意图;
图3为翼型升力曲线图;
图4为旋转桨叶在风中的受力图;
图5为本发明方法的上位机中模块间的关系图;
图6为本发明方法的流程图;
(a)为本发明方法主体流程图;
(b)为本发明方法界面显示模块流程图;
(c)为本发明方法通讯模块流程图;
图7为本发明具体实例的主界面图;
图8为本发明具体实例的风速设置界面图;
图9为本发明具体实例的参数设置界面图;
图10为本发明具体实例的状态显示界面图。
图中1为上位机,2为模拟风轮,3为变频器,4为异步电机,5为试验对象电机,6为风轮旋转面,7为界面显示模块,8为通讯模块,9为数值处理模块。
具体实施方式
下面结合附图对本发明做进一步说明。
当风速为12m/s,转子转速Ω=1680r/min,桨距角β=2度;风机以1.0的参数为例,半径R=30.3,风轮转速N=1680/78.16=21.5RPM,λ=2*3.14*NR/W,环境温度对应空气密度为1.025P时:
一、建立风的数学模型
本实例中建立风的模型为基本风。基本风由风电场测风所得的韦布尔分布参数确定:
K——韦布尔分布的形状参数;Г——伽马函数。
二、根据风轮不同的特性模拟要求,将风的模型特性带入风轮动力特性表达式中,求出其特性方程
(1)风能利用系数
cp=-1.6x2+1.9x+0.02=0.399
(2)风轮机的机械输出功率
(3)风轮机的机械转矩TW
(4)风力机风轮的输出特性方程
求得矩阵A、B、C为
得到状态方程为:
三、求出传递函数
上述方法最终由上位机实现,其过程包括以下步骤:
步骤一、输入命令,当输入命令为模式选取时进入步骤四;
步骤二、开始;
步骤三、初始化;
步骤四、试验模式选取,即选取风的模型;
步骤五、读取数据,数据包括风速、环境温度、风轮转矩、风轮转速、风轮输出功率、异步电机功率;
步骤六、计算风轮转速输出状态方程中矩阵A、B、C;
步骤七、由式(14)计算传递函数;
步骤八、显示模块显示输出;
步骤九、通讯模块进行通讯;
步骤十、D/A输出给变频器,对异步电机进行调整,返回步骤一。
其中显示模块显示输出过程包括以下步骤:
步骤一、输入用户命令;
步骤二、与风力发电机参数数据库连接;
步骤三、调用函数(连接数据库或从数据库读数或刷新界面或写数据);
步骤四、数据处理;
步骤五、显示。
通讯模块通讯过程包括以下步骤:
步骤一、读数据;
步骤二、判断是否为周期数据,是,进入步骤五,否则进入步骤三;
步骤三、发送非周期数据读取命令;
步骤四、读数据;
步骤五、初始化;
步骤六、切换到Sercos2段配置运行参数;
步骤七、运行参数读取;
步骤八、切换到Sercos4段周期通讯;
步骤九、周期通讯;
步骤十、周期运行。
本实例中上位机选用KZC064系列嵌入式一体化工控机,,采用Borland C++Builder 5.0开发的基于Windows的可视化控制软件,它为用户提供了一个友善的人机交互界面。
上位机包括图形界面模块、数值处理模块和上、下位机之间的通讯模块。
图形界面部分提供了人机交互的界面。通过图形界面,用者可以设定实验的方案,也可以观察实验的进程,并在实验结束后形象的显示系统运行的曲线图。
图7是本实例的主界面图。该界面包括三个子窗体的选择按钮,点击风速设置就可以进入风速模型参数设置的界面,点击参数设置就可以进入风机选型和风机参数设置的界面,点击状态显示就可以进入监测风速、异步电机转速、异步电机转矩、异步电机功率的运行曲线。
图8是本实例的风速设置界面。在这个界面中,我们根据四种风速风况模型,分别输入参数,点击开始就可以得到风速时变曲线;当风速曲线不是很理想时,点击停止按钮,重新设定模拟风速的参数,点击下一步则进行下一个参数设置界面。
图9是本实例的参数设置界面。在该界面中根据桨距角的不同而进行选择,由于风轮机的运行方式不同,所以导致所建的模型有所不同,不同的选择操作界面会链接到不同的运算模型。
图10是本实例的状态显示界面。该界面主要做显示和监控用,前几个界面的参数设定好之后,点击开始就可以同步显示模拟风速曲线、异步电机输出转速曲线、异步电机输出转矩曲线、异步电机输出功率曲线。可以进行实时的对比分析,如果发现存在问题或当前曲线非常成功就可以点击暂停按钮,记录当前数据,点击继续,则继续分析监测。
Claims (7)
1、一种风力发电机组风轮功率输出特性模拟试验方法,通过模拟试验装置进行试验,试验装置包括上位机、模拟风轮和试验对象电机,其中模拟风轮包括变频器和异步电机,变频器分别与上位机、异步电机和试验对象电机相连,异步电机和试验对象电机相连,其特征在于该试验方法包括建立风的数学模型、确定风轮动力特性方程和计算传递函数,由上位机通过软件模块实现。
2、如权利要求1所述的一种风力发电机组风轮功率输出特性模拟试验方法,其特征在于所述的建立风的数学模型包括基本风,其数学模型为:
K——韦布尔分布的形状参数,Г——伽马函数;
阵风数学模型为:
式中:
V WG——阵风风速,单位为m/s,
T G——周期,单位为s,t——时间;
T 1G——阵风启动时间,单位为s,max
G——阵风最大值,单位为m/s;
渐变风数学模型为:
式中:
V WR——渐变风速,单位为m/s,max
R——最大值,单位为m/s,
T 1
R ——起始时间,单位为s,
T 2
R ——终止时间,单位为s,
T R ——保持时间,单位为s;
随机噪声风数学模型为:
式中:_i——均匀分布的随机变量,_i=0~2π,KN——地表粗糙系数,一般取KN=0.004,
F——扰动范围,单位为
m 2,
μ——相对高度的平均风速,单位为m/s;
综合上述四种风速成份,模拟实际作用在风轮机上的风速为
式中:VW——作用在风轮机上的风速,单位为m/s。
3、如权利要求1所述的一种风力发电机组风轮功率输出特性模拟试验方法,其特征在于所述的风轮动力特性方程如下:
式中:
T w ——风轮机的机械转矩,W=V——风速,β——桨距角,Ω——转子转速,
Kd——传动链刚度,Irot——风轮转动惯量,Igen——发电机转动惯量,
Cd——阻尼系数。
4、如权利要求1所述的一种风力发电机组风轮功率输出特性模拟试验方法,其特征在于所述的传递函数如下:
G=C(sI-A)-1B
式中:A、B、C——矩阵,s——拉氏变换因子。
5、如权利要求1所述的一种风力发电机组风轮功率输出特性模拟试验方法,其特征在于所述的上位机通过软件模块实现的过程按以下步骤执行:
步骤一、输入命令,当输入命令为模式选取时进入步骤四;
步骤二、开始;
步骤三、初始化;
步骤四、试验模式选取,即选取风的模型;
步骤五、读取数据,数据包括风速、环境温度、风轮转矩、风轮转速、风轮输出功率、异步电机功率;
步骤六、计算风轮转速输出状态方程中矩阵A、B、C;
步骤七、由式G=C(sI-A)-1B计算传递函数;
步骤八、显示模块显示输出;
步骤九、通讯模块进行通讯;
步骤十、D/A输出给变频器,对异步电机进行调整,返回步骤一。
6、如权利要求5所述的一种风力发电机组风轮功率输出特性模拟试验方法,其特征在于步骤八中的显示模块显示输出过程包括以下步骤:
步骤一、输入用户命令;
步骤二、与风力发电机参数数据库连接;
步骤三、调用函数,调用连接数据库或从数据库读数或刷新界面或写数据;
步骤四、数据处理;
步骤五、显示。
7、如权利要求5所述的一种风力发电机组风轮功率输出特性模拟试验方法,其特征在于步骤九中的通讯模块通讯过程包括以下步骤:
步骤一、读数据;
步骤二、判断是否为周期数据,是,进入步骤五,否则进入步骤三;
步骤三、发送非周期数据读取命令;
步骤四、读数据;
步骤五、初始化;
步骤六、切换到Sercos2段配置运行参数;
步骤七、运行参数读取;
步骤八、切换到Sercos4段周期通讯;
步骤九、周期通讯;
步骤十、周期运行。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2006100480432A CN100516810C (zh) | 2006-10-19 | 2006-10-19 | 风力发电机组风轮功率输出特性模拟试验方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2006100480432A CN100516810C (zh) | 2006-10-19 | 2006-10-19 | 风力发电机组风轮功率输出特性模拟试验方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1936531A true CN1936531A (zh) | 2007-03-28 |
CN100516810C CN100516810C (zh) | 2009-07-22 |
Family
ID=37954157
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2006100480432A Expired - Fee Related CN100516810C (zh) | 2006-10-19 | 2006-10-19 | 风力发电机组风轮功率输出特性模拟试验方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100516810C (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101545953B (zh) * | 2009-04-24 | 2011-01-05 | 重庆大学 | 小型风力发电机性能测试方法及其专用测试平台 |
CN102175373A (zh) * | 2011-03-11 | 2011-09-07 | 西安交通大学 | 一种大型风力发电机组全功率测试系统 |
CN102841545A (zh) * | 2012-09-27 | 2012-12-26 | 武汉华大电力自动技术有限责任公司 | 一种风力特性的仿真控制装置及其仿真控制装方法 |
CN105679160A (zh) * | 2016-03-31 | 2016-06-15 | 华北电力大学 | 风电机组运行控制教学模拟装置及其模拟方法 |
CN106405220A (zh) * | 2016-10-20 | 2017-02-15 | 盐城工学院 | 增程式辅助发电系统输出功率测试装置及系统 |
CN110262462A (zh) * | 2019-07-05 | 2019-09-20 | 商飞信息科技(上海)有限公司 | 一种基于实时控制器的风机并网测试和研发平台 |
-
2006
- 2006-10-19 CN CNB2006100480432A patent/CN100516810C/zh not_active Expired - Fee Related
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101545953B (zh) * | 2009-04-24 | 2011-01-05 | 重庆大学 | 小型风力发电机性能测试方法及其专用测试平台 |
CN102175373A (zh) * | 2011-03-11 | 2011-09-07 | 西安交通大学 | 一种大型风力发电机组全功率测试系统 |
CN102841545A (zh) * | 2012-09-27 | 2012-12-26 | 武汉华大电力自动技术有限责任公司 | 一种风力特性的仿真控制装置及其仿真控制装方法 |
CN105679160A (zh) * | 2016-03-31 | 2016-06-15 | 华北电力大学 | 风电机组运行控制教学模拟装置及其模拟方法 |
CN105679160B (zh) * | 2016-03-31 | 2018-11-30 | 华北电力大学 | 风电机组运行控制教学模拟装置及其模拟方法 |
CN106405220A (zh) * | 2016-10-20 | 2017-02-15 | 盐城工学院 | 增程式辅助发电系统输出功率测试装置及系统 |
CN106405220B (zh) * | 2016-10-20 | 2018-12-14 | 盐城工学院 | 增程式辅助发电系统输出功率测试装置及系统 |
CN110262462A (zh) * | 2019-07-05 | 2019-09-20 | 商飞信息科技(上海)有限公司 | 一种基于实时控制器的风机并网测试和研发平台 |
Also Published As
Publication number | Publication date |
---|---|
CN100516810C (zh) | 2009-07-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Adaramola et al. | Experimental investigation of wake effects on wind turbine performance | |
Plourde et al. | Simulations of three-dimensional vertical-axis turbines for communications applications | |
CN103573552B (zh) | 风力涡轮机及其降低转子不平衡的控制方法 | |
CN100516810C (zh) | 风力发电机组风轮功率输出特性模拟试验方法 | |
CN102879199B (zh) | 一种全工况风力机模拟器 | |
CN106777499B (zh) | 一种双馈异步风力发电机组的整机动态建模方法 | |
CN101430246A (zh) | 风力发电模拟试验平台 | |
CN111859650B (zh) | 基于在线式联合仿真的风电机组传动链虚拟地面试验方法 | |
CN104675629A (zh) | 一种变速风力发电机组的最大风能捕获方法 | |
Jiang et al. | Experimental evaluation of a tree-shaped quad-rotor wind turbine on power output controllability and survival shutdown capability | |
Rogowski et al. | Numerical analysis of a small-size vertical-axis wind turbine performance and averaged flow parameters around the rotor | |
Abo-Khalil et al. | Dynamic modeling and control of wind turbines for grid-connected wind generation system | |
Gaonkar et al. | An investigation on the performance and viability of a hybrid twisted blade profile for a horizontal axis micro wind turbine | |
CN201344857Y (zh) | 风力发电模拟试验平台 | |
CN115202238B (zh) | 海上风电机组机械-电气联合硬件在环高精度仿真方法 | |
Corbus et al. | Analysis and comparison of test results from the small wind research turbine test project | |
CN203570507U (zh) | 大型风力发电机组缩比控制系统 | |
Pitance et al. | Experimental validation of Pharwen code using data from Vertical-axis wind turbines | |
CN109611268A (zh) | 一种双叶轮水平轴风力机设计优化方法 | |
Islam et al. | Optimum blade design of pinwheel type horizontal axis wind turbine for low wind speed areas | |
Shi | Model-based detection for ice on wind turbine blades | |
Zhang | Analysis and optimisation of a novel wind turbine | |
Tchakoua et al. | Basis of theoretical formulations for new approach for modelling darrieus-type vertical axis wind turbine rotors using electrical equivalent circuit analogy | |
Fazylova et al. | Analysis of rotor'critical mode of operation to be employed in the design of a wind generation control unit | |
Lindenburg | Comparison of PHATAS Versions and the Wind turbine Module |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20090722 Termination date: 20131019 |