A hybrid finite element-plane wave expansion method is proposed. Phononic crystal plates in contact with acoustic half-spaces are considered. Complex Bloch wavenumbers are obtained from a nonlinear eigenvalue problem. Attenuation band diagrams are shown for 1D and 2D phononic crystal plates.