Applied Physics Letters, 14 January 2008
Appl. Phys. Lett. 92, 026101 (2008) (2 pages)
©2008 American Institute of Physics. All rights reserved. Rightslink - Permissions for ReusePermissions for ReuseAbout Rightslink

Up: Issue Table of Contents
Go to: Previous Article | Next Article
Other formats: HTML (smaller files) | PDF (117 kB)

Comment on “Method to analyze electromechanical stability of dielectric elastomers” [Appl. Phys. Lett. 91, 061921 (2007)]

Andrew N. Norris(a)

Mechanical and Aerospace Engineering, Rutgers University, Piscataway, New Jersey 08854, USA

(Received: 18 September 2007; accepted: 17 December 2007; published online: 14 January 2008)

The model of Zhao and Suo can be readily generalized to predict the critical breakdown electric field Ec of elastomers with arbitrary elastic strain energy function. An explicit expression for Ec is presented for elastomeric thin films under biaxial strain and comparisons are made with experimental data using a two term Ogden rubber elasticity model. Simplified results for uniaxial and for equibiaxial stress provide further insight into the electromechanical stability model. ©2008 American Institute of Physics


Contents

The paper of Zhao and Suo1 describes a fully nonlinear electromechanical model for the phenomenon of electrical breakdown in thin elastomers. The purpose of this comment is to point out some analytical simplifications which provide further insight into their model and to provide explicit formulas useful for elastomer design.

The results here stem from the observation that the determinant of the Hessian H of Eq. (4) in Ref. 1 may be factored, leading to semiexplicit formulas for the critical values of the electrical and mechanical parameters. It may be checked that the determinant reduces to a quadratic in z,

det  <b>H</b> = ((<i>µ</i><sup>2</sup><i>epsilon</i><sup>−1</sup>)/(<i>lambda</i><sub>1</sub><sup>8</sup><i>lambda</i><sub>2</sub><sup>8</sup>))[5 + 3(<i>lambda</i><sub>1</sub><sup>2</sup> + <i>lambda</i><sub>2</sub><sup>2</sup>)<i>lambda</i><sub>1</sub><sup>2</sup><i>lambda</i><sub>2</sub><sup>2</sup> + <i>lambda</i><sub>1</sub><sup>6</sup><i>lambda</i><sub>2</sub><sup>6</sup> + [2 − (<i>lambda</i><sub>1</sub><sup>2</sup> + <i>lambda</i><sub>2</sub><sup>2</sup>)<i>lambda</i><sub>1</sub><sup>2</sup><i>lambda</i><sub>2</sub><sup>2</sup>]<i>z</i> − 3<i>z</i><sup>2</sup>],1

where the nondimensional parameter z is

<i>z</i> = (<i>µ</i> <i>epsilon</i>)<sup>−1</sup><i>D</i>-tilde<sup>2</sup> = <i>µ</i><sup>−1</sup><i>epsilon</i> <i>E</i>-tilde<sup>2</sup><i>lambda</i><sub>1</sub><sup>4</sup><i>lambda</i><sub>2</sub><sup>4</sup> = <i>µ</i><sup>−1</sup><i>epsilon</i> <i>E</i><sup>2</sup><i>lambda</i><sub>1</sub><sup>2</sup><i>lambda</i><sub>2</sub><sup>2</sup>,2

and all other notation follows.1 The roots of the quadratic are real and of opposite sign, so there is a unique positive value of z at which the Hessian is no longer positive definite. It turns out that the same structure of the Hessian is retained for free energy of the form

<i>W</i>(<i>lambda</i><sub>1</sub>,<i>lambda</i><sub>2</sub>,<i>D</i>-tilde) = <i>U</i>(<i>lambda</i><sub>1</sub>,<i>lambda</i><sub>2</sub>)+(<i>D</i>-tilde<sup>2</sup>/(2 <i>epsilon</i>))<i>lambda</i><sub>1</sub><sup>−2</sup><i>lambda</i><sub>2</sub><sup>−2</sup>,3

where U(lambda1,lambda2)=psi(lambda1,lambda2,lambda<sub>1</sub><sup>−1</sup>lambda<sub>2</sub><sup>−1</sup>). An equation similar to Eq. (1) is obtained, and taking the single positive root shows that the critical value of the electric field satisfies

<i>epsilon</i> <i>E</i><sub><i>c</i></sub><sup>2</sup> = (1/6)[4 <i>lambda</i><sub>1</sub><i>lambda</i><sub>2</sub><i>U</i><sub>12</sub> − <i>lambda</i><sub>1</sub><sup>2</sup><i>U</i><sub>11</sub> − <i>lambda</i><sub>2</sub><sup>2</sup><i>U</i><sub>22</sub> + sqrt((<i>lambda</i><sub>1</sub><sup>2</sup><i>U</i><sub>11</sub> + <i>lambda</i><sub>2</sub><sup>2</sup><i>U</i><sub>22</sub> − 4 <i>lambda</i><sub>1</sub><i>lambda</i><sub>2</sub><i>U</i><sub>12</sub>)<sup>2</sup> + 12 <i>lambda</i><sub>1</sub><sup>2</sup><i>lambda</i><sub>2</sub><sup>2</sup>(<i>U</i><sub>11</sub><i>U</i><sub>22</sub> − <i>U</i><sub>12</sub><sup>2</sup>))],4

where Uij=[partial-derivative]2U/[partial-derivative]lambdai[partial-derivative]lambdaj. If the stretches lambda1 and lambda2 are prescribed, then Eq. (4) is sufficient to estimate the critical field strength. Otherwise, if the nominal stresses s1 and s2 are prescribed then the stretches are determined from

<i>s</i><sub><i>j</i></sub> = <i>U</i><sub><i>j</i></sub> − <i>lambda</i><sub><i>j</i></sub><sup>−1</sup><i>epsilon</i> <i>E</i><sub><i>c</i></sub><sup>2</sup>,  <i>j</i> = 1,2.5

Under equibiaxial strain lambda1=lambda2=lambda, Eq. (4) becomes

<i>epsilon</i> <i>E</i><sub><i>c</i></sub><sup>2</sup> = ((<i>lambda</i><sup>2</sup>)/3)(<i>U</i><sub>11</sub> + <i>U</i><sub>12</sub>),6

where the critical value of the stress s1=s2=s is

<i>U</i><sub>1</sub>−(1/(3 <i>lambda</i>))(<i>U</i><sub>11</sub> + <i>U</i><sub>12</sub>) = <i>s</i>.7

Consider the Ogden model for rubber2

<i>psi</i>(<i>lambda</i><sub>1</sub>,<i>lambda</i><sub>2</sub>,<i>lambda</i><sub>3</sub>) = [summation]<sub><i>p</i> = 1</sub><sup><i>N</i></sup>((<i>µ</i><sub><i>p</i></sub>)/(<i>alpha</i><sub><i>p</i></sub>))(<i>lambda</i><sub>1</sub><sup><i>alpha</i><sub><i>p</i></sub></sup> + <i>lambda</i><sub>2</sub><sup><i>alpha</i><sub><i>p</i></sub></sup> + <i>lambda</i><sub>3</sub><sup><i>alpha</i><sub><i>p</i></sub></sup>),8

for which the critical electrical field strength is

<i>epsilon</i> <i>E</i><sub><i>c</i></sub><sup>2</sup> = (1/3)[summation]<sub><i>p</i> = 1</sub><sup><i>N</i></sup><i>µ</i><sub><i>p</i></sub>[(<i>alpha</i><sub><i>p</i></sub> − 1)<i>lambda</i><sup><i>alpha</i><sub><i>p</i></sub></sup> + (2 <i>alpha</i><sub><i>p</i></sub> + 1)<i>lambda</i><sup>−2 <i>alpha</i><sub><i>p</i></sub></sup>].9

If the stress is prescribed then lambda is given by

(1/(3 <i>lambda</i>))[summation]<sub><i>p</i> = 1</sub><sup><i>N</i></sup><i>µ</i><sub><i>p</i></sub>[(4 − <i>alpha</i><sub><i>p</i></sub>)<i>lambda</i><sup><i>alpha</i><sub><i>p</i></sub></sup> − (4 + 2 <i>alpha</i><sub><i>p</i></sub>)<i>lambda</i><sup>−2 <i>alpha</i><sub><i>p</i></sub></sup>] = <i>s</i>.10

These parameterize the critical electrical and mechanical fields in terms of lambda.

Values of the critical breakdown voltage for the elastomer VHB 4905/4910 have been reported by Plante and Dubowsky3 and by Kofod et al.4 Assuming the Ogden model with N=2, Plante and Dubowsky3 measured values of alpha1=1.445(1.450), alpha2=4.248(8.360), µ1=43  560(112,200)  Pa, µ2=117.4(0.1045)  Pa for elastomer films of initial thickness L3=1.5  mm at low (high) stretch rates. Using these values, the critical breakdown voltage Vc=L3<i>E</i>-tildec predicted by Eq. (9) is compared with the data of Refs. 3,4 in Fig. 1. The material dielectric constant was chosen as epsilond=12 to fit the curves with the data, where epsilon=epsilondepsilon0 and epsilon0=8.85×10−12  F/m is the free space permittivity. The agreement is reasonable, given that the experiments were not performed in a state of pure equibiaxial stress. Unlike previous theories, e.g., Ref. 5 the Zhao and Suo model provides an estimate of the breakdown strength that takes into account nonlinear electromechanical effects.

Figure 1.

Some useful explicit results can be determined for the one term Ogden model (N=1,alpha1,µ1,-->alpha,µ). Under equibiaxial stress the critical stretch satisfies lambda>=lambdac where lambdac=[(4+2alpha)/(4−alpha)]1/(3alpha) is the s=0 value. This obviously requires that alpha<4. The critical field Ec has a unique minimum at lambda0=[2(2alpha+1)/(alpha−1)]1/(3alpha) if alpha>1. Zhao and Suo1 considered alpha=2, for which lambdac[approximate]1.26, lambda0[approximate]1.47 and the minimum value of sqrt((<i>epsilon</i> /<i>µ</i>))Ec is 1.038.

Finally, we note that the neo-Hookean constitutive model of Zhao and Suo1 is apparently unique among the N=1 Ogden models in that it yields a simple formula for uniaxial stress. Thus, Eq. (5) with N=1, alpha=2 for j=2 and s2=0 yields the relation lambda<sub>1</sub><sup>2</sup>=3lambda<sub>2</sub><sup>2</sup>/(lambda<sub>2</sub><sup>6</sup>−1) between the stretches. Hence, we can parameterize the critical values in terms of 1<lambda2<=lambdac[approximate]1.26:

sqrt(((<i>epsilon</i>)/(<i>µ</i>)))<i>E</i><sub><i>c</i></sub> = ((2/3)<i>lambda</i><sub>2</sub><sup>2</sup>+(1/3)<i>lambda</i><sub>2</sub><sup>−4</sup>)<sup>1/2</sup>,11<i>a</i>

(<i>s</i><sub>1</sub>/(<i>µ</i>)) = ((<i>lambda</i><sub>2</sub>)/(sqrt(3)))(((4 − <i>lambda</i><sub>2</sub><sup>6</sup>))/(sqrt( <i>lambda</i><sub>2</sub><sup>6</sup> − 1))).11<i>b</i>

In this case, Ec is a monotonically decreasing function of the stress s1, and sqrt((<i>epsilon</i> /<i>µ</i>))Ec-->1 in the limit of large uniaxial stress. Figure 3(b) in Ref. 1 indicates that this is the smallest achievable value of the critical electric field strength. Generalization of the formulas (10) to alpha[not-equal]2 is possible but far more complicated.

In summary, the model of Zhao and Suo readily generalizes to arbitrary elastic strain energy. The explicit results reported here, such as Eq. (4), can be used to compare different elastic constitutive models, and should be helpful in the design of elastomeric actuators.

REFERENCES


  1. X. Zhao and Z. Suo, Appl. Phys. Lett. 91, 061921 (2007). [ISI] first citation in article
  2. R. W. Ogden, Non-Linear Elastic Deformations (Ellis Horwood, Chichester, UK, 1984). first citation in article
  3. J.-S. Plante and S. Dubowsky, Int. J. Solids Struct. 43, 7727 (2006). [Inspec] first citation in article
  4. G. Kofod, R. Kornbluh, R. Pelrine, and P. Sommer–Larsen, in Proceedings of Smart Structures and Materials 2001: Electroactive Polymer Actuators and Devices (EAPAD) (SPIE, San Diego, CA, 2001), Vol. 4329, pp. 141––147. first citation in article
  5. T. Blythe and D. Bloor, Electrical Properties of Polymers (Cambridge University Press, London, 2005). first citation in article

FIGURES


Full figure (17 kB)

Fig. 1. The data show reported critical breakdown voltages as a function of the equibiaxial prestrain lambda for films of VHB 4905/4910 elastomer, from Refs. 3,4. The curves are the predictions of Eq. (9) using the elasticity parameters from Ref. 3 with epsilond=12. First citation in article

FOOTNOTES

aElectronic mail: norris@rutgers.edu.


Up: Issue Table of Contents
Go to: Previous Article | Next Article
Other formats: HTML (smaller files) | PDF (117 kB)