Matched asymptotic solution for crease nucleation in soft solids

A soft solid subjected to a large compression develops sharp self-contacting folds at its free surface, known as creases. Creasing is physically different from structural elastic instabilities, like buckling or wrinkling. Indeed, it is a fully nonlinear material instability, similar to a phase-transformation. This work provides theoretical insights of the physics behind crease nucleation. Creasing is proved to occur after a global bifurcation allowing the co-existence of an outer deformation and an inner solution with localised self-contact at the free surface. The most fundamental result here is the analytic prediction of the nucleation threshold, in excellent agreement with experiments and numerical simulations. A matched asymptotic solution is given within the intermediate region between the two co-existing states. The self-contact acts like the point-wise disturbance in the Oseen's correction for the Stokes flow past a circle. Analytic expressions of the matching solution and its range of validity are also derived

Geometric control by active mechanics of epithelial gap closure

Epithelial wound healing is one of the most important biological processes occurring during the lifetime of an organism. It is a self-repair mechanism closing wounds or gaps within tissues to restore their functional integrity. In this work we derive a new diffuse interface approach for modelling the gap closure by means of a variational principle in the framework of non-equilibrium thermodynamics. We investigate the interplay between the crawling with lamellipodia protrusions and the supracellular tension exerted by the actomyosin cable on the closure dynamics. These active features are modeled as Korteweg forces into a generalised chemical potential. From an asymptotic analysis, we derive a pressure jump across the gap edge in the sharp interface limit. Moreover, the chemical potential diffuses as a Mullins-Sekerka system, and its interfacial value is given by a Gibbs-Thompson relation for its local potential driven by the curvature-dependent purse-string tension. The finite element simulations show an excellent quantitative agreement between the closure dynamics and the morphology of the edge with respect to existing biological experiments. The resulting force patterns are also in good qualitative agreement with existing traction force microscopy measurements. Our results shed light on the geometrical control of the gap closure dynamics resulting from the active forces that are chemically activated around the gap edge.Shedding light on the geometric control of the gap closure dynamics in epithelial wound healing through a novel diffuse interface mathematical model derived by means of a variational principle in the framework of non-equilibrium thermodynamics

Mass transport in morphogenetic processes: a second gradient theory for volumetric growth and material remodeling

International audienceIn this work, we derive a novel thermo-mechanical theory for growth and remodeling of biological materials in morphogenetic processes. This second gradient hyperelastic theory is the first attempt to describe both volumetric growth and mass transport phenomena in a single-phase continuum model, where both stress- and shape-dependent growth regulations can be investigated. The diffusion of biochemical species (e.g. morphogens, growth factors, migration signals) inside the material is driven by configurational forces, enforced in the balance equations and in the set of constitutive relations. Mass transport is found to depend both on first- and on second-order material connections, possibly withstanding a chemotactic behavior with respect to diffusing molecules. We find that the driving forces of mass diffusion can be written in terms of covariant material derivatives reflecting, in a purely geometrical manner, the presence of a (first-order) torsion and a (second-order) curvature. Thermodynamical arguments show that the Eshelby stress and hyperstress tensors drive the rearrangement of the first- and second-order material inhomogeneities, respectively. In particular, an evolution law is proposed for the first-order transplant, extending a well-known result for inelastic materials. Moreover, we define the first stress-driven evolution law of the second-order transplant in function of the completely material Eshelby hyperstress

Pattern selection in growing tubular tissues

International audienceTubular organs display a wide variety of surface morphologies including circumferential and longitudinal folds, square and hexagonal undulations, and finger-type protrusions. Surface morphology is closely correlated to tissue function and serves as a clinical indicator for physiological and pathological conditions, but the regulators of surface morphology remain poorly understood. Here, we explore the role of geometry and elasticity on the formation of surface patterns. We establish morphological phase diagrams for patterns selection and show that increasing the thickness or stiffness ratio between the outer and inner tubular layers induces a gradual transition from circumferential to longitudinal folding. Our results suggest that physical forces act as regulators during organogenesis and give rise to the characteristic circular folds in the esophagus, the longitudinal folds in the valves of Kerckring, the surface networks in villi, and the crypts in the large intestine

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

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

3D-Calibration of the IMU

International audienceA new calibration method for Inertial Measurement Unit (IMU) of strapdown inertial technology was presented. IMU has been composed of accelerometers, gyroscopes and a circuit of signal processing. Normally, a rate transfer test and multi-position tests are used for IMU calibration. The new calibration method is based on whole angle rotation or finite rotation. In fact it is suggested to turn over IMU around three axes simultaneously. In order to solve the equation of calibration, it is necessary to provide an equality of a rank of basic matrix into degree of basic matrix. The results of simulated IMU data presented to demonstrate the performance of the new calibration method