Modeling the constitutive behavior of ferroelectric,ferromagnetic and multiferroic materials by using the condensed method (CM)

Due to the multifunctional applicability, smart materials are of particular interest in the field of material modeling. Most of the developed models, describing the nonlinear behavior, are implemented within the framework of the Finite Element Method (FEM). However, most investigations are restricted to simple boundary value problems (BVP) under uniaxial loading and their goal is the calculation of hysteresis loops. Regarding this circumstance, the so-called condensed method (CM) is introduced to investigate the macroscopic polycrystalline ferroelectric material behavior at a global material point without any kind of discretization scheme. In the presented paper, the CM is extended towards ferromagnetic and multiferroic material behavior. Moreover, numerical results for a pure ferromagnetic behavior and a comparison between the magnetoelectric coupling coefficient calculated by the FEM and the CM are presented. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Linking macroscopic deformation processes to microstructure evolution using an FE-FFT-based micro-macro transition and non-conserved phase-fields

The purpose of this work is the multiscale FE-FFT-based prediction of macroscopic material behavior, micromechanical fields and bulk microstructure evolution in polycrystalline materials subjected to macroscopic mechanical loading. The macroscopic boundary value problem (BVP) is solved using implicit finite element (FE) methods. In each macroscopic integration point, the microscopic BVP is embedded, the solution of which is found employing fast Fourier transform (FFT), fixed-point and Green's function methods. The mean material response is determined by the stress-strain relation at the micro scale or rather the volume average of the micromechanical fields. The evolution of the microstructure is modeled by means of non-conserved phase-fields. As an example, the proposed methodology is applied to the modeling of stress-induced martensitic phase transformations in metal alloys. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Constitutive modeling of ferromagnetic and ferroelectric behaviors and application to multiferroic composites

In this paper, the constitutive modeling of nonlinear multifield behavior as well as the finite element implementation are presented. Nonlinear material models describing the magneto-ferroelectric or electro-ferromagnetic behaviors are presented. Both physically and phenomenologically motivated constitutive models have been developed for the numerical calculation of principally different nonlinear magnetostrictive behaviors. Further, the nonlinear ferroelectric behavior is based on a physically motivated constitutive model. On this basis, the polarization in the ferroelectric and magnetization in the ferromagnetic and magnetostrictive phases, respectively, are simulated and the resulting effects analyzed. Numerical simulations focus on the calculation of magnetoelectric coupling and on the prediction of local domain orientations going along with the poling process, thus supplying information on favorable electric-magnetic loading sequences. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Aspects of a direct homogenization procedure for electro-mechanically coupled problems

The aim of this work is to discuss a micro–macro homogenization procedure for electro–mechanically coupled problems. In this context a two–scale homogenization ansatz for ferroelectric ceramics based on an FE²-approach is presented. The microscopic discretization of the heterogeneous structure of the polycrystalline material allows for the incorporation of microscopic effects, which are necessary to determine the corresponding overall macroscopic material response. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Nonlinear numerical simulation of ferroelectric-ferromagnetic multifunctional composites

The theoretical background of nonlinear constitutive multifield behavior is presented. Nonlinear material models describing the ferroelectric or ferromagnetic behaviors are presented. Both physically and phenomenologically motivated constitutive models have been developed for the numerical calculation of the nonlinear magnetostrictive and ferroelectric behaviors. On this basis, the polarization in the ferroelectric and magnetization in the ferromagnetic respectively magnetostrictive phases are simulated and the resulting effects analyzed. The developed tools enable the prediction of the electromagnetomechanical properties of smart multiferroic composites and supply useful means for their optimization. Goals are to improve the efficiency of ME coupling and to reduce damage associated with the poling process. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A phenomenological constitutive model for magnetostrictive materials and ferroelectric ceramics

The contribution is concerned with a thermodynamic consistent constitutive model for magnetostrictive materials and ferroelectric ceramics. It captures the nonlinear phenomenological behavior which is described by hysteresis effects. Magnetostrictive alloys and ferroelectric ceramics belong to the multifunctional materials. In recent years these materials have become widely–used in actor and sensor applications. They characterize an inherent coupling between deformation and magnetic or electric field. Due to the similarities of the coupled differential equations a uniform approach is applied for both phenomena. The presented three–dimensional material model is thermodynamically motivated. It is based on the definition of a specific free energy function and a switching criterion. Furthermore an additive split of strain and the magnetic or electric field in a reversible and an irreversible part is suggested. The irreversible quantities serve as internal variables, which is analog to plasticity theory. A one–to–one–relation between the two internal variables provides conservation of volume for the irreversible strains. The presented material model can approximate the ferromagnetic or ferroelectric hysteresis curve and the related butterfly hysteresis. Furthermore an extended approach for ferrimagnetic behavior, which occurs in magnetostrictive materials, is presented. Some numerical simulations demonstrate the capability of the presented model. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Modelling of cracks on different scales in ferroelectric materials

In this work we study contributions to the effective fracture toughness of ferroelectric materials arising from effects on macroscopic and mesoscopic scales of the system. On the macroscopic scale, the crack in a ferroelectric material is modeled taking into account an extended theory of stresses at interfaces in dielectric solids [1-3]. We predict several new effects, such as the “poling effect”, “collinear effect” and the coupling of a Mode-II shear loading and the Mode-I SIF. Further, on the mesoscopic scale, we study the influence of polarization switching limited to the fracture process zone (small scale switching) on the fracture toughness. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A hybrid element formulation for electromechanical problems

The difficulty in the modeling of ferroelectric materials is the coverage of the complicated interactions between electrical and mechanical quantities on the macroscale, which are caused by switching processes on the microscale. In the present work we present an electric hybrid element formulation where the stresses and the electric fields are derived by constitutive relations as presented in [1]. Therefore the displacements, the electric potential and the electric displacements are approximated by bilinear ansatz functions. Applying a static condensation procedure we obtain a modified finite element formulation governed by the degrees of freedoms associated to the displacements and the electric potential. The anisotropic material behavior is modeled within a coordinate-invariant formulation [6] for an assumed transversely isotropic material [4]. In this context a general return algorithm is applied to compute the remanent quantities at the actual timestep. Resulting hysteresis loops for the ferroelectric ceramics are presented. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Phase field simulations in ferroelectric materials

The hindering of domain wall movement by defects in ferroelectric materials is closely connected to electric fatigue. A movable domain wall in a ferroelectric material in most cases is modelled as a singular surface which allows the use of configurational forces. In contrast, the present approach treats the polarization as an order parameter, extending the total energy by a phase separation energy and a domain wall energy. The polarization then no longer has a discontinuity at the domain wall but is a continuous vector field (phase field). As an example, a numerical simulation of domain evolution under stress free boundary conditions is presented. (© 2006 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)     

A finite element framework for the numerical implementation of boundary potentials

This contribution deals with the implications of boundary potential energies on deformational mechanics in the framework of the finite element method at finite strains. The common material models in continuum mechanics are taking the bulk into account, nevertheless, neglecting the boundary. However, boundary effects sometimes play a dominant role in the material behavior, e.g. surface tension in fluids. The boundary potentials, in general, are allowed to depend not only on the boundary deformation gradient but also on the spatial surface–normal / curve–tangent, as well. For the finite element implementation, a suitable curvilinear coordinate system attached to the boundary is defined and corresponding geometrical and kinematical derivations are carried out. Afterwards, the discretization of the generalized weak formulation, including boundary potentials, is carried out and finally numerical examples are presented to demonstrate the boundary effects due to the different proposed material behavior. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A hyperelastic biphasic fiber reinforced model for articular cartilage considering the distribution and orientation of collagen fibers

Articular cartilage is a multiphase material consisting of fluids and electrolytes, which is described with the Theory of Porous Media. The mechanical characteristics of articular cartilage are porosity, incompressible material behavior combined with transversely isotropic behavior for solid and fluid phases. There are two central points to model articular cartilage: the poro-viscosity of the porous matrix and the visco elasticity, and orientation of the collagen fibers. A numerical example is presented. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Laminate-based modelling of microstructure and switching in ferroelectrics

The purpose of the work is the thermodynamics-based modelling of the polarisation and the deformation microstructure in the ferroelectric single crystal with the help of a laminate-based approach. The incremental variational-based rate-dependent macroscopic model for dissipative ferroelectric material [1] and the laminate-based microscopic model [2] established in the literature are taken as basis and shall be further extended to a single crystal laminate structure dependent on the loading frequency based on the coupled electromechanical framework taking the effect of polarisation into account. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Computational characterization of ferroelectric solids with a 3D Preisach model based on an orientation distribution function

Ferroelectric or ferromagnetic materials show an interaction between mechanical deformations and polarization or magnetization. A few multiferroic materials possess both ferroic properties and exhibit a magneto-electric (ME) coupling. These ME properties can be achieved in two-phase composites, which combine ferroelectric and ferromagnetic characteristics. To predict a realistic material behavior and a more precise ME coefficient, the application of suitable material models which describe the nonlinear hysteretic behavior is of particular importance. In the present contribution we focus on the characterization of a nonlinear ferroelectric material behavior, in terms of a 3D Preisach model based on an orientation distribution function. (© 2016 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)     

Study on statistically similar RVEs for real microstructures based on different statistical descriptors

A method for constructing statistically similar representative volume elements (SSRVEs) for a real dual-phase (DP) steel microstructure is presented in this contribution. The advantageous material properties of such kind of steels originate from the interaction of the microstructure constituents of the material on the microscale. In order to capture these effects directly in the material modeling, the FE² method is a suitable tool, where an RVE representing the microstructure of a material is used in the microscopic boundary value problem, which is solved at each Gauss-point of the macroscopic boundary value problem. However, RVEs based on real microstructures typically implicate high computational expenses due to the complexity of the underlying microstructure and its discretization. SSRVEs, which have a lower complexity than conventional RVEs but are still able to represent the material, can be used instead. Here, different statistical measures for the construction of SSRVEs and their comparison are the main focus. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Comparison of numerical homogenization techniques in BT-CFO multiferroic particle composites

This work is restricted to linear material behavior, i.e. the structure is considered to be in a perfectly poled state. Different numerical homogenization methods are investigated and used to calculate effective properties of a 0-2 composite modelled in reprensentative volume elements. Bariumtitanate (BT) and cobaltferrite (CFO) are employed in the Finite Element model, where the roles of matrix and inclusion are mutable in principle. Mixed magnetoelectromechanical boundary conditions based on different homogenization theories are applied to the model. The calculated macroscopic behaviors described by the different approaches are compared and presented in the paper. The special focus is on the prediction of coefficients of magnetoelectric coupling with respect to an optimization of the structural arrangement of the composite. (© 2014 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)