Implementation of a constitutive micromechanical model for damage analysis in glass mat reinforced composite structures

A micromechanical model based on a probabilistic approach is implemented in the finite element code CASTEM 2000 to develop numerical simulations that efficiently predict the overall damaged behaviour of random oriented fibre composites. The proposed damage constitutive model is based upon the generalised Mori and Tanaka scheme and Eshelbyʹs equivalence theory. Damage mechanisms occurring at each composite constituent (fibres, matrix and interface) are associated to Weibull probabilistic functions to model their onset and progressive growth at the microscopic scale level. It is obvious that the damaged behaviour of the composite material depends widely on the microscopic material parameters (fibre length, fibre volume fraction, fibre orientation, …). On one hand, the micromechanical model uses homogenisation techniques which enabled us to link these microscopic parameters to the material behaviour and to evaluate explicitly their influences. On the other hand, the implementation of the derived behaviour law into a finite element code enabled us to reflect on the effect of these microscopic parameters on the overall response of a simple composite structure presenting heterogeneous stress fields. In fact, the damage evolution in each constituent (local scale) and the related stiffness reduction are estimated at any material point (integration point) or node of the considered structure subject to a specific loading. Numerical simulations of a composite plate with a hole under in-plane tension were performed to validate the implementation of the behaviour law. Numerical results have been compared to experimental curves and damage evolutions monitored by acoustic emission techniques. Simulations agree well with experimental results in terms of damage onset and growth.