# Reliability Analysis of Key Components of the Pod Based on Grey System Theory

Authors: Dongfang Hu; Hang Su; Zuhui Shen
Journal: Mathematical Problems in Engineering
 (2020)

---

## Abstract

For improving the reliability of key components in an airborne pod, gray system theory is introduced into the virtual design of airborne pod. Based on this theory, key components of airborne pod are analyzed and mathematical models of the gray relational model and absolute and relative gray relational model are established, respectively. The differences and relations of the three methods for evaluating the robustness of products are researched. Considering the influence of relevant factors, the change rate of each evaluation index is calculated, and the reliability analysis of the airborne pod body is realized ultimately. Finally, the No. 2 test pod body is the optimal solution, which provides an effective theoretical basis for determining the optimal structure of the airborne pod body.


---

## Body

## 1. Introduction

For any mechanical system, in order to make the design, development, and manufacturing reliable, that is, to study its reliability, it is necessary to study the reliability of a mechanical system. Mechanical systems are composed of many factors. For the grey system [1], the reliability of a mechanical engineering system describes the ability of the system to accurately, timely, and coordinately complete specified tasks under specified conditions and within a specified time. When this ability is measured with probability, it is reliability [2]. In the service process of a mechanical system, a reliable mechanical system will play a vital role in guaranteeing its working performance.In recent years, scholars have conducted a series of research and exploration in the field of reliability analysis in view of the diversity and complexity of mechanical structure reliability. The most commonly used methods are probability statistics, fuzzy mathematics, evidence theory, and gray system theory. After consulting the relevant literature, the application research of uncertain systems in various engineering fields has made corresponding progress, for example, combining the analysis method and the checkpoint method with the MATLAB software platform to improve the reliability of key component analysis [3]. The gray system is extended to the theory of probability statistics, the gray probability density function and gray reliability calculation method are studied, and the gray reliability calculation model of stress intensity interference is established [4]. A generalized SPE flux sequence and an extreme value SPE flux sequence taking the occurrence time as a sequence are proposed, and a solar proton event prediction method based on gray GM (1,1) and interval is established [5]. An algorithm based on the self-help method and fuzzy mathematics is proposed, and a fuzzy self-help fuzzy mathematics model based on fusion of bad information objects is established. The measured data of bad information characteristic information objects are self-sampled, and the maximum entropy algorithm is used to obtain the measured data. Self-service distribution of data: self-service fusion sequence uses self-service distribution to estimate the true value and interval of the measured pressure value through the fuzzy membership function [6]. Combining the evidence theory analysis method, an effective model and an efficient and reliable mechanical structure are proposed [7].However, according to the relevant literature review, the existing research content rarely sees the grayness in the structural design of aeronautical products, especially airborne equipment. This article is based on this purpose.


## 2. Establishment of the Mathematical Model

### 2.1. Deng’s Grey System Theory

If the behavior factor of a certain aspect of the mechanical system isX0 and X0 is affected by multiple factors xii=1,2,…,n, then this method of using the gray correlation degree of factor xi to factor X0 to express the influence of xi on X0 is called gray correlation analysis. The method in which X0 affects the size is called the gray relational analysis method. Based on Deng’s gray relational analysis theory, the reliability of the outer frame of the airborne pod is regarded as a gray system, and its structural reliability assessment is determined by various performance indicators.Selecting the structural reliability of the pod cabin as the research object and introducing the concept of behavior factors, the cabin structure evaluation index (one of the reliability evaluation indicators described in the technical parameter file of the airborne pod)X0 system characteristic behavior sequence is(1)X0=x01,x02,…,x0k,…,x0n.In the formula, it is the first evaluation data of a certain behavioral factor.According to the relevant technical requirements of airborne pod, the reliability of the pod structure depends on the interaction of multiple behavior factors. Based on this, a gray relational mathematical model is constructed.Theorem 1.
Assume that there are all kinds of factors affecting the reliability of products. By testing and extracting relevant data, the sequence of behavioral indicators of the firsti≤m factor is obtained as follows:(2)Xi=xi1,xi2,…,xik,…,xin.In the formula,I represents the number of the first influencing factors, i=1,2,…,m;K is the serial number of data, k=1,2,…,n; and xik is the K measured value of influencing factor xi.From formulas (1) and (2), the corresponding forms of the sequence of system characteristic behavior X0 and Xi comparison data are deduced, which are, in turn, as follows:(3)X0=x01,x02,…,x0k,…,x0n,Xi=xi1,xi2,…,xik,…,xin,Xm=xm1,xm2,…,xmk,…,xmn.


### 2.2. Determination of Mechanical Reliability Parameters

Mechanical reliability can generally be divided into structural reliability and mechanical reliability. Structural reliability mainly considers the strength of the mechanical structure and the failure caused by fatigue, wear, and fracture under load; the mechanism reliability does not consider the failure caused by the strength, but the failure caused by the kinematics during the action, so the reliability is not bias.As showed in Figure1, the appearance of the pod body frame is generally spherical. According to the design technical parameters of the pod body, such as the azimuth attitude during flight, the symmetrical shape of aerodynamic characteristics, the extreme environmental requirements at high and low temperatures, the air tightness requirements of the internal environment, the antivibration, and anti-impact requirements, 30 closely related evaluation indicators are determined.Figure 1
Solid modeling of the pod frame.Due to the difference of units between data, there are obvious differences in the order of magnitude of each other. It is obviously necessary to integrate data in dimension. Referring to the evaluation criteria of aeronautical product quality parameters, it is artificially divided into 10 grades: when the evaluation coefficient is 8–10, the index is good, when the evaluation coefficient is 6-7, the index is general, and when the index is 4-5, the index is secondary. Based on 30 performance indicators, the performance of the virtual prototype was tested sequentially. According to the test data, 15 indicators for an evaluation coefficient of 6 were selected as the key factors. According to the given technical requirements document for pod design, the 15 evaluation indexes are sorted from important to general according to the key parts embodied in the design document, as shown in Table1.Table 1
Assessment sequence of reliability parameters for airborne pods.
Serial numberSymboltEchnical parameterEvaluation coefficient (≤10)1x01Strength and stiffness of outer spherical shells7 (good)2x02Strength and stiffness of the inner frame9 (perfect)3x03Antivibration characteristics of outer spherical shells6 (good)4x04Vibration resistance of the inner frame8 (perfect)5x05Symmetry characteristics of outer spherical shells7 (good)6x06Symmetry characteristic of the inner frame7 (good)7x07Gas tightness characteristics of the outer spherical shell8 (perfect)8x08Corrosion resistance of outer spherical shells6 (good)9x09Corrosion resistance of the inner frame8 (perfect)10x010Temperature characteristics of outer spherical shell materials8 (perfect)11x011Temperature characteristics of inner frame materials9 (perfect)12x012Vibration resistance of the azimuth motor7 (good)13x013Strength and stiffness of connectors6 (good)14x014Fatigue resistance of the bearing drive7 (good)15x015Impact resistance of the bearing drive8 (perfect)According to the research emphasis, the structural reliability of the pod body is selected as the research object, and the system characteristic behavior sequence composed of 15 evaluation indexes of the airborne pod body is obtained:(4)X0=7,9,6,8,7,7,8,6,8,8,9,7,6,7,8.In the simulation environment, 10 sets of representative virtual test prototypes are selected from a large number of random test prototypes (i.e., random samples) according to different design ideas. According to the abovementioned performance tests, the evaluation tests of the abovementioned 15 indicators were carried out, respectively. According to the corresponding experimental conditions, 10 sets of comparative data sequences of different results are obtained:(5)X1=7.7,9.1,6.5,8.3,7.4,7.7,8.6,6.2,8.6,8.4,9.3,7.9,6.1,7.4,8.5,X2=7.3,9.4,6.1,8.7,7.6,7.3,8.4,6.5,8.3,8.1,9.5,7.5,6.2,7.1,8.6,X3=7.1,9.3,6.7,8.4,7.3,7.1,8.2,6.8,8.1,8.8,9.3,7.3,6.7,7.3,8.7,X4=7.6,9.2,6.4,8.7,7.2,7.7,8.6,6.1,8.2,8.7,9.1,7.9,6.5,7.4,8.2,X5=7.8,9.1,6.2,8.6,7.7,7.3,8.4,6.4,8.6,8.1,9.5,7.4,6.1,7.3,8.1,X6=7.3,9.5,6.6,8.1,7.3,7.6,8.1,6.7,8.9,8.4,9.3,7.2,6.5,7.4,8.6,X7=7.7,9.2,6.3,8.5,7.7,7.3,8.5,6.3,8.2,8.6,9.4,7.5,6.1,7.3,8.4,X8=7.2,9.5,6.6,8.2,7.4,7.8,8.2,6.6,8.7,8.5,9.2,7.3,6.5,7.1,8.7,X9=7.5,9.3,6.8,8.3,7.6,7.1,8.6,6.1,8.3,8.7,9.1,7.1,6.8,7.3,8.3,X10=7.8,9.2,6.1,8.5,7.3,7.8,8.4,6.9,8.7,8.1,9.3,7.5,6.1,7.7,8.6.Thus, a preliminary data sample for evaluating the reliability of the pod structure is obtained.


## 2.1. Deng’s Grey System Theory

If the behavior factor of a certain aspect of the mechanical system isX0 and X0 is affected by multiple factors xii=1,2,…,n, then this method of using the gray correlation degree of factor xi to factor X0 to express the influence of xi on X0 is called gray correlation analysis. The method in which X0 affects the size is called the gray relational analysis method. Based on Deng’s gray relational analysis theory, the reliability of the outer frame of the airborne pod is regarded as a gray system, and its structural reliability assessment is determined by various performance indicators.Selecting the structural reliability of the pod cabin as the research object and introducing the concept of behavior factors, the cabin structure evaluation index (one of the reliability evaluation indicators described in the technical parameter file of the airborne pod)X0 system characteristic behavior sequence is(1)X0=x01,x02,…,x0k,…,x0n.In the formula, it is the first evaluation data of a certain behavioral factor.According to the relevant technical requirements of airborne pod, the reliability of the pod structure depends on the interaction of multiple behavior factors. Based on this, a gray relational mathematical model is constructed.Theorem 1.
Assume that there are all kinds of factors affecting the reliability of products. By testing and extracting relevant data, the sequence of behavioral indicators of the firsti≤m factor is obtained as follows:(2)Xi=xi1,xi2,…,xik,…,xin.In the formula,I represents the number of the first influencing factors, i=1,2,…,m;K is the serial number of data, k=1,2,…,n; and xik is the K measured value of influencing factor xi.From formulas (1) and (2), the corresponding forms of the sequence of system characteristic behavior X0 and Xi comparison data are deduced, which are, in turn, as follows:(3)X0=x01,x02,…,x0k,…,x0n,Xi=xi1,xi2,…,xik,…,xin,Xm=xm1,xm2,…,xmk,…,xmn.


## 2.2. Determination of Mechanical Reliability Parameters

Mechanical reliability can generally be divided into structural reliability and mechanical reliability. Structural reliability mainly considers the strength of the mechanical structure and the failure caused by fatigue, wear, and fracture under load; the mechanism reliability does not consider the failure caused by the strength, but the failure caused by the kinematics during the action, so the reliability is not bias.As showed in Figure1, the appearance of the pod body frame is generally spherical. According to the design technical parameters of the pod body, such as the azimuth attitude during flight, the symmetrical shape of aerodynamic characteristics, the extreme environmental requirements at high and low temperatures, the air tightness requirements of the internal environment, the antivibration, and anti-impact requirements, 30 closely related evaluation indicators are determined.Figure 1
Solid modeling of the pod frame.Due to the difference of units between data, there are obvious differences in the order of magnitude of each other. It is obviously necessary to integrate data in dimension. Referring to the evaluation criteria of aeronautical product quality parameters, it is artificially divided into 10 grades: when the evaluation coefficient is 8–10, the index is good, when the evaluation coefficient is 6-7, the index is general, and when the index is 4-5, the index is secondary. Based on 30 performance indicators, the performance of the virtual prototype was tested sequentially. According to the test data, 15 indicators for an evaluation coefficient of 6 were selected as the key factors. According to the given technical requirements document for pod design, the 15 evaluation indexes are sorted from important to general according to the key parts embodied in the design document, as shown in Table1.Table 1
Assessment sequence of reliability parameters for airborne pods.
Serial numberSymboltEchnical parameterEvaluation coefficient (≤10)1x01Strength and stiffness of outer spherical shells7 (good)2x02Strength and stiffness of the inner frame9 (perfect)3x03Antivibration characteristics of outer spherical shells6 (good)4x04Vibration resistance of the inner frame8 (perfect)5x05Symmetry characteristics of outer spherical shells7 (good)6x06Symmetry characteristic of the inner frame7 (good)7x07Gas tightness characteristics of the outer spherical shell8 (perfect)8x08Corrosion resistance of outer spherical shells6 (good)9x09Corrosion resistance of the inner frame8 (perfect)10x010Temperature characteristics of outer spherical shell materials8 (perfect)11x011Temperature characteristics of inner frame materials9 (perfect)12x012Vibration resistance of the azimuth motor7 (good)13x013Strength and stiffness of connectors6 (good)14x014Fatigue resistance of the bearing drive7 (good)15x015Impact resistance of the bearing drive8 (perfect)According to the research emphasis, the structural reliability of the pod body is selected as the research object, and the system characteristic behavior sequence composed of 15 evaluation indexes of the airborne pod body is obtained:(4)X0=7,9,6,8,7,7,8,6,8,8,9,7,6,7,8.In the simulation environment, 10 sets of representative virtual test prototypes are selected from a large number of random test prototypes (i.e., random samples) according to different design ideas. According to the abovementioned performance tests, the evaluation tests of the abovementioned 15 indicators were carried out, respectively. According to the corresponding experimental conditions, 10 sets of comparative data sequences of different results are obtained:(5)X1=7.7,9.1,6.5,8.3,7.4,7.7,8.6,6.2,8.6,8.4,9.3,7.9,6.1,7.4,8.5,X2=7.3,9.4,6.1,8.7,7.6,7.3,8.4,6.5,8.3,8.1,9.5,7.5,6.2,7.1,8.6,X3=7.1,9.3,6.7,8.4,7.3,7.1,8.2,6.8,8.1,8.8,9.3,7.3,6.7,7.3,8.7,X4=7.6,9.2,6.4,8.7,7.2,7.7,8.6,6.1,8.2,8.7,9.1,7.9,6.5,7.4,8.2,X5=7.8,9.1,6.2,8.6,7.7,7.3,8.4,6.4,8.6,8.1,9.5,7.4,6.1,7.3,8.1,X6=7.3,9.5,6.6,8.1,7.3,7.6,8.1,6.7,8.9,8.4,9.3,7.2,6.5,7.4,8.6,X7=7.7,9.2,6.3,8.5,7.7,7.3,8.5,6.3,8.2,8.6,9.4,7.5,6.1,7.3,8.4,X8=7.2,9.5,6.6,8.2,7.4,7.8,8.2,6.6,8.7,8.5,9.2,7.3,6.5,7.1,8.7,X9=7.5,9.3,6.8,8.3,7.6,7.1,8.6,6.1,8.3,8.7,9.1,7.1,6.8,7.3,8.3,X10=7.8,9.2,6.1,8.5,7.3,7.8,8.4,6.9,8.7,8.1,9.3,7.5,6.1,7.7,8.6.Thus, a preliminary data sample for evaluating the reliability of the pod structure is obtained.


## 3. Establishment of Dun’s Grey Correlation Model

### 3.1. Interval Value Processing

Based on the correlation degree problem in the gray system theory, it is known that the reliability analysis is mainly to characterize the similarity degree of the system characteristic behavior line graphs of each comparison data line, and the representation index is exactly 15 points, selected according to the point coordinates. The similarity on the graph is used to estimate the result. Then, the evaluation range 1∼10 will be the key to the interval value, and the interval value method will be selected for processing. The interval-valued operator will be introduced to obtain the following relationship:(6)Xi=xi1,xi2,…,xik,…,xin,XiD=xi1d,xi2d,…,xikd,…,xind,(7)xikd=xik−minxikmaxxik−minxik,k=1,2,…,n.Then, according to the interval operator quoted by equations (6) and (7), 11 sets of different data (X0 is the feature sequence) X0,X1,X2,…,X10 are interval-valued, and the gray system algorithm program is written by Microsoft Visual C++ software. The interval-valued data sequence is(8)X0′=0.33,1.00,0.00,0.67,0.33,0.33,0.67,0.00,0.67,0.67,1.00,0.33,0.00,0.33,0.67,X1′=0.50,0.94,0.13,0.69,0.41,0.50,0.78,0.03,0.78,0.72,1.00,0.56,0.00,0.41,0.75,X2′=0.35,0.97,0.00,0.76,0.44,0.35,0.68,0.12,0.65,0.59,1.00,0.41,0.03,0.29,0.74,X3′=0.15,1.00,0.00,0.65,0.23,0.15,0.58,0.04,0.54,0.81,1.00,0.23,0.00,0.23,0.77,X4′=0.48,1.00,0.10,0.84,0.35,0.52,0.81,0.00,0.68,0.84,0.97,0.58,0.13,0.42,0.68,X5′=0.50,0.88,0.03,0.74,0.47,0.35,0.68,0.09,0.74,0.59,1.00,0.38,0.00,0.35,0.59,X6′=0.27,1.00,0.03,0.53,0.27,0.37,0.53,0.07,0.80,0.63,0.93,0.23,0.00,0.30,0.70,X7′=0.48,0.94,0.06,0.73,0.48,0.36,0.73,0.06,0.64,0.76,1.00,0.42,0.00,0.36,0.70,X8′=0.23,1.00,0.03,0.57,0.30,0.43,0.57,0.03,0.73,0.67,0.90,0.27,0.00,0.20,0.73,X9′=0.44,1.00,0.22,0.69,0.47,0.31,0.78,0.00,0.69,0.81,0.94,0.31,0.22,0.38,0.69,X10′=0.53,0.97,0.00,0.75,0.38,0.53,0.72,0.25,0.81,0.63,1.00,0.44,0.00,0.50,0.78.


### 3.2. Difference Sequence and Range

The 11 sets of new sequences derived after interval-valued are, respectively, recorded asX0′,X1′,X2′,…,X10′, and the following relations are derived from the four axioms of gray correlation [8]:(9)Δik=x0′k−xi′k,(10)Δi=Δi1,Δi2,…,Δin,i=0,1,2,…,m.According to equations (9) and (10), the difference data sequence (retaining 2 significant digits) is obtained, which essentially calculates the geometric distance between the corresponding two points:(11)Δ1=0.17,0.06,0.13,0.02,0.08,0.17,0.11,0.03,0.11,0.05,0.00,0.23,0.00,0.08,0.08,Δ2=0.02,0.03,0.00,0.09,0.11,0.02,0.01,0.12,0.02,0.08,0.00,0.08,0.03,0.04,0.07,Δ3=0.18,0.00,0.00,0.02,0.10,0.18,0.09,0.04,0.13,0.14,0.00,0.10,0.00,0.10,0.10,Δ4=0.15,0.00,0.10,0.17,0.02,0.19,0.14,0.00,0.01,0.17,0.03,0.25,0.13,0.09,0.01,Δ5=0.17,0.12,0.03,0.07,0.14,0.02,0.01,0.09,0.07,0.08,0.00,0.05,0.00,0.02,0.08,Δ6=0.06,0.00,0.03,0.14,0.06,0.04,0.14,0.07,0.13,0.04,0.07,0.10,0.00,0.03,0.03,Δ7=0.15,0.06,0.06,0.06,0.15,0.03,0.06,0.06,0.03,0.09,0.00,0.09,0.00,0.03,0.03,Δ8=0.10,0.00,0.03,0.10,0.03,0.10,0.10,0.03,0.06.0.06,0.10,0.06,0.00,0.13,0.06,Δ9=0.11,0.00,0.22,0.02,0.14,0.02,0.11,0.00,0.02,0.14,0.06,0.02,0.22,0.05,0.02,Δ10=0.20,0.03,0.00,0.08,0.05,0.20,0.05,0.25,0.14,0.04,0.00,0.11,0.00,0.17,0.11.For each set of different sequences, the data are retrieved according to the written C++ program, and the maximum and minimum values are, respectively, recorded as(12)M=maximaxkΔik,(13)m=miniminkΔik.According to equations (12) and (13), it is necessary to find the extreme difference value, wherein λ0M,m corresponding to the first group Δ0 is the maximum value and the minimum value of the corresponding sequence of the group, and the rest is denoted as λ1,λ2,…,λ10; then, there are(14)λ1=0.23,0,λ2=0.12,0,λ3=0.18,0,λ4=0.25,0,λ5=0.17,0,λ6=0.14,0,λ7=0.15,0,λ8=0.13,0,λ9=0.22,0,λ10=0.25,0.From each limit value, find the extreme difference between the maximum value and the minimum value:(15)M=0.25,m=0.


### 3.3. Correlation Coefficient Sequence and Relevance

Definition 1.
The system behavior sequence and the comparison data sequence are as shown in equation (3), defining the resolution coefficient δ=γ09,γ01,γ010,γ07,γ03,γ08,γ05,γ04,γ06,γ02∈0,1; let(16)γx0k,xik=miniminkx0k−xik+ζmaximaxkx0k−xikx0k−xik+ζmaximaxkx0k−xik,(17)γX0,Xi=1n∑k=1nγx0k,xik.
The given real numberγx0k,xik constitutes a sequence of gray correlation coefficients. If the coefficient γX0,Xi satisfies the gray correlation four axioms, it is called the gray correlation degree of Xi to X0. According to formulas (12), (13), (16), and (17), the following operation is performed, and the formula of the correlation coefficient is rewritten as(18)γ0ik=m+ζMΔik+ζM,ζ∈0,1,k=1,2,…,n;i=1,2,…,m.
According to formula (18), the gray correlation coefficient sequence and Deng’s gray correlation degree are derived:(19)β0i=γ0i1,γ0i2,…,γoik,…,γ0in,γ0i=1n∑k=1nγ0ik;i=1,2,…,m.
Since the weights of each sequence data are similar, the resolution coefficientζ is 0.5, and the abovementioned algorithm is implemented by C++ program to obtain the corresponding data sequence:(20)β01=0.42,0.68,0.49,0.86,0.61,0.42,0.53,0.81,0.53,0.71,1.00,0.35,1.00,0.61,0.61,β02=0.86,0.81,1.00,0.58,0.53,0.86,0.93,0.51,0.86,0.64,1.00,0.64,0.81,0.76,0.64,β03=0.41,1.00,1.00,0.86,0.56,0.41,0.58,0.76,0.49,0.47,1.00,0.56,1.00,0.56,0.56,β04=0.45,1.00,0.56,0.42,0.86,0.40,0.47,1.00,0.93,0.42,0.81,0.33,0.49,0.58,0.93,β05=0.42,0.51,0.81,0.64,0.47,0.86,0.93,0.58,0.64,0.61,1.00,0.71,1.00,0.86,0.61,β06=0.68,1.00,0.81,0.47,0.68,0.76,0.47,0.64,0.49,0.76,0.64,0.56,1.00,0.81,0.81,β07=0.45,0.68,0.68,0.68,0.45,0.81,0.68,0.68,0.81,0.58,1.00,0.58,1.00,0.81,0.81,β08=0.56,1.00,0.81,0.56,0.81,0.56,0.56,0.81,0.68,0.68,0.56,0.68,1.00,0.49,0.68,β09=0.53,1.00,0.36,0.86,0.47,0.86,0.53,1.00,0.86,0.47,0.68,0.86,0.36,0.71,0.86,β10=0.38,0.81,1.00,0.61,0.71,0.38,0.71,0.33,0.47,0.76,1.00,0.53,1.00,0.42,0.53.
The gray correlation degrees corresponding toX0 are(21)γ01=0.642,γ02=0.762,γ03=0.681,γ04=0.643,γ05=0.710,γ06=0.705,γ07=0.713,γ08=0.696,γ09=0.694,γ10=0.643.
In summary,γ02 has the greatest correlation with the behavioral data sequence X0 in the data sequence of 10 sets of virtual prototype performance tests, that is, the test index of the test cabin No. 2 is the most consistent with the assessment criteria.
In order to visually explain the accuracy of the gray correlation analysis, 10 sets of data are, respectively, made into data line diagrams composed of 10 sets of data sequences, as shown in Figure2.
Figure2 shows that the 15 performance indicators extracted from the cabin reliability analysis correspond to 15 coordinate points, respectively. According to the Deng’s gray correlation analysis model, the closeness of the data lines X1,X2,…,X10 and X0 is compared by 15 points. The value is judged by the proximity of the standard value (the corresponding value of the behavioral feature sequence). Therefore, Deng’s gray correlation method is also called the “point correlation method” [9]. Through theoretical calculations, the gray relational order is obtained:(22)δ=γ02,γ07,γ05,γ06,γ08,γ09,γ03,γ04,γ10,γ01.
In Figure2, the correlation between the corresponding 10 sets of comparison data sequence lines and the standard characteristic data sequence lines is consistent with the “point correlation method” analysis trend. According to the gray correlation sequence, the case of the virtual prototype cabins and standard characteristic data sequences of No. 2 and No. 7 is investigated.
As shown in Figure3(a), it can be seen from the comparison between the data display and the test sample cabin No. 1 in Figure 3(b) that although X02 has the greatest correlation with X0, X02 is not the best in shape similar to X0. X0 set of data line graphs, such as X07, has a better degree of similarity. The so-called similarity is the degree to which the graph is consistent with the shape of the standard line graph E. Since the cabin has an economic cost as one of its design requirements, the similarity and proximity of the standard line should be considered separately. Obviously, C is more consistent with B in the similar aspect, but the point correlation does not clearly distinguish this. From the perspective of relevance, as shown in Figure 4, the higher the value, the better the relevance of the indicator.
In the figure, the 15-point correlation index betweenβ02 and β07 is distributed between 0.4 and 1.0. If 8 is taken as the baseline, the virtual prototype of the 7th has 6 indicators that meet the requirements, while the virtual prototype 2 has eight, and if 0.9 is used as the baseline, there are 2 on the 7th and 3 indicators on the 2nd, so the performance of the No. 7 virtual prototype is slightly worse.
The abovementioned results are consistent with Deng’s related analysis results. The main reason is that Deng’s gray correlation analysis model focuses on the similarity of the two curves on the relative displacement of several points on the curve, and the specific analysis does not reflect the characteristics of the system. The relative change rate of the curve and the comparison curve, that is, formula (16), contains only local parameters that reflect the closeness of the two curves and lacks the overall parameters that reflect the closeness of the two curves. The method is judged by the closeness of the limited point, and it cannot fully reflect the specific gray-related requirements such as economic factors, product demand factors, and social factors of the entire mechanical product research object. Therefore, a correlation analysis model based on a macroperspective is needed to ensure that the similarity of similar prototypes is not excluded. The introduction of generalized gray correlation analysis is to consider the relative rate of change factors. The meaning is here.
The biggest innovation in the generalized gray relational analysis system is that the similarity and proximity of the curves can be discussed separately for different engineering actual requirements.Figure 2
Data plot of the performance test of the pod cabin.Figure 3
Data line diagram of the test cabin performance test for No. 2 and No. 7.
(a)(b)Figure 4
Analysis of the test performance of test cabins No. 2 and No. 7.


## 3.1. Interval Value Processing

Based on the correlation degree problem in the gray system theory, it is known that the reliability analysis is mainly to characterize the similarity degree of the system characteristic behavior line graphs of each comparison data line, and the representation index is exactly 15 points, selected according to the point coordinates. The similarity on the graph is used to estimate the result. Then, the evaluation range 1∼10 will be the key to the interval value, and the interval value method will be selected for processing. The interval-valued operator will be introduced to obtain the following relationship:(6)Xi=xi1,xi2,…,xik,…,xin,XiD=xi1d,xi2d,…,xikd,…,xind,(7)xikd=xik−minxikmaxxik−minxik,k=1,2,…,n.Then, according to the interval operator quoted by equations (6) and (7), 11 sets of different data (X0 is the feature sequence) X0,X1,X2,…,X10 are interval-valued, and the gray system algorithm program is written by Microsoft Visual C++ software. The interval-valued data sequence is(8)X0′=0.33,1.00,0.00,0.67,0.33,0.33,0.67,0.00,0.67,0.67,1.00,0.33,0.00,0.33,0.67,X1′=0.50,0.94,0.13,0.69,0.41,0.50,0.78,0.03,0.78,0.72,1.00,0.56,0.00,0.41,0.75,X2′=0.35,0.97,0.00,0.76,0.44,0.35,0.68,0.12,0.65,0.59,1.00,0.41,0.03,0.29,0.74,X3′=0.15,1.00,0.00,0.65,0.23,0.15,0.58,0.04,0.54,0.81,1.00,0.23,0.00,0.23,0.77,X4′=0.48,1.00,0.10,0.84,0.35,0.52,0.81,0.00,0.68,0.84,0.97,0.58,0.13,0.42,0.68,X5′=0.50,0.88,0.03,0.74,0.47,0.35,0.68,0.09,0.74,0.59,1.00,0.38,0.00,0.35,0.59,X6′=0.27,1.00,0.03,0.53,0.27,0.37,0.53,0.07,0.80,0.63,0.93,0.23,0.00,0.30,0.70,X7′=0.48,0.94,0.06,0.73,0.48,0.36,0.73,0.06,0.64,0.76,1.00,0.42,0.00,0.36,0.70,X8′=0.23,1.00,0.03,0.57,0.30,0.43,0.57,0.03,0.73,0.67,0.90,0.27,0.00,0.20,0.73,X9′=0.44,1.00,0.22,0.69,0.47,0.31,0.78,0.00,0.69,0.81,0.94,0.31,0.22,0.38,0.69,X10′=0.53,0.97,0.00,0.75,0.38,0.53,0.72,0.25,0.81,0.63,1.00,0.44,0.00,0.50,0.78.


## 3.2. Difference Sequence and Range

The 11 sets of new sequences derived after interval-valued are, respectively, recorded asX0′,X1′,X2′,…,X10′, and the following relations are derived from the four axioms of gray correlation [8]:(9)Δik=x0′k−xi′k,(10)Δi=Δi1,Δi2,…,Δin,i=0,1,2,…,m.According to equations (9) and (10), the difference data sequence (retaining 2 significant digits) is obtained, which essentially calculates the geometric distance between the corresponding two points:(11)Δ1=0.17,0.06,0.13,0.02,0.08,0.17,0.11,0.03,0.11,0.05,0.00,0.23,0.00,0.08,0.08,Δ2=0.02,0.03,0.00,0.09,0.11,0.02,0.01,0.12,0.02,0.08,0.00,0.08,0.03,0.04,0.07,Δ3=0.18,0.00,0.00,0.02,0.10,0.18,0.09,0.04,0.13,0.14,0.00,0.10,0.00,0.10,0.10,Δ4=0.15,0.00,0.10,0.17,0.02,0.19,0.14,0.00,0.01,0.17,0.03,0.25,0.13,0.09,0.01,Δ5=0.17,0.12,0.03,0.07,0.14,0.02,0.01,0.09,0.07,0.08,0.00,0.05,0.00,0.02,0.08,Δ6=0.06,0.00,0.03,0.14,0.06,0.04,0.14,0.07,0.13,0.04,0.07,0.10,0.00,0.03,0.03,Δ7=0.15,0.06,0.06,0.06,0.15,0.03,0.06,0.06,0.03,0.09,0.00,0.09,0.00,0.03,0.03,Δ8=0.10,0.00,0.03,0.10,0.03,0.10,0.10,0.03,0.06.0.06,0.10,0.06,0.00,0.13,0.06,Δ9=0.11,0.00,0.22,0.02,0.14,0.02,0.11,0.00,0.02,0.14,0.06,0.02,0.22,0.05,0.02,Δ10=0.20,0.03,0.00,0.08,0.05,0.20,0.05,0.25,0.14,0.04,0.00,0.11,0.00,0.17,0.11.For each set of different sequences, the data are retrieved according to the written C++ program, and the maximum and minimum values are, respectively, recorded as(12)M=maximaxkΔik,(13)m=miniminkΔik.According to equations (12) and (13), it is necessary to find the extreme difference value, wherein λ0M,m corresponding to the first group Δ0 is the maximum value and the minimum value of the corresponding sequence of the group, and the rest is denoted as λ1,λ2,…,λ10; then, there are(14)λ1=0.23,0,λ2=0.12,0,λ3=0.18,0,λ4=0.25,0,λ5=0.17,0,λ6=0.14,0,λ7=0.15,0,λ8=0.13,0,λ9=0.22,0,λ10=0.25,0.From each limit value, find the extreme difference between the maximum value and the minimum value:(15)M=0.25,m=0.


## 3.3. Correlation Coefficient Sequence and Relevance

Definition 1.
The system behavior sequence and the comparison data sequence are as shown in equation (3), defining the resolution coefficient δ=γ09,γ01,γ010,γ07,γ03,γ08,γ05,γ04,γ06,γ02∈0,1; let(16)γx0k,xik=miniminkx0k−xik+ζmaximaxkx0k−xikx0k−xik+ζmaximaxkx0k−xik,(17)γX0,Xi=1n∑k=1nγx0k,xik.
The given real numberγx0k,xik constitutes a sequence of gray correlation coefficients. If the coefficient γX0,Xi satisfies the gray correlation four axioms, it is called the gray correlation degree of Xi to X0. According to formulas (12), (13), (16), and (17), the following operation is performed, and the formula of the correlation coefficient is rewritten as(18)γ0ik=m+ζMΔik+ζM,ζ∈0,1,k=1,2,…,n;i=1,2,…,m.
According to formula (18), the gray correlation coefficient sequence and Deng’s gray correlation degree are derived:(19)β0i=γ0i1,γ0i2,…,γoik,…,γ0in,γ0i=1n∑k=1nγ0ik;i=1,2,…,m.
Since the weights of each sequence data are similar, the resolution coefficientζ is 0.5, and the abovementioned algorithm is implemented by C++ program to obtain the corresponding data sequence:(20)β01=0.42,0.68,0.49,0.86,0.61,0.42,0.53,0.81,0.53,0.71,1.00,0.35,1.00,0.61,0.61,β02=0.86,0.81,1.00,0.58,0.53,0.86,0.93,0.51,0.86,0.64,1.00,0.64,0.81,0.76,0.64,β03=0.41,1.00,1.00,0.86,0.56,0.41,0.58,0.76,0.49,0.47,1.00,0.56,1.00,0.56,0.56,β04=0.45,1.00,0.56,0.42,0.86,0.40,0.47,1.00,0.93,0.42,0.81,0.33,0.49,0.58,0.93,β05=0.42,0.51,0.81,0.64,0.47,0.86,0.93,0.58,0.64,0.61,1.00,0.71,1.00,0.86,0.61,β06=0.68,1.00,0.81,0.47,0.68,0.76,0.47,0.64,0.49,0.76,0.64,0.56,1.00,0.81,0.81,β07=0.45,0.68,0.68,0.68,0.45,0.81,0.68,0.68,0.81,0.58,1.00,0.58,1.00,0.81,0.81,β08=0.56,1.00,0.81,0.56,0.81,0.56,0.56,0.81,0.68,0.68,0.56,0.68,1.00,0.49,0.68,β09=0.53,1.00,0.36,0.86,0.47,0.86,0.53,1.00,0.86,0.47,0.68,0.86,0.36,0.71,0.86,β10=0.38,0.81,1.00,0.61,0.71,0.38,0.71,0.33,0.47,0.76,1.00,0.53,1.00,0.42,0.53.
The gray correlation degrees corresponding toX0 are(21)γ01=0.642,γ02=0.762,γ03=0.681,γ04=0.643,γ05=0.710,γ06=0.705,γ07=0.713,γ08=0.696,γ09=0.694,γ10=0.643.
In summary,γ02 has the greatest correlation with the behavioral data sequence X0 in the data sequence of 10 sets of virtual prototype performance tests, that is, the test index of the test cabin No. 2 is the most consistent with the assessment criteria.
In order to visually explain the accuracy of the gray correlation analysis, 10 sets of data are, respectively, made into data line diagrams composed of 10 sets of data sequences, as shown in Figure2.
Figure2 shows that the 15 performance indicators extracted from the cabin reliability analysis correspond to 15 coordinate points, respectively. According to the Deng’s gray correlation analysis model, the closeness of the data lines X1,X2,…,X10 and X0 is compared by 15 points. The value is judged by the proximity of the standard value (the corresponding value of the behavioral feature sequence). Therefore, Deng’s gray correlation method is also called the “point correlation method” [9]. Through theoretical calculations, the gray relational order is obtained:(22)δ=γ02,γ07,γ05,γ06,γ08,γ09,γ03,γ04,γ10,γ01.
In Figure2, the correlation between the corresponding 10 sets of comparison data sequence lines and the standard characteristic data sequence lines is consistent with the “point correlation method” analysis trend. According to the gray correlation sequence, the case of the virtual prototype cabins and standard characteristic data sequences of No. 2 and No. 7 is investigated.
As shown in Figure3(a), it can be seen from the comparison between the data display and the test sample cabin No. 1 in Figure 3(b) that although X02 has the greatest correlation with X0, X02 is not the best in shape similar to X0. X0 set of data line graphs, such as X07, has a better degree of similarity. The so-called similarity is the degree to which the graph is consistent with the shape of the standard line graph E. Since the cabin has an economic cost as one of its design requirements, the similarity and proximity of the standard line should be considered separately. Obviously, C is more consistent with B in the similar aspect, but the point correlation does not clearly distinguish this. From the perspective of relevance, as shown in Figure 4, the higher the value, the better the relevance of the indicator.
In the figure, the 15-point correlation index betweenβ02 and β07 is distributed between 0.4 and 1.0. If 8 is taken as the baseline, the virtual prototype of the 7th has 6 indicators that meet the requirements, while the virtual prototype 2 has eight, and if 0.9 is used as the baseline, there are 2 on the 7th and 3 indicators on the 2nd, so the performance of the No. 7 virtual prototype is slightly worse.
The abovementioned results are consistent with Deng’s related analysis results. The main reason is that Deng’s gray correlation analysis model focuses on the similarity of the two curves on the relative displacement of several points on the curve, and the specific analysis does not reflect the characteristics of the system. The relative change rate of the curve and the comparison curve, that is, formula (16), contains only local parameters that reflect the closeness of the two curves and lacks the overall parameters that reflect the closeness of the two curves. The method is judged by the closeness of the limited point, and it cannot fully reflect the specific gray-related requirements such as economic factors, product demand factors, and social factors of the entire mechanical product research object. Therefore, a correlation analysis model based on a macroperspective is needed to ensure that the similarity of similar prototypes is not excluded. The introduction of generalized gray correlation analysis is to consider the relative rate of change factors. The meaning is here.
The biggest innovation in the generalized gray relational analysis system is that the similarity and proximity of the curves can be discussed separately for different engineering actual requirements.Figure 2
Data plot of the performance test of the pod cabin.Figure 3
Data line diagram of the test cabin performance test for No. 2 and No. 7.
(a)(b)Figure 4
Analysis of the test performance of test cabins No. 2 and No. 7.


## 4. Generalized Grey Relational Analysis Model

### 4.1. Generalized Grey Relational Method

The generalized gray relational analysis model has three objectives for the analysis of pod reliability: first, the curve approach of Deng’s gray correlation analysis further demonstrates the basis of reliability optimal solution selection from a new perspective; second, the data sequence curve. The relative change rate problem between the two is added to the gray correlation analysis, and the original “point correlation analysis” is changed to the “face correlation analysis” between the curves. Starting from the overall situation of the curve distribution, the two curves are obtained by the starting point of zeroing. They are all placed on the coordinate axis, and the size of the area sandwiched between them is used as the criterion for judging the optimal solution of reliability analysis from the macroscopic perspective of the data. Third, the difference between the proximity and similarity of the curve is analyzed accordingly. Also, in the engineering practice, the corresponding factors of economic and design requirements are added. The generalized gray relational analysis is divided into three categories: absolute association method, relative association method, and comprehensive association method, which do not strictly satisfy the gray association four axioms.


### 4.2. Grey Absolute Correlation Method

Theorem 2.
Set the behavior sequenceXi=xi1,xi2,…,xin and record the corresponding polyline Xi−xi1 as Xi−xi1=xi1−xi1,xi2−xi1,…,xin−xi1; then,(23)si=∫1nXi−xi1dt.Definition 2.
Let the behavior sequenceXi=xi1,xi2,…,xin and D be defined as sequence operators and have the following relationship:(24)XiD=xi1d,xi2d,…,xind,(25)xikd=xik=xi1,k=1,2,…,n.
According to equations (24) and (25), D is called the starting point zeroing operator, and the starting point of XiD is zeroed, which is recorded as(26)Xi0=xi01,xi02,…,xi0n.Theorem 3.
We arbitrarily take two sequences of behavior:(27)Xi=xi1,xi2,…,xin,Xj=xj1,xj2,…,xjn.The starting point zero image is(28)Xi0=xi01,xi02,…,xi0n,Xj0=xj01,xj02,…,xj0n.It is referred to as(29)si−sj=∫1nXi0−Xj0dt.IfXi0 is always above Xj0, then si−sj≥0; if Xi0 is always below Xj0, then si−sj≤0; and if Xi0 intersects with Xj0, the si−sj sign is indefinite.Definition 3.
The sequence ofXi0: the sum of the time intervals between the observation data is the length of Xj0. It should be pointed out that the number of observations in two sequences of equal length is not necessarily equal.Theorem 4.
Let sequenceX0 and Xi be of the same length, s0 and si are as described in equation (15) of Theorem 2, and si−s0 is as described in equation (19) of Theorem 2; then,(30)ε0i=1+s0+si1+s0+si+s0−si.Equation (30) is the absolute degree of gray correlation between X0 and Xi.According to the satisfaction condition, for 11 sets of data sequences, since the time intervals between adjacent observation data are the same,X0 and X1,X2,…,X10 are said to constitute an isochronous sequence.Theorem 5.
Let sequenceX0 and Xi be the same length, and both are 1-time-distance sequences, with the corresponding relationship:(31)X00=x001,x002,…,x00n,Xi0=xi01,xi02,…,xi0n.The abovementioned sequence is the starting point zeroing image ofX0 and Xi, respectively, and the following correspondence is given:(32)s0=∑k=2n−1x00k+12x00n,(33)si=∑k=2n−1xi0k+12xi0n,(34)s0−si=∑k=2n−1xi0k−x00k+12xi0n−x00n.From equations (32)–(34), the following absolute correlation algorithm is derived:(35)ε0i=1+∑k=2n−1x00k+1/2x00n+∑k=2n−1xi0k+1/2xi0n1+∑k=2n−1x00k+1/2x00n+∑k=2n−1xi0k+1/2xi0n+∑k=2n−1xi0k−x00k+1/2xi0n−x00n.According to the gray absolute relevance theory, the related algorithm is written based on the C++ programming language, and the starting point zeroing process is performed according to the data sequence obtained byX0 and X1,X2,…,X10:(36)X00=0.00,2.00,−1.00,1.00,0.00,0.00,1.00,−1.00,1.00,1.00,2.00,0.00,−1.00,0.00,1.00,X10=0.00,1.40,−1.20,0.60,−0.30,0.00,0.90,−1.50,0.90,0.70,1.60,0.20,−1.60,−0.30,0.80,X20=0.00,2.10,−1.20,1.40,0.30,0.00,1.10,−0.80,1.00,0.80,2.20,0.20,−1.10,−0.20,1.30,X30=0.00,2.20,−0.40,1.30,0.20,0.00,1.10,−0.30,1.00,1.70,2.20,0.20,−0.40,−0.20,1.60,X40=0.00,1.60,−1.20,1.10,−0.40,0.10,1.00,−1.60,0.60,1.10,1.50,0.30,−1.10,−0.20,0.60,X50=0.00,1.30,−1.60,0.80,−0.10,−0.50,0.60,−1.40,0.80,0.30,1.70,−0.40,−1.70,−0.50,0.30,X60=0.00,2.20,−0.70,0.80,0.00,0.30,0.80,−0.60,1.60,1.10,2.00,−0.10,−0.80,0.10,1.30,X70=0.00,1.50,−1.40,0.80,0.00,−0.40,0.80,−1.40,0.50,0.90,1.70,−0.20,−1.60,−0.40,0.70,X80=0.00,2.30,−0.60,1.00,0.20,0.60,1.00,−0.60,1.50,1.30,2.00,0.10,−0.70,−0.10,1.50,X90=0.00,1.80,−0.70,0.80,0.10,−0.40,1.10,−1.40,0.80,1.20,1.60,−0.40,−0.70,−0.20,0.80,X100=0.00,1.40,−1.70,0.70,−0.50,0.00,0.60,−0.90,0.90,0.30,1.50,−0.30,−1.70,−0.10,0.80.The corresponding absolute relevance is(37)ε01=0.691,ε02=0.932,ε03=0.803,ε04=0.80,ε05=0.588,ε06=0.882,ε07=0.638,ε08=0.824,ε09=0.875,ε010=0.592.According to the theory of gray absolute correlation analysis, the corresponding C++ program is written to calculate the gray absolute correlation order, which is discharged from left to right according to the degree of relevance:(38)τ=ε02,ε06,ε09,ε08,ε03,ε04,ε01,ε07,ε010,ε05.It can be seen from the results of the absolute correlation degree that the results based on the “face correlation method” are numerically different from those based on the “point correlation method,” mainly because the considerations are different because the corresponding curve ofX1,X2,…,X10 and X0 to the absolute degree of correlation reflects the similarity between the curves.Specifically, to evaluate the reliability of the test prototype, as shown in Figures5(a) and 5(b), the area between the data line X2 and the data line X0 in (a) is smaller to the data line X6 in Figure 5(b). The area between the data lines X0 is essentially reflecting the similarity between the curves. After the starting point is zeroed, the starting points are classified into the same point, which is more convincing for the comparison of similarities. With reference to the design requirements of the pod cabin, the optimal solution of the pod structure is studied. The proximity of the standard index will be used as an important reference. In order to meet the requirements of high-quality structural design, No. 2 cabin should be selected in the sixth cabin.Figure 5
Test performance curve of test cabins No. 2 and No. 6 after zeroing at the starting point.
(a)(b)


### 4.3. Grey Relative Correlation Method

The gray relative association method is to characterize the approximate relationship between the behavior sequenceX0 and the rate of change of the comparison sequence Xi relative to the starting point. The closer the rate of change of X0 and Xi is, the larger the relative correlation degree will be.Definition 4.
LetX0 and Xi have the same length, and the initial values are not equal to 0. X0′ and Xi′ are their initial values, respectively, and then, the absolute correlation between X0′ and Xi′ is the relative degree of X0 and Xi.
According to the parameters of the existing data sequencesX0 and X1,X2,…,X10, the initial value processing is performed, that is, the data is divided by the first item number of the sequence, and the new series is obtained, which is consistent with the absolute gray correlation degree algorithm, and is written according to C++. Program operation get its corresponding gray relative relevance:(39)γ06=0.930,γ07=0.752,γ08=0.874,γ09=0.905,γ010=0.698,γ01=0.790,γ02=0.961,γ03=0.844,γ04=0.853,γ05=0.713.
It can be seen from the results of gray relative correlation that the change rate of the curve corresponding toX2 is the closest to the relative starting point of X0. The greater the relative correlation degree, the higher the approximation degree of the curve shape, which is reflected in the correlation degree value, that is, the relative test cabin No. 2 X0 has a high relative correlation. According to the gray relative association analysis theory, the gray relative correlation order is obtained, which is discharged from left to right according to the degree of relevance:(40)ζ=γ02,γ06,γ09,γ08,γ04,γ03,γ01,γ07,γ05,γ010.
According to the parameters of the existing data sequencesX0 and X2, X6, the starting line is zeroed, and the graph is placed on the coordinate axis x=0, so that the correlation between the behavior data sequence line X0 and the comparison data sequence line Xi passes through two. The sequence is characterized by a rate of change from the starting point.
Specifically, the reliability of the test prototype was evaluated. As shown in Figures5(a) and 5(b), relative to the starting point in (a), the rate of change of the data is better than the index in (b), which is essentially in the reflection curve. For the proximity between the two, refer to the design requirements of the pod and study an economical pod structure. It will be an important reference to the standard index of proximity. From the perspective of economy, choose No. 2 pod in cabin 6.


### 4.4. Grey Comprehensive Association Method

The gray comprehensive correlation method combines the curve approximation considered by absolute correlation analysis with the similarity considered by relative correlation analysis, adds weight coefficient, and adjusts according to the different emphasis of the technical parameters of airborne pod body to get the optimal data solution.Definition 5.
LetX0 and Xi have the same length, and the initial values are not equal to 0. ε01 and r01 are the absolute correlation and relative degree of X0 and Xi, respectively, and θ∈0,1 is called(41)ρ0i=θε0i+1−θr0i,that is, the comprehensive degree of association between X0 and Xi.
According to the abovementioned algorithm principle and the foregoing algorithm, the relative correlation operator is added on the basis of the absolute correlation degree, and the derived value is brought into the equation (41).
Due to technical requirements, the proximity and similarity of the reliability analysis curve of the pod cabin are important factors (taking economic factors as an important criterion, and the indicators far beyond the design standards will also be considered unreasonable); according to this reason, the weight coefficient is taken, after the operation, and the following data are obtained:(42)ρ01=0.741,ρ02=0.947,ρ03=0.824,ρ04=0.827,ρ05=0.651,ρ06=0.906,ρ07=0.695,ρ08=0.850,ρ09=0.890,ρ010=0.645.
According to the comprehensive correlation analysis, it can be clearly seen that the test cabin No. 2 can meet the reliability evaluation index better than the test cabin of the sixth test cabin.
However, Deng’s correlation analysis is the same as the abovementioned three correlation degree models. It is based on the proximity of the curve to infer the similarity between the curves. They still have a deficiency in distinguishing the two situations. In order to further distinguish the difference, a kind is given. The grey correlation degree is based on similarity and gray correlation degree model based on proximity [10].


### 4.5. Similarity Gray Correlation Degree and the Proximity Gray Correlation Degree

Definition 6.
LetXi and Xj be of the same length, si−sj, as shown in equation (29); say(43)εij=11+si−sj.
It is the gray correlation degree betweenXi and Xj based on similarity.
The similarity correlation is used to assess the degree of similarity between the sequencesXi and Xj in geometry. The more similar the Xi and Xj are in geometry, the larger the ε09 value, and vice versa.Definition 7.
LetXi and Xj be the same length, si−sj, as shown in equation (27):(44)si−sj=∫1nXi−Xjdt.
Then,(45)ρij=11+si−sj.
This is the proximity-based gray correlation betweenXi and Xj.
Proximity correlation is used to assess the proximity of sequenceXi and Xj in space. The closer the Xi and Xj are in space, the larger the ρij value, and vice versa. The similarity correlation only applies to sequences Xi and Xj. The situation is exactly the same, which is the reason for the quality grading of the 15 performance indicators of the cabin frame.
Based on the abovementioned theory, two correlation analysis are performed on 10 sets of data sequences, and the results are shown in Table2.
According to Table2, it can be seen that when the quality classification of the product of the pod cabin frame is carried out, the similarity gray correlation degree is used as the evaluation principle, that is, the overall geometric shape of the curve is similar and the “parallelism” is maintained within a certain range. Therefore, it is No. 2, No. 9, No. 6, No. 4, No. 8, No. 1, No. 3, No. 7, No. 10, and No. 5, when a certain characteristic of the product is analyzed as described in this section. For the proximity to be focused, the choice is 2, 5, 9, 7, 7, 3, 6, 6, 8, 1, and 10.
It should be noted that the value of proximity is small, but it does not mean that the reliability of the product is low. For example, the corresponding line of No. 6 test prototypes is far higher than the standard line and far away from the standard thread, so the proximity is low, which fully meets the design requirements, but the economy is poor. In addition, all performance test indicators below the technical requirements in the project are unqualified, according to the trend of typical line diagrams shown in Figure6:
According to Figures3(b) and 6(a)–6(d), test capsules No. 4, No. 6, No. 7, No. 9, and No. 10 perform well in terms of proximity and similarity. However, as a design standard prototype, it is obviously unqualified. It is found that some of its evaluation points are lower than the standard index points (characteristic behavior sequence data values), so the scope of consideration is excluded. Referring to the curve of the test chamber No. 2 shown in Figure 3(a), the test cabin No. 2 of the comprehensive correlation degree ρ02=0.947 is selected as the optima l cabin to be assembled with the pod bracket and its components.Table 2
Cabin data curve correlation degree property evaluation ranking.
RankingSymbolSimilarityRankingSymbolProximity1ε020.51281ρ020.17862ε090.40002ρ050.17863ε060.35093ρ090.16954ε040.29414ρ070.16675ε080.23535ρ030.16406ε010.21286ρ040.15397ε030.20417ρ060.15388ε070.18708ρ080.15389ε0100.16959ρ010.149310ε050.168110ρ0100.1429Figure 6
Performance test curve of test cabins No. 4, No. 6, No. 8, and No. 9. (a) Data line 4. (b) Data line 6. (c) Data line 8. (d) Data line diagram 9.
(a)(b)(c)(d)


## 4.1. Generalized Grey Relational Method

The generalized gray relational analysis model has three objectives for the analysis of pod reliability: first, the curve approach of Deng’s gray correlation analysis further demonstrates the basis of reliability optimal solution selection from a new perspective; second, the data sequence curve. The relative change rate problem between the two is added to the gray correlation analysis, and the original “point correlation analysis” is changed to the “face correlation analysis” between the curves. Starting from the overall situation of the curve distribution, the two curves are obtained by the starting point of zeroing. They are all placed on the coordinate axis, and the size of the area sandwiched between them is used as the criterion for judging the optimal solution of reliability analysis from the macroscopic perspective of the data. Third, the difference between the proximity and similarity of the curve is analyzed accordingly. Also, in the engineering practice, the corresponding factors of economic and design requirements are added. The generalized gray relational analysis is divided into three categories: absolute association method, relative association method, and comprehensive association method, which do not strictly satisfy the gray association four axioms.


## 4.2. Grey Absolute Correlation Method

Theorem 2.
Set the behavior sequenceXi=xi1,xi2,…,xin and record the corresponding polyline Xi−xi1 as Xi−xi1=xi1−xi1,xi2−xi1,…,xin−xi1; then,(23)si=∫1nXi−xi1dt.Definition 2.
Let the behavior sequenceXi=xi1,xi2,…,xin and D be defined as sequence operators and have the following relationship:(24)XiD=xi1d,xi2d,…,xind,(25)xikd=xik=xi1,k=1,2,…,n.
According to equations (24) and (25), D is called the starting point zeroing operator, and the starting point of XiD is zeroed, which is recorded as(26)Xi0=xi01,xi02,…,xi0n.Theorem 3.
We arbitrarily take two sequences of behavior:(27)Xi=xi1,xi2,…,xin,Xj=xj1,xj2,…,xjn.The starting point zero image is(28)Xi0=xi01,xi02,…,xi0n,Xj0=xj01,xj02,…,xj0n.It is referred to as(29)si−sj=∫1nXi0−Xj0dt.IfXi0 is always above Xj0, then si−sj≥0; if Xi0 is always below Xj0, then si−sj≤0; and if Xi0 intersects with Xj0, the si−sj sign is indefinite.Definition 3.
The sequence ofXi0: the sum of the time intervals between the observation data is the length of Xj0. It should be pointed out that the number of observations in two sequences of equal length is not necessarily equal.Theorem 4.
Let sequenceX0 and Xi be of the same length, s0 and si are as described in equation (15) of Theorem 2, and si−s0 is as described in equation (19) of Theorem 2; then,(30)ε0i=1+s0+si1+s0+si+s0−si.Equation (30) is the absolute degree of gray correlation between X0 and Xi.According to the satisfaction condition, for 11 sets of data sequences, since the time intervals between adjacent observation data are the same,X0 and X1,X2,…,X10 are said to constitute an isochronous sequence.Theorem 5.
Let sequenceX0 and Xi be the same length, and both are 1-time-distance sequences, with the corresponding relationship:(31)X00=x001,x002,…,x00n,Xi0=xi01,xi02,…,xi0n.The abovementioned sequence is the starting point zeroing image ofX0 and Xi, respectively, and the following correspondence is given:(32)s0=∑k=2n−1x00k+12x00n,(33)si=∑k=2n−1xi0k+12xi0n,(34)s0−si=∑k=2n−1xi0k−x00k+12xi0n−x00n.From equations (32)–(34), the following absolute correlation algorithm is derived:(35)ε0i=1+∑k=2n−1x00k+1/2x00n+∑k=2n−1xi0k+1/2xi0n1+∑k=2n−1x00k+1/2x00n+∑k=2n−1xi0k+1/2xi0n+∑k=2n−1xi0k−x00k+1/2xi0n−x00n.According to the gray absolute relevance theory, the related algorithm is written based on the C++ programming language, and the starting point zeroing process is performed according to the data sequence obtained byX0 and X1,X2,…,X10:(36)X00=0.00,2.00,−1.00,1.00,0.00,0.00,1.00,−1.00,1.00,1.00,2.00,0.00,−1.00,0.00,1.00,X10=0.00,1.40,−1.20,0.60,−0.30,0.00,0.90,−1.50,0.90,0.70,1.60,0.20,−1.60,−0.30,0.80,X20=0.00,2.10,−1.20,1.40,0.30,0.00,1.10,−0.80,1.00,0.80,2.20,0.20,−1.10,−0.20,1.30,X30=0.00,2.20,−0.40,1.30,0.20,0.00,1.10,−0.30,1.00,1.70,2.20,0.20,−0.40,−0.20,1.60,X40=0.00,1.60,−1.20,1.10,−0.40,0.10,1.00,−1.60,0.60,1.10,1.50,0.30,−1.10,−0.20,0.60,X50=0.00,1.30,−1.60,0.80,−0.10,−0.50,0.60,−1.40,0.80,0.30,1.70,−0.40,−1.70,−0.50,0.30,X60=0.00,2.20,−0.70,0.80,0.00,0.30,0.80,−0.60,1.60,1.10,2.00,−0.10,−0.80,0.10,1.30,X70=0.00,1.50,−1.40,0.80,0.00,−0.40,0.80,−1.40,0.50,0.90,1.70,−0.20,−1.60,−0.40,0.70,X80=0.00,2.30,−0.60,1.00,0.20,0.60,1.00,−0.60,1.50,1.30,2.00,0.10,−0.70,−0.10,1.50,X90=0.00,1.80,−0.70,0.80,0.10,−0.40,1.10,−1.40,0.80,1.20,1.60,−0.40,−0.70,−0.20,0.80,X100=0.00,1.40,−1.70,0.70,−0.50,0.00,0.60,−0.90,0.90,0.30,1.50,−0.30,−1.70,−0.10,0.80.The corresponding absolute relevance is(37)ε01=0.691,ε02=0.932,ε03=0.803,ε04=0.80,ε05=0.588,ε06=0.882,ε07=0.638,ε08=0.824,ε09=0.875,ε010=0.592.According to the theory of gray absolute correlation analysis, the corresponding C++ program is written to calculate the gray absolute correlation order, which is discharged from left to right according to the degree of relevance:(38)τ=ε02,ε06,ε09,ε08,ε03,ε04,ε01,ε07,ε010,ε05.It can be seen from the results of the absolute correlation degree that the results based on the “face correlation method” are numerically different from those based on the “point correlation method,” mainly because the considerations are different because the corresponding curve ofX1,X2,…,X10 and X0 to the absolute degree of correlation reflects the similarity between the curves.Specifically, to evaluate the reliability of the test prototype, as shown in Figures5(a) and 5(b), the area between the data line X2 and the data line X0 in (a) is smaller to the data line X6 in Figure 5(b). The area between the data lines X0 is essentially reflecting the similarity between the curves. After the starting point is zeroed, the starting points are classified into the same point, which is more convincing for the comparison of similarities. With reference to the design requirements of the pod cabin, the optimal solution of the pod structure is studied. The proximity of the standard index will be used as an important reference. In order to meet the requirements of high-quality structural design, No. 2 cabin should be selected in the sixth cabin.Figure 5
Test performance curve of test cabins No. 2 and No. 6 after zeroing at the starting point.
(a)(b)


## 4.3. Grey Relative Correlation Method

The gray relative association method is to characterize the approximate relationship between the behavior sequenceX0 and the rate of change of the comparison sequence Xi relative to the starting point. The closer the rate of change of X0 and Xi is, the larger the relative correlation degree will be.Definition 4.
LetX0 and Xi have the same length, and the initial values are not equal to 0. X0′ and Xi′ are their initial values, respectively, and then, the absolute correlation between X0′ and Xi′ is the relative degree of X0 and Xi.
According to the parameters of the existing data sequencesX0 and X1,X2,…,X10, the initial value processing is performed, that is, the data is divided by the first item number of the sequence, and the new series is obtained, which is consistent with the absolute gray correlation degree algorithm, and is written according to C++. Program operation get its corresponding gray relative relevance:(39)γ06=0.930,γ07=0.752,γ08=0.874,γ09=0.905,γ010=0.698,γ01=0.790,γ02=0.961,γ03=0.844,γ04=0.853,γ05=0.713.
It can be seen from the results of gray relative correlation that the change rate of the curve corresponding toX2 is the closest to the relative starting point of X0. The greater the relative correlation degree, the higher the approximation degree of the curve shape, which is reflected in the correlation degree value, that is, the relative test cabin No. 2 X0 has a high relative correlation. According to the gray relative association analysis theory, the gray relative correlation order is obtained, which is discharged from left to right according to the degree of relevance:(40)ζ=γ02,γ06,γ09,γ08,γ04,γ03,γ01,γ07,γ05,γ010.
According to the parameters of the existing data sequencesX0 and X2, X6, the starting line is zeroed, and the graph is placed on the coordinate axis x=0, so that the correlation between the behavior data sequence line X0 and the comparison data sequence line Xi passes through two. The sequence is characterized by a rate of change from the starting point.
Specifically, the reliability of the test prototype was evaluated. As shown in Figures5(a) and 5(b), relative to the starting point in (a), the rate of change of the data is better than the index in (b), which is essentially in the reflection curve. For the proximity between the two, refer to the design requirements of the pod and study an economical pod structure. It will be an important reference to the standard index of proximity. From the perspective of economy, choose No. 2 pod in cabin 6.


## 4.4. Grey Comprehensive Association Method

The gray comprehensive correlation method combines the curve approximation considered by absolute correlation analysis with the similarity considered by relative correlation analysis, adds weight coefficient, and adjusts according to the different emphasis of the technical parameters of airborne pod body to get the optimal data solution.Definition 5.
LetX0 and Xi have the same length, and the initial values are not equal to 0. ε01 and r01 are the absolute correlation and relative degree of X0 and Xi, respectively, and θ∈0,1 is called(41)ρ0i=θε0i+1−θr0i,that is, the comprehensive degree of association between X0 and Xi.
According to the abovementioned algorithm principle and the foregoing algorithm, the relative correlation operator is added on the basis of the absolute correlation degree, and the derived value is brought into the equation (41).
Due to technical requirements, the proximity and similarity of the reliability analysis curve of the pod cabin are important factors (taking economic factors as an important criterion, and the indicators far beyond the design standards will also be considered unreasonable); according to this reason, the weight coefficient is taken, after the operation, and the following data are obtained:(42)ρ01=0.741,ρ02=0.947,ρ03=0.824,ρ04=0.827,ρ05=0.651,ρ06=0.906,ρ07=0.695,ρ08=0.850,ρ09=0.890,ρ010=0.645.
According to the comprehensive correlation analysis, it can be clearly seen that the test cabin No. 2 can meet the reliability evaluation index better than the test cabin of the sixth test cabin.
However, Deng’s correlation analysis is the same as the abovementioned three correlation degree models. It is based on the proximity of the curve to infer the similarity between the curves. They still have a deficiency in distinguishing the two situations. In order to further distinguish the difference, a kind is given. The grey correlation degree is based on similarity and gray correlation degree model based on proximity [10].


## 4.5. Similarity Gray Correlation Degree and the Proximity Gray Correlation Degree

Definition 6.
LetXi and Xj be of the same length, si−sj, as shown in equation (29); say(43)εij=11+si−sj.
It is the gray correlation degree betweenXi and Xj based on similarity.
The similarity correlation is used to assess the degree of similarity between the sequencesXi and Xj in geometry. The more similar the Xi and Xj are in geometry, the larger the ε09 value, and vice versa.Definition 7.
LetXi and Xj be the same length, si−sj, as shown in equation (27):(44)si−sj=∫1nXi−Xjdt.
Then,(45)ρij=11+si−sj.
This is the proximity-based gray correlation betweenXi and Xj.
Proximity correlation is used to assess the proximity of sequenceXi and Xj in space. The closer the Xi and Xj are in space, the larger the ρij value, and vice versa. The similarity correlation only applies to sequences Xi and Xj. The situation is exactly the same, which is the reason for the quality grading of the 15 performance indicators of the cabin frame.
Based on the abovementioned theory, two correlation analysis are performed on 10 sets of data sequences, and the results are shown in Table2.
According to Table2, it can be seen that when the quality classification of the product of the pod cabin frame is carried out, the similarity gray correlation degree is used as the evaluation principle, that is, the overall geometric shape of the curve is similar and the “parallelism” is maintained within a certain range. Therefore, it is No. 2, No. 9, No. 6, No. 4, No. 8, No. 1, No. 3, No. 7, No. 10, and No. 5, when a certain characteristic of the product is analyzed as described in this section. For the proximity to be focused, the choice is 2, 5, 9, 7, 7, 3, 6, 6, 8, 1, and 10.
It should be noted that the value of proximity is small, but it does not mean that the reliability of the product is low. For example, the corresponding line of No. 6 test prototypes is far higher than the standard line and far away from the standard thread, so the proximity is low, which fully meets the design requirements, but the economy is poor. In addition, all performance test indicators below the technical requirements in the project are unqualified, according to the trend of typical line diagrams shown in Figure6:
According to Figures3(b) and 6(a)–6(d), test capsules No. 4, No. 6, No. 7, No. 9, and No. 10 perform well in terms of proximity and similarity. However, as a design standard prototype, it is obviously unqualified. It is found that some of its evaluation points are lower than the standard index points (characteristic behavior sequence data values), so the scope of consideration is excluded. Referring to the curve of the test chamber No. 2 shown in Figure 3(a), the test cabin No. 2 of the comprehensive correlation degree ρ02=0.947 is selected as the optima l cabin to be assembled with the pod bracket and its components.Table 2
Cabin data curve correlation degree property evaluation ranking.
RankingSymbolSimilarityRankingSymbolProximity1ε020.51281ρ020.17862ε090.40002ρ050.17863ε060.35093ρ090.16954ε040.29414ρ070.16675ε080.23535ρ030.16406ε010.21286ρ040.15397ε030.20417ρ060.15388ε070.18708ρ080.15389ε0100.16959ρ010.149310ε050.168110ρ0100.1429Figure 6
Performance test curve of test cabins No. 4, No. 6, No. 8, and No. 9. (a) Data line 4. (b) Data line 6. (c) Data line 8. (d) Data line diagram 9.
(a)(b)(c)(d)


## 5. Conclusions

(1)
According to the gray system theory, designers can quickly and effectively optimize the less relevant cabins and propose improvements. It provides a numerical basis for further research on pod structure optimization and vibration safety inspection.(2)
Both Dun’s gray correlation method and the generalized gray correlation method select the test cabin No. 2 as the optimal solution. However, it can be found that the possibility that the similar test chamber No. 6 is selected as the optimal solution still exists, for the “point correlation.” The conclusions are often more convincing with the analysis of “face associations,” and the considerations are more comprehensive and thorough.(3)
In the generalized gray relational analysis model, the concept of relative change rate between data curves is introduced. Starting from the overall situation of curve distribution, the size of the area between them is used as the criterion for judging the selection of the optimal test prototype for such aviation products, the basis.(4)
In the virtual design of the airborne pod, the reliability analysis of the pod cabin based on the generalized correlation method is introduced. The data curve similarity and proximity are practically applied in the reliability analysis of the pod structure, and the details are as follows. For the difference and connection of the reliability evaluation of mechanical products, the optimal solution for the economic and design requirements is obtained as the No. 2 cabin.(5)
Through the analysis of the generalized correlation method for each evaluation index by the relative gray correlation model, the reliability analysis of the airborne pod cabin frame is finally realized, which provides a specific theoretical basis for determining the optimal structure of the airborne pod cabin frame.



---

*Source: 8763642-2020-08-17.xml*