This paper employs the numerical assembly method (NAM) to determine the βexactβ frequencyβresponse amplitudes of a multiple-span beam carrying a number of various concentrated elements and subjected to a harmonic force, and the exact natural frequencies and mode shapes of the beam for the case of zero harmonic force. First, the coefficient matrices for the intermediate concentrated elements, pinned support, applied force, left-end support and right-end support of a beam are derived. Next, the overall coefficient matrix for the whole vibrating system is obtained using the numerical assembly technique of the conventional finite element method (FEM). Finally, the exact dynamic response amplitude of the forced vibrating system corresponding to each specified exciting frequency of the harmonic force is determined by solving the simultaneous equations associated with the last overall coefficient matrix. The graph of dynamic response amplitudes versus various exciting frequencies gives the frequencyβresponse curve for any point of a multiple-span beam carrying a number of various concentrated elements. For the case of zero harmonic force, the above-mentioned simultaneous equations reduce to an eigenvalue problem so that natural frequencies and mode shapes of the beam can also be obtained.