CN102789539B - 一种飞行器升力面结构优化设计方法 - Google Patents
一种飞行器升力面结构优化设计方法 Download PDFInfo
- Publication number
- CN102789539B CN102789539B CN201210158035.9A CN201210158035A CN102789539B CN 102789539 B CN102789539 B CN 102789539B CN 201210158035 A CN201210158035 A CN 201210158035A CN 102789539 B CN102789539 B CN 102789539B
- Authority
- CN
- China
- Prior art keywords
- space
- design
- lift face
- high accuracy
- aircraft lift
- 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
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种飞行器升力面结构优化设计方法:①采用LHD在设计空间内采取试验设计样本并获得飞行器升力面高精度分析模型响应值;②构造径向基函数代理模型;③在设计空间使用LHD方法随机生成一定数量的网格点,计算其对应的径向基函数代理模型响应值;④给定全局空间缩减率M1%和局部空间缩减率M2%,进行2次伪缩减;⑤在每个聚类空间内利用遗传算法对径向基函数代理模型寻优;⑥得到最优点与上次空间全局近似最优点的空间距离,若满足缩减准则,则删除当前空间得到新空间;⑦获取当前全局近似最优解的真实响应值,判断其是否收敛,如收敛则停止优化,否则反复迭代直至找到最优解。本发明提出的方法提高了优化效率,节约了飞行器升力面结构优化设计成本。
Description
技术领域
本发明涉及一种飞行器升力面结构优化设计方法,属于飞行器优化设计中多学科耦合优化设计技术领域。
背景技术
上世纪80年代起,美国NASA Langley研究中心的Sobieski等一批航空领域的科学家和工程设计人员,提出一种新的飞行器设计方法:多学科设计优化(Multidisciplinary Design Optimization,MDO)。MDO的提出受到了各国研究机构、高校和工业界的普遍重视(如NASA、斯坦福大学、麻省理工学院、Airbus公司和Boeing公司等),掀起了一股MDO的研究热潮。一方面,为了提高分析精度和可信度,高精度分析模型被广泛用于飞行器设计中,例如结构分析中采用的有限元分析(FEA)模型、气动分析中使用的计算流体力学(CFD)分析模型等。高精度分析模型在提高分析精度和可信度的同时也带来了计算耗时的困难,虽然当今计算机软硬件技术已经有了长足的发展,然而,调用高精度分析模型完成一次分析仍然极其耗时,例如使用CFD模型完成一次气动仿真分析需要数小时甚至数十小时。另一方面,飞行器设计问题往往涉及多个相互耦合的学科。譬如,飞行器设计涉及气动、结构、动力、隐身、控制等学科,各学科相互影响,相互制约,飞行器的性能是各学科耦合的综合体现。由于学科之间的耦合关系,飞行器设计问题的系统分析表现为多学科分析。本质上,多学科分析过程是一个典型的非线性求解过程,每次多学科分析都需要进行多次迭代,计算耗时,如果各学科都采用高精度分析模型,其计算量将非常庞大。再次,飞行器优化设计过程中需要经过反复迭代方能收敛到局部或全局最优解,而每次迭代都需要进行多次飞行器设计问题的多学科分析,可见计算成本将进一步增加。因此,直接采用传统的全局优化方法与高精度飞行器分析模型来解决飞行器优化设计会进一步增加设计成本。如何解决飞行器多学科优化设计这个难题困扰许多设计专家。为了解决这样的问题,基于代理模型的飞行器优化设计方法引起了越来越多的注意。
基于代理模型的飞行器升力面多学科耦合设计优化方法本质就是构造分析结果与高精度模型接近(包括气动、结构、隐身等学科分析模型以及全系统分析模型),计算成本更低的代理模型并代替高精度分析模型用于优化。由于高精度分析模型计算一次所需时间的量级为小时或天,而代理模型计算一次所用时间的量级仅为秒甚至毫秒,因此与高精度分析模型的计算时间相比,代理模型以及基于代理模型优化的计算时间往往可以忽略不计。目前常用的代理模型方法包括多项式响应面、Kriging模型、径向基函数、移动最小二乘法以及神经网络等。
多项式响应面能够快速有效地解决低阶相对简单分析模型;Kriging模型能够对一个未知的高阶非线性分析模型进行高精度构造,但是该模型很难获得与使用。径向基函数(Radial Basis Function,RBF)是最常用的代理模型方法之一,其优点在于对于高阶非线性的飞行器高精度分析模型,径向基函数在全局近似精度较高,并且随着试验设计样本的增加,所构造的径向基代理模型的近似精度会提高;在试验设计样本附近,近似精度较高。按照代理模型在优化过程中的使用方式可以分为静态代理模型和动态代理模型。静态代理模型在优化前构造完毕,在优化过程中代理模型保持不变,而动态代理模型在优化过程中根据已知信息逐步进行更新直至优化收敛。与静态代理模型相比,动态代理模型在优化效率和结果精度方面更具有优势。
设计空间是飞行器升力面气动、热、结构多学科耦合设计优化中较重要的一环。一方面,设计空间过小可能使优化设计快速地收敛到局部最优解;另一方面,设计空间过大可能使优化很难收敛,甚至不能收敛。因此,从较大的设计空间放缩到合适的包括全局最优的设计空间是飞行器升力面气动、热、结构多学科耦合设计优化中重要的一步。
径向基函数(Radial Basis Function,RBF)是最常用的代理模型方法之一,其优点在于对于高阶非线性的飞行器高精度分析模型,径向基函数在全局近似精度较高;并且随着试验设计样本的增加,所构造的径向基代理模型的近似精度会提高;在试验设计样本附近,近似精度较高。按照代理模型在优化过程中的使用方式可以分为静态代理模型和自适应代理模型。静态代理模型是通过一次试验设计采取足够多的试验设计样本,然后构造代理模型,在优化过程中代理模型保持不变;而自适应代理模型是在优化设计过程序列采取试验设计样本,然后在每次优化迭代过程中根据已知信息逐步改进和更新代理模型,直至优化收敛。与静态代理模型相比,自适应代理模型在优化效率和结果精度方面更具有优势。
基于代理模型的优化技术在航空航天领域具有很广的应用前景,但是国内研究开展较晚。国外研究机构、工业界和商业软件公司纷纷给予该类算法开发了高效的优化器并用于飞行器优化设计,例如播音公司的Boeing探索器,Altair公司Hyperstudy中的ARSM优化器,以及美国Sandia国家研究实验室的DAKOTA等。
径向基函数(RBF)代理模型的基本形式为:
其中,f(x)是代理模型的响应值;B为权向系数向量,B=[β1,…,βm]T且B应满足公式(2);βi是第i个权向系数,1≤i≤m;m是试验设计样本个数;是基函数向量,φ是基函数,||·||是二范数,x是预测点;xi是第i个试验设计样本。
fi=yi (2)
其中,fi为第i个试验设计样本的径向基函数代理模型预测值;yi为第i个试验设计样本所对应的分析模型响应值。于是有:
AB=Y (3)
其中,Y为试验设计样本所对应的分析模型响应值向量。
由公式(3)可得公式(4):
B=A-1Y (4)
本发明采用的另外一种重要技术是模糊聚类法。
在缩减设计空间时,模糊聚类法用于对网格点的隶属度进行分类,完成数据的模糊聚类划分,通过聚类的放缩与删减,最终得到全局最优解所在的重点区域。模糊c-均值聚类方法(Fuzzy c-means Clustering Method,FCM)是由Bezkek于1981年提出的,它是目前广泛采用的一种模糊聚类算法。考虑一个含有m个网格点的集合ndim代表网格点的维数,此集合依据一定的准则用模糊聚类的方法分成c(2≤c≤cmax)个模糊子集,这里c是指一给定的聚类个数,cmax为聚类个数的最大值,所用的准则是优化一个用来表征聚类的性能指标的目标函数Jcon(U,v),其基本形式为
其中,U为隶属度矩阵,U=[uik];uik表示第k个网格点属于第i个聚类的隶属度;v=(v1,v2,…,vc),vi是第i个聚类中心,1≤i≤c;con是一个大于1的常数(通常取2);dik是欧几里德距离,其数学表达式如公式(6)所示。
dik=||xk-vi|| (6)
其中,xk为第k个网格点。
使用拉格朗日乘子法将式(5)转化为无约束优化问题,如式(7)所示。
其中,为无约束的目标函数;λk为第k个网格点对应的拉格朗日乘子,求公式(7)的极小化问题,得公式(8):
由公式(8)可得:
对于高维设计问题(设计变量大于3),FCM方法在相似度与空间距离方面有很大的缺陷,由Bezdek和Dunn于1979年提出的一种基于最大似然估计距离的模糊聚类方法Gath-Geva clustering method(GGM)用来解决高维设计变量中优化设计的不足。该聚类与FCM聚类方式唯一不同的是距离函数不再是欧几里德距离,而是最大似然估计距离,具体形式为:
其中,det(Fωi)表示对Fωi取行列式;Fωi为第i个聚类的协方差,1≤i≤c,Fωi的值由公式(11)可得;exp()表示底数为e的指数函数。
其中,ω=2为加权幂指数。
发明内容
本发明的目的是为了解决传统的全局优化方法对飞行器高精度分析模型进行优化设计过程中存在计算分析耗时的缺点,提出了一种飞行器升力面结构优化设计方法。
本发明的一种飞行器升力面结构优化设计方法,具体实现步骤如下:
步骤1:给定飞行器升力面高精度分析模型、初始设计空间(用符号S(0)表示)、设计变量。
步骤2:利用拉丁超方试验设计方法(LHD)在初始设计空间S(0)中构造Ninitial个初始试验设计样本。初始试验设计样本的个数Ninitial由公式(12)确定,用符号N记录迭代次数,此时设置迭代次数N=1。
其中,n表示设计变量的个数,即设计空间的维数。
步骤3:通过调用步骤1给定的飞行器升力面高精度分析模型,计算/仿真当前试验设计样本所对应的飞行器升力面高精度分析模型的响应值,并将试验设计样本及其相对应的响应值存储到试验设计样本数据库中;所述试验设计样本数据库中包含的内容包括:通过拉丁超方试验设计方法所得到的试验设计样本及其相对应的飞行器升力面高精度分析模型响应值、每次迭代过程中的全局近似最优解及其相对应的飞行器升力面高精度分析模型响应值。
步骤4:利用步骤3中所述试验设计样本数据库中的所有试验设计样本及其相对应的飞行器升力面高精度分析模型响应值、每次迭代过程中的全局近似最优解及其相对应的飞行器升力面高精度分析模型响应值构造径向基函数代理模型。
步骤5:在第N次迭代(当前)的设计空间S(N-1)中,使用拉丁超方试验设计方法随机生成mgrid_points个网格点,网格点个数mgrid_points由公式(13)确定。用步骤4得到的径向基函数代理模型计算所有网格点对应的径向基函数代理模型响应值。给定整体空间缩减率M1%,将网格点对应的径向基函数代理模型响应值中较小的M1%网格点保留,删除其余网格点,从而完成对设计空间的第一次伪缩减。
mgrid_points=60n (13)
步骤6:给定聚类空间个数c,对经过步骤5的处理后保留下的网格点使用模糊聚类法进行搜索与分类,具体为:使用公式(14)对经过步骤5的处理后保留下的网格点进行聚类,依次得到第k个网格点属于第i个聚类的隶属度uik,根据隶属度将该网格点分成c个聚类空间(i=1,2,…,c),i=1,2,…,c。
其中,did为距离函数,con是一个大于1的常数,vi为第i个聚类的距离中心,m为网格点的个数, 表示向上取整,xk为第k个网格点。
当设计变量小于等于3时,距离函数dik通过公式(15)得到:
dik=||xk-vi|| (15)
当设计变量大于3时,距离函数dik通过公式(16)得到:
其中,det(Fωi)表示对Fωi取行列式;Fωi为第i个聚类的协方差,1≤i≤c,Fωi的值由公式(17)可得;exp()表示底数为e的指数函数。
其中,ω=2为加权幂指数。
步骤7:给定局部空间缩减率M2%,依次将c个聚类空间中径向基函数代理模型响应值较小的M2%个网格点保留,删除其余网格点,完成对设计空间的第二次伪缩减,用(i=1,2,…,c)表示每个聚类空间中删除(1-M2%)个网格点后得到的c个新聚类空间。
步骤8:采用遗传算法在每个新聚类空间对步骤4中所述径向基函数代理模型进行优化,得到第N次迭代中各个新聚类空间的近似最优点(用(i=1,2,…,c)表示)。将对应的径向基函数代理模型响应值中的最小值对应的近似最优点作为寻优过程中的当前全局近似最优解(用表示)。
步骤9:根据步骤8中所述与通过公式(18)定义两种欧氏距离,分别称为第一欧氏距离(用D0表示)和第二欧氏距离(用Di表示),由公式(19)得到第N次迭代的缩减因子(用λ表示);由公式(20)得到第N次迭代的控制因子(用σ表示)。用λ0表示初始缩减因子,其为人为设定值,λ0在区间[0,1]中;用σ0表示初始控制因子;如果λ<λ0并且σ>σ0,表明当前设计空间S(N-1)需要缩减,进行步骤10的操作;否则,进行步骤12的操作。
其中,||·||是二范数;i=1,2,…,c。
其中,Lmax为S(N-1)的最大边长,设计空间S(N-1)被认为是超矩形。
其中,V为空间S(N-1)的超体积,V0为初始空间S(0)的超体积。
步骤10:对公式(18)中距离Di(i=1,2,…,c)从小到大排序,得到新的距离D′i,如式(21)所示;然后给定需删除聚类数(用ndel表示,ndel<c),则距离D′i对应的聚类空间(i=c,c-1,…,c-ndel+1)则被删除。
D′i=sort(Di) (21)
其中,i=1,2,…,c,sort()表示从小到大排序函数。
步骤11:删除当前设计空间S(N-1),重组新的设计空间S(N)。将经过步骤10的操作后剩下的所有网格点的各维的最小值和最大值分别作为新的设计空间S(N)各维的最小值和最大值。
步骤12:获取径向基函数代理模型的当前全局近似最优解对应的飞行器升力面高精度分析模型响应值(用y*(N)表示),并将此全局近似最优解及其飞行器升力面高精度分析模型响应值y*(N)添加到试验设计样本数据库中,并设置迭代次数N自增1。
步骤13:通过公式(22)判断飞行器升力面高精度分析模型响应值y*(N)是否收敛,如果收敛则将步骤12中所述全局近似最优解作为飞行器升力面高精度分析模型的全局最优解输出,结束整个优化过程;否则,重复步骤4至步骤13,直到找到飞行器升力面高精度分析模型的全局最优解。
其中,为第N次全局近似最优解对应的径向基函数代理模型响应值;ε1、ε2为人为确定的小量,ε1取[0.0001,0.01]区间中的值,ε2=10ε1。
有益效果
本发明提出的飞行器升力面结构优化设计方法克服了传统的全局优化方法在求解飞行器升力面高精度分析模型设计问题时存在的计算耗时的缺点,能够有效的降低计算成本,提高优化效率,有助于缩短飞行器升力面高精度分析模型设计的周期。用更少的调用飞行器升力面高精度分析模型的次数得到精度更高的全局最优解。
附图说明
图1为本发明具体实施例中的一种飞行器升力面结构优化设计方法的操作流程图;
图2为本发明具体实施例中的飞行器升力面外形布局图;
图3为本发明具体实施例中的飞行器升力面内部结构图及设计变量分布图。
具体实施方式
为了更好的说明本发明的技术方案,下面结合附图和具体实施例,对本发明做进一步说明,并通过与传统飞行器升力面结构优化设计方法结果比较,对本发明的综合性能进行验证分析。
本实施例优化设计的目标是在满足颤振速度约束条件下飞行器升力面质量最小。对于高速飞行器来讲,颤振速度是一个很重要的非线性指标,当飞行器飞行速度达到颤振临界速度时将会引起非线性自激振动,从而使飞行器遭到破坏。
一种飞行器升力面结构优化设计方法,其操作流程如图1所示,具体为:
步骤1:给定飞行器升力面高精度分析模型、初始设计空间(用符号S(0)表示)、设计变量。
本实施例中,优化模型是低展弦比飞行器升力面,其气动外形如图2所示。其中,翼根长4米,翼梢长1.8米,前缘后掠角为34度,后缘后掠角为18度。内部结构分为蒙皮、翼肋与梁,将它们的厚度分别设为三个设计变量,范围为[0.5,3]毫米,建立有限元模型,共有956个壳单元和540个节点,如图3所示。其中,设计变量1为蒙皮厚度,设计变量2为翼肋厚度,设计变量3为梁厚度。飞行器结构材料使用铝,材料参数如表1示。
表1铝材料参数
结构质量由公式(23)得到:
其中,mass结构质量;为Si指第i个单元的面积;τi为第i个单元的厚度;ρ为材料密度。
颤振速度可由公式(24)得到:
其中,[M*]为模态质量矩阵;{ξ}为模态坐标;[K*]为模态刚度矩阵;ρ为空气密度;V为颤振速度;[A*]为广义的模态坐标气动力系数。假定临界颤振速度为1800m/s。
将飞行器升力面高精度分析模型转化为如下的数学模型:
min mass
st.Vf≥1800m/s (25)
0.5mm≤x1、x2 x3≤3mm
其中,Vf为颤振速度,x1,x2,x3为飞行器升力面蒙皮、翼肋和梁的厚度,为优化模型的设计变量。通过拉格朗日乘子法,将公式(25)转化为无约束优化问题如公式(26)所示。
YI=FI+0.5M((max(0,GI+μ/M))2-(μ/M)2) (26)
其中,FI,GI分别是结构质量mass与颤振速度约束方程(1800-Vf)的归一化。
步骤2:因为设计变量个数n=3,故在初始设计空间S(0)中,使用拉丁超方试验设计方法生成初始样本点个数Ninitial=6,并设置迭代次数N=1。
步骤3:通过调用步骤1给定的飞行器升力面高精度分析模型,计算/仿真当前试验设计样本所对应的飞行器升力面高精度分析模型的响应值,并将试验设计样本及其相对应的响应值存储到试验设计样本数据库中。计算一次该飞行器升力面高精度分析模型大约需要花50秒。
步骤4:利用步骤3中所述试验设计样本数据库中的所有试验设计样本及其相对应的飞行器升力面高精度分析模型响应值、每次迭代过程中的全局近似最优解及其相对应的飞行器升力面高精度分析模型响应值构造径向基函数代理模型。
步骤5:在第N次设计空间S(N-1)中,使用拉丁超方试验设计方法随机生成60×3=180个网格点。用步骤4得到的径向基函数代理模型计算所有网格点对应的径向基函数代理模型响应值。给定整体空间缩减率M1%=70%,将网格点对应的径向基函数代理模型响应值中较小的70%网格点保留,删除其余网格点,从而完成对设计空间的第一次伪缩减。
步骤6:给定聚类空间个数c=3,对经过步骤5的处理后保留下的网格点使用模糊聚类法进行搜索与分类,具体为:使用公式(14)对经过步骤5的处理后保留下的网格点进行聚类,依次得到第k个网格点属于第i个聚类的隶属度uik,根据隶属度将该网格点分成3个聚类空间(i=1,2,3)
步骤7:给定局部空间缩减率M2%=70%,依次将3个聚类空间中径向基函数代理模型响应值较小的70%个网格点保留,删除其余网格点,完成对设计空间的第二次伪缩减,得到3个新的聚类空间(i=1,2,3)。
步骤8:采用遗传算法在每个新聚类空间对步骤4中所述径向基函数代理模型进行优化,得到第N次迭代中各个新聚类空间的近似最优点(i=1,2,3)。将对应的径向基函数代理模型响应值中的最小值对应的近似最优点作为寻优过程中的当前全局近似最优解
步骤9:根据步骤8中所述与通过公式(18)定义两种欧氏距离,分别称为第一欧氏距离D0和第二欧氏距离Di,由公式(19)得到第N次迭代的缩减因子λ表示;由公式(20)得到第N次迭代的控制因子σ表示。给定初始缩减因子λ0=0.1与初始控制因子σ0=0.001,如果λ<λ0并且σ>σ0满足,表明当前设计空间S(N-1)需要缩减,进行步骤10的操作;否则,进行步骤12的操作。
步骤10:对公式(18)中距离Di(i=1,2,…,c)从小到大排序,得到新的距离D′i,如式(21)所示;给定需要删除的聚类个数ndel=1,因此,Ds=max(D1,D2,D3)对应的即被删除。
步骤11:删除当前设计空间S(N-1),重组新的设计空间S(N)。将经过步骤10的操作后剩下的所有网格点的各维的最小值和最大值分别作为新的设计空间S(N)各维的最小值和最大值。
步骤12:获取径向基函数代理模型的当前全局近似最优解对应的飞行器升力面高精度分析模型响应值y*(N),并将此全局近似最优解及其飞行器升力面高精度分析模型响应值y*(N)添加到试验设计样本数据库中,并设置迭代次数N自增1。
步骤13:通过公式(22)判断飞行器升力面高精度分析模型响应值y*(N)是否收敛,如果收敛则将步骤12中所述全局近似最优解作为飞行器升力面高精度分析模型的全局最优解输出,结束整个优化过程;否则,重复步骤4至步骤13,直到找到飞行器升力面高精度分析模型的全局最优解。公式(22)中的ε1=0.01,ε2=0.1。
经过上述步骤的操作,即完成对低展弦比飞行器升力面的优化,为了保证设计问题的鲁棒性,对该设计问题连续进行了11次优化,取11次中调用分析模型次数中的中间优化结果。同时,采用传统的飞行器优化方法序列二次规划(SQP)和遗传算法(GA)对该飞行器升力面进行优化,其优化结果如表2所示。
表2本发明方法、SQP、GA升力面优化结果对比表
由表2可知,由于该模型的特殊性,序列二次规划(SQP)不能对该有限元模型进行优化,不能进行寻优,颤振速度不能满足设计要求;遗传算法(GA)在调用了758次有限元分析模型后得到了目标结构质量为81.70kg,每次调用分析模型大约需要50s,因此遗传算法优化所消耗的总时间为50×758=37900sec=10.53h。而本发明只调用了69次分析模型,目标结构质量为72.38kg,所消耗的总时间为50×69=3450sec=57.5min。显然,本发明用了大约是遗传算法十分之一的时间,得到了目标更小的结构质量,大大提高了优化效率,节约设计成本。
根据飞行器升力面优化实例分析可见,相比于传统的优化设计方法,本发明有助于提高飞行器升力面优化设计结果与设计质量;另一方面,涉及飞行器升力面高精度分析模型的优化问题,本发明还能大大提高的优化效率,降低飞行器优化设计成本,缩短优化设计周期。
以上所述的具体描述,对发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例,用于解释本发明,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种飞行器升力面结构优化设计方法,其特征在于:其实现过程包括步骤1至步骤13,具体为:
步骤1:给定飞行器升力面高精度分析模型、初始设计空间S(0)、设计变量;
步骤2:利用拉丁超方试验设计方法在初始设计空间S(0)中构造初始试验设计样本;用符号N记录迭代次数,此时设置迭代次数N=1;
步骤3:通过调用步骤1给定的飞行器升力面高精度分析模型,计算/仿真当前试验设计样本所对应的飞行器升力面高精度分析模型的响应值,并将试验设计样本及其相对应的响应值存储到试验设计样本数据库中;所述试验设计样本数据库中包含的内容包括:通过拉丁超方试验设计方法所得到的试验设计样本及其相对应的飞行器升力面高精度分析模型响应值、每次迭代过程中的全局近似最优解及其相对应的飞行器升力面高精度分析模型响应值;
步骤4:利用步骤3中所述试验设计样本数据库中的所有试验设计样本及其相对应的飞行器升力面高精度分析模型响应值、每次迭代过程中的全局近似最优解及其相对应的飞行器升力面高精度分析模型响应值构造径向基函数代理模型;
步骤5:在第N次迭代的设计空间S(N-1)中,使用拉丁超方试验设计方法随机生成网格点;用步骤4得到的径向基函数代理模型计算所有网格点对应的径向基函数代理模型响应值;给定整体空间缩减率M1%,将网格点对应的径向基函数代理模型响应值中较小的M1%网格点保留,删除其余网格点;
步骤6:给定聚类空间个数c,对经过步骤5的处理后保留下的网格点使用模糊聚类法进行搜索与分类;
步骤7:给定局部空间缩减率M2%,依次将c个聚类空间中径向基函数代理模型响应值较小的M2%个网格点保留,删除其余网格点,用表示每个聚类空间中删除(1-M2%)个网格点后得到的c个新聚类空间;
步骤8:采用遗传算法在每个新聚类空间对步骤4中所述径向基函数代理模型进行优化,得到第N次迭代中各个新聚类空间的近似最优点将对应的径向基函数代理模型响应值中的最小值对应的近似最优点作为寻优过程中的当前全局近似最优解
步骤9:根据步骤8中所述与通过公式(18)定义两种欧氏距离,分别称为第一欧氏距离D0和第二欧氏距离Di,由公式(19)得到第N次迭代的缩减因子λ;由公式(20)得到第N次迭代的控制因子σ;用λ0表示初始缩减因子,其为人为设定值,λ0在区间[0,1]中;用σ0表示初始控制因子;如果λ<λ0并且σ>σ0,表明当前设计空间S(N-1)需要缩减,进行步骤10的操作;否则,进行步骤12的操作;
其中,||·||是二范数;i=1,2,…,c;
其中,Lmax为S(N-1)的最大边长,设计空间S(N-1)被认为是超矩形;
其中,V为空间S(N-1)的超体积,V0为初始空间S(0)的超体积;
步骤10:对公式(18)中距离Di(i=1,2,…,c)从小到大排序,得到新的距离D′i,如式(21)所示;然后给定需删除聚类数ndel,ndel<c,距离D′i对应的聚类空间被删除;
D′i=sort(Di) (21)
其中,i=1,2,…,c,sort()表示从小到大排序函数;
步骤11:删除当前设计空间S(N-1),重组新的设计空间S(N);将经过步骤10的操作后剩下的所有网格点的各维的最小值和最大值分别作为新的设计空间S(N)各维的最小值和最大值;
步骤12:获取径向基函数代理模型的当前全局近似最优解对应的飞行器升力面高精度分析模型响应值y*(N),并将此全局近似最优解及其飞行器升力面高精度分析模型响应值y*(N)添加到试验设计样本数据库中,并设置迭代次数N自增1;
步骤13:通过公式(22)判断飞行器升力面高精度分析模型响应值y*(N)是否收敛,如果收敛则将步骤12中所述全局近似最优解作为飞行器升力面高精度分析模型的全局最优解输出,结束整个优化过程;否则,重复步骤4至步骤13,直到找到飞行器升力面高精度分析模型的全局最优解;
其中,为第N次全局近似最优解对应的径向基函数代理模型响应值;ε1、ε2为人为确定的小量。
2.如权利要求1所述的一种飞行器升力面结构优化设计方法,其特征在于:其步骤2中所述初始试验设计样本的个数Ninitial由公式(12)确定;
其中,n表示设计变量的个数,即设计空间的维数。
3.如权利要求1或2所述的一种飞行器升力面结构优化设计方法,其特征在于:其步骤5中所述网格点个数mgrid_points由公式(13)确定;
mgrid_points=60n (13)
其中,n表示设计变量的个数。
4.如权利要求1或2所述的一种飞行器升力面结构优化设计方法,其特征在于:其步骤6中所述对经过步骤5的处理后保留下的网格点使用模糊聚类法进行搜索与分类的具体方法为:使用公式(14)对经过步骤5的处理后保留下的网格点进行聚类,依次得到第k个网格点属于第i个聚类的隶属度μik,根据隶属度将该网格点分成c个聚类空间 i=1,2,…,c;
其中,dik为距离函数,con是一个大于1的常数,vi为第i个聚类的距离中心,m为网格点的个数, 表示向上取整,xk为第k个网格点;
当设计变量小于等于3时,距离函数dik通过公式(15)得到:
dik=||xk-vi|| (15)
当设计变量大于3时,距离函数dik通过公式(16)得到:
其中,det(Fωi)表示对Fωi取行列式;Fωi为第i个聚类的协方差,1≤i≤c,Fωi的值由公式(17)可得;exp()表示底数为e的指数函数;
其中,ω=2为加权幂指数。
5.如权利要求1或2所述的一种飞行器升力面结构优化设计方法,其特征在于:其步骤13中所述ε1取[0.0001,0.01]区间中的值,ε2=10ε1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210158035.9A CN102789539B (zh) | 2012-05-21 | 2012-05-21 | 一种飞行器升力面结构优化设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210158035.9A CN102789539B (zh) | 2012-05-21 | 2012-05-21 | 一种飞行器升力面结构优化设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102789539A CN102789539A (zh) | 2012-11-21 |
CN102789539B true CN102789539B (zh) | 2015-05-13 |
Family
ID=47154942
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210158035.9A Active CN102789539B (zh) | 2012-05-21 | 2012-05-21 | 一种飞行器升力面结构优化设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102789539B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103149841A (zh) * | 2013-02-25 | 2013-06-12 | 西北工业大学 | 降低俯仰气动非线性特性的飞行器总体随控优化方法 |
CN104050338A (zh) * | 2014-07-07 | 2014-09-17 | 北京理工大学 | 一种geo卫星小推力器推进剂消耗耦合分析方法 |
CN106777727A (zh) * | 2016-12-26 | 2017-05-31 | 电子科技大学 | 一种可穿戴电子设备系统多学科设计优化的方法 |
CN107092710B (zh) * | 2017-02-27 | 2018-09-14 | 北京航空航天大学 | 一种基于超体积迭代全局优化算法的确定最优结构尺寸的方法 |
CN107451316B (zh) * | 2017-04-01 | 2019-06-18 | 南京航空航天大学 | 一种卫星的构型设计方法 |
CN107357996B (zh) * | 2017-07-17 | 2019-01-08 | 北京航空航天大学 | 一种基于超体积迭代策略的代理模型试验设计方法 |
CN107515994B (zh) * | 2017-09-11 | 2020-04-03 | 西北工业大学 | 一种自主水下航行器的壳体结构多保真度设计优化方法 |
CN108363310B (zh) * | 2018-01-15 | 2020-09-08 | 北京航空航天大学 | 人工智能程序员书写数字飞行器源代码的聚类决策方法 |
CN108459993B (zh) * | 2018-02-02 | 2021-01-05 | 北京理工大学 | 基于快速追峰采样的复杂高维系统优化方法 |
CN110442886B (zh) * | 2018-05-03 | 2023-05-16 | 中国航天科工飞航技术研究院(中国航天海鹰机电技术研究院) | 一种悬浮试验装置及设计方法 |
CN109236258B (zh) * | 2018-10-27 | 2019-07-23 | 中国地质大学(北京) | 一种基于自适应代理模型的致密油藏压裂水平井优化方法 |
CN110263429B (zh) * | 2019-06-18 | 2021-11-19 | 西安电子科技大学 | 基于机器学习的索网反射面天线形面精度调整方法 |
CN114329791B (zh) * | 2021-12-31 | 2024-06-07 | 北京航空航天大学 | 基于模块集成与数据管理的飞行器机翼结构综合优化方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101944141A (zh) * | 2010-08-18 | 2011-01-12 | 北京理工大学 | 一种高效的基于模糊聚类自适应径向基全局优化方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100219296A1 (en) * | 2008-11-01 | 2010-09-02 | Alexander J. Shelman-Cohen | Reduced drag system for windmills, fans, propellers, airfoils, and hydrofoils |
-
2012
- 2012-05-21 CN CN201210158035.9A patent/CN102789539B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101944141A (zh) * | 2010-08-18 | 2011-01-12 | 北京理工大学 | 一种高效的基于模糊聚类自适应径向基全局优化方法 |
Non-Patent Citations (2)
Title |
---|
Yanping LIAO,Li LIU,Teng Long.《Multi-Objective Aerodynamic and Stealth Performance Optimization for Airfoil Using Kriging Surrogate Model》.《2011 IEEE 3rd international Conference on Communication Software and Networks》.2011,全文. * |
彭磊,刘莉,龙腾.《基于动态径向基函数代理模型的优化策略》.《机械工程学报》.2011,第47卷(第7期), * |
Also Published As
Publication number | Publication date |
---|---|
CN102789539A (zh) | 2012-11-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102789539B (zh) | 一种飞行器升力面结构优化设计方法 | |
Wu et al. | Data-driven reduced order model with temporal convolutional neural network | |
CN104866692B (zh) | 一种基于自适应代理模型的飞行器多目标优化方法 | |
Jun et al. | Application of a PCA-DBN-based surrogate model to robust aerodynamic design optimization | |
Yasong et al. | Global aerodynamic design optimization based on data dimensionality reduction | |
CN102682173A (zh) | 基于自适应径向基函数代理模型的飞行器优化设计方法 | |
Han et al. | Surrogate-based aerodynamic shape optimization with application to wind turbine airfoils | |
Li et al. | Vectorial surrogate modeling approach for multi-failure correlated probabilistic evaluation of turbine rotor | |
CN115564114B (zh) | 一种基于图神经网络的空域碳排放短期预测方法及系统 | |
Liu et al. | Prediction and optimization of airfoil aerodynamic performance using deep neural network coupled Bayesian method | |
Scarth et al. | Reliability-based aeroelastic design of composite plate wings using a stability margin | |
CN105718634A (zh) | 一种基于非概率区间分析模型的翼型鲁棒优化设计方法 | |
CN111832839B (zh) | 基于充分增量学习的能耗预测方法 | |
You et al. | A novel mixed uncertainty support vector machine method for structural reliability analysis | |
Sun et al. | Shape optimization of underwater wings with a new multi-fidelity bi-level strategy | |
Tang et al. | Aerodynamic robustness optimization and design exploration of centrifugal compressor impeller under uncertainties | |
Wang et al. | A nacelle inlet design approach with more three-dimensional geometric consideration | |
Peng et al. | Multi-fidelity nonlinear unsteady aerodynamic modeling and uncertainty estimation based on Hierarchical Kriging | |
Luo et al. | Aerodynamic optimization of a transonic fan rotor by blade sweeping using adaptive Gaussian process | |
Lei et al. | An efficient method by nesting adaptive Kriging into Importance Sampling for failure-probability-based global sensitivity analysis | |
Wang et al. | An efficient geometric constraint handling method for surrogate-based aerodynamic shape optimization | |
Tang et al. | Adaptive dynamic surrogate-assisted evolutionary computation for high-fidelity optimization in engineering | |
Liu et al. | Uncertainty optimization design of airfoil based on adaptive point adding strategy | |
Rumpfkeil et al. | Construction of multi-fidelity locally optimized surrogate models for uncertainty quantification | |
Zhu et al. | Design of an RBF surrogate model for low Reynolds number airfoil based on transfer learning |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |