# Dynamic Electromechanical Response of a Viscoelastic Dielectric Elastomer under Cycle Electric Loads

Authors: Junjie Sheng; Yuqing Zhang
Journal: International Journal of Polymer Science
 (2018)

---

## Abstract

Dielectric elastomer (DE) is able to produce large electromechanical deformation which is time-dependent due to the viscoelasticity. In the current study, a thermodynamic model is set up to characterize the influence of viscoelasticity on the electromechanical and dynamic response of a viscoelastic DE. The time-dependent dynamic deformation, the hysteresis, and the dynamic stability undergoing viscoelastic dissipative processes are investigated. The results show that the electromechanical stability has strong frequency dependence; the viscoelastic DE can attain a larger stretch in the dynamic response than the quasistatic actuation. Furthermore, with the decreasing frequency of the applied electric load, the viscoelastic DE system will present dynamic stability evolution from an aperiodic motion to the quasiperiodic motion. The DE system may also experience a stability evolution from a single cycle motion to multicycle motion with the increasing relaxation times. The value and variation trend of the amplitude of the stretch are highly dependent on the excitation frequency and the relaxation time.


---

## Body

## 1. Introduction

Due to their fast response time, soft, lightweight, low-cost, and high energy density, dielectric elastomers (DEs) have been developed to use in high performance applications as artificial muscles, Braille displays, life-like robots, tunable lens, and power generators [1–6]. It consists of a soft elastomeric membrane sandwiched between two compliant electrodes on both sides. Most of the existing studies on DEs have focused on quasistatic deformation [5–8], neglecting the effect of inertia and viscoelasticity.However, to perform as an electromechanical actuator, a Dielectric elastomer (DE) is often subject to transient, time-dependent forces and voltages [9–11] and is mostly expected to deform at high frequencies in applications, where inertia and viscoelasticity can play a significant role in the dynamic application. Applications exploiting the dynamic behavior of DEs have long been realized through experiments including vibrotactile display for mobile applications [12], frequency tuning [13], pumps [14], and acoustic actuator [15]. Recently, researches have been also carried out on modeling the nonlinear vibrations of hyperelastic DE membranes [16–19]. Zhu et al. [16] studied the resonant behavior of a prestretched membrane of a DE and subsequently analyzed the nonlinear oscillations of a DE balloon [17]. Based on simple geometrical and spherical capacitor assumptions, Yong et al. [18] investigated the dynamics of a thick-walled DE spherical shell. Li et al. [19] analyzed the electromechanical and dynamic analyses of tunable pure-shear DE-based resonator and identified the safe operation range for failure prevention while actuating the resonator. Xu et al. [20] obtained an analytical model for the DE by the Euler-Lagrange equation to study the dynamic analysis of a DE with stretching deformation. Jia et al. [21] investigated the response time and dynamic range for a DE actuator. Li et al. [22] presented an analysis of the nonlinear dynamics of a DE as electromechanical resonator (DEER) configured as a pure-shear actuator. An investigation in our previous paper [23] reported the nonlinear dynamic characteristics of a DE membrane undergoing in-plane deformation. Wang et al. [24] investigated the strain-stiffening effect on the nonlinear vibration of a circular DE membrane subjected to electromechanical load.Although the above studies have attempted to model the nonlinear dynamical behavior of the DE to probe its time-dependent performance, they have neglected the viscoelasticity of DE. The time-dependency of a DE can cause dissipation in the system and significantly affect its dynamic performance and coupling efficiency [25, 26]. In particular, experiments have shown that viscoelasticity can significantly influence the electromechanical transduction and its application [4, 9, 10, 25–27].Recently, a few studies have been developed concerning the dynamic performance of viscoelastic DE [28–32]. Zhang et al. [28] studied the coupled nonlinear oscillation and stability evolution of viscoelastic DE under nonequibiaxial tensile forces by utilizing the method of virtual work, developed an analytical model to characterize the dynamic performance of a homogeneously deformed viscoelastic DE under the conditions of equal-biaxial force, uniaxial force, and pure-shear state by using the Euler-Lagrange equation [29], and predicted the damping effect on the dynamic performance of DE based on the standard linear solid rheological model [30]. Wang et al. [31] presented a model describing the nonlinear dynamic visco-hyperelastic behaviors of DE and explained the material’s dynamic energy dissipation mechanism. However, the effect of viscoelasticity and relaxation time on hysteresis and stability evolution of the DE was neglected in previously reported papers. Zhang et al. [32] investigated the static and dynamic performance of a hinge configuration with integrated dielectric elastomers. Our previous paper [33] also investigated the effect of temperature on the dynamic electromechanical performance of viscoelastic dielectric.Based on the theory of nonequilibrium thermodynamics [34–36] and the nonlinear vibration modeling [16–19], we aim to characterize the electromechanical and dynamic response of viscoelastic dielectric elastomer, to predict how viscoelasticity affect its dynamic performance and hysteresis process by comparing with the quasistatic response, and to present a physical interpretation on the instability and stability evolution coupled by viscoelasticity, relaxation time, and dynamic deformation.In this study, the influence of the viscoelasticity effects on the electromechanical and dynamic characteristics of a DE membrane is considered by setting up a thermodynamic model. With the proposed model, the dynamic responses of a DE film subject to a cyclic electric field are investigated and compared with the quasistatic response. Then, the dynamic oscillation, phase diagrams, and Poincaré maps of the viscoelastic DE are studied. A detailed analysis shows the influence of the frequency, viscoelasticity, and relaxation time on the dynamic stability evolution and the hysteresis.


## 2. Thermodynamic Modeling of a Viscoelastic DE

Considering a piece of viscoelastic dielectric, it deforms when subject to both in-plane biaxial force and a voltage-induced force through the thickness. Figure1 illustrates a membrane of a DE, sandwiched between two compliant electrodes. The DE deforms from its original configuration L1×L2×L3 to the current configuration l1×l2×l3 in terms of an equal biaxial expansion. Here, we define λii=1,2,3=li/Li as the stretch ratio. The stretch λi is homogeneous through the DE and is a function of time. Due to the incompressibility of the material, we express the stretch ratio in the thickness direction as λ3=λ1-1λ2-1.Figure 1
Schematics of a viscoelastic dielectric elastomer subject biaxial force in the plane and voltage through the thickness applied via stretchable electrodes.The nominal stresses are defined ass1=P1/L2L3 and s2=P2/L1L3. The nominal electric field is E~=Φ/L3 and the nominal electric displacement is defined by D~=Q/L1L2, the true electric field by E=Φ/l3, and the true electric displacement by D=Q/l1l2. We adopt the model of ideal DE, assuming that the permittivity of the elastomer ε is taken to be constant. That is, the true electric displacement is taken to be linear in the true electric field, D=εE [19, 33, 34]. Based on the assumption of the DE membrane as a parallel capacitor, we obtain the relation between the charge and the voltage as(1)Q=ΦεL1L2L3λ12λ22.When the membrane is subject to forces and voltage, all three quantitiesΦ, λ1, and λ2 can vary, so that the variation of the charge is(2)δQ=εL1L2L3λ12λ22δΦ+ΦεL1L2L32λ1λ22δλ1+2λ2λ12δλ2.When the dimensions of the DE varies slightly byδλ1 and δλ2, the tensile forces do work of P1L1δλ1+P2L2δλ2. Associated with a small amount of charge δQ moving through the conducting wire, the applied voltage does a work of ΦδQ. During actuation, the inertia forces in each material element along the x-direction and y-direction are ρL2L3x2d2λ1/dt2 and ρL1L3y2d2λ2/dt2, respectively [19]. The damping forces along the x-direction and y-direction are cxdλ1/dt and cydλ2/dt, respectively, which are linear with respect to the deformation velocity in these directions, where ρ is the density of the elastomer and c is the damping coefficient [20]. By its displacement derivative along the x-direction and y-direction, it is not difficult to obtain the inertial force does work and the damping force does work as(3)ρL2L3d2λ1dt2δλ1∫0L1x2dx=L13ρL2L33d2λ1dt2δλ1ρL1L3d2λ2dt2δλ2∫0L2y2dy=L23ρL1L33d2λ2dt2δλ2cdλ1dtδλ1∫0L1xdx=12cL12δλ1dλ1dtcdλ2dtδλ2∫0L2ydx=12cL22δλ2dλ2dt.Thermodynamics requires that, for arbitrary variation of the system, the variation of the free energy of the membrane should equal the work done by the voltage, the prestress, the inertia force, and the damping force; namely,(4)L1L2L3δW=ΦδQ+P1L1δλ1+P2L2δλ2-L13ρL2L33d2λ1dt2δλ1-L23ρL1L33d2λ2dt2δλ2-12cL12δλ1dλ1dt-12cL22δλ2dλ2dt.Experiments have verified that the electromechanical responses of DEs are highly rate-dependent, which implies that the deformation and actuation of DEs are highly relying on the rates of mechanical or electrical activation [9, 10]. This rate-dependence is mainly induced by the viscoelasticity of the elastomeric polymer matrix and may consequently influence the electromechanical actuation [27]. Viscoelastic relaxation may be represented by a rheological model of springs and dashpots [34–36]. As the DE exhibits elastic deformation and inelastic deformation, the latter of which is time-dependent, we first model this material by assuming that it is composed of two molecular chain networks, A and B, as sketched in Figure 2. The network A is an ideal hyperelastic chain and deforms reversibly, while the network B relaxes in time and dissipates energy. The viscous deformation is represented by a dashpot.Figure 2
A Schematic illustration of a viscoelastic model for dielectric elastomer. It consists of two parallel elements: one element is a spring A and the other is a spring B which connected a dashpot in series.For the network represented by the spring and the dashpot at the bottom, the stretches(λ1,λ2) are due to both the spring and the dashpot; we adopt the well-established multiplication rule that [35–37](5)λ1=λ1eξ1,λ2=λ2eξ2,where (λ1e,λ2e) are the stretches due to the bottom spring and (ξ1,ξ2) are stretches due to the dashpot.To account for the effect of strain-stiffening, we represent both springs by using the Gent model [36–38]. The free-energy function of the elastomer is the sum of the contributions from the two springs and can be written as(6)W=-μAJA2ln⁡1-λ12+λ22+λ1-2λ2-2-3JA-μBJB2ln⁡1-λ12ξ1-2+λ22ξ2-2+λ1-2λ2-2ξ12ξ22-3JB+ε2ΦL32λ12λ22,where μA and μB are the shear moduli of the two springs and JA and JA are the extension limits.A combination of (2) and (4) gives that(7)L1L2L3δW=Φ2εL1L2L32λ1λ22δλ1+2λ12λ2δλ2+ΦεL1L2L3λ12λ22δΦ+P1L1δλ1+P2L2δλ2-L13ρL2L33d2λ1dt2δλ1-L23ρL1L33d2λ2dt2δλ2-12cL12δλ1dλ1dt-12cL22δλ2dλ2dt.We assume that the system is in mechanical and electrostatic equilibrium, and we obtain from the standard calculus of variation in (7) that(8)∂W∂λ1=2εΦL32λ1λ22+P1L2L3-L12ρ3d2λ1dt2-12cL1L2L3dλ1dt,(9)∂W∂λ2=2εΦL32λ12λ2+P2L1L3-L22ρ3d2λ2dt2-12cL2L1L3dλ2dt,(10)D~=∂W∂E~=ΦεL3λ12λ22=εE~λ12λ22.Inserting (6) into (8) and (9), we obtain that(11)L12ρ3d2λ1dt2+μAλ1-λ1-3λ2-21-λ12+λ22+λ1-2λ2-2-3/JA+μBλ1ξ1-2-λ1-3λ2-2ξ12ξ221-λ12ξ1-2+λ22ξ2-2+λ1-2λ2-2ξ12ξ22-3/JB-εΦL32λ1λ22-P1L2L3+12cL1L2L3dλ1dt=0,L22ρ3d2λ2dt2+μAλ2-λ1-2λ2-31-λ12+λ22+λ1-2λ2-2-3/JA+μBλ2ξ2-2-λ1-2λ2-3ξ12ξ221-λ12ξ1-2+λ22ξ2-2+λ1-2λ2-2ξ12ξ22-3/JB-εΦL32λ12λ2-P2L2L3+12cL2L1L3dλ2dt=0.To simplify mathematical expressions, we introduce the following dimensionless field variabless1=P1/L2L3μ, E-=Φ/L3ε/μ=E~/ε/μ, s1=P2/L1L3μ, and D-=Q/L1L2εμ=D~/εμ.While the method described here is applicable to general cases, in the following discussion, we will look at a special case when the dielectric film is under equibiaxial stress, such thatL1=L2=L, λ1=λ2=λ, ξ1=ξ2=ξ, and s1=s2=s. The symmetry of the problem reduces the governing equations (11) to(12)ρL23μd2λdt2+χλ-λ-51-2λ2+λ-4-3/JA+1-χλξ-2-λ-5ξ41-2λ2ξ-2+λ4ξ4-3/JB-εμΦL32λ3-s+12cμL3dλdT=0,where μ=μA+μB and χ=μA/μ are the ratio between the equilibrium and instantaneous moduli. The rate of deformation in the dashpot is described by ξ-1dξ/dt. We relate the rate of deformation in the dashpot to the stress on the dashpot and write as [36](13)dξdt=ξμB6ηλ2ξ-2-λ-4ξ41-2λ2ξ-2+λ-4ξ4-3/JB,where η is the viscosity of the dashpot.We can see that the qualitative properties of the solutions of (12) and (13) depend on the applied loads and viscoelasticity. Equations (12) and (13) give the equations of motion for a viscoelastic DE, which we use in the following analysis to study the quasistatic response and the dynamic characteristics of a viscoelastic DE.


## 3. Simulation Results and Discussion

### 3.1. Quasistatic Response

When the stresss and voltage Φ are quasistatic, the membrane may reach a state of equilibrium. In this case, both the inertial force and the viscous force are zero; in the state of equilibrium, the equation of motion (12) can be expressed as(14)χλ-λ-51-2λ2+λ-4-3/JA+1-χλξ-2-λ-5ξ41-2λ2ξ-2+λ4ξ4-3/JB-E-2λ3-s=0.The inelastic deformation of the material evolves with relaxation timeη/μB. Here we will use it to normalize time t by introducing a dimensionless time, τ=tμB/η. Using the dimensionless time, we write the kinetic equation (13) as(15)χλ-λ-51-2λ2+λ-4-3/JA+1-χλξ-2-λ-5ξ41-2λ2ξ-2+λ4ξ4-3/JB-E-2λ3-s=0.In the following analysis, we takeJA=110 kPa, JB=55 kPa, χ=0.5 [36] as representative values for the most widely used DE.When the viscoelastic DE actuator subject to a cyclic electric load,E-=E-0sin⁡ωτ, but not to forces s=0, with E-0=0.5 to avoid instantaneous instability within the selected frequency range and with ω being the dimensionless frequency, the total stretch λ and the inelastic stretch ξ can be easily evolved by solving (14) and (15). The voltage is applied at time τ=0. The dashpot does not move instantaneously, so that the initial value of the internal variable is λ=ξ=1. The time-dependent stretch, nominal electric field-stretch, and nominal electric field-nominal electric displacement curves of representative applied voltage frequency are presented in Figure 3.Figure 3
A membrane is subject to the applied cyclic nominal electric field. (a), (d), (g), and (j) are the calculated time-dependent total stretch and inelastic stretch. (b), (e), (h), and (k) are the stretch-electric field response. (c), (f), (i), and (l) are the dimensionless nominal electric displacement-nominal electric field curves. Four different dimensionless frequencies are used in the calculation withω=10 for (a), (b), and (c), ω=1 for (d), (e), and (f), ω=0.034 for (g), (h), and (i), ω=0.033 for (j), (k), and (l).
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(k)
(l)It can be seen that, at a relatively high frequency (ω=10) of the applied electric field, the deformation of the viscoelastic DE could keep in phase with the electric field, as shown in Figure 3(a), but the actuation strain (λ) is relatively small (shown with solid curves in Figure 3(a)). This is because that the voltage was varied on a time-scale much faster than the viscoelastic relaxation time. Subjecting to a cyclic load, as the actuation stretch only depends on the magnitude of the applied voltage, the frequency of the stretch doubles that of the electric field in Figure 3(a), and the mean value of stretch may drift away from the initial equilibrium state. The nominal electric field-stretch response (E--λ) in Figure 3(b) and the nominal electric field-nominal electric displacement response (E--D-) in Figure 3(c) show no clear limit cycle, indicating that there will be insignificant hysteresis in the viscoelastic DE for high frequency.When the frequency of the applied electric field decreases to an intermediate frequency, for example,ω=0.1, the magnitude of the actuation strain (λ) becomes much larger, and the amplitude of the inelastic stretch is almost the same as that of λ, showing almost fully inelastic deformation. Limit cycles appear on E--λ and E--D- plots in Figures 3(e) and 3(f). Both plots show hysteresis which indicate that significant amount of energy is dissipated during the electromechanical conversion process.When the DE is actuated by a relatively low frequency voltage forω=0.034, the magnitude of both the actuation stretch and the inelastic stretch attains peak values and there will be most dissipative energy during this cycle, as shown in Figures 3(g), 3(h), and 3(i). The main reason is the geometric nonlinearity introduced by the finite viscoelastic deformation in (14) and (15) according to Hong [34]. When the input electric field is sinusoidal, the output total stretch (λ) and the inelastic stretch (ξ) are not sinusoidal any more. The reason is that the thermoelectromechanically system under a cycle electric load is highly nonlinear and the asymmetry in the shape of the output cyclic deformation λτ at low frequency would lead to significant irreversible energy loss of the system.However, when the frequency decrease to a critical value, forω≤0.033, the viscoelastic DE will suffer the pull-in instability a little while as plotted in Figures 3(j), 3(k), and 3(l), the instability is marked with “×.” When the voltage is applied at a relatively low frequency, the major part of the total deformation is inelastic as shown in Figure 3(j). This result agrees well with experimental observations [9] that the DE more often suffers the pull-in instability at low stretch rates when inelastic deformation is significant.


### 3.2. Dynamic Response

When the DE with stretching deformation oscillates under cycle voltage, the inertia force will contribute to the electromechanical performance. In the following, numerical simulation results of the model are presented to show the vibration and oscillation behavior of the viscoelastic DE.Here we will use the normalized timet by introducing the dimensionless time t~=t/Lρ/3μ, c~=c/2μL3Lρ/3μ and dimensionless relaxation time τ~=η/μB/Lρ/3μ, and we may reduce (12) and (13) to the following form:(16)d2λdt~2+χλ-λ-51-2λ2+λ-4-3/JA+1-χλξ-2-λ-5ξ41-2λ2ξ-2+λ4ξ4-3/JB-E-2λ3-s+c~dλdT=0,dξdt~=16τ~ξλ2ξ-2-λ-4ξ41-2λ2ξ-2+λ-4ξ4-3/JB.To be specific, the following parameters are used in the dynamic examples:JA=110 KPa, JB=55 KPa, χ=0.5, τ~=1, and c~=1 [20, 34–36].The cyclic electric load,E-=E-0sin⁡Ωt~, is applied at time t~=0, with E-0=0.5 to avoid instantaneous instability within the selected frequency range and with Ω being the dimensionless frequency. The dashpot does not move instantaneously, so that the initial value of the internal variable is λ0=ξ0=1 and λ˙0=0. Under these assumptions, when no prestress is applied, s=0, by solving the governing equations (16), the dynamic stability evolution and performance of the viscoelastic elastomer, such as the oscillation diagram, the hysteresis loop, the phase diagram, and the Poincaré map, are simulated.Due to nonlinearity of the system, the dynamic response of the viscoelastic elastomer under time-dependent electric load can be very complicated. When the representative dimensionless frequency varies from 0 to 10, the viscoelastic DE presents oscillations with varying amplitudes shown in Figures4–7.Figure 4
The total stretchλ the viscous stretch ξ (a), the nominal electric field-stretch response (b), the nominal electric field-nominal electric displacement (c), and the phase diagram (d) for the dimensionless frequency Ω=10.Figure 5
The total stretchλ, the viscous stretch ξ (a), the nominal electric field-stretch response (b), the nominal electric field-nominal electric displacement (c), and the phase diagram (d) for the dimensionless frequency Ω=1.
(a)
(b)
(c)
(d)Figure 6
The total stretchλ, the viscous stretch ξ (a), the nominal electric field-stretch response (b), the nominal electric field-nominal electric displacement (c), and the phase diagram (d) for the dimensionless frequency Ω=0.1.
(a)
(b)
(c)
(d)Figure 7
The total stretchλ, the viscous stretch ξ (a), the nominal electric field-stretch response (b), the nominal electric field-nominal electric displacement (c), and the phase diagram (d) for the dimensionless frequency Ω=0.0213.
(a)
(b)
(c)
(d)It can be seen that voltage applied at a relative high frequency (Ω=10) gives a relative small total stretch (λ) and viscous stretch (ξ) in Figure 4(a). At the same time, the nominal electric field-stretch response (E--λ) in Figure 4(b) and the nominal electric field-nominal electric displacement response (E--D-) in Figure 4(c) show no clear limit cycle, indicating that there will be insignificant hysteresis in the viscoelastic DE for high frequency. The phase diagram in Figure 4(d) is not closed for Ω=10 at the outset, emerging as a tangle of interlaced curves, which means that the system is aperiodic and unstable. However, the phase diagram in Figure 4(d-2) forms a closed loop after few cycles, implying a quasiperiodic vibration. The amplitude of the stretch λ is extremely small; in other words, the quasiperiodic vibration is almost static.Figure5(a) shows the time-dependent behavior of the total stretch and the viscous stretch for Ω=11. The stretch increases and stabilizes after few cycles. E--λ, E--D- response and the phase diagram are plotted in Figures 5(b), 5(c), and 5(d). As can be seen from Figures 5(b) and 5(c), the hysteresis loop can be observed, indicating dissipation of energy. The phase paths in Figure 5(d) are presented in a closed loop when the system has achieved stable vibration after few cycles; this result shows that the system undergoes a nonlinear quasiperiodic oscillation.The dynamic behaviors of the viscoelastic DE forΩ=0.1 are shown in Figure 6. When the DE is actuated by a relatively low frequency voltage for Ω=0.1, the magnitude of the total stretch and the viscous stretch are both much larger and there will be significant dissipative energy during this cycle, as shown in Figures 6(a), 6(b), and 6(c), which is similar to the quasistatic response. The phase diagram is closed, indicating a steady state of oscillation, as indicated in Figure 6(d).In addition, the critical value for the frequency to suppress the pull-in instability can be observed by solving (16). The critical value of the frequency is 0.0213. The dynamic performances of the viscoelastic dielectric elastomer for very low frequency, Ω=0.0213, are plotted in Figure 7. It can be seen that the stretch attains a maximum value which is about 1.5 than the critical value of quasistatic response. What is more, the hysteresis loop in the E--λ and E--D- plots becomes much bigger as compared with the high frequency in Figures 4–6 and the quasistatic response in Figure 3. That is to say, the dissipation of energy in the dynamic response is bigger than the quasistatic response. The reason is that the dynamic response of DE system under a cycle electric load is highly nonlinear and the asymmetry in the shape of the output cyclic deformation λt~ at low frequency for dynamic oscillation, and the nonlinear and asymmetry become more and more obviously, which would lead to more significant irreversible energy loss of the system in dynamic response. When the applied voltage frequency decrease to the critical value, for Ω≤0.0213, the viscoelastic De will suffer the pull-in instability after a certain time.Poincaré map can be used to detect the stability transition of the system [20, 23, 28]. The Poincaré maps for the oscillation with the frequency Ω=10, 1, 0.1, and 0.0213 are shown in Figure 8. It can be concluded that there exists strange attractor in the Poincaré map for the oscillation with the frequency Ω=10 in Figure 8(a) for the initiate cycles, which is the representative of an aperiodic vibration motion; with the creep of time, the oscillation of DE the Poincaré map at the frequency of Ω=10 in Figure 8(a-2) shows one point, indicating single cycle oscillation. That is, the DE experiences a dynamic stability evolution from an aperiodic motion to the quasiperiodic motion for a relative high frequency. With the frequency decreasing, the Poincaré maps at the frequency of Ω=1, 0.1, and 0.0213 all have finite points, which mean that the motions are quasiperiodic and stable. As shown in Figure 8(b), the motion of DE at Ω=1 have two periods at least, compared to single cycle motion for frequencies 0.1 (Figure 8(c)) and 0.213 (Figure 8(d)).Figure 8
Poincaré map for the oscillation under a cyclic electric load with different dimensionless frequencies (a)Ω=10, (b) Ω=1, (c) Ω=0.1, and (d) Ω=0.0213.


### 3.3. Effect of Relaxation Time

What is more, the viscoelastic relaxation time of VHB ranges from hundreds of microseconds to hundreds of seconds [39, 40] and affects the creeping deformation ξ of DE; that is, the enhancement of the prestretch (prestretch provides advantages to the electromechanical actuation to improve the electrical breakdown strength, avoiding instability and reducing the geometric thickness) followed by short relaxation approximately doubles that of the same prestretch but with long relaxation [9, 34]. However, the relaxation time effects are seldom analyzed in the scope of viscoelasticity. When no prestress is applied s=0, we solve (16) under three levels of dimensionless viscoelastic relaxation time at the dimensionless frequency of Ω=1, with τ~=1, 10, and 100, and plot the solution in Figure 9.Figure 9
The total stretchesλ, the viscous stretches ξ ((a), (d), and (g)), the phase diagrams ((c), (e), and (h)), and the Poincaré maps ((c), (f), and (i)) at dimensionless frequency of Ω=1 for different dimensionless relaxation times τ~=1, 10, and 100.The oscillation of stretchλ exhibits a strong nonlinearity with the coupled effects of relaxation times and frequencies. For a relative high relaxation time of τ~=10 in Figure 9(d), the stretch creeps with time and requires longer time to achieve a stable motion. For a higher relaxation time of τ~=100, the stretch of λ will drift to a larger value at the stable state of oscillation after many cycles as shown in Figure 9(g-2). It can also be observed that the creeping value of viscous stretch (ξ) reduces gradually with the increasing value of relaxation times (Figures 9(a), 9(d), and 9(g)), and the creep can be almost eliminated in Figure 9(g) when suitable values of relaxation times are applied. When the relaxation time increases, the viscous stretch becomes smaller, that is, because the viscoelastic deformation of the DE fails to keep up with electric field frequency leading to decreasing the viscous stretch.As illustrated in Figure9, the λ-λ˙ phase diagrams (Figures 9(b), 9(e), and 9(h)) are all presented in closed regions indicating that DE system experiences a nonlinear quasiperiodic oscillation and the dynamic oscillation of the DE system is stable. When the relaxation time increases from τ~=1 to τ~=100, the Poincaré maps (Figures 9(c), 9(f), and 9(i)) all show finite points with 2, 3, and 8, respectively. That is to say, the DE system may experiences a stability evolution from a single cycle motion to multicycle motion with the increasing relaxation times for a fixed frequency.


## 3.1. Quasistatic Response

When the stresss and voltage Φ are quasistatic, the membrane may reach a state of equilibrium. In this case, both the inertial force and the viscous force are zero; in the state of equilibrium, the equation of motion (12) can be expressed as(14)χλ-λ-51-2λ2+λ-4-3/JA+1-χλξ-2-λ-5ξ41-2λ2ξ-2+λ4ξ4-3/JB-E-2λ3-s=0.The inelastic deformation of the material evolves with relaxation timeη/μB. Here we will use it to normalize time t by introducing a dimensionless time, τ=tμB/η. Using the dimensionless time, we write the kinetic equation (13) as(15)χλ-λ-51-2λ2+λ-4-3/JA+1-χλξ-2-λ-5ξ41-2λ2ξ-2+λ4ξ4-3/JB-E-2λ3-s=0.In the following analysis, we takeJA=110 kPa, JB=55 kPa, χ=0.5 [36] as representative values for the most widely used DE.When the viscoelastic DE actuator subject to a cyclic electric load,E-=E-0sin⁡ωτ, but not to forces s=0, with E-0=0.5 to avoid instantaneous instability within the selected frequency range and with ω being the dimensionless frequency, the total stretch λ and the inelastic stretch ξ can be easily evolved by solving (14) and (15). The voltage is applied at time τ=0. The dashpot does not move instantaneously, so that the initial value of the internal variable is λ=ξ=1. The time-dependent stretch, nominal electric field-stretch, and nominal electric field-nominal electric displacement curves of representative applied voltage frequency are presented in Figure 3.Figure 3
A membrane is subject to the applied cyclic nominal electric field. (a), (d), (g), and (j) are the calculated time-dependent total stretch and inelastic stretch. (b), (e), (h), and (k) are the stretch-electric field response. (c), (f), (i), and (l) are the dimensionless nominal electric displacement-nominal electric field curves. Four different dimensionless frequencies are used in the calculation withω=10 for (a), (b), and (c), ω=1 for (d), (e), and (f), ω=0.034 for (g), (h), and (i), ω=0.033 for (j), (k), and (l).
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(k)
(l)It can be seen that, at a relatively high frequency (ω=10) of the applied electric field, the deformation of the viscoelastic DE could keep in phase with the electric field, as shown in Figure 3(a), but the actuation strain (λ) is relatively small (shown with solid curves in Figure 3(a)). This is because that the voltage was varied on a time-scale much faster than the viscoelastic relaxation time. Subjecting to a cyclic load, as the actuation stretch only depends on the magnitude of the applied voltage, the frequency of the stretch doubles that of the electric field in Figure 3(a), and the mean value of stretch may drift away from the initial equilibrium state. The nominal electric field-stretch response (E--λ) in Figure 3(b) and the nominal electric field-nominal electric displacement response (E--D-) in Figure 3(c) show no clear limit cycle, indicating that there will be insignificant hysteresis in the viscoelastic DE for high frequency.When the frequency of the applied electric field decreases to an intermediate frequency, for example,ω=0.1, the magnitude of the actuation strain (λ) becomes much larger, and the amplitude of the inelastic stretch is almost the same as that of λ, showing almost fully inelastic deformation. Limit cycles appear on E--λ and E--D- plots in Figures 3(e) and 3(f). Both plots show hysteresis which indicate that significant amount of energy is dissipated during the electromechanical conversion process.When the DE is actuated by a relatively low frequency voltage forω=0.034, the magnitude of both the actuation stretch and the inelastic stretch attains peak values and there will be most dissipative energy during this cycle, as shown in Figures 3(g), 3(h), and 3(i). The main reason is the geometric nonlinearity introduced by the finite viscoelastic deformation in (14) and (15) according to Hong [34]. When the input electric field is sinusoidal, the output total stretch (λ) and the inelastic stretch (ξ) are not sinusoidal any more. The reason is that the thermoelectromechanically system under a cycle electric load is highly nonlinear and the asymmetry in the shape of the output cyclic deformation λτ at low frequency would lead to significant irreversible energy loss of the system.However, when the frequency decrease to a critical value, forω≤0.033, the viscoelastic DE will suffer the pull-in instability a little while as plotted in Figures 3(j), 3(k), and 3(l), the instability is marked with “×.” When the voltage is applied at a relatively low frequency, the major part of the total deformation is inelastic as shown in Figure 3(j). This result agrees well with experimental observations [9] that the DE more often suffers the pull-in instability at low stretch rates when inelastic deformation is significant.


## 3.2. Dynamic Response

When the DE with stretching deformation oscillates under cycle voltage, the inertia force will contribute to the electromechanical performance. In the following, numerical simulation results of the model are presented to show the vibration and oscillation behavior of the viscoelastic DE.Here we will use the normalized timet by introducing the dimensionless time t~=t/Lρ/3μ, c~=c/2μL3Lρ/3μ and dimensionless relaxation time τ~=η/μB/Lρ/3μ, and we may reduce (12) and (13) to the following form:(16)d2λdt~2+χλ-λ-51-2λ2+λ-4-3/JA+1-χλξ-2-λ-5ξ41-2λ2ξ-2+λ4ξ4-3/JB-E-2λ3-s+c~dλdT=0,dξdt~=16τ~ξλ2ξ-2-λ-4ξ41-2λ2ξ-2+λ-4ξ4-3/JB.To be specific, the following parameters are used in the dynamic examples:JA=110 KPa, JB=55 KPa, χ=0.5, τ~=1, and c~=1 [20, 34–36].The cyclic electric load,E-=E-0sin⁡Ωt~, is applied at time t~=0, with E-0=0.5 to avoid instantaneous instability within the selected frequency range and with Ω being the dimensionless frequency. The dashpot does not move instantaneously, so that the initial value of the internal variable is λ0=ξ0=1 and λ˙0=0. Under these assumptions, when no prestress is applied, s=0, by solving the governing equations (16), the dynamic stability evolution and performance of the viscoelastic elastomer, such as the oscillation diagram, the hysteresis loop, the phase diagram, and the Poincaré map, are simulated.Due to nonlinearity of the system, the dynamic response of the viscoelastic elastomer under time-dependent electric load can be very complicated. When the representative dimensionless frequency varies from 0 to 10, the viscoelastic DE presents oscillations with varying amplitudes shown in Figures4–7.Figure 4
The total stretchλ the viscous stretch ξ (a), the nominal electric field-stretch response (b), the nominal electric field-nominal electric displacement (c), and the phase diagram (d) for the dimensionless frequency Ω=10.Figure 5
The total stretchλ, the viscous stretch ξ (a), the nominal electric field-stretch response (b), the nominal electric field-nominal electric displacement (c), and the phase diagram (d) for the dimensionless frequency Ω=1.
(a)
(b)
(c)
(d)Figure 6
The total stretchλ, the viscous stretch ξ (a), the nominal electric field-stretch response (b), the nominal electric field-nominal electric displacement (c), and the phase diagram (d) for the dimensionless frequency Ω=0.1.
(a)
(b)
(c)
(d)Figure 7
The total stretchλ, the viscous stretch ξ (a), the nominal electric field-stretch response (b), the nominal electric field-nominal electric displacement (c), and the phase diagram (d) for the dimensionless frequency Ω=0.0213.
(a)
(b)
(c)
(d)It can be seen that voltage applied at a relative high frequency (Ω=10) gives a relative small total stretch (λ) and viscous stretch (ξ) in Figure 4(a). At the same time, the nominal electric field-stretch response (E--λ) in Figure 4(b) and the nominal electric field-nominal electric displacement response (E--D-) in Figure 4(c) show no clear limit cycle, indicating that there will be insignificant hysteresis in the viscoelastic DE for high frequency. The phase diagram in Figure 4(d) is not closed for Ω=10 at the outset, emerging as a tangle of interlaced curves, which means that the system is aperiodic and unstable. However, the phase diagram in Figure 4(d-2) forms a closed loop after few cycles, implying a quasiperiodic vibration. The amplitude of the stretch λ is extremely small; in other words, the quasiperiodic vibration is almost static.Figure5(a) shows the time-dependent behavior of the total stretch and the viscous stretch for Ω=11. The stretch increases and stabilizes after few cycles. E--λ, E--D- response and the phase diagram are plotted in Figures 5(b), 5(c), and 5(d). As can be seen from Figures 5(b) and 5(c), the hysteresis loop can be observed, indicating dissipation of energy. The phase paths in Figure 5(d) are presented in a closed loop when the system has achieved stable vibration after few cycles; this result shows that the system undergoes a nonlinear quasiperiodic oscillation.The dynamic behaviors of the viscoelastic DE forΩ=0.1 are shown in Figure 6. When the DE is actuated by a relatively low frequency voltage for Ω=0.1, the magnitude of the total stretch and the viscous stretch are both much larger and there will be significant dissipative energy during this cycle, as shown in Figures 6(a), 6(b), and 6(c), which is similar to the quasistatic response. The phase diagram is closed, indicating a steady state of oscillation, as indicated in Figure 6(d).In addition, the critical value for the frequency to suppress the pull-in instability can be observed by solving (16). The critical value of the frequency is 0.0213. The dynamic performances of the viscoelastic dielectric elastomer for very low frequency, Ω=0.0213, are plotted in Figure 7. It can be seen that the stretch attains a maximum value which is about 1.5 than the critical value of quasistatic response. What is more, the hysteresis loop in the E--λ and E--D- plots becomes much bigger as compared with the high frequency in Figures 4–6 and the quasistatic response in Figure 3. That is to say, the dissipation of energy in the dynamic response is bigger than the quasistatic response. The reason is that the dynamic response of DE system under a cycle electric load is highly nonlinear and the asymmetry in the shape of the output cyclic deformation λt~ at low frequency for dynamic oscillation, and the nonlinear and asymmetry become more and more obviously, which would lead to more significant irreversible energy loss of the system in dynamic response. When the applied voltage frequency decrease to the critical value, for Ω≤0.0213, the viscoelastic De will suffer the pull-in instability after a certain time.Poincaré map can be used to detect the stability transition of the system [20, 23, 28]. The Poincaré maps for the oscillation with the frequency Ω=10, 1, 0.1, and 0.0213 are shown in Figure 8. It can be concluded that there exists strange attractor in the Poincaré map for the oscillation with the frequency Ω=10 in Figure 8(a) for the initiate cycles, which is the representative of an aperiodic vibration motion; with the creep of time, the oscillation of DE the Poincaré map at the frequency of Ω=10 in Figure 8(a-2) shows one point, indicating single cycle oscillation. That is, the DE experiences a dynamic stability evolution from an aperiodic motion to the quasiperiodic motion for a relative high frequency. With the frequency decreasing, the Poincaré maps at the frequency of Ω=1, 0.1, and 0.0213 all have finite points, which mean that the motions are quasiperiodic and stable. As shown in Figure 8(b), the motion of DE at Ω=1 have two periods at least, compared to single cycle motion for frequencies 0.1 (Figure 8(c)) and 0.213 (Figure 8(d)).Figure 8
Poincaré map for the oscillation under a cyclic electric load with different dimensionless frequencies (a)Ω=10, (b) Ω=1, (c) Ω=0.1, and (d) Ω=0.0213.


## 3.3. Effect of Relaxation Time

What is more, the viscoelastic relaxation time of VHB ranges from hundreds of microseconds to hundreds of seconds [39, 40] and affects the creeping deformation ξ of DE; that is, the enhancement of the prestretch (prestretch provides advantages to the electromechanical actuation to improve the electrical breakdown strength, avoiding instability and reducing the geometric thickness) followed by short relaxation approximately doubles that of the same prestretch but with long relaxation [9, 34]. However, the relaxation time effects are seldom analyzed in the scope of viscoelasticity. When no prestress is applied s=0, we solve (16) under three levels of dimensionless viscoelastic relaxation time at the dimensionless frequency of Ω=1, with τ~=1, 10, and 100, and plot the solution in Figure 9.Figure 9
The total stretchesλ, the viscous stretches ξ ((a), (d), and (g)), the phase diagrams ((c), (e), and (h)), and the Poincaré maps ((c), (f), and (i)) at dimensionless frequency of Ω=1 for different dimensionless relaxation times τ~=1, 10, and 100.The oscillation of stretchλ exhibits a strong nonlinearity with the coupled effects of relaxation times and frequencies. For a relative high relaxation time of τ~=10 in Figure 9(d), the stretch creeps with time and requires longer time to achieve a stable motion. For a higher relaxation time of τ~=100, the stretch of λ will drift to a larger value at the stable state of oscillation after many cycles as shown in Figure 9(g-2). It can also be observed that the creeping value of viscous stretch (ξ) reduces gradually with the increasing value of relaxation times (Figures 9(a), 9(d), and 9(g)), and the creep can be almost eliminated in Figure 9(g) when suitable values of relaxation times are applied. When the relaxation time increases, the viscous stretch becomes smaller, that is, because the viscoelastic deformation of the DE fails to keep up with electric field frequency leading to decreasing the viscous stretch.As illustrated in Figure9, the λ-λ˙ phase diagrams (Figures 9(b), 9(e), and 9(h)) are all presented in closed regions indicating that DE system experiences a nonlinear quasiperiodic oscillation and the dynamic oscillation of the DE system is stable. When the relaxation time increases from τ~=1 to τ~=100, the Poincaré maps (Figures 9(c), 9(f), and 9(i)) all show finite points with 2, 3, and 8, respectively. That is to say, the DE system may experiences a stability evolution from a single cycle motion to multicycle motion with the increasing relaxation times for a fixed frequency.


## 4. Conclusions

Although the dynamics of DE are widely studied recently, this study gives a further understanding of the effect of viscoelasticity and frequency on dynamic response of DE and can help to guide the design of dynamic and static applications in DE actuators.The dynamic behavior of a viscoelastic DE under a cycle electric field is presented to compare with the quasistatic behavior. The time-dependent dynamic deformation, the hysteresis process, and the dynamic stability undergoing viscoelastic dissipative processes are investigated by using the theory of nonequilibrium thermodynamics. The results show that when the frequency of the applied voltage decreases, whether the quasistatic or the dynamic performance, the total stretch and the viscous stretch both become larger. In addition, the energy dissipated per cycle increases with the decreasing frequency and there exists a critical value of the frequency at which the pull-in instability takes place. When the applied voltage frequency is below a critical value, the viscoelastic DE cannot attain a steady state of oscillation and may have instability. The system of the viscoelastic DE also presents stability transition from an aperiodic motion to the quasiperiodic motion for a relative high frequency and experiences a stability evolution from single cycle motion to multicycle motion with the increasing relaxation times for a fixed frequency. The creeping value of viscous stretch of the dashpot reduces gradually with the increasing value of relaxation times and the creep can be almost eliminated at a suitable value of relaxation time.



---

*Source: 2803631-2018-01-29.xml*