Open Access
Issue
Mechanics & Industry
Volume 22, 2021
Article Number 34
Number of page(s) 11
DOI https://doi.org/10.1051/meca/2021030
Published online 06 May 2021

© Y. Yu et al., Published by EDP Sciences 2021

Licence Creative CommonsThis is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1 Introduction

Harmonic reducer is widely used in aerospace, medical devices, robots, CNC machine tools, packaging equipment, instruments and other fields. As one of the main working components, the flexible bearing (FB)'s damage is one of the main failure reasons of harmonic drive [13]. It is necessary to study the load distribution in flexible bearing for high efficiency, high precision and long life of harmonic drive.

The harmonic reducer is mainly formed of three components: circular spline (CS), flexspline (FS) and wave generator (WG), where the WG consists of a FB and a wave generator cam (WG cam), shown in Figure 1. When the CS is fixed and the WG is the driving member, the FS becomes the follower. The WG enables the FS to generate a controllable elastic deformation, it forces that the teeth at both ends of the major axis of the CS and the FS are in full engagement, while the teeth at both ends of the minor axis are completely detached. When the WG rotates, the teeth of the FS mesh with the teeth of the CS in turn, repeating four states: engaging-in, engaging, engaging-out and disengaging, shown in Figure 1. Thus, this staggered tooth transmission makes the harmonic reducer drive with great deceleration ratio.

The FB is of great difference compared with the ordinary bearing. In structure, the thickness of FB is much smaller than that of ordinary bearing. After the FB is assembled on the WG cam, the inner and outer rings will have a certain deflection, and the shape depends on that of the WG cam. In loading form, when the harmonic reducer works, the flexible bearing will bear symmetrical external load at both ends of the major axis. Even there is no external load, the internal load distribution of FB is quite different from that of ordinary bearing due to the pre-deformation caused by WG cam.

The calculation and research of load distribution in rolling bearing is the analysis basis of kinematics, lubrication mechanics, contact fatigue life and efficiency of rolling bearing, which is necessary for improving the life and the working performance of rolling bearing. For the ordinary rolling bearing, the analysis and calculation of its load distribution has become the mature theory [4,5]. However, the deformation, load and failure modes of the FB are quite different from those of ordinary bearing. Thus, it isn't accurate to study the FB by the existing theory. In order to research the mechanical characteristics of the FB profoundly, the primary step is to obtain its exact internal load distribution.

Several theoretical methods were applied for studying the load distribution in the FB. Le and Shen [6] obtained the load distribution of the FB using the curved beam theory of elasticity and material mechanics, three-moment theory of continuous beam, and energy method theory. The load decomposition obtained by the model is in good agreement with the experimental results. But this method can only calculate the bearing load distribution of the FB mounted on a four-force action type cam. The deflection of the outer ring of the FB was calculated based on the circular radial deformation under a static equilibrium state and Shao and Wen [7] obtained the load on the rolling element by solving a set of nonlinear deformation coordination equations. However, only the case of odd number balls was discussed in the model. A single-force thin-walled ring superposition model was established to calculate the load distribution of the FB fitted on a cam [8,9]. This method could be used to calculate the load distribution in the FB with random number of balls, but the external loads are not considered in this model. And a two-force thin-walled ring super-position model was established to calculate the load distribution of the FB fitted on a cam [10]. However, this method could only be applied for the FB with even number of balls. Therefore, it is necessary to provide a new theoretical method for calculating the load distribution in the FB with random number of balls, which is in force balance at the same time.

According to the general solution of the deformation of the thin-walled ring deduced by Liu and Chiu [11], Xiong et al. [12] developed a universal static analysis model to analyze the load distribution of different FBs in the harmonic reducer. Not only arbitrary numbers of balls (even or odd numbers) can be considered in this model but also arbitrary radial symmetric external loads.

Apart from theoretical analysis, finite element methods (FEMs) were also employed in the load analysis of the FB [13,14]. Compared with the theoretical methods, this kind of method requires a lot of computing time.

In this paper, a superposition model of three forces ring was proposed based on the theory of thin-walled ring and superposition principle. According to this mechanical model, the internal load distribution in the FB could be calculated. And then the radial deformation and the bending normal stress on the outer ring were also calculated and analyzed. Based on the calculation of this mechanical model, we investigated the influences of the number of balls, the angle of the ball position and the load torque on the ball load. And the influences of the numbers of balls and the load torque on the radial deformation and maximum bending normal stress of outer ring were subsequently studied.

thumbnail Fig. 1

Schematic diagram of the working principle of the harmonic reducer.

2 A mechanical model of three-force ring

2.1 The basic theory

A ring whose cross-section size is much smaller than its radius is called a thin-walled ring. It is assumed that the cross-section of the ring is rectangular and the shape is unchanged along the circumference. And the load is uniformly distributed along the width of the ring. Under these conditions, stress and displacement are constant along the width, so it can be treated as a plane problem. For the simplification of calculation, the mechanical model has the following basic assumptions: ①It is assumed that the load distribution inside the FB is equal to that the outer ring bears multiple radial loads, as shown in Figure 2. ②It is assumed that the deformation of the outer ring is small, so the methods of material mechanics and superposition principle could be used to study the load distribution of outer ring.

In this paper, three main formulas about the theory of thin-walled ring [15,16] are given by: d 2 w d ϕ 2 + w = M R 2 E I (1) w d ϕ + d v = 0 (2) θ = ( 1 R ) ( w d v d ϕ ) (3)

Equation (1) is the elastic equation connecting bending moment M and radial displacement w, where E is the elastic modulus of the material and I is the inertial moment of the ring section. Equation (2) is derived from no elongation hypothesis, and it indicates the relationship between radial displacement w and circumferential displacement v. Equation (3) is the expression of normal angle θ of circular section.

The deformation of outer ring of flexible bearing can be analyzed by the theory of thin-walled ring because of its small thick-ness. In the following analysis and calculations, the shape of the ring is indicated by the neutral layer curve of the ring, and the load and deformation of points on the neutral layer are discussed.

thumbnail Fig. 2

Load simplification of outer ring of the FB.

2.2 Distributed bending moment MAB(ϕ) on ring

The structure diagram of the three-force ring model is shown in Figure 3. If F1 = F is supposed, then other two radial forces F2 and F3 can be known according to the equilibrium relationship of the force system. The ring is cut along the horizontal midline and the upper part of the ring is taken for analysis, as shown in Figure 4.

When the bending moment MA, MB and radial force system act respectively on the half ring shown in Figure 4, the normal angles of section A and B could be calculated by the unit force method [17]. Then the total normal angles of two sections are obtained by superposition. According to deformation compatibility condition, the values of these two normal angles are 0, so the bending moments are given by: M A = F R ( h 2 3 h 1 ) / ( 8 π ) M B = F R ( h 1 3 h 2 ) / ( 8 π ) (4)where h1 and h2 are given by: h 1 = 2 2cos ϕ 0 + ϕ 0 tan ϕ 0 2sin ϕ 0 tan ϕ 0 + πtan ϕ 0 / 2 h 2 = 2 2cos ϕ 0 ϕ 0 tan ϕ 0 2sin ϕ 0 tan ϕ 0 +3 πtan ϕ 0 / 2 (5)

In the same way, when the bending moment MA, MB and radial force system act together on the half ring, the total distributed bending moment on ring could be calculated by superposition. It is given by: M B C ( ϕ ) = F R ( h 3 h 4 cosϕ sinϕ ) / 2 0 ϕ ϕ 0 M A C ( ϕ ) = F R ( h 3 h 5 cosϕ ) / 2 ϕ 0 ϕ π (6)where h3, h4 and h5 are given by: h 3 = ( h 1 + h 2 ) / ( 4 π ) tan ϕ 0 / 2 h 4 = ( h 1 h 2 ) / ( 2 π ) tan ϕ 0 / 2 h 5 = ( h 1 h 2 ) / ( 2 π ) +tan ϕ 0 / 2 (7)

thumbnail Fig. 3

Structure diagram of the ring under the action of three forces.

thumbnail Fig. 4

Internal force of the symmetric section of the ring.

2.3 The expression of radial displacement

Substituting equation (6) into equation (1), differential equations are given by: d 2 w B C d ϕ 2 + w B C = F R 3 2 E I ( h 3 h 4 cosϕ sinϕ ) 0 ϕ ϕ 0 d 2 w A C d ϕ 2 + w A C = F R 3 2 E I ( h 3 h 5 cosϕ ) ϕ 0 ϕ π (8)

The general solution of two differential equations is expressed as: w B C = A 1 sinϕ + A 2 cosϕ F R 3 2 E I ( h 3 1 2 h 4 ϕsinϕ + 1 2 ϕcosϕ ) 0 ϕ ϕ 0 (9) w A C = B 1 sinϕ + B 2 cosϕ F R 3 2 E I ( h 3 1 2 h 5 ϕsinϕ ) ϕ 0 ϕ π (10)where A1, A2, B1 and B2 are integral constants.

Because of the symmetry of structure, the normal angle θ and the circumferential displacement v equal to 0 when ϕ = 0 and ϕ = π. Substituting this symmetric condition into equation (3), A1 and B1 are given by: A 1 = F R 3 4 E I , B 1 = F R 3 4 E I ϕ 0 tan ϕ 0 (11)

Then according to the condition of smooth connection, that is when ϕ = ϕ0, equations (9) and (10) should get the same values and derivatives. The same result is given by: B 2 = ( A 1 B 1 )tan ϕ 0 + A 2 F R 3 ϕ 0 4 E I cos 2 ϕ 0 (12)

It seems that the above conditions are not enough to get a solution, this is because of the lack of horizontal constraints for the structure in Figure 3. To completely determine the four integration constants, we let the radial displacement at point B equal to x0 in Figure 3. According to equations (9), (11) and (12), A2 and B2 denoted by x0 can be given by: A 2 = F R 3 2 E I h 3 + x 0 , B 2 = F R 3 2 E I ( h 3 + 1 2 tan ϕ 0 1 2 ϕ 0 ) + x 0 (13)

The expression of x0 needs to be defined. Comparing the mechanical models of two-force ring and three-force ring, the circumferential displacement at ϕ = π/2 on the ring should be 0 to ensure that the displacement of point B relative to the vertical line passing through the center of the circle is constant. And through equation (2), the circumferential displacement of the ring can be obtained by integrating the radial displacement. By integral calculation, the expression of x0 is given by: x 0 = F R 3 2 E I [ 1 2 ( ϕ 0 tan ϕ 0 h 5 ) ( π 2 + 1 ) h 3 ] (14)

Finally, substituting equations (11), (13) and (14) into equations (9) and (10), the radial displacement of three-force ring can be given by: w A B = F H ( ϕ ) 0 ϕ π (15)where H (ϕ) is the coefficient of bending deformation, which is given by: H(ϕ)={R32EI[12(sinϕ+h4ϕsinϕϕcosϕ)+12(ϕ0tanϕ0h5πh3)cosϕh3]0ϕϕ0R32EI[12ϕ0tanϕ0sinϕ12(h5+πh3)cosϕ+12h5ϕsinϕh3]ϕ0ϕπ(16)

3 Load distribution in the flexible bearing

3.1 Calculation parameters

When the FS is loaded, the loading diagram of outer ring of the FB with N balls (N ≧ 5) is shown in Figure 5. For the simplification of calculation, we have the following assumptions: ①It is assumed that there is no friction between the CS and the FS. ②It is assumed that the rotational speed of the shaft has no influence on the radial force state of the FB. ③It is assumed that the centrifugal load of the balls on the outer ring is much smaller than the ball load itself, which can be ignored. The parameters related to external radial distributed load are given in Table 1, and the main geometric parameters of flexible bearing are given in Table 2. The values of Φ1, Φ2 and Φ3 determine the location of the radial distributed load. In this paper, the values of these three variables refer to a relevant paper [18]. The main variables are the load torque and the number of balls, which can be changed in a certain range.

thumbnail Fig. 5

Loading diagram of outer ring of the FB with N balls (N ≧ 5).

Table 1

Relevant Parameters of radial distributed load.

Table 2

Main parameters of the FB.

3.2 The deformation compatibility equation

In this paper, the cosine-cam is used for analysis of internal load of the FB. The ideal pre-deformation caused by this kind of cam is given by: w c a m = w 0 cos( 2 ϕ ) (17)

when Φ2=Φ3, the radial distributed load in Figure 5 is given by [19,20]: q r = q t max tan τcos[ ( ϕ Φ 1 ) π / ( 2 Φ 2 ) ] q t max = π T / ( 2 Φ 2 d r 2 b w ) (18)where qtmax is the maximal value of tangential distributed load; τ is profile angle; dr is the pitch diameter of the FS; bw is the working width of gear ring.

Under the radial distributed load, the radial deformation of the outer ring is given by [15]: w q = B o R 4 / E I k = 2,4,6 q r k / ( k 2 1 ) 2 cos[ k ( ϕ Φ 1 ) ] q r k = 8 Φ 2 q t max cos( k Φ 2 )tanτ / ( π 2 4 k 2 Φ 2 2 ) (19)where Bo is the width of outer ring; R is the neutral layer radius of outer ring; qrk is the coefficient in the k term of series.

In Figure 5, the radial loads on the steel balls are expressed by Fi and each loading point is marked with a number in an anti-clockwise direction. And the angular position of the ith loading point is given by: ϕ i = 2 π ( i 1 ) / N + α R , i = 1,2,... , N (20)where αR is the rotation angle of balls which determines the angular positions of loading points.

Then the corresponding deformation compatibility equation can be established at each loading point. The equation of the ith loading point is given by: w i = w cam , i w q , i P d / 2 (21)where wcam,i is the pre-deformation caused by cam assembly; wq,i is the deformation caused by the external radial load qr; wi is the deformation caused by the internal radial load Fi; Pd is the radial clearance of the FB.

3.3 The superposition algorithm of three-force system

Unknown radial loads Fi (i = 1, 2, 3, …, N, N ≧ 5) shown in Figure 5 can be decomposed into N unknown three-force systems. Each force system is shown in Figure 3, in which the radial force in the middle is Xi (i = 1, 2, 3, …, N, N ≧ 5) and its angular position is determined by equation (20). The radial deformation of different force system at each loading point can be obtained by equations (15) and (16). After all the force systems are superimposed, the total radial deformations at N loading points can be expressed by matrix: [w1w2wN1wN]=[H(0)H(ϕ0)H((N2)ϕ0)H((N1)ϕ0)H(ϕ0)H(0)H(ϕ0)H((N2)ϕ0)H((N2)ϕ0)H(ϕ0)H(0)H(ϕ0)H((N1)ϕ0)H((N2)ϕ0)H(ϕ0)H(0)]×[X1X2XN1XN](22)

where ϕ0 = 2π/N. The above matrix can be simplified as: [ w ] = [ A ] [ X ] (23)where [w] is the column matrix of radial deformation of outer ring; [A] is the deformation coefficient matrix; [X] is the column matrix of unknown forces.

The deformations at loading points in column matrix [w] can be calculated by equation (21), then we can get a set of linear equations. By solving the equations, [X] is given by: [ X ] = [ A ] 1 [ w ] (24)

The unknown radial loads in Figure 4 can be calculated by superposing unknown forces Xi at each loading point. The superposition process is written in the form of matrix, which is given by: [ F 1 F 2 F N 1 F N ] = [ 1 c 0 0 c c 1 c 0 0 0 c 1 0 c c 0 0 1 c c 0 0 c 1 ] × [ X 1 X 2 X N 1 X N ] (25) where c = 1 / ( 2cos ϕ 0 ) . The above matrix can be simplified as: [ F ] = [ C ] [ X ] (26)where [F] is the column matrix of radial loads; [C] is the transformation matrix.

According to the above linear algorithm, the internal radial load of the FB can be calculated. If there are negative results, it indicates that these steel balls are not bearing load. In order to eliminate the negative load, the force boundary condition Fi = 0 should be applied for the position where the ball is not under load. The above calculations can be programmed, and the diagram is shown in Figure 6.

And it should be noted that the superposition algorithm of three-force system is not suitable for the case that the number of balls is less than 5. When there are three balls, the radial load couldn't be determined because the radial deformations at loading points are the same by equations (15) and (16), it is impossible to establish the linear equation. And when there are two or four balls, the three-force system is not balanced, but the radial load can be calculated by two-force ring model.

thumbnail Fig. 6

The block diagram of load distribution calculations for the FB.

3.4 The verification of the flexible bearing analytic model

In order to compare the load distribution results with Xiong's results, the main parameters of bearing LY-6025 in that paper are taken to the above model [12]. The ball load distribution curves with no external loads are shown in Figure 7. The line with squares represents the ball load distribution result calculated by the three-force ring superposition method, while the line with dots represents that by the static analysis model and the line with triangles represents that by the FEM simulation model in Xiong's paper [12]. It can be seen from the figure that the result of three-force ring superposition method and is similar to that of the compared model. The maximal error between two numerical solutions is within 10%, which occurs around the major axis. And both of the number of no-load balls is 8. The reason for the differences may be the hertz contact deformation and the nonlinear part in the compared model. Therefore, the correctness and validity of the three-force ring superposition model are proved.

thumbnail Fig. 7

The ball load distribution comparison of the FB LY-6025.

3.5 The characteristics of load distribution

3.5.1 The influence of number of balls on maximum ball load

In order to obtain the characteristics of load distribution, we took maximum ball load as study object, similar results are obtained for other ball loads. The relevant calculation parameters are listed in Tables 1 and 2.

According to the algorithm in 3.3, the relationship of maximum ball load and number of balls under load is derived, as shown in Figure 8. It can be seen from figure that the maximum ball load under different load torque always decreases with the increase of number of balls.

thumbnail Fig. 8

Relationship of maximum ball load and number of balls under load.

3.5.2 The influence of load torque on maximum ball load

When αR = 0, the relationship of maximum ball load and load torque under different number of balls is derived, as shown in Figure 9. It can be seen from figure that the maximum ball load increases linearly with the increase of the load torque. And the larger the number of balls, the smaller the slope of the line. This explains that when load torque is large, the maximum ball load of the FB with more balls is far less than that of the FB with less balls.

thumbnail Fig. 9

Relationship of maximum ball load and load torque under different number of balls.

3.5.3 The relationship of ball load and angular position

When the number of balls is 23, the continuous load curves under different load torque can be obtained by changing the rotation angle of balls, as shown in Figure 10. It can be seen from figure that the amplitude of ball load is increased significantly by load torque. When T = 0, the width of no-load area at the end of minor axis is 30.39°; When T = 50 Nm, the width of no-load area increases to 44.1°. Apparently, the load torque narrows the load range of the flexible bearing, resulting in the reduction of the number of loaded balls. For the FB with 23 balls, the number of loaded balls decreases from 19 to 17 after the load torque increases from 0 to 50 Nm. Thus, high load torque on the FB should be avoided, which may cause the noise and the abnormal operation of the FB.

thumbnail Fig. 10

Relationship of ball load and angular position under load.

4 The deformation of outer ring of the flexible bearing

4.1 The expression of the radial deformation

An analytical method to account for the structural deformation of the rings was derived [21]. And it assumes that the single load applied to the thin ring is balanced by a symmetric and tangential shear stress distribution, which is not suitable for this paper. In the former two chapters, the relationship between radial deformation and load of three-force ring was derived and the unknown loads were decomposed into several three-force systems. During the above process, the sizes and positions of these three-force systems were derived. Therefore, the total radial deformation of all points on the ring could be calculated by superimposing the radial deformation caused by each three-force system. For the FB with N balls, the radial loads on the outer ring can be decomposed into N three-force systems, and Xi represents the size of the ith force system. According to the equations (15), (16) and (20), the radial deformation caused by the ith force system is given by: w i ( ϕ ) = X i H ( β i ) (27)where βi is given by:   (28)

Therefore, the total radial deformation of outer ring under N three-force systems and radial distributed load qr is given by: w ( ϕ ) = w q + i = 1 N X i H ( β i ) (29)

where wq is determined by equation (19).

4.2 The characteristics of the radial deformation of outer ring

4.2.1 The influence of number of balls on the radial deformation of outer ring

With the increase of the number of balls, the theoretical deformation will be closer to the ideal pre-deformation. In order to study the influence of number of balls on radial deformation, the difference of the expectation and the theoretical result is calculated, which is given by: Δ w = w ( ϕ ) w c a m ( ϕ ) (30)where w ⁡⁡ (ϕ) is the theoretical radial deformation, which can be calculated by equation (29); wcam(ϕ) is the ideal radial deformation, which is determined by equation (17).

When αR = 0 and T = 0, the relationship curve of the maximum difference Δwmax and the number of balls N is obtained, shown in Figure 11. It can be seen from figure that the maximum difference is negative and it appears at the end of the major axis when N ≤ 10. This will have an adverse effect on the engaging state of teeth here. And when N > 10, the maximum difference is positive and it appears at the end of the minor axis. Overall, the maximum difference decreases with the increase of the number of balls.

thumbnail Fig. 11

Relationship of the maximum difference and number of balls under no-load.

4.2.2 The influence of load torque on the radial deformation of outer ring

When N = 23 and αR = 0, the relationship curve of the difference of radial deformation Δw and angular position ϕ under T = 0 and T = 50 Nm is obtained, as shown in Figure 12. In the figure, the dash dotted line represents the theoretical deformation shape when T = 0 and the solid line represents the theoretical deformation shape when T = 50 Nm. And the expected deformation shape is indicated by the dotted line, which represents zero position. And the dots indicate the position of the balls. It can be seen from figure that the amplitude of Δw can be increased by load torque. When T = 50 Nm, the maximum difference at the end of minor axis is 2.3 times of that under no-load. And the deformation shape is actually wavelike. The theoretical curve coincides with the expected curve at the loading points, and the shapes are similar around these points. But there is an obvious difference between the theoretical and expected deformation shape around the minor axis. And the shape is convex, which results in no-load balls. If the radial deformation here is too large, it may cause that the teeth around the minor axis of the FS mesh with the teeth of steel wheel, which should be avoided in harmonic drive. And it can be found that the number of no-load balls will increase under load torque, which is consistent with the result in 3.5.3.

thumbnail Fig. 12

Relationship of the difference of radial deformation and angular position under load.

5 The bending normal stress of outer ring of flexible bearing

5.1 The bending normal stress distribution on the outer ring

The method for calculating the bending moment caused by ball load is similar to the method in 4.1. According to equations (6) and (7), the bending moment at any section of outer ring MF(ϕ) can be calculated by superimposing the bending moment caused by each three-force system, which is given by: M F ( ϕ ) = i = 1 n X i G M ( β i ) (31)where GM(ϕ) is given by: G M ( ϕ ) = { R ( h 3 h 4 cosϕ sinϕ ) / 2 0 ϕ ϕ 0 R ( h 3 h 5 cosϕ ) / 2 ϕ 0 ϕ π (32)

Substituting equation (19) into equation (1), the bending moment caused by radial distributed load is given by: M q ( ϕ ) = B o R 2 k = 2,4,6... q r k / ( k 2 1 ) cos[ k ( ϕ ϕ 1 ) ] (33)

Consequently, the total bending moment on the outer ring is given by: M ( ϕ ) = M F ( ϕ ) + M q ( ϕ ) (34)

Considering the small wall thickness of outer ring, the bending normal stress of outer ring can be calculated by bending theory of straight beam [22], whose expression is given by: σ M = M ( ϕ ) ( r R ) / I (35)where r is the radius of anywhere on cross section of outer ring; R is the neutral layer radius of outer ring; I is the inertial moment on cross section.

When N = 23 and αR = 0, the bending normal stress of outer ring under no load can be calculated by equation (35), the results are shown in Figure 13. The stress curves are smooth under no-load condition. For the bending normal stress outside the cross section of outer ring, the maximum tensile stress which appears at the end of major axis is 182.7 MPa, and the maximum compressive stress which appears at the end of minor axis is −171.6 MPa. And the distribution of the bending normal stress inside the cross section of outer ring is contrary.

When load torque T = 50 Nm, the curve of the bending normal stress distribution changes. It can be seen from figure that load torque will cause obvious extreme points on the stress curve. And the maximum bending normal stress at the end of the major axis increases apparently, but the maximum bending stress at the end of the minor axis barely change.

thumbnail Fig. 13

Bending normal stress distribution of outer ring of the FB.

5.2 The characteristics of maximum bending normal stress

5.2.1 The influence of number of balls on the maximum bending normal stress

The Figure 14 shows the relationship between the maximum bending normal stress of lateral and medial outer ring and the number of balls under no load. The lateral stress represents the maximum bending normal stress outside the cross section of outer ring, and the medial stress represents that inside the cross section of outer ring. It can be found that the lateral and medial stress have the same trend, that is the stress on major axis decreases with the increase of number of balls but the stress on minor axis stays the same. For example, looking into the curves of the lateral stress, we can find that the maximum tensile stress decreases from 249.1 to 181.8 MPa after the number of balls increases from 7 to 25, which is reduced by 27%. And the maximum compressive stress on minor axis is maintained at about 170 MPa.

thumbnail Fig. 14

The influence of number of balls on the maximum bending normal stress.

5.2.2 The influence of load torque on the maximum bending normal stress

The relationship between the maximum bending normal stress of lateral and medial outer ring and load torque is shown in Figure 15. It can be found that the bending normal stress on major axis increases with the increase of load torque but the bending normal stress on minor axis stays the same. For example, looking into the curves of the lateral stress, we can find that the maximum tensile stress increases from 182.7 to 235.5 MPa after the load torque increases from 0 to 100 Nm, which is augmented by 28.9%. And the maximum compressive stress on minor axis is maintained at about −175 MPa.

thumbnail Fig. 15

The influence of load torque on the maximum bending normal stress.

6 Conclusions

A load analysis model of three-force ring for the FB was proposed in this paper and the main conclusions were as follows:

  • Based on the theory of thin-walled ring, a mechanical model of three-force ring for the FB was established to derive the relationship between radial deformation and load of the outer ring. And then the load distribution in the FB was obtained by superposition of three-force ring. This superposition algorithm of three-force ring is suitable for the case that the number of balls is more than 4. And the calculation speed of this linear algorithm is much higher than FEMs.

  • According to the superposition algorithm of three-force ring, the influence of number of balls, load torque and angular position on the maximum ball load was studied.

The results show that the maximum ball load always increases with the increase of number of balls and it increases linearly with the increase of the load torque. For the FB with 23 balls, the width of no-load area at the end of minor axis increases from 30.39° to 44.1° after the load torque increases from 0 to 50 Nm. Consequently, the load torque narrows the load range of the FB, resulting in the reduction of the number of loaded balls. Thus, the load torque should be limited to avoid the noise and abnormal operation of the FB.

  • The theoretical deformation of outer ring of the FB is closer to the expected deformation with the increase of the number of balls. The load torque will increase the difference of the expectation and the theoretical radial deformation.

The deformation shape is actually wavelike. The theoretical curve coincides with the expected curve at the loading points, and they are similar in shape around these points. But there is an obvious difference between the theoretical and expected deformation shape around the minor axis. And the shape is convex, which results in no-load balls. And it should be noticed that a large radial deformation here should be avoided in harmonic drive.

  • The alternating bending normal stress on outer ring of the FB was studied. On the outer side, the maximum tensile stress appears on major axis and the maximum compressive stress appears on minor axis. On the inner side, the maximum tensile stress appears on minor axis and the maximum compressive stress appears on major axis. The radial distributed load produced by load torque will cause obvious extreme points on the stress curve. In addition, the bending normal stress on major axis will decreases with the increase of the number of balls and increases with the increase of the load torque. But the stress on minor axis is hardly affected by these two factors.

Nomenclature

N : Number of balls

R : Neutral radius of outer ring

M : Bending moment

T : Load torque

Xi : Middle radial force in the ith three-force system

Fi : The ith radial load on outer ring

wAB : Radial deformation of three-force ring

Δwmax : Maximum difference of radial deformation

H (ϕ): Coefficient of bending deformation

Acknowledgments

This project is supported by cooperative foundation of BEIJING CTKM HARMONIC DRIVE CO., LTD and Shanghai university (Grant No: D.71-0109-18-076).

References

  1. I.E. Lyuminarsky, S.E. Lyuminarsky, E.S. Lyu minarskaya, Increased Durability of Flexible Gear of Harmonic Drive Due to Flexible Bearing Modification, Proceedings of the 4th International Conference on Industrial Engineering, ICIE 2018, Moscow, Russian Federation, May 15–18, Springer, Cham, 2018, pp. 1287–1292 [Google Scholar]
  2. Y.W. Shen, Q.T. Ye, Theory and Design of Harmonic Gear Transmission, Machine Industry Press, Beijing, 1985, Chap. 5 [Google Scholar]
  3. W. Ostapski, Analysis of the stress state in the harmonic drive generator-flexspline system in relation to selected structural parameters and manufacturing deviations, Bull. Polish Acad. Sci. Technol. Sci. 58, 683–698 (2010) [Google Scholar]
  4. T.A. Harris, M.N. Kotadlas, Rolling Bearing Analysis, Taylor & Francis, New York, 2006 [Google Scholar]
  5. C. Wan, Analysis Method of Rolling Bearing, China Machine Press, 1987 [Google Scholar]
  6. K.X. Le, Y.W. Shen, Load distribution on the wave generator of harmonic gear transmission, Mech. Sci. Technol. Aerosp. Eng. 34, 38–45 (1990) [Google Scholar]
  7. F.C. Shao et al., Internal load analysis of flexible rolling bearings for harmonic gear transmission, Mach. Des. Manuf. 1, 34–37 (1991) [Google Scholar]
  8. K.X. Le, Y.X. Quan, G.R. Zhou et al., A new method for calculating load distribution on harmonic rolling bearing, J. Zhejiang Univ. 28, 396–404 (1994) [Google Scholar]
  9. L. Li, B.D. You, J.M. Wen et al., Load distribution calculation of flexible ball bearing with elliptical cam wave generator, Abstracts of International Conference on Mechanics, Materials and Structural Engineering, 2016, pp. 14– 15 [Google Scholar]
  10. Y. Jiang, Y.Z. Wang, Q. Tong et al., Theoretical calculation and simulation analysis of internal load of flexible bearing in harmonic drive, Mach. Des. Res. 33, 91–94 (2017) [Google Scholar]
  11. J.Y. Liu, Y.P. Chiu, Analysis of a thin elastic ring under arbitrary loading, J. Eng. Ind. 96, 870–876 (1974) [Google Scholar]
  12. Y. Xiong et al., Load analysis of flexible ball bearing in a harmonic reducer, J. Mech. Des. 142, 1–52 (2020) [Google Scholar]
  13. L. Li, Static and dynamic analysis of flexspline and flexible bearing under the cam with typical profile, M.S. thesis, Harbin Institute of Technology, Harbin, 2016 [Google Scholar]
  14. X.Q. Xiang, W.W. Zhang, X.J. Wang, G.Y. Yu, J. Dong, J. Zhang, Stress characteristics' analysis of flexible bearings at different speeds, Metrol. Meas. Tech. 44, 6–8 (2017) [Google Scholar]
  15. M.H. Ivanov, Harmonic Gear Drives, Visajas Kola Press, Moscow, 1981 [Google Scholar]
  16. I. Lopatukhin, L. Khaskel' berg, Deformation of thin-walled rings under the action of n concentrated and distributed forces, Int. J. Adv. Manuf. Technol. 49, 627–634 (2010) [Google Scholar]
  17. R.C. Hibbeler, Statics and Mechanics of Materials, China Machine Press, 2014 [Google Scholar]
  18. Y. Jiang, Y.Z. Wang, K. Zhao et al., Mechanical characteristics analysis of thin-walled flexible bearing in harmonic reducer, Bearing, 10–14 (2017) [Google Scholar]
  19. Y.N. Fengdai, Simulation analysis of load and test technology research for flexible bearing, Shanghai University, 2017 [Google Scholar]
  20. A.M. Wang, Y.Z. Wang, K. Zhao et al., Mechanical calculation on the flexible bearing of harmonic driver, IOP Conf. Ser. Mater. Sci. Eng. 490, (2019) [Google Scholar]
  21. G. Cavallaro, D. Nelias, F. Bon, Analysis of high-speed intershaft cylindrical roller bearing with flexible rings, Tribol. Lubr. Technol. 48, 154–164 (2005) [Google Scholar]
  22. J.R. Barber, Intermediate mechanics of materials, Mater. Today 10, B104 (2007) [Google Scholar]

Cite this article as: Y. Yu, E. Zhu, X. Chen, Y. Wang, Load analysis and deformation research of the flexible bearing based on a three-force ring superposition method, Mechanics & Industry 22, 34 (2021)

All Tables

Table 1

Relevant Parameters of radial distributed load.

Table 2

Main parameters of the FB.

All Figures

thumbnail Fig. 1

Schematic diagram of the working principle of the harmonic reducer.

In the text
thumbnail Fig. 2

Load simplification of outer ring of the FB.

In the text
thumbnail Fig. 3

Structure diagram of the ring under the action of three forces.

In the text
thumbnail Fig. 4

Internal force of the symmetric section of the ring.

In the text
thumbnail Fig. 5

Loading diagram of outer ring of the FB with N balls (N ≧ 5).

In the text
thumbnail Fig. 6

The block diagram of load distribution calculations for the FB.

In the text
thumbnail Fig. 7

The ball load distribution comparison of the FB LY-6025.

In the text
thumbnail Fig. 8

Relationship of maximum ball load and number of balls under load.

In the text
thumbnail Fig. 9

Relationship of maximum ball load and load torque under different number of balls.

In the text
thumbnail Fig. 10

Relationship of ball load and angular position under load.

In the text
thumbnail Fig. 11

Relationship of the maximum difference and number of balls under no-load.

In the text
thumbnail Fig. 12

Relationship of the difference of radial deformation and angular position under load.

In the text
thumbnail Fig. 13

Bending normal stress distribution of outer ring of the FB.

In the text
thumbnail Fig. 14

The influence of number of balls on the maximum bending normal stress.

In the text
thumbnail Fig. 15

The influence of load torque on the maximum bending normal stress.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.