Jesse L. Silverberg

Using origami design principles to fold reprogrammable mechanical metamaterials

Folding robots and metamaterials The same principles used to make origami art can make self-assembling robots and tunable metamaterials—artificial materials engineered to have properties that may not be found in nature (see the Perspective by You). Felton et al. made complex self-folding robots from flat templates. Such robots could potentially be sent through a collapsed building or tunnels and then assemble themselves autonomously into their final functional form. Silverberg et al. created a mechanical metamaterial that was folded into a tessellated pattern of unit cells. These cells reversibly switched between soft and stiff states, causing large, controllable changes to the way the material responded to being squashed. Science, this issue p. 644, p. 647; see also p. 623 Origami folded sheets can be structurally altered by adding defects to control the mechanical properties. [Also see Perspective by You] Although broadly admired for its aesthetic qualities, the art of origami is now being recognized also as a framework for mechanical metamaterial design. Working with the Miura-ori tessellation, we find that each unit cell of this crease pattern is mechanically bistable, and by switching between states, the compressive modulus of the overall structure can be rationally and reversibly tuned. By virtue of their interactions, these mechanically stable lattice defects also lead to emergent crystallographic structures such as vacancies, dislocations, and grain boundaries. Each of these structures comes from an arrangement of reversible folds, highlighting a connection between mechanical metamaterials and programmable matter. Given origami’s scale-free geometric character, this framework for metamaterial design can be directly transferred to milli-, micro-, and nanometer-size systems.

Origami structures with a critical transition to bistability arising from hidden degrees of freedom

Origami is used beyond purely aesthetic pursuits to design responsive and customizable mechanical metamaterials. However, a generalized physical understanding of origami remains elusive, owing to the challenge of determining whether local kinematic constraints are globally compatible and to an incomplete understanding of how the folded sheets material properties contribute to the overall mechanical response. Here, we show that the traditional square twist, whose crease pattern has zero degrees of freedom (DOF) and therefore should not be foldable, can nevertheless be folded by accessing bending deformations that are not explicit in the crease pattern. These hidden bending DOF are separated from the crease DOF by an energy gap that gives rise to a geometrically driven critical bifurcation between mono- and bistability. Noting its potential utility for fabricating mechanical switches, we use a temperature-responsive polymer-gel version of the square twist to demonstrate hysteretic folding dynamics at the sub-millimetre scale.

Collective Motion of Humans in Mosh and Circle Pits at Heavy Metal Concerts

Human collective behavior can vary from calm to panicked depending on social context. Using videos publicly available online, we study the highly energized collective motion of attendees at heavy metal concerts. We find these extreme social gatherings generate similarly extreme behaviors: a disordered gas-like state called a mosh pit and an ordered vortex-like state called a circle pit. Both phenomena are reproduced in flocking simulations demonstrating that human collective behavior is consistent with the predictions of simplified models.

3D imaging and mechanical modeling of helical buckling in Medicago truncatula plant roots

We study the primary root growth of wild-type Medicago truncatula plants in heterogeneous environments using 3D time-lapse imaging. The growth medium is a transparent hydrogel consisting of a stiff lower layer and a compliant upper layer. We find that the roots deform into a helical shape just above the gel layer interface before penetrating into the lower layer. This geometry is interpreted as a combination of growth-induced mechanical buckling modulated by the growth medium and a simultaneous twisting near the root tip. We study the helical morphology as the modulus of the upper gel layer is varied and demonstrate that the size of the deformation varies with gel stiffness as expected by a mathematical model based on the theory of buckled rods. Moreover, we show that plant-to-plant variations can be accounted for by biomechanically plausible values of the model parameters.

Structure-Function Relations and Rigidity Percolation in the Shear Properties of Articular Cartilage

Among mammalian soft tissues, articular cartilage is particularly interesting because it can endure a lifetime of daily mechanical loading despite having minimal regenerative capacity. This remarkable resilience may be due to the depth-dependent mechanical properties, which have been shown to localize strain and energy dissipation. This paradigm proposes that these properties arise from the depth-dependent collagen fiber orientation. Nevertheless, this structure-function relationship has not yet been quantified. Here, we use confocal elastography, quantitative polarized light microscopy, and Fourier-transform infrared imaging to make same-sample measurements of the depth-dependent shear modulus, collagen fiber organization, and extracellular matrix concentration in neonatal bovine articular cartilage. We find weak correlations between the shear modulus |G(∗)| and both the collagen fiber orientation and polarization. We find a much stronger correlation between |G(∗)| and the concentration of collagen fibers. Interestingly, very small changes in collagen volume fraction vc lead to orders-of-magnitude changes in the modulus with |G(∗)| scaling as (vc - v0)(ξ). Such dependencies are observed in the rheology of other biopolymer networks whose structure exhibits rigidity percolation phase transitions. Along these lines, we propose that the collagen network in articular cartilage is near a percolation threshold that gives rise to these large mechanical variations and localization of strain at the tissues surface.

Lattice mechanics of origami tessellations.

Origami-based design holds promise for developing materials whose mechanical properties are tuned by crease patterns introduced to thin sheets. Although there have been heuristic developments in constructing patterns with desirable qualities, the bridge between origami and physics has yet to be fully developed. To truly consider origami structures as a class of materials, methods akin to solid mechanics need to be developed to understand their long-wavelength behavior. We introduce here a lattice theory for examining the mechanics of origami tessellations in terms of the topology of their crease pattern and the relationship between the folds at each vertex. This formulation provides a general method for associating mechanical properties with periodic folded structures and allows for a concrete connection between more conventional materials and the mechanical metamaterials constructed using origami-based design.

Effects of enzymatic treatments on the depth-dependent viscoelastic shear properties of articular cartilage

Osteoarthritis (OA) is a disease that involves the erosion and structural weakening of articular cartilage. OA is characterized by the degradation of collagen and proteoglycans in the extracellular matrix (ECM), particularly at the articular surface by proteinases including matrix metalloproteinases (MMPs) and a disintegrin and metalloproteinase with thrombospondin motifs (ADAMTSs).1 Degradation of collagen and proteoglycans is known to alter shear mechanical properties of cartilage, but study of this phenomenon has been focused on bulk tissue properties. The purpose of this study was to assess microscale cartilage damage induced by trypsin or collagenase using a technique to measure the local shear viscoelastic properties. Safranin‐O histology revealed a decrease in proteoglycans near the articular surface after collagenase and trypsin digestions, with proteoglycan depletion increasing in time. Similarly, confocal reflectance micrographs showed increasing collagen degradation in collagenase treated samples, although the collagen network remained intact after trypsin treatment. Both treatments induced changes in shear modulus that were confined to a narrow range (∼400µm) near tissue surface. In addition, collagenase altered the total energy dissipation distribution by up to a factor of 100, with longer digestion times corresponding to higher energy dissipation. The ability to detect local mechanical signatures in tissue composition and mechanics is an important tool for understanding the spatially non‐uniform changes that occur in articular cartilage diseases such as OA.

Anatomic variation of depth‐dependent mechanical properties in neonatal bovine articular cartilage

Articular cartilage has well known depth‐dependent structure and has recently been shown to have similarly non‐uniform depth‐dependent mechanical properties. Here, we study anatomic variation of the depth‐dependent shear modulus and energy dissipation rate in neonatal bovine knees. The regions we specifically focus on are the patellofemoral groove, trochlea, femoral condyle, and tibial plateau. In every sample, we find a highly compliant region within the first 500 µm of tissue measured from the articular surface, where the local shear modulus is reduced by up to two orders of magnitude. Comparing measurements taken from different anatomic sites, we find statistically significant differences localized within the first 50 µm. Histological images reveal these anatomic variations are associated with differences in collagen density and fiber organization.

X-ray computed tomography uncovers root–root interactions: quantifying spatial relationships between interacting root systems in three dimensions

Research in the field of plant biology has recently demonstrated that inter- and intra-specific interactions belowground can dramatically alter root growth. Our aim was to answer questions related to the effect of inter- vs. intra-specific interactions on the growth and utilization of undisturbed space by fine roots within three dimensions (3D) using micro X-ray computed tomography. To achieve this, Populus tremuloides (quaking aspen) and Picea mariana (black spruce) seedlings were planted into containers as either solitary individuals, or inter-/intra-specific pairs, allowed to grow for 2 months, and 3D metrics developed in order to quantify their use of belowground space. In both aspen and spruce, inter-specific root interactions produced a shift in the vertical distribution of the root system volume, and deepened the average position of root tips when compared to intra-specifically growing seedlings. Inter-specific interactions also increased the minimum distance between root tips belonging to the same root system. There was no effect of belowground interactions on the radial distribution of roots, or the directionality of lateral root growth for either species. In conclusion, we found that significant differences were observed more often when comparing controls (solitary individuals) and paired seedlings (inter- or intra-specific), than when comparing inter- and intra-specifically growing seedlings. This would indicate that competition between neighboring seedlings was more responsible for shifting fine root growth in both species than was neighbor identity. However, significant inter- vs. intra-specific differences were observed, which further emphasizes the importance of biological interactions in competition studies.

Decoupling local mechanics from large-scale structure in modular metamaterials

Significance The forward-design approach to mechanical metamaterials determines a structure’s properties after it is designed. In contrast, the inverse-design approach specifies parameter domains then optimizes a cost function to achieve a desired property. The former approach guarantees a structure will exist but lacks the ability to prescribe function. The latter approach has prescribed function, but existence is not guaranteed. Here, we trivialize these design challenges for metamaterials by decoupling local mechanical properties from the bulk structure, allowing both to be specified independently. This work introduces a design strategy that substantially advances the capacity to engineer mechanical metamaterials by specifically using modular units with more free design parameters than constraints. The fundamental strategy is explicitly demonstrated with an origami- and kirigami-inspired structure. A defining feature of mechanical metamaterials is that their properties are determined by the organization of internal structure instead of the raw fabrication materials. This shift of attention to engineering internal degrees of freedom has coaxed relatively simple materials into exhibiting a wide range of remarkable mechanical properties. For practical applications to be realized, however, this nascent understanding of metamaterial design must be translated into a capacity for engineering large-scale structures with prescribed mechanical functionality. Thus, the challenge is to systematically map desired functionality of large-scale structures backward into a design scheme while using finite parameter domains. Such “inverse design” is often complicated by the deep coupling between large-scale structure and local mechanical function, which limits the available design space. Here, we introduce a design strategy for constructing 1D, 2D, and 3D mechanical metamaterials inspired by modular origami and kirigami. Our approach is to assemble a number of modules into a voxelized large-scale structure, where the module’s design has a greater number of mechanical design parameters than the number of constraints imposed by bulk assembly. This inequality allows each voxel in the bulk structure to be uniquely assigned mechanical properties independent from its ability to connect and deform with its neighbors. In studying specific examples of large-scale metamaterial structures we show that a decoupling of global structure from local mechanical function allows for a variety of mechanically and topologically complex designs.