Modeling of surface mechanical behaviors of soft elastic solids: theory and examples

From a distance: Shuttleworth revisited

The Shuttleworth equation: a linear stress-strain relation ubiquitously used in modeling the behavior of soft surfaces. Its validity in the realm of materials subject to large deformation is a topic of current debate. Here, we allow for large deformation by deriving the constitutive behavior of the surface from the general framework of finite kinematics. We distinguish cases of finite and infinitesimal surface relaxation preceding an infinitesimal applied deformation. The Shuttleworth equation identifies as the Cauchy stress measure in the fully linearized setting. We show that both in finite and linearized cases, measured elastic constants depend on the utilized stress measure. In addition, we discuss the physical implications of our results and analyze the impact of surface relaxation on the estimation of surface elastic moduli in the light of two different test cases.

Elasticity-to-Capillarity Transition in Soft Substrate Deformation

Whereas capillarity controls fluid dynamics at submillimeter scale and elasticity determines the mechanics of rigid solids, their coupling governs elastocapillary deformations on soft solids. Here, we directly probed the deformations on soft substrates induced by sessile nanodroplets. The wetting ridge created around the contact line and the dimple formed underneath the nanodroplet were imaged with a high spatial resolution using atomic force microscopy. The ridge height nonmonotonically depends on the substrate stiffness, and the dimple depth nonlinearly depends on the droplet size. The capillarity of the substrate overcomes the elasticity of the substrate in dominating the deformations when the elastocapillary length is approximately larger than the droplet contact radius, showing an experimental observation of the elasticity-to-capillarity transition. This study provides an experimental approach to investigate nanoscale elastocapillarity, and the insights have the potential to kick-off future work on the fundamentals of solid mechanics.

A model for 3D deformation and reconstruction of contractile microtissues

Tissue morphogenesis and regeneration are essentially mechanical processes that involve coordination of cellular forces, production and structural remodeling of extracellular matrix (ECM), and cell migration. Discovering the principles of cell-ECM interactions and tissue-scale deformation in mechanically-loaded tissues is instrumental to the development of novel regenerative therapies. The combination of high-throughput three-dimensional (3D) culture systems and experimentally-validated computational models accelerate the study of these principles. In our previous work [E. Mailand, et al., Biophys. J., 2019, 117, 975-986], we showed that prominent surface stresses emerge in constrained fibroblast-populated collagen gels, driving the morphogenesis of fibrous microtissues. Here, we introduce an active material model that allows the embodiment of surface and bulk contractile stresses while maintaining the passive elasticity of the ECM in a 3D setting. Unlike existing models, the stresses are driven by mechanosensing and not by an externally applied signal. The mechanosensing component is incorporated in the model through a direct coupling of the local deformation state with the associated contractile force generation. Further, we propose a finite element implementation to account for large deformations, nonlinear active material response, and surface effects. Simulation results quantitatively capture complex shape changes during tissue formation and as a response to surgical disruption of tissue boundaries, allowing precise calibration of the parameters of the 3D model. The results of this study imply that the organization of the extracellular matrix in the bulk of the tissue may not be a major factor behind the morphogenesis of fibrous tissues at sub-millimeter length scales.

A surface flattening method for characterizing the surface stress, drained Poisson's ratio and diffusivity of poroelastic gels

When a poroelastic gel is released from a patterned mold, surface stress drives deformation and solvent migration in the gel and flattens its surface profile in a time-dependent manner. Specifically, the gel behaves like an incompressible solid immediately after removal from the mold, and becomes compressible as the solvent is able to squeeze out of the polymer network. In this work, we use the finite element method (FEM) to simulate this transient surface flattening process. We assume that the surface stress is isotropic and constant, the polymer network is linearly elastic and isotropic, and that solvent flow obeys Darcy's law. The short-time and long-time surface profiles can be used to determine the surface stress and drained Poisson's ratio of the gel. Our analysis shows that the drained Poisson's ratio and the diffusivity of the gel can be obtained using interferometry and high-speed video microscopy, without mechanical measurement.

Molecular Characterization of Polymer Networks

Polymer networks are complex systems consisting of molecular components. Whereas the properties of the individual components are typically well understood by most chemists, translating that chemical insight into polymer networks themselves is limited by the statistical and poorly defined nature of network structures. As a result, it is challenging, if not currently impossible, to extrapolate from the molecular behavior of components to the full range of performance and properties of the entire polymer network. Polymer networks therefore present an unrealized, important, and interdisciplinary opportunity to exert molecular-level, chemical control on material macroscopic properties. A barrier to sophisticated molecular approaches to polymer networks is that the techniques for characterizing the molecular structure of networks are often unfamiliar to many scientists. Here, we present a critical overview of the current characterization techniques available to understand the relation between the molecular properties and the resulting performance and behavior of polymer networks, in the absence of added fillers. We highlight the methods available to characterize the chemistry and molecular-level properties of individual polymer strands and junctions, the gelation process by which strands form networks, the structure of the resulting network, and the dynamics and mechanics of the final material. The purpose is not to serve as a detailed manual for conducting these measurements but rather to unify the underlying principles, point out remaining challenges, and provide a concise overview by which chemists can plan characterization strategies that suit their research objectives. Because polymer networks cannot often be sufficiently characterized with a single method, strategic combinations of multiple techniques are typically required for their molecular characterization.