A numerical framework for modeling anisotropic dielectric elastomers
Introduction
Dielectric elastomers (DEs) are compliant elastic materials and belong to a class of electro-active polymers (EAPs). DEs undergo large deformation when subjected to an electric field [1]. The ability of dielectric elastomers to achieve large deformation with their lightweight, low cost and fast response time [2] have made them attractive for a wide range of applications, such as energy harvesting devices [3], resonators [4], artificial muscles [5], active lenses [6], among many others. However, requirement of high voltage for achieving large strain together with the associated electromechanical instability has limited the widespread use of DE actuators. Several approaches have been proposed so far to circumvent these problems [[7], [8], [9], [10]]. One of such approaches is to impart anisotropy in the mechanical behavior of DE actuators [[11], [12]], which can be realized by various means and at various scales. A significant volume of previous investigations expounds on the electromechanical analysis of electroactive polymer composites, both theoretically [13] as well as experimentally [14]. In this approach, the deformation of the DE actuator is constrained using stiff fibers in the direction that is perpendicular to the desired direction of actuation [[8], [15], [16]]. Using this approach, Lu et al. [14] achieved 28% unidirectional actuation strain in the flat dielectric membranes. The computational modeling these fiber-stiffened DE actuators is performed by employing the simplistic lumped parameter models based on the assumption of homogeneous deformation [17], or using homogenization approaches [[18], [19]] facilitating a numerical solution through finite element method [20]. In another instance, a few recent experiments demonstrate the enhancement in the actuation response of the dielectric elastomers using the mechanical anisotropy of the electrode material [[21], [22], [23]]. With this approach, Cakmak et al. [21] achieved the linear actuation strain of more than 40% at a relatively low electric field (100 V/) using highly aligned carbon nanotube (CNT) sheet electrodes on acrylate adhesive films. In addition to fiber-reinforcement and imparting anisotropy in the material used for preparing the electrodes, there has been a continued interest in yet another way of enhancing the actuation by including high dielectric constant fillers into the polymer matrix [[24], [25], [26], [27], [28]]. Because of these inclusions, soft elastomers show anisotropy in both mechanical as well as dielectric behaviors. The direction of anisotropy of such composite material depends on the direction of the electric field applied during the curing process. The computational modeling of this dispersion anisotropy in dielectric elastomer actuators has received an increasing attention in the recent past [[26], [29], [30]].
In connection to the modeling and simulation of dielectric elastomers, the basic theory of non-linear electroelasticity was developed by Toupin [31], for investigating the quasistatic electromechanical behavior of elastic dielectrics. Recently, several researchers reformulated the electroelasticity theory of finite deformations and made significant efforts towards the numerical implementation for modeling the soft dielectrics [[32], [33], [34], [35], [36], [37], [38], [39], [40], [41], [42], [43]]. The electromechanical behavior and pull-in instability phenomena in finitely deformed isotropic dielectric elastomers, considering the effects of several factors, such as the level of prestress [[44], [45]], time dependent actuation [[46], [47], [48]], and viscoelasticity [[49], [50]], have been reported analytically as well as numerically in several investigations reported previously in the literature. However, limited efforts have been made to model numerically the electromechanical behavior of the DE actuators with dispersion type anisotropy. In this regard, Yong et al. [51] reported a lumped parameter model to investigate the electromechanical instability in transversely isotropic elastomers. Bustamante [52] presented a constitutive model that considers transverse isotropy stemming from the particle reinforced anisotropy. A homogenization approach for electroactive polymers with random particulate microstructure, in the limit of infinitesimal deformations, is presented by Siboni and Castaneda [18]. Recently, Hossain et al. [26] presented a constitutive framework for the modeling of dispersion type anisotropy in electroactive polymers. In view of the growing interest in exploiting anisotropy in the technology of DE actuators, it is imperative to develop a numerical modeling framework to support the design and development of this class of actuators. However, to the best of the authors’ knowledge, a finite element-based numerical framework for modeling the anisotropic dielectric elastomers has not hitherto been discussed in the literature.
To this end, this paper presents a nonlinear, finite deformation, finite element framework for anisotropic dielectric elastomers based on the existing models for incompressible anisotropic neo-Hookean hyperelastic solids [53] and ideal dielectric elastomers [54]. The analytical expressions are derived for the tangent moduli of the anisotropic materials. This material model has been implemented into an in-house finite element code. We propose a computationally efficient staggered algorithm, in which, we decouple the displacement and the electric potential fields and solve for them separately, to avoid the ill-conditioning that arises in the monolithic scheme of solution. A parametric study is performed for investigating the effect of material anisotropy on the pull-in instability of dielectric elastomers with different levels of pre-stretch. Further, we simulate the effect of material anisotropy on the actuation performance of the bending and buckling actuators.
This paper is organized in six sections. In Section 2, we overview the basics of nonlinear continuum electromechanics and the finite element implementation of the coupled electro-elastic problem. A material model for the anisotropic dielectric elastomers and the analytic expressions for tangent moduli are presented in Section 3. A computationally efficient staggered algorithm for the solution of coupled nonlinear equations is proposed in Section 4. In Section 5, we demonstrate the utility of the numerical framework by simulating the effect of material anisotropy on (1) the pull-in instability of the dielectric elastomers, and (2) actuation performance of bimorph bending and buckling actuators. Subsequently, a comparison of the computational efficiency of the proposed staggered algorithm with that of the monolithic approach is presented. The conclusions drawn from the present study are summarized in Section 6.
Section snippets
Basics of non-linear continuum electro-mechanics and finite element model
In this section, we begin by briefly describing the basic governing equations and finite element formulation pertaining to the electrostatic deformation of dielectric elastomers, based on the nonlinear field theory of deformable dielectrics [[31], [32], [33], [34]].
Let be the undeformed configuration of a continuum dielectric body at time . The cartesian position vector denotes an arbitrary material particle of in the undeformed configuration. At time , the deformed
Material model for anisotropic dielectric elastomers
In this section, we present a material model for investigating the anisotropic deformation of dielectric elastomers following the same approach in isotropic dielectric elastomer [54]. We decompose total free energy density into isotropic, anisotropic and electrical energy parts. In this paper, we restrict our attention to transversely isotropic and nearly incompressible hyperelastic materials. In this case, the free energy density function is represented as
A staggered algorithm for simulation of anisotropic dielectric elastomers
For dielectric elastomers the value of dielectric permittivity is of the order of 10−12 F/m, while the bulk-modulus is of the order of N/m. Because of this, the order of mechanical stiffness matrix, i.e., , in Eq. (9) differs several orders of magnitude from the coupling and electrical stiffness matrices, i.e., and . For better condition numbers, we employ a staggered algorithm, in which we decouple the displacement and the electric potential fields and solve for them
Numerical results
The finite-element method and the staggered algorithm for anisotropic dielectric elastomers discussed in the previous section have been implemented using an in-house MATLAB code. In all FE simulations performed in this section, standard displacement-based eight-node hexahedral elements with four degrees of freedom (3 mechanical and 1 electrical) per node have been used. To enforce the condition of incompressibility, we set the bulk modulus to five orders of magnitude higher than the ground
Conclusions
In summary, we have presented a nonlinear finite element framework of dielectric elastomers that incorporates the effect of material anisotropy. A material model for anisotropic dielectric elastomers is constructed by adopting the nonlinear electroelasticity and incompressible anisotropic neo-Hookean material model. The analytical expressions are presented for the tangent moduli. A staggered algorithm is introduced to solve the coupled fields. The effect of anisotropy parameters on the pull-in
Acknowledgments
The authors gratefully acknowledge the financial assistance provided by the Department of Science and Technology, Government of India through Grant No: EMR/2017/003289. The authors wish to thank Prof. D. M. Joglekar, IIT Roorkee for many useful discussions. The authors are grateful to the anonymous reviewers for their insightful comments.
References (63)
- et al.
High-speed electrically actuated elastomers with strain greater than 100%
Science
(2000) - et al.
Dielectric Elastomers as Electromechanical Transducers: Fundamentals, Materials, Devices, Models and Applications of an Emerging Electroactive Polymer Technology
(2011) - et al.
Dielectric elastomers: Stretching the capabilities of energy harvesting
MRS Bull.
(2012) - et al.
Electromechanical and dynamic analyses of tunable dielectric elastomer resonator
Int. J. Solids Struct.
(2012) - et al.
Advances in dielectric elastomers for actuators and artificial muscles
Macromol. Rapid Commun.
(2010) - et al.
Electromechanically driven variable-focus lens based on transparent dielectric elastomer
Appl. Opt.
(2012) - et al.
A modulated voltage waveform for enhancing the travel range of dielectric elastomer actuators
J. Appl. Mech.
(2018) - et al.
Enhanced electroactive response of unidirectional elastomeric composites with high-dielectric-constant fibers
Adv. Mater.
(2014) - et al.
Interpenetrating polymer networks for high-performance electroelastomer artificial muscles
Adv. Mater.
(2006) - et al.
Bioinspired bicipital muscle with fiber-constrained dielectric elastomer actuator
Extreme Mech. Lett.
(2016)
All-organic dielectric-percolative three-component composite materials with high electromechanical response
Appl. Phys. Lett.
High performance electroactive polymers and nano-composites for artificial muscles
J. Intell. Mater. Syst. Struct.
Analysis of microstructural induced enhancement of electromechanical coupling in soft dielectrics
Appl. Phys. Lett.
Dielectric elastomer actuators under equal-biaxial forces, uniaxial forces, and uniaxial constraint of stiff fibers
Soft Matter
Large, uni-directional actuation in dielectric elastomers achieved by fiber stiffening
Appl. Phys. Lett.
Use of aligned fibers to enhance the performance of dielectric elastomer inchworm robots
Voltage-induced torsion of a fiber-reinforced tubular dielectric elastomer actuator
Compos. Sci. Technol.
A finite-strain constitutive theory for electro-active polymer composites via homogenization
Int. J. Non-Linear Mech.
Micro-electromechanics of soft dielectric matrix composites
Int. J. Solids Struct.
Computational stability analysis of periodic electroactive polymer composites across scales
Comput. Methods Appl. Mech. Engrg.
Carbon nanotube sheet electrodes for anisotropic actuation of dielectric elastomers
Carbon
Surface modification of anisotropic dielectric elastomer actuators with uni-and bi-axially wrinkled carbon electrodes for wettability control
Sci. Rep.
Enhanced anisotropic response of dielectric elastomer actuators with microcombed and etched carbon nanotube sheet electrodes
Carbon
Improvement of electromechanical actuating performances of a silicone dielectric elastomer by dispersion of titanium dioxide powder
IEEE Trans. Dielectr. Electr. Insul.
Dielectric study of poly (styrene-co-butadiene) composites with carbon black, silica, and nanoclay
Macromolecules
Modelling electro-active polymers with a dispersion-type anisotropy
Smart Mater. Struct.
Dielectric properties of silicone rubber-titanium dioxide composites prepared by dielectrophoretic assembly of filler particles
Smart Mater. Struct.
A multi-scale approach to model the curing process in magneto-sensitive polymeric materials
Int. J. Solids Struct.
Modelling of iron-filled magneto-active polymers with a dispersed chain-like microstructure
Eur. J. Mech. A Solids
Two-scale computational homogenization of electro-elasticity at finite strains
Comput. Methods Appl. Mech. Engrg.
The elastic dielectric
J. Ration. Mech. Anal.
Cited by (60)
Inverse dynamic design for motion control of soft machines driven by dielectric elastomer actuators
2024, International Journal of Mechanical SciencesA solid-shell model of hard-magnetic soft materials
2024, International Journal of Mechanical SciencesThermo-electro-mechanical effects on nonlinear dynamics of smart dielectric elastomer minimum energy structures
2024, European Journal of Mechanics, A/SolidsNonlinear electromechanical responses in multi-layered fiber-reinforced dielectric elastomer composites
2024, Thin-Walled StructuresAdversarial deep energy method for solving saddle point problems involving dielectric elastomers
2024, Computer Methods in Applied Mechanics and EngineeringElectro-mechanical actuation modulates fracture performance of soft dielectric elastomers
2024, International Journal of Engineering Science