Computational Multi-Scale Stability Analysis of Periodic Electroactive Polymer Composites at Finite Strains

This work is dedicated to multi-scale stability analysis, especially macroscopic and microscopic stability analysis of periodic electroactive polymer (EAP) composites with embedded fibers. Computational homogenization is considered to determine the response of materials at macro-scale depending on the selected representative volume element (RVE) at micro-scale [4, 5]. The quasi-incompressibility condition is considered by implementing a four-field variational formulation on the RVE, see [7]. Based on the works [1–3, 6, 8] the macroscopic instabilities are determined by the loss of strong ellipticity of homogenized moduli. On the other hand, the bifurcation type microscopic instabilities are treated exploiting the Bloch-Floquet wave analysis in context of finite element discretization, which allows to detect the changed critical size of periodicity of the microstructure and critical macroscopic loading points. Finally, representative numerical examples are given which demonstrate the onset of instability surfaces, the stable macroscopic loading ranges, and further a periodic buckling mode at a microscopic instability point is presented. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Two-scale computational homogenization of magneto-electric composites

This contribution presents a two-scale computational homogenization framework for the micro-macro simulation of magneto-electro-mechanically coupled materials. Energetically consistent micro-macro transition conditions will be derived from a generalized form of the classical Hill-Mandel condition. A focus of the work is on the computation of effective magneto-electric moduli which are derived on the basis of an algorithmically consistent linearization of the macroscopic field equations. The method will be applied to the homogenization of magneto-electric composites which are composed of piezomagnetic and piezoelectric phases. The effective magneto-electric moduli of two-phase composites will be computed. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Microscale modelling and homogenization of fiber structured materials

Various phenomena occurring on the macrosscale result from physical and mechanical behaviour on the microscale [1]. For the mechanical modeling and simulation of the heterogeneous composition of fiber structured material, in addition to the material properties the contact between the fibers has to be taken into account. The material behaviour is strongly influenced by the material properties of the fiber, but also by the geometrical structure. Periodically arranged fibers like woven, knitted or plaited fabrics and randomly oriented ones like fleece can be distinguished in their arrangement. In consideration of different lengthscales the problem involves, it is necessary to introduce a multiscale approach based on the concept of a representative volume element (RVE). The macro-micro scale transition requires a method to impose the deformation gradient on the RVE by suited boundary conditions. The reversing scale transition, based on the HILL-MANDEL condition, requires the equality of the macroscopic average of the variation of work on the RVE and the local variation of the work on the macroscale [2]. For the micro-macro transition the averaged stresses have to be extracted by a homogenization scheme. From these results an effective material law can be derived. Beside the theoretical aspects, we present the stress-strain relation for RVE-models and different boundary conditions. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A Variational Homogenization Approach to Electromechanical Hystereses Caused by Domain Wall Movement

In recent years, increasing interest in so-called smart materials such as ferroelectric polymers and ceramics has been shown. Those materials are used in various actuators, sensors, and also in medical devices. In this paper, we outline a micro-macro approach to the modeling of macroscopic hystereses which directly takes into account the microstructural evolution of electrically poled domains. To this end, an incremental variational formulation for a gradient-type phase field model is developed and exploited for the simulation of electromechanically coupled problems. The formulation determines the hysteretic response of the material in terms of an energy-enthalpy and a dissipation function which both depend on the microscopic remanent polarization treated as an order parameter. The gradient-type balance law for the phase field can be considered as a generalization of Biot's equation for standard dissipative materials and may be related to the classical Ginzburg-Landau equation. Furthermore, the variational formulation serves as natural starting point for a compact and symmetric finite element implementation of the coupled micromechanical problem covering the displacement, the electric potential, and the microscopic polarization vector. For this three-field scenario we develop a variational-based homogenization method which determines the overall macroscopic hysteretic properties of a polycrystalline aggregate. The proposed computational method can be used as a numerical laboratory for the improvement of microstructural properties. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Estimation of material properties of cancellous bone using multiscale FEM

Cancellous bone is a spongy type of bone with voids filled by blood marrow. Without much loss of generality it can be modeled as a material with periodic microstructure where overall parameters can be calculated using homogenization methods. Here the multiscale finite element method is applied and the assumed representative volume element (RVE) is a cube with solid frame and fluid core. From the point of view of the finite element method the RVE is a combination of solid and shell elements. As the acoustic excitation is considered, a complex stiffness matrix and complex displacements appear in the solution of the problem. Calculation of overall properties is repeated for different geometries of the solid frame. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Topological sensitivity analysis of a multi‐scale constitutive model considering a cracked microstructure

This paper deals with the sensitivity analysis of the macroscopic elasticity tensor to topological microstructural changes of the underlying material. In particular, the microstucture is topologicaly perturbed by the nucleation of a small circular inclusion. The derivation of the proposed sensitivity relies on the concept of topological derivative, applied within a variational multi‐scale constitutive framework where the macroscopic strain and stress at each point of the macroscopic continuum are defined as volume averages of their microscopic counterparts over a representative volume element (RVE) of material associated with that point. We consider that the RVE can contain a number of voids, inclusions and/or cracks. It is assumed that non‐penetration conditions are imposed at the crack faces, which do not allow the opposite crack faces to penetrate each other. The derived sensitivity leads to a symmetric fourth‐order tensor field over the unperturbed RVE domain, which measures how the macroscopic elasticity parameters estimated within the multi‐scale framework changes when a small circular inclusion is introduced at the micro‐scale level. Copyright © 2009 John Wiley & Sons, Ltd.     

The heterogeneous multiscale finite element method FE-HMM for the homogenization of linear elastic solids

The objective of the present paper is a formulation of the Heterogeneous Multiscale Finite Element Method (FE-HMM) for the homogenization of linear elastic solids in a geometrical linear frame, and doing so, for the first time, of a vector-valued field problem. The macro stiffness is estimated by stiffness sampling on heterogeneous microdomains in terms of a modified quadrature formula, which implies an equivalence of energy densities of the microscale with the macroscale. Existing a-priori estimates are assessed and used for optimal micro-macro refinement strategies in the H¹-norm and the L²-norm. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

An adaptive approach for modeling a fiber-matrix composite with the FE2 method

In materials with a complicated microstructre [1], the macroscopic material behaviour is unknown. In this work a Fiber-Matrix composite is considered with elasto-plastic fibers. A homogenization of the microscale leads to the macroscopic material properties. In the present work, this is realized in the frame of a FE² formulation. It combines two nested finite element simulations. On the macroscale, the boundary value problem is modelled by finite elements, at each integration point a second finite element simulation on the microscale is employed to calculate the stress response and the material tangent modulus. One huge disadvantage of the approach is the high computational effort. Certainly, an accompanying homogenization is not necessary if the material behaves linear elastic. This motivates the present approach to deal with an adaptive scheme. An indicator, which makes use of the different boundary conditions (BC) of the BVP on microscale, is suggested. The homogenization with the Dirichlet BC overestimates the material tangent modulus whereas the Neumann BC underestimates the modulus [2]. The idea for an adaptive modeling is to use both of the BCs during the loading process of the macrostructure. Starting initially with the Neumann BC leads to an overestimation of the displacement response and thus the strain state of the boundary value problem on the macroscale. An accompanying homogenization is performed after the strain reaches a limit strain. Dirichlet BCs are employed for the accompanying homogenization. Some numerical examples demonstrate the capability of the presented method. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Influence of the nonlinear matrix material behaviour on the effective properties of textile-reinforced thermoplastics

For a consistent lightweight design the consideration of the nonlinear macroscopic material behaviour of composites, which is amongst others driven by damage and strain-rate effects on the mesoscale, is required. Therefore, a modelling approach using numerical homogenization techniques is applied to predict the effective nonlinear material behaviour of the composite based on the finite element simulation of a representative volume element (RVE). In this RVE suitable constitutive relations account for the material behaviour of each constituents. While the reinforcing glass fibres are assumed to remain linear elastic, a viscoplastic constitutive law is applied to represent the strain-rate dependent, inelastic deformation of the matrix material. In order to analyse the influence of the nonlinear matrix material behaviour on the global mechanical response of the composite, effective stress-strain-curves are computed for different load cases and compared to experimental observations. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Effective Dynamic Material Properties of Foam-like Microstructures

The effective material parameters of a microstructured material can be found using homogenization procedures based on calculations of a Representative Volume Element (RVE) of the material. In our approach the RVE is calculated in frequency domain and inertia is taken into account, leading to a frequency dependent behavior of the RVE.With the frequency response of the RVE, effective dynamic properties of the material are calculated using an optimization procedure. Due to the frequency dependent material behavior on the microscale a viscoelastic constitutive equation is applied on the macroscale. An example calculation is presented for an auxetic 2-D foam-like microstructure which is modelled as a plane frame structure. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Designing Statistically Similar RVEs for 3D Dual Phase Steel Microstructures

This contribution presents a method for the construction of three-dimensional Statistically Similar Representative Volume Elements (SSRVEs) for dual phase steels (DP steels). From such kind of advanced high strength steels, enhanced material properties are observed, which originate in the interaction of the individual constituents of the material on the microscale. Our aim is to directly incorporate the microstructure in the material modeling, which can be accomplished by applying i. e. the FE² method. A RVE representing the real material is used in the microscopic boundary value problem, which is solved at each macroscopic integration point. Since such RVEs usually exhibit a high complexity due to the underlying real microstructure, high computational costs are a drawback of the approach. We replace this RVE with a SSRVE, which has a lower complexity but which is still able to represent the mechanical behavior of the RVE and thus of the real microstructure. Virtual experiments show the performance of the method. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical Prediction of Macroscopic Material Failure

The aim of this contribution is the numerical determination of macroscopic material properties based on constitutive relationships characterising the microscale. A macroscopic failure criterion is computed using a three dimensional finite element formulation. The proposed finite element model implements the Strong Discontinuity Approach (SDA) in order to include the localised, fully nonlinear kinematics associated with the failure on the microscale. This numerical application exploits further the Enhanced–Assumed–Strain (EAS) concept to decompose additively the deformation gradient into a conforming part corresponding to a smooth deformation mapping and an enhanced part reflecting the final failure kinematics of the microscale. This finite element formulation is then used for the modelling of the microscale and for the discretisation of a representative volume element (RVE). The macroscopic material behaviour results from numerical computations of the RVE. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Micro-mechanical analysis of creep behavior in a multipass weld

A method of evaluating creep response of the multipass weld based on the micro-macro mechanics approach is introduced. Multipass weld microstructure that consists from columnar, coarse and fine grained zones is considered. Materials of these constituents assumed to be isotropic. Weld metal properties of inelastic behavior have general type of symmetry and are described by the anisotropic creep constitutive model. To model the microstructure of the multipass weld metal the representative volume element (RVE) is created for CAE Abaqus. Material properties of weld metal grain type zones are taken from the experiments. Numerical tests on uniform loading of the RVE are performed. Creep material properties for equivalent weld material are found for welds with different number of passes. The symmetry type of the creep material properties of multi-pass weld are evaluated for the equivalent weld material. As an example of macro model analysis of the welding, the creep calculation of the cylindrical shell with the welding under the uniform inner pressure is performed. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A posteriori error analysis of the heterogeneous multiscale method for homogenization problems

In this Note we derive a posteriori error estimates for a multiscale method, the so-called heterogeneous multiscale method, applied to elliptic homogenization problems. The multiscale method is based on a macro-to-micro formulation. The macroscopic method discretizes the physical problem in a macroscopic finite element space, while the microscopic method recovers the unknown macroscopic data on the fly during the macroscopic stiffness matrix assembly process. We propose a framework for the analysis allowing to take advantage of standard techniques for a posteriori error estimates at the macroscopic level and to derive residual-based indicators in the macroscopic domain for adaptive mesh refinement. To cite this article: A. Abdulle, A. Nonnenmacher, C. R. Acad. Sci. Paris, Ser. I 347 (2009).     

Computational homogenization of non-stationary transport processes in masonry structures

A fully coupled transient heat and moisture transport in a masonry structure is examined in this paper. Supported by several successful applications in civil engineering the nonlinear diffusion model proposed by Künzel (1997) [16] is adopted in the present study. A strong material heterogeneity together with a significant dependence of the model parameters on initial conditions as well as the gradients of heat and moisture fields vindicates the use of a hierarchical modeling strategy to solve the problem of this kind. Attention is limited to the classical first order homogenization in a spatial domain developed here in the framework of a two step (meso–macro) multi-scale computational scheme (FE² problem). Several illustrative examples are presented to investigate the influence of transient flow at the level of constituents (meso-scale) on the macroscopic response including the effect of macro-scale boundary conditions. A two-dimensional section of Charles Bridge subjected to actual climatic conditions is analyzed next to confirm the suitability of algorithmic format of FE² scheme for the parallel computing.     

FE-implementation of plane strain gradient elasticity constitutive equations obtained by homogenization

This paper presents a plane strain gradient elasticity model together with a strength criterion based on the maximum hoop stress. It is developed by higher order homogenization based on a cylindrical RVE containing a void and is implemented as mixed–type finite element formulation into the FE–code Abaqus. Numerical simulations are performed in order to study the size effect related to the onset of failure of specimens of porous elastic material. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On the Computational Homogenization of Transient Diffusion Problems

In this contribution, we address the computational treatment of transient diffusion problems with heterogeneous microstructures using first-order homogenization. There, we treat two different cases, firstly, when the transient part at the microscale can be neglected due to the vanishingly small size of the representative volume element (RVE) and secondly, when no steady state is reached at the microscale. This is the case when dealing with a relaxed version of the scale separation condition. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Influence of the Homogenization on the Transient Behaviour of Size Distributed Polycrystals

The deformation in polycrystals is often heterogenous, e.g. due to grain size dependent hardening. In a semi-analytical representative volume element (RVE), a log-normal distributed grain size is assumed together with a grain size dependent local plastic behavior. The numerical results are well approximated by a simple analytical expression. The effect of the homogenization comparison stiffness on the transient behaviour is explained using a simplified localization equation. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Homogenization for rate-independent systems

This paper is devoted to the homogenization for a class of rate-independent systems described by the energetic formulation. The associated nonlinear partial differential system has periodically oscillating coefficients, but has the form of a standard evolutionary variational inequality. Thus, the model applies to standard linearized elastoplasticity with hardening. Using the recently developed methods of two-scale convergence, periodic unfolding and the new introduced one, periodic folding, we show that the homogenized problem can be represented as a two-scale limit which is again an energetic formulation, but now involving the macroscopic variable in the physical domain as well as the microscopic variable in the periodicity cell. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)