法律状态公告日
法律状态信息
法律状态
2015-03-04
授权
授权
2013-10-09
实质审查的生效 IPC(主分类):F03D1/06 申请日:20130614
实质审查的生效
2013-09-04
公开
公开
技术领域
本发明属于水平轴风力机翼型设计领域,具体涉及适用于大型水平轴风 力机叶片根部的一族厚翼型轮廓及其设计方法。
背景技术
研发出性能优异的风力机专用翼型族是提高水平轴风力机叶片能量捕获 率的关键,而风力机翼型的研发目标与风力机的发展需求密切相关。近年来, 中国风电产业迅猛发展,各种规模的风场在不同地区建立,新装风力机的容 量不断增加,叶片尺寸显著变长。
首先,不同的风资源环境对风力机叶片要求不同。中国南方广大地区风 资源困乏,北方地区和沿海地区属于风资源丰富地区但却存在各种不利的条 件。华北和西北地区风沙严重,极易对风力机叶片表面造成污染和侵蚀。沿 海地区空气湿度较大且热带气旋濒临频繁。一方面温湿的空气容易对叶片表 面造成腐蚀破坏其几何轮廓;另一方面,频繁的热带气旋侵犯对风力机叶片 的结构强度是一个巨大考验。这要求设计出在变工况条件下性能优异且稳定 的翼型系列。
其次,风力机的叶片日趋大型化对叶片的材料和结构特性提出了更高要 求;同时,叶片尺寸的显著增长,风切变效应的影响更加明显,加之在风力 机运行过程中叶片不同部位的风速呈线性分布使得叶片各部位翼型的流动雷 诺数差异很大。这些要求开发出性能优良且随着Re变化稳定的适用于叶片中 部和尾部的大厚度翼型系列。
另一方面,增加翼型厚度及尾缘厚度可以增加翼型的横截面积、改善翼 型的厚度分布,显著提高翼型的结构刚度。同时适当提高翼型的尾缘厚度可 以减小吸力面逆压梯度,从而延迟湍流分离,提高最大升力系数。尽管阻力 系数也会增加,但翼型的升阻比是增加的,所以翼型的气动效率得到了加强。 通常,大厚度翼型适用于叶片根部,该处流动常常处于大攻角的流动分离流 动状态,具有较强的三维流动特性,因此性能优良的大厚度、钝尾缘翼型的 设计,对于研发高性能大型风电叶片具有重要意义。
目前应用于风力机叶片内侧的翼型的相对厚度多为30%~40%,尾缘厚度 也较小。而且现有厚翼型的设计雷诺数单一;设计攻角较小,远远低于风力 机运行时根部区域翼型所处的攻角范围;升力系数失速剧烈,在大攻角下较 小。这些缺陷不利于提高大型风力机的风能捕获率。
发明内容
发明目的:鉴于以上问题,本发明目的在于提出一族在运行攻角下具有 高升力系数且变工况性能稳定的大厚度钝尾缘翼型,可以满足大型风力机叶 片在结构和气动方面对翼型的需求。
本发明是通过以下技术方案实现的。
一族大厚度钝尾缘风力机翼型,在运行攻角下具有高升力系数且变工况 性能稳定,包括四种大厚度翼型,其特征在于,所述四种大厚度翼型的最大 相对厚度依次为45%、50%、55%和60%,尾缘相对厚度依次为7%、9%、12% 和16%;所述四种翼型的最大厚度的弦向相对位置位于30%~35%之间,上表 面最大厚度的弦向相对位置在28%~33%之间,下表面最大厚度的弦向相对位 置在33%~37%之间;所述四种翼型的中弧线,均在弦线的上面,且在靠近前 缘处和靠近尾缘处分别有一个向上的凸起峰;所述四种翼型的最大相对弯度 值在2%~3.2%之间,最大弯度的弦向相对位置在74%~79%之间;所述四种 翼型的前缘半径在0.099~0.125倍弦长之间;其中,所述最大相对厚度是各翼 型上下两个面之间的最大厚度与弦长的比值,各所述弦向相对位置是指弦向 位置相对弦长的比值。
优选地,所述四个翼型的中弧线沿着弦长的分布特征一致,靠近前缘处 的中弧线向上凸起峰,其弦向相对位置在22%~23%之间;靠近尾缘处的中弧 线向上凸起峰,其弦向相对位置在75%~78%之间。
优选地,所述四种大厚度钝尾缘翼型用于风力机叶片的根部区域,在运 行攻角为15°~30°的范围内,具有高升力系数和稳定变工况性能。
优选地,所述四种大厚度钝尾缘翼型的设计雷诺数依次为4.0×106、3.5 ×106、3.0×106、2.5×106。
优选地,45%相对厚度的翼型的设计雷诺数为4×106,该翼型吸力面和 压力面的无量纲几何坐标分别为:
表一45%相对厚度翼型的吸力面坐标(x/c,y/c)
其中,x/c值表示翼型吸力面上某点在弦线方向上相对于前缘的位置,y/c 值表示从弦线到翼型吸力面上某点的高度;
表二45%相对厚度翼型的压力面坐标(x/c,y/c)
其中,x/c值表示翼型压力面上某点在弦线方向上相对于前缘的位置, y/c值表示从弦线到翼型压力面上某点的高度。
进一步地,45%相对厚度的翼型,其失速攻角约为8°,失速后,升力系 数变化平缓,在8°~15°攻角范围内升力系数先略下降后上升,失速特性参 数为Mstall=85.56×10-5。失速特性参数的定义式如下:
其中,cl为升力系数,clmax为最大升力系数,α为攻角,αstall为失速攻角。
进一步地,所述45%相对厚度的翼型,在15°~30°大攻角范围内,翼型 的升力系数从1.6开始略下降,然后开始稳定上升,在30°时达到1.7387。
优选地,50%相对厚度翼型的设计雷诺数为3.5×106,该翼型吸力面和压 力面的无量纲几何坐标分别为
表三50%相对厚度翼型的吸力面坐标(x/c,y/c)
表四50%相对厚度翼型的压力面坐标(x/c,y/c)
进一步地,50%相对厚度翼型,其失速攻角为5°,升力系数第一次达到 最大的数值为1.2258,失速后,升力系数在5°~15°攻角范围内先下降后上 升,整体升力系数值维持在1.08~1.22,失速特性参数为Mstall=101.25×10-5。 进一步地,所述50%相对厚度翼型,在15°~31°大攻角范围内,升力系数从 1.17开始稳定上升,在31°时达到1.6209。
优选地,55%相对厚度的翼型的设计雷诺数为3×106,该翼型吸力面和 压力面的无量纲几何坐标分别为:
表五55%相对厚度翼型的吸力面坐标(x/c,y/c)
表六55%相对厚度翼型的压力面坐标(x/c,y/c)
进一步地,55%相对厚度翼型,其失速攻角为5°,升力系数第一次达到 最大的数值为1.1543,失速后的升力系数在5°~15°攻角范围内下降到约1 后,平稳维持在1.02~1.08水平,失速特性参数为Mstall=46.41×10-5。进一步 地,所述55%相对厚度翼型,在15°~30°大攻角范围内,升力系数从1.0308 以近似线性变化的方式上升,在30°时达到1.6797。
优选地,60%相对厚度翼型的设计雷诺数为2.5×106,该翼型吸力面和压 力面的无量纲几何坐标分别为:
表七60%相对厚度翼型的吸力面坐标(x/c,y/c)
表八60%相对厚度翼型的压力面坐标(x/c,y/c)
进一步地,60%相对厚度的翼型,其失速攻角约为2°,最大升力系数值 为0.8165,失速后,升力系数在2°~15°攻角范围内略下降到0.63左右后稳 定上升,失速平缓特性为Mstall=179.17×10-5。进一步地,所述60%相对厚度 的翼型,在14°~30°大攻角范围内,升力系数从1.0779以近似线性变化的方 式上升,在30°时达到1.783。
根据本发明的另一方面,还提供了一种用于设计上述翼型的设计方法, 包括正设计步骤和反设计步骤,其特征在于,所述正设计步骤中,从原始翼 型出发,首先保持原始翼型的弯度分布和最大厚度位置基本不变,调节翼型 的相对厚度达到目标相对厚度;之后保持翼型的相对厚度不变,调节翼型的 尾缘厚度达到目标尾缘厚度;完成所述正设计步骤后,进行反设计步骤,所 述反设计步骤中,通过调节翼型表面的压力分布,来微调翼型的前缘形状、 相对厚度及其分布、相对弯度及其分布。
优选地,正设计步骤中,所述调节翼型的尾缘厚度达到目标尾缘厚度, 通过保持相对厚度不变,沿中弧线对称放大尾缘的方式实现。优选地,所述 调节翼型表面的压力分布,主要是对翼型的边界层内前缘吸力面压力峰的高 低、压力峰的弦向位置和压力峰回落的梯度等三个主要因素进行调节。
有益结果:本发明提供的一族在运行攻角下具有高升力系数的大厚度钝 尾缘翼型可以弥补大型风力机叶片设计时面临的厚翼型族的匮乏。大的相对 厚度和尾缘厚度可以有效的改善叶片的结构属性。更重要的是,该翼型族的 设计目标明确把叶片实际运行时翼型所处的各自雷诺数条件、大攻角下的升 力系数作为气动目标,而不是像先前把单一雷诺数、较小攻角下的最大升阻 比和最大升力系数做为气动目标,更贴近大型风力机叶片根部叶片的性能需 求。
附图说明
图1为本发明所提出的大厚度钝尾缘翼型族的基本型线;
图2为本发明所提出的各翼型的中弧线分布曲线图;
图3为本发明所提出的各翼型在各自设计雷诺数下的升力系数曲线;
图4为相对厚度为45%的翼型升力系数随Re变化的曲线图;
图5为相对厚度为50%的翼型升力系数随Re变化的曲线图;
图6为相对厚度为55%的翼型升力系数随Re变化的曲线图;
图7为相对厚度为60%的翼型升力系数随Re变化的曲线图;
图中,45%代表相对厚度为45%的翼型;50%代表相对厚度为50%的翼型; 55%代表相对厚度为55%的翼型;60%代表相对厚度为60%的翼型。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚明白,以下参照附图并举 实施例,对发明进一步详细说明。
本发明所提出的翼型是采用正反方法相结合的方式得到。从原始翼型出 发,采用正设计方法,首先保持翼型的弯度分布和最大厚度位置基本不变, 设计得到四种相对厚度的翼型;进一步实现厚尾缘造型设计。厚尾缘的设计 有两种主要的思路:直接截断原翼型的尾部区域或保持相对厚度不变,对称 放大尾缘。直接截断尾缘会显著的改变翼型的厚度和弯度分布,对翼型的气 动性能影响极大,尤其会使得翼型的升力系数退化。在此采用后一种造型方 法,保持相对厚度和最大厚度位置基本不变,对称放大尾缘,可以有效降低 边界逆压梯度,改善翼型的气动性能。这样,经过初始造型设计得到符合相 对厚度、尾缘厚度以及最大相对厚度位置基本预期范围的初始翼型。
在正设计中,保持相对厚度不变,提高翼型的尾缘厚度可以降低某一流 动条件下的压力梯度,延迟边界层的分离,获得较高的最大升力系数。反设 计这一阶段进一步以气动性能的优化为核心,修改翼型的压力分布特征,来 改进翼型的气动参数。大厚度翼型在攻角很小时就会发生转捩,(或许会伴 随着分离泡的再附)形成湍流边界层;湍流边界层较难分离,在尾缘厚度较 大时,分离点更靠近尾缘。对翼型气动性能影响最显著的压力分布特征是边 界层内前缘吸力面压力峰的高低、压力峰的弦向位置和压力峰回落的梯度。 通过对这三者的修改,来间接地细微地调节翼型的前缘形状、相对厚度及其 分布、相对弯度及其分布,以降低翼型的流动分离点随着攻角的增加自尾缘 向前缘移动的速度,从而使得翼型在失速后趋于变化平缓,甚至伴随一个升 力系数继续增加的过程。大厚度翼型在大攻角时通常处于完全湍流状态,受 Re的影响较小,具有较高的Re稳定性。如此反复优化,得到符合设计预期 的四种翼型,如图1所示。
图1为本发明所提出的大厚度钝尾缘翼型族的基本型线,各翼型的基本 几何参数如表九所示,对应的四种翼型的无量纲型线坐标如表一至表八所示. 主要参数为翼型的具体相对厚度、最大相对厚度位置、相对弯度、最大弯度 位置、前缘半径、尾缘厚度等(数值均为相对于弦长的比较值)。其中,X(Y+)、 X(Y-)分别表示上下表面最大厚度位置的弦向位置。
表九一族大厚度钝尾缘翼型的基本几何参数
四种大厚度翼型的最大厚度在弦长方向的相对位置位于30%~35%之间, 具体说来,上表面最大厚度位置在(0.28~0.33)倍弦长间,下表面最大厚度 位置在(0.33~0.37)倍弦长之间;如图2所示,四个翼型的中弧线沿着弦长 的分布特征一致,均在弦长的上面,也就说翼型上表面较厚;中弧线在靠近 前缘出和尾缘处均有一个向上凸起峰,弦向位置分别在(0.22~0.23)倍弦长 和(0.75~0.78)倍弦长之间;最大相对弯度值相近,在(0.020~0.032)弦长 之间,最大弯度的弦向位置在74%~79%之间;四种翼型的前缘半径依次在 (0.099~0.125)弦长之间,尾缘厚度依相对厚度的增加渐次增加;这些有效 保证了翼型族的几何兼容性。
本发明所提出的大厚度钝尾缘翼型是尤其适用于大型风力机叶片。根据 我国风资源分布特征和大型风力机叶片(5MW风力机)各部分雷诺数的分布, 拟定45%到60%四种大厚度翼型的四种设计雷诺数分别为:4.0×106、3.5× 106、3.0×106、2.5×106。
另一方面,对于5MW容量的风力机叶片而言,45%及以上厚度的翼型基 本上适用于叶片展向位置的10%~30%位置处。由于该位置叶片的最大扭角是 受到限制的,所以造成最大弦长位置以内尤其是根部区域截面翼型的实际流 动攻角过大(5MW的风力机叶片根部翼型的攻角范围为12°~25°)。对于 大厚度翼型而言,在这一攻角下,升力曲线通常处于失速后区域,流动边界 层也通常为湍流边界层,阻力系数异常大。因此,本次发明中,以最大升力 系数、大攻角范围(15°~30°)内升力系数为主要气动目标,同时限制升力 系数要随着雷诺数Re的变化稳定。
对本发明所提出的翼型族的气动特征的说明均是基于粘性-无粘迭代的 涡面元方法(如本领域中所熟知的RFOIL)的数值计算结果,计算条件为各 自设计雷诺数下的自然转捩情况。表二给出了翼型族的关键气动参数。其中 最大升力系数为翼型升力系数第一次取得的最大值。翼型的失速平缓特性参 数所考察的攻角范围在这里是从失速攻角到15°攻角。由于在大攻角下,翼 型的升力系数再次以近似线性的方式增加(基于数值计算结果),翼型在大 攻角下的升力特征主要由升力的变化区间来表征。表中给出了15°和30°攻 角下各自的升力系数。由于失速攻角较小,而所考察的是大攻角范围内升力 系数随着Re变化的稳定性,所以在此只用不同雷诺数下升力曲线段的对比 来定性说明。图3为本发明所提出的翼型族数值求解得到的自然转捩条件各 自设计雷诺数下的升力系数曲线。
失速特性参数的定义式如下:
表十各翼型的基本气动参数
结合图3和表十,说明本发明所提出的大厚度钝尾缘翼型族的气动性能 特征。
45%相对厚度翼型的设计雷诺数为4×106。该翼型的失速攻角较高 (8°),因此升力系数第一次达到最大时候的数值在四种翼型中也最大 (1.6427)。失速后,翼型的升力系数变化较为平缓,在(8°~15°)攻角范 围内有一个略下降后上升的过程,但整体升力系数值维持在高水平,失速特 性参数仅为Mstall=85.56×10-5。在大攻角范围内(15°~30°),翼型的升力系 数从1.6开始略下降,然后开始一个稳定的上升过程,在30°时达到1.7387.
50%相对厚度翼型的设计雷诺数为3.5×106。随着相对厚度的增加,该翼 型自然转捩条件下的失速攻角提前为5°,升力系数第一次达到最大的数值 为1.2258.失速后的升力系数在(5°~15°)攻角范围内有一个下降后上升的 过程,但整体升力系数值维持在(1.08~1.22),失速特性参数为 Mstall=101.25×10-5。在大攻角范围内(15°~31°),翼型的升力系数从1.17开 始稳定上升,在31°时达到1.6209.
55%相对厚度翼型的设计雷诺数为3×106。该翼型自然转捩条件下的失 速攻角也为5°,升力系数第一次达到最大的数值为1.1543.失速后的升力系 数在(5°~15°)攻角范围内下降到1左右后变化极其平稳,维持在1.02~1.08 水平,所以失速平缓特性仅为Mstall=46.41×10-5。在大攻角范围内(15°~30°), 翼型的升力系数从1.0308以近似线性变化的方式上升,在30°时达到1.6797.
60%相对厚度翼型的设计雷诺数为2.5×106。由于相对厚度和尾缘厚度过 大,该翼型自然转捩条件下的升力系数在攻角很小时候(2°)就第一次达到 了最大,值为0.8165。失速后的升力系数在(2°~15°)攻角范围内略下降 到0.63(4°)左右后稳定上升,失速平缓特性为Mstall=179.17×10-5。在大攻 角范围内(14°~30°),翼型的升力系数从1.0779以近似线性变化的方式上 升,在30°时达到1.783.
综上,本发明所提出的翼型族第一次达到最大值后经历一个平缓的失速 过程,升力系数继续增加。在大攻角范围内,升力系数持续稳定的上升并处 于较高水平。在15°到30°攻角之间,该翼型族的升力系数大致在1.0到1.7 之间近似线性的变化(45%厚度翼型性能优异,在这一范围内的升力系数为 1.3812~1.7387)。所以达到了在大攻角范围内具有较高升力系数的气动目标。
在翼型的升力系数随Re变化的稳定性方面,依据数值计算结果,图4 展示了45%翼型的升力曲线随着Re的变化,图5展示了50%翼型的升力曲 线随着Re的变化,图6展示了55%翼型的升力曲线随着Re的变化,图7则 展示了60%翼型的升力曲线的变化。从图中可以明显的看出,本发明提出的 这组大厚度钝尾缘翼型在运行攻角范围内具有优异的雷诺数稳定性,有利于 风力发电机获得更稳定的输出。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本 发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在 本发明的范围之内。
机译: 绕射流控制的机翼钝尾缘周围的多重喷射
机译: 飞机起升表面的尾缘和尾缘的支承肋的制造方法
机译: 飞机起升表面的尾缘和尾缘的支承肋的制造方法