1. Introduction
In the field of soft actuators, a fiber-reinforced elastomeric enclosure (FREE) developed by McKibben et al. [
1,
2,
3,
4,
5,
6] has attracted attention because of its simple structure, high force-to-weight ratio, and high pressure. These unique capabilities enable a variety of potential medical device applications [
7,
8]. In addition, soft actuators conform to the surroundings and move with dexterity [
9,
10,
11]. Many soft actuators based on FREEs have been developed. They have a wide range of motions via combinations of programmed actions, including extension, bending, and twisting. Thus, it is essential to understand how to design a soft actuator to achieve a particular spatial motion. To make the design simpler and more systematic, there is a need to develop soft actuators that are easy to manufacture and program with respect to the complex environment.
As actuators based on FREEs can better match the behavior of biological muscles than others, research groups have begun to analyze the characteristics of soft actuators [
12,
13,
14]. Finite-element analysis [
15] is commonly used to optimize and guide the parameters for a soft actuator, and can also suggest designs. However, it can only analyze the movement of the foundation, and the calculations can be very large computationally.
Krishnan et al. [
16,
17] analyzed the mobility of generalized FREEs based on simple geometric relationships that result from the inextensibility of fibers and fluidic actuation. They examined the static behavior, such as extension, contraction, and twisting, by the principle of virtual work. However, this predicted only one-dimensional shape changes of soft actuators, without exploring actual application scenarios. To make it more flexible, some researchers arranged different fiber angles on the soft actuator to enable the single-drive actuator to perform a variety of actions [
18,
19,
20]. These reports were based on geometric modeling, rather than a force analysis of the soft actuator.
Uppalapati et al. [
21,
22] combined multiple modules in parallel to maximize the workspace. They introduced an inverse formulation using Cosserat rods [
23,
24,
25]. However, they did not establish a direct numerical relationship between pressure and deformation of a soft actuator. The actuator length was assumed to be constant, which is effective for bending or winding modules, but not for extension or twisting. Hence, there is a need to design a soft actuator that performs a specific motion, especially in complex environments.
To improve predictions, we used a geometric analysis model based on the principle of virtual work, combined with a mechanical analysis model based on Cosserat rod theory. Together, they predicted deformations of a soft actuator, especially gravity-induced shape deformation.
The advantages of using a combined model are:
(1) The geometric analysis estimates the extension, contraction, and torsion angle at any time. The Cosserat rod theory cannot perform this requirement.
(2) The Cosserat rod theory analyzes soft actuator deformation via mechanical characteristics. External forces can be integrated into the model. The geometric analysis model cannot account for the effects of external forces.
Hence, the combination of these analytical models provides more accurate deformation features. Extension, bending, and twisting modules are fabricated, and deformations are verified. Experimental data agreed well with the proposed model. Finally, we randomly combined modules and performed a series of target capture experiments to demonstrate this approach.
2. Analytical Modeling of Actuator Segments
In this section, we describe in detail the Cosserat rod theory and the geometric modeling of soft actuators.
2.1. The Cosserat Rod Theory
The centerline curve of a rod with
and material orientation
, as functions of a reference parameter
, are the main characteristics, as shown in
Figure 1.
The
and
are kinematic variables that represent the linear and angular rates of change of
and
in a local frame, respectively. Thus, they have the following relationship:
where
represents the skew symmetric matrix of the vector
. This is defined:
From the force and moment balance analysis of the rod in
Figure 1, the following can be obtained:
where
indicates the internal force and
indicates moment vectors.
indicates the tangential direction of a point on the soft actuator.
is the distributed force acting on the actuator and
is the applied moment distribution.
We use linear constitutive laws to map these variables to the internal forces and moments, as shown in Equations (5) and (6):
where
is the stiffness matrix for shear and extension, and is the stiffness matrix for bending and torsion. is the shear modulus, and , , and are the second moments of area for the tube cross sections about the principal axes.
It is necessary to make a brief explanation here. For the bending and twisting modules and the axial force . In other cases, .
The explicit model equations can be obtained by combining Equations (1)–(6) and the following:
where
, and because we assume no body moments,
.
Then, the boundary conditions are determined. At , , and . Matlab ode45 was used with the shooting method to solve the boundary-value problem.
2.2. Geometric Model Analysis
The relationship between soft actuator deformation and the fiber angle is given by:
where
, and
,
indicate axial elongation and radial elongation, respectively, and
,
indicate the helical angles of the fiber, as shown in
Figure 2.
is the net angular deformation of the actuator’s end.
and
are the initial actuator lengths.
,
are the deformed elemental axial length and radius, respectively.
The fibers were inextensible, so the relationship between input pressure and axial elongation can be obtained by the virtual work principle, as follows:
where
is the pressure and
is the volume of the actuator cavity.
is the volume of rubber, and
is Young’s modulus of the material. This was numerically inverted to solve for
, given
from Equations (9)–(13).
Modeling of Extension, Bending and Twisting
To model extension and bending, the helical angles of the fibers were set to
. If the angles are symmetrical, the soft actuator becomes longer under continuous air pressure and will not twist. Therefore, Equations (9)–(10) can be simplified as follows:
The only difference between the bending and extension modules is that the bending module must add another fiber as the binding layer. The added fiber allows the soft actuator to bend towards the fixed side. The process of change with the bending module is shown in
Figure 3.
The twisting module is a combination of the bending and rotating modules, as shown in
Figure 4. The helical angle of the fiber had an asymmetrical arrangement, with a motion similar to the predator action of an octopus arm. In this case, Equation (14) is not applicable and
can only be solved with Equations (9)–(13). Finally, the relationship between pressure and the deformation of the soft actuator can be obtained through the principle of virtual work under ideal conditions.
All of the above cases were obtained for an ideal state, without considering the effects of external forces, such as gravity. In
Figure 5, the gray cylinder on the left depicts the length of an actuator without the effect of gravity. The dark gray cylinder on the right depicts the change in length when considering the effect of gravity. Hence, there is a large difference when gravity is considered for a soft actuator.
To describe the above actions more accurately, the Cosserat rod theory was used to improve the soft actuator model under external loads. Initially, the equivalent axial force generated by internal pressurization of the tube was determined with Equation (16). This equation makes a preliminary hypothesis for the model derivation.
where
is the cross-sectional area of the soft actuator and is equal to
.
3. Manufacturing Soft Actuators
The actuators needed five materials: a hollow latex tube, fibers, rubber cement, a latex tube surface-treatment agent, and a metal rod, as shown in
Figure 6a,b. The latex tube had 15-mm inner and 20-mm outer diameters.
The soft actuators were fabricated by mounting the latex tube on the metal rod, which had a diameter equal to the inner diameter of the latex tube. To paste the fiber and latex tube together, a layer of treatment agent (J-740, DOLDSUN TIME) was applied to the surface of the latex tube. This agent enhanced the fiber restraint effect.
Then, the metal rod was fixed for machine control of the fiber-winding angle on the latex tube. The molding effects are shown in
Figure 7. Finally, we used rubber cement to adhere the fiber to the latex tube.
4. Experimental Validation
To fabricate extending, bending, and twisting actuators, we used a hollow latex tube with
and
. The fiber reinforcement was Kevlar. It was assumed that the fiber was inextensible. For the extension and bending modules, the fiber angles were set to
, and
for the twisting module. Detailed parameters are listed in
Table 1.
4.1. Module Testing
The modules were verified separately with experimental data. The model data of the combined virtual work principle and the Cosserat rod theory were compared and analyzed with experimental data.
4.1.1. Extension Module Test Experiment
The length of the extension module was set to 0.478 m. In
Figure 8, the red line is the experimental data for the actuator length, and the blue line is the length of the extension module predicted by the model. The experimental data were always greater than the data predicted by the model, and the maximum error was 0.0056 m. The reason for the large error may be because the extension module was not an ideal straight line during the fabrication process. There was a certain degree of bending. Therefore, with continuous pressure, the soft actuator was slightly bent along the elongation, which led to the actual length of the soft actuator always being less than the length predicted by the model. In general, the data predicted by the model matched the experimental results well, and were effective enough to predict the shape change of the model.
4.1.2. Bending Module Test Experiment
The length of the bending module was 0.378 m. Then, with input pressures
from 40 kPa to 110 kPa, we used the analytical model to predict the shape of the actuator, as shown in
Figure 9. The direction of gravity was positive along the z-axis. The abscissa represents the deviation of the soft actuator from the y-axis direction. Experimental data under corresponding conditions were also acquired. The position of the soft actuator must be determined via image processing.
In
Figure 9, there is good agreement between the analytical model and the experimental data. To better evaluate the accuracy of the analytical model, we used end-position errors to measure the model. Thus, the prediction model was evaluated by determining the position error of the end point between the model and experiment. The error was calculated by:
where
are the experimental position coordinates of the end of the
z-axis, and
are the position coordinates obtained by analytical modeling.
and
are determined similarly. Using this evaluation method, we obtained position errors of 0.0139 m, 0.012 m, 0.0403 m, and 0.0519 m under various input pressures (40–110 kPa).
The possible reason for the error is that the curvature of the soft actuator was not a constant curvature model with increasing input air pressure. Thus, there may have been deviations in the estimated position of the actuator. Another reason is that the cover at the end of the actuator had weight, which could make the predicted shape deviate from the actual shape.
4.1.3. Twisting Module Test Experiment
The length of the twisting module was 0.7 m. As its motion was three-dimensional, a helical radius was used to measure the model, extract the experimental data, and evaluate the accuracy of the model. The radius could be obtained by the virtual work principle and the Cosserat rod theory. Initially, the
of the actuator under ideal conditions was solved by using the virtual work principle. Then, the Cosserat rod theory was used to obtain the actuator attitude information. The curvature of the actuator centerline could be obtained by this method. The relationship between curvature and the helical radius could be obtained from previous reports [
22]:
where
is the twist per unit length.
represents the angle between the tangent direction at the end of the actuator and the z-axis direction of the global coordinate system.
The change in radius of the twisting module could be obtained by the above method, as shown in
Figure 10. The experimental data were always larger than the predicted data from the model. The maximum error was 0.01927 m, which may be because the cover at the end of the torsion module was made of metal and had a certain weight that affected the helical radius. However, the agreement between the spatial twisting model and experimental data was quite high, indicating the reliability of the proposed model.
4.2. Model Series Motion Experiment
The function of a single module is limited. Therefore, to highlight the operation of the soft actuator, we fabricated the extension, bending, and twisting modules at different scales. The twisting module was used as the end actuator, and it captured the target, as shown in
Figure 11a. Then, we combined the torsion module with the bending and elongation modules to grasp the target. In
Figure 11b,c, the twisting, bending, or extension modules could grasp the target.
5. Conclusions
Although there have been many reports on model analyses for soft actuators, especially for FREEs, few systematically analyzed separate modules, especially considering the effects of external forces. Here, we improved the description of spatial state changes in soft actuators. Firstly, a geometrically exact steady-state model (Cosserat rod theory) and the virtual work principle were used to predict deformation characteristics. The most important contribution is that the combination of these two models increased the design speed of the soft arm and improved the accuracy of shape prediction. Secondly, to verify the effectiveness of the model, we used the model to predict extension, bending, and twisting deformation actuators. Then, a comparison of the model analysis and the experimental data for the three modules verified that the result errors were 0.0056 m, 0.0519 m, and 0.01927 m. Finally, to expand the function of the soft actuator, the extension, bending, and twisting modules were randomly combined to grasp objects in space by controlling the pressure of the three modules.
A soft actuator constructed using FREEs can easily change its shape by means of external forces. Hence, in future work, we plan to change this structure into a module with variable stiffness, then we will optimize the structure of the soft actuator to minimize the volume. Finally, we plan to examine dynamic modeling analysis of the soft actuator because most modeling of the flexible actuator is still based on quasi-static motion.
Author Contributions
X.H. was the supervisor providing funding and administrating the project, and he reviewed and edited the manuscript. Z.C. carried out the investigation, methodology, analysis, and validation, and wrote the original manuscript. A.Z. and T.L. have made great contribution to the application of MATLAB in simulation. Z.Q. built a platform for the verification model and helped to make different soft actuator modules. S.L. helped to test the parameters of the soft actuator, and also helped to perform some verification experiments. All authors have read and agreed to the published version of the manuscript.
Funding
This research was funded by Guangxi Natural Science Foundation Program 2018JJA160218, the National Natural Science Foundation of China under Grant 51965014, 51765014 and 51705094.
Acknowledgments
The authors would like to thank the Peng Cheng Laboratory for providing support of material parameters.
Conflicts of Interest
The authors declare no conflict of interest.
References
- Chou, C.-P.; Hannaford, B. Measurement and modeling of McKibben pneumatic artificial muscles. IEEE Trans. Robot. Autom. 1996, 12, 90–102. [Google Scholar] [CrossRef] [Green Version]
- Carlo Ferraresi, W.F.; Walter Franco, W.; Bertetto, A. Flexible pneumatic actuators: A comparison between the McKibben and the straight fibres muscles. J. Robot. Mechatron. 2001, 13, 56–63. [Google Scholar] [CrossRef]
- Kothera, C.S.; Jangid, M.; Sirohi, J.; Wereley, N.M. Experimental characterization and static modeling of McKibben actuators. J. Mech. Des. 2009, 131, 091010. [Google Scholar] [CrossRef]
- Tondu, B. Robust and accurate closed-loop control of McKibben artificial muscle contraction with a linear single integral action. Actuators 2014, 3, 142–161. [Google Scholar] [CrossRef] [Green Version]
- Walker, J.; Zidek, T.; Harbel, C.; Yoon, S.; Strickland, F.S.; Kumar, S.; Shin, M. Soft robotics: A review of recent developments of pneumatic soft actuators. Actuators 2020, 9, 3. [Google Scholar] [CrossRef] [Green Version]
- Furukawa, S.; Wakimoto, S.; Kanda, T.; Hagihara, H. A Soft Master-Slave Robot Mimicking Octopus Arm Structure Using Thin Artificial Muscles and Wire Encoders. Actuators 2019, 8, 40. [Google Scholar] [CrossRef] [Green Version]
- Roche, E.T.; Wohlfarth, R.; Overvelde, J.T.; Vasilyev, N.V.; Pigula, F.A.; Mooney, D.J.; Bertoldi, K.; Walsh, C.J. A bioinspired soft actuated material. Adv. Mater. 2014, 26, 1200–1206. [Google Scholar] [CrossRef] [PubMed]
- Polygerinos, P.; Wang, Z.; Galloway, K.C.; Wood, R.J.; Walsh, C.J. Soft robotic glove for combined assistance and at-home rehabilitation. Robot. Auton. Syst. 2015, 73, 135–143. [Google Scholar] [CrossRef] [Green Version]
- Zhang, Z.; Philen, M. Modeling, analysis and experiments of inter yarn compaction effects in braided composite actuators. J. Compos. Mater. 2013, 47, 3211–3225. [Google Scholar] [CrossRef]
- Mortier, K. Braided Pneumatic Muscles for Rehabilitation Apparatus; Politecnico di Torino: Turin, Italy, 2014. [Google Scholar]
- Kim, J.Y.; Mazzoleni, N.; Bryant, M. Modeling of Resistive Forces and Buckling Behavior in Variable Recruitment Fluidic Artificial Muscle Bundles. Actuators 2021, 10, 42. [Google Scholar] [CrossRef]
- Marchese, A.D.; Tedrake, R.; Rus, D. Dynamics and trajectory optimization for a soft spatial fluidic elastomer manipulator. Int. J. Robot. Res. 2016, 35, 1000–1019. [Google Scholar] [CrossRef]
- Polygerinos, P.; Wang, Z.; Overvelde, J.T.; Galloway, K.C.; Wood, R.J.; Bertoldi, K.; Walsh, C.J. Modeling of soft fiber-reinforced bending actuators. IEEE Trans. Robot. 2015, 31, 778–789. [Google Scholar] [CrossRef] [Green Version]
- Singh, G.; Krishnan, G. An isoperimetric formulation to predict deformation behavior of pneumatic fiber reinforced elastomeric actuators. In Proceedings of the 2015 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS), Hamburg, Germany, 28 September–2 October 2015; pp. 1738–1743. [Google Scholar]
- Moseley, P.; Florez, J.M.; Sonar, H.A.; Agarwal, G.; Curtin, W.; Paik, J. Modeling, design, and development of soft pneumatic actuators with finite element method. Adv. Eng. Mater. 2016, 18, 978–988. [Google Scholar] [CrossRef]
- Krishnan, G.; Bishop-Moser, J.; Kim, C.; Kota, S. Kinematics of a generalized class of pneumatic artificial muscles. J. Mech. Robot. 2015, 7, 041014. [Google Scholar] [CrossRef]
- Singh, G.; Krishnan, G. A constrained maximization formulation to analyze deformation of fiber reinforced elastomeric actuators. Smart Mater. Struct. 2017, 26. [Google Scholar] [CrossRef]
- Connolly, F.; Polygerinos, P.; Walsh, C.J.; Bertoldi, K.J.S.R. Mechanical programming of soft actuators by varying fiber angle. Soft Robot. 2015, 2, 26–32. [Google Scholar] [CrossRef] [Green Version]
- Connolly, F.; Walsh, C.J.; Bertoldi, K. Automatic design of fiber-reinforced soft actuators for trajectory matching. Proc. Natl. Acad. Sci. USA 2017, 114, 51–56. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Singh, G.; Krishnan, G. Designing Fiber-Reinforced Soft Actuators for Planar Curvilinear Shape Matching. Soft Robot. 2020, 7, 109–121. [Google Scholar] [CrossRef] [PubMed]
- Uppalapati, N.K.; Krishnan, G. Design and Modeling of Soft Continuum Manipulators Using Parallel Asymmetric Combination of Fiber-Reinforced Elastomers. J. Mech. Robot. 2021, 13, 011010. [Google Scholar] [CrossRef]
- Uppalapati, N.K.; Krishnan, G. Towards pneumatic spiral grippers: Modeling and design considerations. Soft Robot. 2018, 5, 695–709. [Google Scholar] [CrossRef]
- Rucker, D.C.; Webster III, R.J. Statics and dynamics of continuum robots with general tendon routing and external loading. IEEE Trans. Robot. 2011, 27, 1033–1044. [Google Scholar] [CrossRef]
- Gazzola, M.; Dudte, L.; McCormick, A.; Mahadevan, L. Forward and inverse problems in the mechanics of soft filaments. R. Soc. Open Sci. 2018, 5, 171628. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Till, J.; Aloi, V.; Rucker, C. Real-time dynamics of soft and continuum robots based on Cosserat rod models. Int. J. Robot. Res. 2019, 38, 723–746. [Google Scholar] [CrossRef]
| Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
© 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).