skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on March 1, 2027

Title: Correlation function-dependent bounds and Mori–Tanaka effective stiffness estimates: A numerical validation study on biphase transversely isotropic composites
This paper investigates the validity of two different analytical homogenization methods: the Mori–Tanaka mean-field theory and Milton’s correlation function-dependent bounds. We focus on biphase linearly elastic transversely isotropic composites. The composites consist of a matrix reinforced with long fibers of either circular or irregular cross section shapes formed by overlapping circles, with different degrees of radius polydispersity. The Mori–Tanaka effective stiffness depends on the phase moduli, volume fractions, and on a few geometric descriptors of the fibers that can be readily evaluated. In contrast, the computation of Milton’s bounds requires finer knowledge of the microstructure, in terms of two and three-point spatial correlation functions, which are not always analytically tractable. We thus consider very specific random microstructure geometries with known correlation functions. The effective moduli estimates of the two methods are validated against the results of numerical homogenization using the finite element method. It is shown that the Mori– Tanaka predictions of the effective transverse bulk modulus are significantly more accurate than those of the transverse and axial shear moduli. In addition, the predictions of the scheme generally deteriorate with an increasing fiber volume fraction. By contrast, the average of Milton’s upper and lower bounds provides a highly accurate estimate for all three independent effective moduli, without any limitation on the fiber concentration. This study highlights the indisputable effect of the spatial correlation functions on the effective properties of composites, and aspires to pave the way towards the development of more predictive, correlation function-dependent homogenization methods.  more » « less
Award ID(s):
2135584 2416344
PAR ID:
10657144
Author(s) / Creator(s):
;
Publisher / Repository:
Elsevier
Date Published:
Journal Name:
Mechanics of Materials
Volume:
214
Issue:
C
ISSN:
0167-6636
Page Range / eLocation ID:
105580
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. This study aims to detect in which microstructure conditions the Mori–Tanaka scheme is inappropriate to calculate the effective stiffness of a two-phase matrix-inclusion system.We analyze the discrepancy between numerical and Mori–Tanaka stiffness estimates in two-dimensional (2D) solids with crack-like flat cavities. The maximum transfer entropy that occurs between a microstructure feature and a stiffness component discrepancy cannot only detect the phase change between a Mori–Tanaka-like cracked solid to a non-Mori–Tanaka-like cracked solid, but also reveal at which load step that phase change first occurs and which microstructure features most affect that phase change. Further analysis with a binary classifier based on a support vector machine (SVM) algorithm shows that the systematic calculation of nine microstructure features based on six statistical crack network descriptors at each step of a loading path can inform the detection of a microstructure transition. The microstructure features identified here could thus be used to trigger the transition from one homogenization scheme to another during incremental stiffness updates, for example, during the simulation of a load path. 
    more » « less
  2. In the present work, we analyze the applicability of two-step homogenization applied to 3D woven composites with high crimp reinforcement. The available micromechanical homogenization approaches (Hashin, Chamis, Hashin-Shtrikman bounds etc.) were developed and validated for unidirectional composites. These formulas have also been used by the community to homogenize tows in 2D and 3D woven composites including reinforcement architectures with high crimp ratios. However, a rigorous study of their applicability to high-crimp geometries is yet to be performed. We utilize Finite Element Analysis (FEA) to calculate the overall engineering constants (Young’s moduli and shear moduli) of tows having various crimp (𝐶𝑅) and wavelength-to-fiber diameter (𝜆/𝑑) ratios. For this analysis, periodic sinusoidal unit cells following shapes of individual fibers are used. Fiber volume fraction is set to 70% and is the same in all cases. Transversely isotropic carbon fiber and isotropic epoxy matrix are used. The results are compared with overall responses of tows modeled using homogenized tow properties obtained from micromechanics and FEA as well as explicitly modeled tows containing multiple parallel fibers. The results of our analysis show dependence of the overall elastic properties on both crimp ratio and the normalized wavelength. Separation of fiber/tow scales is achieved at 𝜆/𝑑 = 50. 
    more » « less
  3. Short carbon fiber-reinforced polymer composites are widely used in polymer extrusion additive manufacturing (AM), including large-area additive manufacturing (LAAM), due to their enhanced mechanical properties as compared to neat polymers. However, the mechanical properties of these composites depend on microstructural characteristics, including fibers and micro-voids, which are determined during processing. In this work, the correlation between fibers and micro-voids within the microstructure of LAAM polymer composites throughout various processing stages of short carbon fiber-reinforced acrylonitrile butadiene styrene (SCF/ABS) is investigated. The processing stages considered here include the incoming pellets, a single freely extruded strand, a single regularly deposited bead, and a single regularly deposited bead pressed by a mechanical roller. A high-resolution X-ray micro-computed tomography (µCT) system is employed to characterize the microstructural features in terms of the fibers (volume fraction, fiber orientation tensor) and micro-voids (volume fraction, sphericity) in the SCF/ABS samples. The results indicate that micro-voids exist within the microstructure of the SCF/ABS composite in all four stages considered here and that the micro-void volume fraction and micro-void sphericity vary among the test samples. Moreover, the results show a considerable variation in fiber orientation and fiber volume fraction within the microstructure throughout all the stages considered; however, all the samples show the highest alignment in the extrusion/print direction. Furthermore, a correlation is identified between the fiber orientation and the micro-void volume fraction within samples from all four stages considered here. This finding suggests that fibers tend to align more in the extrusion/print direction in regions with less micro-void content. 
    more » « less
  4. The transverse strength of fiber-reinforced composites is a matrix-dominated property whose accurate prediction iscrucial to designing and optimizing efficient, lightweight structures. State-of-the-art analytical models for compositestrength predictions do not account for fiber distribution, orientation, and curing-induced residual stress that greatlyinfluence damage initiation and failure propagation at the microscale. This work presents a novel methodology to develop an analytical solution for transverse composite strength based on computational micromechanics that enables the modeling of stress concentration due to representative volume elements (RVE) morphology and residual stress. Finiteelement simulations are used to model statistical samples of composite microstructures, generate stress-strain curves,and correlate statistical descriptors of the microscale to stress concentration factors to predict transverse strength as a function of fiber volume fraction. Tensile tests of thin plies validated this approach for carbon- and glass-reinforced composites showing promise to obtain a generalized analytical model for transverse composite strength prediction. 
    more » « less
  5. null (Ed.)
    The effect of residual stress build-up on the transverse properties of thermoset composites is studied through direct and inverse process modeling approaches. Progressive damage analysis is implemented to characterize composite stiffness and strength of cured composites microstructures. A size effect study is proposed to define the appropriate dimensions of Representative Volume Elements (RVEs). A comparison between periodic (PBCs) and flat (FBCs) boundary conditions during curing is performed on converged RVEs to establish computationally efficient methodologies. Transverse properties are analyzed as a function of the fiber packing through the nearest fiber distance statistical descriptor. A reasonable mechanical equivalence is achieved for RVEs consisting of 40 fibers. It has been found that process-induced residual stresses and fiber packing significantly contribute to the scatter in composites transverse strength. Variation of ±5% in average strength and 18% in standard deviation are observed with respect to ideally cured RVEs that neglect residual stresses. It is established that process modeling is needed to optimize the residual stress state and improve composite performance. 
    more » « less