CN111046555A - Time-varying degradation quality characteristic compensation full life cycle quality robustness optimization method - Google Patents

Time-varying degradation quality characteristic compensation full life cycle quality robustness optimization method Download PDF

Info

Publication number
CN111046555A
CN111046555A CN201911267915.8A CN201911267915A CN111046555A CN 111046555 A CN111046555 A CN 111046555A CN 201911267915 A CN201911267915 A CN 201911267915A CN 111046555 A CN111046555 A CN 111046555A
Authority
CN
China
Prior art keywords
quality
life cycle
characteristic
electromagnetic relay
robustness
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.)
Granted
Application number
CN201911267915.8A
Other languages
Chinese (zh)
Other versions
CN111046555B (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.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
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 Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN201911267915.8A priority Critical patent/CN111046555B/en
Publication of CN111046555A publication Critical patent/CN111046555A/en
Application granted granted Critical
Publication of CN111046555B publication Critical patent/CN111046555B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

A time-varying degradation quality characteristic compensation full-life cycle quality robustness optimization method belongs to the field of electromagnetic relay quality optimization design. Analyzing and determining key design parameters of the electromagnetic relay, carrying out linearity analysis, expanding the whole life cycle of the electromagnetic relay based on a K-L expansion method, dividing the whole life cycle of the electromagnetic relay into a plurality of sub-life intervals, modeling quality robustness characteristic parameters of the sub-life intervals by using a Kriging method, carrying out unified expression, establishing a quality robustness model of the whole life cycle of the electromagnetic relay, selecting an appropriate quality characteristic level for a quality robustness characteristic requirement level, calculating the quality robustness characteristic deviation degree of the electromagnetic relay under the action of time-varying degradation parameters, compensating the quality deviation to the requirement of the quality robustness of the whole life cycle, generating batch samples by using Monte Carlo, calculating output characteristics, and verifying and optimizing effects. A new method is provided for solving the design optimization of the quality robustness of the electromagnetic relay product in the whole life cycle.

Description

Time-varying degradation quality characteristic compensation full life cycle quality robustness optimization method
Technical Field
The invention relates to a quality robustness optimization method, in particular to a time-varying degradation quality characteristic compensation full-life cycle quality robustness optimization method, and belongs to the field of electromagnetic relay quality optimization design.
Background
The electromagnetic relay has the typical advantages of excellent universality, extremely high isolation strength, excellent interference resistance strength and the like, and is widely applied to modern equipment such as industrial operation and control systems, aerospace planes, space shuttles, manned and unmanned spacecrafts and the like to complete the functions of control signal transmission, function control execution, energy system power distribution and the like. With the development of modern equipment towards high sensitivity, high power consumption per unit volume and high reliability of the whole life cycle, especially the requirement of high robustness of the whole life cycle is obvious, and correspondingly, the requirement of robustness of the whole life cycle of the electromagnetic relay is improved to an unprecedented height. The development of the full-life-cycle robustness optimization of the electromagnetic relay is a necessary way for ensuring the application reliability of the electromagnetic relay and designing and optimizing the full-life-cycle robustness of the electromagnetic relay.
The robustness optimization design method with the aim of weakening the influence of parameter fluctuation on the product quality consistency can obviously improve the quality robustness in the industrial application process. The most common method for robust design optimization is to combine an experimental design method with a Taguchi method, which can only optimize in a limited range and has a large amount of calculation. In order to solve the problem of dimension and calculation of distribution of quality characteristic calculation and tolerance of the robust design, a substitution model is provided, the problem of quality output characteristic calculation in the robust design process is solved to a certain extent by introducing the substitution model, but the calculation accuracy of the substitution model and the accuracy of a training sample of the substitution model become the problems restricting the optimization development of the robust design. In addition, aiming at the problem that the Taguchi method is only suitable for single-target optimization, robust design optimization methods based on a sensitive domain, a multi-target interaction entropy, an intelligent algorithm and a hierarchical multi-target are provided, the methods better realize the requirements of multi-target robust design, but the methods mainly pay attention to the mutual influence among output characteristics and often ignore the influence of the correlation among input parameters on the optimal robust solution. It should be noted that under the comprehensive action of multiple parameters such as mechanical-electrical-magnetic-thermal-vibration-impact and the like in the actual working process, a large number of time-varying degradation parameters exist in the parameters of key parts of the product, and the time-varying performance degradation affects the whole life quality robustness of the product, so that the mean value and the variance of the quality robustness characteristic of the permanent magnet actuator are greatly changed. The robustness design method based on initial (time invariant) design is optimized aiming at initial quality design, and the analysis and design optimization of the time variant (dynamic) robustness of the whole life cycle of the permanent magnet actuator are rarely researched.
Aiming at the influence of the performance degradation of parts in the whole life cycle of the electromagnetic relay on the quality robustness, the invention provides a method for modeling the quality robustness characteristic of the whole life cycle of the electromagnetic relay based on K-L expansion and a Krigin model and optimizing the anti-degradation quality robustness of the whole life cycle of the electromagnetic relay based on time-varying degradation quality characteristic compensation, thereby realizing the optimization of the quality robustness of the whole life cycle of the electromagnetic relay.
Disclosure of Invention
The invention aims to solve the problem of missing of a full-life-cycle degradation-resistant optimization design method under the degradation action of an electromagnetic relay, provides a full-life-cycle quality robustness optimization method for time-varying degradation quality characteristic compensation, and particularly relates to how to perform quality robustness characteristic degradation-resistant optimization on the electromagnetic relay under the condition that the performance of parts of the electromagnetic relay is continuously degraded along with the working to obtain the electromagnetic relay with the full-life-cycle quality robustness so as to overcome the problems in the background technology.
In order to achieve the purpose, the invention adopts the following technical scheme: the method for optimizing the quality robustness of the full life cycle of the time-varying degradation quality characteristic compensation comprises the following steps:
the method comprises the following steps: analyzing and determining key design parameters of the electromagnetic relay, and recording the key design parameters as a set U-U (U)1,u2,…,un);
Step two: carrying out linearity analysis on the key design parameters of the electromagnetic relay determined by analysis in the step one, and screening out key nonlinear designsThe parameter is expressed as the set V ═ V (V)1,v2,…,vm);
Step three: the method comprises the steps that the whole life cycle of the electromagnetic relay is unfolded based on a K-L unfolding method, and the whole life cycle is divided into a plurality of sub-life intervals;
step four: modeling the quality robustness characteristic parameters of the sub-life intervals by using a Kriging method, carrying out unified expression on each sub-life interval, and establishing a quality robustness model of the whole life cycle of the electromagnetic relay;
step five: selecting a proper quality characteristic level according to the quality robustness characteristic requirement level of the whole life cycle of the electromagnetic relay;
step six: calculating the deviation degree of the electromagnetic relay quality robustness characteristic under the action of the time-varying degradation parameter based on the quality robustness level requirement obtained by the analysis in the fifth step
Figure BDA0002313388020000031
Wherein b is a quality robustness level;
step seven: compensation for quality robust feature bias using nonlinear parameter combinations, i.e., let
Figure BDA0002313388020000032
Thereby compensating the mass offset to the full life cycle mass robustness requirement;
step eight: and generating batch samples by using Monte Carlo according to the optimized parameter combination of the electromagnetic relay, calculating the output characteristics of the batch samples, performing statistical analysis on the multiple output characteristics, and verifying the optimization effect.
Compared with the prior art, the invention has the beneficial effects that: the method comprises the steps of determining key design parameters of the electromagnetic relay through analysis, carrying out linearity analysis, screening out key nonlinear design parameters, establishing a quality robustness characteristic parameter model of the whole life cycle by using a K-L expansion method and a Krigin method, selecting a proper quality characteristic level according to the quality robustness characteristic requirement level of the whole life cycle of the electromagnetic relay, calculating the quality robustness characteristic deviation degree of the electromagnetic relay under the action of time-varying degradation parameters, and compensating the quality robustness characteristic deviation degree by using a nonlinear parameter combination to realize the optimization requirement of the quality robustness of the whole life cycle.
Drawings
FIG. 1 is a flow chart of a time varying degradation quality feature compensated full life cycle quality robustness optimization method of the present invention;
fig. 2 is a schematic diagram of the compensation effect of the compensation performed in step seven in the embodiment of the present invention.
Detailed Description
The technical solutions in the present invention will be described clearly and completely with reference to the accompanying drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the invention, rather than all embodiments, and all other embodiments obtained by those skilled in the art without any creative work based on the embodiments of the present invention belong to the protection scope of the present invention.
As shown in FIGS. 1-2, the invention discloses a time-varying degradation quality characteristic compensation full life cycle quality robustness optimization method, which comprises
The method comprises the following steps: analyzing the electromagnetic relay according to the working principle, the composition and the manufacturing process parameters of the electromagnetic relay, determining key design parameters (such as armature length, armature thickness, yoke length, yoke thickness and the like) of the electromagnetic relay, and recording the key design parameters as a set of U-U (U-U)1,u2,…,un) Wherein: n is the number of key design parameters;
step two: and (3) carrying out linearity analysis on the key design parameters of the electromagnetic relay analyzed and determined in the step one by using a single parameter or contribution rate analysis method, screening out key nonlinear design parameters from the parameters, and recording the parameters as a set V (V-V)1,v2,…,vm) Wherein: m is the number of non-linear design parameters;
step three: the method comprises the following steps of unfolding the whole life cycle of the electromagnetic relay based on a K-L unfolding method, dividing the whole life cycle into a plurality of sub-life intervals, determining the sub-life intervals to lay a foundation for determining the unity expression of the whole life cycle of the electromagnetic relay, wherein the unfolding process is as follows:
Figure BDA0002313388020000041
wherein: xm(t) and Xr(t) is the life cycle interval process { X (t) epsilon XI(T), T ∈ T } and a radius function, ζj∈ζI=[-1,1]J is 1,2, and K is a standard interval variable which should satisfy
Figure BDA0002313388020000042
λjE [0, ∞) and
Figure BDA0002313388020000051
is composed of
Figure BDA0002313388020000059
According to the Mercer theorem, the characteristic value and the characteristic function of the autocorrelation coefficient can obtain that the correlation coefficient is expressed by the following spectral decomposition:
Figure BDA0002313388020000052
wherein: t, t' are any two moments in the life cycle, the characteristic value lambdajAnd a characteristic function
Figure BDA0002313388020000058
The solution can be solved by Fredholm integration:
Figure BDA0002313388020000053
wherein: the characteristic function should satisfy
Figure BDA0002313388020000054
δijCan be obtained by a Kronecker-delta function;
step four: modeling the quality robustness characteristic parameters of the sub-life intervals by using a Krigin method, carrying out unified expression on each sub-life interval, and establishing an electromagnetic relay full-life cycle quality robustness model, wherein the calculation process is as follows:
the full-life-cycle quality robustness model includes both output characteristic statistics (mean and variance) and a multivariate spatio-temporal model of output characteristic time series relationship, assuming that the full-life-cycle output characteristic is a model including a magnitude of Ns×NtThe matrix of (a) is a regionalized variable in the same time series and space domainWherein
Figure BDA0002313388020000056
Is an Euclidean space dimension plus a time dimension, is a multivariate space-time random field, and:
U(s,t)=[U1(s,t),U2(s,t),K,Up(s,t)]T,p≥3 (4)
wherein: s is as large as NsRepresenting spatial field coordinates, t ∈ NtRepresenting the coordinates of a time field, improving the interpolation precision of a main variable by adding a related covariate,
the interpolation formula of the space-time synergetic kriging is as follows:
Figure BDA0002313388020000057
wherein:
Figure BDA0002313388020000061
is (s, t)0Is output as a characteristic estimate, U11(s,t)1iFor output characteristics being principal variables, U12(s,t)2iIs a co-variable, M1And M2Is the number of interpolation points and M1<M2,ν1iV and v2iThe weight coefficient corresponding to the main variable and the cooperative variable can be determined by a space-time variation functionAnd then the step of determining the number of the first time,
the space-time variation function can be obtained by introducing Lagrange factor
Figure BDA0002313388020000067
And (3) solving:
Figure BDA0002313388020000062
wherein η (g) is a space-time variation function which can be solved by the formula (4),
assuming spatio-temporal random fields
Figure BDA0002313388020000063
Space-time distance d of two positions (d ═ d)s,dt),dsFor vectors that can represent sample spatial distance and direction information, dtThe temporal distance between samples is expressed, so the spatio-temporal variation function can be expressed as:
Figure BDA0002313388020000064
the average value and the variance of the key index output characteristics of the stability of the full life cycle can be obtained through a Kriging model according to the formulas (6) and (7), for the establishment and optimization of the subsequent stable design optimization model of the full life cycle, the converted output characteristics are recorded as f (·), and the average value and the variance are respectively recorded as muf(·)And σf(·)
Step five: selecting a proper quality characteristic level according to the quality robustness characteristic requirement level of the whole life cycle of the electromagnetic relay, wherein the relationship between the quality robustness level and the reliability and the number of the invalid products is shown in the following table;
TABLE I relationship between robust design level and DPMO and reliability
Figure BDA0002313388020000065
Step six: calculating time-varying degradation based on quality robustness level requirements obtained from step five analysisCharacteristic deviation degree of electromagnetic relay quality robustness under parameter action
Figure BDA0002313388020000066
Wherein b is a quality robustness level;
step seven: compensation for quality robust feature bias using nonlinear parameter combinations, i.e., let
Figure BDA0002313388020000071
Thereby compensating the mass shift to the full life cycle mass robustness requirement, the compensation effect is shown in fig. 2;
step eight: and generating batch samples by using Monte Carlo according to the optimized parameter combination of the electromagnetic relay, calculating the output characteristics of the batch samples, performing statistical analysis on the multiple output characteristics, and verifying the optimization effect.
It will be evident to those skilled in the art that the invention is not limited to the details of the foregoing illustrative embodiments, and that the present invention may be embodied in other specific forms without departing from the spirit or essential attributes thereof. The present embodiments are therefore to be considered in all respects as illustrative and not restrictive, the scope of the invention being indicated by the appended claims rather than by the foregoing description, and all changes which come within the meaning and range of equivalency of the claims are therefore intended to be embraced therein. Any reference sign in a claim should not be construed as limiting the claim concerned.
Furthermore, it should be understood that although the present description refers to embodiments, not every embodiment may contain only a single embodiment, and such description is for clarity only, and those skilled in the art should integrate the description, and the embodiments may be combined as appropriate to form other embodiments understood by those skilled in the art.

Claims (3)

1. The method for optimizing the quality robustness of the time-varying degradation quality characteristic compensation in the whole life cycle is characterized by comprising the following steps: the method comprises the following steps:
the method comprises the following steps: analyzing and determining key design parameters of the electromagnetic relay, and recording the key design parameters as a set U-U (U)1,u2,…,un);
Step two: carrying out linearity analysis on the key design parameters of the electromagnetic relay determined by the analysis in the step one, screening out key nonlinear design parameters from the key nonlinear design parameters, and recording the key nonlinear design parameters as a set V (V ═ V)1,v2,…,vm);
Step three: the method comprises the steps that the whole life cycle of the electromagnetic relay is unfolded based on a K-L unfolding method, and the whole life cycle is divided into a plurality of sub-life intervals;
step four: modeling the quality robustness characteristic parameters of the sub-life intervals by using a Kriging method, carrying out unified expression on each sub-life interval, and establishing a quality robustness model of the whole life cycle of the electromagnetic relay;
step five: selecting a proper quality characteristic level according to the quality robustness characteristic requirement level of the whole life cycle of the electromagnetic relay;
step six: calculating the deviation degree of the electromagnetic relay quality robustness characteristic under the action of the time-varying degradation parameter based on the quality robustness level requirement obtained by the analysis in the fifth step
Figure FDA0002313388010000011
Wherein b is a quality robustness level;
step seven: compensation for quality robust feature bias using nonlinear parameter combinations, i.e., let
Figure FDA0002313388010000012
Thereby compensating the mass offset to the full life cycle mass robustness requirement;
step eight: and generating batch samples by using Monte Carlo according to the optimized parameter combination of the electromagnetic relay, calculating the output characteristics of the batch samples, performing statistical analysis on the multiple output characteristics, and verifying the optimization effect.
2. The time-varying degradation quality feature compensated full-life cycle quality robustness optimization method of claim 1, wherein: the unfolding process of the third step is as follows:
Figure FDA0002313388010000013
wherein: xm(t) and Xr(t) is the life cycle interval process { X (t) epsilon XI(T), T ∈ T } and a radius function, ζj∈ζI=[-1,1]J is 1,2, and K is a standard interval variable which should satisfy
Figure FDA0002313388010000021
λjE [0, ∞) and
Figure FDA0002313388010000022
is composed of
Figure FDA0002313388010000023
According to the Mercer theorem, the characteristic value and the characteristic function of the autocorrelation coefficient can obtain that the correlation coefficient is expressed by the following spectral decomposition:
Figure FDA0002313388010000024
wherein: t, t' are any two moments in the life cycle, the characteristic value lambdajAnd a characteristic function
Figure FDA0002313388010000025
The solution can be solved by Fredholm integration:
Figure FDA0002313388010000026
wherein: the characteristic function should satisfy
Figure FDA0002313388010000027
δijFrom KroneThe cker-delta function.
3. The time-varying degradation quality feature compensated full-life cycle quality robustness optimization method of claim 1, wherein: the calculation process of the step four is as follows:
assume that the full life cycle output characteristic is one including a magnitude of Ns×NtThe matrix of (a) is a regionalized variable in the same time series and space domain
Figure FDA0002313388010000028
Wherein
Figure FDA0002313388010000029
Is an Euclidean space dimension plus a time dimension, is a multivariate space-time random field, and:
U(s,t)=[U1(s,t),U2(s,t),K,Up(s,t)]T,p≥3 (4)
wherein: s is as large as NsRepresenting spatial field coordinates, t ∈ NtRepresenting the coordinates of a time field, improving the interpolation precision of a main variable by adding a related covariate,
the interpolation formula of the space-time synergetic kriging is as follows:
Figure FDA00023133880100000210
wherein:
Figure FDA0002313388010000031
is (s, t)0Is output as a characteristic estimate, U11(s,t)1iFor output characteristics being principal variables, U12(s,t)2iIs a co-variable, M1And M2Is the number of interpolation points and M1<M2,ν1iV and v2iThe weight coefficient corresponding to the main variable and the cooperative variable can be determined by a space-time variation function,
the spatio-temporal variation function can be obtained by introducing lagsLangeri factor
Figure FDA0002313388010000032
And (3) solving:
Figure FDA0002313388010000033
wherein η (g) is a space-time variation function which can be solved by the formula (4),
assuming spatio-temporal random fields
Figure FDA0002313388010000034
Space-time distance d of two positions (d ═ d)s,dt),dsFor vectors that can represent sample spatial distance and direction information, dtThe temporal distance between samples is expressed, so the spatio-temporal variation function can be expressed as:
Figure FDA0002313388010000035
the average value and the variance of the key index output characteristics of the stability of the full life cycle can be obtained through a Kriging model according to the formulas (6) and (7), for the establishment and optimization of the subsequent stable design optimization model of the full life cycle, the converted output characteristics are recorded as f (·), and the average value and the variance are respectively recorded as muf(·)And σf(·)
CN201911267915.8A 2019-12-11 2019-12-11 Time-varying degradation quality characteristic compensation full life cycle quality robustness optimization method Active CN111046555B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911267915.8A CN111046555B (en) 2019-12-11 2019-12-11 Time-varying degradation quality characteristic compensation full life cycle quality robustness optimization method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911267915.8A CN111046555B (en) 2019-12-11 2019-12-11 Time-varying degradation quality characteristic compensation full life cycle quality robustness optimization method

Publications (2)

Publication Number Publication Date
CN111046555A true CN111046555A (en) 2020-04-21
CN111046555B CN111046555B (en) 2022-04-08

Family

ID=70235716

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911267915.8A Active CN111046555B (en) 2019-12-11 2019-12-11 Time-varying degradation quality characteristic compensation full life cycle quality robustness optimization method

Country Status (1)

Country Link
CN (1) CN111046555B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115130249A (en) * 2022-08-03 2022-09-30 哈尔滨工业大学 Relay tolerance design method comprehensively considering manufacturing cost and quality loss
CN116415413A (en) * 2023-02-16 2023-07-11 哈尔滨工业大学 Electromagnetic relay parameter robust optimization and reliability design method under uncertain condition

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108229076A (en) * 2018-02-12 2018-06-29 哈尔滨工业大学 A kind of parameters of electromagnetic relay robustness global optimizing method based on response surface
CN108256268A (en) * 2018-02-12 2018-07-06 哈尔滨工业大学 A kind of relay parameter global optimizing method based on K-means radial basis function
CN108256267A (en) * 2018-02-12 2018-07-06 哈尔滨工业大学 A kind of relay quality fluctuation based on radial basis function neural network inhibits design method
CN108268744A (en) * 2018-02-12 2018-07-10 哈尔滨工业大学 A kind of frame circuit breaker Robust Parameter Design method based on Kriging
US20190138926A1 (en) * 2017-12-22 2019-05-09 Beihang University Degradation modeling and lifetime prediction method considering effective shocks

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20190138926A1 (en) * 2017-12-22 2019-05-09 Beihang University Degradation modeling and lifetime prediction method considering effective shocks
CN108229076A (en) * 2018-02-12 2018-06-29 哈尔滨工业大学 A kind of parameters of electromagnetic relay robustness global optimizing method based on response surface
CN108256268A (en) * 2018-02-12 2018-07-06 哈尔滨工业大学 A kind of relay parameter global optimizing method based on K-means radial basis function
CN108256267A (en) * 2018-02-12 2018-07-06 哈尔滨工业大学 A kind of relay quality fluctuation based on radial basis function neural network inhibits design method
CN108268744A (en) * 2018-02-12 2018-07-10 哈尔滨工业大学 A kind of frame circuit breaker Robust Parameter Design method based on Kriging

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
梁慧敏等: "航天继电器电磁系统参数设计方案的稳健性研究", 《低压电器》 *
邓杰等: "直动式继电器电磁系统稳健性参数优化设计", 《电器与能效管理技术》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115130249A (en) * 2022-08-03 2022-09-30 哈尔滨工业大学 Relay tolerance design method comprehensively considering manufacturing cost and quality loss
CN115130249B (en) * 2022-08-03 2024-04-12 哈尔滨工业大学 Relay tolerance design method comprehensively considering manufacturing cost and quality loss
CN116415413A (en) * 2023-02-16 2023-07-11 哈尔滨工业大学 Electromagnetic relay parameter robust optimization and reliability design method under uncertain condition
CN116415413B (en) * 2023-02-16 2023-10-03 哈尔滨工业大学 Electromagnetic relay parameter robust optimization and reliability design method under uncertain condition

Also Published As

Publication number Publication date
CN111046555B (en) 2022-04-08

Similar Documents

Publication Publication Date Title
Genton et al. On the likelihood function of Gaussian max-stable processes
CN111046555B (en) Time-varying degradation quality characteristic compensation full life cycle quality robustness optimization method
Zhang et al. On a nonautonomous SEIRS model in epidemiology
CN110929421B (en) Relay optimization method comprehensively considering permanent magnet degradation and manufacturing uncertainty
Lapenta et al. Particle acceleration and energy conservation in particle in cell simulations
Gao et al. Analysis of gene expression programming for approximation in engineering design
CN112446110A (en) Application method of EOASM algorithm in proxy model construction of robot palletizer driving arm seat
Liu et al. Optimal parameter selection for nonlinear multistage systems with time-delays
CN115470704A (en) Dynamic multi-objective optimization method, device, equipment and computer readable medium
Chen et al. A comparative study of production efficiency in coastal region and non-coastal region in Mainland China: An application of metafrontier model
Zhang et al. Fractional Brownian fields for response surface metamodeling
Xie et al. A novel mutual fractional grey bernoulli model with differential evolution algorithm and its application in education investment forecasting in China
Fang et al. Structural reliability analysis using response surface method with improved genetic algorithm
Wang et al. Recursive quadratic filtering for linear discrete non-Gaussian systems over time-correlated fading channels
Bulut et al. Comparing energy demand estimation using various statistical methods: the case of Turkey
CN111259340B (en) Saturation load prediction method based on logistic regression
CN106647250B (en) Based on offline optimization/online lookup table mode double-layer structure forecast Control Algorithm
CN111158241A (en) Time-lag correlation H-infinity control method of linear singular system with uncertain time lag
CN115130249B (en) Relay tolerance design method comprehensively considering manufacturing cost and quality loss
Zhou et al. Application of the three-parameter discrete direct grey model to forecast China’s natural gas consumption
CN115907177A (en) Urban carbon emission prediction method and system
CN113361176B (en) Nonlinear characteristic value topology optimization method and system considering frequency-dependent material
Wang et al. Asymptotic properties of a stochastic chemostat including species death rate
Liu et al. A novel supervised multi-model modeling method based on k-means clustering
CN111123696B (en) Redundant channel-based networked industrial control system state estimation method

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