A surrogate model for computational homogenization of elastostatics at finite strain using high-dimensional model representation-based neural network

We propose a surrogate model for two-scale computational homogenization of elastostatics at finite strains. The macroscopic constitutive law is made numerically available via an explicit formulation of the associated macroenergy density. This energy density is constructed by using a neural network architecture that mimics a high-dimensional model representation. The database for training this network is assembled through solving a set of microscopic boundary value problems with the prescribed macroscopic deformation gradients (input data) and subsequently retrieving the corresponding averaged energies (output data). Therefore, the two-scale computational procedure for nonlinear elasticity can be broken down into two solvers for microscopic and macroscopic equilibrium equations that work separately in two stages, called the offline and online stages. The finite element method is employed to solve the equilibrium equation at the macroscale. As for microscopic problems, an FFT-based collocation method is applied in tandem with the Newton-Raphson iteration and the conjugate-gradient method. Particularly, we solve the microscopic equilibrium equation in the Lippmann-Schwinger form without resorting to the reference medium. In this manner, the fixed-point iteration that might require quite strict numerical stability conditions in the nonlinear regime is avoided. In addition, we derive the projection operator used in the FFT-based method for homogenization of elasticity at finite strain.     

A nonlocal computational homogenization of softening quasi-brittle materials

In this paper, a computational counterpart of the experimental investigation is presented based on a nonlocal computational homogenization technique for extracting damage model parameters in quasi-brittle materials with softening behavior. The technique is illustrated by introducing the macroscopic nonlocal strain to eliminate the mesh sensitivity in the macroscale level as well as the size dependence of the representative volume element (RVE) in the first-order continuous homogenization. The macroscopic nonlocal strains are computed at each direction, and both the local and nonlocal strains are transferred to the microscale level. Two RVEs with similar geometries and material properties are introduced for each macroscopic Gauss point, in which the microscopic damage variables and the macroscale consistent tangent modulus and its derivatives are obtained by imposing the macroscopic nonlocal strain on the first RVE, and the macroscopic stress is computed by employing the microscopic damage variables and imposing the macroscopic local strain over the second RVE. Finally, numerical examples are solved to illustrate the performance of the proposed nonlocal computational homogenization technique for softening quasi-brittle materials.     

Two-scale homogenization of electromechanically coupled boundary value problems

The contribution addresses a direct micro-macro transition procedure for electromechanically coupled boundary value problems. The two-scale homogenization approach is implemented into a so-called FE²-method which allows for the computation of macroscopic boundary value problems in consideration of microscopic representative volume elements. The resulting formulation is applicable to the computation of linear as well as nonlinear problems. In the present paper, linear piezoelectric as well as nonlinear electrostrictive material behavior are investigated, where the constitutive equations on the microscale are derived from suitable thermodynamic potentials. The proposed direct homogenization procedure can also be applied for the computation of effective elastic, piezoelectric, dielectric, and electrostrictive material properties.     

Inelastic analysis of granular interfaces via computational contact homogenization

The macroscale response of granular contact interfaces is investigated. In order to circumvent the difficulties associated with a direct resolution of such heterogeneous contact problems, where highly mobile particles residing between a deformable body and a rigid surface govern the microscale dynamics, a space–time contact homogenization methodology is developed. The overall approach is based on a separation of spatial as well as temporal scales and proposes an idealized purely frictional macroscale response. The induced macroscale dissipation is directly associated with the microscale dissipation mechanisms due to (i) an inelastic constitutive response for the boundary layer of the deformable body and (ii) frictional interaction among the components of the three‐body contact system. The consequences of a viscoelastic boundary layer that sustains damage due to highly localized deformation in the vicinity of the particles are investigated extensively within a fully nonlinear computational setting that accounts for incompressibility. The effective coefficient of friction that is induced by the homogenization methodology as the fundamental macroscale observable is found to be of a non‐Amontons as well as a non‐Coulomb type. The proposed analysis framework is amenable to a multiscale implementation within a coupled micro–macro approach and yields insight into the macroscopic dynamics of similar heterogeneous interfaces with varying degrees of mobility associated with the roughness features. Copyright © 2010 John Wiley & Sons, Ltd.     

A coupled thermo-chemo-mechanical reduced-order multiscale model for predicting process-induced distortions,residual stresses,and strength

We study residual stresses and part distortion induced by a manufacturing process of a polymer matrix composite and its effect on the component strength. Unlike most of the thermo-chemo-mechanical models in the literature where governing multiphysics equations are directly formulated on the macroscale, we present a multiscale-multiphysics approach. To address the enormous computational complexity involved, a reduced-order homogenization was originally developed for a single physics problem is employed. The proposed reduced-order two-scale thermo-chemo-mechanical model has been validated for predicting part distortion beam strength in three-point bending test. It is shown that while macroscopic stresses are relatively low, and therefore often ignored in practice, stresses at the scale of microconstituents are significant and may have an effect on the overall composite component strength.     

Numerical aspects of non-local modeling of the damage evolution in elastic–plastic materials

The design of mechanical systems in modern industrial plants requires reliable and efficient methods to predict the behavior of structural materials. For complex loading conditions, the behavior of the structural materials is determined by damage evolution, strain rate and temperature. The subject of the article is the modeling of the damage evolution in elastic–plastic materials of structural components, which are utilized at various temperatures. To achieve this goal, a hybrid model of steel cracking is applied. The hybrid model uses a finite element simulation combined with an experimental test realized in the macroscale. By using the hybrid model, the modeling of the damage evolution affords possibilities of determining macroscopic effects of the steel micro-defects. An essence of solving the predicting behavior of structural materials with micro-defects consists in time integration procedures for constitutive equations. In the article a semi-implicit time integration procedure is presented. The semi-implicit time integration procedure is suitable for the inelastic materials (compressible or incompressible) with the combined kinematic–isotropic hardening behavior. Its numerical solutions are stable, namely without the oscillatory behavior. By spatial averaging over a representative volume (RV), the homogenization technique (HT) is used for the defining of non-local variables in the constitutive equations. Evolutionary algorithms (EAs) based on local selections are applied to perform the homogenization technique. Within the framework of the large strain theory, the non-local continuum satisfies the objectivity requirements. Limitations on applicability of the -integral approach to construct crack growth resistance curves are also presented.     

Modeling of Advanced High Strength Steels with the realistic microstructure–strength relationships

The objective of the work is to consider the first-order effects of the realistic microstructure morphology in the macroscale modeling of the multiphase Advanced High Strength Steels (AHSS). Instead of using constitutive equations at macroscale, the strength–microstructure relationship is studied in the forms of micromechanical and multiscale models that do not make considerable simplifications with regard to the microscale geometry and topology. The trade-off between the higher computational time and the higher accuracy has been offset with a stochastic approach in the construction of the microscale models. The multiphase composite effects of AHSS microstructure is considered in realistic microstructural models that are stochastically built from AHSS micrographs. Computational homogenization routines are used to couple micro and macroscale and resultant stress–strain relations are compared for models built with the simplified and idealized geometries of the microstructure. The results from this study show that using a realistic representation of the microstructure, either for DP or TRIP steel, could improve the accuracy of the predicted stress and strain distribution. The resultant globally averaged effective stress and strain fields from realistic microstructure model were able to accurately capture the onset of the plastic instability in the DP steel. It is shown that the macroscale mechanical behavior is directly affected by the level of complexities in the microscale models. Therefore, greater accuracy could be achieved if these stochastic realistic microstructures are used at the microscale models.     

Modeling macroscopic extended continua with the aid of numerical homogenization schemes

Even in the range of small elastic deformations the behavior of foams is not well described by only two elastic constants. Usually the manufacturers give values of the material parameters depending on the loading conditions. This problem is investigated on a microscopic scale by a simple beam model and on the macroscopic scale by an extended continuum model. It has been found that this approach shows the size effect [J. Mater. Sci. 18 (1983) 2572] that cannot be described within the framework of the standard continuum mechanical setting. The existence of the size effect within this model can be explained by independent rotations which do not scale with the displacement field.While macroscopic material parameters are generally unknown for foams the macroscopic properties are derived from the microscale where the parameters are assumed to be known. After evaluation of the microscopic constitutive equations, which are also considered to be known, the quantities are mappped back to the macroscale by a homogenization procedure. This approach is known from literature as FE² model, see e.g. [V. Kouznetsova, Computational homogenization for the multi-scale analysis of multi-phase materials, PhD-thesis, Technical University of Eindhoven, 2002], [Int. J. Numer. Meth. Eng., 54 (2002) 1235] or [Arch. Appl. Mech., 72 (2002) 300]. It is shown that the Cosserat continuum and the FE² model are able to describe the same effects qualitatively.     

Two-scale analysis for deformation-induced anisotropy of polycrystalline metals

The anisotropic macroscopic mechanical behavior of polycrystalline metals is characterized by incorporating the microscopic constitutive model of single crystal plasticity into the two-scale modeling based on the mathematical homogenization theory, which enables us to derive both micro- and macro-scale governing equations. The two-scale simulations are conducted to evaluate the macroscopic anisotropy induced by microscopic plastic deformation histories of the polycrystalline aggregate. In the simulations, the representative volume element (RVE) composed of several crystal grains is uniformly loaded in one direction, unloaded to macroscopically zero stress in a certain stage of deformation and then re-loaded in the different directions. The last re-loading calculations provide different macroscopic responses of the RVE, which can be the appearance of material anisotropy. We then try to examine the effects of the intergranular and intragranular behaviors on the anisotropy by means of various illustrations of microscopic plastic deformation process without referring to the change of crystallographic orientations.     

Multiscale analysis of laminated plates with integrated piezoelectric fiber composite actuators

This study is concerned with the detailed analysis of fiber-reinforced composite plates with integrated piezoceramic fiber composite actuators. A multiscale framework based on the asymptotic expansion homogenization method is used to couple the microscale and macroscale field variables. The microscale fluctuations in the mechanical displacement and electric potential are related to the macroscale deformation and electric fields through 36 distinct characteristic functions. The local mechanical and charge equilibrium equations yield a system of partial differential equations for the characteristic functions that are solved using the finite element method. The homogenized electroelastic properties of a representative material element are computed using the characteristic functions and the material properties of the fiber and matrix. The three-dimensional macroscopic equilibrium equations for a laminated piezoelectric plate are solved analytically using the Eshelby-Stroh formalism. The formulation admits different boundary conditions at the edges and is applicable to thick and thin laminated plates. The microscale stresses and electric displacement in the fibers and matrix are computed from the macroscale fields through interscale transfer operators. The multiscale analysis procedure is illustrated using two model problems. In the first model problem, a simply-supported sandwich plate consisting of a piezoceramic fiber composite shear actuator embedded between two graphite/polymer layers is studied. The second model problem concerns a cantilever graphite/polymer substrate with segmented piezoceramic fiber composite extension actuators attached to its top and bottom surfaces. Results are presented for the homogenized material properties, macroscale deformation, macroscale average stresses and microscale stress distributions.     

An iterative method for coupling of deformation and failure mechanisms on different length scales

An iterative method for coupling of numerical simulations on two length scales is presented. The computations on the microscale and on the macroscale are linked via a suitable macroscopic constitutive law. The parameters of this material law depend on the deformation history and are obtained from simulations using microstructurally representative volume elements (RVEs) subjected to strain paths derived from the associated material points in the macroscopic structure. Thus, different constitutive parameter sets are assigned to different regions of the macrostructure. The microscopic and macroscopic simulations are performed iteratively and interact mutually via the strain paths and the constitutive parameters, respectively. As an example, the strip tension test for a porous material is modelled using the finite element (FE) method. The coupling procedure, the material law and its numerical implementation are described. The method is shown to allow for a detailed simulation of the deformation mechanisms both on the micro‐ and the macroscale as well as for an assessment of their interactions while keeping the computational efforts reasonably low. Copyright © 2000 John Wiley & Sons, Ltd.     

Multiscale prediction of mechanical behavior of ferrite–pearlite steel with numerical material testing

Multiscale mechanical behaviors of ferrite–pearlite steel were predicted using numerical material testing (NMT) based on the finite element method. The microstructure of ferrite–pearlite steel is regarded as a two‐component aggregate of ferrite crystal grains and pearlite colonies. In NMT, the macroscopic stress–strain curve and the deformation state of the microstructure were examined by means of a two‐scale finite element analysis method based on the framework of the mathematical homogenization theory. The microstructure of ferrite–pearlite steel was modeled with finite elements, and constitutive models for ferrite crystal grains and pearlite colonies were prepared to describe their anisotropic mechanical behavior at the microscale level. While the anisotropic linear elasticity and the single crystal plasticity based on representative characteristic length have been employed for the ferrite crystal grains, the constitutive model of a pearlite colony was newly developed in this study. For that reason, the constitutive behavior of the pearlite colony was investigated using NMT on a smaller scale than the scale of the ferrite–pearlite microstructure, with the microstructure of the pearlite colony modeled as a lamellar structure of ferrite and cementite phases with finite elements. On the basis of the numerical results, the anisotropic constitutive model of the pearlite colony was formulated based on the normal vector of the lamella. The components of the anisotropic elasticity were estimated with NMT based on the finite element method, where the elasticity of the cementite phase was numerically evaluated with a first‐principles calculation. Also, an anisotropic plastic constitutive model for the pearlite colony was formulated with two‐surface plasticity consisting of yield functions for the interlamellar shear mode and yielding of the overall lamellar structure. After addressing the microscopic modeling of ferrite–pearlite steel, NMT was performed with the finite element models of the ferrite–pearlite microstructure and with the microscopic constitutive models for each of the components. Finally, the results were compared with the corresponding experimental results on both the macroscopic response and the microscopic deformation state to ascertain the validity of the numerical modeling. Copyright © 2011 John Wiley & Sons, Ltd.     

Multiscale fatigue life prediction model for heterogeneous materials

A multiscale fatigue life prediction model is developed for heterogeneous materials. The proposed model combines a two‐scale asymptotic homogenization approach in time with a ‘block cycle jump’ technique into a unified temporal multiscale framework that can be effectively utilized for arbitrary material architectures and constitutive equations of microphases. The unified temporal multiscale approach in combination with a spatial multiscale approach based on the reduced order homogenization is characterized for high temperature ceramic matrix composites. Copyright © 2012 John Wiley & Sons, Ltd.     

Application of Cosserat Theory to the Modelling of Reinforced Carbon Nananotube Beams

This paper develops a mechanical model for multifunctional reinforced carbon nanotube (CNT) beams. The model is obtained by introducing the couple stresses into the constitutive equations of linear viscoelastic theory. The material functions are determined using the homogenization method.     

A multiscale framework for localizing microstructures towards the onset of macroscopic discontinuity

This paper presents a multiscale computational homogenization model for the post localization behavior of a macroscale domain crossed by a cohesive discontinuity emanating from microstructural damage. The stress–strain and the cohesive macroscopic responses are obtained incorporating the underlying microstructure, in which the damage evolution results in the formation of a strain localization band. The macro structural kinematics entails a discontinuous displacement field and a non-uniform deformation field across the discontinuity. Novel scale transitions are formulated to provide a consistent coupling to the continuous microscale kinematics. From the solution of the micromechanical boundary value problem, the macroscale stress responses at both sides of the discontinuity are recovered, providing automatically the cohesive tractions at the interface. The effective displacement jump and deformation field discontinuity are derived from the same microscale analysis. This contribution focusses on scale transition relations and on the solution procedure at the microlevel; the highlights of the approach are demonstrated on microscale numerical examples. Coupled two-scale solution strategy will be presented in a subsequent paper.     

Framework for deformation induced anisotropy in glassy polymers

In this paper a constitutive model for glassy polymers is developed. Glassy polymers consist of a number of polymer chains that at a microscopic level form a network. If the distribution of the polymer chains shows some preferred direction, the mechanical response at a global macroscopic level will be anisotropic. To incorporate the orientational distribution of the polymer chains, a homogenization procedure involving a chain orientation distribution function was undertaken. When polymers are exposed to external loading, the chains at the microscopic level orient in a certain manner, leading to an evolution of the macroscopic anisotropic properties. This phenomenon was modeled by use of evolution equations for the chains at a microscopic level and are then—by using the orientation distribution function—transformed to the macroscopic level. The theories involved are developed in a large strain setting in which a multiplicative split of the deformation gradient for the elastic-viscoplastic response is adopted. Various numerical experiments were conducted to evaluate the model that was developed.     

基于耦合扩展多尺度有限元方法的功能梯度材料热应力分析

以高效模拟功能梯度材料(FGM)微观非均质性对整体热力学性能的影响为研究目的,通过随机形态描述函数(RMDF)法和体积分数的指数分布建立FGM二维微结构,在此基础上,发展了FGM热应力分析的耦合扩展多尺度有限元方法(CEMsFEM)。该方法基于扩展多尺度有限元方法(EMsFEM)的基本思想,对温度场和位移场构造数值基函数,以把微观非均质材料性质带到宏观响应中。同时为了考虑泊松效应导致的不同方向间的耦合作用,在位移场数值基函数中增加了耦合附加项。通过数值基函数建立宏微观单元信息的映射关系,在宏观尺度求解有效方程,节约计算量。为了更好地考虑微观载荷的影响,把结构的真实响应分解为宏观响应和微观扰动,进一步推导出修正的宏观载荷向量。通过不同体积分数分布的FGM在不同载荷工况下的热应力分析算例验证了本文中方法的正确性和有效性,最后讨论了微结构的尺寸效应对结构热力学响应的影响。