This work proposes a novel multi-step dynamic FE compression method to generate Representative Volume Elements (RVEs) of composites containing a variety of inclusions. This method is actualized through a sequential and four-stage procedure: (1) Sparse and periodic inclusions exhibiting a predefined orientation distribution are generated by implementing a modified random sequential adsorption algorithm, (2) Sparse inclusions undergo biaxial compression into the region of the targeted RVE via a dynamic FE analysis, (3) Periodic inclusions are compressed into the region utilizing a dynamic FE analysis with periodic boundary conditions, and (4) Positions and orientations of the compressed inclusions are extracted and the RVE in computer-aided design format is then generated. The proposed method confers advantages from four distinct perspectives: (1) applicability to various inclusion geometries, including ellipses, lobules, polygons, kidneys and stars, (2) ability to generate periodic RVEs of composites with high inclusion volume fractions (up to 80.0% for circular inclusions), (3) simple and straightforward numerical implementation without explicitly considering inclusion intersection check and (4) capacity to predefine inclusion orientation distribution. Statistical analyses utilizing multiple spatial descriptors, i.e., inclusion orientation angle, local volume fraction, Voronoi cell area, nearest-neighbor distance and orientation, second order intensity function, radial distribution function and two-point probability function, confirm randomness of inclusion distribution in the generated RVEs. The elastic properties and damage behaviors of composites via the FE homogenization method are predicted based on the generated RVEs and are compared with those of available experimental data, the literature and the mean-field homogenization models to demonstrate effectiveness of the proposed multi-step dynamic FE compression method.
Read full abstract