Research Article  Open Access
Luman Zhao, MyungIl Roh, "A Thrust Allocation Method for Efficient Dynamic Positioning of a Semisubmersible Drilling Rig Based on the Hybrid Optimization Algorithm", Mathematical Problems in Engineering, vol. 2015, Article ID 183705, 12 pages, 2015. https://doi.org/10.1155/2015/183705
A Thrust Allocation Method for Efficient Dynamic Positioning of a Semisubmersible Drilling Rig Based on the Hybrid Optimization Algorithm
Abstract
A thrust allocation method was proposed based on a hybrid optimization algorithm to efficiently and dynamically position a semisubmersible drilling rig. That is, the thrust allocation was optimized to produce the generalized forces and moment required while at the same time minimizing the total power consumption under the premise that forbidden zones should be taken into account. An optimization problem was mathematically formulated to provide the optimal thrust allocation by introducing the corresponding design variables, objective function, and constraints. A hybrid optimization algorithm consisting of a genetic algorithm and a sequential quadratic programming (SQP) algorithm was selected and used to solve this problem. The proposed method was evaluated by applying it to a thrust allocation problem for a semisubmersible drilling rig. The results indicate that the proposed method can be used as part of a costeffective strategy for thrust allocation of the rig.
1. Introduction
Research Background. A dynamic positioning (DP) system is especially useful to deepsea working vessel, such as drilling vessel and floating crane. The DP system can maintain the position (with a fixed location or on a predetermined track) and the heading of a working vessel exclusively through the use of active thrusters and propellers. To keep the position means to maintain the desired position within the normal excursion from the desired position and heading in the horizontal plane. For an offshore vessel, which is a semisubmersible drilling rig illustrated in Figure 1 in this case, the DP system generally consists of a position and heading reference system with sensor units, a control system, a power system, and a thruster system [1]. The position and heading reference system, along with the wind sensors and gyrocompasses, can provide the control system with information on the position of the vessel and the magnitude and direction of the environmental forces that affect the position. All of the information obtained is then used to adjust the vessel itself.
The DP system must be very responsive to changes in the conditions due to weather, water depth, and so on in order to properly control the operating parameters of the vessel. Despite the importance of the DP system to the vessel, its use in a real working environment still faces many challenges, such as a power failure and a thruster malfunction, due to its complexity. The fuel expenditure of the DP system is a considerable burden since it requires a large amount of energy. This means that the DP system is overactuated, which renders the thruster to have infinite allocation solutions. Therefore, the thrust allocation can be formulated as an optimization process under numerous constraints in order to minimize power consumption.
Related Works. In the literature, many methods have been proposed to achieve thrust allocation. Johansen and Fossen [2] provided a recent review of thrust allocation methods, but most of these have focused on proposing optimization algorithms intended for thrust allocation for industries, such as automobile, aerospace [3], and robotics, and not for ocean engineering. The relationship between the power consumption and the resulting thrust can usually be approximated by using a quadratic function, so some researchers have formulated optimization problems for thrust allocation as quadratic programming problems and have solved these with some variant of quadratic programming (QP) algorithms [4].
Some of the methods intended for ocean engineering have considered cases with ships or with offshore vessels. Wit [5] proposed a method to achieve the optimal thrust allocation for a ship. The thrust and azimuth direction of the thrusters were selected as design variables for the optimization problem, and the QP algorithm, which is a local optimization algorithm, was used to solve the problem. Similarly, Liang and Cheng [6] and Johansen et al. [7] proposed methods for the ship that include a similar formulation of an optimization problem and a solution through the use of a sequential quadratic programming (SQP) algorithm. Parikshit [8] also proposed a method to achieve the optimal thrust allocation for a drilling rig. The thrust and azimuth direction of thrusters were selected as design variables of an optimization problem, and various global optimization algorithms, such as genetic algorithm (GA) and ITHS (Intelligent Tuned Harmony Search), were used to solve this problem. Zhao et al. [9] also proposed a method that includes a similar formulation of an optimization problem and solution for a drilling rig by using GA.
In this manner, most studies in the field of ocean engineering have determined the thrust and the azimuth direction of the thrusters to be design variables for an optimization problem. In fact, the thrust can be expressed as a function of the speed of rotation, the diameter, and the thrust coefficient of a thruster. In order to reflect the specifications of the given thruster and to more easily control the thruster, the speed of rotation can be used as a design variable instead of the thrust. Meanwhile, a local optimization algorithm can be used as the optimization algorithm, such as the SQP algorithm, to derive an accurate optimum. However, this method is very sensitive to a starting point for optimization and sometimes finds a local optimum. On the other hand, a global optimization algorithm, such as the GA, does not need a starting point for the optimization but can only find a rough optimum [10].
This study proposes an optimal thrust allocation method for an offshore vessel with the aim of achieving the required generalized forces and moment for dynamic positioning while at the same time minimizing total power. Finding such a method is a particularly important issue to operate an offshore vessel while using minimal energy. The method involves mathematically formulating an optimization problem for thrust allocation and selecting and using a suitable optimization algorithm to solve the problem, not developing a new algorithm. When the optimization problem is formulated, the speed of the rotation and the azimuth direction of thrusters are selected as design variables, the total power is selected as an objective function, and the required forces and moment are selected as constraints. A hybrid optimization algorithm that incorporates global and local optimization algorithms was used for the optimization.
Table 1 shows a summary of related studies and a comparison of them with this study.

2. Thrust Allocation for Dynamic Positioning
A thruster system is important for a drilling rig to maintain its position and heading since it can simultaneously provide transverse and longitudinal thrust. An azimuth thruster consists of an electric podded drive that is fitted to the hull and can rotate 360 degrees along the horizontal angle to provide thrust in all directions. The azimuth thruster has been selected in this study because it can provide thrust in any direction in order to act as a propulsor for dynamic positioning and target optimization of a drilling rig. However, the thruster system is usually overactuated in practice, so a thrust allocation problem can be formulated into a constrained optimization problem for the system.
Figure 2 depicts a schematic diagram for the thruster arrangement of the drilling rig and the sign conventions for the forces ( and ) and moment (). In this figure, , , and represent the total forces required in the longitudinal and transverse directions and the total required moment about the vertical direction from the control system, as shown in Figure 1, respectively. Note that the positions of the thrusters are given in a 2D coordinate system (only the horizontal forces and moment are considered). The coordinate system has its origin at the center of gravity (CoG) of the drilling rig.
The environmental forces and the moment tend to move the drilling rig away from its original position during operation at sea. In order to move the drilling rig back to its original or reference position, the control system in the DP system of the drilling rig should first calculate the total forces (thrusts) and moment required. At this time, thrust allocation should be conducted to determine the thrust and azimuth direction of each thruster so that the required forces (and moment) are generated in each of the longitudinal () and transverse () and vertical () directions. The relationship between the thrust and the azimuth direction of each thruster and the total forces and moment that are required can be stated as in (1) to (3). These equations are the governing equations for the thrust allocation and can be used as equality constraints in an optimization problem:
Supposing that the drilling rig is equipped with azimuth thrusters, is the thrust from th thruster and is the azimuth direction (angle) of th thruster. and are the distances of th thruster from the CoG of the drilling rig. Equation (1) indicates that the longitudinal component of the force generated from the thrust and azimuth direction of each thruster should be equal to the total force required in the longitudinal direction. Similarly, (2) indicates that the transverse component of the force from the thrust and azimuth direction of each thruster should be equal to the total force required in the transverse direction. Equation (3) means that the vertical component of the moment generated from the thrust and azimuth direction of each thruster should be equal to the total required moment about the vertical direction.
The goal of achieving optimal thrust allocation is to guarantee efficient power use when dynamically positioning the drilling rig in the specific conditions of a given ocean environment. Meanwhile, the economic impact of minimizing power consumption is also taken into account with the ultimate aim of ensuring the stability and capability of the drilling rig during operation.
3. Thrust Allocation Method
3.1. Optimization Problem for Thrust Allocation
This study mathematically formulated an optimization problem for thrust allocation, and each component of the optimization problem is described below in detail.
3.1.1. Design Variables
Thrust allocation involves determining the thrust () and azimuth direction () of each thruster to generate the forces in longitudinal and transverse directions and the moment about the vertical direction that are required for the DP system of the drilling rig. Thus, the thrust and azimuth direction of each thruster can be used as design variables in the thrust allocation optimization problem. The objective function that is commonly used for this problem is to minimize the total power of the thrusters due to energy savings. In order to use this power as the objective function, it should be mathematically formulated as a function of the design variables. The total power will be proportional to the sum of the thrusts for each thruster. Some researchers [5–9] used the thrust and the azimuth direction of each thruster as design variables and have formulated the objective function (minimization of total power of thrusters) with the corresponding design variables. However, the rotational speed () and the azimuth direction () of each thruster can be used as design variables to reflect the characteristics of the given thruster, to more easily control the thruster, and to more delicately represent the objective function. Thus, the design variables for the thrust allocation optimization problem in this study are as follows: : speed of rotation of each thruster. : azimuth direction of each thruster.
3.1.2. Objective Function
In this study, the speed of the rotation () and the azimuth direction () of each of the thrusters were used as design variables, and thus the objective function (minimization of the total power for the thrusters) was mathematically formulated according to the following procedure where in fact the power consumed by the electric propulsion system depends on the speed of rotation and the azimuth direction of the thrusters.
The thrust and torque of a thruster can be expressed as functions of the speed of rotation ():where is the density of the sea water, is the propeller diameter of the thruster, and are the thrust and torque coefficients, respectively, and is the speed of the rotation of the thruster. The expressions for and can be obtained from the socalled open water test of the thruster. These nondimensional coefficients for thrust and torque were derived through a regression analysis that can be described according to the following parameters [11]:
and coefficients and , , , and terms can be found in a report of the Wageningen BSeries propellers [11] and were used in this study. Here, is the pitch diameter ratio, is the blade area ratio, and is the number of blades of the thruster. The coefficient for the advance is a nondimensional description of the thruster performance that can be expressed according to the speed of rotation and the speed of the advance of the thruster defined as , which can be represented by the speed near each thruster () and the azimuth direction of each thruster (). Since each of the thrusters has a different speed of rotation in practice, and of each thruster can be different:
Now, it is time to formulate the total power consumed by the thrusters of the drilling rig. The power consumption () can be formulated as
By combining (5) and (8), the power consumption of the thrusters can be rewritten aswhere is the number of thrusters.
Therefore, the objective function can be written relating the speed of rotation and the azimuth direction of each thruster. In this study, the total power (power consumption) of the thrusters is minimized by setting the objective function for the thrust allocation optimization problem as
3.1.3. Constraints
To obtain a valid thrust allocation, the governing equations shown in (11) to (13) should be satisfied. The equations represent the thrust and azimuth direction of each thruster that generate the forces required in longitudinal () and transverse () directions as well as the moment required about vertical direction () to dynamically position the vessel. If these equations cannot be satisfied, the drilling rig starts to drift. Thus, these equations are used as equality constraints for the optimization problem.
The governing equation for the longitudinal force can be stated as
Similarly, the governing equation for the transverse force can be stated as
And finally the governing equation for the yawing moment can be stated aswhere , , and are the forces required in the longitudinal and transverse directions and the moment required about the vertical direction, respectively, and are all given from the control system in the dynamic positioning system.
It is essential to minimize the thruster interaction of a semisubmersible drilling rig to ensure effective DP operation [5, 12]. In order to find a solution to provide the most effective thrust for operation, additional constraints should be considered by taking into account the physical locations of the thrusters. Thrusters that are close together have an influence on each other, and if one thruster is in the stream of another, the efficiency would drop significantly. To avoid having this happen, certain angles must be prohibited for the azimuth direction. Each thruster is assumed to have a forbidden angle of for the closest thruster in order to prevent the direct interaction of the thruster, as shown in Figure 3. With the portforward pair, the center of the angle of the outboard thruster is and for the inboard thruster it is The positive angle is measured clockwise while the negative angle is measured anticlockwise from the horizontal axis. Figure 3 shows the ATRs (Attainable Thrust Regions) for the two portforward side thrusters and (). The ATR represents a set of all physically possible surge and sway forces for th thruster.
Considering the ATR, the constraints for the azimuth direction of thrusters and () can be formulated as where represents the forbidden angle of the thrusters. and are the azimuth directions of the thrusters and (), respectively. and are the center of angle for the outboard thruster and for the inboard thruster , respectively.
3.1.4. Summary
The optimization problem for thrust allocation can be summarized as follows:where , (: number of thrusters).
Thus, this problem has one objective function, three equality constraints, and inequality constraints.
3.2. Hybrid Optimization Algorithm for Thrust Allocation
3.2.1. Overview
The optimization algorithms are generally divided into two categories: global and local. Several classes of global optimization algorithms are now available, including the genetic algorithm (GA) [13, 14] and the simulated annealing method. These algorithms are suitable for largescale problems with many local optima. However, these algorithms require further iteration to obtain an acceptable optimum and not an accurate optimum, in contrast to the local optimization algorithms. Several classes of the local optimization algorithms also exist, including sequential linear programming (SLP) [15] and sequential quadratic programming (SQP) [15] and the method of feasible directions (MFD) [16]. Each of these algorithms can effectively determine an accurate optimum. However, in some cases, these algorithms find the local optimum that is closest to the given starting point.
Various attempts have been made to combine global and local optimization algorithms in order to overcome the challenges of using them separately [10, 17–21]. Most of these studies used a combination of a global (e.g., GA) and a local optimization algorithm, which can be referred to as a hybrid optimization algorithm. In this study, the hybrid optimization algorithm was selected and used by incorporating the GA and the SQP algorithm to solve an optimization problem for thrust allocation, as described in Section 3.1. In this case, the SQP algorithm was used to improve the acceptable, global optimum obtained from the GA. That is, at the end of the GA, the SQP algorithm was executed with a starting point that corresponds to the optimum obtained from the GA. Figure 4 shows the general procedure of the hybrid optimization algorithm used in this study.
An experiment on the mathematical optimization problem was performed in order to verify the efficiency, accuracy, and applicability of the hybrid optimization algorithm that was used to optimize the thrust allocation in this study. The selected problem, which is Rastrigin’s problem, one of the benchmark problems, is being widely used to check the efficiency of optimization algorithms [22]. More details about comparative test of the hybrid optimization algorithm can be found in [10].
3.2.2. Optimization Procedure
The following optimization procedure was established to solve the optimization problem that was presented formulated above. First, the initial values were assumed for the design variables. At this time, the values can be randomly generated or can be manually set or extracted from an existing design. Now, these values are transferred to the hybrid optimization algorithm, and the values of an objective function and constraints are then calculated. We then check whether the current values of the design variables are at an optimum or not. If yes, the optimization process finishes and the result will be shown, and if not, the above steps will be repeated until the optimum is found. Figure 5 shows this optimization procedure to find an optimal thrust allocation.
4. Application of the Thrust Allocation Method
4.1. Optimization Target
The optimization target for this study is a semisubmersible drilling rig equipped with eight azimuth thrusters. Figure 5 shows the schematics of this drilling rig with a thruster system and its position. We assumed that each thruster has 4 propeller blades (), a propeller diameter () of 3.6 m, a pitch diameter ratio () of 1.1, a blade area ratio () of 0.7, and a maximum thrust of 600 KN. In addition, we assumed that the speed () near each thruster was of 1 m/s. The longitudinal and transverse positions of each thruster can be seen in Figure 6.
4.2. Mathematical Formulation
This drilling rig has eight thrusters. Thus, the design variables in this problem are as follows:where and are the speed of rotation and the azimuth direction of th thruster, respectively, and thus this problem has 16 design variables.
Now, the objective function of this problem can be stated by using (10):where is the density of sea water, is the propeller diameter of the thruster, and and are the thruster and torque coefficients, respectively. and were calculated from (6) by referring to the report on the Wageningen BSeries propellers [11].
The constraints for the required forces and moment can be stated using (11) to (13):
Equations (18) to (20) ensure that the forces generalized by the eight thrusters are equal to the forces demanded from the control system, and the moment should also be satisfied in the same way to keep the drilling rig in the desired position.
In order to overcome the energy loss due to the thrusterthruster interaction, the ATR should be considered, excluding forbidden zones, by using (14), as shown in Figure 7. For this, each thruster should satisfy the following equations simultaneously:
In addition, there are some limitations on the maximum thrust and azimuth direction for each thruster. Thus, such limitations were also considered as additional constraints in this study. The additional constraints about the maximum thrust and the azimuth direction of each thruster can be stated in the following equations, respectively:where .
Thus, this problem has 3 equality constraints and 24 inequality constraints.
4.3. Optimization Result
The problem in (16) to (23) was solved by using the hybrid optimization algorithm described in Section 3.2, and then the problem was also solved by using the GA in order to compare their relative efficiencies. We assumed that the total required forces () and moment () are given from the control system at 14 time steps, as stated in [8] and as shown in Table 2. For example, the values for the first time step are 50 KN, −600 KN, and −64,000 KN·m. For a total of 14 inputs for the total required forces and moment, the optimization was performed by using the GA and the hybrid optimization algorithm, and the optimization results are shown in Table 2. As mentioned earlier, the hybrid optimization algorithm is a combination of the GA and the SQP algorithm. That is, at the end of the GA, the SQP algorithm was executed with a starting point that corresponds to the optimum obtained from the GA. With respect to the parameters for the GA, the number of individuals in the population was 200, the maximum number of generations was 3,000, the crossover probability was 0.8, and the mutation probability was 0.01 [23]. As shown in Table 2, the hybrid optimization algorithm outperforms the GA although the GA yielded similar results for certain time steps. Figure 8 shows a convergence history of the optimization results for the first time step. This figure shows that the objective function is minimized and that the constraints are satisfied at the same time.

To evaluate the applicability of the proposed method, the optimization results of this study were compared against those obtained in Parikshit’s study [8]. The same inputs were used for the total forces and moment required from the control system for the comparative test. Parikshit used an ITHS (Intelligent Tuned Harmony Search) algorithm to solve this problem. The ITHS algorithm maintains a proper balance between the diversification and intensification throughout the search process by automatically selecting the proper pitch adjustment strategy based on its harmony memory. Table 3 shows a comparison of the optimization results between the existing study and this study. As shown in this table, the results of this study were slightly better than those of Parikshit’s study in terms of a mean value of the ratio.

5. Conclusions and Future Works
In this study, a thrust allocation method was proposed for a semisubmersible drilling rig in order to produce the generalized forces and moment that are required to dynamically position the rig while at the same time minimizing the total power consumed. First, a thrust allocation optimization problem was mathematically formulated with the corresponding design variables, objective function, and constraints. In terms of the design variables, the speed of rotation and the azimuth direction of each thruster were selected. As compared with some studies about optimal thruster allocation, the selection of the design variables can not only represent the objection function more precisely but also control the thruster more easily. The objective function aimed to minimize the total power of the thrusters, and, in terms of the constraints, the governing equations for the thrust and azimuth direction of each thruster were used to generate the forces required in longitudinal and transverse directions as well as the moment about the vertical direction to dynamically position the vessel. Some limitations were also used for each thruster, and additional constraints were introduced by considering the energy loss due to the thrusterthruster interaction and the maximum thrust and azimuth direction of each thruster. The hybrid optimization algorithm was then used to solve the formulated problem. Finally, the proposed method was applied to an example to find the optimal thrust allocation for the semisubmersible drilling rig with 8 thrusters. A comparative test was also performed as part of the current study. In this comparative test, the proposed method was observed to produce slightly better results (about 2% in terms of total power) for thrust allocation than existing study, and thus the proposed method could be used to better determine a strategy to allocate the thruster of the drilling rig.
In the future, the total forces and moment required from the control system, which was the input of this study, will be estimated by considering the current position and the equations of motion of the drilling rig. That is, a more general method to dynamically position the rig will be further studied. In addition, we will improve the present hybrid optimization method and apply the improved method to the formulated problem in the future.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work was partially supported by (a) Brain Korea 21 Plus Program (Education and Research Center for Creative Offshore Plant Engineers of Seoul National University) funded by the Ministry of Education, Republic of Korea, (b) Engineering Research Institute of Seoul National University, Republic of Korea, and (c) Research Institute of Marine Systems Engineering of Seoul National University, Republic of Korea.
References
 A. J. Sørensen, “A survey of dynamic positioning control systems,” Annual Reviews in Control, vol. 35, no. 1, pp. 123–136, 2011. View at: Publisher Site  Google Scholar
 T. A. Johansen and T. I. Fossen, “Control allocation—a survey,” Automatica, vol. 49, no. 5, pp. 1087–1103, 2013. View at: Publisher Site  Google Scholar
 C. Y. Duan, S. J. Zhang, Y. F. Zhao, and X. R. Kong, “Robust control allocation among overactuated spacecraft thrusters under ellipsoidal uncertainty,” Abstract and Applied Analysis, vol. 2014, Article ID 950127, 12 pages, 2014. View at: Publisher Site  Google Scholar
 O. Härkegård, “Dynamic control allocation using constrained quadratic programming,” Journal of Guidance, Control, and Dynamics, vol. 27, no. 6, pp. 1028–1034, 2004. View at: Publisher Site  Google Scholar
 C. D. Wit, Optimal thrust allocation methods for dynamic positioning of ships [M.S. thesis], Delft University of Technology, Delft, The Netherlands, 2009.
 C. C. Liang and W. H. Cheng, “The optimum control of thruster system for dynamically positioned vessels,” Ocean Engineering, vol. 31, no. 1, pp. 97–110, 2004. View at: Publisher Site  Google Scholar
 T. A. Johansen, T. I. Fossen, and S. P. Berge, “Constrained nonlinear control allocation with singularity avoidance using sequential quadratic programming,” IEEE Transactions on Control Systems Technology, vol. 12, no. 1, pp. 211–216, 2004. View at: Publisher Site  Google Scholar
 Y. Parikshit, Analysis, design and optimization of offshore power system network [Ph.D. thesis], National University of Singapore, Singapore, 2013.
 D.W. Zhao, F.G. Ding, J.F. Tan, Y.Q. Liu, and X.Q. Bian, “Optimal thrust allocation based GA for dynamic positioning ship,” in Proceedings of the IEEE International Conference on Mechatronics and Automation (ICMA '10), pp. 1254–1258, Xi’an, China, August 2010. View at: Publisher Site  Google Scholar
 K.Y. Lee, S. Cho, and M.I. Roh, “An efficient globallocal hybrid optimization method using design sensitivity analysis,” International Journal of Vehicle Design, vol. 28, no. 4, pp. 300–317, 2002. View at: Publisher Site  Google Scholar
 M. M. Bernitsas, D. Ray, and P. Kinley, “KT, KQ, and efficiency curves for the Wageningen Bseries propellers,” Tech. Rep., University of Michigan, Ann Arbor, Mich, USA, 1981. View at: Google Scholar
 D. T. Brown and L. Ekstrom, “Vessel thrusterthruster interactions during azimuthing operations,” in Proceedings of the 24th International Conference on Offshore Mechanics and Arctic Engineering, pp. 991–996, Halkidiki, Greece, June 2005. View at: Publisher Site  Google Scholar
 D. E. Goldberg, Genetic Algorithms in Search, Optimization, and Machine Learning, AddisonWesley, Reading, Mass, USA, 1989.
 L. Davis, Handbook of Genetic Algorithms, Van NostrandReinhold, 1991.
 J. S. Arora, Introduction to Optimum Design, Elsevier, 3rd edition, 2012.
 G. Vanderplaats, Numerical Optimization Techniques for Engineering Design, McGrawHill, 1984.
 C. Stork and T. Kusuma, “Hybrid genetic autostatics: new approach for largeamplitude statics with noisy data,” in Proceedings of the 62nd Anniversary International Meeting, pp. 1127–1131, The Society of Exploration Geophysicists, New Orleans, La, USA, 1992. View at: Google Scholar
 M. J. Porsani, P. L. Stoffa, M. K. Sen, R. K. Chunduru, and W. T. Wood, “A combined genetic and linear inversion algorithm for waveform inversion,” in Proceedings of the 63rd Anniversary International Meeting, pp. 692–695, The Society of Exploration Geophysicists, Washington, DC, USA, 1993. View at: Google Scholar
 K. Albarado, R. Hartfield, W. Hurston, and R. Jenkins, “Solid rocket motor design using hybrid optimization,” International Journal of Aerospace Engineering, vol. 2012, Article ID 987402, 9 pages, 2012. View at: Publisher Site  Google Scholar
 E. Deniz Ülker, A. Haydar, and K. Dimililer, “Application of hybrid optimization algorithm in the synthesis of linear antenna array,” Mathematical Problems in Engineering, vol. 2014, Article ID 686730, 7 pages, 2014. View at: Publisher Site  Google Scholar
 F.Z. Oujebbour, A. Habbal, R. Ellaia, and Z. Zhao, “Multicriteria shape design of a sheet contour in stamping,” Journal of Computational Design and Engineering, vol. 1, no. 3, pp. 187–193, 2014. View at: Publisher Site  Google Scholar
 H. Willi and S. Klaus, Test Examples for Nonlinear Programming Codes, Springer, Berlin, Germany, 1981.
 J. J. Grefenstette, “Optimization of control parameters for genetic algorithms,” IEEE Transactions on Systems, Man and Cybernetics, vol. 16, no. 1, pp. 122–128, 1986. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2015 Luman Zhao and MyungIl Roh. 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.