复杂工况下的大型风力机气动性能和尾迹研究Investigation of aerodynamic performance and wake of the large scale wind turbine in complicated conditions

在复杂工况下,大型风力机非定常特性会更严重,导致风力机气动性能变化和尾迹预测更加复杂。本文主要针对稳态偏航、动态偏航、风剪切和随机风速场等复杂工况,基于自由涡尾迹方法,嵌入复杂工况的模块,加入了动态失速模型和三维旋转效应模型修正,实现了复杂工况数值模拟计算,比较了不同复杂工况的气动载荷和尾迹形状。最后,得出了风力机在复杂工况下的气动性能、载荷和尾迹叶尖涡线特性,并计算出风力机在复杂工况下的气动载荷超调量和迟滞时间。对推进自由涡尾迹方法应用于风力机工程的大批工况载荷计算,提高大型风力机的载荷计算精度和设计水平等具有重要意义。     

复杂工况下的大型风力机气动性能和尾迹研究

在复杂工况下,大型风力机非定常特性会更严重,导致风力机气动性能变化和尾迹预测更加复杂。本文主要针对稳态偏航、动态偏航、风剪切和随机风速场等复杂工况,基于自由涡尾迹方法,嵌入复杂工况的模块,加入了动态失速模型和三维旋转效应模型修正,实现了复杂工况数值模拟计算,比较了不同复杂工况的气动载荷和尾迹形状。最后,得出了风力机在复杂工况下的气动性能、载荷和尾迹叶尖涡线特性,并计算出风力机在复杂工况下的气动载荷超调量和迟滞时间。对推进自由涡尾迹方法应用于风力机工程的大批工况载荷计算,提高大型风力机的载荷计算精度和设计水平等具有重要意义。     

台风过境下大型单桩式海上风机结构动力特性研究

风机大型化是我国海上风电技术发展的重要方向. 东南沿海是我国海上风电发展的重要基地, 这一区域频繁发生的台风对海上风机的影响不可忽略. 台风风场与常规大风风场有不同的湍流特性, 同时台风期间较高的风速会引起巨大的台风浪. 本文考虑台风经过期间独特的风场及波浪场, 开展风浪联合作用对大型单桩海上风机影响的研究. 基于DTU 10 MW大型单桩风机, 运用一体化分析软件SIMA建立风浪联合作用下大型单桩风机的耦合数值模型, 研究台风经过不同阶段大型风力机的动力响应特性. 计算结果显示, 叶片变桨能有效降低台风经过时风机叶片所受风载荷, 变桨状态下单桩风机所受风载荷主要来源于塔筒. 在台风经过的不同阶段, 大型单桩海上风机结构表现出不同的动力特性. 台风全过程塔筒运动均受波浪激发一阶频率控制, 塔基上方结构动力载荷以惯性载荷为主, FOVS至FEWS阶段及BOVS阶段至BEWS阶段塔筒运动一阶频率处响应能量增长较小, 响应能量向低频及波频转移. 塔基下方泥面线处剪力响应受波频控制, 弯矩响应受一阶频率控制.      

一种风力机气动计算的全自由涡尾迹模型

采用全自由方式建立风力机尾流场的涡尾迹模型,引入“虚拟周期”的概念,并发展一种自适应松弛因子方法,从而改善了自由尾迹迭代的稳定性,提高了迭代收敛速度。利用建立的自由涡尾迹模型,计算了风力机叶片的尾流场结构、气动性能及叶片载荷,并与实验结果进行了对比分析。结果表明,尖速比越大,自适应松弛因子方法对缩小模型计算时间越有效;全自由涡尾迹模型能准确给出风力机尾流场的结构,包括尾迹的扩张以及叶尖涡和叶根涡的产生、发展和耗散的过程,风轮扭矩与实验数据吻合;叶片载荷分布的计算结果在低风速下与实验值基本一致,但是在大风速下差别较大,说明需要一个准确的失速模型。     

一种磁力滑动式翼型颤振能量俘获器

风致振动是自然界中普遍存在的一种现象,并且蕴藏着巨大的可利用能源.如何充分利用风致振动引起的结构大幅值响应进行能量俘获,为微电子器件供电是能量俘获领域的一个难题.为了高效俘获风致振动能量,文章提出了一种磁力滑动式翼型颤振能量俘获器.基于半经验非线性空气动力学模型并考虑与磁铁位置相关的机电耦合系数,建立了该能量俘获器的动力学模型,搭建了风洞实验平台,制作了实验样机.通过增加风速和降低风速的方式为能量俘获器提供两种不同的初始状态,发现其具有两个临界风速(5.2 m/s和8.3 m/s),降风速实验中在8.3 m/s风速下出现突跳现象.在数值仿真中,在6.8 m/s和8.2 m/s风速下出现了两个突跳点,和一段多解区域.分析了沉浮位移和电压响应,发现沉浮位移以正弦形式响应,输出电压以非正弦形式响应,并出现明显的偶次谐波.仿真的沉浮位移和电压输出波形与实验波形吻合较好,验证了模型的准确性.能量俘获器的均方根电压随电阻的增加而增加,平均功率随电阻增加呈现先增加后降低的趋势.分析了负载电阻对能量俘获性能的影响,在8.6 m/s风速下,实验中能量俘获器的负载电阻接近线圈内阻值时平均功率达到最大值7....     

模态修型减振的旋翼桨毂谐波振动载荷计算

基于模态修型的方法推导了直升机旋翼桨叶根部剪力的计算公式,通过合成桨叶根部载荷推导了旋翼桨毂的谐波振动载荷公式。将本文建立的桨毂谐波载荷计算模型与商用软件CAMRAD计算出的结果进行对比,结果表明:本文所建立模型的旋翼计算频率与CAMRAD计算的频率相比,基阶频率计算误差在2%以内,前十阶频率误差都在8%以内;与已有文献试验测试的桨毂谐波振动载荷相比最大误差在25%以内。说明了本文建立的旋翼桨毂谐波振动载荷计算模型具有可行性和有效性。     

两种不同弦长螺旋桨的风洞试验及分析

本文介绍了两种不同弦长螺旋桨的风洞试验结果。试验在西北工业大学NF-3风洞的三元试验段进行,试验风速分别为V=20m/s、30m/s和40m/s,每个风速下,螺旋桨旋转速度为:900r/min、1200r/min、1500r/min、1800r/min、2100r/min、2400r/min、2700r/min、3000r/min、3300r/min、3600和3900r/min。风洞试验结果表明:当两种螺旋桨的翼型相同、桨叶角沿径向分布相同,但弦长沿径向的分布不同时,它们的拉力、扭矩、功率、效率以及前进比会有显著区别。其中,1#螺旋的弦长大于2#螺旋桨,在相同的试验风速和螺旋桨转速下,1#螺旋桨的拉力、扭矩和功率高于2#螺旋桨,但效率低于2#螺旋桨。     

端面宽度对机械密封极限pcv (pressure×velocity)值的影响

端面宽度是机械密封的重要设计参数,为了研究端面宽度对机械密封极限p_cv值(端面比压p_c×端面平均线速度v)的影响,本文作者选取浸呋喃树脂石墨与微孔常压烧结碳化硅作为密封摩擦副,采用定速度变载荷和变速度定载荷两种不同测试方法获取了五种端面宽度摩擦副的极限p_cv值,分析了不同端面宽度摩擦副在不同测试方法下的摩擦系数时变特性和密封介质温度时变特性,并通过扫描电子显微镜和三维激光形貌仪分析了摩擦副的表面磨损特征. 结果表明:在相同测试方法下,与较宽端面摩擦副相比,较窄端面摩擦副的极限p_cv值较高并且密封失效时摩擦系数较小而密封介质温度较高;在给定端面宽度条件下采用定速度变载荷方法测试时,较高速度15.68 m/s对应摩擦副的极限p_cv值相对于较低速度7.84 m/s时较高,其中极限p_c值较低;达到极限p_cv值时,摩擦副表面的磨损机理主要为黏着磨损,其诱因在于端面摩擦扭矩的突增引起端面过热,导致两密封端面因液膜汽化加剧发生干摩擦. 定速度变载荷和变速度定载荷两种方法均可用于测试机械密封极限p_cv值,但采用变速度定载荷方法测试时密封失效对应的摩擦系数突增幅度相对较不明显,失效的现象较难捕捉,石墨表面被破坏的程度也较轻. 在满足材料强度的要求下,减小端面宽度,有效避免端面过热,有利于机械密封达到更大的极限p_cv值. 本文为机械密封极限p_cv值的测量以及端面宽度的设计和密封的延寿方法提供了指导.      

攒尖四坡屋面风致雪漂移的数值模拟

为获得攒尖四坡屋面的风致雪漂移规律，基于欧拉-欧拉方法和风雪单向耦合假定，运用计算流体动力学软件，选取Mixture模型分别对立方体周边和高低屋面上的风致雪漂移运动进行数值模拟，将模拟结果与两者的实地观测数据对比，探讨分析数值风洞的关键技术和参数设置，验证数值模拟方法的合理性与可靠性。依据攒尖四坡房屋的使用功能要求，设计分析模型与分析工况，在试算的基础上对屋面进行分区。以风速5 m/s，7 m/s，9 m/s，11 m/s，13 m/s和15 m/s，风向角0°，15°，30°和45°以及屋面坡度25°，30°，34°，40°和45°为分析参数，对攒尖四坡房屋的120种工况进行数值模拟，得到屋面各分区侵蚀沉积的基本规律，提出可用于抗雪设计的屋面积雪分布系数。研究表明，风向角的改变会使屋面积雪分布状态发生较大程度的变化，风速和屋面坡度的变化对屋面整体积雪量有较大影响。     

基于动态磁耦合的驰振能量收集器动力学分析

为了提高能量收集系统在低风速下的能量收集效率,将动态磁铁非线性引入到驰振能量收集系统中。在悬臂梁的末端和底座上分别安装一对磁极相斥的磁铁,其中安装在底座上的磁铁与弹簧相连,可随着磁斥力的变化而垂直移动。首先,根据能量法建立了磁耦合驰振能量收集系统的多场耦合振动控制方程。其次,通过Runge-Kutta数值计算方法比较分析了低风速下动态磁耦合驰振能量收集系统(DM-GEH)和固定磁耦合驰振能量收集系统(FM-GEH)的电压输出。DM-GEH系统的切入风速提前了81.82%,在1 m/s～5 m/s风速范围内能量收集效率提高了124.22%。最后,针对弹簧支撑刚度进行参数优化,提升了低风速下的能量收集效率。结果表明,通过改变磁铁支撑方式至弹性支撑将改变系统的振动频率并且降低切入风速,相较于弹簧刚度为1 000 N/m时,弹簧刚度为500 N/m时的系统的切入风速降低了54.55%,能量收集效率提高了15.35%。     

Aeroelastic analysis of wind turbines using a tightly coupled CFD–CSD method

This paper presents aeroelastic analyses of wind turbines, using the compressible flow Helicopter Multi-Block (HMB2) solver of Liverpool University, coupled with a Computational Structural Dynamics method. For this study, the MEXICO and NREL Phase VI wind turbines were employed. A static aeroelastic method was first employed for the analysis of the MEXICO blade and the effect of the torsional stiffness was studied at 10, 15 and 24 m/s axial wind speeds. The torsional deformations showed strong dependency on this parameter and the blade region from mid-span to the tip was the most susceptible to aeroelastic effects. The work progressed by studying both the static and dynamic response on the NREL wind turbine, where the nacelle and the tower were considered. Mean deflections between the static and dynamic methods showed consistency and, due to the structural properties of this blade, flapping modes were dominant. The dynamic aeroelastic method enabled an assessment of the effect of flapping on the blade loads, in conjunction with the effect of tower. Aeroelastic effects were found to be secondary for the MEXICO blade, but had a stronger effect on the larger NREL Phase VI blade.     

Experimental study of the effect of slotted blades on the Savonius wind turbine performance

This study investigates the effect of Reynolds number on the performance of Savonius wind turbine with slotted blades. The turbine performance investigation was based on the torque coefficient( Ct), power coefficient( Cp), and tip speed ratio( TSR). The experiment used two number of blade configuration, blade overlap ratio of 10%, 12.5% and 20%, slotted position of 15%, 20%, 25% and 35%, and also slotted gap width of 3 mm, 5 mm, 7 mm, and 9 mm. The wind speed carried out in this experiment are 5.94 m/s, 6.46 m/s, 6.99 m/s, and 7.27 m/s, which are generated from the fan blowers as a wind source. The Savonius turbine with 10% overlap ratio shows the best performance. The highest Cp obtained is 0.138 by the variation of a 3 mm gap with Re of 1.44 × 10~4 and 0.526 TSR.     

Aeroelastic analysis of large horizontal wind turbine baldes

A nonlinear aeroelastic analysis method for large horizontal wind turbines is described. A vortex wake method and a nonlinear finite element method (FEM) are coupled in the approach. The vortex wake method is used to predict wind turbine aerodynamic loads of a wind turbine, and a three-dimensional (3D) shell model is built for the rotor. Average aerodynamic forces along the azimuth are applied to the structural model, and the nonlinear static aeroelastic behaviors are computed. The wind rotor modes are obtained at the static aeroelastic status by linearizing the coupled equations. The static aeroelastic performance and dynamic aeroelastic responses are calculated for the NH1500 wind turbine. The results show that structural geometrical nonlinearities significantly reduce displacements and vibration amplitudes of the wind turbine blades. Therefore, structural geometrical nonlinearities cannot be neglected both in the static aeroelastic analysis and dynamic aeroelastic analysis.     

Supercritical and subcritical dynamic flow-induced instabilities of a small-scale wind turbine blade placed in uniform flow

There is a growing interest in extracting more power per turbine by increasing the rotor size in offshore wind turbines. As a result, the turbine blades will become longer and therefore more flexible, and a flexible blade is susceptible to flow-induced instabilities. In order to design and build stable large wind turbine blades, the onset of possible flow-induced instabilities should be considered in the design process. Currently, there is a lack of experimental work on flow-induced instabilities of wind turbine blades. In the present study, a series of experiments were conducted and flow-induced instabilities were observed in wind turbine blades. A small-scale flexible blade based on the NREL 5 MW reference wind turbine blade was built using three-dimensional printing technique. The blade was placed in the test section of a wind tunnel and was subjected to uniform oncoming flow, representing the case of a parked wind turbine blade. The blade׳s tip displacement was measured using a non-contacting displacement measurement device as the oncoming wind speed was increased. At a critical wind speed, the blade became unstable and experienced limit cycle oscillations. The amplitude of these oscillations increased with increasing wind speed. Both supercritical and subcritical dynamic instabilities were observed. The instabilities were observed at different angles of attack and for blades both with and without a geometric twist. It was found that the blade twist had a significant influence on the observed instability: a blade without a twist experienced a strong subcritical instability.     

Wake impacts on aerodynamic and aeroelastic behaviors of a horizontal axis wind turbine blade for sheared and turbulent flow conditions

The magnitude and temporal variations of wind speed considerably influence aerodynamic and structural responses of MW-sized horizontal axis wind turbines. Thus, this paper investigates the variations in airloads and blade behavior of a wind turbine blade resulting from operations in sheared and turbulent flow conditions. First, in order to validate the present methods, comparisons of aerodynamic results were made among the blade element momentum method, free-wake method, and numerical results from the previous studies. Then, the validated methods were applied to a national renewable energy laboratory 5 MW reference wind turbine model for fluid–structure interaction analyses. From the numerical simulations, it can be clearly seen that unfavorable airloads and blade deformations occur due to the sheared and turbulent flow conditions. In addition, it is clear that wake impacts are not as substantial at those of high wind speeds; however, the effects obviously affect the aerodynamic and structural behaviors of the blade at lower wind speeds. Therefore, it is concluded that the numerical results markedly indicate the demand for accurate assessment of wake dynamics for accurate estimations of the aerodynamic and structural responses for sheared and turbulent flow environments.     

Criterion of aerodynamic performance of large-scale offshore horizontal axis wind turbines

With the background of offshore wind energy projects, this paper studies aerodynamic performance and geometric characteristics of large capacity wind turbine rotors （1 to 10 MW）, and the main characteristic parameters such as the rated wind speed, blade tip speed, and rotor solidity. We show that the essential criterion of a high- performance wind turbine is a highest possible annual usable energy pattern factor and a smallest possible dimension, capturing the maximum wind energy and producing the maximum annual power. The influence of the above-mentioned three parameters on the pattern factor and rotor geometry of wind turbine operated in China＇s offshore meteoro- logical environment is investigated. The variation patterns of aerodynamic and geometric parameters are obtained, analyzed, and compared with each other. The present method for aerodynamic analysis and its results can form a basis for evaluating aerodynamic performance of large-scale offshore wind turbine rotors.     

Design and verification of airfoil families for large-size wind turbine blades

For the design of wind turbine blades, the use of a family of specially tailored airfoils is particularly important. The dedicated airfoils can dramatically improve the capability of capturing wind power, reduce the structural weight to save the cost of manufacturing and transportation, and lower the inertial loads as well as the loads due to gust. An overview of the world-wide wind turbine airfoil families developed since 1990's is presented, such as the S series, the DU series, the Risø series, and the FFA series. The design and wind-tunnel tests of the Northwestern Polytechnical University (NPU) airfoil family for megawatt-size wind turbines, called the NPU-WA series, are summarized. All tests for the NPU-WA series are carried out in the NF-3 low-speed wind-tunnel with a two-dimensional (2D) test section of 1.6m×0.8m and at the Reynolds number ranging from 1.6×10⁶ to 5×10⁶. The research activities for further improving the NPU-WA airfoils towards lower roughness sensitivity are also reviewed. The development of the new NPUWA series dedicated for multi-megawatt wind turbines is discussed.     

风速不确定性对风力机气动力影响量化研究

风力机气动力学一直是国内外研究的热点课题之一.目前相关研究大都是基于确定性工况条件, 但因风力机常年工作在自然来流复杂环境,风速随机波动致使风电系统呈现不确定性, 对电网稳定性带来巨大挑战,因此进行不确定风速条件下风力机气动力学研究具有重要意义.为揭示不确定性对风力机流场影响机理并明确其对气动力的影响程度,本文提出一种风力机不确定空气动力学分析方法,基于修正叶素动量理论和非嵌入式概率配置点法,建立水平轴风力机不确定性空气动力学响应模型; 以NREL Phase VI S809风力机叶轮为研究对象, 基于该模型提取风力机输出随机响应信息,量化不确定风速对风力机风轮功率、推力、叶片挥舞弯矩和摆振弯矩的影响程度;通过分析流动诱导因子不确定性在叶片展长方向上的分布规律,揭示不确定因素在风力机本体上的传播机制,为风电系统设计及应用提供理论依据和重要参考. 结果表明,风速波动对风力机功率和气动力影响显著,高斯风速标准差由0.05倍增大至0.15倍均值,功率和推力最大波动幅度分别由13.44%和8.00%增大至35.11%和22.02%,叶片挥舞弯矩和摆振弯矩最大波动幅度分别由7.20%和12.84%增大至19.90%和33.49%.来流风速不确定性导致叶片根部位置气流明显波动,可以考虑在该部分采取流动控制措施降低叶片对风速不确定性的敏感程度.