Efficient evaluation of the material response of tissues reinforced by statistically oriented fibres

For several classes of soft biological tissues, modelling complexity is in part due to the arrangement of the collagen fibres. In general, the arrangement of the fibres can be described by defining, at each point in the tissue, the structure tensor (i.e. the tensor product of the unit vector of the local fibre arrangement by itself) and a probability distribution of orientation. In this approach, assuming that the fibres do not interact with each other, the overall contribution of the collagen fibres to a given mechanical property of the tissue can be estimated by means of an averaging integral of the constitutive function describing the mechanical property at study over the set of all possible directions in space. Except for the particular case of fibre constitutive functions that are polynomial in the transversely isotropic invariants of the deformation, the averaging integral cannot be evaluated directly, in a single calculation because, in general, the integrand depends both on deformation and on fibre orientation in a non-separable way. The problem is thus, in a sense, analogous to that of solving the integral of a function of two variables, which cannot be split up into the product of two functions, each depending only on one of the variables. Although numerical schemes can be used to evaluate the integral at each deformation increment, this is computationally expensive. With the purpose of containing computational costs, this work proposes approximation methods that are based on the direct integrability of polynomial functions and that do not require the step-by-step evaluation of the averaging integrals. Three different methods are proposed: (a) a Taylor expansion of the fibre constitutive function in the transversely isotropic invariants of the deformation; (b) a Taylor expansion of the fibre constitutive function in the structure tensor; (c) for the case of a fibre constitutive function having a polynomial argument, an approximation in which the directional average of the constitutive function is replaced by the constitutive function evaluated at the directional average of the argument. Each of the proposed methods approximates the averaged constitutive function in such a way that it is multiplicatively decomposed into the product of a function of the deformation only and a function of the structure tensors only. In order to assess the accuracy of these methods, we evaluate the constitutive functions of the elastic potential and the Cauchy stress, for a biaxial test, under different conditions, i.e. different fibre distributions and different ratios of the nominal strains in the two directions. The results are then compared against those obtained for an averaging method available in the literature, as well as against the integration made at each increment of deformation

The linear elasticity tensor of incompressibile materials

With a universally accepted abuse of terminology, materials having much larger stiffness for volumetric than for shear deformations are called incompressible. This work proposes two approaches for the evaluation of the correct form of the linear elasticity tensor of so-called incompressible materials, both stemming from non-linear theory. In the approach of strict incompressibility, one imposes the kinematical constraint of isochoric deformation. In the approach of quasi-incompressibility, which is often employed to enforce incompressibility in numerical applications such as the Finite Element Method, one instead assumes a decoupled form of the elastic potential (or strain energy), which is written as the sum of a function of the volumetric deformation only and a function of the distortional deformation only, and then imposes that the bulk modulus be much larger than all other moduli. The conditions which the elasticity tensor has to obey for both strict incompressibility and quasi-incompressibility have been derived, regardless of the material symmetry. The representation of the linear elasticity tensor for the quasi-incompressible case differs from that of the strictly incompressible case by one parameter, which can be conveniently chosen to be the bulk modulus. Some important symmetries have been studied in detail, showing that the linear elasticity tensors for the cases of isotropy, transverse isotropy and orthotropy are characterised by one, three and six independent parameters, respectively, for the case of strict incompressibility, and two, four and seven independent parameters, respectively, for the case of quasi-incompressibility, as opposed to the two, five and nine parameters, respectively, of the general compressible cas

Effective properties of hierarchical fiber-reinforced composites via a three-scale asymptotic homogenization approach

The study of the properties of multiscale composites is of great interest in engineering and biology. Particularly, hierarchical composite structures can be found in nature and in engineering. During the past decades, the multiscale asymptotic homogenization technique has shown its potential in the description of such composites by taking advantage of their characteristics at the smaller scales, ciphered in the so-called effective coefficients. Here, we extend previous works by studying the in-plane and out-of-plane effective properties of hierarchical linear elastic solid composites via a three-scale asymptotic homogenization technique. In particular, the approach is adjusted for a multiscale composite with a square-symmetric arrangement of uniaxially aligned cylindrical fibers, and the formulae for computing its effective properties are provided. Finally, we show the potential of the proposed asymptotic homogenization procedure by modeling the effective properties of musculoskeletal mineralized tissues, and we compare the results with theoretical and experimental data for bone and tendon tissues

Synaptic boutons sizes are tuned to best fit their physiological performances

To truly appreciate the myriad of events which relate synaptic function and vesicle dynamics, simulations should be done in a spatially realistic environment. This holds true in particular in order to explain as well the rather astonishing motor patterns which we observed within in vivo recordings which underlie peristaltic contractionsas well as the shape of the EPSPs at different forms of long-term stimulation, presented both here, at a well characterized synapse, the neuromuscular junction (NMJ) of the Drosophila larva (c.f. Figure 1). To this end, we have employed a reductionist approach and generated three dimensional models of single presynaptic boutons at the Drosophila larval NMJ. Vesicle dynamics are described by diffusion-like partial differential equations which are solved numerically on unstructured grids using the uG platform. In our model we varied parameters such as bouton-size, vesicle output probability (Po), stimulation frequency and number of synapses, to observe how altering these parameters effected bouton function. Hence we demonstrate that the morphologic and physiologic specialization maybe a convergent evolutionary adaptation to regulate the trade off between sustained, low output, and short term, high output, synaptic signals. There seems to be a biologically meaningful explanation for the co-existence of the two different bouton types as previously observed at the NMJ (characterized especially by the relation between size and Po), the assigning of two different tasks with respect to short- and long-time behaviour could allow for an optimized interplay of different synapse types. We can present astonishing similar results of experimental and simulation data which could be gained in particular without any data fitting, however based only on biophysical values which could be taken from different experimental results. As a side product, we demonstrate how advanced methods from numerical mathematics could help in future to resolve also other difficult experimental neurobiological issues