AbstractThis paper presents a combined numerical-theoretical study of the macroscopic behavior and local field distributions in a special class of two-dimensional periodic composites with viscoplastic phases. The emphasis is on strongly nonlinear materials containing pores or rigid inclusions. Full-field numerical simulations are carried out using a fast Fourier transform algorithm [Moulinec, H., Suquet, P., 1994. A fast numerical method for computing the linear and nonlinear properties of composites. C. R. Acad. Sci. Paris II 318, 1417–1423.], while the theoretical results are obtained by means of the ‘second-order’ nonlinear homogenization method [Ponte Castañeda, P., 2002. Second-order homogenization estimates for nonlinear composites in...