Abstract
This paper presents a novel constraint handling technique for optimum path generation of four-bar linkages using evolutionary algorithms (EAs). Usually, the design problem is assigned to minimize the error between desired and obtained coupler curves with penalty constraints. It is found that the currently used constraint handling technique is rather inefficient. In this work, we propose a new technique, termed a path repairing technique, to deal with the constraints for both input crank rotation and Grashof criterion. Three traditional path generation test problems are used to test the proposed technique. Metaheuristic algorithms, namely, artificial bee colony optimization (ABC), adaptive differential evolution with optional external archive (JADE), population-based incremental learning (PBIL), teaching-learning-based optimization (TLBO), real-code ant colony optimization (ACOR), a grey wolf optimizer (GWO), and a sine cosine algorithm (SCA), are applied for finding the optimum solutions. The results show that new technique is a superior constraint handling technique while TLBO is the best method for synthesizing four-bar linkages.
1. Introduction
Since the last decade, many researchers have tried to solve the optimization for path generation of four-bar linkages using metaheuristic (MH) algorithms. The objective of path generation problem is to find dimensions of a mechanism, which minimize the target path and the actual path of a point on the coupler link. Path synthesis is one type of kinematic syntheses of four-bar mechanisms [1–13] in which such syntheses are basically classified into two groups. The first category is called dimensional synthesis [1, 2, 4, 5, 7–10]. This synthesis type aimed to find significant link lengths to achieve desirable function, path, and motion generation. The second synthesis is called type synthesis [6, 11] where a designer initially specifies a predefined motion transmission and is supposed not initially to know the mechanism type. This method is analogous to topology design in structural optimization. Having finished synthesizing, a certain mechanism type is received. Position analysis of the four-bar mechanism can be categorized into two groups. The first one is a vector loop or loop closure equation, which is the most traditional method in kinematic analysis, and it is one of the most popular analyses for path synthesis [1, 2, 4, 5, 7–15]. This equation can be solved by using Freudenstein equation. The second analysis technique is a straight forward and a simple method for position analysis involving the use of trigonometric laws for triangles, e.g., the law of cosine [3, 16, 17], whereas the six-bar linkage for steering mechanism also uses the same technique [18]. This work proposes a new computing technique for four-bar linkage position analysis by employing the concept of drawing an arbitrary rectangle using two circles.
The mechanism synthesis can be converted into optimization problem and be solved by using optimizers, where both nongradient- and gradient-based algorithms have been solved this problem. Recently, a nongradient-based optimizer, e.g., evolutionary algorithms (EAs) or metaheuristics (MHs), is a more popular selection in solving such optimization problems. It has been found that the advantages of using MHs are robustness, simplicity of use, and independence of function derivatives; however, they unavoidably lack convergence speed and consistency. At present, many algorithms in this group have been developed, which are expected to enhance in both convergence speed and consistency. Some of the most frequently used MHs for path synthesis are differential evolution (DE) [2, 3, 5, 8, 9, 11, 13], genetic algorithms (GA) [5, 13], particle swarm optimization (PSO) [5, 13], and an imperialist competitive algorithm (ICA) [13], etc. The use of gradient-based method, on the other hand, is somewhat questionable to deal with global optimization and nonsmooth constraints in the path synthesis. Nevertheless, if those aforementioned factors can be alleviated, the advantages of the gradient-based method are better convergence rate and consistency. In the literature, many researchers have combined MHs and a gradient-based optimizer for solving many kinds of real world problems, which is called a hybrid algorithm. Especially for part synthesis, the hybrid optimizers are introduced as the ant gradient [1], hybrid GA [4], and hybrid GA with sequential quadratic programming (SQP) [12]. The hybridization between two or more MHs was also studied such as hybrid GA-DE [7]. Furthermore, the path synthesis is extended to multiobjective optimization, which was solved by using a multiobjective genetic algorithm (MOGA) [10]. From the review literature, it was found that some MHs have been used for solving this task except the work by Sleesongsom and Bureerat [17]; therefore, one of the objectives of this paper is to present the comparative performance of a number of currently used MHs. Those algorithms include the artificial bee colony optimization (ABC), adaptive differential evolution with optional external archive (JADE), population-based incremental learning (PBIL), teaching-learning-based optimization (TLBO), the real-code ant colony optimization (ACOR), a grey wolf optimizer (GWO), a Jaya algorithm (Jaya), and a sine cosine algorithm (SCA).
The path generation is a mechanism synthesis to make a point on a coupler link move along the target path; thus, the objective function is the minimization of the sum square error between the target path and the actual path [5]. The design problem is a constrained optimization problem that comprises two constraints. The first constraint is set for the shortest link in the mechanism to be able to rotate with a complete revolution (crank) in either direction (clockwise or counter clockwise). The second constraint is assigned such that the four link lengths satisfy the Grashof criterion which results in a crank-rocker. From previous work, a simple exterior penalty function technique had used to deal with these constraints [1–5, 7, 8, 10–13]. The new technique proposed by [17] to neglect the first constraint from the optimization problem, which found new technique, provided the better result than the traditional exterior penalty technique. Additionally, the technique had improved in the result, but it increased in time consuming. From the present study can be concluded that the constraint handling technique is an inefficient technique, which needs an improvement [9, 13, 14]. Phukaokaew et al. [14] studied the number of unsuccessful runs from performing MHs for 30 runs, where the path synthesis optimization problems employ the penalty function (PF) technique. This means there is no guarantee that using this technique can promote the good results. The reason is that using the penalty technique leads to an overly narrow feasible region. As a consequence, MHs, which mostly have slow convergence rates, struggle to reach an optimum. Ebrahimi and Payvandy proposed the way to improve the constraint handling technique, which was still based on a penalty function, and they believe the proposed method can enhance the search performance [13].
This paper focuses on two aspects of investigation. Firstly, a new constraint handling technique for path synthesis of a four-bar linkage using MHs is proposed. The method is based on repairing illegitimate design solutions to become feasible solutions. The second investigation is to study the performance comparison of a number of established MHs for solving four-bar linkage path synthesis with the new constraint handling technique, where both convergence rate and consistency of the methods are measured.
The rest of this paper is organized as follows. Section 2 proposes an alternative position analysis of a four-bar linkage. The optimization problem and the constraint handling technique are detailed in Section 3. The numerical experiments are given in Section 4, while the design results are discussed in Section 5. The discussions and conclusions of the study are finally drawn in Sections 6 and 7, respectively.
2. Position Analysis of Four-Bar Linkages
The kinematic diagram of a four-bar linkage is shown in Figure 1. The four-bar linkage is the simplest and most commonly used linkage in many engineering applications. It is composed of a kinematic chain of four binary links connected with four revolute joints (denoted by capital letters) with one link being assigned as a frame. Applications for this mechanism are a window wiper, a door closing mechanism, rock crushers, etc. [9]. Based on the Gruebler equation for planar mechanisms, the mobility or degree-of-freedom of the mechanism is one; thus, it is a constrained mechanism fully operated by one actuator. The path generation for a four-bar linkage is a dimension-based design of the four-bar linkage lengths () and other parameters, which makes the trace point () on the coupler link follow the desire path ().
Let the coordinates of the joints O2 and O4 be () and (), respectively, so the coordinates of points B and O4 can be computed aswhere the angular positions and are positive counter-clockwise. The positions of points and are the intersection points of two circles as illustrated in Figure 1, which is calculated by a vector approach. First, let B and O4 be the centres of the circles with radii and , respectively, and let be the distance between B and O4. Then, generate , the unit vector from B to O4, and , the unit vector perpendicular to . Given that is the vector from B to one of the intersection points, the angle between and denoted by is solved using the law of cosines
The intersection points can then be obtained as
The coupler curve is formed when the crank link rotates. From Figure 1, the coupler point coordinates () in the global coordinate can be expressed aswhere is a unit vector in the direction from B to C and is a unit vector perpendicular to . The distances and can be computed using the given dimensions of , PB, and PC. Also, they can be set as design variables for optimal path synthesis. These vector forms of position analysis are used for four-bar linkage synthesis in this paper.
3. Optimization Problem and Constraint Handling
3.1. Optimization Problem
The path synthesis problem is converted to an optimization problem with an objective function expressed as the sum of square errors between the distances of and their corresponding . A set of the input angles () is assigned as design variables if the prescribed timing is not given. There are two major constraints. The first constraint is set in such a way that the generated mechanism type is a crank-rocker. This leads to two constraints: (i) sum of the shortest and longest links of linkage must be less than the sum of two remaining links and (ii) the shortest link must be an input link where one of its nearby links is set as a frame. The constraint (i) is denoted with Equation (7), while the constraint (ii) is denoted with Equation (6). In cases where the prescribed timing is not promoted, the second constraint is that the input values of must be in either ascending or descending order. The constraint (ii) is added to the design problem, which is denoted by , where . The optimization problem can then be written assubject towhere , is the number of points on prescribed or target curve, and are input link angles. and are lower and upper bounds of a design vector , respectively.
A function evaluation is carried out for the optimization problem (5) by implementing the position analysis detailed in Section 2. To enable the position analysis, the constraints (6)–(8) must be first satisfied. In the previous studies, if the conditions are not met, the objective function will be modified by adding to it a great penalty value. However, in this work, the proposed path repairing algorithm (PRA) is assigned to repair all design solutions to always be feasible before performing position analysis of a four-bar linkage.
3.2. Constraint Handling
Normally, a traditional exterior penalty function technique has been used with the constrained optimization problem (5) [1–5, 7, 8, 10–13], but it was found to be inefficient [9, 13, 14, 17]. Using such a technique is questionable when prespecifying a penalty parameter. If the parameter is too small, the resulting optimum solution may be infeasible, but if it is too large, MH may not be able to find the optimum. This leads us to propose a new strategy to deal with the constraints called the path repairing technique. The technique can be separated into two parts as repair of the input link angle constraint (8) and the Grashof criterion constraint Equations (6) and (7).
3.2.1. Input Link Angle Constraint
For an optimization without prescribed timing, a repairing technique for this phase is shown in Algorithm 1 where the input variables are . Once it is found that the values of are not in either ascending or descending orders, the design variables in the part of will be repaired. Firstly, uniform random number for is generated. The lower bound is set to be 0.0001 in order to avoid repeated values of in cases that some random numbers become zeros. Those random values are then scaled in step 2 so that the sum of them does not exceed . The first angular position of an input link is , its original value. Then, the next value of for is accumulatively added to the next input angle until the last value is obtained as . Then, the new set of input angles is in an ascending direction before returning to the position analysis of a four-bar linkage. This concept was successfully used in a sprayed plate fin heat sink design [19]. The obtained sequence of input angles always obeys the constraint (8).
|
The vector of design variables for path synthesis of a particular four-bar linkage is . Variables are the link lengths of the linkage, and are the angular positions of link , also known as timings of the crank, while other variables are shown in Figure 1. The legitimated set of the timings must obey the condition . During an optimization process, if the set of timings is decoded as, for example, = 0.5, = 0.45, = 0.67, = 1.35, = 4.50, and = 2.10, they violate constraint (8). Then, Algorithm 1 is activated to repair these values. Five (for 6 timings) uniform random numbers are generated, for example, , , , , and . The values of are then scaled down according to step 2 in Algorithm 1 to meet the condition leading to
The output modified values of are then computed as
As a result, by using Algorithm 1, the timings are always feasible. The difference between and will never exceed .
3.2.2. Grashof’s Criterion Constraint
With the same reasons as for the previous constraint, the dimensions of {} must obey the conditions (6) and (7) so that the resulting mechanism is usable. Let the bound constraints of {} in Equation (9) be separately defined as
Then, the relation can be found:
The Grashof criterion repairing technique is shown in Algorithm 2. Firstly, four uniform random numbers {} in the range of [0.0001, 1] are generated. The lower bound is set to be 0.0001 for the same reason as Algorithm 1. Then the auxiliary variables for are computed asWith this computation, it is concluded that is their minimum and is their maximum. These four values fulfil the Grashof criterion if is an input crank. Condition (7) also holds if
|
Therefore, the computational steps 3-4 in Algorithm 2 are applied. Then, the values of {} are all scaled down so that . The values of can then be computed as
Then, set . The values of are assigned to be the lengths of , respectively. With such a computing scheme, the values of {} returned from activating Algorithm 2 will always be feasible. In the search process of MH, when a function evaluation is revoked, feasibility of a design solution will be checked. If it is infeasible, Algorithms 1 and 2 will be used to repair or alter the values of and in and send them back to the main process of a metaheuristic. That means all design solutions in a population of MH are always feasible.
Given that, for example, , , , and , the Grashof criterion is not fulfilled since . These values will then be regenerated by using Algorithm 2. From step 1, given that the values of are randomly generated as , , , and , then the auxiliary variables are computed based on (14) as
Since is greater than 1.00, the values of the auxiliary variables are modified as
Thus, the new link length values are obtained using (16):which results in a crank-rocker four-bar linkage. In cases where the generated value of is greater than that of , their values are swapped. If they are equal, Step 3 in Algorithm 2 is activated ensuring that a crack-rocker is obtained after the repairing process.
4. Numerical Experiment
For evaluating the performance of the proposed path repairing technique, three path synthesis test problems of a four-bar linkage are used, whereas the optimizers are 7 established MHs. To validate the new approach, it will be compared with the exterior penalty function technique, which is traditionally applied in previous work. The path generation problems are detailed as [5, 14]:
Case 1. Path generation without prescribed timing
Design variables are
Target points are .
Limits of the design variables are as follows:
Case 2. Path generation with prescribed timing
Design variables are
Target points are .
Limits of the design variables are
Case 3. Path generation without prescribed timing
Design variables are
Target points are
Limits of the design variables are
The first synthesis problem has a straight line target path without a prescribed timing. The second test problem has a circular prescribed path with prescribed timing while the third test problem has an elliptic path without prescribed timing. The optimizers used to tackle the test problems are 7 well-known and newly developed metaheuristics. Their optimization parameter settings are given below. The population size nP = 100 is used for all algorithms, unless otherwise specified. It should be noted that the terms and variable definitions are from their original sources.(1)Artificial bee colony algorithm (ABC) [20]: the number of food sources for employed bees is set to be nP/2. A trial counter to discard a food source is 100.(2)Real-code ant colony optimization (ACOR) [21]: the parameter settings are q = 0.2 and ξ = 1.(3)Teaching-learning-based optimization (TLBO) [22]: parameter settings are not required.(4)Adaptive differential evolution with optional external archive (JADE) [23]: The parameters are self-adapted during an optimization process.(5)Population-based incremental learning (PBIL) [24]. The learning rate, mutation shift, and mutation rate are set as 0.5, 0.7, and 0.2, respectively.(6)Grey wolf optimizer (GWO) [25]: the GWO has only two main parameters to be adjusted, a and C, where a is decreased from 2 to 0 and C is randomly generated in the range of [0, 2].(7)Sine cosine algorithm (SCA) [26]. The parameter decreases linearly from a = 2 to 0.
It should be noted that the parameter settings used in this paper are either the default values provided or suggested by their authors. Also, they have been used in some previous studies with acceptable results [17, 27–30].
The total number of iterations or generations for each optimization run is set to be 500. Any MH that uses a different population size will be terminated with the total number of function evaluations as 100 × 500. Each MH runs 30 times, so that to measure its convergence rate and consistency. It should be noted that most works in the literature tended to ignore examining the search consistency of MHs. This can be carried out by running MHs many times. In this paper, each MH is run 30 times for each optimization problem. The means of the objective function values obtained from the various MHs are used as their performance indicator. Also, the comparison based on nonparametric statistical Friedman test [17, 31] is employed. Thus, this study would be a proper baseline for the topic of using MHs for four-bar linkage path synthesis in the future.
In conclusion, the path repairing technique is proposed to increase the performance in solving four-bar linkage path synthesis. To evaluate the performance of the proposed technique, three path synthesis test problems and 7 established MHs are used to study. To validate the new constraint handling approach, it is compared with the results that obtained from using the classical exterior penalty function technique, the most popularly used technique for path synthesis. Moreover, statistical results including mean, minimum, maximum, and standard deviation are also reported.
5. Design Results
The optimization of path generation of four-bar mechanisms is to find the optimized link lengths and some other parameters, which minimize the objective function. Three case studies with and without prescribed timing are considered for performance testing of the optimizers and the new constraint handling technique. The mean values of objective functions are used as the main performance indicator, where the lower mean value is more reliable optimizer. The mean objective function value determines MH reliability as MH with lower mean value is likely to be more successful in solving the problem, even with only one optimization run.
The results of Case 1 obtained from the seven optimizers with the novel path repairing technique and the penalty technique are shown in Tables 1 and 2 respectively. In the tables, the number of successful runs (a run that obtains a feasible solution), the mean objective function values from 30 optimization runs (Mean), the worse result (Max), the best result (Min), the standard deviation (Std), and the best linkage that gives the minimum objective function value of each algorithm are given. Each parameter computes based on the objective function defined as the sum of squares of the distances between the desired points, and the actual points, while the term “Error” in the tables is an average distance error of the best found solution computed aswhich has been used as a performance indicator in some previously published papers. is the number of points on the prescribed or target curve. From the results, it is seen that the proposed path repairing technique is far superior to the exterior penalty technique used in the previous studies [1–5, 7, 8, 10–13]. Based on the mean values of the objective function, all seven optimizer performances are greatly improved when implementing the path repairing algorithm as summarized in Table 3. All the optimizers can find the feasible solutions for all 30 runs. The optimizer that gives the best results is TLBO for both the mean value () and the best result (). The second and third best optimizers are ABC and GWO, respectively. The worst performer in this case, according to Mean, is SCA. In Table 2, only ABC can search for feasible solutions for all 30 runs. Based on the mean objective values, the performance of all optimizers deteriorates compared to the results from using PRA in Table 1. The best in cases of using the exterior penalty function technique is ABC, while the second best is GWO as it can search for feasible solutions for 29 from 30 runs with lower mean objective value compared to JADE. Figures 2 and 3 show the path traced by the coupler point of the best solution and the kinematic diagram of best linkage, respectively. It is found that TLBO with the novel design repairing technique gives the best result (Error = 0.0122) better than previous work [5, 14] (Error = 0.1227 (DE), and 0.0166 (DE), respectively). This design result cannot be compared with the previous work [15] due to the number of evaluation of the present work that is lower than ten times when compared with the reference. The best actual path as shown in Figure 3 is closer to the target path than in previous work [5, 14] and comparable with [15]. The four-bar linkage obtained from the best solution as shown in Figure 3 can completely rotate in an ascending direction.
For Case 2 with the number of target points at 6, the results obtained from using the various MHs with PRA and the penalty function technique are shown in Tables 4 and 5, respectively. The coupler curves and the best linkages obtained from TLBO are shown in Figures 4 and 5, respectively. In this case, all MHs are significantly improved by using PRA based on the mean objective function values. They can find feasible solutions for all runs with the given predefined number of function evaluations. Note that PBIL using PRA has higher mean value compared to that using the penalty function technique, but this is computed from only 21 successful runs when using the penalty function technique. Based on using PRA in Table 4, the best method is TLBO while the second and third best algorithms are JADE and PBIL, respectively. The worst optimizer is SCA according to the mean value. In Table 5, based on using the penalty function technique, the top three performers are JADE, GWO, and ABC in that order. In this design case, TLBO produces the best solutions for both constraint handling techniques. For this case, it is proved that using the second repairing technique (Algorithm 2) is better than the penalty function technique. The best optimizer is TLBO as with the first case which gives the best error result as 0.3073, while the results from previous work [5, 13, 14] are error = 2.3496 (DE), 2.5998 (ICA), and 1.6063 (JADE), respectively. The best actual path as shown in Figure 4 is closer with the target path than the previous work [5, 13, 14]. The four-bar linkage obtained from the best solution as shown in Figure 5 can completely rotate in an ascending direction in accordance with the prescribed timing and the circular target path.
In Case 3 with the prescribed curve with 10 points, the results are obtained by those algorithms with PRA, and the penalty function techniques are shown in Tables 6 and 7, respectively. The coupler curve and the best linkage are shown in Figures 6 and 7, respectively. From the results, it is found that TLBO gives the best results for both the mean objective function value and the best result when using PRA. Similarly to the previous two cases, all MHs are greatly improved when using PRA. All the methods can find feasible solutions for all runs with using PRA. The top three performers in cases of using PRA are TLBO, ABC, and JADE in that order while the worst method is SCA based on the mean objective values. In the results of using penalty function technique, ACOR, and SCA cannot find a feasible solution. The best method when using penalty function technique is ABC while the performance of others cannot be evaluated as they can search for feasible solutions for a few optimization runs. The best result from TLBO gives an average distance error of 0.0324 while the previous work [5, 14] had error = 1.9523 (DE), and error = 0.1641 (DE), respectively. It means that the best actual path as shown in Figure 6 is closer with the target path than the previous work [5, 14]. The linkage obtained from the best solution as shown in Figure 7 can completely rotate in an ascending direction in accordance with the elliptic target path.
Table 3 gives a summary of the comparative performance of the various metaheuristics for solving the four-bar linkage path synthesis using the new constraint handling technique and the classical exterior penalty function technique. It shows that the results from using PRA are totally superior to those obtained from using the penalty function.
It is shown that TLBO with PRA is superior to the other MHs using the penalty technique. In this study, the Friedman test and the Tukey–Kramer test are used for comparing the results. Based on the Friedman test, TLBO ranks 1st whereas the second best is JADE at value (0.02) < α (0.05) as shown in Table 8. It can be summarized that TLBO is the best performer for solving the four-bar path synthesis problems cases 1–3. Based on the Tukey–Kramer test, the mean column ranks of TLBO are significantly different from SCA.
The mean values obtained from using all optimizers of the test problems with PRA and the penalty function technique are shown in Table 3. In the statistical study, the results from ACOR and SCA are discarded because they cannot find any feasible solution in Case 3 (without using PRA). This table shows that MHs using the PRA approach give better mean than their counterparts that employ the penalty function technique for all test problems. The of average results of Friedman test are given in Table 9 which shows that MHs with the PRA technique significantly outperforms those using the penalty function technique at value (0.00016) < α (0.05).
Figures 8–10 illustrate the search histories of the best runs of TLBO for Case 1, Case 2, and Case 3, respectively. The search history compares the best runs obtained from using the penalty function technique and the path repairing algorithm. It can be seen that even though TLBO using PRA started with the worse solution (Case 3), the path repairing technique still guided the optimizer to converge considerably faster than when using the penalty function technique. All three figures show the superiority of the proposed path repairing technique for four-bar linkage path synthesis.
6. Discussion
Further discussion is provided in order to investigate the behavior of the proposed constraint handling technique PRA and why it is efficient when used with TLBO. Figure 11 shows the search history of the best runs of TLBO for the case 1 problem using PRA and the penalty function (PF) technique. The figure displays the number of iterations versus the number of infeasible solutions for both runs. The illustration also separates the number of infeasible solutions due to the timing constraint and those caused by the link length constraints. It can be seen that, with the use of PRA, infeasible solutions due to the link length constraints vanished after approximately 50 iterations. The same conclusion is applied in cases of using PF. By using PRA, the numbers of infeasible solutions due to the timing constraint disappear after around 85 iterations while the number of infeasible solutions due to the timing constraint when using PF cannot be suppressed throughout the search process. From this particular comparison, it is shown that the proposed PRA is efficient for dealing with both link length and timing constraints while the penalty function technique failed to solve the timing constraint problem.
Figure 12 shows the search history (objective function versus iterations) of the best runs of the three design cases from using TLBO in combination with PRA. In this figure, the best objective function values obtained from the teaching and learning phases of TLBO are plotted separately. It can be seen from all three path synthesis problems that the best results produced by the learning phase are equal to or better than those obtained from the teaching phase. This implies that the reproduction using the learning operator of TLBO works well with the proposed PRA. In TLBO, the teaching phase is used for exploitation while the learning phase emphasizes more on exploration. That means the proposed technique tends to be efficient with an exploration-based reproduction operator. As a result, further development of TLBO for path synthesis of a four-bar linkage can adapt from the original version by adding a self-adaptive strategy for population sizing.
7. Conclusions
This paper presents a technique to find the optimum parameters of a four-bar linkage for path synthesis using metaheuristics including ABC, JADE, PBIL, TLBO, ACOR, GWO, and SCA. The new technique, called a path repairing technique, is proposed to handle the synthesis constraints effectively. The comparative results of the three case studies show that the new path repairing technique is superior to the penalty function technique traditionally used in four-bar linkage path syntheses in the previous studies. The comparative performance of the metaheuristics shows that the TLBO with PRA is the most efficient method based on both convergence rate and consistency. The results in this work can consider as the baseline for developing an efficient optimizer for four-bar linkage path syntheses in the future. Any optimizer should be tested by running it many times to solve the synthesis problems, and the mean value of an objective function should be used as a performance indicator.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
The authors are grateful for the financial support provided by the King Mongkut’s Institute of Technology Ladkrabang, the Thailand Research Fund (RTA6180010), and the Postdoctoral Program from Research Affairs, Graduate School, Khon Kaen University (58225).