Meso-scale modeling of plasticity in composites

Abstract A two-phase, pixelized representation of a simple composite microstructure is used as a basis for comparing the accuracy of models that extract local information for use on a larger scale. The local microstructure is homogenized on the meso-scale, using a moving window Generalized Method of Cells (GMC) homogenization technique [1]. Effective elastic constitutive behavior is modeled, and several possible representations of the local yield surface of the composite are compared. The goal of this work is to determine a simple and accurate method to approximate a meso-scale yield surface based on the material microstructure, such that the approximation can be used in a macro-scale analysis. In order to determine the accuracy of the homogenized models, each is compared with an analytical study. The problem of a single elasto-plastic circular inclusion in an elasto-plastic matrix under radial loading at an infinite boundary is considered. A benchmark solution to this problem is developed based on the work of Mendelson [2]. Investigation is made into the effect of the variation of several model parameters. The pixel resolution of the digitized circular fiber image, and the moving window size used within the analysis are both varied. The yield surface of the meso-scale element is also varied, according to two different requirements: first yield in the microstructure, or an averaged yield in the matrix microstructure. According to the definition of yield in the material, the yield surface is then fitted with parameters to match the behavior of the material, and to describe the material on the meso-scale. These yield parameters are first defined as those of Hill’s anisotropic yield criterion [3]. The Hill’s yield criterion is chosen for its simplicity as an anisotropic generalization of the Mises yield criterion. A second set of modified Hill’s parameters is also defined, which can be applied to highly anisotropic materials for which the Hill’s criterion is undefined. A finite element analysis is carried out on each meso-scale model, with varying degrees of homogenization. The internal energy and degree of plasticity in these models are compared against the results of the Mendelson benchmark. It is shown that the Hill’s criterion parameters can be successfully modified. A subcell initial yield criterion is shown to be most effective in a material that exhibits high stress concentrations, while a matrix average yield criterion is shown to be more effective in a material with lower stress concentrations.

[1]  J. Aboudi Mechanics of composite materials - A unified micromechanical approach , 1991 .

[2]  G. J. Dvorak,et al.  Generalized Initial Yield Surfaces for Unidirectional Composites , 1974 .

[3]  Sarah C. Baxter,et al.  Analysis of Heterogeneous Composites Based on Moving-Window Techniques , 2003 .

[4]  D. Salinas,et al.  Initial Yield Surface of a Unidirectionally Reinforced Composite , 1972 .

[5]  M. S. Madhava Rao,et al.  Yielding in Unidirectional Composites Under External Loads and Temperature Changes* , 1973 .

[6]  Sarah C. Baxter,et al.  Simulation of local material properties based on moving-window GMC , 2001 .

[7]  Jacob Aboudi,et al.  Micromechanical analysis of yielding of metal matrix composites , 1988 .

[8]  Forrest J. Masters,et al.  Non-Gaussian simulation of local material properties based on a moving-window technique , 2003 .

[9]  Martin Ostoja-Starzewski,et al.  Random field models of heterogeneous materials , 1998 .

[10]  M. P. Kamat,et al.  Three-Dimensional Inelastic Finite Element Analysis of Laminated Composites , 1981 .

[11]  A. Mendelson Elastoplastic analysis of circular cylindrical inclusion in uniformly stressed infinite homogeneous matrix , 1968 .

[12]  M. Paley,et al.  Micromechanical analysis of composites by the generalized cells model , 1992 .

[13]  Marek-Jerzy Pindera,et al.  Plastic Analysis of Complex Microstructure Composites Using the Generalized Method of Cells , 1999 .

[14]  R. Hill A theory of the yielding and plastic flow of anisotropic metals , 1948, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[15]  Carl T. Herakovich,et al.  Mechanics of Fibrous Composites , 1997 .

[16]  Marek-Jerzy Pindera,et al.  An analytical model for the inelastic axial shear response of unidirectional metal matrix composites , 1997 .

[17]  Sarah C. Baxter,et al.  Characterization of Random Composites Using Moving-Window Technique , 2000 .