A Computational Study of Interfacial Stress Distribution in Unidirectional Composites and Its Use in Prediction of Brittle Failure
A computational investigation has been conducted to examine the effect of key microstructural and material parameters, namely the minimum inter-fiber spacing (δ) and the fiber/matrix stiffness ratio (Ef/Em), on the interfacial stress distributions in unidirectional composites subjected to transverse uniaxial tensile load. Representative Volume Elements (RVE’s) containing 144 fibers are constructed using a Monte-Carlo (MC) algorithm, imitating random composite structures. The boundary element method is then used to solve the multi-region elasticity problem on these microstructures. We pay particular attention to the statistics of the distribution of the maximum interfacial stresses computed on each fiber; these are found to follow a Weibull distribution, whose specific shape depends on both microstructural and material parameters. Following the weakest-link theory of Batdorf and Crose [1], we derive a statistical formula for the prediction of brittle failure of composite structures caused by interfacial failure. The limitations and possible extensions to the proposed approach are discussed as well.