Optimization Design of RV Reducer Crankshaft Bearing

: The crankshaft bearing is the key component of a rotate vector (RV) reducer. However, owing to the harsh working load and restricted available space, the bearing often su ﬀ ers from fatigue failure. Therefore, this study proposes a novel optimization method for RV reducer crankshaft bearings. A nonlinear constraint optimization model for the design of the bearing considering the crowned roller proﬁle is formulated and is solved by using a crow search algorithm. The goal of the optimization is to maximize the fatigue life of the bearing. The design variables corresponding to the bearing geometry and crowned roller proﬁle are considered. The load working conditions of the bearing and structure of the RV reducer are analyzed. Various constraints, including geometry, lubrication, strength of the bearing, and structure of the RV reducer, are established. Through the optimization design, the optimum crowned roller proﬁle suitable for the working load of the bearing is obtained, and the stress concentration between the roller and raceway is eliminated. Taking the crankshaft bearing of RV-20E and RV-110E type reducers as examples, the bearings were optimized by the proposed method. After optimization, the bearing life of the RV-20E type reducer is increased by 196%, and the bearing life of the RV-110E type reducer is increased by 168%.


Introduction
The rotate vector (RV) reducer is widely used in industrial robots because of its compact structure, small size, high torsional rigidity, large transmission ratio, and great load capacity [1,2]. Service life is one of the most important performance parameters of a reducer. However, the harsh working environment makes the reducer prone to fail easily [3]. Thus, it is important to investigate the fatigue life of the RV reducer.
For an RV reducer, the service life mainly depends on the rated life of crankshaft bearings [4,5]. Both the RV reducer sample manuals of Nabtesco and standard GB/T37718 take the rated life of the crankshaft bearing as the service life of the reducer. However, the load working condition of the crankshaft bearings are very complex, and the radial size is very small. Therefore, the crankshaft bearings often suffer fatigue failure, which limits the service life of the RV reducer [4]. Therefore, crankshaft bearings use non-standard cylindrical roller bearings [4,5]. For nonstandard bearing, it takes a lot of time to design the bearing with the traditional trial and error method [6]. In addition, the design of the bearing must satisfy various constraints, e.g., geometrical, kinematical, and mechanical constraints. Therefore, it is necessary to use optimization techniques for crankshaft bearings to improve the service life of the reducer.
In recent years, there has been abundant research focused on RV reducers. Li [7] proposed an analytical model for the cycloidal gear considering tooth profile modification. Han [8] analyzed the sensitivity of transmission accuracy of RV reducers. Xu [9] researched the measurement of the Crankshaft bearing is mainly composed of rollers, inner raceways, and outer raceways. Because of the restriction of available space, a cylindrical roller bearing without inner and outer rings is used for the crankshaft bearing. As shown in Figure 2, the outer raceway is integrated into the cycloidal gear. Crankshaft bearing is mainly composed of rollers, inner raceways, and outer raceways. Because of the restriction of available space, a cylindrical roller bearing without inner and outer rings is used for the crankshaft bearing. As shown in Figure 2, the outer raceway is integrated into the cycloidal gear.
Appl. Sci. 2020, 10 The size of the outer raceway is limited by the root circle and the center hole of the cycloidal gear. The inner raceway is integrated into the off-center section of the crankshaft. The radius of the crankshaft bearing distribution circle determines the installation position of the crankshaft bearing in the cycloidal gear. The width of the off-center section is the same as the width of the cycloidal gear, which limits the length of the roller.
Although the components of the crankshaft bearings appear simple, the macro geometry has significant effects on mechanical behavior of bearing. The macro-geometric parameters mainly include the roller length , roller diameter , bearing pitch diameter , and radial clearance P . As shown in Figure 2, the boundary dimensions include the outer diameter ( ), minimum diameter of crankshaft ( ), bearing width (B), diameter of the center hole ( ), and the root circle of the cycloidal gear ( ).

Force Analysis of Crankshaft Bearings
A Cartesian coordinate system, , is employed to analyze the load of the crankshaft bearings. The origin is located at the center of the distribution circle of the pin gears. The direction of the axis coincides with the eccentric direction of the crankshaft. The cycloidal gear bears the forces from the pin gear and crankshaft bearings. The component forces of pin gears on cycloid gear along theaxis and -axis are ∑ and ∑ , respectively. The forces acting on the cycloidal gear by the crankshaft bearing are decomposed into, , , and . The torque caused by the force is balanced with the torque caused by the force ∑ . The force is balanced with force ∑ and is balanced with ∑ . This can be obtained from the force balance on the cycloidal gear: where ∑ = ∑ and = + ln , is the short width coefficient and = / , is the number of pin teeth, is the eccentric distance of the crankshaft, is the number of cycloidal gear teeth, and and is the radius of distribution circle of the pin gear. It is assumed that the output torque of the reducer is . For the torque equilibrium of the cycloidal gear, The size of the outer raceway is limited by the root circle and the center hole of the cycloidal gear. The inner raceway is integrated into the off-center section of the crankshaft. The radius of the crankshaft bearing distribution circle determines the installation position of the crankshaft bearing in the cycloidal gear. The width of the off-center section is the same as the width of the cycloidal gear, which limits the length of the roller.
Although the components of the crankshaft bearings appear simple, the macro geometry has significant effects on mechanical behavior of bearing. The macro-geometric parameters mainly include the roller length L we , roller diameter D we , bearing pitch diameter D m , and radial clearance P d . As shown in Figure 2, the boundary dimensions include the outer diameter (D U ), minimum diameter of crankshaft (D L ), bearing width (B), diameter of the center hole (D h ), and the root circle of the cycloidal gear (D c ).

Force Analysis of Crankshaft Bearings
A Cartesian coordinate system, xoy, is employed to analyze the load of the crankshaft bearings. The origin is located at the center of the distribution circle of the pin gears. The direction of the y axis coincides with the eccentric direction of the crankshaft. The cycloidal gear bears the forces from the pin gear and crankshaft bearings. The component forces of pin gears on cycloid gear along the x-axis and y-axis are P ix and P iy , respectively. The forces acting on the cycloidal gear by the crankshaft bearing are decomposed into, F j1 , F j2 , and F j3 . The torque caused by the force F j1 is balanced with the torque caused by the force p ix . The force F j2 is balanced with force P ix and F j3 is balanced with P iy . This can be obtained from the force balance on the cycloidal gear: where P iy = K y P ix and K y = 2 π 1 k + k 2 −1 2k 2 ln 1+k 1−k , k is the short width coefficient and k = eZ b /R z , z b is the number of pin teeth, e is the eccentric distance of the crankshaft, z g is the number of cycloidal gear teeth, and R z and is the radius of distribution circle of the pin gear.
It is assumed that the output torque of the reducer is T v . For the torque equilibrium of the cycloidal gear, Appl. Sci. 2020, 10, 6520 4 of 20 In Figure 3b, o 1 is the center of the cycloidal gear. o c is the center of the crankshaft bearing. The coordinate system x 1 o c y 1 is rigidly connected with the crankshaft bearing. The angle between the y 1 and y axes is θ. It is assumed that the force between the crankshaft bearing and the cycloidal gear is F. F can be decomposed into the normal force F r and tangential force F t . F r is in the direction of the axis x 1 . F t is in the direction of the axis y 1 . The forces acting on the cycloidal gear by the crankshaft bearing are decomposed into F j1 , F j2 , and F j3 . According to the balance between the acting and reaction forces, Appl. Sci. 2020, 10, x FOR PEER REVIEW 5 of 23 In Figure 3b, is the center of the cycloidal gear. is the center of the crankshaft bearing. The coordinate system is rigidly connected with the crankshaft bearing. The angle between the and axes is . It is assumed that the force between the crankshaft bearing and the cycloidal gear is .
can be decomposed into the normal force and tangential force . is in the direction of the axis .
is in the direction of the axis . The forces acting on the cycloidal gear by the crankshaft bearing are decomposed into , , and . According to the balance between the acting and reaction forces, cos sin sin cos From the above equation, the load of the crankshaft bearing changes periodically as the crank shaft rotates. The load change period was 2 . According to Harris et al [19], the equivalent load of the crankshaft bearing can be calculated by The equivalent load of the bearing obtained through integration is

Fatigue Life of Crankshaft Bearings
In cases where the rolling elements are only lightly profiled or subjected to excessive loads, edge stresses can arise, which need to be considered in the rating life calculation. Therefore, both the roller profile and macro-geometric parameters are considered in this study. First, the load distribution is Therefore, the resultant force acting on the crank bearing is From the above equation, the load of the crankshaft bearing changes periodically as the crank shaft rotates. The load change period was 2π. According to Harris et al. [19], the equivalent load of the crankshaft bearing can be calculated by The equivalent load of the bearing obtained through integration is

Fatigue Life of Crankshaft Bearings
In cases where the rolling elements are only lightly profiled or subjected to excessive loads, edge stresses can arise, which need to be considered in the rating life calculation. Therefore, both the roller profile and macro-geometric parameters are considered in this study. First, the load distribution is calculated using the quasi-static method. Then, the distribution of the contact stress over the length of the rollers is calculated using the non-Hertzian contact method. Finally, the fatigue life of the bearing is evaluated by considering the contact stress of the rollers. As shown in Figure 4, the radial load F m from the cycloidal gear acts on the crankshaft bearings. It is assumed that the relative radial displacement between the crankshaft and the cycloidal gear is δ r . It can be seen from the deformation coordination relationship that the radial deflection at any roller angular position ψ j is where ψ j is the angular position of the jth roller and P d is the radial clearance.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 6 of 23 calculated using the quasi-static method. Then, the distribution of the contact stress over the length of the rollers is calculated using the non-Hertzian contact method. Finally, the fatigue life of the bearing is evaluated by considering the contact stress of the rollers.

Load Distribution of Crankshaft Bearings
As shown in Figure 4, the radial load from the cycloidal gear acts on the crankshaft bearings. It is assumed that the relative radial displacement between the crankshaft and the cycloidal gear is . It can be seen from the deformation coordination relationship that the radial deflection at any roller angular position is where is the angular position of the th roller and is the radial clearance.

Load distribution
Outer raceway Inner raceway Cycloidal gear The total normal contact deformation between two raceways separated by a rolling element under a load is the sum of the deformations between the rolling element and each raceway. Hence, where is the normal contact deformation between inner raceway and the th roller, and is the contact deformation between the outer raceway and the th roller.
The roller at the angular position is subjected to forces from the raceways / . The normal roller loads are related to normal contact deformations as follows [19]: According to Harris [19], for the steel roller-steel raceway contact, The roller is subjected to the forces from the raceways and the centrifugal force . Thus, considering the force equilibrium, The centrifugal force for a steel roller orbiting at speed about a bearing axis is given by The total normal contact deformation between two raceways separated by a rolling element under a load is the sum of the deformations between the rolling element and each raceway. Hence, where δ ij is the normal contact deformation between inner raceway and the jth roller, and δ oj is the contact deformation between the outer raceway and the jth roller. The roller at the angular position ψ j is subjected to forces from the raceways Q oj /Q ij . The normal roller loads are related to normal contact deformations as follows [19]: According to Harris [19], for the steel roller-steel raceway contact, The roller is subjected to the forces from the raceways and the centrifugal force F c . Thus, considering the force equilibrium, The centrifugal force for a steel roller orbiting at speed n m about a bearing axis is given by Appl. Sci. 2020, 10, 6520 6 of 20 Substituting Equations (7)-(10) into Equation (11), Equation (11) can be rewritten as follows: For static equilibrium, the radial load must be equal to the sum of the vertical components of the roller loads: Equations (13) and (14) are a system of simultaneous nonlinear equations. δr and δ ij are unknown and need to be solved. The equations are solved by Newton-Raphson method. Thereafter, roller loads are obtained using Equation (9).

Roller-Raceway Contact Stress
Hertzian contact is often used to solve the contact stress between two cylinders. However, it cannot reveal the stress concentration and the effect of the roller profile. Non-Hertzian contact method can reveal the effect of the roller profile. The distribution of the contact stress over the length of the rollers can be obtained by using non-Hertzian theory. Therefore, it is adopted in this study.
The contact forces between the raceway and roller were calculated using the quasi-static model. When the roller is subjected to a normal load Q, it is assumed that the contact area between the raceway and the roller is A. The contact pressure at coordinate (x, y) is P(x, y). The integral of pressure P(x, y) is equal to the contact force Q.
Contact stress occurs when two surfaces come into contact with each other. Assuming that the contact deformation between the two surfaces at coordinate (x, y) is δ(x, y), then According to the half-space theory, when two bodies come into contact, the elastic deformation δ(x, y) can be expressed by the following equation: where E is the equivalent Young's modulus.
where ν is the Poisson's ratio and E is the elastic modulus. When the roller is subjected to the normal load Q, the contact surfaces are moved toward each other by a displacement δ. For a roller with a crowned profile, z(x, y) is the gap between the roller and raceway at the coordinate (x, y) generated by the crown drop. f (x, y) is the initial distance generated by the geometry of the roller and raceway. Therefore, Equation (17) can be written as follows: Appl. Sci. 2020, 10, 6520 7 of 20 In these equations, the relative displacement δ, pressure distribution P(x, y), and contact area A are unknown. To simplify this problem, the contact area is subdivided into multiple discrete cells. It is assumed that the contact pressure distribution on each cell is constant. Therefore, the integral is replaced by a summation, represented as where f ij is the influence coefficient representing the deflection of cell i because of a uniform pressure over the cell j, which can be obtained according to reference [20]. D i is the contact deformation of cell i. Finally, the equilibrium condition of the applied contact force and the integral of the pressure over the contact yields: Equations (20) and (21) are solved simultaneously for the contact area and contact stress.

Fatigue Life
The basic dynamic load rating of the raceways Q ci /Q co are obtained by 143/108 9/2 2/9 143/108 −9/2 2/9 (22) According to ISO/TR 1281-1, the value of λη is 0.83. γ is the parameter related to the structure of the bearing, γ = D we /d m . The basic dynamic load capacity C r of the bearings is where b m accommodates recent improvements in bearing quality and standard ISO/TR 1281-1 gives values of the factor b m = 1.1. f c is a factor related to the the geometry and the material of the bearings, The contact stress along the roller over the length of the roller is considered. To this end, the rollers were discretized into a number (n s ) of slices. The basic load capacity of a slice of the inner and outer rings, The dynamic equivalent load on the kth slice is: Appl. Sci. 2020, 10, 6520 8 of 20 where P ok j is the contact stress at the kth slice of jth roller of outer ring and P ik j is the contact stress at the kth slice of jth roller of inner ring. The pressures were obtained in the previous section. Finally, the reference bearing fatigue life is obtained as

Mathematical Model of Optimization
Crankshaft bearings are the key components restricting the service life of RV reducers. Therefore, crankshaft bearings are optimized in this study. The objective functions, design variables, and constraints in the optimization design are as follows:

Definition of Design Variables
The mechanical behavior of the crankshaft bearing is significantly affected by the macro-geometric structure of the bearing. The crowned profile gives an axially uniform pressure distribution. These factors significantly affect the fatigue life of the bearing. Therefore, the macro geometry and roller profile are considered in the optimization design. The design parameters of the macro geometry mainly include the roller diameter D we , number of rollers Z, bearing pitch diameter D m , and radial clearance P d . In addition, the radius of the crankshaft distribution circle (r o ) determines the installation position of the crankshaft in the cycloidal gear and the force of the crankshaft bearing.
When the rolling elements are only lightly profiled or subjected to excessive loads, edge stresses can arise. Cylindrical rollers (or raceways) are usually crowned to solve this problem [19]. Fujiwara [21,22] proposed a modified logarithmic shape that allowed for a flat region on the roller, as shown in Figure 5. The crown drop in the axial direction at position x is where Q max is the maximum load acting on the rollers, K 1 is the load coefficient, K 2 is the coefficient of the crown length, and Z m is the crown drop at the edge of the roller end.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 9 of 23 where is the contact stress at the th slice of th roller of outer ring and is the contact stress at the th slice of th roller of inner ring. The pressures were obtained in the previous section.
Finally, the reference bearing fatigue life is obtained as

Mathematical Model of Optimization
Crankshaft bearings are the key components restricting the service life of RV reducers. Therefore, crankshaft bearings are optimized in this study. The objective functions, design variables, and constraints in the optimization design are as follows:

Definition of Design Variables
The mechanical behavior of the crankshaft bearing is significantly affected by the macrogeometric structure of the bearing. The crowned profile gives an axially uniform pressure distribution. These factors significantly affect the fatigue life of the bearing. Therefore, the macro geometry and roller profile are considered in the optimization design. The design parameters of the macro geometry mainly include the roller diameter , number of rollers , bearing pitch diameter , and radial clearance P . In addition, the radius of the crankshaft distribution circle ( ) determines the installation position of the crankshaft in the cycloidal gear and the force of the crankshaft bearing.
When the rolling elements are only lightly profiled or subjected to excessive loads, edge stresses can arise. Cylindrical rollers (or raceways) are usually crowned to solve this problem [19]. Fujiwara [21,22] proposed a modified logarithmic shape that allowed for a flat region on the roller, as shown in Figure 5. The crown drop in the axial direction at position is where is the maximum load acting on the rollers, is the load coefficient, is the coefficient of the crown length, and is the crown drop at the edge of the roller end. Therefore, there are nine mutually independent design variables. , , , , , and define the internal geometric dimensions, and , , and define the profile of the roller. Therefore, there are nine mutually independent design variables. r o , D we , D m , L we , Z, and P d define the internal geometric dimensions, and K 1 , K 2 , and Z m define the profile of the roller.

Objective Function
Crankshaft bearings are the key components limiting the service life of RV reducers. To improve the service life of the reducer, fatigue life of crankshaft bearings is taken as the objective function: where X denotes the vector of design variable, and f (X) is the objective function. The calculation method of L 10r is derived from Section 2.

Constraint Conditions
The design of crankshaft bearing must satisfy various constraints. In practice, the geometry of the crankshaft bearing is limited by the cycloidal gear. The crankshaft bearing should have sufficient load capacity. The crankshaft bearing should be well-lubricated to reduce wear and friction. The assembly conditions of RV reducer must be satisfied. In this section, the constraints conditions are derived from the geometry, lubrication, strength of crankshaft bearings, and structure of the RV reducer.

Constraints on Bearing Pitch Diameter
The outer raceway is located between the root circle and the center hole of the cycloidal gear. Therefore, the bearing pitch diameter should be smaller than the difference between the diameter of the center hole and the diameter of the root circle. The inner raceway of the bearing is integrated in the off-center section of the crankshaft. The eccentric distance of the crankshaft is e. Therefore,

Constraints on Roller Diameter
Geometric constraints of the cycloidal gear determine the upper bound of the roller diameter:

Constraints on Roller Length
For cylindrical roller bearings, the length of the roller should not be longer than 2.5 times the diameter of the rollers [23]. As shown in Figure 2, the roller length is not greater than the width of the crankshaft.
3.3.4. Constraints on Crowned Roller Profile K 1 , K 2 , and Z m , are the parameters that control the profile of rollers. K 1 is the load coefficient, K 2 is the coefficient of the crown length, and Z m is the crown drop at the edge of the roller end. The ranges of the parameters for the logarithmic function are given by [22]

Constraints on Smooth Operation of the Bearing
To guarantee smooth operation of the rollers without collision, an appropriate spacing should be provided between the rollers. The minimum angular spacing between the rollers is assumed to be greater than one in the crankshaft bearing [23].

Constraints on Contact Strength
The contact force is concentrated on a small contact area between the raceway and the roller. The contact pressure is very high, resulting in permanent deformation. Experience has shown that the permanent deformation is less than 0.0001 times the roller diameter to ensure safe operation of the bearing. Alternatively, the maximum contact pressure should be less than 4000 MPa [24]. The contact stress is defined in Section 2.3.2. Thus,

Lubricant Film Thickness Constraints
Proper lubrication prevents severe friction between rollers and raceways. It can avoid metal-to-metal contact between raceway and roller. According to Dowson and Higginson [25], the minimum film thickness for inner and outer raceways is where ζ is the pressure-viscosity coefficient, µ is viscosity of oil, U out is the sliding velocity between the outer raceway and the roller, U in is the sliding velocity between the inner raceway and the roller, and Q max is the maximum load on the rollers. To assess the lubrication appropriateness, the film parameter (Λ) is defined [19] as follows: where h out min and h in min are the minimum lubricant film thickness for the outer and inner raceway, respectively. s R and s m is the root mean square roughness of the roller surface and the raceway surface, respectively. In this study, the value of s m and s R is 1.25 µm.
To ensure sufficient lubrication between the roller raceways, the film parameter value should be less than 3 [6].

Constraints on Structure of RV Reducer
For cylindrical roller bearings, the thickness of the bearing ring should not be less than εD we , where ε is between 0.3 and 0.4 [23,26]. The outer ring is integrated into the cycloidal gear. The thickness between the outer and center holes and the root circle is greater than εD we .
To ensure the strength of the cycloidal gear, the value of ε is 0.4. As shown in Figure 1, the planetary gears are fixed with crankshafts. The radius of the crankshaft bearing distribution circle is equal to the center distance of the sun gear and the planetary gear.
where z 1 , z 2 are the number of teeth of the sun gear and planetary gear, respectively, and m is the gear module. The sun gear meshes with n planet wheels.
According to the theory of planetary transmission, to ensure the assembly conditions of multiple planetary gears, the sum of the number of teeth of the sun and planetary gears should be an integer multiple of the number of planetary gears, Therefore, the radius of the crankshaft bearing distribution circle should be an integer multiple of the product of the number of planetary gears and the gear module:

Crow Search Algorithm
In this study, the optimization of the crankshaft bearing is a typical constraint nonlinear optimization problem. Metaheuristic algorithms have shown promising performance for solving real-world optimization problems that are extremely nonlinear. The crow search algorithm (CSA) is a novel metaheuristic optimization method proposed by Askarzadeh [27]. The CSA has demonstrated its effectiveness and ability to solve many complex engineering optimization problems, such as environmental dispatch, renewable energy, and distribution networks [28]. Therefore, CSA is adopted to solve the optimization problem.
CSA is a population-based metaheuristic algorithm. It is based on the intelligent behavior of crows and searches for the optimal solution through iteration. CSA is initialized with a certain number of crows (flock size) with a random position. The position of each crow represents a solution to the optimization problem. A vector is used to define the position of the crow. The vector dimension represents the number of design variables. During the iteration, the position of each crow is updated. There are two parameters to control the generating of the position. They are flight length and awareness probability respectively. The best solution of each crow is memorized during the iteration. The optimization process stops when the maximum number of iterations is reached. The best position of the crows is select as the optimal solution to the optimization problem.
In Section 3.3, constraints on geometry, lubrication, contact strength of crankshaft bearings, and structure of the RV reducer are defined. To deal with the constraints, the penalty method is adopted in this study. The constrained problem is transformed into an unconstrained problem by the penalty function.
The penalty function is defined as where τ i 1 and υ j ≥ 0. These penalty function parameters should be sufficiently large for this application, but the values depend on the required quality of the solution. M is the number of inequality constraints, and N is the number of equality constraints.
H i [g i (X)] and F j ξ j (X) are index functions: CSA is adopted to solve the optimization problem. Figure 6 shows the flowchart of the proposed optimization design of the bearing. First, the objective function, design variables, and constraints are defined. Second, the adjustable parameters of the CSA are specified. The parameter of the RV reducer and initial parameters of the bearing are given. The position and memory of the crows is initialized. Third, a new position is generated, and the feasibility of the solution is checked through constraints. The process is repeated until the maximum number of iterations is reached.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 13 of 23 In Section 3.3, constraints on geometry, lubrication, contact strength of crankshaft bearings, and structure of the RV reducer are defined. To deal with the constraints, the penalty method is adopted in this study. The constrained problem is transformed into an unconstrained problem by the penalty function.
The penalty function is defined as where ≫ 1 and ≥ 0. These penalty function parameters should be sufficiently large for this application, but the values depend on the required quality of the solution. is the number of inequality constraints, and is the number of equality constraints.
CSA is adopted to solve the optimization problem. Figure 6 shows the flowchart of the proposed optimization design of the bearing. First, the objective function, design variables, and constraints are defined. Second, the adjustable parameters of the CSA are specified. The parameter of the RV reducer and initial parameters of the bearing are given. The position and memory of the crows is initialized. Third, a new position is generated, and the feasibility of the solution is checked through constraints. The process is repeated until the maximum number of iterations is reached.

Optimization Example
The proposed optimization method is applied to the RV 20E and RV 110E-type reducers. The macro geometry and crowned roller profile of crankshaft bearings are optimized. The features of the reducers are listed in Table 1. The grease used is MolywhiteRE00. The output speed of the reducer is 10 rev/min. The optimization is implemented with MATLAB 2017a on a Windows 10 Intel(R) Core(TM) i7-7700 CPU @3.6GHz with 8 GB RAM.

Optimization Design of Crankshaft Bearing of RV-20E
The convergence history of the fatigue life of crankshaft bearing of RV-20E is shown in Figure 7.

Optimization Example
The proposed optimization method is applied to the RV 20E and RV 110E-type reducers. The macro geometry and crowned roller profile of crankshaft bearings are optimized. The features of the reducers are listed in Table 1. The grease used is MolywhiteRE00. The output speed of the reducer is 10 rev/min. The optimization is implemented with MATLAB 2017a on a Windows 10 Intel(R) Core(TM) i7-7700 CPU @3.6GHz with 8 GB RAM.

Optimization Design of Crankshaft Bearing of RV-20E
The convergence history of the fatigue life of crankshaft bearing of RV-20E is shown in Figure 7.   Table 2 shows the comparison of the parameters before and after optimization. Before optimization, the profile of the roller was flat. Therefore, the values of , , and are set to 0. In practical situations, the optimization results must be rounded. The rounded parameters and life are also listed in the table. By substituting the optimization variables into the above constraints, it can be known that all optimization variables can satisfy the constraints. Before optimization, the fatigue life of the crankshaft bearing was 4475.4 million revolutions. After optimization, the bearing life is  Table 2 shows the comparison of the parameters before and after optimization. Before optimization, the profile of the roller was flat. Therefore, the values of k 1 , k 2 , and Z m are set to 0. In practical situations, the optimization results must be rounded. The rounded parameters and life are also listed in the table. By substituting the optimization variables into the above constraints, it can be known that all optimization variables can satisfy the constraints. Before optimization, the fatigue life of the crankshaft bearing was 4475.4 million revolutions. After optimization, the bearing life is increased to 13,254.2 million revolutions. The fatigue life after optimization is increased by 196% compared with the life before optimization. After optimization, the value of the clearance is negative. This is because a proper negative clearance increases the bearing life. The structure of the crankshaft bearing of RV-20E before and after optimization is shown in Figure 8. The space of the cycloid gear is better used after optimization. The structural parameters of the cycloidal gear have not changed during the design optimization. The diameter of outer raceway is increased, while the diameter of the inner raceway is decreased. The roller diameter is increased, but the number of rollers is reduced. This makes the internal structure of the crankshaft bearing more reasonable. In addition, the radius of the crankshaft bearing distribution circle is increased. This is because the load of the crankshaft bearing from the cycloidal gear decreases as the radius of the bearing distribution circle increases.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 15 of 23 increased to 13,254.2 million revolutions. The fatigue life after optimization is increased by 196% compared with the life before optimization. After optimization, the value of the clearance is negative. This is because a proper negative clearance increases the bearing life. The structure of the crankshaft bearing of RV-20E before and after optimization is shown in Figure 8. The space of the cycloid gear is better used after optimization. The structural parameters of the cycloidal gear have not changed during the design optimization. The diameter of outer raceway is increased, while the diameter of the inner raceway is decreased. The roller diameter is increased, but the number of rollers is reduced. This makes the internal structure of the crankshaft bearing more reasonable. In addition, the radius of the crankshaft bearing distribution circle is increased. This is because the load of the crankshaft bearing from the cycloidal gear decreases as the radius of the bearing distribution circle increases. When the minimum thickness of cycloidal gear is too small, the structure strength of the cycloidal gear may be weakened. However, the constraint on the thickness of the cycloidal gear is considered in Section 3.3.8. Therefore, the thickness of the cycloidal gear after optimization will not be too thin. When the inner raceway diameter is too small, the contact stress between the inner raceway and bearings may exceed the strain limit. The constraints on contact strength have been When the minimum thickness of cycloidal gear is too small, the structure strength of the cycloidal gear may be weakened. However, the constraint on the thickness of the cycloidal gear is considered in Section 3.3.8. Therefore, the thickness of the cycloidal gear after optimization will not be too thin. When the inner raceway diameter is too small, the contact stress between the inner raceway and bearings may exceed the strain limit. The constraints on contact strength have been proposed in Section 3.3.6, too. Therefore, the strength of the cycloidal gear and crankshaft bearing is not reduced after optimization.
The optimization roller profile of crankshaft bearing of RV 20E is shown in Figure 9. The roller profile is partially crowned. The crown drop at roller end is 10 µm. The ratio of crowning length to effective contact length is 0.7. The flat region increases the bearing capacity of the roller-raceway. The crown drop eliminates edge stress concentration.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 16 of 23 proposed in Section 3.3.6, too. Therefore, the strength of the cycloidal gear and crankshaft bearing is not reduced after optimization. The optimization roller profile of crankshaft bearing of RV 20E is shown in Figure 9. The roller profile is partially crowned. The crown drop at roller end is 10 μm. The ratio of crowning length to effective contact length is 0.7. The flat region increases the bearing capacity of the roller-raceway. The crown drop eliminates edge stress concentration. The contact stress on the roller that bears the maximum load is shown in Figure 10. The contact stress on the contact region was calculated using the non-Hertzian contact method. Before optimization, the contact stress in the middle of the roller was approximately 1400 MPa. An extremely large edge stress of 2190 MPa occurred at the end of roller. After optimization, the roller profile suitable for the work condition is obtained. The stress concentration at roller end is eliminated. The maximum contact stress of the roller is only 1245 MPa. As the roller diameter is larger, the contact stress between the raceway and roller decreases overall, and the stress in the middle of the roller decreases to 1150 MPa. The contact stress on the roller that bears the maximum load is shown in Figure 10. The contact stress on the contact region was calculated using the non-Hertzian contact method. Before optimization, the contact stress in the middle of the roller was approximately 1400 MPa. An extremely large edge stress of 2190 MPa occurred at the end of roller. After optimization, the roller profile suitable for the work condition is obtained. The stress concentration at roller end is eliminated. The maximum contact stress of the roller is only 1245 MPa. As the roller diameter is larger, the contact stress between the raceway and roller decreases overall, and the stress in the middle of the roller decreases to 1150 MPa.

Optimization Design of Crankshaft Bearing of RV-110E
The convergence history of the bearing fatigue life of RV-110E is shown in Figure 11. After approximately 4 h of running time, optimization results satisfying all the design constraints are obtained. The fatigue life increases to 4840.1 million revolutions after 76 cycles. After this, there are small increases in the fatigue life, which increases to 5064.8 million revolutions after 196 cycles.

Optimization Design of Crankshaft Bearing of RV-110E
The convergence history of the bearing fatigue life of RV-110E is shown in Figure 11. A comparison of the parameters of RV-110E before and after optimization is presented in Table  3. By substituting the optimization variables into the above constraints, it can be known that all optimization variables can satisfy the constraints. Before optimization, the fatigue life of the crankshaft bearing was 1912.2 million revolutions. After optimization, the bearing life is increased to 5125.3 million revolutions. The fatigue life after optimization is increased by 168% compared with that before optimization.  Figure 12 shows the structure of the crankshaft bearing of RV-110E before and after optimization. For RV-110E, the change in the structure of the crankshaft bearings is similar to that of RV-20E. The outer raceway diameter and roller diameter are also increased. The inner raceway diameter is decreased, and the number of rollers is reduced. The radius of the crankshaft bearing distribution circle is also increased. Although these dimensions have been changed, due to the constraints, the amount of change is within a limited range. Therefore, the change will not weaken the cycloidal gear and crankshaft. A comparison of the parameters of RV-110E before and after optimization is presented in Table 3. By substituting the optimization variables into the above constraints, it can be known that all optimization variables can satisfy the constraints. Before optimization, the fatigue life of the crankshaft bearing was 1912.2 million revolutions. After optimization, the bearing life is increased to 5125.3 million revolutions. The fatigue life after optimization is increased by 168% compared with that before optimization.  Figure 12 shows the structure of the crankshaft bearing of RV-110E before and after optimization. For RV-110E, the change in the structure of the crankshaft bearings is similar to that of RV-20E. The outer raceway diameter and roller diameter are also increased. The inner raceway diameter is decreased, and the number of rollers is reduced. The radius of the crankshaft bearing distribution circle is also increased. Although these dimensions have been changed, due to the constraints, the amount of change is within a limited range. Therefore, the change will not weaken the cycloidal gear and crankshaft.  Figure 13 shows the crown drops of crankshaft bearing rollers of RV-110E. The crown drop at the end of the roller is 12 μm. The ratio of the crowning length to the effective contact length is 0.63. The roller diameter is increased after optimization, which helps to reduce the contact stress between the raceway and the roller.  Figure 13 shows the crown drops of crankshaft bearing rollers of RV-110E. The crown drop at the end of the roller is 12 µm. The ratio of the crowning length to the effective contact length is 0.63. The roller diameter is increased after optimization, which helps to reduce the contact stress between the raceway and the roller. The contact stress on the roller that bears the maximum load is shown in Figure 14. Before optimization, the flat profile of the roller resulted in edge stress concentrations at the roller ends. The stress in the middle of the roller was approximately 1500 MPa. The maximum stress was 2660 MPa at the end of the roller. After optimization, the edge stress concentration is effectively eliminated. The maximum stress on the roller is reduced to 1380 MPa. The stress in the middle of the roller is reduced to 1300 MPa. The distribution of the contact stress over the length of the roller is much more uniform. This is beneficial for increasing the bearing life. The contact stress on the roller that bears the maximum load is shown in Figure 14. Before optimization, the flat profile of the roller resulted in edge stress concentrations at the roller ends. The stress in the middle of the roller was approximately 1500 MPa. The maximum stress was 2660 MPa at the end of the roller. After optimization, the edge stress concentration is effectively eliminated. The maximum stress on the roller is reduced to 1380 MPa. The stress in the middle of the roller is reduced to 1300 MPa. The distribution of the contact stress over the length of the roller is much more uniform. This is beneficial for increasing the bearing life.

Conclusions
This study proposes a novel optimization method for RV reducer crankshaft bearings considering the bearing geometry and crowned roller profile. The working load of the bearing and the structure of the reducer are considered. Constraints derived from the geometry, lubrication, strength of crankshaft bearings, and structure of the RV reducer are established. The fatigue life of the bearing is optimized by the proposed optimization method. The research summary of conclusions is as follows: (1) A comprehensive optimization design model of an RV reducer crankshaft bearing is constructed. Based on the CSA algorithm, a nonlinear constraint optimization program is developed to solve the optimization problem. The convergence process shows that the optimization program has high efficiency and ability to solve the optimization problem of RV reducer crankshaft bearings. (2) The roller profile and macro-geometric parameters of the bearing are both considered in this study. Through the optimization design, the optimum crowned roller profile suitable for the working load of the bearing is obtained, and the stress concentration between the roller and the raceway is eliminated. Compared with the situation before optimization, the limited geometric space of the cycloidal gear is used more effectively. (3) Taking the RV-20E and RV-110E reducer crankshaft bearings as examples, the two bearings are optimized by the proposed method. After optimization, all the optimization variables satisfy the constraint conditions. The RV-20E bearing life is increased by 196%, whereas that of the RV110 is increased by 168%.
The optimization design method proposed in this study has an important guiding significance for improving the service life of RV reducers. It can thus provide a theoretical basis for the design of RV reducer crankshaft bearings. The optimization method developed in this study can also be used in the design of other nonstandard roller bearings

Conclusions
This study proposes a novel optimization method for RV reducer crankshaft bearings considering the bearing geometry and crowned roller profile. The working load of the bearing and the structure of the reducer are considered. Constraints derived from the geometry, lubrication, strength of crankshaft bearings, and structure of the RV reducer are established. The fatigue life of the bearing is optimized by the proposed optimization method. The research summary of conclusions is as follows: (1) A comprehensive optimization design model of an RV reducer crankshaft bearing is constructed.
Based on the CSA algorithm, a nonlinear constraint optimization program is developed to solve the optimization problem. The convergence process shows that the optimization program has high efficiency and ability to solve the optimization problem of RV reducer crankshaft bearings. (2) The roller profile and macro-geometric parameters of the bearing are both considered in this study. Through the optimization design, the optimum crowned roller profile suitable for the working load of the bearing is obtained, and the stress concentration between the roller and the raceway is eliminated. Compared with the situation before optimization, the limited geometric space of the cycloidal gear is used more effectively. (3) Taking the RV-20E and RV-110E reducer crankshaft bearings as examples, the two bearings are optimized by the proposed method. After optimization, all the optimization variables satisfy the constraint conditions. The RV-20E bearing life is increased by 196%, whereas that of the RV110 is increased by 168%.
The optimization design method proposed in this study has an important guiding significance for improving the service life of RV reducers. It can thus provide a theoretical basis for the design of RV reducer crankshaft bearings. The optimization method developed in this study can also be used in the design of other nonstandard roller bearings