Integrated Procedure for Homogenization of Particle Random Com- posites Using Virtual Element Method

Marco Pingaro;Patrizia Trovalusci;Emanuele Reccia

2017-01-01

---

Abstract

Particle-based composites, used in many engineering applications or present in nature, exhibit a mi- crostructure made of randomly distributed inclusions (particles) embedded into a dissimilar matrix. A key aspect, recently investigated by many researchers, is the evaluation of appropriate mechanical properties to be adopted for the study of their behavior. Homogenization procedures may be adopted for the definition of equivalent moduli able to take into account at the macroscale the material prop- erties emerging from the internal microstructure. Respect to the classic homogenization approach, in the case of materials with random microstructure it is not possible to a-priori define a Represen- tative Volume Element (RVE), this being an unknown of the problem. A possible way to solve this problem is to approach the RVE using finite–size scaling of intermediate control volume elements, named Statistical Volume Elements (SVEs), and proceed to homogenization. Here homogenization, consistent with a generalized Hill–Mandel condition, is adopted in conjunction with a statistical procedure, by which scale–dependent bounds on classical moduli are obtained using Dirichlet and Neumann boundary conditions for solving boundary value problems. The outlined procedure has provided significant results, also extended to non–classical continuum formulations [1], but with high computational cost which prevents the possibility to perform series of parametric analyses. The procedure here proposed automates all the steps to perform: from the simulations of each ran- dom realization of the microstructure to the solutions of the boundary value problems for the SVEs, up to the evaluation of the final size of the RVE for the homogenization of the random medium. Moreover, the adoption of an innovative computational method, such as the Virtual Element Method (VEM) [2], allow us to reduce the computational burden. The VEM methodology has many compu- tational advantages such as robust stiffness matrix (can be exactly computed in precision machine) and accuracy versus the number of degrees of freedom. For the numerical analysis we adopt a polygonal mesh for the matrix and a single VEM element for the inclusions. The results obtained by adopting this integrated homogenization procedure with VEM are compared with the results pre- viously obtained, by some of the authors, using a standard Finite Elements procedure taking into account two different types of inclusions, either stiffer or softer than the matrix. Several simulations are then performed by modifying the material contrast (ratio between the moduli of the materials components) deriving the size of the RVE for performing homogenization on various kinds of two– phases random composites.