A 1D constitutive law for hysteresis effects in piezoceramics

This paper is concerned with a macroscopic constitutive law for domain switching effects, which occur in piezoelectric ceramics. The thermodynamical framework of the law is based on two scalar valued functions: the Helmholtz free energy and a switching surface. In common usage, the remanent polarization and the remanent strain are employed as internal variables. The novel aspect of the present work is to introduce an irreversible electric field, which serves besides the irreversible strain as internal variable. The irreversible electric field has only theoretical meaning, but leads to advantages within the finite element implementation, where displacement and the electric potential are the nodal degrees of freedoms. A common assumption is a one-to-one relation between the irreversible strain and the polarization. This simplification is not employed in the present paper. To accomplish enough space for the polarization, resulting from an applied electric field, the irreversible strains are additively split and a special hardening function is introduced. This balances the amount of space and the domain switching due to polarization. The constitutive model reproduces the ferroelastic and the ferroelectric hysteresis as well as the butterfly hysteresis for piezoelectric ceramics and it accounts for the mechanical depolarization effect. (© 2006 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)     

Pseudoelasticity of Shape Memory Alloys - a 1D Material Model and Finite Element Implementation

This contribution is concerned with the formulation of a 1D-constitutive model accounting for the pseudoelastic behavior of shape memory alloys. The stress-strain-relationship is idealized by a hysteresis both in the compression as in the tension loading range. It is characterized by an upper loading path, which is to be ascribed to the transformation of the lattice to a martensitic structure. Unloading the material, a lower path is described, because of the reverse transformation into austenitic lattice. The constitutive model is based on a switching criterion which serves as a potential function for the evolution of the internal state variables. The model distinguishes between local and global variables to describe the hysteresis effects for the compression and tension range. A strain driven algorithm which captures the complete nonlinear material behavior is presented. The boundary value problem is solved for a truss element applying the finite element method. A consistent linearization of the nonlinear equations is derived. Simple examples will demonstrate the applicability of the proposed model. For future developments the usage of shape memory alloys within civil engineering structures is aimed. The advantage of the material is the very good damping behavior and the potential to overcome great strains. Both properties are distinguished to be of engineering interest. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Full 3D finite element Cosserat formulation with application in layered structures

In this paper, a full three-dimensional (3D) finite element Cosserat formulation is developed within the principles of continuum mechanics in the small deformation framework. The developed finite element formulation is general; however, the proposed constitutive laws incorporate the effect of the internal length parameter of 3D layered continua. The extension of the existing two-dimensional (2D) Cosserat formulation to the 3D framework is novel and is consistent with plate theory which can be considered as the 3D version of beam theory. The results demonstrate a high level of consistency with the analytical solutions predicted by plate theory as well as predictions by alternative numerical techniques such as the discrete element method.     

Modeling and Simulation of Rate-Dependent Magneto-Active Polymers

In the present work, the magneto-viscoelastic behavior of MAPs is studied by a thermodynamically consistent constitutive model. A finite deformation based framework of nonlinear magneto-viscoelastic coupling is introduced with a multiplicative decomposition of the deformation gradient. The viscosity is captured by evolution equations of the internal variables introduced. We propose energy functions for pure magnetic and magneto-mechanical coupling such that saturation behavior of the magnetostriction and magnetization is captured. After having established the general framework, the model is studied for homogeneous deformations for the purpose of a least-square-based parameter identification from experimental data. The model predictions of non-linear magneto-mechanical responses with strong rate and field dependency are presented. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Mathematical modelling of finite strain magneto-viscoelastic deformations

Large strain magneto-viscoelastic deformations in the presence of a finite magnetic field are modelled in this paper. Internal dissipation mechanisms are proposed using a multiplicative decomposition of the deformation gradient and an additive decomposition of the magnetic induction. Using thermodynamically consistent constitutive and evolution laws, numerical results showing stress relaxation and magnetic field evolution are presented to illustrate the theory. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Modeling and computation of metal powder compaction processes

In this short communication the numerical treatment of a thermomechanical consistent finite strain viscoplasticity model for metal powder compaction is discussed. The convex single surface yield function evolves according to two evolution equations and remains convex under all loading conditions. The very challenging numerical treatment on local level for integrating the constitutive model requires particular globally convergent Newton-like method with inequality constraints so that a stable solution scheme results. This is embedded into a time-adaptive finite element program which makes use of diagonally-implicit Runge-Kutta methods combined with a Multilevel-Newton algorithm. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Theory and Numerics of Rate-Dependent Incremental Variational Formulations in Ferroelectricity

This paper is concerned with macroscopic continuous and discrete variational formulations for domain switching effects at small strains, which occur in ferroelectric ceramics. The developed new three–dimensional model is thermodynamically–consistent and determined by two scalar–valued functions: the energy storage function (Helmholtz free energy) and the dissipation function, which is in particular rate–dependent. The constitutive model successfully reproduces the ferroelastic and the ferroelectric hysteresis as well as the butterfly hysteresis for ferroelectric ceramics. The rate–dependent character of the dissipation function allows us also to reproduce the experimentally observed rate dependency of the above mentioned hysteresis phenomena. An important aspect is the numerical implementation of the coupled problem. The discretization of the two–field problem appears, as a consequence of the proposed incremental variational principle, in a symmetric format. The performance of the proposed methods is demonstrated by means of a benchmark problem. (© 2008 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 gradient model for finite strain elastoplasticity coupled with damage

This paper describes the formulation of an implicit gradient damage model for finite strain elastoplasticity problems including strain softening. The strain softening behavior is modeled through a variant of Lemaitre's damage evolution law. The resulting constitutive equations are intimately coupled with the finite element formulation, in contrast with standard local material models. A 3D finite element including enhanced strains is used with this material model and coupling peculiarities are fully described. The proposed formulation results in an element which possesses spatial position variables, nonlocal damage variables and also enhanced strain variables. Emphasis is put on the exact consistent linearization of the arising discretized equations.

A numerical set of examples comparing the results of local and the gradient formulations relative to the mesh size influence is presented and some examples comparing results from other authors are also presented, illustrating the capabilities of the present proposal.     

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)     

A thermodynamic consistent model for an assumed transversely isotropic piezoceramic with nonlinear ferroelectric hysteresis

A characteristic feature of ferroelectric crystals is the appearance of a spontaneous polarisation, where its direction can be reversed by an applied electric field. This quantity, that has a maximum value at high electric‐fields, depends on the loading history of the material. In this paper we discuss a thermodynamic consistent phenomenological model for an assumed transversely isotropic ferroelectric crystal, where the history dependency is modelled by internal variables. The anisotropic behaviour is governed by isotropic tensor functions, depending on a finite set of invariants, that satisfy automatically the symmetry relationships of the considered body. The main goal of this investigation is to capture some characteristics of nonlinear ferroelectrica, such as the polarisation‐electric‐field and the strain‐electric‐field (butterfly) hysteresis loops. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Hysteresis loops for para–ferromagnetic phase transitions

In this paper, by an extension of the Ginzburg–Landau theory, we propose a mathematical model describing hard magnets within which we are able to explore the para–ferromagnetic transition and by using the Landau–Lifshitz–Gilbert equation, to study the 3D evolution of magnetic field. Finally, the hysteresis loops are obtained and represented by numerical implementations.     

A stereographic projection based phase field model for ferromagnetics

By using the stereographic projection to satisfy the constant magnetization magnitude, a constraint-free phase field model is developed for ferromagnetic materials. Implemented by finite element method, the model is shown to be capable of simulating magnetic domain, vortex, and hysteresis curves. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Dynamic responses of the impact drive mechanism modeled by the distributed parameter system

The piezoelectric actuator has been used for precision positioning from micro-meter down to nano-meter scale. In this paper, the impact drive mechanism (IDM) is designed to achieve a high accuracy and ability in precision positioning motion, where the frictional force is described by the Leuven model combined with the Bounc–Wen model of the hysteresis. The frictional model allows accurate dynamic modeling both in the sliding and the presliding regimes without using switching functions. The governing equations with the hysteresis effects of the distributed parameter system are formulated to obtain the dynamic responses. By using the finite element formulation, numerical solutions due to effects of the piezoelectric element (PE) are provided to compare between the distributed and lumped parameter systems of the IDM. It is shown that the neglect in the mass of the PE will cause the precision errors in the scale of tens nano-meters.     

Normal and tangential continuous elements for least‐squares mixed finite element methods

We consider a finite element discretization of the primal first‐order least‐squares mixed formulation of the second‐order elliptic problem. The unknown variables are displacement and flux, which are approximated by equal‐order elements of the usual continuous element and the normal continuous element, respectively. We show that the error bounds for all variables are optimal. In addition, a field‐based least‐squares finite element method is proposed for the 3D‐magnetostatic problem, where both magnetic field and magnetic flux are taken as two independent variables which are approximated by the tangential continuous and the normal continuous elements, respectively. Coerciveness and optimal error bounds are obtained. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2004.     

Hysteresis loops for para–ferromagnetic phase transitions

In this paper, by an extension of the Ginzburg–Landau theory, we propose a mathematical model describing hard magnets within which we are able to explore the para–ferromagnetic transition and by using the Landau–Lifshitz–Gilbert equation, to study the 3D evolution of magnetic field. Finally, the hysteresis loops are obtained and represented by numerical implementations.     

Finite Element Analysis of General Magnetomechanical Coupling Phenomena

In this paper we consider the derivation of a finite element model for fully–coupled magneto–mechanical boundary value problems based on a vector–valued magnetic potential. The formulation allows for the incorporation of nonlinear, anisotropic and hysteretic constitutive response. A piezomagnetic inclusion problem is computed to test the implementation. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

混凝土材料的粘塑性损伤本构模型

基于不可逆热力学,引入运动硬化、等向硬化和损伤内变量,构造了相应的自由能函数和流动势函数,推导出了混凝土材料的粘塑性损伤本构模型．数值模拟的结果表明,该模型能够避开屈服面和破坏准则的基本假设来描述混凝土材料的以下特性:压缩载荷作用下的体积膨胀现象;应变率敏感性;峰值后由损伤和破坏引起的应力软化和刚度退化现象A·D2由于此模型避开了根据各种变形阶段选择与其相应的本构模型的繁琐计算,因此更便于纳入复杂工况下应力分析有限元程序中．