Objective evaluation of linearization procedures in nonlinear homogenization: A methodology and some implications on the accuracy of micromechanical schemes

A systematic methodology for an accurate evaluation of various existing linearization procedures sustaining mean fields theories for nonlinear composites is proposed and applied to recent homogenization methods. It relies on the analysis of a periodic composite for which an exact resolution of both the original nonlinear homogenization problem and the linear homogenization problems associated with the chosen linear comparison composite (LCC) with an identical microstructure is possible. The effects of the sole linearization scheme can then be evaluated without ambiguity. This methodology is applied to three different two-phase materials in which the constitutive behavior of at least one constituent is nonlinear elastic (or viscoplastic): a reinforced composite, a material in which both phases are nonlinear and a porous material. Comparisons performed on these three materials between the considered homogenization schemes and the reference solution bear out the relevance and the performances of the modified second-order procedure introduced by Ponte Castañeda in terms of prediction of the effective responses. However, under the assumption that the field statistics (first and second moments) are given by the local fields in the LCC, all the recent nonlinear homogenization procedures still fail to provide an accurate enough estimate of the strain statistics, especially for composites with high contrast.     

A nonlinear homogenization procedure for periodic masonry

The present paper deals with the problem of the determination of the in-plane behavior of masonry material. The masonry is considered as a composite material composed by a regular distribution of blocks connected by horizontal and vertical mortar joints. The overall constitutive relationships of the regular masonry are derived by a rational micromechanical and homogenization procedure. Linear elastic constitutive relationship is considered for the blocks, while a new special nonlinear constitutive law is proposed for the mortar joints. In particular, a mortar constitutive law, which accounts for the coupling of the damage and friction phenomena occurring during the loading history, is proposed; the developed model is based on an original micromechanical analysis of the damage process of the mortar joint. Then, an effective nonlinear homogenization procedure, representing the main novelty of the paper, is proposed; it is based on the transformation field analysis, using the technique of the superposition of the effects and the finite element method. The presented methodology is implemented in a numerical code. Finally, numerical applications are performed in order to assess the performances of the proposed procedure in reproducing the mechanical behavior of masonry material.     

Local probabilistic homogenization of two-dimensional model foams accounting for micro structural disorder

The objective of the present study is a numerical analysis of disorder effects in solid structural foams caused by their random irregular micro structure. Using a strain energy based concept, the effective material response is computed in a geometrically non-linear homogenization analysis. The probabilistic homogenization is based on the analysis of a large scale statistically representative volume element. The stochastic information about the scatter in the material response on the lowest possible level is generated by a subsequent division of the representative volume element in substructures consisting of a single cell wall intersection and parts of the adjacent cell walls. For each of the substructures, a homogenization analysis is performed. The results for the local effective stress and strain components are evaluated by means of stochastic methods. The approach is illustrated by a number of examplary studies on the uncertainty of the effective material response of two-dimensional model foams with linear and non-linear elastic material behavior on the cell wall level of structural hierarchy.     

Numerical simulation and homogenization of two-phase flow in heterogeneous porous media

A mathematically rigorous method of homogenization is presented and used to analyze the equivalent behavior of transient flow of two incompressible fluids through heterogeneous media. Asymptotic expansions and H-convergence lead to the definition of a global or effective model of an equivalent homogeneous reservoir. Numerical computations to obtain the homogenized coefficients of the entire reservoir have been carried out via a finite element method. Numerical experiments involving the simulation of incompressible two-phase flow have been performed for each heterogeneous medium and for the homogenized medium as well as for other averaging methods. The results of the simulations are compared in terms of the transient saturation contours, production curves, and pressure distributions. Results obtained from the simulations with the homogenization method presented show good agreement with the heterogeneous simulations.     

Numerical homogenization of elastic and thermal material properties for metal matrix composites (MMC)

A two-scale material modeling approach is adopted in order to determine macroscopic thermal and elastic constitutive laws and the respective parameters for metal matrix composite (MMC). Since the common homogenization framework violates the thermodynamical consistency for non-constant temperature fields, i.e., the dissipation is not conserved through the scale transition, the respective error is calculated numerically in order to prove the applicability of the homogenization method. The thermomechanical homogenization is applied to compute the macroscopic mass density, thermal expansion, elasticity, heat capacity and thermal conductivity for two specific MMCs, i.e., aluminum alloy Al2024 reinforced with 17 or 30 % silicon carbide particles. The temperature dependency of the material properties has been considered in the range from 0 to \(500{\,}^\circ \mathrm {C}\), the melting temperature of the alloy. The numerically determined material properties are validated with experimental data from the literature as far as possible.     

Cosserat model for periodic masonry deduced by nonlinear homogenization

The paper deals with the problem of the determination of the in-plane behavior of periodic masonry material. The macromechanical equivalent Cosserat medium, which naturally accounts for the absolute size of the constituents, is derived by a rational homogenization procedure based on the Transformation Field Analysis. The micromechanical analysis is developed considering a Cauchy model for masonry components. In particular, a linear elastic constitutive relationship is considered for the blocks, while a nonlinear constitutive law is adopted for the mortar joints, accounting for the damage and friction phenomena occurring during the loading history. Some numerical applications are performed on a Representative Volume Element characterized by a selected commonly used texture, without performing at this stage structural analyses. A comparison between the results obtained adopting the proposed procedure and a nonlinear micromechanical Finite Element Analysis is presented. Moreover, the substantial differences in the nonlinear behavior of the homogenized Cosserat material model with respect to the classical Cauchy one, are illustrated.     

岩土介质中局部化变形研究的一些新进展

岩土中的剪切带产生通常是断裂和失稳的前兆, 常引起滑坡等地 质灾害, 对油气或水合物等矿藏的开采容易导致压实带形成. 压实带的产生使岩土孔隙减小, 颗粒破碎, 从而导致沉陷、对油气的挤压等.本文重点对剪切带的宽度和压实带的研究现状作综述, 阐述在该问题的研究上从微观 观测到宏观表述的重要性. 通过分析, 认为还有如下问题值得研究: 怎样 准确得出局部化变形带的宽度和方向?局部化变形带形成后孔隙水和油气的渗流会发生什么 变化?颗粒层次的微结构特性是通过什么方式影响宏观表现的等等.     

Distribution-enhanced homogenization framework and model for heterogeneous elasto-plastic problems

Multi-scale computational models offer tractable means to simulate sufficiently large spatial domains comprised of heterogeneous materials by resolving material behavior at different scales and communicating across these scales. Within the framework of computational multi-scale analyses, hierarchical models enable unidirectional transfer of information from lower to higher scales, usually in the form of effective material properties. Determining explicit forms for the macroscale constitutive relations for complex microstructures and nonlinear processes generally requires numerical homogenization of the microscopic response. Conventional low-order homogenization uses results of simulations of representative microstructural domains to construct appropriate expressions for effective macroscale constitutive parameters written as a function of the microstructural characterization. This paper proposes an alternative novel approach, introduced as the distribution-enhanced homogenization framework or DEHF, in which the macroscale constitutive relations are formulated in a series expansion based on the microscale constitutive relations and moments of arbitrary order of the microscale field variables. The framework does not make any a priori assumption on the macroscale constitutive behavior being represented by a homogeneous effective medium theory. Instead, the evolution of macroscale variables is governed by the moments of microscale distributions of evolving field variables. This approach demonstrates excellent accuracy in representing the microscale fields through their distributions. An approximate characterization of the microscale heterogeneity is accounted for explicitly in the macroscale constitutive behavior. Increasing the order of this approximation results in increased fidelity of the macroscale approximation of the microscale constitutive behavior. By including higher-order moments of the microscale fields in the macroscale problem, micromechanical analyses do not require boundary conditions to ensure satisfaction of the original form of Hill's lemma. A few examples are presented in this paper, in which the macroscale DEHF model is shown to capture the microscale response of the material without re-parametrization of the microscale constitutive relations.     

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.     

Assessment of material uncertainties in solid foams based on local homogenization procedures

The present study is concerned with a numerical prediction and assessment of uncertainties in the macroscopic material properties of solid foams. The material properties are determined by means of a homogenization analysis considering a large scale representative volume element. The microstructure for the representative volume element is determined randomly using a Voronoï tesselation in Laguerre geometry with prescribed cell size distribution. For assessment of the scatter in the effective material response, the homogenization scheme is applied to subsets of the large scale representative volume element. By this means, an interrelation between the local microstructural characteristics and the local mesoscopic material response is established. In a first approach, the individual cells of the foam microstructure are employed as testing volume elements. As an alternative, a moving window technique is applied. The sets of homogenization results obtained by both approaches are evaluated by stochastic methods. For the local effective properties, a distinct scatter is observed. The results in both cases reveal that the local porosity is the most important influence parameter. For the microstructures investigated, only weak local correlations of the effective stiffnesses with a rapid spatial decay of the correlation is observed.     

Thermomechanical Multiscale Constitutive Modeling: Accounting for Microstructural Thermal Effects

In this work we present a thermomechanical multiscale constitutive model for materials with microstructure. In these materials thermal effects at microscale have an impact on the effective macroscopic stress. As a result, it turns out that the homogenized stress depends upon the macroscopic temperature and its gradient. In order to allow this interplay to be thermodynamically valid, we resort to a macroscopic extended thermodynamics whose elements are derived from the microscopic behavior using homogenization concepts. Hence, the thermodynamics implications of this new class of multiscale models are discussed. A variational approach based on the Hill–Mandel Principle of Macro-homogeneity, and which makes use of the volume averaging concept over a local representative volume element (RVE), is employed to derive the thermal and mechanical equilibrium problems at the RVE level and the corresponding homogenization expressions for the effective heat flux and stress. The material behavior at the RVE level is described through standard phenomenological constitutive models. To sum up, the novel contribution of the model presented here is that it allows to include the microscopic temperature fluctuation field, obtained from the multiscale thermal analysis, in the micro-mechanical problem at the RVE level while keeping thermodynamic consistency.     

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.     

Multiscale homogenization of n-component composites with semi-elliptical random interface defects

Effective elastic characteristics of periodic multicomponent composite materials with random interface defects are studied in the paper. The defects are assumed to be semi-elliptical and lying with major semi axes along the interfaces, where minor and major semi-axes as well as the defects number are given as input random variables. The homogenization approach has a multiscale character—some algebraic approximation is used first to calculate effective elastic parameters of the interphase including all defects located at the same interface. Equations for interphase random elastic parameters are obtained using MAPLE symbolic mathematics in conjunction with probabilistic generalized perturbation method. A different homogenization method is applied at the micro scale, where the cell problem is solved numerically using the Finite Element Method (FEM) program. Since the composites considered exhibit random variations of both elastic properties and the interface defects, the overall homogenized characteristics must be obtained as random quantities, which is realized on the micro scale by the Monte-Carlo simulation. The proposed interface defects model obeys the porosity effects resulting from the nature of some matrices in engineering composites as well as the interface cracks appearing as a result of composites ageing during static or fatigue fracture.     

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种典型本构关系的周期性多孔材料平面问题进行了计算, 对比细致模型的应变能、最大位移和等效泊松比, 对程序及迭代方法的准确性进行了验证. 之后对一种三元橡胶基复合材料进行多尺度均质化, 将其分为芯丝尺度和层间尺度. 用线弹性的均质化方法得到了芯丝尺度的等效弹性参数, 并将其作为层间尺度的材料参数. 在层间尺度应用非线性弹性均质化方法对结构进行计算, 得到材料的宏观等效性能, 并以实验结果为基准进行评价.      

A Lagrangian-based continuum homogenization approach applicable to molecular dynamics simulations

The continuum notions of effective mechanical quantities as well as the conditions that give meaningful deformation processes for homogenization problems with large deformations are reviewed. A continuum homogenization model is presented and recast as a Lagrangian-based approach for heterogeneous media that allows for an extension to discrete systems simulated via molecular dynamics (MD). A novel constitutive relation for the effective stress is derived so that the proposed Lagrangian-based approach can be used for the determination of the “stress–deformation” behavior of particle systems. The paper is concluded with a careful comparison between the proposed method and the Parrinello–Rahman approach to the determination of the “stress–deformation” behavior for MD systems. When compared with the Parrinello–Rahman method, the proposed approach clearly delineates under what conditions the Parrinello–Rahman scheme is valid.