Voronoi cell finite element method (VCFEM) is introduced in this paper to describe the elastic-plastic-creep behavior of particle reinforced composites. The interfacial damage is simulated by partly debonding between Matrix and inclusion. A validation of the nonlinear behavior of the cell element has been carry out by comparing VCFEM results with those calculated by the general finite element package MARC and ABAQUS, and good agreements are found. A microstructure with five inclusions is taken as an example to describe the cyclic stress-strain behavior under different particulate orientation condition, and it shows the influence of the topological microstructure of inclusions. Thermomechanical fatigue properties are also investigated and the loops of stress-strain show the great differences of fatigue behavior between the in-phase case and out-of-case.