Structural Optimisation: Biomechanics of the Femur

A preliminary iterative 3D meso-scale structural model of the femur was developed, in which bar and shell elements were used to represent trabecular and cortical bone respectively. The cross-sectional areas of the bar elements and the thickness values of the shell elements were adjusted over successive iterations of the model based on a target strain stimulus, resulting in an optimised construct. The predicted trabecular architecture, and cortical thickness distribution showed good agreement with clinical observations, based on the application of a single leg stance load case during gait. The benefit of using a meso-scale structural approach in comparison to micro or macro-scale continuum approaches to predictive bone modelling was achievement of the symbiotic goals of computational efficiency and structural description of the femur.Comment: Accepted by Engineering and Computational Mechanics (Proceedings of the ICE

New Mechanics of Traumatic Brain Injury

The prediction and prevention of traumatic brain injury is a very important aspect of preventive medical science. This paper proposes a new coupled loading-rate hypothesis for the traumatic brain injury (TBI), which states that the main cause of the TBI is an external Euclidean jolt, or SE(3)-jolt, an impulsive loading that strikes the head in several coupled degrees-of-freedom simultaneously. To show this, based on the previously defined covariant force law, we formulate the coupled Newton-Euler dynamics of brain's micro-motions within the cerebrospinal fluid and derive from it the coupled SE(3)-jolt dynamics. The SE(3)-jolt is a cause of the TBI in two forms of brain's rapid discontinuous deformations: translational dislocations and rotational disclinations. Brain's dislocations and disclinations, caused by the SE(3)-jolt, are described using the Cosserat multipolar viscoelastic continuum brain model. Keywords: Traumatic brain injuries, coupled loading-rate hypothesis, Euclidean jolt, coupled Newton-Euler dynamics, brain's dislocations and disclinationsComment: 18 pages, 1 figure, Late

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

A Muscle Model Based on Feldman's Lambda Model: 3D Finite Element Implementation

This paper presents the introduction of Feldman's muscle model in a three dimensional continuum finite element model of the human face. This model is compared to the classical Hill-type muscle modelComment: CMBBE'2013, Salt Lake City : United States (2013

A multiscale mechanobiological model of bone remodelling predicts site-specific bone loss in the femur during osteoporosis and mechanical disuse

We propose a multiscale mechanobiological model of bone remodelling to investigate the site-specific evolution of bone volume fraction across the midshaft of a femur. The model includes hormonal regulation and biochemical coupling of bone cell populations, the influence of the microstructure on bone turnover rate, and mechanical adaptation of the tissue. Both microscopic and tissue-scale stress/strain states of the tissue are calculated from macroscopic loads by a combination of beam theory and micromechanical homogenisation. This model is applied to simulate the spatio-temporal evolution of a human midshaft femur scan subjected to two deregulating circumstances: (i) osteoporosis and (ii) mechanical disuse. Both simulated deregulations led to endocortical bone loss, cortical wall thinning and expansion of the medullary cavity, in accordance with experimental findings. Our model suggests that these observations are attributable to a large extent to the influence of the microstructure on bone turnover rate. Mechanical adaptation is found to help preserve intracortical bone matrix near the periosteum. Moreover, it leads to non-uniform cortical wall thickness due to the asymmetry of macroscopic loads introduced by the bending moment. The effect of mechanical adaptation near the endosteum can be greatly affected by whether the mechanical stimulus includes stress concentration effects or not.Comment: 25 pages, 10 figure

Remodeling of biological tissue: Mechanically induced reorientation of a transversely isotropic chain network

A new class of micromechanically motivated chain network models for soft biological tissues is presented. On the microlevel, it is based on the statistics of long chain molecules. A wormlike chain model is applied to capture the behavior of the collagen microfibrils. On the macrolevel, the network of collagen chains is represented by a transversely isotropic eight chain unit cell introducing one characteristic material axis. Biomechanically induced remodeling is captured by allowing for a continuous reorientation of the predominant unit cell axis driven by a biomechanical stimulus. To this end, we adopt the gradual alignment of the unit cell axis with the direction of maximum principal strain. The evolution of the unit cell axis' orientation is governed by a first-order rate equation. For the temporal discretization of the remodeling rate equation, we suggest an exponential update scheme of Euler-Rodrigues type. For the spatial discretization, a finite element strategy is applied which introduces the current individual cell orientation as an internal variable on the integration point level. Selected model problems are analyzed to illustrate the basic features of the new model. Finally, the presented approach is applied to the biomechanically relevant boundary value problem of an in vitro engineered functional tendon construct.Comment: LaTeX2e, 19 pages, 9 figure

A FEM-experimental approach for the development of a conceptual linear actuator based on tendril's free coiling

Within the vastness of the plant species, certain living systems show tendril structures whose motion is of particular interest for biomimetic engineers. Tendrils sense and coil around suitable grips, and by shortening in length, they erect the remaining plant body. To achieve contraction, tendrils rotate along their main axis and shift from a linear to a double-spring geometry. This phenomenon is denoted as the free-coiling phase. In this work, with the aim of understanding the fundamentals of the mechanics behind the free coiling, a reverse-engineering approach based on the finite element method was firstly applied. The model consisted of an elongated cylinder with suitable material properties, boundary, and loading conditions, in order to reproduce the kinematics of the tendril. The simulation succeeded in mimicking coiling faithfully and was therefore used to validate a tentative linear actuator model based on the plant’s working principle. More in detail, exploiting shape memory alloy materials to obtain large reversible deformations, the main tendril features were implemented into a nickel-titanium spring-based testing model. The results of the experimental tests confirmed the feasibility of the idea in terms of both functioning principles and actual performance. It can be concluded that the final set-up can be used as a base for a prototype design of a new kind of a linear actuator