Implementation of time-weighted residual method for simulation of flexural waves in multi-span Timoshenko beams subjected to various types of external loads: from stationary loads to accelerating moving masses