《风力发电机运行仿真》由会员分享，可在线阅读，更多相关《风力发电机运行仿真（31页珍藏版）》请在金锄头文库上搜索。
1、基于 MATLAB 的“风力发电机运行仿真” 软件设计摘要关键词A 11 前言1.1 建模仿真的发展现状20世纪 5060 年代, 自动控制领域普遍采用计算机模拟方法研究控制系 统动态过程和性能。“计算机模拟”实质上是数学模型在计算机上的解算运行,当 时的计算机是模拟计算机, 后来发展为数字计算机。1961 年 G.W.Morgenthler 首次对仿真一词作了技术性的解释，认为“仿真”是指在实际系统尚不存在的情 况下，对于系统或活动本质的复现。目前，比较流行于工程技术界的技术定义是 系统仿真是通过对系统模型的实验，研究一个存在的或设计中的系统。仿真的三 要素之间的关系可用三个基本活动来描述。
2、如图1图1系统仿真三要素之间的关系20世纪50年代初连续系统仿真在模拟计算机上进行, 50年代中出现数字仿 真技术, 从此计算机仿真技术沿着模拟仿真和数字仿真两个方面发展。60年代 初出现了混和模拟计算机, 增加了模拟仿真的逻辑控制功能, 解决了偏微分方 程、差分方程、随机过程的仿真问题。从60-70代发展了面向仿真问题的仿真语 言。20世纪80年代末到90年代初, 以计算机技术、通讯技术、智能技术等为 代表的信息技术的迅猛发展, 给计算机仿真技术在可视仿真基础上的进一步发 展带来了契机, 出现了多媒体仿真技术。多媒体仿真技术充分利用了视觉和听觉 媒体的处理和合成技术, 更强调头脑、视觉和听觉
3、的体验, 仿真中人与计算机交 互手段也更加丰富。80 年代初正式提出了“虚拟现实”一词。虚拟现实是一种 由计算机全部或部分生成的多维感觉环境, 给参与者产生视觉、听觉、触觉等 各种感官信息, 使参与者有身临其境的感觉, 同时参与者从定性和定量综合集 成的虚拟环境中可以获得对客观世界中客观事物的感性和理性的认识。图2体现 了仿真科学与技术的发展进程。图2仿真科学与技术的发展以美国为代表的发达国家高度重视仿真技术的发展和应用。美国等西方国家 除军事用途外的其它行业中的仿真技术及应用都居于世界领先水平，如飞行模拟 器、车辆运输仿真、电力系统、石油化工仿真系统等。经过几个五年计划的努力，我国仿真技术得
4、到了快速的发展，并取得了突破 性成果，和长足的进步。在某些方面达到了国际先进水平。但总体水平，特别是 应用水平与发达国家比较还有差距，需要进一步努力，加速发展仿真技术以缩小 差距1.2 本仿真软件简介2 风力发电机各部分数学模型及仿真2.1 风力机风能利用系数(功率系数)Cp是指单位时间内风力机所获得的能量与风能 之比。它是评定风力机气动特性优劣的只要参数，其定义式：1)式中：P为风力机的功率，单位是W；P为空气密度，单位是kg/m3；S 为风轮的扫风面积，单位是 m2V 为来流风速，单位是 m/s在设计 Savonius 风力机时要考虑两个重要的结构参数：一个是重叠比 OL(Overlap
5、ratio), 一个是高径比 AP (Aapect ratio):OL 二 S / d3)AP 二 H / d叶片重叠比对 Savonius 风力机的各种性能影响很大。如图 4 的风洞试验数 据所示，具有不同的叶片重叠比的风力机的最大功率系数相差很大，合理设计叶 片重叠比可以改善风力机的静态启动特性，对风力机的动态力矩变化的战俘和相 位也具有一定的影响。0 了0 Q2 0.4 (Id O.tt0.03尖速比XOL=Q 十 OLxO.2 饥曲图 4 具有不同重叠比的 Savonius 风力机的性能叶片高径比也对风力机的性能影响很大，一般来说叶片高径比越大风力机性 能越好。目前实际应用中的 Sav
6、onius 风力机的叶片高径比一般为 14，准确数值 要根据设计目标、成本和安装地点的风况特点来决定。叶尖速比入是叶片的叶尖圆周速度与风速之比，用来描述风轮在不同风速中 的状态：4)式中：n 为风轮转速，单位是 r/s；3为风轮角频率，单位是rad/s；R 为风轮半径，单位是 m；V 为上游风速，单位是 m/s风力机通过叶片捕获风能，将风能转换为作用的发电机转子上的机械能，将 吸收的叶片转矩为作用在发电机转子上的机械转矩。风力机吸收功率可以表示为风速的函数，其模型表示为：P = 1 pC AV3(5)2 p所以风力机的机械转矩为：P1=一兀pC R32 p其中相关参数的设定会影响风力机的输出效
7、率。对于风力机建模，主要有两种方式，一种是对发电机的实测数据的查表法， 另一种是根据相关的公式进行垂直轴风力发电机输出参数建模。风能利用系数 Cp的函数曲线如图6所示，由图可知，当叶尖速比在1左右时，输出效率最大。 为了简化模型，我们在仿真过程中设定叶尖速比在速度范围之内为恒定值。图 6 最佳 Savonius 型风力机输出效率及转矩效率随叶尖速比的变化通过图6可以拟合曲线，得到风能利用系数Cp的函数：风机输出效率朋转矩做率C = f （九）=0.3656 x （-R ）2 + 0.6505 x （竺）（7）pVV由此可得输出功率为：P = f （九）x * pAV 3 = 0.3656 x
8、（罟）2 + 0.6505 x （罟）x | pAV 3（8）根据数学模型，进行Simulink仿真。输入为风速V，空气密度P和发电机 电磁转矩T_em；叶尖速比入和风力机受风面积A为常数。输出为风力机输出功 率P，输出转矩Te，输出转速3和发电机输入转速3_em。同时我们还设定了风 力机的启动风速和最高风速。模型如图 7 所示，并对其进行了封装。封装界面如 图 8。imdm in41 t mIsadir deniiivPioductlrt grjtorMlnMaxlDividedProd ud1)5吕q电u门|；UA3/IEatu，atifln DivideVO U1：_斤 nirid ip
9、eed Vo ij-i小-hl图 7 风力机仿真模型2K:OKParam et er smni Tjlade r： + io:CaiLcelHelpApplyblaie ar ma l；rr2):maz xTind 弓pwiGn/畧j:m i n vn. ml ip ee d (ti/1 sj :图 8 风力机封装界面设定输入风速为15m/s，空气密度为1.29kg/m3,仿真时间为10秒，得输出 曲线如图 9。横坐标为时间，纵坐标分别为功率、转矩、转速、发电机输入转速 风力机输出功率约为620W。图 9 风力机仿真输出曲线2.2 发电机2.2.1 永磁发电机永磁同步发电机由绕线转子同步发电机
10、发展而来，定子与普通同步发电机基 本相同，转子为永磁体，一般无阻尼绕组，因此不存在励磁绕组的铜损耗，同时 无需外部提供励磁电源，可以提高效率；转子上没有滑轮，可以提高系统的稳定 性。风力机输出的机械转矩带动发电机转子转动。永磁同步发电机的转子为永磁式结构，转子的磁链由永磁体决定。我们将定 子电压在dp0同步旋转坐标系下进行分解，其中，同步旋转坐标系的d轴是转子 磁链的方向。在此基础上建立发电机定子电压的 d 轴和 q 轴分量的表达式：u = R ida du =Riqa qd九.+ d CO Adt e qdA+ q + CO A dt e d9)式中：id和iq分别为发电机的d轴和q轴电流；
11、ud和uq分别为定子电压Eg的d轴和q轴分量;入d和入q分别是d轴和q轴的磁链Ra为定子电阻；a3e为电角频率定义磁链的 d 轴和 q 轴的分量的表达式为：式中：Ld和Lq分别为发电机的d轴和q轴电感；入。为永磁体产生的磁链定义q轴的反电势eq=3e入。，而d轴的反电势ed=0,因为发电机的转子为 对称结构，这里我们可以假设发电机的d轴和q轴的电感相等，即Ld=Lq=L。将 式（8）带入到式（7）并整理化简得到：di R .111)d = a I + CO I + U dt L d e q L d Jdi R ./.1 ） 、1q = a l CO （l + 九）+ U、dt L q edLo
12、 L q12)因为Ld=Lq=L，则永磁同步发电机的电磁转矩表达式为：T = 1.5（L L ）i i + i 九=1.5n i 九ed q d q q opq o式中n为极对数。pMATLAB 中提供的永磁发电机模型当其输入为负时，作为发电机使用。输 入有转矩Tm和转速3两种选择。本仿真软件中，我们选择转速3输入。参数设 置如图13所示，从上到下依次表示：定子相电阻，d轴和q轴定子电感，选择 机器常数，感应磁链，电压常数，转矩常数，极对数，初始条件（转速、角度、 电流）。通过调整感应磁链和极对数得到较为理想的曲线。图 14 和图 15 分别为 仿真模型和仿真结果。图 13 永磁同步发电机参数
13、设置图 15 发电机输出电压曲线2.2.2 电励磁同步发电机风力发电中所用的同步发电机绝大部分是三相同步电机，其输出联接到邻近 的三相电网或输配电线。普通三相同步发电机的原理结构如图 5 所示。在定子铁 心上有若干槽，槽内嵌有均匀分布的在空间彼此相隔 120电角的三相电枢绕组 aa、bb和cc，转子上装有磁极和励磁绕组，当励磁绕组通以直流电流If 后，电机内产生磁场。转子被风力机带动旋转，则磁场与定子三相绕组之间有相 对运动，从而在定子三相绕组中感应出三个幅值相同，彼此相隔120电角的交 流电势。这个交流电势的频率f决定于电机的极对数p和转子转速n,即f=pn/60。图 5 三相同步发电机结构
14、原理图我们将定子电压在dp0同步旋转坐标系下进行分解，则定子方程可表示为:d九+ ddt一九eq11)=R i +a0d九0dt u = R i +q 一九q a q dte d定义磁链的表达式：枕=Li + M i + M id d df fD Dv九=Li + M i/、q q q Q Q(12)九=L i0 0 0式中：Mf为定子绕组和励磁绕组之间的互感系数；MJ为定子绕组和直轴阻尼绕组之间的互感系数；Mq为定子绕组和交轴阻尼绕组之间的互感系数 发电机的电磁转矩为：T = 1.5n (九 i 九 i )ep d q q q(13)MATLAB 中的同步电机有国际标准单位和标幺值两种，我们
15、选择的是国际标准单位值的。输入分功率Pm和转速3两种，为了与永磁电机的输入统一，我们 还是选择转速3输入。参数设置如图16所示，参数依次为：视在功率、线电压、 频率、励磁电流，电子电阻、漏磁电感、d轴和q轴电感，定子的漏磁回路电阻、 漏磁电感，阻尼器的d轴和q轴电阻和漏磁电感，极对数，初始值(初始速度偏 差、电角度、线电流、相角、初始励磁电压)。Conf i eurat i on ; Pa rame t ers : Advaneed.ITominal rowsvoltaEe f reauency, f i eld current Pn(Vk) Vn(ViinE)ifn(Al :800 400 50Stater Rs (ohm) LL Lml. Lma
