Matrix-based Chebyshev spectral approach to dynamic analysis of non-uniform Timoshenko beams