In this study, we developed polynomial shape functions for finite element analysis using poly pyramid elements. A polyhedral domain is divided into poly pyramid elements, and polynomial shape functions are defined based on the tetrahedral subdomains within a poly pyramid element. The present poly pyramid finite elements pass patch tests exactly, and numerical experiments show good performance of the poly pyramid finite elements.