风力机旋转叶片的多体动力学数值分析 multi-body dynamics numerical analysis of rotating blade of horizontal axis wind turbine.pdfVIP

第26卷第4期 太阳能学报 V01．26．No．4 2005年8月 ACrAENERGIAESOLARISsINICA Aug．，2005 风力机旋转叶片的多体动力学数值分析 李德源，叶枝全，陈 严 (汕头大学能源所，汕头515063) 摘要：在现代柔性多体系统动力学理论的框架上，研究了风力机叶片的空间旋转运动与其弹性变形间的耦合 关系；对旋转叶片进行有限元离散，采用凯恩(Kane)方法，建立了一般形式的大型风力机叶片柔性多体动力 学方程；结合大型叶片几何构形特点，导出了旋转叶片的空间梁单元动力矩阵的显式形式；以1．5MW风力机叶 片为例，计算了该叶片的固有频率并与常规有限元数值分析结果进行了比较，体现了旋转叶片的动力刚化效应， 更准确地反映了叶片的动力学特性。 关键词：旋转叶片；多体动力学；l(aIle方程；固有频率；动力刚化 中图分类号：Ⅱ(83 文献标识码：A 应等[6娟]。本文在现代计算柔性多体系统动力学理 O 引 言 论的框架上，考虑风力机叶片的空间旋转运动与其 随着风力机功率的越来越大，为了从风中吸取更 弹性变形问的耦合关系，应用凯恩(Kane)方法，对旋 多的能量，叶片必须设计成足够长的尺寸，这使其柔 转叶片进行有限元离散，用有限单元节点位移表示 性进一步增加，对其运行精度的控制要求越来越高， 叶片上任一点的弹性变形，建立了大型风力机叶片 需要准确地了解叶片工作过程中的动力学响应。风 柔性多体动力学方程。为了方便工程实际应用，结 力机工作时，叶片绕中心轴转动，叶片的变形使其绕 合大型叶片展弦比较大，可以简化为空间梁元的特 转轴的转动陨量及质量分布等发生改变，从而影响叶 点，导出了旋转叶片空间梁单元矩阵的显式形式，以 片的转动形态，而叶片转动形态的改变，必然使叶片 1．5Mw风力机叶片为例，计算了其固有频率，从中 的弹性变形发生改变，叶片的转动与柔性叶片弹性变 可以看出多体动力学理论较常规动力学理论更能准 形是相互耦合的，可见，旋转叶片的动力学分析不同 确反映旋转叶片的动力特性。 于传统的不考虑结构的刚体位移的结构动力学分析， 它需要考虑叶片的转动与其弹性变形问的相互耦合 1旋转叶片的有限元动力学方程 关系及其所导致的动力学效应，该问题属于现代柔性 为了分析叶片的运动，引入如图l所示的坐标 多体动力学的分析范畴。近年来，该理论在风力机领 系，O郴为固定坐标系，o’戈7，，727为固定在旋转叶片 域日益受到重视，已被越来越多的研究者引入风力机 上一点的动坐标系，图中画出了绕定轴z旋转的风 结构动力学分析中¨。j。 力机叶片的一部分。设[A]为动系相对于固定坐标 研究物体变形与其刚性运动的相互作用或耦 系的方向余弦矩阵，它反映弹性叶片的刚体转动， 合，以及这种耦合所导致的动力学效应，正是柔性多 [A]的表达式为： 体系统动力学的核心特征。旋转叶片可视为绕转轴 rcos(茗，戈7)cos(’，，戈7)cos(z，戈7)] 旋转的旋臂梁，国际上对旋转梁的建模大都是对连 [A]=I 续系统采用能量原理等建立偏微分形式的系统方 Lcos(茄，彳7) 程，并用各种试函数求旋转梁的固有频率和位移响 (1) 收稿日期：200B．1l掰 基金项目：国家高技术研究发展计划(863计划)专项经费； 国家自然科学基金(50∞6023)；汕头大学博士后 基金资助 万方数据万方数据 474 太 阳 能 学 报