水平轴风力机叶片自振频率计算方法幻灯片.pdfVIP

2008年第3期 玻璃钢／复合材料 35 水平轴风力机叶片自振频率计算方法研究 陈雨，袁国青 (同济大学航空航天与力学学院，上海200092) 摘要：本文应用商业有限元软件，通过选取不同的单元类型(梁单元和壳体单元)对国内200kw风力机叶片进行模态分 析。结果表明，两种方法精确相符，从而证明在计算叶片自振基频时可把复杂叶片结构简化为一悬臂梁来进行分析。鉴于把 叶片简化为一悬臂梁进行结构动力有限元计算易于实现计算的程序化，基频计算具有足够的精度且计算效率高，因此其对于 叶片的初步结构设计具有一定意义，在分析叶片基频时不必建立复杂的三维叶片有限元模型。 关键词：风力机叶片；振型；自振频率；有限元法 中图分类号：TK83文献标识码：A 文章编号：1003一0999(2008)03—0035一04 1 前言 方法有瑞利法、迭代法、矩阵分析法和有限元法。对 振动是工程结构中经常遇到的现象，风力机也 于迭代法和矩阵分析法【2】，主要求解弯曲振动频率 不例外。当风力机在自然条件下运营时。作用在风 及扭转振动频率，不能够求解耦合振动频率。瑞利 力机叶片上的空气动力、惯性力和弹性力等交变荷 法求解一阶挥舞和摆振方向自振频率具有很高的精 载会使叶片结构产生振动，有可能影响风力机的正 度，但计算叶片弯扭耦合振动频率方法不可靠，误差 常运行，为防止灾难性振动的发生，在设计风力机时 很大【3】。有限元法求解叶片结构的自振特性是一 就必须进行动力学分析，包括对叶片结构自振特性 种可靠而有效的方法，但是风力机叶片的结构设计 的分析。 过程往往是一个逐步调整结构布局和材料铺层方案 风力机叶轮叶片的振动主要有三种振动形式： 反复计算的过程，再加上叶片结构的复杂性，三维有 ①挥舞，是叶片在垂直于旋转平面方向上的弯曲振 限元模型建模相当困难，数值分析的效率较低，因此 动；②摆振，是叶片在旋转平面内的弯曲振动；、③扭 在叶片的初步设计中很少直接应用三维壳单元或实 振，是绕叶片变距轴的扭转振动。此外，在空气动 体单元进行有限元分析，只是到详细设计阶段才会 力、惯性力和弹性力的耦合作用下，这三种形式的振 选用该法进行分析。 动还会发生耦合，引起气动弹性问题。一般风力机 由于叶片的剖面尺寸远较叶片的长度为小，故 叶片气动弹性动态不稳定包括挥舞一扭转不稳定、扭 可以把叶片结构的计算模型简化为一根端固定的悬 转一摆振不稳定、失速颤振和失速诱导的摆振等…。 臂梁。本文拟把复杂的叶片结构简化为一端固定一 端自由的悬臂梁进行结构动力特性分析，从而使创 实际上作用在叶片上的荷载为动荷载，且具有 交变性和随机性，当叶片的固有频率与激振力的频 建叶片结构的有限元模型大大简化。同时，应用商 率相同时，就会产生共振，从而影响叶片的安全。’由 业有限元软件对风力机叶片进行三维有限元分析， 于动荷载很难用数学式来描述，所以工程上一般并 以证明此方法的可靠性。 不对叶片在动荷载作用下的动态响应进行动态分 2计算结构固有频率和振型的有限 析，而是着重于分析计算叶片的自振频率，以防发生 元法 共振。一般叶片自振频率不与叶轮的转速频率的整 由结构动力有限元可知，叶片总体结构有限元 数倍重合，就可以避免发生共振，实践中最重要的是 模型的运动方程为： 避免低阶自振频率与叶轮的转速频率的整数倍重