A methodology for an accurate evaluation of the linearization procedures in nonlinear mean field homogenization

A systematic methodology for the evaluation of the linearization procedures sustaining mean field homogenization theories for nonlinear composite materials is proposed and applied as an illustration to various recently proposed ‘affine’ and ‘second-order’ formulations for nonlinear elasticity. It relies on the analysis of composites for which both the exact nonlinear homogenization problem and the homogenization problem associated with the ‘linear comparison material’ defined by the linearization procedure can be solved numerically with the same accuracy and for the same microstructure. The comparison of the results then provides a rigorous evaluation of the effects of the sole linearization method. To cite this article: A. Rekik et al., C. R. Mecanique 333 (2005).     

Stationary variational estimates for the effective response and field fluctuations in nonlinear composites

This paper presents a variational method for estimating the effective constitutive response of composite materials with nonlinear constitutive behavior. The method is based on a stationary variational principle for the macroscopic potential in terms of the corresponding potential of a linear comparison composite (LCC) whose properties are the trial fields in the variational principle. When used in combination with estimates for the LCC that are exact to second order in the heterogeneity contrast, the resulting estimates for the nonlinear composite are also guaranteed to be exact to second-order in the contrast. In addition, the new method allows full optimization with respect to the properties of the LCC, leading to estimates that are fully stationary and exhibit no duality gaps. As a result, the effective response and field statistics of the nonlinear composite can be estimated directly from the appropriately optimized linear comparison composite. By way of illustration, the method is applied to a porous, isotropic, power-law material, and the results are found to compare favorably with earlier bounds and estimates. However, the basic ideas of the method are expected to work for broad classes of composites materials, whose effective response can be given appropriate variational representations, including more general elasto-plastic and soft hyperelastic composites and polycrystals.     

Strength homogenization of matrix-inclusion composites using the linear comparison composite approach

A homogenization procedure to estimate the macroscopic strength of nonlinear matrix-inclusion composites with different strength characteristics of the matrix and inclusions, respectively, is presented in this paper. The strength up-scaling is formulated within the framework of the yield design theory and the linear comparison composite (LCC) approach, introduced by Ponte Castaneda (2002) and extended to frictional models by Ortega et al. (2011), which estimates the macroscopic strength of composite materials in terms of an optimally chosen linear thermo–elastic comparison composite with a similar underlying microstructure. In the paper various combinations for the underlying material behavior for the individual phases of the composite are considered: The matrix phase can be a quasi frictional material characterized either by a Drucker–Prager-type (hyperbolic) or an elliptical strength criterion, which predicts a strength limit also in hydrostatic compression, while the inclusion phase either may represent empty pores, pore voids filled with a pore fluid, rigid inclusions, or solid inclusions, whose strength characteristics also maybe described by a Drucker–Prager-type or an elliptical strength criterion. For generating the homogenized strength criterion efficiently in such general cases of matrix-inclusion composites, a novel algorithm is proposed in the paper. The validation of the proposed strength homogenization procedure for selected combinations of strength characteristics of the matrix material and the inclusions is conducted by comparisons with experimental results and alternative existing strength homogenization models.     

基于迭代法的非线性弹性均质化研究

微观结构对复合材料的宏观力学性能具有至关重要的影响, 通过合理设计复合材料微观结构可以得到期望的宏观性能. 均质化方法作为一种有效的设计方法, 它从微观结构的角度出发, 利用均匀化的概念, 实现了对复合材料宏观力学性能的预测和设计. 而当考虑非线性因素, 均质化的实现就非常困难. 本文利用双渐近展开方法, 将位移按照宏观位移和微观位移展开, 推导了非线性弹性均质化方程. 通过直接迭代法, 对非线性弹性均质化方程进行了求解, 并给出了具体的迭代方法和实现步骤. 本文基于迭代步骤和非线性弹性均质化方程编写MATLAB 程序, 对3种典型本构关系的周期性多孔材料平面问题进行了计算, 对比细致模型的应变能、最大位移和等效泊松比, 对程序及迭代方法的准确性进行了验证. 之后对一种三元橡胶基复合材料进行多尺度均质化, 将其分为芯丝尺度和层间尺度. 用线弹性的均质化方法得到了芯丝尺度的等效弹性参数, 并将其作为层间尺度的材料参数. 在层间尺度应用非线性弹性均质化方法对结构进行计算, 得到材料的宏观等效性能, 并以实验结果为基准进行评价.      

Estimation of very narrow bounds to the behavior of nonlinear incompressible elastic composites

Variational bounds for the effective behavior of nonlinear composites are improved by incorporating more-detailed morphological information. Such bounds, which are obtained from the generalized Hashin–Shtrikman variational principles, make use of a reference material with the same microstructure as the nonlinear composite. The geometrical information is contained in the effective properties of the reference material, which are explicitly present in the analytical formulae of the nonlinear bounds. In this paper, the variational approach is combined with estimates for the effective properties of the reference composite via the asymptotic homogenization method (AHM), and applied to a hexagonally periodic fiber-reinforced incompressible nonlinear elastic composite, significantly improving some recent results.     

Determination of the effective elastic–plastic response of metal–ceramic composites

The mechanical response of metal–ceramic composites is analysed through a homogenization model accounting for the mechanical behaviour of the constituent materials. In order to achieve this purpose a nonlinear homogenization method based on the phase field approach has been suitably implemented into a numerical code. A prescribed homogenized strain state is applied to a unit volume element of a metal–ceramic composite with proportional loading in which all components of the strain tensor are proportional to one scalar parameter. The mechanical response of the material has been modeled by considering a von Mises plasticity model for the metal phase and a Drucker–Prager associative elastic–plastic material model for the ceramic phase. A two stages plasticity has been obtained in which inelastic strain develops in the metal phase followed by a fully plastic response. A comparison with a finite element model of the stress–strain response of an axisymmetric unit cell has been carried out with the purpose to validate the homogenization based modeling presented in the paper. Plastic parameters of a Drucker–Prager yield surface for the homogenized composite have been calculated at different materials compositions. Associative Drucker–Prager plasticity has been found to be accurate for high ceramic content.     

Material sensitivity analysis in homogenization of linear elastic composites

Summary  The main goal of the paper is to present theoretical aspects and the finite element method (FEM) implementation of the sensitivity analysis in homogenization of composite materials with linear elastic components, using effective modules approach. The deterministic sensitivity analysis of effective material properties is presented in a general form for an n-components periodic composite, and is illustrated by the examples of 1D as well as of 2D heterogeneous structures. The results of the sensitivity analysis presented in the paper confirm the usefulness of the homogenization method in computational analysis of composite materials the method may be applied to computational optimization of engineering composites, to the shape-sensitivity studies and, after some probabilistic extensions, to stochastic sensitivity analysis of random composites. Received 10 November 2000; accepted for publication 24 April 2001     

Limit analysis of composite materials based on an ellipsoid yield criterion

Employing repeating unit cell (RUC) to represent the microstructure of periodic composite materials, this paper develops a numerical technique to calculate the plastic limit loads and failure modes of composites by means of homogenization technique and limit analysis in conjunction with the displacement-based finite element method. With the aid of homogenization theory, the classical kinematic limit theorem is generalized to incorporate the microstructure of composites. Using an associated flow rule, the plastic dissipation power for an ellipsoid yield criterion is expressed in terms of the kinematically admissible velocity. Based on nonlinear mathematical programming techniques, the finite element modelling of kinematic limit analysis is then developed as a nonlinear mathematical programming problem subject to only a small number of equality constraints. The objective function corresponds to the plastic dissipation power which is to be minimized and an upper bound to the limit load of a composite is then obtained. The nonlinear formulation has a very small number of constraints and requires much less computational effort than a linear formulation. An effective, direct iterative algorithm is proposed to solve the resulting nonlinear programming problem. The effectiveness and efficiency of the proposed method have been validated by several numerical examples. The proposed method can provide theoretical foundation and serve as a powerful numerical tool for the engineering design of composite materials.     

Macroscopic behavior and field fluctuations in viscoplastic composites: Second-order estimates versus full-field simulations

This work presents a combined numerical and theoretical study of the effective behavior and statistics of the local fields in random viscoplastic composites. The full-field numerical simulations are based on the fast Fourier transform (FFT) 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 estimates follow from the so-called “second-order” procedure [Ponte Castañeda, P., 2002a. Second-order homogenization estimates for nonlinear composites incorporating field fluctuations: I—Theory. J. Mech. Phys. Solids 50, 737-757]. Two-phase fiber composites with power-law phases are considered in detail, for two different heterogeneity contrasts corresponding to fiber-reinforced and fiber-weakened composites. Both the FFT simulations and the corresponding “second-order” estimates show that the strain-rate fluctuations in these systems increase significantly, becoming progressively more anisotropic, with increasing nonlinearity. In fact, the strain-rate fluctuations tend to become unbounded in the limiting case of ideally plastic composites. This phenomenon is shown to correspond to the localization of the strain field into bands running through the composite along certain preferred orientations determined by the loading conditions. The bands tend to avoid the fibers when they are stronger than the matrix, and to pass through the fibers when they are weaker than the matrix. In general, the “second-order” estimates are found to be in good agreement with the FFT simulations, even for high nonlinearities, and they improve, often in qualitative terms, on earlier nonlinear homogenization estimates. Thus, it is demonstrated that the “second-order” method can be used to extract accurate information not only for the macroscopic behavior, but also for the anisotropic distribution of the local fields in nonlinear composites.     

Nonlinear micromechanical formulation of the high fidelity generalized method of cells

The recent High Fidelity Generalized Method of Cells (HFGMC) micromechnical modeling framework of multiphase composites is formulated in a new form which facilitates its computational efficiency that allows an effective multiscale material–structural analysis. Towards this goal, incremental and total formulations of the governing equations are derived. A new stress update computational method is established to solve for the nonlinear material constituents along with the micromechanical equations. The method is well-suited for multiaxial finite increments of applied average stress or strain fields. Explicit matrix form of the HFGMC model is presented which allows an immediate and convenient computer implementation of the offered method. In particular, the offered derivations provide for the residual field vector (error) in its incremental and total forms along with an explicit expression for the Jacobian matrix. This enables the efficient iterative computational implementation of the HFGMC as a stand alone. Furthermore, the new formulation of the HFGMC is used to generate a nested local-global nonlinear finite element analysis of composite materials and structures. Applications are presented to demonstrate the efficiency of the proposed approach. These include the behavior of multiphase composites with nonlinearly elastic, elastoplastic and viscoplastic constituents.     

Asymptotic homogenization of three-dimensional thermoelectric composites

Thermoelectric composites are promising for high efficiency energy conversion between thermal flows and electric conduction, though their effective behaviors remain poorly understood due to nonlinear thermoelectric coupling. In this paper, we develop an asymptotic homogenization theory to analyze the effective behavior of three-dimensional (3D) thermoelectric composites, built on the observation that the equations governing microscopic field fluctuations in the composite are actually linear instead of nonlinear after separation of length scales. A set of solutions similar to Green's function method are used to construct the unit cell problem, and appropriate interfacial continuity conditions and boundary conditions are derived. The homogenized governing equations are then developed for thermoelectric composites, and they are further reduced for a special case wherein the heat flow and electric conduction in the composite remains one-dimensional (1D) at macroscopic scale, even though the composite itself is 3D in general. The general homogenization theory is implemented using finite element method, and a key constant in the constructed solutions is determined using the reformulated eigenvalue problem. The algorithm is validated, and is applied for a number of case studies for the effective behavior of thermoelectric composites.     

Homogenization of elasto-plastic composites coupled with a nonlinear finite element analysis of the equivalent inclusion problem

This paper deals with the micromechanical modeling of particle reinforced elasto-plastic composites under general non-monotonic loading histories. Incremental mean-field (MF) homogenization models offer an excellent cost-effective solution, however there are cases where their predictions are inaccurate. Here, we assess the applicability of the equivalent inclusion representation, which sustains many homogenization schemes. To this end, MF models are fully coupled with a finite element (FE) solution of the equivalent inclusion problem (EIP). Consequently, Eshelby’s tensor is not used and most (but not all) approximations involved in the generalization of MF models from linear elasticity to the nonlinear regime are avoided. The proposal is implemented for Mori-Tanaka (M-T) and dilute inclusion models and applied to several composite systems with elasto-plastic matrix and spherical or ellipsoidal particles, subjected to various loadings (tension, plane strain, cyclic tension/compression). The predictions are verified against reference full-field FE simulations of multiparticle cells. Results show that the M-T model coupled with the nonlinear FE solution of the EIP is very accurate at the macro level up to 25% volume fraction of reinforcement, while the phase averages remain accurate as long as the volume fraction does not exceed 15%. The strain concentration tensor computed almost exactly from single inclusion FE analysis is compared against approximate expressions assumed by classical MF models. Implications for the development of advanced MF homogenization models are discussed.     

Identification of the elastic stiffness of composites using the virtual fields method and digital image correlation

This paper presents an effective methodology for characterizing the mechanical parameters of composites using digital image correlation combined with the virtual fields method.By using a three-point bending test configuration,this method can identify all mechanical parameters of the material with merely a single test.Successful results verified that this method is especially effective for characterizing composite materials.In this study,the method is applied to measure the orthotropic elastic parameters of fiber-reinforced polymer-matrix composites before and after the hygrothermal aging process.The results indicate that the hygrothermal aging environment significantly influences the mechanical property of a composite.The components of the parameters in the direction of the fiber bundle decreased significantly.From the accuracy analysis,we found that the actual measurement accuracy is sensitive to a shift of the horizontal edges and rotation of the vertical edges.     

A micromechanics-based nonlocal constitutive equation incorporating three-point statistics for random linear elastic composite materials

A variational formulation employing the minimum potential and complementary energy principles is used to derive a micromechanics-based nonlocal constitutive equation for random linear elastic composite materials, relating ensemble averages of stress and strain in the most general situation when mean fields vary spatially. All information contained in the energy principles is retained; we employ stress polarization trial fields utilizing one-point statistics so that the resulting nonlocal constitutive equation incorporates up through three-point statistics. The variational structure is developed first for arbitrary heterogeneous linear elastic materials, then for randomly inhomogeneous materials, then for general n-phase composite materials, and finally for two-phase composite materials, in which case explicit variational upper and lower bounds on the nonlocal effective modulus tensor operator are derived. For statistically uniform infinite-body composites, these bounds are determined even more explicitly in Fourier transform space. We evaluate these in detail in an example case: longitudinal shear of an aligned fiber or void composite. We determine the full permissible ranges of the terms involving two- and three-point statistics in these bounds, and thereby exhibit explicit results that encompass arbitrary isotropic in-plane phase distributions; we also develop a nonlocal “Milton parameter”, the variation of whose eigenvalues throughout the interval [0, 1] describes the full permissible range of the three-point term. Example plots of the new bounds show them to provide substantial improvement over the (two-point) Hashin–Shtrikman bounds on the nonlocal operator tensor, for all permissible values of the two- and three-point parameters. We next discuss further applications of the general nonlocal operator bounds: to any three-dimensional scalar transport problem e.g. conductivity, for which explicit results are given encompassing the full permissible ranges of the two- and three-point statistics terms for arbitrary three-dimensional isotropic phase distributions; and to general three-dimensional composites, where explicit results require future research. Finally, we show how the work just summarized, treating elastostatics, can be generalized to elastodynamics, first in general, then explicitly for the longitudinal shear example.     

基于均匀化理论韧性复合材料塑性极限分析

运用细观力学中的均匀化方法分析了韧性复合材料的塑性极限承载能力.从反映复合材料细观结构的代表性胞元入手,将均匀化理论运用到塑性极限分析中,计算由理想刚塑性、Mises组分材料构成的复合材料的极限承载能力.运用机动极限方法和有限元技术,最终将上述问题归结为求解一组带等式约束的非线性数学规划问题,并采用一种无搜索直接迭代算法求解.为复合材料的强度分析提供了一个有效手段.     

Recoverable strains in composite shape memory alloys

New upper bounds are proposed for a generic problem of geometric compatibility, which covers the problem of bounding the effective recoverable strains in composite shape memory alloys (SMAs), such as polycrystalline SMAs or rigidly reinforced SMAs. Both the finite deformation and infinitesimal strain frameworks are considered. The methodology employed is a generalization of a homogenization approach introduced by Milton and Serkov [2000. Bounding the current in nonlinear conducting composites. J. Mech. Phys. Solids 48, 1295-1324] for nonlinear composites in infinitesimal strains. Some analytical and numerical examples are given to illustrate the method.     

A mesomechanical approach to inhomogeneous particulate composites undergoing localized damage:: part I—a mesodomain simulation

This paper presents a mesoscale homogenization methodology to deal with the damage localization in inhomogeneous particulate-reinforced composites. An effective local particle volume fraction, f ^v_loc,eff, which includes the particle size effect is suggested to characterize the damage-inducing aspect of clusters. The clustering regions in a particulate composite accommodating the saturated local damage at an applied stress amplitude are simulated by mesodomains, subdomains of homogeneous medium and homogeneous damage distribution. The size of the mesodomains is determined by the transition condition from local damage to global damage via macro-mechanics. The mesodomain positions are found in a materials science manner by mapping the area contours of f ^v_loc,eff for the composite. Transformation of a clustering composite into a two-homogeneous phase material enables one to appropriately illustrate the local constitutive behaviours, and paves the way to predict saturated local damage life.     

Homogenization of nonlinear problems in the mechanics of composites

A further development of the homogenization method is proposed to solve the physically nonlinear equilibrium problems for the laminated plates or the plates made of functionally graded materials. In the linear case, according to this method, the corresponding solution is a superposition of the solution to the global problem in the entire domain and the solution to the local problem in a representative domain, e.g., in a periodicity cell. In the nonlinear case, such a superposition is not valid, which complicates the application of the homogenization method. In order to eliminate this difficulty, it is possible to combine the homogenization method and the linearization method when solving a boundary value problem or a variational problem. In the mechanics of deformable solids, the constitutive relations can be considered as equations with respect to velocities or the stress and strain differentials in time or in the loading parameter. When these equations are linear with respect to velocities, it is possible to use the homogenization method. In this paper such an approach is illustrated by the example of a symmetric laminated plate bent under a uniformly distributed time-dependent load.     

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.