风力发电机运行仿真.docx

《风力发电机运行仿真.docx》由会员分享，可在线阅读，更多相关《风力发电机运行仿真.docx（41页珍藏版）》请在第一文库网上搜索。

1、基于MATLAB的“风力发电机运行仿真”软件设计摘要关键词1前言1.1 建模仿真的发展现状20世纪50-60年代，自动控制领域普遍采用计算机模拟方法研究控制系统动态过程和性能。“计算机模拟”实质上是数学模型在计算机上的解算运行，当时的计算机是模拟计算机，后来发展为数字计算机。1961年G.W.Morgenthler首次对仿真一词作了技术性的解释，认为“仿真”是指在实际系统尚不存在的情况下，对于系统或活动本质的复现。目前，比较流行于工程技术界的技术定义是系统仿真是通过对系统模型的实验，研究一个存在的或设计中的系统。仿真的三要素之间的关系可用三个基本活动来描述。如图1图1系统仿真三要素之间的关系2

2、0世纪50年代初连续系统仿真在模拟计算机上进行，50年代中出现数字仿真技术，从此计算机仿真技术沿着模拟仿真和数字仿真两个方面发展。60年代初出现了混和模拟计算机，增加了模拟仿真的逻辑控制功能，解决了偏微分方程、差分方程、随机过程的仿真问题。从60-70代发展了面向仿真问题的仿真语言。20世纪80年代末到90年代初，以计算机技术、通讯技术、智能技术等为代表的信息技术的迅猛发展，给计算机仿真技术在可视仿真基础上的进一步发展带来了契机，出现了多媒体仿真技术。多媒体仿真技术充分利用了视觉和听觉媒体的处理和合成技术，更强调头脑、视觉和听觉的体验，仿真中人与计算机交互手段也更加丰富。80年代初正式提出了

3、“虚拟现实” 一词。虚拟现实是一种由计算机全部或部分生成的多维感觉环境，给参与者产生视觉、听觉、触觉等各种感官信息，使参与者有身临其境的感觉，同时参与者从定性和定量综合集成的虚拟环境中可以获得对客观世界中客观事物的感性和理性的认识。图2体现了仿真科学与技术的发展进程。A仿真对象的发展演化性、自治性、涌现性仿真系统功能、性能仿百控制系统仿真几何、力学仿真力学、物理控制理论系统论复杂性科学仿真研究领域的扩展 一图2仿真科学与技术的发展以美国为代表的发达国家高度重视仿真技术的发展和应用。美国等西方国家除军事用途外的其它行业中的仿真技术及应用都居于世界领先水平，如飞行模拟器、车辆运输仿真、电力系统、石

4、油化工仿真系统等。经过几个五年计划的努力，我国仿真技术得到了快速的发展，并取得了突破性成果，和长足的进步。在某些方面达到了国际先进水平。但总体水平，特别是应用水平与发达国家比较还有差距，需要进一步努力，加速发展仿真技术以缩小差距1.2 本仿真软件简介2风力发电机各部分数学模型及仿真2. 1风力机风能利用系数(功率系数)cP是指单位时间内风力机所获得的能量与风能之比。它是评定风力机气动特性优劣的只要参数，其定义式：。一 式中：P为风力机的功率，单位是W；P为空气密度，单位是kg/m3；S为风轮的扫风面积，单位是n?；V为来流风速，单位是m/s在设计Savonius风力机时要考虑两个重要的结构参数

5、：一个是重叠比0L(Overlap ratio), 一个是高径比 AP (Aapect ratio)：(2)(3)OL = S/dAP=H/d叶片重叠比对Savonius风力机的各种性能影响很大。如图4的风洞试验数据所示，具有不同的叶片重叠比的风力机的最大功率系数相差很大，合理设计叶片重叠比可以改善风力机的静态启动特性，对风力机的动态力矩变化的战俘和相位也具有一定的影响。尖速比h图4具有不同重叠比的Savonius风力机的性能叶片高径比也对风力机的性能影响很大，一般来说叶片高径比越大风力机性能越好。目前实际应用中的Savonius风力机的叶片高径比一般为14,准确数值要根据设计目标、成本和安装

6、地点的风况特点来决定。叶尖速比人是叶片的叶尖圆周速度与风速之比，用来描述风轮在不同风速中的状态：(4)2/rRn coRz =式中：n为风轮转速，单位是r/s；3为风轮角频率，单位是rad/s；R为风轮半径，单位是m；V为上游风速，单位是m/s风力机通过叶片捕获风能，将风能转换为作用的发电机转子上的机械能，将吸收的叶片转矩为作用在发电机转子上的机械转矩。风力机吸收功率可以表示为风速的函数，其模型表示为：P =（5）所以风力机的机械转矩为:(6)T = - = -jipCpRy 60 2Vp %其中相关参数的设定会影响风力机的输出效率。对于风力机建模，主要有两种方式，一种是对发电机的实测数据的查

7、表法,另一种是根据相关的公式进行垂直轴风力发电机输出参数建模。风能利用系数Cp的函数曲线如图6所示，由图可知，当叶尖速比在1左右时，输出效率最大。为了简化模型，我们在仿真过程中设定叶尖速比在速度范围之内为恒定值。风机输出效率及转矩效率图6最佳Savonius型风力机输出效率及转矩效率随叶尖速比的变化通过图6可以拟合曲线，得到风能利用系数Cp的函数：Cp = /(2) = -0.3656x (等 了 + 0.6505 x (竿)(7)由此可得输出功率为：P = /(A)xlp4V3= -0.3656 x ()2 + 0.6505 x ()x-pAV3 (8)22根据数学模型，进行Simulink

8、仿真。输入为风速V,空气密度P和发电机电磁转矩T_em；叶尖速比人和风力机受风面积A为常数。输出为风力机输出功率P，输出转矩Te,输出转速3和发电机输入转速3_em。同时我们还设定了风力机的启动风速和最高风速。模型如图7所示，并对其进行了封装。封装界面如图8。Fcn2Dividel图7风力机仿真模型图8风力机封装界面设定输入风速为15m/s,空气密度为l.29kg/nA仿真时间为10秒，得输出曲线如图9。横坐标为时间，纵坐标分别为功率、转矩、转速、发电机输入转速。风力机输出功率约为620Wo图9风力机仿真输出曲线2.2发电机2. 2.1永磁发电机永磁同步发电机由绕线转子同步发电机发展而来，定子

9、与普通同步发电机基本相同，转子为永磁体，一般无阻尼绕组，因此不存在励磁绕组的铜损耗，同时无需外部提供励磁电源，可以提高效率；转子上没有滑轮，可以提高系统的稳定性。风力机输出的机械转矩带动发电机转子转动。永磁同步发电机的转子为永磁式结构，转子的磁链由永磁体决定。我们将定子电压在dpO同步旋转坐标系下进行分解，其中，同步旋转坐标系的d轴是转子磁链的方向。在此基础上建立发电机定子电压的d轴和q轴分量的表达式：(9)ud= Rjd + 一 叫儿atd儿.%= R/q + /at式中：id和iq分别为发电机的d轴和q轴电流；Ud和uq分别为定子电压Eg的d轴和q轴分量;入d和入q分别是d轴和q轴的磁链R

10、a为定子电阻；3e为电角频率定义磁链的d轴和q轴的分量的表达式为：(10)式中：Ld和Lq分别为发电机的d轴和q轴电感；入o为永磁体产生的磁链定义q轴的反电势eq=3。入。，而d轴的反电势ed=O,因为发电机的转子为对称结构，这里我们可以假设发电机的d轴和q轴的电感相等，即Ld=Lq=L。将式（8）带入到式（7）并整理化简得到：Ra .1dtd_、dt(11)一_ld + (Oelq + UdRa ./.1 i 、1一_+彳) + 7%因为Ld=Lq=L,则永磁同步发电机的电磁转矩表达式为:1=1.5也一4)* +翱= 1.5%翱(12)式中1为极对数。MATLAB中提供的永磁发电机模型当其输

11、入为负时，作为发电机使用。输入有转矩Tm和转速3两种选择。本仿真软件中，我们选择转速输入。参数设置如图13所示，从上到下依次表示：定子相电阻，d轴和q轴定子电感，选择机器常数，感应磁链，电压常数，转矩常数，极对数，初始条件（转速、角度、电流）。通过调整感应磁链和极对数得到较为理想的曲线。图14和图15分别为仿真模型和仿真结果。1000500-500-1000ConfigurationParameters AdvancedStator phase resistance Rs (ohm):2. 8750Inductances Ld(K) Lq(H):8. 5e-3,8. 5e-3Specify：

12、Fluz linkage established by maenets (V. s)Fluz linkage established by macnets (V. s):1. 5Voltace Constant (V peak L-L / krpm):4353. 1185Torque Constant (N. m / A peak):36Pole pairs p ():16Initial conditions wn(rad/s) thetam(det) ia/ ib(A):0j0j 0, 0图13永磁同步发电机参数设置图14永磁同步发电机模型HIIIWiWIMiBllilMliiliN用川川川

13、川川川川图15发电机输出电压曲线2. 2. 2电励磁同步发电机风力发电中所用的同步发电机绝大部分是三相同步电机，其输出联接到邻近的三相电网或输配电线。普通三相同步发电机的原理结构如图5所示。在定子铁心上有若干槽，槽内嵌有均匀分布的在空间彼此相隔12()。电角的三相电枢绕组aa、bb和cc，转子上装有磁极和励磁绕组，当励磁绕组通以直流电流If后，电机内产生磁场。转子被风力机带动旋转，则磁场与定子三相绕组之间有相对运动，从而在定子三相绕组中感应出三个幅值相同，彼此相隔12()。电角的交流电势。这个交流电势的频率f决定于电机的极对数p和转子转速n,即f=pn/6()o还是选择转速3输入。参数设置如图

14、16所示，参数依次为：视在功率、线电压、频率、励磁电流，电子电阻、漏磁电感、d轴和q轴电感，定子的漏磁回路电阻、漏磁电感，阻尼器的d轴和q轴电阻和漏磁电感，极对数，初始值(初始速度偏差、电角度、线电流、相角、初始励磁电压)。图5三相同步发电机结构原理图我们将定子电压在dpO同步旋转坐标系下进行分解，则定子方程可表示为:以尸R/q +(11)“odt定义磁链的表达式：= L,id + M,if + MDiD = L3q + MJq4 = L/o式中：ML为定子绕组和励磁绕组之间的互感系数；MJ为定子绕组和直轴阻尼绕组之间的互感系数；MJ为定子绕组和交轴阻尼绕组之间的互感系数发电机的电磁转矩为：-=L5%(一1X4)(12)(13)MATLAB中的同步电机有国际标准单位和标幺值两种，我们选择的是国际标准单位值的。输入分功率Pm和转速3两种，为了与永磁电机的输入统一，我们Conficuration Parameters AdvancedNominal power, voltaee, frequency, field current Pn(VA) Vn(Vrms) fn(Hz) ifn(A):800 400 50Stator Rs (ohm) LILm。：1.6