大型风力机叶片气动结构一体化没计方法研究.pdfVIP

会议及刊物论文征集I 型风力机叶片最为显著的气动与结构耦合现象。探讨将多学科优化方法应用于风力机叶片的 气动结构一体化设计是必要的。本文采用多学科优化理论，采用叶素动量理论和有限元结构 分析方法，对大型风力机叶片进行了气动／结构一体化设计。 1．气动分析 目前，风力机叶片的气动分析普遍采用叶素动量理论(BEM)，由于计算效率高和讨算精 度也相对较高，在风力机叶片的气动分析和气动外形设计中得到广泛应用。但是根据叶素动 量理论(BEM)建立的叶片气动模型都是刚性的，不考虑叶片运彳了结构的弯曲和扭转对气动外形 的影响，从而也影响了叶片气动力计算的准确度。同时，也要求将各叶素的气动载荷施加给 结构模型，分析气动载荷变化对结构模型的影响。NREL开发的％VT_Perf适用于建立考虑构 变形影响的气动模型，并可给出各叶素的气动载荷和输出功率，具有与结构模型集成的优良 特性。 叶片在气动力作用下发生弯曲变形和扭转变形，弯曲变形改变了叶素剑旋转中心的距离， 降低了叶素切向的速度，扭转变形改变了叶素的攻角，二者都对叶素的气动载荷产生重要影 响。图1为在半径r处的叶素上的气流变化情况。在半径r处叶素的沿展向的方向．；，法向商， ．， 几何扭角为￡，，转角为p。则来流风速K在叶素处的分量圪。为 ‰=Kcosp0一口) 图l叶片弯曲变形对叶素气流的影响 叶素在旋转平面内的速度分量V,o为： 匕。=D．r(1+b) 合成风速％=压巧，雷诺数Re=警空气密度p=1．225姆脯’c是叶素的弦长， ∥是空气粘性系数，a，b分别为轴向和切向诱导因子。 叶素随结构在扭矩的作用下发生顺桨方向的扭转，扭转角度缈叶片r半径处叶素的入流角① 和攻角口可表示为： ①一n掣 口=①一缈一0 431 ■■一W-I舀arV_．1会议及刊物论文征集 合成气流速度K在叶素上的升力L和阻力D分解为法向力只和切向力f， E=￡cos≯+Dsin砂 只=Lsin矽一Dcos≯ 载荷计算过程中采用了Glauert修正、叶尖和叶根损失模型。为了简化气动载荷传递给结 构模型的过程，而又不失一般性，结构节点与叶素的中心位置在展向重合，并且假设叶素上 的气动力可以用叶素中心的集中力替代。 图2叶素上的气流速度分量 2．结构分析 叶片结构模型—般采用梁模型。本研究采用了六自由度Timoshenko梁单元建立悬臂梁模 型，梁节点位于叶片的变矩中心。风力机叶片运行中主要受气动载荷、离心力载荷和重力载 荷，气动力载苟由气动模型提供，为了简化气动载荷在结构节点上施加过程，叶素中心与对 应的结构模型节点重合。离心力和重力通过叶片旋转角速度和重力加速度来实现。同时结构 模型要满足强度和变形要求。 叶片大部分截面不是标准翼型，叶片结构模型的单元属性通过叶片截面的几何性质赋予， 需要的参数包括沿x轴和Y轴的惯性矩，。、，，和惯性积，。，绕叶片变矩中心轴的惯性矩等， 所用材料的弹性模量E和剪变模量G。截面的惯性距和惯性积为： I。=I旷dA j Iy=l妒dA