Chiral topographic instability in shrinking spheres

Nature Computational Science volume 2, pages 632–640 (2022)Cite this article

5125 Accesses

3 Citations

66 Altmetric

Metrics details

Many biological structures exhibit intriguing morphological patterns adapted to environmental cues, which contribute to their important biological functions and also inspire material designs. Here, we report a chiral wrinkling topography in shrinking core–shell spheres, as observed in excessively dehydrated passion fruit and experimentally demonstrated in silicon core–shells under air extraction. Upon shrinkage deformation, the surface initially buckles into a buckyball pattern (periodic hexagons and pentagons) and then transforms into a chiral mode. The neighbouring chiral cellular patterns can further interact with each other, resulting in secondary symmetry breaking and the formation of two types of topological network. We develop a core–shell model and derive a universal scaling law to understand the underlying morphoelastic mechanism and to effectively describe and predict such chiral symmetry breaking far beyond the critical instability threshold. Moreover, we show experimentally that the chiral characteristic adapted to local perturbation can be harnessed to effectively and stably grasp small-sized objects of various shapes and made of different stiff and soft materials. Our results not only reveal chiral instability topographies, providing fundamental insights into the surface morphogenesis of the deformed core–shell spheres that are ubiquitous in the real world, but also demonstrate potential applications of adaptive grasping based on delicate chiral localization.

Morphological pattern formation across length scales is energetically favourable for thin-walled living matter such as fruits1,2, vegetables3, leaves4,5,6, embryos7, organs8, tumours9 and brains10, where spontaneous symmetry breaking during growth or dehydration is normally considered to be a crucial factor in their complex wrinkling topography6,11,12. For example, pollen grains of angiosperm flowers exhibit self-folding when exposed to a dry environment to prevent further desiccation13. Growth-induced residual stress accumulates during tumour progression, driving the global buckling collapse of blood and lymphatic vessels, which makes the vascular delivery of anticancer drugs ineffective9. Symmetry breaking in evolving wrinkling patterns during brain development results in the thickness difference between gyri and sulci, which is closely linked to neurodevelopment disorders such as lissencephaly, polymicrogyria, autism spectrum disorders and schizophrenia14. In terms of its practical use, symmetry breaking in the formation of surface morphology patterns has found ever-increasing applications in various fields, such as micro/nanofabrication of flexible electronic devices15,16, surface self-cleaning and anti-fouling17, synthetic camouflaging skins18, shape-morphing soft actuators19 and adaptive aerodynamic drag control20. The precise prediction, control and manipulation of reversible instability morphologies would be key for relevant applications.

Prior works3,12,21,22,23 on morphological pattern formation in stressed spherical core–shells, a typical structure omnipresent in nature and industrial technologies, have demonstrated a variety of intriguing topographies such as dimple, buckyball and labyrinth modes. Here, we report a chiral instability topography in core–shell spheres. We observed that a drying passion fruit (Passiflora edulia Sims) initially buckles into a periodic buckyball pattern consisting of hexagons and pentagons, evolving into a chiral mode, and forms intriguing chiral topological networks upon excessive shrinkage (Fig. 1). Inspired by this natural phenomenon, we explored, both theoretically and experimentally, the morphological pattern formation and evolution of highly deformed core–shell spheres, especially the emergence of a chiral pattern and chiral ridge networks with symmetry breaking at the advanced bifurcation. We established a mathematical model and a scaling law to capture the chiral instability of core–shell spheres and explored a potential application of perturbation-adaptive chiral localization.

a–h, Natural observations (a–d) and model predictions (e–h) on day 1 (a,e), day 2 (b,f), day 4 (c,g) and day 7 (d,h). Upon shrinkage, the core–shell spheres first buckle into a buckyball pattern (periodic hexagons and pentagons in b and f) and then transform to a chiral ridge (g) and eventually to a ridge network (h) with the coalescence of neighbouring chiral ridges. The core experiences isotropic shrinking (Supplementary Sections I and II and Video 1).

To understand the underlying mechanism and to effectively predict the morphogenesis process, we consider an elastic spherical shell supported by a soft core. Upon shrinkage, the shell buckles elastically to relieve the compressive stress while the core concurrently deforms to maintain perfect bonding at the interface. In shallow shell theory24, the coordinates of the core–shell system can be Cartesian in a tangent plane (or curvilinear and orthogonal). This framework can only describe a part of the spherical geometry (Extended Data Fig. 1), but it is competent here for theoretical analyses. The thickness of the surface layer is denoted by hf, while the radius of the system is represented by R. The Young's modulus and Poisson's ratio of the surface layer are denoted by Ef and νf, respectively, while Es and νs are the corresponding material properties of the soft core. The elastic strain energy Πf in the shell can be written as the sum of the bending energy Πben and membrane energy Πmem thus

where \(D={E}_{\mathrm{f}}{h}_{\mathrm{f}}^{3}/[12(1-{\nu }_{\mathrm{f}}^{2})]\) and \({J}_{\mathrm{f}}={E}_{\mathrm{f}}{h}_{\mathrm{f}}/(1-{\nu }_{\mathrm{f}}^{2})\) stand for, respectively, the flexural and extensional rigidities of the shell, and \({\overline{{{{\mathbf{L}}}}}}_{\mathrm{f}}\) represents the dimensionless elastic matrix. The membrane strain tensor and curvature tensor are denoted by γ and K, respectively. The elastic behaviour of the core can be described by a Winkler-type foundation25,26 as

in which \({K}_{\mathrm{s}}={\overline{E}}_{\mathrm{s}}\sqrt{{p}^{2}+{q}^{2}}/2R\) denotes the stiffness of the core23,27, w stands for deflection, \({\overline{E}}_{\mathrm{s}}={E}_{\mathrm{s}}/(1-{\nu }_{\mathrm{s}}^{2})\), and p and q represent the wavenumbers along the latitude and longitude directions, respectively.

The critical buckling of a core–shell sphere upon shrinkage is analogous to the hydrostatic instability of a spherical shell where an isotropic stress state remains in the pre-buckling stage, that is, σαβδαβ = −σ, in which δαβ is the Kronecker delta, σ denotes the external hydrostatic pressure and the Greek indices α and β take values in {1, 2}. According to Koiter's theory24, elastic stability is primarily determined by the second variation of the total potential energy (Πt = Πf + Πs), and one obtains the equilibrium partial differential equations by using the divergence theorem,

where a comma in a subscript denotes a partial derivative. As an ansatz, we consider the following forms for the displacements in the critical buckling state:

in which A, B and C refer to the amplitudes of waves. Substituting equations (4) into equations (3) and minimizing with respect to k = p2 + q2, one obtains the critical conditions for the onset of wrinkling:

where kcr, σcr and ℓcr denote, respectively, the critical wavenumber, the compressive stress and the wavelength, \(c=\sqrt{3(1-{\nu }_{\mathrm{f}}^{2})}\). Here, we define a key dimensionless parameter \({C}_{\mathrm{s}}=({E}_{\mathrm{s}}/{E}_{\mathrm{f}}){(R/{h}_{\mathrm{f}})}^{3/2}\) that characterizes the stiffness ratio of core–shells and the geometric curvature to classify pattern selection. Once the critical wavenumber kcr is solved, the theoretical buckling stress and wavelength can be calculated (Fig. 2a). During the natural dehydration process of passion fruit, the moduli of both the surface layer and the soft core may become larger (meaning that the surface layer and the core become stiffer), but we observed that the wrinkling wavelength in experiments (Fig. 1 and Supplementary Video 1) remains almost unchanged, and this critical wavelength ℓcr has some inherent (yet implicit) relation with the modulus ratio Es/Ef (equation (5)). Therefore, it is reasonable to approximate in the calculation that the modulus ratio Es/Ef remains relatively constant upon dehydration. Note that, although both natural and numerical observations (Fig. 1b,f) show that the buckyball pattern consisting of hexagons and pentagons covers the whole sphere (non-developable surface), the prevailing buckling mode in core–shell spheres is hexagonal. Also within the shallow shell framework (a part of sphere)24, it is an analytical challenge to apply both hexagons and pentagons to describe the entire spherical surface. Hence, we assume this dominant hexagonal mode (displacement field) in equation (4), and the critical wrinkling condition based on our theory shows good agreement with numerical simulations. Equation (5), in fact, covers the classical buckling case of a spherical shell without a core (Ks = 0), for which there are explicit solutions for the critical threshold, that is, σ0 = Efhf/cR, k0 = 2cR/hf and \({\ell }_{0}=\uppi \sqrt{2R{h}_{\mathrm{f}}/c}\).

a, The critical hexagonal wrinkling wavelength ℓcr as a function of the dimensionless parameter \({C}_{\mathrm{s}}=({E}_{\mathrm{s}}/{E}_{\mathrm{f}}){(R/{h}_{\mathrm{f}})}^{3/2}\) that characterizes the modulus ratio and curvature. b, A scaling law (Methods) for the hexagonal-to-chiral mode transition. Our theoretical predictions agree well with FEM simulations, where C1 denotes the slope.

Source data

Although the critical buckling condition can be predicted analytically by using stability analysis, the secondary bifurcation with the hexagonal-to-chiral mode transition in the post-buckling stage remains a theoretical challenge. Here, we derived a scaling law to provide further insight into such chiral symmetry breaking far beyond the critical threshold (Methods). We assumed that each Y-shaped ridge in the wrinkling hexagons can be regarded as a bilayer system and thus that the chiral ridge instability of core–shell spheres can be simplified as the buckling of bilayered plates under compression. Minimization of the system energy leads to chiral strains that obey the linear relation in Fig. 2b, confirmed by numerical simulations.

To trace the whole post-buckling topographic evolution, we applied the finite element method (FEM) by accounting for various geometric and material parameters (Supplementary Section II). The main challenge lies in the solution of nonlinear equations, since multiple solution branches in the post-buckling regime can be connected via multiple bifurcations. Moreover, for instabilities that are extremely localized (for example, the ridge network shown in Fig. 1c,d), there must exist a local transfer of elastic strain energy from one part of the system to the neighbouring regions, and global solution methods may encounter difficulties in convergence. To solve this difficulty, we implemented a pseudodynamic algorithm by introducing velocity-dependent damping and inertial terms, which can be naturally viewed as a perturbation to allow the calculation to pass through the unstable transitions and to trigger chiral symmetry breaking (Methods). The bifurcation portraits of the dimensionless deflection ∣w∣/hf for various core–shell spheres with different Cs upon shrinkage are plotted in Fig. 3. Periodic buckyball (with hexagons prevailing) wrinkling patterns with supercritical bifurcation emerge initially at the critical thresholds. Upon further shrinkage, hexagonal-to-chiral mode transitions occur, where Y-shaped ridges in the wrinkling hexagons may buckle into chiral ridges. Neighbouring chiral cellular modes can further interact with each other to form two types of topological network. While symmetry is eventually broken with further shrinkage, leading to universal hexagonal-to-chiral mode transitions, different Cs values result in different critical thresholds and wavelengths for the buckyball (with hexagon dominating) buckling mode.

a–f, Diagrams for Cs values of 12.7 (a), 9.09 (b), 7.07 (c), 3.98 (d), 3.18 (e) and 2.55 (f), showing the buckyball pattern (with hexagons prevailing) (i) and chiral ridge networks (ii and iii). Excess shrinkage leads to advanced symmetry breaking of the buckyball mode, transforming into the chiral mode and the chiral ridge network eventually.

Source data

Guided by this theoretical understanding, we next designed a demonstrative experiment to harness such an instability mechanism to achieve pattern tunability, by using liquid silicone that can solidify into any desired shape in a well-designed mould. We made a spherical shell with a hexagonal pattern on the surface, a cavity and a small hole (diameter ~4 mm) for air extraction to induce shrinkage (Methods). Since silicone has a much lower elastic modulus than passion fruit, the smooth shell structure does not buckle into hexagonal patterns (cannot reach the advanced bifurcation range shown in Fig. 3) but exhibits global deformation upon pressure loading condition by air extraction (Methods and Supplementary Video 5). To focus on the chiral bifurcation and to facilitate instability morphology control at this bifurcation, we fabricated artificial hexagonal patterns on the shell surface. We extracted air slowly (~2 mL s−1) from the sample to control the pressure (~10 kPa) so that a state of homogeneous compression could be perfectly achieved. Notably, these well-designed hexagonal networks on the surface of the sample buckle into chiral patterns (Fig. 4a–d and Supplementary Video 2), analogous to the observation of highly dehydrated passion fruits and model predictions (Fig. 1). Furthermore, we can flexibly control the position of local chiral networks by imposing external perturbation as illustrated in Fig. 4e–h (Methods and Supplementary Video 3), consistent with FEM simulations in Fig. 4i–l. These experiments not only demonstrate a hexagonal-to-chiral mode transition, consistent with our theoretical predictions, but also shed light on rational designs of controllable chiral patterns.

a–d, The experimental formation of a chiral ridge network with continuous air extraction, showing the hexagonal-to-chiral mode transition with increasing shrinkage of core–shells (Supplementary Video 2). e–l, The localization of tunable chiral networks on curved surfaces (Supplementary Video 3) triggered by a perturbation (poke by a rod) in experiments (e–h), consistent with numerical simulations (i–l).

Based on these insights, we show that this perturbation-induced chiral instability can be harnessed to effectively and stably grasp small-sized objects with different geometries and made of different stiff or soft materials. The object to grasp acts as a local perturbation when in contact with the hexagonal-patterned shell and is then adaptively locked by the induced local chiral networks. Similar to the aforementioned experimental setup, we fabricated a hemispherical shell with a hexagonal surface pattern as the main body of the gripper. A small hole was made at the bottom of the cap for air extraction. Then, the whole gripper was fixed onto a lifting frame to steadily control the movement. When the curved hemispherical cap touches the target, the contact perturbation-induced symmetry breaking triggers chiral network localization. The chiral pattern and the interface friction spontaneously adapt to the interactions at the contacting areas, which are naturally influenced by the shape and stiffness of the object, so that different objects can be grasped by this smart locking together with air extraction (Fig. 5, Supplementary Fig. 4 and Video 4). When we restored the pressure difference, that is, inflated the cap cavity, the chiral networks elastically reverted back to hexagons, releasing the grasped object. The contrast experiments showed that the hemispherical caps with a smooth surface (no chiral instability) could not grasp those objects at all (Supplementary Video 5), supporting the critical role of the chiral network localization in the grasping process.

a–j, Grasping of different objects: diamond (a,b), nut (c), screw (d), mung bean (e), soya bean (f), blueberry (g), heart-shaped candy (h), irregular shaped glass (i) and glass ball (j). The chiral deformation enables effective, target-adaptive grasping (Supplementary Video 4).

We have unveiled chiral-mode symmetry breaking during excessive shrinkage of core–shell spheres, which can be formulaically described and precisely predicted by our theories and computations, in good agreement with carefully designed experiments. Beyond the critical buckyball wrinkling, chiral ridges emerge on the curved surfaces upon excess deformation, and the neighbouring chiral cellular Y-shaped modes can further interact with each other to form advanced chiral topological networks. The critical buckyball wrinkling conditions can be obtained analytically by using linear stability analysis, while strong nonlinearity (both geometric and material) in the post-buckling regime of shrinking spheres results in considerable difficulties in the theoretical predictions of advanced bifurcations and their associated morphological patterns. Consequently, theoretical analyses on secondary and multiple bifurcations of chiral instability have to resort to dimensional analysis (scaling law) based on certain simplified models. From the computational standpoint, the major challenge in extremely shrinking spheres at large strain is the solution of highly nonlinear equations. The most classical solution method to solve nonlinear static problems is the path-following continuation technique such as that of Riks, while numerical convergence cannot always be ensured for extreme wrinkling problems upon large deformations, since a large number of solution branches can be connected via multiple bifurcations. This fact motivated us to apply the dynamic relaxation method to leap over some localized energy barriers in the nonlinear evolution paths, while the dynamic method cannot straightforwardly predict subcritical bifurcations and hysteresis. Making progress in both theoretical and computational analyses of multiple bifurcations in highly nonlinear evolution paths might require more advanced mathematical approaches.

Inspired by the chiral instability topography induced by local perturbation, we demonstrated an exemplar application of target-adaptive grasping based on chiral localization, while future work may take advantage of smart active materials such as hard-magnetic soft materials and liquid-crystal elastomers to enhance multifunctional designs under multiphysics stimuli. Our results not only provide physical insights into the wrinkling topography of highly deformed core–shell spheres by a universal law but also pave a promising way for realizing multifunctional surfaces by harnessing fruitful topography on curved geometry.

We carried out dimensional analysis to predict the chiral bifurcation of core–shell spheres (Extended Data Fig. 1) upon dehydration (equivalent to thermal shrinkage). Based on the experimental observations and numerical calculations, we assumed that each cellular ridge before chiral instability can be viewed as a layered plate and thus the chiral bifurcation of a cellular ridge can be simplified as the buckling of a bilayer subject to shrinking strain (Extended Data Fig. 1c). Such a plate-like ridge has length L and thickness t and comprises an upper layer of width hf and a lower layer of width hs. Each layer has a Young's modulus Eζ, Poisson's ratio νζ and bending stiffness \({D}_{\zeta }={E}_{\zeta }{t}^{3}/[12(1-{\nu }_{\zeta }^{2})]\), where ζ is ‘f’ or ‘s’.

The bending energies of the upper and lower layers can be expressed as

where uf and us denote, respectively, the out-of-plane deflection of the upper and lower layers, while Ω1 and Ω2 represent the area of the mid surface of the upper and lower layer, respectively.

As an ansatz, we consider the following forms for the deflections in the chiral buckling state:

where the functions Φf(z) and Φs(z) can be expanded into series of exponential decay functions as

where kfi and ksi are coefficients of the following order:

and the displacement continuity condition is satisfied at the interface of upper and lower layers, that is, Φf(hs) = Φs(hs).

According to equations (8) to (12), one obtains

Substituting equation (13) into equations (6) and (7), the bending energies read

in which \({a}_{1}=\iint {\left[{\sum }_{i}{A}_{\mathrm{f}i}\left({k}_{\mathrm{f}i}\tilde{z}{h}_{\mathrm{f}}\right)\sin \left(\uppi \tilde{y}\right)\right]}^{2}{{{\rm{d}}}}\tilde{y} \, {{{\rm{d}}}}\tilde{z}\), \({a}_{2}=\iint {\left[{\sum }_{i}{A}_{\mathrm{s}i}\left({k}_{\mathrm{s}i}\tilde{z}{h}_{\mathrm{s}}\right)\sin \left(\uppi \tilde{y}\right)\right]}^{2}{{{\rm{d}}}}\tilde{y} \, {{{\rm{d}}}}\tilde{z}\), \(\tilde{y}=y/L\) and \(\tilde{z}=z/{h}_{\zeta }\).

The membrane energy can be determined by the in-plane strains given by (note that, for simplicity, the subscript ζ has been omitted)

where εsh is the thermal shrinking strain, and v and w represent the in-plane displacements in the mid surface along the y and z directions, respectively, the order of which can be determined by minimizing the membrane energy. Consequently, the in-plane displacements in the mid surface can be approximated as v = By and w = Cz, in which B and C refer to the slopes of variation.

The membrane energies of the upper and lower layers can be expressed as

According to equations (8) to (12) and (16) to (18), the membrane energies read

Since the upper and lower layers buckle simultaneously, combining equations (14), (15), (21) and (22) leads to

namely,

Note that a1/a2 is a non-negative constant. Based on calculations and equation (24), the scaling law yields the following explicit form for the chiral shrinking strain εc:

where C1 = 0.029 is a fitting coefficient. The scaling law in equation (25) agrees well with finite element simulations for chiral bifurcation (Fig. 2b).

We performed finite element simulations in commercial software Abaqus based on parameters similar to experimental observations. Since the deformation of core–shell spheres can be large (up to 30% shrinking strain), we applied the widely used hyperelastic neo-Hookean (nHk) constitutive law for both the surface layer and the soft core, while more sophisticated hyperelastic constitutions such as the Mooney–Rivlin (MR) model were also examined but showed trivial quantitative differences that did not change the substantial nonlinear mechanism of the instability problem. The elastic strain energy density function of the nHk model is defined as

in which \({C}_{10}=E/4\left(1+\nu \right)\) and \({D}_{1}=6\left(1-2\nu \right)/E\) are material parameters. The volume change reads \(J=\det ({{{\mathbf{F}}}})\), where F is the deformation gradient tensor. The first strain invariant reads \({I}_{1}={{{\rm{tr}}}}({{{{\mathbf{F}}}}}^{\mathrm{T}}\cdot {{{\mathbf{F}}}})\). We coupled eight-node hexahedral volume (C3D8R) elements for the soft core and thin shell (S4R) elements for the surface layer by using a ‘tie’ constraint at the interface. Mesh convergence was carefully examined for all simulations. The main challenge is the solution of the nonlinear equations, as numerous post-buckling solution branches can be connected via multiple bifurcations23,28. Therefore, we applied the dynamic relaxation method to allow the calculation to pass through the unstable transitions, which introduces velocity-dependent damping (C) and artificial inertial (M) terms into the static equilibrium equation (R(U, λ) = 0), leading to

where R is the residual force, U denotes unknown variables and λ represents an incremental loading parameter. Realistic definitions of mass and damping were not necessary; thus, we set these quantities to obtain optimal convergence of t → U(t) for large values of time t (no physical meaning here). When the model is stable (quasi-static), viscous energy dissipation remains quite small such that the artificial damping does not notably perturb the solution. When the system tends to be dynamically unstable, nodal velocities increase, and thus, part of the elastic strain energy released can be dissipated by the damping. A shrinkage load (equivalent to thermal expansion or residual strain) was applied to the core while the surface layer was loading free, which can be expressed as

where α, ΔT and I stand for the thermal expansion coefficient, temperature change and second-order identity tensor, respectively. The shrinkage load εsh can also be characterized by an isotropic residual strain εsh = εres = −λI. In the numerical calculations shown in Fig. 1e–h, we took R/h = 50 and \({C}_{\mathrm{s}}=({E}_{\mathrm{s}}/{E}_{\mathrm{f}}){(R/{h}_{\mathrm{f}})}^{3/2}=9.09\).

To realize flexible tunability of chiral patterns and to further harness the hexagonal-to-chiral mode transition for achieving smart surfaces, we designed demonstrative experiments based on air extraction from silicon core–shell spheres. The simple experimental system consists of two combined hemispherical caps with a channel connecting the internal cavity and an external tube for air extraction. To achieve a hexagonal network on the surface of the hemispherical cap, we designed a mould with a hexagonal network by applying three-dimensional printing technology. Then, we poured in two-part liquid silicone (Hongyejie Technology Co. Ltd.) in 1:1 mass ratio. Liquid silicone needs to stand for 3 hours at 25 °C to cure fully. To create a cavity in the centre of the sample, we applied a hemispherical lid with a diameter slightly smaller than the outer diameter to cover the bottom of the mould when the liquid silicone was curing. After the liquid silicone had cured and was demoulded, we glued two identical hemispherical caps together. The typical parameters of the samples were an outer diameter of 2R = 70 mm, a diameter of the inner cavity of 2r = 58 mm and a hexagonal cellular length of L = 4.33 mm, height of H = 2.61 mm and thickness of t = 0.75 mm. The experimental procedure to realize functional chiral surfaces is illustrated in Extended Data Fig. 2. The inner cavity of the samples was pumped out and depressurized to create a state of homogeneous shrinkage. To demonstrate the effects of shrinkage on the hexagonal-to-chiral mode transition, we slowly exhausted the air in the samples to mimic dehydration-induced shrinkage of passion fruit. When the samples deformed elastically to certain values, the hexagonal network lost stability and buckled into a chiral topography (Fig. 4a–d). Note that this mode transition is reversible when the air re-enters the sample and the pressure difference is restored. To further illustrate the tunability of the chiral localization, we applied a small disturbance (poke by a rod) somewhere on the surface to trigger the hexagonal-to-chiral mode transformation (Fig. 4e–h) while the sample was subjected to homogeneous shrinkage, which was in good agreement with finite element simulations (Fig. 4i–l). This strategy can provide enlightenment for the design of programmable functional surfaces such as adaptive grasping based on chiral localization.

Based on the aforementioned experiment, we present a target-adaptive gripper which can grasp small objects based on a hexagonal-to-chiral mode transformation. Simple structure, easy control, shape adaptation and filterable grasping are prominent advantages of the chiral gripper. The gripper system consists of a hemispherical shell with hexagonal topography, an air channel and a lifting frame that can move up and down (Supplementary Fig. 3). The air channel and the hemispherical part constitute a cavity structure, the former being connected to an external exhaust device to trigger the hexagonal-to-chiral mode transition by air extraction. The lifting frame is combined with the cap to control the motion. The working principle of the gripper is introduced as follows: The lifting frame descends to make the gripper approach a target. When the hexagonal network on the curved surface touches the object, the contact perturbation triggers the hexagonal-to-chiral topographic deformation that can well fit with the targeted shape. Then, the exhaust device begins to pump air. With increasing air extraction, the chiral topography can lock the object tightly to achieve a stable grasp. Finally, the object leaves the desk when raising the lifting frame. When the pressure difference is restored, the chiral topography elastically reverts back to hexagonal networks, releasing the grasped object. We carried out topographic grasping experiments on stiff or soft objects of different shapes and sizes (Fig. 5 and Supplementary Fig. 4). Our experiments showed that the gripper can smartly and stably grasp various small-sized objects. To further demonstrate the crucial role played by the chiral topography in robust grasping, we performed contrast experiments by making a hemispherical cap with a smooth surface. Except for the lack of the initial hexagonal network on the surface, the other parameters of the gripper remained exactly the same as in the aforementioned grasping experiments. With the smooth surface, the targets slid off, leading to failure of effective grasping (Supplementary Video 5). Our experiments not only prove the critical role of the chiral topography in effective, target-adaptive grasping but also shed light on smart gripper designs.

Further information on research design is available in the Nature Research Reporting Summary linked to this article.

Source data for the FEM computations shown in Figs. 2 and 3 are available with this manuscript.

The code used in this study can be obtained from Zenodo29.

Yin, J., Cao, Z., Li, C., Sheinman, I. & Chen, X. Stress-driven buckling patterns in spheroidal core/shell structures. Proc. Natl Acad. Sci. USA 105, 19132–19135 (2008).

Article Google Scholar

Yin, J., Chen, X. & Sheinman, I. Anisotropic buckling patterns in spheroidal film/substrate systems and their implications in some natural and biological systems. J. Mech. Phys. Solids 57, 1470–1484 (2009).

Article Google Scholar

Li, B., Jia, F., Cao, Y. P., Feng, X. Q. & Gao, H. Surface wrinkling patterns on a core–shell soft sphere. Phys. Rev. Lett. 106, 234301 (2011).

Article Google Scholar

Dervaux, J. & Ben Amar, M. Morphogenesis of growing soft tissues. Phys. Rev. Lett. 101, 068101 (2008).

Article Google Scholar

Huang, C., Wang, Z., Quinn, D., Suresh, S. & Hsia, K. J. Differential growth and shape formation in plant organs. Proc. Natl Acad. Sci. USA 115, 12359–12364 (2018).

Article Google Scholar

Xu, F., Fu, C. & Yang, Y. Water affects morphogenesis of growing aquatic plant leaves. Phys. Rev. Lett. 124, 038003 (2020).

Article Google Scholar

Ben Amar, M. & Jia, F. Anisotropic growth shapes intestinal tissues during embryogenesis. Proc. Natl Acad. Sci. USA 110, 10525–10530 (2013).

Article Google Scholar

Ciarletta, P., Balbi, V. & Kuhl, E. Pattern selection in growing tubular tissues. Phys. Rev. Lett. 113, 248101 (2014).

Article Google Scholar

Ciarletta, P. Buckling instability in growing tumor spheroids. Phys. Rev. Lett. 110, 158102 (2013).

Article Google Scholar

Tallinen, T. et al. On the growth and form of cortical convolutions. Nat. Phys. 12, 588–593 (2016).

Article Google Scholar

Brau, F. et al. Multiple-length-scale elastic instability mimics parametric resonance of nonlinear oscillators. Nat. Phys. 7, 56–60 (2011).

Article Google Scholar

Stoop, N., Lagrange, R., Terwagne, D., Reis, P. M. & Dunkel, J. Curvature-induced symmetry breaking determines elastic surface patterns. Nat. Mater. 14, 337–342 (2015).

Article Google Scholar

Katifori, E., Alben, S., Cerda, E., Nelson, D. R. & Dumais, J. Foldable structures and the natural design of pollen grains. Proc. Natl Acad. Sci. USA 107, 7635–7639 (2010).

Article Google Scholar

Holland, M., Budday, S., Goriely, A. & Kuhl, E. Symmetry breaking in wrinkling patterns: Gyri are universally thicker than sulci. Phys. Rev. Lett. 121, 228002 (2018).

Article Google Scholar

Bowden, N., Brittain, S., Evans, A. G., Hutchinson, J. W. & Whitesides, G. M. Spontaneous formation of ordered structures in thin films of metals supported on an elastomeric polymer. Nature 393, 146–149 (1998).

Article Google Scholar

Rogers, J. A., Someya, T. & Huang, Y. Materials and mechanics for stretchable electronics. Science 327, 1603–1607 (2010).

Article Google Scholar

Pocivavsek, L. et al. Topography-driven surface renewal. Nat. Phys. 14, 948–953 (2018).

Article Google Scholar

Pikul, J. H. et al. Stretchable surfaces with programmable 3D texture morphing for synthetic camouflaging skins. Science 358, 210–214 (2017).

Article Google Scholar

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

Article Google Scholar

Terwagne, D., Brojan, M. & Reis, P. M. Smart morphable surfaces for aerodynamic drag control. Adv. Mater. 26, 6608–6611 (2014).

Article Google Scholar

Cao, G., Chen, X., Li, C., Ji, A. & Cao, Z. Self-assembled triangular and labyrinth buckling patterns of thin films on spherical substrates. Phys. Rev. Lett. 100, 036102 (2008).

Article Google Scholar

Breid, D. & Crosby, A. J. Curvature-controlled wrinkle morphologies. Soft Matter 9, 3624–3630 (2013).

Google Scholar

Xu, F., Zhao, S., Lu, C. & Potier-Ferry, M. Pattern selection in core–shell spheres. J. Mech. Phys. Solids 137, 103892 (2020).

Article MathSciNet Google Scholar

van der Heijden, A. M. A. W.T. Koiter's Elastic Stability of Solids and Structures (Cambridge Univ. Press, 2009).

Biot, M. A. Bending of an infinite beam on an elastic foundation. J. Appl. Mech. 4, A1–A7 (1937).

Article Google Scholar

Allen, H. G. Analysis and Design of Structural Sandwich Panels (Pergamon, 1969).

Zhao, Y., Cao, Y., Feng, X. Q. & Ma, K. Axial compression-induced wrinkles on a core-shell soft cylinder: Theoretical analysis, simulations and experiments. J. Mech. Phys. Solids 73, 212–227 (2014).

Article MathSciNet Google Scholar

Groh, R. M. J., Avitabile, D. & Pirrera, A. Generalised path-following for well-behaved nonlinear structures. Comput. Methods Appl. Mech. Eng. 331, 394–426 (2018).

Article MathSciNet Google Scholar

Xu, F., Huang, Y., Zhao, S. & Feng, X. Q. Chiral topographic instability in shrinking spheres. Zenodo https://doi.org/10.5281/zenodo.7025830 (2022).

Download references

This work is supported by the National Natural Science Foundation of China (grants no. 12122204, 11872150 and 11921002), Shanghai Pilot Program for Basic Research-Fudan University (grant no. 21TQ1400100-21TQ010), Shanghai Shuguang Program (grant no. 21SG05), Shanghai Rising-Star Program (grant no. 19QA1400500) and young scientist project of the MOE innovation platform.

Institute of Mechanics and Computational Engineering, Department of Aeronautics and Astronautics, Fudan University, Shanghai, P. R. China

Fan Xu, Yangchao Huang & Shichen Zhao

Institute of Biomechanics and Medical Engineering, AML, Department of Engineering Mechanics, Tsinghua University, Beijing, P. R. China

Xi-Qiao Feng

You can also search for this author in PubMed Google Scholar

F.X. and X.-Q.F. conceived the idea. F.X. designed the research. Y.H. and S.Z. conducted the experiments. F.X. and Y.H. developed the theoretical models and carried out the dimensional analyses. Y.H. and S.Z. performed the numerical simulations. F.X., Y.H. and S.Z. interpreted the results. F.X. and Y.H. wrote the manuscript. All the authors provided helpful discussions.

Correspondence to Fan Xu or Xi-Qiao Feng.

The authors declare that they have no competing interests.

Nature Computational Science thanks Francesco Dal Corso, Ahmer Wadee and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Handling editor: Jie Pan, in collaboration with the Nature Computational Science team. Peer reviewer reports are available.

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

(a) Cross section. (b) Geometry of a core-shell sphere. (c) Schematic of chiral buckling of a Y-shaped cellular representative layered plate.

(a) Pour liquid silicone on a 3D printed mold with hexagonal network on the surface. (b) Create a cavity in the sample by using a hemispherical cover when the liquid silicone is curing. (c) A silicon hemispherical shell with hexagonal pattern on the surface. (d) Two hemispherical shells with hexagonal network are glued and can be separated, with a channel connecting the internal cavity and external tube for air extraction.

Supplementary Figs. 1–4 and Table 1.

Natural dehydration of passion fruit and numerical simulation.

Hexagonal-to-chiral topography transition induced by air extraction.

Chiral topography formation induced by surface disturbance.

Chiral topography for adaptive grasping.

Contrast experiments with smooth surface

FEM source data

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Xu, F., Huang, Y., Zhao, S. et al. Chiral topographic instability in shrinking spheres. Nat Comput Sci 2, 632–640 (2022). https://doi.org/10.1038/s43588-022-00332-y

Download citation

Received: 23 April 2022

Accepted: 09 September 2022

Published: 24 October 2022

Issue Date: October 2022

DOI: https://doi.org/10.1038/s43588-022-00332-y

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

Nature Computational Science (2022)