The main objective of the study is to reduce the dwell time of wagons at stations and to improve the efficiency of shunting locomotive utilization. The problem has a combinatorial nature, since an increase in the number of loading and unloading fronts leads to a sharp growth in the number of feasible service variants. During the research, a mathematical model describing the servicing process of industrial sidings was developed. The study addressed the problem of determining the optimal sequence of wagon deliveries and the optimal distribution of workload among shunting locomotives. Under conditions where two or more shunting locomotives are used, an optimization method based on the indicator of wagon-hours reduction (σ) was proposed for allocating loading and unloading fronts. Using combinatorial properties, it was shown that many possible allocation variants are symmetric, which allowed the development of a mathematical solution that simplifies the search for an optimal solution. Computational results demonstrated that, at the hypothetical railway station “N-1”, applying the optimal service sequence reduces wagon dwell time by 21% compared with an arbitrary sequence. At the hypothetical station “N-2”, distributing wagon groups between two shunting locomotives improves the efficiency of the servicing process by 26% compared with using a single locomotive. Based on the proposed mathematical model and algorithm, a practical software tool was developed that enables the automatic determination of service sequences for loading and unloading fronts. The software allows the identification of optimal servicing orders, analysis of alternative variants, and evaluation of the efficiency of shunting locomotive utilization.