A Novel Hybrid Algorithm for the Forward Kinematics Problem of 6 DOF Based on Neural Networks

The closed kinematic structure of Gough–Stewart platforms causes the kinematic control problem, particularly forward kinematics. In the traditional hybrid algorithm (backpropagation neural network and Newton–Raphson), it is difficult for the neural network part to train different datasets, causing training errors. Moreover, the Newton–Raphson method is unable to operate on a singular Jacobian matrix. In this study, in order to solve the forward kinematics problem of Gough–Stewart platforms, a new hybrid algorithm is proposed based on the combination of an artificial bee colony (ABC)–optimized BP neural network (ABC–BPNN) and a numerical algorithm. ABC greatly improves the prediction ability of neural networks and can provide a superb initial value to numerical algorithms. In the design of numerical algorithms, a modification of Newton’s method (QMn-M) is introduced to solve the problem that the traditional algorithm model cannot be solved when it is trapped in singular matrix. Results show that the maximal improvement in ABC–BPNN error optimization was 46.3%, while the RMSE index decreased by 42.1%. Experiments showed the feasibility of QMn-M in solving singular matrix data, while the percentage improvement in performance for the average number of iterations and required time was 14.4% and 13.9%, respectively.


Introduction
Gough-Stewart platforms (GSPs) are typical six-degree-of-freedom (6-DOF) parallel robots that have the advantages of high rigidity, high precision, and large carrying capacity [1]. GSPs are widely used in various large-scale motion simulation platforms such as flight, automotive, ship, and tank simulators. Generally, GSPs comprise three structures: universal and cylindrical and universal joint (6-UCU), 6-UPS, and 6-SPS, in which 6 represents six identical structures, U represents a universal joint, C represents a cylindrical joint, P represents a prismatic joint, and S represents a spherical joint. GSPs can be divided into 6-6, 3-3, and 6-3 structures according to the connection modes of the upper and lower hinge points. In this paper, the 6-UCU structure and 6-6 connection model are used. GSPs, with its structure and benefits, compensates for the many deficiencies of serial robots. It has a wide variety of applications in the future industry and intelligent manufacturing.
Because it can give mapping between a Cartesian and joint space, the kinematics problem is critical for parallel robotics. Solving the forward kinematics problem (FKP) is an essential step in the modeling and control of parallel robots, particularly for real-time applications. The FKP is a challenging and essential robotics topic in GSPs. Due to the high nonlinearity and the varied closed-loop kinematic architectures of parallel robots [2], there is currently no acknowledged generic solution to solve the FKP. To control GSPs, advanced algorithms with higher computation loads that achieve better performance can be used. Therefore, it is significant to find a more powerful algorithm to reduce the computation time of the FKP. The FKP of GSPs is, therefore, a hot topic for researchers.
Studies on the FKP of GSPs can be classified into the traditional approach and the intelligent algorithm approach. Conventional traditional approaches include the numerical method, analytical method, and adding extra sensors. The numerical method can obtain an iterative solution without suffering the multisolution problem [3]. The Newton-Raphson (NR) method is a prominent numerical approach for analyzing a parallel robot's forward kinematics [4]. The NR method with good initial values is used to numerically solve nonlinear FKP equations for root finding. The method is overly sensitive to initial values, and results diverge if the initial values are inappropriate. The analytical method tends to find a closedform solution [5,6]. This method forms a series of sophisticated derivations and is only useful for certain structures. The method of adding more sensors [7] is a strategy to find a unique approach with the least amount of computing. Unfortunately, due to the high cost of this method, and measurement and assembly faults, its application range is limited.
With the advancement of computer application technology and artificial intelligence, an increasing number of researchers have begun to use intelligent algorithms to solve the FKP, such as artificial neural networks (ANNs) and support vector machines (SVM). Morell et al. [8] solved FKP by using the support vector regression approach, which had a unique notion, but the model training time was too long and would not completely satisfy the requirements in real-time control. ANNs are employed in the field of intelligent algorithms to train the inverse kinematics values of the GSP to produce a set of forward kinematics solutions [9]. Subsequently, the NR method is used to obtain the exact values of the approximations. Using the hybrid strategy to solve FKP has been recognized by some scholars. This hybrid strategy [10] can be combined with a high level of solving speed and accuracy. Moreover, it can leverage the superb initial value provided by a neural network to perform a Newtonian iteration, which compensates for the Newtonian iteration's sensitivity to the initial value. It also overcomes the problem of a neural network's predicted value having a forecasting deviation. The approach is adaptable to parallel robots with a variety of structural properties [11]. Many researchers have used this foundation to improve the algorithm's efficiency and accuracy in real-time motion. To ensure the feasibility of realtime application, Zhu et al. exploited the optimization of numerical iteration efficiency by creating a deviation-driven algorithm [12]. After that, a velocity Jacobian matrix iteration was used to improve the hybrid algorithm's overall speed [13]. The advantages of the neural network and genetic algorithm are combined in the neurogenetic algorithm [14], which maximizes the robot's workspace volume.
For the neural network part, the well-known BP neural network (BPNN) is mostly employed in the traditional hybrid strategy. BPNN has good nonlinear mapping ability, i.e., f {L O1 , L O2 , L O3 , L O4 , L O5 , L O6 } → {α, β, γ, x, y, z}, which is the relationship between joint displacement and end-effector pose. Although the BPNN model can achieve promising performance and has a flexible network structure, it still has some drawbacks: (i) For BPNNs, the training effect depends on the dataset, i.e., limited self-adaptation ability for different datasets. (ii) It is possible to fall into a local optimum for the model system with six inputs and six outputs, which cannot ensure a global training effect. (iii) The model is sensitive to initial parameters. Randomly initializing weights and thresholds impacts the model. Moreover, there is also a flaw with the numerical algorithm part. (iv) A set of inverse of Jacobian matrices J −1 are formed in the process of solving nonlinear equations with the Newton-Raphson method [12]. That is, when the Jacobian matrix is singular, the equations have no solution. On the basis of the above problems, this paper proposes a novel hybrid algorithm that can guarantee efficient and accurate problem solving, so that any value input to the mechanism under the feasible motion space can have a solution, rendering it a universal algorithm for the 6-DOF forward kinematics problem.
The major contributions of our work are as follows: • An optimization method via an artificial bee colony (ABC) to optimize the BPNN (ABC-BPNN) is proposed. Good weights and thresholds are obtained through the process of the ABC's population iteration to prevent the training model from falling into a local optimum. • A modification of Newton's method for solving nonlinear equations with a singular Jacobian matrix is introduced. • We used QMn-M (a modification of Newton's method from Lv et al. [15]) combined with a simplified Newtonian iteration (SNR). • We used the length error threshold to reduce the frequency of ANN calls to improve the efficiency of real-time control.
First, ABC-BPNN [16] was employed to achieve outstanding prediction values, adjust them to the requirements of this hybrid strategy, and improve the overall strategy convergence and computing performance. Second, the QMn-M algorithm [15] can effectively solve the problem of a numerical method for solving forward kinematics suffering from a singular Jacobian matrix. Third, the global Newton-Raphson method with monotonic descent (GNRDM) proposed by Yang et al. [17] was combined with SNR to improve the efficiency of the algorithm [12]. Inspired by this technique, we used QMn-M combined with SNR to ensure the accuracy of numerical solutions. Moreover, it can reduce the computation time and iterations in overall algorithm. Lastly, the length error threshold ε 0 was designed into the overall algorithm process to determine whether the previous set of solutions Q pre satisfied the required value of the current iteration. If this is the case, it is directly used to carry out SNR. As a result, the frequent invocation of ANNs can be reduced. The overall operation efficiency and the problem-solving efficiency of ANNs in real-time motion control can be improved.
The structure of this paper is as follows: In Section 2, we describe the establishment of the kinematics model of GSP, including inverse kinematics problem analysis and the established FKP equations of GSP. The structure of ANNs and the proposed ABC-BPNN applied on the FKP is described in Section 3. Section 4 mainly describes the design and analysis of numerical algorithms. Section 5 outlines a comparative experiment that was conducted on a neural network and numerical algorithm. Lastly, Section 6 concludes the paper.

Inverse Kinematics
Compared with the FKP, the inverse kinematics problem (IKP) is relatively easy. The GSP consists of two platforms that are connected by several links, as shown in Figure 1. The base platform is fixed to the ground, and the mobile platform works as the end effector. The motion of six links causes the movement of the end effector.
According to Figure 1, the coordinates of the six lower hinge points are defined as A i (i = 1, 2, . . . . . . , 6), the coordinates of the six upper hinge points are B i (i = 1, 2, . . . . . . , 6). L i (i = 1, 2, . . . . . . , 6) is the length of the link between the connection point A i and B i at initial state. The radii of the six lower-and higher-hinge points are R 1 and R 2 respectively, as shown in Figure 2. The overall structure of 6-6 GSP and its initial state coordinate system are shown in Figure 2. For the base coordinate system {A}, defined by its base platform and its origin O A in the platform's geometric center. The moving coordinate system {B} is attached to the mobile platform, with its origin O B at the platform's geometric center.
The IKP is to calculate the value of joint variable L Oi by giving the pose of the moving coordinate system. The pose of moving coordinate system {B} is expressed as follows: where x, y, z are the surge, sway, and heave lines of the mobile platform respectively; and α, β, γ are the roll, pitch, and yaw angles of the mobile platform. Equation (3), which is the IKP's formula for the GSP, can be used to solve the displacement of i-th actuator using the geometric technique.
where L Oi denotes the displacement of the i-th link, L li indicates the length of i-th link when Q is the initial state, and R is the rotation matrix in the (Pitch-Roll-Yaw) Euler angle representation, as follows: where cα is cosα, sα is sinα, etc.

FKP Equations
The numerical algorithm's solution necessitates the construction of a set of nonlinear equations for iteration and root finding. To establish the FKP equations of GSP, IKP's formula of GSP, Equation (3) is needed. f i (Q) is the function mapping the pose from Cartesian space to the displacement of the i-th actuator L Oi in joint space. The FKP of GSP is to solve corresponding pose Q by giving L Oi ; we can construct Equation (5) to form a set of nonlinear equations whose constructed equations are Equation (6).
where L ni is the input value of the FKP (the displacement of the i-th actuator).

Artificial Bee Colony-Based BP Neural Network Algorithm
Artificial neural networks (ANNs) have good nonlinear mapping ability, and can realize the mapping relationship between the displacement of several links and the endeffector pose without considering the intermediate operation process. In the ANN-based IKP method, data were obtained by IKP for sample training. L Oi was the input, Q was the output, and only approximate output values could be obtained. The backpropagation neural network (BPNN) is a widely used artificial neural network in GSPs [9].

ABC-BP Neural Network
According to the mapping relation of f {L O1 , L O2 , L O3 , L O4 , L O5 , L O6 } → {α, β, γ, x, y, z}, a BPNN with 6 input nodes, 6 output nodes, and 1 hidden layer was designed as shown in Figure 3. Sigmoid function was chosen to be the activation function.
The BPNN model achieved promising performance and flexible network structure. However, it had several deficiencies that prevented it from completing this mission. To optimize the BPNN, the artificial bee colony (ABC) algorithm was used to ensure the stability of the deviation between the predicted and real values of the IKP's dataset of BPNN training, and to enhance the model training effect, namely, ABC-BPNN [18][19][20][21]. For the BPNN model, the initialization weights and thresholds are given as follows, respectively: Due to its benefits of fast convergence and global search, the ABC method is frequently used to solve optimization issues. Bees and food sources are the two most important components of the ABC algorithm. Scout, employed, and onlookers are the three sorts of bees. Scout bees are entrusted to search for food sources at random, while employed and onlooker bees are in charge of nectar mining [22,23]. Numerous studies showed that the ABC algorithm outperforms the genetic algorithm (GA), particle swarm optimization (PSO), and other methods in terms of accuracy and convergence rate [24,25]. In fact, the principle of the ABC algorithm is dealing with the problem of function optimization by simulating the nectar-gathering mechanism of actual bees. The ABC algorithm iterates continuously, compares the advantages and disadvantages of the problem, keeps the good individuals and eliminates the bad ones, and constantly approaches the global optimal solution. The ABC algorithm replaces the solution to the optimal honey of the colony with the solution to the optimal weights and thresholds of the neural network. The optimized BPNN is trained by the global optimal weights and thresholds provided by ABC. The whole ABC-BPNN algorithm is shown in Algorithm A1. The specific steps are as follows:

2.
Initialize the parameters of the ABC algorithm: self-variable dimension D, population number N (the number of employed bees is N 1 , the number of onlookers is N 2 , and the number of solutions is N 3 ), maximal cycle number G, and threshold of iteration number for scout bees limit.
where N input is the number of nodes in the input layer, N hidden is the number of nodes in the hidden layer, and N output is the number of nodes in the output layer. 3.
Employed-bee search. New solution V i (j) is generated by employed bees searching the field and is compared with old solution X i (j) by using the principle of the greedy algorithm. The strategy compares the size by calculating the fitness value. If the fitness value of the new solution is greater than that of the solution, the solution is replaced by the updated solution. Otherwise, the number of updates is f ailures + 1. The fitness value is calculated with Equation (10), and V i (j) is calculated with Equation (9).
where X i (j) is a random selection of other employed bees, different from themselves, where MSE i is the mean square error generated by training the i-th solution in the BPNN,ŷ im indicates the predicted value during training in the BPNN, y im indicates the actual value, and n indicates the total number of training samples.
where X t i represents the i-th individual in the population of the t generation. 4.
Onlookers search. Onlookers calculate possible value P i with Equation (13), and use the roulette to find a new solution among existing ones.

5.
Scout-bee search. The rapid convergence of employed and scout bees may lead to a decrease in the overall diversity of the population. In order to avoid the population entering the local optimum, the search mechanism of scout bees was designed on the basis of the ABC algorithm. If the number of update f ailures exceeds limit, the solution is discarded and replaced by a new solution generated by Equation (14), and the number of update f ailures is initialized to zero. In this paper, the value of limit is the product of N 3 and D. 6.
If the current cycle number is greater than maximal cycle number G, Step 3 is repeated.
Otherwise, the solution with the maximal fitness is output at the end of the training. 7.
According to the optimal solution provided by the ABC algorithm, the initial weights and thresholds of the BPNN are obtained. Subsequently, the ABC-BPNN model is trained and tested with the sample data to achieve the pose prediction of GSP.

Newton-Raphson Algorithm
In the traditional hybrid algorithm (BPNN and numerical algorithm), Newton-Raphson is used as the main method for numerical solutions. Some other improved Newtonian methods are used for the numerical iteration of optimization algorithms, as shown in Figure 4.
• Newton-Raphson method: The formula of the NR method for solving nonlinear equations of FKP is: where J is a Jacobian matrix, and its formula is: • Newton downhill method: The formula for the Newton downhill (NR-dh) method for solving nonlinear equations of FKP is: where λ is the descent factor; when λ = 1, it is the NR method. When ||F(Q n + 1)|| < ||F(Q n )||, λ remains the same. On the other hand, λ = 2 −t , (t = 0, 1, 2 . . . . . .).

Improved Newton's Method for Nonlinear Singular Equations
Neither the traditional nor the above method can solve the singular Jacobian matrix of FKP nonlinear equations. According to Equation (15), a new set of J is formed continuously in the process of the solving iteration, and its inverse is required. In the process of iteration, when J forms a singular matrix, Rank(J)< 6, its inverse cannot be solved. Therefore, the general solution formula of NR method should be: Q n+1 =Q n − J + Q n F(Q n ) (n = 0, 1, 2, 3 . . . . . .) (18) where continuous iterating pseudoinverse matrix J + slows down the solution process, affects the overall efficiency, and cannot guarantee accuracy of computation. To solve these problems, any value that is input into the feasible motion space of the mechanism can have a solution. Therefore, an improved Newton's method for solving nonlinear singular equations was introduced [15]. Lv et al. [15] combined the methods that were proposed in Kou et al. [26], and Howk et al. [27]; a diagonal matrix containing parameters was added to the singular Jacobian matrix to change the irreversible characteristics of the Jacobian matrix in the iterative process, and the convergence order of this new singular algorithm was second. The equation proposed in Kou et al. [26] is as follows: where µ i (n) ∈ R\{0}. when J is for the singular matrix, the diag µ i (n) F(Q n ) join can solve this problem. The formula of the improved Newton's method for nonlinear singular equations is: when n = 0, the formula is: when n ≥ 1, the formula is: where, when γ = 1 /2, it is the QMn-M algorithm. When γ = 0, it is the PC-M algorithm. When γ = 1, the two-step iterative Equations (20) and (21) change into the single-step iterative Equation (19), which is Newton's improvement on the standard. The recommended algorithm is the QMn-M algorithm at γ = 1 /2, which performed better on the FKP of GSP, with better iterations and total running time than those of the others. A comparison of several algorithms is shown in Table 1. Table 1. Analysis of the advantages and disadvantages of several numerical algorithms.

Numerical Algorithms Advantages Disadvantages
Newton-Raphson Easy to implement and fast convergence.
Sensitive to initial value, and the iteration is easy to diverge.
Newton downhill Reduces sensitivity to initial values.
The strong coupling of the Jacobian matrix reduces the iteration speed in FKP equations.

QMn-M or PC-M
It can solve the nonlinear equations of a singular Jacobian matrix, and only two function values need to be calculated in each iteration step 1 .
The amount of matrix calculation is increased, and more parameter settings are added. 1 With the QMn-M or PC-M algorithm, only two function values need to be calculated in each iteration step because, in Equation (21), F γQ n−1 + (1 − γ)Q * n−1 in the last iteration formula is in Equation (21). The only need is to calculate two functions: F(Q n ) and F (γQ n + (1 − γ)Q * n ).

Simplified Newton Iteration to Optimize QMn-M Algorithm
The QMn-M algorithm solves singular Jacobian matrix J and reduces its problemsolving speed. To ensure the global performance, and its operation speed and accuracy, we used the simplified Newton iteration (SNR) method to optimize the QMn-M algorithm.
Simplified Newton Iteration: where C is a constant Jacobian matrix. When it is not satisfied: ||F(Q n+1 )|| < ||F(Q n )|| (23) to renew the Jacobian matrix. The QMn-M and SNR algorithm can effectively improve the time efficiency of the algorithm and guarantee the accuracy of the value, as shown in Figure 5. Compared with the traditional hybrid algorithm, it can solve the global pose problem.

FKP of Real-Time Control Improvement
In the control of a 6-DOF parallel robot, a real-time solution of FKP is more practical. The control cycle and sampling interval are usually milliseconds. Therefore, the solution of FKP connected between control cycles is approximate. According to the deviation-driven algorithm on real-time FKP [12], and by Equation (24), Q pre (the value of the last solution) is reused, and ABC-BPNN calls are reduced. If the formula is satisfied, Q pre is used to directly carry out a simplified Newton iteration. Otherwise, ABC-BPNN is called to give an approximation.
where ε 0 is the length error threshold that controls whether to call ANNs. The whole hybrid algorithm is shown in Figure 6 and Algorithm A2. Figure 6. Algorithm for forward kinematics problem of 6-DOF. ∆E 1 refers to the difference between Q n+1 and Q n of each iteration, and ε 1 is the threshold for controlling the end of the iteration.

Experiment
The experimental process can be divided into three parts: • Part 1: Inverse kinematics solution. The corresponding data are obtained to form the dataset and normalized. Several experiments were conducted on the GSP (6-UCU) as shown in Figure 1. The positions of hinge points are given in Table 2. The maximal iterations N max were set to 20, and the permissible error was 10 −4 mm. The following experiments were carried out in the environment of MATLAB R2019b, and the operating system was Windows 10. When Q =[0, 0, 0, 0, 0, 0] T , the initial value of each links variable was:

Data Acquisition
A specific number of poses inside the workspace should be chosen before the training process, and the associated actuator displacements should be acquired by solving IKP with the corresponding actuator displacement Equation (3). The sample sets for training ANNs are formed by combining the two groups of data and then reversing their mappings. Essentially, ANNs use calculated data in a joint space as inputs, and corresponding pose data in Cartesian space as predicted outputs.

Neural Network Experiment
To train the ANNs, the network/data management toolbox in MATLAB R2019b was used. In the neural network comparison experiment, 300 groups of continuous data were collected from the machine. We divided the dataset at a ratio of 7:3. The training set consisted of 210 groups, and the test set consisted of 90 groups. We introduced a genetic algorithm-backpropagation neural network (GA-BPNN) to compare with our ABC-BPNN and the traditional BPNN. The BPNN, GA-BPNN, and ABC-BPNN parameter settings were the same, where the learning rate was 0.001, and set epochs were 200. In this experiment, the three neural network models were compared through metrics, as shown in Table 3. The average error of each variable predicted by the neural network model is given in Table 4. In this paper, we calculated the error through the numerical deviation of the corresponding data. The error calculation equation used is as follows: whereQ i represents the predicted value of a pose variable, Q i represents the real value of a pose variable, and m represents the individual number of the same variable inQ mi . Figure 7 shows the predicted values of the three models in a test set. A group of data were selected for comparison, and the pose was the predicted value of the neural network. The experimental results in the figure demonstrate that ABC-BPNN had a better prediction ability than that of the traditional BPNN or GA-BPNN. Figure 8 compares the error of the surge(x) predicted value of the neural network models with 10 groups of sample data. the error calculation of the surge variable was performed with Equation (26). Through the optimized neural network, its error value obviously decreased, as shown in the figure. By using a heuristic algorithm, the network model's prediction performance was enhanced, and ABC-BPNN was more able to adjust to changes in the process.   According to the metrics given in Table 3, the ABC algorithm optimized the initialization weights well, which improved the prediction ability of the BPNN model. Compared with GA-BPNN and BPNN, the RMSE index of ABC-BPNN decreased by 42.1% and 31.4%, respectively, and the R 2 performance of ABC-BPNN improved by 11.0% and 5.6%, respectively. The average error of the three models is given in Table 4; the maximal improvement of ABC-BPNN error optimization was 46.3%. By comparing the average error, the ABC algorithm could reduce the average error of various variables, and its effect was better than that of the genetic algorithm. Experimental results show that ABC-BPNN could reduce the dependence on the dataset and the prediction error of the neural network.

Numerical Algorithm Experiment
We used the data obtained by the neural network to conduct experiments in MATLAB R2019b. On a PC with a 2.30 GHz processor and 8 GB RAM, the FKP that corresponded to the obtained data was solved. An error threshold ε 2 was used to control the QMn-M algorithm in the experiment, and an error threshold ε 1 was subsequently used to control the SNR algorithm. The number of iterations and overall operation time fluctuate depending on ε 1 . The design of the error threshold can be changed with specific data until the value that renders it the most efficient is found. Error threshold ε 1 is used to govern the end of the SNR algorithm and the entire process. Its value is generally accurate to four decimal places (as shown in Table 5, The precise level was 10 −4 ).
After the neural network experiment, the predicted values were iterated numerically. In this section, the average numbers of iterations, required operation time, and different iterative methods are compared. Table 5 shows the comparison of seven iteration methods, giving the average number of iterations and the required time. As shown in Table 5, the required time and average number of iterations of the proposed QMn-M and SNR were reduced, and its performance was better than that of the same second-order NR algorithm after being optimized by SNR. In Table 6, the traditional NR, and NR-dh with QMn-M and PC-M are compared. QMn-M and PC-M could solve nonlinear singular equations, and the advantages of this algorithm in solving the forward kinematics problem are demonstrated. Experiments showed the feasibility of QMn-M and PC-M in solving singular matrix data. (µ i (n) , λ i (n) ) about Q * n and Q n+1 in Equation (21) were set to be the same parameters. In the QMn-M experiment, different (µ i (n) , λ i (n) ) values caused different levels of performance, and unsuitable (µ i (n) , λ i (n) ) values may render it divergent and unsolvable. Figure 9 shows the comparative experiment of the operation time required by four algorithms in 10 groups of data. As described in the Figure 9, QMn-M and Newton-Raphson had the best performance, followed by PC-M, and NR-dh had the worst performance. Figure 10 is a comparison of QMn-M, QMn-M and SNR, and the traditional NR method. It shows that SNR combined with a numerical algorithm could greatly improve the efficiency of the algorithm and better guarantee accuracy. Furthermore, we found the percentage improvement performance of other iterative algorithms by comparing with QMn-M and SNR, which is illustrated in Table 7. Compared to NR-hd, the maximal percentage improvement in performance for the average number of iterations and required time was 68.8% and 70.60%, respectively. Compared to NR, the minimal percentage improvement in performance for the average number of iterations and required time was 14.4% and 13.9%, respectively.

Conclusions
This paper mainly optimized the traditional hybrid algorithm for FKP. The ABC-BPNN model was used to train the IKP's value to obtain a good initial value, and the QMn-M and SNR algorithm was used to iteratively calculate the excellent initial value. The comparison chart in the numerical algorithm experiment shows that the QMn-M and SNR algorithm was superior to other methods in terms of iterations and operation time. It also has the ability to calculate a singular matrix. For continuous samples, its solution does not need to continuously call ABC-BPNN, which greatly reduces the overall operation time required, and makes it convenient to solve on the real-time platform. This new hybrid algorithm can not only ensure that the GSP can have solutions in feasible space, but also meet the speed and accuracy of real-time computation. This algorithm also has portability, which is beneficial to solve the FKP for different structural platforms.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A
The whole ABC-BPNN algorithm for the pose estimation of Gough-Stewart platforms is as follows.
Algorithm A1 ABC-BP neural network for FKP.

Input:
The displacement of current GSP, L Oi ; Output: The pose of current end-effector Q prediction . 1: Establish BP neural network model. 2: Initialize the parameters of the ABC algorithm. 3: Employed bees generate candidate solution V i (j), V i (j) shown in Equation (9). 4: Fitness value comparison in Equations (10) and (12). If New > Old, then Old ← New, else f ailures + 1. 5: Compute the possible value P i of the solution, shown in (13), search according to P i to generate new solutions. 6: Fitness value comparison in Equations (10) and (12). If New > Old then Old ← New, else f ailures + 1; 7: If f ailures > limit, then scout bees initial solution, shown in Equation (14); else, move to Step 8. 8: If current cycles > G, then save the optimal fitness value; else, move to Step 3. 9: Get the optimal weights and thresholds and train the BP neural network.
The whole hybrid algorithm to solve forward kinematics problem for Gough-Stewart platforms is as follows.
Algorithm A2 Hybrid algorithm for forward kinematics problem of 6-DOF.