Articles | Volume 13, issue 1
https://doi.org/10.5194/ms-13-123-2022
https://doi.org/10.5194/ms-13-123-2022
Research article
 | 
25 Feb 2022
Research article |  | 25 Feb 2022

Multi-objective optimization of a redundantly actuated parallel robot mechanism for special machining

Haiqiang Zhang, Jianglong Tang, Qing Gao, Guohua Cui, Kan Shi, and Yan'an Yao
Abstract

In order to improve the accuracy and efficiency of special machining for a complex surface, a 2RPU-2SPR (where R, P, U, and S stand for revolute, prismatic, universal, and spherical joints, respectively) over-constrained redundantly actuated parallel robot mechanism is proposed. And six performance evaluation indexes are established to ensure the working performance including workspace, motion/force transmission efficiency, stiffness, dexterity, energy efficiency, and the inertia coupling index. Furthermore, a collaborative optimal configuration algorithm is conducted based on an orthogonal experimental design algorithm and a multi-objective particle swarm optimization algorithm. On the basis given above, a simulation analysis of a multi-objective optimization is conducted. Compared with two traditional, intelligent optimization algorithms of a multi-objective particle swarm optimization algorithm and an orthogonal experimental design method, the improved collaborative multi-objective optimization algorithm has a better optimization effect.

Dates
1 Introduction

With the continuous development of the manufacturing industry, machine tools play an important role in the processing of large and complex specially shaped work pieces. The requirements of parts that process quality and accuracy are also increasing, and processing equipment is developed in the mode of traditional equipment to advanced equipment and high-level equipment. Zhang et al. (2019a) presented a novel redundantly actuated parallel manipulator that can be employed to form a five-axis hybrid kinematic machine tool for large heterogeneous complex structural component machining in the aerospace field. Yang et al. (2019) invented a new three degrees of freedom (DOF) asymmetric translational parallel manipulator that is adapted to pick-and-place operation under heavy load. Wu et al. (2019) provided formulations of shape singularity and illustrate the application of shape singularity in an adjustable compliance mechanism design. Therefore, the research on redundantly actuated parallel mechanism is great valuable for special machining. Zhang et al. (2018) proposed a novel mobile serial–parallel mechanism with legs for in-pipe use and analyzed the singularities of the parallel mechanism. In addition, the parallel mechanism can achieve a continuous and equivalent rotation.

However, an over-constrained redundantly actuated parallel mechanism is a special parallel mechanism which has the characteristics of multi-variability, strong coupling, and complex nonlinearity. Reveles et al. (2016) presented a method for the joint trajectory planning of kinematically redundant parallel manipulators, and the working modes could be chosen and combined to complete the whole path by applying this approach. Furthermore, an accurate dynamic model of over-constrained redundantly actuated parallel mechanisms is extremely complex because of the coupling relationship between the branches of the parallel mechanism. Zhang et al. (2019c) simplified the dynamics mathematical model of the parallel manipulator by employing D'Alembert's principle and the principle of virtual work. Russo et al. (2018) solved the problems of inverse and forward kinematic for a parallel mechanism with 3-UPR (where R, P, and U, stand for revolute, prismatic, and universal joints, respectively) architecture and computed the Jacobian matrix to evaluate the singular positions of the end-effector. Chen et al. (2012) derived the velocity and the force Jacobian matrix of the parallel mechanism and analyzed the distribution of different dexterity measures in workspace, respectively. Tian et al. (2020) presented a general approach for the error modeling of a lower-mobility parallel mechanism based on screw theory, and then the geometric errors affecting the compensable accuracy of the end-effector were separated using the properties of dual vector space. Zou et al. (2021) proposed a 3 DOF parallel robot without rotation ability and calculated the distributions of maximum joint effect of force and speed. At present, many researchers at home and abroad have carried out abundant related researches on multi-objective performance optimization in different application fields. Enferadi et al. (2017) proposed a kinematics index that is called global workspace conditioning index and adopted a multi-objective evolutionary algorithm to find the optimal solutions and Pareto frontier. Asker et al. (2019) conducted dimensional synthesis based on maximizing the workspace index and the global force transmission and stiffness indices to optimize the device dimensions. Zhang et al. (2017) reflected a competitive mechanism based on multi-objective particle swarm optimization. Zeng et al. (2021) adopted an improved calculation of crowding distance, which substituted a traditional addition operation with multiplication operation. In addition, differential evolution was concatenated to diversify the population, and the key parameters of differential evolution had been assigned different adjusting strategies to further enhance the overall performance. A novel hybrid system was proposed by Xiang et al. (2015) based on an improved gravitational search algorithm which has good performance in solving function optimization problems. In order to accelerate the convergence to the Pareto frontier, Liu et al. (2018) presented a dynamic multiple population particle swarm optimization algorithm based on the decomposition and prediction. Wu et al. (2020) proposed a new algorithm to improve the convergence speed of a dynamic multi-objective optimization algorithm by combining the evolutionary algorithm and the dynamic strategy. Houssein et al. (2022) adopted a multi-objective optimization algorithm, based on the slime mold algorithm, which had the ability to provide better solutions in terms of the Pareto sets proximity and inverted generational distance in decision space indicators compared with the other algorithms.

In short, reasonable structural parameters are the basis of mechanical design. The quality of the mechanism parameters directly determines whether the mechanical performance of the mechanism is excellent or not. However, there are many design variables and solving processes of the objective function when designing the structural parameters of over-constrained redundantly actuated parallel mechanism. The research mechanism parameters cannot be optimized effectively if one simply considers a certain performance index when optimizing the parameters of parallel mechanism. The main contribution of this paper is as follows: a novel, redundantly over-constrained 2PRU-2SPR parallel mechanism is presented, and its structure is introduced in detail. Simultaneously, the kinematics and dynamics analysis of a 2RPU-2SPR parallel mechanism is deduced. What is more, in order to obtain more reasonable design structural parameters and improve the working performance of the 2RPU-2SPR parallel mechanism, six performance evaluation indexes are established to ensure the working performance, including workspace, stiffness, motion/force transmission efficiency, dexterity, energy efficiency, and the inertia coupling index. Furthermore, a multi-objective particle swarm optimization configuration algorithm based on orthogonal experiment response design is proposed in this paper, aiming at special machining complex surfaces. Compared with a traditional multi-objective optimization algorithm, the improved, collaborative multi-objective optimization algorithm has a better optimization effect.

2 2RPU-2SPR over-constrained redundantly actuated parallel mechanism

2.1 Description of each kinematic joint

An over-constrained 2RPU-2SPR redundantly actuated parallel mechanism is proposed in order to meet the machining requirements of complex surfaces and improve the machining quality. As shown in Fig. 1, the mechanism is mainly composed of a fixed platform, moving platform, end-effector, and two joints of actuating branches with the same structure. The fixed platform of the mechanism relates to moving platform through these two joints of actuating branches. One joint of branches is the RPU kinematic branch chain, and the kinematic joints include a revolute (R) joint, a prismatic (P) joint, and a universal (U) joint. The other joint of branches is the SPR kinematic branch, and the kinematic joints include a spherical (S) joint, a prismatic (P) joint, and a revolute (R) joint. It is worth noting that the mechanism has 3 degrees of freedom, namely two rotational degrees of freedom and one translational degree of freedom. However, the number of actuating joints is four, and the number of actuating joints is greater than free degrees. So, the mechanism belongs to a redundant actuated mechanism.

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

Figure 1The five-axis linkage processing equipment.

Download

2.2 Establishment of mechanism coordinate system

With B as the center, a fixed coordinate system Bxyz is established on the mechanism base B1B2B3B4, in which the directions of the x axis and y axis are parallel to BB4 or BB1, and the direction of the z axis is determined according to the right-hand rule. The coordinate system of the moving platform A1A2A3A4 is established with the geometric center of A in the mechanism, and the directions of u axis and v axis on the coordinate system are parallel to AA4 or AA1. Similarly, the w axis is determined according to the right-hand rule (as shown in Fig. 2).

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

Figure 2Structure diagram of a parallel mechanism.

Download

2.3 Kinematics and dynamics analysis of 2RPU-2SPR parallel mechanism

The kinematics problem of a parallel mechanism is to construct the mapping relationship between the actuating component and the output component and to study the regular motion task performed by output components under the action of an actuating force or torque. A kinematics and dynamics analysis of the parallel mechanism is a critical step. It is not only a necessary condition for an integrated performance evaluation, scale synthesis and real-time control research but also a prerequisite for improving the motion accuracy of the mechanism. The introduction of the redundantly actuated structure of the redundantly actuated parallel mechanism means that it can effectively improve the stiffness of the mechanism and workspaces, avoid singularity, and keep the actuating force stable. However, this will make it difficult to analyze the kinematics and dynamics of the parallel mechanism.

As shown in Fig. 2, it can be obtained based on the closed vector method as follows:

(1) b i + d i s i = p + a i , a i = B R A A a i ,

where di represents the linear displacement of the actuator, and bi is the position vector of the joint hinge point in the fixed coordinate system. p=[xyz]T represent the position vector of moving platform coordinate system under the fixed coordinate system, and si is the unit vector of actuator. BRA represents rotation matrix of moving platform relative to the fixed coordinate system, ai represents position vector of Aai in fixed coordinate system, and Aai represents the position vector of the hinge point of moving platform in dynamic coordinate system, i.e., as follows:

(2) b 1 = 0 r b 0 T , b 2 = - r b 0 0 T , b 3 = 0 - r b 0 T , b 4 = r b 0 0 T A a 1 = 0 r a 0 T , A a 2 = - r a 0 0 T , A a 3 = 0 - r a 0 T , A a 4 = r a 0 0 T .

The homogeneous coordinate transformation matrix T can be obtained by three transformations, that is, revolute α around the x axis, moving h along the z axis, and then revolute β around the v axis, and BRA represents the rotation matrix. Then, the homogeneous coordinate transformation matrix can be expressed in Eq. (3) as follows:

(3) T = B R A p 0 1 = c β 0 s β 0 s α s β c α - s α c β - h s α - s β c α s α c α c β h c α 0 0 0 1 ,

where sin and cos are abbreviated as s and c, respectively.

As can be seen from Eq. (3), the position vector p=[xyz]T is a function of independent parameter variable X=[zαβ]T, and the coupled motion of its pose can be expressed in Eq. (4) as follows:

(4) x = 0 , y = - h sin α , z = h cos α , y = - z tan α .

The inverse position solution of Eq. (1) can be expressed in Eq. (5) as follows:

(5) d i = | | p + a i - b i | | , s i = ( p + a i - b i ) / d i ,

where ||•|| represents the Euclidean norm of actuating vector, and si is the unit vector of the actuator.

Equation (5) calculates arithmetic square root, which can be expressed in Eq. (6) as follows:

(6) d 1 = z tan α - r a c α + r b 2 + z + r a s α 2 d 2 = z + r a c α c β 2 + b - r a c α 2 + z tan α + r a s α s β 2 d 3 = z tan α + r a c α - r b 2 + z - r a s α 2 d 4 = z - r a c α c β 2 + b - r a c α 2 + z tan α - r a s α s β 2 .

Given the position of the moving platform, the actuating joint of 2RPU-2SPR over-constrained redundantly actuated parallel mechanism can be obtained from the inverse solution equation of pose, as shown in Eq. (6).

Then, we find the first derivative of time from Eq. (6) as follows:

(7) q ˙ = J 0 X ˙ ,

where q˙=[d˙1,d˙2,d˙3,d˙4]T and J0 are the Jacobian matrix of the parallel mechanism. For the derivation process of the Jacobian matrix, i.e., J0, has been explained by Zhang et al. (2020) in detail. The acceleration of the actuating joint can be obtained by deriving the time on both sides of Eq. (7), as follows:

(8) q ¨ = J ˙ 0 X ˙ + J 0 X ¨ ,

where J˙0 is the derivative of Jacobian matrix.

Due to the coupling relationship between branches of the parallel mechanism, its accurate dynamic model is extremely complex. Therefore, it is necessary to establish an accurate dynamic model as much as possible and build a mapping relationship between the output of the mechanism's moving platform and the force of actuating joint. According to the virtual work principle and D'Alembert's principle, the forces acting on mass center of each system component include gravity, inertial force, and inertial moment. When considering the whole system and ignoring influence of friction in components, the general expression of generalized coordinate dynamic model, based on workspace, can be expressed in Eq. (9) as follows:

(9) M ( X ) X ¨ + C ( X ˙ , X ) X ˙ + N ( X ) = J 0 T τ a ,

where M(X) is a positive definite inertia matrix of the parallel robot mechanism, C(X˙,X) is a centripetal and Coriolis force coefficient matrix, N(X) is the sum of gravity term matrix and external load, and τa is the actuating branch chain control force vector.

3 Establishment of mechanism performance evaluation index

3.1 Workspace evaluation index

Workspace refers to working area that can be reached by the end-effector on the moving platform, which is the direct embodiment of working performance to mechanism. It is also a crucial performance index to measure the comprehensive performance of a mechanism for a redundantly actuated parallel mechanism introduced in this paper. Wang et al. (2020) proposed a variable human–robot workspace based on the user and introduced two kinds of trajectory planning methods. Generally, the analysis methods for the workspace of parallel mechanism mainly include an analytical and numerical method. The method used in this paper is limit boundary search method, which divides the parallel mechanism into four motion branches and then determines the motion range of a single branch chain according to mechanism parameters and motion joint constraints. The envelope surface of the parallel mechanism motion space is obtained, and the workspace range of the mechanism is determined.

The size of a global workspace can be measured by the volume ratio of an effective workspace to search space, namely Eq. (10), as follows:

(10) WSI = N feasible d z d α d β N total d z d α d β .

In the formula, α, β is the attitude angle, Nfeasible indicates the reachable work area, and Ntotal represents entire work area.

WSI is a measurement index of a mechanism workspace. Its value indicates the effective workspace range of a mechanism. According to index definition, the value range of WSI should be between 0 and 1. In a search space, the value of 0 indicates that the end-effector has no pose points that can meet the motion conditions, and the value of 1 indicates that the end-effector can reach all points.

3.2 Motion/force transmission performance index

In order to optimize motion/force transfer of a parallel mechanism, Zhao et al. (2017) proposed the constant power transfer and motion/force transfer criteria of a parallel mechanism under a different position based on the screw theory, establishing a motion/force transfer characteristic index to easily unify the dimension based on the unit index in a calculation process. It is usually used as a main performance evaluation index of the redundantly actuated parallel mechanism, especially for the motion joint of a revolute and prismatic coupling mechanism.

The minimum value of the input and output transmission power is expressed as λ1 and λ2, respectively. The local motion/force transfer performance index can be expressed as the minimum value of the input and output indexes, i.e., χ. Considering that the configuration value changes with the change in position, the average value of the task space is used as global evaluation index to evaluate performance of mechanism which can be expressed in Eq. (11) as follows:

(11) χ = W t χ d W t W t d W t .

3.3 Stiffness evaluation index

The stiffness value of an over-constrained redundantly actuated parallel mechanism will directly affect the motion accuracy and stability when performing surface machining tasks. The greater the stiffness value of the mechanism, the stronger the ability of the mechanism to resist deformation. Sufficient stiffness can ensure that the mechanism maintains high motion accuracy and stability, even when it is disturbed by outside world. In order to optimize the dimension and structure parameters of parallel robotic manipulators, Wang et al. (2017) introduced a new type of stiffness performance evaluation index based on kineto-elastic statics analysis. The stiffness model of a manipulator was expressed as a set of two sets of matrix equations based on a unified mathematical formula, and the elasticity of the connection between the independent connecting rod and connecting rod was described in the form of constraints. Klimchik et al. (2018) proposed a unified and systematic method, which saved the calculation time of stiffness modeling.

The characteristic formula of stiffness matrix can be expressed in Eq. (11) as follows:

(12) τ = J 0 T K 0 J v Δ X , K = J 0 T K 0 J v ,

where, in the following:

(13)K0=Ka04×404×404×4Kc04×402×402×4Kτ,Ka=diagka1ka2ka3ka4,Kc=diagkr1kr2kr3kr4Kτ=kτl10+τ(e1×l10)kτl10+τ(e1×l10)01×201×2kτl30+τ(e1×l30)kτl30+τ(e1×l30)(14)Jv=JaJrcJτc,Jrc=Jc(1:4,1:6),Jτc=03×1Tl10T03×1T(e1×l10)T03×1Tl30T03×1T(e1×l30)TF(15)τi=e1×e2,e1=100T,e2=R100T.

kai is tension compression stiffness coefficient of a branch rod, and kri is bending stiffness coefficient.

In order to evaluate the stiffness distribution of a mechanism in the workspace, the maximum and minimum eigenvalues of K are usually calculated with Eq. (16) to measure stiffness performance of parallel mechanism, i.e., as follows:

(16) K 1 = max { eig ( K ) } .

In order to evaluate the overall stiffness performance of mechanism, the global stiffness index is defined in Eq. (17) as follows:

(17) k s = W K 1 d W W d W .

The larger the value of ks, the stronger the resistance to deformation of mechanism, and the better its stiffness performance, and the higher its motion accuracy when performing work tasks.

3.4 Dexterity index

Dexterity refers to how much posture the end-effector can obtain in a workspace. The better the dexterity of the mechanism, the more posture the end-effector can achieve in the workspace. However, its physical meaning is difficult to express clearly because of the inconsistency of dimensions for a parallel mechanism. Therefore, the characteristic length is used to normalize the Jacobian matrix so that the Jacobian matrix has no dimension. Wei et al. (2018) found the inverse of a robot forward-kinematics Jacobian matrix, which can be used in a singularity analysis, designed for isotropy and optimal control, based on a dexterity index.

Similarly, in order to accurately evaluate the motion flexibility of the mechanism in the workspace, the global condition number index is defined in Eq. (18) as follows:

(18) k g j = W k J d W W d W ,

where, in the following:

(19)kJ=1/cond(J0)=λmin(J0)/λmax(J0)(20)J0=J0(:,1)1LJ0(:,2:3)(21)J0=J0(:,1)J0(:,2:3)(22)L=n1Trace[J0(:,2:3)TJ0(:,2:3)]n2Trace[J0(:,1)TJ0(:,1)].

λmax and λmin represent the maximum and minimum singular values of matrix, and n1 and n2 represent the number of movement and rotation, where n1=1 and n2=2. “Trace” represents the sum of the main diagonal elements of the matrix. The range of kJ is [0,1], and its value changes with the change in the mechanism pose. Its value is closer to 1, which means that the flexibility of the mechanism is better. The closer the value is to 0, the worse the flexibility of the mechanism is, and the closer it is to a singular position, the harder mechanism is to control.

3.5 Energy efficiency index

The energy efficiency is one of the key indicators to evaluate the feasibility of a mechanism. It is a macro embodiment to measure the size parameters of a mechanism and the rationality of structural design. Zhang et al. (2019b) proposed a simple method to quickly solve the maximum virtual power coefficient, which greatly simplifies establishment process of motion force transmission performance evaluation index. At the level of dynamic characteristic analysis, a dimensionless mechanism performance evaluation index independent of coordinate system is established to evaluate the energy conversion efficiency of mechanism, which has important guiding significance for reasonably optimizing mechanical structure and improving dynamic performance of mechanism.

The output energy on the end-effector of the moving platform of a mechanism is regarded as the effective energy, and its proportion in all energy input by parallel mechanism is defined as the energy efficiency, which can be expressed in Eq. (23) as follows:

(23) η = E p / E all ,

where Ep is the energy of the moving platform, and Eall is the energy of a parallel mechanism system.

In order to evaluate the energy transfer efficiency of an organization more accurately, the global energy efficiency index is defined in Eq. (24) as follows:

(24) k ge = W η 2 d W W d W .

3.6 Inertia coupling index

One of the main factors affecting the dynamic response characteristics of parallel mechanisms is the inertial coupling index, and the coupling generally exists in redundant actuated parallel mechanisms. In fact, the coupling characteristics of inertial force/torque will affect the output accuracy of the mechanism. Li et al. (2019) established a dynamic model considering the joint friction effect and gave inertia coupling a characteristic evaluation index. Then, the variation law in the coupling evaluation index with the motion of the mechanism was studied in a workspace of the parallel mechanism. Therefore, in order to improve the control accuracy of mechanism and reduce the error caused by inertial coupling, it is of great significance to design a reasonable method to evaluate the inertial coupling characteristics of a parallel robot mechanism.

Similarly, in order to evaluate the inertial coupling characteristics of a mechanism more reasonably, the global inertial coupling characteristic index is defined in Eq. (25) as follows:

(25) k gc = W k c d W W d W ,

where, in the following:

(26)kc=|Mx12|+|Mx13|+|Mx21|+|Mx23|+|Mx31|+|Mx32||M11|+|M22|+|M33|(27)M(X)=(Mxij)3×3.

The meaning of the global inertial coupling characteristic index kgc is that the smaller kgc and the weaker coupling characteristic of mechanism leads to a higher control accuracy of the mechanism.

4 Research on multi-objective optimization

4.1 Design of initial parameters for multi-objective optimization

First, the initial structural parameters of a parallel mechanism are set on the premise of ensuring the machining of complex surfaces, and then the constraints are determined according to the engineering requirements and geometric relations. The crucial problems of an over-constrained redundantly actuated parallel mechanism are selecting reasonable structural parameters under known constraints for multi-objective optimization, which requires the comprehensive consideration of performance evaluation indexes to make the mechanism achieve a more balanced and excellent performance.

Therefore, according to the above design ideas and purposes, the mathematical model of a multi-objective optimization design of a redundantly actuated parallel mechanism can be described in Eq. (28) as follows:

(28) maximize / minimize  F ( X ) = k 1 ( X ) k 2 ( X ) k m ( X ) subject to  g i ( X ) 0 h i ( X ) = 0 ,

where X=x1x2xn, matrix X represents the parameter design variables, the element xi represents the structural parameters to be designed, and km(X) represents the objective function that is the performance evaluation index. gi(X) represents the inequality constraint equation, namely the design parameter range of the over-constrained parallel mechanism according to the engineering requirements. Simultaneously, hi(X) represents the equality constraint equation, namely the Jacobin matrix of the dynamic platform and actuating joint of the mechanism, the relationship between the branch chain and the dynamic platform, and the position inverse solution equation.

The above six global performance evaluation indexes are taken as the objective function of the 2RPU-2SPR over-constrained redundantly actuated parallel mechanism, namely workspace, stiffness, motion/force transmission efficiency, dexterity, energy efficiency, and the inertia coupling index. In order to facilitate the representation and description, the optimization objective can be expressed as ki(X).

This is where maxF(X)=k1(X)k2(X)k3(X)k4(X)k5(X) k6(X)], k1, k2, k3, k4, k5, k6 represent the global index of workspace, motion/force transmission efficiency, stiffness, dexterity, energy efficiency, and the inertia coupling respectively.

Therefore, the main optimization design variables of the 2RPU-2SPR over-constrained redundantly actuated parallel mechanism can be summarized in Eq. (29) as follows:

(29) X = r a r b h a h 1 h 2 R 1 R 2 ,

where ra, rb, ha, h1, h2, R1, and R2 represents the radius of the moving platform, fixed platform radius, thickness of moving platform, length of actuating lower branch chain, length of actuating upper branch chain, outer circle radius of actuating lower branch chain, and outer circle radius of actuating upper branch chain, respectively.

4.2 Multi-objective particle swarm optimization algorithm and orthogonal experimental design method

The multi-objective particle swarm optimization algorithm is a biological evolutionary computer technology. Its basic principle is to randomly select the population within the range of design parameter variables, and then the initialized population corrects its own parameters through learning. The correction process depends on the previous population, and the population constantly updates to search for local optimization. There is no crossover and mutation, and the setting parameters are simple compared with the genetic algorithm. It does not need a frequent adjustment and is easy to implement. However, it is not conducive to find the global best if the initial population deviates from the Pareto solution set and the search process will be long. In other words, the selection of the initial population will have a great impact on the search ability of the particle swarm optimization algorithm and will directly affect the speed and efficiency of the multi-objective optimization.

Trial design is a crucial step in mechanism parameter design. An appropriate test method can not only reduce the cost and experiment time but also help to obtain more accurate and actual orthogonal experiment results. There are many experimental design methods, including the factor experiment, response surface model, and Taguchi experimental design method. The test method used in this paper is the orthogonal experiment method, which is an efficient method for finding the objective solution of multiple design parameters. It can analyze various changes in input and output parameters under the conditions of a short test cycle and little test data.

4.3 Multi-objective particle swarm optimization based on orthogonal experimental design

Through the above analysis, a multi-objective optimization of the structural parameters of a 2RPU-2SPR over-constrained redundantly actuated parallel mechanism needs to comprehensively measure seven design variables and six objective functions, which makes it difficult to establish a dynamic and kinematic model. Therefore, this paper aims to find a simple and effective multi-objective optimization design method.

An improved multi-objective particle swarm optimization configuration algorithm, based on orthogonal experimental response surface fitting design, is proposed. This method combines the advantages of an orthogonal experimental design method and multi-objective particle swarm optimization algorithm to effectively analyze the influence of the initial population on a multi-objective optimization algorithm. In the process of an optimization calculation, a fuzzy experimental value is generated in advance through the orthogonal experimental design. The experimental value is taken as the initial population in a feasible region, and then the multi-objective particle swarm optimization algorithm is used to search uniformly and carefully in the whole feasible region. The basic steps of the multi-objective particle swarm optimization algorithm, based on an orthogonal experimental design, are as follows:

  1. Establishing the multi-objective optimization design mathematical model of a redundantly actuated 2RPU-2SPR over-constrained redundantly actuated parallel mechanism, involving the design variables, constraints, and objective functions.

  2. Determining the factor of the design parameters and selecting an appropriate orthogonal table to obtain the data of the experiment.

  3. Analyzing the correlation between parameter variables and the objective function and calculating the fitting accuracy of the output results.

  4. Allowing, on this basis, for the population to search in the whole feasible region and, finally, obtaining the optimal Pareto solution and Pareto frontier with dominant indexes.

Table 1Experiment design objective function value.

Download Print Version | Download XLSX

5 Simulation analysis of multi-objective optimization

5.1 Initial optimal Pareto solution and Pareto frontier

In the multi-objective particle swarm optimization design, the population size is 10, maximum number of iterations is 50, inertia weight is 0.9, global increment is 0.9, particle increment is 0.9, maximum speed is 0.1, penalty value for the running failure is 1.0 × 1030, and target value for the running failure is 1.0 × 1030. The parameter range of the design variable factors is shown in Table 1.

Table 2Parameter variable design matrix in the experimental design.

Download Print Version | Download XLSX

According to the parameter range of the design variable factors in Table 1, a 50 × 7 experiment design matrix can be generated, and some data are shown in Table 2.

Table 3Experiment design objective function value.

Download Print Version | Download XLSX

The data points of the experiment in the design matrix will be programmed and calculated to obtain the corresponding objective function values. Some results are shown in Table 3.

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

Figure 3Parametric design variable combination matching.

Download

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

Figure 4Multi-objective optimization design Pareto front target response value.

Download

After the simulation analysis of the experimental design, the ideal Pareto solution and Pareto front of the mechanism can be obtained. The value process of the design parameters h1, h2, ha, r1, r2, ra, and rb is shown in Fig. 3. The blue line represents the best design variable matching combination in a feasible region, and the red dot represents the ideal initial optimal Pareto solution, with values of 0.5061, 0.698, 0.01, 0.02424, 0.1473, 0.4, and 0.469, respectively. The response values of the objective functions k1, k2, k3, k4, k5, and k6 are shown in Fig. 4. The blue line represents function response value corresponding to the combination of design variables in feasible area, and the red dot represents Pareto front, with values of 0.260417, 0.682066, 7.63×108, 0.423911, 0.862877, and 0.30866, respectively.

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

Figure 5Main effect diagram of the different objective functions.

Download

5.2 Mapping relationship between design variables and multi-performance indicators

The corresponding target performance change trend can be obtained according to different levels of design variables. It can be seen from Fig. 5a that the workspace performance k of a parallel mechanism is affected by the thickness ha of moving platform ra and fixed platform rb. According to curve change trend, the larger the values of design variables ra and rb, the lower the average value of performance evaluation index k1. It can be seen from Fig. 5b that the rod length h2 and thickness ha of the moving platform have an impact on the performance evaluation index k2, and the impact needs to be considered in sections. The first half shows that the rod length h2 and the thickness ha of the moving platform have a positive correlation with the flexibility of the mechanism, and the latter part of the conclusions is just the opposite. It also shows that two design variables have a nonlinear relationship with evaluation index k2. Figure 5c shows that, the larger the design variable ra is, the higher the mean value of performance evaluation index k3, while the larger the design variable rb, the lower the mean value of performance evaluation index k3. As shown in Fig. 5d, the performance evaluation index k4 is affected by the design variables ra and rb, which is similar to the performance evaluation index k3. The extreme points of the design variable ra in Fig. 5e and f on performance evaluation indexes k5 and k6 are not uniformly monotonous. In short, the existence of a main effect relationship depends on the trend of the change curve. If the slope of the curve is larger, then it indicates that the design variable has a greater impact on the target performance.

5.3 Correlation analysis of input parameter variables to output targets

Correlation analysis is a linear analysis method based on Pearson correlation coefficient, Spearman correlation coefficient, and Kendall correlation coefficient. For the Pearson correlation coefficient, a calculation formula can be expressed in Eq. (30) as follows:

(30)rXY=(X-X)(Y-Y)(X-X)2(Y-Y)2=xyx2y2=SXYSXSY(31)x=X-X,y=(Y-Y)(32)SXY=xyn-1,SX=x2n-1,SY=y2n-1,

where SXY represents the total variance of the sample, SX represents the sample standard deviation of X, and SY represents the sample standard deviation of Y.

https://ms.copernicus.org/articles/13/123/2022/ms-13-123-2022-f06

Figure 6Correlation matrix of the design variables and objective functions.

Download

The characterization of the correlation can be measured by a positive correlation and negative correlation. If 0<r<1, it means that the dependent variable is positively correlated with independent variable. At this time, the dependent variable will increase with the increase in the independent variable; if -1<r<0, it means that the dependent variable is negatively correlated with the independent variable. At this time, the dependent variable will decrease with the increase in the independent variable. The closer the absolute value of r is to 1, the higher correlation between the independent variables and dependent variables and the stronger the dependence between them will be. On the contrary, the closer the absolute value of r is to 0, the weaker the correlation between independent variable and dependent variable will be. The numerical relationships of the correlation matrix between the design variables and objective function of the 2RPU-2SPR over-constrained parallel mechanism are shown in Fig. 6.

The prediction algorithm of the regression model usually adopts the mean square error, root mean square error, mean absolute error, R2, adjusted R2, etc. In this paper, the determination coefficient method is used to predict the regression model, and its mathematical can be expressed in Eq. (33) as follows:

(33) R 2 = 1 - i ( y ^ i - y i ) 2 i ( y i - y i ) 2 ,

where i(y^i-yi)2 represents the sum of the square differences between real value and predicted value, and i(yi-yi)2 represents the sum of the square differences between real value and mean value.

The numerical range of the fitting accuracy is (0,1). The quality of the regression model can be judged according to the value of the fitting accuracy. When the fitting accuracy of the experiment sample is higher, the value is closer to 1, and when the fitting accuracy is lower, the value is closer to 0. The fitting accuracy between the target functions in experiment design is shown in Fig. 7.

https://ms.copernicus.org/articles/13/123/2022/ms-13-123-2022-f07

Figure 7Fitting accuracy between the objective functions.

Download

Table 4Optimal Pareto solution set for a multi-objective optimization.

Download Print Version | Download XLSX

Because the data samples of the experiment are relatively discrete, the full factor experiment design is not adopted, and there is coupling and competition among the objective functions, making the distribution law not obvious, so it is difficult to achieve a high degree of response surface linear fitting. Through multi-objective particle swarm optimization, 48 sets of Pareto non-inferior solution sets can be obtained because of the orthogonal experimental design, and the response value of the corresponding objective function is the Pareto front, as shown in Table 4.

Table 5Design values of the main parameters of 2RPU-2SPR over a constrained parallel mechanism.

Download Print Version | Download XLSX

Table 6Comparison of objective function values after an improved collaborative multi-objective optimization algorithm.

Download Print Version | Download XLSX

Table 7Comparison of three traditional optimization algorithms.

Download Print Version | Download XLSX

5.4 Comparative analysis of optimization results

After the optimization, there are 48 groups of optimal Pareto solutions. But these objective functions are basically contradictory and conflicting. The increase in one performance index will inevitably damage the other. The Pareto solution cannot continue to improve the performance of the mechanism without weakening other performance indexes. Finally, the optimization system selects the design matching value of the 14th group of parameters as the optimal Pareto solution, according to the minimum distance method. The values of the design parameters before and after optimization are shown in Table 5, and the response values of performance evaluation indexes before and after optimization are shown in Table 6.

It can be seen from Table 6 that the workspace index k1 has been reduced by 27.27 % after optimization, while, for the global motion/force transmission efficiency k2, global stiffness index k3, global dexterity index k4, and global energy efficiency index k5, it can be seen from the numerical value that the improvement rates after optimization are 21.80 %, 16.69 %, 72.98 %, and 38.83 %, respectively. The smaller the value of the global inertia coupling index k6, the better. Its performance index is improved by 14.59 % after optimization, indicating that the improvement of other performance evaluation indexes must be at expense of the workspace, and the optimization design result is feasible.

Table 7 compares the three intelligent optimization algorithms. The comparison results show that the traditional orthogonal experimental design is not very ideal. The global stiffness, global flexibility and global energy transfer efficiency index have decreased. The particle swarm optimization algorithm shows better results compared with the orthogonal experimental design. However, the stiffness performance index is more of a concern for mechanism designers for complex surface machining tasks. So, the results obtained by particle swarm optimization algorithm are also not very ideal. Using the particle swarm optimization algorithm based on the orthogonal experimental design proposed in this paper, the stiffness performance of the mechanism is greatly improved, and there are some improvements in other performance indexes, which are more in line with the application needs of practical machining.

The optimization comparison results show that multi-objective optimization design problem is complex, and the relationship between the performance indexes of the mechanism are influenced and restricted by each other. The improvement in some performance indexes must be based on sacrifice of other indexes. The parameter design should be considered selectively and comprehensively, according to needs of actual processing tasks. The over-constrained redundantly actuated parallel mechanism 2RPU-2SPR is mainly for specific complex surface machining tasks. Figure 8 shows the optimal design combination of parameters obtained by multi-objective particle swarm optimization, which can be selected and decided on for designers.

https://ms.copernicus.org/articles/13/123/2022/ms-13-123-2022-f08

Figure 8Pareto solutions and the optimal solution.

Download

As shown in Fig. 8, the competitive relationship between the performance evaluation indexes k1 and k6 is relatively obvious. The increase in the k1 value will, correspondingly, lead to the increase in the k6 value. However, the performance evaluation indexes of k2 and k6 maintain a negative growth relationship in local areas. Therefore, we should consider the value of each objective function and compromise on the parameter design variables in the multi-objective parameter design, so as to obtain the best mechanism performance and the best design parameter combination as far as possible.

6 Conclusions

In order to complete the machining task of a complex curved surface accurately and efficiently, this paper proposes an over-constrained redundantly actuated parallel mechanism and introduces and analyses six performance evaluation indexes. However, because of difficulties in establishing kinematic and dynamic models, it is laborious to balance the optimization of multi-performance indexes. Therefore, a multi-objective particle swarm optimization configuration algorithm, based on an orthogonal experimental design, is proposed. The orthogonal experimental design is used to obtain the initial population through polynomial function fitting, and then the multi-objective particle swarm optimization algorithm is used to search uniformly and carefully in the whole feasible region space, so as to obtain optimal scale parameters. It is of great significance to improve the potential working space, kinematics, stiffness, and dynamic performance of the parallel mechanism and to meet the requirements of surface processing tasks. Simultaneously, the research results of this paper have an important theoretical basis and engineering application value to guide the development of hybrid configuration processing equipment.

The following conclusions are drawn from the experimental results:

  1. The kinematics and dynamics analysis are the basis of an integrated performance evaluation, scale synthesis, and real-time control researches for an over-constrained redundantly actuated parallel mechanism. This paper analyses the pose inverse solution of the 2RPU-2SPR parallel mechanism, based on the closed vector method, and briefly introduces the kinematics and dynamics of the mechanism.

  2. For the performance requirements of special machining, this paper proposes a performance evaluation system for the 2RPU-2SPR over-constrained redundantly actuated parallel mechanism, which integrates the global index of workspace, motion/force transmission efficiency, stiffness, dexterity, energy efficiency, and inertia coupling. It is conducive to the subsequent comprehensive optimization design of the 2RPU-2SPR over-constrained redundantly actuated parallel mechanism, based on an orthogonal experimental design algorithm and multi-objective particle swarm optimization algorithm.

  3. Compared with the two traditional intelligent optimization algorithms, the results show that particle swarm optimization configuration algorithm, based on an improved orthogonal experiment design, has a better optimization effect in which the workspace, motion/force transmission efficiency, stiffness, dexterity, energy efficiency, and inertia coupling index of optimized 2RPU-2SPR over-constrained redundantly actuated parallel mechanism are improved to a certain extent.

  4. An over-constrained redundantly actuated parallel mechanism is a special kind of parallel mechanism which has the characteristics of multi-variability, strong coupling, and complex nonlinearity. The evaluation indexes of a multi-objective optimization design are often contradictory, and the designer needs to compromise and select an appropriate parameter design value.

Data availability

All the data used in this article can be made available upon reasonable request. Please contact the corresponding author (hqzhang@sdust.edu.cn).

Author contributions

HZ proposed the novel parallel robot and developed the theoretical analysis. JT and QG wrote the paper under the guidance of HZ. GC and YY verified the established the data results and supervised the whole research work.

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 The 30th International Conference on Robotics in Alpe–Adria–Danube Region, RAAD 2021, Futuroscope-Poitiers, France, 21–23 June 2021.

Acknowledgements

The authors would like to acknowledge the financial support from the Fundamental Research Funds for Shanghai Collaborative Innovation Center of Intelligent Manufacturing Robot Technology for Large Components (grant no. ZXP20211101) and the Key Laboratory of Vehicle Advanced Manufacturing, Measuring and Control Technology (Beijing Jiaotong University), Ministry of Education, China.

Financial support

This research has been supported by the Fundamental Research Funds for Shanghai Collaborative Innovation Center of Intelligent Manufacturing Robot Technology for Large Components (grant no. ZXP20211101) and Key Laboratory of Vehicle Advanced Manufacturing, Measuring and Control Technology (Beijing Jiaotong University), Ministry of Education, China.

Review statement

This paper was edited by Juan Sandoval and reviewed by four anonymous referees.

References

Asker, A., Assal, and S. F. M.: Systematic approach for designing a multi-function sit-to-stand mobility assistive device based on performance optimization, Adv. Robotics., 33, 90–105, https://doi.org/10.1080/01691864.2019.1573152, 2019. 

Chen, X. L. and Sun, X. Y.: Dexterity Analysis of a 4-UPS-RPS Parallel Mechanism, Int. J. Adv. Robot. Syst., 

Enferadi, J. and Nikrooz, R.: The Performance Indices Optimization of a Symmetrical Fully Spherical Parallel Mechanism for Dimensional Synthesis, J. Intell. Robot. Syst., 90, 1–17, https://doi.org/10.1007/s10846-017-0675-6, 2017. 9, 591–593, https://doi.org/10.4028/www.scientific.net/AMR.591-593.754, 2012. 

Houssein, E. H., Mahdy, M. A., and Shebl, D.: An efficient slime mould algorithm for solving multi-objective optimization problems, Expert. Syst. Appl., 187, 1–17, https://doi.org/10.1016/j.eswa.2021.115870, 2022. 

Klimchik, A., Pashkevich, A., and Chablat, D.: Fundamentals of manipulator stiffness modeling using matrix structural analysis, Mech. Mach. Theory, 133, 365–394, https://doi.org/10.1016/j.mechmachtheory.2018.11.023, 2018. 

Li, Y. B., Zheng, H., Xu, T. T., Peng, S., Wang, Z. S., and Qin, S. Y.: Dynamic Modeling with Joint Friction and Research on the Inertia Coupling Property of a 5-PSS/UPU Parallel Manipulator, Int. J. Mach. Tool. Manu., 55, 43–52, https://doi.org/10.3901/JME.2019.03.043, 2019. 

Liu, R. C., Li, J., and Jing, F.: A dynamic multiple populations particle swarm optimization algorithm based on decomposition and prediction, Appl. Soft. Comput., 73, 434–459, https://doi.org/10.1016/j.asoc.2018.08.015, 2018. 

Reveles, D., Pamanes, A. G., and Wenger, P.: Trajectory planning of kinematically redundant parallel manipulators by using multiple working modes, Mech. Mach. Theory, 98, 216–230, https://doi.org/10.1016/j.mechmachtheory.2015.09.011, 2016. 

Russo, M., Herrero, S., Altuzarra, O., and Ceccarelli, M.: Kinematic analysis and multi-objective optimization of a 3-UPR parallel mechanism for a robotic leg, Mech. Mach. Theory, 120, 192–202, https://doi.org/10.1016/j.mechmachtheory.2017.10.004, 2018. 

Tian, W., Shen, Z., Lv, D., and Yin, F.: A Systematic Approach for Accuracy Design of Lower-Mobility Parallel Mechanism, Robotica, 38, 2173–2188, https://doi.org/10.1017/S0263574720000028, 2020. 

Wang, H., Zhang, L., Chen, G., and Huang, S.: Parameter optimization of heavy-load parallel manipulator by introducing stiffness distribution evaluation index, Mech. Mach. Theory, 108, 244–259, https://doi.org/10.1016/j.mechmachtheory.2016.10.011, 2017. 

Wang, H. B., Lin, M., Jin, Z., and Yan, H.: A 4-DOF Workspace Lower Limb Rehabilitation Robot: Mechanism Design, Human Joint Analysis and Trajectory Planning, Appl. Sci., 10, 42–45, https://doi.org/10.3390/app10134542, 2020. 

Wei, L. and Angeles, J.: The design of a 3–CPS parallel robot for maximum dexterity, Mech. Mach. Theory, 122, 279–291, https://doi.org/10.1016/j.mechmachtheory.2018.01.003, 2018. 

Wu, X. Y. and Bai, S. P.: Architectural singularities of parallel mechanisms with prismatic joints due to special designs of platform shapes, Mech. Sci., 10, 449–464, https://doi.org/10.5194/ms-10-449-2019, 2019. 

Wu, Y., Shi, L. L., and Liu, X. X.: A new dynamic strategy for dynamic multi-objective optimization, Inform. Sciences, 529, 116–131, https://doi.org/10.1016/j.ins.2020.04.011, 2020. 

Xiang, J., Han, X. H., and Duan, F.: A novel hybrid system for feature selection based on an improved gravitational search algorithm and k-NN method, Appl. Soft. Comput., 31, 293–307, https://doi.org/10.1016/j.asoc.2015.01.043, 2015. 

Yang, Y., Tang, Y., Chen, H., Peng, Y., and Pu, H.: Mechanism design and parameter optimization of a new asymmetric translational parallel manipulator, Mech. Sci., 10, 255–272, https://doi.org/10.5194/ms-10-255-2019, 2019. 

Zeng, N., Song, D., Li, H., and You, Y.: A competitive mechanism Integrated multi-objective whale optimization algorithm with differential evolution, Neurocomputing, 432, 170–182, https://doi.org/10.1016/j.neucom.2020.12.065, 2021. 

Zhang, H., Fang, H., Zhang, D., Luo, X., and Zhao, F.: Kinematics and Dynamics Simulation Analysis of a 3-DOF Parallel Mechanism for Application in Hybrid Machine, in: Recent Advances in Mechanisms, Transmissions and Applications, edited by: Wang, D., Petuya, V., Chen, Y., and Yu, S., MeTrApp 2019, Mechanisms and Machine Science, Vol. 79, Springer, Singapore, https://doi.org/10.1007/978-981-15-0142-5_25, 2020. 

Zhang, H. Q., Fang, H., Zhang, D., Luo, X., and Zou, Q.: Forward kinematics and workspace determination of a novel redundantly actuated parallel manipulator, Int. J. Aerospace Eng., 2019, 1–14, https://doi.org/10.1155/2019/4769174, 2019a. 

Zhang, H. Q., Fang, H. R., and Jiang, B. S.: Motion-force transmissibility characteristic analysis of a redundantly actuated and over constrained parallel machine, Int. J. Autom. Comput., 16, 150-162, https://doi.org/10.1007/s11633-018-1156-5, 2019b. 

Zhang, H. Q., Fang, H. R., and Jiang, B. S.: Dynamic performance evaluation of a redundantly actuated and over-constrained parallel manipulator, Int. J. Autom. Comput., 16, 274–285, https://doi.org/10.1007/s11633-018-1147-6, 2019c. 

Zhang, X. Y. and Zheng, X. T.: A competitive mechanism based multi-objective particle swarm optimizer with fast convergence, Inform. Sciences, 427, 63–76, https://doi.org/10.1016/j.ins.2017.10.037, 2017. 

Zhang, Z. and Meng, G.: Design and analysis of a six degrees of freedom serial–parallel robotic mechanism with multi-degree of freedom legs, Int. J. Adv. Robot. Syst., 15, 1–14, https://doi.org/10.1177/1729881418812643, 2018. 

Zhao, Y. Z., Wang, J. L., Cao, Y. C., Liang, B. W., and Zhao, T. S.: Constant motion/force transmission analysis and synthesis of a class of translational parallel mechanisms, Mech. Mach. Theory, 108, 57–74, https://doi.org/10.1016/j.mechmachtheory.2016.10.008, 2017. 

Zou, Q., Zhang, D., Zhang, S., and Luo, X.: Kinematic and dynamic analysis of a 3-DOF parallel mechanism, Int. J. Mech. Mater. Des., 17, 587–599, https://doi.org/10.1007/s10999-021-09548-8, 2021. 

Download
Short summary
Compared with the two traditional intelligent optimization algorithms, the results show that the particle swarm optimization configuration algorithm, based on an improved orthogonal experiment design, has a better optimization effect, and the results are improved to a certain extent. The evaluation indexes of a multi-objective optimization design are often contradictory. The designer needs to compromise and select an appropriate parameter design value.