Journal of Robotics
Volume 2019 (2019), Article ID 9531387, 10 pages
https://doi.org/10.1155/2019/9531387

Research Article

Construction and Layout of Coordinated Multirobots System Based on Parallel Thought

Wanqiang Xi,1,2 Yaoyao Wang,¹ Yadong Ding,¹ and Bai Chen¹

¹College of Mechanical and Electrical Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
²College of Mechanical Engineering, Jiangsu University of Technology, Changzhou 213001, China

Correspondence should be addressed to Yaoyao Wang; yywang_cmee@nuaa.edu.cn

Received 1 October 2018; Revised 4 December 2018; Accepted 19 December 2018; Published 3 January 2019

Academic Editor: Yangmin Li

Copyright © 2019 Wanqiang Xi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Aiming at the problems that are the construction and the layout optimization of the coordinated multirobots system, the concept of an equivalent parallel robot (EPR) is proposed based on parallel thought in this paper. That is, each robot is regarded as a branched chain of the EPR in coordination. Thus, it can be converted into the problem for the type synthesis and the layout of the parallel robot. Firstly, the definition of robotic characteristic (C) set and the corresponding operational rule are given. Then, the minimum occupied area, cycle time, and condition number are used as the optimized objectives. The weight ratio is used to express the importance of the three objectives to meet the requirements of different scenes. Afterwards, based on the theory of set, the three translational (3T) operational requirements that can be extended to other degree of freedom (DOF) are analyzed. Finally, the feasibility of the proposed methods is verified by the DOF analysis and simulation experiment. Moreover, the multiobjective layout problem is optimized based on the artificial bee colony (ABC) algorithm. The results show that the constructed robot combinations and the layout optimization method are satisfactory.

1. Introduction

The development of robotics has continuously improved the work capabilities of robots, and the field of robotic applications is constantly expanding [1, 2], for example, assembly work in automated factories, deep-sea operations, fault handling in the nuclear industry, and operational tasks in space. On the one hand, due to the complexity of the task, when a single robot is difficult to complete the task, people hope to accomplish this through coordination of multirobots [3]. On the other hand, people also hope to improve the efficiency of the robots system during the operation through the coordination of the multirobots [4]. Furthermore, when the working environment changes or the partial system fails, multirobots can still accomplish their scheduled tasks through their own fault tolerance [5]. In conclusion, compared with a single robot, multirobots system has higher stiffness, accuracy, load capacity, and flexibility [6, 7].

As a new form of robotic application, the coordinated multirobots have increasingly attracted the interest and concern of academic both at home and abroad. For the different operational tasks, different DOF may be required, so the robots need to be chosen to meet the operational DOF. Although the combinations of robots can meet the operational requirements, the unreasonable layout may reduce the efficiency of the system, occupy more areas, and have poor motion performance [8–10]. The previous work of the authors has done some research on the type synthesis of the coordinated multirobots system [11]. However, the set in the description of the position and posture of a robot is not detailed enough, and the layout optimization is still lacking, and it cannot constitute the overall design of the coordinated multirobots system. These problems will be studied in this paper.

The layout of the multirobots system is to optimize its layout position according to different operational tasks after the branched robots are determined. The positions of the branched robots are adjusted for the different requirements of the task to achieve the best layout effect. However, for a parallel robot, it is generally achieved through dimension synthesis to meet the operational requirements; thus its branches will not change once they are laid out. At present, there are few relevant literatures on the layout of coordinated multirobots system. It mainly depends on the actual experience and cannot achieve the best results. Michele Gadaleta et al. proposed the method for the energy-optimal layout design of robotic work cells, which numerically retrieved the industrial robot base position to achieve the minimum energy consumption [12]. Zhen Yang Lim et al. presented a layout optimization method for multirobots cellular manufacturing system and compared the effects of several algorithms on the optimization results [13, 14]. However, either the optimized objectives are too single or the result of optimization is only an approximate areas. It is not possible to give an accurate layout according to the different operational tasks and the importance of the optimized objectives. This paper proposes to use the minimum occupied area, cycle time, and condition number as the optimized objectives. The weight coefficient is used to express the importance of the three objectives to meet the requirements for different scenes of the tasks, and the artificial bee colony (ABC) algorithm is used to optimize the multiobjective layout problem.

The main contributions of this paper are as follows. (a) The problem for the construction and layout optimization of a coordinated multirobots system is transformed into the type synthesis and layout of a parallel mechanism based on the parallel thought. (b) The concept of C set is presented, which can conveniently represent the end kinematic features (for the robot selection of multirobots) and the position and posture of the robot. (c) The multiobjective layout optimization method with weighted ratios is proposed, and the minimum occupied area, cycle time, and condition number are used as the optimized objectives. The remainder of this paper is organized as follows. In Section 2, the basic concept of robotic set is given and corresponding intersection rule is analyzed. In Section 3, the minimum occupied area, cycle time, and condition number are taken as optimized objectives for the layout. The selection and layout of a specific task (3T) are mainly designed in Section 4. Lastly, Section 5 draws the conclusion.

2. Characteristic Set and Its Intersection Rule

The characteristic () set is used to describe the end-feature of a mechanism and the position and the posture of its base [11], as shown in

where represents the translational features of robot’s end-effector in the direction of , , and and the position of its base. represents the rotational features and the posture of its base. They are just as shown in Figure 1.

(a) Transnational features

(b) Rotational features

Figure 1: Description for the end-features of the set.

The meaning of the numerical values of the set’s elements is described in the literature [11]. The following examples are illustrated.

For example, there is a four-DOF robot with two translational features in the and directions and two rotational features rotating around the and axes, which belongs to a 2T2R robot. If the world coordinate system is set as the base coordinate system of the robot, its set can be written as . If the position and posture of the robot’s base coordinate system is not 0 relative to the world coordinate system (neither the coordinates nor the RPY angles are 0), its set can be written as .

The operation of set includes intersection and union, and the intersection is suitable for parallel mechanism [15, 16]. It is can be defined in

where stands for the intersection operator.

3. Layout of Multirobots

The layout optimization of a coordinated multirobots system can be considered as a parallel robot whose branched chains’ positions are changed to achieve the minimum occupied area, cycle time, and condition number.

3.1. Occupied Area

Reasonable layout will reduce the occupied area of the completely coordinated multirobots system in the workshop [13]. For a parallel robot, the objective can be seen as the minimum of the static platform area by changing the position and the posture of its branches to meet the required tasks. The calculation method of the occupied area is that the graphic area formed by the plane where the base of each robot is located. Moreover, it can be calculated based on the value of each robot’s set.

3.2. Cycle Time

In order to be more efficient, the production shop requires less time for the whole process. Therefore, it takes minimal time for EPR to complete the task. The cycle time of the EPR is determined by its rotation angle and the maximum angular velocity. The following equations are specified.

Assume that the angular velocity of the jth joint of No.i robot is . In a cycle time T, the continuous rotation angle of the jth joint is , and the time of the jth joint is

However, the running time of each robot’s joint may be different. Now set the running time of No.i robot in the cycle time as

During the cycle time T, the simultaneous tasks run for

Therefore, the total cycle time is the sum of the time for each independent task; that is,

Similar to the EPR, it is required to arrange the position and the posture of each branch for a parallel robot reasonably. When the task is completed, the branch with the longest consumption time is considered as the minimum cycle time of the parallel robot. In the following sections, the specific example will be given.

3.3. Condition Number

For serial robots, the statics performance index of robot’s flexibility is generally used by Yoshikawa’s manipulability index [17]. That is,

where is the Jacobian matrix of the robot and is the inverse of the Jacobian matrix. However, for a parallel robot, the condition number is often used as its performance index [18]. In a highly coupled multirobots cooperative system, such as multirobots transporting, it can be seen as an EPR, and its condition number is [19]

where is the norm of the Jacobian matrix and for the 2-norm, the condition number can be written as [20]

where and are the maximum and minimum singular values of the Jacobian matrix , respectively, and the range of the condition number is .

When designing a parallel robot, the value of the condition number should be 1 as far as possible [21]. Because when , the singular values are equal, and the arm is in an isotropic configuration with the highest flexibility.

3.4. Weight Fitness

Layout optimization means that the reasonable layout is selected to meet the different requirements. In this paper, the minimum occupied area, cycle time, and condition number are taken as the optimized objectives of the layout. However, in actual situation, one or two combinations may be dominant. Therefore, the proportion of these three evaluation indicators is particularly important, and the weight ratio will be used to express the importance of the three indicators.

The fitness of the area is expressed as

where is the occupied area and is the minimum occupied area. The fitness of the operation time is expressed as

where is the cycle operation time and is the minimum cycle operation time. The fitness of condition number is expressed as

where represents the condition number of the coordinated multirobots system and represents the minimum condition number.

According to different operational requirements, a reasonable weight ratio is designed to represent the importance of the system’s indicators:

where , , and are the weight coefficients and . Their values are determined by their importance. The more important a certain indicator is, the greater its weight coefficient is.

The values of , , and in the above three equations are obtained by ignoring the other two indicators. For example, when is calculated, only the area indicator is considered, and the requirement of the operation time and the condition number are abandoned, that is, in the case where = 1, b = 0, and c = 0.

To optimize the above-mentioned three objectives, the following assumptions are needed:

(1)If there is a nonstatic robot in the system (the base can move), its maximum moving area is set as the static area.(2)The moving point is set as the base coordinate point of the robot.(3)Each joint runs at the maximum speed.

4. Case Study

For example, a procedure of a factory’s production line needs to transport the goods with heavy or bulk from worktable A to worktable B (only involving the three-dimensional translation of the space). However, many small-load robots with the same or different types cannot complete this task alone. Considering the cost savings, no new heavy-load robot will be purchased. At this time, it is necessary to filter and reorganize the existing robots in the factory and build a coordinated multirobots handling system composed of multiple robots.

Next, the above-mentioned type synthesis method and weight multiobjective optimization method will be used in the chain selection and the layout design of the EPR separately.

4.1. The Selection of the 3T EPR’s Chains

The set of the 3T EPR’s end-manipulator is expressed as follows:

The literature [11] shows that the set of the EPR’s end-manipulator is determined by

where , , and are the sets of the EPR’s three-branched robots, respectively, and the set of the branched robots can be expressed as

It can be seen from (15) that the set of the 3T EPR can be obtained by the intersection of the sets for the above three-branched robots. As seen in Table 1, there are sixteen combinations of the branched chains, and the intersections must satisfy (2).

Table 1: The combinations of the branched robots with 3T motion.

Chain1
Chain2
Chain3
Combinations
Structure

symmetry

asymmetry

symmetry

asymmetry

symmetry

The classification and the condition of the combinations from Table 1 are listed in Table 2. According to the Table 1, the intersection just as has fifteen combinations and the proper subset of intersection just as has 9 combinations.

Table 2: Classification for the combinations of the branched robots.

Objectives
Combinations
Conditions

No requirements

Two are in , and

Three ∥ cannot appear in , , and

Three ∥ cannot appear in , , , and

and

Proper subset is

∥∥

Three is ∥ in , , and

∥

Three is ∥ in , , , and

Two is ∥ in , ,

No requirements

Since each branched chain of EPR is an independent robot, it can be a parallel or serial robot. The common robots with 3 to 6 DOF are listed in Table 3.

Table 3: Common robots with 3 to 6 DOF.

Branched chain
Category
Robots

C3T
serial
P₁P₂P₃, R₁R₂R₂

parallel
Delta, 3-PRRR, 3-RPS

C3T1R
serial
R₁R₂R₂R₃

parallel
4-PU₁U₂

C3T2R
serial
R₁R₂R₂R₃R₄

parallel
5-RPUR

C3T3R
serial
Puma560 (6R)

parallel
Stewart (6-UPS)

Then the combinations of the robots are selected to meet the task requirements (3T), and the choice of this paper is the combination of . That is, a parallel robot 3-PRRR and two six DOF serial robots are used as the EPR’s branched chains to transport the goods, as shown in Figure 2. The three pneumatic hands are fixed under the moving platform of the 3-PRRR and the end of the serial robots to transport triangular goods. In this process, the three hands contact the goods closely without any change in position and posture between them. At this moment, the goods can be seen as a moving platform for the EPR, and the three robots are taken as the three-branched chains for the EPR. The screw theory [22] is used to analyze the DOF of EPR in the Equations (20), (21), and (24).

Figure 2: Three robots transporting triangular goods in coordination.

As shown in Figure 2, the parallel robot 3-PRRR is composed of a moving platform and a static platform, which are connected by the same three-branched chains. It is not difficult to see that the kinematic screw system of the parallel robot 3-PRRR is

Its corresponding reciprocal screw is

The set of the parallel robot 3-PRRR can be written as (22). If the base coordinate system of 3-PRRR is set to the world coordinate system o-xyz, its set can be expressed as (23).

Similarly, the kinematic screw system of the six DOF serial robot is

For a 6-DOF serial robot, it is obvious that the set of its end kinematic features can be described as

where the six elements are nonzero. According to (15), the end-features set of EPR can be described as

4.2. The Layout Design of the 3T EPR

When multirobots work in cooperative task, the layout will have a direct impact on the flexibility, operation time, and workspace of the multirobots system. Therefore, a reasonable layout is necessary. The method proposed in this paper is to use the minimum cycle time, occupied area, and condition number of the EPR as the optimized objectives, and their weight coefficients are chosen according to specific requirements.

As shown in Figure 2, the coordinate systems , and are set to the base coordinates of No. 1, No. 2, and No. 3 robots, respectively. The Cartesian coordinate systems , and are set to the moving coordinates of No. 1, No. 2, and No. 3 robots, respectively. is set to the coordinate system for the triangular goods, and is also set to the world coordinate system of the EPR system. Assume that the starting coordinate A of the No.1 robot’s manipulator is (0.2, 0.2, 0.3) and the ending coordinate B of its manipulator after the translation is (0.4, 0.4, 0.5). The DH parameters of the six DOF serial robot are listed in Table 4, and the corresponding coordinate systems of its linkages are shown in Figure 3.

Table 4: Six DOF serial robot parameters of links.

(rad)
(m)
(m)
(rad)

1
θ₁(-π/2)
0.325
0
0

2
0
0.1
0
-π/2

3
θ₃(-π/2)
-0.1
0.300
0

4
0
0.405
0
-π/2

5
0
0
0
π/2

6
0
0
0
-π/2

Figure 3: Six DOF serial robot links coordinate system.

The task requirement is to move the goods from point A to point B and optimize the layout of the EPR system with the minimum occupied area, cycle time, and condition number.

The coordinate system of the EPR system is established and the schematic diagram is shown in Figure 4. The circular modules with colors represent the base of the three robots, and the triangle modules represent the hands of the robots. The yellow pattern indicates the 3-PRRR parallel robot, and the green indicates the 6-DOF serial robot. The coordinate system of EPR system is projected onto the xy plane of the world coordinate system, and the coordinate of each module is as shown in Figure 4. It is worth noting that, in the position layout of the robot, the base coordinate system of the fixed 3-PRRR is set as the world coordinate system, and only the positions of two six DOF serial robots can move.

Figure 4: The coordinate systems diagram of the EPR layout.

The mathematical form of the three optimized objectives is described as follows:

where represents the area optimized function. represents the time optimized function. represents the condition number optimized function. In addition, the circle of radius and the rotation angle θ are the movable range of the serial robots.

The artificial bee colony (ABC) algorithm [23, 24] is used to optimize the multiobjective layout parameters, and its detailed procedure can be depicted in Figure 5. Since the EPR motion form belongs to the type of 3T, the two six-DOF serial robots (No. 2 and No. 3) keep the posture in the layout, as shown in Figure 4.

Figure 5: The flow chart of ABC algorithm.

Assume that the maximum velocity of the parallel robot 3-PRRR’s prismatic joint is v=2.5m/s, the maximum velocity of the serial robot’s rotational joint is ω=2πrad/s, and they run at the maximum velocity. Since the condition number of 3-PRRR is always 1 and the operation status of two serial robots is the same, the condition number of EPR is equivalent to the condition number of a serial robot. As long as the condition number of the serial robot is minimal, the condition number of EPR will be minimal.

Figures 6–8 show the plane distribution diagrams between cycle time, occupied area, and condition number. Figure 9 is a spatial distribution diagram of the three objectives set as spatial coordinate axes. The method of Monte Carlo is used to plot the diagrams in this paper. It can be seen that the scatter plots between the three objectives are not evenly distributed. The cycle time mainly concentrates on the range of 0.35 to 0.55, and the condition number mainly concentrates on the range of 8.5 to 11. In general, the condition number increases when the occupied area increases.

Figure 6: Scatter plot of cycle time and occupied area.

Figure 7: Scatter plot of occupied area and condition number.

Figure 8: Scatter plot of cycle time and condition number.

Figure 9: Scatter plot of cycle time, occupied area and condition number in the space.

Different layout will affect the values of the three objectives, and it is particularly important to choose the layout location that meets the objective requirements. Next, the ABC algorithm will be used to search for the optimal solution according to the importance of its optimized objectives. The results obtained are shown in Table 5, and the corresponding set of each branched robot in the EPR is shown in Table 6. When a=1, b=0, and c=0, the minimum occupied area =1.5500m² can be obtained. Similarly, when a=0, b=1, and c=0, the minimum cycle time =0.0903s, and when a=0, b=0, and c=1, the condition number =8.4253. As seen in Table 5, the values of the objectives decrease as the weight coefficient increases, which is consistent with the results of multiobjective optimization. From Table 6, the C sets of each group of branched robots represent their end kinematic features, the position, and the posture of their layout, which greatly facilitates the selection and layout of the multirobots.

Table 5: Optimization results with different weight coefficients.

Number

area
time
condition

1
1
0
0
0.1500
0
1.5500
-
-

2
0.6
0.2
0.2
0.1283
5.8478
1.5837
0.1281
8.5469

3
0.33
0.34
0.33
0.1185
5.5168
1.6146
0.1154
8.4512

4
0.2
0.6
0.2
0.0773
0.2950
1.6260
0.1349
8.7242

5
0
1
0
0.1500
4.6911
-
0.0903
-

6
0.2
0.2
0.6
0.1263
5.7336
1.5923
0.1239
8.4776

7
0
0
1
0.1500
5.5328
-
-
8.4253

Table 6: The set of each robot under different weight coefficients.

Number
No.1 (C)
No.2 (C)
No.3 (C)

1

2

3

4

5

6

7

5. Conclusions

In this paper, the concept of EPR is proposed based on the parallel thought. The problem for the construction and layout optimization of a coordinated multirobots system is transformed into the selection and layout of the branched chains of a parallel robot. Firstly, the concept of set and its operational rule are presented, which can conveniently represent the end kinematic features and the position and posture of the robot. The selection scheme of EPR satisfying the requirements is constructed based on the intersection rule of the set. After that, the ABC algorithm is used to optimize the selected robots according to the principle of minimum occupied area, cycle time, and condition number. Based on the importance of the three optimized objectives, a multiobjective optimization for the multirobots system with weighted ratios is proposed. The occupied area, cycle time, condition number, and the locations of the robots under different weight coefficient can be obtained by the simulation experiments. The results show that the proposed methods are feasible.

In future, the work may involve the selection and layout of multirobots in more complex situations. For example, consider the situations of more tasks, more operational sequences, and more robots.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This research was supported by the National Natural Science Foundation of China (Grants nos. 51575256, 51705243, and 51706098) and the Natural Science Foundation of Jiangsu Province (BK20170789).

References

L. Ding, X. Li, Q. Li, and Y. Chao, “Nonlinear Friction and Dynamical Identification for a Robot Manipulator with Improved Cuckoo Search Algorithm,” Journal of Robotics, vol. 2018, Article ID 8219123, 10 pages, 2018.

M. Hosseini, R. Meattini, G. Palli, and C. Melchiorri, “A Wearable Robotic Device Based on Twisted String Actuation for Rehabilitation and Assistive Applications,” Journal of Robotics, vol. 2017, Article ID 3036468, 11 pages, 2017.

P. R. Culmer, A. E. Jackson, S. Makower et al., “A control strategy for upper limb robotic rehabilitation with a dual robot system,” IEEE/ASME Transactions on Mechatronics, vol. 15, no. 4, pp. 575–585, 2010.

S. Manzoor, S. Lee, and Y. Choi, “A Coordinated Navigation Strategy for Multi-Robots to Capture a Target Moving with Unknown Speed,” Journal of Intelligent & Robotic Systems, vol. 87, no. 3-4, pp. 627–641, 2017.

D.-H. Lee, H. Park, J.-H. Park, M.-H. Baeg, and J.-H. Bae, “Design of an anthropomorphic dual-arm robot with biologically inspired 8-DOF arms,” Intelligent Service Robotics, vol. 10, no. 2, pp. 137–148, 2017.

S. Erhart and S. Hirche, “Internal Force Analysis and Load Distribution for Cooperative Multi-Robot Manipulation,” IEEE Transactions on Robotics, vol. 31, no. 5, pp. 1238–1243, 2015.

W. Wan and K. Harada, “Developing and Comparing Single-Arm and Dual-Arm Regrasp,” IEEE Robotics & Automation Letters, vol. 1, no. 1, pp. 243–250, 2016.

S. Angra, R. Sehgal, and Z. Samsudeen Noori, “Cellular manufacturing-A time-based analysis to the layout problem,” International Journal of Production Economics, vol. 112, no. 1, pp. 427–438, 2008.

M. Pellicciari, G. Berselli, F. Leali, and A. Vergnano, “A method for reducing the energy consumption of pick-and-place industrial robots,” Mechatronics, vol. 23, no. 3, pp. 326–334, 2013.

I. Suemitsu, K. Izui, T. Yamada, S. Nishiwaki, A. Noda, and T. Nagatani, “Simultaneous optimization of layout and task schedule for robotic cellular manufacturing systems,” Computers & Industrial Engineering, vol. 102, pp. 396–407, 2016.

W. Xi, B. Chen, Y. Wang, and F. Ju, “Type synthesis of coordinated multi-robot system based on parallel thought,” Transactions of the Canadian Society for Mechanical Engineering, vol. 42, no. 2, pp. 164–176, 2018.

M. Gadaleta, G. Berselli, and M. Pellicciari, “Energy-optimal layout design of robotic work cells: Potential assessment on an industrial case study,” Robotics and Computer-Integrated Manufacturing, vol. 47, pp. 102–111, 2017.

Z. Y. Lim, S. G. Ponnambalam, and K. Izui, “Multi-objective hybrid algorithms for layout optimization in multi-robot cellular manufacturing systems,” Knowledge-Based Systems, vol. 120, pp. 87–98, 2017.

Y. L. Zhen, S. G. Ponnambalam, and K. Izui, “Nature inspired algorithms to optimize robot workcell layouts,” Applied Soft Computing, vol. 49, no. C, pp. 570–589, 2016.

X. Meng and F. Gao, “Basic problems and criteria for synthesis of robotics,” International Journal of Mechanisms and Robotic Systems, vol. 1, no. 1, pp. 35–48, 2013.

F. Gao, J. Yang, and Q. J. Ge, “Type Synthesis of Parallel Mechanisms Having the Second Class GF Sets and Two Dimensional Rotations,” Journal of Mechanisms and Robotics, vol. 3, no. 1, pp. 623–632, 2011.

T. Yoshikawa, “Manipulability of robotic mechanisms,” International Journal of Robotics Research, vol. 4, no. 2, pp. 3–9, 1985.

J. K. Salisbury and J. J. Craig, “Articulated hands: force control and kinematic issues,” International Journal of Robotics Research, vol. 1, no. 1, pp. 4–17, 1982.

M. H. Abedinnasab, Y. J. Yoon, and H. Zohoor, “Exploiting Higher Kinematic Performance - Using a 4-Legged Redundant PM Rather than Gough-Stewart Platforms,” Serial and Parallel Robot Manipulators - Kinematics, Dynamics, Control and Optimization, 2012.

B. Zi and S. Qian, Design, Analysis and Control of Cable-suspended Parallel Robots and Its Applications, Springer Nature, 2017.

L. Wang, W. Rong, L. Qi, and Z. Qin, “Optimal design of a 6-DOF parallel robot for precise manipulation,” in Proceedings of the 2009 IEEE International Conference on Mechatronics and Automation, ICMA 2009, pp. 3933–3937, August 2009.

Y. Fang and L.-W. Tsai, “Structure synthesis of a class of 4-DoF and 5-DoF parallel manipulators with identical limb structures,” International Journal of Robotics Research, vol. 21, no. 9, pp. 799–810, 2002.

W. Xi, B. Chen, L. Ding et al., “Dynamic Parameter Identification for Robot Manipulators with Nonlinear Friction Model,” Transactions of the Chinese Society for Agricultural Machinery, vol. 48, no. 2, pp. 393–399, 2017.

L. Ding, H. Wu, and Y. Yao, “Chaotic Artificial Bee Colony Algorithm for System Identification of a Small-Scale Unmanned Helicopter,” International Journal of Aerospace Engineering, vol. 2015, Article ID 801874, 11 pages, 2015.