# A Computing Method to Determine the Performance of an Ionic Liquid Gel Soft Actuator

Authors: Bin He; Chenghong Zhang; Yanmin Zhou; Zhipeng Wang
Journal: Applied Bionics and Biomechanics
 (2018)

---

## Abstract

A new type of soft actuator material—an ionic liquid gel (ILG) that consists of BMIMBF4, HEMA, DEAP, and ZrO2—is polymerized into a gel state under ultraviolet (UV) light irradiation. In this paper, we first propose that the ILG conforms to the assumptions of hyperelastic theory and that the Mooney-Rivlin model can be used to study the properties of the ILG. Under the five-parameter and nine-parameter Mooney-Rivlin models, the formulas for the calculation of the uniaxial tensile stress, plane uniform tensile stress, and 3D directional stress are deduced. The five-parameter and nine-parameter Mooney-Rivlin models of the ILG with a ZrO2 content of 3 wt% were obtained by uniaxial tensile testing, and the parameters are denoted as c10, c01, c20, c11, and c02 and c10, c01, c20, c11, c02, c30, c21, c12, and c03, respectively. Through the analysis and comparison of the uniaxial tensile stress between the calculated and experimental data, the error between the stress data calculated from the five-parameter Mooney-Rivlin model and the experimental data is less than 0.51%, and the error between the stress data calculated from the nine-parameter Mooney-Rivlin model and the experimental data is no more than 8.87%. Hence, our work presents a feasible and credible formula for the calculation of the stress of the ILG. This work opens a new path to assess the performance of a soft actuator composed of an ILG and will contribute to the optimized design of soft robots.


---

## Body

## 1. Introduction

At room temperature, the ionogel is a polymerized gelatinous mixture from the ionic liquid and polymer matrix under UV irradiation [1, 2]. The high conductivity and stability of ionic liquids and the good mechanical properties of polymers make ionic liquid gels an ideal replacement for the traditional electroactive polymers [3]. Due to their high environmental adaptability and low-pressure impedance characteristics, soft robots have shown broad application potential in the fields of biology, medicine, agriculture, and so on. The adoption of electroactive polymer (EAP) materials for soft robots has become a hot research topic in recent years. Progress in electrochemical actuators has been made over the past few decades due to their desirable mechanical properties for intelligent robots, which are an alternative for air- and fluid-derived equipment [4–7]. Ionic liquid gels (ILG) are suitable building blocks for advanced actuators due to their tunable ionic conductivity, chemical stability, thermal stability, and simple ion transport [8-9].Electrochemical actuators have been further developed over the past few decades for their desirable mechanical properties in intelligent robots, which are alternatives to air- and fluid-derived devices. The flexible ionic conductivity of ILG is better suited to the evolution of building blocks due to the simpler ion transport actuators [10].Noncovalent interactions provide the gels with the very high mechanical strength and excellent self-healing ability of supramolecular materials [11, 12]. Based on these studies, we used ZrO2 to fabricate supramolecular nanocomposites with the electrochemical behavior of an ionic liquid and mechanical strength of an ionogel polymer.Among numerical approaches, finite element analysis is one of the most effective and most common information extraction methods to evaluate and optimize robot designs. By using this method, analytical models can be greatly simplified, which greatly increases the computational efficiency. The largest drawback is that ignorance of the nonlinear and constitutive model and the simplification of the computational model often lead to coarse solutions [13, 14]. Numerical simulations offer sufficient insight for each case during general soft robot design.Lee et al. used the finite element method (FEM) to successfully predict the mechanical behavior of an ionic polymer-metal composite (IPMC) actuator [15]. Wang et al. established a model based on the FEM to determine the electromechanical bending behavior of photocurable ionogel actuators (PIA) [16]. He et al. used the common large-scale finite element analysis software ANSYS to simulate an ILG, which is based on the SOLID186 element and the nonlinear hyperelastic Mooney-Rivlin model [17].Because the use of an ILG actuator requires deep understanding of the mechanical properties of the soft robot, it is necessary to establish a theoretical model to obtain its performance index. To reduce the experimental cost and time, we systematically analyzed the ILG via numerical simulations and verified the accuracy of the calculation, which contributed to the development of the ILG in the soft robot. The numerical simulation results matched the corresponding experiments, proving the validity of the model [18, 19].


## 2. Fabrication of the Ionic Liquid Gel

In the experiments, the ILG was composed of 1-butyl-3-methylimidazolium tetrafluoroborate (BMIMBF4), hydroxyethyl methacrylate (HEMA), 2-diethoxyacetophenone (DEAP), and ZrO2, with masses of 900 mg, 68.6 mg, 1.4 mg, and 30 mg, respectively. The mixed solution was then placed into a magnetic stirrer to form a suspension. Following that, the sample was placed on an ML-3500C Maxima-type cold light source for ultra-high-intensity UV radiation curing to induce polymerization. The UV intensity is 90000uW / cm2 (15″/ 380 mm distance) [20].Figure1(a) shows a comparison of the morphology of the ionic liquid before and after gel formation: the left figure shows the liquid state, while the right shows the solid state formed after polymerization. The schematic diagram of the principle behind gel formation is illustrated in Figures 1(b) and 1(c). Under UV light irradiation and the action of the DEAP catalyst, the polymer matrix was cross-linked into a porous network structure.Figure 1
The proposed mechanism of the BMIMBF4-based ionogel under ultraviolet light. (a) ILG solution and ILG solid, respectively. (b) The ILG solution includes BMIMBF4, HEMA, DEAP, and ZrO2. (c) HEMA and ZrO2 are cross-coupled to form a 3D network.


## 3. Material Nonlinearity and Parameters

The morphological analysis of the freeze-dried sample by scanning electron microscopy (SEM) showed that a porous microstructure was present throughout the ionogel. Distilled water was used to replace the internal ionic liquid in the ILG after freeze-drying treatment, and then, an S4800 Hitachi high-resolution field-emission scanning electron microscope was used to scan the sample. Figure2 shows the spatial structure of the ionic liquid carrier HEMA for the typical 3D porous structure with a 5000x magnification, in which its matrix is cross-linked to generate a 3D support skeleton, offering good mechanical strength and self-repairing performance. Due to the porous 3D network in the ionogel, gelated BMIMBF4 retains a relatively high ionic conductivity.Figure 2
SEM image of a freeze-dried BMIMBF4-based gel after replacing the ionic liquid with water [17].


### 3.1. Hyperelastic Hypothesis

Assuming that the ILG is an isotropic, incompressible, hyperelastic body, the following assumptions can be made based on the theory of continuum mechanics to study its mechanical properties.(1)
The strain energy functionW of a unit mass of material is an analytic function of the strain tensor of the natural state, termed the hyperelastic hypothesis. If the rate of change of W is equal to the power of the stress, then the material is a hyperelastic material. The mechanical properties of a hyperelastic material are described by the strain energy density function W, which has many functions.
(2)
Isotropy can be assumed.
(3)
The volume of the material before and after deformation can be assumed to be the same.λ1, λ2, and λ3 are set as the x, y, and z directions of the main (extension) deformation rate, respectively, given by
(1)λ1=xx0,λ2=yy0,λ3=dd0,where x, y, and d are the length, width, and thickness, respectively, and x0, y0, and d0 are the corresponding initial values before deformation.Because the material is incompressible, its volume is the same before and after deformation, giving(2)xyd=x0y0d0,namely,
(3)xydx0y0d0=λ1λ2λ3=1.


### 3.2. Hyperelastic Stress

Considering the mechanical performance requirements for potential applications, three kinds of deformation states should be examined: Figure3(a) shows the uniaxial tensile stress, Figure 3(b) shows uniform prestretching in the X and Y directions, and Figure 3(c) shows the Maxwell stress increased along the thickness.Figure 3
Analyses of the stress in different states.The physical properties are mainly expressed by the strain energy function, and each model is a special form of this function [21–23]. Once the form of the strain energy function W is determined, the Cauchy stress tensor P can be given by
(4)σ=−pI+2∂W∂I1B−2∂W∂I2B−1,where I is the unit tensor, which is the left Gauss deformation tensor, and p is the hydrostatic pressure resulting from the assumption of incompressibility.
(5)I1=B,I2=12I12−trB2,I3=detB,where B is the component of the Green strain tensor. The relationship between the invariants and the principal elongation is of a function of B.
(6)I1=trB=Bii=λ12+λ22+λ32,I2=12trB2−trB2=12BiiBii−BijBji=λ12λ22+λ22λ32+λ32λ12=1λ12+1λ22+1λ32,I3=detB=λ12λ22λ32.Hence, the isotropic and incompressible deformation process of an ILG is given as(7)I3=λ1λ2λ3=1.According to (4) and (6), we can obtain
(8)σi=2λi2∂W∂I1−1λi2∂W∂I2−p,where I1, I2, and I3 are the relative changes in the length, surface area, and volume of the elastomer, respectively.


### 3.3. Mooney-Rivlin Model

The Mooney-Rivlin model was chosen after comparing various hyperelastic constitutive models. The mechanical properties of ionic gel materials can be studied by using the Mooney-Rivlin formula, which is considered to be a nonlinear finite element of ionic gels in this study [24–26].The strain energy function in the Mooney-Rivlin model equation is as follows:(9)W=∑i≠jncijI1−3iI2−3j,where cij is a constant.The Mooney-Rivlin model is the most widely used strain energy function in the finite element method. It assumes that the strain energy density is a first-order function of the principal strain constant. Under large deformations, the mechanical properties of the ILG, as an incompressible hyperelastic material, are described.


### 3.4. Selection of the Constitutive Model for the ILG

Yeoh noted that the practical value of the higher-order strain energy function is small because the reproducibility of the ILG material is not sufficient and does not allow accurate estimation of a large number of parameters [27].Due to its simplicity and practicality, the Mooney-Rivlin model is widely used in finite element analysis. The first-order Mooney-Rivlin model describes the mechanical properties of incompressible hyperelastic materials under large deformation. The higher-order Mooney-Rivlin model can obtain a good approximation for the solution of large strain.The mechanical response of the hyperelastic material model is determined by the strain energy density function. The Mooney-Rivlin constant of the material must be accurately evaluated to obtain a reliable result from the hyperelastic analysis. In finite element analysis, the hyperelastic material is generally assumed to be a homogeneous isotropic material whose elastic modulus (Young’s modulus)E, initial shear modulus G0, and Poisson’s ratio v satisfy the following relationship [28, 29].
(10)E=2G01+v.


## 3.1. Hyperelastic Hypothesis

Assuming that the ILG is an isotropic, incompressible, hyperelastic body, the following assumptions can be made based on the theory of continuum mechanics to study its mechanical properties.(1)
The strain energy functionW of a unit mass of material is an analytic function of the strain tensor of the natural state, termed the hyperelastic hypothesis. If the rate of change of W is equal to the power of the stress, then the material is a hyperelastic material. The mechanical properties of a hyperelastic material are described by the strain energy density function W, which has many functions.
(2)
Isotropy can be assumed.
(3)
The volume of the material before and after deformation can be assumed to be the same.λ1, λ2, and λ3 are set as the x, y, and z directions of the main (extension) deformation rate, respectively, given by
(1)λ1=xx0,λ2=yy0,λ3=dd0,where x, y, and d are the length, width, and thickness, respectively, and x0, y0, and d0 are the corresponding initial values before deformation.Because the material is incompressible, its volume is the same before and after deformation, giving(2)xyd=x0y0d0,namely,
(3)xydx0y0d0=λ1λ2λ3=1.


## 3.2. Hyperelastic Stress

Considering the mechanical performance requirements for potential applications, three kinds of deformation states should be examined: Figure3(a) shows the uniaxial tensile stress, Figure 3(b) shows uniform prestretching in the X and Y directions, and Figure 3(c) shows the Maxwell stress increased along the thickness.Figure 3
Analyses of the stress in different states.The physical properties are mainly expressed by the strain energy function, and each model is a special form of this function [21–23]. Once the form of the strain energy function W is determined, the Cauchy stress tensor P can be given by
(4)σ=−pI+2∂W∂I1B−2∂W∂I2B−1,where I is the unit tensor, which is the left Gauss deformation tensor, and p is the hydrostatic pressure resulting from the assumption of incompressibility.
(5)I1=B,I2=12I12−trB2,I3=detB,where B is the component of the Green strain tensor. The relationship between the invariants and the principal elongation is of a function of B.
(6)I1=trB=Bii=λ12+λ22+λ32,I2=12trB2−trB2=12BiiBii−BijBji=λ12λ22+λ22λ32+λ32λ12=1λ12+1λ22+1λ32,I3=detB=λ12λ22λ32.Hence, the isotropic and incompressible deformation process of an ILG is given as(7)I3=λ1λ2λ3=1.According to (4) and (6), we can obtain
(8)σi=2λi2∂W∂I1−1λi2∂W∂I2−p,where I1, I2, and I3 are the relative changes in the length, surface area, and volume of the elastomer, respectively.


## 3.3. Mooney-Rivlin Model

The Mooney-Rivlin model was chosen after comparing various hyperelastic constitutive models. The mechanical properties of ionic gel materials can be studied by using the Mooney-Rivlin formula, which is considered to be a nonlinear finite element of ionic gels in this study [24–26].The strain energy function in the Mooney-Rivlin model equation is as follows:(9)W=∑i≠jncijI1−3iI2−3j,where cij is a constant.The Mooney-Rivlin model is the most widely used strain energy function in the finite element method. It assumes that the strain energy density is a first-order function of the principal strain constant. Under large deformations, the mechanical properties of the ILG, as an incompressible hyperelastic material, are described.


## 3.4. Selection of the Constitutive Model for the ILG

Yeoh noted that the practical value of the higher-order strain energy function is small because the reproducibility of the ILG material is not sufficient and does not allow accurate estimation of a large number of parameters [27].Due to its simplicity and practicality, the Mooney-Rivlin model is widely used in finite element analysis. The first-order Mooney-Rivlin model describes the mechanical properties of incompressible hyperelastic materials under large deformation. The higher-order Mooney-Rivlin model can obtain a good approximation for the solution of large strain.The mechanical response of the hyperelastic material model is determined by the strain energy density function. The Mooney-Rivlin constant of the material must be accurately evaluated to obtain a reliable result from the hyperelastic analysis. In finite element analysis, the hyperelastic material is generally assumed to be a homogeneous isotropic material whose elastic modulus (Young’s modulus)E, initial shear modulus G0, and Poisson’s ratio v satisfy the following relationship [28, 29].
(10)E=2G01+v.


## 4. Stress Calculation

### 4.1. Uniaxial Tension (State I)

#### 4.1.1. Five-Parameter Mooney-Rivlin Model

From (9), the strain energy equation in the five-parameter Mooney-Rivlin model can be written as
(11)W=c10I1−3+c01I2−3+c20I1−32+c11I1−3I2−3+c02I2−32.From (8), we get
(12)σi=2λi2∂W∂I1−1λi2∂W∂I2−p.The stresses in all directions are(13)σ1=2λ12∂W∂I1−1λ12∂W∂I2−p,(14)σ2=2λ22∂W∂I1−1λ22∂W∂I2−p,(15)σ3=2λ32∂W∂I1−1λ32∂W∂I2−p,where ∂W/∂I1 and ∂W/∂I2 are the partial differentials in the strain energy function W for I1 and I2, respectively, and σ1, σ2, and σ3 are the stresses in the x, y, and z directions, respectively.When uniformly stretched in theX direction, λ2=λ3, and from (3), we get
(16)λ2=λ3=1λ1.From (6) and (16), we obtain
(17)I1=λ12+2λ1,I2=1λ12+2λ1.Because only axial tensile deformation is considered, the stress in the other two directions is zero.(18)σ2=σ3=0.From (14) or (15) and (16), we obtain
(19)p=21λ1∂W∂I1−λ1∂W∂I2.Substituting (19) into (13),
(20)σ1=2λ12−1λ1∂W∂I1+1λ1∂W∂I2.For (11), the partial differentials of the strain energy function W for I1 and I2 are given by
(21)∂W∂I1=c10+2c02I1−3+c11I2−3,∂W∂I2=c01+2c02I2−3+c11I1−3.Therefore,(22)σ1=2λ12−1λ1c10+2c02I1−3+c11I2−3+1λ1c01+2c02I2−3+c11I1−3.


#### 4.1.2. Nine-Parameter Mooney-Rivlin Model

From (9), the strain energy equation in the nine-parameter Mooney-Rivlin model is
(23)W=c10I1−3+c01I2−3+c20I1−32+c11I1−3I2−3+c02I2−32+c30I1−33+c21I1−32I2−3+c12I1−3I2−32+c03I2−33.For (23), the partial differentials of the strain energy function W for I1 and I2 are given by
(24)∂W∂I1=c10+2c20I1−3+c11I2−3+3c30I1−32+2c21I1−3I2−3+c12I2−32=M,∂W∂I2=c01+2c02I2−3+c11I1−3+3c03I2−32+2c12I1−3I2−3+c21I1−32=N.Therefore,(25)σ1=2λ12−1λ1M+1λ1N.


### 4.2. Evenly Stretched in theX and Y Directions (State II)

#### 4.2.1. Five-Parameter Mooney-Rivlin Model

From (9), the strain energy equation in the five-parameter Mooney-Rivlin model is
(26)W=c10I1−3+c01I2−3+c20I1−32+c11I1−3I2−3+c02I2−32.From (8), we get
(27)σi=2λi2∂W∂I1−1λi2∂W∂I2−p.The stresses in all directions are(28)σ1=2λ12∂W∂I1−1λ12∂W∂I2−p,(29)σ2=2λ22∂W∂I1−1λ22∂W∂I2−p,(30)σ3=2λ32∂W∂I1−1λ32∂W∂I2−p.When uniformly stretched in theX and Y directions, λ1=λ2, and from (3), we obtain
(31)λ1=λ2=1λ3.From (6) and (31), we get
(32)I1=2λ12+2λ14,I2=2λ12+λ14.Because tensile deformations are only considered in theX and Y directions, the stress in the thickness direction is zero.
(33)σ1=σ2,σ3=0.From (28) or (29) and (31), we obtain
(34)p=21λ14∂W∂I1−λ14∂W∂I2.Substituting (34) into (28) or (29),
(35)σ1=σ2=2λ12−1λ14∂W∂I1+λ12∂W∂I2.For (26), the partial differentials of the strain energy function W for I1 and I2 are given by
(36)∂W∂I1=c10+2c20I1−3+c11I2−3,∂W∂I2=c01+2c02I2−3+c11I1−3.Therefore,(37)σ1=σ2=2λ12−1λ14c10+2c02I1−3+c11I2−3+λ12c01+2c02I2−3+c11I1−3.


#### 4.2.2. Nine-Parameter Mooney-Rivlin Model

From (9), the strain energy equation of the nine-parameter Mooney-Rivlin model is
(38)W=c10I1−3+c01I2−3+c20I1−32+c11I1−3I2−3+c02I2−32+c30I1−33+c21I1−32I2−3+c12I1−3I2−32+c03I2−33.For (38), the partial differentials of the strain energy function W for I1 and I2 are given by
(39)∂W∂I1=c10+2c20I1−3+c11I2−3+3c30I1−32+2c21I1−3I2−3+c12I2−32=M.∂W∂I2=c01+2c02I2−3+c11I1−3+3c03I2−32+2c12I1−3I2−3+c21I1−32=N.Therefore,(40)σ1=σ2=2λ12−1λ14M+λ12N.


### 4.3. Applied Maxwell Stress Loading (State III)

The Maxwell stress can be written as(41)p=−ςς0V2=−ςς0u2d2,where ς is the insulation constant, ς0 is the vacuum dielectric constant (8.85 × 10−12F/m), V is the electric field intensity, and u is the voltage.


#### 4.3.1. Five-Parameter Mooney-Rivlin Model

From (9), the strain energy equation in the five-parameter Mooney-Rivlin model is
(42)W=c10I1−3+c01I2−3+c20I1−32+c11I1−3I2−3+c02I2−32.From (8), we get
(43)σi=2λi2∂W∂I1−1λi2∂W∂I2−p.The stresses in all directions are(44)σ1=2λ12∂W∂I1−1λ12∂W∂I2−p,(45)σ2=2λ22∂W∂I1−1λ22∂W∂I2−p,(46)σ3=2λ32∂W∂I1−1λ32∂W∂I2−p.When the Maxwell stress is loaded, namely,λ1=λ2, (3) gives
(47)λ1=λ2=1λ3.From (6) and (45), we get
(48)I1=2λ12+2λ14,I2=2λ12+λ14.For (42), the partial differentials of the strain energy function W for I1 and I2 are given by
(49)∂W∂I1=c10+2c20I1−3+c11I2−3,∂W∂I2=c01+2c02I2−3+c11I1−3.Because the materials uniformly stretched in theX and Y directions,
(50)σ1=σ2.From (41), we obtain
(51)σ3=−ςς0V2=−ςς0u2d2=−ςς0u2λ3d02.Substituting (49), (50), and (51) into (44) and (46),
(52)σ1=σ2=2λ12c10+2c20I1−3+c11I2−3−1λ12c01+2c02I2−3+c11I1−3−p,(53)σ3=21λ14c10+2c20I1−3+c11I2−3−λ14c01+2c02I2−3+c11I1−3−p,(54)σ3=−ςς0u2λ3d02=−ςς0λ14u2d02.From (53) and (54),
(55)p=21λ14c10+2c20I1−3+c11I2−3−λ14c01+2c02I2−3+c11I1−3+ςς0λ14u2d02.Therefore,(56)σ1=σ2=2λ12c10+2c20I1−3+c11I2−3−1λ12c01+2c02I2−3+c11I1−3−21λ14c10+2c20I1−3+c11I2−3−λ14c01+2c02I2−3+c11I1−3−ςς0λ14u2d02,σ3=−ςς0λ14u2d02.


#### 4.3.2. Nine-Parameter Mooney-Rivlin Model

From (9), the strain energy equation in the nine-parameter Mooney-Rivlin model is
(57)W=c10I1−3+c01I2−3+c20I1−32+c11I1−3I2−3+c02I2−32+c30I1−33+c21I1−32I2−3+c12I1−3I2−32+c03I2−33.For (57), the partial differentials of the strain energy function W for I1 and I2 are given by
(58)∂W∂I1=c10+2c20I1−3+c11I2−3+3c30I1−32+2c21I1−3I2−3+c12I2−32=M,∂W∂I2=c01+2c02I2−3+c11I1−3+3c03I2−32+2c12I1−3I2−3+c21I1−32=N.Because the materials uniformly stretched in theX and Y directions,
(59)σ1=σ2.From (41),
(60)σ3=−ςς0V2=−ςς0u2d2=−ςς0u2λ3d02.Substituting (58), (59), and (60) into (44) and (46),
(61)σ1=σ2=2λ12M−1λ12N−p,(62)σ3=21λ14M−λ14N−p=−ςς0u2λ3d02=−ςς0λ14u2d02.From (61) and (62),
(63)p=21λ14M−λ14N+ςς0λ14u2d02.Therefore,(64)σ1=σ2=2λ12M−1λ12N−21λ14M−λ14N−ςς0λ14u2d02,σ3=−ςς0λ14u2d02.


## 4.1. Uniaxial Tension (State I)

### 4.1.1. Five-Parameter Mooney-Rivlin Model

From (9), the strain energy equation in the five-parameter Mooney-Rivlin model can be written as
(11)W=c10I1−3+c01I2−3+c20I1−32+c11I1−3I2−3+c02I2−32.From (8), we get
(12)σi=2λi2∂W∂I1−1λi2∂W∂I2−p.The stresses in all directions are(13)σ1=2λ12∂W∂I1−1λ12∂W∂I2−p,(14)σ2=2λ22∂W∂I1−1λ22∂W∂I2−p,(15)σ3=2λ32∂W∂I1−1λ32∂W∂I2−p,where ∂W/∂I1 and ∂W/∂I2 are the partial differentials in the strain energy function W for I1 and I2, respectively, and σ1, σ2, and σ3 are the stresses in the x, y, and z directions, respectively.When uniformly stretched in theX direction, λ2=λ3, and from (3), we get
(16)λ2=λ3=1λ1.From (6) and (16), we obtain
(17)I1=λ12+2λ1,I2=1λ12+2λ1.Because only axial tensile deformation is considered, the stress in the other two directions is zero.(18)σ2=σ3=0.From (14) or (15) and (16), we obtain
(19)p=21λ1∂W∂I1−λ1∂W∂I2.Substituting (19) into (13),
(20)σ1=2λ12−1λ1∂W∂I1+1λ1∂W∂I2.For (11), the partial differentials of the strain energy function W for I1 and I2 are given by
(21)∂W∂I1=c10+2c02I1−3+c11I2−3,∂W∂I2=c01+2c02I2−3+c11I1−3.Therefore,(22)σ1=2λ12−1λ1c10+2c02I1−3+c11I2−3+1λ1c01+2c02I2−3+c11I1−3.


### 4.1.2. Nine-Parameter Mooney-Rivlin Model

From (9), the strain energy equation in the nine-parameter Mooney-Rivlin model is
(23)W=c10I1−3+c01I2−3+c20I1−32+c11I1−3I2−3+c02I2−32+c30I1−33+c21I1−32I2−3+c12I1−3I2−32+c03I2−33.For (23), the partial differentials of the strain energy function W for I1 and I2 are given by
(24)∂W∂I1=c10+2c20I1−3+c11I2−3+3c30I1−32+2c21I1−3I2−3+c12I2−32=M,∂W∂I2=c01+2c02I2−3+c11I1−3+3c03I2−32+2c12I1−3I2−3+c21I1−32=N.Therefore,(25)σ1=2λ12−1λ1M+1λ1N.


## 4.1.1. Five-Parameter Mooney-Rivlin Model

From (9), the strain energy equation in the five-parameter Mooney-Rivlin model can be written as
(11)W=c10I1−3+c01I2−3+c20I1−32+c11I1−3I2−3+c02I2−32.From (8), we get
(12)σi=2λi2∂W∂I1−1λi2∂W∂I2−p.The stresses in all directions are(13)σ1=2λ12∂W∂I1−1λ12∂W∂I2−p,(14)σ2=2λ22∂W∂I1−1λ22∂W∂I2−p,(15)σ3=2λ32∂W∂I1−1λ32∂W∂I2−p,where ∂W/∂I1 and ∂W/∂I2 are the partial differentials in the strain energy function W for I1 and I2, respectively, and σ1, σ2, and σ3 are the stresses in the x, y, and z directions, respectively.When uniformly stretched in theX direction, λ2=λ3, and from (3), we get
(16)λ2=λ3=1λ1.From (6) and (16), we obtain
(17)I1=λ12+2λ1,I2=1λ12+2λ1.Because only axial tensile deformation is considered, the stress in the other two directions is zero.(18)σ2=σ3=0.From (14) or (15) and (16), we obtain
(19)p=21λ1∂W∂I1−λ1∂W∂I2.Substituting (19) into (13),
(20)σ1=2λ12−1λ1∂W∂I1+1λ1∂W∂I2.For (11), the partial differentials of the strain energy function W for I1 and I2 are given by
(21)∂W∂I1=c10+2c02I1−3+c11I2−3,∂W∂I2=c01+2c02I2−3+c11I1−3.Therefore,(22)σ1=2λ12−1λ1c10+2c02I1−3+c11I2−3+1λ1c01+2c02I2−3+c11I1−3.


## 4.1.2. Nine-Parameter Mooney-Rivlin Model

From (9), the strain energy equation in the nine-parameter Mooney-Rivlin model is
(23)W=c10I1−3+c01I2−3+c20I1−32+c11I1−3I2−3+c02I2−32+c30I1−33+c21I1−32I2−3+c12I1−3I2−32+c03I2−33.For (23), the partial differentials of the strain energy function W for I1 and I2 are given by
(24)∂W∂I1=c10+2c20I1−3+c11I2−3+3c30I1−32+2c21I1−3I2−3+c12I2−32=M,∂W∂I2=c01+2c02I2−3+c11I1−3+3c03I2−32+2c12I1−3I2−3+c21I1−32=N.Therefore,(25)σ1=2λ12−1λ1M+1λ1N.


## 4.2. Evenly Stretched in theX and Y Directions (State II)

### 4.2.1. Five-Parameter Mooney-Rivlin Model

From (9), the strain energy equation in the five-parameter Mooney-Rivlin model is
(26)W=c10I1−3+c01I2−3+c20I1−32+c11I1−3I2−3+c02I2−32.From (8), we get
(27)σi=2λi2∂W∂I1−1λi2∂W∂I2−p.The stresses in all directions are(28)σ1=2λ12∂W∂I1−1λ12∂W∂I2−p,(29)σ2=2λ22∂W∂I1−1λ22∂W∂I2−p,(30)σ3=2λ32∂W∂I1−1λ32∂W∂I2−p.When uniformly stretched in theX and Y directions, λ1=λ2, and from (3), we obtain
(31)λ1=λ2=1λ3.From (6) and (31), we get
(32)I1=2λ12+2λ14,I2=2λ12+λ14.Because tensile deformations are only considered in theX and Y directions, the stress in the thickness direction is zero.
(33)σ1=σ2,σ3=0.From (28) or (29) and (31), we obtain
(34)p=21λ14∂W∂I1−λ14∂W∂I2.Substituting (34) into (28) or (29),
(35)σ1=σ2=2λ12−1λ14∂W∂I1+λ12∂W∂I2.For (26), the partial differentials of the strain energy function W for I1 and I2 are given by
(36)∂W∂I1=c10+2c20I1−3+c11I2−3,∂W∂I2=c01+2c02I2−3+c11I1−3.Therefore,(37)σ1=σ2=2λ12−1λ14c10+2c02I1−3+c11I2−3+λ12c01+2c02I2−3+c11I1−3.


### 4.2.2. Nine-Parameter Mooney-Rivlin Model

From (9), the strain energy equation of the nine-parameter Mooney-Rivlin model is
(38)W=c10I1−3+c01I2−3+c20I1−32+c11I1−3I2−3+c02I2−32+c30I1−33+c21I1−32I2−3+c12I1−3I2−32+c03I2−33.For (38), the partial differentials of the strain energy function W for I1 and I2 are given by
(39)∂W∂I1=c10+2c20I1−3+c11I2−3+3c30I1−32+2c21I1−3I2−3+c12I2−32=M.∂W∂I2=c01+2c02I2−3+c11I1−3+3c03I2−32+2c12I1−3I2−3+c21I1−32=N.Therefore,(40)σ1=σ2=2λ12−1λ14M+λ12N.


## 4.2.1. Five-Parameter Mooney-Rivlin Model

From (9), the strain energy equation in the five-parameter Mooney-Rivlin model is
(26)W=c10I1−3+c01I2−3+c20I1−32+c11I1−3I2−3+c02I2−32.From (8), we get
(27)σi=2λi2∂W∂I1−1λi2∂W∂I2−p.The stresses in all directions are(28)σ1=2λ12∂W∂I1−1λ12∂W∂I2−p,(29)σ2=2λ22∂W∂I1−1λ22∂W∂I2−p,(30)σ3=2λ32∂W∂I1−1λ32∂W∂I2−p.When uniformly stretched in theX and Y directions, λ1=λ2, and from (3), we obtain
(31)λ1=λ2=1λ3.From (6) and (31), we get
(32)I1=2λ12+2λ14,I2=2λ12+λ14.Because tensile deformations are only considered in theX and Y directions, the stress in the thickness direction is zero.
(33)σ1=σ2,σ3=0.From (28) or (29) and (31), we obtain
(34)p=21λ14∂W∂I1−λ14∂W∂I2.Substituting (34) into (28) or (29),
(35)σ1=σ2=2λ12−1λ14∂W∂I1+λ12∂W∂I2.For (26), the partial differentials of the strain energy function W for I1 and I2 are given by
(36)∂W∂I1=c10+2c20I1−3+c11I2−3,∂W∂I2=c01+2c02I2−3+c11I1−3.Therefore,(37)σ1=σ2=2λ12−1λ14c10+2c02I1−3+c11I2−3+λ12c01+2c02I2−3+c11I1−3.


## 4.2.2. Nine-Parameter Mooney-Rivlin Model

From (9), the strain energy equation of the nine-parameter Mooney-Rivlin model is
(38)W=c10I1−3+c01I2−3+c20I1−32+c11I1−3I2−3+c02I2−32+c30I1−33+c21I1−32I2−3+c12I1−3I2−32+c03I2−33.For (38), the partial differentials of the strain energy function W for I1 and I2 are given by
(39)∂W∂I1=c10+2c20I1−3+c11I2−3+3c30I1−32+2c21I1−3I2−3+c12I2−32=M.∂W∂I2=c01+2c02I2−3+c11I1−3+3c03I2−32+2c12I1−3I2−3+c21I1−32=N.Therefore,(40)σ1=σ2=2λ12−1λ14M+λ12N.


## 4.3. Applied Maxwell Stress Loading (State III)

The Maxwell stress can be written as(41)p=−ςς0V2=−ςς0u2d2,where ς is the insulation constant, ς0 is the vacuum dielectric constant (8.85 × 10−12F/m), V is the electric field intensity, and u is the voltage.


### 4.3.1. Five-Parameter Mooney-Rivlin Model

From (9), the strain energy equation in the five-parameter Mooney-Rivlin model is
(42)W=c10I1−3+c01I2−3+c20I1−32+c11I1−3I2−3+c02I2−32.From (8), we get
(43)σi=2λi2∂W∂I1−1λi2∂W∂I2−p.The stresses in all directions are(44)σ1=2λ12∂W∂I1−1λ12∂W∂I2−p,(45)σ2=2λ22∂W∂I1−1λ22∂W∂I2−p,(46)σ3=2λ32∂W∂I1−1λ32∂W∂I2−p.When the Maxwell stress is loaded, namely,λ1=λ2, (3) gives
(47)λ1=λ2=1λ3.From (6) and (45), we get
(48)I1=2λ12+2λ14,I2=2λ12+λ14.For (42), the partial differentials of the strain energy function W for I1 and I2 are given by
(49)∂W∂I1=c10+2c20I1−3+c11I2−3,∂W∂I2=c01+2c02I2−3+c11I1−3.Because the materials uniformly stretched in theX and Y directions,
(50)σ1=σ2.From (41), we obtain
(51)σ3=−ςς0V2=−ςς0u2d2=−ςς0u2λ3d02.Substituting (49), (50), and (51) into (44) and (46),
(52)σ1=σ2=2λ12c10+2c20I1−3+c11I2−3−1λ12c01+2c02I2−3+c11I1−3−p,(53)σ3=21λ14c10+2c20I1−3+c11I2−3−λ14c01+2c02I2−3+c11I1−3−p,(54)σ3=−ςς0u2λ3d02=−ςς0λ14u2d02.From (53) and (54),
(55)p=21λ14c10+2c20I1−3+c11I2−3−λ14c01+2c02I2−3+c11I1−3+ςς0λ14u2d02.Therefore,(56)σ1=σ2=2λ12c10+2c20I1−3+c11I2−3−1λ12c01+2c02I2−3+c11I1−3−21λ14c10+2c20I1−3+c11I2−3−λ14c01+2c02I2−3+c11I1−3−ςς0λ14u2d02,σ3=−ςς0λ14u2d02.


### 4.3.2. Nine-Parameter Mooney-Rivlin Model

From (9), the strain energy equation in the nine-parameter Mooney-Rivlin model is
(57)W=c10I1−3+c01I2−3+c20I1−32+c11I1−3I2−3+c02I2−32+c30I1−33+c21I1−32I2−3+c12I1−3I2−32+c03I2−33.For (57), the partial differentials of the strain energy function W for I1 and I2 are given by
(58)∂W∂I1=c10+2c20I1−3+c11I2−3+3c30I1−32+2c21I1−3I2−3+c12I2−32=M,∂W∂I2=c01+2c02I2−3+c11I1−3+3c03I2−32+2c12I1−3I2−3+c21I1−32=N.Because the materials uniformly stretched in theX and Y directions,
(59)σ1=σ2.From (41),
(60)σ3=−ςς0V2=−ςς0u2d2=−ςς0u2λ3d02.Substituting (58), (59), and (60) into (44) and (46),
(61)σ1=σ2=2λ12M−1λ12N−p,(62)σ3=21λ14M−λ14N−p=−ςς0u2λ3d02=−ςς0λ14u2d02.From (61) and (62),
(63)p=21λ14M−λ14N+ςς0λ14u2d02.Therefore,(64)σ1=σ2=2λ12M−1λ12N−21λ14M−λ14N−ςς0λ14u2d02,σ3=−ςς0λ14u2d02.


## 4.3.1. Five-Parameter Mooney-Rivlin Model

From (9), the strain energy equation in the five-parameter Mooney-Rivlin model is
(42)W=c10I1−3+c01I2−3+c20I1−32+c11I1−3I2−3+c02I2−32.From (8), we get
(43)σi=2λi2∂W∂I1−1λi2∂W∂I2−p.The stresses in all directions are(44)σ1=2λ12∂W∂I1−1λ12∂W∂I2−p,(45)σ2=2λ22∂W∂I1−1λ22∂W∂I2−p,(46)σ3=2λ32∂W∂I1−1λ32∂W∂I2−p.When the Maxwell stress is loaded, namely,λ1=λ2, (3) gives
(47)λ1=λ2=1λ3.From (6) and (45), we get
(48)I1=2λ12+2λ14,I2=2λ12+λ14.For (42), the partial differentials of the strain energy function W for I1 and I2 are given by
(49)∂W∂I1=c10+2c20I1−3+c11I2−3,∂W∂I2=c01+2c02I2−3+c11I1−3.Because the materials uniformly stretched in theX and Y directions,
(50)σ1=σ2.From (41), we obtain
(51)σ3=−ςς0V2=−ςς0u2d2=−ςς0u2λ3d02.Substituting (49), (50), and (51) into (44) and (46),
(52)σ1=σ2=2λ12c10+2c20I1−3+c11I2−3−1λ12c01+2c02I2−3+c11I1−3−p,(53)σ3=21λ14c10+2c20I1−3+c11I2−3−λ14c01+2c02I2−3+c11I1−3−p,(54)σ3=−ςς0u2λ3d02=−ςς0λ14u2d02.From (53) and (54),
(55)p=21λ14c10+2c20I1−3+c11I2−3−λ14c01+2c02I2−3+c11I1−3+ςς0λ14u2d02.Therefore,(56)σ1=σ2=2λ12c10+2c20I1−3+c11I2−3−1λ12c01+2c02I2−3+c11I1−3−21λ14c10+2c20I1−3+c11I2−3−λ14c01+2c02I2−3+c11I1−3−ςς0λ14u2d02,σ3=−ςς0λ14u2d02.


## 4.3.2. Nine-Parameter Mooney-Rivlin Model

From (9), the strain energy equation in the nine-parameter Mooney-Rivlin model is
(57)W=c10I1−3+c01I2−3+c20I1−32+c11I1−3I2−3+c02I2−32+c30I1−33+c21I1−32I2−3+c12I1−3I2−32+c03I2−33.For (57), the partial differentials of the strain energy function W for I1 and I2 are given by
(58)∂W∂I1=c10+2c20I1−3+c11I2−3+3c30I1−32+2c21I1−3I2−3+c12I2−32=M,∂W∂I2=c01+2c02I2−3+c11I1−3+3c03I2−32+2c12I1−3I2−3+c21I1−32=N.Because the materials uniformly stretched in theX and Y directions,
(59)σ1=σ2.From (41),
(60)σ3=−ςς0V2=−ςς0u2d2=−ςς0u2λ3d02.Substituting (58), (59), and (60) into (44) and (46),
(61)σ1=σ2=2λ12M−1λ12N−p,(62)σ3=21λ14M−λ14N−p=−ςς0u2λ3d02=−ςς0λ14u2d02.From (61) and (62),
(63)p=21λ14M−λ14N+ςς0λ14u2d02.Therefore,(64)σ1=σ2=2λ12M−1λ12N−21λ14M−λ14N−ςς0λ14u2d02,σ3=−ςς0λ14u2d02.


## 5. Experimental Analysis

### 5.1. Production of Tensile Test Sample

The prepared solution was poured into a mold and polymerized into a gel under the irradiation of a UV lamp. The gel is white in color and is filled in the transparent glass mold as shown in Figure4. Then, the sample was cut to the size shown in Figure 5, with a thickness of 3.4 mm. As shown in Figure 6, the experimental instrument was a UTM2502 electronic universal testing machine. The mechanical sensor on the testing machine can achieve a precision of 0.1 mN, the displacement sensor has a precision of 0.001 mm, and the stretching rate is 500 mm/min.Figure 4
The gel after polymerization under a UV lamp.Figure 5
Test sample size (in mm).Figure 6
Uniaxial tensile test.


### 5.2. Experimental Results

As seen from Figure6, the ILG becomes longer and thinner as the load increases, which is consistent with the assumption that the material is incompressible.The tensile stress-strain curve of the ILG is shown in Figure7, and the average tensile strength (Young’s modulus) of the material obtained from the tensile stress-strain curve is 7.6 kPa.Figure 7
Stress-strain curve of the ionic liquid gel polymer.In the ionogel, BMIMBF4 exhibits a high level of hyperelastic toughness when the tensile deformation reaches 360%. The tensile test showed that the tensile properties of the gel increased with an increase in ZrO2 content. The increase in ZrO2 content should generate more cross-linking sites and higher conversion rates (shown in Figure 2), which contributes to the overall mechanical properties. As the content of ZrO2 increases, the tensile strength of the ILG increases, while the elongation rate decreases. Considering the above data, the optimum amount of ZrO2 to provide a large tensile strain and tensile strength is 3 wt%.


### 5.3. Stress Calculation

The relationship between the real stressσi and engineering stress σE is
(65)σi=σEλ,where σi is the true stress and σE is the engineering stress. The calculated stress given in Table 1 is the engineering stress.Table 1
Comparison of the calculated tensile stress with the experimental data from the uniaxial tensile test.
Measurement data
Experimental results
Five-parameter Mooney-Rivlin model
Nine-parameter Mooney-Rivlin model
Force(mN)
Deformation (mm)
Strain (%)
Stress (kPa)
Calculated stress (σE) (kPa)
Error (%)
Calculated stress (σE) (kPa)
Error (%)
4.97
2.0
5.7
0.50
0.50
0.46
0.50
0.00
8.21
3.5
10
0.86
0.86
0.51
0.86
0.07
14.88
7.0
20
1.70
1.70
0.26
1.70
0.03
20.52
10.5
30
2.54
2.54
0.06
2.55
0.26
25.35
14.0
40
3.38
3.37
0.24
3.39
0.29
29.47
17.5
50
4.21
4.21
0.13
4.23
0.44
33.06
21.0
60
5.04
5.04
0.09
5.07
0.49
36.22
24.5
70
5.86
5.86
0.05
5.90
0.68
38.94
28.0
80
6.68
6.69
0.07
6.73
0.8
41.48
31.5
90
7.50
7.50
0.01
7.57
0.87
43.58
35.0
100
8.30
8.30
0.05
8.39
1.10
45.50
38.5
110
9.10
9.10
0.01
9.21
1.21
47.13
42.0
120
9.88
9.88
0.02
10.02
1.41
48.67
45.5
130
10.65
10.65
0.01
10.82
1.57
49.98
49.0
140
11.41
11.41
0.04
11.60
1.67
51.03
52.5
150
12.15
12.14
0.06
12.37
1.80
51.95
56.0
160
12.86
12.86
0.01
13.12
2.03
52.75
59.5
170
13.56
13.56
0.01
13.86
2.19
53.40
63.0
180
14.24
14.24
0.0
14.58
2.35
53.90
66.5
190
14.89
14.90
0.05
15.27
2.58
54.32
70.0
200
15.52
15.53
0.06
15.95
2.79
54.65
73.5
210
16.12
16.14
0.11
16.61
3.05
54.78
77.0
220
16.70
16.72
0.10
17.25
3.29
54.86
80.5
230
17.25
17.27
0.11
17.86
3.55
54.94
84.0
240
17.78
17.79
0.07
18.45
3.78
54.84
87.5
250
18.28
18.28
0.01
19.02
4.02
54.72
91.0
260
18.74
18.74
0.00
19.55
4.32
54.47
94.5
270
19.18
19.16
0.08
20.05
4.56
54.01
98.0
280
19.57
19.55
0.09
20.53
4.89
53.61
101.5
290
19.93
19.90
0.14
20.97
5.19
53.23
105.0
300
20.24
20.21
0.13
21.37
5.57
52.48
108.5
310
20.50
20.49
0.07
21.73
6.00
51.83
112.0
320
20.73
20.71
0.08
22.05
6.38
51.02
115.5
330
20.91
20.90
0.05
22.34
6.81
50.29
119.0
340
21.04
21.04
0.01
22.57
7.29
49.21
122.5
350
21.12
21.13
0.06
22.77
7.80
48.25
126.0
360
21.16
21.18
0.08
22.92
8.31
47.17
129.5
370
21.15
21.17
0.10
23.03
8.87We next tested whether the above-derived engineering stress expression for the ILG is accurate and whether it can be used for the design of an ILG actuator or sensor. To verify the accuracy and practicability of the induced stress expression of the ILG, the deformation rate at each point in the uniaxial tensile test was substituted into the engineering stress expression to calculate the engineering stress. A comparison of the calculated uniaxial tensile stress with the experimental data is given in Table1.For the five-parameter Mooney-Rivlin model, the parametersc10, c01, c20, c11, and c02 are as follows.
(66)c10=2.3582,c01=−0.87482,c20=−0.066187,c11=0.19663,c02=0.21939.For the nine-parameter Mooney-Rivlin model, the parametersc10, c01, c20, c11, c02, c30, c21, c12, and c03 are as follows.
(67)c10=0.55361,c01=1.0009,c20=130.08,c11=−270.06,c02=142.91,c30=0.014438,c21=−0.17711,c12=−31.71,c03=17.028.As seen from Table1 in the comparisons of the calculated tensile stress with the experimental data, the relative error of the five-parameter Mooney-Rivlin model is less than 0.51% and the relative error of the nine-parameter Mooney-Rivlin model is no more than 8.87%.Figure8 shows that the stress-strain curve calculated by the five-parameter model almost coincides with the experimental curve. The first half of the stress-strain curve calculated by the nine-parameter model almost coincides with the experimental curve, while in the second half, the error between the stress-strain curve calculated by the nine-parameter model and the experimental curve becomes increasingly larger, but the error remains small.Figure 8
Comparison of the calculated tensile stress with the experimental data from the uniaxial tensile test.The above analyses indicate that the simulated values are consistent with the experimental values, making our derivation a feasible and credible stress formula for the calculation of the ILG properties.Due to our current laboratory conditions, only the uniaxial tensile test was performed. The next step is to improve the laboratory conditions in order to carry out the plane uniform tensile test and the Maxwell stress experiments. Further studies of the five-parameter and nine-parameter Mooney-Rivlin models will be carried out.


## 5.1. Production of Tensile Test Sample

The prepared solution was poured into a mold and polymerized into a gel under the irradiation of a UV lamp. The gel is white in color and is filled in the transparent glass mold as shown in Figure4. Then, the sample was cut to the size shown in Figure 5, with a thickness of 3.4 mm. As shown in Figure 6, the experimental instrument was a UTM2502 electronic universal testing machine. The mechanical sensor on the testing machine can achieve a precision of 0.1 mN, the displacement sensor has a precision of 0.001 mm, and the stretching rate is 500 mm/min.Figure 4
The gel after polymerization under a UV lamp.Figure 5
Test sample size (in mm).Figure 6
Uniaxial tensile test.


## 5.2. Experimental Results

As seen from Figure6, the ILG becomes longer and thinner as the load increases, which is consistent with the assumption that the material is incompressible.The tensile stress-strain curve of the ILG is shown in Figure7, and the average tensile strength (Young’s modulus) of the material obtained from the tensile stress-strain curve is 7.6 kPa.Figure 7
Stress-strain curve of the ionic liquid gel polymer.In the ionogel, BMIMBF4 exhibits a high level of hyperelastic toughness when the tensile deformation reaches 360%. The tensile test showed that the tensile properties of the gel increased with an increase in ZrO2 content. The increase in ZrO2 content should generate more cross-linking sites and higher conversion rates (shown in Figure 2), which contributes to the overall mechanical properties. As the content of ZrO2 increases, the tensile strength of the ILG increases, while the elongation rate decreases. Considering the above data, the optimum amount of ZrO2 to provide a large tensile strain and tensile strength is 3 wt%.


## 5.3. Stress Calculation

The relationship between the real stressσi and engineering stress σE is
(65)σi=σEλ,where σi is the true stress and σE is the engineering stress. The calculated stress given in Table 1 is the engineering stress.Table 1
Comparison of the calculated tensile stress with the experimental data from the uniaxial tensile test.
Measurement data
Experimental results
Five-parameter Mooney-Rivlin model
Nine-parameter Mooney-Rivlin model
Force(mN)
Deformation (mm)
Strain (%)
Stress (kPa)
Calculated stress (σE) (kPa)
Error (%)
Calculated stress (σE) (kPa)
Error (%)
4.97
2.0
5.7
0.50
0.50
0.46
0.50
0.00
8.21
3.5
10
0.86
0.86
0.51
0.86
0.07
14.88
7.0
20
1.70
1.70
0.26
1.70
0.03
20.52
10.5
30
2.54
2.54
0.06
2.55
0.26
25.35
14.0
40
3.38
3.37
0.24
3.39
0.29
29.47
17.5
50
4.21
4.21
0.13
4.23
0.44
33.06
21.0
60
5.04
5.04
0.09
5.07
0.49
36.22
24.5
70
5.86
5.86
0.05
5.90
0.68
38.94
28.0
80
6.68
6.69
0.07
6.73
0.8
41.48
31.5
90
7.50
7.50
0.01
7.57
0.87
43.58
35.0
100
8.30
8.30
0.05
8.39
1.10
45.50
38.5
110
9.10
9.10
0.01
9.21
1.21
47.13
42.0
120
9.88
9.88
0.02
10.02
1.41
48.67
45.5
130
10.65
10.65
0.01
10.82
1.57
49.98
49.0
140
11.41
11.41
0.04
11.60
1.67
51.03
52.5
150
12.15
12.14
0.06
12.37
1.80
51.95
56.0
160
12.86
12.86
0.01
13.12
2.03
52.75
59.5
170
13.56
13.56
0.01
13.86
2.19
53.40
63.0
180
14.24
14.24
0.0
14.58
2.35
53.90
66.5
190
14.89
14.90
0.05
15.27
2.58
54.32
70.0
200
15.52
15.53
0.06
15.95
2.79
54.65
73.5
210
16.12
16.14
0.11
16.61
3.05
54.78
77.0
220
16.70
16.72
0.10
17.25
3.29
54.86
80.5
230
17.25
17.27
0.11
17.86
3.55
54.94
84.0
240
17.78
17.79
0.07
18.45
3.78
54.84
87.5
250
18.28
18.28
0.01
19.02
4.02
54.72
91.0
260
18.74
18.74
0.00
19.55
4.32
54.47
94.5
270
19.18
19.16
0.08
20.05
4.56
54.01
98.0
280
19.57
19.55
0.09
20.53
4.89
53.61
101.5
290
19.93
19.90
0.14
20.97
5.19
53.23
105.0
300
20.24
20.21
0.13
21.37
5.57
52.48
108.5
310
20.50
20.49
0.07
21.73
6.00
51.83
112.0
320
20.73
20.71
0.08
22.05
6.38
51.02
115.5
330
20.91
20.90
0.05
22.34
6.81
50.29
119.0
340
21.04
21.04
0.01
22.57
7.29
49.21
122.5
350
21.12
21.13
0.06
22.77
7.80
48.25
126.0
360
21.16
21.18
0.08
22.92
8.31
47.17
129.5
370
21.15
21.17
0.10
23.03
8.87We next tested whether the above-derived engineering stress expression for the ILG is accurate and whether it can be used for the design of an ILG actuator or sensor. To verify the accuracy and practicability of the induced stress expression of the ILG, the deformation rate at each point in the uniaxial tensile test was substituted into the engineering stress expression to calculate the engineering stress. A comparison of the calculated uniaxial tensile stress with the experimental data is given in Table1.For the five-parameter Mooney-Rivlin model, the parametersc10, c01, c20, c11, and c02 are as follows.
(66)c10=2.3582,c01=−0.87482,c20=−0.066187,c11=0.19663,c02=0.21939.For the nine-parameter Mooney-Rivlin model, the parametersc10, c01, c20, c11, c02, c30, c21, c12, and c03 are as follows.
(67)c10=0.55361,c01=1.0009,c20=130.08,c11=−270.06,c02=142.91,c30=0.014438,c21=−0.17711,c12=−31.71,c03=17.028.As seen from Table1 in the comparisons of the calculated tensile stress with the experimental data, the relative error of the five-parameter Mooney-Rivlin model is less than 0.51% and the relative error of the nine-parameter Mooney-Rivlin model is no more than 8.87%.Figure8 shows that the stress-strain curve calculated by the five-parameter model almost coincides with the experimental curve. The first half of the stress-strain curve calculated by the nine-parameter model almost coincides with the experimental curve, while in the second half, the error between the stress-strain curve calculated by the nine-parameter model and the experimental curve becomes increasingly larger, but the error remains small.Figure 8
Comparison of the calculated tensile stress with the experimental data from the uniaxial tensile test.The above analyses indicate that the simulated values are consistent with the experimental values, making our derivation a feasible and credible stress formula for the calculation of the ILG properties.Due to our current laboratory conditions, only the uniaxial tensile test was performed. The next step is to improve the laboratory conditions in order to carry out the plane uniform tensile test and the Maxwell stress experiments. Further studies of the five-parameter and nine-parameter Mooney-Rivlin models will be carried out.


## 6. Conclusions

In this paper, an ILG is modeled by the hyperelastic nonlinear finite element model. The simulation results show that the Mooney-Rivlin model can well adapt to the constitutive relation of the material [30, 31]. The main advantage of the ILG is that the stress-strain curve can be obtained by the performance parameters of the material in a relatively short time, which provides a theoretical basis for the optimal design of a soft robot.The simulation results show that the average error between the calculated data and the experimental data is small, and the model has a good correlation with the experimental data. The model requires the input of the ILG material parameters. A standard uniaxial stretching method is used to obtain the desired ILG material parameters.In the future, we will seek a generalized algorithm for identifying the ILG mechanical properties. Notably, all the results of this study show that there is a good correlation between the 3D theoretical assumptions and the experimental conditions, which proves that our method can be used to optimize the design of a soft robot. This work opens a new path to study the performance of ILG soft actuators, which will be the direction of future work.



---

*Source: 8327867-2018-05-02.xml*