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

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.


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. ese factors limit the application of rigid robots in unstructured and complex environments [1][2][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][5][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.

Design of the ILG Soft Robot
e ILG soft actuator consists of a 5-layer structure, as shown in Figure 1. e middle layer is an electroactive layer composed of the ILG material, which functions to store ionic liquids. e 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. e soft robot consists of three ILG actuators, a common copper electrode, and three independent copper electrodes, as shown in Figure 2. e input unit consists of a common electrode and an independent electrode. e common electrode is a copper column, with the upper end fixed with a fixture, and the copper electrode is connected to a power source. e soft robot consists of three ILG actuators evenly distributed around the common electrode at 120°. Two copper electrodes sandwich one ILG actuator. Figure 3 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.

Preparation of the Soft Actuator.
e 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]. BMIMBF 4 was selected as the ionic liquid. e carriers of ionic liquids need to have interconnected porous structures, large specific surface area and porosity, and good mechanical strength and electrochemical stability. erefore, HEMA was selected to prepare the carrier. is 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. e mixed solution consists of BMIMBF 4 , HEMA, and ZrO 2 . Under UV light, the ILG is generated by the polymerization reaction induced by DEAP. e polymer matrix is cross-linked into a porous network structure. ZrO 2 enhances the mechanical strength of the gel, and as the ZrO 2 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 BMIMBF 4 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).

Soft Actuator Deformation Mechanism.
For the ILG actuator, the application of a voltage on the electrode can   2 Advances in Materials Science and Engineering cause bending deformation, which is the result of the coupling of the electric field, chemical field, flow field, and field force. e 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). erefore, BMIM + and BF 4 − 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]. e entire actuator therefore bends toward the anode side, as shown in Figure 5.
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]: where i is the ion type (i � 1 represents the cation, and i � 2 represents the anion), j i is the flow of ions in the PHEMA pores, d i 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 z i is the number of ionic charges. e 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. e continuous equation of concentration changing with time is where ∇j i is the flow change of the i-th ion in the PHEMA pores and t is the time. e internal electric field balance equation of the ILG actuator can be expressed as follows: where D is the electric displacement, E is the electric field, ρ is the net charge density, and k e is the effective dielectric constant of the ILG. e 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: where ε c is the induced strain, α(s) is the induced strain coupling coefficient, and s is the Laplace domain. erefore, the axial induced stress is where E t is the elastic modulus of the electrode layer and d(s) is the axial-induced stress coupling coefficient. e bending moment of the actuator is where b and b 1 are half the thickness of the actuator and ILG layer, respectively, and w is the width of the actuator. e 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. e cantilever beam flexibility matrix is as follows [14]:

Advances in Materials Science and Engineering
where R θx , R θy , and R θz represent rotational flexibility around the x, y, and z axes, respectively, and R x , R y , and R z 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.

Hyperelastic Arruda-Boyce Model
Material nonlinearity is caused by the nonlinearity of the constitutive relationship of the material. e 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 Figure 6.
e 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 function W, which is used to define the strain energy density in the material.
e Arruda-Boyce strain energy function W is as follows [15][16][17]: with c 1 � (1/2), c 2 � (1/20), c 3 � (11/1050), c 4 � (19/7000), and c 5 � (519/673750), where W is the strain energy density, G is the initial shear modulus, λ m is the deformation rate, J el is the elastic volume ratio, and D is the temperature-dependent material parameter associated with the bulk modulus. A material is completely incompressible when J el � 1. I 1 is the first strain invariant of the deviatoric strain, and the relationship between I 1 and the main tensile strain rates λ 1 , λ 2 , and λ 3 is as follows: where λ 1 , λ 2 , and λ 3 are the main (extension) deformation rates along the x, y, and z directions, respectively. I 1 , I 2 , and I 3 are the relative changes of the length, the surface area, and the volume of the elastomer, respectively. e material stress is obtained by differentiating the stretching: where σ n is the normal stress and λ n is the stretching in the loading direction. e isotropic and incompressible deformation process of an ILG is given by �� When the ILG is uniformly stretched, λ 2 � λ 3 ; equation (11) can be used to calculate e true stress is obtained by Differentiating equation (8), we obtain

Model for the Soft Robot.
To establish a gripping model scheme as shown in Figure 7, 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 k g is the curvature profile when the object is gripped [18,19].

Deformation Analysis for Soft Actuators.
At s � s * , the position of any material point can be represented by the following formula [20], as shown in Figure 8.
where, in Cartesian coordinates, X and Y are represented as follows: where η is a dummy variable and θ is the angle between the unit tangent vector and the horizontal direction.
where the prime (′) represents the partial derivative with respect to s. Assuming that r ⇀ is continuous, θ and r ′ ⇀ are also continuous.
At point s � ς, the jump in any function χ � χ(s, θ(s), θ ′ (s)) is expressed by a compact notation: where e jump will be associated with the force applied at the discrete point. e bending moment M is described by a constitutive equation: where k 0 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. e governing equation can be derived from the equilibrium of the linear and angular momentum [21][22][23]:

Deformation Analysis Governing Formula for the Soft
Robot. e corner of the object is assumed to be in contact with the soft robot at s � ξ. Two new unit vectors are defined as follows, as shown in Figure 9: 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.
where f 1 is the normal force and f 2 is the corresponding friction force. Equation (21) is applied to segment s ∈ (ξ, l), and we obtain the following equation: e expression of the potential energy W of the soft robot is established: where Equation (21) can be used to establish boundary value problems.  where the solution θ � θ * satisfies the conditions e soft robot is divided into K segments. e total potential energy expression is approximated as follows: where the constraint function W c is

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.

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 Figure 10. (x 1 , y 1 , z 1 ) and (x 2 , y 2 , z 2 ) are the position co-  Advances in Materials Science and Engineering ordinates of points P 1 and P 2 before deformation, respectively. e object is deformed to a new position under the action of external forces. (u 1 , v 1 , w 1 ) and (u 2 , v 2 , w 2 ) are the deformation coordinates of points P 1 and P 2 after deformation, respectively. e formula for the undeformed length of a truss element is as follows: e expression of the total Lagrangian strain ε along the deformation axis is According to the change of the elastic energy, the product of the stiffness matrix and the displacement vector is expressed as follows [2,24]: where E is Young's modulus, A is the cross sectional area, and l is the length of the truss element. e tangent stiffness matrix is obtained by differentiating equation (33) with respect to the displacement vector: e mass matrix is the same as that for a linear truss element [25,26]:

Nonlinear FEM for Beam Elements.
e two famous beam theories are the Euler-Bernoulli and Timoshenko beam theories. e 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. erefore, 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]: 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: X (x 1 , y 1 , z 1 ) (x 2 , y 2 , z 2 ) (x 1 + u 1 , y 1 + v 1 , z 1 + w 1 ) Figure 10: Nondeformed and deformed geometry of a truss element.
where [B] z [N] zs , Assume that By applying Taylor expansion without higher order terms, the tangential sti ness matrix is expressed as follows: where e mass matrix is the same as that of a linear beam element.

Nonlinear Convergence Algorithms.
Although the Newton-Raphson algorithm is used for nonlinear static analysis, the tangent sti ness 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. e Cris eld method uses arcs instead of vertical lines to search for solutions. e increment of the load factor becomes an unknown problem to be solved in the iterative process.

Assume that
where λ is the load parameter and F ⇀ is the preselected load vector. e residual force vector of the i-th step is expressed as follows: Figure 11: Scheme of the Cris eld method.

Advances in Materials Science and Engineering
where U ⇀ is the displacement vector.
where δ U ⇀ i 1 is the first increment of the load parameter in the i-th step and δλ i 1 is the first increment of the displacement vector in the i-th step.
Substituting (45) and (46) into equation (44), we obtain When the second step incremental search path is perpendicular to the normal of the first step incremental path, we obtain 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. e 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 e convergent load parameters and displacement vectors are as follows: (50)

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 parameters G and λ m are as follows: (51) A comparison of the uniaxial tensile stress calculation data with the experimental data is given in Table 1. Table 1, 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. is 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. Figure 12 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. e above analyses indicate that the calculation formula of the ILG stress is feasible and reliable.

Experiment on the Soft Actuator.
e square wave of the input signal has a magnitude of 4 V, a period of 30 s, and a duty cycle of 50%. e 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. e interval of the displacement reading is 3 s. Figure 13 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. e 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.
In Figure 14, 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. e displacement calculation results for the different input voltages are shown in Table 2.
According to Figure 15, 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. e relative error between the calculated and experimental displacement curves is small. erefore, the calculation formula of the soft actuator is feasible and credible.

Working Performance Test of the Soft Actuator.
e longterm 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]. e 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.  e 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 Figure 16.

Advances in Materials Science and Engineering
e experimental results show that the current peaks remain unchanged over 30 minutes, and the current response of the ILG actuator is stable. erefore, the internal electrolyte is not substantially decomposed.
From the above experiments, the ILG actuator can maintain stable current response and bending displacement in air. ese results show that the ILG actuator can effectively avoid the defects of traditional EAP actuators and has excellent driving performance.

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.    Table 3. e 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.
A set of screenshots of the soft robot manipulator's experimental process is shown in Figure 17. e 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. e experiments verify the effectiveness of ILG soft robots.
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. e ILG material avoids the defects of IPMCs, and the soft robot proposed in this paper can theoretically work in air for a long time.

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. e soft robot can grab a 105 mg object under a 3.5 V voltage. e results prove that the ILG bionic soft robot has good development prospects.