Kumar Gaurav, Rajiv Tiwari and Twinkle Mandawat

(Department of Mechanical Engineering, Indian Institute of Technology Guwahati, Guwahati 781039, India)

Abstract: Optimization of cylindrical roller bearings (CRBs) has been performed using a robust design. It ensures that the changes in the objective function, even in the case of variations in design variables during manufacturing, have a minimum possible value and do not exceed the upper limit of a desired range of percentage variation. Also, it checks the feasibility of design outcome in presence of manufacturing tolerances in design variables. For any rolling element bearing, a long life indicates a satisfactory performance. In the present study, the dynamic load carrying capacity Cd, which relates to fatigue life, has been optimized using the robust design. In roller bearings, boundary dimensions (i.e., bearing outer diameter, bore diameter and width) are standard. Hence, the performance is mainly affected by the internal dimensions and not the bearing boundary dimensions mentioned formerly. In spite of this, besides internal dimensions and their tolerances, the tolerances in boundary dimensions have also been taken into consideration for the robust optimization. The problem has been solved with the elitist non-dominating sorting genetic algorithm (NSGA-II). Finally, for the visualization and to ensure manufacturability of CRB using obtained values, radial dimensions drawing of one of the optimized CRB has been made. To check the robustness of obtained design after optimization, a sensitivity analysis has also been carried out to find out how much the variation in the objective function will be in case of variation in optimized value of design variables. Optimized bearings have been found to have improved life as compared with standard ones.

Keywords: cylindrical roller bearing; optimization; robust design; elitist non-dominating sorting genetic algorithm (NSGA-II); fatigue life; dynamic load carrying capacity

0 Introduction

Cylindrical roller bearings(CRBs) are used for fast operations because of their low-friction torque characteristics. They have the capacity to sustain large radial loads and a few thrust loads at moderate speeds and are being used in several industries like the aerospace, automobiles, and marines. If equipped with few lubrication holes, a special retainer and a specially contoured raceway and roller profiles, it can operate satisfactorily even under the high-speed and high-load conditions for applications in aerospace industry. There are many different causes of rolling bearing failure, for example, the fatigue, wear, running marks, corrosion/erosion, over-hitting, and overload. However, fatigue is often considered the classic failure mode of rolling bearings, which is related with the dynamic capacity of the bearing. To have a minimum variation of the objective function (e.g., fatigue life of the bearing), in case of presence of manufacturing tolerances in the design variables, the optimization must be a robust design. The robust design also ensures that the design is still feasible, i.e., the dimensions obtained after robust optimization are such that even if there is variation in their value during manufacturing, the resulting variation in the objective function will be within a desirable range. In the present study, a robust design optimization method has been used, which is based on the flexible physical programming approach. Flexibility of robust design optimization means that the designer is allowed to specify the various levels of each design variable as well as of the objective function.

Messac and Ismail Yahaya[1]came forth with a methodology for the robust design optimization, the methodology being dependent on a flexible physical programming approach and ensuring that the objective function varies by minimum possible value in case of uncertainties or variations in design variables as per designer’s choice. Chen et al.[2]examined how effective the physical programming approach is in ensuring that the design outcome does turn out to be a robust one and established that the physical programming approach is superior over the other conventional methods. Giassi et al.[3]used robust design techniques to develop a multi-disciplinary optimization method for use in concurrent engineering. The Taguchi’s design[4]model was mainly based on the loss function. With this model, the quality indeed got improved, but there were several inefficiencies in the method, and for problems with non-linearity of higher orders, it would often give less than satisfactory results. Hence, for the present study, the robust design has been achieved using flexible programming approach, which treats each objective function and its desired tolerance, and the design variables and their maximum possible tolerances independently, while expresses designer preference, hence eliminating all the above stated deficiencies.

Changsen[5]examined the rolling bearings to explain their load distribution, deformation, kinematics, geometry, and stresses, hence calculating the dynamic load rating, static load rating, life and Elasto-hydrodynamic lubricant film thickness. Harris[6]analyzed the rolling element bearings comprehensively. Seireg[7]studied the use of different methods of optimization in designing the mechanical devices, and also analyzed a few of its illustrative examples. Deb et al.[8]developed the NSGA-II algorithm, which is implemented in the current study. NSGA-II algorithm possesses less computational complexity, gives better convergence and better spread for Pareto optimal fronts as compared with its contemporaries. Branke et al.[9]studied multi-objective evolutionary algorithm, but in the modified form, which can concentrate the search on knee regions. Hirani[10]optimized the hydrodynamic bearings to minimize the power loss, temperature rise, and oil feed flow.

Link et al.[11]performed optimization on a suction muffler and bearing to reduce the bearing power losses and enhance the suction muffler energy efficiency using the genetic algorithm (GA). Saruhan et al.[12]developed a method for optimum design of tilting pad bearings, which is based on GA. Hirani and Suh[13]performed the optimization of fluid-film steadily loaded journal bearings to minimize the power loss and oil flow. Wang[14]studied the usefulness of parallel optimization for the lubrication analysis with the help of GA. Kim et al.[15]developed a method for identification of bearing parameters, which was based on the scheme of global optimization making use of measured unbalance response of rotor-bearing system. Gupta et al.[16]performed single and multi-objective optimization of ball bearings taking the dynamic capacity, static capacity, and elasto-hydrodynamic minimum film thickness as the objective functions. The study also deals with sensitivity analysis of various design variables to observe their impact on the objective functions. Xie et al.[17]did a study on grease lubricated ball bearing’s Elasto-hydrodynamic lubrication and performed a numerical calculation on it. Tiwari et al.[18]performed the optimization of dynamic capacity for a tapered roller bearing and then carried out a sensitivity analysis to find out the effect of design variables on dynamic capacity. Solanki et al.[19]performed bearing design optimization to maximize its load carrying capacity with respect to clearance, and to minimize the film thickness and temperature rise. Panda et al.[20]took help of the particle-swarm optimization (PSO) algorithm for a radial ball bearing to maximize its fatigue life. Dragoni[21]performed optimization of tapered roller bearings to maximize its static load capacity taking the internal dimensions as design variables. Tiwari and Chandran[22]performed the multi-objective optimization of CRBs taking the dynamic capacity, Elasto-hydrodynamic minimum film thickness, and the maximum temperature as objective functions. The study also deals with comparison of life of the customary and optimized bearings. Dewangan et al.[23]performed bearing optimization to minimize the power loss using GA. Kim et al.[24]utilized a micro-genetic algorithm to optimize an angular contact ball bearing, thereby enhancing the values of its radial and axial stiffness. Tudose et al.[25]did a study to find out the optimal bench bearing arrangement in order to maximize its life. Kalyan and Tiwari[26]performed the multi-objective optimization of needle roller bearings taking dynamic capacity and elasto-hydrodynamic minimum film thickness as the objective functions. The study also deals with sensitivity analysis to find out the variations in the two objective functions with specific variation in design variables. Verma and Tiwari[27]performed robust single-objective optimization of tapered roller bearing taking dynamic capacity as objective function. It gives a robust optimal design by taking variability into account. Salunkhe et al.[28]utilized different optimization techniques to find out the usefulness of the modified design of rolling contact bearing. In addition to the above mentioned literature, standard texts[29-30]have also been referred to to get a basic understanding of bearings and optimization algorithms. Hardly any study has yet been done on the robust optimization of rolling bearings, as there are always some manufacturing tolerances in the bearing design drawing. With these objectives in mind, the robust design optimization of CRBs is done using the NSGA-II algorithm in this study.

In the present paper, a constrained non-linear robust optimization problem for the design of CRB has been proposed and a detailed description of the procedure to solve it is described. The problem consists of nine design variables, four of which are basic design variables and the remaining are constrain parameters. The basic design variables are basic geometries of the bearing that directly affect the dynamic capacity. The five constraint parameters are used to restrict the feasible design space. The dynamic capacity is optimized, while due to the tolerances in the bearing dimensions, its variation has been minimized and the tolerances have been taken from the bearing catalogues which help in finding the upper and lower bounds of the design variables. The main work has been done on the mathematical problem formulation of the robust optimization. Constraints based on geometry and strength criterion have been formulated. A total of nineteen constraints for deterministic optimization and twenty-one constraints for robust optimization have been included. The design problem has been optimized by real-coded GA and the dynamic capacity obtained after optimization is more than that available in the standard catalogue.

Hence, in the present section, in addition to the brief introduction, a summary of the related works done previously on bearing design has also been stated and motivation for the present work is derived. In Section 1, the geometrical aspects of CRBs, on which the present study is focused, have been explained in brief. Besides, mathematical formulations have been derived based on studies done previously in similar fields, i.e., the robust optimization, and the single and multi-objective optimizations of rolling-element bearings. The section has been further divided into sub-sections that deal individually with the formulations of objective functions, constraints, and variable bounds required for the optimization algorithm. In Section 2, implementation of NSGA-II to rolling element bearing using different parameters in the optimization algorithm NSGA-II have been discussed. In Section 3, the results obtained have been discussed quantitatively and qualitatively, and in Section 4 conclusions have been made regarding the nature of obtained results.

1 Geometrical Aspects and Problems Formulation for Robust Optimization of CRBs

CRB is identified using standard boundary dimensions, namely, the outer diameter (D), bore diameter (d), and bearing width (B).To avoid sharp edges, chamfering is provided on the faces of both raceways. As shown in Fig.1, the chamfering dimensions for the roller bearing are the outer ring chamfer heightr1, the outer ring chamfer widthr2, the inner ring chamfer heightr3, and the inner ring chamfer widthr4.The minimum value of these chamfering radii is provided in the bearing catalogues. The internal geometry of CRB is described well using the following: roller diameterDr, pitch diameterDm, number of rollersZ, and effective roller lengthle.For the present problem, the boundary dimensions and minimum chamfer radius are consulted from the SKF general catalogue[31].

In the current study, a robust optimization (which is non-linear in nature) will be performed for CRBs, for which an evolutionary algorithm NSGA-II has been chosen because of its elite-preserving strategy, explicit diversity-preserving mechanism, and overall ability to preserve a better solution spread. First, a brief on solving a deterministic optimization is provided.

A deterministic optimization problem comprises 3 parts: objectives, constraints, and variable bounds.

1)Objectives:

Minimize/Maximize

f(x)

(1)

Fig.1 Typical CRB cross-section

wherexrefers to the design variable vector to be optimized andfis the objective function. Here, the design variable vectorxconsists of the variables to be used for the deterministic optimization, which in case of CRB is given as

x=[DmDrZleKDminKDmaxεeβ]T

(2)

whereKDmaxis the maximum roller diameter limiter,KDminis the minimum roller diameter limiter,εis the parameter for the outer ring strength consideration,eis the parameter for mobility condition andβis the parameter for the effective length of the roller. The parametersDm,Dr,Zandlerepresent the bearing internal geometry. The remaining five parameters give the bearing a feasible design space in the form of constraints to make the dimensions obtained from optimization possible for manufacture. The range of values needs to be input for each of these nine variables, out of which a set of variables is obtained that gives the optimum objective function. The range for each of these variables will be discussed in Section 1.3.

2)Constraints:

Subject to

gj(x,p)≥0,j=1,2,……,J

(3)

whereprefers to vector of those parameters in the constraints that are not present in the objective function:

p=[DdB]T

(4)

whose value remains constant during optimization, but may undergo variations during manufacturing, hence their variation must be considered, which will be discussed in the robust optimization;grefers to the design constraints for deterministic optimization;jis the number of design constraints used in deterministic optimization.

The values for these parameters will be discussed in Section 1.3.

3)Variable Bounds:

xl(L)≤xl≤xl(U),l=1,2,……,L

(5)

Decision variablexlis thel-th element of the design variable vectorx, and will take a value within the lower and upper boundsxl(L)andxl(U), respectively; superscript (L) refers to lower limit and superscript (U) refers to upper limit. The procedure to solve robust optimization is explained as follows.

As has been explained in the work by Messac and Ismail Yahaya[1]on the robust design using physical programming approach, certain modifications need to be made in the objective function, constraints, and variable bounds in order to make the design robust. The new objective functions, constraints, and design space are formulated as follows. Hence, the physical programming problem model can be expressed as

a)Objective function:

(6)

The only change from Eq.(1) in the case of deterministic optimization is the inclusion of term Δf, which comprises terms from the vector Δx; although the design variable vector to be used for optimization is stillxonly. It must be noted that Δ may be the variation in design variable, constraint or objective function. The vector Δxcomprises

Δx=[ΔDmΔDrΔle]T

(7)

where ΔDmis the bearing pitch diameter tolerance, ΔDris the roller diameter tolerance, and Δleis the roller length tolerance. The value of each of these tolerances has to be input by the user. Decision of the value for each of these tolerances will be discussed in Section 1.3.

b)Subject to the constraints

gj(x,p)≥Δgj(x,p,Δx,Δp)

j=1,2,……,J

(8)

The change from Eq.(3) in the case of deterministic optimization is the inclusion of term Δg, which comprises of terms from the vector Δxand Δp.The vector Δpcomprises

Δp=[ΔDΔdΔB]T

(9)

where ΔDis the bearing outer diameter tolerance, Δdis the bearing bore diameter tolerance, and ΔBis the bearing width tolerance. Here the value for each of these tolerances must be input by the user as well. Decision of the value for each of these tolerances is discussed in Section 1.3.

The constraints from Eq. (9) can be further stated as

hk(x,p)=gk(x,p)-Δgk(x,p, Δx, Δp)

k=1,2,...,K

(10)

whereKis the number of robust design constraints used in robust optimization. Eq.(10) must be written in the inequality form as

hk(x,p)≥0

(11)

wherehis the design constraint in inequality form. It must be noted that the number of deterministic optimization constraints is different from the number of robust design constraints. Hence, they must be represented by different values of “J” and “K”, respectively. Eqs.(8) and (10) show the conversion of each deterministic optimization constraint to the robust optimization constraint. For Eq. (8), only the deterministic constraint is utilized, hence “j” is used for this equation. While for Eq. (10), it is already converted to robust design constraint, hence it is shown by “k”, wherekis the number of robust design constraints used. But it must be noted that for all those constraints that are getting converted from deterministic to robust, the values of “J” and “K” are the same, although here the values differ for “j” and “k” due to the addition of two new constraints in the robust optimization prescribing the bounds of variation in the objective function. The presence of these two extra constraints in case of the robust optimization shall be explained in details in Section 1.2.

The bounds for design variables undergoing manufacturing variations will change from those in Eq.(5) and the new variable bounds are given as

xl(L)+Δxl≤xl≤xl(U)-Δxl,l=1,2,……,L

(12)

whereLis the number of design variablesx.To find out the variation in objective function, as seen in Eq.(6), it can be calculated as

(13)

Similarly, the changes made in constraints of deterministic optimization to convert it to a robust optimization as seen in Eq.(8) are as follows:

(14)

whereMis the number of parametersp.

1.1 Objective Functions

The function to maximize/minimize is designed here. The dynamic capacityCdis optimized, which is defined as “constant radial load that can be endured by a group of evidently similar bearings for a rating life of one million revolutions of the inner ring for a static load and fixed outer ring”. The expression for dynamic capacity of CRB is mathematically defined as[6]

(15)

with

and

(16)

wherebm(=1.36) is a component inCdto take into account the latest developments in technology,λ(=0.61) is the reduction component for the discrepancies that may occur during the manufacturing and mounting of roller bearing,ν(=1.36) is a component to justify the edge loading andiis number of rows of roller[32]. For the fatigue life calculation, an equivalent radial load is obtained based on actual applied radial and axial loads on the bearing, and the dynamic capacity obtained from Eq. (15) is used in conjunction with the equivalent radial load for the life estimation.

max[f(x)]=max[Cd]=min[-Cd]

(17)

Since the robust optimization is performed, the objective is modified according to Eq. (6) as

min{f+Δf},f=-Cd

(18)

As is already known, the variation in the objective function will only be due to the design variables that do have a tolerance and are present in the mathematical formula for the objective function. Therefore, in case ofCd, the variation is not due to the bearing outer diameterD, bore diameterdor widthB.Using Eq. (13) and (18), there is

(19)

Here, only positive sign is taken so as to find out the magnitude of variation in objective function to reduce this magnitude to a minimum possible value using Eq.(18) and make it lie within a desirable range with the constraints in the form of Eq.(11). Detailed explanation of how to obtain

has been given in Appendix A.

Hence, with Eqs. (20)-(22), the revised values of objective function due to tolerances in design variables can be obtained. ThenCdcan be maximized and the variation in its magnitude can be minimized with NSGA-II algorithm[8, 28].

(20)

(21)

(22)

1.2 Constraints

The constraints affecting the optimization are defined here. Conversion of the constraints from deterministic optimization to robust optimization is also discussed. Constraints help in defining a feasible domain of values for design parameters and variation in the objective function. For the deterministic optimization[26], 19 constraints are taken. But since the robust optimization is performed, two more constraints are in the present optimization problem. They have been explained as follows. The complete list of constraints used in deterministic as well as robust optimization is given in Table B1 of Appendix B. All the constraints are explained for deterministic optimization.

1)Constraints 1 and 2:

The pitch diameter of the bearing must lie between the bearing inner diameter and bearing outer diameter, thereby ensuring enough chamfering space for the inner and outer rings’ sharp corners, as given in Table B1 of Appendix B. The values ofr1minandr3minare given as input. Decision of their values is discussed in Section 1.3.

2)Constraints 3 and 4:

The range of values that can be taken for mean diameter of the roller is based on strength as well as geometric considerations, while its lower limit came with the aid of contact stress expression, and the upper limit came with bearing geometry. It must have a minimum value such that the roller can endure the contact stress arising due to external loads. The contact stresses have been explained in details in Constraint 14. The maximum value for rollerDris limited by boundary dimensions, i.e., bearing outer diameterDand inner diameterd; and it can be calculated by taking into consideration the radii for chamfering of the inner and outer ring corners, so that chamfered corners can be obtained. Therefore, from Ref. [6], the range of values forDrcan be taken as given in Table B1 of Appendix B.Qmaxis the maximum normal load acting between the roller and the raceway, andσCmaxis the maximum contact stress on the roller.Qmaxcan be found out using the Stribeck’s expression[6], which is given in Table B1 of Appendix B.Fris the radial load acting on the roller,α(= 00) is the angle between the cylindrical roller and the outer raceway. For Constraint 3, the maximum contact stressσCmaxis taken as the safe contact stressσCsafe.The values for the parametersσCmaxandσCsafeneed to be input. Decision of their values andFris discussed in Section 1.3.

3)Constraints 5 and 6:

For the number of rollers, the lower limit is calculated from the bearing pitch diameter’s lower limit and the roller mean diameter’s upper limit. Similarly, the upper limit is calculated from the bearing pitch diameter’s upper limit and roller mean diameter’s lower limit. Hence,as indicated in Ref. [6], the range of values forZis listed in Table B1 of Appendix B.

4)Constraints 7 and 8:

The factorsKDminandKDmaxdecide the feasible lower and upper limits for the diameter of rolling element as given in Table B1 of Appendix B. The limits to these parameters have been chosen based on the parametric study by using initial optimization runs as Ref. [33].

5)Constraints 9 and 10:

In order to take care of the running mobility of bearings, the difference between the pitch diameter and bearing average diameter must not exceed a certain specified value(Table B1 of Appendix B). Possible choices for the parameter have been chosen from the parametric study[33].

6)Constraint 11:

The bearing ring thickness at the bottom of the outer raceway must exceed or be equal toεDr, whereεis an unknown constraint constant[22], as given in Table B1 of Appendix B.

7)Constraint 12:

The outer ring thickness must be sufficiently large to allow for the corners to be chamfered (Table B1 of Appendix B).

8)Constraint 13:

It is observed that the stresses developed in the inner ring are always greater than those in the outer ring. That is why there must also be a constraint for the thickness of ring. In CRBs, the design for inner ring must be such that it has greater endurance for stresses than the outer ring since the former has to take more stresses because of its convex curvature (Table B1 of Appendix B).

9)Constraint 14:

Ref.[6] accords to the ISO standards. The value for this maximum contact stress shall vary only if the bearing is a ball bearing in place of roller bearing, or else it shall always be the same for all roller bearings, as per the ISO standards (Table B1 of Appendix B).

10)Constraint 15:

Similarly, the outer ring must also not be stressed beyond the permissible value, as has been shown in Table B1 of Appendix B.

11)Constraint 16:

In roller bearings, to make sure that the rollers are able to roll effortlessly on the inner and outer ring raceways, and not collide among themselves, they must be suitably spaced. In standard bearings, there must at least be an angular spacing of one degree between two rollers. The constraint is explained in Table B1 of Appendix B.

Considering Fig.1, let the angular gap between two consecutive rollers (i.e., angle subtended between the ending of one roller and beginning of immediately next roller) beψ.Let the angle subtended between the endings of two consecutive rollers beφ.Let 2θbe the angle subtended by each of the roller at the center of bearing bore (i.e., angle subtended between the ending and beginning of the same roller). Hence, as can be seen in Fig. 1, there is

ψ=φ-2θ

(23)

(24)

and

(25)

Hence, using Eqs. (23), (24) and (25), there is

(26)

To maintain the angular spacing to always be at least 1°,

ψ≥1°

(27)

Hence, using Eqs. (26) and (27), there is

(28)

Now, Eq. (28) is written as

(29)

12)Constraints 17 and 18:

In the roller bearing design, the roller length also plays an important role, hence its constraining must also be done with utmost care. In the present case,leis considered analogous toDr.leandBput a check on the upper limit ofle(i.e.,le<βB).Moreover, the overall roller length must not exceed a constraint value so that it does not extend from the chamfering of the rings. The same can be seen in Table B1 of Appendix B. The limits to these parameters have been chosen based on the parametric study by using initial optimization runs as in Ref. [33].

13)Constraint 19:

Thickness of the outer ring must be sufficient to ensure the occurrence of the maximum dynamic shear stress at the center of the ring even for the most adverse case. The maximum dynamic shear stress must arise at thrice theZstatic(Table B1 of Appendix B).Zstaticis the depth of occurrence of maximum static shear stress.

14)Constraint 20 & 21:

This constraint is not applicable for deterministic optimization and applies only for robust optimization. In physical programming for the robust design, the designer needs to specify the preference of variation in the objective function due to variation in input variables. Different ranges for percentage variations are prescribed in the objective functionCdfor different CRBs. The reason for prescribing different ranges is elaborated in Section 3. The objective function variations at the design variable tolerance values shall always be within this range and towards the lower limit since the objective function variation is minimized, as can be seen in Eqs. (18) and (19). As is known, each design variable tolerance value input in this study is the maximum possible tolerance for the particular design variable for that CRB, hence the objective function variation corresponding to all those design variable tolerance inputs collectively will also be the most possible for that CRB and objective function. The total objective function variation is directly proportional to each design variable tolerance, as can be seen in Eqs. (13) and (19). Hence, the variation in that objective function for that CRB due to tolerances in the design variables (which shall always be less than the tolerances input by us) shall never be greater than the upper limit of the range of percentage variation prescribed in the constraints. The prescribed ranges of variation in the objective function for each of the CRBs are mentioned in Table 1. The detailed reason behind choosing different values for percentage variation in the objective function for different CRBs is expounded in Section 3. For CRB NU202, as can be seen in Table 1, a 1.2%-1.5% variation is prescribed in the objective functionCd.Hence,

(30)

For CRB NU202, to express the constraint in the form as given in Table B1 of Appendix B, there is

As is known, the above constraint is only to define robustness, hence no modification is needed as done previously for conversion from the deterministic optimization to the robust optimization. To express it in the form as given in Eq. (11), it is

(31)

(32)

As has been studied by Messac and IsmailYahaya[1], the constraints used in the deterministic optimization have been modified in a suitable manner to use it for the robust optimization using Eqs.(8), (10), and (14). All the conversions are in Table B1 of Appendix B (along with their derivations given below the table).

Table 1 Prescribed ranges of variation in objective function for each of the CRBs

1.3 Bounds for Design Variables

Decision of the range for design variables are present in the vector

which can be used as input for any CRB in case of deterministic optimization and the robust optimization. But before this, it must be known how to decideD,d,B,r1min,r2min,r3min, andr4min,Cdand the limiting speed to be input for each bearing. The value for each of these parameters has been taken from Ref.[31] and listed in Table 2. These parameters are used in deciding the range for the basic design variablesDm,Dr,Zandle.The bounds help to define a more concise feasible search space to be used in GAs for faster convergence of the solution. To find the bounds in the case of deterministic optimization for basic design variables for CRB NU202, inputs are boundary dimensions (Table 2), the radial load equals to 10% ofCd(Table 2), the speed of shaft equals to the limiting speed for NU202, i.e., 26000 r/min (Table 2), and the maximum contact stress is as mentioned in Table 3 for the type of bearing (i.e. roller or ball). As is known, the rated life (orL10) for CRB is mathematically given as

(33)

whereL10is in millions of revolutions,Cdis inNandFr(radial load) is inN.Hence, it is preferred to take a lower value of the radial load in order to achieve greater rated life for CRB considering the fact that the aim of this paper is to maximize the fatigue life. Hence, the radial load is chosen as 10% ofCd.However, a different value for the radial load may also be chosen for the present study, while considering the fact that it should be small enough to give a good rated life. It must be noted that the speed of shaft equal to the limiting speed is taken for each of the bearing, as it is unsafe to operate any bearing beyond its limiting speed. The bounds of the basic design variables in case of deterministic optimization have been calculated with the help of the constraints used in the deterministic optimization.

Table 2 Standard CRB parameters used for the design

Table 3 Material properties of the bearing steel

The detailed explanation for how to obtain the bounds for each of the basic design parameters from the constraints (Table B1 of Appendix B) is given in Appendix C. The bounds for these variables shall remain the same even in case of robust optimization when they are used for the deterministic optimization[22], because the variable bounds differ for any CRB for the robust and deterministic optimization only for those variables which have certain tolerances during manufacturing. The bounds for these design variables would also be the same irrespective of the type of CRB.

The bounds for the design variables of the CRB NU202 have been mentioned for the deterministic optimization in Table C1 of Appendix C.

The values of tolerances for the basic design variables and bearing boundary dimensions for each of the CRB have been mentioned in Table D1 of Appendix D. Detailed explanation of how to determine the tolerances for each dimension of the CRBs has been given in Appendix D. While in the case of bearing boundary dimension, the tolerances are always taken as negative. In case of basic design variables, the tolerances may be taken as negative and positive. The input for tolerances in bearing boundary dimensions is taken as negative because they have been conventionally accepted as negative values only. But it must be noted that robust design constraints and the robust design objective function treat negative and positive tolerances (whether of basic design variables or bearing boundary dimensions) in the same manner, as can be seen from Eq.(20) and Eqs.(B.8)-(B.109) in Appendix B.

With the help of robust optimization, a set of design variables shall be achieved, which, even if in the worst case, varies by a value equal to the tolerances (Table D1 in Appendix D). The variation in that objective function for that CRB shall never be greater than the upper limit of the range of percentage variation prescribed in Constraints 20 and 21. To convert the range of bounds for basic design parameters from the deterministic optimization to the robust optimization, Eq. (12) can be used. The pitch diameterDmfor CRB NU202 is considered. For the deterministic optimization and referring to Table C1 in Appendix C, there is

15.6≤Dm≤33.8

(34)

As is given in Table D1 of Appendix D, the variation in the pitch diameter ΔDm(mm) for CRB NU202 is

|ΔDm|=0.010?ΔDm=±0.010

(35)

Hence,with Eqs. (34) and (35), for the robust optimization, the range for pitch diameter is

15.6+ΔDm≤Dm≤33.8-ΔDm

(36)

Although the positive and negative tolerances are treated the same by the robust design constraints and objective functions, different values of optimized design variables and optimized objective functions are obtained on taking a design variable tolerance as negative or positive. This is due to the fact that the sign of design variable tolerance will cause the new robust design bounds of basic design variables to be different, even if the tolerance magnitude remains the same. Hence, although the final optimized objective functions and design variables for the cases of positive and negative tolerances shall have very close values, the probability of them being exactly the same is very low.

However, the final optimized objective functions as well as design variables are not at all affected by the sign of tolerances in bearing boundary dimensions and only affected by their magnitude.

Consider negative tolerances for basic design variables:

Using Eq. (12) and data for CRB NU202 basic design variables from Table C1 and design variable tolerances from Table D1 of Appendix D for ΔDm=-0.010, there is

15.6-0.010≤Dm≤33.8+0.010?

15.59≤Dm≤33.81

(37)

For ΔDr=-0.008, there is

0.2659-0.008≤Dm≤9.1+0.008?

0.2579≤Dm≤9.108

(38)

For Δle=-0.120, there is

0.0106-0.120≤le≤9.35+0.120?

-0.1094≤le≤9.47

(39)

Neglecting the negative value,there is

0≤le≤9.47

(40)

Hence, the bounds for the basic design variables are obtained for the robust optimization of CRB NU202 for negative tolerances in this way. Similarly, the bounds for the robust optimization for the other CRBs are obtained as well, for both negative and positive tolerances. The bounds for the design variables of CRB NU202 have been mentioned for the robust optimization for negative tolerances in Table 4.

Table 4 Bounds of design variables for NU 202 for the robust optimization for negative tolerances

Since the bounds for the design variablesKDmin,KDmax,ε,eandβremain the same irrespective of the type of optimization (whether deterministic or robust) and the type of CRBs, the bounds for the basic design variablesDm,Drandlefor the robust optimization of each of the CRBs for the prescribed maximum positive and negative tolerances (Table D1 of Appendix D) have been mentioned in Table C2 of Appendix C. It must be noted that the bounds for the design variableZremain the same irrespective of the type of optimization (whether deterministic or robust), but they do vary with the type of CRB being considered.

1.4 Constraint Handling

Constrained optimization problems involve inequality and/or equality constraints. Solving such problems using GA requires a constraint handling method. There are several ways to handle constraints, but owing to the simplicity and ease of operation, both the penalty function[34]and penalty parameter-less[35]approaches are frequently applied, although the method best suited for optimization algorithms working with population is the penalty parameter-less method. The aforesaid method incurs a summation of penalty terms 〈gj(X)〉 and 〈hk(X)〉 corresponding to the constraint violation functionCV(X).It calculates the overall normalized constraint violation of the solutionXas follows:

(41)

2 Implementation of NSGA-II to Rolling Element Bearings

A flow chart of NSGA-II[8,29]is shown in Fig.2, which gives the procedure for solving the present optimization problem. It must be noted that although the different multi-objective optimization GAs work on slightly different paths of approach when dealing with a multi-objective optimization problem, they behave in exactly the same manner when dealing with a single objective optimization problem as the objective function. Hence, whenever a multi-objective optimization GA is used for solving single-objective optimization problems, the usual entire procedure, i.e., initialization, evaluation, selection, reproduction, re-evaluation, and so on, stays the same, as it would be for a simple real-coded GA. But what makes NSGA-II the best multi-objective optimization algorithm is its low computational complexity and diversity preservation. It also has the advantage of elitism, which is also the advantage of the most other multi-objective optimization procedure. So, when NSGA-II is used for single objective optimization, the benefits it brings along are low computational complexity and diversity preservation, which are not present in the other potential single-objective optimization methods, namely, evolution strategies, evolutionary programming, genetic programming, etc. This is the reason why NSGA-II is used but not any other available multi-objective optimization methods in place of the other available potential single-objective optimization methods.

Fig.2 Flowchart of NSGA-II algorithm

In the present problem of CRB, nine design variables have been considered and encoded in real-coded chromosomes, thereafter appropriate modifications are made to the NSGA-II code to take the number of rollersZas integer value. Table 3 lists the material properties like Young’s modulus, Poisson’s ratio and acceptable contact stressσCsafe.For the bearing, the alloy steel 52100 grade bearing steel is chosen as the material due to its features of high wear resistance, rolling fatigue strength, cost-effectiveness, and long working life. These values are used as inputs to decide the range for the basic design variablesDm,Dr,Z, andlefor the deterministic optimization (Appendix C), and also as inputs for the optimization constraints (Table B1 of Appendix B).

The bounds for the design variables of the CRB NU202 are mentioned for the deterministic optimization in Table D1. With Eq. (12), the bounds of the basic design variablesDm,Drandlefor deterministic optimization (Table D1) along with data for design variable tolerances from Table D1 are taken as inputs to decide the range for the basic design variables for the robust optimization in the similar manner as that in Eqs. (37) - (40). The values of tolerances for the basic design variablesDm,Dr, andle, and bearing boundary dimensions for each of the CRB mentioned in Table D1 are taken as inputs for robust optimization constraints, as listed in Table B1 of Appendix B.

The bounds for the design variables of the CRB NU202 are mentioned for the robust optimization for negative tolerances in Table 4. The bounds for the basic design variablesDm,Dr, andlefor the robust optimization for each of the CRBs for the prescribed maximum positive and negative tolerances (Table D1 of Appendix D) are mentioned in Table C2 of Appendix C.

3 Results and Observations

The results obtained after the robust optimization of NU202 for negative tolerances (Table D1 of Appendix D) and prescribed variation of 1.2%-1.5% in the objective function (Table 1) for the population size 1000 and different values of number of generations are given in Table 5. In the case of the CRB NU202, for a population size of 1000, the number of generations used for obtaining the final solution is 125. The reason for using this value for the number of generations corresponding to population size 1000 is the convergence of values of the objective function towards a higher magnitude and different design variables towards a particular value as the chosen value of number of generation approaches (Table 5). Convergence of the different design variables and the objective function towards increasing value of number of generations can be seen in Figs. 3, 4 and 5. Even higher values of number of generations can be used for the given population size, but it is a futile exercise, as the results obtained are almost the same as those obtained now, and it is also computationally exhaustive (Table 5 and Fig. 6) because the execution time keeps on increasing by increasing number of generations. Hence, to minimize the execution time as much as possible, a minimum possible number of generations are kept and these values are strictly stuck to. It must be noted that for a different value of population size, the value for the number of generations at which convergence occurs might be different. After finalizing the population size and the number of generations for each run, GA parameters (mutation probability, crossover probability, mutation index, crossover index and random seed) are chosen. The lower value of GA parameters gives the maximum spread of solution points and hence diverse results. Therefore, lower values are chosen for GA parameters for the present study[26]. As mentioned earlier, a random seed has also been used as the GA parameter and used to generate initial population for each run. GA parameter variations and their favoured values have been given in Table 6. Similarly, NSGA-II can be implemented for other standard CRBs for the negative as well as positive tolerances for any prescribed range of variation in the objective function with the data from Tables 2, 3, 6, 7 and 8. A certain value is used for population size and the number of generations for the robust optimization ofCd.In this study, a population size of 1000 is kept and the number of generations is varied up to the point of convergence in design variables and objective functions. The final optimized dimensions for the robust optimization ofCdfor each of the CRB are mentioned in Table 2. The negative as well as positive tolerances after prescribing suitable ranges of variation in the objective function have been provided in Table 7.

Fig.3 Convergence of Dm(mm), Dr (mm), Z and le (mm) with number of generations for robust optimization of CRB NU202 for negative tolerances

Fig.4 Convergence of KDmin, KDmax and ε with number of generations for robust optimization of CRB NU202 for negative tolerances

Fig.5 Convergence of e, β and Cd(N) with number of generations for robust optimization of CRB NU202 for negative tolerances

Table 5 Robust optimization results for Cd of CRB NU202 for different values of number of generations

Table 6 GA parameter variations and favoured values

Fig.6 Increase in execution time with number of generations for robust optimization of CRB NU202 for negative tolerances

Table 7 Robust optimization results for dynamic capacity of different CRBs for negative and positive tolerances

Table 8 Robust optimization results for Cd of CRB NU202 for negative tolerances and different values of prescribed variation in objective function

It can be seen from Table 8 that when NSGA-II is implemented for the robust optimization ofCdfor NU202 for negative tolerances, and different ranges of variation are prescribed in the objective function and using data from Tables 1, 2, 3 and 4. The final optimized values of design variables and objective functions have decreased from the previously prescribed range of variation 1.2%-1.5%. This is because while the tolerances input for CRB NU202 remains the same, an increase in prescribed range of variation demands a decrease in one or all of the basic design parametersDm,Dr, andle, as can be seen in Eqs. (20)-(22). Due to this decrease in these basic design parameters,Cdhas decreased. Although a highly flexible design of the bearing is preferred (i.e., acceptability of high percentages of variations in the objective function), the importance of a high value of the dynamic capacity cannot be ignored. Hence, a balance needs to be maintained between the flexibility in design and load carrying capacity. Hence, by observing the results from Table 8, the final range of variation was prescribed to be 1.2%-1.5% for the robust optimization ofCdfor NU202 for negative tolerances, considering the fact that the final optimized objective function after performing the robust optimization for any objective for any CRB using positive/negative tolerance must be comparable to the value obtained on the deterministic optimization for that objective function for that CRB[22]. For example, Ref. [22] proposes that for NU202, for single objective optimization ofCd, (Cd)optimum=22160 N. From Table 8, an optimized value 22081 N is obtained for the robust optimization ofCdfor NU202 for negative tolerances, which lies near the above value and is only for prescribed range of variation 1.2%-1.5%. Hence, 1.2%-1.5% is chosen as the prescribed range of variation.

Similarly, different ranges of variation in the objective function can be prescribed for the robust optimization ofCdfor each of the CRB mentioned in Table 2, for both negative and positive tolerances. The prescribed ranges of variation in the objective function for each of the CRBs have been mentioned in Table 1. Now, after deciding a prescribed range of variation of 1.2%-1.5% for the robust optimization ofCdfor NU202 for negative tolerances, NSGA-II can be implemented for the robust optimization ofCdfor NU202 for positive tolerances, prescribing the same range of variation in the objective function as the case of negative tolerance and using data from Tables 1, 2, 3 and 4. Now, the bounds of basic design variablesDm,Dr, andlehave changed from the case of negative tolerance for NU202 (Table C2 of Appendix C). Hence, although the final optimized values of design variables and objective functions are very close to the case of negative tolerance for NU202, they are surely different as can be seen in Table 9. But the idea of still using the same prescribed range of variation as negative tolerance in case of positive tolerance does turn out to be useful, as the optimized value ofCd=21990 N (Table 9) is still comparable to the value (Cd)optimum=22160 N[22]. Hence, the same range of variation is prescribed in the objective function for the robust optimization ofCdfor each of the CRB for negative and positive tolerances.

To check whether the optimized dimensions are geometrically feasible or not, NU202 bearing is drawn for dimensions obtained using the robust optimization ofCdof NU202 for negative tolerances, as given in Table 9, which uses a solid modelling software and is shown in Fig. 7, where the cross-sectional view of the optimized bearing NU-202 is shown. From the figures, it can be inferred that no interference is there between the dimensions.

The objective function variations of the design variable tolerance values input in this study shall always be within the range prescribed by us in Constraints 20 and 21.

The problem that the objective functionCdvaries with each of the basic design variableDm,Drandleundergoing manufacturing tolerances is investigated. Hence, using results from Table 9 for the negative tolerances of CRB NU202 and Eqs. (20) - (22), there is

(42)

Table 9 Robust optimization results for Cd of CRB NU202 for negative and positive tolerances for a prescribed variation of 1.2%-1.5% in objective function

Fig.7 Radial dimensions of optimized bearing for robust optimization of NU202 with Cd as objective function

Hence, the change rate ofCdwith each of the design variablesDm,Drandleis positive. It means that for a positive tolerance in any of the dimensionsDm,Drandleduring the manufacturing of bearings, the objective functionCdwill increase, while a negative tolerance in the same will cause a decrease in theCd.The objective functionCdalso changes at the maximum rate forDrand it is least affected by the changes inDm.Although a positive tolerance for the design variablesDm,Dr, andledoes not seem to be a problem, it causes an increase in the objective functionCd, which is good for the bearing. But in that case, if the final optimized dimensions are obtained using the robust optimization and only the deterministic optimization is performed, the tolerances may cause the final manufactured dimensions to be violating one or more than one of the Constraints 1-19.

The total percentage variation in the objective functionCdare due to the tolerances in each of the design variablesDm,Dr, andle.Hence, using the results from Eq. (42) and with the data from Table D1 for the negative tolerances of CRB NU202, there is

(43)

Hence,with the results obtained in Eq. (43), there is

(44)

The negative sign in Eq. (44) has been mentioned to show that the net variation inCdof the CRB NU202 for negative tolerances in the design variablesDm,Dr, andlewill be negative. As is observed above, with Eq. (43), a negative tolerance in the design variablesDm,Dr, andlewill cause a decrease in theCd.

As is mentioned, the objective function variation magnitudes at the design variable tolerance values input for any CRB for any objective function shall always be within the range prescribed in this study for the CRBs in Table 1 and towards the lower limit that has been observed in Eq. (44). From Table 1, the prescribed range of variation for CRB NU202 for robust optimization ofCdis 1.2%-1.5% and a net variation of 1.2% is obtained in the objective functionCd, as can be seen in Eq. (44).

But it must be insured that the optimum results obtained for the robust optimization satisfies each of the constraints. To this end, the values for the left side of each of the robust optimization constraints are calculated, which is given in Table B1 of Appendix B, corresponding to the results obtained for the robust optimization of NU202 for negative tolerances and prescribed variation in the objective function as mentioned in Table 9. The obtained values for left side corresponding to each constraint have been mentioned in Appendix E, where all the values are greater than zero, hence it is certain that the optimum solution obtained for robust optimization of NU202 for negative tolerances and prescribed variation in the objective function does satisfy all the constraints. Hence, it can be concluded that the results in Table 7 obtained on robust optimization of each of the CRBs for negative and positive tolerances and prescribed variation in the objective function do satisfy the 21 robust optimization constraints, as mentioned in Table B1.

Consider Constraint 1 and the value of its left side in Appendix E. The value of left side for Constraint 1 is 10.4498, which indicates that if only Constraint 1 needs to be followed to find out the value ofDm, it (obtained through optimization) can be still decreased by 10.4498 mm. Similarly, the value of left side for Constraint 2 is 7.7112, which means that if only Constraint 2 needs to be followed to find out the value ofDm, it (obtained through optimization) can still be increased by 7.7112 mm. But in order to determine the value of optimumDm, both Constraints 1 and 2 have to be followed in addition to several other constraints out of the total 21 available constraints. Hence, the optimum value is obtained forDm.The condition that each of the constraint containsDmas a design variable has been satisfied. It must be noted that each of the nine design variables considered in this study may be the only design variable in a constraint, or it may be present as one of the several design variables in a constraint.

Similarly, consider Constraint 7, the value of its left side in Appendix E is 1.0961, which indicates that if only Constraint 7 has to be followed to find out the optimum values of the design variablesDrandKDmin, both values (obtained through optimization) can still be changed simultaneously or individually, such that the value of the expression

{2Dr-KDmin(D-d)-KDmin(|ΔD|+|Δd|)}

does not decrease by more than 1.0961 mm. But in order to determine the value of optimumDrandKDminin addition to Constraint 7, several other constraints out of total 21 constraints have to be followed. Hence, the optimum value for each of the design variables,DrandKDminare obtained. The condition that each of the constraint containsDrandKDminas the design variable(s) in combined form or separately has been satisfied.

For any of the constraint (whether consisting of just one design variable or a set of design variables), the nearer the value of LHS to zero is, the more critical the constraint would be. This criticality of a constraint can be seen in reference to all the constraints or in reference to only those constraints, sharing at least one common design variable. Consider the criticality of a constraint in reference to all the constraints. For the set of design variables present in this most critical constraint, the convergence to a stable value shall occur more quickly than convergence to a stable value for the set of design variables present in less critical constraints. The explanation for it has been given in the next paragraph. Consider the criticality of a constraint in reference to the constraints sharing at least one common design variable. Then out of all the constraints consisting of a design variable, the final optimized value for that design variable shall depend on the constraint that is most critical for that design variable.

In Appendix E for Constraints 11 and 16, the value of left-hand side corresponding to the optimum results for the robust optimization of NU202 for negative tolerances, as given in Table 9, lies very close to 0, which indicates that these are the most critical constraints (in reference to all the constraints). Being a critical constraint means that for each run, these constraints will limit the range of design variables (out of the range that has been input) to optimize the objective function. Hence, the algorithm will keep only those populations as favourable parent populations in further generations, which do satisfy these constraints (i.e., the most critical ones). Hence, for the design variables that are present in these critical constraints, parent populations are carried forward for many generations. Hence for those design variables, the convergence to a stable value can be said to have occurred very quickly, i.e., in the initial generations itself. Hence, in the present case of robust optimization of the dynamic capacity of NU202 for negative tolerances (Table 9) and also for the robust optimization of each of the CRBs for positive/negative tolerance (Table 7), design variablesDm,Dr,εandZpresent in Constraints 11 and 16 shall converge to a stable value more quickly than the other design variables.

As can be observed in the final results for the optimization ofCdfor CRB NU202 (Tables 7 and 9), while some of the design variables likele,KDmin,KDmax,ε,eandβare converging to the boundary values of their prescribed range, the others likeDm,DrandZconverge to a value lying somewhere between the bounds of prescribed range. First of all, the point of convergence of basic design variables,Dm,Dr,Zandlerepresenting the bearing internal geometry are discussed. As can be seen in Section 1.3, the upper and lower bounds for these basic design variables for robust optimization have been derived from the bounds for these design variables for the deterministic optimization, which themselves are based on a few of the constraints from Table B1 of Appendix B, as can be seen in Appendix C. As has been mentioned previously, if the criticality of a constraint is considered in reference to the constraints having at least one common design variable, out of all the constraints consisting of a design variable, the final optimized value for that design variable shall depend on the constraint that is most critical for that design variable. Hence, for any of the basic design variables, convergence to a boundary value of the prescribed range shall occur, only if the critical constraint for that design variable is the same as the constraint that had been used to find its upper and lower bound for the deterministic optimization in Appendix C. The point of convergence of the remaining five design variables are discussed, namely,KDmin,KDmax,ε,e, andβthat give the bearing a feasible design space in the form of constraints, as given in Table B1 of Appendix B.

To understand the convergence ofKDminandKDmax, consider Constraints 7 and 8, as can be seen in Table B1 of Appendix B.

(45)

As can be seen in Eq. (43),the change rate ofCdwith respect to the design variableDris positive, which means that in order to maximizeCd(which is the purpose of this study),Drmust also be high. As can be seen in Eq. (45), in order to make the search space forDrcorresponding to Constraints 7 and 8, it only consists of higher values forDr.The upper and lower bounds for this constraint must also be high, for which high values ofKDminandKDmaxare required. Hence, the optimization algorithm tries to obtain values for these design variablesKDminandKDmaxas high as possible, while always takes into consideration the fact that the values to which the design variablesDr,KDminandKDmaxare constrained through this constraint do not violate any of the remaining 19 constraints. The optimization algorithm succeeds in increasing the values for these design variables,KDminandKDmax, up to the maximum possible value (i.e., the upper bound of prescribed range) without violating any of the remaining 19 constraints, i.e., the values to whichKDminandKDmaxconverge to are the upper bounds for their prescribed range.

To understand the convergence ofe, consider Constraints 9 and 10, as can be seen in Table B1 of Appendix B.

(0.5-e)(D+d)≤Dm≤(0.5+e)(D+d)

(46)

With the help of Eq.(43),Dmmust also be high in order to maximize dynamic capacity. As can be seen in Eq. (46), a high value ofewill give a higher value for the upper bound for search space ofDmcorresponding to Constraints 9 and 10. As can be seen in Eq. (46), a higher value ofeshall cause the lower bound for the search space ofDmcorresponding to Constraints 9 and 10 to be lower, which is not a problem, as the optimization algorithm can still look for higher values forDm.Hence, a highest possible value foreis searched for, such that the values to whichDmandeare constrained through this constraint do not violate any of the remaining 19 constraints. Hence, the value to which the design variableeconverges to is the upper bound for its prescribed range.

To understand the convergence ofε, consider Constraint 11, as can be seen in Table B1 of Appendix B.

0.5(D-Dm-Dr)-εDr≥0?

Dm+(1+2ε)Dr≤D

(47)

As can be seen in Eq.(47), a smaller value ofεwill increase the upper bound for search space of bothDmandDr, corresponding to Constraint 11, to an even higher value. Hence, the optimization algorithm converges the design variableεto the least possible value (i.e., lower bound for its prescribed range), meanwhile ensures that the values to which the design variablesDm,Drandeare constrained through this constraint do not violate any of the remaining 20 constraints.

To understand the convergence ofβ, consider Constraint 17, as can be seen in Table B1 of Appendix B.

βB-le≥0 orle≤βB

(48)

As can be seen in Eq. (43),Cdcan be maximized by increasingle.As can be seen in Eq. (48), a higher value ofβis well suited to make the upper bound of the search space forlehigher. Hence, the value of the design variableβis increased as much as possible, while always takes care of the fact that the values to which the design variablesleandβare constrained through this constraint do not violate any of the remaining 20 constraints. Hence, the design variableβconverges to the upper bound of its prescribed range.

4 Conclusions

The present study is primarily focused on having an optimum design of cylindrical bearing that not only gives us the maximum load carrying capacity but also remains least affected by the possible practical error that can occur in the form of tolerances during its manufacturing. The study presented some mathematical formulations that make it easier to achieve the objective of having a robust design. The formulations have shown that in order to achieve the robustness, adjustments need to be made in all three components of any optimization process, i.e., the objective function, constraints and variable bounds. Hence, based on the results of the present study, the following conclusions can be drawn:

1)It is possible to optimize a bearing design in such a manner that its variations due to tolerances in design variables (whose maximum value has been prescribed by us in both positive and negative direction) are never greater than the upper limit of the range of percentage variation prescribed for it during manufacturing. This is possible to be manufactured, as solid modelling software can be used to draw such a bearing with geometrical parameters as obtained by the robust optimization.

2)The change rate for the objective functionCdcorresponding to each of the design variablesDm,Drandleundergoing variations during manufacturing is positive.

3)Out of all the parameters, which have been given tolerance and affect the objective function, the roller diameter has the maximum effect on dynamic capacity and the mean diameter has the least.

In this work, the robust design of CRBs has been carried out with a single objective function, namely, the dynamic capacityCd.However, bearings may have different application requirements, so a multi-objective optimization of robust design for CRBs can be carried out by considering simultaneously and in pairs of objective functions asCd, maximum bearing temperatureTmax, and the elasto-hydrodynamic minimum film thicknesshminas objective functions. Apart from this, the vibration performance and fatigue life of a rolling element bearing can be significantly affected by variable compliances of the inherent bearings, the rotor, and housing compliances, and these can be considered as a future work.

Appendices

Note: The electronic version of Appendix A,B,C,D,E is available at http://hit.alljournals.cn/jhit_cn/ch/reader/view_abstract.aspx?file_no=202211070000001&edit_id=20221107103135001& flag=2.