This paper presents a novel algorithm for solving meshed DC railway networks with nonlinear load characteristics. The work addresses voltage constraints and power flow analysis in complex traction systems, providing an efficient solution for railway electrification infrastructure analysis and optimization.