一种基于混合准则的机载分布式POS布局优化方法技术

本公开提供了一种机载分布式POS布局优化方法。首先，基于Fisher信息矩阵和模态应变能准则分别确定机翼结构待监控振型阶数范围和监控振型；然后，基于监控振型可观准则和模态置信度准则分别确定子IMU数量上限和下限；其次，确定布局方案编码并生成初始种群；之后，改进模态应变能准则，并将其与模态置信度准则和熵值法联合起来建立多目标函数模型；最后，运行遗传算法，确定分布式POS最优布局方案。该方法先确定子IMU安装数量的范围，避免了优化算法在全部布局方案中寻找最优解，减小了计算量，提高了优化算法的寻优效率；熵值法优化了改进模态应变能准则和模态置信度准则在多目标函数中的权重，提高了优化算法的寻优精度。

【技术实现步骤摘要】
一种基于混合准则的机载分布式POS布局优化方法
本公开涉及航空航天
，具体而言，涉及机载分布式位置姿态系统(POS)布局优化方法。可用于确定机载分布式POS中多个子系统的最优安装数量和位置，特别适合于机载阵列天线合成孔径雷达等多载荷同时对地观测需要确定位置姿态测量系统最优布局方式的情况。
技术介绍
集成多个或多种遥感载荷的机载分布式对地观测系统已经成为目前对地观测系统的主要发展趋势之一。这些多任务遥感载荷实现高精度成像的前提是获取各遥感载荷安置点的高精度运动参数。机载分布式位置姿态系统(PositionandOrientationSystem，POS)是多任务遥感载荷机载对地观测系统获取各遥感载荷安置点的位置、速度、姿态等运动参数的主要手段。机载分布式POS通常由一个或少量高精度主POS(也称为主系统)、多个中/低精度子惯性测量单元(InertialMeasurementUnit，IMU)(也称为子系统)、导航计算机系统(POSComputerSystem，PCS)和后处理软件四部分组成。高精度主POS通过传递对准将高精度的位置、速度、姿态等运动信息传递给各子IMU，从而实现各子IMU处运动信息的精确测量。其中，高精度主POS一般位于机舱内，提供时空基准信号，各子IMU安装在各遥感载荷安置点附近，为各遥感载荷提供精确的运动参数信息。对于机载多任务遥感载荷对地观测系统，分布式POS需要获取所有载荷安置点的高精度运动参数。由于载机对机载分布式POS的体积、重量均有严格的限制，特别是观测载荷分布在机翼下的情况，又因为成本等系统工程方面的考虑，在各观测载荷处均安装一个子IMU的方式既不合理也不现实。如何获取分布式POS子IMU的最佳安装数量和安装位置，即实现分布式POS的最优布局，是机载分布式POS的关键技术之一。分布式POS的优化布局问题属于工程中的最优化问题，需要用到各类的最优化算法。在运用优化算法进行最优布局时，首先需要确定候选布局方案组成的解空间。其中，候选布局方案包括不同安装数量下不同安装位置组成的所有方案。优化算法将在该解空间内寻找最优布局方案。目前，主要有两种方式确定候选布局方案的解空间：第一种方式是在全部候选布局方案的范围中寻找最优解，其中全部候选方案包括安装数量过少或过多的这些明显次优的布局方案，例如专利申请号为202010489022.4的专利采用全部候选布局方案作为解空间，导致求解计算量较大。例如有10个遥感载荷安置点，分布式POS子IMU可能的布局方案就多达个。即便优化算法具有快速寻优的特点，也没有必要损失时间对安装数量过少或过多的这些次优方案进行求解。第二种方式是先确定安装的确切数量，然后在这一确定的数量对应的不同安装位置的解空间内寻找最优解。该方式的优点是极大的缩小了解空间，提高了优化算法的寻优效率。但是确定安装数量是目前传感器布局优化的难点。现有确定传感器安装数量的方法各有优缺点，很难保证确定的数量是最优数量，这将导致在此基础上经过最优位置求解后最终确定的布局方案并非最优。例如，目前确定安装数量的方法主要有三类：经验法、基于结构损伤灵敏度分析或考虑模态频率重频的方法、基于模态识别领域的准则。第一类方法依赖经验，往往根据经验直接确定安装数量，这类方法缺乏理论指导，精度不高；第二类方法中，结构损伤灵敏度分析适于监测结构损伤而非分布式POS中测量点如何布设问题；此外，分布式POS中机翼结构也没有模态频率重频的情况，因此第二类方法对于分布式POS的布局来说不适用。第三类方法中，常用的有监控振型可观准则和模态置信度准则(ModalAssuranceCriterion，MAC)。第一种准则考虑了系统的可观性，根据系统可观性理论，要让包含n阶监控振型的系统可观，最少需要n个传感器。其中监控振型具体指包含系统模态振型主要阶次的振型，监控振型需要相关准则事先确定。(刘伟.基于有效独立的改进传感器优化布置方法研究[J].振动与冲击,2013,32(06):54-62.)。该准则的优点是用较少的传感器就能保证监控振型被辨识出来，但是辨识精度不高，导致后面使用多目标函数模型确定最优布局方案的精度也不高。另外，传感器最优安装数量应该既能保证较高的布局方案的精度，又能保证较低的布局成本。但是该准则确定的安装数量不能保证较高的布局方案的精度，因此并非最优的安装数量。而在MAC中，首先逐个增加传感器的安装数量，然后计算MAC矩阵最大的非对角元素。最大的非对角元素随着传感器数量的增加越小越好，这样各监控振型之间的相关性越小，监控振型的空间交角越大。当继续增加传感器安装数量而MAC矩阵最大的非对角元素不再减小时，此时的安装数量就是所求传感器的安装数量(YiT.H.,LiH.N..Optimalsensorplacementforstructuralhealthmonitoringbasedonmultipleoptimizationstrategies[J].TheStructuralDesignofTallandSpecialBuildings,2011,20(7):881-900.)。MAC能够保证获取的监控振型的空间交角大，辨识精度高，但是由于安装数量随辨识精度的提高而迅速增加，通过该准则得到的安装数量通常较多，导致布局成本较高，因此这样确定出的安装数量也并非最优。综上，在机载分布式POS最优布局时，如果能先确定出子IMU安装数量的范围，然后在该数量范围对应的那部分解空间内寻找最优解，将在适当缩小解空间、提高优化算法寻优效率的同时保证最优解的质量。在确定了候选布局方案组成的解空间之后，需要建立多目标函数模型，然后采用优化算法求解约束条件下该多目标函数的最优解，即可确定出最优布局方案。多目标函数模型建立主要基于优化准则。常用的优化准则包括最小均方根误差准则、MAC、模态应变能准则、Gramian矩阵准则等。上述优化准则各有优缺点。其中，最小均方根误差准则是在确定最优布局方案时，以载体上所有载荷安置点处运动参数的均方根误差为准则，优点是能直观反应布局方案优劣，但是计算量过大，计算过程耗时过长。MAC有利于避免监控振型间的空间交角过小而丢失重要模态，缺点是会选择一些模态应变能较低的位置来安装传感器，即将部分传感器布置在振动信号弱的位置上，也就是结构变形较小的位置。这样一来，对于机载分布式POS而言，意味着未安装子IMU的载荷处于变形较大的区域。由于载机变形复杂，无法精确建模，因此利用已安装的子IMU处的运动参数，采用最小二乘拟合等方法获取的未安装子IMU的载荷点运动参数的精度就不高，降低了机载分布式POS对所有载荷安置点运动参数获取的整体精度。而模态应变能准则和Gramian矩阵准则都是将结构振动较大处作为优先安装传感器的位置。其中Gramian矩阵准则仅考虑了模态位移对结构振动幅值的影响，忽略了模态应变能对结构振动幅值的影响。而模态应变能准则同时考虑了模态位移和模态应变能对结构振动幅值的影响，因此相比Gramian矩阵准则，模态应变能准则在原理上更具优势。但是模态应变能准则无法保证监控振型间的较大的空间交角，容本文档来自技高网...

【技术保护点】
1.一种基于混合准则的机载分布式POS布局优化方法，具体步骤为：/n1.1以Fisher信息矩阵迹的变化率为准则确定机翼结构待监控振型阶数的范围；/n1.2使用模态应变能准则，确定机翼结构的监控振型；/n1.3使用监控振型可观准则确定子IMU安装数量的下限，并使用MAC确定子IMU安装数量上限；/n1.4确定布局方案的编码并生成初始种群；/n1.5改进模态应变能准则，并将其与模态置信度准则和熵值法联合起来建立多目标函数模型；/n1.6运行遗传算法，求解多目标函数，将适应度值最大个体对应的布局方案作为机载分布式POS的最优布局方案。/n

【技术特征摘要】
1.一种基于混合准则的机载分布式POS布局优化方法，具体步骤为：
1.1以Fisher信息矩阵迹的变化率为准则确定机翼结构待监控振型阶数的范围；
1.2使用模态应变能准则，确定机翼结构的监控振型；
1.3使用监控振型可观准则确定子IMU安装数量的下限，并使用MAC确定子IMU安装数量上限；
1.4确定布局方案的编码并生成初始种群；
1.5改进模态应变能准则，并将其与模态置信度准则和熵值法联合起来建立多目标函数模型；
1.6运行遗传算法，求解多目标函数，将适应度值最大个体对应的布局方案作为机载分布式POS的最优布局方案。


2.根据权利要求1所述的基于混合准则的机载分布式POS布局优化方法，其特征在于：所述的步骤1.1中针对以Fisher信息矩阵迹的变化率为准则确定机翼结构待监控振型阶数的范围，具体步骤为：
2.1获取机翼前n阶模态振型矩阵Φn
从有限元模型提取m个载荷安置点对应位置的前n阶模态数据，得到的模态振型矩阵Φn的表达式为：



其中n为模态振型的阶数；m为载荷安置点的个数；φj为第j阶模态振型向量j＝1,2,…,n；φij为第i个载荷安置点的第j阶模态振型向量的模态位移，i＝1,2,…,m，j＝1,2,…,n；
2.2计算Fisher信息矩阵的迹
Fisher信息矩阵Q(n)可以表示为：



由于模态振型具有正交性，则上式可以写成：



于是Fisher信息矩阵Q(n)的迹tr(Q(n))可以表示为：



2.3计算Fisher信息矩阵迹的变化率ROC，确定待监控振型阶数的范围1～nc
Fisher信息矩阵Q(n)的迹在前n阶次和前n+1阶次上的变化率的表达式为：



ROC随着振型阶数n的增加而变化；当n达到某一数值nc，ROC值趋于0。可以确定待监控振型阶数的范围为1～nc。


3.根据权利要求2所述的基于混合准则的机载分布式POS布局优化方法，其特征在于：所述的步骤1.2中针对使用模态应变能准则确定机翼结构的监控振型，具体步骤为：
3.1获取机翼结构整体刚度矩阵K
提取m个载荷安置点的结构整体刚度矩阵K，其表达式为：



其中，矩阵中的元素kab为刚度系数，a＝1,2,…,m，b＝1,2,…,m，m为载荷安置点的个数；
3.2使用模态应变能准则选出监控振型
在划定了待监控振型阶数范围1～nc的基础上，根据模态应变能大小，选取应变能贡献较大的模态作为监控振型；模态应变能MSE是机翼结构整体刚度矩阵和模态振型共同作用的结果，机翼结构第j阶模态应变能的表达式为：



根据上式计算各阶模态的模态应变能，将最大模态应变能对应阶次的模态振型作为监控振型，监控振型的阶数为nd。


4.根据权利要求2所述的基于混合准则的机载分布式POS布局优化方法，其特征在于：所述的步骤1.3中针对使用监控振型可观准则确定子IMU安装数量的下限N1，并使用MAC确定子IMU安装数量上限N2，具体步骤为：
4.1确定子IMU安装数量的下限N1
子IMU安装数量的下限由监控振型阶次的数量nd决定；监控振型的可观性决定布置子IMU数量至少为监控振型阶次的数量，如果子IMU数量少于监控振型阶次的数量，系统不可观测到所有监控振型，从而导致丢失重要模态；因此子IMU安装数量的下限N1等于监控振型的阶数nd。
4.2使用MAC确定子IMU安装数量的上限N2
MAC矩阵的表达式为：



其中MACef为MAC矩阵的第e行、第f列元素，e＝1,2,…,nd，f＝1,2,…,nd，nd为监控振型阶数；
MAC矩阵的元素MACef的表达式为：



其中e＝1,2,…,nd，f＝1,2,…,nd，nd为监控振型阶数；
确定子IMU安装数量上限N2步骤为：
a)利用列主元QR分解法对模态振型矩阵Φn进行列主元QR分解，得到子IMU的N个初始安置点，其余m-N个安置点称为候选安置点，令P＝m-N；
b)求初始安置点的模态置信度矩阵MACN，并求MACN的最大非对角元素maxN，其中MACN表示安置点个数为N的MAC矩阵，maxN是MACN的最大非对角元素；
c)从P个候选安置点选择一个安置点添加到初始安置点集合中，得到包含N+1个初始安置点的集合，求其模态置信度矩阵MACN+1，并求其最大非对角元素maxN+1；
d)计算MACN+1和MACN的最大非对角元素之差F＝maxN+1-maxN；
e)从N+1个初始安置点中移除第c)步添加的候选安置点，从P个候选安置点重新选择一个添加到初始安置点集合中，并计算其最大非对角元素之差F；重复移除、添加候选安置点和计算F的步骤，得到F(l)，l＝1,2,…,P，并将|F(l)|中的最小值|F(l)|min对应的安置点加入到初始安置点中，初始安置点扩充到N+1个，候选安置点缩减到P-1个；...

【专利技术属性】
技术研发人员：宫晓琳，符倚伦，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：北京;11