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)     

An Adaptive Finite Element Approach for Brittle Fracture

The extension of the finite element method to take discrete fracture and failure modes into account is a current field of research. In recent times, first results in terms of cohesive element formulations have been introduced into commercial applications. Such element formulations are able to cover the discrete behaviour of interfaces between different materials or the mechanical processes of thin layers. These approaches are not suitable for simulations with unknown crack paths in homogeneous materials, due to the initial elastic phase of the material formulation and the necessity to define potential crack paths a priori. The presented strategy starts with an unextended model and modifies the structure during the computations in terms of an adaptive procedure. The idea is to generate additional elements, based on the cohesive element formulation, to approximate arbitrary crack paths. For this purpose, a failure criterion is introduced. For nodes where the limiting value is reached, cohesive elements are introduced between the volume element boundaries of accordingly facets and corresponding nodes are duplicated. Necessary modifications for this application on system level as well as the element and the material formulation are introduced. By means of some numerical examples, the functionality of the presented procedure is demonstrated. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Determination of Effective Properties of MMC by Computational Homogenization

On the macrolevel Metal Matrix Composites (MMCs) resemble a homogeneous material. However, on the microlevel they show an inhomogeneous microstructure. This paper will have how heterogeneities affect the overall properties and the behaviour of a material (i. e. the effective properties). This is done using computational homogenization techniques. Finite element (FE) simulations were conducted in ABAQUS in connection with MATLAB, using material parameters for aluminium alloy AA2124 and SiC to develop a representative volume element (RVE) of the MMC AMC217xe. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Optimal Experimental Design to Identify the Average Stress-Strain Response in Short Fiber-Reinforced Plastics

We show how to use optimal experimental design methods for the parameter identification of short fiber reinforced plastic (SFRP) materials. The experimental data is given by computer simulations of representative volume elements (RVE) of the SFRP material. The experiments are designed such that a minimal number of RVE simulations is required and that the model response attains a minimal variance for a class of strains and fiber orientations. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Comparison of Statistical Descriptors for the Construction of Statistically Similar RVEs

A method for the construction of three-dimensional statistically similar representative volume elements (SSRVEs) is presented. Since the beneficial material properties of microheterogeneous materials originate in the microstructure, its incorporation in material modeling is desired. The FE₂ method is a suitable tool to accomplish this step. Applying this method, a microscopic boundary value problem, given by a representative volume element (RVE) is attached to every integration point of the macroscale. Such RVEs taken from real microstructures exhibit a high level of complexity, thus raise computational costs. The use of a SSRVE with lower complexity is able to decrease these costs. Here, different statistical measures are compared in view of the performance in the SSRVE construction procedure. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Local probabilistic homogenization schemes for assessment of material uncertainties in solid foams

The present study is concerned with a numerical scheme for prediction of the effective properties of solid foams considering their uncertainty. The approach is based on an analysis of a large-scale, statistically representative volume element which is subdivided into small-scale testing volume elements. Application of a standard homogenization scheme to the testing volume elements together with a stochastic evaluation yields a complete probabilistic characterization of the material which may be used for a random field definition of the material behaviour in a macroscopic effective field analysis of foam structures. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Effective parameters of cancellous bone

Cancellous bone is a two–component structure consisting of the bone frame and interstitial blood marrow. In the scope of this presentation, the multiscale finite element method is used for its modeling. This method results from a combination of homogenization theory and the theory of finite elements and is based on the calculation of effective material parameters by investigating representative volume elements (RVEs). For the particular kind of material considered here, a cubic two–phase RVE is assumed where the dry skeleton is modeled in different ways. Apart from the variations of the geometry, the influence of the usage of different types of finite elements is studied in this context. Note that the presence of a liquid phase requires dynamic investigation including the viscous phenomena. To this end, acoustic excitation and an analysis in the complex domain are chosen. The method permits calculation of the effective material parameters such as Young's modulus, bulk modulus and Poisson's ratio and furthermore the simulation of the behaviour of the complete bone or of its parts. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Iterative versus direct parallel substructuring methods in semiconductor device modelling

The numerical simulation of semiconductor devices is extremely demanding in term of computational time because it involves complex embedded numerical schemes. At the kernel of these schemes is the solution of very ill‐conditioned large linear systems. In this paper, we present the various ingredients of some hybrid iterative schemes that play a central role in the robustness of these solvers when they are embedded in other numerical procedures. On a set of two‐dimensional unstructured mixed finite element problems representative of semiconductor simulation, we perform a fair and detailed comparison between parallel iterative and direct linear solution techniques. We show that iterative solvers can be robust enough to solve the very challenging linear systems that arise in those simulations. Copyright © 2004 John Wiley & Sons, Ltd.     

Investigation of microcracks in ferroelectric materials by application of a grain-boundary-motivated cohesive law

Ferroelectric materials exhibit a huge potential for engineering applications – ranging from electrical actuators (inverse piezoelectric effect) to sensor technology (direct piezoelectric effect). To give an example, lead zirconate titanate (PZT) is a typical perovskite ion crystal possessing ferroelectric properties. In this contribution, we are particularly interested in the modelling of microcracking effects in ferroelectric materials. In view of Finite-Element-based simulations, the geometry of a natural grain structure, as observed on the so-called micro-level, is represented by an appropriate mesh. While the response on the grains themselves is approximated by coupled continuum elements, grain boundaries are numerically incorporated via so-called cohesive-type elements. For the sake of simplicity, switching effects in the bulk material will be neglected. The behaviour of the grain boundaries is modelled by means of cohesive-type laws. Identifying grain boundaries as potential failure zones leading to microcracking, cohesive-type elements consequently offer a great potential for numerical simulations. As an advantage, in the case of failure they do not a priori result in ill-conditioned systems of equations as compared with the application of standard continuum elements to localised deformations. Finally, representative constitutive relations for both the bulk material and the grain boundaries, enable two-dimensional studies of low-cycle-fatigue motivated benchmark boundary value problems. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A contribution to the modeling of metal foams on the mesoscale

Due to their useful properties in lightweight construction and due to their excellent behavior in energy absorption for example in crash mechanics, metal foams became an interesting, often utilized and investigated material. For the determination of the mechanical properties of foams without the help of expensive experiments, a way for computing these properties is searched. The problem in doing so is that foams can be composed out of randomly distributed edges and faces with varying thickness and of other inhomogeneities on the mesoscale like imperfections. The goal in this paper is, to investigate the influence of these irregularities on the mechanical, linear elastic properties of a metal foam on the macroscale and to determine the size of a representative volume element, for which the irregularities on the mesoscale do not have a great influence on the linear elastic properties. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Modeling of microscopic failure in heterogeneous materials

The proper modeling of state-of-the-art engineering materials requires a profound understanding of the nonlinear macroscopic material behavior. Especially for heterogeneous materials the effective macroscopic response is amongst others driven by damage effects and the inelastic material behavior of the individual constituents [1]. Since the macroscopic length scale of such materials is significantly larger than the fine-scale structure, a direct modeling of the local structure in a component model is not convenient. Multiscale techniques can be used to predict the effective material behavior. To this end, the authors developed a modeling technique based on representative volume elements (RVE) to predict the effective material behavior on different length scales. The extended finite element method (XFEM) is used to model discontinuities within the material structure independent of the underlying FE mesh. A dual enrichment strategy allows for the combined modeling of kinks (material interfaces) and jumps (cracks) within the displacement field [2]. The gradual degradation of the interface is thereby controlled by a cohesive zone model. In addition to interface failure, a non-local strain driven continuum damage model has been formulated to efficiently detect localization zones within the material phases. An integral formulation introduces a characteristic length scale and assures the convergence of the approach upon mesh refinement [3]. The proposed method allows for an efficient modeling of substantial failure mechanisms within a heterogeneous structure without the need of remeshing or element substitution. Due to the generality of the approach it can be used on different length scales. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Cosserat Parameter Identification within the Frame of the Discrete Element Method

One of the main challenges using the Discrete Element Method is that there is no direct compliance to the well known continuum parameters such as elastic moduli. In this article we show how homogenization procedures using representative volume elements composed of discrete particles lead to Cosserat continua. Simulating a shear test with discrete elements it becomes obvious, that the evolving microstructure is mainly composed of contact chains that form triangles and quadrilaterals. For these contact chains we set up contact energies in normal and shear directions and combine those to derive the effective energy of the material. By comparison of this energy to a Cosserat energy we can derive formulas for the Lamé and Cosserat parameters. They are now only dependent on the interaction energies and radii of the particles. To show the validity of our assumptions and derivations we present some discrete element simulations of shear tests. (© 2012 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)     

A computational two scaled model for the simulation of micro-heterogeneous low-alloyed TRIP steels

In transformation induced plasticity (TRIP) steel a diffusionless austenitic-martensitic phase transformation induced by plastic deformation can be observed, resulting in excellent macroscopic properties. In particular low-alloyed TRIP steels, which can be obtained at lower production costs than high-alloyed TRIP steel, combine this mechanism with a heterogeneous arrangement of different phases at the microscale, namely ferrite, bainite, and retained austenite. The macroscopic behavior is governed by a complex interaction of the phases at the micro-level and the inelastic phase transformation from retained austenite to martensite. A reliable model for low-alloyed TRIP steel should therefore account for these microstructural processes to achieve an accurate macroscopic prediction. To enable this, we focus on a multiscale method often referred to as FE² approach, see [6]. In order to obtain a reasonable representative volume element, a three-dimensional statistically similar representative volume element (SSRVE) [1] can be used. Thereby, also computational costs associated with FE² calculations can be significantly reduced at a comparable prediction quality. The material model used here to capture the above mentioned microstructural phase transformation is based on [3] which was proposed for high alloyed TRIP steels, see also e.g. [8]. Computations based on the proposed two-scale approach are presented here for a three dimensional boundary value problem to show the evolution of phase transformation at the microscale and its effects on the macroscopic properties. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Multiscale Modelling of the Effective Viscoplastic Material Behaviour of Textile-Reinforced Polymers Using XFEM

In this contribution a modelling approach using numerical homogenisation techniques is applied to predict the effective nonlinear material behaviour of composites from simulations of a representative volume element (RVE). Numerical models of the heterogeneous material structure in the RVE are generated using the eXtended Finite Element Method (XFEM) which allows for a regular mesh. Suitable constitutive relations account for the material behaviour of the constituents. The influence of the nonlinear matrix material behaviour on the composite is studied in a physically nonlinear FE simulation of the local material behaviour in the RVE ­ effective stress-strain curves are computed and compared to experimental observations. The approach is currently augmented by a damage model for the fibre bundle. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A phase field model for martensitic transformations with a temperature-dependent separation potential

Metallic materials often exhibit a complex microstructure with varying material properties in the different phases. Of major importance in mechanical engineering is the evolution of the austenitic and martensitic phases in steel. The martensitic transformation can be induced by heat treatment or by plastic surface deformation at low temperatures. A two dimensional elastic phase field model for martensitic transformations considering several martensitic orientation variants to simulate the phase change at the surface is introduced in [1]. However here, only one martensitic orientation variant is considered for the sake of simplicity. The separation potential is temperature dependent. Therefore, the coefficients of the Landau polynomial are identified by results of molecular dynamics (MD) simulations for pure iron [1]. The resulting separation potential is applied to analyse the mean interface velocity with respect to temperature and load. The interface velocity is computed by use of the dissipative part to the configurational forces balance as suggested in [3]. The model is implemented in the finite element code FEAP using standard 4-node elements with bi-linear shape functions. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Micro-mechanical studies on the effect of various stress-states on ductile damage and failure

The paper deals with the effect of different stress states on damage and failure behavior of ductile materials. To be able to model these effects a continuum damage model has been proposed taking into account the dependence of the stress intensity, the stress triaxiality and the Lode parameter on the constitutive equations. The model is based on the introduction of damaged and fictitious undamaged configurations. Only experiments are not adequate enough to determine all constitutive parameters. Therefore, additional three-dimensional micro-mechanical simulations of representative volume elements have been performed to get more insight in the complex damage mechanisms. These simulations cover a wide range of different void sizes, void shapes and void distributions. After all, the results from the micro-mechanical simulations are used to propose the damage equations and to identify corresponding parameters. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Computational Homogenization of Micromechanically Resolved Textile Materials

This contribution deals with textile materials. On the macroscopic level textiles are characterized by a large area-to-thickness ratio, such that it is numerically efficient to treat the textile structure as a shell. To capture the contact behavior, fibers within a representative volume element are explicitly modeled by means of one dimensional beam elements on the microscopic level. A suitable, shell-specific homogenization method is developed, which connects the homogeneous shell specific macro level to a fiber structured micro level. This contribution investigates the determination of the nonlinear constitutive behavior of textile materials. Selected examples for the macroscopic behavior of microscopic heterogeneous fiber structured materials are presented. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Computation of the effective nonlinear material behaviour of composites using X-FEM - Analysis of the matrix material behaviour

For a consequent 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, the modelling approach using numerical homogenization techniques based on the simulation of representative volume elements which are modelled by the extended finite element method (X–FEM) is currently extended to nonlinear material behaviour. While the glass fibres are assumed to remain linear elastic, a viscoplastic constitutive law accounts for strain–rate dependence and inelastic deformation of the matrix material. This paper describes the process of finding suitable constitutive relations for the polymeric matrix material Polypropylene in the small–strain regime. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)