水平轴风力机叶片气弹建模与响应分析 Modeling and analysis of the aeroelastic response of the horizontal axis wind turbine blade.pdfVIP

第30卷第2期 空气动力学学报 V01．30．No．2 2012年4月 ACTAAERODYNAMICASINICA Apr．，2012 文章编号：0258—1825(2012)02—0192一06 水平轴风力机叶片气弹建模与响应分析 陆 洋， 周桂林 (南京航空航天大学旋翼动力学国防科技重点实验室，江苏南京210016) 摘 要：基于Hamilton原理在惯性系下建立叶片动力学方程。针对风力机预弯叶片，引入上反坐标系描述叶片的 流模型采用Suzuki广义动态尾迹人流模型。利用Newmark数值积分方法获得叶片气弹响应的稳态周期解。分别 以美国国家可再生能源实验室PhaseVI非定常空气动力学实验及其公开的1．5Mw风力机叶片为算例计算了有／ 无预弯叶片的气弹响应，验证了本文所建模型的正确性。 关键词：风力机；气弹响应；惯性系；预弯叶片；建模 中图分类号：TK83 文献标识码：A 片的动态响应。 IJ 引 吾 可见，为了保证叶片响应的计算精度，有限元方 大型水平轴风力机系统是强非线性流一刚柔耦合 法已逐渐替代模态法成为风力机叶片结构动力分析 的周期时变多体系统，结构和运动关系非常复杂，对 的主流方法，气动力建模则仍以动量叶素理论和二维 风力机叶片进行动力学建模时必须考虑叶片的几何 准定常气动力模型为主。 非线性以及非定常气动载荷等因素，理论推导和数值 然而，现有的叶片气弹建模方法仍存在以下局限 计算都比较困难。另一方面，随着风力机额定功率不 性：1)基于有限元方法的风力机动力学方程通常建 断增加，叶片尺寸越来越大，细长叶片在交变载荷下 立在非惯性系下，即处理动能部分时将其分为旋转坐 的挥舞变形更大。为防止叶尖与塔架碰撞，出现了预 标系的动能和牵连运动引起的动能两部分，破坏了推 弯叶片设计概念，增加了风力机叶片气弹建模的难 导和编程的一致性；2)建模时，没有考虑叶片预弯对 度。 叶片气动弹性动力响应的影响；3)气动力模型多采用 风力机叶片气弹建模包括叶片结构动力学建模 二维准定常空气动力计算模型，没有考虑叶片动态人 和空气动力学建模。文献[1]采用模态分析法进行叶 流以及动态失速影响。这些局限性，使得传统建模方 法的计算精度和适用性受到了制约。 片结构动力学建模，采用Greenberg’s二维准定常气 动力模型计算气动力，研究了风力机气弹稳定性；文 针对前述叶片建模方法的局限性，本文开展了如 献E23采用5节点18自由度梁单元对叶片进行有限 元离散，采用二维准定常气动模型，利用拟线性法计 阵[5]，通过递归计算的方式得到惯性坐标系下的位 算了水平轴风机工作状态下的气动弹性响应；文献 置矢量，从而在惯性系下建立了风力机叶片动力学方 程，将动能部分合为一体，避免了传统建模中的繁琐； E3]以Marc有限元软件为基础平台，采用动量叶素 理论二维气动力模型，计算叶片的气弹响应；文献[4] 2)引入上反坐标系方法描述风力机叶片的预弯变形， 采用二节点梁单元对叶片进行有限元离散，采用 的针对风机叶片的广义动态尾迹人流模型(GDW)[6]， Leishman—Beddoes非定常模型计算气动力，分析了叶 09 *收稿日期：201109；修订日期：2011-11-22 万方数据 第2期