Next Article in Journal
Global Optimization of UAV Area Coverage Path Planning Based on Good Point Set and Genetic Algorithm
Next Article in Special Issue
Design and Simulation Analysis of an Electromagnetic Damper for Reducing Shimmy in Electrically Actuated Nose Wheel Steering Systems
Previous Article in Journal
Numerical Simulation of Aircraft Icing under Local Thermal Protection State
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Reliability-Oriented Configuration Optimization of More Electrical Control Systems

1
School of Automation Science and Electrical Engineering, Beihang University, Beijing 100191, China
2
Beijing Advanced Innovation Center for Big-Data Based Precision Medicine, Beihang University, Beijing 100191, China
3
Shenyuan Honors College, Beihang University, Beijing 100191, China
*
Author to whom correspondence should be addressed.
Aerospace 2022, 9(2), 85; https://doi.org/10.3390/aerospace9020085
Submission received: 24 November 2021 / Revised: 20 January 2022 / Accepted: 28 January 2022 / Published: 6 February 2022

Abstract

:
More electrical vehicles adopt dissimilar redundant control systems with dissimilar power supplies and dissimilar actuators to achieve high reliability and safety, but this introduces more intricacy into the configuration design. Currently, it is difficult to identify the optimum configuration via the conventional trial-and-error approach within an acceptable timeframe. Hence, it is imperative to discover novel methods for the configuration design of more electrical vehicles. This paper introduced the design specification of more electric vehicles and investigated the contribution of different kinds of actuators, presenting a new multi-objective configuration optimization approach on the foundation of system reliability, weight, power, and cost. By adopting the non-dominated sorting genetic algorithm-II (NSGA-II), the Pareto optimization design set was obtained. Then, the analytic hierarchy process (AHP) was introduced to make a comprehensive decision on the schemes in the Pareto set and determine the optimal system configuration. Eventually, numerical results indicated that the reliability of our designed configuration increased by 5.89% and 55.34%, respectively, compared with dual redundancies and single redundancy configurations, which verified the effectiveness and practicability of the proposed method.

1. Introduction

Safety critical systems in aircraft [1], carrier rockets [2], ships [3] and large machines [4] usually adopt redundant actuation systems to guarantee high reliability and safety [5]. In such designs, the rest systems can complete the work when one or more actuators fail. To reduce the number of common-cause faults and common-mode faults, dissimilar redundancy systems are adopted consequentially [6], in which different actuators with the same performances are used. With the rapid development of more electrical technologies, and more electric vehicles based on different power supply systems, actuation systems and computers have become increasingly common. In terms of actuation systems, hydraulic actuators (HAs), electro-hydrostatic actuators (EHAs) and electromechanical actuators (EMAs) [7] are leveraged in more electric vehicles. For instance, more electrical aircrafts are adopting heterogeneous actuation systems to move control surfaces based on HAs and EHAs [8]. Since HAs utilize hydraulic power supply systems while EHAs and EMAs adopt electrical power supply systems, the configuration design of more electric vehicles (MEVs) faces substantial challenges. The first one is how to design more electric control systems to maintain high mission reliability without permissive cost and weight. The second one is how to strike a balance between power supply systems and actuator distribution. The last one is how to determine the types of actuators for more electric control systems, which include HAs, EHAs, and EMAs. Given the aforementioned factors, the design of safety-critical system configurations becomes a daunting task. Therefore, an appropriate and effective approach should be adopted to solve the combination-explosion problem of system configurations, thus achieving the optimization design of a more electric control system (MECS) [9,10].
The optimization design of an electric control system needs to realize multi-objective optimization and comprehensively consider various contradictory indicators. Traditional methods include the weighting method [11,12], the constraint method [13,14], and goal programming [15,16,17]. However, these methods require the designer to master corresponding background knowledge to determine weights, energy dissipation, and other indicators, and the explosive configuration is another challenge that arises from those methods. In addition, the program needs to be run independently many times, which may lead to inconsistent results each time and render it more difficult for designers to make eventual decisions. In recent decades, various multi-objective intelligent optimization algorithms have emerged and vigorously developed, for example, the genetic algorithm (GA) [18], differential evolution algorithm (DE) [19], particle swarm optimization algorithm (PSO) [20,21], and non-dominated sorting genetic algorithm (NSGA) [22]. Among them, NSGA is famous for its good performance due to the adoption of non-dominated sorting and fitness-sharing strategy. The former increases the possibility of superior characteristics inherited by the next generation, while the latter sustains population diversity, overcoming the overmuch reproduction of super individuals and preventing prematurity convergence [23]. Nevertheless, there are also disadvantages such as high computational complexity in NSGA. Thus, in 2002, Deb et al. proposed NSGA-II [24]. In contrast to NSGA, the rapid non-dominated sorting method, elitism preservation strategy and congestion comparison were introduced in NSGA-II, which greatly reduced the computational complexity of the algorithm. Moreover, NSGA-II also expanded the distribution space of the solution set in the Pareto frontier, thus maintaining population diversity [25]. At present, NSGA-II has been widely applied to tackle multi-objective optimization problems in various fields. Xia et al. [26] investigated the multi-objective optimization problem for AUV conceptual design, where NSGA-II was applied to find the optimal Pareto frontier, with a comprehensive consideration of cost, effectiveness, and risks. The result verified the effectiveness of the algorithm. Alam et al. [27] studied the problem of AUV design and construction and employed NSGA-II to determine the optimum design of a torpedo-like AUV with a total length of 1.3 m. In addition, a kind of heavier-than-water underwater vehicle (HUV) was regarded as a research object by Liu et al. [28], and NSGA-II was adopted to establish a global approximation model, thus assisting the eventual optimal design of HUV. Nevertheless, these methods ignored the optimization designs of the global architectures of dissimilar control systems. As mentioned above, the optimization design of MECS is to solve the multi-objective optimization issue, where some contradictory indicators require all-round considerations. Therefore, the NSGA-II, combined with the AHP [29,30,31,32] method, was adopted in this paper to solve this problem, in which the binary encoding mode was employed and the decision variables were confined in discrete space. Moreover, the mutation process was modified to better match the characteristics of the problem. Eventually, the case study verified the validity of our approach.
The other sections of the paper are arranged as follows: Section 2 demonstrates the system descriptions. Section 3 presents the optimization approach on the foundation of NSGA-II and AHP method. Section 4 displays the case study, while the conclusions are presented in Section 5.

2. Mathematical Modeling of More Electric Control System (MECS)

According to the design specification of safety critical systems, more electric control systems require high reliability and safety. For example, the possibility of mission failure per flight owing to certain material damages in the flight control system should not surpass the upper limit. Normally, the failure rate of a flight control system is 10 9 / h to 10 10 / h for a commercial aircraft. Hence, it is necessary to utilize a dissimilar flight control system to maintain high reliability.

2.1. Single Control System Structure

A schematic diagram of the basic structure of a more electric control system, in which the controller, actuator, and sensor are main components, is illustrated in Figure 1. Since the electric power supply and hydraulic power supply are provided simultaneously, the actuator can be selected from HAs, EHAs, or EMAs, as shown in Figure 1. A HA is powered by a hydraulic supply system while an EHA or EMA is powered by an electric supply system. Since a centralized hydraulic power supply system has the characteristics of high-power density and fast response, HAs are widely used in control systems. However, the control system is heavy and easily exposed to oil contamination when HAs are adopted as they need pipeline transmission from the centralized power supply to the actuator. An EHA consists of the motor, pump, and cylinder that replaces the pipeline transmission by wire transmission. Hence, the weight of an EHA is light, whereas its heat dissipation is poor. An EMA consists of the motor and the ball screws. The command is transmitted from the wire controlling the motor and ball screws to drive the load. Therefore, an EMA is simple and light weight, but it gets stuck easy. Hence, to determine which type of actuator should be used, one has to comprehensively consider the weight, performance, and reliability.

2.2. Redundant Configuration of Actuation System

The more electric control of safety critical systems not only requires the high-precision performance but also requires extremely high reliability and safety. Thus, redundant designs are often adopted, such as dual-redundant actuators and triple-redundant actuators. In order to integrate the advantages of hydraulic power and electric power, the typical isomorphic actuation system and heterogeneous actuation system shown in Figure 2, Figure 3 and Figure 4 are used in real-world applications. In an aircraft, the HA/HA, HA/EHA, and HA/EMA are classic dual-redundant actuation systems, as shown in Figure 2, Figure 3 and Figure 4.
Generally, the regulation of a dissimilar redundant actuator is to design the different actuator powered by a different power supply for one surface. In such design, when any actuator or power supply system fails, the other can fulfill the task through fault switching. Table 1 presents the typical configuration of a dissimilar redundant actuation system.
Remark 1.
Triple redundancy was considered in this research because when the number of redundancies increases, the performance indicators of an MECS will increase accordingly, which is not beneficial for system operation. Besides, an actuation system with more than three redundancies is rarely used in practice.

2.3. Redundant Configuration of More Electric Control System

As mentioned above, dissimilar redundancy technology is widely used in safety critical system design in order to improve system reliability. In a commercial aircraft, various controllers, actuators, and power supplies are applied in the aircraft actuation system. A typical redundant actuation system configuration based on high reliability for a commercial aircraft is illustrated in Figure 5. Although multiple control computers, actuators, and power supply systems can achieve high reliability, weights and costs will increase correspondingly, with maintainability decreased as well. Therefore, it is imperative to optimize the quantity of redundant control computers, power supply systems, and actuators. At same time, designers have to solve the common faults of dissimilar power supplies and dissimilar actuators as shown in Figure 2, Figure 3 and Figure 4.
In Figure 5, I expresses the information nodes of the control computer; P describes the power supply nodes; A is the actuator nodes; and S shows the surface of the aircraft. Herein, the actuator includes an HA, EHA, and EMA, as shown in Figure 2, Figure 3 and Figure 4. Both the fonts of actuation nodes and rudder surfaces are bold since they are our main considerations.
Definition 1.
Power supply module set is  P = { P i } ( i = 1 , 2 , , n ) . For aircraft shown in Figure 5,  P 1 = H means the power supply;  P 1 is the hydraulic power; and  P 2 = E means that the power supply  P 2 is electric power.
Definition 2.
Actuator module set is  A = { A j } ( j = 1 , 2 , , m ) . For aircraft shown in Figure 5,  A 1 = H A means that the actuator  A 1 is a hydraulic actuator.  A 2 = E H A means that the actuator  A 2 is an electro-hydrostatic actuator, and  A 3 = E M A means that the actuator A 3 is a mechatronic actuator.
Remark 2.
Here, the commercial aircraft redundant actuation system is introduced as a typical MESC merely. In fact, our method is extensible and can be applied to more kinds of systems. In Section 4, we abstract a typical MECS with 3 hydraulic power supplies and 2 electric power supplies as the research object, which can be also extended to more combinations.

3. Multi-Objective Optimization of MECS Based on NSGA-II and AHP

Though an MECS can improve system reliability and ensure safety, when the system structure is more sophisticated, the corresponding structural indicators such as weight and cost will increase enormously. The design of an MECS may also increase the system volume and faults derived from the heavy weight. Therefore, we need to optimize factors such as the weight, power dissipation, and reliability of an MECS to ensure a compromise between these indicators.

3.1. Multi-Objective Optimization Modeling of MECS

The primary task for optimization design is to transfer the optimization problem into mathematical descriptions. Generally, an optimization problem is composed of three elements: objective functions, decision variables, and constraint conditions. The optimization problem of an MECS can be expressed as
J = max R s ( t ) R s ( t ) = f ( H i , E i , H A j , E H A j , E M A j ) s . t . 0 m ( H i , E i , H A j , E H A j , E M A j ) M 0 ψ ( H i , E i , H A j , E H A j , E M A j ) Ψ 0 c ( H i , E i , H A j , E H A j , E M A j ) C 0 g ( H i , E i , H A j , E H A j , E M A j ) G i = = { 1 , , m } j = = { 1 , , n }
where function J is the maximum reliability; m ( H i , E i , H A j , E H A j , E M A j ) is the actual evaluation mass of the MECS; M is the superior limit of mass; ψ ( H i , E i , H A j , E H A j , E M A j ) expresses the actual evaluation power dissipation of the MECS; Ψ is the superior limit of power dissipation; c ( H i , E i , H A j , E H A j , E M A j ) describes the actual evaluation cost of the MECS; c represents the superior limit of cost; g ( H i , E i , H A j , E H A j , E M A j ) represents the other actual evaluation indicators of MECS; and G represents the superior limit of other indicators. The specific evaluation methods of objective functions and constraints are given below.

3.1.1. Objective Function

According to Equation (1), reliability is the most essential indicator and the objective function. System reliability refers to the ability of the system to meet the specified functions within the specified time and under the specified conditions. Only when specific functions of an MECS are controlled effectively at the same time can the overall function of the system be guaranteed. Therefore, the overall functional reliability of an MECS is defined as the ability to effectively control the specific functions simultaneously within the specified time and under the specified conditions.
Based on the definition above and by considering an MECS with eight actuation functions, the functional reliability of an MECS can be expressed as
R S = Pr F 1 F 2 F 8 = Pr S
where F 1 , , F 8 represent the pivotal functions that MECS should accomplish.
By considering the multiple control surface combinations involved in the realization of system functions, Equation (2) can be rewritten as
R S = Pr x = 1 p F 1 y = 1 q F 2 z = 1 m F 8
where x = 1 p F 1 y = 1 q F 2 z = 1 m F 8 represents a certain combination that realizes the pivotal functions.
To simplify the calculation, the program first calculated the overall reliability of electric power supplies E 1 , E 2 and then calculated the overall reliability of hydraulic power supplies H 1 , H 2 . Next, we calculated the minimum path and the program performing disjoint operation. Through decoupling analysis, nine minimal paths were obtained. Consequently, the functional reliability of an MECS can be calculated by
R S = Pr i = 1 9 S i = q E 1 q E 2 q E 5 q E 6 p E B p E D p E H p E J p E N p E O p E P p E Q p E R p E S p E T p E U p E V p E W p E X p H 1 p H 3 p P 1 p P 2 p K 1 p K 2 p D D p Y Y p C 1 p C 2 p C 3 + q E 1 p E 2 q E 5 q E 6 p E B p E D p E H p E J p E N p E O p E P p E Q p E R p E S p E T p E U p E V p E W p E X p H 3 p P 1 p P 2 p K 1 p K 2 p D D p Y Y p C 1 p C 2 p C 3 + + q E 1 p E 2 p E 5 p E B p E D p E H p E J p E N p E O p E P p E Q p E R p E S p E T p E U p E V p E W p E X p H 1 p H 3 p P 1 p P 2 p K 1 p K 2 p D D p C 1 p C 2 p C 3 + p E 1 p E 5 p E B p E D p E H p E J p E N p E O p E P p E Q p E R p E S p E T p E U p E V p E W p E X p H 1 p H 3 p P 1 p P 2 p K 1 p K 2 p D D p C 1 p C 2 p C 3
Remark 3.
Equation (4) gives the relationship between system reliability and component system. In fact, system reliability is determined by each component reliability, that is, the actuator reliability, while actuator reliability is associated with weight, power consumption, and cost. These indicators will eventually have an influence on system reliability.

3.1.2. Constraint Conditions

The aim of optimization design is to identify the optimal solution from the practical solutions. The optimal solution can meet the goal of design as far as possible. The elevation of weight will restrict the MECS installation power; hence, weight is a vital index for system property. Power dissipation is another vital property index. In the optimization design of an MECS, particularly for a highly-efficient MECS, how to decrease the power dissipation with the existent technology and apparatus is a challenge which need to be tackled by designers. Meanwhile, as an MECS has to harbor high performance and reliability, so the reliability of an MECS should be considered. For practical application, the cost is an indispensable indicator that must be evaluated. Eventually, weight, power dissipation, cost, and reliability are determined as constraint conditions.
  • Weight
Weight is a vital evaluation index and decisive factor of an MECS. When we evaluate weight, the difficulties usually include various components, many of which exhibit a nonlinear growth relationship with design requirements, and serious coupling with other systems with different types of actuators. For an EHA, since it consists of the integration block, cylinder, pump, and motor, its weight evaluation can be presented as
M E H A = M b l o c k + M c y l i n d e r + M p u m p + M m o t o r
where M b l o c k , M c y l i n d e r , M p u m p refer to the weight of the integration block, cylinder, and pump, respectively. M m o t o r denotes the weight of the brushless direct-current motor (BLDC). All of them have the same unit, kilograms.
A permanent magnet brushless direct-current motor (BLDC) is utilized in the EHA due to its satisfactory control property and great power-to-mass ratio. The weight forecast of the BLDC is expressed by
M m o t o r = 0.628 T 3 / 3.5 + 0.783
in which T represents the torque of BLDC.
The pump weight is proportionate to pump output. Thus, M p u m p can be expressed as [33]
M p u m p = 0.339 D + 2.038
The integration block is the EHA frame and involves the indispensable parts such as the nonreturn valve, filter, and accumulator. The overall weight is speculated by the EHA power, which is expressed as [33]
M b l o c k = 0.105 P E H A + 2
where P E H A denotes the maximal power of EHA.
The fluid cylinder weight includes the four parts stated below:
M c y l i n d e r = M c o v e r + M s h e l l + M p i s t o n + M r o d
where M c o v e r , M s h e l l , M p i s t o n , M r o d are the weights of the cylinder cover, shell, piston, and plunger rod, respectively. All of them can be calculated by
M r o d   = π 4 × d r o d   2 × L r o d   × ρ s t e e l   M p i s t o n   = A × t p i s t o n   × ρ c o p p e r   M s h e l l   = π 4 × d s h e l l   2 4 A π × L s h e l l   × ρ s t e e l   M c o v e r   = 2 × π 4 × d s h e l l   2 × t c o v e r   × ρ s t e e l  
Thus, the total weight can be evaluated as
M = i = 1 n H A M H A + i = 1 n E H A M E H A + i = 1 n E M A M E M A + M e l e c + M p i p e + M w i r e
where M H A , M E H A , M E M A represent the weights of HA, EHA, and EMA, respectively, and M e l e c , M p i p e , M w i r e represent the weights of electric equipment, pipe and wire, respectively.
  • Power efficiency
For an EHA, the motor converts electric energy to hydraulic power, pushing the cylinder to realize the movement of the control surface. For an EHA system, when its control surface load is certain, it saves more energy, resulting in lower power dissipation. The control surface torque T S , control surface velocity ω S , and pressure P S are known design specifications. The aim of decreasing the power dissipation of an EHA is to decrease the power output of the motor P , which is written as
P = T × W
The motor is directly connected to the pump, so the torque and speed of the motor and pump are the same. The torque of the pump is expressed as
T = J p m × d ω d t + K f r i c × ω + ( p f + 2 p p i p e ) × D 2 π
where J p m is the pump rotational inertia; K f r i c denotes the viscosity factor, J p m × d ω / d t = 0 ; the pump rotation velocity is constant; p f denotes the differential pressure between the 2 cylinder cavities; and p p i p e denotes the pressure consumption within the pipe. Hence, Equation (13) is written as
T = K f r i c × ω + ( p f + 2 p p i p e ) × D 2 π
in which p f = F / A , F is the loading force of an EHA and A denotes the piston effective area, which can be calculated by
A = k F P S
in which k > 1 denotes a logical excess margin, and F can be described by
F = T S × sin ( θ + 30 ° ) L
For an EMA, the power dissipation is determined by the ball screw. Therefore, Equation (17) gives the calculation method for the drive torque of the ball screw:
T a = F a × I 2 π × n 1
where T a represents the drive torque; I is the lead screw; n 1 is the positive efficiency of the feed screw; and F a refers to the axial load, which can be expressed by
F a = F + μ m g
where F means the axial cutting force of the lead screw, μ is the comprehensive friction coefficient, and m refers the weight of the worktable and workpiece.
After the calculation of the drive torque T a , the motor power can be determined accordingly; thus, the corresponding power dissipation of the motor can be obtained.
For the power dissipation evaluation of a HA, the volumetric efficiency of the pump η 0 is the main affecting factor, which can be described by
η 0 = Q Q 0 × 100 %
where Q 0 means the theoretical flow of the pump, and Q is the actual pump flow, which is written as
Q = Q p u m p + Q p i p e + Q c y l i n d e r
where Q p u m p denotes the pump leak, Q p i p e denotes the loss of flow within the hydraulic tube, and Q c y l i n d e r denotes the cylinder flow.
Q p u m p = ξ × ( p 1 p 2 ) = ξ × F / A Q c y l i n d e r = A × v Q p i p e = 2 × ξ × p p i p e
where ξ is the pump leakance coefficient, p p i p e denotes the pipe pressure drop, and v denotes the cylinder speed.
  • Cost
For an MECS, cost is another index worthy of consideration. In this paper, the total costs are divided into two parts:
C = C m a n u + C o p e r
where C represents the total cost of an MECS, C m a n u refers to the manufacturing cost, and C o p e r represents the operation cost.
However, when we evaluated the total cost, the component cost exhibited a nonlinear growth relationship with the requirement. To simplify this process, this paper utilized different configurations to obtain the costs of all components. In the evaluation of manufacturing costs, the main components considered included the HA, EHA, EMA, oil boxes, pipelines, cables, engine-driven pump (EDP), and electric motor driven pump (EMP). The evaluation was conducted mainly based on the similarity principle. Thus, the manufacturing costs is expressed as
C m a n u = i = 1 n e d p C e d p + i = 1 n e m p C e m p + i = 1 n t a n k C t a n k + i = 1 n m o t o r C m o t o r + i = 1 n a c t C a c t + i = 1 n p i p e C p i p e + i = 1 n w i r e C w i r e
where C e d p represents the cost of the EDP; C e m p refers to the cost of the EMP; C t a n k is the cost of the tank; C m o t o r is the cost of the motor; C a c t describes the cost of the actuator; C p i p e expresses the cost of the pipe; and C w i r e denotes the wire cost.
The operation expense primarily denotes the Direct Operating Costs (DOCs), which are changeable costs directly derived from operating the aircraft [34]. The fuel expense and maintenance expenditure of DOCs are changeable. Hence, the DOCs are predominantly related to fuel expense and maintenance expenditure herein. The operation expense speculation is written as
C o p e r = C o i l + C m a i n
where C o i l and C m a i n represent the fuel cost and maintenance cost, respectively.
Based on the statistical proportion of main flight control maintenance cost in the total maintenance expenditure, the fuel expense speculation primarily considers the fuel consumption during the life of the commercial aircraft. The function of fuel cost with respect to time is expressed as
C o i l ( t ) = 50 × C 0 × W e i g h t × ( 1.02 t 1 )
where C 0 means the average annual fuel cost of the aircraft and W e i g h t represents the fuel weight.
The maintenance cost of the aircraft in the whole life cycle mainly has two parts [35]. One is called the cyclic cost C C , which is related to the take-off and landing of the aircraft, such as the maintenance of its braking device, flap system, and landing gear. The other cost is related to the flight time of the aircraft, such as the regular loss and replacement of parts caused by each hour of flight. This part of the cost is called hourly cost C H . The evaluation of the maintenance cost is usually accomplished through the statistic of man-hour cost. Equation (26) provides the method for man-hour cost calculation
C m a i n = C H × t + C C × t + M × ( 1 / α ) × R × t

3.1.3. Design Variables

The expression of a goal suggests that there are substantial that which ought to be identified in configuration designs. Nevertheless, considering all those variables during optimization makes convergence difficult and is quite time-consuming. Hence, merely vital variables having a remarkable influence on MECS property ought to be optimized. Overall, the selective principles of the design variates are stated below:
  • The quantity of design variates ought to be decreased to the utmost extent. Overall, the quantity of design variates in mechanical optimization design should not surpass 5.
  • The variables ought to exert a remarkable impact on the goal function. Indexes affecting the constraint and property directly ought to be chosen as design variates.
  • The chosen variates ought to be independent.
  • The variates ought to be chosen as per the optimization goal.
According to the aforesaid principles, the optimization variates are chosen via:
x = x 1 x 2 x 3 = n H A n E H A n E M A

3.2. Multi-Objective Optimization Based on NSGA-II

The optimization issue herein aimed at the configuration design of an MECS. NSGA-II is proper for optimizing multi-objective issues as it exhibits strong distributed capability and rapid convergence. The flow chart of NSGA-II is shown in Figure 6. Subsequently, the encoding/decoding mode and fast non-dominated sorting process are elaborated.

3.2.1. Encoding and Decoding

In our research, the binary encoding mode was adopted in NSGA-II, which was simple, efficient, and easy to design and use for implementing crossover and mutation operations. Two types of energy information were encoded through this mode, while the configuration of each actuation system was also expressed in binary. Figure 7 illustrates this coding mode.

3.2.2. Fast Non-Dominated Sorting

Compared with non-dominated sorting approach utilized in NSGA arithmetic, the rapid non-dominated sorting utilized in NSGA-II algorithm requires lower computational complexity, which is only about O ( M N 2 ) . The overall algorithm flow was presented as follow:
Algorithm 1: Fast non-dominated sorting
1:    Fast-non-dominated-sort ( P )
2:     for each  p P
3:      for each  q P
4:       if  p n q , then     # if p is dominated by q , then add q to S p
5:         S p = S p q
6:       else if  q n p , then
7:         n p = n p + 1
8:      if  n p = 0 , then
9:         p r a n k = 1 , F 1 = F 1 P     # when n p of the individual is 0, then this individual is the first level of Pareto
10:    The comparison of dominating relationships between individuals, S p and n p , are introduced for storage and records, respectively; n represents the comparison of dominating relationships. The solution of n p = 0 is stored in the records of level 1, and the solution of level 1 has higher priority than that of level 2.
11:    i = 1.
12:    while  F i =  do
13:    H=
14:    for each p F i
15:     for each q S p         # Sort all the individuals in S p
16:       n q = n q 1
17:       if  n q = 0 , then       # when n q of the individual is 0, it is a non-dominated individual
18:          q r a n k = i + 2 , H = H q # The Pareto level of this individual is the current highest level plus 1. At this moment, the initial value of i was 0, so we added 2.
19:    end while
20:     i = i + 1 , F i = H
21:    Loop the program to obtain level 2, level 3… The computational complexity is O ( M N 2 )

3.2.3. Crowding Degree Calculation

Crowing degree I d refers to the local crowding distance between any two adjacent points whose level are the same in the target space. The purpose of introducing crowding degree is to improve the distribution uniformity of the Pareto solution set. Furthermore, it could not only boost the diversity of population but also enhance system robustness. The crowding degree I d could be expressed as the length of the maximum rectangle of individual i on both sides, where the rectangle only includes i itself.
The main procedure to determine the degree of individual congestion involves three steps:
  • Define the crowding degree of every individual in population i as I d = 0 ;
  • Define the crowding degree 0 d and i d of boundary individuals as according to each evaluation indicator;
  • Define the crowding degree of marginal individuals as a larger number L ( 1 ) d = L ( l ) d = M to prioritize individuals on the sorting edge; thus, the crowding degree of any other individual I d can be expressed as
I d = j = 1 m f j i + 1 f j i 1
where j refers to each evaluation indicator; f j i + 1 represents the j t h evaluation indicator value of individual i + 1 ; and f j i + 1 represents the j t h evaluation indicator value of individual i 1 .
Through the aforesaid steps and calculation, every individual was endowed with two attributes, i.e., the crowding distance and the rank. This laid the foundation for the follow-up processes.

3.2.4. Optimal Selection

The optimal selection avoids the loss of effective factors and ensures the survival rate of high-performance individuals, so it can ensure that the Pareto optimal solution is continuously optimized. The optimal selection can not only improve the efficiency and convergence of the optimization but also ensure the uniformity of the optimization process. The selection process is completed by comparing the results of fast non-dominated sorting and crowding calculation, and better individuals are selected after comparison. Two steps are involved during this process.
  • The first step is the rank comparison. Select two individuals a and b randomly and make comparison between A r a n k (the non-dominated rank of individual a ) and B r a n k (the non-dominated rank of individual b ). When A r a n k < B r a n k , a is better than b and vice versa. Moreover, the crowding degree requires to be compared when A r a n k = B r a n k ;
  • The second step is the crowding degree comparison. When condition A d > B d is satisfied, it indicates that individual a is better; otherwise, individual b is better. Then, the better individual is selected to continue the following optimal processes.

3.2.5. Crossover

The main function of crossover is to simulate gene recombination in the process of heredity and evolution. There are various approaches to implement it, such as uniform crossover, multi-point crossover, and binary crossover. In our research, the simulated binary crossover is adopted for crossover operation, which is expressed as
y 1 j = 0.5 × 1 + γ j x 1 j + 1 γ j x 2 j y 2 j = 0.5 × 1 γ j x 1 j + 1 + γ j x 2 j 2 u j 1 η + 1 i f   u j 0.5 1 2 1 u j 1 η + 1 e l s e
where u j 0 , 1 ; x 1 j and x 2 j are parent individuals; y 1 j and y 2 j are offspring individuals; and η > 0 refers to the cross-distribution index. Generally, η = 20 is the best value of the cross-distribution index by default.

3.2.6. Mutation

Mutation is widely applied to simulate variation links in the process of biological heredity and evolution, which refers to the replacement between genes, that is, substitute the other gene values with the alleles on the locus to create new individuals. Through mutation, not only could the local searchability be improved but also the population diversity is guaranteed.
This paper introduced the polynomial mutation method into the research on multi-objective optimization problems. The form of mutation operator is expressed as
V k = V k + δ u k L k   δ = 2 u + ( 1 2 u ) 1 δ 1 η m + 1 1 η m + 1 1   i f   u 0.5 1 2 ( 1 u ) + 2 ( u 0.5 ) 1 δ 2 η m + 1 1 η m + 1 e l s e   i f   u > 0.5
where δ 1 = V k L k u k L k ,   δ 2 = u k V k u k L k , u denotes a stochastic number in interval [0,1], η m denotes the distribution index, and V k is a parent individual.

3.3. Comprehensive Evaluation of System Configuration Based on AHP

The analytic hierarchy process (AHP) is a multi-objective decision analysis approach that combines qualitation and quantitation analyses, where elements associated with task decisions are divided into the object level, criterion level and scheme level. The AHP mathematizes the decision-making via a few quantitation data based on deep analyses, such as the influential factors and inner association of the decision-making problems. AHP can offer an easy decision approach for intricate decision issues with several standards and no evident structure features. Questions to be evaluated in MECS optimization design ought to be methodized and layered, and the structure model of hierarchy analysis ought to be constructed. A three-level hierarchy structure of AHP was presented by Figure 8. Level 1 is the general objective of decision-making, which denotes the optimum design of MECS. Level 2 denotes the criterion layer, where the weight, power dissipation, expense, and dependability are utilized as evaluation standards of system multi-objective optimization design. Level 3 denotes the scheme layer, where the solutions are determined in the Pareto frontier acquired by NSGA-II. Besides, the process of AHP mainly includes the judgment matrix structure, weight coefficient calculation, consistency of judgment matrix verification, and weight coefficients of goal level calculation. The flow chart of AHP was shown in Figure 9.
  • Construct the decision-making model for AHP according to Figure 8.
  • Structure the judgement matrix J . Judgement matrix is established as per the association between the goals in the criterion layer.
J = 1 3 9 1 1 3 1 1 3 1 3 1 9 1 3 1 1 9 1 3 9 1
3.
Validate the judgement matrix coherence. Coherence index is computed via:
C I = λ max n n 1 C R = C I R I
where C I denotes the coherence index of the judgment matrix, R I denotes the average stochastic coherence index of the matrix (specific values were presented by Table 2), C R denotes the stochastic coherence ratio of the matrix, λ max denotes the maximal value of the matrix characteristic value, and n denotes the matrix order.
By computation, if λ max = 3.7638 and CI = (3.7638 − 4)/(4 − 1) = −0.0787, then CR = −0.0787/0.90 = −0.0875 < 0.1, which suggests that the weight matrix is coherent. When the consistence test is not met, return to step 2 and reconstruct the judgement matrix.
4.
Compute the weighted coefficient between the contrasted elements with the relevant standards. Compute the continued product M i of each row element in A , the product of every row element, and its n-th root w ¯ 1 .
     M i = i = 1 n x i j = 27 0.0370 0.0041 27 T , i = 1 , , n w ¯ 1 = M i n = 2.2795 0.4387 0.2533 2.2795 T , i = 1 , , n
Normalize w ¯ 1 as W i = w ¯ 1 / i = 1 n w ¯ 1 , where W i denotes the weighted coefficient of every factor.
W = 0.4341 0.0835 0.0482 0.4341
5.
Speculate the design in the Pareto frontier as per the weighted coefficients of every standard, and afterwards get the optimum design of MECS.

4. Case Study and Discussion

The diagram of an MECS with three fluid power supplies and two electric power supplies (3H2E) was presented by Figure 10. The tendency toward more electricity causes the transformation from hydraulic power into standby energy, or the replacement of hydraulic power with electrical energy, which produces more options for all control surfaces. For every actuator, it can be an HA linked to hydraulic source or an EHA linked to electrical source. In addition, each actuator is controlled by at least one controller. Therefore, the number of alternative configurations of the MECS is exceedingly high.
In this section, we introduced the proposed optimization method into the case study of the MECS illustrated in Figure 10. By setting the maximum optimization iterations of NSGA-II as 200, and population size as 500, we obtained the Pareto frontier of the optimization result illustrated in Figure 11. W refers to the weight of the MECS, P is the power consumption of the MECS, and C represents the expense of the MECS. The red-round-dot-curved surface is a 3D Pareto optimum front curved surface. Every dot denotes an actuation design. The green line, marked by hollow pentagram; blue line, marked by cross symbol; and yellow line, marked by triangle were projected on three planes. The Pareto front was utilized to realize multi-objective optimization and to achieve a certain decrease in configuration quantity.
Moreover, the graph in Figure 11 presents the relationship among the three objectives. The line in x, y plane presents that weight and power have the relation of promotion. The line in x, z plane presents that weight and expense are contradictory. The line in y, z plane displays contradictory power and expense of an MECS. When the power decreases, the expense is elevated. The design objective is to minimize power and expense; hence, those two functions conflict with one another. Due to the conflicts between design indexes, it is hard to make each index optimum simultaneously.
In addition, to better illustrate the relationship between these three objectives and system reliability, we projected Figure 11 onto three 2-D planes, respectively, and obtained the 2D-pareto frontiers in Figure 12, Figure 13 and Figure 14. Reliability was used as the scale parameter to color different system configurations, i.e., the brighter the color, the higher the reliability.
Through multi-objective optimization, the Pareto optimization frontier was acquired. Nevertheless, finding the best way to identify the optimal solution in the Pareto frontier remains a daunting challenge. Thus, AHP was introduced to solve this problem. According to Section 3.3, diverse optimum solutions can be realized via structuring diverse judgment matrices when we use AHP to perform decision analyses. Thus, to better observe the association between the goals more, the judgment matrix was presented by Table 3, and the eventual result was illustrated in Figure 15 and Figure 16.
Through a comprehensive analysis of AHP, the score of each solution on the Pareto frontier was acquired. The solution with the greatest score is the best one. As presented by Figure 15 and Figure 16, the optimal design strategy is number 96 of the Pareto frontiers, with triple redundancy configuration used in the MCES. The reliability of our designed configuration increased by 5.89% and 55.34% respectively, compared with the dual redundancy and single redundancy configurations. The results indicated that our approach combining multi-objective optimization and decision-making could realize the multi-objective optimization design of an MECS. The solutions on the Pareto frontier acquired by NSGA-II arithmetic were uniformly distributed, which suggested that we could offer designers more diversified choices. Subsequently, AHP was employed for the eventual decision analyses, which enabled designers to introduce predilections and experiences for different goals.

5. Conclusions and Future Work

The present research demonstrated the multi-objective optimization for an MECS. By designing the level length, the objectives are optimized. However, it is hard to minimize every goal simultaneously due to the conflicts between these objectives. Thus, an appropriate multi-objective optimization algorithm, which can compromise multiple objectives, is required to acquire the Pareto optimum solution set. Besides, since the Pareto frontier solution is not a solution but a solution set, an approach of multi-objective decision-making analysis is required to identify the optimum solution in the Pareto frontier. Herein, a combination of multi-objective optimization and multi-objective decision-making was utilized. The NSGA-II was employed to determine the optimum solution set, i.e., the Pareto frontier. Subsequently, the optimum design was acquired via an AHP. The outcome revealed the practicability of our optimization approach. In addition, the present research discovered that the multi-objective optimization and multi-objective decision-making approaches could be utilized in the optimization design of an MECS, and that this approach suited other components and systems as well.
Future work should encompass the improvement of real-time performance for an optimization algorithm, with more evaluation indicators considered.

Author Contributions

Conceptualization, Z.L. and R.C.; methodology, Z.L. and R.C.; software, D.L.; validation, Z.L. and D.L.; writing—original draft preparation, Z.L.; writing—review and editing, S.W.; supervision, S.W. and J.S.; project administration, J.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (No. 51620105010) and the program of China Scholarship Council (No. 202106020106).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

This study does not report any data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sun, K.; Gebre-Egziabher, D. Air data fault detection and isolation for small UAS using integrity monitoring framework. Navig.-J. Inst. Navig. 2021, 68, 577–600. [Google Scholar] [CrossRef]
  2. Kotyegov, V.; Kutovoy, S. Transfer of aerospace technologies for ensuring safety of industrial and civil objects. In Proceedings of the Second IEEE International Workshop on Intelligent Data Acquisition and Advanced Computing Systems: Technology and Applications, 2003 Proceedings, Lviv, Ukraine, 8–10 September 2003; pp. 136–138. [Google Scholar]
  3. Cossentino, M.; Lopes, S.; Renda, G.; Sabatucci, L.; Zaffora, F. Smartness and autonomy for shipboard power systems reconfiguration. In Proceedings of the International Conference on Modelling and Simulation for Autonomous Systems, Palermo, Italy, 29–31 October 2019; pp. 317–333. [Google Scholar]
  4. Joung, E.; Lee, H.; Oh, S.; Kim, G. Derivation of railway software safety criteria and management procedure. In Proceedings of the 2009 International Conference on Information and Multimedia Technology, Jeju, Korea, 16–18 December 2009; pp. 64–67. [Google Scholar]
  5. Cai, B.G.; Jin, C.M.; Ma, L.C.; Cao, Y.; Nakamura, H. Analysis on the application of on-chip redundancy in the safety-critical system. IEICE Electron. Express 2014, 11, 20140153. [Google Scholar] [CrossRef] [Green Version]
  6. Wang, S.P.; Cui, X.Y.; Jian, S.; Tomovic, M.M.; Jiao, Z.X. Modeling of reliability and performance assessment of a dissimilar redundancy actuation system with failure monitoring. Chin. J. Aeronaut. 2016, 29, 799–813. [Google Scholar] [CrossRef] [Green Version]
  7. Li, T.Y.; Chen, T.T.; Ye, K.; Jiang, X.X. Research and simulation of marine electro-hydraulic steering gear system based on EASY5. In Proceedings of the IEEE 3rd Information Technology, Networking, Electronic and Automation Control Conference (ITNEC), Chengdu, China, 15–17 March 2019; pp. 1972–1979. [Google Scholar]
  8. Jiao, Z.; Yu, B.; Wu, S.; Shang, Y.; Huang, H.; Tang, Z.; Wei, R.; Li, C. An intelligent design method for actuation system architecture optimization for more electrical aircraft. Aerosp. Sci. Technol. 2019, 93, 105079. [Google Scholar] [CrossRef]
  9. Sarlioglu, B.; Morris, C.T. More electric aircraft: Review, challenges, and opportunities for commercial transport aircraft. IEEE Trans. Transp. Electrif. 2015, 1, 54–64. [Google Scholar] [CrossRef]
  10. Benzaquen, J.; He, J.; Mirafzal, B. Toward more electric powertrains in aircraft: Technical challenges and advancements. CES Trans. Electr. Mach. Syst. 2021, 5, 177–193. [Google Scholar] [CrossRef]
  11. Arbaiy, N. Weighted linear fractional programming for possibilistic multi-objective problem. In Computational Intelligence in Information Systems; Springer: Cham, Switzerland, 2015; pp. 85–94. [Google Scholar]
  12. Marler, R.T.; Arora, J.S. The weighted sum method for multi-objective optimization: New insights. Struct. Multidiscip. Optim. 2010, 41, 853–862. [Google Scholar] [CrossRef]
  13. Du, Y.; Xie, L.; Liu, J.; Wang, Y.; Xu, Y.; Wang, S. Multi-objective optimization of reverse osmosis networks by lexicographic optimization and augmented epsilon constraint method. Desalination 2014, 333, 66–81. [Google Scholar] [CrossRef]
  14. Yang, X.; Leng, Z.; Xu, S.; Yang, C.; Yang, L.; Liu, K.; Song, Y.; Zhang, L. Multi-objective optimal scheduling for CCHP microgrids considering peak-load reduction by augmented ε-constraint method. Renew. Energy 2021, 172, 408–423. [Google Scholar] [CrossRef]
  15. Anastasopoulos, P.C.; Sarwar, M.T.; Shankar, V.N. Safety-oriented pavement performance thresholds: Accounting for unobserved heterogeneity in a multi-objective optimization and goal programming approach. Anal. Methods Accid. Res. 2016, 12, 35–47. [Google Scholar] [CrossRef]
  16. De, P.; Deb, M. Using goal programming approach to solve fuzzy multi-objective linear fractional programming problems. In Proceedings of the 2016 IEEE International Conference on Computational Intelligence and Computing Research (ICCIC), Chennai, India, 15–17 December 2016; pp. 1–5. [Google Scholar]
  17. Anukokila, P.; Radhakrishnan, B.; Anju, A. Goal programming approach for solving multi-objective fractional transportation problem with fuzzy parameters. RAIRO-Oper. Res. 2019, 53, 157–178. [Google Scholar] [CrossRef] [Green Version]
  18. Chang, Y.F.; Xie, H.; Huang, W.C.; Zeng, L. An improved multi-objective quantum genetic algorithm based on cellular automaton. In Proceedings of the 9th IEEE International Conference on Software Engineering and Service Science (ICSESS), Beijing, China, 23–25 November 2018; pp. 342–345. [Google Scholar]
  19. Altay, E.V.; Alatas, B. Differential evolution and sine cosine algorithm based novel hybrid multi-objective approaches for numerical association rule mining. Inf. Sci. 2021, 554, 198–221. [Google Scholar] [CrossRef]
  20. Tang, Q.R.; Li, Y.H.; Deng, Z.Q.; Chen, D.; Guo, R.Q.; Huang, H. Optimal shape design of an autonomous underwater vehicle based on multi-objective particle swarm optimization. Nat. Comput. 2020, 19, 733–742. [Google Scholar] [CrossRef]
  21. Su, S.J.; Han, J.; Xiong, Y.P. Optimization of unmanned ship’s parametric subdivision based on improved multi-objective PSO. Ocean Eng. 2019, 194, 106617. [Google Scholar] [CrossRef]
  22. Guria, C.; Bhattacharya, P.K.; Gupta, S.K. Multi-objective optimization of reverse osmosis desalination units using different adaptations of the non-dominated sorting genetic algorithm (NSGA). Comput. Chem. Eng. 2005, 29, 1977–1995. [Google Scholar] [CrossRef]
  23. Srinivas, N.; Deb, K. Muiltiobjective optimization using nondominated sorting in genetic algorithms. Evol. Comput. 1994, 2, 221–248. [Google Scholar] [CrossRef]
  24. Deb, K.; Pratap, A.; Agarwal, S.; Meyarivan, T. A fast and elitist multiobjective genetic algorithm: NSGA-II. IEEE Trans. Evol. Comput. 2002, 6, 182–197. [Google Scholar] [CrossRef] [Green Version]
  25. Yu, H.L.; Castelli-Dezza, F.; Cheli, F.; Tang, X.L.; Hu, X.S.; Lin, X.K. Dimensioning and power management of hybrid energy storage systems for electric vehicles with multiple optimization criteria. IEEE Trans. Power Electron. 2021, 36, 5545–5556. [Google Scholar] [CrossRef]
  26. Xia, G.; Liu, C.; Chen, X. Multi-objective optimization for AUV conceptual design based on NSGA-II. In Proceedings of the OCEANS 2016-Shanghai, Shanghai, China, 10–13 April 2016; pp. 1–6. [Google Scholar]
  27. Alam, K.; Ray, T.; Anavatti, S. Design and construction of an autonomous underwater vehicle. Neurocomputing 2014, 142, 16–29. [Google Scholar] [CrossRef]
  28. Liu, X.Y.; Yuan, Q.Q.; Zhao, M.; Cui, W.C.; Ge, T. Multiple objective multidisciplinary design optimization of heavier-than-water underwater vehicle using CFD and approximation model. J. Mar. Sci. Technol. 2017, 22, 135–148. [Google Scholar] [CrossRef]
  29. Yu, B.; Wu, S.; Jiao, Z.X.; Shang, Y.X. Multi-objective optimization design of an electrohydrostatic actuator based on a particle swarm optimization algorithm and an analytic hierarchy process. Energies 2018, 11, 2426. [Google Scholar] [CrossRef] [Green Version]
  30. Cheng, X.; He, X.Y.; Dong, M.; Xu, Y.F. Transformer fault probability calculation based on analytic hierarchy process. In Proceedings of the 2021 International Conference on Sensing, Measurement & Data Analytics in the era of Artificial Intelligence (ICSMD), Nanjing, China, 21–23 October 2021; pp. 1–6. [Google Scholar]
  31. Medineckiene, M.; Zavadskas, E.K.; Björk, F.; Turskis, Z. Multi-criteria decision-making system for sustainable building assessment/certification. Arch. Civ. Mech. Eng. 2015, 15, 11–18. [Google Scholar] [CrossRef]
  32. Altuzarra, A.; Gargallo, P.; Moreno-Jiménez, J.M.; Salvador, M. Influence, relevance and discordance of criteria in AHP-Global Bayesian prioritization. Int. J. Inf. Technol. Decis. Mak. 2013, 12, 837–861. [Google Scholar] [CrossRef]
  33. Wu, S.; Yu, B.; Jiao, Z.X.; Shang, Y.X. Preliminary design and multi-objective optimization of electro-hydrostatic actuator. Proc. Inst. Mech. Eng. Part G J. Aerosp. Eng. 2017, 231, 1258–1268. [Google Scholar] [CrossRef]
  34. Lampl, T.; Knigsberger, R.; Hornung, M. Design and evaluation of distributed electric drive architectures for high-lift control systems. In Deutsche Luft- und Raumfahrtkongress; Technische Universität München: Munich, Germany, 2017. [Google Scholar]
  35. Zhao, T. Acquisition Cost Estimating Methodology for Aircraft Conceptual Design. Ph.D. Thesis, Cranfield University, Cranfield, UK, 2004. [Google Scholar]
Figure 1. Structure of typical actuation system.
Figure 1. Structure of typical actuation system.
Aerospace 09 00085 g001
Figure 2. Configuration of actuation system with two HAs.
Figure 2. Configuration of actuation system with two HAs.
Aerospace 09 00085 g002
Figure 3. Configuration of actuation system with one HA and one EMA.
Figure 3. Configuration of actuation system with one HA and one EMA.
Aerospace 09 00085 g003
Figure 4. Configuration of actuation system with one HA and one EHA.
Figure 4. Configuration of actuation system with one HA and one EHA.
Aerospace 09 00085 g004
Figure 5. Configuration of commercial aircraft redundant actuation system.
Figure 5. Configuration of commercial aircraft redundant actuation system.
Aerospace 09 00085 g005
Figure 6. Flow chart of NSGA-II.
Figure 6. Flow chart of NSGA-II.
Aerospace 09 00085 g006
Figure 7. Binary encoding mode adopted in NSGA-I.
Figure 7. Binary encoding mode adopted in NSGA-I.
Aerospace 09 00085 g007
Figure 8. Three-level hierarchical framework of AHP.
Figure 8. Three-level hierarchical framework of AHP.
Aerospace 09 00085 g008
Figure 9. Flowchart of AHP.
Figure 9. Flowchart of AHP.
Aerospace 09 00085 g009
Figure 10. A typical 3H2E more electric control system (MECS).
Figure 10. A typical 3H2E more electric control system (MECS).
Aerospace 09 00085 g010
Figure 11. 3D-Pareto frontier of multi-objective optimization of MECS.
Figure 11. 3D-Pareto frontier of multi-objective optimization of MECS.
Aerospace 09 00085 g011
Figure 12. 2D Pareto frontier between weight, power consumption, and reliability.
Figure 12. 2D Pareto frontier between weight, power consumption, and reliability.
Aerospace 09 00085 g012
Figure 13. 2D Pareto frontier among weight, cost, and reliability.
Figure 13. 2D Pareto frontier among weight, cost, and reliability.
Aerospace 09 00085 g013
Figure 14. 2D Pareto frontier among power consumption, cost, and reliability.
Figure 14. 2D Pareto frontier among power consumption, cost, and reliability.
Aerospace 09 00085 g014
Figure 15. AHP assessment outcomes of Pareto frontier.
Figure 15. AHP assessment outcomes of Pareto frontier.
Aerospace 09 00085 g015
Figure 16. Reliability comparison among different configurations.
Figure 16. Reliability comparison among different configurations.
Aerospace 09 00085 g016
Table 1. Typical configuration of dissimilar redundant actuation system.
Table 1. Typical configuration of dissimilar redundant actuation system.
RedundancyActuator TypePower Supply
Dual redundanciesHA, HAHydraulic power
HA, EMAHydraulic and electric power
HA, EHAHydraulic and electric power
Triple redundanciesHA, EHA, and EMAHydraulic and electric power
Table 2. Values of the average stochastic coherence index.
Table 2. Values of the average stochastic coherence index.
n123456789
RI0.000.000.580.901.121.241.321.411.45
Table 3. The association between the objectives.
Table 3. The association between the objectives.
ObjectivesWeightPower DissipationCostReliability
Weight1391
Power dissipation1/311/31/3
Cost1/91/311/9
Reliability1391
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liao, Z.; Wang, S.; Shi, J.; Liu, D.; Chen, R. Reliability-Oriented Configuration Optimization of More Electrical Control Systems. Aerospace 2022, 9, 85. https://doi.org/10.3390/aerospace9020085

AMA Style

Liao Z, Wang S, Shi J, Liu D, Chen R. Reliability-Oriented Configuration Optimization of More Electrical Control Systems. Aerospace. 2022; 9(2):85. https://doi.org/10.3390/aerospace9020085

Chicago/Turabian Style

Liao, Zirui, Shaoping Wang, Jian Shi, Dong Liu, and Rentong Chen. 2022. "Reliability-Oriented Configuration Optimization of More Electrical Control Systems" Aerospace 9, no. 2: 85. https://doi.org/10.3390/aerospace9020085

APA Style

Liao, Z., Wang, S., Shi, J., Liu, D., & Chen, R. (2022). Reliability-Oriented Configuration Optimization of More Electrical Control Systems. Aerospace, 9(2), 85. https://doi.org/10.3390/aerospace9020085

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop