Simulation of planar soft tissues using a structural constitutive model: Finite element implementation and validation

Computational implementation of physical and physiologically realistic constitutive models is critical for numerical simulation of soft biological tissues in a variety of biomedical applications. It is well established that the highly nonlinear and anisotropic mechanical behaviors of soft tissues are an emergent behavior of the underlying tissue microstructure. In the present study, we have implemented a structural constitutive model into a finite element framework specialized for membrane tissues. We noted that starting with a single element subjected to uniaxial tension, the non-fibrous tissue matrix must be present to prevent unrealistic tissue deformations. Flexural simulations were used to set the non-fibrous matrix modulus because fibers have little effects on tissue deformation under three-point bending. Multiple deformation modes were simulated, including strip biaxial, planar biaxial with two attachment methods, and membrane inflation. Detailed comparisons with experimental data were undertaken to insure faithful simulations of both the macro-level stress–strain insights into adaptations of the fiber architecture under stress, such as fiber reorientation and fiber recruitment. Results indicated a high degree of fidelity and demonstrated interesting microstructural adaptions to stress and the important role of the underlying tissue matrix. Moreover, we apparently resolve a discrepancy in our 1997 study (Billiar and Sacks, 1997. J. Biomech. 30 (7), 753–756) where we observed that under strip biaxial stretch the simulated fiber splay responses were not in good agreement with the experimental results, suggesting non-affine deformations may have occurred. However, by correctly accounting for the isotropic phase of the measured fiber splay, good agreement was obtained. While not the final word, these simulations suggest that affine fiber kinematics for planar collagenous tissues is a reasonable assumption at the macro level. Simulation tools such as these are imperative in the design and simulation of native and engineered tissues.     

Finite element modeling of impulsive excitation and shear wave propagation in an incompressible,transversely isotropic medium

Elastic properties of materials can be measured by observing shear wave propagation following localized, impulsive excitations and relating the propagation velocity to a model of the material. However, characterization of anisotropic materials is difficult because of the number of elasticity constants in the material model and the complex dependence of propagation velocity relative to the excitation axis, material symmetries, and propagation directions. In this study, we develop a model of wave propagation following impulsive excitation in an incompressible, transversely isotropic (TI) material such as muscle. Wave motion is described in terms of three propagation modes identified by their polarization relative to the material symmetry axis and propagation direction. Phase velocities for these propagation modes are expressed in terms of five elasticity constants needed to describe a general TI material, and also in terms of three constants after the application of two constraints that hold in the limit of an incompressible material. Group propagation velocities are derived from the phase velocities to describe the propagation of wave packets away from the excitation region following localized excitation. The theoretical model is compared to the results of finite element (FE) simulations performed using a nearly incompressible material model with the five elasticity constants chosen to preserve the essential properties of the material in the incompressible limit. Propagation velocities calculated from the FE displacement data show complex structure that agrees quantitatively with the theoretical model and demonstrates the possibility of measuring all three elasticity constants needed to characterize an incompressible, TI material.     

A new cell-based FE model for the mechanics of embryonic epithelia

In order to overcome a significant stiffening artefact associated with current finite element (FE) models for the mechanics of embryonic epithelia, two new FE formulations were developed. Cell–cell interfacial tensions γ are represented by constant-force rod elements as in previous models. However, the viscosity of the cytoplasm with its embedded organelles and filament networks is modeled using viscous triangular elements, it is modeled using either radial and circumferential dashpots or an orthogonal dashpot system rather than the viscous triangular elements typical of previous two-dimensional FE models. The models are tested against tissue (epithelium) stretching because it gives rise to significant changes in cell shape and against cell sorting because it involves high rates of cell rearrangement. The orthogonal dashpot system is found to capture cell size and shape effects well, give the model cells characteristics that are consistent with those of real cells, provide high computational efficiency and hold promise for future three-dimensional analyses.     

Soft-cuticle biomechanics: A constitutive model of anisotropy for caterpillar integument

The mechanical properties of soft tissues are important for the control of motion in many invertebrates. Pressurized cylindrical animals such as worms have circumferential reinforcement of the body wall; however, no experimental characterization of comparable anisotropy has been reported for climbing larvae such as caterpillars. Using uniaxial, real-time fluorescence extensometry on millimeter scale cuticle specimens we have quantified differences in the mechanical properties of cuticle to circumferentially and longitudinally applied forces. Based on these results and the composite matrix-fiber structure of cuticle, a pseudo-elastic transversely isotropic constitutive material model was constructed with circumferential reinforcement realized as a Horgan-Saccomandi strain energy function. This model was then used numerically to describe the anisotropic material properties of Manduca cuticle. The constitutive material model will be used in a detailed finite-element analysis to improve our understanding of the mechanics of caterpillar crawling.     

A comparison of cortical elastic properties in the craniofacial skeletons of three primate species and its relevance to the study of human evolution

When a force is applied to an object, the resulting pattern of strain is a function of both the object's geometry and its elastic properties. Thus, knowledge of elastic properties in craniofacial cortical bone is indispensable for exploring the biomechanics and adaptation of primate skulls. However, elastic properties, such as density and stiffness, cannot be measured in all species, particularly extinct species known only from fossils. In order for advanced engineering techniques such as finite element analysis (FEA) to be applied to questions of primate and hominid craniofacial functional morphology, it is important to understand interspecific patterns of variation in elastic properties. We hypothesized that closely related species would have similar patterns of bone elastic properties, and that similarities with extant species should allow reasonable predictions of elastic properties in the skeletons of extinct primate species. In this study, we tested this hypothesis by measuring elastic properties in five areas of the external cortex of the baboon craniofacial skeleton using an ultrasonic technique, and by comparing the results to existing data from macaque and human crania. Results showed that cortical density, thickness, elastic and shear moduli, and anisotropy varied among areas in the baboon cranium. Similar variation had previously been found in rhesus and human crania, suggesting area-specific elastic patterns in the skulls of each species. Comparison among species showed differences, suggesting species-specific patterns. These patterns were more similar between macaques and baboons for density, maximum elastic and shear stiffness, and anisotropy than between either of these and humans. This finding demonstrates that patterns of cortical elastic properties are generally similar in closely related primate species with similar craniofacial morphology. Thus, reasonable estimates of cortical bone elastic properties should be possible for extinct species through the study of phylogenetically related and functionally similar modern forms. For example, reasonable elastic property estimates of cortical bone from fossil hominid skulls should be possible once adequate information about such properties in extant great apes is added to our current data from humans, macaques, and baboons. Such data should eventually allow FEA of craniofacial function in fossil hominids.     

A study of the viscoelastic effect in a bone remodeling model

The present paper addresses the following question can a simple regulatory bone remodeling model predict effects of viscosity on the trabecular morphology? For that, we propose an extension of a previous bone remodeling model by taking into account the viscosity properties of the tissue. Zener’s law is used to describe the mechanical behavior of the bone and a specific law of the apparent bone density rate is proposed. Based on stability analysis, numerical simulations are then performed to investigate the viscosity role on simulations of the bone remodeling process. We show that the viscous contribution affects the evolution of the apparent bone density, by slowing down the adaptation process, which seems to be confirmed by simulations with real data obtained from rat tibia.     

A computationally efficient strategy to estimate muscle forces in a finite element musculoskeletal model of the lower limb

Concurrent multiscale simulation strategies are required in computational biomechanics to study the interdependence between body scales. However, detailed finite element models rarely include muscle recruitment due to the computational burden of both the finite element method and the optimization strategies widely used to estimate muscle forces. The aim of this study was twofold: first, to develop a computationally efficient muscle force prediction strategy based on proportional-integral-derivative (PID) controllers to track gait and chair rise experimental joint motion with a finite element musculoskeletal model of the lower limb, including a deformable knee representation with 12 degrees of freedom; and, second, to demonstrate that the inclusion of joint-level deformability affects muscle force estimation by using two different knee models and comparing muscle forces between the two solutions. The PID control strategy tracked experimental hip, knee, and ankle flexion/extension with root mean square errors below 1°, and estimated muscle, contact and ligament forces in good agreement with previous results and electromyography signals. Differences up to 11% and 20% in the vasti and biceps femoris forces, respectively, were observed between the two knee models, which might be attributed to a combination of differing joint contact geometry, ligament behavior, joint kinematics, and muscle moment arms. The tracking strategy developed in this study addressed the inevitable tradeoff between computational cost and model detail in musculoskeletal simulations and can be used with finite element musculoskeletal models to efficiently estimate the interdependence between muscle forces and tissue deformation.     

New method of three-dimensional analysis of bipedal locomotion for the study of displacements of the body and body-parts centers of mass in man and non-human primates: Evolutionary framework

The current biomechanical interpretation of the chimpanzee's bipedal walking argues that larger lateral and vertical displacements of the body center of mass occur in the chimpanzee's “side-to-side” gait than in the human striding gait. The evolutionary hypothesis underlying this study is the following: during the evolution of human bipedalism one of the necessary changes could have been the progressive reduction of these displacements of the body center of mass. In order to quantitatively test this hypothesis, it is necessary to obtain simultaneously the trajectories of the centers of mass of the whole body and of the different body parts. To solve this problem, a new method of three-dimensional analysis of walking, associated with a volumetric modelling of the body, has been developed based on finite-element modelling. An orthogonal synchrophotographic device yielding four synchronous pictures of the walking subject allows a qualitative analysis of the photographic sequences together with the results of their quantitative analysis. This method was applied to an adult man, a 3-year-old girl and a 9-year-old male chimpanzee. Our results suggest that the trajectory of the body center of mass of the human is distinguished from that of the chimpanzee not by a lower movement amplitude but by the synchronization of the transverse and vertical displacements into two periodic curves in phase with one another. The non-human primate uses its repertoire of arboreal movements in its bipedal terrestrial gait, provisionally referred to as a “rope-walker” gait. We show that the interpretation of a “side-to-side” gait is not applicable to the chimpanzee. We argue that similarly this interpretation and the initial hypothesis presuppose a basic symmetric structure of the gait, in relation to the sagittal plane of progression, similar to the human one. This lateral symmetry of the right and left displacements of the center of gravity, in phase with the right and left single supports of walking, is probably a very derived feature of the human gait. We suggest that low lateral and vertical displacements of the body center of mass are not indicative of a progressive bipedal gait and we discuss the new evolutionary implications of our results. © 1993 Wiley-Liss, Inc.     

A model for the neuronal substrate of dead reckoning and memory in arthropods: a comparative computational and behavioral study

Returning to the point of departure after exploring the environment is a key capability for most animals. In the absence of landmarks, this task will be solved by integrating direction and distance traveled over time. This is referred to as path integration or dead reckoning. An important question is how the nervous systems of navigating animals such as the 1 mm³ brain of ants can integrate local information in order to make global decision. In this article we propose a neurobiologically plausible system of storing and retrieving direction and distance information. The path memory of our model builds on the well established concept of population codes, moreover our system does not rely on trigonometric functions or other complex non-linear operations such as multiplication, but only uses biologically plausible operations such as integration and thresholding. We test our model in two paradigms; in the first paradigm the system receives input from a simulated compass, in the second paradigm, the model is tested against behavioral data recorded from 17 ants. We were able to show that our path memory system was able to reliably encode and compute the angle of the vector pointing to the start location, and that the system stores the total length of the trajectory in a dependable way. From the structure and behavior of our model, we derive testable predictions both at the level of observable behavior as well as on the anatomy and physiology of its underlying neuronal substrate.