Next Article in Journal
A Study on the Influence of Nonlinear Vibration on Fretting Damage of Involute Spline Pairs in Aero-Engines
Next Article in Special Issue
Coarse-Grained Molecular Dynamics Simulations of Organic Friction Modifier Adsorption on Rough Surfaces under Shear
Previous Article in Journal
Numerical Simulation and Wear Resistance Property of Ni-Based Alloy Coating on the Surface of Ti-6Al-4V Substrate
Previous Article in Special Issue
A Methodological Approach to Assessing the Tribological Properties of Lubricants Using a Four-Ball Tribometer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Predicting Wear under Boundary Lubrication: A Decisive Statistical Study

1
ASC Goerlach, Robert-Bosch-Straße 60/1, 72810 Gomaringen, Germany
2
Operational Research Group, School of Mathematical Sciences, University of Southampton, Southampton SO17 1BJ, UK
3
Business School, Faculty of Social Sciences, University of Southampton, Southampton SO17 1BJ, UK
*
Author to whom correspondence should be addressed.
Lubricants 2023, 11(12), 514; https://doi.org/10.3390/lubricants11120514
Submission received: 15 September 2023 / Revised: 20 October 2023 / Accepted: 22 November 2023 / Published: 4 December 2023
(This article belongs to the Special Issue Advances in Boundary Lubrication)

Abstract

:
The forthcoming revolution in mobility and the use of lubricants to ensure ecological friendliness intensifies the pressure on tribology for predictors in new life cycles, mainly addressing wear. The current paper aims to obtain such predictors by studying how the wear processes that occur in a standard FE8 bearing test rig under thin film lubrication are conducted by the properties of the lubricant rather than simple viscosity parameters. Assuming that the activity of a lubricant with respect to the temperature, surface, and chemicals is a matter of its chemical potential, the results show that the nature of the base oil is a key parameter, apart from the chemical structure of the additives. Moreover, it becomes clear that chemical predictors are changing by altering the conditions they are exposed to. As an important result, the lubricant is effective in the prevention of wear if it has the capacity to uptake and transmit electrical charges due to its polarisability during a wear process.

1. Introduction

1.1. A Few Aspects

The discipline of how power converts in a system of mating components with relative motion is the central topic of tribology. In a general view, at least two solid components are moved relatively against each other under force (the power input) and in the presence of a lubricant as a soft material in a randomly defined environment (see Figure 1).
Tribology appears as a scattered cross-discipline, merging physical sciences and chemistry while adapting to self-generated theories. However, a concise definition is hard to state, and it is commonly established (see, for example, [1,2,3,4,5] and many more). The friction and lubrication of solids (solid-state bodies) were studied previously by Bowden and Tabor [6]. Lubrication as a soft matter medium, separating the solid-state bodies, has been extensively studied within EHL (elastohydrodynamic lubrication) (see, for example, [7,8,9,10,11,12,13] and many more). Basically, EHL is the topic of how a lubricant, as a body of soft matter, enters and passes a converging gap between two solid bodies. It undertakes the effort to define the role of the hydrostatic pressure that lifts off the moving solid parts. This effect reduces the direct touch of the solid bodies might lead to the consequence of wear. Hence, EHL is related to the Reynolds theory of fluid transportation in pipes (see, for example, [14,15,16] and many more). Similar to modern drive train technology (e.g., automotive), lubrication films narrow down in order to reduce the friction caused by splash losses, and thin films in the region of sub-microns and nanoscales appear. Thin film lubrication (TFL) became a field of intense research (see, for example, [17,18,19,20] and more). While within EHL and TFL, a lubricant is treated as a continuum, the continuous reduction in the tribofilm into a molecular dimension (a few nanometers) means that their chemical properties, e.g., electron transfer and photon emission processes, emerge. These processes are the domain of Tribocatalysis (TC) (see, for example, [21] and more). Chemical processes related to catalysis are discussed as significant in understanding the role of functional additives in TFL and boundary lubrication. Plenty of studies have been carried out in order to investigate how the chemical structure is related to the function and the life of a drive train component (see, for example, [22] and more). Irregularities by the appearance of unexpectedly high wear and other effects caused by the use of additives that should prevent wear are reported (see, for example, [23,24] and more). Materials exposed to wear, where the solid bodies are touching each other, act under transient and non-equilibrium conditions (see, for example, [25,26,27] and more), emitting and exchanging photons, electrons, charged and non-charged particles (see, for example, [28,29,30] and more), denoted commonly as fracto-emission (FE). Lowering the film thickness in tribological contact brings both solid-state and lubricants under perpetual stress, causing the excitation of the molecular and solid-state structure. As the input frequency of stress is high compared to the relaxation and repairs the capability of the partner, electron and photon exchange processes appear, but also chemical reactions. The engineering of the transition from TFL to TC and FE is hard to parametrise by the components tested due to the many apparent scales (from a few nanometers to meters). Both engineering molecular processes and trans-scale modeling have given the motivation to undertake multi-scale modeling (see, for example, [17,31,32,33] and more).

1.2. White Etching Cracks (WEC) and Lubricants

WEC came into the interest of industry and research within more than a decade [34]. The failure has been attributed to material degradation processes (see, for example, [35,36] and more) in relation to lubricants and parametrised conditions (see, for example, [37,38] and more). Even though some patterns with respect to the compositions of lubricants are reported (see, for example, [39,40] and more), a clear pattern in relation to the structure of any chemical and stoichiometry is missing.

1.3. A Distinct Perspective

The continuous reduction in the tribogap (schematical sketch, see Figure 2) puts the lubricant constituents under stress. A (see Figure 2) reflects the situation where the thermal motion and mutual electronic excitation by neighbouring molecules are in balance with relaxation to the ground state. This basically reflects EHL and viscous flow. As their collision rates increase and the space becomes narrower, the molecules react with density and increase in symmetry, and their narrowing could lead to a superstructure (see Figure 2B). This self-organised structure acts as a new species, different in physics (e.g., A), solubility, viscosity, friction, and functional release. Narrowing the tribogap further (see Figure 2C) finally causes the disruption of the molecular system as the collision rates and the mutual touchdown do not allow structural relaxation.
We argue that the capability of a lubricant to undergo the transition from A → B → C is a key parameter in lubrication.

2. Materials and Methods

With reference to the basic investigations, as described in the previous section, we tried to set up an empirical-based concept on how the molecular transition from EHL toward TFL and TC could be described. As lubricants in the life cycle are a) constituted of many components and b) exposed to complex loading, we took the concept presented earlier [41,42,43]. Basically, the mutual interaction therein is described as a matter of a non-dimensional factor Di,j, which arises from the chemical potential of a pair i,j of constituents. Di,j could be interpreted as a lubricant stress factor caused by impact stress from the tribological gap (see Figure 2). This concept neglects the equilibrium that each molecular system tries to reach under stress but describes instead the proximity or distance to the equilibrium. Ideally, the equilibrium is reached if Di,j = 1 is the chemical potential of a set of equal components i,j, which, per definition, means Di,j = 1. As an activity parameter, we took, according to [41,42,43], the dipolemoment divided by the sum of all dipolemoments as well as the induced dipolemonent divided by the sum of all induced dipolemoments as non-dimensional parameters to calculate the stress response factor, Dij, between the components i and j. If Di,j < 1, it stays within {0,1} and denotes it’s probability to be present in the system. If, in contrast, Di,j > 1, we define the appearance of clusters by dividing Di,j by a number n (n ∈ Z), so that Di,j comes back into {0,1}. The number n is interpreted as degeneration, as the factor leading to clusters within the system of constituents. Notably, Di,j and the clusters appear simply as transient probabilities independent of how far they are from a thermodynamic equilibrium. Within a real lubricant, all combinations i,j have to be calculated by means of increasing activity. This leads to a matrix of cluster and non-cluster presence for each activity level on all possible combinations i,j related to the given activity. As all combinations are allowed per activity, the individual percentages of the dipole (cluster/non-cluster) and inducible dipoles (polarisabiliy) (cluster/non-cluster) are calculated. The dipole moments and polarisabilites are calculated using the semiempirical PM3 method (Hyperchen 8.0™). The chemical structures are derived either by brand name or by retro-analysing the lubricant with infrared and XRF analyses. All calculations described here are automised using the system software tool Solu42TM 1.1.

3. Experimental Section

In order to evaluate the system, an FE8 test rig (according) to the DIN 51819 standard was used (see Figure 3 and Figure 4 (taken from [39]).
A set of test lubricants divided into Extreme Pressure (EP) and Anti-Wear (AW) candidates given by the chemical structure of the additives was used (see [43]) and Appendix A. All permutations for any component pair i,j as a function of a virtual activity were calculated using Software Package Solu42TM. For the driving factor for the interaction of the components, their activity coefficient taking the relative dipolemoment was finally given by the dipolemoment of an individual divided by the sum of all components used. The same investigation was used to calculate the relative polarisabilities (the induced dipole moment). The dipolemoments and polarisabilities of all components are retrieved via a calculation using the semiempirical parametrised method 3 (PM3) with Hyperchem(8.0)™. As explicitly explained in [43], the factor Di,j is obtained as follows:
D i , j =   ln ( a i ) / ln ( a j )
(ai): Dipole activity for the relative dipolemoments and polarisability activity for the relative polarisability. The term Di,j is calculated using means in relation to increasing tribological activity (aj), which is basically the activity with respect to the surface (S) and the “inner” structure of the fluid (L).
If Di,j > 1, it is interpreted as a cluster (CL). Calculating all combinations, it comes down to the following predictors:
CL_Di_Add_S: Dipolar cluster formation of the additive at the surface
CL_Po_Add_S: Polarisability cluster formation of the additive at the surface
Di_Add_S: Dipolar formation of the additive at the surface
Po_Add_S: Polarisability formation of the additive at the surface
CL_Di_Add_L: Dipolar cluster formation of the additive within the liquid (L)
CL_Po_Add_L: Polarisability cluster formation of the additive within the liquid (L)
Di_Add_S_L: Dipolar formation of the additive within the liquid (L)
Po_Add_S_L: Polarisability formation of the additive within the liquid (L)
CL_Di_L_S: Dipolar cluster formation of the liquid at surface S
CL_Po_L_S: Polarisability cluster formation of the liquid at surface S
If more than one Baseoil component is present, e.g., Li,j, the combinations are extended for each i,j.
CL_Di_Add_Li,j: Dipolar cluster formation of the additive within the liquid Li,j
CL_Po_Add_Li,j: Polarisability cluster formation of the additive within the liquid Li,j
Di_Add_S Li,j: Dipolar formation of the additive within the liquid Li,j
Po_Add_S_ Li,j: Polarisability formation of the additive within the liquid Li,j
CL_Di_Li,j: Dipolar cluster formation of the additive within the liquid Li,j
CL_Po_Li,j: Polarisability Cluster formation of the Additive within the liquid Li,j
Di_Li,j_S: Dipolar formation of the liquid Li,j at surface S
Po_Li,j_S: Polarisability formation of the liquid Li,j at surface S
CL_Di_Li,j_S: Dipolar cluster formation of the liquid Li,j at surface S
CL_Po_Li,j_S: Polarisability cluster formation of the liquid Li,j at surface S

4. Results

The results of the component testing versus the variants in lubricants and their chemical predictors are exemplified in Table 1 for activity parameter 1, which represents the lowest activity. For the computation of the entirety of the chemical predictors, we used Hyperchem 8.0™ as a base in order to calculate the dipolemoment and the polarisability of each component. The computation of the chemical predictors as a function of imposed surface activity was carried out by the software package Solu42TM. The results were split into non-clusters and clusters (CL) for the various combinations given by the base oil and the additives. For computation, an increase in surface activity was used (for details, please contact the authors). In order to obtain the main predictors explaining the individual wear of the test rig components, e.g., cage, rollers, and washers, but also the average wear and the wear of single components, a variety of statistical methods were used.
The boxplot is a standardised graphical representation that effectively illustrates the distribution of data by utilising its five-number summary, which includes the minimum value, the first quartile (Q1), the median, the third quartile (Q3), and the maximum value. Boxplots provide valuable insights into the presence of outliers and their corresponding values, as well as the symmetry of data, the degree of data clustering, and the skewness of data [44]. We calculated the average wear for each lubricant to better represent its corresponding wear. Figure 5 illustrates that the average wear of AW has outliers (2 dots in red rectangle shape). After removing these outliers, a clearer representation of the distribution of the average wear for both EP and AW is shown in Figure 6. The spread of average wear for EP is noticeably wider than AW. Despite their medians being closely matched, EP displays a right-skewed distribution, while AW leans toward a left-skewed pattern.

4.1. Comparison of Wear Characteristics between EP and AW Lubricants

From Figure 7, overall, the median lines on the box plots suggest that EP lubricants tend to have higher wear attributes than AW across all four wear types. For the roller and cage types, AW data cluster tightly around the median and show much lower wear compared to EP. By contrast, for the housing and motorside types, the wear data for AW are more spread out, with the medians slightly different from EP.

4.1.1. EP Lubricant Wear Profile

Under EP lubrication, the wear values associated with the rollers, housing, and motorside mostly fall within the 1 to 200 range, while the wear cage displays range from 10 to approximately 500. This range has enhanced data dispersion that is considerably higher than other wear types within the EP category (Figure 8).

4.1.2. AW Lubricant Wear Profile

The wear characteristics of rollers, housing, and the motorside in AW lubricants typically fall between 0 and about 150. On the other hand, the wear cage varies from 30 to 350; however, this variation is not as pronounced as with EP lubricants. Notably, the wear roller in AW exhibits lower wear than the other types, as 50% of the data fall between 0 to approximately 70. The chemical codes AW6a and AW10a stand out as outliers since they fall beyond the average range indicated by the box plot’s whiskers, which is another important point to make (Figure 9).

4.2. Correlation between Average Wear and Chemical Predictors of EP and AW Lubricants

We used correlation analysis in our study to evaluate the linear relationships between chemical predictors and average wear. The correlation coefficients (r) were calculated to provide an indication of the strength and direction of these interactions. An r value near 1 or −1 indicates a strong positive or negative connection, respectively, whereas values closer to 0 imply a weaker link. Positive correlation indicates that one variable increases as the other increases, and negative correlation means that, as one variable increases, the other decreases.
A significance level (α) of 0.05 was used to examine the statistical significance of the observed connections. This criterion indicates that there is a 5% chance we could reject the null hypothesis of no correlation when it is true (Type I error). In this analysis, correlations with p-values less than 0.05 were considered statistically significant, implying strong evidence against the null hypothesis.
To emphasis the most significant connections discovered in our data, we selected highlighted correlations that exceeded a magnitude of 0.6, representing a moderately strong correlation.
For activities ranging from 1 to 10, both CLPo_Add_L1 and Po_L1_S in AW indicate a positive and moderately strong correlation with average wear. Meanwhile, other chemical predictors display only weak correlations. However, they are not statistically significant (p-value < 0.05) (Figure 10). While the correlation coefficient indicates a relationship, the absence of statistical significance shows that conclusions should be drawn with caution. This means that while the correlation value suggests a relationship, the results are not strong enough (at the 5% significance level) to conclude that it is not simply caused by random variability in the data. AW’s limited dataset size may have contributed to the high p-value.
In the EP group, activity 4 is the only activity where Po_L1_S exhibits a positive and moderately strong relationship. Activity 2 shows that CLDi_L1_S and Di_Add_L1 have a negative, moderately strong correlation with average wear. A positive, moderately strong correlation is evident between average wear and both Di_Add_S (for activities 15–80) and Po_L1_S1 (for activities 4 and 15–10,000). Conversely, both Di_Add_L1 and Di_L1_S show a negative, moderately strong correlation between the average wear in activities 2, 4, and 15–10,000 and 15–10,000, respectively. Importantly, all these observed correlations are statistically significant (Figure 11).
A positive relationship with average wear suggests that as the values of these chemical predictors increase, the average wear also tends to increase. On the other hand, a negative relationship with average wear suggests that, as these chemical predictors increase, the average wear tends to decrease.
Furthermore, the intermediate strength of these correlations shows that these relationships are notable but not exactly linear.
Finally, the fact that these correlations are statistically significant suggests that the observed relationships are not the result of random chance. The data show sufficient evidence against the null hypothesis.

4.3. Relationship of Average Wear and Mean of Cluster and Non-Cluster Chemical Predictors of EP and AW Lubricants:

In the EP category in Figure 12, we computed the average for the attributes of both cluster and non-cluster chemical predictors for each chemical code. When plotted against their average wear, the following distinct patterns emerged:
  • The cluster and non-cluster groups scatter in contrasting directions, almost as if they mirror each other.
  • There is a dramatic shift in patterns from activity 1 to activity 15. Beyond that, from activity 20 to 10,000, the changes are subtle.
  • In activity 1, the non-cluster group stays on the left side of the graph. As activities progress, it moves gradually to the right. Conversely, the cluster group follows mirrored behaviour, initially starting on the right and moving leftward. The side switch can be observed from activity 7.
  • Low-wear and high-wear groups can be observed in both cluster and non-cluster groups where data points above the 100 mean wear are considered high wear and below this is low wear.
In the AW category in Figure 13, we applied the same method as in EP and obtained the same pattern. However, there were some differences in AW.
  • There is a dramatic shift in patterns in activities 1, 10 and 15.
  • The side switch starts in activity 10, which is one activity after EP.

4.4. Relationship of Average Wear and Chemical Predictors of EP and AW Lubricants:

In Figure 14, Di_Add_L1 exclusively resides in the low-wear zone and exhibits a negative correlation with average wear. By contrast, other chemical predictors are present in both the high- and low-wear zones. Notably, ‘CLDi_Add_S’ is observed in activities 1, 2, 4, and 7, while both ‘CLPo_Add_S’ and ‘CLDi_L1_S’ appear in activities 1, 2, 4, 7, and 10.
For AW (Figure 15), Di_Add_L1 had a similar trend to EP, which exclusively resides in the low-wear zone with a negative relationship with the average wear except in activities 2,15,25, and 30, which show a positive relationship. Additionally, CLDi_Add_L1 only appears in the high-wear and AW group.
The differences in oil codes between the low- and high-wear groups are explained in Figure 16. The blue and green oil codes, which stand for Trimethylolepropaneoctoate (TMP) and PAO8 + TMP, respectively, are most frequently found in the high-wear category. The presence of the orange, purple, and red oil codes, which stand for PAO8, PAO40 + PAO8, and PAO40 + TMP, respectively, distinguished the low-wear group by contrast. Notably, CLDi_L1_S is primarily found in the low-wear zone, while Di_Add_L1 is only present there. Nevertheless, both wear zones are covered by chemical predictors, such as CLDi_Add_S, CLPo_Add_S, CLDi_Add_L1, CLPo_Add_L1, Di_Add_S, Po_Add_S, Di_L1_S, and Po_L1_S.
Figure 17 shows that for AW, TMP and PAO 8 are the only two oil codes in AW. In particular, the connection between Di_Add_L1 and CLDi_L1_S is constant and is only observed in the low-wear category. This suggests that increasing the predictor has no effect on the average wear; instead, it stays the same.

4.5. Relationship of Average Wear and Chemical Predictors of AW Lubricants with outliers:

On the other hand, if outliers are taken into account, Figure 18 shows that Di_Add_L1 and CLDi_L1_S appear in the high-wear range with the outliers, although they still show the same consistent association.
The correlation coefficient for the AW group, as shown in Figure 19, highlights a difference when outliers are taken into account: the correlations between average wear and chemical predictors for AW and EP shift in the reverse directions, but AW still shows no statistically significant result.

4.6. Logistics Regression to Observe Relationship of Wear Category and Chemical Predictors

As we classified the wear into low and high categories, we utilised g logistic regression to determine how chemical predictors influence these distinct wear levels. Notably, the logistic regression analysis is exclusively applied to the EP group, as AW did not yield statistically significant results.
To address issues with multicollinearity, we chose to include only three predictors (CLPo_Add_L1, Di_L1_S, and Di_Add_L1) in our model. These three predictors were chosen due to their low variance inflation factors (VIFs) (Figure 20), which indicates insignificant collinearity issues. In addition, their correlation coefficients with the wear category are the highest compared to those of the other variables (Figure 21).
This logistic regression model shows that “CLPo_Add_L1” and “Di_Add_L1” have statistically significant positive and negative influences, respectively, on the probability of being in the high wear category, while “Di_L1_S” has a statistically significant negative influence (Figure 22).
For each unit increase in “CLPo_Add_L1,” log-odds of being in the high wear category increased by 1.9086 units (p < 0.001). An increase in one unit in “Di_L1_S” was associated with a decrease of approximately 0.0725 units in log-odds of the high wear category (p < 0.05). Conversely, for every unit increase in “Di_Add_L1,” log-odds of being in the high wear category decreased by approximately 0.3294 units (p < 0.001).
These findings illuminate the relationships between these chemical predictors and wear patterns. The statistically significant coefficients indicate that “CLPo_Add_L1” positively influences high wear, whereas “Di_L1_S” and “Di_Add_L1” negatively influence high wear in EP lubricants. “CLPo_Add_L1” had a statistically significant positive coefficient, showing that higher quantities of this chemical predictor are related to an increased chance of high wear in EP lubricants. By contrast, the predictors “Di_L1_S” and “Di_Add_L1” had statistically significant negative coefficients, indicating that higher quantities of these chemical predictors are associated with a reduced likelihood of high wear.
In conclusion, our logistic regression analysis reveals the significant influence of these chemical predictors on wear categories, casting light on the factors that contribute to wear patterns in our study.

4.7. Other Significant Findings

The AW group, without outliers, only shows a significant correlation in the wear roller. CLPo_Add_L1 shows a strong positive relationship with the wear roller for all activities. By contrast, ‘Po_L1_S’ demonstrates this strength of association only in activities 1–15, whereas ‘CLDi_Add_S’ demonstrates a similar strength of association only in activity 10 (Figure 23).
Given the high correlation observed with CLPo_Add_L1, implying a strong linear relationship, we proceeded to conduct a linear regression analysis with this predictor (see Figure 24). Additionally, we included Po_L1_S as it contributed to an improved R-squared value. The results of this regression analysis yielded an adjusted R-squared of 0.67. This figure indicates that the selected independent variables in the model, such as “CLPo_Add_L1” and “Po_L1_S,” explain approximately 66% of the variation in the dependent variable, “Wear_Roller.” Notably, some prospective predictors were excluded from the model on purpose. Their lack of statistically significant relationships with the dependent variable and their inability to enhance the adjusted R-squared value led to this conclusion.
Interpreting the coefficients, for every one unit of change in CLPo_Add_L1, there was an estimated change of 17.0164 units in the wear roller. Similarly, a one-unit change in Po_L1_S corresponded to an estimated change of 4.33 units in the wear roller. These findings illuminate the impact of both CLPo_Add_L1 and Po_L1_S on the wear roller in the AW group.

5. Conclusions

Commonly accepted hypotheses in tribology state that, in the case of boundary lubrication, functional additives prevent direct metal–metal contact in the contact zone of the machinery element. Traditionally, boundary lubrication is associated with the relation of the film height in the contact toward the height of the asperities. The film height is calculated traditionally using the theory of elastohydrodynamics (EHL), stating that the viscosity (e.g., film thickness) locally present in the contact zone is the prime predictor defining a boundary situation. Classical EHL states that viscosity, as the relevant predictor, transforms using the local energy impact in the contact zone given by temperature, temperature flashes, and shear impact, both leading to a decay in the viscosity. By contrast, the pressure assumed to be locally present in the contact zone compresses the lubricant, which is associated with an increase in viscosity in the lubricant by means of a viscosity-pressure coefficient. A lot of computational effort is carried out in order to obtain an insight into these processes and the structure–property relationship given by the chemical structure of the lubricant and their ingredients. However, up until now, no specific rules have been seen regarding how wear processes are assumed to prevent, but in some cases, are also found to promote and are predictable. The transit of EHL toward thin film lubrication (TFL) and tribocatalysis (TC) is accompanied by various physical transitions within the lubricant, notified as their “inner” structure herein. The exposure to the impact of tribological stress causes the lubricant to transform its spatial solute structure before it starts to react chemically. These structural changes are described within this paper by calculating the chemical potential, differentiated in the dipolar and induced dipolar (polarisability) activity appearing at the surface (S) and in the solvent (L). Solutes in solvents and at surfaces may create clusters under stress in order to dissipate it. This study permuted all different structural rearrangements of the lubricants in relation to imposed external activity. This activity could arise from temperature, pressure, and surfaces becoming active due to the tribological impact. The study here applies the assumption comprising the formation of “inner” structures, clusters (CL), and non-clusters to a variety of lubricants with well-known chemical structures and their ability to prevent or induce wear using the bearing test rig FE8 (DIN 51819). The lubricants were divided into two different base oils (Poly-a-Olefine (PAO) and an Ester (Trimethylpropane Octoate (TMP) and were placed under different extreme pressure (EP) and anti-wear (AW) additives. The wear within this test is scattered over the single constituents, e.g., cage, rollers, and washer. For all the different lubricants, the “inner” structure was calculated as dipolar and induced dipolar activity and their likelihood to form clusters at the surface and within the base oil. A detailed statistical investigation, searching for relations of the “inner” structure with individuals, but also the average wear led to the following results:
The base oil, rather than the additives, significantly guides the appearance of high or low wear.
The lubricant predictors change abruptly due to the increase in the surface activity. This leads to the conclusion that determining the activity of a lubricant during any test requires a careful definition of how the load capacity of the component is related to the load capacity of the lubricant. This is to say, the result of any load-related tribological test may abruptly change under different test conditions.
The encapsulation of induced dipole moments (clustered polarisability) and additives in a solvent are relevant as clusters promote significant wear. We are interpreting this result by the fact that a wear process is highly transient by exchanging electrons, charged and uncharged particles. These transient exchange processes are scavenged if the lubricant uptakes them via charge dissipation due to its polarisation ability. If the “inner” structure of the lubricant clusters the polarisability within the solvent, this functionality is blocked, and the wear process is not dissipated in the lubricant.
The presence of outliers in the AW dataset, combined with its small size, suggests the unfavourable effect on statistical metrics in smaller datasets, thus misleading the interpretation of the results. The situation emphasises the importance of dataset size and quality in tribological research. It raises the delicate problem of how to approach outliers, particularly in smaller datasets, where their influence can be clearer. The topic for discussion is whether outliers should be thoroughly accounted for in order to maintain the dataset’s accuracy or whether they should be set aside because they may represent unique but real-world conditions. This opens the door to a broader discussion about improving data-gathering procedures in tribology. Future research could focus on increasing the AW dataset to either confirm or refine the current study’s conclusions while carefully examining the influence of outliers in the analysis.

6. Summary

The transition from EHL toward TFL and TC requires a detailed understanding of how lubricants uptake the exposition to contact stress. Apart from viscosity-related parameters, the “inner” structure of a lubricant comes into the picture. Despite the huge variety of chemicals, the simple calculation of the dipolar and induced dipolar structure leads to powerful (e.g., robust) predictors for wear under thin film lubrication (TFL) conditions. Despite the fact that the chemical structures had to be assumed due to their technical grade, the prediction is seemingly very robust if the extreme pressure additives are taken. This insight is even more important as new lubricants with respect to the prediction of wear have to be found for E-mobility and ecology purposes.

Author Contributions

Conceptualization, W.H. and L.K.; methodology, W.H. and L.K.; software, W.H. and L.K.; validation, W.H. and B.G.; formal analysis, L.K.; investigation, W.H.; resources, W.H.; data curation, L.K.; writing—original draft preparation, L.K.; writing—review and editing, W.H. and L.K.; visualization, W.H. and L.K.; supervision, J.F.; project administration, J.F.; funding acquisition, B.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy restrictions.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

List of Additives
EP
  • Triphenylphosphorothionate
  • Tridodecyltrisulfide
  • Sulfurated Isobutene S2Cl2 Route
  • Sulfurated Isobutene NaSx Route
  • Sulfurated Rapseed Oil
  • Dithiophosphoric Acid Isopropylester
AW
  • Methylenebisdithiocarbamate
  • Molybdenumbisdithiocarbamate
  • Zinc-2-ethylhexyldithiophosphate
  • Molybdenum-2-ethylhexyldithiophosphate
  • Tricresylphosphate
  • 2-Ethylhexyldithiophosphate-Amine neutralised
  • Bis-C13-alkylammonium -dihexylphosphate
  • 2.6-Dimercaptodithiazole

References

  1. Bushan, B. Introduction to Tribology; Wiley: Hoboken, NJ, USA, 2013; ISBN 978-1-119-94453-9. [Google Scholar]
  2. Moore, D.F. Principals and Applications of Tribology; Pergamon Press: Elmsford, UK, 1975; ISBN 0-08-0179-020. [Google Scholar]
  3. Shizhu, W.; Ping, H. Principles of Tribology; John Wiley & Sons: Hoboken, NJ, USA; Tsinghua University Press: Beijing, China, 2012. [Google Scholar]
  4. Deters, L. Reibung, Verschleiß und Schmierung. In Konstruktionselemente des Maschinenbaus 2; Sauer, B., Ed.; Springer: Berlin/Heidelberg, Germany, 2008; pp. 1–68. ISBN 978-3-540-76653-7. [Google Scholar]
  5. Shi, J.; Zhao, B.; He, T.; Tu, L.; Lu, X.; Xu, H. Tribology and dynamic characteristics of textured journal-thrust coupled bearing considering thermal and pressure coupled effects. Tribol. Int. 2023, 180, 108292. [Google Scholar] [CrossRef]
  6. Bowden, F.P. The Friction And Lubrication of Solids. In Oxford Classsic Texts in the Physical Sciences; Oxford University Press: Oxford, UK, 2001; ISBN 13-978-0198507772. [Google Scholar]
  7. Priest, M.; Ehret, P.; Flamand, L.; Dalmaz, G.; Childs TH, C.; Dowson, D.; Berthier, Y.; Taylor, C.M.; Lubrecht, A.A. Lubrication at the Frontier: The Role of the Interface and Surface Layers in the Thin Film and Boundary Regime, 1st ed.; Elsevier Science: Amsterdam, The Netherlands, 1999; ISBN 9780080535661. [Google Scholar]
  8. Zhu, D.; Wang, Q.J. Elastohydrodynamic Lubrication: A Gateway to Interfacial Mechanics—Review and Prospect. J. Tribol. 2011, 133, 041001. [Google Scholar] [CrossRef]
  9. Wills, G. Lubrication Fundamentals (Mechanical Engineering); Marcel Dekker: New York, NY, USA, 1980. [Google Scholar]
  10. Ginsberg, J.H. Advanced Engineering Dynamics; Cambridge University Press: Cambridge, UK, 1995; ISBN 13 978-0521470216. [Google Scholar]
  11. Zhang, Y.; Biboulet, N.; Venner, C.H.; Lubrecht, A. Prediction of the Stribeck curve under full-film Elastohydrodynamic Lubrication. Tribol. Int. 2020, 149, 105569. [Google Scholar] [CrossRef]
  12. Bi, Z.; Mueller, D.W.; Zhang, W. State of the art of friction modelling at interfaces subjected to elastohydrodynamic lubrication (EHL). Friction 2020, 9, 207–227. [Google Scholar] [CrossRef]
  13. Ewen, J.P.; Spikes, H.A.; Dini, D. Contributions of Molecular Dynamics Simulations to elastohydrodynamic Lubrication. Tribol. Lett. 2021, 69, 24. [Google Scholar] [CrossRef]
  14. Hamrock, B.J.; Steven, R.; Schmid, S.R.; Jacobson, B.O. Fundamentals of Fluid Film Lubrication, 2nd ed.; CRC Press Inc.: Boca Raton, FL, USA, 2004; ISBN 13 978-0824753719/10 0824753712. [Google Scholar]
  15. Szeri, A.Z. Fluid Film Lubrication: Theory and Design, 2nd ed.; Cambridge University Press: Cambridge, UK, 2011; ISBN 13 978-0521898232/10 0521898234. [Google Scholar]
  16. Zakharov, S. Hydrodynamic lubrication research: Current situation and future prospects. J. Frict. Wear 2010, 31, 56–67. [Google Scholar] [CrossRef]
  17. Tan, Y.; Huang, W.; Wang, X. Tribochemistry of ZDDP in molecular orbital calculations. Tribol. Int. 2004, 37, 447–450. [Google Scholar] [CrossRef]
  18. Khonsari, M.M.; Booser, R.E. Applied Tribology: Bearing Design and Lubrication; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2008; ISBN 978-0-470-05711-7. [Google Scholar]
  19. Kudish, I.I. Elastohydrodynamic Lubrication for Line and Point Contacts; CRC Press: Boca Raton, FL, USA, 2013; ISBN 978-1-4665-8390-0. [Google Scholar]
  20. Spikes, H.A. Triboelectrochemistry: Influence of applied electrical potentials on friction and wear of lubricated contacts. Tribol. Lett. 2020, 68, 90. [Google Scholar] [CrossRef]
  21. Li, X.; Tong, W.; Shi, J.; Chen, Y.; Zhang, Y.; An, Q. Tribocatalysis mechanisms: Electron transfer and transition. J. Mater. Chem. A 2023, 11, 4458–4472. [Google Scholar] [CrossRef]
  22. Holweger, W. Novel Predictors for Friction and Wear in Drivetrain Applications. In Friction, Lubrication and Wear; Chowdhury, M.A., Ed.; Intechopen: London, UK, 2019. [Google Scholar]
  23. Sommer, K.; Heinz, R.; Schöfer, J. Verschleiß Metallischer Werkstoffe; Springer Vieweg: Wiebaden, Germany, 2014; ISBN 978-3-8348-2464-6. [Google Scholar]
  24. Hostis, B.L.; Minfray, C.; Frégonèse, M.; Verdu, C.; Ter-Ovanessian, B.; Vacher, B.; Le Mogne, T.; Jarnias, F.; D’Ambros, A.D.-C. Influence of lubricant formulation on rolling contact fatigue of gears—Interaction of lubricant additives with fatigue cracks. Wear 2017, 382–383, 113–122. [Google Scholar] [CrossRef]
  25. Kovalchenko, M.V.; Korchak, S.N.; Kuznetsov, O.L.; Erdemir, A.G.V.; Samsonov, G.V. The effects of mechanical stresses on catalytic processes. J. Catal. 2002, 212, 232–235. [Google Scholar]
  26. Erdemir, A. Tribocatalysis: A new field at the frontier of tribology and catalysis. Surf. Sci. 2014, 631, 2–8. [Google Scholar]
  27. Chen, Z.; Zhang, R.; Gao, Y.; Yang, Q.; Liu, J. The role of catalysts in friction and wear. Friction 2019, 7, 271–293. [Google Scholar]
  28. Dickinson, T.; Jensen, L.C.; Lee, S.; Scudiero, L.; Langford, S.C. Fracto-emission and electrical transients due to interfacial failure. J. Adhes. Sci. Technol. 1994, 8, 1285–1309. [Google Scholar] [CrossRef]
  29. Dickinson, J.T.; Jensen, L.C.; Jahan-Latibari, A. Fracto-emission: The role of charge separation. J. Vac. Sci. Technol. A 1984, 2, 1112–1116. [Google Scholar] [CrossRef]
  30. Kajdas, C.K. Importance of the triboemission process for tribochemical reaction. Tribol. Int. 2005, 38, 337–353. [Google Scholar] [CrossRef]
  31. Huang, W.; Tan, Y.; Chen, B.; Dong, J.; Wang, X. The binding of antiwear additives to iron surfaces: Quantum chemical calculations and tribological tests. Tribol. Int. 2003, 36, 163–168. [Google Scholar] [CrossRef]
  32. Jayadas, N.H.; Prabhakaran Nair, P. Study of the Anti-Wear Properties of Coconut Oil Using Quantum Chemical Calculations and Tribological Tests. ASME J. Tribol. 2006, 128, 654–659. [Google Scholar] [CrossRef]
  33. Vakis, A.; Yastrebov, V.; Scheibert, J.; Nicola, L.; Dini, D.; Minfray, C.; Almqvist, A.; Paggi, M.; Lee, S.; Limbert, G.; et al. Modeling and simulation in tribology across scales: An overview. Tribol. Int. 2018, 125, 169–199. [Google Scholar] [CrossRef]
  34. Evans, M.-H. An updated review: White etching cracks (WECs) and axial cracks in wind turbine gearbox bearings, Materials Science and Technology. Mater. Sci. Technol. 2016, 32, 1133–1169. [Google Scholar] [CrossRef]
  35. Šmeļova, V.; Schwedt, A.; Wang, L.; Holweger, W.; Mayer, J. Microstructural changes in White Etching Cracks (WECs) and their relationship with those in Dark Etching Region (DER) and White Etching Bands (WEBs) due to Rolling Contact Fatigue (RCF). Int. J. Fatigue 2017, 100, 148–158. [Google Scholar] [CrossRef]
  36. Spille, J.; Wranik, J.; Barteldes, S.; Mayer, J.; Schwedt, A.; Zürcher, M.; Lutz, T.; Wang, L.; Holweger, W. A study on the initiation processes of white etching cracks (WECs) in AISI 52100 bearing steel. Wear 2021, 477, 203864. [Google Scholar] [CrossRef]
  37. Holweger, W.; Schwedt, A.; Rumpf, V.; Mayer, J.; Bohnert, C.; Wranik, J.; Spille, J.; Wang, L. A Study on Early Stages of White Etching Crack Formation under Full Lubrication Conditions. Lubricants 2022, 10, 24. [Google Scholar] [CrossRef]
  38. Esmaeili, K.; Wang, L.; Harvey, T.J.; White, N.M.; Holweger, W. A Study on the Influence of Electrical Discharges on the Formation of White Etching Cracks in Oil-Lubricated Rolling Contacts and Their Detection Using Electrostatic Sensing Technique. Lubricants 2023, 11, 164. [Google Scholar] [CrossRef]
  39. Wranik, J.; Holweger, W.; Lutz, T.; Albrecht, P.; Reichel, B.; Wang, L. A Study on Decisive Early Stages in White Etching Crack Formation Induced by Lubrication. Lubricants 2022, 10, 96. [Google Scholar] [CrossRef]
  40. Gould, B.; Demas, N.; Erck, R.; Lorenzo-Martin, M.C.; Ajayi, O.; Greco, A. The Effect of Electrical Current on Premature Failures and Microstructural Degradation in Bearing Steel. Int. J. Fatigue 2021, 145, 106078. [Google Scholar] [CrossRef]
  41. Shivakumar, D.; Williams, J.; Wu, Y.; Damm, W.; Shelley, J.; Sherman, W. Prediction of Absolute Solvation Free Energies using Molecular Dynamics Free Energy Perturbation and the OPLS Force Field. J. Chem. Theory Comput. 2010, 6, 1509–1519. [Google Scholar] [CrossRef]
  42. Mezei, M. The finite difference thermodynamic integration, tested on calculating the hydration free energy difference between acetone and dimethylamine in water. J. Chem. Phys. 1987, 86, 7084–7088. [Google Scholar] [CrossRef]
  43. Holweger, W.; Bobbio, L.; Mo, Z.; Fliege, J.; Goerlach, B.; Simon, B. A Computational Study on the Role of Lubricants under Boundary Lubrication. Lubricants 2023, 11, 80. [Google Scholar] [CrossRef]
  44. Albright, S.C.; Winston, W.L. Business Analytics: Data analysis & decision making. In Describing the Distribution of a Single Variable; Cengage Learning, Inc.: Beijing, China, 2016; pp. 49–51. [Google Scholar]
Figure 1. Abstract representation of a Tribosystem.
Figure 1. Abstract representation of a Tribosystem.
Lubricants 11 00514 g001
Figure 2. Reducing the tribogap from (A) (EHL) to (B) (TFL) to (C) (TC) as a scheme. Red dots shall mark cluster formation as consequence of spatial stress.
Figure 2. Reducing the tribogap from (A) (EHL) to (B) (TFL) to (C) (TC) as a scheme. Red dots shall mark cluster formation as consequence of spatial stress.
Lubricants 11 00514 g002
Figure 3. Structure of the FE8 test rig: 1. Frame and crane for the test head. 2. Torque transducer 3. Thermocouples for each bearing. 4. Test head. 5. Acceleration sensor 6. Driving unit 7. Control cabinet for terminals for cooling and fan heating. (Source [39]).
Figure 3. Structure of the FE8 test rig: 1. Frame and crane for the test head. 2. Torque transducer 3. Thermocouples for each bearing. 4. Test head. 5. Acceleration sensor 6. Driving unit 7. Control cabinet for terminals for cooling and fan heating. (Source [39]).
Lubricants 11 00514 g003
Figure 4. Test head with axial cylindrical roller bearings. 1. Housing. 2. Test bearing 2 (test head-sided bearing). 3. Spacer 4. Test bearing 1 (motor sided bearing). 5. Shaft. 6. Clamping bolt. 7. Bearing seat. 8. Drainpipe. 9. Cap. 10. Bearing support with screwed-on pilot pin. 11. Lid cup of spring package. 12. Lid. 13. Auxiliary bearing. Both test bearings consist of stationary housing (HW) and a rotating shaft washer (SW). Figure adapted from Source [39] and [43]).
Figure 4. Test head with axial cylindrical roller bearings. 1. Housing. 2. Test bearing 2 (test head-sided bearing). 3. Spacer 4. Test bearing 1 (motor sided bearing). 5. Shaft. 6. Clamping bolt. 7. Bearing seat. 8. Drainpipe. 9. Cap. 10. Bearing support with screwed-on pilot pin. 11. Lid cup of spring package. 12. Lid. 13. Auxiliary bearing. Both test bearings consist of stationary housing (HW) and a rotating shaft washer (SW). Figure adapted from Source [39] and [43]).
Lubricants 11 00514 g004
Figure 5. Box plot comparing the average wear between EP and AW (excluding outliers).
Figure 5. Box plot comparing the average wear between EP and AW (excluding outliers).
Lubricants 11 00514 g005
Figure 6. Box plot comparing the average wear between EP and AW.
Figure 6. Box plot comparing the average wear between EP and AW.
Lubricants 11 00514 g006
Figure 7. Box plot comparing the wear components between EP and AW.
Figure 7. Box plot comparing the wear components between EP and AW.
Lubricants 11 00514 g007
Figure 8. Box plot comparing the wear components of EP.
Figure 8. Box plot comparing the wear components of EP.
Lubricants 11 00514 g008
Figure 9. Box plot comparing the wear components of AW with outliers (dots in red rectangle shape) (left) and without outliers (right).
Figure 9. Box plot comparing the wear components of AW with outliers (dots in red rectangle shape) (left) and without outliers (right).
Lubricants 11 00514 g009
Figure 10. Correlation coefficient between average wear and chemical predictors of AW across different chemical activities (excluding outliers).
Figure 10. Correlation coefficient between average wear and chemical predictors of AW across different chemical activities (excluding outliers).
Lubricants 11 00514 g010
Figure 11. Correlation coefficient between the average wear and chemical predictors of EP across different chemical activities.
Figure 11. Correlation coefficient between the average wear and chemical predictors of EP across different chemical activities.
Lubricants 11 00514 g011
Figure 12. Scatter plot of average mean and average attributes of cluster and non-cluster chemical predictors across chemical activity 1–40.
Figure 12. Scatter plot of average mean and average attributes of cluster and non-cluster chemical predictors across chemical activity 1–40.
Lubricants 11 00514 g012
Figure 13. Scatter plot of average mean and average attributes of cluster and non-cluster chemical predictors in each chemical activity (AW—excluding outliers).
Figure 13. Scatter plot of average mean and average attributes of cluster and non-cluster chemical predictors in each chemical activity (AW—excluding outliers).
Lubricants 11 00514 g013
Figure 14. Scatter plot for average wear and chemical predictors in each chemical activity (EP).
Figure 14. Scatter plot for average wear and chemical predictors in each chemical activity (EP).
Lubricants 11 00514 g014
Figure 15. Scatter plot for average wear and chemical predictors in each chemical activity (AW) excluding outliers.
Figure 15. Scatter plot for average wear and chemical predictors in each chemical activity (AW) excluding outliers.
Lubricants 11 00514 g015
Figure 16. Scatter plot of average wear and chemical predictors grouped by oil code (EP).
Figure 16. Scatter plot of average wear and chemical predictors grouped by oil code (EP).
Lubricants 11 00514 g016
Figure 17. Scatter plot of average wear and each chemical predictor (AW) grouped by oil code (excluding outliers).
Figure 17. Scatter plot of average wear and each chemical predictor (AW) grouped by oil code (excluding outliers).
Lubricants 11 00514 g017
Figure 18. Scatter plot of average wear and each chemical predictor (AW) grouped by oil code (including outliers).
Figure 18. Scatter plot of average wear and each chemical predictor (AW) grouped by oil code (including outliers).
Lubricants 11 00514 g018
Figure 19. Correlation coefficient of average wear and chemical predictors of EP (left) and AW (right) across chemical activity (including outliers).
Figure 19. Correlation coefficient of average wear and chemical predictors of EP (left) and AW (right) across chemical activity (including outliers).
Lubricants 11 00514 g019
Figure 20. VIF factor.
Figure 20. VIF factor.
Lubricants 11 00514 g020
Figure 21. Correlation matrix (red for negative and green for positive, the intensity represents the magnitude).
Figure 21. Correlation matrix (red for negative and green for positive, the intensity represents the magnitude).
Lubricants 11 00514 g021
Figure 22. Logistic regression result (EP).
Figure 22. Logistic regression result (EP).
Lubricants 11 00514 g022
Figure 23. Correlation of wear cage and EP predictors (left) and wear roller and AW predictors excluding outliers (right).
Figure 23. Correlation of wear cage and EP predictors (left) and wear roller and AW predictors excluding outliers (right).
Lubricants 11 00514 g023
Figure 24. Linear regression wear roller—AW excluding outliers.
Figure 24. Linear regression wear roller—AW excluding outliers.
Lubricants 11 00514 g024
Table 1. Testing results (example from virtual activity 1).
Table 1. Testing results (example from virtual activity 1).
Lubricant
Code
FrictionWear ComponentsPredictors
Wear
Roller
(mg.)
Wear
Cage
(mg.)
Wear Housing
(mg.)
Wear Motorside
(mg.)
Di_Add_SCLDi_Add_SPo_Add_SCLPo_Add_SDi_L1_SCLDi_L1_SPo_L1_SDi_Add_L1CLDi_Add_L1CLPo_Add_L1
AW131.50.5340.51014.45016.8737.01010.64016.554.48
AW1024.2530.5221103.25100.75013.63015.0635.38014.62015.425.89
AW10a57.2520951030.751252.751402.25031.59011.98014.5712.4823.2106.17
AW1118131.5138.25141140.25015.79016.929.21016.92012.189
AW11a24.2536.75149.536.2550.5031.17011.78012.8112.6425.6905.9
AW1a33035600033.39013.37015.329.1224.6604.15
AW631.568.595.75162.25168.25014.07015.9135.5012.62016.335.57
AW6a65.75624.25386963.5859.5032.81012.76014.9810.924.3604.19
EP 14c19.5462430.525.5025.05010.1927.5607.06026.613.52
EP 14d19.54.7510.53.252.75028.94012.11028.478.63017.764.09
EP 1d2428.558.251919.5019.86023.6021.938.2321.7904.6
EP 2b31.32201.25415.75216219.75022.94026.72014.398.99021.65.37
EP 624.8122.75265.25153.75150029.41029.2613.96014.4306.596.35
EP 6a 9378.258.54.75022.35025.1609.6313.2421.0108.6
EP 6e22.5109.75275.25108.25134.25018.05017.9721.57022.308.2311.88
EP 6f261.75255.2531.25029.52010.87010.5914.3927.8406.79
EP122.4173.5387.5111.5112.75026.63026.7414.56014.1908.819.08
EP1420.06278259137.25249.25170.75172.5030.93031.1812.38011.808.045.67
EP14a30.18226616727.755.755.75023.07026.8709.4410.9122.8106.9
EP1a29.6035229897.2553.55762.25022.44025.04010.0914.1520.2308.05
EP1b23.78140.75480.75174.5150.75021.63023.81017.3813.36016.137.7
EP1c18.594149.5119.2599.75022.97026.81015.29.11020.575.34
EP1e2191.25368.75132148.5016.83016.8823.13022.53010.1710.47
EP1f25.25138.7557.569.2589.75029.61010.82011.0815.3726.7706.35
EP220.53146015171.75383.75212.5186033.15033.4611.1010.5306.695.08
EP2a34.00044352132.2541.255960.5019.83023.69023.078.0320.6604.73
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Goerlach, B.; Holweger, W.; Kitirach, L.; Fliege, J. Predicting Wear under Boundary Lubrication: A Decisive Statistical Study. Lubricants 2023, 11, 514. https://doi.org/10.3390/lubricants11120514

AMA Style

Goerlach B, Holweger W, Kitirach L, Fliege J. Predicting Wear under Boundary Lubrication: A Decisive Statistical Study. Lubricants. 2023; 11(12):514. https://doi.org/10.3390/lubricants11120514

Chicago/Turabian Style

Goerlach, Bernd, Walter Holweger, Lalita Kitirach, and Joerg Fliege. 2023. "Predicting Wear under Boundary Lubrication: A Decisive Statistical Study" Lubricants 11, no. 12: 514. https://doi.org/10.3390/lubricants11120514

APA Style

Goerlach, B., Holweger, W., Kitirach, L., & Fliege, J. (2023). Predicting Wear under Boundary Lubrication: A Decisive Statistical Study. Lubricants, 11(12), 514. https://doi.org/10.3390/lubricants11120514

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop