0

### Research Papers

J. Tribol. 2013;135(2):021001-021001-11. doi:10.1115/1.4023410.

The effect of single-sided and double-sided harmonic surface waviness on the film thickness, pressure, and temperature oscillations in an elastohydrodynamically lubricated eccentric-tappet pair has been investigated in relation to the eccentricity and the waviness wavelength. The results show that, during one working cycle, the waviness causes significant fluctuations of the oil film, pressure, and temperature, as well as a reduction in minimum film thickness. Smaller wavelength causes more dramatic variations in oil film. The fluctuations of the pressure, film thickness, temperature, and traction coefficient caused by double-sided waviness are nearly the same compared with the single-sided waviness, but the variations are less intense.

Commentary by Dr. Valentin Fuster

### Research Papers: Applications

J. Tribol. 2013;135(2):021101-021101-12. doi:10.1115/1.4023147.

Novel Al-6061 hybrid composite materials with varying range of SiC particulate and constant weight percentage of B4C particulate and 1% of magnesium alloy were fabricated by the stir casting technique (liquid metallurgy route). The mechanical and tribological properties of the hybrid composites and that of Al-6061 unreinforced alloy were examined by a Rockwell hardness test machine, pin-on-disc test machine, and Optical Microscope and Scanning Electron Microscope (SEM). A plan of experiment generated through Taguchi's technique was used to conduct experiments based on the L27 orthogonal array. The developed analysis of variance (ANOVA) and the regression equations were used to investigate the influence of parameters like sliding speed, applied load, sliding time, and percentage of reinforcement and their interactions on the dry sliding wear and friction coefficient of the composites. The motto of the present study is ‘the smaller the better’ to identify the optimum conditions for dry sliding wear and the friction coefficient. The results indicate that wear rate and friction coefficient were highly influenced by applied load, sliding speed, percentage of reinforcement and sliding time whereas the interaction between these parameters show only a minor influence in Hybrid Metal Matrix Composites (HMMCs). The wear surface morphology and wear mechanism of the pins were investigated using SEM and were correlated with wear test results. Finally, confirmation tests were carried out to verify the experimental results. It is concluded that Al-6061 hybrid composite can replace the conventional material used in the brake disc of automobiles owing to improved hardness and strength and reduced wear rate.

Commentary by Dr. Valentin Fuster

### Research Papers: Contact Mechanics

J. Tribol. 2013;135(2):021401-021401-10. doi:10.1115/1.4023409.

The evolution of the maximum contact stresses in amorphous carbon coated silicon during sliding wear against a Si3N4 ball was investigated. Amorphous carbon coating was prepared on a silicon substrate by the electron cyclotron resonance (ECR) plasma sputtering method. Surface morphologies of the coating and counterpart were measured by an atomic force microscope (AFM). The friction and wear behavior of the coating was studied by a ball-on-disk tribometer. The cross-sections of the wear tracks at different wear stages were observed with a scanning electron microscope (SEM). Maximum contact stresses with different coating thicknesses were calculated by the three-dimensional semi-analytical method (SAM). The results demonstrated that when taking surface asperities into consideration, maximum shear stress at the bonding interface and adjacent substrate showed a dramatic increase during wear and should be responsible for the initiation and propagation of the cracks observed at the final stage of sliding.

Commentary by Dr. Valentin Fuster
J. Tribol. 2013;135(2):021402-021402-11. doi:10.1115/1.4023207.

Preclinical wear evaluation is extremely important in hip replacements, wear being one of the main causes of failure. Experimental tests are attractive but highly cost demanding; thus predictive models have been proposed in the literature, mainly based on finite element simulations. In such simulations, the effect of friction is usually disregarded, as it is considered not to affect the contact pressure distribution. However, a frictional contact could also result in a shift of the location of the nominal contact area, which can thus modify the wear maps. The aim of this study is to investigate this effect in wear prediction for metal-on-metal implants. Wear assessment was based on a purpose-developed mathematical model, extension of a previous one proposed by the same authors for metal-on-plastic implants. The innovative aspect of the present study consists in the implementation of a modified location of the nominal contact point due to friction, which takes advantage of the analytical formulation of the wear model. Simulations were carried out aimed at comparing total and resurfacing hip replacements under several gait conditions. The results highlighted that the adoption of a frictional contact yields lower linear wear rates and wider worn areas, while for the adopted friction coefficient ($f=0.2$), the total wear volume remains almost unchanged. The comparison between total and resurfacing replacements showed higher scaled wear volumes (wear volume divided by wear factor) for the latter, in agreement with the literature. The effect of the boundary conditions (in vivo versus in vitro) was also investigated remarking their influence on implant wear and the need to apply more physiological-like conditions in hip simulators. In conclusion although friction is usually neglected in numerical wear predictions, as it does not affect markedly the contact pressure distribution, its effect in the location of the theoretical contact point was observed to influence wear maps. This achievement could be useful for increasing the correlation between numerical and experimental simulations, usually based on the total wear volume. In order to improve the model reliability, future studies will be devoted to implement the geometry update by combining the present model to finite element analyses. On the other hand, further experimental investigations are required to get out from the wide dispersion of wear factors reported in the literature.

Commentary by Dr. Valentin Fuster

### Research Papers: Elastohydrodynamic Lubrication

J. Tribol. 2012;135(2):021501-021501-8. doi:10.1115/1.4023084.

The dynamic behaviors of a single elastohydrodynamic lubricated (EHL) contact between a rolling element and raceways under wider load and speed ranges are analyzed numerically based on the transient EHL model and the free vibration model. The discrete convolution and fast fourier transform method is implemented in order to increase the computational efficiency associated with elastic deformations and the semisystem approach is applied to improve the solution convergence under severe conditions. The change of mutual approach is selected as the standard of bearing vibrations and the inlet length and dimensionless natural frequency corresponding to the working load and speed are determined. The numerical results demonstrate that the stiffness increases with the increasing load and decreases with speed. However, the changes of damping are different in various working conditions, especially under heavier load and higher speed. It is also indicated that the stiffness and damping increase with the increase in ambient viscosity and the decrease in pressure-viscosity coefficient.

Commentary by Dr. Valentin Fuster
J. Tribol. 2013;135(2):021502-021502-13. doi:10.1115/1.4023085.

Under elastohydrodynamic (EHL) conditions, the temperature in the contact zone determines the load resistance of the lubricant film. Therefore, an efficient assessment of the scuffing risks requires accurate contact temperature predictions. The work presented in this paper develops a simple and prompt temperature model exploitable within any thermal EHL modeling approach. The model incorporates a multilayer lubricant film representation for the heat equation solution. The developments also include an original heat repartition factor expression and a simple formula for handling intermediate values of the Peclet number. The rolling/sliding conditions in the inlet cause temperature rises that also affect the film resistance in the contact zone; this study proposes an inlet temperature rise equation. This formula offers temperature predictions in agreement with reference values. The contact zone temperature predictions for the line of contact problem under sliding/rolling conditions agree remarkably well with published numerical results; for the evaluations presenting the higher absolute difference, the correspondence remained over 94% and 95% for the maximum and mean temperatures, respectively. Both line and elliptical contact conditions were tested and compared to experimental data available in the literature. The analysis evidenced the precision of the estimates; thus attesting to the accuracy of the model under any contact conditions. Finally, the results indicate that, depending on the pressure and speed combination, the shearing zone may occupy around 30% of the film thickness.

Commentary by Dr. Valentin Fuster
J. Tribol. 2013;135(2):021503-021503-11. doi:10.1115/1.4023206.

It has been recognized that purposefully designed surface texturing can contribute to the improvement of tribological performance of elements and friction reduction. However, its optimal parameters may depend on the operating conditions. This paper investigated the effect of a triangle-shaped dimples array on the tribological performance of the lubricated point-contacts under different lubrication regimes, based on the rotational sliding experiment of a patterned steel disk against smooth steel balls. The dimples arrays were produced by laser process and characterized by the 3D profilometer. A series of tests were conducted with different dimple parameters including depth, coverage ratio, size, and direction. Stribecklike curves were obtained to depict the transition of lubrication regimes, and the electrical contact resistance was utilized to qualitatively characterize the lubrication status. The test results showed that the dimples arrays with different sizes, depths and coverage ratios had a distinct effect on the friction behaviors. Compared with the nontextured surfaces, when the dimple depth decreased from 30μm to zero with fixed coverage ratio and size, the friction coefficient firstly decreased, and then increased. The friction coefficient finally approached that of the nontextured surface, during which the lowest value appeared at the dimple depth of approximately 10∼15μm. The coverage ratio of texture showed the similar effect on the friction coefficient. Usually, the coverage ratio of approximately 10% resulted in the lowest friction coefficient. The dimple size and direction also had obvious effects on the friction coefficient. Thus, it can be concluded that there exists a set of optimal values for the dimple depth, coverage ratio, size, and direction to realize the friction reduction.

Commentary by Dr. Valentin Fuster
J. Tribol. 2013;135(2):021504-021504-9. doi:10.1115/1.4023148.

Surface temperature measurements within sliding contacts are useful since interfacial heat dissipation is closely linked to tribological behavior. One of the most powerful techniques for such measurements is in-contact temperature mapping whereby a sliding contact is located beneath an infrared microscope. In this approach, one of the specimens must be transparent to infrared and coated such that radiation components can be distinguished and isolated from background values. Despite its effectiveness, a number of practical constraints prevent this technique from being applied to rough surfaces—a research area where temperature maps could provide much needed two-dimension input data to inform mixed and boundary friction models. The research described in this paper is aimed at improving the infrared temperature mapping technique in terms of validity, robustness, and spatial resolution, so that measurements of rough surfaces contacts can be made. First, Planck's law is applied in order to validate the use of surface coating as a means of removing background radiation. Second, a refined method of calibration is put forward and tested, which negates the need for a soft aluminum coating and hence enables rough surfaces to be measured. Finally, the use of super-resolution algorithms is assessed in order extend spatial resolution beyond the current limit of 6 $μm$.

Commentary by Dr. Valentin Fuster

### Research Papers: Hydrodynamic Lubrication

J. Tribol. 2012;135(2):021701-021701-9. doi:10.1115/1.4007807.

A comprehensive analysis method is proposed to resolve the problem of simulating the complex thermoflow with two kinds of distinct characteristic lengths in a dry gas seal. A conjugated simulation of the complicated heat transfer and the gas film flow is carried out by using the commercial computational fluid dynamics (CFD) software CFX. By using the proposed method, three-dimensional velocity and pressure fields in the gas film flow and the temperature distribution within the sealing rings are investigated for three kinds of film thickness, respectively. A comparison of thermohydrodynamic characteristics of the dry gas seal is conducted between the sealed gas of air and helium. The latter one is used in a helium compressor for a high-temperature gas-cooled reactor (HTGR). From comparisons and discussions of a series of simulation results, it will be found that the comprehensive proposal is effective and simulation results are reasonable. Even under a hypothetical accidental condition, the maximum temperature rise in the dry gas seal is within the acceptable range of HTGR safety requirements.

Commentary by Dr. Valentin Fuster
J. Tribol. 2012;135(2):021702-021702-11. doi:10.1115/1.4023086.

Part of the heat generated by the shearing of the lubricating film during operation of a hydrodynamic bearing is transferred to the bearing components. In the case of the pad, which is usually fully submerged in the lubricating oil, heat is further transferred at the pad free walls to the oil by convection. This mechanism causes a thermal gradient in a pad and, consequently, its thermal deflection. In large hydrodynamic thrust bearings, thermal deflection of the pads is an important phenomenon influencing bearing performance. For such bearings, pad distortion can reach the level of hydrodynamic film thickness and can significantly change the bearing's properties. In this paper, the study of the influence of the heat convection coefficient on the predicted performance of a large hydrodynamic thrust bearing is presented. Two sets of convection coefficients at the pad free surfaces are investigated with the use of thermo-elasto-hydrodynamic (TEHD) calculations. An analysis is carried out for the Itaipu hydro turbine thrust bearing with the outer diameter equal to 5.2 m, which is one of the biggest hydro power plants in the world. The results of the theoretical predictions are compared to the measured data collected during bearing operation.

Commentary by Dr. Valentin Fuster

### Technical Briefs

J. Tribol. 2012;135(2):024501-024501-5. doi:10.1115/1.4023083.

This paper presents a hybrid mobility solution approach to the analysis of dynamically loaded misaligned journal bearings. Mobility data obtained for misaligned bearings (calculated from a finite element representation of the Reynolds equation) are compared with existing curve-fitted mobility maps representative of a perfectly aligned bearing. A relative error analysis of mobility magnitude and direction provides a set of misaligned journal bearing configurations (midplane eccentricity ratio and normalized misalignment angle), where existing curve-fitted mobility map components based on aligned bearings can be used to calculate the resulting journal motion. For bearing configurations where these mobility maps are not applicable, the numerical simulation process proceeds using a complete finite element solution of the Reynolds equation. A numerical example representing a misaligned main bearing in a four-stroke automotive engine illustrates the hybrid solution method. Substantial savings in computational time are obtained using the hybrid approach over the complete finite element solution method without loss of computational accuracy.

Commentary by Dr. Valentin Fuster