In this paper, by transforming the isotropy elastodynamic equations into the first-order hyperbolic ones, we have improved the O(2,4) finite difference time domain method with the staggered grid. Using the approach in this paper, we calculated the band gap in two-dimensional phononic crystals. The calculated and measured results match well.