The lubrication behavior and mechanical characteristics of the main bearing area of an aero-engine main shaft bearing determine the reliability and life of the main shaft bearing. In aero-engine main shaft bearings, the lubricant not only plays the role of lubrication but also affects the dynamic characteristics of the bearing; therefore, if the lubricant drag force is insufficient, it will lead to rolling body slipping. Slipping not only affects the reliability of the bearing operation but also will make the temperature of the contact area instantaneously increase, leading to the occurrence of gluing, scraping and other lubrication failure phenomena in the main bearing area. A lubricant under the shear conditions of traction characteristics is actually the external manifestation of rheological properties. Rheological properties are one of the elastic fluid power lubrication theories and are an important part of the study. Elasto-hydrodynamic lubrication theory of the oil film pressure, film thickness and temperature and solid domains interact to form a thermal–fluid–solid coupling relationship; this coupling relationship affects the main bearing area of the lubrication behavior and mechanical properties, thus affecting the lubrication state of the bearings and dynamic characteristics. With the continuous improvement of aero-engine performance requirements for main shaft bearings, it is of great significance to carry out a coupling study of the lubrication behavior and mechanical properties of the bearing contact zone under heavy load, high speed and high temperature conditions to improve the service performance, reliability and life of the bearings.

By optimizing the carburizing heat treatment process, the grain size of the carburized layer of M50NiL steel was successfully refined to the sub-micron level. The mechanism for the generation of a large number of sub-micron crystal regions (SMCR) is that dislocations are entangled and linked due to the pinning effect of nanometer-sized carbides. In this study, a stacking fault energy (SFE) model for austenite in M50NiL steel was established. First-principles calculations were employed to investigate the effects of alloying elements, as well as the position and quantity of carbon (C) atoms, on the generalized stacking fault energy (GSFE). The variations in SFE were further analyzed in combination with differential charge density calculations. The simulation results revealed that the addition of alloying elements excluding nickel led to a reduction in the unstable stacking fault energy. Differential charge density analysis indicated that this decrease was associated with the weakening of Fe–Fe bonds in the L0 layer, where stacking faults occurred. When C atoms are interstitially dissolved near the L0 layer, the Fe–Fe bonds near the L0 layer are enhanced, and the unstable stacking fault energy is correspondingly increased. Compared with the pure iron system, the combined effect of alloying elements and C atoms in M50NiL steel maintained a relatively low level of both the unstable stacking fault energy and the stacking fault formation barrier, provided that C atoms were not dissolved in the L1 layer. This condition was favorable for dislocation slip. Meanwhile, the stable stacking fault energy significantly increased, enhancing the stability of austenite. Based on these simulation results, the relationship between the GSFE of austenite in M50NiL steel and the formation of subgrains and twins within the submicron crystalline regions of the carburized layer was discussed.