Forced vibration of two-dimensional functionally graded beams is studied in this paper by a developed
meshfree boundary-domain integral equation method. Material properties of the functionally graded
beams are assumed varying continuously either in longitudinal or transvers direction following the
exponential function. The boundary-domain integral equations are derived by using the elastostatic
fundamental solutions based on the two-dimensional elastic theory. Radial integral method (RIM) is
employed to transform the domain integrals into boundary integrals. A meshfree scheme is achieved
through assuming the displacements and accelerations in the domain integrals by a combination of the
radial basis function and polynomials with time dependent coefficients. Wilson-θ, Houbolt as well as two
kinds of damped Newmark's algorithms are applied to accomplish the time integration. The forced vibration
of the functionally graded beam subjected by the harmonic loading and transient loading are
investigated in detail. Numerical examples demonstrate that the four mentioned time integral schemes
are all adapted well to the developed meshfree boundary-domain integral equation method for analyzing
the forced vibration of homogeneous structures. For the analysis of FG structures, it is shown that the
damped Newmark's algorithm can achieve more stable and accurate results