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.     

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.     

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.     

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.     

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.     

Domain switching effect on fracture of piezoelectric solids

Rational design of smart sensors and actuators that consist of piezoelectric solids requires a thorough understanding of the constitutive behavior of this material under mechanical and electrical loading. Domain switching is the cause of significant nonlinearity in the constitutive behavior of piezoelectric solids, which may be enhanced in the presence of cracks. In this paper, the response of piezoelectric solids is formulated by coupling thermal, electrical, and mechanical effects. The corresponding finite element equations are derived and applied in the solution of the piezoelectric center crack problems. The effects of domain switching are evaluated on the near tip stress intensity factors.     

Two-scale micromechanics-based probabilistic modeling of domain switching in ferroelectric ceramics

A two-scale micromechanics model is developed in this paper to analyze domain switching in ferroelectric ceramics, using a probabilistic domain switching criterion based on energetic analysis. The microstructure of ferroelectric ceramics at two distinct length scales, domains and grains, has been carefully analyzed. The interaction at domain level is accounted for by energy minimization theory, while the fluctuation at grain level is analyzed using ellipsoidal two-point correlation function. The model has been implemented by Monte Carlo method, and applied to simulate the electric poling and mechanical depoling of Pb(Zr_xTi_1-x)O₃ (PZT) ceramics across morphotropic phase boundary (MPB). The drastically different switching characteristics of PZT ceramics across MPB has been captured, and good agreement with experiments has been observed. The effects of the transformation strains and spontaneous polarizations are highlighted, confirming the proposition of Li et al. [2005. Domain switching in polycrystalline ferroelectric ceramics. Nature Materials 4, 776–781] that the strain compatibility plays a dominant role in domain switching in ferroelectric ceramics.     

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.     

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.     

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.     

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.     

Finite element simulation of a polycrystalline ferroelectric based on a multidomain single crystal switching model

In this paper, we compute the constitutive behavior of a ferroelectric ceramic by a plane strain finite element model, where each element represents a single grain in the polycrystal. The properties of a grain are described by the microscopic model for switching in multidomain single crystals of ferroelectric materials presented by Huber et al. [J. Mech. Phys. Solids 47 (1999) 1663]. The poling behavior of the polycrystal is obtained by employing the finite element formulation for electromechanical boundary value problems developed by Landis [Int. J. Numer. Meth. Eng. 55 (2002) 613]. In particular, we address the influence of the single grain properties and the interaction between grains, respectively.     

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.     

An I-integral method for crack-tip intensity factor variation due to domain switching in ferroelectric single-crystals

In the present study, an I-integral method is established for solving the crack-tip intensity factors of ferroelectric single-crystals. The I-integral combined with the phase field model is successfully used to investigate crack-tip intensity factor variations due to domain switching in ferroelectricity subjected to electromechanical loadings, which exhibits several advantages over previous methods based on small-scale switching. First, the shape of the switching zone around a crack tip is predicted by the time-dependent Ginzburg–Landau equation, which does not require preset energy-based switching criterion. Second, the I-integral can directly solve the crack-tip intensity factors and decouple the crack-tip intensity factors of different modes based on superimposing an auxiliary state onto an actual state. Third, the I-integral is area-independent, namely, the I-integral is not affected by the integral area size, the polarization distributions, or domain walls. This makes the I-integral applicable to large-scale domain switching. To this end, the electro-elastic field intensity factors of an impermeable crack in PbTiO₃ ferroelectric single crystals are evaluated under electrical, mechanical, and combined loading. The intensity factors obtained by the I-integral agree well with those obtained by the extrapolation technique. From numerical results, the following conclusions can be drawn with respect to fracture behavior of ferroelectrics under large-scale switching. Under displacement controlled mechanical loading, the stress intensity factors (SIFs) decrease monotonically due to the domain switching process, which means a crack tip shielding or effective switching-induced toughening occurs. If an external electric field is applied, the electric displacement intensity factor (EDIF) increases in all cases, i.e., the formed domain patterns enhance the electric crack tip loading. The energy release rate, expressed by the crack-tip J-integral, is reduced by the domain switching in all examples, which underlines the switching-induced-toughening effect. In contrast, under stress controlled load, the SIF evolves due to large-scale switching to a stable value, which is higher than the non-switching initial value, i.e., fracture is promoted in this case.     

基于RBF插值技术的CFD/CSD非线性耦合分析方法研究

基于非结构混合网格的N-S方程求解器和结构柔度影响系数法,发展了一种考虑气动、结构非线性的基于RBF插值技术CFD/CSD耦合分析方法,适用于解决现代大展弦比飞机的非线性静气动弹性问题。该方法采用时间相关法（即求解非定常方程组,用长时间的渐近解趋于定常状态）求解静气弹分析时的定常流动。考虑大展弦比飞机结构变形问题为大变形小应力问题,在利用柔度系数法求解结构方程时,假设每次求解结构方程时应力与应变为线性关系,整体静气弹分析过程为非线性关系,因此每次求解结构方程时要更新柔度影响系数矩阵。在非定常N-S方程每求解一个时间步耦合一次结构有限元分析,由于结构有限元分析的时间相对于气动分析时间是很短的,所以这种方法实际上近似使用了一次求解非定常气动力的时间完成了整个静气动弹性分析的过程。对于气动网格与结构有限元网格不一致性,本文采用径向基函数（RBF）插值方法中的TPS方法进行结构弹性变形和气动载荷插值,采用虚功原理完成气动载荷数据交换。为了节省气弹分析时间,采用动网格方法对气动网格进行更新,本文基于RBF插值方法发展一种适用于混合网格（四面体、三棱柱、金字塔和六面体）变形的动网格方法,可以保证附面层网格的质量与分布从而准确模拟其流动。利用该方法对M6机翼、DLR-F6翼身组合体和某大型客机机翼进行了静气动弹性特性分析,结果验证了本文开发的非线性CFD/CSD耦合分析方法的可行性、精确性和高效性。     

A micromechanics method to study the effect of domain switching on fracture behavior of polycrystalline ferroelectric ceramics

ＩｎｔｒｏｄｕｃｔｉｏｎＲｅｃｅｎｔｌｙ ,ｔｈｅｆｅｒｒｏｅｌｅｃｔｒｉｃｃｅｒａｍｉｃｓｈａｓｓｕｃｈｅｘｃｅｌｌｅｎｔｃｈａｒａｃｔｅｒｉｓｔｉｃｓｏｆｐｉｅｚｏｅｌｅｃｔｒｉｃｉｔｙａｎｄｐｙｒｏｅｌｅｃｔｒｉｃｉｔｙｅｔｃ .ｔｈａｔｉｔｂｅｃｏｍｅｓｏｎｅｏｆｔｈｅｍｏｓｔｉｍｐｏｒｔａｎｔｆｕｎｃｔｉｏｎａｌｍａｔｅｒｉａｌｓ.Ｆｏｒｉｎｓｔａｎｃｅ ,ｔｈｅｗｉｄｅｌｙａｐｐｌｉｅｄｓｅｎｓｏｒｓ,ｔｒａｎｓｄｕｃｅｒｓａｎｄａｃｔｕａｔｏｒｓｅｔｃ .ａｒｅｍａｄｅｏｆｔｈｅｆｅｒｒｏｅ…     

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.