A new domain-boundary element method is developed for elastodynamic analysis of functionally graded Timoshenko beams. Three governing partial differential equations of motion are derived by considering through-the-thickness variations of the physical properties. Weighted-residual forms are imposed utilizing the static fundamental solutions. These forms are then reduced to three integral equations containing domain integrals with time derivatives of unknown functions. Through domain discretization and shape function approximation, integral equations are converted to a system of ordinary differential equations in time. Forced dynamic response is revealed by solving the system of equations via Houbolt method. Comparison of dynamic responses generated for homogeneous beams to those calculated through an analytical solution and finite difference method verify the developed procedures. Further parametric analyses are performed for functionally graded Timoshenko beams under step, harmonic, and impulsive loadings. The numerical results presented illustrate the influence of material inhomogeneity on time histories of deflection and stress. Domain-boundary element method is demonstrated to be an effective technique for elastodynamic analysis of functionally graded structures. (C) 2017 Elsevier Ltd. All rights reserved.