The paper proposes a novel approach to modeling electrified transportation systems. The proposed solution reflects the mechanical dynamics of vehicles as well as the distribution and losses of electric supply. Moreover, energy conversion losses between the mechanical and electrical subsystems and their bilateral influences are included. Such a complete model makes it possible to replicate, e.g., the impact of voltage drops on vehicle acceleration or the necessity of partial disposal of regenerative braking energy due to temporary lack of power transmission capability. The modeling methodology uses a flexible twin data-bus structure, which poses no limitation on the number of vehicles and enables modeling complex traction power supply structures. The proposed solution is suitable for various electrified transportation systems including suburban and urban systems. The modeling methodology is applicable i.a. to Matlab/Simulink, which makes it broadly available and customizable, and provides short computation time. The applicability and accuracy of the method were verified by comparing simulation and measurement results on an exemplary trolleybus system operating in Pilsen, Czech Republic. Simulation of daily operation of an area including four supply sections and maximal simultaneous number of nine vehicles showed a good conformance with the measured data, with the difference in the total consumed energy not exceeding 5%.