Abstract: The piecewise interpolating polynomial is employed to approximate arbitrary dynamic loads in the Duhamel integration for the solution of dynamic response of multiple-degree of-freedom systems. The relavant formulas are derived. Because the Duhamel integration of piecewise polynomial is exact, the proposed solution is more accurate and computationally effort saving than the traditional numerical integration schemes.