In this paper, a Heaviside operator method for solving the special solution of partial differential equations with nonhomogeneous constant coefficients is given. When the nonhomogeneous term is a polynomial, the solution is concise and straightforward. When the nonhomogeneous term is a quasi-polynomial, it can be transformed into a polynomial case by exponential transformation, so as to obtain a special solution.