This study aims to deeply understand the effect of contact stress and slip ratio on wear performances of bainitic rail steels. The results showed that the wear loss increased as the contact stress and slip ratio increased. Based on the surface damage morphology and microstructural analyses, it revealed that the rolling contact fatigue wear mechanism played a significant role under the low slip ratio, but the dominant wear mechanism transferred to the abrasive wear at the high slip ratio. Meanwhile, the bainitic steel specifically presented worse wear resistance under the abrasive wear mode. Compared with the influence of a slip ratio, the increase in contact stress led to severer plastic flows and contributed to the propagation of cracks. In addition, the contact stress and slip ratio had the opposite effect on the friction coefficient, that is, the friction coefficient of bainitic steels behaved the inverse proportion with the contact stress, but positive proportion with the slip ratio. At last, the increase in slip ratio had more significant effect on the reduction of retained austenite (RA) than the enlargement of contact stress due to the fact that the RA would probably be removed before the martensitic transformation occurred under the abrasive wear mechanism.

The effects of cyclic loading on the surface microstructure evolution of different contact locations in a used pearlitic rail were studied. Microstructures were analyzed using Scanning Electron Microscopy (SEM). Meanwhile, grain boundaries and crystallographic orientations were explored via Electron Backscatter Diffraction (EBSD). At last, wheel–rail contact probabilities and forces were calculated using rail profiles. The results indicate that the side wear region located in the gauge face was 71.5% in the high-angle grain boundaries (HAGBs) fraction, 0.88 in the Kernel Average Misorientation (KAM) value, 36% in the recrystallization (REX) fraction, and had a predominant orientation in grains. The rolling contact fatigue (RCF) region situated at the gauge corner was 66.3% in the HAGBs fraction, 0.92 in the KAM value, 33% in the REX fraction, and was mis-orientated in grains. The region located at the edge of the running band was 60.7% in the low-angle grain boundaries (LAGBs) fraction, 0.97 in the KAM value, 12% in the REX fraction, and was mis-orientated in grains. Continuous dynamic recrystallization (cDRX) took place in wear and RCF regions during the cyclic rolling contact loading, creating ultra-fine grains with a transformation from LAGBs to HAGBs, lower KAM values, and more REX. Grains oriented along [111] parallel to the vertical direction in the wear region were influenced by the dominant normal force, while grains in the RCF region were non-oriented, which was attributed to large lateral and vertical forces of similar magnitudes.

A novel methodology is proposed for accelerating the calculation efficiency of the boundary element modeling of rolling contact. This methodology involves the implementation of an initial solution estimate. The method is to provide the initial estimate value by means of simplifying the method, which is used for the iterative calculation of the boundary element method to solve the normal and tangential contact problems. In the normal contact problem, the normal pressure and contact patch are provided as the initial values for the iterative calculation of the boundary element method. In the tangential contact problem, the initial values for tangential stress and adhesion-slip distribution are provided. A first novel aspect is breaking the initial iteration setting of the traditional BEM, which can significantly reduce the iterations. A second novelty is that a method for determining the potential contact area is proposed to ensure the correctness of boundary element modeling without increasing the computational cost. In the following section, an analysis is conducted to ascertain the impact of the initial solution estimate method on the efficacy of boundary element modeling. The result is that efficiency increased by 69.1% of normal contact calculations, with the initial contact patch exerting the most significant influence. The efficiency of the process under investigation increased by 56.9%, and the most pronounced effect is the distribution of adhesion-slip.