Next Article in Journal
Evaluation of a Serrated Edge to Mitigate the Adverse Effects of a Backward-Facing Step on an Airfoil
Next Article in Special Issue
An Approach for Using a Tensor-Based Method for Mobility-User Pattern Determining
Previous Article in Journal / Special Issue
Global Navigation Satellite System Spoofing Detection in Inertial Satellite Navigation Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Using Functionally Redundant Inertial Measurement Units to Increase Reliability and Ensure Fault Tolerance

by
Ivan M. Kuznetsov
*,
Konstantin K. Veremeenko
,
Maxim V. Zharkov
and
Andrey N. Pronkin
Flight-Navigation and Information-Measuring Systems Department, Moscow Aviation Institute (National Research University), 125993 Moscow, Russia
*
Author to whom correspondence should be addressed.
Inventions 2023, 8(6), 159; https://doi.org/10.3390/inventions8060159
Submission received: 31 October 2023 / Revised: 9 December 2023 / Accepted: 11 December 2023 / Published: 18 December 2023
(This article belongs to the Special Issue Recent Advances and New Trends in Signal Processing)

Abstract

:
This paper aims to assess the possibility of using functionally redundant inertial units to solve problems of increasing reliability and ensuring the fault tolerance of the various classes and purposes of aircraft navigation systems. We present the results of studying failure detection methods to improve the accuracy and reliability of a strapdown functionally redundant inertial unit. The resulting structural redundancy of the strapdown inertial measurement unit is designed to increase the fault tolerance and accuracy of strapdown inertial navigation systems. The methods for detecting sensor failures in functionally redundant inertial units are based on the use of the equations of functionally redundant inertial unit compliance to nominal requirements for the accuracy of measuring the input action vector. To describe the methods for detecting and eliminating failed sensor and algorithm designs based on them, we gave the mathematical models of the measurement vector of functionally redundant inertial units concerning the measured vector and the error identification condition, including the residual of the matching equations with the size due to the level of redundancy, determining the total number of matching equations. The main criterion for determining a failed sensor is non-compliance with the nominal value of the residual included in a certain number of matching equations of the information received from such meters. The developed algorithms are examined using simulation methods. The study of the selected structure of the functionally redundant inertial units shows that the proposed approaches are efficient. Also, we manage to identify the main characteristics of the algorithms for detecting sensor failures that are structurally a part of the functionally redundant inertial units.

1. Introduction

The current requirements toward maintaining the operability of navigation systems of aircraft performing various tasks to improve reliability characteristics in terms of ensuring a certain level of fault tolerance due to possible malfunctions during operation limit the use of traditional navigation systems. The performance degradation of onboard navigation information measurers structurally included in aircraft onboard navigation systems (not to mention their failure) can significantly affect the navigation system’s performance in terms of the required accuracy and availability of the navigation information. The use of information redundancy has traditionally been one of the ways to solve such problems.
The aircraft control systems of various classes and purposes, including unmanned aerial vehicles and small spacecrafts, must meet the requirements toward ensuring reliability and operability in the case of possible failures. One of the ways to increase reliability and ensure fault tolerance has traditionally been the use of information redundancy [1,2,3,4,5].
One of the essential advantages of strapdown inertial navigation systems (SINSs) is the possibility of increasing the reliability of measuring the input parameter vector at lower total costs compared to traditional inertial navigation systems based on gyrostabilized platforms when the goal is achieved via the redundancy at the system level [4,5]. In the case of the strapdown version, one can solve this problem by using an excessive number of measuring elements as part of an inertial measurement unit (IMU) [6,7,8]. In this case, the IMU is called functionally redundant because it can perform its function using various combinations, for example, accelerometers measuring the components of the specific force vector and angular rate sensors (ARSs) that provide information about the motion parameters of a moving object relative to the center of mass. Also, both the hardware (the failure of the most sensitive element) and information (as a measure of the accuracy of the output information) reliability are monitored [5].
Redundancy refers to the difference between the number of primary information sensors used in the IMU and its number, the minimum necessary for setting a measuring basis that provides the measurements of the specific force vector, and the absolute angular velocity vector. In general, the measuring system can be non-orthogonal [4,9], while the optimality of the functionally redundant inertial unit (FRIU) technical solution is still achieved, for example, in terms of the weight and size and economic indicators.
Many researchers have studied methods of fault detection and isolation of inertial measurement systems with redundant sensors in their structure. The most common of them include such methods as generalized likelihood testing (GLT), which is not able to detect and isolate the faults of several autonomous inertial measurement units (SIMUs) along a certain axial direction. This is because, when solving the decoupling matrix, row vectors of the orthogonal projection matrix are used, which do not satisfy the conditions of linear correlation. The GLT method is used to detect a system fault, while the linear prediction method is used to estimate the inertial sensor values. A faulty inertial sensor signal can be isolated by comparing the calculated value with changes in the inertial sensor signal, and information about the fault is recorded for system recovery. The work [10] analyzes the disadvantages of the GLT method and proposes a methodology for constructing the decoupling matrix by selecting maximal linearly independent systems from an array of orthogonal projections and then orthogonalizing and combining the maximal uncorrelated set.
The article [11] proposes a method for detecting and isolating multiple faults based on a GLT and a linear prediction approach. The work [12] shows the approach in terms of the fault identification accuracy of tiny faults with an acceptable computational complexity.
There is also a known group of methods, such as singular value decomposition (SVD) [13], which use the null space component of faults with respect to a configuration matrix to obtain the geometric fault detection and isolation (FDI) technique. As well as SVD and GLT, the most frequent FDI techniques imply the use of an optimal parity vector [14,15,16] approach usually intended to recognize false alarms and wrong isolations.
One of the most recent studies [4] reviews the most common schemes of the FRIU configuration and its ensuring fault tolerance methods; a review of redundant inertial navigation technology can be found in [5].
Summarizing the analytical review, it can be concluded that in the reviewed works devoted to the consideration of ensuring FRIU-based fault tolerance methods, the effectiveness of a single method is usually considered in terms of the basic criteria for the effectiveness of its work, such as the time of failure detection [10,11,12,13,14,15,16]; its fault identification accuracy (i.e., depending on the FRIU working state [17]); the FRIU configuration and, related with this, the lever arm effects [18]; and the frequency of false alarms and wrong isolations, together with a study of general inertial system reliability [19] and its accuracy [20,21,22,23].
The research object is a fault-tolerant spacecraft SINS, configured with another one to be functionally and structurally redundant, so the research is quite specific, requiring special investigation and incurring problems that were not considered before. The methods of detecting failures to improve the accuracy and reliability of the strapdown FRIU are considered. In this case, the arising structural redundancy of the strapdown IMU (SIMU) is intended to increase the SINS fault tolerance and accuracy.
This paper is organized as follows. In Section 2, the classification of the sensor failure types as well as the concept of a fault tolerance system design providing the two-level procedures with its mathematical models are given. The methodology of the failure detection algorithm operability confirmation, including several scenarios and simulation results, is presented in Section 3. A discussion of the simulation results, including a comparison to the results of other recent studies, is given in Section 4. Section 5 is devoted to the research conclusions.

2. Materials and Methods

This paper aims to analyze the possibility of using FRIUs to increase the reliability and ensure the fault tolerance of aircraft navigation systems of various classes and purposes. To achieve the research goal, we set the following tasks:
  • Classify the types of sensor failures included in the FRIU;
  • Develop a mathematical model for implementing algorithms for detecting and eliminating failed sensors in terms of the established classification;
  • Design a research methodology for the proposed algorithms for detecting and eliminating failed sensors;
  • Perform a simulation in accordance with the developed methodology of research on algorithms for detecting and eliminating failed sensors;
  • Analyze the research results and estimate the operability of the proposed technical solutions and the main performance parameters of the developed algorithms.

2.1. Classification of Sensor Failure Types

The inertial measurement fault tolerance system is designed to detect a failure in the measurements of an IMU and identify a failed sensor. The failure of an inertial sensor means exceeding the estimate of the measurement error limits in the acceptable range of values.
Depending on the effect on the navigation parameter calculation error using the SINS algorithm, we suggest dividing the inertial sensors into three types:
-
Instantaneous failure, which is when the measurement error estimate with the inertial sensor exceeds the limits of the tolerance range of values, allowing one to identify the failure in one measurement cycle unambiguously;
-
Medium-level failure, which can be detected only after the accumulation and analysis of the statistical measurement indicators;
-
Low-level failure, which is when using measurements of a failed inertial sensor leads to an increase in the errors of the SINS at a speed higher than that guaranteed by the passport characteristics.
The concept of a fault tolerance system design for inertial measurements based on two inertial information units identical in type and characteristics relies on a redundant number of measuring elements as part of a single (consisting of two inertial information units) IMU. In this case, the IMU is called functionally redundant because it can perform its function using various combinations of accelerometers that measure the components of the specific force vector and ARSs that provide information about the movement parameters of a moving object relative to the mass center.
Redundancy refers to the difference between the number of primary information sensors used in the IMU and their number, which is minimally required for creating a measuring basis that provides the specific force vector and the absolute angular velocity vector measurements.
The geometry of the FRIU measuring axes’ location relative to each other is represented by two conditional groups in orthogonal and asymmetrically arranged directions. In the latter group, the sensors’ measuring axes can be directed along the cone generatrix with different half-angles [2,3,7] or perpendicular to the regular polyhedra faces, depending on the required redundancy level [9]. There are algorithms based on which failed sensors are detected and eliminated [24,25].
For sensors measuring the projections of the input action vector on their measuring axes (accelerometers and absolute angular velocity sensors), the non-redundant measuring basis (at a redundancy level equal to (in Equation (1)) m = k 3 , that is, at k = 3 , m = 0 ) is designed using 3 sensors whose measuring axes are non-collinear and non-planar. The total number of different measurement bases in a unit of k sensors with a redundancy level of m = k 3 is calculated using Equation (1):
n = C k 3 = C m + 3 3 ,
and the total number of different measuring bases in a unit of 6 sensors reaches 20.
The concept of a fault tolerance system design provides for a three-level failure detection procedure:
-
Failure detection in one measurement cycle (Level I procedure);
-
Failure detection based on an analysis of the measurement residual statistical characteristics (Level II procedure);
-
Failure detection based on complex information processing from the SINS built according to the FRIU measurements (Level III procedure).
We discuss the procedures of the first two levels and describe their main relationships in this section.

2.2. Failure Detection in One Measurement Cycle

As mentioned above, the total number of different measuring bases in a unit of 6 sensors reaches 20. This step allowed us to identify the failure of two sensors in one measurement cycle and detect the failure of the third. The basis for an algorithm for the Level I procedure design is the information about the attitude of the measuring axes of the sensors in the axes of the coordinate system (CS) “O” associated with the installation base on which both inertial information units are located. Such a characteristic of the IMU is called an alignment one, and for a unit of 6 sensors, it is presented as a direct cosine matrix (DCM) of the size (6 × 3).
Let the unit vector characterize the attitude of the measuring axis of the i-th sensor in the CS “O”. Then, the transposed alignment matrix of a 6-sensor unit can be represented as follows (Equation (2)):
E O T = [ e ¯ o 1 , , e ¯ o k ] ,
where E O T —the transposed alignment matrix; e ¯ o i —the vectors of the direct cosines of the size (1 × 3) (i = 1 … k), which make up the DCM of the size (i × 3), describing the attitude of the measuring axes of the i sensors in the CS “O”, and the mathematical model of the FRIU measurement vector ( m ¯ ) concerning the measured vector V ¯ O takes the following form (Equation (3)):
m ^ = E O V ¯ O + Δ m ¯ ,
where m ^ —the measured value of the vector components of the FRIU; E O —the alignment matrix of the 6-sensor unit; V ¯ O —the measured vector, Δ m ¯ —a vector of the size (6 × 3) of the sensor instrumental errors.
Using the normal distribution law of the unit sensor instrumental errors (typical for measuring elements made using unified technology), the estimation of the maximum likelihood of the measured vector is achieved using the weighted least squares algorithm (Equations (4)–(6)):
V ^ = H m ^ ,
V ^ = ( E O T Q E O ) 1 E O T Q ,
Q = K m 1 ,
where V ^ —the estimated vector; H—the transformation matrix of the size (3 × 6) of the measurement vector of the single-component unit sensors ( m ¯ ) into the measured vector estimate; m ^ —the FRIU sensor measurement vector of the size (6 × 1); E O —the direct cosine matrix (alignment) of the size (6 × 3), describing the attitude of the measuring axes of the sensors in the CS “O”; and K m —the measurement error covariance matrix of the size (6 × 6).
The algorithm (Equations (4)–(6)) provides the smallest variance of the vector V ¯ error estimate on ( σ v 2 ) in the asymptotic sense and the finite domain.
The measured vector in the expression (Equation (2)) can be the specific force vector ( V ¯ = n ¯ ) or the absolute angular velocity vector ( V ¯ = ω ¯ ).
With equally accurate measurements, when the weight matrix is the identity one, the optimal algorithm is transformed into the least squares method (LSM) algorithm and reduced to multiplying the vector m ¯ by a constant pseudo-inverted matrix, the use of which is justified in the general case, with the lack of reliable information about the statistical characteristics of the measuring element errors:
H = ( E O T E O ) 1 E O T ,
An algorithm to be designed for detecting and identifying failures in one measurement cycle should use the IMU matching equations of the input action vector’s measurement accuracy nominal requirements.
Let ρ ¯ be the s-dimensional residual matching equation vector (the dimension of the vector is due to the redundancy level that determines the matching equation number), and let ε ¯ t o l be these s-dimensional residual tolerance values’ vector. Then, the conditions for the IMU error identification exceeding the nominal values can be expressed by the following vector equation:
ρ ¯ = [ C ] m ^ ε ¯ t o l ,
where ρ ¯ —the s-dimensional residual matching equations, s—the number of matching equations, [ C ] —the matching matrix of the size (s × k) corresponding to the condition (Equation (9)), m ^ —the k-dimensional unit sensor real measurement vector, and ε ¯ t o l —the s-dimensional vector of these residual tolerance values.
[ C ] E O = 0 ,
When the condition (Equation (9)) is met, the vector ε ¯ t o l can be found from the expression:
[ C ] Δ m ¯ ε ¯ t o l ,
The validity of this expression can be established by substituting (Equation (4)) into the main theoretical expression (Equation (8)). Indeed, when the condition (Equation (9)) is met, the useful signal is excluded from the expression (Equation (8)), and only the residual due to the measurement error remains.
To find the [ C ] matrix elements, a system of algebraic equations corresponding to the measurements of four sensors should be used. Using the condition (Equation (9)) for a unit of 6 sensors, 15 equations have to be made, each responsible only for the sensors that the signals are part of.
With known values of the [ C ] matrix elements for the 6-sensor IMU, 15 matching equations can be formed, each of which includes the data of only 4 measuring elements (mi):
C 1.1 m 1 + C 1.2 m 2 + C 1.3 m 3 + C 1.4 m 4 = ρ 1 ε ¯ 1 t o l ; C 2.1 m 1 + C 2.2 m 2 + C 2.3 m 3 + C 2.4 m 4 = ρ 2 ε ¯ 2 t o l ; C 3.1 m 1 + C 3.2 m 2 + C 3.3 m 3 + C 3.6 m 6 = ρ 3 ε ¯ 3 t o l ; C 15.3 m 3 + C 15.4 m 4 + C 15.5 m 5 + C 15.6 m 6 = ρ 15 ε ¯ 15 t o l .
The information from each measurer is included in 10 specific equations, which, in case of failure, do not correspond to the nominal value of the residual. The failed sensors are determined by the number of these equations.
When establishing the tolerance for the residual of the matching equations, we assumed that the instrumental errors of the SINS’s measuring elements made using unified technology are independent and normally distributed. Then, according to Equations (8)–(10), the residual vector component SD can be calculated as follows:
σ ε i t o l = σ Δ m max ( j = 1 n C i j 2 ) 1 / 2 ,
where σ ε i t o l —the residual vector component SD, and σ Δ m max —the unit first measurement norm SD of i = 1, …, 15; j = 1, …, 6.
In (Equation (12)), the Δ m max first norm of the measurement vector of the unit is calculated as the maximum modulo value of the error component of the measurement vector at the time of calculating the residual vector.
However, the residual vector determination, which ensures the specified navigation accuracy and reliable identification of failed sensors, is a responsible engineering task, which should be solved by using reasoned simulation and hardware-in-the-loop simulation.

2.3. Failure Detection Based on the Analysis of Measurement Residual Statistical Characteristics

To detect faulty measurements caused, for example, by the failure of inertial sensors (in the form of an abrupt change in readings or an increasing error (medium and low-level failures)), an algorithm to test the null hypothesis about the non-bias of the matching equation residuals should be used, an alternative to which is the hypothesis about the presence of the average value in the residuals above the thresholds, set according to the manufacturers of inertial sensors and confirmed by the results of experiments (simulation and hardware-in-the-loop simulation). The hypothesis that the mean values of the differences are equal to zero is tested based on the normalized difference criterion (Student’s criterion [26]) between the actually obtained statistical characteristics and the predicted ones.
To calculate the statistical characteristics of the measurements based on the last l cycles, one should calculate the estimates of the moving averages of the residuals of the matching equations M ρ i ( t k ) and the variances of the D ρ i ( t k ) residuals for each correspondence equation i at the time t k :
M ρ i ( t k ) = ( j = 1 n ρ i ( t j ) ) / l ,
D ρ i ( t k ) = ( j = 1 n ρ i ( t j ) M ρ i ( t k ) ) 2 ) / l ,
where M ρ i ( t k ) —the matching equations on the l cycles residual estimates, ρ i ( t j ) —the i-th residual of the matching equations at the time t j , and l—the number of cycles used to calculate the moving average.
The hypothesis of the non-bias of the residuals of the matching equations at a time t k should be tested as follows:
q ρ i ( t k ) = ( M ρ i ( t k ) M i ρ ) / Y ρ i ( t k ) ,
where q ρ i ( t k ) —the calculated value of the criterion for testing the hypothesis, Y ρ i ( t k ) = ( D ρ i ( t k ) + D i ρ / d ) 1 / 2 , where M i ρ and D i ρ are determined according to the data from the inertial information unit manufacturer.
The condition for making a decision on the presence of errors above the threshold set by the manufacturer in the measurements of the sensors included in the i-th characteristic equation, with r the number of the freedom level, has the following form:
If q ρ i ( t k ) > q c r i t , then the hypothesis is rejected: the presence of faulty measurements in the characteristic equation are diagnosed. The critical value of the criterion q c r i t is selected based on the condition J ( q c r i t ) = 1 α r . α r —the significance level of the criterion (probability of a type-1 error), and J—the Laplace function.
Similar to Section 2.2, the failed sensor should be determined according to the matching equation numbers.
It is advisable to organize fast (n 10), medium (n 30), and slow cycles (n 60 s accordingly) for testing the hypothesis of the non-bias of the residuals, which will ensure the rapid detection of failures not detected by the Level I procedure and, at the same time, detect medium-level failures.

3. Results

To analyze the algorithms for detecting and eliminating failures of FRIU sensors, we provided an analytical assessment of the fault tolerance system performance in terms of identifying the possibilities for detecting and identifying failures with a different number of simultaneously failed inertial sensors.
The Level I and Level II procedures of failure detection and identification are based on the use of the so-called matching equations (Equation (10)). In the case of using two inertial measurement units for the FRIU design, the number of such equations is equal to 15.
If there is a failure of one sensor for 10 matching equations, the condition will not be met (see Equation (8)), and the hypothesis will not be confirmed (see Equation (14)); however, for the 5 remaining equations, this condition will be met and, respectively, the hypothesis will be confirmed. If there is a failure of two sensors for the 14 matching equations, the considered condition will not be met, and the hypothesis will also not be confirmed. At the same time, for the one remaining equation, the condition will be met, and the hypothesis will be confirmed. Finally, if there are failures of three or more sensors for all 15 matching equations, the verification condition will not be met, and the hypothesis will not be confirmed.
In this regard, it can be concluded that the Level I and Level II procedures make it possible to identify the failure of two sensors and detect the failure of the third one.
To verify the operability of the considered algorithms for ensuring the fault tolerance of the SINS measurements with an excessive number of measurers and identify their main characteristics, we performed a full-fledged simulation, which involved considering the functionality of the proposed system for ensuring the fault tolerance of the FRIU measurements in various scenarios. The key features of such scenarios can be divided into several groups:
  • The study of the sensitivity of the considered system for ensuring the fault tolerance of SINS measurements with an excessive number of meters at different levels of sensor accuracy (systematic and noise components) of the FRIU (Level I and Level II procedures);
  • An assessment of the requirements for the accuracy of the FRIU alignment based on two inertial information units;
  • An assessment of the influence of the relative angular oscillations of two inertial information units.

3.1. Sensitivity Study with the Use of Various Accuracy Classes’ Sensors

In the first group of scenarios, we tried to identify the dependence of the main parameters of the algorithm for detecting sensor failure (residuals of the matching equations and their standard deviation (SD) from the values of the main errors of the sensors under study), with various levels of systematic and noise components in the range from 10−3 to 10 deg/hour for the ARSs and from 10−5g to 10−1g m/s2 for the accelerometers for the same trajectory conditions in which the FRIU operates. Table 1 provides examples of the calculations of the tolerance values of the residuals of the matching equations and their SD systems when operating algorithms of the Level I procedure (Figure 1 and Figure 2 clearly display it). Table 2 and Table 3 display an example of calculating the tolerance values of the residuals of the matching equations according to the Student’s criterion for fast, medium, and slow cycles (Level II procedure) (Figure 3 and Figure 4).
The results obtained within the framework of this study allow for identifying the levels of sensitivity of the proposed approach. For the Level I procedure, for example, the ARS threshold values and their SD of the matching equations almost fit the linear function (see Figure 1) of its accuracy class, and it is the most sensitive for the consumer grade (10 deg/hour) and equals approximately 20% for its mean value. At the same time, the threshold levels for the tactical (10−1 deg/hour) and navigation (10−3 deg/hour ARS drift rate, respectively (i.e., IMU 1000 [27])), on average, are two times bigger than its accuracy level. The accelerometer threshold values and their SD of the matching equations also almost fit the linear function (see Figure 2) of its accuracy class and are also the most sensitive for the consumer grade (10−1g m/s2) and equal approximately 37% for its mean value. At the same time, the threshold level for tactical (10−3g m/s2) is six times and navigation (10−5g m/s2 accelerometer zero drifts, respectively) is two orders of magnitude greater than its accuracy level.
For the Level II procedure, the ARS threshold values and their SD of the matching equations have an even more linear function fit for the accuracy class, compared to the Level I procedure (see Figure 3), and the sensitivity increases within the sensor accuracy class logarithmically from the “fast” to “slow” cycle simulation test conditions. For the consumer grade (10 deg/hour), the threshold levels have the range of 116% for “fast”, 21% for “medium”, and 7% for “slow” cycles, respectively; at the same time, for the tactical grade (10−1 deg/hour), it has the range of 1115% for “fast”, 210% for “medium”, and 73% for “slow” cycles, respectively. Last, the threshold levels for the navigation grade (10−3 deg/hour ARS drift rate, respectively) are in the range of 126% for “fast”, 24% for “medium”, and 21% for “slow” cycles, respectively. The accelerometer threshold values obtained within the Level II procedure simulation and their SD of the matching equations also almost fit the linear function (see Figure 4) of its accuracy class, and as well as for the ARS, its sensitivity increases within the sensor accuracy class logarithmically from the “fast” to “slow” cycle simulation test conditions. For the consumer grade (10−1g m/s2), the threshold levels have the range of 116% for “fast”, 21% for “medium”, and 7% for “slow” cycles, respectively; at the same time, for the tactical grade (10−3g m/s2), it has the range of 101% for “fast”, 20% for “medium”, and 6.8% for “slow” cycles, respectively. Last, the threshold levels for the navigation grade (10−5g m/s2 accelerometer zero drifts, respectively) are in the range of 1141% for “fast”, 199% for “medium”, and 62% for “slow” cycles, respectively.
The obtained results allow for concluding that the combination of various Level II procedure cycles can lead to a preferred FDI accuracy level and can be customized and optimized due to the overall required system reliability level.

3.2. Assessment of Adjustment Accuracy Requirements

The realization of the potential advantages of functionally redundant meters is associated with the complexity of the measuring element unit (MEU) design. Generally, this circumstance leads to technical difficulties in ensuring the alignment characteristics of the unit and requires a reasonable assignment of tolerances to ensure the nominal attitude of its measuring axes.
In this regard, in the second group of the simulation modeling, the requirements for the accuracy of the MEU adjustment were evaluated, and the unit output information accuracy was investigated.
For a class of functionally redundant accelerometer units (AUs) of a conical structure with a nominal alignment characteristic of the form (Equation (16)):
H = [ sin χ cos β i , cos χ , sin χ sin β i ] T ,
where H—the transformation matrix, and where i = 1, …, k, βi = 2 π i 1 / k , the expression for the functional transformation matrix will take the form (Equation (17)):
H = 2 cos β 1 k sin χ , , 2 cos β 1 k sin χ , , 2 cos β k k sin χ 1 k cos χ , , 1 k cos χ , , 1 k cos χ 2 sin β 1 k sin χ , , 2 sin β i k sin χ , , 2 sin β k k sin χ ,
where χ —the half-angle of the cone, along the generatrix of which the measuring axes of the sensors are uniformly located (with the βi angles).
For the MEU design chosen within the modeling, the first norm of the H transformation matrix is determined by the second row:
H I = 1 cos χ ,
The value of the angle χ in the selected configuration is 54.74° [4,9], which determines the values of the expression’s corresponding norms. This step allowed us to estimate the maximum error of the MEU output depending on the adjustment errors, configuration of the unit, and flight conditions:
δ V ^ I = 2 H I r ¯ i I V ^ I Θ i j max ,
H I = 1.7013, r ¯ i I = 0.89, ( V ^ I = 1 g–cruise flight mode).
As a result, the following calculated ratios can be obtained: Θ i j max = δ n i / 2.75 .
As part of the simulation, the requirements for the accuracy of the alignment were studied. The results from calculating the tolerance value of the adjustment error (3σ) are shown in Table 4. During the calculation, we assumed that the error of the AU output information ( δ n i ) should not exceed the accelerometer zero-offset error. Table 4 also shows the results for the orthogonal characteristics of the MEU structure and the conical structure with an angle χ equal to 78°.
The results obtained within the framework of the FRIU accuracy alignment requirements study based on two inertial information units shows that it depends on its structure and increases significantly better (up to 27% for conical ( χ = 54.74 ° ) and 80% for conical ( χ = 78 ° )) as well as for the precise adjustment (10 ang. s. error) and for the rough one (80 ang. s. error) and allows for stating the fact of more strict requirements for the more non-orthogonal unit structure.

3.3. Evaluation of the Relative Angular Oscillation Effect

To ensure the design of the FRIU based on already existing orthogonal IMUs, providing the necessary technological design of their relative angular position, the developers of the FRIU should estimate the effect of angular fluctuations of the IMUs that will lead to a change in the values of the matrix elements E O (Equation (2)). The analytical expression for calculating the matching equations’ residual error can be obtained by formally varying the expression (Equation (20))
ρ ¯ = [ C ] E O V O + [ C ] Δ m ,
which can be obtained by substituting (Equation (3)) with (Equation (8)) in the resulting expression
ρ = [ C ] E O V ¯ O
E O = [ δ e ¯ O 1 , δ e ¯ O i ] T ,
where E O —the matrix of the size (6 × 3) of the unit measuring axes’ adjustment errors.
The row of the matrix E O describes the measuring axis adjustment error of the i-th sensor; it is expressed in terms of the vector product components δ e ¯ O i = Θ ¯ i × e ¯ O i , one of the multipliers of which is the small rotation vector ( Θ ¯ i ), which characterizes the difference between the real measuring axis attitude of the i-th sensor from the nominal one caused by angular fluctuations of two IMUs in their own housing.
The error in calculating the residual of the matching equations caused by the relative angular fluctuations of the two IMUs can be estimated using a simulation, presenting a model of the change in the components of the small rotation vector ( Θ ¯ i ), for example, in the form of a harmonic oscillation with specified parameters (the circular oscillation frequency of 100 rad/s; the oscillation amplitude of 1°; and the number of cycles–105, respectively).
Table 5 shows the results of the simulation with the ARSs as the measuring elements used.
Table 6 shows the results of the simulation with accelerometers as the measuring elements used.
The results obtained within the framework of the relative angular oscillation effect study shows the level of its contribution to the overall matching equation residual calculating error. Depending on the specific construction and increased dynamics, this issue has the level up to percent fractions in its mean value but adds considerable noise, which increases the probability of type-1 errors.
The simulation showed that the proposed approach to detect failures of inertial sensors structurally included in the FRIU performs well. The approach considered is a set of functional algorithms that implement the concept of a system design that is ensuring the fault tolerance of the functioning of both the sensors and systems that operate using the measuring information received from such sensors, primarily SINSs. In addition, this approach makes it possible to detect failures in one measurement cycle and to identify failures based on the analysis of the statistical characteristics of the measurement residuals.
The study of the sensitivity of failure detection algorithms when using sensors of various accuracy classes in the first group of simulation scenarios made it possible to determine the tolerance values of residuals and their SD, as well as to establish thresholds for making a decision to exclude failed sensors from the SINS’s functional algorithms. The given simulation results allowed us to select sensitive elements for the FRIU design for the configuration scheme under consideration with the selected redundancy characteristics.
The assessment of the requirements for the FRIU accuracy alignment based on two inertial information units in the second group of simulation scenarios shows that the requirements for the accuracy of the alignment depend on the structure of the AU and increase significantly as the non-orthogonality of the unit structure increases. In more dynamic flight modes, the alignment accuracy requirements may exceed the design and technological capabilities of their implementation. In this case, analytical alignment is the only way to solve the problem.
The third group of simulation scenarios made it possible to determine the assessment of the influence of the relative angular oscillations of the two inertial information units and to obtain numerical values of error estimates for calculating the residual of the matching equations for the selected model of the external disturbance of the considered FRIU.

4. Discussion

The analysis of the simulation results allowed us to conclude that the proposed approach to detecting failures of inertial sensors structurally included in the FRIU operates well. The numerical values obtained in Table 1, Table 2 and Table 3 reflect the dependences of the main parameters of the algorithm for detecting sensor failure, the residual of the matching equations, and its SD on the values of the main errors of the sensors under study in various levels of the systematic and noise components of the ARSs and accelerometers. In this case, the tolerance values of the residuals of the matching equations can be obtained directly during the operation of the algorithm. The approach is somewhat universal because there is no need to work with a useful signal; it is enough to have an a priori assessment of the characteristics of the sensors and work on the residuals due to the error of its measurements.
A feature of the proposed approach that distinguishes it from those presented in other works [10,11,12,24,25] is the ability of the algorithm to detect so-called medium-level failures using the statistical hypothesis-testing technique. Among the other things, the algorithm provides the smallest variance of the vector error of estimation in the asymptotic sense and the finite domain. The way it is designed is that it works specifically with the residuals—information about all the types of failures is contained in it and a useful signal can be excluded from consideration, thereby avoiding errors affecting the occurrence of false alarms [21,22]. One of the advantages of the proposed approach is the high value of fault identification accuracy with a high probability due to the use of a large number of correspondence equations, which greatly simplifies the procedure for identifying failed sensors, eliminates false alarms with an appropriately pre-performed procedure for factory calibration of sensors, and, in addition, unlike the other works, is a customizable parameter [17,19,21,22].
Moreover, the algorithm takes into account the influence of tolerances on ensuring the nominal attitude of the FRIU measuring axes, which allows for solving the problems associated with alignment that inevitably arise, for example, when installing ready-made measuring units in a structure that provides a redundancy of measurements [17,19,22,23]. The assessment of the requirements for the accuracy of the alignment of the MEU (Table 4), with the parameter chosen as the main criterion (the accuracy of the output information of the unit), allowed us to establish the dependence of the structure of the MEU and the dynamics of the movement of the object on these requirements. In this case, analytical accounting of the adjustment parameters may be the only way to improve the accuracy of the meters.
Finally, due to the conditions of the numerical experiment, in contrast to published works in the subject area [3,7,10,11,12,17,19,22,23], the algorithm provides for the possibility of considering the influence of mutual angular oscillations of the FRIU based on already existing orthogonal IMUs, for which the necessary technological design of their mutual angular position is provided. The numerical values obtained during the experiment of the estimates of the errors in calculating the residual of the matching equations for the selected model of the external disturbance of the considered FRIU (Table 5 and Table 6) make it possible to increase the final accuracy of the algorithm for detecting and eliminating sensor failures in the presence of such disturbances.

5. Conclusions

As a result of this study, we suggested ways to detect failures to improve the accuracy and reliability of the strapdown FRIU. In this case, the arising structural redundancy of the SIMU is intended to increase the SINS’s fault tolerance and accuracy levels.
When substantiating the concept of fault tolerance, we classified the types of failures. Thus, we showed the concept of a fault tolerance system design that provides a three-level failure detection procedure to detect failures of all types.
The possibilities of detecting and identifying failures using the proposed approach were studied using simulation-modeling methods. We considered a wide range of modeling scenarios, providing the most comprehensive evaluation of the detecting sensor failure algorithms’ main characteristics that are structurally an FRIU part.
The research results allowed us to conclude that the proposed approach to the fault tolerance system design is working. Also, an analytical assessment of the fault tolerance system characteristics in terms of an evaluation of the capabilities of detecting and identifying failures with a different number of simultaneously failed inertial sensors revealed that Level I and Level II procedures allowed us to identify the failure of two sensors and to detect the failure of the third one.
The research novelty lies in the ability of the algorithm to detect medium-level failures with the help of the statistical hypothesis-testing technique, taking into account the influence of tolerances on ensuring the nominal attitude of the measuring axes of functionally redundant inertial units and the relative angular oscillation of the inertial measurement unit structurally included in the strapdown inertial navigation system.
Finally, the use of first- and second-level procedure sets makes it possible to provide the so-called FRIU predictive diagnostics and to warn the consumer in advance of the need to exclude such measurements from the navigation solution.
A further direction of the subject area can involve studying a Level III procedure based on the complex processing of information from parallel-executed SINS algorithms constructed from measurements of a functionally redundant measuring unit containing six identical ARSs and accelerometers included in two inertial information units. This procedure is assumed to allow for identifying and eliminating failures of low-level sensors.

Author Contributions

Conceptualization, K.K.V.; methodology, K.K.V.; software, M.V.Z. and A.N.P.; validation, I.M.K. and A.N.P.; formal analysis, I.M.K.; investigation, I.M.K.; resources, I.M.K.; data curation, A.N.P.; writing—original draft preparation, I.M.K.; writing—review and editing, I.M.K.; supervision, M.V.Z.; project administration, K.K.V.; funding acquisition, K.K.V. All authors have read and agreed to the published version of the manuscript.

Funding

The research was carried out within the state assignment of the Ministry of Science and Higher Education of the Russian Federation (theme no. FSFF-2023-0005).

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
ARSAngular Rate Sensors
AU Accelerometer Unit
CSCoordinate System
DCMDirect Cosine Matrix
FDIFault Detection and Isolation
FRIUFunctionally Redundant Inertial Unit
GLTGeneralized Likelihood Testing
IMUInertial Measurement Unit
LSMLeast Squares Method
MEMathematical Expectation
MEUMeasuring Element Unit
SDStandard Deviation
SIMUStrapdown Inertial Measurement Unit
SINSStrapdown Inertial Navigation System
SVDSingular Value Decomposition
NotationDescription
k, l, m, n, sNumber of sensors, number of cycles used to calculate the moving average, redundancy level, total number of different measuring bases, number of matching equations
E O Alignment DCM
e ¯ o i Vectors of the direct cosines of the size (1 × 3), (i = 1 … k)
m ^ Measured value of the FRIU vector components of the size (6 × 1) (rad/m/s2)
V ¯ O Measured vector of the size (6 × 1) (rad/m/s2)
Δ m ¯ Vector of the instrumental errors of the sensors (rad/m/s2)
V ^ Estimated measurement vector of the size (3 × 1) (rad/m/s2)
HTransformation matrix of the size (3 × 6)
K m The measurement error covariance matrix of the size (6 × 6)
σ ν 2 Measurement vector variance (rad2/ m2/s4)
n ¯ The specific force vector of the size (3 × 1) (m/s2)
ω ¯ The absolute angular velocity vector of the size (3 × 1) (rad)
ρ ¯ The vector of the residuals of the matching equations
[ C ] The correspondence matrix of the size (s × k)
ε ¯ t o l s-dimensional vector of residual tolerance values
σ ε i t o l SD of the residual vector components
σ Δ m max SD of the measurement unit first measurement norm
M ρ i ( t k ) Moving averages of the matching equation residuals at the time t k
ρ i ( t j ) i-th residual of the matching equations at the time t j , within the cycle of moving average calculation (from j = 1…l)
q ρ i ( t k ) Calculated value of the criterion for testing the hypothesis
D i ρ Matching equation residual variance
q c r i t The critical value of the criterion for testing the hypothesis
JLaplace function
α r Level of significance of the criterion (probability of type-1 error)
βiFRIU configuration coordinate angles (rad)
χ Cone half-angle (sensor measurement axis oriented along its generatrix) (rad)
ρ The error vector of the matching equation residuals
E O Measuring axis adjustment error matrix
δ e ¯ O i Adjustment error direct cosine vector of the size (1 × 3), (i = 1 … k)
Θ ¯ i Small rotation vector (rad)

References

  1. Allerton, D.J.; Jia, H. An Error Compensation Method for Skewed Redundant Inertial Configuration. In Proceedings of the 58th Annual Meeting of The Institute of Navigation and CIGTF 21st Guidance Test Symposium 2002, Albuquerque, NM, USA, 24–26 June 2002; pp. 142–147. [Google Scholar]
  2. Osman, A.; Wright, B.; Nassar, S. Multi-sensor inertial navigation systems employing skewed redundant inertial sensors. In Proceedings of the ION GNSS 19th International Technical Meeting of the Satellite Division, Fort Worth, TX, USA, 26–29 September 2006; Volume 35, pp. 69–87. [Google Scholar]
  3. Sturza, M.A. Navigation system integrity monitoring using redundant measurements. J. Inst. Navig. 1988, 35, 483–501. [Google Scholar] [CrossRef]
  4. Li, X.; Xu, F.; Xu, G. Redundant Inertial Measurement Unit Reconfiguration and Trajectory Replanning of Launch Vehicle; Navigation: Science and Technology; Springer Nature: Singapore, 2023. [Google Scholar] [CrossRef]
  5. Zhai, X.; Ren, Y.; Wang, L.; Zhu, T.; He, Y.; Lv, B. A Review of Redundant Inertial Navigation Technology. In Proceedings of the 2021 International Conference on Computer, Control and Robotics (ICCCR), Shanghai, China, 8–10 January 2021; pp. 272–278. [Google Scholar] [CrossRef]
  6. Colomina, I.; Giménez, M.; Rosales, J.J. Redundant IMUs for precise trajectory determination. Int. Arch. Photogramm. Remote Sens. 2004, 35-B1, 159–165. [Google Scholar]
  7. Pejsa, A.J. Optimum skewed redundant inertial navigators. AIAA J. 1974, 12, 899–902. [Google Scholar] [CrossRef]
  8. Zheng, X. Reliability study of redundant configuration of IMU fiber optic gyroscope inertial measurement unit. In Proceedings of the 2009 8th International Conference on Reliability, Maintainability and Safety (ICRMS), Chengdu, China, 20–24 July 2009. [Google Scholar]
  9. Sukkarieh, S.; Gibbens, P.; Brocholsky, B. A low-cost redundant inertial measurement unit for unmanned air vehicles. Int. J. Robot. Res. 2000, 19, 1089–1103. [Google Scholar] [CrossRef]
  10. Ye, S.; Zhang, M. Analysis of Generalized Likelihood Test Method Deficiencies and Improvements for Fault Detection and Isolation of Redundant Strapdown Inertial Measurement Unit. Mechatron. Autom. Eng. 2017, 327–335. [Google Scholar] [CrossRef]
  11. Cheng, J.; Sun, X.; Chen, D.; Cheng, C.; Mou, H.; Liu, P. Multi-fault detection and isolation for redundant strapdown inertial navigation system. In Proceedings of the 2018 IEEE/ION Position, Location and Navigation Symposium (PLANS), Monterey, CA, USA, 23–26 April 2018; pp. 1562–1569. [Google Scholar] [CrossRef]
  12. Hu, X.; Zhang, X.; Peng, X.; Yang, D. A Novel Algorithm for the Fault Diagnosis of a Redundant Inertial Measurement Unit. IEEE Access 2020, 8, 46081–46091. [Google Scholar] [CrossRef]
  13. Shim, D.S.; Yang, C.K. Geometric FDI based on SVD for redundant inertial sensor systems. In Proceedings of the 2004 5th Asian Control Conference, Melbourne, Australia, 20–23 July 2004; pp. 1094–1100. [Google Scholar]
  14. Yang, J.; Zhang, H.Y. Satellite fault diagnosis based on optimal parity vector method. IFAC Proc. 2001, 34, 437–441. [Google Scholar] [CrossRef]
  15. Jing, Y.; Zhang, H. Optimal parity vector method for satellite fault diagnosis. Chin. J. Aeronaut. 2002, 23, 183–186. [Google Scholar]
  16. Li, B.; Zhang, Z. Optimal parity vector technology applied to fault detection of redundant inertial measurement unit. Aerosp. Control 2016, 34, 86–90. [Google Scholar]
  17. Xu, B.; Wen, Z.; Zhang, G.; Li, P. A Novel Fault Detection and Isolation Technology for Redundant Inertial Measurement Unit. IEEE Trans. Instrum. Meas. 2023, 72, 3519711. [Google Scholar] [CrossRef]
  18. Kim, E.J.; Kim, Y.H.; Choi, M.J.; Song, J.W. Optimal Configuration Based on Lever Arm Effect Rejection Performance Index for Redundant Inertial Sensors. IEEE Trans. Instrum. Meas. 2022, 71, 8501111. [Google Scholar] [CrossRef]
  19. Liang, H.; Guo, Y.; Zhao, X. Fault Detection and Isolation of the Multi-Sensor Inertial System. Micromachines 2021, 12, 593. [Google Scholar] [CrossRef]
  20. Zhu, T.; Ren, Y.; Wang, L.; Zhai, X.; Lu, F.; Zou, T. Accuracy Improvement of a Redundant Inertial Measurement Unit Brought about by the Dual-Axis Rotational Motion. In Proceedings of the 2021 International Conference on Computer, Control and Robotics (ICCCR), Shanghai, China, 8–10 January 2021; pp. 138–145. [Google Scholar] [CrossRef]
  21. Cheng, J.; Sun, X.; Liu, P.; Mou, H. An Improved Residual Chi-Square Test Fault Isolation Approach in Four-Gyro SINS. IEEE Access 2019, 7, 174400–174411. [Google Scholar] [CrossRef]
  22. Cheng, J.; Sun, X.; Mou, H. A Modified GLT Double Faults Isolation Approach Based on MLE and RPV for Six-Gyro Redundant SINS. IEEE Access 2019, 7, 5312–5332. [Google Scholar] [CrossRef]
  23. Wang, F.F.; Zhang, T.; Fu, W.X. An improved principal component analysis algorithm on FDI of redundant inertial measurement unit. In Proceedings of the 37th Chinese Control Conference (CCC), Wuhan, China, 25–27 July 2018; pp. 6082–6086. [Google Scholar]
  24. Gai, E.; Harrison, J.V.; Daly, K.C. Generalized likelihood test for FDI in redundant sensor configurations. J. Guid. Control 1979, 2, 9–17. [Google Scholar]
  25. Medvedev, A. Fault detection and isolation by a continuous parity space method. Automatica 1995, 31, 1039–1044. [Google Scholar] [CrossRef]
  26. Pustylnik, E.I. Statistical Analysis Techniques and Observation Processing; Nauka: Moscow, Russia, 1968. [Google Scholar]
  27. Inertial Measurement Units Catalogue 2022 (LLC RPC «Optolink»). Available online: http://optolink.ru/documents/Catalog2022/Optolink_Catalog_2022_p2.pdf (accessed on 25 November 2023).
Figure 1. Matching equations’ acceptable residuals and their SD (for ARSs) (Level I procedure).
Figure 1. Matching equations’ acceptable residuals and their SD (for ARSs) (Level I procedure).
Inventions 08 00159 g001
Figure 2. Matching equations’ acceptable residuals and their SD (for accelerometers) (Level I procedure).
Figure 2. Matching equations’ acceptable residuals and their SD (for accelerometers) (Level I procedure).
Inventions 08 00159 g002
Figure 3. Matching equations’ acceptable residuals and their SD (for ARSs) (Level II procedure).
Figure 3. Matching equations’ acceptable residuals and their SD (for ARSs) (Level II procedure).
Inventions 08 00159 g003
Figure 4. Matching equations’ acceptable residuals and their SD (for accelerometers) (Level II procedure).
Figure 4. Matching equations’ acceptable residuals and their SD (for accelerometers) (Level II procedure).
Inventions 08 00159 g004
Table 1. Acceptable values of residuals and their SD matching equations (ARSs) (Level I procedure).
Table 1. Acceptable values of residuals and their SD matching equations (ARSs) (Level I procedure).
ARSs Systematic Error Value, deg/Hour
1010−110−3
εtolσεtolεtolσεtolεtolσεtol
1.894.010.210.451.76 × 10−33.74 × 10−3
Accelerometers’ Zero Drift Value, m/s2
10−1 g10−3 g10−5 g
εtolσεtolεtolσεtolεtolσεtol
0.370.780.060.350.040.18
Table 2. Acceptable values of residuals and their SD matching equations (ARSs) (Level II procedure).
Table 2. Acceptable values of residuals and their SD matching equations (ARSs) (Level II procedure).
ARSs Systematic Error Value, deg/Hour
1010−110−3
εtolσεtolεtolσεtolεtolσεtol
“fast” cycle (n = 10)
11.6137.750.1150.3741.26 × 10−34.1 × 10−3
“average” cycle (n = 30)
2.1665.9562.1 × 10−25.9 × 10−22.4 × 10−46.5 × 10−4
“slow” cycle (n = 60)
0.73231.8877.3 × 10−31.87 × 10−22.1 × 10−42.1 × 10−4
Table 3. Acceptable values of residuals and their SD matching equations (accelerometers) (Level II procedure).
Table 3. Acceptable values of residuals and their SD matching equations (accelerometers) (Level II procedure).
Accelerometers’ Zero Drift Value, m/s2
10−1g10−3g10−5g
εtolσεtolεtolσεtolεtolσεtol
“fast” cycle (n = 10)
1.1393.70.010.0331.12 × 10−33.65 × 10−3
“average” cycle (n = 30)
0.2120.5840.0025.3 × 10−31.96 × 10−45.43 × 10−4
“slow” cycle (n = 60)
7.18 × 10−20.1856.7 × 10−41.6 × 10−36.16 × 10−51.61 × 10−4
Table 4. Identification of failed measuring elements in an IMU of six sensors.
Table 4. Identification of failed measuring elements in an IMU of six sensors.
Characteristic
of the Unit Structure
Tolerance Value of the Adjustment Error, ang. s.
δni = 1 × 10−4 gδni = 3 × 10−4 gδni = 6 × 10−4 gδni = 8 × 10−4 g
Orthogonal, coinciding with the orts of the
body frame
10306080
Conical χ = 54.74 ° 7.3224458
Conical χ = 78 ° 26.41317
Table 5. Estimation of the matching equations’ residual calculating error for ARSs (deg./hour).
Table 5. Estimation of the matching equations’ residual calculating error for ARSs (deg./hour).
θ m a x , deg. V O , deg/Hour
51020
ME *SDMESDMESD
1/36010−60.4410−60.8810−61.7
1/6010−52.610−55.210−510.5
110−41580.0013160.002633
* ME—mathematical expectation (expected value).
Table 6. Estimation of the matching equations’ residual calculating for accelerometers (m/s2).
Table 6. Estimation of the matching equations’ residual calculating for accelerometers (m/s2).
θ m a x , deg. V O , m/s2
9.812 × 9.813 × 9.81
ME *SDMESDMESD
1/36010−910−410−910−410−910−4
1/6010−910−310−80.0210−80.004
110−70.0810−70.1710−60.25
* ME—mathematical expectation (expected value).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kuznetsov, I.M.; Veremeenko, K.K.; Zharkov, M.V.; Pronkin, A.N. Using Functionally Redundant Inertial Measurement Units to Increase Reliability and Ensure Fault Tolerance. Inventions 2023, 8, 159. https://doi.org/10.3390/inventions8060159

AMA Style

Kuznetsov IM, Veremeenko KK, Zharkov MV, Pronkin AN. Using Functionally Redundant Inertial Measurement Units to Increase Reliability and Ensure Fault Tolerance. Inventions. 2023; 8(6):159. https://doi.org/10.3390/inventions8060159

Chicago/Turabian Style

Kuznetsov, Ivan M., Konstantin K. Veremeenko, Maxim V. Zharkov, and Andrey N. Pronkin. 2023. "Using Functionally Redundant Inertial Measurement Units to Increase Reliability and Ensure Fault Tolerance" Inventions 8, no. 6: 159. https://doi.org/10.3390/inventions8060159

APA Style

Kuznetsov, I. M., Veremeenko, K. K., Zharkov, M. V., & Pronkin, A. N. (2023). Using Functionally Redundant Inertial Measurement Units to Increase Reliability and Ensure Fault Tolerance. Inventions, 8(6), 159. https://doi.org/10.3390/inventions8060159

Article Metrics

Back to TopTop