What makes an accurate and reliable subject-specific finite element model? A case study of an elephant femur

The Tarsometatarsus of the Ostrich Struthio camelus: Anatomy, Bone Densities, and Structural Mechanics

Background

The ostrich Struthio camelus reaches the highest speeds of any extant biped, and has been an extraordinary subject for studies of soft-tissue anatomy and dynamics of locomotion. An elongate tarsometatarsus in adult ostriches contributes to their speed. The internal osteology of the tarsometatarsus, and its mechanical response to forces of running, are potentially revealing about ostrich foot function.

Methods/Principal Findings

Computed tomography (CT) reveals anatomy and bone densities in tarsometatarsi of an adult and a young juvenile ostrich. A finite element (FE) model for the adult was constructed with properties of compact and cancellous bone where these respective tissues predominate in the original specimen. The model was subjected to a quasi-static analysis under the midstance ground reaction and muscular forces of a fast run. Anatomy–Metatarsals are divided proximally and distally and unify around a single internal cavity in most adult tarsometatarsus shafts, but the juvenile retains an internal three-part division of metatarsals throughout the element. The juvenile has a sparsely ossified hypotarsus for insertion of the m. fibularis longus, as part of a proximally separate third metatarsal. Bone is denser in all regions of the adult tarsometatarsus, with cancellous bone concentrated at proximal and distal articulations, and highly dense compact bone throughout the shaft. Biomechanics–FE simulations show stress and strain are much greater at midshaft than at force applications, suggesting that shaft bending is the most important stressor of the tarsometatarsus. Contraction of digital flexors, inducing a posterior force at the TMT distal condyles, likely reduces buildup of tensile stresses in the bone by inducing compression at these locations, and counteracts bending loads. Safety factors are high for von Mises stress, consistent with faster running speeds known for ostriches.

Conclusions/Significance

High safety factors suggest that bone densities and anatomy of the ostrich tarsometatarsus confer strength for selectively critical activities, such as fleeing and kicking predators. Anatomical results and FE modeling of the ostrich tarsometatarsus are a useful baseline for testing the structure’s capabilities and constraints for locomotion, through ontogeny and the full step cycle. With this foundation, future analyses can incorporate behaviorally realistic strain rates and distal joint forces, experimental validation, and proximal elements of the ostrich hind limb.

Biomechanical implications of intraspecific shape variation in chimpanzee crania: moving toward an integration of geometric morphometrics and finite element analysis

In a broad range of evolutionary studies, an understanding of intraspecific variation is needed in order to contextualize and interpret the meaning of variation between species. However, mechanical analyses of primate crania using experimental or modeling methods typically encounter logistical constraints that force them to rely on data gathered from only one or a few individuals. This results in a lack of knowledge concerning the mechanical significance of intraspecific shape variation that limits our ability to infer the significance of interspecific differences. This study uses geometric morphometric methods (GM) and finite element analysis (FEA) to examine the biomechanical implications of shape variation in chimpanzee crania, thereby providing a comparative context in which to interpret shape-related mechanical variation between hominin species. Six finite element models (FEMs) of chimpanzee crania were constructed from CT scans following shape-space Principal Component Analysis (PCA) of a matrix of 709 Procrustes coordinates (digitized onto 21 specimens) to identify the individuals at the extremes of the first three principal components. The FEMs were assigned the material properties of bone and were loaded and constrained to simulate maximal bites on the P(3) and M(2) . Resulting strains indicate that intraspecific cranial variation in morphology is associated with quantitatively high levels of variation in strain magnitudes, but qualitatively little variation in the distribution of strain concentrations. Thus, interspecific comparisons should include considerations of the spatial patterning of strains rather than focus only on their magnitudes.