Search

We report surprising morphological changes of suspension droplets (containing class II hydrophobin protein HFBI from Trichoderma reesei in water) as they evaporate with a contact line pinned on a rigid solid substrate. Both pendant and sessile droplets display the formation of an encapsulating elastic film as the bulk concentration of solute reaches a critical value during evaporation, but the morphology of the droplet varies significantly: for sessile droplets, the elastic film ultimately crumples in a nearly flattened area close to the apex while in pendant droplets, circumferential wrinkling occurs close to the contact line. These different morphologies are understood through a gravito-elastocapillary model that predicts the droplet morphology and the onset of shape changes, as well as showing that the influence of the direction of gravity remains crucial even for very small droplets (where the effect of gravity can normally be neglected). The results pave the way to control droplet shape in several engineering and biomedical applications.

Initially, stressed plates are widely used in modern fabrication techniques, such as additive manufacturing and UV lithography, for their tunable morphology by application of external stimuli. In this work, we propose a formal asymptotic derivation of the Föppl–von Kármán equations for an elastic plate with initial stresses, using the constitutive theory of nonlinear elastic solids with initial stresses under the assumptions of incompressibility and material isotropy. Compared to existing works, our approach allows us to determine the morphological transitions of the elastic plate without prescribing the underlying target metric of the unstressed state of the elastic body. We explicitly solve the derived FvK equations in some physical problems of engineering interest, discussing how the initial stress distribution drives the emergence of spontaneous curvatures within the deformed plate. The proposed mathematical framework can be used to tailor shape on demand, with applications in several engineering fields ranging from soft robotics to four-dimensional printing.

In several pathological conditions, such as coronavirus infections, multiple sclerosis, Alzheimer's and Parkinson's diseases, the physiological shape of axons is altered and a periodic sequence of bulges appears. Experimental evidences suggest that such morphological changes are caused by the disruption of the microtubules composing the cytoskeleton of the axon. In this paper, we develop a mathematical model of damaged axons based on the theory of continuum mechanics and nonlinear elasticity. The axon is described as a cylinder composed of an inner passive part, called axoplasm, and an outer active cortex, composed mainly of F-actin and able to contract thanks to myosin-II motors. Through a linear stability analysis we show that, as the shear modulus of the axoplasm diminishes due to the disruption of the cytoskeleton, the active contraction of the cortex makes the cylindrical configuration unstable to axisymmetric perturbations, leading to a beading pattern. Finally, the non-linear evolution of the bifurcated branches is investigated through finite element simulations.

Mechanical instabilities can be exploited to design innovative structures, able to change their shape in the presence of external stimuli. In this work, we derive a mathematical model of an elastic beam subjected to an axial force and constrained to smoothly slide along a rigid support, where the distance between the rod midline and the constraint is fixed and finite. Using both theoretical and computational techniques, we characterize the bifurcations of such a mechanical system, in which the axial force and the natural curvature of the beam are used as control parameters. We show that, in the presence of a straight support, the rod can deform into shapes exhibiting helices and perversions, namely transition zones connecting together two helices with opposite chirality. The mathematical predictions of the proposed model are also compared with some experiments, showing a good quantitative agreement. In particular, we find that the buckled configurations may exhibit multiple perversions and we propose a possible explanation for this phenomenon based on the energy landscape of the mechanical system.

© 2020 Elsevier. This manuscript version is made available under the CC-BY-NC-ND 4.0 license http://creativecommons.org/licenses/by-nc-nd/4.0/ In this work, we study the mechanics of metamaterial sheets inspired by the pellicle of Euglenids. They are composed of interlocking elastic rods which can freely slide along their edges. We characterize the kinematics and the mechanics of these structures using the special Cosserat theory of rods and by assuming axisymmetric deformations of the tubular assembly. Through an asymptotic expansion, we investigate both structures that comprise a discrete number of rods and the limit case of a sheet composed by infinitely many rods. We apply our theoretical framework to investigate the stability of these structures in the presence of an axial load. Through a linear analysis, we compute the critical buckling force for both the discrete and the continuous case. For the latter, we also perform a numerical post-buckling analysis, studying the non-linear evolution of the bifurcation through finite elements simulations.

The mathematical modeling of the contraction of a muscle is a crucial problem in biomechanics. Several different models of muscle activation exist in literature. A possible approach to contractility is the so-called active strain: it is based on a multiplicative decomposition of the deformation gradient into an active contribution, accounting for the muscle activation, and an elastic one, due to the passive deformation of the body. We show that the active strain approach does not allow to recover the experimental stress–stretch curve corresponding to a uniaxial deformation of a skeletal muscle, whatever the functional form of the strain energy. To overcome such difficulty, we introduce an alternative model, that we call mixture active strain approach, where the muscle is composed of two different solid phases and only one of them actively contributes to the active behavior of the muscle.

Active materials are media for which deformations can occur in absence of loads, given an external stimulus. Two approaches to the modeling of such materials are mainly used in literature, both based on the introduction of a new tensor: an additive stress $\mathsf{P}_{\text{act}}$ in the active stress case and a multiplicative strain $\mathsf{F}_{a}$ in the active strain one. Aim of this paper is the comparison between the two approaches on simple shears. Considering an incompressible and transversely isotropic material, we design constitutive relations for $\mathsf{P}_{\text{act}}$ and $\mathsf{F}_{a}$ so that they produce the same results for a uniaxial deformation along the symmetry axis. We then study the two approaches in the case of a simple shear deformation. In a hyperelastic setting, we show that the two approaches produce different stress components along a simple shear, unless some necessary conditions on the strain energy density are fulfilled. However, such conditions are very restrictive and rule out the usual elastic strain energy functionals. Active stress and active strain therefore produce different results in shear, even if they both fit uniaxial data. Our results show that experimental data on the stress-stretch response on uniaxial deformations are not enough to establish which activation approach can capture better the mechanics of active materials. We conclude that other types of deformations, beyond the uniaxial one, should be taken into consideration in the modeling of such materials.

Understanding the mechanics of brain embryogenesis can provide insights on pathologies related to brain development, such as lissencephaly, a genetic disease which causes a reduction of the number of cerebral sulci. Recent experiments on brain organoids have confirmed that gyrification, i.e. the formation of the folded structures of the brain, is triggered by the inhomogeneous growth of the peripheral region. However, the rheology of these cellular aggregates and the mechanics of lissencephaly are still matter of debate. In this work, we develop a mathematical model of brain organoids based on the theory of morpho-elasticity. We describe them as non-linear elastic bodies, composed of a disk surrounded by a growing layer called cortex. The external boundary is subjected to a tissue surface tension due the intercellular adhesion forces. We show that the resulting surface energy is relevant at the small length scales of brain organoids and affects the mechanics of cellular aggregates. We perform a linear stability analysis of the radially symmetric configuration and we study the post-buckling behaviour through finite element simulations. We find that the process of gyrification is triggered by the cortex growth and modulated by the competition between two length scales: the radius of the organoid and the capillary length generated by surface tension. We show that a solid model can reproduce the results of the in-vitro experiments. Furthermore, we prove that the lack of brain sulci in lissencephaly is caused by a reduction of the cell stiffness: the softening of the organoid strengthens the role of surface tension, delaying or even inhibiting the onset of a mechanical instability at the free boundary.

A soft solid is said to be initially stressed if it is subjected to a state of internal stress in its unloaded reference configuration. In physical terms, its stored elastic energy may not vanish in the absence of an elastic deformation, being also dependent on the spatial distribution of the underlying material inhomogeneities. Developing a sound mathematical framework to model initially stressed solids in nonlinear elasticity is key for many applications in engineering and biology. This work investigates the links between the existence of elastic minimizers and the constitutive restrictions for initially stressed materials subjected to finite deformations. In particular, we consider a subclass of constitutive responses in which the strain energy density is taken as a scalar-valued function of both the deformation gradient and the initial stress tensor. The main advantage of this approach is that the initial stress tensor belongs to the group of divergence-free symmetric tensors satisfying the boundary conditions in any given reference configuration. However, it is still unclear which physical restrictions must be imposed for the well-posedness of this elastic problem. Assuming that the constitutive response depends on the choice of the reference configuration only through the initial stress tensor, under given conditions we prove the local existence of a relaxed state given by an implicit tensor function of the initial stress distribution. This tensor function is generally not unique, and can be transformed according to the symmetry group of the material at fixed initial stresses. These results allow one to extend Ball's existence theorem of elastic minimizers for the proposed constitutive choice of initially stressed materials. This article is part of the theme issue ‘Rivlin's legacy in continuum mechanics and applied mathematics’.

Solid tumours have the ability to assemble their own vascular network for optimizing their access to the vital nutrients. These new capillaries are morphologically different from normal physiological vessels. In particular, they have a much higher spatial tortuosity forcing an impaired flow within the peritumoral area. This is a major obstacle for the efficient delivery of antitumoral drugs. This work proposes a morpho-elastic model of the tumour vessels. A tumour capillary is considered as a growing hyperelastic tube that is spatially constrained by a linear elastic environment, representing the interstitial matter. We assume that the capillary is an incompressible neo-Hookean material, whose growth is modelled using a multiplicative decomposition of the deformation gradient. We study the morphological stability of the capillary by means of the method of incremental deformations superposed on finite strains, solving the corresponding incremental problem using the Stroh formulation and the impedance matrix method. The incompatible axial growth of the straight capillary is found to control the onset of a bifurcation towards a tortuous shape. The post-buckling morphology is studied using a mixed finite element formulation in the fully nonlinear regime. The proposed model highlights how the geometrical and the elastic properties of the capillary and the surrounding medium concur to trigger the loss of marginal stability of the straight capillary and the nonlinear development of its spatial tortuosity.

Residual stresses may appear in elastic bodies, owing to the formation of misfits in the microstructure, driven by plastic deformations or thermal or growth processes. They are especially widespread in living matter, resulting from dynamic remodelling processes aimed at optimizing the overall structural response to environmental physical forces. From a mechanical viewpoint, residual stresses are classically modelled through the introduction of a virtual incompatible state that collects the local relaxed states around each material point. In this work, we employ an alternative approach based on a strain energy function that constitutively depends only on the deformation gradient and the residual stress tensor. In particular, our objective is to study the morphological stability of an incompressible sphere, made of a neo-Hookean material, and subjected to given distributions of residual stresses. The boundary value elastic problem is studied with analytic and numerical tools. Firstly, we perform a linear stability analysis on the prestressed solid sphere using the method of incremental deformations. The marginal stability conditions are given as a function of a control parameter, which is the dimensionless variable that represents the characteristic intensity of the residual stresses. Secondly, we perform finite-element simulations using a mixed formulation in order to investigate the postbuckling morphology in the fully nonlinear regime. Considering different initial distributions of the residual stresses, we find that different morphological transitions occur around the material domain, where the hoop residual stress reaches its maximum compressive value. The loss of spherical symmetry is found to be controlled by the mechanical and geometrical properties of the sphere, as well as the spatial distribution of the residual stress. The results provide useful guidelines for designing morphable soft spheres, for example by controlling residual stresses through active deformations. They finally suggest a viable solution for the nondestructive characterization of residual stresses in soft tissues, such as solid tumours.

This work investigates the morphological stability of a soft body composed of two heavy elastic layers attached to a rigid surface and subjected only to the bulk gravity force. Using theoretical and computational tools, we characterize the selection of different patterns as well as their nonlinear evolution, unveiling the interplay between elastic and geometric effects for their formation. Unlike similar gravity-induced shape transitions in fluids, such as the Rayleigh–Taylor instability, we prove that the nonlinear elastic effects saturate the dynamic instability of the bifurcated solutions, displaying a rich morphological diagram where both digitations and stable wrinkling can emerge. The results of this work provide important guidelines for the design of novel soft systems with tunable shapes, with several applications in engineering sciences. This article is part of the themed issue ‘Patterning through instabilities in complex media: theory and applications.’