Next Article in Journal
Compressive and Thermal Properties of Non-Structural Lightweight Concrete Containing Industrial Byproduct Aggregates
Previous Article in Journal
Numerical Simulation Study of Gas-Solid Heat Transfer and Decomposition Processes of Limestone Calcined with Blast Furnace Gas in a Parallel Flow Regenerative Lime Kiln
Previous Article in Special Issue
Experimental Investigation on the Influence of Fiber Path Curvature on the Mechanical Properties of Composites
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Accuracy Comparison of Micromechanics Models of Particulate Composites against Microstructure-Free Finite Element Modeling

Department of Mechanical Engineering, University of Manitoba, Winnipeg, MB R3T 2N2, Canada
Materials 2022, 15(11), 4021; https://doi.org/10.3390/ma15114021
Submission received: 28 April 2022 / Revised: 31 May 2022 / Accepted: 3 June 2022 / Published: 6 June 2022
(This article belongs to the Special Issue Experimental and Numerical Investigation of Composite Materials)

Abstract

:
Micromechanics models of composite materials are preferred in the analysis and design of composites for their high computational efficiency. However, the accuracy of the micromechanics models varies widely, depending on the volume fraction of inclusions and the contrast of phase properties, which have not been thoroughly studied, primarily due to the lack of complete and representative experimental data. The recently developed microstructure-free finite element modeling (MF-FEM) is based on the fact that, for a particulate-reinforced composite, if the characteristic size of the inclusions is much smaller than the composite representative volume element (RVE), the elastic properties of the RVE are independent of inclusion shape and size. MF-FEM has a number of advantages over the conventional microstructure-based finite element modeling. MF-FEM predictions have good to excellent agreement with the reported experiment results. In this study, predictions produced by MF-FEM are used in replace of experimental data to compare the accuracy of selected micromechanics models of particulate composites. The results indicate that, only if both the contrasts in phase Young’s moduli and phase Poisson’s ratios are small, the micromechanics models are able to produce accurate predictions. In other cases, they are more or less inaccurate. This study may serve as a guide for the appropriate use of the micromechanics models.

1. Introduction

The analytical or semi-analytical solutions developed from micromechanics models of composite materials provide an efficient method for the prediction of effective elastic properties. Compared with the experimental methods and the computer modeling approaches, micromechanics-based analytical solutions are preferred by material engineers in the analysis and design of new composites. The design process is usually an iterative process and often requires a large number of elasticity characterizations for the intermediate designs, which would be very costly and time-consuming if the experimental methods or computer modeling were to be used. Most micromechanics models have been developed for particulate or short-fiber reinforced composites [1,2,3,4,5,6], primarily due to their wide applications in engineering structures and industrial products. These micromechanics models are also important for the description of mechanical behavior of functionally graded materials and structures [7,8]. However, the accuracy of these micromechanics models varies widely, and it is affected by factors such as the volume fraction of inclusions and the contrast of phase properties. Information regarding the accuracy variation of the micromechanics models is necessary for material engineers to appropriately apply them in composite analysis and design. However, a thorough study on the accuracy variation of micromechanics models has not yet been conducted, primarily due to the lack of complete and representative experimental data. Although mechanical testing is the ideal method for evaluating the accuracy of micromechanics models, mechanical testing is expensive and time-consuming. Experimental data reported in the literature do not cover the entire range of inclusion volume-fraction and are unable to represent the various contrasts of phase properties, e.g., [9,10,11] among others. Furthermore, experimental data are often adversely affected by factors such as experimental setup, defects in test specimens, and random measurement errors.
Conventional microstructure-based finite element modeling (MB-FEM) is a more efficient approach than the experimental methods, but it has a number of deficiencies. For example, it is time-consuming to generate a workable geometric model for the composite microstructure; it is challenging to achieve a high volume-fraction of inclusions; and the quality of the finite element mesh is affected by the geometric shape of inclusions and the minimum distance between neighboring inclusions. The recently developed microstructure-free finite element modeling (MF-FEM) [12] may eliminate these deficiencies. It has been demonstrated that the elastic properties predicted by the MF-FEM have good to excellent agreement with the reported experimental data [12].
In this paper, predictions of MF-FEM are used in replacement of experimental data to study how the accuracy of micromechanics models varies with the volume fraction of inclusions and the contrasts of phase properties. The objective of this study is to provide a guide for the appropriate use of the micromechanics models.

2. Materials and Methods

Elastic properties of representative two-phase composites are predicted by the MF-FEM, which are used as references in the comparison of selected micromechanics models. The four particulate composites described in Table 1 are used for the study. They represent different natural or engineering composite materials, with typical combinations of contrasts in the phase Young’s moduli and phase Poisson’s ratios. Generally, a phase material that has larger Young’s modulus would have smaller Poisson’s ratio, and vice versa. The contrast of a phase property is measured by the ratio of the larger value to the smaller one.
The elasticity of a composite material and its constituents are often described by four elastic constants, i.e., Young’s modulus, shear modulus, bulk modulus, and Poisson’s ratio. These properties are not independent from each other but are related via the general elasticity relations. If two of them are known, the rest can be determined from the elasticity relations. For the above reason, the MF-FEM and the selected micromechanics models are required to predict at least two of the four elastic properties.
For consistency, the symbols listed in Table 2 are used in the description of elastic properties of composites and their phases.

2.1. Microstructure-Free Finite Element Modeling (MF-FEM) of Composite Representative Volume Element (RVE)

Microstructure-free finite element modeling (MF-FEM) is a recently developed approach for the study of particulate composites [12]. The fundamental difference between the MF-FEM and the conventional microstructure-based finite element modeling (MB-FEM) is in the way that the inclusions are represented in the finite element model. In the MB-FEM, a geometric model must be constructed to describe the composite microstructure that consists of a matrix and inclusions, which involves a great amount of work but it is not necessary. Based on the continuum mechanics of composites [13,14,15], for a particulate composite, if (1) the inclusions have small aspect ratio, (2) their characteristic size is much smaller than the representative volume element (RVE), and (3) the inclusions have a statistically uniform random distribution in the RVE, the elastic properties of the RVE are independent of the shape and size of the inclusions. The above assumptions are also adopted in almost all the micromechanics models of particulate composites. Previous studies [5,16,17,18] indicate that if the maximum characteristic size of inclusions is smaller than 0.04 times of the RVE dimensions, the RVE elastic properties are almost independent of inclusion shape and size. Now that the shape and size of small inclusions have no effect, the MF-FEM use elements to represent the inclusions. In the version of MF-FEM implemented previously [12], first a uniform mesh of brick elements is generated for the RVE that contains no inclusion; then, a number of the elements are randomly selected and assigned with the inclusion properties, and the rest of the elements are set to have the matrix properties. There are different methods to implement the random selection of elements to represent inclusions. In this study, the MATLAB function p = randperm(N) is used to generate a random permutation of the integers from 1 to N without repeating numbers, where N is the total number of elements in the RVE finite element mesh; the first n elements are selected as the inclusions, the total volume of the selected n elements is required equal to the volume of inclusions, the latter is calculated from the desired volume-fraction of inclusions. No periodicity is required in the distribution of inclusions. One sample of RVE constructed in such a way is displayed in Figure 1.
The RVE is in a cubic shape and the dimensionless side length is L = 100 , and the element size is 2. Therefore, the size ratio between the elements and the RVE is 1:50, which is smaller than the threshold, such that the shape and size of inclusions have no effect on the elastic properties of the RVE. The task of the MF-FEM is to characterize the effective Young’s modulus and effective Poisson’s ratio of the RVE in the three axial directions displayed in Figure 1.
Theoretically, if the assumptions described in the above are satisfied, the composite RVE should be isotropic in its elastic properties. However, the actual distribution of inclusions may not be strictly random and uniform. Therefore, the RVE properties may still be slightly anisotropic [12]. To eliminate the residual anisotropy, the properties in the x , y , and z directions are averaged. The boundary conditions used in the characterization of RVE properties are listed in Table 3.
Based on the mean-field homogenization theory [19,20], the effective Young’s moduli ( E ¯ i ) and effective Poisson’s ratios ( ν ¯ i j ) of the RVE in the three axial directions are determined from the average stresses ( σ ¯ i ) and average strains ( ε ¯ i ), i.e.,
E ¯ i = σ ¯ i ε ¯ i   ,   ( i = x , y , z )
ν ¯ i j = ε ¯ j ε ¯ i   ,   ( i , j = x , y , z )
The average stresses and strains are calculated from the finite element stresses ( σ i ) and strains ( ε i ) via
σ ¯ i = 1 V V σ i   d V ,       ε ¯ i = 1 V V ε i   d V   ( i = x , y , z )
where V is the volume of the RVE.
The effective Young’s modulus and effective Poisson’s ratio of the RVE are obtained as
E ¯ = E ¯ x + E ¯ y + E ¯ z 3
ν ¯ = ν ¯ x y + ν ¯ y x + ν ¯ y z + ν ¯ z y + ν ¯ z x + ν ¯ x z 6

2.2. The Selected Micromechanics Models of Particulate Composites

A large number of micromechanics models have been developed for the prediction of composite elastic properties [3,5,6] but not all of them will be compared in this study, due to the capability of the MF-FEM. The selected micromechanics models should satisfy the following criteria:
  • The micromechanics model was explicitly developed for particulate or short-fiber reinforced composites, where the composites can be considered a homogeneous and isotropic material at the length scale of the RVE.
  • The micromechanics model produces explicit analytical solutions for at least two of the four elastic properties.
  • The analytical solutions do not require an empirical coefficient.
Based on the above criteria, seven micromechanics models were selected. They are the Voigt and Reuss bounds, the Hashin–Shtrikman bounds, the Voigt–Reuss–Hill average, the Mori–Tanaka method, the generalized self-consistent scheme, the isotropized Voigt–Reuss model, and the product of exponential functions. To be referred later in this study, the analytical solutions produced by the selected micromechanics models are presented below with the symbols listed in Table 2.
 (1)
The Voigt and Reuss (VR) bounds [21,22]
The Voigt and Reuss bounds are also referred to as the rule of mixtures and the inverse rule of mixtures in the literature when they are applied to predict the elastic properties of unidirectional fiber-reinforced composites. They are included in this comparison study due to their capability of providing the upper and lower bounds for the elastic properties of particulate composites. The Voigt and Reuss (VR) formulas can be generically presented as
P ¯ V = f 1 P 1 + f 2 P 2
P ¯ R = P 1   P 2 f 1 P 2 + f 2 P 1
where, P and P ¯ represent a generic elastic property, which can be any one of the four elasticity constants; the subscripts V and R indicate Voigt and Reuss models, respectively. The formulas are often used to estimate bounds for Young’s modulus, sometimes also for shear modulus, Poisson’s ratio, and bulk modulus.
 (2)
The Hashin-Shtrikman (HS) bounds [23]
The Hashin-Shtrikman (HS) bound formulas in Equation (7a–d) are derived from the variational principle of minimum potential energy with the concept of stress polarization, which have a stricter mathematical base than the VR bounds.
K ¯ L = K 1 + α 1 ,   α 1 = f 2 ( K 2 K 1 ) ( 3 K 1 + 4   G 1 ) 3 f 2 K 1 + 3   f 1 K 2 + 4 G 1
K ¯ U = K 2 α 2 ,   α 2 = f 1 ( K 2 K 1 ) ( 3 K 2 + 4 G 2 ) 3 f 1 K 2 + 3 f 2 K 1 + 4 G 2
G ¯ L = G 1 + β 1 ,   β 1 = 5 f 2 ( G 2 G 1 ) G 1 ( 3 K 1 + 4 G 1 ) ( 12 f 2 + 8 ) ( G 1 ) 2 + ( 12 f 1 G 2 + ( 6 f 2 + 9 ) K 1 ) G 1 + 6 G 2 K 1 f 1
G ¯ U = G 2 β 2 ,   β 2 = 5 f 1 ( G 2 G 1 ) G 2 ( 3 K 2 + 4 G 2 ) ( 12 f 1 + 8 ) ( G 2 ) 2 + ( 12 f 2 G 1 + ( 6 f 1 + 9 ) K 2 ) G 2 + 6 G 1 K 2 f 2
The subscripts, L and U , represent the lower and the upper bound of effective bulk modulus ( K ¯ ) or effective shear modulus ( G ¯ ).
 (3)
The Voigt–Reuss–Hill (VRH) average [24]
The Voigt–Reuss–Hill average in Equation (8) provides a simple way to estimate the elastic constants of particulate composites.
P ¯ = P ¯ V + P ¯ R 2
where P ¯ V and P ¯ R are determined by Equation (6a,b).
 (4)
The Mori–Tanaka (MT) model [25]
The Mori–Tanaka model is based on Eshelby’s elasticity theory for inhomogeneity in an infinite medium [20,26]. The average internal stress in the matrix is calculated with the transformation of the strain. Interaction among inclusions and the presence of a free boundary are considered in the average elastic energy. The MT formulas for the determination of effective shear modulus ( G ¯ ) and effective bulk modulus ( K ¯ ) are presented in Equation (9a,b).
G ¯ = G 1 + f 2 ( G 1 G 1 ) 1 + f 1 ( G 2 G 1 ) G 1 + G 1 ( 9 K 1 + 8 G 1 ) 6 ( K 1 + 2 G 1 )
K ¯ = K 1 + f 2 ( K 2 K 1 ) 1 + f 1 ( K 2 K 1 ) K 1 + 4 3 G 1
 (5)
The Generalized self-consistent (GSC) model [3]
In the generalized self-consistent model, a spherical inclusion is embedded in a concentric spherical annulus of the matrix material, which is in turn embedded in an infinite medium of unknown effective properties. This three-phase model is an improved version of the original tow-phase self-consistent model [27], where the spherical inclusion is directly embedded in the infinite medium of unknown effective properties. Although the concept is simple, the obtained analytical solutions are complex. The effective bulk modulus in Equation (10b) has the same expression as that in Equation (9b) of the MT model, but the solution of effective shear modulus is much more complex, as displayed in Equation (10a).
G ¯ = ( B + B 2 A C A ) G 1
K ¯ = K 1 + f 2 ( K 2 K 1 ) 1 + f 1 ( K 2 K 1 ) K 1 + 4 3 G 1
The three coefficients A , B , and C in Equation (10a) are expressed by the phase shear moduli ( G 1 , G 2 ), phase Poisson’s ratios ( ν 1 , ν 2 ), phase volume fractions ( f 1 , f 2 ), and another three coefficients ( η 1 , η 2 and η 3 ).
    A = 8 ( G 2 G 1 1 ) ( 4 5 ν 1 ) η 1 f 2 10 3 2 ( 63 ( G 2 G 1 1 ) η 2 + 2 η 1 η 3 ) f 2 7 3 + 252 ( G 2 G 1 1 ) η 2 f 2 5 3 50 ( G 2 G 1 1 ) ( 7 12 ν 1 + 8 ν 1 2 ) η 2 f 2 + 4 ( 7 10 ν 1 ) η 2 η 3 B = 2 ( G 2 G 1 1 ) ( 1 5 ν 1 ) η 1 f 2 10 3 + 2 ( 63 ( G 2 G 1 1 ) η 2 + 2 η 1 η 3 ) f 2 7 3 252 ( G 2 G 1 1 ) η 2 f 2 5 3 + 75 ( G 2 G 1 1 ) ( 3 ν 1 ) η 2 ν 1 f 2 + 3 2 ( 15 ν 1 7 ) η 2 η 3 C = 4 ( G 2 G 1 1 ) ( 5 ν 1 7 ) η 1 f 2 10 3 2 ( 63 ( G 2 G 1 1 ) η 2 + 2 η 1 η 3 ) f 2 7 3 + 252 ( G 2 G 1 1 ) η 2 f 2 5 3 + 25 ( G 2 G 1 1 ) ( ν 1 2 7 ) η 2 f 2 ( 7 + 5 ν 1 ) η 2 η 3
The expressions of η 1 , η 2 , and η 3 are provided below.
η 1 = ( G 2 G 1 1 ) ( 7 10 ν 1 ) ( 7 + 5 ν 2 ) + 105 ( ν 2 ν 1 ) η 2 = ( G 2 G 1 1 ) ( 7 + 5 ν 2 ) + 35 ( 1 ν 2 ) η 3 = ( G 2 G 1 1 ) ( 8 10 ν 1 ) + 15 ( 1 ν 1 )
 (6)
The Isotropized Voigt-Reuss (Iso-VR) model [28]
The Iso-VR model is the result of isotropization of the conventional Voigt and Reuss model, which together represent a transversely isotropic material model. The isotropization assumes that the normal and the shear strain energy in the transversely isotropic model are respectively equivalent to those in the isotropized model. The resulting formula is provided in the equation below.
P ¯ = 2 1 P ¯ V + 1 P ¯ R
The above formula is applicable to effective Young’s modulus and effective shear modulus [28].
 (7)
The product of exponential functions (PEF) [29]
The product of exponential functions in Equation (12) was established from experimental data of Young’s moduli measured from bovine bone specimens. The formula is also applicable to effective shear modulus.
P ¯ = P 1 f 1 P 2 f 2
In the above equation, P 1 and P 2 are originally the Young’s moduli of the organic and the inorganic phase in the bone. Bone is a typical composite material in nature, the formula is thus included in this comparison study.
Two of the four effective properties that are predicted by the MF-FEM and by each of the selected micromechanics models are marked in Table 4. The rest of effective properties are determined via the elasticity relations.
In this comparison study, the predictions of MF-FEM serve as the references in replace of experimental data to compare the accuracy of the selected micromechanics models. For the micromechanics models that predict specific values of the effective properties, the accuracy is measured by the relative error defined in Equation (13).
δ = | P ¯ MF FEM P ¯ MM | P ¯ MF FEM × 100 %
where P ¯ MF FEM and P ¯ MM are, respectively, the effective property predicted by the MF-FEM and one of the selected micromechanics models.
For the micromechanics models that estimate bounds, i.e., VR and HS, their performance is evaluated by how well the MF-FEM predictions are enclosed within the estimated bounds.

3. Results

Figure A1, Figure A2, Figure A3 and Figure A4 in the Appendix A present the effective properties of the composites in Table 1, predicted by the MF-FEM and the selected micromechanics models with the volume fraction of the stiffer varying from 0.0 to 1.0. Based on the above results, relative errors in the effective properties predicted by VRH, GSC, MT, Iso-VR, and PEF are calculated using Equation (13). The obtained relative errors are displayed in Figure A5, Figure A6, Figure A7 and Figure A8 in the Appendix A. The relative errors are inhomogeneous and demonstrate fluctuation over the range of volume fraction; it is inconvenient to use them to measure the overall accuracy. Therefore, averages of the errors, calculated using the following equation, are used to measure the overall performance of the models.
δ ¯ = i n δ i n
where i indicates one of the volume fractions f 2 = { 0.1 ,   0.2 ,   ,   0.9 } and n = 9 .
The average relative errors in the effective Young’s modulus, shear modulus, bulk modulus, and Poisson’s ratio of the four composites (Table 1) are displayed in Figure 2.
The following observations can be made from the results regarding the accuracy of the five models, i.e., VRH, GSC, MT, Iso-VR, and PEF.
(1)
The accuracy of the models is inhomogeneous over the range of volume fraction; see Figure A5, Figure A6, Figure A7 and Figure A8 in the Appendix A. For GSC and MT models, the maximum relative error usually occurs in the second half of the range. For VRH, Iso-VR, and PEF, the accuracy fluctuates over the range.
(2)
Only for Composite #1, which has small contrasts in both its phase Young’s moduli and phase Poisson’s ratios, all the models have reasonable accuracy in all four effective properties. The maximum relative error is below 1%; see Figure A5a, Figure A6a, Figure A7a and Figure A8a, in addition to Figure 2.
(3)
For Composite #2, which has a small contrast of phase Young’s moduli but a large contrast in phase Poisson’s ratios, the models have acceptable accuracy for effective Young’s modulus and effective shear modulus, with the maximum relative error below 6%. But the accuracy for effective bulk modulus and effective Poisson’s ratio is quite poor, with average error above 15%.
(4)
For Composites #3 and #4, both have a large contrast in phase Young’s moduli, only if the volume fraction of the stiffer phase is low (<0.2), the models have acceptable accuracy, cf. Figure A1, Figure A2, Figure A3 and Figure A4. None of the models have acceptable overall accuracy, see Figure 2.
(5)
Generally, PEF has better performance than the other models; see Figure A5, Figure A6, Figure A7 and Figure A8 and Figure 2.
The above observations indicate that contrasts in both phase Young’s moduli and phase Poisson’s ratios affect the accuracy of the models, but the effect from phase Young’s modulus is dominant.
Results related to the VR and HS bounds are presented in Figure 3, Figure 4, Figure 5 and Figure 6. A number of phenomena can be detected from the results.
  • Only if the contrasts of phase Young’s moduli and phase Poisson’s ratios are small, the VR and HS bounds are able to enclose the MF-FEM predictions.
  • The gap between the upper and the lower bound of either HS or VR model is primarily dependent upon the contrast of phase Young’s moduli. If the contrast of phase Young’s moduli is small, the bounds are tight; otherwise, the bounds are loose. The contrast of phase Poisson’s ratios has a much lower significant effect on the gap.
  • Contrary to the observations reported in some of the previous studies, the HS bounds may not be always enclosed by the VR bounds, e.g., the effective Young’s moduli in Figure 3b, the effective bulk moduli in Figure 5b, and the effective Poisson’s ratios in Figure 6b. This phenomenon is related to the large contrast of phase Poisson’s ratios.
  • MF-FEM predictions may be out of both the HS and VR bounds, e.g., the effective Young’s moduli in Figure 3b and the effective Poisson’s ratios in Figure 6c.

4. Discussion

The results of this comparison study indicate that the differences between the predictions of MF-FEM and those of the selected micromechanics models are affected by both the contrasts of phase properties and the volume fractions; the largest differences occur in Composites #3 and #4, which have large contrasts in phase properties, particularly in phase Young’s moduli, and when the volume-fraction of the stiffer phase is high. Although previous experimental studies produced a large volume of data for composites that have small contrasts of phase properties and low volume-fractions of inclusions, e.g., [9,10,11,30,31] among many others, experimental data for composites that have large contrast and high volume-fraction of inclusions are rare in the literature for unknown reasons. Nevertheless, composites that have large phase contrast and high volume-fraction of inclusions indeed exist, bone is a good example. If bone is considered as an organic-inorganic composite [29], the organic and the inorganic phases have material properties of extraordinary disparity; nevertheless, the bone has superior mechanical properties over many engineering materials.
The rationality of using MF-FEM predictions in replacement of experimental data for the evaluation of micromechanics model accuracy is supported by a number of facts. First, MF-FEM predictions have good to excellent agreement with the experimental data produced from composites that have small phase contrast or from composites that have large phase contrast but low volume-fraction of inclusions [12]. Second, the MF-FEM is based on the same fundamental assumptions as the selected micromechanics models, that is, the composite is assumed statistically homogeneous and isotropic, and the material properties of the composite are independent of the shape and size of the inclusions. Theoretically, micromechanics models are more comparable with the MF-FEM than with mechanical testing, considering that the accuracy of both the MF-FEM and micromechanics models is only affected by their assumptions, while mechanical testing data may contain errors from different sources. To make analytical solution possible and simple, various special assumptions are adopted in the selected micromechanics models, which narrow the applicable scope of the models and result in low accuracy when the assumptions are not satisfied. The MF-FEM does not introduce any special assumptions. Therefore, it can be reasonably stated that the MF-FEM may not be able to completely replace mechanical testing, but the MF-FEM is more accurate than the micromechanics models.
The differences between the MF-FEM and the selected micromechanics models are induced mainly by the special assumptions adopted in the models. For example, the assumption of dilute dispersion is explicitly or implicitly adopted in many micromechanics models including GSC and MT. By assuming that an inclusion is embedded in an infinite medium, dilute dispersion disregards the effect of interaction among inclusions. However, the rationality of this assumption disappears when the volume-fraction of inclusions is high. For composites that have small contrast of phase properties, the effect of inclusion interaction is trivial; for composites that have large phase contrast, the effect becomes much more significant. The dilute dispersion assumption is implicitly adopted in GSC and MT [32,33]. The limitations of conventional micromechanics models have been theoretically discussed in the literature [34,35,36,37]. The assumption of dilute dispersion is identified as the fundamental source for the limitations. The effect of dilute dispersion can only be fully revealed when the volume fraction is close to full packing of inclusions [3]. The Mori–Tanaka method is found unacceptably inaccurate when it is applied far beyond the usual dilute range [34]. Both VRH and Iso-VR are based on the iso-strain and iso-stress assumptions, which represent the two extreme scenarios that phase materials work together. Under the iso-strain condition, the phase materials work in parallel to achieve maximum stiffness, while under the iso-stress condition, the phase materials work in series to have maximum flexibility. The actual situation is somewhere between the two extreme scenarios, which may not have been accurately represented by either VRH or Iso-VR. The PEF formula was developed from bone mechanical testing data using statistical methods [29]. Although the formula is sometimes more accurate than the other models in the tested cases, the physical meaning of the formula is not as clear as the others. Another common issue for VRH, Iso-VR, and PEF is that they treat the four elasticity constants as independent parameters, but they are actually dependent upon each other.

5. Conclusions

In this comparison study, the accuracy of the selected micromechanics models, including some of the most commonly used ones, is evaluated against the MF-FEM. The study results indicate that, only if the particulate composite has small contrasts in both of its phase Young’s moduli and phase Poisson’s ratios, the selected models have reasonable accuracy. When applied to composites that have large contrasts of phase properties, particularly the large contrast in phase Young’s moduli, the models may have extremely large error when the volume-fraction of the stiffer phase is high. The main limitation of the current study is that it is purely numerical based; an experimental study with the emerging 3D-printing technique would demonstrate the prospective application of the MF-FEM in engineering.

Funding

This research was funded by the Natural Sciences and Engineering Research Council (NSERC) with grant number RGPIN-2019-05372 and the APC was funded by Research Manitoba.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The funds from the Natural Sciences and Engineering Research Council (NSERC) of Canada and the Research Manitoba are gratefully acknowledged.

Conflicts of Interest

The author declares no conflict of interest.

Appendix A

Figure A1, Figure A2, Figure A3 and Figure A4 present predictions of the four effective properties, i.e., Young’s modulus, shear modulus, bulk modulus, and Poisson’s ratio, by the MF-FEM and by the selected micromechanics models for the four representative composites in Table 1, with the volume fraction of the stiffer phase varying from 0.0 to 1.0.
Figure A1. Composite #1—small contrast in both phase Young’s modulus and phase Poisson’s ratio. (a) Young’s modulus; (b) shear modulus; (c) bulk modulus; (d) Poisson’s ratio.
Figure A1. Composite #1—small contrast in both phase Young’s modulus and phase Poisson’s ratio. (a) Young’s modulus; (b) shear modulus; (c) bulk modulus; (d) Poisson’s ratio.
Materials 15 04021 g0a1aMaterials 15 04021 g0a1b
Figure A2. Composite #2—small contrast in phase Young’s modulus and large contrast in phase Poisson’s ratio. (a) Young’s modulus; (b) shear modulus; (c) bulk modulus; (d) Poisson’s ratio.
Figure A2. Composite #2—small contrast in phase Young’s modulus and large contrast in phase Poisson’s ratio. (a) Young’s modulus; (b) shear modulus; (c) bulk modulus; (d) Poisson’s ratio.
Materials 15 04021 g0a2
Figure A3. Composite #3—large contrast in phase Young’s modulus and small contrast in phase Poisson’s ratio. (a) Young’s modulus; (b) shear modulus; (c) bulk modulus; (d) Poisson’s ratio.
Figure A3. Composite #3—large contrast in phase Young’s modulus and small contrast in phase Poisson’s ratio. (a) Young’s modulus; (b) shear modulus; (c) bulk modulus; (d) Poisson’s ratio.
Materials 15 04021 g0a3
Figure A4. Composite #4—large contrast in both phase Young’s modulus and phase Poisson’s ratio. (a) Young’s modulus; (b) shear modulus; (c) bulk modulus; (d) Poisson’s ratio.
Figure A4. Composite #4—large contrast in both phase Young’s modulus and phase Poisson’s ratio. (a) Young’s modulus; (b) shear modulus; (c) bulk modulus; (d) Poisson’s ratio.
Materials 15 04021 g0a4
Figure A5, Figure A6, Figure A7 and Figure A8 show relative errors measured by Equation (13), in the effective Young’s modulus, shear modulus, bulk modulus, and Poisson’s ratio.
Figure A5. Relative error in effective Young’s modulus.
Figure A5. Relative error in effective Young’s modulus.
Materials 15 04021 g0a5
Figure A6. Relative error in effective shear modulus.
Figure A6. Relative error in effective shear modulus.
Materials 15 04021 g0a6aMaterials 15 04021 g0a6b
Figure A7. Relative error in effective bulk modulus. [Note: the pink line is coincident with the blue one].
Figure A7. Relative error in effective bulk modulus. [Note: the pink line is coincident with the blue one].
Materials 15 04021 g0a7
Figure A8. Relative error in effective Poisson’s ratio.
Figure A8. Relative error in effective Poisson’s ratio.
Materials 15 04021 g0a8

References

  1. Huang, H.-B.; Huang, Z.M. Micromechanical prediction of elastic-plastic behavior of a short fiber or particle reinforced composite. Compos. Part A Appl. Sci. Manuf. 2020, 134, 105889. [Google Scholar] [CrossRef]
  2. Jagadeesh, G.; Setti, S.G. A review on micromechanical methods for evaluation of mechanical behavior of particulate reinforced metal matrix composites. J. Mater. Sci. 2020, 55, 9848–9882. [Google Scholar] [CrossRef]
  3. Christensen, R.M. A critical evaluation for a class of micromechanics models. J. Mech. Phys. Solids 1990, 38, 379–404. [Google Scholar] [CrossRef]
  4. Kundalwal, S.I.; Ray, M.C. Effective properties of a novel composite reinforced with short carbon fibers and radially aligned carbon nanotubes. Mech. Mater. 2012, 53, 47–60. [Google Scholar] [CrossRef]
  5. Pindera, M.-J.; Khatam, H.; Drago, A.S.; Bansal, Y. Micromechanics of spatially uniform heterogeneous media: A critical review and emerging approaches. Compos. Part B 2009, 40, 349–378. [Google Scholar] [CrossRef]
  6. Raju, B.; Hiremath, S.R.; Roy Mahapatra, D. A review of micromechanics based models for effective elastic properties of reinforced polymer matrix composites. Compos. Struct. 2018, 204, 607–619. [Google Scholar] [CrossRef]
  7. Singh, R.K.; Rastogi, V. A review on solid state fabrication methods and property characterization of functionally graded materials. Mater. Today: Proc. 2021, 47, 3930–3935. [Google Scholar] [CrossRef]
  8. Barbaros, I.; Yang, Y.; Safaei, B.; Yang, Z.; Qin, Z.; Asmael, M. State-of-the-art review of fabrication, application, and mechanical properties of functionally graded porous nanocomposite materials. Nanotechnol. Rev. 2022, 11, 321–371. [Google Scholar] [CrossRef]
  9. Smith, J.C. Experimental values for the elastic constants of a particulate-filled glassy polymer. J. Res. Natl. Bur. Stand. A Phys. Chem. 1976, 80, 45–49. [Google Scholar] [CrossRef]
  10. Doi, H.; Fujiwara, Y.; Miyake, K.; Oosawa, Y. A systematic investigation of elastic moduli of WC-Co alloys. Metall. Mater. Trans. B 1970, 1, 1417–1425. [Google Scholar] [CrossRef]
  11. Richard, T.G. The mechanical behavior of a solid microsphere filled composite. J. Compos. Mater. 1975, 9, 108–113. [Google Scholar] [CrossRef]
  12. Luo, Y. Microstructure-free finite element modeling for elasticity characterization and design of fine-particulate composites. J. Compos. Sci. 2022, 6, 35. [Google Scholar] [CrossRef]
  13. Dvorak, G. Micromechanics of Composite Materials; Springer: Berlin/Heidelberg, Germany, 2013. [Google Scholar]
  14. Christensen, R.M. Mechanics of Composite Materials; Dover Publications: Mineola, NY, USA, 2012. [Google Scholar]
  15. Sharma, S. Mechanics of Particle- and Fiber-Reinforced Polymer Nanocomposites: From Nanoscale to Continuum Simulations; Wiley: Hoboken, NJ, USA, 2021. [Google Scholar]
  16. Trias, D.; Costa, J.; Turon, A.; Hurtado, J. Determination of the critical size of a statistical representative volume element (SRVE) for carbon reinforced polymers. Acta Mater. 2006, 54, 3471–3484. [Google Scholar] [CrossRef]
  17. Harper, L.T.; Qian, C.; Turner, T.A.; Li, S.; Warrior, N.A. Representative volume elements for discontinuous carbon fibre composites—Part 1, Boundary conditions. Compos. Sci. Technol. 2012, 72, 225–234. [Google Scholar] [CrossRef]
  18. Pelissou, C.; Baccou, J.; Monerie, Y.; Perales, F. Determination of the size of the representative volume element for random quasi-brittle composites. Int. J. Solids Struct. 2009, 46, 2482–2855. [Google Scholar] [CrossRef] [Green Version]
  19. Pierard, O.; Friebel, C.; Doghri, I. Mean-field homogenization of multi-phase thermo-elastic composites: A general framework and its validation. Compos. Sci. Technol. 2004, 64, 1587–1603. [Google Scholar] [CrossRef]
  20. Eshelby, J.D. The determination of the elastic field of an ellipsoidal inclusion, and related problems. Proc. Royal Soc. Lond. Ser. A 1957, 241, 376. [Google Scholar]
  21. Voigt, W. Uber die Beziehung zwischen den beiden Elastizitatskonstanten Isotroper Korper. Wied. Ann. 1889, 38, 573–587. [Google Scholar] [CrossRef] [Green Version]
  22. Reuss, A. Berechnung der Fließgrense von Mischkristallen auf Grund der Plastizitätsbedingung für Einkristalle. Z. Angew. Math. Und Mech. 1929, 9, 49–58. [Google Scholar] [CrossRef]
  23. Hashin, Z.; Shtrikman, S. A variational approach to the theory of the elastic behaviour of multiphase materials. J. Mech. Phys. Solids 1963, 11, 127–140. [Google Scholar] [CrossRef]
  24. Chuang, D.H.; Buessem, W.R. The Voigt-Reuss-Hill approximation and elastic moduli of polycrystalline MgO, CaF2, β-ZnS, ZnSe, and CdTe. J. Appl. Phys. 1967, 38, 2535. [Google Scholar] [CrossRef]
  25. Mori, T.; Tanaka, K. Average stress in matrix and average elastic energy of materials with misfitting inclusions. Acta Metall. 1973, 21, 571–574. [Google Scholar] [CrossRef]
  26. Eshelby, J.D. Elastic inclusions and inhomogeneities. Prog. Solid Mech. 1961, 2, 89–140. [Google Scholar]
  27. Cristensen, R.M.; Lo, K.H. Solutions for effective shear properties in three phase sphere and cylinder models. J. Mech. Phys. Solids 1979, 27, 315–330. [Google Scholar] [CrossRef]
  28. Luo, Y. Isotropized Voigt-Reuss model for prediction of elastic properties of particulate composites. Mech. Adv. Mater. Struct. 2021. [Google Scholar] [CrossRef]
  29. Luo, Y.; Wu, X. Bone quality is dependent on the quantity and quality of organic-inorganic phases. J. Med. Biol. Eng. 2020, 40, 273–281. [Google Scholar] [CrossRef]
  30. Ravichandran, K.S. Elastic properties of two-phase composites. J. Am. Ceram. Soc. 1994, 77, 1178–1184. [Google Scholar] [CrossRef]
  31. Ibrahim, I.A.; Mohamed, F.A.; Lavernia, E.J. Particulate reinforced metal matrix composites—A review. J. Mater. Sci. 1991, 26, 1137–1156. [Google Scholar] [CrossRef]
  32. Hill, R. A self-consistent mechanics of composite materials. J. Mech. Phys. Solids 1965, 13, 213–222. [Google Scholar] [CrossRef]
  33. Benveniste, Y. Revisiting the generalized self-consistent scheme in composites: Clarification of some aspects and a new formulation. J. Mech. Phys. Solids 2008, 56, 2984–3002. [Google Scholar] [CrossRef]
  34. Christensen, R.M.; Schantz, H.G.; Shapiro, J. On the range of validity of the Mori-Tanaka method. J. Mech. Phys. Solids 1992, 40, 69–73. [Google Scholar] [CrossRef]
  35. Ferrari, M. Asymmetry and the high concentration limit of the Mori-Tanaka effective medium theory. Mech. Mater. 1991, 11, 251–256. [Google Scholar] [CrossRef]
  36. Gluzman, S.; Mityushev, V.; Nawalaniec, W. Computational Analysis of Structured Media; Academic Press: Cambridge, MA, USA, 2017. [Google Scholar]
  37. Mityushev, V.; Drygas, P. Effective properties of fibrous composites and cluster convergence. Multiscale Model. Simul. 2019, 17, 696–715. [Google Scholar] [CrossRef]
Figure 1. Representative volume element used in MF-FEM.
Figure 1. Representative volume element used in MF-FEM.
Materials 15 04021 g001
Figure 2. The average relative errors in the effective elastic properties.
Figure 2. The average relative errors in the effective elastic properties.
Materials 15 04021 g002
Figure 3. MF-FEM predictions and VR/HS bounds of effective Young’s modulus.
Figure 3. MF-FEM predictions and VR/HS bounds of effective Young’s modulus.
Materials 15 04021 g003aMaterials 15 04021 g003b
Figure 4. Relation between MF-FEM predictions and VR/HS bounds of effective shear modulus.
Figure 4. Relation between MF-FEM predictions and VR/HS bounds of effective shear modulus.
Materials 15 04021 g004
Figure 5. Relation between MF-FEM predictions and VR/HS bounds of effective bulk modulus.
Figure 5. Relation between MF-FEM predictions and VR/HS bounds of effective bulk modulus.
Materials 15 04021 g005
Figure 6. Relation between MF-FEM predictions and VR/HS bounds of effective Poisson’s ratio.
Figure 6. Relation between MF-FEM predictions and VR/HS bounds of effective Poisson’s ratio.
Materials 15 04021 g006
Table 1. Two-phase particulate composites and their phase properties.
Table 1. Two-phase particulate composites and their phase properties.
Composite #Softer PhaseStiffer PhasePhase Contrast
Young’s Modulus (MPa)Poisson’s RatioYoung’s Modulus (MPa)Poisson’s RatioYoung’s ModulusPoisson’s Ratio
180.00.20120.00.15Small Small
280.00.45120.00.15Small Large
380.00.2012,000.00.15LargeSmall
480.00.4512,000.00.15LargeLarge
Table 2. Symbols of composite and phase properties.
Table 2. Symbols of composite and phase properties.
E i Young’s modulus of phase i E ¯ Effective Young’s modulus of the composite
G i Shear modulus of phase i G ¯ Effective shear modulus of the composite
K i Bulk modulus of phase i K ¯ Effective bulk modulus of the composite
ν i Poisson’s ratio of phase i ν ¯ Effective Poisson’s ratio of the composite
f i Volume fraction of phase i P ¯ ,   P i Generic property of the composite and phase i
Table 3. RVE boundary conditions for the characterization of composite effective properties.
Table 3. RVE boundary conditions for the characterization of composite effective properties.
RVE SurfaceYoung’s Modulus ( E ¯ i ,   i = x , y , z ) and Poisson’s Ratio ( ν ¯ i j ,   i , j = x , y , z )
E ¯ x ,   ν ¯ x y ,   ν ¯ x z E ¯ y ,   ν ¯ y x ,   ν ¯ y z E ¯ z ,   ν ¯ z x ,   ν ¯ z y
x = 0 u x = 0 u x = 0 u x = 0
y = 0 u y = 0 u y = 0 u y = 0
z = 0 u y = 0 u y = 0 u z = 0
x = 100 u x = 1 Homogeneous u x Homogeneous u x
y = 100 Homogeneous u y u y = 1 Homogeneous u y
z = 100 Homogeneous u z Homogeneous u z u z = 1
Table 4. Effective properties that are directly determined by MF-FEM and micromechanics models.
Table 4. Effective properties that are directly determined by MF-FEM and micromechanics models.
Young’s Modulus ( E ¯ ) Shear Modulus ( G ¯ ) Bulk Modulus ( K ¯ ) Poisson’s Ratio ( ν ¯ )
MF-FEMx x
VR boundsxx
HS bounds xx
VRHxx
MT xx
GSC xx
Iso-VRxx
PEFxx
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Luo, Y. An Accuracy Comparison of Micromechanics Models of Particulate Composites against Microstructure-Free Finite Element Modeling. Materials 2022, 15, 4021. https://doi.org/10.3390/ma15114021

AMA Style

Luo Y. An Accuracy Comparison of Micromechanics Models of Particulate Composites against Microstructure-Free Finite Element Modeling. Materials. 2022; 15(11):4021. https://doi.org/10.3390/ma15114021

Chicago/Turabian Style

Luo, Yunhua. 2022. "An Accuracy Comparison of Micromechanics Models of Particulate Composites against Microstructure-Free Finite Element Modeling" Materials 15, no. 11: 4021. https://doi.org/10.3390/ma15114021

APA Style

Luo, Y. (2022). An Accuracy Comparison of Micromechanics Models of Particulate Composites against Microstructure-Free Finite Element Modeling. Materials, 15(11), 4021. https://doi.org/10.3390/ma15114021

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop