Continuum soft tissue models from upscaling of arrays of hyperelastic cells

Abstract

Constitutive models for soft tissue mechanics are typically constructed by fitting phenomenological models to experimental measurements. However, a significant challenge is to rationally construct soft tissue models that encode the properties of the constituent cells and their extracellular matrix. This work presents a framework to derive multiscale soft tissue models that incorporate properties of individual cells without assuming homogeneity or periodicity at the cell level. We consider a viscoelastic model for each cell (which can deform, grow and divide), that we couple to form a network description of a one-dimensional line of cells. We use a discrete-to-continuum approach to form (nonlinear) continuum partial differential equation models for the tissue. These models elucidate the contrasting role of the two forms of dissipation: substrate dissipation localizes the deformation to the neighbourhood of the free boundary and inhibits morphoelastic growth, whereas internal cell dissipation promotes spatial uniformity and does not influence the elongation length. Furthermore, cell division is shown to increase the rate of elongation of the array compared with growth alone, provided the substrate dissipation is proportional to the cell surface area.