Validity of the Cauchy-Born rule applied to discrete cellular-scale models of biological tissues.

The development of new models of biological tissues that consider cells in a discrete manner is becoming increasingly popular as an alternative to continuum methods based on partial differential equations, although formal relationships between the discrete and continuum frameworks remain to be established. For crystal mechanics, the discrete-to-continuum bridge is often made by assuming that local atom displacements can be mapped homogeneously from the mesoscale deformation gradient, an assumption known as the Cauchy-Born rule (CBR). Although the CBR does not hold exactly for noncrystalline materials, it may still be used as a first-order approximation for analytic calculations of effective stresses or strain energies. In this work, our goal is to investigate numerically the applicability of the CBR to two-dimensional cellular-scale models by assessing the mechanical behavior of model biological tissues, including crystalline (honeycomb) and noncrystalline reference states. The numerical procedure involves applying an affine deformation to the boundary cells and computing the quasistatic position of internal cells. The position of internal cells is then compared with the prediction of the CBR and an average deviation is calculated in the strain domain. For center-based cell models, we show that the CBR holds exactly when the deformation gradient is relatively small and the reference stress-free configuration is defined by a honeycomb lattice. We show further that the CBR may be used approximately when the reference state is perturbed from the honeycomb configuration. By contrast, for vertex-based cell models, a similar analysis reveals that the CBR does not provide a good representation of the tissue mechanics, even when the reference configuration is defined by a honeycomb lattice. The paper concludes with a discussion of the implications of these results for concurrent discrete and continuous modeling, adaptation of atom-to-continuum techniques to biological tissues, and model classification

A 3D discrete model of the diaphragm and human trunk

In this paper, a 3D discrete model is presented to model the movements of the trunk during breathing. In this model, objects are represented by physical particles on their contours. A simple notion of force generated by a linear actuator allows the model to create forces on each particle by way of a geometrical attractor. Tissue elasticity and contractility are modeled by local shape memory and muscular fibers attractors. A specific dynamic MRI study was used to build a simple trunk model comprised of by three compartments: lungs, diaphragm and abdomen. This model was registered on the real geometry. Simulation results were compared qualitatively as well as quantitatively to the experimental data, in terms of volume and geometry. A good correlation was obtained between the model and the real data. Thanks to this model, pathology such as hemidiaphragm paralysis can also be simulated.Comment: published in: "Lung Modelling", France (2006

Colloquium: Mechanical formalisms for tissue dynamics

The understanding of morphogenesis in living organisms has been renewed by tremendous progressin experimental techniques that provide access to cell-scale, quantitative information both on theshapes of cells within tissues and on the genes being expressed. This information suggests that ourunderstanding of the respective contributions of gene expression and mechanics, and of their crucialentanglement, will soon leap forward. Biomechanics increasingly benefits from models, which assistthe design and interpretation of experiments, point out the main ingredients and assumptions, andultimately lead to predictions. The newly accessible local information thus calls for a reflectionon how to select suitable classes of mechanical models. We review both mechanical ingredientssuggested by the current knowledge of tissue behaviour, and modelling methods that can helpgenerate a rheological diagram or a constitutive equation. We distinguish cell scale ("intra-cell")and tissue scale ("inter-cell") contributions. We recall the mathematical framework developpedfor continuum materials and explain how to transform a constitutive equation into a set of partialdifferential equations amenable to numerical resolution. We show that when plastic behaviour isrelevant, the dissipation function formalism appears appropriate to generate constitutive equations;its variational nature facilitates numerical implementation, and we discuss adaptations needed in thecase of large deformations. The present article gathers theoretical methods that can readily enhancethe significance of the data to be extracted from recent or future high throughput biomechanicalexperiments.Comment: 33 pages, 20 figures. This version (26 Sept. 2015) contains a few corrections to the published version, all in Appendix D.2 devoted to large deformation

A novel model for one-dimensional morphoelasticity. Part I - Theoretical foundations

While classical continuum theories of elasticity and viscoelasticity have long been used to describe the mechanical behaviour of solid biological tissues, they are of limited use for the description of biological tissues that undergo continuous remodelling. The structural changes to a soft tissue associated with growth and remodelling require a mathematical theory of ‘morphoelasticity’ that is more akin to plasticity than elasticity. However, previously-derived mathematical models for plasticity are difficult to apply and interpret in the context of growth and remodelling: many important concepts from the theory of plasticity do not have simple analogues in biomechanics.\ud \ud In this work, we describe a novel mathematical model that combines the simplicity and interpretability of classical viscoelastic models with the versatility of plasticity theory. While our focus here is on one-dimensional problems, our model builds on earlier work based on the multiplicative decomposition of the deformation gradient and can be adapted to develop a three-dimensional theory. The foundation of this work is the concept of ‘effective strain’, a measure of the difference between the current state and a hypothetical state where the tissue is mechanically relaxed. We develop one-dimensional equations for the evolution of effective strain, and discuss a number of potential applications of this theory. One significant application is the description of a contracting fibroblast-populated collagen lattice, which we further investigate in Part II

Validity of the Cauchy-Born rule applied to discrete cellular-scale models of biological tissues

The development of new models of biological tissues that consider cells in a discrete manner is becoming increasingly popular as an alternative to PDE-based continuum methods, although formal relationships between the discrete and continuum frameworks remain to be established. For crystal mechanics, the discrete-to-continuum bridge is often made by assuming that local atom displacements can be mapped homogeneously from the mesoscale deformation gradient, an assumption known as the Cauchy-Born rule (CBR). Although the CBR does not hold exactly for non-crystalline materials, it may still be used as a first order approximation for analytic calculations of effective stresses or strain energies. In this work, our goal is to investigate numerically the applicability of the CBR to 2-D cellular-scale models by assessing the mechanical behaviour of model biological tissues, including crystalline (honeycomb) and non-crystalline reference states. The numerical procedure consists in precribing an affine deformation on the boundary cells and computing the position of internal cells. The position of internal cells is then compared with the prediction of the CBR and an average deviation is calculated in the strain domain. For centre-based models, we show that the CBR holds exactly when the deformation gradient is relatively small and the reference stress-free configuration is defined by a honeycomb lattice. We show further that the CBR may be used approximately when the reference state is perturbed from the honeycomb configuration. By contrast, for vertex-based models, a similar analysis reveals that the CBR does not provide a good representation of the tissue mechanics, even when the reference configuration is defined by a honeycomb lattice. The paper concludes with a discussion of the implications of these results for concurrent discrete/continuous modelling, adaptation of atom-to-continuum (AtC) techniques to biological tissues and model classification