# Growth patterns for shape-shifting elastic bilayers

^{a}John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA 02138;^{b}Institute for Computational Engineering and Sciences, University of Texas at Austin, Austin, TX 78712;^{c}Wyss Institute for Biologically Inspired Engineering, Harvard University, Cambridge, MA 02138;^{d}Department of Physics, Harvard University, Cambridge, MA 02138;^{e}Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA 02138;^{f}Kavli Institute for NanoBio Science and Technology, Harvard University, Cambridge, MA 02138

See allHide authors and affiliations

Edited by John A. Rogers, Northwestern University, Evanston, IL, and approved September 6, 2017 (received for review May 31, 2017)

## Significance

Many biological forms, such as leaves, flowers, and faces, are shaped by complex growth patterns. How can we prescribe the rules of growth on a simple surface so that it will morph into a flower or a face? Here, we solve this inverse problem of designing the growth patterns for an anisotropically growing elastic bilayer structure and prove that it can be used to achieve any target surface shape from any reference shape. We demonstrate the applicability of this result via the computational design of growth patterns for animal, vegetable, and mineral surfaces—a face, a flower, and a canyon. Our solution provides algorithms for engineering complex functional shapes in tissues, and actuation systems in soft robotics, and elsewhere.

## Abstract

Inspired by the differential-growth-driven morphogenesis of leaves, flowers, and other tissues, there is increasing interest in artificial analogs of these shape-shifting thin sheets made of active materials that respond to environmental stimuli such as heat, light, and humidity. But how can we determine the growth patterns to achieve a given shape from another shape? We solve this geometric inverse problem of determining the growth factors and directions (the metric tensors) for a given isotropic elastic bilayer to grow into a target shape by posing and solving an elastic energy minimization problem. A mathematical equivalence between bilayers and curved monolayers simplifies the inverse problem considerably by providing algebraic expressions for the growth metric tensors in terms of those of the final shape. This approach also allows us to prove that we can grow any target surface from any reference surface using orthotropically growing bilayers. We demonstrate this by numerically simulating the growth of a flat sheet into a face, a cylindrical sheet into a flower, and a flat sheet into a complex canyon-like structure.

Nonuniform in-plane growth of thin sheets generically leads to metric frustration that is relieved by out-of-plane buckling. This mechanism lies at the heart of many morphogenetic processes in botany, such as the shaping of a leaf (1), the blooming of a flower (2), or the explosive dispersal of seeds from certain pods (3). From an engineering perspective, these examples raise the possibility of biomimetic design: programming shape (morphogramming) into matter that can be actuated with environmental signals such as light, temperature, or concentration.

The theory of non-Euclidean plates and shells that links the elastic response of materials into the correctly invariant framework provided by differential geometry (1, 4⇓⇓⇓–8) is a natural starting place to analyze the growth and form of sheets and shells. Growing a thin structure by changing the in-plane intrinsic distances and angles between material elements makes its metric non-Euclidean; generically, this implies that the strain-free reference configuration may not be physically realizable in 3D space. Therefore, the system settles into a residually strained equilibrium configuration that is determined by a local minimum of the energetic cost of stretching and bending the sheet. Generally, this state might not be unique, and, typically, there will be a range of metastable configurations accessible to the system. This raises the natural question of the inverse problem: How should one program growth patterns into a sheet so that it morphs into a prescribed target shape? Recent attempts to solve this question have focused on theoretical designs of optimal growth patterns for the weakly nonlinear deformations of thin shape-shifting isotropic elastic sheets (9), axisymmetric growth patterns for morphable shells (10), design of director fields into deformable nematic elastic sheets (11), or a 4D phytomimetic printing approach based on a linearized elastic analysis to derive the print paths of an anisotropic bilayer made of a responsive ink (12). However, there is no general theoretical or computational framework to solve this inverse problem.

Here, we address this question in the context of a growing elastic bilayer, inspired by the growth and form of plant organs such as leaves and flowers that are usually made of two cell layers that adhere to each other and can grow independently. This bilayer geometry may be naturally described in terms of an infinitesimally-thin “midsurface” and a “thickness” *i*) occurs only in the in-plane directions (tangent to the midsurface) and that (*ii*) the shell can be divided into several “layers,” with growth constant through the thickness direction for each layer. These assumptions allow us to represent growth within each layer as a tensor field specified at each point of the midsurface. “Isotropic growth” consists of an equal growth factor in all in-plane directions at every point, and is thus encoded by one independent degree of freedom (the isotropic scaling factor) per layer at each point on the midsurface. With “orthotropic growth,” the growth factor of the material at each point is a function of the in-plane direction, providing three degrees of freedom per layer for each location on the midsurface: the two growth factors in the orthogonal “principal growth directions” and the planar rotation angle of this axis. Given the ubiquity of the bilayer geometry in plant organs, and the ease of additive manufacturing techniques that allow us to approach this possibility, we will focus on this case from now on.

Such a bilayer where each of the two layers can experience independent orthotropic growth (Fig. 1, *Left*) has six degrees of freedom at every location along the common midsurface. Any choice of these degrees of freedom constitutes a growth field for the shell, giving rise to an equilibrium shape in

Classical differential geometry of 2D surfaces (13) tells us that a surface is uniquely defined by two symmetric quadratic forms (the first and second fundamental form), which consist of six quantities at every location along it. Naive counting of the degrees of freedom suggests that since we have three degrees of freedom associated with the in-plane growth of each of the layers in the composite bilayer, solutions to the inverse problem should be possible. In this work, we show that a solution does indeed always exist for the inverse problem when regularized as an elastic energy minimization problem. Furthermore, while this solution typically does not yield a resultant shape with zero residual strain (due to incompatibility when crossing from one layer to the other through the midsurface), we show that the residual strain is constant (i.e., independent of the realization of the bilayer). We provide a simple algebraic expression for the growth factors and angles that achieve this solution, allowing any bilayer to grow into any target shape.

## Geometry and Elasticity

We parameterize the midsurface of a shell using curvilinear coordinates *Right*. By using the Kirchhoff–Love assumption, any material point

Similarly, we can describe the “growth” of

### Elastic Energy of a Curved Monolayer.

Given the “rest quantities” *SI Appendix* for details) the depth-integrated elastic energy as*SI Appendix* for details of this equivalence).

### Elastic Energy of a Bilayer.

A bilayer is made of two monolayers of thickness *Left*). Here, ‘bottom’ and ‘top’ are used with the convention of midsurface normal vectors pointing upward (i.e., into the layer whose metric is

The elastic energy of a bilayer made of two layers with metrics given by *SI Appendix*)**1**.

### Energy Equivalence Between Monolayers and Bilayers.

A natural question that is raised by the geometric and mechanical description of the composite bilayer is whether we can relate it to an equivalent monolayer with appropriate first and second fundamental forms. In this work we show that, indeed, if

To see this equivalence, consider the ansatz *SI Appendix* for details)**3** and **4** become functions of *SI Appendix* for details).

## Forward Problem of Growth

The result of Eqs. **3** and **4** can be used to compute the quadratic forms *SI Appendix* for the mathematical details of these constraints.

To give a sense of how this plays out practically, we consider two examples of the forward problem. First, consider the case of Fig. 1, *Left*, where we prescribed orthotropic growth to each of the layers, with the top layer growing in one principal direction and the bottom layer growing in the orthogonal direction. Both layers experience 1D expansion with the same, constant factor **4**. It is easy to see that these two forms, being spatially homogeneous, automatically satisfy both of the Peterson–Mainardi–Codazzi equations. Gauss’ equation, however, cannot be satisfied: The Gauss curvature *Lower Left* is therefore characterized by residual strain (see *SI Appendix* for details).

Second, consider a modification of the example of Fig. 1, where the two layers instead grow in the same direction, but with different growth factors *SI Appendix* for more details). Our solution, summarized in Eqs. **3** and **4**, therefore not only generalizes Timoshenko’s theory to nonlinear elasticity and arbitrary growth factors, but can also be seen as its extension that allows us to proceed from strips to surfaces.

## Inverse Problem of Growth

As noted above, inhomogeneous orthotropic growth for each layer in the bilayer can be represented in terms of three degrees of freedom at every point—two growth factors corresponding to the growth in orthogonal principal directions and one angle to define the axes. Together, these three values define the symmetric positive definite first fundamental form at a given location. For orthotropic growth of a bilayer, we therefore have independent control of the first fundamental forms for the bottom and top layer, which provides six degrees of freedom for each point on the surface. This leads to the following theorem:

### Theorem 1.

*Given any surface* *and any shaped planar region* *(topologically compatible with* *) with embedding* *, so that the immersion of* *is equal to* *, there exists a bilayer* *, with sufficiently small desired thickness* *so that* *defines the energetically equilibrated configuration of the bilayer’s midsurface. This bilayer is defined by**where* *and* *are the first and second fundamental forms of* *Moreover, this embedding is free of residual strain, with the exception of a deformation-independent strain due to incompatibility of the bilayers at their common boundary. Finally, the metrics* *and* *as defined in* Eq. **5** can *be decomposed as orthotropic growth of a bilayer on*

#### Proof.

Eq. **3** states that the minimum energy embedding of a curved monolayer, with **1** we can immediately see that the minimum energy configuration of a curved monolayer with embedding **4**, resulting in Eq. **5**. From Eq. **3**, we can see that the resulting bilayer embedding is free of residual strain except for the already-mentioned deformation-independent term. Finally, by performing a spectral decomposition of the metric, we can write both

With this theorem, we have the capability to grow any initial bilayer structure into any target shape. However, we need to specify two practical considerations. First of all, the metrics **5** have to be positive-definite to be admissible (13). We show in the *SI Appendix* that this results in the constraint

Secondly, the trajectory of growth that is followed to transition from initial to final bilayer metrics might pose complications in practice. Although the embedding of the geometry that globally minimizes the energy given the final bilayer metrics is always unique and identical to the target configuration, some growth trajectories might result in a metastable shape different from the target shape. Since any metastable shape can always be “snapped” into the desired target shape, this issue is of secondary importance to our main contribution, yet provides some interesting questions. In particular, to guarantee that a growth process always results in the target shape, we anticipate two numerical and/or physical difficulties that could arise during interpolation of the metrics. First, it is important for the symmetric positive definite (SPD) matrices encoding the bilayer metrics to vary smoothly over time and remain SPD throughout the interpolation, and, moreover, for the eigenvectors and eigenvalues (encoding the anisotropy amount and direction) to vary as smoothly as possible during interpolation. Second, interpolated bilayer metrics at intermediate growth stages should remain as compatible as possible; by this, we mean that the first and second fundamental form corresponding to the interpolated bilayer metrics should deviate as little as possible from satisfying the Gauss and Peterson–Mainardi–Codazzi compatibility relations. Otherwise, residual strain accumulates during the interpolation and could pose an energetic barrier between the final realized shape and the target solution. Again, the target solution would still be the global energetic minimum, but physically or numerically, we could find ourselves stuck at a local minimum. A particular numerical example of this is further detailed in the *SI Appendix*, for the case of the snapdragon flower growth case shown below.

A systematic approach to completely avoiding metastable states during the growth process is a very interesting direction for future work, and we can now pose the problem in a crisp way: Is it possible to interpolate two pairs of bilayer metrics, so that the equivalent minimum-energy midsurface fundamental forms are always compatible? Even more interesting is the prospect to harness incompatibility at intermediate growth stages to control the final grown shape, or exploit multiple solutions depending on the spatiotemporal distribution of growth, for instance, by snap-through (18).

## Results

For parametrized surfaces, we can use Eq. **5** directly to solve the inverse problem algebraically. In simple situations associated with parametric surfaces that have explicit fundamental forms (e.g., a hemisphere, catenoid, and a saddle), we can carry out these computations analytically (*SI Appendix*). However, to truly demonstrate the usefulness of our approach, we need to show how to design growth patterns for complex shapes with multiple spatial scales. We do this by using numerical methods to minimize Eq. **2**. For a given 3D target surface, we first triangulate it using a mesh that can capture the smallest length scale of interest, and then compute the current first and second fundamental forms **5**, and decompose each of them spectrally into local orthotropic growth rules for a given reference mesh and its corresponding bilayer metrics. We then interpolate the growth factors on each point in the initial mesh into a set of discrete values between those of the reference and target surfaces and solve a sequence of problems to determine the intermediate equilibrium configuration at each of the discrete growth steps. This provides a way to visualize the transition from initial to final configuration in a quasi-static manner, and further allows us to design arbitrary way-points between the initial and final state (*SI Appendix*).

Inspired by recent work on floral morphogenesis (19⇓–21), we first show how we can grow a cylinder made of two thin sheets of the same uniform thickness into a snapdragon flower. In this case, we assume the initial cylinder, with **5**. The snapdragon flower, with *SI Appendix* for more details). In Fig. 2, *Left*, and Movie S1, we show a sequence of intermediate shapes viewed from two different angles and, for comparison, the actual snapdragon flower.

To highlight the ability of our approach to capture complex surface geometries with features on multiple scales, we turn to the human face. In Fig. 2, *Center*, and Movie S2, we show that we can grow an initially flat bilayer disk into a 3D model of a human face, that of the physicist Max Planck. In this case, we linearly interpolate the growth factors from unity, on the initial disk, to their final values computed from Eq. **5**, and show the result for energetic equilibria at two intermediate stages.

Finally, we use our inverse-design theory to grow a simulacrum of a complex inanimate surface, a horseshoe bend in the Colorado River in Arizona. Using United States Geological Survey (USGS) elevation data, we create our 3D target shape and grow it from a rectangular sheet. As in the snapdragon example, we use anchor points to guide the interpolation (see *SI Appendix* for more details). In Fig. 2, *Right*, and Movie S3, we show a sequence of intermediate shapes obtained during the growth process. These numerical results demonstrate the practical validity of our theoretical framework that allows us to capture the shapes of complex absolute-scale-independent surfaces from the animal, vegetable, and mineral world.

## Discussion

This study poses and solves the inverse-design problem of designing growth patterns for creating complex shapes from a uniformly thin isotropic elastic bilayer capable of sustaining orthotropic growth. It opens the way for formulating and solving other variants of the inverse-design problem for growth-metric tensors encountered in such cases as a single growing monolayer, isotropically growing bilayers (23), or the most general case of orthotropically growing bilayers with incompatible metrics, all of which will generally require numerical approaches for both the forward and inverse problems. The generalization to account for situations where the thicknesses of the two layers *SI Appendix*) provides yet another perspective, as this can be exploited for applications such as artificial lenses (24) or controlled actuation of the curvature in the presence of constraints. It is worth noting that in all of these situations, there are specific instances where the final residually strained state may constitute an orbit of connected minima associated with a Goldstone mode, e.g., saddle-like solutions in heated or swollen plates (refs. 9 and 15, p. 158), but the general nature of these states and how to design them remains open.

Our approach is agnostic to the actual mechanism that induces this growth, whether it is heat, light, swelling, or biological growth, as long as it is mathematically equivalent to changing the metric of the constituent layers. While the experimental realization of bilayers is easy using additive manufacturing, achieving general orthotropic growth is currently more challenging. An existing technique for orthotropic growth relies on locally embedding fibrils in an isotropically growing base material, providing control over the angle between two fixed orthotropic swelling factors through the print direction (12). To also vary the other two degrees of freedom would require further control over either the density of the fibrils in two orthogonal directions, or the density in one direction together with the isotropic growth factor of the base material, for each point in each layer. Alternatively, one could consider discrete lattice or origami-type surfaces that approximate the features of a continuous surface at a larger scale. This simplifies the problem by replacing the exact local control of growth by approximate control of an appropriate nonlocal average and is likely to be the first to be realizable experimentally.

## Acknowledgments

We thank E. Coen for the digital surface meshes and image of the snapdragon flower in Fig. 2. This work was supported in part by the Swiss National Science Foundation (to W.M.v.R.), National Science Foundation Grant DMS-1304211 (to E.V.), National Science Foundation Grant DMR 14-20570 (to L.M.), DMREF 15-33985 (to L.M.), and Army Research Office Grant W911NF-15-1-0166 (to L.M.).

## Footnotes

↵

^{1}Present address: Department of Mechanical Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139.- ↵
^{2}To whom correspondence should be addressed. Email: lmahadev{at}g.harvard.edu.

Author contributions: L.M. conceived the study; W.M.v.R., E.V., and L.M. designed research; W.M.v.R. and E.V. performed research; W.M.v.R., E.V., and L.M. analyzed data; and W.M.v.R., E.V., and L.M. wrote the paper.

The authors declare no conflict of interest.

This article is a PNAS Direct Submission.

This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1709025114/-/DCSupplemental.

## References

- ↵.
- Liang H,
- Mahadevan L

- ↵.
- Liang H,
- Mahadevan L

- ↵.
- Armon S,
- Efrati E,
- Kupferman R,
- Sharon E

- ↵
- ↵.
- Klein Y,
- Efrati E,
- Sharon E

- ↵
- ↵
- ↵.
- Lewicka M,
- Mahadevan L,
- Pakzad MR

- ↵.
- Jones GW,
- Mahadevan L

- ↵.
- Kim J,
- Hanna JA,
- Byun M,
- Santangelo CD,
- Hayward RC

- ↵
- ↵
- ↵.
- Struik DJ

- ↵.
- Weischedel C,
- Tuganov A,
- Hermansson T,
- Linn J,
- Wardetzky M

- ↵.
- Mansfield E

- ↵.
- do Carmo MP

- ↵
- ↵
- ↵
- ↵
- ↵.
- Coen E,
- Rebocho AB

- ↵
- ↵.
- Pezzulla M,
- Smith GP,
- Nardinocchi P,
- Holmes DP

- ↵.
- Lee LP

## Citation Manager Formats

## Article Classifications

- Physical Sciences
- Applied Physical Sciences