Modified Pso Method for Robust Control of 3rps Parallel Manipulators

Copyright q 2010 T.-S. Zhan and C.-C. Kao. 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. We propose an effective method to design a modified particle swarm optimization MPSO singularity control method for a fully parallel robot manipulator. By adopting MPSO to obtain simple and effective estimated damping values, the result automatically adjusts the damping value around a singular point and greatly improves the accuracy of system responses. This method works by damping accelerations of the end effector, so that accelerations along the degenerated directions are zero at a singular point. These velocities, however, may not be zero in some situations, in which case, fluctuations will be encountered around a singular point. To overcome this drawback, we propose a control scheme that uses both damped acceleration and damped velocity, called the hybrid damped resolved-acceleration control HDRAC scheme. The MPSO optimization method can immediately provide optimal damping factors when used in an online application. Our proposed approach offers such superior features as easy implementation, stable convergence characteristics, and good computational efficiency. The main advantage of the HDRAC with MPSO in the 3RPS parallel manipulator control system is that it is not necessary for the system to plan its path for avoiding the singular point; thus, the workspace can be improved. Illustrative examples are provided to show the effectiveness of this HDRAC in practical applications, and experimental results verifying the utility of the proposed control scheme are presented.


Introduction
Parallel robotics offers the advantages of heavy load capacity, high rigidity, and low inertia.The triaxial RPS parallel manipulator, which is equipped with revolute joints on the base platform and spherical joints on the moving platform, has been widely applied to electronic gaming machines, virtual reality movies, simulations of aircraft and vehicles, and so on.In their research on 3RPS parallel robotics, Lee and Shah 1 used loop closure equations to deduce kinematic equations for the relationship between link lengths and the position of the moving platform.In a separate study, they also deduced the dynamic equations that included the drive forces and reacting forces of 3RPS parallel robotics in the case of inverse dynamics 2, 3 .Following their work, Fang and Huang 4 simplified the mathematic models and established the relationship among three input-output variables.Li and Xu 5, 6 deal with the kinematic characteristics in terms of workspace and dexterity vary with differences in the arrangement of actuators are investigated in detail.Generally, there are three types of singularity of parallel mechanisms, that is, inverse, direct, and combined singularities.In this study, we deal with the direct kinematic singularities.Therefore, the author divided the procedure of finding "direct kinematic singularities" of 3RPS platform mechanism.In horizontal position of the moving platform, to begin with, he assumed the direction of the moving platform and determined its horizontal position at this moment, that is, its projection on the fixed base.Then, since the direction of the moving platform has been assumed, also, a moving platform remaining in this direction can only translate vertically, it is possible to locate the limiting position while this direction has always been kept; this is the boundary of workspace, which makes the determinant of Jacobian matrix equal to zero and thus is the "direct kinematic singularities".
There may exist singular configurations emerging from the motions of mechanisms, and resolved-rate control and resolved-acceleration control 7-9 are the two general approaches used to control robots.In the neighborhood of singularities, very high speeds are required to produce even small changes in the position or orientation of the end effector.In this paper, singular value decomposition SVD is used in order to show that fluctuations around the singular point are due to unnecessary nonzero end-effector velocities along the degenerated directions when the manipulator is at a singularity.This paper proposes a hybrid damped resolved-acceleration control HDRAC scheme, which considers both damped acceleration and velocity as a way to remove these unnecessary velocities.
Particle swarm optimization PSO is a stochastic population-based optimization approach and was first published by Kennedy and Eberhart in 1995 10, 11 .PSO has been shown to be an efficient, roust, and simple optimization algorithm.Most studies of the PSO are empirical with only a few theoretical analyses that concentrate on understanding particle trajectories.A large amount of research has been proposed to improve the performance of PSO.From these studies, much effort has been invested to obtain a better understanding of its convergence properties.These studies concentrated mostly on a better understanding of the basic PSO control parameters, namely, the acceleration coefficients, inertia weight, velocity clamping, and swarm size 12-14 .In this paper, a modified PSO algorithm is proposed to improve the searching ability and prevent from being trapped in a local optimum.The main difference of the MPSO from the PSO is its fitness function which considers the "distance" to avoid converging to a local optimum.From these empirical studies, it can be concluded that the MPSO is sensitive to control parameter choices, specifically the inertia weight, acceleration coefficients, and velocity clamping 15 .However, wrong initialization of these parameters may lead to divergent or cyclic behavior.
In this study, singular value decomposition SVD is used to show that the abovementioned fluctuation is due to unnecessary nonzero end-effector velocities along the degenerated directions when the manipulator is at a singularity.The hybrid damped resolved-acceleration control scheme HDRAC , which considers both damped acceleration and velocity, is proposed to remove these unnecessary velocities.Numerical simulations and the results of experiments are compared in order to demonstrate the applicability and effectiveness of the proposed algorithms at a singularity for a 3RPS parallel manipulator.

Singularity Analysis of the 3RPS Parallel Manipulator
As shown in Figure 1, one end of the 3RPS parallel manipulator is a moving platform and the other end is a fixed base.The moving platform is connected to the fixed base by three limbs via three spherical joints B 1 , B 2 , and B 3 and three revolute joints A 1 , A 2 , and A 3 .
Each limb consists of two links connected in a series by a prismatic joint.Figure 2 shows that the origin P of the coordinate system x y z T is located at the mass center of the base moving platform.The Cartesian coordinate vector X Y Z α β γ T denotes the fixed coordinates, where X Y Z T and α β γ T denote, respectively, the vectors of position and orientation.The fixed and moving platforms of the 3 RPS are circles with radius a and b, respectively.The distance between the three spherical joints on the moving platform is w.

Dynamics of the 3PRS Parallel Manipulator
Let us focus on the mass center O i of the ith link, as shown in Figure 2.There are three components of the weight force of the moving platform that act on the ith link as i F pix , i F piy , and i F piz .The weight force of the ith link is W i , the acting force of the separately excited DC servomotor on the revolution parts is i F 1motor , and the acting force on each vertical leg is i F 2motor .The 3RPS parallel manipulator dynamic motion of the links and moving platform is described by using Newton-Euler equations, which are formulated by separately considering the free-body diagrams of the links and the moving platform 16-19 .Through reduction and incorporation, we can calculate the acting force on each vertical leg by using the following force equation for the moving vertical leg 19 : where W li is the weight of the ith link, i F pix D −1 Q, and matrices D and Q are given in Appendix A.
With regard to the revolution parts, we have the equation for the first motor as follows: where K t1 is the torque coefficient, K b1 is the back emf of the motor, l a1 is the armature inductance, r a1 is the armature resistance, V T is the terminal voltage, ω is the rotor speed, and i a1 is the armature current.

Deduction of Direct Kinetic Singularities
In this section, we adopt a method 16-19 to formulate the direct singular position of the 3RPS parallel manipulator.The distance between any two points of B 1 , B 2 , and B 3 is always w, which gives us the following relationships:

2.4
By differentiating 2.4 with respect to time, and writing them in matrix form, we get

Particle Swarm Optimization
Birds particles flocking optimizes a certain objective function in a PSO system.Each agent knows its best value so far pbest and its position.This information is an analogy of personal experiences of each agent.Moreover, each agent knows the best value so far in the group gbest among pbests.This information is an analogy of knowledge of how the other agents around them have performed.The PSO concept 10, 11 consists of changing the velocity of each particle toward its pbest and gbest locations.In the PSO, each particle moves to a new position according to new velocity and the previous positions of the particle.This is compared with the best position generated by previous particles in the fitness function, and the best one is kept; so each particle accelerates in the direction of not only the local best solution but also the global best position.If a particle discovers a new probable solution, other particles will move closer to it to explore the region more completely in the process.
In general, there are three attributes, current position x j , current velocity v j , and past best position pbest j , for particles in the search space to present their features.Each particle in the swarm is iteratively updated according to the aforementioned attributes.For example, 10-14 , the jth particle is represented as x j x j,1 , x j,2 , . . ., x j,g in the gdimensional space.The best previous position of the jth particle is recorded and represented as pbest j pbest j,1 , pbest j,2 , . . ., pbest j,g .The index of best particle among all particles in the group is represented by the gbest g .The rate of the position change velocity for particle j is represented as v j v j,1 , v j,2 , . . ., v j,g .The modified velocity and position of each particle can be calculated using the current velocity and distance from pbest j,g to gbest j,g as shown in the following formulas 14 : where n is the number of particles in a group; m is the number of members in a particle; t is the pointer of iterations generations ; v t j,g is the velocity of the particle j at iteration t, V min g ≤ v t j,g ≤ V max g ; w is the inertia weighting factor; c 1 , c 2 are the acceleration constants; Rand , Rand * are random numbers between 0 and 1; x t j,g is the current position of particle j at iteration t; pbest j is the pbest of particle j; gbest g is the gbest of the group g.
In the above procedures, the parameter V max g determines the resolution or fitness, with which regions are searched between the present position and the target position.If V max g is too high, particles might fly past good solutions.If V max g is too low, particles may not explore sufficiently beyond local solutions.
The constants c 1 and c 2 represent the weighting of the stochastic acceleration terms that pull each particle toward pbest and gbest positions.Low values allow particles to roam far from the target regions before being tugged back.On the other hand, high values result in an abrupt movement toward or past target regions.
Suitable selection of inertia weighting factor w provides a balance between global and local explorations, thus requiring less iteration on average to find a sufficiently optimal solution.As originally developed, w often decreases linearly from about 0.9 to 0.4 during a run.In general, the inertia weighting factor w is set according to the following equation 11, 12 : where iter max is the maximum number of iterations generations , and iter is the current number of iterations.

Modified Particle Swarm Optimization
The main point of the MPSO differs from the PSO to consider the "distance" in its fitness function to avoid converging to a local optimum.Assign a rank i.e., the number place 1, 2, 3, . . ., etc. RE k to the calculated error of each new individual, v k , k 1, . . ., PS, PS is the population size.A combined population with 2 × PS individuals is formed.Unlike previously developed statistic methods, the concept of "distance" is added to the fitness function to prevent from being trapped in a local minimum.The fitness score of the kth individual is modified by 20, 21 where ρ is an adaptive decay scale, ρ max is set as 0.7, and ρ min is set as 0.005 in this paper.RD k is the rank of D k assigned to the kth individual, where D k is the distance from the individual to the current best solution vector and is given by where v k is the vector of the kth individual in the combined population, and v best is the current best solution vector.An adaptive scheme is defined as 20

Mathematical Problems in Engineering
where ρ Δ is the step size; F min is the minimum value of fitness functions; R is the regulating scale and is set as 1.25 in this paper, and g max is the maximum allowable number of iterations.Individuals will be ranked ascending according to their fitness scores by a sorting algorithm.The PS individuals are transcribed along with their fitness for the next generation.If the new population does not include the current best solution, the best solution must be replaced with the last individual in the new population.In addition, a gradually decreased decay scale can satisfy a successive statistic searching process by first using the diversification bigger ρ to explore more regions and then the intensification smaller ρ to exploit the neighborhood of an elite solution.The current best solution point A for a minimum fitness problem as shown in the Figure 3 may not reach the global optimum 20, 21 , and there are three electable solutions.Generally, solutions with slightly better fitness point C or B prevailed, so the solution was trapped into the valley prematurely.The more attractive solution point G is relatively far away from point A, but it nears the global optimal.To prevent prematurity, point G with slightly worse fitness than C, it needs a higher rank to be selected.That is, a higher RD k is awarded to a longer D k .

Stopping Criteria
Stopping criteria is given in the following order: 1 maximum allowable number of iterations reached.
2 number of iterations reached without improving the current best solution.
Figure 4 shows the flow chart of the proposed algorithm.

Resolved-Acceleration Control
The relationship between the velocities of the mass center P of the moving platform and the joint velocity q for parallel manipulators can be represented as where ṙ and ω are the linear and angular velocities, respectively.The Cartesian coordinate vector q X, Y, Z, α, β, γ T with six variables is chosen to describe the position and orientation of the moving base, that is, q is the velocity and J is the Jacobian matrix.Differentiating 4.1 gives a ≡ r α J q J q, 4.2 where r and α are the linear and angular accelerations, respectively.Now, by expanding 2.1 -2.3 , we can obtain the dynamic equation of the manipulator in the matrix form of τ M q q N q, q G q M q q f q, q , 4.3 where M q is the inertia matrix of the manipulator, N q, q contains the Coriolis and centrifugal force/torque components, G q contains the gravity force/torque components, and τ is the vector of actuator forces.The second-order nonlinear coupled dynamic equation 4.3 can be linearized by inputting the inverse dynamics as follows: where qd is the vector of the desired accelerations, so that q qd .The resolved-acceleration control scheme 19 is  where K D and K P are gain matrices, subscript d denotes the desired value, ε r ε d − r is the positional error, and ε e is the orientation error.Unfortunately, this control scheme 4.5a breaks down when J −1 does not exist, which occurs at a singular configuration.

Damped Least-Square Method
In this section, the damped least-square method for the inverse problem of 4. where which always exists for ρ / 0. This solution is a compromise between the residual error, J q − ẋ, and the velocity q of the point P .Singular value decomposition SVD can provide insight into the singularities of the inverse Jacobian 19,22,23 .This study deals with nonredundant manipulators, and there are two orthogonal matrices U where Σ diag σ 1 , . . ., σ 6 and σ i are the singular values of J.The vectors u i and v i are the ith left and the ith right singular vectors, respectively.Substituting 4.9 into 4.8 gives where When the manipulator is far away from a singular point i.e., σ i ρ , we have σ i / σ 2 i ρ 2 ≈ 1/σ i , which implies that J * ρ ≈ J −1 .When the manipulator is near a singular point, it can be seen from 4.10 that the solutions for the velocity q have finite values.

The Hybrid Damped Resolved-Acceleration Control Scheme
In this paper, we propose a hybrid damped resolved-acceleration control scheme HDRAC that restricts velocity as well as acceleration.The HDRAC is used to solve the following optimization problem: where a * d is defined in 4.2 , ρ is the damping factor, and ρ hd is an additional damping factor which is dependent on the smallest singular value of the Jacobian matrix.The joint acceleration commands qhd have to minimize 4.11 .When the manipulator is far away from the singular point, there is no unnecessary joint velocity, so ρ hd will be zero.If the manipulator is at a singular point, the unnecessary joint velocity will have to be removed by joint acceleration commands, so ρ hd must be specified.The redundant components of the joint acceleration commands qhd will be equal to the unnecessary components of q multiplied by ρ hd to minimize 4.11 .where I 6 denotes a 6 × 6 identity matrix.In order to obtain the solutions for this optimization problem in a closed form, the discrete-time technique will be used.With discrete-time control, the back difference is always used to approximate the derivative, so that

Mathematical Problems in Engineering
where Δt is the sampling time.Through reduction and incorporation, the solution of the optimization problem is 19 Meanwhile, we have denoted the second term as qhdc .
The SVD theory can be reformulated as qhdc ρ hd The size of the deceleration region of the HDRAC can be adjusted by a linear function ρ hd , which is chosen as where σ min is the smallest singular value or singular parameter 19 of the Jacobian matrix.δ hd is a user-defined threshold value that represents the deceleration region.Unnecessary joint velocities will be partially removed when σ min is in this region and wholly removed at the next sampling interval when ρ hd 1/Δt.Note that the maximum value of ρ hd is 1/Δt.If the manipulator is at a singular point of σ i ρ > 0, i 1, . . ., k, and j k 1, . . ., 6, then 4.15 can be reduced to qhdc ≈ −1 Δt 4.17 We will next prove that this term removes unnecessary joint velocities at the next sampling interval.Multiplying both sides of 4.13 by v i v T i yields

4.18
Mathematical Problems in Engineering 13 By using an ideal computed-torque scheme, we force joint velocity v i i.e., v T i q along the degenerated direction to zero at the next sampling interval.The joint acceleration command along v i i.e., the product of v T i qhd should then be −1/Δt v T i q , which follows from 4.18 .This implies that 4.17 will remove unnecessary joint velocities along all degenerated directions v i , i k 1, . . ., 6.According to 4.18 , the desired joint velocities along v i , i 1, . . ., k are also decreased by the amount of Δtv T i qhd or ρ 2 / σ 2 i ρ 2 v T i q .Note that the maximum coefficient of v T i q v i in 4.15 is −1/Δt when σ i 0.

Stability Analysis
In this section, we will show that the HDRAC is asymptotically stable and discuss its convergent property.Let the orientation error be ε e f θ e u e , where u e and θ e are the rotational axis and the angle between the current orientation and the desired orientation, respectively, and f θ e θ e for very small values of θ e .Suppose that ρ > 0 when the manipulator is at a singular point.A sufficient condition for the global asymptotic convergence of the HDRAC in the whole workspace of the manipulator in a situation, where rd ṙd 0 and α d ω d 0, is that K D is a positive-definite matrix and where k pr and k pe are positive.The equilibrium is at the point where q 0. Let a Lyapunov function be 24 where ε r r d − r and E θ e 0 f φ dφ 1/2 θ 2 e , for f φ φ.Clearly, L q, q ≥ 0. It is recognized that θe u e ω d − ω .Differentiating E with respect to time, we obtain Ė ∂E ∂θ e θe ε e ω d − ω .

5.3
Evaluating ∂L/∂t in terms of the solutions of 5.2 yields If the ideal computed-torque control is used, q qhd , 4.14 can be substituted into 5.4 , which gives us 24

Mathematical Problems in Engineering
If rd ṙd 0 and α d ω d 0, then ∂L ∂t − qT J T K D J q − ρ hd ρ 2 qT q. 5.6 Since K D is a positive-definite matrix, ∂L/∂t ≤ 0 for all q, and ∂L/∂t 0 only for q 0. Combining the above results and L q, q ≥ 0, it is shown by the Lyapunov theorem that q is asymptotically stable, and the equilibrium point is q 0. However, for the proposed control scheme, ∂L/∂t has an additional term −ρ hd ρ 2 qT q more than that of the HDRAC, to ensure that ∂L/∂t is zero only for q 0, that is, the equilibrium point is q 0.

Numerical Simulations and Experimental Results
The purpose of inverse dynamics analysis is to determine the required torque of the separately-excited DC servomotor for a given motion trajectory of the moving platform.The method 19 is adopted to analyze the 3RPS parallel manipulator dynamics of the links and moving platform by using Newton-Euler equations of motion, which are formulated by considering the free-body diagrams of the links and moving platform separately.
Let us focus on the mass center O i of the ith link as shown in Figure 5, where i F pix , i F piy , and i F piz are three components of the weight force of the moving platform acting on the ith link, W i is the weight force of the ith link, i F 1motor is the acting force of the separatelyexcited DC servomotor on the revolution parts, and the i F 2motor is the acting force on each vertical leg.Each link is placed along the X i axis; therefore, its inertia tensor about X i , Y i and Z i is a diagonal matrix with diagonal elements I ixx ,I iyy , and I izz .The angular momentum 19 of the ith link is The rate change of H i is given by Ḣi Ḣi xyzi ω i × H i .

6.2
It follows from Figures 2 and 5 that the external moment about the mass center O i of the ith link is given by where i r pi is the vector from O i to B i with respect to the link frame with unit vector i i j i k i .The equations of rotational motion of the ith link can be obtained by equating the rate of change of H i with the external moment.By equating 6.2 with 6.3 , we obtain

Mathematical Problems in Engineering
İiyy ω iy I iyy ωiy − i r pix i F piz − i r piz i F pix , I iyy ω iy ω ix − I ixx ω ix ω iy İizz ω iz I iz ωiz i r pix i F piy − i r piy i F pix .

6.4
Since İixx İiyy İizz 0, i F pix passing through the mass center O i do not generate torques, and i r pi only has values in the X i -axis and couples with the separately-excited DC motor; the above equations 6.4 can be simplified as follows: where J m1 is the moment of inertia; B m1 is the damped coefficient; K T 1 and K T 2 are the torque constants; i a1 and i a2 are the armature currents.We now process to find the total force acting on the mass center O i as follows:

6.6
Once the force terms i F piy and i F piz have been obtained in 6.6 of link motion, one should take into consideration the force term i F pix .The angular momentum 19 of the moving platform is described by H p I pxx ω px i p I pyy ω py j p I pzz ω pz k p , 6.7 where I pxx , I pyy , and I pzz are, respectively, the mass moments of inertia of the moving platform about the X, Y , and Z axes.The rate change of angular momentum is The total torque acting on mass center P of the moving platform is given by where P P i is the vector from P to B i , i 1, 2, 3.

Numerical Simulations
In this section, numerical simulations are offered to demonstrate the applicability and effectiveness of the proposed algorithms of the 3RPS manipulator platform.The physical model is shown in Figure 6.The moving platform of the manipulator is a circle with radius b 22.08 cm, and the distance between any two of the three spherical joints on the moving platform is w 36.16cm.The fixed base is an inscribed triangle within a circle with radius a 41.72 cm.In the physical model, the constraints of the lengths d 1 , d 2 , and d 3 are ±10 cm, and the revolution angles θ 1 , θ 2 , and θ 3 are ±20 • .A program has been written in order to calculate the singular point.It was found that when d 1 0.0717 cm, d 2 9.7488 cm, and d 3 0.0771 cm, and θ 1 19.9834 • , θ 2 0.2438 • , and θ 3 19.9074• , then one singular point exists in the workspace.That is, the center coordinate of the moving platform with respect to the fixed base is 1.6144, −2.7933, 0.0308 , and the unit is cm.In this study, point P is required to move at a constant speed from the starting point 0, 0, 0 , to then pass through the singular point 1.6144, −2.7933, 0.0308 , and then it arrives at its final point 5, −5, 5 .
In this study, the HDRAC scheme is used to solve the damping factors ρ by means of the MPSO optimization method.Damping factors ρ are not fixed values and are used to damp the accelerations of point P when the manipulator is in the neighborhood of singularities.The dependence of the damping factors on the accelerations can change.For on-line applications, the optimization method can immediately provide optimal damping factors.
Figures 7 a -7 f show numerical simulations of the trajectories and velocities along the X-axis, Y -axis, and Z-axis for the controller incorporating HDRAC with fixed damping values, and for the controller incorporating HDRAC with MPSO.For the controller incorporating HDRAC with MPSO, when the manipulator is near the singular point, the system functions are more smooth and there is no fluctuation.

Experimental Setup
In order to demonstrate the proposed control rules, PC-based experimental equipment was devised.The parallel manipulator was coupled with a DC motor.The translation position was measured by a sensor.A photo of the experimental equipment is shown in Figure 6.
In general, such parallel robotic is controlled by hydraulic servo, that is, to control the extension or retraction of the limb through driving with hydraulic system.Limiter sensor or measuring tools are generally installed on the link of the limb to enable accurate changes in its length when controlling positions; the characteristic of position limiter sensor is shown in the Table 1.The controller is based on a PC with Pentium-586 CPU.In order to measure the angle, angular speed, position, and velocity of the manipulator, the interface device was implemented by the motion control card PCI-7344.The graphical software of Simulink and the control software LABVIEW were used to implement the proposed control rule.At the same time, this software also determines the linear converts between the physical scales and measured voltages from sensors.A block diagram of the overall system with the controller incorporating HDRAC with MPSO and acceleration feedback is illustrated in Figure 8  the 3RPS parallel manipulator was divided into three parts: an actuator, parallel manipulator, and controller.The actuator part consisted of a DC motor and a driver.The parameters of the DC motor were K t 0.21 Nm/A, J m 0.000037 Nms 2 , B m 0.000013 Nms/rad, l a 2.4 mH, and r a 4.8 Ω.

Experimental Results
Figures 9 a -9 f show the experimental result of trajectory and velocity along the Xaxis, Y -axis, and Z-axis with the DARAC controller that was obtained by using the 3RPS parallel manipulator, in which the end effector passed through a singular point.The DARAC controller functioned, such that joint acceleration along the degenerated direction was damped when the end effector was close to the singular point.It was found that when the manipulator approached the singular point, the system was unstable for the DARAC controller and displayed great fluctuations.Figures 10 a -10 f show the experimental result of trajectory and velocity along the X-axis, Y -axis, and Z-axis for the controller incorporating HDRAC with fixed damping values and for the controller incorporating HDRAC with MPSO.It was found that, by using the controller incorporating HDRAC with MPSO, a very good performance could be achieved under the same conditions and that the system was generally stable and had very slight fluctuations.These experimental results demonstrate that the HDRAC overcomes the problem of fluctuation when the manipulator arrives at the singular point.It also proves that the proposed controller incorporating HDRAC with MPSO has great robustness.

Conclusion
This paper proposes a simple HDRAC scheme with an MPSO method, which has the merits of offering a high-quality solution and better computation efficiency in the neighborhood of a singular point.The proposed method directly damps velocities and accelerations, such that unnecessary joint velocities are removed at a singular point.The proposed HDRAC scheme with MPSO control is asymptotically stable and remedies the problem of fluctuation.In this study, HDRAC with MPSO is used to solve the damping factors ρ by the optimization method.The damping factors ρ are not fixed values for on-line applications, and the optimization method can immediately provide optimal damping factors.Because it shortens the detection processes, the MPSO method can satisfactorily predict optimal damping values.
Comparisons between our numerical simulations and experimental results showed that the proposed method offers greater stability and efficiency and solves the problems of instability and fluctuation more easily and quickly.The main advantage of using a controller incorporating HDRAC with MPSO is that the 3RPS parallel manipulator control system does not need to plan its path to avoid any singular point, because the controller will command the manipulator to move near the singular point with minimum trajectory errors.Moreover, in comparison with a controller incorporating HDRAC with fixed damping values, the proposed approach makes it easier and simpler to improve the performance of the motion in the neighborhood of singular points.

A.
The items contained in matrix D are as follows: The items contained in matrix Q are as follows: A.2

B.
The items contained in matrix A are as follows:

C.
The coefficients are as follows:

Figure 3 :
Figure 3: The concept of distances.
population and evaluate the fitness function Create the new feasible population, repair strategy, and evaluate the fitness function Does the calculated fitness function value improve future?Update the best solution Sort and select Are the stopping rules satisfied?Output the optimal solution

Figure 4 :
Figure 4: The flow chart of the MPSO.

Figure 5 :
Figure 5: Free-body diagram of acting forces on the ith link.

Figure 6 :
Figure 6: Physical model of the 3RPS parallel manipulator.

Table 1 :
The characteristic of position limiter sensor.Static platform error due to arm position sensor errors mm , this error represents the vector conversion error from joint coordinates to the platform reference coordinates Translation

Figure 7 :Figure 8 :
Figure 7: Numerical simulations for the controller incorporating HDRAC with fixed damping values and the controller incorporating HDRAC with MPSO. a Translation responses along the X-axis.b Velocities along the X-axis.c Translation responses along the Y -axis.d Velocities along the Y -axis.e Translation responses along the Z-axis.f Velocities along the Z-axis.

Figure 9 :
Figure 9: Experimental results for the controller incorporating DARAC with fixed damping values.a Translation responses along the X-axis.b Velocities along the X-axis.c Translation responses along the Y -axis.d Velocities along the Y -axis.e Translation responses along the Z-axis.f Velocities along the Z-axis.

Figure 10 :
Figure 10: Experimental results for the controller incorporating HDRAC with fixed damping values and the controller incorporating HDRAC with MPSO. a Translation responses along the X-axis.b Velocities along the X-axis.c Translation responses along the Y -axis.d Velocities along the Y -axis.e Translation responses along the Z-axis.f Velocities along the Z-axis.