On the simulation of textile reinforced composites and structures

This paper presents experimental and numerical methods to perform simulations of the mechanical behavior of textile reinforced composites and structures. The first aspect considered refers to the meso-to-macro transition in the framework of the finite element (FE) method. Regarding an effective modelling strategy the Binary Model is used to represent the discretized complex architecture of the composite. To simulate the local response and to compute the macroscopic stress and stiffness undergoing small strain a user routine is developed. The results are transfered to the macroscopic model during the solution process. The second aspect concerns the configuration of the fiber orientation and textile shear deformation in complex structural components. To take these deformations which affect the macroscopic material properties into account they are regarded in a macroscopic FE model. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On the simulation of textile reinforced concrete using a multi scale approach

Textile reinforced concrete (TRC) is a composite of textile structures made of multi-filament yarns (rovings) within a cementitious matrix. Experimental investigations of textile reinforced concrete specimen show very complex failure mechanisms on different length scales. Therefore mechanical models on the micro, meso and macro scale are introduced. The paper presents a hierarchical material model of TRC on three scales. While on the micro scale the individual filaments of the fiber bundles are distinguished to determine an effective roving behavior, models on the meso scale are used to predict the macroscopic response of the composite material. (© 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)     

Extended finite element modeling of deformable porous media with arbitrary interfaces

In this paper, an enriched finite element method is presented for numerical simulation of saturated porous media. The arbitrary discontinuities, such as material interfaces, are encountered via the extended finite element method (X-FEM) by enhancing the standard FEM displacements. The X-FEM technique is applied to the governing equations of porous media for the spatial discretization, followed by a generalized Newmark scheme used for the time domain discretization. In X-FEM, the material interfaces are represented independently of element boundaries and the process is accomplished by partitioning the domain with some triangular sub-elements whose Gauss points are used for integration of the domain of elements. Finally, several numerical examples are analyzed, including the dynamic analysis of the failure of lower San Fernando dam, to demonstrate the efficiency of the X-FEM technique in saturated porous soils.     

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)     

Finite element (FE) simulation of textile-reinforced composites using a COSSERAT Element

The efficient simulation of textile reinforced composites requires reliable material parameters describing the macroscopic material behavior. Due to the micro–structure, decoupled tensile and bending stiffnesses are observed in experimental investigations. This motivates the use of higher order continuum theories. For numerical simulation of this material behavior, the formulation and application of a volume element based on the Cosserat continuum theory is presented. (© 2008 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)     

Modeling of the effective viscoelastic material behavior of textile reinforced composites using a multi-scale approach

The increasing importance of constructive lightweight in modern engineering science involves the use of advanced materials like textile reinforced composites. In order to reduce development costs, efficient numerical simulations are needed to model the macroscopic behavior of the final product. Focussing on long term phenomena, which are important when parts made of composites with rate-dependent material behavior are assembled by bolted or screwed joints, a two-step homogenization procedure is used to obtain an effective homogeneous equivalent material at the macroscopic scale. (© 2014 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)     

On the macroscopic material behaviour of textile reinforced concrete undertensile loading

This paper concerns with the development of the macroscopic material behaviour of textile reinforced concrete (TRC) using an analytical approach. Therefore the heterogeneous structure of TRC is modelled on the mesoscopic level. The overall material behaviour on the macroscopic level is obtained by means of the homogenisation technique. The analytical approach is based on the micro mechanical solution for a single inclusion according to Eshelby . In extension of this solution for multidirectional reinforced concrete an effective field approximation is used. This approach considers the interactions between the different orientated rovings and the micro cracks in an average sense. For the mechanical modelling of the bond behaviour between roving and matrix after initiating of the macro cracking a slip based bond model with a multiple linear shear stress-slip relation is used. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A state-of-the-art review of the X-FEM for computational fracture mechanics

This paper presents a review of the extended finite element method X-FEM for computational fracture mechanics. The work is dedicated to discussing the basic ideas and formulation for the newly developed X-FEM method. The advantage of the method is that the element topology need not conform to the surfaces of the cracks. Moreover, X-FEM coupled with LSM makes possible the accurate solution of engineering problems in complex domains, which may be practically impossible to solve using the standard finite element method.     

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)     

A constitutive model for shape memory alloy fibres and its applicability to prestressed composites

This contribution is concerned with a constitutive model for shape memory fibres. The 1D-constitutive model accounts for the pseudoplastic and shape memory effect (SME). The macroscopic answer of the material is determined by the evolution from a twinned martensitic lattice into a deformed and detwinned one. On the macroscopic scale these effects are responsible for the upper boundary of the hysteresis which is situated around the origin of the stress-strain-diagram. During the phase transition process inelastic strains arise. When the lattice is fully detwinned, a linear elastic branch at the end of the hysteresis is observed. The initial state of the material is recovered by unloading and heating the material subsequently. The constitutive model is derived from the Helmholtz' free energy and fulfils the 2nd law of thermodynamics. For the present model five internal state variables are employed. Two of them are used to describe the inelastic strain and a backstress. The others represent the martensitic volume fraction and are necessary to describe the SME. The latter variables are depending on the deformation state as well as on temperature. A change on temperature goes along with a reduction of the inelastic strain. The model is incorporated in a fibre matrix discretization to prestress the surrounding structure. The boundary value problem is solved for a truss element applying the finite element method. Examples will demonstrate the applicability in engineering structures. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A micromechanics model for FRC composites

A multiscale model for FRC composite structures taking into consideration the complex interactions at the scales of the fiber and microcracks is proposed. At the scale of the single fiber, a semi-analytical model characterizes the microslip behavior at the interface between the matrix and the fiber in terms of the overall composite stresses. The influence of fiber bundles on microcrack bridging and arrest is taken into account within the framework of linear elastic fracture mechanics. Upscaling to the macroscopic level using continuum micromechanics shows that the macroscopic deformation of the FRC composite is governed by a ’TERZAGHI’ like effective stress. For the finite element analyses of failure behavior at the scale of the composite structure, an ’interface solid element’ technique is used to consider localized cracking. Selected numerical and semi-analytical results together with experimental validations are provided. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

About the Microstructural Effects of Polycrystalline Materials and their Macroscopic Representation at Finite Deformation

In the sheet bulk metal forming field, the strict geometrical requirements of the workpieces result in a need of a precise prediction of the material behaviour. The simulation of such forming processes requires a valid material model, performing well for a huge variety of different geometrical characteristics and finite deformation. Because of the crystalline nature of metals, anisotropies have to be taken into account. Macroscopically observable plastic deformation is traced back to dislocations within considered slip systems in the crystals causing plastic anisotropy on the microscopic and the macroscopic level. A finite crystal plasticity model is used to model polycrystalline materials in representative volume elements (RVEs) of the microstructure. A multiplicative decomposition of the deformation gradient into elastic and plastic parts is performed, as well as a volumetric-deviatoric split of the elastic contribution. In order to circumvent singularities stemming from the linear dependency of the slip system vectors, a viscoplastic power-law is introduced providing the evolution of the plastic slips and slip resistances. The model is validated with experimental microstructural data under deformation. The validation on the macroscopic scale is performed through the reproduction of the experimentally calculated initial yield surface. Additionally, homogenised stress-strain curves from the microstructure build the outcome for a suitable effective material model. Through optimisation techniques, effective material parameters can be determined and compared to results from real forming processes. (© 2012 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)     

Modeling of bainitic phase transformation

In this work, we present a macroscopic material model for simulation of austenite to bainite and of austenite to martensite transformations accompanied by transformation-induced plasticity (TRIP), which is an important phenomenon in metal forming processes. In order to account for the incubation time the model considers nucleation of the bainite phase. When this quantity attains a barrier term, growth of bainite volume fraction is started. The model formulation allows for individual evolutions of upper and lower bainite. Furthermore, the numerical implementation of the constitutive equations into a finite element program is described. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical Simulation of Crack Growth

In recent years the X-FEM based on the partition of unity method and the strong discontinuity approach (SDA) have shown to be powerful tools to model crack growth. Both methods model the crack surface by introducing additional d.o.f.. In the X-FEM the nodes in the mesh around a crack are globally enhanced with new d.o.f. while in in the SDA the new d.o.f. are commonly introduced as internal ones. Thus the jump displacement fields are constant across elements. Therefore the d.o.f. can be condensed on element level which results in jumps in the displacement field at element edges. In this contribution the strong discontinuity approach is used approximating the displacement jump linearly across the crack length similar as e.g. in [3]. New additional nodes of the cracked elements that lie on the element edges are introduced but are not considered as internal nodes but remain global. Thus crack path continuity is automatically given. These global d.o.f. approximate the discontinuous part of the displacement field. The sum of the aforementioned part and the continuous displacement field represent the total displacement field including a possible jump. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)