The numerical simulation of particulate reinforced composites by using a two-dimensional VCFEM formulated with plastic, thermal, and creep strain

The working environment of industrial fields where the Particulate Reinforced Composites (PRCs) are commonly used is complex. In these harsh environments, PRCs produce some nonlinear strains, such as plastic, thermal, and creep strains, which affect the performance of PRCs to some extent. In addition, the actual structure of PRCs contains a massive amount of inclusions, and the number of inclusions could reach the magnitude of millions or even tens of millions. For the large-scale numerical simulation of PRCs, ordinary displacement finite element requires a large number of meshes, which increases their computational complexity and unsolvability. For real PRCs with random particle distribution, homogenisation theory does not accurately capture the true state of stress concentration in PRCs. The Voronoi cell finite element method (VCFEM) can be a good solution to the drawbacks of the above two numerical simulation methods. The subject of the present work explores a modified complementary energy generalized function that takes the constitutive relationship of nonlinear strain into account, proposing a two-dimensional VCFEM, formulated with plastic, thermal, and creep strain. The simulation results are compared with the computed results of a commercial finite element software Marc to verify the validity and accuracy of the VCFEM. Results of the comparison show that the VCFEM has the advantages of simple meshing and faster computation under the same precision. The VCFEM, considering plastic, thermal, and creep strains proposed in this paper, lays the foundation for future large-scale numerical simulations of PRCs containing a large quantity of inclusion structures and can realistically capture the stress concentration state of PRCs. At the end of the article, the impact of the order of stress function on the calculation results is discussed.