Homogenization estimates for the average behavior and field fluctuations in cubic and hexagonal viscoplastic polycrystals

The “second-order” homogenization procedure (J. Mech. Phys. Solids 50 (2002) 737) is used to compute estimates of the self-consistent type for the effective response of cubic and hexagonal viscoplastic polycrystals with isotropic textures. The method, which requires the computation of the averages of the stress field and the covariances of its fluctuations over the various grain orientations in an optimally selected “linear comparison polycrystal,” is also used to generate information on the heterogeneity of the stress and strain-rate fields within the polycrystals. In contrast with earlier estimates of the self-consistent type, such as those arising from the “incremental” and “tangent” schemes, the new estimates for the effective behavior are found to satisfy all known bounds, even in the strongly nonlinear, rate-insensitive limit. In addition, they are found to satisfy a recently proposed scaling law at large grain anisotropy. The fluctuations of the stresses and strain rates, which are nonzero for all grain orientations, are found to generally increase with decreasing strain-rate sensitivity (i.e., increasing nonlinearity) and with increasing grain anisotropy (which is typically higher for lower-symmetry systems).     

Second-order theory for the effective behavior and field fluctuations in viscoplastic polycrystals

A recently developed “second-order” homogenization procedure (Ponte Castañeda (J. Mech. Phys. Solids 50 (2002a, b) 737, 759)) is extended to viscoplastic polycrystals and applied to compute the effective response of a certain special class of isotropic polycrystals. The method itself reduces to a simple expression requiring the computation of the averages of the stress field and the covariances of its fluctuations over the various grain orientations in an optimally selected “linear comparison polycrystal”. Therefore, the method not only allows the determination of the effective behavior of the polycrystal, but as a byproduct also yields information on the heterogeneity of the stress and strain-rate fields within the polycrystal. An application is given for a model 2-dimensional, isotropic polycrystal with power-law behavior for the constituent grains. The resulting predictions for the effective behavior are found to satisfy sharp bounds available from the literature and to be consistent with the results of recent numerical simulations. The associated averages and fluctuations of the stresses and strain rates are found to depend strongly on the strain-rate sensitivity (i.e., nonlinearity) and grain anisotropy. In particular, the stress and strain-rate fluctuations were found to grow and become strongly anisotropic with increasing values of the nonlinearity and grain anisotropy parameters.     

Second-order homogenization estimates for nonlinear composites incorporating field fluctuations: II—applications

In Part I of this work, an improved “second-order” homogenization theory was developed. This new theory makes use of generalized secant moduli that are intermediate between the standard secant and tangent moduli of the nonlinear phases, and that depend not only on the averages, or first-moments of the fields in the phases, but also on the second-moments of the field fluctuations, or phase covariance tensors. In this article, the theory, which is known to be exact to second-order in the heterogeneity contrast, is applied to the special cases of rigidly reinforced and porous materials. These are cases corresponding to infinite contrast where fairly explicit analytical expressions of the Hashin-Shtrikman and self-consistent-type may be generated for nonlinear composites. The results show that the new theory improves on the earlier theory (Ponte Castañeda, J. Mech. Phys. Solids 44 (1996) 827) in at least two ways. First, the new theory satisfies rigorous bounds, even near the percolation limit, where field fluctuations become important, and the earlier second-order theory had been found to fail. Second, the new theory predicts fully compressible behavior for porous materials with an incompressible matrix phase, where the earlier theory had also been found to fail. In addition, the new estimates are found to be in better agreement with numerical simulations available from the literature.     

Infinite-contrast periodic composites with strongly nonlinear behavior: Effective-medium theory versus full-field simulations

This paper presents a combined numerical-theoretical study of the macroscopic behavior and local field distributions in a special class of two-dimensional periodic composites with viscoplastic phases. The emphasis is on strongly nonlinear materials containing pores or rigid inclusions. Full-field numerical simulations are carried out using a fast Fourier transform algorithm [Moulinec, H., Suquet, P., 1994. A fast numerical method for computing the linear and nonlinear properties of composites. C. R. Acad. Sci. Paris II 318, 1417–1423.], while the theoretical results are obtained by means of the ‘second-order’ nonlinear homogenization method [Ponte Castañeda, P., 2002. Second-order homogenization estimates for nonlinear composites incorporating field fluctuations. I. Theory. J. Mech. Phys. Solids 50, 737–757]. The effect of nonlinearity and inclusion concentration is investigated in the context of power-law (with strain-rate sensitivity m) behavior for the matrix phase under in-plane shear loadings. Overall, the ‘second-order’ estimates are found to be in good agreement with the numerical simulations, with the best agreement for the rigidly reinforced materials. For the porous systems, as the nonlinearity increases (m decreases), the strain field is found to localize along shear bands passing through the voids (the strain fluctuations becoming unbounded) and the effective stress exhibits a singular behavior in the dilute limit. More specifically, for small porosities and fixed nonlinearity $m>0$, the effective stress decreases linearly with increasing porosity. However, for ideally plastic behavior $(m=0)$, the dependence on porosity becomes non-analytic. On the other hand, for rigidly-reinforced composites, the strain field adopts a tile pattern with bounded strain fluctuations, and no singular behavior is observed (to leading order) in the dilute limit.     

Localization of elastic deformation in strongly anisotropic, porous, linear materials with periodic microstructures: Exact solutions and dilute expansions

Exact solutions are derived for the problem of a two-dimensional, infinitely anisotropic, linear-elastic medium containing a periodic lattice of voids. The matrix material possesses either one infinitely soft, or one infinitely hard loading direction, which induces localized (singular) field configurations. The effective elastic moduli are computed as functions of the porosity in each case. Their dilute expansions feature half-integer powers of the porosity, which can be correlated to the localized field patterns. Statistical characterizations of the fields, such as their first moments and their histograms are provided, with particular emphasis on the singularities of the latter. The behavior of the system near the void close-packing fraction is also investigated. The results of this work shed light on corresponding results for strongly non-linear porous media, which have been obtained recently by means of the “second-order” homogenization method, and where the dilute estimates also exhibit fractional powers of the porosity.     

Second-order homogenization estimates for nonlinear composites incorporating field fluctuations: I—theory

This paper is concerned with the development of an improved second-order homogenization method incorporating field fluctuations for nonlinear composite materials. The idea is to combine the desirable features of two different, earlier methods making use of “linear comparison composites”, the properties of which are chosen optimally from suitably designed variational principles. The first method (Ponte Castañeda, J. Mech. Phys. Solids 39 (1991) 45) makes use of the “secant” moduli of the phases, evaluated at the second moments of the strain field over the phases, and delivers bounds, but these bounds are only exact to first-order in the heterogeneity contrast. The second method (Ponte Castañeda, J. Mech. Phys. Solids 44 (1996) 827) makes use of the “tangent” moduli, evaluated at the phase averages (or first moments) of the strain field, and yields estimates that are exact to second-order in the contrast, but that can violate the bounds in some special cases. These special cases turn out to correspond to situations, such as percolation phenomena, where field fluctuations, which are captured less accurately by the second-order method than by the bounds, become important. The new method delivers estimates that are exact to second-order in the contrast, making use of generalized secant moduli incorporating both first- and second-moment information, in such a way that the bounds are never violated. Some simple applications of the new theory are given in Part II of this work.     

On the overall behavior, microstructure evolution, and macroscopic stability in reinforced rubbers at large deformations: I—Theory

This work presents an analytical framework for determining the overall constitutive response of elastomers that are reinforced by rigid or compliant fibers, and are subjected to finite deformations. The framework accounts for the evolution of the underlying microstructure, including particle rotation, which results from the finite changes in geometry that are induced by the applied loading. In turn, the evolution of the microstructure can have a significant geometric softening (or hardening) effect on the overall response, leading to the possible development of macroscopic instabilities through loss of strong ellipticity of the homogenized incremental moduli. The theory is based on a recently developed “second-order” homogenization method, which makes use of information on both the first and second moments of the fields in a suitably chosen “linear comparison composite,” and generates fairly explicit estimates—linearizing properly—for the large-deformation effective response of the reinforced elastomers. More specific applications of the results developed in this paper will be presented in Part II.     

A general hyperelastic model for incompressible fiber-reinforced elastomers

This work presents a new constitutive model for the effective response of fiber-reinforced elastomers at finite strains. The matrix and fiber phases are assumed to be incompressible, isotropic, hyperelastic solids. Furthermore, the fibers are taken to be perfectly aligned and distributed randomly and isotropically in the transverse plane, leading to overall transversely isotropic behavior for the composite. The model is derived by means of the “second-order” homogenization theory, which makes use of suitably designed variational principles utilizing the idea of a “linear comparison composite.” Compared to other constitutive models that have been proposed thus far for this class of materials, the present model has the distinguishing feature that it allows consideration of behaviors for the constituent phases that are more general than Neo-Hookean, while still being able to account directly for the shape, orientation, and distribution of the fibers. In addition, the proposed model has the merit that it recovers a known exact solution for the special case of incompressible Neo-Hookean phases, as well as some other known exact solutions for more general constituents under special loading conditions.     

Modeling the macroscopic behavior of two-phase nonlinear composites by infinite-rank laminates

A new approach is proposed for estimating the macroscopic behavior of two-phase nonlinear composites with random, particulate microstructures. The central idea is to model composites by sequentially laminated constructions of infinite rank whose macroscopic behavior can be determined exactly. The resulting estimates incorporate microstructural information up to the two-point correlation functions, and require the solution to a Hamilton–Jacobi equation with the inclusion concentration and the macroscopic fields playing the role of ‘time’ and ‘spatial’ variables, respectively. Because they are realizable, by construction, these estimates are guaranteed to be convex, to satisfy all pertinent bounds, to exhibit no duality gap, and to be exact to second order in the heterogeneity contrast. Sample results are provided for two- and three-dimensional power-law composites, and are compared with other homogenization estimates, as well as with numerical simulations available from the literature. The estimates are found to give physically sensible predictions for all the cases considered, even for extreme values of the nonlinearity and heterogeneity contrast. Interestingly, in the case of isotropic porous materials under hydrostatic loadings, the estimates agree exactly with standard Gurson-type models for viscoplastic porous media.     

Numerical simulation of elasto-plastic deformation of composites: evolution of stress microfields and implications for homogenization models

The deformation of a composite made up of a random and homogeneous dispersion of elastic spheres in an elasto-plastic matrix was simulated by the finite element analysis of three-dimensional multiparticle cubic cells with periodic boundary conditions. “Exact” results (to a few percent) in tension and shear were determined by averaging 12 stress-strain curves obtained from cells containing 30 spheres, and they were compared with the predictions of secant homogenization models. In addition, the numerical simulations supplied detailed information of the stress microfields, which was used to ascertain the accuracy and the limitations of the homogenization models to include the nonlinear deformation of the matrix. It was found that secant approximations based on the volume-averaged second-order moment of the matrix stress tensor, combined with a highly accurate linear homogenization model, provided excellent predictions of the composite response when the matrix strain hardening rate was high. This was not the case, however, in composites which exhibited marked plastic strain localization in the matrix. The analysis of the evolution of the matrix stresses revealed that better predictions of the composite behavior can be obtained with new homogenization models which capture the essential differences in the stress carried by the elastic and plastic regions in the matrix at the onset of plastic deformation.     

Homogenization-based constitutive models for porous elastomers and implications for macroscopic instabilities: I—Analysis

The purpose of this paper is to provide homogenization-based constitutive models for the overall, finite-deformation response of isotropic porous rubbers with random microstructures. The proposed model is generated by means of the “second-order” homogenization method, which makes use of suitably designed variational principles utilizing the idea of a “linear comparison composite.” The constitutive model takes into account the evolution of the size, shape, orientation, and distribution of the underlying pores in the material, resulting from the finite changes in geometry that are induced by the applied loading. This point is key, as the evolution of the microstructure provides geometric softening/stiffening mechanisms that can have a very significant effect on the overall behavior and stability of porous rubbers. In this work, explicit results are generated for porous elastomers with isotropic, (in)compressible, strongly elliptic matrix phases. In spite of the strong ellipticity of the matrix phases, the derived constitutive model may lose strong ellipticity, indicating the possible development of shear/compaction band-type instabilities. The general model developed in this paper will be applied in Part II of this work to a special, but representative, class of isotropic porous elastomers with the objective of exploring the complex interplay between geometric and constitutive softening/stiffening in these materials.     

Micromechanics of particle-reinforced elasto-viscoplastic composites: Finite element simulations versus affine homogenization

The micromechanics of elasto-viscoplastic composites made up of a random and homogeneous dispersion of spherical inclusions in a continuous matrix was studied with two methods. The first one is an affine homogenization approach, which transforms the local constitutive laws into fictitious linear thermo-elastic relations in the Laplace–Carson domain so that corresponding homogenization schemes can apply, and the temporal response is computed after a numerical inversion of Laplace transform. The second method is the direct numerical simulation by finite elements of a three-dimensional representative volume element of the composite microstructure. The numerical simulations carried out over different realizations of the composite microstructure showed very little scatter and thus provided – for the first time – “exact” results in the elasto-viscoplastic regime that can be used as benchmarks to check the accuracy of other models. Overall, the predictions of the affine homogenization model were excellent, regardless of the volume fraction of spheres, of the loading paths (shear, uniaxial tension and biaxial tension as well as monotonic and cyclic deformation), particularly at low strain rates. It was found, however, that the accuracy decreased systematically as the strain rate increased. The detailed information of the stress and strain microfields given by the finite element simulations was used to analyze the source of this difference, so that better homogenization methods can be developed.     

Field fluctuations and macroscopic properties for nonlinear composites

A recently introduced nonlinear homogenization method [J. Mech. Phys. Solids 50 ( 2002) 737–757] is used to estimate the effective behavior and the associated strain and stress fluctuations in two-phase, power-law composites with aligned-fiber microstructures, subjected to anti-plane strain, or in-plane strain loading. Using the Hashin–Shtrikman estimates for the relevant “linear comparison composite,” results are generated for two-phase systems, including fiber-reinforced and fiber-weakened composites. These results, which are known to be exact to second-order in the heterogeneity contrast, are found to satisfy all known bounds. Explicit analytical expressions are obtained for the special case of rigid-ideally plastic composites, including results for arbitrary contrast and fiber concentration. The effective properties, as well as the phase averages and fluctuations predicted for these strongly nonlinear composites appear to be consistent with deformation mechanisms involving shear bands. More specifically, for the case where the fibers are stronger than the matrix, the predictions appear to be consistent with the shear bands tending to avoid the fibers, while the opposite would be true for the case where the fibers are weaker.     

A new variational estimate for the effective response of hyperelastic composites

A variational method for predicting the effective properties of hyperelastic composites in terms of available estimates for “hyperelastic comparison composites” is proposed. In some cases this estimate can produce a lower bound on the effective energy-density function. This nonlinear-comparison variational procedure is specialized to classes of fiber and statistically isotropic composites with the aid of appropriate choices of comparison composites with neo-Hookean phases. The end results are given in terms of closed-form expressions for the effective strain energy-density functions, from which the stress-strain relations can be extracted analytically. Explicit analytical estimates for the overall responses of composites whose phases behaviors are governed by the Gent model are obtained. The results for the fiber composites are compared with corresponding finite element simulations of periodic fiber composites as well as with other available estimates. A fine agreement between the predictions obtained via the various estimates is revealed even in the limit of infinite contrast between the properties of the phases.     

Finite deformation of incompressible fiber-reinforced elastomers: A computational micromechanics approach

The in-plane finite deformation of incompressible fiber-reinforced elastomers was studied using computational micromechanics. Composite microstructure was made up of a random and homogeneous dispersion of aligned rigid fibers within a hyperelastic matrix. Different matrices (Neo-Hookean and Gent), fibers (monodisperse or polydisperse, circular or elliptical section) and reinforcement volume fractions (10–40%) were analyzed through the finite element simulation of a representative volume element of the microstructure. A successive remeshing strategy was employed when necessary to reach the large deformation regime in which the evolution of the microstructure influences the effective properties. The simulations provided for the first time “quasi-exact” results of the in-plane finite deformation for this class of composites, which were used to assess the accuracy of the available homogenization estimates for incompressible hyperelastic composites.     

Homogenization-based constitutive models for magnetorheological elastomers at finite strain

This paper proposes a new homogenization framework for magnetoelastic composites accounting for the effect of magnetic dipole interactions, as well as finite strains. In addition, it provides an application for magnetorheological elastomers via a “partial decoupling” approximation splitting the magnetoelastic energy into a purely mechanical component, together with a magnetostatic component evaluated in the deformed configuration of the composite, as estimated by means of the purely mechanical solution of the problem. It is argued that the resulting constitutive model for the material, which can account for the initial volume fraction, average shape, orientation and distribution of the magnetically anisotropic, non-spherical particles, should be quite accurate at least for perfectly aligned magnetic and mechanical loadings. The theory predicts the existence of certain “extra” stresses—arising in the composite beyond the purely mechanical and magnetic (Maxwell) stresses—which can be directly linked to deformation-induced changes in the microstructure. For the special case of isotropic distributions of magnetically isotropic, spherical particles, the extra stresses are due to changes in the particle two-point distribution function with the deformation, and are of order volume fraction squared, while the corresponding extra stresses for the case of aligned, ellipsoidal particles can be of order volume fraction, when changes are induced by the deformation in the orientation of the particles. The theory is capable of handling the strongly nonlinear effects associated with finite strains and magnetic saturation of the particles at sufficiently high deformations and magnetic fields, respectively.     

Onset of macroscopic instabilities in fiber-reinforced elastomers at finite strain

In this paper, we investigate theoretically the possible development of instabilities in fiber-reinforced elastomers (and other soft materials) when they are subjected to finite-strain loading conditions. We focus on the physically relevant class of “macroscopic” instabilities, i.e., instabilities with wavelengths that are much larger than the characteristic size of the underlying microstructure. To this end, we make use of recently developed homogenization estimates, together with a fundamental result of Geymonat, Müller and Triantafyllidis linking the development of these instabilities to the loss of strong ellipticity of the homogenized constitutive relations. For the important class of material systems with very stiff fibers and random microstructures, we derive a closed-form formula for the critical macroscopic deformation at which instabilities may develop under general loading conditions, and we show that this critical deformation is quite sensitive to the loading orientation relative to the fiber direction. The result is also confronted with classical estimates (including those of Rosen) for laminates, which have commonly been used as two-dimensional (2-D) approximations for actual fiber-reinforced composites. We find that while predictions based on laminate models are qualitatively correct for certain loadings, they can be significantly off for other more general 3-D loadings. Finally, we provide a parametric analysis of the effects of the matrix and fiber properties and of the fiber volume fraction on the onset of instabilities for various loading conditions.     

A fully automated numerical tool for a comprehensive validation of homogenization models and its application to spherical particles reinforced composites

This paper presents a fully automated numerical tool for computing the accurate effective properties of two-phase linearly elastic composites reinforced by randomly distributed spherical particles. Virtual microstructures were randomly generated by an algorithm based on molecular dynamics. Composites effective properties were computed using a technique based on Fast Fourier Transforms (FFT). The predictions of the numerical tool were compared to those of analytical homogenization models for a broad range of phases mechanical properties contrasts and spheres volume fractions. It is found that none of the tested analytical models provides accurate estimates for the whole range of contrasts and volume fractions tested. Furthermore, no analytical homogenization models stands out of the others as being more accurate for the investigated range of volume fractions and contrasts. The new fully automated tool provides a unique means for computing, once and for all, the accurate properties of composites over a broad range of microstructures. In due course, the database generated with this tool might replace analytical homogenization models.