CN101458735A - 一种具有高升阻比的翼型 - Google Patents

一种具有高升阻比的翼型 Download PDF

Info

Publication number
CN101458735A
CN101458735A CNA2008102370169A CN200810237016A CN101458735A CN 101458735 A CN101458735 A CN 101458735A CN A2008102370169 A CNA2008102370169 A CN A2008102370169A CN 200810237016 A CN200810237016 A CN 200810237016A CN 101458735 A CN101458735 A CN 101458735A
Authority
CN
China
Prior art keywords
aerofoil profile
airfoil
molded lines
drag ratio
profile
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CNA2008102370169A
Other languages
English (en)
Inventor
陈进
王旭东
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chongqing University
Original Assignee
Chongqing University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Chongqing University filed Critical Chongqing University
Priority to CNA2008102370169A priority Critical patent/CN101458735A/zh
Publication of CN101458735A publication Critical patent/CN101458735A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Wind Motors (AREA)

Abstract

本发明涉及基于泛函的通用翼型型线形状设计的翼型,通过建立翼型型线的泛函集成方程,然后对该集成方程建立优化设计模型,根据所要求的目标函数选择不同的方程项数和系数来完成目标翼型的设计。本发明提供的基于泛函的通用翼型型线的高升阻比翼型与传统的风力机翼型相比,在正常工作功角范围内,不仅升力系数高,而且具有很高的升阻比,失速较晚。它的构成型线不但能够发挥通用翼型型线的优势,其构成的新型型线还能满足单一型线不能达到的技术指标,从而提高了风力机的输出功率,降低了风力发电成本。

Description

一种具有高升阻比的翼型
技术领域
本发明涉及风力机叶片翼型的设计,具体涉及一种具有高升阻比的翼型。它除了适应于风力机叶片外,也于适用于飞机机翼的设计制造。
背景技术
风力机叶片翼型设计理论是决定风力机功率特性和载荷特性的根本因素,一直是各国学者研究的热点所在。目前叶片翼型的研究和设计都是基于特定的几何轮廓曲面,来研究其空气动力学特性,并进行参数的设计和优化。风力机专门翼型是在航空翼型的基础上发展起来的,已经研究出的航空叶片翼型有美国航空航天局的NACA230系列、NACA44系列、NACA63-2系列、NACALS系列和美国贝尔公司的FX系列等;专门的风力机翼型有美国可再生能源实验室的NREL S系列、瑞典FFA-W系列、丹麦Ris
Figure A200810237016D0004092706QIETU
系列和荷兰DU系列等。
翼型的性能对风力机的气动性能具有决定性的影响,高性能翼型的研究是风力机发展的一项基础性研究。高升力、低阻力的翼型一直是翼型设计中所追求的目标。对于翼型的气动外形设计目前通常采用反设计方法,这种方法是根据给定希望达到的气动特性,压力分布以及初始的基本翼型,通过几何和流动控制方程,逐步逼近给定的气动特性。这种设计方法不够直接,设计过程也比较耗时,而且往往逼近的结果不够准确。因此,直接从翼型的形状表达出发,寻找一种更为直接高效的设计方法,就显得尤为必要。
此外,由于特定的几何曲线数学模型一经确定,其固有的几何特性和数学特性无法变更,因而性能受到根本制约。对于翼型的设计如果只基于给定的翼型形状,进行参数的优化、改进,并不能对表征叶片翼型型线本质特征的根本因素—翼型型线的形函数本身进行优选,因而难以在叶片翼型型线的研究和创新中取得根本的突破。
因此,研究风力机叶片翼型型线应具有怎样的函数特性,该函数的数学模型如何表达,该数学模型表征的几何曲线具有怎样的空气动力学特性,如何对该几何曲线进行形状优化等根本问题,对于建立风力机叶片翼型型线的完整理论,并根据风电场的实际情况构造出各种具有最佳性能的叶片翼型型线,具有重大的理论意义和工业应用价值。
综上所述,本专利通过研究翼型的通用特性,建立其泛函集成方程,进而提出了一种全新的、高效的翼型优化设计方法。翼型采用泛函集成方程来表达,其几何形状及解析特性可通过级数系数的调整和优化加以控制;形状优化的泛函极值问题也可以通过对可变项数级数系数的参数优化来逐次逼近。进而实现各种翼型形状的构造。
发明内容
本发明的目的就是提供了一种基于泛函的通用翼型型线形状设计的翼型。
本发明以上目的是通过以下技术方案来实现的,首先建立翼型型线的泛函集成方程,然后对该集成方程建立优化设计模型,根据所要求的目标函数选择不同的方程项数和系数来完成目标翼型的设计。
建立笛卡儿坐标系下z平面上的一个圆,通过改变圆心的位置,利用儒科夫斯基变换式:
ζ=f(z)=z+a2/z             (1)
就能够变换成ζ平面上的一个翼型,其中,a为1/4翼型弦长。
此时,ζ平面上的翼型在笛卡儿坐标系下就可以表示为:
ξ = ( r + a 2 / r ) cos θ η = ( r - a 2 / r ) sin θ - - - ( 2 )
其中,r为翼型的矢径长度。
为了表征翼型的通用性,r其实是一个关于θ的函数,表示为:
Figure A200810237016D00052
通过选取不同的就可以变换出各种不同性能的翼型。
翼型型线是由具有翼型特性的曲线组成的。根据以上的儒科夫斯基变换理论,基于Taylor级数思想的可取函数
Figure A200810237016D00054
的广义泛函方程可以集成为
Figure A200810237016D00055
k=1,2,3,...n         (4)
可以看出该式满足
Figure A200810237016D00056
保证了翼型的尖后缘特性。
对风力机翼型而言,其升力与阻力比是最重要的性能指标,对提高风轮效率和对整个风力机组工作性能的影响关系极大。在翼型设计过程中,一般期望得到较大的升力同时阻力又保持很小,所以将翼型的升阻比作为目标函数,建立广义翼型型线泛函极值模型:
f(x)=max(Cl/Cd)          (5)
其中Cl位翼型的升力系数,Cd为翼型的阻力系数。
对于翼型型线的泛函集成方程,选择函数
Figure A200810237016D00057
的前6项系数为优化设计变量,记作
X=(x1,x2,x3,x4,x5,x6)        (6)
式中变量的变化范围根据翼型的形状来确定。当各个变量的值超过一定范围时,该表征式将不再具有翼型的形状特征,所以对变量的范围进行了约束:
Xmin≤X≤Xmax。(7)
除了对形状的要求以外,其它的约束主要针对翼型的结构和强度。通常情况下,叶片主要产生功率区翼型的最大相对厚度为20%左右,选择翼型的最大相对厚度为X,即施加约束:
t/c=X            (8)
其中t为翼型的最大厚度,c为翼型的弦长。
翼型的结构和强度特性除了和最大相对厚度有关以外,还与翼型最大相对厚度所处的弦长位置有关,一般翼型的转矩中心在其弦长的1/4处左右,所以翼型的最大相对厚度在弦长所处的位置控制在
0.2≤x/c≤0.3。             (9)
基于以上的翼型型线优化设计模型,以型线集成方程系数为变量,通过改变、舍取级数项数,从而使型线得到优化,设计出具有高升阻比的翼型。
有效结果:本发明提供的基于泛函的通用翼型型线的高升阻比翼型与传统的风力机翼型相比,在正常工作功角范围内,不仅升力系数高,而且具有很高的升阻比,失速较晚。它的构成型线不但能够发挥通用翼型型线的优势,其构成的新型型线还能满足单一型线不能达到的技术指标,从而提高了风力机的输出功率,降低了风力发电成本,必将促进风电产业的大力发展。
附图说明
下面结合附图给出的一个非限定性的实施例对本发明作进一步说明:
图1——相对厚度18%的翼型形状;
图2——相对厚度18%的翼型升力和功角关系图;
图3——相对厚度18%的翼型升阻比关系图;
图1为最大相对厚度为18%的翼型,该翼型最大相对厚度在弦长的27%处。
具体实施方式
在实施例给出的翼型的设计过程中,根据TAYLOR级数思想和级数函数收敛原则取前6项,选择一定的相对厚度来进行翼型形状设计。这里给出了最大相对厚度为18%的翼型,表1给出了这种翼型泛函集成方程的系数,将这些系数代入到集成方程中就能得到翼型的形状,该翼型的最大相对厚度为18%,最大相对厚度在弦长的27%处。见图1。
表1 翼型泛函集成方程系数
 
最大相对厚度 X1 x2 x3 x4 x5 x6
18% 0.1132 0.1328 0.0072 0.0261 0.0001 0.0189
采用XFOIL软件对于本发明提出的这种翼型的气动性能进行了计算分析。雷诺数Re=1.6×106时,最大相对厚度为18%的本发明翼型的升阻比见表2,当功角为5°时,拥有最大升阻比161.8465,而且该翼型在风力机叶片正常工作功角0°~10°内均具有很大的升阻比。
表2雷诺数为Re=1.6×106的最大相对厚度为18%翼型的升阻比
 
角度(°) 0 1 2 3 4 5 6 7
升阻比 126.8564 138.9669 148.8447 156.1282 160.9059 161.8465 157.6797 147.1154
角度(°) 8 9 10 11 12 13 14 15
升阻比 138.2119 119.0821 98.9927 80.6389 63.9088 50.6272 40.6724 33.3490

Claims (2)

1、一种具有高升阻比的翼型,由具有翼型特性的型线平滑连接构成,其特征在于,具有翼型特性的型线按一下方式构建:
建立在笛卡儿坐标系下z平面上的一个圆,通过改变圆心的位置,利用儒科夫斯基变换式:
Figure A200810237016C0002085642QIETU
变换成
Figure A200810237016C0002161335QIETU
平面上的一个翼型,其中,a为1/4翼型弦长;
此时,
Figure A200810237016C0002161335QIETU
平面上的翼型在笛卡儿坐标系下表示为:
ξ = ( r + a 2 / r ) cos θ η = ( r - a 2 / r ) sin θ - - - ( 2 )
其中,r为翼型的矢径长度,是一个关于θ的函数,表示为:
Figure A200810237016C00022
取函数
Figure A200810237016C00023
的广义泛函方程可以集成为
Figure A200810237016C00024
k = 1,2,3 , . . . n - - - ( 4 )
该式满足
Figure A200810237016C00026
保证了翼型的尖后缘特性;
翼型族的形状由式(4)表示的通用翼型型线来控制,上述泛函表征式集成了所有具有翼型特性的光滑曲线;集成型线系数参数通过利用Taylor级数及其收敛原则,以该系数为变量,通过改变、舍取级数项数,最终优化设计得到了具有高升阻比的集成翼型系列。
2、根据权利要求1所述的具有高升阻比的翼型,其特征是:将翼型的升阻比作为目标函数,建立广义翼型型线泛函极值模型:
f(x)=max(Cl/Cd)              (5)
其中Cl位翼型的升力系数,Cd为翼型的阻力系数;
对于翼型型线的泛函集成方程(4),选择函数
Figure A200810237016C0002161406QIETU
(θ)的前6项系数为优化设计变量,记作
X=(x1,x2,x3,x4,x5,x6)         (6)
式中变量的变化范围根据翼型的形状来确定,当各个变量的值超过一定范围时,该表征式将不再具有翼型的形状特征,所以对变量的范围进行了约束:
Xmin≤X≤Xmax              (7)
其它的约束主要针对翼型的结构和强度,通常情况下,叶片主要产生功率区翼型的最大相对厚度为20%左右,选择翼型的最大相对厚度为X,即施加约束:
t/c=X                        (8)
其中t为翼型的最大厚度,c为翼型的弦长;
翼型的转矩中心在其弦长的1/4处左右,所以翼型的最大相对厚度在弦长所处的位置控制在
0.2≤x/c≤0.3。             (9)
基于以上的翼型型线优化设计模型,以型线集成方程系数为变量,通过改变、舍取级数项数,从而使型线得到优化,设计出具有高升阻比的翼型。
CNA2008102370169A 2008-12-31 2008-12-31 一种具有高升阻比的翼型 Pending CN101458735A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNA2008102370169A CN101458735A (zh) 2008-12-31 2008-12-31 一种具有高升阻比的翼型

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNA2008102370169A CN101458735A (zh) 2008-12-31 2008-12-31 一种具有高升阻比的翼型

Publications (1)

Publication Number Publication Date
CN101458735A true CN101458735A (zh) 2009-06-17

Family

ID=40769592

Family Applications (1)

Application Number Title Priority Date Filing Date
CNA2008102370169A Pending CN101458735A (zh) 2008-12-31 2008-12-31 一种具有高升阻比的翼型

Country Status (1)

Country Link
CN (1) CN101458735A (zh)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101813074A (zh) * 2010-03-25 2010-08-25 张玉福 简化的环形大气动力机
CN101615216B (zh) * 2009-07-28 2011-01-05 重庆大学 对翼型尾缘进行流线形增厚的方法
CN101923584B (zh) * 2009-10-30 2012-08-15 重庆大学 风力机专用翼型设计方法及风力机专用翼型
CN103136422A (zh) * 2013-01-11 2013-06-05 重庆大学 翼型集成与b样条结合的中等厚度翼型设计方法
CN104317992A (zh) * 2014-10-11 2015-01-28 重庆工商大学 风力机翼型正设计方法及风力机翼型族
CN104354875A (zh) * 2014-10-31 2015-02-18 中国航天空气动力技术研究院 一种考虑运动轨迹的增升装置优化设计方法
CN104401504A (zh) * 2014-11-19 2015-03-11 中国地质大学(武汉) 一种固定翼航测型无人机设计方法
CN104612892A (zh) * 2014-12-30 2015-05-13 中国科学院工程热物理研究所 一种风力机翼型的多学科优化设计方法
CN104834772A (zh) * 2015-04-22 2015-08-12 复旦大学 基于人工神经网络的飞机翼型/机翼反设计方法
CN105822510A (zh) * 2016-04-12 2016-08-03 南京航空航天大学 一种叶尖系列化分段叶片及其设计方法
CN108331708A (zh) * 2018-01-25 2018-07-27 华北电力大学 风力发电叶片翼型线的确定方法
CN111859801A (zh) * 2020-07-16 2020-10-30 湖北工业大学 失速型风力机翼型的设计方法

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101615216B (zh) * 2009-07-28 2011-01-05 重庆大学 对翼型尾缘进行流线形增厚的方法
CN101923584B (zh) * 2009-10-30 2012-08-15 重庆大学 风力机专用翼型设计方法及风力机专用翼型
CN101813074A (zh) * 2010-03-25 2010-08-25 张玉福 简化的环形大气动力机
CN103136422B (zh) * 2013-01-11 2015-12-02 重庆大学 翼型集成与b样条结合的中等厚度翼型设计方法
CN103136422A (zh) * 2013-01-11 2013-06-05 重庆大学 翼型集成与b样条结合的中等厚度翼型设计方法
CN104317992A (zh) * 2014-10-11 2015-01-28 重庆工商大学 风力机翼型正设计方法及风力机翼型族
CN104354875A (zh) * 2014-10-31 2015-02-18 中国航天空气动力技术研究院 一种考虑运动轨迹的增升装置优化设计方法
CN104354875B (zh) * 2014-10-31 2016-08-24 中国航天空气动力技术研究院 一种考虑运动轨迹的增升装置优化设计方法
CN104401504A (zh) * 2014-11-19 2015-03-11 中国地质大学(武汉) 一种固定翼航测型无人机设计方法
CN104612892B (zh) * 2014-12-30 2017-03-29 中国科学院工程热物理研究所 一种风力机翼型的多学科优化设计方法
CN104612892A (zh) * 2014-12-30 2015-05-13 中国科学院工程热物理研究所 一种风力机翼型的多学科优化设计方法
CN104834772A (zh) * 2015-04-22 2015-08-12 复旦大学 基于人工神经网络的飞机翼型/机翼反设计方法
CN104834772B (zh) * 2015-04-22 2018-01-12 复旦大学 基于人工神经网络的飞机翼型/机翼反设计方法
CN105822510A (zh) * 2016-04-12 2016-08-03 南京航空航天大学 一种叶尖系列化分段叶片及其设计方法
CN108331708A (zh) * 2018-01-25 2018-07-27 华北电力大学 风力发电叶片翼型线的确定方法
CN111859801A (zh) * 2020-07-16 2020-10-30 湖北工业大学 失速型风力机翼型的设计方法
CN111859801B (zh) * 2020-07-16 2022-07-19 湖北工业大学 失速型风力机翼型的设计方法

Similar Documents

Publication Publication Date Title
CN101458735A (zh) 一种具有高升阻比的翼型
CN101923584B (zh) 风力机专用翼型设计方法及风力机专用翼型
CN103136422B (zh) 翼型集成与b样条结合的中等厚度翼型设计方法
CN101615216B (zh) 对翼型尾缘进行流线形增厚的方法
CN109902384B (zh) 一种基于气弹模型的风力机柔性叶片预弯预扭设计方法
CN104405596A (zh) 一种风力发电机组低风速翼型族
CN106777567A (zh) 一种风力机钝尾缘翼型优化设计方法
Gupta et al. Modeling and aerodynamic analysis of small scale, mixed airfoil horizontal axis wind turbine blade
CN105404743A (zh) B样条与曲率光滑连续性结合的风力机翼型设计方法
CN104317992A (zh) 风力机翼型正设计方法及风力机翼型族
CN102094767A (zh) 一族用于兆瓦级风力机叶片的翼型
CN109409019A (zh) 一种风力机叶片气动外形的预弯优化方法
CN102444540A (zh) 一种水平轴风力发电机的风轮叶片翼型
CN106777525A (zh) 考虑设计叶尖速比风轮静、动态影响的风力机气动设计方法
CN104018998B (zh) 一种用于兆瓦级风力机叶片的21%厚度主翼型
CN112211779A (zh) 一种基于生物特征的仿生风力发电机叶片
CN103939283A (zh) 垂直轴风力机专用高效叶片
Khakpour Nejadkhaki et al. A Design Methodology for a Flexible Wind Turbine Blade With an Actively Variable Twist Distribution to Increase Region 2 Efficiency
CN101886619A (zh) 风力发电机叶尖专用翼型
CN106741924A (zh) 一种全工况条件下的高升阻比低力矩特性9%厚度旋翼翼型
Halder et al. An Improved Aerodynamic Modeling of the Wind Turbine for the Wind Power Harvesting
El-Okda et al. Design of a small horizontal axis wind turbine
Méndez et al. Airfoil family design for large offshore wind turbine blades
CN202348553U (zh) 一种水平轴风力发电机的风轮叶片翼型
CN112065651B (zh) 一种用于风力发电机组的风轮叶片层的翼型

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C12 Rejection of a patent application after its publication
RJ01 Rejection of invention patent application after publication

Open date: 20090617