Micromechanics of ferroelectric polymer-based electrostrictive composites

Electrostriction refers to the strain induced in a dielectric by electric polarization, which is usually very small for practical application. In this paper, we present a micromechanical analysis on the effective electrostriction of a ferroelectric polyvinylidene fluoride trifluoroethylene [P(VDF-TrFE)] polymer-based composite, where the exact connections between the effective electrostrictive coefficients and effective elastic moduli are established, and numerical algorithm for the prediction of the effective electrostrictive coefficients of the composite in terms of its microstructural information is developed. From our calculations, enhanced electrostriction in the composite has been demonstrated, and optimal microstructure for electrostriction enhancement has been identified. Our analysis provides a mechanism for the electrostriction enhancement, where the electrostrictive strain several times higher than that of polymer matrix can be obtained, if the microstructure of the composites can be carefully tailored.     

A study of electromechanical switching in ferroelectric single crystals

This article documents both modeling and experimental studies developed to investigate the switching behavior of ferroelectric single crystals. The theoretical model makes a priori ansatz that switching follows the evolution of a particular domain pattern. The choice of this configuration is dictated by the requirement that domains remain compatible during evolution, giving rise to a low-energy path for the overall switching. The construction of this pattern is achieved using multirank laminates. It offers an advantage of specifying different types of domain wall movements, leading to a distinction for the switching types. A loading experiment is performed on a barium titanate (BaTiO₃) single crystal with a constant compressive stress and a cyclic electric field. Both 180 and 90 coercive fields are measured as input parameters required for the theoretical framework. The simulation results show good agreement with the observed strains measured by the present and other available experiments. It is found that depolarization has a non-trivial influence on attainable actuation strains.     

On the coupling of non-linear normal modes

The phenomenon of internal resonance is known as the exchange of energy between the modes and the existence of coupled-mode response under a single-mode excitation. This phenomenon is observed whenever a non-linear normal mode loses its stability, called the modal coupling. The details of modal coupling are formulated in the free vibrations of two-degree-of-freedom systems, and compared with internal resonance. The theory is based on the structural change in Poincaré map due to the stability change of normal modes. It is shown that every change in stability of normal modes gives rise to a pitchfork or a period-doubling bifurcation. The functional form is derived to compute the coupled modes by the method of harmonic balance. Examples are given to describe the procedure of stability analysis of non-linear normal modes, to compute the coupled modes, and then to demonstrate that results of internal resonances can be derived by model coupling. Other examples are given to demonstrate that the results of some modal couplings cannot be obtained by internal resonances.     

Numerical simulation of domain switching in multilayer ferroelectric actuators

Micromechanical finite element methods are developed based on a nonlinear constitutive model of ferroelectric polycrystals. Electromechanical behaviors ahead of an internal electrode tip are numerically simulated in multilayer ferroelectric actuators. Around the electrode edge, the nonuniform electric field generates a concentration of stress due to the incompatible strain as well as spontaneous strain. The preferred domain switching enhances the concentration of residual stress and may cause the actuators to crack. An electrically permeable crack emanating from an internal electrode is analyzed. A large scale domain switching zone is found in the vicinity of crack tips. The larger the actuating strain and electric field are, the larger the switching zone will be. The size of switching zone even reaches the scale of crack length with increasing electromechanical loading.     

Development of inter- and intragranular stresses during switching of ferroelectric polycrystals

Ferroelectrics are crystalline inorganic materials consisting of domains with different directions of spontaneous polarization. By application of sufficiently high electric fields, these domains can switch into a common direction, thus making the material piezoelectric. Due to ferroelasticity, the domains can be also switched into different states by the application of mechanical stress. In polycrystalline materials, as used in most applications, electric and stress fields interact so as to maintain compatibility. We study the influence of grain-to-grain interactions on the overall and local switching behavior and in particular the induced stresses inside grains and across grain boundaries. The behavior inside each grain is represented by the single-crystal model of [Huber, J.E., Fleck, N.A., Landis, C.M., McMeeking, R.M., 1999. A constitutive model for ferroelectric polycrystals. Journal of the Mechanics and Physics of Solids 47 (8), 1663–1697] and the polycrystal response is obtained through a two-dimensional multi-grain model in which grains are represented individually. We investigate the effect of random grain orientations, both in the plane of consideration and in three directions, and compare plane strain with plane stress conditions. It is found that the overall piezoelectric response under electric loading is not dependent only on the intra- and intergranular stresses in the plane but is also significantly affected by stresses in through-thickness direction.     

A switching fractional calculus-based controller for normal non-linear dynamical systems

In this paper, a fractional calculus-based terminal sliding mode controller is introduced for finite-time control of non-autonomous non-linear dynamical systems in the canonical form. A fractional terminal switching manifold which is appropriate for canonical integer-order systems is firstly designed. Then some conditions are provided to avoid the inherent singularities of the conventional terminal sliding manifolds. A non-smooth Lyapunov function is adopted to prove the finite time stability and convergence of the sliding mode dynamics. Afterward, based on the sliding mode control theory, an equivalent control and a discontinuous control law are designed to guarantee the occurrence of the sliding motion in finite time. The proposed control scheme uses only one control input to stabilize the system. The proposed controller is also robust against system uncertainties and external disturbances. Two illustrative examples show the effectiveness and applicability of the proposed fractional finite-time control strategy. It is worth noting that the proposed sliding mode controller can be applied for control and stabilization of a large class of non-autonomous non-linear uncertain canonical systems.     

Micromechanics of ferroelectric domain switching behavior: Part II: Constitutive relations and hysteresis

A constitutive relation is developed to describe the nonlinear behavior of ferroelectric ceramics subjected to external stress and electric field. The theoretical development considers each domain as an inclusion. The Helmholtz and Gibbs free energy of the constituent element are derived by using a micromechanics approach. They are functionals of the orientation distribution function (ODF) that represents the domain distribution patterns. By applying the internal variable theory and expanding ODF in Fourier series, the yield condition, evolution of ODF, and constitutive relation are obtained. Theoretical results agree with experiments.     

An optical method to assess electromechanical coupling in ferroelectric ceramics

A new technique is described, which allows the assessment of elastic and inelastic regions around a macroscopic defect in ferroelectric-ferroelastic ceramics. The accuracy and robustness of the method are demonstrated on a PZT plate with a centered hole subjected to uni-axial compressive stresses. From the electrical potential distribution on the sample surface, the mechanical response of the material is obtained at different load levels.     

Macroscopical non-linear material model for ferroelectric materials inside a hybrid finite element formulation

A new approach for modeling hysteretic non-linear ferroelectric ceramics is presented, based on a fully ferroelectric/ferroelastic coupled macroscopic material model. The material behavior is described by a set of yield functions and the history dependence is stored in internal state variables representing the remanent polarization and the remanent strain. For the solution of the electromechanical coupled boundary value problem, a hybrid finite element formulation is used. Inside this formulation the electric displacement is available as nodal quantity (i.e. degree of freedom) which is used instead of the electric field to determine the evolution of remanent polarization. This involves naturally the electromechanical coupling. A highly efficient integration technique of the constitutive equations, defining a system of ordinary differential equations, is obtained by a customized return mapping algorithm. Due to some simplifications of the algorithm, an analytical solution can be calculated. The automatic differentiation technique is used to obtain the consistent tangent operator. Altogether this has been implemented into the finite element code FEAP via a user element. Extensive verification tests are performed in this work to evaluate the behavior of the material model under pure electrical and mechanical as well as coupled and multi-axial loading conditions.     

The influence of mechanical constraint upon the switching of a ferroelectric memory capacitor

The role of mechanical constraint upon the switching response of a ferroelectric thin film memory capacitor is explored. The memory capacitor is represented by a two dimensional ferroelectric island whose non-linear behaviour is modelled by a crystal plasticity constitutive law within the finite element method. The switching response of the device, in terms of remnant charge storage, is determined as a function of geometry and constraint. Various types of constraint on the ferroelectric capacitor are considered, including the presence of a silicon dioxide passivation layer, a silicon substrate and metallic electrodes. The effect of the relative resistance to 90 degree switching and 180 degree switching is also explored in a tetragonal ferroelectric device. Throughout the study, the finite element calculations are compared with the behaviour of a material element subjected to various degrees of mechanical constraint.     

Linear and non-linear viscoelastic properties of ethylene vinyl acetate/nano-crystalline cellulose composites

This paper reports on the melt rheological properties of ethylene vinyl acetate containing between 0 and 10 wt.% of nano-crystalline cellulose (NCC). A complete set of rheological tests including frequency sweeps, shear transients, and uniaxial elongations was performed. Frequency sweeps showed that at low frequencies, a pseudo solid-like behavior was obtained for NCC concentrations higher than 5%. This behavior was related to hydrogen bonding between NCC particles and the creation of particle networks as the result of particle–particle interactions. For transient shear tests, all compositions presented a stress overshoot at high shear rates before reaching a steady state. It was found that the amplitude of this overshoot depends on both NCC content and shear rate. On the other hand, the time to reach the maximum was found to be highly shear rate dependent but concentration dependence was rather weak. For uniaxial extensional flow, higher extensional viscosity was observed with increasing NCC content. On the other hand, strain hardening was found to decrease with increasing NCC content.     

Tuning of non-uniform switch toughening in ferroelectric composites by an electric field

This paper deals with a mode III interfacial crack subject to anti-plane stress and in-plane electric fields. The analysis concentrates on the tuning of fracture toughness from non-uniform ferroelectric-ferroelastic domain switch-ing by an electric field. The electric loading changes the size of the asymmetric switching zone. Employing the weight function method, we obtain the electrically-dependent switch toughening for stationary and quasi-static growing interfacial cracks, respectively. Multi-domain solutions are derived for non-poled and fully-poled ferroelectric composites. Numer-ical results are presented on the electric field tuning of the critical applied stress intensity factor. The research provides ways to optimize fracture properties of ferroelectric compos-ites by altering the electric field.     

The non-linear dynamic response of the Euler-Bernoulli beam with an arbitrary number of switching cracks

In this study the non-linear dynamic response of the Euler-Bernoulli beam in presence of multiple concentrated switching cracks (i.e. cracks that are either fully open or fully closed) is addressed. The overall behaviour of such a beam is non-linear due to the opening and closing of the cracks during the dynamic response; however, it can be regarded as a sequence of linear phases each of them characterised by different number and positions of the cracks in open state. In the paper the non-linear response of the beam with switching cracks is evaluated by determining the exact modal properties of the beam in each linear phase and evaluating the corresponding time history linear response through modal superposition analysis. Appropriate initial conditions at the instant of transition between two successive linear phases have been considered and an energy control has been enforced with the aim of establishing the minimum number of linear modes that must be taken into account in order to obtain accurate results. Some numerical applications are presented in order to illustrate the efficiency of the proposed approach for the evaluation of the non-linear dynamic response of beams with multiple switching cracks. In particular, the behaviour under different boundary conditions both for harmonic loading and free vibrations has been investigated.     

Nonlinear electromechanical fields and localized polarization switching of 1-3 piezoelectric/polymer composites

This paper examines the nonlinear electromechanical response of 1-3 piezoelectric/polymer composites. The piezocomposites contain square or circular piezoelectric rods in an epoxy matrix. Experiments were conducted to measure the displacement versus electric field curves, using the device specimen of the 1-3 piezocomposites. Three dimensional finite element analysis was also carried out to study the electromechanical fields in the 1-3 piezocomposites by introducing a model for polarization switching. Comparison was then made between simulation and experiment.     

Electroelastic field concentrations and polarization switching induced by circular electrode at the interface of piezoelectric disk composites

The electromechanical field concentrations due to circular electrode at the interface of piezoelectric disk composites have been discussed. This paper consists of two parts. In the first part, the problem of an internal electrode embedded at the interface of two dissimilar semi-infinite piezoelectric solids was formulated by means of Hankel transforms, and the solution was solved exactly. In the second part, finite element analysis was carried out to study electromechanical response in piezoelectric disk composites containing a circular electrode at the interface by introducing a model for polarization switching in local areas of field concentrations. A nonlinear behavior induced by localized polarization switching was observed between the strain and the voltage applied to the electrode.     

Effects of interface contacts on the magneto electro-elastic coupling for fiber reinforced composites

Imperfect bonding between constituents is studied where displacements, electric and magnetic static potentials are considered to have a jump proportional to the normal component of the mechanical traction, electric displacement and magnetic flux. This condition may model various interface damages or the thin glue layer between two adjacent phases. They are termed as the mechanically compliant, dielectrically weakly capacitance and magnetically weakly inductance at the interface. It is shown that while the more imperfect the interface is, the overall properties become weaker, such as longitudinal shear stiffness, out-of-plane piezoelectric coupling, and magnetoelectric coupling. Out-of-plane piezomagnetic coupling, transverse dielectric permittivity and transverse dielectric permeability exhibit no influence by imperfect bonding. The imperfect interface proposed is mimicked by the springs, capacitors and inductances for the mechanical, electric and magnetic interaction between the phases and are highly sensitive to the interphase properties. The results are compared mainly with the self consistent model reported in the literature and good agreements are shown.     

Linear and non-linear nature of the flow of polypropylene filled with ferrite particles: from low to concentrated composites

The flow behavior of a filled suspension consisting of ferrite particles suspended in a polypropylene matrix with and without the addition of a commercial dispersant (Solplus DP310) was studied. The composites were filled with 10, 20, 30, and 40 vol.%. Both capillary and parallel disk rotational flows were employed. On the one hand, dynamic results confirm general trends found for highly concentrated systems. The higher is the filler level, the lower is the linear viscoelastic domain. When adding the dispersant agent, it was shown a larger linear viscoelastic domain, lower moduli values and thus, lower viscosity. Also, the critical strain, G′ and G′′ showed a power law dependency on the volume fraction. On the other hand, the capillary results showed no dependency of the flow properties on the die. Thus, no slip of the suspension at the wall was observed. Actually, this experimental finding elucidated that the significant decrease on viscosity produced by the addition of the dispersant agent at 40 vol.% is principally due to lubricant effects and not at all to slip contributions. The results also reveal three distinct flow regimes. Low, moderate, and high shear rates lead to different microstructure under flow.     

The shift of Curie temperature and evolution of ferroelectric domain in ferroelectric crystals

A micromechanics-based thermodynamic model for the phase transition of ferroelectric crystals is developed and, with it, the shift of Curie temperature and evolution of ferroelectric phase upon cooling are examined. This approach differs from the classical phenomenological one in that the evolution of new domain concentration can be predicted. We start out by formulating the Gibbs free energy of a generic, two-phase crystal consisting of the parent paraelectric phase and the transformed ferroelectric phase, at a given level of temperature, stress, and electric field. The thermodynamic driving force for domain growth is then derived and, together with the resistance force, a kinetic equation is established. The derived driving force is found to arise from three different sources of Gibbs free energy: (i) the interaction energy due to the heterogeneity of electromechanical moduli of the parent and product phases, (ii) the energy dissipation due to spontaneous polarization, and (iii) the self-energy of the dual-phase system due to the existence of polarization strain and electric polarization. For a BaTiO₃ crystal the electromechanical heterogeneity is found to play a rather significant role that seems not to have been recognized before. The derived shift recovers to the Clausius-Clapeyron relation if such heterogeneity disappears. We have examined in detail several factors that affect the shift of Curie temperature, and calculated the evolution of overall polarization and dielectric constant of a BaTiO₃ crystal. The results are found to be consistent with available test data.     

On the thermomechanical coupling in finite strain plasticity theory with non-linear kinematic hardening by means of incremental energy minimization

The thermomechanical coupling in finite strain plasticity theory with non-linear kinematic hardening is analyzed within the present paper. This coupling is of utmost importance in many applications, e.g., in those showing low cycle fatigue (LCF) under large strain amplitudes. Since the by now classical thermomechanical coupling originally proposed by Taylor and Quinney cannot be used directly in case of kinematic hardening, the change in heat as a result of plastic deformation is computed by applying the first law of thermodynamics. Based on this balance law, together with a finite strain plasticity model, a novel variationally consistent method is elaborated. Within this method and following Stainier and Ortiz (2010), all unknown variables are jointly and conveniently computed by minimizing an incrementally defined potential. In sharp contrast to previously published works, the evolution equations are a priori enforced by employing a suitable parameterization of the flow rule and the evolution equations. The advantages of this parameterization are, at least, twofold. First, it leads eventually to an unconstrained stationarity problem which can be directly applied to any yield function being positively homogeneous of degree one, i.e., the approach shows a broad range of application. Secondly, the parameterization provides enough flexibility even for a broad range of non-associative models such as kinematic hardening of Armstrong–Frederick-type. Different to Stainier and Ortiz (2010), the continuous variational problem is approximated by a standard, fully-implicit time integration. The applicability of the resulting numerical implementation is finally demonstrated by analyzing the thermodynamically coupled response for a loading cycle.     

On the lower bounds for critical loads under large deformations in non-linear hyperelastic composites with imperfect interlaminar adhesion

The present paper investigates a mechanism of compressive fracture for heterogeneous incompressible non-linear materials with special kinds of defects of interfacial adhesion under large deformations. The analysis finds the lower bounds for the critical load. In order to calculate the bounds, the problem of the internal instability is considered within the scope of the exact statement based on the application of the model of a piecewise-homogeneous medium and the equations of the 3-D stability theory. The solution of the 3-D problem is found for the most general case accounting for large deformations and the biaxiality of compressive loads. The characteristic determinants are derived for the first four modes, which are more commonly observed. Special attention is given to the calculation of critical loads for hyperelastic layers described by a simplified version of Mooney's potential, namely the neo-Hookean potential.