Load Calculation of the Most Loaded Rolling Element for a Rolling Bearing with Internal Radial Clearance

This paper presents a new model for calculation of load for the most loaded rolling element in a rolling bearing with internal radial clearance. The calculation is based on a so-called load factor. By multiplying this factor by the value of the external radial load, the load transferred by the most loaded rolling element of the bearing is obtained. The values of the load factor are shown in the tables and diagrams, which makes the model very suitable for practical use. The load factors are shown for a ball bearing as well as for a roller bearing. The model considers two support positions of the inner ring on an even and odd number of rolling elements. The new model was compared with the most commonly used models up to now. The results showed greater accuracy of the studied model.


Introduction
One of the most important tasks of rolling bearings is to transfer load from the rotary elements to the supports. This takes place through a certain number of rolling elements, from the rotating ring to the stationary one. The load that rolling elements transfer, i.e., the internal load distribution, is an essential feature of rolling bearing, in particular for the determination of static and dynamic load rating, service life, the characteristics of stiffness, the accuracy of rotation, vibration characteristics, the losses due to friction, heat calculation, etc.
It is known that the calculation of the load distribution on the rolling elements is very complicated. It is nonlinear and a statically indeterminate problem, which requires that in addition to the static equilibrium equation, additional conditions have to be introduced, based on the connection between the contact stress and deformation, according to Hertz's contact theory [1]. Furthermore, in the case of rolling bearings, the internal load distribution is a time-varying value, which is dependent on the position of the rolling elements on the raceways. Time variability, static indeterminacy, and nonlinearity have meant that the calculation of the internal load distribution can only be solved numerically, through many iterations, using Newton's iteration or a similar method. However, because of its significant importance, especially for engineering practice, a large number of authors are dealing with methods to simplify the calculation of the internal load distribution of rolling bearings, and are trying to develop as simple but sufficiently accurate methods as possible.
The first papers that dealt with these problems are related to the study by Stribeck [2]. Stribeck used the Hertzian contact theory and derived an equation for calculation of the load of the most loaded rolling element for the radial ball bearing with zero clearance. The resulting solution for rolling bearing (ball or roller) with nominal clearances was later expanded by Palmgren [3]. Stribeck's expressions are still in use for the calculation of the static load rating of rolling bearing [4,5]. However, the Stribeck's expressions have one serious drawback. They do not consider the size of the internal radial clearance, but only its presence, which is described through the so-called Stribeck's number. That number is 4.37 for the ball bearing without radial clearance, and 4.08 for a roller bearing. For bearings with clearance greater than 0, Stribeck's number is 5. Later in the calculation process, Sjövall introduced the so-called Sjövall's integral [6].
The first comprehensive mathematical model for determining the load distribution for rolling bearing was developed by Jones [7]. He observed an ideal rolling bearing with an internal radial clearance and proposed a mathematical model to determine the load distribution as a nonlinear equation system. The model of Jones was later developed by Kovalev and Narodecki [8], as well as Harris [1]. Harris [1], using Sjövall's integrals, describes the size of the contact zone through a load distribution factor (ε) and solves the load distribution problem by the method of trial and error. Based on the assumed distribution factor, the deflection of the bearing is determined. In the next step, equilibrium equations are solved, and the process is repeated until the given accuracy is satisfied. Harris's model was later used for different purposes and has been further developed by a large number of authors, such as Grünberg [9,10] Houpert [11], Ricci [12] Oswald [13], Rusu and Tudose [14,15], Xiao et al. [16], Fang et al. [17], etc. However, the mathematical calculation using the Harris model implies many iterations of computers, which often causes complications in concrete applications in practice.
The great practical significance comes from a paper by a research group gathered around Professor Mitrovic [18][19][20][21]. Mitrovic and associates explored the impact of structural and tribological parameters on the load distribution and the working capacity of rolling bearing. In addition, Mitrovic explored the impact of load distribution on the performance of the rolling bearing radial rigidity, the static load ratings, and the service life [22]. Mitrovic used the load distribution factor as an indicator of the load distribution nonuniformity. Mitrovic's research was later expanded by Zeljkovic and Zivkovic [23]. The analysis presented in the resulting papers [18][19][20][21] shows that among the most important factors determining the degree of load distribution non-uniformity are a radial clearance and a bearing deflection. A higher level of a bearing deflection will provide a load transfer across a larger number of rolling elements and more even load distribution. However, a higher level of the bearing deflection requires higher values of the external radial loads. In further papers [24][25][26][27][28][29], these phenomena were studied, and it was concluded that for a given bearing, the internal load distribution depends mostly on the size of the external radial load and the internal radial clearance. In a paper by Tomović [27] are shown diagrams from which it is possible to directly determine the number of active rolling elements for a given bearing type, depending on the load applied and the internal radial clearance.
As a continuation of research that has been shown previously [24][25][26][27][28][29], in this paper a new mathematical model for calculating the load of the most loaded rolling element in a bearing with an internal radial clearance is derived and presented. The model considers a radially loaded rolling bearing. Inspired by previous papers [18][19][20][21][22][23], the proposed model for the calculation uses a load factor. Multiplying this factor with the value of the external radial load, the load carried by the most loaded rolling element of a bearing was obtained. For limit values of the bearing deflection, the values of the load distribution factor were derived and presented through the appropriate tables and diagrams. Factors of load distribution were derived for the rolling ball bearing and rolling bearing with cylindrical rollers. The presented mathematical model was compared with some of the so far most widely used models. The results of this analysis are presented in this paper.

Starting Assumptions
The starting assumptions were as follows: -A single-row, radially loaded rolling bearing, ball, or roller is being analyzed; - The contact surfaces of the bearing are ideally made; - The angle between the rolling elements is constant, and it can be calculated using the equation Appl. Sci. 2020, 10, 6934 The cage excludes the interplay between the rolling elements; - The direction of the external radial load passes through the center of the bearing rings; -All elements of the bearing are rigid, and there is no flexion. There are only contact deformations and corresponding to Hertz's theory of elasticity [1]; - The inner ring of the bearing is movable, while the outer is fixed; - The bearing works in the isothermal conditions.

Boundary Positions of Rolling Bearing
One of the most common simplifications used in the mathematical analysis of rolling bearing is the observation of the bearing in just one position, when the center of one of the rolling elements coincides with the direction of action of the external radial load (Figure 1a). Although from the point of load distribution, this is the most unfavorable position of support, analyzing only one position is insufficient and does not provide complete information about the behavior of the bearing in operation. Namely, during working, the set of rolling elements is constantly oscillating between two boundary positions, as shown in Figure 1. This phenomenon has been detailed in previous papers [25,28,30]. It is clear that, due to the different position of the sets of the rolling elements, the load distribution in the two boundary positions will be different. Therefore, the value of the contact forces and contact deformations between elements of the bearing will be different in these boundary positions. Differences in the level of contact force and contact deformation will inevitably cause the differences in other running characteristics of the bearing, such as the level of deflection, stiffness, vibration characteristics, rotation accuracy, and so on.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 3 of 18 -The cage excludes the interplay between the rolling elements; - The direction of the external radial load passes through the center of the bearing rings; -All elements of the bearing are rigid, and there is no flexion. There are only contact deformations and corresponding to Hertz's theory of elasticity [1]; - The inner ring of the bearing is movable, while the outer is fixed; - The bearing works in the isothermal conditions.

Boundary Positions of Rolling Bearing
One of the most common simplifications used in the mathematical analysis of rolling bearing is the observation of the bearing in just one position, when the center of one of the rolling elements coincides with the direction of action of the external radial load (Figure 1a). Although from the point of load distribution, this is the most unfavorable position of support, analyzing only one position is insufficient and does not provide complete information about the behavior of the bearing in operation. Namely, during working, the set of rolling elements is constantly oscillating between two boundary positions, as shown in Figure 1. This phenomenon has been detailed in previous papers [25,28,30]. It is clear that, due to the different position of the sets of the rolling elements, the load distribution in the two boundary positions will be different. Therefore, the value of the contact forces and contact deformations between elements of the bearing will be different in these boundary positions. Differences in the level of contact force and contact deformation will inevitably cause the differences in other running characteristics of the bearing, such as the level of deflection, stiffness, vibration characteristics, rotation accuracy, and so on.
If we want to describe as accurately as possible the behavior of the bearing in operation, then we need to introduce this second boundary support position in the analysis. In previous studies [24][25][26][27][28], the following terms were adopted for these positions of support: The first boundary position of the support occurs when the center of one of the rolling elements coincides with the direction of external radial load, and the other when the rollers are symmetrically arranged relative to the direction of action of external radial load. In the following section, the value If we want to describe as accurately as possible the behavior of the bearing in operation, then we need to introduce this second boundary support position in the analysis.
In previous studies [24][25][26][27][28], the following terms were adopted for these positions of support: The first boundary position of the support occurs when the center of one of the rolling elements coincides with the direction of external radial load, and the other when the rollers are symmetrically arranged relative to the direction of action of external radial load. In the following section, the value Appl. Sci. 2020, 10, 6934 4 of 18 relating to the boundary position of support on an odd number of rolling elements will have index o, and the value relating to the boundary position of support on an even number rolling elements will have index e.

Active Rolling Elements
Active rolling elements are those elements involved in transferring the external radial load. The area bounded by the active rolling elements is called the load zone [1].
In a radially loaded rolling bearing, the load can be transferred only by rolling elements that are under the meridian plane, and only those elements that have contact with both bearing rings. This is often not the case in bearings with radial internal clearance. In bearings with radial internal clearance, due to the negative impact of clearance, it is necessary to achieve extra relative movements between the bearing rings, so that the individual rolling elements come into contact with both rings and begin to transfer external radial load. This is only possible if there is an adequate level of contact deformations between the elements of the bearing, i.e., if the bearing is loaded with the adequate intensity of external radial force.
The number of active rolling elements for a given type of bearing depends solely on the internal radial clearance and external radial loads. In [24,25], this phenomenon is explained in detail. In [26,27], the procedure for quickly and easily determining the number of active elements of the bearing is presented. The size of the necessary deflection and the necessary external radial load for entry of individual rolling elements in contact with the bearing rings can be determined on the basis of the equations shown in [25,26], respectively.
Since these values are very significant for the calculation of the internal load distribution, below the basic equations for their determination will be briefly presented.

The Maximum Number of Rolling Elements in the Load Zone of the Bearing
The maximum number of rolling elements that can be found below the meridian plain exclusively depends on the total number of rolling elements in the bearing, and according to [25], can be calculated from the following equations:

•
For the BSO position: • For the BSE position:

Boundary Bearing Deflection
Tomović [25] presented forms for the calculation of the boundary values of bearing deflection, which is required to achieve load transfer over q rolling elements, as where q is the number of active rolling elements of the bearing, e is the radial clearance, and t q is the coefficient of boundary bearing deflection, obtained by the following equations [25]: • For the BSO position: Appl. Sci. 2020, 10, 6934

of 18
• For the BSE position: The coefficient t q can be found in tables from the previous study for reference [25].

Boundary Load of Bearing
The boundary external load of bearing that is required for the q of the active rolling elements can be calculated as [26] where s q is the coefficient of the boundary bearing load, K is the effective stiffness coefficient, and n is the exponent (n = 3/2 for a ball bearing and n = 10/9 for a bearing with the rollers). The coefficient s q for a BSO or BSE position can be calculated with the following equations: The coefficient s q can be read from the tables that are shown in [26]. Boundary loads that are described by Equation (7) correspond to the bearing deflection from the Equation (4).

Number of Active Rolling Elements
One of the biggest problems in the calculation of load distribution in the rolling bearing is to define the number of active rolling elements (q). This problem was solved using the trial and error method [1], and more iterations were required for a satisfactory solution. Using Equation (7) and the table given in a paper by Tomović [26], this number can be very easily determined. By simply comparing the applied external load to the calculated boundary values, it is possible to conclude how many rolling elements are engaged in load transfer.

Load Factor of the Most Loaded Rolling Element of Bearing
The most loaded rolling element of bearing is the one element that stands closest to the external radial force direction. In boundary positions of support are as shown in Figure 2. These are the rolling elements which carry the label zero. From the engineering point of view, the boundary support position that is shown in Figure 2a is particularly significant. In this position, the center of the most loaded rolling element lies in the direction of the external radial force. This is the most disadvantaged case for support, because the rolling elements are most loaded at the time when they find themselves on the external radial load direction, i.e., transfer the largest proportion of the external load. The maximum contact stress and deformation occur on the place of the most loaded rolling element at this position [29]. Therefore, the calculation of the largest number of bearing running performance is calculated concerning this position, such as static and dynamic load capacity, service life, etc.
In engineering practice and scientific research, only this support position is generally considered in the analyses. However, in many cases, this approach is inadequate, because it does not give enough information about the behavior of the bearing in operation. Namely, a set of rolling elements constantly changes its position during operation of the bearing, which causes a temporal change of the internal load distribution, as well as the temporal change of the intensity of contact stress and deformation and the stiffness of the bearing [31,32]. Therefore, the research presented in [24][25][26][27][28] proposed the introduction of another one boundary position, shown in Figure 2b. The analysis of load distribution, which is shown in this paper, will be made regarding these two boundary positions.
the internal load distribution, as well as the temporal change of the intensity of contact stress and deformation and the stiffness of the bearing [31,32]. Therefore, the research presented in [24][25][26][27][28] proposed the introduction of another one boundary position, shown in Figure 2b. The analysis of load distribution, which is shown in this paper, will be made regarding these two boundary positions. In both of the positions shown above, the rolling elements are arranged symmetrically relative to the vertical axis, representing the direction of action of the external radial force. On this basis, it can be assumed that the load distribution is symmetrical to the external load direction and the vertical bearing axis. This creates the possibility of adopting a simple labeling system, as shown in Figure 2. The most loaded rolling element is marked with 0. The other rolling elements have been given labels in ascending order relative to their position relative to the vertical bearing axis.

BSO Position of Support
The equation of equilibrium for the BSO support position can be written as (Figure 2a) where Q is the external radial load, F0 is the load of the most loaded rolling element, Fi is the load of the i-th rolling element, γ is the angular spacing between rolling bearing elements, and q is the number of active rolling elements. The load transferred by the i-th rolling elements is equal to where δi is the contact deformation at the position of the i-th rolling element. Relative to Equation (11) the load of the 0-th rolling element is equal to: By dividing Equation (11) by Equation (12), the following can be obtained: According to [26], the contact deformation δi is obtainable as In both of the positions shown above, the rolling elements are arranged symmetrically relative to the vertical axis, representing the direction of action of the external radial force. On this basis, it can be assumed that the load distribution is symmetrical to the external load direction and the vertical bearing axis. This creates the possibility of adopting a simple labeling system, as shown in Figure 2. The most loaded rolling element is marked with 0. The other rolling elements have been given labels in ascending order relative to their position relative to the vertical bearing axis.

BSO Position of Support
The equation of equilibrium for the BSO support position can be written as (Figure 2a) where Q is the external radial load, F 0 is the load of the most loaded rolling element, F i is the load of the i-th rolling element, γ is the angular spacing between rolling bearing elements, and q is the number of active rolling elements. The load transferred by the i-th rolling elements is equal to where δ i is the contact deformation at the position of the i-th rolling element. Relative to Equation (11) the load of the 0-th rolling element is equal to: By dividing Equation (11) by Equation (12), the following can be obtained: According to [26], the contact deformation δ i is obtainable as where a 2·i+1 is the boundary-bearing deflection until entering of i-th rolling element in contact with bearing rings. Furthermore, if Equation (14) is included in Equation (13), the relationship between the load of the most loaded and the arbitrary rolling element is obtained: Since this is for the BSO support position, the contact deformation of the most loaded rolling element is equal to total deflection of the bearing, i.e., δ 0 = w o ; thus, Equation (15) can be written as If the above equation is inserted in Equation (10), after arranging, the relationship between the external load and the load of the most loaded rolling element was obtained, is If it is further taken that The force which is transferred by the most loaded rolling element is obtained, as: In the equation above k max,o represents a load factor of the most loaded rolling element. This factor is, in fact, a ratio of the size of the external load and the load transferred by the most loaded rolling element of bearing.
In the boundary case, when switching from one support system to the other [25], the bearing rings' relative movement is equal to the boundary deflection, which is required for supporting the inner ring on the q of the rolling elements, i.e., w o = a q . Then the expression for the load factor of the most loaded rolling element is obtained as Further taking Equation (4) into account, after arranging, the above equation is obtained in its final form: Appl. Sci. 2020, 10, 6934

of 18
Respectively, if the expressions for t q from Equation (5) are taken into account, then, after arranging, we can obtain

BSE Position of Support
The labeling system for the BSE position of support is shown on Figure 2b. The condition of force equilibrium for this boundary position can be expressed by the following equation: Again, the most loaded rolling elements carry the 0 marks. Their contact deformations (δ 0 ), according to [19], can be calculated as follows: According to [26], contact deformation δ i can be obtained as If Equation (25) is now inserted into (13), the relation between the loads of the most loaded rolling element and on the arbitrary rolling element is obtained, after arranging, as Respectively, based on Equation (25), if w e is included instead of δ 0 , the following is obtained: By inserting the above equation into the equation of equilibrium (Equation (23)), the relationship between the most loaded rolling element and the external radial load is obtained, after arranging: where Factor k emax is the load factor of the most loaded rolling element. In the boundary moment, when switching from one to the other the support system [25], If this is included in Equation (29), it is obtained that the k emax factor is equal to If Equations (4) and (6) are further taken into account, the final form of the equation for the BSE position of support is obtained, as The above equations refer to the boundary moment when the i-th rolling element engages with the bearing rings.

The Load Factor of the Most Loaded Rolling Element: Tables and Diagrams
Factors k max and k maxq show which part of the external load is transferred through the most loaded rolling element. This load is relevant when determining the static load rating of the bearing.
Factor k maxq represents the load factor of the most loaded rolling element at the boundary moment when switching from one support system to the other [25]. At this moment, the inner ring of the bearing begins to support on q the rolling elements of bearing. This factor, in fact, represents a relationship between the boundary external radial load (Q q ) and the load transferred by the most loaded rolling element. Thus, the load that is transferred by the most loaded rolling element, in the case when the inner ring begins to support the q of rolling elements, is equal to The value Q q presents the minimal or the boundary external radial load, which is necessary for the support of the inner ring on the q rolling elements, and can be calculated using Equations (7)- (9).
The values of factor k maxq for the radial ball bearing are given in Table 1, and in Table 2 for the roller bearing. The mentioned factors are obtained using Equations (22) and (33) for a bearing with a total number of rolling elements (z), which varies in limits from 5 to 20. In the tables, given vertically, are the maximum possible number of active rolling elements. This number is shown separately for a case of BSO (z so ) and BSE (z se ) position of support, and was obtained by Equations (2) and (3). In addition, the dependence of the coefficient k maxq on the total number of rolling elements in a bearing is shown in Figures 3 and 4. The different curves in the diagrams are given in relation to the different number of active rolling elements. The diagrams shown correspond to the boundary values of the bearing deflection, at which the q-th rolling element engages into contact with both rings, i.e., when the end two rolling elements from the zone of load begin to engage in the transfer of the external load. At that moment, these two elements have a 0-th contact deformation, i.e., the 0-th load.

Discussion of Equations, Tables, and Diagrams
From Equations (22) and (33), it can be seen that factor kmaxq solely depends on the total number and number of active rolling elements. At first glance, one gets the impression that factor kmaxq is not dependent on the internal radial clearance. However, for a given external load, the number of active rolling elements is directly dependent on the internal radial clearance [27], and through this parameter, the internal radial clearance directly influences the size of the factor kmaxq. Figures 3 and 4 clearly show that with the increase of the total number of rolling elements, the value of kmaxq factor declines. However, this decline is inconsiderable, and it is getting smaller as the total number of rolling elements in the bearing increases. This means that the total number of rolling

Discussion of Equations, Tables, and Diagrams
From Equations (22) and (33), it can be seen that factor kmaxq solely depends on the total number and number of active rolling elements. At first glance, one gets the impression that factor kmaxq is not dependent on the internal radial clearance. However, for a given external load, the number of active rolling elements is directly dependent on the internal radial clearance [27], and through this parameter, the internal radial clearance directly influences the size of the factor kmaxq. Figures 3 and 4 clearly show that with the increase of the total number of rolling elements, the value of kmaxq factor declines. However, this decline is inconsiderable, and it is getting smaller as the total number of rolling elements in the bearing increases. This means that the total number of rolling

Discussion of Equations, Tables, and Diagrams
From Equations (22) and (33), it can be seen that factor k maxq solely depends on the total number and number of active rolling elements. At first glance, one gets the impression that factor k maxq is not dependent on the internal radial clearance. However, for a given external load, the number of active rolling elements is directly dependent on the internal radial clearance [27], and through this parameter, the internal radial clearance directly influences the size of the factor k maxq . Figures 3 and 4 clearly show that with the increase of the total number of rolling elements, the value of k maxq factor declines. However, this decline is inconsiderable, and it is getting smaller as the total number of rolling elements in the bearing increases. This means that the total number of rolling elements does not have a significant impact on the value of factor k maxq . With regard to the value of factor k maxq , the biggest impact is the number of active rolling elements. Besides that, the image in Figures 3 and 4 show that with the increase of the total number of rolling elements in the bearing, the value of the load factor asymptotically leans toward a constant value. These statements offer the possibility of approximation and wide-open opportunities for further research and the adoption of general values for the k maxq factor.
From Figures 3 and 4, it can also be noted that there is a great similarity between the diagrams that show the value of factors k maxq for the ball and for the bearing with rollers. However, the values of the factor k maxq at the bearings with rollers are slightly lower when supporting on q = 5-10 rolling elements. This difference in the size of k maxq factor occurs due to the different values of the exponent n, which in the ball bearing is n = 3/2, and in the roller bearing is n = 10/9.

The Validity Analysis of the Proposed Model
In order to verify the proposed model, a calculation of the load that transfers the most loaded roller element for ball bearing 6008 was made. The calculation was carried out in relation to the following equation: Here is k 0 , the load factor of the 0-th rolling element. The values of the coefficient k 0 were adopted as follows: • In cases where the assumed external load was less than the boundary load, i.e., when Q ≤ Q 3 , the load transfer inside the bearing occurs via one or two rolling elements, i.e., as per support system 1 [25,30]. In these cases, for the load calculation is used the coefficient k 0 = k max3 = 1. By increasing the load above Q 3 , new rolling elements engage in contact with the bearing rings, and the part of the load that transfers the most loaded rolling element decreases.

•
For the loads which are in the interval Q q < Q ≤ Q q+1 , the value of the load factor is obtained by interpolation, using the next equation: The load factor values selected in this way give the most accurate calculation results.

•
For the loads that are above the boundary load, corresponding to the maximum possible number of the active rolling elements, the value of the load factor k 0 = k 0zs is taken, and respectively for the bearing 6008 k 0 = k 06 .
The basic parameters of the ball bearing 6008, which are necessary for the calculation, are shown in Table 3. The coefficients t q and s q are taken from tables shown in [25,26]. The values of the load factor of the most loaded rolling element (k 0q = k maxq ) are taken from Table 1. The total number of rolling elements: z = 12 Maximum possible number of active rolling elements: z s = 6 The effective coefficient of the bearing: K = 2.88 × 10 5 N/mm 3/2 In Table 4 are shown the boundary values of external radial load Q q , corresponding to values of the internal radial clearance, which are considered in this analysis and which range from 0 to 50 µm. In this way, the entire range of recommended internal clearance values for bearing 6008 is covered, for all classes (C2 to C5) [33]. The values of the boundary loads Q q were obtained using Equation (7). The calculation results are shown in Table S1 (Supplementary Materials). Results are shown in relation to the proposed model (index F 0 ), Stribeck's method (S), and by solving the system of equilibrium equations inside the bearing (ER: exact result). The equilibrium equations were derived based on the procedure shown in [25]. The solutions obtained in this way were considered accurate. In relation to these results, the error of the proposed method and the Stribeck method is obtained in the following procedure. The calculation errors, expressed as a percentage, are shown below the calculation results.
The load varied in the range from 0 to 10,000 N, which is close to the static load rating of the ball bearing 6008, which amounts C 0 = 11,600 N. For the boundary loads shown in Table 4, the calculation was also made. These results are also shown in Table S1.
For efficient analysis, the results shown in Table S1 are shown in Figures 5-12. The left side of the figure shows loads which transfer the 0-th roller bearing element. In order to compare the calculation results, the calculation error according to the proposed method and the Stribeck method is shown on the right side.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 13 of 18 The calculation results are shown in Table S1 (Supplementary Materials). Results are shown in relation to the proposed model (index F0), Stribeck's method (S), and by solving the system of equilibrium equations inside the bearing (ER: exact result). The equilibrium equations were derived based on the procedure shown in [25]. The solutions obtained in this way were considered accurate. In relation to these results, the error of the proposed method and the Stribeck method is obtained in the following procedure. The calculation errors, expressed as a percentage, are shown below the calculation results.
The load varied in the range from 0 to 10,000 N, which is close to the static load rating of the ball bearing 6008, which amounts C0 = 11,600 N. For the boundary loads shown in Table 4, the calculation was also made. These results are also shown in Table S1.
For efficient analysis, the results shown in Table S1 are shown in Figures 5-12. The left side of the figure shows loads which transfer the 0-th roller bearing element. In order to compare the calculation results, the calculation error according to the proposed method and the Stribeck method is shown on the right side.

Conclusions
This paper presents a new mathematical model for calculating the load transferred by the most loaded rolling element of the bearing. The model is based on so-called load distribution factors. The load factors show which part of the external radial load transfers the most loaded rolling element of the bearing.
The model is characterized by great simplicity, which makes it very suitable for practical use. By simply multiplying the load factor by the size of the external radial load, the size of the load transferred over the most loaded rolling element is obtained.
The load factors are shown in the form of tables, in particular for ball bearings and in particular, for a roller bearing, and correspond to all types of a radial loaded ball and roller bearings.
The factors in the tables can be easily read based on only two parameters: the total number of rolling elements in the bearing and the number of active rolling elements.
The total number of rolling elements in the bearing does not have a significant impact on the size of the load factors. With regard to the size of the load factor, the most significant influence is of the number of active rolling elements. This fact creates a presumption for the adoption of general load factor values for different types of bearings. General factors for all types of bearings would only be dependent on the number of active rolling elements in the bearing.
Because of its simplicity, the calculation proposed in this paper can be compared to a calculation

Conclusions
This paper presents a new mathematical model for calculating the load transferred by the most loaded rolling element of the bearing. The model is based on so-called load distribution factors. The load factors show which part of the external radial load transfers the most loaded rolling element of the bearing.
The model is characterized by great simplicity, which makes it very suitable for practical use. By simply multiplying the load factor by the size of the external radial load, the size of the load transferred over the most loaded rolling element is obtained.
The load factors are shown in the form of tables, in particular for ball bearings and in particular, for a roller bearing, and correspond to all types of a radial loaded ball and roller bearings.
The factors in the tables can be easily read based on only two parameters: the total number of rolling elements in the bearing and the number of active rolling elements.
The total number of rolling elements in the bearing does not have a significant impact on the size of the load factors. With regard to the size of the load factor, the most significant influence is of the number of active rolling elements. This fact creates a presumption for the adoption of general load factor values for different types of bearings. General factors for all types of bearings would only be dependent on the number of active rolling elements in the bearing.
Because of its simplicity, the calculation proposed in this paper can be compared to a calculation

Conclusions
This paper presents a new mathematical model for calculating the load transferred by the most loaded rolling element of the bearing. The model is based on so-called load distribution factors. The load factors show which part of the external radial load transfers the most loaded rolling element of the bearing.
The model is characterized by great simplicity, which makes it very suitable for practical use. By simply multiplying the load factor by the size of the external radial load, the size of the load transferred over the most loaded rolling element is obtained.
The load factors are shown in the form of tables, in particular for ball bearings and in particular, for a roller bearing, and correspond to all types of a radial loaded ball and roller bearings.
The factors in the tables can be easily read based on only two parameters: the total number of rolling elements in the bearing and the number of active rolling elements.
The total number of rolling elements in the bearing does not have a significant impact on the size of the load factors. With regard to the size of the load factor, the most significant influence is of the number of active rolling elements. This fact creates a presumption for the adoption of general load factor values for different types of bearings. General factors for all types of bearings would only be dependent on the number of active rolling elements in the bearing.
Because of its simplicity, the calculation proposed in this paper can be compared to a calculation based on Stribeck's numbers. Both calculations are based on simple mathematical operations, do not require the usage of computers, and they are very useful for the fast calculations and fieldwork.
The analysis showed that for bearings with a radial clearance greater than zero, the results obtained by the proposed method have significantly higher precision. For bearings with a radial clearance equal to zero, the Stribeck method has greater precision.
Stribeck's number for ball bearings with zero clearance is S = 4.37, and for bearings with internal clearance other than zero S = 5. These values of Stribeck's numbers have been in use for many years, and represent the generally adopted value for calculating the load of the most loaded ball. For bearings with zero radial clearance, the use of a Stribeck number is fully justified. However, the analysis made in this paper indicates the need for the correction of the Stribeck's number for bearings with a radial clearance greater than zero. Specifically, Stribeck's calculation only differentiates whether or not the clearance is present in the bearing, regardless of its size. The Stribeck's number for all bearings with a clearance is the same, regardless of the size of the clearance.
The model proposed in this paper considers the impact of the clearance size on the load distribution within the bearing, through the number of active rolling elements of the bearing. In this regard, it is suggested that a correction of the Stribeck's number be made concerning the number of active rolling elements in the bearing.
Funding: This research received no external funding.

Conflicts of Interest:
The authors declare no conflict of interest.

Nomenclature γ
angle between rolling elements z total number of rolling elements BSO boundary position of support on an odd number of rolling elements BSE boundary position of support on an even number of rolling elements z s the maximum number of rolling elements that can be found in the loaded zone q number of active rolling elements e internal radial clearance a q deflection of bearing t q coefficient of boundary deflection of bearing Q bearing radial load Q q boundary load of bearing K Hertz's stiffness coefficient s q coefficient of boundary radial load F rolling element load δ contact deformation n exponent w bearing deflection k load factor of the rolling element k q load factor of the rolling element in the boundary moment k max load factor of the most loaded rolling element F max load of the most loaded rolling element