风力发电机运行仿真.docx
《风力发电机运行仿真.docx》由会员分享,可在线阅读,更多相关《风力发电机运行仿真.docx(32页珍藏版)》请在咨信网上搜索。
1、基于 MATLAB 的“风力发电机运行仿真” 软件设计摘要关键词1 前言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)PC = p 1 p SV32式中:P 为风力机的功率,单位是W;为空气密度,单位是kg/m3;S 为风轮的扫风面积,单位是 m2;V
5、 为来流风速,单位是m/s在设计 Savonius 风力机时要考虑两个重要的结构参数:一个是重叠比 OL (Overlap ratio),一个是高径比 AP (Aapect ratio):(4)OL = S / d ()2AP = H / d ()3叶片重叠比对 Savonius 风力机的各种性能影响很大。如图 4 的风洞试验数 据所示, 具有不同的叶片重叠比的风力机的最大功率系数相差很大, 合理设计叶 片重叠比可以改善风力机的静态启动特性, 对风力机的动态力矩变化的战俘和相位也具有一定的影响。图 4 具有不同重叠比的 Savonius 风力机的性能叶片高径比也对风力机的性能影响很大, 一般来
6、说叶片高径比越大风力机性 能越好。 目前实际应用中的 Savonius 风力机的叶片高径比一般为 14, 准确数值 要根据设计目标、成本和安装地点的风况特点来决定。叶尖速比 是叶片的叶尖圆周速度与风速之比, 用来描述风轮在不同风速中的状态:入 = 2Rn = R V V式中:n 为风轮转速,单位是r/s; 为风轮角频率,单位是rad/s;R 为风轮半径,单位是m;V 为上游风速,单位是m/s风力机通过叶片捕获风能, 将风能转换为作用的发电机转子上的机械能, 将 吸收的叶片转矩为作用在发电机转子上的机械转矩。风力机吸收功率可以表示为风速的函数,其模型表示为:1(5)P = p C AV 32 p
7、所以风力机的机械转矩为:P 1 V 2T = = pC R3 2 p 入(6)其中相关参数的设定会影响风力机的输出效率。对于风力机建模,主要有两种方式,一种是对发电机的实测数据的查表法, 另一种是根据相关的公式进行垂直轴风力发电机输出参数建模。风能利用系数pC 的函数曲线如图 6 所示,由图可知,当叶尖速比在 1 左右时,输出效率最大。为了简化模型,我们在仿真过程中设定叶尖速比在速度范围之内为恒定值。图 6 最佳 Savonius 型风力机输出效率及转矩效率随叶尖速比的变化通过图 6 可以拟合曲线,得到风能利用系数 C 的函数:C = f (入) = 一0.3 6 5 ( R)2 + 0.6
8、5 0 (pR) ()7p V V由此可得输出功率为:P = f (入) 1 pAV 3 = 一0.3656 (R)2 + 0.6505 (R) 1 pAV 3 (8)2 V V 2根据数学模型,进行 Simulink 仿真。输入为风速 V,空气密度 和发电机 电磁转矩 T_em;叶尖速比 和风力机受风面积 A 为常数。输出为风力机输出功 率 P,输出转矩 Te,输出转速 和发电机输入转速 _em。同时我们还设定了风 力机的启动风速和最高风速。 模型如图 7 所示, 并对其进行了封装。 封装界面如图 8。图 7 风力机仿真模型图 8 风力机封装界面设定输入风速为 15m/s,空气密度为 1.2
9、9kg/m3,仿真时间为 10 秒,得输出 曲线如图 9。 横坐标为时间, 纵坐标分别为功率、 转矩、 转速、 发电机输入转速。 风力机输出功率约为 620W。图 9 风力机仿真输出曲线2.2 发电机2.2.1 永磁发电机永磁同步发电机由绕线转子同步发电机发展而来, 定子与普通同步发电机基 本相同,转子为永磁体,一般无阻尼绕组,因此不存在励磁绕组的铜损耗,同时 无需外部提供励磁电源, 可以提高效率; 转子上没有滑轮, 可以提高系统的稳定 性。风力机输出的机械转矩带动发电机转子转动。永磁同步发电机的转子为永磁式结构, 转子的磁链由永磁体决定。 我们将定 子电压在 dp0 同步旋转坐标系下进行分解
10、, 其中, 同步旋转坐标系的 d 轴是转子 磁链的方向。在此基础上建立发电机定子电压的 d 轴和 q 轴分量的表达式: (9)ud = Ra id + dt(d入)d - e入quq = Ra iq + d入qdt + e入d式中:id 和 iq 分别为发电机的 d 轴和 q 轴电流;ud 和 uq 分别为定子电压 Eg 的 d 轴和 q 轴分量; d 和 q 分别是 d 轴和 q 轴的磁链a 为电角频率R 为定子电阻;e定义磁链的 d 轴和 q 轴的分量的表达式为:(入 = L i + 入入q(d) = Lq(d)i q(d) o (10)式中:Ld 和 Lq 分别为发电机的 d 轴和 q
11、轴电感;o 为永磁体产生的磁链定义 q 轴的反电势 eq= e o ,而 d 轴的反电势 ed=0,因为发电机的转子为 对称结构,这里我们可以假设发电机的 d 轴和 q 轴的电感相等,即Ld=Lq=L。将 式(8)带入到式(7)并整理化简得到:d(di)t(d) = - id + oe iq + L(1) ud = - iq - oe (id + L(1) 入o ) + L(1) uq因为 Ld=Lq=L,则永磁同步发电机的电磁转矩表达式为:T = 1.5(L - L )i i + i 入 = 1.5n i 入e d q d q q o p q o(11)(12)式中 n 为极对数。pMATL
12、AB 中提供的永磁发电机模型当其输入为负时,作为发电机使用。输 入有转矩 Tm 和转速 两种选择。本仿真软件中,我们选择转速 输入。参数设 置如图 13 所示,从上到下依次表示:定子相电阻, d 轴和 q 轴定子电感,选择 机器常数,感应磁链,电压常数,转矩常数,极对数,初始条件(转速、角度、 电流)。通过调整感应磁链和极对数得到较为理想的曲线。图 14 和图 15 分别为仿真模型和仿真结果。图 13 永磁同步发电机参数设置图 14 永磁同步发电机模型图 15 发电机输出电压曲线2.2.2 电励磁同步发电机l入 = L i入 = L i + M , i | q q q Q Q风力发电中所用的同
13、步发电机绝大部分是三相同步电机, 其输出联接到邻近 的三相电网或输配电线。 普通三相同步发电机的原理结构如图 5 所示。 在定子铁 心上有若干槽, 槽内嵌有均匀分布的在空间彼此相隔 120电角的三相电枢绕组 aa、 bb和 cc,转子上装有磁极和励磁绕组,当励磁绕组通以直流电流 If 后, 电机内产生磁场。 转子被风力机带动旋转, 则磁场与定子三相绕组之间有相 对运动,从而在定子三相绕组中感应出三个幅值相同,彼此相隔120电角的交流电势。这个交流电势的频率 f 决定于电机的极对数 p 和转子转速n,即 f=pn/60。图 5 三相同步发电机结构原理图我们将定子电压在 dp0 同步旋转坐标系下进
14、行分解,则定子方程可表示为:u id + dt(d入)d - Oe入q|u = R i + d入q - O 入 | q a q dt e d u0 = Ra i0 + d入0dt定义磁链的表达式:(入 = L i + M , i + M , i| d d d f f D D0 0 0式中:fdM 为定子绕组和励磁绕组之间的互感系数;M 为定子绕组和直轴阻尼绕组之间的互感系数;M 为定子绕组和交轴阻尼绕组之间的互感系数Q(11)(12)发电机的电磁转矩为:T = 1.5n (入 i - 入 i )e p d q q q(13)MATLAB 中的同步电机有国际标准单位和标幺值两种,我们选择的是国际
15、标准单位值的。输入分功率 Pm 和转速两种,为了与永磁电机的输入统一,我们还是选择转速输入。参数设置如图16 所示,参数依次为:视在功率、线电压、 频率、 励磁电流, 电子电阻、 漏磁电感、 d 轴和 q 轴电感, 定子的漏磁回路电阻、 漏磁电感,阻尼器的 d 轴和 q 轴电阻和漏磁电感,极对数,初始值(初始速度偏 差、电角度、线电流、相角、初始励磁电压)。图 16 电励磁同步发电机参数设置这一发电机模块需要配合励磁模块一起使用,在这里需要说明一下。图 17 为励磁模块的基础模型,从左到右依次是:低通滤波器,超前滞后补偿,主调节 器,饱和度,励磁模型,阻尼器。通过改变主调节器的增益、时间常数和
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 风力发电机 运行 仿真
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【鱼**】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【鱼**】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。