1. 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 backlash of RV reducers. Hsieh [
10] established an analytical model of system dynamics in cycloidal pin gear reducers with small tooth differences. Jin [
11] analyzed the transmission error of the RV reducer by using virtual machine technology. Wang [
12] proposed a new contact model of cycloidal meshing pair and analyzed the transmission accuracy of RV reducer. The influence of tooth profile modification was analyzed. Li [
13] analyzed the meshing characteristics of RV reducers considering the manufacturing error. Chen [
14] analyzed the vibration characteristics of the RV reducer. Xu [
4] established a dynamic model of the RV reducers. The transmission mechanisms of cycloidal gear were analyzed. Wei [
15] analyzed the lubrication characteristics of RV reducers. Zhang [
5] analyzed the mixed lubrication of crankshaft bearings of an RV reducer. Wang [
16] used genetic algorithm to optimize the design of cycloid reducer. Chu [
17] proposed a selective assembly method for RV reducers using genetic algorithms. Wang [
18] proposed a multi-objective optimization method for cycloidal speed reducer. The objective functions are the reducer volume, loads on the crankshaft bearings, and stress of cycloidal gear.
The above review shows that most of the published works focus on the meshing characteristics of the cycloidal pair, tooth profile modification, and performance characteristics of the RV reducer. However, few scholars have focused on optimizing the fatigue life of the reducers. Although studies [
16,
18] have optimized the fatigue life of cycloidal reducers by using a multi-objective optimization method, the geometric parameters and roller profile of the bearing have not been comprehensively considered during the optimization design.
Therefore, an optimization method for RV reducer crankshaft bearings considering the bearing geometry and crowned roller profile is presented. First, the performance of the bearing considering the load working condition and roller profile is analyzed. Second, a mathematical model for the optimization design is developed. Third, two types of crankshaft bearings are optimized by the proposed optimization method. A comparison of the parameters of the bearing before and after optimization is presented. Finally, a research summary of the study is provided.
2. Performance of Crankshaft Bearings of Rotate Vector (RV) Reducers
Crankshaft bearings are the key components restricting the service life of RV reducers. The performance of the bearing is analyzed in this section. First, the structure of the RV reducers is analyzed. Second, the loads on the crankshaft bearing in the reducer are calculated. Finally, the fatigue life of the crankshaft bearing considering the clearance and crowned roller profile is evaluated.
2.1. Structure of RV Reducer
The RV reducer is a two-stage planetary transmission mechanism. The first stage is involute gear transmission. The second stage is cycloid gear transmission.
Figure 1 is the structure of the RV reducer. In a cycloidal–pinwheel transmission, the crankshaft and cycloidal gear are connected by the crankshaft bearing. The off-center section of the crankshafts causes the cycloidal gear to rotate in eccentricity. The cycloidal gear meshes with the pin gear. Therefore, the cycloid gear also rotates around its own axis. The crankshaft transmits the spinning motion of the cycloidal gear to the output mechanism of the reducer.
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.
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
. As shown in
Figure 2, the boundary dimensions include the outer diameter (
), minimum diameter of crankshaft (
), bearing width (
), diameter of the center hole (
), and the root circle of the cycloidal gear (
).
2.2. 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 the
-axis 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
,
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,
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,
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
. 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
2.3. 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.
2.3.1. 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.
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
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).
2.3.2. 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
, it is assumed that the contact area between the raceway and the roller is
. The contact pressure at coordinate
is
. The integral of pressure
is equal to the contact force
.
Contact stress occurs when two surfaces come into contact with each other. Assuming that the contact deformation between the two surfaces at coordinate
is
, then
According to the half-space theory, when two bodies come into contact, the elastic deformation
can be expressed by the following equation:
where
is the equivalent Young’s modulus.
where
is the Poisson’s ratio and
is the elastic modulus.
When the roller is subjected to the normal load
, the contact surfaces are moved toward each other by a displacement
. For a roller with a crowned profile,
is the gap between the roller and raceway at the coordinate
generated by the crown drop.
is the initial distance generated by the geometry of the roller and raceway. Therefore, Equation (17) can be written as follows:
In these equations, the relative displacement
, pressure distribution
, and contact area
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
is the influence coefficient representing the deflection of cell
because of a uniform pressure over the cell
, which can be obtained according to reference [
20].
is the contact deformation of cell
.
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.
2.3.3. Fatigue Life
The basic dynamic load rating of the raceways
are obtained by
According to ISO/TR 1281-1, the value of
is 0.83.
is the parameter related to the structure of the bearing,
. The basic dynamic load capacity
of the bearings is
where
accommodates recent improvements in bearing quality and standard ISO/TR 1281-1 gives values of the factor
.
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
of slices. The basic load capacity of a slice of the inner and outer rings,
, is
The dynamic equivalent load on the
th slice is:
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
3. 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:
3.1. 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 number of rollers , bearing pitch diameter , and radial clearance . 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.
3.2. 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
denotes the vector of design variable, and
is the objective function. The calculation method of
is derived from
Section 2.
3.3. 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.
3.3.1. 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
. Therefore,
3.3.2. Constraints on Roller Diameter
Geometric constraints of the cycloidal gear determine the upper bound of the roller diameter:
3.3.3. 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
,
, and
, are the parameters that control the profile of rollers.
is the load coefficient,
is the coefficient of the crown length, and
is the crown drop at the edge of the roller end. The ranges of the parameters for the logarithmic function are given by [
22]
3.3.5. 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].
3.3.6. 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
[
24]. The contact stress is defined in
Section 2.3.2. Thus,
3.3.7. 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,
is the sliding velocity between the outer raceway and the roller,
is the sliding velocity between the inner raceway and the roller, and
is the maximum load on the rollers.
To assess the lubrication appropriateness, the film parameter (Λ) is defined [
19] as follows:
where
and
are the minimum lubricant film thickness for the outer and inner raceway, respectively.
and
is the root mean square roughness of the roller surface and the raceway surface, respectively. In this study, the value of
and
is 1.25 μm.
To ensure sufficient lubrication between the roller raceways, the film parameter value should be less than 3 [
6].
3.3.8. Constraints on Structure of RV Reducer
For cylindrical roller bearings, the thickness of the bearing ring should not be less than
, 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
.
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
are the number of teeth of the sun gear and planetary gear, respectively, and
is the gear module. The sun gear meshes with
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:
3.4. 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
and
. 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.
and
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.
5. 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