Articles | Volume 13, issue 1
Mech. Sci., 13, 535–542, 2022
https://doi.org/10.5194/ms-13-535-2022

Special issue: Advances in Service and Industrial Robotics – RAAD2021

Mech. Sci., 13, 535–542, 2022
https://doi.org/10.5194/ms-13-535-2022
Research article
15 Jun 2022
Research article | 15 Jun 2022

The singular value decomposition method of improved incidence matrix for isomorphism identification of epicyclic gear trains

The singular value decomposition method of improved incidence matrix for isomorphism identification of epicyclic gear trains
Mingshuai Zhou1,2, Wei Sun1,2, Rongxuan Wu3, and Wenxu Fu4 Mingshuai Zhou et al.
  • 1Key Laboratory of Metallurgical Equipment and Control Technology, Ministry of Education, Wuhan University of Science and Technology, Wuhan 430081, China
  • 2Hubei Key Laboratory of Mechanical Transmission and Manufacturing Engineering, Wuhan University of Science and Technology, Wuhan 430081, China
  • 3College of Machinery and Automation, Wuhan University of Science and Technology, Wuhan 430081, China
  • 4International College, Wuhan University of Science and Technology, Wuhan 430081, China

Correspondence: Wei Sun (sw@wust.edu.cn)

Abstract

The epicyclic gear train (EGT) is an advanced gear transmission mechanism, which is widely applied in drive systems. It is of great significance to eliminate the same structure in the type synthesis of EGTs. In this paper, an isomorphism identification method of EGTs based on the singular value decomposition of improved incidence matrix is proposed. Firstly, the improved incidence matrix is used to describe the structure of EGTs. Then, the degree sequence of links and kinematic pairs can be extracted and used as a basic information for preliminary screening. The improved incidence matrix can effectively distinguish multiple joints. Next, the topological connection relationship between the links and the kinematic pairs is extracted by the singular value decomposition of the matrix, which is used as the final value for isomorphism identification of EGTs. Finally, the effectiveness of this method is verified through a lot of examples.

1 Introduction

The epicyclic gear train (EGT) has been commonly used in the transmission of various high-end mechanical equipment due to its advantages of large transmission ratio range, compact structure, large bearing capacity, etc. The configuration synthesis is an effective method to ensure the discovery of potential high-performance EGT structures. The purpose of configuration synthesis is to generate all feasible configurations required by a given number of links and degrees of freedom. In the process of configuration synthesis, it is inevitable to produce numerous EGTs with the same structure, so the isomorphism identification of EGTs shows its important significance. However, due to the large number of configurations generated during the synthesis process, manual screening alone is not practical. Therefore, scholars are committed to developing an algorithm that can automatically distinguish repetitive structures. Since the 1950s and 1960s, the isomorphism identification of kinematic chains has always been a theoretical conundrum in the process of their configuration synthesis, and there is still no completed solution to solve this problem (the EGT is a special kinematic chain with gear pairs). To solve the problem of isomorphism identification of kinematic chains, the proposed methods are mainly divided into the following two categories:

  1. The first one is to find the characteristic constant of kinematic chains. Yan and Hall (1982, 1981), Tsai (1987), and Tsai and Lin (1989) utilized methods of the characteristic polynomial based on adjacency matrix for isomorphism identification of kinematic chains, but there are cases of failure. The standard code method (Shin and Krishnamurty, 1992, 1993, 1994) was proposed by Shin et al., but its detection efficiency is not high. The Hamming number (Rao and Rao, 1993a, b; Rao and Raju, 1991) was also used for isomorphism identification of kinematic chains. Its accuracy of the identification is superior, which can effectively differentiate the structural information of the matrix, but there are still some failure cases. Sun et al. (2017) used an improved Hamming number method to solve the problem of isomorphism identification of kinematic chains with multiple joints and have achieved good results, but their method does not have the ability to decode. The eigenvalues and eigenvectors of the adjacency matrix (Cubillo and Wan, 2005; Zhou et al., 2018; Liu and Yu, 2012) were utilized for isomorphism identification, with a certain accuracy, but there are also counterexamples. The loop method (Yang and Ding, 2018; Prasad Raju Pathapati and Rao, 2002) was applied for isomorphism identification of kinematic chains, but the number of loops will increase sharply with the increase of the number of vertices in the topological graph, which will cause a large number of operations. The distance matrix (Yadav et al., 1995; Yadav et al., 1996) was used to describe the connection relationship between the links of the kinematic chain, which represents the shortest distance between the links. However, the structural information contained in the matrix will be destroyed due to some highly symmetrical structures. In recent years, Sun et al. (2020) proposed a set of methods based on the power of the adjacency matrix to extract the feature invariants of kinematic chains, which can be used for the kinematic chain or the EGT, and have achieved good results. Deng et al. (2020) proposed a method for isomorphism identification using topological index, which can also be applied to the EGT. It has achieved good results.

  2. The second one is to detect isomorphism from the isomorphic definition of graphs. It is common to start from the topological graph of the kinematic chain and quickly search for the bijective relationship between the vertices and edges of the isomorphic graph through optimization algorithms. Genetic algorithms (Rao, 2003; Feng and Chen, 2001) and an ant algorithm (Yang et al., 2007) were used to identify isomorphism. The advantage of these methods is the fact that they will not easily fall into the local optimal solutions but easily find the global optimal solutions. However, their convergence speed is slow, and they cannot guarantee a high enough correct rate. The neural network algorithm (Rao, 2002) also has similar characteristics.

The research method of pure-rotating-pair kinematic chains is relatively mature, but the problem of EGTs has become more complicated due to the addition of gear pairs in ordinary kinematic chains, which leads to the inapplicability of many existing isomorphism identification methods. Therefore, people have to explore new solutions to expand the scope of application of algorithms.

For the isomorphic identification of the EGT, it is most common to use an adjacent matrix to describe the EGT, but few scholars use the incidence matrix to describe it, mainly because the incidence matrix is not necessarily a square matrix, and the subsequent data processing is difficult. Starting from improving the incidence matrix, the singular value decomposition of the matrix is applied to extract the incidence matrix invariants, which is used as a judgment for the isomorphic identification of the EGT. The method has certain advantages in accuracy and efficiency.

2 Topological descriptions of the EGT

In this paper, the structure of the EGT is described in a topological graph. The topological graph can clearly show its structural characteristics, and it is convenient to convert it into a matrix or array to quantify the structural information of the EGT. Figure 1a shows a three-dimensional model of a Simpson planetary gear train, and Fig. 1b shows its structural graph. This paper utilizes the bicolor topological graph by Yang and Ding (2018) to represent the EGT. As shown in Fig. 1c, the solid vertex represents the link, the hollow vertex represents the multiple joints, the solid line edge represents the rotating pair, and the dotted line edge represents the gear pair. The edges associated with the hollow vertex can be regarded as the same kinematic pair. The serial numbers of links and kinematic pairs are marked as “ i ” and “ i ”, respectively, in Fig. 1d.

https://ms.copernicus.org/articles/13/535/2022/ms-13-535-2022-f01

Figure 1Topological description of the EGT.

Download

3 The improved incidence matrix

In order to facilitate analysis and processing, the structural information of the topology can be stored in a digital matrix. Consequently, this paper proposes an improved incidence matrix to describe the topological graph of the EGT.

Definition 1. The definition of the improved incidence matrix A is as follows:

(1) A = a i j m × n ,

where i (1, 2, …, m) represents the ith vertex (i.e., the ith link), and j (1, 2, …, n) represents the jth edge (i.e., the jth kinematic pair); then,

(2) a i j = 1 Link i is  adjacent  to  Rotating   Pair j , 2 Link i is  adjacent  to  Gear  Pair j , 0 Link i is  not  adjacent  to  Pair j .

For example, the following equation (Eq. 3) is the improved incidence matrix corresponding to the EGT topological graph in Fig. 1, in which the second row represents Link 2, which is adjacent to Gear Pair 1 in the first column and Rotating Pair 7 in the seventh column, respectively.

(3) A 6 × 7 = 0 0 0 0 1 2 1 2 0 0 0 0 0 1 2 2 0 0 1 0 0 0 2 2 0 0 0 1 0 0 2 1 0 2 0 0 0 0 1 0 0 1

The improved incidence matrix Am×n (i.e., the link–kinematic-pair matrix) can clearly express the information about links and kinematic pairs of the EGT. Therefore, the basic structural constants obtained from the matrix are as follows:

  1. The number of rows of the matrix represents the number of links of the EGT, and the number of columns of the matrix represents the number of kinematic pairs of the EGT.

  2. Nonzero elements in Row i of the matrix represent kinematic pairs adjacent to Link i, and different numbers represent different types of kinematic pairs.

  3. The number of nonzero elements in Column j of the matrix is the number of links adjacent to Kinematic Pair j. When it is more than two, it is with multiple joints, i.e., the compound kinematic pair.

Definition 2. The code of link (CL) is the kinematic pair information connected with the corresponding link; the sequence of the CL from large to small is called the CLS, which is expressed as CL in the paper. It is defined as follows:

(4) C L = [ c i l ] 1 × m ,

where

(5)cil=rpi+gpi10-p;(6)10p-1max{gpi}<10p,pN+.

Here, rpi represents the number of rotating pairs associated with Link i, and gpi represents the number of gear pairs associated with Link i.

Definition 3. The code of kinematic pair (CKP) is the number of links connected with the corresponding kinematic pair; the sequence of the CKP from large to small is called the CKPS, which is expressed as Ckp in the paper. It is defined as follows:

(7) C kp = c j kp 1 × n ,

where

(8) c j kp = n l j .

Here, nlj indicates the number of links adjacent to Kinematic Pair j.

https://ms.copernicus.org/articles/13/535/2022/ms-13-535-2022-f02

Figure 2The flowchart of isomorphism identification.

Download

https://ms.copernicus.org/articles/13/535/2022/ms-13-535-2022-f03

Figure 3The topological graphs of 6-link EGTs a1 and a2.

Download

4 Isomorphism identification

Definition 4. (Feng and Chen, 2001) Given undirected graphs G1=(V1,E1) and G2=(V2,E2) and if there is bijection f such that u,vV1, [u,v]E1[f(u),f(v)]E2, then G1 and G2 are isomorphic and denoted as G1G2.

According to the definition of isomorphism, the nodes of the topological graph with the same structure have a one-to-one mapping relationship. It is deduced that the topological adjacency relationship between nodes of isomorphic topological graph is also the same. There is then the following definition:

Definition 5. If the number and type of links and kinematic pairs and the topological constraints between links in two EGTs are identical, then these two EGTs could be considered to be isomorphic.

For the above basic structural constants, the CLS and the CKPS can be obtained directly from the improved incidence matrix. However, it is difficult to extract directly the topological constraints between links. Generally, the eigenvalue and eigenvector method is used to extract the characteristic structural constants of matrices, but the precondition for this method is that the matrix must be a square matrix. The numbers of rows and columns of the improved incidence matrix are related to the numbers of links and kinematic pairs, so it is not necessarily a square matrix. The singular value does not change with the exchange of rows and columns of the matrix, so it is a characteristic invariant of the matrix. Therefore, the singular value of the improved incidence matrix can be used as the structural characteristic information of the topological connection between the EGT's links.

https://ms.copernicus.org/articles/13/535/2022/ms-13-535-2022-f04

Figure 4The topological graphs of 10-link EGTs b1 and b2.

Download

For the matrix ARm×n, where m and n are any positive integers, the following matrix decomposition theorem exists:

Theorem 1. If ARrm×n, where r>0, then there is the orthogonal matrix Um×m, Vn×n, which makes

(9) A = UDV T ,

where Dm×n is a rectangular diagonal matrix, and r is the number of nonzero elements on the diagonal of the matrix.

It has the following definition: if ARm×n, then the nonnegative square root of the nonzero eigenvalue of the matrix ATA is the singular value of matrix A, and all of them are denoted as σ(A). Equation (9) is called the singular value decomposition (SVD). The calculation formula of SVD is as follows:

(10) A T A x = λ x ,

where λ is the eigenvalue of the matrix ATA, and x is the eigenvector corresponding to λ.

(11) | A T A - λ E | = 0

By solving Eq. (11), we can obtain the eigenvalue set Λ=[λ]1×r of the matrix ATA. Thus, the singular value set σ(A) of matrix A is

(12) σ ( A ) = λ 1 × r .

Therefore, the algorithmic steps of isomorphism identification are as follows:

  • Step 1. One should attain the bicolor topological graphs of the two EGTs to be compared and construct the improved incidence matrices of the two EGTs.

  • Step 2. One should compare the basic structure constants of the two matrices, including the sequence of the CL and the sequence of the CKP. If they are the same, then go to Step 3; if not, they are isomeric, so go to Step 4.

  • Step 3. Calculate the singular values of the two incidence matrices. If they are the same, then they are considered to be isomorphic; otherwise, they are isomeric, so go to Step 4.

  • Step 4. The algorithm is over.

The isomorphism identification flowchart is given in Fig. 2.

https://ms.copernicus.org/articles/13/535/2022/ms-13-535-2022-f05

Figure 5The topological graphs of 10-link EGTs c1 and c2.

Download

5 Numerical examples

Case 1. Figure 3 shows the bicolor topological graphs of two 6-link EGTs with multiple joints. The method in this paper is utilized in the isomorphism identification of the two EGTs, and the specific process is as follows.

The two improved incidence matrices of the two EGTs are as follows:

(13)Aa1=200010021000020120100002201000021020000110,(14)Aa2=200010021000200120100002100200011000000122,(15)Ca1L=Ca2L=[2.1,2,1.2,1.2,1.2,1.1],(16)Ca1kp=Ca2kp=[4,2,2,2,2,2,2].

That is, the basic structure constants of the two EGTs are the same.

The singular value decompositions of the two matrices are performed:

(17)σ(Aa1)=[4.0851,3.3098,2.5041,1.6807,1.5164,0.9811],(18)σ(Aa2)=[4.0851,3.3098,2.5041,1.6807,1.5164,0.9811].

Here, σ(Aa1)=σ(Aa2), which proves that the two EGTs are isomorphic. By changing the positions of vertex 6 and vertex 5 in Fig. 3, right, we can know that they are isomorphic. The result of isomorphism identification is consistent with reality.

Case 2. Figure 4 shows the topological graphs of two 8-link EGTs. It is observed that the two are isomeric. The method in this paper is used in the isomorphism identification of the two EGTs, and the specific process is as follows.

The two improved incidence matrices of the two EGTs are as follows:

(19)Ab1=20021000002200010000021000002200120110000000001220000001020000001100000000010002,(20)Ab2=20021000002200010010021000200200120100000000012000000011000000000102000000000212,(21)Cb1L=[3.1,2,1.3,1.2,1.2,1.2,1.1,1.1],(22)Cb2L=[2.2,2.1,2.0,1.3,1.2,1.2,1.1,1.1],(23)Cb1kp=[5,2,2,2,2,2,2,2,2,2],(24)Cb2kp=[5,2,2,2,2,2,2,2,2,2].

Here, Cb1kp=Cb2kp but Cb1LCb2L, which shows that the two EGTs are isomeric. The result of isomorphism identification is consistent with reality. The judgment process is over in advance.

Case 3. Figure 5 shows the topological graphs of two 10-link EGTs. It is carefully observed that the two EGTs are isomorphic. Then the method in this paper is employed in the isomorphism identification of the two EGTs.

(25)Ac1=2000100020000220000000100002101000000000012000000200000210000102200001200000000000021000002000010110000000000001200200000100000200,(26)Ac2=2000001200000220100000000002200000100000021000001000000120000000000002100000020001011000000000000000122000010000002000001000210022,(27)Cc1L=Cc2L=[3,2.3,2.1,1.2,1.2,1.2,1.2,1.2,1.1,1.1],(28)Cc1kp=Cc2kp=[6,2,2,2,2,2,2,2,2,2,2,2,2],(29)σ(Ac1)=σ(Ac2)=[4.5835,3.7688,3.3779,3.2899,2.7814,2.1596,1.7020,1.4714,1.3280,1.1526].

These two EGTs are isomorphic, which means that the result of isomorphism identification is consistent with reality.

The three above cases show two different results, namely, isomorphism and isomerism. Since the number of links of the EGTs in the first two cases is small, we can easily observe whether the two EGTs are isomorphic manually. However, in the last case, the EGTs have a large number of links and complex topology, so it will be difficult to manually judge whether they are isomorphic. Moreover, in the process of configuration synthesis, a large number of EGT structures will be generated, and manual comparison will be very time-consuming. Therefore, the automatic algorithm of isomorphism identification shows its main significance.

6 Conclusions
  1. In this paper, the improved incidence matrix is used to describe the EGT structure. The improved incidence matrix is convenient to retain the information of links and kinematic pairs of the EGT. Compared with the adjacency matrix, it can effectively and intuitively determine multiple joints.

  2. Compared with some mainstream methods, the SVD method is very simple in implementation and judgment efficiency is high. For example, the loop method needs to find all loops (or loops with given requirements) and then classify and code these. As the number of loops increases, the amount of computation increases exponentially. The adjacency matrix eigenvalue and eigenvector method utilizes the link-to-link adjacency matrix to store structural information. The method is relatively simple, but there are many counterexamples and low reliability. The standard code is used for isomorphism identification, but the power of the matrix needs to be used many times to compare the symmetrical links in the mechanism. It is computationally expensive and has many steps. Relatively speaking, the SVD method only needs to compare CLS and CKPS for preliminary discrimination, and the two can be obtained directly from the improved incidence matrix. Calculating the singular value of the improved incidence matrix, the topological connection relationship between EGTs' links can be obtained. At present, there have been no failure cases.

Data availability

All data included in this study are available upon request by contacting the corresponding author.

Author contributions

MZ conceived the original ideas and drafted the manuscript. WS supervised the project and reviewed the writing. RW produced the figures in this article. WF helped with the writing and language.

Competing interests

The contact author has declared that neither they nor their co-authors have any competing interests.

Disclaimer

Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Special issue statement

This article is part of the special issue “Advances in Service and Industrial Robotics – RAAD2021”. It is a result of the 30th International Conference on Robotics in Alpe-Adria-Danube Region, RAAD 2021, Futuroscope-Poitiers, France, 21–23 June 2021.

Acknowledgements

This work is supported in part by the National College Students Innovation and Entrepreneurship Training Program (China) (grant no. S202110488043), the National Natural Science Foundation of China (grant no. 51875418), Youth Fund of Hubei Science and Technology Department (grant no. 2021CFB135), Foundation of Hubei Provincial Education Department (grant no. B2020011), and Wuhan University of Science and Technology (WUST) National Defense Pre-research Foundation (grant no. GF202008).

Financial support

This research has been supported by the Youth Fund of Hubei Science and Technology Department (grant no. 2021CFB135), and the National College Students Innovation and Entrepreneurship Training Program (China) (grant no. S202110488043), National Natural Science Foundation of China (grant no. 51875418), Foundation of Hubei Provincial Education Department (grant no. B2020011), Wuhan University of Science and Technology (WUST) National Defense Pre-research Foundation (grant no. GF202008).

Review statement

This paper was edited by Mohamed Amine Laribi and reviewed by Joonhong Park and one anonymous referee.

References

Cubillo, J. P. and Wan, J.: Comments on mechanism kinematic chain isomorphism identification using adjacent matrices, Mech. Mach. Theory, 40, 131–139, https://doi.org/10.1016/j.mechmachtheory.2004.07.004, 2005. 

Deng, T., Xu, H., Tang, P., Liu, P., and Yan, L.: A novel algorithm for the isomorphism detection of various kinematic chains using topological index, Mech. Mach. Theory, 146, 103740, https://doi.org/10.1016/j.mechmachtheory.2019.103740, 2020. 

Feng, C. and Chen, Y.: Mechanism kinematics chain isomorphism identification based on genetic algorithms, Chin. J. Mech. Eng., 37, 27–30, https://doi.org/10.3901/JME.2001.10.027, 2001. 

Liu, J. and Yu, D.: Representations & isomorphism identification of planar kinematic chains with multiple joints based on the converted adjacent matrix, J. Mech. Eng., 48, 15–21, https://doi.org/10.3901/JME.2012.05.015, 2012. 

Prasad Raju Pathapati, V. V. N. R. and Rao, A. C.: A new technique based on loops to investigate displacement isomorphism in planetary gear trains, J. Mech. Design, 124, 662–675, https://doi.org/10.1115/1.1503373, 2002. 

Rao, A. C.: An artificial neural network approach to mechanism kinematic chain isomorphism identification, Mech. Mach. Theory, 37, 549–551, https://doi.org/10.1016/S0094-114X(02)00004-6, 2002. 

Rao, A. C.: A genetic algorithm for epicyclic gear trains, Mech. Mach. Theory, 38, 135–147, https://doi.org/10.1016/S0094-114X(02)00094-0, 2003. 

Rao, A. C. and Raju, D. V.: Application of the hamming number technique to detect isomorphism among kinematic chains and inversions, Mech. Mach. Theory, 26, 55–75, https://doi.org/10.1016/0094-114X(91)90022-V, 1991. 

Rao, A. C. and Rao, C. N.: Loop based pseudo hamming values–I testing isomorphism and rating kinematic chains, Mech. Mach. Theory, 28, 113–127, https://doi.org/10.1016/0094-114X(93)90051-V, 1993a. 

Rao, A. C. and Rao, C. N.: Loop based pseudo hamming values–II inversions, preferred frames and actuators, Mech. Mach. Theory, 28, 129–143, https://doi.org/10.1016/0094-114X(93)90052-W, 1993b. 

Shin, J. K. and Krishnamurty, S.: Development of a standard code for colored graphs and its application to kinematic chains, Proceedings of the ASME 1992 Design Technical Conferences, 22nd Biennial Mechanisms Conference: Flexible Mechanisms, Dynamics, and Analysis, Scottsdale, Arizona, USA, 13–16 September 1992, DETC1992-0387, https://doi.org/10.1115/DETC1992-0387, 1992. 

Shin, J. K. and Krishnamurty, S.: Standard code technique in the enumeration of epicyclic gear trains, Mech. Mach. Theory, 28, 347–355, https://doi.org/10.1016/0094-114X(93)90075-7, 1993. 

Shin, J. K. and Krishnamurty, S.: On identification and canonical numbering of pin-jointed kinematic chains, J. Mech. Design, 116, 182–188, https://doi.org/10.1115/1.2919344, 1994. 

Sun, L., Cui, R., Ye, Z., Zhou, Y., Xu, Y., and Wu, C.: Similarity recognition and isomorphism identification of planar kinematic chains, Mech. Mach. Theory, 145, 103678, https://doi.org/10.1016/j.mechmachtheory.2019.103678, 2020. 

Sun, W., Kong, J., and Sun, L.: The improved hamming number method to detect isomorphism for kinematic chain with multiple joints, J. Adv. Mech. Des. Syst., 11, JAMDSM0061, https://doi.org/10.1299/jamdsm.2017jamdsm0061, 2017.  

Tsai, L. W.: An application of the linkage characteristic polynomial to the topological synthesis of epicyclic gear trains, J. Mech. Trans. and Automation, 109, 329–336, https://doi.org/10.1115/1.3258798, 1987. 

Tsai, L. W. and Lin, C. C.: The creation of nonfractionated, two-degree-of-freedom epicyclic gear trains, J. Mech. Trans. and Automation, 111, 524–529, https://doi.org/10.1115/1.3259033, 1989. 

Yadav, J. N., Pratap, C. R., and Agrawal, V. P.: Computer aided detection of isomorphism among kinematic chains and mechanisms using the concept of modified distance, Mech. Mach. Theory, 31, 439–444, https://doi.org/10.1016/0094-114X(95)00078-D, 1996. 

Yadav, J. N., Pratap, C. R., and Agrawal, V. P.: Detection of isomorphism among kinematic chains using the distance concept, J. Mech. Design, 117, 607–611, https://doi.org/10.1115/1.2826728, 1995. 

Yan, H. S. and Hall, A. S.: Linkage characteristic polynomials: definitions, coefficients by inspection, J. Mech. Design, 103, 578–584, https://doi.org/10.1115/1.3254957, 1981. 

Yan, H. S. and Hall, A. S.: Linkage characteristic polynomials: assembly theorems, uniqueness, J. Mech. Design, 104, 11–20, https://doi.org/10.1115/1.3256301, 1982. 

Yang, P., Pei, Z., Liao, N., and Yang, B.: Isomorphism identification for epicyclic gear mechanism based on mapping property and ant algorithm, Eng. Comput., 23, 49–54, https://doi.org/10.1007/s00366-006-0041-y, 2007. 

Yang, W. and Ding, H.: The perimeter loop-based method for the automatic isomorphism detection in planetary gear trains, J. Mech. Design, 140, 123302, https://doi.org/10.1115/1.4041572, 2018. 

Zhou, Z., Kong, J., Sun, L., and Xiao, Q.: Representation and isomorphism identification of planar kinematic chains with multiple joints on value-added matrix, Mech. Sci. Tech. Aer. Eng., 37, 657–662, https://doi.org/10.13433/j.cnki.1003-8728.2018.0501, 2018. 

Download
Short summary
The epicyclic gear train (EGT) is a gear transmission mechanism. When the number of designed structures is increasing and structures are becoming ever more complex, manual screening is very difficult. Therefore, scholars are committed to developing an algorithm that can distinguish repetitive structures. In this paper, we use a method to transform the structural information of EGT into a matrix and extract the hidden information in the matrix to judge whether there is a repeated structure.