<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Robot. AI</journal-id>
<journal-title>Frontiers in Robotics and AI</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Robot. AI</abbrev-journal-title>
<issn pub-type="epub">2296-9144</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/frobt.2018.00137</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Robotics and AI</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Toward a Dielectric Elastomer Resonator Driven Flapping Wing Micro Air Vehicle</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name><surname>Cao</surname> <given-names>Chongjing</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/670545/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Burgess</surname> <given-names>Stuart</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name><surname>Conn</surname> <given-names>Andrew T.</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
<xref ref-type="corresp" rid="c001"><sup>&#x0002A;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/290512/overview"/>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>Bristol Robotics Laboratory</institution>, <addr-line>Bristol</addr-line>, <country>United Kingdom</country></aff>
<aff id="aff2"><sup>2</sup><institution>Department of Aerospace Engineering, University of Bristol</institution>, <addr-line>Bristol</addr-line>, <country>United Kingdom</country></aff>
<aff id="aff3"><sup>3</sup><institution>Department of Mechanical Engineering, University of Bristol</institution>, <addr-line>Bristol</addr-line>, <country>United Kingdom</country></aff>
<author-notes>
<fn fn-type="edited-by"><p>Edited by: Perla Maiolino, University of Cambridge, United Kingdom</p></fn>
<fn fn-type="edited-by"><p>Reviewed by: Claudio Rossi, Polytechnic University of Madrid, Spain; Massimiliano Gei, Cardiff University, United Kingdom</p></fn>
<corresp id="c001">&#x0002A;Correspondence: Andrew T. Conn <email>a.conn&#x00040;bristol.ac.uk</email></corresp>
<fn fn-type="other" id="fn001"><p>This article was submitted to Soft Robotics, a section of the journal Frontiers in Robotics and AI</p></fn></author-notes>
<pub-date pub-type="epub">
<day>23</day>
<month>01</month>
<year>2019</year>
</pub-date>
<pub-date pub-type="collection">
<year>2018</year>
</pub-date>
<volume>5</volume>
<elocation-id>137</elocation-id>
<history>
<date date-type="received">
<day>21</day>
<month>09</month>
<year>2018</year>
</date>
<date date-type="accepted">
<day>17</day>
<month>12</month>
<year>2018</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x000A9; 2019 Cao, Burgess and Conn.</copyright-statement>
<copyright-year>2019</copyright-year>
<copyright-holder>Cao, Burgess and Conn</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p></license>
</permissions>
<abstract><p>In the last two decades, insect-inspired flapping wing micro air vehicles (MAVs) have attracted great attention for their potential for highly agile flight. Insects flap their wings at the resonant frequencies of their flapping mechanisms. Resonant actuation is highly advantageous as it amplifies the flapping amplitude and reduces the inertial power demand. Emerging soft actuators, such as dielectric elastomer actuators (DEAs) have large actuation strains and thanks to their inherent elasticity, DEAs have been shown a promising candidate for resonant actuation. In this work a double cone DEA configuration is presented, a mathematic model is developed to characterize its quasi-static and dynamic performance. We compare the high frequency performance of two most common dielectric elastomers: silicone elastomer and polyacrylate tape VHB. The mechanical power output of the DEA is experimentally analyzed as a DEA-mass oscillator. Then a flapping wing mechanism actuated by this elastic actuator is demonstrated, this design is able to provide a peak flapping amplitude of 63&#x000B0; at the frequency of 18 Hz.</p></abstract>
<kwd-group>
<kwd>flying robots</kwd>
<kwd>bio-inspired robotics</kwd>
<kwd>insect flight</kwd>
<kwd>electro-active polymer</kwd>
<kwd>dielectric elastomer actuator</kwd>
<kwd>resonance</kwd>
</kwd-group>
<counts>
<fig-count count="13"/>
<table-count count="3"/>
<equation-count count="6"/>
<ref-count count="38"/>
<page-count count="11"/>
<word-count count="7617"/>
</counts>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="s1">
<title>Introduction</title>
<p>In the last two decades, insect-inspired flapping wing micro air vehicles (MAVs) have attracted significant research interest for their potential for highly agile flight. Many MAV designs have been developed, such as Microrobotic Fly , DelFly (De Croon et al., <xref ref-type="bibr" rid="B13">2009</xref>) and Robotic hummingbird (Keennon et al., <xref ref-type="bibr" rid="B20">2012</xref>). One challenge all flapping wing MAV researchers have been facing is the extremely high power demands required for autonomous flight at micro scales. In nature, insects solve this problem by taking advantage of their elastic thorax and muscle system as a damped oscillator and flap their wings at its resonant frequency . Most species of flying insects utilize indirect flight muscles called dorsoventral (dvm) and dorsolateral (dlm) muscles to elevate and depress their wings, respectively, as illustrated in Figures <xref ref-type="fig" rid="F1">1A,B</xref>. The muscles are described as indirect as instead of directly driving the wings, they deform the highly elastic notum, the top plate of the thorax, which then drives the wings through the pleural wing process. Insects naturally excite this natural oscillator at its resonant frequency which amplifies the flapping stroke and greatly reduces inertial power demands . In conventional robotic technologies, actuators and power transmission mechanisms are all rigid, in order to incorporate elasticity into the MAV designs, additional elastic elements have to be added, which increase the weight and complexity of the system. To date, only few studies have utilized resonant excitation in flapping wing MAV designs [see e.g., (Baek et al., <xref ref-type="bibr" rid="B1">2009</xref>; Bolsman et al., <xref ref-type="bibr" rid="B2">2009</xref>; Zhang and Deng, <xref ref-type="bibr" rid="B38">2017</xref>)]. By adding a spring to a motor driven flapping wing MAV, Baek et al. (<xref ref-type="bibr" rid="B1">2009</xref>) showed a 30% average power reduction by driving the MAV at its resonance. Microrobotic Fly utilized the inherent elasticity of piezoelectric actuator and achieved a tethered flight at the resonant frequency of this system of 110 Hz. A comprehensive review on flapping wing MAV designs with integrated elastic elements can be found in Zhang and Rossi (<xref ref-type="bibr" rid="B37">2017</xref>).</p>
<fig id="F1" position="float">
<label>Figure 1</label>
<caption><p>Transverse view of the thorax of a typical indirect flight insect. The dorsoventral (dvm) and dorsolateral (dlm) muscles elevating <bold>(A)</bold> and depressing <bold>(B)</bold> the wings indirectly via the notum and pleural wing process (redrawn from Pringle, <xref ref-type="bibr" rid="B26">1957</xref>).</p></caption>
<graphic xlink:href="frobt-05-00137-g0001.tif"/>
</fig>
<p>Emerging soft actuators, such as dielectric elastomer actuators (DEAs) offer an alternative paradigm for designing flapping wing MAVs. Firstly, as an actuator, a DEA can generate large linear actuation strains that are comparable and even greater than muscles. Secondly, as an elastomer, the inherent elasticity of this material makes it capable of storing and releasing elastic energy, which is like the highly efficient elastomeric protein (resilin) found in the thorax of insects (Dudley, <xref ref-type="bibr" rid="B14">2002</xref>). The large linear stroke and the inherent elasticity make DEA ideal for a novel flapping wing MAV design that mimics the resonant actuation found in insects without the complex transmission mechanisms and any additional elastic elements found in conventional flapping wing MAV designs. Another advantage of DEAs for MAV applications is the good scaling capability. Miniaturizing DEAs can be fairly easy, DEA membranes with the thickness of hundreds of nanometers have been reported (T&#x000F6;pper et al., <xref ref-type="bibr" rid="B30">2015</xref>; Weiss et al., <xref ref-type="bibr" rid="B34">2016</xref>), and multiple layers of DEA membranes can be stacked to generate a lager force and power output (Carpi et al., <xref ref-type="bibr" rid="B10">2007</xref>; Kovacs et al., <xref ref-type="bibr" rid="B21">2009</xref>). To date, only few DEA actuated flapping wing MAVs have been developed [see e.g., (Burgess et al., <xref ref-type="bibr" rid="B5">2009</xref>; Henke et al., <xref ref-type="bibr" rid="B17">2017</xref>; Lau et al., <xref ref-type="bibr" rid="B22">2017</xref>)]. Lau and his co-workers (Lau et al., <xref ref-type="bibr" rid="B23">2014</xref>, <xref ref-type="bibr" rid="B22">2017</xref>) have developed rolled and stacked DEAs to flap the wings and these two designs demonstrated a flapping stroke of 10 and 2&#x000B0;, respectively at 1 Hz. Henke et al. (<xref ref-type="bibr" rid="B17">2017</xref>) have demonstrated a dragonfly-like flapping robot using a minimum energy structure mechanism and the flapping stroke is estimated to be about 15&#x000B0; at a few hertz. It should be noted that no aforementioned DEA-driven flapping mechanisms has utilized resonant actuation and the generated stroke and frequency are far from being high enough to enable flight. In contrast, insects have a typical flapping stroke of 120&#x000B0; and frequencies of 20 to over 200 Hz depending on the size of the species (Brodsky, <xref ref-type="bibr" rid="B4">1994</xref>). In this work, we present a DEA-driven flapping wing MAV design which seeks to achieve a larger stroke amplitude and wingbeat frequency by utilizing the resonance actuation of DEAs. The concept of this design is illustrated in Figures <xref ref-type="fig" rid="F2">2A,B</xref> which is clearly inspired by the insect elastic thorax.</p>
<fig id="F2" position="float">
<label>Figure 2</label>
<caption><p><bold>(A)</bold> Proposed flapping mechanism where the antagonistic linear double cone DEA is attached to the frame and drives the wings via a transmission mechanism. <bold>(B)</bold> Schematic design of a double cone actuator with the inner frame radius, <italic>b</italic>, and spacer height, <italic>L</italic>, labeled.</p></caption>
<graphic xlink:href="frobt-05-00137-g0002.tif"/>
</fig>
<p>This paper is organized as follows. In section Double Cone Linear DEA Design, the design of the double cone linear actuator is introduced. The performance of two types of elastomers for DEAs in high frequency domain at which insect inspired MAV operate is compared. Analytical model is developed to predict the quasi-static performance of the double cone DEA and a spring-mass model is used to characterize its natural frequency. In section Power Output of the Double Cone DEA, we investigate the mechanical power output of the DEA as a function of actuation frequency. Then in section Flapping Wing MAV Design the resonant actuation of the DEA on a flapping mechanism is investigated. Finally, conclusions and future work are discussed.</p>
</sec>
<sec id="s2">
<title>Double Cone Linear DEA Design</title>
<sec>
<title>Design Overview</title>
<p>In this work, a double cone DEA configuration is employed as was developed in Choi et al. (<xref ref-type="bibr" rid="B11">2003</xref>) and Conn and Rossiter (<xref ref-type="bibr" rid="B12">2012</xref>). The general concept involves two dielectric elastomer membranes being separately bonded to two circular frames with a protrusion from a strut deforms the center of the membranes out-of-plane to form a double cone shape, as shown in Figure <xref ref-type="fig" rid="F2">2</xref>. The two membranes can be actuated separately to achieve bidirectional actuation. Apart from the advantage of natural agonist-antagonist configuration, the double cone DEA design also has a good compactness, high mass-specific energy density and can be fabricated consistently. The actuation principle of a double cone DEA is explained as follows. Due to the similar geometry and stretch ratio, top and bottom membranes exert an equal reaction force on the strut, which makes the strut balance in the middle (assuming the weight of the strut is negligible). However, when an electric field is applied on one membrane, the generated Maxwell stress reduces the force exerted by this membrane, the imbalance in forces drives the strut toward the actuated membrane side until another force equilibrium is achieved. As the Maxwell stress is directly related to the electric field applied to the DEA, electric fields in the range of 50&#x02013;100 V/&#x003BC;m are usually used.</p>
</sec>
<sec>
<title>Actuator Material Characterization</title>
<p>Polyacrylate tape VHB (3M) and silicone elastomers are the two predominant DEA materials. VHB has the advantages of large energy density and actuation strain (Carpi et al., <xref ref-type="bibr" rid="B9">2011</xref>), its inherent adhesiveness and wide commercial availability also make it ideal for DEA prototyping. A potential drawback of VHB elastomer is its high viscoelasticity which limits the actuation bandwidth. On the other hand, silicone elastomers typically have a lower viscoelasticity, which makes them potentially more suitable for high frequency applications, such as MAVs where the frequency is usually over 10 Hz. The low viscoelasticity of silicone elastomers also makes them more likely to achieve resonance than VHB. To determine which material is more suitable for MAV applications, in this work, we compare the performance of VHB 4905 against an off-the-shelf silicone membrane (PARKER EAP 40 &#x003BC;m).</p>
<p>The fabrication process of the double cone DEA is described as follows. For VHB based DEA, a 0.5 mm thick VHB 4905 membrane was stretched biaxially by 4 &#x000D7; 4 and then bonded to a rigid acrylic frame with an inner radius <italic>b</italic>. For the silicone DEA, the elastomer membrane has an initial thickness of 40 &#x003BC;m, it was bonded to a 0.1 mm thickness Mylar ring using silicone adhesive (Smooth-On Sil-Poxy) without any pre-stretch and then attached to the acrylic frame. No pre-stretch was introduced to the Parker silicone elastomer as the out-of-plane deformation can introduce significant stretch on the membrane. In our previous study, the Parker silicone double cone DEAs with no pre-stretch showed the highest stroke and work output (Cao and Conn, <xref ref-type="bibr" rid="B6">2018</xref>). Second, carbon conductive grease (M.G. Chemicals Ltd) was used as the compliant electrodes and was applied using a fine brush. Copper tape was used to connect the compliant electrodes and high voltage cables. Finally, the two circular frames were connected together using nylon fasteners. One nylon spacer with a height <italic>L</italic> was used as the support strut to deform the DEA membranes. The DEA components are shown in Figure <xref ref-type="fig" rid="F2">2</xref>.</p>
<p>In this work, the dynamic performance of the double cone DEAs is characterized by the free stroke as a function of the excitation frequency. The experiment followed the protocol suggested in Carpi et al. (<xref ref-type="bibr" rid="B8">2015</xref>) and is described as follows. Two double cone DEA specimens were fabricated, one using Parker silicone elastomer and the other one using VHB 4905. The geometries of the two specimens are identical, where the radius of the actuator membrane <italic>b</italic> &#x0003D; 20 mm and the strut height <italic>L</italic> &#x0003D; 30 mm. The frames of the actuator were fixed to the testing rig, while leaving the nylon spacer free to move. For each DEA, the two membranes were driven by two 180&#x000B0; out-of-phase sinusoidal voltage actuation waves with the amplitude of 1,920 V for Parker silicone membranes and 2,500 V for VHB (which equates to a nominal electric field of 60 and 100 V/&#x003BC;m, respectively). A frequency sweep from 1 to 100 Hz with a step of 5 Hz (step was reduced to 2 Hz near resonance) were conducted and at each frequency, 50 cycles were repeated to allow the actuator to reach a steady state. The actuation signals were generated using MATLAB (MathWorks) and sent to a high voltage amplifier (Ultravolt 5HVA23-BP1) via a DAQ (National Instruments BNC-2111). A laser displacement sensor (LK-G152 and LKGD500, Keyence) was used to measure the displacement of the DEA.</p>
<p>The displacement amplitudes of the DEAs and the phase differences between the input (actuation voltage) and output (displacement) signals against excitation frequency are shown in Figures <xref ref-type="fig" rid="F3">3A,B</xref>. As can be seen in Figure <xref ref-type="fig" rid="F3">3A</xref>, at low frequencies, both silicone and VHB specimens have a stroke close to 2 mm. Two peaks can be observed for silicone DEA specimen at 20 and 52 Hz, respectively. The first peak has a lower amplitude of 3 mm while the second peak has the amplitude of 7.1 mm, and is the mechanical resonant frequency of this specimen. The two amplitude peaks for silicone cone DEAs also have been demonstrated by Rizzello et al. (<xref ref-type="bibr" rid="B27">2015a</xref>,<xref ref-type="bibr" rid="B28">b</xref>) and is believed to due to the non-linearity of the Maxwell pressure (Maxwell pressure is proportional to the square of the actuation voltage) and this complex electro-mechanical dynamic phenomenon will be investigated further in the future. The amplitude at the resonance of this silicone DEA is over 355% than that at low frequencies (&#x0003C;10 Hz). The amplitude of the VHB specimen reduces continuously with the increasing actuation frequency and no peak was observed, this is due to the significant viscoelasticity of the material. As pointed out by Carpi et al. (<xref ref-type="bibr" rid="B8">2015</xref>) that, for high viscous damping materials which have no amplitude peak, the resonant frequency can be determined by the value that causes a phase shift of 90&#x000B0;. From Figure <xref ref-type="fig" rid="F3">3B</xref>, a 90&#x000B0; phase shift is found at 52 Hz for both silicone and VHB specimens. Which indicates that 52 Hz is the mechanical resonant frequency for both Parker silicone and VHB samples. This experiment clearly shows that silicone elastomer is superior to VHB for this double cone DEA configuration at in terms of bandwidth and resonant actuation thanks to its low viscous damping. In later sections, Parker silicone elastomer will be adopted for double cone DEA design.</p>
<fig id="F3" position="float">
<label>Figure 3</label>
<caption><p>Frequency sweep experiment for Parker silicone and VHB double cone DEA. <bold>(A)</bold> Amplitudes of the silicone and VHB specimen against actuation frequency. <bold>(B)</bold> Phase difference of the input and output signals of the silicone and VHB specimen against actuation frequency.</p></caption>
<graphic xlink:href="frobt-05-00137-g0003.tif"/>
</fig>
</sec>
<sec>
<title>Double Cone DEA Model</title>
<p>In this section we present a double cone DEA mathematic model. Hodgins et al. (<xref ref-type="bibr" rid="B18">2014</xref>) and Rizzello et al. (<xref ref-type="bibr" rid="B27">2015a</xref>,<xref ref-type="bibr" rid="B28">b</xref>) have developed approximated mathematical models to characterize the performance of conical DEAs with biasing springs and biasing mass. In their simulations, two important simplification assumptions were made, the first being a truncated cone shape approximation and the second involves homogeneous stress distribution on the membrane. However, as have been shown by an analytical model developed by He et al. (<xref ref-type="bibr" rid="B16">2010</xref>), Wang et al. (<xref ref-type="bibr" rid="B32">2016</xref>), and Wang (<xref ref-type="bibr" rid="B33">2018</xref>), the membrane deformation is in fact non-truncated and the stress distribution on the membrane is inhomogeneous. By utilizing this analytical approach, Bortot and Gei (<xref ref-type="bibr" rid="B3">2015</xref>) attempted to maximum the harvested energy of a cone DE generator by tuning the pre-stretch, geometrical ratio and intensity of maximum external load. Here we adopt this approach which is able to characterize the complex deformation and stress distribution on a conical DEA. We extend this to characterize double cone DEAs quasi-statically, and the simulation is then verified against experiments. Due to the extremely high computational cost to use this approach in a dynamic manner, we then use a simplified mass-spring oscillator model to predict the dynamic performance of this actuator.</p>
<sec>
<title>Quasi-Static Model</title>
<p>As can be seen from Figure <xref ref-type="fig" rid="F2">2</xref>, the double cone configuration is in fact two pre-deformed single cones with an inertial load (Nylon spacer in this case). Due to the symmetry of the top and bottom cones, we will analyse a single cone first. A piece of elastomer membrane with an initial thickness <italic>T</italic> is first pre-stretched biaxially by a stretch ratio of &#x003BB;<sub><italic>p</italic></sub>. Then this pre-stretched membrane is bounded to a rigid ring with the radius <italic>b</italic> and an inner disk of radius <italic>a</italic>, as shown in Figure <xref ref-type="fig" rid="F4">4A</xref>. Compliant electrodes are coated on both sides of the membrane. A force <italic>F</italic> moves the inner disk out-of-plane relative to the outer ring by a distance <italic>h</italic> with an actuation voltage V applied across the electrodes. The membrane is deformed into a conical shape, as illustrated in Figure <xref ref-type="fig" rid="F4">4B</xref>. Due to this out-of-plane deformation, a particle on the membrane at position <italic>R</italic> (as marked red Figure <xref ref-type="fig" rid="F4">4A</xref>) now moves to the position of [<italic>r</italic>(<italic>R</italic>), <italic>z</italic>(<italic>R</italic>)] in Figure <xref ref-type="fig" rid="F4">4B</xref>, where <italic>r</italic> is the current radius and <italic>z</italic> is the distance to the undeformed plane. The coordinates of (<italic>r, z</italic>) for <italic>R</italic> &#x0003D; [<italic>a, b</italic>] describe the geometry of the conical DEA, and a summarize of this model is presented as follows [after (He et al., <xref ref-type="bibr" rid="B16">2010</xref>) and adopted from Cao and Conn (<xref ref-type="bibr" rid="B6">2018</xref>)].</p>
<fig id="F4" position="float">
<label>Figure 4</label>
<caption><p>The conceptual illustrations of a single cone DEA. <bold>(A)</bold> Pre-stretched membrane is bonded to a rigid ring and a central disk. <bold>(B)</bold> DEA membrane deforms out-of-plane with the existence of a force <italic>F</italic> and a voltage <italic>V</italic>.</p></caption>
<graphic xlink:href="frobt-05-00137-g0004.tif"/>
</fig>
<p>The state of the point <italic>R</italic> (<italic>r</italic> ,<italic>z</italic> ,&#x003B8;) can be expressed as</p>
<disp-formula id="E1"><label>(1)</label><mml:math id="M1"><mml:mtable class="eqnarray" columnalign="right center left"><mml:mtr><mml:mtd><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:mi>r</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>R</mml:mi></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:msubsup><mml:mrow><mml:mo>&#x003BB;</mml:mo></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>&#x02032;</mml:mi></mml:mrow></mml:msubsup><mml:mo class="qopname">cos</mml:mo><mml:mi>&#x003B8;</mml:mi><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<disp-formula id="E2"><label>(2)</label><mml:math id="M2"><mml:mtable class="eqnarray" columnalign="right center left"><mml:mtr><mml:mtd><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:mi>z</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>R</mml:mi></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:msubsup><mml:mrow><mml:mo>&#x003BB;</mml:mo></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>&#x02032;</mml:mi></mml:mrow></mml:msubsup><mml:mo class="qopname">sin</mml:mo><mml:mi>&#x003B8;</mml:mi><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<disp-formula id="E3"><label>(3)</label><mml:math id="M3"><mml:mtable class="eqnarray" columnalign="right center left"><mml:mtr><mml:mtd><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>R</mml:mi></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mrow><mml:mi>s</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mi>R</mml:mi><mml:msub><mml:mrow><mml:mi>s</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mfrac><mml:mo class="qopname">sin</mml:mo><mml:mi>&#x003B8;</mml:mi><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where <italic>s</italic><sub>1</sub> and <italic>s</italic><sub>2</sub> are the nominal longitudinal and radial stress, respectively and <inline-formula><mml:math id="M4"><mml:msubsup><mml:mrow><mml:mo>&#x003BB;</mml:mo></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>&#x02032;</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula> is the radial stretch due to the out-of-plane deformation.</p>
<p>In quasi-static state, the external force <italic>F</italic> and the reaction force exerted by the deformed membrane are balanced, and this relationship gives</p>
<disp-formula id="E4"><label>(4)</label><mml:math id="M5"><mml:mtable class="eqnarray" columnalign="right center left"><mml:mtr><mml:mtd><mml:mn>2</mml:mn><mml:mi>&#x003C0;</mml:mi><mml:mfrac><mml:mrow><mml:mi>T</mml:mi></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mo>&#x003BB;</mml:mo></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:msub><mml:mrow><mml:mo>&#x003BB;</mml:mo></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mfrac><mml:mi>R</mml:mi><mml:msub><mml:mrow><mml:mi>&#x003C3;</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo class="qopname">sin</mml:mo><mml:mi>&#x003B8;</mml:mi><mml:mo>=</mml:mo><mml:mi>F</mml:mi><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where &#x003C3;<sub>1</sub> &#x0003D; &#x003BB;<sub>1</sub><italic>s</italic><sub>1</sub>.</p>
<p>The three differential equations (1&#x02013;3) and an algebraic equation (4) govern the state of the DEA and are numerically solved in MATLAB using shooting method and &#x0201C;<italic>ode15</italic>&#x0201D; function.</p>
<p>This model can be used in two ways such that if the force <italic>F</italic> and voltage <italic>V</italic> are given, this model can predict the out-of-plane deformation <italic>h</italic> of the single cone DEA, or if the out-of-plane deformation <italic>h</italic> and the voltage <italic>V</italic> are known, this model can predict the force <italic>F</italic>. This model can be extended to a double cone DEA by the constrain: the sum of the out-of-plane deformation of the two membranes <italic>h</italic><sub>1</sub> and <italic>h</italic><sub>2</sub> is a constant value <italic>L</italic>. The physical meaning of this constrain is that as the out-of-plane deformation for a double cone DEA is caused by the spacer (as shown in Figure <xref ref-type="fig" rid="F2">2</xref>), the total deformation of the two membranes should always be the length of the spacer <italic>L</italic>. For a given deformation of the double cone DEA, the force exerted by the DEA can be solved using this model as -<italic>F</italic><sub>1</sub>&#x0002B;<italic>F</italic><sub>2</sub> where <italic>F</italic><sub>1</sub> and <italic>F</italic><sub>2</sub> the forces exerted by top and bottom membranes, respectively.</p>
</sec>
<sec>
<title>Quasi-Static Model Validation</title>
<p>The stress&#x02013;strain relationship of the silicone elastomer is described using the Ogden hyperelastic model (Ogden, <xref ref-type="bibr" rid="B25">1972</xref>) in this work, as written in equation (5). In order to obtain the parameters in the Ogden model, a pure-shear pull test was conducted on a sample of the elastomer. Figure <xref ref-type="fig" rid="F5">5</xref> shows the experimental result and the fitted Ogden model. The model agrees very well with the measured data with the parameters of <inline-formula><mml:math id="M6"><mml:msub><mml:mrow><mml:mi>&#x003BC;</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>1</mml:mn><mml:mo>.</mml:mo><mml:mn>287</mml:mn><mml:mo>&#x000D7;</mml:mo><mml:mn>1</mml:mn><mml:msup><mml:mrow><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>5</mml:mn></mml:mrow></mml:msup><mml:mtext>&#x000A0;</mml:mtext><mml:mi>P</mml:mi><mml:mi>a</mml:mi></mml:math></inline-formula>, <inline-formula><mml:math id="M7"><mml:msub><mml:mrow><mml:mi>&#x003BC;</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>2</mml:mn><mml:mo>.</mml:mo><mml:mn>6</mml:mn><mml:mo>&#x000D7;</mml:mo><mml:mn>1</mml:mn><mml:msup><mml:mrow><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>4</mml:mn></mml:mrow></mml:msup><mml:mtext>&#x000A0;</mml:mtext><mml:mi>P</mml:mi><mml:mi>a</mml:mi></mml:math></inline-formula>, &#x003B1;<sub>1</sub> &#x0003D; 2 and &#x003B1;<sub>2</sub> &#x0003D; 4.</p>
<disp-formula id="E5"><label>(5)</label><mml:math id="M8"><mml:mtable columnalign='left'><mml:mtr><mml:mtd><mml:msub><mml:mi>&#x003C3;</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>&#x003BC;</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:msub><mml:mo>&#x003BB;</mml:mo><mml:mn>1</mml:mn></mml:msub></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi>&#x003B1;</mml:mi><mml:mn>1</mml:mn></mml:mrow></mml:msup><mml:mo>&#x02212;</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:msub><mml:mo>&#x003BB;</mml:mo><mml:mn>1</mml:mn></mml:msub></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi>&#x003B1;</mml:mi><mml:mn>1</mml:mn></mml:mrow></mml:msup><mml:msup><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:msub><mml:mo>&#x003BB;</mml:mo><mml:mn>2</mml:mn></mml:msub></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi>&#x003B1;</mml:mi><mml:mn>1</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:mfrac></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mtext>&#x02009;&#x02009;&#x02009;&#x02009;&#x02009;&#x02009;&#x02009;&#x02009;&#x02009;&#x02009;&#x02009;&#x02009;&#x02009;</mml:mtext><mml:mo>+</mml:mo><mml:msub><mml:mi>&#x003BC;</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:msub><mml:mo>&#x003BB;</mml:mo><mml:mn>1</mml:mn></mml:msub></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi>&#x003B1;</mml:mi><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mo>&#x02212;</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:msub><mml:mo>&#x003BB;</mml:mo><mml:mn>1</mml:mn></mml:msub></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi>&#x003B1;</mml:mi><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:msup><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:msub><mml:mo>&#x003BB;</mml:mo><mml:mn>2</mml:mn></mml:msub></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi>&#x003B1;</mml:mi><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:mfrac></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>To verify the quasi-static model developed, three double cone DEAs with the spacer heights <italic>L</italic> &#x0003D; 20, 25 and 30 mm were used. The pre-stretch ratio is 1.0 &#x000D7; 1.0, <italic>a</italic> &#x0003D; 4 mm and <italic>b</italic> &#x0003D; 20 mm. The DEA was fixed to the testing rig and was deformed from its equilibrium position by a linear rail at a speed of 0.01 mm/s to eliminate the effect of viscosity. The experimental setup is illustrated in Figure <xref ref-type="fig" rid="F6">6</xref> and the experimental results are shown in Figure <xref ref-type="fig" rid="F7">7</xref>. The model agrees very well with the experimental results. Despite the non-linearity of the elastomer, the double cone DEA demonstrates an approximated linear passive force-displacement relationship, which means the stiffness of the actuator is approximately constant.</p>
<fig id="F5" position="float">
<label>Figure 5</label>
<caption><p>Pure shear test of the silicone elastomer and the Ogden model fit.</p></caption>
<graphic xlink:href="frobt-05-00137-g0005.tif"/>
</fig>
<fig id="F6" position="float">
<label>Figure 6</label>
<caption><p>Double cone DEA passive quasi-static test setup.</p></caption>
<graphic xlink:href="frobt-05-00137-g0006.tif"/>
</fig>
<fig id="F7" position="float">
<label>Figure 7</label>
<caption><p>Passive quasi-static experimental results and model predictions. <italic>L</italic> &#x0003D; 20, 25, 30 mm.</p></caption>
<graphic xlink:href="frobt-05-00137-g0007.tif"/>
</fig>
</sec>
<sec>
<title>Double Cone DEA Resonance Characterization</title>
<p>In the last section, we have shown that a double cone DEA has an approximated linear force-displacement relationship, and can be treated as a linear spring. In dynamic actuation, this DEA can be simplified as a linear spring with an inertial load (mass of the spacer and any additional inertial load connected to the spacer). This configuration is very similar to that of a mechanical undamped oscillator, hence for dynamic characterization, in this work we simplify the DEA-mass oscillator as a linear undamped oscillator. As a result, the natural frequency <italic>f</italic> <sub>0</sub> of the DEA-mass oscillator can be expressed as <inline-formula><mml:math id="M9"><mml:msub><mml:mrow><mml:mi>f</mml:mi></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msqrt><mml:mrow><mml:mfrac><mml:mrow><mml:mi>K</mml:mi></mml:mrow><mml:mrow><mml:mi>M</mml:mi></mml:mrow></mml:mfrac></mml:mrow></mml:msqrt></mml:mrow><mml:mrow><mml:mn>2</mml:mn><mml:mi>&#x003C0;</mml:mi></mml:mrow></mml:mfrac></mml:math></inline-formula>, where <italic>K</italic> (N/m) is the approximated stiffness of the DEA and <italic>M</italic> (kg) is the mass of the spacer and any payload in the undamped oscillator (assuming the mass of the DEA membranes is negligible). The damping ratio of the DEA is neglected in the natural frequency characterization as the Parker silicone elastomer shows very low viscoelasticity thus we assume effect of damping from the elastomer can be negligible and the good agreement between the model prediction and experimental results suggests that this assumption holds.</p>
<p>To verify the simplified dynamic model, five DEA-mass oscillator specimens were fabricated with different stiffness and mass values. To adjust the stiffness <italic>K</italic> of these specimens, different spacer length <italic>L</italic> were used, as the spacer becomes longer, the two DEA membranes deform out-of-plane further and causing an increase in the stiffness <italic>K</italic>. The stiffness was also adjusted by adding a second layer of elastomer on the DEA. By adding an additional layer of membrane with the same pre-stretch to both top and bottom cones, the stiffness is doubled. The payload was adjusted by adding nuts to the spacer. Table <xref ref-type="table" rid="T1">1</xref> lists the detailed values of the stiffness and mass of the five specimens, as well as their model predicted and measured resonant frequencies. The predicted resonant frequencies show a very good agreement with the experimental results with an average relative error of 1.98%. The promising results confirm that the linear spring simplification and the negligible damping ratio assumption hold. By comparing the five specimens, it can be noticed that the DEAs with heavier payload and lower stiffness have larger amplitudes at their resonance.</p>
<table-wrap position="float" id="T1">
<label>Table 1</label>
<caption><p>Comparison of the measured and model predicted resonant frequencies for the five DEA-mass oscillators.</p></caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th/>
<th valign="top" align="center"><bold>DEA 1</bold></th>
<th valign="top" align="center"><bold>DEA 2</bold></th>
<th valign="top" align="center"><bold>DEA 3</bold></th>
<th valign="top" align="center"><bold>DEA 4</bold></th>
<th valign="top" align="center"><bold>DEA 5</bold></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left"><italic>K</italic> (N/m)</td>
<td valign="top" align="center">63.5</td>
<td valign="top" align="center">63.5</td>
<td valign="top" align="center">63.5</td>
<td valign="top" align="center">41.4</td>
<td valign="top" align="center">127</td>
</tr>
<tr>
<td valign="top" align="left"><italic>M</italic> (g)</td>
<td valign="top" align="center">0.603</td>
<td valign="top" align="center">1.26</td>
<td valign="top" align="center">2.538</td>
<td valign="top" align="center">1.26</td>
<td valign="top" align="center">1.26</td>
</tr>
<tr>
<td valign="top" align="left"><italic>d</italic> (mm)</td>
<td valign="top" align="center">7.2</td>
<td valign="top" align="center">8.8</td>
<td valign="top" align="center">9.35</td>
<td valign="top" align="center">10.85</td>
<td valign="top" align="center">8.12</td>
</tr>
<tr>
<td valign="top" align="left"><italic>f<sub>0</sub></italic> predicted (Hz)</td>
<td valign="top" align="center">51.8</td>
<td valign="top" align="center">35.7</td>
<td valign="top" align="center">25.2</td>
<td valign="top" align="center">28.8</td>
<td valign="top" align="center">50.5</td>
</tr>
<tr>
<td valign="top" align="left"><italic>f<sub>0</sub></italic> measured (Hz)</td>
<td valign="top" align="center">52</td>
<td valign="top" align="center">37</td>
<td valign="top" align="center">26</td>
<td valign="top" align="center">28</td>
<td valign="top" align="center">51</td>
</tr>
<tr>
<td valign="top" align="left">Relative error of <italic>f<sub>0</sub></italic></td>
<td valign="top" align="center">0.39%</td>
<td valign="top" align="center">3.64%</td>
<td valign="top" align="center">1.98%</td>
<td valign="top" align="center">2.9%</td>
<td valign="top" align="center">0.98%</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
</sec>
</sec>
<sec id="s3">
<title>Power Output of the Double Cone DEA</title>
<sec>
<title>Mechanical Power Output Against Excitation Frequency</title>
<p>High power output from the actuator is crucial for flapping wing MAV designs. The large amplitude of the silicone double cone DEA at its resonance could lead to a high power output. In this section, we investigate the mechanical power output of the silicone double cone DEA against its excitation frequency. Stacking multiple layers of DEA membranes have been shown to be able to amplify the force output of a DEA (Carpi et al., <xref ref-type="bibr" rid="B10">2007</xref>; Kovacs et al., <xref ref-type="bibr" rid="B21">2009</xref>). So in addition to exploiting resonance, increasing the number of layers of the DEA membrane may also increase the power output of the DEA. In this section, we also analyze the effect of layer numbers to the power output of the actuator. As increasing the number of layers will also increase the stiffness and the resonant frequency, in order to the keep the resonant frequency constant and ensures a consistent comparison, the payload also increases together with the layer number (recall that <inline-formula><mml:math id="M10"><mml:msub><mml:mrow><mml:mi>f</mml:mi></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msqrt><mml:mrow><mml:mfrac><mml:mrow><mml:mi>K</mml:mi></mml:mrow><mml:mrow><mml:mi>M</mml:mi></mml:mrow></mml:mfrac></mml:mrow></mml:msqrt></mml:mrow><mml:mrow><mml:mn>2</mml:mn><mml:mi>&#x003C0;</mml:mi></mml:mrow></mml:mfrac></mml:math></inline-formula>, as <italic>K</italic> increases, the mass <italic>M</italic> should increase with the same ratio to ensure a constant <italic>f</italic> <sub>0</sub>). As the DEA-mass system oscillates, the instant mechanical power output is written as <inline-formula><mml:math id="M11"><mml:msub><mml:mrow><mml:mi>P</mml:mi></mml:mrow><mml:mrow><mml:mi>m</mml:mi><mml:mi>e</mml:mi><mml:mi>c</mml:mi><mml:mi>h</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mi>M</mml:mi><mml:mfrac><mml:mrow><mml:msup><mml:mrow><mml:mi>d</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:msup><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:mfrac><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mfrac></mml:math></inline-formula>, where <italic>x</italic> is the displacement of the payload and <italic>t</italic> is time. The average mechanical power output is simply <inline-formula><mml:math id="M12"><mml:msub><mml:mrow><mml:mi>P</mml:mi></mml:mrow><mml:mrow><mml:mi>m</mml:mi><mml:mi>e</mml:mi><mml:mi>c</mml:mi><mml:mi>h</mml:mi><mml:mstyle class="text"><mml:mtext>_</mml:mtext></mml:mstyle><mml:mi>a</mml:mi><mml:mi>v</mml:mi><mml:mi>g</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>T</mml:mi></mml:mrow></mml:mfrac><mml:msubsup><mml:mrow><mml:mo>&#x0222B;</mml:mo></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mi>T</mml:mi></mml:mrow></mml:msubsup><mml:msub><mml:mrow><mml:mi>P</mml:mi></mml:mrow><mml:mrow><mml:mi>m</mml:mi><mml:mi>e</mml:mi><mml:mi>c</mml:mi><mml:mi>h</mml:mi></mml:mrow></mml:msub><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:math></inline-formula>, where <italic>T</italic> is the period of one actuation cycle.</p>
</sec>
<sec>
<title>Experiment and Results</title>
<p>Four silicone double cone DEA specimens with different numbers of layers of membranes were fabricated, each specimen has the same frame radius <italic>b</italic> of 13 mm and the spacer height <italic>L</italic> of 13 mm. As discussed in previous section, to keep the resonant frequency constant, the payload has to increase with the layer numbers. Hence for the four specimens with 2, 4, 6 and 8 layers of silicone membranes, the payload for each specimen is 1, 2, 3 and 4 grams, respectively, which results in a constant resonant frequency of 45.8 Hz predicted by the dynamic model. Each layer includes Parker silicone elastomer, Mylar ring frame, carbon grease electrode and copper tape connection, the total mass is 0.12 g. Same frequency sweep from 1 to 100 Hz as described in Section Double Cone Linear DEA Design was conducted, and the experimental data was analyzed in MATLAB.</p>
<p>Figure <xref ref-type="fig" rid="F8">8</xref> shows the measured average mechanical power output of each specimen against the excitation frequency. It can be noticed that at low frequencies (&#x0003C;35 Hz), the DEA outputs very low mechanical power. Then the average mechanical power output increases rapidly as the excitation frequency increases and peaks at the resonant frequency. The power output then drops when the excitation frequency is greater than the resonant frequency. Figure <xref ref-type="fig" rid="F9">9</xref> shows the time series of the displacement of the payload, the kinetic energy of the payload and the elastic energy stored in the elastomer at the resonance of the 2 layers DEA specimen. It is clear that at resonance, a close to 90&#x000B0; phase shift can be observed between the elastic energy stored in elastomer and the kinetic energy of the payload, which means that when each half stroke starts, the elastic energy of the elastomer is converted into the kinetic energy of the payload, this gives the mass a higher acceleration and faster velocity, hence a larger mechanical power output. When the mass approaches the end of a half stroke, its kinetic energy is converted into the elastic potential energy in the membranes, which is then released in the next half cycle. Apart from increasing the mechanical power output of the DEA, elastic energy recovery can also increase the efficiency, a detailed elastic energy recovery study on double cone DEAs can be found in Cao and Conn (<xref ref-type="bibr" rid="B7">2017</xref>).</p>
<fig id="F8" position="float">
<label>Figure 8</label>
<caption><p>Average mechanical power against excitation frequency for 2, 4, 6, and 8 layers of double cone silicone DEAs. Resonant frequency is fixed at about 46 Hz.</p></caption>
<graphic xlink:href="frobt-05-00137-g0008.tif"/>
</fig>
<fig id="F9" position="float">
<label>Figure 9</label>
<caption><p>Time series of the displacement, kinetic energy of the mass and elastic potential energy in the DEA membrane at its resonance. Kinetic energy and elastic energy are 90&#x000B0; out-of-phase.</p></caption>
<graphic xlink:href="frobt-05-00137-g0009.tif"/>
</fig>
<p>From Figure <xref ref-type="fig" rid="F8">8</xref>, by comparing the mechanical power output of the four specimens, increasing the number of layers can increase the mechanical power output of a DEA. It should be noted that increasing the layer numbers also increases the capacitance of the DEA which increases the payload on the high voltage amplifier hence the maximum number of layers can be added is restricted by the high voltage amplifier (current amplifier has a maximum power output of 1 W). Table <xref ref-type="table" rid="T2">2</xref> compares the average mechanical power output of the DEAs of 2, 4, 6, and 8 layers of membranes and the mass-specific powers at their resonant frequency, the mechanical power scaled up with the layer number as expected. However, maybe due to the increase in capacitance and the connection issue between DEA electrodes and high voltage channels, <italic>n</italic> layer specimens (<italic>n</italic> &#x0003E; 2) did not generate as much mechanical powers as <italic>n</italic>/(<italic>n</italic> &#x0003D; 2) times than the 2-layer one. The mass-specific power, which is the mechanical power output divided by the total mass of the actuator, is about 100 mW/g at resonance. By comparison, the mass-specific power of insect flight muscle is estimated to be between 80 to 83 mW/g (Lehmann and Dickinson, <xref ref-type="bibr" rid="B24">1997</xref>; Tu and Daniel, <xref ref-type="bibr" rid="B31">2004</xref>). Piezoelectric actuators can have a mass-specific power as high as 400 mW/g (Steltz and Fearing, <xref ref-type="bibr" rid="B29">2007</xref>), however, the actuation strokes is very low (&#x0003C;1%) and only this high mass-specific power is only achieved at high frequencies (say &#x0003E; 100 Hz) (Huber et al., <xref ref-type="bibr" rid="B19">1997</xref>). On the contrary, the double cone DEA presented in this work has a very large stroke at its resonance (e.g., 99% stroke relative to the height of the DEA in Figure <xref ref-type="fig" rid="F9">9</xref>), it also has a mass-specific power close to insect flight muscle and an actuation bandwidth close to medium to large insect species, which makes it ideal for insect-inspired flapping wing MAV applications.</p>
<table-wrap position="float" id="T2">
<label>Table 2</label>
<caption><p>Average mechanical power at resonance and the corresponding mass-specific power of 2, 4, 6, and 8 layers of double cone silicone DEA.</p></caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th valign="top" align="left"><bold><italic>Number of layers n</italic></bold></th>
<th valign="top" align="center"><bold>2</bold></th>
<th valign="top" align="center"><bold>4</bold></th>
<th valign="top" align="center"><bold>6</bold></th>
<th valign="top" align="center"><bold>8</bold></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left"><italic>P<sub><italic>max</italic></sub></italic> (mW)</td>
<td valign="top" align="center">64.49</td>
<td valign="top" align="center">88.38</td>
<td valign="top" align="center">186.6</td>
<td valign="top" align="center">209.1</td>
</tr>
<tr>
<td valign="top" align="left"><italic>P<sub><italic>max</italic></sub></italic>/<italic>P<sub><italic>max</italic></sub></italic>(<italic>n</italic> &#x0003D; 2)</td>
<td valign="top" align="center">100%</td>
<td valign="top" align="center">137%</td>
<td valign="top" align="center">289%</td>
<td valign="top" align="center">324%</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Mass-specific power</italic> (mW/g)</td>
<td valign="top" align="center">134.35</td>
<td valign="top" align="center">92.06</td>
<td valign="top" align="center">129.58</td>
<td valign="top" align="center">108.91</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
</sec>
<sec id="s4">
<title>Flapping Wing MAV Design</title>
<sec>
<title>Design Overview</title>
<p>In this section we present a double cone DEA driven flapping wing MAV design with no additional elastic element added since the elasticity is embedded in the actuator itself. A schematic illustration of this design is shown in Figure <xref ref-type="fig" rid="F10">10</xref> with all the parameters listed in Table <xref ref-type="table" rid="T3">3</xref>. In the last section we have shown that increasing the number of layers of DEA membranes can increase its mechanical power output, as a result, in this MAV design, eight layers of silicone elastomer were used to ensure a good power output from the actuator. The flapping mechanism includes a slider crank mechanism which converts the linear motion of the DEA into a reciprocal flapping motion and a double rocker mechanism to amplify the flapping stroke. Singularities are avoided by including two mechanical stops the slider. It should be pointed out that since the aim of this design is to demonstrate the concept of a DE resonator driven flapping wing MAV in benchtop tests and hence no optimization has been done to the size and weight of this mechanism to achieve free flight.</p>
<fig id="F10" position="float">
<label>Figure 10</label>
<caption><p>Schematic diagram of the DE resonator driven flapping wing MAV design (not to scale).</p></caption>
<graphic xlink:href="frobt-05-00137-g0010.tif"/>
</fig>
<table-wrap position="float" id="T3">
<label>Table 3</label>
<caption><p>Design parameters and values for the flapping wing MAV design.</p></caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th valign="top" align="left"><bold>Parameter</bold></th>
<th valign="top" align="center"><bold><italic>L</italic></bold></th>
<th valign="top" align="center"><bold><italic>b</italic></bold></th>
<th valign="top" align="center"><bold><italic>L<sub><bold>1</bold></sub></italic></bold></th>
<th valign="top" align="center"><bold><italic>L<sub><bold>2</bold></sub></italic></bold></th>
<th valign="top" align="center"><bold><italic>L<sub><bold>3</bold></sub></italic></bold></th>
<th valign="top" align="center"><bold><italic>L<sub><bold>4</bold></sub></italic></bold></th>
<th valign="top" align="center"><bold><italic>L<sub><bold>5</bold></sub></italic></bold></th>
<th valign="top" align="center"><bold><italic>L<sub><bold><italic>wing</italic></bold></sub></italic></bold></th>
<th valign="top" align="center"><bold><italic>m<sub><bold><italic>wing</italic></bold></sub></italic></bold></th>
<th valign="top" align="center"><bold><italic>AR</italic></bold></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">Value</td>
<td valign="top" align="center">13 mm</td>
<td valign="top" align="center">13 mm</td>
<td valign="top" align="center">8 mm</td>
<td valign="top" align="center">4 mm</td>
<td valign="top" align="center">8.5 mm</td>
<td valign="top" align="center">18 mm</td>
<td valign="top" align="center">5 mm</td>
<td valign="top" align="center">40 mm</td>
<td valign="top" align="center">0.14 g</td>
<td valign="top" align="center">6</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec>
<title>Flapping Tests and Results</title>
<p>The experimental setup is similar to that of the DEA-mass dynamic test. The flapping wing MAV was fixed to the testing rig and two square waves with an amplitude of 1,920 V and 180&#x000B0; phase shift were used to drive the antagonistic DEA membranes. An actuation frequency from 1 to 50 Hz with the steps of 5 Hz (step was reduced to 2 Hz near amplitude peak) was tested. The experiments were recorded by a camera (GoPro Session 5) at 100 frames-per-second. The stroke of the DEA was measured by the laser displacement sensor and the flapping stroke was estimated based on the DEA stroke and the DEA to flapping transmission ratio of this mechanism and verified against video footage (example video in <xref ref-type="supplementary-material" rid="SM1">Supplementary Materials</xref>).</p>
<p>In the first test, wings were not included to the mechanism, hence the DEA actuated the mechanism only. As is shown as the blue curve in Figure <xref ref-type="fig" rid="F11">11</xref>, the rotational stroke of this mechanism increases with the increasing actuation frequency first and peaks at a frequency of 34 Hz. The stroke then begins to reduce as the frequency increases further above 34 Hz, this is due to the inertia of the linkages and frictional loss.</p>
<fig id="F11" position="float">
<label>Figure 11</label>
<caption><p>Flapping strokes of the MAV against excitation frequency with no wing connected, wing frames only and full wings attached.</p></caption>
<graphic xlink:href="frobt-05-00137-g0011.tif"/>
</fig>
<p>In the second test, only the wing frames were connected to the flapping mechanism, which introduced the inertia of the wings to this system. No aerodynamic force was included in this set of experiment. As is shown in red dash curve in Figure <xref ref-type="fig" rid="F11">11</xref>, the inertia of the wings lowers the resonant frequency from 34 to 22 Hz, the flapping stroke is also reduced from 96 to 83&#x000B0;. Above the resonant frequency, the flapping stroke drops sharply as the frequency increases.</p>
<p>In the last test, wing membranes were included, which introduces aerodynamic drag to the system. The aerodynamic force reduces the peak flapping stroke to 63&#x000B0; at 18 Hz. Figure <xref ref-type="fig" rid="F12">12</xref> shows the flapping wing MAV with full wings attached at its top and bottom stop at 18 Hz. The measured displacement of the DEA and the estimated flapping angle of the mechanism at 18 Hz are shown in Figure <xref ref-type="fig" rid="F13">13A</xref>, as can be seen, despite that a square wave was applied to the DEA, the displacement is approximately sinusoidal.</p>
<fig id="F12" position="float">
<label>Figure 12</label>
<caption><p>DEA driven flapping wing MAV with full wings attached at 18 Hz, showing upper <bold>(a)</bold> and lower <bold>(b)</bold> stroke reversal.</p></caption>
<graphic xlink:href="frobt-05-00137-g0012.tif"/>
</fig>
<fig id="F13" position="float">
<label>Figure 13</label>
<caption><p><bold>(A)</bold> The displacement of the DEA (measured) and the flapping angle (estimated) at 18 Hz. <bold>(B)</bold> the speed ratio of the four-bar mechanism.</p></caption>
<graphic xlink:href="frobt-05-00137-g0013.tif"/>
</fig>
<p>An empirical equation that describes the mass that can be supported during hovering is given by Ellington (<xref ref-type="bibr" rid="B15">1999</xref>) in Equation 6, based on this equation, the current prototype can generate a lift of about 133 mg.</p>
<disp-formula id="E7"><label>(6)</label><mml:math id="M13"><mml:mtable class="eqnarray" columnalign="right center left"><mml:mtr><mml:mtd><mml:mi>m</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn><mml:mo>.</mml:mo><mml:mn>387</mml:mn><mml:mfrac><mml:mrow><mml:msup><mml:mrow><mml:mo class="MathClass-ord">&#x02205;</mml:mo></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:msup><mml:mrow><mml:mi>n</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:msup><mml:mrow><mml:mi>R</mml:mi></mml:mrow><mml:mrow><mml:mn>4</mml:mn></mml:mrow></mml:msup><mml:msub><mml:mrow><mml:mi>C</mml:mi></mml:mrow><mml:mrow><mml:mi>L</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mi>A</mml:mi><mml:mi>R</mml:mi></mml:mrow></mml:mfrac><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where <italic>m</italic> is mass (kg), &#x02205; is flapping stroke (rad), <italic>n</italic> is frequency (Hz), <italic>R</italic> is the wing length (m), <italic>C</italic><sub><italic>L</italic></sub> is the lift coefficient (<italic>C</italic><sub><italic>L</italic></sub> &#x0003D; 2 to 3 for insect hovering) and <italic>AR</italic> is the aspect ratio of the wing.</p>
</sec>
<sec>
<title>Discussion</title>
<p>The MAV test shows promising initial results for DEA driven flapping wing MAVs, however, it also suggests that to achieve hovering, significant improvements have to be made in both DEA designs and flapping mechanism optimizations. As the empirical aerodynamic equation (equation 6) suggests, to generate sufficient lift for a 10 g flapping wing MAV to achieve hovering flight, and by assuming the stroke amplitude &#x003C6; &#x0003D; 120&#x000B0; and <italic>AR</italic> &#x0003D; 7 [commonly seen in flying insects ], the MAV has to flap its wings at 100 Hz for a wing length of 40 mm or at 20 Hz for a wing length of 80 mm. For DEAs with multiple layers stacked together, due to its large RC time constant (&#x003C4; &#x0003D; <italic>RC</italic>, where <italic>R</italic> is the resistance and <italic>C</italic> is the capacitance), it is favorable to actuate the DEAs at relatively low frequencies which will allow sufficient charging and discharging periods. As a result, a DEA actuated MAV with a pair long wings and a low resonant frequency possibly are preferred. Based on this discussion, further development can be made several key areas. First is to improve the electrode conductivity using alternative materials, such as carbon nanotubes [such as in Yuan et al. (<xref ref-type="bibr" rid="B36">2008</xref>)]. The reduction in the surface resistant of the DEA will improve its RC time constant, which means the actuator can respond faster to the input signal thus increasing the mechanical power output. Second is to optimize the coupling between the DEA and flapping mechanism so that long wings are used, and the resonance frequency is tuned at a relatively low value (say 15&#x02013;25 Hz). Last but not least, the flapping mechanism can be optimized in terms of speed ratios by adjusting the link ratios. In Figure <xref ref-type="fig" rid="F13">13B</xref>, the speed ratio of the flapping wing mechanism against the displacement of the DEA is plotted. It can be noted that the two reversals (i.e., upper and lower end stops) the speed ratios are the lowest, which suggest that when the DEA is trying the accelerate the wings, it has the poorest mechanical advantage. The low speed ratios at two reversals can be a significant limitation of this mechanism design. In the future work, the mechanism design should ensure that the speed ratios at the two reversals are the highest in one complete cycle. Compliant transmission mechanisms, such as active hinges as adopted in Wood (<xref ref-type="bibr" rid="B35">2008</xref>) can be used in the future mechanism designs to eliminate the effect of frictions and the elastic energy stored in the hinges can help accelerate the wings in the beginning of each stroke.</p>
</sec>
</sec>
<sec sec-type="conclusions" id="s5">
<title>Conclusion</title>
<p>In this work, a double cone DEA driven flapping wing MAV was presented. First, two predominant dielectric elastomer materials polyacrylate tape VHB (3M) and silicone elastomers were compared. Frequency sweep test showed that the DEA made with silicone elastomer demonstrates a high peak in amplitude at its resonance while the DEA made with VHB has no peak at resonance due to its high viscoelasticity. This result suggested that the low viscoelasticity of silicone elastomers makes them a better candidate for flapping wing MAV actuations. An analytical model and a simplified mass-spring model were used to characterize the quasi-static and dynamic performance of the silicone double cone DEA, respectively. The analytical model adopted in this work was able to predict the quasi-static performance of the actuator accurately and both the model and experimental results showed that double cone DEA has an approximated linear force-displacement relationship, which suggested that the DEA can be simplified as a linear spring. This finding then allowed us to use a classic spring-mass oscillator model to estimate the natural frequency of the DEA-mass oscillator in an extremely simple way yet high accuracy. The mechanical power output of the DEA was analyzed against the excitation frequency and the number of layers of membranes stacked to the DEA. The mechanical power output was shown to have a peak at the predicted resonant frequency of the DEA-mass system and the power scaled up with the number of layers of DEA membranes added to the DEA. A 209.1 mW average mechanical power output and an equivalent of 108.9 W/kg mass-specific mechanical power density from an 8-layer DEA was demonstrated. Subsequently, by using the same 8-layer DEA to drive the proposed insect-inspired flapping wing MAV, a peak flapping stroke of 63&#x000B0; at 18 Hz was demonstrated. While this flapping performance requires further optimization toward a hovering flight, this prototype far outperforms previous DEA-based flapping mechanisms without exploiting resonant actuations [e.g., a flapping amplitude of 15&#x000B0; in Henke et al. (<xref ref-type="bibr" rid="B17">2017</xref>) and 10&#x000B0; at a frequency of 5 Hz in Lau et al. (<xref ref-type="bibr" rid="B23">2014</xref>)]. For future work, we will continue working on optimizing the actuator and maximizing its mass-specific power, also, new flapping mechanism designs and fabrication techniques will be attempted to optimize the power transmission and efficiency.</p>
</sec>
<sec id="s6">
<title>Data Availability Statement</title>
<p>Data are available at the University of Bristol data repository, data.bris, at <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.5523/bris.1pcnxotj9hpfr2wnjnhf1x5xst">https://doi.org/10.5523/bris.1pcnxotj9hpfr2wnjnhf1x5xst</ext-link>.</p>
</sec>
<sec id="s7">
<title>Author Contributions</title>
<p>CC: data curation. CC, SB, and AC: formal analysis. CC, SB, and AC: investigation. CC, SB, and AC: methodology. SB and AC: supervision. CC: writing&#x02013;original draft. AC: writing&#x02013;review and editing.</p>
<sec>
<title>Conflict of Interest Statement</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
</sec>
</body>
<back>
<ack><p>This work was supported by the EPSRC Centre for doctoral training in Future Autonomous and Robotic Systems (FARSCOPE) at the Bristol Robotics Laboratory where CC is a PhD student. AC was supported by EPSRC grant EP/P025846/1.</p>
</ack>
<sec sec-type="supplementary-material" id="s8">
<title>Supplementary Material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/frobt.2018.00137/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/frobt.2018.00137/full#supplementary-material</ext-link></p>
<supplementary-material xlink:href="Video_1.MP4" id="SM1" mimetype="video/mp4" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Baek</surname> <given-names>S. S.</given-names></name> <name><surname>Ma</surname> <given-names>K. Y.</given-names></name> <name><surname>Fearing</surname> <given-names>R. S.</given-names></name></person-group> (<year>2009</year>). <article-title>Efficient resonant drive of flapping-wing robots</article-title>, in <source>Intelligent Robots and Systems, IROS 2009. IEEE/RSJ International Conference</source> (<publisher-loc>St. Louis, MO</publisher-loc>).</citation></ref>
<ref id="B2">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bolsman</surname> <given-names>C.</given-names></name> <name><surname>Goosen</surname> <given-names>J.</given-names></name> <name><surname>Van Keulen</surname> <given-names>F.</given-names></name></person-group> (<year>2009</year>). <article-title>Design overview of a resonant wing actuation mechanism for application in flapping wing MAVs</article-title>. <source>Int. J. Micro. Air Vehicles</source> <volume>1</volume>, <fpage>263</fpage>&#x02013;<lpage>272</lpage>. <pub-id pub-id-type="doi">10.1260/175682909790291500</pub-id></citation></ref>
<ref id="B3">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bortot</surname> <given-names>E.</given-names></name> <name><surname>Gei</surname> <given-names>M.</given-names></name></person-group> (<year>2015</year>). <article-title>Harvesting energy with load-driven dielectric elastomer annular membranes deforming out-of-plane</article-title>. <source>Extreme Mech. Lett.</source> <volume>5</volume>, <fpage>62</fpage>&#x02013;<lpage>73</lpage>. <pub-id pub-id-type="doi">10.1016/j.eml.2015.09.009</pub-id></citation></ref>
<ref id="B4">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Brodsky</surname> <given-names>A. K.</given-names></name></person-group> (<year>1994</year>). <source>The Evolution of Insect Flight.</source> <publisher-loc>Oxford</publisher-loc>: <publisher-name>Oxford University Press</publisher-name>.</citation></ref>
<ref id="B5">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Burgess</surname> <given-names>S.</given-names></name> <name><surname>Ling</surname> <given-names>C.</given-names></name> <name><surname>Conn</surname> <given-names>A.</given-names></name> <name><surname>Araromi</surname> <given-names>S.</given-names></name> <name><surname>Wang</surname> <given-names>J.</given-names></name> <name><surname>Vaidyanathan</surname> <given-names>R.</given-names></name></person-group> (<year>2009</year>). <article-title>Development of a novel Electro Active Polymer (EAP) actuator for driving the wings of flapping micro air vehicle</article-title>, in <source>Computer Aided Optimum Design in Engineering</source>, eds <person-group person-group-type="editor"><name><surname>Hern&#x000E1;ndez</surname> <given-names>S.</given-names></name> <name><surname>Brebbia</surname> <given-names>C. A</given-names></name></person-group>. (<publisher-loc>Southampton, UK</publisher-loc>: <publisher-name>WIT Press</publisher-name>), <fpage>207</fpage>&#x02013;<lpage>217</lpage>.</citation></ref>
<ref id="B6">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Cao</surname> <given-names>C.</given-names></name> <name><surname>Conn</surname> <given-names>A.</given-names></name></person-group> (<year>2018</year>). <article-title>Performance optimization of a conical dielectric elastomer actuator</article-title>. <source>Actuators</source> <volume>7</volume>:<fpage>32</fpage>. <pub-id pub-id-type="doi">10.3390/act7020032</pub-id></citation></ref>
<ref id="B7">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Cao</surname> <given-names>C.</given-names></name> <name><surname>Conn</surname> <given-names>A. T.</given-names></name></person-group> (<year>2017</year>). <article-title>Elastic actuation for legged locomotion</article-title>, in <source>Electroactive Polymer Actuators and Devices (EAPAD) 2017</source> (<publisher-loc>Portland</publisher-loc>).</citation></ref>
<ref id="B8">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Carpi</surname> <given-names>F.</given-names></name> <name><surname>Anderson</surname> <given-names>I.</given-names></name> <name><surname>Bauer</surname> <given-names>S.</given-names></name> <name><surname>Frediani</surname> <given-names>G.</given-names></name> <name><surname>Gallone</surname> <given-names>G.</given-names></name> <name><surname>Gei</surname> <given-names>M.</given-names></name> <etal/></person-group>. (<year>2015</year>). <article-title>Standards for dielectric elastomer transducers</article-title>. <source>Smart Materials Struct.</source> <volume>24</volume>:<fpage>105025</fpage>. <pub-id pub-id-type="doi">10.1088/0964-1726/24/10/105025</pub-id></citation></ref>
<ref id="B9">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Carpi</surname> <given-names>F.</given-names></name> <name><surname>Kornbluh</surname> <given-names>R.</given-names></name> <name><surname>Sommer-Larsen</surname> <given-names>P.</given-names></name> <name><surname>Alici</surname> <given-names>G.</given-names></name></person-group> (<year>2011</year>). <article-title>Electroactive polymer actuators as artificial muscles: are they ready for bioinspired applications?</article-title> <source>Bioinspir. Biomim.</source> <volume>4</volume>:<fpage>045006</fpage>. <pub-id pub-id-type="doi">10.1088/1748-3182/6/4/045006</pub-id></citation></ref>
<ref id="B10">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Carpi</surname> <given-names>F.</given-names></name> <name><surname>Salaris</surname> <given-names>C.</given-names></name> <name><surname>De Rossi</surname> <given-names>D.</given-names></name></person-group> (<year>2007</year>). <article-title>Folded dielectric elastomer actuators</article-title>. <source>Smart Materials Struct.</source> <volume>16</volume>:<fpage>S300</fpage>. <pub-id pub-id-type="doi">10.1088/0964-1726/16/2/S15</pub-id></citation></ref>
<ref id="B11">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Choi</surname> <given-names>H. R.</given-names></name> <name><surname>Jung</surname> <given-names>K. M.</given-names></name> <name><surname>Kwak</surname> <given-names>J. W.</given-names></name> <name><surname>Lee</surname> <given-names>S. W.</given-names></name> <name><surname>Kim</surname> <given-names>H. M.</given-names></name> <name><surname>Jeon</surname> <given-names>J. W.</given-names></name> <etal/></person-group>. (<year>2003</year>). <article-title>Digital polymer motor for robotic applications</article-title>, in <source>Robotics and Automation, 2003. 2003 IEEE International Conference on Robotics and Automation (ICRA&#x00027;03)</source> (<publisher-loc>Taipei</publisher-loc>).</citation></ref>
<ref id="B12">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Conn</surname> <given-names>A. T.</given-names></name> <name><surname>Rossiter</surname> <given-names>J.</given-names></name></person-group> (<year>2012</year>). <article-title>Towards holonomic electro-elastomer actuators with six degrees of freedom</article-title>. <source>Smart Mater. Struct.</source> <volume>3</volume>:<fpage>035012</fpage>. <pub-id pub-id-type="doi">10.1088/0964-1726/21/3/035012</pub-id></citation></ref>
<ref id="B13">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>De Croon</surname> <given-names>G.</given-names></name> <name><surname>De Clercq</surname> <given-names>K.</given-names></name> <name><surname>Ruijsink</surname> <given-names>R.</given-names></name> <name><surname>Remes</surname> <given-names>B.</given-names></name> <name><surname>de Wagter</surname> <given-names>C.</given-names></name></person-group> (<year>2009</year>). <article-title>Design, aerodynamics, and vision-based control of the DelFly</article-title>. <source>Int. J. Micro Air Vehicles</source> <volume>1</volume>, <fpage>71</fpage>&#x02013;<lpage>97</lpage>. <pub-id pub-id-type="doi">10.1260/175682909789498288</pub-id></citation></ref>
<ref id="B14">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Dudley</surname> <given-names>R.</given-names></name></person-group> (<year>2002</year>). <source>The Biomechanics of Insect Flight: Form, Function, Evolution</source>. <publisher-loc>Princeton, NJ</publisher-loc>: <publisher-name>Princeton University Press</publisher-name>.</citation></ref>
<ref id="B15">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ellington</surname> <given-names>C. P.</given-names></name></person-group> (<year>1999</year>). <article-title>The novel aerodynamics of insect flight: applications to micro-air vehicles</article-title>. <source>J. Exp. Biol.</source> <volume>23</volume>, <fpage>3439</fpage>&#x02013;<lpage>3448</lpage>.</citation></ref>
<ref id="B16">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>He</surname> <given-names>T.</given-names></name> <name><surname>Cui</surname> <given-names>L.</given-names></name> <name><surname>Chen</surname> <given-names>C.</given-names></name> <name><surname>Suo</surname> <given-names>Z.</given-names></name></person-group> (<year>2010</year>). <article-title>Nonlinear deformation analysis of a dielectric elastomer membrane&#x02013;spring system</article-title>. <source>Smart Mater. Struct.</source> <volume>8</volume>:<fpage>085017</fpage>. <pub-id pub-id-type="doi">10.1088/0964-1726/19/8/085017</pub-id></citation></ref>
<ref id="B17">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Henke</surname> <given-names>E.</given-names></name> <name><surname>Wilson</surname> <given-names>K.</given-names></name> <name><surname>Anderson</surname> <given-names>I.</given-names></name></person-group> (<year>2017</year>). <article-title>Entirely soft dielectric elastomer robots</article-title>, in <source>Electroactive Polymer Actuators and Devices (EAPAD) 2017</source>, (Portland).</citation></ref>
<ref id="B18">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hodgins</surname> <given-names>M.</given-names></name> <name><surname>Rizzello</surname> <given-names>G.</given-names></name> <name><surname>Naso</surname> <given-names>D.</given-names></name> <name><surname>York</surname> <given-names>A.</given-names></name> <name><surname>Seelecke</surname> <given-names>S.</given-names></name></person-group> (<year>2014</year>). <article-title>An electro-mechanically coupled model for the dynamic behavior of a dielectric electro-active polymer actuator</article-title>. <source>Smart Mater. Struct.</source> <volume>10</volume>:<fpage>104006</fpage>. <pub-id pub-id-type="doi">10.1088/0964-1726/23/10/104006</pub-id></citation></ref>
<ref id="B19">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Huber</surname> <given-names>J.</given-names></name> <name><surname>Fleck</surname> <given-names>N.</given-names></name> <name><surname>Ashby</surname> <given-names>M.</given-names></name></person-group> (<year>1997</year>). <article-title>The selection of mechanical actuators based on performance indices</article-title>. <source>Royal Soc.</source> <volume>453</volume>, <fpage>2185</fpage>&#x02013;<lpage>2205</lpage>. <pub-id pub-id-type="doi">10.1098/rspa.1997.0117</pub-id></citation></ref>
<ref id="B20">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Keennon</surname> <given-names>M.</given-names></name> <name><surname>Klingebiel</surname> <given-names>K.</given-names></name> <name><surname>Won</surname> <given-names>H.</given-names></name> <name><surname>Andriukov</surname> <given-names>A.</given-names></name></person-group> (<year>2012</year>). <article-title>Development of the nano hummingbird: a tailless flapping wing micro air vehicle</article-title>, in <source>AIAA Aerospace Sciences Meeting</source> (<publisher-loc>Reston, VA</publisher-loc>).</citation></ref>
<ref id="B21">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kovacs</surname> <given-names>G.</given-names></name> <name><surname>D&#x000FC;ring</surname> <given-names>L.</given-names></name> <name><surname>Michel</surname> <given-names>S.</given-names></name> <name><surname>Terrasi</surname> <given-names>G.</given-names></name></person-group> (<year>2009</year>). <article-title>Stacked dielectric elastomer actuator for tensile force transmission</article-title>. <source>Sensors Actuators A Phys.</source> <volume>2</volume>, <fpage>299</fpage>&#x02013;<lpage>307</lpage>. <pub-id pub-id-type="doi">10.1016/j.sna.2009.08.027</pub-id></citation></ref>
<ref id="B22">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Lau</surname> <given-names>G. K.</given-names></name> <name><surname>Chin</surname> <given-names>Y. W.</given-names></name> <name><surname>La</surname> <given-names>T. G.</given-names></name></person-group> (<year>2017</year>). <article-title>Development of elastomeric flight muscles for flapping wing micro air vehicles</article-title>, in <source>Electroactive Polymer Actuators and Devices (EAPAD) 2017</source> (<publisher-loc>Portland, OR</publisher-loc>).</citation></ref>
<ref id="B23">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lau</surname> <given-names>G. K.</given-names></name> <name><surname>Lim</surname> <given-names>H. T.</given-names></name> <name><surname>Teo</surname> <given-names>J. Y.</given-names></name> <name><surname>Chin</surname> <given-names>Y. W.</given-names></name></person-group> (<year>2014</year>). <article-title>Lightweight mechanical amplifiers for rolled dielectric elastomer actuators and their integration with bio-inspired wing flappers</article-title>. <source>Smart Mater. Struct.</source> <volume>2</volume>:<fpage>025021</fpage>. <pub-id pub-id-type="doi">10.1039/9781782625209</pub-id></citation></ref>
<ref id="B24">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lehmann</surname> <given-names>F.</given-names></name> <name><surname>Dickinson</surname> <given-names>M.</given-names></name></person-group> (<year>1997</year>). <article-title>The changes in power requirements and muscle efficiency during elevated force production in the fruit fly Drosophila Melanogaster</article-title>. <source>J. Exp. Biol.</source> <volume>7</volume>, <fpage>1133</fpage>&#x02013;<lpage>1143</lpage>.</citation></ref>
<ref id="B25">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ogden</surname> <given-names>W.</given-names></name></person-group> (<year>1972</year>). <article-title>Large deformation isotropic elasticity-on the correlation of theory and experiment for incompressible rubberlike solids</article-title>. <source>Proc. R. Soc. A</source> <volume>326</volume>, <fpage>565</fpage>&#x02013;<lpage>584</lpage>. <pub-id pub-id-type="doi">10.1098/rspa.1972.0026</pub-id></citation></ref>
<ref id="B26">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Pringle</surname> <given-names>J.</given-names></name></person-group> (<year>1957</year>). <source>Insect Flight.</source> <publisher-loc>Cambridge</publisher-loc>: <publisher-name>Cambridge University Press</publisher-name>.</citation></ref>
<ref id="B27">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rizzello</surname> <given-names>G.</given-names></name> <name><surname>Hodgins</surname> <given-names>M.</given-names></name> <name><surname>Naso</surname> <given-names>D.</given-names></name> <name><surname>York</surname> <given-names>A.</given-names></name> <name><surname>Seelecke</surname> <given-names>S.</given-names></name></person-group> (<year>2015a</year>). <article-title>Modeling of the effects of the electrical dynamics on the electromechanical response of a DEAP circular actuator with a mass&#x02013;spring load</article-title>. <source>Smart Materials Struct.</source> <volume>9</volume>:<fpage>094003</fpage>. <pub-id pub-id-type="doi">10.1088/0964-1726/24/9/094003</pub-id></citation></ref>
<ref id="B28">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rizzello</surname> <given-names>G.</given-names></name> <name><surname>Hodgins</surname> <given-names>M.</given-names></name> <name><surname>Naso</surname> <given-names>D.</given-names></name> <name><surname>York</surname> <given-names>A.</given-names></name> <name><surname>Seelecke</surname> <given-names>S.</given-names></name></person-group> (<year>2015b</year>). <article-title>Dynamic modeling and experimental validation of an annular dielectric elastomer actuator with a biasing mass</article-title>. <source>J. Vibration Acoustics</source> <volume>1</volume>:<fpage>011005</fpage>. <pub-id pub-id-type="doi">10.1115/1.4028456</pub-id></citation></ref>
<ref id="B29">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Steltz</surname> <given-names>E.</given-names></name> <name><surname>Fearing</surname> <given-names>R.</given-names></name></person-group> (<year>2007</year>). <article-title>Dynamometer power output measurements of piezoelectric actuators</article-title>, in <source>2007 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS)</source> (<publisher-loc>Rome</publisher-loc>).</citation></ref>
<ref id="B30">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>T&#x000F6;pper</surname> <given-names>T.</given-names></name> <name><surname>Weiss</surname> <given-names>F.</given-names></name> <name><surname>Osmani</surname> <given-names>B.</given-names></name> <name><surname>Bippes</surname> <given-names>C.</given-names></name> <name><surname>Leung</surname> <given-names>V.</given-names></name> <name><surname>M&#x000FC;ller</surname> <given-names>B.</given-names></name></person-group> (<year>2015</year>). <article-title>Siloxane-based thin films for biomimetic low-voltage dielectric actuators</article-title>. <source>Sensors Actuators A Phys.</source> <volume>233</volume>, <fpage>32</fpage>&#x02013;<lpage>41</lpage>. <pub-id pub-id-type="doi">10.1016/j.sna.2015.06.014</pub-id></citation></ref>
<ref id="B31">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Tu</surname> <given-names>M.</given-names></name> <name><surname>Daniel</surname> <given-names>T.</given-names></name></person-group> (<year>2004</year>). <article-title>Submaximal power output from the dorsolongitudinal flight muscles of the hawkmoth Manduca sexta</article-title>. <source>J. Exp. Biol.</source> <volume>26</volume>, <fpage>4651</fpage>&#x02013;<lpage>4662</lpage>. <pub-id pub-id-type="doi">10.1242/jeb.01321</pub-id></citation></ref>
<ref id="B32">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname> <given-names>B.</given-names></name> <name><surname>Wang</surname> <given-names>Z.</given-names></name> <name><surname>He</surname> <given-names>T.</given-names></name></person-group> (<year>2016</year>). <article-title>Investigation on the viscoelastic behaviors of a circular dielectric elastomer membrane undergoing large deformation</article-title>. <source>AIP Adv.</source> <volume>12</volume>:<fpage>125127</fpage>. <pub-id pub-id-type="doi">10.1063/1.4973639</pub-id></citation></ref>
<ref id="B33">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname> <given-names>H.</given-names></name></person-group> (<year>2018</year>). <article-title>Viscoelastic analysis of a spring-connected dielectric elastomer actuator undergoing large inhomogeneous deformation</article-title>. <source>Int. J. Mech. Sci.</source> <fpage>17</fpage>&#x02013;<lpage>23</lpage>. <pub-id pub-id-type="doi">10.1016/j.ijmecsci.2017.12.004</pub-id></citation></ref>
<ref id="B34">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Weiss</surname> <given-names>F.</given-names></name> <name><surname>T&#x000F6;pper</surname> <given-names>T.</given-names></name> <name><surname>Osmani</surname> <given-names>B.</given-names></name> <name><surname>Peters</surname> <given-names>S.</given-names></name> <name><surname>Kovacs</surname> <given-names>G.</given-names></name> <name><surname>M&#x000FC;ller</surname> <given-names>B.</given-names></name></person-group> (<year>2016</year>). <article-title>Electrospraying nanometer-thin elastomer films for low-voltage dielectric actuators</article-title>. <source>Adv Electron. Mater.</source> <volume>2</volume>:<fpage>1500476</fpage>. <pub-id pub-id-type="doi">10.1002/aelm.201500476</pub-id></citation></ref>
<ref id="B35">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wood</surname> <given-names>R. J.</given-names></name></person-group> (<year>2008</year>). <article-title>The first take off of a biologically inspired at-scale robotic insect</article-title>, in <source>IEEE Transcations on Robotics</source> <volume>24</volume>, <fpage>341</fpage>&#x02013;<lpage>347</lpage>. <pub-id pub-id-type="doi">10.1109/TRO.2008.916997</pub-id></citation></ref>
<ref id="B36">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yuan</surname> <given-names>W. H. L.</given-names></name> <name><surname>Yu</surname> <given-names>Z. B.</given-names></name> <name><surname>Lam</surname> <given-names>T.</given-names></name> <name><surname>Biggs</surname> <given-names>J.</given-names></name> <name><surname>Ha</surname> <given-names>S. M.</given-names></name> <etal/></person-group>. (<year>2008</year>). <article-title>Fault-tolerant dielectric elastomer actuators using single-walled carbon nanotube electrodes</article-title>. <source>Adv. Materials</source> <volume>3</volume>, <fpage>612</fpage>&#x02013;<lpage>625</lpage>. <pub-id pub-id-type="doi">10.1002/adma.200701018</pub-id></citation></ref>
<ref id="B37">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname> <given-names>C.</given-names></name> <name><surname>Rossi</surname> <given-names>C.</given-names></name></person-group> (<year>2017</year>). <article-title>A review of compliant transmission mechanisms for bio-inspired flapping-wing micro air vehicles</article-title>. <source>Bioinspir. Biomim.</source> <volume>12</volume>:<fpage>p025005</fpage>. <pub-id pub-id-type="doi">10.1088/1748-3190/aa58d3</pub-id><pub-id pub-id-type="pmid">28079026</pub-id></citation></ref>
<ref id="B38">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname> <given-names>J.</given-names></name> <name><surname>Deng</surname> <given-names>X.</given-names></name></person-group> (<year>2017</year>). <article-title>Resonance principle for the design of flapping wing micro air vehicles</article-title>. <source>IEEE Trans. Robot.</source> <volume>33</volume>, <fpage>183</fpage>&#x02013;<lpage>197</lpage>. <pub-id pub-id-type="doi">10.1109/TRO.2016.2626457</pub-id></citation></ref>
</ref-list>
</back>
</article>