2.1. Basic slicing algorithm

To compactly cover the target surface, path points should be generated by processing corresponding points on the target surface. This is commonly implemented by slicing the model with a set of planes to extract the intersection points. It imitates the slicing process in additive manufacturing but with a much longer spacing between slice planes. In spray painting, most of the path is in a raster pattern, which ensures a better uniformity of coating thickness [Reference Kao and Prinz14, Reference Gleeson, Jakobsson, Salman, Ekstedt, Sandgren, Edelvik, Carlson and Lennartson15]. Some key parameters of the raster pattern path could be presented by these slice planes. For example, spacing distance between two adjacent passes is expressed by that between two adjacent planes.

By importing the STL model of the target surface and transforming it to fit the XZ plane of the world coordinate system, a group of slice planes could be defined by their Z coordinates. These planes are perpendicular to the XZ plane and parallel to each other. Based on the plane, the intersection points on STL facets’ edges could be calculated by:

(1a) \begin{align} x_{\textrm{inter}} & =x_1 + \frac{({z_{\textrm{sli}}} -{z_1})({x_2} -{x_1})}{{z_2} -{z_1}} \end{align}

(1b) \begin{align} y_{\textrm{inter}} & = y_1 + \frac{({z_{\textrm{sli}}} -{z_1})({y_2} -{y_1})}{{z_2} -{z_1}} \end{align}

(1c) \begin{align} z_{\textrm{inter}} & = z_{\textrm{sli}} \end{align}

where $({x_{\textrm{inter}}},{y_{\textrm{inter}}},{z_{\textrm{inter}}})$ is the coordinate of the resulting intersection point, $({x_1},{y_1},{z_1})$ and $({x_2},{y_2},{z_2})$ are coordinates of the edge’s two vertices, and $z_{\textrm{sli}}$ is the Z coordinate of the current slice plane.

In the STL model, a vertex is shared with several adjacent facets. For each vertex, let us calculate the average normal vector among its adjacent facets. This vector is defined as the normal vector of the vertex. Moreover, the normal vector of the intersection point could be constructed by the average normal vector between the sliced edge’s two vertices. By offsetting each intersection point along the normal vector with a certain spray distance, a path pass for spray painting could be generated.

This method works when the slice plane is within the upper and lower boundaries. However, when the plane slices some special places around the boundary, for instance, the highest or lowest vertex of the model, the basic method could only generate limited path points compared with the middle passes. This phenomenon is called “boundary effect.” Besides, the basic method lacks further processing to boost path’s smoothness, making path points discrete and non-uniform. The above problems would cause an unsatisfactory thickness distribution and weaker robot control [Reference Hegels, Wiederkehr and Müller16]. To ameliorate these two situations, the boundary fitting demand and path smoothness demand should be met in the proposed method. Figure 1 illustrates the outline of the basic slicing algorithm.

Figure 1.

Outline of basic slicing algorithm.

2.2. Thickness modeling

Based on the generated pass, the following work is calculating the resulting coating thickness. This is a preliminary step for further thickness uniformity optimization. To simulate it, two parts should be considered here: the instant model and the dynamic model.

The instant model aims to compute the coating thickness formed by a spray gun targeted at an arbitrary surface in an instant. To model such a scenario, the first step is finding the standard thickness distribution. It is formed by a spray gun targeted at a standard plane with a perpendicular spray direction and a standard spray distance. To implement it, an experimental test should be conducted beforehand to collect related thickness data. Afterwards, a distribution function is adopted to fit the data. In this study, the ellipse dual- $\beta$ distribution model is selected owing to its flexibility in describing the shape of the thickness distribution [Reference Yu, Cheng, Zhang and Ou7, Reference Zhang, Huang, Gao and Wang17]. It could be given by:

(2) \begin{equation} \begin{aligned} T_{\textrm{std}}\!\left ( p \right ) & = T_{\textrm{std}}\!\left ({{x_p},{y_p}} \right ) \\[5pt] & ={T_{\textrm{max}}}{{\left ({1 - \frac{{x_p^2}}{{{a^2}}}} \right )}^{{\beta _1} - 1}}{{\left [{1 - \frac{{y_p^2}}{{{b^2}\!\left ({1 - \frac{{x_p^2}}{{{a^2}}}} \right )}}} \right ]}^{{\beta _2} - 1}},\; & \frac{{x_p^2}}{{{a^2}}} + \frac{{y_p^2}}{{{b^2}}}\le 1 \end{aligned} \end{equation}

where $T_{\textrm{max}}$ is the maximum coating thickness in an instant, it is related to the paint flow rate and paint transfer efficiency of the spray gun. $a$ is the semi-major axis of the elliptical coverage area, and $b$ is the semi-minor axis. They could describe the cross-section contour of the distribution. $\beta _1$ is the exponent of the distribution on the Y-axis section and $\beta _2$ is that on the X-axis section. They could delineate the contour on the longitudinal section.

To make the model more general, the next step is projecting the distribution function from the standard plane onto an arbitrary surface point. This course is based on differential geometry [Reference Conner, Atkar, Rizzi and Choset18]:

(3) \begin{equation} T\!\left ( p_1 \right ) = \begin{cases} T_{\textrm{std}}\!\left ({{p}} \right )\frac{{{h^2}cos\alpha }}{{{{{\overrightarrow{gp_1} }}^2}{{\cos }^3}\theta }} &,{\rm{ }}\alpha \le 90^\circ \\[5pt] 0 &,{\rm{ }}\alpha \ge 90^\circ \end{cases} \end{equation}

where $g$ is the path point, $p_1$ is the target point, $\;\overrightarrow{gp_1} \;$ is the vector from $g$ to $p_1$ , ${P}\;$ is the standard plane mentioned above, ${p}\;$ is the intersection point between $\;\overrightarrow{gp_1} \;$ and $\;{P}$ , $h$ is the standard spray distance between $g$ and $\;{P}$ , which has been determined in the experimental test, $\theta$ is the angle between spray direction and $\;\overrightarrow{gp_1}$ , and $\alpha$ is the angle between the normal vector of $p_1$ and the opposite direction of $\overrightarrow{gp_1}$ . Hence, an instant model could be constructed as discussed above. The instant model is illustrated in Fig. 2.

Figure 2.

Instant model. The left figure illustrates the standard thickness distribution and the right figure depicts the projection from the standard plane to an arbitrary surface point.

In contrast to the instant model, the dynamic model is targeted at finding the coating thickness formed by a spray gun moving along a certain pass in a certain time interval. Dynamic model during time interval $[0,t]$ could be configured by integrating instant model $\;T\;$ over time.

(4a) \begin{equation} T_{\textrm{dyn}}=\int _0^tTdt \end{equation}

To calculate it, a continuous path is divided into enough discrete segments where each segment is a line segment bounded by two adjacent path points [Reference Arıkan and Balkan19]. Hence, an instant time could be approximated as a time interval as follows:

(4b) \begin{equation} dt = \frac{|\overrightarrow{g_ig_{i+1}}|}{v} \end{equation}

where $\; |{\overrightarrow{{g_i}{g_{i + 1}}} } |\;$ is the length of the path segment bounded by path points $\;{g_i}\;$ and $\;{g_{i + 1}}$ , and $v$ is the average spray gun motion speed in this segment. In terms of $p_1$ , the cumulative thickness ${T_{\textrm{dyn}}} ( p_1 )$ could be computed by summing all the products between the instant thickness $\;{T_{i + 1}} ( p_1 )\;$ gained in the discrete path segment $\;{g_i}{g_{i + 1}}\;$ and the corresponding time interval.

(4c) \begin{equation} {{T_{\textrm{dyn}}}\!\left ( p_1 \right ) = \mathop \sum \limits _{i = 1}^{n - 1}{\rm{ }}{T_{i + 1}}\!\left ( p_1 \right )\frac{{\left |{\overrightarrow{{g_i}{g_{i + 1}}} } \right |}}{v}} \end{equation}

where $i$ and $n$ are the index and the total number of the path points, respectively. This model is depicted in Fig. 3.

Figure 3.

Dynamic model.

3.1. Boundary point detection

To meet the demand for boundary fitting, the first step is identifying the boundary points of the target surface, which are evident for human beings but a little tricky for computers. In the STL model, an edge is defined as a line segment bounded by two vertices on a certain facet. According to the topology of a meshed model, each boundary edge is exclusive to only one facet while other edges are shared with different facets. Based on this principle, the whole loop of boundary edges could be searched and boundary points on these edges are extracted. A boundary point is shared with two boundary edges and the other two boundary points on these edges are its neighboring boundary points. Figure 5(a) depicts this course.

Figure 5.

Boundary point detection. (a) Boundary point extraction. (b) Segmentation.

Generally, a boundary loop could be divided into four independent segments by four corner points. To identify the corner points, let us connect each boundary point’s two boundary edges and compare the angle between these two edges as shown in Fig. 5(b). The point with one of the largest four angles is a corner point. Afterwards, to initialize the segmentation of boundary points, let us randomly choose one corner point as the initial point on the initial segment and then randomly choose one of its neighboring boundary points as the second point. Afterwards, by excluding the former point on the segment, the remaining neighboring boundary point of the current point is selected as the next point. This step repeats until another corner point is searched. The next segment begins with the new corner point and a new step continues in the same way. This loop ends when the initial corner point is searched again. Figure 5(b) also illustrates this course.

Thus far, the four independent boundary segments could be distinguished. Among them, the segment with the highest average Z coordinate is defined as the upper boundary, whereas the segment with the lowest one is defined as the lower boundary.

3.2. Sample point generation

Based on the intersection points mentioned in Section 2.1, several typical points should be selected from them to reflect the general geometric information and facilitate further processing. The term “key points” is used here to describe such points. They are selected at a uniform length interval between the leftmost and the rightmost intersection point of a pass.

As shown in Fig. 6, the search starts from the leftmost intersection point, which is selected as the first key point. During each step, the length of the line segment adds up one by one. When the accumulative length of segments exceeds the length interval $invl$ , the two vertices of the current segment remain to be selected. If the current accumulative length is closer to $invl$ , the right vertex is selected as the next key point. Otherwise, the left vertex is selected. In the next step, the accumulative length is cleared and the searching continues from the new key point. At the end of this loop, the rightmost intersection point is selected as the last key point. By offsetting each key point along its normal vector with a value of spray distance, a corresponding sample point is formed, whose spray direction is opposite to the normal vector:

(5a) \begin{equation} \begin{bmatrix} x_{\textrm{smp}} \\[5pt] y_{\textrm{smp}} \\[5pt] z_{\textrm{smp}} \\[5pt] 1 \\[5pt] \end{bmatrix} = \begin{bmatrix} I_{3\times 3} & h_1\cdot \overrightarrow{n_{\textrm{key}}} \\[5pt] 0\;0\;0 & 1 \end{bmatrix} \begin{bmatrix} x_{\textrm{key}} \\[5pt] y_{\textrm{key}} \\[5pt] z_{\textrm{key}} \\[5pt] 1 \\[5pt] \end{bmatrix} \end{equation}

(5b) \begin{equation} \overrightarrow{{n_{\textrm{smp}}}} = - \overrightarrow{{n_{\textrm{key}}}} \end{equation}

where $({x_{\textrm{key}}},{y_{\textrm{key}}},{z_{\textrm{key}}})$ and $\overrightarrow{{n_{\textrm{key}}}}$ are the coordinate and normal vector of the key point, $({x_{\textrm{smp}}},{y_{\textrm{smp}}},{z_{\textrm{smp}}})$ and $\overrightarrow{{n_{\textrm{smp}}}}$ are those of the sample point, $I_{3 \times 3}$ is a $3\times 3$ identity matrix, and $h_1$ is spray distance. Figure 6 depicts the generation of sample points from key points. By offsetting different key points with different spray distances, various shapes of the pass could be generated based on the resulting sample points.

Figure 6.

Sample point generation. Key points are searched with a uniform length interval. The right vertex of Segment 8 is selected as the key point because the accumulative length until Segment 8 just exceeds the length interval and is closer to it than the former one.

3.3. Curve fitting

Fitting boundary points and sample points extracted above, the pass could be expressed by a polynomial curve [Reference Noakes and Popiel20]. Generally, a space curve could be defined by a parametric equation:

(6a) \begin{equation} \overrightarrow r = (x(t),{\rm{ }}y(t),{\rm{ }}z(t)) \end{equation}

Firstly, let $x$ be a free parameter:

(6b) \begin{equation} x(t) = t \end{equation}

Regarding the $y$ - $t$ relationship, a fourth-order polynomial is used to fit the sample points. It could be approximated by:

(6c) \begin{equation} y(t) ={c_0} +{c_1}t +{c_2}{t^2} +{c_3}{t^3} \end{equation}

where $c$ is the coefficient of the $y$ - $t$ polynomial.

In terms of the $z$ - $t$ relationship, another fourth-order polynomial should be introduced to meet the demand for boundary fitting. Therefore, coordinates of the upper and lower boundary points are utilized in the curve fitting. Firstly, let us assume that these two groups of points are translated with the same distance along the z axis to the place near the slice plane:

(6d) \begin{equation}{z_{\textrm{upr}}} - ({{\bar{z}_{\textrm{upr}}} -{z_{\textrm{sli}}}})\,=\,{a_0} +{a_1}{x_{\textrm{upr}}} +{a_2}x_{\textrm{upr}}^2 +{a_3}x_{\textrm{upr}}^3+{a_4}x_{\textrm{upr}}^4 \end{equation}

(6e) \begin{equation} z_{\textrm{lwr}} + ({{z_{\textrm{sli}}} -{\bar{z}_{\textrm{lwr}}}})\,=\,{b_0} +{b_1}{x_{\textrm{lwr}}} +{b_2}x_{\textrm{lwr}}^2 +{b_3}x_{\textrm{lwr}}^3 +{b_4}x_{\textrm{lwr}}^4 \end{equation}

where $x_{\textrm{upr}}$ , $x_{\textrm{lwr}}$ , and $z_{\textrm{upr}}$ , $z_{\textrm{lwr}}$ are each X and Z coordinates of the upper and lower boundary points, $\bar{z}_{\textrm{upr}}$ and $\bar{z}_{\textrm{lwr}}$ are the average Z coordinates of the upper and lower boundary points, $z_{\textrm{sli}}$ is the Z coordinates of the slice plane indicating the rough Z-axis position of the pass, $a$ is the coefficient of the translated upper boundary polynomial, and $b$ is the coefficient of the lower one. To offer more possibilities of the generated path, the first or the last $z_{\textrm{sli}}$ could exceed the maximum or minimum Z coordinate of the model. In this case, its y−t relationship is constructed by fitting corresponding upper or lower boundary points instead of sample points.

Besides, the shape of the curve should vary gradually from pass to pass. It ensures an evener overlap and brings about a more uniform thickness distribution. To fulfill it, a compound polynomial is constructed by a linear interpolation between two coefficients with the same order from the two translated polynomials:

(6f) \begin{equation} {d_i}=\frac{{{z_{\textrm{sli}}} -{\bar{z}_{\textrm{lwr}}}}}{{{\bar{z}_{\textrm{upr}}} -{\bar{z}_{\textrm{lwr}}}}}{a_i} + \frac{{{\bar{z}_{\textrm{upr}}} -{z_{\textrm{sli}}}}}{{{\bar{z}_{\textrm{upr}}} -{\bar{z}_{\textrm{lwr}}}}}{b_i},\quad i \in \left \{{0,1,2,3,4} \right \} \end{equation}

(6g) \begin{equation} z(t) ={d_0} +{d_1}t +{d_2}{t^2} +{d_3}{t^3} +{d_4}{t^4} \end{equation}

where $i$ is the index of the polynomial coefficient, and $d$ is the compound coefficient. Eq. (6g) denotes the final $z$ - $t$ polynomial. In the curve fitting, the fourth-order polynomial is selected to describe the path. As shown in Eq. (6f), it could facilitate further construction by processing its five explicit and dominating coefficients. Besides, the fourth order is sufficient to approximate various path passes. Meanwhile, it could also avoid Runge’s phenomenon where the generated path could be highly oscillatory when the degree comes to be higher. Figure 7 depicts the $z$ - $t$ and $y$ - $t$ polynomial curves.

Figure 7.

$z$ - $t$ and $y$ - $t$ polynomial curves. $x(t)=t$ .

Thus far, with an arbitrary $t$ , the coordinates of path points could be obtained based on Eq. (6a). Owing to the varying curvature of the path curve, points should be interpolated with a uniform length interval as explained in Section 3.2, instead of directly plugging a series of uniform $t$ into the polynomial. It ensures a uniform distribution of the path points.

3.4. Path coordinate system definition

Up to now, a series of consecutive path points with certain positions have been generated. However, orientation, another significant parameter to define a robot path, has not been expressed. In robotics, a common way to describe these two parameters is attaching a coordinate system called {Path} to each path point. In this scenario, the origin of {Path} is the position of the path point, the Y axis of {Path} is along the spray gun’s motion direction, the Z axis of {Path} is along the spray direction, and the X axis is orthogonal to both of the Y axis and Z axis. In this case, the motion direction is defined as the direction of a connecting line between the current path point and the next one. Meanwhile, the spray direction is mapped from the opposite normal vector of the closest intersection point to the path point.

In the definition of the coordinate system, the orthogonality between the motion direction and the spray direction, the two independent axes, is required to be ensured. It could be fulfilled by rotating the spray direction to be perpendicular to the motion direction. Let $\overrightarrow{{p_1}{p_2}}$ represent the vector from path point $p_1$ to its adjacent path point $\;{p_2}$ , $\;\overrightarrow{{n_1}}$ represents the original spray direction of $p_1$ , $\;\overrightarrow{{i_{\{\textrm{path}\}}}} \;$ , $\;\overrightarrow{{j_{\{\textrm{path}\}}}} \;$ $\;\overrightarrow{{k_{\{\textrm{path}\}}}} \;$ are the standard bases of the X axis, Y axis, and Z axis in {Path} attached to $p_1$ , and $\theta$ is the angle between the original spray direction and the adjusted one:

(7a) \begin{equation}{\theta = \arccos \!\left ({\overrightarrow{{n_1}} \cdot \overrightarrow{{j_{\{\textrm{path}\}}}} } \right ) - \frac{\pi }{2}} \end{equation}

hence $\overrightarrow{{k_{\{\textrm{path}\}}}}$ is given by:

(7b) \begin{equation}{\overrightarrow{{k_{\{\textrm{path}\}}}} = \frac{{{R_{\overrightarrow{ i_{\textrm{\{Path\}}}} }}\left ( \theta \right )\overrightarrow{{n_1}} }}{{\left |{\overrightarrow{{n_1}} } \right |}}} \end{equation}

where $R_{\overrightarrow{i_{\textrm{\{Path\}}}} }(\theta )$ is the desired rotation matrix that rotates $\;\overrightarrow{{n_1}} \;$ about $\overrightarrow{{i_{\{\textrm{path}\}}}}$ to be orthogonal to $\;\overrightarrow{{j_{\{\textrm{path}\}}}}$ . It could be obtained by:

(7c) \begin{equation} R_{\overrightarrow{i_{\textrm{\{Path\}}}} }(\theta )={{I}_{3\times 3}}+\textrm{sin}\theta \overrightarrow{i_{\textrm{\{Path\}}}}+\left ( 1-\textrm{cos}\theta \right ){{{\overrightarrow{i_{\textrm{\{Path\}}}}}}\,^{2}} \end{equation}

Therefore, {Path} could be completely defined by a 4 $\times$ 4 transformation matrix [Reference Liu, Xu, Zhou and Fang21]:

(7d) \begin{equation}{T = \left [{\begin{array}{c@{\quad}c@{\quad}c@{\quad}c}{\overrightarrow{{i_{\{\textrm{path}\}}}} } &{\overrightarrow{{j_{\{\textrm{path}\}}}} } &{\overrightarrow{{k_{\{\textrm{path}\}}}} } &{{P_{\{\textrm{path}\}}}} \\[5pt] 0 & 0 & 0 & 1 \end{array}} \right ]} \end{equation}

where $P_{\{\textrm{path}\}}$ is the origin of {Path}. These coordinates are expressed from {W}, the world coordinate system. Figure 8 illustrates the definition of {Path}.

Figure 8.

Definition of {Path}.

4.1. Influencing factors

With parameters of the generated path and the model to calculate the resulting coating thickness, an optimization could be carried out to obtain the optimal uniformity of coating thickness. In spray painting, a variety of parameters could impact uniformity. Among them, spray distance, spray speed, and spacing distance are three of the most dominating factors related to path generation.

As discussed in Section 3.2, spray distance is introduced there to express the offset distance from the key point to the sample point along the normal vector. Hence, spray distance affects the shape of the pass, which fits the sample points. Besides, a higher spray distance would bring about a lower coating thickness and a wider coverage area, and vice versa. To gain a uniform coating thickness, it is workable to adjust spray distance for each key point.

Spray speed directly affects the dynamic model of the thickness simulation as shown in Eq. (4c). The higher the spray speed, the lower the coating thickness, and vice versa. In this study, the spray gun is assumed to move at a constant speed along a certain pass.

Differing from the above two parameters, spacing distance impacts the uniformity indirectly. A standard thickness distribution is in a volcano shape where the thickness value is higher in the middle while lower on the two tails. It makes the thickness uneven. To alleviate it, coating thickness distributions resulting from two adjacent passes should overlap each other as shown in Fig. 9. It would accumulate a higher coating thickness on the side of the coverage area, making the compositive thickness more uniform. In this scenario, spacing distance impacts the degree of overlap. Therefore, it also has an influence on the thickness uniformity.

Figure 9.

Two overlapped thickness distributions.

4.2. Objective function

To evaluate the uniformity of coating thickness, the coefficient of variation is selected to be the indicator [Reference Guan and Chen11]:

(8a) \begin{equation} {{c_v} = \frac{{{\sigma _{\textrm{thk}}}}}{{{\mu _{\textrm{thk}}}}}}\times 100\% \end{equation}

where $\;{c_v}\;$ is the coefficient of variation, $\;{\sigma _{\textrm{thk}}}\;$ and $\;{\mu _{\textrm{thk}}}\;$ are the standard deviation and the mean of the coating thickness, respectively. As $c_v$ is a dimensionless number, it is feasible to evaluate the uniformity among thickness datasets with different units and different desired values. If a specific desired thickness $T_{\textrm{desired}}$ is determined, $\;{\mu _{\textrm{thk}}}\;$ is replaced by this value and $\;{\sigma _{\textrm{thk}}}\;$ is revised to be the standard deviation from $T_{\textrm{desired}}$ .

In spray painting, the overall coating thickness is accumulated pass by pass and the global uniformity of the target surface is a result of the local uniformity between each neighboring pass pair. To consider the relationship between the global uniformity and the local one, the objective function here is constructed as a combination of them. For the first pass and the last pass, the global uniformity is selected as the objective function. It considers the thickness of the whole surface points. In terms of the pass within them, the corresponding objective function is the local uniformity, which only includes the surface points between the current pass and the last neighboring one. Hence, the final objective function is given by

(8b) \begin{equation} \min f = \begin{cases} c_{v_{\textrm{global}}} &,\textrm{the first or the last pass} \\[5pt] c_{v_{\textrm{local}}} &,\textrm{otherwise.} \end{cases} \end{equation}

Considering the parameters mentioned in Section 4.1, the optimization problem could be defined as:

(8c) \begin{equation}{\min f\left ({H,v,d} \right )} \end{equation}

where $H$ is an array stores spray distance of each sample point on the current pass, $v$ is the spray speed of the current pass, and $d$ is the spacing distance between the current pass and the last pass. Figure 10 outlines the factors in the optimization.

Figure 10.

Factors in the optimization. The current pass is Pass 2. $H$ is an array stores spray distance of each sample point on the current pass, $v$ is the spray speed of the current pass, and $d$ is the spacing distance between the current pass and the last pass.

4.3. Optimization method: PSO

Even though the problem has been defined in Eq. (8a), the objective function $min f$ is non-continuous and not specific. Therefore, the original optimization method requires the derivative of the objective function is not workable in this case. To deal with it, a global search algorithm needs to be adopted, which only requires calculating the value of the objective function during running. Among global search algorithms, the PSO [Reference Wang, Tan and Liu22, Reference Bouraine and Azouaoui23] is preferred in spray painting problems owing to its convergence efficiency and feasibility for a large number of optimization variables [Reference Guo, Chen, Zhang and Wei24]. Thus, it is selected to solve the optimization problem here.

PSO is a bio-inspired algorithm. It heuristically imitates a flock of birds when they forage in groups. In PSO, a particle is a position consisting of the optimization variables. In terms of the initialization of PSO, the algorithm would randomly form a swarm of particles within the search space. The number of particles in the swarm is defined as “size of swarm.” For each particle, the position with the minimal value of the objective function so far is called “particle’s best-known position.” For the whole swarm, the position with the minimal value so far among all of the particles is called “swarm’s best-known position.”

At each step of PSO, a particle would move to a new position toward a certain direction. This direction is constructed by the vector toward the particle’s best-known position, the vector toward the swarm’s best-known position, and the vector from the last position to the current position. The influence degrees of these three factors are expressed by the cognitive coefficient (between [0,4]), social coefficient (between [0,4]), and inertia weight (between [0.1,1.1]). By this means, the position of each particles is iteratively updated. Finally, these particles would gradually converge around where the global minimum of the objective function is.