Skip to main content

Thank you for visiting You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Bubble casting soft robotics


Inspired by living organisms, soft robots are developed from intrinsically compliant materials, enabling continuous motions that mimic animal and vegetal movement1. In soft robots, the canonical hinges and bolts are replaced by elastomers assembled into actuators programmed to change shape following the application of stimuli, for example pneumatic inflation2,3,4,5. The morphing information is typically directly embedded within the shape of these actuators, whose assembly is facilitated by recent advances in rapid prototyping techniques6,7,8,9,10,11. Yet, these manufacturing processes have limitations in scalability, design flexibility and robustness. Here we demonstrate a new all-in-one methodology for the fabrication and the programming of soft machines. Instead of relying on the assembly of individual parts, our approach harnesses interfacial flows in elastomers that progressively cure to robustly produce monolithic pneumatic actuators whose shape can easily be tailored to suit applications ranging from artificial muscles to grippers. We rationalize the fluid mechanics at play in the assembly of our actuators and model their subsequent morphing. We leverage this quantitative knowledge to program these soft machines and produce complex functionalities, for example sequential motion obtained from a monotonic stimulus. We expect that the flexibility, robustness and predictive nature of our methodology will accelerate the proliferation of soft robotics by enabling the assembly of complex actuators, for example long, tortuous or vascular structures, thereby paving the way towards new functionalities stemming from geometric and material nonlinearities.

This is a preview of subscription content

Access options

Buy article

Get time limited or full article access on ReadCube.


All prices are NET prices.

Fig. 1: From flow to programmed actuation.
Fig. 2: Flow determines form.
Fig. 3: Pneumatic actuation produces curvature.
Fig. 4: Programmed mechanics from curvature.

Data availability

The data used in this study are available from the corresponding authors upon reasonable request. Source data are provided with this paper.

Code availability

The codes that support the findings of this study are available from the corresponding authors upon reasonable request.


  1. 1.

    Whitesides, G. M. Soft robotics. Angew. Chem. Int. Ed. 57, 4258–4273 (2018).

    CAS  Google Scholar 

  2. 2.

    Polygerinos, P., Wang, Z., Galloway, K. C., Wood, R. J. & Walsh, C. J. Soft robotic glove for combined assistance and at-home rehabilitation. Rob. Autom. Syst. 73, 135–143 (2015).

    Google Scholar 

  3. 3.

    Shepherd, R. F. et al. Multigait soft robot. Proc. Natl Acad. Sci. USA 108, 20400–20403 (2011).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  4. 4.

    Martinez, R. V. et al. Robotic tentacles with three-dimensional mobility based on flexible elastomers. Adv. Mater. 25, 205–212 (2013).

    CAS  PubMed  Google Scholar 

  5. 5.

    Roche, E. T. et al. A bioinspired soft actuated material. Adv. Mater. 26, 1200–1206 (2014).

    CAS  PubMed  Google Scholar 

  6. 6.

    Mosadegh, B. et al. Pneumatic networks for soft robotics that actuate rapidly. Adv. Funct. Mater. 24, 2163–2170 (2014).

    CAS  Google Scholar 

  7. 7.

    Yap, H. K., Ng, H. Y. & Yeow, C.-H. High-force soft printable pneumatics for soft robotic applications. Soft Robot. 3, 144–158 (2016).

    Google Scholar 

  8. 8.

    Wehner, M. et al. An integrated design and fabrication strategy for entirely soft, autonomous robots. Nature 536, 451–455 (2016).

    ADS  CAS  PubMed  Google Scholar 

  9. 9.

    Zhao, H., Li, Y., Elsamadisi, A. & Shepherd, R. Scalable manufacturing of high force wearable soft actuators. Extreme Mech. Lett. 3, 89–104 (2015).

    Google Scholar 

  10. 10.

    Paek, J., Cho, I. & Kim, J. Microrobotic tentacles with spiral bending capability based on shape-engineered elastomeric microtubes. Sci. Rep. 5, 10768 (2015).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  11. 11.

    Jin, L., Forte, A. E., Deng, B., Rafsanjani, A. & Bertoldi, K. Kirigami-inspired inflatables with programmable shapes. Adv. Mater. 32, 2001863 (2020).

    CAS  Google Scholar 

  12. 12.

    Ilievski, F., Mazzeo, A. D., Shepherd, R. F., Chen, X. & Whitesides, G. M. Soft robotics for chemists. Angew. Chem. Int. Ed. 50, 1890–1895 (2011).

    CAS  Google Scholar 

  13. 13.

    Marchese, A. D., Onal, C. D. & Rus, D. Autonomous soft robotic fish capable of escape maneuvers using fluidic elastomer actuators. Soft Robot. 1, 75–87 (2014).

    PubMed  PubMed Central  Google Scholar 

  14. 14.

    Yang, D. et al. Buckling pneumatic linear actuators inspired by muscle. Adv. Mater. Technol. 1, 1600055 (2016).

    Google Scholar 

  15. 15.

    Overvelde, J. T. B., Kloek, T., D’haen, J. J. A. & Bertoldi, K. Amplifying the response of soft actuators by harnessing snap-through instabilities. Proc. Natl Acad. Sci. USA 112, 10863–10868 (2015).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  16. 16.

    Hawkes, E. W., Blumenschein, L. H., Greer, J. D. & Okamura, A. M. A soft robot that navigates its environment through growth. Sci. Robot. 2, aan3028 (2017).

    Google Scholar 

  17. 17.

    Majidi, C. Soft robotics: a perspective—current trends and prospects for the future. Soft Robot. 1, 5–11 (2014).

    Google Scholar 

  18. 18.

    Gladman, A. S., Matsumoto, E. A., Nuzzo, R. G., Mahadevan, L. & Lewis, J. A. Biomimetic 4D printing. Nat. Mater. 15, 413–418 (2016).

    ADS  PubMed  Google Scholar 

  19. 19.

    Kanik, M. et al. Strain-programmable fiber-based artificial muscle. Science 365, 145–150 (2019).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  20. 20.

    William Boley, J. W. et al. Shape-shifting structured lattices via multimaterial 4D printing. Proc. Natl Acad. Sci. USA 116, 20856–20862 (2019).

    PubMed  PubMed Central  Google Scholar 

  21. 21.

    Guseinov, R., McMahan, C., Pérez, J., Daraio, C. & Bickel, B. Programming temporal morphing of self-actuated shells. Nat. Commun. 11, 237 (2020).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  22. 22.

    Kim, Y., Yuk, H., Zhao, R., Chester, S. A. & Zhao, X. Printing ferromagnetic domains for untethered fast-transforming soft materials. Nature 558, 274–279 (2018).

    ADS  CAS  PubMed  Google Scholar 

  23. 23.

    Hu, W., Lum, G. Z., Mastrangeli, M. & Sitti, M. Small-scale soft-bodied robot with multimodal locomotion. Nature 554, 81–85 (2018).

    ADS  CAS  PubMed  Google Scholar 

  24. 24.

    Acome, E. et al. Hydraulically amplified self-healing electrostatic actuators with muscle-like performance. Science 359, 61–65 (2018).

    ADS  CAS  PubMed  Google Scholar 

  25. 25.

    Gorissen, B. et al. Elastic inflatable actuators for soft robotic applications. Adv. Mater. 29, 1604977 (2017).

    Google Scholar 

  26. 26.

    Polygerinos, P. et al. Soft robotics: review of fluid-driven intrinsically soft devices; manufacturing, sensing, control, and applications in human-robot interaction. Adv. Eng. Mater. 19, 1700016 (2017).

    Google Scholar 

  27. 27.

    Siéfert, E., Reyssat, E., Bico, J. & Roman, B. Bio-inspired pneumatic shape-morphing elastomers. Nat. Mater. 18, 24–28 (2019).

    ADS  PubMed  Google Scholar 

  28. 28.

    Vasios, N., Gross, A. J., Soifer, S., Overvelde, J. T. B. & Bertoldi, K. Harnessing viscous flow to simplify the actuation of fluidic soft robots. Soft Robot. 7, 1–9 (2020).

    PubMed  Google Scholar 

  29. 29.

    Jambon-Puillet, E., Piéchaud, M. R. & Brun, P.-T. Elastic amplification of the Rayleigh–Taylor instability in solidifying melts. Proc. Natl Acad. Sci. USA 118, e2020701118 (2021).

    MathSciNet  CAS  PubMed  PubMed Central  Google Scholar 

  30. 30.

    Bretherton, F. P. The motion of long bubbles in tubes. J. Fluid Mech. 10, 166–188 (1961).

    ADS  MathSciNet  MATH  Google Scholar 

  31. 31.

    Aussillous, P. & Quéré, D. Quick deposition of a fluid on the wall of a tube. Phys. Fluids 12, 2367–2371 (2000).

    ADS  CAS  MATH  Google Scholar 

  32. 32.

    Klaseboer, E., Gupta, R. & Manica, R. An extended Bretherton model for long Taylor bubbles at moderate capillary numbers. Phys. Fluids 26, 032107 (2014).

    ADS  Google Scholar 

  33. 33.

    Takagi, D. & Huppert, H. E. Flow and instability of thin films on a cylinder and sphere. J. Fluid Mech. 647, 221–238 (2010).

    ADS  MathSciNet  MATH  Google Scholar 

  34. 34.

    Lee, A. et al. Fabrication of slender elastic shells by the coating of curved surfaces. Nat. Commun. 7, 11155 (2016).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  35. 35.

    Audoly, B. & Pomeau, Y. Elasticity and Geometry: From Hair Curls to the Non-Linear Response of Shells (Oxford Univ. Press, 2010).

  36. 36.

    Gerbode, S. J., Puzey, J. R., McCormick, A. G. & Mahadevan, L. How the cucumber tendril coils and overwinds. Science 337, 1087–1091 (2012).

    ADS  CAS  PubMed  Google Scholar 

  37. 37.

    Miller, J. T., Lazarus, A., Audoly, B. & Reis, P. M. Shapes of a suspended curly hair. Phys. Rev. Lett. 112, 068103 (2014).

    ADS  CAS  PubMed  Google Scholar 

  38. 38.

    Yang, H. D., Greczek, B. T. & Asbeck, A. T. Modeling and analysis of a high-displacement pneumatic artificial muscle with integrated sensing. Front. Robot. AI 5, 136 (2019).

    PubMed  PubMed Central  Google Scholar 

  39. 39.

    Kaufman, J. J. et al. Structured spheres generated by an in-fibre fluid instability. Nature 487, 463–467 (2012).

    ADS  CAS  PubMed  Google Scholar 

  40. 40.

    Marthelot, J., Strong, E. F., Reis, P. M. & Brun, P.-T. Designing soft materials with interfacial instabilities in liquid films. Nat. Commun. 9, 4477 (2018).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  41. 41.

    Zhao, C. et al. Layered nanocomposites by shear-flow-induced alignment of nanosheets. Nature 580, 210–215 (2020).

    ADS  CAS  PubMed  Google Scholar 

  42. 42.

    Nerger, B. A., Brun, P.-T. & Nelson, C. M. Marangoni flows drive the alignment of fibrillar cell-laden hydrogels. Sci. Adv. 6, eaaz7748 (2020).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  43. 43.

    Cybulski, O. et al. Concentric liquid reactors for chemical synthesis and separation. Nature 586, 57–63 (2020).

    ADS  CAS  PubMed  Google Scholar 

Download references


This work was supported by NSF CAREER award (CBET 2042930) and through the Princeton University Materials Research Science and Engineering Center (NSF Grant DMR-1420541) and the Princeton Yang Family and David T. Wilkinson Innovation Funds.

Author information




T.J.J., J.M. and P.-T.B. conceived the project. T.J.J., E.J.-P. and P.-T.B. conducted the experiments and analysed the data. T.J.J. performed the Kirchhoff rod simulations. E.J.-P. performed the finite element simulations. All authors wrote the manuscript.

Corresponding author

Correspondence to P.-T. Brun.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Peer review information Nature thanks Detlef Lohse, Orlin Velev and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Extended data figures and tables

Extended Data Fig. 1 Characterization of elastomers.

a, Oscillatory shear rheology measurements for VPS-08, 16 and 32 (strain 1%, frequencies {1, 2.5, 2.5} Hz). The cross indicates τc. b, Rescaled time varying viscosity for t < τc. The dashed lines are fits of equation (5). c, Dogbone uniaxial tensile tests data for multiple samples of VPS-32, VPS-16 and VPS-08. The black curves are fits of the various constitutive models. d, Table showing rheological and elastic material constants for VPS

Source data.

Extended Data Fig. 2 Cross-section shape.

a, Schematic for the computation of the shape of the liquid bath. b, Meniscus shape profiles obtained from solving equation (9) for various values of \(R/{{\mathscr{l}}}_{{\rm{c}}}\) and hi/R. c, Model shape profiles obtained from matching of the bath solution (equation (9)) to the film thickness solution (equation (2)) for R = 1.6 mm. The solid lines represent the membrane and the dashed lines represent the bath. The color codes the average membrane thickness hf. d, e, Void fraction ϕ =  (1 − hi/R)2 (d) and average membrane thickness hf (e) as a function of the axial position for a metre-long actuator (sample 1). Error bars represent propagated measurement uncertainty. f, Void fraction for three different samples. Uncertainty is shown as box (middle 50%) and whisker (full range) plots overlaid on the markers. g, Membrane thickness for three different samples. The upward (respectively downward) facing triangles represent the maximum thickness \(h(0)\) (minimum thickness \(h({\psi }_{0})\)), while the diamonds indicates the average thickness hf. Uncertainty is shown as box (middle 50%) and whisker (full range) plots overlaid on the markers

Source data.

Extended Data Fig. 3 Finite element simulations.

a, Deformation of the simulated actuator shown in b for \(P\approx 23\) kPa (experimental shape, G = 0.36 MPa, \({J}_{m}=14\)). Inset: view of the deformed cross-section. The color codes the von Mises stress while the black lines show the edges of the undeformed configuration. b, Curvature κ as a function of the inflation pressure P. Circles represent experimental data (VPS-32, R = 1.57 mm, \(h(0)=52\,{\rm{\mu }}{\rm{m}}\), hi/R = 0.34), the solid curve is the finite element simulation using the experimental cross-sectional shape extracted by image analysis, the dashed curve is the simulation using the modeled cross-section using the experimental value of \(h(0)\). Inset: comparison between the experimental and model shapes. c, Curvature κ as a function of the applied pressure P for our model shapes as we vary hf and R (see legend in d). d, Same data as c rescaled according to equation (17), the dashed line is a power law fit

Source data.

Extended Data Fig. 4 Pneumatic bending of actuator.

a, Schematic for the pneumatic inflation of a bubble cast actuator cross-section. b, Schematic for the bending of a bubble cast actuator. c, Actuator curvature \(\kappa \) as a function of applied pressure P for varied length L (i), thickness hf (ii), radius R (iii), and shear modulus G (iv; see legend in Fig. 3a, b). d, Three-dimensional plot of the rescaled curvature \(\kappa R\) as a function of the rescaled pressure \(P/G\) and the rescaled membrane thickness R/hf. Markers show the experimental data while the blue surface is our model (equation (16)) using a Gent constitutive law with \({\psi }_{0}={\rm{\pi }}/4,\,\varphi =0.4225,\,\chi =0.45,\,{J}_{m}=14\)

Source data.

Extended Data Fig. 5 Blocking force elastica.

a, Series of images (i) and calculated elastic curve (ii) for a bubble cast actuator inflated with one end clamped and the other end blocked by a wall (see Extended Data Fig. 7c; scale bar, 1 cm). b, Blocking force F as a function of the inflation pressure P for various actuators (R, G, hf and L have been varied)

Source data.

Extended Data Fig. 6 Fabrication with a twist.

a, c, Two tubes are kept tightly together with a connector (not represented here) that allows the rotation of one tube with respect to the other. Applying the rotation at the very end of the drainage step, that is, around the gelation point where the polymer is still deformable, yet would not flow significantly anymore under the action of gravity results in the actuator shown in b and d. Colours are a guide to eye; only one polymer is used. Rotating the cylinder (as shown in a (respectively c)) rotates the thin membrane and thus the direction of actuation under pressure. In this example, the 90° (180°) rotation produces two curvatures equal in magnitude but in orthogonal planes (resp. with opposite signs).

Extended Data Fig. 7 Experimental set-ups.

a, Experimental set-up for Bretherton-like flow during fabrication. The camera captures the bubble front velocity U to determine the Capillary number Ca. b, Experimental set-up for inflation-bending experiments. The pressure sensor records the internal–external pressure difference P while the camera records the shape of the actuator. c, Experimental set-up for blocking force–pressure experiments. The pressure sensor records the internal–external pressure difference P while the load cell measures the force F. d, Experimental set-up for force-elongation experiments. The pressure sensor records the internal–external pressure difference P while the Instron measures the force F and displacement \({\ell }\).

Supplementary information

Supplementary Information

Supplementary text, Supplementary equations and Supplementary references

Supplementary Video 1

Bubble casting fabrication method. A tortuous channel in an acrylic mould is filled with a curing VPS melt. An elongated bubble is then injected through the channel. The VPS film is then given time to cure as gravity causes drainage. The resulting soft actuator bends when inflated, demonstrating the ability to grasp a blackberry.

Supplementary Video 2

High aspect ratio pneumatic muscles. Comparison between an inflated actuator at pressures P and a simulation of a Kirchhoff rod with changing natural curvature κ according to equation (3) of the main text. A human-scale actuator demonstrates coiling behaviour when inflated with air. An artificial limb mimics a human arm with a contractile force. Actuators with diameters of 12.8 mm and 1.0 mm are inflated to lift a water bottle and paper clip. A high-aspect-ratio curved actuator wraps itself around a raspberry and lifts it without any damage.

Supplementary Video 3

Sequential actuation of digits. The programming logic is depicted via schematic representation. A bubble is injected through the channel digits (1,2,3,4) in a sequential manner at various waiting times τw,i. The evolving viscosity increases as the polymer melt cross-links and diverges at the curing time τc as the melt cures into an elastic solid. Injecting air into the resulting actuator causes a sequential bending of the actuator digits.

Supplementary Video 4

Morphing to complex shapes. An actuator with a tortuous path is inflated in water. The actuator bends out of plane to a curve that lives on the surface of a sphere. A spiral-shaped actuator is attached to a circular membrane, transforming local curvature change of the rod to a global bending of the membrane. An actuator is attached to the edge of a membrane to perform self-folding of a five-sided box. Two branched actuators are bonded to either side of a membrane in the shape of a caudal fin. The fin actuators are alternatively inflated in a water bath with tracer particles. During fabrication a tubular mould is rotated 180 degrees at the gelation point. The resulting actuator bends in two different directions. During fabrication a three sectioned tubular mould is rotated at 180 degrees and 90 degrees at the gelation point. The resulting actuator bends in three directions and in two planes.

Supplementary Video 5

A soft machine to fetch in a constrained space. A modular soft machine made of two actuators is programmed to grip and then retrieve a ball from a cylindrical vessel using a single pressure source.

Supplementary Video 6

Programming curvature variation. A rectilinear actuator, built by bubble casting, is programmed to have three sections of varying curvatures (curvature ratio 4:2:1). The coding is accomplished by injecting the bubble through the tube in a step wise sequence such that the final film thickness is predicted. The actuator is inflated on a water bath forming a spiral with appropriate curvatures.

Source data

Rights and permissions

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Jones, T.J., Jambon-Puillet, E., Marthelot, J. et al. Bubble casting soft robotics. Nature 599, 229–233 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI:


By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.


Quick links

Nature Briefing

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing