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)     

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)     

Effective elastic properties of particulate-reinforced compossite materials

A numerical approach for determination of the effective properties of particulate composite materials has been developed. A representative volume element (RVE) of the composite material is analyzed with help of the finite-element method. Uniform boundary displacements or tractions are applied on the boundaries of the RVE for introducing the known average strain in the RVE. Local stress and strain distributions in the RVE are calculated using the finite-element method. Different effective elastic constants can be calculated by averaging the local fields corresponding to different sets of boundary conditions. The present approach allows us to determine the effective properties of particle-reinforced composites with acceptable accuracy. The calculated effective properties of the composite are between the upper and lower Hashin—Shtrikman bounds. The results based on the present approach lead to higher stiffness of composites in comparison with analytical approaches.Institute fur Werkstoffwissenschaften, Fachberech Werkseoffwissenschaften, Martin-Luther-Universität Halle-Wittenberg, D-06099 Halle, Germany. Published in Mekhanika Kompozitnykh Materialov, Vol. 33, No. 4, pp. 450–459, July–August, 1997.     

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)     

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)     

A RVE-based micro mechanical model for short fiber reinforced plastic materials including matrix damage and fiber fracture

The representative volume element (RVE) method is applied to a fiber reinforced polymer material undergoing matrix damage and fiber fracture. Results of RVE computations are compared to uniaxial tensile tests performed with the composite material. It is shown that the macroscopic behavior of the composite material can accurately be predicted by RVE computations. (© 2011 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)     

The effective heat conduction of short fibre reinforced materials considering the fibre distribution

This contribution focuses on the effective heat conductivity of short fibre reinforced materials. For this purpose, a representative volume element (RVE), which is able to represent all possible fibre orientation distributions, is introduced and modelled in ABAQUS. Subsequently, the effective heat conductivity of the RVE is derived, employing a numerical homogenisation scheme, and a phenomenological material model is fitted to those results. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

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)     

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)     

Micromechanical Approach of Martensitic Transformation Problem in Steels

To determine the mechanical behavior of material involving the martensitic phase transformation (for example, steels like 100Cr6), a representative volume element (RVE) model including phase transformation criterion is desireable at micromechanical approach. A framework combining the Eshelby's inclusion theory as well as continuum mechanics with phase-transformation (PT) critical condition at RVE model is presented briefly. And application of this model to estimate the critical aspect ratio of martensitic plate or lath inside homogeneneous stress field is also included, where the RVE can be under uniaxial tension/compression or pure shear loading case. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Multiscale modeling of continuous fiber and fabric reinforced rubber components

Fabric or continuous fiber reinforced rubber components (e.g. tires, air springs, industrial hoses, conveyor belts or membranes) are underlying high deformations in application and show a complex, nonlinear material behavior. A particular challenge depicts the simulation of these composites. In this contribution we show the identification of the stress and strain distributions by using an uncoupled multiscale modeling method, see [1]. Within this method, two representation levels are described: One, the meso level, where all constituents of the composite are shown in a discrete manner by a representative volume element (RVE) and secondly, the macro level, where the structural behavior of the component is defined by a smeared anisotropic hyperelastic constitutive law. Uncoupled means that the RVE does not drive the macroscopic material behavior directly as in a coupled approach, where a RVE boundary value problem has to be solved at every integration point of the macro level. Thus an uncoupled approach leads to a tremendous reduction in numerical effort because the boundary value problem of a RVE just has to be solved at a point of interest, see [1]. However, the uncoupled scale transition has to fulfill the HILL–MANDEL condition of energetic equivalence of both scales. We show the calibration of material parameters for a given constitutive model for fiber reinforced rubber by fitting experimental data on the macro level. Additionally, we demonstrate the determination of effective properties of the yarns. Finally, we compare the energies of both scales in terms of compliance with the HILL–MANDEL condition by using the example of a biaxial loaded sample and discuss the consequences for the mesoscopic level. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Analysis of an aggregation‐based algebraic two‐grid method for a rotated anisotropic diffusion problem

A two‐grid convergence analysis based on the paper [Algebraic analysis of aggregation‐based multigrid, by A. Napov and Y. Notay, Numer. Lin. Alg. Appl. 18 (2011), pp. 539–564] is derived for various aggregation schemes applied to a finite element discretization of a rotated anisotropic diffusion equation. As expected, it is shown that the best aggregation scheme is one in which aggregates are aligned with the anisotropy. In practice, however, this is not what automatic aggregation procedures do. We suggest approaches for determining appropriate aggregates based on eigenvectors associated with small eigenvalues of a block splitting matrix or based on minimizing a quantity related to the spectral radius of the iteration matrix. Copyright © 2015 John Wiley & Sons, Ltd.     

Modelling additive manufactured materials using a crystal plasticity model

In this contribution the macroscopic behaviour of additively manufactured Inconel 718 will be identified using a crystal plasticity model on the meso-level. The experimentally observed grain structures will be used to generate a representative volume element. Based on this RVE macroscopic mechanical parameters will be identified and compared with experimental results. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

界面弧形裂纹对混凝土开裂强度的影响研究

混凝土由于水分蒸发、干缩、泌水以及骨料与砂浆变形不一致等原因会导致骨料与砂浆的界面层中产生弧形裂纹,从而对混凝土开裂强度产生很大影响.从细观角度将混凝土视作由粗骨料和水泥砂浆组成的两相复合材料,并将界面层视为粗骨料与水泥砂浆的接触层进行分析.首先基于相互作用直推估计(interaction direct derivative, IDD)法,考虑混凝土中骨料颗粒的相互作用,将施加在混凝土表征体积元的远场外荷载等效为无限大基体中含单一骨料的等效外荷载.然后,将等效外荷载转化为最大和最小主应力,基于断裂力学理论得到界面层中弧形裂纹的应力强度因子,并根据复合型裂纹幂准则判断弧形裂纹是否发生开裂,进而来研究混凝土开裂强度的变化规律.通过与数值模拟结果的比较,验证了界面弧形裂纹应力强度因子解析解的有效性,参数分析结果表明,当裂纹与最大主应力垂直或与最小主应力呈45°夹角时,骨料周围弧形裂纹最易发生开裂破坏.随着裂纹长度增加,混凝土受拉和受压开裂强度先减小后增大,且均存在最不利的裂纹长度.混凝土开裂强度随着骨料体积分数的增加而增大,随着骨料粒径的增大而减小.在裂纹长度较小时,增大骨料的弹性模量有利于提高混凝土开裂强度.骨料周围承受同号应力可以提高混凝土的开裂强度,反之,异号应力会降低开裂强度.     

Numerical homogenization of composite structures with rhombic fiber arrangements

A numerical procedure is developed to determine effective material properties of unidirectional fiber reinforced composites with rhombic fiber arrangements. With the assumption of a periodic micro structure a representative volume element (RVE) is considered, where the phases have isotropic or transversely isotropic material characterizations. The interface between the phases is treated as perfect. The procedure handles the primary non-rectangular periodicity with homogenization techniques based on finite element models. Due to appropriate boundary conditions applied to the RVE elastic effective coefficients are derived. Six different boundary condition states are required to get all coefficients of the stiffness tensor. Results are listed and compared with other publications and good agreements are shown. Furthermore new results are presented, which exhibit the orthotropic behavior of such composites caused by the rhombic fiber arrangement. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Spatial Random Material Property Model for Vibration of Composites

Investigation of vibration and buckling of thin walled composite structures is very sensitive to parameters like uncertain material properties and thickness imperfections. Because of the manufacturing process and others, thin walled composite and other structures show uncertainties in material properties, and other parameters which cannot be reduced by refined discretization. These parameters are mostly spatial distributed in nature. Here I introduce a semivariogram type material property model to predict the spatial distributed material property (like young's modulus) over the structure. The computation of semivariogram parameters needs the local material properties over a prespecified gird. The material properties at each grid have been obtained by considering a statistically homogeneous representative volume element (RVE) at each gird. According to random nature of the spatial arrangement of fibers, the statistically homogeneous RVE is obtained using image processing. The effective material properties of the RVE have been obtained numerically with the help of periodic boundary condition. The methodology is applied to a composite panel model and modal analysis has been carried. The results of the modal analysis (eigen values and mode shapes) are compared with experimental modal analysis results which are in good agreement. Using the presented material property model we can better predict the vibration characteristics of the thin walled composite structures with the inherent uncertainties. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Multiscale simulations of concrete mechanical tests

In civil engineering, computational modeling is widely used in the design process at the structural level. In contrast to that, an automated support for the selection or design of construction materials is currently not available. Specification of material properties and model parameters has a strong influence on the results. Therefore, an uncoupled two-step approach is employed to provide relatively quick and reliable simulations of concrete (mortar) tests. First, the Mori–Tanaka method is utilized to include the majority of small aggregates and air voids. The strain incremental form of MT approach serves for the prediction of material properties subsequently used in the finite element simulations of mechanical tests.     

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)     

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)