CN111814360B - Method for evaluating use reliability of airborne electromechanical equipment - Google Patents
Method for evaluating use reliability of airborne electromechanical equipment Download PDFInfo
- Publication number
- CN111814360B CN111814360B CN202010859077.XA CN202010859077A CN111814360B CN 111814360 B CN111814360 B CN 111814360B CN 202010859077 A CN202010859077 A CN 202010859077A CN 111814360 B CN111814360 B CN 111814360B
- Authority
- CN
- China
- Prior art keywords
- fault
- mechanical
- follows
- electronic
- distribution
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/20—Administration of product repair or maintenance
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/12—Timing analysis or timing optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- General Physics & Mathematics (AREA)
- Entrepreneurship & Innovation (AREA)
- Quality & Reliability (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Operations Research (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Economics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
The invention discloses a reliability evaluation method for use of airborne electromechanical equipment, which comprises the steps of analyzing a fault mode of a specific product, counting fault data of mechanical and electronic components according to the fault mode, further fitting a product life distribution model, estimating distribution model parameters, respectively calculating fault weights of mechanical and electronic components of a system, then constructing a mixed distribution model of the system, finally drawing a system reliability curve of target equipment, then accurately reflecting the reliability level of the components, searching weak links of the system, accurately reflecting the overall reliability level of the system, and providing a basis for reducing the risk rate of the system, designing preventive maintenance cycles and optimizing maintenance cost.
Description
Technical Field
The invention relates to the technical field of reliability engineering, in particular to a method for evaluating the use reliability of airborne electromechanical equipment.
Background
The airborne electromechanical product is assisted by the severe environment, is easy to be subjected to the comprehensive action of various environmental stresses such as high temperature, high humidity and acid gas and various mechanical stresses such as vibration and impact, and has complex fault behaviors. In the design, development and production stages, due to various constraints such as cost, environment, test and calculation methods and the like, the reliability prediction and identification of the airborne product cannot completely reflect the reliability level of the product in the real environment, and huge hidden danger is caused to the reliability of the airborne electromechanical product after being put into use.
Currently, the reliability evaluation of the airborne electromechanical products is usually performed on the assumption that the airborne electromechanical products are subjected to a certain distribution, but the airborne products are generally formed by combining mechanical components and electronic components, and the fault distribution models of the mechanical components and the electronic components often do not belong to the same distribution type. Meanwhile, the existing reliability comprehensive method is often to simplify the system into a series model according to a reliability logic method, but the traditional reliability mathematical model of the series system has the defect that the reliability of the system is more conservative when being calculated, so that the reliability analysis result is not accurate enough.
Disclosure of Invention
The invention aims to provide a method for evaluating the use reliability of airborne electromechanical equipment, which can accurately reflect the reliability level of components, search weak links of a system, accurately reflect the overall reliability level of the system and provide a basis for reducing the risk rate of the system, designing preventive maintenance cycles and optimizing the maintenance cost.
The technical scheme adopted by the invention is as follows:
a method for evaluating the use reliability of airborne electromechanical equipment comprises the following steps:
s1: establishing a target equipment sample set; the method comprises the following specific steps:
s 1.1: counting fault data of target equipment, and constructing a fault sample database X ═ X 1 ,…,x M M is the number of fault samples;
s 1.2: analyzing the fault mode of the fault data in the fault sample database and classifying the fault data to form a mechanical sample database X 1 And an electronic sample database X 2 (ii) a The method specifically comprises the following steps:
s2: fitting a mechanical component service life distribution model and calculating a mechanical failure weight omega 1 :
S3: fitting an electronic component service life distribution model and calculating an electronic fault weight omega 2 ;
S4: constructing a mixed distribution model of target equipment, and acquiring the system reliability R (t) of the model;
in formula (1), R (t) represents system reliability, F (t) represents system fault distribution function, F 1 (t) represents a mechanical assembly failure distribution function, F 2 (t) represents the electronic component failure distribution function, Φ 1 Normal distribution function is expressed, t represents operation time, mu and sigma are mean and variance of mechanical assembly life distribution model, and alpha is electronic assembly lifeThe shape parameter of the distribution model, beta is the scale parameter of the electronic component service life distribution model;
s5: and drawing a system reliability curve of the target equipment.
Further, the specific process of step s1.2 is as follows:
s 1.2.1: counting fault data of mechanical parts in target equipment, and constructing a mechanical sample databaseM 1 Number of mechanical samples;
s 1.2.2: counting fault data of electronic components in target equipment, and constructing an electronic sample databaseM 2 Number of electronic samples, M 2 =M-M 1 。
Further, the step S2 specifically includes the following steps:
s 2.1: calculating distribution parameters of a service life distribution model of the mechanical assembly;
the service life of the mechanical assembly follows the log normal distribution, and the fault function of the mechanical assembly is as follows:
in the formula (2), phi 1 Representing a normal distribution function, t representing the operation time, and mu and sigma being the mean and variance of the life distribution function of the mechanical component;
the solution equations for μ and σ are as follows:
in the formula (3), lnL represents a log-likelihood function, t (i) For fault sample X 1 N is the total amount of observation samples, Z 0 =(lnt 0 - μ)/σ as statistic, Φ (-Z) 0 ) Is a standard normal distribution functionNumber phi (Z) 0 ) Is a standard normal distribution density function;
s 2.2: calculating a mechanical failure weight ω 1 ;
In the formula (4), M is the total number of system faults, M 1 The number of mechanical samples.
Further, the step S3 specifically includes the following steps:
s 3.1: calculating distribution parameters of the electronic component service life distribution model;
the electronic component life follows a weibull distribution, and the electronic component failure function is as follows:
in the formula (5), alpha is a shape parameter, and beta is a scale parameter;
the shape parameter alpha and the scale parameter beta are solved through a likelihood equation, and the solving equation is as follows:
in the formula (6), t (i) For fault sample X 2 N is the total amount of observation samples, M 2 The number of failures of the electronic component;
s 3.2: calculating an electronic fault weight ω 2 ;
ω 2 =1-ω 1 (7);
In the formula (7), ω 1 Is the mechanical failure weight.
Further, the fault sample database is established based on fault time.
The invention has the following beneficial effects:
the method comprises the steps of counting fault data of mechanical and electronic components according to fault modes by analyzing the fault modes of specific products, further fitting a product life distribution model, calculating distribution model parameters, calculating fault weights of mechanical and electronic components of a system respectively, constructing a mixed distribution model of the system, finally drawing a system reliability curve of target equipment, accurately reflecting the reliability level of the components, searching weak links of the system, accurately reflecting the overall reliability level of the system, and providing a basis for reducing the risk rate of the system, designing preventive maintenance periods and optimizing maintenance cost.
Drawings
FIG. 1 is a system reliability curve in accordance with an embodiment.
Detailed Description
The method for evaluating the use reliability of the airborne electromechanical equipment comprises the following steps: the method comprises the following steps:
s1: a sample set of target devices is established.
The method comprises the following specific steps:
s 1.1: counting fault data of target equipment, and constructing a fault sample database X ═ X 1 ,…,x M M is the number of fault samples;
s 1.2: analyzing the fault mode of the fault data in the fault sample database and classifying the fault data to form a mechanical sample database X 1 And an electronic sample database X 2 (ii) a The method comprises the following specific steps:
s 1.2.1: counting fault data of mechanical parts in target equipment, and constructing a mechanical sample databaseM 1 Number of mechanical samples;
s 1.2.2: counting fault data of electronic components in target equipment, and constructing an electronic sample databaseNumber of electronic samples, M 2 =M-M 1 。
S2: service life distribution die for fitting mechanical assemblyType, calculating the mechanical failure weight ω 1 。
The specific process is as follows:
s 2.1: calculating distribution parameters of a service life distribution model of the mechanical assembly;
the service life of the mechanical assembly follows the lognormal distribution, and the fault function of the mechanical assembly is as follows:
in the formula (1), phi 1 Representing a normal distribution function, t representing the running time, and mu and sigma representing the mean and variance of a mechanical component life distribution model;
the solution equations for the distribution parameters μ and σ are as follows:
in the formula (2), lnL is a log-likelihood function, t (i) For fault sample X 1 N is the total amount of observation samples, Z 0 =(lnt 0 - μ)/σ as statistic, Φ (-Z) 0 ) Is a standard normal distribution function, phi (Z) 0 ) Is a standard normal distribution density function;
s 2.2: calculating a mechanical failure weight ω 1 ;
In formula (3), M is the total number of system faults, M 1 The number of mechanical samples.
S3: fitting an electronic component service life distribution model and calculating an electronic fault weight omega 2 。
The specific process is as follows:
s 3.1: calculating distribution parameters of the electronic component service life distribution model;
the electronic component life follows a weibull distribution, and the electronic component failure function is as follows:
in the formula (4), alpha is a shape parameter, and beta is a scale parameter;
the shape parameter alpha and the scale parameter beta are solved through a likelihood equation, and the solving equation is as follows:
in the formula (5), t (i) For fault sample X 2 N is the total amount of observation samples, M 2 Is the number of failures of the electronic component;
s 3.2: calculating an electronic fault weight ω 2 ;
ω 2 =1-ω 1 (6);
In the formula (6), ω 1 Is the mechanical failure weight.
S4: constructing a mixed distribution model of the target equipment;
in the formula (7), R (t) represents the system reliability, F (t) represents the system fault distribution function, and F 1 (t) represents a mechanical assembly failure distribution function, F 2 (t) represents the electronic component failure distribution function, Φ 1 Representing a normal distribution function, t representing the running time, mu and sigma representing the mean value and variance of the mechanical assembly service life distribution model, alpha representing the shape parameter of the electronic assembly Weibull distribution model, and beta representing the scale parameter of the electronic assembly service life distribution model;
s5: and drawing a system reliability curve of the target equipment.
For a better understanding of the present invention, the following embodiments are provided to further explain the technical solutions of the present invention.
The invention comprises the following steps:
s1: and (4) sorting and classifying system fault data in a service environment, and establishing a target equipment sample set.
The method comprises the following specific steps:
s 1.1: counting fault data of target equipment, and constructing a fault sample database X ═ X 1 ,…,x M M is the number of fault samples;
s 1.2: analyzing the fault mode of the fault data in the fault sample database and classifying the fault data to form a mechanical sample database X 1 And an electronic sample database X 2 (ii) a The method specifically comprises the following steps:
s 1.2.1: counting fault data of mechanical parts in target equipment based on fault time, and constructing a mechanical sample databaseM 1 Number of mechanical samples;
s 1.2.2: counting fault data of electronic components in target equipment based on fault time, and constructing an electronic sample databaseM 2 Number of electronic samples, M 2 =M-M 1 。
In this embodiment, a system sample 25 is observed, wherein a fault data sample 20 is present. The fault sample data is shown in table 1:
table 1:
s2: fitting a mechanical component service life distribution model and calculating a mechanical fault weight omega 1 。
The specific process is as follows:
s 2.1: and calculating distribution parameters of the mechanical component service life distribution model.
The service life of the mechanical assembly follows the log normal distribution, and the fault function of the mechanical assembly is as follows:
in the formula (2), phi 1 Normal distribution function is expressed, t represents the running time, and μ and σ are the mean and variance of the mechanical component life distribution model.
The solution equations for μ and σ are as follows:
in equation (3), lnL is a log-likelihood function, t (i) For fault sample X 1 N is the total amount of observation samples, Z 0 =(lnt 0 - μ)/σ as statistic, Φ (-Z) 0 ) Is a standard normal distribution function, phi (Z) 0 ) Is a standard normal distribution density function. s 2.2: calculating a mechanical failure weight ω 1 。
In the formula (4), M is the total number of system faults, M 1 The number of mechanical samples.
The calculation result of this embodiment is as follows:
s3: fitting an electronic component service life distribution model and calculating an electronic fault weight omega 2 。
The specific process is as follows:
s 3.1: calculating distribution parameters of the electronic component service life distribution model;
the service life of the electronic component follows Weibull distribution, and the fault function of the electronic component is as follows:
in the formula (5), alpha is a shape parameter, and beta is a scale parameter;
the shape parameter alpha and the scale parameter beta are solved through a likelihood equation, and the solving equation is as follows:
in the formula (6), t (i) For fault sample X 2 N is the total amount of observation samples, M 2 The number of failures of the electronic component;
s 3.2: calculating an electronic fault weight ω 2 ;
ω 2 =1-ω 1 (7);
In the formula (7), ω 1 Is the mechanical failure weight.
The calculation results in this example are as follows:
Shape parameter alpha is 0.931, and scale parameter beta is 9.304
S4: and constructing a mixed distribution model of the target equipment.
In formula (1), r (t) represents the system reliability, and t represents the operating time.
S5: and calculating the mean time between failures of the system, and drawing a system reliability curve of the target equipment as shown in figure 1.
According to the method, service life distribution models of the mechanical component and the electronic component are respectively fitted based on fault data generated by the airborne electromechanical equipment in the service process, then a new system reliability model is constructed based on the mixed distribution model on the basis of considering the fault weight of each component, and finally a system reliability curve of the target equipment is drawn, so that the reliability level of the component is accurately reflected, weak links of the system are searched, the overall reliability level of the system can be accurately reflected, and a basis is provided for reducing the risk rate of the system, designing a preventive maintenance period and optimizing the maintenance cost.
Claims (5)
1. The use reliability assessment method of the airborne electromechanical equipment is characterized by comprising the following steps: the method comprises the following steps:
s1: establishing a target equipment sample set; the method comprises the following specific steps:
s 1.1: counting fault data of target equipment, and constructing a fault sample database X ═ X 1 ,…,x M M is the number of fault samples;
s 1.2: analyzing the fault mode of the fault data in the fault sample database and classifying the fault data to form a mechanical sample database X 1 And an electronic sample database X 2 (ii) a The method specifically comprises the following steps:
s2: fitting a mechanical component service life distribution model and calculating a mechanical fault weight omega 1 :
S3: fitting an electronic component service life distribution model and calculating an electronic fault weight omega 2 ;
S4: constructing a mixed distribution model of target equipment, and acquiring the system reliability R (t) of the model;
in formula (1), R (t) represents system reliability, F (t) represents system fault distribution function, F 1 (t) represents a mechanical assembly failure distribution function, F 2 (t) represents the electronic component failure distribution function, Φ 1 Representing a normal distribution function, t representing the running time, mu and sigma representing the mean value and variance of the mechanical assembly service life distribution model, alpha representing the shape parameter of the electronic assembly service life distribution model, and beta representing the scale parameter of the electronic assembly service life distribution model;
s5: and drawing a system reliability curve of the target equipment.
2. The airborne electromechanical device use reliability assessment method of claim 1, wherein: the specific process of step s1.2 is as follows:
s 1.2.1: counting fault data of mechanical parts in target equipment, and constructing a mechanical sample databaseM 1 Number of mechanical samples;
3. The airborne electromechanical device usage reliability assessment method according to claim 1, characterized in that: the specific process of step S2 is as follows:
s 2.1: calculating distribution parameters of a service life distribution model of the mechanical assembly;
the service life of the mechanical assembly follows the log normal distribution, and the fault function of the mechanical assembly is as follows:
in the formula (2), phi 1 Representing a normal distribution function, t representing the running time, and mu and sigma representing the mean value and the variance of the service life distribution function of the mechanical component;
the solution equations for μ and σ are as follows:
in the formula (3), lnL represents a log-likelihood function, t (i) For fault sample X 1 N is the total amount of observation samples, Z 0 =(lnt 0 - μ)/σ as statistic, Φ (-Z) 0 ) Is a standard normal distribution function, phi (Z) 0 ) Is a standard normal distribution density function;
s 2.2: calculating a mechanical failure weight ω 1 ;
In the formula (4), M is the total number of system faults, M 1 The number of mechanical samples.
4. The airborne electromechanical device usage reliability assessment method according to claim 1, characterized in that: the specific process of step S3 is as follows:
s 3.1: calculating distribution parameters of the electronic component service life distribution model;
the electronic component life follows a weibull distribution, and the electronic component failure function is as follows:
in the formula (5), alpha is a shape parameter, and beta is a scale parameter;
the shape parameter alpha and the scale parameter beta are solved through a likelihood equation, and the solving equation is as follows:
in the formula (6), t (i) For fault sample X 2 N is the total amount of observation samples, M 2 The number of failures of the electronic component;
s 3.2: calculating an electronic fault weight ω 2 ;
ω 2 =1-ω 1 (7);
In the formula (7), ω 1 Is the mechanical failure weight.
5. The airborne electromechanical device usage reliability assessment method according to claim 1, characterized in that: the specific process of step S3 is as follows: the fault sample database is established based on fault time.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010859077.XA CN111814360B (en) | 2020-08-24 | 2020-08-24 | Method for evaluating use reliability of airborne electromechanical equipment |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010859077.XA CN111814360B (en) | 2020-08-24 | 2020-08-24 | Method for evaluating use reliability of airborne electromechanical equipment |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111814360A CN111814360A (en) | 2020-10-23 |
CN111814360B true CN111814360B (en) | 2022-09-30 |
Family
ID=72859675
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010859077.XA Active CN111814360B (en) | 2020-08-24 | 2020-08-24 | Method for evaluating use reliability of airborne electromechanical equipment |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111814360B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112395692B (en) * | 2020-11-30 | 2022-05-24 | 福建省特种设备检验研究院泉州分院 | Disassembly-free bridge crane reliability assessment method |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103971024A (en) * | 2014-05-26 | 2014-08-06 | 华北电力大学(保定) | Method for evaluating reliability of relaying protection systems under small sample failure data |
CN106547265B (en) * | 2016-11-01 | 2019-01-29 | 中车株洲电力机车研究所有限公司 | A kind of live reliability estimation method and system of rail traffic electronic control unit |
CN108459991B (en) * | 2018-03-12 | 2021-12-21 | 哈尔滨工程大学 | Method for obtaining equipment reliability data |
CN111414674A (en) * | 2020-02-24 | 2020-07-14 | 常州工学院 | Reliability evaluation method for logging instrument |
-
2020
- 2020-08-24 CN CN202010859077.XA patent/CN111814360B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN111814360A (en) | 2020-10-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111461555B (en) | Production line quality monitoring method, device and system | |
CN107357275B (en) | Non-gaussian industrial process fault detection method and system | |
CN109613428A (en) | It is a kind of can be as system and its application in motor device fault detection method | |
CN112083244B (en) | Integrated intelligent diagnosis system for faults of avionic equipment | |
CN108898311A (en) | A kind of data quality checking method towards intelligent distribution network repairing dispatching platform | |
CN110174878A (en) | A kind of unmanned intelligence equipment is healthy and ensures integrated management general-purpose system | |
Filz et al. | Simulation-based assessment of quality inspection strategies on manufacturing systems | |
CN111126759A (en) | Electric energy meter state evaluation method based on abnormal event fault correlation degree | |
CN109916634A (en) | A kind of aviation turbofan engine remaining life prediction technique and system | |
CN111814360B (en) | Method for evaluating use reliability of airborne electromechanical equipment | |
CN109145319A (en) | Key equipment cutting tool method for predicting residual useful life based on deep neural network | |
CN106022636A (en) | Astronavigation component PID implementation effect evaluation method | |
CN117131437A (en) | Isolation method for vibration faults of engine | |
CN116773239A (en) | Intelligent gas meter controller reliability life prediction method | |
CN109213119B (en) | Complex industry key component fault prediction method and system based on online learning | |
CN111965995B (en) | Three-in-one aircraft avionics system health assessment and prediction system and method thereof | |
CN111474905B (en) | Parameter drift fault diagnosis method in manufacturing process of electromechanical product | |
CN115455647A (en) | Moonlet service life evaluation method and system based on random process | |
CN109388829B (en) | Electronic product service life measuring and calculating method | |
Li et al. | Simulation of Reliability Prediction Based on Multiple Factors for Spinning Machine | |
Song et al. | An empirical study of comparison of code metric aggregation methods–on embedded software | |
CN111506998B (en) | Method for constructing parameter drift fault feature sample library in manufacturing process of electromechanical product | |
CN111209973B (en) | Process monitoring method based on mixed kernel PCA-CCA and kernel density estimation | |
CN115566997B (en) | Photovoltaic module attenuation test system and method for air compression station | |
Xiaomei et al. | A software intelligence level evaluation method based on analytic hierarchy process and principal component analysis |
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 |