Interface stress distributions in transversely loaded continuous fiber composites: parallel computation in multi-fiber RVEs using the boundary element method
Interface stress distributions in transversely loaded continuous fiber composites: parallel computation in multi-fiber RVEs using the boundary element method
复制标题
DOI：
10.1016/j.compscitech.2003.07.006
复制
发表时间：
2004-07-01
影响因子：
9.1
通讯作者：
Papathanasiou, TD
中科院分区：
文献类型：
--
作者：
Chen, XM;Papathanasiou, TD
本文介绍了横向加载连续纤维增强复合材料中微观结构对界面应力（空间和统计分布）影响的计算研究结果。我们使用边界元方法（BEM）的并行实现来代替传统的单元模型，使用蒙特卡罗（MC）算法来分析包含多达144个单个光纤横截面的大型代表性体积单元（RVEs），这些截面位于感兴趣的域内。边界元法非常适合这种多包含建模，因为它的相对精度和容易网格复杂的几何形状。通过控制最小允许的纤维间距，MC算法允许生成具有不同程度均匀性（或无序性）的微结构。在这些多光纤rve上计算空间相关函数和总体属性，以确保它们的大小允许推断出适当的统计信息。在验证了计算机程序并给出数值算例的基础上，研究了关键的微观结构和材料参数，即最小纤维间距（delta）、纤维/基体刚度比（E-f/E-m）和纤维体积分数（phi）对应力分布的影响。特别注意的是在每根纤维上计算的最大界面应力分布的统计；它们遵循威布尔分布，其具体形状取决于材料和微观结构参数。这些信息可以很容易地合并到预测复合失效的统计模型中。(C) 2003 Elsevier Ltd.版权所有。
We present the results of a computational investigation on the effect of microstructure on the (spatial and statistical distribution) of interface stresses in transversely loaded continuous fiber-reinforced composites. Instead of traditional unit-cell models, we use a parallel implementation of the boundary element method (BEM) to analyze large representative volume elements (RVEs) containing up to 144 individual fiber cross-sections placed within the domain of interest using a Monte Carlo (MC) algorithm. The BEM is well suited for this multi-inclusion modeling because of its relative accuracy and the ease of meshing complex geometries. By controlling the minimum allowable inter-fiber spacing, the MC algorithm allows for the generation of microstructures showing various degrees of homogeneity (or, disorder). Spatial correlation functions and overall properties are computed on these multi-fiber RVEs to ensure that their size allows for appropriate statistics to be inferred. After validating the computer code and giving some numerical examples, the effects of key microstructural and material parameters, namely the minimum inter-fiber spacing (delta), the fiber/matrix stiffness ratio (E-f/E-m) and the fiber volume fraction (phi) on the stress distributions are examined. Particular attention is placed on the statistics of the distribution of the maximum interface stresses computed on each fiber; these are found to follow a Weibull-like distribution, whose specific shape depends on material and microstructural parameters. This information can be readily incorporated into statistical models for the prediction of composite failure. (C) 2003 Elsevier Ltd. All rights reserved.