Modeling the formation of double rolls from heterogeneously patterned gels

Biomimetic growth in polymer gels

By modeling gels growing in confined environments, we uncover a biomimetic feedback mechanism between the evolving gel and confining walls that enables significant control over the properties of the grown gel. Our new model describes the monomer adsorption, polymerization and cross-linking involved in forming new networks and the resultant morphology and mechanical behavior of the grown gel. Confined between two hard walls, a thin, flat "parent" gel undergoes buckling; removal of the walls returns the gel to the flat structure. Polymerization and cross-linking in the confined parent generates the next stage of growth, forming a random copolymer network (RCN). When the walls are removed, the RCN remains in the buckled state, simultaneously "locking in" these patterns and increasing the Young's modulus by two orders of magnitude. Confinement of thicker gels between harder or softer 3D walls leads to controllable mechanical heterogeneities, where the Young's modulus between specific domains can differ by three orders of magnitude. These systems effectively replicate the feedback between mechanics and morphology in biological growth, where mechanical forces guide the structure formation throughout stages of growth. The findings provide new guidelines for shaping "growing materials" and introducing new approaches to matching form and function in synthetic systems.

Pure measures of bending for soft plates

This paper, originally motivated by a question raised by Wood and Hanna [Soft Matter, 2019, 15, 2411], shows that pure measures of bending for soft plates can be defined by introducing the class of bending-neutral deformations, which represent finite incremental changes in the plate's shape that do not induce any additional bending. This class of deformations is subject to a geometric compatibility condition, which is fully characterized. A tensorial pure measure of bending, which is invariant under bending-neutral deformations, is described in detail. As shown by an illustrative class of examples, the general notion of a pure measure of bending could be useful in formulating direct theories for soft plates, where stretching and bending energies are treated separately.

Formation of rolls from liquid crystal elastomer bistrips

Formation of desired three-dimensional (3D) shapes from flat thin sheets with programmed non-uniform deformation profiles is an effective strategy to create functional 3D structures. Liquid crystal elastomers (LCEs) are of particular use in programmable shape morphing due to their ability to undergo large, reversible, and anisotropic deformation in response to a stimulus. Here we consider a rectangular monodomain LCE thin sheet divided into one high- and one low-temperature strip, which we dub a 'bistrip'. Upon activation, a discontinuously patterned, anisotropic in-plane stretch profile is generated, and induces buckling of the bistrip into a rolled shape with a transitional bottle neck. Based on the non-Euclidean plate theory, we derive an analytical model to quantitatively capture the formation of the rolled shapes from a flat bistrip with finite thickness by minimizing the total elastic energy involving both stretching and bending energies. Using this analytical model, we identify the critical thickness at which the transition from the unbuckled to buckled configuration occurs. We further study the influence of the anisotropy of the stretch profile on the rolled shapes by first converting prescribed metric tensors with different anisotropy to a unified metric tensor embedded in a bistrip of modified geometry, and then investigating the effect of each parameter in this unified metric tensor on the rolled shapes. Our analysis sheds light on designing shape morphing of LCE thin sheets, and provides quantitative predictions on the 3D shapes that programmed LCE sheets can form upon activation for various applications.

Volume-constrained deformation of a thin sheet as a route to harvest elastic energy

Thin sheets exhibit rich morphological structures when subjected to external constraints. These structures store elastic energy that can be released on demand when one of the constraints is suddenly removed. Therefore, when adequately controlled, shape changes in thin bodies can be utilized to harvest elastic energy. In this paper, we propose a mechanical setup that converts the deformation of the thin body into a hydrodynamic pressure that potentially can induce a flow. We consider a closed chamber that is filled with an incompressible fluid and is partitioned symmetrically by a long and thin sheet. Then, we allow the fluid to exchange freely between the two parts of the chamber, such that its total volume is conserved. We characterize the slow, quasistatic, evolution of the sheet under this exchange of fluid, and derive an analytical model that predicts the subsequent pressure drop in the chamber. We show that this evolution is governed by two different branches of solutions. In the limit of a small lateral confinement we obtain approximated solutions for the two branches and characterize the transition between them. Notably, the transition occurs when the pressure drop in the chamber is maximized. Furthermore, we solve our model numerically and show that this maximum pressure behaves nonmonotonically as a function of the lateral compression.

Buckling-induced interaction between circular inclusions in an infinite thin plate

Design of slender artificial materials and morphogenesis of thin biological tissues typically involve stimulation of isolated regions (inclusions) in the growing body. These inclusions apply internal stresses on their surrounding areas that are ultimately relaxed by out-of-plane deformation (buckling). We utilize the Föppl-von Kármán model to analyze the interaction between two circular inclusions in an infinite plate that their centers are separated a distance of 2ℓ. In particular, we investigate a region in phase space where buckling occurs at a narrow transition layer of length ℓ_{D} around the radius of the inclusion, R (ℓ_{D}≪R). We show that the latter length scale defines two regions within the system, the close separation region, ℓ-R∼ℓ_{D}, where the transition layers of the two inclusions approximately coalesce, and the far separation region, ℓ-R≫ℓ_{D}. While the interaction energy decays exponentially in the latter region, E_{int}∝e^{-(ℓ-R)/ℓ_{D}}, it presents nonmonotonic behavior in the former region. While this exponential decay is predicted by our analytical analysis and agrees with the numerical observations, the close separation region is treated only numerically. In particular, we utilize the numerical investigation to explore two different scenarios within the final configuration: The first where the two inclusions buckle in the same direction (up-up solution) and the second where the two inclusions buckle in opposite directions (up-down solution). We show that the up-down solution is always energetically favorable over the up-up solution. In addition, we point to a curious symmetry breaking within the up-down scenario; we show that this solution becomes asymmetric in the close separation region.

Modeling the behavior of inclusions in circular plates undergoing shape changes from two to three dimensions

Growth of biological tissues and shape changes of thin synthetic sheets are commonly induced by stimulation of isolated regions (inclusions) in the system. These inclusions apply internal forces on their surroundings that, in turn, promote 2D layers to acquire complex 3D configurations. We focus on a fundamental building block of these systems, and consider a circular plate that contains an inclusion with dilative strains. Based on the Föppl-von Kármán (FvK) theory, we derive an analytical model that predicts the 2D-to-3D shape transitions in the system. Our findings are summarized in a phase diagram that reveals two distinct configurations in the post-buckling region. One is an extensive profile that holds close to the threshold of the instability, and the second is a localized profile, which preempts the extensive solution beyond the buckling threshold. While the former solution is derived as a perturbation around the flat configuration, assuming infinitesimal amplitudes, the latter solution is derived around a buckled state that is highly localized. We show that up to vanishingly small corrections that scale with the thickness, this localized configuration is equivalent to that expected for ultra-thin sheets, which completely relax compressive stresses. Our findings agree quantitatively with direct numerical minimization of the FvK energy. Furthermore, we extend the theory to describe shape transitions in polymeric gels, and compare the results with numerical simulations that account for the complete elastodynamic behavior of the gels. The agreement between the theory and these simulations indicates that our results are observable experimentally. Notably, our findings can provide guidelines to the analysis of more complicated systems that encompass interaction between several buckled inclusions.