CN106845019A - 一种自适应机翼翼型设计方法 - Google Patents
一种自适应机翼翼型设计方法 Download PDFInfo
- Publication number
- CN106845019A CN106845019A CN201710109049.4A CN201710109049A CN106845019A CN 106845019 A CN106845019 A CN 106845019A CN 201710109049 A CN201710109049 A CN 201710109049A CN 106845019 A CN106845019 A CN 106845019A
- Authority
- CN
- China
- Prior art keywords
- aerofoil profile
- centerdot
- variable
- optimization
- camber
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/12—Computing arrangements based on biological models using genetic models
- G06N3/126—Evolutionary algorithms, e.g. genetic algorithms or genetic programming
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Biophysics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Optimization (AREA)
- Artificial Intelligence (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Aviation & Aerospace Engineering (AREA)
- Automation & Control Theory (AREA)
- Physiology (AREA)
- Genetics & Genomics (AREA)
- Computer Hardware Design (AREA)
- Biomedical Technology (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Structures Of Non-Positive Displacement Pumps (AREA)
Abstract
本发明公开了一种自适应机翼翼型设计方法,首先采用CST参数化方法获得所选择的初始翼型的参数化描述,再确定局部变形的优化设计变量并选择优化算法初始化种群,然后以在升力系数CL=0.9时实现升阻比最大、前缘弯度改变后的最大弯曲应力最小为优化目标,以无塑形变形为约束对种群进行占优排序,最后选择出满足终止条件的优化翼型。本发明以翼型期望的气动特性作为设计目标,建立起自适应机翼翼型同原始翼型的约束关系,且能为变形驱动技术设计提供合理的目标翼型,提高了工程可实现性,进而提高了优化设计效率。
Description
技术领域
本发明属于机翼翼型设计技术领域,具体涉及一种自适应机翼翼型设计方法。
背景技术
机翼对飞机的飞行性能影响极大,与机体的结构和飞机的总体布置也有关系,在全面考虑机翼的参数选择问题时,重点是其剖面形状即翼型和其平面形状几何参数的选择。机翼翼型以及其在机翼上的配置情况,对气动特性影响极大,显然,只有选用良好的翼型并进行正确的配置,才可能保证机翼具有良好的气动特性。
在目前的机翼翼型设计过程中,翼型形状通常是根据某一特定飞行状态、飞行环境优化得到的,在满足该状态飞行器性能最优的前提下,通过舵面、增升装置等辅助机构,实现该翼型在其他典型状态下的使用需求,这就意味着飞行器在非设计点使用时性能有所损失。由于上述传统翼型设计方法的局限性以及人们对飞行器性能、经济性等的要求也是越来越提高,在上世纪八十年代,美国学者提出任务自适应概念,成为解决上述矛盾的一个重要研究方向。所谓自适应机翼翼型就是通过改变原始翼型几何特征,如弯度,以此来满足机翼在不同工况下设计要求。
从目前国内关于自适应机翼翼型的公开专利来看,部分工作专注于开发变形驱动技术,在其设计过程中仅能以原始机翼翼型及其采用的驱动技术为约束,而没有合理的能够满足气动特性要求的自适应机翼翼型作为校核指标,因而设计较为盲目;也有部分文献开展自适应机翼翼型气动特性设计,但其方案仍采用传统翼型优化手段,仅考虑翼型改变后的气动特性指标,设计过程未能提出合理的约束,使得设计出来的翼型在工程上并不实用。如果简单的将上述两方面工作结合以弥补各自不足,那么整个设计流程效率较低并不适用。
发明内容
针对现有技术的问题,本发明提供了一种自适应机翼翼型设计方法,该方法以翼型期望的气动特性作为设计目标,建立起自适应机翼翼型同原始翼型的约束关系,且能为变形驱动技术设计提供合理的目标翼型,提高了工程可实现性。
为了达到上述的目的,本发明所采用的技术方案是:
一种自适应机翼翼型设计方法,该方法包括以下步骤:
步骤1:选择初始翼型,采用CST参数化方法获得初始翼型的参数化描述;
步骤2:确定局部变形的优化设计变量,将影响弯度部分的主要权重系数作为优化变量;
步骤3:选择优化算法初始化种群;
步骤4:将初始种群进行处理,以在升力系数CL=0.9时实现升阻比最大、前缘弯度改变后的最大弯曲应力最小时为优化目标,以前缘变弯后不产生轴向应力为约束对种群进行占优排序;
步骤5:判断优化目标是否满足终止条件:
若优化目标满足终止条件,则结束优化过程,并输出优化的新翼型结果;
若优化目标不满足终止条件,则对初始种群进行修改,并回到步骤3。
进一步的是,采用CST参数化方法对初始翼型进行参数化描述,该方法具有良好的形状控制能力,同时隐含连续性约束,能够保证在翼型上下型线变形后依然能够实现二阶导数连续,该翼型型线数据点参数化定义如下:
其中y表示翼型数据点纵坐标,x翼型数据点横坐标,c为翼型弦长,zle为前缘纵坐标,zte尾缘纵坐标,的公式如下:
在描述翼型时,N1和N2为常量,的具体公式如下:
其中bi为控制点权重,通过zle可以调整前缘位置,通过bi可以调整曲线形状。
进一步的是,确定局部变形的优化设计变量并进行初始化种群的方法如下:
由步骤1可知,该方法的可变参数为(zle,zte,b0low,...,bnlow,b0up,...,bnup),下标up和low分别代表翼型上下型线,在众多可变参数中,zte可由翼型尾缘数据确定;b0low,b0up可以合并用前缘半径Rle代替,b0low,b0up合并后可以保证前缘曲率半径连续。在其余控制点权重参数中,可以根据可变弯度部分在翼型中的占比(A)确定下述权重为待优化变量:
其余非优化变量权重系数(b(k+1)low,...,bnlow,b(k+1)up,...,bnup)可根据原始翼型数据型线数据求得,由于上下型线权重系数求法相同,下面以上型线权重系数求解方法进行说明,非优化变量权重系数应满足下式:
其中下标m代表有m个已知点,x1,…,xm代表翼固定段数据点,将其整理可知其形式为Xm×(n-k)b(n-k)×1=Ym×1,其中Xm×(n-k)为其系数矩阵,b(n-k)×1为其待求权重系数矩阵,由于待拟合点x1,…,xm数量远多于(n-k),即上式未知量个数远大于方程组个数,因此该线性方程组没有传统意义下的解,通过求解其法方程组的方法求得其余权重系数,即b=(XTX)-1XTY,用同样方法可以求得下型线非优化变量权重系数;因此确定优化变量有(zle,Rle,b1low,...,bklow,b1up,...,bkup),根据选用的优化算法初始化种群,优化算法主要有遗传算法、蚁群算法、禁忌搜索和粒子群算法等。
进一步的是,将种群内每个样本通过CST参数化方法生成翼型数据文件以Pointwise支持的Segment格式文件输出,通过PowerShell及Pointwise软件支持的Glyph语言调用Pointwise并实现动网格方法建立变形后翼型的计算网格,此后通过PowerShell及Ansys Fluent的jou脚本语言对CFD求解过程中的控制方程、离散格式、求解方法、边界条件、流场初始化等参数进行设置,最终迭代求解得到弯度改变后的翼型在不同攻角下的升阻比,通过拟合得到在CL=0.9时升阻比,CL=0.9时的升阻比最大设计指标是本方法的一个设计工况点。
进一步的是,由于在翼型设计一般先于结构设计,在无准确结构的情况下,无法通过有限元精确分析受力,因此采用纯弯无剪力的梁模型来估算弯度改变后的最大弯曲应力,在气动设计时,假设翼型型线为等厚度均质截面梁,则变弯度后弯曲梁应力计算公式如下:
其中E为材料杨氏模量,z为梁中性面到表面距离,为翼型型线弯度,可知最大应力发生在弯度变化最大处;由于CST参数化方法建立翼型型线描述二阶处处连续,因此可以通过变弯度后对应弯度变化量,进而计算得到最大弯曲应力σmax。
进一步的是,自适应翼型在弯度变化后,翼型型线长度不变成为一个重要约束,该约束满足应通过以下过程保证:
通过几何关系计算自适应机翼翼型变弯度后的上下型线长度,应等于原始翼型型线长度,即满足下式:
其中Loriginal为初始翼型变形区前缘弧长,yup为翼型上型线斜率,ylow为下型线斜率。
进一步的是,采用基于Pareto多目标遗传算法进行占优排序,在多目标优化过程中,种群个体需要根据优化结果及约束被赋予Pareto秩,以区别每一个种群优劣,进而控制该个体对下一代种群产生影响的能力大小,在本例中Pareto秩根据多目标遗传优化算法中占优排序的方式定义,占优的种群个体Pareto秩越小,代表符合预期性能越优,在有约束时,将种群中i个体对j个体占优定义如下:
(1)单约束时,i、j不满足约束,i接近约束边界;
(2)多约束时,i、j不满足约束,i约束对j约束占优;
(3)i满足约束、j不满足约束;
(4)i与j都满足约束,且i对j占优。
在本方法中,是单约束多目标的优化过程,在升力系数CL=0.9时升阻比大为优,最大弯曲应力最小为优,约束按下式进行计算:
当得到满足约束及优化目标的Pareto前沿时,在从Pareto前沿上综合选择性能最优结果。
进一步的是,若优化目标不满足终止条件,通过多目标遗传优化算法交叉、变异等方法对初始种群进行修改,并回到步骤3。
进一步的是,所述终止条件指优化迭代次数达到预设次数,该预设次数根据优化时间限制进行选择,优化过程中,优化迭代次数越多,一般优化结果越佳,但优化时间相应也会增加。
本发明的有益效果是:
在本发明中提出了一种自适应机翼翼型设计方法,该方法通过在翼型优化设计过程中考虑升阻比、最大弯曲应力等几何结构特征参数,并且建立起自适应机翼翼型同原始翼型的约束关系即无塑形变形为约束,同时又以期望的气动特性作为设计目标,使得通过该设计方法得到的翼型能够为变形驱动技术设计提供合理的较高性能的目标翼型,也能提高自适应机翼翼型的工程可实现性,进而提高优化设计效率。
附图说明
图1是本发明的设计流程图;
图2是原始翼型与自适应机翼翼型几何形状图;
图3是原始翼型与自适应机翼翼型的极曲线图。
具体实施方式
为了使发明的目的、技术方案和优点更加清楚,下面结合附图对本发明作进一步阐述。
在本实施例中,如图1所示,一种自适应机翼翼型设计方法,该方法包括以下步骤:
步骤1:选择RAE2822翼型作为初始翼型,采用CST参数化方法获得初始翼型的参数化描述;
步骤2:确定局部变形的优化设计变量,将影响弯度部分的主要权重系数作为优化变量;
步骤3:选择优化算法初始化种群;
步骤4:将初始种群进行处理,以在升力系数CL=0.9时实现升阻比最大、前缘弯度改变后的最大弯曲应力最小时为优化目标,以前缘变弯后不产生轴向应力为约束对种群进行占优排序;
步骤5:判断优化目标是否满足终止条件:
若优化目标满足终止条件,则结束优化过程,并输出优化的新翼型结果;
若优化目标不满足终止条件,则对初始种群进行修改,并回到步骤3。
作为一种优化方案,采用CST参数化方法对初始翼型进行参数化描述,该方法具有良好的形状控制能力,同时隐含连续性约束,能够保证在翼型上下型线变形后依然能够实现二阶导数连续,该翼型型线数据点参数化定义如下:
其中y表示翼型数据点纵坐标,x翼型数据点横坐标,c为翼型弦长,zle为前缘纵坐标,zte尾缘纵坐标,的公式如下:
在描述翼型时,N1和N2为常量,其中,N1=0.5,N2=1,的具体公式如下:
其中bi为控制点权重,n为伯恩斯坦多项式阶数,通过zle可以调整前缘位置,通过bi可以调整曲线形状。
作为一种优化方案,确定局部变形的优化设计变量并进行初始化种群的方法如下:
由步骤1可知,该方法的可变参数为(zle,zte,b0low,...,bnlow,b0up,...,bnup),下标up和low分别代表翼型上下型线,在众多可变参数中,zte可由翼型尾缘数据确定;b0low,b0up可以合并用前缘半径Rle代替,b0low,b0up合并后可以保证前缘曲率半径连续。在其余控制点权重参数中,可以根据可变弯度部分在翼型中的占比(A)确定下述权重为待优化变量:
在本实施例中,A取20%,
该式意义是将伯恩斯坦函数中影响变弯度部分的主要权重系数作为优化变量,其值由优化算法随机赋值初始化,其余非优化变量权重系数(b(k+1)low,...,bnlow,b(k+1)up,...,bnup)可根据原始翼型数据型线数据求得,由于上下型线权重系数求法相同,下面以上型线权重系数求解方法进行说明,非优化变量权重系数应满足下式:
其中下标m代表有m个已知点,x1,…,xm代表翼固定段数据点,将其整理可知其形式为Xm×(n-k)b(n-k)×1=Ym×1,其中Xm×(n-k)为其系数矩阵,b(n-k)×1为其待求权重系数矩阵,由于待拟合点x1,…,xm数量远多于(n-k),即上式未知量个数远大于方程组个数,因此该线性方程组没有传统意义下的解,通过求解其法方程组的方法求得其余权重系数,即b=(XTX)-1XTY,用同样方法可以求得下型线非优化变量权重系数;因此确定优化变量有(zle,Rle,b1low,...,bklow,b1up,...,bkup),根据选用的优化算法初始化种群,优化算法主要有遗传算法、蚁群算法、禁忌搜索和粒子群算法等,在本实施例中,种群的生成采用十进制编码,确定种群样本数量N,根据待优化变量取值范围及待优化变量数量,随机生成N串十进制数字序列。
作为一种优化方案,将种群内每个样本通过CST参数化方法生成翼型数据文件以Pointwise支持的Segment格式文件输出,通过PowerShell及Pointwise软件支持的Glyph语言调用Pointwise并实现动网格方法建立变形后翼型的计算网格,此后通过PowerShell及AnsysFluent的jou脚本语言对CFD求解过程中的控制方程、离散格式、求解方法、边界条件、流场初始化等参数进行设置,最终迭代求解得到弯度改变后的翼型在不同攻角下的升阻比,通过拟合得到在CL=0.9时升阻比,CL=0.9时的升阻比最大设计指标是本方法的一个设计工况点。
作为一种优化方案,由于在翼型设计一般先于结构设计,在无准确结构的情况下,无法通过有限元精确分析受力,因此采用纯弯无剪力的梁模型来估算弯度改变后的最大弯曲应力,在气动设计时,假设翼型型线为等厚度均质截面梁,则变弯度后弯曲梁应力计算公式如下:
其中E为材料杨氏模量,z为梁中性面到表面距离,为翼型型线弯度,可知最大应力发生在弯度变化最大处;由于CST参数化方法建立翼型型线描述二阶处处连续,因此可以通过变弯度后对应弯度变化量,进而计算得到最大弯曲应力σmax。
作为一种优化方案,自适应翼型在弯度变化后,翼型型线长度不变成为一个重要约束,该约束满足应通过以下过程保证:
通过几何关系计算自适应机翼翼型变弯度后的上下型线长度,应等于原始翼型型线长度,即满足下式:
其中Loriginal为初始翼型变形区前缘弧长,yup为翼型上型线斜率,ylow为下型线斜率。
作为一种优化方案,采用基于Pareto多目标遗传算法进行占优排序,在多目标优化过程中,种群个体需要根据优化结果及约束被赋予Pareto秩,以区别每一个种群优劣,进而控制该个体对下一代种群产生影响的能力大小,在本例中Pareto秩根据多目标遗传优化算法中占优排序的方式定义,占优的种群个体Pareto秩越小,代表符合预期性能越优,在有约束时,将种群中i个体对j个体占优定义如下:
(1)单约束时,i、j不满足约束,i接近约束边界;
(2)多约束时,i、j不满足约束,i约束对j约束占优;
(3)i满足约束、j不满足约束;
(4)i与j都满足约束,且i对j占优。
在本方法中,是单约束多目标的优化过程,在升力系数CL=0.9时升阻比大为优,最大弯曲应力最小为优,约束按下式进行计算:
作为一种优化方案,若优化目标不满足终止条件,通过多目标遗传优化算法交叉、变异等方法对初始种群进行修改,并回到步骤3。
作为一种优化方案,所述终止条件指优化迭代次数达到预设次数,该预设次数根据优化时间限制进行选择,优化过程中,优化迭代次数越多,一般优化结果越佳,但优化时间相应也会增加,本例中,优选优化迭代次数为20次。
当得到满足约束及优化目标的Pareto前沿时,在从Pareto前沿上综合选择性能最优的结果。
本发明中的方法通过在翼型优化设计过程中考虑升阻比、最大弯曲应力等几何结构特征参数,建立起自适应机翼翼型同原始翼型的约束关系,同时又以期望的气动特性作为设计目标,使得通过该设计方法得到的翼型能够为变形驱动技术设计提供合理目标翼型,也能提高自适应机翼翼型的工程可实现性,进而提高优化设计效率,相对传统的机翼翼型设计方法更为全面,优化后的机翼翼型与优化前机翼翼型几何形状对比图如图2所示,优化后的机翼翼型与优化前机翼翼型极曲线对比图如图3所示。
以上显示和描述了本发明的基本原理、主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
Claims (9)
1.一种自适应机翼翼型设计方法,其特征在于:该方法包括以下步骤:
步骤1:选择初始翼型,采用CST参数化方法获得初始翼型的参数化描述;
步骤2:确定局部变形的优化设计变量,将影响弯度部分的主要权重系数作为优化变量;
步骤3:选择优化算法初始化种群;
步骤4:将初始种群进行处理,以在升力系数CL=0.9时实现升阻比最大、前缘弯度改变后的最大弯曲应力最小时为优化目标,以前缘变弯后不产生轴向应力为约束对种群进行占优排序;
步骤5:判断优化目标是否满足终止条件:
若优化目标满足终止条件,则结束优化过程,并输出优化的新翼型结果;
若优化目标不满足终止条件,则对初始种群进行修改,并回到步骤3。
2.根据权利要求1所述的一种自适应机翼翼型设计方法,其特征在于:采用CST参数化方法对初始翼型进行参数化描述,该翼型型线数据点参数化定义如下:
其中y表示翼型数据点纵坐标,x翼型数据点横坐标,c为翼型弦长,zle为前缘纵坐标,zte尾缘纵坐标,的公式如下:
在描述翼型时,N1和N2为常量,的具体公式如下:
其中bi为控制点权重,通过zle可以调整前缘位置,通过bi可以调整曲线形状。
3.根据权利要求2所述的一种自适应机翼翼型设计方法,其特征在于:确定局部变形的优化设计变量并进行初始化种群的方法如下:
由步骤1可知,该方法的可变参数为(zle,zte,b0low,...,bnlow,b0up,...,bnup),下标up和low分别代表翼型上下型线,在众多可变参数中,zte可由翼型尾缘数据确定;b0low,b0up可以合并用前缘半径Rle代替,在其余控制点权重参数中,可以根据可变弯度部分在翼型中的占比(A)确定下述权重为待优化变量:
其余非优化变量权重系数(b(k+1)low,...,bnlow,b(k+1)up,...,bnup)可根据原始翼型数据型线数据求得,由于上下型线权重系数求法相同,下面以上型线权重系数求解方法进行说明,非优化变量权重系数应满足下式:
其中下标m代表有m个已知点,x1,...,xm代表翼固定段数据点,将其整理可知其形式为Xm×(n-k)b(n-k)×1=Ym×1,其中Xm×(n-k)为其系数矩阵,b(n-k)×1为其待求权重系数矩阵,由于待拟合点x1,...,xm数量远多于(n-k),即上式未知量个数远大于方程组个数,因此该线性方程组没有传统意义下的解,通过求解其法方程组的方法求得其余权重系数,即b=(XTX)-1XTY,用同样方法可以求得下型线非优化变量权重系数;因此确定优化变量有(zle,Rle,b1low,...,bklow,b1up,...,bkup),根据选用的优化算法初始化种群。
4.根据权利要求3所述的一种自适应机翼翼型设计方法,其特征在于:
将种群内每个样本通过CST参数化方法生成翼型数据文件以Pointwise支持的Segment格式文件输出,通过PowerShell及Pointwise软件支持的Glyph语言调用Pointwise并实现动网格方法建立变形后翼型的计算网格,此后通过PowerShell及Ansys Fluent的jou脚本语言对CFD求解过程中的控制方程、离散格式、求解方法、边界条件、流场初始化等参数进行设置,最终迭代求解得到弯度改变后的翼型在不同攻角下的升阻比,通过拟合得到在CL=0.9时升阻比。
5.根据权利要求4所述的一种自适应机翼翼型设计方法,其特征在于:
在气动设计时,假设翼型型线为等厚度均质截面梁,则变弯度后弯曲梁应力计算公式如下:
其中E为材料杨氏模量,z为梁中性面到表面距离,为翼型型线弯度,可知最大应力发生在弯度变化最大处;由于CST参数化方法建立翼型型线描述二阶处处连续,因此可以通过变弯度后对应弯度变化量,进而计算得到最大弯曲应力σmax。
6.根据权利要求5所述的一种自适应机翼翼型设计方法,其特征在于:
自适应翼型在弯度变化后,翼型型线长度不变成为一个重要约束,该约束满足应通过以下过程保证:
通过几何关系计算自适应机翼翼型变弯度后的上下型线长度,应等于原始翼型型线长度,即满足下式:
其中Loriginal为初始翼型变形区前缘弧长,yup为翼型上型线斜率,ylow为下型线斜率。
7.根据权利要求6所述的一种自适应机翼翼型设计方法,其特征在于:采用基于Pareto多目标遗传算法进行占优排序。
8.根据权利要求7所述的一种自适应机翼翼型设计方法,其特征在于:若优化目标不满足终止条件,通过多目标遗传优化算法交叉、变异等方法对初始种群进行修改,并回到步骤3。
9.根据权利要求8所述的一种自适应机翼翼型设计方法,其特征在于:所述终止条件指优化迭代次数达到预设次数,该预设次数根据优化时间限制进行选择。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710109049.4A CN106845019B (zh) | 2017-02-27 | 2017-02-27 | 一种自适应机翼翼型设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710109049.4A CN106845019B (zh) | 2017-02-27 | 2017-02-27 | 一种自适应机翼翼型设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106845019A true CN106845019A (zh) | 2017-06-13 |
CN106845019B CN106845019B (zh) | 2020-04-10 |
Family
ID=59134155
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710109049.4A Active CN106845019B (zh) | 2017-02-27 | 2017-02-27 | 一种自适应机翼翼型设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106845019B (zh) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109670271A (zh) * | 2019-01-15 | 2019-04-23 | 上海交通大学 | 基于敏感点多点增量成形的大型薄壁加筋板的校形方法 |
CN109783858A (zh) * | 2018-12-13 | 2019-05-21 | 航天神舟飞行器有限公司 | 一种用于低雷诺数翼型优化的方法 |
CN110110384A (zh) * | 2019-04-12 | 2019-08-09 | 电子科技大学 | 一种基于浓度策略的非支配遗传算法在飞机翼形设计优化中的应用方法 |
CN110704944A (zh) * | 2019-09-12 | 2020-01-17 | 北京航空航天大学 | 一种面向变弯度翼型的参数化建模方法 |
CN110750839A (zh) * | 2019-08-26 | 2020-02-04 | 深圳市道通智能航空技术有限公司 | 一种螺旋桨翼型的设计方法及终端设备 |
CN110816874A (zh) * | 2019-10-11 | 2020-02-21 | 成都飞机工业(集团)有限责任公司 | 一种通过地面滑行试验辨识双发飞机平衡极曲线的方法 |
CN111967090A (zh) * | 2020-08-04 | 2020-11-20 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种优化设计空间的动态改进方法 |
CN112001033A (zh) * | 2020-09-03 | 2020-11-27 | 哈尔滨工程大学 | 基于联合cst算法的仿生蟹翼型优化设计方法 |
CN112580151A (zh) * | 2020-12-27 | 2021-03-30 | 复旦大学 | 基于目标空间的自然层流翼型的优化问题指标获取方法 |
CN112632703A (zh) * | 2020-12-24 | 2021-04-09 | 中国航空工业集团公司沈阳空气动力研究所 | 满足结构约束的机翼翼型前后缘变形外形参数化方法 |
CN113195994A (zh) * | 2019-03-20 | 2021-07-30 | 日挥环球株式会社 | 热交换器群的设计方法及处理工厂 |
CN113343347A (zh) * | 2021-02-01 | 2021-09-03 | 复旦大学 | 一种翼型前缘cst的垂直补偿修正方法及装置 |
CN113361017A (zh) * | 2021-07-03 | 2021-09-07 | 西北工业大学 | 一种考虑机构约束的管射无人机二次折叠翼气动设计方法 |
CN116796387A (zh) * | 2023-06-16 | 2023-09-22 | 安徽工业大学 | 基于全信息自适应合作博弈理论的风力机翼型优化方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102060101A (zh) * | 2010-12-21 | 2011-05-18 | 南京航空航天大学 | 一种用于变体机翼的蒙皮 |
KR20130139089A (ko) * | 2012-06-12 | 2013-12-20 | 금오공과대학교 산학협력단 | 잠수함 최적화 설계 탐색시스템 |
-
2017
- 2017-02-27 CN CN201710109049.4A patent/CN106845019B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102060101A (zh) * | 2010-12-21 | 2011-05-18 | 南京航空航天大学 | 一种用于变体机翼的蒙皮 |
KR20130139089A (ko) * | 2012-06-12 | 2013-12-20 | 금오공과대학교 산학협력단 | 잠수함 최적화 설계 탐색시스템 |
Non-Patent Citations (3)
Title |
---|
CHUANZHENLIU 等: "Application of the 3D multi-block CST method to hypersonic aircraft optimization", 《AEROSPACE SCIENCE AND TECHNOLOGY》 * |
付鸿雁 等: "简化自适应机翼的气动外形优化设计", 《空气动力学学报》 * |
刘航 等: "自适应翼型的气动外形优化设计", 《航空学报》 * |
Cited By (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109783858A (zh) * | 2018-12-13 | 2019-05-21 | 航天神舟飞行器有限公司 | 一种用于低雷诺数翼型优化的方法 |
CN109783858B (zh) * | 2018-12-13 | 2023-09-29 | 航天神舟飞行器有限公司 | 一种用于低雷诺数翼型优化的方法 |
CN109670271B (zh) * | 2019-01-15 | 2023-03-07 | 上海交通大学 | 基于敏感点多点增量成形的大型薄壁加筋板的校形方法 |
CN109670271A (zh) * | 2019-01-15 | 2019-04-23 | 上海交通大学 | 基于敏感点多点增量成形的大型薄壁加筋板的校形方法 |
CN113195994B (zh) * | 2019-03-20 | 2022-05-31 | 日挥环球株式会社 | 热交换器群的设计方法及处理工厂 |
CN113195994A (zh) * | 2019-03-20 | 2021-07-30 | 日挥环球株式会社 | 热交换器群的设计方法及处理工厂 |
US11879687B2 (en) | 2019-03-20 | 2024-01-23 | Jgc Corporation | Method of designing heat exchanger group and processing plant |
CN110110384A (zh) * | 2019-04-12 | 2019-08-09 | 电子科技大学 | 一种基于浓度策略的非支配遗传算法在飞机翼形设计优化中的应用方法 |
CN110750839A (zh) * | 2019-08-26 | 2020-02-04 | 深圳市道通智能航空技术有限公司 | 一种螺旋桨翼型的设计方法及终端设备 |
CN110704944A (zh) * | 2019-09-12 | 2020-01-17 | 北京航空航天大学 | 一种面向变弯度翼型的参数化建模方法 |
CN110704944B (zh) * | 2019-09-12 | 2021-10-01 | 北京航空航天大学 | 一种面向变弯度翼型的参数化建模方法 |
CN110816874A (zh) * | 2019-10-11 | 2020-02-21 | 成都飞机工业(集团)有限责任公司 | 一种通过地面滑行试验辨识双发飞机平衡极曲线的方法 |
CN110816874B (zh) * | 2019-10-11 | 2022-04-08 | 成都飞机工业(集团)有限责任公司 | 一种通过地面滑行试验辨识双发飞机平衡极曲线的方法 |
CN111967090A (zh) * | 2020-08-04 | 2020-11-20 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种优化设计空间的动态改进方法 |
CN112001033B (zh) * | 2020-09-03 | 2022-08-26 | 哈尔滨工程大学 | 基于联合cst算法的仿生蟹翼型优化设计方法 |
CN112001033A (zh) * | 2020-09-03 | 2020-11-27 | 哈尔滨工程大学 | 基于联合cst算法的仿生蟹翼型优化设计方法 |
CN112632703A (zh) * | 2020-12-24 | 2021-04-09 | 中国航空工业集团公司沈阳空气动力研究所 | 满足结构约束的机翼翼型前后缘变形外形参数化方法 |
CN112632703B (zh) * | 2020-12-24 | 2024-05-07 | 中国航空工业集团公司沈阳空气动力研究所 | 满足结构约束的机翼翼型前后缘变形外形参数化方法 |
CN112580151B (zh) * | 2020-12-27 | 2022-03-18 | 复旦大学 | 基于目标空间的自然层流翼型的优化问题指标获取方法 |
CN112580151A (zh) * | 2020-12-27 | 2021-03-30 | 复旦大学 | 基于目标空间的自然层流翼型的优化问题指标获取方法 |
CN113343347B (zh) * | 2021-02-01 | 2022-06-17 | 复旦大学 | 一种翼型前缘cst的垂直补偿修正方法及装置 |
CN113343347A (zh) * | 2021-02-01 | 2021-09-03 | 复旦大学 | 一种翼型前缘cst的垂直补偿修正方法及装置 |
CN113361017A (zh) * | 2021-07-03 | 2021-09-07 | 西北工业大学 | 一种考虑机构约束的管射无人机二次折叠翼气动设计方法 |
CN113361017B (zh) * | 2021-07-03 | 2022-09-16 | 西北工业大学 | 一种考虑机构约束的管射无人机二次折叠翼气动设计方法 |
CN116796387A (zh) * | 2023-06-16 | 2023-09-22 | 安徽工业大学 | 基于全信息自适应合作博弈理论的风力机翼型优化方法 |
CN116796387B (zh) * | 2023-06-16 | 2024-03-26 | 安徽工业大学 | 基于全信息自适应合作博弈理论的风力机翼型优化方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106845019B (zh) | 2020-04-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106845019A (zh) | 一种自适应机翼翼型设计方法 | |
EP3926507A1 (en) | Helicopter rotor airfoil determination method and system | |
CN104408260B (zh) | 一种潮流能水轮机叶片翼型设计方法 | |
CN109484623B (zh) | 宽速域大升力线斜率对称翼型设计方法及翼型 | |
CN109190283A (zh) | 一种考虑高湍流自由来流效应的风力机翼型气动优化方法 | |
CN106126860B (zh) | 一种考虑加工误差的高超声速机翼鲁棒优化设计方法 | |
Li et al. | Aerodynamic optimization of wind turbine airfoils using response surface techniques | |
CN109460566A (zh) | 一种风力机叶片内侧厚翼型的气动稳健优化设计方法 | |
CN101923584B (zh) | 风力机专用翼型设计方法及风力机专用翼型 | |
CN106055791B (zh) | 基于预估校正算法的飞行器全局气动优化方法 | |
CN101615216B (zh) | 对翼型尾缘进行流线形增厚的方法 | |
CN104834772A (zh) | 基于人工神经网络的飞机翼型/机翼反设计方法 | |
CN107180134A (zh) | 一种可重复使用的天地往返飞行器外形设计方法 | |
CN102235325A (zh) | 基于翼型集成和混合尾缘改型的风力机叶尖翼型设计方法 | |
CN108733914A (zh) | 基于人工神经网络的跨音速翼型自然层流延迟转捩设计方法 | |
Rogalsky et al. | Hybridization of differential evolution for aerodynamic design | |
CN104863799A (zh) | 一种利用贝塞尔函数曲线的风力机翼型设计方法 | |
CN107194122A (zh) | 一种变速风机叶片的多工况气动优化的改进方法 | |
CN112926132B (zh) | 一种考虑三维效应影响下的固定翼翼型气动外形设计方法 | |
Brown et al. | Rapidly recovering wind turbine wakes with dynamic pitch and rotor speed control | |
CN113742861A (zh) | 一种适用于风力发电机风洞试验的叶片模型优化设计方法 | |
CN111310282A (zh) | 一种适用于高原环境的直升机旋翼翼型生成方法及系统 | |
CN106741924A (zh) | 一种全工况条件下的高升阻比低力矩特性9%厚度旋翼翼型 | |
CN105787217B (zh) | 一种飞机用波纹翼型的优化设计方法 | |
CN106227985A (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 |