In this work, a non‐conforming three‐dimensional finite element coupled with direct methods and homogenization technique is presented for the limit analysis of periodic metal‐matrix composites. Using this element, which is constructed from bilinear shape functions and enriched by internal second‐order polynomials, limit analysis of composite material can be efficiently carried out. Accuracy and overall performance are illustrated through comparison with different structural solid elements in the context of direct as well as incremental methods. It is shown that the limit domain of periodic composites for different fiber distributions and volume fractions provides a foundation for the structural design. (© 2010 Wiley‐VCH Verlag GmbH & Co. KGaA, Weinheim)