弯扭耦合对风力机叶片气弹模态特性的影响.pdf
《弯扭耦合对风力机叶片气弹模态特性的影响.pdf》由会员分享,可在线阅读,更多相关《弯扭耦合对风力机叶片气弹模态特性的影响.pdf(7页珍藏版)》请在咨信网上搜索。
1、第 11 卷 第 4 期 新 能 源 进 展 Vol.11 No.4 2023 年 8 月 ADVANCES IN NEW AND RENEWABLE ENERGY Aug.2023 *收稿日期:2023-03-08 修订日期:2023-04-12 基金项目:国家自然科学基金项目(51776044,51105079);广东省自然科学基金项目(2020A1515010844)通信作者:夏鸿建,E-mail: 引用本文:曲珍壮,夏鸿建,李德源,等.弯扭耦合对风力机叶片气弹模态特性的影响J.新能源进展,2023,11(4):333-339.Citation:QU Zhenzhuang,XIA Hon
2、gjian,LI Deyuan,et al.Effect of bending-torsion coupling on aeroelastic modal characteristics of wind turbine bladesJ.Advances in new and renewable energy,2023,11(4):333-339.文章编号:2095-560X(2023)04-0333-07 弯扭耦合对风力机叶片气弹模态特性的影响*曲珍壮,夏鸿建,李德源,林伟豪(广东工业大学 机电工程学院,广州 510006)摘 要:弯扭耦合对风力机叶片气动弹性响应有重要影响。为研究弯扭耦合对风
3、力机叶片气弹模态特性的影响,采用几何精确梁理论和叶素动量定理,建立风力机叶片的气弹耦合方程,在 Beddoes-Leishman 动态失速模型的基础上对气弹耦合方程线性化,构建特征平衡方程。由复特征值分析法计算 DTU 10 MW 叶片的气弹频率和气弹阻尼比。通过在截面刚度矩阵中引入耦合系数,考虑弯扭耦合对气弹模态的影响。挥舞顺桨耦合会增加挥舞模态气弹频率,降低挥舞模态气弹阻尼比。挥舞失速耦合会降低挥舞模态气弹频率,增大挥舞模态气弹阻尼比。挥舞顺桨耦合会降低风力机叶片的颤振风速。关键词:风力机叶片;几何精确梁;弯扭耦合;颤振分析 中图分类号:TK83;O313.7 文献标志码:A DOI:10
4、.3969/j.issn.2095-560X.2023.04.005 Effect of Bending-Torsion Coupling on Aeroelastic Modal Characteristics of Wind Turbine Blades QU Zhenzhuang,XIA Hongjian,LI Deyuan,LIN Weihao(School of Electromechanical Engineering,Guangdong University of Technology,Guangzhou 510006,China)Abstract:The bending-tor
5、sion coupling has an important influence on the aeroelastic response of wind turbine blades.To study the effect of bending-torsion coupling on the aeroelastic modal characteristics of wind turbine blades,the geometrically exact beam theory and blade element momentum theory are used to establish the
6、aeroelastic coupling equation for wind turbine blades.Based on the Beddoes-Leishman dynamic stall model,the aeroelastic coupling equation is linearized and the characteristic equilibrium equation is constructed.The DTU 10 MW blades aeroelastic frequencies and damping ratio are calculated using compl
7、ex eigenvalue analysis.The influence of bending-torsion coupling on the aeroelastic modes is considered by introducing coupling coefficients into the section stiffness matrix.The flap-twist to feather coupling increases the flapwise aeroelastic frequency and decreases the flapwise aeroelastic dampin
8、g ratio,while the flap-twist to stall coupling decreases the flapwise aeroelastic frequency and increases the flapwise aeroelastic damping ratio.The flap-twist to feather can reduce the flutter speed of wind turbine blades.Keywords:wind turbine blade;geometrically exact beam;bending-torsion coupling
9、;flutter analysis 0 引 言 弯扭耦合对风力机叶片气动弹性响应有重要影响。叶片在气动力作用下弯曲后,会改变攻角的大小,进而改变叶片的气动力,形成反馈循环。弯扭耦合可以由叶片的几何形状(几何耦合)或各向异性叶片材料(材料耦合)引起。几何耦合是叶片几何形状弯曲的结果(预弯曲、后掠或负载偏移),从而引起额外的扭转。材料耦合是叶片表皮或梁帽的纤维方向偏离主方向,从而引起弯曲和扭转的耦合。弯扭耦合可以用来减轻风力机叶片的负载。根据 LOBITZ1、BOTTASSO2等的研究,挥舞扭转顺桨耦合会使风力机叶片的疲劳载荷降低10%20%。然而弯扭耦合除了降低叶片负载外,还会改变叶片334 新
10、 能 源 进 展 第 11 卷 的气弹模态特性。LOBITZ3研究了弯扭耦合叶片与非弯扭耦合叶片发生颤振失稳临界风速,发现弯扭耦合叶片发生颤振失稳的临界风速会比非弯扭耦合叶片降低 12%左右。LARWOOD4研究后掠叶片的颤振特性,发现后掠叶片颤振失稳的模态是二阶挥舞与一阶扭转的耦合模态。随着风力机叶片长度和柔性的增加,叶片的非线性效应越加显著。KALLESE5采用非线性二阶 Bernoulli-Euler 梁理论,研究了 5 MW 叶片弯曲变形对气弹稳定性的影响,表明叶片的几何非线性对其气弹稳定性有明显的影响。传统采用欧拉梁结构模型6或模态法7构建叶片结构模型,难以表达柔性叶片的几何非线性变
11、形。由于风力机叶片气弹耦合状态方程的建立比较复杂,基于频域的颤振分析研究相对较少,黄俊东等8采用多体动力学方法和叶素动量理论的方法建立叶片气弹耦合方程,研究了美国可再生能源实验室(National Renewable Energy Laboratory,NREL)发布的 5 MW 风力机气弹稳定性问题。和其他非线性模型相比,几何精确梁模型可以更为准确方便地计算叶片的变形位置以及叶片扭转的角度。通过在截面刚度矩阵中加入耦合项,研究由叶片材料引起的弯扭耦合问题。在结构动力学求解时,对结构动力学方程进行线性化处理,可以方便地与气动力线性化结合起来,进行特征值分析。本文首先根据几何精确梁理论,结合叶素
12、动量定理,建立风力机叶片气弹耦合方程,找出叶片在不同工况下的平衡位置。然后对气动力进行线性化,并引入 Beddoes-Leishman(B-L)动态失速模型,构建特征平衡方程,采用特征值分析方法,研究风力机叶片的气弹模态特性。通过在截面刚度矩阵中引入耦合系数9,以丹麦技术大学(Technical University of Denmark,DTU)10 MW RWT10叶片为参考叶片,研究了挥舞扭转耦合对气弹模态特性的影响。1 几何精确梁变形的描述 几何精确梁采用平截面假设,梁的位形由梁轴线的位移和相应截面的转动决定。如图 1 有全局坐标系(不考虑叶片预弯预扭以及截面偏心等因素)、不变形位置坐
13、标系(考虑叶片真实几何形状)和变形位置坐标系(考虑气弹耦合后叶片的变形)。图 1 梁的转动坐标系 Fig.1 The rotating coordinate system of the beam 在变形位置,梁上任意一点 p 的位置向量可以写为,ppt z x yt zx yuuR (1)式中:pu为剪切中心在全局坐标系下的坐标;z 为截面在梁轴线上的坐标;为截面转动矩阵;pR为截面坐标系下的位置向量。叶片坐标系的建立为,在叶片根部以各个截面的剪心为原点,建立局部坐标系,将局部坐标系平移和旋转后,得到叶片的真实几何形状。三维旋转矩阵 有 9 个元素,包含 6 个约束关系,难以直接作为姿态变量,
14、欧拉四元数是描述截面转动非奇异的最小参数表达,本文采用欧拉四元数 q描述截面的转动11。T00123,qq q q qqq (2)式中:四元数之间满足222201231qqqq;T()表示向量的转置。四元数 q和转动矩阵 的关系为 2T0021q2q2I+q+qq (3)2 结构动力学方程的建立 2.1 几何精确梁本构关系 几何精确梁理论基于 Reinssner 应变,根据参考文献12-13,定义应变 为平动应变向量,应变 K为转动应变向量,截面等效应力 N、M 和应变、K之间的关系为 NSMK (4)第 4 期 曲珍壮等:弯扭耦合对风力机叶片气弹模态特性的影响 335 式中:S 为 6 6
15、截面刚度矩阵。不考虑耦合时,S矩阵只包含对角线元素。通过在截面刚度矩阵 S 中设置耦合项46K可以引入挥舞扭转耦合。464466KK K (5)式中:为挥舞扭转耦合系数9;44K和66K分别为挥舞刚度和扭转刚度。应变与变形的关系为 T0u (6)T*12300,2K KKqqK (7)式中:0和0K为梁的初始位置和曲率参数;()为在参考构形中对截面在 z 轴上坐标的导数;为四元数乘法。2.2 动力学平衡方程 首先建立风力机叶片弱形式动力学平衡方程14,然后通过数值积分和微分方法离散求解。对内力沿单元轴线积分,得到内力虚功为15 eeTTint0dlZWNK M (8)对气动力沿轴线积分,重力进
16、行体积积分,得到单元的外力虚功 eeTTTextaaaap0ddlVZVWu fmug (9)式中:au为轴线气动节点的虚位移;a 为轴线节点气动点的虚转动;af为轴线节点的气动力向量;am为轴线节点的气动弯矩向量;g 为重力加速度向量;为叶片体密度;el为单元长度;下标 e 表示叶片单元。各点的速度进行体积分,得到单元的惯性力虚功 eTinertdppVVWuu (10)根据虚功原理,单元内力、外力和惯性力虚功满足 eeeintextinert0WWW (11)通过 Gauss-Lobatto 积分和微分求积法14离散平衡方程,得到基于节点虚位移表达的内力、外力和惯性力虚功15。eintin
17、teextexteinertinertWdRWdRWdReTeTeT (12)单元节点 k 虚位移表示为 kkkkxkykzkxyzd=(13)单元位移e12 Kdddd。把式(12)代入式(11),并在全局坐标系下进行组装,得到等效节点内力、外力和惯性力的离散平衡方程 intextinert0RRR (14)3 动力学方程线性化 为了构建叶片气弹线性平衡方程,形成特征值求解,需要对动力学平衡方程进行线性化处理。分为两个部分,第一部分是结构动力学方程线性化,第二部分是气动力线性化。3.1 结构动力学线性化 分别对节点内力、外力和惯性力矩阵线性化15得到 intintRKd (15)Taexte
18、xtRKdBFB d (16)inertinertinertinertRKdCdMd (17)将式(15)式(17)代入平衡方程(14)可以得到 Ta0M dC dK dBFB d (18)式中:intextinertKKKK;inertMM;inertCC;为弱形式求积长度系数16;为Gauss-Lobatto积分系数;B为微分求积定位矩阵。3.2 气动力线性化 根据叶素动量定理,作用在叶片第k个截面的气动力可以写为 a2a12kkkkc WfC (19)式中:c为截面弦长;a为空气密度;W为相对来流风速;TLDM,CCCC,其中LC、DC、MC分别为升力系数、阻力系数和扭矩系数。336 新
19、 能 源 进 展 第 11 卷 考虑叶片振动和入流角的影响,相对来流风速W和攻角可以写为 WVnxRy (20)1VnxtgRy (21)式中:Vn为垂直旋转平面的来流风速;为入流角;为桨距角;为叶片扭角;y 和x 分别为摆振和挥舞方向振动速度;R为翼型截面剪心到轮毂中心的距离;为叶片旋转速度,见图2。图 2 叶片翼型截面 Fig.2 Blade airfoil section 当来流风速和叶片旋转速度均已确定时,叶片在某一瞬时受到的气动力是振动速度的函数。为研究在某一工况下气动力对叶片振动特性的影响,首先将af进行坐标转换,得到全局坐标下的气动力aF,然后采用变分原理对气动力进行线性化处理。
- 配套讲稿:
如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。