FACTA UNIVERSITATIS Series: Mechanical Engineering Vol. 18, N o 1, 2020, pp. 121 - 134 https://doi.org/10.22190/FUME191129013T © 2020 by University of Niš, Serbia | Creative Commons License: CC BY-NC-ND Original scientific paper THE SELECTION OF OPTIMAL REVERSIBLE TWO-SPEED PLANETARY GEAR TRAINS FOR MACHINE TOOL GEARBOXES Sanjin Troha 1 , ŽeljkoVrcan 1 , Dimitar Karaivanov 2 , Madina Isametova 3 1 University of Rijeka, Faculty of Engineering, Croatia 2 University of Chemical Technology and Metallurgy, Sofia, Bulgaria 3 Satbayev University, Almaty, Kazakhstan Abstract. The application of multi-criteria optimization to two-carrier two-speed planetary gear trains is outlined in this paper. In order to determine the mathematical model of multi-criteria optimization variables, the objective functions and conditions must be determined first. Two-carrier two-speed planetary gear trains with brakes on coupled shafts are analyzed in this paper. The mathematical model covers the determination of the set of the Pareto optimal solutions as well as the method for selecting an optimal solution from this set. A numerical example is provided to illustrate the procedure in which the optimal two-speed planetary gear train is selected and defined by design parameters. Key Words: Multi-criteria Optimization, Two-speed Planetary Gear Trains, Pareto Optimal Solutions, Coupled Shafts 1. INTRODUCTION Multi-criteria optimization problems are very common in many scientific and technical solutions. The optimization of gear trains as complete technical systems implies a complex mathematical model that has to describe actual system operation in actual circumstances. Planetary gear trains (PGT)s are a type of geared transmission which offers many advantages in comparison to conventional gearboxes. Therefore, the area of application of single-stage and multi-stage PGTs in mechanical engineering is increasing. Multi-stage PGTs are obtained by linking the shafts of one or two single stage PGTs. A special multi-stage PGT is a two-speed, two-carrier PGT consisting of two coupling shafts and four external shafts. This type of compound gear train has many important characteristics, the most notable being the ability to change the transmission ratio and the direction of rotation of the output Received November 29, 2019 / Accepted March 03, 2020 Corresponding author: Sanjin Troha University of Rijeka, Faculty of Engineering, Vukovarska ul. 58, 51000, Rijeka, Croatia E-mail: sanjin.troha@riteh.hr 122 S. TROHA, Ž. VRCAN, D. KARAIVANOV, M. ISAMETOVA shaft under load on demand. Therefore, they are particularly suited for applications as main drives in machinery, e.g. machine tools, cranes, etc. 1.1. State of the art The application of multi-criteria optimization to gear transmissions, especially planetary gear transmissions has not been the subject of many research papers. A population-based evolutionary multi-objective optimization approach, based on the concept of Pareto optimality, is proposed in paper [1] for the design of helical gears. Paper [2] deals with the selection of the best parameters in order to obtain the required gear quality and with the optimization of the design process itself. An analytical and computer aided procedure for the multi-criteria design optimization of multi-stage gear transmission is presented in paper [3]. The process of planetary gear transmission optimization is shown in paper [4] as a method which leads to the optimal solution.On the other hand, there were very few research efforts dealing with two-speed, two-carrier PGTs until 2003 [5]. Two-carrier PGTs consisting of two coupled and four external shafts, which enable two-speed transmissions, have significant application as gearboxes [5]. The possible schemes of these transmissions are presented in [5-10], while possible transmission structures with convenient brake layouts which could be used as two-speed transmissions are examined in [5]. A method for investigating the transmission ratio, the internal power flows and the efficiency of complex multi-carrier gearings is presented in [7]. An optimization of the two-carrier two-speed PGTs with brakes on single shafts is provided in [11]. In this example, a fishing boat transmission was chosen as input data for the numerical example of multi-criteria gearbox optimization. This paper provides an optimization of the two-carrier two-speed PGTs with brakes on the coupled shafts, in continuation of the optimal selection choice methodology application. The characteristics of a machine tool transmission have been used as input data for the numerical example of multi-criteria optimization application. Apart from the determination of the set of the Pareto optimal solutions, the weighted coefficient method was applied in order to determine the optimal solution. 2. MATHEMATICAL MODEL FOR PLANETARY GEAR TRAIN OPTIMIZATION The two-carrier two-speed PGTs with brakes on coupled shafts are built from basic types of PGT. The basic type of PGT (type 2k-h, variant A) is an arrangement using a central sun gear with external gearing (1), external ring gear with internal gearing (3), planet gears with external gearing(2) and planet carrier (h), as shown in Fig.1. The planets are in simultaneous mesh with the sun gear and the ring gear. Also, a Wolf- Arnaudov’s symbol can be used, indicating the torque on the main elements as a function of basic transmission ratio i0. The equations for the basic transmission ratio and the ideal torque ratio calculation are also pointed out in Fig. 1. The carrier shaft is the summary element of the basic PGT, as a negative transmission ratio is obtained by stopping the planet carrier, indicating a change of the direction of rotation of the output element. The Choice of Optimal Reversible Two-Speed Planetary Gear Train for Machine Tools Gearboxes 123 Basic transmission ratio: 0 130 zzi . Ideal torque ratio: 3 3 0 1 1 | | 1 T z t i T z . Fig. 1 Basic type of PGT and Wolf-Arnaudov’s symbol with torque ratios (1 – sun gear; 2 – planet; 3 – ring gear; h – planet carrier) The process of finding the optimal solution starts with the definition of a mathematical model, as stated in [12]. The complete mathematical model of the basic type of PGT was described in the aforementioned paper and a brief summary will be also presented in this section. It is necessary to define variables, objective functions and functional constraints in order to define a mathematical model. 2.1. Variables The following variables are considered by this model: the number of teeth of sun gear z1, the number of teeth of planet gears z2, the number of teeth of ring gear z3, the number of planets w n , gear module n m and gear face width b. The optimization variables are of the mixed type: the gear tooth numbers are positive and negative integers, the number of planets is a discrete value, the module is a discrete standard value (acc. to ISO 54), while the face width is a continuous variable. The gear tooth numbers and the number of planets are non-dimensional values while the module and the face width are given in millimeters. 2.2. Objective functions The characteristics used by the model to determine the objective functions are the volume, mass, efficiency and manufacturing cost of PGTs. The volume of the PGTs is used as an overall dimension expression, and the gears are approximated with a cylinder volume with the diameter equal to the pitch diameter and the height equal to the face width. The fact that the planets are inside the ring gear makes it possible for the PGT volume to be expressed by Eq. (1) 2 23 3 cos cos cos4 wt tn zm bV (1) where t is the transverse pressure angle, wt23 is the working transverse pressure angle for the pair 2-3 and β is the helix angle at the pitch diameter. 124 S. TROHA, Ž. VRCAN, D. KARAIVANOV, M. ISAMETOVA Since the mass of a particular gear is determined as gear volume multiplied by the density of gear material and fact that mass is determined as the sum of all gear masses in a PGT, this criterion has been expressed as Eq. (2): 23 2 2 2 33 12 2 2 2 22 12 2 2 2 112 2 cos cos cos cos cos cos cos 25.0 wt t wt t w wt tn zkzknzk m bm (2) Efficiency is one of the most important criteria for the design and evaluation of the arrangement quality. The calculation of the gear transmission efficiency is generally confined to losses depending on the friction on tooth flanks in contact while neglecting the losses in bearings and losses due to oil viscosity, i.e. restricted to the calculation of contact power losses [12-14]. The model, followed by the developed computer program, is adjusted to the most commonly used variant with the sun gear as the input element, and the carrier as the output element while the ring gear remains stationary. Basic PGT efficiency in this case is given by Eq. (3) [12]: 0 00 1 1 i i (3) where 0 is the efficiency with the planet carrier stationary, as expressed by Eq. (4) 32113 3 0 20.035.015.0 1 zzzzz z (4) The economic demands must be also taken into consideration in the techno- economical optimization, as these demands are directly related to production costs. The time needed for the manufacture of gears is taken as a measure of the production costs and as an economic factor. This function is then determined as a sum of the time periods needed for the manufacturing of sun gear (TP1), planets (TP2) and ring gear (TP3), i.e. 321 TTnTF wT (5) The production times are determined according to Fette, Lorenc and Höfler [15]. 2.3. Functional constraints The functional constraints are the conditions required for the proper operation of a system. There are numerous exceptions that need to be taken into consideration for PGTs to operate correctly in comparison to conventional gear transmissions. The exceptions presented in this model are related to assembly conditions, geometrical conditions and strength conditions.The assembly conditions include the conditions of coaxiality, adjacency and conjunction [16]. The geometrical conditions are related to the undercutting and profile interference, the ratio of the pressure angle to the working transverse pressure angle, the tooth thickness and the tooth space width, the transverse contact ratio value, the sliding speeds at the point of contact, the ratio of the pinion face width to the pinion reference diameter, etc. These conditions have been ensured in accordance with the actual standards (ISO TC 60 list of Standards 090915). The strength conditions, safety factors for bending strength and surface durability of each gear, are checked according to ISO 6336-1 to ISO 6336-3 [17]. The Choice of Optimal Reversible Two-Speed Planetary Gear Train for Machine Tools Gearboxes 125 2.4. Steps in the optimization process The optimization process begins by generating all solutions for the assigned input data. All 6-tuples of design parameters (z1, z2, z3, nw, mn, b) satisfying the functional constraints are generated for the given input data (transmission ratio, input number of revolution, input torque, service life in hours, application factor, accuracy grade Q(DIN 3961), and the values of the objective functions for every 6-tuple are computed. These 6- tuples form a set of feasible solutions. An optimal solution is then selected, based on the established objective functions and constraints, and determined by variables. The mathematical model of nonlinear multi-criteria problem can be formulated as follows: 1 2max ( ), ( ), , ( ) subject to k f x f x f x x S (6) Here, f1(x),..., fk(x) are objective functions, x = (x1,...,xn) is the vector of decision variables and S is the set of feasible solutions. Every point x S is mapped to the point (f1(x), f2(x),..., fk(x)) in k dimensional objective space. Therefore, one can observe the objective set: 1 2 {(( ( ), ( ), , ( ) | ) k F f x f x f x x S (7) The notation „maxˮ determines a simultaneous maximization of all the objective functions. If any objective function has to be minimized, the minimization of function fi(x) is performed by maximization of function fi(x). According to the structure of feasible set S, discrete multi-criteria optimization problems do exist. In this PGTs problem, six decision variables exist, corresponding to the basic design parameters: x = (x1,x2,x3,x4,x5,x6) = (z1,z2,z3,nw,mn,b). Furthermore, there are four objective functions: volume V(x), mass m(x), efficiency (x) and production costs T(x): 1 2 3 4 ( ) ( ), ( ) ( ), ( ) ( ), ( ) ( ) p f x V x f x m x f x x f x T x (8) Therefore, the mathematical model of nonlinear multi-criteria problem in concrete task, can be formulated as follows: 1 2 3 4 max{ ( ), ( ), ( ), ( )} subject to f x f x f x f x x S (9) As multi-criteria optimization problems are mathematically ill-defined which can be seen from the definition, a criterion for selecting the optimal solution must be defined. The most important criterion for selecting these „equally goodˮ solutions is the Pareto optimality concept: The solution x S is Pareto optimal if no solution y S exists which maintains fi(x) fi(y) for all i = 1,...,n and maintains strict inequality, i.e. fi(x) < fi(y) for at least one index i. Determination of the Pareto optimal solutions set is the first stage in optimal solution finding. The optimal solution is selected in the next stage, where the weighted coefficients method is applied to select the optimal solution from the Pareto solutions set. 126 S. TROHA, Ž. VRCAN, D. KARAIVANOV, M. ISAMETOVA 2.5. Weighted coefficients method For this method, the following scalarized problem must be set up: 0 0 1 1 max ( ) ( ) ( ) . . M m m f x w f x w f x s t x S (10) Weighted coefficients (or weights) iw are positive real numbers and 0 0 1 ( ) ( ) ( ) i i i f x f f x are normalized objective functions where 0 i f are normalizing coefficients [12]. All solutions obtained by using this method are Pareto optimal [12]. This model may be used regardless of the existence of priority functions or not [15]. The complete optimization procedure is implemented in the PLANGEARS software. 3. TWO-SPEED TWO-CARRIER PLANETARY GEAR TRAINS 3.1. Two-carrier planetary gear trains structures and labeling method In cases where two-speed transmissions are required, a mechanism obtained by connecting two basic PGTs shown in Fig. 1 is one of the best suited design solutions. By joining two shafts of one PGT with two shafts of another PGT a mechanism is formed with four external shafts in total, Fig. 2. Fig. 2 Symbolic representation of a compound planetary gear train with four external shafts The two component trains can be joined in in 12 different ways, resulting in a PGT with four external shafts [18]. An alphanumerical label (S11…S56) is attached to each of the 12 structural schemes, indicating the ways of connection between the shafts of the main elements of both component trains (Fig. 3). In every presented scheme it is also possible to place the brakes as well as the driving and the driven machine on external shafts in 12 different ways (V1…V12), corresponding to layout variants (Fig. 4). The compound trains in consideration can be classified into three different groups according to whether the brakes are placed on the coupled shafts, on the single shafts or both on the coupled and the single shaft. The Choice of Optimal Reversible Two-Speed Planetary Gear Train for Machine Tools Gearboxes 127 Fig. 3 Systematization of all schemes of two-carrier planetary gear train with four external shafts Fig. 4 Systematization of all layout variants (A-input shaft, B-output shaft) 3.2. Operations of planetary gear trains with different layout variants By placing the brakes on two shafts, a braking system is obtained in which the alternating activation of the brakes shifts the power flow through the PGT, causing a change of the transmission ratio. Some PGTs of this type are described in [5,7,18,19,20]. The possible power flow paths for PGTs are analyzed, and functions of the transmission ratio for some trains of this type are deduced in [5,18,19]. 15 kinematic schemes of the considered type are presented in [6], and achievable values of transmission ratios and efficiencies are given. A computer program DVOBRZ for the selection of an optimal variant of similar multi-speed PGTs is described in [5,18,19], and charts of shifting capabilities for all possible two-speed PGTs are given in [5]. Each variant has its own characteristics that determine the possibilities of transmission ratio changes. Some variants can be presumed to work in both transmission ratios as reducers and multipliers, while other variants work like a reducer with one ration and like a multiplier with the other ratio. Also, some variants change the direction of rotation when after a transmission ratio change, while other variants keep the direction of rotation after changing the transmission ratio. The transmission ratio of each PGT stage depends only on its basic transmission ratio (ideal torque ratio). The compound train with brakes on the coupled shafts (layout variant V1 in the Fig. 4) with power flows when some of the brakes are active is symbolically shown in Fig. 5 by means of a Wolf-Arnaudov symbol. The green dotted line represents the power of relative motion. There are two possible directions: from the sun gear to the ring gear or from the ring gear to the sun gear. Also, the expressions for transmission ratios and efficiency obtained by using torque method when brake Br1 is activated are given on the left side, while the expressions for brake Br2 activated are given on the right side [5]. 128 S. TROHA, Ž. VRCAN, D. KARAIVANOV, M. ISAMETOVA 1 (1 ) 1 I I B II I I Br I II A II II t t T t t t i t t T t t ; 2 (1 ) (1 )(1 ) 1 B I II I II Br I II A T t t t t i t t T 0 0 0with losses 1 without losses 1 ( ) ( ) (1 ) I II I II II IIB Br IB II II t t tT tT t t ; with losses 0 0 2 without losses ( ) (1 )(1 ) ( ) (1 )(1 ) B I I II II Br B I II T t t T t t Fig. 5 Power flows on the Wolf-Arnaudov’s symbol through the train with brakes on the coupled shafts Regardless of which brake is applied, both the component trains operate actively, as seen in Fig. 5. The power input and output are on the single shafts. Also, the direction of the power flow is the same in both variants. When the upper brake (Br1) is applied, the input element is the sun gear of the first stage. The power is transmitted through the ring gear of the first stage and the ring gear of the second stage to the output element - carrier of the second stage. In the other case, when the lower brake (Br2) is applied, the path from the input element (sun gear of the first stage) to the output element (carrier of the second stage) includes the carrier of the first stage and the sun gear of the second stage. As this transmission changes the direction of rotation with the transmission ratio, it is suitable for application in the machine tools which have a working motion with considerable load at low speed and a return motion to the initial position at high speed and light load. The Choice of Optimal Reversible Two-Speed Planetary Gear Train for Machine Tools Gearboxes 129 4. RESULTS AND DISCUSSION Variant S15V1(Figs. 6 and 7) was chosen to demonstrate the procedure of multi-criteria optimization application. A symbolic review of the transmission composition with kinematic scheme is shown in Fig. 6, and the path of the power flow through the transmission is shown in Fig. 7. a) b) Fig. 6 Symbolic review of transmission composition (a), kinematic scheme (b) a) b) Fig. 7 Power flow through the transmission on the conceptual scheme with brake Br1 applied (a) and with brake Br2 applied (b) The type of transmission is selected according to the transmission requirements of the machine tool concerned. The necessary transmission ratios are: iBr1 = 6 in one direction (with Br1activated) and iBr2 = 40 in the other (with Br2 activated). In this case, the ideal torque ratios for both the planetary gear stages can be defined from the shifting capabilities diagram, Fig. 8 [5]. 130 S. TROHA, Ž. VRCAN, D. KARAIVANOV, M. ISAMETOVA Fig. 8 Shifting capabilities diagram of compound trains S15V1 The transmission ratios are defined as functions of ideal torque ratios, 1 ( / ) Br I II i t t (1 ) II t and 1 (1 )(1 ) Br I II i t t , so the ideal torque ratios are defined by using equations from Fig. 8: 09.5It and 57.5IIt [5]. 4.1. The first compound gear train stage (I) With brake Br1 applied and brake Br2 turned off, the carrier of the first stage and sun gear of the second stage are immovable. The first stage is determined in this mode as the torque of the first stage is greater in this mode. The input of the system (A) is the sun gear of the first stage (Figs. 6 and 7a). The output of the system (B) is the planet carrier of the second stage (Figs. 6 and 7). The power is transmitted to the ring gear of the first stage and then to the ring gear of the second stage and finally to the carrier of the second stage. The input data required for the multi-criteria optimization application is: i0 = 5.09, nin = 2850min 1 , Tin = 33.5Nm (P = 10kW), L = 8000 h, KA = 1.25, IT7 for all gears, material z1/material z2/material z3= 20MoCr4/20MoCr4/34CrNiMo6, SH min = 1.1, SF min = 1.2, i = 3%, z1 = 1530. The feasible set consists of 834 solutions, from which 43 Pareto solutions are selected. By application of the weighted coefficient method with weighted coefficients: w1 = 0.5, w2 = 0.0, w3 = 0.0, w4 = 0.5 the solution shown in Table 1 is obtained, with a set of objective functions values shown in Table 2. The weighted coefficients are chosen in accordance with techno-economic optimization requirements. The Choice of Optimal Reversible Two-Speed Planetary Gear Train for Machine Tools Gearboxes 131 Table 1 Optimal solution of the first stage in the first case Variable values x1 = z1 x2 = z2 x3 = z3 x4 = nw x5 = mn x6 = b 15 31 -78 3 2 16 Table 2 Objective functions for solution shown in Table 1 f1 in mm 3 f2 in kg f3 f4 in min 305815.19 1.4529 0.9865 88.523 By prioritizing other objective functions, i.e. by choosing other values for weighted coefficients, other solutions would become optimal. The differences between solutions obtained in this way are logical. For example, the efficiency can be increased only with a large number of teeth considering Eq. (3) and the fact that efficiency is the only function that has to be maximized. 4.2. The second compound gear train stage (II) With brake Br2 applied and brake Br1 turned off, the ring gears of both the stages are immovable. The second stage is determined in this mode as the torque at the sun gear of the second stage will be greater in this mode. As in the previous mode, the input of the system (A) is the sun gear of the first stage and the output of the system (B) is the carrier of the second stage (Figs. 6 and 7b). The power is transmitted to the carrier of the first stage and then to the sun gear of the second stage and finally to the carrier of the second stage. Since the ideal torque ratio of this stage is tII = 5.57, the basic transmission ratio is i0 = t = 5.57 and the transmission ratio of the second stage in this mode is i = 1 i0 = 6.57. The input parameters of the second stage are equal to the output parameters of the first stage. Because of that, with brake Br2 applied, the ring gear of the first stage is also immovable, and the corresponding input number of revolutions at the second stage is calculated by means of the already defined first stage: 1- min677.459 2.6 2850 i n nn outIIinII The torque on the sun gear of the second stage is calculated using the following equation: 1 (1 ) 33.5 (1 5.2) 207.7 Nm II A I T T t The other input data required for the multi-criteria optimization application in stage is equal to the one in the first stage: L = 800 h, KA = 1.25, IT7 for all gears, material z1/material z2/material z3=20MoCr4/20MoCr4/34CrNiMo6, SH min = 1.1, SF min = 1.2, i = 3%, z1 = 1530. The feasible set consists of 1778 solutions, from which 43 Pareto solutions have been isolated. By application of the weighted coefficient method with weighted coefficients: w1 = 0.5, w2 = 0.0, w3 = 0.0, w4 = 0.5 the solution shown in Table 3 was obtained, using a set of objective function values shown in Table 4. 132 S. TROHA, Ž. VRCAN, D. KARAIVANOV, M. ISAMETOVA Table 3 Optimal solution of the second stage Variable values x1 = z1 x2 = z2 x3 = z3 x4 = nw x5 = mn x6 = b 15 32 -81 3 2.75 27 Table 4 Objective function for solution shown in Table 3 f1 in mm 3 f2 in kg f3 f4 in min 1040497.69 4.985 0.9865 121.85 By comparing the design parameters of both the stages it can be concluded that the design parameters enable a compact transmission design that is very important for installation in a machine tool. Also, the Pareto optimality concept as the criterion for selecting an equally good solution can be applied to compound PGT according to these criteria. Furthermore, the weighted coefficient method can be used for selecting the optimal solution from a Pareto set, and it can be adjusted to varying impacts of individual criteria functions. Since the numbers of teeth of all gear are known, it is now possible to determine the realized transmission ratios and efficiency with brake Br1 applied and with brake Br2 applied, according to equations given in Fig. 5. Ideal torque ratio in the first stage: 3 1 | | | 78 | 5.2 15 I z t z . Ideal torque ratio in the second stage: 3 1 | | | 81 | 5.4 15 II z t z . Transmission ratio with Br1 activated: 1 5.2 (1 ) (1 5.4) 6.16 5.4 I Br II II t i t t . Transmission ratio with Br2 activated: 2 (1 ) (1 ) (1 5.2) (1 5.4) 39.68 Br I II i t t . The deviations of the actual and required transmission ratios are in the permissible range. Basic efficiency of the first stage: 3 0 3 1 1 2 3 0.15 0.35 0.20 ( 78) 0.15 0.35 0.20 1 1 0.984 ( 78) 15 15 31 ( 78) I z z z z z z . Basic efficiency of the second stage: 3 0 3 1 1 2 3 0.15 0.35 0.20 ( 81) 0.15 0.35 0.20 1 1 0.98442 ( 81) 15 15 32 ( 81) II z z z z z z . The Choice of Optimal Reversible Two-Speed Planetary Gear Train for Machine Tools Gearboxes 133 Total efficiency with Br1 applied: 0 0 0 1 1 0.9816 (1 ) I II I II II II Br I II II t t t t t t . Total efficiency with Br2 applied: 0 0 2 (1 )(1 ) 0.97361 (1 )(1 ) I I II II Br I II t t t t . It can be noticed that both efficiencies are very high (97…98%), and that the efficiency with brake Br1 on is slightly higher than the efficiency with brake Br2 on. The second step can be performed by including the efficiency determined by defined gear tooth numbers. The calculation of input torque is shown here: 1 0 (1 ) 33.5 (1 0.984 5.2) 204.91 Nm II A I I T T t . By applying the procedure of the optimal solution where only the value of the input torque has been changed, the solution shown in Table 5 is obtained using the objective function shown in Table 6. The feasible set consists of 1767 solutions, from which 44 Pareto solutions can be deducted. The weighted coefficients have the following values: w1 = 0.5, w2 = 0.0, w3 = 0.0, w4 = 0.5. Table 5 Optimal solution of the second stage in the second step Variable values x1 = z1 x2 = z2 x3 = z3 x4 = nw x5 = mn x6 = b 15 32 -81 3 2.75 26 Table 6 Objective function for solution shown in Table 5 f1 in mm 3 f2 in kg f3 f4 in min 1001960.7375 4.80 0.984 119.611 A difference is noticed in one variable only – face width. As the input torque is slightly lower, the expected face width will be smaller too. The distinction is negligible; therefore, it is not necessary to carry out the procedure including the efficiency determined by a defined number of teeth of the first stage. 5. CONCLUSION An original method that combines two computer programs (DVOBRZ and PLANGEARS) for multi-criteria optimization of two-carrier two-speed PGTs with brakes on coupled shafts has been presented in this paper. These compound gear trains consist of two basic type of PGTs and have considerable application in systems which need different transmission ratios and direction changes (e.g. as machine tool gearboxes which work with 134 S. TROHA, Ž. VRCAN, D. KARAIVANOV, M. ISAMETOVA a considerably greater transmission ratio in one direction and direction changing with a smaller transmission ratio in the other). The same procedure was successfully implemented in the optimal solution choice of two-carrier two-speed PGTs with brakes on single shafts, leading to a universal method of compound PGT optimization. The optimal solution is determined by considering design parameters, such as mass and production cost as objective functions, and by using multi-criteria optimization and the weight coefficient method for choosing the optimal solution from the Pareto optimal solution set. This approach can be successfully used for the basic PGT type and compound gear trains assembled from basic types, as shown in this paper. The results obtained using this procedure are in accordance with the literature on technical system optimization and indicate a good choice of optimization methods. Furthermore, this approach indicates a possibility for application to other PGT types. REFERENCES 1. Tudose, L., Buiga, O., Jucan, D., Stefanache, C., 2008, Multi-objective optimization in helical gears design, The Fifth International Symposium about Design in Mechanical Engineering-KOD 2008, Novi Sad, Serbia, pp. 77-84. 2. Tkachev, A., Goldfarb, V., 2009, The concept of optimal design for spur and helical gears, The 3 rd International Conference Power Transmissions '09, October 2009, Chalkidiki, Greece, pp.59-62. 3. Rosić, B., 2001, Multicriterion optimization of multistage gear train transmission, Facta Universitatis- Series Mechanical Engineering, 1(8), 2001, pp. 1107-1115. 4. Brüser, P., Grüschow, G., 1989, Otimierung von Planetengetrieben,Antiebstechnik, 2, pp. 64-67. 5. Troha, S., 2011, Analysis of a planetary change gear train’s variants, (in Croatian), PhD Thesis, University of Rijeka, Engineering Faculty, Rijeka, Croatia. 6. Kudrjavtsev, V. N., Kirdyiashev, I. N., 1977, Planetary gears, (in Russian), Handbook, Mashinostroenie, Leningrad. 7. Arnaudow, K., Karaivanov, D., 2005, Systematik, Eigenschaften und Möglichkeiten von zusammengesetzten Mehrsteg-Planetengetrieben, Antriebstechnik, 5, pp. 58-65. 8. Ivanov, A.N., 1990, Evaluation of diametric dimensions of planetary gearboxes in the design phase, (in Russian),Vestnikmashinostroenie, 7, pp. 16-19. 9. Lechner, G., Naunheimer, H., 1999, Automotive Transmissions, Springer-Verlag, Heidelberg. 10. Jelaska, D., 2012, Gears and Gear Drives, University of Split, Croatia. 11. Stefanović-Marinović, J., Troha, S., Milovančević, M., 2017, An application of multicriteria optimization to the two-carrier two-speed planetary gear trains, Facta Universitatis-Series Mechanical Engineering, 15(1), pp. 85-95. 12. Stefanović-Marinović, J., Petković, M., Stanimirović, I., Milovančević, M., 2011, A model of planetary gear multicriteria optimization, Transactions of FAMENA, 35(4), pp. 21-34. 13. Sriatih, A., Yedukondalu, G., Jagadeesh, A., 2011, Mechanical efficiency of planetary gear trains: An Estimate, Mechanical Engineering Research, 1(1), pp. 97-102. 14. Del Castillo, J.M., 2002, Theanalytical expression of the efficiency of planetary gear trains, Mechanism and Machine Theory, 37, pp. 197-214. 15. Stefanović-Marinović, J., 2008, Multicriterionoptimization of planetary power transmission gear pairs, (in Serbian), PhD Thesis, University of Niš, Faculty of Mechanical Engineering, Niš, Serbia. 16. Niemann, G., Winter, H., 1989, Maschinenelemente, Band II, zweite völlig neubearbeite Auflage, Springer-Verlag Berlin. 17. International Organization for Standardization, 2006, Calculation of load capacity of spur and helical gears, International standard ISO 6336-2. 18. Troha, S., Žigulić, R., Karaivanov, D., 2014, Kinematic operating modes of two-speed two-carrier planetary gear trains with four external shafts, Transactions of FAMENA, 38(1), pp. 63-76. 19. Troha, S., Lovrin, N., Milovančević, M., 2012, Selection of the two–carrier shifting planetary gear train controlled by clutches and brakes, Transactions of FAMENA, 36(3), pp. 1-12. 20. Troha, S., Petrov, P., Karaivanov, D., 2009, Regarding the optimization of coupled two carrier planetary gears with two coupled and four external shafts, Machine building and electrical engineering, 1, pp. 49–56.