Power Loss Analysis of an Oil-Jet Lubricated Angular Contact Ball Bearing: Theoretical and Experimental Investigations

: This study presents a theoretical and experimental analysis to quantify the power losses generated by an oil jet lubricated angular contact ball bearing. The analysis is conducted for a moderate speed range (N · dm product less than 10 6 ) and a limited applied load (<5% of the static capacity). The lubrication regime of each ball is studied through a theoretical model and varies from Iso-Viscous Rigid to Elasto-Hydrodynamic. Therefore, the hydrodynamic effects are considered in the power loss calculation. An experimental campaign is carried out and the influence of several parameters (applied load, oil injection temperature, speed, etc.) is studied. A good agreement is found between the developed model and the measurements. It is shown that the radial applied load has no influence on power losses, unlike speed and axial load. This can be explained by the load distribution and the hydrodynamic rolling contribution on the low loaded balls


Introduction
Rolling Element Bearings (REBs) are important components in mechanical transmission.Understanding the physics involved in these parts is of primary interest to improve the design and modeling of overall systems.REB power losses are still a topic under investigation, as can be seen in several recent articles [1][2][3][4][5].The different sources of dissipation in a bearing are known and presented in [6], and the main contributions are: -Sliding at the ball/ring contact surface area, which is generated by a speed difference between the two surfaces.This contribution is highly load-dependent and has been studied extensively [7][8][9][10].-Hydrodynamic rolling at the ball/ring contact inlet.It was first considered by Tevaarwerk [11], who considered this force as the result of pressure along the contact.More recent works [12][13][14] studied this force as a result of traction and sliding effects.-Oil film shearing at the cage and rings or ball interface.It is modeled as a journal bearing [15,16] and is generally negligible [15].-Drag and churning losses between the ball and the oil/air mixture.These contributions have been studied at high rotational speeds [15][16][17] and are predominant when the product N•dm is greater than 10 6 [15].However, for low to moderate speeds, this contribution is limited.
There are two main approaches to estimate REB power losses.The first one is based on a global approach.The idea is to predict the REB power losses without specific knowledge of each contribution and without precise information on the internal geometry and kinematics.In this field, two global models are generally found: Harris-Palmgren [18,19] and SKF [20] models.Previous studies have shown that experimental recalibration is necessary for accurate prediction [21][22][23].Therefore, this approach is not considered in the present study.
The second approach is a local one.The idea is to calculate power losses from each contribution.It requires precise information about the bearing's internal geometry (ball diameter, osculation).In the literature, not all the above-mentioned contributions are estimated by similar relationships (e.g., in [1,3,24]).As far as the influence of hydrodynamic rolling on deep groove ball bearings (DGBB) or angular contact ball bearings (ACBB) is concerned, most of the studies estimate this force in the elastohydrodynamic lubricated (EHL) regime for loaded balls.This consideration makes sense as ball bearing (BB) contacts are known to have high contact pressures.A few studies [25,26] also considered the isoviscous rigid (IVR) and piezo-viscous rigid regime (PVR) to estimate the hydrodynamic rolling force on loaded balls.These studies were performed at a very low N•d m product (∼ 10 4 and for limited contact loads.On the other hand, it has been shown by Darul et al. [6] that unloaded BB in DGBB must also be considered in the IVR regime to accurately predict REB power losses.This study also showed that, for moderate speeds (N•dm product less than 10 6 ) and low applied loads (Fr < 5% C 0 ), most of the power losses are due to hydrodynamic rolling [21].This study was carried out on radially loaded DGBB.The aim of this paper is to investigate if this approach is suitable for ACBB.To this end, an oil-jet lubricated ACBB is studied for moderate speeds and for limited axial and radial loads (equivalent load ≤ 5% of the static capacity).
The first part of the article deals with a static model, which is necessary to access contact information.Subsequently, the different lubrication regimes observed in the contacts are described and compared to the ones which may occur in a DGBB.As contact temperature is very influential on the lubrication regime, the mechanical model is coupled to a thermal one which aims to simulate the bearing thermal behavior.Next, the power loss calculations are presented, and the model is compared to experimental measurements for different operating conditions: speed, oil injection temperature, oil flow rate, and radial and axial load.Finally, a specific load case is analyzed to study load-independent power losses in the ACBB.

Static Model
A static model, according to Jones' analysis [7] and detailed in [19], has been implemented to evaluate power dissipation at each contact point within the REB.This section provides a concise summary of the fundamental procedures employed.

REB Geometry
The macro-geometry (outer and inner diameters, etc.) of the REB is assumed to be known.However, other parameters must be evaluated to compute the static method: (a) The free contact angle (when no load is applied) α 0 (Figure 1a) is defined as: where P d is the diametral clearance of the REB and A is the distance between raceway groove curvature centers ( r o , r i ), defined as: The free contact angle is given by manufacturers for ACBB.(b) The geometrical ratio γ: where α is the contact angle under loading, depending on α 0 .D is the ball diameter and d m is the REB mean diameter. 2

Distribution of Internal Loading
The aim is to calculate the load applied to each ball (Qnj) located at a position ψj (Figure 1b).In the present study, the speed range is moderate.Neither centrifugal force nor gyroscopic moment are considered, resulting in similar contact angles and loads on the inner and outer rings.From this hypothesis, the global REB equilibrium can be written as: and the load  is related to the REB radial displacement  and axial displacement  : where  is a load-deflection factor defined in [19].Finally, the global equilibrium is reached for an optimal value  and  .A set of simulations is performed for a DGBB and an ACBB.The bearing and the material properties are given in Tables 1 and 2. A radial load of 1 kN is applied on the DGBB.An axial load of 400 N is applied on the ACBB.The load distribution is shown in Figure 2.For the radially loaded DGBB, there is a loaded zone and an unloaded one.For the axially loaded ACBB, there is only a loaded zone.It is worth noting that the contact load applied to the loaded balls in the DGBB is much greater than the load applied to each RE in the ACBB.

Distribution of Internal Loading
The aim is to calculate the load applied to each ball (Q nj ) located at a position ψ j (Figure 1b).In the present study, the speed range is moderate.Neither centrifugal force nor gyroscopic moment are considered, resulting in similar contact angles and loads on the inner and outer rings.From this hypothesis, the global REB equilibrium can be written as: and the load Q n j is related to the REB radial displacement δ r and axial displacement δ a : where K n is a load-deflection factor defined in [19].Finally, the global equilibrium is reached for an optimal value δ r and δ a .A set of simulations is performed for a DGBB and an ACBB.The bearing and the material properties are given in Tables 1 and 2. A radial load of 1 kN is applied on the DGBB.An axial load of 400 N is applied on the ACBB.The load distribution is shown in Figure 2.For the radially loaded DGBB, there is a loaded zone and an unloaded one.For the axially loaded ACBB, there is only a loaded zone.It is worth noting that the contact load applied to the loaded balls in the DGBB is much greater than the load applied to each RE in the ACBB.Then, one can wonder if the contact load in the ACBB would be sufficient to reach the EHL regime.This question can be answered by examining the lubrication regime of each ball.

Characteristics Value
Young modulus (GPa) 210 Poisson coefficient (−) 0.3 Then, one can wonder if the contact load in the ACBB would be sufficient to reach the EHL regime.This question can be answered by examining the lubrication regime of each ball.

Lubrication Regime in Full Film
To investigate the lubrication regime on each contact, Moes' parameters [27] are used.

Moes' Parameters
Depending on the operating conditions, represented by L and N parameters, the dimensionless oil film thickness H can be estimated: with  = and U*, G* or W* as the Dowson's parameters [28]:

Lubrication Regime in Full Film
To investigate the lubrication regime on each contact, Moes' parameters [27] are used.

Moes' Parameters
Depending on the operating conditions, represented by L and N parameters, the dimensionless oil film thickness H can be estimated: with λ = R x R y and U*, G* or W* as the Dowson's parameters [28]: The rolling velocity v r is calculated as the average point of contact velocity of both the ring and the rolling element: Lubricants 2024, 12, 14 5 of 20 The ball rotation speed ω R and the cage rotational speed ω c are defined in [19].

Mapping the Lubrication Regime
It is necessary to reconstruct the Moes' map since it relies on the ellipticity of the contact.The map is created using the asymptotic solutions of each lubrication regime, as thoroughly explained in [27].The outcomes of this mapping are illustrated in Figure 3a, where the different lubrication regimes are also represented.Simulations were run for the above-mentioned DGBB and ACBB, using the loading conditions defined in Section 2, a rotational speed of 6500 rpm, and a mineral oil with characteristics given in Table 3.The dimensionless oil film thickness for each ball was plotted for different contact temperatures ranging from 30 • C to 80 • C. For the DGBB (a), it is evident that the lubrication regime for loaded RE is in the EHL zone (as an example, regimes associated with two different loaded balls are plotted in Figure 3a).For the unloaded balls, the lubrication regime is IVR.However, for the ACBB (b), all RE are loaded and the points are not located in the EHL zone but in the piezo-viscous rigid (PVR) regime.Based on this behavior, two points can be noticed here: - The ball lubrication regime may change from one ball to another, due to the load distribution.

-
The ball lubrication regime may change according to contact temperature.Another way to visualize the regime change is proposed by the authors.The idea is to plot the evolution of the lubrication regime for each ball (contact load constant) as a function of the contact temperature.This visualization, shown in Figure 3c for DGBB and Figure 3d for ACBB, makes it possible to highlight the influence of contact temperature on the lubrication regime.
Lubricants 2024, 12, x FOR PEER REVIEW 5 of 20 The rolling velocity  is calculated as the average point of contact velocity of both the ring and the rolling element: The ball rotation speed  and the cage rotational speed  are defined in [19].

Mapping the Lubrication Regime
It is necessary to reconstruct the Moes' map since it relies on the ellipticity of the contact.The map is created using the asymptotic solutions of each lubrication regime, as thoroughly explained in [27].The outcomes of this mapping are illustrated in Figure 3a, where the different lubrication regimes are also represented.Simulations were run for the above-mentioned DGBB and ACBB, using the loading conditions defined in Section 2, a rotational speed of 6500 rpm, and a mineral oil with characteristics given in Table 3.The dimensionless oil film thickness for each ball was plotted for different contact temperatures ranging from 30 °C to 80 °C.For the DGBB (a), it is evident that the lubrication regime for loaded RE is in the EHL zone (as an example, regimes associated with two different loaded balls are plotted in Figure 3a).For the unloaded balls, the lubrication regime is IVR.However, for the ACBB (b), all RE are loaded and the points are not located in the EHL zone but in the piezo-viscous rigid (PVR) regime.Based on this behavior, two points can be noticed here: - The ball lubrication regime may change from one ball to another, due to the load distribution.

-
The ball lubrication regime may change according to contact temperature.Another way to visualize the regime change is proposed by the authors.The idea is to plot the evolution of the lubrication regime for each ball (contact load constant) as a function of the contact temperature.This visualization, shown in Figure 3c for DGBB and Figure 3d for ACBB, makes it possible to highlight the influence of contact temperature on the lubrication regime.

Characteristics
A first conclusion can be drawn from this study: the lubrication regime of ACBB RE is not as binary as for radially loaded DGBB.For ACBB, RE are neither in hydrodynamics nor in EHL, but between these two regimes.This must be considered when calculating the power losses.

Power Loss Modeling
To estimate REB power losses, three contributions are considered: sliding, hydrodynamic rolling, and drag losses.Other contributions, such as elastic rolling resistance or cage interactions, are neglected.Sliding and hydrodynamic rolling occur in the contact area, with an ellipse for BB.The dimensions of this area can be calculated by estimating the load distribution.

Ball/Ring Contact Area
For a given ball n°j,  is the ellipse dimension in the transverse rolling direction ( ⃗): and b the ellipse dimension in the rolling direction ( ⃗): with ′ the equivalent Young's modulus;  the curvature sum; and  and  the ellipse dimensionless parameters, given in [19].
Moreover, the maximal Hertz pressure on the contact can be estimated by: A first conclusion can be drawn from this study: the lubrication regime of ACBB RE is not as binary as for radially loaded DGBB.For ACBB, RE are neither in hydrodynamics nor in EHL, but between these two regimes.This must be considered when calculating the power losses.

Power Loss Modeling
To estimate REB power losses, three contributions are considered: sliding, hydrodynamic rolling, and drag losses.Other contributions, such as elastic rolling resistance or cage interactions, are neglected.Sliding and hydrodynamic rolling occur in the contact area, with an ellipse for BB.The dimensions of this area can be calculated by estimating the load distribution.

Ball/Ring Contact Area
For a given ball n • j, a is the ellipse dimension in the transverse rolling direction ( → y ): and b the ellipse dimension in the rolling direction ( → x ): with E ′ the equivalent Young's modulus; ρ equiv the curvature sum; and a ′ and b ′ the ellipse dimensionless parameters, given in [19].Moreover, the maximal Hertz pressure on the contact can be estimated by:

Sliding
Assuming a constant friction coefficient µ sl , the sliding force on each contact can be estimated from [29] as: As a first approximation, µ sl is assumed to be 0.05, consistent with global models found in the literature [20].
The sliding velocity v sl can be estimated from Jones [7] and Harris analyses [19].Finally, power losses due to sliding on inner and outer rings are:

Hydrodynamic Rolling
Hydrodynamic rolling is due to an asymmetric pressure distribution.This is directly linked to the oil film thickness and so to the lubrication regime.As has been shown previously, for low load applications, the lubrication regimes covered are IVR, PVR, and EHL.Hence, the hydrodynamic rolling contribution is not the same for each regime.
-For the IVR regime, hydrodynamic rolling can be estimated using Houpert's work [30]: -For the EHL regime, the works of Tevaarwerk [11] or Houpert [31] can be used: with k = b a -For PVR regime, Biboulet's work [31] can be implemented.It enables a smooth transition from IVR to EHL regimes.This formula has also been used in [26]: with c Bib = 0.85 1 + 4.72M −0.39 k 0.484N 0.165 Finally, power losses from hydrodynamic rolling are calculated by using the rolling velocity defined in Equation ( 11):

Drag Losses
It is assumed that the rolling elements are moving through a fluid that is a mixture of air and oil.Then, the estimation of drag losses is achieved using Pouly's work [15]: with ρ mixture as the mixture density, c d as the drag coefficient, and A c as the cross-sectional area.These parameters are defined in [15].
The three contributions considered in this study have been described.These models enable the total REB power loss to be estimated, as shown in Figure 4.As it has been underlined that the ball lubrication regime may change with contact temperature, this dissipation calculation has to be associated with a thermal model.To this end, a thermal network is used.
with  as the mixture density,  as the drag coefficient, and  as the cross-sectional area.These parameters are defined in [15].
The three contributions considered in this study have been described.These models enable the total REB power loss to be estimated, as shown in Figure 4.As it has been underlined that the ball lubrication regime may change with contact temperature, this dissipation calculation has to be associated with a thermal model.To this end, a thermal network is used.

Thermal Model
The thermal network implemented is based on Pouly's work [15] for the bearing modelling and on Brossier's work [23] for the oil injection modelling.Concerning the oil flow injection rate , only a part of the injected oil passes through the bearing.The coefficient of penetration  , which is approximately 0.4 [23], indicates this.The thermal model, in Figure 5, requires the bearing inner ring (IR) and the bearing outer ring (OR) measurements.From that, power losses (rolling and sliding) are injected at the IR/RE and OR/RE contacts, and at the air/oil mixture inside the REB (drag losses).Solving the thermal equilibrium results in the temperature of the ball.Therefore, it is assumed that the temperature at the point of contact is the average temperature of the ring and the RE.

Test Rig and Test Procedure
The test rig used is presented in [21,23].Only the most relevant details are presented here.The test rig is modular.It is composed of a motor, a torquemeter, two identical oil-

Thermal Model
The thermal network implemented is based on Pouly's work [15] for the bearing modelling and on Brossier's work [23] for the oil injection modelling.Concerning the oil flow injection rate Q, only a part of the injected oil passes through the bearing.The coefficient of penetration β, which is approximately 0.4 [23], indicates this.The thermal model, in Figure 5, requires the bearing inner ring (IR) and the bearing outer ring (OR) measurements.From that, power losses (rolling and sliding) are injected at the IR/RE and OR/RE contacts, and at the air/oil mixture inside the REB (drag losses).Solving the thermal equilibrium results in the temperature of the ball.Therefore, it is assumed that the temperature at the point of contact is the average temperature of the ring and the RE.
with  as the mixture density,  as the drag coefficient, and  as the cross-sectional area.These parameters are defined in [15].
The three contributions considered in this study have been described.These models enable the total REB power loss to be estimated, as shown in Figure 4.As it has been underlined that the ball lubrication regime may change with contact temperature, this dissipation calculation has to be associated with a thermal model.To this end, a thermal network is used.

Thermal Model
The thermal network implemented is based on Pouly's work [15] for the bearing modelling and on Brossier's work [23] for the oil injection modelling.Concerning the oil flow injection rate , only a part of the injected oil passes through the bearing.The coefficient of penetration  , which is approximately 0.4 [23], indicates this.The thermal model, in Figure 5, requires the bearing inner ring (IR) and the bearing outer ring (OR) measurements.From that, power losses (rolling and sliding) are injected at the IR/RE and OR/RE contacts, and at the air/oil mixture inside the REB (drag losses).Solving the thermal equilibrium results in the temperature of the ball.Therefore, it is assumed that the temperature at the point of contact is the average temperature of the ring and the RE.

Test Rig and Test Procedure
The test rig used is presented in [21,23].Only the most relevant details are presented here.The test rig is modular.It is composed of a motor, a torquemeter, two identical oil-

Test Rig and Test Procedure
The test rig used is presented in [21,23].Only the most relevant details are presented here.The test rig is modular.It is composed of a motor, a torquemeter, two identical oil-bath lubricated "support blocks", and a "measurement block".The test rig architecture is presented in Figure 6 and a visualization is shown in Figure 7.
Lubricants 2024, 12, x FOR PEER REVIEW 9 of 20 bath lubricated "support blocks", and a "measurement block".The test rig architecture is presented in Figure 6 and a visualization is shown in Figure 7.

Limits and Measurements Performed during This Test Campaign
For this study, the measurement block contains two 7210 REBs (see Table 1).The tested operating conditions are associated with a maximum rotational speed equal to 8070 rpm, a maximum radial load equal to 1.5 kN, and a maximum axial load equal to 650 N. bath lubricated "support blocks", and a "measurement block".The test rig architecture is presented in Figure 6 and a visualization is shown in Figure 7.

Limits and Measurements Performed during This Test Campaign
For this study, the measurement block contains two 7210 REBs (see Table 1).The tested operating conditions are associated with a maximum rotational speed equal to 8070 rpm, a maximum radial load equal to 1.5 kN, and a maximum axial load equal to 650 N.

Limits and Measurements Performed during This Test Campaign
For this study, the measurement block contains two 7210 REBs (see Table 1).The tested operating conditions are associated with a maximum rotational speed equal to 8070 rpm, a maximum radial load equal to 1.5 kN, and a maximum axial load equal to 650 N. The OR and IR temperatures are measured with thermocouples (a telemetry device is used for the IR).Other temperatures can be measured, such as for ambient air, block surface, oil injection temperature, etc.The radial load is applied with a hydraulic cylinder.Axial load is applied with corrugated lock washers.Oil is injected using a peristaltic pump.The total shaft torque loss is measured with the torquemeter.The REB torque loss is obtained by subtracting the REB support block torque.A calibration phase is therefore required.

Calibration
This phase is necessary to isolate the power losses of the tested REB from the total measured torque.The objective is to estimate the power dissipation of the support block as a function of speed, loads and outer ring temperatures.As explained in [21], a calibration block is used instead of the measurement one.This block has two lip seals and two oil-bath lubricated REBs which are identical to the ones used in the support blocks (DGBB).The calibration process takes place in two steps.The first one consists of the measurement of the lip seal power losses.The second one refers to the measure of support REB resistive torque.This procedure is repeated for several speeds, loads, and temperatures.Therefore, the support blocks are characterized, and tests can be performed.

Range and Accuracy
The torquemeter used is a non-contact one.Its full scale is 10 N•m and its accuracy is about 0.02 N•m.Temperatures are measured with type T thermocouples.Their accuracy is 0.5 • C within the range of −40 • C to +125 • C. The force sensor has a nominal load of 20 kN and an accuracy of 0.4% of the nominal value.The oil flow accuracy is approximately 5%.

Oil Lubrication
A mineral oil is used to perform tests.Its physical properties are presented in Table 2.The oil flow rate is controlled with a peristaltic pump (two oil flow rates were tested: 10 L/h and 15 L/h) and the oil is heated at 50 • C or 70 • C. Oil temperature injection is measured just before injection.

Rotational Speed
The impact of the rotational speed has been tested from 3200 rpm to 8070 rpm.Two intermediate speeds (4800 rpm and 6500 rpm) have been investigated.

Load
As this study aims to focus on limited applied load, the radial load has been changed from 50 N to 1.5 kN (< 5% of the static capacity).Regarding the axial load, a preload is applied to the bearing from 225 N to 675 N.

Test Matrix
Finally, the test matrix is indicated in Table 4.

Power Loss Measurement Results
This part investigates the power loss measurement results obtained with the

Evolution with Speed
Figure 8 shows the evolution of the inner ring temperature for different speeds.All the other parameters are identical for these experiments referring to test n • 1 (Table 4).It can be observed that the higher the speed is, the faster the IR heats up.

Power Loss Measurement Results
This part investigates the power loss measurement results obtained with the 7210 ACBB.The influence of each parameter is studied.

Evolution with Speed
Figure 8 shows the evolution of the inner ring temperature for different speeds.All the other parameters are identical for these experiments referring to test n°1 (Table 4).It can be observed that the higher the speed is, the faster the IR heats up.The corresponding power losses are plotted in Figure 9a.First of all, it can be seen that power losses decrease with time.This behavior can be explained by the oil viscosity decrease as the bearing heats up.Moreover, power losses increase with speed.But, since the REB thermal behavior is not identical for the different curves (for example, at 1000 s, the IR temperature is around 90 °C at 8070 rpm and 55 °C at 3200 rpm), it seems more relevant to compare these evolutions at the same thermal equilibrium.For this reason, in Figure 9b power losses are plotted with respect to the OR-IR mean temperature.It can be concluded that, at an identical mean temperature, power losses are multiplied by 3 when the speed is twice higher.Furthermore, the OR-IR mean temperature studied here is included in the following range: (25-70 °C).The corresponding power losses are plotted in Figure 9a.First of all, it can be seen that power losses decrease with time.This behavior can be explained by the oil viscosity decrease as the bearing heats up.Moreover, power losses increase with speed.But, since the REB thermal behavior is not identical for the different curves (for example, at 1000 s, the IR temperature is around 90 • C at 8070 rpm and 55 • C at 3200 rpm), it seems more relevant to compare these evolutions at the same thermal equilibrium.For this reason, in Figure 9b power losses are plotted with respect to the OR-IR mean temperature.It can be concluded that, at an identical mean temperature, power losses are multiplied by 3 when the speed is twice higher.Furthermore, the OR-IR mean temperature studied here is included in the following range: (25-70 • C).

Evolution with Oil Injection Temperature
Figure 10 shows the evolution of the power losses for two different oil injection temperatures (25 • C and 50 • C) at 6500 rpm (referring to test n • 2).This figure highlights that power losses are not modified by the oil injection temperature.As explain by Dowson and Higginson [28], these findings confirm that the controlling viscosity is associated with the bulk temperature of the REB elements and not with the oil injection one.Therefore, changing the oil temperature can affect the thermal behavior of the bearing, but not the losses directly.

Evolution with Oil Injection Temperature
Figure 10 shows the evolution of the power losses for two different oil injection temperatures (25 °C and 50 °C) at 6500 rpm (referring to test n°2).This figure highlights that power losses are not modified by the oil injection temperature.As explain by Dowson and Higginson [28], these findings confirm that the controlling viscosity is associated with the bulk temperature of the REB elements and not with the oil injection one.Therefore, changing the oil temperature can affect the thermal behavior of the bearing, but not the losses directly.

Evolution with Oil Flow Rate
In Figure 11, the evolution of the power losses is plotted for two different oil flow rates (10 L/h and 15 L/h) at 8070 rpm (referring to test n°3).The power losses are slightly modified by the oil flow rate.This is due to the speed regime, which is moderate in this

Evolution with Oil Flow Rate
In Figure 11, the evolution of the power losses is plotted for two different oil flow rates (10 L/h and 15 L/h) at 8070 rpm (referring to test n • 3).The power losses are slightly modified by the oil flow rate.This is due to the speed regime, which is moderate in this case study.Indeed, the drag contribution, which is the only one that depends on the oil flow rate, is not predominant at moderate rotational speed.case study.Indeed, the drag contribution, which is the only one that depends on the oil flow rate, is not predominant at moderate rotational speed.

Evolution with Axial Preload
The evolution of power losses for different axial preloads is shown in Figure 12.In Figure 12a, the radial load is equal to 0 N, and three axial preloads are tested at 4800 rpm (referring to test n°4).In Figure 12b, the radial load is equal to 650 N, and two axial preloads are tested at 6500 rpm (referring to test n°5).For all cases, increasing the axial pre-

Evolution with Axial Preload
The evolution of power losses for different axial preloads is shown in Figure 12.In Figure 12a, the radial load is equal to 0 N, and three axial preloads are tested at 4800 rpm (referring to test n • 4).In Figure 12b, the radial load is equal to 650 N, and two axial preloads are tested at 6500 rpm (referring to test n • 5).For all cases, increasing the axial preload increases power losses.These observations, for different operating conditions, are consistent with global power loss models [18][19][20].

Evolution with Axial Preload
The evolution of power losses for different axial preloads is shown in Figure 12.In Figure 12a, the radial load is equal to 0 N, and three axial preloads are tested at 4800 rpm (referring to test n°4).In Figure 12b, the radial load is equal to 650 N, and two axial preloads are tested at 6500 rpm (referring to test n°5).For all cases, increasing the axial preload increases power losses.These observations, for different operating conditions, are consistent with global power loss models [18][19][20].

Evolution with Radial Load
The evolution of power losses for different radial loads is shown in Figure 13.In Figure 13a, the axial preload is equal to 450 N and three radial loads are tested at 4800 rpm (referring to test n°6).In Figure 13b, the axial preload is equal to 450 N and two radial

Evolution with Radial Load
The evolution of power losses for different radial loads is shown in Figure 13.In Figure 13a, the axial preload is equal to 450 N and three radial loads are tested at 4800 rpm (referring to test n • 6).In Figure 13b, the axial preload is equal to 450 N and two radial loads are tested at 6500 rpm (referring to test n • 7).These figures underline that power losses are almost uninfluenced by the radial load.This result is quite surprising and needs further investigation in the next part of the paper through the use of the developed model.

Comparison According to Different Parameters
The experiments are compared with the presented model.Figure 14 shows the comparison for different parameters that have been investigated through the experimental campaign.The black curves represent the power losses calculated with the developed

Comparison between Measurements and Power Loss Calculations 6.1. Comparison According to Different Parameters
The experiments are compared with the presented model.Figure 14 shows the comparison for different parameters that have been investigated through the experimental campaign.The black curves represent the power losses calculated with the developed model.The gray shaded areas represent the measurement uncertainty.Several points can be highlighted from these comparisons:

Comparison According to Different Parameters
The experiments are compared with the presented model.Figure 14 shows the comparison for different parameters that have been investigated through the experimental campaign.The black curves represent the power losses calculated with the developed model.The gray shaded areas represent the measurement uncertainty.Several points can be highlighted from these comparisons:  -In Figure 14a, the model shows a very good agreement with experiments.Speed appears to be correctly taken into account in the model.The power loss distribution at 50 °C is shown in Figure 15 for different rotational speeds.Hydrodynamic rolling is the major contributor (about 95% of the losses).Sliding is very small (<1%), as well as drag, especially at low speeds.-In Figure 14a, the model shows a very good agreement with experiments.Speed appears to be correctly taken into account in the model.The power loss distribution at 50 • C is shown in Figure 15 for different rotational speeds.Hydrodynamic rolling is the major contributor (about 95% of the losses).Sliding is very small (<1%), as well as drag, especially at low speeds.
-In Figure 14a, the model shows a very good agreement with experiments.Speed appears to be correctly taken into account in the model.The power loss distribution at 50 °C is shown in Figure 15 for different rotational speeds.Hydrodynamic rolling is the major contributor (about 95% of the losses).Sliding is very small (<1%), as well as drag, especially at low speeds.-In Figure 14b, the influence of the oil temperature injection is similar for the experiments and for the developed model.It emphasizes that the oil temperature which influences power losses is the one at the contact interface between the balls and the rings.This oil temperature is directly related to the temperature of the balls and the rings, and not to the oil injection one.

-
In Figure 14c, the calculated power losses are slightly modified by the oil flow rate.The oil flow rate only changes the drag contribution through the intermediary of the oil volume fraction.This fraction modifies the mixture properties, especially its density.Multiplying the oil flow rate by 1.5 leads to an increase of drag losses equal to 16%.However, as shown in Figure 15, the drag contribution represents only 5% of the total power losses.Therefore, the influence of the oil flow rate is limited.It should -In Figure 14b, the influence of the oil temperature injection is similar for the experiments and for the developed model.It emphasizes that the oil temperature which influences power losses is the one at the contact interface between the balls and the rings.This oil temperature is directly related to the temperature of the balls and the rings, and not to the oil injection one.-In Figure 14c, the calculated power losses are slightly modified by the oil flow rate.The oil flow rate only changes the drag contribution through the intermediary of the oil volume fraction.This fraction modifies the mixture properties, especially its density.Multiplying the oil flow rate by 1.5 leads to an increase of drag losses equal to 16%.However, as shown in Figure 15, the drag contribution represents only 5% of the total power losses.Therefore, the influence of the oil flow rate is limited.It should be noted that the oil flow is still necessary to avoid starvation effects, which could have an impact on the hydrodynamic rolling contribution.- In Figure 14d, a good agreement is found between the model and the experiments, especially at high temperatures for a stabilized thermal behavior.One can note a constant increase when the axial preload changes from 225 N to 450 N, whatever the temperature is.
If the influence of speed, preload, oil temperature, and flow rate seem well understood and predicted, the impact of radial load needs a more in-depth analysis.

Study of the Radial Load Influence
The study of radial load influence is performed on test n • 6.The rotational speed is equal to 4800 rpm, the axial preload to 450 N, and three values of the radial load are investigated: 50 N, 650 N, and 1500 N.
First, the load distribution is determined for these different load cases, as shown in Figure 16: From this distribution and the estimation of the ball temperature through the thermal network, the lubrication regime of each ball is estimated as a function of the contact temperature (see Figure 17).It can be observed that the PVR regime is present when the contact load is not very high (up to 60 N).Beyond this value, the contact load is sufficient to reach the EHL regime.When the contact load is very low, the IVR regime is found.
From the regime distribution, the power losses are calculated and compared with measurements.At 4800 rpm, Figure 18 shows that power losses calculated with the model are slightly higher than the measured ones.However, the values are within the measurement uncertainty.Moreover, theoretical calculations confirm that the radial load has almost no influence on power losses (in the tested range).In order to understand the influence of lubrication regime on hydrodynamic rolling power losses, Figure 19 shows the evolution of hydrodynamic rolling distribution at 50 • C with radial load.In this figure, 'Rolling EHL' corresponds to power losses calculated by just considering the EHL regime on all the balls, whereas 'Rolling PVR & IVR' corresponds to power losses added by considering the PVR and IVR regimes on low loaded balls.These results highlight that hydrodynamic rolling in the EHL regime decreases as the radial load decreases.But this decrease is compensated by hydrodynamic rolling in PVR & IVR regimes, which leads to an almost constant value of power losses (similar to experimental results). - In Figure 14d, a good agreement is found between the model and the experiments, especially at high temperatures for a stabilized thermal behavior.One can note a constant increase when the axial preload changes from 225 N to 450 N, whatever the temperature is.
If the influence of speed, preload, oil temperature, and flow rate seem well understood and predicted, the impact of radial load needs a more in-depth analysis.

Study of the Radial Load Influence
The study of radial load influence is performed on test n°6.The rotational speed is equal to 4800 rpm, the axial preload to 450 N, and three values of the radial load are investigated: 50 N, 650 N, and 1500 N.
First, the load distribution is determined for these different load cases, as shown in Figure 16: From this distribution and the estimation of the ball temperature through the thermal network, the lubrication regime of each ball is estimated as a function of the contact temperature (see Figure 17).It can be observed that the PVR regime is present when the contact load is not very high (up to 60 N).Beyond this value, the contact load is sufficient to reach the EHL regime.When the contact load is very low, the IVR regime is found.From the regime distribution, the power losses are calculated and compared with measurements.At 4800 rpm, Figure 18 shows that power losses calculated with the model stood and predicted, the impact of radial load needs a more in-depth analysis.

Study of the Radial Load Influence
The study of radial load influence is performed on test n°6.The rotational speed is equal to 4800 rpm, the axial preload to 450 N, and three values of the radial load are investigated: 50 N, 650 N, and 1500 N.
First, the load distribution is determined for these different load cases, as shown in Figure 16: From this distribution and the estimation of the ball temperature through the thermal network, the lubrication regime of each ball is estimated as a function of the contact temperature (see Figure 17).It can be observed that the PVR regime is present when the contact load is not very high (up to 60 N).Beyond this value, the contact load is sufficient to reach the EHL regime.When the contact load is very low, the IVR regime is found.From the regime distribution, the power losses are calculated and compared with measurements.At 4800 rpm, Figure 18 shows that power losses calculated with the model     Using the same methodology, the power dissipations were calculated for test n°7 at 6500 rpm and for two radial loads: 650 N and 1500 N. The results are shown in Figure 20.Once again, the model is perfectly in line with the experiments and confirms the low influence of radial load on power dissipation.

Conclusions
A theoretical power loss model is developed for oil-jet-lubricated ACBB.The model is based on a static analysis which gives the load distribution on each ball and accurate information on the contact conditions.Then, the lubrication regime of each rolling element is investigated using Moes' parameters.It is shown that this lubrication regime can vary from IVR to EHL depending on the loads and contact temperatures.Sliding, hydrodynamic rolling, and drag contributions are calculated at each contact of each ball, resulting in a global power loss generation value.Power losses are injected into a thermal network to account for heat dissipation inside the bearing.The temperature of the balls is calculated and used to estimate the oil viscosity in the contacts.There is therefore a strong coupling between heat generation and heat dissipation, which can affect the lubrication regime.
The model is then compared with measurements.The experimental study is conducted to investigate the influence of several parameters on power losses: rotational speed, oil injection temperature and flow rate, and axial and radial loads.For the range studied, it is shown that speed and axial preload are the two main influential parameters.Surprisingly, the radial load has almost no effect on the power losses.This behavior can

Conclusions
A theoretical power loss model is developed for oil-jet-lubricated ACBB.The model is based on a static analysis which gives the load distribution on each ball and accurate information on the contact conditions.Then, the lubrication regime of each rolling element is investigated using Moes' parameters.It is shown that this lubrication regime can vary from IVR to EHL depending on the loads and contact temperatures.Sliding, hydrodynamic rolling, and drag contributions are calculated at each contact of each ball, resulting in a global power loss generation value.Power losses are injected into a thermal network to account for heat dissipation inside the bearing.The temperature of the balls is calculated and used to estimate the oil viscosity in the contacts.There is therefore a strong coupling between heat generation and heat dissipation, which can affect the lubrication regime.
The model is then compared with measurements.The experimental study is conducted to investigate the influence of several parameters on power losses: rotational speed, oil injection temperature and flow rate, and axial and radial loads.For the range studied, it is shown that speed and axial preload are the two main influential parameters.Surprisingly, the radial load has almost no effect on the power losses.This behavior can be explained by variations associated with the load distribution and the lubrication regime.
As a very good agreement is found between the theoretical model and the measurements, it is highlighted that more than 90% of the power losses are due to the hydrodynamic rolling contribution.Furthermore, the load-independent power losses in ACBB can be explained by the change in the lubrication regime on the contacts.

Figure 1 .
Figure 1.(a) Free contact angle.(b) Ball location and applied loads.

Figure 1 .
Figure 1.(a) Free contact angle.(b) Ball location and applied loads.

Figure 4 .
Figure 4. Diagram of power loss calculation.

Figure 5 .
Figure 5. Thermal network of a bearing.

Figure 4 .
Figure 4. Diagram of power loss calculation.

Figure 4 .
Figure 4. Diagram of power loss calculation.

Figure 5 .
Figure 5. Thermal network of a bearing.

Figure 5 .
Figure 5. Thermal network of a bearing.

Figure 7 .
Figure 7. Visualization of the test rig.

Figure 7 .
Figure 7. Visualization of the test rig.

Figure 7 .
Figure 7. Visualization of the test rig.

Figure 8 .
Figure 8. Evolution of inner ring (IR) temperature for 4 different speeds.

Figure 8 .
Figure 8. Evolution of inner ring (IR) temperature for 4 different speeds.

Figure 9 . 5 . 3 . 2 .
Figure 9. (a) Evolution of power losses for 4 different speeds with respect to time.(b) Evolution of power losses for 4 different speeds with respect to outer ring (OR)-inner ring (IR) mean temperature.5.3.2.Evolution with Oil Injection Temperature Figure 10 shows the evolution of the power losses for two different oil injection temperatures (25 °C and 50 °C) at 6500 rpm (referring to test n°2).This figure highlights that power losses are not modified by the oil injection temperature.As explain by Dowson and

Figure 9 .
Figure 9. (a) Evolution of power losses for 4 different speeds with respect to time.(b) Evolution of power losses for 4 different speeds with respect to outer ring (OR)-inner ring (IR) mean temperature.

Figure 9 .
Figure 9. (a) Evolution of power losses for 4 different speeds with respect to time.(b) Evolution of power losses for 4 different speeds with respect to outer ring (OR)-inner ring (IR) mean temperature.

Figure 10 .
Figure 10.Evolution of power losses for 2 different oil injection temperatures with respect to OR-IR mean temperature.

Figure 10 .
Figure 10.Evolution of power losses for 2 different oil injection temperatures with respect to OR-IR mean temperature.

Figure 11 .
Figure 11.Evolution of power losses for 2 different oil flow rates with respect to OR-IR mean temperature.

Figure 11 .
Figure 11.Evolution of power losses for 2 different oil flow rates with respect to OR-IR mean temperature.

Figure 11 .
Figure 11.Evolution of power losses for 2 different oil flow rates with respect to OR-IR mean temperature.

Figure 12 .
Figure 12.Evolution of power losses for different axial preloads with respect to OR-IR mean temperature: (a) 4800 rpm, (b) 6500 rpm.

Figure 12 .
Figure 12.Evolution of power losses for different axial preloads with respect to OR-IR mean temperature: (a) 4800 rpm, (b) 6500 rpm.

Lubricants 2024 ,Figure 13 .
Figure 13.Evolution of power losses for different radial loads with respect to OR-IR mean temperature: (a) 4800 rpm, (b) 6500 rpm.

Figure 13 .
Figure 13.Evolution of power losses for different radial loads with respect to OR-IR mean temperature: (a) 4800 rpm, (b) 6500 rpm.

Figure 13 .
Figure 13.Evolution of power losses for different radial loads with respect to OR-IR mean temperature: (a) 4800 rpm, (b) 6500 rpm.

Figure 14 .
Figure 14.Comparison between model and experiments for the different parameters investigated.(a) Influence of speed (b) Influence of oil injection temperature.(c) Influence of oil flow rate.(d) Influence of axial preload.

Figure 14 .
Figure 14.Comparison between model and experiments for the different parameters investigated.(a) Influence of speed (b) Influence of oil injection temperature.(c) Influence of oil flow rate.(d) Influence of axial preload.

Figure 15 .
Figure 15.Power loss distribution for different speeds.

Figure 15 .
Figure 15.Power loss distribution for different speeds.

Figure 18 .
Figure 18.Comparison between model and experiments for three different radial loads, speed = 4800 rpm.

Figure 18 .
Figure 18.Comparison between model and experiments for three different radial loads, speed = 4800 rpm.

Figure 18 .
Figure 18.Comparison between model and experiments for three different radial loads, speed = 4800 rpm.

Figure 19 .
Figure 19.Hydrodynamic rolling distribution, at 50 • C, for three different radial loads, speed = 4800 rpm.Using the same methodology, the power dissipations were calculated for test n • 7 at 6500 rpm and for two radial loads: 650 N and 1500 N. The results are shown in Figure20.Once again, the model is perfectly in line with the experiments and confirms the low influence of radial load on power dissipation.Lubricants 2024, 12, x FOR PEER REVIEW 18 of 20

Figure 20 .
Figure 20.Comparison between model and experiments for two different radial loads, speed = 6500 rpm.

Figure 20 .
Figure 20.Comparison between model and experiments for two different radial loads, speed = 6500 rpm.

λ
Dimensions of the contact ellipse [m]A Distance between raceway groove curvature centers[m] C 0 Static capacity [N] d m Mean diameter [m] D Rolling element diameter [m] materials parameter = α E ′ [−] H Dimensionless film thickness [m] K Load-deflection factor [ N•mm −1.5 N, L Moes parameters, [rpm] P d Diametral clearance [m] P Power losses, [W] Q Oil flow rate [m 3 •s −1 ] Q n Load on a ball [N] r Raceway groove curvature radius [m] R x , R y Equivalent radii [m] U * Dimensionless speed = v r η E ′ R x [−] v r Rolling velocity [m•s −1 ] v sl Sliding velocity [m•s −1 ] Radii ratio = R x R y [−] ρ equiv Curvature sum [m −1 ] δ Deflectionor contact deformation [m] η Dynamic viscosity [Pa•s] µ sl Friction coefficient [−] ψ Location angle of a ball [rad] ω Rotational speed [rad•s −1 ] i Refers to inner ring o Refers to outer ring j Refers to ball number