Next Article in Journal
An Efficient Optimal Derivative-Free Fourth-Order Method and Its Memory Variant for Non-Linear Models and Their Dynamics
Next Article in Special Issue
Numerical Analysis of the Effect of the Vortex Finder on the Hydrocyclone’s Split Water Ratio and Separation Performance
Previous Article in Journal / Special Issue
An Analysis of Numerical Homogenisation Methods Applied on Corrugated Paperboard
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Hierarchical Design Framework for the Design of Soft Robots

by
Philip Frederik Ligthart
* and
Martin Philip Venter
Department Mechanical and Mechatronic Engineering, Stellenbosch University, Stellenbosch 7600, South Africa
*
Author to whom correspondence should be addressed.
Math. Comput. Appl. 2023, 28(2), 47; https://doi.org/10.3390/mca28020047
Submission received: 15 February 2023 / Revised: 8 March 2023 / Accepted: 16 March 2023 / Published: 21 March 2023
(This article belongs to the Special Issue Current Problems and Advances in Computational and Applied Mechanics)

Abstract

:
This paper demonstrates the effectiveness of a hierarchical design framework in developing environment-specific behaviour for fluid-actuated soft robots. Our proposed framework employs multi-step optimisation and reduced-order modelling to reduce the computational expense associated with simulating non-linear materials used in the design process. Specifically, our framework requires the designer to make high-level decisions to simplify the optimisations, targeting simple objectives in earlier steps and more complex objectives in later steps. We present a case study, where our proposed framework is compared to a conventional direct design approach for a simple 2D design. A soft pneumatic bending actuator was designed that is able to perform asymmetrical motion when actuated cyclically. Our results show that the hierarchical framework can find almost 2.5 times better solutions in less than 3% of the time when compared to a direct design approach.

1. Introduction

Soft robots have a design space that is complex and co-dependent. As a result of the highly compliant materials, successful designs require consideration of both shape and behaviour [1]. Soft robots are more biologically compatible than traditional robots, with potential applications in fields such as human-to-machine interfaces, search and rescue, and other delicate operations [2]. Despite the complex shape-dependant control, high degree of deformation, and non-linear material behaviour characteristic, it is still the convention to use a build-and-test development strategy [3]. However, the traditional build-and-test design paradigm is expensive and does not allow rapid exploration of the design space. Attempts have been made to use computational resources to aid in the design process, but these also have drawbacks, such as high computational expense. Researchers have made recent advancements in modelling [4,5] and optimising [6,7] techniques for soft robots that improve accuracy and computational efficiency [8]. Nonetheless, these methods often exist in isolation and provide limited benefits for designing a full robot [9]. For this reason, Pinskier and Howard [10] propose a hierarchical design framework to make the design process more efficient and allow for future modelling techniques’ integration. In this paper, we demonstrate a hierarchical design framework to enable the efficient design of soft robots. Using a hierarchical framework, we aim to improve the design process and make it easier to incorporate new modelling methods.

2. Background

Soft robotics is a promising field of robotics that offers unique benefits not achievable with traditional rigid robots. Soft robots possess both active and passive compliance due to their complex material properties, which are costly and difficult to achieve with rigid robots [11]. This compliance renders them safer for human–machine interactions and ideal for delicate handling operations such as fruit picking [2,12]. Additionally, their material properties allow them to deform and adapt to their environment, making them suitable for navigating unstructured terrains. Hence, they hold great potential for search-and-rescue missions [13]. The advent of 3D printing technology has enabled the complete manufacture of soft robots, including the incorporation of sensory equipment within their morphology, allowing for embodied intelligence and control [14,15]. Given these advantages, soft robotics is poised to have a significant impact on the future of robotics.
In the field of soft robots pneumatic actuation is the dominant actuation technology [16]. For this reason, we focus on pneumatic networks (PNs) in this paper. There are two main types of PNs, namely the slow pneumatic network (sPN) [17] and the fast pneumatic network (fPN) [18]. Many soft robots which have been developed are created by assembling multiple of these actuators together. Marchese et al. [19] created a soft robotic fish tail capable of escape manoeuvres using two sPNs, Schiller et al. [20] designed a gecko-inspired soft robot using six fPNs connected together, and Ilievski et al. [17] created a soft robotic gripper using six connected sPNs. In many of the designs found in the literature, a trial-and-error design approach is followed [3]. Both physical testing and FE simulations are used during the trial-and-error approach. This direct design process is feasible when the problem is simple and small; however, when a more elaborate problem is to be solved, it is difficult to manually iterate a design to achieve an acceptable solution, so to this end, authors have used computational optimisation techniques, such as gradient-based optimisation or genetic algorithms for the iteration of the designs [6,7,21].
Using computational resources to optimise the design process is advantageous for the designer, as it alleviates the need for manual iterations. Frameworks, such as those developed by Runge et al. [22], have been developed which incorporate both FE simulations and kinematic modelling techniques to allow for the automated design of soft robotics. However, such frameworks remain largely parametric and hence are limited in terms of design space exploration. Evolutionary design methodologies have also been explored with good success [10,21,23]. Many works make use of simulated environments to evolve the soft bodies; this is pragmatic but leads to a simulation-to-reality gap [10]. Other design frameworks have been developed that use topology optimisation and evolutionary methods to automatically design soft actuators and components [6,7]; however, when considering material nonlinearities, these methods are also plagued with a high computational burden. For an in-depth review of the autonomous design of soft robots, the reader is referred to Pinskier and Howard [10].

3. Materials and Methods

The hierarchical framework described here builds on the work of Pinskier and Howard [10]. The hierarchical framework presented here is aimed at condensing the large design problem into smaller sub-problems that are simpler or computationally more efficient to solve. The smaller sub-problems can then be assembled back into the large design problem. In this approach, the designer makes high-level decisions, and computational resources are utilised to optimise smaller sub-problems.
The hierarchical framework is comprised of multiple levels, shown in Figure 1. Level 1 constitutes the desired outcome of the entire design process. High-level decisions are made to simplify the design. The simplification of level 1 becomes level 2. The same holds when moving from level 2 to level 3. Simplifications can comprise splitting the design into sub-components, making use of symmetry, and repetition, among others. Each simplification is represented as a hierarchical step (HS) within the level. There can be multiple HSs in a single level, as in the case where a desired robot is split into multiple sub-components. For example, a rehabilitation glove can potentially be split into five sub-components, each representing a different finger. Once the design has been sufficiently reduced or simplified in the forward pass of the framework, each HS is then solved, starting from the last level. Each next HS that is solved is assembled back into the complete design in the first level. This is the backward pass through the framework.
Within the framework, each HS is treated as a separate design problem. A single HS is shown in Figure 2. It must have a clearly defined target, properties, and constraints that are enforced by the designer based on application specifications and availability. These constraints can also help to reduce the design space. The design problem can often be framed as an optimisation problem, but it is important to make sure the optimisation is well defined. Optimisations can lead to unexpected or undesirable outcomes if the target properties or constraints are inaccurate as noted by Lehman et al. [24]. Clearly defining properties at each level of the hierarchy helps keep the design pragmatic and avoids unnecessary optimisation. It is important to remember that each HS in the design framework has its own distinct target; this is opposed to a single optimisation step, where a single objective function is to be achieved. When reducing the large problem into smaller sub-problems, the designer must carefully consider the effects of the simplifications. Using repeating units or incorporating symmetry can have consequences that restrict the design space.
Within the design framework, FE simulations can be utilised. However, when the FE simulation requires a non-linear analysis, the computational demands can be substantial. To address this challenge, reduced-order models can be employed to enhance the efficiency of optimising sub-problems. Many such models can be found in the literature, for example, piecewise constant curvature [4], Cosserat rods [5], topology optimisation [7], or a bespoke reduced model specifically tailored for the application at hand can be developed.

4. Design of Asymmetrical Bending Motion Actuator

To demonstrate the above-described framework, a soft pneumatic actuator capable of performing asymmetrical bending motion is designed. Asymmetric motion is the objective for two reasons. Firstly, asymmetrical motion is essential for soft robots that locomote or swim. During forward motion, a walking or swimming robot requires its legs to follow an asymmetrical path. Furthermore, such an actuator presents a more complex design than a conventional soft pneumatic bending actuator (SPBA). Figure 3 provides an overview of the hierarchical framework implemented in the case study. The framework starts with the desired design in the first level, progressively breaking it down into simpler problems. The optimisation commences from the last level, working back to the first level. A summary of each level’s function is presented in Table 1 and is discussed below.

4.1. Forward Pass

When transitioning from the first level to the second level, a simplification is made by using repeating elements instead of a single continuous morphology. This simplification limits the possible deformations, but it is deemed acceptable for the desired outcome. A custom, simplified model was created to represent the actuator and to optimise the problem. The model is a 2 × 2 transformation of points, representing the movement of an element being actuated. Despite the non-linear path followed by a soft actuator, the start and end points of each actuation are only considered at this point, allowing for the use of a linear transformation. This simplification represents the model well enough to allow for rapid exploration with a low computation expense. While this approach represents a simplification of the actual model, it is permissible since later stages in the framework allow us to fully describe the model and perform more complete optimisations [25]. When the transformation is applied to a point (vector), it moves in space. Figure 4a shows an example of the two right-hand nodes of a unit being transformed. Figure 4b shows the result of when the two right-hand side nodes of each element in a cascade of connected elements are transformed. The figure shows that the model provides a good representation of an sPN during actuation. Mathematically, a point x 0 , y 0 is transformed into the point x t , y t using the equation
x t y t = x s c l x s h y s h y s c l Transformation x 0 y 0
where s c l represents scale and s h represents shear in the respective coordinate directions.
We used the reduced model to explore the design space. When examining the path traced by the actuator’s tip, we observed that if there is only a starting point and one actuation point, the resulting path would be a line connecting the start and actuation point, i.e., start–actuation point–start. To create a closed area, at least one additional point is required. Thus, it was found that asymmetrical motion can be achieved by applying two consecutive transformations followed by a return transformation to the original state. Specifically, the path would be, start–actuation point 1–actuation point 2–start. These findings show that asymmetrical motion is possible, and to achieve it, the actuator will require at least two discrete actuations. Using the knowledge gained from level 2, a single unit will be designed that has two discrete chambers within it. The two chambers will be actuated independently. The properties and constraints for the forward pass are shown in Table 2.

4.2. Backward Pass

The backward pass solves the simplified problems generated in the forward pass. This starts with the chamber layout optimisation in level 3. The chosen chamber layout is to have two independent chambers placed above and below one another within the repeating element. Two chambers were used to allow for two distinct actuations. This is graphically represented in Figure 5a. The two chambers are actuated according to Figure 5b. The x and y coordinates of each chamber corner node point are the design variables which are optimised. The design variables are represented as
x = x 1 , y 1 , x 2 , y 2 , x 3 , y 3 , x 4 , y 4 , x 5 , y 5 , x 6 , y 6 , x 7 , y 7 , x 8 , y 8
The constraints from Table 2 were added to the optimiser as well as additional bounds to each node of the chambers to enforce the minimum wall thickness constraint, as seen in Figure 5a. An objective function is constructed where asymmetrical bending motion is desired. The asymmetry is quantified as the area enclosed by the tip of the actuator during a complete actuation cycle. The coordinates of the actuator tip are extracted during actuation, and the surveyor’s area formula [28] is used to calculate the enclosed area. The formula makes use of the Cartesian coordinates of each vertex of a simple polygon. The formula is
A = 1 2 i = 1 n y i + y i + 1 x i x i + 1
As the optimisation is a minimisation problem, the objective function is the negative of the calculated area A.
An additional simplification of reducing the number of elements in the model was made in order to reduce the computational burden of simulating the actuator. A collection of three units is the smallest possible representation that describes the actuator. When the repeating unit is cascaded to form the full-length actuator, it is in effect only the second element that is repeated, as the first and last elements have unique boundary conditions [29]. This is more clearly represented in Figure 6.
A non-linear finite element solver was used to simulate the units when the pressures were applied to the internal voids according to Figure 5b. The optimiser was then used to minimise the objective function. An optimisation was initialised using random parameters, and multiple restarts were used. The properties of the non-linear solver and the optimiser can be found in Appendix A.
Figure 7 shows an optimised three-element actuator. The red line shows the path followed by the tip of the actuator during actuation. When the optimised element is cascaded to the full nine elements, the resulting actuation is shown in Figure 8. The figures show that the goal of achieving asymmetrical motion is achieved.

4.3. Results and Discussion

To evaluate the effectiveness of the hierarchical design framework, or hierarchical framework method, we also used a direct design approach as a comparison. The direct method was not designed using the hierarchical framework. Instead, it used all the nodes of all the chambers as design variables, and the entire actuator was simulated. This approach is referred to as the direct method. Figure 9 illustrates the graphical representation of the two different methods.
In Figure 10, we compare the run times and solutions of the two methods. As shown, the hierarchical framework method is capable of finding suitable solutions faster than the direct method. The lower solution time of the hierarchical framework method enables the efficient execution of multiple restarts, which is advantageous when dealing with numerous local optima scattered within the design space. Due to the enormous size of the direct method optimisation problem, solution times are far longer, making multiple restarts infeasible. While direct optimisation may be able to provide better solutions, it is unlikely that it will find such an optimal solution in a comparable time frame. Table 3 summarises the properties of the optimisation problems in the two scenarios.
From the framework above, the possibility arises of using a different HS at each level. By adding a new HS to one level, subsequent levels can leverage its results to uncover diverse design solutions. This occurs as the effects of the HS spread throughout the framework and enable the exploration of new areas of the design space. The capability to add new HSs highlights the scalability of the framework. For instance, the framework can incorporate changes such as recombining units in level 2, or altering the chamber layout in level 3, as shown within the hierarchical framework in Figure 11.

5. Conclusions

From the case study presented above, it was shown that using the proposed framework significantly reduced the development time when compared to a direct design approach. The reduced development time also leads to better overall performance, as the reduced computational time can allow for multiple restarts, which aid in more rapidly exploring the design space. The framework allows the designer to make high-level decisions that reduce the computational burden placed on the optimiser by introducing simplifications and creating reduced model representations. Because simplifications inherently reduce the design space, it aids in rapidly finding solutions but can also restrict the final design. The restrictions can be mitigated by using the final solutions of previous simpler designs as a starting point for more direct optimisations. This leverages the reduced design space to obtain a good starting point for a more direct optimisation. While the framework in this study did not incorporate force feedback, it is an essential element for the functionality of soft actuators that interact with the environment. These robots may need to propel themselves forward or manipulate objects, and force feedback is essential to achieving these goals. However, integrating force feedback into the design problem can add additional complexities. Therefore, the next stage in the development of the framework should involve the incorporation of force feedback along with the output displacement to achieve functional soft robots.

Author Contributions

Conceptualisation, P.F.L. and M.P.V.; methodology, P.F.L.; software, P.F.L.; formal analysis, P.F.L.; investigation, P.F.L.; resources, M.P.V.; data curation, P.F.L.; writing—original draft preparation, P.F.L.; writing—review and editing, M.P.V.; visualisation, P.F.L.; supervision, M.P.V.; project administration, M.P.V.; funding acquisition, M.P.V. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded in part by the National Research Foundation of South Africa (NRF) grant number 129381.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to this research forming part of a larger study.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
MDPIMultidisciplinary Digital Publishing Institute
FEFinite element
HSHierarchical step
SPBASoft pneumatic bending actuator
PNPneumatic network
sPNSlow pneumatic network
fPNFast pneumatic network
MMFDModified method of feasible directions
HFmHierarchical framework method
DmDirect method

Appendix A. Optimisation Loop Details

Appendix A.1. Non-Linear Solver

The non-linear solver used was Marc, and the pre- and post-processor used was Mentat [30]. A planar plane strain analysis was performed, using quad4 elements with a maximum size of 0.9 mm. The elements used were Marc element type 118, which are reduced integration and Herman formulation elements. The pressures were applied as edge loads, and self-contact was activated. Convergence properties used in the solver are shown in Table A1.
Table A1. Non-linear solver convergence properties.
Table A1. Non-linear solver convergence properties.
PropertyValue
Iterative procedureFull Newton–Raphson
Relative residual force tolerance0.1

Appendix A.2. Optimiser

The design optimiser tool (DOT) from Van der Plaats Research & Development was used as the optimiser [31]. Method 1 of the optimiser, which represents the modified method of feasible directions (MMFD), was used. The optimiser properties are shown in Table A2. The optimiser was not tuned to suit any of the different design methods; rather, the default optimiser parameters were used to eliminate any bias.
Table A2. DOT optimiser properties.
Table A2. DOT optimiser properties.
PropertyValue
Absolute Convergence Criteria m a x 0.0001 × A B S F 0 , 1.0 × 10 20
Relative convergence criteria0.001
Maximum optimisation iterations100
Gradient calculationForward difference
Relative finite difference step0.001

Appendix A.3. Function Calls and Evaluation Time

When running an optimisation, the efficiency of the optimiser can be determined by the number of function calls required to converge to an optimum. In this work, this was not the approach that was taken. The function call is a non-linear FE simulation of an sPN, which has changing geometry. Due to the changing geometry, the problem and, hence, the run time of the simulation change.In the case where the internal voids are small and there is little deformation, the simulation run time can be in the order of 30 s. The run time can increase by a factor 4 or more if the internal voids are larger and there is self contact which occurs. Hence, function calls are not equivalent in this optimisation pipeline, and the optimisation time was rather used as a comparison metric. To ensure that all tests run times could be compared, a single computer was used for all simulations.

Appendix A.4. Computation Hardware

The computer that was used for all simulations is detailed in Table A3 and Table A4.
Table A3. Hardware specifications.
Table A3. Hardware specifications.
ProcessorIntel(R) Core(TM) i5-6200U CPU @ 2.30 GHz 2.40 GHz
Installed RAM16.0 GB
System type64-bit operating system, ×64 based processor
Table A4. Operating system.
Table A4. Operating system.
Operating SystemWindows
Edition10 Pro
Version22H2
OS build19,045.2311

References

  1. Rus, D.; Tolley, M.T. Design, fabrication and control of soft robots. Nature 2015, 521, 467–475. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Tawk, C.; Alici, G. A Review of 3D-Printable Soft Pneumatic Actuators and Sensors: Research Challenges and Opportunities. Adv. Intell. Syst. 2021, 3, 2000223. [Google Scholar] [CrossRef]
  3. Ma, L.K.; Zhang, Y.; Liu, Y.; Zhou, K.; Tong, X. Computational Design and Fabrication of Soft Pneumatic Objects with Desired Deformations. ACM Trans. Graph. 2017, 36, 1–12. [Google Scholar] [CrossRef]
  4. Webster III, R.J.; Jones, B.A. Design and kinematic modeling of constant curvature continuum robots: A review. Int. J. Robot. Res. 2010, 29, 1661–1683. [Google Scholar] [CrossRef]
  5. Trivedi, D.; Lotfi, A.; Rahn, C.D. Geometrically exact models for soft robotic manipulators. IEEE Trans. Robot. 2008, 24, 773–780. [Google Scholar] [CrossRef]
  6. Chen, Y.; Xia, Z.; Zhao, Q. Optimal Design of Soft Pneumatic Bending Actuators Subjected to Design-Dependent Pressure Loads. IEEE/ASME Trans. Mechatronics 2019, 24, 2873–2884. [Google Scholar] [CrossRef]
  7. Caasenbrood, B.; Pogromsky, A.; Nijmeijer, H. A Computational Design Framework for Pressure-Driven Soft Robots through Nonlinear Topology Optimization. In Proceedings of the 2020 3rd IEEE International Conference on Soft Robotics (RoboSoft), New Haven, CT, USA, 15 May–15 July 2020; pp. 633–638. [Google Scholar] [CrossRef]
  8. Grazioso, S.; Di Gironimo, G.; Siciliano, B. A Geometrically Exact Model for Soft Continuum Robots: The Finite Element Deformation Space Formulation. Soft Robot. 2019, 6, 790–811. [Google Scholar] [CrossRef]
  9. Runge, G.; Raatz, A. A framework for the automated design and modelling of soft robotic systems. CIRP Ann. 2017, 66, 9–12. [Google Scholar] [CrossRef]
  10. Pinskier, J.; Howard, D. From Bioinspiration to Computer Generation: Developments in Autonomous Soft Robot Design. Adv. Intell. Syst. 2022, 4, 2100086. [Google Scholar] [CrossRef]
  11. Wang, Z.; Chen, M.Z.Q.; Yi, J. Soft robotics for engineers. HKIE Trans. 2015, 22, 88–97. [Google Scholar] [CrossRef]
  12. Navas, E.; Fernandez, R.; Sepulveda, D.; Armada, M.; Gonzalez-De-santos, P. Soft grippers for automatic crop harvesting: A review. Sensors 2021, 21, 2689. [Google Scholar] [CrossRef]
  13. Tolley, M.T.; Shepherd, R.F.; Mosadegh, B.; Galloway, K.C.; Wehner, M.; Karpelson, M.; Wood, R.J.; Whitesides, G.M. A Resilient, Untethered Soft Robot. Soft Robot. 2014, 1, 213–223. [Google Scholar] [CrossRef] [Green Version]
  14. Paul, C. Morphological computation: A basis for the analysis of morphology and control requirements. Robot. Auton. Syst. 2006, 54, 619–630. [Google Scholar] [CrossRef]
  15. Howard, D.; Eiben, A.E.; Kennedy, D.F.; Mouret, J.B.; Valencia, P.; Winkler, D. Evolving embodied intelligence from materials to machines. Nat. Mach. Intell. 2019, 1, 12–19. [Google Scholar] [CrossRef] [Green Version]
  16. Xavier, M.S.; Tawk, C.D.; Zolfagharian, A.; Pinskier, J.; Howard, D.; Young, T.; Lai, J.; Harrison, S.M.; Yong, Y.K.; Bodaghi, M.; et al. Soft Pneumatic Actuators: A Review of Design, Fabrication, Modeling, Sensing, Control and Applications. IEEE Access 2022, 10, 59442–59485. [Google Scholar] [CrossRef]
  17. Ilievski, F.; Mazzeo, A.D.; Shepherd, R.F.; Chen, X.; Whitesides, G.M. Soft Robotics for Chemists. Angew. Chem. 2011, 50, 1890–1895. [Google Scholar] [CrossRef]
  18. Mosadegh, B.; Polygerinos, P.; Keplinger, C.; Wennstedt, S.; Shepherd, R.F.; Gupta, U.; Shim, J.; Bertoldi, K.; Walsh, C.J.; Whitesides, G.M. Pneumatic networks for soft robotics that actuate rapidly. Adv. Funct. Mater. 2014, 24, 2163–2170. [Google Scholar] [CrossRef] [Green Version]
  19. Marchese, A.D.; Onal, C.D.; Rus, D. Autonomous Soft Robotic Fish Capable of Escape Maneuvers Using Fluidic Elastomer Actuators. Soft Robot. 2014, 1, 75–87. [Google Scholar] [CrossRef] [Green Version]
  20. Schiller, L.; Seibel, A.; Schlattmann, J. Toward a Gecko-Inspired, Climbing Soft Robot. Front. Neurorobot. 2019, 13, 106. [Google Scholar] [CrossRef] [Green Version]
  21. Cheney, N.; MacCurdy, R.; Clune, J.; Lipson, H. Unshackling Evolution: Evolving Soft Robots with Multiple Materials and a Powerful Generative Encoding. SIGEVOlution 2014, 7, 11–23. [Google Scholar] [CrossRef]
  22. Runge, G.; Wiese, M.; Günther, L.; Raatz, A. A framework for the kinematic modeling of soft material robots combining finite element analysis and piecewise constant curvature kinematics. In Proceedings of the 2017 3rd International Conference on Control, Automation and Robotics (ICCAR), Nagoya, Japan, 22–24 April 2017; pp. 7–14. [Google Scholar] [CrossRef]
  23. Hiller, J.; Lipson, H. Automatic Design and Manufacture of Soft Robots. IEEE Trans. Robot. 2012, 28, 457–466. [Google Scholar] [CrossRef]
  24. Lehman, J.; Clune, J.; Misevic, D.; Adami, C.; Altenberg, L.; Beaulieu, J.; Bentley, P.J.; Bernard, S.; Beslon, G.; Bryson, D.M.; et al. The Surprising Creativity of Digital Evolution: A Collection of Anecdotes from the Evolutionary Computation and Artificial Life Research Communities. Artif. Life 2020, 26, 274–306. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Hamza, K.; Saitou, K. Design Optimization of Vehicle Structures for Crashworthiness Using Equivalent Mechanism Approximations. J. Mech. Des. 2004, 127, 485–492. [Google Scholar] [CrossRef]
  26. Conradie, N.T. A Scale-Invariant Generative Design Process for 2D Soft Robot Actuators. Masters Thesis, Stellenbosch University, Stellenbosch, South Africa, 2021. [Google Scholar]
  27. Ellis, D.R.; Venter, M.P.; Venter, G. Generative Design Procedure for Embedding Specified Planar Behavior in Modular Soft Pneumatic Actuators. Soft Robot. 2022, 9, 552–561. [Google Scholar] [CrossRef]
  28. Braden, B. The Surveyor’s Area Formula. Coll. Math. J. 1986, 17, 326–337. [Google Scholar] [CrossRef]
  29. Joubert, I.J. Generative Design Using Lindenmayer-Systems and Numerical Optimization. Masters Thesis, Stellenbosch University, Stellenbosch, South Africa, 2020. [Google Scholar]
  30. Hexagon, A.B. Marc 2021.4, Nexus Portal, Marc 2021.4 Volume A: Theory and User Information; Hexagon AB: Irvine, CA, USA, 2023. [Google Scholar]
  31. Vanderplaats Research & Development, Inc. DOT Design Optimisation Tools Users Manual, Version 7; Vanderplaats Research & Development, Inc.: Colorado Springs, CO, USA, 2017. [Google Scholar]
Figure 1. Hierarchical framework made up of multiple hierarchical steps. The complete design is broken down into simpler design problems in the forward pass. The simpler design problems are then solved and reassembled in the backward pass. The directions of the forward and backward passes are represented by the dashed arrows.
Figure 1. Hierarchical framework made up of multiple hierarchical steps. The complete design is broken down into simpler design problems in the forward pass. The simpler design problems are then solved and reassembled in the backward pass. The directions of the forward and backward passes are represented by the dashed arrows.
Mca 28 00047 g001
Figure 2. Single hierarchical step representing an optimisation problem. Multiple candidate solutions can be generated, which may be used in succeeding HSs.
Figure 2. Single hierarchical step representing an optimisation problem. Multiple candidate solutions can be generated, which may be used in succeeding HSs.
Mca 28 00047 g002
Figure 3. A complete hierarchical framework with figure representations of each hierarchical step in the framework.
Figure 3. A complete hierarchical framework with figure representations of each hierarchical step in the framework.
Mca 28 00047 g003
Figure 4. Custom reduced order representation of an sPN undergoing actuation.
Figure 4. Custom reduced order representation of an sPN undergoing actuation.
Mca 28 00047 g004
Figure 5. Unit definitions. (a) The chamber layout within a unit. The dotted lines represent the boundaries of each chamber node. (b) Actuation pattern. Both chambers start and end at the off position, and this allows for cyclic actuation. Each coloured line corresponds to the actuation of the matching coloured chamber.
Figure 5. Unit definitions. (a) The chamber layout within a unit. The dotted lines represent the boundaries of each chamber node. (b) Actuation pattern. Both chambers start and end at the off position, and this allows for cyclic actuation. Each coloured line corresponds to the actuation of the matching coloured chamber.
Mca 28 00047 g005
Figure 6. Reduction of a full-length actuator to its simplest representation. All units with unique boundary conditions remain in the reduced representation. Units with similar boundary conditions are condensed into a single unit representative of the group [29].
Figure 6. Reduction of a full-length actuator to its simplest representation. All units with unique boundary conditions remain in the reduced representation. Units with similar boundary conditions are condensed into a single unit representative of the group [29].
Mca 28 00047 g006
Figure 7. Three-unit finite element result showing the optimised chamber designs. The red line shows the enclosed area, which was maximised during the optimisation. For the three-element method, the midpoint of the right-hand side of the middle element is used as the tip to avoid the bulging of the tip having too large of an impact on the optimisation.
Figure 7. Three-unit finite element result showing the optimised chamber designs. The red line shows the enclosed area, which was maximised during the optimisation. For the three-element method, the midpoint of the right-hand side of the middle element is used as the tip to avoid the bulging of the tip having too large of an impact on the optimisation.
Mca 28 00047 g007
Figure 8. Result from the three-element optimisation design expanded to the full nine-unit actuator. FE mesh is shown in black with the tip displacements traced in red.
Figure 8. Result from the three-element optimisation design expanded to the full nine-unit actuator. FE mesh is shown in black with the tip displacements traced in red.
Mca 28 00047 g008
Figure 9. The different design method representations. Red dots represent design points, blue squares represent repeated points, the grey background shows FE simulated, and the wavy background shows non-FE simulated. (a) A complete domain where every chamber node is a design point and the entire actuator is simulated. This is a direct optimisation. (b) The hierarchical method is used where repetition is applied and only 3 elements are simulated.
Figure 9. The different design method representations. Red dots represent design points, blue squares represent repeated points, the grey background shows FE simulated, and the wavy background shows non-FE simulated. (a) A complete domain where every chamber node is a design point and the entire actuator is simulated. This is a direct optimisation. (b) The hierarchical method is used where repetition is applied and only 3 elements are simulated.
Mca 28 00047 g009
Figure 10. Run time and objective function comparison between the hierarchical framework method (HFm) and the direct method (Dm). R represents Restart. The objective function is defined in Equation (3). The solid line shows the results when both methods have the same initial conditions which were a single 1 × 1 mm void at the centre of each chamber domain. All other optimisations (shown as dotted lines) have random starting points. Best solution values, A = (1.07, 843 ), and B = (50.65, 342 ).
Figure 10. Run time and objective function comparison between the hierarchical framework method (HFm) and the direct method (Dm). R represents Restart. The objective function is defined in Equation (3). The solid line shows the results when both methods have the same initial conditions which were a single 1 × 1 mm void at the centre of each chamber domain. All other optimisations (shown as dotted lines) have random starting points. Best solution values, A = (1.07, 843 ), and B = (50.65, 342 ).
Mca 28 00047 g010
Figure 11. An expanded hierarchical framework with examples of additional HSs in level 2 and level 3.
Figure 11. An expanded hierarchical framework with examples of additional HSs in level 2 and level 3.
Mca 28 00047 g011
Table 1. Hierarchical case study summary.
Table 1. Hierarchical case study summary.
Level 1Level 2Level 3
Forward PassDesired designCreate reduced modelChamber design
Backward PassCascade to full lengthOptimise 3 element modelChamber design
Table 2. Hierarchical framework forward pass properties. Properties were chosen based on similar studies. Material [16], Geometry [6,26,27], Pressure [6,18,27].
Table 2. Hierarchical framework forward pass properties. Properties were chosen based on similar studies. Material [16], Geometry [6,26,27], Pressure [6,18,27].
PropertyValue/Specification
MaterialMold star 30
Pressure0 to 15 kPa
Internal Chambers2
Unit dimensions80 × 80 × 2 mm
Number of units9
Minimum sidewall thickness2 mm
Table 3. Breakdown of different optimisation problems in each method.
Table 3. Breakdown of different optimisation problems in each method.
MethodDesign VariablesConstraintsSide ConstraintsFE Size
Hierarchical Framework16232 1 3
Direct method1442288Full
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ligthart, P.F.; Venter, M.P. A Hierarchical Design Framework for the Design of Soft Robots. Math. Comput. Appl. 2023, 28, 47. https://doi.org/10.3390/mca28020047

AMA Style

Ligthart PF, Venter MP. A Hierarchical Design Framework for the Design of Soft Robots. Mathematical and Computational Applications. 2023; 28(2):47. https://doi.org/10.3390/mca28020047

Chicago/Turabian Style

Ligthart, Philip Frederik, and Martin Philip Venter. 2023. "A Hierarchical Design Framework for the Design of Soft Robots" Mathematical and Computational Applications 28, no. 2: 47. https://doi.org/10.3390/mca28020047

APA Style

Ligthart, P. F., & Venter, M. P. (2023). A Hierarchical Design Framework for the Design of Soft Robots. Mathematical and Computational Applications, 28(2), 47. https://doi.org/10.3390/mca28020047

Article Metrics

Back to TopTop