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

Previous section: TITLE PAGE
Next section: REFERENCES

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.


Previous section: TITLE PAGE
Next section: REFERENCES