[go: up one dir, main page]

CN114925441B - A method for airborne distributed PHM computational modeling - Google Patents

A method for airborne distributed PHM computational modeling Download PDF

Info

Publication number
CN114925441B
CN114925441B CN202210115309.XA CN202210115309A CN114925441B CN 114925441 B CN114925441 B CN 114925441B CN 202210115309 A CN202210115309 A CN 202210115309A CN 114925441 B CN114925441 B CN 114925441B
Authority
CN
China
Prior art keywords
algorithm
phm
computing
airborne
component
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
CN202210115309.XA
Other languages
Chinese (zh)
Other versions
CN114925441A (en
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.)
Air Force Engineering University of PLA
Original Assignee
Air Force Engineering University of PLA
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 Air Force Engineering University of PLA filed Critical Air Force Engineering University of PLA
Priority to CN202210115309.XA priority Critical patent/CN114925441B/en
Publication of CN114925441A publication Critical patent/CN114925441A/en
Application granted granted Critical
Publication of CN114925441B publication Critical patent/CN114925441B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/02Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06Power analysis or power optimisation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/12Timing analysis or timing optimisation
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D10/00Energy efficient computing, e.g. low power processors, power management or thermal management

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses an airborne distributed PHM calculation modeling method, which specifically comprises the following steps: modeling an airborne PHM operation platform; PHM calculates task allocation optimization conditions. According to the method, the computing resources are classified and graded, the performance of the computing resources is defined, and a model of the computing resources is constructed; classifying the algorithms by distinguishing online algorithms from offline algorithms and the accuracy and cost of the algorithms, and constructing a PHM algorithm model; the importance of the airborne equipment is divided, the monitoring requirements of the components are determined, and an operation requirement model of the airborne PHM system is constructed corresponding to the corresponding algorithm.

Description

一种机载分布式PHM计算建模方法A method for airborne distributed PHM computational modeling

技术领域Technical Field

本发明涉及面向故障预测与健康管理(Prognostic and Health Management,PHM)计算的建模方法,具体涉及一种机载分布式PHM计算建模方法。The present invention relates to a modeling method for fault prediction and health management (PHM) calculation, and in particular to an airborne distributed PHM calculation modeling method.

背景技术Background technique

随着航空技术的不断发展,飞机的复杂度不断提高,PHM技术成为先进飞 机必备的技术之一。PHM技术的实现依赖强大的数据处理能力,目前的PHM 算法大多基于通用处理器平台进行开发和运行,然而飞机上的能量和空间都有 限,基于通用处理器平台开发的算法在运行时会消耗大量的CPU资源和能量, 使得基于通用处理器的PHM计算在机载平台上的实现受到严重制约。随着 PHM技术的不断发展,机载PHM系统监测的参数越来越多,机载PHM计算 需要处理的数据的数量和类型快速增长,PHM算法也越来越复杂,因此分布式 嵌入式系统开始成为机载PHM系统进行计算的主要平台。随着分布式系统通 信技术的不断发展,机载设备间数据和任务的界限已经被打破,各个机载设备 之间可以实现计算资源和数据的共享,由此为分布式计算提供了可能。但是现 有的针对机载PHM的分布式嵌入式系统计算任务都固化在对应的计算资源内, 并且运行状态固定,导致计算资源的严重浪费。With the continuous development of aviation technology, the complexity of aircraft is constantly increasing, and PHM technology has become one of the essential technologies for advanced aircraft. The implementation of PHM technology relies on powerful data processing capabilities. Most of the current PHM algorithms are developed and run based on general-purpose processor platforms. However, the energy and space on the aircraft are limited. The algorithms developed based on the general-purpose processor platform will consume a lot of CPU resources and energy when running, which seriously restricts the implementation of PHM calculations based on general-purpose processors on airborne platforms. With the continuous development of PHM technology, the number of parameters monitored by airborne PHM systems is increasing, the number and type of data that need to be processed by airborne PHM calculations are growing rapidly, and PHM algorithms are becoming more and more complex. Therefore, distributed embedded systems have begun to become the main platform for airborne PHM systems to perform calculations. With the continuous development of distributed system communication technology, the boundaries of data and tasks between airborne devices have been broken, and computing resources and data can be shared between various airborne devices, which provides the possibility for distributed computing. However, the existing distributed embedded system computing tasks for airborne PHM are solidified in the corresponding computing resources, and the running state is fixed, resulting in a serious waste of computing resources.

针对机载PHM系统中分布式计算任务分配和调度的优化成为优化机载 PHM系统运行能力的当务之急。Optimizing the allocation and scheduling of distributed computing tasks in airborne PHM systems has become a top priority for optimizing the operational capabilities of airborne PHM systems.

发明内容Summary of the invention

针对现有技术存在的问题,本发明提供一种机载分布式PHM计算建模方 法,具体包括下列步骤:In view of the problems existing in the prior art, the present invention provides an airborne distributed PHM computational modeling method, which specifically includes the following steps:

步骤一、机载PHM运行平台建模Step 1: Modeling the airborne PHM operation platform

STEP1.机载分布式计算资源建模STEP 1. Modeling of airborne distributed computing resources

机载PHM系统由许多能够运行PHM算法的计算资源组成,假设某一机载 PHM系统共有M个可用于运行PHM算法的计算资源;The airborne PHM system consists of many computing resources that can run the PHM algorithm. Assume that a certain airborne PHM system has a total of M computing resources that can be used to run the PHM algorithm.

(1)计算资源分类CType (1) Computing Resource Classification C Type

这些计算资源是不同的嵌入式平台,具有不同的处理器架构,每个计算资 源是同构系统或异构系统;根据机载PHM系统的计算资源所应用的嵌入式平 台以及系统结构,对计算资源进行分类,假设全机PHM系统中的计算资源共 有n类,为每一类资源分配一个标号i=1,2,…n,如表1所示;These computing resources are different embedded platforms with different processor architectures. Each computing resource is a homogeneous system or a heterogeneous system. The computing resources are classified according to the embedded platform and system structure used by the computing resources of the airborne PHM system. Assuming that there are n types of computing resources in the whole PHM system, a label is assigned to each type of resource. i = 1, 2, ... n, as shown in Table 1;

表1机载PHM计算资源分类Table 1 Classification of airborne PHM computing resources

其中,指代单ARM架构的计算资源,/>指代双DSP架构的计算资 源,指代DSP+FPGA异构架构的计算资源,/>指代 ARM+DSP+FPGA异构架构的计算资源;目前嵌入式平台发展的趋势是异构平 台,在不同的功能需求下经常有新的嵌入式异构系统架构出现,无法穷举所有 的架构类型,此处下标的规律是按照单核架构、多核同构架构、简单异构架构、 复杂异构架构的类型排序,当有新的嵌入式系统架构应用到机载PHM系统中 时,在类代号中继续补充新的下标序号即可;in, Refers to computing resources of a single ARM architecture, /> Refers to the computing resources of the dual DSP architecture. Refers to the computing resources of DSP+FPGA heterogeneous architecture,/> Refers to the computing resources of ARM+DSP+FPGA heterogeneous architecture; the current trend of embedded platform development is heterogeneous platform. New embedded heterogeneous system architectures often appear under different functional requirements. It is impossible to list all architecture types. The subscript rules here are sorted according to the types of single-core architecture, multi-core homogeneous architecture, simple heterogeneous architecture, and complex heterogeneous architecture. When a new embedded system architecture is applied to the airborne PHM system, the new subscript number can be added to the class code;

(2)算力分级CAbility (2) Computing Power Classification C Ability

在对计算资源分类的基础上,针对每一类计算资源,按照其计算能力进行 分级,规定第i号计算资源的算力等级为i=1,2,…,M;假设全机 PHM系统中吞吐量最小的计算资源为第l号计算资源,则规定其算力等级 />并规定/>为单位计算能力;假设全机PHM系统中吞吐 量最大的计算资源为第k号计算资源,参照单位计算能力的定义,则规定其算 力等级/>则CAbilitymax即为全机PHM系统中计算资源的算 力等级的最大值;计算资源参照单位计算能力向下取整,得到全机机载PHM 系统所有PHM算法模型运行平台对应的算力等级/>i=1,2,…,M;On the basis of the classification of computing resources, each type of computing resource is graded according to its computing power, and the computing power level of the i-th computing resource is specified as i = 1, 2, ..., M; Assuming that the computing resource with the smallest throughput in the whole PHM system is the computing resource No. 1, then its computing power level is specified/> And stipulate/> is the unit computing power; assuming that the computing resource with the largest throughput in the whole PHM system is the kth computing resource, refer to the definition of unit computing power to define its computing power level/> Then C Abilitymax is the maximum value of the computing power level of the computing resources in the whole aircraft PHM system; the computing resources are rounded down with reference to the unit computing power to obtain the computing power level corresponding to all PHM algorithm model running platforms of the whole aircraft PHM system/> i = 1, 2, ..., M;

依据表1以及算力分级,机载PHM系统中的每一个计算资源都可以对应 i=1,2,…,n,j=1,2,…,M,由此得到表2;According to Table 1 and the computing power classification, each computing resource in the airborne PHM system can correspond to and i=1,2,…,n,j=1,2,…,M, thus we get Table 2;

表2机载PHM计算资源性能Table 2 Performance of airborne PHM computing resources

STEP2.PHM算法性能建模STEP 2. PHM algorithm performance modeling

对于某个机载部件,需要对PHM算法的运行需求进行建模;For a certain airborne component, the operational requirements of the PHM algorithm need to be modeled;

(1)确定实时性:PHM工程人员在进行PHM算法开发时,根据算法处 理数据的时效性以及算法是否需要使用历史数据,将PHM算法分为在线算法 和离线算法;处理的数据具有时效性且不需要使用历史数据的算法模型即为具 有实时性的在线算法;处理的数据反应的是部件一段时间内的状态参数,或者 需要使用历史数据的算法模型即为不具有实时性的离线算法;(1) Determine real-time performance: When developing PHM algorithms, PHM engineers divide PHM algorithms into online algorithms and offline algorithms based on the timeliness of the data processed by the algorithm and whether the algorithm needs to use historical data. Algorithm models that process time-sensitive data and do not need to use historical data are online algorithms with real-time performance. Algorithm models that process data that reflects the status parameters of components over a period of time or need to use historical data are offline algorithms that do not have real-time performance.

(2)确定计算资源需求:对于需要在机载PHM系统上运行的算法α, PHM工程人员在进行PHM算法开发时,需要针对算法的复杂度、逻辑结构、 数据长度以及数据缓存等需求,通过算法在各类计算资源、各个算力等级的计 算资源上进行测试,根据测试结果,确定各类资源上该算法在类资源上 运行时,最少需要/>α代表不同的算法,i=1,2,…,n,算力等级的计 算资源才能满足算法运行的需求,由此能够为每个算法可以得到一个运行需求, 即为表3;(2) Determine the computing resource requirements: For the algorithm α that needs to be run on the airborne PHM system, PHM engineers need to test the algorithm on various computing resources and computing resources of various computing power levels according to the algorithm complexity, logical structure, data length, and data cache requirements when developing the PHM algorithm. Based on the test results, determine the algorithm’s performance on various resources. When running on a class resource, at least/> α represents different algorithms, i = 1, 2, ..., n, and computing resources of a certain computing power level can meet the requirements of algorithm operation. Thus, an operation requirement can be obtained for each algorithm, which is shown in Table 3;

表3算法α运行需求Table 3 Algorithm α operation requirements

STEP3.部件监测需求建模STEP 3. Component monitoring requirements modeling

假设:飞机可执行m种飞行任务,飞机的机载PHM系统一共对x个机载部 件进行状态监测,机载PHM系统获取其健康状态参数以表征部件的健康状态; 首先构建参数以表征部件的监测需求;Assumption: An aircraft can perform m types of flight missions. The aircraft's onboard PHM system monitors the status of x onboard components. The onboard PHM system obtains its health status parameters to characterize the health status of the components. First, construct parameters to characterize the monitoring requirements of the components.

(1)部件重要度i=1,2,…,x:(1) Component Importance i = 1, 2, ..., x:

根据部件失效对飞行安全产生的危害度对部件进行重要度的划分,部件失 效对飞行安全危害越大则越大;The importance of components is divided according to the degree of harm caused by component failure to flight safety. The greater the harm caused by component failure to flight safety, the greater the importance of components. The bigger;

(2)任务重要度j=1,2,…,m:(2) Task importance j = 1, 2, ..., m:

飞机的每一飞行任务都有一个编号j,j=1,2,…,m,其中m为飞机可执行 飞行任务种类,根据部件失效对飞行任务完成度的危害对部件进行部件的任务 重要度进行划分,部件失效对某一飞行任务的完成度危害越大则越大, 执行不同的飞行任务时同一个部件的/>可能相同;Each flight mission of the aircraft has a number j, j = 1, 2, ..., m, where m is the type of flight mission that the aircraft can perform. The components are divided into task importance according to the harm of component failure to the completion of the flight mission. The greater the harm of component failure to the completion of a certain flight mission, the greater the The larger the size, the more likely it is that the same component will perform different flight missions. Probably the same;

(3)部件健康等级NHealth(t):(3) Component health level N Health (t):

t表示NHealth(t)是时间的函数,随着时间的推移,部件的性能逐渐退化,NHealth(t)改变;t represents N Health (t) as a function of time. As time goes by, the performance of the component gradually degrades and N Health (t) changes;

机载部件的失效是逐渐退化的过程,在部件退化的初期,其退化过程往往 是平稳的,基于当前对部件健康状态等级的评估对部件进行健康等级的划分, 部件失效的可能性越大则NHealth(t)越大;The failure of airborne components is a gradual degradation process. In the early stage of component degradation, the degradation process is often stable. The components are divided into health levels based on the current assessment of the component health status level. The greater the possibility of component failure, the greater N Health (t);

(4)监测需求NMonitor(t):(4) Monitoring requirement N Monitor (t):

所构建的监测需求NMonitor(t),表征某一个部件对于监测的需求迫切度;The constructed monitoring requirement N Monitor (t) represents the urgency of a component’s need for monitoring;

STEP4.设置PHM算法运行方式STEP 4. Set the PHM algorithm operation mode

对于某个部件,工程人员在其PHM算法的开发阶段,需要根据其退化模 型为其选择若干种算法,在不同的监测需求下对应不同的监测采样率和故障诊 断与寿命预测算法;For a certain component, engineers need to select several algorithms for it according to its degradation model during the development phase of its PHM algorithm, and use different monitoring sampling rates and fault diagnosis and life prediction algorithms according to different monitoring requirements;

算法、采样率、监测间隔的设置具有以下规律:The settings of the algorithm, sampling rate, and monitoring interval follow the following rules:

①监测需求较低时,所选择的算法具备低功耗、低算力需求的特点,算法 可以是离线算法或者是在线算法,传感器的采样率较低,监测间隔时间 长;算法的目标是监测并发现部件的早期故障特征;① When the monitoring demand is low, the selected algorithm has the characteristics of low power consumption and low computing power requirements. The algorithm can be an offline algorithm or an online algorithm. The sampling rate of the sensor is low and the monitoring interval is long. The goal of the algorithm is to monitor and discover the early fault characteristics of the components.

②随着监测需求的增大,所选择的算法的实时性、计算精度提高,以在线 算法为主,传感器的采样率提高,监测间隔时间缩短甚至为持续监测; 算法的目标是在部件出现早期故障特征之后,准确识别部件的加速退化 过程;② As the monitoring demand increases, the real-time performance and calculation accuracy of the selected algorithm are improved, with the online algorithm as the main one, the sampling rate of the sensor is increased, the monitoring interval is shortened or even continuous monitoring is adopted; The goal of the algorithm is to accurately identify the accelerated degradation process of the component after the early fault characteristics of the component appear;

③当监测需求继续增大,应当选择具备高精度和实时性特点的算法,传感 器采样率达到最高,监测方式为实时监测;算法的目的是在部件的加速 退化阶段准确预测部件的剩余使用寿命,在部件发生故障时能够立即对 故障进行准确的识别和定位;③ When the monitoring demand continues to increase, an algorithm with high precision and real-time characteristics should be selected, the sensor sampling rate should be the highest, and the monitoring method should be real-time monitoring; the purpose of the algorithm is to accurately predict the remaining service life of the component during the accelerated degradation stage of the component, and to be able to immediately and accurately identify and locate the fault when the component fails;

在所有x个机载部件中,最高的监测需求为PHM开发人员为 每个部件在每个确定的NMonitor(t)下找到一种PHM算法,算法(1,0)表示部件 1在NMonitor(t)=0时运行的算法,算法(1,1)表示部件1在NMonitor(t)=1时运 行的算法,以此类推,算法(x,N)表示部件x在NMonitor(t)=N时运行的算法;PHM开发人员为每个部件在每个确定的NMonitor(t)下都设定PHM算法运行时数 据的采样率和算法两次连续运行之间的时间间隔;得到表4;Among all x airborne components, the highest monitoring requirement is The PHM developer finds a PHM algorithm for each component under each determined N Monitor (t). Algorithm (1, 0) represents the algorithm that component 1 runs when N Monitor (t) = 0, and algorithm (1, 1) represents the algorithm that component 1 runs when N Monitor (t) = 1. Similarly, algorithm (x, N) represents the algorithm that component x runs when N Monitor (t) = N. The PHM developer sets the sampling rate of the PHM algorithm runtime data and the time interval between two consecutive runs of the algorithm for each component under each determined N Monitor (t). Table 4 is obtained.

表4PHM算法运行方式Table 4 PHM algorithm operation mode

STEP5.确定算法运行性能STEP 5. Determine the algorithm performance

根据表4中的采样率和监测间隔,PHM工程人员在进行PHM算法开发时, 需要对PHM算法在不同的计算资源下的性能进行测试,测试内容为:对于算 法α,当其运行在类,计算能力为/>的计算资源下时,可以确定其 单次执行的计算延迟/>和功耗/>得到表5;According to the sampling rate and monitoring interval in Table 4, PHM engineers need to test the performance of the PHM algorithm under different computing resources when developing the PHM algorithm. The test content is: for algorithm α, when it runs on Class, computing power is/> When the computing resources are sufficient, the computing delay of a single execution can be determined./> and power consumption/> Table 5 is obtained;

表5 PHM算法α运行性能Table 5 PHM algorithm α running performance

STEP6.确定PHM系统运行需求STEP 6. Determine the PHM system operation requirements

当机载PHM系统运行时,能够根据每个机载部件的当前时刻的监测需求 确定当前时刻整个PHM系统中需要运行的所有算法及其采样率、运行间隔和 算法运行性能,且假设每个PHM算法都只在一个计算资源上运行;假设t时刻 一共有l个需要运行的PHM算法,则这些需要运行的PHM算法及其采样率、 运行间隔和算法运行性能构成一个集合;就确定出t时刻机载PHM系统的运行 需求和运行性能,即机载PHM系统运行的需求及任务分配的空间;When the airborne PHM system is running, it can determine all algorithms that need to be run in the entire PHM system at the current moment and their sampling rates, running intervals and algorithm running performance according to the monitoring requirements of each airborne component at the current moment, and assume that each PHM algorithm is only run on one computing resource; assume that there are a total of l PHM algorithms that need to be run at time t, then these PHM algorithms that need to be run and their sampling rates, running intervals and algorithm running performance constitute a set; the running requirements and running performance of the airborne PHM system at time t are determined, that is, the running requirements of the airborne PHM system and the space for task allocation;

步骤二、PHM计算任务分配优化空间构造Step 2: PHM calculation task allocation and optimization of spatial structure

假设某一时刻,机载PHM系统需要同时运行K个算法;Assume that at a certain moment, the airborne PHM system needs to run K algorithms simultaneously;

在获得PHM系统运行的需求及任务分配的空间的基础上,在任务分配空 间内对PHM计算任务分配优化;Based on the requirements of PHM system operation and the space of task allocation, the PHM computing task allocation is optimized within the task allocation space;

以机载PHM系统中所有算法执行的延迟和功耗作为评价指标构造性能指 标函数J(D):The performance index function J(D) is constructed by taking the delay and power consumption of all algorithms executed in the airborne PHM system as evaluation indicators:

其中,in,

为K维列向量,表示一共有K个需要运行的算法,分别分配到 n1、n2、…、nK号计算资源上运行; is a K-dimensional column vector, indicating that there are K algorithms that need to be run, which are respectively assigned to n 1 , n 2 , …, n K computing resources for execution;

和/>分别表示第i个算法在当前分配的计算资源中执行的延迟时间 和功耗; and/> They represent the delay time and power consumption of the i-th algorithm executed in the currently allocated computing resources respectively;

k1、k2为权重因子,表示算法执行延迟和计算功耗在性能指标中的权重, 满足k1+k2=1,k1、k2≥0;k 1 and k 2 are weight factors, indicating the weights of algorithm execution delay and computing power consumption in performance indicators, satisfying k 1 +k 2 =1, k 1 , k 2 ≥0;

约束条件为:算法需要在能够在满足最低算力等级需求的计算资源上运行;The constraints are: the algorithm needs to run on computing resources that can meet the minimum computing power level requirements;

在上述分配优化条件的基础上,再利用寻优算法,将每个需要运行的 PHM算法分配到一个计算资源上,使得上述分配优化条件中的性能指标函数 J(D)最小,即达到PHM系统算法运行的总体性能最优。On the basis of the above allocation optimization conditions, the optimization algorithm is used to allocate each PHM algorithm that needs to be run to a computing resource, so that the performance indicator function J(D) in the above allocation optimization conditions is minimized, that is, the overall performance of the PHM system algorithm operation is optimal.

本发明通过将计算资源进行分类、分级,明确计算资源的性能,构建计算 资源的模型。通过区分在线、离线算法以及算法的精度、代价对算法进行分类, 构建PHM算法模型。通过对机载设备的重要度进行划分,确定部件的监测需 求,并对应相应的算法,构建机载PHM系统的运行需求模型。The present invention classifies and grades computing resources, clarifies the performance of computing resources, and constructs a computing resource model. The algorithm is classified by distinguishing online and offline algorithms and the accuracy and cost of the algorithm to construct a PHM algorithm model. The monitoring requirements of components are determined by dividing the importance of airborne equipment, and the corresponding algorithms are used to construct an operation requirement model of the airborne PHM system.

具体实施方式Detailed ways

本发明的机载分布式PHM计算建模方法具体包括下列步骤:The airborne distributed PHM computational modeling method of the present invention specifically includes the following steps:

步骤一、机载PHM运行平台建模Step 1: Modeling the airborne PHM operation platform

首先对针对飞机PHM系统的计算需求与计算能力进行建模;First, the computational requirements and computing power for the aircraft PHM system are modeled;

STEP1.机载分布式计算资源建模STEP 1. Modeling of airborne distributed computing resources

机载PHM系统由许多能够运行PHM算法的计算资源组成,假设某一机载 PHM系统共有M个可用于运行PHM算法的计算资源;The airborne PHM system consists of many computing resources that can run the PHM algorithm. Assume that a certain airborne PHM system has a total of M computing resources that can be used to run the PHM algorithm.

(1)计算资源分类CType (1) Computing Resource Classification C Type

这些计算资源是不同的嵌入式平台,具有不同的处理器架构,每个计算资 源是同构系统或异构系统;根据机载PHM系统的计算资源所应用的嵌入式平 台以及系统结构,对计算资源进行分类,假设全机PHM系统中的计算资源共 有n类,为每一类资源分配一个标号如表1所示;These computing resources are different embedded platforms with different processor architectures. Each computing resource is a homogeneous system or a heterogeneous system. The computing resources are classified according to the embedded platform and system structure used by the computing resources of the airborne PHM system. Assuming that there are n types of computing resources in the whole PHM system, a label is assigned to each type of resource. As shown in Table 1;

表1机载PHM计算资源分类Table 1 Classification of airborne PHM computing resources

其中,指代单ARM架构的计算资源,/>指代双DSP架构的计算资 源,指代DSP+FPGA异构架构的计算资源,/>指代 ARM+DSP+FPGA异构架构的计算资源……目前嵌入式平台发展的趋势是异构 平台,在不同的功能需求下经常有新的嵌入式异构系统架构出现,无法穷举所 有的架构类型,此处下标的规律是按照单核架构、多核同构架构、简单异构架 构、复杂异构架构的类型排序,当有新的嵌入式系统架构应用到机载PHM系 统中时,在类代号中继续补充新的下标序号即可。in, Refers to computing resources of a single ARM architecture, /> Refers to the computing resources of the dual DSP architecture. Refers to the computing resources of DSP+FPGA heterogeneous architecture,/> Refers to the computing resources of ARM+DSP+FPGA heterogeneous architecture... The current development trend of embedded platforms is heterogeneous platforms. New embedded heterogeneous system architectures often appear under different functional requirements. It is impossible to list all architecture types. The subscript rule here is to sort by single-core architecture, multi-core homogeneous architecture, simple heterogeneous architecture, and complex heterogeneous architecture. When a new embedded system architecture is applied to the airborne PHM system, just continue to add new subscript numbers in the class code.

(2)算力分级CAbility (2) Computing Power Classification C Ability

在对计算资源分类的基础上,针对每一类计算资源,按照其计算能力进行 分级,规定第i号计算资源的算力等级为假设全机 PHM系统中吞吐量最小的计算资源为第l号计算资源(例如某ARM),则规定 其算力等级/>并规定为单位计算能力。假设全机PHM 系统中吞吐量最大的计算资源为第k号计算资源(例如某FPGA),参照单位计 算能力的定义,则规定其算力等级/>则CAbilitymax即为全 机PHM系统中计算资源的算力等级的最大值。例如:假设上述l号计算资源 (某ARM)的吞吐量为2Mb/s,则/>假设上述k号计算资源(吞 吐量最大的计算资源,某FPGA)的吞吐量为1234.5Mb/s,则规定其算力等级 为/>(1234.5/2=617.25向下取整)为全机PHM系 统中计算资源的算力等级的最大值;假设a号计算资源的吞吐量为666.6Mb/s, 则其算力等级/>(666.6/2=333.3向下取整)。其余的计算资源参 照单位计算能力向下取整,得到全机机载PHM系统所有PHM算法模型运行平 台对应的算力等级/> On the basis of the classification of computing resources, each type of computing resource is graded according to its computing power, and the computing power level of the i-th computing resource is specified as Assuming that the computing resource with the smallest throughput in the whole PHM system is the No. 1 computing resource (for example, an ARM), its computing power level is specified/> And stipulate Assume that the computing resource with the largest throughput in the whole PHM system is the kth computing resource (e.g., a certain FPGA). Referring to the definition of unit computing power, its computing power level is specified. Then C Abilitymax is the maximum value of the computing power level of computing resources in the whole PHM system. For example, assuming that the throughput of the above computing resource No. 1 (a certain ARM) is 2Mb/s, then/> Assuming that the throughput of the computing resource No. k (the computing resource with the highest throughput, a certain FPGA) is 1234.5Mb/s, its computing power level is defined as/> (1234.5/2=617.25 rounded down) is the maximum value of the computing power level of computing resources in the whole PHM system; assuming that the throughput of computing resource a is 666.6Mb/s, its computing power level/> (666.6/2=333.3 rounded down). The remaining computing resources are rounded down with reference to the unit computing power to obtain the computing power level corresponding to all PHM algorithm model operation platforms of the entire aircraft's onboard PHM system/>

依据表1以及算力分级,机载PHM系统中的每一个计算资源都可以对应和/>由此得到表2。(例如,1号资 源为单ARM架构,为/>类计算资源,算力等级为1;2号资源为双DSP同 构架构,为/>类计算资源,算力等级为3;M号资源为ARM+DSP+FPGA 异构架构,为/>类计算资源,算力等级为/>)。According to Table 1 and the computing power classification, each computing resource in the airborne PHM system can correspond to and/> Thus, Table 2 is obtained. (For example, resource No. 1 is a single ARM architecture, which is/> Class computing resources, computing power level 1; resource No. 2 is a dual DSP homogeneous architecture, for/> Class computing resources, computing power level 3; M resources are ARM+DSP+FPGA heterogeneous architecture, for/> Class computing resources, computing power level is/> ).

表2机载PHM计算资源性能Table 2 Performance of airborne PHM computing resources

STEP2.PHM算法性能建模STEP 2. PHM algorithm performance modeling

对于某个机载部件,需要对PHM算法的运行需求进行建模;For a certain airborne component, the operational requirements of the PHM algorithm need to be modeled;

(1)确定实时性:PHM工程人员在进行PHM算法开发时,根据算法处 理数据的时效性以及算法是否需要使用历史数据,将PHM算法分为在线算法 和离线算法。处理的数据具有时效性且不需要使用历史数据的算法模型即为具 有实时性的在线算法;处理的数据反应的是部件一段时间内的状态参数,或者 需要使用历史数据的算法模型即为不具有实时性的离线算法。(1) Determine real-time performance: When developing PHM algorithms, PHM engineers divide PHM algorithms into online algorithms and offline algorithms based on the timeliness of the data processed by the algorithm and whether the algorithm needs to use historical data. An online algorithm is one that processes data that is time-sensitive and does not require the use of historical data; an offline algorithm is one that processes data that reflects the state parameters of a component over a period of time or requires the use of historical data.

(2)确定计算资源需求:对于需要在机载PHM系统上运行的算法α(例 如,α代表滑动平均滤波算法),PHM工程人员在进行PHM算法开发时,需要 针对算法的复杂度、逻辑结构、数据长度以及数据缓存等需求,通过算法在各 类计算资源、各个算力等级的计算资源上进行测试,根据测试结果,确定各类 资源上该算法在类资源上运行时,最少需要(α代表不同的算法,i=1,2,…,n)算力等级的计算资源才能满足算法运 行的需求,由此能够为每个算法可以得到一个运行需求,即为表3。(2) Determine the computing resource requirements: For the algorithm α (for example, α represents the sliding average filter algorithm) that needs to be run on the airborne PHM system, PHM engineers need to test the algorithm on various computing resources and computing resources of various computing power levels according to the algorithm complexity, logical structure, data length, and data cache requirements when developing the PHM algorithm. Based on the test results, determine the algorithm's performance on various resources. When running on a class resource, at least (α represents different algorithms, i=1, 2, ..., n) Computing resources of a certain computing power level can meet the requirements of algorithm operation, so an operation requirement can be obtained for each algorithm, which is shown in Table 3.

表3算法α运行需求Table 3 Algorithm α operation requirements

例如,表3表示算法α需要在类、算力等级大于/>的计算资源 上;或者在/>类、算力等级大于/>的计算资源上;或者在在/>类、 算力等级大于的计算资源上…;或者在在/>类、算力等级大于 />的计算资源上能够运行。For example, Table 3 shows that algorithm α needs to Class, computing power level greater than/> computing resources; or in/> Class, computing power level greater than/> on computing resources; or on/> Class, computing power level greater than on computing resources of ...; or on /> Class, computing power level greater than/> can run on computing resources.

STEP3.部件监测需求建模STEP 3. Component monitoring requirements modeling

假设:飞机可执行m种飞行任务,飞机的机载PHM系统一共对x个机载部 件进行状态监测,机载PHM系统获取其健康状态参数以表征部件的健康状态; 首先构建参数以表征部件的监测需求;Assumption: An aircraft can perform m types of flight missions. The aircraft's onboard PHM system monitors the status of x onboard components. The onboard PHM system obtains its health status parameters to characterize the health status of the components. First, construct parameters to characterize the monitoring requirements of the components.

(1)部件重要度 (1) Component Importance

根据部件失效对飞行安全产生的危害度对部件进行重要度的划分,部件失 效对飞行安全危害越大则越大;The importance of components is divided according to the degree of harm caused by component failure to flight safety. The greater the harm caused by component failure to flight safety, the greater the importance of components. The bigger;

(2)任务重要度 (2) Task importance

飞机的每一飞行任务都有一个编号j(j=1,2,…,m),其中m为飞机可执行 飞行任务种类,根据部件失效对飞行任务完成度的危害对部件进行部件的任务 重要度进行划分,部件失效对某一飞行任务的完成度危害越大则越大, 执行不同的飞行任务时同一个部件的/>可能相同;Each flight mission of the aircraft has a number j (j = 1, 2, ..., m), where m is the type of flight mission that the aircraft can perform. The components are divided into task importance according to the harm of component failure to the completion of the flight mission. The greater the harm of component failure to the completion of a certain flight mission, the greater the The larger the size, the more likely it is that the same component will perform different flight missions. Probably the same;

(3)部件健康等级NHealth(t):(3) Component health level N Health (t):

t表示NHealth(t)是时间的函数,随着时间的推移,部件的性能逐渐退化, NHealth(t)改变。t represents N Health (t) as a function of time. As time goes by, the performance of the component gradually degrades and N Health (t) changes.

机载部件的失效是逐渐退化的过程,在部件退化的初期,其退化过程往往 是平稳的,基于当前对部件健康状态等级的评估对部件进行健康等级的划分, 部件失效的可能性越大则NHealth(t)越大;The failure of airborne components is a gradual degradation process. In the early stage of component degradation, the degradation process is often stable. The components are divided into health levels based on the current assessment of the component health status level. The greater the possibility of component failure, the greater N Health (t);

(4)监测需求NMonitor(t):(4) Monitoring requirement N Monitor (t):

所构建的监测需求NMonitor(t),表征某一个部件对于监测的需求迫切度;The constructed monitoring requirement N Monitor (t) represents the urgency of a component’s need for monitoring;

STEP4.设置PHM算法运行方式STEP 4. Set the PHM algorithm operation mode

对于某个部件,工程人员在其PHM算法的开发阶段,需要根据其退化模 型为其选择若干种算法,在不同的监测需求下对应不同的监测采样率和故障诊 断与寿命预测算法;For a certain component, engineers need to select several algorithms for it according to its degradation model during the development phase of its PHM algorithm, and use different monitoring sampling rates and fault diagnosis and life prediction algorithms according to different monitoring requirements;

算法、采样率、监测间隔的设置具有以下规律:The settings of the algorithm, sampling rate, and monitoring interval follow the following rules:

④监测需求较低时,所选择的算法具备低功耗、低算力需求的特点,算法 可以是离线算法或者是在线算法,传感器的采样率较低,监测间隔时间 长;算法的目标是监测并发现部件的早期故障特征;④ When the monitoring demand is low, the selected algorithm has the characteristics of low power consumption and low computing power requirements. The algorithm can be an offline algorithm or an online algorithm. The sampling rate of the sensor is low and the monitoring interval is long. The goal of the algorithm is to monitor and discover the early fault characteristics of the components.

⑤随着监测需求的增大,所选择的算法的实时性、计算精度提高,以在线 算法为主,传感器的采样率提高,监测间隔时间缩短甚至为持续监测; 算法的目标是在部件出现早期故障特征之后,准确识别部件的加速退化 过程;⑤ With the increase of monitoring demand, the real-time performance and calculation accuracy of the selected algorithm are improved, with the online algorithm as the main one, the sampling rate of the sensor is increased, the monitoring interval is shortened or even continuous monitoring; The goal of the algorithm is to accurately identify the accelerated degradation process of the component after the early fault characteristics of the component appear;

⑥当监测需求继续增大,应当选择具备高精度和实时性特点的算法,传感 器采样率达到最高,监测方式为实时监测;算法的目的是在部件的加速 退化阶段准确预测部件的剩余使用寿命,在部件发生故障时能够立即对 故障进行准确的识别和定位;⑥ When the monitoring demand continues to increase, an algorithm with high precision and real-time characteristics should be selected, the sensor sampling rate should be the highest, and the monitoring method should be real-time monitoring; the purpose of the algorithm is to accurately predict the remaining service life of the component during the accelerated degradation stage of the component, and to be able to immediately and accurately identify and locate the fault when the component fails;

在所有x个机载部件中,最高的监测需求为PHM开发人员为 每个部件在每个确定的NMonitor(t)下找到一种PHM算法,算法(1,0)表示部件1在NMonitor(t)=0时运行的算法(/>三个指 标中,仅NMissionj的取值可以为0,当/>为0时NMonitor(t)取0),算法 (1,1)表示部件1在NMonitor(t)=1时运行的算法,……,算法(x,N)表示部件x 在NMonitor(t)=N时运行的算法(NMonitor(t)为大于等于0的整数,因为 />三个指标都是大于等于0的整数,在前面指 标的定义中已经明确);PHM开发人员为每个部件在每个确定的NMonitor(t)下 都设定PHM算法运行时数据的采样率和算法两次连续运行之间的时间间隔;采 样率(1,0)表示部件1在NMonitor(t)=0时设置的数据采样率,监测间隔(1,0)表 示部件1在NMonitor(t)=0时算法连续两次运行之间的时间间隔,采样率(1,1) 表示部件1在NMonitor(t)=1时设置的数据采样率,监测间隔(1,1)表示部件1 在NMonitor(t)=1时算法连续两次运行之间的时间间隔,……,采样率(x,N)表 示部件x在NMonitor(t)=N时设置的数据采样率,监测间隔(x,N)表示部件x在 NMonitor(t)=N时算法连续两次运行之间的时间间隔;由此可以得到表4;Among all x airborne components, the highest monitoring requirement is The PHM developer finds a PHM algorithm for each component under each determined N Monitor (t). Algorithm (1, 0) represents the algorithm that component 1 runs when N Monitor (t) = 0 (/> Among the three indicators, only NMissionj can be 0. When N Monitor (t) is 0, N Monitor (t) is 0). Algorithm (1, 1) represents the algorithm that component 1 runs when N Monitor (t) = 1. ..., Algorithm (x, N) represents the algorithm that component x runs when N Monitor (t) = N (N Monitor (t) is an integer greater than or equal to 0, because/> The three indicators are all integers greater than or equal to 0, which have been clearly defined in the previous indicator definitions); PHM developers set the data sampling rate and the time interval between two consecutive runs of the PHM algorithm for each component under each determined N Monitor (t); sampling rate (1, 0) represents the data sampling rate set by component 1 when N Monitor (t) = 0, monitoring interval (1, 0) represents the time interval between two consecutive runs of the algorithm by component 1 when N Monitor (t) = 0, sampling rate (1, 1) represents the data sampling rate set by component 1 when N Monitor (t) = 1, monitoring interval (1, 1) represents the time interval between two consecutive runs of the algorithm by component 1 when N Monitor (t) = 1, ..., sampling rate (x, N) represents the data sampling rate set by component x when N Monitor (t) = N, monitoring interval (x, N) represents the time interval between two consecutive runs of the algorithm by component x when N Monitor (t) = N; thus, Table 4 can be obtained;

表4PHM算法运行方式Table 4 PHM algorithm operation mode

STEPS.确定算法运行性能STEPS. Determine the algorithm performance

根据表4中的采样率和监测间隔,PHM工程人员在进行PHM算法开发时, 需要对PHM算法在不同的计算资源下的性能进行测试,测试内容为:对于算 法α,当其运行在类,计算能力为/>的计算资源下时,可以确定其 单次执行的计算延迟/>和功耗/>假设对于算法α而言:算法α在 />类资源下最低计算能力需求为3,在类资源下最低计算能力需求为 2,在/>类资源下最低计算能力需求为1;在类、算力为3的计算资 源下单次执行的计算延迟为/>功耗为/>类、算力为 2的计算资源下单次执行的计算延迟为/>功耗为/>类、 算力为3的计算资源下单次执行的计算延迟为/>功耗为/>类、算力为1的计算资源下单次执行的计算延迟为/>功耗为 />类、算力为2的计算资源下单次执行的计算延迟为 />功耗为/>在/>类、算力为3的计算资源下单次执行的 计算延迟为/>功耗为/>规定算法无法运行的计算资源的 />由此得到表5。According to the sampling rate and monitoring interval in Table 4, PHM engineers need to test the performance of the PHM algorithm under different computing resources when developing the PHM algorithm. The test content is: for algorithm α, when it runs on Class, computing power is/> When the computing resources are sufficient, the computing delay of a single execution can be determined./> and power consumption/> Assume that for algorithm α: Algorithm α in/> The minimum computing power requirement for this resource type is 3. The minimum computing power requirement for class resources is 2. The minimum computing power requirement for this resource type is 1; The computational delay of a single execution under the computing resource of class and computing power 3 is/> Power consumption is/> exist The computational delay of a single execution under the computing resource of class and computing power 2 is/> Power consumption is/> exist The computing delay of a single execution under the computing resource of class and computing power 3 is/> Power consumption is/> exist The computational delay of a single execution under a computing resource of class 1 is/> Power consumption is/> exist The computational delay of a single execution under the computing resource class with a computing power of 2 is/> Power consumption is/> In/> The computational delay of a single execution under the computing resource of class and computing power 3 is/> Power consumption is/> The specified algorithm cannot run on computing resources/> This resulted in Table 5.

表5 PHM算法α运行性能Table 5 PHM algorithm α running performance

同理,根据工程人员进行的性能测试,每个算法都能够得到一个表征该算 法运行性能的表;Similarly, based on the performance tests conducted by engineers, each algorithm can obtain a table that characterizes the algorithm's operating performance;

STEP6.确定PHM系统运行需求STEP 6. Determine the PHM system operation requirements

当机载PHM系统运行时,能够根据每个机载部件的当前时刻的监测需求 确定当前时刻整个PHM系统中需要运行的所有算法及其采样率、运行间隔和 算法运行性能,且假设每个PHM算法都只在一个计算资源上运行;假设t时刻 一共有l个需要运行的PHM算法,则这些需要运行的PHM算法及其采样率、 运行间隔和算法运行性能构成一个集合。就确定出t时刻机载PHM系统的运行 需求和运行性能,即机载PHM系统运行的需求及任务分配的空间;When the airborne PHM system is running, it can determine all the algorithms that need to be run in the entire PHM system at the current moment and their sampling rates, running intervals, and algorithm running performance according to the monitoring requirements of each airborne component at the current moment, and assume that each PHM algorithm is only run on one computing resource; assume that there are l PHM algorithms that need to be run at time t, then these PHM algorithms that need to be run and their sampling rates, running intervals, and algorithm running performance constitute a set. The running requirements and running performance of the airborne PHM system at time t are determined, that is, the running requirements of the airborne PHM system and the space for task allocation;

步骤二、PHM计算任务分配优化空间构造Step 2: PHM calculation task allocation and optimization of spatial structure

假设某一时刻,机载PHM系统需要同时运行K个算法。Assume that at a certain moment, the airborne PHM system needs to run K algorithms simultaneously.

在获得PHM系统运行的需求及任务分配的空间的基础上(步骤一STEP6 中确定的,t时刻需要运行的PHM算法及其采样率、运行间隔和算法运行性 能),在任务分配空间内对PHM计算任务分配优化;Based on the requirements of PHM system operation and the space of task allocation (the PHM algorithm to be run at time t and its sampling rate, operation interval and algorithm operation performance determined in step 1, STEP 6), optimize the PHM calculation task allocation within the task allocation space;

PHM算法执行的延迟时间和计算产生的功耗是评价整个PHM系统运行平 台性能的重要指标,而PHM计算任务分配优化的目的就是使所有计算资源上 运行的计算任务总体的延迟最低且功耗最小,从而以低功耗提高对系统监测数 据处理的实时性;The delay time of PHM algorithm execution and the power consumption generated by calculation are important indicators for evaluating the performance of the entire PHM system operation platform. The purpose of PHM calculation task allocation optimization is to minimize the overall delay and power consumption of the calculation tasks running on all computing resources, thereby improving the real-time performance of system monitoring data processing with low power consumption.

以机载PHM系统中所有算法执行的延迟和功耗作为评价指标构造性能指 标函数J(D):The performance index function J(D) is constructed by taking the delay and power consumption of all algorithms executed in the airborne PHM system as evaluation indicators:

其中,in,

为K维列向量,表示一共有K个需要运行的算法,分别分配到 n1、n2、…、nK号计算资源上运行; is a K-dimensional column vector, indicating that there are K algorithms that need to be run, which are respectively assigned to n 1 , n 2 , …, n K computing resources for execution;

和/>分别表示第i个算法在当前分配的计算资源中执行的延迟时间 和功耗; and/> They represent the delay time and power consumption of the i-th algorithm executed in the currently allocated computing resources respectively;

k1、k2为权重因子,表示算法执行延迟和计算功耗在性能指标中的权重, 满足k1+k2=1,k1、k2≥0;k 1 and k 2 are weight factors, indicating the weights of algorithm execution delay and computing power consumption in performance indicators, satisfying k 1 +k 2 =1, k 1 , k 2 ≥0;

约束条件为:算法需要在能够在满足最低算力等级需求的计算资源上运行;The constraints are: the algorithm needs to run on computing resources that can meet the minimum computing power level requirements;

在上述分配优化条件的基础上,再利用各种寻优算法,如群粒子群算法 (GarcíagonzaloE,Fernándezmartínez JL.A brief historical review of particle swarmoptimization(PSO)[J].Journal of Bioinformatics&Intelligent Control,2012,1(1):3-16.), 灰狼算法(MIRJALILI S,MIRJALILI S M,LEWIS A.Grey wolf optimization[J].Advances in Engineering Software,2014,69(7):46-61.),将每个需要 运行的PHM算法分配到一个计算资源上,使得上述分配优化条件中的性能指 标函数J(D)最小,即达到PHM系统算法运行的总体性能最优(性能指标函数 的构建参考最优控制理论中的性能指标函数)。On the basis of the above allocation optimization conditions, various optimization algorithms are used, such as particle swarm optimization (Garcíagonzalo E, Fernándezmartínez JL. A brief historical review of particle swarm optimization (PSO) [J]. Journal of Bioinformatics & Intelligent Control, 2012, 1 (1): 3-16.), grey wolf algorithm (MIRJALILI S, MIRJALILI S M, LEWIS A. Grey wolf optimization [J]. Advances in Engineering Software, 2014, 69 (7): 46-61.), to allocate each PHM algorithm to be run to a computing resource, so that the performance index function J (D) in the above allocation optimization conditions is minimized, that is, the overall performance of the PHM system algorithm operation is optimal (the construction of the performance index function refers to the performance index function in the optimal control theory).

Claims (1)

1.一种机载分布式PHM计算建模方法,其特征在于,具体包括下列步骤:1. An airborne distributed PHM computational modeling method, characterized in that it specifically includes the following steps: 步骤一、机载PHM运行平台建模Step 1: Modeling the airborne PHM operation platform STEP1.机载分布式计算资源建模STEP 1. Modeling of airborne distributed computing resources 机载PHM系统由许多能够运行PHM算法的计算资源组成,假设某一机载PHM系统共有M个可用于运行PHM算法的计算资源;The airborne PHM system consists of many computing resources that can run the PHM algorithm. Assume that a certain airborne PHM system has a total of M computing resources that can be used to run the PHM algorithm; (1)计算资源分类CType (1) Computing Resource Classification C Type 这些计算资源是不同的嵌入式平台,具有不同的处理器架构,每个计算资源是同构系统或异构系统;根据机载PHM系统的计算资源所应用的嵌入式平台以及系统结构,对计算资源进行分类,假设全机PHM系统中的计算资源共有n类,为每一类资源分配一个标号i=1,2,…n,如表1所示;These computing resources are different embedded platforms with different processor architectures. Each computing resource is a homogeneous system or a heterogeneous system. The computing resources are classified according to the embedded platform and system structure used by the computing resources of the airborne PHM system. Assuming that there are n types of computing resources in the whole PHM system, a label is assigned to each type of resource. i = 1, 2, ... n, as shown in Table 1; 表1机载PHM计算资源分类Table 1 Classification of airborne PHM computing resources 其中,指代单ARM架构的计算资源,/>指代双DSP架构的计算资源,/>指代DSP+FPGA异构架构的计算资源,/>指代ARM+DSP+FPGA异构架构的计算资源;目前嵌入式平台发展的趋势是异构平台,在不同的功能需求下经常有新的嵌入式异构系统架构出现,无法穷举所有的架构类型,此处下标的规律是按照单核架构、多核同构架构、简单异构架构、复杂异构架构的类型排序,当有新的嵌入式系统架构应用到机载PHM系统中时,在类代号中继续补充新的下标序号即可;in, Refers to computing resources of a single ARM architecture, /> Refers to the computing resources of the dual DSP architecture, /> Refers to the computing resources of DSP+FPGA heterogeneous architecture,/> Refers to the computing resources of ARM+DSP+FPGA heterogeneous architecture; the current trend of embedded platform development is heterogeneous platform. New embedded heterogeneous system architectures often appear under different functional requirements. It is impossible to list all architecture types. The subscript rules here are sorted according to the types of single-core architecture, multi-core homogeneous architecture, simple heterogeneous architecture, and complex heterogeneous architecture. When a new embedded system architecture is applied to the airborne PHM system, the new subscript number can be added to the class code; (2)算力分级CAbility (2) Computing Power Classification C Ability 在对计算资源分类的基础上,针对每一类计算资源,按照其计算能力进行分级,规定第i号计算资源的算力等级为i=1,2,…,M;假设全机PHM系统中吞吐量最小的计算资源为第l号计算资源,则规定其算力等级/>并规定/>为单位计算能力;假设全机PHM系统中吞吐量最大的计算资源为第k号计算资源,参照单位计算能力的定义,则规定其算力等级/>则CAbilitymax即为全机PHM系统中计算资源的算力等级的最大值;计算资源参照单位计算能力向下取整,得到全机机载PHM系统所有PHM算法模型运行平台对应的算力等级/>i=1,2,…,M;On the basis of the classification of computing resources, each type of computing resource is graded according to its computing power, and the computing power level of the i-th computing resource is specified as i = 1, 2, ..., M; Assuming that the computing resource with the smallest throughput in the whole PHM system is the computing resource No. 1, then its computing power level is specified/> And stipulate/> is the unit computing power; assuming that the computing resource with the largest throughput in the whole PHM system is the kth computing resource, refer to the definition of unit computing power to define its computing power level/> Then C Abilitymax is the maximum value of the computing power level of the computing resources in the whole aircraft PHM system; the computing resources are rounded down with reference to the unit computing power to obtain the computing power level corresponding to all PHM algorithm model running platforms of the whole aircraft PHM system/> i = 1, 2, ..., M; 依据表1以及算力分级,机载PHM系统中的每一个计算资源都可以对应i=1,2,…,n,j=1,2,…,M,由此得到表2;According to Table 1 and the computing power classification, each computing resource in the airborne PHM system can correspond to and i=1,2,…,n,j=1,2,…,M, thus we get Table 2; 表2机载PHM计算资源性能Table 2 Performance of airborne PHM computing resources STEP2.PHM算法性能建模STEP 2. PHM algorithm performance modeling 对于某个机载部件,需要对PHM算法的运行需求进行建模;For a certain airborne component, the operational requirements of the PHM algorithm need to be modeled; (1)确定实时性:PHM工程人员在进行PHM算法开发时,根据算法处理数据的时效性以及算法是否需要使用历史数据,将PHM算法分为在线算法和离线算法;处理的数据具有时效性且不需要使用历史数据的算法模型即为具有实时性的在线算法;处理的数据反应的是部件一段时间内的状态参数,或者需要使用历史数据的算法模型即为不具有实时性的离线算法;(1) Determine real-time performance: When developing PHM algorithms, PHM engineers divide PHM algorithms into online algorithms and offline algorithms based on the timeliness of the data processed by the algorithm and whether the algorithm needs to use historical data. The algorithm model that processes data with timeliness and does not need to use historical data is an online algorithm with real-time performance; the algorithm model that processes data that reflects the status parameters of the component over a period of time or needs to use historical data is an offline algorithm that does not have real-time performance. (2)确定计算资源需求:对于需要在机载PHM系统上运行的算法α,PHM工程人员在进行PHM算法开发时,需要针对算法的复杂度、逻辑结构、数据长度以及数据缓存等需求,通过算法在各类计算资源、各个算力等级的计算资源上进行测试,根据测试结果,确定各类资源上该算法在类资源上运行时,最少需要/>α代表不同的算法,i=1,2,…,n,算力等级的计算资源才能满足算法运行的需求,由此能够为每个算法可以得到一个运行需求,即为表3;(2) Determine the computing resource requirements: For the algorithm α that needs to be run on the airborne PHM system, PHM engineers need to test the algorithm on various computing resources and computing resources of various computing power levels based on the algorithm complexity, logical structure, data length, and data cache requirements when developing the PHM algorithm. Based on the test results, determine the algorithm's performance on various resources. When running on a class resource, at least/> α represents different algorithms, i = 1, 2, ..., n, and computing resources of a certain computing power level can meet the requirements of algorithm operation. Therefore, an operation requirement can be obtained for each algorithm, which is shown in Table 3; 表3算法α运行需求Table 3 Algorithm α operation requirements STEP3.部件监测需求建模STEP 3. Component monitoring requirements modeling 假设:飞机可执行m种飞行任务,飞机的机载PHM系统一共对x个机载部件进行状态监测,机载PHM系统获取其健康状态参数以表征部件的健康状态;首先构建参数以表征部件的监测需求;Assumption: An aircraft can perform m types of flight missions. The aircraft's onboard PHM system monitors the condition of x onboard components in total. The onboard PHM system obtains its health status parameters to characterize the health status of the components. First, the parameters are constructed to characterize the monitoring requirements of the components. (1)部件重要度i=1,2,…,x:(1) Component Importance i = 1, 2, ..., x: 根据部件失效对飞行安全产生的危害度对部件进行重要度的划分,部件失效对飞行安全危害越大则越大;The importance of components is divided according to the degree of harm caused by component failure to flight safety. The greater the harm caused by component failure to flight safety, the greater the importance of components. The bigger; (2)任务重要度j=1,2,…,m:(2) Task importance j = 1, 2, ..., m: 飞机的每一飞行任务都有一个编号j,j=1,2,…,m,其中m为飞机可执行飞行任务种类,根据部件失效对飞行任务完成度的危害对部件进行部件的任务重要度进行划分,部件失效对某一飞行任务的完成度危害越大则越大,执行不同的飞行任务时同一个部件的/>可能相同;Each flight mission of the aircraft has a number j, j = 1, 2, ..., m, where m is the type of flight mission that the aircraft can perform. The components are divided into task importance according to the harm of component failure to the completion of the flight mission. The greater the harm of component failure to the completion of a certain flight mission, the greater the The larger the size, the more likely it is that the same component will perform different flight missions. Probably the same; (3)部件健康等级NHealth(t):(3) Component health level N Health (t): t表示NHealth(t)是时间的函数,随着时间的推移,部件的性能逐渐退化,NHealth(t)改变;t represents N Health (t) as a function of time. As time goes by, the performance of the component gradually degrades and N Health (t) changes; 机载部件的失效是逐渐退化的过程,在部件退化的初期,其退化过程往往是平稳的,基于当前对部件健康状态等级的评估对部件进行健康等级的划分,部件失效的可能性越大则NHealth(t)越大;The failure of airborne components is a gradual degradation process. In the early stage of component degradation, the degradation process is often stable. The components are divided into health levels based on the current assessment of the component health status level. The greater the possibility of component failure, the greater N Health (t); (4)监测需求NMonitor(t):(4) Monitoring requirement N Monitor (t): 所构建的监测需求NMonitor(t),表征某一个部件对于监测的需求迫切度;The constructed monitoring requirement N Monitor (t) represents the urgency of a component’s need for monitoring; STEP4.设置PHM算法运行方式STEP 4. Set the PHM algorithm operation mode 对于某个部件,工程人员在其PHM算法的开发阶段,需要根据其退化模型为其选择若干种算法,在不同的监测需求下对应不同的监测采样率和故障诊断与寿命预测算法;For a certain component, engineers need to select several algorithms for it according to its degradation model during the development phase of its PHM algorithm, corresponding to different monitoring sampling rates and fault diagnosis and life prediction algorithms under different monitoring requirements; 算法、采样率、监测间隔的设置具有以下规律:The settings of the algorithm, sampling rate, and monitoring interval follow the following rules: ⑦监测需求较低时,所选择的算法具备低功耗、低算力需求的特点,算法可以是离线算法或者是在线算法,传感器的采样率较低,监测间隔时间长;算法的目标是监测并发现部件的早期故障特征;⑦ When the monitoring demand is low, the selected algorithm has the characteristics of low power consumption and low computing power requirements. The algorithm can be an offline algorithm or an online algorithm. The sampling rate of the sensor is low and the monitoring interval is long. The goal of the algorithm is to monitor and discover the early fault characteristics of the components. ⑧随着监测需求的增大,所选择的算法的实时性、计算精度提高,以在线算法为主,传感器的采样率提高,监测间隔时间缩短甚至为持续监测;算法的目标是在部件出现早期故障特征之后,准确识别部件的加速退化过程;⑧ As the monitoring demand increases, the real-time performance and calculation accuracy of the selected algorithm are improved, with the online algorithm as the main one, the sampling rate of the sensor is increased, the monitoring interval is shortened or even continuous monitoring; the goal of the algorithm is to accurately identify the accelerated degradation process of the component after the early fault characteristics of the component appear; ⑨当监测需求继续增大,应当选择具备高精度和实时性特点的算法,传感器采样率达到最高,监测方式为实时监测;算法的目的是在部件的加速退化阶段准确预测部件的剩余使用寿命,在部件发生故障时能够立即对故障进行准确的识别和定位;⑨ When the monitoring demand continues to increase, an algorithm with high precision and real-time characteristics should be selected, the sensor sampling rate should be the highest, and the monitoring method should be real-time monitoring; the purpose of the algorithm is to accurately predict the remaining service life of the component during the accelerated degradation stage of the component, and to be able to immediately and accurately identify and locate the fault when the component fails; 在所有x个机载部件中,最高的监测需求为PHM开发人员为每个部件在每个确定的NMonitor(t)下找到一种PHM算法,算法(1,0)表示部件1在NMonitor(t)=0时运行的算法,算法(1,1)表示部件1在NMonitor(t)=1时运行的算法,以此类推,算法(x,N)表示部件x在NMonitor(t)=N时运行的算法;PHM开发人员为每个部件在每个确定的NMonitor(t)下都设定PHM算法运行时数据的采样率和算法两次连续运行之间的时间间隔;得到表4;Among all x airborne components, the highest monitoring requirement is The PHM developer finds a PHM algorithm for each component under each determined N Monitor (t). Algorithm (1, 0) represents the algorithm that component 1 runs when N Monitor (t) = 0, and algorithm (1, 1) represents the algorithm that component 1 runs when N Monitor (t) = 1. Similarly, algorithm (x, N) represents the algorithm that component x runs when N Monitor (t) = N. The PHM developer sets the sampling rate of the PHM algorithm runtime data and the time interval between two consecutive runs of the algorithm for each component under each determined N Monitor (t). Table 4 is obtained. 表4 PHM算法运行方式Table 4 PHM algorithm operation mode STEP5.确定算法运行性能STEP 5. Determine the algorithm performance 根据表4中的采样率和监测间隔,PHM工程人员在进行PHM算法开发时,需要对PHM算法在不同的计算资源下的性能进行测试,测试内容为:对于算法α,当其运行在类,计算能力为/>的计算资源下时,可以确定其单次执行的计算延迟/>和功耗得到表5;According to the sampling rate and monitoring interval in Table 4, PHM engineers need to test the performance of the PHM algorithm under different computing resources when developing the PHM algorithm. The test content is: for algorithm α, when it runs on Class, computing power is/> When the computing resources are sufficient, the computing delay of a single execution can be determined./> and power consumption Table 5 is obtained; 表5 PHM算法α运行性能Table 5 PHM algorithm α running performance STEP6.确定PHM系统运行需求STEP 6. Determine the PHM system operation requirements 当机载PHM系统运行时,能够根据每个机载部件的当前时刻的监测需求确定当前时刻整个PHM系统中需要运行的所有算法及其采样率、运行间隔和算法运行性能,且假设每个PHM算法都只在一个计算资源上运行;假设t时刻一共有l个需要运行的PHM算法,则这些需要运行的PHM算法及其采样率、运行间隔和算法运行性能构成一个集合;就确定出t时刻机载PHM系统的运行需求和运行性能,即机载PHM系统运行的需求及任务分配的空间;When the airborne PHM system is running, all algorithms that need to be run in the entire PHM system at the current moment and their sampling rates, running intervals and algorithm running performance can be determined according to the monitoring requirements of each airborne component at the current moment, and it is assumed that each PHM algorithm is only run on one computing resource; assuming that there are a total of l PHM algorithms that need to be run at time t, then these PHM algorithms that need to be run and their sampling rates, running intervals and algorithm running performance constitute a set; the running requirements and running performance of the airborne PHM system at time t are determined, that is, the running requirements of the airborne PHM system and the space for task allocation; 步骤二、PHM计算任务分配优化空间构造Step 2: PHM calculation task allocation and optimization of spatial structure 假设某一时刻,机载PHM系统需要同时运行K个算法;Assume that at a certain moment, the airborne PHM system needs to run K algorithms simultaneously; 在获得PHM系统运行的需求及任务分配的空间的基础上,在任务分配空间内对PHM计算任务分配优化;Based on the requirements of PHM system operation and the space of task allocation, the PHM computing task allocation is optimized within the task allocation space; 以机载PHM系统中所有算法执行的延迟和功耗作为评价指标构造性能指标函数J(D):The performance index function J(D) is constructed by taking the delay and power consumption of all algorithms executed in the airborne PHM system as evaluation indicators: 其中,in, 为K维列向量,表示一共有K个需要运行的算法,分别分配到n1、n2、…、nK号计算资源上运行; is a K-dimensional column vector, indicating that there are K algorithms that need to be run, which are respectively assigned to n 1 , n 2 , …, n K computing resources for execution; 和/>分别表示第i个算法在当前分配的计算资源中执行的延迟时间和功耗; and/> They represent the delay time and power consumption of the i-th algorithm executed in the currently allocated computing resources respectively; k1、k2为权重因子,表示算法执行延迟和计算功耗在性能指标中的权重,满足k1+k2=1,k1、k2≥0;k 1 and k 2 are weight factors, indicating the weights of algorithm execution delay and computing power consumption in the performance index, satisfying k 1 +k 2 =1, k 1 , k 2 ≥0; 约束条件为:算法需要在能够在满足最低算力等级需求的计算资源上运行;The constraints are: the algorithm needs to run on computing resources that can meet the minimum computing power level requirements; 在上述分配优化条件的基础上,再利用寻优算法,将每个需要运行的PHM算法分配到一个计算资源上,使得上述分配优化条件中的性能指标函数J(D)最小,即达到PHM系统算法运行的总体性能最优。On the basis of the above allocation optimization conditions, the optimization algorithm is used to allocate each PHM algorithm that needs to be run to a computing resource, so that the performance indicator function J(D) in the above allocation optimization conditions is minimized, that is, the overall performance of the PHM system algorithm operation is optimal.
CN202210115309.XA 2022-01-30 2022-01-30 A method for airborne distributed PHM computational modeling Active CN114925441B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210115309.XA CN114925441B (en) 2022-01-30 2022-01-30 A method for airborne distributed PHM computational modeling

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210115309.XA CN114925441B (en) 2022-01-30 2022-01-30 A method for airborne distributed PHM computational modeling

Publications (2)

Publication Number Publication Date
CN114925441A CN114925441A (en) 2022-08-19
CN114925441B true CN114925441B (en) 2024-06-07

Family

ID=82805435

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210115309.XA Active CN114925441B (en) 2022-01-30 2022-01-30 A method for airborne distributed PHM computational modeling

Country Status (1)

Country Link
CN (1) CN114925441B (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA3058203A1 (en) * 2017-06-15 2018-12-20 Aurora Flight Sciences Corporation Autonomuos aircraft health systems and methods
CN109917772A (en) * 2017-12-13 2019-06-21 北京航空航天大学 A PHM Rapid Prototyping System for Remote Online Evaluation of Equipment Status
CN110544328A (en) * 2019-09-02 2019-12-06 哈尔滨工业大学 An embedded high-efficiency computing platform and method for airborne fault prediction and health management
CN112839382A (en) * 2020-12-30 2021-05-25 北京邮电大学 A joint allocation method of communication and computing resources driven by video semantics in the Internet of Vehicles
WO2021218003A1 (en) * 2020-04-27 2021-11-04 中国电子科技集团公司第十四研究所 Radar embedded health management system

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9424693B2 (en) * 2014-03-10 2016-08-23 Embraer S.A. Maintenance planning optimization for repairable items based on prognostics and health monitoring data

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA3058203A1 (en) * 2017-06-15 2018-12-20 Aurora Flight Sciences Corporation Autonomuos aircraft health systems and methods
CN109917772A (en) * 2017-12-13 2019-06-21 北京航空航天大学 A PHM Rapid Prototyping System for Remote Online Evaluation of Equipment Status
CN110544328A (en) * 2019-09-02 2019-12-06 哈尔滨工业大学 An embedded high-efficiency computing platform and method for airborne fault prediction and health management
WO2021218003A1 (en) * 2020-04-27 2021-11-04 中国电子科技集团公司第十四研究所 Radar embedded health management system
CN112839382A (en) * 2020-12-30 2021-05-25 北京邮电大学 A joint allocation method of communication and computing resources driven by video semantics in the Internet of Vehicles

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于云计算的飞机PHM体系架构研究;李耀华;尚金秋;;计算机工程;20171231(12);12-16 *
数据驱动视角下飞机故障预测与健康管理系统设计及验证;刘亮;周博;于涛;张宁;;计算机测量与控制;20170725(07);21-25 *

Also Published As

Publication number Publication date
CN114925441A (en) 2022-08-19

Similar Documents

Publication Publication Date Title
EP3847549B1 (en) Minimizing impact of migrating virtual services
CN103729248B (en) A kind of method and apparatus of determination based on cache perception task to be migrated
CN108595301B (en) A method and system for predicting server energy consumption based on machine learning
US20080320269A1 (en) Method and apparatus for ranking of target server partitions for virtual server mobility operations
CN103279392B (en) A kind of load sorting technique run on virtual machine under cloud computing environment
CN118550711A (en) Method and system for improving calculation efficiency
CN103955398B (en) Virtual machine coexisting scheduling method based on processor performance monitoring
US20080320123A1 (en) Process and methodology for generic analysis of metrics related to resource utilization and performance
CN103559077B (en) Optimized automatic transfer method and system for virtual machines
US11042209B2 (en) Control of the energy consumption of a server cluster
CN118153516B (en) Processor function simulation verification method and system
CN115904666A (en) Deep learning training task scheduling system for GPU clusters
CN114925441B (en) A method for airborne distributed PHM computational modeling
Swain et al. Efficient straggler task management in cloud environment using stochastic gradient descent with momentum learning-driven neural networks
KR20190078453A (en) Migration System and Method by Fuzzy Value Rebalance in Distributed Cloud Environment
CN116185584A (en) Multi-tenant database resource planning and scheduling method based on deep reinforcement learning
US20110055831A1 (en) Program execution with improved power efficiency
WO2024046283A1 (en) Task scheduling method, model generation method, and electronic device
CN116360921A (en) Cloud platform resource optimal scheduling method and system for electric power Internet of things
KR100547625B1 (en) Intelligent Monitoring System and Method for Grid Information Service
Wang et al. HARRD: Real-time software rejuvenation decision based on hierarchical analysis under weibull distribution
CN114326475A (en) Reliability optimization method of unmanned domain controller based on competitive failure mode
Huang et al. Is the powersave governor really saving power?
CN114048013B (en) Offline task scheduling optimization method and device
CN119576589B (en) A control method, device and medium for a server device

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