CN115774900B - 一种不确定条件下的变构型飞行器指令鲁棒优化设计方法 - Google Patents
一种不确定条件下的变构型飞行器指令鲁棒优化设计方法 Download PDFInfo
- Publication number
- CN115774900B CN115774900B CN202211444860.5A CN202211444860A CN115774900B CN 115774900 B CN115774900 B CN 115774900B CN 202211444860 A CN202211444860 A CN 202211444860A CN 115774900 B CN115774900 B CN 115774900B
- Authority
- CN
- China
- Prior art keywords
- aircraft
- variable
- model
- uncertainty
- configuration
- 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
- 238000000034 method Methods 0.000 title claims abstract description 84
- 238000005457 optimization Methods 0.000 title claims abstract description 56
- 238000013461 design Methods 0.000 title claims abstract description 25
- 230000008569 process Effects 0.000 claims abstract description 21
- 238000010206 sensitivity analysis Methods 0.000 claims abstract description 8
- 238000004088 simulation Methods 0.000 claims abstract description 8
- 238000005070 sampling Methods 0.000 claims description 13
- 239000011159 matrix material Substances 0.000 claims description 11
- 238000004458 analytical method Methods 0.000 claims description 10
- 230000000739 chaotic effect Effects 0.000 claims description 9
- 230000008878 coupling Effects 0.000 claims description 8
- 238000010168 coupling process Methods 0.000 claims description 8
- 238000005859 coupling reaction Methods 0.000 claims description 8
- 230000014509 gene expression Effects 0.000 claims description 6
- 230000004044 response Effects 0.000 claims description 6
- 230000035945 sensitivity Effects 0.000 claims description 6
- 230000003416 augmentation Effects 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 5
- 239000003795 chemical substances by application Substances 0.000 claims description 5
- 238000005315 distribution function Methods 0.000 claims description 5
- 230000001419 dependent effect Effects 0.000 claims description 4
- 101150036841 minJ gene Proteins 0.000 claims description 4
- 238000007476 Maximum Likelihood Methods 0.000 claims description 3
- 238000000342 Monte Carlo simulation Methods 0.000 claims description 3
- 238000012098 association analyses Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 238000013178 mathematical model Methods 0.000 claims description 3
- 238000005381 potential energy Methods 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 3
- 230000003190 augmentative effect Effects 0.000 claims description 2
- 230000001186 cumulative effect Effects 0.000 claims description 2
- 230000005684 electric field Effects 0.000 claims description 2
- 238000011156 evaluation Methods 0.000 claims description 2
- 238000013507 mapping Methods 0.000 claims description 2
- 230000035939 shock Effects 0.000 claims description 2
- 238000010200 validation analysis Methods 0.000 claims description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims 1
- 125000003275 alpha amino acid group Chemical group 0.000 claims 1
- 150000001875 compounds Chemical class 0.000 claims 1
- 238000001228 spectrum Methods 0.000 claims 1
- 239000000126 substance Substances 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 230000035899 viability Effects 0.000 description 1
Classifications
-
- 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
- Feedback Control In General (AREA)
Abstract
本发明公开了一种不确定条件下的变构型飞行器指令鲁棒优化设计方法,包括如下步骤:设计飞行器变构型调节策略,构建飞行器动力学模型,采用模型敏感度分析方法,确定飞行器变构型调节的代理模型;量化飞行过程中的不确定项,表征飞行器模型不确定性的传播方式,推演出变构型飞行器不确定性的具体形式;确定变构型飞行器优化性能指标,构建出变构型飞行器模型特性对闭环极限性能的约束关系,设计变构型飞行器指令的分层迭代策略,解算出不确定性条件下变构型飞行器的任务指令;采用数值仿真验证变构型飞行器所获得任务指令的有效性。本发明能够量化不确定性的影响,满足变构型飞行器不确定条件下的复杂任务要求。
Description
技术领域
本发明涉及变构型飞行器设计技术领域,尤其是一种不确定条件下的变构型飞行器指令鲁棒优化设计方法。
背景技术
变构型飞行器运行环境更为复杂,难以精确描述其动力学特性,加之各种随机干扰使得模型具有大的不确定性。更为重要的是,变构型飞行器受到结构承载和推进效率的限制,需要对飞行状态尤其是攻角、侧滑角等进行严格约束,同时还应满足舵面饱和、热壅塞等多重约束条件,亟待突破不确定条件下变构型飞行器任务指令的鲁棒优化技术,为工程应用提供有价值的设计方法。
变构型飞行器任务包线更大,系统存在明显的多约束强耦合关系,同时飞行器受到外界强不确定扰动的影响,模型参数易变,标称情况下的优化结果往往不满足于实际的性能需求。
发明内容
本发明所要解决的技术问题在于,提供一种不确定条件下的变构型飞行器指令鲁棒优化设计方法,量化不确定性的影响,将变构模型参数的不确定性融入到优化设计中形成迭代,获得飞行器的任务指令,满足变构型飞行器不确定条件下的复杂任务要求。
为解决上述技术问题,本发明提供一种不确定条件下的变构型飞行器指令鲁棒优化设计方法,包括如下步骤:
步骤1、设计飞行器变构型调节策略,构建飞行器动力学模型,采用模型敏感度分析方法,确定飞行器变构型调节的代理模型;
步骤2、量化飞行过程中的不确定项,表征飞行器模型不确定性的传播方式,推演出变构型飞行器不确定性的具体形式;
步骤3、确定变构型飞行器优化性能指标,构建出变构型飞行器模型特性对闭环极限性能的约束关系,设计变构型飞行器指令的分层迭代策略,解算出不确定性条件下变构型飞行器的任务指令;
步骤4、采用数值仿真验证变构型飞行器所获得任务指令的有效性。
优选的,步骤1中,设计飞行器变构型调节策略,构建飞行器动力学模型,采用模型敏感度分析方法,确定飞行器变构型调节的代理模型具体包括如下步骤:
步骤11、针对飞行器构型,给出飞行过程中的动压、过载和热流多约束条件,采用工程估算方法估算飞行器的气动力,对飞行器几何外形进行网格划分,生成若干面元,据飞行状态,面元参数以及激波/膨胀波理论获得的各面元的来流参数,以及面元距前缘距离和非稳态速度面元数据,得到每个面元上总气动力和总压,对所有面元的力和力矩进行求和,计算出最终的气动力和力矩;
确定飞行器飞行过程的受力与力矩同变构型参数与飞行状态之间的函数关系,结合拉格朗日方程、虚功原理和密歇尔斯基等式,构建变构型飞行器的动力学模型,
式中,V和ω分别为飞行器的线速度和角速度,m和J分别表示质量和转动惯量矩阵,F和M为飞行器的气动力和力矩,δ是构型参数;
步骤12、针对构建的变构型飞行器模型数据库和非线性数学模型,采用拉丁超方采样方法,优选模型数据库中合理的样本点,通过改变数据排列顺序提高采样性能,进而对样本势能进行优化,引入物理中电场电势的含义,通过整个样本空间内的电势大小Ee进行评价
其中:Ns为样本空间的容量,si,sj分别表示第i个和第j个样本点,基于电势大小的输出值,实现最优正交拉丁超方采样;
步骤13、针对采样出的变构型飞行器模型数据,采用灵敏度分析方法,确定特征向量灵敏度矩阵SE,根据灵敏度矩阵中的元素幅值判断表征系统状态与运动模态间的关系,提炼与飞行器飞行过程相关的主要影响参数,采用极大似然辨识方法,选取使似然函数L(θ)达到最大的参数值θ,作为参数θ的估计值,确定变构型飞行器代理模型结构,将飞行器复杂的气动力、力矩和推力表达式表示成关于飞行状态与变构型调节参数的函数,推演出变构型飞行器代理模型形式,以升力系数CL和阻力系数CD为例,表示为:
其中为升力相关的模型系数,/>为阻力相关的模型系数,这些模型系数是构型参数δ的函数;
步骤14、评估变构型飞行器不同速域范围内模型动态特性之间的相似性,依据相似性准则整合不同区域的模型,将飞行器算力模型表达式的待定系数表示成变构型调节参数的函数,确定融合变构调节参数的线性变参模型,
其中,x为状态变量;y为输出变量;u为输入变量;状态矩阵A,B以及输入矩阵C,D会随着高度h与马赫数Ma而发生变化;
步骤15、采用模态分析和闭环极限理论和方法,当特征向量灵敏度矩阵一列中的元素小于该列最大元素的10%时,便将该元素对应的状态视为无关状态,对变构型飞行器动力学模型的平衡状态和动态特性进行折衷分析,静态特性包括配平状态和配平输入等,动态特性包括系统稳定性、零极点位置,估算飞行器的可控安全边界。
优选的,步骤2中,量化飞行过程中的不确定项,表征飞行器模型不确定性的传播方式,推演出变构型飞行器不确定性的具体形式具体包括如下步骤:
步骤21、采用广义混沌多项式量化飞行过程中的不确定性,基于伽辽金理论将混沌多项式表示的不确定扰动引入到飞行器模型中,得到不确定条件下的混沌多项式展开
其中,i=1,2,…,NY,NY是输出响应Yi中变量的个数,ε为随机气动摄动变量(ε=[ε1,…,εd],d为随机变量的维数),为第i个输出响应的多项式展开系数,而Ψj为第i个输出响应的混沌多项式;对于存在多维和混合分布的不确定项,如果它们是相互独立的,则概率密度分布为每个随机变量概率密度分布函数的乘积,如果不确定项存在两个相互关联的分量,采用纳塔夫方法进行转换,以随机变量的边缘累积分布函数和协方差矩阵为基本信息,将相关非正态随机变量直接转换为独立标准正态随机变量;
步骤22、采用截断误差作为指标,利用谱分解等理论自适应更新混沌多项展开式的基函数,当截断误差超过阈值时,混沌多项展开式使用的基将自适应地更新,将变构型飞行器考虑不确定性的模型组转换为等价的增广模型形式,对对随机变量进行采样得到有
其中,x为飞行状态,s=1,2,…,ns,ns为采样点个数,从而将不确定性建模问题转变为确定的增广建模问题;
步骤23、将变构型飞行段的多约束条件转换为增广模型的等价约束,采用能量状态法获得多约束条件下变构型飞行器的可行区域,将速度、高度、轨迹角和质量作为状态变量,将迎角和舵面作为控制变量,总能量可以看作是系统连续的状态变量,以飞行器可行区域边界对应的是飞行器能量状态变化率为0,通过计算能量状态变化率来确定飞行器可行区域。
优选的,步骤3中,确定变构型飞行器优化性能指标,构建出变构型飞行器模型特性对闭环极限性能的约束关系,设计变构型飞行器指令的分层迭代策略,解算出不确定性条件下变构型飞行器的任务指令具体包括如下步骤:
步骤31、根据最小期望/协方差,结合指标加权策略,确立多约束条件下变构型飞行器指令优化的性能指标;具体为:分析变构型飞行器模型耦合特性,综合稳定性、控制精度、鲁棒性要求,根据最小期望/协方差,结合指标加权策略,确立多约束条件下飞行器指令优化的性能指标,描述如下:
式(1)中,是优化的目标函数,/>由状态变量组成的向量,/>和/>给出了状态变量的约束范围,p是需要优化的指令,考虑到不确定性会存在于状态变量/>中或目标函数的输出/>中,采用鲁棒优化方法,构造关于原优化目标函数f的适应度函数F,即
式(2)中,Fn为基于适应度函数鲁棒性能指标,N是不确定性参数的数量,设计变量的约束不变,而和/>描述了不确定性的约束,/>是不确定性参数p的所有可能情况的集合,由pub和plb分别定义了上界和下界,/>是一个在/>中满足一定分布规律Λ(plbn,pubn)的随机参数,采用鲁棒正则化方法,构造如下的适应度函数:
式(3)中,δ是不确定性影响下指令关于标称位置的偏移量,而φ(δ)是飞行器指令关于标称位置偏移量的概率密度函数;
步骤32、根据飞行器不确定性与性能指标的关联分析,选取性能指标的权重系数,采用概率密度函数衡量不同类型不确定项对指标的影响程度,确定出的性能指标为
minJ=ωμJμ+ωσJσ (12)
其中,Jμ和Jσ为性能指标的均值和方差对应的权重值;
步骤33、通过分析变构型飞行器性能要求之间的耦合特性,构建指令的分层优化策略,优化问题可以表示如下:
minJ(P,K,r) (13)
其中,J表示目标函数,优化的目标是使得目标函数值取得最小。P,K,r为设计变量,即在优化设计过程中可以被用户控制的变量:P为构型参数,K为输入指令参数,r为状态参数,分层优化是将优化问题进行拆分若干个子系统,首先完成对一个子系统的优化,之后将最优的参数代入到下一步进行优化,获得不确定性条件下满足多约束条件的飞行器任务指令。
优选的,步骤4中,采用数值仿真验证变构型飞行器所获得任务指令的有效性具体包括如下步骤:
步骤41、针对给定标称状态下,对任务指令进行分析,评估标称状态下任务优化的有效性,构建飞行器构型参数与闭环性能可达边界的映射关系,通过改变构型参数,分析闭环性能,验证标称状态下变构型设计的有效性;
步骤42、在强不确定的严苛条件下,对变构型飞行器模型参数进行拉偏,通过蒙脱卡洛仿真分析,依靠大量重复的仿真实验来,分析不同输入条件下,系统输出的变化范围,解算飞行性能是否满足期望的要求,从而验证控制性能增强设计结果的有效性。
本发明的有益效果为:采用鲁棒优化方法可使飞行器在执行不同的使命,以及某一使命的各个任务段都具有满意的飞行效能,因而可有效增加航程和节省能耗,同时能够执行多种使命,具有更高的环境适应、应变、攻击和生存能力。
附图说明
图1为本发明的变构型飞行器建模流程示意图。
图2为本发明的变构型飞行器不确定量化与包线走廊确定流程示意图。
图3为本发明的变构型飞行器指令的鲁棒优化设计流程示意图。
具体实施方式
一种不确定条件下的变构型飞行器指令鲁棒优化设计方法,包括如下步骤:
步骤1、设计飞行器变构型调节策略,构建飞行器动力学模型,采用模型敏感度分析方法,确定飞行器变构型调节的代理模型。如图1所示,具体包括如下步骤:
(1)针对飞行器典型构型,根据任务需求,设计典型任务的变构型策略,给出飞行过程中的动压、过载和热流等多约束条件,采用工程估算方法构建飞行器的算力数据库,确定飞行器飞行过程中的受力与力矩同变构型参数与飞行状态之间的函数关系,结合拉格朗日方程、虚功原理和密歇尔斯基等式,构建变构型飞行器的动力学模型。
(2)针对构建的变构型飞行器模型数据库和非线性数学模型,采用拉丁超方采样方法优选模型数据库中合理的样本点通过改变数据排列顺序提高采样性能,进而对样本势能进行优化,实现最优正交拉丁超方采样。
(3)针对采样出的变构型飞行器模型数据,采用灵敏度分析方法,提炼与飞行器飞行过程相关的主要影响参数。采用极大似然辨识方法,确定变构型飞行器代理模型结构,将飞行器复杂的气动力、力矩和推力表达式表示成关于飞行状态与变结构调节参数的函数,推演出变构型飞行器代理模型形式。
(4)评估变构型飞行器不同包线范围内模型动态特性之间的相似性,依据相似性准则整合不同区域的模型,将飞行器算力模型表达式的待定系数表示成变构型调节参数的函数,确定融合变构型调节参数的代理形式。
(5)采用模态分析和闭环极限理论和方法,对变构型飞行器动力学模型的平衡状态和动态特性进行折衷分析,考虑飞行条件、变构型调节和不确定性因素对静态平衡特性(飞行包线、配平舵面、迎角和推力)及动态特性(不稳定极点和右半平面零点)影响,结合对象模型的输入和状态约束条件,对系统特性随飞行状态的变化规律以及闭环性能极限进行了分析,得到飞行器的可控安全边界。
步骤2、量化飞行过程中的不确定项,表征飞行器模型不确定性的传播方式,推演出变构型飞行器不确定性的具体形式。
该步骤结合混沌多项式扩展法和投影法,对变构型飞行器运动过程存在的不确定性进行量化,应用多元混沌多项式,将飞行器模型的不确定性转化为其等价确定性矩阵,进而采用能量状态法确定变构型飞行器的包线走廊,具体的实施过程如图2所示,包括如下步骤:
(1)采用广义混沌多项式量化飞行过程中的不确定性,基于伽辽金理论将混沌多项式表示的不确定扰动引入到变构型飞行器模型中,对于存在多维和混合分布的不确定项,如果它们是相互独立的,则概率密度分布为每个随机变量概率密度分布函数的乘积,如果不确定项存在两个相互关联的分量,可以采用纳塔夫方法进行转换。
(2)采用截断误差作为指标,利用谱分解等理论自适应更新混沌多项展开式的基函数,将变构型飞行器考虑不确定性的模型组转换为等价的增广模型形式。
(3)将变构型飞行过程中的多约束条件转换为增广模型的等价约束,采用能量状态法获得不确定条件下变构型飞行器的包线走廊。
步骤3、确定变构型飞行器优化性能指标,构建出变构型飞行器模型特性对闭环极限性能的约束关系,设计变构型飞行器指令的分层迭代策略,解算出不确定性条件下变构型飞行器的任务指令。
具体的实施过程如图3所示,包括如下步骤:
(1)分析变构型飞行器模型耦合特性,综合稳定性、控制精度、鲁棒性等要求,根据最小期望/协方差,结合指标加权策略,确立多约束条件下飞行器指令优化的性能指标,可以描述如下:
式(1)中,是优化的目标函数,/>由状态变量组成的向量,/>和/>给出了状态变量的约束范围,p是需要优化的指令。考虑到不确定性会存在于状态变量/>中或目标函数的输出/>中,采用鲁棒优化方法,构造关于原优化目标函数f的适应度函数F,即
式(2)中,Fn为基于适应度函数鲁棒性能指标,N是不确定性参数的数量。设计变量的约束不变,而和/>描述了不确定性的约束。/>是不确定性参数p的所有可能情况的集合,由pub和plb分别定义了上界和下界,/>是一个在/>中满足一定分布规律Λ(plbn,pubn)的随机参数。采用鲁棒正则化方法,构造如下的适应度函数:
式(3)中,δ是不确定性影响下指令关于标称位置的偏移量,而φ(δ)是飞行器指令关于标称位置偏移量的概率密度函数。
(2)根据飞行器运动过程中的不确定性与性能指标的关联分析,选取性能指标的权重系数,采用概率密度函数衡量不同类型不确定项对指标的影响程度。概率密度函数是对不确定性传播的一种量化,它反映了不确定项在当前时刻的整体分布情况,预测不确定性的传播。
(3)通过分析飞行器性能要求之间的耦合特性,构建指令的分层优化策略,获得不确定性条件下满足多约束条件的飞行器任务指令。变构型飞行器任务指令的鲁棒优化策略,通过耦合特性分析与优化策略的选取,降低变构型飞行器多变量多目标优化问题的复杂性,解决优化过程中的迭代效率低甚至无解的问题,提高不确定性条件下飞行器指令获取效率。
步骤4、采用数值仿真验证变构型飞行器所获得任务指令的有效性。具体包括如下步骤:
(1)针对给定标称状态下,对任务指令进行分析,评估标称状态下任务优化的有效性。
(2)在强不确定的严苛条件下,对变构型飞行器模型参数进行拉偏,通过蒙特卡洛仿真分析,验证所获得结果的有效性。
本发明给出一种不确定条件下变构型飞行器的鲁棒优化设计方法,量化不确定性的影响,将变构模型参数的不确定性融入到优化设计中形成迭代,获得飞行器的任务指令,满足变构型飞行器不确定条件下的复杂任务要求。
Claims (2)
1.一种不确定条件下的变构型飞行器指令鲁棒优化设计方法,其特征在于,包括如下步骤:
步骤1、设计飞行器变构型调节策略,构建飞行器动力学模型,采用模型敏感度分析方法,确定飞行器变构型调节的代理模型;具体包括如下步骤:
步骤11、针对飞行器构型,给出飞行过程中的动压、过载和热流多约束条件,采用工程估算方法估算飞行器的气动力,对飞行器几何外形进行网格划分,生成若干面元,据飞行状态、面元参数以及激波/膨胀波理论获得的各面元的来流参数,以及面元距前缘距离和非稳态速度面元数据,得到每个面元上总气动力和总压,对所有面元的力和力矩进行求和,计算出最终的气动力和力矩;
确定飞行器飞行过程的受力与力矩同变构型参数与飞行状态之间的函数关系,结合拉格朗日方程、虚功原理和密歇尔斯基等式,构建变构型飞行器的动力学模型,
式中,V和ω分别为飞行器的线速度和角速度,m和J分别表示质量和转动惯量矩阵,F和M为飞行器的气动力和力矩,δ是构型参数;
步骤12、针对构建的变构型飞行器模型数据库和非线性数学模型,采用拉丁超方采样方法,选取模型数据库中合理的样本点,通过改变数据排列顺序提高采样性能,进而对样本势能进行优化,引入物理中电场电势的含义,通过整个样本空间内的电势大小Ee进行评价
其中:Ns为样本空间的容量,si,sj分别表示第i个和第j个样本点,基于电势大小的输出值,实现最优正交拉丁超方采样;
步骤13、针对采样出的变构型飞行器模型数据,采用灵敏度分析方法,确定特征向量灵敏度矩阵SE,根据灵敏度矩阵中的元素幅值判断表征系统状态与运动模态间的关系,提炼与飞行器飞行过程相关的主要影响参数,采用极大似然辨识方法,选取使似然函数L(θ)达到最大的参数值θ,作为参数θ的估计值,确定变构型飞行器代理模型结构,将飞行器复杂的气动力、力矩和推力表达式表示成关于飞行状态与变构型调节参数的函数,推演出变构型飞行器代理模型形式,升力系数CL和阻力系数CD表示为:
其中为升力相关的模型系数,/>为阻力相关的模型系数,这些模型系数是构型参数δ的函数;
步骤14、评估变构型飞行器不同速域范围内模型动态特性之间的相似性,依据相似性准则整合不同区域的模型,将飞行器算力模型表达式的待定系数表示成变构型调节参数的函数,确定融合变构调节参数的线性变参模型,
其中,x为状态变量;y为输出变量;u为输入变量;状态矩阵A,B以及输入矩阵C,D会随着高度h与马赫数Ma而发生变化;
步骤15、采用模态分析和闭环极限理论和方法,当特征向量灵敏度矩阵一列中的元素小于该列最大元素的10%时,便将该元素对应的状态视为无关状态,对变构型飞行器动力学模型的平衡状态和动态特性进行折衷分析,静态特性包括配平状态和配平输入,动态特性包括系统稳定性、零极点位置,估算飞行器的可控安全边界;
步骤2、量化飞行过程中的不确定项,表征飞行器模型不确定性的传播方式,推演出变构型飞行器不确定性的具体形式;具体包括如下步骤:
步骤21、采用广义混沌多项式量化飞行过程中的不确定性,基于伽辽金理论将混沌多项式表示的不确定扰动引入到飞行器模型中,得到不确定条件下的混沌多项式展开
其中,i=1,2,…,NY,NY是输出响应Yi中变量的个数,ε为随机气动摄动变量,ε=[ε1,…,εd],d为随机变量的维数,为第i个输出响应的多项式展开系数,而Ψj为第i个输出响应的混沌多项式;对于存在多维和混合分布的不确定项,如果它们是相互独立的,则概率密度分布为每个随机变量概率密度分布函数的乘积,如果不确定项存在两个相互关联的分量,采用纳塔夫方法进行转换,以随机变量的边缘累积分布函数和协方差矩阵为基本信息,将相关非正态随机变量直接转换为独立标准正态随机变量;
步骤22、采用截断误差作为指标,利用谱分解理论自适应更新混沌多项展开式的基函数,当截断误差超过阈值时,混沌多项展开式使用的基将自适应地更新,将变构型飞行器考虑不确定性的模型组转换为等价的增广模型形式,对随机变量进行采样得到有
其中,x为飞行状态,s=1,2,…,ns,ns为采样点个数,从而将不确定性建模问题转变为确定的增广建模问题;
步骤23、将变构型飞行段的多约束条件转换为增广模型的等价约束,采用能量状态法获得多约束条件下变构型飞行器的可行区域,将速度、高度、轨迹角和质量作为状态变量,将迎角和舵面作为控制变量,总能量看作是系统连续的状态变量,以飞行器可行区域边界对应的是飞行器能量状态变化率为0,通过计算能量状态变化率来确定飞行器可行区域;
步骤3、确定变构型飞行器优化性能指标,构建出变构型飞行器模型特性对闭环极限性能的约束关系,设计变构型飞行器指令的分层迭代策略,解算出不确定性条件下变构型飞行器的任务指令;具体包括如下步骤:
步骤31、根据最小期望/协方差,结合指标加权策略,确立多约束条件下变构型飞行器指令优化的性能指标;具体为:分析变构型飞行器模型耦合特性,综合稳定性、控制精度、鲁棒性要求,根据最小期望/协方差,结合指标加权策略,确立多约束条件下飞行器指令优化的性能指标,描述如下:
式(9)中,是优化的目标函数,/>由状态变量组成的向量,/>和/>给出了状态变量的约束范围,p是需要优化的指令,考虑到不确定性会存在于状态变量/>中或目标函数的输出/>中,采用鲁棒优化方法,构造关于原优化目标函数f的适应度函数F,即
式(10)中,Fn为基于适应度函数鲁棒性能指标,N是不确定性参数的数量,设计变量的约束不变,而和/>描述了不确定性的约束,/>是不确定性参数p的所有可能情况的集合,由pub和plb分别定义了上界和下界,/>是一个在/>中满足一定分布规律Λ(plbn,pubn)的随机参数,采用鲁棒正则化方法,构造如下的适应度函数:
式(11)中,δ是不确定性影响下指令关于标称位置的偏移量,而φ(δ)是飞行器指令关于标称位置偏移量的概率密度函数;
步骤32、根据飞行器不确定性与性能指标的关联分析,选取性能指标的权重系数,采用概率密度函数衡量不同类型不确定项对指标的影响程度,确定出的性能指标为
minJ=ωμJμ+ωσJσ (12)
其中,Jμ和Jσ为性能指标的均值和方差对应的权重值;
步骤33、通过分析变构型飞行器性能要求之间的耦合特性,构建指令的分层优化策略,优化问题表示如下:
minJ(P,K,r) (13)
其中,J表示目标函数,优化的目标是使得目标函数值取得最小,P,K,r为设计变量,即在优化设计过程中被用户控制的变量:P为构型参数,K为输入指令参数,r为状态参数,分层优化是将优化问题进行拆分若干个子系统,首先完成对一个子系统的优化,之后将最优的参数代入到下一步进行优化,获得不确定性条件下满足多约束条件的飞行器任务指令;
步骤4、采用数值仿真验证变构型飞行器所获得任务指令的有效性。
2.如权利要求1所述的不确定条件下的变构型飞行器指令鲁棒优化设计方法,其特征在于,步骤4中,采用数值仿真验证变构型飞行器所获得任务指令的有效性具体包括如下步骤:
步骤41、针对给定标称状态下,对任务指令进行分析,评估标称状态下任务优化的有效性,构建飞行器构型参数与闭环性能可达边界的映射关系,通过改变构型参数,分析闭环性能,验证标称状态下变构型设计的有效性;
步骤42、在强不确定的严苛条件下,对变构型飞行器模型参数进行拉偏,通过蒙脱卡洛仿真分析,依靠大量重复的仿真实验来,分析不同输入条件下,系统输出的变化范围,解算飞行性能是否满足期望的要求,从而验证控制性能增强设计结果的有效性。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211444860.5A CN115774900B (zh) | 2022-11-18 | 2022-11-18 | 一种不确定条件下的变构型飞行器指令鲁棒优化设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211444860.5A CN115774900B (zh) | 2022-11-18 | 2022-11-18 | 一种不确定条件下的变构型飞行器指令鲁棒优化设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115774900A CN115774900A (zh) | 2023-03-10 |
CN115774900B true CN115774900B (zh) | 2023-12-15 |
Family
ID=85389353
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211444860.5A Active CN115774900B (zh) | 2022-11-18 | 2022-11-18 | 一种不确定条件下的变构型飞行器指令鲁棒优化设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115774900B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117252130B (zh) * | 2023-11-17 | 2024-02-06 | 中国空气动力研究与发展中心计算空气动力研究所 | 飞行器近场压力分布预测方法、系统及相关设备 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107066640A (zh) * | 2016-11-22 | 2017-08-18 | 南京航空航天大学 | 基于保护映射的变构型飞行器模型迭代设计方法 |
CN107480335A (zh) * | 2017-07-12 | 2017-12-15 | 南京航空航天大学 | 一种高超声速飞行器模型迭代设计方法 |
CN107985626A (zh) * | 2017-10-27 | 2018-05-04 | 南京航空航天大学 | 一种基于可变构型空天飞行器的气动布局设计方法 |
CN109976154A (zh) * | 2019-03-04 | 2019-07-05 | 北京理工大学 | 一种基于混沌多项式和序列凸优化的飞行器轨迹优化方法 |
CN114115354A (zh) * | 2021-12-13 | 2022-03-01 | 北京航空航天大学 | 一种异质平台协同路径规划方法 |
CN114595510A (zh) * | 2022-01-27 | 2022-06-07 | 南京航空航天大学 | 一种空天飞行器融合结构外形参数的综合建模方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11132476B2 (en) * | 2017-12-13 | 2021-09-28 | Facebook, Inc. | Automatic aircraft design optimization based on joint aerodynamic, structural, and energy performance |
-
2022
- 2022-11-18 CN CN202211444860.5A patent/CN115774900B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107066640A (zh) * | 2016-11-22 | 2017-08-18 | 南京航空航天大学 | 基于保护映射的变构型飞行器模型迭代设计方法 |
CN107480335A (zh) * | 2017-07-12 | 2017-12-15 | 南京航空航天大学 | 一种高超声速飞行器模型迭代设计方法 |
CN107985626A (zh) * | 2017-10-27 | 2018-05-04 | 南京航空航天大学 | 一种基于可变构型空天飞行器的气动布局设计方法 |
CN109976154A (zh) * | 2019-03-04 | 2019-07-05 | 北京理工大学 | 一种基于混沌多项式和序列凸优化的飞行器轨迹优化方法 |
CN114115354A (zh) * | 2021-12-13 | 2022-03-01 | 北京航空航天大学 | 一种异质平台协同路径规划方法 |
CN114595510A (zh) * | 2022-01-27 | 2022-06-07 | 南京航空航天大学 | 一种空天飞行器融合结构外形参数的综合建模方法 |
Non-Patent Citations (3)
Title |
---|
Nonlinear Aircraft Closed Loop Control System Simulation Based on Dynamic Inversion;Jing Zhou等;《ACSE 2014》;全文 * |
Robust Multiple Model Predictive Control for Ascent Trajectory Tracking of Aerospace Vehicles;Rui Cao等;《IEEE Transactions on Aerospace and Electronic Systems》;第58卷(第2期);全文 * |
吸气式高超声速飞行器爬升段关键任务点的鲁棒优化;汤佳骏等;《宇航学报》;第41卷(第5期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN115774900A (zh) | 2023-03-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Houtzager et al. | Recursive predictor-based subspace identification with application to the real-time closed-loop tracking of flutter | |
Han et al. | Surrogate-based optimization | |
CN108121856B (zh) | 一种全飞行域飞行器动稳定性分析方法 | |
CN115774900B (zh) | 一种不确定条件下的变构型飞行器指令鲁棒优化设计方法 | |
CN113836634B (zh) | 一种大差异性气动数据的深度神经网络建模方法 | |
Cheung et al. | Local Lagrangian reduced-order modeling for the Rayleigh-Taylor instability by solution manifold decomposition | |
CN107065897A (zh) | 三自由度直升机显式模型预测控制方法 | |
White et al. | Neural networks predict fluid dynamics solutions from tiny datasets | |
CN111651940A (zh) | 基于pso-narx的涡轴发动机稳态模型辨识方法 | |
CN115495965A (zh) | 复杂航空结构在混合不确定性下时变可靠性的分析方法 | |
CN114722655B (zh) | 一种基于局部有限寿命疲劳约束条件的结构拓扑优化方法 | |
CN109446557B (zh) | 一种基于概率密度演化的随机气动弹性系统稳定性分析方法 | |
CN114595510A (zh) | 一种空天飞行器融合结构外形参数的综合建模方法 | |
CN108009320B (zh) | 一种高超声速飞行器面向控制的多系统关联建模方法 | |
Ge et al. | Genetic algorithm-based sensor allocation with nonlinear centralized fusion observable degree | |
CN116306246A (zh) | 一种大壁板铆接变形预测与优化方法 | |
CN115933381B (zh) | 一种多约束条件下的空天飞行器控制性能增强设计方法 | |
CN113158353B (zh) | 一种基于频域无参数辨识的倒立摆系统模型参数辨识方法 | |
Amrit et al. | Efficient multi-objective aerodynamic optimization by design space dimension reduction and co-kriging | |
Schirrer et al. | A comprehensive robust control design and optimization methodology for complex flexible-structure systems | |
Paw et al. | Uncertainty modeling, analysis and robust flight control design for a small UAV system | |
Campa et al. | Robust control and analysis of a towed underwater vehicle | |
Iannelli et al. | A balanced mode decomposition approach for equation-free reduced-order modeling of LPV aeroservoelastic systems | |
CN110990940B (zh) | 基于msvr的机翼装配定位布局设计方法 | |
Hmaiddouch et al. | Observer design for a class of discrete-time Takagi-Sugeno models with unmeasurable premise variables: Application to an asynchronous motor |
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 |