Coupled VEM–BEM Approach for Isotropic Damage Modelling in Composite Materials

Numerical prediction of composite damage behaviour at the microscopic level is still a challenging engineering issue for the analysis and design of modern materials. In this work, we document the application of a recently developed numerical technique based on the coupling between the Virtual Element Method (VEM) and the Boundary Element Method (BEM) within the framework of continuum damage mechanics (CDM) to model the in-plane damage evolution characteristics of composite materials. BEM is a widely adopted and efficient numerical technique that reduces the problem dimensionality due to its underlying formulation. It substantially simplifies the pre-processing stage and decreases the computational effort without affecting the solution's accuracy. VEM is a recent generalization to general polygonal mesh elements of the Finite Element Method that ensures noticeable simplification in the data preparation stage of the analysis, notably for computational micro-mechanics problems, whose analysis domain often features complex geometries. The numerical technique has been applied to artificial microstructures, representing the transverse section of composite material with stiffer circular-shaped inclusions embedded in a softer matrix. BEM is used to model the inclusions that are supposed to behave within the linear elastic range, while VEM is used to model the surrounding matrix material, developing nonlinear behaviours. Numerical results are reported and discussed to validate the proposed method.