This paper presents a numerical analysis method to determine flange contact at variable wheel positions. The shapes of the wheel and rail surface functions with surface parameters. The Newton-Rhapson method for wheel-rail contact can provide fast solutions, but may not yield true values at optimization process with the condition that minimum distance is zero can time-consuming. A compound method, combining the Newton-Rhapson methods the optimization process method is proposed to provide exact solutions efficiently.