# Application and Analysis of an Ionic Liquid Gel in a Soft Robot

Authors: Chenghong Zhang; Bin He; Zhipeng Wang; Yanmin Zhou; Aiguo Ming
Journal: Advances in Materials Science and Engineering
 (2019)

---

## Abstract

Due to their light weight, flexibility, and low energy consumption, ionic electroactive polymers have become a hotspot for bionic soft robotics and are ideal materials for the preparation of soft actuators. Because the traditional ionic electroactive polymers, such as ionic polymer-metal composites (IPMCs), contain water ions, a soft actuator does not work properly upon the evaporation of water ions. An ionic liquid polymer gel is a new type of ionic electroactive polymer that does not contain water ions, and ionic liquids are more thermally and electrochemically stable than water. These liquids, with a low melting point and a high ionic conductivity, can be used in ionic electroactive polymer soft actuators. An ionic liquid gel (ILG), a new type of soft actuator material, was obtained by mixing 1-butyl-3-methylimidazolium tetrafluoroborate (BMIMBF4), hydroxyethyl methacrylate (HEMA), diethoxyacetophenone (DEAP) and ZrO2 and then polymerizing this mixture into a gel state under ultraviolet (UV) light irradiation. An ILG soft actuator was designed, the material preparation principle was expounded, and the design method of the soft robot mechanism was discussed. Based on nonlinear finite element theory, the deformation mechanism of the ILG actuator was deeply analyzed and the deformation of the soft robot when grabbing an object was also analyzed. A soft robot was designed with the soft actuator as the basic module. The experimental results show that the ILG soft robot has good driving performance, and the soft robot can grab a 105 mg object at an input voltage of 3.5 V.


---

## Body

## 1. Introduction

Traditional robots are typically constructed of rigid motion joints based on hard materials and can perform tasks quickly and accurately. However, such robots have limited movement flexibility and low environmental adaptability and cannot work under space-constrained conditions. These factors limit the application of rigid robots in unstructured and complex environments [1–3].In nature, mollusks are widely distributed in oceans, rivers, and lakes and on land. After hundreds of millions of years of evolution, this animal has gradually developed the characteristics of a large deformation ability, a light weight, and a high power density ratio and can achieve efficient movement under complex natural environment conditions by changing its body shape. In recent years, researchers have attempted to apply the biological principles of mollusks to the research and design of robotics. A soft robot is composed of a soft material that can withstand large strains, has the capabilities for continuous deformation and drive-structure integration and can arbitrarily change its shape and size over a wide range. Such a robot has strong adaptability to unstructured environments and broad application prospects in military reconnaissance, rescue, and medical surgery [4–6].A flexible actuator based on an ionic electroactive polymer (EAP) has the advantages of a low driving voltage, a large displacement ability, a light weight, etc. Such actuators have become a research hotspot in the field of bionic robots. In the past few decades, electrochemical actuators, which are substitutes for air- and fluid-derived devices, have been further developed due to their ideal mechanical properties for use in intelligent robots. However, traditional ionic EAP actuators, such as those based on IPMCs, have the characteristics of a short working time in nonwater media, a complex manufacturing process, and a high cost. An ILG is a new type of ionic EAP that can work in air for a long time. Because an ILG offers chemical stability, thermal stability, and simpler ion transport, it is more suitable for the production of soft robot actuators than other EAPs.In this paper, a soft robot actuator based on an ILG material is proposed and its preparation, driving mechanism, and design method are deeply analyzed. Finally, the validity and rationality of the soft robot are verified by driving performance and grabbing experiments.


## 2. Design of the ILG Soft Robot

The ILG soft actuator consists of a 5-layer structure, as shown in Figure1. The middle layer is an electroactive layer composed of the ILG material, which functions to store ionic liquids. The outer two layers that wrap the middle layer are electrode layers, which consist of activated carbon. Activated carbon has a high specific surface area, a high electrical conductivity, a high density, a strong adsorption force, etc., making it very suitable for use as an electrode material [7]. In addition, gold foil is selected as a current collector to cover the surface of the activated carbon layers. When the actuator is working, one end is attached to an external metal electrode, and a wire is led from the external electrode to connect to a power supply.Figure 1
Five-layer structure of the soft actuator.The soft robot consists of three ILG actuators, a common copper electrode, and three independent copper electrodes, as shown in Figure2. The input unit consists of a common electrode and an independent electrode. The common electrode is a copper column, with the upper end fixed with a fixture, and the copper electrode is connected to a power source. The soft robot consists of three ILG actuators evenly distributed around the common electrode at 120°. Two copper electrodes sandwich one ILG actuator.Figure 2
Configuration of the soft robot.Figure3 shows the motion cycle of the soft robot to grab an object, which can be described as follows. First, the ILG actuator is naturally stretched, and by applying an input voltage of 3.5 V to the electrodes, the actuator quickly bends outward. Second, when the actuator is close to the target object, the input voltage direction is changed, and the actuator bends inward to clamp the object. Finally, the input signal is changed, and the actuator bends outward to release the object.Figure 3
Motion cycle of the soft robot manipulator to grab an object.


### 2.1. Preparation of the Soft Actuator

The preparation process of an ILG is a process of ionic liquid loading. Ionic liquid loading refers to the immobilization of an ionic liquid to form a solid carrier by physical or chemical means [8]. BMIMBF4 was selected as the ionic liquid. The carriers of ionic liquids need to have interconnected porous structures, large specific surface area and porosity, and good mechanical strength and electrochemical stability. Therefore, HEMA was selected to prepare the carrier. This material is a colorless and transparent liquid that is soluble in water and has low toxicity. HEMA was polymerized under ultraviolet light to form polyhydroxyethyl methacrylate (PHEMA) with a porous structure.The mixed solution consists of BMIMBF4, HEMA, and ZrO2. Under UV light, the ILG is generated by the polymerization reaction induced by DEAP. The polymer matrix is cross-linked into a porous network structure. ZrO2 enhances the mechanical strength of the gel, and as the ZrO2 content increases, the tensile strength of the ILG increases [9, 10]. Figure 4 shows the soft actuator fabrication process. Figure 4(a) shows the ILG sample, which is covered with a layer of BMIMBF4 due to osmosis. Because the ionic liquid has a certain viscosity, the activated carbon electrode can be directly affixed to both sides of the gel, as shown in Figure 4(b). Finally, the gold foil is evenly affixed to the surface of the activated carbon layers to obtain an ILG soft actuator, as shown in Figure 4(c).Figure 4
Fabrication process of the soft actuator: (a) ILG; (b) the activated carbon layer is affixed; (c) the gold foil layer is affixed.
(a)
(b)
(c)


### 2.2. Soft Actuator Deformation Mechanism

For the ILG actuator, the application of a voltage on the electrode can cause bending deformation, which is the result of the coupling of the electric field, chemical field, flow field, and field force. The average pore diameter of the structure is greater than 10 microns, which is much larger than the diameter of free ions in the ionic liquids (nanoscale). Therefore, BMIM+ and BF4− can move freely inside the carrier. Under the action of the electric field, positive and negative ions accumulate in the electrode layers on both sides of the ILG. Because the cation volume is much larger than the anion volume, the volume of the ILG cathode side expands and the volume of the ILG anode side shrinks [11, 12]. The entire actuator therefore bends toward the anode side, as shown in Figure 5.Figure 5
Working principle of the actuator.When an electric field is applied, the movement of the internal ions of the ILG actuator in the porous medium can be described by the Nernst–Planck equation [13]:(1)ji=−di∇ϖi+FRTϖizi∇Θ,where i is the ion type (i = 1 represents the cation, and i = 2 represents the anion), ji is the flow of ions in the PHEMA pores, di is the diffusion coefficient of the ions, ϖi is the ion concentration, ∇ϖi is the ion concentration gradient, Θ is the potential, ∇Θ is the potential gradient, F is the Faraday constant, R is the gas constant, T is the absolute temperature, and zi is the number of ionic charges. The first term on the right-hand side of the formula indicates the contribution of the cation or anion concentration diffusion; the second term indicates the contribution of electromigration.The continuous equation of concentration changing with time is(2)∂ϖi∂t=−∇ji,where ∇ji is the flow change of the i‐th ion in the PHEMA pores and t is the time.The internal electric field balance equation of the ILG actuator can be expressed as follows:(3)E=Dke=−∇Θ,∇D=ρ=Fϖ1−ϖ2,where D is the electric displacement, E is the electric field, ρ is the net charge density, and ke is the effective dielectric constant of the ILG.The effect of the electrode layer on the bending deformation of the actuator along the width and the thickness is ignored. Considering only the relationship between the axial deformation of the electrode layer and the net charge density, the axial induced strain can be expressed as follows:(4)εcs=ρsαs,where εc is the induced strain, αs is the induced strain coupling coefficient, and s is the Laplace domain.Therefore, the axial induced stress is(5)σcs=εcEt=ρsds,ds=αsEt,where Et is the elastic modulus of the electrode layer and ds is the axial-induced stress coupling coefficient.The bending moment of the actuator is(6)Ms=2∫b1bσcsxwdx=2∫b1bρsdsxwdx,where b and b1 are half the thickness of the actuator and ILG layer, respectively, and w is the width of the actuator.The ILG actuator is the basic unit of a soft robot and is the source of its deformation. For flexible mechanisms, the degrees of freedom and constraints are determined by the flexibility matrix. The cantilever beam flexibility matrix is as follows [14]:(7)R=diagRθx,Rθy,Rθz,Rx,Ry,Rz=diaglEaIx,lEaIy,lGaJ,l312EaIy,l312EaIx,lEaA,where Rθx, Rθy, and Rθz represent rotational flexibility around the x, y, and z axes, respectively, and Rx, Ry, and Rz represent translational flexibility along the x, y,and z axes, respectively. If the order of magnitude of an element in one direction is significantly greater than that of the element in the other directions, then the mechanism has a degree of freedom in this direction. Using this method, the equivalent constraints or degrees of freedom of the soft robot can be established.


## 2.1. Preparation of the Soft Actuator

The preparation process of an ILG is a process of ionic liquid loading. Ionic liquid loading refers to the immobilization of an ionic liquid to form a solid carrier by physical or chemical means [8]. BMIMBF4 was selected as the ionic liquid. The carriers of ionic liquids need to have interconnected porous structures, large specific surface area and porosity, and good mechanical strength and electrochemical stability. Therefore, HEMA was selected to prepare the carrier. This material is a colorless and transparent liquid that is soluble in water and has low toxicity. HEMA was polymerized under ultraviolet light to form polyhydroxyethyl methacrylate (PHEMA) with a porous structure.The mixed solution consists of BMIMBF4, HEMA, and ZrO2. Under UV light, the ILG is generated by the polymerization reaction induced by DEAP. The polymer matrix is cross-linked into a porous network structure. ZrO2 enhances the mechanical strength of the gel, and as the ZrO2 content increases, the tensile strength of the ILG increases [9, 10]. Figure 4 shows the soft actuator fabrication process. Figure 4(a) shows the ILG sample, which is covered with a layer of BMIMBF4 due to osmosis. Because the ionic liquid has a certain viscosity, the activated carbon electrode can be directly affixed to both sides of the gel, as shown in Figure 4(b). Finally, the gold foil is evenly affixed to the surface of the activated carbon layers to obtain an ILG soft actuator, as shown in Figure 4(c).Figure 4
Fabrication process of the soft actuator: (a) ILG; (b) the activated carbon layer is affixed; (c) the gold foil layer is affixed.
(a)
(b)
(c)


## 2.2. Soft Actuator Deformation Mechanism

For the ILG actuator, the application of a voltage on the electrode can cause bending deformation, which is the result of the coupling of the electric field, chemical field, flow field, and field force. The average pore diameter of the structure is greater than 10 microns, which is much larger than the diameter of free ions in the ionic liquids (nanoscale). Therefore, BMIM+ and BF4− can move freely inside the carrier. Under the action of the electric field, positive and negative ions accumulate in the electrode layers on both sides of the ILG. Because the cation volume is much larger than the anion volume, the volume of the ILG cathode side expands and the volume of the ILG anode side shrinks [11, 12]. The entire actuator therefore bends toward the anode side, as shown in Figure 5.Figure 5
Working principle of the actuator.When an electric field is applied, the movement of the internal ions of the ILG actuator in the porous medium can be described by the Nernst–Planck equation [13]:(1)ji=−di∇ϖi+FRTϖizi∇Θ,where i is the ion type (i = 1 represents the cation, and i = 2 represents the anion), ji is the flow of ions in the PHEMA pores, di is the diffusion coefficient of the ions, ϖi is the ion concentration, ∇ϖi is the ion concentration gradient, Θ is the potential, ∇Θ is the potential gradient, F is the Faraday constant, R is the gas constant, T is the absolute temperature, and zi is the number of ionic charges. The first term on the right-hand side of the formula indicates the contribution of the cation or anion concentration diffusion; the second term indicates the contribution of electromigration.The continuous equation of concentration changing with time is(2)∂ϖi∂t=−∇ji,where ∇ji is the flow change of the i‐th ion in the PHEMA pores and t is the time.The internal electric field balance equation of the ILG actuator can be expressed as follows:(3)E=Dke=−∇Θ,∇D=ρ=Fϖ1−ϖ2,where D is the electric displacement, E is the electric field, ρ is the net charge density, and ke is the effective dielectric constant of the ILG.The effect of the electrode layer on the bending deformation of the actuator along the width and the thickness is ignored. Considering only the relationship between the axial deformation of the electrode layer and the net charge density, the axial induced strain can be expressed as follows:(4)εcs=ρsαs,where εc is the induced strain, αs is the induced strain coupling coefficient, and s is the Laplace domain.Therefore, the axial induced stress is(5)σcs=εcEt=ρsds,ds=αsEt,where Et is the elastic modulus of the electrode layer and ds is the axial-induced stress coupling coefficient.The bending moment of the actuator is(6)Ms=2∫b1bσcsxwdx=2∫b1bρsdsxwdx,where b and b1 are half the thickness of the actuator and ILG layer, respectively, and w is the width of the actuator.The ILG actuator is the basic unit of a soft robot and is the source of its deformation. For flexible mechanisms, the degrees of freedom and constraints are determined by the flexibility matrix. The cantilever beam flexibility matrix is as follows [14]:(7)R=diagRθx,Rθy,Rθz,Rx,Ry,Rz=diaglEaIx,lEaIy,lGaJ,l312EaIy,l312EaIx,lEaA,where Rθx, Rθy, and Rθz represent rotational flexibility around the x, y, and z axes, respectively, and Rx, Ry, and Rz represent translational flexibility along the x, y,and z axes, respectively. If the order of magnitude of an element in one direction is significantly greater than that of the element in the other directions, then the mechanism has a degree of freedom in this direction. Using this method, the equivalent constraints or degrees of freedom of the soft robot can be established.


## 3. Hyperelastic Arruda–Boyce Model

Material nonlinearity is caused by the nonlinearity of the constitutive relationship of the material. The constitutive relationship of a material includes the stress, strain, strain rate, load duration, temperature, etc. After a load is unloaded, residual strain will appear in a nonlinear material, as shown in Figure6.Figure 6
Stress-strain curve of nonlinear materials.The hyperelastic Arruda–Boyce model, typically used for relatively largely deformed materials, is used to simulate ILGs. Hyperelastic materials can be described by the strain energy functionW, which is used to define the strain energy density in the material. The Arruda–Boyce strain energy function W is as follows [15–17]:(8)W=G∑k=15ckλm2k−1I1k−3k+1DJel2−12−lnJel,with c1=1/2, c2=1/20, c3=11/1050, c4=19/7000, and c5=519/673750, where W is the strain energy density, G is the initial shear modulus, λm is the deformation rate, Jel is the elastic volume ratio, and D is the temperature-dependent material parameter associated with the bulk modulus. A material is completely incompressible when Jel = 1.I
1 is the first strain invariant of the deviatoric strain, and the relationship between I1 and the main tensile strain rates λ1, λ2, and λ3 is as follows:(9)I1=λ12+λ22+λ32,I2=1λ12+1λ22+1λ32,I3=λ12×λ22×λ32,where λ1, λ2, and λ3 are the main (extension) deformation rates along the x, y, and z directions, respectively. I1, I2, and I3 are the relative changes of the length, the surface area, and the volume of the elastomer, respectively.The material stress is obtained by differentiating the stretching:(10)σn=∂W∂λn,where σn is the normal stress and λn is the stretching in the loading direction.The isotropic and incompressible deformation process of an ILG is given by(11)I3=λ1λ2λ3=1.When the ILG is uniformly stretched,λ2=λ3; equation (11) can be used to calculate(12)λ2=λ3=1λ1.The true stress is obtained by(13)σ0=σnλn.Differentiating equation (8), we obtain(14)σ0=Gλn2−λn−1+15λm2λn4−λn−2+11175λm4λn6−λn−3+19875λm6λn8−λn−4+51967375λm8λn10−λn−5.


### 4. Model for the Soft Robot

To establish a gripping model scheme as shown in Figure7, assume that the soft robot is an inextensible, fully flexible elastic rod. Assume that l is the length of the evenly distributed rod, EI is the flexural rigidity, ρ is the mass per unit length, and kg is the curvature profile when the object is gripped [18, 19].Figure 7
Lifting of an object, considering the normal force and the static friction force at the contact point of the soft actuator.


### 4.1. Deformation Analysis for Soft Actuators

Ats=s∗, the position of any material point can be represented by the following formula [20], as shown in Figure 8.(15)r⇀s=s∗=Xs=s∗E⇀1+Ys=s∗E⇀2,where, in Cartesian coordinates, X and Y are represented as follows:(16)Xs=s∗=Xs=0+∫0s∗cosθηdη,Ys=s∗=Ys=0+∫0s∗sinθηdη,where η is a dummy variable and θ is the angle between the unit tangent vector and the horizontal direction.(17)r′⇀=cosθsE⇀1+sinθsE⇀2,where the prime (′) represents the partial derivative with respect to s. Assuming that r⇀ is continuous, θ and r′⇀ are also continuous.Figure 8
Ats = 0, F0 and M0 are the terminal force and terminal moment, respectively; at s = l, Fl and Ml are the terminal force and terminal moment; and at point s = ξ, Fξ is the force.At points=ς, the jump in any function χ=χs,θs,θ′s is expressed by a compact notation:(18)χς=χς+−χς−,where(19)χς−=lims↗ςχs,θs,θ′s,χς+=lims↘ςχs,θs,θ′s.The jump will be associated with the force applied at the discrete point.The bending momentM is described by a constitutive equation:(20)M⇀=EIθ′−k0E⇀3,where k0 represents the curvature. In addition to the gravity per unit length and end load of the rod, the singular force Fξ⇀ at s=ξ must be considered. Figure 7 shows the force that simulates the contact of an object with the soft robot. The governing equation can be derived from the equilibrium of the linear and angular momentum [21–23]:(21)m′⇀−ρgE2⇀=0,∂∂sEIθ′−k0+m2cosθ−m1sinθ=0,m⇀ξ+Fξ⇀=0,where m⇀=m1E⇀1+m2E⇀2.


### 4.2. Deformation Analysis Governing Formula for the Soft Robot

The corner of the object is assumed to be in contact with the soft robot ats = ξ. Two new unit vectors are defined as follows, as shown in Figure 9:(22)H⇀1=cosθξE⇀1+sinθξE⇀2,H⇀1=−sinθξE⇀1+cosθξE⇀2,where H⇀1 is the tangent to the centerline of the soft robot and H⇀2 is perpendicular to the centerline of the soft robot.(23)F⇀ξ=f1H⇀1+f2H⇀2=N1E⇀1+N2E⇀2,where f1 is the normal force and f2 is the corresponding friction force.Figure 9
At points = ξ, the object is in contact with the soft robot. The object is assumed to remain stationary due to the static friction force f1H⇀1 and the normal force f2H⇀2.Equation (21) is applied to segment s∈ξ,l, and we obtain the following equation:(24)m⇀ξ+=−ρgl−ξE⇀2,m⇀ξ−=N1E⇀1+N2E⇀2−ρgl−ξE⇀2.The expression of the potential energyW of the soft robot is established:(25)W=∫0ξEI2θ′−k02+ρgYs−m⇀ξ−⋅r′⇀ds+∫ξlEI2θ′−k02+ρgYsds,where(26)Ys=∫0ssinθξdξ.Equation (21) can be used to establish boundary value problems.(27)EIθ″−k0′−ρgl−scosθ=0,s∈0,ξ,EIθ″−k0′−ρgl−scosθ+N2cosθ−N1sinθ=0,s∈ξ,l,where the solution θ=θ∗ satisfies the conditions(28)θ0=0,θ′l=k0l,θξ=0,θ′ξ=0,∫0ξcosθsds=c2.The soft robot is divided intoK segments. The total potential energy expression is approximated as follows:(29)Wd=EI2∑i=1Kθi−θi−1ds−k02ds+ρgl−s0sinθ0ds2+∑i=1K−1ρgl−sisinθids+ρgl−sKsinθnds2−N1cosθ0ds2−N1∑i=1I−1cosθids−N1cosθIds2−N2sinθ0ds2−N2∑i=1I−1sinθids−N2sinθIds2,where the constraint function Wc is(30)Wc=ds2+∑i=1I−1cosθids+cosθIds2−c2.


## 4. Model for the Soft Robot

To establish a gripping model scheme as shown in Figure7, assume that the soft robot is an inextensible, fully flexible elastic rod. Assume that l is the length of the evenly distributed rod, EI is the flexural rigidity, ρ is the mass per unit length, and kg is the curvature profile when the object is gripped [18, 19].Figure 7
Lifting of an object, considering the normal force and the static friction force at the contact point of the soft actuator.


## 4.1. Deformation Analysis for Soft Actuators

Ats=s∗, the position of any material point can be represented by the following formula [20], as shown in Figure 8.(15)r⇀s=s∗=Xs=s∗E⇀1+Ys=s∗E⇀2,where, in Cartesian coordinates, X and Y are represented as follows:(16)Xs=s∗=Xs=0+∫0s∗cosθηdη,Ys=s∗=Ys=0+∫0s∗sinθηdη,where η is a dummy variable and θ is the angle between the unit tangent vector and the horizontal direction.(17)r′⇀=cosθsE⇀1+sinθsE⇀2,where the prime (′) represents the partial derivative with respect to s. Assuming that r⇀ is continuous, θ and r′⇀ are also continuous.Figure 8
Ats = 0, F0 and M0 are the terminal force and terminal moment, respectively; at s = l, Fl and Ml are the terminal force and terminal moment; and at point s = ξ, Fξ is the force.At points=ς, the jump in any function χ=χs,θs,θ′s is expressed by a compact notation:(18)χς=χς+−χς−,where(19)χς−=lims↗ςχs,θs,θ′s,χς+=lims↘ςχs,θs,θ′s.The jump will be associated with the force applied at the discrete point.The bending momentM is described by a constitutive equation:(20)M⇀=EIθ′−k0E⇀3,where k0 represents the curvature. In addition to the gravity per unit length and end load of the rod, the singular force Fξ⇀ at s=ξ must be considered. Figure 7 shows the force that simulates the contact of an object with the soft robot. The governing equation can be derived from the equilibrium of the linear and angular momentum [21–23]:(21)m′⇀−ρgE2⇀=0,∂∂sEIθ′−k0+m2cosθ−m1sinθ=0,m⇀ξ+Fξ⇀=0,where m⇀=m1E⇀1+m2E⇀2.


## 4.2. Deformation Analysis Governing Formula for the Soft Robot

The corner of the object is assumed to be in contact with the soft robot ats = ξ. Two new unit vectors are defined as follows, as shown in Figure 9:(22)H⇀1=cosθξE⇀1+sinθξE⇀2,H⇀1=−sinθξE⇀1+cosθξE⇀2,where H⇀1 is the tangent to the centerline of the soft robot and H⇀2 is perpendicular to the centerline of the soft robot.(23)F⇀ξ=f1H⇀1+f2H⇀2=N1E⇀1+N2E⇀2,where f1 is the normal force and f2 is the corresponding friction force.Figure 9
At points = ξ, the object is in contact with the soft robot. The object is assumed to remain stationary due to the static friction force f1H⇀1 and the normal force f2H⇀2.Equation (21) is applied to segment s∈ξ,l, and we obtain the following equation:(24)m⇀ξ+=−ρgl−ξE⇀2,m⇀ξ−=N1E⇀1+N2E⇀2−ρgl−ξE⇀2.The expression of the potential energyW of the soft robot is established:(25)W=∫0ξEI2θ′−k02+ρgYs−m⇀ξ−⋅r′⇀ds+∫ξlEI2θ′−k02+ρgYsds,where(26)Ys=∫0ssinθξdξ.Equation (21) can be used to establish boundary value problems.(27)EIθ″−k0′−ρgl−scosθ=0,s∈0,ξ,EIθ″−k0′−ρgl−scosθ+N2cosθ−N1sinθ=0,s∈ξ,l,where the solution θ=θ∗ satisfies the conditions(28)θ0=0,θ′l=k0l,θξ=0,θ′ξ=0,∫0ξcosθsds=c2.The soft robot is divided intoK segments. The total potential energy expression is approximated as follows:(29)Wd=EI2∑i=1Kθi−θi−1ds−k02ds+ρgl−s0sinθ0ds2+∑i=1K−1ρgl−sisinθids+ρgl−sKsinθnds2−N1cosθ0ds2−N1∑i=1I−1cosθids−N1cosθIds2−N2sinθ0ds2−N2∑i=1I−1sinθids−N2sinθIds2,where the constraint function Wc is(30)Wc=ds2+∑i=1I−1cosθids+cosθIds2−c2.


## 5. Nonlinear Finite Element Method

Geometric nonlinearity arises from the nonlinear relationship between the strain and displacement. Currently, research on geometric nonlinearity mainly focuses on three types of problems: (1) a large displacement with a small strain, (2) a small displacement with a large strain, and (3) a large displacement with a large strain. A geometric nonlinear problem has two main characteristics. First, due to the large deformation of the structure, the strain and displacement of the structure are nonlinear. Second, a balance equation is established at the position after deformation. In the analysis of large deformation, the displacement of the structure changes continuously and appropriate strain, stress, and constitutive relationships should be adopted.To capture the nonlinear behavior of the structure, the full nonlinear finite element formulas of the truss elements and beam elements are studied.


### 5.1. Nonlinear FEM for Truss Elements

In the Cartesian coordinate system, the object is displaced to a certain position under the action of external forces, as shown in Figure10. (x1, y1, z1) and (x2, y2, z2) are the position coordinates of points P1 and P2 before deformation, respectively. The object is deformed to a new position under the action of external forces. (u1, v1, w1) and (u2, v2, w2) are the deformation coordinates of points P1 and P2 after deformation, respectively.Figure 10
Nondeformed and deformed geometry of a truss element.The formula for the undeformed length of a truss element is as follows:(31)l=x2−x12+y2−y12+z2−z12.The expression of the total Lagrangian strainε along the deformation axis is(32)ε=u2+x2−u1−x1+v2+y2−v1−y1+w2+z2−w1−z1−ll.According to the change of the elastic energy, the product of the stiffness matrix and the displacement vector is expressed as follows [2, 24]:(33)ku=AlEεφ,φ=∂ε∂u1,∂ε∂v1,∂ε∂w1,∂ε∂u2,∂ε∂v2,∂ε∂w2T,where E is Young’s modulus, A is the cross sectional area, and l is the length of the truss element.The tangent stiffness matrix is obtained by differentiating equation (33) with respect to the displacement vector:(34)k⌣=∂ku∂u=EAlφ∂ε∂u+EAlε∂φ∂u=EAlφφT+EAlε∂2φ∂u2.The mass matrix is the same as that for a linear truss element [25, 26]:(35)m=ρAl6200100020010002001100200010020001002.


### 5.2. Nonlinear FEM for Beam Elements

The two famous beam theories are the Euler–Bernoulli and Timoshenko beam theories. The Euler–Bernoulli beam theory assumes that the cross section remains planar and normal to the reference line after bending and its stiffness is higher than the actual stiffness. Timoshenko’s beam theory overcomes this problem by introducing shear deformation into the model, which obtains accurate results for thick beam calculations. Because the soft actuators are slender, the shear deformation is negligible. Therefore, nonlinear Euler–Bernoulli beams with von Karman nonlinearity can be used for modeling and analysis of the soft actuators.In addition to the bending effect, the finite element formula should also include the torsion and tensile effects to reflect the large deformation effect [27, 28]:(36)ε11=e+zρ2−yρ3,ε12=−zρ1,ε13=yρ1,e=u′+v′22+w′22,ρ1=φ′,ρ2=−w″1−w′2,ρ3=−v1−v2,where εij is the engineering strain tensor, e is the axial strain, y and z are the coordinates on the cross section, and ρi is the deformation curvature. u, v, and w are the displacements on the cross section.When a change of energy is applied, the product of the stiffness matrix and the displacement vector is expressed as follows:(37)kh=∫lBTΨTΦϕdx,where(38)B=∂N∂s,Ψ=1v′w′000000100002w′w″00−1−w′20−2v′v″001−v′20,Φ=EA0000GI110000EI220000EI33,ϕ=eρ1ρ2ρ3.Assume that(39)hi=h¯+Δhi.By applying Taylor expansion without higher order terms, the tangential stiffness matrix is expressed as follows:(40)k⌣=∫lBTΨ¯TΦΨ+ΓBdx,where(41)Γ=0000000EAe−2EI33v″ρ300−2EI33v′ρ3000EAe+2EI22w″ρ2002EI33w′ρ2000000000000000000.The mass matrix is the same as that of a linear beam element.


### 5.3. Nonlinear Convergence Algorithms

Although the Newton–Raphson algorithm is used for nonlinear static analysis, the tangent stiffness matrix becomes singular at some points, which makes a global equilibrium solution impossible. Riks proposed a process to track the intersection of the normal to the tangent line with the equilibrium path to solve this problem [29], as shown in Figure 11. The Crisfield method uses arcs instead of vertical lines to search for solutions. The increment of the load factor becomes an unknown problem to be solved in the iterative process.Figure 11
Scheme of the Crisfield method.Assume that(42)F=λF⇀,where λ is the load parameter and F⇀ is the preselected load vector.The residual force vector of thei‐th step is expressed as follows:(43)RU⇀i,λi=RU⇀i−1,λi−1+∂R∂U⇀U⇀i−1,λi−1δU⇀i+∂R∂λU⇀i−1,λi−1δλ1i,where U⇀ is the displacement vector.(44)kU⇀U⇀i−1,λi−1−λi−1F⇀+kTU⇀i−1,λi−1δU⇀1i+F⇀δλ1i=0,where δU⇀1i is the first increment of the load parameter in the i‐th step and δλ1i is the first increment of the displacement vector in the i‐th step.(45)U⇀1i=U⇀i−1+δU⇀1i,(46)λ1i=λi−1+δλ1i.Substituting (45) and (46) into equation (44), we obtain(47)kU⇀U⇀1i,λ1i−λ1iF⇀+kTU⇀1i,λ1iδU⇀2i+F⇀δλ2i=0.When the second step incremental search path is perpendicular to the normal of the first step incremental path, we obtain(48)δU⇀1i⋅δU⇀2i+F⇀δλ1i⋅F⇀δλ2i=0.When solving nonlinear finite element equations by an iterative method, choosing the appropriate convergence criterion is necessary to ensure that the iteration can be terminated. The convergence criterion will directly affect the speed and accuracy of the solution. If the convergence criterion is not appropriately chosen, the calculation will fail.From equations (47) and (48), we obtain(49)δU⇀2i=−kTU⇀1i,λ1i−1kU⇀U⇀1i,λ1i−λ1iF⇀+F⇀δλ2i,δλ2i=δU⇀1ikTU⇀1i,λ1i−1kU⇀U⇀1i,λ1i−λ1iF⇀F⇀δλ1iF⇀−δU⇀1ikTU⇀1i,λ1i−1F⇀.The convergent load parameters and displacement vectors are as follows:(50)λi=λi−1+δλ1i+δλ2i+δλ3i+⋯,U⇀i=U⇀i−1+δU⇀1i+δU⇀2i+δU⇀3i+⋯.


## 5.1. Nonlinear FEM for Truss Elements

In the Cartesian coordinate system, the object is displaced to a certain position under the action of external forces, as shown in Figure10. (x1, y1, z1) and (x2, y2, z2) are the position coordinates of points P1 and P2 before deformation, respectively. The object is deformed to a new position under the action of external forces. (u1, v1, w1) and (u2, v2, w2) are the deformation coordinates of points P1 and P2 after deformation, respectively.Figure 10
Nondeformed and deformed geometry of a truss element.The formula for the undeformed length of a truss element is as follows:(31)l=x2−x12+y2−y12+z2−z12.The expression of the total Lagrangian strainε along the deformation axis is(32)ε=u2+x2−u1−x1+v2+y2−v1−y1+w2+z2−w1−z1−ll.According to the change of the elastic energy, the product of the stiffness matrix and the displacement vector is expressed as follows [2, 24]:(33)ku=AlEεφ,φ=∂ε∂u1,∂ε∂v1,∂ε∂w1,∂ε∂u2,∂ε∂v2,∂ε∂w2T,where E is Young’s modulus, A is the cross sectional area, and l is the length of the truss element.The tangent stiffness matrix is obtained by differentiating equation (33) with respect to the displacement vector:(34)k⌣=∂ku∂u=EAlφ∂ε∂u+EAlε∂φ∂u=EAlφφT+EAlε∂2φ∂u2.The mass matrix is the same as that for a linear truss element [25, 26]:(35)m=ρAl6200100020010002001100200010020001002.


## 5.2. Nonlinear FEM for Beam Elements

The two famous beam theories are the Euler–Bernoulli and Timoshenko beam theories. The Euler–Bernoulli beam theory assumes that the cross section remains planar and normal to the reference line after bending and its stiffness is higher than the actual stiffness. Timoshenko’s beam theory overcomes this problem by introducing shear deformation into the model, which obtains accurate results for thick beam calculations. Because the soft actuators are slender, the shear deformation is negligible. Therefore, nonlinear Euler–Bernoulli beams with von Karman nonlinearity can be used for modeling and analysis of the soft actuators.In addition to the bending effect, the finite element formula should also include the torsion and tensile effects to reflect the large deformation effect [27, 28]:(36)ε11=e+zρ2−yρ3,ε12=−zρ1,ε13=yρ1,e=u′+v′22+w′22,ρ1=φ′,ρ2=−w″1−w′2,ρ3=−v1−v2,where εij is the engineering strain tensor, e is the axial strain, y and z are the coordinates on the cross section, and ρi is the deformation curvature. u, v, and w are the displacements on the cross section.When a change of energy is applied, the product of the stiffness matrix and the displacement vector is expressed as follows:(37)kh=∫lBTΨTΦϕdx,where(38)B=∂N∂s,Ψ=1v′w′000000100002w′w″00−1−w′20−2v′v″001−v′20,Φ=EA0000GI110000EI220000EI33,ϕ=eρ1ρ2ρ3.Assume that(39)hi=h¯+Δhi.By applying Taylor expansion without higher order terms, the tangential stiffness matrix is expressed as follows:(40)k⌣=∫lBTΨ¯TΦΨ+ΓBdx,where(41)Γ=0000000EAe−2EI33v″ρ300−2EI33v′ρ3000EAe+2EI22w″ρ2002EI33w′ρ2000000000000000000.The mass matrix is the same as that of a linear beam element.


## 5.3. Nonlinear Convergence Algorithms

Although the Newton–Raphson algorithm is used for nonlinear static analysis, the tangent stiffness matrix becomes singular at some points, which makes a global equilibrium solution impossible. Riks proposed a process to track the intersection of the normal to the tangent line with the equilibrium path to solve this problem [29], as shown in Figure 11. The Crisfield method uses arcs instead of vertical lines to search for solutions. The increment of the load factor becomes an unknown problem to be solved in the iterative process.Figure 11
Scheme of the Crisfield method.Assume that(42)F=λF⇀,where λ is the load parameter and F⇀ is the preselected load vector.The residual force vector of thei‐th step is expressed as follows:(43)RU⇀i,λi=RU⇀i−1,λi−1+∂R∂U⇀U⇀i−1,λi−1δU⇀i+∂R∂λU⇀i−1,λi−1δλ1i,where U⇀ is the displacement vector.(44)kU⇀U⇀i−1,λi−1−λi−1F⇀+kTU⇀i−1,λi−1δU⇀1i+F⇀δλ1i=0,where δU⇀1i is the first increment of the load parameter in the i‐th step and δλ1i is the first increment of the displacement vector in the i‐th step.(45)U⇀1i=U⇀i−1+δU⇀1i,(46)λ1i=λi−1+δλ1i.Substituting (45) and (46) into equation (44), we obtain(47)kU⇀U⇀1i,λ1i−λ1iF⇀+kTU⇀1i,λ1iδU⇀2i+F⇀δλ2i=0.When the second step incremental search path is perpendicular to the normal of the first step incremental path, we obtain(48)δU⇀1i⋅δU⇀2i+F⇀δλ1i⋅F⇀δλ2i=0.When solving nonlinear finite element equations by an iterative method, choosing the appropriate convergence criterion is necessary to ensure that the iteration can be terminated. The convergence criterion will directly affect the speed and accuracy of the solution. If the convergence criterion is not appropriately chosen, the calculation will fail.From equations (47) and (48), we obtain(49)δU⇀2i=−kTU⇀1i,λ1i−1kU⇀U⇀1i,λ1i−λ1iF⇀+F⇀δλ2i,δλ2i=δU⇀1ikTU⇀1i,λ1i−1kU⇀U⇀1i,λ1i−λ1iF⇀F⇀δλ1iF⇀−δU⇀1ikTU⇀1i,λ1i−1F⇀.The convergent load parameters and displacement vectors are as follows:(50)λi=λi−1+δλ1i+δλ2i+δλ3i+⋯,U⇀i=U⇀i−1+δU⇀1i+δU⇀2i+δU⇀3i+⋯.


## 6. Experimental Analysis

### 6.1. Stress Calculation of the ILG

To verify the accuracy and practicability of the induced stress expression of the ILG, the deformation rate of each point in the uniaxial tensile test is substituted into the stress expression to calculate the stress.For the Arruda–Boyce model, the parametersG and λm are as follows:(51)G=4.10736,λm=4.62529.A comparison of the uniaxial tensile stress calculation data with the experimental data is given in Table1.Table 1
Comparison of the tensile stress calculation data with experimental data obtained in the uniaxial tensile test.
Measurement data
Experimental results
Arruda–Boyce model
Force (mN)
Deformation (mm)
Strain (%)
Stress (kPa)
Calculated stress (σE) (kPa)
Error (%)
4.97
2.0
5.7
0.50
0.54
7.42
8.21
3.5
10
0.86
0.92
6.64
14.88
7.0
20
1.70
1.78
4.70
20.52
10.5
30
2.54
2.64
4.08
25.35
14.0
40
3.38
3.53
4.38
29.47
17.5
50
4.21
4.43
5.21
33.06
21.0
60
5.04
5.34
6.02
36.22
24.5
70
5.86
6.26
6.69
38.94
28.0
80
6.68
7.16
7.31
41.48
31.5
90
7.50
8.06
7.40
43.58
35.0
100
8.30
8.94
7.70
45.50
38.5
110
9.10
9.80
7.65
47.13
42.0
120
9.88
10.63
7.64
48.67
45.5
130
10.65
11.44
7.26
49.98
49.0
140
11.41
12.21
6.93
51.03
52.5
150
12.15
12.96
6.71
51.95
56.0
160
12.86
13.69
6.39
52.75
59.5
170
13.56
14.38
6.02
53.40
63.0
180
14.24
15.05
5.66
53.90
66.5
190
14.89
15.69
5.36
54.32
70.0
200
15.52
16.30
5.01
54.65
73.5
210
16.12
16.88
4.65
54.78
77.0
220
16.70
17.44
4.50
54.86
80.5
230
17.25
17.98
4.31
54.94
84.0
240
17.78
18.50
3.97
54.84
87.5
250
18.28
18.99
3.88
54.72
91.0
260
18.74
19.46
3.72
54.47
94.5
270
19.18
19.91
3.73
54.01
98.0
280
19.57
20.34
4.06
53.61
101.5
290
19.93
20.75
4.22
53.23
105.0
300
20.24
21.15
4.29
52.48
108.5
310
20.50
21.53
5.05
51.83
112.0
320
20.73
21.89
5.59
51.02
115.5
330
20.91
22.23
6.42
50.29
119.0
340
21.04
22.57
7.09
49.21
122.5
350
21.12
22.88
8.51Table1, which presents comparisons of the tensile stress calculation data with experimental data, shows that the minimum relative error of the Arruda–Boyce model is 3.72 and the maximum relative error is 8.51. This result shows that the calculation accuracy of the Arruda–Boyce model is high and satisfies the requirement of material performance analysis in soft robot design.Figure12 shows that the first half of the stress-strain curve calculated by the Arruda–Boyce model almost coincides with the experimental curve, while for the second half, the relative error between the calculated and experimental stress-strain curves becomes increasingly larger but remains small. The above analyses indicate that the calculation formula of the ILG stress is feasible and reliable.Figure 12
Comparison of the tensile stress calculation data with experimental data obtained in the uniaxial tensile test.


### 6.2. Experiment on the Soft Actuator

The square wave of the input signal has a magnitude of 4 V, a period of 30 s, and a duty cycle of 50%. The actuator size is 30 mm × 5 mm × 0.5 mm, the free segment length is 25 mm, and the average thickness of the ILG layer is approximately 0.4 mm.The interval of the displacement reading is 3 s. Figure13 shows that the maximum amplitude of the positive axis of the actuator is 4.8 mm, and the maximum amplitude of the negative axis is 5 mm. The positive and negative displacements of the actuator are basically symmetrical and stable, and the attenuation of the long-term working drive performance is very small.Figure 13
Soft actuator motion screenshots when the input voltage is 4 V.In Figure14, the simulation results are analyzed: when the input voltage is 3.0 V, the maximum displacement of the soft actuator is 3.7 mm; when the input voltage is 3.5 V, the maximum displacement is 4.3 mm; and when the input voltage is 4.0 V, the maximum displacement is 5.0 mm. The displacement calculation results for the different input voltages are shown in Table 2.Figure 14
(a–c) Displacement calculation results for the input voltages of 3 V, 3.5 V, and 4 V, respectively.
(a)
(b)
(c)Table 2
Displacement calculation results for the different input voltages.
Force (V)
3.0
3.5
4.0
Deformation (mm)
3.7
4.3
5.0According to Figure15, the simulation value curve is almost consistent with the experimental value curve, and the change trends of the displacement curves from the simulation and experiment are almost identical. The relative error between the calculated and experimental displacement curves is small. Therefore, the calculation formula of the soft actuator is feasible and credible.Figure 15
Comparison of the displacement calculation data with experimental data.


### 6.3. Working Performance Test of the Soft Actuator

The long-term stable working performance is an important factor in ionic EAP soft actuators. Since the traditional IPMC actuator contains a water-based electrolyte, maintaining a stable driving performance after operation in air for dozens of seconds is difficult due to the evaporation of water and electrolysis [30]. The ILG soft robot does not contain a water-based electrolyte, which avoids the problem of water molecule loss, and can work stably in air for a long time.The current response stability of the soft actuator is tested by applying a square wave signal with a magnitude of 4 V, a period of 30 s, and a duty cycle of 50%, and the current response curve is shown in Figure16. The experimental results show that the current peaks remain unchanged over 30 minutes, and the current response of the ILG actuator is stable. Therefore, the internal electrolyte is not substantially decomposed.Figure 16
Current response curve of the soft actuator.From the above experiments, the ILG actuator can maintain stable current response and bending displacement in air. These results show that the ILG actuator can effectively avoid the defects of traditional EAP actuators and has excellent driving performance.


### 6.4. Experiment of the Soft Robot Grabbing an Object

First, an input voltage of 3.5 V is applied across the electrodes, and the soft robot gradually opens outward. Next, the operating platform is moved close to the target object. Finally, by changing the direction of the supply voltage, the soft robot reversely bends to clamp the object.The parameters of the soft robot are shown in Table3. The values 1∼3 in the first column are the ILG actuator number, corresponding the three sets of execution components on the manipulator, with a total mass of 211 mg.Table 3
Parameters of the soft manipulator.
Number
Total length (mm)
Mass (g)
Free segment length (mm)
Width (mm)
Total thickness (mm)
ILG layer thickness (mm)
1
35
0.068
30
5
0.5
0.4
2
35
0.070
30
5
0.5
0.4
3
35
0.073
30
5
0.5
0.4A set of screenshots of the soft robot manipulator’s experimental process is shown in Figure17. The process of the soft robot grabbing an object is shown in Figure 3. From Figures 17(a) and 17(b), when no object is grabbed, the deformation of the soft actuator is larger, and the deformation is similar to that in Figure 13. In this case, the mass of the object clamped by the soft robot at an input voltage of 3.5 V is approximately 105 mg. The experiments verify the effectiveness of ILG soft robots.Figure 17
Screenshots of the soft robot grabbing an object. (a) Natural stretching. (b) Bending inward. (c) Bending outward. (d) Grabbing an object.
(a)
(b)
(c)
(d)Compared with the soft robot designed by Saito et al., which has a load capacity of approximately 3 mg, the ILG soft robot designed in this paper has a greatly increased load capacity. Since the soft robot consists of three actuators, this ILG soft robot has better adaptability for grabbing complex shape objects than Saito’s soft robot, which consists of two actuators [31].In addition, the service life is an important performance indicator for EAP robots. The ILG material avoids the defects of IPMCs, and the soft robot proposed in this paper can theoretically work in air for a long time.


## 6.1. Stress Calculation of the ILG

To verify the accuracy and practicability of the induced stress expression of the ILG, the deformation rate of each point in the uniaxial tensile test is substituted into the stress expression to calculate the stress.For the Arruda–Boyce model, the parametersG and λm are as follows:(51)G=4.10736,λm=4.62529.A comparison of the uniaxial tensile stress calculation data with the experimental data is given in Table1.Table 1
Comparison of the tensile stress calculation data with experimental data obtained in the uniaxial tensile test.
Measurement data
Experimental results
Arruda–Boyce model
Force (mN)
Deformation (mm)
Strain (%)
Stress (kPa)
Calculated stress (σE) (kPa)
Error (%)
4.97
2.0
5.7
0.50
0.54
7.42
8.21
3.5
10
0.86
0.92
6.64
14.88
7.0
20
1.70
1.78
4.70
20.52
10.5
30
2.54
2.64
4.08
25.35
14.0
40
3.38
3.53
4.38
29.47
17.5
50
4.21
4.43
5.21
33.06
21.0
60
5.04
5.34
6.02
36.22
24.5
70
5.86
6.26
6.69
38.94
28.0
80
6.68
7.16
7.31
41.48
31.5
90
7.50
8.06
7.40
43.58
35.0
100
8.30
8.94
7.70
45.50
38.5
110
9.10
9.80
7.65
47.13
42.0
120
9.88
10.63
7.64
48.67
45.5
130
10.65
11.44
7.26
49.98
49.0
140
11.41
12.21
6.93
51.03
52.5
150
12.15
12.96
6.71
51.95
56.0
160
12.86
13.69
6.39
52.75
59.5
170
13.56
14.38
6.02
53.40
63.0
180
14.24
15.05
5.66
53.90
66.5
190
14.89
15.69
5.36
54.32
70.0
200
15.52
16.30
5.01
54.65
73.5
210
16.12
16.88
4.65
54.78
77.0
220
16.70
17.44
4.50
54.86
80.5
230
17.25
17.98
4.31
54.94
84.0
240
17.78
18.50
3.97
54.84
87.5
250
18.28
18.99
3.88
54.72
91.0
260
18.74
19.46
3.72
54.47
94.5
270
19.18
19.91
3.73
54.01
98.0
280
19.57
20.34
4.06
53.61
101.5
290
19.93
20.75
4.22
53.23
105.0
300
20.24
21.15
4.29
52.48
108.5
310
20.50
21.53
5.05
51.83
112.0
320
20.73
21.89
5.59
51.02
115.5
330
20.91
22.23
6.42
50.29
119.0
340
21.04
22.57
7.09
49.21
122.5
350
21.12
22.88
8.51Table1, which presents comparisons of the tensile stress calculation data with experimental data, shows that the minimum relative error of the Arruda–Boyce model is 3.72 and the maximum relative error is 8.51. This result shows that the calculation accuracy of the Arruda–Boyce model is high and satisfies the requirement of material performance analysis in soft robot design.Figure12 shows that the first half of the stress-strain curve calculated by the Arruda–Boyce model almost coincides with the experimental curve, while for the second half, the relative error between the calculated and experimental stress-strain curves becomes increasingly larger but remains small. The above analyses indicate that the calculation formula of the ILG stress is feasible and reliable.Figure 12
Comparison of the tensile stress calculation data with experimental data obtained in the uniaxial tensile test.


## 6.2. Experiment on the Soft Actuator

The square wave of the input signal has a magnitude of 4 V, a period of 30 s, and a duty cycle of 50%. The actuator size is 30 mm × 5 mm × 0.5 mm, the free segment length is 25 mm, and the average thickness of the ILG layer is approximately 0.4 mm.The interval of the displacement reading is 3 s. Figure13 shows that the maximum amplitude of the positive axis of the actuator is 4.8 mm, and the maximum amplitude of the negative axis is 5 mm. The positive and negative displacements of the actuator are basically symmetrical and stable, and the attenuation of the long-term working drive performance is very small.Figure 13
Soft actuator motion screenshots when the input voltage is 4 V.In Figure14, the simulation results are analyzed: when the input voltage is 3.0 V, the maximum displacement of the soft actuator is 3.7 mm; when the input voltage is 3.5 V, the maximum displacement is 4.3 mm; and when the input voltage is 4.0 V, the maximum displacement is 5.0 mm. The displacement calculation results for the different input voltages are shown in Table 2.Figure 14
(a–c) Displacement calculation results for the input voltages of 3 V, 3.5 V, and 4 V, respectively.
(a)
(b)
(c)Table 2
Displacement calculation results for the different input voltages.
Force (V)
3.0
3.5
4.0
Deformation (mm)
3.7
4.3
5.0According to Figure15, the simulation value curve is almost consistent with the experimental value curve, and the change trends of the displacement curves from the simulation and experiment are almost identical. The relative error between the calculated and experimental displacement curves is small. Therefore, the calculation formula of the soft actuator is feasible and credible.Figure 15
Comparison of the displacement calculation data with experimental data.


## 6.3. Working Performance Test of the Soft Actuator

The long-term stable working performance is an important factor in ionic EAP soft actuators. Since the traditional IPMC actuator contains a water-based electrolyte, maintaining a stable driving performance after operation in air for dozens of seconds is difficult due to the evaporation of water and electrolysis [30]. The ILG soft robot does not contain a water-based electrolyte, which avoids the problem of water molecule loss, and can work stably in air for a long time.The current response stability of the soft actuator is tested by applying a square wave signal with a magnitude of 4 V, a period of 30 s, and a duty cycle of 50%, and the current response curve is shown in Figure16. The experimental results show that the current peaks remain unchanged over 30 minutes, and the current response of the ILG actuator is stable. Therefore, the internal electrolyte is not substantially decomposed.Figure 16
Current response curve of the soft actuator.From the above experiments, the ILG actuator can maintain stable current response and bending displacement in air. These results show that the ILG actuator can effectively avoid the defects of traditional EAP actuators and has excellent driving performance.


## 6.4. Experiment of the Soft Robot Grabbing an Object

First, an input voltage of 3.5 V is applied across the electrodes, and the soft robot gradually opens outward. Next, the operating platform is moved close to the target object. Finally, by changing the direction of the supply voltage, the soft robot reversely bends to clamp the object.The parameters of the soft robot are shown in Table3. The values 1∼3 in the first column are the ILG actuator number, corresponding the three sets of execution components on the manipulator, with a total mass of 211 mg.Table 3
Parameters of the soft manipulator.
Number
Total length (mm)
Mass (g)
Free segment length (mm)
Width (mm)
Total thickness (mm)
ILG layer thickness (mm)
1
35
0.068
30
5
0.5
0.4
2
35
0.070
30
5
0.5
0.4
3
35
0.073
30
5
0.5
0.4A set of screenshots of the soft robot manipulator’s experimental process is shown in Figure17. The process of the soft robot grabbing an object is shown in Figure 3. From Figures 17(a) and 17(b), when no object is grabbed, the deformation of the soft actuator is larger, and the deformation is similar to that in Figure 13. In this case, the mass of the object clamped by the soft robot at an input voltage of 3.5 V is approximately 105 mg. The experiments verify the effectiveness of ILG soft robots.Figure 17
Screenshots of the soft robot grabbing an object. (a) Natural stretching. (b) Bending inward. (c) Bending outward. (d) Grabbing an object.
(a)
(b)
(c)
(d)Compared with the soft robot designed by Saito et al., which has a load capacity of approximately 3 mg, the ILG soft robot designed in this paper has a greatly increased load capacity. Since the soft robot consists of three actuators, this ILG soft robot has better adaptability for grabbing complex shape objects than Saito’s soft robot, which consists of two actuators [31].In addition, the service life is an important performance indicator for EAP robots. The ILG material avoids the defects of IPMCs, and the soft robot proposed in this paper can theoretically work in air for a long time.


## 7. Conclusions

A soft robot is a new flexible structure, and the materials and mechanism design methods of a soft robot are quite different from those of traditional robots. To avoid the shortcomings of traditional rigid manipulators, such as large energy consumption and complex mechanisms, a soft robot based on a new ILG material is proposed in this paper.A high-performance ILG actuator based on the principle of light curing was developed, and a modular design method based on motion and constraints was discussed. An ion transport model of porous media based on the Nernst–Planck equation was established, and the deformation mechanism of the ILG actuator was deeply analyzed in combination with the theory of cantilever beam deformation. Based on nonlinear finite element theory, the deformation of the soft robot when grabbing an object was analyzed. A three-finger type soft robot was designed, and the maximum displacement was 5 mm at a 4 V square wave voltage. The soft robot can grab a 105 mg object under a 3.5 V voltage. The results prove that the ILG bionic soft robot has good development prospects.



---

*Source: 2857282-2019-05-02.xml*