Voltage-controlled non-axisymmetric vibrations
of soft electro-active tubes with strain-stiffening effect 111Dedicated to Yibin Fu, in esteem and friendship
Abstract
Material properties of soft electro-active (SEA) structures are significantly sensitive to external electro-mechanical biasing fields (such as pre-stretch and electric stimuli), which generate remarkable knock-on effects on their dynamic characteristics. In this work, we analyze the electrostatically tunable non-axisymmetric vibrations of an incompressible SEA cylindrical tube under the combination of a radially applied electric voltage and an axial pre-stretch. Following the theory of nonlinear electro-elasticity and the associated linearized theory for superimposed perturbations, we derive the nonlinear static response of the SEA tube to the inhomogeneous biasing fields for the Gent ideal dielectric model. Using the State Space Method, we efficiently obtain the frequency equations for voltage-controlled small-amplitude three-dimensional non-axisymmetric vibrations, covering a wide range of behaviors, from the purely radial breathing mode to torsional modes, axisymmetric longitudinal modes, and prismatic diffuse modes. We also perform an exhaustive numerical analysis to validate the proposed approach compared with the conventional displacement method, as well as to elucidate the influences of the applied voltage, axial pre-stretch, and strain-stiffening effect on the nonlinear static response and vibration behaviors of the SEA tube. The present study clearly indicates that manipulating electro-mechanical biasing fields is a feasible way to tune the small-amplitude vibration characteristics of an SEA tube. The results should benefit experimental work on, and design of, voltage-controlled resonant devices made of SEA tubes.
keywords:
Soft electro-active tube , Non-axisymmetric vibrations , Biasing fields , State Space Method , Strain-stiffening effect , Active resonator1 Introduction
Among the ever-increasing range of smart materials being currently developed, soft electro-active (SEA) materials such as dielectric elastomers can convert or transduce electrical energy to or from mechanical energy (Pelrine et al., 1998). SEA structures demonstrate superior electro-mechanical coupling properties, as they reduce remarkably in thickness and expand in area when exposed to applied electro-mechanical biasing fields, in contrast to typical piezoelectric materials which are too brittle to undergo large deformations. Other excellent features, such as reversible large deformation, rapid response, low weight and low cost, high elastic energy density, and high conversion efficiency, have attracted wide academic and industrial interest and led to various broad practical applications in soft robotics, biomedical devices, flexible electronics, tunable resonators as well as active waveguides, phononic crystals and metamaterials (Carpi et al., 2011; Anderson et al., 2012; Zhao and Wang, 2014; Lu et al., 2020; Zhu et al., 2010b; Zhao et al., 2016; Wang et al., 2020; Chen et al., 2022; Zhao et al., 2023).
Strong nonlinearity and electro-mechanical coupling make it quite difficult to establish a general continuum mechanics framework. Pioneering works on nonlinear theory of electro-elasticity were conducted by Toupin (1956, 1963) more than half a century ago for static and dynamic analyses of finitely deformed elastic dielectrics. Since the 1980s, there have been numerous reformulations of a general nonlinear continuum theory for electro-magneto-mechanical couplings (Maugin, 1988; Eringen and Maugin, 1990), paralleled with the development of various smart materials and structures with wide-ranging applications. Furthermore, the emergence of SEA materials in recent decades has encouraged new interpretations, advancements, and applications of nonlinear electro-elasticity theory (McMeeking and Landis, 2005; Dorfmann and Ogden, 2006; Suo et al., 2008; Liu, 2013; Dorfmann and Ogden, 2014). A nonlinear continuum framework, accounting for the nonlinear interaction between mechanical and electromagnetic fields, as documented in the monograph by Dorfmann and Ogden (2014), has now successfully been applied to the analysis of electro-active and magneto-active materials undergoing significant deformations (Rudykh and Bertoldi, 2013; Xie et al., 2016; Jandron and Henann, 2018; Fu et al., 2018; Psarra et al., 2019; Su et al., 2019, 2020; Wu and Destrade, 2021).
The linearized incremental theory based on the nonlinear electro-elasticity theory (Baumhauer and Tiersten, 1973; Tiersten, 1981; Maugin, 1988; Eringen and Maugin, 1990; Baesu et al., 2003; Dorfmann and Ogden, 2010) is commonly employed to investigate how biasing fields (induced by, for example, prestretch, internal pressure, and electric stimuli) affect the superimposed small-amplitude dynamic properties of SEA structures. Dorfmann and Ogden (2010, 2014) developed a compact version of the linearized incremental theory in both the Lagrangian and updated Lagrangian descriptions to examine the small-amplitude motions superimposed on finite biasing fields, paying particular attention to SEA materials. We recommend the detailed review by Wu et al. (2016) for a comprehensive comparison of different versions of nonlinear electro-elasticity theories and relevant linearized incremental theories, showing ultimately that ostensibly different theories in the literature on this topic are actually equivalent with no substantive differences.
The analysis of electro-mechanical instabilities has significant theoretical and practical implications. Practically speaking, SEA materials and structures could experience multiple failure mechanisms, such as pull-in or snap-through instability (Zhao and Suo, 2007; Su et al., 2018), electric breakdown (Zurlo et al., 2017), macroscopic and microscopic buckling instability (Bertoldi and Gei, 2011; Rudykh et al., 2014; Goshkoderia and Rudykh, 2017), localized necking of SEA membranes (Fu et al., 2018), bending instability of SEA slabs or bilayers (Su et al., 2019, 2020), barrelling axisymmetric instability of SEA tubes (Melnikov and Ogden, 2018), prismatic instability of SEA tubes (Bortot and Shmuel, 2018), bulging instability of SEA tubes or balloons (Lu et al., 2015a; Wang et al., 2017), and pear-shaped bifurcation from SEA spherical balloons (Xie et al., 2016). In particular, Bertoldi and Gei (2011) analytically identified four different instability criteria for multilayered soft dielectrics: (i) loss of positive definiteness of the tangent electro-elastic constitutive operator, (ii) existence of diffuse modes of bifurcation or microscopic instability modes, (iii) loss of strong ellipticity of the homogenized continuum or macroscopic instability modes, and (iv) electric breakdown. In another work, Rudykh et al. (2014) extensively investigated the multiscale instabilities in layered dielectric elastomers and explained the crucial effect of microstructures on the onset of instabilities: (i) macroscopic instabilities predominate when the stiffer phase’s volume fraction is moderate, (ii) interfacial instabilities start to show up at low stiffer phase volume fractions, and (iii) instabilities of a finite scale, comparable to the microstructure size, appear at high volume fractions of the stiffer phase.
Investigations on the small-amplitude dynamic behaviors of smart systems made of SEA materials subject to biasing fields induced by pre-stretch, internal pressure and electric stimuli, also have theoretical as well as practical significance, see the recent review article by Zhao et al. (2023) on vibrations and waves in soft dielectric structures. Based on Dorfmann and Ogden’s linearized incremental theory of nonlinear electro-elasticity (Dorfmann and Ogden, 2010, 2014), considerable efforts have been devoted to studying small-amplitude elastic waves propagating in finitely deformed SEA structures, such as bulk waves in compressible dielectrics (Galich and Rudykh, 2016), surface waves in a deformed SEA half-space (Dorfmann and Ogden, 2010), Rayleigh-Lamb waves in a deformed ideal dielectric plate (Shmuel et al., 2012; Ziser and Shmuel, 2017; Broderick et al., 2020), torsional, axisymmetric, non-axisymmetric and circumferential waves in a pre-stretched SEA tube subject to an axial or a radial electrical biasing field (Shmuel and DeBotton, 2013; Shmuel, 2015; Su et al., 2016; Wu et al., 2017, 2020; Dorfmann and Ogden, 2020), bulk waves in extremely deformed soft auxetic materials (Galich and Rudykh, 2015), shear or longitudinal wave propagation and tunable band gaps in periodic dielectric laminates (Galich and Rudykh, 2017; Chen et al., 2020), and electrostatically tunable band gaps in finitely deformed SEA fiber composites (Shmuel, 2013). In particular, to overcome the difficulty in solving the inhomogeneous biasing fields generated by the application of a radial electric voltage in an SEA tube, Wu et al. (2017) proposed the State Space Method (SSM), which combines the state-space formalism with the approximate laminate technique, to effectively investigate the elastic waves propagating in SEA tubes.
Electro-mechanical biasing fields lead to changes in the effective material characteristics and geometry, and also allow for active tuning of the vibration behaviors to acquire desired operating performance. Thus, tunable SEA resonators have a wide range of potential applications, including the design of tunable SEA loudspeakers for sound generation, the use of adaptive acoustic absorbers for noise reduction, and the development of active and adaptive vibration isolators and dampers that take advantage of viscoelasticity and stiffness tunability. For example, Dubois et al. (2008) examined experimentally and theoretically the voltage controllability of the resonance frequency for SEA polymer membranes, and found that tuning the voltage might reduce the resonance frequency by up to 77% from its initial value. Zhu et al. (2010b) analytically and experimentally found that the natural frequencies of a circular dielectric membrane can be tuned by varying the in-plane pre-stretch, out-of-plane pressure and voltage. Sugimoto et al. (2013) proposed a lightweight push-pull acoustic transducer using dielectric films for sound generation in advanced audio systems, and their experiments showed that push-pull driving can effectively suppress harmonic distortion. Furthermore, Hosoya et al. (2015) constructed, examined, and evaluated a hemispherical breathing mode loudspeaker driven by a dielectric actuator to determine its repeatability, sound pressure, vibration mode profiles, and acoustic radiation patterns. To absorb sound energy, Lu et al. (2015b) developed an electronically tunable duct silencer, which is formed with dielectric membranes and back cavity, and uses external control signals. With the Space State Method, Zhu et al. (2020) studied axisymmetric torsional and longitudinal vibrations in an SEA tube subject to inhomogeneous biasing fields induced by the combined action of axial pre-stretch and radial voltage, and established electrostatically tunable axisymmetric vibration characteristics. By applying alternating voltages with opposite phases to a dielectric actuator, Zhang et al. (2015) proposed a vibration damper to reduce vibration. Sarban et al. (2011) effectively achieved active vibration isolation after fabricating a core-free rolled tubular SEA actuator and studying its dynamic properties. As a biomedical application, Son and Goulbourne (2012) suggested coupling a SEA tube sensor to an artery segment to give it structural support while also keeping track of its local condition data. Zhu et al. (2010a) demonstrated theoretically the tunability of the natural frequency of breathing modes in a dielectric balloon by varying the pressure or voltage. Mao et al. (2019) used the SSM to investigate the feasibility of tuning the three-dimensional (3D) and small-amplitude torsional and spheroidal vibrations in a dielectric spherical balloon by means of varying internal pressure and radial electric voltage. Recently, the SSM was utilized by Cao et al. (2024) to explore the influences of electro-mechanical biasing fields and fluid added mass effect on the linearized axisymmetric vibration of multilayered SEA circular plates in contact with fluid.
The first objective of this paper is to investigate the strain-stiffening effect on axisymmetric torsional and longitudinal vibrations (hereafter abbreviated as T vibrations and L vibrations) of an SEA tube as a continuation of our previous work (Zhu et al., 2020) wherein the strain-stiffening effect was not taken into consideration. The second purpose of this paper is to clarify how the inhomogeneous biasing fields induced by the radial electric voltage and axial pre-stretch (see Figs. 1(a) and 1(b)) and the strain-stiffening effect impact the superimposed non-axisymmetric small-amplitude vibrations in the SEA tube.
The structure of this paper is as follows. Section 2 derives the basic formulations governing the axisymmetric static deformations and induced inhomogeneous biasing fields of an SEA tube characterized by the Gent ideal dielectric model. In Section 3, we employ the state-space formalism combined with an approximate laminate technique to obtain the frequency equations of the non-axisymmetric vibrations and prismatic vibrations of the deformed SEA tube with generalized rigidly supported boundary condition. Numerical calculations are presented in Section 4 to examine the nonlinear static response to the applied radial voltage of the SEA tube, to validate the excellent convergence rate and accuracy of the SSM, and to demonstrate the influences of the applied voltage, axial pre-stretch, and strain-stiffening effect on the axisymmetric and non-axisymmetric vibration characteristics. Finally, we give a conclusive summary in Section 5 and some relevant mathematical derivations are provided in A.
2 Equations of nonlinear electro-elasticity
The basic equations governing the finite electro-elastic deformations of an incompressible soft electro-elastic body are presented first in Sec. 2.1, as established by Dorfmann and Ogden (2006, 2014). Then we specialize the basic equations of nonlinear electro-elasticity to the finite static axisymmetric deformations of an SEA tube subject to a radial electric field and an axial pre-stretch in Sec. 2.2.
2.1 Finite electro-elasticity theory
Here we consider a soft deformable continuous electro-elastic body undergoing a static finite deformation. The underformed stress-free reference configuration at time is denoted by with its boundary and outward unit normal vector denoted as and , respectively. An arbitrary material point in the stress-free reference configuration is labelled by a position vector . At time , stimulated by external electro-mechanical loads, the electro-elastic body deforms to the deformed or current configuration . Naturally, the material point occupies a new position , and the deformation is described by the mapping where is a continuous and twice differentiable vector function. The boundary and outward unit normal of the current configuration are denoted as and , respectively, and is the deformation gradient tensor, where ‘’ is the gradient operator with respect to in the reference configuration . The quantity measures local volume changes, and is equal to one identically for an incompressible material.
The equation of motion, Gauss’s law, and Faraday’s law can be expressed as follows, respectively, within the ‘quasi-electrostatic approximation’ and in the absence of mechanical body forces as well as free charges and electric currents:
(1) |
where ‘’ and ‘’ are the divergence and curl operators with respect to in the current configuration , is the mass density (which is unchanged during the motion due to the material incompressibility), and , and are the total Cauchy stress tensor, Eulerian electric displacement vector and electric field vector in , respectively.
Following nonlinear electro-elasticity theory (Dorfmann and Ogden, 2014), it is convenient to present the nonlinear constitutive relations for incompressible SEA materials in terms of a total energy function or modified free energy function, per unit reference volume in as
(2) |
where is the total nominal stress tensor, and are the nominal electric displacement vector and nominal electric field vector, respectively, which are the Lagrangian counterparts of and , and is a Lagrange multiplier due to the incompressibility constraint. For an incompressible isotropic SEA material (), the energy density function depends on five scalar quantities, for example the following five invariants:
(3) |
where is the right Cauchy-Green deformation tensor, with the superscript T indicating the transpose operator.
Therefore, combination of Eqs. (2) and (3) provides the total Cauchy stress tensor and the Eulerian electric field vector as
(4) |
where is the identity tensor in , is the left Cauchy-Green deformation tensor, and the shorthand notation is adopted throughout this paper.
There is no electric field in the surrounding vacuum when an electric voltage is applied to the surfaces of the SEA body covered with flexible electrodes. Thus, the mechanical and electric boundary conditions to be satisfied on the boundary are written in Eulerian form as
(5) |
where is the applied mechanical traction vector per unit area of and is the free surface charge density on .
2.2 Finite axisymmetric deformations of an SEA tube
The problem of nonlinear axisymmetric deformations of an SEA tube subject to a radial electric field, internal/external pressures, and an axial pre-stretch has previously been discussed by Zhu et al. (2010c); Shmuel and DeBotton (2013); Zhou et al. (2014); Shmuel (2015); Melnikov and Ogden (2016); Su et al. (2016); Wu et al. (2017); Bortot and Shmuel (2018). In this section, we briefly outline for completeness the formulations governing the nonlinear axisymmetric deformations for arbitrary energy function when the SEA tube is subject to a radial voltage combined with an axial pre-stretch (see Fig. 1(b)). We then specialize the results to the Gent ideal dielectric model and obtain explicit expressions for both the nonlinear static response and the radially inhomogeneous biasing fields.
The schematic diagrams of an SEA tube with flexible electrodes before and after activation are illustrated in Figs. 1(a) and 1(b), respectively. For convenience, we use the cylindrical coordinate systems and to describe the undeformed and deformed configurations, respectively. In the undeformed configuration, the inner and outer radii and the thickness of the tube are specified as , , and , respectively, with the tube length denoted by . To realize electro-mechanically tunable vibration characteristics of the SEA tube, a radial electric voltage is applied to the electrodes of the tube, which is also simultaneously subject to a uniform axial pre-stretch . As a result, the inner and outer radii, the thickness and the length of the deformed tube are , , , and , respectively. The inner-to-outer radius ratios in the undeformed and deformed configurations are defined as and , respectively.
For an incompressible, initially isotropic tube, cylindrically axisymmetric deformations are described by
(6) |
where . Thus, the deformation gradient tensor can be represented as
(7) |
where ‘diag’ denotes the diagonal matrix, is the circumferential principal stretch, and is the radial principal stretch. Thus, the left and right Cauchy-Green tensors can be written as in their respective eigenvector bases.
For cylindrically axisymmetric deformations and an applied radial electric field, the biasing Eulerian electric displacement vector only has a radial component and the non-zero component of its Lagrangian counterpart, , is . As a result, the five independent scalar invariants in Eq. (3) can be written now in the form
(8) |
which, when substituted into the initial constitutive relations (4), yields the non-zero components of the total stress tensor and the Eulerian electric field vector as
(9) |
From Eq. (8), we see that the five invariants can be written in terms of three independent variables only, for instance: , and . Then we define a new reduced energy density function and obtain from Eqs. (8) and (9) the following relations:
(10) |
where , , and .
For axisymmetric deformations invariant along the axis, all the initial physical quantities depend only on . In this case, Faraday’s law (1)3 is then fulfilled automatically. With the help of Eq. (10)1, the equation of motion (1)1 and Gauss’s law (1)2 simplify to
(11) |
respectively.
The inner and outer electrode-coated surfaces of the deformed tube carry equal free charges with opposite sign (i.e. ). The electric field vanishes in the surrounding vacuum based on Gauss’s theorem and neglecting edge effects. Then integrating Eq. (11)2, we can obtain the solution of the radial electric displacement as
(12) |
The electric field vector is curl-free, and hence an electrostatic potential is introduced to write . Inserting Eq. (12) into Eq. (10)3 and integrating the resulting equation from to results in the electric potential difference between the inner and outer surfaces as
(13) |
Integrating Eq. (11)1 from to and using the change of variable (which is derived from Eq. (6) with the definition ), we obtain
(14) |
where and are circumferential stretches of the inner and outer surfaces of the SEA tube, respectively. Here, we assume that both the inner and outer surfaces of the tube are traction-free, i.e., . Thus, Eq. (14)2 can be written as
(15) |
which establishes a general expression of nonlinear axisymmetric response of the circumferential stretch to the electrical variable or (which is included in ) for an arbitrary energy function. Similarly, we can obtain from Eq. (14)1 the radial normal stress as
(16) |
The formulations obtained above are completely universal, valid for any isotropic SEA tube. We now specialize the preceding results to the Gent ideal dielectric model, which is characterized by the following (reduced) total energy density function:
(17) |
where is the shear modulus of the SEA tube in the absence of biasing fields, is the permittivity of an ideal dielectric material, and the dimensionless Gent constant reflects the limiting chain extensibility of rubber networks (Gent, 1996), accounting for the strain-stiffening effect. Moreover, the Gent ideal dielectric model (17) reduces to the neo-Hookean model in the limit of .
Substituting Eq. (17)2 into Eq. (13) and integrating, we obtain the explicit expression between the dimensionless electric voltage and dimensionless surface free charge as
(18) |
Substitution of Eq. (17)2 into Eq. (15) gives
(19) |
Making use of the relation and Eq. (12), we can obtain the second integration part in Eq. (19) as
(20) |
The first integration part in Eq. (19) can be derived as follows:
(21) |
where are the two roots of the following quadratic polynomial equation of :
(22) |
Therefore, integrating Eq. (21) yields
(23) |
where
(24) |
Substituting Eqs. (20), (21) and (23) into Eq. (19), we obtain the nonlinear axisymmetric deformation relation between and as
(25) |
By means of Eqs. (18) and (25), we finally get the nonlinear axisymmetric response between and for the Gent ideal dielectric model as
(26) |
Moreover, the radially inhomogeneous biasing fields (including the circumferential stretch , the radial electric displacement , the radial normal stress , and the Lagrange multiplier ) needed to calculate the natural frequency of electrostatically tunable non-axisymmetric vibrations can be derived, using Eqs. (6)1, (12), (16), and (9)1, as
(27) |
where is the dimensionless radial coordinate in the deformed configuration, is the dimensionless radial electric displacement, and and are the dimensionless radial normal stress and dimensionless Lagrange multiplier, respectively.
3 Incremental fields and non-axisymmetric vibration analysis
To analyze the non-axisymmetric vibrations of a finitely deformed SEA tube, we employ the linearized incremental theory of electro-elasticity (Dorfmann and Ogden, 2010, 2014), the main parts of which are first summarized for sake of completeness. Then, the incremental governing equations in cylindrical coordinates are recast into the state-space formalism, which we use, combined with the approximate laminate technique, to obtain the frequency equations.
3.1 Linearized incremental theory
We superimpose a time-dependent infinitesimal incremental motion and an infinitesimal incremental electric displacement upon a finitely deformed configuration (with the boundary and outward unit ). Here and henceforth, a superposed dot indicates the increment in the quantity concerned. According to the incremental theory of electro-elasticity (Dorfmann and Ogden, 2010, 2014), we can write the linearized incremental incompressibility condition, incremental governing equations and incremental constitutive relations for incompressible SEA materials in updated Lagrangian form as
(28) |
(29) |
respectively, where
(30) |
Here, is the incremental mechanical displacement vector, is the incremental displacement gradient tensor, is the incremental Lagrange multiplier, and , , and are the push-forward counterparts of the increments of the total nominal stress and of the Lagrangian electric displacement and electric field, respectively. The subscript ‘’ is used to indicate the resultant push-forward variables.
In Eq. (30), , and are fourth-, third- and second-order tensors, respectively, which are referred to as the instantaneous electro-elastic moduli tensors. Their component forms satisfy
(31) |
where , and are the referential electro-elastic moduli tensors related to the total energy function , with components
(32) |
With no external fields in the surrounding vacuum, the updated Lagrangian forms of the incremental mechanical and electric boundary conditions satisfied on take the following simple forms,
(33) |
where is the updated Lagrangian incremental traction vector per unit area of and is the incremental surface charge density on .
3.2 Incremental equations and state-space formalism in cylindrical coordinates
In this subsection, the incremental governing equations (28)-(30) are first specialized to the cylindrical coordinates . The basic incremental governing equations of the deformed SEA tube are the incremental incompressibility constraint,
(34) |
together with the incremental equations of motion and incremental Gauss’s law,
(35) |
(36) |
and the incremental constitutive equations,
(37) |
where , and are the effective material parameters associated with the instantaneous electro-elastic moduli , and (their explicit expressions are provided in Eq. (41) in the paper by Wu et al. (2017)). In the process of deriving Eq. (37), the incremental displacement gradient tensor was specialized to the cylindrical coordinates. In addition, an incremental electric potential , defined by , was introduced so that the incremental Faraday law (29)3 is satisfied automatically.
Note that the non-zero components of the instantaneous electro-elastic moduli tensors , and were derived by Wu et al. (2017) for the nonlinear axisymmetric deformations of incompressible isotropic SEA tubes subject to a combination of radial electric displacement and axial pre-stretch (see their Appendix B for specific expressions). It is worth emphasizing that instantaneous physical properties of the SEA tubes are significantly sensitive to applied electro-mechanical biasing fields, which generate remarkable knock-on influences on the dynamic characteristics of superimposed small-amplitude motions, as demonstrated below.
As discussed in Sec. 2.2, the physical quantities for cylindrically axisymmetric deformations in the SEA tube subject to the radial electric voltage are radially inhomogeneous, leading to the -dependence of the instantaneous electro-elastic moduli. Thus, the resultant incremental displacement equations are, in general, a set of coupled partial differential equations with variable coefficients, which are likely intractable to solve analytically and difficult to solve numerically. Here we employ the state-space method (SSM) (Wu et al., 2017; Zhu et al., 2020), combining the incremental state-space formalism with the approximate laminate technique, to obtain the frequency equations. Following a standard derivation procedure (omitted here for simplicity), we transform the original incremental governing equations (34)-(37) into the following incremental state equation:
(38) |
which is a first-order system of differential equations with respect to , where , is the incremental state vector, with , , , , , , and being the state variables, and is an 8 8 system matrix, which depends on the radial coordinate, the instantaneous electro-elastic moduli, and partial derivatives with respect to , and . The specific expressions for the elements of the system matrix can be found in the Appendix C in the paper of Wu et al. (2017); they are not reproduced here for brevity. We emphasize that the state equation (38) is applicable for any form of energy function.
3.3 Approximate laminate technique
We assume that generalized rigidly supported conditions (Ding and Chen, 2001) are applied to the two ends of the deformed SEA tube, so that the end cross-sections are in smooth contact with rigid plattens at both ends. In addition, the electric inductions in the surrounding vacuum near the tube ends can be neglected, so that the incremental electric displacement at both ends is zero. Therefore, the incremental mechanical and electric boundary conditions are
(39) |
To satisfy the incremental boundary conditions (39), we assume the following formal solutions for the harmonic non-axisymmetric free vibrations of the deformed SEA tube:
(40) |
where and are the dimensionless radial and axial coordinates, respectively; and denote the circumferential mode number and axial mode number, respectively; = is the imaginary unit, and is the circular frequency of free vibrations. Note that all the unknown functions (, etc.) in Eq. (40) are dimensionless. Then, inserting Eq. (40) into the state equation (38) and non-dimensionalizing related variables, we obtain the dimensionless form of the incremental state equation as
(41) |
where is the dimensionless incremental state vector and is the 8 8 dimensionless system matrix, with its four partitioned sub-matrices given by
(42) |
in which the dimensionless quantities are defined as
(43) |
and is the dimensionless circular frequency.
Note that different combinations of mode numbers and result in different types of vibrations. It is obvious from Eqs. (41) and (42) that the eight unknown incremental state variables are fully coupled for non-axisymmetric vibrations with and (see Fig. 1(c)).
For the incremental axisymmetric vibrations with and , the incremental fields are independent of the coordinate (i.e. ) and the state equation (41) reduces to two uncoupled classes of incremental axisymmetric vibrations: the axisymmetric longitudinal vibrations (L vibrations) where the mechanical displacement components and coupled with the incremental electrical quantities are non-zero (see Fig. 1(d)); and the purely torsional vibrations (T vibrations) with the sole displacement component (see Fig. 1(e)). The cylindrical breathing mode with , characterized by the sole radial displacement component , is a special mode of the L vibrations. The state equations governing the incremental axisymmetric vibrations of deformed SEA tubes were obtained by Zhu et al. (2020) (see their Eqs. (30), (31), (41) and (42)) and are omitted here for brevity.
Additionally, we consider incremental fields independent of the coordinate (i.e. and ), so that holds. To satisfy the incremental boundary conditions (39), the harmonic solutions for the incremental vibrations independent of are assumed as
(44) |
and hence, the incremental state equation (41) specializes to
(45) |
where is the dimensionless incremental state vector, and the dimensionless system matrix is obtained from Eq. (42) as
(46) |
Therefore, the incremental vibrations independent of ( and ) are described by Eqs. (44)-(46), which are identified as the prismatic vibrations, where the non-zero mechanical displacement components and are coupled with the incremental electrical quantities and the SEA tube remains prismatic while its cross-section loses its circular shape (see Fig. 1(f)). Furthermore, prismatic vibrations do not depend on the length-to-thickness ratio . Note that the nomenclature of prismatic vibrations corresponds to the elastic counterpart of prismatic diffuse modes in the realm of instability (Haughton and Ogden, 1979; Bortot and Shmuel, 2018). When , Eqs. (44)-(46) governing the prismatic vibrations reduce to those of the cylindric breathing mode.
Now we adopt the approximate laminate technique (Fan and Zhang, 1992; Chen and Ding, 2002) and divide the tube into a laminate with equal thin sublayers. The thickness of each sublayer is , which is sufficiently small for the system matrices and within each sublayer to be approximately regarded as constant. Without loss of generality, the dimensionless deformed radial coordinate itself and the effective material parameters take the values at the mid-plane of each sublayer.
Here, we use , , and to denote the deformed radial coordinates of the inner, outer and middle surfaces of the -th sublayer. Their dimensionless forms of the deformed radial coordinate are
(47) |
respectively. Applied to each sublayer, the incremental state equation (41) becomes
(48) |
where is the approximated constant system matrix within the -th thin sublayer, obtained by setting . Consequently, the formal solution to Eq. (48) in the -th sublayer is
(49) |
which yields the transfer relation between the incremental state vectors at the inner and outer surfaces of the -th sublayer as
(50) |
Making use of the continuity conditions at the fictitious interfaces between equally divided sublayers that require the eight state variables be continuous, we can derive from Eq. (50) the following transfer relation connecting the incremental state vectors and at the outer and inner surfaces of the deformed SEA tube:
(51) |
where is the global transfer matrix. Similar derivations are applicable to the case of the prismatic vibrations described by and .
3.4 Frequency equations of non-axisymmetric vibrations
To proceed, we assume that the inner and outer surfaces of the SEA tube are traction-free and that the applied electric voltage between these two surfaces is kept fixed. Thus, the incremental mechanical and electric boundary conditions (33) are:
(52) |
Combining Eq. (52) with Eq. (51) results in a set of independent linear algebraic equations:
(53) |
where denote the elements of the global transfer matrix . For non-trivial solutions to exist, the determinant of the coefficient matrix in Eq. (53) must vanish. Thus, we have
(54) |
which is the characteristic frequency equation for small-amplitude non-axisymmetric vibrations of the activated SEA tube subject to radially inhomogeneous biasing fields for mode numbers and .
Following the same derivation procedure, we can also acquire the characteristic frequency equation of the incremental prismatic vibrations () as:
(55) |
where are the elements of the global transfer matrix for the prismatic vibrations.
For the incremental axisymmetric vibrations including the L vibrations, the T vibrations and the breathing vibrations, we refer to the frequency equations (40) and (43) in the paper by Zhu et al. (2020).
Note that the frequency equations derived above are applicable for any form of the energy function of the incompressible isotropic SEA tubes. To calculate numerically the natural frequency of electrostatically tunable non-axisymmetric vibrations, we specialize the analysis to the Gent ideal dielectric model (17). Then the dimensionless effective material parameters appearing in Eqs. (42) and (43) are
(56) |
where = and = . In the limit of , we have and , and the effective material parameters listed in Eq. (56) reduce to their equivalent forms for the neo-Hookean ideal dielectric model (see Eq. (44) in Zhu et al. (2020)).
4 Numerical results and discussions
In this section, we investigate the influences of the electro-mechanical biasing fields (i.e., the combined action of radial voltage and axial pre-stretch ) and strain-stiffening effect on the nonlinear response and non-axisymmetric vibration characteristics of the SEA tube characterized by the Gent ideal dielectric model (17).
4.1 Axisymmetric nonlinear static response
We first investigate the axisymmetric nonlinear static response under the combined action of electric voltage and axial pre-stretch.
Based on the nonlinear response equations described in Section. 2.2 for the neo-Hookean and Gent ideal dielectric models, we plot in Fig. 2 the nonlinear static response variation curves (i.e. the dimensionless electric voltage versus the circumferential stretch ) of the SEA tube under different electro-mechanical biasing fields. In Fig. 2, the Gent constant models the strain-stiffening effect of the SEA tube near a limiting stretch and the ‘N-H model’ legend denotes the neo-Hookean ideal dielectric model. It should be emphasized that the electro-mechanical limit-point instability (i.e., the reaches a maximum) occurs in neo-Hookean SEA tubes (Shmuel and DeBotton, 2013; Shmuel, 2015; Wu et al., 2017; Zhu et al., 2020) but not in Gent SEA tubes because of the strain-stiffening effect (Zhou et al., 2014; Bortot and Shmuel, 2018). Fig. 2 displays the nonlinear response curves for different values of with a fixed pre-stretch , and the variation curve for the neo-Hookean model is also included. It shows that increases monotonically as the radial electric voltage increases, which means physically that the tube reduces in thickness and expands in the radial direction with increasing . When the applied voltage keeps growing, the tube reaches its limiting stretch due to the effect of strain-stiffening for a Gent SEA tube. The smaller is, the smaller the limiting stretch is, and the stronger the strain-stiffening effect is, while for large values (e.g. ), the curve tends to that of the neo-Hookean model.
In Fig. 2, we additionally take the influence of axial pre-stretch into account. We mark out the electro-mechanical instability voltages (i.e., for , for , and for ) for the neo-Hookean model under three different pre-stretches for reference. For a fixed value, when subject to larger axial pre-stretch, lower radial voltage is required to obtain equal circumferential stretch , which means physically that the SEA tube is easier to deform for a larger axial pre-stretch.
4.2 Validation of the State Space Method (SSM)
The effectiveness of the SSM is first verified in terms of its convergence and accuracy for the 3D free vibration analysis of the Gent SEA tube subject to biasing fields in this section.
NoL | 20 | 40 | 60 | 80 | 100 | 120 | 140 | 160 |
---|---|---|---|---|---|---|---|---|
0.30495 | 0.30495 | 0.30495 | 0.30495 | 0.30495 | 0.30495 | 0.30495 | 0.30495 | |
0.35255 | 0.35255 | 0.35255 | 0.35255 | 0.35255 | 0.35255 | 0.35255 | 0.35255 | |
0.44378 | 0.44378 | 0.44378 | 0.44378 | 0.44378 | 0.44378 | 0.44378 | 0.44378 |
NoL | 20 | 40 | 60 | 80 | 100 | 120 | 140 | 160 |
---|---|---|---|---|---|---|---|---|
1.20181 | 1.20180 | 1.20180 | 1.20180 | 1.20180 | 1.20180 | 1.20180 | 1.20180 | |
1.97680 | 1.97698 | 1.97702 | 1.97703 | 1.97704 | 1.97704 | 1.97704 | 1.97704 | |
2.85347 | 2.85284 | 2.85273 | 2.85269 | 2.85267 | 2.85266 | 2.85265 | 2.85265 |
For the convergence analysis, we list in Tables 1 and 2 the first three natural frequencies of the non-axisymmetric vibration calculated by the SSM with different numbers of discretized layers (NoL) for a thin and slender SEA tube and a thick and short SEA tube of the Gent ideal dielectric model, respectively. Moreover, the circumferential mode number and axial mode number are all equal to one. Clearly, we can find the natural frequencies with arbitrary precision using the current SSM because the results clearly indicate an excellent convergence rate with increasing NoL. The calculated frequencies based on the SSM are close to those of the original SEA tube that is subject to radially inhomogeneous biasing fields. Thus, 120 discretized layers will be chosen hereafter, which is assumed to have high accuracy.
The radially inhomogeneous biasing fields are induced by the applied radial electric voltage. However, when there is no applied radial electric voltage, the deformation of the pre-stretched SEA tube is homogeneous, which makes it feasible to obtain the exact frequency solutions to the non-axisymmetric vibrations through the conventional displacement method. The detailed derivations for the analytical frequency equations of the non-axisymmetric vibrations were provided in Appendix C of our previous work (Zhu et al., 2020) for an arbitrary energy function. A in this paper additionally gives the effective material parameters and analytical frequency equations of three kinds of special vibrations (including the breathing mode, T vibrations and prismatic vibrations) in a pre-stretched hyperelastic tube characterized by the Gent model. Therefore, the accuracy of SSM can be verified by making a comparison to the results obtained from the exact solutions.
Based on the SSM and the exact solutions, we carry out an accuracy analysis of the first two dimensionless vibration frequencies versus the axial mode number or the circumferential mode number for the pre-stretched () thin and slender SEA tube and for the thick and short SEA tube with no applied voltage, for the following four kinds of vibrations: Fig. 3 for L vibrations including the breathing mode , Fig. 3 for T vibrations , Fig. 3 for non-axisymmetric vibrations , and Fig. 3 for prismatic vibrations . For the Gent material model, the dimensionless Gent constant is chosen to be 97.2.
In Fig. 3, regardless of geometric sizes, there is only one radial vibration frequency for the breathing mode as a result of the tube incompressibility. Moreover, the natural frequency of breathing mode in a thin and slender SEA tube is lower than those of other first-order L vibration modes , while in the case of thick and short tube, the natural frequency of the breathing mode is larger than those of the three L vibration modes of the first order. In addition, we make the link with the data in the case of a thick and short neo-Hookean SEA tube from Zhu et al. (2020) and that of the Gent model with a large () in both Figs. 3 and 3, denoted by the black-dashed line and green diamond, respectively. As expected, the natural frequencies of the Gent tube with a large coincide with those of the neo-Hookean tube. For the four kinds of vibrations shown in Fig. 3, it is clear that the vibration frequencies calculated by the SSM agree with those predicted from the exact solutions in the entire axial and circumferential mode number range.
In conclusion, the SSM-based numerical results are highly accurate for the 3D free vibration analysis because of the great convergence rate and the excellent agreement with the exact solutions.
4.3 Strain-stiffening effect on axisymmetric and prismatic vibrations
In this subsection, we investigate how the strain-stiffening effect affects the natural frequencies of axisymmetric vibrations (including the L vibrations and T vibrations) and prismatic vibrations. If not otherwise stated, we consider a thin and slender tube with the geometric sizes set as and .
First, for a pre-stretched () thin and slender Gent ideal SEA tube with , the variations of the first-order dimensionless natural frequency with the axial mode number are displayed in Figs. 4 and 4 under different radial electric voltages for the L and T vibrations, respectively. In Fig. 4, we exhibit the variation trend of the first-order natural frequency with the circumferential mode number for the prismatic vibrations.
It is seen from Fig. 4 that the first-order frequency of L vibrations goes up monotonously with increasing . As the voltage increases, the natural frequency decreases in the entire axial mode number range. Specifically, we find that a relatively low voltage (e.g., ) barely affects the L vibration frequency, but as continues growing, there exists a big gap between the two frequency spectra stimulated by and . Such a dramatic change occurs due to the tube’s rapid expansion, resulting in the sharp decrease of the global stiffness when .
For the T vibrations shown in Fig. 4, the vibration frequency increases monotonously and linearly with the axial mode number for an arbitrary applied radial voltage . Taking the electro-mechanical instability voltage exclusive to the neo-Hookean SEA tube as the demarcation point, the change of voltage makes no difference to the vibration frequency for this kind of thin and slender pre-stretched SEA tube for . However, for , the natural frequency has a remarkable rise with the applied voltage because of the strain-stiffening effect.
For the prismatic vibrations depicted in Fig. 4, it is found that the circumferential mode number corresponding to the lowest natural frequency depends on the applied voltage when . For example, the lowest vibration frequency is obtained at for , while for , it occurs at . For (such as the cases and ), the natural frequency increases monotonically with the increasing circumferential mode number and the lowest frequency is achieved at . Moreover, the natural frequency for goes up monotonically with the applied voltage in the entire circumferential mode number range.
Then, to clearly reveal the influences of voltage and strain-stiffening effect on the vibration behaviors, we exhibit in Fig. 5 the variations of the natural frequency with the dimensionless radial voltage for an SEA tube characterized by the Gent and neo-Hookean models at three different axial pre-stretches. Four specific cases are included in Fig. 5: breathing mode in Fig. 5; L vibration mode in Fig. 5; T vibration mode in Fig. 5 and prismatic vibration mode in Fig. 5.
For the breathing mode (which can be seen as a special case of L vibrations) shown in Fig. 5, the lines correspond to the Gent model while the symbols denote the neo-Hookean model, and different line and symbol styles represent different axial pre-stretches. Clearly, the natural frequencies for the neo-Hookean SEA tube decrease nonlinearly to zero with the increase in voltage. The voltage corresponding to is called the critical voltage and this threshold is referred to as the breathing mode instability. Axisymmetric barreling instabilities (Simpson and Spector, 1984; Chen et al., 2017; Wu et al., 2020) in the SEA tube arise from the global stiffness rapidly decreasing when the applied voltage gets closer to the critical value . The critical voltages of the breathing mode for the neo-Hookean SEA tube are , , at the axial pre-stretches , , , respectively, which are identical to the electro-mechanical instability voltages mentioned in Section 4.1. Meanwhile, we find that a larger axial pre-stretch results in a lower critical voltage , which means the axial pre-stretch destabilizes the neo-Hookean SEA tube. Moreover, the curves for the neo-Hookean model all start from an identical point when , because the natural frequency of the breathing mode is independent of axial pre-stretch in the absence of voltage (Zhu et al., 2020). For the Gent SEA tube, the natural frequency of the breathing mode goes down at first due to the decrease of the global stiffness and then increases conversely and rapidly because of the strain-stiffening effect when the voltage gradually grows, and the breathing mode instability is eliminated. In addition, the analytical solution Eq. (62) to the breathing mode of the Gent SEA tube when demonstrates that although its natural frequency varies with the axial pre-stretch , the chosen value of Gent constant is large enough to eliminate the effect of the pre-stretch considered; that is why the curves for the Gent model almost start from an identical point when , which is similar to the result of the neo-Hookean model.
For the L vibration mode depicted in Fig. 5, the variation curves of the first-order vibration frequency with the applied voltage are presented for the Gent and neo-Hookean models () at three cases of axial pre-stretch. Generally, the natural frequencies of this L vibration mode for both neo-Hookean and Gent models all decline nonlinearly to zero with an increase in . The phenomenon that a larger axial pre-stretch results in a lower critical voltage exists for both the neo-Hookean model and Gent model with , which is qualitatively similar to that of the breathing mode for the neo-Hookean model. However, for the Gent model with , the critical voltage of an SEA tube subject to either axial pre-extension or pre-compression is larger than that without axial pre-stretch because of the earlier strain-stiffening effect.
For the T vibration mode , we plot the first-order natural frequency of a thin and slender or thick and short SEA tube for three different axial pre-stretches in Fig. 5. The Gent model and neo-Hookean model are all presented for comparison. The neo-Hookean numerical data for a thick and short SEA tube is taken from Zhu et al. (2020). Lines represent the Gent model while symbols denote the neo-Hookean model. For the neo-Hookean model, it is obvious that the horizontal first-order natural frequencies for the SEA tubes in two geometries are independent of the applied voltage and axial pre-stretch within the range of , as explained in Appendix D of Zhu et al. (2020). As for the Gent model, the natural frequency curves of the thin and slender tube are also horizontal when , but as exceeds , its natural frequency curves rise with the increase of voltage. However, when it comes to the thick and short Gent tube, it is apparent that the frequency curves change nonlinearly and finally reduce to zero when the applied voltage gradually increases from zero. Thus, we can find that the increase of the Gent tube thickness strengthens its sensitivity of natural frequency to the voltage. Moreover, although the analytical solution Eq. (65) to the first-order T vibration mode of the Gent SEA tube shows that its natural frequency at relies on the axial pre-stretch , the Gent constant value is large enough to weaken the effect from , making the natural frequencies at different axial pre-stretches almost start from an identical point when .
In Fig. 5, we display the variation curves of prismatic vibration mode in a thin and slender SEA tube for the first-order vibration frequency in response to the applied voltage. A similar variation trend to the breathing mode in Fig. 5 is observed for both the neo-Hookean and Gent models. For example, the natural frequency for the neo-Hookean SEA tube reduces nonlinearly to zero when the voltage increases. The point where for the neo-Hookean model in Fig. 5 corresponds to the prismatic diffuse instability (Haughton and Ogden, 1979; Bortot and Shmuel, 2018). We observe from Fig. 7(a) that the shape of tube cross-section of the prismatic vibration mode for is analogous to that of the breathing mode; at the same time, we are considering the case of a thin-walled tube, which explains why the critical voltages of prismatic vibration mode for the neo-Hookean model in Fig. 5 are very close to those of breathing mode for the neo-Hookean model in Fig. 5. For the Gent SEA tube, when the applied voltage increases, the natural frequency first goes down to a small value because of the decrease of the global stiffness and then increases due to the strain-stiffening effect. Therefore, for this specific prismatic vibration mode , the Gent SEA tube with strain-stiffening effect eliminates the prismatic diffuse instability that can occur in the neo-Hookean tube. However, for other prismatic vibration modes with , the prismatic diffuse instabilities might occur in the Gent SEA tube (Bortot and Shmuel, 2018), as described below.
We further plot the curves of prismatic vibration frequency as functions of the radial voltage for in Fig. 6 for both the neo-Hookean and Gent SEA tubes. We observe from Fig. 6 for the neo-Hookean model that all of the prismatic vibration frequencies decrease nonlinearly to zero in response to the continuously increasing radial voltage. The critical voltage for is higher than those for and the critical voltage goes up with for . Specifically, the critical voltages of prismatic vibrations for different circumferential mode numbers are: for , for , for , for , and for . However, we notice from Fig. 6 that prismatic instability of for the neo-Hookean model does not exist for the Gent model. Similar to Fig. 6, as the radial voltage gradually increases, the first-order vibration frequency decreases to zero with for , for , for , and for , respectively. The critical voltage for the Gent model increases monotonically when increasing the circumferential mode number for . Moreover, the critical voltage of the same circumferential mode number for the Gent model is slightly larger than that for the neo-Hookean model. These phenomena agree well with the results obtained by Bortot and Shmuel (2018) and their Figs. 4(a-c) illustrated the critical voltage for the prismatic diffuse instability as a function of thickness-to-mean radius ratio of the SEA tubes clamped without axial pre-stretch or clamped after pre-stretch, where cases for both the Gent and neo-Hookean models were presented for various circumferential mode numbers.
For a better illustration, we present in Fig. 7 four prismatic vibration mode shapes for in a pre-stretched () thin and slender Gent SEA tube subject to a radial voltage . The intensity of circumferential gridlines represents the distribution of radial displacement, while the radially scattered gridlines denote the distribution of circumferential displacement. According to Eq. (44), the prismatic vibrations only have two coupled components: the radial displacement component and the circumferential displacement component . Thus, the SEA tube remains prismatic in axial direction while its cross-section loses its circular shape for . The circumferential mode number plays an important role in shaping the prismatic-vibrating mode shape. It can be seen from Fig. 7 that vibration crests and troughs are distributed in the prismatic mode shapes for the th order circumferential mode number.
To better understand how strain stiffening affects axisymmetric and prismatic vibrations, we present variation curves of the first-order natural frequency with the Gent parameter in Fig. 8 for a thin and slender SEA tube under different axial pre-stretches and voltages. We examine the strain-stiffening effect under two circumstances: no radial voltage and high radial voltage close to the of the neo-Hookean model. We see in Fig. 8 that as the Gent parameter increases to infinity, the frequency tends to converge to that of the neo-Hookean model, as expected. Recalling Figs. 5(a), 5(b) and 5(d), the natural frequencies approach zero as the applied radial voltage approaches in a neo-Hookean SEA tube. Therefore, for the high radial voltage close to the , the vibration frequency decreases monotonically to zero with the increase of , as shown in the bottom subplots of Figs. 8(a), 8(b) and 8(d). However, the natural frequency for the T vibration mode in Fig. 8(c) decreases to a finite value when increasing to infinity, because the first-order T vibration frequency in a thin and slender SEA tube is independent of radial voltage and axial pre-stretch within the range of , regardless of whether the neo-Hookean or Gent model is considered, as displayed in Fig. 5(c). When there is no radial voltage in each top subplot of Figs. 8(a)-(d), the SEA tube subject to pre-extension or pre-compression gets stiffened rapidly with the strain as becomes relatively small, which explains that the natural frequency increases drastically when decreases further from the small value. In other words, the stronger the strain-stiffening effect is, the easier it is to achieve a strain-stiffened state under a small deformation. Additionally, it can be seen from the top subplots of Figs. 8(a)-(d) that for a fixed pre-stretch and zero radial voltage, the natural frequencies for all vibration modes are independent of the Gent parameter , which physically means that the strain-stiffening effect has no influence on the vibration frequency of a SEA tube without electro-mechanical biasing fields.
4.4 Effect of electro-mechanical biasing fields on non-axisymmetric vibrations.
In this subsection, we expound on how the electro-mechanical biasing fields influence the non-axisymmetric vibration characteristics in a thin and slender SEA tube.
First, we display in Fig. 9 the variation of the first-order dimensionless frequency with the axial mode number in a pre-stretched () Gent SEA tube for different values of radial voltage. Two representative circumferential mode numbers and are selected in Figs. 9 and 9, respectively. It can be seen that the vibration frequency increases almost linearly with the increase of axial mode number . For , a low radial voltage (e.g., ) barely affects the natural frequency, whereas the natural frequency decreases rapidly when the SEA tube is subject to a relatively higher radial voltage (e.g., ). Moreover, when comparing Fig. 9 with to Fig. 9, we observe a more noticeable frequency change within the range of small (e.g., ) when increasing the radial voltage with .
Then, to investigate the influence of circumferential mode number , we plot in Fig. 10 the dimensionless frequency spectra ( versus ) in a pre-strecthed () Gent SEA tube for different values of axial mode number and radial voltage. Specifically, in Fig. 10, we illustrate the first two vibration frequencies without radial voltage for three axial mode numbers . We observe that the first-order frequency change is very small when , and that the frequency then goes up gradually with the increase of . In addition, the second-order frequency increases significantly and monotonically with the axial mode number . In Fig. 10, we consider two axial mode numbers and four radial voltages . We see that the overall frequency curves of non-axisymmetric vibration decline when increasing the radial voltage across the entire range of circumferential mode number. In particular, the voltage-induced frequency reduction of a large circumferential mode number is more obvious than that of a small one.
In Fig. 11, we depict the variation curves of the natural frequency of non-axisymmetric vibration mode with the radial voltage for four representative axial pre-stretches . Both the Gent model and neo-Hookean model are shown for comparison. Lines represent the results based on the Gent model, while symbols indicate those predicted by the neo-Hookean model.
Fig. 11 shows the variation curves of versus with the SEA tube subject to a fixed pre-extention . The first-order frequency goes down to zero monotonically when increasing the voltage for both the neo-Hookean and Gent models. The point where for the non-axisymmetric vibrations corresponds to the 3D non-axisymmetric buckling instability (Haughton and Ogden, 1979; Su, 2020). In the context of voltage-controlled non-axisymmetric instabilities, Su (2020) illustrated the critical stretch and pattern shapes associated with the buckling of SEA tubes in his Fig. 10 and considered specific radius ratios and length aspect ratios, under various applied voltages and torsions. We note from Fig. 11 that the veering phenomenon occurs between the second-order mode (Mode 2) and third-order mode (Mode 3) for the Gent model, while mode crossing happens for the neo-Hookean model. The term ‘veering’ here refers to the situation where two branches approach each other and then veer away and diverge instead of crossing. Mathematically, veering is accompanied by rapid changes in the eigenvectors (Mace and Manconi, 2012; Wu et al., 2017). The sequence of mode order needs to be clarified here: for the mode veering phenomenon, since there is no mode crossing among branches, the sequence is still defined based on the frequency, from low to high. However, for the mode-crossing phenomenon, the modes that remain the same before and after the crossing point are defined as modes with the same order, rather than being based on their frequency.
To proceed, we provide the mode veering validation for the Gent model with top view of the vibration modes in the veering progress in Fig. 12. Four representative radial voltages () are chosen. The second- and third-order vibration modes at are illustrated in Figs. 12(a) and 12(e) for reference as the initial vibration modes. A circular end surface and evenly distributed thickness in Fig. 12(a) demonstrate that the circumferential variations of radial displacement are hardly noticeable, while large circumferential changes in are illustrated in Fig. 12(e), where the tube thickness is not distributed uniformly. As the voltage grows to 0.4 (before the mode veering point), the second- and third-order vibration mode shapes keep consistent with the initial vibration modes with . However, when the voltage is further increased to 0.45 or 0.5, the mode shapes between Mode 2 and Mode 3 switch. In Fig. 12(d), the circumferential change in is exhibited, which is consistent with the mode shape in Fig. 12(e). In Fig. 12(h), the tube end surface remains circular and its thickness is evenly distributed, which corresponds to the mode shape shown in Fig. 12(a). Therefore, when the applied voltage is close to triggering the veering, a slight change in voltage can significantly alter the vibration mode shape. Usually, it leads to a flow of mechanical energy between adjacent vibration modes, which physically expresses a modal transformation. The validation of mode crossing predicted by the neo-Hookean model follows the results in Figs. 13 and 14 for and , respectively.
In Fig. 11, there is no axial pre-stretch, . As the radial voltage increases, the first- and third-order natural frequencies decrease monotonically toward zero (i.e., the 3D non-axisymmetric buckling instability occurs (Haughton and Ogden, 1979; Su, 2020)), while the second-order natural frequency remains independent of the voltage within the range of for both the neo-Hookean and Gent models. Moreover, a mode-crossing phenomenon appears between Mode 2 and Mode 3 for both energy models, and its validation for the Gent model is illustrated in Fig. 13, where we provide the vibration mode shapes of Mode 2 and Mode 3 stimulated by two selected radial voltages and before and after the mode-crossing point. When examining the second-order mode shape at in Fig. 13(a), the top view reveals that the circumferential changes in are barely noticeable and the cross-section remains circular, while the sharp spiral gridlines observed from the front view indicate that undergoes significant changes along the axial direction. Then, considering the third-order mode shape at in Fig. 13(b), the top view reveals significant circumferential change in leading to uneven thickness distribution in the cross-section, while from the front view, only slight change in along the axial direction can be observed. For (after the mode-crossing point), the vibration mode shapes of Mode 2 in Fig. 13(c) and Mode 3 in Fig. 13(d) remain the same as those at . Therefore, we refer to the branch crossing between Mode 2 and Mode 3 in Fig. 11 as the mode crossing, where it is important to emphasize again that the definition of the same order modes is based on the similarity of modes before and after the crossing point, rather than the frequency values of modes. Numerical calculations (not reproduced here) show that a similar phenomenon of mode crossing also exists in the neo-Hookean model.
In Figs. 11 and 11, the axial pre-stretch is compressive, as and , respectively. Clearly, the vibration frequencies of the second-order mode in Fig. 11 and of the first-order mode in Fig. 11 are all independent of the radial voltage within the range of . However, as the radial voltage increases, the frequencies of other vibration modes decrease to zero monotonically, where the 3D non-axisymmetric buckling instabilities occur. In particular, the critical voltage of the first vibration mode initially increases and then decreases until it reaches zero as the axial pre-stretch changes from axial pre-extension to axial pre-compression. Similar to Fig. 11, the mode-crossing phenomenon happens between Mode 2 and Mode 3 in Fig. 11 and between Mode 1 and Mode 2 in Fig. 11 for both the neo-Hookean and Gent models. This mode-crossing phenomenon can be validated by comparing the vibration modes before and after the crossing point. As shown in Fig. 14, the mode shape variations of Mode 1 and Mode 2 for the Gent model and , before () and after () the crossing point, are essentially the same as those depicted in Fig. 13.
To demonstrate the diversity of modes, we present in Fig. 15 the 3D mode shapes of non-axisymmetric vibrations with different combinations of circumferential and axial mode numbers for and . The illustrated 3D mode shapes are composed of three coupled displacement components: the radial displacement , circumferential displacement , and axial displacement . The SEA tube vibrates sinusoidally in the axial and circumferential directions following Eq. (40). Moreover, the axial and circumferential mode numbers are integer multiples of the half-wavelength in the axial direction and the wavelength in the circumferential direction, respectively.
To clearly reveal the effect of axial pre-stretch on the non-axisymmetric vibration characteristics, the variation curves of the first three natural frequencies with the axial pre-stretch in a Gent SEA tube are plotted in Fig. 16 for and under four different radial voltages , , , and . Taking Fig. 11 into consideration, it is clear from Fig. 16 that the second-order natural frequency is barely affected by the axial pre-stretch when . But when (e.g., ), a noticeable frequency veering phenomenon occurs near between Mode 2 and Mode 3, which explains that the second-order frequency for decreases with the axial pre-stretch for , as shown in Fig. 16. Additionally, we observe that the first-order frequency for also decreases when the axial pre-stretch increases from 1.7 to 2. This is because the geometric sizes of the tube increase and the overall stiffness declines as a result of the high voltage close to . Nevertheless, the first- and third-order natural frequency curves diverge from the second-order one in opposite directions when decreasing the axial pre-stretch for . To be specific, the third-order natural frequency continues to increase as decreases, while the first-order frequency goes down gradually as decreases from axial pre-extension to axial pre-compression. In particular, the first-order frequency curve disappears (i.e., the first-order frequency drops to and the 3D non-axisymmetric buckling instability happens) when the axial pre-stretch reaches a critical value . For example, the critical axial pre-stretches are , , , and for , , , and , respectively. Thus, the critical axial pre-stretch has a rise when increasing the radial voltage, which destabilizes the SEA tube. This disappearance can also be observed in Fig. 11, where the first-order frequency curve gradually disappears when decreases from axial pre-extension to axial pre-compression.
Moreover, in order to study the strain-stiffening effect on the non-axisymmetric vibration behaviors, particularly when the SEA tube is subject to low and high radial voltages, the first three natural frequencies of the non-axisymmetric vibration mode are depicted in Fig. 17 as functions of the Gent parameter for a pre-stretched SEA tube. It can be seen that for both low and high radial voltages, the first three natural frequencies decrease rapidly at first with the increase of , and then gradually tend to the straight lines representing the results predicted by the neo-Hookean model. Furthermore, comparing Fig. 17(a) with Fig. 17(b), the corresponding Gent parameter , for which the frequency predictions based on the Gent model start to deviate from those by the neo-Hookean model, increases from approximately for to for . In other words, the strain-stiffening effect manifests itself earlier at relatively high voltages. We note that, as the Gent parameter further decreases from the deviation point between the Gent and neo-Hookean models, the strain-stiffening effect becomes increasingly pronounced, resulting in a rapid increase in frequency. Similar phenomena can be observed for axial pre-compression and are not shown here for brevity.
When there is no axial pre-stretch (i.e., ), it can be found from Fig. 11 that the 3D non-axisymmetric buckling instability voltage () of the first mode for the Gent model is far less than the electro-mechanical instability voltage . According to the results when in Fig. 2, the Gent SEA tube does not get strain-stiffened for , even when the Gent parameter decreases to 10. Thus, the first-order natural frequency of the SEA tube for the mode does not change with the Gent parameter unless it is chosen to be extremely small (e.g., ).
In order to more clearly demonstrate the voltage-controlled vibration frequency of a specific material and discuss the failure possibility of electric breakdown (EB) under the combination of extreme deformation and high voltage conditions for the Gent model, we have selected Silicone CF19-2186 (Shmuel and Pernas-Salomón, 2016) by the manufacturer Nusil as a specific numerical example, with its material properties being , , , and , where and are the relative permittivity and dielectric strength, beyond which the EB phenomenon occurs. Using the material properties of Silicone CF19-2186, the dimensionless breakdown electric field is calculated as , where is the material permittivity with . For the static axisymmetric deformation of SEA tubes subject to axial pre-stretch and radial voltage, the radial electric field can be derived as based on the ideal dielectric energy model (17)2 and constitutive equation (10)3. Thus, based on Eq. (18) the radial electric field of the SEA tube is expressed as , which is inversely proportional to the radial coordinate . Clearly, the maximum value of the radial electric field is obtained at inner surface of the SEA tube and its dimensionless form is .
For a pre-stretched thin and slender Gent SEA tube made of Silicone CF19-2186, we plotted in Fig. 18 the lowest dimensionless vibration frequencies of five different vibration modes as functions of dimensionless radial voltage. It is seen that similar to the previous results, the vibration frequency of the breathing mode first decreases and then increases because of the strain-stiffening effect, while that of the T vibration mode remains unchanged at first and then rises quickly. The frequencies of the L vibration, prismatic vibration and non-axisymmetric vibration modes decrease monotonically to zero, corresponding to the relevant instabilities. The dimensionless instability voltages are marked in Fig. 18, which are for the prismatic vibration mode , for the L vibration mode , and for the non-axisymmetric vibration mode . The EB failure possibility can be determined by comparing with for different applied voltages, as shown in Fig. 18. It is apparent from Fig. 18 that the dimensionless EB voltage is larger than the instability voltages of the relevant vibration modes, which means that the buckling instability or prismatic diffuse instability of a SEA tube made of Silicone CF19-2186 occur before its electric breakdown. In addition, Fig. 18 demonstrates that electrostatically tunable vibration characteristics of SEA tubes made of Silicone CF19-2186 is feasible over a wide voltage range.
5 Conclusions
We investigated the electrostatically tunable non-axisymmetric vibration characteristics of an SEA cylindrical tube with strain-stiffening effect under inhomogeneous biasing fields induced by an axial pre-stretch and a radial electric voltage. First, we used the finite electro-elasticity theory to derive the governing equations of nonlinear axisymmetric static response and radially inhomogeneous biasing fields of an incompressible SEA tube characterized by the Gent ideal dielectric model. Next, based on the relevant linearized theory for incremental fields proposed by Dorfmann and Ogden, we employed the State-Space Method (SSM) to tackle the inhomogeneous biasing fields and obtain the frequency equations for small-amplitude prismatic vibrations and non-axisymmetric vibrations of the activated SEA tube. Finally, we conducted numerical calculations to verify the accuracy and convergence of the SSM in dealing with the non-axisymmetric vibrations, and to thoroughly study the influence of the strain-stiffening effect on the axisymmetric and prismatic vibrations, as well as the influences of electro-mechanical biasing fields and strain-stiffening effect on the non-axisymmetric vibrations. We can summarize our main results as follows:
-
(1)
The SSM is a highly efficient and accurate method to study the superimposed non-axisymmetric vibrations of SEA tubes subject to inhomogeneous biasing fields.
-
(2)
As the radial voltage increases, the vibration frequency of the breathing mode predicted by the neo-Hookean model monotonically decreases to zero due to the reduction in global stiffness, while that anticipated by the Gent model first reduces to a small value owing to the attenuation of global stiffness, and then increases conversely because of the strain-stiffening effect.
-
(3)
For the L vibrations, the vibration frequencies predicted by both the neo-Hookean and Gent models exhibit a nonlinear decline towards zero with the increase of radial voltage and the axisymmetric barreling instabilities occur.
-
(4)
For the T vibrations of thin SEA tubes characterized by the Gent model, the lowest vibration frequency does not change with the radial voltage when it is less than the electro-mechanical instability voltage . However, as the voltage exceeds , the lowest vibration frequency starts to increase significantly as a result of the strain-stiffening effect.
-
(5)
For the prismatic vibration with a circumferential mode number , its lowest vibration frequencies for both the neo-Hookean and Gent models decrease nonlinearly to zero when continuously increasing the voltage and prismatic diffuse instabilities appear in the SEA tube. The critical instability voltage increases monotonically with the circumferential mode number.
-
(6)
For the non-axisymmetric vibrations, the lowest vibration frequencies monotonically go down to zero when increasing the voltage for both the neo-Hookean and Gent models and 3D non-axisymmetric buckling instabilities happen. The lowest frequency curve gradually disappears when the axial pre-stretch goes from pre-extension to pre-compression. Mode veering or crossing phenomena exist in two higher-order vibration modes resulting in modal conversion, which can be judged by analyzing the modal shape evolution.
-
(7)
As the axial pre-stretch reduces from axial pre-extension to axial pre-compression, the lowest vibration frequency of non-axisymmetric vibrations continues to go down and eventually reaches zero, indicating the 3D non-axisymmetric buckling instability. The critical axial pre-stretch increases with the radial voltage, which destabilizes the SEA tube.
-
(8)
Regardless of the vibration type, the vibration frequencies at first decrease rapidly with an increase of the value of the Gent parameter, and then gradually tend to the results predicted by the neo-Hookean model. The strain-stiffening effect manifests itself earlier at relatively high voltages.
Our results indicate the possibility of exploiting the electro-mechanical biasing fields to realize on-demand tunability of small-amplitude vibration behaviors of SEA tubes exhibiting the strain-stiffening effect. The present study provides guidelines for further experimental research and design of SEA tube-based electrostatically tunable resonant systems, which may find a wide range of potential applications, including active vibration isolators, tunable energy harvesters, tunable sound generators, as well as biomedical actuators and sensors.
It is emphasized that the main focus of this work is to investigate the influence of strain-stiffening effect on the vibration behaviors of SEA tubes. Note that Arora et al. (2022) explored the influence of strain-stiffening effect on the stability of soft periodic laminates and achieved deformation-activated negative group velocity by customizing the stiffening behavior of the non-Gaussian soft phases. Thus, how the strain-stiffening effect modulates the buckling instability and wave characteristics of negative group velocity for SEA structures is a topic worthy of further investigation.
Acknowledgments
The work is supported by the National Natural Science Foundation of China (Nos. 11872329, 12192210, 12192211, 12272339, and 12072315), the 111 Project, PR China (No. B21034), the Natural Science Foundation of Zhejiang Province, PR China (No. LD21A020001), the research project from Huanjiang Laboratory (Zhuji, Zhejiang Province), and the China Scholarship Council (No. 202006320363). WB gratefully acknowledges the support of the European Union Horizon 2020 Research and Innovation Programme under the Marie Skłodowska-Curie Actions (Grant No. 896229).
Appendix A Frequency equations of free vibrations in a pre-stretched hyperelastic tube characterized by the Gent model
In this appendix, we use the conventional displacement method to derive the frequency equations of four kinds of vibrations (i.e. non-axisymmetric vibrations, L vibrations (including the breathing mode), T vibrations and prismatic vibrations) in a pre-stretched hyperelastic tube characterized by the Gent model. For arbitrary energy function models, the detailed derivations of the frequency equations of the first three kinds of vibrations in a pre-stretched hyperelastic tube have been provided in our previous work (Zhu et al., 2020).
Three displacement functions , , and to express the displacement components are introduced as:
(57) |
where , , and have the following assumed formal solutions:
(58) |
The specific derivation process for the frequency equation of the non-axisymmetric and axisymmetric vibrations can be found in the Appendix C in the paper of Zhu et al. (2020) and is omitted here for brevity. Specifically, we can refer to the frequency equations (C.11) and (C.12) for the non-axisymmetric vibrations and (C.13) for the axisymmetric vibrations (including the L vibrations and T vibrations) in that paper.
Now consider an axially pre-stretched hyperelastic tube characterized by the incompressible Gent model with strain-energy function Eq. (17) with . We obtain the required effective material parameters as:
(59) |
where , , and for the homogeneous deformation in the hyperelastic tube without the electro-mechanical coupling. According to Eq. (9)1, the Lagrange multiplier can be obtained from as
(60) |
A.1 Breathing mode
For the breathing mode with and , we recall the frequency equation from Eq. (C.20) in the paper of Zhu et al. (2020),
(61) |
which, when combined with Eqs. (59) and (60), results in the following frequency equation of the breathing mode for the pre-stretched Gent hyperelastic tube:
(62) |
where and . Thus, the vibration frequency of the breathing mode of the Gent hyperelastic tube depends on the inner-to-outer radius ratio and axial pre-stretch appearing in , but is independent of the length-to-thickness ratio . When is large enough, the effect of axial pre-stretch can be counteracted, as with the result of the neo-Hookean model (Zhu et al., 2020).
A.2 Purely torsional vibrations (T vibrations)
For the T vibrations with and , the only non-zero displacement component is . Inserting Eq. (59) into Eq. (C.13) of Zhu et al. (2020), the elements (, , , ) related to the T vibrations can be rewritten as
(63) |
where , with and , and and are the Bessel functions of the first and second kinds of order , respectively. Thus, we obtain the frequency equation of the T vibrations for the Gent hyperelastic tube as
(64) |
It is clear from Eq. (64) that the natural frequency of the T vibrations depends on the axial pre-stretch in , the length-to-thickness ratio in and the inner-to-outer radius ratio . Furthermore, we note from Eqs. (64) and (59) that is one of the solutions to the frequency equation of the T vibrations, which yields
(65) |
The solution (65) represents the torsional displacement proportional to the radius (i.e., each cross-section of the tube rotates around its center during the vibration), and we see that its natural frequency depends on the axial pre-stretch and length-to-thickness ratio, but is independent of the inner-to-outer radius ratio.
A.3 Prismatic vibrations
For the case of prismatic vibrations with and , we have and its corresponding assumed solutions can be rewritten from Eq. (58) as
(66) |
Moreover, its corresponding incremental governing equations (see Eq. (C.4) of Zhu et al. (2020)) reduce to
(67) |
where and . Obviously, Eq. (67)2 is a homogeneous Euler equation, with general solution
(68) |
where and are the undetermined constants. Substituting Eq. (67)2 into Eq. (67)3 and using Eq. (68), we obtain the solution of as
(69) |
In addition, Eq. (67)1 is clearly a Bessel equation of order , and its solution is
(70) |
where and are arbitrary constants to be determined. Substituting Eqs. (68)-(70) into Eqs. (57) and (37)1,4,8, we obtain the non-zero incremental transverse stress components as
(71) |
where
(72) |
in which the prime denotes differentiation with respect to .
To satisfy the imposed mechanical boundary conditions (52)2 and ensure that the determinantal condition for non-trivial solutions exists, we obtain the frequency equation of the prismatic vibrations as
(73) |
where the first two rows of the determinant’s elements that correspond to the boundary conditions on the inner surface are written as
(74) |
For the mechanical boundary conditions on the outer surface , we replace the inner radius with the outer radius in Eq. (74) to obtain the elements of the final two rows of the determinant Eq. (73).
References
- Anderson et al. (2012) Anderson, I.A., Gisby, T.A., McKay, T.G., O’Brien, B.M., Calius, E.P., 2012. Multi-functional dielectric elastomer artificial muscles for soft and smart machines. J. Appl. Phys. 112, 041101.
- Arora et al. (2022) Arora, N., Yao, Q., Rudykh, S., 2022. Deformation activated negative group velocity state in soft laminates. Extreme Mech. Lett. 51, 101592.
- Baesu et al. (2003) Baesu, E., Fortune, D., Soós, E., 2003. Incremental behaviour of hyperelastic dielectrics and piezoelectric crystals. Z. Angew. Math. Phys. 54, 160–178.
- Baumhauer and Tiersten (1973) Baumhauer, J., Tiersten, H., 1973. Nonlinear electroelastic equations for small fields superposed on a bias. J. Acoust. Soc. Am. 54, 1017–1034.
- Bertoldi and Gei (2011) Bertoldi, K., Gei, M., 2011. Instabilities in multilayered soft dielectrics. J. Mech. Phys. Solids 59, 18–42.
- Bortot and Shmuel (2018) Bortot, E., Shmuel, G., 2018. Prismatic bifurcations of soft dielectric tubes. Int. J. Eng. Sci. 124, 104–114.
- Broderick et al. (2020) Broderick, H.C., Dorfmann, L., Destrade, M., 2020. Electro-elastic Lamb waves in dielectric plates. Extreme Mech. Lett. 39, 100782.
- Cao et al. (2024) Cao, Y.Z., Wu, B., Carrera, E., Rudykh, S., Chen, W.Q., 2024. Axisymmetric vibration of multilayered electroactive circular plates in contact with fluid. J. Sound Vib. 573, 118189.
- Carpi et al. (2011) Carpi, F., De Rossi, D., Kornbluh, R., Pelrine, R.E., Sommer-Larsen, P., 2011. Dielectric Elastomers as Electromechanical Transducers: Fundamentals, Materials, Devices, Models and Applications of an Emerging Electroactive Polymer Technology. Elsevier.
- Chen and Ding (2002) Chen, W.Q., Ding, H.J., 2002. On free vibration of a functionally graded piezoelectric rectangular plate. Acta Mech. 153, 207–216.
- Chen et al. (2017) Chen, W.Q., Liu, D.Y., Kitipornchai, S., Yang, J., 2017. Bifurcation of pressurized functionally graded elastomeric hollow cylinders. Compos. B: Eng. 109, 259–276.
- Chen et al. (2022) Chen, Y.J., Wu, B., Destrade, M., Chen, W.Q., 2022. Voltage-controlled topological interface states for bending waves in soft dielectric phononic crystal plates. Int. J. Solids Struct. 259, 112013.
- Chen et al. (2020) Chen, Y.J., Wu, B., Su, Y.P., Chen, W.Q., 2020. Effects of strain stiffening and electrostriction on tunable elastic waves in compressible dielectric elastomer laminates. Int. J. Mech. Sci. 176, 105572.
- Ding and Chen (2001) Ding, H.J., Chen, W.Q., 2001. Three Dimensional Problems of Piezoelasticity. Nova Science Publishers.
- Dorfmann and Ogden (2006) Dorfmann, A., Ogden, R.W., 2006. Nonlinear electroelastic deformations. J. Elast. 82, 99–127.
- Dorfmann and Ogden (2010) Dorfmann, A., Ogden, R.W., 2010. Electroelastic waves in a finitely deformed electroactive material. IMA J. Appl. Math. 75, 603–636.
- Dorfmann and Ogden (2014) Dorfmann, L., Ogden, R.W., 2014. Nonlinear Theory of Electroelastic and Magnetoelastic Interactions. Springer, New York.
- Dorfmann and Ogden (2020) Dorfmann, L., Ogden, R.W., 2020. Waves and vibrations in a finitely deformed electroelastic circular cylindrical tube. Proc. R. Soc. A: Math. Phys. Eng. Sci. 476, 20190701.
- Dubois et al. (2008) Dubois, P., Rosset, S., Niklaus, M., Dadras, M., Shea, H., 2008. Voltage control of the resonance frequency of dielectric electroactive polymer (DEAP) membranes. J. Microelectromech. Syst. 17, 1072–1081.
- Eringen and Maugin (1990) Eringen, A.C., Maugin, G.A., 1990. Electrodynamics of Continua I: Foundations and Solid Media. Springer-Verlag, New York.
- Fan and Zhang (1992) Fan, J.R., Zhang, J.Y., 1992. Exact solutions for thick laminated shells. Sci. China A 35, 1343–1355.
- Fu et al. (2018) Fu, Y.B., Dorfmann, L., Xie, Y.X., 2018. Localized necking of a dielectric membrane. Extreme Mech. Lett. 21, 44–48.
- Galich and Rudykh (2015) Galich, P.I., Rudykh, S., 2015. Influence of stiffening on elastic wave propagation in extremely deformed soft matter: from nearly incompressible to auxetic materials. Extreme Mech. Lett. 4, 156–161.
- Galich and Rudykh (2016) Galich, P.I., Rudykh, S., 2016. Manipulating pressure and shear waves in dielectric elastomers via external electric stimuli. Int. J. Solids Struct. 91, 18–25.
- Galich and Rudykh (2017) Galich, P.I., Rudykh, S., 2017. Shear wave propagation and band gaps in finitely deformed dielectric elastomer laminates: long wave estimates and exact solution. J. Appl. Mech. 84, 091002.
- Gent (1996) Gent, A.N., 1996. A new constitutive relation for rubber. Rubber Chem. Technol. 69, 59–61.
- Goshkoderia and Rudykh (2017) Goshkoderia, A., Rudykh, S., 2017. Electromechanical macroscopic instabilities in soft dielectric elastomer composites with periodic microstructures. Eur. J. Mech. A/Solids 65, 243–256.
- Haughton and Ogden (1979) Haughton, D.M., Ogden, R.W., 1979. Bifurcation of inflated circular cylinders of elastic material under axial loading–II. Exact theory for thick-walled tubes. J. Mech. Phys. Solids 27, 489–512.
- Hosoya et al. (2015) Hosoya, N., Baba, S., Maeda, S., 2015. Hemispherical breathing mode speaker using a dielectric elastomer actuator. J. Acoust. Soc. Am. 138, EL424–EL428.
- Jandron and Henann (2018) Jandron, M., Henann, D.L., 2018. A numerical simulation capability for electroelastic wave propagation in dielectric elastomer composites: Application to tunable soft phononic crystals. Int. J. Solids Struct. 150, 1–21.
- Liu (2013) Liu, L.P., 2013. On energy formulations of electrostatics for continuum media. J. Mech. Phys. Solids 61, 968–990.
- Lu et al. (2015a) Lu, T.Q., An, L., Li, J.G., Yuan, C., Wang, T.J., 2015a. Electro-mechanical coupling bifurcation and bulging propagation in a cylindrical dielectric elastomer tube. J. Mech. Phys. Solids 85, 160–175.
- Lu et al. (2020) Lu, T.Q., Ma, C., Wang, T.J., 2020. Mechanics of dielectric elastomer structures: A review. Extreme Mech. Lett. 38, 100752.
- Lu et al. (2015b) Lu, Z.B., Godaba, H., Cui, Y.D., Foo, C.C., Debiasi, M., Zhu, J., 2015b. An electronically tunable duct silencer using dielectric elastomer actuators. J. Acoust. Soc. Am. 138, EL236–EL241.
- Mace and Manconi (2012) Mace, B.R., Manconi, E., 2012. Wave motion and dispersion phenomena: Veering, locking and strong coupling effects. J. Acoust. Soc. Am. 131, 1015–1028.
- Mao et al. (2019) Mao, R.W., Wu, B., Carrera, E., Chen, W.Q., 2019. Electrostatically tunable small-amplitude free vibrations of pressurized electro-active spherical balloons. Int. J. Non-Linear Mech. 117, 103237.
- Maugin (1988) Maugin, G.A., 1988. Continuum Mechanics of Electromagnetic Solids. North-Holland, Amsterdam.
- McMeeking and Landis (2005) McMeeking, R.M., Landis, C.M., 2005. Electrostatic forces and stored energy for deformable dielectric materials. J. Appl. Mech. 72, 581–590.
- Melnikov and Ogden (2016) Melnikov, A., Ogden, R.W., 2016. Finite deformations of an electroelastic circular cylindrical tube. Z. Angew. Math. Phys. 67, 1–20.
- Melnikov and Ogden (2018) Melnikov, A., Ogden, R.W., 2018. Bifurcation of finitely deformed thick-walled electroelastic cylindrical tubes subject to a radial electric field. Z. Angew. Math. Phys. 69, 1–27.
- Pelrine et al. (1998) Pelrine, R.E., Kornbluh, R.D., Joseph, J.P., 1998. Electrostriction of polymer dielectrics with compliant electrodes as a means of actuation. Sens. Actuator A: Phys. 64, 77–85.
- Psarra et al. (2019) Psarra, E., Bodelot, L., Danas, K., 2019. Wrinkling to crinkling transitions and curvature localization in a magnetoelastic film bonded to a non-magnetic substrate. J. Mech. Phys. Solids 133, 103734.
- Rudykh and Bertoldi (2013) Rudykh, S., Bertoldi, K., 2013. Stability of anisotropic magnetorheological elastomers in finite deformations: A micromechanical approach. J. Mech. Phys. Solids 61, 949–967.
- Rudykh et al. (2014) Rudykh, S., Bhattacharya, K., DeBotton, G., 2014. Multiscale instabilities in soft heterogeneous dielectric elastomers. Proc. R. Soc. A: Math. Phys. Eng. Sci. 470, 20130618.
- Sarban et al. (2011) Sarban, R., Jones, R., Mace, B., Rustighi, E., 2011. A tubular dielectric elastomer actuator: Fabrication, characterization and active vibration isolation. Mech. Syst. Signal Process. 25, 2879–2891.
- Shmuel (2013) Shmuel, G., 2013. Electrostatically tunable band gaps in finitely extensible dielectric elastomer fiber composites. Int. J. Solids Struct. 50, 680–686.
- Shmuel (2015) Shmuel, G., 2015. Manipulating torsional motions of soft dielectric tubes. J. Appl. Phys. 117, 174902.
- Shmuel and DeBotton (2013) Shmuel, G., DeBotton, G., 2013. Axisymmetric wave propagation in finitely deformed dielectric elastomer tubes. Proc. R. Soc. A: Math. Phys. Eng. Sci. 469, 20130071.
- Shmuel et al. (2012) Shmuel, G., Gei, M., Debotton, G., 2012. The Rayleigh–Lamb wave propagation in dielectric elastomer layers subjected to large deformations. Int. J. Non-Linear Mech. 47, 307–316.
- Shmuel and Pernas-Salomón (2016) Shmuel, G., Pernas-Salomón, R., 2016. Manipulating motions of elastomer films by electrostatically-controlled aperiodicity. Smart Mater. Struct. 25, 125012.
- Simpson and Spector (1984) Simpson, H.C., Spector, S.J., 1984. On barrelling instabilities in finite elasticity. J. Elast. 14, 103–125.
- Son and Goulbourne (2012) Son, S., Goulbourne, N.C., 2012. Large strain analysis of a soft polymer electromechanical sensor coupled to an arterial segment. J. Intell. Mater. Syst. Struct. 23, 575–586.
- Su (2020) Su, Y.P., 2020. Voltage-controlled instability transitions and competitions in a finitely deformed dielectric elastomer tube. Int. J. Eng. Sci. 157, 103380.
- Su et al. (2018) Su, Y.P., Broderick, H.C., Chen, W.Q., Destrade, M., 2018. Wrinkles in soft dielectric plates. J. Mech. Phys. Solids 119, 298–318.
- Su et al. (2016) Su, Y.P., Wang, H.M., Zhang, C.L., Chen, W.Q., 2016. Propagation of non-axisymmetric waves in an infinite soft electroactive hollow cylinder under uniform biasing fields. Int. J. Solids Struct. 81, 262–273.
- Su et al. (2019) Su, Y.P., Wu, B., Chen, W.Q., Destrade, M., 2019. Finite bending and pattern evolution of the associated instability for a dielectric elastomer slab. Int. J. Solids Struct. 158, 191–209.
- Su et al. (2020) Su, Y.P., Wu, B., Chen, W.Q., Destrade, M., 2020. Pattern evolution in bending dielectric-elastomeric bilayers. J. Mech. Phys. Solids 136, 103670.
- Sugimoto et al. (2013) Sugimoto, T., Ando, A., Ono, K., Morita, Y., Hosoda, K., Ishii, D., Nakamura, K., 2013. A lightweight push-pull acoustic transducer composed of a pair of dielectric elastomer films. J. Acoust. Soc. Am. 134, EL432–EL437.
- Suo et al. (2008) Suo, Z.G., Zhao, X.H., Greene, W.H., 2008. A nonlinear field theory of deformable dielectrics. J. Mech. Phys. Solids 56, 467–486.
- Tiersten (1981) Tiersten, H.F., 1981. Electroelastic interactions and the piezoelectric equations. J. Acoust. Soc. Am. 70, 1567–1576.
- Toupin (1956) Toupin, R.A., 1956. The elastic dielectric. J. Ration. Mech. Anal. 5, 849–915.
- Toupin (1963) Toupin, R.A., 1963. A dynamical theory of elastic dielectrics. Int. J. Eng. Sci. 1, 101–126.
- Wang et al. (2017) Wang, F.F., Yuan, C., Lu, T.Q., Wang, T.J., 2017. Anomalous bulging behaviors of a dielectric elastomer balloon under internal pressure and electric actuation. J. Mech. Phys. Solids 102, 1–16.
- Wang et al. (2020) Wang, Y.F., Wang, Y.Z., Wu, B., Chen, W.Q., Wang, Y.S., 2020. Tunable and active phononic crystals and metamaterials. Appl. Mech. Rev. 72, 040801.
- Wu and Destrade (2021) Wu, B., Destrade, M., 2021. Wrinkling of soft magneto-active plates. Int. J. Solids Struct. 208, 13–30.
- Wu et al. (2020) Wu, B., Destrade, M., Chen, W.Q., 2020. Nonlinear response and axisymmetric wave propagation in functionally graded soft electro-active tubes. Int. J. Mech. Sci. 187, 106006.
- Wu et al. (2017) Wu, B., Su, Y.P., Chen, W.Q., Zhang, C.Z., 2017. On guided circumferential waves in soft electroactive tubes under radially inhomogeneous biasing fields. J. Mech. Phys. Solids 99, 116–145.
- Wu et al. (2016) Wu, B., Zhang, C.L., Zhang, C.Z., Chen, W.Q., 2016. Theory of electroelasticity accounting for biasing fields: Retrospect, comparison and perspective. Adv. Mech. 46, 201601.
- Xie et al. (2016) Xie, Y.X., Liu, J.C., Fu, Y.B., 2016. Bifurcation of a dielectric elastomer balloon under pressurized inflation and electric actuation. Int. J. Solids Struct. 78, 182–188.
- Zhang et al. (2015) Zhang, J., Chen, H., Li, B., McCoul, D., Pei, Q., 2015. Tunable active vibration attenuation using highly deformable dielectric elastomers. Smart Mater. Struct. 24, 115033.
- Zhao and Suo (2007) Zhao, X.H., Suo, Z.G., 2007. Method to analyze electromechanical stability of dielectric elastomers. Appl. Phys. Lett. 91, 061921.
- Zhao and Wang (2014) Zhao, X.H., Wang, Q.M., 2014. Harnessing large deformation and instabilities of soft dielectrics: Theory, experiment, and application. Appl. Phys. Rev. 1, 021304.
- Zhao et al. (2016) Zhao, Z.H., Shuai, C.G., Gao, Y., Rustighi, E., Xuan, Y., 2016. An application review of dielectric electroactive polymer actuators in acoustics and vibration control. J. Phys. Conf. Ser. 744, 012162.
- Zhao et al. (2023) Zhao, Z.N., Chen, Y.J., Hu, X.Y., Bao, R.H., Wu, B., Chen, W.Q., 2023. Vibrations and waves in soft dielectric elastomer structures. Int. J. Mech. Sci. 239, 107885.
- Zhou et al. (2014) Zhou, J.Y., Jiang, L.Y., Khayat, R.E., 2014. Electromechanical response and failure modes of a dielectric elastomer tube actuator with boundary constraints. Smart Mater. Struct. 23, 045028.
- Zhu et al. (2020) Zhu, F.Z., Wu, B., Destrade, M., Chen, W.Q., 2020. Electrostatically tunable axisymmetric vibrations of soft electro-active tubes. J. Sound Vib. 483, 115467.
- Zhu et al. (2010a) Zhu, J., Cai, S.Q., Suo, Z.G., 2010a. Nonlinear oscillation of a dielectric elastomer balloon. Polym. Int. 59, 378–383.
- Zhu et al. (2010b) Zhu, J., Cai, S.Q., Suo, Z.G., 2010b. Resonant behavior of a membrane of a dielectric elastomer. Int. J. Solids Struct. 47, 3254–3262.
- Zhu et al. (2010c) Zhu, J., Stoyanov, H., Kofod, G., Suo, Z.G., 2010c. Large deformation and electromechanical instability of a dielectric elastomer tube actuator. J. Appl. Phys. 108, 074113.
- Ziser and Shmuel (2017) Ziser, Y., Shmuel, G., 2017. Experimental slowing of flexural waves in dielectric elastomer films by voltage. Mech. Res. Commun. 85, 64–68.
- Zurlo et al. (2017) Zurlo, G., Destrade, M., DeTommasi, D., Puglisi, G., 2017. Catastrophic thinning of dielectric elastomers. Phys. Rev. Lett. 118, 078001.