The paper presents a method of computing electrical and mechanical variables of BLDC motors. It takes into account electrical, magnetic and mechanical phenomena in the power supply-converter-BLDC motor-load machine system. The solution to the problem is the so-called circuit-field method. The results determined with the use of time stepping finite element method were used as the parameters of equations of the developed mathematical model. Losses in the motor, losses in transistors and diodes of the converter as well as the actual back EMF waveforms, variable moment of inertia and variable load torque are accounted for. The designed laboratory stand and the test results are presented in the paper. The experimental verification shows the correctness of the developed method, algorithm and program. The developed computational method is universal with respect to different electromechanical systems with cylindrical BLDC motors. It can be applied to electromechanical systems with BLDC motors operating at constant but also variable load torque and moment of inertia.