The Föppl–von Kármán equations of elastic plates with initial stress

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\"{o}ppl-von K\'{a}rm\'{a}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 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 4D printing

Rods coiling about a rigid constraint: helices and perversions

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

Mechanics of axisymmetric sheets of interlocking and slidable rods

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

Elastic Instability behind Brittle Fracture

International audienceWe argue that nucleation of brittle cracks in initially flawless soft elastic solids is preceded by a nonlinear elastic instability, which cannot be captured without accounting for geometrical precise description of finite elastic deformation. As a prototypical problem we consider a homogeneous elastic body subjected to tension and assume that it is weakened by the presence of a free surface which then serves as a site of crack nucleation. We show that in this maximally simplified setting, brittle fracture emerges from a symmetry breaking elastic instability activated by softening and involving large elastic rotations. The implied bifurcation of the homogeneous elastic equilibrium is highly unconventional for nonlinear elasticity as it exhibits an extraordinary sensitivity to geometry, reminiscent of the transition to turbulence in fluids. We trace the post-bifurcational development of this instability beyond the limits of applicability of scale free continuum elasticity and use a phase-field approach to capture the scale dependent sub-continuum strain localization, signaling the formation of actual cracks

Mathematical modelling of axonal cortex contractility

The axonal cortex is composed of a regular structure of F-actin and spectrin able to contract thanks to myosin II motors. Such an active tension is of fundamental importance in controlling the physiological shape of axons. Recent experiments show that axons modulate the contraction of the cortex when subject to mechanical deformations, exhibiting a non-trivial coupling between the hoop and the axial active tension. However, the underlying mechanisms are still poorly understood. In this paper, we propose a continuum model of the axon based on the active strain theory. By using the Coleman–Noll procedure, we shed light on the coupling between the hoop and the axial active strain through the Mandel stress tensor. We propose a qualitative analysis of the system under the simplifying assumption of incompressibility, showing the existence of a stable equilibrium solution. In particular, our results show that the axon regulates the active contraction to maintain a homeostatic stress state. Finally, we propose numerical simulations of the model, using a more suitable compressible constitutive law. The results are compared with experimental data, showing an excellent quantitative agreement

Active elasticity drives the formation of periodic beading in damaged axons

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 nonlinear evolution of the bifurcated branches is investigated through finite element simulations

Rayleigh-Taylor instability in soft elastic layers

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'

Surface tension controls the onset of gyrification in brain organoids

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