风力发电机动力学仿真研究
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
。．拦王三Ｅ三丑 ＴｕｒｂＳｉｍ程序生成全域湍流风模型，用ＩＥＣＫＡＩ
萋言． 一｛１ Ｅ王三Ｅ王１－丑１１
姿＃＝王三二＝Ｆ＝；习 一ｌＬ———＿－＿ｊ———＿—＿一Ｌ————－—上——————Ｌ＿————Ｊ
＿１Ｌ——玉广—矗忑—１赉广■袁矿—喇．０
叶片径向位置（ｍ）
Ｆｉｇ．２
圈２挥舞方向五阶振型 Ｆｌａｐｗｉｓｅ ｆｉｖｅ ｍｏｄｅ ｓｈａｐｅｓ ｏｆ ｔｈｅ ｂｌａｄｅ
Ａｂｓｔｒａｃｔ：Ｔｈｅ ＦＡＳＴ ｃｏｄｅ ｗａｓ Ｉｎｔｒｏｄｕｃｅｄ ａｎｄ ｅｍｐｌｏｙｅｄ ｔＯ ｓｉｍｕｌａｔｅ ｔｈｅ ｄｙｎａｍｉｃ ｒｅｓｐｏｎｓｅ ｏｆ ｗｉｎｄ ｔｕｒ— ｂｉｎｅ；ｔｈｅ ｗｉｎｄ ｔｕｒｂｉｎｅ ｍｏｄｅｌ ｗａｓ ｂｕｉｌｔ ｂｙ ｔｈｅ ｓｕｂｒｏｕｔｉｎｅ ＴｕｒｂＳｉｍ．Ａｆｔｅｒ ｌｏａｄｉｎｇ ｔｈｅ ａｅｒｏｄｙｎａｍｉｃｓ ｌｏａｄｓ ｃａｌｃｕｌａｔｅｄ ｆｒｏｍ ＡｅｒｏＤｙｎ ｏｎ ｔｈｅ ｗｉｎｄ ｔｕｒｂｉｎｅ ｍｏｄｅｌ，ｔｈｅ ｄｙｎａｍｉｃ ｒｅｓｐｏｎｓｅ ｏｆ ｗｉｎｄ ｔｕｒｂｉｎｅ ｗａｓ ａｃｈｉｅｖｅｄ． Ｔｈｅ ＡＤＡＭＳ ｗｉｎｄ ｔｕｒｂｉｎｅ ｍｏｄｅｌ ｗａｓ ｃｒｅａｔｅｄ ｂｙ ＦＡＳＴ ｃｏｄｅ．Ｔｈｅ ｒｅｓｐｏｎｓｅ ａｎａｌｙｓｉｓ ０ｎ ｄｙｎａｍｉｃ ｓｉｍｕｌａｔｉｏｎ ｏｆ ｗｉｎｄ ｔｕｒｂｉｎｅ ｗｅｒｅ ｄｉｓｃｕｓｓｅｄ ｂｙ ＦＡＳＴ ｃｏｄｅ． Ｋｅｙ ｗｏｒｄｓ： ＦＡＳＴ；ｗｉｎｄ ｔｕｒｂｉｎｅ；ｄｙｎａｍｉｃ ｒｅｓｐｏｎｓｅ
应影响了其余叶片的叶尖位移响应。 由图９知，通过ＦＡＳＴ程序，叶片１的叶尖在
挥舞方向上的加速度响应较大，在摆振方向上却影 响了其余叶片叶尖的加速度响应，三支叶片叶尖摆 振方向的加速度响应总体增加，叶片２和叶片３的 增加量大于叶片１，可得叶片２和叶片３受叶片１ 的快速变桨在摆振方向上的影响较大。
叶片１挥舞方向的弯矩（Ｒｏｏｔ Ｍｙｂｌ）和摆振方 向的弯矩（Ｒｏｏｔ Ｍｘｂｌ）见图１０，可知摆振方向的弯 矩值变化不大，而挥舞方向的弯矩值急剧变化，这是 因为叶片变桨后，气动推力减小（图１１）。同时叶片 １叶根的变桨力矩也明显减小（图１２）。
本研究采用ＦＡＳＴ程序进行风力发电机的动 力学仿真分析，通过子程序ＴｕｒｂＳｉｍ生成风文件， 调用子程序ＡｅｒｏＤｙｎ计算气动载荷并加载到风力 发电机的仿真模型上，实现风力发电机的动力学响 应计算。
序中包含叶片参数，风模型及采用的计算模型。 ＴｕｒｂＳｉｍ程序是用来生成全域湍流风模型。Ｍｏｄｅｓ 程序是用来计算柔性体的固有频率和振型，为 ＦＡＳＴ输入参数作准备。
但上述仿真研究的重点是考察叶轮的气动弹性
收稿日期：２００９—１２—２０ 基金项目：新疆金风科技股份有限公司项目（２００６ＢＡＡ０１Ａ０２） 通讯作者：崔新维，Ｅ－ｍａｉｌ：ｃｕｉｘｉｎｗｅｉ＠１２６．ｃｏｒｎ
万方数据
１７８
新疆农业大学学报
２０１０年
稳定性问题，并没有进行风力发电机的动力学仿真 分析。目前国内对风力发电机性能的研究主要借助 于ＧＨ Ｂｌａｄｅｄ软件来完成。
湍流风模型，ＩＥＣ规定的Ａ级湍流强度，轮毂高度 平均风速为１２ ｍ／ｓ，分析时间为９．９～１９．９ ｓ，生成 模型见图４和图５。
襄１ 叶片挥舞方向和摆振方向前五阶固有频率
Ｔａｂｌｅ １ Ｆｉｖｅ ｆｒｅｑｕｅｎｃｉｅｓ ｏｆ ｆｌａｐｗｉｓｅ ａｎｄ ｅｄｇｅｗｉｓｅ ｏｆ ｂｌａｄｅ Ｈｚ
ｏ．ｊＥ王三巨王丑
新疆农业大学学报２０１０，３３（２）：１７７＂－１８１ Ｊｏｕｒｎａｌ ｏ｜Ｘｉｎｊｉａｎｇ Ａｇｒｉｃｕｌｔｕｒａｌ Ｕｎｉｖｅｒｓｉｔｙ
文章编号：１００７—８６１４（２０１０）０２—０１７７—０５
风力发电机动力学仿真研究
付长江，崔新维
（新疆农业大学机械交通学院。鸟鲁木齐８３００５２）
摘要： 采用ＦＡＳＴ程序进行风力发电机的动力学仿真分析，通过子程序ＴｕｒｂＳｉｍ生成湍流风文件，调用子程序 ＡｅｒｏＤｙｎ计算气动载荷并加载到风力发电机的仿真模型上，实现风力发电机的动力学响应计算。通过ＦＡＳＴ程序 建立了ＡＤＡＭＳ的风力发电机模型，对以ＦＡＳＴ实现动力学响应分析进行了探索性的研究。
ｉ
ｌ
一｛曼硝匿
１
１
ｌ
言１ 昌ｌ 硝 匿 姐
．
时间（ｓ） 圈４轮毂高度三个方向的风速 Ｆｉｇ．４ ３Ｄ ｗｉｎｄ ｓｐｅｅｄ ｏｆ ｔｈｅ ｈｕｂ ｈｅｉｇｈｔ
时间（ｓ）
Ｆ嘻５
圈５轮敦高度合风速
Ｔｏｔａｌ ｗｉｎｄ ｓｐｅｅｄ ｏｆ ｔｈｅ ｈｕｂ ｈｅｉｇｈｔ
２结果与分析
ＦＡＳＴ程序可以作风力发电机的动力学响应仿 真分析，可以仿真ＩＥＣ中规定的工况。风力发电机 的起动、正常顺桨停机、叶片顺桨失败停机、叶片突 然顺桨、停机、空转等工况。本研究所仿真的工况为 在１３ ｓ时叶片１停止正常的变桨，在１５ ｓ到达指定 角度。通过ＦＡＳＴ计算，叶尖出平面位移量的瞬时 响应见图６，叶尖平面内位移量的瞬时响应见图７。 从分析结果来看，在此工况下，叶片１叶尖出平面和 平面内的位移响应方向较其余叶片总体上不同，叶 尖位移响应量大，从图６可知，叶片１的叶尖位移响
风力发电机组系统是一个非常复杂的强非线性 流一刚一柔耦合的周期时变多体系统。建立精确的动 态模型是困难的。由于风载荷具有交变性和随机 性，叶片作为弹性结构发生振动是必然的。在正常 工作时，叶片在绕转轴作大范围的空间旋转运动时， 会产生出旋转平面的挥舞振动、旋转平面内的摆振 振动及绕叶片轴线的扭转振动。叶片振动引起的非 周期动载荷通过轮毂直接作用于机组系统。而机组 系统本身就是一个复杂的耦合非线性振动系统，影 响系统振动因素较多，包括传动轴的柔性以及电机 特性等。同时轴系在运转过程中将轴系动载荷通过 轮毂反作用于叶片，影响叶片的运动形态，进而又影 响叶片的受载。
ＡｅｒｏＤｙｎ程序是计算气动载荷的子程序。程
图１ ＦＡＳＴ程序输入、输出文件 Ｆｉｇ．１ ＦＡＳＴ ｉｎｐｕｔ ａｎｄ ｏｕｔｐｕｔ ｆｉｌｅｓ
１．２．２ 风力发电机的ＡＤＡＭＳ模型 ＦＡＳＴ可以生成ＡＤＡＭＳ模型数据，导人数据
即可生成风力发电机的模型。ＦＡＳＴ的功能类似于 ＡＤＡＭＳ的前处理器，此前处理器可以通过简单的 指定ＦＡＳＴ输入文件来建立模型，方便和简化了使 用ＡＤＡＭＳ建立风力发电机模型。但ＦＡＳＴ只能 建立中等复杂的风力发电机模型，一旦ＦＡＳＴ的风 力发电机模型建成，只需再添加一些工作，就可建立 复杂的ＡＤＡＭＳ风力发电机模型。从ＦＡＳＴ传递 到ＡＤＡＭＳ的风力发电机模型包含了所有ＦＡＳＴ 中的参数ｉ也夹带了ＦＡＳＴ中的一些缺陷，如叶片 和塔架的模态截断近似。而许多ＦＡＳＴ模型中无 法计算的量，在ＡＤＡＭＳ模型中仿真却是常规的， 如叶片和塔架的扭转和拉伸自由度、预弯叶片、叶片 的质量和弹性中心的偏置、塔架质量中心的偏置、变 桨电机的动力学响应分析及ＡＤＡＭＳ拥有强大的 图形处理能力。 １．３叶片模型
万方数据
第２期
付长江，等：风力发电机动力学仿真研究
１７９
叶轮转速和桨距角也影响叶片模态。在发生振动 时，叶片的低阶模态响应的贡献率要高于高阶模态 响应。Ｍｏｄｅｓ程序计算叶片模态时将叶片分成柔 性段和刚性段，不仅考虑以上的一些因子，同时也把 叶尖质量也作为一个考虑因子。通过程序Ｍｏｄｅｓ 计算叶片挥舞方向和摆振方向的前五阶固有频率见 表１。其振型见图２和图３。 １．４风模型
ＦＡＳＴ采用线性模态方法来描述风力发电机组 的柔性单元，如塔筒和叶片。ＦＡＳＴ使用两阶的挥 舞模态和一阶的摆振模态，模态受叶片的结构影响。 在建立变速风力发电机叶片模态模型时，选择一个 转速设计点下求解叶片模态振型。叶轮转速对叶片 的模态振型的影响没有对振动频率的影响严重。
影响叶片模态的因素有叶根的约束形式，单位 长度质量，挥舞方向和摆振方向的刚度。柔性较大 的叶片在空间做大范围的运动，运动产生的变形使 刚度发生变化，轮毂对叶片模态也存在影响。同时
关键词：ＦＡＳＴ；风力wenku.baidu.com电视ｉ动力学响应
中图分类号：ＴＫ８３
文献标识码：Ａ
Ｓｔｕｄｙ ｏｎ Ｄｙｎａｍｉｃ Ｓｉｍｕｌａｔｉｏｎ ｏｆ Ｗｉｎｄ Ｔｕｒｂｉｎｅ
Ｆｕ Ｃｈａｎｇ－ｊｉａｎｇ．ＣＵＩ Ｘｉｎ－ｗｅｉ
（Ｃｏｌｌｅｇｅ ｏｆ Ｍｅｃｈａｎｉｃａｌ Ｅｎｇｉｎｅｅｒｉｎｇ ａｎｄ Ｃｏｍｍｕｎｉｃａｔｉｏｎ，Ｘｉｎｊｉａｎｇ Ａｇｒｉｃｕｌｔｕｒａｌ Ｕｎｉｖｅｒｓｉｔｙ，Ｕｒｕ－ ｍｑｉ ８３００５２，Ｃｈｉｎａ）
１ 模型构建
１．１风力发电机模型 风力发电机模型包括风模型、空气动力模型、结
构动力学模型和控制模型。空气动力模型计算叶轮 气动性能和气动载荷，结构动力学模型包括叶轮、塔 架和传动系统的结构动力学：控制模型包括变桨控 制、偏航控制、发电机变速控制。空气动力模型根据 风模型输入计算出叶片载荷作用到结构动力学模型 上，同时结构动力学模型又将结构的变形反馈到空 气动力模型上；控制模型获取当前各控制参数信息， 根据来流情况及其他需要，对机组进行控制，使风力 发电机最大限度捕获风能［６］。本研究采用ＮＲＥＬ 的ＷｉｎｄＰＡＣＴｌ．５ ＭＷ三叶片变速风力发电机模 型，叶轮直径７０ ｍ，叶片是ＮＲＥＬ采用Ｓ系列翼型 设计，长度为３３．２５ ｍ。塔筒高度８２．３９ ｍ，机舱质 量５ １１７０ ｋｇ，叶轮质量１５ １４８ ｋｇ。 １．２采用ＦＡＳＴ建立风力发电机模型
言一÷Ｅ壬，—ｔ—＇同ｉ－－ｉｌｉ ｉ
羲一；；［ｉ ＩＩ二［ｉ ＩＩ
Ｉ
１—Ｉ－！，Ｉ—ｉ 萋一ｊ］目ｉ ｉ干ｉ≠ｉ嗣ｉ ｉ寻ｉ干ｉ；ｉ翮ｉ ｌ
ｉ一／ｉ０
一ｌＬ——￡广—－ｒ药『一—嚣万—丐蠡广吲．０
叶片径向位置（ｍ）
Ｈ￡．３
图３摆振方向五阶振型 ＦＡ２ｅｗｉｓｅ ｆｉｖｅ ｄｉｒｅｃｔｉｏｎ ｓｈａｐｅｓ ｏｆ ｔｈｅ ｂｌａｄｅ
因此，建立风力发电机模型时需要考虑如何简 化及求解。仿真技术广泛吸收了数值计算、力学等
学科的理论基础，借助于计算机技术而发展起来的 学科。利用这项技术，可对许多工程中的实际问题 进行数值建模仿真。Ｗｅｎｄｅｈ推导了适用于风轮桨 叶的气动载荷，用非耦合的非旋转模态研究了其气 动弹性稳定性问题［１］。Ｃｈｏｐｒａ［２］用非线性半刚性模 型研究了桨叶气动弹性响应和稳定性问题。李立本 等［３］建立风力发电机转子叶片的非线性运动方程， 采用模态法求解挥舞、摆振、扭转微分方程并应用了 数值结果对风力发电机的气动弹性稳定性进行了分 析。李德源和叶枝全等¨’５１对旋转叶片进行有限元 离散，采用凯恩（Ｋａｎｅ）方法，建立了一般形式的大 型风力机叶片柔性多体动力学方程，导出了旋转叶 片的空间梁单元动力矩阵的显式形式。计算了叶片 的固有频率。
ＦＡＳＴ程序的输入、输出文件总体结构示意见 图１。ＦＡＳＴ用主程序描述风力发电机运行参数和 基本几何尺寸参数。这些参数包括仿真控制、风力 发电机控制、重力环境条件、自由度选择、风力发电 机初始条件、风力发电机配置、各部件的质量和转动 惯量、传动系、发电机模型、基础模型、塔架模型、机 头偏航动力学参数、叶片模型、气动模型、ＡＤＡＭＳ 数据接口、线性化控制和输出参数说明。
万方数据
１８０
罾
嘲 静 ｇ 暄 斗 茁
；～ 新疆农业大学学报 ．一叶片 ——叶片： ——叶片：
２０１０年
吕
咖
簿
趟 墨
｜
喧
斗
时间（ａ）
圈６叶尖出平面的位移响应 Ｆｉｇ．６ Ｄｉｓｐｌａｃｅｍｅｎｔ ｒｅｓｐｏｎｅ ｏｕｔ－ｏｆ－ｐｌａｎｅ ｂｌａｄｅ ｔｉｐ
ＦＡＳＴ（Ｆａｔｉｇｕｅ，Ａｅｒｏｄｙｎａｍｉｃｓ，Ｓｔｒｕｃｔｕｒｅｓ，
Ｔｕｒｂｕｌｅｎｃｅ）程序现由ＮＲＥＬ开发和维护，是一个 复杂的气弹仿真器，具备计算水平轴两叶片、三叶片 风力发电机的极限载荷和疲劳载荷［７］。ＦＡＳＴ程序 由两、三叶片水平轴风力发电机和气动子程序 ＡｅｒｏＤｙｎ组成。同时开发了与ＧＨ Ｂｌａｄｅｄ的数据 接口［８］。２００５年，ＦＡＳＴ与ＡｅｒｏＤｙｎ通过Ｇｅｒｍａ－ ｎｉｓｃｈｅｒ Ｌｌｏｙｄ Ｗｉｎｄ Ｅｎｅｒｇｉｅ的评估，认为适合“为 设计和认证陆上风力发电机计算载荷”［９１。 １．２．１ ＦＡＳＴ程序输入、输出文件说明
萋言． 一｛１ Ｅ王三Ｅ王１－丑１１
姿＃＝王三二＝Ｆ＝；习 一ｌＬ———＿－＿ｊ———＿—＿一Ｌ————－—上——————Ｌ＿————Ｊ
＿１Ｌ——玉广—矗忑—１赉广■袁矿—喇．０
叶片径向位置（ｍ）
Ｆｉｇ．２
圈２挥舞方向五阶振型 Ｆｌａｐｗｉｓｅ ｆｉｖｅ ｍｏｄｅ ｓｈａｐｅｓ ｏｆ ｔｈｅ ｂｌａｄｅ
Ａｂｓｔｒａｃｔ：Ｔｈｅ ＦＡＳＴ ｃｏｄｅ ｗａｓ Ｉｎｔｒｏｄｕｃｅｄ ａｎｄ ｅｍｐｌｏｙｅｄ ｔＯ ｓｉｍｕｌａｔｅ ｔｈｅ ｄｙｎａｍｉｃ ｒｅｓｐｏｎｓｅ ｏｆ ｗｉｎｄ ｔｕｒ— ｂｉｎｅ；ｔｈｅ ｗｉｎｄ ｔｕｒｂｉｎｅ ｍｏｄｅｌ ｗａｓ ｂｕｉｌｔ ｂｙ ｔｈｅ ｓｕｂｒｏｕｔｉｎｅ ＴｕｒｂＳｉｍ．Ａｆｔｅｒ ｌｏａｄｉｎｇ ｔｈｅ ａｅｒｏｄｙｎａｍｉｃｓ ｌｏａｄｓ ｃａｌｃｕｌａｔｅｄ ｆｒｏｍ ＡｅｒｏＤｙｎ ｏｎ ｔｈｅ ｗｉｎｄ ｔｕｒｂｉｎｅ ｍｏｄｅｌ，ｔｈｅ ｄｙｎａｍｉｃ ｒｅｓｐｏｎｓｅ ｏｆ ｗｉｎｄ ｔｕｒｂｉｎｅ ｗａｓ ａｃｈｉｅｖｅｄ． Ｔｈｅ ＡＤＡＭＳ ｗｉｎｄ ｔｕｒｂｉｎｅ ｍｏｄｅｌ ｗａｓ ｃｒｅａｔｅｄ ｂｙ ＦＡＳＴ ｃｏｄｅ．Ｔｈｅ ｒｅｓｐｏｎｓｅ ａｎａｌｙｓｉｓ ０ｎ ｄｙｎａｍｉｃ ｓｉｍｕｌａｔｉｏｎ ｏｆ ｗｉｎｄ ｔｕｒｂｉｎｅ ｗｅｒｅ ｄｉｓｃｕｓｓｅｄ ｂｙ ＦＡＳＴ ｃｏｄｅ． Ｋｅｙ ｗｏｒｄｓ： ＦＡＳＴ；ｗｉｎｄ ｔｕｒｂｉｎｅ；ｄｙｎａｍｉｃ ｒｅｓｐｏｎｓｅ
应影响了其余叶片的叶尖位移响应。 由图９知，通过ＦＡＳＴ程序，叶片１的叶尖在
挥舞方向上的加速度响应较大，在摆振方向上却影 响了其余叶片叶尖的加速度响应，三支叶片叶尖摆 振方向的加速度响应总体增加，叶片２和叶片３的 增加量大于叶片１，可得叶片２和叶片３受叶片１ 的快速变桨在摆振方向上的影响较大。
叶片１挥舞方向的弯矩（Ｒｏｏｔ Ｍｙｂｌ）和摆振方 向的弯矩（Ｒｏｏｔ Ｍｘｂｌ）见图１０，可知摆振方向的弯 矩值变化不大，而挥舞方向的弯矩值急剧变化，这是 因为叶片变桨后，气动推力减小（图１１）。同时叶片 １叶根的变桨力矩也明显减小（图１２）。
本研究采用ＦＡＳＴ程序进行风力发电机的动 力学仿真分析，通过子程序ＴｕｒｂＳｉｍ生成风文件， 调用子程序ＡｅｒｏＤｙｎ计算气动载荷并加载到风力 发电机的仿真模型上，实现风力发电机的动力学响 应计算。
序中包含叶片参数，风模型及采用的计算模型。 ＴｕｒｂＳｉｍ程序是用来生成全域湍流风模型。Ｍｏｄｅｓ 程序是用来计算柔性体的固有频率和振型，为 ＦＡＳＴ输入参数作准备。
但上述仿真研究的重点是考察叶轮的气动弹性
收稿日期：２００９—１２—２０ 基金项目：新疆金风科技股份有限公司项目（２００６ＢＡＡ０１Ａ０２） 通讯作者：崔新维，Ｅ－ｍａｉｌ：ｃｕｉｘｉｎｗｅｉ＠１２６．ｃｏｒｎ
万方数据
１７８
新疆农业大学学报
２０１０年
稳定性问题，并没有进行风力发电机的动力学仿真 分析。目前国内对风力发电机性能的研究主要借助 于ＧＨ Ｂｌａｄｅｄ软件来完成。
湍流风模型，ＩＥＣ规定的Ａ级湍流强度，轮毂高度 平均风速为１２ ｍ／ｓ，分析时间为９．９～１９．９ ｓ，生成 模型见图４和图５。
襄１ 叶片挥舞方向和摆振方向前五阶固有频率
Ｔａｂｌｅ １ Ｆｉｖｅ ｆｒｅｑｕｅｎｃｉｅｓ ｏｆ ｆｌａｐｗｉｓｅ ａｎｄ ｅｄｇｅｗｉｓｅ ｏｆ ｂｌａｄｅ Ｈｚ
ｏ．ｊＥ王三巨王丑
新疆农业大学学报２０１０，３３（２）：１７７＂－１８１ Ｊｏｕｒｎａｌ ｏ｜Ｘｉｎｊｉａｎｇ Ａｇｒｉｃｕｌｔｕｒａｌ Ｕｎｉｖｅｒｓｉｔｙ
文章编号：１００７—８６１４（２０１０）０２—０１７７—０５
风力发电机动力学仿真研究
付长江，崔新维
（新疆农业大学机械交通学院。鸟鲁木齐８３００５２）
摘要： 采用ＦＡＳＴ程序进行风力发电机的动力学仿真分析，通过子程序ＴｕｒｂＳｉｍ生成湍流风文件，调用子程序 ＡｅｒｏＤｙｎ计算气动载荷并加载到风力发电机的仿真模型上，实现风力发电机的动力学响应计算。通过ＦＡＳＴ程序 建立了ＡＤＡＭＳ的风力发电机模型，对以ＦＡＳＴ实现动力学响应分析进行了探索性的研究。
ｉ
ｌ
一｛曼硝匿
１
１
ｌ
言１ 昌ｌ 硝 匿 姐
．
时间（ｓ） 圈４轮毂高度三个方向的风速 Ｆｉｇ．４ ３Ｄ ｗｉｎｄ ｓｐｅｅｄ ｏｆ ｔｈｅ ｈｕｂ ｈｅｉｇｈｔ
时间（ｓ）
Ｆ嘻５
圈５轮敦高度合风速
Ｔｏｔａｌ ｗｉｎｄ ｓｐｅｅｄ ｏｆ ｔｈｅ ｈｕｂ ｈｅｉｇｈｔ
２结果与分析
ＦＡＳＴ程序可以作风力发电机的动力学响应仿 真分析，可以仿真ＩＥＣ中规定的工况。风力发电机 的起动、正常顺桨停机、叶片顺桨失败停机、叶片突 然顺桨、停机、空转等工况。本研究所仿真的工况为 在１３ ｓ时叶片１停止正常的变桨，在１５ ｓ到达指定 角度。通过ＦＡＳＴ计算，叶尖出平面位移量的瞬时 响应见图６，叶尖平面内位移量的瞬时响应见图７。 从分析结果来看，在此工况下，叶片１叶尖出平面和 平面内的位移响应方向较其余叶片总体上不同，叶 尖位移响应量大，从图６可知，叶片１的叶尖位移响
风力发电机组系统是一个非常复杂的强非线性 流一刚一柔耦合的周期时变多体系统。建立精确的动 态模型是困难的。由于风载荷具有交变性和随机 性，叶片作为弹性结构发生振动是必然的。在正常 工作时，叶片在绕转轴作大范围的空间旋转运动时， 会产生出旋转平面的挥舞振动、旋转平面内的摆振 振动及绕叶片轴线的扭转振动。叶片振动引起的非 周期动载荷通过轮毂直接作用于机组系统。而机组 系统本身就是一个复杂的耦合非线性振动系统，影 响系统振动因素较多，包括传动轴的柔性以及电机 特性等。同时轴系在运转过程中将轴系动载荷通过 轮毂反作用于叶片，影响叶片的运动形态，进而又影 响叶片的受载。
ＡｅｒｏＤｙｎ程序是计算气动载荷的子程序。程
图１ ＦＡＳＴ程序输入、输出文件 Ｆｉｇ．１ ＦＡＳＴ ｉｎｐｕｔ ａｎｄ ｏｕｔｐｕｔ ｆｉｌｅｓ
１．２．２ 风力发电机的ＡＤＡＭＳ模型 ＦＡＳＴ可以生成ＡＤＡＭＳ模型数据，导人数据
即可生成风力发电机的模型。ＦＡＳＴ的功能类似于 ＡＤＡＭＳ的前处理器，此前处理器可以通过简单的 指定ＦＡＳＴ输入文件来建立模型，方便和简化了使 用ＡＤＡＭＳ建立风力发电机模型。但ＦＡＳＴ只能 建立中等复杂的风力发电机模型，一旦ＦＡＳＴ的风 力发电机模型建成，只需再添加一些工作，就可建立 复杂的ＡＤＡＭＳ风力发电机模型。从ＦＡＳＴ传递 到ＡＤＡＭＳ的风力发电机模型包含了所有ＦＡＳＴ 中的参数ｉ也夹带了ＦＡＳＴ中的一些缺陷，如叶片 和塔架的模态截断近似。而许多ＦＡＳＴ模型中无 法计算的量，在ＡＤＡＭＳ模型中仿真却是常规的， 如叶片和塔架的扭转和拉伸自由度、预弯叶片、叶片 的质量和弹性中心的偏置、塔架质量中心的偏置、变 桨电机的动力学响应分析及ＡＤＡＭＳ拥有强大的 图形处理能力。 １．３叶片模型
万方数据
第２期
付长江，等：风力发电机动力学仿真研究
１７９
叶轮转速和桨距角也影响叶片模态。在发生振动 时，叶片的低阶模态响应的贡献率要高于高阶模态 响应。Ｍｏｄｅｓ程序计算叶片模态时将叶片分成柔 性段和刚性段，不仅考虑以上的一些因子，同时也把 叶尖质量也作为一个考虑因子。通过程序Ｍｏｄｅｓ 计算叶片挥舞方向和摆振方向的前五阶固有频率见 表１。其振型见图２和图３。 １．４风模型
ＦＡＳＴ采用线性模态方法来描述风力发电机组 的柔性单元，如塔筒和叶片。ＦＡＳＴ使用两阶的挥 舞模态和一阶的摆振模态，模态受叶片的结构影响。 在建立变速风力发电机叶片模态模型时，选择一个 转速设计点下求解叶片模态振型。叶轮转速对叶片 的模态振型的影响没有对振动频率的影响严重。
影响叶片模态的因素有叶根的约束形式，单位 长度质量，挥舞方向和摆振方向的刚度。柔性较大 的叶片在空间做大范围的运动，运动产生的变形使 刚度发生变化，轮毂对叶片模态也存在影响。同时
关键词：ＦＡＳＴ；风力wenku.baidu.com电视ｉ动力学响应
中图分类号：ＴＫ８３
文献标识码：Ａ
Ｓｔｕｄｙ ｏｎ Ｄｙｎａｍｉｃ Ｓｉｍｕｌａｔｉｏｎ ｏｆ Ｗｉｎｄ Ｔｕｒｂｉｎｅ
Ｆｕ Ｃｈａｎｇ－ｊｉａｎｇ．ＣＵＩ Ｘｉｎ－ｗｅｉ
（Ｃｏｌｌｅｇｅ ｏｆ Ｍｅｃｈａｎｉｃａｌ Ｅｎｇｉｎｅｅｒｉｎｇ ａｎｄ Ｃｏｍｍｕｎｉｃａｔｉｏｎ，Ｘｉｎｊｉａｎｇ Ａｇｒｉｃｕｌｔｕｒａｌ Ｕｎｉｖｅｒｓｉｔｙ，Ｕｒｕ－ ｍｑｉ ８３００５２，Ｃｈｉｎａ）
１ 模型构建
１．１风力发电机模型 风力发电机模型包括风模型、空气动力模型、结
构动力学模型和控制模型。空气动力模型计算叶轮 气动性能和气动载荷，结构动力学模型包括叶轮、塔 架和传动系统的结构动力学：控制模型包括变桨控 制、偏航控制、发电机变速控制。空气动力模型根据 风模型输入计算出叶片载荷作用到结构动力学模型 上，同时结构动力学模型又将结构的变形反馈到空 气动力模型上；控制模型获取当前各控制参数信息， 根据来流情况及其他需要，对机组进行控制，使风力 发电机最大限度捕获风能［６］。本研究采用ＮＲＥＬ 的ＷｉｎｄＰＡＣＴｌ．５ ＭＷ三叶片变速风力发电机模 型，叶轮直径７０ ｍ，叶片是ＮＲＥＬ采用Ｓ系列翼型 设计，长度为３３．２５ ｍ。塔筒高度８２．３９ ｍ，机舱质 量５ １１７０ ｋｇ，叶轮质量１５ １４８ ｋｇ。 １．２采用ＦＡＳＴ建立风力发电机模型
言一÷Ｅ壬，—ｔ—＇同ｉ－－ｉｌｉ ｉ
羲一；；［ｉ ＩＩ二［ｉ ＩＩ
Ｉ
１—Ｉ－！，Ｉ—ｉ 萋一ｊ］目ｉ ｉ干ｉ≠ｉ嗣ｉ ｉ寻ｉ干ｉ；ｉ翮ｉ ｌ
ｉ一／ｉ０
一ｌＬ——￡广—－ｒ药『一—嚣万—丐蠡广吲．０
叶片径向位置（ｍ）
Ｈ￡．３
图３摆振方向五阶振型 ＦＡ２ｅｗｉｓｅ ｆｉｖｅ ｄｉｒｅｃｔｉｏｎ ｓｈａｐｅｓ ｏｆ ｔｈｅ ｂｌａｄｅ
因此，建立风力发电机模型时需要考虑如何简 化及求解。仿真技术广泛吸收了数值计算、力学等
学科的理论基础，借助于计算机技术而发展起来的 学科。利用这项技术，可对许多工程中的实际问题 进行数值建模仿真。Ｗｅｎｄｅｈ推导了适用于风轮桨 叶的气动载荷，用非耦合的非旋转模态研究了其气 动弹性稳定性问题［１］。Ｃｈｏｐｒａ［２］用非线性半刚性模 型研究了桨叶气动弹性响应和稳定性问题。李立本 等［３］建立风力发电机转子叶片的非线性运动方程， 采用模态法求解挥舞、摆振、扭转微分方程并应用了 数值结果对风力发电机的气动弹性稳定性进行了分 析。李德源和叶枝全等¨’５１对旋转叶片进行有限元 离散，采用凯恩（Ｋａｎｅ）方法，建立了一般形式的大 型风力机叶片柔性多体动力学方程，导出了旋转叶 片的空间梁单元动力矩阵的显式形式。计算了叶片 的固有频率。
ＦＡＳＴ程序的输入、输出文件总体结构示意见 图１。ＦＡＳＴ用主程序描述风力发电机运行参数和 基本几何尺寸参数。这些参数包括仿真控制、风力 发电机控制、重力环境条件、自由度选择、风力发电 机初始条件、风力发电机配置、各部件的质量和转动 惯量、传动系、发电机模型、基础模型、塔架模型、机 头偏航动力学参数、叶片模型、气动模型、ＡＤＡＭＳ 数据接口、线性化控制和输出参数说明。
万方数据
１８０
罾
嘲 静 ｇ 暄 斗 茁
；～ 新疆农业大学学报 ．一叶片 ——叶片： ——叶片：
２０１０年
吕
咖
簿
趟 墨
｜
喧
斗
时间（ａ）
圈６叶尖出平面的位移响应 Ｆｉｇ．６ Ｄｉｓｐｌａｃｅｍｅｎｔ ｒｅｓｐｏｎｅ ｏｕｔ－ｏｆ－ｐｌａｎｅ ｂｌａｄｅ ｔｉｐ
ＦＡＳＴ（Ｆａｔｉｇｕｅ，Ａｅｒｏｄｙｎａｍｉｃｓ，Ｓｔｒｕｃｔｕｒｅｓ，
Ｔｕｒｂｕｌｅｎｃｅ）程序现由ＮＲＥＬ开发和维护，是一个 复杂的气弹仿真器，具备计算水平轴两叶片、三叶片 风力发电机的极限载荷和疲劳载荷［７］。ＦＡＳＴ程序 由两、三叶片水平轴风力发电机和气动子程序 ＡｅｒｏＤｙｎ组成。同时开发了与ＧＨ Ｂｌａｄｅｄ的数据 接口［８］。２００５年，ＦＡＳＴ与ＡｅｒｏＤｙｎ通过Ｇｅｒｍａ－ ｎｉｓｃｈｅｒ Ｌｌｏｙｄ Ｗｉｎｄ Ｅｎｅｒｇｉｅ的评估，认为适合“为 设计和认证陆上风力发电机计算载荷”［９１。 １．２．１ ＦＡＳＴ程序输入、输出文件说明