A phenomenological constitutive model for ferroelastic and ferroelectric hysteresis effects in ferroelectric ceramics

This paper is concerned with a macroscopic constitutive law for domain switching effects, which occur in ferroelectric ceramics. The three-dimensional model is thermodynamically consistent and is determined by two scalar valued functions: the Helmholtz free energy and a switching surface. In a kinematic hardening process the movement of the center of the switching surface is controlled by internal variables. In common usage, the remanent polarization and the irreversible strain are employed as internal variables. The novel aspect of the present work is to introduce an irreversible electric field, which serves instead of the remanent polarization as internal variable. The irreversible electric field has only theoretical meaning, but it makes the formulation very suitable for a finite element implementation, where displacements and the electric potential are the nodal degrees of freedom. The paper presents an appropriate implementation into a hexahedral finite brick element. The uni-axial constitutive model successfully reproduces the ferroelastic and the ferroelectric hysteresis as well as the butterfly hysteresis for ferroelectric ceramics. Furthermore it accounts for the mechanical depolarization effect, which occurs if the polarized ferroelectric ceramic is subjected to a compression stress.     

A unified framework for modeling hysteresis in ferroic materials

This paper addresses the development of a unified framework for quantifying hysteresis and constitutive nonlinearities inherent to ferroelectric, ferromagnetic and ferroelastic materials. Because the mechanisms which produce hysteresis vary substantially at the microscopic level, it is more natural to initiate model development at the mesoscopic, or lattice, level where the materials share common energy properties along with analogous domain structures. In the first step of the model development, Helmholtz and Gibbs energy relations are combined with Boltzmann theory to construct mesoscopic models which quantify the local average polarization, magnetization and strains in ferroelectric, ferromagnetic and ferroelastic materials. In the second step of the development, stochastic homogenization techniques are invoked to construct unified macroscopic models for nonhomogeneous, polycrystalline compounds exhibiting nonuniform effective fields. The combination of energy analysis and homogenization techniques produces low-order models in which a number of parameters can be correlated with physical attributes of measured data. Furthermore, the development of a unified modeling framework applicable to a broad range of ferroic compounds facilitates material characterization, transducer development, and model-based control design. Attributes of the models are illustrated through comparison with piezoceramic, magnetostrictive and shape memory alloy data and prediction of material behavior.     

A SIMPLE CONSTITUTIVE MODEL FOR FERROELECTRIC CERAMICS UNDER ELECTRICAL/MECHANICAL LOADING

A simple phenomenological model is developed for describing the macroscopic constitutive response of ferroelectric materials based on consideration of the fact that domain switching is a progressive evolution process with loading. The volume fraction of domain switching is taken as an internal variable, which is derived from the domain nucleation theory. The proposed theory can simulate the dielectric hysteresis, reversed butterfly hysteresis, nonlinear strain-stress hysteresis, as well as electric displacement-stress relation of ferroelectric materials. Its comparison with experimental results and two other theoretical models reveals that the model presented can well predict the nonlinear hysteresis of ferroelectrics under electrical or mechanical loading.     

Non-linear constitutive relations for magnetostrictive materials

In this paper, non-linear deformation behavior of magnetostrictive materials is studied and three magnetoelastic coupling constitutive models are developed. The standard square (SS) constitutive model is developed by means of truncating the polynomial expansion of the Gibbs free energy. The hyperbolic tangent (HT) constitutive equations, which involve a hyperbolic tangent magnetic-field dependence, are proposed to model the magnetic-field-induced strain saturation of magnetostrictive materials in the region of intense magnetic fields. A new model based on density of domain switching (DDS) is established in terms of the basic truth that magnetic domain switching underlies magnetostrictive deformation. In this model, it is assumed that the relation between density of domain switching, defined by the quantity of magnetic domains switched by per unit magnetic field and magnetic field can be described by a density function with normal distribution. The moduli in these constitutive models can be determined by a material function that is proposed to describe the dependence of the peak piezomagnetic coefficient on the compressive pre-stress for one-dimensional cases based on the experimental results published. The accuracy of the non-linear constitutive relations is evaluated by comparing the theoretical values with experimental results of a Terfenol-D rod operated under both compressive pre-stress and bias magnetic field. Results indicate that the SS constitutive equations can accurately predict the experimental results under a low or moderate magnetic field while the HT model can, to some extent, reflect the trend of saturation of magnetostrictive strain under a high magnetic field. The model based on DDS, which is more effective in simulating the experimental curves, can capture the main characteristics of the mechanism of magnetoelastic coupling deformation of a Terfenol-D rod, such as the notable dependence of magnetoelastic response on external stress and the saturation of magnetostrictive strain under intense magnetic fields. In addition, the SS constitutive relation for a general three-dimensional problem is discussed and an approach to characterize the modulus tensors is proposed.     

Three-dimensional finite element simulations of ferroelectric polycrystals under electrical and mechanical loading

Complex, non-linear, irreversible, hysteretic behavior of polycrystalline ferroelectric materials under a combined electro-mechanical loading is a result of domain wall motion, causing simultaneous expansion and contraction of unlike domains, grain sub-divisions that have distinct spontaneous polarization and strain. In this paper, a 3-dimensional finite element method is used to simulate such a polycrystalline ferroelectric under electrical and mechanical loading. A constitutive law due to Huber et al. [1999. A constitutive model for ferroelectric polycrystals. J. Mech. Phys. Solids 47, 1663-1697] for switching by domain wall motion in multidomain ferroelectric single crystals is employed in our model to represent each grain, and the finite element method is used to solve the governing conditions of mechanical equilibrium and Gauss's law. The results provide the average behavior for the polycrystalline ceramic. We compare the outcomes predicted by this model with the available experimental data for various electromechanical loading conditions. The qualitative features of ferroelectric switching are predicted well, including hysteresis and butterfly loops, the effect on them of mechanical compression, and the response of the polycrystal to non-proportional electrical loading.     

Phenomenological model for the macroscopical material behavior of ferroelectric ceramics

A thermodynamically consistent phenomenological model for the simulation of the macroscopic behavior of ferroelectric polycrystalline ceramics is presented. It is based on the choice of microscopically motivated internal state variables, which describe the texture and the polarization state of the polycrystal. Saturation states are defined for the internal state variables. The linear material behavior is modelled by a transversely isotropic piezoelectric constitutive law, where the anisotropy is history dependent. For non-linear irreversible processes, a switching function and associated evolution rules are applied, satisfying the principle of maximum ferroelectric dissipation. Saturation is modelled by the use of energy-barrier functions in the electric enthalpy density function. Numerical examples demonstrate the capability of the proposed model, to predict the typical experimental results.     

磁致伸缩材料的非线性本构关系

给出了磁致伸缩材料的两个非线性本构关系,即标准平方型和双曲正切型。在确定一维问题的本构系数时,基于已有的实验结果,引进一个材料函数,用来描述磁致伸缩材料的压磁系数随预应力变化的关系。将非 线性本构关系的理论模型计算结果与实验曲线对比,结果表明标准平方型本构关系在中低磁场下能精确地模拟实验曲线,而双曲正切型本构关系在高磁场时能反映材料的磁致应变饱和现象。讨论了在标准平方型本构的一般三维情形,给出了确定本构系数的方法。     

Damage model for flexural strength variation of ferroelectric materials induced by electric field

Ferroelectric material is being used widely in aerospace related applications, namely active control of large flexible space trusses, structure acoustics and helicopter rotary blades. Degradation of material properties affects the safety and reliability of these structural members. In particular, the flexural strength of ferroelectric materials tend to reduce when they are exposed to electric field. Experimentation cannot always explain the underlying reasons of these physical phenomena. To that end, a damage model is presented to study the switching process induced by applied electric field that might be the cause of damage. Furthermore, evolution of internal damage is considered using irreversible thermodynamics. Two relations are assumed to associate the flexural strength with the damage variable. Available test data for three-point-bend test are used to test the model for a PZT material exposed to electric field. The results reveal the validity of the present damage model while the critical strain energy density criterion led to a more reasonable correlation.     

Time dependent behavior of ferroelectric materials undergoing changes in their material properties with electric field and temperature

This paper presents a time dependent polarization constitutive model suitable for predicting nonlinear polarization and electro-mechanical strain responses of ferroelectric materials subject to various histories of electric fields. The constitutive model is derived based on a single integral form with nonlinear (electric field and temperature dependent) integrand. The total polarization consists of the time-dependent and residual components. The residual component of the polarization is due to polarization switching in the ferroelectric materials. We use an ‘internal clock’ concept to incorporate the effect of electric field on the rate of polarization. The corresponding strain response is determined through the use of third order piezoelectric constant and/or fourth order electrostrictive constant that vary with polarization stage. It is assumed that in absence of polarization, both piezoelectric and electrostrictive constants are zero. To incorporate the effect of temperature on the overall polarization behavior all material parameters in the constitutive model are allowed to change with the ambient temperature. We present numerical studies on the effect of time, temperature, and electric field on the response of ferroelectric material followed by verification of the constitutive model. Experimental data on lead zirconate titanate (PZT) materials available in the literature are used to verify the model.     

A nonlinear constitutive model for magnetostrictive materials

A general nonlinear constitutive model is proposed for magnetostrictive materials, based on the important physical fact that a nonlinear part of the elastic strain produced by a pre-stress is related to the magnetic domain rotation or movement and is responsible for the change of the maximum magnetostrictive strain with the pre-stress. To avoid the complicity of determining the tensor function describing the nonlinear elastic strain part, this paper proposes a simplified model by means of linearizing the nonlinear function. For the convenience of engineering applications, the expressions of the 3-D (bulk), 2-D (film) and 1-D (rod) models are, respectively, given for an isotropic material and their applicable ranges are also discussed. By comparison with the experimental data of a Terfenol-D rod, it is found that the proposed model can accurately predict the magnetostrictive strain curves in low, moderate and high magnetic field regions for various compressive pre-stress levels. The numerical simulation further illustrates that, for either magnetostrictive rods or thin films, the proposed model can effectively describe the effects of the pre-stress or residual stress on the magnetization and magnetostrictive strain curves, while none of the known models can capture all of them. Therefore, the proposed model enjoys higher precision and wider applicability than the previous models, especially in the region of the high field.The project supported by the National Natural Science Foundation of China (10132010 and 90405005)     

铁电陶瓷宏观单轴力电行为的双面模型

铁电陶瓷以其优越的力电耦合性作为新型的智能材料使用. 提出基于弹塑性双面理论的宏观铁电本构模型. 根据铁电陶瓷内部电畴在外电场和机械场作用下的微观运动，在宏观上除引入材料的畴变面外，还首次引入饱和面，并考虑以畴变面与饱和面之间的广义距离来表征铁电陶瓷的非线性行为. 数值计算结果与实验数据的比较表明所提出的初步理论可适当地反映力电加载下铁电陶瓷的宏观非线性行为.     

A ferroelectric and ferroelastic 3D hexahedral curvilinear finite element

An isoparametric 3D electromechanical hexahedral finite element integrating a 3D phenomenological ferroelectric and ferroelastic constitutive law for domain switching effects is proposed. The model presents two internal variables which are the ferroelectric polarization (related to the electric field) and the ferroelastic strain (related to the mechanical stress). An implicit integration technique of the constitutive equations based on the return-mapping algorithm is developed. The mechanical strain tensor and the electric field vector are expressed in a curvilinear coordinate system in order to handle the transverse isotropy behavior of ferroelectric ceramics. The hexahedral finite element is implemented into the commercial finite element code Abaqus® via the subroutine user element. Some linear (piezoelectric) and non linear (ferroelectric and ferroelastic) benchmarks are considered as validation tests.     

基于深度强化学习算法的颗粒材料应力?应变关系数据驱动模拟研究

颗粒材料的宏观力学行为受颗粒组分等材料参数, 孔隙率、配位数等状态参数的影响, 同时又具备复杂的加载路径和加载历史相关性, 建立包含多个内变量以及各变量间相互关联的颗粒材料本构模型是一个重要的科学难题. 不同于传统的基于屈服面、流动法则和硬化函数框架下的唯象本构模型, 本文基于颗粒物质力学的研究基础, 以颗粒材料平均孔隙率、细观组构参数和弹性刚度参数作为内变量, 结合深度学习方法建立以有向图表征的数据本构模型. 有向图中以不同的链接网络表示不同的内变量信息流动方向, 各个内变量间的映射关系采用循环神经网络来建立, 将各个神经网络相互组合, 形成包含不同内变量且具有不同预测能力的本构模型. 该本构模型的建立过程等价于在众多可能的内变量链接关系空间中寻找最能描述实际材料宏观应力应变行为的优化问题. 因此, 可将有向图本构模型的建立过程看作“马尔可夫决策过程”, 采用深度强化学习算法构建有向图的内变量链接组合优化过程, 具体采用AlphaGo Zero算法自动寻找最优的颗粒材料数据驱动本构模型建模路径. 研究结果表明, 采用有向图和深度强化学习算法可建立起完全依靠“数据驱动”的颗粒材料应力?应变关系. 此外, 本方法提供了一种将不同理论模型从数据角度统一起来, 且基于人工智能算法发展更优模型的研究思路, 可为相似问题的研究提供借鉴.      

Experimental and theoretical study of the nonlinear response of agiant magnetostrictive rod

In this paper, a magnetomechanical coupling constitutive relation of the giant magnetostrictive material was investigated experimentally and theoretically. A grain-oriented magnetostrictive rod of iron and rare earth was tested under a combined magnetomechanical loading. Two types of experimental curves were obtained, i.e., the magnetostrictive curve of the extensional strain vs the magnetic field, and the curve of the magnetic polarization intensity vs the pre-stress. A new theoretical constitutive model, based on the density of domain switching, is developed. Comparison of the theoretical predictions with the experimental results indicates that this model can capture the main characteristics of the magnetoelastic coupling deformation of a giant magnetostrictive rod. The project supported by the National Natural Science Foundation of China (10025209, 10132010, 10102007)     

基于唯象相变理论的磁致伸缩材料耦合滞回动力学模拟

本文基于改进的Landau唯象相变理论,构造一个耦合的非线性常微分方程模型来模拟一维磁致伸缩材料的磁滞动态特性。模型的构造通过引入一个非凸的自由能函数来模拟磁致伸缩磁材料中不可逆的磁极化翻转与磁致应变,该自由能函数的每一个局部极小值都对应材料的一个磁化方向。通过热力学平衡条件建立能刻画磁致伸缩效应的非线性本构关系。所构造的模型成功地模拟出了磁场与弹性场之间的磁滞曲线和蝶形曲线,并采用实验结果对模型进行了验证。     

Continuous symmetries and constitutive laws of dissipative materials within a thermodynamic framework of relaxation: Part I: Formal aspects

An analysis of the continuous symmetries of the constitutive laws of inelastic materials written within a thermodynamical framework of relaxation is performed. This framework relies on the generalization of Gibb’s relationship outside the equilibrium of a uniform system, and the use of the fluctuation theory to model the material dissipation due to its internal microstructure change [Cunat, C., 2001. The DNLR approach and relaxation phenomena. Part I – Historical account and DNLR formalism. Mech. Time-depend. Mater. 5, 39–65]. The approach leads to a viscoelastic like formulation for small deformations, and changes gradually for finite strains towards elastoviscoplasticity (with or without damage) via a dependence of characteristic times with the loading path, in a way similar to the endochronic approach developed by Valanis [Valanis, K.C., 1975. On the fundations of the endochronic theory of viscoplasticity. Arch. Mech. 27, 857–868]. The present thermodynamic framework has been previously applied to elastoviscoplastic materials under cyclic and non-proportional loadings [Dieng, L., Abdul-Latif, A., Haboussi, M., Cunat, C., 2005b. Cyclic plasticity modeling with the distribution of non-linear relaxations approach. Int. J. Plasticity 21, 353–379]. The constitutive laws split into the state laws relating intensive variables (thermodynamics forces) to extensive-like variables, and the complementary evolution laws of the internal variables associated to the dissipative mechanisms. An interpretation of a non-equilibrium thermodynamic approach of irreversible processes in terms of an extremum principle is proposed, associated to a Lagrangian functional. It is shown that one possible choice for the Lagrangian kernel is the material derivative of the internal energy density, augmented by a complementary term that accounts for the evolution laws of the internal variables. Interpreting the material behavior during the non-equilibrium evolution as the Euler–Lagrange equations of the resulting action integral, a differential condition expressing both the local and variational symmetries encapsulated into the Lagrangian formulation is formulated. It is further shown that both symmetry conditions are fully equivalent along the optimal path corresponding to the satisfaction of the constitutive laws. In terms of both practical and methodological aspects, the predictive nature of the symmetry analysis is highlighted, as a systematic tool for the exploitation of the constitutive response. Its performance and utility are exemplified by the construction of a time–temperature equivalence principle for a dry viscous polymer (PA66); the calculated shift factor is shown to well agree with the empirical shift factor given by Williams–Landel–Ferry (WLF) expression. A systematic interpretation of the calculated symmetry groups of the constitutive laws in terms of master curves for various plastic and viscoplastic materials shall be presented in a forthcoming contribution.     

Load-induced oriented damage and anisotropy of rock-like materials

Theoretical model for deformability of brittle rock-like materials in the presence of an oriented damage of their internal structure is formulated and verified experimentally. This model is based on the assumption that non-linearity of the stress–strain curves of these materials is a result of irreversible process of oriented damage growth. It was also assumed that a material response, represented by the strain tensor, is a function of two tensorial variables: the stress tensor and the damage effect tensor that is responsible for the current state of the internal structure of the material. The explicit form of the respective non-linear stress–strain relations that account for the appropriate damage evolution equation was obtained by employing the theory of tensor function representations and by using the results of own experiments on damage growth. Such an oriented damage that grows in the material, described by the second order symmetric damage effect tensor, results in gradual development of the material anisotropy. The validity of the constitutive equations proposed was verified by using the available experimental results for concrete subjected to the plane state of stress. The relevant experimental data for sandstone and concrete subjected to tri-axial state of stress were also used.     

A 3D mechanistic model for brittle materials containing evolving flaw distributions under dynamic multiaxial loading

We present a validated fully 3D mechanism-based micromechanical constitutive model for brittle solids under dynamic multiaxial loading conditions. Flaw statistics are explicitly incorporated through a defect density, and evolving flaw distributions in both orientation and size. Interactions among cracks are modeled by means of a crack-matrix-effective-medium approach. A tensorial damage parameter is defined based upon the crack length and orientation development under local effective stress fields. At low confining stresses, the wing-cracking mechanism dominates, leading to the degradation of the modulus and peak strength of the material, whereas at high enough confining stresses, the cracking mechanism is completely shut-down and dislocation mechanisms become dominant. The model handles general multiaxial stress states, accounts for evolving internal variables in the form of evolving flaw size and orientation distributions, includes evolving anisotropic damage and irreversible damage strains in a thermodynamically consistent fashion, incorporates rate-dependence through the micromechanics, and includes dynamic bulking based on independent experimental data. Simulation results are discussed and compared with experimental results on one specific structural ceramic, aluminum nitride. We demonstrate that this 3D constitutive model is capable of capturing the general constitutive response of structural ceramics.