A Constitutive Model for Magnetostrictive Materials - Theory and Finite Element Implementation

A 3D macroscopic constitutive law for hysteresis effects in magnetostrictive materials is presented and a finite element implementation is provided. The novel aspect of the thermodynamically consistent model is an additive decomposition of the magnetic and the strain field in a reversible and an irreversible part. Employing the irreversible magnetic field is advantageous for a finite element implementation, where the displacements and magnetic scalar potential are the nodal degrees of freedom. To consider the correlation between the irreversible magnetic field and the irreversible strains a one-to-one relation is assumed. The irreversible magnetic field determines as internal variable the movement of the center of a switching surface. This controls the motion of the domain walls during the magnetization process. The evolution of the internal variables is derived from the magnetic enthalpy function by the postulate of maximum dissipation, where the switching surface serves as constraint. The evolution equations are integrated using the backward Euler implicit integration scheme. The constitutive model is implemented in a 3D hexahedral element which provides an algorithmic consistent tangent stiffness matrix. A numerical example demonstrates the capability of the proposed model to reproduce the ferromagnetic hysteresis loops of a Terfenol-D sample. (© 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)     

On domain switching in deformable ferroelectrics, seen as continua with microstructure

We obtain a three-dimensional continuum model for deformable ferroelectric bodies in their polar phase characterized by a spontaneous polarization. This is accomplished by assuming the body as comprised of a continuum with vectorial microstructure: in each point of the body therefore a gross and a fine structure are superposed, the gross structure representing a non linear polarizable elastic body and the vectorial fine structure describing the spontaneous polarization.¶Among the distinctive features of ferroelectric materials, the most interesting is represented by the organization of spontaneous polarization into a domain structure, which minimizes electrostatic energy and which can be modified by the application of electric and deformation fields. This process, called "polarization reversal" or "domain switching", is associated with various hysteresis loops, the most typical being those between spontaneous polarization and electric field (dielectric hysteresis), and between strain and electric field ("butterfly" loop).¶From the balance laws of continua with vectorial microstructure and dissipation inequality we arrive at the evolution equation for the spontaneous polarization which allows for both inertial and dissipative terms and describes domain switching. We postulate a simple interaction mechanism between the spontaneous polarization and the pair electric field, deformation and arrive at, in the quasi-static case, to a minimization problem for a functional which reminds the micromagnetic energy of deformable ferromagnetics.¶For linearized kinematics we study, in the one-dimensional case, stable relative minimizers and give a simple justification for dielectric hysteresis and butterfly loops: under the hypothesis that the domain walls are sharp interfaces, the solutions we find explain the banded twin domains morphology which is typical of many ferroelectrics.     

A finite element modeling of piezoelectric materials with a superimposed stress

Based on the mechanism of domain switching, a three dimensional nonlinear finite element model for piezoelectric materials subjected to electromechancial loading is developed in this contribution. The finally considered model problem deals with differently oriented grains whereby uni-axial, quasi-static cyclic loading is applied. It is assumed that a crystal orientation switches if the reduction in free energy of the grain exceeds a critical energy barrier. The nonlinearity in the small electromechanical loading range is addresses via a polynomial probability function for domain switching. Hysteresis behavior is discussed taking the influence of a superimposed compression state into account. It is observed that the hysteresis loop flattens under the axial compression but elongates under the transverse compression. Irrespective of how the compression is applied, the remnant polarization and as well as the coercive electric field decrease. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Configurational forces in a multiscale approach to piezoelectric materials

Due to the growing interest in determining the macroscopic material response of inhomogeneous materials, computational methods are becoming increasingly concerned with the application of homogenization techniques. In this work, a two-scale classical homogenization of an electro-mechanically coupled material using a FE²-approach is discussed. We explicitly formulated the homogenized coefficients of the elastic, piezoelectric and dielectric tensors for small strain as well as the homogenized remanent strain and remanent polarization. In the homogenization different representative volume elements (RVEs), which capture the micro-structure of the inhomogeneous material, are used to represent the macroscopic material response. Two different schemes are considered. In the first case, domain wall movement is not allowed, but in the second case the movement of the domain walls is taken into account using thermodynamic considerations. Later this technique is used to determine the macroscopic and microscopic configurational forces on defects [2]. These defect situations include the driving force on a crack tip. The effect of the applied electric field on configurational forces at the crack tip is investigated. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On The Formulatio and Numerical Implementation of Dissipative Electro-Mechanics at Large Strains

In the last year an increasing interest in functional materials such as ferroelectric polymers and ceramics has been shown. For those materials viscous effects or electric polarizations cause hysteresis phenomena accompanied with possibly large remanent strains. This paper outlines aspects of the formulation and numerical implementation of dissipative electro-mechanics at large strains. In the first part, we focus on the geometric nature of dissipative electro-mechanics. In a second part, we discuss constitutive assumptions which account for specific problems arising in the geometric nonlinear setting. This concerns the definition of objective energy storage and dissipation functions with suitable symmetries and (weak) convexity properties. With regard to the choice of the internal variables entering these functions, a critical point are kinematic assumptions. Here, we investigate the multiplicative decomposition of the local deformation gradient into reversible and remanent parts as well as the introduction of a remanent metric. In a third part, we summarize details of the constitutive updates as well as the finite element formulations, both on the basis of compact incremental variational principles. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A Finite Element Formulation for the Nonlinear Analysis of Piezoelectric Three-Dimensional Beam Structures

A finite element formulation for a three-dimensional piezoelectric beam which includes geometrical and material nonlinearities is presented. To account for the piezoelectric effect, the coupling between the mechanical stress and the electrical displacement is considered. Based on the Timoshenko theory, an eccentric beam formulation is introduced which provides an efficient model to analyze piezoelectric structures. The geometrically nonlinear assumption allows the calculation of large deformations including buckling analysis. A quadratic approximation of the electric potential through the cross section of the beam ensures the fulfilment of the charge conservation law exactly. This assumption leads to a finite element formulation with six mechanical and five electrical degrees of freedom per node. To take into account the typical ferroelectric hysteresis phenomena, a nonlinear material model is essential. For this purpose, the phenomenological Preisach model is implemented into the beam formulation which provides an efficient determination of the remanent part of the polarization. The applicability of the introduced beam formulation is discussed with respect to available data from literature. (© 2005 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)     

畴极化转动对多晶铁电材料断裂特性的影响

主要基于细观力学方法揭示了畴极化转动对多晶铁电陶瓷的各向异性断裂特性的平均影响。首先,用Eshelby-Mori-Tanaka理论和统计模型分析了无穷大铁电材料体中一椭球夹杂的内、外电弹性场,得到畴极化转动对电弹性场的平均影响;其次,推导了等效多晶铁电陶瓷中含一钱币状裂纹的裂纹扩展力(能量释放率)G_ext,并用它估计了畴极化转动对多晶铁电陶瓷断裂特性的影响。对BaTiO₃陶瓷中裂纹扩展力的计算结果表明,对多晶铁电材料断裂特性分析必须考虑畴极化转动的影响。计算结果得出了与实验相一致的结论:在受较小的力时,外加电场对裂纹扩展产生较大的影响,而且在某种程度上能促进了裂纹扩展。     

Advanced finite element formulations for modeling thin piezoelectric structures

This contribution is concerned with mixed finite element formulations for modeling piezoelectric beam and shell structures. Due to the electromechanical coupling, specific deformation modes are joined with electric field components. In bending dominated problems incompatible approximation functions of these fields cause incorrect results. These effects occur in standard finite element formulations, where interpolation functions of lowest order are used. A mixed variational approach is introduced to overcome these problems. The mixed formulation allows for a consistent approximation of the electromechanical coupled problem. It utilizes six independent fields and could be derived from a Hu-Washizu variational principle. Displacements, rotations and the electric potential are employed as nodal degrees of freedom. According to the Timoshenko theory (beam) and the Reissner-Mindlin theory (shell), the formulations account for constant transversal shear strains. To incorporate three dimensional constitutive relations all transversal components of the electric field and the strain field are enriched by mixed finite element interpolations. Thus the complete piezoelectric coupling is appropriately captured. The common assumption of vanishing transversal stress and dielectric displacement components is enforced in an integral sense. Some numerical examples will demonstrate the capability of the presented finite element formulation. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

An energetic material model for time‐dependent ferroelectric behaviour: existence and uniqueness

We discuss rate‐independent engineering models for the multi‐dimensional behaviour of ferroelectric materials. These models capture the non‐linear and hysteretic behaviour of such materials. We show that these models can be formulated in an energetic framework which is based on the elastic and the electric displacements as reversible variables and on interior, irreversible variables like the remanent polarization. We provide quite general conditions on the constitutive laws which guarantee the existence of a solution. Under more restrictive assumptions we are also able to establish uniqueness results. Copyright © 2006 John Wiley & Sons, Ltd.     

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)     

Relaxation model of piezoelectric materials

Polarization switching inside grains is time dependent. When external applied loading is not quasi-static, macroscopic properties of piezoelectric materials changes with the rate of loading. In this paper, a 2-D micromechanical model is proposed in order to simulate the rate dependent properties of certain perovskite type tetragonal piezoelectric materials based on linear constitutive, nonlinear domain switching, intergranular effects and kinetics models. The material is electrically loaded with an alternating voltage of various frequencies. For the onset of domain switching, energy equation is implemented. Propagation of the domain wall during domain switching in grains is modeled by means of exponential kinetics relation after domain nucleation. Mechanical strain butterfly loops under different frequencies (0.01Hz–1Hz) are simulated. The model gives important insights into the rate dependency of the piezoelectric materials that have been observed in some experiments reported in the literature. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Ⅰ阶梯度损伤理论

从热力学基本定律出发,将应变张量、标量损伤变量、损伤梯度作为Helmholtz自由能函数的状态变量,利用本构泛函展开法在自然状态附近作自由能函数的Taylor展开,未引入附加假设,推导出Ⅰ阶梯度损伤本构方程的一般形式．该形式在损伤为0时可退化为线弹性应力-应变本构方程,在损伤梯度为0时可退化为基于应变等效假设给出的线弹性局部损伤本构方程．一维解析解表明,随着应力增大,损伤场逐步由空间非周期解变为关于空间的类周期解,类周期解的峰值区域形成局部化带．局部化带内的损伤变量将不同于局部化带外的损伤变量,由此可以反映出介质的局部化特征．损伤局部化并不是与损伤同时发生,而是在损伤发生后逐渐显现出来,模型的局部化机制开始启动;损伤局部化的宽度同内部特征长度成正比．     

Simulation of Atomic Force Microscopy for investigating BaTiO3 and LiMn2O4 nanostructures based on Phase Field Approach

Atomic Force Microscopy (AFM) probes the surface features of specimens using an extremely sharp tip scanning the sample surface while the force is applied. AFM is also widely used for investigating the electrically non-conductive materials by applying an electric potential on the tip. Piezoresponse Force Microscopy (PFM) and Electrochemical Strain Microscopy (ESM) are variants of AFM for different materials. Both PFM and ESM signals are obtained by observing the displacement of the tip when applying electric fields during the scanning process. The PFM technique is based on converse piezoelectric effect of ferroelectrics and the ESM technique is based on electrochemical coupling in solid ionic conductors. In this work, two continuum-mechanical formulations for simulation of PFM and ESM are discussed. In the first model, for PFM simulation, a phase field approach based on the Allen-Cahn equation for non-conserved order parameters is employed for ferroelectrics. Here, the polarization vector is chosen as order parameter. Since ferroelectrics have highly anisotropic properties, this model accounts for transversely isotropic symmetry using an invariant formulation. The polarization switching behavior under the electric field will be discussed with some numerical examples. In the simulation of ESM, we employ a constitutive model based on the work of Bohn et al. [8] for the modeling of lithium manganese dioxide LiMn₂O₄ (LMO). It simulates the deformation of the LMO particle according to an applied voltage and the evolution of lithium concentration after removing a DC pulse. The modeling results are compared to experimental data. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Experimental and Numerical Investigation of Nonlinear Piezoelectric Effects Including Minor Hystereses

Piezoelectric ceramics are often used in active structures for shape and vibration control. Since the operation range is not limited to small signals the nonlinear behaviour of the actuator under high electric loads has to be known. There are several approaches in literature to model the hysteretic effects, each having its assets and drawbacks. When a model is able to reproduce the minor loops of the strain - electric field hysteresis, it often lacks the consideration of stress dependence which is fundamental for actuators attached to elastic structures. On the other hand constitutive models which take into account all ferroelectric and ferroelastic effects are not capable of representing the minor hystereses in acceptable calculation times. In this work a phenomenological constitutive model is verified using the experimental data of an active plate structure. Therefore, the ceramic is characterised under mechanically unconstrained conditions and afterwards attached onto a steel plate. The bonding to the substructure leads to a mechanical stress depending on the actuation state of the ceramic. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Characterisation and Modelling of Piezoceramic Actuator Hystereses

Piezoelectric ceramics are often used as actuators in smart structures technology. In the vast majority of papers dealing with this topic only linear constitutive relations are used. However, the electric field-strain relations of such actuators show hysteretic behaviour, which means that the piezoelectric coupling coefficient is not constant. In this study the hysteresis of a mechanically unconstrained actuator is obtained using the Michelson interferometry. The hysteretic behaviour is modelled by a Preisach model. Using these experimental data, for the modelling of an active structure with embedded piezoelectric actuators the actual coupling coefficient can then be determined with the help of the Preisach model. With this procedure the actuation strain of an embedded actuator, including the physical nonlinearities, can be calculated using the material characteristics obtained for an unconstrained actuator. For an experimental validation of the method outlined above, a Lead Zirconate Titanate (PZT) actuator is characterised experimentally and then glued to a cantilever beam. Then, the tip displacement of the actuated beam is determined experimentally and simulated numerically using the above method. The experimental and numerical results agree reasonably well if the shear lag due to the bonding layer between the actuator and the structure is taken into consideration. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A mixed finite element formulation for piezoelectric shells

A finite element formulation to analyze piezoelectric shell problems is presented. A reference surface of the shell is modelled with a four node element. Each node possesses six mechanical degrees of freedom, three displacements and three rotations, and one electric degree of freedom, which is the difference of the electric potential in thickness direction. The formulation is based on the mixed field variational principle of Hu-Washizu. The independent fields are displacements u , electric potential φ, strains E , electric field E , stresses S and dielectric displacements D . The mixed formulation allows an interpolation of the strains and the electric field in thickness direction. Accordingly a three-dimensional material law is incorporated in the variational formulation. It is remarked that no simplification regarding the constitutive law is assumed. The formulation allows the consideration of arbitrary constitutive relations. The normal zero stress condition and the normal zero dielectric displacement condition are enforced by the independent stress and dielectric displacement fields. They are defined as zero in thickness direction. The present shell element fulfills the important patch tests: the in-plane, bending and shear test. Some numerical examples demonstrate the applicability of the present piezoelectric shell element. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Formulation and Numerical Implementation of Ferroelectrics at Large Strains

In recent years, an increasing interest has been shown in functional materials such as ferroelectric polymers. For such materials, viscous effects and electric polarizations cause hysteresis phenomena accompanied with possibly large remanent strains and rotations. Ferroelectric polymers have many attractive characteristics. They are light, inexpensive, fracture tolerant, and pliable. Furthermore, they can be manufactured into almost any conceivable shape and their properties can be tailored to suit a broad range of requirements. In this work, continuous and discrete variational formulations are exploited for the treatment of the non-linear dissipative response of ferroelectric polymers under electrical loading. The point of departure is a general internal variable formulation that determines the hysteretic response of this class of materials in terms of an energy storage and a rate-dependent dissipation function. The ferroelectric constitutive assumptions, which account for specific problems arising in the geometric nonlinear setting, are discussed. With regard to the choice of the internal variables, a critical factor is the kinematic assumptions. Here, we propose the multiplicative decomposition of the local deformation gradient into reversible and remanent parts, where the latter is characterized by a metric tensor. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)