Mathematical Problems in Engineering

Volume 2013 (2013), Article ID 419043, 7 pages

http://dx.doi.org/10.1155/2013/419043

Research Article

Design Optimization of a Speed Reducer Using Deterministic Techniques

Ming-Hua Lin,¹ Jung-Fa Tsai,² Nian-Ze Hu,³ and Shu-Chuan Chang2,4

¹Department of Information Technology and Management, Shih Chien University, No. 70, Dazhi Street, Taipei 10462, Taiwan

²Department of Business Management, National Taipei University of Technology, No. 1, Section 3, Chung-Hsiao E. Road, Taipei 10608, Taiwan

³Department of Information Management, No. 64, Wunhua Road, Huwei Township, Yunlin County 632, Taiwan

⁴Department of Information Management, St. John’s University, No. 499, Section 4, Tam King Road, Tamsui District, New Taipei City 25135, Taiwan

Received 26 June 2013; Accepted 10 September 2013

Academic Editor: Yi-Chung Hu

Copyright © 2013 Ming-Hua Lin et al. 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.

Abstract.

The optimal design problem of minimizing the total weight of a speed reducer under constraints is a generalized geometric programming problem. Since the metaheuristic approaches cannot guarantee to find the global optimum of a generalized geometric programming problem, this paper applies an efficient deterministic approach to globally solve speed reducer design problems. The original problem is converted by variable transformations and piecewise linearization techniques. The reformulated problem is a convex mixed-integer nonlinear programming problem solvable to reach an approximate global solution within an acceptable error. Experiment results from solving a practical speed reducer design problem indicate that this study obtains a better solution comparing with the other existing methods.

1. Introduction

Many engineering design problems are formulated as mathematical programming models. In last few decades, these nonlinear engineering problems have been investigated in much research that solved the formulated problems by different methods. The methods can be generally categorized into metaheuristic and deterministic approaches. To compare the performance of different optimization algorithms, several structural engineering applications are often solved to validate or test the suitability of the optimization algorithms. The speed reducer problem is one of the benchmark problems in structural optimization. The problem represents the design of a simple gear box used in a light airplane between the engine and propeller to allow each to rotate at its most efficient speed.

A large number of algorithms have been developed to solve different engineering optimization problems. In order to overcome the computational drawbacks of existing numerical methods, many metaheuristic algorithms that combine rules and randomness to imitate natural phenomena [1] have been developed. The most general metaheuristic methods include evolutionary computation (EC), tabu search (TS), simulated annealing (SA), ant colony optimization (ACO), and particle swarm (PS) [2]. The surveys of applications and algorithmic advances for metaheuristic algorithms are provided by Glover and Kochenberger [3], Lee and Geem [1], and Bianchi et al. [4]. Li and Papalambros [5] used the global optimization knowledge that is incorporated in several types of rules concerning constraint activity, redundancy, and dominance to solve the speed reducer problem. Ku et al. [6] solved the speed reducer problem using the Taguchi method that emphasizes the design of a robust product insensitive to disturbances. Akhtar et al. [7] developed an optimization algorithm based on a sociobehavioural concept of society and civilization to solve the same problem. The essence of their methodology is derived from the concept that the behaviour of an individual changes and improves due to social interaction with the society leaders. Rao and Xiong [8] proposed a hybrid genetic algorithm that combines the advantages of random search and deterministic search methods to improve the convergence speed and computational efficiency for solving mixed-discrete nonlinear design optimization problems. They also solved the speed reducer problem to demonstrate the effectiveness and robustness of their approach. Cagnina et al. [9] developed a particle swarm optimization algorithm to solve constrained engineering optimization problems and used four standard engineering design problems including a speed reducer problem to validate their algorithm. Jaberipour and Khorram [2] proposed two new harmony search (HS) metaheuristic algorithms for engineering optimization problems with continuous design variables and applied their method to solve the speed reducer problem.

Although the metaheuristic algorithms have the advantages of broad applicability, easy implementation, and robustness, these methods cannot guarantee global optimality of the solution. Several deterministic approaches based on mathematical programming techniques have been developed to solve engineering design problems. Tosserams et al. [10] proposed a decomposed problem formulation based on the augmented Lagrangian penalty function and the block coordinate descent algorithm for quasiseparable multidisciplinary design optimization problems. They solved the speed reducer design problem by the proposed decomposition algorithms. Lu and Kim [11] proposed a decomposition algorithm for the multidisciplinary design optimization problems with complementarity constraints based on the regularization technique and inexact penalty decomposition. They also solved the design problem of a speed reducer. One major deterministic approach to globally solve generalized geometric programming problems is to reformulate the original problems into convex mixed-integer nonlinear programming (MINLP) problems. Some transformation techniques have been developed to convexify the nonconvex functions. Pörn et al. [12] introduced different convexification strategies to deal with posynomial and negative binomial terms. Floudas and Pardalos [13] and Maranas and Floudas [14] proposed exponential transformations to treat nonconvex terms. Lundell et al. [15] proposed some transformation techniques to solve optimization problems including signomial functions to global optimality. Li and Lu [16] applied convexification strategies and piecewise linearization techniques to solve generalized geometric programming problems with free discrete/continuous variables. Tsai and Lin [17] proposed an efficient method to solve a posynomial geometric programming problem with separable functions by applying an appropriate variable transformation and an efficient piecewise linearization formulation. Lin et al. [18] used convexification strategies and piecewise linearization techniques to solve engineering optimization problems including the speed reducer design problem. Lu [19] proposed a convexification transformation method (beta method) based on the concept of 1-convex functions to improve the efficiency of solving generalized geometric programming problems. Huang [20] proposed a deterministic optimization approach to solve geometric programming problems including the speed reducer design problem. His method converts all signomial terms into convex and concave terms, and then the concave terms are further treated with a piecewise linearization method. Lin and Tsai [21] also used convexification strategies and piecewise linearization techniques to solve the speed reducer design problem. This study applies an efficient optimization approach to globally solve speed reducer design problems based on deterministic techniques. In addition to convexification strategies and piecewise linearization techniques, this study applies optimization-based range reduction techniques [22, 23] to improve computational efficiency in globally solving the speed reducer design problem. Compared with existing methods, the proposed method is capable of obtaining a better solution.

The rest of the paper is organized as follows. Section 2 describes the process of globally solving a speed reducer design problem. A practical speed reducer problem is solved in Section 3 to demonstrate the effectiveness of the proposed method. After that, conclusion remarks are made in Section 4.

2. Global Optimization Approach of a Speed Reducer Design Problem

A speed reducer is part of the gear box of mechanical system, and it is used in many other types of applications. The design of the speed reducer is a more challenging benchmark, because it involves seven design variables [25]. As shown in Figure 1, the design of the speed reducer is considered with the face width (), the module of the teeth (), the number of teeth on pinion (), the length of the first shaft between bearings (), the length of the second shaft between bearings (), diameter of the first shaft (), and the diameter of the second shaft (). Another schematic of the speed reducer is presented in Figure 2 with its design variables being labeled.

Figure 1: Speed reducer design [24].

Figure 2: A schematic of the speed reducer [11].

This problem is taken from Golinski [26]. The objective is to minimize the total weight of the speed reducer while satisfying eleven constraints. The constraints include the limits on the bending stress of the gear teeth, surface stress, transverse deflections of shafts 1 and 2 due to transmitted force, and stresses in shafts 1 and 2. The mathematical programming model of a speed reducer problem considered in this study is expressed as follows.

The original speed reducer problem described previously can be simplified as the following generalized geometric programming problem.

The simplified problem above is a nonconvex program. Based on the deterministic techniques, this study transforms the problem into a convex MINLP problem by the convexification strategies and piecewise linearization methods. Then the reformulated problem can be solved by convex MINLP solvers to obtain a global optimal solution.

First, we determine that certain classes of signomial terms in the above simplified problem are convex and do not necessitate any transformations. Consequently, the number of concave functions requiring to be piecewise linearized decreases, and the resulting problem is a computationally efficient model. and are convex terms. According to Maranas and Floudas [27], and are also convex terms. Then, the nonconvex monomials are transformed. By taking exponential transformation [28, 29] on the variable with a positive exponent, the positive monomial terms , , , , , , , , and are transformed into convex terms , , , , , , , , and , respectively, where , . By taking power transformations [28–32] on the variables to make the sum of the exponents not greater than one, the negative monomial terms , , , , and are transformed into convex terms , , , , and , respectively, where , , , .

The nonconvex problem can be convexified and underestimated by the convexification strategies mentioned previously if the inverse transformations (, , , , and , ) are approximated by piecewise linear functions. The efficiency of the piecewise linearization technique has a critical impact on the computational efficiency in solving the reformulated problems. Vielma and Nemhauser [33] proposed a linearization approach that has favorable tightness properties. Their experimental results showed that the Vielma and Nemhauser [33] method significantly outperforms other models. Tsai and Lin [17] employed the Vielma and Nemhauser [33] method to solve posynomial geometric programming problems. This study also adopts the Vielma and Nemhauser [33] method to linearly approximate the inverse transformations. Compared with the Lin and Tsai [21] method, this study utilizes range reduction techniques [22, 23] to further improve the computational efficiency.

Adding more break points can construct a tighter underestimator of the original problem, and the obtained solution is more closer to the real global solution. If is the th constraint and is the solution derived from the reformulated model, then the number of break points does not need to increase until , where is the feasibility tolerance.

3. Computational Experiments

By using the deterministic approach introduced above, this study reformulates the original speed reducer problem as a convex MINLP problem as follows.

In the transformation process, the piecewise linearization technique introduced by Vielma and Nemhauser [33] is utilized to approximate (), (), and (). The Vielma and Nemhauser [33] method represents a piecewise linear function with break points by binary variables. By using binary variables, respectively, we convert this program to a convex MINLP problem with break points, respectively, used in linearly approximating the inverse transformations. The reformulated problems are solved by LINGO [34]. Table 1 lists the reported solutions from LINGO, objective values on the reported solutions, and errors in constraint under different numbers of break points. Using more break points derives a solution with a lower error in constraint. Figure 3 indicates the objective value obtained from the proposed method under different numbers of break points. We observe that the objective value approximates the real global objective value better as the number of break points increases. Figure 4 indicates the CPU time required to solve the speed reducer problem under different numbers of break points. The required CPU time to solve the reformulated model tends to grow exponentially as the number of break points becomes large. To enhance computational efficiency, this study applies the range reduction techniques to effectively tighten variable bounds. Table 2 lists variable bound, solution, objective value, accumulated CPU time, and error in constraint in each iteration to solve this problem by the proposed method with range reduction. Thirty-two break points are used in the piecewise linearization process in each iteration. The accumulated CPU time consists of the CPU time to update variable bounds and solve the reformulated models iteratively. The global solution (3.5, 0.7, 17.0, 7.3, 7.7153190, 3.350282, and 5.286654) with objective 2994.471921 and an error in constraint below 10⁻⁶ can be obtained within 8 minutes. If no range reduction is adopted and 512 line segments are used in the piecewise linearization process, the global solution with an error in constraint below is obtained within 16 minutes as shown in Table 1.

Table 1: Experiment results of the speed reducer problem by the proposed method under different numbers of break points.

No. of break points
No. of binary variables
Reported solution (
CPU time (mm:ss)
Objective value
Error in constraint

8
3
(3.5, 0.7, 17.0, 7.3, 7.714826, 3.347398, 5.286205)
00:07
2993.457535
0.002526473

16
4
(3.5, 0.7, 17.0, 7.3, 7.715248, 3.349748, 5.286589)
00:09
2994.308994
0.000418000

32
5
(3.5, 0.7, 17.0, 7.3, 7.715291, 3.350039, 5.286628)
00:10
2994.408856
0.000157319

64
6
(3.5, 0.7, 17.0, 7.3, 7.715313, 3.350187, 5.286648)
00:42
2994.459757
0.000024774

128
7
(3.5, 0.7, 17.0, 7.3, 7.715318, 3.353125, 5.286653)
01:27
2994.467375
0.000009551

256
8
(3.5, 0.7, 17.0, 7.3, 7.715320, 3.350213, 5.286654)
05:23
2994.470348
0.000001492

512
9
(3.5, 0.7, 17.0, 7.3, 7.715320, 3.350214, 5.286654)
15:09
2994.470603
0.000000596

Table 2: Experiment results of the speed reducer problem by the proposed method with range reduction.

Iteration
Variable bound
Reported solution (
Accumulated CPU time (mm:ss)
Objective value
Error in constraint

1

(3.5, 0.7, 17.0, 7.3, 7.715291, 3.350039, 5.286628)
00:10
2994.408856
0.000157319

2

(3.5, 0.7, 17.0, 7.3, 7.7153190, 3.350282, 5.286654)
07:56
2994.471921
0.000000264

Figure 3: Objective values of the speed reducer problem under different numbers of break points.

Figure 4: CPU time required to solve the speed reducer problem under different numbers of break points.

Table 3 displays the comparison of results with existing methods. The solutions listed in the table are reported from their original research, and the objective values are computed from the objective function = − + + on the reported solutions. Tosserams et al. [10] and Lu and Kim [11] obtained optimal solutions with higher errors and higher objective values than those of the proposed method. Although Ku et al. [6], Lin et al. [18], and Huang [20] obtained optimal solutions with lower objective values than those of the proposed method, the errors in constraint in these three methods are higher than those in the proposed method. Compared with the solutions of Akhtar et al. [7], Rao and Xiong [8], Cagnina et al. [9], Jaberipour and Khorram [2], and Li and Papalambros [5], the proposed method results in a lower objective value under the same feasibility tolerance .

Table 3: Comparisons of optimal solutions of the speed reducer problem by different methods.

Methods
Method type
Reported solution (
Objective value
Error in constraint

Ku et al. [6]
Metaheuristic
(3.6, 0.7, 17, 7.3, 7.8, 3.4, 5.0)
2876.219475
<10⁰

Akhtar et al. [7]
Metaheuristic
(3.506122, 0.700006, 17, 7.549126, 7.859330, 3.365576, 5.289773)
3008.197440
<10⁻⁶

Rao and Xiong [8]
Metaheuristic
(3.5, 0.7, 17, 7.3, 7.8, 3.36, 5.29)
3000.959715
<10⁻⁶

Cagnina et al. [9]
Metaheuristic
(3.5, 0.7, 17, 7.3, 7.8, 3.350214, 5.286683)
2996.347849
<10⁻⁶

Jaberipour and Khorram [2]
Metaheuristic
(3.5, 0.7, 17, 7.3, 7.71533233833903, 3.35021510925684, 5.28666403545462)
2994.477531
<10⁻⁶

Li and Papalambros [5]
Metaheuristic
(3.5, 0.7, 17, 7.299999, 7.715317, 3.350541, 5.286654)
2994.553869
<10⁻⁶

Tosserams et al. [10]
Deterministic
(3.5, 0.7, 17, 7.3, 7.72, 3.35, 5.29)
2996.645783
<10⁻³

Lu and Kim [11]
Deterministic
(3.5, 0.7, 17, 7.3, 7.670396, 3.542421, 5.245814)
3019.583365
<10⁻¹

Lin et al. [18]
Deterministic
(3.5, 0.7, 17, 7.3, 7.715, 3.349, 5.286)
2993.738921
<10⁻²

Huang [20]
Deterministic
(3.495652, 0.7000002, 17, 7.30000007, 7.7120386, 3.343372, 5.285352)
2990.124384
<10⁻²

Proposed method
Deterministic
(3.5, 0.7, 17.0, 7.3, 7.7153190, 3.350282, 5.286654)
2994.471921
<10⁻⁶

4. Conclusions

The speed reducer problem is used in many other types of applications and is one of the standard benchmark problems in structural optimization. Although many metaheuristic algorithms have been developed to solve this problem, these methods cannot guarantee global optimality of the solution. This study applies a deterministic approach based on convexification strategies and piecewise linearization methods to globally solve a speed reducer design problem. However, numerous break points are utilized in the linearization process for reaching an approximate global solution with a low error, and much CPU time is required to solve the reformulated model. Therefore, this study also adopts the range reduction technique to enhance the computational efficiency. Compared with metaheuristic methods, this study guarantees the global optimality of the solution. Compared with other deterministic methods, this study obtains a better solution with a lower error in constraint.

Acknowledgment

The research is supported by the Taiwan NSC Grants NSC 101-2410-H-158-002-MY2 and NSC 102-2410-H-027-012-MY3, and Shih Chien University Grant 102-05-04002.

References

K. S. Lee and Z. W. Geem, “A new meta-heuristic algorithm for continuous engineering optimization: harmony search theory and practice,” Computer Methods in Applied Mechanics and Engineering, vol. 194, no. 36–38, pp. 3902–3933, 2005.

M. Jaberipour and E. Khorram, “Two improved harmony search algorithms for solving engineering optimization problems,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 11, pp. 3316–3331, 2010.

F. W. Glover and G. A. Kochenberger, Handbook of Metaheuristics, International Series in Operations Research & Management Science, Kluwer Academic Publishers, Boston, Mass, USA, 2003.

L. Bianchi, M. Dorigo, L. M. Gambardella, and W. J. Gutjahr, “A survey on metaheuristics for stochastic combinatorial optimization,” Natural Computing, vol. 8, no. 2, pp. 239–287, 2009.

H. L. Li and P. Papalambros, “A production system for use of global optimization knowledge,” Journal of Mechanisms, Transmissions, and Automation in Design, vol. 107, no. 2, pp. 277–284, 1985.

K. J. Ku, S. S. Rao, and L. Chen, “Taguchi-aided search method for design optimization of engineering systems,” Engineering Optimization, vol. 30, no. 1, pp. 1–23, 1998.

S. Akhtar, K. Tai, and T. Ray, “A socio-behavioral simulation model for engineering design optimization,” Engineering Optimization, vol. 34, no. 4, pp. 341–354, 2002.

S. S. Rao and Y. Xiong, “A hybrid genetic algorithm for mixed-discrete design optimization,” Journal of Mechanical Design, vol. 127, no. 6, pp. 1100–1112, 2005.

L. C. Cagnina, S. C. Esquivel, and C. A. Coello Coello, “Solving engineering optimization problems with the simple constrained particle swarm optimizer,” Informatica, vol. 32, no. 3, pp. 319–326, 2008.

S. Tosserams, L. F. P. Etman, and J. E. Rooda, “An augmented Lagrangian decomposition method for quasi-separable problems in MDO,” Structural and Multidisciplinary Optimization, vol. 34, no. 3, pp. 211–227, 2007.

S. Lu and H. M. Kim, “A regularized inexact penalty decomposition algorithm for multidisciplinary design optimization problems with complementarity constraints,” Journal of Mechanical Design, vol. 132, no. 4, Article ID 041005, 12 pages, 2010.

R. Pörn, I. Harjunkoski, and T. Westerlund, “Convexification of different classes of non-convex MINLP problems,” Computers and Chemical Engineering, vol. 23, no. 3, pp. 439–448, 1999.

C. A. Floudas and P. M. Pardalos, State of the Art in Global Optimization: Computational Methods and Applications, vol. 7 of Nonconvex Optimization and Its Applications, Kluwer Academic Publishers, Dordrecht, The Netherlands, 1996.

C. D. Maranas and C. A. Floudas, “Global optimization in generalized geometric programming,” Computers and Chemical Engineering, vol. 21, no. 4, pp. 351–369, 1997.

A. Lundell, J. Westerlund, and T. Westerlund, “Some transformation techniques with applications in global optimization,” Journal of Global Optimization, vol. 43, no. 2-3, pp. 391–405, 2009.

H.-L. Li and H.-C. Lu, “Global optimization for generalized geometric programs with mixed free-sign variables,” Operations Research, vol. 57, no. 3, pp. 701–713, 2009.

J.-F. Tsai and M.-H. Lin, “An efficient global approach for posynomial geometric programming problems,” INFORMS Journal on Computing, vol. 23, no. 3, pp. 483–492, 2011.

M.-H. Lin, J.-F. Tsai, and P.-C. Wang, “Solving engineering optimization problems by a deterministic global optimization approach,” Applied Mathematics & Information Sciences, vol. 6, no. 3, supplement, pp. 1101–1107, 2012.

H.-C. Lu, “An efficient convexification method for solving generalized geometric problems,” Journal of Industrial and Management Optimization, vol. 8, no. 2, pp. 429–455, 2012.

C. H. Huang, “Engineering design by geometric programming,” Mathematical Problems in Engineering, vol. 2013, Article ID 568098, 8 pages, 2013.

M. H. Lin and J. F. Tsai, “Optimal design of a speed reducer,” Applied Mechanics and Materials, vol. 376, pp. 327–330, 2013.

M.-H. Lin and J.-F. Tsai, “Range reduction techniques for improving computational efficiency in global optimization of signomial geometric programming problems,” European Journal of Operational Research, vol. 216, no. 1, pp. 17–25, 2012.

C. S. Adjiman, I. P. Androulakis, and C. A. Floudas, “Global optimization of mixedinteger nonlinear problems,” AIChE Journal, vol. 46, no. 9, pp. 1769–1797, 2000.

R. Hassan, B. Cohanim, O. de Weck, and G. Venter, “A comparison of particle swarm optimization and the genetic algorithm,” in Proceedings of the 1st AIAA Multidisciplinary Design Optimization Specialist Conference, pp. 18–21, April 2005.

A. H. Gandomi and X.-S. Yang, “Benchmark problems in structural optimization,” in Computational Optimization, Methods and Algorithms, vol. 356 of Studies in Computational Intelligence, pp. 259–281, Springer, Berlin, Germany, 2011.

J. Golinski, “Optimal synthesis problems solved by means of nonlinear programming and random methods,” Journal of Mechanisms, vol. 5, no. 3, pp. 287–309, 1970.

C. D. Maranas and C. A. Floudas, “Finding all solutions of nonlinearly constrained systems of equations,” Journal of Global Optimization, vol. 7, no. 2, pp. 143–182, 1995.

R. Pörn, K.-M. Björk, and T. Westerlund, “Global solution of optimization problems with signomial parts,” Discrete Optimization, vol. 5, no. 1, pp. 108–120, 2008.

A. Lundell and T. Westerlund, “Convex underestimation strategies for signomial functions,” Optimization Methods and Software, vol. 24, no. 4-5, pp. 505–522, 2009.

T. Westerlund, “Some transformation techniques in global optimization,” in Global Optimization: From Theory to Implementation, L. Liberti and N. Maculan, Eds., vol. 84, pp. 45–74, Springer, New York, NY, USA, 2007.

J.-F. Tsai and M.-H. Lin, “Global optimization of signomial mixed-integer nonlinear programming problems with free variables,” Journal of Global Optimization, vol. 42, no. 1, pp. 39–49, 2008.

A. Lundell, J. Westerlund, and T. Westerlund, “Some transformation techniques with applications in global optimization,” Journal of Global Optimization, vol. 43, no. 2-3, pp. 391–405, 2009.

J. P. Vielma and G. L. Nemhauser, “Modeling disjunctive constraints with a logarithmic number of binary variables and constraints,” Mathematical Programming, vol. 128, no. 1-2, pp. 49–72, 2011.

LINGO, Release. 11.0, Lindo System Inc., Chicago, Ill, USA, 2004.