CN107145677B - 一种改进的几何参数翼型设计方法 - Google Patents

一种改进的几何参数翼型设计方法 Download PDF

Info

Publication number
CN107145677B
CN107145677B CN201710362161.9A CN201710362161A CN107145677B CN 107145677 B CN107145677 B CN 107145677B CN 201710362161 A CN201710362161 A CN 201710362161A CN 107145677 B CN107145677 B CN 107145677B
Authority
CN
China
Prior art keywords
airfoil
camber
thickness
expression
parameters
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.)
Active
Application number
CN201710362161.9A
Other languages
English (en)
Other versions
CN107145677A (zh
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.)
Beihang University
Original Assignee
Beihang 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 Beihang University filed Critical Beihang University
Priority to CN201710362161.9A priority Critical patent/CN107145677B/zh
Publication of CN107145677A publication Critical patent/CN107145677A/zh
Application granted granted Critical
Publication of CN107145677B publication Critical patent/CN107145677B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Automation & Control Theory (AREA)
  • Structures Of Non-Positive Displacement Pumps (AREA)

Abstract

本发明提出了一种改进的几何参数翼型设计方法,将表征翼型轮廓的弯度、厚度分开表达,并最终由厚度表达式和弯度表达式叠加而构建翼型函数。本发明的设计方法,采用较少的几何参数表征翼型,能够极大的减小计算机硬件系统的计算量,加快优化进程;同时,该设计方法尽量保持了优化空间的连续性,进而显著提高了翼型设计的精度和广度。

Description

一种改进的几何参数翼型设计方法
技术领域
本发明属于航空飞行器气动外形设计技术,特别是涉及一种航空飞行器翼型参数化设计方法。
背景技术
在飞机的各种飞行状态下,机翼是飞机承受升力的主要部件。一般飞机都有对称面,如果平行于对称面在机翼展向任意位置切一刀,切下来的机翼剖面称作为翼剖面或翼型。翼型是机翼和尾翼成形重要组成部分,其直接影响到飞机的气动性能和飞行品质。因此,在飞机的设计过程中,翼型的参数化设计显得尤为重要,其是完成气动优化的基础。
所谓参数化方法,即用有限的参数描述一个特定的外形,参数化方法的好坏将直接导致优化结果的优劣。目前,翼型参数化方法主要分为两大类,即变形法与描述法。对于某一特定形状的翼型而言,应用变形方法能比描述法获得更为准确的拟合效果。但当翼型设计空间较大时,描述法能够使用更少的控制参数以描述更多翼型。对于飞机概念设计阶段进行相对较粗的飞机外形设计优化而言,描述法显然具有更大的优势。
现有的翼型构建方面较常用的描述法主要有:基于特征的PARSEC几何参数法和基于函数的正交基函数法、CST法。基于特征的PARSEC几何参数法使用11个具有物理意义的特征参数描述翼型,该方法虽然有利于直观表示翼型的特征,但是其中涉及参数众多,大部分参数并没有实际意义,无法应用到翼型性能的分析中去,且难以覆盖更广的形状空间;基于函数的方法通过基函数的线性组合来表述翼型,如CST法等,其虽然可以覆盖更广的形状空间,但函数方法所获的参数数量在参数化阶段即已确定,在优化过程中无法动态调整参数对几何控制的精细化程度,而且多项式函数在特定的参数组合下会出现病态解,即该情况下的函数图形完全无法作为翼型使用。
发明内容
针对上述现有技术的不足,本发明提供了一种改进的几何参数翼型设计方法(Improved Geometric Parameters,下称IGP方法)。该方法相较于现有的上述描述法,采用8个优化参数表征翼型,并将翼型的弯度、厚度分开表达,极大地缩减了翼型设计空间的大小、减小了计算量、加快了优化进程。
本发明提供的改进的几何参数翼型设计方法,主要技术方案包括如下步骤:1)、获取与翼型轮廓相关的8个几何参数;2)、建立翼型的弯度表达式:
Figure GDA0002487426930000021
其中,xC和yC分别为翼型弯度上任一点的横坐标和纵坐标,c1、c2、c3、c4分别是三次贝塞尔曲线两个控制点的横、纵坐标,k为控制参数,取值范围为[0,1];3)、建立翼型的厚度表达式:t=t1x0.5+t2x+t3x2+t4x3+t5x4,其中,t为翼型在横坐标x处的厚度,t1、t2、t3、t4、t5为控制参数;4)、构建翼型表达式:该翼型表达式由所述弯度表达式和厚度表达式叠加而成,其采用厚度t和弯度xc,yc的表示式为:
Figure GDA0002487426930000022
其中,xu与yu为翼型上翼面曲线的横纵坐标,xl与yl为翼型下翼面曲线的横纵坐标。
本发明提供的改进的几何参数翼型设计方法还采用如下附属技术方案:
所述的8个几何参数为:相对弯度C,最大弯度所在的弦向位置XC,中弧线后缘夹角αTE,最大弯度处中弧线曲率bXC,相对厚度T,最大厚度所在的弦向位置XT,前缘半径ρ0,后缘夹角βTE
所述翼型表达式中的参数的优化设计空间为:
Figure GDA0002487426930000031
其中,
Figure GDA0002487426930000032
Figure GDA0002487426930000033
分别为ρ0与βTE的无量纲量。
所述控制参数t1、t2、t3、t4、t5可以通过以下方程组获得:
Figure GDA0002487426930000034
对于后缘厚度为0的标准化翼型,满足t(1)=0。
与现有技术相比,本发明至少存在以下技术效果:
第一、采用较少的几何参数构建翼型,能够以几何级数缩减翼型设计空间的大小,加快了概念设计阶段的计算速度;
第二、在拟合翼型过程中,将弯度、厚度拆开,减小了计算量,加快了优化进程;
第三、在优化过程中,保持了优化空间的连续性,提高了设计精度;
第四、设计中不需基础翼型,其控制参数本身也可直接与空气动力学理论中常用的翼型外形参数对应。
附图说明
图1为翼型函数几何参数示意图。
附图标记:附图中的标记说明,C-相对弯度,XC-最大弯度所在的弦向位置,αTE-中弧线后缘夹角,bXC-最大弯度处中弧线曲率,T-相对厚度,XT-最大厚度所在的弦向位置,ρ0-前缘半径,βTE-后缘夹角。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明中的改进的几何参数翼型设计方法(IGP方法)能够用于航空飞行器的翼型设计,该IGP方法能够使得本领域技术人员在概念设计阶段,运用位势流理论进行气动分析时简化计算。
参见附图1,本发明将翼型分解为有厚度的对称翼型和弯板翼型,首先获取翼型的8个几何参数:相对弯度C,最大弯度所在的弦向位置XC,中弧线后缘夹角αTE,最大弯度处中弧线曲率bXC,相对厚度T,最大厚度所在的弦向位置XT,前缘半径ρ0,后缘夹角βTE。例如可以通过计算机等具有采集、运算能力的硬件装置,或本领域其它熟知的手段获取上述几何参数。
弯度基于贝塞尔多项式进行表达,构建弯度表达式如式(1),其中,c1、c2、c3、c4分别是三次贝塞尔曲线两个控制点的横、纵坐标,k为控制参数,取值范围为[0,1],xc,yc分别为翼型模型弯度上任一点的横、纵坐标。
Figure GDA0002487426930000041
厚度则基于多项式基函数进行表达,构建厚度表达式如式(2),其中,t1、t2、t3、t4、t5为控制参数,x为翼型厚度横坐标,对于后缘厚度为0的标准化翼型,满足式(3)。
t=t1x0.5+t2x+t3x2+t4x3+t5x4 (2)
t(1)=0 (3)
构建翼型表达式:该翼型表达式由所述弯度表达式和厚度表达式叠加而成,其采用厚度t和弯度xc,yc的表示式为:
Figure GDA0002487426930000051
其中,xu与yu为翼型上翼面曲线的横纵坐标,xl与yl为翼型下翼面曲线的横纵坐标。
在已知翼型的上述几何参数的情况下,若要推得翼型的具体表达式,则需以上述8个几何参数作为约束,去求解8个优化参数。由于约束与未知数的个数相同,故方程组有唯一解。在实际计算中,由于弯度方程为参数方程,需引入最大弯度处的控制参数kc,同时存在控制方程
Figure GDA0002487426930000052
对其取值进行约束,使方程组的解仍唯一。
最终,经过变换,得到如下反推方程组:弯度反推方程组如式(4),厚度反推方程组如式(5):
Figure GDA0002487426930000053
Figure GDA0002487426930000061
而在实际的翼型构建过程中,若设计空间不连续,就会影响到随后优化流程中的绘图、气动计算等环节,使的优化过程中断。因此为了保证设计空间的连续性,本发明对部分优化参数进行了无量纲化处理,将“危险点”(即设计空间中存在“对应曲线形状过于奇怪”的点)置于取值范围的边缘处。
结合上述反推方程组(4)、(5),以及上述的无量纲化处理,求解出表征翼型的8个优化参数c1、c2、c3、c4、XT、T、
Figure GDA0002487426930000062
最终得到的翼型的参数设计空间,即:
Figure GDA0002487426930000063
Figure GDA0002487426930000064
其中,
Figure GDA0002487426930000065
Figure GDA0002487426930000066
分别为ρ0与βTE的无量纲量,具有如下关系:
Figure GDA0002487426930000067
Figure GDA0002487426930000068
根据上面的计算过程可知,根据方程组(5)、(7)可以将厚度表达式(2)中的控制参数t1、t2、t3、t4、t5采用优化参数XT、T、
Figure GDA0002487426930000071
表示,至此,将翼型表达式:
Figure GDA0002487426930000072
采用上述8个优化参数c1、c2、c3、c4、XT、T、
Figure GDA0002487426930000073
表征出来。
本发明的优点和有益技术效果:
1、总所周知,在翼型优化过程中,变量的个数增长给计算机带来的计算量增长是几何级数的,从而优化所需时间也以几何级数增长,这就要求在保证设计空间能够覆盖设计点的前提下,变量的个数尽量少。本发明在构建翼型函数表达式中采用的控制参数个数少于PARSEC几何参数法、正交基函数法、CST法三种描述法,能够以几何级数缩减翼型设计空间大小,极大地加快了概念设计阶段的计算机计算速度;
2、本发明的IGP方法在构建翼型过程中,将弯度、厚度拆开,对基于薄翼理论进行气动分析的优化问题而言,该IGP方法仅需使用4个参数构建翼型弯度,而其他描述法由于厚度弯度耦合,至少需要10个参数,因此本发明的该IGP方法大大减小了计算机系统的计算量,加快了优化进程;
3、本发明在优化过程中,保持了优化空间的连续性。而传统的多项式函数在特定的参数组合下会出现不连续点的病态解,使得构建的函数图形完全无法作为翼型使用。
尽管己描述了本发明的优选实施例,但本领域内的技术人员一旦得知了基本创造性概念,则可对这些实施例作出另外的变更和修改。所以,所附权利要求意欲解释为包括优选实施例以及落入本发明范围的所有变更和修改。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。

Claims (4)

1.一种用于航空飞行器的改进的几何参数翼型设计方法,其特征在于,包括如下步骤:
1)、获取与翼型轮廓相关的8个几何参数;
2)、建立翼型弯度表达式:
Figure FDA0002538875360000011
其中,xC和yC分别为翼型弯度上任一点的横坐标和纵坐标,c1、c2、c3、c4分别是三次贝塞尔曲线两个控制点的横、纵坐标,k为控制参数,取值范围为[0,1];
3)、建立翼型厚度表达式:t=t1x0.5+t2x+t3x2+t4x3+t5x4,其中,t为厚度,x为翼型厚度横坐标,t1、t2、t3、t4、t5为控制参数;
4)、构建翼型表达式:该翼型表达式由所述弯度表达式和厚度表达式叠加而成,其采用厚度t和弯度xC,yC的表示式为:
Figure FDA0002538875360000012
其中,xu与yu为翼型上翼面曲线的横纵坐标,xl与yl为翼型下翼面曲线的横纵坐标;
其中,所述的8个几何参数为:相对弯度C,最大弯度所在的弦向位置XC,中弧线后缘夹角αTE,最大弯度处中弧线曲率bXC,相对厚度T,最大厚度所在的弦向位置XT,前缘半径ρ0,后缘夹角βTE
2.根据权利要求1所述的设计方法,其特征在于:
在已知翼型的上述几何参数的情况下,若要推得翼型的具体表达式,则需以上述8个几何参数作为约束,去求解8个优化参数,由于约束与未知数的个数相同,故方程组有唯一解,由于弯度方程为参数方程,需引入最大弯度处的控制参数kc,同时存在控制方程
Figure FDA0002538875360000028
对其取值进行约束,使方程组的解仍唯一;
最终,经过变换,得到弯度反推方程组(4)和厚度反推方程组(5):
Figure FDA0002538875360000022
Figure FDA0002538875360000023
将前缘半径ρ0和后缘夹角βTE分别进行无量纲化处理,得到前缘半径的无量纲量
Figure FDA0002538875360000024
和后缘夹角的无量纲量
Figure FDA0002538875360000025
其中
Figure FDA0002538875360000026
结合上述反推方程组(4)和(5),以及上述的无量纲化处理,求解出表征翼型的8个优化参数c1、c2、c3、c4、XT
Figure FDA0002538875360000027
得到所述翼型表达式中的参数的优化设计空间为:
Figure FDA0002538875360000031
3.根据权利要求2所述的设计方法,其特征在于:所述控制参数t1、t2、t3、t4、t5可以通过以下方程组获得:
Figure FDA0002538875360000032
4.根据权利要求1-3任一项所述的设计方法,其特征在于:对于后缘厚度为0的标准化翼型,满足t(1)=0。
CN201710362161.9A 2017-05-22 2017-05-22 一种改进的几何参数翼型设计方法 Active CN107145677B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710362161.9A CN107145677B (zh) 2017-05-22 2017-05-22 一种改进的几何参数翼型设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710362161.9A CN107145677B (zh) 2017-05-22 2017-05-22 一种改进的几何参数翼型设计方法

Publications (2)

Publication Number Publication Date
CN107145677A CN107145677A (zh) 2017-09-08
CN107145677B true CN107145677B (zh) 2020-09-04

Family

ID=59777289

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710362161.9A Active CN107145677B (zh) 2017-05-22 2017-05-22 一种改进的几何参数翼型设计方法

Country Status (1)

Country Link
CN (1) CN107145677B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108518295A (zh) * 2018-03-28 2018-09-11 中国海洋大学 潮流能水轮机叶片层流翼型设计方法
CN110704944B (zh) * 2019-09-12 2021-10-01 北京航空航天大学 一种面向变弯度翼型的参数化建模方法
CN110737949B (zh) * 2019-11-12 2023-05-30 中国工程物理研究院总体工程研究所 一种筒式发射无人机折叠机翼的发射受力分析方法
CN112001033B (zh) * 2020-09-03 2022-08-26 哈尔滨工程大学 基于联合cst算法的仿生蟹翼型优化设计方法
CN113673032A (zh) * 2021-08-19 2021-11-19 中国人民解放军国防科技大学 一种面向气动外形优化的形状空间实验设计方法
CN113656895A (zh) * 2021-08-20 2021-11-16 中国空气动力研究与发展中心高速空气动力研究所 一种描述超临界翼型的改进Parsec参数化方法
CN118246157B (zh) * 2024-05-29 2024-08-23 中国人民解放军陆军指挥学院 一种旋翼无人机的三维桨叶设计方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7454321B1 (en) * 2002-01-07 2008-11-18 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Robust, optimal subsonic airfoil shapes
CN104392075A (zh) * 2014-12-15 2015-03-04 中国飞机强度研究所 一种翼型参数化建模方法
CN104778327A (zh) * 2015-04-23 2015-07-15 复旦大学 基于人工神经网络的飞机翼型优化设计方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7454321B1 (en) * 2002-01-07 2008-11-18 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Robust, optimal subsonic airfoil shapes
CN104392075A (zh) * 2014-12-15 2015-03-04 中国飞机强度研究所 一种翼型参数化建模方法
CN104778327A (zh) * 2015-04-23 2015-07-15 复旦大学 基于人工神经网络的飞机翼型优化设计方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Perez-Arribas 等.Automatic modelling of airfoil data points.《AEROSPACE SCIENCE AND TECHNOLOGY》.2016,第55卷449-457. *
基于CST参数化方法的飞机翼型快速设计研究;徐亚峰;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑 (月刊 )》;20130415(第04期);C031-14 *
基于实用参数控制的翼型参数化方法;王文剑 等;《飞行设计》;20130630;第33卷(第3期);1-4 *

Also Published As

Publication number Publication date
CN107145677A (zh) 2017-09-08

Similar Documents

Publication Publication Date Title
CN107145677B (zh) 一种改进的几何参数翼型设计方法
US9031817B2 (en) Method for designing a natural laminar flow wing of a supersonic aircraft
CN106050739A (zh) 一种冷却风扇用高性能翼型
CN111125910B (zh) 一种风洞试验气动配平损失评估方法
CN105760635A (zh) 一种基于鸮翼的仿生缝翼的设计方法
CN108038259B (zh) 基于曲率生成气动部件外形的方法
CN112926132B (zh) 一种考虑三维效应影响下的固定翼翼型气动外形设计方法
CN111278728A (zh) 用于飞行器的圆顶整流罩及其制造方法
CN111859545A (zh) 一种考虑升力匹配的宽速域高升阻比机翼优化设计方法
CN112623254B (zh) 一种混合层流机翼吸气能量损耗工程计算方法
CN112632703A (zh) 满足结构约束的机翼翼型前后缘变形外形参数化方法
CN105117541B (zh) 一种正向型架外形优化设计方法
CN103043224B (zh) 生成后缘襟翼舵面翼型前缘曲线的双圆法
CN116702309B (zh) 基于柯西曲线的机翼前缘设计方法及系统
CN109455236B (zh) 一种汽车尾翼结构及其设计方法
CN113353285B (zh) 一种旋翼翼型的确定方法及系统
CN113626935B (zh) 一种高巡航效率跨音速月形机翼设计方法
CN115649417A (zh) 一种高亚声速自配平高隐身翼型
CN109598025B (zh) 一种v形尾翼倾角与面积优化设计方法
CN105260498B (zh) 一种大型民机机翼变弯度设计方法
CN112985753A (zh) 一种用于风洞试验的冰形参数化方法
CN110750836A (zh) 一种基于稳态锥动的飞行器俯仰阻尼力矩计算方法
Zhou Design of F1 Race Car Rear Wing Airfoil: Optimizing the Lift to Drag Ratio through Numerical Simulation
CN113343347B (zh) 一种翼型前缘cst的垂直补偿修正方法及装置
CN115593612B (zh) 一种自配平抗失速高性能翼型

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant