As a new type of composite material, particle reinforced composite materials, which has good mechanical properties and secondary machining, have been widely used in mechanical, biological, aerospace, military, motor and other important industrial areas. With the development of science and technology lots of research and numerical simulation have been carried on at home and aboard. Because of the reinforcements, the overall mechanical properties have been significantly improved. At the same time, fracture properties and fatigue characteristics are lower. This paper, based on the VCFEM, lead in the traditional FEM to research particle reinforced composite materials, tending to get a better result in simulating the mechanical property. The principle of the voronoi unit Based on the particle composites as the research object, combine four node isoparametric element with the voronoi cell mesh together to complete the structure calculation.to make the description of the distribution of high stress of the interface of inclusion particles more accurate. As we know that due to the reinforcement, the original features of the material have changed. To a certain extent, reduce its applicability. The interface layer is the important reason of the damage. As Fig1-1, Contact interface between particle and matrix cracking cracks. This makes a third crack boundary. means crack boundary, means inclusion particles, crack boundary, means substrate crack boundary, so we can find that at it still satisfy: (1) Fig 1 Containing inclusions, consider interface debonded voronoi cell At and because of the cracking, the two boundaries have no restraint, so satisfy the following boundary conditions (on border) (2) (on border) (3)