Computational Physics Group
Cohesive Modeling of Dewetting in Particulate Composites: Micromechanics vs. Multiscale Finite Element Analysis
H.M Inglis1, P.H. Geubelle2, K. Matous2,3, H. Tan1 and Y. Huang1
1Department of Mechanical Engineering
2Department of Aerospace Engineering
3Center for Simulation of Advanced Rockets
University of Illinois at Urbana-Champaign
Urbana, IL 61801, USA.
The effect of damage due to particle
debonding on the constitutive response of highly filled
composites is investigated using two multiscale homogenization schemes:
one based on a closed-form micromechanics solution, and the other on
the finite element implementation of the mathematical theory of
homogenization. In both cases, the particle debonding process is
modeled using a bilinear cohesive law which relates cohesive tractions
to displacement jumps along the particle matrix interface. A detailed
comparative assessment between the two homogenization schemes is
presented, with emphasis on the effect of volume fraction, particle
size and particle-to-particle interaction.
We have performed a detailed
comparative assessment of micromechanics and finite element based
homogenization schemes for the problem of debonding damage in a 2-D
particulate composite. Special emphasis has been placed on the ability
of the two schemes to capture particle-to-particle interactions and
the effect of dissimilar particle sizes.
The 2-D micromechanics model developed by Tan et al. (2005a) is effective at capturing key features of the macroscopic stress-strain response for minimal computational effort. The shortcomings of this model are that it cannot capture the instability inherent in the system or the heterogeneous stress and strain fields. In a highly filled composite, interactions between particles are a significant contributor to failure through local stress concentrations and the occurrence of localization. The micromechanics model does not capture these interactions during the failure process, and hence ceases to be predictive under high volume fraction, f > 0:5, when stress concentrations begin to play a significant role in the solution, or when the particle distribution is random, resulting in localization. Both models demonstrate that, for large differences in particle diameters, it is unnecessary to model the debonding of smaller particles, but is sufficient to represent their contribution to damage nucleation and to the stiffness of the matrix.
The ability of the MTHFE code to function as a direct numerical simulation for validation of simpler models has been demonstrated. The MTH-based code has an ability to capture a richness of detail about the physical response of the system. The method is capable of solving more complex loading cases, and can be extended to include different material models.
This work was supported by the Center for Simulation of Advanced Rockets (CSAR) under contract number B341494 by the U.S. Department of Energy. K. Matous and P. H. Geubelle also acknowledge support from ATK/Thiokol (Program Managers, J. Thompson and Dr. I. L. Davis). H. Tan and Y. Huang acknowledge additional support from ONR Composites for Marine Structures Program (Grant N00014-01-1-0205, Program Manager Dr. Y.D.S. Rajapakse).
© 2009 Notre Dame and Dr. Karel Matous