- Control Systems and Robotics Laboratory (CSRL), Department of Mechanical Engineering, School of Engineering, Hellenic Mediterranean University, Heraklion, Greece
Introduction: Soft pneumatic actuators (SPAs) play a pivotal role in soft robotics due to their unique characteristics of compliance, flexibility, and adaptability. There are plenty of approaches that examine the modeling parameters of SPAs, aiming to optimize their design and, thus, achieve the most advantageous responses. Current optimization methods applied to SPAs are usually performed individually for each design parameter without considering the simultaneous effect all parameters can have on the output performance. This modeling shortcoming is essential to be addressed since customized SPAs are used in a variety of applications, each with different output requirements.
Methods: This study provides a generalized design optimization framework for modeling the SPA performance for their motion profiles, the produced strain energy while being deformed, and their stiffness characteristics. Utilizing experimentally validated finite element methods, all geometrical and material parameters of the models are investigated in response surface methodology optimization using the central composite design approach.
Results: The results showcase the entire design space of omnidirectional SPAs, along with their output performance, providing guidelines to the end user for design optimization.
Discussion: The offering of this modeling process for SPAs can be adapted to the demands of any potential application and ensure the best performance with respect to the required output responses.
1 Introduction
For the past couple of decades, soft robotic technology has been oriented toward setting up robotic mechanisms made of highly deformable materials. These mechanisms are characterized by high dexterity and safety, exhibit a variety of motions that offer applicability to a wide range of tasks including dynamic tasks, and are sensitive to the physical interaction, allowing usage in difficult-to-access environments. The effectiveness of these soft robotic mechanisms hinges on the performance of their actuators. To date, soft pneumatic actuators (SPAs) are the most common and widely used soft actuators. SPAs have internal cavities, i.e., chambers, for pneumatic actuation, which can be deformed using negative or positive pressures. The pneumatic actuation approach is particularly favored due to its numerous advantages, including lightweight construction, rapid response times, and straightforward and cost-effective implementation, as well as robust pressure resistance (Xavier et al., 2022). SPAs perform motions such as bending, elongation, contraction, twisting, and combinations of the aforementioned. They can be manufactured by silicone rubber casting approaches, directly by 3D-printing methods using flexible filaments or elastomeric resins, or by thermal bonding of thermoplastics (Niiyama et al., 2015). SPAs find utilization in a range of applications, including minimally invasive surgery (Shi et al., 2022; Decroly et al., 2021; Kwok et al., 2022), rehabilitation (Polygerinos et al., 2013; Cianchetti et al., 2013; Fidinillah et al., 2023), safe human–robot interaction scenarios (Noritsugu, 2005; Ku et al., 2024), biomedical devices (Greef and Delchambre, 2009), and handling of delicate objects (Shintake et al., 2018; Deimel and Brock, 2013; Zhou et al., 2021).
The most popular categories of SPAs are the reinforced and pneumatic network (PneuNet) actuators. A well-known subcategory of SPAs comprises omnidirectional actuators, which were proposed by Suzumori et al. (1991a) and Suzumori et al. (1991b) and continue to attract considerable scholarly attention (Xavier et al., 2021b; Lee et al., 2024; Luo et al., 2023). In the existing literature, the most common omnidirectional actuators are monolithic structures with multiple internal pneumatic/hydraulic chambers (Abidi et al., 2018; Decroly et al., 2021). However, they can also be formed through assemblies of multiple single SPAs connected in series or parallel (O’Brien and Lane, 2001; Kalisky et al., 2017) or in the form of polyhedral networks (Liu et al., 2024).
The design and the modeling of SPAs used in various applications and how it is feasible to optimize their design to achieve certain responses have been extensively discussed by a great number of researchers. In particular, numerous studies have investigated the effect of design parameters, either geometrical or structural, on responses such as the bending performance and the force output of SPAs (Elsayed et al., 2014b; Polygerinos et al., 2015; Shi et al., 2022). In these studies, the most frequent module geometry is chosen as a cylinder with internal chambers (Elsayed et al., 2014a) acting as the pressurization vessels. Aside from this type of SPAs, PneuNets are also examined extensively (Hu et al., 2018; Zhu et al., 2022; Jing et al., 2022; Fidinillah et al., 2023). In the design optimization stages, this is implemented in most cases by using finite element analysis (Yan et al., 2016; Hu et al., 2018; Fidinillah et al., 2023). In there, each geometrical parameter is examined sequentially with respect to the desired output response by testing a range of values while preserving fixed values for the rest. However, this approach may lead to different output results depending on the sequence of investigation being followed. Therefore, this approach fails to consider all the factors that affect the SPA behavior simultaneously. Additionally, every application carries its unique requirements for SPAs, which need to be satisfied under optimum conditions. In this context, there is a lack of a generalized framework regarding the design methodology that could be easy to customize to meet application-specific requirements and achieve the desired responses as well. In this study, we aim to address this need and take the initial stride toward an overall fundamental design methodology of SPAs for optimizing responses such as bending ability, strain energy consumption, and stiffness, which constitute major key elements for evaluating and characterizing the SPA performance and profiles (Fidinillah et al., 2023; Zhu et al., 2022; Jing et al., 2022), based on each application demand, so far lacking in the scientific literature.
In this work, we contribute to a comprehensive design framework for omnidirectional SPAs, and the paper is structured as follows: Section 2.1 focuses on exploring a wide parameter space associated with the modeling of such actuators and providing broader design, geometrical, and structural recommendations for SPAs that the end user could take under consideration; in Section 2.2, guided by the central composite design approach, we simulate omnidirectional actuators with different ranges of values for their constituent geometrical parameters, as well as shore hardness material properties, by using finite element method (FEM) analysis under fully parameterized models. These responses are evaluated and optimized using a response surface methodology (RSM), after experimentally validating the FEM simulation results in Section 2.3; and in Section 3, we use analysis of variance (ANOVA) to examine their responses regarding bending angle performance, actuating pressure-produced strain energy, and the stiffness variation that can be delivered with respect to the parameters under investigation. In this way, we provide a complete guide to the end user for optimizing the behavior of their soft omnidirectional actuators.
2 Materials and methods
2.1 SPA concept design–parameter space
As a starting point, the geometrical parameter space is set up for the desired configuration of the actuator. In this framework study, an exemplary SPA design is adopted, but their parameters can be modified according to specific needs and applications. Here, it is considered as a monolithic module with internal chambers and two distal caps attached to both sides. The module has the following geometrical parameters to consider, as shown in Figure 1: 1) the module length (
2.1.1 Fixed parameters
The module length
Additionally, in this study, several
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1418484/frobt-11-1418484-HTML/image_m/frobt-11-1418484-g002.jpg)
Figure 2. The SPA module diameter-to-module length ratio is tested using FEM simulations for five ratios (0.1, 0.3, 0.5, 0.7, and 1) to investigate (A) their bending ability at constant pressure (displacement values closer to 1 indicate a SPA that bends easier), (B) their energy consumption profile (strain energy values closer to 1 indicate a SPA that generates more internal strains to bend), and (C) their stiffness response when subjected to external pushing force (stiffness values closer to 1 indicate a SPA that deflects less). Ratios of 0.7 and 1 cannot be simulated reaching 90° and, therefore, are omitted in (B,C). The gray and red zones correspond to ranges where there is difficulty for SPAs due to the small dimensions to be fabricated via conventional fabrication methods or bend due to similar dimensions in the two orthogonal axes, respectively. Indicatively, FEM of SPAs with 0.3 and 0.5 ratios is visualized below the graphs for each case.
2.1.2 Parameters under investigation
Utilizing FEM analysis, various module cross sections are explored to identify the best module cross section to utilize. All experimented actuator shapes have a single internal concentric chamber, which is identical in shape to their outer shape, and a cross-sectional area equal between all of them. These single-chamber actuators are simulated, assuming hyperelastic properties for their material, and pressure is applied to their cavities. This results in the transformation of all their unpressurized cross sections to a circular cross section when pressurized, as shown in Figure 3A. From an energy perspective, the achieved circular shape geometrically has no edges or corners, so the distribution of stresses and strains is more uniform and does not concentrate in these regions. With only this in mind, any initial cross section can be utilized if one aims to create a single-chamber SPA that, when pressurized, elongates along its main axis. In this instantiation of the SPA, the main criterion for initial geometrical shape selection should therefore be the ease of fabrication each shape could offer, as well as any application-specific limitations. To achieve more than one direction of motion for the SPA, one is required to increase the number of internally pressurizable chambers. As discussed, a single chamber within the actuated SPA makes it elongate along the main axis as the strains are distributed uniformly around the periphery. In the case of two chambers, the SPA can either elongate along its main axis, if both chambers are equally pressurized, or can bend, along a single plane, toward the opposite direction of the pressurized chamber with the highest pressure. The degrees of freedom (DOFs) therefore include elongation and roll. Introducing three or more chambers, the SPA, apart from elongation and roll, gains pitch as a DOF. It is noted that a SPA becomes steerable in the entire three-dimensional space at a minimum of three chambers, but one can opt for additional chambers to maximize the control resolution of the SPA. The correlation between the number of chambers and the DOFs is shown in Figure 3B. Nonetheless, it is significant that aside from these actuation chambers, there is usually the need for an additional chamber that serves the purpose of enabling tools and other components to pass through the entire length of the SPA, such as cameras and grippers in surgical laparoscopic tools (Yan et al., 2016; Cianchetti et al., 2013; Kwok et al., 2022).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1418484/frobt-11-1418484-HTML/image_m/frobt-11-1418484-g003.jpg)
Figure 3. (A) FEM analysis of various SPAs with different cross-sectional shapes (circular, semicircular, square, and triangular). When the internal chamber is pressurized, a circular cross-sectional shape is reached in all cases; (B) pressurized and unpressurized motion states of SPAs with a different number of chambers. To achieve an omnidirectional SPA design, a minimum of three internal pressurizable chambers are required so that it is steerable in the three-dimensional space offering elongation, roll, and pitch capabilities; (C) FEM investigation on various cross sections for omnidirectional SPAs with three chambers (hexagon, oval, and annulus sector cross sections are presented). The cross section of a chamber with an annulus sector design presents the maximum displacement when pressurized. This displacement is associated with the maximum bending ability of the omnidirectional SPA when a single chamber is pressurized; (D) FEM investigation as in (C) for the amount of strain energy produced when pressurizing a single chamber of the omnidirectional SPA to reach a 90° bending angle. The annulus sector design is shown to generate the least strain in the material that comprises that SPA.
For the cross-sectional shape of the chamber(s), the selection can be made based on the following requirements: i) the wall thickness to be of equal size across its span to allow uniform distribution of the stresses toward the outer surface of the module while being pressurized and ii) less material to be used so that the volume of the pressurizable chambers is increased, and thus, the pressure area required to perform bending, or any other combined motion pattern, can be decreased, consequently limiting the strain exerted on the deformed areas of the SPA. Taking these requirements into consideration, for the case of the one-chamber SPA, the optimal chamber cross-section shape derived is circular. In the case of two symmetrical chambers, the best resulting chamber shape is semicircular. Lastly, in the case of three chambers, these are best to be at a 120° symmetrically spaced arrangement with an annulus sector shape. This geometrical distribution of area in the chambers is the logical derivation to occur when one splits the SPA in continuously more chambers while preserving a uniform outer wall thickness. However, in this study, the hypothesis of the annulus sector being the optimal chamber shape for omnidirectional SPAs with three or more chambers is tested against a number of different chamber cross sections of the same area using FEM analysis. Each cross section is tested in two scenarios that, in subsequent sections, are analyzed in detail. The first scenario aims to investigate the maximum allowable bending angle the SPA can reach while only one chamber is pressurized at constant pressure. The second is to examine the energy cost, i.e., the produced strain energy for bending to a 90° angle, when actuating a single chamber. The results of these tests are shown normalized in Figure 3C and Figure 3D, respectively, and highlight the optimal cross section chamber shape within the context of the maximum bending angle and the least produced strain energy, revealing the superiority of the annular sector over all other major shapes. It is also observed that the annulus sector preserves its superiority as a chamber shape for omnidirectional SPAs with a minimum of three pneumatic chambers irrespective of the scale/size of the chamber area. This validation is performed through FEM analysis of the same SPAs by doubling the chamber area, and the results are given in Supplementary Material.
Design of the selected configuration: The SPA performs best when in a cylindrical form, with at least three pneumatic chambers of an annulus sector cross-sectional shape that are radially arranged around the central axis. Furthermore, the SPA carries one main circular central channel so that various tools can pass through the module and reach the distal end with minimum influence in bending performance, as shown in Figure 4.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1418484/frobt-11-1418484-HTML/image_m/frobt-11-1418484-g004.jpg)
Figure 4. Omnidirectional SPA configuration along with the cross section of the module. Design variables A, B, C, and D represent the wall thickness, inner thickness, central channel radius, and arc angle, respectively.
2.1.3 Cross-section shape-driving parameters
Adding the additional cavity in the center of the SPA generates the selected cross-section shape with a section of annulus. This leads us to further explore the effect of the driving parameters that define the particular cross section. The outer radius of the annulus sector is associated with the thickness between the outer surface of the module and the outer surface of the chamber of the SPA, defined henceforth as wall thickness. The inner radius of the annulus sector is linked to the thickness between the interior surface of the chamber and the surface of the central channel, defined henceforth as inner thickness. The radius of the central channel is also examined. The final variable to consider is the angle of the arc of the annulus sector. Therefore, a total of four variables are selected to investigate their effect on SPA responses, which are depicted as
2.2 Response surface methodology for design optimization
The RSM comprises a set of mathematical and statistical tools suitable for modeling and analyzing complex relationships between multiple independent input variables and the responses they produce on the dependent variable
where
On the other hand, simple linear and interaction models are sometimes adequate to be considered the approximation relationship between factors (input variables) and responses (output variables) and, therefore, efficiently predict the optimum values. In this study, the process of optimization includes the estimation of coefficients, prediction of responses, and evaluation of the developed model and is performed using the statistical analysis suite Design-Expert 13 (Stat-Ease Inc., Minneapolis, MN, United States).
2.2.1 Central composite design
CCD is used as the experimental design methodology as it facilitates the efficient exploration of the factor space, maintaining a balance between precision and resource requirements. It is a fractional factorial design with a center point attached to a group of stars, or axial points, which allows the estimation of the curvature (Njoku and Otisi, 2023; Bhattacharya, 2021). In particular, four independent/input variables (i.e., factors) are considered, namely, wall thickness (A), inner thickness (B), radius of the central channel (C), and arc angle (D), each with three levels, as shown in Figure 4. The examined Y responses of the designated experiment are i) the bending angle at a fixed pressure value; ii) the produced strain energy for the SPA to bend at a 90° bending angle; and iii) the stiffness of the SPA. For the latter response, the aim is to examine the deflection of the SPA, when it has already bent 90°, and then, a force is exerted on its distal top surface. Measuring the displacement of the SPA and given the force applied, one can calculate the stiffness that is equal to the ratio of the force to the provoked deformation. The number of experimental runs (i.e., the series of tests an experiment encompasses) required for CCD is given by
2.2.2 Finite element method analysis
The parametric modeling of the SPA using Python scripting and the conduction of simulations to investigate its behavior in the three aforementioned scenarios are performed using SIMULIA Abaqus FEA. The first case inspects the produced strain energy while the SPA deforms to a bending angle of 90°. In the second case, a force of one unit is applied at the top end-cap surface while the SPA is pressurized at a 90° angle. Thus, it is feasible to measure the deflection that is created and, as a result, extend our study in estimating the stiffness characteristics of the SPA. Lastly, the third case aims to study, given a fixed pressure value, the maximum bending angle the SPA reaches. Consequently, 25 models, which correspond to the number of runs that the CCD generated, are modeled, and their resulting data from the FEM simulations are extracted to complete the corresponding responses in the CCD. Each SPA model is meshed using solid quadratic tetrahedral elements (Abaqus element type C3D10H) with a seed size of 1.5. As for the material, the SPA is a monolithic structure assigned with silicone rubber, which is defined by mass density, hyperelasticity, isotropy, and incompressibility, regarding the requisite material properties. In the examples used to construct the design framework of this study, the Yeoh model is used as the hyperelastic model to capture the highly extensible and elastic behavior of silicone rubber due to its ability to approach large strain problems (Xavier et al., 2021a). The strain energy function of the Yeoh model
where
2.3 Experimental validation of FEM simulations
Physical experiments with the SPAs are performed to validate the accuracy and reliability of the computational models, and their experimental output responses are compared with the corresponding outputs from the FEM simulations. An additional quasistatic investigation is implemented where the values for the A, B, C, and D variables (Figure 4) are preserved to a constant value while testing two different shore hardness hyperelastic materials, namely, Dragon Skin 20A (Smooth-On, Inc.) and Dragon Skin 30A (Smooth-On, Inc.) with elongation at breaks 620% and 364%, respectively. Both have a mass density of
2.3.1 SPA fabrication process
Two hyperelastic materials, material A (Dragon Skin 30A) and material B (Dragon Skin 20A), are used to fabricate the SPA configuration. Several different SPA fabrication processes are reported in the literature, and for this study, the general approach suggested by Xie et al. (2018); Decroly et al. (2021); Cianchetti et al. (2013); and Yan et al. (2016) is adopted with a few custom modifications. Both silicone elastomer materials require a catalyst to obtain their final solid state. To achieve this, the two parts are mixed together at a 1:1 ratio. The mixture is subjected to vacuum to degas and is then casted in the appropriate additive manufactured mold. Once the silicone is demolded, the SPA module is constructed. To be complete, the module is sealed on both ends with distal end caps. The module is integrated with the distal end cap by placing the module vertically inside a small cylindrical container and filled with silicone at a depth of the desired cap length, as shown in Figure 5. The curing of the silicone is expedited in a convection oven set at 45°C.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1418484/frobt-11-1418484-HTML/image_m/frobt-11-1418484-g005.jpg)
Figure 5. Omnidirectional SPA fabrication process that utilizes multi-part molds for silicone casting.
2.3.2 FEM simulations vs experiments
For the bending-angle performance evaluation, a motion capture analysis system (Impulse X2E, PhaseSpace, CA, United States) is used to track the SPA motion in the physical world. Three active motion markers are placed at the base end cap where the SPA is held upright so as to generate a local reference coordinate system. Another single motion marker is placed at the surface of one side of the distal cap of the SPA. The experiment involves a single chamber of the SPA being pressurized until it reaches a 90° bending angle. A comparison between the FEM model and the experimental data is shown in Figure 6A, where the top end-cap trajectories during the particular motions are recorded. The results indicate the good agreement of the FEM model with the physical prototypes for both materials A and B, indicatively, 14% and 5%, respectively. In the SPA case prepared with material A, the deviation concerning the requisite actuation pressure between theoretical and experimental behaviors observed is 9.6%, while in material B, a 11% deviation from the experimental value is spotted in FEM simulation.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1418484/frobt-11-1418484-HTML/image_m/frobt-11-1418484-g006.jpg)
Figure 6. Physical testing of two SPAs with similar geometry but of different shore hardness materials bending at 90°. (A) SPA bending trajectories are experimentally obtained using active LED motion markers placed at the distal end and base of the actuators for the ground plane reference. The trajectories are compared with those extracted from FEM simulations presenting a small deviation in both cases. (B) Experimental stiffness setup that utilizes a uniaxial force sensor and a linear translation stage. Stiffness values are based on the measurement of force and displacement by the force sensor. (C,D) Comparison between experimental and computational results for materials A and B, which illustrate the SPA stiffness and displacement profiles. The values are normalized across the two materials for ease of comparison.
For the stiffness investigation, a contact interaction setup, including a universal testing machine, is used to measure the force exerted on the top surface of the SPA and vice versa, as shown in Figure 6B. The apparatus contains a digitally fabricated fixture where the SPA rear end cap is secured, while a force sensor can come in contact with the distal end cap. In this way, the force and displacement measurements are recorded by the sensor by gradually compressing the actuated SPA. As shown in Figure 6C, for material A, the root mean square error (RMSE) between the experimental measurements and the FEM simulation results for the stiffness is 13% and 4.4% for the displacement due to force. Similarly, as shown in Figure 6D, for material B, the RMSE is 4% and 3.5% for both stiffness and displacement occurring due to force application, respectively. Comparatively, the displacement of material B, due to its lower shore hardness value, as expected, is greater than the displacement of material A. The stiffness responses present a magnitude-increasing pattern up to a certain force threshold and then starts decreasing. This behavior is observed in both materials A and B, in both FEM simulations and experiments, and can be attributed to the high compliance of the system occurring at the early stages of the contacted force that finds equilibrium as the force increases. Simultaneously, the resistivity at specific local areas of the SPA is also increased due to the compression of its silicone body. As a result, it is feasible to inspect how material properties affect the SPA deflection in the case of contact interactions. However, it needs to be clarified that in the proposed design framework, the optimization process attempts to minimize the displacement due to contact by adjusting all of the investigated parameters while also incorporating material selection as a design parameter according to the desired stiffness profile demands.
3 Results
3.1 ANOVA and model fit statistics
The significance of the variables and their interactions is determined using ANOVA with a significance threshold of 95% and a p-value of 0.050. The mathematical models are derived using the ANOVA table and used for optimization purposes. In particular, for the response of displacement after the application of force, a quadratic model should be considered initially for model fitting. However, there is the need for some statistical model reduction, indicated by the difference between the predicted R-squared and the adjusted R-squared value being more than 0.2. In the statistical analysis, if the particular difference is less than 0.2, then it is stated that the model can fit the data and can be assuredly used to interpolate. Therefore, the final CCD obtained for this response is a reduced quadratic. For the responses of strain energy and bending angle, the suggested models are linear. The F-value for every model, which is the ratio of the mean square for the individual term to the mean square for the residual, indicates that all the obtained results are statistically significant with also p-value variables of less than 0.05 (Saha et al., 2010; Sadhukhan et al., 2016). Moreover, adequate precision is associated with the signal-to-noise ratio (SNR). If the SNR is greater than 4, the model has a strong enough signal to be used for optimization. All the statistical measures used for the examination of the performance of this design framework, i.e., predicted and adjusted R-squared values, model F-values, p-values, and adequate precision, are given in a detailed table in Supplementary Material. Apart from that, the model selection along with the statistical significant variables for each response is summarized in Table 1 as they are required to provide a better understanding and justify Figure 7, which is explained in Section 3.2.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1418484/frobt-11-1418484-HTML/image_m/frobt-11-1418484-t001.jpg)
Table 1. Selection of models for each response and their statistical significant terms based on their p-value
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1418484/frobt-11-1418484-HTML/image_m/frobt-11-1418484-g007.jpg)
Figure 7. Response surface methodology (RSM) results using a central composite design (CCD) as an optimization approach. (A) Correlation of design of experiment (DoE) curves: design variables are presented with respect to each response, i.e., displacement after force applied, strain energy, and bending angle. Axes are normalized between 0 and 1, so 0 indicates the minimum value of the desired range based on the demands of the application an SPA will be used for, and 1 indicates the maximum value correspondingly. (B) Interactions of factors (design variables) and their impact on responses. Surface colors represent the gradient range from the lowest, corresponding to blue shades, to the greatest response value, corresponding to red shades. Red dots correspond to response values above the surface, and the pink ones correspond to those below the surface. Red-outlined areas highlight the surface location where the pair of the corresponding variables leads to the optimized, that is, minimization of the displacement of force, minimization of the produced strain energy, and maximization of the bending angle.
3.2 Correlation curves
Figure 7 encapsulates the effect of the investigated independent variables on the aforementioned responses individually and interactively. Analyzing the results shown in Figure 7A, starting from the response of displacement after force, variations in the central channel radius exhibit negligible influence. However, when the wall thickness or inner thickness increases, the displacement decreases. Reversely, when the arc angle increases, the displacement increases as well. As for the response of strain energy, it is observed that as the wall thickness and arc angle increase, the strain energy also increases. The central channel radius and inner thickness do not influence this response significantly. Lastly, concerning the response of bending, when increasing the wall thickness and central channel radius, the bending angle is decreased, while the arc angle, when increasing, leads to its augmentation. The inner thickness, on the other side, does not affect the response. In Figure 7B, two main points are highlighted: 1) the correlations between the independent variables refer only to the statistical significant model terms of the present ANOVA results, depicted in Table 1. The statistical insignificant terms do not influence the responses, so the interaction with them is omitted. This implies that the latter terms may have an impact on SPA behavior; however, with the present DoE, these terms are considered insignificant; 2) for the responses of strain energy and bending angle, where the suggested models are linear, an adjustment from a linear to a two-factor interaction (2FI) model is necessary, including only the products of statistically significant model terms. This adjustment, notably incorporating
4 Discussion
In this work, a design methodology and an optimization framework were presented that provide guidelines and insights to the end user working toward the modeling aspects of a SPA. The methodology starts agnostic of cross-sectional shapes and the number of chambers required and provided the geometrical parameter space one must consider according to a specific application in mind for their SPAs. The framework showcased, among many other characteristics, the small importance in the performance of the external shape during pressurization, while the choice of number of pressurizable chambers was found highly correlated with the DOFs one desires to achieve. In the case of those omnidirectional SPAs, which are highly reported in the literature, the minimum number of chambers was determined along with the optimal chamber shape from the perspective of energy and deformation during pressurization. The best performing chamber cross section was found to be the annulus sector. Examining the cross section of an omnidirectional SPA, comprising three pneumatic chambers with an annulus sector shape and a central circular channel that allows application-specific passage through the SPA body, a number of variables were further investigated for their effect on output responses. The variables and their responses focused on the bending performance with a fixed pressure value and the produced strain energy while bending at a fixed 90° angle. Furthermore, the stiffness of the SPA in combination with its deformed shape after pressurization was a topic of interest for many researchers (Shi et al., 2024; Gao et al., 2023; Lamping and de Payrebrune, 2023; Luo et al., 2023), and in this study, the FEM analysis was used to illustrate the stiffness variation across the design spectrum. The investigation used data obtained from the FEM simulations and was implemented using a central composite design, which appertains to the response surface methodology. These design methods enabled us, for the first time, to consider all the variables of the SPA at once, which was underexplored in the existing scientific literature. Yet, in order to consider the FEM approach trustworthy, a series of experimental studies were conducted with physical prototypes, and the performance results were compared to the simulation results, thus validating that the FEM models adequately capture the real-word conditions.
In the experiments, different shore-hardness SPAs were examined to investigate the effect of material properties in these simulation cases. From the existing literature, an important observation was that material model coefficients often vary and cover a wide range (Marechal et al., 2020; Xavier et al., 2021a). Thus, we note that the characterization of hyperelastic materials depends on several factors, such as the fabrication quality of the test specimens, i.e., dog bones, and the conditions of the experiment. In particular, weather conditions, unequal or insufficient mixing of the silicone elastomer with the catalyst, air trapped in the silicone mixture, overcuring, and examination of different ranges of strain and fitted return different material coefficients for the same hyperelastic material models and can lead to these notable variations. In this context, the coefficients in hyperelastic models can occur in two ways. One way is the conventional one, where the coefficients are determined using uniaxial, biaxial, and shear test data. The test data originate from test specimens being deformed (in most cases elongated) by a tensile tester. Therefore, considering all the before-mentioned factors that can affect the coefficient values, the results will satisfy only the specific conditions of the end user’s specific study and application. An alternative option, in case there is no available material testing equipment, one can opt for fine-tuning the existing literature model coefficients, visualizing the SPA behavior via FEM models while obtaining experimental feedback concerning the approximation of the appropriate values.
The ability to optimize the responses according to the requirements of the application should be emphasized while analyzing the impact of these variables on the responses using CCD. Optimization was defined by either maximizing, minimizing, keeping in range, or targeting a value of a response. Regarding the response of the produced strain energy, it is meaningful to minimize first since less strain energy makes the SPA more energy efficient, durable as an overall system, and facilitates more accurate control. One of the advantages of the SPAs is their high compliance. However, the same compliance becomes an unwanted feature once one requires the SPA to interact with the environment. We demonstrated in this study how to passively alter the design parameter/variable to achieve the minimum displacement response when forces act at a SPA. The model provided the desired output parameters for the SPA to preserve high resistance when accepting an external load, something that is desired by many applications at certain stages of SPA utilization. To achieve greater stiffness, the wall thickness and inner thickness of the SPA should be kept close to the upper bound, whereas the arc angle should be kept close to the lower bound. Finally, with respect to the bending motion in free space of the SPA, it is found useful in most cases to maximize the bending angle under a fixed pressure value. Hence, for maximum bending, the wall thickness and central channel radius should be kept low and the arc angle high. Altogether, this optimization process provided an essential component for modern soft robotics, where SPAs and more specifically omnidirectional SPAs find utilization in more and more domains, and thus, the needs and requirements spread across a wide parameter range.
The present study considered exclusively the design, implying that active stiffness control, neither internal nor external, e.g., through reinforcements/constraints, has been applied. Although reinforcements could also be incorporated in this design framework guide, we believe that as a more generalized methodology for SPAs is being developed, it is crucial to deplete all the possible adjustments that a plain hyperelastic material actuator and its geometry can provide in the first place and then proceed to any additional modifications. Future work could seek to include such passive and active reinforcements to further improve the SPA responses and also examine them in dynamic scenarios.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.
Author contributions
MM: writing–original draft, conceptualization, formal analysis, investigation, methodology, software, validation, visualization, and writing–review and editing. PP: writing–review and editing, conceptualization, funding acquisition, methodology, project administration, resources, supervision, and writing–original draft.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. The work was supported by the European Union through the PALPABLE Project under Grant No. 101092518.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/frobt.2024.1418484/full#supplementary-material
References
Abidi, H., Gerboni, G., Brancadoro, M., Fras, J., Diodato, A., Cianchetti, M., et al. (2018). Highly dexterous 2-module soft robot for intra-organ navigation in minimally invasive surgery. Int. J. Med. Robotics Comput. Assisted Surg. 14, e1875. doi:10.1002/rcs.1875
Arezzo, A., Mintz, Y., Allaix, M. E., Arolfo, S., Bonino, M. A., Gerboni, G., et al. (2016). Total mesorectal excision using a soft and flexible robotic arm: a feasibility study in cadaver models. Surg. Endosc. 31, 264–273. doi:10.1007/s00464-016-4967-x
Bhattacharya, S. (2021). “Central composite design for response surface methodology and its application in pharmacy,” in Response surface methodology in engineering science. Editor P. Kayaroganam (Rijeka: IntechOpen).
Cianchetti, M., Ranzani, T., Gerboni, G., De Falco, I., Laschi, C., and Menciassi, A. (2013). “Stiff-flop surgical manipulator: mechanical design and experimental characterization of the single module,” in 2013 IEEE/RSJ international conference on intelligent robots and systems, Tokyo, Japan, November 03–07, 2013.
Decroly, G., Lambert, P., and Delchambre, A. (2021). A soft pneumatic two-degree-of-freedom actuator for endoscopy. Front. Robotics AI 8, 768236. doi:10.3389/frobt.2021.768236
Deimel, R., and Brock, O. (2013). A compliant hand based on a novel pneumatic actuator. IEEE Int. Conf. Robotics Automation, 2047–2053. doi:10.1109/ICRA.2013.6630851
Diodato, A., Brancadoro, M., Rossi, G. D., Abidi, H., Dall’Alba, D., Muradore, R., et al. (2018). Soft robotic manipulator for improving dexterity in minimally invasive surgery. Surg. Innov. 25, 69–76. doi:10.1177/1553350617745953
Elsayed, Y., Lekakou, C., Geng, T., and Saaj, C. (2014a). Design optimisation of soft silicone pneumatic actuators using finite element analysis. , 44–49. doi:10.1109/AIM.2014.6878044
Elsayed, Y., Vincensi, A., Lekakou, C., Geng, T., Saaj, C. M., Ranzani, T., et al. (2014b). Finite element analysis and design optimization of a pneumatically actuating silicone module for robotic surgery applications. Soft Robot. 1, 255–262. doi:10.1089/soro.2014.0016
Fidinillah, T., Risangtuni, A., Putra, N., Virdyawan, V., Suprijanto, d., and Putri, S. (2023). Design of soft pneumatic actuator in hand rehabilitation robot with fem-based modeling. 109–113. doi:10.1109/ICA58538.2023.10273084
Gao, W., Kang, J., Wang, G., Ma, H., Chen, X., Kadic, M., et al. (2023). Programmable and variable-stiffness robotic skins for pneumatic actuation. Adv. Intell. Syst. 5, 2300285. doi:10.1002/aisy.202300285
Gerboni, G., Ranzani, T., Diodato, A., Ciuti, G., Cianchetti, M., and Menciassi, A. (2015). Modular soft mechatronic manipulator for minimally invasive surgery (mis): overall architecture and development of a fully integrated soft module. Meccanica 50, 2865–2878. doi:10.1007/s11012-015-0267-0
Greef, A., and Delchambre, A. (2009). Towards flexible medical instruments: review of flexible fluidic actuators. Precis. Eng. 33, 311–321. doi:10.1016/j.precisioneng.2008.10.004
Hu, W., Mutlu, R., Li, W., and Alici, G. (2018). A structural optimisation method for a soft pneumatic actuator. Robotics 7, 24. doi:10.3390/robotics7020024
Jing, X., Chen, S., Zhang, C., and Xie, F. (2022). Increasing bending performance of soft actuator by silicon rubbers of multiple hardness. Machines 10, 272. doi:10.3390/machines10040272
Kalisky, T., Wang, Y., Shih, B., Drotman, D., Jadhav, S., Aronoff-Spencer, E., et al. (2017). Differential pressure control of 3d printed soft fluidic actuators. , 6207–6213. doi:10.1109/IROS.2017.8206523
Ku, S., Song, B.-H., Park, T., Lee, Y., and Park, Y.-L. (2024). Soft modularized robotic arm for safe human–robot interaction based on visual and proprioceptive feedback. Int. J. Robotics Res. 0, 1128–1150. doi:10.1177/02783649241227249
Kwok, K.-W., Wurdemann, H., Arezzo, A., Menciassi, A., and Althoefer, K. (2022). Soft robot-assisted minimally invasive surgery and interventions: advances and outlook. Proc. IEEE 110, 871–892. doi:10.1109/JPROC.2022.3167931
Lamidi, S., Olaleye, N., Bankole, Y., Obalola, A., Aribike, E., and Adigun, I. (2022). “Applications of response surface methodology (rsm) in product design, development, and process optimization,” in Response surface methodology (Rijeka: IntechOpen).
Lamping, F., and de Payrebrune, K. M. (2023). A novel and practicable approach for determining the beam parameters of soft pneumatic multi-chamber bending actuators. Appl. Sci. 13, 2822. doi:10.3390/app13052822
Lee, K., Bayarsaikhan, K., Aguilar, G., Realmuto, J., and Sheng, J. (2024). Design and characterization of soft fabric omnidirectional bending actuators. Actuators 13, 112. doi:10.3390/act13030112
Liu, X., Han, X., Hong, W., Wan, F., and Song, C. (2024). Proprioceptive learning with soft polyhedral networks. Int. J. Robotics Res. 0, 02783649241238765. doi:10.1177/02783649241238765
Luo, X., Song, D., Zhang, Z., Wang, S., and Shi, C. (2023). A novel distal hybrid pneumatic/cable-driven continuum joint with variable stiffness capacity for flexible gastrointestinal endoscopy. Adv. Intell. Syst. 5, 2200403. doi:10.1002/aisy.202200403
Marechal, L., Balland, P., Lindenroth, L., Petrou, F., Kontovounisios, C., and Bello, F. (2020). Toward a common framework and database of materials for soft robotics. Soft Robot. 8, 284–297. doi:10.1089/soro.2019.0115
Niiyama, R., Sun, X., Sung, C., An, B., Rus, D., and Kim, S. (2015). Pouch motors: printable soft actuators integrated with computational design. Soft Robot. 2, 59–70. doi:10.1089/soro.2014.0023
Njoku, C. N., and Otisi, S. K. (2023). “Application of central composite design with design expert v13 in process optimization,” in Response surface methodology. Editor P. Kayarogannam (Rijeka: IntechOpen).
Noritsugu, T. (2005). “Pneumatic soft actuator for human assist technology,” in Symposium on fluid power.
O’Brien, D., and Lane, D. (2001). “3d force control system design for a hydraulic parallel bellows continuum actuator,” in Proceedings 2001 ICRA. IEEE International Conference on Robotics and Automation, Seoul, Korea, May 21–26, 2001.
Polygerinos, P., Lyne, S., Wang, Z., Nicolini, L., Mosadegh, B., Whitesides, G., et al. (2013). Towards a soft pneumatic glove for hand rehabilitation. 1512–1517. doi:10.1109/IROS.2013.6696549
Polygerinos, P., Wang, Z., Overvelde, J. T. B., Galloway, K. C., Wood, R. J., Bertoldi, K., et al. (2015). Modeling of soft fiber-reinforced bending actuators. IEEE Trans. Robotics 31, 778–789. doi:10.1109/TRO.2015.2428504
Rezaee, R., Maleki, A., Jafari, A., Mazloomi, S., Zandsalimi, Y., and Mahvi, A. (2014). Application of response surface methodology for optimization of natural organic matter degradation by uv/h2o2 advanced oxidation process. J. Environ. health Sci. Eng. 12, 67. doi:10.1186/2052-336X-12-67
Runciman, M., Darzi, A., and Mylonas, G. P. (2019). Soft robotics in minimally invasive surgery. Soft Robot. 6, 423–443. doi:10.1089/soro.2018.0136
Sadhukhan, B., Mondal, N. K., and Chattoraj, S. (2016). Optimisation using central composite design (ccd) and the desirability function for sorption of methylene blue from aqueous solution onto lemna major. Karbala Int. J. Mod. Sci. 2, 145–155. doi:10.1016/j.kijoms.2016.03.005
Saha, P., Chowdhury, S., Gupta, S., Kumar, I., and Kumar, R. (2010). Assessment on the removal of malachite green using tamarind fruit shell as biosorbent. Clean. - Soil Air Water 38, 437–445. doi:10.1002/clen.200900234
Shi, J., Gaozhang, W., and Wurdemann, H. (2022). Design and characterisation of cross-sectional geometries for soft robotic manipulators with fibre-reinforced chambers. doi:10.1109/RoboSoft54090.2022.9762157
Shi, J., Shariati, A., Abad, S.-A., Liu, Y., Dai, J. S., and Wurdemann, H. A. (2024). Stiffness modelling and analysis of soft fluidic-driven robots using lie theory. Int. J. Robotics Res. 43, 354–384. doi:10.1177/02783649231200595
Shintake, J., Cacucciolo, V., Floreano, D., and Shea, H. (2018). Soft robotic grippers. Adv. Mater. 30, 1707035. doi:10.1002/adma.201707035
Suzumori, K., Iikura, S., and Tanaka, H. (1991a). “Development of flexible microactuator and its applications to robotic mechanisms,” in Proceedings. 1991 IEEE International Conference on Robotics and Automation, Sacramento, CA, April 09–11, 1991.
Suzumori, K., Iikura, S., and Tanaka, H. (1991b). “Flexible microactuator for miniature robots,” in Proceedings. IEEE Micro Electro Mechanical Systems, Nara, Japan, January 30–02, 1991.
Xavier, M. S., Fleming, A. J., and Yong, Y. K. (2021a). Finite element modeling of soft fluidic actuators: overview and recent developments. Adv. Intell. Syst. 3, 2000187. doi:10.1002/aisy.202000187
Xavier, M. S., Tawk, C. D., Yong, Y. K., and Fleming, A. J. (2021b). 3d-printed omnidirectional soft pneumatic actuators: design, modeling and characterization. Sensors Actuators A Phys. 332, 113199. doi:10.1016/j.sna.2021.113199
Xavier, M. S., Tawk, C. D., Zolfagharian, A., Pinskier, J., Howard, D., Young, T., et al. (2022). Soft pneumatic actuators: a review of design, fabrication, modeling, sensing, control and applications. IEEE Access 10, 59442–59485. doi:10.1109/ACCESS.2022.3179589
Xie, R., Su, M., Zhang, Y., and Guan, Y. (2018). “3d-psa: a 3d pneumatic soft actuator with extending and omnidirectional bending motion,” in 2018 IEEE international conference on robotics and biomimetics (ROBIO), Kuala Lumpur, Malaysia, December 12–15, 2018.
Yan, J., Dong, H., Zhang, X., and Zhao, J. (2016). “A three-chambed soft actuator module with omnidirectional bending motion,” in 2016 IEEE international conference on real-time computing and robotics (RCAR), Angkor Wat, Cambodia, June 06–10, 2016.
Ye, G., Ma, L., Li, L., Liu, J., Yuan, S., and Huang, G. (2017). Application of Box–Behnken design and response surface methodology for modeling and optimization of batch flotation of coal. Int. J. Coal Prep. Util. 40, 131–145. doi:10.1080/19392699.2017.1350657
Zhou, L., Ren, L., Chen, Y., Niu, S., Han, Z., and Ren, L. (2021). Bio-inspired soft grippers based on impactive gripping. Adv. Sci. 8, 2002017. doi:10.1002/advs.202002017
Keywords: design framework, design optimization, response surface methodology, finite element analysis, central composite design, soft pneumatic actuators, omnidirectional soft actuators
Citation: Moutousi M and Polygerinos P (2024) Omnidirectional soft pneumatic actuators: a design and optimization framework. Front. Robot. AI 11:1418484. doi: 10.3389/frobt.2024.1418484
Received: 16 April 2024; Accepted: 22 August 2024;
Published: 10 September 2024.
Edited by:
Cecilia Laschi, National University of Singapore, SingaporeReviewed by:
Chaoyang Song, Southern University of Science and Technology, ChinaNicholas Kellaris, Artimus Robotics, United States
Copyright © 2024 Moutousi and Polygerinos. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Panagiotis Polygerinos, polygerinos@hmu.gr