The contact areas between rail and wheel are critical in determining the performance of rail vehicles. Analyzing the mechanics, materials, tribology, and geometry of the interface is important to understand the contact damage mechanisms, fatigue parameter and fatigue life. Developing a simplified and novel approach to predict crack initiation and fatigue life is key for investigating contact fatigue of rail wheel interaction. This study aims to investigate rolling contact fatigue (RCF) cracks and fatigue life in railroad wheels using Finite Element Model (FEM). Using a proper material model that accounts for plastic shakedown, ratchetting, and cyclic hardening for both the rail and wheel materials, the FEM model is analyzed with the help of Abaqus 2020 tool. The submodelling technique was implemented to reduce computation demand. Stress and strain history outputs extracted from the submodel (at 0 mm, 1.3 mm, 2.6 mm, and 3.6 mm depths) were used to predict critical plane orientation, crack initiation and fatigue life using critical plane models. The orientation of critical (crack) plane was found using a novel MATLAB algorithm by searching for a plane with maximum fatigue parameter (FPmax). The investigation reveals that the location of maximum damage is the subsurface and the crack angles fall within the range of (0°, 40°), (62°, 123°) and (143°, 180°) with a 5 % margin and at a depth of 2.6 mm below the wheel surface. The determined fatigue crack initiation and life at this specific location were 6178 and 7174 cycles, respectively. Predicting the critical zone which is susceptible to shorter crack initiation and fatigue life is a crucial to prevent wheel failure.